Constraints on fluid flow processes in the Hellenic Accretionary Complex (eastern Mediterranean Sea) from numerical modeling

The dynamics of accretionary convergent margins are severely influenced by intense deformation and fluid expulsion. To quantify the fluid pressure and fluid flow velocities in the Hellenic subduction system, we set up 2-D hydrogeological numerical models following two seismic reflection lines across the Mediterranean Ridge. These profiles bracket the along-strike variation in wedge geometry: moderate compression and a >4 km thick underthrust sequence in the west versus enhanced compression and <1 km of downgoing sediment in the center. Input parameters were obtained from preexisting geophysical data, drill cores, and new geotechnical laboratory experiments. A permeability-porosity relationship was determined by a sensitivity analysis, indicating that porosity and intrinsic permeability are small. This hampers the expulsion of fluids and leads to the build up of fluid overpressure in the deeper portion of the wedge and in the underthrust sediment. The loci of maximum fluid pressure are mainly controlled by the compactional fluid source, which generally decreases toward the backstop. However, pore pressure is still high at the decollement level at distances <100 km from the deformation front, either by the incorporation of low permeability evaporites or additional compaction of the wedge sediments in the two profiles. In the west, however, formation of a wide accretionary complex is facilitated by high pore pressure zones. When compared to other large accretionary complexes such as Nankai or Barbados, our results not only show broad similarities but also that near-lithostatic pore pressures may be easier to maintain in the Hellenic Arc because of accentuated collision, some underthrust evaporates, and a thicker underthrust sequence

The Devastating 2022 M6.2 Afghanistan Earthquake: Challenges, Processes, and Implications

On June 21st, a Mw6.2 earthquake struck the Afghan-Pakistan-border-region, situated within the India-Asia collision. Thousand thirty-nine deaths were reported, making the earthquake the deadliest of 2022. We investigate the event\u27s rupture processes by combining seismological and geodetic observations, aiming to understand what made it that fatal. Our Interferometric Synthetic Aperture Radar-constrained slip-model and regional moment-tensor inversion, confirmed through field observations, reveal a sinistral rupture with maximum slip of 1.8 m at 5 km depth on a N20°E striking, sub-vertical fault. We suggest that not only external factors (event-time, building stock) but fault-specific factors made the event excessively destructive. Surface rupture was favored by the rock foliation, coinciding with the fault strike. The distribution of Peak-Ground-Velocity was governed by the sub-vertical fault. Maximum slip was large compared to other events globally and might have resulted in peak-frequencies coinciding with resonance-frequencies of the local buildings and demonstrates the devastating impact of moderate-size earthquakes

Highly variable friction and slip observed at Antarctic ice stream bed

The slip of glaciers over the underlying bed is the dominant mechanism governing the migration of ice from land into the oceans, with accelerating slip contributing to sea-level rise. Yet glacier slip remains poorly understood, and observational constraints are sparse. Here we use passive seismic observations to measure both frictional shear stress and slip at the bed of the Rutford Ice Stream in Antarctica using 100,000 repetitive stick-slip icequakes. We find that basal shear stresses and slip rates vary from 10$^4$ to 10$^7$ Pa and 0.2 to 1.5 m per day, respectively. Friction and slip vary temporally over the order of hours, and spatially over 10s of metres, due to corresponding variations in effective normal stress and ice–bed interface material. Our findings suggest that the bed is substantially more complex than currently assumed in ice stream models and that basal effective normal stresses may be significantly higher than previously thought. Our observations can provide constraints on the basal boundary conditions for ice-dynamics models. This is critical for constraining the primary contribution of ice mass loss in Antarctica and hence for reducing uncertainty in sea-level rise projections

The crust in the pamir: Insights from receiver functions

The Cenozoic convergence between India and Asia has created Earth's thickest crust in the Pamir‐Tibet Plateau by extreme crustal shortening. Here we study the crustal structure of the Pamir and western Tian Shan, the adjacent margins of the Tajik, Tarim, and Ferghana Basins, and the Hindu Kush, using data collected by temporary seismic experiments. We derive, compare, and combine independent observations from P and S receiver functions. The obtained Moho depth varies from ~40 km below the basins to a double‐normal thickness of 65–75 km underneath the Pamir and Hindu Kush. A Moho doublet—with the deeper interface down to a depth of ~90 km—coincides with the arc of intermediate‐depth seismicity underneath the Pamir, where Asian continental lower crust delaminates and rolls back. The crust beneath most of the Central and South Pamir has a low Vp/Vs ratio (<1.70), suggesting a dominantly felsic composition, probably a result of delamination/foundering of the mafic rocks of the lower crust. Beneath the Cenozoic gneiss domes of the Central and South Pamir, which represent extensional core complexes, the Vp/Vs ratios are moderate to high (~1.75), consistent with the previously observed, midcrustal low‐velocity zones, implying the presence of crustal partial melts. Even higher crustal average Vp/Vs ratios up to 1.90 are found in the sedimentary basins and along the Main Pamir Thrust. The ratios along the latter—the active thrust front of the Pamir—may reflect fluid accumulations within a strongly fractured fault system

