Regional estimation of catchment‐scale soil properties by means of streamflow recession analysis for use in distributed hydrological models

The estimation of catchment‐scale soil properties, such as water storage capacity and hydraulic conductivity, is of primary interest for the implementation of distributed hydrological models at the regional scale. This estimation is generally performed on the basis of information provided by soil databases. However, such databases are often established for agronomic uses and generally do not document deep‐weathered rock horizons (i.e. pedologic horizons of type C and deeper), which can play a major role in water transfer and storages. Here, we define the Drainable Storage Capacity Index (DSCI), an indicator that relies on the comparison between cumulated streamflow and precipitation to assess catchment‐scale storage capacities. DSCI is found to be reliable to detect underestimation of soil storage capacities in soil databases. We also use the streamflow recession analysis methodology defined by Brutsaert and Nieber in 1977 to estimate water storage capacities and lateral saturated hydraulic conductivities of the nondocumented deep horizons. The analysis is applied to a sample of 23 catchments (0.2–291 km²) located in the Cévennes‐Vivarais region (south of France). For regionalization purposes, the obtained results are compared with the dominant catchment geology and present a clear hierarchy between the different geologies of the area. Hard crystalline rocks are found to be associated with the thickest and less conductive deep soil horizons. Schist rocks present intermediate values of thickness and of saturated hydraulic conductivity, whereas sedimentary rocks and alluvium are found to be less thick and most conductive. These results are of primary interest in view of the future set‐up of distributed hydrological models over the Cévennes‐Vivarais region. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling water fluxes for the analysis of diffuse pollution at the river basin scale

Diffuse pollution is a significant and sometimes even major component of surface water pollution. Diffuse inputs of pollutants to the surface water are related to runoff of precipitation. This means that the analysis of diffuse pollutant fluxes from the land surface to the surface water requires an analysis of water fluxes. In this paper we have modelled the average long‐term total runoff, groundwater recharge index and groundwater residence times for two large European river basins (Rhine and Elbe). We applied and compared two independently developed and recently published methods. We found that with the available large‐scale databases and methods we could simulate successfully the regional patterns of the average long‐term total runoff. The reported groundwater recharge indices and groundwater residence times should be interpreted as estimates based on available knowledge and databases. They do not represent absolute values, but illustrate the possible travel times and spatial patterns of the different runoff components that have to be taken into account for the analysis of diffuse pollution at large regional and temporal scales. Copyright © 2000 John Wiley & Sons, Ltd.     

Data‐Driven Approach for Analyzing Hydrogeology and Groundwater Quality Across Multiple Scales

Recent trends of assimilating water well records into statewide databases provide a new opportunity for evaluating spatial dynamics of groundwater quality and quantity. However, these datasets are scarcely rigorously analyzed to address larger scientific problems because they are of lower quality and massive. We develop an approach for utilizing well databases to analyze physical and geochemical aspects of groundwater systems, and apply it to a multiscale investigation of the sources and dynamics of chloride (Cl^?) in the near‐surface groundwater of the Lower Peninsula of Michigan. Nearly 500,000 static water levels (SWLs) were critically evaluated, extracted, and analyzed to delineate long‐term, average groundwater flow patterns using a nonstationary kriging technique at the basin‐scale (i.e., across the entire peninsula). Two regions identified as major basin‐scale discharge zones—the Michigan and Saginaw Lowlands—were further analyzed with regional‐ and local‐scale SWL models. Groundwater valleys (“discharge” zones) and mounds (“recharge” zones) were identified for all models, and the proportions of wells with elevated Cl^? concentrations in each zone were calculated, visualized, and compared. Concentrations in discharge zones, where groundwater is expected to flow primarily upwards, are consistently and significantly higher than those in recharge zones. A synoptic sampling campaign in the Michigan Lowlands revealed concentrations generally increase with depth, a trend noted in previous studies of the Saginaw Lowlands. These strong, consistent SWL and Cl^? distribution patterns across multiple scales suggest that a deep source (i.e., Michigan brines) is the primary cause for the elevated chloride concentrations observed in discharge areas across the peninsula.     

Bayesian Chance‐Constrained Hydraulic Barrier Design under Geological Structure Uncertainty

