THE RESPONSE OF UK VEGETATION TO ELEVATED TEMPERATURES IN 2006: COUPLING ENVISAT MERIS TERRESTRIAL CHLOROPHYLL INDEX (MTCI) AND MEAN AIR TEMPERATURE

Similar documents
Permanent City Research Online URL:

THE RELATIONSHIP BETWEEN THE MERIS TERRESTRIAL CHLOROPHYLL INDEX AND CHLOROPHYLL CONTENT

SATellite MOnitoring of VEgetation Seasonality SAT-MOVES. Brian O Connor, Ned Dwyer, Fiona Cawkwell

A Study of the Uncertainty in Future Caribbean Climate Using the PRECIS Regional Climate Model

DROUGHT ASSESSMENT USING SATELLITE DERIVED METEOROLOGICAL PARAMETERS AND NDVI IN POTOHAR REGION

SHORT COMMUNICATION EXPLORING THE RELATIONSHIP BETWEEN THE NORTH ATLANTIC OSCILLATION AND RAINFALL PATTERNS IN BARBADOS

Spatial Drought Assessment Using Remote Sensing and GIS techniques in Northwest region of Liaoning, China

Seasonal and interannual relations between precipitation, soil moisture and vegetation in the North American monsoon region

CLIMATE CHANGE IMPACTS ON HYDROMETEOROLOGICAL VARIABLES AT LAKE KARLA WATERSHED

Satellite-derived environmental drivers for top predator hotspots

How does the physical environment influence communities and ecosystems? Hoodoos in Cappadocia, Turkey

8.2 GLOBALLY DESCRIBING THE CURRENT DAY LAND SURFACE AND HISTORICAL LAND COVER CHANGE IN CCSM 3.0 USING AVHRR AND MODIS DATA AT FINE SCALES

Shift in Onset of First Inter Monsoon Rain in Coconut Growing Areas in Sri Lanka

A Time Series of Photo-synthetically Available Radiation at the Ocean Surface from SeaWiFS and MODIS Data

Name ECOLOGY TEST #1 Fall, 2014

Using MERIS and MODIS for Land Cover Mapping in the Netherlands

Assessment of Snow Cover Vulnerability over the Qinghai-Tibetan Plateau

Using digital cameras to monitoring vegetation phenology: Insights from PhenoCam. Andrew D. Richardson Harvard University

Full Version with References: Future Climate of the European Alps

Analysis of meteorological measurements made over three rainy seasons in Sinazongwe District, Zambia.

ASSESSMENT OF DIFFERENT WATER STRESS INDICATORS BASED ON EUMETSAT LSA SAF PRODUCTS FOR DROUGHT MONITORING IN EUROPE

Assimilation of satellite fapar data within the ORCHIDEE biosphere model and its impacts on land surface carbon and energy fluxes

ANALYSIS AND VALIDATION OF A METHODOLOGY TO EVALUATE LAND COVER CHANGE IN THE MEDITERRANEAN BASIN USING MULTITEMPORAL MODIS DATA

RESOLUTION ERRORS ASSOCIATED WITH GRIDDED PRECIPITATION FIELDS

Proceedings of the Second Working Meeting on MERIS and AATSR Calibration and Geophysical Validation (MAVT-2006)

Ryan P. Shadbolt * Central Michigan University, Mt. Pleasant, Michigan

Short Communication Shifting of frozen ground boundary in response to temperature variations at northern China and Mongolia,

Drought Bulletin for the Greater Horn of Africa: Situation in June 2011

J8.4 TRENDS OF U.S. SNOWFALL AND SNOW COVER IN A WARMING WORLD,

Assimilating terrestrial remote sensing data into carbon models: Some issues

SEASONAL VARIABILITY AND PERSISTENCE IN TEMPERATURE SCENARIOS FOR ICELAND

World Geography Chapter 3

Chapter outline. Reference 12/13/2016

