PALADYN v1.0, a comprehensive land surface-vegetation-carbon cycle model of intermediate complexity

PALADYN is presented; it is a new comprehensive and computationally efficient land surface-vegetation-carbon cycle model designed to be used in Earth system models of intermediate complexity for long-term simulations and paleoclimate studies. The model treats in a consistent manner the interaction between atmosphere, terrestrial vegetation and soil through the fluxes of energy, water and carbon. Energy, water and carbon are conserved. PALADYN explicitly treats permafrost, both in physical processes and as an important carbon pool. It distinguishes nine surface types: five different vegetation types, bare soil, land ice, lake and ocean shelf. Including the ocean shelf allows the treatment of continuous changes in sea level and shelf area associated with glacial cycles. Over each surface type, the model solves the surface energy balance and computes the fluxes of sensible, latent and ground heat and upward shortwave and longwave radiation. The model includes a single snow layer. Vegetation and bare soil share a single soil column. The soil is vertically discretized into five layers where prognostic equations for temperature, water and carbon are consistently solved. Phase changes of water in the soil are explicitly considered. A surface hydrology module computes precipitation interception by vegetation, surface runoff and soil infiltration. The soil water equation is based on Darcy's law. Given soil water content, the wetland fraction is computed based on a topographic index. The temperature profile is also computed in the upper part of ice sheets and in the ocean shelf soil. Photosynthesis is computed using a light use efficiency model. Carbon assimilation by vegetation is coupled to the transpiration of water through stomatal conductance. PALADYN includes a dynamic vegetation module with five plant functional types competing for the grid cell share with their respective net primary productivity. PALADYN distinguishes between mineral soil carbon, peat carbon, buried carbon and shelf carbon. Each soil carbon type has its own soil carbon pools generally represented by a litter, a fast and a slow carbon pool in each soil layer. Carbon can be redistributed between the layers by vertical diffusion and advection. For the vegetated macro surface type, decomposition is a function of soil temperature and soil moisture. Carbon in permanently frozen layers is assigned a long turnover time which effectively locks carbon in permafrost. Carbon buried below ice sheets and on flooded ocean shelves is treated differently. The model also includes a dynamic peat module. PALADYN includes carbon isotopes C-13 and C-14, which are tracked through all carbon pools. Isotopic discrimination is modelled only during photosynthesis. A simple methane module is implemented to represent methane emissions from anaerobic carbon decomposition in wetlands (including peatlands) and flooded ocean shelf. The model description is accompanied by a thorough model evaluation in offline mode for the present day and the historical period.