Toeplitz operators with radial symbols on weighted holomorphic Orlicz space

Peer reviewe

The devastating 2022 M6.2 Afghanistan earthquake: Challenges, processes, and implications

On June 21st, a Mw6.2 earthquake struck the Afghan-Pakistan-border-region, situated within the India-Asia collision. Thousand thirty-nine deaths were reported, making the earthquake the deadliest of 2022. We investigate the event's rupture processes by combining seismological and geodetic observations, aiming to understand what made it that fatal. Our Interferometric Synthetic Aperture Radar-constrained slip-model and regional moment-tensor inversion, confirmed through field observations, reveal a sinistral rupture with maximum slip of 1.8 m at 5 km depth on a N20°E striking, sub-vertical fault. We suggest that not only external factors (event-time, building stock) but fault-specific factors made the event excessively destructive. Surface rupture was favored by the rock foliation, coinciding with the fault strike. The distribution of Peak-Ground-Velocity was governed by the sub-vertical fault. Maximum slip was large compared to other events globally and might have resulted in peak-frequencies coinciding with resonance-frequencies of the local buildings and demonstrates the devastating impact of moderate-size earthquakes

Lithospheric delamination beneath the southern Puna plateau resolved by local earthquake tomography

We present a local earthquake tomography to illuminate the crustal and uppermost mantle structure beneath the southern Puna plateau and to test the delamination hypothesis. Vp and Vp/Vs ratios were obtained using travel time variations recorded by 75 temporary seismic stations between 2007 and 2009. In the upper crust, prominent low Vp anomalies are found beneath the main volcanic centers, indicating the presence of magma and melt beneath the southern Puna plateau. Beneath the Moho at around 90-km depth, a strong high Vp anomaly is detected just west of the giant backarc Cerro Galan ignimbrite caldera. This high Vp anomaly is only resolved if earthquakes with an azimuthal gap up to 300 degrees are included in the inversion. However, we show through data subset and synthetic tests that the anomaly is robust due to our specific station-event geometry and interpret it as a delaminated block of lower crust and uppermost mantle lithosphere under the southern Puna plateau. The low velocities in the crust are interpreted as a product of the delamination event that triggered the rise of fluids and melts into the crust and induced the high topography in this part of the plateau. The tomography also reveals the existence of low-velocity anomalies that link arc magmatism at the Ojos del Salado volcanic center with slab seismicity clusters at depths of about 100 and 150 km and support fluid transport in the mantle wedge due to dehydration reaction within the subducted slab

Distributed Acoustic Sensing (DAS) for natural microseismicity studies: A case study from Antarctica

Icequakes, microseismic earthquakes at glaciers, offer insights into the dynamics of ice sheets. For the first time in the Antarctic, we explore the use of fiber optic cables as Distributed Acoustic Sensors (DAS) as a new approach for monitoring basal icequakes. We present the use of DAS for studying icequakes as a case study for the application of DAS to microseismic datasets in other geological settings. Fiber was deployed on the ice surface at Rutford Ice Stream in two different configurations. We compare the performance of DAS with a conventional geophone network for: microseismic detection and location; resolving source and noise spectra; source mechanism inversion; and measuring anisotropic shear-wave splitting parameters. Both DAS array geometries detect fewer events than the geophone array. However, DAS is superior to geophones for recording the microseism signal, suggesting the applicability of DAS for ambient noise interferometry. We also present the first full-waveform source mechanism inversions using DAS anywhere, successfully showing the horizontal stick-slip nature of the icequakes. In addition, we develop an approach to use a 2D DAS array geometry as an effective multi-component sensor capable of accurately characterising shear-wave splitting due to the anisotropic ice fabric. Although our observations originate from a glacial environment, the methodology and implications of this work are relevant for employing DAS in other microseismic environments