Stable isotope characteristics of precipitation of Pamba River basin,Kerala, India

Stable isotope composition of precipitation from Pamba River basin, Kerala, India, is evaluated to understand the role of spatial and temporal variations on rainwater isotope characteristics. Physiographically different locations in the basin showed strong spatial and temporal variations. δ ¹⁸O varied from ?7.63 to ?1.75 ‰ in the lowlands; from ?9.32 to ?1.94 ‰ in the midlands and from ?11.6 to ?4.00 ‰ in the highlands. Local Meteoric Water Lines (LMWL) for the three regions were determined separately and an overall LMWL for the whole of the basin was found to be δ ²H = 6.6 (±0.4) δ ¹⁸ O+10.4 (±2.0). Altitude effect was evident for the basin (0.1 ‰ for δ ¹⁸O and 0.8 ‰ for δ ²H per 100 m elevation), while the amount effect was weak. The precipitation formed from the marine moisture supplied at a steady rate, without much isotopic evolution in this period may have masked the possible depletion of heavier isotopes with increasing rainfall. Consistently high d-excess values showed the influence of recycled vapour, despite the prevailing high relative humidity. The oceanic and continental vapour source origins for the south-west and north-east monsoons were clearly noted in the precipitation in the basin. Rayleigh distillation model showed about 30% rainout of the monsoon vapour mass in the basin.     

Petrogenesis of peralkaline granite dykes of the Straumsvola complex,western Dronning Maud Land,Antarctica

Peralkaline syenite and granite dykes cut the Straumsvola nepheline syenite pluton in Western Dronning Maud Land, Antarctica. The average peralkalinity index (PI?=?molecular Al/[Na?+?K]) of the dykes is 1.20 (n?=?29) and manifests itself in the presence of the Zr silicates eudialyte, dalyite and vlasovite, and the Na–Ti silicate, narsarsukite. The dykes appear to have intruded during slow cooling of the nepheline syenite pluton, and the petrogenetic relationship of the dykes and the pluton cannot be related to closed-system processes at low pressure, given the thermal divide that exists between silica-undersaturated and oversaturated magmas. Major and trace element variations in the dykes are consistent with a combination of fractional crystallization of parental peralkaline magma of quartz trachyte composition, and internal mineral segregation prior to final solidification. The distribution of accessory minerals is consistent with late-stage crystallization of isolated melt pockets. The dykes give an Rb–Sr isochron age of 171?±?4.4 Ma, with variable initial ⁸⁷Sr/⁸⁶Sr ratio (0.7075?±?0.0032), and have an average ε _Nd of ? 12.0. Quartz phenocrysts have δ¹⁸O values of 8.4–9.2‰, which are generally in O-isotope equilibrium with bulk rock. Differences in the δ¹⁸O values of quartz and aegirine (average Δ_{quartz?aegirine} = 3.5‰) suggest aegirine formation temperatures around 500 °C, lower than expected for a felsic magma, but consistent with poikilitic aegirine that indicates subsolidus growth. The negative ε _Nd (< ? 10) and magma δ¹⁸O values averaging 8.6‰ (assuming Δ_quartz?magma = 0.6‰) are inconsistent with a magma produced by closed-system fractional crystallization of a mantle-derived magma. By contrast, the nepheline syenite magma had mantle-like δ¹⁸O values and much less negative ε _Nd (average ??3.1, n?=?3). The country rock has similar δ¹⁸O values to the granite dykes (average 8.0‰, n?=?108); this means that models for the petrogenesis of the granites by assimilation are unfeasible, unless an unexposed high-δ¹⁸O contaminant is invoked. Instead, it is proposed that the peralkaline syenite and granite dykes formed by partial melting of alkali-metasomatised gneiss that surrounds the nepheline syenite, followed by fractional crystallization.     

Hydrochemical and isotopic investigation of a tropical wetland system in the Indian subcontinent

The study reported here is a part of an attempt to establish a comprehensive hydrochemical and isotopic baseline for a tropical wetland system as background data for a range of applications. Surface water samples of Vembanad Lake were collected from 20 stations in three seasons during the period 2007–2009. The analytical results were subjected to different chemical classification techniques to understand processes affecting the chemical concentration of waters. The Piper diagram classified the water samples as 100% alkali group in pre-monsoon followed by 15% in monsoon and 85% in post-monsoon, and for anions 100% samples were of strong acids followed by 90% in monsoon and 100% in post-monsoon season. The plot to decipher the mechanism controlling water chemistry placed the Vembanad Lake in the region of precipitation and rock dominance in the monsoon season and in the field of saline water dominance in pre-monsoon and post-monsoon season. The positive values for the chloro-alkaline indices in pre and post-monsoon season promoted cation exchange in the system. The stable isotopes of water samples ranged from ?20.21 to +17.0‰ and ?5.6 to +3.34‰ for δ ²H and δ ¹⁸O, respectively. The most depleted δ values observed in the monsoon are due to the amount effect. The high enrichment observed in pre-monsoon is primarily due to evaporation and salinity mixing. The variation of isotopes in the whole system point toward the fact that salinity mixing can be indicated by the δ ¹⁸O variation and δ ²H indicates the evaporation effect. The plot of δ ¹⁸O with chloride concentration showed precipitation dominance in the monsoon season, mixing of saline water and evaporation in pre-monsoon season, whereas the post-monsoon samples plot in both fresh and saline region.     

