Hydrochemical characteristics of surface and groundwater and suitability for drinking and agricultural use in the Upper Tigris River Basin,Diyarbakır–Batman,Turkey

Investigations were undertaken into the quality of surface water and groundwater bodies within the Upper Tigris Basin in Turkey to determine their suitability for potable and agricultural use. In the study area, the majority of the groundwater and surface water samples belong to the calcium–magnesium–bicarbonate type (Ca–Mg–HCO₃) or magnesium–calcium–bicarbonate type (Mg–Ca–HCO₃). Chemical analysis of all water samples shows that the mean cation concentrations (in mg/L) were in the order Ca²⁺ > Mg²⁺ > Na⁺ > K⁺ and that of anions are in the order \( \text{HCO}_{3}^{ - } \) > \( \text{SO}_{4}^{2 - } \) > Cl^? > \( \text{CO}_{3}^{ - } \) for all groundwater and surface water samples. The Mg²⁺/Ca²⁺ ratio ranges from 0.21 to 1.30 with most of the values greater than 0.5, indicating that weathering of dolomites is dominant in groundwater. The analysis reveals that all of the samples are neutral to slightly alkaline (pH 7.0–8.7). Groundwater and surface water suitability for drinking usage was evaluated according to the World Health Organization and Turkish Standards (TSE-266) and suggests that most of the samples are suitable for drinking. Various determinants such as sodium absorption ratio, percent sodium (Na %), residual sodium carbonate and soluble sodium percentage revealed that most of the samples are suitable for irrigation. According to MH values, all of the well water samples were suitable for irrigation purposes, but 80 and 81.82% of Zillek springs and surface water samples were unsuitable. As per the PI values, the water samples from the study area are classified as Class I and Class II and are considered to be suitable for irrigation.     

Appraisal of groundwater quality in a crystalline aquifer: a chemometric approach

The purpose of this study is to assess the groundwater quality and identify the processes that control the groundwater chemistry in a crystalline aquifer. A total of 72 groundwater samples were collected during pre- and post-monsoon seasons in the year 2014 in a semi-arid region of Gooty Mandal, Anantapur district, Andhra Pradesh, India. The study utilized chemometric analysis like basic statistics, Pearson’s correlation coefficient (r), principal component analysis (PCA), Gibbs ratio, and index of base exchange to understand the mechanism of controlling the groundwater chemistry in the study area. The results reveal that groundwater in the study area is neutral to slightly alkaline in nature. The order of dominance of cations is Na⁺ > Ca²⁺ > Mg²⁺ > K⁺ while for anions, it is \( {\mathrm{HCO}}_3^{-}>{\mathrm{Cl}}^{-} \)>\( {\mathrm{NO}}_3^{-} \)>\( {\mathrm{SO}}_4^{2-} \)>\( {\mathrm{CO}}_3^{2-}>{\mathrm{F}}^{-} \) in both seasons. Based on the Piper classification, most of the groundwater samples are identified as of sodium bicarbonate (\( {\mathrm{Na}}^{+}-{\mathrm{HCO}}_3^{-}\Big) \) type. According to the results of the principal component analysis (PCA), three factors and two factors were identified pre and post monsoon, respectively. The present study indicates that the groundwater chemistry is mostly controlled by geogenic processes (weathering, dissolution, and ion exchange) and some extent of anthropogenic activities.     

(Ca,Mg)2Si2O6 clinopyroxenes: A solution model based on nonconvergent site-disorder

