Structure, thermodynamic and transport properties of CaAl₂Si₂O₈ liquid. Part I: Molecular dynamics simulations     

Frank J. Spera  Dean Nevins  Ian Cutler 《Geochimica et cosmochimica acta》2009,73(22):6918-2191

  相似文献   

Water diffusion in dacitic melt     

Huaiwei Ni  Harald Behrens 《Geochimica et cosmochimica acta》2009,73(12):3642-3655

H₂O diffusion in dacitic melt was investigated at 0.48-0.95 GPa and 786-893 K in a piston-cylinder apparatus. The diffusion couple design was used, in which a nominally dry dacitic glass makes one half and is juxtaposed with a hydrous dacitic glass containing up to ∼8 wt.% total water (H₂O_t). H₂O concentration profiles were measured on quenched glasses with infrared microspectroscopy. The H₂O diffusivity in dacite increases rapidly with water content under experimental conditions, similar to previous measurements at the same temperature but at pressure <0.15 GPa. However, compared with the low-pressure data, H₂O diffusion at high pressure is systematically slower. H₂O diffusion profiles in dacite can be modeled by assuming molecular H₂O (H₂O_m) is the diffusing species. Total H₂O diffusivity D_H2Ot within 786-1798 K, 0-1 GPa, and 0-8 wt.% H₂O_t can be expressed as: where D_H2Ot is in m²/s, T is temperature in K, P is pressure in GPa, K = exp(1.49 − 2634/T) is the equilibrium constant of speciation reaction (H₂O_m+O?2OH) in the melt, X = C/18.015/[C/18.015 + (100 − C)/33.82], C is wt.% of H₂O_t, and 18.015 and 33.82 g/mol correspond to the molar masses of H₂O and anhydrous dacite on a single oxygen basis. Compared to H₂O diffusion in rhyolite, diffusivity in dacite is lower at intermediate temperatures but higher at superliquidus temperatures. This general H₂O diffusivity expression can be applied to a broad range of geological conditions, including both magma chamber processes and volcanic eruption dynamics from conduit to the surface.  相似文献   

Structure, thermodynamic and transport properties of liquid MgSiO₃: Comparison of molecular models and laboratory results     

Frank J. Spera  Mark S. Ghiorso 《Geochimica et cosmochimica acta》2011,75(5):1272-1296

Liquid MgSiO₃ is a model for the Earth’s magma ocean and of remnant melt present near the core-mantle boundary. Here, models for molten MgSiO₃ are computed employing empirical potential molecular dynamics (EPMD) and results are compared to published results including two EPMD studies and three first-principles molecular dynamics (FPMD) models and to laboratory data. The EPMD results derived from the Oganov (OG) potential come closest to the density of MgSiO₃ liquid at the 1-bar melting point inferred from the melting curve. At higher P, EPMD densities calculated from the OG potential and FPMD broadly match shock wave studies, with the OG potential yielding the better comparison. Matsui (M) potential results deviate from other studies above ∼50 GPa. Overall, results based on the OG potential compare best to experimental densities over the P-T range of the mantle. Isothermally, upon increasing P the mean coordination numbers () of oxygen around Si and Mg monotonically increase with pressure. Tetrahedral Si and octahedral Si monotonically increase and decrease, respectively, whereas pentahedral Si maximizes at 10-20 GPa. Tetrahedral Mg decreases monotonically as P increases whereas pentahedral, octahedral and higher coordination polyhedra each show similar behavior first increasing and then decreasing after attaining a maximum; the P of the maximum for each polyhedra type migrates to higher P as the CN increases. Free oxygen and oxygen with one nearest neighbor of either Si or Mg decreases whereas Si or Mg with two or three nearest oxygens (i.e., tricluster oxygen) increases with increasing P isothermally. The increase of tricluster oxygen is consistent with spectroscopy on MgSiO₃ glass quenched from 2000 K and 0-40 GPa and high-energy X-ray studies constraining the coordination of O around Mg and around Si at 2300 K and 1 bar. Coordination statistics from FPMD studies for O around Si and Si around O are in agreement with the EPMD results based on the M and OG potentials. Mg self-diffusivity is greater than O and Si self-diffusivities for both the M and OG potentials. All D values monotonically decrease with increasing pressure isothermally and all atoms are more diffusive in the M liquid compared to the OG liquid except at T > ∼5000 K and P > 100 GPa. Previously published EPMD diffusivities fall between values given by the M and OG potentials, at least up to 45 GPa. The M liquid is generally less viscous than the OG liquid except at P > ∼80 GPa. Activation energy and volume are around 96 kJ/mol and 1.5 cm³/mol, respectively. The FPMD viscosity results at 120 GPa and 4000 and 4500 K are essentially identical to the values from the M and OG potentials. FPMD viscosity results are similar to the OG results for P < 60 GPa; at higher P, the FPMD viscosities are higher. At 4000 K and 100 GPa the shear viscosity of liquid MgSiO₃ is ∼0.1 Pa s. More extensive laboratory results are required to better define the thermodynamic, transport and structural properties of MgSiO₃ liquids and for comparison with computational studies.  相似文献   