Drought and its effect on vegetation, comparison of NDVI for drought and non-drought years related to Land use classifications

Vegetation Remote Sensing

8.2 Tropospheric ozone

Leveraging Sentinel-1 time-series data for mapping agricultural land cover and land use in the tropics

Global Biogeography. Natural Vegetation. Structure and Life-Forms of Plants. Terrestrial Ecosystems-The Biomes

Drought Estimation Maps by Means of Multidate Landsat Fused Images

MERIS SMILE EFFECT CHARACTERISATION AND CORRECTION DOCUMENT. document title/ titre du document. prepared by/préparé par MERIS ESL

Impact of the 2002 stratospheric warming in the southern hemisphere on the tropical cirrus clouds and convective activity

ESTIMATING SNOWMELT CONTRIBUTION FROM THE GANGOTRI GLACIER CATCHMENT INTO THE BHAGIRATHI RIVER, INDIA ABSTRACT INTRODUCTION

3. Carbon Dioxide (CO 2 )

2007: The Netherlands in a drought again (2 May 2007)

Module 11: Meteorology Topic 3 Content: Climate Zones Notes

Seasonal Climate Forecast August October 2013 Verification (Issued: November 17, 2013)

Analysis of Relative Humidity in Iraq for the Period

Remote detection of giant reed invasions in riparian habitats: challenges and opportunities for management planning

Changes Of Glacier Extent: Using Gis To Integration Maps, Remote Sensing Material And Field Surveying. Petuaniabukta, Central Spitsbergen Case Study

Impact of interannual variability of meteorological parameters on NDVI over Mongolian

N Management in Potato Production. David Mulla, Carl Rosen, Tyler Nigonand Brian Bohman Dept. Soil, Water & Climate University of Minnesota

International Journal of Scientific and Research Publications, Volume 3, Issue 5, May ISSN

A high spectral resolution global land surface infrared emissivity database

Greening of Arctic: Knowledge and Uncertainties

FUTURE PROJECTIONS OF PRECIPITATION CHARACTERISTICS IN ASIA

Statistical Analysis of Long Term Temporal Trends of Precipitation and Temperature in Wainganga Sub-Basin, India

Defining microclimates on Long Island using interannual surface temperature records from satellite imagery

Attribution of Estonian phyto-, ornitho- and ichtyophenological trends with parameters of changing climate

Lecture Topics. 1. Vegetation Indices 2. Global NDVI data sets 3. Analysis of temporal NDVI trends

Spatial Downscaling of TRMM Precipitation Using DEM. and NDVI in the Yarlung Zangbo River Basin

Precipitation patterns alter growth of temperate vegetation

Mozambique. General Climate. UNDP Climate Change Country Profiles. C. McSweeney 1, M. New 1,2 and G. Lizcano 1

Scaling photosynthetic light-use efficiency from canopies to. landscapes

SHORT COMMUNICATION HOW WELL DO THE CENTRAL ENGLAND TEMPERATURE AND THE ENGLAND AND WALES PRECIPITATION SERIES REPRESENT THE CLIMATE OF THE UK?

Satellite-based Lake Surface Temperature (LST) Homa Kheyrollah Pour Claude Duguay

ISPRS Archives XXXVIII-8/W3 Workshop Proceedings: Impact of Climate Change on Agriculture

Climate briefing. Wellington region, May Alex Pezza and Mike Thompson Environmental Science Department

EXTRACTION OF REMOTE SENSING INFORMATION OF BANANA UNDER SUPPORT OF 3S TECHNOLOGY IN GUANGXI PROVINCE

Season Spotter: Using Citizen Science to Validate and Scale Plant Phenology from Near-Surface Remote Sensing

AN ACCURACY ASSESSMENT OF AATSR LST DATA USING EMPIRICAL AND THEORETICAL METHODS

Land Cover Project ESA Climate Change Initiative. Processing chain for land cover maps dedicated to climate modellers.

