Accounting for forest age in the tile-based dynamic global vegetation model JSBACH4 (4.20p7; git feature/forests) – a land surface model for the ICON-ESM

Abstract

Abstract. Natural and anthropogenic disturbances, in particular forest management, affect forest age structures all around the globe. Forest age structures in turn influence key land surface processes, such as photosynthesis and thus the carbon cycle. Yet, many dynamic global vegetation models (DGVMs), including those used as land surface models (LSMs) in Earth system models (ESMs), do not account for subgrid forest age structures, despite being used to investigate land-use effects on the global carbon budget or simulating biogeochemical responses to climate change. In this paper we present a new scheme to introduce forest age classes in hierarchical tile-based DGVMs combining benefits of recently applied approaches the first being a computationally efficient age-dependent simulation of all relevant processes, such as photosynthesis and respiration, using a restricted number of age classes and the second being the tracking of the exact forest age, which is a prerequisite for any implementation of age-based forest management. This combination is achieved by using the tile hierarchy to track the area fraction for each age on an aggregated plant functional type level, whilst simulating the relevant processes for a set of age classes. We describe how we implemented this scheme in JSBACH4, the LSM of the ICOsahedral Non-hydrostatic Earth system model (ICON-ESM). Subsequently, we compare simulation output to global observation-based products for gross primary production, leaf area index, and above-ground biomass to assess the ability of simulations with and without age classes to reproduce the annual cycle and large-scale spatial patterns of these variables. The comparisons show decreasing differences and increasing computation costs with an increasing number of distinguished age classes. The results demonstrate the benefit of the introduction of age classes, with the optimal number of age classes being a compromise between computation costs and error reduction.