Do geographically isolated wetlands influence landscape functions?

Geographically isolated wetlands (GIWs), those surrounded by uplands, exchange materials, energy, and organisms with other elements in hydrological and habitat networks, contributing to landscape functions, such as flow generation, nutrient and sediment retention, and biodiversity support. GIWs constitute most of the wetlands in many North American landscapes, provide a disproportionately large fraction of wetland edges where many functions are enhanced, and form complexes with other water bodies to create spatial and temporal heterogeneity in the timing, flow paths, and magnitude of network connectivity. These attributes signal a critical role for GIWs in sustaining a portfolio of landscape functions, but legal protections remain weak despite preferential loss from many landscapes. GIWs lack persistent surface water connections, but this condition does not imply the absence of hydrological, biogeochemical, and biological exchanges with nearby and downstream waters. Although hydrological and biogeochemical connectivity is often episodic or slow (e.g., via groundwater), hydrologic continuity and limited evaporative solute enrichment suggest both flow generation and solute and sediment retention. Similarly, whereas biological connectivity usually requires overland dispersal, numerous organisms, including many rare or threatened species, use both GIWs and downstream waters at different times or life stages, suggesting that GIWs are critical elements of landscape habitat mosaics. Indeed, weaker hydrologic connectivity with downstream waters and constrained biological connectivity with other landscape elements are precisely what enhances some GIW functions and enables others. Based on analysis of wetland geography and synthesis of wetland functions, we argue that sustaining landscape functions requires conserving the entire continuum of wetland connectivity, including GIWs.

Modeling Water and Nutrient Movement in Sandy Soils Using HYDRUS-2D.

Models help to describe and predict complex processes and scenarios that are difficult to understand or measure in environmental management systems. Thus, model simulations were performed (i) to calibrate HYDRUS-2D for water and solute movement as a possible decision support system for Candler and Immokalee fine sand using data from microsprinkler and drip irrigation methods, (ii) to validate the performance of HYDRUS-2D using field data of microsprinkler and drip irrigation methods, and (iii) to investigate Br, NO, and water movement using annual or seasonal weather data and variable fertigation scenarios. The model showed reasonably good agreement between measured and simulated values for soil water content ( = 0.87-1.00), Br ( = 0.63-0.96), NO-N ( = 0.66-0.98), P ( = 0.25-0.78), and K ( = 0.44-0.99) movement. The model could be successfully used for scheduling irrigation and predicting nutrient leaching for both microsprinkler and drip irrigation systems on Florida's sandy soils.

Back diffusion from thin low permeability zones.

Aquitards can serve as long-term contaminant sources to aquifers when contaminant mass diffuses from the aquitard following aquifer source mass depletion. This study describes analytical and experimental approaches to understand reactive and nonreactive solute transport in a thin aquitard bounded by an adjacent aquifer. A series of well-controlled laboratory experiments were conducted in a two-dimensional flow chamber to quantify solute diffusion from a high-permeability sand into and subsequently out of kaolinite clay layers of vertical thickness 15 mm, 20 mm, and 60 mm. One-dimensional analytical solutions were developed for diffusion in a finite aquitard with mass exchange with an adjacent aquifer using the method of images. The analytical solutions showed very good agreement with measured breakthrough curves and aquitard concentration distributions measured in situ by light reflection visualization. Solutes with low retardation accumulated more stored mass with greater penetration distance in the aquitard compared to high-retardation solutes. However, because the duration of aquitard mass release was much longer, high-retardation solutes have a greater long-term back diffusion risk. The error associated with applying a semi-infinite domain analytical solution to a finite diffusion domain increases as a function of the system relative diffusion length scale, suggesting that the solutions using image sources should be applied in cases with rapid solute diffusion and/or thin clay layers. The solutions presented here can be extended to multilayer aquifer/low-permeability systems to assess the significance of back diffusion from thin layers.

Why wastewater treatment fails to protect stream ecosystems in Europe.

