The role of heterogeneous strain in the development and preservation of a polymetamorphic record in high-P granulites, western Canadian Shield

Mafic rocks in the Chipman domain of the Athabasca granulite terrane, western Canadian Shield, provide the first well‐documented record of two distinct high‐P granulite facies events in the same domain in this region. Textural relations and the results of petrological modelling (NCFMASHT system) of mafic granulites are interpreted in terms of a three‐stage tectonometamorphic history. Stage 1 involved development of the assemblage Grt + Cpx + Qtz ± Pl (M₁) from a primary Opx‐bearing igneous precursor at conditions of 1.3 GPa, 850–900 °C. Field and microstructural observations suggest that M₁ developed synchronously with an early S₁ gneissic fabric. Stage 2 is characterized by heterogeneous deformation (D₂) and synkinematic partial retrogression of the peak assemblage to an amphibole‐bearing assemblage (M₂). Stage 3 involved a third phase of deformation and a return to granulite facies conditions marked by the prograde breakdown of amphibole (Amph₂) to produce matrix garnet (Grt_3a) and the coronitic assemblage Cpx_3b + Opx_3b + Ilm_3b + Pl_3b (M_3b) at 1.0 GPa, 800–900 °C. M₁ and M_3b are correlated with 2.55 and 1.9 Ga metamorphic generations of zircon, respectively, which were dated in a separate study. Heterogeneous strain played a crucial role in both the development and preservation of these rare examples of multiple granulite facies events within single samples. Without this fortuitous set of circumstances, the apparent reaction history could have incorrectly led to an interpretation involving a single‐cycle high‐grade event. The detailed P–T–t–D history constructed for these rocks provides the best evidence to date that much of the east Lake Athabasca region experienced long‐term lower crustal residence from 2.55 to 1.9 Ga, and thus the region represents a rare window into the reactivation and ultimate stabilization processes of cratonic lithosphere.     

First U–Pb SHRIMP age of the Hauterivian stage, Neuquén Basin, Argentina

A high-resolution ion-microprobe (SHRIMP) U–Pb zircon age from a tuff layer intercalated in the ammonoid bearing sedimentary succession of the Neuquén Basin in Argentina provides a robust geochronologic date to add to the absolute ages and to improve the relative chronology of the Early Cretaceous Hauterivian stage. The tuff layer appears interbedded between shales of the upper member (Agua de la Mula) of the Agrio Formation within the Spitidiscus riccardii ammonoid zone (base of the Late Hauterivian) yielding a date of 132.5 ± 1.3 Ma. This date confirms and supports an accurate correlation between the ammonoid biostratigraphy of the Neuquén Basin with the Western Mediterranean Province of the Tethys during the Early Cretaceous and matches with the most recently published time scale. It also casts doubts on the validity of K–Ar ages on glauconite-grains recently reported from the Lower Cretaceous of the Vocontian Basin of France.     

含煤层地质环境下地震波场的数值模拟

本文对含低速煤层地质环境下弹性波场多波多分量地震资料进行了二维数值模拟研究,对人工边界反射进行了有效处理,频散效应得到了有效的压制,对几种不同激发与观测排列方式下的弹性波资料进行了模型计算与分析。     

Continental crust and lithospheric mantle interaction beneath North China: isotopic evidence from granulite xenoliths in Hannuoba, Sino-Korean craton

Mafic granulite and pyroxenite xenoliths from Cenozoic alkaline basalts at Hannuoba, Hebei Province, North China have been selected for a systematic geochemical and Sr–Nd–Pb isotopic study, which provides a unique opportunity to explore nature of the lower crust and the interaction between the continental crust and lithospheric mantle beneath an Archean craton. The major, compatible and incompatible elements and radiogenic isotopes of these xenoliths suggest great chemical heterogeneity of the lower crust beneath the Hannuoba region. Petrological and geochemical evidences indicate a clear cumulate origin, and most likely, they are related to basaltic underplating in different geological episodes. However, the Sr–Nd–Pb isotopic compositions of the xenoliths reveal a profound enriched source signature (EM I) with some influence of EM II, which implies that some portion of pre-existing, old metasomatized subcontinental lithospheric mantle could have played an important role in their genesis. It is suggested that the interaction between continental crust and subcontinental mantle as manifested by basaltic underplating would be closely related to regional tectonic episodes and geodynamic processes in the deep part of subcontinental lithospheric mantle.     

有限弹性变形地壳中的地震波

地壳由半无限大的基岩上一层厚度为H^-的表土层组成，入射地震波为垂直的SH波，产生水平地面运动。当浅源大地震发生时，在极震区以外行波传播产生地面运动将使地壳介质有非线性的有限弹性变形。用小参数摄动法使非线性控制方程为线性化的小参数各阶控制方程，得出头两阶线性控制方程的解析解。     

Lithium and lithium isotope profiles through the upper oceanic crust: a study of seawater-basalt exchange at ODP Sites 504B and 896A

