Greening of the land surface in the world's cold regions consistent with recent warming

Temporal variation characteristics in the association between climate and vegetation in Northwest China

Northwest China has undergone notable alterations in climate and vegetation growth in recent decades. Nevertheless, uncertainties persist concerning the response of different vegetation types to climate change and the underlying mechanisms. This study utilized the Normalized Difference Vegetation Index (NDVI) and three sets of meteorological data to investigate the interannual variations in the association between vegetation and climate (specifically precipitation and temperature) from 1982 to 2015. Several conclusions were drawn. (1) RNDVI-GP (relationship between Growing Season NDVI and precipitation) decreased significantly across all vegetation, while RNDVI-GT (relationship between Growing Season NDVI and temperature) showed an insignificant increase. (2) Trends of RNDVI-GP and RNDVI-GT exhibited great variations across various types of vegetation, with forests displaying notable downward trends in both indices. The grassland exhibited a declining trend in RNDVI-GP but an insignificant increase in RNDVI-GT, while no significant temporal changes in RNDVI-GP or RNDVI-GT were observed in the barren land. (3) The fluctuations in RNDVI-GP and RNDVI-GT closely aligned with variations in drought conditions. Specifically, in regions characterized by VPD (vapor pressure deficit) trends less than 0.02 hpa/yr, which are predominantly grasslands, a rise in SWV (soil water volume) tended to cause a reduction in RNDVI-GP but an increase in RNDVI-GT. However, a more negative trend in SWV was associated with a more negative trend in both RNDVI-GP and RNDVI-GT when the VPD trend exceeded 0.02 hPa/yr, primarily in forests. Our results underscore the variability in the relationship between climate change and vegetation across different vegetation types, as well as the role of drought in modulating these associations.

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.

Spatiotemporal patterns, sustainability, and primary drivers of NDVI-derived vegetation dynamics (2003-2022) in Nepal

Understanding the vegetation dynamics and their drivers in Nepal has significant scientific reference value for implementing sustainable ecological policies. This study provides a comprehensive analysis of the spatio-temporal variations in vegetation cover in Nepal from 2003 to 2022 using MODIS NDVI data and explores the effects of climatic factors and anthropogenic activities on vegetation. Mann-Kendall test was used to assess the significant trend in NDVI and was integrated with the Hurst exponent to predict future trends. The driving factors of NDVI dynamics were analyzed using Pearson's correlation, partial derivative, and residual analysis methods. The results indicate that over the last 20 years, Nepal has experienced an increasing trend in NDVI at 0.0013 year-1, with 80% of the surface area (vegetation cover) showing an increasing vegetation trend (~ 28% with a significant increase in vegetation). Temperature influenced vegetation dynamics in the higher elevation areas, while precipitation and human interventions influenced the lower elevation areas. The Hurst exponent analysis predicts an improvement in the vegetation cover (greening) for a larger area compared to vegetation degradation (browning). A significantly increased area of NDVI residuals indicates a positive anthropogenic influence on vegetation cover. Anthropogenic activities have a higher relative contribution to NDVI variation followed by temperature and then precipitation. The results of residual trend and Hurst analysis in different regions of Nepal help identify degraded areas, both in the present and future. This information can assist relevant authorities in implementing appropriate policies for a sustainable ecological environment.

Trends and drivers of Arctic-boreal fire intensity between 2003 and 2022

Climate change has disproportional effects on Arctic-boreal ecosystems, as the increase of air temperatures in these northern regions is several times higher than the global average. Ongoing warming and drying have resulted in recent record-breaking fire years in Arctic-boreal ecosystems, resulting in substantial carbon emissions that might accelerate climate change. While recent trends in Arctic-boreal burned area have been well documented, it is still unclear how fire intensity has changed. Fire intensity relates to the energy release from combustion and to a large extent drives the impact of a fire on the vegetation and soils, the emission of various gasses and the combustion completeness of different fuels. Here, we used the active fire product from the Moderate Resolution Imaging Spectroradiometer (MODIS) to examine trends in fire radiative power (FRP) over the entire Arctic-boreal region. We found a significant increase in annual median fire intensity between 2003 and 2022 in Eurasian boreal forests, for both daytime (increase of 0.392 MW/km2 per year, R2 = 0.56, p < 0.001) and nighttime fires (increase of 0.175 MW/km2 per year, R2 = 0.47, p < 0.001), while no general trend in FRP was observed in boreal North America. This increase in FRP in Eurasian boreal forests was strongly associated with simultaneous increases in air temperature, vapour pressure deficit, fire weather and fuel availability. We estimated that for Eurasia with each degree increase in air temperature, annual median daytime FRP increases with 1.58 MW/km2 in the tundra and 0.94 MW/km2 in the taiga. Climate change has thus resulted in a widespread and clear increase in fire intensity in central and eastern Eurasia while we could not discern clear trends in Arctic-boreal North America. Arctic-boreal fire intensity may further increase with climate change, with potentially major consequences for fire regimes, carbon emissions and society.

Alpine grassland greening on the Northern Tibetan Plateau driven by climate change and human activities considering extreme temperature and soil moisture

Alpine grassland is among the world's most vulnerable ecosystems, characterized by a high sensitivity to climate change (CC) and human activities (HA). Quantifying the relative contributions of CC and HA to grassland change plays a crucial role in safeguarding grassland ecological security and devising sustainable grassland management strategies. Although there were adequate studies focusing on the separate impacts of CC and HA on alpine ecosystem, insufficient attention has been given to investigating the effects of extreme temperatures and soil moisture. In this study, the spatiotemporal variations of alpine grassland were analyzed based on MODIS NDVI during the growing season from 2000 to 2020 in Naqu, using partial least squares regression and residual analysis methods to analyze the importance of climate factors and the impacts of CC and HA on grassland change. The results show that the NDVI during the growing season in Naqu exhibited an increasing trend of 0.0046/10a. At the biome scale, the most significant and rapid increase was observed in alpine desert and alpine desert grassland. Extreme temperature and soil moisture (SM) exerted a more significant importance on alpine grassland at whole scale. SM always showed a significant importance at biome and grid scale. The contributions of CC and HA to the change during the growing season were calculated as 0.0032/10a and 0.0015/10a, respectively, accounting for 68.05 % and 31.05 %. CC dominated the increase in NDVI during the growing season; HA contributed positively to NDVI in most areas of Naqu. The results are expected to enhance our understanding of grassland variations under CC and HA and provide a scientific basis for future ecological conservation in alpine regions.

An emergent constraint on the thermal sensitivity of photosynthesis and greenness in the high latitude northern forests

Despite the general consensus that the warming over the high latitudes northern forests (HLNF) has led to enhanced photosynthetic activity and contributed to the greening trend, isolating the impact of temperature increase on photosynthesis and greenness has been difficult due to the concurring influence of the CO2 fertilization effect. Here, using an ensemble of simulations from biogeochemical models that have contributed to the Trends in Net Land Atmosphere Carbon Exchange project (TRENDY), we identify an emergent relationship between the simulation of the climate-driven temporal changes in both gross primary productivity (GPP) and greenness (Leaf Area Index, LAI) and the model's spatial sensitivity of these quantities to growing-season (GS) temperature. Combined with spatially-resolved observations of LAI and GPP, we estimate that GS-LAI and GS-GPP increase by 17.0 ± 2.4% and 24.0 ± 3.0% per degree of warming, respectively. The observationally-derived sensitivities of LAI and GPP to temperature are about 40% and 71% higher, respectively, than the mean of the ensemble of simulations from TRENDY, primarily due to the model underestimation of the sensitivity of light use efficiency to temperature. We estimate that the regional mean GS-GPP increased 28.2 ± 5.1% between 1983-1986 and 2013-2016, much larger than the 5.8 ± 1.4% increase from the CO2 fertilization effect implied by Wenzel et al. This suggests that warming, not CO2 fertilization, is primarily responsible for the observed dramatic changes in the HLNF biosphere over the last century.

