Search for: All records

o what extent mechanical anisotropy is required to explain the dynamics of the lithosphere is an important yet unresolved question. If anisotropy affects stress and deformation, and hence processes such as fault loading, how can we quantify its role from observations? Here, we derive analytical solutions and build a theoretical framework to explore how a shear zone with linear anisotropic viscosity can lead to deviatoric stress heterogeneity, strain-rate enhancement, as well as non-coaxial principal stress and strain rate. We develop an open-source finite-element software based on FEniCS for more complicated scenarios in both 2-D and 3-D. Mechanics of shear zones with transversely isotropic and orthorhombic anisotropy subjected to misoriented shortening and simple shearing are explored. A simple regional example for potential non-coaxiality for the Leech River Schist above the Cascadia subduction zone is presented. Our findings and these tools may help to better understand, detect and evaluate mechanical anisotropy in natural settings, with potential implications including the transfer of lithospheric stress and deformation through fault loading. 

El Niño events, the warm phase of the El Niño–Southern Oscillation (ENSO) phenomenon, amplify climate variability throughout the world. Uncertain climate model predictions limit our ability to assess whether these climatic events could become more extreme under anthropogenic greenhouse warming. Palaeoclimate records provide estimates of past changes, but it is unclear if they can constrain mechanisms underlying future predictions. Here we uncover a mechanism using numerical simulations that drives consistent changes in response to past and future forcings, allowing model validation against palaeoclimate data. The simulated mechanism is consistent with the dynamics of observed extreme El Niño events, which develop when western Pacific warm pool waters expand rapidly eastwards because of strongly coupled ocean currents and winds. These coupled interactions weaken under glacial conditions because of a deeper mixed layer driven by a stronger Walker circulation. The resulting decrease in ENSO variability and extreme El Niño occurrence is supported by a series of tropical Pacific palaeoceanographic records showing reduced glacial temperature variability within key ENSO-sensitive oceanic regions, including new data from the central equatorial Pacific. The model–data agreement on past variability, together with the consistent mechanism across climatic states, supports the prediction of a shallower mixed layer and weaker Walker circulation driving more frequent extreme El Niño genesis under greenhouse warming. 

Rock strength has long been linked to lithospheric deformation and seismicity. However, independent constraints on the related elastic heterogeneity are missing, yet could provide key information for solid Earth dynamics. Using coseismic Global Navigation Satellite Systems (GNSS) data for the 2011 M9 Tohoku-oki earthquake in Japan, we apply an inverse method to infer elastic structure and fault slip simultaneously. We find compliant material beneath the volcanic arc and in the mantle wedge within the partial melt generation zone inferred to lie above ~100 km slab depth. We also identify low-rigidity material closer to the trench matching seismicity patterns, likely associated with accretionary wedge structure. Along with traditional seismic and electromagnetic methods, our approach opens up avenues for multiphysics inversions. Those have the potential to advance earthquake and volcano science, and in particular once expanded to InSAR type constraints, may lead to a better understanding of transient lithospheric deformation across scales. 

Abstract Understanding mechanical conditions that lead to complexity in earthquakes is important to seismic hazard analysis. In this study, we simulate physics‐based multicycle dynamic models of the San Andreas fault (Carrizo through San Bernardino sections) and the San Jacinto fault (Claremont and Clark strands). We focus on a complex fault geometry based on the Southern California Earthquake Center Community Fault Model and its effect over multiple earthquake cycles. Using geodetically derived strain rates, we validate the models against geologic slip rates and recurrence intervals at various paleoseismic sites. We find that the interactions among fault geometry, dynamic rupture and interseismic stress accumulation produce stress heterogeneities, leading to rupture segmentation and variability in earthquake recurrence. Our models produce earthquakes with rupture extents similar to a recent comprehensive paleoseismic catalog. The “earthquake gates” of the Big Bend and the Cajon Pass occasionally impede dynamic ruptures. The angle of compression, which is the subtraction of the maximum shear strain rate direction from the local fault strike, can better determine the likelihood of the impedance of restraining bends to dynamic ruptures. Because the Big Bend has an angle of compression of ∼20°, ruptures that traverse the Big Bend, like the 1857 Fort Tejon earthquake, are more frequent than expected based on empirical relations which predict the ∼40° restraining bend to terminate most ruptures. Our models indicate that large ruptures tend to initiate north of the Big Bend and propagate southwards, similar to the 1857 earthquake, providing critical information for ground shaking assessment in the region. 

Abstract Dynamic modeling of sequences of earthquakes and aseismic slip (SEAS) provides a self‐consistent, physics‐based framework to connect, interpret, and predict diverse geophysical observations across spatial and temporal scales. Amid growing applications of SEAS models, numerical code verification is essential to ensure reliable simulation results but is often infeasible due to the lack of analytical solutions. Here, we develop two benchmarks for three‐dimensional (3D) SEAS problems to compare and verify numerical codes based on boundary‐element, finite‐element, and finite‐difference methods, in a community initiative. Our benchmarks consider a planar vertical strike‐slip fault obeying a rate‐ and state‐dependent friction law, in a 3D homogeneous, linear elastic whole‐space or half‐space, where spontaneous earthquakes and slow slip arise due to tectonic‐like loading. We use a suite of quasi‐dynamic simulations from 10 modeling groups to assess the agreement during all phases of multiple seismic cycles. We find excellent quantitative agreement among simulated outputs for sufficiently large model domains and small grid spacings. However, discrepancies in rupture fronts of the initial event are influenced by the free surface and various computational factors. The recurrence intervals and nucleation phase of later earthquakes are particularly sensitive to numerical resolution and domain‐size‐dependent loading. Despite such variability, key properties of individual earthquakes, including rupture style, duration, total slip, peak slip rate, and stress drop, are comparable among even marginally resolved simulations. Our benchmark efforts offer a community‐based example to improve numerical simulations and reveal sensitivities of model observables, which are important for advancing SEAS models to better understand earthquake system dynamics.