There is significant debate about why less than half of European rivers and streams are in good ecological status, despite decades of intense regulatory efforts. Of the multiple stressors that are recognized as potential contributors to stream degradation, we focus on discharge from 26,500 European wastewater treatment plants (WWTPs). We tested the hypothesis that stream ecological status degradation across Europe is related to the local intensity of wastewater discharge, with an expected stream-order (ω) dependence based on the scaling laws that govern receiving stream networks. We found that ecological status in streams (ω≤3) declined consistently with increasing urban wastewater discharge fraction of stream flow (UDF) across river types and basins. In contrast, ecological status in larger rivers (ω≥4) was not related to UDF. From a continental-scale logistic regression model (accuracy 86%) we identified an ecologically critical threshold UDF = 6.5% ± 0.5. This is exceeded by more than one third of WWTPs in Europe, mostly discharging into smaller streams. Our results suggest that new receiving water-specific strategies for wastewater management are needed to achieve good ecological status in smaller streams.

Seasonal dynamics of terrestrially sourced nitrogen influenced Karenia brevis blooms off Florida's southern Gulf Coast.

Harmful algal blooms (HABs) threaten coastal ecological systems, public health, and local economies, but the complex physical, chemical, and biological processes that culminate in HABs vary by locale and are often poorly understood. Despite broad recognition that cultural eutrophication may exacerbate nearshore bloom events, the association is typically not linear and is often difficult to quantify. Off the Gulf Coast of Florida, Karenia brevis blooms initiate in the open waters of the Gulf of Mexico, and advection of cells supplies nearshore blooms. However, past work has struggled to describe the relationship between terrestrial nutrient runoff and bloom maintenance near the Gulf Coast. This study applied a novel nonlinear time series (NLTS) analytical framework to investigate whether nearshore bloom dynamics observed near Charlotte Harbor, FL were causally and systematically driven by terrestrially sourced inputs of nitrogen, phosphorus, and freshwater between 2012 and 2018. Singular spectrum analysis (SSA) isolated low-dimensional, deterministic signals in K. brevis log10-density dynamics and in the dynamics of nine of 10 candidate drivers. The predominantly seasonal K. brevis signal was strong, explaining 77.6% of the total variance in the observed time series. Causal tests with convergent cross-mapping provided evidence that nitrogen concentrations measured at the discharge point of the Caloosahatchee River systematically influenced K. brevis bloom dynamics. However, further causal testing failed to link these nitrogen dynamics to an upstream basin, possibly due to data limitations. The results support the hypothesis that anthropogenic nitrogen runoff facilitated the growth of K. brevis blooms near Charlotte Harbor and suggest that bloom events would be mitigated by nitrogen source and transport controls within the Caloosahatchee and/or Kissimmee River basins. More broadly, this work demonstrates that management-relevant causal inferences into the drivers of HABs may be drawn from available monitoring records.

Emergent dispersal networks in dynamic wetlandscapes.

The connectivity among distributed wetlands is critical for aquatic habitat integrity and to maintain metapopulation biodiversity. Here, we investigated the spatiotemporal fluctuations of wetlandscape connectivity driven by stochastic hydroclimatic forcing, conceptualizing wetlands as dynamic habitat nodes in dispersal networks. We hypothesized that spatiotemporal hydrologic variability influences the heterogeneity in wetland attributes (e.g., size and shape distributions) and wetland spatial organization (e.g., gap distances), in turn altering the variance of the dispersal network topology and the patterns of ecological connectivity. We tested our hypotheses by employing a DEM-based, depth-censoring approach to assess the eco-hydrological dynamics in a synthetically generated landscape and three representative wetlandscapes in the United States. Network topology was examined for two end-member connectivity measures: centroid-to-centroid (C2C), and perimeter-to-perimeter (P2P), representing the full range of within-patch habitat preferences. Exponentially tempered Pareto node-degree distributions well described the observed structural connectivity of both types of networks. High wetland clustering and attribute heterogeneity exacerbated the differences between C2C and P2P networks, with Pareto node-degree distributions emerging only for a limited range of P2P configuration. Wetlandscape network topology and dispersal strategies condition species survival and biodiversity.

The evolution of human population distance to water in the USA from 1790 to 2010.

