Investigation of surface water/groundwater relationships in a rural watershed in Southern Honduras

Despite failed attempts at obtaining a potable water system, the village of El Caracol in Southern Honduras remains committed to improving access to water resources. To assist in this endeavor, an investigation of the hydrogeological characteristics of the local watershed was conducted. Daily precipitation was recorded to examine the relationship between precipitation and approximated river and spring discharges. A Thornthwaite Mather Water Balance Model was used to predict monthly discharges for comparison with observed values, and to infer the percentage of topographic watersheds contributing to the respective discharges. As aquifer porosity in this region is thought to be primarily secondary (i.e., fractures), field observed lineaments were compared with those interpreted from remote sensing imagery in an attempt to determine the usefulness of these interpretations in locating potential water sources for a future project.

Groundwater Under the Butte Hill - The East Camp Mining System and the Berkeley Pit

This lecture discusses monitoring activities of the Berkeley Pit for the past 31 years at the Montana Bureau of Mines and Geology in Butte, Montana.

Reading Our Land Use Legacy in Groundwater: Nitrate Dynamics in the Judith River Watershed

Elevated nitrate in groundwater is common is agricultural areas where fertilizer has been added at high rates for decades. Within the Judith River Wastershed, high native soil fertility allowed for dryland wheat production without N fertilization until the 1980s, yet elevated nitrate levels were frequently observed in shallow aquifers. Dr. Stephanie Ewing presents results for soil, groundwater and surface water analyses from a hydrologically isolated strath terrace near Moccasin, MT. In context of this uniquely well constrained field setting, these observed data, along with land use history and a simple mass balance model, revel the long term development and perturbation of native soil fertility with cultivation.

Geochemistry and Stable Isotopes of Surface Water and Groundwater in the Continental Pit in Butte, Montana, USA

The Continental porphyry Cu‐Mo mine, located 2 km east of the famous Berkeley Pit lake of Butte, Montana, contains two small lakes that vary in size depending on mining activity. In contrast to the acidic Berkeley Pit lake, the Continental Pit waters have near-neutral pH and relatively low metal concentrations. The main reason is geological: whereas the Berkeley Pit mined highly‐altered granite rich in pyrite with no neutralizing potential, the Continental Pit is mining weakly‐altered granite with lower pyrite concentrations and up to 1‐2% hydrothermal calcite. The purpose of this study was to gather and interpret information that bears on the chemistry of surface water and groundwater in the active Continental Pit. Pre‐existing chemistry data from sampling of the Continental Pit were compiled from the Montana Bureau of Mines and Geology and Montana Department of Environmental Quality records. In addition, in March of 2013, new water samples were collected from the mine’s main dewatering well, the Sarsfield well, and a nearby acidic seep (Pavilion Seep) and analyzed for trace metals and several stable isotopes, including dD and d18O of water, d13C of dissolved inorganic carbon, and d34S of dissolved sulfate. In December 2013, several soil samples were collected from the shore of the frozen pit lake and surrounding area. The soil samples were analyzed using X‐ray diffraction to determine mineral content. Based on Visual Minteq modeling, water in the Continental Pit lake is near equilibrium with a number of carbonate, sulfate, and molybdate minerals, including calcite, dolomite, rhodochrosite (MnCO3), brochantite (CuSO4·3Cu(OH)2), malachite (Cu2CO3(OH)2), hydrozincite (Zn5(CO3)2(OH)6), gypsum, and powellite (CaMoO4). The fact that these minerals are close to equilibrium suggests that they are present on the weathered mine walls and/or in the sediment of the surface water ponds. X‐Ray Diffraction (XRD) analysis of the pond “beach” sample failed to show any discrete metal‐bearing phases. One of the soil samples collected higher in the mine, near an area of active weathering of chalcocite‐rich ore, contained over 50% chalcanthite (CuSO4·5H2O). This water‐soluble copper salt is easily dissolved in water, and is probably a major source of copper to the pond and underlying groundwater system. However, concentrations of copper in the latter are probably controlled by other, less‐soluble minerals, such as brochantite or malachite. Although the acidity of the Pavilion Seep is high (~ 11 meq/L), the flow is much less than the Sarsfield Well at the current time. Thus, the pH, major and minor element chemistry in the Continental Pit lakes are buffered by calcite and other carbonate minerals. For the Continental Pit waters to become acidic, the influx of acidic seepage (e.g., Pavilion Seep) would need to increase substantially over its present volume.

