Small-aperture Seismic Array Research Articles

Monitoring induced microseismicity in an urban context using very small seismic arrays: The case study of the Vendenheim EGS project

Monitoring the seismicity induced by fluid injections in deep geothermal reservoirs is often limited by the significant anthropogenic ambient noise level inherent to most enhanced geothermal system (EGS) projects in urban contexts. We report on the performance of a monitoring network made of three small aperture (72 m) seismic arrays composed of 21 3D nodes. We test the setup for four months (from 12 December 2020 to 8 April 2021) of the Strasbourg sequence of induced earthquakes (2019–2022) related to the EGS Georhin project (Vendenheim, France). The deployment starts a few days after the MLv = 3.6 induced earthquake of 4 December 2020 and covers the early shut-in period of the wells. We use a beamforming technique to characterize the main noise sources, which consist of slow apparent velocities of surface waves emitted from mobile anthropogenic sources (motorway and railway traffic). We detect events with a phase-weighted stacking method, which is efficient when wavefronts illuminate the arrays with a high apparent velocity. Earthquakes associated with these detections are located using a matched field processing technique. The obtained catalog includes 216 seismic events, which represent four times more events than the reference catalog from Le Bureau Central Sismologique Français–Renass (the national academic agency in charge of seismicity monitoring in France), and a reduction of the completeness magnitude from 0.3 to −0.5. The clustering of the seismicity is analyzed using waveform correlation. The enriched catalog reveals intermittent seismic activity during the slow and continuous decrease in fluid pressure after shut-in.

Acoustic full waveform inversion for 2-D ambient noise source imaging

SUMMARYWe present a method for estimating seismic ambient noise sources by acoustic full waveform inversion (FWI) of interstation cross-correlations. The method is valid at local scales for laterally heterogeneous media, and ambient noise sources confined to the Earth’s surface. Synthetic tests performed using an actual field array geometry, are used to illustrate three unique aspects of our work. First: the method is able to recover noise sources of arbitrary spatial distribution, both within and outside the receiver array, with high fidelity. This holds true for complex velocity models and does not require a good initial guess for inversion, thereby addressing an outstanding issue in the existing research literature. Second: we analyse the extent of biases in source inversion that arise due to inaccurate velocity models. Our findings indicate that source inversion using simplified (e.g. homogeneous) velocity models may work reliably when lateral variations in velocity structure are limited to 5 or 10 per cent in magnitude, but is vitiated by strong variations of 20 per cent or higher, wherein the effect of strong scattering and/or phase distortions become significant. Finally, our technique is implemented without the adjoint method, which is usually inextricably linked to FWI. Inversions are performed using source kernels computed for each receiver pair, and this approach is computationally tractable for real-world problems with small aperture seismic arrays.

Seismic Monitoring of Novaya Zemlya: Progress, Challenges, and Prospects

AbstractThe permanent seismic stations of the European Arctic maintain a detection threshold of around magnitude 2 for events on and around Novaya Zemlya. Events above magnitude 3 are clearly observed by multiple stations at regional and far-regional distances and, with improved travel-time models, can be located with high accuracy. The monitoring capability for smaller magnitude events is dominated by the small aperture seismic arrays ARCES and SPITS. We review the properties of Novaya Zemlya seismic signals on key stations and discuss how empirical signal processing may enhance detection and interpretation of future events in the region. We present a joint probabilistic location for 21 low-magnitude events between 1986 and 2020 in which the joint probability distribution for all events simultaneously exploits both constraints on earlier events from stations no longer in operation and constraints on newer events from more recently deployed stations. Advances in signal processing, enhanced exploitation of archive data, new permanent stations, and comparative multiple event analysis will all contribute both to a more robust and sensitive detection capability and higher confidence in signal interpretation.

One-Year Seismic Survey of the Tectonic CO2-Rich Site of Mefite d'Ansanto (Southern Italy): Preliminary Insights in the Seismic Noise Wavefield.

