ABSTRACT
Subsidence induced by groundwater depletion is a grave problem in many regions around the world, leading to a permanent loss of groundwater storage within an aquifer and even producing structural damage at the Earth's surface. California's Tulare Basin is no exception, experiencing about a meter of subsidence between 2015 and 2020. However, understanding the relationship between changes in groundwater volumes and ground deformation has proven difficult. We employ surface displacement measurements from Interferometric Synthetic Aperture Radar (InSAR) and gravimetric estimates of terrestrial water storage from the Gravity Recovery and Climate Experiment (GRACE) satellite pair to characterize the hydrological dynamics within the Tulare basin. The removal of the long-term aquifer compaction from the InSAR time series reveals coherent short-term variations that correlate with hydrological features. For example, in the winter of 2018-2019 uplift is observed at the confluence of several rivers and streams that drain into the southeastern edge of the basin. These observations, combined with estimates of mass changes obtained from the orbiting GRACE satellites, form the basis for imaging the monthly spatial variations in water volumes. This approach facilitates the quick and effective synthesis of InSAR and gravimetric datasets and will aid efforts to improve our understanding and management of groundwater resources around the world.
ABSTRACT
Range change data, obtained from Synthetic Aperture Radar satellites, form the basis for estimates of aquifer volume change in California's Central Valley. The estimation algorithm incorporates a function penalizing changes far from known well locations, linking the aquifer volume changes to agricultural, industrial, and municipal pumping within the Tulare basin. We show that the range changes are compatible with the hypothesis that the source of aquifer volume changes are variations in effective pressure around documented wells. Specifically, inclusion of the well distance penalty does not degrade the fit to the observations, inversions with and without it both give variance reductions of 99.6%. The patterns of aquifer volume change vary significantly from the drought year, between October 2015 and October 2016, to a wet year in 2017, and into 2018, a year with near average rainfall. The 2.3 million acre-feet of estimated volume reduction, a lower bound on the amount of water extracted from the basin between October 2015 and 2016, agrees with independent estimates of 1.8 and 2.3 million acre-feet. The aquifer volume reduction is also compatible with a loss of 3.1 km3 (2.5 million acre-feet) in groundwater volume derived from Gravity Recovery and Climate Experiment (GRACE) satellite data.
ABSTRACT
The macroscopic laws controlling the advection and diffusion of solute at the scale of the porous continuum are derived in a general manner that does not place limitations on the geometry and time evolution of the pore space. Special focus is given to the definition and symmetry of the dispersion tensor that is controlling how a solute plume spreads out. We show that the dispersion tensor is not symmetric and that the asymmetry derives from the advective derivative in the pore-scale advection-diffusion equation. When flow is spatially variable across a voxel, such as in the presence of a permeability gradient, the amount of asymmetry can be large. As first shown by Auriault [J.-L. Auriault et al. Transp. Porous Med. 85, 771 (2010)TPMEEI0169-391310.1007/s11242-010-9591-y] in the limit of low Péclet number, we show that at any Péclet number, the dispersion tensor D_{ij} satisfies the flow-reversal symmetry D_{ij}(+q)=D_{ji}(-q) where q is the mean flow in the voxel under analysis; however, Reynold's number must be sufficiently small that the flow is reversible when the force driving the flow changes sign. We also demonstrate these symmetries using lattice-Boltzmann simulations and discuss some subtle aspects of how to measure the dispersion tensor numerically. In particular, the numerical experiments demonstrate that the off-diagonal components of the dispersion tensor are antisymmetric which is consistent with the analytical dependence on the average flow gradients that we propose for these off-diagonal components.
