Lithologic controls on mineralization at the Lagunas Norte high-sulfidation epithermal gold deposit, northern Peru

The 13.1-Moz high-sulfidation epithermal gold deposit of Lagunas Norte, Alto Chicama District, northern Peru, is hosted in weakly metamorphosed quartzites of the Upper Jurassic to Lower Cretaceous Chimú Formation and in overlying Miocene volcanic rocks of dacitic to rhyolitic composition. The Dafne and Josefa diatremes crosscut the quartzites and are interpreted to be sources of the pyroclastic volcanic rocks. Hydrothermal activity was centered on the diatremes and four hydrothermal stages have been defined, three of which introduced Au ± Ag mineralization. The first hydrothermal stage is restricted to the quartzites of the Chimú Formation and is characterized by silice parda, a tan-colored aggregate of quartz-auriferous pyrite–rutile ± digenite infilling fractures and faults, partially replacing silty beds and forming cement of small hydraulic breccia bodies. The δ³⁴S values for pyrite (1.7–2.2?‰) and digenite (2.1?‰) indicate a magmatic source for the sulfur. The second hydrothermal stage resulted in the emplacement of diatremes and the related volcanic rocks. The Dafne diatreme features a relatively impermeable core dominated by milled slate from the Chicama Formation, whereas the Josefa diatreme only contains Chimú Formation quartzite clasts. The third hydrothermal stage introduced the bulk of the mineralization and affected the volcanic rocks, the diatremes, and the Chimú Formation. In the volcanic rocks, classic high-sulfidation epithermal alteration zonation exhibiting vuggy quartz surrounded by a quartz–alunite and a quartz–alunite–kaolinite zone is observed. Company data suggest that gold is present in solid solution or micro inclusions in pyrite. In the quartzite, the alteration is subtle and is manifested by the presence of pyrophyllite or kaolinite in the silty beds, the former resulting from relatively high silica activities in the fluid. In the quartzite, gold mineralization is hosted in a fracture network filled with coarse alunite, auriferous pyrite, and enargite. Alteration and mineralization in the breccias were controlled by permeability, which depends on the type and composition of the matrix, cement, and clast abundance. Coarse alunite from the main mineralization stage in textural equilibrium with pyrite and enargite has δ³⁴S values of 24.8–29.4?‰ and $ {\delta^{18 }}{{\mathrm{O}}_{{\mathrm{S}{{\mathrm{O}}_4}}}} $ values of 6.8–13.9?‰, consistent with H₂S as the dominant sulfur species in the mostly magmatic fluid and constraining the fluid composition to low pH (0–2) and logfO₂ of ?28 to ?30. Alunite–pyrite sulfur isotope thermometry records temperatures of 190–260 °C; the highest temperatures corresponding to samples from near the diatremes. Alunite of the third hydrothermal stage has been dated by ⁴⁰Ar/³⁹Ar at 17.0?±?0.22 Ma. The fourth hydrothermal stage introduced only modest amounts of gold and is characterized by the presence of massive alunite–pyrite in fractures, whereas barite, drusy quartz, and native sulfur were deposited in the volcanic rocks. The $ {\delta^{18 }}{{\mathrm{O}}_{{\mathrm{S}{{\mathrm{O}}_4}}}} $ values of stage IV alunite vary between 11.5 and 11.7?‰ and indicate that the fluid was magmatic, an interpretation also supported by the isotopic composition of barite (δ³⁴S?=?27.1 to 33.8?‰ and $ {\delta^{18 }}{{\mathrm{O}}_{{\mathrm{S}{{\mathrm{O}}_4}}}} $ ?=?8.1 to 12.7?‰). The Δ³⁴S_py–alu isotope thermometry records temperatures of 210 to 280 °C with the highest values concentrated around the Josefa diatreme. The Lagunas Norte deposit was oxidized to a depth of about 80 m below the current surface making exploitation by heap leach methods viable.     

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.     

Hydrochemical and isotopic characteristics of groundwater in the Ndop plain,northwest Cameroon: resilience to seasonal climatic changes

Shallow groundwater (>30 mbgl) is an essential source of drinking water to rural communities in the Ndop plain, northwest Cameroon. As a contribution to water management, the effect of seasonal variation on the groundwater chemistry, hydrochemical controls, drinking quality and recharge were investigated during the peaks of the dry (January) and rainy (September) seasons. Field measurements of physical parameters were preceded by sampling 58 groundwater samples during both seasons for major ions and stable isotope analyses. The groundwater, which was barely acidic (mean pH of 6) and less mineralised (TDS < 272 mg/l), showed no significant seasonal variation in temperature, pH and TDS during the two seasons. The order of cation abundance (meq/l) was Na⁺ > Ca²⁺ > Mg²⁺ > K⁺ and Na⁺ > Mg²⁺ > Ca²⁺ > K⁺ in the dry and rainy seasons, respectively, but that of anions ( \( {\text{HCO}}_{3}^{ - } \)  >  \( {\text{NO}}_{3}^{ - } \)  > Cl^? >  \( {\text{SO}}_{4}^{2 - } \)  > F^?) was similar in both seasons. This suggests a negligible effect of seasonal variations on groundwater chemistry. The groundwater, which was CaMgHCO₃ and NaHCO₃, is chemically evolved rainfall (CaMgSO₄Cl) in the area. Silicate mineral dissolution and cation-exchange were the main controls on groundwater chemistry while there was little anthropogenic influence. The major ions and TDS concentrations classified the water as suitable for human consumption as per WHO guidelines. The narrow cluster of δ¹⁸O and δD of same groundwater from both seasons between the δ¹⁸O and δD values of May–June precipitation along the Ndop Meteoric Water Line indicates meteoric origin, rapid recharge (after precipitation) and timing of recharge between May and June rainfall. Diffuse groundwater recharge mainly occurs at low altitudes (<1,400 m asl) within the plain. Besides major ions and TDS, the similar δ¹⁸O and δD of groundwater from both seasons indicate a consistent groundwater recharge and flow pattern throughout the year and resilience to present day short-term seasonal climatic variations. However, controlled groundwater abstraction is recommended given the increasing demand.     

