Comparison of phenology trends by land cover class: a case study in the Great Basin, USA

Direct impacts of human land use and indirect impacts of anthropogenic climate change may alter land cover and associated ecosystem function, affecting ecological goods and services. Considerable work has been done to identify long‐term global trends in vegetation greenness, which is associated with primary productivity, using remote sensing. Trend analysis of satellite observations is subject to error, and ecosystem change can be confused with interannual variability. However, the relative trends of land cover classes may hold clues about differential ecosystem response to environmental forcing. Our aim was to identify phenological variability and 10‐year trends for the major land cover classes in the Great Basin. This case study involved two steps: a regional, phenology‐based land cover classification and an identification of phenological variability and 10‐year trends stratified by land cover class. The analysis used a 10‐year time series of Advanced Very High Resolution Radiometer satellite data to assess regional scale land cover variability and identify change. The phenology‐based regional classification was more detailed and accurate than national or global products. Phenological variability over the 10‐year period was high, with substantial shifts in timing of start of season of up to 9 weeks. The mean long‐term trends of montane land cover classes were significantly different from valley land cover classes due to a poor response of montane shrubland and pinyon‐juniper woodland to the early 1990s drought. The differential response during the 1990s suggests that valley ecosystems may be more resilient and montane ecosystems more susceptible to prolonged drought. This type of regional‐scale land cover analysis is necessary to characterize current patterns of land cover phenology, distinguish between anthropogenically driven land cover change and interannual variability, and identify ecosystems potentially susceptible to regional and global change.     

Variability and evolution of global land surface phenology over the past three decades (1982–2012)

Monitoring land surface phenology (LSP) is important for understanding both the responses and feedbacks of ecosystems to the climate system, and for representing these accurately in terrestrial biosphere models. Moreover, by shedding light on phenological trends at a variety of scales, LSP provides the potential to fill the gap between traditional phenological (field) observations and the large‐scale view of global models. In this study, we review and evaluate the variability and evolution of satellite‐derived growing season length (GSL) globally and over the past three decades. We used the longest continuous record of Normalized Difference Vegetation Index data available to date at global scale to derive LSP metrics consistently over all vegetated land areas and for the period 1982–2012. We tested GSL, start‐ and end‐of‐season metrics (SOS and EOS, respectively) for linear trends as well as for significant trend shifts over the study period. We evaluated trends using global environmental stratification information in place of commonly used land cover maps to avoid circular findings. Our results confirmed an average lengthening of the growing season globally during 1982–2012 – averaging 0.22–0.34 days yr^?1, but with spatially heterogeneous trends. About 13–19% of global land areas displayed significant GSL change, and over 30% of trends occurred in the boreal/alpine biome of the Northern Hemisphere, which showed diverging GSL evolution over the past three decades. Within this biome, the ‘Cold and Mesic’ environmental zone appeared as an LSP change hotspot. We also examined the relative contribution of SOS and EOS to the overall changes, finding that EOS trends were generally stronger and more prevalent than SOS trends. These findings constitute a step towards the identification of large‐scale phenological drivers of vegetated land surfaces, necessary for improving phenological representation in terrestrial biosphere models.     

Satellite chlorophyll fluorescence measurements reveal large‐scale decoupling of photosynthesis and greenness dynamics in boreal evergreen forests

