Local Earthquakes Research Articles

Integrating earthquake-based passive seismic methods in mineral exploration: Case study from the Gerolekas bauxite mining area, Greece

As the global need for aluminum constantly rises, bauxite is considered to be a critical mineral, and the mining industry is in search of new and effective exploration solutions. In this context, we design and implement a purely earthquake-based passive seismic survey at the Gerolekas bauxite mining site in Greece. It is a very difficult exploration setting, characterized by rough topography, limited accessibility, and a very complex geotectonic regime. We gather a passive seismic data set consisting of four months of continuous recordings (May to August 2018) from 129 stand-alone 3C seismological stations. We then analyze this data set and extract 848 microearthquakes that will serve as sources for the application of local earthquake tomography (LET) and transient-source seismic interferometry (TSI) by autocorrelation. We apply LET to estimate the 3D P- and S-wave velocity models of the subsurface below the study area and TSI by autocorrelation to retrieve the zero-offset virtual reflection responses below each of the recording stations. The velocity models provide a relatively coarse image of a previously completely unexplored part of the mining concession, whereas the higher-resolution virtual reflection imaging illuminates in detail the different interfaces. We also reprocess three lines of legacy active seismic data that were shot in 2003, using the LET P-wave velocity model for depth migration, and confirm the improvement of seismic imaging. Finally, we evaluate the obtained results using well data and jointly interpret them, extracting useful information on the expected target depths and indicating that earthquake-based passive seismic techniques can be an innovative and environmentally friendly option for mineral exploration.

Fluid-assisted intra-plate seismicity at the edge of the Gawler Craton, South Australia

The Australian continent, being void of active plate boundaries, is often perceived as seismically quiescent. However, earthquakes of moderate magnitude (M6+) occur on the continent around once per decade. Such intra-plate activity can pose a significant risk as these earthquakes can occur along fault lines that are either unknown or considered inactive, are often non-periodic, and poorly understood. Within Australia, the spatial distribution of intra-plate seismicity is non-uniform, instead tending to concentrate along certain weak zones of increased activity. One such region is the eastern margin of the Gawler Craton in South Australia, one of the oldest building blocks of the continent. Recently, several new temporary seismic arrays have been deployed in the region, transforming data coverage across southern Australia. In total, 139 new local earthquakes have been recorded, most of which went undetected by the national seismic network. Following relocation, the pattern of earthquakes becomes more localised and appears to coincide with the edge of the Gawler Craton. Further, a spatial association was found between earthquakes and mound springs, which act as the discharge point for groundwater migrating from the Great Artesian Basin. Enhanced fluid pressures (mantle-degassing) within permeable crustal scale fault systems, responding to a regional contractional strain field, appear to be a key driver of seismicity in the area.

Village Design Guidelines as a Tool for Defining Post-Disaster Earthquake Resistant Local Architectural Organization Policies

Due to the destruction of disasters re-establishment of the settlement required to be constructed in a limited time. In the restructuring process of built environment, the primary expectation is to provide risk management against earthquake effects. However, planning decisions developed solely for this purpose may cause the loss of the local characteristics of rural settlements within the borders of the same province vary due to natural and social factors. As a solution, in the scope of the study the earthquake-safe construction model that will develop on reading the potentials of the place is discussed through village design guides. The content approach consists of pre-disaster documenting the local architectural identity data and local earthquake resistant design values determined by inventory analysis. The next step is based on the collection of damage indicators in the post-disaster period to determine the issues requiring improvement. It is envisaged that the parallel development of this content with the legislation will ensure that the guidelines become a design policy in the process of earthquake safe rural construction. Thus, rural planning that is economically, ecologically, and culturally earthquake-resilient can be possible with reference to local architectural identity values.

BPMF: A Backprojection and Matched-Filtering Workflow for Automated Earthquake Detection and Location