Human societies evolved alongside rivers, but how has the relationship between human settlement locations and water resources evolved over time? We conducted a dynamic analysis in the conterminous US to assess the coevolution of humans and water resources from 1790 to 2010. Here we show that humans moved closer to major rivers in pre-industrial periods but have moved farther from major rivers after 1870, demonstrating the dynamics of human reliance on rivers for trade and transport. We show that humans were preferentially attracted to areas overlying major aquifers since industrialization due to the emergent accessibility of groundwater in the 20th century. Regional heterogeneity resulted in diverse trajectories of settlement proximity to major rivers, with the attractiveness of rivers increasing in arid regions and decreasing in humid areas. Our results reveal a historical coevolution of human-water systems, which could inform water management and contribute to societal adaptation to future climate change.

In Situ Measurement of Nitrate Flux and Attenuation Using a Soil Passive Flux Meter.

This work enhances our understanding of catchment-scale N budgets by demonstrating the modification and application of a simple method for direct in situ measurements of vadose zone nitrate leaching and attenuation. We developed a soil passive flux meter (SPFM) to measure solute leaching based on a modified design of ion-exchange resin columns, and we tested the design in numerical simulations, laboratory experiments, plot-scale field experiments, and a catchment-scale field deployment. Our design minimized flow divergence around the resin column to attain nearly 100% capture of surface applied tracers in plot- and catchment-scale deployments. We found that mixing resin with native soil and extending the column height 10 cm above the resin layer minimized divergence of soil water around the column, resulting in a field-measured convergence factor (χ) of 1.3 that was consistent with numerical simulations. For catchment-scale testing, SPFMs were used at nine sites in three dominant land uses (crop, pasture, and turf) with known N inputs in two deployments, one during the 4-mo wet season and an additional set during the 8-mo dry season, to obtain integral annual measures of soil nitrate fluxes. In situ measured nitrate leaching determined from the SPFMs was positively correlated with known N inputs ( = 0.55, < 0.05) and attenuation averaged 67% (± 24% SD) of inputs across all sites. Although N inputs explain a large portion of the variability, our results emphasize the importance of both inter- and intra-land use variability in landscape-scale N budgets.

Diffusion of solutes from depleting sources into and out of finite low-permeability zones.

Two important factors that affect groundwater contaminant persistence are the temporal pattern of contaminant source depletion and solute diffusion into and out of aquitards. This study provides a framework to evaluate the relative importance of these effects on contaminant persistence, with emphasis on the importance of thin aquitards. We developed one-dimensional (1D) analytical solutions for forward and back diffusion in a finite domain with a no flux boundary using the method of images and demonstrated their applicability to measured data from three well-controlled laboratory diffusion experiments with exponentially depleting sources. We used both in situ aquitard solute concentrations and aquifer breakthrough curves for sorbing and non-sorbing solutes. The finite-domain no flux boundary solutions showed better agreement with measured data than was available with semi-infinite approaches, with increasing discrepancy for dimensionless relative diffusion length scale beyond a critical threshold value (Zd > 0.7). We also used a mass balance to demonstrate that the temporal pattern of contaminant source depletion controls the duration of solute mass accumulation in the aquitard, as well as the total solute mass release back into the aquifer. Lower rates of source depletion result in a longer period of mass accumulation in the aquitard and later back diffusion initiation time. The amount of solute mass stored in the aquitard increases with longer loading duration, thereby contributing to overall longer contaminant persistence in aquifers. This study entails widespread implications for anthropogenic waste and contamination sites, which are all dependent on efficient and cost-effective contaminant management strategies.

Multi-decadal trajectories of phosphorus loading, export, and instream retention along a catchment gradient.

