Search for: All records

Sanitary sewer systems are critical urban water infrastructure that protect both human and environmental health. Their design, operation, and monitoring require novel modeling techniques that capture dominant processes while allowing for computationally efficient simulations. Open water flow in sewers and rivers are intrinsically similar processes. With this in mind, we formulated a new parsimonious model inspired by the Width Function Instantaneous Unit Hydrograph (WFIUH) approach, widely used to predict rainfallrunoff relationships in watersheds, to a sanitary sewer system consisting of nearly 10,000 sewer conduits and 120,000 residential and commercial sewage connections in Northern Virginia, U.S.A. Model predictions for the three primary components of sanitary flow, including Base Wastewater Flow (BWF), Groundwater Infiltration (GWI), and Runoff Derived Infiltration and Inflow (RDII), compare favorably with the more computationally demanding industry-standard Storm Water Management Model (SWMM). This novel application of the WFIUH modeling framework should support a number of critical water quality endpoints, including (i) sewer hydrograph separation through the quantification of BWF, GWI, and RDII outflows, (ii) evaluation of the impact of new urban developments on sewage flow dynamics, (iii) monitoring and mitigation of sanitary sewer overflows, and (iv) design and interpretation of wastewater surveillance studies. 

Abstract. Groundwater table dynamics extensively modify the volume of the hyporheic zoneand the rate of hyporheic exchange processes. Understanding the effects ofdaily groundwater table fluctuations on the tightly coupled flow and heattransport within hyporheic zones is crucial for water resourcesmanagement. With this aim in mind, a physically based model is used to explorehyporheic responses to varying groundwater table fluctuationscenarios. The effects of different timing and amplitude of groundwater tabledaily drawdowns under gaining and losing conditions are explored in hyporheiczones influenced by natural flood events and diel river temperaturefluctuations. We find that both diel river temperature fluctuations and dailygroundwater table drawdowns play important roles in determining thespatiotemporal variability of hyporheic exchange rates, temperature ofexfiltrating hyporheic fluxes, mean residence times, and hyporheicdenitrification potentials. Groundwater table dynamics present substantiallydistinct impacts on hyporheic exchange under gaining or losing conditions. Thetiming of groundwater table drawdown has a direct influence on hyporheicexchange rates and hyporheic buffering capacity on thermaldisturbances. Consequently, the selection of aquifer pumping regimes hassignificant impacts on the dispersal of pollutants in the aquifer and thermalheterogeneity in the sediment. 

Abstract Analytical solutions for the three‐dimensional groundwater flow equation have been widely used to gain insight about subsurface flow structure and as an alternative to computationally expensive numerical models. Of particular interest are solutions that decompose prescribed hydraulic head boundaries (e.g., Dirichlet boundary condition) into a collection of harmonic functions. Previous studies estimate the frequencies and amplitudes of these harmonics with a least‐square approach where the amplitudes are fitted given a pre‐assigned set of frequencies. In these studies, an ad hoc and structured discretization of the frequency domain is typically used, excluding dominant frequencies while assigning importance to spurious frequencies, with significant consequences for estimating the fluxes and residence times. This study demonstrates the advantages of using a pre‐assigned frequency spectrum that targets the dominant frequencies based on rigorous statistical analysis with predefined significance levels. The new approach is tested for three hydrologic conceptualizations: (a) a synthetic periodic basin, (b) synthetic bedforms, and (c) a natural mountainous watershed. The performance of the frequency spectrum selection is compared with exact analytical or approximate numerical solutions. We found that the new approach better describes the fluxes and residence times for Dirichlet boundaries with well‐defined characteristics spatial scales (e.g., periodic basins and bedforms). For more complex scenarios, such as natural mountainous watersheds, both pre‐assigned frequency spectrums present similar performance. The spectral solutions presented here can play a central role in developing reduced‐complexity models for assessing regional water and solute fluxes within mountain watersheds and hyporheic zones. 

Abstract In this paper, we develop and validate a rigorous modeling framework, based on Duhamel's Theorem, for the unsteady one‐dimensional vertical transport of a solute across a flat sediment‐water interface (SWI) and through the benthic biolayer of a turbulent stream. The modeling framework is novel in capturing the two‐way coupling between evolving solute concentrations above and below the SWI and in allowing for a depth‐varying diffusivity. Three diffusivity profiles within the sediment (constant, exponentially decaying, and a hybrid model) are evaluated against an extensive set of previously published laboratory measurements of turbulent mass transfer across the SWI. The exponential diffusivity profile best represents experimental observations and its reference diffusivity scales with the permeability Reynolds number, a dimensionless measure of turbulence at the SWI. The depth over which turbulence‐enhanced diffusivity decays is of the order of centimeters and comparable to the thickness of the benthic biolayer. Thus, turbulent mixing across the SWI may serve as a universal transport mechanism, supplying the nutrient and energy fluxes needed to sustain microbial growth, and nutrient processing, in the benthic biolayer of stream and coastal sediments. 

