Crustal structure of the Eastern Alps along the TRANSALP profile from wide-angle seismic tomography

The objective of the TRANSALP project is an investigation of the Eastern Alps with regard to their deep structure and dynamic evolution. The core of the project is a 340-km-long seismic profile at 12°E between Munich and Venice. This paper deals with the P-wave velocity distribution as derived from active source travel time tomography. Our database consists of Vibroseis and explosion seismic travel times recorded at up to 100 seismological stations distributed in a 30-km-wide corridor along the profile. In order to derive a velocity and reflector model, we simultaneously inverted refractions and reflections using a derivative of a damped least squares approach for local earthquake tomography. 8000 travel time picks from dense Vibroseis recordings provide the basis for high resolution in the upper crust. Explosion seismic wide-angle reflection travel times constrain both deeper crustal velocities and structure of the crust–mantle boundary with low resolution. In the resulting model, the Adriatic crust shows significantly higher P-wave velocities than the European crust. The European Moho is dipping south at an angle of 7°. The Adriatic Moho dips north with a gentle inclination at shallower depths. This geometry suggests S-directed subduction. Azimuthal variations of the first-break velocities as well as observations of shear wave splitting reveal strong anisotropy in the Tauern Window. We explain this finding by foliations and laminations generated by lateral extrusion. Based on the P-wave model we also localized almost 100 local earthquakes recorded during the 2-month acquisition campaign in 1999. Seismicity patterns in the North seem related to the Inn valley shear zone, and to thrusting of Austroalpine units over European basement. The alignment of deep seismicity in the Trento-Vicenza region with the top of the Adriatic lower crust corroborates the suggestion of a deep thrust fault in the Southern Alps.     

3-D seismic structure of the source area of the 1993 Latur, India, earthquake and its implications for rupture nucleations

The Latur earthquake (M_w 6.1) of 29 September 1993 is a rare stable continental region (SCR) earthquake that occurred on a previously unknown blind fault. In this study, we determined detailed three-dimensional (3-D) P- and S-wave velocity (V_p, V_s) and Poisson's ratio (σ) structures by inverting the first P- and S-wave high-quality arrival time data from 142 aftershocks that were recorded by a network of temporary seismic stations. The source zone of the Latur earthquake shows strong lateral heterogeneities in V_p, V_s and σ structures, extending in a volume of about 90 × 90 × 15 km³. The mainshock occurred within, but near the boundary, of a low-V_p, high-V_s and low-σ zone. This suggests that the structural asperities at the mainshock hypocenter are associated with a partially fluid-saturated fractured rock in a previously unknown source zone with intersecting fault surfaces. This might have triggered the 1993 Latur mainshock and its aftershock sequence. Our results are in good agreement with other geophysical studies that suggest high conductivity and high concentration of radiogenic helium gas beneath the source zone of the Latur earthquake. Our study provides an additional evidence for the presence of fluid related anomaly at the hidden source zone of the Latur earthquake in the SCR and helps us understand the genesis of damaging earthquakes in the SCR of the world.     

The Tornquist Zone, a north east inclining lithospheric transition at the south western margin of the Baltic Shield: Revealed through a nonlinear teleseismic tomographic inversion

From July 1996 to August 1997 the TOR project operated 130 seismographs in North Germany, Denmark and South Sweden, with the aim of collecting signals from local, regional and teleseismic earthquakes. This data set is particularly interesting since the seismic antenna crosses the most significant geological boundary in Europe, the Tornquist Zone, which in the northern part is the border between the Baltic Shield and the younger European lithosphere. Previous studies have shown significant physical changes in the crust and upper mantle across this transition zone, including two independent teleseismic tomographic studies of the TOR data set. But these two studies disagree on the orientation of the slope of the transition. Both studies used an iterative linearized inversion method. We will in this work Preprint submitted to Elsevier Science 27 July 2005 present an inversion based on Bayesian statistics, where the solution space is examined in order to study a very large number of tomographic solutions and to examine the solution uniqueness and uncertainty. The method is applied to measurements of 3345 relative teleseismic P-phase travel times from 48 teleseismic earthquakes with good azimuthal coverage with respect to the great circle arc of the TOR array. We find the lithospheric transition to be a north east inclination of around 30° to 45° off vertical.     

