A constraint on historic growth in global photosynthesis due to rising CO2

Theory predicts that rising CO2 increases global photosynthesis, a process known as CO2 fertilization, and that this is responsible for much of the current terrestrial carbon sink. The estimated magnitude of the historic CO2 fertilization, however, differs by an order of magnitude between long-term proxies, remote sensing-based estimates and terrestrial biosphere models. Here we constrain the likely historic effect of CO2 on global photosynthesis by combining terrestrial biosphere models, ecological optimality theory, remote sensing approaches and an emergent constraint based on global carbon budget estimates. Our analysis suggests that CO2 fertilization increased global annual terrestrial photosynthesis by 13.5 ± 3.5% or 15.9 ± 2.9 PgC (mean ± s.d.) between 1981 and 2020. Our results help resolve conflicting estimates of the historic sensitivity of global terrestrial photosynthesis to CO2 and highlight the large impact anthropogenic emissions have had on ecosystems worldwide. The authors incorporate terrestrial biosphere models with ecological optimality theory, remote sensing and global carbon budget estimates to constrain the historical effects of CO2 on photosynthesis. They show that CO2 fertilization likely increased global photosynthesis by 13.5% between 1981 and 2020.

We reconciled the apparent difference between the TBM-inferred and satellite-based estimates of the sensitivity of GPP to CO 2 (Fig. 1) by using first-principles theory to incorporate the direct effect of increasing CO 2 on C 3 LUE in the satellite-based estimates.We refer to RS estimates that incorporate theory of the direct effect of CO 2 on LUE as the modified RS-based methods hereafter.The direct effect of CO 2 on LUE reflects the increasing competitiveness of CO 2 relative to O 2 for the active sites of the ribulose-1,5-bisphosphate carboxylase-oxygenase (RuBisCO) enzyme and the increasing competitiveness of CO 2 as atmospheric concentrations rise (Methods).To do so we considered two distinct classes or satellite-based estimates.The first is a commonly used LUE approach based on the Moderate Resolution Imaging Spectroradiometer (MODIS) algorithm (the MA approach) and the second is an ML approach that integrates both satellite and ground observations of ecosystem carbon fluxes.The direct effect of increasing CO 2 on the LUE of canopy photosynthesis 13 was roughly twice as large as the indirect effect of increasing canopy leaf area and thus increasing the fraction of absorbed photosynthetically active radiation(fAPAR), represented in the ML and MA approaches (Fig. 1a,b).The long-term sensitivity of the RS-based estimates of GPP modified to account for both the direct (β LUE R ) and indirect ( β fAPAR R ) effect of increasing CO 2 ( β GPP R ; equation (1)) was 0.50 ± 0.1 (mean ± s.d.) and 0.46 ± 0.1 for the ML and MA approaches, respectively (Fig. 1b), compared to 0.16 ± 0.05 and 0.16 ± 0.06 for the original ML and MA-based estimates, respectively (Fig. 1b).The long-term increase in GPP from the updated RS-based estimates thus more closely approximated that of the TBM ensemble mean ( β GPP R = 0.59 ± 0.16) (Fig. 1b).The modified RS-based methods predict a 7.27 ± 0.7% (ML) and 6.72 ± 0.9% (MA) increase in global annual GPP for a 14.5% increase in atmospheric CO 2 between 1982 and 2012.
Despite the agreement between the updated satellite methods and the TBM model ensemble (Fig. 1b), there is a large spread in individual TBM sensitivities and the true sensitivity is uncertain because of the lack of a comparable observational record.To address this issue, we proposed a constraint on the historic response of surface reflectance, from which vegetation solar energy absorption can be derived.As they integrate land surface observations, they are often regarded as a benchmark to which TBMs should be compared 11 .Such comparisons generally suggest that TBMs overestimate the change in global terrestrial photosynthesis due to too high a sensitivity of photosynthesis to increasing CO 2 (refs.6,11).However, satellite-TBM comparisons are mired by the fact that most satellite-based estimates, be they machine learning (ML) or algorithmically based, do not incorporate the universally observed direct effect of increasing CO 2 on the light-use efficiency (LUE) of leaves of C 3 vegetation 13 .This is because the direct effect of increasing CO 2 on LUE is not directly observable from space 14 .In contrast, observation-based proxies, based on ice-core records of COS 8,15 , eddy-covariance networks 12 and herbarium and field-based deuterium isotopomers 9 , suggest that TBMs may underestimate the sensitivity of global photosynthesis to CO 2 .TBMs themselves show a large range of sensitivities of global terrestrial photosynthesis to CO 2 (refs.10,16,17), though few demonstrate sensitivities as low as the average satellite-inferred values 6,14 or typically as high as those derived from the COS or deuterium proxies 8,9,17 .The spread in estimates of the sensitivity of global terrestrial photosynthesis to CO 2 and the lack of a global constraint, constitutes a large source of uncertainty in future projections of the Earth system 18 and hinders attribution of the various processes responsible for long-term changes in the global terrestrial carbon cycle.
Here, we use remote sensing (RS) observations informed with ecological optimality theory to help constrain the historic response of photosynthesis to rising CO 2. We develop a method to incorporate the direct effect of CO 2 on the rate of canopy-level terrestrial gross primary photosynthesis (GPP) in established satellite-based approaches.We do so using first-principles theory of photosynthetic carbon fixation 19,20 and generate 30-year global datasets of satellite-derived GPP.In addition, we identify an emergent multimodel relationship [21][22][23] between the modelled terrestrial carbon sink and the sensitivity of photosynthesis to CO 2 from the Trends in Net Land-Atmosphere Carbon Exchanges project (TRENDY 24 ).When combined, these approaches constrain the range of plausible estimates of the historic effect of CO 2 on global GPP, resolving the large apparent difference between Fig.  photosynthesis to rising CO 2 by combining the TRENDY modelled sensitivity of global GPP, for which no direct observations exist, with the magnitude of the cumulative global terrestrial residual carbon sink between 1982 and 2016 (S LAND ), for which there are constrained estimates 25,26 .Using cumulative S LAND estimates as opposed to shorter time periods reduces the influence of annual or decadal random error.Model sensitivities of photosynthesis to CO 2 (Supplementary Table 2) were positively correlated (r = 0.70, P < 0.01) with the magnitude of the modelled cumulative terrestrial sink on a multidecadal scale (Fig. 2a), with a stronger CO 2 fertilization effect leading to a larger modelled cumulative sink.This emergent relationship [21][22][23] , provides an opportunity to constrain the wide range in estimates of the sensitivity of GPP to CO 2 with the observed magnitude of S LAND , particularly when combined with the results of the modified RS-based estimates.The full distribution, which includes the TBMs (Supplementary Table 2) and the original RS-based estimates, provides an estimate of β GPP R of 0.54 ± 0.21 (mean ± s.d.; Fig. 2b), which is lower than that derived from the distribution of TBMs ( β GPP R = 0.59 ± 0.16).The posterior TBM distribution, formed by bootstrapping the cumulative land-sink emergent constraint relationship (Fig. 2a) provides a constrained estimate of β GPP R of 0.62 ± 0.14 (Fig. 2b).This is 30% lower than the maximum unconstrained estimate and over 200% higher than that of the original RS-based approaches.The constrained distribution represents a 33% reduction in uncertainty compared to the full distribution of β GPP R (Fig. 2b) and a 13.5% reduction compared to the unconstrained TBM ensemble (Fig. 2b).
Results from the updated RS estimates and the emergent constraint provide a point of comparison for other reported estimates of the sensitivity of global terrestrial photosynthesis to CO 2 .A long-term COS proxy has been proposed 8 , which simulates photosynthetic change on the basis of a mass balance of global COS sources and sinks from 1900 to 2013 and suggests an increase in photosynthesis equivalent to an effective β GPP R of 0.94 (Supplementary Table 1).This is comparable to the highest sensitivity of the TBM models used here 17 .The COS estimate, however, integrates over a longer time period and therefore potentially captures changes in the land surface unrelated to CO 2 , such as reforestation and the agricultural green revolution 27 , and is thus not directly comparable to the emergent constraint and updated RS estimates presented here.Another proxy, based on deuterium isotopomers gathered from herbarium specimens and field trials 9 , suggests an historic change equivalent to a β GPP R of 1.03 (Supplementary Table 1).Although higher than that derived from COS, the deuterium isotopomer estimate reflects the effect of increasing CO 2 on photosynthesis for leaves in full sunlight.As shaded leaves experience stronger light limitation, which results in a lower sensitivity to CO 2 , COS-based estimates could thus reasonably be expected to be higher than the canopy integrated sensitivity.Our results indicate that such larger implied sensitivities are probably overestimates (Fig. 2).
The closer agreement between the updated RS approaches and the TBMs (Fig. 1) allows for their response to CO 2 to be probed more deeply.The sensitivity of C 3 photosynthesis to CO 2 increases strongly with temperature 28 (Fig. 3a; equations ( 2)-( 7)) because the suppression of oxygenation by RuBisCO with increasing CO 2 is greater at higher temperatures.Reduced RuBisCO oxygenation reduces photorespiration at high temperatures, as represented by the temperature dependence of the photosynthetic CO 2 compensation point (Γ*, equation ( 3)).The resulting latitudinal gradient is reproduced by both the TBMs examined (Fig. 3b) and the updated RS approaches (Fig. 3c-e).The results indicate that the influence of CO 2 on photosynthesis at high latitudes is limited due to low temperatures.Estimates of the long-term change in GPP from the updated RS approaches show large changes, particularly in areas of intensive agriculture such as the midwestern United States, central and northern Europe and India (Fig. 3c,d).Compared to the RS approaches (Fig. 3d), the TBMs predict smaller increases in arid mid-and low-latitude regions, particularly in Australia and South Africa but much larger increases in the productive croplands and tropical and temperate forests (Fig. 3d).The lower TBM sensitivity, in particular of shrublands (Fig. 3f), is potentially due to poorly represented TBM processes such as the positive relationship between CO 2 and woody shrub expansion 29 .The lower TBM sensitivities could be from inaccurate representation of greening trends that arise from changes in land management practices such as reforestation 30 .The relatively higher TBM sensitivity regions, particularly tropical forests (Fig. 3), may be due to insufficient TBM representation of nutrient constraints 31 , or the saturation of RS vegetation indices at high leaf area 32 , reflecting large uncertainty about the response of tropical forest photosynthesis to CO 2 (ref.33).In general, the magnitude of the TBM and updated satellite β GPP R suggests that the global terrestrial photosynthetic response to CO 2 is consistent with the response of the light-limited photosynthetic rate which has also been suggested by observations of photosynthesis and biomass changes at the ecosystem scale [34][35][36] , theoretical models 37,38 and by model results showing that    electron-transport limited leaves are responsible for most global carbon assimilated through photosynthesis 39 .
As with any application of the emergent constraint technique, it is important to highlight that many factors could lead to biases and undermine the robustness of the derived constraint.Of primary concern is the potential for emergent constraints to rely on spurious cross-model correlations that are not based on a clear physical relationship 40 .The constraint we identify is based on the relationship between CO 2 and the land sink, for which there is ample observational and theoretical support 3,4 .Although CO 2 fertilization is by no means the sole likely reason for an increasing land sink 4 (other contributions arise from forest regrowth, nitrogen fertilization, growing season extensions and release from cold limitations), such processes are included in the models examined and contribute to the scatter in the relationship between β GPP R and S LAND presented in Fig. 2a.That said, there are many processes inadequately represented in both TBMs and the satellite approaches that could lead to biases in the derived β GPP R .For example, models have been shown to poorly reproduce changes in the seasonal cycle of atmospheric CO 2 (ref.41) and demonstrate a range of responses when compared to results from experimental manipulation 42 .Nutrient limitation, thermal temperature acclimation, water stress, disturbances (including land-use change) and leaf area dynamics are all poorly represented in TBMs 42,43 .Future implementations of new process representations or model structures may lead to updated inference on the response of photosynthesis to CO 2 .
A further source of uncertainty relates to the degree of structural similarity between models and the potential for systematic cross-model biases.For example, if all models in the ensemble had the same missing or biased process representation, which led to systematic bias in the modelled relationship between the sensitivity of photosynthesis to CO 2 and the land sink across models, that could bias the emergent constraint reported here.Systematic cross-model biases with shared structural similarity could also lead to an underestimation of the uncertainty associated with the values derived from the emergent constraint 40,44 .
The models we examine represent the state-of-the-science for land surface modelling and have substantial diversity of process representations and responses to forcings 45 , even for well-studied processes such as photosynthesis.Because of this diversity, there are outlier models with high or low CO 2 sensitivities or S LAND estimates and such 'wrong' models are necessary for the formation of an emergent constraint.Indeed, removal of some outlier models, in particular CABLE and VEGAS, degrades the derived relationship between β GPP R and S LAND presented in Fig. 2a (to r = 0.5, 0.64; P = 0.03, 0.01, respectively) and removal of both models leads to no statistically significant relationship between β GPP R and S LAND (P = 0.15).If future versions of current outlier models are more consistent with the ensemble, the constraint identified here may no longer be evident.
Global photosynthesis is the largest flux of CO 2 in the global carbon cycle and small changes in terrestrial photosynthesis over time can lead to large changes in the net carbon sink.The resulting feedback from the ) effect of increasing CO 2 on GPP, in addition to the effect of land use and climate changes on the fraction of absorbed radiation (fAPAR) (c); the difference between the data presented in b and c (d); the latitudinal distribution of long-term changes in gross primary photosynthesis (∆GPP, PgC) from 1982 to 2012, from the TBM ensemble (orange shaded area, mean, s.d.across models) and ∆GPP predicted from RS approaches with (black, mean, s.d. between MODIS and ML approaches) and without (red) a direct effect of CO 2 on LUE (Methods) (e); and long-term changes in ∆GPP, separated by plant functional types (f).EBF, evergreen broadleaved forest; SAV, savanna; DBF, deciduous broadleaved forests; CRO, croplands; SH, shrublands; ENF, evergreen needleleaf forests; GRA, grasslands; WET, wetlands.