High-pressure phase equilibria and element partitioning experiments on Apollo 15 green C picritic glass: Implications for the role of garnet in the deep lunar interior     

David S. Draper  S. Andrew duFrane  Rachel E. Dwarzski 《Geochimica et cosmochimica acta》2006,70(9):2400-2416

We report results of nominally anhydrous near-liquidus experiments on a synthetic analog to very low-titanium Apollo 15 green C lunar picritic glass from ∼2 to 5 GPa. Apollo 15 green C glass (A15C) is saturated with garnet and pyroxene on the liquidus at ∼3 GPa. However, such an assemblage is unlikely to represent the lunar-mantle source region for this glass, and instead an olivine + orthopyroxene-dominated source is favored, in accord with earlier lower-pressure experiments on A15C. Near-liquidus garnet has a slight but significant majorite component at ∼5 GPa in this iron-rich bulk composition, as expected from our previous work in ordinary-chondritic bulk compositions. Ion microprobe measurements of partitioning of Sr, Ba, Sc, Nd, Sm, Dy, Yb, Y, Zr, Hf, and Th between garnet and coexisting melt in these experiments are the first garnet partition coefficients (D values) available that are directly relevant to lunar compositions. D values for these garnets differ significantly compared to D values for garnets grown in more magnesian, terrestrial bulk compositions, which until now are all that have been available in modeling the possible role of garnet in the lunar interior. For example, D values for heavy rare earth elements are lower than are those from terrestrial basaltic systems. These partitioning values are well-described by the lattice-strain partitioning model, but predictive relationships for garnet partitioning using that model fail to match the measured values, as was the case in our earlier work on chondritic compositions. Using our new D values in place of the “terrestrial” values in a variety of models of lunar petrogenesis, we suggest that garnet is unlikely to be present in the source regions for very titanium-poor lunar liquids despite its appearance on the liquidus of A15C.  相似文献   

Titanium in subduction zone fluids: First insights from abinitio molecular metadynamics simulations     

Jelle van Sijl  Neil L. Allan  Wim van Westrenen 《Geochimica et cosmochimica acta》2010,74(9):2797-7377

The preferential incorporation of High-Field-Strength Elements (HFSE) in rutile (TiO₂), combined with its supposed stability in subduction zone settings, make it an ideal candidate to explain the low HFSE concentrations in subduction-derived magmas. The solubility behaviour of rutile is key to these arguments, but at present experimental and field-based evidence are contradictory.We have used abinitio molecular (meta)dynamics to investigate the coordination environment of Ti(IV) in pure water at 300 and 1000 K and densities ranging from 900-1260 kg m⁻³ (approximate pressures 0.9-3.6 GPa). In all high temperature simulations, the long-range structure of the solvent indicates a breakdown of the hydrogen bonding network as expected for supercritical water. The five-fold coordination of titanium to water is energetically most favourable in aqueous fluids at room temperature and pressure, separated from four and six-fold configurations by ∼175 and ∼200 kJ mol⁻¹, respectively. The average first shell Ti-O distance is 2.00 Å, in excellent agreement with bond lengths obtained from experiments. At similar densities and 1000 K, titanium is on average six-fold coordinate with water, and shows some degree of water dissociation in the first hydration shell. This coordination environment is remarkably persistent with increasing density from 1021 to 1260 kg m⁻³ at constant temperature (1000 K). At lower densities, however, (900 kg m⁻³ at 1000 K), the coordination with first shell water molecules is less than five. The observed coordination changes could promote association of titanium with peralkaline or peraluminous domains in the aqueous fluid, and thereby explain field-and laboratory based evidence of enhanced HFSE concentrations.This study demonstrates that abinitio molecular dynamics has considerable potential to access details of element behaviour in aqueous fluids at geologically relevant conditions that are impossible to examine otherwise. Changes in the solvent structure due to density variations lead to differences in solvent behaviour allowing access to new domains for fluid-solid interaction. Moreover, changes in the solvent structure are strongly linked to the effectiveness of element solvation.  相似文献   

