Historical tree phenology data reveal the seasonal rhythms of the Congo Basin rainforest.

Tropical forest phenology directly affects regional carbon cycles, but the relation between species-specific and whole-canopy phenology remains largely uncharacterized. We present a unique analysis of historical tropical tree phenology collected in the central Congo Basin, before large-scale impacts of human-induced climate change. Ground-based long-term (1937-1956) phenological observations of 140 tropical tree species are recovered, species-specific phenological patterns analyzed and related to historical meteorological records, and scaled to characterize stand-level canopy dynamics. High phenological variability within and across species and in climate-phenology relationships is observed. The onset of leaf phenophases in deciduous species was triggered by drought and light availability for a subset of species and showed a species-specific decoupling in time along a bi-modal seasonality. The majority of the species remain evergreen, although central African forests experience relatively low rainfall. Annually a maximum of 1.5% of the canopy is in leaf senescence or leaf turnover, with overall phenological variability dominated by a few deciduous species, while substantial variability is attributed to asynchronous events of large and/or abundant trees. Our results underscore the importance of accounting for constituent signals in canopy-wide scaling and the interpretation of remotely sensed phenology signals.

Trade-offs in biodiversity and ecosystem services between edges and interiors in European forests.

Forest biodiversity and ecosystem services are hitherto predominantly quantified in forest interiors, well away from edges. However, these edges also represent a substantial proportion of the global forest cover. Here we quantified plant biodiversity and ecosystem service indicators in 225 plots along forest edge-to-interior transects across Europe. We found strong trade-offs: phylogenetic diversity (evolutionary measure of biodiversity), proportion of forest specialists, decomposition and heatwave buffering increased towards the interior, whereas species richness, nectar production potential, stemwood biomass and tree regeneration decreased. These trade-offs were mainly driven by edge-to-interior structural differences. As fragmentation continues, recognizing the role of forest edges is crucial for integrating biodiversity and ecosystem service considerations into sustainable forest management and policy.

Damage to living trees contributes to almost half of the biomass losses in tropical forests.

Accurate estimates of forest biomass stocks and fluxes are needed to quantify global carbon budgets and assess the response of forests to climate change. However, most forest inventories consider tree mortality as the only aboveground biomass (AGB) loss without accounting for losses via damage to living trees: branchfall, trunk breakage, and wood decay. Here, we use ~151,000 annual records of tree survival and structural completeness to compare AGB loss via damage to living trees to total AGB loss (mortality + damage) in seven tropical forests widely distributed across environmental conditions. We find that 42% (3.62 Mg ha-1 year-1 ; 95% confidence interval [CI] 2.36-5.25) of total AGB loss (8.72 Mg ha-1 year-1 ; CI 5.57-12.86) is due to damage to living trees. Total AGB loss was highly variable among forests, but these differences were mainly caused by site variability in damage-related AGB losses rather than by mortality-related AGB losses. We show that conventional forest inventories overestimate stand-level AGB stocks by 4% (1%-17% range across forests) because assume structurally complete trees, underestimate total AGB loss by 29% (6%-57% range across forests) due to overlooked damage-related AGB losses, and overestimate AGB loss via mortality by 22% (7%-80% range across forests) because of the assumption that trees are undamaged before dying. Our results indicate that forest carbon fluxes are higher than previously thought. Damage on living trees is an underappreciated component of the forest carbon cycle that is likely to become even more important as the frequency and severity of forest disturbances increase.

Atmospheric phosphorus deposition amplifies carbon sinks in simulations of a tropical forest in Central Africa.

Spatial redistribution of nutrients by atmospheric transport and deposition could theoretically act as a continental-scale mechanism which counteracts declines in soil fertility caused by nutrient lock-up in accumulating biomass in tropical forests in Central Africa. However, to what extent it affects carbon sinks in forests remains elusive. Here we use a terrestrial biosphere model to quantify the impact of changes in atmospheric nitrogen and phosphorus deposition on plant nutrition and biomass carbon sink at a typical lowland forest site in Central Africa. We find that the increase in nutrient deposition since the 1980s could have contributed to the carbon sink over the past four decades up to an extent which is similar to that from the combined effects of increasing atmospheric carbon dioxide and climate change. Furthermore, we find that the modelled carbon sink responds to changes in phosphorus deposition, but less so to nitrogen deposition. The pronounced response of ecosystem productivity to changes in nutrient deposition illustrates a potential mechanism that could control carbon sinks in Central Africa. Monitoring the quantity and quality of nutrient deposition is needed in this region, given the changes in nutrient deposition due to human land use.

