Phase relation of CaSO<Subscript>4</Subscript> at high pressure and temperature up to 90 GPa and 2300 K

Calcium sulfate (CaSO₄), one of the major sulfate minerals in the Earth’s crust, is expected to play a major role in sulfur recycling into the deep mantle. Here, we investigated the crystal structure and phase relation of CaSO₄ up to ~90 GPa and 2300 K through a series of high-pressure experiments combined with in situ X-ray diffraction. CaSO₄ forms three thermodynamically stable polymorphs: anhydrite (stable below 3 GPa), monazite-type phase (stable between 3 and ~13 GPa) and barite-type phase (stable up to at least 93 GPa). Anhydrite to monazite-type phase transition is induced by pressure even at room temperature, while monazite- to barite-type transition requires heating at least to 1500 K at ~20 GPa. The barite-type phase cannot always be quenched from high temperature and is distorted to metastable AgMnO₄-type structure or another modified barite structure depending on pressure. We obtained the pressure–volume data and density of anhydrite, monazite- and barite-type phases and found that their densities are lower than those calculated from the PREM model in the studied P–T conditions. This suggests that CaSO₄ is gravitationally unstable in the mantle and fluid/melt phase into which sulfur dissolves and/or sulfate–sulfide speciation may play a major role in the sulfur recycling into the deep Earth.     

Geochronological and isotopic records of crustal storage and assimilation in the Wolverine Creek–Conant Creek system,Heise eruptive centre,Snake River Plain

Understanding the processes of differentiation of the Yellowstone–Snake River Plain (YSRP) rhyolites is typically impeded by the apparent lack of erupted intermediate compositions as well as the complex nature of their shallow interaction with the surrounding crust responsible for their typically low O isotopic ratios. A pair of normal-δ¹⁸O rhyolitic eruptions from the Heise eruptive centre in eastern Idaho, the Wolverine Creek Tuff and the Conant Creek Tuff, represent unique magmatic products of the Yellowstone hotspot preserving abundant vestiges of the intermediate differentiation steps leading to rhyolite generation. We address both shallow and deep processes of magma generation and storage in the two units by combining high-precision ID–TIMS U–Pb zircon geochronology, trace element, O and Hf isotopic studies of zircon, and Sr isotopic analyses of individual high-Mg# pyroxenes inherited from lower- to mid-crustal differentiation stages. The zircon geochronology confirms the derivation of both tuffs from the same rhyolitic magma reservoir erupted at 5.5941 ± 0.0097 Ma, preceded by at least 92 ± 14 ky of continuous or intermittent zircon saturation approximating the length of pre-eruptive magma accumulation in the upper crust. Some low-Mg# pyroxenes enclosing zircons predate the eruption by at least 45 ± 27 ky, illustrating the co-crystallisation of major and accessory phases in the near-liquidus rhyolitic melts of the YSRP over a significant period of time. Coeval zircon crystals are isotopically heterogeneous (two populations at εHf ~?5 and ?13), requiring the assembly of isotopically distinct melt pockets directly prior to, or during, the eruption. The primitive Mg# 60–90 pyroxenes are out of isotopic equilibrium with the host rhyolitic melt (⁸⁷Sr/⁸⁶Sr_i = 0.70889), covering a range of ⁸⁷Sr/⁸⁶Sr_i = 0.70705–0.70883 corresponding to ratios typical of the most radiogenic YSRP basalts to the least radiogenic YSRP rhyolites. Together with the low εHf in zircon, the Sr isotopic ratios illustrate limited assimilation dominated by radiogenic Archean crustal source materials incorporated into variably evolved YSRP melts as they progress towards rhyolitic compositions by assimilation–fractional crystallisation.     

Characteristics of deuterium excess parameters for geothermal water in Beijing

This study investigates the characteristics of geothermal water in 10 geothermal fields in Beijing. The relationships between the deuterium excess parameter (d) and temperature, depth, age of geothermal groundwater, groundwater flow field, and Eh were investigated using geothermal groundwater samples. Results showed that (1) the average d value of geothermal water is 5.4, whereas that of the groundwater in normal temperature is 6.04. The differences are induced by the oxygen isotope exchange during the water–rock interaction, which may be more easily completed in geothermal water than in cold groundwater. (2) The d value increases remarkably with the age of the geothermal groundwater. The d value increases from 11.2 to 14.6 when the age of the geothermal water is 12,760 ± 130 a and 38,960 ± 630 a, respectively. Moreover, the isotope heat exchange for composition of the hydrogen and oxygen isotopes in the geothermal groundwater proceeds sufficiently with time. (3) The d value decreases from 5.72 to 3.03 when the depth increases from 125.13 to 3221 m. Generally, in the same area, the d value decreases with depth because the temperature is increasing. (4) The d value of the groundwater gradually reduces from the northern recharge area to the southern discharge area. The average d value is 7.31 in the northern recharge area and 5.68 in the middle Beijing Depression, whereas the d value in the southern area of Fengheying is ?9.20. The larger difference in d values between the recharge and discharge areas is due to the slower velocity of underwater flow, which induces longer time for oxygen exchange. (5) The relationship between the d and Eh is complex. When Eh is <200 mV, the d value of the geothermal water decreases with the decrease in Eh. When Eh is higher than 200 mV, the d value increases slightly with the decrease in Eh. The study of the characteristics of deuterium excess parameters for geothermal water could provide a scientific isotopic evidence for assessment and exploitation measures in geothermal groundwater systems.     

