Palaeoseismic records in lacustrine sediments—A case study of the Daqingshan piedmont fault and Hasuhai Lake in Inner Mongolia,China

As a traditional method for palaeoseismic studies, trenching can be combined with dating techniques to identify palaeoseismic events and the earthquake recurrence interval. However, when using trenches to study palaeoearthquakes, factors such as the active tectonic background of the earthquake‐caused structure, the lithology on both sides of the fault, the geomorphology location and type and the samples and methods for dating will affect the location of the trench. Thus, trenches should be carefully selected and used to identify the impact of ancient earthquakes. The results have substantial uncertainties and limitations. In recent years, scholars have made considerable progress in using other methods to reveal the palaeoseismic information of faults. Moreover, the history of fault activity may have been recorded in the lacustrine sediment adjacent to the fault. Hasuhai Lake is adjacent to the middle segment of the Daqingshan piedmont fault in Inner Mongolia. Since the Holocene, the region has experienced a temperate continental semi‐arid climate with little interference, and Hasuhai Lake and peripheral waters present weak hydrodynamic conditions that provide an ideal location for the study of palaeoseismic records in lacustrine sediments. Sediment samples and samples for dating were collected from three trenches excavated on the periphery of the Hasuhai Lake. Their variations in grain size and magnetic susceptibility revealed that wind and flowing water jointly produced the sedimentary conditions of Hasuhai sediments. The ¹⁴C dating results and variations in the grain size distribution, grain size components and magnetic susceptibility of sediments caused by seismic events obtained in this study were compared with those caused by a series of palaeoseismic events at the middle segment of the Daqingshan piedmont fault reported by previous studies using trenches, knickpoints and palaeosol records. The results identified seven palaeoseismic events recorded near Hasuhai Lake since 12,000 years. The combined use of lacustrine sediment variation characteristics and dating techniques is an effective method for studying palaeoseismic events.     

Amplification of ground motion and waveform complexity in fault zones: examples from the San Andreas and Calaveras Faults

Summary. P -wave seismograms at ranges less than 10 km are synthesized by asymptotic ray theory and by summation of Gaussian beams for point sources located in a low-velocity wedge surrounding a fault. The computations are performed using models of the wedge inferred from the analysis of reflection and refraction experiments across the San Andreas and Hayward-Calaveras faults. Calculations in these models show that the 10–20Hz vertical displacements of earthquakes located at 3–10km depth are amplified by up to an order of magnitude in a 1–2km wide region centred on the fault trace compared to displacements predicted by laterally homogeneous models of the crust. This amplification is not cancelled by high attentuation in the fault zone and compensates for the reduction in amplitudes directly above the source predicted from the radiation pattern of a strike-slip earthquake. Depending on the source depth of the earthquake and the structure and velocity contrast of the wedge, multiple triplications in the travel-time curve of direct P - and S -waves will occur at stations in the fault zone. A wedge model successfully predicts the triplications observed in the P waveforms of aftershocks of the Coyote Lake earthquake recorded in the fault zone, showing that body waves from microearthquakes can be used to determine the three-dimensional velocity structure of the fault zone. The amplification, waveform complexity, and distortion of ray paths introduced by the low- velocity wedge suggest that its effects should be included in the interpretation of strong ground motions and travel times observed in the fault zone. For realistic models of the wedge, asymptotically approximate methods of calculating the body waveforms are strictly valid for frequencies greater than 20Hz. Numerical methods may be necessary to calculate accurately the wavefield at lower frequencies.     

A multidisciplinary study of a slow-slipping fault for seismic hazard assessment: the example of the Middle Durance Fault (SE France)