Phosphorus inputs to many rivers have been reduced in recent decades to mitigate the damaging effects of eutrophication. However, reductions in total phosphorus (TP) inputs rarely correspond with ecological improvements of the river ecosystem. We analyzed a unique weekly long-term data set ranging from 1966 to 2013, covering seven monitoring sites in the Ruhr River in Germany. We identified the relative importance of different TP sources, quantified long-term trajectories of degradation and recovery, including the dynamics of TP retention, and assessed the response of chlorophyll-a (Chl-a) to increasing and decreasing TP concentrations along the whole river gradient. We found that the decline of TP loads at the beginning of the 1980s was dominantly triggered by a reduction of point sources. The cumulative TP retention capacity increased along the river gradient, increasing from effectively zero in the upstream section, to 26% and 36% of TP input in the upper midstream and lower downstream section. This pattern is consistent with higher prevalence of impoundments and weirs downstream, indicating that TP retention is likely associated with sedimentation posing a potential risk due to remobilization of legacy phosphorus. Degradation and recovery pathways differ from upstream to downstream. Along the river continuum we found three distinct types of reversible trajectories: 1. instream storage only during the recovery phase (upstream only); 2. instream storage in both degradation and recovery phases, but with significantly different characteristics depending on TP input load (midstream only); 3. higher instream storage during the recovery phase (downstream only). While in-stream TP loads may recover rapidly, the ecological response to altered nutrient inputs can be associated with considerable time-lags and decouplings between Chl-a and TP concentrations. Therefore, river systems may not return to historically good ecological status solely from massive nutrient reduction, but may also require other management activities.

River network connectivity and fish diversity.

Frequent and severe disruptions of natural river flows associated with human activities significantly alter hydrological connectivity in large river networks, with deleterious effects on fish diversity. Understanding the relationship between fish diversity and river network connectivity is fundamental to ensuring species persistence, ecosystem integrity, and human well-being. Here, we provide a review of the mechanisms by which river network connectivity (RNC) affects fish diversity. We review the relationships between forms, systems and types of RNC and fish diversity, based on more than 100 previous studies. In summary, sustaining RNC promotes fish diversity in longitudinal and lateral axes, and species sorting, dispersal dynamics, and habitat availability are the main factors driving the distribution of fish diversity, followed by nutrition and trophic dynamics. Our work highlights the effects of RNC on fish diversity, and provides a mechanistic understanding of how RNC affects fish diversity across river basins, thus providing scientific guidance for protecting fish biodiversity and improving the health of river network ecosystems.

Stochastic dynamics of wetlandscapes: Ecohydrological implications of shifts in hydro-climatic forcing and landscape configuration.

Wetlands are embedded in landscapes in fractal spatial patterns, and are characterized by highly dynamic, interlinked hydrological, biogeochemical, and ecological functions. We propose here a stochastic approach to evaluate and predict the spatiotemporal hydrologic variability of wetlands at landscape scale (100 km2). Stochastic hydro-climatic forcing (daily rainfall and evapotranspiration) and the landscape topographic setting (spatial structure of wetlands within the landscape) are key drivers of wetland eco-hydrologic functionality. The novelty of our approach lies in the quantification of the hydrological dynamics for all wetlands distributed in a given landscape, and in linking stochasticity of hydroclimatic forcing and ecologically meaningful wetland network metrics. We applied the modeling framework to investigate daily hydrologic dynamics in six landscapes across the U.S. that span gradients of hydroclimate and abundance of wetlands. We assess landscape-scale patterns using four key wetland hydrological attributes that have significance in terms of aquatic habitat suitability and dispersal: (1) Abundance (2) Diversity (3) Persistence, and (4) Accessibility. We observe that the hydrologic responses of each of the six landscapes are driven by the interactions between regional stochastic hydro-climatic forcing and landscape topographic setting. Despite differences in these features, similar scaling relations define diversity (area distributions) and accessibility (separation-distance distributions). Persistence of hydrologic regimes, defined by duration of inundation above thresholds, was least in more-arid settings, and higher in humid settings, consistent with intuitive understanding. These results can support assessments of the spatiotemporal variability of ecohydrological attributes in diverse wetlandscapes, including aquatic species dispersal and habitat suitability for unique flora and fauna.

Spatial patterns of water quality impairments from point source nutrient loads in Germany's largest national River Basin (Weser River).