Global prediction of gross primary productivity under future climate change

The ecosystem gross primary productivity (GPP) is crucial to land-atmosphere carbon exchanges, and changes in global GPP as well as its influencing factors have been well studied in recent years. However, identifying the spatio-temporal variations of global GPP under future climate changes is still a challenging issue. This study aims to develop data-driven approach for predicting the global GPP as well as its monthly and annual variations up to the year 2100 under changing climate. Specifically, Catboost was employed to examine the potential relationship between the GPP and environmental factors, with climate variables, CO2 concentration and terrain attributes being selected as environmental factors. The predicted monthly and annual GPP from Coupled Model Intercomparison Project phase 6 (CMIP6) under future SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5 scenarios were analyzed. The results indicate that the global GPP is predicted to increase under the future climate change in the 21st century. The annual GPP is expected to be 115.122 Pg C, 116.537 Pg C, 117.626 Pg C, and 120.097 Pg C in 2100 under four future scenarios, and the predicted monthly GPP shows seasonal difference. Meanwhile, GPP tends to increase in the northern mid-high latitude regions and decrease in the equatorial regions. For the climate zones form Köppen-Geiger classification, the arid, cold, and polar zones present increased GPP, while GPP in the tropical zone will decrease in the future. Moreover, the high importance of climate variables in GPP prediction illustrates that the future climate change is the main driver of the global GPP dynamics. This study provides a basis for predicting how global GPP responds to future climate change in the coming decades, which contribute to understanding the interactions between vegetation and climate.

Comparative transcriptomics of the chilling stress response in two Asian mangrove species, Bruguiera gymnorhiza and Rhizophora apiculata

Low temperatures largely determine the geographic limits of plant species by reducing survival and growth. Inter-specific differences in the geographic distribution of mangrove species have been associated with cold tolerance, with exclusively tropical species being highly cold-sensitive and subtropical species being relatively cold-tolerant. To identify species-specific adaptations to low temperatures, we compared the chilling stress response of two widespread Indo-West Pacific mangrove species from Rhizophoraceae with differing latitudinal range limits-Bruguiera gymnorhiza (L.) Lam. ex Savigny (subtropical range limit) and Rhizophora apiculata Blume (tropical range limit). For both species, we measured the maximum photochemical efficiency of photosystem II (Fv/Fm) as a proxy for the physiological condition of the plants and examined gene expression profiles during chilling at 15 and 5 °C. At 15 °C, B. gymnorhiza maintained a significantly higher Fv/Fm than R. apiculata. However, at 5 °C, both species displayed equivalent Fv/Fm values. Thus, species-specific differences in chilling tolerance were only found at 15 °C, and both species were sensitive to chilling at 5 °C. At 15 °C, B. gymnorhiza downregulated genes related to the light reactions of photosynthesis and upregulated a gene involved in cyclic electron flow regulation, whereas R. apiculata downregulated more RuBisCo-related genes. At 5 °C, both species repressed genes related to CO2 assimilation. The downregulation of genes related to light absorption and upregulation of genes related to cyclic electron flow regulation are photoprotective mechanisms that likely contributed to the greater photosystem II photochemical efficiency of B. gymnorhiza at 15 °C. The results of this study provide evidence that the distributional range limits and potentially the expansion rates of plant species are associated with differences in the regulation of photosynthesis and photoprotective mechanisms under low temperatures.

Warming effects on grassland soil microbial communities are amplified in cool months

Global warming modulates soil respiration (RS) via microbial decomposition, which is seasonally dependent. Yet, the magnitude and direction of this modulation remain unclear, partly owing to the lack of knowledge on how microorganisms respond to seasonal changes. Here, we investigated the temporal dynamics of soil microbial communities over 12 consecutive months under experimental warming in a tallgrass prairie ecosystem. The interplay between warming and time altered (P < 0.05) the taxonomic and functional compositions of microbial communities. During the cool months (January to February and October to December), warming induced a soil microbiome with a higher genomic potential for carbon decomposition, community-level ribosomal RNA operon (rrn) copy numbers, and microbial metabolic quotients, suggesting that warming stimulated fast-growing microorganisms that enhanced carbon decomposition. Modeling analyses further showed that warming reduced the temperature sensitivity of microbial carbon use efficiency (CUE) by 28.7% when monthly average temperature was low, resulting in lower microbial CUE and higher heterotrophic respiration (Rh) potentials. Structural equation modeling showed that warming modulated both Rh and RS directly by altering soil temperature and indirectly by influencing microbial community traits, soil moisture, nitrate content, soil pH, and gross primary productivity. The modulation of Rh by warming was more pronounced in cooler months compared to warmer ones. Together, our findings reveal distinct warming-induced effects on microbial functional traits in cool months, challenging the norm of soil sampling only in the peak growing season, and advancing our mechanistic understanding of the seasonal pattern of RS and Rh sensitivity to warming.

Climatic limitations on grassland photosynthesis over the Tibetan Plateau shifted from temperature to water

Plant photosynthesis plays an essential role in regulating the global carbon cycle. Therefore, it is essential to understand the limitations imposed by climate on plant photosynthesis to comprehend the impacts of climate change on land carbon dynamics. In this study, taking gross primary productivity as a direct representation of photosynthesis, we employed a light use efficiency model (i.e., the revised EC-LUE) and factorial analysis method to quantify the spatiotemporal variation of temperature- and water-limitations on plant photosynthesis over the Tibetan Plateau (TP) grasslands during growing season (May to October) in 1983-2018. Results revealed a clear spatiotemporal pattern of the temperature- and water-limitations: temperature is the primary climatic limiting factor in the eastern TP, while water is the primary climatic limiting factor in the western TP; the water- and temperature-limitations prevail in summer and spring/autumn, respectively. The water- and temperature-limitations intensified and alleviated, respectively, during 1983 through 2018. There also was a widespread shift from temperature-limitation to water-limitation in the TP, particularly in midsummer (August). Our findings demonstrated the shifting relative importance of climatic limitations on plant photosynthesis under changing climate, which is crucial for predicting future terrestrial carbon cycle dynamics.

Spatiotemporal variations and its driving factors of NDVI in Northwest China during 2000-2021

Northwest China (WTL) is an essential ecological barrier zone of China, an important node of the "Silk Road Economic Belt," and a crucial bearing area for China's execution of the "One Road and One Belt" and "Going Global" strategies. However, its ecology is exceedingly fragile and particularly vulnerable to climate change and human interference. This study explored the spatiotemporal evolution characteristics of vegetation in WTL using NDVI data and investigated its drive mechanisms by geodetector, partial correlation analysis, and residual trend analysis methods. As well as forecasting the trend for vegetation changes. The findings demonstrated that (1) the change in NDVI manifested an overall improvement trend and the distribution in space of NDVI rose from the center to the periphery. 57.07% of the area had a sparse cover of vegetation (NDVI between 0 and 0.2). In addition, about 49% of regions had deterioration tendencies, which were mainly aggregated in HX, QCXDB, QCXDN, and the eastern of QCXQN and QCXXB. (2) The NDVI's shifting trend was unsustainability, and the region of uncertain future accounted for 57.45% of the total, with apparent unsustainability features. (3) The key parameters influencing NDVI spatial distribution were Pre (precipitation), vegetation type, land use type, and soil type. The interaction between two factors enhanced the influence of any single element, which appeared as bivariate and nonlinear enhancements. (4) Both climate variations and human activities have been recognized as key variables affecting NDVI growth. NDVI variance in 73.02% of areas was influenced by the combined effects of climate variations and human activities. However, human activities were the most influential element in NDVI growth, with the relative contributions of 80.28% (19.72% of which was caused by climate variations). These results can be conducive to deepening insights into the local vegetation status, identifying the mechanisms driving vegetation change, and providing scientific recommendations for WTL's ecosystem restoration measures based on actual situations.

