Evidence for on-going inflation of the Socorro Magma Body, New Mexico, from interferometric synthetic aperture radar imaging

Interferometric synthetic aperture radar (InSAR) imaging of the central Rio Grande rift (New Mexico, USA) during 1992–1999 reveals a crustal uplift of several centimeters that spatially coincides with the seismologically determined outline of the Socorro magma body, one of the largest currently active magma intrusions in the Earth's continental crust. Modeling of interferograms shows that the observed deformation may be due to elastic opening of a sill-like intrusion at a rate of a few millimeters per year. Despite an apparent constancy of the geodetically determined uplift rate, thermodynamic arguments suggest that it is unlikely that the Socorro magma body has formed via steady state elastic inflation

A unified continuum representation of post-seismic relaxation mechanisms: semi-analytic models of afterslip, poroelastic rebound and viscoelastic flow

We present a unified continuum mechanics representation of the mechanisms believed to be commonly involved in post-seismic transients such as viscoelasticity, fault creep and poroelasticity. The time-dependent relaxation that follows an earthquake, or any other static stress perturbation, is considered in a framework of a generalized viscoelastoplastic rheology whereby some inelastic strain relaxes a physical quantity in the material. The relaxed quantity is the deviatoric stress in case of viscoelastic relaxation, the shear stress in case of creep on a fault plane and the trace of the stress tensor in case of poroelastic rebound. In this framework, the instantaneous velocity field satisfies the linear inhomogeneous Navier's equation with sources parametrized as equivalent body forces and surface tractions. We evaluate the velocity field using the Fourier-domain Green's function for an elastic half-space with surface buoyancy boundary condition. The accuracy of the proposed method is demonstrated by comparisons with finite-element simulations of viscoelastic relaxation following strike-slip and dip-slip ruptures for linear and power-law rheologies. We also present comparisons with analytic solutions for afterslip driven by coseismic stress changes. Finally, we demonstrate that the proposed method can be used to model time-dependent poroelastic rebound by adopting a viscoelastic rheology with bulk viscosity and work hardening. The proposed method allows one to model post-seismic transients that involve multiple mechanisms (afterslip, poroelastic rebound, ductile flow) with an account for the effects of gravity, non-linear rheologies and arbitrary spatial variations in inelastic properties of rocks (e.g. the effective viscosity, rate-and-state frictional parameters and poroelastic properties)

The complete (3-D) surface displacement field in the epicentral area of the 1999 M_W7.1 Hector Mine Earthquake, California, from space geodetic observations

We use Interferometric Synthetic Aperture Radar (InSAR) data to derive continuous maps for three orthogonal components of the co-seismic surface displacement field due to the 1999 M_w7.1 Hector Mine earthquake in southern California. Vertical and horizontal displacements are both predominantly antisymmetric with respect to the fault plane, consistent with predictions of linear elastic models of deformation for a strike-slip fault. Some deviations from symmetry apparent in the surface displacement data may result from complexity in the fault geometry

Coseismic Deformation from the 1999 M_w 7.1 Hector Mine, California, Earthquake as Inferred from InSAR and GPS Observations

We use interferometric synthetic aperture radar (InSAR) and Global Positioning System (GPS) observations to investigate static deformation due to the 1999 M_w 7.1 Hector Mine earthquake, that occurred in the eastern California shear zone. Interferometric decorrelation, phase, and azimuth offset measurements indicate regions of surface and near-surface slip, which we use to constrain the geometry of surface rupture. The inferred geometry is spatially complex, with multiple strands. The southern third of the rupture zone consists of three subparallel segments extending about 20 km in length in a N45°W direction. The central segment is the simplest, with a single strand crossing the Bullion Mountains and a strike of N10°W. The northern third of the rupture zone is characterized by multiple splays, with directions subparallel to strikes in the southern and central. The average strike for the entire rupture is about N30°W. The interferograms indicate significant along-strike variations in strain which are consistent with variations in the ground-based slip measurements. Using a variable resolution data sampling routine to reduce the computational burden, we invert the InSAR and GPS data for the fault geometry and distribution of slip. We compare results from assuming an elastic half-space and a layered elastic space. Results from these two elastic models are similar, although the layered-space model predicts more slip at depth than does the half-space model. The layered model predicts a maximum coseismic slip of more than 5 m at a depth of 3 to 6 km. Contrary to preliminary reports, the northern part of the Hector Mine rupture accommodates the maximum slip. Our model predictions for the surface fault offset and total seismic moment agree with both field mapping results and recent seismic models. The inferred shallow slip deficit is enigmatic and may suggest that distributed inelastic yielding occurred in the uppermost few kilometers of the crust during or soon after the earthquake

Finite source modelling of magmatic unrest in Socorro, New Mexico, and Long Valley, California