Mid‐to‐high latitude forests play an important role in the terrestrial carbon cycle, but the representation of photosynthesis in boreal forests by current modelling and observational methods is still challenging. In particular, the applicability of existing satellite‐based proxies of greenness to indicate photosynthetic activity is hindered by small annual changes in green biomass of the often evergreen tree population and by the confounding effects of background materials such as snow. As an alternative, satellite measurements of sun‐induced chlorophyll fluorescence (SIF) can be used as a direct proxy of photosynthetic activity. In this study, the start and end of the photosynthetically active season of the main boreal forests are analysed using spaceborne SIF measurements retrieved from the GOME‐2 instrument and compared to that of green biomass, proxied by vegetation indices including the Enhanced Vegetation Index (EVI) derived from MODIS data. We find that photosynthesis and greenness show a similar seasonality in deciduous forests. In high‐latitude evergreen needleleaf forests, however, the length of the photosynthetically active period indicated by SIF is up to 6 weeks longer than the green biomass changing period proxied by EVI, with SIF showing a start‐of‐season of approximately 1 month earlier than EVI. On average, the photosynthetic spring recovery as signalled by SIF occurs as soon as air temperatures exceed the freezing point (2–3 °C) and when the snow on the ground has not yet completely melted. These findings are supported by model data of gross primary production and a number of other studies which evaluated in situ observations of CO₂ fluxes, meteorology and the physiological state of the needles. Our results demonstrate the sensitivity of space‐based SIF measurements to light‐use efficiency of boreal forests and their potential for an unbiased detection of photosynthetic activity even under the challenging conditions interposed by evergreen boreal ecosystems.     

Terrestrial biosphere models need better representation of vegetation phenology: results from the North American Carbon Program Site Synthesis

Phenology, by controlling the seasonal activity of vegetation on the land surface, plays a fundamental role in regulating photosynthesis and other ecosystem processes, as well as competitive interactions and feedbacks to the climate system. We conducted an analysis to evaluate the representation of phenology, and the associated seasonality of ecosystem‐scale CO₂ exchange, in 14 models participating in the North American Carbon Program Site Synthesis. Model predictions were evaluated using long‐term measurements (emphasizing the period 2000–2006) from 10 forested sites within the AmeriFlux and Fluxnet‐Canada networks. In deciduous forests, almost all models consistently predicted that the growing season started earlier, and ended later, than was actually observed; biases of 2 weeks or more were typical. For these sites, most models were also unable to explain more than a small fraction of the observed interannual variability in phenological transition dates. Finally, for deciduous forests, misrepresentation of the seasonal cycle resulted in over‐prediction of gross ecosystem photosynthesis by +160 ± 145 g C m^?2 yr^?1 during the spring transition period and +75 ± 130 g C m^?2 yr^?1 during the autumn transition period (13% and 8% annual productivity, respectively) compensating for the tendency of most models to under‐predict the magnitude of peak summertime photosynthetic rates. Models did a better job of predicting the seasonality of CO₂ exchange for evergreen forests. These results highlight the need for improved understanding of the environmental controls on vegetation phenology and incorporation of this knowledge into better phenological models. Existing models are unlikely to predict future responses of phenology to climate change accurately and therefore will misrepresent the seasonality and interannual variability of key biosphere–atmosphere feedbacks and interactions in coupled global climate models.     

Land Surface Phenology in the Tropics: The Role of Climate and Topography in a Snow-Free Mountain

Leaf phenology represents a major temporal component of ecosystem functioning, and understanding the drivers of seasonal variation in phenology is essential to understand plant responses to climate change. We assessed the patterns and drivers of land surface phenology, a proxy for leafing phenology, for the meridional Espinhaço Range, a South American tropical mountain comprising a mosaic of savannas, dry woodlands, montane vegetation and moist forests. We used a 14-year time series of MODIS/NDVI satellite images, acquired between 2001 and 2015, and extracted phenological indicators using the TIMESAT algorithm. We obtained precipitation data from the Tropical Rainfall Measuring Mission, land surface temperature from the MODIS MOD11A2 product, and cloud cover frequency from the MODIS MOD09GA product. We also calculated the topographic wetness index and simulated clear-sky radiation budgets based on the SRTM elevation model. The relationship between phenology and environmental drivers was assessed using general linear models. Temporal displacement in the start date of the annual growth season was more evident than variations in season length among vegetation types, indicating a possible temporal separation in the use of resources. Season length was inversely proportional to elevation, decreasing 1.58 days per 100 m. Green-up and senescence rates were faster where annual temperature amplitude was higher. We found that water and light availability, modulated by topography, are the most likely drivers of land surface phenology in the region, determining the start, end and length of the growing season. Temperature had an important role in determining the rates of leaf development and the strength of vegetation seasonality, suggesting that tropical vegetation is also sensitive to latitudinal temperature changes, regardless of the elevational gradient. Our work improves the current understanding of phenological strategies in the seasonal tropics and emphasizes the importance of topography in shaping light and water availability for leaf development in snow-free mountains.     

