Major Slip Research Articles

Magnitude Clustering During Stick‐Slip Dynamics on Laboratory Faults

AbstractWe present an analysis of magnitude clustering of microfractures inferred from acoustic emissions (AEs) during stick‐slip (SS) dynamics of faulted Westerly granite samples in frictional sliding experiments, with and without fluids, under triaxial loading with constant displacement rate. We investigate magnitude clustering in time across periods during, preceding and after macroscopic slip events on laboratory faults. Our findings reveal that magnitude clustering exists such that subsequent AEs tend to have more similar magnitudes than expected. Yet, this clustering only exists during macroscopic slip events and is strongest during major slip events in fluid‐saturated and dry samples. We demonstrate that robust magnitude clustering arises from variations in frequency‐magnitude distributions of AE events during macroscopic slip events. These temporal variations indicate a prevalence of larger AE events right after (0.3–3 s) the SS onset. Hence, magnitude clustering is a consequence of non‐stationarities.

Intelligent Tire Prototype in Longitudinal Slip Operating Conditions.

With the recent advances in autonomous vehicles, there is an increasing need for sensors that can help monitor tire-road conditions and the forces that are applied to the tire. The footprint area of a tire that makes direct contact with the road surface, known as the contact patch, is a key parameter for determining a vehicle's effectiveness in accelerating, braking, and steering at various velocities. Road unevenness from features such as potholes and cracks results in large fluctuations in the contact patch surface area. Such conditions can eventually require the driver to perform driving maneuvers unorthodox to normal traffic patterns, such as excessive pedal depressions or large steering inputs, which can escalate to hazards such as the loss of control or impact. The integration of sensors into the inner liner of a tire has proven to be a promising method for extracting real-time tire-to-road contact patch interface data. In this research, a tire model is developed using Abaqus/CAE and analyzed using Abaqus/Explicit to study the nonlinear behavior of a rolling tire. Strain variations are investigated at the contact patch in three major longitudinal slip driving scenarios, including acceleration, braking, and free-rolling. Multiple vertical loading conditions on the tire are applied and studied. An intelligent tire prototype called KU-iTire is developed and tested to validate the strain results obtained from the simulations. Similar operating and loading conditions are applied to the physical prototype and the simulation model such that valid comparisons can be made. The experimental investigation focuses on the effectiveness of providing usable and reliable tire-to-road contact patch strain variation data under several longitudinal slip operating conditions. In this research, a correlation between FEA and experimental testing was observed between strain shape for free-rolling, acceleration, and braking conditions. A relationship between peak longitudinal strain and vertical load in free-rolling driving conditions was also observed and a correlation was observed between FEA and physical testing.

Phenotypic analysis of ataxia in spinocerebellar ataxia type 6 mice using DeepLabCut

This study emphasizes the benefits of open-source software such as DeepLabCut (DLC) and R to automate, customize and enhance data analysis of motor behavior. We recorded 2 different spinocerebellar ataxia type 6 mouse models while performing the classic beamwalk test, tracked multiple body parts using the markerless pose-estimation software DLC and analyzed the tracked data using self-written scripts in the programming language R. The beamwalk analysis script (BAS) counts and classifies minor and major hindpaw slips with an 83% accuracy compared to manual scoring. Nose, belly and tail positions relative to the beam, as well as the angle at the tail base relative to the nose and tail tip were determined to characterize motor deficits in greater detail. Our results found distinct ataxic abnormalities such as an increase in major left hindpaw slips and a lower belly and tail position in both SCA6 ataxic mouse models compared to control mice at 18 months of age. Furthermore, a more detailed analysis of various body parts relative to the beam revealed an overall lower body position in the SCA684Q compared to the CT-longQ27PC mouse line at 18 months of age, indicating a more severe ataxic deficit in the SCA684Q group.

From Chingford to Carsington: impact of soil mechanics on dam construction, 1937–1987