Abstract We introduce BPMF (backprojection and matched filtering)—a complete and fully automated workflow designed for earthquake detection and location, and distributed in a Python package. This workflow enables the creation of comprehensive earthquake catalogs with low magnitudes of completeness using no or little prior knowledge of the study region. BPMF uses the seismic wavefield backprojection method to construct an initial earthquake catalog that is then densified with matched filtering. BPMF integrates recent machine learning tools to complement physics-based techniques, and improve the detection and location of earthquakes. In particular, BPMF offers a flexible framework in which machine learning detectors and backprojection can be harmoniously combined, effectively transforming single-station detectors into multistation detectors. The modularity of BPMF grants users the ability to control the contribution of machine learning tools within the workflow. The computation-intensive tasks (backprojection and matched filtering) are executed with C and CUDA-C routines wrapped in Python code. This leveraging of low-level, fast programming languages and graphic processing unit acceleration enables BPMF to efficiently handle large datasets. Here, we first summarize the methodology and describe the application programming interface. We then illustrate BPMF’s capabilities to characterize microseismicity with a 10 yr long application in the Ridgecrest, California area. Finally, we discuss the workflow’s runtime scaling with numerical resources and its versatility across various tectonic environments and different problems.

Preliminary Result of Swarm Activities in Toba Region Using Dense Temporary Network

For the first time, an earthquake swarm occurred in 2021 in Lake Toba (Indonesia), one of the largest calderas in the world. Although the earthquakes were located in the volcanic environment, the swarm activities could be related to the tectonic activities from the Sumatra fault system. The swarm activities occurred at a very shallow depth and were felt with the intensity of II–IV MMI by local people. The objective of this research was to investigate the characteristics of the earthquake swarm in the Toba Caldera from the waveforms recorded by temporary stations in March – June 2021. A total of 26 seismic stations were deployed for two months in May and June 2021. Automatic detection of seismic waves is important because it can be used for earthquake early warning systems and speeds up analysis of origin time, P-, and S arrival times, and earthquake locations. We found characteristics of spectrum waveform from low to high frequency of a non-swarm earthquake, while high frequency seismograms were generated by the swarm earthquake. Furthermore, we successfully captured the hypocentre that was distributed in the SW – NE direction and may figure out the unknown geological system that is responsible for the swarm activities.

New Insights Into the Active Tectonics of the Northern Canadian Cordillera From an Enhanced Earthquake Catalog

AbstractSeismic activity in the Northern Canadian Cordillera is characterized by diffuse earthquakes that extend hundreds of km northwest from the Yakutat collision zone. We use 25 months of broadband seismic data from Mackenzie Mountain Earthscope Project (MMEP), USArray Transportable Array (TA), and permanent Canadian National Seismic Network stations to present a local earthquake catalog with high sensitivity to small regional events. Deep learning techniques are adopted for both seismic phase detection and association. Event relocations are performed to provide well constrained estimates of earthquake depth distributions. Clusters of seismicity spanning the upper crust are located in the central Richardson Mountains, along the Tintina fault, and in the northeast Selwyn Basin. These clusters suggest that the core of the Richardson Anticlinorium is tectonically active and that the Tintina fault is a locus for low levels of active deformation. We interpret seismicity in the northeast Selwyn Basin as primarily occurring in the hanging wall of the Plateau thrust fault and suggest that some combination of localized duplex structures and lithological strength contrasts both within the Selwyn Basin and between abutting Paleozoic shelf sequences may be responsible for seismicity in the Mackenzie Mountain foreland.

Local earthquake tomography of the Cachar fold and thrust belt in lower Assam, Northeast India

Local earthquake tomography of the Cachar fold and thrust belt in lower Assam, Northeast India

Months-long seismicity transients preceding the 2023 MW 7.8 Kahramanmaraş earthquake, Türkiye