Accurate phenology analyses require bud traits and energy budgets.

Spring phenology is mainly driven by temperature in extratropical ecosystems. Recent evidence highlighted the key role of micrometeorology and bud temperature on delaying or advancing leaf unfolding. Yet, phenology studies, either using ground-based or remote sensing observations, always substitute plant tissue temperature by air temperature. In fact, temperatures differ substantially between plant tissues and the air because plants absorb and lose energy. Here, we build on recent observations and well-established energy balance theories to discuss how solar radiation, wind and bud traits might affect our interpretation of spring phenology sensitivity to warming. We show that air temperature might be an imprecise and biased predictor of bud temperature. Better characterizing the plants' phenological response to warming will require new observations of bud traits and temperature for accurately quantifying their energy budget. As consistent micrometeorology datasets are still scarce, new approaches coupling energy budget modelling and plant traits could help to improve phenology analyses across scales.

FunAndes - A functional trait database of Andean plants.

We introduce the FunAndes database, a compilation of functional trait data for the Andean flora spanning six countries. FunAndes contains data on 24 traits across 2,694 taxa, for a total of 105,466 entries. The database features plant-morphological attributes including growth form, and leaf, stem, and wood traits measured at the species or individual level, together with geographic metadata (i.e., coordinates and elevation). FunAndes follows the field names, trait descriptions and units of measurement of the TRY database. It is currently available in open access in the FIGSHARE data repository, and will be part of TRY's next release. Open access trait data from Andean plants will contribute to ecological research in the region, the most species rich terrestrial biodiversity hotspot.

Liana optical traits increase tropical forest albedo and reduce ecosystem productivity.

Lianas are a key growth form in tropical forests. Their lack of self-supporting tissues and their vertical position on top of the canopy make them strong competitors of resources. A few pioneer studies have shown that liana optical traits differ on average from those of colocated trees. Those trait discrepancies were hypothesized to be responsible for the competitive advantage of lianas over trees. Yet, in the absence of reliable modelling tools, it is impossible to unravel their impact on the forest energy balance, light competition, and on the liana success in Neotropical forests. To bridge this gap, we performed a meta-analysis of the literature to gather all published liana leaf optical spectra, as well as all canopy spectra measured over different levels of liana infestation. We then used a Bayesian data assimilation framework applied to two radiative transfer models (RTMs) covering the leaf and canopy scales to derive tropical tree and liana trait distributions, which finally informed a full dynamic vegetation model. According to the RTMs inversion, lianas grew thinner, more horizontal leaves with lower pigment concentrations. Those traits made the lianas very efficient at light interception and significantly modified the forest energy balance and its carbon cycle. While forest albedo increased by 14% in the shortwave, light availability was reduced in the understorey (-30% of the PAR radiation) and soil temperature decreased by 0.5°C. Those liana-specific traits were also responsible for a significant reduction of tree (-19%) and ecosystem (-7%) gross primary productivity (GPP) while lianas benefited from them (their GPP increased by +27%). This study provides a novel mechanistic explanation to the increase in liana abundance, new evidence of the impact of lianas on forest functioning, and paves the way for the evaluation of the large-scale impacts of lianas on forest biogeochemical cycles.

Forest above-ground volume assessments with terrestrial laser scanning: a ground-truth validation experiment in temperate, managed forests.

BACKGROUND AND AIMS: Quantifying the Earth's forest above-ground biomass (AGB) is indispensable for effective climate action and developing forest policy. Yet, current allometric scaling models (ASMs) to estimate AGB suffer several drawbacks related to model selection and uncertainties about calibration data traceability. Terrestrial laser scanning (TLS) offers a promising non-destructive alternative. Tree volume is reconstructed from TLS point clouds with quantitative structure models (QSMs) and converted to AGB with wood basic density. Earlier studies have found overall TLS-derived forest volume estimates to be accurate, but highlighted problems for reconstructing finer branches. Our objective was to evaluate TLS for estimating tree volumes by comparison with reference volumes and volumes from ASMs. METHODS: We quantified the woody volume of 65 trees in Belgium (from 77 to 2800 L; Pinus sylvestris, Fagus sylvatica, Larix decidua, and Fraxinus excelsior) with QSMs and destructive reference measurements. We tested a volume expansion factor (VEF) approach by multiplying the solid and merchantable volume from QSMs by literature VEF values. KEY RESULTS: Stem volume was reliably estimated with TLS. Total volume was overestimated by +21 % using original QSMs, by +9 % and -12 % using two sets of VEF-augmented QSMs, and by -7.3 % using best-available ASMs. The most accurate method differed per site, and the prediction errors for each method varied considerably between sites. CONCLUSIONS: VEF-augmented QSMs were only slightly better than original QSMs for estimating tree volume for common species in temperate forests. Despite satisfying estimates with ASMs, the model choice was a large source of uncertainty, and species-specific models did not always exist. Therefore, we advocate for further improving tree volume reconstructions with QSMs, especially for fine branches, instead of collecting more ground-truth data to calibrate VEF and allometric models. Promising developments such as improved co-registration and smarter filtering approaches are ongoing to further constrain volumetric errors in TLS-derived estimates.

