Continuous GPS measurements across the Western Alps, 1996–1998

Global Positioning System (GPS) data are analysed from three permanent sites, two in the Western Alps (Grasse, France, and Zimmerwald, Switzerland) and the third in the Po Basin (Torino, Italy), for the 2.5-year period from 1 January 1996 to 1 July 1998. An analysis of the stochastic properties of the position time series reveals a significant amount of spatially and temporally correlated noise, which best fits a model combining white noise and flicker noise. The coloured noise is drastically reduced by spatially filtering the time series, suggesting that it is not due to site-specific effects such as monument motion, but rather to noise sources common to the three sites, such as reference frame, satellite orbit or Earth orientation errors. We find velocity uncertainties (95 per cent confidence interval) of 2.7  mm  yr⁻¹ at GRAS, 2.3  mm  yr⁻¹ at ZIMM, and 2.9  mm  yr⁻¹ at TORI. The residual velocity in the Eurasian reference frame is statistically greater than zero at ZIMM only, with 2.4 ± 2.3  mm  yr⁻¹ of motion in a NW direction. These uncertainties place an upper bound on the expected deformation in the Western Alps. Velocities relative to stable Eurasia do not exceed 3  mm  yr⁻¹, corresponding to a maximum strain rate of 0.03  μstrain  yr⁻¹ over the GRAS–TORI–ZIMM triangle.     

Seismic anisotropy within the uppermost mantle of southern Germany

This paper presents an updated interpretation of seismic anisotropy within the uppermost mantle of southern Germany. The dense network of reversed and crossing refraction profiles in this area made it possible to observe almost 900 traveltimes of the P_n phase that could be effectively used in a time-term analysis to determine horizontal velocity distribution immediately below the Moho. For 12 crossing profiles, amplitude ratios of the P_n phase compared to the dominant crustal phase were utilized to resolve azimuthally dependent velocity gradients with depth. A P -wave anisotropy of 3–4 per cent in a horizontal plane immediately below the Moho at a depth of 30 km, increasing to 11 per cent at a depth of 40 km, was determined. For the axis of the highest velocity of about 8.03 km s⁻¹ at a depth of 30 km a direction of N31°F was obtained. The azimuthal dependence of the observed P_n amplitude is explained by an azimuth-dependent sub-Moho velocity gradient decreasing from 0.06 s⁻¹ in the fast direction to 0 s⁻¹ in the slow direction of horizontal P -wave velocity. From the seismic results in this study a petrological model suggesting a change of modal composition and percentage of oriented olivine with depth was derived.     

On the consistency of earthquake moment release and space geodetic strain rates: Europe

In this paper, approximately 100 VLBI/SLR/GPS velocities map European strain rates from <0.09 × 10⁻⁸ to >9.0 × 10⁻⁸ yr⁻¹ with regional uncertainties of 20 to 40 per cent. Kostrov's formula translates these strain-rate values into regional geodetic moment rates M¯˙ _geodetic . Two other moment rates, M¯˙ _seismic , extracted from a 100-year historical catalogue and M¯˙ _plate , taken from plate-tectonic models, contrast the geodetic rates. In Mediterranean Europe, the ratios of M¯˙ _seismic to M¯˙ _geodetic are between 0.50 and 0.71. In Turkey the ratio falls to 0.22. Although aseismic deformation may contribute to the earthquake deficit ( M¯˙ _seismic values less than M¯˙ _geodetic ), the evidence is not compelling because the magnitudes of the observed shortfalls coincide with the random variations expected in a 100-year catalogue. If the lack of aseismic deformation inferred from the 100-year catalogue holds true for longer periods, then much of Europe's strain budget would have to be accommodated by more frequent or larger earthquakes than have been experienced this century to raise the ratios of M¯˙ _seismic to M¯˙ _geodetic to unity. Improved geological fault data bases, longer historical earthquake catalogues, and densification of the continent's space geodetic network will clarify the roles of aseismic deformation versus statistical quiescence.     

Geodetically determined strain across the southern end of the Tonga-Kermadec-Hikurangi subduction zone