UK Contribution to the European CORINE Land Cover

CALIBRATION OF THERMAL INFRARED CHANNELS OF FENG YUN SATELLITE DATA. Wan Hazli Wan Kadir, Abd Wahid Rasib, Ab. Latif Ibrahim and Zawatil Azhani Zaini

Biology and the hierarchies of life. finer scale. coarser scale. individual. populations. metapopulation. community. ecosystem. biome.

Physiological (Ecology of North American Plant Communities

The Relationship between Vegetation Changes and Cut-offs in the Lower Yellow River Based on Satellite and Ground Data

Data Fusion and Multi-Resolution Data

RESEARCH METHODOLOGY

Analysis of Historical Pattern of Rainfall in the Western Region of Bangladesh

Assessing Drought in Agricultural Area of central U.S. with the MODIS sensor

Unseasonable weather conditions in Japan in August 2014

Landuse and Landcover change analysis in Selaiyur village, Tambaram taluk, Chennai

Jay Lawrimore NOAA National Climatic Data Center 9 October 2013

Energy Systems, Structures and Processes Essential Standard: Analyze patterns of global climate change over time Learning Objective: Differentiate

Table 13 Maximum dewpoint temperature (degrees Celsius), primary meteorological station, H. J. Andrews Experimental Forest, 5/10/72 through 12/31/84

remain on the trees all year long) Example: Beaverlodge, Alberta, Canada

Introduction. Populus trichocarpa TORR. and GRAY. By M. G. R. CANNELL and S. C. WILLETT

The Wide Dynamic Range Vegetation Index and its Potential Utility for Gap Analysis

What is a Biome? An Overview of Biomes. The Holdridge Life Zones. Tundra 9/14/2010. In the following slides, you ll

2015: A YEAR IN REVIEW F.S. ANSLOW

Climate variability rather than overstocking causes recent large scale cover changes of Tibetan pastures

Climate Change and Vegetation Phenology

Climate Modeling Research & Applications in Wales. John Houghton. C 3 W conference, Aberystwyth

Earth s Major Terrerstrial Biomes. *Wetlands (found all over Earth)

SEASONAL AND ANNUAL TRENDS OF AUSTRALIAN MINIMUM/MAXIMUM DAILY TEMPERATURES DURING

WHAT CAN MAPS TELL US ABOUT THE GEOGRAPHY OF ANCIENT GREECE? MAP TYPE 1: CLIMATE MAPS

Regionalized and application-specific compositing - a review of requirements, opportunities and challenges

CHAPTER 1: INTRODUCTION

VIIRS narrowband to broadband land surface albedo conversion: formula and validation

Monitoring Growing Season Length of Deciduous Broad Leaf Forest Derived From Satellite Data in Iran

Transcription:

