Search NASA⌕ Search

Engineering topics

Stephen Sitch

Publications and source records attributed to Stephen Sitch.

18 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Soil Respiration–Driven CO 2 Pulses Dominate Australia’s Flux Variability

The Australian continent contributes substantially to the year-to-year variability of the global terrestrial carbon dioxide (CO 2 ) sink. However, the scarcity of in situ observations in remote areas prevents the deciphering of processes that force the CO 2 flux variability. In this study, by examining atmospheric CO 2 measurements from satellites in the period 2009–2018, we find recurrent end-of-dry-season CO 2 pulses over the Australian continent. These pulses largely control the year-to-year variability of Australia’s CO 2 balance. They cause two to three times larger seasonal variations compared with previous top-down inversions and bottom-up estimates. The pulses occur shortly after the onset of rainfall and are driven by enhanced soil respiration preceding photosynthetic uptake in Australia’s semiarid regions. The suggested continental-scale relevance of soil-rewetting processes has substantial implications for our understanding and modeling of global climate–carbon cycle feedbacks.

Eva-Marie Metz↗
NASA NTRS2023

Harmonising the Land-Use Flux Estimates of Global Models and National Inventories for 2000–2020

As the focus of climate policy shifts from pledges to implementation, there is a growing need to track progress on climate change mitigation at the country level, particularly for the land-use sector. Despite new tools and models providing unprecedented monitoring opportunities, striking differences remain in estimations of anthropogenic land-use CO 2 fluxes between, on the one hand, the national greenhouse gas inventories (NGHGIs) used to assess compliance with national climate targets under the Paris Agreement and, on the other hand, the Global Carbon Budget and Intergovernmental Panel on Climate Change (IPCC) assessment reports, both based on global bookkeeping models (BMs). Recent studies have shown that these differences are mainly due to inconsistent definitions of anthropogenic CO 2 fluxes in managed forests. Countries assume larger areas of forest to be managed than BMs do, due to a broader definition of managed land in NGHGIs. Additionally, the fraction of the land sink caused by indirect effects of human-induced environmental change (e.g. fertilisation effect on vegetation growth due to increased atmospheric CO 2 concentration) on managed lands is treated as non-anthropogenic by BMs but as anthropogenic in most NGHGIs. We implement an approach that adds the CO 2 sink caused by environmental change in countries' managed forests (estimated by 16 dynamic global vegetation models, DGVMs) to the land-use fluxes from three BMs. This sum is conceptually more comparable to NGHGIs and is thus expected to be quantitatively more similar. Our analysis uses updated and more comprehensive data from NGHGIs than previous studies and provides model results at a greater level of disaggregation in terms of regions, countries and land categories (i.e. forest land, deforestation, organic soils, other land uses). Our results confirm a large difference (6.7 GtCO 2 yr −1 ) in global land-use CO 2 fluxes between the ensemble mean of the BMs, which estimate a source of 4.8 GtCO 2 yr −1 for the period 2000–2020, and NGHGIs, which estimate a sink of −1.9 GtCO 2 yr −1 in the same period. Most of the gap is found on forest land (3.5 GtCO 2 yr −1 ), with differences also for deforestation (2.4 GtCO 2 yr −1 ), for fluxes from other land uses (1.0 GtCO 2 yr −1 ) and to a lesser extent for fluxes from organic soils (0.2 GtCO2 yr−1). By adding the DGVM ensemble mean sink arising from environmental change in managed forests (−6.4 GtCO 2 yr −1 ) to BM estimates, the gap between BMs and NGHGIs becomes substantially smaller both globally (residual gap: 0.3 GtCO 2 yr −1 ) and in most regions and countries. However, some discrepancies remain and deserve further investigation. For example, the BMs generally provide higher emissions from deforestation than NGHGIs and, when adjusted with the sink in managed forests estimated by DGVMs, yield a sink that is often greater than NGHGIs. In summary, this study provides a blueprint for harmonising the estimations of anthropogenic land-use fluxes, allowing for detailed comparisons between global models and national inventories at global, regional and country levels. This is crucial to increase confidence in land-use emissions estimates, support investments in land-based mitigation strategies and assess the countries' collective progress under the Global Stocktake of the Paris Agreement.

