Climate-induced tree-mortality pulses are obscured by broad-scale and long-term greening

Vegetation greening has been suggested to be a dominant trend over recent decades, but severe pulses of tree mortality in forests after droughts and heatwaves have also been extensively reported. These observations raise the question of to what extent the observed severe pulses of tree mortality induced by climate could affect overall vegetation greenness across spatial grains and temporal extents. To address this issue, here we analyse three satellite-based datasets of detrended growing-season normalized difference vegetation index (NDVIGS) with spatial resolutions ranging from 30 m to 8 km for 1,303 field-documented sites experiencing severe drought- or heat-induced tree-mortality events around the globe. We find that severe tree-mortality events have distinctive but localized imprints on vegetation greenness over annual timescales, which are obscured by broad-scale and long-term greening. Specifically, although anomalies in NDVIGS (ΔNDVI) are negative during tree-mortality years, this reduction diminishes at coarser spatial resolutions (that is, 250 m and 8 km). Notably, tree-mortality-induced reductions in NDVIGS (|ΔNDVI|) at 30-m resolution are negatively related to native plant species richness and forest height, whereas topographic heterogeneity is the major factor affecting ΔNDVI differences across various spatial grain sizes. Over time periods of a decade or longer, greening consistently dominates all spatial resolutions. The findings underscore the fundamental importance of spatio-temporal scales for cohesively understanding the effects of climate change on forest productivity and tree mortality under both gradual and abrupt changes.

Climate change increases carbon allocation to leaves in early leaf green-up

Global greening, characterized by an increase in leaf area index (LAI), implies an increase in foliar carbon (C). Whether this increase in foliar C under climate change is due to higher photosynthesis or to higher allocation of C to leaves remains unknown. Here, we explored the trends in foliar C accumulation and allocation during leaf green-up from 2000 to 2017 using satellite-derived LAI and solar-induced chlorophyll fluorescence (SIF) across the Northern Hemisphere. The accumulation of foliar C accelerated in the early green-up period due to both increased photosynthesis and higher foliar C allocation driven by climate change. In the late stage of green-up, however, we detected decreasing trends in foliar C accumulation and foliar C allocation. Such stage-dependent trends in the accumulation and allocation of foliar C are not represented in current terrestrial biosphere models. Our results highlight that a better representation of C allocation should be incorporated into models.

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 CO₂, 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"

[Figure: see text].

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.

New generation geostationary satellite observations support seasonality in greenness of the Amazon evergreen forests

Assessing the seasonal patterns of the Amazon rainforests has been difficult because of the paucity of ground observations and persistent cloud cover over these forests obscuring optical remote sensing observations. Here, we use data from a new generation of geostationary satellites that carry the Advanced Baseline Imager (ABI) to study the Amazon canopy. ABI is similar to the widely used polar orbiting sensor, the Moderate Resolution Imaging Spectroradiometer (MODIS), but provides observations every 10-15 min. Our analysis of NDVI data collected over the Amazon during 2018-19 shows that ABI provides 21-35 times more cloud-free observations in a month than MODIS. The analyses show statistically significant changes in seasonality over 85% of Amazon forest pixels, an area about three times greater than previously reported using MODIS data. Though additional work is needed in converting the observed changes in seasonality into meaningful changes in canopy dynamics, our results highlight the potential of the new generation geostationary satellites to help us better understand tropical ecosystems, which has been a challenge with only polar orbiting satellites.

Biophysical impacts of Earth greening largely controlled by aerodynamic resistance

Satellite observations show widespread increasing trends of leaf area index (LAI), known as the Earth greening. However, the biophysical impacts of this greening on land surface temperature (LST) remain unclear. Here, we quantify the biophysical impacts of Earth greening on LST from 2000 to 2014 and disentangle the contributions of different factors using a physically based attribution model. We find that 93% of the global vegetated area shows negative sensitivity of LST to LAI increase at the annual scale, especially for semiarid woody vegetation. Further considering the LAI trends (P ≤ 0.1), 30% of the global vegetated area is cooled by these trends and 5% is warmed. Aerodynamic resistance is the dominant factor in controlling Earth greening's biophysical impacts: The increase in LAI produces a decrease in aerodynamic resistance, thereby favoring increased turbulent heat transfer between the land and the atmosphere, especially latent heat flux.

