A Bayesian hierarchical model that quantifies long-term annual land surface phenology from sparse 30 m Landsat time series.