Frictional heating on a fault zone with finite thickness

Summary. The problem of frictional heating on a fault of finite thickness is considered. Temperature distributions during and after faulting are obtained for faults of various thicknesses. A thick fault approximation is shown to be valid if a fault has a thickness greater than about 1 cm. The thickness of the melted zone is predicted for various frictional stress levels and displacements. The predicted thicknesses are shown to be in reasonable agreement with field observations.     

Authors' reply to a comment on 'An objective method for smoothing palaeomagnetic data'

Summary. King & Rees (1979), in commenting on our paper (Clark & Thompson 1978), have raised a number of subjects beyond the main theme of our paper. We described a statistical method of smoothing noisy data and illustrated its application using a set of previously published palaeomagnetic declination data. King & Rees acknowledge our success in achieving the main aim of our paper, namely the construction of objective confidence limits for the smooth curve presumed to underlie the data. However, they seem to question the value of our general technique of data-analysis when used to estimate past variations in the geomagnetic field. As they interpose in their comments on our statistical analyses references to 'techniques of sedimentary analysis', 'inclination errors', 'mode of acquisition of remanence' and 'stability tests', we take the opportunity to discuss these additional topics as well as those relating to our original paper.     

Simulation of seismicity due to fluid migration in a fault zone

Spatio-temporal variation of rupture activity is modelled assuming fluid migration in a narrow, porous fault zone formed along a vertical strike-slip fault in a semi-infinite elastic medium. The principle of the effective stress coupled to the Coulomb failure criterion introduces mechanical coupling between fault slip and the pore fluid. The fluid is assumed to flow out of a localized high-pressure fluid compartment in the fault at the onset of earthquake rupture. The duration of the earthquake sequence is assumed to be much shorter than the recurrence period of characteristic events on the fault. Both an earthquake swarm and a foreshock–main-shock sequence can be simulated by changing the relative magnitudes of the initial tectonic stress, pore fluid pressure, fracture strength and so on. When an inhomogeneity is introduced into the spatial distribution of fracture strength, high complexity is observed in the spatio-temporal variation of rupture activity. For example, the time interval between two successive events is highly irregular, and a relatively long quiescence of activity is sometimes observed in a foreshock–main-shock sequence. The quiescence is caused by the temporary arresting of rupture extension, due to an encounter with fault segments having locally high strengths. The frequency–magnitude statistics of intermediate-size events obey the Gutenberg–Richter relation. The calculations show the temporal variation of the b value during some foreshock sequences, and the degree of the change seems to depend on the statistical distribution of the fracture strength.     

Detailed 3D seismic imaging of a fault zone beneath Lake Geneva, Switzerland

An efficient high‐resolution (HR) three‐dimensional (3D) seismic reflection system for small‐scale targets in lacustrine settings was developed. In Lake Geneva, near the city of Lausanne, Switzerland, the offshore extension of a complex fault zone well mapped on land was chosen for testing our system. A preliminary two‐dimensional seismic survey indicated structures that include a thin (<40 m) layer of subhorizontal Quaternary sediments that unconformably overlie south‐east‐dipping Tertiary Molasse beds and a major fault zone (Paudèze Fault Zone) that separates Plateau and Subalpine Molasse (SM) units. A 3D survey was conducted over this test site using a newly developed three‐streamer system. It provided high‐quality data with a penetration to depths of 300 m below the water bottom of non‐aliased signal for dips up to 30° and with a maximum vertical resolution of 1.1 m. The data were subjected to a conventional 3D processing sequence that included post‐stack time migration. Tests with 3D pre‐stack depth migration showed that such techniques can be applied to HR seismic surveys. Delineation of several horizons and fault surfaces reveals the potential for small‐scale geologic and tectonic interpretation in three dimensions. Five major seismic facies and their detailed 3D geometries can be distinguished. Three fault surfaces and the top of a molasse surface were mapped in 3D. Analysis of the geometry of these surfaces and their relative orientation suggests that pre‐existing structures within the Plateau Molasse (PM) unit influenced later faulting between the Plateau and SM. In particular, a change in strike of the PM bed dip may indicate a fold formed by a regional stress regime, the orientation of which was different from the one responsible for the creation of the Paudèze Fault Zone. This structure might have later influenced the local stress regime and caused the curved shape of the Paudèze Fault in our surveyed area.     

Quaternary landscape evolution in the San Jacinto fault zone,Peninsular Ranges of Southern California: Transient response to strike-slip fault initiation