Hydrogeochemical and isotopic characteristics of the Ilica geothermal system (Erzurum,Turkey)

In this paper, the hydrochemical isotopic characteristics of samples collected from geothermal springs in the Ilica geothermal field, Eastern Anatolia of Turkey, are examined and described. Low-temperature geothermal system of Ilica (Erzurum, Turkey) located along the Eastern Anatolian fault zone was investigated for hydrogeochemical and isotopic characteristics. The study of ionic and isotopic contents shows that the thermal water of Ilica is mainly, locally fed by groundwater, which changes chemically and isotopically during its circulation within the major fault zone reaching depths. The thermal spring has a temperature of 29–39 °C, with electrical conductivity ranging from 4,000 to 7,510 µS/cm and the thermal water is of Na–HCO₃–Cl water type. The chemical geothermometers applied in the Ilica geothermal waters yielded a maximum reservoir temperature of 142 °C according to the silica geothermometers. The thermal waters are undersaturated with respect to gypsum, anhydrite and halite, and oversaturated with respect to dolomite. The dolomite mineral possibly caused scaling when obtaining the thermal waters in the study area. According to the enthalpy chloride-mixing model, cold water to the thermal water-mixing ratio is changing between 69.8 and 75 %. The δ¹⁸O–δ²H compositions obviously indicate meteoric origin of the waters. Thermal water springs derived from continental precipitation falling on to higher elevations in the study area. The δ¹³C ratio for dissolved inorganic carbonate in the waters lies between 4.63 and 6.48 ‰. In low-temperature waters carbon is considered as originating from volcanic (mantle) CO₂.     

The crystal structure of galgenbergite-(Ce), CaCe2(CO3)4∙H2O

Galgenbergite-(Ce) from the type locality, the railroad tunnel Galgenberg between Leoben and St. Michael, Styria, Austria, was investigated. There it occurs in small fissures of an albite-chlorite schist as very thin tabular crystals building rosette-shaped aggregates associated with siderite, ancylite-(Ce), pyrite and calcite. Electron microprobe analyses gave CaO 9.49, Ce₂O₃ 28.95, La₂O₃ 11.70, Nd₂O₃ 11.86, Pr₂O₃ 3.48, CO₂ 30.00, H₂O 3.07, total 98.55 wt.%. CO₂ and H₂O calculated by stoichiometry. The empirical formula (based on Ca + REE ∑3.0) is $ \mathrm{C}{{\mathrm{a}}_{1.00 }}{{\left( {\mathrm{C}{{\mathrm{e}}_{1.04 }}\mathrm{L}{{\mathrm{a}}_{0.42 }}\mathrm{N}{{\mathrm{d}}_{0.42 }}\mathrm{P}{{\mathrm{r}}_{0.12 }}} \right)}_{2.00 }}{{\left( {\mathrm{C}{{\mathrm{O}}_3}} \right)}_4}\cdot {{\mathrm{H}}_2}\mathrm{O} $ , and the simplified formula is $ \mathrm{CaC}{{\mathrm{e}}_2}{{\left( {\mathrm{C}{{\mathrm{O}}_3}} \right)}_4}\cdot {{\mathrm{H}}_2}\mathrm{O} $ . According to X-ray single crystal diffraction galgenbergite-(Ce) is triclinic, space group $ P\overline{1},a=6.3916(5) $ , b?=?6.4005(4), c?=?12.3898(9) Å, α?=?100.884(4), β?=?96.525(4), γ?=?100.492(4)°, V?=?483.64(6) ų, Z?=?2. The eight strongest lines in the powder X-ray diffraction pattern are [d _calc in Å/(I)/hkl]: 5.052/(100)/011; 3.011/(70)/0-22; 3.006/(66)/004; 5.899/(59)/-101; 3.900/(51)/1-12; 3.125/(46)/-201; 2.526/(42)/022; 4.694/(38)/-102. The infrared absorption spectrum reveals H₂O (OH-stretching mode at 3,489 cm^?1, HOH bending mode at 1,607 cm^?1) and indicates the presence of distinctly non-equivalent CO₃-groups by double and quadruple peaks of their ν1, ν2, ν3 and ν4 modes. The crystal structure of galgenbergite-(Ce) was refined with X-ray single crystal data to R1?=?0.019 for 2,448 unique reflections (I?>?2σ(I)) and 193 parameters. The three cation sites of the structure Ca(1), Ce(2) and Ce(3) have a modest mixed site occupation by Ca and small amount of REE (Ce, La, Pr, Nd) and vice versa. The structure is based on double layers parallel to (001), which are composed of Ca(1)Ce(2)(CO₃)₂ single layers with an ordered chessboard like arrangement of Ca and Ce, and with a roof tile-like stacking of the CO₃ groups. Perpendicular to (001) the double layers are connected to a triclinic framework structure with good cleavage parallel to (001) by a differently organized and more open part of the structure formed by Ce(3)(CO₃)₂(H₂O). Based on the topology of the CaCe(CO₃)₂ single layer in galgenbergite-(Ce), structural relationships to rutherfordine, to aragonite and ancylite type minerals, and to lanthanite are outlined.     

