Coseismic slip model of the 2007 August Pisco earthquake (Peru) as constrained by Wide Swath radar observations

The Pisco earthquake ( M _w 8.0; 2007 August 15) occurred offshore of Peru's southern coast at the subduction interface between the Nazca and South American plates. It ruptured a previously identified seismic gap along the Peruvian margin. We use Wide Swath InSAR observations acquired by the Envisat satellite in descending and ascending orbits to constrain coseismic slip distribution of this subduction earthquake. The data show movement of the coastal regions by as much as 85 cm in the line-of-sight of the satellite. Distributed-slip model indicates that the coseismic slip reaches values of about 5.5 m at a depth of ∼18–20 km. The slip is confined to less than 40 km depth, with most of the moment release located on the shallow parts of the interface above 30 km depth. The region with maximum coseismic slip in the InSAR model is located offshore, close to the seismic moment centroid location. The geodetic estimate of seismic moment is 1.23 × 10²¹ Nm ( M _w 8.06), consistent with seismic estimates. The slip model inferred from the InSAR observations suggests that the Pisco earthquake ruptured only a portion of the seismic gap zone in Peru between 13.5° S and 14.5° S, hence there is still a significant seismic gap to the south of the 2007 event that has not experienced a large earthquake since at least 1687.     

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.     

基于改进证据权重法的北海道地震同震滑坡易发性评价

多模型耦合方法被广泛应用于同震滑坡易发性评价研究,但耦合证据权重法和其他方法的改进型证据权重法鲜有涉及。该文以北海道地震震中区为研究区,基于震前和震后高精度遥感影像解译出5977处同震滑坡(以中小型为主),选取高程、坡度、距断层距离、距水系距离等8个滑坡易发性评价因子,对同震滑坡的空间分布规律和控制因子进行分析,发现同震滑坡集中发育在距断层小于5 km、距水系小于800 m、峰值地面加速度(PGA)0.5 g、高程90~250 m、坡度20°~40°、坡向东至东南向、岩性以页岩和砾岩为主的区域内;基于改进型证据权重法,综合考虑地形、地质和水文条件对北海道同震滑坡易发性进行评价,经ROC曲线评价模型验证,发现改进模型曲线下面积(AUC)为0.916,较原模型(AUC为0.870)精度有明显提高。改进模型能为同震滑坡易发性评价提供新的研究思路,并为震后区域防灾减灾提供科学依据。     

Dynamic rupture and stress change in a normal faulting earthquake in the subducting Cocos plate

A large nearly vertical, normal faulting earthquake ( M _w = 7.1) took place in 1997 in the Cocos plate, just beneath the ruptured fault zone of the great 1985 Michoacan thrust event ( M _w = 8.1). Dynamic rupture and resultant stress change during the 1997 earthquake have been investigated on the basis of near-source strong-motion records together with a 3-D dynamic model.
Dynamically consistent waveform inversion reveals a highly heterogeneous distribution of stress drop, including patch-like asperities and negative stress-drop zones. Zones of high stress drop are mainly confined to the deeper, southeastern section of the vertical fault, where the maximum dynamic stress drop reaches 280 bars (28 MPa). The dynamically generated source time function varies with location on the fault, and yields a short slip duration, which is caused by a short scalelength of stress-drop heterogeneities. The synthetic seismograms calculated from the dynamic model are generally consistent with the strong-motion velocity records in the frequency range lower than 0.5 Hz.
The pattern of stress-drop distribution appears, in some sense, to be consistent with that of coseismic changes in shear stress resulting from the 1985 thrust event. This consistency suggests that the stress transfer from the 1985 event to the subducting plate could be one of the possible mechanisms that increased the chance of the occurrence of the 1997 earthquake.     

Temporal and spatial variations of post-seismic deformation following the 1999 Chi-Chi, Taiwan earthquake