Quantifying effects of cold acclimation and delayed springtime photosynthesis resumption in northern ecosystems

Land carbon dynamics in temperate and boreal ecosystems are sensitive to environmental change. Accurately simulating gross primary productivity (GPP) and its seasonality is key for reliable carbon cycle projections. However, significant biases have been found in early spring GPP simulations of northern forests, where observations often suggest a later resumption of photosynthetic activity than predicted by models. Here, we used eddy covariance-based GPP estimates from 39 forest sites that differ by their climate and dominant plant functional types. We used a mechanistic and an empirical light use efficiency (LUE) model to investigate the magnitude and environmental controls of delayed springtime photosynthesis resumption (DSPR) across sites. We found DSPR reduced ecosystem LUE by 30-70% at many, but not all site-years during spring. A significant depression of LUE was found not only in coniferous but also at deciduous forests and was related to combined high radiation and low minimum temperatures. By embedding cold-acclimation effects on LUE that considers the delayed effects of minimum temperatures, initial model bias in simulated springtime GPP was effectively resolved. This provides an approach to improve GPP estimates by considering physiological acclimation and enables more reliable simulations of photosynthesis in northern forests and projections in a warming climate.

Increased impact of the El Niño-Southern Oscillation on global vegetation under future warming environment

There are broad effects of vegetation changes on regional climate, carbon budget, the water cycle, and ecosystems' productivity. Therefore, further knowledge of the drivers of future vegetation changes is critical to mitigate the influences of global warming. The El Niño-Southern Oscillation (ENSO) is a major mode of interannual climate variability and is likely to affect vegetation on the global scale. Nonetheless, little is known about the causal impacts of ENSO on future vegetation cover with changes in land use and a warming environment. Here, we examined the connections between ENSO and vegetation using leaf area index (LAI) data over the period 2015-2100 from Coupled Modeling Intercomparison Project Phase 6. Our findings indicate that, compared with the historical period 1915-2000, the vegetated areas influenced by ENSO are projected to rise by approximately 55.2% and 20.7% during the twenty-first century of the scenarios SSP2-4.5 and SSP5-8.5, respectively. Though uncertainty for the causal link between ENSO and vegetation changes remains in several regions (i.e., parts of North America, southern Australia, and western Asia), ENSO signature on LAI variations is robust over northern Australia, Amazonia, and parts of Southeast Asia. These results indicate that the influences of ENSO on global vegetation may strengthen in the future.

CO2 fertilization contributed more than half of the observed forest biomass increase in northern extra-tropical land

The existence of a large-biomass carbon (C) sink in Northern Hemisphere extra-tropical ecosystems (NHee) is well-established, but the relative contribution of different potential drivers remains highly uncertain. Here we isolated the historical role of carbon dioxide (CO2 ) fertilization by integrating estimates from 24 CO2 -enrichment experiments, an ensemble of 10 dynamic global vegetation models (DGVMs) and two observation-based biomass datasets. Application of the emergent constraint technique revealed that DGVMs underestimated the historical response of plant biomass to increasing [CO2 ] in forests (β Forest Mod ) but overestimated the response in grasslands (β Grass Mod ) since the 1850s. Combining the constrainedβ Forest Mod (0.86 ± 0.28 kg C m-2  [100 ppm]-1 ) with observed forest biomass changes derived from inventories and satellites, we identified that CO2 fertilization alone accounted for more than half (54 ± 18% and 64 ± 21%, respectively) of the increase in biomass C storage since the 1990s. Our results indicate that CO2 fertilization dominated the forest biomass C sink over the past decades, and provide an essential step toward better understanding the key role of forests in land-based policies for mitigating climate change.

Effects of climate change on vegetation dynamics of the Qinghai-Tibet Plateau, a causality analysis using empirical dynamic modeling

Given the vital role of the Qinghai-Tibet Plateau (QTP) as water tower in Asia and regulator for regional and even global climate, the relationship between climate change and vegetation dynamics on it has received considerable focused attention. Climate change may influence the vegetation growth on the plateau, but clear empirical evidence of such causal linkages is sparse. Herein, using datasets CRU-TS v4.04 and AVHHR NDVI from 1981 to 2019, we quantify causal effects of climate factors on vegetation dynamics with an empirical dynamical model (EDM) -- a nonlinear dynamical systems analysis approach based on state-space reconstruction rather than correlation. Results showed the following: (1) climate change promotes the growth of vegetation on the QTP, and specifically, this favorable influence of temperature is stronger than precipitation's; (2) the direction and strength of climate effects on vegetation varied over time, and the effects are seasonally different; (3) a significant increase in temperature and a slight increase in precipitation are beneficial to vegetation growth, specifically, NDVI will increase within 2% in the next 40 years with the climate trend of warming and humidity. Besides the above results, another interesting finding is that the two seasons in which precipitation strongly influence vegetation in the Three-River Source region (part of the QTP) are spring and winter. This study provides insights into the mechanisms by which climate change affects vegetation growth on the QTP, aiding in the modeling of vegetation dynamics in future scenarios.

Increasing temperature regulates the advance of peak photosynthesis timing in the boreal ecosystem

The shift in vegetation phenology is an essential indicator of global climate change. Numerous researches based on reflectance-based vegetation index data have explored the changes in the start (SOS) and end (EOS) of vegetation life events at long time scales, while a huge discrepancy existed between the phenological metrics of vegetation structure and function. The peak photosynthesis timing (PPT), which is crucial in regulating terrestrial ecosystem carbon balance, has not received much attention. Using two global reconstructed solar-induced chlorophyll fluorescence data (CSIF and GOSIF) directly associated with vegetation photosynthesis, the spatio-temporal dynamics in PPT as well as the key environmental controls across the boreal ecosystem during 2001-2019 were systematically explored. Multi-year mean pattern showed that PPT mainly appeared in the first half of July. Compared to the northern Eurasia, later PPT appeared in the northern North America continent for about 4-5 days. Meanwhile, spatial trend in PPT exhibited an advanced trend during the last two decades. Especially, shrubland and grassland were obvious among all biomes. Spatial partial correlation analysis revealed that preseason temperature was the dominant environmental driver of PPT trends, occupying 81.32% and 78.04% of the total pixels of PPT_CSIF and PPT_GOSIF, respectively. Attribution analysis by ridge regression again emphasized the largest contribution of temperature to PPT dynamics in the boreal ecosystem by 52.22% (PPT_CSIF) and 46.59% (PPT_GOSIF), followed by radiation (PPT_CSIF: 24.44%; PPT_GOSIF: 28.66%) and precipitation (PPT_CSIF: 23.34%; PPT_GOSIF: 24.75%). These results have significant implications for deepening our understanding between vegetation photosynthetic phenology and carbon cycling with respect to future climate change in the boreal ecosystem.

Himalayan alpine ecohydrology: An urgent scientific concern in a changing climate

Climate change is projected to have important impacts on snow and vegetation distribution in global mountains. Despite this, the coupling of ecological shifts and hydrological processes within alpine zones has not attracted significant scientific attention. As the largest and one of the most climatically sensitive mountain systems, we argue that Himalayan alpine ecohydrological processes require urgent scientific attention because up to 1.6 billion people rely on water supplies from the mountains. We review studies from global mountain systems to highlight the importance of considering ecohydrological impacts within Himalayan alpine zones (4100-6000 m.a.s.l), explaining mechanisms for interactions between snow and dwarf plants. Our findings highlight the paucity of monitoring stations within Himalayan alpine systems. We suggest that it is likely that alpine ecological shifts will impact hydrological processes, but we found that specific mechanisms and functional relationships are missing for Himalayan systems, so the strength and direction of ecohydrological relationships is currently unknown. We advocate for more purposeful and widespread monitoring efforts below glaciers and above the treeline, calling for new experiments to query the role of small plants within the Himalayan alpine hydrological system. We outline the need for community engagement with alpine ecohydrological experiments, and we explain how new snow and vegetation products derived from remote sensing observations have the potential to improve scientific understanding of the interacting effects of warming and ecohydrological factors in this sensitive region.

