Allometric relationships and trade-offs in 11 common Mediterranean-climate grasses.

Biomass allocation in plants is the foundation for understanding dynamics in ecosystem carbon balance, species competition, and plant-environment interactions. However, existing work on plant allometry has mainly focused on trees, with fewer studies having developed allometric equations for grasses. Grasses with different life histories can vary in their carbon investment by prioritizing the growth of specific organs to survive, outcompete co-occurring plants, and ensure population persistence. Further, because grasses are important fuels for wildfire, the lack of grass allocation data adds uncertainty to process-based models that relate plant physiology to wildfire dynamics. To fill this gap, we conducted a greenhouse experiment with 11 common California grasses varying in photosynthetic pathway and growth form. We measured plant sizes and harvested above- and belowground biomass throughout the life cycle of annual species, while for the establishment stage of perennial grasses to quantify allometric relationships for leaf, stem, and root biomass, as well as plant height and canopy area. We used basal diameter as a reference measure of plant size. Overall, basal diameter is the best predictor for leaf and stem biomass, height, and canopy area. Including height as another predictor can improve model accuracy in predicting leaf and stem biomass and canopy area. Fine root biomass is a function of leaf biomass alone. Species vary in their allometric relationships, with most variation occurring for plant height, canopy area, and stem biomass. We further explored potential trade-offs in biomass allocation across species between leaf and fine root, leaf and stem, and allocation to reproduction. Consistent with our expectation, we found that fast-growing plants allocated a greater fraction to reproduction. Additionally, plant height and specific leaf area negatively influenced the leaf-to-stem ratio. However, contrary to our hypothesis, there were no differences in root-to-leaf ratio between perennial and annual or C4 and C3 plants. Our study provides species-specific and functional-type-specific allometry equations for both above- and belowground organs of 11 common California grass species, enabling nondestructive biomass assessment in California grasslands. These allometric relationships and trade-offs in carbon allocation across species can improve ecosystem model predictions of grassland species interactions and environmental responses through differences in morphology.

The effect of the vertical gradients of photosynthetic parameters on the CO2 assimilation and transpiration of a Panamanian tropical forest.

Terrestrial biosphere models (TBMs) include the representation of vertical gradients in leaf traits associated with modeling photosynthesis, respiration, and stomatal conductance. However, model assumptions associated with these gradients have not been tested in complex tropical forest canopies. We compared TBM representation of the vertical gradients of key leaf traits with measurements made in a tropical forest in Panama and then quantified the impact of the observed gradients on simulated canopy-scale CO2 and water fluxes. Comparison between observed and TBM trait gradients showed divergence that impacted canopy-scale simulations of water vapor and CO2 exchange. Notably, the ratio between the dark respiration rate and the maximum carboxylation rate was lower near the ground than at the top-of-canopy, leaf-level water-use efficiency was markedly higher at the top-of-canopy, and the decrease in maximum carboxylation rate from the top-of-canopy to the ground was less than TBM assumptions. The representation of the gradients of leaf traits in TBMs is typically derived from measurements made within-individual plants, or, for some traits, assumed constant due to a lack of experimental data. Our work shows that these assumptions are not representative of the trait gradients observed in species-rich, complex tropical forests.

The drivers and impacts of Amazon forest degradation.

Approximately 2.5 × 106 square kilometers of the Amazon forest are currently degraded by fire, edge effects, timber extraction, and/or extreme drought, representing 38% of all remaining forests in the region. Carbon emissions from this degradation total up to 0.2 petagrams of carbon per year (Pg C year-1), which is equivalent to, if not greater than, the emissions from Amazon deforestation (0.06 to 0.21 Pg C year-1). Amazon forest degradation can reduce dry-season evapotranspiration by up to 34% and cause as much biodiversity loss as deforestation in human-modified landscapes, generating uneven socioeconomic burdens, mainly to forest dwellers. Projections indicate that degradation will remain a dominant source of carbon emissions independent of deforestation rates. Policies to tackle degradation should be integrated with efforts to curb deforestation and complemented with innovative measures addressing the disturbances that degrade the Amazon forest.

Earlier snowmelt may lead to late season declines in plant productivity and carbon sequestration in Arctic tundra ecosystems.