Nature of serpentinization and carbonation of ophiolitic peridotites (Eastern Desert,Egypt): constrains from stable isotopes and whole-rock geochemistry

Serpentinites are widespread in the Arabian-Nubian Shield (ANS) of the Eastern Desert of Egypt and usually enclose a tremendous carbonate alteration. Combined investigation of the stable isotope compositions of both O-H in serpentines and O-C in the whole-rock and the chemistry of the fluid-mobile elements (FMEs) in whole-rock serpentinites from Wadi (W.) Alam, Gabal (G.) El-Maiyit, and W. Atalla (Eastern Desert of Egypt) allowed to better understand the subsequent fluid sources of serpentinization and carbonation, as well as impact of these processes on the geochemistry of protolith ultramafic rocks. δ ¹⁸O values of W. Alam and W. Atalla serpentine minerals are close to the unaltered mantle and propose a lower temperature serpentinization if compared with those of G. El-Maiyit rocks. Moreover, δD values of W. Alam and W. Atalla serpentines (? 94 to ? 65‰) correspond to an igneous source that might be hydrothermal solutions mixed with the seawater in the mid-ocean ridge-arc transition setting. On the other hand, G. El-Maiyit serpentine is more depleted in ¹⁸O (with lower δ ¹⁸O values = 4.08–4.85‰), and its δD values (? 73 to 56 ‰) are most probably caused by an interaction with metamorphic fluids, acquired during on-land emplacement of oceanic peridotites or during burial in fore-arc setting. In addition, the oceanic oxygen isotope composition of most studied ophiolitic serpentinites points to the preservation of the pre-obduction δ ¹⁸O signatures and thus local-scale fluid flow at low water/rock ratios. Serpentinization fluids were CO₂-poor and the carbonation of the serpentinites resulted from infiltration of externally derived fluids. δ ¹⁸O_VSMOW values of carbonates in the studied serpentinites vary between heavier oxygen isotope composition in G. El-Maiyit samples (av. = 25.32‰) to lighter composition in W. Alam samples (av. = 19.43‰). However, δ ¹³C values of all serpentinites point mantle source of carbon. This source might have been evolved in mid-ocean ridge (W. Atalla) and subduction zone (W. Alam and G. El-Maiyit) settings. The studied serpentinites are usually enriched in FMEs, particularly Pb, Sr, Cs, and U. These enrichments were most probably the result of serpentinization and/or carbonation.     

Nickelpicromerite,K<Subscript>2</Subscript>Ni(SO<Subscript>4</Subscript>)<Subscript>2</Subscript>•6H<Subscript>2</Subscript>O,a new picromerite-group mineral from Slyudorudnik,South Urals,Russia

A new picromerite-group mineral, nickelpicromerite, K₂Ni(SO₄)₂?·?6H₂O (IMA 2012–053), was found at the Vein #169 of the Ufaley quartz deposit, near the town of Slyudorudnik, Kyshtym District, Chelyabinsk area, South Urals, Russia. It is a supergene mineral that occurs, with gypsum and goethite, in the fractures of slightly weathered actinolite-talc schist containing partially vermiculitized biotite and partially altered sulfides: pyrrhotite, pentlandite, millerite, pyrite and marcasite. Nickelpicromerite forms equant to short prismatic or tabular crystals up to 0.07 mm in size and anhedral grains up to 0.5 mm across, their clusters or crusts up to 1 mm. Nickelpicromerite is light greenish blue. Lustre is vitreous. Mohs hardness is 2–2½. Cleavage is distinct, parallel to {10–2}. D _meas is 2.20(2), D _calc is 2.22 g cm^?3. Nickelpicromerite is optically biaxial (+), α?=?1.486(2), β?=?1.489(2), γ?=?1.494(2), 2V_meas =75(10)°, 2V_calc =76°. The chemical composition (wt.%, electron-microprobe data) is: K₂O 20.93, MgO 0.38, FeO 0.07, NiO 16.76, SO₃ 37.20, H₂O (calc.) 24.66, total 100.00. The empirical formula, calculated based on 14 O, is: K_1.93Mg_0.04Ni_0.98S_2.02O_8.05(H₂O)_5.95. Nickelpicromerite is monoclinic, P2₁/c, a?=?6.1310(7), b?=?12.1863(14), c?=?9.0076(10) Å, β?=?105.045(2)°, V?=?649.9(1) ų, Z?=?2. Eight strongest reflections of the powder XRD pattern are [d,Å-I(hkl)]: 5.386–34(110); 4.312–46(002); 4.240–33(120); 4.085–100(012, 10–2); 3.685–85(031), 3.041–45(040, 112), 2.808–31(013, 20–2, 122), 2.368–34(13–3, 21–3, 033). Nickelpicromerite (single-crystal X-ray data, R?=?0.028) is isostructural to other picromerite-group minerals and synthetic Tutton’s salts. Its crystal structure consists of [Ni(H₂O)₆]²⁺ octahedra linked to (SO₄)^2? tetrahedra via hydrogen bonds. K⁺ cations are coordinated by eight anions. Nickelpicromerite is the product of alteration of primary sulfide minerals and the reaction of the acid Ni-sulfate solutions with biotite.     