Fire-derived phosphorus fertilization of African tropical forests.

Central African tropical forests face increasing anthropogenic pressures, particularly in the form of deforestation and land-use conversion to agriculture. The long-term effects of this transformation of pristine forests to fallow-based agroecosystems and secondary forests on biogeochemical cycles that drive forest functioning are poorly understood. Here, we show that biomass burning on the African continent results in high phosphorus (P) deposition on an equatorial forest via fire-derived atmospheric emissions. Furthermore, we show that deposition loads increase with forest regrowth age, likely due to increasing canopy complexity, ranging from 0.4 kg P ha-1 yr-1 on agricultural fields to 3.1 kg P ha-1 yr-1 on old secondary forests. In forest systems, canopy wash-off of dry P deposition increases with rainfall amount, highlighting how tropical forest canopies act as dynamic reservoirs for enhanced addition of this essential plant nutrient. Overall, the observed P deposition load at the study site is substantial and demonstrates the importance of canopy trapping as a pathway for nutrient input into forest ecosystems.

Understanding 3D structural complexity of individual Scots pine trees with different management history.

Tree functional traits together with processes such as forest regeneration, growth, and mortality affect forest and tree structure. Forest management inherently impacts these processes. Moreover, forest structure, biodiversity, resilience, and carbon uptake can be sustained and enhanced with forest management activities. To assess structural complexity of individual trees, comprehensive and quantitative measures are needed, and they are often lacking for current forest management practices. Here, we utilized 3D information from individual Scots pine (Pinus sylvestris L.) trees obtained with terrestrial laser scanning to, first, assess effects of forest management on structural complexity of individual trees and, second, understand relationship between several tree attributes and structural complexity. We studied structural complexity of individual trees represented by a single scale-independent metric called "box dimension." This study aimed at identifying drivers affecting structural complexity of individual Scots pine trees in boreal forest conditions. The results showed that thinning increased structural complexity of individual Scots pine trees. Furthermore, we found a relationship between structural complexity and stem and crown size and shape as well as tree growth. Thus, it can be concluded that forest management affected structural complexity of individual Scots pine trees in managed boreal forests, and stem, crown, and growth attributes were identified as drivers of it.

Unraveling the relative role of light and water competition between lianas and trees in tropical forests: A vegetation model analysis.

Despite their low contribution to forest carbon stocks, lianas (woody vines) play an important role in the carbon dynamics of tropical forests. As structural parasites, they hinder tree survival, growth and fecundity; hence, they negatively impact net ecosystem productivity and long-term carbon sequestration.Competition (for water and light) drives various forest processes and depends on the local abundance of resources over time. However, evaluating the relative role of resource availability on the interactions between lianas and trees from empirical observations is particularly challenging. Previous approaches have used labour-intensive and ecosystem-scale manipulation experiments, which are infeasible in most situations.We propose to circumvent this challenge by evaluating the uncertainty of water and light capture processes of a process-based vegetation model (ED2) including the liana growth form. We further developed the liana plant functional type in ED2 to mechanistically simulate water uptake and transport from roots to leaves, and start the model from prescribed initial conditions. We then used the PEcAn bioinformatics platform to constrain liana parameters and run uncertainty analyses.Baseline runs successfully reproduced ecosystem gas exchange fluxes (gross primary productivity and latent heat) and forest structural features (leaf area index, aboveground biomass) in two sites (Barro Colorado Island, Panama and Paracou, French Guiana) characterized by different rainfall regimes and levels of liana abundance.Model uncertainty analyses revealed that water limitation was the factor driving the competition between trees and lianas at the drier site (BCI), and during the relatively short dry season of the wetter site (Paracou). In young patches, light competition dominated in Paracou but alternated with water competition between the wet and the dry season on BCI according to the model simulations.The modelling workflow also identified key liana traits (photosynthetic quantum efficiency, stomatal regulation parameters, allometric relationships) and processes (water use, respiration, climbing) driving the model uncertainty. They should be considered as priorities for future data acquisition and model development to improve predictions of the carbon dynamics of liana-infested forests. Synthesis. Competition for water plays a larger role in the interaction between lianas and trees than previously hypothesized, as demonstrated by simulations from a process-based vegetation model.