Article
https://doi.org/10.1038/s41558-023-01867-2effect of increasing CO 2 on photosynthesis (the carbon-concentration feedback) has been estimated to be over four times larger and more uncertain, than the direct carbon-climate feedback 46 .The large differences between estimates of historic changes in GPP [6][7][8][9][10]15 are therefore disconcerting and could potentially lead to incorrect inference about biases in current TBMs 6,14 and long-term changes in related components of the global carbon cycle such as soil respiration 11,47 . Th confluence of approaches we use bounds the plausible range of the historic effect of CO 2 on global terrestrial photosynthesis to a β GPP R of 0.62 ± 0.14 (mean, s.d.; Fig. 2b) and helps to reconcile differences in previous estimates.The results also show that widely used RS-based estimates of global terrestrial photosynthesis need to incorporate the effect of increasing CO 2 on photosynthetic LUE and provide a globally applicable approach that is broadly consistent with the TBMs examined.Together, our results suggest that increases in atmospheric CO 2 have led to a large increase in global photosynthesis since 1982, representing a strong carbon-concentration feedback that has helped to slow down the accumulation of anthropogenic emissions in the atmosphere.

Methods
The β metric of CO 2 sensitivity We quantified the apparent sensitivity of global terrestrial GPP to CO 2 in the RS, TBM and independent proxy estimates using two approaches: (1) the percentage change in GPP with respect to GPP at the start of the time period (equation ( 5) below) and (2) a β metric defined as the response ratio (R) of GPP with respect to CO 2 : where GPP(t) is the value of GPP at time t and Ca(t) is the value of atmospheric [CO 2 ] at time t.Although other methods to calculate the β-factor have been proposed (for example, ref. 48), we use equation ( 1) for ease of interpretation.A β of 1 represents direct proportionality between the GPP CO 2 response and the change in CO 2 .Note that to avoid undue influence of year-to-year variability in GPP, we estimated GPP(t) and GPP(t 0 ) on the basis of a linear regression fit to the GPP time series.