A passive seismic experiment is carried out at the non-volcanic highly degassing site of Mefite d'Ansanto located at the northern tip of the Irpinia region (southern Italy), where the 1980 MS 6.9 destructive earthquake occurred. Between 2020 and 2021, background seismic noise was recorded by deploying a broadband seismic station and a seismic array composed of seven 1 Hz three-component sensors. Using two different array configurations, we were allowed to explore in detail the 1-20 Hz frequency band of the seismic noise wavefield as well as Rayleigh wave phase velocities in the 400-800 m/s range. Spectral analyses and array techniques were applied to one year of data showing that the frequency content of the signal is very stable in time. High frequency peaks are likely linked to the emission source, whereas at low frequencies seismic noise is clearly correlated to meteorological parameters. The results of this study show that small aperture seismic arrays probe the subsurface of tectonic CO2-rich emission areas and contribute to the understanding of the link between fluid circulation and seismogenesis in seismically active regions.

Long-term analysis of cryoseismic events and associated ground thermal stress in Adventdalen, Svalbard

Abstract. The small-aperture Spitsbergen seismic array (SPITS) has been in continuous operation at Janssonhaugen for decades. The high-Arctic location in the Svalbard archipelago makes SPITS an ideal laboratory for the study of cryoseisms, a nontectonic class of seismic events caused by freeze processes in ice, ice–soil, and ice–rock materials. We extracted a catalog of > 100 000 events from the nearly continuous observation period between 2004 and 2021, characterized by short-duration ground shaking of just a few seconds. This catalog contains two main subclasses where one subclass is related to underground coal-mining activities and the other is inferred to be dominated by frost quakes resulting from thermal-contraction cracking of ice wedges and crack-filling vein ice. This inference is supported by the correspondence between peaks in observed seismicity with peaks in modeled ground thermal stress, based on a Maxwellian thermo-viscoelastic model constrained by borehole observations of ground temperature. The inferred frost quakes appear to be dominated by surface wave energy and SPITS proximal source positions, with three main areas that are associated with dynamic geomorphological features, i.e., erosional scarps and a frozen-debris/solifluction lobe. Seismic stations providing year-round, high-temporal-resolution measurements of ground motion may be highly complementary to satellite remote sensing methods, such as InSAR (interferometric synthetic aperture radar), for studying the dynamics of periglacial environments. The long-term observational record presented in this study, containing tens of thousands of cryoseismic events, in combination with a detailed record of borehole ground temperature observations, provides a unique insight into the spatiotemporal patterns of cryoseisms. The observed patterns may guide the development of models that can be used to understand future changes to cryoseismicity based on projected temperatures.

Monitoring the Microseismicity through a Dense Seismic Array and a Similarity Search Detection Technique: Application to the Seismic Monitoring of Collalto Gas-Storage, North Italy

Seismic monitoring in areas where induced earthquakes could occur is a challenging topic for seismologists due to the generally very low signal to noise ratio. Therefore, the seismological community is devoting several efforts to the development of high-quality networks around the areas where fluid injection and storage and geothermal activities take place, also following the national induced seismicity monitoring guidelines. The use of advanced data mining strategies, such as template matching filters, auto-similarity search, and deep-learning approaches, has recently further fostered such monitoring, enhancing the seismic catalogs and lowering the magnitude of completeness of these areas. In this framework, we carried out an experiment where a small-aperture seismic array was installed within the dense seismic network used for monitoring the gas reservoir of Collalto, in North Italy. The continuous velocimetric data, acquired for 25 days, were analysed through the application of the optimized auto-similarity search technique FAST. The array was conceived as a cost-effective network, aimed at integrating, right above the gas storage site, the permanent high-resolution Collalto Seismic Network. The analysis allowed to detect micro-events down to magnitude Ml = −0.4 within a distance of ~15 km from the array. Our results confirmed that the system based on the array installation and the FAST data analysis might contribute to lowering the magnitude of completeness around the site of about 0.7 units.

The Deployment of the Seismometer to Investigate Ice and Ocean Structure (SIIOS) in Northwest Greenland: An Analog Experiment for Icy Ocean World Seismic Deployments

