Characterization of roof runoff microbial quality in four U.S. cities with varying climate and land use characteristics.

Roof runoff has the potential to serve as an important local water source in regions with growing populations and limited water supply. Given the scarcity of guidance regulating the use of roof runoff, a need exists to characterize the microbial quality of roof runoff. The objective of this 2-year research effort was to examine roof runoff microbial quality in four U.S. cities: Fort Collins, CO; Tucson, AZ; Baltimore, MD; and Miami, FL. Seven participants, i.e., homeowners and schools, were recruited in each city to collect roof runoff samples across 13 precipitation events. Sample collection was done as part of a citizen science approach. The presence and concentrations of indicator organisms and potentially human-infectious pathogens in roof runoff were determined using culture methods and digital droplet polymerase chain reaction (ddPCR), respectively. The analyzed pathogens included Salmonella spp., Campylobacter spp., Giardia duodenalis, and Cryptosporidium parvum. Several factors were evaluated to study their influence on the presence of potentially human-infectious pathogens including the physicochemical characteristics (total suspended solids, volatile suspended solids, total dissolved solids, chemical oxygen demand, and turbidity) of roof runoff, concentrations of indicator organisms, presence/absence of trees, storm properties (rainfall depth and antecedent dry period), percent of impervious cover surrounding each sampling location, seasonality, and geographical location. E. coli and enterococci were detected in 73.4% and 96.2% of the analyzed samples, respectively. Concentrations of both E. coli and enterococci ranged from <0 log10 to >3.38 log10 MPN/100 mL. Salmonella spp. invA, Campylobacter spp. ceuE, and G. duodenalis ß - giardin gene targets were detected in 8.9%, 2.5%, and 5.1% of the analyzed samples, respectively. Campylobacter spp. mapA and C. parvum 18S rRNA gene targets were not detected in any of the analyzed samples. The detection of Salmonella spp. invA was influenced by the geographical location of the sampling site (Chi-square p-value < 0.001) as well as the number of antecedent dry days prior to a rain event (p-value = 0.002, negative correlation). The antecedent dry period was negatively correlated with the occurrence of Campylobacter spp. ceuE as well (p-value = 0.07). On the other hand, the presence of G. duodenalis ß-giardin in roof runoff was positively correlated with rainfall depth (p-value = 0.05). While physicochemical parameters and impervious area were not found to be correlated with the presence/absence of potentially human-infectious pathogens, significant correlations were found between meteorological parameters and the presence/absence of potentially human-infectious pathogens. Additionally, a weak, yet significant positive correlation, was found only between the concentrations of E. coli and those of Giardia duodenalis ß-giardin. This dataset represents the largest-scale study to date of enteric pathogens in U.S. roof runoff collections and will inform treatment targets for different non-potable end uses for roof runoff. However, the dataset is limited by the low percent detection of bacterial and protozoan pathogens, an issue that is likely to persist challenging the characterization of roof runoff microbial quality given sampling limitations related to the volume and number of samples.

Single-component multiphase lattice Boltzmann simulation of free bubble and crevice heterogeneous cavitation nucleation.

This work serves as an important extension of previous work on cavitation simulation [Sukop and Or, Phys. Rev. E 71, 046703 (2005)10.1103/PhysRevE.71.046703]. A modified Shan-Chen single-component multiphase lattice Boltzmann method is used to simulate two different heterogeneous cavitation nucleation mechanisms, the free gas bubble model and the crevice nucleation model. Improvements include the use of a real-gas equation of state, a redefined effective mass function, and the exact difference method forcing scheme. As a result, much larger density ratios, better thermodynamic consistency, and improved numerical accuracy are achieved. In addition, the crevice nucleation model is numerically investigated using the lattice Boltzmann method. The simulations show excellent qualitative and quantitative agreement with the heterogeneous nucleation theories.

Facilitating Integration in Interdisciplinary Research: Lessons from a South Florida Water, Sustainability, and Climate Project.