The groundwater community has widely recognized geological structure uncertainty as a major source of model structure uncertainty. Previous studies in aquifer remediation design, however, rarely discuss the impact of geological structure uncertainty. This study combines chance‐constrained (CC) programming with Bayesian model averaging (BMA) as a BMA‐CC framework to assess the impact of geological structure uncertainty in remediation design. To pursue this goal, the BMA‐CC method is compared with traditional CC programming that only considers model parameter uncertainty. The BMA‐CC method is employed to design a hydraulic barrier to protect public supply wells of the Government St. pump station from salt water intrusion in the “1500‐foot” sand and the “1700‐foot” sand of the Baton Rouge area, southeastern Louisiana. To address geological structure uncertainty, three groundwater models based on three different hydrostratigraphic architectures are developed. The results show that using traditional CC programming overestimates design reliability. The results also show that at least five additional connector wells are needed to achieve more than 90% design reliability level. The total amount of injected water from the connector wells is higher than the total pumpage of the protected public supply wells. While reducing the injection rate can be achieved by reducing the reliability level, the study finds that the hydraulic barrier design to protect the Government St. pump station may not be economically attractive.     

Quantifying large‐scale historical formation of accommodation in the Mississippi Delta

Large volumes of new accommodation have formed within the Mississippi Delta plain since the mid‐1950s in association with rapid conversion of coastal wetlands to open water. The three‐dimensional aspects and processes responsible for accommodation formation were quantified by comparing surface elevations, water depths, and vertical displacements of stratigraphic contacts that were correlated between short sediment cores. Integration of data from remotely sensed images, sediment cores, and water‐depth surveys at 10 geologically diverse areas in the delta plain provided a basis for estimating the total volume of accommodation formed by interior‐wetland subsidence and subsequent erosion. Results indicate that at most of the study areas subsidence was a greater contributor than erosion to the formation of accommodation associated with wetland loss. Tens of millions of cubic meters of accommodation formed rapidly at each of the large open‐water bodies that were formerly continuous interior delta‐plain marsh. Together the individual study areas account for more than 440 × 10⁶ × m³ of new accommodation that formed as holes in the Mississippi River delta‐plain fabric between 1956 and 2004. This large volume provides an estimate of the new sediment that would be needed just at the study areas to restore the delta‐plain wetlands to their pre‐1956 areal extent and elevations. Published 2010. This article is a US Government work and is in the public domain in the USA.     

VS2DRTI: Simulating Heat and Reactive Solute Transport in Variably Saturated Porous Media

Variably saturated groundwater flow, heat transport, and solute transport are important processes in environmental phenomena, such as the natural evolution of water chemistry of aquifers and streams, the storage of radioactive waste in a geologic repository, the contamination of water resources from acid‐rock drainage, and the geologic sequestration of carbon dioxide. Up to now, our ability to simulate these processes simultaneously with fully coupled reactive transport models has been limited to complex and often difficult‐to‐use models. To address the need for a simple and easy‐to‐use model, the VS2DRTI software package has been developed for simulating water flow, heat transport, and reactive solute transport through variably saturated porous media. The underlying numerical model, VS2DRT, was created by coupling the flow and transport capabilities of the VS2DT and VS2DH models with the equilibrium and kinetic reaction capabilities of PhreeqcRM. Flow capabilities include two‐dimensional, constant‐density, variably saturated flow; transport capabilities include both heat and multicomponent solute transport; and the reaction capabilities are a complete implementation of geochemical reactions of PHREEQC. The graphical user interface includes a preprocessor for building simulations and a postprocessor for visual display of simulation results. To demonstrate the simulation of multiple processes, the model is applied to a hypothetical example of injection of heated waste water to an aquifer with temperature‐dependent cation exchange. VS2DRTI is freely available public domain software.     

Community‐Based Groundwater Monitoring Network Using a Citizen‐Science Approach