Assessing seismic hazard in continental interiors is difficult because these regions are characterized by low strain rates and may be struck by infrequent destructive earthquakes. In this paper, we provide an example showing that interpretations of seismic cross sections combined with other kinds of studies such as analysis of microseismicity allow the whole seismogenic source area to be imaged in this type of region. The Middle Durance Fault (MDF) is an 80-km-long fault system located southeastern France that has a moderate but regular seismicity and some palaeoseismic evidence for larger events. It behaves as an oblique ramp with a left-lateral-reverse fault slip and has a low strain rate. MDF is one of the rare slow active fault system monitored by a dedicated dense velocimetric short period network. This study showed a fault system segmented in map and cross section views which consists of staircase basement faults topped by listric faults ramping off Triassic evaporitic beds. Seismic sections allowed the construction of a 3-D structural model used for accurate location of microseismicity. Southern part of MDF is mainly active in the sedimentary cover. In its northern part and in Alpine foreland, seismicity deeper than 8 km was also recorded meaning active faults within the crust cannot be excluded. Seismogenic potential of MDF was roughly assessed. Resulting source sizes and estimated slip rates imply that the magnitude upper limit ranges from 6.0 to 6.5 with a return period of a few thousand years. The present study shows that the coupling between 3-D fault geometry imaging and accurate location of microseismicity provides a robust approach to analyse active fault sources and consequently a more refined seismic hazard assessment.     

Estimating rupture scenario likelihood based on dynamic rupture simulations: the example of the segmented Middle Durance fault, southeastern France

The Middle Durance fault system, southeastern France, is a slow active fault that produced moderate-size historical seismic events and shows evidence of at least one   M _w ≳ 6.5  event in the last 29 000 yr. Based on dynamic rupture simulation, we propose earthquake scenarios that are constrained by knowledge of both the tectonic stress field and of the 3-D geometry of the Durance fault system. We simulate dynamic rupture interaction among several fault segmentations of different strikes, dips and rakes, using a 3-D boundary integral equation method. 50 combinations of reasonable stress field orientations, stress field amplitudes and hypocentre locations are tested. The probability of different rupture evolutions is then computed. Each segment ruptures mainly as a single event (44 per cent of the 50 simulations test in this paper). However, the probability that an event triggers simultaneously along three segments is high (26 per cent), leading to a potential rupture length of 45 km. Finally, 2 per cent of the simulations occur along four adjacent segments, producing the greatest total rupture length of 55 km. The simulation results show that the southernmost segment is most easily ruptured (40 per cent), because of its favourable orientation with respect to the tectonic stress and of its favourable location for interaction with the other segments. South-bound unilateral propagation is slightly preferable (41 per cent), compared to north-bound unilateral and bilateral propagation modes. Although, these rupture scenarios cannot be directly translated into probabilities of occurrence, they do provide a better insight as to which rupture scenarios are more likely, an important element to better estimate near-field strong ground motion and seismic hazard.     

Minimum work, fault activity and the growth of critical wedges in fold and thrust belts

Many studies of critical wedges treat the interior of the wedge as continuous and do not address the manner in which it grows from the undeformed state to a typical imbricate wedge. In this paper we present a 2D kinematic–mechanical model which attempts to explain the development of a critical wedge in a fold and thrust belt in terms of both gravitational and frictional work. In the undeformed model a series of thrust faults are defined which have the potential to take up an external displacement. The active fault at a given time is that which minimizes gravitational and frictional work as a result of displacement. Displacement on the active fault causes a change in topography and deformation of other faults which may favour an alternative fault at the next time step. The model is a mixed Lagrangian–Eulerian scheme in which the upper surface, in addition to being deformed, is also subject to erosion, transport and sedimentation. The model predicts propagation of thrust fault activity towards the foreland through time as a result of increasing topographic (gravitational) loads and frictional work on deformed hinterland faults. As the zone of fault activity progresses through the developing critical wedge several faults are active over time-scales of ≈1 Myr. However, a simple chronology or sequence of fault activity cannot be assumed as out-of-sequence thrusting occurs during this overall foreland propagation. The detailed spatial and temporal activity of faults is complex and reflects the interaction between the development of topography, the contrast between basal (décollement) and internal coefficients of friction and the effects of erosion and sedimentation. In particular, rates of erosion and sedimentation are found to be important controls on fault activity both spatially and temporally. Erosion, by locally removing topography above a fault, reduces gravitational and frictional work enabling continued fault activity or reactivation. Sedimentation, conversely, acts to increase gravitational and frictional work on a fault, and therefore has the potential to blanket faults and render them inactive. Model results illustrate the complex feedbacks that can exist between tectonic and surficial mass transport processes.     