Assessing the CO 2 -sensitivity of satellite-based GPP
Recent reports have highlighted that the most commonly used satellite-based estimates of GPP have a much lower CO 2 -sensitivity than that derived from TBMs 6,11 .However, most satellite-based estimates do not incorporate the universally observed direct effect of increasing CO 2 on the LUE of leaves of C 3 vegetation 13 , which is not observable from space 14 .The effect of increasing CO 2 on global terrestrial C 3 photosynthesis that we examine here manifests through two primary pathways: though increasing the biochemical rate of photosynthesis on a leaf area basis 49 , which we refer to as the direct effect and through increases in leaf area on a ground area basis, allowing for the interception of greater amounts of light 50,51 , which we refer to as the indirect effect.The former, direct response, arises because CO 2 is a substrate for the photosynthetic enzyme, RuBisCO.Both CO 2 and O 2 compete at the active site of RuBisCO, so changes in the concentration of either affect the rate at which CO 2 is assimilated, effectively changing the LUE of photosynthesis on a leaf area basis at a given light level.The latter, indirect response of increasing leaf area index (LAI 51 ) and the resulting increase in the (fAPAR), reflects both the increased carbon available to invest in structural growth under elevated CO 2 and potential changes in the hydrological equilibrium due to elevated CO 2 -induced increases in water-use efficiency, which can lead to increased leaf area in water-limited ecosystems [52][53][54] .Both response pathways are incorporated in TBMs 24 and long-term proxies account for each to differing degrees.Most satellite-based estimates, however, do not account for the direct effect of increasing CO 2 on the biochemical rate of photosynthesis 14,55 .
We assessed whether incorporating a CO 2 sensitivity of LUE in RS-based approaches for estimating GPP reconciled the difference between the sensitivity of RS-based GPP to increasing CO 2 and that implied by the emergent constraint.To do so, we develop a CO 2 sensitivity function for incorporating the effect of increasing CO 2 on the LUE of photosynthesis into satellite GPP estimates, based on the conservative assumption that the ecosystem-scale CO 2 sensitivity is consistent with that of the electron-transport limited rate of photosynthesis (A j ).This is supported by reports that the observed CO 2 response of photosynthesis and biomass closely corresponds to the CO 2 -sensitivity of A j (ref.35).In addition, it has been suggested that shaded, and thus primarily electron-transport limited, leaves contribute the most canopy 36,56 and global terrestrial photosynthesis 39 .The assumption is further supported by optimal coordination theory, which posits that photosynthesis under typical daytime field conditions is close to the point where RuBisCO-limited (A c ) and A j are colimiting.The colimitation of A c and A j has been shown to hold across a range of ecosystems 57 , as has the downregulation of the maximum velocity of carboxylation (V cmax ) under elevated CO 2 to maintain coordination 58 .Given that the sensitivity of A j to CO 2 is much smaller than that of A c (ref.59), the sensitivity of A j to CO 2 therefore represents a conservative approach to incorporate a CO 2 sensitivity of LUE 37 in RS estimates of photosynthesis.Note that we also make the conservative assumption that C 4 plants operate at or near CO 2 saturation 60 .
The mechanistic photosynthesis model proposed by ref. 49 captures the biochemical controls of leaf photosynthesis and responses to variations in temperature, light and CO 2 concentration.According to the model, the gross photosynthesis rate, A, is limited by either the capacity of the RuBisCO enzyme for the carboxylation of ribulose-1,5-bisphosphate (RuBP), the electron-transport capacity for RuBP regeneration.In the case of the limitation by the electron-transport capacity for RuBP regeneration, the photosynthetic rate (A j , µmol m −2 s −1 ) is given by: where φ 0 is the intrinsic quantum efficiency, I is the absorbed light (µmol m −2 s −1 ), c i (Pa) is the leaf-internal CO 2 concentration and Γ * (Pa) is the CO 2 compensation point.Parameter Γ * depends on temperature, as estimated through a biochemical rate parameter (r) 61 : where R is the molar gas constant (8.314J mol −1 K −1 ), r 25 = 4.22 Pa, is the photorespiratory point at 25 ˚C, ΔH is the activation energy for Γ * (37.83 kJ mol −1 ) and T is the temperature in K. Assuming the CO 2 sensitivity of light-limited photosynthesis allows for the development of an index of the effect of CO 2 on photosynthetic LUE 37 , which can be incorporated in any RS-based LUE model or empirical upscaling estimate of GPP.By rewriting equation ( 2), substituting c i by the product of atmospheric CO 2 (c a ) and the ratio of leaf-internal to leaf-ambient CO 2 (χ = c i /c a ), the sensitivity of GPP and LUE to CO 2 can be described as: where ϕ CO 2 = c a χ−Γ * c a χ+2Γ * and LUE = GPP/φ 0 I.Note that the indirect effect of CO 2 on GPP through φ 0 I, is explicitly accounted for in satellite-based methods through changes in the fAPAR and considered here as an independent effect.However, the direct effect, through changes in LUE, (ϕ CO 2 ), is not.We used equation ( 4) to derive a scalar, f (CO 2 ), to account for the direct effect of CO 2 in any LUE-based estimate of GPP (for example, satellite or empirical upscaling approaches).To do so, we calculated ∆GPP in year t due to the effect of CO 2 on LUE as GPP(t = 0) × f (CO 2 ), where: f (CO 2 ) thus represents the fractional increase in LUE due to the direct effect of CO 2 relative to a baseline period (here 1982, the start of the time series for the satellite-based methods considered).
The sensitivity of LUE to CO 2 thus depends on both Γ * , which is calculated by means of equation ( 3), and χ.We estimated χ using the least-cost hypothesis 19,62 .This states that an optimal long-term effective value of χ can be predicted as a result of plants minimizing their total https://doi.org/10.1038/s41558-023-01867-2carbon costs associated with photosynthetic carbon gain and explicitly expressed with the following model: where D is vapour pressure deficit and η * is the viscosity of water relative to its value at 25 °C (ref.63) and b is the ratio of the cost of maintaining carboxylation relative to that of maintaining transpiration 19 .The Michaelis-Menten coefficient of RuBisCO (K) is given by: where  60 .We thus make the conservative assumption of no direct CO 2 effect on LUE in the C 4 proportion of each pixel.