Arctic warming is affecting snow cover and soil hydrology, with consequences for carbon sequestration in tundra ecosystems. The scarcity of observations in the Arctic has limited our understanding of the impact of covarying environmental drivers on the carbon balance of tundra ecosystems. In this study, we address some of these uncertainties through a novel record of 119 site-years of summer data from eddy covariance towers representing dominant tundra vegetation types located on continuous permafrost in the Arctic. Here we found that earlier snowmelt was associated with more tundra net CO2 sequestration and higher gross primary productivity (GPP) only in June and July, but with lower net carbon sequestration and lower GPP in August. Although higher evapotranspiration (ET) can result in soil drying with the progression of the summer, we did not find significantly lower soil moisture with earlier snowmelt, nor evidence that water stress affected GPP in the late growing season. Our results suggest that the expected increased CO2 sequestration arising from Arctic warming and the associated increase in growing season length may not materialize if tundra ecosystems are not able to continue sequestering CO2 later in the season.

Demographic composition, not demographic diversity, predicts biomass and turnover across temperate and tropical forests.

The growth and survival of individual trees determine the physical structure of a forest with important consequences for forest function. However, given the diversity of tree species and forest biomes, quantifying the multitude of demographic strategies within and across forests and the way that they translate into forest structure and function remains a significant challenge. Here, we quantify the demographic rates of 1961 tree species from temperate and tropical forests and evaluate how demographic diversity (DD) and demographic composition (DC) differ across forests, and how these differences in demography relate to species richness, aboveground biomass (AGB), and carbon residence time. We find wide variation in DD and DC across forest plots, patterns that are not explained by species richness or climate variables alone. There is no evidence that DD has an effect on either AGB or carbon residence time. Rather, the DC of forests, specifically the relative abundance of large statured species, predicted both biomass and carbon residence time. Our results demonstrate the distinct DCs of globally distributed forests, reflecting biogeography, recent history, and current plot conditions. Linking the DC of forests to resilience or vulnerability to climate change, will improve the precision and accuracy of predictions of future forest composition, structure, and function.

Warming increased bark beetle-induced tree mortality by 30% during an extreme drought in California.

Quantifying the responses of forest disturbances to climate warming is critical to our understanding of carbon cycles and energy balances of the Earth system. The impact of warming on bark beetle outbreaks is complex as multiple drivers of these events may respond differently to warming. Using a novel model of bark beetle biology and host tree interactions, we assessed how contemporary warming affected western pine beetle (Dendroctonus brevicomis) populations and mortality of its host, ponderosa pine (Pinus ponderosa), during an extreme drought in the Sierra Nevada, California, United States. When compared with the field data, our model captured the western pine beetle flight timing and rates of ponderosa pine mortality observed during the drought. In assessing the influence of temperature on western pine beetles, we found that contemporary warming increased the development rate of the western pine beetle and decreased the overwinter mortality rate of western pine beetle larvae leading to increased population growth during periods of lowered tree defense. We attribute a 29.9% (95% CI: 29.4%-30.2%) increase in ponderosa pine mortality during drought directly to increases in western pine beetle voltinism (i.e., associated with increased development rates of western pine beetle) and, to a much lesser extent, reductions in overwintering mortality. These findings, along with other studies, suggest each degree (°C) increase in temperature may have increased the number of ponderosa pine killed by upwards of 35%-40% °C-1 if the effects of compromised tree defenses (15%-20%) and increased western pine beetle populations (20%) are additive. Due to the warming ability to considerably increase mortality through the mechanism of bark beetle populations, models need to consider climate's influence on both host tree stress and the bark beetle population dynamics when determining future levels of tree mortality.

Hydraulically-vulnerable trees survive on deep-water access during droughts in a tropical forest.

Deep-water access is arguably the most effective, but under-studied, mechanism that plants employ to survive during drought. Vulnerability to embolism and hydraulic safety margins can predict mortality risk at given levels of dehydration, but deep-water access may delay plant dehydration. Here, we tested the role of deep-water access in enabling survival within a diverse tropical forest community in Panama using a novel data-model approach. We inversely estimated the effective rooting depth (ERD, as the average depth of water extraction), for 29 canopy species by linking diameter growth dynamics (1990-2015) to vapor pressure deficit, water potentials in the whole-soil column, and leaf hydraulic vulnerability curves. We validated ERD estimates against existing isotopic data of potential water-access depths. Across species, deeper ERD was associated with higher maximum stem hydraulic conductivity, greater vulnerability to xylem embolism, narrower safety margins, and lower mortality rates during extreme droughts over 35 years (1981-2015) among evergreen species. Species exposure to water stress declined with deeper ERD indicating that trees compensate for water stress-related mortality risk through deep-water access. The role of deep-water access in mitigating mortality of hydraulically-vulnerable trees has important implications for our predictive understanding of forest dynamics under current and future climates.