In 1937 the investigation of the failure during construction of Chingford number 2 embankment dam brought the concepts of modern soil mechanics to bear on the construction of British embankment dams. Effective stress theory highlighted the crucial significance of pore-water pressure in understanding embankment stability and it became standard practice to monitor pore pressures during embankment construction using twin-tube hydraulic piezometers. By the 1980s, embankment stability might have seemed assured, with embankment design utilising limit equilibrium stability analyses that incorporated strength parameters of fill and foundation materials determined by laboratory testing. However, in 1984 the upstream slope of Carsington dam suffered a major slip in the final stage of construction. Recently developed advanced finite-element analyses were able to demonstrate that progressive failure was a major factor in the slope instability, but the failure was also a reminder of the continuing need for experience-based judgement in design.

Rapid rupture characterization for the 2023 MS 6.2 Jishishan earthquake

On December 18, 2023, the MS 6.2 Jishishan earthquake occurred in the northeastern region of the Qinghai-Xizang Plateau, causing heavy casualties and property damage in Gansu and Qinghai Provinces. In this study, we integrate space imaging geodesy, finite fault inversion, and back-projection methods to decipher its rupture property, including fault geometry, coseismic slip distribution, rupture direction, and propagation speed. The results reveal that the seismogenic fault dips to the southwest at an angle of 29°. The major slip asperity is dominated by reverse slip and is concentrated within a depth range of 7–16 ​km, which explains the significant uplift near the epicenter observed by both the Sentinel-1 ascending and descending InSAR data. Moreover, the teleseismic array waveforms indicate a northwest propagating rupture with an overall slow rupture velocity of ∼1.91 ​km/s (AK array) or 1.01 ​km/s (AU array).

Energy landscape of perylenediimide chromophoric aggregates.

Understanding the self-assembly of conjugated organic materials at the molecular level is crucial in their potential applications as active components in electronic and optoelectronic devices. The type of aggregation significantly influences the intriguing electronic and optical characteristics differing from their constituent molecules. Perylenediimides (PDIs), electron-deficient molecules exhibiting remarkable n-type semiconducting properties, are among the most explored organic fluorescent materials due to their high fluorescence efficiency, photostability, and optoelectronic properties. PDI derivatives are reported to form well-tailored supramolecular architectures: cofacial with minor slip (H-aggregates), staggered with major slip (J-aggregates), magic angle stacking (M-aggregates), rotated (X-aggregates), rotated orthogonal ((+)-aggregates), etc. H*-aggregates are defined here as an ideal case of H-aggregate with an eclipsed configuration. Although numerous reports regarding the formation and optical properties of various PDI aggregates are known, the key driving force within the PDI units guiding the self-assembly to form distinct aggregate systems remains elusive. To unravel the molecular-level mechanisms behind the self-assembly of PDI units by probing the intermolecular interactions, symmetry-adapted perturbation theory-based energy decomposition, potential energy surface scans, and non-covalent interaction index analyses were employed on PDI dimer models. Quantum theory of atoms in molecules and frontier molecular orbital analyses were implemented on the dimer models to comprehend the effect of heteroatoms and orbital interactions in stabilising the X-aggregates over the other PDI aggregate systems. Competition between the attractive and repulsive non-covalent interactions dictates a stability order of X > H > J > M > (+) > H* for the PDI aggregate system, while in the parent perylene system, the stability order was found to be X > (+) > H > M > J > H*.

A Bayesian Source Model for the 2022 Mw6.6 Luding Earthquake, Sichuan Province, China, Constrained by GPS and InSAR Observations

