Partial derivatives of surface wave phase velocities for flat anisotropic models

Summary. Surface wave behaviour in flat anisotropic structures is first illustrated by performing an exact computation on a simple two-layer model. The variational procedure of Smith & Dahlen is then used to compute the partial derivatives of surface wave phase velocities with respect to the elastic parameters in more realistic earth models. Linear relationships between the partial derivatives for a general anisotropic structure and those for a transversely isotropic structure are derived. When considering waves propagating in a fixed direction, there are only four independent derivatives for Rayleigh waves, and two for Love waves. To avoid the lack of resolution in an inverse method, we propose to use physically constrained models. These results are illustrated by using a model with hexagonal symmetry and a symmetry axis oriented either vertically or horizontally. Quasi-Love- and quasi-Rayleigh-wave partial derivatives are computed for both axis orientations. Modes up to the second overtone and periods ranging between 45 and 130 s have been considered. Finally, anomalies of phase velocity are computed in an oceanic model made of 1/6 oriented olivine crystals with horizontal or vertical preferred orientations of the a -axis.     

Inversion of multimode surface wave data: evidence for sub-lithospheric anisotropy

Summary. Multiparameter inversions of multimode dispersion data are performed for two large regions: the Pacific Ocean and North America. Anisotropy is taken into account by considering transversely isotropic structures with a vertical axis of symmetry. Two fundamental questions are studied in detail: (1) how to make the inverted models consistent when using different sets of parameters, (2) what is the significance of transversely isotropic inversion for the actual Earth's structure? It is proved that full consistency of the inverted models can be achieved by properly taking into account some a priori informations on the model and it is shown that the use of transversely isotropic models with vertical axis of symmetry does not cause severe limitations when interpreting the data. The models we have obtained are discussed in the light of these investigations. Considering an olivine-rich upper mantle, we make a tentative interpretation of these models in terms of preferred orientation of the a -axis of the crystals in one fixed horizontal direction.     

SANICLAY: simple anisotropic clay plasticity model

SANICLAY is a new simple anisotropic clay plasticity model that builds on a modification of an earlier model with an associated flow rule, in order to include simulations of softening response under undrained compression following K_o consolidation. Non‐associativity is introduced by adopting a yield surface different than the plastic potential surface. Besides, the isotropic hardening of the yield surface both surfaces evolve according to a combined distortional and rotational hardening rule, simulating the evolving anisotropy. Although built on the general premises of critical state soil mechanics, the model induces a critical state line in the void ratio–mean effective stress space, which is a function of anisotropy. To ease interpretation, the model formulation is presented firstly in the triaxial stress space and subsequently, its multiaxial generalization is developed systematically, in a form appropriate for implementation in numerical codes. The SANICLAY is shown to provide successful simulation of both undrained and drained rate‐independent behaviour of normally consolidated sensitive clays, and to a satisfactory degree of accuracy of overconsolidated clays. The new model requires merely three constants more than those of the modified Cam clay model, all of which are easily calibrated from well‐established laboratory tests following a meticulously presented procedure. Copyright © 2006 John Wiley & Sons, Ltd.     

Crustal motion and deformation in Greece from a decade of GPS measurements, 1993–2003

The Hellenic plate boundary region, located in the collision zone between the Nubian/Arabian and Eurasian lithospheric plates, is one of the seismo-tectonically most active areas of Europe. During the last 15 years, GPS measurements have been used to determine the crustal motion in the area of Greece with the aim to better understand the geodynamical processes of this region. An extended reoccupation network covering whole Greece has been measured periodically in numerous GPS campaigns since the late eighties, and a continuous GPS network has been operated in the region of the Ionian Sea since 1995. In this paper, we present a new detailed high-quality solution of continuous and campaign-type measurements acquired between 1993 and 2003. During the GPS processing, a special effort was made to obtain consistent results with highest possible accuracies and reliabilities. Data of 54 mainly European IGS and EUREF sites were included in the GPS processing in order to obtain results which are internally consistent with the European kinematic field and order to allow for a regional interpretation. After an overview of the results of the IGS/EUREF sites, the results from more than 80 stations in Greece are presented in terms of velocities, time series, trajectories and strain rates. Previous geodetic, geological and seismological findings are generally confirmed and substantially refined. New important results include the observation of deformation zones to the north and to the south of the North Aegean Trough and in the West Hellenic arc region, arc-parallel extension of about 19 mm/yr along the Hellenic arc, and compression between the Ionian islands and the Greek mainland. Due to continuous long-term observations of 4–8 years, it was possible to extract height changes from the GPS time series. In Greece, we observe a differential subsidence of the order of 2 mm/yr between the northern and central Ionian islands across the Kefalonia fault zone. The differential subsidence of the central Ionian islands with respect to the northwestern Greek mainland amounts to 4 mm/yr.     

Partial crustal melting beneath the Betic Cordillera (SE Spain): The case study of Mar Menor volcanic suite

The Neogene Volcanic Province (NVP) within the Betic Cordillera (SE Spain) consists of three main metapelitic enclave suites (from SW to NE: El Hoyazo, Mazarrón and Mar Menor). Since the NVP represents a singular place in the world where crustal enclaves were immediately quenched after melting, their microstructures provide a “photograph” of the conditions at depth just after the moment of the melting.