Short term prediction of earthquake magnitude, time, and location is currently not possible. In some cases, however, documented observations have been retrospectively considered as precursory. Here we present seismicity transients starting approx. 8 months before the 2023 MW 7.8 Kahramanmaraş earthquake on the East Anatolian Fault Zone. Seismicity is composed of isolated spatio-temporal clusters within 65 km of future epicentre, displaying non-Poissonian inter-event time statistics, magnitude correlations and low Gutenberg-Richter b-values. Local comparable seismic transients have not been observed, at least since 2014. Close to epicentre and during the weeks prior to its rupture, only scarce seismic activity was observed. The trends of seismic preparatory attributes for this earthquake follow those previously documented in both laboratory stick-slip tests and numerical models of heterogeneous earthquake rupture affecting multiple fault segments. More comprehensive earthquake monitoring together with long-term seismic records may facilitate recognizing earthquake preparation processes from other regional deformation transients.

Machine learning versus manual earthquake location workflow: testing LOC-FLOW on an unusually productive microseismic sequence in northeastern Italy

It is an open question whether machine-learning (ML) methods can be trusted in areas where dense and localized seismic networks are in operation, and prompt and accurate detection and location of earthquakes are essential to guide decision-making processes that contribute to seismic-risk-mitigation-strategies, even for very low-magnitude events. To address these concerns, we compare the performance of a widely-used ML phase picker, PhaseNet, integrated with several popular earthquake location methods (included in LOC-FLOW), with the results obtained by the workflow adopted since 2012 by the Collalto Seismic Network, installed to monitor natural and potentially induced microearthquakes nearby an underground gas storage. The tested dataset concerns the most populated microseismic sequence observed so far (374 events, ML⩽2.5, August 2021, Refrontolo, NE-Italy), as its unusual productivity raised some criticalities in the combination of automatic routines, and time-consuming manual revision of phase picks adopted by the standard workflow. LOC-FLOW is able to detect the majority of the events listed in the manually revised catalog, demonstrating its ability to efficiently and accurately build earthquake catalogs from continuous seismic data. We highlight both the advantages and limitations of the ML-picker and recommend the use of template-matching-techniques in the final stage of processing to increase the number of events.

Reference seismic crustal model of the Dinarides

Abstract. Continental collision zones are structurally one of the most heterogeneous areas intermixing various different units within a relatively small space. A good example of this is the Dinarides, a mountain chain situated in the central Mediterranean, where thick carbonates cover older crystalline basement units and remnants of subducted oceanic crust. This is further complicated by the highly variable crustal thickness ranging from 20 to almost 50 km. In terms of spatial extension, this area is relatively small but covers tectonically differentiated domains making, any seismic or geological analysis complex, with significant challenges in areas that lack seismic information on crustal structure. Presently there is no comprehensive 3D crustal model of the Dinarides (and surrounding areas). Using the compilations of previous studies and employing kriging interpolation, we created a vertically and laterally varying crustal model defined on a regular grid for the wider area of the Dinarides, also covering parts of Adriatic Sea and the SW part of the Pannonian Basin. The model is divided by three interfaces, Neogene deposit bottom, carbonate rock complex bottom and Moho discontinuity, with seismic velocities (P and S waves) and density defined at each grid point. To validate the newly derived model, we calculated travel times for an earthquake recorded on several seismic stations in the Dinarides area. The calculated travel times show significant improvement when compared to the simple 1D model used for routine earthquake location in Croatia. The model derived in this work represents the first step towards improving our knowledge of the crustal structure in the complex area of the Dinarides. We hope that the newly assembled model will be useful for all forthcoming studies (e.g., as a starting model for seismic tomography, as a model for earthquake simulations) which require knowledge of the crustal structure.

Weighing Geophysical Data With Trans‐Dimensional Algorithms: An Earthquake Location Case Study

AbstractIn geophysical inverse problems, the distribution of physical properties in an Earth model is inferred from a set of measured data. A necessary step is to select data that are best suited to the problem at hand. This step is performed ahead of solving the inverse problem, generally on the basis of expert knowledge. However, expert‐opinion can introduce bias based on pre‐conceptions. Here we apply a trans‐dimensional algorithm to automatically weigh data on the basis of how consistent they are with the fundamental hypotheses made to solve the inverse problem. We demonstrate this approach by inverting arrival times for the location of a seismic source in an elastic half‐space, assuming a point‐source and uniform weights in concentric shells. The key advantage is that the data do no longer need to be selected by an expert, but they are assigned varying weights during the inversion procedure.

