Over the past several decades, the land modeling community has recognized the importance of nutrient regulation on the global terrestrial carbon cycle. Implementations of nutrient limitation in land models are diverse, varying from applying simple empirical down-regulation of potential gross primary productivity under nutrient deficit conditions to more mechanistic treatments. In this study, we introduce a new approach to model multinutrient (nitrogen [N] and phosphorus [P]) limitations in the Energy Exascale Earth System Model (E3SM) Land Model version 1 (ELMv1-ECA). The development is grounded on (1) advances in representing multiple-consumer, multiple-nutrient competition; (2) a generic dynamic allocation scheme based on water, N, P, and light availability; (3) flexible plant CNP stoichiometry; (4) prognostic treatment of N and P constraints on several carbon cycle processes; and (5) global data sets of plant physiological traits. Through benchmarking the model against best knowledge of global plant and soil carbon pools and fluxes, we show that our implementation of nutrient constraints on the present-day carbon cycle is robust at the global scale. Compared with predecessor versions, ELMv1-ECA better predicts global-scale gross primary productivity, ecosystem respiration, leaf area index, vegetation biomass, soil carbon stocks, evapotranspiration, N2O emissions, and NO3- leaching. Factorial experiments indicate that representing the phosphorus cycle improves modeled carbon fluxes, while considering dynamic allocation improves modeled carbon stock density. We also highlight the value of using the International Land Model Benchmarking (ILAMB) package to evaluate and document performance during model development.

Soil erosion delivers enormous amounts of macro-nutrients including nitrogen (N) and phosphorus (P) from land to rivers, potentially sustaining water column bioavailable nutrient levels for decades. In this study, we represent erosional N and P fluxes in the Energy Exascale Earth System Model (E3SM) and apply the model to the continental United States. We estimate that during 1991-2019 soil erosion delivers 775 Gg yr-1 (1 Gg = 109 g) of particulate N (PN) and 328 Gg yr-1 of particulate P (PP) on average to the drainage basins of the northern Gulf of Mexico, including the Mississippi/Atchafalaya River and other rivers draining to the Texas Gulf and the Eastern Gulf. Our model simulation shows that in these rivers PP is the dominant P constituent and over 55% of P exported by erosion comes from soil P pools that could become bioavailable within decades. More importantly, we find that during 1991-2019 erosional N and P fluxes increase at rates of about 15 Gg N yr-1 and 6 Gg P yr-1, respectively, due to increased extreme rains in the Mississippi/Atchafalaya river basin, and this intensification of erosional N and P fluxes drive the significant increase of riverine PN and PP yields to the northern Gulf of Mexico. With extreme rains projected to increase with warming, erosional nutrient fluxes in the region would likely continue to rise in the future, thus complicating the effort of reducing eutrophication in the inland and coastal waters.

We present a representation of nitrogen and phosphorus cycling in the Functionally Assembled Terrestrial Ecosystem Simulator, a demographic vegetation model within the Energy Exascale Earth System land model. This representation is modular, and designed to allow testing of multiple hypothetical approaches for carbon-nutrient coupling in plants. Novel model hypotheses introduced in this work include, (a) the controls on plant acquisition of aqueous mineralized nutrients in the soil and (b) fairly straight forward methods of allocating nutrients to specific plant organs and their losses through live plant turnover as well as litter fluxes generated through plant mortality. This combines the new with pre-existing hypotheses (such as nitrogen fixation and soil decomposition) into a system that can accommodate plant-soil dynamics for a large number of size- and functional-type-resolved plant cohorts within a time-since-disturbance-resolved ecosystem. Root uptake of nutrients is governed by fine root biomass, and plants vary in their fine root biomass allocation in order to balance carbon and nutrient limitations to growth. We test the sensitivity of the model to a wide range of parameter variations and structural representations, and in the context of observations at Barro Colorado Island, Panama. A key model prediction is that plants in the high-light-availability canopy positions allocate more carbon to fine roots than plants in low-light understory environments, given the widely different carbon versus nutrient constraints of these two niches within a given ecosystem. This model provides a basis for exploring carbon-nutrient coupling with vegetation demography within Earth system models.

The maximum photosynthetic carboxylation rate (V_cmax ) is an influential plant trait that has multiple scaling hypotheses, which is a source of uncertainty in predictive understanding of global gross primary production (GPP). Four trait-scaling hypotheses (plant functional type, nutrient limitation, environmental filtering, and plant plasticity) with nine specific implementations were used to predict global V_cmax distributions and their impact on global GPP in the Sheffield Dynamic Global Vegetation Model (SDGVM). Global GPP varied from 108.1 to 128.2 PgC yr^-1 , 65% of the range of a recent model intercomparison of global GPP. The variation in GPP propagated through to a 27% coefficient of variation in net biome productivity (NBP). All hypotheses produced global GPP that was highly correlated (r = 0.85-0.91) with three proxies of global GPP. Plant functional type-based nutrient limitation, underpinned by a core SDGVM hypothesis that plant nitrogen (N) status is inversely related to increasing costs of N acquisition with increasing soil carbon, adequately reproduced global GPP distributions. Further improvement could be achieved with accurate representation of water sensitivity and agriculture in SDGVM. Mismatch between environmental filtering (the most data-driven hypothesis) and GPP suggested that greater effort is needed understand V_cmax variation in the field, particularly in northern latitudes.

Global terrestrial models currently predict that the Amazon rainforest will continue to act as a carbon sink in the future, primarily owing to the rising atmospheric carbon dioxide (CO2) concentration. Soil phosphorus impoverishment in parts of the Amazon basin largely controls its functioning, but the role of phosphorus availability has not been considered in global model ensembles—for example, during the Fifth Climate Model Intercomparison Project. Here we simulate the planned free-air CO2 enrichment experiment AmazonFACE with an ensemble of 14 terrestrial ecosystem models. We show that phosphorus availability reduces the projected CO2-induced biomass carbon growth by about 50% to 79 ± 63 g C m−2 yr−1 over 15 years compared to estimates from carbon and carbon–nitrogen models. Our results suggest that the resilience of the region to climate change may be much less than previously assumed. Variation in the biomass carbon response among the phosphorus-enabled models is considerable, ranging from 5 to 140 g C m−2 yr−1, owing to the contrasting plant phosphorus use and acquisition strategies considered among the models. The Amazon forest response thus depends on the interactions and relative contributions of the phosphorus acquisition and use strategies across individuals, and to what extent these processes can be upregulated under elevated CO2.