Giacomo Grassi↗
NASA NTRS2022

Global Carbon Budget 2022

Accurate assessment of anthropogenic carbon dioxide (CO 2 ) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO 2 emissions (E FOS ) are based on energy statistics and cement production data, while emissions from land-use change (E LUC ), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO 2 concentration is measured directly, and its growth rate (G ATM ) is computed from the annual changes in concentration. The ocean CO 2 sink (S OCEAN ) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO 2 sink (S LAND ) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (B IM ), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2021, E FOS increased by 5.1 % relative to 2020, with fossil emissions at 10.1 ± 0.5 GtC yr −1 (9.9 ± 0.5 GtC yr −1 when the cement carbonation sink is included), and E LUC was 1.1 ± 0.7 GtC yr −1 , for a total anthropogenic CO 2 emission (including the cement carbonation sink) of 10.9 ± 0.8 GtC yr −1 (40.0 ± 2.9 GtCO 2 ). Also, for 2021, G ATM was 5.2 ± 0.2 GtC yr −1 (2.5 ± 0.1 ppm yr −1 ), S OCEAN was 2.9 ± 0.4 GtC yr −1 , and S LAND was 3.5 ± 0.9 GtC yr −1 , with a B IM of −0.6 GtC yr −1 (i.e. the total estimated sources were too low or sinks were too high). The global atmospheric CO 2 concentration averaged over 2021 reached 414.71 ± 0.1 ppm. Preliminary data for 2022 suggest an increase in E FOS relative to 2021 of +1.0 % (0.1 % to 1.9 %) globally and atmospheric CO 2 concentration reaching 417.2 ppm, more than 50 % above pre-industrial levels (around 278 ppm). Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2021, but discrepancies of up to 1 GtC yr −1 persist for the representation of annual to semi-decadal variability in CO 2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) a low agreement between the different methods on the magnitude of the land CO 2 flux in the northern extratropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set. The data presented in this work are available at https://doi.org/10.18160/GCP-2022 (Friedlingstein et al., 2022b).

Pierre Friedlingstein↗
NASA NTRS2022

The Ozone–Climate Penalty Over South America and Africa By 2100

Climate change has the potential to increase surface ozone (O3) concentrations, known as the “ozone–climate penalty”, through changes to atmospheric chemistry, transport and dry deposition. In the tropics, the response of surface O3 to changing climate is relatively understudied but has important consequences for air pollution and human and ecosystem health. In this study, we evaluate the change in surface O3 due to climate change over South America and Africa using three state-of-the-art Earth system models that follow the Shared Socioeconomic Pathway 3-7.0 emission scenario from CMIP6. In order to quantify changes due to climate change alone, we evaluate the difference between simulations including climate change and simulations with a fixed present-day climate. We find that by 2100, models predict an ozone–climate penalty in areas where O3 is already predicted to be high due to the impacts of precursor emissions, namely urban and biomass burning areas, although on average, models predict a decrease in surface O3 due to climate change. We identify a small but robust positive trend in annual mean surface O3 over polluted areas. Additionally, during biomass burning seasons, seasonal mean O3 concentrations increase by 15 ppb (model range 12 to 18 ppb) in areas with substantial biomass burning such as the arc of deforestation in the Amazon. The ozone–climate penalty in polluted areas is shown to be driven by an increased rate of O3 chemical production, which is strongly influenced by NOx concentrations and is therefore specific to the emission pathway chosen. Multiple linear regression finds the change in NOx concentration to be a strong predictor of the change in O3 production, whereas increased isoprene emission rate is positively correlated with increased O3 destruction, suggesting NOx-limited conditions over the majority of tropical Africa and South America. However, models disagree on the role of climate change in remote, low-NOx regions, partly because of significant differences in NOx concentrations produced by each model. We also find that the magnitude and location of the ozone–climate penalty in the Congo Basin has greater inter-model variation than that in the Amazon, so further model development and validation are needed to constrain the response in central Africa. We conclude that if the climate were to change according to the emission scenario used here, models predict that forested areas in biomass burning locations and urban populations will be at increasing risk of high O3 exposure, irrespective of any direct impacts on O3 via the prescribed emission scenario.