North American vegetation patterns observed with the NOAA-7 advanced very high resolution radiometer

Spectral vegetation index measurements derived from remotely sensed observations show great promise as a means to improve knowledge of land vegetation patterns. The daily, global observations acquired by the Advanced Very High Resolution Radiometer, a sensor on the current series of U.S. National Oceanic and Atmospheric Administration meteorological satellites, may be particularly well suited for global studies of vegetation. Preliminary results from analysis of North American observations, extending from April to November 1982, show that the vegetation index patterns observed correspond to the known seasonality of North American natural and cultivated vegetation. Integration of the observations over the growing season produced measurements that are related to net primary productivity patterns of the major North American natural vegetation formations. Regions of intense cultivation were observed as anomalous areas in the integrated growing season measurements. These anomalies can be explained by contrasts between cultivation practices and natural vegetation phenology. Major new information on seasonality, annual extent and interannual variability of vegetation photosynthetic activity at continental and global scales can be derived from these satellite observations.This research is in part supported through NASA Cooperative Agreement NCC 5–26 from the NASA/Goddard Space Flight Center.     

Landscape controls on the timing of spring,autumn, and growing season length in mid‐Atlantic forests

The timing of spring leaf development, trajectories of summer leaf area, and the timing of autumn senescence have profound impacts to the water, carbon, and energy balance of ecosystems, and are likely influenced by global climate change. Limited field‐based and remote‐sensing observations have suggested complex spatial patterns related to geographic features that influence climate. However, much of this variability occurs at spatial scales that inhibit a detailed understanding of even the dominant drivers. Recognizing these limitations, we used nonlinear inverse modeling of medium‐resolution remote sensing data, organized by day of year, to explore the influence of climate‐related landscape factors on the timing of spring and autumn leaf‐area trajectories in mid‐Atlantic, USA forests. We also examined the extent to which declining summer greenness (greendown) degrades the precision and accuracy of observations of autumn offset of greenness. Of the dominant drivers of landscape phenology, elevation was the strongest, explaining up to 70% of the spatial variation in the onset of greenness. Urban land cover was second in importance, influencing spring onset and autumn offset to a distance of 32 km from large cities. Distance to tidal water also influenced phenological timing, but only within ~5 km of shorelines. Additionally, we observed that (i) growing season length unexpectedly increases with increasing elevation at elevations below 275 m; (ii) along gradients in urban land cover, timing of autumn offset has a stronger effect on growing season length than does timing of spring onset; and (iii) summer greendown introduces bias and uncertainty into observations of the autumn offset of greenness. These results demonstrate the power of medium grain analyses of landscape‐scale phenology for understanding environmental controls on growing season length, and predicting how these might be affected by climate change.     

Disentangling the mechanisms behind winter snow impact on vegetation activity in northern ecosystems

Although seasonal snow is recognized as an important component in the global climate system, the ability of snow to affect plant production remains an important unknown for assessing climate change impacts on vegetation dynamics at high‐latitude ecosystems. Here, we compile data on satellite observation of vegetation greenness and spring onset date, satellite‐based soil moisture, passive microwave snow water equivalent (SWE) and climate data to show that winter SWE can significantly influence vegetation greenness during the early growing season (the period between spring onset date and peak photosynthesis timing) over nearly one‐fifth of the land surface in the region north of 30 degrees, but the magnitude and sign of correlation exhibits large spatial heterogeneity. We then apply an assembled path model to disentangle the two main processes (via changing early growing‐season soil moisture, and via changing the growth period) in controlling the impact of winter SWE on vegetation greenness, and suggest that the “moisture” and “growth period” effect, to a larger extent, result in positive and negative snow–productivity associations, respectively. The magnitude and sign of snow–productivity association is then dependent upon the relative dominance of these two processes, with the “moisture” effect and positive association predominating in Central, western North America and Greater Himalaya, and the “growth period” effect and negative association in Central Europe. We also indicate that current state‐of‐the‐art models in general reproduce satellite‐based snow–productivity relationship in the region north of 30 degrees, and do a relatively better job of capturing the “moisture” effect than the “growth period” effect. Our results therefore work towards an improved understanding of winter snow impact on vegetation greenness in northern ecosystems, and provide a mechanistic basis for more realistic terrestrial carbon cycle models that consider the impacts of winter snow processes.     