Summary. A method of simultaneous reduction is presented for determining strain rates from multiple triangulation surveys where common triangulation stations have been used, but the angles of the old survey have not necessarily been reobserved. This method is applied to triangulation in the northern South Island, at the southern end of the Tonga-Kermadec-Hikurangi subduction zone. From a profile of shear strain across the Indian-Pacific plate boundary, the displacement of the Indian plate relative to the Pacific is calculated to be 54 ± 9 mm yr⁻¹ at an azimuth of 84°± 10°, in remarkable agreement with the motion predicted by global plate tectonic models. Most of this motion occurs within a 150 km wide zone bounded on the east by the Hikurangi Trough. Within this zone the motion is partitioned: near the Hikurangi Trough no slip is occurring at the upper surface of the subducting Pacific plate (the subduction thrust) and motion is predominantly thrusting normal to the trough axis: to the west is a region of predominantly dextral strike slip faulting. This pattern is consistent with Fitch's model of oblique subduction. To the south of the profile, a change is observed in the azimuth of the faulting along a line which marks the southern extent of the subduction slab, indicating the end of the partitioned motion.     

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.     

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.     

Fault propagation and climatic control of sedimentation on the Ghoubbet Rift Floor: insights from the Tadjouraden cruise in the western Gulf of Aden

A detailed geophysical survey of the Ghoubbet Al Kharab (Djibouti) clarifies the small-scale morphology of the last submerged rift segment of the propagating Aden ridge before it enters the Afar depression. The bathymetry reveals a system of antithetic normal faults striking N130°E, roughly aligned with those active along the Asal rift. The 3.5 kHz sub-bottom profiler shows how the faults cut distinct layers within the recent, up to 60 m thick, sediment cover on the floor of the basin. A large volcanic structure, in the centre of the basin, the 'Ghoubbet' volcano, separates two sedimentary flats. The organization of volcanism and the planform of faulting, with en echelon subrifts along the entire Asal–Ghoubbet rift, appear to confirm the westward propagation of this segment of the plate boundary. Faults throughout the rift have been active continuously for the last 8400 yr, but certain sediment layers show different offsets. The varying offsets of these layers, dated from cores previously retrieved in the southern basin, imply Holocene vertical slip rates of 0.3–1.4 mm yr⁻¹ and indicate a major decrease in sedimentation rate after about 6000 yr BP, and a redistribution of sediments in the deepest troughs during the period that preceded that change.     

Active tectonics of an apparently aseismic region: distributed active strike-slip faulting in the Hangay Mountains of central Mongolia

We identify and describe a series of east–west left-lateral strike-slip faults (named the Songino-Margats, the Hag Nuur, the Uliastay and the South Hangay fault systems) in the Hangay mountains of central Mongolia: an area that has little in the way of recorded seismicity and which is often considered as a rigid block within the India–Eurasia collision zone. The strike-slip faults of central Mongolia constitute a previously unrecognized hazard in this part of Mongolia. Each of the strike-slip faults show indications of late Quaternary activity in the form of aligned sequences of sag-ponds and pressure-ridges developed in alluvial deposits. Total bed-rock displacements of ∼3 km are measured on both the Songino-Margats and South Hangay fault systems. Bed-rock displacements of 11 km are observed across the Hag Nuur fault. Cumulative offset across the Uliastay fault systems are unknown but are unlikely to be large. We have no quantitative constraint on the age of faulting in the Hangay. The ≤20 km of cumulative slip on the Hangay faults might, at least in part, be inherited from earlier tectonic movements. Our observations show that, despite the almost complete absence of instrumentally recorded seismicity in the Hangay, this part of Mongolia is cut through by numerous distributed strike-slip faults that accommodate regional left-lateral shear between Siberia and China. Central Mongolia is thus an important component of the India–Eurasia collision that would be overlooked in models of the active tectonics based on the distribution of seismicity. We suggest that active faults such as those identified in the Hangay of Mongolia might exist in other, apparently aseismic, regions within continental collision zones.     

2001 August earthquake swarm at Shadwan Island, Gulf of Suez, Egypt