Thermodynamic mixing properties of Mg(Al,Cr)2O4 spinel crystalline solution at high temperatures and pressures

Three Al-Cr exchange isotherms at 1,250°, 1,050°, and 796° between Mg(Al, Cr)₂O₄ spinel and (Al, Cr)₂O₃ corundum crystalline solutions have been studied experimentally at 25 kbar pressure. Starting from gels of suitable bulk compositions, close approach to equilibrium has been demonstrated in each case by time studies. Using the equation of state for (Al, Cr)₂O₃ crystalline solution (Chatterjee et al. 1982a) and assuming that the Mg(Al, Cr)₂O₄ can be treated in terms of the asymmetric Margules relation, the exchange isotherms were solved for Δ G ^*, and . The best constrained data set from the 1,250° C isotherm clearly shows that the latter two quantities do not overlap within three standard deviations, justifying the choice of asymmetric Margules relation for describing the excess mixing properties of Mg(Al, Cr)₂O₄ spinels. Based on these experiments, the following polybaric-polythermal equation of state can be formulated: , P expressed in bars, T in K, G _m ^ex and W _G,i ^Sp in joules/mol. Temperature-dependence of G _m ^ex is best constrained in the range 796–1,250° C; extrapolation beyond that range would have to be done with caution. Such extrapolation to lower temperature shows tentatively that at 1 bar pressure the critical temperature, T _c, of the spinel solvus is 427° C, with dT_c/dP≈1.3 K/kbar. The critical composition, X _c, is 0.42 , and changes barely with pressure. Substantial error in calculated phase diagrams will result if the significant positive deviation from ideality is ignored for Al-Cr mixing in such spinels.     

Weathering Processes in the Min Jiang: Major Elements, {}^{87}{\text{Sr/}}{}^{86}{\text{Sr}},\;\delta {}^{34}{\text{S}}_{{\text{SO}}_{\text{4}} } ,\;{\text{and}}\;\delta {}^{18}{\text{O}}_{{\text{SO}}_{\text{4}} }

We investigated the dissolved major elements, $ {}^{87}{\text{Sr/}}{}^{86}{\text{Sr}},\;\delta {}^{34}{\text{S}}_{{\text{SO}}_{\text{4}} } ,\;{\text{and}}\;\delta {}^{18}{\text{O}}_{{\text{SO}}_{\text{4}} } $ composition of the Min Jiang, a headwater tributary of the Chang Jiang (Yangtze River). A forward calculation method was applied to quantify the relative contribution to the dissolved load from rain, evaporite, carbonate, and silicate reservoirs. Input from carbonate weathering dominated the major element composition (58–93%) and that from silicate weathering ranged from 2 to 18% in unperturbed Min Jiang watersheds. Most samples were supersaturated with respect to calcite, and the CO₂ partial pressures were similar to or up to ～5 times higher than atmospheric levels. The Sr concentrations in our samples were low (1.3–2.5 μM) with isotopic composition ranging from 0.7108 to 0.7127, suggesting some contribution from felsic silicates. The Si/(Na^* + K) ratios ranged from 0.5 to 2.5, which indicate low to moderate silicate weathering intensity. The $ \delta {}^{34}{\text{S}}_{{\text{SO}}_{\text{4}} } \;{\text{and}}\;\delta {}^{18}{\text{O}}_{{\text{SO}}_{\text{4}} } $ for five select samples showed that the source of dissolved sulfate was combustion of locally consumed coal. The silicate weathering rates were 23–181 × 10³ mol/km²/year, and the CO₂ consumption rates were 31–246 × 10³ mol/km²/year, which are moderate on a global basis. Upon testing various climatic and geomorphic factors for correlation with the CO₂ consumption rate, the best correlation coefficients found were with water temperature (r ² = 0.284, p = 0.009), water discharge (r ² = 0.253, p = 0.014), and relief (r ² = 0.230, p = 0.019).     

A hydrochemical study of the Hammam Righa geothermal waters in north-central Algeria

This study focuses on the hydrochemical characteristics of 47 water samples collected from thermal and cold springs that emerge from the Hammam Righa geothermal field, located in north-central Algeria. The aquifer that feeds these springs is mainly situated in the deeply fractured Jurassic limestone and dolomite of the Zaccar Mount. Measured discharge temperatures of the cold waters range from 16.0 to 26.5 °C and the hot waters from 32.1 to 68.2 °C. All waters exhibited a near-neutral pH of 6.0–7.6. The thermal waters had a high total dissolved solids (TDS) content of up to 2527 mg/l, while the TDS for cold waters was 659.0–852.0 mg/l. Chemical analyses suggest that two main types of water exist: hot waters in the upflow area of the Ca–Na–SO₄ type (Hammam Righa) and cold waters in the recharge zone of the Ca–Na–HCO₃ type (Zaccar Mount). Reservoir temperatures were estimated using silica geothermometers and fluid/mineral equilibria at 78, 92, and 95 °C for HR4, HR2, and HR1, respectively. Stable isotopic analyses of the δ¹⁸O and δD composition of the waters suggest that the thermal waters of Hammam Righa are of meteoric origin. We conclude that meteoric recharge infiltrates through the fractured dolomitic limestones of the Zaccar Mount and is conductively heated at a depth of 2.1–2.2 km. The hot waters then interact at depth with Triassic evaporites located in the hydrothermal conduit (fault), giving rise to the Ca–Na–SO₄ water type. As they ascend to the surface, the thermal waters mix with shallower Mg-rich groundwater, resulting in waters that plot in the immature water field in the Na–K–Mg diagram. The mixing trend between cold groundwaters from the recharge zone area (Zaccar Mount) and hot waters in the upflow area (Hammam Righa) is apparent via a chloride-enthalpy diagram that shows a mixing ratio of 22.6 < R < 29.2 %. We summarize these results with a geothermal conceptual model of the Hammam Righa geothermal field.     