The thermobarometric information provided by the different microstructural assemblages has been integrated with the geophysical and geodynamical published data into a model of the petrologic evolution of the Mar Menor enclaves. They were equilibrated at 2–3 kbar, 850–900 °C, and followed a sequence of heating melt producing reactions. A local cooling event evidenced by minor melt crystallization preceded the eruption.

The lower crustal studies presented in this work contribute to the knowledge of: (i) the partial melting event beneath the Mar Menor volcanic suite through a petrologic detailed study of the enclaves; (ii) how the microstructures of fast cooled anatectic rocks play an important role in tracing the magma evolution in a chamber up to the eruption, and how they can be used as pseudothermobarometers; (iii) the past and current evolution of the Alborán Domain (Betic Cordillera) and Mediterranean Sea, and how the base of a metapelitic crust has melted within an active geodynamic setting.     

Crustal structure at the easternmost termination of the Variscan belt based on CELEBRATION 2000 and ALP 2002 data

The eastern margin of the Variscan belt in Europe comprises plate boundaries between continental blocks and terranes formed during different tectonic events. The crustal structure of that complicated area was studied using the data of the international refraction experiments CELEBRATION 2000 and ALP 2002. The seismic data were acquired along SW–NE oriented refraction and wide-angle reflection profiles CEL10 and ALP04 starting in the Eastern Alps, passing through the Moravo-Silesian zone of the Bohemian Massif and the Fore-Sudetic Monocline, and terminating in the TESZ in Poland. The data were interpreted by seismic tomographic inversion and by 2-D trial-and-error forward modelling of the P waves. Velocity models determine different types of the crust–mantle transition, reflecting variable crustal thickness and delimiting contacts of tectonic units in depth. In the Alpine area, few km thick LVZ with the Vp of 5.1 km s^− 1 dipping to the SW and outcropping at the surface represents the Molasse and Helvetic Flysch sediments overthrust by the Northern Calcareous Alps with higher velocities. In the Bohemian Massif, lower velocities in the range of 5.0–5.6 km s^− 1 down to a depth of 5 km might represent the SE termination of the Elbe Fault Zone. The Fore-Sudetic Monocline and the TESZ are covered by sediments with the velocities in the range of 3.6–5.5 km s^− 1 to the maximum depth of 15 km beneath the Mid-Polish Trough. The Moho in the Eastern Alps is dipping to the SW reaching the depth of 43–45 km. The lower crust at the eastern margin of the Bohemian Massif is characterized by elevated velocities and high Vp gradient, which seems to be a characteristic feature of the Moravo-Silesian. Slightly different properties in the Moravian and Silesian units might be attributed to varying distances of the profile from the Moldanubian Thrust front as well as a different type of contact of the Brunia with the Moldanubian and its northern root sector. The Moho beneath the Fore-Sudetic Monocline is the most pronounced and is interpreted as the first-order discontinuity at a depth of 30 km.     

Zircon U–Pb ages, Hf and O isotopes constrain the crustal architecture of the ultrahigh-pressure Dabie orogen in China

The crustal structure of the Dabie orogen was reconstructed by a combined study of U–Pb ages, Hf and O isotope compositions of zircons from granitic gneiss from North Dabie, the largest lithotectonic unit in the orogen. The results were deciphered from metamorphic history to protolith origin with respect to continental subduction and exhumation. Zircon U–Pb dating provides consistent ages of 751 ± 7 Ma for protolith crystallization, and two group ages of 213 ± 4 to 245 ± 17 Ma and 126 ± 4 to 131 ± 36 Ma for regional metamorphism. Majority of zircon Hf isotope analyses displays negative ε_Hf(t) values of − 5.1 to − 2.9 with crust Hf model ages of 1.84 to 1.99 Ga, indicating protolith origin from reworking of middle Paleoproterozoic crust. The remaining analyses exhibit positive ε_Hf(t) values of 5.3 to 14.5 with mantle Hf model ages of 0.74 to 1.11 Ga, suggesting prompt reworking of Late Mesoproterozoic to Early Neoproterozoic juvenile crust. Zircon O isotope analyses yield δ¹⁸O values of − 3.26 to 2.79‰, indicating differential involvement of meteoric water in protolith magma by remelting of hydrothermally altered low δ¹⁸O rocks. North Dabie shares the same age of Neoproterozoic low δ¹⁸O protolith with Central Dabie experiencing the Triassic UHP metamorphism, but it was significantly reworked at Early Cretaceous in association with contemporaneous magma emplacement. The Rodinia breakup at about 750 Ma would lead to not only the reworking of juvenile crust in an active rift zone for bimodal protolith of Central Dabie, but also reworking of ancient crust in an arc-continent collision zone for the North Dabie protolith. The spatial difference in the metamorphic age (Triassic vs. Cretaceous) between the northern and southern parts of North Dabie suggests intra-crustal detachment during the continental subduction. Furthermore, the Dabie orogen would have a three-layer structure prior to the Early Cretaceous magmatism: Central Dabie in the upper, North Dabie in the middle, and the source region of Cretaceous magmas in the lower.