Optimality principles explaining divergent responses of alpine vegetation to environmental change

Recent increases in vegetation greenness over much of the world reflect increasing CO₂ globally and warming in cold areas. However, the strength of the response to both CO₂ and warming in those areas appears to be declining for unclear reasons, contributing to large uncertainties in predicting how vegetation will respond to future global changes. Here, we investigated the changes of satellite-observed peak season absorbed photosynthetically active radiation (F_max ) on the Tibetan Plateau between 1982 and 2016. Although climate trends are similar across the Plateau, we identified robust divergent responses (a greening of 0.31 ± 0.14% year^-1 in drier regions and a browning of 0.12 ± 0.08% year^-1 in wetter regions). Using an eco-evolutionary optimality (EEO) concept of plant acclimation/adaptation, we propose a parsimonious modelling framework that quantitatively explains these changes in terms of water and energy limitations. Our model captured the variations in F_max with a correlation coefficient (r) of .76 and a root mean squared error of .12 and predicted the divergent trends of greening (0.32 ± 0.19% year^-1 ) and browning (0.07 ± 0.06% year^-1 ). We also predicted the observed reduced sensitivities of F_max to precipitation and temperature. The model allows us to explain these changes: Enhanced growing season cumulative radiation has opposite effects on water use and energy uptake. Increased precipitation has an overwhelmingly positive effect in drier regions, whereas warming reduces F_max in wetter regions by increasing the cost of building and maintaining leaf area. Rising CO₂ stimulates vegetation growth by enhancing water-use efficiency, but its effect on photosynthesis saturates. The large decrease in the sensitivity of vegetation to climate reflects a shift from water to energy limitation. Our study demonstrates the potential of EEO approaches to reveal the mechanisms underlying recent trends in vegetation greenness and provides further insight into the response of alpine ecosystems to ongoing climate change.

Quantifying the direct effects of long-term dynamic land use intensity on vegetation change and its interacted effects with economic development and climate change in jiangsu, China

Vegetation change reflects sensitive responses of ecosystem environment to global climate change as well as land use. It is well known that land use type and its transformation affect vegetation change. However, how the changes in land use intensity (LUI) within different land use types impact vegetation and the interactions with other drivers remain poorly understood. We measured the LUI of Jiangsu Province, China, within the main land use types in 1995, 2000, 2005, 2010, 2015 and 2018 by combining remote sensing-based land use data with representative county scale economic and social indicators. Structural equation models (SEMs) were built to quantify the influences of long term LUI on vegetation change interacting with economic development, climate change and topographical conditions in transformed land, cropland, rural settlements and urbanized land, respectively. Seventy percent of significant vegetation change existed in non-transformed land use types. Although the area with a vegetation greening trend is larger than that with a vegetation browning trend, the vegetation browning areas is prominent in urbanized lands and some croplands in south basins. The constructed SEMs suggested the dominant negative effect of fast economic development regardless of land use types, while LUI played important and different direct and indirect effects on affecting vegetation change significantly interacting with economic development and climate change in different land use types. The LUI increasing led a vegetation greening in cropland, and stronger than climate warming with both positive direct and indirect effects for influencing climate change. The LUI change took negative effects on vegetation change in rural and urban areas, while a positive indirect effect of LUI increasing in urbanized land signaled the positive results of human managements. We then provided some land use-specific suggestions on basin scale for land management in Jiangsu. Our results highlight the necessity of long-term LUI quantification and promote the understanding of its effects on vegetation change interacted with other drivers within different land use types. This can be very helpful for sustainable land use and managements in regions with fast economic development.

Remote sensing hydrological indication: Responses of hydrological processes to vegetation cover change in mid-latitude mountainous regions

Hydrological processes in mid-latitude mountainous regions are greatly affected by changes in vegetation cover that induced by the climate change. However, studies on hydrological processes in mountainous regions are limited, because of difficulties in building and maintaining basin-wide representative hydrological stations. In this study, a new method, remote sensing technology for monitoring river discharge by combining satellite remote sensing, unmanned aerial vehicles and hydrological surveying, was used for evaluating the runoff processes in the Changbai Mountains, one of the mid-latitude mountainous regions in the eastern part of Northeast China. Based on this method, the impact of vegetation cover change on hydrological processes was revealed by combining the data of hydrological processes, meteorology, and vegetation cover. The results showed a decreasing trend in the monitored river discharge from 2000 to 2021, with an average rate of -5.13 × 10⁵ m³ yr^-1. At the monitoring section mainly influenced by precipitation, the precipitation-induced proportion of changes in river discharge to annual average river discharge and its change significance was only 6.5 % and 0.23, respectively, showing the precipitation change was not the cause for the decrease in river discharge. A negative impact of evapotranspiration on river discharge was found, and the decrease in river discharge was proven to be caused by the increasing evapotranspiration, which was induced by the drastically increased vegetation cover under a warming climate. Our findings suggested that increases in vegetation cover due to climate change could reshape hydrological processes in mid-latitude mountainous regions, leading to an increase in evapotranspiration and a subsequent decrease in river discharge.

Climate Controls on the Spatial Variability of Vegetation Greenup Rate across Ecosystems in Northern Hemisphere

Variations in individual phenological events in response to global change have received considerable attentions. However, the development of phenological stages is relatively neglected, especially based on in situ observation data. In this study, the rate of vegetation greenup (Vgreenup) across the Northern Hemisphere was examined for different plant functional types (PFTs) by using eddy covariance flux data from 40 sites (417 site-years). Then, the controls of climatic variables on the spatial distribution of Vgreenup across PFTs were further investigated. The mean Vgreenup was 0.22 ± 0.11 g C m^-2 day^-2 across all sites, with the largest and lowest values observed in cropland and evergreen needle-leaf forest, respectively. A strong latitude dependence by Vgreenup was observed in both Europe and North America. The spatial variations of Vgreenup were jointly regulated by the duration of greenup (Dgreenup) and the amplitude of greenup (Agreenup). However, the predominant factor was Dgreenup in Europe, which changed to Agreenup in North America. Spring climatic factors exerted significant influences on the spatial distribution of Vgreenup across PFTs. Specifically, increasing temperature tended to shorten Dgreenup and promote Agreenup simultaneously, resulting in an acceleration of Vgreenup. Dryness had a depression effect on Vgreenup for the whole study area, as exhibited by a lower Vgreenup with increasing vapor pressure deficit or decreasing soil moisture. However, Vgreenup in North America was only significantly and positively correlated with temperature. Without the limitation of other climatic factors, the temperature sensitivity of Vgreenup was higher in North America (0.021 g C m^-2 day^-2 °C^-1) than in Europe (0.015 g C m^-2 day^-2 °C^-1). This study provides new cognitions for Vgreenup dynamics from in situ observations in complement to satellite observations, which can improve our understanding of terrestrial carbon cycles.

Greening or browning? The macro variation and drivers of different vegetation types on the Qinghai-Tibetan Plateau from 2000 to 2021