Abstract In anticipation of future spacecraft missions to icy ocean worlds, the Seismometer to Investigate Ice and Ocean Structure (SIIOS) was funded by National Aeronautics and Space Administration, to prepare for seismologic investigations of these worlds. During the summer of 2018, the SIIOS team deployed a seismic experiment on the Greenland ice sheet situated, approximately, 80 km north of Qaanaaq, Greenland. The seismometers deployed included one Trillium 120 s Posthole (TPH) broadband seismometer, 13 Silicon Audio flight-candidate seismometers, and five Sercel L28 4.5 Hz geophones. Seismometers were buried 1 m deep in the firn in a cross-shaped array centered on a collocated TPH and Silicon Audio instrument. One part of the array consisted of Silicon Audio and Sercel geophones situated 1 m from the center of the array in the ordinal directions. A second set of four Silicon Audio instruments was situated 1 km from the center of the array in the cardinal directions. A mock-lander spacecraft was placed at the array center and instrumented with four Silicon Audio seismometers. We performed an active-source experiment and a passive-listening experiment that lasted for, approximately, 12 days. The active–source experiment consisted of 9–12 sledgehammer strikes to an aluminum plate at 10 separate locations up to 100 m from the array center. The passive experiment recorded the ice-sheet ambient background noise, as well as local and regional events. Both datasets will be used to quantify differences in spacecraft instrumentation deployment strategies, and for evaluating science capabilities for single-station and small-aperture seismic arrays in future geophysical missions. Our initial results indicate that the flight-candidate seismometer performs comparably to the TPH at frequencies above 0.1 Hz and that instruments coupled to the mock-lander perform comparably to ground-based instrumentation in the frequency band of 0.1–10 Hz. For future icy ocean world missions, a deck-coupled seismometer would perform similarly to a ground-based deployment across the most frequency bands.

Rotaphone-CY: The Newest Rotaphone Model Design and Preliminary Results from Performance Tests with Active Seismic Sources.

Rotaphone-CY is a six-component short-period seismograph that is capable of the co-located recording of three translational (ground velocity) components along three orthogonal axes and three rotational (rotation rate) components around the three axes in one device. It is a mechanical sensor system utilizing records from elemental sensors (geophones) arranged in parallel pairs to derive differential motions in the pairs. The pairs are attached to a rigid frame that is anchored to the ground. The model design, the latest one among various Rotaphone designs based on the same principle and presented elsewhere, is briefly introduced. The upgrades of the new model are a 32-bit A/D converter, a more precise placing of the geophones to parallel pairs and a better housing, which protects the instrument from external electromagnetic noise. The instrument is still in a developmental stage. It was tested in a field experiment that took place at the Geophysical Observatory in Fürstenfeldbruck (Germany) in November 2019. Four Rotaphones-CY underwent the huddle-testing phase of the experiment as well as the field-deployment phase, in which the instruments were installed in a small-aperture seismic array of a triangular shape. The preliminary results from this active-source experiment are shown. Rotaphone-CY data are verified, in part, by various approaches: mutual comparison of records from four independent Rotaphone-CY instruments, waveform matching according to rotation-to-translation relations, and comparison to array-derived rotations when applicable. The preliminary results are very promising and they suggest the good functionality of the Rotaphone-CY design. It has been proved that the present Rotaphone-CY model is a reliable instrument for measuring short-period seismic rotations of the amplitudes as small as 10 rad/s.

Identification of the Nature of Seismic Events that Occurred in the East European Platform as Recorded by the Mikhnevo Small-Aperture Seismic Array at Regional Distances

This article summarizes the observations of the Mikhnevo small-aperture seismic array (SSA) for seismic events within the East European Platform based on analysis of a large number of experimental data. A method for recognizing earthquakes and short-lived group blasts, which is based on discriminating the nature of slope for the curve of distribution of the ratio between amplitude spectra of P and Lg waves log(AP/ALg) relative to the frequency of oscillation, is proposed. It is shown that the curve has a certain trend, whose slope is different for explosions and earthquakes. A new criterion, that is, the sign of the coefficient in the equation that connects the values of log(AP/ALg) and the frequency, is proposed. A negative trend with the increase in frequency is typical of short shots, while a positive trend is typical of earthquakes. The method has been tested based on the records of 405 industrial explosions and 24 earthquakes, with epicenters located 350–900 km from the Mikhnevo SSA.

Identification of Induced Seismicity in the Fault Zone of the Korobkovskoye Deposits (Kursk Magnetic Anomaly) Based on Observations of a Small-Aperture Seismic Array

