Eric Hetland

The 2011 Magnitude 9.0 Tohoku-Oki Earthquake: Mosaicking the Megathrust from Seconds to Centuries

Detailed geophysical measurements reveal features of the 2011 Tohoku-Oki megathrust earthquake. Geophysical observations from the 2011 moment magnitude (Mw) 9.0 Tohoku-Oki, Japan earthquake allow exploration of a rare large event along a subduction megathrust. Models for this event indicate that the distribution of coseismic fault slip exceeded 50 meters in places. Sources of high-frequency seismic waves delineate the edges of the deepest portions of coseismic slip and do not simply correlate with the locations of peak slip. Relative to the Mw 8.8 2010 Maule, Chile earthquake, the Tohoku-Oki earthquake was deficient in high-frequency seismic radiation—a difference that we attribute to its relatively shallow depth. Estimates of total fault slip and surface secular strain accumulation on millennial time scales suggest the need to consider the potential for a future large earthquake just south of this event.

Southern San Andreas‐San Jacinto fault system slip rates estimated from earthquake cycle models constrained by GPS and interferometric synthetic aperture radar observations

We use ground geodetic and interferometric synthetic aperture radar satellite observations across the southern San Andreas (SAF)-San Jacinto (SJF) fault systems to constrain their slip rates and the viscosity structure of the lower crust and upper mantle on the basis of periodic earthquake cycle, Maxwell viscoelastic, finite element models. Key questions for this system are the SAF and SJF slip rates, the slip partitioning between the two main branches of the SJF, and the dip of the SAF. The best-fitting models generally have a high-viscosity lower crust (η = 10^(21) Pa s) overlying a lower-viscosity upper mantle (η = 10^(19) Pa s). We find considerable trade-offs between the relative time into the current earthquake cycle of the San Jacinto fault and the upper mantle viscosity. With reasonable assumptions for the relative time in the earthquake cycle, the partition of slip is fairly robust at around 24–26 mm/a for the San Jacinto fault system and 16–18 mm/a for the San Andreas fault. Models for two subprofiles across the SAF-SJF systems suggest that slip may transfer from the western (Coyote Creek) branch to the eastern (Clark-Superstition hills) branch of the SJF from NW to SE. Across the entire system our best-fitting model gives slip rates of 2 ± 3, 12 ± 9, 12 ± 9, and 17 ± 3 mm/a for the Elsinore, Coyote Creek, Clark, and San Andreas faults, respectively, where the large uncertainties in the slip rates for the SJF branches reflect the large uncertainty in the slip rate partitioning within the SJF system.

A multiscale approach to estimating topographically correlated propagation delays in radar interferograms

When targeting small amplitude surface deformation, using repeat orbit Interferometric Synthetic Aperture Radar (InSAR) observations can be plagued by propagation delays, some of which correlate with topographic variations. These topographically-correlated delays result from temporal variations in vertical stratification of the troposphere. An approximate model assuming a linear relationship between topography and interferometric phase has been used to correct observations with success in a few studies. Here, we present a robust approach to estimating the transfer function, K, between topography and phase that is relatively insensitive to confounding processes (earthquake deformation, phase ramps from orbital errors, tidal loading, etc.). Our approach takes advantage of a multiscale perspective by using a band-pass decomposition of both topography and observed phase. This decomposition into several spatial scales allows us to determine the bands wherein correlation between topography and phase is significant and stable. When possible, our approach also takes advantage of any inherent redundancy provided by multiple interferograms constructed with common scenes. We define a unique set of component time intervals for a given suite of interferometric pairs. We estimate an internally consistent transfer function for each component time interval, which can then be recombined to correct any arbitrary interferometric pair. We demonstrate our approach on a synthetic example and on data from two locations: Long Valley Caldera, California, which experienced prolonged periods of surface deformation from pressurization of a deep magma chamber, and one coseismic interferogram from the 2007 Mw 7.8 Tocapilla earthquake in northern Chile. In both examples, the corrected interferograms show improvements in regions of high relief, independent of whether or not we pre-correct the data for a source model. We believe that most of the remaining signals are predominately due to heterogeneous water vapor distribution that requires more sophisticated correction methods than those described here.