Satellite Evidence of Phenological Differences Between Urbanized and Rural Areas of the Eastern United States Deciduous Broadleaf Forest

We used a 10-year record (1990–99) of composited and cloud-screened reflectances from the Advanced Very High Resolution Radiometer (AVHRR) to test for phenological differences between urban and rural areas in the eastern United States deciduous broadleaf forest (DBF). We hypothesized that well-documented urban heat island effects would be associated with alterations in temperature-sensitive vegetation phenology. Our objectives were thus (a) to investigate possible differences in the start of the growing season (SOS) and end of the growing season (EOS) between the urban and DBF land covers, (b) to investigate related differences in greenness amplitude and fractional cover, and (c) to develop a generalized additive model (GAM) to predict the spatial variation of observed differences. By analyzing individual 1° latitude by 1° longitude blocks, we found that, on average, urbanization is associated with a growing season expansion of 7.6 days. Most of this effect is caused by an earlier SOS in urban areas. In all cases, urban regions had lower fractional cover and greenness amplitude. The GAM model failed to produce a viable model for differences in EOS, probably because it is dominated by photoperiod controls with only a minor temperature impact. SOS differences were predicted with an accuracy of about 2.4 days, with a GAM consisting of smoothed functions of mean annual average temperature, urban fractional cover, and the urban vs DBF greenness amplitude difference. We speculate that evidence of a phenological response to warming indicates that global warming, without reduction in DBF vegetation cover and greenness amplitude, may increase carbon sequestration in mesic deciduous forests. Received 6 June 2001; accepted 23 October 2001.     

Changes in the structure and function of northern Alaskan ecosystems when considering variable leaf‐out times across groupings of species in a dynamic vegetation model

The phenology of arctic ecosystems is driven primarily by abiotic forces, with temperature acting as the main determinant of growing season onset and leaf budburst in the spring. However, while the plant species in arctic ecosystems require differing amounts of accumulated heat for leaf‐out, dynamic vegetation models simulated over regional to global scales typically assume some average leaf‐out for all of the species within an ecosystem. Here, we make use of air temperature records and observations of spring leaf phenology collected across dominant groupings of species (dwarf birch shrubs, willow shrubs, other deciduous shrubs, grasses, sedges, and forbs) in arctic and boreal ecosystems in Alaska. We then parameterize a dynamic vegetation model based on these data for four types of tundra ecosystems (heath tundra, shrub tundra, wet sedge tundra, and tussock tundra), as well as ecotonal boreal white spruce forest, and perform model simulations for the years 1970–2100. Over the course of the model simulations, we found changes in ecosystem composition under this new phenology algorithm compared with simulations with the previous phenology algorithm. These changes were the result of the differential timing of leaf‐out, as well as the ability for the groupings of species to compete for nitrogen and light availability. Regionally, there were differences in the trends of the carbon pools and fluxes between the new phenology algorithm and the previous phenology algorithm, although these differences depended on the future climate scenario. These findings indicate the importance of leaf phenology data collection by species and across the various ecosystem types within the highly heterogeneous Arctic landscape, and that dynamic vegetation models should consider variation in leaf‐out by groupings of species within these ecosystems to make more accurate projections of future plant distributions and carbon cycling in Arctic regions.     

Seasonal carbon dynamics and water fluxes in an Amazon rainforest