Until the Mw 6.6 Luding earthquake ruptured the Moxi section of the Xianshuihe fault (XSHF) on 5 September 2022, the region had not experienced an Mw >6 earthquake since instrumental records began. We used Global Positioning System (GPS) and Sentinel-1 interferometric synthetic aperture radar (InSAR) observations to image the coseismic deformation and constrain the location and geometry of the seismogenic fault using a Bayesian method We then present a distributed slip model of the 2022 Mw6.6 Luding earthquake, a left-lateral strike-slip earthquake that occurred on the Moxi section of the Xianshuihe fault in the southwest Sichuan basin, China. Two tracks (T26 and T135) of the InSAR data captured a part of the coseismic surface deformation with the line-of-sight displacements range from ∼−0.16 m to ~0.14 m in the ascending track and from ~−0.12 m to ~0.10 m in the descending track. The inverted best-fitting fault model shows a pure sinistral strike-slip motion on a west-dipping fault plane with a strike of 164.3°. We adopt a variational Bayesian approach and account for the uncertainties in the fault geometry to retrieve the distributed slip model. The inverted result shows that the maximum slip of ~1.82 m occurred at a depth of 5.3 km, with the major slip concentrated within depths ranging from 0.9–11 km. The InSAR-determined moment is 1.3 × 1019 Nm, with a shear modulus of 30 GPa, equivalent to Mw 6.7. The published coseismic slip models of the 2022 Luding earthquake show apparent differences despite the use of similar geodetic or seismic observations. These variations underscore the uncertainty associated with routinely performed source inversions and their interpretations for the underlying fault model.

Complex conjugate rupture of the 2014 Mw 6.2 Ludian (Yunnan, China) earthquake

The Mw 6.2 Ludian earthquake, struck Yunnan, China on August 3rd, 2014, leading to a complex aftershock distribution and enormous disaster. The Ludian earthquake is well recorded by the China Earthquake Network. The rupture process of the Ludian earthquake is constrained by the more local strong motion data including the LLT station with epicentral distance of ∼2 km and teleseismic body waves and surface waves in this study. The local strong motion data constrained the epicenter’s location and rupture process better. Waveform fit of key local strong-motion stations supported that the Ludian earthquake ruptured two conjugate faults with strike angles of 77° and 160°, respectively. The rupture simultaneously took place on two conjugate faults and lasted approximately 10 s. The peak slip, rise time, and average slip rate are 0.63 m, 0.69 s, and 2.12 m/s, respectively. The result shows that the main rupture zone was on the fault with a strike of 160°. The depth of the major slip patch was 0–4 km with a centroid depth of 2 km. It indicates that the shallow slip is the main reason for the enormous disaster from the earthquake source perspective.

Exploring the effects of high pressure on hydrogen bonding in pharmaceutical cocrystals: A systematic study of pyridine dicarboxylic acid systems using synchrotron and neutron diffraction

Pharmaceutical cocrystals use common robust hydrogen bonding synthons to create novel materials with different physicochemical properties. In this systematic study of a series of cocrystals, we explore the effect of high pressure on one of these commonly used motifs, the acid-pyridine motif, to assess the commonality of behaviour under extreme conditions. We have surveyed five pyridine dicarboxylic acid systems using both synchrotron and neutron diffraction methods to elucidate the changes in structure. We observe that the hydrogen bonding in these systems compress at a similar rate despite the changes to the molecular make-up of the solids and that on compression the changes in structure are indicative that the layers move along the major slip planes in the structure. We have observed two phase transitions to new forms of the pyrazine:malonic acid system, one for each stoichiometric ratio. This study demonstrates that the combination of two complementary diffraction approaches is key to understanding polymorphic behaviour at high pressure.

Numerical Modeling of Earthquake Cycles Based On Navier‐Stokes Equations With Viscoelastic‐Plasticity Rheology

