Various thermodynamic properties of H2O that are defined as pressure or temperature derivatives of some other variable, such as isothermal compressibility (β, pressure derivative of density), ...isobaric thermal expansion (α, temperature derivative of density), and specific isobaric heat capacity (cf, temperature derivative of enthalpy), all show large magnitudes near the critical point, reflecting large variations in fluid density and specific enthalpy with small changes in temperature and pressure. As a result, mass (related to fluid density) and energy (related to fluid enthalpy) transport in this PT region are sensitive to changing PT conditions. Addition of NaCl to H2O causes the region of anomalous behavior, here defined as the critical region, to migrate to higher temperatures and pressures. The critical region is defined as that region of PT space in which the dimensionless reduced susceptibility χ~ ≥ 0.5. When NaCl is added to H2O, the critical region migrates to higher temperature and pressure. However, the absolute magnitudes of thermodynamic properties that are defined as temperature and/or pressure derivatives (α, β, and cf) all decrease with increasing salinity. Thus, the mass and energy transporting capacities of hydrothermal fluids in the critical region become less sensitive to changing PT conditions as the salinity increases. For example, quartz solubility can be described as a function of fluid density, and because density becomes less sensitive to changing PT conditions as salinity increases, quartz solubility also becomes less sensitive to changing PT conditions as fluid salinity increases. Similarly, fluxibility describes the ability of a fluid to transport heat by buoyancy‐driven convection, and fluxibility decreases with increasing salinity. Results of this study show that the mass and energy transport capacity of fluids in the Earth's crust are maximized in the critical region and that the sensitivity to changing PT conditions decreases with increasing salinity.
Pure H2O exhibits anomalous behavior in the vicinity of the critical point, reflected by large variations in density and specific enthalpy with small changes in temperature and pressure. Mass and energy transport properties that are temperature or pressure derivatives of density and specific enthalpy thus show large variability near the critical point. Addition of NaCl to H2O causes the critical region in which fluid properties exhibit anomalous behavior to migrate to higher temperatures and pressures.
Melt inclusions (MI) are considered the best tool available for determining the pre-eruptive volatile contents of magmas. H2O and CO2 concentrations of the glass phase in MI are commonly used both as ...a barometer and to track magma degassing behavior during ascent due to the strong pressure dependence of H2O and CO2 solubilities in silicate melts. The often unstated and sometimes overlooked requirement for this method to be valid is that the glass phase in the MI must represent the composition of the melt that was trapped at depth in the volcanic plumbing system. However, melt inclusions commonly contain a vapor bubble that formed after trapping owing to differential shrinkage of the melt compared to the host crystal, and/or crystallization at the inclusion-host interface. Such bubbles may contain a substantial portion of volatiles, such as CO2, that were originally dissolved in the melt. In this study, we determined the contribution of CO2 in the vapor bubble to the overall CO2 content of MI based on quantitative Raman analysis of the vapor bubbles in MI from the 1959 Kilauea Iki (Hawaii), 1960 Kapoho (Hawaii), 1974 Fuego volcano (Guatemala), and 1977 Seguam Island (Alaska) eruptions. We found that the bubbles typically contain 40 to 90% of the total CO2 in the MI. Reconstructing the original CO2 content by adding the CO2 in the bubble back into the melt results in an increase in CO2 concentration by as much as an order of magnitude (thousands of parts per million). Reconstructed CO2 concentrations correspond to trapping pressures that are significantly greater than one would predict based on analysis of the volatiles in the glass alone. Trapping depths can be as much as 10 km deeper than estimates that ignore the CO2 in the bubble. In addition to CO2 in the vapor bubbles, many MI showed the presence of a carbonate mineral phase. Failure to recognize the carbonate during petrographic examination or analysis of the glass and to include its contained CO2 when reconstructing the CO2 content of the originally trapped melt will introduce additional errors into the calculated volatile budget. Our results emphasize that accurate determination of the pre-eruptive volatile content of melts based on analysis of melt inclusions must consider the volatiles contained in the bubble (and carbonates, if present). This can be accomplished either by analysis of the bubble and the glass followed by mass-balance reconstruction of the original volatile content of the melt, or by re-homogenization of the MI prior to conducting microanalysis of the quenched, glassy MI.
Melt inclusions (MI) represent the best source of information concerning the pre-eruptive volatile contents of magmas. If the trapped melt is enriched in volatile species, following trapping the MI ...may generate a vapor bubble containing volatiles that have exsolved from the melt. Thermodynamic modeling of vapor-saturated albitic composition (NaAlSi3O8) MI shows that the CO2 content of the melt phase in the MI is sensitive to small amounts of post-entrapment crystallization (PEC), whereas the H2O content of the melt is less sensitive to PEC. During PEC, CO2 is transferred from the melt to the vapor phase and the vapor bubble may contain a significant amount, if not most, of the CO2 in the MI. The contrasting behaviors of H2O and CO2 during PEC lead to H2O-CO2 trends that are similar to those predicted for open-system degassing during magma ascent and decompression. Thus, similar H2O-CO2 trends may be produced if (1) vapor-saturated MI are trapped at various depths along a magmatic ascent path, or (2) MI having the same volatile content are all trapped at the same depth, but undergo different amounts of PEC following trapping. It is not possible to distinguish between these two contrasting interpretations based on MI volatile data alone. However, by examining the volatile trends within the context of other geochemical monitors of crystallization or magma evolution progress, it may be possible to determine whether the volatile trends were generated along a degassing path or if they reflect various amounts of PEC in an originally homogeneous melt inclusion assemblage. The volatile trends resulting from PEC of MI described in this study are directly applicable to silica-rich (granitic) MI trapped in non-ferromagnesian host phases, and are only qualitatively applicable to more mafic melt compositions and/or host phases owing to modifications resulting from Fe exchange with the host and to post-entrapment re-equilibration processes.
The 1959 Kilauea Iki eruption provides a unique opportunity to investigate the process of shallow magma mixing, its impact on the magmatic volatile budget and its role in triggering and driving ...episodes of Hawaiian fountaining. Melt inclusions hosted by olivine record a continuous decrease in H2O concentration through the 17 episodes of the eruption, while CO2 concentrations correlate with the degree of post-entrapment crystallization of olivine on the inclusion walls. Geochemical data, when combined with the magma budget and with contemporaneous eruption observations, show complex mixing between episodes involving hot, geochemically heterogeneous melts from depth, likely carrying exsolved vapor, and melts which had erupted at the surface, degassed and drained-back into the vent. The drained-back melts acted as a coolant, inducing rapid cooling of the more primitive melts and their olivines at shallow depths and inducing crystallization and vesiculation and triggering renewed fountaining. A consequence of the mixing is that the melts became vapor-undersaturated, so equilibration pressures cannot be inferred from them using saturation models. After the melt inclusions were trapped, continued growth of vapor bubbles, caused by enhanced post-entrapment crystallization, sequestered a large fraction of CO2 from the melt within the inclusions. This study, while cautioning against accepting melt inclusion CO2 concentrations "as measured" in mixed magmas, also illustrates that careful analysis and interpretation of post-entrapment modifications can turn this apparent challenge into a way to yield novel useful insights into the geochemical controls on eruption intensity.
The shifts in wavenumber of the ν3(SiO4) (approximately 1008 cm-1) Raman band of fully crystalline synthetic zircon with changing pressure (P) and temperature (T) were calibrated for application as a ...Raman spectroscopic pressure sensor in optical cells to about 1000 °C and 10 GPa. The relationship between wavenumber (ν) of this band and T from 22 to 950 °C is described by the equation ν (cm-1) = 7.54·10-9·T3 - 1.61·10-5·T2 - 2.89·10-2·T + 1008.9, where T is given in °C. The pressure dependence is nearly linear over the studied range in P. At approximately 25 °C, the θν/θP slope to 6.6 GPa is 5.69 cm-1/GPa, and that to 2 GPa is 5.77 cm-1/GPa. The θν/θP slope does not significantly change with temperature, as determined from experiments conducted along isotherms up to 700 °C. Therefore, this pressure sensor has the advantage that a constant θν/θP slope of 5.8 ± 0.1 cm-1/GPa can be applied in experiments to pressures of at least about 6.6 GPa without introducing a significant error. The pressure sensor was tested to determine isochores in experiments with H2O+Na2Si3O7 and H2O+NaAlSi3O8 fluids to 803 °C and 1.65 GPa. These pressures were compared to pressures calculated from the equation of state (EoS) of H2O based on the measured vapor dissolution or ice melting temperature for the same experiment. Pressures determined from the zircon sensor in runs in which NaAlSi3O8 melt dissolved in aqueous fluid were close to or lower than the pressure calculated from the EoS of H2O using the vapor dissolution or ice melting temperature. In experiments with H2O+Na2O+SiO2 fluids, however, the pressure obtained from the Raman spectrum of zircon was often significantly higher than that estimated from the EoS of H2O. This suggests that the pressures along some critical curves of water-silicate melt pseudobinary systems should be revised.
Various thermodynamic properties of H sub(2)O that are defined as pressure or temperature derivatives of some other variable, such as isothermal compressibility ( beta , pressure derivative of ...density), isobaric thermal expansion ( alpha , temperature derivative of density), and specific isobaric heat capacity (cf, temperature derivative of enthalpy), all show large magnitudes near the critical point, reflecting large variations in fluid density and specific enthalpy with small changes in temperature and pressure. As a result, mass (related to fluid density) and energy (related to fluid enthalpy) transport in this PT region are sensitive to changing PT conditions. Addition of NaCl to H sub(2)O causes the region of anomalous behavior, here defined as the critical region, to migrate to higher temperatures and pressures. The critical region is defined as that region of PT space in which the dimensionless reduced susceptibility.
Abstract
Various thermodynamic properties of H
2
O that are defined as pressure or temperature derivatives of some other variable, such as isothermal compressibility (β, pressure derivative of ...density), isobaric thermal expansion (α, temperature derivative of density), and specific isobaric heat capacity (
c
f
, temperature derivative of enthalpy), all show large magnitudes near the critical point, reflecting large variations in fluid density and specific enthalpy with small changes in temperature and pressure. As a result, mass (related to fluid density) and energy (related to fluid enthalpy) transport in this
PT
region are sensitive to changing
PT
conditions. Addition of NaCl to H
2
O causes the region of anomalous behavior, here defined as the critical region, to migrate to higher temperatures and pressures. The critical region is defined as that region of
PT
space in which the dimensionless reduced susceptibility
≥ 0.5. When NaCl is added to H
2
O, the critical region migrates to higher temperature and pressure. However, the absolute magnitudes of thermodynamic properties that are defined as temperature and/or pressure derivatives (α, β, and
c
f
) all decrease with increasing salinity. Thus, the mass and energy transporting capacities of hydrothermal fluids in the critical region become less sensitive to changing
PT
conditions as the salinity increases. For example, quartz solubility can be described as a function of fluid density, and because density becomes less sensitive to changing
PT
conditions as salinity increases, quartz solubility also becomes less sensitive to changing
PT
conditions as fluid salinity increases. Similarly, fluxibility describes the ability of a fluid to transport heat by buoyancy‐driven convection, and fluxibility decreases with increasing salinity. Results of this study show that the mass and energy transport capacity of fluids in the Earth's crust are maximized in the critical region and that the sensitivity to changing
PT
conditions decreases with increasing salinity.