Journal cover Journal topic
Hydrology and Earth System Sciences An interactive open-access journal of the European Geosciences Union

Journal metrics

  • IF value: 3.535 IF 3.535
  • IF 5-year<br/> value: 4.292 IF 5-year
  • SNIP value: 1.523 SNIP 1.523
  • IPP value: 3.478 IPP 3.478
  • SJR value: 1.859 SJR 1.859
  • h5-index value: 54 h5-index 54
Hydrol. Earth Syst. Sci. Discuss., 11, 8443-8492, 2014
© Author(s) 2014. This work is distributed
under the Creative Commons Attribution 3.0 License.
Research article
23 Jul 2014
Review status
This discussion paper is under review for the journal Hydrology and Earth System Sciences (HESS).
Resolving terrestrial ecosystem processes along a subgrid topographic gradient for an earth-system model
Z. M. Subin1,2, P. C. D. Milly2,3, B. N. Sulman2,4, S. Malyshev1,2, and E. Shevliakova1,2
1Princeton Environmental Institute, Princeton University, Princeton, New Jersey, USA
2Geophysical Fluid Dynamics Laboratory, National Oceanic and Atmospheric Administration, Princeton, New Jersey, USA
3US Geological Survey, Princeton, New Jersey, USA
4School of Public and Environmental Affairs, Indiana University, Bloomington, Indiana, USA

Abstract. Soil moisture is a crucial control on surface water and energy fluxes, vegetation, and soil carbon cycling. Earth-system models (ESMs) generally represent an areal-average soil-moisture state in gridcells at scales of 50–200 km and as a result are not able to capture the nonlinear effects of topographically-controlled subgrid heterogeneity in soil moisture, in particular where wetlands are present. We addressed this deficiency by building a subgrid representation of hillslope-scale topographic gradients, TiHy (Tiled-hillslope Hydrology), into the Geophysical Fluid Dynamics Laboratory (GFDL) land model (LM3). LM3-TiHy models one or more representative hillslope geometries for each gridcell by discretizing them into land model tiles hydrologically coupled along an upland-to-lowland gradient. Each tile has its own surface fluxes, vegetation, and vertically-resolved state variables for soil physics and biogeochemistry. LM3-TiHy simulates a gradient in soil moisture and water-table depth between uplands and lowlands in each gridcell. Three hillslope hydrological regimes appear in non-permafrost regions in the model: wet and poorly-drained, wet and well-drained, and dry; with large, small, and zero wetland area predicted, respectively. Compared to the untiled LM3 in stand-alone experiments, LM3-TiHy simulates similar surface energy and water fluxes in the gridcell-mean. However, in marginally wet regions around the globe, LM3-TiHy simulates shallow groundwater in lowlands, leading to higher evapotranspiration, lower surface temperature, and higher leaf area compared to uplands in the same gridcells. Moreover, more than four-fold larger soil carbon concentrations are simulated globally in lowlands as compared with uplands. We compared water-table depths to those simulated by a recent global model-observational synthesis, and we compared wetland and inundated areas diagnosed from the model to observational datasets. The comparisons demonstrate that LM3-TiHy has the capability to represent some of the controls of these hydrological variables, but also that improvement in parameterization and input datasets are needed for more realistic simulations. We found large sensitivity in model-diagnosed wetland and inundated area to the depth of conductive soil and the parameterization of macroporosity. With improved parameterization and inclusion of peatland biogeochemical processes, the model could provide a new approach to investigating the vulnerability of Boreal peatland carbon to climate change in ESMs.

Citation: Subin, Z. M., Milly, P. C. D., Sulman, B. N., Malyshev, S., and Shevliakova, E.: Resolving terrestrial ecosystem processes along a subgrid topographic gradient for an earth-system model, Hydrol. Earth Syst. Sci. Discuss., 11, 8443-8492, doi:10.5194/hessd-11-8443-2014, 2014.
Search HESSD
Discussion paper