AbstractVisco‐elastic‐plastic modeling approaches for long‐term tectonic deformation assume that co‐seismic fault displacement can be integrated over 1000s–10,000s years (tens of earthquake cycles) with the appropriate failure law, and that short‐timescale fluctuations in the stress field due to individual earthquakes have no effect on long‐term behavior. Models of the earthquake rupture process generally assume that the tectonic (long‐range) stress field or kinematic boundary conditions are steady over the course of multiple earthquake cycles. This study is aimed to fill the gap between long‐term and short‐term deformations by modeling earthquake cycles with the rate‐and‐state frictional (RSF) relationship in Navier‐Stokes equations. We reproduce benchmarks at the earthquake timescale to demonstrate the effectiveness of our approach. We then discuss how these high‐resolution models degrade if the time‐step cannot capture the rupture process accurately and, from this, infer when it is important to consider coupling of the two timescales and the level of accuracy required. To build upon these benchmarks, we undertake a generic study of a thrust fault in the crust with a prescribed geometry. It is found that lower crustal rheology affects the periodic time of characteristic earthquake cycles and the inter‐seismic, free‐surface deformation rate. In particular, the relaxation of the surface of a cratonic region (with a relatively strong lower crust) has a characteristic double‐peaked uplift profile that persists for thousands of years after a major slip event. This pattern might be diagnostic of active faults in cratonic regions.

SLIP SOURCE MODEL OF THE 1995 NEFTEGORSK EARTHQUAKE (NORTH SAKHALIN) FROM GEODETIC DATA

The May 27, 1995 Mw=7.0 Neftegorsk earthquake occurred in the north of Sakhalin Island, rupturing the Upper Piltun fault, a secondary feature of the main Hokkaido-Sakhalin regional fault zone. The fault geometry, coseismic slip model, and Coulomb stress changes in the earthquake focal area were calculated based on a finite fault modeling. We used near-field coseismic offsets at 24 points obtained by comparison between predating triangulation and GPS observations, which were collected before and after the earthquake. Our slip distribution model shows two major slip patches. Larger slip asperity (amplitude up to 6.36 m) was characterized by right-lateral strike-slip movements, which correspond to focal mechanism of the earthquake, whereas the northern segment has reverse fault mechanism with maximum slip of 2.64 m. The fault length and width, average slip and stress drop values are estimated at 78 km, 28 km, 1.91 m and 11.3 MPa, respectively. The estimated release moment is approximately 7.49×1019 N∙m equal to Mw=7.2, which is larger than that reported by the USGS and GCMT but consistent with the values reported by other researchers. The coseismic Coulomb stress changes enhanced the stress by more than 10 MPa on the southern segment of the Gyrgylaninsky fault and middle section of the Hokkaido–Sakhalin fault. Seismic risks on the nearest faults cannot be ignored in the future despite the fact that the earthquake with a magnitude of 5.8 occurred in 2010 near the Gyrgylaninsky fault. The recent GPS rates in the surroundings of the Neftegorsk surface rupture mean that the recurrence interval for similar earthquakes may be more than a thousand years.

Re‐Examining Temporal Variations in Intermediate‐Depth Seismicity

AbstractChanges in the frequency of intermediate‐depth (60–300 km) earthquakes in response to static stress transfer can provide insights into the mechanisms of earthquake generation within subducting slabs. In this study, we use the most up‐to‐date global and regional earthquake catalogs to show that both the aftershock productivity of large earthquakes, and the changes in the frequency of intermediate‐depth earthquakes around the timing of major megathrust slip, support the view that faults within the slab are relatively insensitive to static stress transfer on the order of earthquake stress drops. We interpret these results to suggest the population of faults within the slab are much further from their failure stress than is typical for shallow fault systems. We also find that aftershock productivity varies within slabs over small spatial scales, indicating that the mechanism that enables faults to rupture at intermediate depths is likely to be spatially heterogeneous over length‐scales of a few tens of kilometres. We suggest dehydration‐related weakening mechanisms can best account for this heterogeneity.

The 2021 Ms 6.0 Luxian (China) Earthquake: Blind Reverse‐Fault Rupture in Deep Sedimentary Formations Likely Induced by Pressure Perturbation From Hydraulic Fracturing

