Abstract
The terrestrial ecosystem in East Asia mainly consists of semi-arid regions that are sensitive to climate change. Therefore, gross primary productivity (GPP) in East Asia could be highly variable and vulnerable to climate change, which can significantly affect the local carbon budget. Here, we examine the spatial and temporal characteristics of GPP variability in East Asia and its relationship with climate factors over the last three decades. We detect an abrupt decrease in GPP over Eastern China-Mongolia region around the year 2000. This is attributed to an abrupt decrease in precipitation associated with the phase shift of the Pacific decadal oscillation (PDO). We also evaluate the reproducibility of offline land surface models to simulate these abrupt changes. Of the twelve models, eight were able to simulate this abrupt response, while the others failed due to the combination of an exaggerated CO2 fertilization effect and an underrated climate impact. For accurate prediction, it is necessary to improve the sensitivity of the GPP to changes in CO2 concentrations and the climate system.
Similar content being viewed by others
Introduction
Land carbon sink is a major driver of the interannual variability of atmospheric CO2 concentrations1,2,3. Gross primary productivity (GPP), the organic carbon sequestrated by terrestrial plants through photosynthesis4, has a significant impact on the variability of the land carbon sink at regional and global scales5,6. GPP is largely modulated by global environmental changes, such as CO2 concentrations and meteorological factors5. For instance, GPP increases in response to rising atmospheric CO2 concentrations, which is referred to as the CO2 fertilization effect7,8. In addition, the productivity of terrestrial ecosystems significantly depends on climate conditions such as precipitation and temperature3,9,10. Thus, the amount of land carbon sequestration varies with these factors and could exhibit strong interannual variability11. Moreover, strong meteorological forcing can cause abrupt changes in GPP, which are non-linear and non-stationary12. For example, lower absolute soil moisture leads to a prominent abrupt decrease in GPP in arid regions of Europe13. Specifically, these abrupt decreases in terrestrial carbon fluxes can lead to a deterioration in local carbon uptake.
Previous studies have reported abrupt changes in vegetation productivity and the underlying mechanisms in semi-arid ecosystems. Ma et al.14 have examined how climate extremes impact semi-arid ecosystems in southeastern Australia, with abrupt shifts in phenology and vegetation productivity in drought years. Berdugo et al.15 has shown the ubiquity of abrupt productivity loss in global drylands due to climatic, edaphic, and human factors from 2000 to 2019. They show that abrupt decreases in Normalized Difference Vegetation Index (NDVI) have been detected mainly in semi-arid forest regions15. As semi-arid ecosystems are highly sensitive to water availability15,16, drought or negative trend in precipitation could lead to abrupt decreases in vegetation productivity or loss of ecosystem resilience14,17.
East Asia is known to be a significant contributor to the global carbon cycle7,18,19. Fossil fuel CO2 emissions in East Asia are about 1.5 PgC year−1, and about 13–27% of these emissions have been offset by its terrestrial ecosystem from 1990 to 20097. Based on the prevailing semi-arid regions and historical extreme climate events in Inner East Asia20, it is possible that there have been abrupt changes in GPP over East Asia20. However, there is a limited number of studies on the abrupt changes in GPP, particularly the decrease in this region. Although several studies have detected the change points of GPP in East Asia4,5,16, they have mostly focused on greening and monotonic increasing trend of GPP. Furthermore, it is unclear whether current land models can reproduce or simulate the abrupt change in this region.
The primary objective of this study is to identify the major variability and occurrence of abrupt changes in GPP in East Asia over the last three decades. In detail, we examine the dominant spatio-temporal characteristics of GPP variability, detect the significant abrupt changes in GPP, and identify their main drivers. We utilize multiple GPP datasets to validate the consistency of the results. Furthermore, we assess the capabilities of land models to simulate the variability and abrupt changes in GPP.
Results
The abrupt shift in GPP over Eastern China-Mongolia
To investigate the dominant spatio-temporal characteristics of GPP variability, we examine the first leading mode of GPP in East Asia using the Empirical Orthogonal Function (EOF) (see “Data and methods” section, Fig. 1). The first eigenvectors and principal components (PCs) of GPPFLUXCOM and GPPNIRv explain 21.4% and 16.4% of the total variance, respectively. The maximum variation of GPPFLUXCOM and GPPNIRv is consistently located in Eastern China-Mongolia region (40°–52°N and 110°–124°E). Furthermore, the PCs of both first leading modes consistently show a distinctive phase shift around the year 2000. Specifically, there is a positive phase until the 2000s and afterward a negative phase, suggesting a decrease in vegetation productivity around the year 2000 in Eastern China-Mongolia region. In this region, the first eigenvectors and PCs of GPPFLUXCOM and GPPNIRv explain 41.2% and 31.2% of the total variance, respectively (Supplementary Fig. 1). Note that the phase shift in GPP with decadal time scale is comparable in magnitude to their interannual variability, and even larger for GPPNIRv in particular.
(a, b) Eigenvectors of the first leading mode and (c, d) its principal components (PCs) time series from the Empirical Orthogonal Function (EOF) analysis of JJA (June–July–August) GPPFLUXCOM (1980–2016) and GPPNIRv (1982–2018) over East Asia (24°–52°N and 100°–149°E). The values of explained variance for GPPFLUXCOM and GPPNIRv are 21.4% and 16.4% of the total variance, respectively. The red box in a and b indicates Eastern China-Mongolia region (40°–52°N and 110°–124°E).
Eastern China-Mongolia is a semi-arid temperate grassland region with limited precipitation21,22. Previous studies have reported that vegetation productivity is strongly dependent on water availability in semi-arid regions across Northern China16,22, suggesting the possibility that changes in GPP are caused by changes in precipitation over Eastern China-Mongolia region. Therefore, we further examined the changes in GPP and precipitation in this region and their relationship (Fig. 2). There is a prominent decreasing shift in GPP in the late 1990s; a positive anomaly before the late 1990s and a negative anomaly afterward. This result is consistent with the PCs of the first leading mode shown in Fig. 1c and d. The pronounced decreasing shift in precipitation in the late 1990s is also observed. This is consistent with previous studies that found a decreasing shift in summer precipitation in Northeastern and North China in the late 1990s23,24. From 1999 to 2007, the average GPPNIRv (GPPFLUXCOM) decreased by 11% (6%) compared to the period of 1990 to 1998, and precipitation also decreased by 28% during the same period. The decreases of GPPNIRv and GPPFLUXCOM in 1998 and 2001 are 36.3 and 12.0 TgC month−1, respectively. These changes are about four times larger than the standard deviation of GPPNIRv (1982–2018) and GPPFLUXCOM (1980–2016).
(a) Time series of JJA mean GPPFLUXCOM anomaly (1980–2016, black solid line), CRUNCEP precipitation anomaly (1980–2016, blue dashed line), and the Pacific decadal oscillation (PDO) index (1980–2018, gray bar). (b) Same as (a), but for GPPNIRv anomaly (1982–2018, black solid line). GPP and precipitation anomalies are averaged over Eastern China-Mongolia region (40°–52°N and 110°–124°E) as shown in the red box in Fig. 1.
The GPPFLUXCOM and GPPNIRv anomalies are positively correlated with the precipitation anomaly in Eastern China-Mongolia (r = 0.87 and 0.63, P < 0.01). These findings exhibit similarity across another precipitation dataset (Supplementary Fig. 2). These results are consistent with a previous study that showed a positive relationship between decreasing precipitation and GPP over Northern China from 1999 to 201116. On the contrary, temperature and solar radiation have relatively low correlation coefficients (r = − 0.54, P < 0.01 and r = − 0.38, P < 0.05) with GPPFLUXCOM anomaly and no statistically significant correlation with GPPNIRv anomaly at the 95% confidence level. These results suggest that the decreasing shift in precipitation is probably responsible for that of GPP in the late 1990s. Moreover, considering the positive correlation between precipitation and evapotranspiration (r = 0.68, P < 0.01) observed in Eastern China-Mongolia, it is conceivable that other processes may contribute to the dramatic shift. For example, changes in the strength of the land–atmosphere coupling can lead to rapid changes in GPP25, which will be discussed in more detail in the “Discussion” section.
We further examine whether these decreasing shifts of GPP and precipitation are statistically significant based on the Lepage test (see “Data and methods” section). The phase shifts of GPPFLUXCOM, GPPNIRv, and precipitation around 1999 are all statistically significant (Fig. 3). GPPFLUXCOM and GPPNIRv show the most significant abrupt shifts in 1999–2000 and 2000–2001, respectively (P < 0.05). Precipitation also shows the most significant abrupt shift in 1999–2000 (P < 0.05). This result exhibits similarity across another precipitation dataset. (Supplementary Fig. 3). This is consistent with the results of previous studies23,24, which showed the abrupt change point of summer precipitation around 1999 over Northeastern China. These highest HK values are largely determined by the standardized Wilcoxon rank sum statistic, indicating that these statistically significant abrupt shifts mostly come from the large changes in the mean state. The results of the Lepage test, conducted with different lengths of the moving window, indicate the consistent abrupt shift year for GPP and precipitation (Supplementary Fig. 4).
Time series of the Lepage statistic (HK) values of JJA mean GPPFLUXCOM (black solid line), GPPNIRv (black dashed line), and CRUNCEP precipitation (blue solid line) for a window length of 9 years. If the HK value is higher than 5.99, the difference between the means of the two samples is significant at the 95% confidence level. GPP and precipitation anomalies are averaged over Eastern China-Mongolia region (40°–52°N and 110°–124°E) as shown in the red box in Fig. 1.
Vegetation productivity is significantly regulated by temperature and CO2 concentrations, as well as precipitation10. For example, the abrupt changes in NDVI in Southwest China were mainly driven by changes in temperature during 1982–201512. Therefore, to confirm the quantitative contributions of climate factors and CO2 to the regime shift in GPP around 1999, we evaluate the relative contributions of changes in temperature, precipitation, and CO2 concentrations to the GPP shift using the multiple linear regression method (see “Data and methods” section). With this method, we reconstructed GPP using climate factors and CO2, estimating the contributions of the climate and CO2 factors to the GPP shift between period 1 (P1: 1990 to 1998) and period 2 (P2: 1999 to 2007). It is evident that there is a decrease in the reconstructed GPPFLUXCOM (− 1.42) and GPPNIRv (− 0.93) between P1 and P2, indicating a significant decrease in GPP during P2 (Fig. 4). It is clear that the decrease in GPP is mostly due to the change in precipitation (GPPFLUXCOM: − 1.31, GPPNIRv: − 1.19). Temperature also contributes to the decrease, but its impact is minor (GPPFLUXCOM: − 0.29, GPPNIRv: − 0.064). Increasing CO2 leads to an increase in GPP in P2 due to the fertilization effect7,26. These results indicate that the change in precipitation is mainly responsible for the statistically significant regime shift in GPP around 1999. In summary, there is a significant abrupt change in GPP in Eastern China-Mongolia over the past three decades, which can be attributed to the abrupt decrease in precipitation.
Differences in JJA mean of normalized GPPFLUXCOM and GPPNIRv between P1 (1990 to 1998) and P2 (1999 to 2007), reconstructed from multiple linear regression based on normalized climate factors (precipitation and temperature) and CO2. Each bar is the quantitative contribution of precipitation (blue), temperature (red), and CO2 (green) to the total difference in GPP (gray). A negative (positive) value indicates a lower (higher) mean state of GPP during P2 compared to P1. All variables are averaged over Eastern China-Mongolia region (40°–52°N and 110°–124°E) as shown in the red box in Fig. 1.
The inter-model diversity of the abrupt shift in GPP over Eastern China-Mongolia
To support the robustness of the regime shift in GPP, we further analyze the output of offline land surface models participating in TRENDY. Using the simulations provided by TRENDY, we isolate the time-varying GPPclimate-forcing and the GPPCO2-forcing (see “Data and methods” section). As shown in Fig. 5a, an abrupt shift in GPP around the year 2000 is well captured in the multi-model ensemble (MME) mean of the TRENDY models. This change is statistically significant based on the Lepage test (P < 0.05) (Fig. 5b). However, there is significant diversity and spread among the individual TRENDY models, especially in their HK values (Fig. 5b). Of the twelve models, only eight models show a significant abrupt shift in GPP in 2000–2001 (Strong shift group: CLASS-CTEM, CLM5.0, JULES-ES, JSBACH, ORCHIDEE-CNP, ORCHIDEE, SDGVM, and VISIT), while the other models do not simulate abrupt shifts (Weak shift group: LPJ-GUESS, CABLE-POP, DLEM, and ISBA-CTRIP).
(a) Time series of JJA mean GPP anomalies: multi-model ensemble (MME) mean of GPPTRENDY (black dashed line), the mean of the strong shift group (red solid line), and the weak shift group (blue solid line). The red (blue) horizontal bars show the mean values of GPP during P1′ (1991 to 1999) and P2′ (2000 to 2008) of the strong shift group (the weak shift group). (Strong shift group: CLASS-CTEM, CLM5.0, JULES-ES, JSBACH, ORCHIDEE-CNP, ORCHIDEE, SDGVM, and VISIT; Weak shift group: LPJ-GUESS, CABLE-POP, DLEM, and ISBA-CTRIP). (b) Time series of Lepage statistics (HK) values of MME mean of JJA GPPTRENDY (black dashed line) and individual TRENDY models (colored solid lines) for a window length of 9 years. If the HK value is higher than 5.99, the difference between the means of the two samples is significant at the 95% confidence level. GPP anomalies are averaged over Eastern China-Mongolia region (40°–52°N and 110°–124°E) as shown in the red box in Fig. 1.
We further investigate the differences in GPP between the strong and weak shift groups (Fig. 5a). The red and blue horizontal bars show the mean values of GPP for the two groups during period 1 (P1′: 1991 to 1999) and period 2 (P2′: 2000 to 2008). Although the two groups show quite similar interannual variability, their decadal variability is different. The weak shift group shows a smaller mean state change (− 5.64 TgC month−1) in GPP between the two periods than the strong shift group (− 17.4 TgC month−1).
To examine the driver of the abrupt shift in GPP and the cause of the differences between the groups, we quantify the contributions of the climate and CO2 forcings to the mean state change in GPP from the TRENDY models. We calculated the differences in normalized GPP separated by each forcing by subtracting the mean of P1′ from P2′ (Fig. 6). In all models, GPP decreases from P1′ to P2′ and is mostly driven by the GPPclimate-forcing (Fig. 6a). The decrease in GPPclimate-forcing is mostly driven by changes in precipitation rather than temperature and solar radiation (Supplementary Fig. 5). Therefore, in terms of GPPclimate-forcing, the TRENDY models also present the aforementioned precipitation-induced regime shift in GPP, which is consistent with satellite and reanalysis data.
(a) Differences in JJA mean of normalized GPP between P1′ and P2′ for the individual TRENDY models. Each bar is the quantitative contribution of climate (red) and CO2 forcing (blue) to the total difference in normalized GPP (gray). The models are arranged in order of the magnitude of the HK values. (b) Difference in JJA mean of normalized GPP between P1′ and P2′ from the strong and weak shift groups. A negative (positive) value indicates a lower (higher) mean state of GPP during P2′ compared to P1′. GPP anomalies are averaged over Eastern China-Mongolia region (40°–52°N and 110°–124°E) as shown in the red box in Fig. 1.
The effect of CO2 forcing is opposite to that of climate forcing, but the impact is relatively small: the MME mean of GPPCO2-forcing differences (0.26) is smaller in magnitude compared to that of GPPclimate-forcing (− 1.40). In particular, the four models in the weak shift group (LPJ-GUESS, CABLE-POP, DLEM, and ISBA-CTRIP) have a higher sensitivity of GPP to CO2 and a lower sensitivity to climate compared to the other group (Fig. 6a). These characteristics are more obvious in the group means (Fig. 6b): the MME means of the GPPCO2-forcing and GPPclimate-forcing differences for the weak shift group (the strong shift group) are 0.42 (0.18) and − 0.98 (− 1.60), respectively. Therefore, the differences in the sensitivity of GPP to CO2 and climate forcing between the models lead to significant differences in the changes in GPP between the two groups.
Discussion
We have investigated the first EOF mode of GPP variability and found significant abrupt changes in GPP in East Asia over the past three decades. The PCs of the first leading modes of GPPFLUXCOM and GPPNIRv consistently show a distinctive phase shift around the year 2000 in Eastern China-Mongolia. The abrupt change in GPP around the year 2000 is consistently detected by satellite, reanalysis datasets, and land surface models. The high correlation between GPP and precipitation, and the abrupt decrease in precipitation indicate that the abrupt decrease in GPP resulted from the changes in precipitation. Quantitative analysis using the land surface models shows that the shift in GPP is mostly driven by changes in precipitation, supporting the GPP-precipitation mechanism.
It has been reported that the Pacific decadal oscillation (PDO) phase change can affect the regional precipitation in East Asia through atmospheric teleconnections24,28,29,30,31,32. After the late 1990s, the Pacific entered a negative phase of the PDO (Fig. 2), and the upper-level atmospheric circulation changed accordingly. The East Asian westerly jet stream (EAWJS) is weakened and poleward shifted. The poleward-shifted EAWJS changes the jet-related secondary meridional-vertical circulation29, making an anomalous descending motion to the northern parts of East Asia. Subsequently, the anticyclonic circulation dominates the target region in a negative PDO phase, resulting in less precipitation in the late 1990s. Therefore, the precipitation shifts around the year 2000 are related to the phase shift of the PDO. Our results suggest that the abrupt changes in GPP in Eastern China-Mongolia region depend on the phase change of the PDO.
Of the twelve offline land surface models, eight were able to simulate this abrupt response, while the others failed to capture it. The differences between the TRENDY models can be attributed to the diverse sensitivities of CO2 and climate forcing. This implies that differences in the sensitivity of GPP to CO2 and climate forcing between the models resulted in the inter-model diversity in simulating the abrupt changes. One factor contributing to the differences in sensitivity is that terrestrial photosynthetic CO2 assimilation, which is the basis for model estimation of GPP, is represented and parameterized differently across terrestrial biosphere models33. This leads to different responses of GPP to key environmental variables, including climate and CO2 forcing33. Therefore, this leads to structural differences (e.g., sensitivities) and large uncertainties in the models34. Therefore, improving the response and parameterization of photosynthesis in each model would be helpful to reduce the uncertainties in the TRENDY models and investigate the more accurate variability of GPP in Eastern China-Mongolia.
It is possible that land–atmosphere coupling, a positive feedback between vegetation and precipitation recycling25, is involved in the abrupt change in GPP over Eastern China-Mongolia. In this feedback loop, reduced precipitation in water-limited areas leads to a decline in vegetation, which can decrease evapotranspiration. A decrease in precipitation due to the reduced evapotranspiration is expected to lead to a decrease in vegetation. The Eastern China-Mongolia region is known to have a high precipitation recycling ratio35; in other words, precipitation is highly contributed by evapotranspiration from the same region. Thus, precipitation and evapotranspiration have a positive correlation (r = 0.68, P < 0.01) from 1990 to 2007. Accordingly, evapotranspiration also shows a significant abrupt shift in 1999–2000 in the Eastern China-Mongolia region (P < 0.05). Given these, it is conceivable that the positive feedback could contribute to the dramatic shift in GPP and its persistence. Due to the significant effects of the feedback, further studies are needed to understand the role of land–atmosphere coupling in abrupt changes in vegetation productivity.
Semiarid carbon flux is known to be influenced not only by precipitation, but also by temperature, as warming-induced water deficit also leads to a decrease in photosynthesis1. In this study, compared to precipitation, temperature exhibits a relatively low correlation coefficient (r = − 0.54, P < 0.01) with GPPFLUXCOM anomaly in 1980–2016 and no statistically significant correlation with GPPNIRv anomaly at the 95% confidence level. However, over a shorter time period (1998–2016), GPPFLUXCOM and GPPNIRv anomalies have a stronger negative correlation with temperature anomaly (r = − 0.68, P < 0.01 and − 0.56, P < 0.05). However, this relationship was not considered in this study. Indeed, GPP showed a delayed response to the recovery of precipitation in 2003–2004 (Fig. 2). Based on the negative correlation between GPP and temperature, it can be suggested that the positive temperature anomaly in the early 2000s may have contributed to the delay in GPP recovery in 2003–2004 (figure not shown). Therefore, it is necessary to consider temperature to understand the overall variability of GPP in East Asia comprehensively.
The GPP time span (Fluxcom, 1950–2016; NIRv, 1982–2018; Trendy, 1700–2018) does not account for recent changes in GPP. A recent abrupt shift towards a hotter and drier climate in inner East Asia20 could indicate the possibility of frequent abrupt changes in GPP. Thus, it is necessary to conduct an extended study to investigate the impacts of a hotter and drier climate on GPP in inner East Asia. In addition, future abrupt changes in vegetation productivity in East Asia need to be investigated. Drought risks in East Asia (30°–55°N and 110°–140°E) would increase due to rapidly rising evaporative demand27. Moreover, the projections based on the representative concentration pathways (RCPs) RCP8.5 and RCP4.5 indicate that East Asia (30°–70°N and 80°–180°E) is one of the regions that will experience the most significant expansion of drylands by the end of the twenty-first century36,37. The high correlation between water availability and GPP suggests that abrupt changes in vegetation productivity may become more widespread and frequent in East Asia in the future. Furthermore, the drylands of East Asia will play an important role in regulating the interannual variability of global dryland GPP37. This suggests that the intense occurrence of abrupt changes could exacerbate the interannual variation of the global land sink in the future. In this regard, it is essential to conduct further research on future abrupt changes in vegetation productivity and terrestrial carbon flux. Further research should focus on the northeastern part of East Asia, where dryland expansion is expected27, and the western part of East Asia, where future drought risks are expected to be higher37. This further research could contribute to developing carbon–neutral strategies by enhancing the understanding of the interannual variability of the land sink.
Data and methods
GPP datasets
Terrestrial photosynthesis is a highly uncertain process due to the lack of direct observations of GPP on a global scale. To obtain robust results on GPP responses to climate variability, we used both data-driven and process-based GPP datasets with over 30 years of records for the East Asian region.
FLUXCOM
FLUXCOM (version RS + METEO) provides ensembles of upscaled GPP products on a global scale (www.fluxcom.org) for the period from 1950 to 201638,39. FLUXCOM GPP was developed using machine learning techniques to merge FLUXNET site-level observations, satellite remote sensing, and meteorological data to scale these fluxes to the global scale34,40,41. This study used the FLUXCOM GPP, which is calculated by averaging GPP ensembles generated from three machine learning algorithms (Random Forest, Artificial Neural Network, and Multivariate Adaptive Regression Splines) and two flux partitioning methods (daytime and nighttime) for the period from 1980 to 201642.
Near-infrared reflectance (NIRv)
Previous studies have shown that both NIRv and SIF effectively capture global-scale changes in terrestrial GPP43,44,45. This study used NIRv due to its broader data coverage and higher resolution compared to satellite-based SIF46 to investigate long-term trends in GPP. The NIRv data are derived from the long-term satellite datasets of the Advanced Very High Resolution Radiometer (AVHRR) reflectance from the Land Long Term Data Record v4 (LTDR) product for the period from 1982 to 2018. This study utilizes the dataset for all available time periods. NIRv is calculated as a function of monthly NDVI and near-infrared reflection of the total pixel (NIRT)43,44. A subtraction of 0.08 from the NDVI is applied to eliminate the influence of bare soils46 [Eq. (1)].
TRENDY v8
The model intercomparison project TRENDY v8 provides GPP output from Dynamic Global Vegetation Models (DGVMs) for the period from 1700 to 201847. These models are forced with the same meteorological conditions, atmospheric CO2 concentrations, and land use datasets34,46,48. TRENDY v8 provides a suite of simulations that allow the response of land to climate, CO2, and land use forcing to be separated. The “S0” simulation is a baseline simulation with time-invariant pre-industrial CO2 concentrations, climate, and land use. The “S1” simulation is forced by time-varying atmospheric CO2 concentrations, but with fixed climate and land use information. The “S2” simulation is forced by time-varying atmospheric CO2 concentrations and climate with constant land use change.
This study used the GPP output from 12 DGVMs from 1980 to 2018, which was regridded to a spatial resolution of 1° × 1° for consistent analysis in the multi-model domain (Table 1). Based on the experimental designs, we utilize the experimental simulations to evaluate the relative contributions of CO2 fertilization and climate forcing to the change in GPP. For each TRENDY model, the CO2 fertilization effect-driven change in GPP (GPPCO2-forcing) is calculated from the difference between S1 and S0, while the climate change-driven change in GPP (GPPclimate-forcing) is calculated from the difference between S2 and S1. In addition, the response of GPP to total forcing (GPPTRENDY) is calculated as the difference between S2 and S0, which is the sum of CO2 fertilization and climate change-driven GPP. Note that this study excludes the effect of land use change on GPP, as land use change does not appear to be a significant driver of GPP change over the study period (Supplementary Fig. 6).
We focused on averaged GPP in the boreal summer months of June–July–August (JJA) to examine spatio-temporal variations across the East Asia region (24°–52°N and 100°–149°E). The JJA means of GPPFLUXCOM and GPPNIRv have the highest values and variability, accounting for 78.3% and 78.7% of annual vegetation productivity, respectively (Supplementary Fig. 7).
Meteorological data
We used precipitation, temperature, and solar radiation data from CRUNCEPv7 for the period 1980 to 201649, along with evapotranspiration data from ERA5 for the period 1980 to 201851. We also utilized precipitation data from GPCC for the period 1982 to 2018 to verify the consistency of the results52. Atmospheric CO2 concentration data from the NOAA GML Carbon Cycle Cooperative Global Air Sampling Network were also used for the period 1980 to 201650. Table 2 summarizes the meteorological data used in this study.
EOF-based detection of dominant spatial and temporal variability
The EOF analysis is considered to be an effective approach to investigate the temporal and spatial coherence for each orthogonal component of complex time-varying spatial patterns53. It performs a linear transformation, defined in terms of the eigenvectors of their covariance matrix, which provides dimension reduction54. This analysis allows us to examine the structural features of a long-time series variable and extract its representation of variability. Therefore, it effectively captures representations of climate variability and identifies the impacts of large-scale climate change and abnormal climate on ecosystems55. Thus, the EOF analysis has been used to characterize the dominant modes of variability in GPP and vegetation indexes in previous studies53,56,57.
Lepage test for abrupt change detection
There are several methods for detecting change points in adjacent data on a decadal scale. Among these methods, the Lepage test58 has few underlying assumptions, which makes it adopted for detecting various types of climate changes, including linear trends, cyclical variations, step-like changes, and discontinuous changes23. This method has been used extensively in numerous studies to detect abrupt changes23,24,59,60. The Lepage statistic (HK) is calculated as the sum of the squares of the standardized Wilcoxon (W) and Ansari-Bradley (A) statistics [Eq. (2)]61, which represent the comparison of the mean and variance between two samples, respectively58.
where E(W), E(A), V(W), and V(A) stand for the expected value and variance of W and A, respectively. If the HK is greater than 5.99, the difference between the two samples is significant at the 95% confidence level. A detailed description of the methods can be found in an earlier study23. In this study, the Lepage test with a 9-year moving window is used to identify the inter-decadal abrupt change point. We defined the year with the largest HK value above 5.99 as the year in which an abrupt change occurred, based on previous studies24,62.
Quantifying the contribution of climate and CO2 to the change in GPP
The multiple linear regression method was used to analyze the linear relationship between multiple variables63. This method was used to investigate the relative contributions of climatic (temperature and precipitation) and CO2 factors to vegetation productivity [Eq. (3)].
where Prec, Temp, and CO2 stand for precipitation, temperature, and CO2 concentration, respectively. \(\alpha\), \(\upbeta\), and \(\gamma\) represent the partial regression coefficients for each variable. \(\delta\) represents the residuals. We used normalized GPP, precipitation, temperature, and CO2 concentrations.
Data availability
All data used in this study are publicly available and can be downloaded from the corresponding websites (FLUXCOM: https://fluxcom.org/; NIRv: https://data.tpdc.ac.cn/; Access to TRENDY v.8 data can be obtained by contacting Stephen Sitch: (S.A.Sitch@exeter.ac.uk); CRUNCEPv7 datasets (precipitation, temperature, and solar radiation): https://rda.ucar.edu/datasets/ds314.3/citation/; ERA5 datasets (evapotranspiration): https://cds.climate.copernicus.eu/; GPCC (precipitation): https://psl.noaa.gov/; CO2 concentration: https://gml.noaa.gov/).
Code availability
The computer codes that support the analysis within this paper are available from the corresponding author on request.
References
Zhao, H., Jia, G., Wang, H., Zhang, A. & Xu, X. Seasonal and interannual variations in carbon fluxes in East Asia semi-arid grasslands. Sci. Total Environ. 668, 1128–1138 (2019).
Kim, J.-S., Kug, J.-S., Yoon, J.-H. & Jeong, S.-J. Increased atmospheric CO2 growth rate during El Niño driven by reduced terrestrial productivity in the CMIP5 ESMs. J. Clim. 29, 8783–8805 (2016).
Park, S.-W. et al. Two aspects of decadal ENSO variability modulating the long-term global carbon cycle. Geophys. Res. Lett. 47, e086390 (2020).
He, P. et al. Spatiotemporal evolutionary and mechanism analysis of grassland GPP in China. Ecol. Indic. 143, 109323 (2022).
Ma, J. et al. Trends and controls of terrestrial gross primary productivity of China during 2000–2016. Environ. Res. Lett. 14, 084032 (2019).
Zhang, Y., Dannenberg, M. P., Hwang, T. & Song, C. E. Niño-Southern Oscillation-Induced Variability of Terrestrial Gross Primary Production During the Satellite era. J. Geophys. Res. Biogeosci. 124, 2419–2431 (2019).
Piao, S. L. et al. The carbon budget of terrestrial ecosystems in East Asia over the last two decades. Biogeosciences 9, 3571–3586 (2012).
Park, S.-W., Kug, J.-S., Jun, S.-Y., Jeong, S.-J. & Kim, J.-S. Role of cloud feedback in continental warming response to CO2 physiological forcing. J. Clim. 34, 8813–8828 (2021).
Sakalli, A., Cescatti, A., Dosio, A. & Gücel, M. U. Impacts of 2°C global warming on primary production and soil carbon storage capacity at pan-European level. Clim. Serv. 7, 64–77 (2017).
Liu, Z. et al. Global divergent responses of primary productivity to water, energy, and CO2. Environ. Res. Lett. 14, 124044 (2019).
Friedlingstein, P. et al. Global carbon budget 2022. Earth Syst. Sci. Data 14, 4811–4900 (2022).
Xu, X., Liu, H., Lin, Z., Jiao, F. & Gong, H. Relationship of abrupt vegetation change to climate change and ecological engineering with multi-timescale analysis in the Karst region, Southwest China. Remote Sens. 11, 1564 (2019).
Sungmin, O. & Park, S. K. Flash drought drives rapid vegetation stress in arid regions in Europe. Environ. Res. Lett. 18, 014028 (2023).
Ma, X., Huete, A., Moran, S., Ponce-Campos, G. & Eamus, D. Abrupt shifts in phenology and vegetation productivity under climate extremes. J. Geophys. Res. Biogeosci. 120, 2036–2052 (2015).
Berdugo, M., Gaitán, J. J., Delgado-Baquerizo, M., Crowther, T. W. & Dakos, V. Prevalence and drivers of abrupt vegetation shifts in global drylands. Proc. Natl. Acad. Sci. 119, e2123393119 (2022).
Yuan, W. et al. Multiyear precipitation reduction strongly decreases carbon uptake over northern China. J. Geophys. Res. Biogeosci. 119, 881–896 (2014).
Poulter, B. et al. Contribution of semi-arid ecosystems to interannual variability of the global carbon cycle. Nature 509, 600–603 (2014).
Piao, S. et al. Contribution of climate change and rising CO2 to terrestrial carbon balance in East Asia: A multi-model analysis. Glob. Planet. Change 75, 133–142 (2011).
Kim, J.-S. et al. Arctic warming-induced cold damage to East Asian terrestrial ecosystems. Commun. Earth Environ. 3, 16 (2022).
Zhang, P. et al. Abrupt shift to hotter and drier climate over inner East Asia beyond the tipping point. Science 370, 1095–1099 (2020).
Kang, X. et al. Modeling impacts of climate change on carbon dynamics in a steppe ecosystem in Inner Mongolia. China. J. Soils Sedim. 11, 562–576 (2011).
Xu, H.-J., Wang, X.-P., Zhao, C.-Y. & Yang, X.-M. Diverse responses of vegetation growth to meteorological drought across climate zones and land biomes in northern China from 1981 to 2014. Agric. For. Meteorol. 262, 1–13 (2018).
Liu, Y., Huang, G. & Huang, R. Inter-decadal variability of summer rainfall in Eastern China detected by the Lepage test. Theor. Appl. Climatol. 106, 481–488 (2011).
Xu, Z., Fan, K. & Wang, H. Decadal variation of summer precipitation over China and associated atmospheric circulation after the late 1990s. J. Clim. 28, 4086–4106 (2015).
D’Odorico, P., Bhattachan, A., Davis, K. F., Ravi, S. & Runyan, C. W. Global desertification: Drivers and feedbacks. Adv. Water Resourc. 51, 326–344 (2013).
Park, S.-W., Kim, J.-S. & Kug, J.-S. The intensification of Arctic warming as a result of CO2 physiological forcing. Nat. Commun. 11, 2098 (2020).
Ha, K.-J., Moon, S., Timmermann, A. & Kim, D. Future changes of summer monsoon characteristics and evaporative demand over Asia in CMIP6 simulations. Geophys. Res. Lett. 47, e2020087492 (2020).
Qian, C. & Zhou, T. Multidecadal variability of North China aridity and its relationship to PDO during 1900–2010. J. Clim. 27, 1210–1222 (2014).
Zhu, Y., Wang, H., Ma, J., Wang, T. & Sun, J. Contribution of the phase transition of Pacific decadal oscillation to the late 1990s’ shift in East China summer rainfall. J. Geophys. Res. Atmos. 120, 8817–8827 (2015).
Piao, J., Chen, W. & Chen, S. Water vapour transport changes associated with the interdecadal decrease in the summer rainfall over Northeast Asia around the late-1990s. Int. J. Climatol. 41, E1469–E1482 (2021).
Dong, X. Influences of the Pacific decadal oscillation on the east Asian summer monsoon in non-ENSO years. Atmos. Sci. Lett. 17, 115–120 (2016).
Jiang, J., Zhou, T., Chen, X. & Wu, B. Central Asian precipitation shaped by the tropical pacific decadal variability and the atlantic multidecadal variability. J. Clim. 34, 7541–7553 (2021).
Rogers, A. et al. A roadmap for improving the representation of photosynthesis in Earth system models. New Phytol. 213, 22–42 (2017).
O’Sullivan, M. et al. Process-oriented analysis of dominant sources of uncertainty in the land carbon sink. Nat. Commun. 13, 4781 (2022).
van der Ent, R. J., Wang-Erlandsson, L., Keys, P. W. & Savenije, H. H. G. Contrasting roles of interception and transpiration in the hydrological cycle. Part 2: Moisture recycling. Earth Syst. Dynam. 5, 471–489 (2014).
Huang, J., Yu, H., Guan, X., Wang, G. & Guo, R. Accelerated dryland expansion under climate change. Nat. Clim. Change 6, 166–171 (2016).
Yao, J. et al. Accelerated dryland expansion regulates future variability in dryland gross primary production. Nat. Commun. 11, 1665 (2020).
Tramontana, G. et al. Predicting carbon dioxide and energy fluxes across global FLUXNET sites with regression algorithms. Biogeosciences 13, 4291–4313 (2016).
Jung, M. et al. Compensatory water effects link yearly global land CO2 sink changes to temperature. Nature 541, 516–520 (2017).
Jung, M. et al. The FLUXCOM ensemble of global land-atmosphere energy fluxes. Sci. Data 6, 74 (2019).
Deng, M. et al. The response of vegetation to regional climate change on the Tibetan plateau based on remote sensing products and the dynamic global vegetation model. Remote Sens. 14, 3337 (2022).
Jung, M. et al. Scaling carbon fluxes from eddy covariance sites to globe: Synthesis and evaluation of the FLUXCOM approach. Biogeosciences 17, 1343–1365 (2020).
Badgley, G., Field, C. B. & Berry, J. A. Canopy near-infrared reflectance and terrestrial photosynthesis. Sci. Adv. 3, e1602244 (2017).
Yang, R. et al. Divergent historical GPP trends among state-of-the-art multi-model simulations and satellite-based products. Earth Syst. Dyn. 13, 833–849 (2022).
Guanter, L. et al. Global and time-resolved monitoring of crop photosynthesis with chlorophyll fluorescence. Proc. Natl. Acad. Sci. 111, E1327–E1333 (2014).
Wang, S., Zhang, Y., Ju, W., Qiu, B. & Zhang, Z. Tracking the seasonal and inter-annual variations of global gross primary production during last four decades using satellite near-infrared reflectance data. Sci. Total Environ. 755, 142569 (2021).
Friedlingstein, P. et al. Global carbon budget 2019. Earth Syst. Sci. Data 11, 1783–1838 (2019).
Sitch, S. et al. Recent trends and drivers of regional sources and sinks of carbon dioxide. Biogeosciences 12, 653–679 (2015).
Viovy, N. CRUNCEP Version 7: Atmospheric Forcing Data for the Community Land Model. Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory (2018). Accessed 23 Aug 2021.
Dlugokencky et al. Atmospheric Carbon Dioxide Dry Air Mole Fractions from the NOAA GML Carbon Cycle Cooperative Global Air Sampling Network, 1968–2020. In (2021). Acessed 1 September 2023.
Hersbach, H. et al. ERA5 Monthly Averaged Data on Single Levels from 1940 to Present, Copernicus Climate Change Service (C3S) Climate Data Store (CDS) (2023). Accessed 1 Sep 2023.
Schneider, U. et al. GPCC Full Data Reanalysis Version 6.0 at 1.0°: Monthly Land-Surface Precipitation from Rain-Gauges built on GTS-based and Historic Data (2011). Accessed 1 Sep 2023.
Sun, Y. et al. OCO-2 advances photosynthesis observation from space via solar-induced chlorophyll fluorescence. Science 358, 5747 (2017).
Wilks, D. S. Chapter 12-Principal Component (EOF) Analysis. International Geophysics Vol. 100, 519–562 (Academic Press, 2011).
Hannachi, A., Jolliffe, I. T. & Stephenson, D. B. Empirical orthogonal functions and related techniques in atmospheric science: A review. Int. J. Climatol. 27, 1119–1152 (2007).
Hao, L. et al. Detection of the coupling between vegetation leaf area and climate in a multifunctional watershed, Northwestern China. Remote Sens. 8, 1032 (2016).
He, K., Li, W. & He, R. Variability of remotely sensed solar-induced chlorophyll fluorescence in relation to climate indices. Environments 9, 121 (2022).
Lepage, Y. A combination of Wilcoxon’s and Ansari-Bradley’s statistics. Biometrika 58, 213–217 (1971).
Yonetani, T. Detection of long term trend, cyclic variation and step-like change by the Lepage test. J. Meteorol. Soc. Japan. Ser. II(71), 415–418 (1993).
Kwon, M., Jhun, J.-G., Wang, B., An, S.-I. & Kug, J.-S. Decadal change in relationship between east Asian and WNP summer monsoons. Geophys. Res. Lett. 32, 16 (2005).
Yoo, Y.-E., Son, S.-W., Lee, J.-H. & Min, S.-K. Abrupt decrease of wintertime cold nights in Korea in the late 1980s. Asia-Pac. J. Atmos. Sci. 55, 31–39 (2019).
Liu, Y. et al. Changes in lake area in the Inner Mongolian plateau under climate change: The role of the Atlantic multidecadal oscillation and arctic sea ice. J. Clim. 33, 1335–1349 (2020).
Zhang, W. et al. Quantifying the contribution of climate change and human activities to biophysical parameters in an arid region. Ecol. Indic. 129, 107996 (2021).
Mauritsen, T. et al. Developments in the MPI-M Earth system model version 1.2 (MPI-ESM1.2) and its response to increasing CO2. J. Adv. Model. Earth Syst. 11, 998–1038 (2019).
Haverd, V. et al. A new version of the CABLE land surface model (Subversion revision r4601) incorporating land use and land cover change, woody vegetation demography, and a novel optimisation-based approach to plant coordination of photosynthesis. Geosci. Model Dev. 11, 2995–3026 (2018).
Melton, J. R. & Arora, V. K. Competition between plant functional types in the Canadian terrestrial ecosystem model (CTEM) v.2.0. Geosci. Model Dev. 9, 323–361 (2016).
Lawrence, D. M. et al. The community land model version 5: Description of new features, benchmarking, and impact of forcing uncertainty. J. Adv. Model. Earth Syst. 11, 4245–4287 (2019).
Tian, H. et al. North American terrestrial CO2 uptake largely offset by CH4 and N2O emissions: Toward a full accounting of the greenhouse gas budget. Clim. Change 129, 413–426 (2015).
Decharme, B. et al. Recent changes in the ISBA-CTRIP land surface system for use in the CNRM-CM6 climate model and in global off-line hydrological applications. J. Adv. Model. Earth Syst. 11, 1207–1252 (2019).
Sellar, A. A. et al. UKESM1: Description and evaluation of the UK earth system model. J. Adv. Model. Earth Syst. 11, 4513–4558 (2019).
Smith, B. et al. Implications of incorporating N cycling and N limitations on primary production in an individual-based dynamic vegetation model. Biogeosciences 11, 2027–2054 (2014).
Krinner, G. et al. A dynamic global vegetation model for studies of the coupled atmosphere-biosphere system. Glob. Biogeochem. Cycles 19, 2199 (2005).
Goll, D. S. et al. Carbon–nitrogen interactions in idealized simulations with JSBACH (version 310). Geosci. Model Dev. 10, 2009–2030 (2017).
Walker, A. P. et al. The impact of alternative trait-scaling hypotheses for the maximum photosynthetic carboxylation rate (Vcmax) on global gross primary production. New Phytol. 215, 1370–1386 (2017).
Kato, E., Kinoshita, T., Ito, A., Kawamiya, M. & Yamagata, Y. Evaluation of spatially explicit emission scenario of land-use change and biomass burning using a process-based biogeochemical model. J. Land Use Sci. 8, 104–122 (2013).
Acknowledgements
This research was supported by National Research Foundation of Korea (NRF) grant funded by the Korean government (NRF-2022R1A3B1077622 and NRF-2021M3I6A1086808), and Korea Environment Industry & Technology Institute (KEITI) through ‘Project for developing an observation-based GHG emissions geospatial information map', funded by Korea Ministry of Environment (MOE) (RS-2023-00232066). Jin-Soo Kim was supported by Korea Environmental Industry & Technology Institute (KEITI) through ‘Project for developing an observation-based GHG emissions geospatial information map’, funded by Korea Ministry of Environment (MOE) (RS-2023-00232066).
Author information
Authors and Affiliations
Contributions
D.-B.L. compiled the data, conducted analyses, prepared the figures, and wrote the manuscript. J.-S.Kim, S.-W.P., and J.-S.Kug designed the research and wrote the majority of the manuscript content. All of the authors discussed the study results and reviewed the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Lee, D., Kim, JS., Park, SW. et al. An abrupt shift in gross primary productivity over Eastern China-Mongolia and its inter-model diversity in land surface models. Sci Rep 13, 22971 (2023). https://doi.org/10.1038/s41598-023-49763-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-49763-1
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.