Mineral and fluid inclusions in zircon of UHP metamorphic rocks from the CCSD-main drill hole: A record of metamorphism and fluid activity

The Chinese Continental Scientific Drilling (CCSD) main drill hole (0–3000 m) in Donghai, southern Sulu orogen, consists of eclogite, paragneiss, orthogneiss, schist and garnet peridotite. Detailed investigations of Raman, cathodoluminescence, and microprobe analyses show that zircons from most eclogites, gneisses and schists have oscillatory zoned magmatic cores with low-pressure mineral inclusions of Qtz, Pl, Kf and Ap, and a metamorphic rim with relatively uniform luminescence and eclogite-facies mineral inclusions of Grt, Omp, Phn, Coe and Rt. The chemical compositions of the UHP metamorphic mineral inclusions in zircon are similar to those from the matrix of the host rocks. Similar UHP metamorphic P–T conditions of about 770 °C and 32 kbar were estimated from coexisting minerals in zircon and in the matrix. These observations suggest that all investigated lithologies experienced a joint in situ UHP metamorphism during continental deep subduction. In rare cases, magmatic cores of zircon contain coesite and omphacite inclusions and show patchy and irregular luminescence, implying that the cores have been largely altered possibly by fluid–mineral interaction during UHP metamorphism.

Abundant H₂O–CO₂, H₂O- or CO₂-dominated fluid inclusions with low to medium salinities occur isolated or clustered in the magmatic cores of some zircons, coexisting with low-P mineral inclusions. These fluid inclusions should have been trapped during magmatic crystallization and thus as primary. Only few H₂O- and/or CO₂-dominated fluid inclusions were found to occur together with UHP mineral inclusions in zircons of metamorphic origin, indicating that UHP metamorphism occurred under relatively dry conditions. The diversity in fluid inclusion populations in UHP rocks from different depths suggests a closed fluid system, without large-scale fluid migration during subduction and exhumation.     

Site scale modeling of slow-moving landslides,a 3D viscoplastic finite element modeling approach

This paper presents an advanced 3D numerical methodology to reproduce the kinematics of slow active landslides, more precisely, to reproduce the nearly constant strain rate (secondary creep) and the acceleration/deceleration of the moving mass due to hydrological changes. For this purpose, finite element analyses are performed in a large area covering a long time-span (12 years), in order to exhibit different interacting slope movements. First, we perform a stability analysis using the shear strength reduction (SSR) technique with a Mohr-Coulomb failure criteria. It is done in order to compute factors of safety (FS) and to identify two different scenarios, the first one being stable (FS > 1) and the second one being unstable (FS < 1). In the studied test case, the Portalet landslide (Central Spanish Pyrenees), the first scenario corresponds to an initial stable configuration of the slope and the second one to an unstable excavated configuration. Second, taking the first scenario as an initial condition, a time-dependent analysis is performed using a coupled formulation to model solid skeleton and pore fluids interaction, and a simplified ground water model that takes into account daily rainfall intensity. In this case, a viscoplastic constitutive model based on Perzyna’s theory is applied to reproduce soil viscous behavior and the delayed creep deformation due to the excavation. The fluidity parameter is calibrated to reproduce displacements measured by the monitoring systems. Our results demonstrate that 3D analyses are preferable to 2D ones for reproducing in a more realistic way the slide behavior. After calibration, the proposed model is able to simulate successfully short- and medium-term predictions during stages of primary and secondary creep.     

Sealing efficiency analysis for shallow-layer caprocks in CO<Subscript>2</Subscript> geological storage

The CO₂ migrated from deeper to shallower layers may change its phase state from supercritical state to gaseous state (called phase transition). This phase transition makes both viscosity and density of CO₂ experience a sharp variation, which may induce the CO₂ further penetration into shallow layers. This is a critical and dangerous situation for the security of CO₂ geological storage. However, the assessment of caprock sealing efficiency with a fully coupled multi-physical model is still missing on this phase transition effect. This study extends our previous fully coupled multi-physical model to include this phase transition effect. The dramatic changes of CO₂ viscosity and density are incorporated into the model. The impacts of temperature and pressure on caprock sealing efficiency (expressed by CO₂ penetration depth) are then numerically investigated for a caprock layer at the depth of 800 m. The changes of CO₂ physical properties with gas partial pressure and formation temperature in the phase transition zone are explored. It is observed that phase transition revises the linear relationship of CO₂ penetration depth and time square root as well as penetration depth. The real physical properties of CO₂ in the phase transition zone are critical to the safety of CO₂ sequestration. Pressure and temperature have different impact mechanisms on the security of CO₂ geological storage.     