Experiments at high pressure and temperature indicate that excess Ca may be dissolved in diopside. If the (Ca, Mg)₂Si₂O₆ clinopyroxene solution extends to more Ca-rich compositions than CaMgSi₂O₆, macroscopic regular solution models cannot strictly be applied to this system. A nonconvergent site-disorder model, such as that proposed by Thompson (1969, 1970), may be more appropriate. We have modified Thompson's model to include asymmetric excess parameters and have used a linear least-squares technique to fit the available experimental data for Ca-Mg orthopyroxene-clinopyroxene equilibria and Fe-free pigeonite stability to this model. The model expressions for equilibrium conditions \(\mu _{{\text{Mg}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{opx}}} = \mu _{{\text{Mg}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{cpx}}} \) (reaction A) and \(\mu _{{\text{Ca}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{opx}}} = \mu _{{\text{Ca}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{cpx}}} \) (reaction B) are given by: 1 $$\begin{gathered} \Delta \mu _{\text{A}}^{\text{O}} = {\text{RT 1n}}\left[ {\frac{{(X_{{\text{Mg}}}^{{\text{opx}}} )^2 }}{{X_{{\text{Mg}}}^{{\text{M1}}} \cdot X_{{\text{Mg}}}^{{\text{M2}}} }}} \right] - \frac{1}{2}\{ W_{21} [2(X_{{\text{Ca}}}^{{\text{M2}}} )^3 - (X_{{\text{Ca}}}^{{\text{M2}}} ] \hfill \\ {\text{ + 2W}}_{{\text{22}}} [X_{{\text{Ca}}}^{{\text{M2}}} )^2 - (X_{{\text{Ca}}}^{{\text{M2}}} )^3 + \Delta {\text{G}}_{\text{*}}^{\text{0}} (X_{{\text{Ca}}}^{{\text{M1}}} \cdot X_{{\text{Ca}}}^{{\text{M2}}} )\} \hfill \\ {\text{ + W}}^{{\text{opx}}} (X_{{\text{Wo}}}^{{\text{opx}}} )^2 \hfill \\ \Delta \mu _{\text{B}}^{\text{O}} = {\text{RT 1n}}\left[ {\frac{{(X_{{\text{Ca}}}^{{\text{opx}}} )^2 }}{{X_{{\text{Ca}}}^{{\text{M1}}} \cdot X_{{\text{Ca}}}^{{\text{M2}}} }}} \right] - \frac{1}{2}\{ 2W_{21} [2(X_{{\text{Mg}}}^{{\text{M2}}} )^2 - (X_{{\text{Mg}}}^{{\text{M2}}} )^3 ] \hfill \\ {\text{ + W}}_{{\text{22}}} [2(X_{{\text{Mg}}}^{{\text{M2}}} )^3 - (X_{{\text{Mg}}}^{{\text{M2}}} )^2 + \Delta {\text{G}}_{\text{*}}^{\text{0}} (X_{{\text{Mg}}}^{{\text{M1}}} \cdot X_{{\text{Mg}}}^{{\text{M2}}} )\} \hfill \\ {\text{ + W}}^{{\text{opx}}} (X_{{\text{En}}}^{{\text{opx}}} )^2 \hfill \\ \hfill \\ \end{gathered} $$ where 1 $$\begin{gathered} \Delta \mu _{\text{A}}^{\text{O}} = 2.953 + 0.0602{\text{P}} - 0.00179{\text{T}} \hfill \\ \Delta \mu _{\text{B}}^{\text{O}} = 24.64 + 0.958{\text{P}} - (0.0286){\text{T}} \hfill \\ {\text{W}}_{{\text{21}}} = 47.12 + 0.273{\text{P}} \hfill \\ {\text{W}}_{{\text{22}}} = 66.11 + ( - 0.249){\text{P}} \hfill \\ {\text{W}}^{{\text{opx}}} = 40 \hfill \\ \Delta {\text{G}}_*^0 = 155{\text{ (all values are in kJ/gfw)}}{\text{.}} \hfill \\ \end{gathered} $$ . Site occupancies in clinopyroxene were determined from the internal equilibrium condition 1 $$\begin{gathered} \Delta G_{\text{E}}^{\text{O}} = - {\text{RT 1n}}\left[ {\frac{{X_{{\text{Ca}}}^{{\text{M1}}} \cdot X_{{\text{Mg}}}^{{\text{M2}}} }}{{X_{{\text{Ca}}}^{{\text{M2}}} \cdot X_{{\text{Mg}}}^{{\text{M1}}} }}} \right] + \tfrac{1}{2}[(2{\text{W}}_{{\text{21}}} - {\text{W}}_{{\text{22}}} )(2{\text{X}}_{{\text{Ca}}}^{{\text{M2}}} - 1) \hfill \\ {\text{ + }}\Delta G_*^0 (X_{{\text{Ca}}}^{{\text{M1}}} - X_{{\text{Ca}}}^{{\text{M2}}} ) + \tfrac{3}{2}(2{\text{W}}_{{\text{21}}} - {\text{W}}_{{\text{22}}} ) \hfill \\ {\text{ (1}} - 2X_{{\text{Ca}}}^{{\text{M1}}} )(X_{{\text{Ca}}}^{{\text{M1}}} + \tfrac{1}{2})] \hfill \\ \end{gathered} $$ where δG _E ⁰ =153+0.023T+1.2P. The predicted concentrations of Ca on the clinopyroxene Ml site are low enough to be compatible with crystallographic studies. Temperatures calculated from the model for coexisting ortho- and clinopyroxene pairs fit the experimental data to within 10° in most cases; the worst discrepancy is 30°. Phase relations for clinopyroxene, orthopyroxene and pigeonite are successfully described by this model at temperatures up to 1,600° C and pressures from 0.001 to 40 kbar. Predicted enthalpies of solution agree well with the calorimetric measurements of Newton et al. (1979). The nonconvergent site disorder model affords good approximations to both the free energy and enthalpy of clinopyroxenes, and, therefore, the configurational entropy as well. This approach may provide an example for Febearing pyroxenes in which cation site exchange has an even more profound effect on the thermodynamic properties.     

Energetics of hydration of cordierite and water barometry in cordierite-granulites

The Gibbs free energy and volume changes attendant upon hydration of cordierites in the system magnesian cordierite-water have been extracted from the published high pressure experimental data at \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) =P _total, assuming an ideal one site model for H₂O in cordierite. Incorporating the dependence of ΔG and ΔV on temperature, which was found to be linear within the experimental conditions of 500°–1,000°C and 1–10,000 bars, the relation between the water content of cordierite and P, T and \(f_{{\text{H}}_{\text{2}} {\text{O}}} \) has been formulated as $$\begin{gathered} X_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{crd}}} = \hfill \\ \frac{{f_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{P, T}}} }}{{\left[ {{\text{exp}}\frac{1}{{RT}}\left\{ {64,775 - 32.26T + G_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{1, }}T} - P\left( {9 \times 10^{ - 4} T - 0.5142} \right)} \right\}} \right] + f_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{P, T}}} }} \hfill \\ \end{gathered} $$ The equation can be used to compute H₂O in cordierites at \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) <1. Our results at different P, T and partial pressure of water, assuming ideal mixing of H₂O and CO₂ in the vapour phase, are in very good agreement with the experimental data of Johannes and Schreyer (1977, 1981). Applying the formulation to determine \(X_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{crd}}} \) in the garnet-cordierite-sillimanite-plagioclase-quartz granulites of Finnish Lapland as a test case, good agreement with the gravimetrically determined water contents of cordierite was obtained. Pressure estimates, from a thermodynamic modelling of the Fe-cordierite — almandine — sillimanite — quartz equilibrium at \(P_{{\text{H}}_{\text{2}} {\text{O}}} = 0\) and \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) =P_total, for assemblages from South India, Scottish Caledonides, Daly Bay and Hara Lake areas are compatible with those derived from the garnetplagioclase-sillimanite-quartz geobarometer.     