We employed the well-established Horton-Strahler, hierarchical, stream-order (ω) scheme to investigate scaling of nutrient loads (P and N) from ~845 wastewater treatment plants (WWTPs) distributed along the river network in urbanized Weser River, the largest national basin in Germany (~46K km2; ~8.4 million population). We estimated hydrologic and water quality impacts at the reach- and basin-scales, at two steady river discharge conditions (median flow, QR50; low-flow, QR90). Of the five WWTPs class-sizes (1 ≤ k ≤ 5), ~68% discharge to small low-order streams (ω < 3). We found large variations in capacity to dilute WWTP nutrient loads because of variability in (1) treated wastewater discharge (QU) within and among different class-sizes, and (2) river discharge (QR) within low-order streams (ω < 3) resulting from differences in drainage areas. For QR50, reach-scale water quality impairment assessed by nutrient concentration was likely at 136 (~16%) locations for P and 15 locations (~2%) for N. About 90% of these locations were lower-order streams (ω < 3). At QR50 and only with dilution, basin-scale cumulative nutrient loads from multiple upstream WWTPs increase impaired locations to 266 (~32% of total) for P. Considering in-stream uptake decreased P-impaired streams to 225 (~27%), suggesting the dominant role of dilution in the Weser River basin. Role of in-stream uptake diminished along the flow paths, while dilution in larger streams (4 ≤ ω ≤ 7) minimizes the impact of WWTP loads. Under QR90 conditions [(QR50/QR90) ~ 2.5], water quality impaired locations will likely double for the basin-scale analyses. Long-term water quality data suggested that diffuse sources are the primary contributors for water quality impairments in large streams. Our data-modeling synthesis approach is transferable to other urbanized river basins and extends understanding of point source impacts on water quality across spatial scales.

Locally-calibrated light transmission visualization methods to quantify nonaqueous phase liquid mass in porous media.

Five locally-calibrated light transmission visualization (LTV) methods were tested to quantify nonaqueous phase liquid (NAPL) mass and mass reduction in porous media. Tetrachloroethylene (PCE) was released into a two-dimensional laboratory flow chamber packed with water-saturated sand which was then flushed with a surfactant solution (2% Tween 80) until all of the PCE had been dissolved. In all the LTV methods employed here, the water phase was dyed, rather than the more common approach of dyeing the NAPL phase, such that the light adsorption characteristics of NAPL did not change as dissolution progressed. Also, none of the methods used here required the use of external calibration chambers. The five visualization approaches evaluated included three methods developed from previously published models, a binary method, and a novel multiple wavelength method that has the advantage of not requiring any assumptions about the intra-pore interface structure between the various phases (sand/water/NAPL). The new multiple wavelength method is also expected to be applicable to any translucent porous media containing two immiscible fluids (e.g., water-air, water-NAPL). Results from the sand-water-PCE system evaluated here showed that the model that assumes wetting media of uniform pore size (Model C of Niemet and Selker, 2001) and the multiple wavelength model with no interface structure assumptions were able to accurately quantify PCE mass reduction during surfactant flushing. The average mass recoveries from these two imaging methods were greater than 95% for domain-average NAPL saturations of approximately 2.6x10(-2), and were approximately 90% during seven cycles of surfactant flushing that sequentially reduced the average NAPL saturation to 7.5x10(-4).

Simplified contaminant source depletion models as analogs of multiphase simulators.

Four simplified dense non-aqueous phase liquid (DNAPL) source depletion models recently introduced in the literature are evaluated for the prediction of long-term effects of source depletion under natural gradient flow. These models are simple in form (a power function equation is an example) but are shown here to serve as mathematical analogs to complex multiphase flow and transport simulators. The spill and subsequent dissolution of DNAPLs was simulated in domains having different hydrologic characteristics (variance of the log conductivity field=0.2, 1 and 3) using the multiphase flow and transport simulator UTCHEM. The dissolution profiles were fitted using four analytical models: the equilibrium streamtube model (ESM), the advection dispersion model (ADM), the power law model (PLM) and the Damkohler number model (DaM). All four models, though very different in their conceptualization, include two basic parameters that describe the mean DNAPL mass and the joint variability in the velocity and DNAPL distributions. The variability parameter was observed to be strongly correlated with the variance of the log conductivity field in the ESM and ADM but weakly correlated in the PLM and DaM. The DaM also includes a third parameter that describes the effect of rate-limited dissolution, but here this parameter was held constant as the numerical simulations were found to be insensitive to local-scale mass transfer. All four models were able to emulate the characteristics of the dissolution profiles generated from the complex numerical simulator, but the one-parameter PLM fits were the poorest, especially for the low heterogeneity case.