Spatial and temporal variations in global soil respiration and their relationships with climate and land cover

Soil respiration (R _s) represents the largest flux of CO₂ from terrestrial ecosystems to the atmosphere, but its spatial and temporal changes as well as the driving forces are not well understood. We derived a product of annual global R _s from 2000 to 2014 at 1 km by 1 km spatial resolution using remote sensing data and biome-specific statistical models. Different from the existing view that climate change dominated changes in R _s, we showed that land-cover change played a more important role in regulating R _s changes in temperate and boreal regions during 2000-2014. Significant changes in R _s occurred more frequently in areas with significant changes in short vegetation cover (i.e., all vegetation shorter than 5 m in height) than in areas with significant climate change. These results contribute to our understanding of global R _s patterns and highlight the importance of land-cover change in driving global and regional R _s changes.

Summer soil drying exacerbated by earlier spring greening of northern vegetation

Earlier vegetation greening under climate change raises evapotranspiration and thus lowers spring soil moisture, yet the extent and magnitude of this water deficit persistence into the following summer remain elusive. We provide observational evidence that increased foliage cover over the Northern Hemisphere, during 1982-2011, triggers an additional soil moisture deficit that is further carried over into summer. Climate model simulations independently support this and attribute the driving process to be larger increases in evapotranspiration than in precipitation. This extra soil drying is projected to amplify the frequency and intensity of summer heatwaves. Most feedbacks operate locally, except for a notable teleconnection where extra moisture transpired over Europe is transported to central Siberia. Model results illustrate that this teleconnection offsets Siberian soil moisture losses from local spring greening. Our results highlight that climate change adaptation planning must account for the extra summer water and heatwave stress inherited from warming-induced earlier greening.

Satellite-observed pantropical carbon dynamics