Quaternary normal faulting in southeastern Sicily (Italy):a seismic source for the 1693 large earthquake

We present geological and morphological data, combined with an analysis of seismic reflection lines across the Ionian offshore zone and information on historical earthquakes, in order to yield new constraints on active faulting in southeastern Sicily. This region, one of the most seismically active of the Mediterranean, is affected by WNW–ESE regional extension producing normal faulting of the southern edge of the Siculo–Calabrian rift zone. Our data describe two systems of Quaternary normal faults, characterized by different ages and related to distinct tectonic processes. The older NW–SE-trending normal fault segments developed up to ≈400  kyr ago and, striking perpendicular to the main front of the Maghrebian thrust belt, bound the small basins occurring along the eastern coast of the Hyblean Plateau. The younger fault system is represented by prominent NNW–SSE-trending normal fault segments and extends along the Ionian offshore zone following the NE–SW-trending Avola and Rosolini–Ispica normal faults. These faults are characterized by vertical slip rates of 0.7–3.3  mm  yr ⁻¹ and might be associated with the large seismic events of January 1693. We suggest that the main shock of the January 1693 earthquakes ( M ~ 7) could be related to a 45  km long normal fault with a right-lateral component of motion. A long-term net slip rate of about 3.7  mm  yr ⁻¹ is calculated, and a recurrence interval of about 550 ± 50  yr is proposed for large events similar to that of January 1693.     

Fault architecture, basin structure and evolution of the Gulf of Corinth Rift, central Greece

The style of extension and strain distribution during the early stages of intra-continental rifting is important for understanding rift-margin development and can provide constraints for lithospheric deformation mechanisms. The Corinth rift in central Greece is one of the few rifts to have experienced a short extensional history without subsequent overprinting. We synthesise existing seismic reflection data throughout the active offshore Gulf of Corinth Basin to investigate fault activity history and the spatio-temporal evolution of the basin, producing for the first time basement depth and syn-rift sediment isopachs throughout the offshore rift. A major basin-wide unconformity surface with an age estimated from sea-level cycles at ca . 0.4 Ma separates distinct seismic stratigraphic units. Assuming that sedimentation rates are on average consistent, the present rift formed at 1–2 Ma, with no clear evidence for along-strike propagation of the rift axis. The rift has undergone major changes in relative fault activity and basin geometry during its short history. The basement depth is greatest in the central rift (maximum ∼3 km) and decreases to the east and west. In detail however, two separated depocentres 20–50 km long were created controlled by N- and S-dipping faults before 0.4 Ma, while since ca . 0.4 Ma a single depocentre (80 km long) has been controlled by several connected N-dipping faults, with maximum subsidence focused between the two older depocentres. Thus isolated but nearby faults can persist for timescales ca . 1 Ma and form major basins before becoming linked. There is a general evolution towards a dominance of N-dipping faults; however, in the western Gulf strain is distributed across several active N- and S-dipping faults throughout rift history, producing a more complex basin geometry.     

Post-seismic reloading and temporal clustering on a single fault

Geological studies show evidence for temporal clustering of large earthquakes on individual fault systems. Since post-seismic deformation due to the inelastic rheology of the lithosphere may result in a variable loading rate on a fault throughout the interseismic period, it is reasonable to expect that the rheology of the non-seismogenic lower crust and mantle lithosphere may play a role in controlling earthquake recurrence times. We study this phenomenon using a 2-D, finite element method continuum model of the lithosphere containing a single strike-slip fault. This model builds on a previous study using a 1-D spring-dashpot-slider analogue of a single fault system to study the role of Maxwell viscoelastic relaxation in producing non-periodic earthquakes. In our 2-D model, the seismogenic portion of the fault slips when a predetermined yield stress is exceeded; stress accumulated on the seismogenic fault is shed to the viscoelastic layers below and recycled back to the seismogenic fault through viscoelastic relaxation. We find that random variation of the fault yield stress from one earthquake to the next can cause the earthquake sequence to be clustered; the amount of clustering depends on a non-dimensional number, W , called the Wallace number defined as the standard deviation of the randomly varied fault yield stress divided by the effective viscosity of the system times the tectonic loading rate. A new clustering metric based on the bimodal distribution of interseismic intervals allows us to investigate clustering behaviour of systems over a wide range of model parameters and those with multiple viscoelastic layers. For models with   W ≥ 1  clustering increases with increasing W , while those with   W ≤ 1  are unclustered, or quasi-periodic.     