Geochemical and isotopic characteristics of fluids in the Niutuozhen geothermal field,North China

Niutuozhen geothermal field is located in the Jizhong graben, belonging to the northern part of Bohai Bay Basin in North China. Chemical and isotopic analyses were carried out on 14 samples of the geothermal fluids discharged from Neogene Minghuazhen (Nm), Guantao (Ng), and Jixianian Wumishan (Jxw) formations. The δ²H and δ¹⁸O in water, δ¹³C in CH₄, δ¹³C in CO₂, and ³He/⁴He ratio in the gases were analyzed in combination with chemical analyses on the fluids in the Niutuozhen geothermal field. The chemical and isotopic compositions indicate a meteoric origin of the thermal waters. The reservoir temperatures estimated by chemical geothermometry are in the range between 60 and 108 °C. The results show that the gases are made up mainly by N₂ (18.20–97.42 vol%), CH₄ (0.02–60.95 vol%), and CO₂ (0.17–25.14 vol%), with relatively high He composition (up to 0.52 vol%). The chemical and isotopic compositions of the gas samples suggest the meteoric origin of N₂, predominant crustal origins of CH₄, CO₂, and He. The mantle-derived He contributions are calculated to be from 5 to 8% based on a crust–mantle binary mixing model. The deep temperatures in the Jxw reservoir were evaluated based on gas isotope geothermometry to be in the range from 141 to 165 °C. The mantle-derived heat fraction in the surface heat flow is estimated to be in the range of 48–51% based on ³He/⁴He ratios.     

Isotope–geochemical characterization and geothermometrical modeling of Uttarakhand geothermal field,India

Uttarakhand geothermal area, located in the central belt of the Himalayan geothermal province, is one of the important high temperature geothermal fields in India. In this study, the chemical characteristics of the thermal waters are investigated to identify the main geochemical processes affecting the composition of thermal waters during its ascent toward the surface as well as to determine the subsurface temperature of the feeding reservoir. The thermal waters are mainly Ca–Mg–HCO₃ type with moderate silica and TDS concentrations. Mineral saturation states calculated from PHREEQC geochemical code indicate that thermal waters are supersaturated with respect to calcite, dolomite, aragonite, chalcedony, quartz (SI > 0), and undersaturated with respect to gypsum, anhydrite, and amorphous silica (SI < 0). XRD study of the spring deposit samples fairly corroborates the predicted mineral saturation state of the thermal waters. Stable isotopes (δ¹⁸O, δ²H) data confirm the meteoric origin of the thermal waters with no oxygen-18 shift. The mixing phenomenon between thermal water with shallow ground water is substantiated using tritium (³H) and chemical data. The extent of dilution is quantified using tritium content of thermal springs and non-thermal waters. Classical geothermometers, mixing model, and multicomponent fluid geothermometry modeling (GeoT) have been applied to estimate the subsurface reservoir temperature. Among different classical geothermometers, only quartz geothermometer provide somewhat reliable estimation (96–140 °C) of the reservoir temperature. GeoT modeling results suggest that thermal waters have attained simultaneous equilibrium with respect to minerals like calcite, quartz, chalcedony, brucite, tridymite, cristobalite, talc, at the temperature 130 ± 5 °C which is in good agreement with the result obtained from the mixing model.     

Later stages of evolution of an epithermal system: Au–Ag mineralizations at Apigania Bay, Tinos Island, Cyclades, Hellas, Greece

Precious metals accompany all types of epithermal deposits. In general, the largest of these deposits occur in intrusive or extrusive rocks of alkaline or calc-alkaline affinity. The Apigania Bay vein system and Au–Ag mineralization is hosted in Mesozoic marbles and schists, and is composed primarily of five nearly parallel, high-angle quartz veins that extend for at least 200 m. Gold–silver mineralization, in association with more than thirty ore and vein minerals, is developed in three stages and occurs at the contact of marbles and schists. Zones of epidote–chlorite–calcite and sericite–albite alteration are associated with precious metal-bearing milky and clear quartz veins. Fluid inclusion studies suggest that hydrothermal mineralization was deposited under hydrostatic pressures of ~100 bars, at temperature of 120–235°C, from low to moderate, calcium-bearing, saline fluids of 0.2 to 6.8 equiv. wt.% NaCl. Calculated isotope compositions (δ¹⁸O?=??4.7‰ to 1.7‰ and δD?=??120‰ to ?80‰) for waters in equilibrium with milky and clear quartz are consistent with mixing with dilute, low temperature meteoric ore fluids. Calculated δ ¹³C_CO2 (0.6‰ to 1.1‰) and δ ³⁴S_H2S (?7.3 to ?0.3‰) compositions of the ore fluids indicate exchange, in an open system, with a metasedimentary source. Gold and silver deposition was associated with degassing of hydrogen due to intense uplift of the mineralizing area. The physicochemical conditions of mineralization stages I to III range between 200°C and 150°C, $f_{{\text{S}}_2 } = 10^{ - 18.1} $ to 10^?16.8, $f_{{\text{O}}_2 } = 10^{ - 44.0} $ to 10^?41.5, pH?=?6.9 to7.6, $f_{{\text{H}}_{\text{2}} {\text{S}}} = 10^{ - 3.4} $ to 10^?2.6 and $a_{{\text{H}}_{\text{2}} {\text{S}}} = 10^{ - 2.7} $ to 10^?2.6. Apigania Bay could be possibly considered the latest evolutional phase of Tinos hydrothermal system.     