Zinc isotope fractionation under vaporization processes and in aqueous solutions

Equilibrium Zn isotope fractionation was investigated using first-principles quantum chemistry methods at the B3LYP/6-311G* level. The volume variable cluster model method was used to calculate isotope fractionation factors of sphalerite, smithsonite, calcite, anorthite, forsterite, and enstatite. The water-droplet method was used to calculate Zn isotope fractionation factors of Zn²⁺-bearing aqueous species; their reduced partition function ratio factors decreased in the order \(\left[ {{\text{Zn}}\left( {{\text{H}}_{2} {\text{O}}} \right)_{6} } \right]^{2 + } > \left[ {{\text{ZnCl}}\left( {{\text{H}}_{2} {\text{O}}} \right)_{5} } \right]^{ + } > \left[ {{\text{ZnCl}}_{2} \left( {{\text{H}}_{2} {\text{O}}} \right)_{4} } \right] > \left[ {{\text{ZnCl}}_{3} \left( {{\text{H}}_{2} {\text{O}}} \right)_{2} } \right]^{ - } > {\text{ZnCl}}_{4} ]^{2 - }\). Gaseous ZnCl₂ was also calculated for vaporization processes. Kinetic isotope fractionation of diffusional processes in a vacuum was directly calculated using formulas provided by Richter and co-workers. Our calculations show that in addition to the kinetic isotope effect of diffusional processes, equilibrium isotope fractionation also contributed nontrivially to observed Zn isotope fractionation of vaporization processes. The calculated net Zn isotope fractionation of vaporization processes was 7–7.5‰, with ZnCl₂ as the gaseous species. This matches experimental observations of the range of Zn isotope distribution of lunar samples. Therefore, vaporization processes may be the cause of the large distribution of Zn isotope signals found on the Moon. However, we cannot further distinguish the origin of such vaporization processes; it might be due either to igneous rock melting in meteorite bombardments or to a giant impact event. Furthermore, isotope fractionation between Zn-bearing aqueous species and minerals that we have provided helps explain Zn isotope data in the fields of ore deposits and petrology.     

Geochemical assessment of fluoride enrichment and nitrate contamination in groundwater in hard-rock aquifer by using graphical and statistical methods

This systematic study was carried out with objective to delineate the various sources responsible for \(\hbox {NO}_{3}^{-}\) contamination and \(\hbox {F}^{-}\) enrichment by utilizing statistical and graphical methods. Since Central Ground Water Board, India, indicated susceptibility of \(\hbox {NO}_{3}^{-}\) contamination and \(\hbox {F}^{-}\) enrichment, in most of the groundwater, \(\hbox {NO}_{3}^{-}\) and \(\hbox {F}^{-}\) concentration primarily observed \({>}45\) and \({>}1.5~\hbox {mg/L}\), respectively, i.e., higher than the permissible limit for drinking water. Water Quality Index (WQI) indicates \({\sim }22.81\%\) groundwater are good-water, \({\sim }71.14\%\) groundwater poor-water, \({\sim }5.37\%\) very poor-water and 0.67% unsuitable for drinking purpose. Piper diagram indicates \({\sim }59.73\%\) groundwater hydrogeochemical facies are Ca–Mg–\(\hbox {HCO}_{3 }\) water-types, \({\sim }28.19\%\) Ca–Mg–\(\hbox {SO}_{4}\)–Cl water-types, \({\sim }8.72\%\) Na–K–\(\hbox {SO}_{4}\)–Cl water-types and 3.36% Na–K–\(\hbox {HCO}_{3 }\) water-types. This classification indicates dissolution and mixing are mainly controlling groundwater chemistry. Salinity diagram indicate \({\sim }44.30\%\) groundwater under in low sodium and medium salinity hazard, \({\sim }49.66\%\) groundwater fall under low sodium and high salinity hazard, \({\sim }3.36\%\) groundwater fall under very-high salinity hazard. Sodium adsorption ratio indicates \({\sim }97\%\) groundwater are in excellent condition for irrigation. The spatial distribution of \(\hbox {NO}_{3}^{-}\) indicates significant contribution of fertilizer from agriculture lands. Fluoride enrichment occurs in groundwater through the dissolution of fluoride-rich minerals. By reducing the consumption of fertilizer and stress over groundwater, the water quality can be improved.     

Fluvial geochemistry of Subarnarekha River basin,India