Raman spectroscopic measurements of CO₂ density: Experimental calibration with high-pressure optical cell (HPOC) and fused silica capillary capsule (FSCC) with application to fluid inclusion observations   总被引：3，自引：0，他引：3  

Xiaolin Wang  I.-Ming Chou  Robert C. Burruss  Yucai Song 《Geochimica et cosmochimica acta》2011,75(14):4080-4093

Raman spectroscopy is a powerful method for the determination of CO₂ densities in fluid inclusions, especially for those with small size and/or low fluid density. The relationship between CO₂ Fermi diad split (Δ, cm⁻¹) and CO₂ density (ρ, g/cm³) has been documented by several previous studies. However, significant discrepancies exist among these studies mainly because of inconsistent calibration procedures and lack of measurements for CO₂ fluids having densities between 0.21 and 0.75 g/cm³, where liquid and vapor phases coexist near room temperature.In this study, a high-pressure optical cell and fused silica capillary capsules were used to prepare pure CO₂ samples with densities between 0.0472 and 1.0060 g/cm³. The measured CO₂ Fermi diad splits were calibrated with two well established Raman bands of benzonitrile at 1192.6 and 1598.9 cm⁻¹. The relationship between the CO₂ Fermi diad split and density can be represented by: ρ = 47513.64243 − 1374.824414 × Δ + 13.25586152 × Δ² − 0.04258891551 × Δ³ (r² = 0.99835, σ = 0.0253 g/cm³), and this relationship was tested by synthetic fluid inclusions and natural CO₂-rich fluid inclusions. The effects of temperature and the presence of H₂O and CH₄ on this relationship were also examined.  相似文献   

Complexation of cadmium to sulfur and oxygen functional groups in an organic soil   总被引：7，自引：0，他引：7  

Torbjörn Karlsson  Erik Björn 《Geochimica et cosmochimica acta》2007,71(3):604-614

Cadmium (Cd) is a toxic trace element and due to human activities soils and waters are contaminated by Cd both on a local and global scale. It is widely accepted that chemical interactions with functional groups of natural organic matter (NOM) is vital for the bioavailability and mobility of trace elements. In this study the binding strength of cadmium (Cd) to soil organic matter (SOM) was determined in an organic (49% organic C) soil as a function of reaction time, pH and Cd concentration. In experiments conducted at native Cd concentrations in soil (0.23 μg g⁻¹ dry soil), halides (Cl, Br) were used as competing ligands to functional groups in SOM. The concentration of Cd in the aqueous phase was determined by isotope-dilution (ID) inductively-coupled-plasma-mass-spectrometry (ICP-MS), and the activity of Cd²⁺ was calculated from the well-established Cd-halide constants. At higher Cd loading (500-54,000 μg g⁻¹), the Cd²⁺ activity was directly determined by an ion-selective electrode (ISE). On the basis of results from extended X-ray absorption fine structure (EXAFS) spectroscopy, a model with one thiolate group (RS⁻) was used to describe the complexation (Cd²⁺ + RS⁻ ? CdSR⁺; log K_CdSR) at native Cd concentrations. The concentration of thiols (RSH; 0.047 mol kg⁻¹ C) was independently determined by X-ray absorption near-edge structure (XANES) spectroscopy. Log K_CdSR values of 11.2-11.6 (pK_a for RSH = 9.96), determined in the pH range 3.1-4.6, compare favorably with stability constants for the association between Cd and well-defined thiolates like glutathione. In the concentration range 500-54,000 μg Cd g⁻¹, a model consisting of one thiolate and one carboxylate (RCOO⁻) gave the best fit to data, indicating an increasing role for RCOOH groups as RSH groups become saturated. The determined log K_CdOOCR of 3.2 (Cd²⁺ +  RCOO⁻ ? CdOOCR⁺; log K_CdOOCR; pK_a for RCOOH = 4.5) is in accordance with stability constants determined for the association between Cd and well-defined carboxylates. Given a concentration of reduced sulfur groups of 0.2% or higher in NOM, we conclude that the complexation to organic RSH groups may control the speciation of Cd in soils, and most likely also in surface waters, with a total concentration less than 5 mg Cd g⁻¹ organic C.  相似文献   

