Application of the numerical manifold method for stress wave propagation across rock masses

In this paper, the numerical manifold method (NMM) is extended to study wave propagation across rock masses. First, improvements to the system equations, contact treatment, and boundary conditions of the NMM are performed, where new system equations are derived based on the Newmark assumption of the space–time relationship, the edge‐to‐edge contact treatment is further developed for the NMM to handle stress wave propagation across discontinuities, and the viscous non‐reflection boundary condition is derived based on the energy minimisation principle. After the modification, numerical comparisons between the original and improved NMM are presented. The results show that the original system equations result in artificial numerical damping, which can be overcome by the Newmark system equations. Meanwhile, the original contact scheme suffers some calculation problems when modelling stress wave propagation across a discontinuity, which can be solved by the proposed edge‐to‐edge contact scheme. Subsequently, the influence of the mesh size and time step on the improved NMM for stress wave propagation is studied. Finally, 2D wave propagation is modelled, and the model's results are in good agreement with the analytical solution. Copyright © 2013 John Wiley & Sons, Ltd.     

On the relationship between historical land‐use change and water availability: the case of the lower Tarim River region in northwestern China

Natural ecosystems in the region of the lower Tarim River in northwestern China strongly deteriorated since the 1950s due to an expanding desertification. As a result, the downstream Tarim River reaches became permanently dry land. This historical evolution in land‐use change is typically the result of the anthropogenic impact on natural ecosystems. On the basis of a spatially distributed hydrological catchment model bidirectionally linked with a fully hydrodynamic MIKE11 river model, land‐use changes characterized by historical changes in leaf area index (LAI) of vegetation, as well as the evolution of irrigated surface areas, can be causally related to changes in water resources (groundwater storage and surface water resources). An increased surface area of irrigated (agricultural) land, together with a majority of inefficient irrigation methods, did lead to a strong increase of water resources consumption of the farmlands located in the upper Tarim River area. Evidently, this evolution influenced available water resources downstream in the Tarim basin. As a result, farmland has been gradually relocated to the upstream regions. This has led to reduced flows from the upper Tarim stream, which subsequently accelerated the dropping of the groundwater level downstream in the basin. This study moreover demonstrates that land surface biomass changes (cumulative LAI) along the lower Tarim River are strongly related to the changes in groundwater storage. Copyright © 2012 John Wiley & Sons, Ltd.     

Isotope geochemistry and Re–Os geochronology of the Yanjiagou Mo deposit in the central North China Craton

The Yanjiagou deposit, located in the central North China Craton (NCC), is a newly found porphyry‐type Mo deposit. The Mo mineralization here is spatially associated with the Mapeng batholith. In this study, we identify four stages of ore formation in this deposit: pyrite phyllic stage (I), quartz–pyrite stage (II), quartz–pyrite–molybdenite stage (III), which is the main mineralization stage, and quartz–carbonate stage (IV). We present sulphur and lead isotope data on pyrite, and rhenium and osmium isotopes of molybdenite from the porphyry deposit and evaluate the timing and origin of ore formation. The δ³⁴S values of the pyrite range from ‐1.1‰ to −0.6‰, with an average of −0.875‰, suggesting origin from a mixture of magmatic/mantle sources and the basement rocks. The Pb isotope compositions of the pyrite show a range of 16.369 to 17.079 for ²⁰⁶Pb/²⁰⁴Pb, 15.201 to 15.355 for ²⁰⁷Pb/²⁰⁴Pb, and 36.696 to 37.380 for ²⁰⁸Pb/²⁰⁴Pb, indicating that the ore‐forming materials were derived from a mixture of lower crust (or basement rocks) and mantle. Rhenium contents in molybdenite samples from the main ore stage are between 74.73 to 254.43 ppm, with an average of 147.9 ppm, indicating a mixed crustal‐mantle source for the metal. Eight molybdenite separates yield model ages ranging from 124.17 to 130.80 Ma and a mean model age of 128.46 Ma. An isochron age of 126.7 ± 1.1 Ma (MSWD = 2.1, initial ¹⁸⁷Os = 0.0032 ± 0.0012 ppb) is computed, which reveals a close link between the Mo mineralization and the magmatism that generated the Mapeng batholith. The age is close to the zircon U–Pb age of ca. 130 Ma from the batholith reported in a recent study. The age is also consistent with the timing of mineralization in the Fuping ore cluster in the central NCC, as well as the peak time of lithosphere thinning and destruction of the NCC. We evaluate the spatio‐temporal distribution of the Mo deposits in the NCC and identify three important molybdenum provinces along the northern and southern margins of the craton formed during three distinct episodes: Middle to Late Triassic (240–220 Ma), Early Jurassic (190–175 Ma), and Late Jurassic to Early Cretaceous (150–125 Ma). The third period is considered to mark the most important metallogenic event, coinciding with the peak of lithosphere thinning and craton destruction in the NCC. Copyright © 2014 John Wiley & Sons, Ltd.     

Iron isotopic systematics of UHP eclogites respond to oxidizing fluid during exhumation