We investigate surface deformation associated with currently active crustal magma bodies in Socorro, New Mexico, and Long Valley, California, USA. We invert available geodetic data from these locations to constrain the overall geometry and dynamics of the inferred deformation sources at depth. Our best-fitting model for the Socorro magma body is a sill with a depth of 19 km, an effective diameter of 70 km and a rate of increase in the excess magma pressure of 0.6 kPa yr^(−1). We show that the corresponding volumetric inflation rate is ∼6×10^(−3) km^3 yr^(−1), which is considerably less than previously suggested. The measured inflation rate of the Socorro magma body may result from a steady influx of magma from a deep source, or a volume increase associated with melting of the magma chamber roof (i.e. crustal anatexis). In the latter case, the most recent major injection of mantle-derived melts into the middle crust beneath Socorro may have occurred within the last several tens to several hundreds of years. The Synthetic Interferometric Aperture Radar (InSAR) data collected in the area of the Long Valley caldera, CA, between June 1996 and July 1998 reveal an intracaldera uplift with a maximum amplitude of ∼11 cm and a volume of 3.5×10^(−2) km^3. Modelling of the InSAR data suggests that the observed deformation might be due to either a sill-like magma body at a depth of ∼12 km or a pluton-like magma body at a depth of ∼8 km beneath the resurgent dome. Assuming that the caldera fill deforms as an isotropic linear elastic solid, a joint inversion of the InSAR data and two-colour laser geodimeter data (which provide independent constraints on horizontal displacements at the surface) suggests that the inferred magma chamber is a steeply dipping prolate spheroid with a depth of 7–9 km and an aspect ratio in excess of 2:1. Our results highlight the need for large radar look angles and multiple look directions in future InSAR missions

Estimate of differential stress in the upper crust from variations in topography and strike along the San Andreas fault

The major bends of the San Andreas fault in California are associated with significant variations in the along-fault topography. The topography-induced perturbations in the intermediate principal stress may result in the rotation of the fault with respect to the maximum compression axis provided that the fault is non-vertical, and the slip is horizontal. The progressive fault rotation may produce additional topography via thrust faulting in the adjacent crust, resulting in a positive feedback. The observed rotation of the fault plane due to the along-fault variations in topography is used to infer the magnitude of the in situ differential stress. Our results suggest that the average differential stress in the upper crust around the San Andreas fault is of the order of 50 MPa, implying that the effective fault strength is about a factor of two lower than predictions based on Byerlee's law and the assumption of hydrostatic pore pressure

Deformation due to a pressurized horizontal circular crack in an elastic half-space, with applications to volcano geodesy

We consider deformation due to sill-like magma intrusions using a model of a horizontal circular crack in a semi-infinite elastic solid. We present exact expressions for vertical and horizontal displacements of the free surface of a half-space, and calculate surface displacements for a special case of a uniformly pressurized crack. We derive expressions for other observable geophysical parameters, such as the volume of a surface uplift/subsidence, and the corresponding volume change due to fluid injection/withdrawal at depth. We demonstrate that for essentially oblate (i.e. sill-like) source geometries the volume change at the source always equals the volume of the displaced material at the surface of a half-space. Our solutions compare favourably to a number of previously published approximate models. Surface deformation due to a ‘point’ crack (that is, a crack with a large depth-to-radius ratio) differs appreciably from that due to an isotropic point source (‘Mogi model’). Geodetic inversions that employ only one component of deformation (either vertical or horizontal) are unlikely to resolve the overall geometry of subsurface deformation sources even in a simplest case of axisymmetric deformation. Measurements of a complete vector displacement field at the Earth's surface may help to constrain the depth and morphology of active magma reservoirs. However, our results indicate that differences in surface displacements due to various axisymmetric sources may be subtle. In particular, the sill-like and pluton-like magma chambers may give rise to differences in the ratio of maximum horizontal displacements to maximum vertical displacements (a parameter that is most indicative of the source geometry) that are less than 30 per cent. Given measurement errors in geodetic data, such differences may be hard to distinguish

Deformation on Nearby Faults Induced by the 1999 Hector Mine Earthquake

Interferometric Synthetic Aperture Radar observations of surface deformation due to the 1999 Hector Mine earthquake reveal motion on several nearby faults of the eastern California shear zone.We document both vertical and horizontal displacements of several millimeters to several centimeters across kilometerwide zones centered on pre-existing faults.Portions of some faults experienced retrograde (that is, opposite to their long-term geologic slip) motion during or shortly after the earthquake.The observed deformation likely represents elastic response of compliant fault zones to the permanent co-seismic stress changes. The induced fault displacements imply decreases in the effective shear modulus within the kilometer-wide fault zones, indicating that the latter are mechanically distinct from the ambient crustal rocks