Multi-interferogram method for measuring interseismic deformation: Denali Fault, Alaska

Studies of interseismic strain accumulation are crucial to our understanding of continental deformation, the earthquake cycle and seismic hazard. By mapping small amounts of ground deformation over large spatial areas, InSAR has the potential to produce continental-scale maps of strain accumulation on active faults. However, most InSAR studies to date have focused on areas where the coherence is relatively good (e.g. California, Tibet and Turkey) and most analysis techniques (stacking, small baseline subset algorithm, permanent scatterers, etc.) only include information from pixels which are coherent throughout the time-span of the study. In some areas, such as Alaska, where the deformation rate is small and coherence very variable, it is necessary to include information from pixels which are coherent in some but not all interferograms. We use a three-stage iterative algorithm based on distributed scatterer interferometry. We validate our method using synthetic data created using realistic parameters from a test site on the Denali Fault, Alaska, and present a preliminary result of  10.5 ± 5.0  mm yr⁻¹ for the slip rate on the Denali Fault based on a single track of radar data from ERS1/2.     

Dynamic stress field of a kinematic earthquake source model with k-squared slip distribution

Source models such as the k -squared stochastic source model with k -dependent rise time are able to reproduce source complexity commonly observed in earthquake slip inversions. An analysis of the dynamic stress field associated with the slip history prescribed in these kinematic models can indicate possible inconsistencies with physics of faulting. The static stress drop, the strength excess, the breakdown stress drop and critical slip weakening distance D _c distributions are determined in this study for the kinematic k -squared source model with k -dependent rise time. Several studied k -squared models are found to be consistent with the slip weakening friction law along a substantial part of the fault. A new quantity, the stress delay, is introduced to map areas where the yielding criterion of the slip weakening friction is violated. Hisada's slip velocity function is found to be more consistent with the source dynamics than Boxcar, Brune's and Dirac's slip velocity functions. Constant rupture velocities close to the Rayleigh velocity are inconsistent with the k -squared model, because they break the yielding criterion of the slip weakening friction law. The bimodal character of D _c / D _tot frequency–magnitude distribution was found. D _c approaches the final slip D _tot near the edge of both the fault and asperity. We emphasize that both filtering and smoothing routinely applied in slip inversions may have a strong effect on the space–time pattern of the inferred stress field, leading potentially to an oversimplified view of earthquake source dynamics.     

Coseismic deformation revealed by inversion of strong motion and GPS data: the 2003 Chengkung earthquake in eastern Taiwan

A moderate earthquake of   M _w= 6.8  occurred on 2003 December 10. It ruptured the Chihshang Fault in eastern Taiwan which is the most active segment of the Longitudinal fault as a plate suture fault between the Luzon arc of the Philippine Sea plate and the Eurasian plate. The largest coseismic displacements were 13 cm (horizontal) and 26 cm (vertical). We analyse 40 strong motion and 91 GPS data to model the fault geometry and coseismic dislocations. The most realistic shape of the Chihshang fault surface is listric in type. The dipping angle of the seismic zone is steep (about 60°–70°) at depths shallower than 10 km and then gradually decreases to 40°–50° at depths of 20–30 km. Thus the polygonal elements in Poly3D are well suited for modelling complex surfaces with curving boundaries. Using the strong motion data, the displacement reaches 1.2 m dip-slip on the Chihshang Fault and decreases to 0.1 m near surface. The slip averages 0.34 m, releasing a scalar moment of 1.6E26 dyne-cm. For GPS data, our model reveals that the maximal dislocation is 1.8 m dip-slip. The dislocations decrease to 0.1 m near the surface. The average slip is 0.48 m, giving a scalar moment of 2.2E26 dyne-cm. Regarding post-seismic deformation, a displacements of 0.5 m were observed near the Chihshang Fault, indicating the strain had not been totally released, as a probable result of near-surface locking of the fault zone.