Satellite‐based observations indicate that seasonal patterns in canopy greenness and productivity in the Amazon are negatively correlated with precipitation, with increased greenness occurring during the dry months. Flux tower measurements indicate that the canopy greening that occurs during the dry season is associated with increases in net ecosystem productivity (NEP) and evapotranspiration (ET). Land surface and terrestrial biosphere model simulations for the region have predicted the opposite of these observed patterns, with significant declines in greenness, NEP, and ET during the dry season. In this study, we address this issue mainly by developing an empirically constrained, light‐controlled phenology submodel within the Ecosystem Demography model version 2 (ED2). The constrained ED2 model with a suite of field observations shows markedly improved predictions of seasonal ecosystem dynamics, more accurately capturing the observed patterns of seasonality in water, carbon, and litter fluxes seen at the Tapajos National Forest, Brazil (2.86°S, 54.96°W). Long‐term simulations indicate that this light‐controlled phenology increases the resilience of Amazon forest NEP to interannual variability in climate forcing.     

Successional and seasonal variations in soil and litter microbial community structure and function during tropical postagricultural forest regeneration: a multiyear study

Soil microorganisms regulate fundamental biochemical processes in plant litter decomposition and soil organic matter (SOM) transformations. Understanding how microbial communities respond to changes in vegetation is critical for improving predictions of how land‐cover change affects belowground carbon storage and nutrient availability. We measured intra‐ and interannual variability in soil and forest litter microbial community composition and activity via phospholipid fatty acid analysis (PLFA) and extracellular enzyme activity across a well‐replicated, long‐term chronosequence of secondary forests growing on abandoned pastures in the wet subtropical forest life zone of Puerto Rico. Microbial community PLFA structure differed between young secondary forests and older secondary and primary forests, following successional shifts in tree species composition. These successional patterns held across seasons, but the microbial groups driving these patterns differed over time. Microbial community composition from the forest litter differed greatly from those in the soil, but did not show the same successional trends. Extracellular enzyme activity did not differ with forest succession, but varied by season with greater rates of potential activity in the dry seasons. We found few robust significant relationships among microbial community parameters and soil pH, moisture, carbon, and nitrogen concentrations. Observed inter‐ and intrannual variability in microbial community structure and activity reveal the importance of a multiple, temporal sampling strategy when investigating microbial community dynamics with land‐use change. Successional control over microbial composition with forest recovery suggests strong links between above and belowground communities.     

Global vegetation productivity response to climatic oscillations during the satellite era

Climate control on global vegetation productivity patterns has intensified in response to recent global warming. Yet, the contributions of the leading internal climatic variations to global vegetation productivity are poorly understood. Here, we use 30 years of global satellite observations to study climatic variations controls on continental and global vegetation productivity patterns. El Niño‐Southern Oscillation (ENSO) phases (La Niña, neutral, and El Niño years) appear to be a weaker control on global‐scale vegetation productivity than previously thought, although continental‐scale responses are substantial. There is also clear evidence that other non‐ENSO climatic variations have a strong control on spatial patterns of vegetation productivity mainly through their influence on temperature. Among the eight leading internal climatic variations, the East Atlantic/West Russia Pattern extensively controls the ensuing year vegetation productivity of the most productive tropical and temperate forest ecosystems of the Earth's vegetated surface through directionally consistent influence on vegetation greenness. The Community Climate System Model (CCSM4) simulations do not capture the observed patterns of vegetation productivity responses to internal climatic variations. Our analyses show the ubiquitous control of climatic variations on vegetation productivity and can further guide CCSM and other Earth system models developments to represent vegetation response patterns to unforced variability. Several winter time internal climatic variation indices show strong potentials on predicting growing season vegetation productivity two to six seasons ahead which enables national governments and farmers forecast crop yield to ensure supplies of affordable food, famine early warning, and plan management options to minimize yield losses ahead of time.     

Measuring phenological variability from satellite imagery