We use GPS displacements collected in the 15 months after the 1999 Chi-Chi, Taiwan earthquake  ( M _w 7.6)  to evaluate whether post-seismic deformation is better explained by afterslip or viscoelastic relaxation of the lower crust and upper mantle. We find that all viscoelastic models tested fail to fit the general features in the post-seismic GPS displacements, in contrast to the satisfactory fit obtained with afterslip models. We conclude that afterslip is the dominant mechanism in the 15-month period, and invert for the space–time distribution of afterslip, using the Extended Network Inversion Filter. Our results show high slip rates surrounding the region of greatest coseismic slip. The slip-rate distribution remains roughly stationary over the 15-month period. In contrast to the limited coseismic slip on the décollement, afterslip is prominent there. Maximum afterslip of 0.57 m occurs downdip and to the east of the hypocentral region. Afterslip at hypocentral depths is limited to the southern part of the main shock rupture, with little or no slip on the northern section where coseismic slip was greatest. Whether this results from along strike variations in frictional properties or dynamic conditions that locally favour stable sliding is not clear. In general, afterslip surrounds the area of greatest coseismic slip, consistent with post-seismic slip driven by the main shock stress change. The total accumulated geodetic afterslip moment is  3.8 × 10¹⁹ N m  , significantly more than the seismic moment released by aftershocks,  6.6 × 10¹⁸ N m  . Afterslip and aftershocks appear to have different temporal evolutions and some spatial correlations, suggesting that aftershock rates may not be completely controlled by the rate of afterslip.     

The fractal geometry of the surface ruptures of the 1999 Chi-Chi earthquake, Taiwan

The 1999 September 20 Chi-Chi earthquake is the largest seismic event which occurred in the island during the twentieth century. Available seismic data relative to this earthquake are of high quality, and surface ruptures identified as features associated to the Chelungpu fault can be clearly observed at the surface and precisely mapped. We calculated the fractal dimension ( D ) and b value of Gutenberg–Richter law for 6-month aftershocks of the Chi-Chi earthquake for the fault area, and find that the surface ruptures exhibit self-similar geometry only within specific ruler intervals. The D values of the surface ruptures reflect the fault slip and geometry at depth. More importantly, the small-size aftershocks seem more likely to occur within high D value and high b value areas, whereas small D value and small b value areas have a high potential for medium- and large-size aftershocks.     

A note on the surface volume change of shallow earthquakes

Summary. This note presents an exact analytical formula for determining the magnitude of coseismic surface volume change (δ V ) of earthquake faults in a half-space. For a Poisson solid, the formula is remarkably simple; δ V = M _zz |8μ, where M _zz is one of the moment tensor elements of the source. Maximum δ V values derive from dip slip on faults plunging 45°. For these events, surface volume changes of 0.0001 and 4.3 km³ can be expected for magnitude 5 and 8 earthquakes respectively. All of the coseismic surface volume change is recovered in the interseismic period through relaxation of the Earth and rebound of the surface. A useful rule of thumb for estimating the magnitude of vertical rebound in 45° dip slip events is δ h _p=Δ s /24, where Δ s is the coseismic slip on the fault.     

A kinematic viscoelastic model of the San Francisco earthquake of 1906

Summary. We construct a model of the San Andreas fault zone based on a rectangular fault in an elastic layer overlying a viscoelastic half-space. We alllow both steady and episodic aseismic slip at depth on the fault as well as a large-scale relative plate driving force. We use the model to explain the aseismic changes in geodetic triangulation angles observed during the 40 years following the 1906 San Francisco earthquake. The most important results are that viscoelastic relaxation can explain the data very well, and that the driving force of relative plate motion can be characterized by a horizontal distance scale perpendicular to the plate boundary of hundreds of kilometres.     

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.     

A Bayesian approach to estimating tectonic stress from seismological data

Earthquakes are conspicuous manifestations of tectonic stress, but the non-linear relationships between the stresses acting on a fault plane, its frictional slip, and the ensuing seismic radiation are such that a single earthquake by itself provides little information about the ambient state of stress. Moreover, observational uncertainties and inherent ambiguities in the nodal planes of earthquake focal mechanisms preclude straightforward inferences about stress being drawn on the basis of individual focal mechanism observations. However, by assuming that each earthquake in a small volume of the crust represents a single, uniform state of stress, the combined constraints imposed on that stress by a suite of focal mechanism observations can be estimated. Here, we outline a probabilistic (Bayesian) technique for estimating tectonic stress directions from primary seismological observations. The Bayesian formulation combines a geologically motivated prior model of the state of stress with an observation model that implements the physical relationship between the stresses acting on a fault and the resultant seismological observation. We show our Bayesian formulation to be equivalent to a well-known analytical solution for a single, errorless focal mechanism observation. The new approach has the distinct advantage, however, of including (1) multiple earthquakes, (2) fault plane ambiguities, (3) observational errors and (4) any prior knowledge of the stress field. Our approach, while computationally demanding in some cases, is intended to yield reliable tectonic stress estimates that can be confidently compared with other tectonic parameters, such as seismic anisotropy and geodetic strain rate observations, and used to investigate spatial and temporal variations in stress associated with major faults and coseismic stress perturbations.     