Water level monitoring provides essential information about the condition of aquifers and their responses to water extraction, land‐use change, and climatic variability. It is important to have a spatially distributed, long‐term monitoring well network for sustainable groundwater resource management. Community‐based monitoring involving citizen scientists provides an approach to complement existing government‐run monitoring programs. This article demonstrates the feasibility of establishing a large‐scale water level monitoring network of private water supply wells using an example from Rocky View County (3900 km²) in Alberta, Canada. In this network, community volunteers measure the water level in their wells, and enter these data through a web‐based data portal, which allows the public to view and download these data. The close collaboration among the university researchers, county staff members, and community volunteers enabled the successful implementation and operation of the network for a 5‐year pilot period, which generated valuable data sets. The monitoring program was accompanied by education and outreach programs, in which the educational materials on groundwater were developed in collaboration with science teachers from local schools. The methodology used in this study can be easily adopted by other municipalities and watershed stewardship groups interested in groundwater monitoring. As governments are starting to rely increasingly on local municipalities and conservation authorities for watershed management and planning, community‐based groundwater monitoring provides an effective and affordable tool for sustainable water resources management.     

Reverse Water‐Level Fluctuations Associated with Fracture Connectivity

Reverse water‐level fluctuations (RWFs), a phenomenon in which water levels rise briefly in response to pumping, were detected in monitoring wells in a fractured siliciclastic aquifer system near a deep public supply well. The magnitude and timing of RWFs provide important information that can help interpret aquifer hydraulics near pumping wells. A RWF in a well is normally attributed to poroelastic coupling between the solid and fluid components in an aquifer system. In addition to revealing classical pumping‐induced poroelastic RWFs, data from pressure transducers located at varying depths and distances from the public supply well suggest that the RWFs propagate rapidly through fractures to influence wells hundreds of meters from the pumping well. The rate and cycling frequency of pumping is an important factor in the magnitude of RWFs. The pattern of RWF propagation can be used to better define fracture connectivity in an aquifer system. Rapid, cyclic head changes due to RWFs may also serve as a mechanism for contaminant transport.     

The relationship between dissolved organic carbon in stream water and soil organic carbon pools at different spatial scales

The relationship between stream water DOC concentrations and soil organic C pools was investigated at a range of spatial scales in subcatchments of the River Dee system in north‐east Scotland. Catchment percentage peat cover and soil C pools, calculated using local, national and international soils databases, were related to mean DOC concentrations in streams draining small‐ (<5 km²), medium‐ (12–38 km²) and large‐scale (56–150 km²) catchments. The results show that, whilst soil C pool is a good predictor of stream water DOC concentration at all three scales, the strongest relationships were found in the small‐scale catchments. In addition, in both the small‐ and large‐scale catchments, percentage peat cover was as a good predictor of stream water DOC concentration as catchment soil C pool. The data also showed that, for a given soil C pool, streams draining lowland (<700 m) catchments had higher DOC concentrations than those draining upland (>700 m) catchments, suggesting that disturbance and land use may have a small effect on DOC concentration. Our results therefore suggest that the relationship between stream water DOC concentration and catchment soil C pools exists at a range of spatial scales and this relationship appears to be sufficiently robust to be used to predict the effects of changes in catchment soil C storage on stream water DOC concentration. Copyright © 1999 John Wiley & Sons, Ltd.     

Surfactant effects on the water‐stable aggregation of wettable soils from the continental USA

Surfactants may affect soil structure differently depending upon the soil or the quality of rainfall or irrigation water. This study examined whether the water‐stable aggregation of 11 wettable soils was affected by surfactants and the water in which the soils were sieved. The study also examined whether the wettable soils' water drop penetration time (WDPT) was affected by surfactants, water drop quality, and elapsed time since the surfactants were applied. Two nonionic surfactants and a surfactant‐free water control were sprayed (by misting) upon air‐dry soil, then WDPT was measured 1 and 72 h thereafter. Subsequently, this treated soil was slowly wetted with an aerosol to its water content at a matric potential of ?3 kPa, then immediately sieved for 600 s in water that contained either appreciable or few electrolytes. Water‐stable aggregation, quantified as mean weight diameter (MWD), varied widely among soils, ranging from 0.10 to 1.36 mm. The MWDs were affected (at p = 0.06) by surfactant treatments, depending upon the soil but not sieving water quality. Surfactants affected the MWD of an Adkins loamy sand and Feltham sand, two of the three coarsest‐textured soils. Although WDPTs never exceeded 5 s, depending upon the soil WDPTs were affected by surfactant treatments but not by water drop quality. After surfactant application, WDPTs generally decreased with time for three soils but increased with time for one soil. Findings suggested that surfactants interacted (1) with clay mineralogy to affect MWD and (2) with soluble calcium to affect WDPT for certain soils. Surfactant treatments but not water quality affected both MWD and WDPT for some but not all of 11 wettable, US soils. Published 2012. This article is a US Government work and is in the public domain in the USA.     

