Temperature increase reduces global yields of major crops in four independent estimates.

Wheat, rice, maize, and soybean provide two-thirds of human caloric intake. Assessing the impact of global temperature increase on production of these crops is therefore critical to maintaining global food supply, but different studies have yielded different results. Here, we investigated the impacts of temperature on yields of the four crops by compiling extensive published results from four analytical methods: global grid-based and local point-based models, statistical regressions, and field-warming experiments. Results from the different methods consistently showed negative temperature impacts on crop yield at the global scale, generally underpinned by similar impacts at country and site scales. Without CO2 fertilization, effective adaptation, and genetic improvement, each degree-Celsius increase in global mean temperature would, on average, reduce global yields of wheat by 6.0%, rice by 3.2%, maize by 7.4%, and soybean by 3.1%. Results are highly heterogeneous across crops and geographical areas, with some positive impact estimates. Multimethod analyses improved the confidence in assessments of future climate impacts on global major crops and suggest crop- and region-specific adaptation strategies to ensure food security for an increasing world population.

Dominant regions and drivers of the variability of the global land carbon sink across timescales.

Net biome productivity (NBP) dominates the observed large variation of atmospheric CO2 annual increase over the last five decades. However, the dominant regions controlling inter-annual to multi-decadal variability of global NBP are still controversial (semi-arid regions vs. temperate or tropical forests). By developing a theory for partitioning the variance of NBP into the contributions of net primary production (NPP) and heterotrophic respiration (Rh ) at different timescales, and using both observation-based atmospheric CO2 inversion product and the outputs of 10 process-based terrestrial ecosystem models forced by 110-year observational climate, we tried to reconcile the controversy by showing that semi-arid lands dominate the variability of global NBP at inter-annual (<10 years) and tropical forests dominate at multi-decadal scales (>30 years). Results further indicate that global NBP variability is dominated by the NPP component at inter-annual timescales, and is progressively controlled by Rh with increasing timescale. Multi-decadal NBP variations of tropical rainforests are modulated by the Pacific Decadal Oscillation (PDO) through its significant influences on both temperature and precipitation. This study calls for long-term observations for the decadal or longer fluctuations in carbon fluxes to gain insights on the future evolution of global NBP, particularly in the tropical forests that dominate the decadal variability of land carbon uptake and are more effective for climate mitigation.

Inconsistencies of interannual variability and trends in long-term satellite leaf area index products.

Understanding the long-term performance of global satellite leaf area index (LAI) products is important for global change research. However, few effort has been devoted to evaluating the long-term time-series consistencies of LAI products. This study compared four long-term LAI products (GLASS, GLOBMAP, LAI3g, and TCDR) in terms of trends, interannual variabilities, and uncertainty variations from 1982 through 2011. This study also used four ancillary LAI products (GEOV1, MERIS, MODIS C5, and MODIS C6) from 2003 through 2011 to help clarify the performances of the four long-term LAI products. In general, there were marked discrepancies between the four long-term LAI products. During the pre-MODIS period (1982-1999), both linear trends and interannual variabilities of global mean LAI followed the order GLASS>LAI3g>TCDR>GLOBMAP. The GLASS linear trend and interannual variability were almost 4.5 times those of GLOBMAP. During the overlap period (2003-2011), GLASS and GLOBMAP exhibited a decreasing trend, TCDR no trend, and LAI3g an increasing trend. GEOV1, MERIS, and MODIS C6 also exhibited an increasing trend, but to a much smaller extent than that from LAI3g. During both periods, the R2 of detrended anomalies between the four long-term LAI products was smaller than 0.4 for most regions. Interannual variabilities of the four long-term LAI products were considerably different over the two periods, and the differences followed the order GLASS>LAI3g>TCDR>GLOBMAP. Uncertainty variations quantified by a collocation error model followed the same order. Our results indicate that the four long-term LAI products were neither intraconsistent over time nor interconsistent with each other. These inconsistencies may be due to NOAA satellite orbit changes and MODIS sensor degradation. Caution should be used in the interpretation of global changes derived from the four long-term LAI products.

Attribution of seasonal leaf area index trends in the northern latitudes with "optimally" integrated ecosystem models.