Well-constrained case studies of transient landscape response to external forcing are needed to improve our understanding of erosion processes in tectonically active mountain belts. The Peninsular Ranges portion of the San Jacinto fault zone (SJFZ) is an excellent location for such a study because it displays pronounced geomorphic disequilibrium resulting from initiation of a major strike-slip fault in the past 1.0 to 2.5 million years. We recognize two geomorphic domains in this region: (1) a relict low-relief upland domain consisting of broad flat valleys and low-gradient streams and (2) very steep, rough topography with deeply incised canyons and retreating erosional knickpoints. Pleistocene sediments exposed along and near the SJFZ include fluvial conglomerate, sandstone, and mudstone, with weak paleosols and west- to NW-directed paleocurrents. These sediments accumulated in a low-gradient stream system (represented by domain 1) during an early phase of slip in the SJFZ, prior to the modern phase of erosion and degradation (domain 2). Late Pliocene or early Pleistocene initiation of the SJFZ triggered a wave of headward erosion and stream capture that is still migrating NW along the fault zone. Using the total distance that capture points have migrated along the fault zone and a range of possible ages for fault initiation, the rate of knickpoint retreat is estimated at ∼ 12 to 44 km/my.To explore the signal of transient geomorphic response to fault initiation, we analyzed 23 tributaries along an ∼ 20-km portion of the main fault valley within domain 2. The analysis reveals three zones with distinctive morphologies: (1) strongly convex longitudinal profiles in the NW, (2) a large (ca. 5–6 km²) landslide in the central zone, and (3) concave tributaries in the SE with profile complexity decreasing and catchment area increasing from NW to SE. The distribution of these zones suggests close spatial and temporal association of active fault slip, bedrock incision, deep-seated landslides, and erosional modification. The fundamental driving force behind these processes is profound geomorphic disequilibrium resulting from initiation of the SJFZ. We suggest that landslides may have played a significant role in shaping the morphology of this fault zone, and that the influence of landslides may be underestimated in areas where characteristic landforms and deposits are obscured by later erosion and faulting.     

Reply to comment by A. Douglas on 'Systematic determination of earthquake rupture directivity and fault planes from analysis of long-period P-wave spectra'

We respond to the comments by Douglas regarding our earlier paper by emphasizing that our automated method was intended to distinguish between the primary and auxiliary fault planes in earthquake focal mechanisms and does not always produce reliable results for rupture velocity and rupture length.     

Relay zone evolution: a history of repeated fault propagation and linkage,central Corinth rift,Greece

Established models indicate that, before being breached, relay zones along rift borders can evolve either by lengthening and rotating during progressive overlap of growing fault segments (isolated fault model), or, by simply rotating without lengthening before breaching (coherent fault model). The spatio‐temporal distribution of vertical motions in a relay zone can thus be used to distinguish fault growth mechanisms. Depositional relay zones that develop at sea level and accommodate both deposition on the ramp itself as well as transfer of sediments from the uplifting footwall into the hangingwall depocentres and provide the most complete record of vertical motions. We examine the development of a depositional relay ramp on the border of the active Corinth rift, Greece to reconstruct fault interaction in time and space using both onshore and offshore (2D seismic lines) data. The Akrata relay zone developed over a period of ca. 0.5 Myr since the Middle Pleistocene between the newly forming East Helike Fault (EHF) that propagated towards the older, more established Derveni Fault (DF). The relay zone captured the Krathis River, which deposited prograding Gilbert‐type deltas on the sub‐horizontal ramp. Successive oblique faults record progressive linkage and basinward migration of accommodation along the ramp axis, whereas marine terraces record diachronous uplift in their footwalls. Although early linkage of the relay zone occurs, continuous propagation and linkage of the EHF onto the static DF is recorded before final beaching. Rotation on forced folds above the upward and laterally propagating normal faults at the borders of the relay zone represents the ramp hinges. The Akrata relay zone cannot be compared directly to a simple fault growth model because (1) the relay zone connects two fault segments of different generations; (2) multiple linkages during propagation was facilitated by the presence of pre‐existing crustal structures, inherited from the Hellenide fold and thrust belt. The linkage of the EHF to the DF contributed to the westward and northward propagation of the southern rift border.     

Earthquake sequences on a frictional fault model with non-uniform strengths and relaxation times

Summary. The space and time characteristics of earthquake sequences, including a main shock, aftershocks and the recurrence of major shocks in a long time range, are investigated on a frictional fault model with non-uniform strengths and relaxation times, which is subjected to a time-dependent shear stress. Aftershocks with low stress drop take place successively in spaced regions so as to fill the gaps which have not yet been ruptured since the main shock, while those with high stress drop occur in and around the regions left unruptured during the main faulting. The frequency decay of aftershocks with time follows a hyperbolic law with the rates p consistent with observations. There are good linear relations in logarithmic scales for source area versus frequency and seismic moment versus frequency of the generated aftershocks. The b -value obtained in the present experiments appears slightly larger than that for observations. It was found that more heterogeneous distribution of the fault strength give smaller p and larger b -values. The recurrence of major shocks, particularly of very large shocks with high stress drop, is often preceded by a completely silent period of activity or very low activity with a small number of foreshocks. The major shocks take place successively in adjacent unruptured regions and sometimes show slow-speed migrations. These results provide explanations to various observations of earthquake sequences.