Sub-Saharan African Ground Water Protection—Building on International Experience

Sub-Saharan Africa faces significant challenges in dealing with ground water pollution. These countries can look to successes and missteps on other continents to help choose their own individual paths to ensuring reliable and clean supplies of ground water. In the large view, sub-Saharan Africa can define specific levels of acceptable risk in water quality that drive cleanup efforts and are amenable to acceptance across national and geographic boundaries. Ground water quality databases must be expanded, and data must be available in an electronic form that is flexible, expandable, and uniform, and that can be used over wide geographic areas. Guidance from other continents is available on well construction, sampling and monitoring, interim remediation, technical impracticability, monitored natural attenuation, and many specific issues such as how to deal with small waste generators and septic contamination of water supply wells. It is important to establish a common African view on the appropriateness of other nations’ ground water quality guidance for African issues, economic conditions, and community circumstances. Establishing numerical, concentration-based, water quality action levels for pollutants in ground water, which many neighboring African nations could hold comparable, would set the stage for risk-based remediation of contaminated sites. Efforts to gain public, grass-roots understanding and support for stable and balanced enforcement of standards are also key. Finally, effective capacity building in the region could be an eventual solution to ground water quality problems; with increased numbers of trained environmental professionals, ground water throughout the region can be protected and contaminated sites cleaned up.     

Availability and Content of Domestic Well Records in the United States

Ground water professionals within public and private sectors use well records as data sources. Both the availability and the technical content of domestic well records in the 50 states are of interest to them. Well record availability is dependent on legal requirements, filing systems, databases, and storage places. Forty-six states have statewide regulations or other legal requirements for filing completion reports for domestic wells. Fifty-one agencies across the country maintain domestic well records. Filing systems involve location, chronology, date, or number information. Thirty-one of the 51 agencies maintain varying types of databases containing record information or data related to the records. Overall, records are kept in central offices in 43 states and in regional offices in three states.
The technical content of the records was evaluated for general, location, hydrogeology, and well construction information to assess the relative value of the records for use in national pesticide surveys. Technical information tabulated from the well records collected for this paper included nine items in a general category, 28 items in a well-construction category, eight items in a hydrogeology category, and six items in a location category. Items in the general and location categories identified the well location and ownership. Construction category items include those describing well-construction parameters such as grout, casing, and screen. Hydrogeology category items include static water level, aquifer media, and estimated yield. The three items always requested were owner's name, driller's name, and static water level. The three least-requested items, ranging from 16 percent to 10 percent, were packers, drilling fluid, and geologic formation.     

Water and Wastewater Treatment in Africa – Current Practices and Challenges

Sustainable access to safe drinking water and basic sanitation is an important part of the millennium development goals (MDGs). For most African countries, an extensive effort is needed for the last three remaining years for the achievement of the MDGs, especially in Sub‐Saharan Africa. Current practices for water and wastewater treatment in Africa are insufficient to ensure safe water and basic sanitation. To address this challenge, joint efforts are needed, including transforming to green economy, innovating technologies, improving operation and maintenance, harvesting energy, improving governance and management, promoting public participation, and establishing water quality standards.     

Information support for studying water resources of Lake Ladoga and its basin

A regional geographic information system for water resources of Lake Ladoga basin was developed as an information support for works aimed to study and practically use the water resources of Lake Ladoga and its basin. The system incorporates databases with the results of hydrochemical, hydrological, and meteorological measurements; information on water use; and an atlas of digital maps. The developed inquiry and communication system is based on web-interfaces, provides for remote access to databases, and enables interactive data exchange.     

Lateral and subsurface flows impact arctic coastal plain lake water budgets

