Land surface phenology from VEGETATION and PROBA-V data. Assessment over deciduous forests

Land surface phenology has been widely retrieved although no consensus exists on the optimal satellite dataset and the method to extract phenology metrics. This study is the first comprehensive comparison of vegetation variables and methods to retrieve land surface phenology for 1999-2017 time series of Copernicus Global Land products derived from SPOT-VEGETATION and PROBA-V data. We investigated the sensitivity of phenology to (I) the input vegetation variable: normalized difference vegetation index (NDVI), leaf area index (LAI), fraction of absorbed photosynthetically active radiation (FAPAR), and fraction of vegetation cover (FCOVER); (II) the smoothing and gap filling method for deriving seasonal trajectories; and (III) the method to extract phenological metrics: thresholds based on a percentile of the annual amplitude of the vegetation variable, autoregressive moving averages, logistic function fitting, and first derivative methods. We validated the derived satellite phenological metrics (start of the season (SoS) and end of the season (EoS)) using available ground observations of Betula pendula, B. alleghaniensis, Acer rubrum, Fagus grandifolia, and Quercus rubra in Europe (Pan-European PEP725 network) and the USA (National Phenology Network, USA-NPN). The threshold-based method applied to the smoothed and gap-filled LAI V2 time series agreed best with the ground phenology, with root mean square errors of ˜10 d and ˜25 d for the timing of SoS and EoS respectively. This research is expected to contribute for the operational retrieval of land surface phenology within the Copernicus Global Land Servic

Mapping Crop Cycles in China Using MODIS-EVI Time Series

As the Earth’s population continues to grow and demand for food increases, the need for improved and timely information related to the properties and dynamics of global agricultural systems is becoming increasingly important. Global land cover maps derived from satellite data provide indispensable information regarding the geographic distribution and areal extent of global croplands. However, land use information, such as cropping intensity (defined here as the number of cropping cycles per year), is not routinely available over large areas because mapping this information from remote sensing is challenging. In this study, we present a simple but efficient algorithm for automated mapping of cropping intensity based on data from NASA’s (NASA: The National Aeronautics and Space Administration) MODerate Resolution Imaging Spectroradiometer (MODIS). The proposed algorithm first applies an adaptive Savitzky-Golay filter to smooth Enhanced Vegetation Index (EVI) time series derived from MODIS surface reflectance data. It then uses an iterative moving-window methodology to identify cropping cycles from the smoothed EVI time series. Comparison of results from our algorithm with national survey data at both the provincial and prefectural level in China show that the algorithm provides estimates of gross sown area that agree well with inventory data. Accuracy assessment comparing visually interpreted time series with algorithm results for a random sample of agricultural areas in China indicates an overall accuracy of 91.0% for three classes defined based on the number of cycles observed in EVI time series. The algorithm therefore appears to provide a straightforward and efficient method for mapping cropping intensity from MODIS time series data

A Review of Vegetation Phenological Metrics Extraction Using Time-Series, Multispectral Satellite Data

Vegetation dynamics and phenology play an important role in inter-annual vegetation changes in terrestrial ecosystems and are key indicators of climate-vegetation interactions, land use/land cover changes, and variation in year-to-year vegetation productivity. Satellite remote sensing data have been widely used for vegetation phenology monitoring over large geographic domains using various types of observations and methods over the past several decades. The goal of this paper is to present a detailed review of existing methods for phenology detection and emerging new techniques based on the analysis of time-series, multispectral remote sensing imagery. This paper summarizes the objective and applications of detecting general vegetation phenology stages (e.g., green onset, time or peak greenness, and growing season length) often termed “land surface phenology,” as well as more advanced methods that estimate species-specific phenological stages (e.g., silking stage of maize). Common data-processing methods, such as data smoothing, applied to prepare the time-series remote sensing observations to be applied to phenological detection methods are presented. Specific land surface phenology detection methods as well as species-specific phenology detection methods based on multispectral satellite data are then discussed. The impact of different error sources in the data on remote-sensing based phenology detection are also discussed in detail, as well as ways to reduce these uncertainties and errors. Joint analysis of multiscale observations ranging from satellite to more recent ground-based sensors is helpful for us to understand satellite-based phenology detection mechanism and extent phenology detection to regional scale in the future. Finally, emerging opportunities to further advance remote sensing of phenology is presented that includes observations from Cubesats, near-surface observations such as PhenoCams, and image data fusion techniques to improve the spatial resolution of time-series image data sets needed for phenological characterization