The fluvial geochemistry of the Subarnarekha River and its major tributaries has been studied on a seasonal basis in order to assess the geochemical processes that explain the water composition and estimate solute fluxes. The analytical results show the mildly acidic to alkaline nature of the Subarnarekha River water and the dominance of \(\hbox {Ca}^{2+}\) and \(\hbox {Na}^{+}\) in cationic and \(\hbox {HCO}_{3}^{-}\) and \({\hbox {Cl}}^{-}\) in anionic composition. Minimum ionic concentration during the monsoon and maximum concentration in the pre-monsoon seasons reflect concentrating effects due to decrease in the river discharge and increase in the base flow contribution during the pre-monsoon and dilution effects of atmospheric precipitation in the monsoon season. The solute acquisition processes are mainly controlled by weathering of rocks, with minor contribution from marine and anthropogenic sources. Higher contribution of alkaline earth \((\hbox {Ca}^{2+}{+}\,\hbox {Mg}^{2+})\) to the total cations \((\hbox {TZ}^{+})\) and high \((\hbox {Na}^{+}+\hbox {K}^{+})/\hbox {Cl}^{-}\), \((\hbox {Na}^{+}+\hbox {K}^{+})/\hbox {TZ}^{+}\), \(\hbox {HCO}_{3}^{-}/(\hbox {SO}_{4}^{2-}+\hbox {Cl}^{-})\) and low \((\hbox {Ca}^{2+}+\hbox {Mg}^{2+})/(\hbox {Na}^{+}+\hbox {K}^{+})\) equivalent ratios suggest that the Subarnarekha River water is under the combined influence of carbonate and silicate weathering. The river water is undersaturated with respect to dolomite and calcite during the post-monsoon and monsoon seasons and oversaturated in the pre-monsoon season. The pH–log \(\hbox {H}_{4}\hbox {SiO}_{4}\) stability diagram demonstrates that the water chemistry is in equilibrium with the kaolinite. The Subarnarekha River annually delivered \(1.477\times 10^{6}\) ton of dissolved loads to the Bay of Bengal, with an estimated chemical denudation rate of \(77\hbox { ton km}^{-2}\hbox { yr}^{-1}\). Sodium adsorption ratio, residual sodium carbonate and per cent sodium values placed the studied river water in the ‘excellent to good quality’ category and it can be safely used for irrigation.     

Chemistry, Raman and infrared spectroscopic characterization of the phosphate mineral reddingite: (MnFe)3(PO4)2(H2O,OH)3, a mineral found in lithium-bearing pegmatite

Detailed investigation of an intermediate member of the reddingite–phosphoferrite series, using infrared and Raman spectroscopy, scanning electron microcopy and electron microprobe analysis, has been carried out on a homogeneous sample from a lithium-bearing pegmatite named Cigana mine, near Conselheiro Pena, Minas Gerais, Brazil. The determined formula is $ ({\text{Mn}}_{1.60} {\text{Fe}}_{1.21} {\text{Ca}}_{0.01} {\text{Mg}}_{0.01} )_{\sum 2.83} ({\text{PO}}_{4} )_{2.12} \cdot ({\text{H}}_{2} {\text{O}}_{2.85} {\text{F}}_{0.01} )_{\sum 2.86} $ , indicating predominance in the reddingite member. Raman spectroscopy coupled with infrared spectroscopy supports the concept of phosphate, hydrogen phosphate and dihydrogen phosphate units in the structure of reddingite-phosphoferrite. Infrared and Raman bands attributed to water and hydroxyl stretching modes are identified. Vibrational spectroscopy adds useful information to the molecular structure of reddingite–phosphoferrite.     

Adsorption kinetics of chromium(III) removal from aqueous solutions using natural red earth

Laboratory-scale-simulated experiments were carried out using Cr(III) solutions to identify the Cr(III) retention behavior of natural red earth (NRE), a natural soil available in the northwestern coastal belt of Sri Lanka. The effects of solution pH, initial Cr(III) concentration and the contact time were examined. The NRE showed almost 100 % Cr(III) adsorption within the first 90 min. [initial [Cr(III)] = 0.0092–0.192 mM; initial pH 4.0–9.0]. At pH 2 (298 K), when particle size ranged from 125 to 180 μm the Cr(III) adsorption data were modeled according to Langmuir convention assuming site homogeneity. The pH-dependent Cr(III) adsorption data were quantified by diffused layer model assuming following reaction stoichiometries: $$ \begin{aligned} 2\, {>}{\text{AlOH}}_{{({\text{s}})}} + {\text{ Cr }}\left( {\text{OH}} \right)_{{ 2\,({\text{aq}})}}^{ + } \, \to \, \left( { {>}{\text{AlO}}} \right)_{ 2} {\text{Cr}}_{{({\text{s}})}}^{ + } + {\text{ 2H}}_{ 2} {\text{O}} \quad {\text{log K 15}}. 5 6\\ 2\, {>}{\text{FeOH}}_{{({\text{s}})}} + {\text{ Cr}}\left( {\text{OH}} \right)_{{ 2\,({\text{aq}})}}^{ + } \, \to \, \left( { {>}{\text{FeO}}} \right)_{ 2} {\text{Cr}}_{{({\text{s}})}}^{ + } + {\text{ 2H}}_{ 2} {\text{O}}\quad {\text{log K 5}}.0 8.\\ \end{aligned} $$ The present data showed that NRE can effectively be used to mitigate Cr(III) from aqueous solutions and this method is found to be simple, effective, economical and environmentally benign.     

Assessment of groundwater quality for drinking and irrigation purposes using hydrochemical studies in Nalbari district of Assam,India