Cumulative strain and boundary effects in the gravity flow of a viscous slab

A simple mathematical model based on a well-known flow equation permits determination of progressive strains at all points within a viscous slab flowing in an inclined rectangular channel. Differences in strain patterns between central and lateral domains of the slab may aid in determining boundary effects within natural structures formed by gravity-sliding.     

A model of the quiet solar atmosphere

The solar atmosphere may be divided into a number of isolated active components and a quiet residue. On the largest scale the latter is dominated by a general dipole magnetic field of strength 1–2 G; its observable components are flux concentrations in supergranule boundary regions (SBRs), spicules, mottles and polar plumes. The velocity field in the SBRs is discussed. There are continuous gas streaming motions up and down between the photosphere and the corona; spicules may be mainly downward moving gas.A unifying model is developed of these various components, as well as the heating mechanism of the whole quiet atmosphere. Highly ordered velocity fields of the cell, together with a gravitational wave, cause a vertical magnetic force tube to collapse below a critical level; the result is an upward eruption of a vortex ring at the Alfvén velocity. The complex mass velocity pattern may explain spicules, mottles and plumes, as well as unobservable streaming motions.The quiet atmosphere is divided into regions above SBRs and those above the inner parts of the cells. Hydromagnetic eruptions from the former may account for the entire heat requirement of the atmosphere. The model atmosphere has a chromosphere-corona transition layer which bulges upwards above the SBRs and so conforms with EUV data. The energy and mass balances in this solar atmosphere are considered, and it is also shown to be consistent with the radio data.     

Field Investigation of a New Recharge Approach for ASR Projects in Near‐Surface Aquifers

Aquifer storage and recovery (ASR) is the artificial recharge and temporary storage of water in an aquifer when water is abundant, and recovery of all or a portion of that water when it is needed. One key limiting factor that still hinders the effectiveness of ASR is the high costs of constructing, maintaining, and operating the artificial recharge systems. Here we investigate a new recharge method for ASR in near‐surface unconsolidated aquifers that uses small‐diameter, low‐cost wells installed with direct‐push (DP) technology. The effectiveness of a DP well for ASR recharge is compared with that of a surface infiltration basin at a field site in north‐central Kansas. The performance of the surface basin was poor at the site due to the presence of a shallow continuous clay layer, identified with DP profiling methods, that constrained the downward movement of infiltrated water and significantly reduced the basin recharge capacity. The DP well penetrated through this clay layer and was able to recharge water by gravity alone at a much higher rate. Most importantly, the costs of the DP well, including both the construction and land costs, were only a small fraction of those for the infiltration basin. This low‐cost approach could significantly expand the applicability of ASR as a water resources management tool to entities with limited fiscal resources, such as many small municipalities and rural communities. The results of this investigation demonstrate the great potential of DP wells as a new recharge option for ASR projects in near‐surface unconsolidated aquifers.     

The suspended sediment concentration distribution in the Bohai Sea, Yellow Sea and East China Sea

The distribution of the suspended sediment concentration (SSC) in the Bohai Sea, Yellow Sea and East China Sea (BYECS) is studied based on the observed turbidity data and model simulation results. The observed turbidity results show that (i) the highest SSC is found in the coastal areas while in the outer shelf sea areas turbid water is much more difficult to observe, (ii) the surface layer SSC is much lower than the bottom layer SSC and (iii) the winter SSC is higher than the summer SSC. The Regional Ocean Modeling System (ROMS) is used to simulate the SSC distribution in the BYECS. A comparison between the modeled SSC and the observed SSC in the BYECS shows that the modeled SSC can reproduce the principal features of the SSC distribution in the BYECS. The dynamic mechanisms of the sediment erosion and transport processes are studied based on the modeled results. The horizontal distribution of the SSC in the BYECS is mainly determined by the current-wave induced bottom stress and the fine-grain sediment distribution. The current-induced bottom stress is much higher than the wave-induced bottom stress, which means the tidal currents play a more significant role in the sediment resuspension than the wind waves. The vertical mixing strength is studied based on the mixed layer depth and the turbulent kinetic energy distribution in the BYECS. The strong winter time vertical mixing, which is mainly caused by the strong wind stress and surface cooling, leads to high surface layer SSC in winter. High surface layer SSC in summer is restricted in the coastal areas.