Significant increases in remotely sensed vegetation indices in the northern latitudes since the 1980s have been detected and attributed at annual and growing season scales. However, we presently lack a systematic understanding of how vegetation responds to asymmetric seasonal environmental changes. In this study, we first investigated trends in the seasonal mean leaf area index (LAI) at northern latitudes (north of 30°N) between 1982 and 2009 using three remotely sensed long-term LAI data sets. The most significant LAI increases occurred in summer (0.009 m2  m-2  year-1 , p < .01), followed by autumn (0.005 m2  m-2  year-1 , p < .01) and spring (0.003 m2 m-2  year-1 , p < .01). We then quantified the contribution of elevating atmospheric CO2 concentration (eCO2 ), climate change, nitrogen deposition, and land cover change to seasonal LAI increases based on factorial simulations from 10 state-of-the-art ecosystem models. Unlike previous studies that used multimodel ensemble mean (MME), we used the Bayesian model averaging (BMA) to optimize the integration of model ensemble. The optimally integrated ensemble LAI changes are significantly closer to the observed seasonal LAI changes than the traditional MME results. The BMA factorial simulations suggest that eCO2 provides the greatest contribution to increasing LAI trends in all seasons (0.003-0.007 m2  m-2  year-1 ), and is the main factor driving asymmetric seasonal LAI trends. Climate change controls the spatial pattern of seasonal LAI trends and dominates the increase in seasonal LAI in the northern high latitudes. The effects of nitrogen deposition and land use change are relatively small in all seasons (around 0.0002 m2  m-2  year-1 and 0.0001-0.001 m2  m-2  year-1 , respectively). Our analysis of the seasonal LAI responses to the interactions between seasonal changes in environmental factors offers a new perspective on the response of global vegetation to environmental changes.

Nonlinear variations of forest leaf area index over China during 1982-2010 based on EEMD method.

Variations in leaf area index (LAI) are critical to research on forest ecosystem structure and function, especially carbon and water cycle, and their responses to climate change. Using the ensemble empirical mode decomposition (EEMD) method and global inventory modeling and mapping studies (GIMMS) LAI3g dataset from 1982 to 2010, we analyzed the nonlinear feature and spatial difference of forest LAI variability over China for the past 29 years in this paper. Results indicated that the national-averaged forest LAI was characterized by quasi-3- and quasi-7-year oscillations, which generally exhibited a rising trend with an increasing rate. When compared with 1982, forest LAI change by 2010 was more evident than that by 1990 and 2000. The largest increment of forest LAI occurred in Central and South China, while along the southeastern coastal areas LAI increased at the fastest pace. During the study period, forest LAI experienced from decrease to increase or vice versa across much of China and varied monotonically for only a few areas. Focusing on regional-averaged trend processes, almost all eco-geographical regions showed continuously increasing trends in forest LAI with different magnitudes and speeds, other than tropical humid region and temperate humid/subhumid region, where LAI decreased initially and increased afterwards.

Delayed autumn phenology in the Northern Hemisphere is related to change in both climate and spring phenology.

The timing of the end of the vegetation growing season (EOS) plays a key role in terrestrial ecosystem carbon and nutrient cycles. Autumn phenology is, however, still poorly understood, and previous studies generally focused on few species or were very limited in scale. In this study, we applied four methods to extract EOS dates from NDVI records between 1982 and 2011 for the Northern Hemisphere, and determined the temporal correlations between EOS and environmental factors (i.e., temperature, precipitation and insolation), as well as the correlation between spring and autumn phenology, using partial correlation analyses. Overall, we observed a trend toward later EOS in ~70% of the pixels in Northern Hemisphere, with a mean rate of 0.18 ± 0.38 days yr-1 . Warming preseason temperature was positively associated with the rate of EOS in most of our study area, except for arid/semi-arid regions, where the precipitation sum played a dominant positive role. Interestingly, increased preseason insolation sum might also lead to a later date of EOS. In addition to the climatic effects on EOS, we found an influence of spring vegetation green-up dates on EOS, albeit biome dependent. Our study, therefore, suggests that both environmental factors and spring phenology should be included in the modeling of EOS to improve the predictions of autumn phenology as well as our understanding of the global carbon and nutrient balances.

Global impacts of the 1980s regime shift.