The earthquake swarm that struck Shadwan Island at the entrance of the Gulf of Suez in 2001 August included 408 events. Almost all of these events (94 per cent) were microearthquakes and only 6 per cent had small measurable magnitudes  (5.0 > M _L≥ 3.0)  . Most of the earthquakes were weak and followed each other so closely in time that they could not be identified at more distant stations. The fault plane solutions of the strongest events of the swarm show almost identical focal mechanisms, predominately normal faulting with a significant sinistral strike-slip component for nodal planes trending NW–SE. A comparison with the mechanisms of the 1969 and 1972 events which took place 20 km north of the swarm epicentral region shows similarities in faulting type and orientation of nodal planes. The azimuths of T -axes determined from focal mechanisms in this study are oriented in the NNE–SSW direction. This direction is consistent with the present-day stress field derived from borehole breakouts in the southern Gulf of Suez and the last phase of stress field changes in the Late Pleistocene, as well as with recent GPS results.
The source parameters of the largest  ( M _L≥ 3.0)  events of the 2001 August Shadwan swarm have been estimated from the P -wave spectra of the Egyptian National Seismograph Network (ENSN). Averaging of the values obtained at different stations shows relatively similar source parameters, including a fault length of  0.65 ≤ L ≤ 2 km  , a seismic moment of  7.1 × 10¹²≤ Mo ≤ 3.0 × 10¹⁴ N m  and a stress drop of  0.4 ≤Δσ≤ 10  bar.     

Strain accumulation along an oblique plate boundary: the Reykjanes Peninsula, southwest Iceland

We use annual GPS observations on the Reykjanes Peninsula (RP) from 2000 to 2006 to generate maps of surface velocities and strain rates across the active plate boundary. We find that the surface deformation on the RP is consistent with oblique plate boundary motion on a regional scale, although considerable temporal and spatial strain rate variations are observed within the plate boundary zone. A small, but consistent increase in eastward velocity is observed at several stations on the southern part of the peninsula, compared to the 1993–1998 time period. The 2000–2006 velocities can be modelled by approximating the plate boundary as a series of vertical dislocations with left-lateral motion and opening. For the RP plate boundary we estimate left-lateral motion  18⁺⁴₋₃ mm yr⁻¹  and opening of  7⁺³₋₂ mm yr⁻¹  below a locking depth of  7⁺¹₋₂ km  . The resulting deep motion of  20⁺⁴₋₃ mm yr⁻¹  in the direction of  N(100⁺⁸₋₆)°E  agrees well with the predicted relative North America–Eurasia rate. We calculate the areal and shear strain rates using velocities from two periods: 1993–1998 and 2000–2006. The deep motion along the plate boundary results in left-lateral shear strain rates, which are perturbed by shallow deformation due to the 1994–1998 inflation and elevated seismicity in the Hengill–Hrómundartindur volcanic system, geothermal fluid extraction at the Svartsengi power plant, and possibly earthquake activity on the central part of the peninsula.     

Large-scale velocity field and strain tensor in Iran inferred from GPS measurements: new insight for the present-day deformation pattern within NE Iran

A network of 26 GPS sites was implemented in Iran and Northern Oman to measure displacements in this part of the Arabia–Eurasia collision zone. We present the GPS velocity field obtained from three surveys performed in 1999 September, 2001 October and 2005 September and the deduced strain tensor. This study refines previous studies inferred from only the two first surveys. Improvements are significant in NE Iran. The present-day shortening rate across the mountain belts of NE Iran is estimated to 5 ± 1 mm yr⁻¹ at about N11°, 2 ± 1 mm yr⁻¹ of NS shortening across the eastern Kopet Dag and 3 ± 1 mm yr⁻¹ of NS shortening across Binalud and Kuh-e-Sorkh. Our GPS measurements emphasize the varying character of the Kopet Dag deformation between its southeastern part with prevailing thrusting at low rates and its northwestern part with dominant strike-slip activity at increasing rates. The principal axes of the horizontal strain tensor appears very homogeneous from the Zagros to the Alborz and the Kopet-Dag (N20°) and in eastern Iran (Makran and Lut block: N30°). Only NW Iran suffers a variable strain pattern which seems to wrap the Caspian basin. The strain tensor map underlines the existence of large homogeneous tectonic provinces in terms of style and amplitude of the deformation.     