Finite-frequency traveltime tomography of San Francisco Bay region crustal velocity structure

Seismic velocity structure of the San Francisco Bay region crust is derived using measurements of finite-frequency traveltimes. A total of 57 801 relative traveltimes are measured by cross-correlation over the frequency range 0.5–1.5 Hz. From these are derived 4862 'summary' traveltimes, which are used to derive 3-D P -wave velocity structure over a 341 × 140 km² area from the surface to 25 km depth. The seismic tomography is based on sensitivity kernels calculated on a spherically symmetric reference model. Robust elements of the derived P -wave velocity structure are: a pronounced velocity contrast across the San Andreas fault in the south Bay region (west side faster); a moderate velocity contrast across the Hayward fault (west side faster); moderately low velocity crust around the Quien Sabe volcanic field and the Sacramento River delta; very low velocity crust around Lake Berryessa. These features are generally explicable with surface rock types being extrapolated to depth ∼10 km in the upper crust. Generally high mid-lower crust velocity and high inferred Poisson's ratio suggest a mafic lower crust.     

Analogue modelling of inverted domino‐style basement fault systems

Inversion of pre‐existing extensional fault systems is common in rift systems, back‐arc basins and passive margins. It can significantly influence the development of structural traps in hydrocarbon basins. The analogue models of domino‐style basement fault systems shown in this paper produced, on extension, characteristic hangingwall growth stratal wedges that, when contracted and inverted, formed classic inversion harpoon geometries and asymmetric hangingwall contractional fault‐propagation folds. Segmented footwall shortcut faults formed as the basement faults were progressively back‐rotated and steepened. The pre‐existing extensional fault architectures, basement fault geometries and the relative hangingwall and footwall block rotations exerted fundamental controls on the inversion styles. Digital image correlation (DIC) strain monitoring illustrated complex vertical fault segmentation and linkage during inversion as the major faults were reactivated and strain was progressively transferred onto footwall shortcut faults. Hangingwall deformation during inversion was dominated by significant back‐rotation as the inversion progressed. The mechanical stratigraphy of the cover sequences strongly influenced the fold and fault evolution of the reactivated fault systems. The implications of the experimental results for the interpretation and analysis of inversion structures are discussed and are compared with natural examples of inverted basement‐involved extensional faults observed in seismic datasets.     

Detection of fluid migration pathways in seismic data: implications for fault seal analysis

A new and efficient method for fault seal analysis using seismic data is presented. It uses multiple seismic attributes and neural networks to enhance fluid migration pathways, including subtle features that are not detectable using single attributes only. The method may be used as a first estimate of fault seal or to calibrate results from other techniques. The results provide information about which faults and fault segments are sealing or leaking. Fluid flow along individual faults appears to be focused along zones of weakness, and fault seal research should thus be focused on finding such weak locations within fault zones, a task that is best done using three‐dimensional (3D) seismic data. Under certain conditions, it is suggested that fluids migrate along fault planes by a diapiric fluid flow mechanism. The results assist in calibrating the bulk hydraulic properties of faults and rock formations and can be used in basin modelling.     

Evidence for different scaling of earthquake source parameters for large earthquakes depending on faulting mechanism

