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.     

Characteristics of extreme dust events observed over two urban areas in Iran

Determination of dust loading in the atmosphere is important not only from the public health point of view, but also for regional climate changes. The present study focuses on the characteristics of two major dust events for two urban areas in Iran, Kermanshah and Tehran, over the period of 4 years from 2006 to 2009. To detect extreme dust outbreaks, various datasets including synoptic data, dust concentration, reanalysis data and numerical results of WRF and HYSPLIT models were used. The weather maps demonstrate that for these events dusts are mainly generated when wind velocity is high and humidity is low in the lower troposphere and the region is under the influence of a thermal low. The event lasts until the atmospheric stability prevails and the surface wind speed weakens. The thermal low nature of the synoptic conditions of these major events is also responsible for deep boundary layer development with its thermals affecting the vertical dust flux over the region. Trajectory studies show that the dust events originated from deserts in Iraq and Syria and transported towards Iran. The main distinction between the two types of mobilizations seems to affect the dust concentrations in the Tehran urban area.     

Consistent seasonal snow cover depth and duration variability over the Western Himalayas (WH)

Precipitation in solid form, i.e., snow, during winter season over the Western Himalayas (WH) leads to the build-up of seasonal snow cover. Seasonal snow cover build-up (snow cover depth and duration) largely depends on atmospheric variables such as temperature, precipitation, radiation, wind, etc. Integrated (combined) influence of atmospheric variables on seasonal snow cover gets reflected in terms of spatial and temporal variability in seasonal snow cover build-up pattern. Hence spatial and temporal variability of seasonal snow cover build-up can serve as a good indicator of climate change in high altitude mountainous regions like the WH. Consistent seasonal snow cover depth and duration, delay days and early melt days of consistent seasonal snow cover at 11 stations spread across different mountain ranges over the WH were analyzed. Mean, maximum and percentiles (25th, 50th, 75th, 90th and 95th) of consistent seasonal snow cover depth and duration show decline over the WH in the recent past 2–3 decades. Consistent seasonal snow cover is found to melt early and snow cover build-up pattern is found to show changes over the WH. Decline in consistent seasonal snow cover depth, duration and changing snow cover build-up pattern over the WH in recent decades indicate that WH has undergone considerable climate change and winter weather patterns are changing in the WH.     

The interaction between surface water and groundwater and its effect on water quality in the Second Songhua River basin,northeast China

The relationship between surface water and groundwater not only influences the water quantity, but also affects the water quality. The stable isotopes (δD, δ ¹⁸O) and hydrochemical compositions in water samples were analysed in the Second Songhua River basin. The deep groundwater is mainly recharged from shallow groundwater in the middle and upper reaches. The shallow groundwater is discharged to rivers in the downstream. The runoff from upper reaches mainly contributed the river flow in the downstream. The CCME WQI indicated that the quality of surface water and groundwater was ‘Fair’. The mixing process between surface water and groundwater was simulated by the PHREEQC code with the results from the stable isotopes. The interaction between surface water and groundwater influences the composition of ions in the mixing water, and further affects the water quality with other factors.     

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.     

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.     

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.