Late Quaternary activity of faults and recurrence interval of earthquakes in the eastern Hokuriku region, northern central Japan, on the basis of precise cryptotephra analysis of fluvial terrace sequences

About 2000 active faults are known to exist within the land area of Japan. Most of these active faults have deformed the topographic surfaces which were formed in the late Quaternary, including fluvial terraces; and the formative ages of these terraces are estimated mainly by tephrochronology. Fluvial terraces in the eastern Hokuriku region, comprising the Toyama, Tonami, and Kanazawa Plains, northern central Japan, are widely distributed and have been deformed by reverse active faults. The formative age of terraces in this area has not been reported, as volcanic ash deposits are rarely visible within terrace deposits and the overlying loamy soil, and outcrops of fluvial terraces are quite scarce in this area. In the present study, we carried out a drilling survey on these terraces to obtain samples of the overlying loamy soil and upper part of terrace deposits. From these samples, we extracted some well-known widespread volcanic ash, from which we were able to estimate the approximate age of the terraces and the vertical slip rate of the active faults. Late Quaternary fluvial terraces in eastern Hokuriku are divided into 12 levels: Terraces 1 to 12 in descending order. Widespread tephras such as the Kikai-Tozurahara Tephra (K-Tz: 95 ka) are contained in the lowest part of the loamy soil in Terrace 4 and the Daisen-Kurayoshi Pumice (DKP: 55 ka) is present in the lowest part of the loamy soil in Terrace 6. From the ages and the vertical displacements of the fluvial terraces, the late Quaternary average vertical slip rates of active faults in eastern Hokuriku are estimated to be 0.2–0.9 mm/year (Uozu fault), 0.1–0.4 mm/year (Kurehayama fault), 0.1–0.3 mm/year (Takashozu fault), 0.1–0.4 mm/year (Hohrinji fault), and 0.5–0.8 mm/year (Morimoto-Togashi fault). We also estimated the recurrence interval of earthquakes related to active faults from displacement per event and ages of terraces and no significant difference in vertical displacement per single earthquake for different active faults, and recurrence intervals tend to be inversely proportional to vertical displacement rates. This study demonstrates that a combination of drilling of loamy soil and precise cryptotephra analysis of fluvial terraces can be used to estimate the formative age of the terraces and the average slip rate of active faults in areas where volcanic ash deposits are rare.     

Stress transfer relations among the earthquakes that occurred in Kerman province, southern Iran since 1981

We explore the possible stress triggering relationship of the   M ≥ 6.4  earthquakes that occurred in Kerman Province, southern Iran since 1981. We calculated stress changes due to both coseismic sudden movement in the upper crust and the time-dependent viscous relaxation of the lower crust and/or upper mantle following the event. Four events of   M ≥ 6.4  occurred between 1981 and 2005, on and close to the Gowk fault, show a clear Coulomb stress load to failure relationship. The  2003 M = 6.5  Bam earthquake, however, which occurred approximately 95 km SW of the closest Gowk event, shows a very weak stress relation to preceding earthquakes. The coseismic static stress change at the hypocentre of the Bam earthquake is quite small (∼0.006 bars). The time-dependent post-seismic stress change could be 26 times larger or 7 times lower than that of coseismic static stress alone depending on the choice of viscoelastic crustal model and the effective coefficient of friction. Given the uncertainties in the viscoelastic earth models and the effective coefficient of friction, we cannot confidently conclude that the 2003 Bam event was brought closer to failure through coseismic or post-seismic stress loading. Interestingly, the southern Gowk segment with a similar strike to that of the Bam fault, experienced a stress load of up to 8.3 bars between 1981 and 2003, and is yet to have a damaging earthquake.     

昌马断裂带活动构造地貌之研究