Scaling relationships between seismic moment, rupture length, and rupture width have been examined. For this purpose, the data from several previous studies have been merged into a database containing more than 550 events. For large earthquakes, a dependence of scaling on faulting mechanism has been found. Whereas small and large dip-slip earthquakes scale in the same way, the self-similarity of earthquakes breaks down for large strike-slip events. Furthermore, no significant differences in scaling could be found between normal and reverse earthquakes and between earthquakes from different regions. Since the thickness of the seismogenic layer limits fault widths, most strike-slip earthquakes are limited to rupture widths of between 15 and 30 km while the rupture length is not limited. The aspect ratio of dip-slip earthquakes is similar for all earthquake sizes. Hence, the limitation in rupture width seems to control the maximum possible rupture length for these events. The different behaviour of strike-slip and dip-slip earthquakes can be explained by rupture dynamics and geological fault growth. If faults are segmented, with the thickness of the seismogenic layer controlling the length of each segment, strike-slip earthquakes might rupture connected segments more easily than dip-slip events, and thus could produce longer ruptures than dip-slip events of the same width     

Morphostructure,tectono‐sedimentary evolution and seismic potential of the Horseshoe Fault,SW Iberian Margin

High‐resolution acoustic and seismic data acquired 100 km offshore Cape São Vicente, image with unprecedented detail one of the largest active reverse faults of the SW Iberian Margin, the Horseshoe Fault (HF). The HF region is an area seismogenically active, source of the largest magnitude instrumental and historical earthquake (M_w > 6) occurred in the SW Iberian Margin. The HF corresponds to a N40 trending, 110 km long, and NW‐verging active thrust that affects the whole sedimentary sequence and reaches up to the seafloor, generating a relief of more than 1 km. The along‐strike structural variability as well as fault trend suggests that the HF is composed by three main sub‐segments: North (N25), Central (N50) and South (N45). Swath‐bathymetry, TOBI sidescan sonar backscatter and parametric echosounder TOPAS profiles reveal the surface morphology of the HF block, characterized by several, steep (20°) small scarps located on the hangingwall, and a succession of mass transport deposits (i.e. turbidites) on its footwall, located in the Horseshoe Abyssal Plain. A succession of pre‐stack depth‐migrated multichannel seismic reflection profiles across the HF and neighbouring areas allowed us to constrain their seismo‐stratigraphy, structural geometry, tectono‐sedimentary evolution from Upper Jurassic to present‐day, and to calculate their fault parameters. Finally, on the basis of segment length, surface fault area and seismogenic depth we evaluated the seismic potential of the HF, which in the worst‐case scenario may generate an earthquake of magnitude M_w 7.8 ± 0.1. Thus, considering the tectonic behaviour and near‐shore location, the HF should be recognized in seismic and tsunami hazard assessment models of Western Europe and North Africa.     

Reactivation of basement faults: interplay of ice‐sheet advance,glacial lake formation and sediment loading

The Emme Delta is a small glacilacustrine delta, which developed on the southern flank of the Wesergebirge Mountains in NW Germany. Shallow shear‐wave seismic surveys allow a detailed assessment of the structural style of the delta body. Two different fault systems are developed within the delta, both showing syn‐sedimentary activity. The faults have planar to slightly listric geometries and show vertical offsets in a range of 2–15 m. They form small graben and half‐graben systems, which locally show roll‐over structures. The fill of the half‐grabens has a wedge‐shaped geometry, with the greatest sediment thickness close to the fault. The fault system in the upper portion of the Emme Delta is restricted to the delta body and probably gravity induced. In the lower portion of the delta, normal faults occur that originate in the underlying Jurassic basement rocks and penetrate into the delta deposits. The grid of seismic lines shows that the normal faults are trending E–W. This fits to a late Triassic–early Jurassic deformation phase in the Central European Basin System. We hypothese that these faults were reactivated during the Pleistocene by the advancing ice‐sheet, water and sediment loading. Based on the seismic data set, an overall model for the reactivation of the basement fault was developed. The advancing ice‐sheet caused far field extension, which might have reactivated pre‐existing normal faults. Later, the fault activity was enhanced due to sediment and water loading. In addition, high pore pressure due to lake formation might have supported the slip processes along the faults. After glacial unloading and lake drainage, the fault activity stopped.     

Structural style and stratigraphic architecture of fault propagation folding in extensional settings: a seismic example from the Smørbukk area, Halten Terrace, Mid-Norway