High photosynthetic capacity of Sahelian C3 and C4 plants.

The semi-arid ecosystems of the African Sahel play an important role in the global carbon cycle and are among the most sensitive ecosystems to future environmental pressures. Still, basic data of photosynthetic characteristics of Sahelian vegetation are very limited, preventing us to properly understand these ecosystems and to project their response to future global changes. Here, we aim to study and quantify key leaf traits, including photosynthetic parameters and leaf nutrients (Nleaf and Pleaf), of common C3 and C4 Sahelian plants (trees, lianas, and grasses) at the Dahra field site (Senegal). Dahra is a reference site for grazed semi-arid Sahelian savannah ecosystems in carbon cycle studies. Within the studied species, we found that photosynthetic parameters varied considerably between functional types. We also found significant relationships between and within photosynthetic parameters and leaf traits which mostly differed in their slopes from C3 to C4 plants. In agreement with the leaf economic spectrum, strong relationships (R2 = 0.71) were found between SLA and Nleaf whereby C3 and C4 plants showed very similar relationships. By comparing our data to a global dataset of plant traits, we show that measured Sahelian plants exhibit higher photosynthetic capacity (Asat) compared to the non-Sahelian vegetation, with values that are on average a fourfold of the global average. Moreover, Sahelian C3 plants showed photosynthetic nutrient use efficiencies that were on average roughly twice as high as global averages. We interpreted these results as the potential adaptation of Sahelian plants to short growing season lengths via an efficient nutrient allocation to optimize photosynthesis during this period. Our study provides robust estimates of key functional traits, but also traits relationships that will help to calibrate and validate vegetation models over this data-poor region.

Comparable canopy and soil free-living nitrogen fixation rates in a lowland tropical forest.

Biological nitrogen fixation (BNF) is a fundamental part of nitrogen cycling in tropical forests, yet little is known about the contribution made by free-living nitrogen fixers inhabiting the often-extensive forest canopy. We used the acetylene reduction assay, calibrated with 15N2, to measure free-living BNF on forest canopy leaves, vascular epiphytes, bryophytes and canopy soil, as well as on the forest floor in leaf litter and soil. We used a combination of calculated and published component densities to upscale free-living BNF rates to the forest level. We found that bryophytes and leaves situated in the canopy in particular displayed high mass-based rates of free-living BNF. Additionally, we calculated that nearly 2 kg of nitrogen enters the forest ecosystem through free-living BNF every year, 40% of which was fixed by the various canopy components. Our results reveal that in the studied tropical lowland forest a large part of the nitrogen input through free-living BNF stems from the canopy, but also that the total nitrogen inputs by free-living BNF are lower than previously thought and comparable to the inputs of reactive nitrogen by atmospheric deposition.

Drivers of carbon stocks in forest edges across Europe.

Forests play a key role in global carbon cycling and sequestration. However, the potential for carbon drawdown is affected by forest fragmentation and resulting changes in microclimate, nutrient inputs, disturbance and productivity near edges. Up to 20% of the global forested area lies within 100 m of an edge and, even in temperate forests, knowledge on how edge conditions affect carbon stocks and how far this influence penetrates into forest interiors is scarce. Here we studied carbon stocks in the aboveground biomass, forest floor and the mineral topsoil in 225 plots in deciduous forest edges across Europe and tested the impact of macroclimate, nitrogen deposition and smaller-grained drivers (e.g. microclimate) on these stocks. Total carbon and carbon in the aboveground biomass stock were on average 39% and 95% higher at the forest edge than 100 m into the interior. The increase in the aboveground biomass stock close to the edge was mainly related to enhanced nitrogen deposition. No edge influence was found for stocks in the mineral topsoil. Edge-to-interior gradients in forest floor carbon changed across latitude: carbon stocks in the forest floor were higher near the edge in southern Europe. Forest floor carbon decreased with increasing litter quality (i.e. high decomposition rate) and decreasing plant area index, whereas higher soil temperatures negatively affected the mineral topsoil carbon. Based on high-resolution forest fragmentation maps, we estimate that the additional carbon stored in deciduous forest edges across Europe amounts to not less than 183 Tg carbon, which is equivalent to the storage capacity of 1 million ha of additional forest. This study underpins the importance of including edge influences when quantifying the carbon stocks in temperate forests and stresses the importance of preserving natural forest edges and small forest patches with a high edge-to-interior surface area.