The present work was carried out in Nalbari district of Assam (India) with an objective to assess the quality of groundwater and to check its suitability for drinking and irrigation purposes. Groundwater samples were collected from 50 different locations during pre- and post-monsoon seasons of 2016. Results of chemical analysis revealed that mean concentration of cations varied in the order Ca²⁺?>?Na⁺?>?Mg²⁺?>?K⁺, while for anions the order was HCO₃ ^??>?Cl^??>?SO₄^2??>?NO₃^2??>?F^? during both pre- and post-monsoon seasons. The suitability of groundwater samples for drinking purpose was assessed by comparing the results of physico-chemical analysis of groundwater with Indian Standards. Further, its suitability for irrigation purpose was assessed by evaluating several parameters like sodium adsorption ratio (SAR), sodium percentage (Na%), magnesium ratio, Kelly’s ratio and residual sodium carbonate (RSC). The SAR values obtained for all the samples were plotted against EC values in the US Salinity Laboratory diagram, and it was revealed that the most of the samples fall under water type C2-S1 indicating medium salinity and low SAR. Further, it was found that the majority of the samples belong to Ca–Mg–HCO₃ hydrochemical facies followed by Ca–Mg–Cl–SO₄, whereas only a few samples belong to Na–K–HCO₃ hydrochemical facies.     

The crystal structure of waylandite from Wheal Remfry, Cornwall, United Kingdom

Sr- and Ca-rich waylandite, $ {\left( {{\hbox{B}}{{\hbox{i}}_{0.{54}}}{\hbox{S}}{{\hbox{r}}_{0.{31}}}{\hbox{C}}{{\hbox{a}}_{0.{25}}}{{\hbox{K}}_{0.0{1}}}{\hbox{B}}{{\hbox{a}}_{0.0{1}}}} \right)_{\Sigma 1.12}}{{\hbox{H}}_{0.{18}}}{\left( {{\hbox{A}}{{\hbox{l}}_{{2}.{96}}}{\hbox{C}}{{\hbox{u}}_{0.0{2}}}} \right)_{\Sigma 2.98}}{\left[ {{{\left( {{{\hbox{P}}_{0.{97}}}{{\hbox{S}}_{0.0{3}}}{\hbox{S}}{{\hbox{i}}_{0.0{1}}}} \right)}_{\Sigma 1.00}}{{\hbox{O}}_4}} \right]_2}{\left( {\hbox{OH}} \right)_6} $ , from Wheal Remfry, Cornwall, United Kingdom has been investigated by single-crystal X-ray diffraction and electron microprobe analyses. Waylandite crystallises in space group R $ \overline 3 $ ? m, with the cell parameters: a?=?7.0059(7) Å, c?=?16.3431(12) Å and V?=?694.69(11) ų. The crystal structure has been refined to R ₁?=?3.76%. Waylandite has an alunite-type structure comprised of a rhombohedral stacking of (001) composite layers of corner-shared AlO₆ octahedra and PO₄ tetrahedra, with (Bi,Sr,Ca) atoms occupying icosahedrally coordinated sites between the layers.     

Influence of ignimbrite on the chemistry of river water in Shirasu plateau,Japan

River, rain and spring water samples from a region covered in “Shirasu” ignimbrite were collected on Kyushu Island, Japan. The analytical results were subjected to multivariate statistical analysis and stoichiometric calculation to understand the geographical distribution of chemical components in water and to extract geochemical underlying factors. The multivariate statistical analysis showed that the river-water chemistry is only slightly influenced by hot springs or polluted waters, but is highly controlled by weathering of ignimbrite. On the basis of the stoichiometric calculation based on water–rock interaction, the water chemistry was successfully estimated by a simple equation:\({\left[ {{\text{Si}}} \right]}{\text{ = 2}}{\left[ {{\text{Na}}^{{\text{ + }}} } \right]}{\text{ + }}{\left[ {{\text{Mg}}^{{{\text{2 + }}}} } \right]}\) in the upstream area, complemented by \({\left[ {{\text{Si}}} \right]}{\text{ = }}{\left[ {{\text{Na}}^{{\text{ + }}} } \right]}{\text{ - 3}}{\left[ {{\text{K}}^{{\text{ + }}} } \right]}{\text{ + }}{\left[ {{\text{Mg}}^{{{\text{2 + }}}} } \right]}{\text{ - 2}}{\left[ {{\text{Ca}}^{{{\text{2 + }}}} } \right]}\) in the downstream area.     

Majoritic garnet grains within shock-induced melt veins in amphibolites from the Ries impact crater suggest ultrahigh crystallization pressures between 18 and 9 GPa