Abstract. Vegetation phenological phenomena are closely related to seasonal dynamics of the lower atmosphere and are therefore important elements in global models and vegetation monitoring. Normalized difference vegetation index (NDVI) data derived from the National Oceanic and Atmospheric Administration's Advanced Very High Resolution Radiometer (AVHRR) satellite sensor offer a means of efficiently and objectively evaluating phenological characteristics over large areas. Twelve metrics linked to key phenological events were computed based on time-series NDVI data collected from 1989 to 1992 over the conterminous United States. These measures include the onset of greenness, time of peak NDVI, maximum NDVI, rate of greenup, rate of senescence, and integrated NDVI. Measures of central tendency and variability of the measures were computed and analyzed for various land cover types. Results from the analysis showed strong coincidence between the satellite-derived metrics and predicted phenological characteristics. In particular, the metrics identified interannual variability of spring wheat in North Dakota, characterized the phenology of four types of grasslands, and established the phenological consistency of deciduous and coniferous forests. These results have implications for large-area land cover mapping and monitoring. The utility of remotely sensed data as input to vegetation mapping is demonstrated by showing the distinct phenology of several land cover types. More stable information contained in ancillary data should be incorporated into the mapping process, particularly in areas with high phenological variability. In a regional or global monitoring system, an increase in variability in a region may serve as a signal to perform more detailed land cover analysis with higher resolution imagery.     

Satellite microwave detection of boreal forest recovery from the extreme 2004 wildfires in Alaska and Canada

The rate of vegetation recovery from boreal wildfire influences terrestrial carbon cycle processes and climate feedbacks by affecting the surface energy budget and land‐atmosphere carbon exchange. Previous forest recovery assessments using satellite optical‐infrared normalized difference vegetation index (NDVI) and tower CO₂ eddy covariance techniques indicate rapid vegetation recovery within 5–10 years, but these techniques are not directly sensitive to changes in vegetation biomass. Alternatively, the vegetation optical depth (VOD) parameter from satellite passive microwave remote sensing can detect changes in canopy biomass structure and may provide a useful metric of post‐fire vegetation response to inform regional recovery assessments. We analyzed a multi‐year (2003–2010) satellite VOD record from the NASA AMSR‐E (Advanced Microwave Scanning Radiometer for EOS) sensor to estimate forest recovery trajectories for 14 large boreal fires from 2004 in Alaska and Canada. The VOD record indicated initial post‐fire canopy biomass recovery within 3–7 years, lagging NDVI recovery by 1–5 years. The VOD lag was attributed to slower non‐photosynthetic (woody) and photosynthetic (foliar) canopy biomass recovery, relative to the faster canopy greenness response indicated from the NDVI. The duration of VOD recovery to pre‐burn conditions was also directly proportional (P < 0.01) to satellite (moderate resolution imaging spectroradiometer) estimated tree cover loss used as a metric of fire severity. Our results indicate that vegetation biomass recovery from boreal fire disturbance is generally slower than reported from previous assessments based solely on satellite optical‐infrared remote sensing, while the VOD parameter enables more comprehensive assessments of boreal forest recovery.     

The greening and browning of Alaska based on 1982–2003 satellite data

Aim To examine the trends of 1982–2003 satellite‐derived normalized difference vegetation index (NDVI) values at several spatial scales within tundra and boreal forest areas of Alaska. Location Arctic and subarctic Alaska. Methods Annual maximum NDVI data from the twice monthly Global Inventory Modelling and Mapping Studies (GIMMS) NDVI 1982–2003 data set with 64‐km² pixels were extracted from a spatial hierarchy including three large regions: ecoregion polygons within regions, ecozone polygons within boreal ecoregions and 100‐km climate station buffers. The 1982–2003 trends of mean annual maximum NDVI values within each area, and within individual pixels, were computed using simple linear regression. The relationship between NDVI and temperature and precipitation was investigated within climate station buffers. Results At the largest spatial scale of polar, boreal and maritime regions, the strongest trend was a negative trend in NDVI within the boreal region. At a finer scale of ecoregion polygons, there was a strong positive NDVI trend in cold arctic tundra areas, and a strong negative trend in interior boreal forest areas. Within boreal ecozone polygons, the weakest negative trends were from areas with a maritime climate or colder mountainous ecozones, while the strongest negative trends were from warmer basin ecozones. The trends from climate station buffers were similar to ecoregion trends, with no significant trends from Bering tundra buffers, significant increasing trends among arctic tundra buffers and significant decreasing trends among interior boreal forest buffers. The interannual variability of NDVI among the arctic tundra buffers was related to the previous summer warmth index. The spatial pattern of increasing tundra NDVI at the pixel level was related to the west‐to‐east spatial pattern in changing climate across arctic Alaska. There was no significant relationship between interannual NDVI and precipitation or temperature among the boreal forest buffers. The decreasing NDVI trend in interior boreal forests may be due to several factors including increased insect/disease infestations, reduced photosynthesis and a change in root/leaf carbon allocation in response to warmer and drier growing season climate. Main conclusions There was a contrast in trends of 1982–2003 annual maximum NDVI, with cold arctic tundra significantly increasing in NDVI and relatively warm and dry interior boreal forest areas consistently decreasing in NDVI. The annual maximum NDVI from arctic tundra areas was strongly related to a summer warmth index, while there were no significant relationships in boreal areas between annual maximum NDVI and precipitation or temperature. Annual maximum NDVI was not related to spring NDVI in either arctic tundra or boreal buffers.     

