Dynamics and Dislodgment from Pore Constrictions of a Trapped Nonwetting Droplet Stimulated by Seismic Waves

Seismic waves affect fluid flow and transport processes in porous media. Therefore, quantitative understanding of the role of seismic waves in subsurface hydrodynamics is important for the development of practical applications and prediction of natural phenomena. We present a theoretical fluid dynamics model to describe how low-frequency elastic waves mobilize isolated droplets trapped in pores by capillary resistance. The ability of the theoretical model to predict the critical mobilization amplitudes (Ac) and the displacement dynamics of the nonwetting droplet are validated against computational fluid dynamics (CFD) simulations. Our theory has the advantage of rapid calculation of Ac for various scenarios. Both theory and CFD simulations show that the Ac increases with increasing wave frequency. The theoretical and computational models agree well in the low-frequency range both in terms of predicting the displacement history of the droplet and its eventual dislodgment, but their results begin to diverge with increasing wave frequency since the Hagen-Poiseuille flow approximation in the model becomes invalid. Relative to a previous viscous seismic model, our model compares more favorably to experimental observations. The model is thus appropriate for predicting trapped nonwetting droplet dynamics in and dislodgement from pore constrictions by low-frequency elastic waves

Analysis of the Effects of Dam Release Properties and Ambient Groundwater Flow on Surface Water‐Groundwater Exchange Over a 100‐km‐Long Reach

Hydroelectric dams often create highly dynamic downstream flows that promote surface water‐groundwater (SW‐GW) interactions including bank storage, the temporary storage of river water in the riverbank. Previous research on SW‐GW exchanges in dammed rivers has primarily been at single study sites, which has limited the understanding of how these exchanges evolve as dam releases travel downstream. This study evaluates how dam releases affect SW‐GW exchange continuously over a 100‐km distance. This is accomplished by longitudinally routing water releases through a synthetic river and modeling bed and bank fluid and solute exchange across transverse transects spaced along the reach. Peak and square dam release hydrograph shapes with three magnitudes (0.5, 1.0, and 1.5 m) were considered. The effect of four ambient groundwater flow conditions (very slightly losing, neutral, and two gaining from the perspective of the river) was evaluated for each dam release scenario. Both types of dam release shapes cause SW‐GW interaction over the entire 100‐km distance, and our results show that square type releases cause bank storage exchange well beyond this distance. Strongly gaining conditions reduce the amount of exchange and allow flushing of river‐sourced solute out of the bank after the dam pulse has passed. Both neutral and losing conditions have larger fluid and solute flux into the bank and limit the amount of solute that returns to the river. Our results support that river corridors downstream of dams have increased river‐aquifer connectivity and that this enhanced connectivity can extend at least 100 km downstream

The Role of Eddies inside Pores in the Transition from Darcy to Forchheimer Flows

We studied the role of intra-pore eddies, from viscous to inertial flows, in modifying continuum-scale flow inside pores. Flow regimes spanning Reynolds Number Re ∼ 0 to 1350 are divided into three zones - one zone follows Darcy flow, and the other two zones describe non-Darcy or Forchheimer flow. During viscous flows, i.e., Re \u3c 1, stationary eddies occupy about 1/5 of the pore volume. Eddies grow when Re \u3e 1, and their growth leads to the deviation from Darcy\u27s law and the emergence of Forchheimer flow manifested as a characteristic reduction in the apparent hydraulic conductivity Ka. The reduction in Ka is due to the narrowing of the flow channel which is a consequence of the growth in eddies. The two zones of Forchheimer flow correspond to the changes in rate of reduction in Ka, which in turn are due to the changes in eddy growth rate. Since the characteristics of Forchheimer flow are specific to pore geometry, our results partly explain why a variety of Forchheimer models are expected and needed for different porous media

Theory for Dynamic Longitudinal Dispersion in Fractures and Rivers with Poiseuille Flow

We present a theory for dynamic longitudinal dispersion coefficient (D) for transport by Poiseuille flow, the foundation for models of many natural systems, such as in fractures or rivers. Our theory describes the mixing and spreading process from molecular diffusion, through anomalous transport, and until Taylor dispersion. D is a sixth order function of fracture aperture (b) or river width (W). The time (T) and length (L) scales that separate preasymptotic and asymptotic dispersive transport behavior are T = b2/(4D m), where Dm is the molecular diffusion coefficient, and L = b4 / 48μDm ∂p / ∂x, where p is pressure and μ is viscosity. In the case of some major rivers, we found that L is ∼150W. Therefore, transport has to occur over a relatively long domain or long time for the classical advection-dispersion equation to be valid

Pore Geometry Effects on Intrapore Viscous to Inertial Flows and on Effective Hydraulic Parameters