Rates of active deformation in the Aegean Sea and surrounding regions

Abstract Average strain rates are calculated from earthquakes in the period 1908-81 that occurred in the Aegean Sea extensional region, and in the convergent zone associated with the Hellenic Trench. In spite of large uncertainties resulting from the use of an M_S : M_o relationship, seismic N-S extensional rates in the Aegean are in the region 20–60 mm yr^-1 whereas seismic shortening rates in the Hellenic Trench are less than about 15 mm yr^-1. This is surprising because Africa and Eurasia are known to be converging, not separating. This apparent anomaly is caused by most of the convergence in the Hellenic Trench occurring aseismically. By contrast, the seismic extensional rates in the Aegean agree quite well with those expected from other arguments. The present day extensional rates are sufficiently high for McKenzie's instantaneous stretching model to be applicable. There is some evidence that these high extensional rates have operated throughout the last 5 Myr.     

Monitoring overburden layer changes and fault movements from time-lapse seismic data on the Valhall Field

A method developed by Røste et al. , which discriminates between layer thickness and velocity changes, is tested on pre-stack time-lapse seismic ocean bottom cable (OBC) data from the Valhall Field. A key parameter in this discrimination process is the dilation factor, α, which is the relative velocity change divided by the relative thickness change within a given layer. The high quality and good repeatability of the OBC data enables us to estimate α with a reasonable accuracy for α-values between 0 and −5. For α-values below −5, complementary information, like for instance geomechanics is required. For the top reservoir horizon we estimate a maximum subsidence of 0.50 m ± 0.29 m and a corresponding velocity decrease for the sequence from the seabed to the top reservoir of 2.0 m s⁻¹± 0.40 m s⁻¹.
Time-lapse distortion patterns with characteristic time-shift versus offset signatures are observed. The positions and vertical extents of the distortion zones are determined from ray path considerations and modelling. The distortion zones correlate with buried faults, indicating that a (time-lapse) distortion zone might be produced by a localized slip in a fault zone. We present an extended method which allows for vertical (in addition to lateral) variations in the relative thickness and velocity changes. This method can be viewed as a simplified version of time-lapse tomography, but involving fewer unknown parameters, giving more stability to the estimated changes in thickness and velocity. Using this technique, we are able to estimate α for positions with localized time-lapse distortions.     

Stacking for the frequencies and Qs of 0S0 and 1S0

Summary. Using nine IDA records for the Indonesian earthquake of 1977 August 19, we have formed an optimal linear combination of the records and have measured the frequency and Q of ₀ S ₀ and ₁ S ₀. The frequency was measured using the moment ratio method. The attenuation was measured by the minimum width method and by the time-lapse method. The frequency and attenuation were measured simultaneously by varying them to obtain a best fit to the data. A 2000-hr stack, the sum of nine individual records, for ₀ S ₀ gave a frequency of 0.814664 mHz±4 ppm. The values for the Q of ₀ S ₀ for the three different methods of measurement were 5600,5833 and 5700, respectively. The error in the estimates of Q ^-1 is about 5 per cent for the minimum power method. For ₁ S ₀ a 300-hr stack yielded a frequency of 1.63151 mHz±30 ppm. The values of Q for this mode were 1960, 1800 and 1850, respectively, with an error in Q ^-1 of about 12 per cent for the minimum power method.     

Dynamical rupture process on a three-dimensional fault with non-uniform frictions and near-field seismic waves