Ni(II) sorption on natural chlorite     

Åsa Zazzi  Anna-Maria Jakobsson  Susanna Wold 《Applied Geochemistry》2012

Sorption of Ni(II) onto chlorite surfaces was studied as a function of pH (5–10), ionic strength (0.01–0.5 M) and Ni concentration (10⁻⁸–10⁻⁶ M) in an Ar atmosphere using batch sorption with radioactive ⁶³Ni as tracer. Such studies are important since Ni(II) is one of the major activation products in spent nuclear fuel and sorption data on minerals such as chlorite are lacking. The sorption of Ni(II) onto chlorite was dependent on pH but not ionic strength, which indicates that the process primarily comprises sorption by surface complexation. The maximum sorption was at pH ∼ 8 (K_d = ∼10⁻³ cm³/g). Desorption studies over a period of 1–2 weeks involving replacement of the aqueous solution indicated a low degree of desorption. The acid–base properties of the chlorite mineral were determined by titration and described using a non-electrostatic surface complexation model in FITEQL. A 2-pK NEM model and three surface complexes, Chl_OHNi²⁺, Chl_OHNi(OH)⁺ and Chl_OHNi(OH)₂, gave the best fit to the sorption results using FITEQL. The high K_d values and low degree of desorption observed indicate that under expected groundwater conditions, a large fraction of Ni(II) that is potentially leachable from spent nuclear fuel may be prevented from migrating by sorption onto chlorite surfaces.  相似文献   

Hydration of rhyolitic glass during weathering as characterized by IR microspectroscopy   总被引：1，自引：0，他引：1  

Tadashi Yokoyama  Satoshi Okumura  Satoru Nakashima 《Geochimica et cosmochimica acta》2008,72(1):117-125

The mechanism and rate of hydration of rhyolitic glass during weathering were studied. Doubly polished thin sections of two rhyolites with different duration of weathering (Ohsawa lava: 26,000 yr, Awanomikoto lava: 52,000 yr) were prepared. Optical microscope observation showed that altered layers had developed along the glass surfaces. IR spectral line profile analysis was conducted on the glass sections from the surface to the interior for a length of 250 μm and the contents of molecular H₂O (H₂O_m), OH species (OH) and total water (H₂O_t) were determined. The diffusion profile of H₂O_m in Ohsawa lava extends beyond the layer observed by optical microscope. The content of H₂O_m in the hydrated region is much higher than that of OH species. Thus, the reaction from H₂O_m to OH appears to be little and H₂O_m is the dominant water species moving into the glass during weathering. Based on the concentration profiles, the diffusion coefficients of H₂O_m(D_H2Om) and H₂O_t(D_H2Ot) were determined to be 2.8 × 10⁻¹⁰ and 3.4 × 10⁻¹⁰ μm² s⁻¹ for Ohsawa lava, and 5.2 × 10⁻¹¹ and 4.1 × 10⁻¹¹ μm² s⁻¹ for Awanomikoto lava, respectively. The obtained D_H2Om during weathering are more than 2-3 orders of magnitude larger than the diffusion coefficient at ∼20 °C that is extrapolated from the diffusivity data for >400 °C. This might suggest that the mechanism of water transport is different at weathering conditions and >400 °C.  相似文献   

Spectroscopic analysis (FTIR, Raman) of water in mafic and intermediate glasses and glass inclusions   总被引：1，自引：0，他引：1  

Maxime Mercier  Andrea Di Muro  Nicole Métrich  Olfa Belhadj 《Geochimica et cosmochimica acta》2010,74(19):5641-5656