Shock-induced melt veins in amphibolites from the Nördlinger Ries often have chemical compositions that are similar to bulk rock (i.e., basaltic), but there are other veins that are confined to chlorite-rich cracks that formed before the impact and these are poor in Ca and Na. Majoritic garnets within the shock veins show a broad chemical variation between three endmembers: (1) \({}^{\text{VIII}}{{\text{M}^{2+}}_3} {}^{\text{VI}}{\text{Al}}_{2} ({}^{\text{IV}}{\text{SiO}}_{4} )_{3}\) (normal garnet, Grt), (2) \({}^{\text{VIII}}{{\text{M}^{2+}}_3} {}^{\text{VI}}[{\text{M}}^{2 + } ({\text{Si,Ti}})]({}^{\text{IV}}{\text{SiO}}_{4} )_{3}\)  (majorite, Maj), and (3) \({}^{\text{VIII}}({{\text {Na} {\text M}^{2+}}_2}) {}^{\text{VI}}[ ({\text{Si,Ti}}){\text {Al}}]({}^{\text{IV}}{\text{SiO}}_{4} )_{3}\) (Na-majorite₅₀Grt₅₀), whereby M²⁺ = Mg²⁺, Fe²⁺, Mn²⁺, Ca²⁺. In particular, we observed a broad variation in ^VI(Si,Ti) which ranges from 0.12 to 0.58 cations per formula unit (cpfu). All these majoritic garnets crystallized during shock pressure release at different ultrahigh pressures. Those with high ^VI(Si,Ti) (0.36–0.58 cpfu) formed at high pressures and temperatures from amphibole-rich melts, while majoritic garnets with lower ^VI(Si,Ti) of 0.12–0.27 cpfu formed at lower pressures and temperatures from chlorite-rich melts. Furthermore, majoritic garnets with intermediate values of ^VI(Si,Ti) (0.24–0.39) crystallized from melts with intermediate contents of Ca and Na. To the best of our knowledge the ‘MORB-type’ Ca–Na-rich majoritic garnets with maximum contents of 2.99 wt% Na₂O and calculated crystallisation pressures of 16–18 GPa are the most extreme representatives ever found in terrestrial shocked materials. At the Ries, the duration of the initial contact and compression stage at the central location of impact is estimated to only ~ 0.1 s. We used a ~ 200-µm-thick shock-induced vein in a moderately shocked amphibolite to model its pressure–temperature–time (P–T–t) path. The graphic model manifests a peak temperature of ~ 2600 °C for the vein, continuum pressure lasting for ~ 0.02 s, a quench duration of ~ 0.02 s and a shock pulse of ~ 0.038 s. The small difference between the continuum pressure and the pressure of majoritic garnet crystallization underlines the usefulness of applying crystallisation pressures of majoritic garnets from metabasites for calculation of dynamic shock pressures of host rocks. Majoritic garnets of chlorite provenance, however, are not suitable for the determination of continuum pressure since they crystallized relatively late during shock release. An extraordinary glass- and majorite-bearing amphibole fragment in a shock-vein of one amphibolite documents the whole unloading path.     

Genesis of secondary uranium minerals associated with jasperoid veins,El Erediya area,Eastern Desert,Egypt

Uranium mineralization in the El Erediya area, Egyptian Eastern Desert, has been affected by both high temperature and low temperature fluids. Mineralization is structurally controlled and is associated with jasperoid veins that are hosted by a granitic pluton. This granite exhibits extensive alteration, including silicification, argillization, sericitization, chloritization, carbonatization, and hematization. The primary uranium mineral is pitchblende, whereas uranpyrochlore, uranophane, kasolite, and an unidentified hydrated uranium niobate mineral are the most abundant secondary uranium minerals. Uranpyrochlore and the unidentified hydrated uranium niobate mineral are interpreted as alteration products of petscheckite. The chemical formula of the uranpyrochlore based upon the Electron Probe Micro Analyzer (EPMA) is . It is characterized by a relatively high Zr content (average ZrO₂ = 6.6 wt%). The average composition of the unidentified hydrated uranium niobate mineral is , where U and Nb represent the dominant cations in the U and Nb site, respectively. Uranophane is the dominant U⁶⁺ silicate phase in oxidized zones of the jasperoid veins. Kasolite is less abundant than uranophane and contains major U, Pb, and Si but only minor Ca, Fe, P, and Zr. A two-stage metallogenetic model is proposed for the alteration processes and uranium mineralization at El Erediya. The primary uranium minerals were formed during the first stage of the hydrothermal activity that formed jasperoid veins in El Eradiya granite (130–160 Ma). This stage is related to the Late Jurassic–Early Cretaceous phase of the final Pan-African tectono-thermal event in Egypt. After initial formation of El Erediya jasperoid veins, a late stage of hydrothermal alteration includes argillization, dissolution of iron-bearing sulfide minerals, formation of iron-oxy hydroxides, and corrosion of primary uranium minerals, resulting in enrichment of U, Ca, Pb, Zr, and Si. During this stage, petscheckite was altered to uranpyrochlore and oxy-petscheckite. Uranium was likely transported as uranyl carbonate and uranyl fluoride complexes. With change of temperature and pH, these complexes became unstable and combined with silica, calcium, and lead to form uranophane and kasolite. Finally, at a later stage of low-temperature supergene alteration, oxy-petscheckite was altered to an unidentified hydrated uranium niobate mineral by removal of Fe.     

A multi-isotope approach for understanding sources of water,carbon and sulfur in natural springs of the Central Appalachian region

Natural springs have been reliable sources of domestic water and have allowed for the development of recreational facilities and resorts in the Central Appalachians. The structural history of this area is complex and it is unknown whether these natural springs receive significant recharge from modern precipitation or whether they discharge old water recharged over geological times scales. The main objective of this study was to use stable isotopes of water ( $\delta^{18} {\text{O}}_{{{\text{H}}_{2} {\text{O}}}}$ and $\delta^{2} {\text{H}}_{{{\text{H}}_{2} {\text{O}}}}$ ), dissolved inorganic carbon ( $\delta^{13} {\text{C}}_{\text{DIC}}$ ) and dissolved sulfate ( $\delta^{34} {\text{S}}_{{{\text{SO}}_{4} }}$ and $\delta^{18} {\text{O}}_{{{\text{SO}}_{4} }}$ ) to delineate sources of water, carbon and sulfur in several natural springs of the region. Our preliminary isotope data indicate that all springs are being recharged by modern precipitation. The oxygen isotope composition indicates that waters in thermal springs did not encounter the high temperatures required for O isotope exchange between the water and silicate/carbonate minerals, and/or the residence time of water in the aquifers was short due to high flow rates. The carbon isotopic composition of dissolved inorganic carbon and sulfur/oxygen isotopic composition of dissolved sulfate provide evidence of low-temperature water–rock interactions and various biogeochemical transformations these waters have undergone along their flow path.     

The H<Subscript>2</Subscript>O solubility of alkali basaltic melts: an experimental study