Surface ozone↗
NASA NTRS2022

A Process-Model Perspective on Recent Changes in the Carbon Cycle of North America

Continental North America has been found to be a carbon (C) sink over recent decades by multiple studies employing a variety of estimation approaches. However, several key questions and uncertainties remain with these assessments. Here we used results from an ensemble of 19 state-of-the-art dynamic global vegetation models from the TRENDYv9 project to improve these estimates and study the drivers of its interannual variability. Our results show that North America has been a C sink with a magnitude of 0.37 ± 0.38 (mean and one standard deviation) PgC year−1 for the period 2000–2019 (0.31 and 0.44 PgC year−1 in each decade); split into 0.18 ± 0.12 PgC year−1 in Canada (0.15 and 0.20), 0.16 ± 0.17 in the United States (0.14 and 0.17), 0.02 ± 0.05 PgC year−1 in Mexico (0.02 and 0.02) and 0.01 ± 0.02 in Central America and the Caribbean (0.01 and 0.01). About 57% of the new C assimilated by terrestrial ecosystems is allocated into vegetation, 30% into soils, and 13% into litter. Losses of C due to fire account for 41% of the interannual variability of the mean net biome productivity for all North America in the model ensemble. Finally, we show that drought years (e.g., 2002) have the potential to shift the region to a small net C source in the simulations (−0.02 ± 0.46 PgC year−1). Our results highlight the importance of identifying the major drivers of the interannual variability of the continental-scale land C cycle along with the spatial distribution of local sink-source dynamics.

Guillermo Murray-Tortarolo↗
NASA NTRS2022

Process-Oriented Analysis of Dominant Sources of Uncertainty in the Land Carbon Sink

The observed global net land carbon sink is captured by current land models. All models agree that atmospheric CO 2 and nitrogen deposition driven gains in carbon stocks are partially offset by climate and land-use and land-cover change (LULCC) losses. However, there is a lack of consensus in the partitioning of the sink between vegetation and soil, where models do not even agree on the direction of change in carbon stocks over the past 60 years. This uncertainty is driven by plant productivity, allocation, and turnover response to atmospheric CO 2 (and to a smaller extent to LULCC), and the response of soil to LULCC (and to a lesser extent climate). Overall, differences in turnover explain ~70% of model spread in both vegetation and soil carbon changes. Further analysis of internal plant and soil (individual pools) cycling is needed to reduce uncertainty in the controlling processes behind the global land carbon sink.

Attribution↗
NASA NTRS2021

Are Land-Use Change Emissions in Southeast Asia Decreasing or Increasing?

Southeast Asia is a region known for active land-use changes (LUC) over the past 60 years; yet, how trends in net CO2 uptake and release resulting from LUC activities (net LUC flux) have changed through past decades remains uncertain. The level of uncertainty in net LUC flux from process-based models is so high that it cannot be concluded that newer estimates are necessarily more reliable than older ones. Here, we examined net LUC flux estimates of Southeast Asia for the 1980s−2010s from older and newer sets of Dynamic Global Vegetation Model simulations (TRENDY v2 and v7, respectively), and forcing data used for running those simulations, along with two book-keeping estimates (H&N and BLUE). These estimates yielded two contrasting historical LUC transitions, such that TRENDY v2 and H&N showed a transition from increased emissions from the 1980s to 1990s to declining emissions in the 2000s, while TRENDY v7 and BLUE showed the opposite transition. We found that these contrasting transitions originated in the update of LUC forcing data, which reduced the loss of forest area during the 1990s. Further evaluation of remote sensing studies, atmospheric inversions, and the history of forestry and environmental policies in Southeast Asia supported the occurrence of peak emissions in the 1990s and declining thereafter. However, whether LUC emissions continue to decline in Southeast Asia remains uncertain as key processes in recent years, such as conversion of peat forest to oil-palm plantation, are yet to be represented in the forcing data, suggesting a need for further revision.