Interannual variations and trends in global land surface phenology derived from enhanced vegetation index during 1982–2010

Land surface phenology is widely retrieved from satellite observations at regional and global scales, and its long-term record has been demonstrated to be a valuable tool for reconstructing past climate variations, monitoring the dynamics of terrestrial ecosystems in response to climate impacts, and predicting biological responses to future climate scenarios. This study detected global land surface phenology from the advanced very high resolution radiometer (AVHRR) and the Moderate Resolution Imaging Spectroradiometer (MODIS) data from 1982 to 2010. Based on daily enhanced vegetation index at a spatial resolution of 0.05 degrees, we simulated the seasonal vegetative trajectory for each individual pixel using piecewise logistic models, which was then used to detect the onset of greenness increase (OGI) and the length of vegetation growing season (GSL). Further, both overall interannual variations and pixel-based trends were examined across Koeppen’s climate regions for the periods of 1982–1999 and 2000–2010, respectively. The results show that OGI and GSL varied considerably during 1982–2010 across the globe. Generally, the interannual variation could be more than a month in precipitation-controlled tropical and dry climates while it was mainly less than 15 days in temperature-controlled temperate, cold, and polar climates. OGI, overall, shifted early, and GSL was prolonged from 1982 to 2010 in most climate regions in North America and Asia while the consistently significant trends only occurred in cold climate and polar climate in North America. The overall trends in Europe were generally insignificant. Over South America, late OGI was consistent (particularly from 1982 to 1999) while either positive or negative GSL trends in a climate region were mostly reversed between the periods of 1982–1999 and 2000–2010. In the Northern Hemisphere of Africa, OGI trends were mostly insignificant, but prolonged GSL was evident over individual climate regions during the last 3 decades. OGI mainly showed late trends in the Southern Hemisphere of Africa while GSL was reversed from reduced GSL trends (1982–1999) to prolonged trends (2000–2010). In Australia, GSL exhibited considerable interannual variation, but the consistent trend lacked presence in most regions. Finally, the proportion of pixels with significant trends was less than 1 % in most of climate regions although it could be as large as 10 %.     

日光诱导叶绿素荧光对亚热带常绿针叶林物候的追踪

