Long-period ocean sound waves constrain shallow slip and tsunamis in megathrust ruptures

Great earthquakes along subduction-zone plate boundaries, like the magnitude 9.0 Tohoku-Oki, Japan, event, deform the seafloor to generate massive tsunamis. Tsunami wave heights near shore are greatest when excitation occurs far offshore near the trench, where water depths are greatest and fault slip is shallow. Unfortunately the rupture process there is poorly constrained with land-based geodetic and even seafloor deformation measurements. Here we demonstrate, through dynamic rupture simulations of the Tohoku event, that long-period sound waves in the ocean, observable with ocean-bottom pressure sensors and/or seismometers, can resolve the shallow rupture process and tsunami excitation near the trench. These waves could potentially be used to improve local tsunami early warning systems

Earthquake Slip Between Dissimilar Poroelastic Materials

A mismatch of elastic properties across a fault induces normal stress changes during spatially nonuniform in-plane slip. Recently, Rudnicki and Rice showed that similar effects follow from a mismatch of poroelastic properties (e.g., permeability) within fluid-saturated fringes of damaged material along the fault walls; in this case, it is pore pressure on the slip plane and hence effective normal stress that is altered during slip. The sign of both changes can be either positive or negative, and they need not agree. Both signs reverse when rupture propagates in the opposite direction. When both elastic and poroelastic properties are discontinuous across the fault, steady sliding at a constant friction coefficient, f, is unstable for arbitrarily small f if the elastic mismatch permits the existence of a generalized Rayleigh wave. Spontaneous earthquake rupture simulations on regularized slip-weakening faults confirm that the two effects have comparable magnitudes and that the sign of the effective normal stress change cannot always be predicted solely from the contrast in elastic properties across the fault. For opposing effects, the sign of effective normal stress change reverses from that predicted by the poroelastic mismatch to that predicted by the elastic mismatch as the rupture accelerates, provided that the wave speed contrast exceeds about 5–10% (the precise value depends on the poroelastic contrast and Skempton's coefficient). For faults separating more elastically similar materials, there exists a minimum poroelastic contrast above which the poroelastic effect always determines the sign of the effective normal stress change, no matter the rupture speed.Earth and Planetary SciencesEngineering and Applied Science

Adjoint-based inversion for stress and frictional parameters in earthquake modeling

We present an adjoint-based optimization method to invert for stress and frictional parameters used in earthquake modeling. The forward problem is linear elastodynamics with nonlinear rate-and-state frictional faults. The misfit functional quantifies the difference between simulated and measured particle displacements or velocities at receiver locations. The misfit may include windowing or filtering operators. We derive the corresponding adjoint problem, which is linear elasticity with linearized rate-and-state friction with time-dependent coefficients derived from the forward solution. The gradient of the misfit is efficiently computed by convolving forward and adjoint variables on the fault. The method thus extends the framework of full-waveform inversion to include frictional faults with rate-and-state friction. In addition, we present a space-time dual-consistent discretization of a dynamic rupture problem with a rough fault in antiplane shear, using high-order accurate summation-by-parts finite differences in combination with explicit Runge--Kutta time integration. The dual consistency of the discretization ensures that the discrete adjoint-based gradient is the exact gradient of the discrete misfit functional as well as a consistent approximation of the continuous gradient. Our theoretical results are corroborated by inversions with synthetic data. We anticipate that adjoint-based inversion of seismic and/or geodetic data will be a powerful tool for studying earthquake source processes; it can also be used to interpret laboratory friction experiments.Comment: Updated title, added additional references, provided additional details in sections 1 and 5, fixed typo

Earthquake ruptures with thermal weakening and the operation of major faults at low overall stress levels

We model ruptures on faults that weaken in response to flash heating of microscopic asperity contacts (within a rate-and-state framework) and thermal pressurization of pore fluid. These are arguably the primary weakening mechanisms on mature faults at coseismic slip rates, at least prior to large slip accumulation. Ruptures on strongly rate-weakening faults take the form of slip pulses or cracks, depending on the background stress. Self-sustaining slip pulses exist within a narrow range of stresses: below this range, artificially nucleated ruptures arrest; above this range, ruptures are crack-like. Natural earthquakes will occur as slip pulses if faults operate at the minimum stress required for propagation. Using laboratory-based flash heating parameters, propagation is permitted when the ratio of shear to effective normal stress on the fault is 0.2–0.3; this is mildly influenced by reasonable choices of hydrothermal properties. The San Andreas and other major faults are thought to operate at such stress levels. While the overall stress level is quite small, the peak stress at the rupture front is consistent with static friction coefficients of 0.6–0.9. Growing slip pulses have stress drops of ∼3 MPa; slip and the length of the slip pulse increase linearly with propagation distance at ∼0.14 and ∼30 m/km, respectively. These values are consistent with seismic and geologic observations. In contrast, cracks on faults of the same rheology have stress drops exceeding 20 MPa, and slip at the hypocenter increases with distance at ∼1 m/km

Thermo- and Hydro-mechanical Processes along Faults during Rapid Slip

Field observations of maturely slipped faults show a generally broad zone of damage by cracking and granulation. Nevertheless, large shear deformation, and therefore heat generation, in individual earthquakes takes place with extreme localization to a zone <1–5 mm wide within a finely granulated fault core. Relevant fault weakening processes during large crustal events are therefore likely to be thermal. Further, given the porosity of the damage zones, it seems reasonable to assume groundwater presence. It is suggested that the two primary dynamic weak- ening mechanisms during seismic slip, both of which are expected to be active in at least the early phases of nearly all crustal events, are then as follows: (1) Flash heating at highly stressed frictional micro-contacts, and (2) Thermal pressurization of fault-zone pore fluid. Both have characteristics which promote extreme localization of shear. Macroscopic fault melting will occur only in cases for which those processes, or others which may sometimes become active at large enough slip (e.g., thermal decomposition, silica gelation), have not sufficiently reduced heat generation and thus limited temperature rise. Spontaneous dynamic rupture modeling, using procedures that embody mechanisms (1) and (2), shows how faults can be statically strong yet dynamically weak, and oper- ate under low overall driving stress, in a manner that generates negligible heat and meets major seismic constraints on slip, stress drop, and self-healing rupture mode.Earth and Planetary SciencesEngineering and Applied Science

A Finite Difference Method for Off-fault Plasticity throughout the Earthquake Cycle

We have developed an efficient computational framework for simulating multiple earthquake cycles with off-fault plasticity. The method is developed for the classical antiplane problem of a vertical strike-slip fault governed by rate-and-state friction, with inertial effects captured through the radiationdamping approximation. Both rate-independent plasticity and viscoplasticity are considered, where stresses are constrained by a Drucker-Prager yield condition. The off-fault volume is discretized using finite differences and tectonic loading is imposed by displacing the remote side boundaries at a constant rate. Time-stepping combines an adaptive Runge-Kutta method with an incremental solution process which makes use of an elastoplastic tangent stiffness tensor and the return-mapping algorithm. Solutions are verified by convergence tests and comparison to a finite element solution. We quantify how viscosity, isotropic hardening, and cohesion affect the magnitude and off-fault extent of plastic strain that develops over many ruptures. If hardening is included, plastic strain saturates after the first event and the response during subsequent ruptures is effectively elastic. For viscoplasticity without hardening, however, successive ruptures continue to generate additional plastic strain. In all cases, coseismic slip in the shallow sub-surface is diminished compared to slip accumulated at depth during interseismic loading. The evolution of this slip deficit with each subsequent event, however, is dictated by the plasticity model. Integration of the off-fault plastic strain from the viscoplastic model reveals that a significant amount of tectonic off-set is accommodated by inelastic deformation (~0.1 m per rupture, or ~10% of the tectonic deformation budget)

"Внутренняя картина болезни" и ее динамика у больных желчнокаменной болезнью

ЖЕЛЧНОКАМЕННАЯ БОЛЕЗНЬВНУТРЕННИЕ БОЛЕЗНИАСТЕНИЯКАЧЕСТВО ЖИЗНИПСИХОТЕРАПИ