Lianas in silico, ecological insights from a model of structural parasitism.

Tropical forests are a critical component of the Earth system, storing half of the global forest carbon stocks and accounting for a third of terrestrial photosynthesis. Lianas are structural parasites that can substantially reduce the carbon sequestration capacity of these forests. Simulations of this peculiar growth form have only recently started and a single vegetation model included lianas so far. In this work we present a new liana implementation within the individual based model Formind. Initial tests indicate high structural realism both horizontal and vertical. In particular, we benchmarked the model against empirical observations of size distribution, mean liana cluster size and vertical leaf distribution for the Paracou site in French Guiana. Our model predicted a reduction of above-ground biomass between 10% for mature stands to 45% for secondary plots upon inclusion of lianas in the simulations. The reduced biomass was the result of a lower productivity due to a combination of lower tree photosynthesis and high liana respiration. We evaluated structural metrics (LAI, basal area, mean tree-height) and carbon fluxes (GPP, respiration) by comparing simulations with and without lianas. At the equilibrium, liana productivity was 1.9t C ha - 1 y - 1 , or 23% of the total GPP and the forest carbon stocks were between 5% and 11% lower in simulations with lianas. We also highlight the main strengths and limitations of this new approach and propose new field measurements to further the understanding of liana ecology in a modelling framework.

Century-long apparent decrease in intrinsic water-use efficiency with no evidence of progressive nutrient limitation in African tropical forests.

Forests exhibit leaf- and ecosystem-level responses to environmental changes. Specifically, rising carbon dioxide (CO2 ) levels over the past century are expected to have increased the intrinsic water-use efficiency (iWUE) of tropical trees while the ecosystem is gradually pushed into progressive nutrient limitation. Due to the long-term character of these changes, however, observational datasets to validate both paradigms are limited in space and time. In this study, we used a unique herbarium record to go back nearly a century and show that despite the rise in CO2 concentrations, iWUE has decreased in central African tropical trees in the Congo Basin. Although we find evidence that points to leaf-level adaptation to increasing CO2 -that is, increasing photosynthesis-related nutrients and decreasing maximum stomatal conductance, a decrease in leaf δ13 C clearly indicates a decreasing iWUE over time. Additionally, the stoichiometric carbon to nitrogen and nitrogen to phosphorus ratios in the leaves show no sign of progressive nutrient limitation as they have remained constant since 1938, which suggests that nutrients have not increasingly limited productivity in this biome. Altogether, the data suggest that other environmental factors, such as increasing temperature, might have negatively affected net photosynthesis and consequently downregulated the iWUE. Results from this study reveal that the second largest tropical forest on Earth has responded differently to recent environmental changes than expected, highlighting the need for further on-ground monitoring in the Congo Basin.

Semi-automatic extraction of liana stems from terrestrial LiDAR point clouds of tropical rainforests.

Lianas are key structural elements of tropical forests having a large impact on the global carbon cycle by reducing tree growth and increasing tree mortality. Despite the reported increasing abundance of lianas across neotropics, very few studies have attempted to quantify the impact of lianas on tree and forest structure. Recent advances in high resolution terrestrial laser scanning (TLS) systems have enabled us to quantify the forest structure, in an unprecedented detail. However, the uptake of TLS technology to study lianas has not kept up with the same pace as it has for trees. The slower technological adoption of TLS to study lianas is due to the lack of methods to study these complex growth forms. In this study, we present a semi-automatic method to extract liana woody components from plot-level TLS data of a tropical rainforest. We tested the method in eight plots from two different tropical rainforest sites (two in Gigante Peninsula, Panama and six in Nouragues, French Guiana) along an increasing gradient of liana infestation (from plots with low liana density to plots with very high liana density). Our method uses a machine learning model based on the Random Forest (RF) algorithm. The RF algorithm is trained on the eigen features extracted from the points in 3D at multiple spatial scales. The RF based liana stem extraction method successfully extracts on average 58% of liana woody points in our dataset with a high precision of 88%. We also present simple post-processing steps that increase the percentage of extracted liana stems from 54% to 90% in Nouragues and 65% to 70% in Gigante Peninsula without compromising on the precision. We provide the entire processing pipeline as an open source python package. Our method will facilitate new research to study lianas as it enables the monitoring of liana abundance, growth and biomass in forest plots. In addition, the method facilitates the easier processing of 3D data to study tree structure from a liana-infested forest.