植被物候期(春季返青和秋季衰老)是表征生物响应和陆地碳循环的基础信息。由于常绿针叶林冠层绿度的季节变动较弱,遥感提取常绿针叶林的物候信息存在着较大的不确定性,是目前区域物候监测中的难点。利用MODIS植被指数(归一化植被指数NDVI和增强型植被指数EVI)、GOME-2日光诱导叶绿素荧光(SIF)和通量数据(总初级生产力GPP)估算2007—2011年亚热带常绿针叶林物候期,用来比较三类遥感指数估算常绿针叶林物候的差异。结果表明:基于表征光合作用物候的通量GPP数据估算得到5年内亚热带常绿针叶林生长季开始时间(SOS_(GPP))为第63天,生长季结束时间(EOS_(GPP))为第324天,生长季长度为272天;基于反映植被光合作用特征的SIF曲线获得物候信息要滞后GPP物候期,其中生长季开始时间滞后19天,生长季结束时间滞后2天;基于传统植被指数NDVI和EVI的物候期滞后GPP物候期的时间要大于SIF滞后期,其中植被指数SOS滞后SOS_(GPP)31天,植被指数EOS滞后EOS_(GPP)10—17天。虽然基于3种遥感指数估算的春季和秋季物候都滞后于通量GPP的物候期,但是卫星SIF的物候信息能够更好地捕捉常绿针叶林的生长阶段。同时,春季温度是影响森林生长季开始时间的最重要因素;秋季水分和辐射是影响生长季结束时间的关键因素。由此可见,SIF估算的亚热带常绿针叶林的春季和秋季物候的滞后时间要短于传统植被指数,能更好地追踪常绿林光合作用的季节性,为深入研究陆地生态系统碳循环及其对气候变化的响应提供重要的基础。     

Strong contribution of autumn phenology to changes in satellite‐derived growing season length estimates across Europe (1982–2011)

Land Surface Phenology (LSP) is the most direct representation of intra‐annual dynamics of vegetated land surfaces as observed from satellite imagery. LSP plays a key role in characterizing land‐surface fluxes, and is central to accurately parameterizing terrestrial biosphere–atmosphere interactions, as well as climate models. In this article, we present an evaluation of Pan‐European LSP and its changes over the past 30 years, using the longest continuous record of Normalized Difference Vegetation Index (NDVI) available to date in combination with a landscape‐based aggregation scheme. We used indicators of Start‐Of‐Season, End‐Of‐Season and Growing Season Length (SOS, EOS and GSL, respectively) for the period 1982–2011 to test for temporal trends in activity of terrestrial vegetation and their spatial distribution. We aggregated pixels into ecologically representative spatial units using the European Landscape Classification (LANMAP) and assessed the relative contribution of spring and autumn phenology. GSL increased significantly by 18–24 days decade^?1 over 18–30% of the land area of Europe, depending on methodology. This trend varied extensively within and between climatic zones and landscape classes. The areas of greatest growing‐season lengthening were the Continental and Boreal zones, with hotspots concentrated in southern Fennoscandia, Western Russia and pockets of continental Europe. For the Atlantic and Steppic zones, we found an average shortening of the growing season with hotspots in Western France, the Po valley, and around the Caspian Sea. In many zones, changes in the NDVI‐derived end‐of‐season contributed more to the GSL trend than changes in spring green‐up, resulting in asymmetric trends. This underlines the importance of investigating senescence and its underlying processes more closely as a driver of LSP and global change.     

Comparison of Phenology Models for Predicting the Onset of Growing Season over the Northern Hemisphere

Vegetation phenology models are important for examining the impact of climate change on the length of the growing season and carbon cycles in terrestrial ecosystems. However, large uncertainties in present phenology models make accurate assessment of the beginning of the growing season (BGS) a challenge. In this study, based on the satellite-based phenology product (i.e. the V005 MODIS Land Cover Dynamics (MCD12Q2) product), we calibrated four phenology models, compared their relative strength to predict vegetation phenology; and assessed the spatial pattern and interannual variability of BGS in the Northern Hemisphere. The results indicated that parameter calibration significantly influences the models'' accuracy. All models showed good performance in cool regions but poor performance in warm regions. On average, they explained about 67% (the Growing Degree Day model), 79% (the Biome-BGC phenology model), 73% (the Number of Growing Days model) and 68% (the Number of Chilling Days-Growing Degree Day model) of the BGS variations over the Northern Hemisphere. There were substantial differences in BGS simulations among the four phenology models. Overall, the Biome-BGC phenology model performed best in predicting the BGS, and showed low biases in most boreal and cool regions. Compared with the other three models, the two-phase phenology model (NCD-GDD) showed the lowest correlation and largest biases with the MODIS phenology product, although it could catch the interannual variations well for some vegetation types. Our study highlights the need for further improvements by integrating the effects of water availability, especially for plants growing in low latitudes, and the physiological adaptation of plants into phenology models.