Vegetation greenness is one of the main indicators to characterize changes in terrestrial ecosystems. China has implemented a few large-scale ecological restoration programs on the Qinghai-Tibetan Plateau (QTP) to reverse the trend of ecosystem degradation. Although the effectiveness of these programs is beginning to show, the mechanisms of vegetation degradation under climate change and human activities are still controversial. Existing studies have mostly focused on changes in overall vegetation change, with less attention on the drivers of change in different vegetation types. In this study, earth satellite observation records were used to robustly map changes in vegetation greenness on the QTP from 2000 to 2021. The random forest (RF) algorithm was further used to detect the drivers of greenness browning on the QTP as a whole and in seven different vegetation types. The results show that an overall trend of greening in all seven vegetation types on the QTP over a 21-year period. The area of greening was 46.54×104 km2, and browning was 5.32×104 km2, representing a quarter and 2.86% of the natural vegetation area, respectively. The results of the browning driver analysis show that areas with high altitude, reduced annual precipitation, high intensity of human activity, average annual maximum and average annual minimum precipitation of approximately 500 mm are most susceptible to browning on the QTP. For the seven different vegetation types, their top 6 most important browning drivers and the ranking of drivers differed. DEM and precipitation changes are important drivers of browning for seven vegetation types. These results reflect the latest spatial and temporal dynamics of vegetation on the QTP and highlight the common and characteristic browning drivers of vegetation ecosystems. They provide support for understanding the response of different vegetation to natural and human impacts and for further implementation of site-specific restoration measures.

Quantifying the contributions of climate change and human activities to vegetation dynamic in China based on multiple indices

Distinguishing the respective roles of climate change and anthropogenic activities can provide crucial information for sustainable management of the environment. Here, using the residual trend method (RESTREND), which measures the residue of the actual and potential trends of vegetation, we quantified the relative contributions of human activities (e.g., ecological restoration, overgrazing, and urbanization) and climate change (the warmer and wetter trend) to vegetation dynamics in China during 1988-2018 based on multiple vegetation indices, including the vegetation optical depth (Ku-VOD, C-VOD), normalized difference vegetation index (NDVI), and gross primary productivity (GPP). The results showed that the VOD, NDVI, and GPP exhibited overall increasing trends during 1988-2018. Human activities contributed >70% to the increases in NDVI and GPP in China, whereas a counterbalanced contribution of human activities and climate change was identified for the VOD dynamics (51% vs. 49%). Regions with high contributions from human activities to NDVI, GPP, and VOD were located in northeastern, southern, central, and northwestern China. In northern China, the positive impacts of human activities on NDVI (78%) and BEPS-GPP (83%) were greater than those of climate change. In contrast, human activities contributed 96% to the decrease in Ku-VOD over the same period. Before 2000, climate change promoted increases in GPP and NDVI in most regions of southern China. The increasing rates of GPP and NDVI accelerated after 2000 due to afforestation. However, human activities like overgrazing and urbanization have led to decreases in Ku-VOD in northern and southwestern China, and in C-VOD in northeastern, eastern, central, southwestern, and southern China. In all, the relative roles of climate and human factors varied in different regions when NDVI, GPP, or VOD were individually considered. Our results highlighted that the regional-scale vegetation conditions should be taken into full account to achieve sustainable management of ecosystems.

Precipitation and soil nutrients determine the spatial variability of grassland productivity at large scales in China

Changes in net primary productivity (NPP) to global change have been studied, yet the relative impacts of global change on grassland productivity at large scales remain poorly understood. Using 182 grassland samples established in 17 alpine meadows (AM) and 21 desert steppes (DS) in China, we show that NPP of AM was significantly higher than that of DS. NPP increased significantly with increasing leaf nitrogen content (LN) and leaf phosphorus content (LP) but decreased significantly with increasing leaf dry matter content (LDMC). Among all abiotic factors, soil nutrient factor was the dominant factor affecting the variation of NPP of AM, while the NPP of DS was mainly influenced by the changing of precipitation. All abiotic factors accounted for 62.4% of the spatial variation in the NPP of AM, which was higher than the ability to explain the spatial variation in the NPP of DS (43.5%). Leaf traits together with soil nutrients and climatic factors determined the changes of the grassland productivity, but the relative contributions varied somewhat among different grassland types. We quantified the effects of biotic and abiotic factors on grassland NPP, and provided theoretical guidance for predicting the impacts of global change on the NPP of grasslands.

Quantifying Northern High Latitude Gross Primary Productivity (GPP) Using Carbonyl Sulfide (OCS)