Micro-Raman spectroscopy, even though a very promising technique, is not still routinely applied to analyse H₂O in silicate glasses. The accuracy of Raman water determinations critically depends on the capability to predict and take into account both the matrix effects (bulk glass composition) and the analytical conditions on band intensities. On the other hand, micro-Fourier transform infrared spectroscopy is commonly used to measure the hydrous absorbing species (e.g., hydroxyl OH⁻ and molecular H₂O) in natural glasses, but requires critical assumptions for the study of crystal-hosted glasses. Here, we quantify for the first time the matrix effect of Raman external calibration procedures for the quantification of the total H₂O content (H₂O_T = OH⁻ + H₂O_m) in natural silicate glasses. The procedures are based on the calibration of either the absolute (external calibration) or scaled (parameterisation) intensity of the 3550 cm⁻¹ band. A total of 67 mafic (basanite, basalt) and intermediate (andesite) glasses hosted in olivines, having between 0.2 and 4.8 wt% of H₂O, was analysed. Our new dataset demonstrates, for given water content, the height (intensity) of Raman H₂O_T band depends on glass density, reflectance and water environment. Hence this matrix effect must be considered in the quantification of H₂O by Raman spectroscopy irrespective of the procedure, whereas the parameterisation mainly helps to predict and verify the self-consistency of the Raman results. In addition, to validate the capability of the micro-Raman to accurately determine the H₂O content of multicomponent aluminosilicate glasses, a subset of 23 glasses was analysed by both micro-Raman and micro-FTIR spectroscopy using the band at 3550 cm⁻¹. We provide new FTIR absorptivity coefficients (ε₃₅₅₀) for basalt (62.80 ± 0.8 L mol⁻¹ cm⁻¹) and basanite (43.96 ± 0.6 L mol⁻¹ cm⁻¹). These values, together with an exhaustive review of literature data, confirm the non-linear decline of the FTIR absorptivity coefficient (ε₃₅₅₀) as the glass depolymerisation increases. We demonstrate the good agreement between micro-FTIR and micro-Raman determination of H₂O in silicate glasses when the matrix effects are properly considered.  相似文献   

A surface structural model for ferrihydrite I: Sites related to primary charge, molar mass, and mass density     

Tjisse Hiemstra  Willem H. Van Riemsdijk 《Geochimica et cosmochimica acta》2009,73(15):4423-455

A multisite surface complexation (MUSIC) model for ferrihydrite (Fh) has been developed. The surface structure and composition of Fh nanoparticles are described in relation to ion binding and surface charge development. The site densities of the various reactive surface groups, the molar mass, the mass density, the specific surface area, and the particle size are quantified. As derived theoretically, molecular mass and mass density of nanoparticles will depend on the types of surface groups and the corresponding site densities and will vary with particle size and surface area because of a relatively large contribution of the surface groups in comparison to the mineral core of nanoparticles. The nano-sized (∼2.6 nm) particles of freshly prepared 2-line Fh as a whole have an increased molar mass of M ∼ 101 ± 2 g/mol Fe, a reduced mass density of ∼3.5 ± 0.1 g/cm³, both relatively to the mineral core. The specific surface area is ∼650 m²/g. Six-line Fh (5-6 nm) has a molar mass of M ∼ 94 ± 2 g/mol, a mass density of ∼3.9 ± 0.1 g/cm³, and a surface area of ∼280 ± 30 m²/g. Data analysis shows that the mineral core of Fh has an average chemical composition very close to FeOOH with M ∼ 89 g/mol. The mineral core has a mass density around ∼4.15 ± 0.1 g/cm³, which is between that of feroxyhyte, goethite, and lepidocrocite. These results can be used to constrain structural models for Fh. Singly-coordinated surface groups dominate the surface of ferrihydrite (∼6.0 ± 0.5 nm⁻²). These groups can be present in two structural configurations. In pairs, the groups either form the edge of a single Fe-octahedron (∼2.5 nm⁻²) or are present at a single corner (∼3.5 nm⁻²) of two adjacent Fe octahedra. These configurations can form bidentate surface complexes by edge- and double-corner sharing, respectively, and may therefore respond differently to the binding of ions such as uranyl, carbonate, arsenite, phosphate, and others. The relatively low PZC of ferrihydrite can be rationalized based on the estimated proton affinity constant for singly-coordinated surface groups. Nanoparticles have an enhanced surface charge. The charging behavior of Fh nanoparticles can be described satisfactory using the capacitance of a spherical Stern layer condenser in combination with a diffuse double layer for flat plates.  相似文献   

CO₂ solubility in Martian basalts and Martian atmospheric evolution     

Ben D. Stanley  Marc M. Hirschmann  Anthony C. Withers 《Geochimica et cosmochimica acta》2011,75(20):5987-6003