The Influence of Benthic Macrofauna on the Erodibility of Intertidal Sediments with Varying mud Content in Three New Zealand Estuaries

Fine sediment inputs can alter estuarine ecosystem structure and function. However, natural variations in the processes that regulate sediment transport make it difficult to predict their fate. In this study, sediments were sampled at different times (2011–2012) from 45 points across intertidal sandflat transects in three New Zealand estuaries (Whitford, Whangamata, and Kawhia) encompassing a wide range in mud (≤63 μm) content (0–56 %) and macrofaunal community structure. Using a core-based erosion measurement device (EROMES), we calculated three distinct measures of sediment erosion potential: erosion threshold (? _c ; N m^?2), erosion rate (ER; g m^?2 s^?1), and change in erosion rate with increasing bed shear stress (m _e ; g N^?1 s^?1). Collectively, these measures characterized surface (? _c and ER) and sub-surface (m _e ) erosion. Benthic macrofauna were grouped by functional traits (size and motility) and data pooled across estuaries to determine relationships between abiotic (mud content, mean grain size) and biotic (benthic macrofauna, microbial biomass) variables and erosion measures. Results indicated that small bioturbating macrofauna (predominantly freely motile species <5 mm in size) destabilized surface sediments, explaining 23 % of the variation in ? _c (p ≤ 0.01) and 59 % of the variation in ER (p ≤ 0.01). Alternatively, mud content and mean grain size cumulatively explained 61 % of the variation in m _e (p ≤ 0.01), where increasing mud and grain size stabilized sub-surface sediments. These results highlight that the importance of biotic and abiotic predictors vary with erosion stage and that functional group classifications are a useful way to determine the impact of benthic macrofauna on sediment erodibility across communities with different species composition.     

Behavior of 10-Å phase at low temperatures

The thermal evolution of 10-Å phase Mg₃Si₄O₁₀(OH)₂·H₂O, a phyllosilicate which may have an important role in the storage/release of water in subducting slabs, was studied by X-ray single-crystal diffraction in the temperature range 116–293 K. The lattice parameters were measured at several intervals both on cooling and heating. The structural model was refined with intensity data collected at 116 K and compared to the model refined at room temperature. As expected for a layer silicate on cooling in this temperature range, the a and b lattice parameters undergo a small linear decrease, α _a  = 1.7(4) 10^?6 K^?1 and α _b  = 1.9(4) 10^?6 K^?1, where α is the linear thermal expansion coefficient. The greater variation is along the c axis and can be modeled with the second order polynomial c _T  = c ₂₉₃(1 + 6.7(4)10^?5 K^?1ΔT + 9.5(2.5)10^?8 K^?2(ΔT)²) where ΔT = T ? 293 K; the monoclinic angle β slightly increased. The cell volume thermal expansion can be modeled with the polynomial V _T  = V ₂₉₃ (1 + 8.0 10^?5 K^?1 ΔT + 1.4 10^?7 K^?2 (ΔT)²) where ΔT = T ? 293 is in K and V in ų. These variations were similar to those expected for a pressure increase, indicating that T and P effects are approximately inverse. The least-squares refinement with intensity data measured at 116 K shows that the volume of the SiO₄ tetrahedra does not change significantly, whereas the volume of the Mg octahedra slightly decreases. To adjust for the increased misfit between the tetrahedral and octahedral sheets, the tetrahedral rotation angle α changes from 0.58° to 1.38°, increasing the ditrigonalization of the silicate sheet. This deformation has implications on the H-bonds between the water molecule and the basal oxygen atoms. Furthermore, the highly anisotropic thermal ellipsoid of the H₂O oxygen indicates positional disorder, similar to the disorder observed at room temperature. The low-temperature results support the hypothesis that the disorder is static. It can be modeled with a splitting of the interlayer oxygen site with a statistical distribution of the H₂O molecules into two positions, 0.6 Å apart. The resulting shortest O_bas–O_W distances are 2.97 Å, with a significant shortening with respect to the value at room temperature. The low-temperature behavior of the H-bond system is consistent with that hypothesized at high pressure on the basis of the Raman spectra evolution with P.     

Crystal chemistry of natural and synthetic lead oxyhalides. Part I. Crystal structure of Pb<Subscript>13</Subscript>O<Subscript>10</Subscript>Cl<Subscript>6</Subscript>