Multiscale InSAR Time Series (MInTS) analysis of surface deformation

[1] We present a new approach to extracting spatially and temporally continuous ground deformation fields from interferometric synthetic aperture radar (InSAR) data. We focus on unwrapped interferograms from a single viewing geometry, estimating ground deformation along the line-of-sight. Our approach is based on a wavelet decomposition in space and a general parametrization in time. We refer to this approach as MInTS (Multiscale InSAR Time Series). The wavelet decomposition efficiently deals with commonly seen spatial covariances in repeat-pass InSAR measurements, since the coefficients of the wavelets are essentially spatially uncorrelated. Our time-dependent parametrization is capable of capturing both recognized and unrecognized processes, and is not arbitrarily tied to the times of the SAR acquisitions. We estimate deformation in the wavelet-domain, using a cross-validated, regularized least squares inversion. We include a model-resolution-based regularization, in order to more heavily damp the model during periods of sparse SAR acquisitions, compared to during times of dense acquisitions. To illustrate the application of MInTS, we consider a catalog of 92 ERS and Envisat interferograms, spanning 16 years, in the Long Valley caldera, CA, region. MInTS analysis captures the ground deformation with high spatial density over the Long Valley region.

Coseismic fault slip of the 2008 Mw 7.9 Wenchuan earthquake estimated from InSAR and GPS measurements

[1] We infer co-seismic fault slip during the 2008 Mw 7.9 Wenchuan earthquake from interferometric synthetic aperture radar (InSAR) and GPS observations of ground deformation. We use ALOS/PALSAR data from ascending orbits on six tracks, and we do not use data that are strongly affected by ionospheric perturbations. We use a fault model composed of three planar fault segments of the Beichuan fault, and one planar segment representing the parallel Pengguan fault. Maximum thrust-slip is up to 6.7 m near the surface, and occurs in two locations, near Yingxiu in the south and Beichuan in the center of the rupture. Maximum strike-slip is over 4 m, and occurs near Pingtong and Nanba along the northern end of the rupture. We find that the ratio of coseismic thrust- to strike-slip on the Beichuan fault decreases from 1.5 to 0.7 from the SW to the NE. Citation: Feng, G., E. A. Hetland, X. Ding, Z. Li, and L. Zhang (2010), Coseismic fault slip of the 2008 Mw 7.9 Wenchuan earthquake estimated from InSAR and GPS measurements, Geophys. Res. Lett., 37, L01302, doi:10.1029/2009GL041213.

Inference of Multiple Earthquake-Cycle Relaxation Timescales from Irregular Geodetic Sampling of Interseismic Deformation

Abstract Characterizing surface deformation throughout a full earthquake cycle is a challenge due to the lack of high‐resolution geodetic observations of duration comparable to that of characteristic earthquake recurrence intervals (250–10,000 years). Here we approach this problem by comparing long‐term geologic slip rates with geodetically derived fault slip rates by sampling only a short fraction (0.001%–0.1%) of a complete earthquake cycle along 15 continental strike‐slip faults. Geodetic observations provide snapshots of surface deformation from different times through the earthquake cycle. The timing of the last earthquake on many of these faults is poorly known, and may vary greatly from fault to fault. Assuming that the underlying mechanics of the seismic cycle are similar for all faults, geodetic observations from different faults may be interpreted as samples over a significantly larger fraction of the earthquake cycle than could be obtained from the geodetic record along any one fault alone. As an ensemble, we find that geologically and geodetically inferred slip rates agree well with a linear relation of 0.94±0.09. To simultaneously explain both the ensemble agreement between geologic and geodetic slip‐rate estimates with observations of rapid postseismic deformation, we consider the predictions from simple two‐layer earthquake‐cycle models with both Maxwell and Burgers viscoelastic rheologies. We find that a two‐layer Burgers model, with two relaxation timescales, is consistent with observations of deformation throughout the earthquake cycle, whereas the widely used two‐layer Maxwell model with a single relaxation timescale, is not, suggesting that the earthquake cycle is effectively characterized by a largely stress‐recoverable rapid postseismic stage and a much more slowly varying interseismic stage.