本文根据野外地质调查资料和人工探槽的剖折,并结合¹⁴C年龄数据,讨论了昌马活动断裂带晚更新世晚期以来形成的构造地貌的类型、特征及分期问题.     

Geodetic inversion for space–time distribution of fault slip with time-varying smoothing regularization

We have developed a new geodetic inversion method for space–time distribution of fault slip velocity with time-varying smoothing regularization in order to reconstruct accurate time histories of aseismic fault slip transients. We introduce a temporal smoothing regularization on slip and slip velocity through a Bayesian state space approach in which the strength of regularization (temporal smoothness of slip velocity) is controlled by a hyperparameter. The time-varying smoothing regularization is realized by treating the hyperparameter as a time-dependent stochastic variable and adopting a hierarchical Bayesian state space model, in which a prior distribution on the hyperparameter is introduced in addition to a conventional Bayesian state space model. We have tested this inversion method on two synthetic data sets generated by simulated aseismic slip transients. Results show that our method reproduces well both rapid changes of slip velocity and steady-state velocity without significant oversmoothing and undersmoothing, which has been hard to overcome by the conventional Bayesian approach with time-independent smoothing regularization. Application of this method to transient deformation in 2002 caused by a silent earthquake off the Boso peninsula, Japan, also shows similar advantages of this method over the conventional approach.     

High-frequency radiation from crack (stress drop) models of earthquake faulting

We present a theory for the radiation of high-frequency waves by earthquake faults. We model the fault as a planar region in which the stress drops to the kinematic friction during slip. This model is entirely equivalent to a shear crack. For two-dimensional fault models we show that the high frequencies originate from the stress and slip velocity concentrations in the vicinity of the fault's edges. These stress concentrations radiate when the crack expands with accelerated motion. The most efficient generation of high-frequency waves occurs when the rupture velocity changes abruptly. In this case, the displacement spectrum has an ω^-2 behaviour at high frequencies. The excitation is proportional to the intensity of the stress concentration near the crack tips and to the change in the focusing factor due to rupture velocity. We extend these two-dimensional results to more general three-dimensional fault models in the case when the rupture velocity changes simultaneously on the rupture front. Results are similar to those described for two-dimensional faults. We apply the theory to the case of a circular fault that grows at constant velocity and stops suddenly. The present theory is in excellent agreement with a numerical solution of the same problem.
Our results provide upper bounds to the high-frequency radiation from more realistic models in which rupture velocity does not change suddenly. The ω^-2 is the minimum possible decay at high frequencies for any crack model of the source.     

Tectonic geomorphology of the Ash Hill fault, Panamint Valley, California

ABSTRACT
Panamint Valley, in eastern California, is an extensional basin currently bounded by active, dextral-normal oblique-slip faults. There is considerable debate over the tectonic and topographic evolution of the valley. The least-studied structure, the Ash Hill fault, runs for some 50 km along the valley's western edge, and active strands of the fault continue south into the neighbouring Slate Range. Vertical displacement on the fault is valley-side up, creating topography that conflicts with the gross morphology of the valley itself. We use this topography, along with kinematic and geological markers, to constrain the Quaternary slip rate and orientation of the Ash Hill fault. The fault offsets all but the active channel deposits in the valley, and slickenlines indicate a strike-slip to dip-slip ratio of 3.5:1. An offset volcanic unit dated at 4 Ma provides a minimum slip rate of 0.3±0.1 mm yr⁻¹, and a long-term strike-slip to dip-slip ratio of 5.2:1. Slip on the fault has warped a palaeolake shoreline within the valley. Simple elastic dislocation modelling of the vertical deformation of the shoreline suggests total fault slip of ≈60 m, valley-side up. The shoreline probably dates to 120–150 ka, implying a late Quaternary slip rate of 0.4–0.5 mm yr⁻¹. We suggest two possible mechanisms for the apparently anomalous slip behaviour of the Ash Hill fault. The fault may be a listric structure related to the proposed low-angle fault underlying Panamint Valley. Alternatively, the Ash Hill fault is a high-angle fault, implying that the valley is currently bounded by high-angle dextral-slip faults. Lack of detailed subsurface information precludes any knowledge of the true relationships between the presently active faults.     