The results of recording and determining the locations of seismic events triggered by industrial explosions conducted by OJSC Kombinat KMAruda at the Korobkovskoe iron ore deposit (Belgorod Region) are presented. Seismic signals were recorded by a small-aperture seismic array located on the surface in the immediate vicinity of the mine workings. The confinement of induced seismicity to the fault zone and the migration of seismic sources along its strike have been established. Such studies on the territory of the Korobkovskoe iron ore deposit have been performed for the first time.

The Multiscale Structure of the Longmen Shan Central Fault Zone from Local and Teleseismic Data Recorded by Short-Period Dense Arrays

ABSTRACTThe Longmen Shan fault zone (FZ), which consists of the back-range, the central, and the front-range faults, acts as the boundary between the Sichuan basin and eastern Tibet. In this study, local and teleseismic waveforms recorded by a 2D small aperture seismic array (176 temporary short-period seismometers) deployed by China University of Geosciences (Beijing) from 22 October to 20 November 2017 and a dense linear seismic array of 16 stations deployed by Geophysical Exploration Center, China Earthquake Administration during July 2008 are used to study the FZ structure by analyzing FZ-trapped waves (FZTWs), the radial-to-vertical amplitude ratio, and travel-time delays. Based on power density spectra analysis, FZTWs from local events with larger amplitudes and longer wavetrains are clearly observed at stations 6002–6003, 6013–6025, and W025–W032. The dispersion measured from trapped waves is quite weak. The near-surface shear velocity structure estimated from the radial-to-vertical amplitude ratios of local initial P waves shows a low-velocity zone around the surface rupture trace. The slight time delay of direct P waves examined from local and teleseismic events indicates a relatively shallow slow structure beneath the arrays. Through the comprehensive analysis of the central FZ, our results suggest a shallow low-velocity zone with a width of ∼150–160 m along the surface rupture trace. Moreover, our P-wave receiver functions reveal that the Moho depth beneath the Longmen Shan FZ is approximately 45 km, and receiver functions at stations located within the surface rupture zone show more complicated waveforms than those off the surface rupture.

Installation and Performance of the Albuquerque Seismological Laboratory Small-Aperture Posthole Array

AbstractThe Global Seismographic Network (GSN) has been used extensively by seismologists to characterize large earthquakes and image deep earth structure. Although the network’s original design goals have been met, the seismological community has suggested that the incorporation of small-aperture seismic arrays at select sites may improve performance of the network and enable new observations. As a pilot study for this concept, we have created a 500 m aperture, nine-element broadband seismic array around the GSN station ANMO (Albuquerque, New Mexico) at the U.S. Geological Survey Albuquerque Seismological Laboratory (ASL). The array was formed by supplementing the secondary borehole seismometer (90 m depth) at ANMO with eight additional 2.6 m posthole sites. Each station’s seismometer was oriented using a fiber optic gyroscope to within 2.0° of north. Data quality, particularly on the vertical components, is excellent with median power levels closely tracking the secondary sensor at ANMO at frequencies lower than 1 Hz. Horizontal component data are more variable at low frequencies (<0.02 Hz), with the type of installation and local geography appearing to strongly influence the amount of tilt-induced noise. Throughout the article, we pose several fundamental questions related to the variability and precision of seismic wavefield measurements that we seek to address with data from this array. In addition, we calculate the array response and show a few examples of using the array to obtain back azimuths of a local event and a continuous narrowband noise source. The apparent velocity of the event across the array is then used to infer the local P-wave velocity at the ASL. Near-real-time data collected from the array along with collocated meteorological, magnetic, and infrasound data are freely available in near-real time from the Incorporated Research Institutions for Seismology Data Management Center.

Recognizing of Microseismic Source Types Based on Small-Aperture Seismic Array Data

Abstract—In this paper, we consider the problem of recording and processing the microseismic data with their subsequent interpretation and making the decisions concerning the probable threats associated with the violation of integrity of terrestrial environment due to mining activity. These decisions are made based on the detection of microseismic sources and recognition of their types by processing the multichannel seismograms from small-aperture seismic arrays. The procedures of seismogram processing include the algorithms for detecting the signals of microseismic sources in the noisy seismograms from the array sensors and the procedures for estimating various parameters of these sources. In this paper, we propose a new detection algorithm which is based on evaluating multiple coherence of the array seismograms and is capable of identifying the waves radiated from microseismic sources with complex focal mechanisms against the diffuse seismic noise. In contrast to the widely used STA/LTA array detector which is aimed at separating plane seismic waves, the proposed detector is intended for the recognition of waves from local seismic events in the complexly structured subsurface environment. Based on the numerical experiments with model and real seismic array data containing signals from the sources with complex focal mechanisms, it is shown that the types of microseismic sources can be recognized with the use of various algorithms of seismic array data analysis: with new algorithms for estimating the coordinates of microseismic sources and with traditional algorithms of space-time spectral analysis.