Arctic thaw lakes are an important source of water for aquatic ecosystems, wildlife, and humans. Many recent studies have observed changes in Arctic surface waters related to climate warming and permafrost thaw; however, explaining the trends and predicting future responses to warming is difficult without a stronger fundamental understanding of Arctic lake water budgets. By measuring and simulating surface and subsurface hydrologic fluxes, this work quantified the water budgets of three lakes with varying levels of seasonal drainage, and tested the hypothesis that lateral and subsurface flows are a major component of the post‐snowmelt water budgets. A water budget focused only on post‐snowmelt surface water fluxes (stream discharge, precipitation, and evaporation) could not close the budget for two of three lakes, even when uncertainty in input parameters was rigorously considered using a Monte Carlo approach. The water budgets indicated large, positive residuals, consistent with up to 70% of mid‐summer inflows entering lakes from lateral fluxes. Lateral inflows and outflows were simulated based on three processes; supra‐permafrost subsurface inflows from basin‐edge polygonal ground, and exchange between seasonally drained lakes and their drained margins through runoff and evapotranspiration. Measurements and simulations indicate that rapid subsurface flow through highly conductive flowpaths in the polygonal ground can explain the majority of the inflow. Drained lakes were hydrologically connected to marshy areas on the lake margins, receiving water from runoff following precipitation and losing up to 38% of lake efflux to drained margin evapotranspiration. Lateral fluxes can be a major part of Arctic thaw lake water budgets and a major control on summertime lake water levels. Incorporating these dynamics into models will improve our ability to predict lake volume changes, solute fluxes, and habitat availability in the changing Arctic. Published 2016. This article is a U.S. Government work and is in the public domain in the USA.     

New Attempt to Reduce the Harm of Smoking: Reducing the Nitrosamines Level in Tobacco Smoke by Microwave Irradiation

In order to protect the environment and public health, microwave irradiation was investigated as a means of reducing the nitrosamines level of tobacco smoke under mild conditions where the microwave energy was 1 kW and the irradiation time was shorter than 2 min. The microwave‐induced elimination of nitrosamines such as N‐nitrosopyrrolidine (NPYR) and N‐nitrosodiphenylamine (NDPA) was investigated and the impact of the presence of water and salt in the medium on the reduction was assessed in detail. The existence of water in the medium was crucial for the microwave‐induced reduction of NDPA, because the rapid movement of water molecules under the influence of microwaves strongly affects the adjacent nitrosamines leading to their decomposition. This was the first time that microwaves have been used to directly irradiate cigarettes sealed in their packaging, and the nitrosamines content of tobacco and smoke of these irradiated cigarettes were analyzed with care. Irradiation for 90 s with microwave energy of 1kW was seen to reduce 25 to 30% of the nitrosamines in tobacco and 50 to 60% of the levels present in the mainstream smoke of the cigarette. Moreover, lower mutagenic activity and higher CHO cell livability were found for the irradiated cigarettes in comparison with the control according to the results of in vitro tests.     

Estimation of Generic Subslab Attenuation Factors for Vapor Intrusion Investigations

Generic indoor air:subslab soil gas attenuation factors (SSAFs) are important for rapid screening of potential vapor intrusion risks in buildings that overlie soil and groundwater contaminated with volatile chemicals. Insufficiently conservative SSAFs can allow high‐risk sites to be prematurely excluded from further investigation. Excessively conservative SSAFs can lead to costly, time‐consuming, and often inconclusive actions at an inordinate number of low‐risk sites. This paper reviews two of the most commonly used approaches to develop SSAFs: (1) comparison of paired, indoor air and subslab soil gas data in empirical databases and (2) comparison of estimated subslab vapor entry rates and indoor air exchange rates (IAERs). Potential error associated with databases includes interference from indoor and outdoor sources, reliance on data from basements, and seasonal variability. Heterogeneity in subsurface vapor plumes combined with uncertainty regarding vapor entry points calls into question the representativeness of limited subslab data and diminishes the technical defensibility of SSAFs extracted from databases. The use of reasonably conservative vapor entry rates and IAERs offers a more technically defensible approach for the development of generic SSAF values for screening. Consideration of seasonal variability in building leakage rates, air exchange rates, and interpolated vapor entry rates allows for the development of generic SSAFs at both local and regional scales. Limitations include applicability of the default IAERs and vapor entry rates to site‐specific vapor intrusion investigations and uncertainty regarding applicability of generic SSAFs to assess potential short‐term (e.g., intraday) variability of impacts to indoor air.     

The Daily Erosion Project – daily estimates of water runoff,soil detachment,and erosion