A Novel Strategy to Reconstruct NDVI Time-Series with High Temporal Resolution from MODIS Multi-Temporal Composite Products

Vegetation indices (VIs) data derived from satellite imageries play a vital role in land surface vegetation and dynamic monitoring. Due to the excessive noises (e.g., cloud cover, atmospheric contamination) in daily VI data, temporal compositing methods are commonly used to produce composite data to minimize the negative influence of noise over a given compositing time interval. However, VI time series with high temporal resolution were preferred by many applications such as vegetation phenology and land change detections. This study presents a novel strategy named DAVIR-MUTCOP (DAily Vegetation Index Reconstruction based on MUlti-Temporal COmposite Products) method for normalized difference vegetation index (NDVI) time-series reconstruction with high temporal resolution. The core of the DAVIR-MUTCOP method is a combination of the advantages of both original daily and temporally composite products, and selecting more daily observations with high quality through the temporal variation of temporally corrected composite data. The DAVIR-MUTCOP method was applied to reconstruct high-quality NDVI time-series using MODIS multi-temporal products in two study areas in the continental United States (CONUS), i.e., three field experimental sites near Mead, Nebraska from 2001 to 2012 and forty-six AmeriFlux sites evenly distributed across CONUS from 2006 to 2010. In these two study areas, the DAVIR-MUTCOP method was also compared to several commonly used methods, i.e., the Harmonic Analysis of Time- Series (HANTS) method using original daily observations, Savitzky–Golay (SG) filtering using daily observations with cloud mask products as auxiliary data, and SG filtering using temporally corrected composite data. The results showed that the DAVIR-MUTCOP method significantly improved the temporal resolution of the reconstructed NDVI time series. It performed the best in reconstructing NDVI time-series across time and space (coefficient of determination (R2 = 0.93 ~ 0.94) between reconstructed NDVI and ground-observed LAI). DAVIR-MUTCOP method presented the highest robustness and accuracy with the change of the filtering parameter (R2 = 0.99 ~ 1.00, bias = 0.001, root mean square error (RMSE) = 0.020). Only MODIS data were used in this study; nevertheless, the DAVIR-MUTCOP method proposed a universal and potential way to reconstruct daily time series of other VIs or from other operational sensors, e.g., AVHRR and VIIRS

Detecting Inter-Annual Variations in the Phenology of Evergreen Conifers Using Long-Term MODIS Vegetation Index Time Series

Long-term observations of vegetation phenology can be used to monitor the response of terrestrial ecosystems to climate change. Satellite remote sensing provides the most efficient means to observe phenological events through time series analysis of vegetation indices such as the Normalized Difference Vegetation Index (NDVI). This study investigates the potential of a Photochemical Reflectance Index (PRI), which has been linked to vegetation light use efficiency, to improve the accuracy of MODIS-based estimates of phenology in an evergreen conifer forest. Timings of the start and end of the growing season (SGS and EGS) were derived from a 13-year-long time series of PRI and NDVI based on a MAIAC (multi-angle implementation of atmospheric correction) processed MODIS dataset and standard MODIS NDVI product data. The derived dates were validated with phenology estimates from ground-based flux tower measurements of ecosystem productivity. Significant correlations were found between the MAIAC time series and ground-estimated SGS (R-2 = 0.36-0.8), which is remarkable since previous studies have found it difficult to observe inter-annual phenological variations in evergreen vegetation from satellite data. The considerably noisier NDVI product could not accurately predict SGS, and EGS could not be derived successfully from any of the time series. While the strongest relationship overall was found between SGS derived from the ground data and PRI, MAIAC NDVI exhibited high correlations with SGS more consistently (R-2 > 0.6 in all cases). The results suggest that PRI can serve as an effective indicator of spring seasonal transitions, however, additional work is necessary to confirm the relationships observed and to further explore the usefulness of MODIS PRI for detecting phenology.Peer reviewe

