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.

Borehole geophysics and horizontal loop electromagnetic profiling (Max-Min) were integrated with regional and site-scale geological and geochemical data to investigate the occurrence of, and possible pathways for, saltwater intrusion near fracture zones on a small island in British Columbia, Canada. An island-wide geochemical study identified a number of coastal wells that are contaminated by seawater; however, the occurrence of high salinity groundwater is spatially irregular due to variable fracturing of the bedrock. To investigate the influence of fracturing on the presence of high salinity groundwaters, geophysical investigations were undertaken at several sites. The nature of the bedrock permeability at these sites, with respect to lithology and fracture zone proximity, is described from geologic and hydrogeologic investigations and supported using surface EM profiling. Fractures and bedding contacts within boreholes, which were suspected to dominate bedrock permeability on the basis of outcrop studies, were identified using borehole video camera in conjunction with normal resistivity, spontaneous potential and natural gamma logs. Flow meter logs, acquired under a variety of aquifer stress conditions including static, tidal and pumping are used to identify potential water transmitting fractures and the locations of entry points for fresh and saline groundwater. The low flow rates measured under the various stress conditions confirm that groundwater flow is minimal and is restricted to mudstone units and single, generally sub-vertical fractures. The low natural groundwater discharge rates near the coast, even at close distances to fracture zones, may be key to the occurrence of saltwater intrusion on many parts of the island.

Soil and groundwater contamination are often managed by establishing on-site cleanup targets within the context of risk assessment or risk management measures. Decision-makers rely on modeling tools to provide insight; however, it is recognized that all models are subject to uncertainty. This case study compares suggested remediation requirements using a site-specific numerical model and a standardized analytical tool to evaluate risk to a downgradient wetland receptor posed by on-site chloride impacts. The base case model, calibrated to observed non-pumping and pumping conditions, predicts a peak concentration well above regulatory criteria. Remediation scenarios are iteratively evaluated to determine a remediation design that adheres to practical site constraints, while minimizing the potential for risk to the downgradient receptor. A nonlinear uncertainty analysis is applied to each remediation scenario to stochastically evaluate the risk and find a solution that meets the site-owner risk tolerance, which in this case required a risk-averse solution. This approach, which couples nonlinear uncertainty analysis with a site-specific numerical model provides an enhanced level of knowledge to foster informed decision-making (i.e., risk-of-success) and also increases stakeholder confidence in the remediation design.

No abstract