LM4-SHARC v1.0: Resolving the Catchment-scale Soil-Hillslope Aquifer-River Continuum for the GFDL Earth System Modeling Framework

Abstract. Catchment-scale representation of the groundwater and its interaction with other parts of the hydrologic cycle is crucial for accurately depicting the land water-energy balance in Earth system models (ESMs). Despite existing efforts to describe the groundwater in the land component of ESMs, most ESMs still need a prognostic framework for describing catchment-scale groundwater based on its emergent properties to understand its implication to the broader Earth system. To fill this gap, we developed a new parameterization scheme for resolving the groundwater and its two-way interactions with the unsaturated soil and stream at the catchment-scale. We implemented this new parameterization scheme (SHARC, or Soil-Hillslope Aquifer-River Continuum), in the Geophysical Fluid Dynamics Laboratory land model (i.e., LM4-SHARC) and evaluated its performance. By bridging the gap between hydraulic groundwater theory and ESMs' land hydrology, the new LM4-SHARC has been applied to the Providence headwater catchment at Southern Sierra, NV, and tested against in-situ observations. We found that LM4-SHARC leads to noticeable improvements in representing key hydrologic variables such as streamflow, near-surface soil moisture, and soil temperature. In addition to enhancing the representation of the water and energy balance, our analysis showed that accounting for groundwater convergence can induce a more significant hydrologic contrast with higher sensitivity of soil water storage to groundwater properties in the riparian zone. Our findings indicate the feasibility of incorporating two-way interactions among groundwater, unsaturated soil, and streams into the hydrological components of ESMs and further need to explore the implications of these interactions in the context of Earth system dynamics.