Temporal evolution of DNAPL source and contaminant flux distribution: impacts of source mass depletion.

We investigated, using model simulations, the changes occurring in the distribution of dense non-aqueous phase liquid (DNAPL) mass (Sn) within the source zone during depletion through dissolution, and the resulting changes in the contaminant flux distribution (J) at the source control plane (CP). Two numerical codes (ISCO3D and T2VOC) were used to simulate selected scenarios of DNAPL dissolution and transport in three-dimensional, heterogeneous, spatially correlated, random permeability fields with emplaced sources. Data from the model simulations were interpreted based on population statistics (mean, standard deviation, coefficient of variation) and spatial statistics (centroid, second moments, variograms). The mean and standard deviation of the Sn and J distributions decreased with source mass depletion by dissolution. The decrease in mean and standard deviation was proportional for the J distribution resulting in a constant coefficient of variation (CV), while for the Sn distribution, the mean decreased faster than the standard deviation. The spatial distributions exhibited similar behavior as the population distribution, i.e., the CP flux distribution was more stable (defined by temporally constant second moments and range of variograms) than the Sn distribution. These observations appeared to be independent of the heterogeneity of the permeability (k) field (variance of the log permeability field=1 and 2.45), correlation structure (positive vs. negative correlation between the k and Sn domains) and the DNAPL dissolution model (equilibrium vs. rate-limited), for the cases studied. Analysis of data from a flux monitoring field study (Hill Air Force Base, Utah) at a DNAPL source CP before and after source remediation also revealed temporal invariance of the contaminant flux distribution. These modeling and field observations suggest that the temporal evolution of the contaminant flux distribution can be estimated if the initial distribution is known. However, the findings are preliminary and broader implications to sampling strategies for remediation performance assessment need to be evaluated in additional modeling and experimental studies.

Laboratory investigation of flux reduction from dense non-aqueous phase liquid (DNAPL) partial source zone remediation by enhanced dissolution.

This study investigated the benefits of partial removal of dense nonaqueous phase liquid (DNAPL) source zones using enhanced dissolution in eight laboratory scale experiments. The benefits were assessed by characterizing the relationship between reductions in DNAPL mass and the corresponding reduction in contaminant mass flux. Four flushing agents were evaluated in eight controlled laboratory experiments to examine the effects of displacement fluid property contrasts and associated override and underride on contaminant flux reduction (R(j)) vs. mass reduction (R(m)) relationships (R(j)(R(m))): 1) 50% ethanol/50% water (less dense than water), 2) 40% ethyl-lactate/60% water (more dense than water), 3) 18% ethanol/26% ethyl-lactate/56% water (neutrally buoyant), and 4) 2% Tween-80 surfactant (also neutrally buoyant). For each DNAPL architecture evaluated, replicate experiments were conducted where source zone dissolution was conducted with a single flushing event to remove most of the DNAPL from the system, and with multiple shorter-duration floods to determine the path of the R(j)(R(m)) relationship. All of the single-flushing experiments exhibited similar R(j)(R(m)) relationships indicating that override and underride effects associated with cosolvents did not significantly affect the remediation performance of the agents. The R(j)(R(m)) relationship of the multiple injection experiments for the cosolvents with a density contrast with water tended to be less desirable in the sense that there was less R(j) for a given R(m). UTCHEM simulations supported the observations from the laboratory experiments and demonstrated the capability of this model to predict R(j)(R(m)) relationships for non-uniformly distributed NAPL sources.

Evaluation of simplified mass transfer models to simulate the impacts of source zone architecture on nonaqueous phase liquid dissolution in heterogeneous porous media.

