The quantification of recharge and trans-valley underflow is needed in arid regions to estimate the impacts of new water withdrawals on the water table. However, for mountainous desert areas, such estimates are highly challenging, due to data scarcity, heterogeneous soils, and long residence times. Conventional assessment employs isolated groundwater models configured with simplified uniform estimates of recharge. Here, we employed a data-constrained surface-subsurface process model to provide an ensemble of spatially distributed recharge and underflow estimates using perturbed parameters. Then, the Model-Independent Parameter Estimation and Uncertainty Quantification (PEST) package was used to calibrate MODFLOW aquifer hydraulic conductivity for this ensemble and reject implausible recharge values. This novel dual-model approach, broadly applicable to mountainous arid regions, was designed to maximally exploit available data sources. It can assimilate groundwater head observations, reject unrealistic parameters, and narrow the range of estimated drawdowns due to pumping. We applied this approach to the Chuckwalla basin in California, USA to determine natural recharge. Simulated recharge concentrates along alluvial fans at the mountain fronts and ephemeral washes where runoff water infiltrates. If an evenly distributed recharge was employed as in conventional studies, it would result in regional biases in estimated drawdown and larger uncertainty bounds. We also note that the speed of groundwater recovery does not guarantee sustainability: heavy pumping induces large hydraulic gradients that initially recover quickly when pumping is halted, but the system may not ultimately recover to pre-pumping levels.

The DOE and BLM identified 285,000 acres of desert land in the Chuckwalla valley in the western U.S., for solar energy development. In addition to several approved solar projects, a pumped storage project was recently proposed to pump nearly 8000 acre-ft-yr of groundwater to store and stabilize solar energy output. This study aims at providing estimates of the amount of naturally-occurring recharge, and to estimate the impact of the pumping on the water table. To better provide the locations and intensity of natural recharge, this study employs an integrated, physically-based hydrologic model, PAWS+CLM, to calculate recharge. Then, the simulated recharge is used in a parameter estimation package to calibrate spatially-distributed K field. This design is to incorporate all available observational data, including soil moisture monitoring stations, groundwater head, and estimates of groundwater conductivity, to constrain the modeling. To address the uncertainty of the soil parameters, an ensemble of simulations are conducted, and the resulting recharges are either rejected or accepted based on calibrated groundwater head and local variation of the K field. The results indicate that the natural total inflow to the study domain is between 7107 and 12,772 afy. During the initial-fill phase of pumped storage project, the total outflow exceeds the upper bound estimate of the inflow. If the initial-fill is annualized to 20 years, the average pumping is more than the lower bound of inflows. The results indicate after adding the pumped storage project, the system will nearing, if not exceeding, its maximum renewable pumping capacity. The accepted recharges lead to a drawdown range of 24 to 45 ft for an assumed specific yield of 0.05. However, the drawdown is sensitive to this parameter, whereas there is insufficient data to adequately constrain this parameter.