Seeing the Canopy for the Branches: Improved Within Canopy Scaling of Leaf Nitrogen

Transitioning across biological scales is a central challenge in land surface models. Processes that operate at the scale of individual leaves must be scaled to canopies, and this is done using dedicated submodels. Here, we focus on a submodel that prescribes how light and nitrogen are distributed through plant canopies. We found a mathematical inconsistency in a submodel implemented in the Community and Energy Land Models (CLM and ELM), which incorporates twigs, branches, stems, and dead leaves in nitrogen scaling from leaf to canopy. The inconsistency leads to unrealistic (physically impossible) values of the nitrogen scaling coefficient. The mathematical inconsistency is a general mistake, that is, would occur in any model adopting this particular submodel. We resolve the inconsistency by allowing distinct profiles of stems and branches versus living leaves. We implemented the updated scheme in the ELM and find that the correction reduces global mean gross primary production (GPP) by 3.9 Pg C (3%). Further, when stems and branches are removed from the canopy in the updated model (akin to models that ignore shading from stems), global GPP increases by 4.1 Pg C (3.2%), because of reduced shading. Hence, models that entirely ignore stem shading also introduce errors in the global spatial distribution of GPP estimates, with a strong signal in the tropics, increasing GPP there by over 200 g C m⁻² yr⁻¹. Appropriately incorporating stems and other nonphotosynthesizing material into the light and nitrogen scaling routines of global land models, will improve their biological realism and accuracy.