THE RESPONSE OF UK VEGETATION TO ELEVATED TEMPERATURES IN 2006: COUPLING ENVISAT MERIS TERRESTRIAL CHLOROPHYLL INDEX (MTCI) AND MEAN AIR TEMPERATURE Samuel Almond a, Doreen S. Boyd b, Paul J. Curran c and Jadunandan Dash d a School of Conservation Sciences, Bournemouth University, Talbot Campus, Bournemouth BH12 5BB, UK salmond@bournemouth.ac.uk b School of Geography, University of Nottingham, University Park, Nottingham NG7 2RD, UK c Office of the Vice-Chancellor, Bournemouth University, Talbot Campus, Bournemouth BH12 5BB, UK d School of Geography, University of Southampton, Southampton SO17 1BJ, UK KEY WORDS: Envisat MTCI, index, vegetation growing season, phenology, senescence, chlorophyll, air temperature ABSTRACT: Climate is one of the key variables driving ecosystems at local to global scales. How and to what extent vegetation responds to climate variability is a challenging topic for global change analysis. Autumn 2006 was the warmest in the UK since records began in 1659 (Meteorological Office, 16th October 2006) and this delayed the seasonal senescence of vegetation. Earth observation data provides an opportunity to study this phenomenon and put it into the context of ecosystem dynamics over a longer time period. The European Space Agency (ESA) uses data recorded by the Medium Resolution Imaging Spectrometer (MERIS) in the red / NIR bands to produce an operational product called the MERIS Terrestrial Chlorophyll Index (MTCI). The MTCI is related to the position of the red edge in vegetation spectra and can be used to determine the chlorophyll content of vegetation and thus its phenology. The 2006 growing season was investigated using MTCI data. Weekly MTCI composites for 2003-2006 were used to produce land-cover specific MTCI time series which were related to temporally-coincident temperature data. Analysis demonstrated that for forest and heath and grassland (i) the MTCI was sensitive to seasonal variation in inferred chlorophyll content, (ii) the high autumnal temperatures were related to an extended growing season and (ii) there was a strong relationship between seasonal temperature variation and MTCI. 1.1 Background 1. INTRODUCTION The most recent assessment report of the Intergovernmental Panel on Climate Change (IPCC 2001) predicted further increases in global mean temperature, climate variability and extreme events. Phenology is one of the main bio-indicators of climate change impacts on ecosystems (Schleip et al., 2006). Seasonal temperature change in temperate forests is known to start both tree growth and leafing (Sparks et al., 2005, Fisher et al., 2007) and senescence. The forecasted change in climate is likely to have consequences for land use practices (Thomas, 2006) and vegetation phenology is an important indicator of ecosystem response to climate change. It has been established that air temperature is the most important controlling factor related to trees phenology in the temperate latitudes, whereas precipitation and photoperiod play a less pronounced role in phenological development (Chen et al., 2005). Vegetation response to variation in seasonal temperature is likely to affect photosynthetic rates (Bassow and Bazzaz, 1998) however; autumn senescence in deciduous vegetation is triggered both by environmental conditions and leaf specific maturation processes. The effects of temperature on leaf senescence is difficult to identify because factors, such as air temperature, day length and rainfall, often change at the same time (Rosenthal and Camm, 1997). The ability to couple observations of vegetation phenology with climatic factors over large scales is vital to monitor, predict and manage impacts of climatic change on ecosystems. Remote sensing data have been used in the study of vegetation condition and seasonal vegetation dynamics for many years (Reed and Bradley, 2006). The ability to observe vegetation phenology remotely over large areas provides a unique opportunity to monitor the effects of climatic change on vegetated canopies at local to global scales. The vast majority of studies have used data from the Advanced Very High Resolution Radiometer (Zhang et al., 2006). A new generation of remote sensing data sources are now available and they improve greatly our ability to identify changes in ecosystem phenology. 1.2 Envisat MTCI One new remote sensing data source is the Envisat MERIS Terrestrial Chlorophyll Index (MTCI). The MTCI was designed to monitor vegetation condition via an estimation of chlorophyll content. Vegetation growth cycles can be characterised through changes in chlorophyll concentration and leaf area index, which determine chlorophyll content (Curran et al. 2007). Methods used to determine chlorophyll content from remote sensing data have focused upon locating the red edge position (REP) between the red absorption feature and the near infrared reflectance maxima of a vegetation reflectance spectrum (Curran et al. 2007). However, such methods are time-