The 21 March 2022 Mw 5.1 Tampen Spur Earthquake, North Sea: Location, Moment Tensor, and Context

ABSTRACT The 21 March 2022 Mw 5.1 Tampen Spur earthquake is the largest event in the North Sea over the last 33 yr. The earthquake was recorded by a 10,708 sensors permanent reservoir monitoring system (PRM) deployed at the Snorre field only a few kilometers away from the epicenter. The event was also recorded on regional networks, as well as at teleseismic distances. Here, we take advantage of this large number of observations to relocate this seismic event, estimate a moment tensor, and refine the depth of the earthquake. The different data sets and location methods yield similar locations, suggesting a well-constrained earthquake location with a depth of 22 ± 5 km. Furthermore, a moment tensor is estimated, which indicates a thrust fault mechanism. Deep crustal earthquakes (>20 km) have been documented within several continental rifts, for example, in Basin and Range in the United States and in northern Germany. Studying deep intraplate earthquakes such as the event observed here can offer insights into the rheology and structures in the deeper portions of the crust. This thrust earthquake implies the reactivation of an older fault in the lower crust consistent with the ridge-push forces from the mid-Atlantic ridge. Such brittle failure at depth indicates a lower crust with little to no fluids present at this depth within the Tampen Spur region.

Structure of shallow magma sources beneath Augustine Volcano (Alaska) inferred from local earthquake tomography

Augustine is an active stratovolcano occupying almost entirely a small island in the Cook Inlet in Alaska. In this study, we use the arrival time data from local seismicity recorded by 15 permanent seismic stations deployed on the flanks of Augustine Volcano by the Alaska Volcano Observatory and invert them by the body wave local earthquake tomography. The resulting model includes the 3D distributions of the P and S wave velocities (Vp and Vs), Vp/Vs ratio and the relocated seismic events. The resolved area extends down to the depth of 2 km below sea level (b.s.l). At shallow depths, within the edifice, we observe a contrasted anomaly with very high Vp, low Vs and very high Vp/Vs, which represents the rigid volcano core composed of igneous rocks (high Vp), but strongly fractured and saturated with fluids and melts (low Vs). At larger depths, below sea level, we retrieve a columnar anomaly of high Vp/Vs ratio representing the magma conduit. On the top of this anomaly, there is a zone with low Vp/Vs coinciding with the strong seismicity cluster, which is associated with the degassing of deep fluids.

Study on VLF Electric Field Anomalies Caused by Seismic Activity on the Western Coast of the Pacific Rim

In order to explore the correlation between earthquakes and ionospheric very low-frequency (VLF) electric field disturbances, this article uses VLF data observed by the China Earthquake Electromagnetic Satellite (CSES) to analyze very low-frequency signals before and after earthquakes from January 2019 to March 2023 in terms of the amplitude and signal-to-noise ratio of electric field power spectrum disturbances. Taking 73 earthquakes with a magnitude of 6.0 or higher occurring in the Circum-Pacific seismic belt as an example, comprehensive research on the VLF electric field disturbance phenomenon caused by strong earthquakes is conducted, considering both the earthquake location and source mechanism. The research results indicate the following: (1) there is a strong correlation between earthquakes with a magnitude of 6.0 or above and abnormal disturbances in the VLF electric field, which often occur within 20 days before the earthquake and within 800 km from the epicenter. (2) From the perspective of earthquake-prone areas, the VLF electric field anomalies observed before earthquakes in the Ryukyu Islands of the Taiwan region exhibit small and concentrated field fluctuations, while the Taiwan Philippines region exhibits larger field fluctuations and more dispersed fluctuations. The discovery of this correlation between seismic ionospheric phenomena and seismic activity provides a new and effective approach to earthquake monitoring, which can be used for earthquake prediction, early warning, and disaster prevention and reduction work.

Co-seismic surface displacement of the June 21, 2022 MW6 Khōst, Afghanistan earthquake from InSAR observations

