Spatial and temporal variability of groundwater recharge in a sandstone aquifer in a semiarid region

Abstract

Abstract. With the aim to understand the spatial and temporal variability of groundwater recharge, a high-resolution, spatially distributed numerical model (MIKE SHE) representing surface water and groundwater was used to simulate responses to precipitation in a 2.16 km2 upland catchment on fractured sandstone near Los Angeles, California. Exceptionally high temporal and spatial resolution was used for this catchment modeling: hourly climate data, a 20 m×20 m grid in the horizontal plane, and 240 numerical layers distributed vertically within the thick vadose zone and in the upper part of the groundwater zone. The finest practical spatial and temporal resolutions were selected to accommodate the large degree of surface and subsurface variability of catchment features. Physical property values for the different lithologies were assigned based on previous on-site investigations, whereas the parameters controlling streamflow and evapotranspiration were derived from calibration to continuous streamflow at the outfall and to average hydraulic heads from 17 wells. Confidence in the calibrated model was enhanced by validation through (i) comparison of simulated average recharge to estimates based on the applications of the chloride mass-balance method to data from the groundwater and vadose zones within and beyond the catchment, (ii) comparison of the water isotope signature (18O and 2H) in shallow groundwater to the variability of isotope signatures for precipitation events over an annual cycle, and (iii) comparison of simulated recharge time series and observed fluctuation of water levels. The average simulated recharge across the catchment for the period 1995–2014 is 16 mm yr−1 (4 % of the average annual precipitation), which is consistent with previous estimates obtained by using the chloride mass balance method (4.2 % of the average precipitation). However, one of the most unexpected results was that local recharge was simulated to vary from 0 to >1000 mm yr−1 due to episodic precipitation and overland runoff effects. This recharge occurs episodically with the major flux events at the bottom of the evapotranspiration zone, as simulated by MIKE SHE and confirmed by the isotope signatures, occurring only at the end of the rainy season. This is the first study that combines MIKE SHE simulations with the analysis of water isotopes in groundwater and rainfall to determine the timing of recharge in a sedimentary bedrock aquifer in a semiarid region. The study advances the understanding of recharge and unsaturated flow processes and enhances our ability to predict the effects of surface and subsurface features on recharge rates. This is crucial in highly heterogeneous contaminated sites because different contaminant source areas have widely varying recharge and, hence, groundwater fluxes impacting their mobility.