[1] This work investigates the role of terrain and subsurface heterogeneity on the interactions between groundwater dynamics and land surface energy fluxes using idealized simulations. A three-dimensional variably saturated groundwater code (ParFlow) coupled to a land surface model (Common Land Model) is used to account for both vertical and lateral water and pressure movement. This creates a fully integrated approach, coupling overland and subsurface flow while having an explicit representation of the water table and all land surface processes forced by atmospheric data. Because the water table is explicitly represented in these simulations, regions with stronger interaction between water table depth and the land surface energy balance (known as critical zones) can be identified. This study uses simple terrain and geologic configurations to demonstrate the importance of lateral surface and subsurface flows in determining land surface heat and moisture fluxes. Strong correlations are found between the land surface fluxes and water table depth across all cases, including terrain shape, subsurface heterogeneity, vegetation type, and climatological region. Results show that different land forms and subsurface heterogeneities produce very different water table dynamics and land surface flux responses to atmospheric forcing. Subsurface formation and properties have the greatest effect on the coupling between the water table and surface heat and moisture fluxes. Changes in landform and land surface slope also have an effect on these interactions by influencing the fraction of rainfall contributing to overland flow versus infiltration. This directly affects the extent of the critical zone with highest coupling strength along the hillside. Vegetative land cover, as seen in these simulations, has a large effect on the energy balance at the land surface but a small effect on streamflow and water table dynamics and thus a limited impact on the land surface-subsurface interactions. Although climate forcing has a direct effect on water table dynamics and feedbacks to the land surface, in this study it does not overcome that of subsurface heterogeneity and terrain.Citation: Rihani, J. F., R. M. Maxwell, and F. K. Chow (2010), Coupling groundwater and land surface processes: Idealized simulations to identify effects of terrain and subsurface heterogeneity on land surface energy fluxes, Water Resour. Res., 46, W12523,

Abstract. Continental-scale hyper-resolution simulations constitute a grand challenge in characterizing nonlinear feedbacks of states and fluxes of the coupled water, energy, and biogeochemical cycles of terrestrial systems. Tackling this challenge requires advanced coupling and supercomputing technologies for earth system models that are discussed in this study, utilizing the example of the implementation of the newly developed Terrestrial Systems Modeling Platform (TerrSysMP v1.0) on JUQUEEN (IBM Blue Gene/Q) of the Jülich Supercomputing Centre, Germany. The applied coupling strategies rely on the Multiple Program Multiple Data (MPMD) paradigm using the OASIS suite of external couplers, and require memory and load balancing considerations in the exchange of the coupling fields between different component models and the allocation of computational resources, respectively. Using the advanced profiling and tracing tool Scalasca to determine an optimum load balancing leads to a 19 % speedup. In massively parallel supercomputer environments, the coupler OASIS-MCT is recommended, which resolves memory limitations that may be significant in case of very large computational domains and exchange fields as they occur in these specific test cases and in many applications in terrestrial research. However, model I/O and initialization in the petascale range still require major attention, as they constitute true big data challenges in light of future exascale computing resources. Based on a factor-two speedup due to compiler optimizations, a refactored coupling interface using OASIS-MCT and an optimum load balancing, the problem size in a weak scaling study can be increased by a factor of 64 from 512 to 32 768 processes while maintaining parallel efficiencies above 80 % for the component models.

No abstract

The effects of terrain, soil moisture heterogeneity, subsurface properties, and water table dynamics on the development and behavior of the atmospheric boundary layer are studied through a set of idealized numerical experiments. The mesoscale atmospheric model Advanced Regional Prediction System (ARPS) is used to isolate the effects of subsurface heterogeneity, terrain, and soil moisture initialization. The simulations are initialized with detailed soil moisture distributions obtained from offline spin-ups using a coupled surface-subsurface model (ParFlow-CLM). In these idealized simulations, we observe that terrain effects dominate the planetary boundary layer (PBL) development during early morning hours, while the soil moisture signature overcomes that of terrain during the afternoon. Water table and subsurface properties produce a similar effect as that of soil moisture as their signatures (reflected in soil moisture profiles, energy fluxes, and evaporation at the land surface) can also overcome that of terrain during afternoon hours. This is mostly clear for land surface energy fluxes and evaporation at the land surface. We also observe the coupling between water table depth and planetary boundary layer depth in our cases is strongest within wet-to-dry transition zones. This extends the findings of previous studies which demonstrate the subsurface connection to surface energy fluxes is strongest in such transition zones. We investigate how this connection extends into the atmosphere and can affect the structure and development of the convective boundary layer.