Where you drink water: An assessment of the Tennessee, USA public water supply.

Monitoring drinking water quality is essential to protect people's health and wellbeing. In the United States, the Safe Drinking Water Information System (SDWIS) database records the occurrence of a drinking water violation regulation in public water systems. A notable shortcoming of SDWIS is the lack of the contaminant concentration level about the allowable maximum contaminant threshold. In this study, we take advantage of both the SDWIS violation database and the contaminants sampling database at the state level to examine the drinking water quality of all kinds of drinking water systems in detail. We obtained sampling data (i.e., the concentration level of contaminants) of public water systems (PWSs) in Tennessee and explored the statistical distribution of contaminant concentration data in relation to the enforceable maximum regulatory contaminant level). We use both SDWIS violation records and actual concentrations of contaminants from the sampling data to study the factors that influence the drinking water quality of PWSs. We find that different types of violations were more frequent in (1) specific geological regions, (2) counties with PWSs that serve a larger population (10,000 to 100,000 people), and (3) places with abundant surface water, such as near a lake or major river. Additionally, the distribution of measured concentrations for many contaminants was not smooth but was punctuated by discontinuities at selected levels, such as at 50% of the maximum contaminant level. Such anomalies in the sampling data do not indicate violations, but more investigation is needed to determine the reasons behind the punctuated changes.

A Search for Freshwater in the Saline Aquifer of Coastal Bangladesh.

In the polder region of coastal Bangladesh, shallow groundwater is primarily brackish with unpredictable occurrence of freshwater pockets. Delta building processes, including the codeposition of fresh-to-saline porewater and sediments, have formed the shallow aquifer. Impermeable clay facies and the lack of a topographical gradient limit the flow of groundwater and its mixing with surface water so controls on spatial variability of salinity are not obvious. By characterizing groundwater-surface water (GW-SW) interactions, this study attempted to identify areas of potable groundwater for the polder communities. We used transects of piezometers, cores, electromagnetic induction, and water chemistry surveys to explore two sources of potential fresh groundwater: (1) tidal channel-aquifer exchange and (2) meteoric recharge. Fresh groundwater proved difficult to find due to heterogeneous subsurface lithology, asymmetrical tidal dynamics, extreme seasonal fluctuations in rainfall, and limited field data. Geophysical observations suggest substantial lateral variability in shallow subsurface conductivity profiles. Piezometers show varying degrees of tidal pressure attenuation away from the channels. Nevertheless, the active exchange of freshwater appears to be limited due to low permeability of banks and surface sediments. Results indicate that pockets of fresh groundwater cannot be identified using readily available hydrogeological methods, so alternative drinking water sources should be pursued. By better understanding the hydrogeology of the system, however, communities will be better equipped to redirect water management resources to more feasible and sustainable drinking water options.

Minimizing irrigation water demand: An evaluation of shifting planting dates in Sri Lanka.

Climate change coupled with increasing demands for water necessitates an improved understanding of the water-food nexus at a scale local enough to inform farmer adaptations. Such assessments are particularly important for nations with significant small-scale farming and high spatial variability in climate, such as Sri Lanka. By comparing historical patterns of irrigation water requirements (IWRs) to rice planting records, we estimate that shifting rice planting dates to earlier in the season could yield water savings of up to 6%. Our findings demonstrate the potential of low-cost adaptation strategies to help meet crop production demands in water-scarce environments. This local-scale assessment of IWRs in Sri Lanka highlights the value of using historical data to inform agricultural management of water resources when high-skilled forecasts are not available. Given national policies prioritizing in-country production and farmers' sensitivities to water stress, decision-makers should consider local degrees of climate variability in institutional design of irrigation management structures.

Sources of salinity and arsenic in groundwater in southwest Bangladesh.