A robust estimation of the earthquake location, seismic moment, and fault geometry is essential for objective seismic hazard assessment. Seismic events in a remote location, specifically in the absence of seismic and GNSS networks, can be investigated effectively using the InSAR-based technique. This study adopts the Differential Interferometric SAR (DInSAR) technique to quantify the co-seismic surface displacement caused by the June 21, 2022, Khōst MW6 earthquake that occurred along the western plate boundary between the Indian and Eurasian plate. The interferograms show that the maximum surface deformation occurred on the northwest and southwest of the fault line. From coherence, the Line of Sight (LOS) displacement, and the co-seismic surface displacement analysis, it has been observed that surface deformation was most pronounced in the southwest region of the fault line, and the surface has moved to the opposite direction along the fault line, which indicates a sinistral slightly oblique strike-slip movement. This InSAR-based observation appears consistent with the seismic waveforms derived from co-seismic surface displacements. Further, it has been argued that the slip deficit accumulated during the period of the last about 48 years along the frontal region of the northward extension of the Suleiman range and associated fault zone is qualitatively estimated at about 1.5 m, which is consistent with the seismic waveforms derived finite slip model.

The micro‐seismicity of Co. Donegal (Ireland): Defining baseline seismicity in a region of slow lithospheric deformation

AbstractA catalogue of precisely located micro‐seismicity is fundamental for investigating seismicity and rock physical properties in active tectonic and volcanic regions and for the definition of a ‘baseline’ seismicity, required for a safe future exploitation of georesource areas. In this study, we produce the first manually revised catalogue of micro‐seismicity for Co. Donegal region (Ireland), an area of about 50K M2 of on‐going deformation, aimed at localizing natural micro‐seismic events occurred between 2012 and 2015. We develop a stochastic method based on a Markov chain Monte Carlo (McMC) sampling approach to compute earthquake hypocentral location parameters. Our results indicates that micro‐seismicity is present with magnitudes lower than 2 (the highest magnitude is 2.8).The recorded seismicity is almost clustered along previously mapped NE‐SW trending, steeply dipping faults and confined within the upper crust (focal depth less than 10 km). We also recorded anthropogenic seismicity mostly related to quarries' activity in the study area.

Velocity structure (1D) and earthquakes relocation in the flat-slab to normal plate transition zone in the Argentinean transarc

The Argentinean Andean foreland region between latitudes 31°S to 33.5°S presents an intense crustal and upper mantle seismic activity that was recorded by a local seismological network called CHARSME (Chile Argentina SeisMic Experiment), and was used to obtain new velocity determinations. In a first moment, from 4 P-wave velocity (VP) models proposed by other authors, the detail in the layers was increased, resulting in better models used as initial model in the investment. Finally, by carefully selecting a set of earthquakes, a simultaneous velocity inversion, hypocentre determination and station correction were performed using the VELEST code. In order to obtain better velocity models (VP and VS) and VP/VS ratio, different tests were carried out in all processes.The final VP model presented in this paper incorporates considerations made in other seismological studies. This model accurately determined velocities between 10 and 65 km depth. The upper mantle velocity (VP) was defined in 8.05 km/s. From the Wadatti diagram, which represents a line whose slope is related to the wave velocities and the Poisson's coefficient, a VP/VS value of 1.732 was obtained. Alternative VP/VS values were also considered to invert VP, VS and VP/VS simultaneously. The minimum detailed model of VP, VS and VP/VS decreased the RMS by 25% relative to the initial models, implying an improvement in the uncertainty of locations. The new and precise locations of crustal earthquakes allowed defining clusters within Cuyania, a tectonic terrain that collided to Gondwana in the Early Paleozoic. Although it is known that foreland seismicity is mainly concentrated in this terrain, this work shows that there are well-defined zones within Cuyania that are more seismically active. Earthquakes of intermediate depth (upper mantle) are strongly agglomerated in the flat-slab zone of the subducting Nazca Plate.The discontinuities present in our model suggest upper, middle and lower crustal structuring, as well as a marked contrast between the sediments and the underlying basement. It was possible to interpret an upper crust consisting of sediments and basement to a depth of 25 km, a middle crust between 25 and 45 km and a lower crust between 45 and 55 km. We also defined the Mohorovičić (Moho) discontinuity at 55 km depth. The crustal velocities are in agreement with previous seismological studies.