Spatial heterogeneity and environmental predictors of permafrost region soil organic carbon stocks.

Large stocks of soil organic carbon (SOC) have accumulated in the Northern Hemisphere permafrost region, but their current amounts and future fate remain uncertain. By analyzing dataset combining >2700 soil profiles with environmental variables in a geospatial framework, we generated spatially explicit estimates of permafrost-region SOC stocks, quantified spatial heterogeneity, and identified key environmental predictors. We estimated that Pg C are stored in the top 3 m of permafrost region soils. The greatest uncertainties occurred in circumpolar toe-slope positions and in flat areas of the Tibetan region. We found that soil wetness index and elevation are the dominant topographic controllers and surface air temperature (circumpolar region) and precipitation (Tibetan region) are significant climatic controllers of SOC stocks. Our results provide first high-resolution geospatial assessment of permafrost region SOC stocks and their relationships with environmental factors, which are crucial for modeling the response of permafrost affected soils to changing climate.

A spatial emergent constraint on the sensitivity of soil carbon turnover to global warming.

Carbon cycle feedbacks represent large uncertainties in climate change projections, and the response of soil carbon to climate change contributes the greatest uncertainty to this. Future changes in soil carbon depend on changes in litter and root inputs from plants and especially on reductions in the turnover time of soil carbon (τs) with warming. An approximation to the latter term for the top one metre of soil (ΔCs,τ) can be diagnosed from projections made with the CMIP6 and CMIP5 Earth System Models (ESMs), and is found to span a large range even at 2 °C of global warming (-196 ± 117 PgC). Here, we present a constraint on ΔCs,τ, which makes use of current heterotrophic respiration and the spatial variability of τs inferred from observations. This spatial emergent constraint allows us to halve the uncertainty in ΔCs,τ at 2 °C to -232 ± 52 PgC.

Forest responses to simulated elevated CO2 under alternate hypotheses of size- and age-dependent mortality.

Elevated atmospheric carbon dioxide (eCO2 ) is predicted to increase growth rates of forest trees. The extent to which increased growth translates to changes in biomass is dependent on the turnover time of the carbon, and thus tree mortality rates. Size- or age-dependent mortality combined with increased growth rates could result in either decreased carbon turnover from a speeding up of tree life cycles, or increased biomass from trees reaching larger sizes, respectively. However, most vegetation models currently lack any representation of size- or age-dependent mortality and the effect of eCO2 on changes in biomass and carbon turnover times is thus a major source of uncertainty in predictions of future vegetation dynamics. Using a reduced-complexity form of the vegetation demographic model the Functionally Assembled Terrestrial Ecosystem Simulator to simulate an idealised tropical forest, we find increases in biomass despite reductions in carbon turnover time in both size- and age-dependent mortality scenarios in response to a hypothetical eCO2 -driven 25% increase in woody net primary productivity (wNPP). Carbon turnover times decreased by 9.6% in size-dependent mortality scenarios due to a speeding up of tree life cycles, but also by 2.0% when mortality was age-dependent, as larger crowns led to increased light competition. Increases in aboveground biomass (AGB) were much larger when mortality was age-dependent (24.3%) compared with size-dependent (13.4%) as trees reached larger sizes before death. In simulations with a constant background mortality rate, carbon turnover time decreased by 2.1% and AGB increased by 24.0%, however, absolute values of AGB and carbon turnover were higher than in either size- or age-dependent mortality scenario. The extent to which AGB increases and carbon turnover decreases will thus depend on the mechanisms of large tree mortality: if increased size itself results in elevated mortality rates, then this could reduce by about half the increase in AGB relative to the increase in wNPP.

Species-Specific Shifts in Diurnal Sap Velocity Dynamics and Hysteretic Behavior of Ecophysiological Variables During the 2015-2016 El Niño Event in the Amazon Forest.