2D probabilistic prediction of sparsely measured earth properties constrained by geophysical imaging fully accounting for tomographic reconstruction ambiguity

Many hydrological, environmental, or engineering exploration tasks require predicting spatially continuous scenarios of sparsely measured borehole logging data. We present a methodology to probabilistically predict such scenarios constrained by ill-posed geophysical tomography. Our approach allows for transducing tomographic reconstruction ambiguity into the probabilistic prediction of spatially continuous target parameter scenarios. It is even applicable to data sets where petrophysical relations in the survey area are non-unique, i.e., different facies related petrophysical relations may be present. We employ static two-layer artificial neural networks (ANNs) for prediction and additionally evaluate, whether the training performance of the ANNs can be used to rank geophysical tomograms, which are mathematically equal reconstructions of physical parameter distributions in the ground. We illustrate our methodology using a realistic synthetic database for maximal control about the prediction performance and ranking potential of the approach. For doing so, we try to link geophysical radar and seismic tomography as input parameters to porosity of the ground as target parameter of ANN. However, the approach is flexible and can cope with any combination of geophysical tomograms and hydrologic, environmental or engineering target parameters. Ranking of equivalent geophysical tomograms based on additional borehole logging data is found to be generally possible, but risks remain that the ranking based on the ANN training performance does not fully coincide with the closeness of geophysical tomograms to ground truth. Since geophysical field data sets do usually not offer control options similar to those used in our synthetic database, we do not recommend the utilization of recurrent ANNs to learn weights for the individual geophysical tomograms used in the prediction procedure.     

Using cloud water path and cloud top temperature for estimating convective and stratiform rainfall from SEVIRI daytime data

A satellite rainfall retrieval technique is proposed here. The relationships of rain rate with each of cloud water path (CWP) and cloud top temperature (CTT) are investigated. The CWP and CTT are retrieved from SEVIRI data (spinning enhanced visible and infrared imager), and corresponding rain rates are measured by weather radar. The rain rates are compared to corresponding CWP and then to corresponding CTT. The investigation demonstrates an exponential functional dependency between rain rates and CWP for low and moderate rain rates (stratiform rainfall). Conversely, the rain rates are more closely related to CTT for high rain rates (convective rainfall). Therefore, two separate relationships are established for rain rate retrievals. The results show rain rates estimated by the developed scheme are in good correlation with those observed by weather radar.     

Mineralogy and geochemistry of granitoids from Kinnaur region,Himachal Higher Himalaya,India: Implication on the nature of felsic magmatism in the collision tectonics

Felsic magmatism in the southern part of Himachal Higher Himalaya is constituted by Neoproterozoic granite gneiss (GGn), Early Palaeozoic granitoids (EPG) and Tertiary tourmaline-bearing leucogranite (TLg). Magnetic susceptibility values (<3 ×10^?3 SI), molar Al₂ O ₃/(CaO + Na₂ O + K ₂O) (≥1.1), mineral assemblage (bt–ms–pl–kf–qtz ± tur ± ap), and the presence of normative corundum relate these granitoids to peraluminous S-type, ilmenite series (reduced type) granites formed in a syncollisional tectonic setting. Plagioclase from GGn (An₁₀–An₃₁) and EPG (An₁₅–An₃₃) represents oligoclase to andesine and TLg (An₂–An₁₅) represents albite to oligoclase, whereas compositional ranges of K-feldspar are more-or-less similar (Or₈₈ to Or₉₅ in GGn, Or₈₆ to Or₉₇ in EPG and Or₈₇ to Or₉₄ in TLg). Biotites in GGn (Mg/Mg + Fe^t= 0.34–0.45), EPG (Mg/Mg + Fe^t= 0.27–0.47), and TLg (Mg/Mg + Fe^t= 0.25–0.30) are ferribiotites enriched in siderophyllite, which stabilised between FMQ and HM buffers and are characterised by dominant 3Fe\(\rightleftharpoons \)2Al, 3Mg\(\rightleftharpoons \)2Al substitutions typical of peraluminous (S-type), reducing felsic melts. Muscovite in GGn (Mg/Mg + Fe^t=0.58–0.66), EPG (Mg/Mg + Fe^t=0.31?0.59), and TLg (Mg/Mg + Fe^t=0.29–0.42) represent celadonite and paragonite solid solutions, and the tourmaline from EPG and TLg belongs to the schorl-elbaite series, which are characteristics of peraluminous, Li-poor, biotite-tourmaline granites. Geochemical features reveal that the GGn and EPG precursor melts were most likely derived from melting of biotite-rich metapelite and metagraywacke sources, whereas TLg melt appears to have formed from biotite-muscovite rich metapelite and metagraywacke sources. Major and trace elements modelling suggest that the GGn, EPG and TLg parental melts have experienced low degrees (～13, ～17 and ～13%, respectively) of kf–pl–bt fractionation, respectively, subsequent to partial melting. The GGn and EPG melts are the results of a pre-Himalayan, syn-collisional Pan-African felsic magmatic event, whereas the TLg is a magmatic product of Himalayan collision tectonics.     

