Groundwaters at Risk: Wetland Loss Changes Sources,Lengthens Pathways,and Decelerates Rejuvenation of Groundwater Resources

Wetland loss alters the hydrology of wetlandscapes in poorly understood ways. To quantify the effects of wetland loss on subsurface hydrology, a physically based hydrologic model that simulates the timing and pathways of subsurface hydrologic connections was coupled with wetland inventories over a 50‐year period during which substantial wetland loss occurred. The model revealed, based on vertical variations in saturated hydraulic conductivities, wetland loss of different degrees led to a contraction of catchment contributing areas to local surface waters but an expansion of contributing areas to the regional surface water body. This shift in groundwater contributing areas reflected (1) a decrease in baseflow contribution to the local surface water bodies, and (2) an increase in the transit time and length of subsurface hydrologic connections with an associated increase in the magnitude and age of baseflow discharging to the regional surface water body. The model also showed regions with thick permeable aquifers were particularly sensitive to the loss of wetlands. Our ability to predict these changes in hydrology of the watershed provides important support for designing science‐based policies to promote sustainable water resource management.     

土地利用对流域氮素记忆效应的影响研究

"记忆效应"是一些流域实行了多年的氮污染控制实践却仍未见成效的主要原因之一.然而,目前对流域氮素记忆效应的认识还很有限.本研究利用温州市珊溪水源地13个集水区近5年的河流TN浓度逐月监测数据,运用重标极差(R/S)分析和Spearman秩相关分析方法,探讨了土地利用对流域氮素记忆效应的影响.结果表明,河流TN平均浓度与耕地(r=0.495,p=0.086)和建设用地(r=0.692,p=0.009)面积比例正相关,与林地(r=-0.604,p=0.029)面积比例显著负相关,与园地(r=-0.039,p=0.900)、未利用地(r=-0.176,p=0.566)、集水区面积(r=-0.335,p=0.263)和河网密度(r=0.148,p=0.629)的相关性不显著.13条入库支流TN浓度序列的Hurst指数变异范围为0.33~0.72,意味着其中11个集水区具有显著的氮素记忆效应(0.5Hurst指数1.0),而其余2个集水区不显著(0Hurst指数0.5).Hurst指数与耕地(r=-0.482,p=0.095)和建设用地(r=-0.311,p=0.301)面积比例均负相关,而与园地(r=0.479,p=0.098)和林地(r=0.510,p=0.075)面积比例正相关.同时,Hurst指数与集水区面积大小(r=0.118,p=0.700)和未利用地面积比例(r=0.032,p=0.917)的相关性不显著,而与河网密度负相关(r=-0.529,p=0.063).流域内"源"功能为主土地利用面积的增加和氮输移效率的提高加剧了河流氮污染,而对氮素记忆效应的作用则相反.因此,尽管增加"汇"景观面积、降低氮输移效率能缓解河流氮污染,但"汇"景观中截留的遗留氮对河流氮浓度会造成更长时间的影响.     

Managing Uncertainty in Runoff Estimation with the U.S. Environmental Protection Agency National Stormwater Calculator

The U.S. Environmental Protection Agency National Stormwater Calculator (NSWC) simplifies the task of estimating runoff through a straightforward simulation process based on the EPA Stormwater Management Model. The NSWC accesses localized climate and soil hydrology data, and options to experiment with low‐impact development (LID) features for parcels up to 5 ha in size. We discuss how the NSWC treats the urban hydrologic cycle and focus on the estimation uncertainty in soil hydrology and its impact on runoff simulation by comparing field‐measured soil hydrologic data from 12 cities to corresponding NSWC estimates in three case studies. The default NSWC hydraulic conductivity is 10.1 mm/h, which underestimates conductivity measurements for New Orleans, Louisiana (95 ± 27 mm/h) and overestimates that for Omaha, Nebraska (3.0 ± 1.0 mm/h). Across all cities, the NSWC prediction, on average, underestimated hydraulic conductivity by 10.5 mm/h compared to corresponding measured values. In evaluating how LID interact with soil hydrology and runoff response, we found direct hydrologic interaction with pre‐existing soil shows high sensitivity in runoff prediction, whereas LID isolated from soils show less impact. Simulations with LID on higher permeability soils indicate that nearly all of pre‐LID runoff is treated; while features interacting with less‐permeable soils treat only 50%. We highlight the NSWC as a screening‐level tool for site runoff dynamics and its suitability in stormwater management.     

An Investigation of Errors in Distributed Models' Stream Discharge Prediction Due to Channel Routing