Current climate change scenarios indicate warmer temperatures and the potential for more extreme droughts in the tropics, such that a mechanistic understanding of the water cycle from individual trees to landscapes is needed to adequately predict future changes in forest structure and function. In this study, we contrasted physiological responses of tropical trees during a normal dry season with the extreme dry season due to the 2015-2016 El Niño-Southern Oscillation (ENSO) event. We quantified high resolution temporal dynamics of sap velocity (Vs), stomatal conductance (gs) and leaf water potential (ΨL) of multiple canopy trees, and their correlations with leaf temperature (Tleaf) and environmental conditions [direct solar radiation, air temperature (Tair) and vapor pressure deficit (VPD)]. The experiment leveraged canopy access towers to measure adjacent trees at the ZF2 and Tapajós tropical forest research (near the cities of Manaus and Santarém). The temporal difference between the peak of gs (late morning) and the peak of VPD (early afternoon) is one of the major regulators of sap velocity hysteresis patterns. Sap velocity displayed species-specific diurnal hysteresis patterns reflected by changes in Tleaf. In the morning, Tleaf and sap velocity displayed a sigmoidal relationship. In the afternoon, stomatal conductance declined as Tleaf approached a daily peak, allowing ΨL to begin recovery, while sap velocity declined with an exponential relationship with Tleaf. In Manaus, hysteresis indices of the variables Tleaf-Tair and ΨL-Tleaf were calculated for different species and a significant difference (p < 0.01, α = 0.05) was observed when the 2015 dry season (ENSO period) was compared with the 2017 dry season ("control scenario"). In some days during the 2015 ENSO event, Tleaf approached 40°C for all studied species and the differences between Tleaf and Tair reached as high at 8°C (average difference: 1.65 ± 1.07°C). Generally, Tleaf was higher than Tair during the middle morning to early afternoon, and lower than Tair during the early morning, late afternoon and night. Our results support the hypothesis that partial stomatal closure allows for a recovery in ΨL during the afternoon period giving an observed counterclockwise hysteresis pattern between ΨL and Tleaf.

Beyond Static Benchmarking: Using Experimental Manipulations to Evaluate Land Model Assumptions.

Land models are often used to simulate terrestrial responses to future environmental changes, but these models are not commonly evaluated with data from experimental manipulations. Results from experimental manipulations can identify and evaluate model assumptions that are consistent with appropriate ecosystem responses to future environmental change. We conducted simulations using three coupled carbon-nitrogen versions of the Community Land Model (CLM, versions 4, 4.5, and-the newly developed-5), and compared the simulated response to nitrogen (N) and atmospheric carbon dioxide (CO2) enrichment with meta-analyses of observations from similar experimental manipulations. In control simulations, successive versions of CLM showed a poleward increase in gross primary productivity and an overall bias reduction, compared to FLUXNET-MTE observations. Simulations with N and CO2 enrichment demonstrate that CLM transitioned from a model that exhibited strong nitrogen limitation of the terrestrial carbon cycle (CLM4) to a model that showed greater responsiveness to elevated concentrations of CO2 in the atmosphere (CLM5). Overall, CLM5 simulations showed better agreement with observed ecosystem responses to experimental N and CO2 enrichment than previous versions of the model. These simulations also exposed shortcomings in structural assumptions and parameterizations. Specifically, no version of CLM captures changes in plant physiology, allocation, and nutrient uptake that are likely important aspects of terrestrial ecosystems' responses to environmental change. These highlight priority areas that should be addressed in future model developments. Moving forward, incorporating results from experimental manipulations into model benchmarking tools that are used to evaluate model performance will help increase confidence in terrestrial carbon cycle projections.

Variation in hydroclimate sustains tropical forest biomass and promotes functional diversity.

The fate of tropical forests under climate change is unclear as a result, in part, of the uncertainty in projected changes in precipitation and in the ability of vegetation models to capture the effects of drought-induced mortality on aboveground biomass (AGB). We evaluated the ability of a terrestrial biosphere model with demography and hydrodynamics (Ecosystem Demography, ED2-hydro) to simulate AGB and mortality of four tropical tree plant functional types (PFTs) that operate along light- and water-use axes. Model predictions were compared with observations of canopy trees at Barro Colorado Island (BCI), Panama. We then assessed the implications of eight hypothetical precipitation scenarios, including increased annual precipitation, reduced inter-annual variation, El Niño-related droughts and drier wet or dry seasons, on AGB and functional diversity of the model forest. When forced with observed meteorology, ED2-hydro predictions capture multiple BCI benchmarks. ED2-hydro predicts that AGB will be sustained under lower rainfall via shifts in the functional composition of the forest, except under the drier dry-season scenario. These results support the hypothesis that inter-annual variation in mean and seasonal precipitation promotes the coexistence of functionally diverse PFTs because of the relative differences in mortality rates. If the hydroclimate becomes chronically drier or wetter, functional evenness related to drought tolerance may decline.