A number of recent papers have stressed the importance of lateral and vertical fault propagation on sediment geometries in active rift settings. However, the majority of these studies have been based on outcrop data. This contribution addresses the evolution of a single, major normal fault and its interaction with adjacent active faults using high-resolution 3D seismic data from the Smørbukk and Smørbukk South hydrocarbon fields, Halten Terrace, Mid-Norway. The major fault dividing the two fields, the Trestakk–Smørbukk fault, evolves from a southern segment with a well-defined set of rift wedges in its hangingwall to a northern segment where the fault tip is buried and a fault-tip fold is developed. Isochore maps of three Jurassic intervals illustrate a south to north evolution where, initially, Early Jurassic fault activity is limited to the southern part of the study area. Middle to Upper Jurassic intervals display a northwards migration in activity and linkage with two other major faults in the study area. This northwards migration had a profound effect on sediment geometries and depocentres in an area where previously only Late Jurassic rift activity has been recognized.     

Earthquake deformations in the Lisan deposits and seismotectonic implications

Summary. Earthquake deformations and induced sedimentary structures preserved in Quaternary sediments include faults, folds, fissures, slumps, sand boils and other effects of liquefaction. Such deformations and structures are well preserved in the Lisan deposits of the Dead Sea. Of most importance are the fold-type deformations known as décollement structures which are present all along the eastern side of the Lisan and seem to decrease gradually westwards to disappear in the middle of the Lisan. These may indicate that palaeoearthquakes originating along the Araba fault have triggered such structures due to shaking of elastoplastic unconsolidated sediments over gentle slopes dipping to the west.
Preliminary results from studies on décollement structures preserved in a section representing some 1733 years of continuous deposition in the uppermost? Pleistocene, in the vicinity of Wadi Araba, indicate that: (1) seismic activity has fluctuated with time. Average recurrence period is about 340 ± 20yr for earthquakes with magnitudes greater than or equal to 6.5, Earthquakes with magnitude greater than 7 seem to have occurred along the Araba fault. (2) Deduced earthquake magnitudes conform to the frequency–magnitude relationship: log N = 5.24–0.68 M . (3) The deduced seismic slip rate along the Araba fault seems to be not less than 0.64 ± 0.04 cm yr⁻¹.     

Late Quaternary systematic stream offsets caused by repeated large seismic events along the Kunlun fault, northern Tibet

The Kunlun fault is one of the largest strike-slip faults in northern Tibet, China. In this paper, we focus upon the Kusai Lake–Kunlun Pass segment of the fault to understand the geomorphic development of offset streams caused by repeated large seismic events, based on tectono-geomorphic analysis of high-resolution satellite remote sensing images combined with field studies. The results indicate that systematic left-lateral stream offsets appear at various scales across the fault zone: Lateral offsets of small gullies caused by the 2001 M_w 7.8 Kunlun earthquake vary typically from 3 m to 6 m, meanwhile streams with cumulative offsets of 10 m, 25–30 m, 50–70 m, 250–300 m and 750–1400 m have resulted from repeated large seismic events during the late Quaternary. An average slip rate of 10 ± 1 mm/year has been estimated from the lateral stream offsets and ¹⁴C ages of alluvial fan surfaces incised by the streams. A three-dimensional model showing tectono-geomorphic features along a left-lateral strike-slip fault is also presented. The Kusai Lake–Kunlun Pass segment provides an opportunity to understand the relationship between geomorphic features produced by individual large seismic events and long-term geomorphic development caused by repeated large seismic events along a major strike-slip fault.     

Relocation of M ≥ 2 events of the 1989 Dôbi seismic sequence in Afar: evidence for earthquake migration