Tracking Montana’s Groundwater: 20 Years of Lessons Learned

Long-term water-level measurements provide a fundamental indicator of the health of Montana's groundwater resources. For more than 20 years the Groundwater Assessment Program has monitored groundwater levels across the state. This lecture reviews the State's groundwater use and shows how climate variability, groundwater development, and land-use has impacted different aquifers.

Footprint of recycled water subsidies downwind of Lake Michigan

Continental evaporation is a significant and dynamic flux within the atmospheric water budget, but few methods provide robust observational constraints on the large-scale hydroclimatological and hydroecological impacts of this ‘recycled-water' flux. We demonstrate a geospatial analysis that provides such information, using stable isotope data to map the distribution of recycled water in shallow aquifers downwind from Lake Michigan. The δ2H and δ18O values of groundwater in the study region decrease from south to north, as expected based on meridional gradients in climate and precipitation isotope ratios. In contrast, deuterium excess (d = δ2H − 8 × δ18O) values exhibit a significant zonal gradient and finer-scale spatially patterned variation. Local d maxima occur in the northwest and southwest corners of the Lower Peninsula of Michigan, where ‘lake-effect' precipitation events are abundant. We apply a published model that describes the effect of recycling from lakes on atmospheric vapor d values to estimate that up to 32% of recharge into individual aquifers may be derived from recycled Lake Michigan water. Applying the model to geostatistical surfaces representing mean d values, we estimate that between 10% and 18% of the vapor evaporated from Lake Michigan is re-precipitated within downwind areas of the Lake Michigan drainage basin. Our approach provides previously unavailable observational constraints on regional land-atmosphere water fluxes in the Great Lakes Basin and elucidates patterns in recycled-water fluxes that may influence the biogeography of the region. As new instruments and networks facilitate enhanced spatial monitoring of environmental water isotopes, similar analyses can be widely applied to calibrate and validate water cycle models and improve projections of regional hydroecological change involving the coupled lake-atmosphere-land system. Read More: http://www.esajournals.org/doi/abs/10.1890/ES12-00062.1

Tree species and diversity effects on soil water seepage in a tropical plantation

Plant diversity has been shown to influence the water cycle of forest ecosystems by differences in water consumption and the associated effects on groundwater recharge. However, the effects of biodiversity on soil water fluxes remain poorly understood for native tree species plantations in the tropics. Therefore, we estimated soil water fluxes and assessed the effects of tree species and diversity on these fluxes in an experimental native tree species plantation in Sardinilla (Panama). The study was conducted during the wet season 2008 on plots of monocultures and mixtures of three or six tree species. Rainfall and soil water content were measured and evapotranspiration was estimated with the Penman-Monteith equation. Soil water fluxes were estimated using a simple soil water budget model considering water input, output, and soil water and groundwater storage changes and in addition, were simulated using the physically based one-dimensional water flow model Hydrus-1D. In general, the Hydrus simulation did not reflect the observed pressure heads, in that modeled pressure heads were higher compared to measured ones. On the other hand, the results of the water balance equation (WBE) reproduced observed water use patterns well. In monocultures, the downward fluxes through the 200 cm-depth plane were highest below Hura crepitans (6.13 mm day−1) and lowest below Luehea seemannii (5.18 mm day−1). The average seepage rate in monocultures (±SE) was 5.66 ± 0.18 mm day−1, and therefore, significantly higher than below six-species mixtures (5.49 ± 0.04 mm day−1) according to overyielding analyses. The three-species mixtures had an average seepage rate of 5.63 ± 0.12 mm day−1 and their values did not differ significantly from the average values of the corresponding species in monocultures. Seepage rates were driven by the transpiration of the varying biomass among the plots (r = 0.61, p = 0.017). Thus, a mixture of trees with different growth rates resulted in moderate seepage rates compared to monocultures of either fast growing or slow growing tree species. Our results demonstrate that tree-species specific biomass production and tree diversity are important controls of seepage rates in the Sardinilla plantation during the wet season.