Dependence of the evolution of carbon dynamics in the northern permafrost region on the trajectory of climate change.

We conducted a model-based assessment of changes in permafrost area and carbon storage for simulations driven by RCP4.5 and RCP8.5 projections between 2010 and 2299 for the northern permafrost region. All models simulating carbon represented soil with depth, a critical structural feature needed to represent the permafrost carbon-climate feedback, but that is not a universal feature of all climate models. Between 2010 and 2299, simulations indicated losses of permafrost between 3 and 5 million km2 for the RCP4.5 climate and between 6 and 16 million km2 for the RCP8.5 climate. For the RCP4.5 projection, cumulative change in soil carbon varied between 66-Pg C (1015-g carbon) loss to 70-Pg C gain. For the RCP8.5 projection, losses in soil carbon varied between 74 and 652 Pg C (mean loss, 341 Pg C). For the RCP4.5 projection, gains in vegetation carbon were largely responsible for the overall projected net gains in ecosystem carbon by 2299 (8- to 244-Pg C gains). In contrast, for the RCP8.5 projection, gains in vegetation carbon were not great enough to compensate for the losses of carbon projected by four of the five models; changes in ecosystem carbon ranged from a 641-Pg C loss to a 167-Pg C gain (mean, 208-Pg C loss). The models indicate that substantial net losses of ecosystem carbon would not occur until after 2100. This assessment suggests that effective mitigation efforts during the remainder of this century could attenuate the negative consequences of the permafrost carbon-climate feedback.

Land use change and El Niño-Southern Oscillation drive decadal carbon balance shifts in Southeast Asia.

An integrated understanding of the biogeochemical consequences of climate extremes and land use changes is needed to constrain land-surface feedbacks to atmospheric CO2 from associated climate change. Past assessments of the global carbon balance have shown particularly high uncertainty in Southeast Asia. Here, we use a combination of model ensembles to show that intensified land use change made Southeast Asia a strong source of CO2 from the 1980s to 1990s, whereas the region was close to carbon neutral in the 2000s due to an enhanced CO2 fertilization effect and absence of moderate-to-strong El Niño events. Our findings suggest that despite ongoing deforestation, CO2 emissions were substantially decreased during the 2000s, largely owing to milder climate that restores photosynthetic capacity and suppresses peat and deforestation fire emissions. The occurrence of strong El Niño events after 2009 suggests that the region has returned to conditions of increased vulnerability of carbon stocks.

Vegetation demographics in Earth System Models: A review of progress and priorities.

Numerous current efforts seek to improve the representation of ecosystem ecology and vegetation demographic processes within Earth System Models (ESMs). These developments are widely viewed as an important step in developing greater realism in predictions of future ecosystem states and fluxes. Increased realism, however, leads to increased model complexity, with new features raising a suite of ecological questions that require empirical constraints. Here, we review the developments that permit the representation of plant demographics in ESMs, and identify issues raised by these developments that highlight important gaps in ecological understanding. These issues inevitably translate into uncertainty in model projections but also allow models to be applied to new processes and questions concerning the dynamics of real-world ecosystems. We argue that stronger and more innovative connections to data, across the range of scales considered, are required to address these gaps in understanding. The development of first-generation land surface models as a unifying framework for ecophysiological understanding stimulated much research into plant physiological traits and gas exchange. Constraining predictions at ecologically relevant spatial and temporal scales will require a similar investment of effort and intensified inter-disciplinary communication.

Matrix approach to land carbon cycle modeling: A case study with the Community Land Model.