To understand possible volcanogenic fluxes of CO₂ to the Martian atmosphere, we investigated experimentally carbonate solubility in a synthetic melt based on the Adirondack-class Humphrey basalt at 1-2.5 GPa and 1400-1625 °C. Starting materials included both oxidized and reduced compositions, allowing a test of the effect of iron oxidation state on CO₂ solubility. CO₂ contents in experimental glasses were determined using Fourier transform infrared spectroscopy (FTIR) and Fe³⁺/Fe^T was measured by Mössbauer spectroscopy. The CO₂ contents of glasses show no dependence on Fe³⁺/Fe^T and range from 0.34 to 2.12 wt.%. For Humphrey basalt, analysis of glasses with gravimetrically-determined CO₂ contents allowed calibration of an integrated molar absorptivity of 81,500 ± 1500 L mol⁻¹ cm⁻² for the integrated area under the carbonate doublet at 1430 and 1520 cm⁻¹. The experimentally determined CO₂ solubilities allow calibration of the thermodynamic parameters governing dissolution of CO₂ vapor as carbonate in silicate melt, K_II, (Stolper and Holloway, 1988) as follows: , ΔV⁰ = 20.85 ± 0.91 cm³ mol⁻¹, and ΔH⁰ = −17.96 ± 10.2 kJ mol⁻¹. This relation, combined with the known thermodynamics of graphite oxidation, facilitates calculation of the CO₂ dissolved in magmas derived from graphite-saturated Martian basalt source regions as a function of P, T, and f_O2. For the source region for Humphrey, constrained by phase equilibria to be near 1350 °C and 1.2 GPa, the resulting CO₂ contents are 51 ppm at the iron-wüstite buffer (IW), and 510 ppm at one order of magnitude above IW (IW + 1). However, solubilities are expected to be greater for depolymerized partial melts similar to primitive shergottite Yamato 980459 (Y 980459). This, combined with hotter source temperatures (1540 °C and 1.2 GPa) could allow hot plume-like magmas similar to Y 980459 to dissolve 240 ppm CO₂ at IW and 0.24 wt.% of CO₂ at IW + 1. For expected magmatic fluxes over the last 4.5 Ga of Martian history, magmas similar to Humphrey would only produce 0.03 and 0.26 bars from sources at IW and IW + 1, respectively. On the other hand, more primitive magmas like Y 980459 could plausibly produce 0.12 and 1.2 bars at IW and IW + 1, respectively. Thus, if typical Martian volcanic activity was reduced and the melting conditions cool, then degassing of CO₂ to the atmosphere may not be sufficient to create greenhouse conditions required by observations of liquid surface water. However, if a significant fraction of Martian magmas derive from hot and primitive sources, as may have been true during the formation of Tharsis in the late Noachian, that are also slightly oxidized (IW + 1.2), then significant contribution of volcanogenic CO₂ to an early Martian greenhouse is plausible.  相似文献   

Surface heat flow and CO₂ emissions within the Ohaaki hydrothermal field,Taupo Volcanic Zone,New Zealand     

Clinton Rissmann  Bruce Christenson  Cynthia Werner  Matthew Leybourne  Jim Cole  Darren Gravley 《Applied Geochemistry》2012