Pelagic molybdenum concentration anomalies and the impact of sediment resuspension on the molybdenum budget in two tidal systems of the North Sea

The seasonal dynamics of molybdenum (Mo) were studied in the water column of two tidal basins of the German Wadden Sea (Sylt-Rømø and Spiekeroog) between 2007 and 2011. In contrast to its conservative behaviour in the open ocean, both, losses of more than 50% of the usual concentration level of Mo in seawater and enrichments up to 20% were observed repeatedly in the water column of the study areas. During early summer, Mo removal by adsorption on algae-derived organic matter (e.g. after Phaeocystis blooms) is postulated to be a possible mechanism. Mo bound to organic aggregates is likely transferred to the surface sediment where microbial decomposition enriches Mo in the pore water. First δ98/95Mo data of the study area disclose residual Mo in the open water column being isotopically heavier than MOMo (Mean Ocean Molybdenum) during a negative Mo concentration anomaly, whereas suspended particulate matter shows distinctly lighter values. Based on field observations a Mo isotope enrichment factor of ε = −0.3‰ has been determined which was used to argue against sorption on metal oxide surfaces. It is suggested here that isotope fractionation is caused by biological activity and association to organic matter. Pelagic Mo concentration anomalies exceeding the theoretical salinity-based concentration level, on the other hand, cannot be explained by replenishment via North Sea waters alone and require a supply of excess Mo. Laboratory experiments with natural anoxic tidal flat sediments and modelled sediment displacement during storm events suggest fast and effective Mo release during the resuspension of anoxic sediments in oxic seawater as an important process for a recycling of sedimentary sulphide bound Mo into the water column.

Groundwater age distributions at a public drinking water supply well field derived from multiple age tracers (⁸⁵Kr, ³H/³He, and ³⁹Ar)

Groundwater age is a key aspect of production well vulnerability. Public drinking water supply wells typically have long screens and are expected to produce a mixture of groundwater ages. The groundwater age distributions of seven production wells of the Holten well field (Netherlands) were estimated from tritium-helium (3H/3He), krypton-85 (85Kr), and argon-39 (39Ar), using a new application of a discrete age distribution model and existing mathematical models, by minimizing the uncertainty-weighted squared differences of modeled and measured tracer concentrations. The observed tracer concentrations fitted well to a 4-bin discrete age distribution model or a dispersion model with a fraction of old groundwater. Our results show that more than 75 of the water pumped by four shallow production wells has a groundwater age of less than 20 years and these wells are very vulnerable to recent surface contamination. More than 50 of the water pumped by three deep production wells is older than 60 years. 3H/3He samples from short screened monitoring wells surrounding the well field constrained the age stratification in the aquifer. The discrepancy between the age stratification with depth and the groundwater age distribution of the production wells showed that the well field preferentially pumps from the shallow part of the aquifer. The discrete groundwater age distribution model appears to be a suitable approach in settings where the shape of the age distribution cannot be assumed to follow a simple mathematical model, such as a production well field where wells compete for capture area.

Dating groundwater in the Bohemian Cretaceous Basin: Understanding tracer variations in the subsurface

