International audienceGeophysical 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

[1] We present an analytic solution for the deformation near an infinite strike-slip fault in an elastic layer overlying a linear viscoelastic half-space. The theory is valid for any linear viscoelastic rheology and any earthquake sequence. This is a generalization of the work of J. C. Savage and colleagues, which only holds for models with a uniform shear modulus, Maxwell viscoelastic lower half-space, and periodic rupture recurrence. We demonstrate the theory for models of an elastic layer over a Maxwell, standard linear solid, Burgers, and triviscous half-space. For each of these models, we calculate example postseismic displacements, interseismic displacements in a periodic earthquake sequence, and interseismic displacements for a nonperiodic earthquake sequence. Our solution is for a simple geometry; however, the model presents an elegant tool to explore the evolution of displacements for relatively complex rheologies and rupture recurrence histories.

.[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.

[1] 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.

18 Pa s can explain the apparent compression regardless of the assumed mantle viscosity. The geodetic velocities that we determined after removing the postseismic contribution indicate that within current GPS detection limits the BR is stable east of the CNSB and is undergoing rapid right-lateral shear west of the CNSB.