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.     

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.     

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.     

Ionic conductivity in gem-quality single-crystal alkali feldspar from the Eifel: temperature,orientation and composition dependence

We measured the ion conductivity of single-crystal alkali feldspar originating from two different locations in the Eifel/Germany, named Volkesfeld and Rockeskyller sanidine and having potassium site fractions \(C_\mathrm{K}\) of 0.83 and 0.71, respectively. The dc conductivities resulting from electrochemical impedance spectroscopy over the temperature range of 300–900 \(^{\circ }\hbox {C}\) show a weak composition dependence but pronounced differences between the b-direction [\(\perp (010)\)] and \(c^{*}\)-direction [\(\perp (001)\)] of the monoclinic feldspar structure. Conductivity activation energies obtained from the observed linear Arrhenius plots are close to 1.2 eV in all cases, which is closely similar to the activation energies of the \(^{22}\mathrm{Na}\) tracer diffusivity in the same crystals. Taking into account literature data on K tracer diffusion and diffusion correlation effects, the present results point to a predominance of the interstitialcy mechanism over the vacancy mechanism in mass and charge transport on the alkali sublattice in potassium-rich alkali feldspar.     

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.     

Anisotropic mechanical properties of zircon and the effect of radiation damage

This study provides new insights into the relationship between radiation-dose-dependent structural damage due to natural U and Th impurities and the anisotropic mechanical properties (Poisson’s ratio, elastic modulus and hardness) of zircon. Natural zircon samples from Sri Lanka (see Muarakami et al. in Am Mineral 76:1510–1532, 1991) and synthetic samples, covering a dose range of zero up to 6.8 × 10¹⁸ α-decays/g, have been studied by nanoindentation. Measurements along the [100] crystallographic direction and calculations, based on elastic stiffness constants determined by Özkan (J Appl Phys 47:4772–4779, 1976), revealed a general radiation-induced decrease in stiffness (~54 %) and hardness (~48 %) and an increase in the Poisson’s ratio (~54 %) with increasing dose. Additional indentations on selected samples along the [001] allowed one to follow the amorphization process to the point that the mechanical properties are isotropic. This work shows that the radiation-dose-dependent changes of the mechanical properties of zircon can be directly correlated with the amorphous fraction as determined by previous investigations with local and global probes (Ríos et al. in J Phys Condens Matter 12:2401–2412, 2000a; Farnan and Salje in J Appl Phys 89:2084–2090, 2001; Zhang and Salje in J Phys Condens Matter 13:3057–3071, 2001). The excellent agreement, revealed by the different methods, indicates a large influence of structural and even local phenomena on the macroscopic mechanical properties. Therefore, this study indicates the importance of acquiring better knowledge about the mechanical long-term stability of radiation-damaged materials.