Search

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. Plain Language Summary: This work introduces a new set of nutrient cycling hypotheses incorporated into a terrestrial biosphere model. This includes the cycling of carbon, nitrogen and phosphorus, and focuses mainly on plant acquisition, allocation, and turnover. An analysis shows the model offers reasonable responses to perturbations in parameter constants and variability in climate forcing, considering its design balance between process complexity and parameterization burden. Key Points: The nutrient enabled ELM‐FATES model represents expected pattern responses to nutrient availability and parameter perturbationsThe model has been designed to introduce a reasonably small parameterization burdenThese hypothesis can capture some but not all elements of carbon‐nutrient dynamics and can be further inter‐compared with other hypotheses [ABSTRACT FROM AUTHOR]

Topographic heterogeneity and lateral subsurface flow at the hillslope scale of ≤1 km may have outsized impacts on tropical forest through their impacts on water available to plants under water-stressed conditions. However, vegetation dynamics and finer-scale hydrologic processes are not concurrently represented in Earth system models. In this study, we integrate the Energy Exascale Earth System Model (E3SM) land model (ELM) that includes the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), with a three-dimensional hydrology model (ParFlow) to explicitly resolve hillslope topography and subsurface flow and perform numerical experiments to understand how hillslope-scale hydrologic processes modulate vegetation along water availability gradients at Barro Colorado Island (BCI), Panama. Our simulations show that groundwater table depth (WTD) can play a large role in governing aboveground biomass (AGB) when drought-induced tree mortality is triggered by hydraulic failure. Analyzing the simulations using random forest (RF) models, we find that the domain-wide simulated AGB and WTD can be well predicted by static topographic attributes, including surface elevation, slope, and convexity, and adding soil moisture or groundwater table depth as predictors further improves the RF models. Different model representations of mortality due to hydraulic failure can change the dominant topographic driver for the simulated AGB. Contrary to the simulations, the observed AGB in the well-drained 50 ha forest census plot within BCI cannot be well predicted by the RF models using topographic attributes and observed soil moisture as predictors, suggesting other factors such as nutrient status may have a larger influence on the observed AGB. The new coupled model may be useful for understanding the diverse impact of local heterogeneity by isolating the water availability and nutrient availability from the other external and internal factors in ecosystem modeling. [ABSTRACT FROM AUTHOR]