In this article, the effects of different diverging-converging pore geometries were investigated, and the microscale fluid flow and effective hydraulic properties from these pores were compared with that of a pipe from viscous to inertial laminar flow regimes. The flow fields are obtained using computational fluid dynamics, and the comparative analysis is based on a new dimensionless hydraulic shape factor β, which is the specific surface scaled by the length of pores. Results from all diverging-converging pores show an inverse pattern in velocity and vorticity distributions relative to the pipe flow. The hydraulic conductivity K of all pores is dependent on and can be predicted from β with a power function with an exponent of 3/2. The differences in K are due to the differences in distribution of local friction drag on the pore walls. At Reynolds number (Re) ∼ 0 flows, viscous eddies are found to exist almost in all pores in different sizes, but not in the pipe. Eddies grow when Re →1 and leads to the failure of Darcy\u27s law. During non-Darcy or Forchheimer flows, the apparent hydraulic conductivity Ka decreases due to the growth of eddies, which constricts the bulk flow region. At Re \u3e 1, the rate of decrease in Ka increases, and at Re \u3e\u3e 1, it decreases to where the change in Ka ≈ 0, and flows once again exhibits a Darcy-type relationship. The degree of nonlinearity during non-Darcy flow decreases for pores with increasing β. The nonlinear flow behavior becomes weaker as β increases to its maximum value in the pipe, which shows no nonlinearity in the flow; in essence, Darcy\u27s law stays valid in the pipe at all laminar flow conditions. The diverging-converging geometry in pores plays a critical role in modifying the intrapore fluid flow, implying that this property should be incorporated in effective larger-scale models, e.g., pore-network models

Ground-Based Thermography of Fluvial Systems at Low and High Discharge Reveals Potential Complex Thermal Heterogeneity Driven by Flow Variation and Bioroughness

Temperature is a primary physical and biogeochemical variable in aquatic systems. Field-based measurement of temperature at discrete sampling points has revealed temperature variability in fluvial systems, but traditional techniques do not readily allow for synoptic sampling schemes that can address temperature-related questions with broad, yet detailed, coverage. We present results of thermal infrared imaging at different stream discharge (base flow and peak flood) conditions using a handheld IR camera. Remotely sensed temperatures compare well with those measured with a digital thermometer. The thermal images show that periphyton, wood, and sandbars induce significant thermal heterogeneity during low stages. Moreover, the images indicate temperature variability within the periphyton community and within the partially submerged bars. The thermal heterogeneity was diminished during flood inundation, when the areas of more slowly moving water to the side of the stream differed in their temperature. The results have consequences for thermally sensitive hydroecological processes and implications for models of those processes, especially those that assume an effective stream temperature

Active Layer Groundwater Flow: The Interrelated Effects of Stratigraphy, Thaw, and Topography

The external drivers and internal controls of groundwater flow in the thawed “active layer” above permafrost are poorly constrained because they are dynamic and spatially variable. Understanding these controls is critical because groundwater can supply solutes such as dissolved organic matter to surface water bodies. We calculated steady‐state three‐dimensional suprapermafrost groundwater flow through the active layer using measurements of aquifer geometry, saturated thickness, and hydraulic properties collected from two major landscape types over time within a first‐order Arctic watershed. The depth position and thickness of the saturated zone is the dominant control of groundwater flow variability between sites and during different times of year. The effect of water table depth on groundwater flow dwarfs the effect of thaw depth. In landscapes with low land‐surface slopes (2–4%), a combination of higher water tables and thicker, permeable peat deposits cause relatively constant groundwater flows between the early and late thawed seasons. Landscapes with larger land‐surface slopes (4–10%) have both deeper water tables and thinner peat deposits; here the commonly observed permeability decrease with depth is more pronounced than in flatter areas, and groundwater flows decrease significantly between early and late summer as the water table drops. Groundwater flows are also affected by microtopographic features that retain groundwater that could otherwise be released as the active layer deepens. The dominant sources of groundwater, and thus dissolved organic matter, are likely wet, flatter regions with thick organic layers. This finding informs fluid flow and solute transport dynamics for the present and future Arctic.Plain Language SummaryGroundwater flow in permafrost watersheds is potentially a key component of global carbon budgets because permafrost soil stores vast amounts of carbon that could be mobilized due to a warming climate and the corresponding increase in soil thaw. In addition to carrying carbon, groundwater can supply important nutrients and solutes to surface waters. However, we do not yet understand the factors that control groundwater flow in soils above permafrost because saturation changes rapidly and continuously, and soil hydraulic properties are largely unknown. We created measurement‐informed calculations of groundwater flow from areas of permafrost with different characteristics and found that soil types, which vary based on the slope of the land surface, are the most important control. Near‐surface soils were identical in hillslopes and valleys, whereas deeper soils in hillslopes allowed for less groundwater flow than in valleys. In early summer, when only the near‐surface soils were thawed, groundwater flows in the hillslopes and valley were similar. In late summer, when the deeper soil was thawed, groundwater flow in the valley remained high, but flow in the hillslope was negligible. Our observations also showed that small mounds on the land surface caused groundwater to be trapped behind underground ice dams.Key PointsDetailed measurements of hydraulic head, hydraulic conductivity, and saturated thicknesses in active layers were made over time and spaceThree main soil layers consistently comprise the stratigraphy of the active layer across the studied Arctic watershedGroundwater flow depends most on the depth of the water table and the subsurface stratigraphy, which varies based on landscape typePeer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/151862/1/wrcr24085_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/151862/2/wrcr24085.pd