Water runoff and sediment transport from agricultural uplands are substantial threats to water quality and sustained crop production. To improve soil and water resources, farmers, conservationists, and policy‐makers must understand how landforms, soil types, farming practices, and rainfall interact with water runoff and soil erosion processes. To that end, the Iowa Daily Erosion Project (IDEP) was designed and implemented in 2003 to inventory these factors across Iowa in the United States. IDEP utilized the Water Erosion Prediction Project (WEPP) soil erosion model along with radar‐derived precipitation data and government‐provided slope, soil, and management information to produce daily estimates of soil erosion and runoff at the township scale (93 km² [36 mi²]). Improved national databases and evolving remote sensing technology now permit the derivation of slope, soil, and field‐level management inputs for WEPP. These remotely sensed parameters, along with more detailed meteorological data, now drive daily WEPP hillslope soil erosion and water runoff estimates at the small watershed scale, approximately 90 km² (35 mi²), across sections of multiple Midwest states. The revisions constitute a substantial improvement as more realistic field conditions are reflected, more detailed weather data are utilized, hill slope sampling density is an order of magnitude greater, and results are aggregated based on surface hydrology enabling further watershed research and analysis. Considering these improvements and the expansion of the project beyond Iowa it was renamed the Daily Erosion Project (DEP). Statistical and comparative evaluations of soil erosion simulations indicate that the sampling density is adequate and the results are defendable. The modeling framework developed is readily adaptable to other regions given suitable inputs. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Flexural deformation capacity of rectangular RC columns determined by the CAE method

A non‐parametric empirical approach, called the conditional average estimator (CAE) method, has been implemented for the estimation of the flexural deformation capacity of reinforced concrete rectangular columns expressed in terms of the ultimate (‘near collapse’) drift. Two databases (PEER and Fardis), which represent subsets of the original databases, were used. Four input parameters were employed in the basic model: axial load index, index related to confinement, shear span index, and concrete compressive strength. The results of analyses suggest that, in general, ultimate drift decreases with increasing axial load index, and increases with better confinement. An increase in the shear span‐to‐depth ratio has a beneficial effect until a turning point is reached. After that the opposite trend can be observed, i.e. a decrease in the ultimate drift with further increasing of the shear span‐to‐depth ratio. No clear trend is observed in the case of concrete compressive strength. The predictions, obtained by using the Fardis database are in general somewhat larger than the predictions from the PEER database, due to the difference in the definition of ultimate drift. The scatter of results is large. The local coefficient of variation, which is a measure for dispersion, amounts to about 0.2–0.5. The ultimate drifts obtained by using the two databases, were compared with the values predicted by the Eurocode 8 empirical formula. Copyright © 2006 John Wiley & Sons, Ltd.     

A linear systems approach to watershed transport simulation

Models that simulate loadings of pollutants from agricultural landscapes to surface waters often operate at time scales that are relatively coarse (e.g. daily) compared with how fast water moves in streams, suggesting a commensurate physical scale that is substantially larger than typical agricultural fields. In general, as pollutants enter water and move downstream, longitudinal dispersive effects and travel time de‐synchronization tend to cause flattening and broadening of concentration peaks—an effect with implications for potential impacts on ecological and human health, and for which adequate representation is thus important for risk assessment. In‐stream transport is often approximated in practice using numerical implementation of the one‐dimensional advection–dispersion equation (ADE), with streams discretized into linked homogeneous segments. However, when a daily time step is employed, limitations inherent in the finite difference methodology may constrain simulated dispersion in lotic waters to unrepresentative or unrealistic magnitudes. In this paper, a convolution‐based approach to surface water transport is suggested as an alternative to the ADE, for use in combination with daily input loading models. This approach offers the advantage of greater flexibility in representing longitudinal mixing by using impulse response functions (IRF) to represent inter‐segment transport. Networks of stream segments are represented using nested convolutions, implemented using forward and inverse discrete Fourier transform to simplify calculations. Enhanced representational flexibility arises from the freedom afforded the modeller in selecting each segment's IRF, which may be chosen to represent dispersive regimes ranging from pure advection (plug flow) to compete mixing, and beyond to the sort of long‐tailed mixing characterized by fractal inverse frequency power‐law scaling. The approach is explored in proof‐of‐concept exercises that make use of atrazine monitoring data sets collected over common time periods from upstream and downstream locations within the same watersheds. Published 2014. This article is a U.S. Government work and is in the public domain in the USA.