Crystals of lead oxychloride Pb₁₃O₁₀Cl₆ have been synthesized on the basis of high-temperature solid-state reactions. The Pb₁₃O₁₀Cl₆ structure was studied using X-ray single-crystal diffraction analysis. The compound is monoclinic, and the space group is C2/c; the unit-cell dimensions are a = 16.1699(14), b = 7.0086(6), c = 23.578(2) Å, β = 97.75°, and V = 2647.6(4) ų. The structure has been solved by direct methods and refined to R ₁ = 0.0505 for 2671 observed unique reflections. The structure is a 3D framework consisting of OPb₄ tetrahedrons. Chlorine atoms are located in the framework channels. The structure contains seven symmetrically independent Pb atoms, which are coordinated by 2 to 4 O^2? and 2 to 4 Cl^? anions. The synthesized compound is compared with other natural and synthetic lead oxyhalides.     

Stability at high-pressure,elastic behaviour and pressure-induced structural evolution of CsAlSi<Subscript>5</Subscript>O<Subscript>12</Subscript>, a potential host for nuclear waste

The elastic and structural behaviour of the synthetic zeolite CsAlSi₅O₁₂ (a = 16.753(4), b = 13.797(3) and c = 5.0235(17) Å, space group Ama2, Z = 2) were investigated up to 8.5 GPa by in situ single-crystal X-ray diffraction with a diamond anvil cell under hydrostatic conditions. No phase-transition occurs within the P-range investigated. Fitting the volume data with a third-order Birch–Murnaghan equation-of-state gives: V ₀ = 1,155(4) ų, K _T0 = 20(1) GPa and K′ = 6.5(7). The “axial moduli” were calculated with a third-order “linearized” BM-EoS, substituting the cube of the individual lattice parameter (a ³, b ³, c ³) for the volume. The refined axial-EoS parameters are: a ₀ = 16.701(44) Å, K _T0a = 14(2) GPa (β_a = 0.024(3) GPa^?1), K′ _a = 6.2(8) for the a-axis; b ₀ = 13.778(20) Å, K _T0b = 21(3) GPa (β_b = 0.016(2) GPa^?1), K′ _b = 10(2) for the b-axis; c ₀ = 5.018(7) Å, K _T0c = 33(3) GPa (β_c = 0.010(1) GPa^?1), K′ _c = 3.2(8) for the c-axis (K _T0a:K _T0b:K _T0c = 1:1.50:2.36). The HP-crystal structure evolution was studied on the basis of several structural refinements at different pressures: 0.0001 GPa (with crystal in DAC without any pressure medium), 1.58(3), 1.75(4), 1.94(6), 3.25(4), 4.69(5), 7.36(6), 8.45(5) and 0.0001 GPa (after decompression). The main deformation mechanisms at high-pressure are basically driven by tetrahedral tilting, the tetrahedra behaving as rigid-units. A change in the compressional mechanisms was observed at P ≤ 2 GPa. The P-induced structural rearrangement up to 8.5 GPa is completely reversible. The high thermo-elastic stability of CsAlSi₅O_12, the immobility of Cs at HT/HP-conditions, the preservation of crystallinity at least up to 8.5 GPa and 1,000°C in elastic regime and the extremely low leaching rate of Cs from CsAlSi₅O₁₂ allow to consider this open-framework silicate as functional material potentially usable for fixation and deposition of Cs radioisotopes.     

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.     

CO<Subscript>2</Subscript> content of andesitic melts at graphite-saturated upper mantle conditions with implications for redox state of oceanic basalt source regions and remobilization of reduced carbon from subducted eclogite

We have performed experiments to determine the effects of pressure, temperature and oxygen fugacity on the CO₂ contents in nominally anhydrous andesitic melts at graphite saturation. The andesite composition was specifically chosen to match a low-degree partial melt composition that is generated from MORB-like eclogite in the convective, oceanic upper mantle. Experiments were performed at 1–3 GPa, 1375–1550?°C, and fO₂ of FMQ ?3.2 to FMQ ?2.3 and the resulting experimental glasses were analyzed for CO₂ and H₂O contents using FTIR and SIMS. Experimental results were used to develop a thermodynamic model to predict CO₂ content of nominally anhydrous andesitic melts at graphite saturation. Fitting of experimental data returned thermodynamic parameters for dissolution of CO₂ as molecular CO₂: ln(K ⁰) = ?21.79?±?0.04, ΔV ⁰?=?32.91?±?0.65 cm³mol^?1, ΔH ⁰?=?107?±?21 kJ mol^?1, and dissolution of CO₂ as CO₃ ^2?: ln(K ⁰ ) = ?21.38?±?0.08, ΔV ⁰?=?30.66?±?1.33 cm³ mol^?1, ΔH ⁰?=?42?±?37 kJ mol^?1, where K ⁰ is the equilibrium constant at some reference pressure and temperature, ΔV ⁰ is the volume change of reaction, and ΔH ⁰ is the enthalpy change of reaction. The thermodynamic model was used along with trace element partition coefficients to calculate the CO₂ contents and CO₂/Nb ratios resulting from the mixing of a depleted MORB and the partial melt of a graphite-saturated eclogite. Comparison with natural MORB and OIB data suggests that the CO₂ contents and CO₂/Nb ratios of CO₂-enriched oceanic basalts cannot be produced by mixing with partial melts of graphite-saturated eclogite. Instead, they must be produced by melting of a source containing carbonate. This result places a lower bound on the oxygen fugacity for the source region of these CO₂-enriched basalts, and suggests that fO₂ measurements made on cratonic xenoliths may not be applicable to the convecting upper mantle. CO₂-depleted basalts, on the other hand, are consistent with mixing between depleted MORB and partial melts of a graphite-saturated eclogite. Furthermore, calculations suggest that eclogite can remain saturated in graphite in the convecting upper mantle, acting as a reservoir for C.     

