Refine
Year of publication
Language
- English (77)
Is part of the Bibliography
- yes (77)
Keywords
- Earthquake source observations (11)
- Seismicity and tectonics (6)
- Induced seismicity (5)
- Time-series analysis (5)
- Statistical seismology (4)
- Earthquake dynamics (3)
- Fracture and flow (3)
- fault (3)
- induced seismicity (3)
- Body waves (2)
From June to August 2021, we deployed a dense seismic nodal network across the Hengill geothermal area in southwest Iceland to image and characterize faults and high-temperature zones at high resolution.
The nodal network comprised 498 geophone nodes spread across the northern Nesjavellir and southern Hverahlio geothermal fields and was complemented by an existing permanent and temporary backbone seismic network of a total of 44 short-period and broadband stations.
In addition, we recorded distributed acoustic sensing data along two fiber optic telecommunication cables near the Nesjavellir geothermal power plant with commercial interrogators.
During the time of deployment, a vibroseis survey took place around the Nesjavellir power plant.
Here, we describe the network and the recorded datasets.
Furthermore, we showsome initial results that indicate a high data quality and highlight the potential of the seismic records for various follow up studies, such as high-resolution event location to delineate faults and body- and surface-wave tomographies to image the subsurface velocity structure in great detail.
Natural gas can be temporarily stored in a variety of underground facilities, such as depleted gas and oil fields, natural aquifers and caverns in salt rocks. Being extensively monitored during operations, these systems provide a favourable opportunity to investigate how pressure varies in time and space and possibly induces/triggers earthquakes on nearby faults. Elaborate and detailed numerical modelling techniques are often applied to study gas reservoirs. Here we show the possibilities and discuss the limitations of a flexible and easily formulated tool that can be straightforwardly applied to simulate temporal pore-pressure variations and study the relation with recorded microseismic events. We use the software POEL (POroELastic diffusion and deformation) which computes the poroelastic response to fluid injection/extraction in a horizontally layered poroelastic structure. We further develop its application to address the presence of vertical impermeable faults bounding the reservoir and of multiple injection/extraction sources. Exploiting available information on the reservoir geometry and physical parameters, and records of injection/extraction rates for a gas reservoir in southern Europe, we perform an extensive parametric study considering different model configurations. Comparing modelled spatiotemporal pore-pressure variations with in situ measurements, we show that the inclusion of vertical impermeable faults provides an improvement in reproducing the observations and results in pore-pressure accumulation near the faults and in a variation of the temporal pore-pressure diffusion pattern. To study the relation between gas storage activity and recorded local microseismicity, we applied different seismicity models based on the estimated porepressure distribution. This analysis helps to understand the spatial distribution of seismicity and its temporal modulation. The results show that the observed microseismicity could be partly linked to the storage activity, but the contribution of tectonic background seismicity cannot be excluded.
Reservoir-triggered seismicity has been observed near dams during construction, impoundment, and cyclic filling in many parts of the earth. In Turkey, the number of dams has increased substantially over the last decade, with Ataturk Dam being the largest dam in Turkey with a total water capacity of 48.7 billion m(3). After the construction of the dam, the monitoring network has improved. Considering earthquakes above the long-term completeness magnitude of M-C = 3.5, the local seismicity rate has substantially increased after the filling of the reservoir. Recently, two damaging earthquakes of M-w 5.5 and M-w 5.1 occurred in the town of Samsat near the Ataturk Reservoir in 2017 and 2018, respectively. In this study, we analyze the spatio-temporal evolution of seismicity and its source properties in relation to the temporal water-level variations and the stresses resulting from surface loading and pore-pressure diffusion. We find that water-level and seismicity rate are anti-correlated, which is explained by the stabilization effect of the gravitational induced stress imposed by water loading on the local faults. On the other hand, we find that the overall effective stress in the seismogenic zone increased over decades due to pore-pressure diffusion, explaining the enhanced background seismicity during recent years. Additionally, we observe a progressive decrease of the Gutenberg-Richter b-value. Our results indicate that the stressing rate finally focused on the region where the two damaging earthquakes occurred in 2017 and 2018.
The Alpine mountains in central Europe are characterized by a heterogeneous crust accumulating different tectonic units and blocks in close proximity to sedimentary foreland basins. Centroid moment tensor inversion provides insight into the faulting mechanisms of earthquakes and related tectonic processes but is significantly aggravated in such an environment. Thanks to the dense AlpArray seismic network and our flexible bootstrap-based inversion tool Grond, we are able to test different setups with respect to the uncertainties of the obtained moment tensors and centroid locations. We evaluate the influence of frequency bands, azimuthal gaps, input data types, and distance ranges and study the occurrence and reliability of non-double-couple (DC) components. We infer that for most earthquakes (M-w >= 3.3) a combination of time domain full waveforms and frequency domain amplitude spectra in a frequency band of 0.02-0.07 Hz is suitable. Relying on the results of our methodological tests, we perform deviatoric moment tensor (MT) inversions for events with M-w > 3.0. Here, we present 75 solutions for earthquakes between January 2016 and December 2019 and analyze our results in the seismotectonic context of historical earthquakes, seismic activity of the last 3 decades, and GNSS deformation data. We study regions of comparably high seismic activity during the last decades, namely the Western Alps, the region around Lake Garda, and the eastern Southern Alps, as well as clusters further from the study region, i.e., in the northern Dinarides and the Apennines. Seismicity is particularly low in the Eastern Alps and in parts of the Central Alps. We apply a clustering algorithm to focal mechanisms, considering additional mechanisms from existing catalogs. Related to the N-S compressional regime, E-W-to-ENE-WSW-striking thrust faulting is mainly observed in the Friuli area in the eastern Southern Alps. Strike-slip faulting with a similarly oriented pressure axis is observed along the northern margin of the Central Alps and in the northern Dinarides. NW-SE-striking normal faulting is observed in the NW Alps, showing a similar strike direction to normal faulting earthquakes in the Apennines. Both our centroid depths and hypocentral depths in existing catalogs indicate that Alpine seismicity is predominantly very shallow; about 80% of the studied events have depths shallower than 10 km.
On 12 August 2021, a > 220 s lasting complex earthquake with Mw > 8.2 hit the South Sandwich Trench. Due to its remote location and short interevent times, reported earthquake parameters varied significantly between different international agencies. We studied the complex rupture by combining different seismic source characterization techniques sensitive to different frequency ranges based on teleseismic broadband recordings from 0.001 to 2 Hz, including point and finite fault inversions and the back-projection of high-frequency signals.
We also determined moment tensor solutions for 88 aftershocks. The rupture initiated simultaneously with a rupture equivalent to a M-w 7.6 thrust earthquake in the deep part of the seismogenic zone in the central subduction interface and a shallow megathrust rupture, which propagated unilaterally to the south with a very slow rupture velocity of 1.2 km/s and varying strike following the curvature of the trench.
The slow rupture covered nearly two-thirds of the entire subduction zone length, and with M-w 8.2 released the bulk of the total moment of the whole earthquake.
Tsunami modeling indicates the inferred shallow rupture can explain the tsunami records. The southern segment of the shallow rupture overlaps with another activation of the deeper part of the megathrust equivalent to M-w 7.6. The aftershock distribution confirms the extent and curvature of the rupture. Some mechanisms are consistent with the mainshocks, but many indicate also activation of secondary faults. Rupture velocities and radiated frequencies varied strongly between different stages of the rupture, which might explain the variability of published source parameters.
Plain Language Summary
The earthquake of 12 August 2021 along the deep-sea trench of the South Sandwich Islands in the South Atlantic reached a magnitude of 8.2 and triggered a tsunami. The automatic earthquake parameter determination of different agencies showed very different results shortly after the earthquake and partially underestimated the tsunami potential of the earthquake.
A possible reason was the complex rupture process and that the tsunami was generated by a long and shallow slow slip rupture sandwiched between more conventional fast slip subevents at its northern and southern ends. In addition, the fault surface, which extended over 450 km, was highly curved striking 150 degrees-220 degrees.
We investigated the different components of the seismic wavefields in different frequency ranges and with different methods.
The analysis shows how even complex earthquakes can be deciphered by combining analyzing methods. The comparison with aftershocks and the triggered tsunami waves confirms our model that explains the South Sandwich rupture by four subevents in the plate boundary along the curved deep-sea trench. Here, the depth, rupture velocities, and slip on each segment of the rupture vary considerably.
The method can also be applied to other megathrust earthquakes and help to further improve tsunami warnings in the future.
With the present study, we introduce a fast and robust method to calculate the source displacement spectra of small earthquakes on a local to regional scale. The work is based on the publicly available Qopen method of full envelope inversion, which is further tuned for the given purpose. Important source parameters-seismic moment, moment magnitude, corner frequency, and high-frequency fall off-are determined from the source spectra by fitting a simple earthquake source model. The method is demonstrated by means of a data set comprising the 2018 West Bohemia earthquake swarm. We report moment magnitudes, corner frequencies, and centroid moment tensors inverted from short-period body waves with the Grond package for all earthquakes with a local magnitude larger than 1.8. Moment magnitudes calculated by envelope inversion show a very good agreement to moment magnitudes resulting from the probabilisitc moment tensor inversion. Furthermore, source displacement spectra from envelope inversion show a good agreement with spectra obtained by multiple taper analysis of the direct onsets of body waves but are not affected by the large scatter of the second. The seismic moments obtained with the envelope inversion scale with corner frequencies according to M-0 proportional to f(c)(-4.7). Earthquakes of the present data set result in a smaller stress drop for smaller magnitudes. Self-similarity of earthquake rupture is not observed. In addition, we report frequency-dependent site amplification at the used stations.
Ground subsidence caused by natural or anthropogenic processes affects major urban areas worldwide. Sinkhole formation and infrastructure fractures have intensified in the federal capital of Maceio (Alagoas, Brazil) since early 2018, forcing authorities to relocate affected residents and place buildings under demolition. In this study, we present a 16-year history (2004-2020) of surface displacement, which shows precursory deformations in 2004-2005, reaching a maximum cumulative subsidence of approximately 200 cm near the Mundau Lagoon coast in November 2020. By integrating the displacement observations with numerical source modelling, we suggest that extensive subsidence can be primarily associated with the removal of localized, deep-seated material at the location and depth where salt is mined. We discuss the accelerating subsidence rates, influence of severe precipitation events on the aforementioned geological instability, and related hazards. This study suggests that feedback destabilization mechanisms may arise in evaporite systems due to anthropogenic activities, fostering enhanced and complex superficial ground deformation.
Seismicity models are probabilistic forecasts of earthquake rates to support seismic hazard assessment.
Physics-based models allow extrapolating previously unsampled parameter ranges and enable conclusions on underlying tectonic or human-induced processes.
The Coulomb Failure (CF) and the rate-and-state (RS) models are two widely used physics-based seismicity models both assuming pre-existing populations of faults responding to Coulomb stress changes.
The CF model depends on the absolute Coulomb stress and assumes instantaneous triggering if stress exceeds a threshold, while the RS model only depends on stress changes.
Both models can predict background earthquake rates and time-dependent stress effects, but the RS model with its three independent parameters can additionally explain delayed aftershock triggering.
This study introduces a modified CF model where the instantaneous triggering is replaced by a mean time-to-failure depending on the absolute stress value.
For the specific choice of an exponential dependence on stress and a stationary initial seismicity rate, we show that the model leads to identical results as the RS model and reproduces the Omori-Utsu relation for aftershock decays as well stress-shadowing effects.
Thus, both CF and RS models can be seen as special cases of the new model. However, the new stress response model can also account for subcritical initial stress conditions and alternative functions of the mean time-to-failure depending on the problem and fracture mode.
On October 9, 2014, a Mw 7.1-6.7 seismic doublet occurred at the Juan Fernandez microplate, close to the triple junction with Pacific and Nazca plates. The Mw 7.1 earthquake is the largest earthquake ever to have been recorded in the region. Its thrust focal mechanism is also unusual for the region, although the northern part of the microplate is expected to undergo compression. The region is remote and seismological data is limited to a seismic station at similar to 600 km distance on Easter Island and teleseismic observations for the largest events. We use a combination of advanced seismological techniques to overcome the lack of local data and resolve earthquake source parameters for the doublet and its aftershock sequence, being able to reconstruct the chronology of the sequence and the geometry of affected fault segments. Our results depict a complex seismic sequence characterized by the interplay of thrust and strike-slip earthquakes along different structures, including a second, reversed strike slip-thrust seismic doublet in November 2014. Seismicity occurred within the microplate and only in the late part of the sequence migrated northward, towards the microplate boundary. The first largest doublet, whose rupture kinematic is well explained by stress changes imparted by the first subevent on the second one, may have activated unmapped E-W and NE-SW faults or an internal curved pseudofault, attributed to the longterm rotation of the microplate. Few large, thrust earthquakes are observed within the sequence, taking place in the vicinity of mapped compressional ridges. We suggest that compressional stresses in the northern part of the microplate and at its boundary are partially accommodated aseismically. However, the occasional occurrence of large, impulsive thrust earthquakes, with a considerable tsunamigenic potential, poses a relevant hazard for islands in the South Pacific region.