land-use change↗
NASA NTRS2021

Response to Comments on “Recent global decline of CO2 fertilization effects on vegetation photosynthesis”

Our study suggests that the global CO2 fertilization effect (CFE) on vegetation photosynthesis has declined during the past four decades. The Comments suggest that the temporal inconsistency in AVHRR data and the attribution method undermine the results’ robustness. Here, we provide additional evidence that these arguments did not affect our finding and that the global decline in CFE is robust.

evapotranspiration↗
NASA NTRS2021

Modelled land use and land cover change emissions – a spatio-temporal comparison of different approaches

Quantifying the net carbon flux from land use and land cover changes (fLULCC) is critical for understanding the global carbon cycle and, hence, to support climate change mitigation. However, large-scale fLULCC is not directly measurable and has to be inferred from models instead, such as semi-empirical bookkeeping models and process-based dynamic global vegetation models (DGVMs). By definition, fLULCC estimates are not directly comparable between these two different model types. As an important example, DGVM-based fLULCC in the annual global carbon budgets is estimated under transient environmental forcing and includes the so-called loss of additional sink capacity (LASC). The LASC results from the impact of environmental changes on land carbon storage potential of managed land compared to potential vegetation and accumulates over time, which is not captured in bookkeeping models. The fLULCC from transient DGVM simulations, thus, strongly depends on the timing of land use and land cover changes mainly because LASC accumulation is cut off at the end of the simulated period. To estimate the LASC, the fLULCC from pre-industrial DGVM simulations, which is independent of changing environmental conditions, can be used. Additionally, DGVMs using constant present-day environmental forcing enable an approximation of bookkeeping estimates. Here, we analyse these three DGVM-derived fLULCC estimations (under transient, pre-industrial, and present-day forcing) for 12 models within 18 regions and quantify their differences as well as climate- and CO2-induced components and compare them to bookkeeping estimates. Averaged across the models, we find a global fLULCC (under transient conditions) of 2.0±0.6 PgC/yr for 2009–2018, of which ∼40 % are attributable to the LASC (0.8±0.3 PgC/yr). From 1850 onward, the fLULCC accumulated to 189±56 PgC with 40±15 PgC from the LASC. Around 1960, the accumulating nature of the LASC causes global transient fLULCC estimates to exceed estimates under present-day conditions, despite generally increased carbon stocks in the latter. Regional hotspots of high cumulative and annual LASC values are found in the USA, China, Brazil, equatorial Africa, and Southeast Asia, mainly due to deforestation for cropland. Distinct negative LASC estimates in Europe (early reforestation) and from 2000 onward in the Ukraine (recultivation of post-Soviet abandoned agricultural land), indicate that fLULCC estimates in these regions are lower in transient DGVM compared to bookkeeping approaches. Our study unravels the strong dependence of fLULCC estimates on the time a certain land use and land cover change event happened to occur and on the chosen time period for the forcing of environmental conditions in the underlying simulations. We argue for an approach that provides an accounting of the fLULCC that is more robust against these choices, for example by estimating a mean DGVM ensemble fLULCC and LASC for a defined reference period and homogeneous environmental changes (CO2 only).

land use and land cover change↗
NASA NTRS2020

Causes of slowing‐down seasonal CO2 amplitude at Mauna Loa

Changing amplitude of the seasonal cycle of atmospheric CO2 (SCA) in the northern hemisphere is an emerging carbon cycle property. Mauna Loa (MLO) station (20°N, 156°W), which has the longest continuous northern hemisphere CO2 record, shows an increasing SCA before the 1980s (p < .01), followed by no significant change thereafter. We analyzed the potential driving factors of SCA slowing‐down, with an ensemble of dynamic global vegetation models (DGVMs) coupled with an atmospheric transport model. We found that slowing‐down of SCA at MLO is primarily explained by response of net biome productivity (NBP) to climate change, and by changes in atmospheric circulations. Through NBP, climate change increases SCA at MLO before the 1980s and decreases it afterwards. The effect of climate change on the slowing‐down of SCA at MLO is mainly exerted by intensified drought stress acting to offset the acceleration driven by CO2 fertilization. This challenges the view that CO2 fertilization is the dominant cause of emergent SCA trends at northern sites south of 40°N. The contribution of agricultural intensification on the deceleration of SCA at MLO was elusive according to land–atmosphere CO2 flux estimated by DGVMs and atmospheric inversions. Our results also show the necessity to adequately account for changing circulation patterns in understanding carbon cycle dynamics observed from atmospheric observations and in using these observations to benchmark DGVMs.