The terrestrial carbon (C) cycle has been commonly represented by a series of C balance equations to track C influxes into and effluxes out of individual pools in earth system models (ESMs). This representation matches our understanding of C cycle processes well but makes it difficult to track model behaviors. It is also computationally expensive, limiting the ability to conduct comprehensive parametric sensitivity analyses. To overcome these challenges, we have developed a matrix approach, which reorganizes the C balance equations in the original ESM into one matrix equation without changing any modeled C cycle processes and mechanisms. We applied the matrix approach to the Community Land Model (CLM4.5) with vertically-resolved biogeochemistry. The matrix equation exactly reproduces litter and soil organic carbon (SOC) dynamics of the standard CLM4.5 across different spatial-temporal scales. The matrix approach enables effective diagnosis of system properties such as C residence time and attribution of global change impacts to relevant processes. We illustrated, for example, the impacts of CO2 fertilization on litter and SOC dynamics can be easily decomposed into the relative contributions from C input, allocation of external C into different C pools, nitrogen regulation, altered soil environmental conditions, and vertical mixing along the soil profile. In addition, the matrix tool can accelerate model spin-up, permit thorough parametric sensitivity tests, enable pool-based data assimilation, and facilitate tracking and benchmarking of model behaviors. Overall, the matrix approach can make a broad range of future modeling activities more efficient and effective.

Carbon cycle confidence and uncertainty: Exploring variation among soil biogeochemical models.

Emerging insights into factors responsible for soil organic matter stabilization and decomposition are being applied in a variety of contexts, but new tools are needed to facilitate the understanding, evaluation, and improvement of soil biogeochemical theory and models at regional to global scales. To isolate the effects of model structural uncertainty on the global distribution of soil carbon stocks and turnover times we developed a soil biogeochemical testbed that forces three different soil models with consistent climate and plant productivity inputs. The models tested here include a first-order, microbial implicit approach (CASA-CNP), and two recently developed microbially explicit models that can be run at global scales (MIMICS and CORPSE). When forced with common environmental drivers, the soil models generated similar estimates of initial soil carbon stocks (roughly 1,400 Pg C globally, 0-100 cm), but each model shows a different functional relationship between mean annual temperature and inferred turnover times. Subsequently, the models made divergent projections about the fate of these soil carbon stocks over the 20th century, with models either gaining or losing over 20 Pg C globally between 1901 and 2010. Single-forcing experiments with changed inputs, temperature, and moisture suggest that uncertainty associated with freeze-thaw processes as well as soil textural effects on soil carbon stabilization were larger than direct temperature uncertainties among models. Finally, the models generated distinct projections about the timing and magnitude of seasonal heterotrophic respiration rates, again reflecting structural uncertainties that were related to environmental sensitivities and assumptions about physicochemical stabilization of soil organic matter. By providing a computationally tractable and numerically consistent framework to evaluate models we aim to better understand uncertainties among models and generate insights about factors regulating the turnover of soil organic matter.

A global trait-based approach to estimate leaf nitrogen functional allocation from observations.

Nitrogen is one of the most important nutrients for plant growth and a major constituent of proteins that regulate photosynthetic and respiratory processes. However, a comprehensive global analysis of nitrogen allocation in leaves for major processes with respect to different plant functional types (PFTs) is currently lacking. This study integrated observations from global databases with photosynthesis and respiration models to determine plant-functional-type-specific allocation patterns of leaf nitrogen for photosynthesis (Rubisco, electron transport, light absorption) and respiration (growth and maintenance), and by difference from observed total leaf nitrogen, an unexplained "residual" nitrogen pool. Based on our analysis, crops partition the largest fraction of nitrogen to photosynthesis (57%) and respiration (5%) followed by herbaceous plants (44% and 4%). Tropical broadleaf evergreen trees partition the least to photosynthesis (25%) and respiration (2%) followed by needle-leaved evergreen trees (28% and 3%). In trees (especially needle-leaved evergreen and tropical broadleaf evergreen trees) a large fraction (70% and 73%, respectively) of nitrogen was not explained by photosynthetic or respiratory functions. Compared to crops and herbaceous plants, this large residual pool is hypothesized to emerge from larger investments in cell wall proteins, lipids, amino acids, nucleic acid, CO2 fixation proteins (other than Rubisco), secondary compounds, and other proteins. Our estimates are different from previous studies due to differences in methodology and assumptions used in deriving nitrogen allocation estimates. Unlike previous studies, we integrate and infer nitrogen allocation estimates across multiple PFTs, and report substantial differences in nitrogen allocation across different PFTs. The resulting pattern of nitrogen allocation provides insights on mechanisms that operate at a cellular scale within leaves, and can be integrated with ecosystem models to derive emergent properties of ecosystem productivity at local, regional, and global scales.