Simulations of strong ground motion in SW Iberia for the 1969 February 28 (Ms= 8.0) and the 1755 November 1 (M∼ 8.5) earthquakes – II. Strong ground motion simulations

This is the second paper of a series of two concerning strong ground motion in SW Iberia due to earthquakes originating from the adjacent Atlantic area. The aim of this paper is to use the velocity model that was proposed and validated in the companion paper for seismic intensity modelling of the 1969 ( M _s= 8.0) and 1755 ( M = 8.5–8.7) earthquakes.
First, we propose a regression to convert simulated values of Peak Ground Velocity (PGV) into Modified Mercalli Intensity (MMI) in SW Iberia, and using this regression, we build synthetic isoseismal maps for a large ( M _s= 8.0) earthquake that occurred in 1969. Based on information on the seismic source provided by various authors, we show that the velocity model effectively reproduces macroseismic observations in the whole region. We also confirm that seismic intensity distribution is very sensitive to a small number of source parameters: rupture directivity, fault strike and fault dimensions. Then, we extrapolate the method to the case of the great ( M = 8.5–8.7) 1755 earthquake, for a series of hypotheses recently proposed by three authors about the location of the epicentral region. The model involving a subduction-related rupture in the Gulf of Cádiz results in excessive ground motion in northern Morocco, suggesting that the source of the 1755 earthquake should be located further west. A rupture along the western coast of Portugal, compatible with an activation of the passive western Iberian margin, would imply a relatively low average slip, which, alone, would could not account for the large tsunami observed in the whole northern Atlantic ocean. A seismic source located below the Gorringe Bank seems the most likely since it is more efficient in reproducing the distribution of high intensities in SW Iberia due to the 1755 earthquake.     

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.     

Source investigation of a small event using empirical Green's functions and simulated annealing

We propose a two-step inversion of three-component seismograms that (1) recovers the far-field source time function at each station and (2) estimates the distribution of co-seismic slip on the fault plane for small earthquakes (magnitude 3 to 4). The empirical Green's function (EGF) method consists of finding a small earthquake located near the one we wish to study and then performing a deconvolution to remove the path, site, and instrumental effects from the main-event signal.
The deconvolution between the two earthquakes is an unstable procedure: we have therefore developed a simulated annealing technique to recover a stable and positive source time function (STF) in the time domain at each station with an estimation of uncertainties. Given a good azimuthal coverage, we can obtain information on the directivity effect as well as on the rupture process. We propose an inversion method by simulated annealing using the STF to recover the distribution of slip on the fault plane with a constant rupture-velocity model. This method permits estimation of physical quantities on the fault plane, as well as possible identification of the real fault plane.
We apply this two-step procedure for an event of magnitude 3 recorded in the Gulf of Corinth in August 1991. A nearby event of magnitude 2 provides us with empirical Green's functions for each station. We estimate an active fault area of 0.02 to 0.15 km² and deduce a stress-drop value of 1 to 30 bar and an average slip of 0.1 to 1.6 cm. The selected fault of the main event is in good agreement with the existence of a detachment surface inferred from the tectonics of this half-graben.     

Mechanics of normal fault networks subject to slip-weakening friction

We seek to understand how the stress interactions and the slip-weakening process combine within a non-coplanar, normal fault network to allow a slip instability to develop, and shape the final slip distribution on the system. In a first part, we perform a non-linear spectral analysis to investigate the conditions of stability and the process of slip initiation in an antiplane non-coplanar fault system subject to a slip-dependent friction law. That numerical model allows determining the zones that are able to slip within a fault network, as well as the location of the stress singularities. The resulting slip profiles on the faults show only a few different shapes, some of them with long, linear sections. This leads to formulate a general classification of slip profiles that can be used to infer the degree of fault interaction within any non-coplanar system. In a second part of work, we use our modelling to try reproducing the cumulative slip profiles measured on three real normal interacting faults forming a large-scale en echelon system. For that, we assume that cumulative slip profiles can be compared to the first static modal solution of our conceptual model. We succeed reproducing the profiles quite well using a variable weakening along the faults. Overall, the weakening rate decreases in the direction of propagation of the fault system. Yet, modelling the slip along the propagating, isolated termination segment of the system requires an unlikely distribution of weakening. This suggests that factors not considered in our analysis may contribute to slip profile shaping on isolated, propagating faults.