Modeling the impact of liana infestation on the demography and carbon cycle of tropical forests.

There is mounting empirical evidence that lianas affect the carbon cycle of tropical forests. However, no single vegetation model takes into account this growth form, although such efforts could greatly improve the predictions of carbon dynamics in tropical forests. In this study, we incorporated a novel mechanistic representation of lianas in a dynamic global vegetation model (the Ecosystem Demography Model). We developed a liana-specific plant functional type and mechanisms representing liana-tree interactions (such as light competition, liana-specific allometries, and attachment to host trees) and parameterized them according to a comprehensive literature meta-analysis. We tested the model for an old-growth forest (Paracou, French Guiana) and a secondary forest (Gigante Peninsula, Panama). The resulting model simulations captured many features of the two forests characterized by different levels of liana infestation as revealed by a systematic comparison of the model outputs with empirical data, including local census data from forest inventories, eddy flux tower data, and terrestrial laser scanner-derived forest vertical structure. The inclusion of lianas in the simulations reduced the secondary forest net productivity by up to 0.46 tC  ha-1  year-1 , which corresponds to a limited relative reduction of 2.6% in comparison with a reference simulation without lianas. However, this resulted in significantly reduced accumulated above-ground biomass after 70 years of regrowth by up to 20 tC /ha (19% of the reference simulation). Ultimately, the simulated negative impact of lianas on the total biomass was almost completely cancelled out when the forest reached an old-growth successional stage. Our findings suggest that lianas negatively influence the forest potential carbon sink strength, especially for young, disturbed, liana-rich sites. In light of the critical role that lianas play in the profound changes currently experienced by tropical forests, this new model provides a robust numerical tool to forecast the impact of lianas on tropical forest carbon sinks.

Large-sized rare tree species contribute disproportionately to functional diversity in resource acquisition in African tropical forest.

Increasing evidence is available for a positive effect of biodiversity on ecosystem productivity and standing biomass, also in highly diverse systems as tropical forests. Biodiversity conservation could therefore be a critical aspect of climate mitigation policies. There is, however, limited understanding of the role of individual species for this relationship, which could aid in focusing conservation efforts and forest management planning. This study characterizes the functional specialization and redundancy for 95% of all tree species (basal area weighted percentage) in a diverse tropical forest in the central Congo Basin and relates this to species' abundance, contribution to aboveground carbon, and maximum size. Functional characterization is based on a set of traits related to resource acquisition (wood density, specific leaf area, leaf carbon, nitrogen and phosphorus content, and leaf stable carbon isotope composition). We show that within both mixed and monodominant tropical forest ecosystems, the highest functional specialization and lowest functional redundancy are solely found in rare tree species and significantly more in rare species holding large-sized individuals. Rare species cover the entire range of low and high functional redundancy, contributing both unique and redundant functions. Loss of species supporting functional redundancy could be buffered by other species in the community, including more abundant species. This is not the case for species supporting high functional specialization and low functional redundancy, which would need specific conservation attention. In terms of tropical forest management planning, we argue that specific conservation of large-sized trees is imperative for long-term maintenance of ecosystem functioning.

Air temperature optima of vegetation productivity across global biomes.

The global distribution of the optimum air temperature for ecosystem-level gross primary productivity ([Formula: see text]) is poorly understood, despite its importance for ecosystem carbon uptake under future warming. We provide empirical evidence for the existence of such an optimum, using measurements of in situ eddy covariance and satellite-derived proxies, and report its global distribution. [Formula: see text] is consistently lower than the physiological optimum temperature of leaf-level photosynthetic capacity, which typically exceeds 30 °C. The global average [Formula: see text] is estimated to be 23 ± 6 °C, with warmer regions having higher [Formula: see text] values than colder regions. In tropical forests in particular, [Formula: see text] is close to growing-season air temperature and is projected to fall below it under all scenarios of future climate, suggesting a limited safe operating space for these ecosystems under future warming.