AbstractTo resolve the occurrence mechanism of the Luxian earthquake occurred near an active hydraulic fracturing well pad on 15 September 2021, we first use seismic waveforms to invert the focal mechanism solution and centroid depth, and then utilize Sentinel‐1 Synthetic Aperture Radar images and Global Navigation Satellite System observations to determine the seismogenic fault and slip distribution. Our results show that the Luxian event ruptured a previously‐unmapped southwest‐dipping reverse fault intersecting with one horizontal well of the well pad. Major slip occurred below the shale gas bed (∼4 km deep) but above the crystalline basement (∼7 km deep). Further analyses on preseismic surface deformation initiated in the northeast of the event reveal that pore pressure near the hypocenter was increased by ∼4.5 MPa. Consequently, fracking fluid injected through the horizontal well and other wells of nearby pads likely flowed directly into the fault zones and prompted the occurrence of the Luxian earthquake.

Hydromagnetic flow of magnetite–water nanofluid utilizing adapted Buongiorno model

The hydromagnetic flow of magnetite–water nanofluid due to a rotating stretchable disk has been numerically assessed. The nanofluid flow has been modeled utilizing the adapted Buongiorno model that considers the volume fraction-dependent effective nanofluid properties and the major slip mechanisms. In addition, experimentally gleaned functions of effective dynamic viscosity and effective thermal conductivity are deployed. The modeled equations are transformed into a first-order ODEs scheme employing Von Kármán’s similarity conversions and then resolved via the Runge–Kutta algorithm through the shooting technique. The impact of pertinent terms over the physical quantities, nanoliquid temperature and nanoliquid concentration is explained with the support of graphs. Results show that rising volume fraction of magnetite nanoparticles (NPs) and magnetic field term enhance the drag force. Mass transport rate is demoted with augmenting values of magnetic field parameter whereas is promoted with increase in Schmidt number. Further, it is detected that the changes in stretching strength parameter are directly proportional to Nusselt number and inversely proportional to the thermal field. The findings of this numerical analysis have applications in spin coating, rotating disk reactors, storage devices for computers, food processing, and rotating heat exchangers.

Numerical Study on the Heat Transfer Characteristics of Cu-Water and TiO2-Water Nanofluid in a Circular Horizontal Tube

A numerical simulation of convective heat transfer coefficient (hconv) was studied with Cu-Water and TiO2-Water nanofluids flowing through a circular tube subjected to uniform wall heat flux boundary conditions under laminar and turbulent regimes. Four different concentrations of nanofluids (ɸ = 0.5, 1, 1.5 and 2%) were used for the analysis and the Reynolds number (Re) was varied from laminar (500 to 2000) to turbulent flow regime (5000 to 20,000). The dependence of hconv on Re and ɸ was investigated using a single-phase Newtonian approach. In comparison to base fluid, average hconv enhancements of 10.4% and 7.3% were noted, respectively, for the maximum concentration (ɸ = 2%) and Re = 2000 for Cu-Water and TiO2—water nanofluids in the laminar regime. For the same ɸ under the turbulent regime (Re = 20,000), the enhancements were noted to be 14.6% and 13.2% for both the nanofluids, respectively. The random motion (Brownian motion) and heat diffusion (thermophoresis) by nanosized particles are the two major slip mechanisms that have more influence on the enhancement of hconv. In addition, the Nusselt number (Nu) of the present work was validated for water with the Shah and Dittus Boelter equation and found to have good agreement for both the regimes.

The 2021 Loyalty Islands Earthquake (Mw 7.7): Tsunami Waveform Inversion and Implications for Tsunami Forecasting for New Zealand