The iron stable isotope compositions (δ⁵⁶Fe) and iron valence states of ultrahigh‐pressure eclogites from Bixiling in the Dabie orogen belt, China, were measured to trace the changes of geochemical conditions during vertical transportation of earth materials, for example, oxygen fugacity. The bulk Fe³⁺/ΣFe ratios of retrograde eclogites, determined by Mössbauer spectroscopy, are consistently higher than those of fresh eclogites, suggesting oxidation during retrograde metamorphism and fluid infiltration. The studied eclogites (five samples) display limited mid‐ocean ridge basalts (MORB)‐like (~0.10‰) δ⁵⁶Fe values, which are indistinguishable from their protoliths, that is, gabbro cumulates formed through differentiation of mantle‐derived basaltic magma. This suggests that Fe isotope fractionation during continental subduction is limited. Garnet separates display limited δ⁵⁶Fe variation ranging from ?0.08 ± 0.07 ‰ to 0.02 ± 0.07‰, whereas coexisting omphacite displays a large variation of δ⁵⁶Fe values from 0.15 ± 0.07‰ to 0.47 ± 0.07‰. Omphacite also has highly variable Fe³⁺/ΣFe ratios from 0.367 ± 0.025 to 0.598 ± 0.024, indicating modification after peak metamorphism. Omphacite from retrograde eclogites has elevated Fe³⁺/ΣFe ratios (0.54–0.60) compared to that from fresh eclogites (~0.37), whereas garnet displays a narrow range of ferric iron content with Fe³⁺/ΣFe ratios from 0.039 ± 0.013 to 0.065 ± 0.022. The homogenous δ⁵⁶Fe values and Fe³⁺/ΣFe ratios of garnet suggest that it survived the retrograde metamorphism and preserved its Fe‐isotopic features and ferric contents of peak metamorphism. Because of similar diffusion rates of Fe and Mg in garnet and omphacite, and constant Δ²⁶Mg_{omphacite‐garnet} values (1.14 ± 0.04‰), equilibrium iron isotope fractionation between garnet and omphacite was probably achieved during peak metamorphism. Elevated Fe³⁺/ΣFe ratios of omphacite from retrograde eclogites and variant Δ⁵⁶Fe_{omphacite‐garnet} values of the studied eclogites (0.13 ± 0.10‰ to 0.48 ± 0.10‰) indicate that oxidized geofluid infiltration resulted in the elevation of δ⁵⁶Fe values of omphacite during retrograde metamorphism.     

Effects of low- to high-angle normal faults on sedimentary architectures in the Eocene Wenchang Formation,Enping Sag,Pearl River Mouth Basin,South China Sea

The effects of low- to high-angle (>30°) normal faults on sedimentary architectural units in the Eocene Wenchang Formation, Enping Sag, Pearl River Mouth Basin (PRMB), South China Sea were investigated utilising a high-quality 3D seismic data set and restored paleogeomorphology. It has been shown that sequence stratigraphic units and sedimentary architecture are significantly controlled by the low- to high-angle normal faults. The Wenchang Formation, a second-order sequence, can be subdivided into two para-second-sequences (the Lower and Upper Wenchang sequences, E₂W^L and E₂W^U) and seven third-order sequences (from base to top: SQ1～SQ7). The low-angle fault confined sequence architecture of the Wenchang Formation is mainly characterised by lateral stacking with the ratio of the vertical subsidence (V) to horizontal slip (H) being reduced from 1/2 for E₂W^L to 1/6 for E₂W^U. In contrast, the high-angle fault confined sequence is characterised by vertical stacking with the ratio of V/H close to 1 for sequences SQ1 to SQ7. In the 3D seismic area, the features of sediment-dispersal pattern were interpreted based on an integrated analysis of paleogeomorphology, seismic reflection characteristics, stratal thickness distribution and multiple attribute clustering. The results show that the large-scale fan delta, belt-shape lacustrine deposit and bird-foot braided delta systems mainly developed in the low-angle fault confined sequences, whereas small-scale fan delta, rhombus-shaped lacustrine deposit and lobe-shaped braided delta systems inherited tectono-sedimentary architectures in the high-angle fault confined sequences.     

The Mesozoic structural evolution of the Gorgon Platform,North Carnarvon Basin,Australia

The Gorgon Platform is located on the southeastern edge of the Exmouth Plateau in the North Carnarvon Basin, North West Shelf, Australia. A structural analysis using three-dimensional (3D) seismic data has revealed four major sets of extensional faults, namely, (1) the Exmouth Plateau extensional fault system, (2) the basin bounding fault system (Exmouth Plateau–Gorgon Platform Boundary Fault), (3) an intra-rift fault system in the graben between the Exmouth Plateau and the Gorgon Platform and (4) an intra-rift fault system within the graben between the Exmouth Plateau and the Exmouth Sub-basin. Fault throw-length analyses imply that the initial fault segments, which formed the Exmouth Plateau–Gorgon Platform Boundary Fault (EG Boundary Fault), were subsequently connected vertically and laterally by both soft- and hard-linked structures. These major extensional fault systems were controlled by three different extensional events during the Early and Middle Jurassic, Late Jurassic and Early Cretaceous, and illustrate the strong role of structural inheritance in determining fault orientation and linkage. The Lower and Middle Jurassic and Upper Jurassic to Lower Cretaceous syn-kinematic sequences are separated by unconformities.