Nonaqueous phase liquid (NAPL) dissolution was studied in three-dimensional (3D) heterogeneous experimental aquifers (25.5 cm x 9 cm x 8.5 cm) with two different longitudinal correlation lengths (2.1 cm and 1.1 cm) and initial spill volumes (22.5 ml and 10.5 ml). Spatial and temporal distributions of NAPL during dissolution were measured using magnetic resonance imaging (MRI). At high NAPL spill volume, average effluent concentrations initially increased during dissolution, as NAPL pools transitioned to NAPL ganglia, and then decreased as the total NAPL-water interfacial area decreased over time. Experimental results were used to test six dissolution models: (i and ii) a one-dimensional (1D) model using either specific NAPL-water interfacial area values estimated from MR images at each time step (i.e., 1D quasi-steady state model), or an empirical mass transfer (Sh') correlation (i.e., 1D transient model), (iii and iv) a multiple analytical source superposition technique (MASST) using either the NAPL distribution determined from MR images at each time step (i.e., MASST steady state model), or the NAPL distribution determined from mass balance calculations (i.e., MASST transient model), (v) an equilibrium streamtube model, and (vi) a 3D grid-scale pool dissolution model (PDM) with a dispersive mass flux term. The 1D quasi-steady state model and 3D PDM captured effluent concentration values most closely, including some concentration fluctuations due to changes in the extent of flow reduction. The 1D transient, MASST steady state and transient, and streamtube models all showed a monotonic decrease in effluent concentration values over time, and the streamtube model was the most computationally efficient. Changes during dissolution of the effective NAPL-water interfacial area estimated from imaging data are similar to changes in effluent concentration values. The 1D steady state model incorporates estimates of the effective NAPL-water interfacial area directly at each time point; the 3D PDM does so indirectly through mass balance and a relative permeability function, which causes reduced water flow through high saturation NAPL regions. Hence, when model accuracy is required, the results indicate that a surrogate of this effective interfacial area is required. Approaches to include this surrogate in the MASST and streamtube models are recommended.

High-resolution reconstruction of the United States human population distribution, 1790 to 2010.

Where do people live, and how has this changed over timescales of centuries? High-resolution spatial information on historical human population distribution is of great significance to understand human-environment interactions and their temporal dynamics. However, the complex relationship between population distribution and various influencing factors coupled with limited data availability make it a challenge to reconstruct human population distribution over timescales of centuries. This study generated 1-km decadal population maps for the conterminous US from 1790 to 2010 using parsimonious models based on natural suitability, socioeconomic desirability, and inhabitability. Five models of increasing complexity were evaluated. The models were validated with census tract and county subdivision population data in 2000 and were applied to generate five sets of 22 historical population maps from 1790-2010. Separating urban and rural areas and excluding non-inhabitable areas were the most important factors for improving the overall accuracy. The generated gridded population datasets and the production and validation methods are described here.

Passive flux meter measurement of water and nutrient flux in saturated porous media: bench-scale laboratory tests.

The passive nutrient flux meter (PNFM) is introduced for simultaneous measurement of both water and nutrient flux through saturated porous media. The PNFM comprises a porous sorbent pre-equilibrated with a suite of alcohol tracers, which have different partitioning coefficients. Water flux was estimated based on the loss of loaded resident tracers during deployment, while nutrient flux was quantified based on the nutrient solute mass captured on the sorbent. An anionic resin, Lewatit 6328 A, was used as a permeable sorbent and phosphate (PO4(3-)) was the nutrient studied. The phosphate sorption capacity of the resin was measured in batch equilibration tests as 56 mg PO4(3-) g(-1), which was determined to be adequate capacity to retain PO4(3-) loads intercepted over typical PNFM deployment periods in most natural systems. The PNFM design was validated with bench-scale laboratory tests for a range of 9.8 to 28.3 cm d(-1) Darcy velocities and 6 to 43 h deployment durations. Nutrient and water fluxes measured by the PNFM averaged within 6 and 12% of the applied values, respectively, indicating that the PNFM shows promise as a tool for simultaneous measurement of water and nutrient fluxes.