AbstractA tsunamigenic earthquake with thrust faulting mechanism occurred southeast of the Loyalty Islands, New Caledonia, in the Southern Vanuatu subduction zone on the 10th of February 2021. The tsunami was observed at coastal gauges in the surrounding islands and in New Zealand. The tsunami was also recorded at a new DART network designed to enhance the tsunami forecasting capability of the Southwestern Pacific. We used the tsunami waveforms in an inversion to estimate the fault slip distribution. The estimated major slip region is located near the trench with maximum slip of 4 m. This source model with an assumed rupture velocity of 1.0 km/s can reproduce the observed seismic waves. We evaluated two tsunami forecasting approaches for coastal regions in New Zealand: selecting a pre‐computed scenario, and interpolating between two pre‐computed scenarios. For the evaluation, we made a reference map of tsunami threat levels in New Zealand using the estimated source model. The results show that the threat level maps from the pre‐computed Mw 7.7 scenario located closest to the epicenter, and from an interpolation of two scenarios, match the reference threat levels in most coastal regions. Further improvements to enhance the system toward more robust warnings include expansion of scenario database and incorporation of tsunami observation around tsunami source regions. We also report on utilization of the coastal gauge and DART station data for updating forecasts in real‐time during the event and discuss the differences between the rapid‐response forecast and post‐event retrospective forecasts.

Continuous estimation of coseismic and early postseismic slip phenomena via the GNSS carrier phase to fault slip approach: a case study of the 2011 Tohoku-Oki sequence

The transition process from coseismic to early postseismic phenomena within a half-day remains a significant topic for understanding the slip budget and friction properties of the fault. However, the investigation of this phase has undergone limited advancement. This is mainly due to the lack of precision pertaining to the Global Navigation Satellite System (GNSS) analysis caused by difficulty of separation between fault slip and other unknown parameters such as tropospheric delay. Therefore, we propose an alternative approach (phase-to-slip; PTS) that detects fault slip directly from GNSS carrier phase variation without conventional positioning. Since the PTS simultaneously estimates fault slip and other unknowns, we can quantitatively evaluate their separation accuracy and contribution in the estimation. This study attempts to continuously estimate the coseismic and early postseismic slip of the 2011 Tohoku-Oki sequence using the PTS method. We analyzed 1-Hz carrier phase data for approximately 2 h before and after the mainshock origin (2011/03/11/05:46UTC). As a result, we successfully obtained the coseismic slip of the Mw9.0 mainshock and two major aftershocks in the off-Ibaraki (Mw7.8) and off Iwate (Mw7.4) regions. For all three events, the estimated slip distribution, equivalent moment magnitude, and calculated displacement field well agreed with those from the conventional positioning. Additionally, our results suggest postseismic slip mainly in the downdip area adjacent to the mainshock rupture. We obtained three major slip areas around the downdip region near Aomori, Iwate, and Miyagi. The estimated slip amount reaches approximately 0.1–0.2 m in 34 min immediately after the mainshock. The equivalent seismic moment of the slip areas near Iwate and Miyagi were approximately Mw7.4. This amount is similar or slightly greater than the estimations based on conventional positioning. Significant overlap is observed in the locations of the slip areas. In this manner, PTS continuously detected a series of coseismic and early postseismic slips with timescale ranging from a few minutes to an hour. Our results demonstrate the capability of the PTS method for broadband fault slip monitoring and its potential contribution to the investigation of early postseismic phenomena.

Present-day crustal strain and major fault slip rates in North China determined using GNSS observations