BACKGROUND: High salinity and arsenic (As) concentrations in groundwater are widespread problems in the tidal deltaplain of southwest Bangladesh. To identify the sources of dissolved salts and As, groundwater samples from the regional shallow Holocene aquifer were collected from tubewells during the dry (May) and wet (October) seasons in 2012-2013. Thirteen drill cores were logged and 27 radiocarbon ages measured on wood fragments to characterize subsurface stratigraphy. RESULTS: Drill cuttings, exposures in pits and regional studies reveal a >5 m thick surface mud cap overlying a ~30 m thick upper unit of interbedded mud and fine sand layers, and a coarser lower unit up to 60 m thick dominated by clean sands, all with significant horizontal variation in bed continuity and thickness. This thick lower unit accreted at rates of ~2 cm/year through the early Holocene, with local subsidence or compaction rates of 1-3 mm/year. Most tubewells are screened at depths of 15-52 m in sediments deposited 8000-9000 YBP. Compositions of groundwater samples from tubewells show high spatial variability, suggesting limited mixing and low and spatially variable recharge rates and flow velocities. Groundwaters are Na-Cl type and predominantly sulfate-reducing, with specific conductivity (SpC) from 3 to 29 mS/cm, high dissolved organic carbon (DOC) 11-57 mg/L and As 2-258 ug/L, and low sulfur (S) 2-33 mg/L. CONCLUSIONS: Groundwater compositions can be explained by burial of tidal channel water and subsequent reaction with dissolved organic matter, resulting in anoxia, hydrous ferric oxide (HFO) reduction, As mobilization, and sulfate (SO4) reduction and removal in the shallow aquifer. Introduction of labile organic carbon in the wet season as rice paddy fertilizer may also cause HFO reduction and As mobilization. Variable modern recharge occurred in areas where the clay cap pinches out or is breached by tidal channels, which would explain previously measured (14)C groundwater ages being less than depositional ages. Of samples collected from the shallow aquifer, Bangladesh Government guidelines are exceeded in 46 % for As and 100 % for salinity.

Impacts of management alternatives on rice yield and nitrogen losses to the environment: A case study in rural Sri Lanka.

Maintaining crop yields is vital as populations increase, but environmental degradation resulting from cultivation must be prevented. In particular, freshwater resources are at risk of nitrate leaching from superfluous fertilization. This research explores the tradeoffs between maximizing yield and limiting environmental impacts of rice production in Sri Lanka. The DeNitrification-DeComposition (DNDC) model was used to examine how various combinations of fertilization and irrigation management affect yield, nitrous oxide (N2O) emissions, and nitrogen (N) leaching in paddy systems under climate and soil conditions in the dry zone of Sri Lanka from 1991 to 2010. Simulated fertilizer application rates ranged from zero to 2700 kgN/ha and simulated irrigation schemes were continuously flooded, marginally flooded, and rain-fed. Increasing fertilizer levels from zero to 300 kgN/ha per year increased yield but application of fertilizer beyond that amount ceased to affect yield for any of the three irrigation schemes. The combination of management options for obtaining the maximum grain yield, near 9000 kgC/ha, with the greatest amount of N uptake and relatively low nitrate leaching was using 225 kgN/ha under a continuously flooded regime. This research explores how cultivation in rice-growing regions in south Asia affects the environment and the N cycle, and demonstrates how informed management of these systems can reduce external inputs of N fertilizer without impacting yield.

The effects of household management practices on the global warming potential of urban lawns.

Nitrous oxide (N2O) emissions are an important component of the greenhouse gas (GHG) budget for urban turfgrasses. A biogeochemical model DNDC successfully captured the magnitudes and patterns of N2O emissions observed at an urban turfgrass system at the Richland Creek Watershed in Nashville, TN. The model was then used to study the long-term (i.e. 75 years) impacts of lawn management practice (LMP) on soil organic carbon sequestration rate (dSOC), soil N2O emissions, and net Global Warming Potentials (net GWPs). The model simulated N2O emissions and net GWP from the three management intensity levels over 75 years ranged from 0.75 to 3.57 kg N ha(-1)yr(-1) and 697 to 2443 kg CO2-eq ha(-1)yr(-1), respectively, which suggested that turfgrasses act as a net carbon emitter. Reduction of fertilization is most effective to mitigate the global warming potentials of turfgrasses. Compared to the baseline scenario, halving fertilization rate and clipping recycle as an alternative to synthetic fertilizer can reduce net GWPs by 17% and 12%, respectively. In addition, reducing irrigation and mowing are also effective in lowering net GWPs. The minimum-maintenance LMP without irrigation and fertilization can reduce annual N2O emissions and net GWPs by approximately 53% and 70%, respectively, with the price of gradual depletion of soil organic carbon, when compared to the intensive-maintenance LMP. A lawn age-dependent best management practice is recommended: a high dose fertilizer input at the initial stage of lawn establishment to enhance SOC sequestration, followed by decreasing fertilization rate when the lawn ages to minimize N2O emissions. A minimum-maintained LMP with clipping recycling, and minimum irrigation and mowing, is recommended to mitigate global warming effects from urban turfgrass systems. Among all practices, clipping recycle may be a relatively malleable behavior and, therefore, a good target for interventions seeking to reduce the environmental impacts of lawn management through public education. Our results suggest that a long-term or a chronosequence study of turfgrasses with varying ages is warranted to capture the complete dynamics of contribution of turfgrasses to global warming.