Geophysical methods applied to fault characterization and earthquake potential assessment in the Lower Tagus Valley, Portugal

The study region is located in the Lower Tagus Valley, central Portugal, and includes a large portion of the densely populated area of Lisbon. It is characterized by a moderate seismicity with a diffuse pattern, with historical earthquakes causing many casualties, serious damage and economic losses. Occurrence of earthquakes in the area indicates the presence of seismogenic structures at depth that are deficiently known due to a thick Cenozoic sedimentary cover. The hidden character of many of the faults in the Lower Tagus Valley requires the use of indirect methodologies for their study. This paper focuses on the application of high-resolution seismic reflection method for the detection of near-surface faulting on two major tectonic structures that are hidden under the recent alluvial cover of the Tagus Valley, and that have been recognized on deep oil-industry seismic reflection profiles and/or inferred from the surface geology. These are a WNW–ESE-trending fault zone located within the Lower Tagus Cenozoic basin, across the Tagus River estuary (Porto Alto fault), and a NNE–SSW-trending reverse fault zone that borders the Cenozoic Basin at the W (Vila Franca de Xira–Lisbon fault). Vertical electrical soundings were also acquired over the seismic profiles and the refraction interpretation of the reflection data was carried out. According to the interpretation of the collected data, a complex fault pattern disrupts the near surface (first 400 m) at Porto Alto, affecting the Upper Neogene and (at least for one fault) the Quaternary, with a normal offset component. The consistency with the previous oil-industry profiles interpretation supports the location and geometry of this fault zone. Concerning the second structure, two major faults were detected north of Vila Franca de Xira, supporting the extension of the Vila Franca de Xira–Lisbon fault zone northwards. One of these faults presents a reverse geometry apparently displacing Holocene alluvium. Vertical offsets of the Holocene sediments detected in the studied geophysical data of Porto Alto and Vila Franca de Xira–Lisbon faults imply minimum slip rates of 0.15–0.30 mm/year, three times larger than previously inferred for active faults in the Lower Tagus Valley and maximum estimates of average return periods of 2000–5000 years for M 6.5–7 co-seismic ruptures.     

Computation of Large Anisotropic Seismic Heterogeneities (CLASH)

A general tomographic technique is designed in order (i) to operate in anisotropic media; (ii) to account for the uneven seismic sampling and (iii) to handle massive data sets in a reasonable computing time. One modus operandi to compute a 3-D body wave velocity model relies on surface wave phase velocity measurements. An intermediate step, shared by other approaches, consists in translating, for each period of a given mode branch, the phase velocities integrated along ray paths into local velocity perturbations. To this end, we develop a method, which accounts for the azimuthal anisotropy in its comprehensive form. The weakly non-linear forward problem allows to use a conjugate gradient optimization. The Earth's surface is regularly discretized and the partial derivatives are assigned to the individual grid points. Possible lack of lateral resolution, due to the inescapable uneven ray path coverage, is taken into account through the a priori covariances on parameters with laterally variable correlation lengths. This method allows to efficiently separate the 2ψ and the 4ψ anisotropic effects from the isotropic perturbations. Fundamental mode and overtone phase velocity maps, derived with real Rayleigh wave data sets, are presented and compared with previous maps. The isotropic models concur well with the results of Trampert & Woodhouse. Large 4ψ heterogeneities are located in the tectonically active regions and over the continental lithospheres such as North America, Antarctica or Australia. At various periods, a significant 4ψ signature is correlated with the Hawaii hotspot track. Finally, concurring with the conclusions of Trampert & Woodhouse, our phase velocity maps show that Rayleigh wave data sets do need both 2ψ and 4ψ anisotropic terms.