Heat capacity of hydrous trachybasalt from Mt Etna: comparison with CaAl<Subscript>2</Subscript>Si<Subscript>2</Subscript>O<Subscript>8</Subscript> (An)–CaMgSi<Subscript>2</Subscript>O<Subscript>6</Subscript> (Di) as basaltic proxy compositions

The specific heat capacity (C _p) of six variably hydrated (~3.5 wt% H₂O) iron-bearing Etna trachybasaltic glasses and liquids has been measured using differential scanning calorimetry from room temperature across the glass transition region. These data are compared to heat capacity measurements on thirteen melt compositions in the iron-free anorthite (An)–diopside (Di) system over a similar range of H₂O contents. These data extend considerably the published C _p measurements for hydrous melts and glasses. The results for the Etna trachybasalts show nonlinear variations in, both, the heat capacity of the glass at the onset of the glass transition (i.e., C _p ^g ) and the fully relaxed liquid (i.e., C _p ^l ) with increasing H₂O content. Similarly, the “configurational heat capacity” (i.e., C _p ^c  = C _p ^l  ? C _p ^g ) varies nonlinearly with H₂O content. The An–Di hydrous compositions investigated show similar trends, with C _p values varying as a function of melt composition and H₂O content. The results show that values in hydrous C _p ^g , C _p ^l and C _p ^c in the depolymerized glasses and liquids are substantially different from those observed for more polymerized hydrous albitic, leucogranitic, trachytic and phonolitic multicomponent compositions previously investigated. Polymerized melts have lower C _p ^l and C _p ^c and higher C _p ^g with respect to more depolymerized compositions. The covariation between C _p values and the degree of polymerization in glasses and melts is well described in terms of SM_hydrous and NBO/T _hydrous. Values of C _p ^c increase sharply with increasing depolymerization up to SM_hydrous ~ 30–35 mol% (NBO/T _hydrous ~ 0.5) and then stabilize to an almost constant value. The partial molar heat capacity of H₂O for both glasses (\( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{g}} \)) and liquids (\( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{l}} \)) appears to be independent of composition and, assuming ideal mixing, we obtain a value for \( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{l}} \) of 79 J mol^?1 K^?1. However, we note that a range of values for \( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{l}} \) (i.e., ~78–87 J mol^?1 K^?1) proposed by previous workers will reproduce the extended data to within experimental uncertainty. Our analysis suggests that more data are required in order to ascribe a compositional dependence (i.e., nonideal mixing) to \( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{l}} \).     

Petrogenesis of low-δ<Superscript>18</Superscript>O quartz porphyry dykes,Koegel Fontein complex,South Africa

This paper investigates the origin of low-δ¹⁸O quartz porphyry dykes associated with the 144–133 Ma Koegel Fontein Igneous Complex, which was intruded during the initial phase of breakup of Africa and South America. The 25-km diameter Rietpoort Granite is the largest and youngest phase of activity, and is roofed by a 10-km diameter pendant of gneiss. Quartz porphyry (QP) dykes, up to 15 m in width, strike NW–SE across the complex. The QP dykes that intruded outside the granite have similar quartz phenocryst δ¹⁸O values (average 8.0‰, ± 0.7, n?=?33) to the granite (average 8.3?±?1.0, n?=?7). The QP dykes that intruded the roof pendant have quartz phenocrysts with more variable δ¹⁸O values (average 1.6‰, ± 2.1, n?=?55). In some cases quartz phenocrysts have δ¹⁸O values as low as ? 2.5‰. The variation in δ¹⁸O value within the quartz crystal population of individual dykes is small relative to the overall range, and core and rim material from individual quartz phenocrysts in three samples are identical within error. There is no evidence that quartz phenocryst δ¹⁸O values have been affected by fluid–rock interaction. Based on a ?_quartz?magma value of 0.6‰, magma δ¹⁸O values must have been as low as ? 3.1‰. Samples collected along the length of the two main QP dykes that traverse the roof pendant have quartz phenocryst δ¹⁸O values that range from +?1.1 to +?4.6‰, and ? 2.3 to +?5.6‰, respectively. These δ¹⁸O values correlate negatively (r = ? 0.96) with initial ⁸⁷Sr/⁸⁶Sr, which can be explained by the event that lowered δ¹⁸O values of the source being older than the dykes. We suggest that the QP dykes were fed by magma produced by partial melting of gneiss, which had been variably altered at high temperature by ¹⁸O-depleted meteoric water during global glaciation at ~?550 Ma. The early melts had variable δ¹⁸O value but as melt pockets interconnected during melting, the δ¹⁸O values approached that of average gneiss. Variable quartz phenocryst δ¹⁸O values in the same dyke can be explained by vertical emplacement, at variable rates of ascent along the dyke. The lateral variation in quartz, and hence magma δ¹⁸O value at a particular point along a single dyke would depend on the rate of ascent of magma at that point along the dyke, and the ‘age’ of the particular magma batch.     