Drought planning and management: using high spatial resolution as part of the solution.

Water scarcity is intensified by drought, a phenomenon that impacts many sectors of society and affects virtually all climate zones. The Palmer drought indices are widely used by scientists and policy makers to understand drought and model its components. Despite the spatial heterogeneity and variability in variables required by the Palmer model, regional index values are most commonly used for real-time drought assessment. Local stakeholders charged with developing flexible and tailored water management policies have articulated the need for drought indices calculated at finer spatial resolutions than a regional scale. We use the Pacific Northwest United States (U.S.) as a study area to demonstrate the differences between drought indices calculated for U.S. climate divisions with those calculated at a 0.5° by 0.5° latitude/longitude resolution. Our results indicate that regional values of the two cumulative Palmer drought indices do not represent finer-resolution values well. For half of the study area, the pictures of drought (as determined by regional and finer-resolution values) are drastically different more than 30% of the time. Thus, quite often water managers do not have a clear understanding of the relative severity of drought in their area, which can have serious implications for drought mitigation and adaptation.

Gaining perspective on the water-energy nexus at the community scale.

Water and energy resources are interrelated but their influence on each other is rarely considered. To quantify the water and energy portfolios associated with a community's water-energy nexus (WEN) and the influence of geographic location on resources, we present the WEN tool. The WEN tool quantifies a community's transport (consumed for or lost before delivery) and nexus (energy for water and water for energy) resources so communities can assess their resource flows. In addition, to provide insight into the full range of impacts of water and energy resource acquisition and to frame the influence of geography on resources, we coin the term "urban resource islands". The concept of urban resource islands provides a framework for considering the implication of geography on a community's water and energy resource acquisition and use. The WEN tool and the concept of resource islands can promote communities to think about their hidden resources and integrate such concepts into their sustainability trade-off analyses and policy decisions. In this paper, we use Tucson, Arizona, United States as a case study.

Modeling effects of multinode wells on solute transport.

Long-screen wells or long open boreholes with intraborehole flow potentially provide pathways for contaminants to move from one location to another in a ground water flow system. Such wells also can perturb a flow field so that the well will not provide water samples that are representative of ground water quality a short distance away from the well. A methodology is presented to accurately and efficiently simulate solute transport in ground water systems that include wells longer than the grid spacing used in a simulation model of the system and hence are connected to multiple nodes of the grid. The methods are implemented in a MODFLOW-compatible solute-transport model and use MODFLOW's Multi-Node Well Package but are generic and can be readily implemented in other solute-transport models. For nonpumping multinode wells (used to simulate open boreholes or observation wells, for example) and for low-rate pumping wells (in which the flow between the well and the ground water system is not unidirectional), a simple routing and local mixing model was developed to calculate nodal concentrations within the borehole. For high-rate pumping multinode wells (either withdrawal or injection, in which flow between the well and the ground water system is in the same direction at all well nodes), complete and instantaneous mixing in the wellbore of all inflows is assumed.

Simulation of future stream alkalinity under changing deposition and climate scenarios.

Models of soil and stream water acidification have typically been applied under scenarios of changing acidic deposition, however, climate change is usually ignored. Soil air CO2 concentrations have potential to increase as climate warms and becomes wetter, thus affecting soil and stream water chemistry by initially increasing stream alkalinity at the expense of reducing base saturation levels on soil exchange sites. We simulate this change by applying a series of physically based coupled models capable of predicting soil air CO2 and stream water chemistry. We predict daily stream water alkalinity for a small catchment in the Virginia Blue Ridge for 60 years into the future given stochastically generated daily climate values. This is done for nine different combinations of climate and deposition. The scenarios for both climate and deposition include a static scenario, a scenario of gradual change, and a scenario of abrupt change. We find that stream water alkalinity continues to decline for all scenarios (average decrease of 14.4 microeq L-1) except where climate is gradually warming and becoming more moist (average increase of 13 microeq L-1). In all other scenarios, base cation removal from catchment soils is responsible for limited alkalinity increase resulting from climate change. This has implications given the extent that acidification models are used to establish policy and legislation concerning deposition and emissions.

Simulation of solute transport across low-permeability barrier walls.