Experimental determination of the reaction paragonite=jadeite+kyanite+H2O,and internally consistent thermodynamic data for part of the system Na2O−Al2O3−SiO2−H2O,with applications to eclogites and blueschists

Hydrothermal reversal experiments have been performed on the upper pressure stability of paragonite in the temperature range 550–740 ° C. The reaction $$\begin{gathered} {\text{NaAl}}_{\text{3}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{1 0}}} ({\text{OH)}}_{\text{2}} \hfill \\ {\text{ paragonite}} \hfill \\ {\text{ = NaAlSi}}_{\text{2}} {\text{O}}_{\text{6}} + {\text{Al}}_{\text{2}} {\text{SiO}}_{\text{5}} + {\text{H}}_{\text{2}} {\text{O}} \hfill \\ {\text{ jadeite kyanite vapour}} \hfill \\ \end{gathered}$$ has been bracketed at 550 ° C, 600 ° C, 650 ° C, and 700 ° C, at pressures 24–26 kb, 24–25.5 kb, 24–25 kb, and 23–24.5 kb respectively. The reaction has a shallow negative slope (? 10 bar °C^?1) and is of geobarometric significance to the stability of the eclogite assemblage, omphacite+kyanite. The experimental brackets are thermodynamically consistent with the lower pressure reversals of Chatterjee (1970, 1972), and a set of thermodynamic data is presented which satisfies all the reversal brackets for six reactions in the system Na₂O-Al₂O₃-SiO₂-H₂O. The Modified Redlich Kwong equation for H₂O (Holloway, 1977) predicts fugacities which are too high to satisfy the reversals of this study. The P-T stabilities of important eclogite and blueschist assemblages involving omphacite, kyanite, lawsonite, Jadeite, albite, chloritoid, and almandine with paragonite have been calculated using thermodynamic data derived from this study.     

Semi-empirical model for retrieval of soil moisture using RISAT-1 C-Band SAR data over a sub-tropical semi-arid area of Rewari district,Haryana (India)

We have estimated soil moisture (SM) by using circular horizontal polarization backscattering coefficient (\(\sigma ^{\mathrm{o}}_{\mathrm{RH}}\)), differences of circular vertical and horizontal \(\sigma ^{\mathrm{o}} \, (\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}})\) from FRS-1 data of Radar Imaging Satellite (RISAT-1) and surface roughness in terms of RMS height (\({\hbox {RMS}}_{\mathrm{height}}\)). We examined the performance of FRS-1 in retrieving SM under wheat crop at tillering stage. Results revealed that it is possible to develop a good semi-empirical model (SEM) to estimate SM of the upper soil layer using RISAT-1 SAR data rather than using existing empirical model based on only single parameter, i.e., \(\sigma ^{\mathrm{o}}\). Near surface SM measurements were related to \(\sigma ^{\mathrm{o}}_{\mathrm{RH}}\), \(\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}}\) derived using 5.35 GHz (C-band) image of RISAT-1 and \({\hbox {RMS}}_{\mathrm{height}}\). The roughness component derived in terms of \({\hbox {RMS}}_{\mathrm{height}}\) showed a good positive correlation with \(\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}} \, (R^{2} = 0.65)\). By considering all the major influencing factors (\(\sigma ^{\mathrm{o}}_{\mathrm{RH}}\), \(\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}}\), and \({\hbox {RMS}}_{\mathrm{height}}\)), an SEM was developed where SM (volumetric) predicted values depend on \(\sigma ^{\mathrm{o}}_{\mathrm{RH}}\), \(\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}}\), and \({\hbox {RMS}}_{\mathrm{height}}\). This SEM showed \(R^{2}\) of 0.87 and adjusted \(R^{2}\) of 0.85, multiple R=0.94 and with standard error of 0.05 at 95% confidence level. Validation of the SM derived from semi-empirical model with observed measurement (\({\hbox {SM}}_{\mathrm{Observed}}\)) showed root mean square error (RMSE) = 0.06, relative-RMSE (R-RMSE) = 0.18, mean absolute error (MAE) = 0.04, normalized RMSE (NRMSE) = 0.17, Nash–Sutcliffe efficiency (NSE) = 0.91 (\({\approx } 1\)), index of agreement (d) = 1, coefficient of determination \((R^{2}) = 0.87\), mean bias error (MBE) = 0.04, standard error of estimate (SEE) = 0.10, volume error (VE) = 0.15, variance of the distribution of differences \(({\hbox {S}}_{\mathrm{d}}^{2}) = 0.004\). The developed SEM showed better performance in estimating SM than Topp empirical model which is based only on \(\sigma ^{\mathrm{o}}\). By using the developed SEM, top soil SM can be estimated with low mean absolute percent error (MAPE) = 1.39 and can be used for operational applications.     

Ephesite,Na(LiAl2) [Al2Si2O10] (OH)2: I. thermal stability and standard state thermodynamic properties