High-pressure displacive phase transition of a natural Mg-rich pigeonite

A high-pressure single-crystal X-ray diffraction study has been carried out on a P2₁/c natural Mg-rich pigeonite sample with composition ca. Wo₆En₇₆Fs₁₈ using a diamond anvil-cell. The unit-cell parameters were determined at 14 different pressures to 7.14 GPa. The sudden disappearance of the b-type reflections (h + k = odd) and a strong discontinuity (about 2.8%) in the unit-cell volume indicated a first-order P2₁/c–C2/c phase transition between 4.66 and 4.88 GPa. The P(V) data of the P2₁/c phase were fitted to 4.66 GPa by a third-order Birch–Murnaghan equation of state (BM3 EoS), whereas the limited number of experimental data collected within the C2/c phase between 4.88 and 7.14 GPa were fitted using the same equation of state but with K′ constrained to the value obtained for the P2₁/c fitting. The equation of state coefficients are V ₀ = 424.66(6) ų, K _T0 = 104(2) GPa and K′ = 8(1) for the P2₁/c phase, and V ₀ = 423.6(1) ų, K _T0 = 112.4(8) GPa, and K′ fixed to 8(1) for the C2/c phase. The axial moduli for a, b, and c for the P2₁/c phase were obtained using also a BM3-EoS, while for the C2/c phase only a linear calculation could be performed, and therefore the same approach was applied for comparison also to the P2₁/c phase. In general the C2/c phase exhibits axial compressibilities (β _c > β _a >> β _b) lower than those of the P2₁/c phase (β _b > β _c ≈ β _a; similar to those found in previous studies in clinopyroxenes and orthopyroxenes). The lower compressibility of the C2/c phase compared with that of the P2₁/c could be ascribed to the greater stiffness along the b direction. A previously published relationship between P _c and M2 average cation radius (i.r.) has been updated using all the literature data on P2₁/c clinopyroxene containing large cations at M2 site and our new data. The following weighted regression was obtained: P _c (GPa) = 26(4) ? 28(5) ×  i.r (Å), R ² = 0.97. This improved equation can be used to predict the critical pressure of natural P2₁/c clinopyroxene samples just knowing the composition at M2 site.     

Crystal structure and refined formula of garyansellite Mg<Subscript>2</Subscript>Fe<Superscript>3+</Superscript>(PO<Subscript>4</Subscript>)<Subscript>2</Subscript>(OH) · 2H<Subscript>2</Subscript>O

Single-crystal study of the structure (R = 0.0268) was performed for garyansellite from Rapid Creek, Yukon, Canada. The mineral is orthorhombic, Pbna, a = 9.44738(18), b = 9.85976(19), c = 8.14154(18) Å, V = 758.38(3) ų, Z = 4. An idealized formula of garyansellite is Mg₂Fe³⁺(PO₄)₂(OH) · 2H₂O. Structurally the mineral is close to other members of the phosphoferrite–reddingite group. The structure contains layers of chains of M(2)O₄(OH)(H₂O) octahedra which share edges to form dimers and connected by common edges with isolated from each other M(1)O₄(H₂O)₂ octahedra. The neighboring chains are connected to the layer through the common vertices of M(2) octahedra and octaahedral layers are linked through PO₄ tetrahedra.     

Influence of the <Emphasis Type="Italic">X</Emphasis>-site composition on tourmaline’s crystal structure: investigation of synthetic K-dravite,dravite, oxy-uvite,and magnesio-foitite using SREF and Raman spectroscopy