Abstract Physics‐based distributed hydrological models that include groundwater are widely used to understand and predict physical and biogeochemical processes within watersheds. Typically, due to computational limitations, watershed modelers minimize the number of elements used in domain discretization, smoothing or even ignoring critical topographic features. We use an idealized model to investigate the implications of mesh refinement along streams and ridges for modeling three‐dimensional groundwater flow and transport in mountainous watersheds. For varying degrees of topographic complexity level (TCL), which increases with the level of mesh refinement, and geological heterogeneity, we estimate and compare steady state baseflow discharge, mean age, and concentration of subsurface weathering products. Results show that ignoring lower‐order streams or ridges diminishes flow through local flow paths and biases higher the contribution of intermediate and regional flow paths, and biases baseflow older. The magnitude of the bias increases for systems where permeability rapidly decreases with depth and is dominated by shallow flow paths. Based on a simple geochemical model, the concentration of weathering products is less sensitive to the TCL, partially due to the thermodynamic constraints on chemical reactions. Our idealized model also reproduces the observed emergent scaling relationship between the groundwater contribution to streamflow and drainage area, and finds that this scaling relationship is not sensitive to mesh TCL. The bias effects have important implications for the use of hydrological models in the interpretation of environmental tracer data and the prediction of biogeochemical evolution of stream water in mountainous watersheds. 

Abstract Many water quality and ecosystem functions performed by streams occur in the benthic biolayer, the biologically active upper (~5 cm) layer of the streambed. Solute transport through the benthic biolayer is facilitated by bedform pumping, a physical process in which dynamic and static pressure variations over the surface of stationary bedforms (e.g., ripples and dunes) drive flow across the sediment‐water interface. In this paper we derive two predictive modeling frameworks, one advective and the other diffusive, for solute transport through the benthic biolayer by bedform pumping. Both frameworks closely reproduce patterns and rates of bedform pumping previously measured in the laboratory, provided that the diffusion model's dispersion coefficient declines exponentially with depth. They are also functionally equivalent, such that parameter sets inferred from the 2D advective model can be applied to the 1D diffusive model, and vice versa. The functional equivalence and complementary strengths of these two models expand the range of questions that can be answered, for example, by adopting the 2D advective model to study the effects of geomorphic processes (such as bedform adjustments to land use change) on flow‐dependent processes and the 1D diffusive model to study problems where multiple transport mechanisms combine (such as bedform pumping and turbulent diffusion). By unifying 2D advective and 1D diffusive descriptions of bedform pumping, our analytical results provide a straightforward and computationally efficient approach for predicting, and better understanding, solute transport in the benthic biolayer of streams and coastal sediments. 

Abstract Coupled groundwater flow and heat transport within hyporheic zones extensively affect water, energy, and solute exchange with surrounding sediments. The local and cumulative implications of this tightly coupled process strongly depend on characteristics of drivers (i.e., discharge and temperature of the water column) and modulators (i.e., hydraulic and thermal properties of the sediment). With this in mind, we perform a systematic numerical analysis of hyporheic responses to understand how the temporal variability of river discharge and temperature affect flow and heat transport within hyporheic zones. We identify typical time series of river discharge and temperature from gauging stations along the headwater region of Mississippi River Basin, which are characterized by different degrees of flow alteration, to drive a physics‐based model of the hyporheic exchange process. Our modeling results indicate that coupled groundwater flow and heat transport significantly affects the dynamic response of hyporheic zones, resulting in substantial differences in exchange rates and characteristic time scales of hyporheic exchange processes. We also find that the hyporheic zone dampens river temperature fluctuations increasingly with higher frequency of temperature fluctuations. This dampening effect depends on the system transport time scale and characteristics of river discharge and temperature variability. Furthermore, our results reveal that the flow alteration reduces the potential of hyporheic zones to act as a temperature buffer and hinders denitrification within hyporheic zones. These results have significant implications for understanding the drivers of local variability in hyporheic exchange and the implications for the development of thermal refugia and ecosystem functioning in hyporheic zones.