Interdisciplinary research is increasingly called upon to find solutions to complex sustainability problems, yet co-creating usable knowledge can be challenging. This article offers broad lessons for conducting interdisciplinary science from the South Florida Water, Sustainability, and Climate Project (SFWSC), a 5-year project funded by the U.S. National Science Foundation (NSF). The goal was to develop a holistic decision-making framework to improve understanding of the complex natural-social system of South Florida water allocation and its threats from climate change, including sea level rise, using a water resources optimization model as an integration mechanism. The SFWSC project faced several challenges, including uncertainty with tasks, high task interdependence, and ensuring communication among geographically dispersed members. Our hypothesis was that adaptive techniques would help overcome these challenges and maintain scientific rigor as research evolved. By systematically evaluating the interdisciplinary management approach throughout the project, we learned that integration can be supported by a three-pronged approach: (1) Build a well-defined team and leadership structure for collaboration across geographic distance and disciplines, ensuring adequate coordination funding, encouraging cross-pollination, and allowing team structure to adapt; (2) intentionally design a process and structure for facilitating collaboration, creating mechanisms for routine analysis, and incorporating collaboration tools that foster communication; and (3) support integration within the scientific framework, by using a shared research output, and encouraging team members to adapt when facing unanticipated constraints. These lessons contribute to the international body of knowledge on interdisciplinary research and can assist teams attempting to develop sustainable solutions in complex natural-social systems.

A hydro-economic model of South Florida water resources system.

South Florida's water infrastructure and ecosystems are under pressure from socio-economic growth. Understanding the region's water resources management tradeoffs is essential for developing effective adaptation strategies to cope with emerging challenges such as climate change and sea level rise, which are expected to affect many other regions in the future. We describe a network-based hydro-economic optimization model of the system to investigate the tradeoffs, incorporating the economic value of water in urban and agricultural sectors and economic damages due to urban flooding while also accounting for water supply to sustain fragile ecosystems such as the Everglades and coastal estuaries. Results illustrate that maintaining high reliability of urban water supply under scenarios of reduced water availability (i.e., drier climate conditions) may trigger economic losses to the Everglades Agricultural Area, which will likely become more vulnerable as competition over scarce water resources increases. More pronounced economic losses are expected in urban and agricultural areas when flows to the Everglades are prioritized. Flow targets for coastal estuaries are occasionally exceeded under optimal flow allocations to various demand nodes, indicating that additional storage may be needed to maintain the environmental integrity of the estuarine ecosystems. Wetter climate conditions, on the other hand, generally lead to increased flows throughout the system with positive effects on meeting water demands, although flood mitigation efforts will necessitate additional releases to the estuaries. Strengths and limitations of the hydro-economic model are discussed.

Value of irrigation water usage in South Florida agriculture.

This study estimates economic loss from South Florida croplands when usage of agricultural irrigation water is altered. In South Florida, 78% of the total value of farm products sold is comprised of cropland products. The majority of Florida citrus and sugarcane are produced in the area, and agricultural irrigation was the largest sector of water use in 2010, followed by public water supply. The Florida Department of Environmental Protection announced in December 2012 that traditional sources of fresh groundwater will have difficulty meeting all of the additional demands by 2030. A shortage of water will impose significant damage to the rural and agriculture economy in Florida, which may lead to higher prices and costs for consumers to purchase citrus or other Florida agriculture products. This paper presents a methodology for estimating economic loss when usage of irrigation water is altered, and examines economic values of irrigation water use for South Florida cropland. The efficient allocation of irrigation water across South Florida cropland is also investigated in order to reduce economic cost to the South Florida agricultural sector.

Economic impacts of urban flooding in South Florida: Potential consequences of managing groundwater to prevent salt water intrusion.

High-value urban zones in coastal South Florida are considered particularly vulnerable to salt water intrusion into the groundwater-based, public water supplies caused by sea level rise (SLR) in combination with the low topography, existing high water table, and permeable karst substrate. Managers in the region closely regulate water depths in the extensive South Florida canal network to control closely coupled groundwater levels and thereby reduce the risk of saltwater intrusion into the karst aquifer. Potential SLR adaptation strategies developed by local managers suggest canal and groundwater levels may have to be increased over time to prevent the increased salt water intrusion risk to groundwater resources. However, higher canal and groundwater levels cause the loss of unsaturated zone storage and lead to an increased risk of inland flooding when the recharge from rainfall exceeds the capacity of the unsaturated zone to absorb it and the water table reaches the surface. Consequently, higher canal and groundwater levels are also associated with increased risk of economic losses, especially during the annual wet seasons. To help water managers and urban planners in this region better understand this trade-off, this study models the relationships between flood insurance claims and groundwater levels in Miami-Dade County. Via regression analyses, we relate the incurred number of monthly flood claims in 16 Miami-Dade County watersheds to monthly groundwater levels over the period from 1996 to 2010. We utilize these estimated statistical relationships to further illustrate various monthly flood loss scenarios that could plausibly result, thereby providing an economic quantification of a "too much water" trade-off. Importantly, this understanding is the first of its kind in South Florida and is exceedingly useful for regional-scale hydro-economic optimization models analyzing trade-offs associated with high water levels.