The crystal structures of synthetic K-dravite [^XK^YMg ₃ ^Z Al ₆ ^T Si₆O₁₈(BO₃) ₃ ^V (OH) ₃ ^W (OH)], dravite [^XNa^YMg ₃ ^Z Al ₆ ^T Si₆O₁₈(BO₃) ₃ ^V (OH) ₃ ^W (OH)], oxy-uvite [^XCa^YMg ₃ ^Z Al ₆ ^T Si₆O₁₈(BO₃) ₃ ^V (OH) ₃ ^W O], and magnesio-foitite [^X?^Y(Mg₂Al)^ZAl ₆ ^T Si₆O₁₈(BO₃) ₃ ^V (OH) ₃ ^W (OH)] are investigated by polarized Raman spectroscopy, single-crystal structure refinement (SREF), and powder X-ray diffraction. The use of compositionally simple tourmalines characterized by electron microprobe analysis facilitates the determination of site occupancy in the SREF and band assignment in the Raman spectra. The synthesized K-dravite, oxy-uvite, and magnesio-foitite have significant Mg–Al disorder between their octahedral sites indicated by their respective average 〈Y–O〉 and 〈Z–O〉 bond lengths. The Y- and Z-site compositions of oxy-uvite (^YMg_1.52Al_1.48(10) and ^ZAl_4.90Mg_1.10(15)) and magnesio-foitite (^YAl_1.62Mg_1.38(18) and ^ZAl_4.92Mg_1.08(24)) are refined from the electron densities at each site. The Mg–Al ratio of the Y and Z sites is also determined from the relative integrated peak intensities of the Raman bands in the O–H stretching vibrational range (3250–3850 cm^?1), producing values in good agreement with the SREF data. The unit cell volume of tourmaline increases from magnesio-foitite (1558.4(3) ų) to dravite (1569.5(4)–1571.7(3) ų) to oxy-uvite (1572.4(2) ų) to K-dravite (1588.1(2) ų), mainly due to lengthening of the crystallographic c-axis. The increase in the size of the X-site coordination polyhedron from dravite (Na) to K-dravite (K) is accommodated locally in the crystal structure, resulting in the shortening of the neighboring O1–H1 bond. In oxy-uvite, Ca²⁺ is locally associated with a deprotonated W (O1) site, whereas vacant X sites are neighbored by protonated W (O1) sites. Increasing the size of the X-site-occupying ion does not detectably affect bonding between the other sites; however, the higher charge of Ca and the deprotonated W (O1) site in oxy-uvite are correlated to changes in the lattice vibration Raman spectrum (100–1200 cm^?1), particularly for bands assigned to the T ₆O₁₈ ring. The Raman spectrum of magnesio-foitite shows significant deviations from those of K-dravite, dravite, and oxy-uvite in both the lattice and O–H stretching vibrational ranges (100–1200 and 3250–3850 cm^?1, respectively). The vacant X site is correlated with long- and short-range changes in the crystal structure, i.e., deformation of the T ₆O₁₈ ring and lengthening of the O1–H1 and O3–H3 bonds. However, X-site vacancies in K-dravite, dravite, and oxy-uvite result only in the lengthening of the neighboring O1–H1 bond and do not result in identifiable changes in the lattice-bonding environment.     

Sulfate Reduction and Sulfur Cycles at Two Seagrass Beds Inhabited by Cold Affinity <Emphasis Type="Italic">Zostera marina</Emphasis> and Warm Affinity <Emphasis Type="Italic">Halophila nipponica</Emphasis> in Temperate Coastal Waters

To evaluate the impact of invading seagrass on biogeochemical processes associated with sulfur cycles, we investigated the geochemical properties and sulfate reduction rates (SRRs) in sediments inhabited by invasive warm affinity Halophila nipponica and indigenous cold affinity Zostera marina. A more positive relationship between SRR and below-ground biomass (BGB) was observed at the H. nipponica bed (SRR = 0.6809 × BGB ? 4.3162, r ² = 0.9878, p = 0.0006) than at the Z. marina bed (SRR = 0.3470 × BGB ? 4.0341, r ² = 0.7082, p = 0.0357). These results suggested that SR was more stimulated by the dissolved organic carbon (DOC) exuded from the roots of H. nipponica than by the DOC released from the roots of Z. marina. Despite the enhanced SR in spring-summer, the relatively lower proportion (average, 20%) of acid-volatile sulfur (AVS) in total reduced sulfur and the strong correlation between total oxalate-extractable Fe (Fe_(oxal)) and chromium-reducible sulfur (CRS = 0.2321 × total Fe_(oxal) + 1.8180, r ² = 0.3344, p = 0.0076) in the sediments suggested the rapid re-oxidation of sulfide and precipitation of sulfide with Fe. The turnover rate of the AVS at the H. nipponica bed (0.13 day^?1) was 2.5 times lower than that at the Z. marina bed (0.33 day^?1). Together with lower AVS turnover, the stronger correlation of SRR to BGB in the H. nipponica bed suggests that the extension of H. nipponica resulting from the warming of seawater might provoke more sulfide accumulation in coastal sediments.     

Synthesis,TEM characterization and thermal behaviour of LiNiSi<Subscript>2</Subscript>O<Subscript>6</Subscript> pyroxene