Despite evidence from a number of Earth systems that abrupt temporal changes known as regime shifts are important, their nature, scale and mechanisms remain poorly documented and understood. Applying principal component analysis, change-point analysis and a sequential t-test analysis of regime shifts to 72 time series, we confirm that the 1980s regime shift represented a major change in the Earth's biophysical systems from the upper atmosphere to the depths of the ocean and from the Arctic to the Antarctic, and occurred at slightly different times around the world. Using historical climate model simulations from the Coupled Model Intercomparison Project Phase 5 (CMIP5) and statistical modelling of historical temperatures, we then demonstrate that this event was triggered by rapid global warming from anthropogenic plus natural forcing, the latter associated with the recovery from the El Chichón volcanic eruption. The shift in temperature that occurred at this time is hypothesized as the main forcing for a cascade of abrupt environmental changes. Within the context of the last century or more, the 1980s event was unique in terms of its global scope and scale; our observed consequences imply that if unavoidable natural events such as major volcanic eruptions interact with anthropogenic warming unforeseen multiplier effects may occur.

Vegetation canopy structure mediates the response of gross primary production to environmental drivers across multiple temporal scales.

Gross primary production (GPP) is a critical component of the global carbon cycle and plays a significant role in the terrestrial carbon budget. The impact of environmental factors on GPP can occur through both direct (by influencing photosynthetic efficiency) and indirect (through the modulation of vegetation structure) pathways, but the extent to which these mechanisms contribute has been seldom quantified. In this study, we used structural equation modeling and observations from the FLUXNET network to investigate the direct and indirect effects of environmental factors on terrestrial ecosystem GPP at multiple temporal scales. We found that canopy structure, represented by leaf area index (LAI), is a crucial intermediate factor in the GPP response to environmental drivers. Environmental factors affect GPP indirectly by altering canopy structure, and the relative proportion of indirect effects decreased with increasing LAI. The study also identified different effects of environmental factors on GPP across time scales. At the half-hourly time scale, radiation was the primary driver of GPP. In contrast, the influences of temperature and vapor pressure deficit took on greater prominence at longer time scales. About half of the total effect of temperature on GPP was indirect through the regulation of canopy structure, and the indirect effect increased with increasing time scale (GPPNT-based models: 0.135 (half-hourly) vs. 0.171 (daily) vs. 0.189 (weekly) vs. 0.217 (monthly); GPPDT-based models: 0.139 vs. 0.170 vs. 0.187 vs. 0.215; all values were reported in gC m-2 d-1 °C-1, P < 0.001); while the indirect effect of radiation on GPP was comparatively lower, accounting for less than a quarter of the total effect. Furthermore, we observed a direct, negative-to-positive impact of precipitation on GPP across timescales. These findings provide crucial information on the interplay between environmental factors and LAI on GPP and enable a deeper understanding of the driving mechanisms of GPP.

A global dataset of the fraction of absorbed photosynthetically active radiation for 1982-2022.

The fraction of absorbed photosynthetically active radiation (FPAR) is an essential biophysical parameter that characterizes the structure and function of terrestrial ecosystems. Despite the extensive utilization of several satellite-derived FPAR products, notable temporal inconsistencies within each product have been underscored. Here, the new generation of the GIMMS FPAR product, GIMMS FPAR4g, was developed using a combination of a machine learning algorithm and a pixel-wise multi-sensor records integration approach. PKU GIMMS NDVI, which eliminates the orbital drift and sensor degradation issues, was used as the data source. Comparisons with ground-based measurements indicate root mean square errors ranging from 0.10 to 0.14 with R-squared ranging from 0.73 to 0.87. More importantly, our product demonstrates remarkable spatiotemporal coherence and continuity, revealing a persistent terrestrial darkening over the past four decades (0.0004 yr-1, p < 0.001). The GIMMS FPAR4g, available for half-month intervals at a spatial resolution of 1/12° from 1982 to 2022, promises to be a valuable asset for in-depth analyses of vegetation structures and functions spanning the last 40 years.

Global critical soil moisture thresholds of plant water stress.

During extensive periods without rain, known as dry-downs, decreasing soil moisture (SM) induces plant water stress at the point when it limits evapotranspiration, defining a critical SM threshold (θcrit). Better quantification of θcrit is needed for improving future projections of climate and water resources, food production, and ecosystem vulnerability. Here, we combine systematic satellite observations of the diurnal amplitude of land surface temperature (dLST) and SM during dry-downs, corroborated by in-situ data from flux towers, to generate the observation-based global map of θcrit. We find an average global θcrit of 0.19 m3/m3, varying from 0.12 m3/m3 in arid ecosystems to 0.26 m3/m3 in humid ecosystems. θcrit simulated by Earth System Models is overestimated in dry areas and underestimated in wet areas. The global observed pattern of θcrit reflects plant adaptation to soil available water and atmospheric demand. Using explainable machine learning, we show that aridity index, leaf area and soil texture are the most influential drivers. Moreover, we show that the annual fraction of days with water stress, when SM stays below θcrit, has increased in the past four decades. Our results have important implications for understanding the inception of water stress in models and identifying SM tipping points.

