The land sink of anthropogenic carbon emissions, a crucial component of mitigating climate change, is primarily attributed to the CO₂ fertilization effect on global gross primary productivity (GPP). However, direct observational evidence of this effect remains scarce, hampered by challenges in disentangling the CO₂ fertilization effect from other long-term drivers, particularly climatic changes. Here, we introduce a novel statistical approach to separate the CO₂ fertilization effect on GPP and daily maximum net ecosystem production (NEPmax) using eddy covariance records across 38 extratropical forest sites. We find the median stimulation rate of GPP and NEPmax to be 16.4 ± 4% and 17.2 ± 4% per 100 ppm increase in atmospheric CO₂ across these sites, respectively. To validate the robustness of our findings, we test our statistical method using factorial simulations of an ensemble of process-based land surface models. We acknowledge that additional factors, including nitrogen deposition and land management, may impact plant productivity, potentially confounding the attribution to the CO₂ fertilization effect. Assuming these site-specific effects offset to some extent across sites as random factors, the estimated median value still reflects the strength of the CO₂ fertilization effect. However, disentanglement of these long-term effects, often inseparable by timescale, requires further causal research. Our study provides direct evidence that the photosynthetic stimulation is maintained under long-term CO₂ fertilization across multiple eddy covariance sites. Such observation-based quantification is key to constraining the long-standing uncertainties in the land carbon cycle under rising CO₂ concentrations.

Elevated atmospheric CO2 (eCO2) influences the carbon assimilation rate and stomatal conductance of plants, and thereby can affect the global cycles of carbon and water. However, the extent to which these physiological effects of eCO2 influence the land-atmosphere exchange of carbon and water is uncertain. In this study, we aim at developing a method to detect the emergence of the physiological CO2 effects on various variables related to carbon and water fluxes. We use a comprehensive process-based land surface model QUINCY (QUantifying Interactions between terrestrial Nutrient CYcles and the climate system) to simulate the leaf-level effects of increasing atmospheric CO2 concentrations and their century-long propagation through the terrestrial carbon and water cycles across different climate regimes and biomes. We then develop a statistical method based on the signal-to-noise ratio to detect the emergence of the eCO2 effects. The signal in gross primary production (GPP) emerges at relatively low eCO2 (Δ[CO2] ~ 20 ppm) where the leaf area index (LAI) is relatively high. Compared to GPP, the eCO2 effect causing reduced 28 transpiration water flux (normalized to leaf area) emerges only at relatively high CO2 increase (Δ[CO2] >> 40 ppm), due to the high sensitivity to climate variability and thus lower signal-to-noise ratio. In general, the response to eCO2 is detectable earlier for variables of the carbon cycle than the water cycle, when plant productivity is not limited by climatic constraints, and stronger in forest-dominated rather than in grass- dominated ecosystems. Our results provide a step towards when and where we expect to detect physiological CO2 effects in in-situ flux measurements, how to detect them and encourage future efforts to improve the understanding and quantification of these effects in observations of terrestrial carbon and water dynamics.