Ocean Drilling Program (ODP) Hole 504B near the Costa Rica Rift is the deepest hole drilled in the ocean crust, penetrating a volcanic section, a transition zone and a sheeted dike complex. The distribution of Li and its isotopes through this 1.8-km section of oceanic crust reflects the varying conditions of seawater alteration with depth. The upper volcanic rocks, altered at low temperatures, are enriched in Li (5.6-27.3 ppm) and have heavier isotopic compositions (δ⁷Li=6.6-20.8‰) relative to fresh mid-ocean ridge basalt (MORB) due to uptake of seawater Li into alteration clays. The Li content and isotopic compositions of the deeper volcanic rocks are similar to MORB, reflecting restricted seawater circulation in this section. The transition zone is a region of mixing of seawater with upwelling hydrothermal fluids and sulfide mineralization. Li enrichment in this zone is accompanied by relatively light isotopic compositions (−0.8-2.1‰) which signify influence of basalt-derived Li during mineralization and alteration. Li decreases with depth to 0.6 ppm in the sheeted dike complex as a result of increasing hydrothermal extraction in the high-temperature reaction zone. Rocks in the dike complex have variable isotopic values that range from −1.7 to 7.9‰, depending on the extent of hydrothermal recrystallization and off-axis low-temperature alteration. Hydrothermally altered rocks are isotopically light because ⁶Li is preferentially retained in greenschist and amphibolite facies minerals. The δ⁷Li values of the highly altered rocks of the dike complex are complementary to those of high-temperature mid-ocean ridge vent fluids and compatible to equilibrium control by the alteration mineral assemblage. The inventory of Li in basement rocks permits a reevaluation of the role of oceanic crust in the budget of Li in the ocean. On balance, the upper 1.8 km of oceanic crusts remains a sink for oceanic Li. The observations at 504B and an estimated flux from the underlying 0.5 km of gabbro suggest that the global hydrothermal flux is at most 8×10⁹ mol/yr, compatible with geophysical thermal models. This work defines the distribution of Li and its isotopes in the upper ocean crust and provides a basis to interpret the contribution of subducted lithosphere to arc magmas and cycling of crustal material in the deep mantle.     

Wide-angle seismic velocities in heterogeneous crust

Seismic velocities measured by wide-angle surveys are commonly used to constrain material composition in the deep crust. Therefore, it is important to understand how these velocities are affected by the presence of multiscale heterogeneities. The effects may be characterised by the scale of the heterogeneity relative to the dominant seismic wavelength (λ); what is clear is that heterogeneities of all scales and strengths bias wide-angle velocities to some degree. Waveform modelling was used to investigate the apparent wide-angle P -wave velocities of different heterogeneous lower crusts. A constant composition (50 per cent felsic and 50 per cent ultramafic) was formed into a variety of 1- and 2-D heterogeneous arrangements and the resulting wide-angle seismic velocity was estimated. Elastic, 1-D models produced the largest velocity shift relative to the true average velocity of the medium (which is the velocity of an isotropic mixture of the two components). Thick (width > λ) horizontal layers, as a result of Fermat's Principle, provided the largest increase in velocity; thin (width ≪λ) vertical layers produced the largest decrease in velocity. Acoustic 2-D algorithms were shown to be inadequate for modelling the kinematics of waves in bodies with multiscale heterogeneities. Elastic, 2-D modelling found velocity shifts (both positive and negative) that were of a smaller magnitude than those produced by 1-D models. The key to the magnitude of the velocity shift appears to be the connectivity of the fast (and/or slow) components. Thus, the models with the highest apparent levels of connectivity between the fast phases, the 1-D layers, produced the highest-magnitude velocity shifts. To understand the relationship between measured seismic velocities and petrology in the deep crust it is clear that high-resolution structural information (which describes such connectivity) must be included in any modelling.     

构造层次与大陆壳动力学机制转变关系

针对构造层次研究现状和存在的问题，把大陆壳划分为深部，中部，浅部和浅－表部四个构造层次。依据各自所处特定的构造位置、组成构造岩类型、形成的制约因素和地质时代等方面的区别，各构造层次分别是壳－幔间滑动，大陆张裂、隆－滑构造和变质核杂岩构造等多种大陆壳动力学机制转变过程中的产物。构造层次与大陆壳动力学机制转变关系的确定，更有利于古老板构造连续性、整体性的研究以及多期、多层次、多旋回大陆地壳演化模式的建     

Attenuation differences in layer 2A in intermediate- and slow-spreading oceanic crust

In situ seismic attenuationQ⁻¹logs are derived from borehole velocity profiles and reveal sharp boundaries between morphologies of the extrusive volcanic layers in intermediate- and slow-spreading oceanic crust.Q⁻¹logs are calculated from the scattering attenuation associated with vertical velocity heterogeneity in Ocean Drilling Program Holes 504B and 896A and in Hole 395A, located in 5.9–7.3 Ma crust on the Pacific and Atlantic plates, respectively. Our results strongly tie crustal properties to seismic measurables and observed geological structures: we find that the scattering attenuation can be used to identify the extrusive volcanic sequence because it is closely related to changes in the degree of vertical heterogeneity. We interpret a distinct decrease in the Q⁻¹log at the transition below the extrusive volcanic layer to correspond with the seismic layer 2A/2B boundary. The boundary is located at 465 m depth below the sea floor in both Hole 395A and 504B, although this is likely to be a coincidence of the sediment thickness at these sites. Layer 2A is estimated to be approximately 150 m thick in Hole 504B and > 300 m thick in Hole 395A. Cyclic sequences of high-porosity pillows and low-porosity massive units in the uppermost 100 m of volcanics in Hole 395A result in large velocity heterogeneities which cause > 5 times more attenuation in this layer than in Hole 504B. In Hole 896A, by contrast, fewer pillows, more massive flows, and a greater volume of carbonate veins decrease the velocity heterogeneity and attenuation significantly over only 1 km distance from Hole 504B. We conclude that the attenuation in the extrusive volcanics of the ocean crust is largely controlled by variation in local heterogeneity and morphology as well as by subsequent hydrothermal alteration. The observed differences inQ⁻¹profiles and layer 2A thickness at these sites may be attributed to variations in the volume and duration of volcanic activity at mid-ocean spreading centers for these Pacific and Atlantic ridge segments.