consuming and are not accurate indicators of chlorophyll content at high chlorophyll contents. The MTCI overcomes the potential limitations of using the REP to derive chlorophyll content. It is calculated using the ratio of the difference in reflectance between band 10 and band 9 and the difference in reflectance between band 9 and band 8 of the MERIS standard band setting: Site New Forest National Park, Hampshire Vegetation type Ancient deciduous forest heath and grassland area of Bodmin Moor, Cornwall Wildwood, Kent heathland MTCI = R753.75 - R708.75 / R708.75 - R681.25 where R753.75, R708.75, R681.25 are reflectance in the centre wavelengths (nm) of the MERIS standard band setting. The MTCI is very simple to calculate yet it is sensitive to all and notably, high values of chlorophyll content. This coupled with the virtues of the MERIS sensor (e.g., radiometrically it is the most accurate imaging spectrometer in space (Curran and Steele, 2004)); fine spectral resolution, moderate spatial resolution (300m) and three-day repeat cycle) (Delwart et al., 2007) has lead to the adoption of the MTCI as an ESA Level 2 land product. Given that the MTCI is the only available chlorophyll index from a spaceborne sensor there is now real opportunity for monitoring vegetation function and condition systemically and reliably and thus the MTCI has now been used in several applications with success (Curran and Dash, 2006; Dash and Curran, 2007). Initial evaluation of the MTCI demonstrates a strong positive relationship with chlorophyll content as measured in the field (r2 of 0.8) (Dash and Curran 2007). The Meteorological Office and the Royal Netherlands Meteorological Institute recorded May to October 2006 as the warmest on record in the UK since 1659. These observed temperatures led to reports of an extended growing season and thus delayed senescence (BBC News accessed 30th October 2006). This paper aims to explore the response of natural vegetation at selected sites in England to this anomaly using a 2003 to 2006 time-series of the MTCI. 2. STUDY SITES AND AUXILIARY DATA Bringwood, Hereforeshire Area 571km2, of which 223km2 is and ancient, 164.5km2 is designated heath and grassland 14.0km2 14.2km2 14.5km2 Table 1. Locations of natural vegetation used in this study 2.1 Satellite sensor data High spatial resolution Terra Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and aerial imagery (25cm true colour composites) were used to confirm both the accuracy of the GIS and that no change in land cover had taken place during 2003 2006. Envisat MERIS reduced resolution (1.2km) data was used to produce 8 day MTCI composites. Radiometric and geometric corrections using BEAM software (the UK-MM-PAF) preceded the derivation of a level 2 product. An arithmetic mean composite was produced from the imagery for each 8 day period during the growing season (March to November) of the vegetation for 2003 2006. The composites were then stacked chronologically for each year. Areas of interest were derived from the GIS layers and used to locate pixels located inside the chosen sites. The pixels were then averaged to produce a single value for each 8 day composite in the layer stack and a temporal profile produced for each site. The selected study sites were located in Hampshire, Cornwall, Herefordshire and Kent, England (Figure 1). Sites were relatively flat with near homogeneous vegetation cover and were larger than 3.6km2 (approximately 3x3 MERIS pixel in reduced resolution) and at geographical locations throughout England. GIS vector layers produced by English Nature were used to determine cover type and site size (Table 1). Figure 1. Location of the study sites within England Figure 2. A summary of the methodology used to monitor vegetation phenology

