This Supplementary Material includes a detailed model description with equations. Section 1 summarises the general structure and vertical discretisation of vegetation and soil, and introduces general parameters (Tab 1). Section 2 describes the canopy processes, such as photosynthesis and stomatal coupling, with parameters in Tab. 2. Section 3 introduces vegetation growth, turnover and dynamics and the corresponding parameters are in Tab. 3. The soil biochemistry is described in Section 4, and its parameters are in Tab. 4. Section 5 describes the implementation of the isotope code, with parameters in Tab. 5. Section 6 describes the radiation scheme, surface energy balance and soil hydrology, with parameters described in Tab. 6. The PFT-specific parameters are listed in Tab. 7. 1 General model structure and discretisation Each gridcell of the model is subdivided into nested tiles, each of which is occupied by one specific vegetation-type, representing a plant functional type (PFT). The number of tiles per gridcell is flexible, making it is easy to implement more/different PFTs in the future. In the model, vegetation is represented by an average individual composed of a range of structural pools (leaves, sapwood, heartwood, coarse roots, fine roots, and fruit), a fast overturning, respiring non-structural pool (labile), as well as a seasonal, non-respiring, and non-structural storage pool (reserve). Tree vegetation types are furthermore characterised by their height (m), diameter (m), and stand density (m −2). Soil biogeochemistry is represented using five organic pools: metabolic (met), structural (str) and and woody (wl) litter, as well as fast (f) and slow (s) overturning soil organic matter. Each of these pools contains carbon (C), nitrogen (N) and phosphorus (P), as well as 13 C, 14 C, and 15 N. The unit of the pools is mol X m −2 for vegetation and mol X m −3 for soil biogeochemical pools, where X represents any of these elements. In addition, the model represents the following soil biogeochemical pools (NH 4 , NO 3 , NO y , N 2 O, N 2 , and PO 4), with equivalent units. The model operates on a half-hourly timescale (denoted as dt). Vegetation processes are assumed to respond to these instantaneous conditions and associated fluxes with a process-specific lag time (τ process mavg , see Tab. 1), representing a form of