A pyroxene with composition LiNiSi₂O₆ was synthesized at T = 1,473 K and P = 2.0 GPa; the cell parameters at T = 298 K are a = 9.4169(6) Å, b = 8.4465(7) Å, c = 5.2464(3) Å, β = 110.534(6)°, V = 390.78(3) Å³. TEM examination of the LiNiSi₂O₆ pyroxene showed the presence of h + k odd reflections indicative of a primitive lattice, and of antiphase domains obtained by dark field imaging of the h + k odd reflections. A HT in situ investigation was performed by examining TEM selected area diffraction patterns collected at high temperature and synchrotron radiation powder diffraction. In HTTEM the LiNiSi₂O₆ was examined together with LiCrSi₂O₆ pyroxene. In LiCrSi₂O₆ the h + k odd critical reflections disappear at about 340 K; they are sharp up to the transition temperature and do not change their shape until they disappear. In LiNiSi₂O₆ the h + k odd reflections are present up to sample deterioration at 650 K. A high temperature synchrotron radiation powder diffraction investigation was performed on LiNiSi₂O₆ between 298 and 773 K. The analysis of critical reflections and of changes in cell parameters shows that the space group is P-centred up to the highest temperature. The comparative analysis of the thermal and spontaneous strain contributions in P2₁/c and C2/c pyroxenes indicates that the high temperature strain in P-LiNiSi₂O₆ is very similar to that due to thermal strain only in C2/c spodumene and that a spontaneous strain contribution related to pre-transition features is not apparent in LiNiSi₂O₆. A different high-temperature behaviour in LiNiSi₂O₆ with respect to other pyroxenes is suggested, possibly in relation with the presence of Jahn–Teller distortion of the M1 polyhedron centred by low-spin Ni³⁺.     

Theoretical calculation of equilibrium Mg isotope fractionation between silicate melt and its vapor

Isotope fractionation during the evaporation of silicate melt and condensation of vapor has been widely used to explain various isotope signals observed in lunar soils, cosmic spherules, calcium–aluminum-rich inclusions, and bulk compositions of planetary materials. During evaporation and condensation, the equilibrium isotope fractionation factor (α) between high-temperature silicate melt and vapor is a fundamental parameter that can constrain the melt’s isotopic compositions. However, equilibrium α is difficult to calibrate experimentally. Here we used Mg as an example and calculated equilibrium Mg isotope fractionation in MgSiO₃ and Mg₂SiO₄ melt–vapor systems based on first-principles molecular dynamics and the high-temperature approximation of the Bigeleisen–Mayer equation. We found that, at 2500 K, δ²⁵Mg values in the MgSiO₃ and Mg₂SiO₄ melts were 0.141?±?0.004 and 0.143?±?0.003‰ more positive than in their respective vapors. The corresponding δ²⁶Mg values were 0.270?±?0.008 and 0.274?±?0.006‰ more positive than in vapors, respectively. The general \(\alpha - T\) equations describing the equilibrium Mg α in MgSiO₃ and Mg₂SiO₄ melt–vapor systems were: \(\alpha_{{{\text{Mg}}\left( {\text{l}} \right) - {\text{Mg}}\left( {\text{g}} \right)}} = 1 + \frac{{5.264 \times 10^{5} }}{{T^{2} }}\left( {\frac{1}{m} - \frac{1}{{m^{\prime}}}} \right)\) and \(\alpha_{{{\text{Mg}}\left( {\text{l}} \right) - {\text{Mg}}\left( {\text{g}} \right)}} = 1 + \frac{{5.340 \times 10^{5} }}{{T^{2} }}\left( {\frac{1}{m} - \frac{1}{{m^{\prime}}}} \right)\), respectively, where m is the mass of light isotope ²⁴Mg and m′ is the mass of the heavier isotope, ²⁵Mg or ²⁶Mg. These results offer a necessary parameter for mechanistic understanding of Mg isotope fractionation during evaporation and condensation that commonly occurs during the early stages of planetary formation and evolution.     

Influence of oxygen fugacity and temperature on the redox state of iron in natural silicic aluminosilicate melts

The influence of oxygen fugacity (fO₂) and temperature on the valence and structural state of iron was experimentally studied in glasses quenched from natural aluminosilicate melts of granite and pantellerite compositions exposed to various T-fO₂ conditions (1100–1420°C and 10^?12–10^?0.68 bar) at a total pressure of 1 atm. The quenched glasses were investigated by Mössbauer spectroscopy. It was shown that the effect of oxygen fugacity on the redox state of iron at 1320–1420°C can be described by the equation log(Fe³⁺/Fe²⁺) = k log(fO₂) + q, where k and q are constants depending on melt composition and temperature. The Fe³⁺/Fe²⁺ ratio decreases with decreasing fO₂ (T = const) and increasing temperature (fO₂ = const). The structural state of Fe³⁺ depends on the degree of iron oxidation. With increasing Fe³⁺/Fe²⁺ ≥ 1, the dominant coordination of Fe³⁺ changes from octahedral to tetrahedral. Ferrous iron ions occur in octahedral (and/or five-coordinated) sites independent of Fe³⁺/Fe²⁺.