Experiments were conducted to determine the water solubility of alkali basalts from Etna, Stromboli and Vesuvius volcanoes, Italy. The basaltic melts were equilibrated at 1,200°C with pure water, under oxidized conditions, and at pressures ranging from 163 to 3,842 bars. Our results show that at pressures above 1 kbar, alkali basalts dissolve more water than typical mid-ocean ridge basalts (MORB). Combination of our data with those from previous studies allows the following simple empirical model for the water solubility of basalts of varying alkalinity and fO₂ to be derived: \textH ₂ \textO( \textwt% ) = \text H ₂ \textO_\textMORB ( \textwt% ) + ( 5.84 ×10 ^{- 5} *\textP - 2.29 ×10 ^{- 2} ) ×( \textNa₂ \textO + \textK₂ \textO )( \textwt% ) + 4.67 ×10 ^{- 2} ×\Updelta \textNNO - 2.29 ×10 ^{- 1} {\text{H}}_{ 2} {\text{O}}\left( {{\text{wt}}\% } \right) = {\text{ H}}_{ 2} {\text{O}}_{\text{MORB}} \left( {{\text{wt}}\% } \right) + \left( {5.84 \times 10^{ - 5} *{\text{P}} - 2.29 \times 10^{ - 2} } \right) \times \left( {{\text{Na}}_{2} {\text{O}} + {\text{K}}_{2} {\text{O}}} \right)\left( {{\text{wt}}\% } \right) + 4.67 \times 10^{ - 2} \times \Updelta {\text{NNO}} - 2.29 \times 10^{ - 1} where H₂O_MORB is the water solubility at the calculated P, using the model of Dixon et al. (1995). This equation reproduces the existing database on water solubilities in basaltic melts to within 5%. Interpretation of the speciation data in the context of the glass transition theory shows that water speciation in basalt melts is severely modified during quench. At magmatic temperatures, more than 90% of dissolved water forms hydroxyl groups at all water contents, whilst in natural or synthetic glasses, the amount of molecular water is much larger. A regular solution model with an explicit temperature dependence reproduces well-observed water species. Derivation of the partial molar volume of molecular water using standard thermodynamic considerations yields values close to previous findings if room temperature water species are used. When high temperature species proportions are used, a negative partial molar volume is obtained for molecular water. Calculation of the partial molar volume of total water using H₂O solubility data on basaltic melts at pressures above 1 kbar yields a value of 19 cm³/mol in reasonable agreement with estimates obtained from density measurements.     

Paragenesis of unusual Fe-cordierite (sekaninaite)-bearing paralava and clinker from the Kuznetsk coal basin,Siberia, Russia

Sekaninaite (X_Fe > 0.5)-bearing paralava and clinker are the products of ancient combustion metamorphism in the western part of the Kuznetsk coal basin, Siberia. The combustion metamorphic rocks typically occur as clinker beds and breccias consisting of vitrified sandstone–siltstone clinker fragments cemented by paralava, resulting from hanging-wall collapse above burning coal seams and quenching. Sekaninaite–Fe-cordierite (X_Fe = 95–45) is associated with tridymite, fayalite, magnetite, ± clinoferrosilite and ±mullite in paralava and with tridymite and mullite in clinker. Unmelted grains of detrital quartz occur in both rocks (<3 vol% in paralavas and up to 30 vol% in some clinkers). Compositionally variable siliceous, K-rich peraluminous glass is <30% in paralavas and up to 85% in clinkers. The paralavas resulted from extensive fusion of sandstone–siltstone (clinker), and sideritic/Fe-hydroxide material contained within them, with the proportion of clastic sediments ≫ ferruginous component. Calculated dry liquidus temperatures of the paralavas are 1,120–1,050°C and 920–1,050°C for clinkers, with calculated viscosities at liquidus temperatures of 10^1.6–7.0 and 10^7.0–9.8 Pa s, respectively. Dry liquidus temperatures of glass compositions range between 920 and 1,120°C (paralava) and 920–960°C (clinker), and viscosities at these temperatures are 10^9.7–5.5 and 10^8.8–9.7 Pa s, respectively. Compared with worldwide occurrences of cordierite–sekaninaite in pyrometamorphic rocks, sekaninaite occurs in rocks with X_Fe (mol% FeO/(FeO + MgO)) > 0.8; sekaninaite and Fe-cordierite occur in rocks with X_Fe 0.6–0.8, and cordierite (X_Fe < 0.5) is restricted to rocks with X_Fe < 0.6. The crystal-chemical formula of an anhydrous sekaninaite based on the refined structure is | \textK_0.02 |(\textFe_1.54^{2 +} \textMg_0.40 \textMn_0.06 )_{\Upsigma 2.00}^M [(\textAl_1.98 \textFe_0.02^{2 +} \textSi_1.00 )_{\Upsigma 3.00}^T1 (\textSi_3.94 \textAl_2.04 \textFe_0.02^{2 +} )_{\Upsigma 6.00}^T2 \textO₁₈ ]. \left| {{\text{K}}_{0.02} } \right|({\text{Fe}}_{1.54}^{2 + } {\text{Mg}}_{0.40} {\text{Mn}}_{0.06} )_{\Upsigma 2.00}^{M} [({\text{Al}}_{1.98} {\text{Fe}}_{0.02}^{2 + } {\text{Si}}_{1.00} )_{\Upsigma 3.00}^{T1} ({\text{Si}}_{3.94} {\text{Al}}_{2.04} {\text{Fe}}_{0.02}^{2 + } )_{\Upsigma 6.00}^{T2} {\text{O}}_{18} ].     