Seismogenic ancient structures of the central and northern part of the East European platform

The analysis of the location of the epicenters of earthquakes that occurred in the central and northern part of the East European platform in 2009-2016, recorded by the seismic stations of the GS RAS and the small aperture seismic array of IGD RAS Mikhnevo was performed. The results obtained indirectly indicate the seismic activity of the Riphean structures of the region, disturbing the surface of the basement, and their possible activation at the present time. Available data on historical earthquakes also confirm their relevance to paleorifts. It seems important to take into account the position of the ancient aulacogens in assessing the seismic hazard of the East European platform.

Assessment of the Seismic Impact of Industrial Explosions in the Central Part of the East European Platform

The attenuation coefficients of longitudinal and shear waves propagating in the Earth’s upper crust below the central region of the East European Platform have been determined using seismic observations along a profile from the Novogurovskii limestone quarry to the Mikhnevo small-aperture seismic array. Assuming that dissipative properties of the Earth’s crust are invariable, we drew a contour map characterizing the propagation of seismic energy generated by open pit mine blasting in the central part of the East European Platform; the map reflects both historical and recent earthquakes as well. The distribution of seismic impacts of industrial explosions is compared with the distribution of zones of tectonic dislocations revealed by lineament analysis. The map also contains estimates of the intensity of open pit blasting in 2013–2017.

Small-aperture array as a tool to monitor fluid injection- and extraction-induced microseismicity: applications and recommendations

The monitoring of microseismicity during temporary human activities such as fluid injections for hydrofracturing, hydrothermal stimulations or wastewater disposal is a difficult task. The seismic stations often cannot be installed on hard rock, and at quiet places, noise is strongly increased during the operation itself and the installation of sensors in deep wells is costly and often not feasible. The combination of small-aperture seismic arrays with shallow borehole sensors offers a solution. We tested this monitoring approach at two different sites: (1) accompanying a fracking experiment in sedimentary shale at 4 km depth and (2) above a gas field under depletion. The small-aperture arrays were planned according to theoretical wavenumber studies combined with simulations considering the local noise conditions. We compared array recordings with recordings available from shallow borehole sensors and give examples of detection and location performance. Although the high-frequency noise on the 50-m-deep borehole sensors was smaller compared to the surface noise before the injection experiment, the signals were highly contaminated during injection by the pumping activities. Therefore, a set of three small-aperture arrays at different azimuths was more suited to detect small events, since noise recorded on these arrays is uncorrelated with each other. Further, we developed recommendations for the adaptation of the monitoring concept to other sites experiencing induced seismicity.

Velocity Model of Western Volgo-Uralia from Receiver Functions

The results of studying the deep structure of the Earth’s crust and upper mantle in the central part of the Russian platform from receiver functions are presented. The records of teleseismic waves by the Monakovo small-aperture seismic array in the region of the northwestern slope of the Tokmovskii Arch of the Volga–Kama anteclise are used. The modification of the P-receiver function method (Vinnik, 1977) suggested in (Sanina et al., 2014) for analyzing the receiver functions in the regions with a complexly structured upper part of the section and the presence of a thick sedimentary cover is applied. The method is based on separating the high- and low-frequency components of the seismic record and successive reconstruction of the V-s velocity section in the upper part of the crust, which is performed first and, next, the entire deep section of the crust and the mantle down to a depth of ~300 km. The positions of the seismic conversion boundaries in the crust and upper mantle beneath the Monakovo array are determined. The upper mantle velocity section constructed based on the observations at the Mikhnevo array (Sanina et al., 2014) is compared with the world data on the ancient Precambrian platform.

Identification of Local Seismic Events Using the Mikhnevo Seismic Array

