Dynamic spatio-temporal patterns of metapopulation occupancy in patchy habitats.

Spatio-temporal dynamics in habitat suitability and connectivity among mosaics of heterogeneous wetlands are critical for biological diversity and species persistence in aquatic patchy landscapes. Despite the recognized importance of stochastic hydroclimatic forcing in driving wetlandscape hydrological dynamics, linking such effects to emergent dynamics of metapopulation poses significant challenges. To fill this gap, we propose here a dynamic stochastic patch occupancy model (SPOM), which links parsimonious hydrological and ecological models to simulate spatio-temporal patterns in species occupancy in wetlandscapes. Our work aims to place ecological studies of patchy habitats into a proper hydrologic and climatic framework to improve the knowledge about metapopulation shifts in response to climate-driven changes in wetlandscapes. We applied the dynamic version of the SPOM (D-SPOM) framework in two wetlandscapes in the US with contrasting landscape and climate properties. Our results illustrate that explicit consideration of the temporal dimension proposed in the D-SPOM is important to interpret local- and landscape-scale patterns of habitat suitability and metapopulation occupancy. Our analyses show that spatio-temporal dynamics of patch suitability and accessibility, driven by the stochasticity in hydroclimatic forcing, influence metapopulation occupancy and the topological metrics of the emergent wetlandscape dispersal network. D-SPOM simulations also reveal that the extinction risk in dynamic wetlandscapes is exacerbated by extended dry periods when suitable habitat decreases, hence limiting successful patch colonization and exacerbating metapopulation extinction risks. The proposed framework is not restricted only to wetland studies but could also be applied to examine metapopulation dynamics in other types of patchy habitats subjected to stochastic external disturbances.

Strong hydroclimatic controls on vulnerability to subsurface nitrate contamination across Europe.

Subsurface contamination due to excessive nutrient surpluses is a persistent and widespread problem in agricultural areas across Europe. The vulnerability of a particular location to pollution from reactive solutes, such as nitrate, is determined by the interplay between hydrologic transport and biogeochemical transformations. Current studies on the controls of subsurface vulnerability do not consider the transient behaviour of transport dynamics in the root zone. Here, using state-of-the-art hydrologic simulations driven by observed hydroclimatic forcing, we demonstrate the strong spatiotemporal heterogeneity of hydrologic transport dynamics and reveal that these dynamics are primarily controlled by the hydroclimatic gradient of the aridity index across Europe. Contrasting the space-time dynamics of transport times with reactive timescales of denitrification in soil indicate that ~75% of the cultivated areas across Europe are potentially vulnerable to nitrate leaching for at least one-third of the year. We find that neglecting the transient nature of transport and reaction timescale results in a great underestimation of the extent of vulnerable regions by almost 50%. Therefore, future vulnerability and risk assessment studies must account for the transient behaviour of transport and biogeochemical transformation processes.

Convergence of DNAPL source strength functions with site age.

Dissolution of dense nonaqueous phase liquid (DNAPL) source zones can be accurately predicted based on appropriate characterization of the source zone architecture, which controls the rate of mass discharge or source strength function. However, the architecture changes temporally as the source zone mass is depleted by dissolution. To generalize comparisons between contaminated sites with different porewater velocities or contaminant solubilities, site age is defined in terms of the fraction of contaminant mass that has been eluted from the source zone by aqueous dissolution. Here changes in DNAPL architecture during dissolution of a source zone were measured by light transmission visualization in laboratory flow chambers. Architectures measured at ages corresponding to initial conditions, 20, 50, and 90% mass removal were used in an equilibrium streamtube (EST) model to accurately predict subsequent dissolution. It is shown both experimentally and theoretically that as DNAPL contaminated sites age, fractional reductions in contaminant discharge and mass converge to become equal, regardless of the initial architecture. This behavior is a consequence of convergence from log-normal to exponential behavior. Analysis of errors in dissolution predictions suggests that the age of many contaminated sites is likely sufficient that architecture and source strength function characterization may not be necessary as it can be assumed with reasonable accuracy that future dissolution will follow an exponential decay model.

Reactive tracer tests to predict dense nonaqueous phase liquid dissolution dynamics in laboratory flow chambers.

Reactive tracer tests were conducted to evaluate the relationship between contaminant mass reduction, Rm, and flux reduction, Rj, in laboratory experiments with porous media contaminated with a dense nonaqueous phase liquid (DNAPL). The reduction in groundwater contaminant flux resulting from partial mass removal was obtained from continuous and pulsed cosolvent and surfactant flushing dissolution tests in laboratory flow chambers packed with heterogeneous porous media. Using the streamtubes concept a Lagrangian analytical solution was applied to study the contaminant dissolution. The analytical solution was independently parametrized using nonreactive and reactive tracertests and the predicted dissolution was compared to the observed data. Analytical solution parameters related to aquifer hydrodynamic heterogeneities were determined from a nonreactive tracer, while those related to DNAPL spatial distribution heterogeneity were obtained from a reactive tracer. Reactive travel time variance, derived from this combination of tracers, was used to predict the relationship between Rm and Rj. Predictions based on the tracer tests closely matched measured dissolution data, suggesting that tracers can be used to characterize the DNAPL spatial distribution heterogeneity controlling the dissolution behavior. Experimental results demonstrated that increased reactive travel time variance led to greater flux reduction for a given partial mass removal.

Fluid and porous media property effects on dense nonaqueous phase liquid migration and contaminant mass flux.

The effects of fluid and porous media properties on dense nonaqueous phase liquid (DNAPL) migration and associated contaminant mass flux generation were evaluated. Relationships between DNAPL mass and solute mass flux were generated by measuring steady-state mass flux following stepwise injection of perchloroethylene (PCE) into flow chambers packed with homogeneous porous media. The effects of fluid properties including density and interfacial tension (IFT), and media properties including grain size and wettability were evaluated by varying the density contrast and interfacial tension properties between PCE and water, and by varying the porous media mean grain diameter and wettability characteristics. Contaminant mass flux was found to increase as grain size decreased, suggesting enhanced lateral and vertical DNAPL spreading with higher fluid entry pressure. Mass flux showed a slight increase as the DNAPL approached neutral buoyancy, likely due to enhanced vertical spreading above the injection point. DNAPL spatial distribution and contaminant mass flux were only minimally affected by IFT and by intermediate-level wettability changes, but were dramatically affected by wettability reversal. The relationship between DNAPL loading and flux generation became more linear as grain size decreased and density contrast between fluids decreased. These results imply that capillary flow characteristics of the porous media and fluid properties will control mass flux generation from source zones.

Evaluation of remediation performance and cost for field-scale single-phase microemulsion (SPME) flushing.

Pump-and-treat strategies employed to contain contaminant plumes have been demonstrated to be inefficient, and often ineffective, for remediation of aquifers contaminated by nonaqueous phase liquids (NAPLs). This has prompted the development of alternative technologies to provide enhanced remediation of NAPL source zones. The results reported here are from a solubilization study wherein the objective was to achieve high NAPL remediation efficiencies using a low concentration of chemical additives in the flushing solution. A surfactant/alcohol mixture was used to generate a Winsor Type I system, where the NAPL was solubilized and transported as a single-phase microemulsion (SPME). For the SPME process, only 5.5 wt% of the flushing solution comprised chemical additives. The costs associated with this project are discussed and approaches aimed at reducing costs are suggested.