Incorporating a CO 2 sensitivity into satellite-based GPP
The approach for incorporating a CO 2 sensitivity we outline above (equation ( 5)) can be incorporated into any satellite-based photosynthesis product.Here, we test it on two broadly used approaches.The first, the MODIS MOD17 algorithm (GPP MODIS (ref.65)) and the second an empirical upscaling method based on a model tree ensemble (GPP MTE (ref.66)).We applied the MODIS MOD17 GPP algorithm driven by 30-year (1982-2012) Global Inventory Modeling and Mapping Studies (GIMMS3g) fAPAR data 67 , to calculate a new 30-year global monthly gridded (0.5°) dataset of MODIS-derived GPP: where  66 ), which does not account for the direct effect of CO 2 on LUE, we followed the approach outlined for the MODIS GPP product.Specifically, we applied the CO 2 function (equation ( 5) to spatially distributed GPP MTE , as: Early RS GPP models 37,70 advocated for including a CO 2 effect on LUE, though primarily used the larger, light-saturated, sensitivity.A recent review 7 found that the most widely used modern RS GPP approaches 65,66 do not include a CO 2 effect on LUE and of the 3 that did (out of 14 assessed) 2 are enzyme kinetics, not LUE, models (BESS 72 and BEPS 70 ).The third (cFix 71 ) assumes the light-saturated CO 2 sensitivity, which is not suitable for global application given the large contribution of RuBP regeneration-limited leaves 36,73 .Some recent studies 12,74,75 incorporated a CO 2 effect on LUE but the approach taken typically requires the reparameterization of the LUE model and is thus not easily applicable to other RS GPP products.The approach proposed here provides a generic and conservative method for incorporating CO 2 effects on LUE in any RS GPP product, which allows us to quantify the relative importance of incorporating a CO 2 effect in RS GPP products and reconciles the large difference between RS and TBM-derived sensitivities to CO 2 .

Constraining terrestrial photosynthesis CO 2 sensitivity
Emergent constraints have gained prominence in recent years as a means by which to infer unobserved quantities of interest in land surface and climate models [21][22][23] .The underlying core concept is that, although there is a large spread in the model estimates of an observed variable X and an unobserved variable Y across models, the relationship linking the two is sometimes tightly constrained across models.Given the existence of a strong and robust relationship across models between X and Y, observations of X can be used to generate a probabilistic inference, or constraint, on Y.This approach has been termed 'emergent' because the functional relationship cannot be diagnosed from a single model but rather emerges from examining the model spread [21][22][23] .
The emergent constraint identified in this study links the sensitivity of GPP to CO 2 (β GPP R , see definition below) to the magnitude of the cumulative residual terrestrial sink (S LAND ) between 1982 and 2016.It is derived from a linear regression across an ensemble of TBMs between the modelled cumulative S LAND and the sensitivity of GPP to CO 2 .We use global simulations from 15 TBMs (Supplementary Table 2) run as part of the Trends in Net Land-Atmosphere Exchange (TRENDY v.6) initiative (https://sites.exeter.ac.uk/trendy) (v.6 data are reported in ref. 76).In TRENDY, common input forcing data were prescribed for a series of model experiments from 1901 to 2015.Here we use both the results of the TRENDY v.6 scenario S3 simulations (temporally dynamic climate, CO 2 and land use) as reported in the Global Carbon Project (GCP 76 ) and the TRENDY v.6 scenario S1 simulations (CO 2 -only: temporally dynamic CO 2 , time-invariant climate; pre-industrial land-use mask).For more details on the TRENDY project see ref. 24 and for details of the TRENDY v.6 simulations used here see ref. 76.
We estimated β GPP R for each TRENDY v.6 TBM from annual GPP from the S1 (CO 2 -only) simulations, performed by 15 models (Supplementary Table 2), using equation (1) over the 1982-2012 period (to maintain consistency with the RS methods assessed).Cumulative S LAND (PgC) is calculated from the annual S LAND (PgC yr −1 ) reported by the GCP 76 for each TRENDY v.6 TBM, which represents the annual total net biome productivity plus emissions from land-use change.
The emergent constraint approach relies on a statistical relationship between a model predicted variable for which an observational constraint exists and one for which there is no observational constraint available [21][22][23] .In the case of the relationship between β GPP R and S LAND , estimates of S LAND are made annually by the Global Carbon Project, along with the associated uncertainties 25 .The S LAND values we use as the constraint are the cumulative reported annual values of the residual land sink from the Global Carbon Project 25 over the period from 1982 to 2016.Note that the period we used was chosen to both coincide with the satellite observations we use and to be sufficiently long so as to minimize the effect of macroclimatic events such as strong El Nino periods and volcanic eruptions.
The Global Carbon Project reports S LAND uncertainty both on an annual, decadal and a cumulative basis, with an average uncertainty of 0.9 PgC yr −1 for each of the four decades included in this study.For the https://doi.org/10.1038/s41558-023-01867-21960-2020 period where direct atmospheric CO 2 measurements are available, the Global Carbon Project estimates residual land carbon uptake of 135 ± 25 GtC (mean ± s.d.; ref. 25, Table 8), with a near-zero unattributed budget imbalance.The budget closure is interpreted as evidence of a coherent community understanding of the emissions and global sinks for this period 25 .This provides a cumulative S LAND reference uncertainty of 18.5%, which we apply to the cumulative fluxes of the period examined (1982 to 2016).It should be noted that the uncertainty on cumulative S LAND is itself uncertain and is estimated by the Global Carbon Project based on the most up-to-date versions of the land surface models they use 25 .Any future reduction in cumulative S LAND uncertainty would decrease the uncertainty of β GPP R reported here.

Fig. 2 |
Fig. 2 | A constraint on the sensitivity of global terrestrial photosynthesis to CO 2 .a, The relationship between the modelled sensitivity of GPP to CO 2 ( β GPP R , TRENDY experiment S1: dynamic CO 2 only) and the modelled cumulative terrestrial carbon sink (PgC, TRENDY experiment S3: dynamic CO 2 , climate and land use).Individual TRENDY model details and β GPP R values are listed in Supplementary Table2.The red line and shaded area show the best linear fit across models and the associated prediction standard error (dashed red) and standard deviation of prediction error (dashed grey) intervals.The vertical Fig. 2 | A constraint on the sensitivity of global terrestrial photosynthesis to CO 2 .a, The relationship between the modelled sensitivity of GPP to CO 2 ( β GPP R , TRENDY experiment S1: dynamic CO 2 only) and the modelled cumulative terrestrial carbon sink (PgC, TRENDY experiment S3: dynamic CO 2 , climate and land use).Individual TRENDY model details and β GPP R values are listed in Supplementary Table2.The red line and shaded area show the best linear fit across models and the associated prediction standard error (dashed red) and standard deviation of prediction error (dashed grey) intervals.The vertical Fig. 2 | A constraint on the sensitivity of global terrestrial photosynthesis to CO 2 .a, The relationship between the modelled sensitivity of GPP to CO 2 ( β GPP R , TRENDY experiment S1: dynamic CO 2 only) and the modelled cumulative terrestrial carbon sink (PgC, TRENDY experiment S3: dynamic CO 2 , climate and land use).Individual TRENDY model details and β GPP R values are listed in Supplementary Table2.The red line and shaded area show the best linear fit across models and the associated prediction standard error (dashed red) and standard deviation of prediction error (dashed grey) intervals.The vertical Article https://doi.org/10.1038/s41558-023-01867-2

Fig. 3 |
Fig. 3 | Spatial differences in the estimated long-term changes in global terrestrial photosynthesis from LUE theory, TBMs and satellite observations combined with theory.a-f, The global distribution of: the sensitivity of terrestrial photosynthesis on a leaf area basis to CO 2 (β LUE R ) due to changes in LUE (a); CO 2 -induced changes in terrestrial photosynthesis (∆GPP, gC m −2 yr −1 ) from 1982 to 2012 from an ensemble of TBMs (TBMs; TRENDY-S1) (b); mean long-term changes in GPP from the two updated satellite methods, which includes a modelled direct (β LUE R ) and measured indirect ( β fAPAR R

1 | Long-term changes in global annual photosynthesis from TBMs and multiple
based on simulations from process-based models in the TRENDY project model ensemble (orange, mean ± s.d.) and two different satellite approaches (empirical MODIS algorithm (MA, solid lines); a machine learning method (ML, dashed lines)).Estimates from the satellite approaches were obtained allowing for an effect of increasing CO 2 on either the fAPAR (red lines, dots), the LUE of 64c and K o are the Michaelis-Menten coefficient of RuBisCO for carboxylation and oxygenation, respectively, expressed in partial pressure units and P o is the partial pressure of O 2 .K responds to temperature through K c and K o , the temperature responses for which are described using a temperature response function described by equation (3) with specific parameters: ΔH is 79.43 kJ mol −1 for K c and 36.38 kJ mol −1 for K o , r 25 is 39.97 kPa for K c and 27.48 kPa for K o (ref.61).We applied this derived sensitivity to the RS approaches detailed below, on a per-pixel basis in proportion to the percentage of C 3 plants in a given pixel64, as C 4 plants operate at or near CO 2 saturation LUE max represents biome-specific maximum LUE, f(D) represents a water stress reduction scalar based on the atmospheric vapour pressure deficit and f (T min ) represents a low-temperature stress reduction scalar.LUE max , f(D) and f (T min ) are parameterized according to ref.68.Value f (CO 2 ) is estimated on a per-pixel basis using equation(5).