Mauna Loa↗
NASA NTRS2020

Evaluation of global terrestrial evapotranspiration using state-of-the-art approaches in remote sensing, machine learning and land surface modeling

Evapotranspiration (ET) is critical in linking global water, carbon and energy cycles. However, direct measurement of global terrestrial ET is not feasible. Here, we first reviewed the basic theory and state-of-the-art approaches for estimating global terrestrial ET, including remote-sensing-based physical models, machine-learning algorithms and land surface models (LSMs). We then utilized 4 remote-sensing-based physical models, 2 machine-learning algorithms and 14 LSMs to analyze the spatial and temporal variations in global terrestrial ET. The results showed that the ensemble means of annual global terrestrial ET estimated by these three categories of approaches agreed well, with values ranging from 589.6 mm/yr (6.56×10^4 cu.km/yr) to 617.1 mm/yr (6.87×10^4 cu.km/yr). For the period from 1982 to 2011, both the ensembles of remote-sensing-based physical models and machine-learning algorithms suggested increasing trends in global terrestrial ET (0.62 mm/sq.yr with a significance level of p<0.05 and 0.38 mm yr−2 with a significance level of p<0.05, respectively). In contrast, the ensemble mean of the LSMs showed no statistically significant change (0.23 mm/sq.yr, p>0.05), although many of the individual LSMs reproduced an increasing trend. Nevertheless, all 20 models used in this study showed that anthropogenic Earth greening had a positive role in increasing terrestrial ET. The concurrent small interannual variability, i.e., relative stability, found in all estimates of global terrestrial ET, suggests that a potential planetary boundary exists in regulating global terrestrial ET, with the value of this boundary being around 600 mm/yr. Uncertainties among approaches were identified in specific regions, particularly in the Amazon Basin and arid/semiarid regions. Improvements in parameterizing water stress and canopy dynamics, the utilization of new available satellite retrievals and deep-learning methods, and model–data fusion will advance our predictive understanding of global terrestrial ET.

surface modeling↗
NASA NTRS2019

State of the science in reconciling top‐down and bottom‐up approaches for terrestrial CO2 budget

Robust estimates of CO2 budget, CO2 exchanged between the atmosphere and terrestrial biosphere, are necessary to better understand the role of the terrestrial biosphere in mitigating anthropogenic CO2 emissions. Over the past decade, this field of research has advanced through understanding of the differences and similarities of two fundamentally different approaches: “top‐down” atmospheric inversions and “bottom‐up” biosphere models. Since the first studies were undertaken, these approaches have shown an increasing level of agreement, but disagreements in some regions still persist, in part because they do not estimate the same quantity of atmosphere–biosphere CO2 exchange. Here, we conducted a thorough comparison of CO2 budgets at multiple scales and from multiple methods to assess the current state of the science in estimating CO2 budgets. Our set of atmospheric inversions and biosphere models, which were adjusted for a consistent flux definition, showed a high level of agreement for global and hemispheric CO2 budgets in the 2000s. Regionally, improved agreement in CO2 budgets was notable for North America and Southeast Asia. However, large gaps between the two methods remained in East Asia and South America. In other regions, Europe, boreal Asia, Africa, South Asia, and Oceania, it was difficult to determine whether those regions act as a net sink or source because of the large spread in estimates from atmospheric inversions. These results highlight two research directions to improve the robustness of CO2 budgets: (a) to increase representation of processes in biosphere models that could contribute to fill the budget gaps, such as forest regrowth and forest degradation; and (b) to reduce sink–source compensation between regions (dipoles) in atmospheric inversion so that their estimates become more comparable. Advancements on both research areas will increase the level of agreement between the top‐down and bottom‐up approaches and yield more robust knowledge of regional CO2 budgets.

