The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and Key Points: • Updated Community Land Model has more hydrological and ecological process fidelity and more comprehensive representation of land management. • The model is systematically evaluated using International Land Model Benchmarking system and shows marked improvement over prior versions. parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5. Plain Language Summary The Community Land Model (CLM) is the land component of the widely used Community Earth System Model (CESM). Here, we introduce model developments included in CLM version 5 (CLM5), the default land component for CESM2 which will be used for the Coupled Model Intercomparison Project (CMIP6). CLM5 includes many new and updated processes including (1) hydrology and snow features such as spatially explicit soil depth, canopy snow processes, a simple firn model, and a more mechanistic river model, (2) plant hydraulics and hydraulic redistribution, (3) revised nitrogen cycling with flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake, (4) expansion to six crop types (global) and time-evolving irrigated areas and fertilization rates, (5) improved urban building energy model, and (6) carbon isotopes. New optional features include a demographically structured dynamic vegetat...

Earth System Models (ESMs) are essential tools for understanding and predicting global change, but they cannot explicitly resolve hillslope-scale terrain structures that fundamentally organize water, energy, and biogeochemical stores and fluxes at subgrid scales. Here we bring together hydrologists, Critical Zone scientists, and ESM developers, to explore how hillslope structures may modulate ESM grid-level water, energy, and biogeochemical fluxes. In contrast to the one-dimensional (1-D), 2-to 3-m deep, and free-draining soil hydrology in most ESM land models, we hypothesize that 3-D, lateral ridge-to-valley flow through shallow and deep paths and insolation contrasts between sunny and shady slopes are the top two globally quantifiable organizers of water and energy (and vegetation) within an ESM grid cell. We hypothesize that these two processes are likely to impact ESM predictions where (and when) water and/or energy are limiting. We further hypothesize that, if implemented in ESM land models, these processes will increase simulated continental water storage and residence time, buffering terrestrial ecosystems against seasonal and interannual droughts. We explore efficient ways to capture these mechanisms in ESMs and identify critical knowledge gaps preventing us from scaling up hillslope to global processes. One such gap is our extremely limited knowledge of the subsurface, where water is stored (supporting vegetation) and released to stream baseflow (supporting aquatic ecosystems). We conclude with a set of organizing

Earth's terrestrial near-subsurface environment can be divided into relatively porous layers of soil, intact regolith, and sedimentary deposits above unweathered bedrock. Variations in the thicknesses of these layers control the hydrologic and biogeochemical responses of landscapes. Currently, Earth System Models approximate the thickness of these relatively permeable layers above bedrock as uniform globally, despite the fact that their thicknesses vary systematically with topography, climate, and geology. To meet the need for more realistic input data for models, we developed a high-resolution gridded global data set of the average thicknesses of soil, intact regolith, and sedimentary deposits within each 30 arcsec (1 km) pixel using the best available data for topography, climate, and geology as input. Our data set partitions the global land surface into upland hillslope, upland valley bottom, and lowland landscape components and uses models optimized for each landform type to estimate the thicknesses of each subsurface layer. On hillslopes, the data set is calibrated and validated using independent data sets of measured soil thicknesses from the U.S. and Europe and on lowlands using depth to bedrock observations from groundwater wells in the U.S. We anticipate that the data set will prove useful as an input to regional and global hydrological and ecosystems models.

No abstract

[1] In this article, a method for drainage network extraction from high-resolution digital elevation models (DEMs; e.g., those derived from airborne laser swath mapping) is presented, which requires just two user-defined parameters and is capable of handling discontinuous valley networks. The accuracy and robustness of the method are illustrated using synthetic valley networks that mimic the complexities of real landscapes and for which the true drainage network is known exactly by construction. The method involves six principal steps: optimal Wiener filtering to remove microtopographic noise, mapping of the contour curvature, identification of valley heads using a user-defined contour-curvature threshold criterion, routing of a unit discharge of water from each valley head using a multiple-flow-direction routing algorithm, removal of discontinuous reaches from the drainage network using a user-defined discharge-per-upstream-valley-head threshold criterion, and thinning of the valley network to a single pixel width. The method yields accurate results using the same user-defined parameters for the two field sites considered in this study, suggesting that for DEMs with resolution of approximately 1 m, the method has the ability to produce accurate results for a variety of landscapes by using the same parameter values used in this study.