New Radar Interferometric Time Series Analysis Toolbox Released

Interferometric synthetic aperture radar (InSAR) has become an important geodetic tool for measuring deformation of Earth’s surface due to various geophysical phenomena, including slip on earthquake faults, subsurface migration of magma, slow‐moving landslides, movement of shallow crustal fluids (e.g., water and oil), and glacier flow. Airborne and spaceborne synthetic aperture radar (SAR) instruments transmit microwaves toward Earth’s surface and detect the returning reflected waves. The phase of the returned wave depends on the distance between the satellite and the surface, but it is also altered by atmospheric and other effects. InSAR provides measurements of surface deformation by combining amplitude and phase information from two SAR images of the same location taken at different times to create an interferogram. Several existing open‐source analysis tools [Rosen et al., 2004; Rosen et al., 2011; Kampes et al., 2003 ; Sandwell et al., 2011] enable scientists to exploit observations from radar satellites acquired at two different epochs to produce a surface displacement map.

An open source Bayesian Monte Carlo isotope mixing model with applications in Earth surface processes

The implementation of isotopic tracers as constraints on source contributions has become increasingly relevant to understanding Earth surface processes. Interpretation of these isotopic tracers has become more accessible with the development of Bayesian Monte Carlo (BMC) mixing models, which allow uncertainty in mixing end-members and provide methodology for systems with multicomponent mixing. This study presents an open source multiple isotope BMC mixing model that is applicable to Earth surface environments with sources exhibiting distinct end-member isotopic signatures. Our model is first applied to new δ18O and δD measurements from the Athabasca Glacier, which showed expected seasonal melt evolution trends and vigorously assessed the statistical relevance of the resulting fraction estimations. To highlight the broad applicability of our model to a variety of Earth surface environments and relevant isotopic systems, we expand our model to two additional case studies: deriving melt sources from δ18O, δD, and 222Rn measurements of Greenland Ice Sheet bulk water samples and assessing nutrient sources from ɛNd and 87Sr/86Sr measurements of Hawaiian soil cores. The model produces results for the Greenland Ice Sheet and Hawaiian soil data sets that are consistent with the originally published fractional contribution estimates. The advantage of this method is that it quantifies the error induced by variability in the end-member compositions, unrealized by the models previously applied to the above case studies. Results from all three case studies demonstrate the broad applicability of this statistical BMC isotopic mixing model for estimating source contribution fractions in a variety of Earth surface systems.

The weight of the mountains: Constraints on tectonic stress, friction, and fluid pressure in the 2008 Wenchuan earthquake from estimates of topographic loading

Though it is widely recognized that large mountain ranges produce significant stresses in the Earths crust, these stresses are not commonly quantified. Nonetheless, near large mountains topography may affect fault activity by changing the stress balance on the faults. In this work, we calculate the stress field from topography in the Longmen Shan (Sichuan, China) and resolve those stresses on several models of the faults that ruptured in the 2008 Mw 7.9 Wenchuan earthquake. We find that the topography results in shear stresses up to 20 MPa and normal stresses up to 80 MPa on the faults, with significant variability across the faults. Topographic stresses generally load the fault in a normal and left-lateral shear sense, opposite to the inferred coseismic slip sense, and thus inhibit the coseismic slip. We estimate the tectonic stress needed to overcome topographic and lithostatic stresses by assuming that the direction of maximum shear accumulated on the faults is roughly collinear with the inferred coseismic slip. We further estimate the static friction and pore fluid pressure assuming that the fault was, on average, at Mohr-Coulomb failure at the time of the Wenchuan earthquake. We use a Bayesian inversion strategy, yielding posterior probability distributions for the estimated parameters. We find most likely estimates of maximum tectonic compressive stress near 0.6 ρgz and oriented ∼E-W, and minimum tectonic stress near 0.2 ρgz. Static friction on the fault is near 0.2, and pore fluid pressure is between 0 and 0.4 of the total pressure.