The combined effect of F and H2O on interdiffusion between peralkaline dacitic and rhyolitic melts

The effective binary diffusion coefficient (EBDC) of silicon has been measured during the interdiffusion of peralkaline, fluorine-bearing (1.3 wt% F), hydrous (3.3 and 6 wt% H₂O), dacitic and rhyolitic melts at 1.0 GPa and temperatures between 1100°C and 1400°C. From Boltzmann-Matano analysis of diffusion profiles the diffusivity of silicon at 68 wt% SiO₂ can be described by the following Arrhenius equations (with standard errors): $$\begin{gathered} {\text{with 1}}{\text{.3 wt\% F and 3}}{\text{.3\% H}}_{\text{2}} {\text{O:}} \hfill \\ {\text{D}}_{{\text{Si}}} = \begin{array}{*{20}c} { + {\text{3}}{\text{.59}}} \\ {{\text{3}}{\text{.66}} \times {\text{10}}^{ - {\text{9}}} } \\ { - {\text{1}}{\text{.86}}} \\ \end{array} {\text{exp}}\left( {{{ - {\text{86}}{\text{.1}} \pm {\text{8}}{\text{.9}}} \mathord{\left/ {\vphantom {{ - {\text{86}}{\text{.1}} \pm {\text{8}}{\text{.9}}} {{\text{RT}}}}} \right. \kern-\nulldelimiterspace} {{\text{RT}}}}} \right) \hfill \\ {\text{with 1}}{\text{.3 wt\% F and 6}}{\text{.0\% H}}_{\text{2}} {\text{O:}} \hfill \\ {\text{D}}_{{\text{Si}}} = \begin{array}{*{20}c} { + {\text{3}}{\text{.59}}} \\ {{\text{3}}{\text{.51}} \times {\text{10}}^{ - {\text{8}}} } \\ { - {\text{1}}{\text{.77}}} \\ \end{array} {\text{exp}}\left( {{{ - {\text{109}}{\text{.5}} \pm {\text{8}}{\text{.9}}} \mathord{\left/ {\vphantom {{ - {\text{109}}{\text{.5}} \pm {\text{8}}{\text{.9}}} {{\text{RT}}}}} \right. \kern-\nulldelimiterspace} {{\text{RT}}}}} \right) \hfill \\ \end{gathered} $$ where D is in m²s^?1 and activation energies are in kJ/mol. Diffusivities measured at 64 and 72 wt% SiO₂ are only slightly different from those at 68 wt% SiO₂ and frequently all measurements are within error of each other. Silicon, aluminum, iron, magnesium, and calcium EBDCs were also calculated from diffusion profiles by error function inversion techniques assuming constant diffusivity. With one exception, silicon EBDCs calculated by error function techniques are within error of Boltzmann-Matano EBDCs. Average diffusivities of Fe, Mg, and Ca were within a factor of 2.5 of silicon diffusivities whereas Al diffusivities were approximately half those of silicon. Alkalies diffused much more rapidly than silicon and non-alkalies, however their diffusivities were not quantitatively determined. Low activation energies for silicon EBDCs result in rapid diffusion at magmatic temperatures. Assuming that water and fluorine exert similar effects on melt viscosity at high temperatures, the viscosity can be calculated and used in the Eyring equation used to determine diffusivities, typically to within a factor of three of those measured in this study. This correlation between viscosity and diffusivity can be inverted to calculate viscosities of fluorine- and water-bearing granitic melts at magmatic temperatures; these viscosities are orders of magnitude below those of hydrous granitic melts and result in more rapid and effective separation of granitic magmas from partially molten source rocks. Comparison of Arrhenius parameters for diffusion measured in this study with Arrhenius parameters determined for diffusion in similar compositions at the same pressure demonstrates simple relationships between Arrhenius parameters, activation energy-E_a, kJ/mol, pre-exponential factor-D_o, m²s^?1, and the volatile, X=F or OH^?, to oxygen, O, ratio of the melt {(X/X+O)}: $$\begin{gathered} {\text{E}}a = - {\text{1533\{ }}{{\text{X}} \mathord{\left/ {\vphantom {{\text{X}} {\left( {{\text{X}} + {\text{O}}} \right)}}} \right. \kern-\nulldelimiterspace} {\left( {{\text{X}} + {\text{O}}} \right)}}{\text{\} }} + {\text{213}}{\text{.3}} \hfill \\ {\text{D}}_{\text{O}} = {\text{2}}{\text{.13}} \times {\text{10}}^{ - {\text{6}}} {\text{exp}}\left[ { - {\text{6}}{\text{.5\{ }}{{\text{X}} \mathord{\left/ {\vphantom {{\text{X}} {\left( {{\text{X}} + {\text{O}}} \right)}}} \right. \kern-\nulldelimiterspace} {\left( {{\text{X}} + {\text{O}}} \right)}}{\text{\} }}} \right] \hfill \\ \end{gathered} $$ These relationships can be used to estimate diffusion in various melts of dacitic to rhyolitic composition containing both fluorine and water. Calculations for the contamination of rhyolitic melts by dacitic enclaves at 800°C and 700°C provide evidence for the virtual inevitability of diffusive contamination in hydrous and fluorine-bearing magmas if they undergo magma mixing of any form.     

Reversed experimental calibration of the garnet-clinopyroxene Fe — Mg exchange thermometer

The temperature-sensitive Fe,Mg exchange equilibrium,