A large number of events with sources in the immediate vicinity of an array are usually detected during seismological observations with seismic arrays. These events should be detected and correctly interpreted during processing of seismic array records in order to avoid clogging up the event catalog. This problem can be solved by classifying records of local events by genetic features and creating a databank with the most representative samples. The present paper considers local events recorded using a unique scientific setup, the Mikhnevo small aperture seismic array. Epicenters of local seismic events are located less than 5 km from the center of the array. Seismic responses of acoustic shock waves are also examined. Seismic events caused by anthropogenic sources are identified and classified using cluster, cross-correlation, and wavelet analysis. Events accompanied only by the arrival of surface waves, as well as events represented by body, surface, and acoustic waves, are identified. Shock wave events are classified as a separate category. A small group of supposedly natural weak events is also found. As a result, a databank of waveforms of local seismic events for the Mikhnevo seismic array is established. In the future, this will make it possible to automate their identification when investigating the seismicity of the East European Platform.

Dynamical links between small- and large-scale mantle heterogeneity: Seismological evidence

We identify PKP•PKP scattered waves (also known as P′•P′) from earthquakes recorded at small-aperture seismic arrays at distances less than 65°. P′•P′ energy travels as a PKP wave through the core, up into the mantle, then scatters back down through the core to the receiver as a second PKP. P′•P′ waves are unique in that they allow scattering heterogeneities throughout the mantle to be imaged. We use array-processing methods to amplify low amplitude, coherent scattered energy signals and resolve their incoming direction. We deterministically map scattering heterogeneity locations from the core–mantle boundary to the surface. We use an extensive dataset with sensitivity to a large volume of the mantle and a location method allowing us to resolve and map more heterogeneities than have previously been possible, representing a significant increase in our understanding of small-scale structure within the mantle. Our results demonstrate that the distribution of scattering heterogeneities varies both radially and laterally. Scattering is most abundant in the uppermost and lowermost mantle, and a minimum in the mid-mantle, resembling the radial distribution of tomographically derived whole-mantle velocity heterogeneity. We investigate the spatial correlation of scattering heterogeneities with large-scale tomographic velocities, lateral velocity gradients, the locations of deep-seated hotspots and subducted slabs. In the lowermost 1500 km of the mantle, small-scale heterogeneities correlate with regions of low seismic velocity, high lateral seismic gradient, and proximity to hotspots. In the upper 1000 km of the mantle there is no significant correlation between scattering heterogeneity location and subducted slabs. Between 600 and 900 km depth, scattering heterogeneities are more common in the regions most remote from slabs, and close to hotspots. Scattering heterogeneities show an affinity for regions close to slabs within the upper 200 km of the mantle. The similarity between the distribution of large-scale and small-scale mantle structures suggests a dynamic connection across scales, whereby mantle heterogeneities of all sizes may be directed in similar ways by large-scale convective currents.

Ambient vibration characterization and monitoring of a rock slope close to collapse

We analyse the ambient vibration response of Alpe di Roscioro (AdR), an incipient rock slope failure located above the village Preonzo in southern Switzerland. Following a major failure in May 2012 (volume ∼210 000 m3), the remaining unstable rock mass (∼140 000 m3) remains highly fractured and disrupted, and has been the subject of intensive monitoring. We deployed a small-aperture seismic array at the site shortly after the 2012 failure. The measured seismic response exhibited strong directional amplification (factors up to 35 at 3.5 Hz), higher than previously recorded on rock slopes. The dominant direction of ground motion was found to be parallel to the predominant direction of deformation and perpendicular to open fractures, reflecting subsurface structure of the slope. We then equipped the site with two semi-permanent seismic stations to monitor the seismic response with the goal of identifying changes caused by internal damage that may precede subsequent failure. Although failure has not yet occurred, our data reveal important variations in the seismic response. Amplification factors and resonant frequencies exhibit seasonal trends related (both directly and inversely) to temperature changes and are sensitive to freezing periods (resonant frequencies increase with temperature and during freezing). We attribute these effects to thermal expansion driving microcrack closure, in addition to ice formation, which increase fracture and bulk rock stiffness. We find the site response at AdR is linear over the measured range of weak input motions spanning two orders of magnitude. Our results further develop and refine ambient vibration methods used in rock slope hazard assessment.