The Ethiopian side of central Afar was struck in August 1989 by the largest seismic sequence (three 6.1 ≤ M _s ≤ 6.3 events, 15 with M _s or m _b ≥ 5.0) since that of Serdo in 1969. Using the Djibouti seismological network, we relocated 297 of the events of that sequence. As most of the large events took place outside the network, we assessed the accuracy and stability of earthquake relocations by using three different velocity models and two relocation codes to try to relate individual shocks to distinct faults and surface breaks. A majority of the events apparently occurred underneath the floor of the Dôbi graben, an area about 45  km long and 15  km wide, rupturing boundary and inner floor faults, in agreement with the surface cracks and scarps that we mapped in the area. The relocation shows that the principal events propagated about 50  km northwestwards along the graben in the first 40  hr. A day and a half after the beginning of the sequence, smaller events ( M ≤ 4) started to propagate more than 55  km eastwards, towards Asal Lake. Using two three-component stations installed near the Ethiopian border, we could determine reliable depths for 21 events. The depths are compatible with a seismogenic crust about 14  km thick in the Dôbi and Hanle graben area. Although the Dôbi sequence ruptured about 50  km of the fault array extending from Serdo to Asal, where the regional stress was released by earthquakes in 1969 and 1978, respectively, a seismic gap about 50  km long still subsists along the northern part of the Gaggade region (Der'êla half-graben).     

On the use of satellite altimetry to infer the earthquake rupture characteristics: application to the 2004 Sumatra event

Standard data and methods, such as the inversion of seismic and GPS data, have been used extensively to infer the details of the 2004 December 26 earthquake. The unprecedented large size of this event gave the opportunity to modern altimeters to provide the first clear records of a tsunami in deep ocean, therefore allowing us to study the rupture history from an independent perspective. We invert the Jason-1 and Topex–Poseidon altimetry records, considering the new constraints available on the geometry of the fault plane, and taking them into account in a 3-D rupture model. The data are corrected for the non-negligible effect of satellite motion during measurements. Our results show that the rupture propagated over the 1500 km of subduction zone initially identified by the aftershock distribution, with a magnitude of   M _w= 9.1  . Our solution compares well with the latitudinal distribution of slip inferred from other data sets, with a maximum of energy release north of Sumatra, and two other slip patches near the Nicobar and Andaman islands. Based on waveform comparison, we assert that the shallow portion of the megathrust offshore Banda Aceh had slip amplitudes of more than 20 m. Also, we find that significant amounts of slip (about 10 m) concentrated below the Andaman islands and did not propagate on the shallow portion of the interface. Although synthetic tests tend to show less resolution in the northern part of the rupture, this solution is compatible with the near-field data (GPS, coral heads and imagery), and would allow one to explain the apparent paradox between the large local displacements and the moderate tsunami observed locally. Finally, we demonstrate the rapidly dominating effect of propagation and slip distribution over the rupture velocity, and how it precludes the direct estimate of this latter parameter.     

Source history of the 1905 great Mongolian earthquakes (Tsetserleg, Bolnay)

Two great Mongolian earthquakes, Tsetserleg and Bolnay, occurred on 1905 July 9 and 23. We determined the source history of these events using body waveform inversion. The Tsetserleg rupture (azimuth N60°) correspond to a N60° oriented branch of the long EW oriented Bolnay fault.
Historical seismograms recorded by Wiechert instruments are digitized and corrected for the geometrical deformation due to the recording system. We use predictive filters to recover the signals lost at the minute marks.
The total rupture length for the Tsetserleg earthquake may reach up to 190 km, in order to explain the width of the recorded body waves. This implies adding 60 km to the previously mapped fault. The rupture propagation is mainly eastward. It starts at the southwest of the central subsegment, showing a left lateral strike-slip with a reverse component. The total duration of the modelled source function is 65 s. The seismic moment deduced from the inversion is 10²¹ N m, giving a magnitude   M _w = 8  .
The nucleation of the Bolnay earthquake was at the intersection between the main fault (375 km left lateral strike-slip) and the Teregtiin fault (N160°, 80 km long right lateral strike-slip with a vertical component near the main fault). The rupture was bilateral along the main fault: 100 km to the west and 275 km to east. It also propagated 80 km to the southeast along the Teregtiin fault. The source duration was 115 s. The moment magnitude Mw varies between 8.3 and 8.5.
The nucleation and rupture depths remain uncertain. We tested three cases: (1) nucleation and rupture depth limited to the seismogenic zone; (2) nucleation in the seismogenic zone and rupture propagation going to the base of the crust and (3) nucleation within the crust–upper mantle interface and rupture propagation within the upper mantle.