atmospheric inversion↗
NASA NTRS2019

Global Carbon Budget 2019

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (E(FF)) are based on energy statistics and cement production data, while emissions from land use change (E(LUC)), mainly deforestation, are based on land use and land use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (G(ATM)) is computed from the annual changes in concentration. The ocean CO2 sink (S(OCEAN)) and terrestrial CO2 sink (S(LAND)) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (B(IM)), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2009–2018), E(FF) was 9.5±0.5 GtC/yr, E(LUC) 1.5±0.7 GtC/yr, G(ATM) 4.9±0.02 GtC/yr (2.3±0.01 ppm/yr), S(OCEAN) 2.5±0.6 GtC/yr, and S(LAND) 3.2±0.6 GtC/yr, with a budget imbalance B(IM) of 0.4 GtC/yr indicating overestimated emissions and/or underestimated sinks. For the year 2018 alone, the growth in E(FF) was about 2.1 % and fossil emissions increased to 10.0±0.5 GtC/yr, reaching 10 GtC/yr for the first time in history, E(LUC) was 1.5±0.7 GtC/yr, for total anthropogenic CO2 emissions of 11.5±0.9 GtC/yr (42.5±3.3 GtCO2). Also for 2018, G(ATM) was 5.1±0.2 GtC/yr (2.4±0.1 ppm/yr), S(OCEAN) was 2.6±0.6 GtC/yr, and S(LAND) was 3.5±0.7 GtC/yr, with a B(IM) of 0.3 GtC. The global atmospheric CO2 concentration reached 407.38±0.1 ppm averaged over 2018. For 2019, preliminary data for the first 6–10 months indicate a reduced growth in E(FF) of +0.6 % (range of −0.2 % to 1.5 %) based on national emissions projections for China, the USA, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. Overall, the mean and trend in the five components of the global carbon budget are consistently estimated over the period 1959–2018, but discrepancies of up to 1 GtC/yr persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations shows (1) no consensus in the mean and trend in land use change emissions over the last decade, (2) a persistent low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018a, b, 2016, 2015a, b, 2014, 2013).

global carbon budget 2019↗
NASA NTRS2019

Negative extreme events in gross primary productivity and their drivers in China during the past three decades

Climate extremes have remarkable impacts on ecosystems and are expected to increase with future global warming. However, only few studies have focused on the ecological extreme events and their drivers in China. In this study, we carried out an analysis of negative extreme events in gross primary productivity (GPP) in China and the sub-regions during 1982–2015, using monthly GPP simulated by 12 process-based models (TRENDYv6) and an observation-based model (Yao-GPP). Extremes were defined as the negative 5th percentile of GPP anomalies, which were further merged into individual extreme events using a three-dimensional contiguous algorithm. Spatio-temporal patterns of negative GPP anomalies were analyzed by taking the 1000 largest extreme events into consideration. Results showed that the effects of extreme events decreased annual GPP by 2.8% (i.e. 208 TgC/year) in TRENDY models and 2.3% (i.e. 151 TgC/year) in Yao-GPP. Hotspots of extreme GPP deficits were mainly observed in North China (−53 gC/sq.m year) in TRENDY models and Northeast China (−42 gC/sq.m year) in Yao-GPP. For China as a whole, attribution analyses suggested that extreme low precipitation was associated with 40%–50% of extreme negative GPP events. Most events in northern and western China could be explained by meteorological droughts (i.e. low precipitation) while GPP extreme events in southern China were more associated with temperature extremes, in particular with cold spells. GPP was revealed to be much more sensitive to heat/drought than to cold/wet extreme events. Combined with projected changes in climate extremes in China, GPP negative anomalies caused by drought events in northern China and by temperature extremes in southern China might be more prominent in the future.

Climate change↗