Summary. Dynamical rupture process on the fault is investigated in a quasi-three-dimensional faulting model with non-uniform distributions of static frictions or the fracture strength under a finite shearing pre-stress. The displacement and stress time functions on the fault are obtained by solving numerically the equations of motion with a finite stress—fracture criterion, using the finite difference method.
If static frictions are homogeneous or weakly non-uniform, the rupture propagates nearly elliptically with a velocity close to that of P waves along the direction of pre-stress and with a nearly S wave velocity in the direction perpendicular to it. The rise time of the source function and the final displacements are larger around the centre of the fault. In the case when the static frictions are heavily non-uniform and depend on the location, the rupture propagation becomes quite irregular with appreciably decreased velocities, indicating remarkable stick-slip phenomena. In some cases, there remain unruptured regions where fault slip does not take place, and high stresses remain concentrated up to the final stage. These regions could be the source of aftershocks at a next stage.
The stick—slip faulting and irregular rupture propagation radiate high-frequency seismic waves, and the near-field spectral amplitudes tend to show an inversely linear frequency dependence over high frequencies for heavily non-uniform frictional faults.     

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.     

Experimental study of the anisotropy of longitudinal and transverse waves from local earthquake records

Summary. The paper gives the results of a study of the anisotropy of seismic wave velocities within the Ashkhabad test field in Central Asia. The anisotropy was studied by analysing variations in the values of apparent velocities of first arrivals for epicentral distances ranging from 30 to 130 km and by analysing the delays (Δ t_s1-_s2 ) between the arrival times of shear waves with different polarizations.
The velocities of P -waves vary with azimuth from 5.3 to 6.27 km s^-1 and the velocities of S -waves vary from 3.15 to 3.5 km s^-1.
The delay times Δ t_S1 - _S2 depend on the direction of the propagation. The character of the variation of the propagation velocity of the longitudinal wave, the presence of two differently polarized shear waves S 1 and S 2 propagating at different velocities, and the character of the distribution of Δ t_S1 - _S2 on the stereogram suggest that the symmetry of the anisotropic medium is close to hexagonal with a nearly horizontal symmetry axis coinciding with the direction of maximal velocity. The azimuth of the symmetry axis of the medium is 140° and coincides with the direction of geological faults.     

Uplift and heat flow following the injection of magmas into the lithosphere

Summary. The thermal effect of a rapid injection of hot magmas into the lower part of the lithosphere is modelled as an increase in heat production through the invaded region. The change in surface heat flow and the uplift resulting from the thermal expansion are determined in three-dimensional axially symmetric geometry: they are expressed as the space time convolutions of a Green's function with the anomalous heat production.
The anomalies with shorter wavelength (compared to the lithospheric thickness) are attenuated. This filtering affects the surface uplift more than the heat flow anomaly; the attenuation effect is larger when only the lower part of the lithosphere is invaded.
The uplift time constant is of the same order as the heat conduction time if the lower lithosphere is invaded by magmas at a moderate rate (i.e. the rate of injection does not exceed the equivalent of 0.1 per cent of the lithospheric volume in 10⁶yr). Fifty per cent of the total uplift takes place in about 80 × 10⁶yr for a lithosphere 100 km thick. The uplift is slightly faster when the whole lithosphere is invaded. The heat flow anomaly is delayed when the lower part of the lithosphere is invaded.
The spatial extent and the timing of the uplift and heat flow anomalies are critical in determining the mechanism's feasibility. Magma injections explain rapid uplifts [> 100 m (10⁶ yr)⁻¹] only if the magma is supplied at a very high rate (i.e. at least 10 per cent of the lithosphere volume per 10⁶yr). It is a feasible mechanism for uplifts that occur over longer periods of time (≊ 30 × 10⁶yr) such as those that seem to have occurred when the African plate came to rest with respect to the mantle.     

On the quality of marine magnetic anomaly sources and sea-floor topography

Summary. The quality of marine magnetic anomaly sources is described with the power-density representation of a stochastic model of random temporal and spatial emplacement of the marine magnetic anomaly source in the oceanic crust. Typical values of sea-floor spreading and emplacement parameters define a high-fidelity process of recording and sea-surface detection of the palaeomagnetic field reversals for spreading rates over 20 mm yr^-1. An analogous stochastic model is developed for the formation of sea-floor topography by normal faulting. It is shown that the random process of normal faulting observed in the inner walls of the FAMOUS rift valley can account for the quality of the adjacent West Rift Mountains topography.