Changes in terrestrial tropical carbon stocks have an important role in the global carbon budget. However, current observational tools do not allow accurate and large-scale monitoring of the spatial distribution and dynamics of carbon stocks¹. Here, we used low-frequency L-band passive microwave observations to compute a direct and spatially explicit quantification of annual aboveground carbon (AGC) fluxes and show that the tropical net AGC budget was approximately in balance during 2010 to 2017, the net budget being composed of gross losses of -2.86 PgC yr^-1 offset by gross gains of -2.97 PgC yr^-1 between continents. Large interannual and spatial fluctuations of tropical AGC were quantified during the wet 2011 La Niña year and throughout the extreme dry and warm 2015-2016 El Niño episode. These interannual fluctuations, controlled predominantly by semiarid biomes, were shown to be closely related to independent global atmospheric CO₂ growth-rate anomalies (Pearson's r = 0.86), highlighting the pivotal role of tropical AGC in the global carbon budget.

Changes in timing of seasonal peak photosynthetic activity in northern ecosystems

Seasonality in photosynthetic activity is a critical component of seasonal carbon, water, and energy cycles in the Earth system. This characteristic is a consequence of plant's adaptive evolutionary processes to a given set of environmental conditions. Changing climate in northern lands (>30°N) alters the state of climatic constraints on plant growth, and therefore, changes in the seasonality and carbon accumulation are anticipated. However, how photosynthetic seasonality evolved to its current state, and what role climatic constraints and their variability played in this process and ultimately in carbon cycle is still poorly understood due to its complexity. Here, we take the "laws of minimum" as a basis and introduce a new framework where the timing (day of year) of peak photosynthetic activity (DOY_Pmax ) acts as a proxy for plant's adaptive state to climatic constraints on its growth. Our analyses confirm that spatial variations in DOY_Pmax reflect spatial gradients in climatic constraints as well as seasonal maximum and total productivity. We find a widespread warming-induced advance in DOY_Pmax (-1.66 ± 0.30 days/decade, p < 0.001) across northern lands, indicating a spatiotemporal dynamism of climatic constraints to plant growth. We show that the observed changes in DOY_Pmax are associated with an increase in total gross primary productivity through enhanced carbon assimilation early in the growing season, which leads to an earlier phase shift in land-atmosphere carbon fluxes and an increase in their amplitude. Such changes are expected to continue in the future based on our analysis of earth system model projections. Our study provides a simplified, yet realistic framework based on first principles for the complex mechanisms by which various climatic factors constrain plant growth in northern ecosystems.

Air temperature optima of vegetation productivity across global biomes

The global distribution of the optimum air temperature for ecosystem-level gross primary productivity ([Formula: see text]) is poorly understood, despite its importance for ecosystem carbon uptake under future warming. We provide empirical evidence for the existence of such an optimum, using measurements of in situ eddy covariance and satellite-derived proxies, and report its global distribution. [Formula: see text] is consistently lower than the physiological optimum temperature of leaf-level photosynthetic capacity, which typically exceeds 30 °C. The global average [Formula: see text] is estimated to be 23 ± 6 °C, with warmer regions having higher [Formula: see text] values than colder regions. In tropical forests in particular, [Formula: see text] is close to growing-season air temperature and is projected to fall below it under all scenarios of future climate, suggesting a limited safe operating space for these ecosystems under future warming.

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.

Coupling of ecosystem-scale plant water storage and leaf phenology observed by satellite

Plant water storage is fundamental to the functioning of terrestrial ecosystems by participating in plant metabolism, nutrient and sugar transport, and maintenance of the integrity of the hydraulic system of the plant. However, a global view of the size and dynamics of the water pools stored in plant tissues is still lacking. Here, we report global patterns of seasonal variations in ecosystem-scale plant water storage and their relationship with leaf phenology, based on space-borne measurements of L-band vegetation optical depth. We find that seasonal variations in plant water storage are highly synchronous with leaf phenology for the boreal and temperate forests, but asynchronous for the tropical woodlands, where the seasonal development of plant water storage lags behind leaf area by up to 180 days. Contrasting patterns of the time lag between plant water storage and terrestrial groundwater storage are also evident in these ecosystems. A comparison of the water cycle components in seasonally dry tropical woodlands highlights the buffering effect of plant water storage on the seasonal dynamics of water supply and demand. Our results offer insights into ecosystem-scale plant water relations globally and provide a basis for an improved parameterization of eco-hydrological and Earth system models.

Post-drought decline of the Amazon carbon sink

Amazon forests have experienced frequent and severe droughts in the past two decades. However, little is known about the large-scale legacy of droughts on carbon stocks and dynamics of forests. Using systematic sampling of forest structure measured by LiDAR waveforms from 2003 to 2008, here we show a significant loss of carbon over the entire Amazon basin at a rate of 0.3 ± 0.2 (95% CI) PgC yr⁻¹ after the 2005 mega-drought, which continued persistently over the next 3 years (2005–2008). The changes in forest structure, captured by average LiDAR forest height and converted to above ground biomass carbon density, show an average loss of 2.35 ± 1.80 MgC ha⁻¹ a year after (2006) in the epicenter of the drought. With more frequent droughts expected in future, forests of Amazon may lose their role as a robust sink of carbon, leading to a significant positive climate feedback and exacerbating warming trends.

Forests of the Amazon Basin have experienced frequent and severe droughts in recent years with significant impacts on their carbon cycling. Here, using satellite LiDAR samples from 2003 to 2008, the authors show the long-term legacy of these droughts with persistent loss of carbon stocks after the 2005 drought.

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 m²  m^-2  year^-1 , p < .01), followed by autumn (0.005 m²  m^-2  year^-1 , p < .01) and spring (0.003 m² m^-2  year^-1 , p < .01). We then quantified the contribution of elevating atmospheric CO₂ concentration (eCO₂ ), 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 eCO₂ provides the greatest contribution to increasing LAI trends in all seasons (0.003-0.007 m²  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 m²  m^-2  year^-1 and 0.0001-0.001 m²  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.

Estimation of leaf area index and its sunlit portion from DSCOVR EPIC data: Theoretical basis

This paper presents the theoretical basis of the algorithm designed for the generation of leaf area index and diurnal course of its sunlit portion from NASA's Earth Polychromatic Imaging Camera (EPIC) onboard NOAA's Deep Space Climate Observatory (DSCOVR). The Look-up-Table (LUT) approach implemented in the MODIS operational LAI/FPAR algorithm is adopted. The LUT, which is the heart of the approach, has been significantly modified. First, its parameterization incorporates the canopy hot spot phenomenon and recent advances in the theory of canopy spectral invariants. This allows more accurate decoupling of the structural and radiometric components of the measured Bidirectional Reflectance Factor (BRF), improves scaling properties of the LUT and consequently simplifies adjustments of the algorithm for data spatial resolution and spectral band compositions. Second, the stochastic radiative transfer equations are used to generate the LUT for all biome types. The equations naturally account for radiative effects of the three-dimensional canopy structure on the BRF and allow for an accurate discrimination between sunlit and shaded leaf areas. Third, the LUT entries are measurable, i.e., they can be independently derived from both below canopy measurements of the transmitted and above canopy measurements of reflected radiation fields. This feature makes possible direct validation of the LUT, facilitates identification of its deficiencies and development of refinements. Analyses of field data on canopy structure and leaf optics collected at 18 sites in the Hyytiälä forest in southern boreal zone in Finland and hyperspectral images acquired by the EO-1 Hyperion sensor support the theoretical basis.

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.

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.

Evaporative cooling over the Tibetan Plateau induced by vegetation growth

In the Arctic, climate warming enhances vegetation activity by extending the length of the growing season and intensifying maximum rates of productivity. In turn, increased vegetation productivity reduces albedo, which causes a positive feedback on temperature. Over the Tibetan Plateau (TP), regional vegetation greening has also been observed in response to recent warming. Here, we show that in contrast to arctic regions, increased growing season vegetation activity over the TP may have attenuated surface warming. This negative feedback on growing season vegetation temperature is attributed to enhanced evapotranspiration (ET). The extra energy available at the surface, which results from lower albedo, is efficiently dissipated by evaporative cooling. The net effect is a decrease in daily maximum temperature and the diurnal temperature range, which is supported by statistical analyses of in situ observations and by decomposition of the surface energy budget. A daytime cooling effect from increased vegetation activity is also modeled from a set of regional weather research and forecasting (WRF) mesoscale model simulations, but with a magnitude smaller than observed, likely because the WRF model simulates a weaker ET enhancement. Our results suggest that actions to restore native grasslands in degraded areas, roughly one-third of the plateau, will both facilitate a sustainable ecological development in this region and have local climate cobenefits. More accurate simulations of the biophysical coupling between the land surface and the atmosphere are needed to help understand regional climate change over the TP, and possible larger scale feedbacks between climate in the TP and the Asian monsoon system.

Detection and attribution of vegetation greening trend in China over the last 30 years

The reliable detection and attribution of changes in vegetation growth is a prerequisite for the development of strategies for the sustainable management of ecosystems. This is an extraordinary challenge. To our knowledge, this study is the first to comprehensively detect and attribute a greening trend in China over the last three decades. We use three different satellite-derived Leaf Area Index (LAI) datasets for detection as well as five different process-based ecosystem models for attribution. Rising atmospheric CO2 concentration and nitrogen deposition are identified as the most likely causes of the greening trend in China, explaining 85% and 41% of the average growing-season LAI trend (LAIGS) estimated by satellite datasets (average trend of 0.0070 yr(-1), ranging from 0.0035 yr(-1) to 0.0127 yr(-1)), respectively. The contribution of nitrogen deposition is more clearly seen in southern China than in the north of the country. Models disagree about the contribution of climate change alone to the trend in LAIGS at the country scale (one model shows a significant increasing trend, whereas two others show significant decreasing trends). However, the models generally agree on the negative impacts of climate change in north China and Inner Mongolia and the positive impact in the Qinghai-Xizang plateau. Provincial forest area change tends to be significantly correlated with the trend of LAIGS (P < 0.05), and marginally significantly (P = 0.07) correlated with the residual of LAIGS trend, calculated as the trend observed by satellite minus that estimated by models through considering the effects of climate change, rising CO2 concentration and nitrogen deposition, across different provinces. This result highlights the important role of China's afforestation program in explaining the spatial patterns of trend in vegetation growth.

Seasonally different response of photosynthetic activity to daytime and night-time warming in the Northern Hemisphere

Over the last century the Northern Hemisphere has experienced rapid climate warming, but this warming has not been evenly distributed seasonally, as well as diurnally. The implications of such seasonal and diurnal heterogeneous warming on regional and global vegetation photosynthetic activity, however, are still poorly understood. Here, we investigated for different seasons how photosynthetic activity of vegetation correlates with changes in seasonal daytime and night-time temperature across the Northern Hemisphere (>30°N), using Normalized Difference Vegetation Index (NDVI) data from 1982 to 2011 obtained from the Advanced Very High Resolution Radiometer (AVHRR). Our analysis revealed some striking seasonal differences in the response of NDVI to changes in day- vs. night-time temperatures. For instance, while higher daytime temperature (Tmax) is generally associated with higher NDVI values across the boreal zone, the area exhibiting a statistically significant positive correlation between Tmax and NDVI is much larger in spring (41% of area in boreal zone--total area 12.6×10(6) km2) than in summer and autumn (14% and 9%, respectively). In contrast to the predominantly positive response of boreal ecosystems to changes in Tmax, increases in Tmax tended to negatively influence vegetation growth in temperate dry regions, particularly during summer. Changes in night-time temperature (Tmin) correlated negatively with autumnal NDVI in most of the Northern Hemisphere, but had a positive effect on spring and summer NDVI in most temperate regions (e.g., Central North America and Central Asia). Such divergent covariance between the photosynthetic activity of Northern Hemispheric vegetation and day- and night-time temperature changes among different seasons and climate zones suggests a changing dominance of ecophysiological processes across time and space. Understanding the seasonally different responses of vegetation photosynthetic activity to diurnal temperature changes, which have not been captured by current land surface models, is important for improving the performance of next generation regional and global coupled vegetation-climate models.

Contribution of semi-arid ecosystems to interannual variability of the global carbon cycle

The land and ocean act as a sink for fossil-fuel emissions, thereby slowing the rise of atmospheric carbon dioxide concentrations. Although the uptake of carbon by oceanic and terrestrial processes has kept pace with accelerating carbon dioxide emissions until now, atmospheric carbon dioxide concentrations exhibit a large variability on interannual timescales, considered to be driven primarily by terrestrial ecosystem processes dominated by tropical rainforests. We use a terrestrial biogeochemical model, atmospheric carbon dioxide inversion and global carbon budget accounting methods to investigate the evolution of the terrestrial carbon sink over the past 30 years, with a focus on the underlying mechanisms responsible for the exceptionally large land carbon sink reported in 2011 (ref. 2). Here we show that our three terrestrial carbon sink estimates are in good agreement and support the finding of a 2011 record land carbon sink. Surprisingly, we find that the global carbon sink anomaly was driven by growth of semi-arid vegetation in the Southern Hemisphere, with almost 60 per cent of carbon uptake attributed to Australian ecosystems, where prevalent La Niña conditions caused up to six consecutive seasons of increased precipitation. In addition, since 1981, a six per cent expansion of vegetation cover over Australia was associated with a fourfold increase in the sensitivity of continental net carbon uptake to precipitation. Our findings suggest that the higher turnover rates of carbon pools in semi-arid biomes are an increasingly important driver of global carbon cycle inter-annual variability and that tropical rainforests may become less relevant drivers in the future. More research is needed to identify to what extent the carbon stocks accumulated during wet years are vulnerable to rapid decomposition or loss through fire in subsequent years.

A two-fold increase of carbon cycle sensitivity to tropical temperature variations

Earth system models project that the tropical land carbon sink will decrease in size in response to an increase in warming and drought during this century, probably causing a positive climate feedback. But available data are too limited at present to test the predicted changes in the tropical carbon balance in response to climate change. Long-term atmospheric carbon dioxide data provide a global record that integrates the interannual variability of the global carbon balance. Multiple lines of evidence demonstrate that most of this variability originates in the terrestrial biosphere. In particular, the year-to-year variations in the atmospheric carbon dioxide growth rate (CGR) are thought to be the result of fluctuations in the carbon fluxes of tropical land areas. Recently, the response of CGR to tropical climate interannual variability was used to put a constraint on the sensitivity of tropical land carbon to climate change. Here we use the long-term CGR record from Mauna Loa and the South Pole to show that the sensitivity of CGR to tropical temperature interannual variability has increased by a factor of 1.9 ± 0.3 in the past five decades. We find that this sensitivity was greater when tropical land regions experienced drier conditions. This suggests that the sensitivity of CGR to interannual temperature variations is regulated by moisture conditions, even though the direct correlation between CGR and tropical precipitation is weak. We also find that present terrestrial carbon cycle models do not capture the observed enhancement in CGR sensitivity in the past five decades. More realistic model predictions of future carbon cycle and climate feedbacks require a better understanding of the processes driving the response of tropical ecosystems to drought and warming.

Large-scale variations in the vegetation growing season and annual cycle of atmospheric CO2 at high northern latitudes from 1950 to 2011

We combine satellite and ground observations during 1950-2011 to study the long-term links between multiple climate (air temperature and cryospheric dynamics) and vegetation (greenness and atmospheric CO(2) concentrations) indicators of the growing season of northern ecosystems (>45°N) and their connection with the carbon cycle. During the last three decades, the thermal potential growing season has lengthened by about 10.5 days (P < 0.01, 1982-2011), which is unprecedented in the context of the past 60 years. The overall lengthening has been stronger and more significant in Eurasia (12.6 days, P < 0.01) than North America (6.2 days, P > 0.05). The photosynthetic growing season has closely tracked the pace of warming and extension of the potential growing season in spring, but not in autumn when factors such as light and moisture limitation may constrain photosynthesis. The autumnal extension of the photosynthetic growing season since 1982 appears to be about half that of the thermal potential growing season, yielding a smaller lengthening of the photosynthetic growing season (6.7 days at the circumpolar scale, P < 0.01). Nevertheless, when integrated over the growing season, photosynthetic activity has closely followed the interannual variations and warming trend in cumulative growing season temperatures. This lengthening and intensification of the photosynthetic growing season, manifested principally over Eurasia rather than North America, is associated with a long-term increase (22.2% since 1972, P < 0.01) in the amplitude of the CO(2) annual cycle at northern latitudes. The springtime extension of the photosynthetic and potential growing seasons has apparently stimulated earlier and stronger net CO(2) uptake by northern ecosystems, while the autumnal extension is associated with an earlier net release of CO(2) to the atmosphere. These contrasting responses may be critical in determining the impact of continued warming on northern terrestrial ecosystems and the carbon cycle.

Evaluation of terrestrial carbon cycle models for their response to climate variability and to CO2 trends

The purpose of this study was to evaluate 10 process-based terrestrial biosphere models that were used for the IPCC fifth Assessment Report. The simulated gross primary productivity (GPP) is compared with flux-tower-based estimates by Jung et al. [Journal of Geophysical Research 116 (2011) G00J07] (JU11). The net primary productivity (NPP) apparent sensitivity to climate variability and atmospheric CO2 trends is diagnosed from each model output, using statistical functions. The temperature sensitivity is compared against ecosystem field warming experiments results. The CO2 sensitivity of NPP is compared to the results from four Free-Air CO2 Enrichment (FACE) experiments. The simulated global net biome productivity (NBP) is compared with the residual land sink (RLS) of the global carbon budget from Friedlingstein et al. [Nature Geoscience 3 (2010) 811] (FR10). We found that models produce a higher GPP (133 ± 15 Pg C yr(-1) ) than JU11 (118 ± 6 Pg C yr(-1) ). In response to rising atmospheric CO2 concentration, modeled NPP increases on average by 16% (5-20%) per 100 ppm, a slightly larger apparent sensitivity of NPP to CO2 than that measured at the FACE experiment locations (13% per 100 ppm). Global NBP differs markedly among individual models, although the mean value of 2.0 ± 0.8 Pg C yr(-1) is remarkably close to the mean value of RLS (2.1 ± 1.2 Pg C yr(-1) ). The interannual variability in modeled NBP is significantly correlated with that of RLS for the period 1980-2009. Both model-to-model and interannual variation in model GPP is larger than that in model NBP due to the strong coupling causing a positive correlation between ecosystem respiration and GPP in the model. The average linear regression slope of global NBP vs. temperature across the 10 models is -3.0 ± 1.5 Pg C yr(-1) °C(-1) , within the uncertainty of what derived from RLS (-3.9 ± 1.1 Pg C yr(-1) °C(-1) ). However, 9 of 10 models overestimate the regression slope of NBP vs. precipitation, compared with the slope of the observed RLS vs. precipitation. With most models lacking processes that control GPP and NBP in addition to CO2 and climate, the agreement between modeled and observation-based GPP and NBP can be fortuitous. Carbon-nitrogen interactions (only separable in one model) significantly influence the simulated response of carbon cycle to temperature and atmospheric CO2 concentration, suggesting that nutrients limitations should be included in the next generation of terrestrial biosphere models.

Changes in satellite-derived spring vegetation green-up date and its linkage to climate in China from 1982 to 2010: a multimethod analysis

The change in spring phenology is recognized to exert a major influence on carbon balance dynamics in temperate ecosystems. Over the past several decades, several studies focused on shifts in spring phenology; however, large uncertainties still exist, and one understudied source could be the method implemented in retrieving satellite-derived spring phenology. To account for this potential uncertainty, we conducted a multimethod investigation to quantify changes in vegetation green-up date from 1982 to 2010 over temperate China, and to characterize climatic controls on spring phenology. Over temperate China, the five methods estimated that the vegetation green-up onset date advanced, on average, at a rate of 1.3 ± 0.6 days per decade (ranging from 0.4 to 1.9 days per decade) over the last 29 years. Moreover, the sign of the trends in vegetation green-up date derived from the five methods were broadly consistent spatially and for different vegetation types, but with large differences in the magnitude of the trend. The large intermethod variance was notably observed in arid and semiarid vegetation types. Our results also showed that change in vegetation green-up date is more closely correlated with temperature than with precipitation. However, the temperature sensitivity of spring vegetation green-up date became higher as precipitation increased, implying that precipitation is an important regulator of the response of vegetation spring phenology to change in temperature. This intricate linkage between spring phenology and precipitation must be taken into account in current phenological models which are mostly driven by temperature.

Surface urban heat island across 419 global big cities

Urban heat island is among the most evident aspects of human impacts on the earth system. Here we assess the diurnal and seasonal variation of surface urban heat island intensity (SUHII) defined as the surface temperature difference between urban area and suburban area measured from the MODIS. Differences in SUHII are analyzed across 419 global big cities, and we assess several potential biophysical and socio-economic driving factors. Across the big cities, we show that the average annual daytime SUHII (1.5 ± 1.2 °C) is higher than the annual nighttime SUHII (1.1 ± 0.5 °C) (P < 0.001). But no correlation is found between daytime and nighttime SUHII across big cities (P = 0.84), suggesting different driving mechanisms between day and night. The distribution of nighttime SUHII correlates positively with the difference in albedo and nighttime light between urban area and suburban area, while the distribution of daytime SUHII correlates negatively across cities with the difference of vegetation cover and activity between urban and suburban areas. Our results emphasize the key role of vegetation feedbacks in attenuating SUHII of big cities during the day, in particular during the growing season, further highlighting that increasing urban vegetation cover could be one effective way to mitigate the urban heat island effect.

Comment on "Drought-induced reduction in global terrestrial net primary production from 2000 through 2009"

Zhao and Running (Reports, 20 August 2010, p. 940) reported a reduction in global terrestrial net primary production (NPP) from 2000 through 2009. We argue that the small trends, regional patterns, and interannual variations that they describe are artifacts of their NPP model. Satellite observations of vegetation activity show no statistically significant changes in more than 85% of the vegetated lands south of 70°N during the same 2000 to 2009 period.

Large seasonal swings in leaf area of Amazon rainforests

Despite early speculation to the contrary, all tropical forests studied to date display seasonal variations in the presence of new leaves, flowers, and fruits. Past studies were focused on the timing of phenological events and their cues but not on the accompanying changes in leaf area that regulate vegetation-atmosphere exchanges of energy, momentum, and mass. Here we report, from analysis of 5 years of recent satellite data, seasonal swings in green leaf area of approximately 25% in a majority of the Amazon rainforests. This seasonal cycle is timed to the seasonality of solar radiation in a manner that is suggestive of anticipatory and opportunistic patterns of net leaf flushing during the early to mid part of the light-rich dry season and net leaf abscission during the cloudy wet season. These seasonal swings in leaf area may be critical to initiation of the transition from dry to wet season, seasonal carbon balance between photosynthetic gains and respiratory losses, and litterfall nutrient cycling in moist tropical forests.

Evidence for a significant urbanization effect on climate in China

China has experienced rapid urbanization and dramatic economic growth since its reform process started in late 1978. In this article, we present evidence for a significant urbanization effect on climate based on analysis of impacts of land-use changes on surface temperature in southeast China, where rapid urbanization has occurred. Our estimated warming of mean surface temperature of 0.05 degrees C per decade attributable to urbanization is much larger than previous estimates for other periods and locations. The spatial pattern and magnitude of our estimate are consistent with those of urbanization characterized by changes in the percentage of urban population and in satellite-measured greenness.

Climate-driven increases in global terrestrial net primary production from 1982 to 1999

Recent climatic changes have enhanced plant growth in northern mid-latitudes and high latitudes. However, a comprehensive analysis of the impact of global climatic changes on vegetation productivity has not before been expressed in the context of variable limiting factors to plant growth. We present a global investigation of vegetation responses to climatic changes by analyzing 18 years (1982 to 1999) of both climatic data and satellite observations of vegetation activity. Our results indicate that global changes in climate have eased several critical climatic constraints to plant growth, such that net primary production increased 6% (3.4 petagrams of carbon over 18 years) globally. The largest increase was in tropical ecosystems. Amazon rain forests accounted for 42% of the global increase in net primary production, owing mainly to decreased cloud cover and the resulting increase in solar radiation.

Climatic control of the high-latitude vegetation greening trend and Pinatubo effect

A biogeochemical model of vegetation using observed climate data predicts the high northern latitude greening trend over the past two decades observed by satellites and a marked setback in this trend after the Mount Pinatubo volcano eruption in 1991. The observed trend toward earlier spring budburst and increased maximum leaf area is produced by the model as a consequence of biogeochemical vegetation responses mainly to changes in temperature. The post-Pinatubo decline in vegetation in 1992-1993 is apparent as the effect of temporary cooling caused by the eruption. High-latitude CO(2) uptake during these years is predicted as a consequence of the differential response of heterotrophic respiration and net primary production.

A large carbon sink in the woody biomass of Northern forests

The terrestrial carbon sink, as of yet unidentified, represents 15-30% of annual global emissions of carbon from fossil fuels and industrial activities. Some of the missing carbon is sequestered in vegetation biomass and, under the Kyoto Protocol of the United Nations Framework Convention on Climate Change, industrialized nations can use certain forest biomass sinks to meet their greenhouse gas emissions reduction commitments. Therefore, we analyzed 19 years of data from remote-sensing spacecraft and forest inventories to identify the size and location of such sinks. The results, which cover the years 1981-1999, reveal a picture of biomass carbon gains in Eurasian boreal and North American temperate forests and losses in some Canadian boreal forests. For the 1.42 billion hectares of Northern forests, roughly above the 30th parallel, we estimate the biomass sink to be 0.68 +/- 0.34 billion tons carbon per year, of which nearly 70% is in Eurasia, in proportion to its forest area and in disproportion to its biomass carbon pool. The relatively high spatial resolution of these estimates permits direct validation with ground data and contributes to a monitoring program of forest biomass sinks under the Kyoto protocol.

Higher northern latitude normalized difference vegetation index and growing season trends from 1982 to 1999

Normalized difference vegetation index data from the polar-orbiting National Oceanic and Atmospheric Administration meteorological satellites from 1982 to 1999 show significant variations in photosynthetic activity and growing season length at latitudes above 35 degrees N. Two distinct periods of increasing plant growth are apparent: 1982-1991 and 1992-1999, separated by a reduction from 1991 to 1992 associated with global cooling resulting from the volcanic eruption of Mt. Pinatubo in June 1991. The average May to September normalized difference vegetation index from 45 degrees N to 75 degrees N increased by 9% from 1982 to 1991, decreased by 5% from 1991 to 1992, and increased by 8% from 1992 to 1999. Variations in the normalized difference vegetation index were associated with variations in the start of the growing season of -5.6, +3.9, and -1.7 days respectively, for the three time periods. Our results support surface temperature increases within the same period at higher northern latitudes where temperature limits plant growth.