Active Layer Groundwater Flow: The Interrelated Effects of Stratigraphy, Thaw, and Topography

The external drivers and internal controls of groundwater flow in the thawed “active layer” above permafrost are poorly constrained because they are dynamic and spatially variable. Understanding these controls is critical because groundwater can supply solutes such as dissolved organic matter to surface water bodies. We calculated steady‐state three‐dimensional suprapermafrost groundwater flow through the active layer using measurements of aquifer geometry, saturated thickness, and hydraulic properties collected from two major landscape types over time within a first‐order Arctic watershed. The depth position and thickness of the saturated zone is the dominant control of groundwater flow variability between sites and during different times of year. The effect of water table depth on groundwater flow dwarfs the effect of thaw depth. In landscapes with low land‐surface slopes (2–4%), a combination of higher water tables and thicker, permeable peat deposits cause relatively constant groundwater flows between the early and late thawed seasons. Landscapes with larger land‐surface slopes (4–10%) have both deeper water tables and thinner peat deposits; here the commonly observed permeability decrease with depth is more pronounced than in flatter areas, and groundwater flows decrease significantly between early and late summer as the water table drops. Groundwater flows are also affected by microtopographic features that retain groundwater that could otherwise be released as the active layer deepens. The dominant sources of groundwater, and thus dissolved organic matter, are likely wet, flatter regions with thick organic layers. This finding informs fluid flow and solute transport dynamics for the present and future Arctic

Dynamics of Hyporheic Flow and Heat Transport Across a Bed-to-Bank Continuum in a Large Regulated River

The lower Colorado River (LCR) near Austin, Texas is heavily regulated for hydropower generation. Daily water releases from a dam located 23 km upstream of our study site in the LCR caused the stage to fluctuate by more than 1.5 m about a mean depth of 1.3 m. As a result, the river switches from gaining to losing over a dam storage-release cycle, driving exchange between river water and groundwater. We assessed the hydrologic impacts of this by simultaneous temperature and head monitoring across a bed-to-bank transect. River-groundwater exchange flux is largest close to the bank and decreases away from the bank. Correspondingly, both the depth of the hyporheic zone and the exchange time are largest close to the bank. Adjacent to the bank, the streambed head response is hysteretic, with the hysteresis disappearing with distance from the bank, indicating that transient bank storage affects the magnitude and direction of vertical exchange close to the bank. Pronounced changes in streambed temperature are observed down to a meter. When the river stage is high, which coincides with when the river is coldest, downward advection of heat from a previous cycles\u27 warm-water pulse warms the streambed. When the river is at its lowest stage but warmest temperature, upwelling groundwater cools the streambed. Future research should consider and focus on a more thorough understanding of the impacts of dam regulation on the hydrologic, thermal, biogeochemical, and ecologic dynamics of rivers and their hyporheic and riparian zones

Groundwater Flow and Exchange Across the Land Surface Explain Carbon Export Patterns in Continuous Permafrost Watersheds

Groundwater flow regimes in the seasonally thawed soils in areas of continuous permafrost are relatively unknown despite their potential role in delivering water, carbon, and nutrients to streams. Using numerical groundwater flow models informed by observations from a headwater catchment in arctic Alaska, United States, we identify several mechanisms that result in substantial surface‐subsurface water exchanges across the land surface during downslope transport and create a primary control on dissolved organic carbon loading to streams and rivers. The models indicate that surface water flowing downslope has a substantial groundwater component due to rapid surface‐subsurface exchanges across a range of hydrologic states, from unsaturated to flooded. Field‐based measurements corroborate the high groundwater contributions, and river dissolved organic carbon concentrations are similar to that of groundwater across large discharge ranges. The persistence of these groundwater contributions in arctic watersheds will influence carbon export to rivers as thaw depth increases in a warmer climate