• A strain model and an elastic block model of North China were established. • The western Ordos and the northern North China exhibited larger strain rates. • The NNE- and NWW-trending faults showed dextral and sinistral slip respectively. • The fault deformation in North China was controlled by neighboring blocks. North China is characterized by numerous faults, where many earthquakes were recorded. However, the fault slip and crustal deformation mechanism underlying the effects of neighboring block motions remain poorly understood. In this study, we first constructed a strain model using global navigation satellite system (GNSS) velocities to analyze the relationship between seismicity and strain. We then constructed a block model based on the fault distribution and GNSS observations to obtain fault slip. The strain and fault slip rates were then used to determine the crustal deformation mechanism in North China. The strain model suggests the western border of the Ordos Block and the northern North China Plain had larger maximum shear strain rates. These areas are located in a transition between positive and negative dilatation rates, which indicate intense crustal deformation. Based on the block model, we found the NNE-trending faults were dominated by dextral slip. The northern and southern borders of the Ordos Block indicate left-lateral motion (0.4–2.5 mm/a), while its western and eastern margins exhibited right-lateral slip (1.0–3.0 mm/a). The periphery of the Ordos Block exhibited normal slip (0.2–3.2 mm/a), except at the southwestern tip. These results suggest that the eastward extrusion of the Tibetan Plateau caused left-lateral shear in North China, which led to sinistral slip of NNE-trending faults. The northward expansion of the northeastern Tibetan Plateau also contributes to sinistral slip along the northern and southern margins of the Ordos Block, as well as dextral slip along its western and eastern margins.

Late Pleistocene Boulaaba travertine and calcareous tufa, Kasserine, Central Tunisia: Implications for North African spring-fed Ca-carbonates

Late Pleistocene (≈21000 BP) travertine and calcareous tufa deposits have been recognized in the region of Boulaaba, Kasserine in Central Tunisia with a variety of morphological aspects and depositional facies. Morphological aspects of these deposits consist of two elongated continental Ca-carbonate masses defined as fissure-ridges built-up around old fissure networks. Distinction between travertine and tufa lithofacies, based on textural, compositional and topographic settings regarding the entire carbonate masses, have led for travertine to 5 lithofacies: (1) banded travertine fissure infill, (2) bacterial shrubs, (3) stromatolitic laminae, (4) pisoids and (5) speleothem-like travertine growth habits. Tufa is composed of autochtonous tufa lithofacies with (1) phytoherme framestone tufa, (2) calcified insect larval constructions, and clastic alllochtonous tufa lithofacies with (1) matrix-supported microdetrital tufa and (2) grain-supported macrodetrital phytoclast tufa. Both lithofacies show variable diagenetic evolutions with an overall lithification. These deposits are concentrated along the Foussana trough-related normal faults and fractures that belong to the extensional regime of the major strike slip fault of Kasserine that enhanced and facilitated fluid upwelling along with movements of faults which damage and fracture the rocks and increase the overall permeability. During this phase (Late Pleistocene), the deposition of Boulaaba travertine and tufa has implied some paleo-hydroclimatic conditions marked in most North African travertine and tufa depositing sites by (1) a higher water table with abundant effective rainfall on the recharge areas, (2) dense plant cover and soils giving acidic water infiltrations, (3) karstic dissolution of Cretaceous limestones (especially the limestone aquifers of Bireno and Abiod), and (4) a strong water mineralization of the groundwater.

Could thermal pressurization have induced the frequency-dependent rupture during the 2019 Mw8.0 Peru intermediate-depth earthquake?

SUMMARY The rupture process of earthquakes at intermediate-depth (∼70–300 km) have rarely been illuminated by a joint analysis of geodetic and seismic observations, hindering our understanding of their dynamic rupture mechanisms. Here we present detailed rupture process of the 2019 Mw8.0 Peru earthquake at the depth of 122 km depth, derived with a holistic approach reconciling InSAR and broad-band seismic waveform data. The joint inversion of InSAR observations and teleseismic body waves results in a finite rupture model that extends ∼200 km along strike, with unilateral rupture towards north that lasted for ∼60 s. There are four major slip patches in the finite fault model which are well corresponding to the position and timing of the sources in back-projection and multiple points source results. The largest slip patch, which occurred ∼40 s after the rupture initiation, had a longer and smoother rise time, and radiated much weaker high-frequency seismic waves compared to other smaller slip patches. This distinct frequency-dependent rupture could be explained by a strong dynamic weakening mechanism. We question whether thermal pressurization of pore free water rather than thermal run away could be such a mechanism. Our frequency content analysis could be generalized to study other earthquakes including those deeper than 300 km.