2.2 Metrological data Climate data within 20 km of each study site (local station) together with the Central England Temperature (CET) series was provided by the Metrological Office (UK). Average daily temperatures (T mean ) were calculated as the mean of the daily maximum (T max ) and the daily minimum (T min ) for both the local station and CET series and used in this study. The CET is often used in broad UK studies to typify the English climate as it is strongly correlated (p < 0.001) with local station observations throughout the UK (Croxton et al., 2006). 3. RESULTS AND DISCUSSION Temporal MTCI profile for all sites indicated an increase MTCI (and thus inferred chlorophyll content) for the autumnal period of 2006 compared to the same period for each of the earlier three years. Late October 2006 MTCI values were comparable with mid September values for the years 2003 and 2004 indicating a delay in senescence of the vegetation canopies. As an example Figure 3 illustrates the temporal MTCI profiles for the (a) lowland heath and grassland and (b) lowland in the New Forest National Park. Similar trends were observed in all three forest sites; there was a very strong positive correlation between MTCI for all forest sites in each year (e.g., r 2 = 0.97 in 2006). Indicating that changes in MTCI is a response to regional and national rather than local factors. Figure 3. Seasonal MTCI derived from layer stacks for the New Forest study site (a) lowland heath and grassland (b) Both the local T mean and CET T mean data were correlated with MTCI for each study site. A strong positive correlation between local T mean and MTCI was observed at all sites. Figure 4a shows this relationship for forest and heath and grassland for the New Forest site for the year of 2006 and Figure 4b shows this relationship for the period 2003 2006 for forest at the New Forest site. Correlation analysis demonstrated a strong positive relationship between T mean and MTCI for all study sites (mean r 2 for period 2003 2006; Bodmin moor 0.92, Bringwood 0.90 and Wildwood 0.88), indicating that the inferred chlorophyll canopy content is correlated strongly to mean monthly temperature, carbon dioxide, when water/nutrients and light are not limiting growth of vegetation is controlled by temperature (Deng et al. 2007). T mean is potentially a limiting factor to vegetation growth. Recent climate change has been found to affect species phenology, especially in the earlier onset of spring events in mid- and higher latitudes, whilst there is not enough information on autumn events to draw firm conclusions (Menzel et al. 2006). Findings in this paper support the notion of delayed senescence due to favourable growing conditions. Figure 4. (a) Relationship between MTCI and local T mean for woodland and grass / heath land sites in the New Forest using 2006 data; (b) the correlation between MTCI and local T mean for woodland sites in the New Forest for years 2003-2006, demonstrating the consistent relationship between MTCI and local T mean. Investigation of the difference between both MTCI and CET mean for 2006 and MTCI and CET mean for the previous three years

revealed a marked increase in temperature and inferred chlorophyll content for all sites. From the period July November 2006 mean inferred chlorophyll content was higher, delaying seasonal senescence. Statistical analysis needs to determine whether inferred chlorophyll content was significantly higher throughout the autumn period. However, during mid-october the mean inferred chlorophyll content for all sites was 54% greater than in previous years. This trend (Figure 5) follows that of T mean, however a response lag of approximately two weeks between temperature and MTCI was observed. Figure 5 also suggests that a seasonal shift may have occurred. Rather than experiencing only an extended growing season in 2006, the vegetation canopies may also have experienced a delay in the greening-up period. This is suggested by the below average inferred chlorophyll content until the month of July. Figures 3a&b also indicates the same trend of a slower greening up period in 2006. A longer running average of yearly observations would be needed to confirm this observation. 5. ACKNOWLEDGEMENTS We thank all those who helped with the data acquisition and processing; in particular Infoterra for providing the majority of the MTCI weekly composites used, the Meteorological Office (UK) for climate data and Bournemouth University for a studentship to SA. 6. REFERENCES Bassow, S.L. and Bazzaz, F.A., 1998. How environmental conditions affect canopy leaf-level photosynthesis in four deciduous tree species. Ecology, 79, pp. 2660-2675. BBC News, 2006, Britain blooms in extended summer, Available at http://news.bbc.co.uk/1/hi/uk/6066354.stm (accessed 30th October 2006). Case, M.J., Peterson, D.L. and David, L., 2007. Growth climate relations of Lodgepole Pine in the North Cascades National Park, Washington. Northwest Science, 81, pp. 62-75. Chen, X., Q., Hu, B. and Yu, R., 2005. Spatial and temporal variation of phenological growing season and climate change impacts in temperate eastern China. Global Change Biology, 11, pp. 1118-1130. Croxton, P.J., Huber, K., Collinson, N. and Sparks, T., 2006. How well do the Central England temperature and the England and Wales precipitation series represent the climate of the UK? International Journal of Climatology, 26, pp. 2287-2292. Curran, P.J. and Steele, C.M., 2005. MERIS: The re-branding of an ocean sensor. International Journal of Remote Sensing, 26, 1781-1798. Figure 5. The variation in MTCI and temperature in 2006 compared to the running average 2003 2005 for all sites. MTCI was compared to the CET mean. The monthly mean temperature was above average for the duration of the autumnal 2006 period 4. CONCLUSIONS The MTCI was sensitive to seasonal variations in inferred canopy chlorophyll content and could be a valuable tool for phenological studies. Higher than average autumnal temperatures in 2006 delayed observed senescence and resulted in an extended growing season at all four sites. There was a strong positive correlation between temperature and MTCI, from which we can infer that chlorophyll content was sensitive to temperature. Long term monitoring is required to determine whether climatic warming is leading to changing phenology and an increase in the growing season. However, a continual increase in mean monthly temperatures may not consequentially lead to an increase in canopy chlorophyll content as other growth conditions will become limiting. The use of the MERIS MTCI affords the ability to carry out long term monitoring. Curran, P.J., Dash, J. and Llewellyn, G.M., 2007. Indian Ocean tsunami: The use of MERIS (MTCI) data to infer salt stress in coastal vegetation. International Journal of Remote Sensing, 28, pp. 729-735. Dash, J. and Curran, P.J., 2004. The MERIS Terrestrial Chlorophyll Index. International Journal of Remote Sensing, 25, pp. 5003-5013. Dash, J. and Curran, P.J., 2007. Evaluation of the MERIS terrestrial chlorophyll index (MTCI). Advances in Space Research, 39, pp. 100-104. Delward, S., Preusker, R., Bourg, L., Santer, R., Ramon, D. and Fischer, J., 2007. MERIS in flight calibration. International Journal of Remote Sensing, 28, pp. 479-496. Deng, F.P., Su, G.L. and Liu, C., 2007. Seasonal variation of MODIS vegetation indexes and their statistical relationship with climate over the subtropic evergreen forest in Zhejiang, China. IEEE Geoscience and Remote Sensing Letters, 4, pp. 236-240. Houghton, J.T., Ding, Y., Griggs, D.J., Noguer, M., van der Linden, P.J., Dai, X., Maskell, K., Johnson, (eds) 2001. Intergovernmental Panel on Climate Change, The Scientific Basis. Contribution of Working Group I to the Third assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge UK.