High temporal resolution modeling of the impact of rain, tides, and sea level rise on water table flooding in the Arch Creek basin, Miami-Dade County Florida USA.

Modeling of groundwater levels in a portion of the low-lying coastal Arch Creek basin in northern Miami-Dade County in Southeast Florida USA, which is subject to repetitive flooding, reveals that rain-induced short-term water table rises can be viewed as a primary driver of flooding events under current conditions. Areas below 0.9m North American Vertical Datum (NAVD) elevation are particularly vulnerable and areas below 1.5m NAVD are vulnerable to exceptionally large rainfall events. Long-term water table rise is evident in the groundwater data, and the rate appears to be consistent with local rates of sea level rise. Linear extrapolation of long-term observed groundwater levels to 2060 suggest roughly a doubling of the number of days when groundwater levels exceed 0.9m NAVD and a threefold increase in the number of days when levels exceed 1.5m NAVD. Projected sea level rise of 0.61m by 2060 together with increased rainfall lead to a model prediction of frequent groundwater-related flooding in areas<0.9m NAVD. However, current simulations do not consider the range of rainfall events that have led to water table elevations>1.5m NAVD and widespread flooding of the area in the past. Tidal fluctuations in the water table are predicted to be more pronounced within 600m of a tidally influenced water control structure that is hydrodynamically connected to Biscayne Bay. The inland influence of tidal fluctuations appears to increase with increased sea level, but the principal driver of high groundwater levels under the 2060 scenario conditions remains groundwater recharge due to rainfall events.

Developing Benthic Class Specific, Chlorophyll-a Retrieving Algorithms for Optically-Shallow Water Using SeaWiFS.

This study evaluated the ability to improve Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) chl-a retrieval from optically shallow coastal waters by applying algorithms specific to the pixels' benthic class. The form of the Ocean Color (OC) algorithm was assumed for this study. The operational atmospheric correction producing Level 2 SeaWiFS data was retained since the focus of this study was on establishing the benefit from the alternative specification of the bio-optical algorithm. Benthic class was determined through satellite image-based classification methods. Accuracy of the chl-a algorithms evaluated was determined through comparison with coincident in situ measurements of chl-a. The regionally-tuned models that were allowed to vary by benthic class produced more accurate estimates of chl-a than the single, unified regionally-tuned model. Mean absolute percent difference was approximately 70% for the regionally-tuned, benthic class-specific algorithms. Evaluation of the residuals indicated the potential for further improvement to chl-a estimation through finer characterization of benthic environments. Atmospheric correction procedures specialized to coastal environments were recognized as areas for future improvement as these procedures would improve both classification and algorithm tuning.

Improved lattice Boltzmann modeling of binary flow based on the conservative Allen-Cahn equation.

The primary and key task of binary fluid flow modeling is to track the interface with good accuracy, which is usually challenging due to the sharp-interface limit and numerical dispersion. This article concentrates on further development of the conservative Allen-Cahn equation (ACE) [Geier et al., Phys. Rev. E 91, 063309 (2015)10.1103/PhysRevE.91.063309] under the framework of the lattice Boltzmann method (LBM), with incorporation of the incompressible hydrodynamic equations [Liang et al., Phys. Rev. E 89, 053320 (2014)10.1103/PhysRevE.89.053320]. Utilizing a modified equilibrium distribution function and an additional source term, this model is capable of correctly recovering the conservative ACE through the Chapman-Enskog analysis. We also simulate four phase-tracking benchmark cases, including one three-dimensional case; all show good accuracy as well as low numerical dispersion. By coupling the incompressible hydrodynamic equations, we also simulate layered Poiseuille flow and the Rayleigh-Taylor instability, illustrating satisfying performance in dealing with complex flow problems, e.g., high viscosity ratio, high density ratio, and high Reynolds number situations. The present work provides a reliable and efficient solution for binary flow modeling.

Granular encapsulation of light hydrophobic liquids (LHL) in LHL-salt water systems: Particle induced densification with quartz sand.