Over the summer of 2015, the National Water Center hosted the National Flood Interoperability Experiment (NFIE) Summer Institute. The NFIE organizers introduced a national‐scale distributed hydrologic modeling framework that can provide flow estimates at around 2.67 million reaches within the continental United States. The framework generates discharges by coupling a given Land Surface Model (LSM) with the Routing Application for Parallel Computation of Discharge (RAPID). These discharges are then accumulated through the National Hydrography Dataset Plus stream network. The framework can utilize a variety of LSMs to provide the runoff maps to the routing component. The results obtained from this framework suggested that there still exists room for further enhancements to its performance, especially in the area of peak timing and magnitude. The goal of our study was to investigate a single source of the errors in the framework's discharge estimates, which is the routing component. The authors substitute RAPID which is based on the simplified linear Muskingum routing method by the nonlinear routing component the Iowa Flood Center have incorporated in their full hydrologic Hillslope‐Link Model. Our results show improvement in model performance across scales due to incorporating new routing methodology.     

Trait differences between naturalized and invasive plant species independent of residence time and phylogeny

The ability to predict which alien plants will transition from naturalized to invasive prior to their introduction to novel regions is a key goal for conservation and has the potential to increase the efficacy of weed risk assessment (WRA). However, multiple factors contribute to plant invasion success (e.g., functional traits, range characteristics, residence time, phylogeny), and they all must be taken into account simultaneously in order to identify meaningful correlates of invasion success. We compiled 146 pairs of phylogenetically paired (congeneric) naturalized and invasive plant species in Australia with similar minimum residence times (i.e., time since introduction in years). These pairs were used to test for differences in 5 functional traits (flowering duration, leaf size, maximum height, specific leaf area [SLA], seed mass) and 3 characteristics of species’ native ranges (biome occupancy, mean annual temperature, and rainfall breadth) between naturalized and invasive species. Invasive species, on average, had larger SLA, longer flowering periods, and were taller than their congeneric naturalized relatives. Invaders also exhibited greater tolerance for different environmental conditions in the native range, where they occupied more biomes and a wider breadth of rainfall and temperature conditions than naturalized congeners. However, neither seed mass nor leaf size differed between pairs of naturalized and invasive species. A key finding was the role of SLA in distinguishing between naturalized and invasive pairs. Species with high SLA values were typically associated with faster growth rates, more rapid turnover of leaf material, and shorter lifespans than those species with low SLA. This suite of characteristics may contribute to the ability of a species to transition from naturalized to invasive across a wide range of environmental contexts and disturbance regimes. Our findings will help in the refinement of WRA protocols, and we advocate the inclusion of quantitative traits, in particular SLA, into the WRA schemes.     

Modeling Long‐Term Trends of Chlorinated Ethene Contamination at a Public Supply Well

A mass‐balance solute‐transport modeling approach was used to investigate the effects of dense nonaqueous phase liquid (DNAPL) volume, composition, and generation of daughter products on simulated and measured long‐term trends of chlorinated ethene (CE) concentrations at a public supply well. The model was built by telescoping a calibrated regional three‐dimensional MODFLOW model to the capture zone of a public supply well that has a history of CE contamination. The local model was then used to simulate the interactions between naturally occurring organic carbon that acts as an electron donor, and dissolved oxygen (DO), CEs, ferric iron, and sulfate that act as electron acceptors using the Sequential Electron Acceptor Model in three dimensions (SEAM3D) code. The modeling results indicate that asymmetry between rapidly rising and more gradual falling concentration trends over time suggests a DNAPL rather than a dissolved source of CEs. Peak concentrations of CEs are proportional to the volume and composition of the DNAPL source. The persistence of contamination, which can vary from a few years to centuries, is proportional to DNAPL volume, but is unaffected by DNAPL composition. These results show that monitoring CE concentrations in raw water produced by impacted public supply wells over time can provide useful information concerning the nature of contaminant sources and the likely future persistence of contamination.     

Variable Source Area Hydrology Modeling with the Water Erosion Prediction Project Model

In nondegraded watersheds of humid climates, subsurface flow patterns determine where the soil saturates and where surface runoff is occurring. Most models necessarily use infiltration‐excess (i.e., Hortonian) runoff for predicting runoff and associated constituents because subsurface flow algorithms are not included in the model. In this article, we modify the Water Erosion Prediction Project (WEPP) model to simulate subsurface flow correctly and to predict the spatial and temporal location of saturation, the associated lateral flow and surface runoff, and the location where the water can re‐infiltrate. The modified model, called WEPP‐UI, correctly simulated the hillslope drainage data from the Coweeta Hydrologic Laboratory hillslope plot. We applied WEPP‐UI to convex, concave, and S‐shaped hillslope profiles, and found that multiple overland flow elements are needed to simulate distributed lateral flow and runoff well. Concave slopes had the greatest runoff, while convex slopes had the least. Our findings concur with observations in watersheds with saturation‐excess overland flow that most surface runoff is generated on lower concave slopes, whereas on convex slopes runoff infiltrates before reaching the stream. Since the WEPP model is capable of simulating both saturation‐excess and infiltration‐excess runoff, we expect that this model will be a powerful tool in the future for managing water quality.     