Meteorological Office, 2006. Exceptionally warm extended summer 2006. Available at www.metoffice.gov.uk (accessed 18th October 2006). Menzel, A., Sparks, T. H., Estrella, N. and Roy, D.B., 2006, Altered geographic and temporal variability in phenology in response to climate change. Global Ecology and Biogeography, 15, pp. 498-504. Reed, B.C., 2006. Trend analysis of time series phenology of North America derived from satellite data. GIScience and Remote Sensing, 43, pp. 24-38. Rosenthal, S.I. and Camm, E.L., 1997. Photosynthetic decline and pigment loss during autumn foliar senescence in western larch (Larix occidentalis). Tree Physiology,17, pp. 767-775. Schleip, C., Menzel, A., Estrella, N. and Dose, V., 2006. The use of Bayesian analysis to detect recent changes in phenological events throughout the year. Agricultural and Forest Meteorology, 141, pp. 179-191. Sparks, T.H., Croxton, P.J., Collinson, N. and Taylor, P., 2005. Example of phonological change, past and present, in UK farming, Annals of Applied Biology, 146, pp. 531-537. Thomas, A., 2006, Climatic change and potential agricultural productivity in China, Erdkunde, 60, 157-172. van Oldenborgh, G.J. 2006, Extraordinarily mild European autumn 2006 due to global warming. Global Change Newsletter, 67, pp. 18-20. Zhang, X.Y., Friedl, M.A., Schaaf, C.B., 2006. Global vegetation phenology from Moderate Resolution Imaging Spectroradiometer (MODIS): Evaluation of global patterns and comparison with in situ measurements. Journal of Geophysical Research, 111: Article number G04017.