Mapping the Spatial Distribution of Winter Crops at Sub-Pixel Level Using AVHRR NDVI Time Series and Neural Nets

For large areas, it is difficult to assess the spatial distribution and inter-annual variation of crop acreages through field surveys. Such information, however, is of great value for governments, land managers, planning authorities, commodity traders and environmental scientists. Time series of coarse resolution imagery offer the advantage of global coverage at low costs, and are therefore suitable for large-scale crop type mapping. Due to their coarse spatial resolution, however, the problem of mixed pixels has to be addressed. Traditional hard classification approaches cannot be applied because of sub-pixel heterogeneity. We evaluate neural networks as a modeling tool for sub-pixel crop acreage estimation. The proposed methodology is based on the assumption that different cover type proportions within coarse pixels prompt changes in time profiles of remotely sensed vegetation indices like the Normalized Difference Vegetation Index (NDVI). Neural networks can learn the relation between temporal NDVI signatures and the sought crop acreage information. This learning step permits a non-linear unmixing of the temporal information provided by coarse resolution satellite sensors. For assessing the feasibility and accuracy of the approach, a study region in central Italy (Tuscany) was selected. The task consisted of mapping the spatial distribution of winter crops abundances within 1 km AVHRR pixels between 1988 and 2001. Reference crop acreage information for network training and validation was derived from high resolution Thematic Mapper/Enhanced Thematic Mapper (TM/ETM+) images and official agricultural statistics. Encouraging results were obtained demonstrating the potential of the proposed approach. For example, the spatial distribution of winter crop acreage at sub-pixel level was mapped with a cross-validated coefficient of determination of 0.8 with respect to the reference information from high resolution imagery. For the eight years for which reference information was available, the root mean squared error (RMSE) of winter crop acreage at sub-pixel level was 10%. When combined with current and future sensors, such as MODIS and Sentinel-3, the unmixing of AVHRR data can help in the building of an extended time series of crop distributions and cropping patterns dating back to the 80s.JRC.H.4-Monitoring Agricultural Resource

A threshold method for robust and fast estimation of land-surface phenology using Google Earth Engine

Cloud-based platforms are changing the way of analyzing remotely sensed data by providing high computational power and rapid access to massive volumes of data. Several types of studies use cloud-based platforms for global-scale analyses, but the number of land-surface phenology (LSP) studies that use cloud-based platforms is low. We analyzed the performance of the state-of-the-art LSP algorithms and propose a new threshold-based method that we implemented in Google Earth Engine (GEE). This new LSP method, called maximum separation (MS) method, applies a moving window that estimates the ratio of observations that exceed a given threshold before and after the central day. The start and end of the growing season are the days of the year when the difference between the ratios before and after the central day are minimal and maximal. The MODIS phenology metrics estimated with the MS method showed similar performances as traditional threshold methods when compared with ground estimations derived from the PhenoCam dataset, a network of digital cameras that provides near-surface remotely sensed observations of vegetation phenology. The main advantage of the MS method is that it can be directly applied to daily nonsmoothed time series without any additional preprocessing steps. The implementation of the proposed method in GEE allowed the processing of global phenological maps derived from MODIS. The distribution of code in GEE allows the reproducibility of results and the rapid processing of LSP metrics by the scientific community

Doctor of Philosophy