A Water Yield‐Oriented Practical Approach for Multifunctional Forest Management and its Application in Dryland Regions of China

Mountainous forest areas are vitally important for water supply in dryland regions which suffer from high erosion risk and severe water shortage. Massive afforestation, mainly for erosion control, may reduce the water yield and threaten local water supply security. Moreover, many over‐dense forests due to a strict logging ban policy have produced remarkably negative impacts for both forests (e.g., low timber quality, restricted natural regeneration, and high stand instability) and water yield. To satisfy the rapidly increasing demands on water supply and other services, a practical approach for managing forest stands in a multifunctional way, which particularly addresses water yielding, is urgently required. For this purpose, we integrated the existing knowledge and experience, designed an “ideal” stand structure to represent multifunctional forest (MFF) and determined its key parameters (a ground coverage of >0.7, a canopy density around 0.7, and an H/DBH ratio (tree height [m] to the diameter at breast height [cm]) of <0.7). Moreover, a decision process for MFF stand management was recommended as: (1) investigating the site quality; (2) identifying the site‐specific main forest functions; (3) quantifying the stand structure; (4) diagnosing the stand structure by comparing with the “ideal” one; and (5) arranging the functions/structure‐oriented management measures. In this way, the water‐yielding function can be improved and meanwhile other forest functions can be promoted.     

Drainage Impacts on Surficial Water Retention Capacity of a Prairie Pothole Watershed

Wetland restoration has been proposed as a tool to mitigate excess runoff and associated nonpoint source pollution in the Upper Midwestern United States. This study quantified the surficial water retention capacity of existing and drained wetlands for the Greater Blue Earth River Basin (GBERB), an intensively drained agricultural watershed. Using airborne light detection and ranging, the historic depressional storage was determined to be 152 mm. Individual depression analysis suggested that the restoration of most drained areas would have little impact on the storage capacity of the GBERB because the majority (53%) of retention capacity was in large depressions (>40 ha) which comprised only a small proportion (<1.0) of the observed depressions. Accounting for change in storage and the difference in annual evapotranspiration (ET) between wetlands and the croplands that replaced them, restoration of all depressions in the Minnesota portion of GBERB would provide a maximum of 131 mm additional capacity over and above the modern day capacity (193 mm; 56 mm depressional storage; 60 mm wetland ET; and 77 mm cropland ET). Considering that depressional depths in smaller areas are within the range of uncertainty of the lidar digital elevation models and larger depressions have the most storage, we conclude that efforts to increase the surficial water‐holding capacity of the GBERB would be best served in the restoration of large (>40 ha) depressions.     

Empirical Estimation of Stream Discharge Using Channel Geometry in Low‐Gradient,Sand‐Bed Streams of the Southeastern Plains

Manning's equation is used widely to predict stream discharge (Q) from hydraulic variables when logistics constrain empirical measurements of in‐bank flow events. Uncertainty in Manning's roughness (n_M) is the major source of error in natural channels, and sand‐bed streams pose difficulties because flow resistance is affected by flow‐dependent bed configuration. Our study was designed to develop and validate models for estimating Q from channel geometry easily derived from cross‐sectional surveys and available GIS data. A database was compiled consisting of 484 Q measurements from 75 sand‐bed streams in Alabama, Georgia, South Carolina, North Carolina (Southeastern Plains), and Florida (Southern Coastal Plain), with six New Zealand streams included to develop statistical models to predict Q from hydraulic variables. Model error characteristics were estimated with leave‐one‐site‐out jackknifing. Independent data of 317 Q measurements from 55 Southeastern Plains streams indicated the model (Q = A_cR_H^0.6906S^0.1216; where A_c is the channel area, R_H is the hydraulic radius, and S is the bed slope) best predicted Q, based on Akaike's information criterion and root mean square error. Models also were developed from smaller Q range subsets to explore if subsets increased predictive ability, but error fit statistics suggested that these were not reasonable alternatives to the above equation. Thus, we recommend the above equation for predicting in‐bank Q of unbraided, sandy streams of the Southeastern Plains.