Impact of 3D Var GSI-ENKF hybrid data assimilation system

The hybrid two-way coupled 3DEnsVar assimilation system was tested with the NCMRWF global data assimilation forecasting system. At present, this system consists of T574L64 deterministic model and the grid-point statistical interpolation analysis scheme. In this experiment, the analysis system is modified with a two-way coupling with an 80 member Ensemble Kalman Filter of T254L64 resolution and runs are carried out in parallel to the operational system for the Indian summer monsoon season (June–September) for the year 2015 to study its impact. Both the assimilation systems are based on NCEP GFS system. It is found that hybrid assimilation marginally improved the quality of the forecasts of all variables over the deterministic 3D Var system, in terms of statistical skill scores and also in terms of circulation features. The impact of the hybrid system in prediction of extreme rainfall and cyclone track is discussed.     

Role of land state in a high resolution mesoscale model for simulating the Uttarakhand heavy rainfall event over India

In 2013, Indian summer monsoon witnessed a very heavy rainfall event (>30 cm/day) over Uttarakhand in north India, claiming more than 5000 lives and property damage worth approximately 40 billion USD. This event was associated with the interaction of two synoptic systems, i.e., intensified subtropical westerly trough over north India and north-westward moving monsoon depression formed over the Bay of Bengal. The event had occurred over highly variable terrain and land surface characteristics. Although global models predicted the large scale event, they failed to predict realistic location, timing, amount, intensity and distribution of rainfall over the region. The goal of this study is to assess the impact of land state conditions in simulating this severe event using a high resolution mesoscale model. The land conditions such as multi-layer soil moisture and soil temperature fields were generated from High Resolution Land Data Assimilation (HRLDAS) modelling system. Two experiments were conducted namely, (1) CNTL (Control, without land data assimilation) and (2) LDAS, with land data assimilation (i.e., with HRLDAS-based soil moisture and temperature fields) using Weather Research and Forecasting (WRF) modelling system. Initial soil moisture correlation and root mean square error for LDAS is 0.73 and 0.05, whereas for CNTL it is 0.63 and 0.053 respectively, with a stronger heat low in LDAS. The differences in wind and moisture transport in LDAS favoured increased moisture transport from Arabian Sea through a convectively unstable region embedded within two low pressure centers over Arabian Sea and Bay of Bengal. The improvement in rainfall is significantly correlated to the persistent generation of potential vorticity (PV) in LDAS. Further, PV tendency analysis confirmed that the increased generation of PV is due to the enhanced horizontal PV advection component rather than the diabatic heating terms due to modified flow fields. These results suggest that, two different synoptic systems merged by the strong interaction of moving PV columns resulted in the strengthening and further amplification of the system over the region in LDAS. This study highlights the importance of better representation of the land surface fields for improved prediction of localized anomalous weather event over India.     

Assessing the vegetation canopy influences on wind flow using wind tunnel experiments with artificial plants

Wind erosion causes serious problems and considerable threat in most regions of the world. Vegetation on the ground has an important role in controlling wind erosion by covering soil surface and absorbing wind momentum. A set of wind tunnel experiments was performed to quantitatively examine the effect of canopy structure on wind movement. Artificial plastic vegetations with different porosity and canopy shape were introduced as the model canopy. Normalized roughness length (Z ₀/H) and shear velocity ratio (R) were analyzed as a function of roughness density (λ). Experiments showed that Z ₀/H increases and R decreases as λ reaches a maximum value, λ _max, while the values of Z ₀/H and R showed little change with λ value beyond as λ _max.     

新疆西天山松湖铁矿床火山岩年代学和岩石地球化学