Addition of granular materials to floating crude oil slicks can be effective in capturing and densifying the floating hydrophobic phase, which settles by gravity. Interaction of light hydrophobic liquids (LHL) with quartz sand was investigated in LHL-salt water systems. The LHLs studied were decane, tetradecane, hexadecane, benzene, toluene, ethylbenzene, m-xylene, and 2-cholorotoluene. Experiments were conducted with fine quartz sand (passing sieve No. 40 with openings 0.425 mm). Each LHL was dyed with few crystals of Sudan IV dye for ease of visual observation. A volume of 0.5 mL of each LHL was added to 100 mL salt water (34 g/L). Addition of one gram of quartz sand to the floating hydrophobic liquid layer resulted in formation of sand-encapsulated globules, which settled due to increased density. All LHLs (except for a few globules of decane) formed globules covered with fine sand particles that were heavy enough to settle by gravity. The encapsulated globules were stable and retained their shape upon settling. Polarity of hydrophobic liquids as the main factor of aggregation with minerals was found to be insufficient to explain LHL aggregation with sand. Contact angle measurements were made by submerging a large quartz crystal with the LHL drop on its surface into salt water. A positive correlation was observed between the wetting angle of LHL and the LHL volume captured (r = 0.75). The dependence of the globule density on globule radius was analyzed in relation to the coverage (%) of globule surface (LHL-salt water interface) by fine quartz particles.

Geostatistical Borehole Image-Based Mapping of Karst-Carbonate Aquifer Pores.

Quantification of the character and spatial distribution of porosity in carbonate aquifers is important as input into computer models used in the calculation of intrinsic permeability and for next-generation, high-resolution groundwater flow simulations. Digital, optical, borehole-wall image data from three closely spaced boreholes in the karst-carbonate Biscayne aquifer in southeastern Florida are used in geostatistical experiments to assess the capabilities of various methods to create realistic two-dimensional models of vuggy megaporosity and matrix-porosity distribution in the limestone that composes the aquifer. When the borehole image data alone were used as the model training image, multiple-point geostatistics failed to detect the known spatial autocorrelation of vuggy megaporosity and matrix porosity among the three boreholes, which were only 10 m apart. Variogram analysis and subsequent Gaussian simulation produced results that showed a realistic conceptualization of horizontal continuity of strata dominated by vuggy megaporosity and matrix porosity among the three boreholes.

Modeling the impact of restoration efforts on phosphorus loading and transport through Everglades National Park, FL, USA.

Ecosystems of Florida Everglades are highly sensitive to phosphorus loading. Future restoration efforts, which focus on restoring Everglades water flows, may pose a threat to the health of these ecosystems. To determine the fate and transport of total phosphorus and evaluate proposed Everglades restoration, a water quality model has been developed using the hydrodynamic results from the M3ENP (Mike Marsh Model of Everglades National Park)--a physically-based hydrological numerical model which uses MIKE SHE/MIKE 11 software. Using advection-dispersion with reactive transport for the model, parameters were optimized and phosphorus loading in the overland water column was modeled with good accuracy (60%). The calibrated M3ENP-AD model was then modified to include future bridge construction and canal water level changes, which have shown to increase flows into ENP. These bridge additions increased total dissolved phosphorus (TP) load downstream in Shark Slough and decreased TP load in downstream Taylor Slough. However, there was a general decrease in TP concentration and TP mass per area over the entire model domain. The M3ENP-AD model has determined the mechanisms for TP transport and quantified the impacts of ENP restoration efforts on the spatial-temporal distribution of phosphorus transport. This tool can be used to guide future Everglades restoration decisions.

Quantifying data worth toward reducing predictive uncertainty.

The present study demonstrates a methodology for optimization of environmental data acquisition. Based on the premise that the worth of data increases in proportion to its ability to reduce the uncertainty of key model predictions, the methodology can be used to compare the worth of different data types, gathered at different locations within study areas of arbitrary complexity. The method is applied to a hypothetical nonlinear, variable density numerical model of salt and heat transport. The relative utilities of temperature and concentration measurements at different locations within the model domain are assessed in terms of their ability to reduce the uncertainty associated with predictions of movement of the salt water interface in response to a decrease in fresh water recharge. In order to test the sensitivity of the method to nonlinear model behavior, analyses were repeated for multiple realizations of system properties. Rankings of observation worth were similar for all realizations, indicating robust performance of the methodology when employed in conjunction with a highly nonlinear model. The analysis showed that while concentration and temperature measurements can both aid in the prediction of interface movement, concentration measurements, especially when taken in proximity to the interface at locations where the interface is expected to move, are of greater worth than temperature measurements. Nevertheless, it was also demonstrated that pairs of temperature measurements, taken in strategic locations with respect to the interface, can also lead to more precise predictions of interface movement.

Solute and heat transport model of the Henry and hilleke laboratory experiment.