Nitrogen addition delays the emergence of an aridity-induced threshold for plant biomass.

Crossing certain aridity thresholds in global drylands can lead to abrupt decays of ecosystem attributes such as plant productivity, potentially causing land degradation and desertification. It is largely unknown, however, whether these thresholds can be altered by other key global change drivers known to affect the water-use efficiency and productivity of vegetation, such as elevated CO2 and nitrogen (N). Using >5000 empirical measurements of plant biomass, we showed that crossing an aridity (1-precipitation/potential evapotranspiration) threshold of ∼0.50, which marks the transition from dry sub-humid to semi-arid climates, led to abrupt declines in aboveground biomass (AGB) and progressive increases in root:shoot ratios, thus importantly affecting carbon stocks and their distribution. N addition significantly increased AGB and delayed the emergence of its aridity threshold from 0.49 to 0.55 (P < 0.05). By coupling remote sensing estimates of leaf area index with simulations from multiple models, we found that CO2 enrichment did not alter the observed aridity threshold. By 2100, and under the RCP 8.5 scenario, we forecast a 0.3% net increase in the global land area exceeding the aridity threshold detected under a scenario that includes N deposition, in comparison to a 2.9% net increase if the N effect is not considered. Our study thus indicates that N addition could mitigate to a great extent the negative impact of increasing aridity on plant biomass in drylands. These findings are critical for improving forecasts of abrupt vegetation changes in response to ongoing global environmental change.

[Division of winter wheat yield estimation by remote sensing based on MODIS EVI time series data and spectral angle clustering].

Crop yield estimation division is the basis of crop yield estimation; it provides an important scientific basis for estimation research and practice. In the paper, MODIS EVI time-series data during winter wheat growth period is selected as the division data; JiangSu province is study area; A division method combined of advanced spectral angle mapping(SVM) and K-means clustering is presented, and tested in winter wheat yield estimation by remote sensing. The results show that: division method of spectral angle clustering can take full advantage of crop growth process that is reflected by MODIS time series data, and can fully reflect region differences of winter wheat that is brought by climate difference. Compared with the traditional division method, yield estimation result based on division result of spectral angle clustering has higher R2 (0.702 6 than 0.624 8) and lower RMSE (343.34 than 381.34 kg x hm(-2)), reflecting the advantages of the new division method in the winter wheat yield estimation. The division method in the paper only use convenient obtaining time-series remote sensing data of low-resolution as division data, can divide winter wheat into similar and well characterized region, accuracy and stability of yield estimation model is also very good, which provides an efficient way for winter wheat estimation by remote sensing, and is conducive to winter wheat yield estimation.

Seasonal peak photosynthesis is hindered by late canopy development in northern ecosystems.

The seasonal dynamics of the vegetation canopy strongly regulate the surface energy balance and terrestrial carbon fluxes, providing feedbacks to climate change. Whether the seasonal timing of maximum canopy structure was optimized to achieve a maximum photosynthetic carbon uptake is still not clear due to the complex interactions between abiotic and biotic factors. We used two solar-induced chlorophyll fluorescence datasets as proxies for photosynthesis and the normalized difference vegetation index and leaf area index products derived from the moderate resolution imaging spectroradiometer as proxies for canopy structure, to characterize the connection between their seasonal peak timings from 2000 to 2018. We found that the seasonal peak was earlier for photosynthesis than for canopy structure in >87.5% of the northern vegetated area, probably leading to a suboptimal maximum seasonal photosynthesis. This mismatch in peak timing significantly increased during the study period, mainly due to the increasing atmospheric CO2, and its spatial variation was mainly explained by climatic variables (43.7%) and nutrient limitations (29.6%). State-of-the-art ecosystem models overestimated this mismatch in peak timing by simulating a delayed seasonal peak of canopy development. These results highlight the importance of incorporating the mechanisms of vegetation canopy dynamics to accurately predict the maximum potential terrestrial uptake of carbon under global environmental change.

Comment on "Recent global decline of CO2 fertilization effects on vegetation photosynthesis".

Wang et al. (Research Articles, 11 December 2020, p. 1295) reported a large decrease in CO2 fertilization effect (CFE) across the globe during the period 1982­2015 and suggested that ecosystem models underestimate the rate of CFE decline. We find that their claims are artifacts of incorrect processing of satellite data and problematic methods for deriving and comparing CFE between satellite data and model simulations.