新疆西天山阿吾拉勒铁成矿带内众多铁矿床的围岩火山岩时代、岩石成因和构造背景迄今尚未得到很好地约束。文章对该铁成矿带内的松湖铁矿床围岩粗面英安岩和流纹岩进行了锆石LA_ICP_MS U_Pb测年和地球化学研究。粗面英安岩和流纹岩的锆石206Pb/238U加权平均年龄分别为(326.8±2.7)Ma和(327.3±1.7)Ma。粗面英安岩和流纹岩准铝质,均富集轻稀土元素、Rb、K、Zr、Hf,显示中等Eu负异常,亏损Ba、Sr、P、Ta、Nb和Ti。锆石饱和温度计算结果显示,粗面英安岩的母岩浆温度较高(774~812℃,平均792℃),流纹岩母岩浆的温度较低(713~790℃,平均750℃)。粗面英安岩具有非常低Sr高Yb的特征,Nb/Ta比值为10.8~11.4,具有较高的w(Th)(≥8.1×10-6)和高的Th/Ce比值(≥0.31);流纹岩具有低Sr低Yb特征,Nb/Ta比值为8.5~9.7,w(Th)(≥5.3×10-6)较高,Th/Ce比值为0.14~0.75,据此推测粗面英安岩母岩浆可能是上地壳部分熔融的产物,流纹岩母岩浆可能源于比粗面英安岩母岩浆更深的地壳的部分熔融。松湖铁矿区粗面英安岩和流纹岩具有弧火山岩的地球化学特征,结合构造环境判别图解和区域地质情况,推断其形成于大陆弧环境。     

Si and O self-diffusion in hydrous forsterite and iron-bearing olivine from the perspective of defect chemistry

We discuss the experimental results of silicon and oxygen self-diffusion coefficients in forsterite and iron-bearing olivine from the perspective of defect chemistry. Silicon diffusion is dominated by V_O ^··-associated V_Si″″, whereas oxygen diffusion is dominated by hopping of V_O ^·· under anhydrous conditions, and by (OH)_O ^· under hydrous conditions. By considering the charge neutrality condition of [(OH)_O ^·] = 2[V_Me″] in hydrous forsterite and iron-bearing olivine, we get D _Si ∝ (\(C_{{{\text{H}}_{2} {\text{O}}}}\))^1/3 and D _O ∝ (\(C_{{{\text{H}}_{2} {\text{O}}}}\))⁰, which explains the experimental results of water effects on oxygen and silicon self-diffusion rates (Fei et al. in Nature 498:213–215, 2013; J Geophys Res 119:7598–7606, 2014). The \(C_{{{\text{H}}_{2} {\text{O}}}}\) dependence of creep rate in the Earth’s mantle should be close to that given by Si and O self-diffusion coefficients obtained under water unsaturated conditions.     

Hydrogen sites in the dense hydrous magnesian silicate phase E: a pulsed neutron powder diffraction study

Hydrogen site positions and occupancy in the crystal structure of dense hydrous magnesium silicate (DHMS) phase E were determined for the first time by pulsed neutron powder diffraction. A fully deuterated pure phase E powder sample, which had space group \(R\overline{3} m\) and lattice parameters of a = 2.97065(8) Å and c = 13.9033(4) Å, was synthesized at 15 GPa and 1100 °C. Through quantitative evaluation of refined structure parameters obtained with sufficient spatial resolution and very high signal-to-background ratio, we conclude that the O–D dipoles in the refined phase E structure are tilted by 24° from the direction normal to the layers of edge-shared MgO₆ octahedra (octahedral layers). The tilted dipole structure of phase E is in remarkable contrast to that of brucite, Mg(OH)₂, which has dipoles exactly normal to the octahedral layer. This contrast exists because the O–Si–O bonding unique in the phase E structure connects two adjacent octahedral layers and thereby reduces the interlayer O···O distance. This shrinkage of the interlayer distance induces the tilting of the O–D dipole and also generates unique O–D···O hydrogen bonding connecting all the layers in the phase E structure.     

Thermal equation of state of natural tourmaline at high pressure and temperature

Synchrotron-based in situ angle-dispersive X-ray diffraction experiments were conducted on a natural uvite-dominated tourmaline sample by using an external-heating diamond anvil cell at simultaneously high pressures and temperatures up to 18 GPa and 723 K, respectively. The angle-dispersive X-ray diffraction data reveal no indication of a structural phase transition over the P–T range of the current experiment in this study. The pressure–volume–temperature data were fitted by the high-temperature Birch–Murnaghan equation of state. Isothermal bulk modulus of K ₀ = 96.6 (9) GPa, pressure derivative of the bulk modulus of \(K_{0}^{\prime } = 12.5 \;(4)\), thermal expansion coefficient of α ₀ = 4.39 (27) × 10^?5 K^?1 and temperature derivative of the bulk modulus (?K/?T) _P  = ?0.009 (6) GPa K^?1 were obtained. The axial thermoelastic properties were also obtained with K _a0 = 139 (2) GPa, \(K_{a0}^{\prime }\) = 11.5 (7) and α _a0 = 1.00 (11) × 10^?5 K^?1 for the a-axis, and K _c0 = 59 (1) GPa, \(K_{c0}^{\prime }\) = 11.4 (5) and α _c0 = 2.41 (24) × 10^?5 K^?1 for the c-axis. Both of axial compression and thermal expansion exhibit large anisotropic behavior. Thermoelastic parameters of tourmaline in this study were also compared with that of the other two ring silicates of beryl and cordierite.     