SEAWAT is a coupled version of MODFLOW and MT3DMS designed to simulate variable-density ground water flow and solute transport. The most recent version of SEAWAT, called SEAWAT Version 4, includes new capabilities to represent simultaneous multispecies solute and heat transport. To test the new features in SEAWAT, the laboratory experiment of Henry and Hilleke (1972) was simulated. Henry and Hilleke used warm fresh water to recharge a large sand-filled glass tank. A cold salt water boundary was represented on one side. Adjustable heating pads were used to heat the bottom and left sides of the tank. In the laboratory experiment, Henry and Hilleke observed both salt water and fresh water flow systems separated by a narrow transition zone. After minor tuning of several input parameters with a parameter estimation program, results from the SEAWAT simulation show good agreement with the experiment. SEAWAT results suggest that heat loss to the room was more than expected by Henry and Hilleke, and that multiple thermal convection cells are the likely cause of the widened transition zone near the hot end of the tank. Other computer programs with similar capabilities may benefit from benchmark testing with the Henry and Hilleke laboratory experiment.

Lattice Boltzmann models for flow and transport in saturated karst.

Flow and transport simulation in karst aquifers remains a significant challenge for the ground water modeling community. Darcy's law-based models cannot simulate the inertial flows characteristic of many karst aquifers. Eddies in these flows can strongly affect solute transport. The simple two-region conduit/matrix paradigm is inadequate for many purposes because it considers only a capacitance rather than a physical domain. Relatively new lattice Boltzmann methods (LBMs) are capable of solving inertial flows and associated solute transport in geometrically complex domains involving karst conduits and heterogeneous matrix rock. LBMs for flow and transport in heterogeneous porous media, which are needed to make the models applicable to large-scale problems, are still under development. Here we explore aspects of these future LBMs, present simple examples illustrating some of the processes that can be simulated, and compare the results with available analytical solutions. Simulations are contrived to mimic simple capacitance-based two-region models involving conduit (mobile) and matrix (immobile) regions and are compared against the analytical solution. There is a high correlation between LBM simulations and the analytical solution for two different mobile region fractions. In more realistic conduit/matrix simulation, the breakthrough curve showed classic features and the two-region model fit slightly better than the advection-dispersion equation (ADE). An LBM-based anisotropic dispersion solver is applied to simulate breakthrough curves from a heterogeneous porous medium, which fit the ADE solution. Finally, breakthrough from a karst-like system consisting of a conduit with inertial regime flow in a heterogeneous aquifer is compared with the advection-dispersion and two-region analytical solutions.

Distribution of multiphase fluids in porous media: comparison between lattice Boltzmann modeling and micro-x-ray tomography.

A parallel implementation of the three-dimensional Shan-and-Chen multicomponent, multiphase lattice Boltzmann method (LBM) was used to simulate the equilibrium distributions of two immiscible fluids in porous media. The simulations were successfully validated against cone-beam x-ray microtomographic data on the distribution of oil (decane), water, and air phases in a 5-mm cube of porous medium composed of packed quartz sand grains. The results confirm that LBM models allow for the straightforward incorporation of complex pore space geometry determined from x-ray microtomography measurements and that simulated wetting and nonwetting phase distributions are consistent with x-ray observations on both macroscopic and microscopic scales.

Proposed approximation for contact angles in Shan-and-Chen-type multicomponent multiphase lattice Boltzmann models.

We propose a method for approximating the adhesion parameters in the Shan and Chen multicomponent, multiphase lattice Boltzmann model that leads to the desired fluid-solid contact angle. The method is a straightforward application of Young's equation with substitution of the Shan and Chen cohesion parameter and a density factor for the fluid-fluid interfacial tension, and the adhesion parameters for the corresponding fluid-solid interfacial tensions.

Simulation of gaseous diffusion in partially saturated porous media under variable gravity with lattice Boltzmann methods.

Liquid distributions in unsaturated porous media under different gravitational accelerations and corresponding macroscopic gaseous diffusion coefficients were investigated to enhance understanding of plant growth conditions in microgravity. We used a single-component, multiphase lattice Boltzmann code to simulate liquid configurations in two-dimensional porous media at varying water contents for different gravity conditions and measured gas diffusion through the media using a multicomponent lattice Boltzmann code. The relative diffusion coefficients (D rel) for simulations with and without gravity as functions of air-filled porosity were in good agreement with measured data and established models. We found significant differences in liquid configuration in porous media, leading to reductions in D rel of up to 25% under zero gravity. The study highlights potential applications of the lattice Boltzmann method for rapid and cost-effective evaluation of alternative plant growth media designs under variable gravity.

Lattice Boltzmann method for homogeneous and heterogeneous cavitation.

The Shan and Chen single component multiphase lattice Boltzmann model is used to simulate homogeneous and heterogeneous cavitation. The simulations show excellent agreement with the equation of state for homogeneous cavitation and with energy considerations based on interface formation and bubble expansion for heterogeneous cavitation.