Journal metrics

Journal metrics

  • IF value: 4.252 IF 4.252
  • IF 5-year value: 4.890 IF 5-year 4.890
  • CiteScore value: 4.49 CiteScore 4.49
  • SNIP value: 1.539 SNIP 1.539
  • SJR value: 2.404 SJR 2.404
  • IPP value: 4.28 IPP 4.28
  • h5-index value: 40 h5-index 40
  • Scimago H index value: 51 Scimago H index 51
Volume 11, issue 10 | Copyright
Geosci. Model Dev., 11, 4085-4102, 2018
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.

Model description paper 11 Oct 2018

Model description paper | 11 Oct 2018

Development and evaluation of a variably saturated flow model in the global E3SM Land Model (ELM) version 1.0

Gautam Bisht1, William J. Riley1, Glenn E. Hammond2, and David M. Lorenzetti3 Gautam Bisht et al.
  • 1Climate & Ecosystem Sciences Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA
  • 2Applied Systems Analysis and Research Department, Sandia National Laboratories, Albuquerque, NM 87185-0747, USA
  • 3Sustainable Energy Systems Group, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA

Abstract. Improving global-scale model representations of near-surface soil moisture and groundwater hydrology is important for accurately simulating terrestrial processes and predicting climate change effects on water resources. Most existing land surface models, including the default E3SM Land Model (ELMv0), which we modify here, routinely employ different formulations for water transport in the vadose and phreatic zones. Clark et al. (2015) identified a variably saturated Richards equation flow model as an important capability for improving simulation of coupled soil moisture and shallow groundwater dynamics. In this work, we developed the Variably Saturated Flow Model (VSFM) in ELMv1 to unify the treatment of soil hydrologic processes in the unsaturated and saturated zones. VSFM was tested on three benchmark problems and results were evaluated against observations and an existing benchmark model (PFLOTRAN). The ELMv1-VSFM's subsurface drainage parameter, fd, was calibrated to match an observationally constrained and spatially explicit global water table depth (WTD) product. Optimal spatially explicit fd values were obtained for 79% of global 1.9° × 2.5° grid cells, while the remaining 21% of global grid cells had predicted WTD deeper than the observationally constrained estimate. Comparison with predictions using the default fd value demonstrated that calibration significantly improved predictions, primarily by allowing much deeper WTDs. Model evaluation using the International Land Model Benchmarking package (ILAMB) showed that improvements in WTD predictions did not degrade model skill for any other metrics. We evaluated the computational performance of the VSFM model and found that the model is about 30% more expensive than the default ELMv0 with an optimal processor layout. The modular software design of VSFM not only provides flexibility to configure the model for a range of problem setups but also allows for building the model independently of the ELM code, thus enabling straightforward testing of the model's physics against other models.

Publications Copernicus
Short summary
Most existing global land surface models used to study impacts of climate change on water resources routinely use different models for near-surface unsaturated soil and the deeper groundwater table. We developed a model that uses a unified treatment of soil hydrologic processes throughout the entire soil column. Using a calibrated drainage parameter, the new model is able to correctly predict deep water table depth as reported in an observationally constrained global dataset.
Most existing global land surface models used to study impacts of climate change on water...