Ephesite, Na(LiAl₂) [Al₂Si₂O₁₀] (OH)₂, has been synthesized for the first time by hydrothermal treatment of a gel of requisite composition at 300≦T(° C)≦700 and \(P_{H_2 O}\) upto 35 kbar. At \(P_{H_2 O}\) between 7 and 35 kbar and above 500° C, only the 2M₁ polytype is obtained. At lower temperatures and pressures, the 1M polytype crystallizes first, which then inverts to the 2M₁ polytype with increasing run duration. The X-ray diffraction patterns of the 1M and 2M₁ poly types can be indexed unambiguously on the basis of the space groups C2 and Cc, respectively. At its upper thermal stability limit, 2M₁ ephesite decomposes according to the reaction (1) $$\begin{gathered} {\text{Na(LiAl}}_{\text{2}} {\text{) [Al}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{{\text{10}}} {\text{] (OH)}}_{\text{2}} \hfill \\ {\text{ephesite}} \hfill \\ {\text{ = Na[AlSiO}}_{\text{4}} {\text{] + LiAl[SiO}}_{\text{4}} {\text{] + }}\alpha {\text{ - Al}}_{\text{2}} {\text{O}}_{\text{3}} {\text{ + H}}_{\text{2}} {\text{O}} \hfill \\ {\text{nepheline }}\alpha {\text{ - eucryptite corundum}} \hfill \\ \end{gathered}$$ Five reversal brackets for (1) have been established experimentally in the temperature range 590–750° C, at \(P_{H_2 O}\) between 400 and 2500 bars. The equilibrium constant, K, for this reaction may be expressed as (2) $$log K{\text{ = }}log f_{{\text{H}}_{\text{2}} O}^* = 7.5217 - 4388/T + 0.0234 (P - 1)T$$ where \(f_{H_2 O}^* = f_{H_2 O} (P,T)/f_{H_2 O}^0\) (1,T), with T given in degrees K, and P in bars. Combining these experimental data with known thermodynamic properties of the decomposition products in (1), the following standard state (1 bar, 298.15 K) thermodynamic data for ephesite were calculated: H _f,298.15 ⁰ =-6237372 J/mol, S _298.15 ⁰ =300.455 J/K·mol, G _298.15 ⁰ =-5851994 J/mol, and V _298.15 ⁰ =13.1468 J/bar·mol.     

Some geobarometers involving cordierite in the FeO-Al2O3-SiO2(±H2O) system: refinements,thermodynamic calibration,and applicability in granulite facies rocks

Five geobarometers involving cordierite have been formulated for quantitative pressure sensing in high grade metapelites. The relevant reactions in the FeO-Al₂O₃-SiO₂ (±H₂O) system are based on the assemblages (A) cordierite-garnet-sillimanite-quartz, (B) cordierite-spinel-quartz, (C) cordierite-garnet-spinel-sillimanite, (D) cordierite-garnet-orthopyroxene-quartz and (E) cordierite-orthopyroxene-sillimanite-quartz. Application of the barometric formulations to a large number of granulite grade rocks indicates that the cordierite-garnet-sillimanite-quartz equilibrium is widely applicable and registers pressures which are in good agreement with the “consensus” pressure estimates. The dispersion in the computed P values, expressed as one standard deviation, is within ±1.2 kbar. The geobarometers (B) and (C) also yield pressures which are reasonable and compare well with those computed from equilibrium (A). The estimated pressures from (D) and (E), both involving orthopyroxene, are at variance with these estimates. It has been argued that the discrepancy in pressures obtained from these geobarometers stems from an inadequate knowledge of activity-composition relations and/or errors in input thermodynamic data of aluminous orthopyroxene. The convergence of pressure values estimated from the barometric formulations, especially (A), (B) and (C), implies that the present formulations are more dependable than the existing formulations and are also capable of setting limits on P values in response to varying $$\begin{gathered} {\text{1/2Fe}}_{\text{2}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} \hfill \\ {\text{ = 1/3Fe}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} {\text{ + 2/3Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + 5/6SiO}}_{\text{2}} {\text{. (A)}} \hfill \\ {\text{1/2Fe}}_{\text{2}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} {\text{ = FeAl}}_{\text{2}} {\text{O}}_{\text{4}} {\text{ + 5/2SiO}}_{\text{2}} {\text{. (B)}} \hfill \\ {\text{Fe}}_{\text{2}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} {\text{ + FeAl}}_{\text{2}} {\text{O}}_{\text{4}} \hfill \\ = {\text{Fe}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} {\text{ + 2Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{. (C)}} \hfill \\ {\text{1/2Fe}}_{\text{2}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} {\text{ + Fe}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} \hfill \\ = {\text{Fe}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} {\text{ + 3/2SiO}}_{\text{2}} .{\text{ (D)}} \hfill \\ {\text{1/2Fe}}_{\text{2}} {\text{Al}}{}_{\text{4}}{\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} \hfill \\ = 1/2{\text{Fe}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} {\text{ + Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + 1/2SiO}}_{\text{2}} .{\text{ (E)}} \hfill \\ \end{gathered}$$ . The present communication addresses the calibration, applicability and reliability of these barometers with reference to granulite facies metapelites.     

H-bonding scheme and cation partitioning in axinite: a single-crystal neutron diffraction and Mössbauer spectroscopic study

The crystal chemistry of a ferroaxinite from Colebrook Hill, Rosebery district, Tasmania, Australia, was investigated by electron microprobe analysis in wavelength-dispersive mode, inductively coupled plasma–atomic emission spectroscopy (ICP–AES), ⁵⁷Fe Mössbauer spectroscopy and single-crystal neutron diffraction at 293 K. The chemical formula obtained on the basis of the ICP–AES data is the following: \( ^{X1,X2} {\text{Ca}}_{4.03} \,^{Y} \left( {{\text{Mn}}_{0.42} {\text{Mg}}_{0.23} {\text{Fe}}^{2 + }_{1.39} } \right)_{\varSigma 2.04} \,^{Z1,Z2} \left( {{\text{Fe}}^{3 + }_{0.15} {\text{Al}}_{3.55} {\text{Ti}}_{0.12} } \right)_{\varSigma 3.82} \,^{T1,T2,T3,T4} \left( {{\text{Ti}}_{0.03} {\text{Si}}_{7.97} } \right)_{\varSigma 8} \,^{T5} {\text{B}}_{1.96} {\text{O}}_{30} \left( {\text{OH}} \right)_{2.18} \). The ⁵⁷Fe Mössbauer spectrum shows unambiguously the occurrence of Fe²⁺ and Fe³⁺ in octahedral coordination only, with Fe²⁺/Fe³⁺ = 9:1. The neutron structure refinement provides a structure model in general agreement with the previous experimental findings: the tetrahedral T1, T2, T3 and T4 sites are fully occupied by Si, whereas the T5 site is fully occupied by B, with no evidence of Si at the T5, or Al or Fe³⁺ at the T1–T5 sites. The structural and chemical data of this study suggest that the amount of B in ferroaxinite is that expected from the ideal stoichiometry: 2 a.p.f.u. (for 32 O). The atomic distribution among the X1, X2, Y, Z1 and Z2 sites obtained by neutron structure refinement is in good agreement with that based on the ICP–AES data. For the first time, an unambiguous localization of the H site is obtained, which forms a hydroxyl group with the oxygen atom at the O16 site as donor. The H-bonding scheme in axinite structure is now fully described: the O16–H distance (corrected for riding motion effect) is 0.991(1) Å and an asymmetric bifurcated bonding configuration occurs, with O5 and O13 as acceptors [i.e. with O16···O5 = 3.096(1) Å, H···O5 = 2.450(1) Å and O16–H···O5 = 123.9(1)°; O16···O13 = 2.777(1) Å, H···O13 = 1.914(1) Å and O16–H···O13 = 146.9(1)°].     

Chemistry of Neutral and Alkaline Waters with Low Al<Superscript>3+</Superscript> Activity Against Hydroxyaluminosilicate HAS<Subscript>B</Subscript> Solubility. The Evidence from Ground and Surface Waters of the Sudetes Mts. (SW Poland)

A wide set of aqueous chemistry data (574 water analyses) from natural environments has been used to testify and validate of the solubility of synthetic hydroxyaluminosilicate (HAS_B)_, Al₂Si₂O₅(OH)₄. The ground and surface waters represent regolith and/or fissure aquifers in various (magmatic, sedimentary and metamorphic) bedrocks in the Sudetes Mts. (SW Poland). The solubility of HAS_B in natural waters was calculated using the method proposed by Schneider et al. (Polyhedron 23:3185–3191, 2004). Results confirm usefulness and validity of this method. The HAS_B solubility obtained from the field data (logK_sp = −44.7 ± 0.58) is lower than it was estimated (logK_sp = −40.6 ± 0.15) experimentally (Schneider et al. Polyhedron 23:3185–3191, 2004). In the waters studied the equilibrium with HAS_B is maintained at pH above 6.7 and at [Al³⁺] ≤ 10⁻¹⁰. Silicon activity (log[H₄SiO₄]) ranges between −4.2 and −3.4. Due to the calculation method used, the K_sp mentioned above cannot be considered as a classical solubility constant. However, it can be used in the interpretation of aluminium solubility in natural waters. The HAS_B has solubility lower than amorphous Al(OH)₃, and higher than proto-imogolite. From water samples that are in equilibrium with respect to HAS_B, the solubility product described by the reaction, is calculated to be logK_sp = 14.0 (±0.7) at 7°C.     

Sources and Processes Affecting Nitrate in a Dam-Controlled Subtropical River,Southwest China

Excess nutrient (N and P) loads are recognized as the major cause of serious water quality problems in China. River systems play a very important role in nitrate (NO₃ ^?) transportation and transformation in the aquatic environment. To understand and clarify the sources and processes affecting NO₃ ^? in river basins, we have examined spatial and temporal variations of concentration and dual-isotopic composition of NO₃ ^? in the dam-controlled Jialing River, a major tributary of the Yangtze River where land use is dominated by agriculture. Water samples were collected in July 2008 and February 2009 from the main channel of the Jialing River and its major tributaries. The δ¹⁵N and δ¹⁸O of NO₃ ^? range from 1.5 to 11.0 ‰ (average 6.2 ‰) and ?5.0 to 11.1 ‰ (average, 1.6 ‰), respectively. NO₃ ^? isotope data and δ¹⁸O of water interpreted in combination with hydrological and chemical data suggest that most of the NO₃ ^? input is from nitrification during the rainy season, and discharge of sewage and manure in the upper course and from cities accounts for much of the NO₃ ^? load during the dry season. The construction of cascade dams has led to retention of Si and a decrease in the Si/N ratio, implying that assimilation and/or denitrification may significantly affect NO₃ ^? in the dam area, as demonstrated by NO₃ ^? and dissolved Si concentrations, and \(\updelta^{ 1 5} {\text{N}}_{{{\text{NO}}_{3} }}\) and \(\updelta^{ 1 8} {\text{O}}_{{{\text{NO}}_{3} }}\) values. This study indicates that dual-isotopic data can be used to identify NO₃ ^? pollution sources and the processes NO₃ ^? has undergone during its retention and transport in the watershed of the dam-controlled Jialing River.     

Stabilitätsbeziehungen zwischen Chlorit,Cordierit und Almandin bei der Metamorphose

The stability relations between cordierite and almandite in rocks, having a composition of CaO poor argillaceous rocks, were experimentally investigated. The starting material consisted of a mixture of chlorite, muscovite, and quartz. Systems with widely varying Fe²⁺/Fe²⁺+Mg ratios were investigated by using two different chlorites, thuringite or ripidolite, in the starting mixture. Cordierite is formed according to the following reaction: $${\text{Chlorite + muscovite + quartz}} \rightleftharpoons {\text{cordierite + biotite + Al}}_{\text{2}} {\text{SiO}}_{\text{5}} + {\text{H}}_{\text{2}} {\text{O}}$$ . At low pressures this reaction characterizes the facies boundary between the albite-epidotehornfels facies and the hornblende-hornfels facies, at medium pressures the beginning of the cordierite-amphibolite facies. Experiments were carried out reversibly and gave the following equilibrium data: 505±10°C at 500 bars H₂O pressure, 513±10°C at 1000 bars H₂O pressure, 527±10°C at 2000 bars H₂O pressure, and 557±10°C at 4000 bars H₂O pressure. These equilibrium data are valid for the Fe-rich starting material, using thuringite as the chlorite, as well as for the Mg-rich starting mixture with ripidolite. At 6000 bars the equilibrium temperature for the Mg-rich mixture is 587±10°C. In the Fe-rich mixture almandite was formed instead of cordierite at 6000 bars. The following reaction was observed: $${\text{Thuringite + muscovite + quartz}} \rightleftharpoons {\text{almandite + biotite + Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + H}}_{\text{2}} {\text{O}}$$ . Experiments with the Fe-rich mixture, containing Fe²⁺/Fe²⁺+Mg in the ratio 8∶10, yielded three stability fields in a P,T-diagram (Fig.1):

1. Above 600°C/5.25 kb and 700°C/6.5 kb almandite+biotite+Al₂SiO₅ coexist stably, cordierite being unstable.
2. The field, in which almandite, biotite and Al₂SiO₅ are stable together with cordierite, is restricted by two curves, passing through the following points:
   1. 625°C/5.5 kb and 700°C/6.5 kb,
   2. 625°C/5.5 kb and 700°C/4.0 kb.
3. At conditions below curves 1 and 2b, cordierite, biotite, and Al₂SiO₅ are formed, but no garnet.

An appreciable MnO-content in the system lowers the pressures needed for the formation of almandite garnet, but the quantitative influence of the spessartite-component on the formation of almandite could not yet be determined. the Mg-rich system with Fe²⁺/Fe²⁺+Mg=0.4 garnet did not form at pressures up to 7 kb in the temperature range investigated. Experiments at unspecified higher pressures (in a simple squeezer-type apparatus) yielded the reaction: $${\text{Ripidolite + muscovite + quartz}} \rightleftharpoons {\text{almandite + biotite + Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + H}}_{\text{2}} {\text{O}}$$ . Further experiments are needed to determine the equilibrium data. The occurence of garnet in metamorphic rocks is discussed in the light of the experimental results.     

Water diffusion in Mount Changbai peralkaline rhyolitic melt

Diffusion couple experiments with wet half (up to 4.6 wt%) and dry half were carried out at 789–1,516 K and 0.47–1.42 GPa to investigate water diffusion in a peralkaline rhyolitic melt with major oxide concentrations matching Mount Changbai rhyolite. Combining data from this work and a related study, total water diffusivity in peralkaline rhyolitic melt can be expressed as:

$ D_{{{\text{H}}_{ 2} {\text{O}}_{\text{t}} }} = D_{{{\text{H}}_{ 2} {\text{O}}_{\text{m}} }} \left( {1 - \frac{0.5 - X}{{\sqrt {[4\exp (3110/T - 1.876) - 1](X - X^{2} ) + 0.25} }}} \right), $

$ {\text{with}}\;D_{{{\text{H}}_{ 2} {\text{O}}_{\text{m}} }} = \exp \left[ { - 1 2. 7 8 9- \frac{13939}{T} - 1229.6\frac{P}{T} + ( - 27.867 + \frac{60559}{T})X} \right], $

where D is in m² s^?1, T is the temperature in K, P is the pressure in GPa, and X is the mole fraction of water and calculated as X = (C/18.015)/(C/18.015 + (100 ? C)/33.14), where C is water content in wt%. We recommend this equation in modeling bubble growth and volcanic eruption dynamics in peralkaline rhyolitic eruptions, such as the ~1,000-ad eruption of Mount Changbai in North East China. Water diffusivities in peralkaline and metaluminous rhyolitic melts are comparable within a factor of 2, in contrast with the 1.0–2.6 orders of magnitude difference in viscosities. The decoupling of diffusivity of neutral molecular species from melt viscosity, i.e., the deviation from the inversely proportional relationship predicted by the Stokes–Einstein equation, might be attributed to the small size of H₂O molecules. With distinct viscosities but similar diffusivity, bubble growth controlled by diffusion in peralkaline and metaluminous rhyolitic melts follows similar parabolic curves. However, at low confining pressure or low water content, viscosity plays a larger role and bubble growth rate in peralkaline rhyolitic melt is much faster than that in metaluminous rhyolite.