Corrosion behavior of anodic oxidized TiO<Subscript>2</Subscript> film in seawater

TiO₂ films were formed on metallic titanium substrates by the anodic oxidation method in H₂SO₄ solution under the 80V D.C.. Phase component and microstructure were characterized by X-ray diffraction (XRD) and scanning electron microscopy (SEM). Water contact angles on titanium oxide film surface were measured under both dark and sunlight illumination conditions. Corrosion tests were carried out in seawater under different illumination conditions by electrochemistry impedance spectrum (EIS) and polarization curves. The result showed that the TiO₂ film prepared by the anodic oxidation method was anatase with a uniform structure and without obvious pores or cracks on its surface. The average water contact angle of the film was 116.4° in dark, in contrast to an angle of 42.7° under the UV illumination for 2 hours, which demonstrates good hydrophobic property. The anti-corrosion behavior of the TiO₂ film was declining with the extended immersion time. Under dark conditions, however, the hydrophobic TiO₂ film retarded the water infiltrating into the substrate. The impedance changed slowly and the corrosion current density was 2 orders of magnitude lower than that with the film illuminated by sunlight. All of those mentioned above indicate that the TiO₂ film possesses much better performance under dark condition, and it can be applied as an engineering material under dark seawater environment.     

Fine-resolution simulation of surface current and sea ice in the Arctic Mediterranean Seas

A fine-resolution model is developed for ocean circulation simulation in the National Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG), Chinese Academy of Sciences, and is applied to simulate surface current and sea ice variations in the Arctic Mediterranean Seas. A dynamic sea ice model in elastic-viscous-plastic rheology and a thermodynamic sea ice model are employed. A 200-year simulation is performed and a dimatological average of a 10-year period (141st-150th) is presented with focus on sea ice concentration and surface current variations in the Arctic Mediterranean Seas. The model is able to simulate well the East Greenland Current, Beaufort Gyre and the Transpolar Drift, but the simulated West Spitsbergen Current is small and weak. In the March climatology, the sea ice coverage can be simulated well except for a bit more ice in east of Spitsbergen Island. The result is also good for the September scenario except for less ice concentration east of Greenland and greater ice concentration near the ice margin. The extra ice east of Spitsbergen Island is caused by sea ice current convergence forced by atmospheric wind stress.     

Assessment on agricultural drought risk based on variable fuzzy sets model

Drought is one of the major natural disasters causing huge agricultural losses annually. Regional agricultural drought risk assessment has great significance for reducing regional disaster and agricultural drought losses. Based on the fuzzy characteristics of agricultural drought risk, variable fuzzy sets model was used for comprehensively assessing agricultural drought risk of Liaoning Province in China. A multi-layers and multi-indices assessment model was established according to variable fuzzy sets theory, and agricultural drought risk of all 14 prefecture-level cities was respectively estimated in terms of dangerousness, vulnerability, exposure and drought-resistibility. By calculating the combination weights of four drought risk factors, agricultural drought risk grade of each city was obtained. Based on the assessment results, the spatial distribution maps of agricultural drought risk were drawn. The results shows that eastern cities have lower drought dangerousness than western cities in Liaoning Province totally. Most cities are located in low drought vulnerability region and high drought exposure region. Because of frequent and severe drought since 2000, most cities are located in lower drought-resistibility region. Comprehensive agricultural drought risk presents apparent spatial characteristics, escalating from the east to the west. Drought dangerousness is the most important factor influencing comprehensive agricultural drought risk. Through the spatial distribution maps of drought risk, decision makers could find out drought situation and make decisions on drought resistance conveniently.     

Fault type analysis along the San Andreas Fault zone: A numerical approach

Finite Element (FE) modeling under plane stress condition is used to analyze the fault type variation with depth along and around the San Andreas Fault (SAF) zone. In this simulation elastic rheology was used and was thought justifiable as the variation in depth from 0.5 km to 20 km was considered. Series of calculations were performed with the variation in domain properties. Three types of models were created based on simple geological map of California, namely, 1) single domain model considering whole California as one homogeneous domain, 2) three domains model including the North American plate, Pacific plate, and SAF zone as separate domains, and 3) Four domains model including the three above plus the Garlock Fault zone. Mohr-Coulomb failure criterion and Byerlee's law were used for the calculation of failure state. All the models were driven by displacement boundary condition imposing the fixed North American plate and Pacific plate motion along N34°W vector up to the northern terminus of SAF and N50°E vector motion for the subducting the Gorda and Juan de Fuca plates. Our simulated results revealed that as the depth increased, the fault types were generally normal, and at shallow depth greater strike slip and some thrust faults were formed. It is concluded that SAF may be terminated as normal fault at depth although the surface expression is clearly strike slip.