globalchange  > 气候变化事实与影响
DOI: doi:10.1038/nclimate2533
论文题名:
Three decades of multi-dimensional change in global leaf phenology
作者: Robert Buitenwerf
刊名: Nature Climate Change
ISSN: 1758-997X
EISSN: 1758-7117
出版年: 2015-03-02
卷: Volume:5, 页码:Pages:364;368 (2015)
语种: 英语
英文关键词: Phenology
英文摘要:

Changes in the phenology of vegetation activity may accelerate or dampen rates of climate change by altering energy exchanges between the land surface and the atmosphere1, 2 and can threaten species with synchronized life cycles3, 4, 5. Current knowledge of long-term changes in vegetation activity is regional6, 7, 8, or restricted to highly integrated measures of change such as net primary productivity9, 10, 11, 12, 13, which mask details that are relevant for Earth system dynamics. Such details can be revealed by measuring changes in the phenology of vegetation activity. Here we undertake a comprehensive global assessment of changes in vegetation phenology. We show that the phenology of vegetation activity changed severely (by more than 2 standard deviations in one or more dimensions of phenological change) on 54% of the global land surface between 1981 and 2012. Our analysis confirms previously detected changes in the boreal and northern temperate regions6, 7, 8. The adverse consequences of these northern phenological shifts for land-surface–climate feedbacks1, ecosystems4 and species3 are well known. Our study reveals equally severe phenological changes in the southern hemisphere, where consequences for the energy budget and the likelihood of phenological mismatches are unknown. Our analysis provides a sensitive and direct measurement of ecosystem functioning, making it useful both for monitoring change and for testing the reliability of early warning signals of change14.

Recent climate change has shifted species distributions15, 16 and leaf phenology17, 18 around the world, leading to mismatches in previously synchronized phenological cycles3, 4. Such mismatches greatly increase the risk of extinction for affected species, and ongoing climatic and phenological change is expected to further increase this risk5. Despite documenting and predicting effects of climate change on many organisms, these previous studies do not provide an easy way of inferring how widespread such changes are or where they are most severe. In addition to being a symptom of climate change, vegetation change also feeds back to the climate system by forcing rates of energy exchange between the land surface and the atmosphere. Changes in the vigour and timing of vegetation activity can therefore accelerate or slow down rates of climate change1. Yet, the extent to which changes in vegetation phenology will impact the climate system by modifying albedo, transpiration, partitioning between latent and sensible heat in the atmosphere, and cloud formation, has been identified as a major source of uncertainty in climate change projections2, 19.

To quantify changes in global vegetation activity, previous studies have used remotely sensed data to quantify changes in primary productivity9, 10, 11, 12. These studies have indicated an overall increase in net primary productivity (NPP) during the 1980s and 1990s (ref. 9), whereas evidence for a decrease during the 2000s (ref. 11) has been debated20. Although estimating NPP is important for describing carbon sequestration, it is a highly integrated metric that masks important details of the nature of change. For example, it provides no information on the likelihood of phenological mismatches and limited information on consequences for the land-surface energy budget. Consequently, constant NPP does not guarantee that vegetation is not responding to changing climates and increased atmospheric CO2 in ways that affect the functioning of the Earth system. To quantify intra-annual shifts in the timing and vigour of vegetation activity, remotely sensed absorption of photosynthetically active radiation by the land surface can be used to directly infer photosynthetic activity.

We present a global analysis of change in the seasonal pattern of photosynthetically active radiation absorbed by the land surface as measured by the normalized difference vegetation index (NDVI). We improve on previous analyses that have used NDVI to infer phenological change in two important ways. First, previous studies on long-term changes in leaf phenology have been regional6, 7, 8. We analyse the GIMMS3g data, a global record from 1981 to 2012, at 0.083° and 15-day resolution. Second, a problem that has prevented global analyses of phenology is that the information content of phenological metrics is not universal. For example, the onset of the growing season is an informative metric in deciduous forests, but less useful in evergreen forests. We use an improved method to estimate 21 ecologically interpretable metrics of the phenological cycle from the data (Fig. 1) and evaluate the magnitude of change within 83 phenologically similar zones (hereafter called phenomes, shown in Supplementary Fig. 1) to account for the fact that the information content of these metrics differs between ecosystem types. These phenomes were identified using a cluster analysis of the phenological data. To obtain a spatially comparable measure of change, the change per pixel was scaled by the variance of change for the phenome to which it was assigned. Change is therefore reported in standard deviations (s.d.), which can be interpreted as a measure of the severity of change.

Figure 1: Calculated phenological metrics.
Calculated phenological metrics.

NDVI data for a single phenological year of an illustrative pixel, showing the 21 phenological metrics used in this study. Each labelled point represents a date and an associated NDVI value. The integral of the curve is also calculated. The fitted spline (f), first and second derivatives (f′,f”) were used to calculate phenological metrics. gsl, growing season length.

Data.