Carbon dioxide emissions and heat flow have been determined from the Ohaaki hydrothermal field, Taupo Volcanic Zone (TVZ), New Zealand following 20 a of production (116 MW_e). Soil CO₂ degassing was quantified with 2663 CO₂ flux measurements using the accumulation chamber method, and 2563 soil temperatures were measured and converted to equivalent heat flow (W m⁻²) using published soil temperature heat flow functions. Both CO₂ flux and heat flow were analysed statistically and then modelled using 500 sequential Gaussian simulations. Forty subsoil CO₂ gas samples were also analysed for stable C isotopes. Following 20 a of production, current CO₂ emissions equated to 111 ± 6.7 T/d. Observed heat flow was 70 ± 6.4 MW, compared with a pre-production value of 122 MW. This 52 MW reduction in surface heat flow is due to production-induced drying up of all alkali–Cl outflows (61.5 MW) and steam-heated pools (8.6 MW) within the Ohaaki West thermal area (OHW). The drying up of all alkali–Cl outflows at Ohaaki means that the soil zone is now the major natural pathway of heat release from the high-temperature reservoir. On the other hand, a net gain in thermal ground heat flow of 18 MW (from 25 MW to 43.3 ± 5 MW) at OHW is associated with permeability increases resulting from surface unit fracturing by production-induced ground subsidence. The Ohaaki East (OHE) thermal area showed no change in distribution of shallow and deep soil temperature contours despite 20 a of production, with an observed heat flow of 26.7 ± 3 MW and a CO₂ emission rate of 39 ± 3 T/d. The negligible change in the thermal status of the OHE thermal area is attributed to the low permeability of the reservoir beneath this area, which has limited production (mass extraction) and sheltered the area from the pressure decline within the main reservoir. Chemistry suggests that although alkali–Cl outflows once contributed significantly to the natural surface heat flow (∼50%) they contributed little (<1%) to pre-production CO₂ emissions due to the loss of >99% of the original CO₂ content due to depressurisation and boiling as the fluids ascended to the surface. Consequently, the soil has persisted as the major (99%) pathway of CO₂ release to the atmosphere from the high temperature reservoir at Ohaaki. The CO₂ flux and heat flow surveys indicate that despite 20 a of production the variability in location, spatial extent and magnitude of CO₂ flux remains consistent with established geochemical and geophysical models of the Ohaaki Field. At both OHW and OHE carbon isotopic analyses of soil gas indicate a two-stage fractionation process for moderate-flux (>60 g m⁻² d⁻¹) sites; boiling during fluid ascent within the underlying reservoir and isotopic enrichment as CO₂ diffuses through porous media of the soil zone. For high-flux sites (>300 g m⁻² d⁻¹), the δ¹³CO₂ signature (−7.4 ± 0.3‰ OHW and −6.5 ± 0.6‰ OHE) is unaffected by near-surface (soil zone) fractionation processes and reflects the composition of the boiled magmatic CO₂ source for each respective upflow. Flux thresholds of <30 g m⁻² d⁻¹ for purely diffusive gas transport, between 30 and 300 g m⁻² d⁻¹ for combined diffusive–advective transport, and ?300 g m⁻² d⁻¹ for purely advective gas transport at Ohaaki were assigned. δ¹³CO₂ values and cumulative probability plots of CO₂ flux data both identified a threshold of ∼15 g m⁻² d⁻¹ by which background (atmospheric and soil respired) CO₂ may be differentiated from hydrothermal CO₂.  相似文献   

Influence of glass polymerisation and oxidation on micro-Raman water analysis in alumino-silicate glasses   总被引：1，自引：0，他引：1  

Maxime Mercier  Andrea Di Muro  Daniele Giordano  Priscille Lesne  Bruno Scaillet  Gilles Montagnac 《Geochimica et cosmochimica acta》2009,73(1):197-217

The development of an accurate analytical procedure for determination of dissolved water in complex alumino-silicate glasses via micro-Raman analysis requires the assessment of the spectra topology dependence on glass composition. We report here a detailed study of the respective influence of bulk composition, iron oxidation state and total water content on the absolute and relative intensities of the main Raman bands related to glass network vibrations (LF: ∼490 cm⁻¹; HF: ∼960 cm⁻¹) and total water stretching (H₂O_T: ∼3550 cm⁻¹) in natural glasses. The evolution of spectra topology was examined in (i) 33 anhydrous glasses produced by the re-melting of natural rock samples, which span a very large range of polymerisation degree (NBO/T from 0.00 to 1.16), (ii) 2 sets of synthetic anhydrous basaltic glasses with variable iron oxidation state (Fe³⁺/Fe_T from 0.05 to 0.87), and (iii) 6 sets of natural hydrous glasses (C_H2OT from 0.4 to 7.0 wt%) with NBO/T varying from 0.01 to 0.76.In the explored domain of water concentration, external calibration procedure based on the H₂O_T band height is matrix-independent but its accuracy relies on precise control of the focusing depth and beam energy on the sample. Matrix-dependence strongly affects the internal calibrations based on H₂O_T height scaled to that of LF or HF bands but its effect decreases from acid (low NBO/T, SM) to basic (high NBO/T, SM) glasses. Structural parameters such as NBO/T (non-bridging oxygen per tetrahedron) and SM (sum of structural modifiers) describe the matrix-dependence better than simple compositional parameters (e.g. SiO₂, Na₂O + K₂O). Iron oxidation state has only a minor influence on band topology in basalts and is thus not expected to significantly affect the Raman determinations of water in mafic (e.g. low SiO₂, iron-rich) glasses. Modelling the evolution of the relative band height with polymerisation degree allows us to propose a general equation to predict the dissolved water content in natural glasses: