Jan Verbesselt1, Rob J Hyndman2, Achim Zeilis3, Darius Culvenor1

  1. Remote sensing team, CSIRO Sustainable Ecosystems, Private Bag 10, Melbourne VIC 3169, Australia
  2. Department of Econometrics and Business Statistics, Monash University, Melbourne VIC 3800, Australia
  3. Institute for Statistics, Leopold-Franzens-Universitt Innsbruck, 6020 Innsbruck, Austria

Remote Sensing of Environment, 114(12), 2970-2980.

A challenge in phenology studies is understanding what constitutes significant phenological change amidst background variation (e.g. noise) and ecosystem disturbances (e.g. fires). The majority of phenological studies have focussed on extracting critical points in the seasonal growth cycle (e.g. Start-of-spring), without exploiting the full temporal detail. Moreover, the high degree of phenological variability between years demonstrates the necessity of distinguishing long term phenological change from temporal variability. Here, we evaluate the phenological change detection ability of a method for detecting Breaks For Additive Seasonal and Trend (BFAST). BFAST integrates the decomposition of time series into trend, seasonal, and noise components with methods for detecting change within time series. BFAST detects significant phenological changes within time series by exploiting the full time series without needing to derive phenological metrics. The times and numbers of trend and phenological changes are iteratively estimated by fitting piecewise robust linear models, of which the parameters are used to characterize change by its magnitude and direction. We tested BFAST by simulating 16-day Normalized Difference Vegetation Index (NDVI) time series with varying amounts of seasonality and noise, containing abrupt disturbances (e.g. fires) and long term phenological changes. This revealed that BFAST is able to accurately detect the number and timing of phenological changes within time series while accounting for disturbances (e.g. fires) and noise. The simulation study also showed that the phenological change detection is influenced by the signal to noise ratio of the time series. Application of the method on 16-day NDVI MODIS images from 2000 until 2009 for a forested study area in south eastern Australia confirmed these results. Phenological change is more easily detected in grasslands where the seasonal amplitude is larger than 0.3 NDVI when compared to evergreen forests where the seasonal amplitude is approximately 0.1 NDVI while noise levels were the same. BFAST present a novel approach for the detection of significant long term phenological changes within full time series which is necessary to study spatio-temporal patterns in land cover phenology, distinguish change from interannual variability in a global change context. The method can be applied to other disciplines dealing with seasonal time series data, such as biology, hydrology, and climatology to detect and characterize change within time series. The methods described in this study are available in the BFAST package for R (R Development Core Team, 2009).

Keywords: seasonal change, phenology, change detection, time series, disturbance,climate change, remote sensing, NDVI, MODIS.

Working paper

Online paper


1 2 3 5
March 7th, 2017

Coherent Probabilistic Forecasts for Hierarchical Time Series

February 28th, 2017

Forecasting with temporal hierarchies

February 14th, 2017

The Australian Macro Database: An online resource for macroeconomic research in Australia

February 14th, 2017

Macroeconomic forecasting for Australia using a large number of predictors

January 13th, 2017

Visualising forecasting algorithm performance using time series instance spaces

January 1st, 2017

Associations between outdoor fungal spores and childhood and adolescent asthma hospitalisations

January 1st, 2017

Grouped functional time series forecasting: an application to age-specific mortality rates

December 7th, 2016

Exploring the influence of short-term temperature patterns on temperature-related mortality: a case-study of Melbourne, Australia

October 13th, 2016

Reconciling forecasts: the hts and thief packages

September 20th, 2016

smoothAPC package for R

September 20th, 2016

stR package for R

September 15th, 2016

Forecasting large collections of related time series

August 22nd, 2016

thief package for R

June 21st, 2016

Exploring time series collections used for forecast evaluation

May 25th, 2016

ISCRR time series workshop

May 6th, 2016

Automatic foRecasting using R

February 29th, 2016

On sampling methods for costly multi-objective black-box optimization

February 19th, 2016

Dynamic Algorithm Selection for Pareto Optimal Set Approximation

February 4th, 2016

Forecasting uncertainty in electricity smart meter data by boosting additive quantile regression

January 30th, 2016

Bayesian rank selection in multivariate regressions

January 25th, 2016

Probabilistic Energy Forecasting: Global Energy Forecasting Competition 2014 and Beyond

January 24th, 2016

Long-term forecasts of age-specific participation rates with functional data models

January 1st, 2016

Bagging exponential smoothing methods using STL decomposition and Box-Cox transformation

January 1st, 2016

Fast computation of reconciled forecasts for hierarchical and grouped time series

December 31st, 2015

Measuring forecast accuracy

November 26th, 2015

Forecasting hierarchical and grouped time series through trace minimization

November 2nd, 2015

Forecasting big time series data using R

October 7th, 2015

Optimal forecast reconciliation for big time series data

October 5th, 2015

Google workshop: Forecasting and visualizing big time series data

September 16th, 2015


August 25th, 2015

New IJF editors

August 17th, 2015

Machine learning bootcamp

August 7th, 2015

Statistical issues with using herbarium data for the estimation of invasion lag-phases

June 30th, 2015

Exploring the feature space of large collections of time series

June 26th, 2015

Exploring the boundaries of predictability: what can we forecast, and when should we give up?

June 25th, 2015

Automatic algorithms for time series forecasting

June 23rd, 2015

MEFM: An R package for long-term probabilistic forecasting of electricity demand

June 19th, 2015

Probabilistic forecasting of peak electricity demand

June 10th, 2015

Do human rhinovirus infections and food allergy modify grass pollen–induced asthma hospital admissions in children?

June 8th, 2015

STR: A Seasonal-Trend Decomposition Procedure Based on Regression

June 4th, 2015

Probabilistic time series forecasting with boosted additive models: an application to smart meter data

June 1st, 2015

Large-scale unusual time series detection

May 26th, 2015

Visualization of big time series data

May 22nd, 2015

Probabilistic forecasting of long-term peak electricity demand

April 20th, 2015

A note on the validity of cross-validation for evaluating time series prediction

April 4th, 2015

Discussion of “High-dimensional autocovariance matrices and optimal linear prediction”

April 1st, 2015

Change to the IJF editors

February 23rd, 2015

Visualization and forecasting of big time series data

January 12th, 2015

Visualizing and forecasting big time series data

December 24th, 2014

Bivariate data with ridges: two-dimensional smoothing of mortality rates