The northern section of the Bohemian Cretaceous Basin has been the site of intensive U exploitation with harmful impacts on groundwater quality. The understanding of groundwater flow and age distribution is crucial for the prediction of the future dispersion and impact of the contamination. State of the art tracer methods (3H, 3He, 4He, 85Kr, 39Ar and 14C) were, therefore, used to obtain insights to ageing and mixing processes of groundwater along a north–south flow line in the centre of the two most important aquifers of Cenomanian and middle Turonian age. Dating of groundwater is particularly complex in this area as: (i) groundwater in the Cenomanian aquifer is locally affected by fluxes of geogenic and biogenic gases (e.g. CO2, CH4, He) and by fossil brines in basement rocks rich in Cl and SO4; (ii) a thick unsaturated zone overlays the Turonian aquifer; (iii) a periglacial climate and permafrost conditions prevailed during the Last Glacial Maximum (LGM), and iv) the wells are mostly screened over large depth intervals. Large disagreements in 85Kr and 3H/3He ages indicate that processes other than ageing have affected the tracer data in the Turonian aquifer. Mixing with older waters (>50 a) was confirmed by 39Ar activities. An inverse modelling approach, which included time lags for tracer transport throughout the unsaturated zone and degassing of 3He, was used to estimate the age of groundwater. Best fits between model and field results were obtained for mean residence times varying from modern up to a few hundred years. The presence of modern water in this aquifer is correlated with the occurrence of elevated pollution (e.g. nitrates). An increase of reactive geochemical indicators (e.g. Na) and radiogenic 4He, and a decrease in 14C along the flow direction confirmed groundwater ageing in the deeper confined Cenomanian aquifer. Radiocarbon ages varied from a few hundred years to more than 20 ka. Initial 14C activity for radiocarbon dating was calibrated by means of 39Ar measurements. The 14C age of a sample recharged during the LGM was further confirmed by depleted stable isotope signatures and near freezing point noble gas temperature. Radiogenic 4He accumulated in groundwater with concentrations increasing linearly with 14C ages. This enabled the use of 4He to validate the dating range of 14C and extend it to other parts of this aquifer. In the proximity of faults, 39Ar in excess of modern concentrations and 14C dead CO2 sources, elevated 3He/4He ratios and volcanic activity in Oligocene to Quaternary demonstrate the influence of gas of deeper origin and impeded the application of 4He, 39Ar and 14C for groundwater dating.

The global carbon budget 1959–2011

Accurate assessments of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the climate policy process, and project future climate change. Present-day analysis requires the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. Here we describe datasets and a methodology developed by the global carbon cycle science community to quantify all major components of the global carbon budget, including their uncertainties. We discuss changes compared to previous estimates, consistency within and among components, and methodology and data limitations. CO2 emissions from fossil fuel combustion and cement production (EFF) are based on energy statistics, while emissions from Land-Use Change (ELUC), including deforestation, are based on combined evidence from land cover change data, fire activity in regions undergoing deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. Finally, the global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms. For the last decade available (2002–2011), EFF was 8.3 ± 0.4 PgC yr−1, ELUC 1.0 ± 0.5 PgC yr−1, GATM 4.3 ± 0.1 PgC yr−1, SOCEAN 2.5 ± 0.5 PgC yr−1, and SLAND 2.6 ± 0.8 PgC yr−1. For year 2011 alone, EFF was 9.5 ± 0.5 PgC yr−1, 3.0 percent above 2010, reflecting a continued trend in these emissions; ELUC was 0.9 ± 0.5 PgC yr−1, approximately constant throughout the decade; GATM was 3.6 ± 0.2 PgC yr−1, SOCEAN was 2.7 ± 0.5 PgC yr−1, and SLAND was 4.1 ± 0.9 PgC yr−1. GATM was low in 2011 compared to the 2002–2011 average because of a high uptake by the land probably in response to natural climate variability associated to La Niña conditions in the Pacific Ocean. The global atmospheric CO2 concentration reached 391.31 ± 0.13 ppm at the end of year 2011. We estimate that EFF will have increased by 2.6% (1.9–3.5%) in 2012 based on projections of gross world product and recent changes in the carbon intensity of the economy. All uncertainties are reported as ±1 sigma (68% confidence assuming Gaussian error distributions that the real value lies within the given interval), reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. This paper is intended to provide a baseline to keep track of annual carbon budgets in the future.