We used the GIMMS3g data set (downloaded from http://ecocast.arc.nasa.gov/data/pub/gimms/3g), which provides Advanced Very High Resolution Radiometer (AVHRR) NDVI data. The data are 15-day maximum-value composites from 1981 to 2012 at 0.083° resolution. Compositing reduces atmospheric effects (clouds, aerosols) and data have further been processed to reduce effects of navigation errors, major volcanic eruptions and orbital drift of older satellites, and have been subjected to several sensor calibration steps. Finally, a rigorous analysis was carried out to resolve remaining discrepancies between data from different generations of AVHRR sensors30. Data quality scores were provided for each observation.

Defining phenological metrics.

Data for each pixel were smoothed using a cubic spline function, which was weighted by data quality scores. A spline, unlike a parametric function, maintains a high fidelity to the data while retaining continuous first and second derivatives. We designed a two-step method to flexibly define the start of the time series for each pixel. First, the average Julian day of minimum NDVI (trough day) was calculated for each pixel using the 31-year time series. In the second step, the exact trough day for each phenological year was determined within a 180-day window (90 days on either side) around the 31-year mean trough day.

The period between two consecutive trough days constitutes a phenological year. For each phenological year 20 additional metrics were extracted (shown in Fig. 1). Instead of defining the ‘start of the season’ and ‘end of season’ as a single date, we calculated three dates for start and end of season, which improved our ability to detect change. For example, in some deciduous forests, invasive understorey plants extend the growing season by retaining leaves longer than native forest trees31. Such change is difficult to describe with a single metric—that is, a large decrease in NDVI still occurs when the trees drop their leaves while the understorey plants will cause low NDVI to persist for longer. ‘Leaf-on 1’, ‘leaf-on 2’, and ‘leaf-on 3’ represent the first increase, the fastest increase, and end of increasing (that is, the start of peak) photosynthetic activity, respectively. These dates can be interpreted as leaf emergence, rapid leaf expansion, and attainment of a full canopy. Similarly, three leaf-off dates were calculated and the number of days between leaf-on and leaf-off dates represented three measures of growing season length: ‘gsl-long’, ‘gsl’, and ‘gsl-peak’.

We calculate the three leaf-on and leaf-off dates using the first and second derivatives of the spline function. For smooth functions—for example, a sine—the maximum of the second derivative indicates when the function starts going concave up, and could hence be used to define the start of the growing season. However, NDVI data is not always smooth and sometimes produces a spline with multiple peaks within a phenological year. Instead of sacrificing accuracy by smoothing the data further, we centred and scaled both the first and second derivatives, summed them, and calculated the days at which this sum is maximal and minimal (shown in Fig. 1). This method selects for the maximum (or minimum) in the second derivative that is closest to the maximum (or minimum) in the first derivative and provided a robust way of defining the three leaf-on and leaf-off metrics.

This procedure yielded a 31-year time series of 21 phenological metrics for 2,075,445 pixels.

Phenological change over time.

For every pixel, change for each of the 21 metrics was defined as the difference between averages of the beginning and the end of the time series—that is, between 1981–1990 and 2003–2012. Trend analysis was avoided as there was no reason to assume a priori that change should be linear or monotonic. Comparing ten-year averages buffers against effects of anomalous climatic events in a single year on the change estimate.

Comparing change spatially.

As explained in the text, we scaled the change of each pixel by the variance of phenologically similar regions to obtain a spatially comparable measure of change. Because existing land cover classifications based on the biome concept did not adequately summarize phenological variation, we used a model-based clustering method to identify the optimal number of phenologically similar regions (phenomes). The procedure resulted in 83 phenomes (clustering details are given in the Supplementary Information).

Summarizing change.

The final step in our analysis was to summarize change globally, by grouping phenomes with similar type and magnitude of change into nine syndromes (Fig. 3a) using unsupervised hierarchical clustering (details are given in the Supplementary Information). For each syndrome, average change in s.d. was plotted along two axes, one representing change in NDVI, the other representing phenological change (that is, change in dates or number of days; Fig. 3b). The full workflow of our analysis is summarized in Supplementary Fig. 4.

  1. Bonan, G. B. Forests and climate change: Forcings, feedbacks, and the climate benefits of forests. Science 320, 14441449 (2008).
  2. Myhre, G. et al. in Climate Change 2013: The Physical Science Basis (eds Stocker, T. F. et al.) 659740 (IPCC, Cambridge Univ. Press, 2013).
  3. Both, C., Bouwhuis, S., Lessells, C. M. & Visser, M. E. Climate change and population declines in a long-distance migratory bird. Nature 441, 8183 (2006).
  4. Post, E. et al. Ecological dynamics across the Arctic associated with recent climate change. Science 325, 13551358 (2009).
  5. Thomas, C. D. et al. Extinction risk from climate change. Nature 427, 145148 (2004).
  6. Zhu, W. et al. Extension of the growing season due to delayed autumn over mid and high latitudes in North America during 1982–2006. Glob. Ecol. Biogeogr. 21, 260271 (2012). URL:
http://www.nature.com/nclimate/journal/v5/n4/full/nclimate2533.html
Citation statistics:
资源类型: 期刊论文
标识符: http://119.78.100.158/handle/2HF3EXSE/4825
Appears in Collections:气候变化事实与影响
科学计划与规划
气候变化与战略

Files in This Item:
File Name/ File Size Content Type Version Access License
nclimate2533.pdf(1077KB)期刊论文作者接受稿开放获取View Download

Recommended Citation:
Robert Buitenwerf. Three decades of multi-dimensional change in global leaf phenology[J]. Nature Climate Change,2015-03-02,Volume:5:Pages:364;368 (2015).
Service
Recommend this item
Sava as my favorate item
Show this item's statistics
Export Endnote File
Google Scholar
Similar articles in Google Scholar
[Robert Buitenwerf]'s Articles
百度学术
Similar articles in Baidu Scholar
[Robert Buitenwerf]'s Articles
CSDL cross search
Similar articles in CSDL Cross Search
[Robert Buitenwerf]‘s Articles
Related Copyright Policies
Null
收藏/分享
文件名: nclimate2533.pdf
格式: Adobe PDF
此文件暂不支持浏览
所有评论 (0)
暂无评论
 

Items in IR are protected by copyright, with all rights reserved, unless otherwise indicated.