Earthquake Early Warning Starting From 3 s of Records on a Single Station With Machine Learning

AbstractWe introduce the Ensemble Earthquake Early Warning System (E3WS), a set of Machine Learning (ML) algorithms designed to detect, locate, and estimate the magnitude of an earthquake starting from 3 s of P‐waves recorded by a single station. The system is made of six Ensemble ML algorithms trained on attributes computed from ground acceleration time series in the temporal, spectral, and cepstral domains. The training set comprises data sets from Peru, Chile, Japan, and the STEAD global data set. E3WS consists of three sequential stages: detection, P‐phase picking, and source characterization. The latter involves magnitude, epicentral distance, depth, and back azimuth estimation. E3WS achieves an overall success rate in the discrimination between earthquakes and noise of 99.9%, with no false positive (noise mis‐classified as earthquakes) and very few false negatives (earthquakes mis‐classified as noise). All false negatives correspond to M ≤ 4.3 earthquakes, which are unlikely to cause any damage. For P‐phase picking, the Mean Absolute Error is 0.14 s, small enough for earthquake early warning purposes. For source characterization, the E3WS estimates are virtually unbiased, have better accuracy for magnitude estimation than existing single‐station algorithms, and slightly better accuracy for earthquake location. By updating estimates every second, the approach gives time‐dependent magnitude estimates that follow the earthquake source time function. E3WS gives faster estimates than present alert systems relying on multiple stations, providing additional valuable seconds for potential protective actions.

Constructing a 3‐D Radially Anisotropic Crustal Velocity Model for Oklahoma Using Full Waveform Inversion

AbstractOver the past decade, the seismicity rate in the state of Oklahoma has increased significantly, which has been linked to industrial operations, such as saltwater injection and hydraulic fracturing. Taking advantage of induced earthquakes and recently deployed seismometers, we construct a 3‐D radially anisotropic seismic velocity model for the crust of Oklahoma by using full waveform inversion. To mitigate the well‐known cycle‐skipping problem, we use misfit functions based on phase and waveform differences in several frequency bands. Relative velocity perturbations in the inverted model allow us to delineate major geological provinces in Oklahoma, such as the Anadarko Basin and the Cherokee Platform/Shelf. In addition, radial anisotropy in the inverted model reflects deformation within the crust of Oklahoma, which might correlate with sedimentary layering, microcracks/fractures, as well as dominant orientation of anisotropic minerals. The crystalline basement beneath Oklahoma can be inferred from the new velocity model, which enables us to better classify induced seismicity in current earthquake catalogs. Furthermore, synthetic experiments suggest that the new velocity model enables us to better constrain earthquake locations in Oklahoma, especially for determining their depths, which are important for investigating induced seismicity.

PLATE TECTONICS, LARGEST EARTHQUAKES AND RESULTS OF GEODETIC MEASUREMENTS

The article presents the hypothesis of lithospheric plate tectonics and its significance in the location of earthquake and volcanic eruptions on Earth. Information is given on the largest earthquakes in terms of energy and the number of victims of earthquakes on our planet with reference to plate tectonics, including earthquakes in Turkey and Syria on February 6, 2023. An increase in the role of geodesy in the study of geodynamic phenomena in recent decades is shown in connection with the introduction of satellite measuring instruments, which make it possible to determine distances of hundreds and thousands of kilometers with an accuracy of the order of 10-8–10-9 parameters of the Earth, as well as to determine the propagation of deformations from strong earthquakes on the earth's surface. It is concluded that the strongest earthquakes provoke deformations of the earth's crust, covering vast territories of the Earth with an area of hundreds of thousands and millions of square km.