Volatile emissions and gas geochemistry of Hot Spring Basin,Yellowstone National Park,USA

We characterize and quantify volatile emissions at Hot Spring Basin (HSB), a large acid-sulfate region that lies just outside the northeastern edge of the 640 ka Yellowstone Caldera. Relative to other thermal areas in Yellowstone, HSB gases are rich in He and H₂, and mildly enriched in CH₄ and H₂S. Gas compositions are consistent with boiling directly off a deep geothermal liquid at depth as it migrates toward the surface. This fluid, and the gases evolved from it, carries geochemical signatures of magmatic volatiles and water–rock reactions with multiple crustal sources, including limestones or quartz-rich sediments with low K/U (or ^40?Ar/^4?He). Variations in gas chemistry across the region reflect reservoir heterogeneity and variable degrees of boiling. Gas-geothermometer temperatures approach 300 °C and suggest that the reservoir feeding HSB is one of the hottest at Yellowstone. Diffuse CO₂ flux in the western basin of HSB, as measured by accumulation-chamber methods, is similar in magnitude to other acid-sulfate areas of Yellowstone and is well correlated to shallow soil temperatures. The extrapolation of diffuse CO₂ fluxes across all the thermal/altered area suggests that 410 ± 140 t d^− 1 CO₂ are emitted at HSB (vent emissions not included). Diffuse fluxes of H₂S were measured in Yellowstone for the first time and likely exceed 2.4 t d^− 1 at HSB. Comparing estimates of the total estimated diffuse H₂S emission to the amount of sulfur as SO₄²⁻ in streams indicates ~ 50% of the original H₂S in the gas emission is lost into shallow groundwater, precipitated as native sulfur, or vented through fumaroles. We estimate the heat output of HSB as ~ 140–370 MW using CO₂ as a tracer for steam condensate, but not including the contribution from fumaroles and hydrothermal vents. Overall, the diffuse heat and volatile fluxes of HSB are as great as some active volcanoes, but they are a small fraction (1–3% for CO₂, 2–8% for heat) of that estimated for the entire Yellowstone system.     

A Hybrid Finite-Difference and Analytic Element Groundwater Model

Regional finite-difference models tend to have large cell sizes, often on the order of 1–2 km on a side. Although the regional flow patterns in deeper formations may be adequately represented by such a model, the intricate surface water and groundwater interactions in the shallower layers are not. Several stream reaches and nearby wells may occur in a single cell, precluding any meaningful modeling of the surface water and groundwater interactions between the individual features. We propose to replace the upper MODFLOW layer or layers, in which the surface water and groundwater interactions occur, by an analytic element model (GFLOW) that does not employ a model grid; instead, it represents wells and surface waters directly by the use of point-sinks and line-sinks. For many practical cases it suffices to provide GFLOW with the vertical leakage rates calculated in the original coarse MODFLOW model in order to obtain a good representation of surface water and groundwater interactions. However, when the combined transmissivities in the deeper (MODFLOW) layers dominate, the accuracy of the GFLOW solution diminishes. For those cases, an iterative coupling procedure, whereby the leakages between the GFLOW and MODFLOW model are updated, appreciably improves the overall solution, albeit at considerable computational cost. The coupled GFLOW–MODFLOW model is applicable to relatively large areas, in many cases to the entire model domain, thus forming an attractive alternative to local grid refinement or inset models.     

Estimating recharge rates with analytic element models and parameter estimation

Quantifying the spatial and temporal distribution of recharge is usually a prerequisite for effective ground water flow modeling. In this study, an analytic element (AE) code (GFLOW) was used with a nonlinear parameter estimation code (UCODE) to quantify the spatial and temporal distribution of recharge using measured base flows as calibration targets. The ease and flexibility of AE model construction and evaluation make this approach well suited for recharge estimation. An AE flow model of an undeveloped watershed in northern Wisconsin was optimized to match median annual base flows at four stream gages for 1996 to 2000 to demonstrate the approach. Initial optimizations that assumed a constant distributed recharge rate provided good matches (within 5%) to most of the annual base flow estimates, but discrepancies of >12% at certain gages suggested that a single value of recharge for the entire watershed is inappropriate. Subsequent optimizations that allowed for spatially distributed recharge zones based on the distribution of vegetation types improved the fit and confirmed that vegetation can influence spatial recharge variability in this watershed. Temporally, the annual recharge values varied >2.5-fold between 1996 and 2000 during which there was an observed 1.7-fold difference in annual precipitation, underscoring the influence of nonclimatic factors on interannual recharge variability for regional flow modeling. The final recharge values compared favorably with more labor-intensive field measurements of recharge and results from studies, supporting the utility of using linked AE-parameter estimation codes for recharge estimation.     

Are toxic contaminants accumulating in Massachusetts coastal sediments following startup of the Massachusetts Bay outfall: a comprehensive comparison of baseline and post-diversion periods

The Massachusetts Water Resources Authority (MWRA) conducts monitoring to address concerns related to the 2000 diversion of secondarily treated effluent discharge into Massachusetts Bay. Baseline data (1992-2000) showed multiple regions defined by physical and chemical composition. Near the Massachusetts Bay outfall, there is a series of heterogeneous sediments in relatively close proximity to the primary historic source of contaminants (Boston Harbor). Farfield sediments exhibited greater compositional definition from one another, which was attributed to the greater spatial separation of the sampling locations. Factors that influence contaminant variability include local and distributed sources, and are primarily related to gradients in depositional environments. Post-diversion sediment data suggest that 4 years of treated effluent discharge has not increased contaminant concentrations to the bay system. However, abundance of the sewage tracer, Clostridium perfringens, has increased variably in sediments located within 2 km of the outfall, providing a distinct effluent signal near the outfall.     

On Constraining Pilot Point Calibration with Regularization in PEST

Ground water model calibration has made great advances in recent years with practical tools such as PEST being instrumental for making the latest techniques available to practitioners. As models and calibration tools get more sophisticated, however, the power of these tools can be misapplied, resulting in poor parameter estimates and/or nonoptimally calibrated models that do not suit their intended purpose. Here, we focus on an increasingly common technique for calibrating highly parameterized numerical models—pilot point parameterization with Tikhonov regularization. Pilot points are a popular method for spatially parameterizing complex hydrogeologic systems; however, additional flexibility offered by pilot points can become problematic if not constrained by Tikhonov regularization. The objective of this work is to explain and illustrate the specific roles played by control variables in the PEST software for Tikhonov regularization applied to pilot points. A recent study encountered difficulties implementing this approach, but through examination of that analysis, insight into underlying sources of potential misapplication can be gained and some guidelines for overcoming them developed.     

Comparison of TL multiple aliquot, single grain GLSL SAR and C-14 ages for the Puritjarra, Australia, rock shelter

We report single-grain quartz luminescence ages for the Puritjarra rock shelter, with the aim of resolving an apparent discrepancy between ages obtained by ¹⁴C and a variety of luminescence methods, previously reported. Ages now found at all depths to 75 cm (ages to 30 ka) can be interpreted as largely resolving the differences. This implied caveat arises because single-grain methods are statistically inefficient. As a consequence, a degree of interpretation is inevitable in analysing the data. The emphasis in the present paper is an analysis making use of weighted histograms. The measurements by single grain OSL and ¹⁴C, including ABOX–SC, taken together, can be regarded as compatible. They indicate human occupation of the Puritjarra rock shelter at least as early as 30 ka BP.