dissertationCheatgrass, or Downy Brome, is an exotic winter annual weed native to the Mediterranean region. Since its introduction to the U.S., it has become a significant weed and aggressive invader of sagebrush, pinion-juniper, and other shrub communities, where it can completely out-compete native grasses and shrubs. In this research, remotely sensed data combined with field collected data are used to investigate the distribution of the cheatgrass in Central Utah, to characterize the trend of the NDVI time-series of cheatgrass, and to construct a spatially explicit populationbased model to simulate the spatial-temporal dynamics of the cheatgrass. This research proposes a method for mapping the canopy closure of invasive species using remotely sensed data acquired at different dates. Different invasive species have their own distinguished phenologies and the satellite images in different dates could be used to capture the phenology. The results of cheatgrass abundance prediction have a good fit with the field data for both linear regression and regression tree models, although the regression tree model has better performance than the linear regression model. To characterize the trend of NDVI time-series of cheatgrass, a novel smoothing algorithm named RMMEH is presented in this research to overcome some drawbacks of many other algorithms. By comparing the performance of RMMEH in smoothing a 16-day composite of the MODIS NDVI time-series with that of two other methods, which are the 4253EH, twice and the MVI, we have found that RMMEH not only keeps the original valid NDVI points, but also effectively removes the spurious spikes. The reconstructed NDVI time-series of different land covers are of higher quality and have smoother temporal trend. To simulate the spatial-temporal dynamics of cheatgrass, a spatially explicit population-based model is built applying remotely sensed data. The comparison between the model output and the ground truth of cheatgrass closure demonstrates that the model could successfully simulate the spatial-temporal dynamics of cheatgrass in a simple cheatgrass-dominant environment. The simulation of the functional response of different prescribed fire rates also shows that this model is helpful to answer management questions like, "What are the effects of prescribed fire to invasive species?" It demonstrates that a medium fire rate of 10% can successfully prevent cheatgrass invasion

Fernerkundung der Vegetationsphänologie über MODIS NDVI Daten - Herausforderungen bei der Datenverarbeitung und -validierung mittels Bodenbeobachtungen zahlreicher Arten und LiDAR

Phenology, the cyclic events in living organisms is triggered by climatic conditions and indicators of climate change. They are important factors influencing species interactions and ecosystem functioning. This thesis deals with the estimation of phenological metrics (Land Surface Phenology or LSP) from MODIS based time series NDVI data. Results of data analysis emphasises the role of ground observations, topography and LiDAR characteristics of forest stand in describing the variability in LSP.Phänologie, die zyklischen Stadien von lebenden Organismen werden über klimatische Verhältnisse gesteuert und dienen als Indikatoren des Klimawandels. Diese Faktoren beeinflussen maßgeblich die Interaktionen zwischen Arten und sind für das Funktionieren von Ökosystemen ausschlaggebend. Diese Arbeit behandelt die Bestimmung von phänologischen Metriken (Phänologie der Landoberfläche oder LSP) unter Verwendung von MODIS basierten NDVI Zeitreihen. Die Ergebnisse der Datenanalyse hebt die Wichtigkeit von Bodenbeobachtungen, Topographie und LiDAR Merkmalen von Waldbeständen bei der Beschreibung der LSP Variabilität hervor

The potential of satellite-observed crop phenology to enhance yield gap assessments in smallholder landscapes

Many of the undernourished people on the planet obtain their entitlements to food via agricultural-based livelihood strategies, often on underperforming croplands and smallholdings. In this context, expanding cropland extent is not a viable strategy for smallholders to meet their food needs. Therefore, attention must shift to increasing productivity on existing plots and ensuring yield gaps do not widen. Thus, supporting smallholder farmers to sustainably increase the productivity of their lands is one part of a complex solution to realizing universal food security. However, the information (e.g., location and causes of cropland underperformance) required to support measures to close yield gaps in smallholder landscapes are often not available. This paper reviews the potential of crop phenology, observed from satellites carrying remote sensing sensors, to fill this information gap. It is suggested that on a theoretical level phenological approaches can reveal greater intra-cropland thematic detail, and increase the accuracy of crop extent maps and crop yield estimates. However, on a practical level the spatial mismatch between the resolution at which crop phenology can be estimated from satellite remote sensing data and the scale of yield variability in smallholder croplands inhibits its use in this context. Similarly, the spatial coverage of remote sensing-derived phenology offers potential for integration with ancillary spatial datasets to identify causes of yield gaps. To reflect the complexity of smallholder cropping systems requires ancillary datasets at fine spatial resolutions which, often, are not available. This further precludes the use of crop phenology in attempts to unpick the causes of yield gaps. Research agendas should focus on generating fine spatial resolution crop phenology, either via data fusion or through new sensors (e.g., Sentinel-2) in smallholder croplands. This has potential to transform the applied use of remote sensing in this context