Groundwater flow, transport, and residence times through topography‐driven basins with exponentially decreasing permeability and porosity

Abstract

[1] In this paper, we investigate the effects of systematic and local heterogeneity on groundwater flow, transport, and residence time distributions (RTDs) of basins where groundwater flow is topography driven. Systematic heterogeneity is represented by an exponentially depth-decreasing hydraulic conductivity and porosity, and local heterogeneity is represented by the dispersivity. The RTDs for both a simple basin with one flow system and a basin with nested local and regional systems gradually evolve to a power law RTD with more pronounced systematic heterogeneity. Exponential decrease of poromechanical properties enhances shallow circulation and subdues deep and regional flows leading to longer flushing times for the large part of the domain, while the shallower portions flush solutes rapidly. Therefore, deeper basins lead to more persistent and pronounced power law RTDs when the poromechanical properties systematically decrease with depth. Separate contributions to the RTD due to stagnation zones associated with local flow cells and due to deeper immobile zones were identified; each leads to a different tailing behavior. Local heterogeneity slightly enhances the power law RTD by causing the tailing to begin earlier but does not affect the late time portion of the RTD. Systematic depth-dependent heterogeneity is an important factor controlling the circulation and associated RTDs of subsurface fluids. It contributes significantly to generation of power law RTDs.Citation: Cardenas, M. B., and X.-W. Jiang (2010), Groundwater flow, transport, and residence times through topography-driven basins with exponentially decreasing permeability and porosity, Water Resour.

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

“…Cardenas and Jiang, 2010). But the implementation of rapid changes in K s with depth in free boundary models with a priori unknown location of water table may cause high computational cost and numerical difficulties (such as numerical instability, dispersion and artificial oscillation) for the particle transport solution, which are compounded as the rate of K s change with depth increases.…”

Legacy Nutrient Dynamics at the Watershed Scale: Principles, Modeling, and Implications

“…For areas having thick unsaturated zones, water and solute movement through the unsaturated zone can also be protracted, resulting in long travel time between the land surface and underlying groundwater (McMahon and B€ ohlke, 2006). The movement of water through subsurface flow systems is better characterized by a distribution of transit times and is typically skewed with long tails, approximating a power-law distribution (Cardenas and Jiang, 2010;Godsey et al, 2010;Stewart et al, 2010). In other words, there is typically longer travel time of dissolved nutrient than that of the groundwater mean transit time.…”

Dépendance à la profondeur et modèles exponentiels de perméabilité dans des cônes de déjection graveleux

“…Vertical permeability variations in each aquifer provide key information for quantifying groundwater flow and solute transport; one example is the typical information aids in realizing geological heterogeneity in geostatistical methods, which often require a fundamental assumption of stationarity (de Marsily 1986;Koltermann and Gorelick 1996;Wackernagel 2010). Analytical and numerical solutions indicate the importance of depth-dependent hydraulic conductivity K in groundwater systems of various scales (Saar and Manga 2004;Marion et al 2008;Cardenas and Jiang 2010;Zlotnik et al 2011). A systematic decreasing trend in either porous media or fractured media can also be understood by using semi-empirical models based on simplifications and well-established relations among permeability, porosity, fracture aperture, and effective stress .…”