Low-permeability, non-reactive barrier walls are often used to contain contaminants in an aquifer. Rates of solute transport through such barriers are typically many orders of magnitude slower than rates through the aquifer. Nevertheless, the success of remedial actions may be sensitive to these low rates of transport. Two numerical simulation methods for representing low-permeability barriers in a finite-difference groundwater-flow and transport model were tested. In the first method, the hydraulic properties of the barrier were represented directly on grid cells and in the second method, the intercell hydraulic-conductance values were adjusted to approximate the reduction in horizontal flow, allowing use of a coarser and computationally efficient grid. The alternative methods were tested and evaluated on the basis of hypothetical test problems and a field case involving tetrachloroethylene (PCE) contamination at a Superfund site in New Hampshire. For all cases, advective transport across the barrier was negligible, but preexisting numerical approaches to calculate dispersion yielded dispersive fluxes that were greater than expected. A transport model (MODFLOW-GWT) was modified to (1) allow different dispersive and diffusive properties to be assigned to the barrier than the adjacent aquifer and (2) more accurately calculate dispersion from concentration gradients and solute fluxes near barriers. The new approach yields reasonable and accurate concentrations for the test cases.

A finite-volume ELLAM for three-dimensional solute-transport modeling.

A three-dimensional finite-volume ELLAM method has been developed, tested, and successfully implemented as part of the U.S. Geological Survey (USGS) MODFLOW-2000 ground water modeling package. It is included as a solver option for the Ground Water Transport process. The FVELLAM uses space-time finite volumes oriented along the streamlines of the flow field to solve an integral form of the solute-transport equation, thus combining local and global mass conservation with the advantages of Eulerian-Lagrangian characteristic methods. The USGS FVELLAM code simulates solute transport in flowing ground water for a single dissolved solute constituent and represents the processes of advective transport, hydrodynamic dispersion, mixing from fluid sources, retardation, and decay. Implicit time discretization of the dispersive and source/sink terms is combined with a Lagrangian treatment of advection, in which forward tracking moves mass to the new time level, distributing mass among destination cells using approximate indicator functions. This allows the use of large transport time increments (large Courant numbers) with accurate results, even for advection-dominated systems (large Peclet numbers). Four test cases, including comparisons with analytical solutions and benchmarking against other numerical codes, are presented that indicate that the FVELLAM can usually yield excellent results, even if relatively few transport time steps are used, although the quality of the results is problem-dependent.

Anthropogenic sources of arsenic and copper to sediments in a suburban lake, Northern Virginia.

Mass balances of total arsenic and copper for a suburban lake in densely populated northern Virginia were calculated using date collected during 1998. Mass-balance terms were precipitation; stream inflow, including road runoff; stream outflow; and contributions from leaching of pressure-treated lumber. More mass of arsenic and copper was input to the lake than was output the 1998 lake-retention rates were 70% for arsenic and 20% for copper. The arsenic mass balance compared well with a calculated annual mass accumulation in the top 1 cm of the lake sediments; however, the calculated contribution of copper to the lake was insufficient to account for the amount of copper in this zone. Leaching experiments were conducted on lumber treated with chromated copper arsenate (CCA) to quantify approximate amounts of arsenic and copper contributed by this source. Sources to lake sediments included leaching of CCA-treated lumber (arsenic, 50%; copper, 4%), streamwater (arsenic, 50%; copper, 90%), and atmospheric deposition (arsenic, 1%; copper, 3%). Results of this study suggest that CCA-treated lumber and road runoff could be significant nonpoint sources of arsenic and copper, respectively, in suburban catchments.

Aggregating Fine-Scale Ecological Knowledge to Model Coarser-Scale Attributes of Ecosystems.

As regional and global scales become more important to ecologists, methods must be developed for the application of existing fine-scale knowledge to predict coarser-scale ecosystem properties. This generally involves some form of model in which fine-scale components are aggregated. This aggregation is necessary to avoid the cumulative error associated with the estimation of a large number of parameters. However, aggregation can itself produce errors that arise because of the variation among the aggregated components. The statistical expectation operator can be used as a rigorous method for translating fine-scale relationships to coarser scales without aggregation errors. Unfortunately this method is too cumbersome to be applied in most cases, and alternative methods must be used. These alternative methods are typically partial corrections for the variation in only a few of the fine-scale attributes. Therefore, for these methods to be effective, the attributes that are the most severe sources of error must be identified a priori. We present a procedure for making these identifications based on a Monte Carlo sampling of the fine-scale attributes. We then present four methods of translating fine-scale knowledge so it can be applied at coarser scales: (1) partial transformations using the expectation operator, (2) moment expansions, (3) partitioning, and (4) calibration. These methods should make it possible to apply the vast store of fine-scale ecological knowledge to model coarser-scale attributes of ecosystems.