The northern high latitude (NHL, 40°N to 90°N) is where the second peak region of gross primary productivity (GPP) other than the tropics. The summer NHL GPP is about 80% of the tropical peak, but both regions are still highly uncertain (Norton et al. 2019, https://doi.org/10.5194/bg-16-3069-2019). Carbonyl sulfide (OCS) provides an important proxy for photosynthetic carbon uptake. Here we optimize the OCS plant uptake fluxes across the NHL by fitting atmospheric concentration simulation with the GEOS-CHEM global transport model to the aircraft profiles acquired over Alaska during NASA's Carbon in Arctic Reservoirs Vulnerability Experiment (2012-2015). We use the empirical biome-specific linear relationship between OCS plant uptake flux and GPP to derive the six plant uptake OCS fluxes from different GPP data. Such GPP-based fluxes are used to drive the concentration simulations. We evaluate the simulations against the independent observations at two ground sites of Alaska. The optimized OCS fluxes suggest the NHL plant uptake OCS flux of -247 Gg S year^-1, about 25% stronger than the ensemble mean of the six GPP-based OCS fluxes. GPP-based OCS fluxes systematically underestimate the peak growing season across the NHL, while a subset of models predict early start of season in Alaska, consistent with previous studies of net ecosystem exchange. The OCS optimized GPP of 34 PgC yr^-1 for NHL is also about 25% more than the ensembles mean from six GPP data. Further work is needed to fully understand the environmental and biotic drivers and quantify their rate of photosynthetic carbon uptake in Arctic ecosystems.

Assessment of climate change effects on vegetation and river hydrology in a semi-arid river basin

Climate change plays a key role in changing vegetation productivity dynamics, which ultimately affect the hydrological cycle of a watershed through evapotranspiration (ET). Trends and correlation analysis were conducted to investigate vegetation responses across the whole Upper Jhelum River Basin (UJRB) in the northeast of Pakistan using the normalized difference vegetation index (NDVI), climate variables, and river flow data at inter-annual/monthly scales between 1982 and 2015. The spatial variability in trends calculated with the Mann-Kendall (MK) trend test on NDVI and climate data was assessed considering five dominant land use/cover types. The inter-annual NDVI in four out of five vegetation types showed a consistent increase over the 34-year study period; the exception was for herbaceous vegetation (HV), which increased until the end of the 1990s and then decreased slightly in subsequent years. In spring, significant (p<0.05) increasing trends were found in the NDVI of all vegetation types. Minimum temperature (Tmin) showed a significant increase during spring, while maximum temperature (Tmax) decreased significantly during summer. Average annual increase in Tmin (1.54°C) was much higher than Tmax (0.37°C) over 34 years in the UJRB. Hence, Tmin appears to have an enhancing effect on vegetation productivity over the UJRB. A significant increase in NDVI, Tmin and Tmax during spring may have contributed to reductions in spring river flow by enhancing evapotranspiration observed in the watershed of UJRB. These findings provide valuable information to improve our knowledge and understanding about the interlinkages between vegetation, climate and river flow at a watershed scale.

Revealing four decades of snow cover dynamics in the Hindu Kush Himalaya

Knowledge about the distribution and dynamics of seasonal snow cover (SSC) is of high importance for climate studies, hydrology or hazards assessment. SSC varies considerably across the Hindu Kush Himalaya both in space and time. Previous studies focused on regional investigations or the influence of snow melt on the local hydrological system. Here, we present a systematic assessment of metrics to evaluate SSC dynamics for the entire HKH at regional and basin scale based on AVHRR GAC data at a 0.05° spatial and daily temporal resolution. Our findings are based on a unique four-decade satellite-based time series of snow cover information. We reveal strong variability of SSC at all time scales. We find significantly decreasing SSC trends in individual summer and winter months and a declining tendency from mid-spring to mid-fall, indicating a shift in seasonality. Thanks to this uniquely spatio-temporally resolved long-term data basis, we can particularly highlight the unique temporally variable character of seasonal snow cover and its cross-disciplinary importance for mountain ecosystems and downstream regions.

Direct and indirect impacts of urbanization on vegetation growth across the world's cities

Urban environments, regarded as "harbingers" of future global change, may exert positive or negative impacts on urban vegetation growth. Because of limited ground-based experiments, the responses of vegetation to urbanization and its associated controlling factors at the global scale remain poorly understood. Here, we use satellite observations from 2001 to 2018 to quantify direct and indirect impacts of urbanization on vegetation growth in 672 worldwide cities. After controlling for the negative direct impact of urbanization on vegetation growth, we find a widespread positive indirect effect that has been increasing over time. These indirect effects depend on urban development intensity, population density, and background climate, with more pronounced positive effects in cities with cold and arid environments. We further show that vegetation responses to urbanization are modulated by a cities' developmental status. Our findings have important implications for understanding urbanization-induced impacts on vegetation and future sustainable urban development.

Ecological Dynamics and Regeneration Expansion of Treeline Ecotones in Response to Climate Change in Northern Bhutan Himalayas

The alpine treeline ecotones are an early indicator of vegetation’s response to changes in climate, and the advancement of diffuse treeline ecotones has been associated with mean annual warming temperatures. However, the knowledge of how tree demographic size, age and population distribution, and regeneration decrease with increasing elevation and mean annual temperature remain fragmentary in Bhutan. There was no explanation of how treelines migrate in response to the climate. Therefore, the objectives of this study were to investigate tree demographic size and age and population distribution, as well as the regeneration expansion of treeline ecotones of Abies densa trees in response to climate change. Demographic data from thirty transect bands from treeline ecotones and reconstructed mean annual temperatures from tree-rings were used. Regression analysis was used to establish a relationship between elevation/temperature and demographic tree size and age, as well as to determine recruitment frequency distributions and whether these could be driven by climate change. The tree demography indicated that the treeline ecotone in our sampling site is temperature limited. Hence, cooler temperatures at higher elevations should drive decreases in basal diameter, age and recruitment frequencies. From the dendroecological analysis, the diffuse treeline ecotones appear to be climbing on average 1.00 m per year in Northern Bhutan. We also found that the recruitment frequency has increased over recent years (1850–2017), as temperatures continue to rise. The thermal treeline ecotones will be likely to serve as a line of bioclimatic reference against which other zones of bioclimate can be defined. With documented responses of treeline ecotones toward mean annual temperatures, the expectation is that additional warming will continue to influence regeneration expansion in the future. This dynamic response of treeline ecotones towards the climate acts as an indicator of climate change. Information about climbing treelines and altered ecotones should be a vital part of the material for decision makers to consider, to assess impacts and threats to Himalayan alpine biota.

Dispersal and fire limit Arctic shrub expansion

Arctic shrub expansion alters carbon budgets, albedo, and warming rates in high latitudes but remains challenging to predict due to unclear underlying controls. Observational studies and models typically use relationships between observed shrub presence and current environmental suitability (bioclimate and topography) to predict shrub expansion, while omitting shrub demographic processes and non-stationary response to changing climate. Here, we use high-resolution satellite imagery across Alaska and western Canada to show that observed shrub expansion has not been controlled by environmental suitability during 1984-2014, but can only be explained by considering seed dispersal and fire. These findings provide the impetus for better observations of recruitment and for incorporating currently underrepresented processes of seed dispersal and fire in land models to project shrub expansion and climate feedbacks. Integrating these dynamic processes with projected fire extent and climate, we estimate shrubs will expand into 25% of the non-shrub tundra by 2100, in contrast to 39% predicted based on increasing environmental suitability alone. Thus, using environmental suitability alone likely overestimates and misrepresents shrub expansion pattern and its associated carbon sink.

Satellite observations document trends consistent with a boreal forest biome shift

The boreal forest biome is a major component of Earth's biosphere and climate system that is projected to shift northward due to continued climate change over the coming century. Indicators of a biome shift will likely first be evident along the climatic margins of the boreal forest and include changes in vegetation productivity, mortality, and recruitment, as well as overall vegetation greenness. However, the extent to which a biome shift is already underway remains unclear because of the local nature of most field studies, sparsity of systematic ground-based ecological monitoring, and reliance on coarse resolution satellite observations. Here, we evaluated early indicators of a boreal forest biome shift using four decades of moderate resolution (30 m) satellite observations and biogeoclimatic spatial datasets. Specifically, we quantified interannual trends in annual maximum vegetation greenness using an ensemble of vegetation indices derived from Landsat observations at 100,000 sample sites in areas without signs of recent disturbance. We found vegetation greenness increased (greened) at 38 [29, 42] % and 22 [15, 26] % of sample sites from 1985 to 2019 and 2000 to 2019, whereas vegetation greenness decreased (browned) at 13 [9, 15] % and 15 [13, 19] % of sample sites during these respective periods [95% Monte Carlo confidence intervals]. Greening was thus 3.0 [2.6, 3.5] and 1.5 [0.8, 2.0] times more common than browning and primarily occurred in cold sparsely treed areas with high soil nitrogen and moderate summer warming. Conversely, browning primarily occurred in the climatically warmest margins of both the boreal forest biome and major forest types (e.g., evergreen conifer forests), especially in densely treed areas where summers became warmer and drier. These macroecological trends reflect underlying shifts in vegetation productivity, mortality, and recruitment that are consistent with early stages of a boreal biome shift.

Phenology-based seasonal terrestrial vegetation growth response to climate variability with consideration of cumulative effect and biological carryover

Vegetation growth is influenced not only by climate variability but also by its past states. However, the differences in the degree of the climate variability and past states affecting vegetation growth over seasons are still poorly understood, particularly given the cumulative climate effects. Relying on the Normalized Difference Vegetation Index (NDVI) data from 1982 to 2014, the vegetation growing season was decomposed into three periods (sub-seasons) - green-up (GSgp), maturity (GSmp), and senescence (GSsp) - following a phenology-based definition. A distributed lag model was then utilized to analyze the time-lag effect of vegetation growth response to climatic factors including precipitation, temperature, and solar radiation during each sub-season. On this basis, the relative importance of climatic factors and vegetation growth carryover (VGC) effect on vegetation growth was quantified at the phenology-based seasonal scale. Results showed that the longest peak lag of precipitation, temperature, and solar radiation occurred in the GSmp, GSsp, and GSgp, with 1.27 (1.13 SD), 0.89 (1.02 SD), and 0.80 (1.04 SD) months, respectively. The influence of climate variability was strongest in the GSgp, and diminished over the season, while the opposite for the VGC effect. The relative influence of each climatic factor also varied between sub-seasons. Vegetation in more than 58% of areas was more affected by temperature in the GSgp, and the proportion decreased to 34.00% and 31.78% in the GSmp and GSsp, respectively. Precipitation and solar radiation acted as the dominant climatic factors in only 28.80% and 20.88% of vegetation areas in the GSgp, but they increased to 35.21%, 32.61% in the GSmp, and 38.20%, 30.02% in the GSsp, respectively. The increased regions influenced by precipitation were mainly in dry areas especially for the boreal and cool temperate climate zones, while increased regions influenced by solar radiation were primarily located in moist areas of mid-high latitudes of the Northern Hemisphere. By introducing the cumulative climate effect, our findings highlight seasonal patterns of vegetation growth affected by climate variability and the VGC effect. The results provide a more comprehensive perspective on climate-vegetation interactions, which may help us to accurately forecast future vegetation growth under accelerating global warming.

An earlier start of the thermal growing season enhances tree growth in cold humid areas but not in dry areas

Climatic warming alters the onset, duration and cessation of the vegetative season. While previous studies have shown a tight link between thermal conditions and leaf phenology, less is known about the impacts of phenological changes on tree growth. Here, we assessed the relationships between the start of the thermal growing season and tree growth across the extratropical Northern Hemisphere using 3,451 tree-ring chronologies and daily climatic data for 1948-2014. An earlier start of the thermal growing season promoted growth in regions with high ratios of precipitation to temperature but limited growth in cold-dry regions. Path analyses indicated that an earlier start of the thermal growing season enhanced growth primarily by alleviating thermal limitations on wood formation in boreal forests and by lengthening the period of growth in temperate and Mediterranean forests. Semi-arid and dry subalpine forests, however, did not benefit from an earlier onset of growth and a longer growing season, presumably due to associated water loss and/or more frequent early spring frosts. These emergent patterns of how climatic impacts on wood phenology affect tree growth at regional to hemispheric scales hint at how future phenological changes may affect the carbon sequestration capacity of extratropical forest ecosystems.

Assessing dynamic vegetation model parameter uncertainty across Alaskan arctic tundra plant communities

As the Arctic region moves into uncharted territory under a warming climate, it is important to refine the terrestrial biosphere models (TBMs) that help us understand and predict change. One fundamental uncertainty in TBMs relates to model parameters, configuration variables internal to the model whose value can be estimated from data. We incorporate a version of the Terrestrial Ecosystem Model (TEM) developed for arctic ecosystems into the Predictive Ecosystem Analyzer (PEcAn) framework. PEcAn treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We examined how variation in 21 parameters in the equation for gross primary production influenced model sensitivity and uncertainty in terms of two carbon fluxes (net primary productivity and heterotrophic respiration) and two carbon (C) pools (vegetation C and soil C). We set up different parameterizations of TEM across a range of tundra types (tussock tundra, heath tundra, wet sedge tundra, and shrub tundra) in northern Alaska, along a latitudinal transect extending from the coastal plain near Utqiaġvik to the southern foothills of the Brooks Range, to the Seward Peninsula. TEM was most sensitive to parameters related to the temperature regulation of photosynthesis. Model uncertainty was mostly due to parameters related to leaf area, temperature regulation of photosynthesis, and the stomatal responses to ambient light conditions. Our analysis also showed that sensitivity and uncertainty to a given parameter varied spatially. At some sites, model sensitivity and uncertainty tended to be connected to a wider range of parameters, underlining the importance of assessing tundra community processes across environmental gradients or geographic locations. Generally, across sites, the flux of net primary productivity (NPP) and pool of vegetation C had about equal uncertainty, while heterotrophic respiration had higher uncertainty than the pool of soil C. Our study illustrates the complexity inherent in evaluating parameter uncertainty across highly heterogeneous arctic tundra plant communities. It also provides a framework for iteratively testing how newly collected field data related to key parameters may result in more effective forecasting of Arctic change.

The Human Impact on All Soil-Forming Factors during the Anthropocene

Soil-the thin outer skin of the Earth's land-is a critical and fragile natural resource. Soil is the basis for almost all global agriculture and the medium in which most terrestrial biological activity occurs. Here, we reconsider the five forming factors of soil originally suggested more than a century ago (parent material, time, climate, topography, and organisms) and updated over the years to add human activity as the sixth forming factor. We demonstrate how present anthropogenic activity has become the leading component influencing each one of the original forming factors. We thus propose that, starting from the Anthropocene, human activity should no longer be considered as a separate forming factor but rather a main driving force of each of the five original ones. We suggest that the importance of soil and the strong direct and indirect effects of anthropogenic factors on soil-forming factors should be considered together to ensure sustainability of this critical resource.

A Pixel-Based Vegetation Greenness Trend Analysis over the Russian Tundra with All Available Landsat Data from 1984 to 2018

As Arctic warming continues, its impact on vegetation greenness is complex, variable and inherently scale-dependent. Studies with multiple spatial resolution satellite observations, with 30 m resolution included, on tundra greenness have been implemented all over the North American tundra. However, finer resolution studies on the greenness trends in the Russian tundra have only been carried out at a limited local or regional scale and the spatial heterogeneity of the trend remains unclear. Here, we analyzed the fine spatial resolution dataset Landsat archive from 1984 to 2018 over the entire Russian tundra and produced pixel-by-pixel greenness trend maps with the support of Google Earth Engine (GEE). The entire Russian tundra was divided into six geographical regions based on World Wildlife Fund (WWF) ecoregions. A Theil–Sen regression (TSR) was used for the trend identification and the changed pixels with a significance level p < 0.05 were retained in the final results for a subsequent greening/browning trend analysis. Our results indicated that: (1) the number of valid Landsat observations was spatially varied. The Western and Eastern European Tundras (WET and EET) had denser observations than other regions, which enabled a trend analysis during the whole study period from 1984 to 2018; (2) the most significant greening occurred in the Yamal-Gydan tundra (WET), Bering tundra and Chukchi Peninsula tundra (CT) during 1984–2018. The EET had a greening trend of 2.3% and 6.6% and the WET of 3.4% and 18% during 1984–1999 and 2000–2018, respectively. The area of browning trend was relatively low when we first masked the surface water bodies out before the trend analysis; and (3) the Landsat-based greenness trend was broadly similar to the AVHRR-based trend over the entire region but AVHRR retrieved more browning areas due to spectral mixing adjacent effects. Higher resolution images and field measurement studies are strongly needed to understand the vegetation trend over the Russian tundra ecosystem.

The tempo of greening in the European Alps: Spatial variations on a common theme

The long-term increase in satellite-based proxies of vegetation cover is a well-documented response of seasonally snow-covered ecosystems to climate warming. However, observed greening trends are far from uniform, and substantial uncertainty remains concerning the underlying causes of this spatial variability. Here, we processed surface reflectance of the moderate resolution imaging spectroradiometer (MODIS) to investigate trends and drivers of changes in the annual peak values of the Normalized Difference Vegetation Index (NDVI). Our study focuses on above-treeline ecosystems in the European Alps. NDVI changes in these ecosystems are highly sensitive to land cover and biomass changes and are marginally affected by anthropogenic disturbances. We observed widespread greening for the 2000-2020 period, a pattern that is consistent with the overall increase in summer temperature. At the local scale, the spatial variability of greening was mainly due to the preferential response of north-facing slopes between 1900 and 2400 m. Using high-resolution imagery, we noticed that the presence of screes and outcrops locally magnified this response. At the regional scale, we identified hotspots of greening where vegetation cover is sparser than expected given the elevation and exposure. Most of these hotspots experienced delayed snow melt and green-up dates in recent years. We conclude that the ongoing greening in the Alps primarily reflects the high responsiveness of sparsely vegetated ecosystems that are able to benefit the most from temperature and water-related habitat amelioration above treeline.

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

[Figure: see text].

Inconsistency of Global Vegetation Dynamics Driven by Climate Change: Evidences from Spatial Regression

Global greening over the past 30 years since 1980s has been confirmed by numerous studies. However, a single-dimensional indicator and non-spatial modelling approaches might exacerbate uncertainties in our understanding of global change. Thus, comprehensive monitoring for vegetation’s various properties and spatially explicit models are required. In this study, we used the newest enhanced vegetation index (EVI) products of Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 6 to detect the inconsistency trend of annual peak and average global vegetation growth using the Mann–Kendall test method. We explored the climatic factors that affect vegetation growth change from 2001 to 2018 using the spatial lag model (SLM), spatial error model (SEM) and geographically weighted regression model (GWR). The results showed that EVImax and EVImean in global vegetated areas consistently showed linear increasing trends during 2001–2018, with the global averaged trend of 0.0022 yr−1 (p < 0.05) and 0.0030 yr−1 (p < 0.05). Greening mainly occurred in the croplands and forests of China, India, North America and Europe, while browning was almost in the grasslands of Brazil and Africa (18.16% vs. 3.08% and 40.73% vs. 2.45%). In addition, 32.47% of the global vegetated area experienced inconsistent trends in EVImax and EVImean. Overall, precipitation and mean temperature had positive impacts on vegetation variation, while potential evapotranspiration and vapour pressure had negative impacts. The GWR revealed that the responses of EVI to climate change were inconsistent in an arid or humid area, in cropland or grassland. Climate change could affect vegetation characteristics by changing plant phenology, consequently rendering the inconsistency between peak and mean greening. In addition, anthropogenic activities, including land cover change and land use management, also could lead to the differences between annual peak and mean vegetation variations.

BVOC Emissions From a Subarctic Ecosystem, as Controlled by Insect Herbivore Pressure and Temperature

Abstract

The biogenic volatile organic compounds, BVOCs have a central role in ecosystem–atmosphere interactions. High-latitude ecosystems are facing increasing temperatures and insect herbivore pressure, which may affect their BVOC emission rates, but evidence and predictions of changes remain scattered. We studied the long-term effects of + 3 °C warming and reduced insect herbivory (achieved through insecticide sprayings) on mid- and late summer BVOC emissions from field layer vegetation, supplemented with birch saplings, and the underlying soil in Subarctic mountain birch forest in Finland in 2017–2018. Reduced insect herbivory decreased leaf damage by 58–67% and total ecosystem BVOC emissions by 44–72%. Of the BVOC groups, total sesquiterpenes had 70–80% lower emissions with reduced herbivory, and in 2017 the decrease was greater in warmed plots (89% decrease) than in ambient plots (34% decrease). While non-standardized total BVOC, monoterpene, sesquiterpene and GLV emissions showed instant positive responses to increasing chamber air temperature in midsummer samplings, the long-term warming treatment effects on standardized emissions mainly appeared as changes in the compound structure of BVOC blends and varied with compounds and sampling times. Our results suggest that the effects of climate warming on the total quantity of BVOC emissions will in Subarctic ecosystems be, over and above the instant temperature effects, mediated through changes in insect herbivore pressure rather than plant growth. If insect herbivore numbers will increase as predicted under climate warming, our results forecast herbivory-induced increases in the quantity of Subarctic BVOC emissions.

Graphic Abstract

Modeling vegetation greenness and its climate sensitivity with deep-learning technology

Climate sensitivity of vegetation has long been explored using statistical or process-based models. However, great uncertainties still remain due to the methodologies' deficiency in capturing the complex interactions between climate and vegetation. Here, we developed global gridded climate-vegetation models based on long short-term memory (LSTM) network, which is a powerful deep-learning algorithm for long-time series modeling, to achieve accurate vegetation monitoring and investigate the complex relationship between climate and vegetation. We selected the normalized difference vegetation index (NDVI) that represents vegetation greenness as model outputs. The climate data (monthly temperature and precipitation) were used as inputs. We trained the networks with data from 1982 to 2003, and the data from 2004 to 2015 were used to validate the models. Error analysis and sensitivity analysis were performed to assess the model errors and investigate the sensitivity of global vegetation to climate change. Results show that models based on deep learning are very effective in simulating and predicting the vegetation greenness dynamics. For models training, the root mean square error (RMSE) is <0.01. Model validation also assure the accuracy of our models. Furthermore, sensitivity analysis of models revealed a spatial pattern of global vegetation to climate, which provides us a new way to investigate the climate sensitivity of vegetation. Our study suggests that it is a good way to integrate deep-learning method to monitor the vegetation change under global change. In the future, we can explore more complex climatic and ecological systems with deep learning and coupling with certain physical process to better understand the nature.

Increased association between climate change and vegetation index variation promotes the coupling of dominant factors and vegetation growth

Vegetation productivity dynamics are closely related to climate change, and water availability determines vegetation growth in water-limited ecosystems. Nevertheless, how changes in the interactions between climatic factors and vegetation activity variation regulate the relationship between their trends remains unclear. The Normalized Difference Vegetation Index (NDVI) is an effective proxy of vegetation growth. First, we investigated the NDVI trends, and the results revealed a vegetation activity with weaker greening and greater spatial heterogeneity after an obvious land-cover breakpoint in 1999 compared with that before 1999 in northwest China. Notably, the Loess Plateau greatly led the greenness trends, but the Tibet Plateau showed mean browning after 1999, which implied that the coupling of climate change and vegetation trends varied with spatio-temporal changes. Subsequently, using the Geographical Detector Method (GDM), we quantified and compared the association between climate change and the interannual variability of NDVI in the two stages. Vegetation productivity variation is more closely related to changes in climatic factors after 1999 compared with that before 1999. Precipitation (PPT) and vapor pressure deficit (VPD) are the primary constraints to vegetation growth in both stages. Patterns in NDVI trend increases are consistent with those of increased PPT and decreased VPD and vice versa after 1999. However, the same patterns were not observed before 1999 because of the weak association between climate change and NDVI variation. This implicated a great significance of the association between climate change and changes in vegetation activity for the prediction of potential carbon sequestration due to the shift of dominant factors and their trends under future climate change.

Integrating the evidence for a terrestrial carbon sink caused by increasing atmospheric CO2

Atmospheric carbon dioxide concentration ([CO₂ ]) is increasing, which increases leaf-scale photosynthesis and intrinsic water-use efficiency. These direct responses have the potential to increase plant growth, vegetation biomass, and soil organic matter; transferring carbon from the atmosphere into terrestrial ecosystems (a carbon sink). A substantial global terrestrial carbon sink would slow the rate of [CO₂ ] increase and thus climate change. However, ecosystem CO₂ responses are complex or confounded by concurrent changes in multiple agents of global change and evidence for a [CO₂ ]-driven terrestrial carbon sink can appear contradictory. Here we synthesize theory and broad, multidisciplinary evidence for the effects of increasing [CO₂ ] (iCO₂ ) on the global terrestrial carbon sink. Evidence suggests a substantial increase in global photosynthesis since pre-industrial times. Established theory, supported by experiments, indicates that iCO₂ is likely responsible for about half of the increase. Global carbon budgeting, atmospheric data, and forest inventories indicate a historical carbon sink, and these apparent iCO₂ responses are high in comparison to experiments and predictions from theory. Plant mortality and soil carbon iCO₂ responses are highly uncertain. In conclusion, a range of evidence supports a positive terrestrial carbon sink in response to iCO₂ , albeit with uncertain magnitude and strong suggestion of a role for additional agents of global change.

Seasonal biological carryover dominates northern vegetation growth

The state of ecosystems is influenced strongly by their past, and describing this carryover effect is important to accurately forecast their future behaviors. However, the strength and persistence of this carryover effect on ecosystem dynamics in comparison to that of simultaneous environmental drivers are still poorly understood. Here, we show that vegetation growth carryover (VGC), defined as the effect of present states of vegetation on subsequent growth, exerts strong positive impacts on seasonal vegetation growth over the Northern Hemisphere. In particular, this VGC of early growing-season vegetation growth is even stronger than past and co-occurring climate on determining peak-to-late season vegetation growth, and is the primary contributor to the recently observed annual greening trend. The effect of seasonal VGC persists into the subsequent year but not further. Current process-based ecosystem models greatly underestimate the VGC effect, and may therefore underestimate the CO2 sequestration potential of northern vegetation under future warming.