Shift from feeding to sustainably nourishing urban China: A crossing-disciplinary methodology for global environment-food-health nexus.

Dietary change is a win-win opportunity to address the nexus of health and the environment. To prevent city dwellers from developing non-communicable diseases, in 2013, China updated the 2000 version of nutrition-based dietary reference intake (DRI) guidelines. However, whether the DRI guidelines have a positive effect on the environment is not well understood. Here, we explored the systematic effects of urbanization on China's health and environmental nexus based on survey data. Then, we optimized the diets of 18 age-gender groups to reduce carbon emissions, water consumption, and land use while meeting the healthy nutrition goals of both DRI guidelines. The results showed that the optimal diets based on the DRI 2013 outperformed these on DRI 2000 in improving China's environmental sustainability, although these diets did not always perform better at an individual scale. Our findings suggest that dietary changes can reduce carbon, water, and ecological footprints by 24%, 15%, and 22% in 2050, respectively; however, the differences in age-specific and gender-specific health goals cannot be neglected.

China and India lead in greening of the world through land-use management.

Satellite data show increasing leaf area of vegetation due to direct (human land-use management) and indirect factors (climate change, CO2 fertilization, nitrogen deposition, recovery from natural disturbances, etc.). Among these, climate change and CO2 fertilization effect seem to be the dominant drivers. However, recent satellite data (2000-2017) reveal a greening pattern that is strikingly prominent in China and India, and overlapping with croplands world-wide. China alone accounts for 25% of the global net increase in leaf area with only 6.6% of global vegetated area. The greening in China is from forests (42%) and croplands (32%), but in India is mostly from croplands (82%) with minor contribution from forests (4.4%). China is engineering ambitious programs to conserve and expand forests with the goal of mitigating land degradation, air pollution and climate change. Food production in China and India has increased by over 35% since 2000 mostly due to increasing harvested area through multiple cropping facilitated by fertilizer use and surface/ground-water irrigation. Our results indicate that the direct factor is a key driver of the "Greening Earth", accounting for over a third, and likely more, of the observed net increase in green leaf area. They highlight the need for realistic representation of human land-use practices in Earth system models.

Velocity of change in vegetation productivity over northern high latitudes.

Warming is projected to increase the productivity of northern ecosystems. However, knowledge on whether the northward displacement of vegetation productivity isolines matches that of temperature isolines is still limited. Here we compared changes in the spatial patterns of vegetation productivity and temperature using the velocity of change concept, which expresses these two variables in the same unit of displacement per time. We show that across northern regions (>50° N), the average velocity of change in growing-season normalized difference vegetation index (NDVIGS, an indicator of vegetation productivity; 2.8 ± 1.1 km yr-1) is lower than that of growing-season mean temperature (T GS; 5.4 ± 1.0 km yr-1). In fact, the NDVIGS velocity was less than half of the T GS velocity in more than half of the study area, indicating that the northward movement of productivity isolines is much slower than that of temperature isolines across the majority of northern regions (about 80% of the area showed faster changes in temperature than productivity isolines). We tentatively attribute this mismatch between the velocities of productivity and temperature to the effects of limited resource availability and vegetation acclimation mechanisms. Analyses of ecosystem model simulations further suggested that limited nitrogen availability is a crucial obstacle for vegetation to track the warming trend.

The Effect of Afforestation on Soil Moisture Content in Northeastern China.

Widespread afforestation programs sequester carbon from the atmosphere and mitigate the rising of atmospheric carbon dioxide (CO2). Meanwhile, afforestation carbon sequestration may cost soil water. However, changes in soil moisture content (SMC) after large-scale afforestation or reforestation have rarely been quantified. In this study, we measured changes in SMC following afforestation using a paired plots method with data from 757 plots in Northeastern China. We found a marginally significant decline in soil moisture content of the top 1-m soil (SMC0-1m) after afforestation (P = 0.08) at the regional scale. The SMC responses to afforestation also vary across species. For example, significant SMC decrease are found for Populus spp. plantations (P < 0.05) and plantations of Pinus sylvestris var. mongolica (P < 0.05). Splitting the first meter of the soil profile into different depth intervals revealed that SMC declined significantly in shallow layers (0-30 cm) for Populus spp. and Pinus sylvestris var. mongolica. We also found that when SMC in the control exceeded a specific threshold, SMC for all five tree species considered tended to decrease, suggesting that the effects of afforestation on soil hydrology vary across different regions.