河南老湾金矿床上上河矿段矿床地质和成矿流体地球化学

河南桐柏老湾金矿床是桐柏-大别山(北坡)金银成矿带内大型造山带型金矿床之一。文章对该矿床的上上河矿段进行了矿床地质和成矿流体地球化学研究,旨在查明该矿段的流体成矿过程。根据矿脉穿插关系、矿石结构构造、矿物共生组合以及黄铁矿的粒度和晶形,将老湾金矿上上河矿段成矿过程划分为:石英粗粒自形黄铁矿(Ⅰ)、石英细粒半自形-他形黄铁矿(Ⅱ)、石英多金属硫化物(Ⅲ)及石英碳酸盐(Ⅳ)4个阶段。镜下观察显示,矿床中的包裹体类型有含CO_2包裹体(LH_2O+LCO_2+VCO_2)、纯CO_2包裹体(LCO_2+VCO_2)、液相包裹体(LH_2O+VH_2O)及少量含子晶包裹体(LH_2O+VH_2O+S)。第Ⅰ阶段、第Ⅱ阶段和第Ⅲ阶段均可见含CO_2包裹体、纯CO_2包裹体和液相包裹体,有时可见含CO_2包裹体与液相包裹体共存。流体包裹体显微测温结果表明,成矿流体可近似看做中温、低盐度、富CO_2的NaCl-H_2O-CO_2体系,纯CO_2包裹体和液相包裹体所代表的流体可能是由含CO_2包裹体所代表的Na Cl-H_2O-CO_2流体经不混溶形成的,三者在寄主矿物沉淀时,被同时捕获而共存。从第Ⅰ阶段到第Ⅳ阶段,成矿流体温度从303~379℃逐渐降低到138~195℃,盐度w(Na Cleq)从4.07%~9.59%逐渐降低到1.06%~2.74%。在成矿的第Ⅰ阶段成矿流体发生了不混溶作用,而在第Ⅱ阶段和第Ⅲ阶段流体中的CO_2起泡分离再次引发了不混溶作用。从第Ⅰ阶段到第Ⅲ阶段,成矿流体的δ18OH_2O从6.56‰~9.71‰经1.89‰~4.01‰变化到0.08‰,δDH_2O从-78.1‰~-64.2‰经-79.5‰~-76.3‰变化到-72.6‰,表明老湾金矿第Ⅰ阶段成矿流体主要为岩浆热液,第Ⅱ阶段成矿流体中有少量大气降水加入,第Ⅲ阶段成矿流体中大气降水的比例明显加大。     

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.     

In situ high-temperature X-ray diffraction and spectroscopic study of fibroferrite,FeOH(SO<Subscript>4</Subscript>)·5H<Subscript>2</Subscript>O

The thermal dehydration process of fibroferrite, FeOH(SO₄)·5H₂O, a secondary iron-bearing hydrous sulfate, was investigated by in situ high-temperature synchrotron X-ray powder diffraction (HT-XRPD), in situ high-temperature Fourier transform infrared spectroscopy (HT-FTIR) and thermal analysis (TGA-DTA) combined with evolved gas mass spectrometry. The data analysis allowed the determination of the stability fields and the reaction paths for this mineral as well as characterization of its high-temperature products. Five main endothermic peaks are observed in the DTA curve collected from room T up to 800 °C. Mass spectrometry of gases evolved during thermogravimetric analysis confirms that the first four mass loss steps are due to water emission, while the fifth is due to a dehydroxylation process; the final step is due to the decomposition of the remaining sulfate ion. The temperature behavior of the different phases occurring during the heating process was analyzed, and the induced structural changes are discussed. In particular, the crystal structure of a new phase, FeOH(SO₄)·4H₂O, appearing at about 80 °C due to release of one interstitial H₂O molecule, was solved by ab initio real-space and reciprocal-space methods. This study contributes to further understanding of the dehydration mechanism and thermal stability of secondary sulfate minerals.