Classification of Satellite Time Series-derived Land Surface Phenology Focused on the Northern Fertile Crescent

Land surface phenology describes events in a seasonal vegetation cycle and can be used in a variety of applications from predicting onset of future drought conditions, to revealing potential limits of historical dry farming, to guiding more accurate dating of archeological sites. Traditional methods of monitoring vegetation phenology use data collected in situ. However, vegetation health indices derived from satellite remote sensor data, such as the normalized difference vegetation index (NDVI), have been used as proxy for vegetation phenology due to their repeated acquisition and broad area coverage. Land surface phenology is accessible in the NDVI satellite record when images are processed to be intercomparable over time and temporally ordered to create a time series. This study utilized NDVI time series to classify areas of similar vegetation phenology in the northern Fertile Crescent, an area from the middle Mediterranean coast to southern/south-eastern Turkey to western Iran and northern Iraq. Phenological monitoring of the northern Fertile Crescent is critical due to the area\u27s minimal water resources, susceptibility to drought, and understanding ancient historical reliance on precipitation for subsistence dry farming. Delineation of phenological classes provides areal and temporal synopsis of vegetation productivity time series. Phenological classes were developed from NDVI time series calculated from NOAA\u27s Advanced Very High Resolution Radiometer (AVHRR) imagery with 8 × 8 km spatial resolution over twenty-five years, and by NASA\u27s Moderate Resolution Imaging Spectroradiometer (MODIS) with 250 × 250 m spatial resolution over twelve years. Both AVHRR and MODIS time series were subjected to data reduction techniques in spatial and temporal dimensions. Optimized ISODATA clusters were developed for both of these data reduction techniques in order to compare the effects of spatial versus temporal aggregation. Within the northern Fertile Crescent study area, the spatial reduction technique showed increased cluster cohesion over the temporal reduction method. The latter technique showed an increase in temporal smoothing over the spatial reduction technique. Each technique has advantages depending on the desired spatial or temporal granularity. Additional work is required to determine optimal scale size for the spatial data reduction technique

The CACAO Method for Smoothing, Gap Filling, and Characterizing Seasonal Anomalies in Satellite Time Series

Consistent, continuous, and long time series of global biophysical variables derived from satellite data are required for global change research. A novel climatology fitting approach called CACAO (Consistent Adjustment of the Climatology to Actual Observations) is proposed to reduce noise and fill gaps in time series by scaling and shifting the seasonal climatological patterns to the actual observations. The shift and scale CACAO parameters adjusted for each season allow quantifying shifts in the timing of seasonal phenology and inter-annual variations in magnitude as compared to the average climatology. CACAO was assessed first over simulated daily Leaf Area Index (LAI) time series with varying fractions of missing data and noise. Then, performances were analyzed over actual satellite LAI products derived from AVHRR Long-Term Data Record for the 1981-2000 period over the BELMANIP2 globally representative sample of sites. Comparison with two widely used temporal filtering methods-the asymmetric Gaussian (AG) model and the Savitzky-Golay (SG) filter as implemented in TIMESAT-revealed that CACAO achieved better performances for smoothing AVHRR time series characterized by high level of noise and frequent missing observations. The resulting smoothed time series captures well the vegetation dynamics and shows no gaps as compared to the 50-60% of still missing data after AG or SG reconstructions. Results of simulation experiments as well as confrontation with actual AVHRR time series indicate that the proposed CACAO method is more robust to noise and missing data than AG and SG methods for phenology extraction

Operational Large-Area Land-Cover Mapping: An Ethiopia Case Study

Knowledge of land cover and land use nationally is a prerequisite of many studies on drivers of land change, impacts on climate, carbon storage and other ecosystem services, and allows for sufficient planning and management. Despite this, many regions globally do not have accurate and consistent coverage at the national scale. This is certainly true for Ethiopia. Large-area land-cover characterization (LALCC), at a national scale is thus an essential first step in many studies of land-cover change, and yet is itself problematic. Such LALCC based on remote-sensing image classification is associated with a spectrum of technical challenges such as data availability, radiometric inconsistencies within/between images, and big data processing. Radiometric inconsistencies could be exacerbated for areas, such as Ethiopia, with a high frequency of cloud cover, diverse ecosystem and climate patterns, and large variations in elevation and topography. Obtaining explanatory variables that are more robust can improve classification accuracy. To create a base map for the future study of large-scale agricultural land transactions, we produced a recent land-cover map of Ethiopia. Of key importance was the creation of a methodology that was accurate and repeatable and, as such, could be used to create earlier, comparable land-cover classifications in the future for the same region. We examined the effects of band normalization and different time-series image compositing methods on classification accuracy. Both top of atmosphere and surface reflectance products from the Landsat 8 Operational Land Imager (OLI) were tested for single-time classification independently, where the latter resulted in 1.1% greater classification overall accuracy. Substitution of the original spectral bands with normalized difference spectral indices resulted in an additional improvement of 1.0% in overall accuracy. Three approaches for multi-temporal image compositing, using Landsat 8 OLI and Moderate Resolution Imaging Spectroradiometer (MODIS) data, were tested including sequential compositing, i.e., per-pixel summary measures based on predefined periods, probability density function compositing, i.e., per-pixel characterization of distribution of spectral values, and per-pixel sinusoidal models. Multi-temporal composites improved classification overall accuracy up to 4.1%, with respect to single-time classification with an advantage of the Landsat OLI-driven composites over MODIS-driven composites. Additionally, night-time light and elevation data were used to improve the classification. The elevation data and its derivatives improved classification accuracy by 1.7%. The night-time light data improve producer’s accuracy of the Urban/Built class with the cost of decreasing its user’s accuracy. Results from this research can aid map producers with decisions related to operational large-area land-cover mapping, especially with selecting input explanatory variables and multi-temporal image compositing, to allow for the creation of accurate and repeatable national-level land-cover products in a timely fashion

Evaluation of the impact of climate and human induced changes on the Nigerian forest using remote sensing

The majority of the impact of climate and human induced changes on forest are related to climate variability and deforestation. Similarly, changes in forest phenology due to climate variability and deforestation has been recognized as being among the most important early indicators of the impact of environmental change on forest ecosystem functioning. Comprehensive data on baseline forest cover changes including deforestation is required to provide background information needed for governments to make decision on Reducing Emissions from Deforestation and Forest Degradation (REED). Despite the fact that Nigeria ranks among the countries with highest deforestation rates based on Food and Agricultural Organization estimates, only a few studies have aimed at mapping forest cover changes at country scales. However, recent attempts to map baseline forest cover and deforestation in Nigeria has been based on global scale remote sensing techniques which do not confirm with ground based observations at country level. The aim of this study is two-fold: firstly, baseline forest cover was estimated using an ‘adaptive’ remote sensing model that classified forest cover with high accuracies at country level for the savanna and rainforest zones. The first part of this study also compared the potentials of different MODIS data in detecting forest cover changes at regional (cluster level) scale. The second part of this study explores the trends and response of forest phenology to rainfall across four forest clusters from 2002 to 2012 using vegetation index data from the MODIS and rainfall data obtained from the TRMM.Tertiary Education Trust Fund, Nigeri

REMOTE SENSING DATA ANALYSIS FOR ENVIRONMENTAL AND HUMANITARIAN PURPOSES. The automation of information extraction from free satellite data.

This work is aimed at investigating technical possibilities to provide information on environmental parameters that can be used for risk management. The World food Program (WFP) is the United Nations Agency which is involved in risk management for fighting hunger in least-developed and low-income countries, where victims of natural and manmade disasters, refugees, displaced people and the hungry poor suffer from severe food shortages. Risk management includes three different phases (pre-disaster, response and post disaster) to be managed through different activities and actions. Pre disaster activities are meant to develop and deliver risk assessment, establish prevention actions and prepare the operative structures for managing an eventual emergency or disaster. In response and post disaster phase actions planned in the pre-disaster phase are executed focusing on saving lives and secondly, on social economic recovery. In order to optimally manage its operations in the response and post disaster phases, WFP needs to know, in order to estimate the impact an event will have on future food security as soon as possible, the areas affected by the natural disaster, the number of affected people, and the effects that the event can cause to vegetation. For this, providing easy-to-consult thematic maps about the affected areas and population, with adequate spatial resolution, time frequency and regular updating can result determining. Satellite remote sensed data have increasingly been used in the last decades in order to provide updated information about land surface with an acceptable time frequency. Furthermore, satellite images can be managed by automatic procedures in order to extract synthetic information about the ground condition in a very short time and can be easily shared in the web. The work of thesis, focused on the analysis and processing of satellite data, was carried out in cooperation with the association ITHACA (Information Technology for Humanitarian Assistance, Cooperation and Action), a center of research which works in cooperation with the WFP in order to provide IT products and tools for the management of food emergencies caused by natural disasters. These products should be able to facilitate the forecasting of the effects of catastrophic events, the estimation of the extension and location of the areas hit by the event, of the affected population and thereby the planning of interventions on the area that could be affected by food insecurity. The requested features of the instruments are: • Regular updating • Spatial resolution suitable for a synoptic analysis • Low cost • Easy consultation Ithaca is developing different activities to provide georeferenced thematic data to WFP users, such a spatial data infrastructure for storing, querying and manipulating large amounts of global geographic information, and for sharing it between a large and differentiated community; a system of early warning for floods, a drought monitoring tool, procedures for rapid mapping in the response phase in a case of natural disaster, web GIS tools to distribute and share georeferenced information, that can be consulted only by means of a web browser. The work of thesis is aimed at providing applications for the automatic production of base georeferenced thematic data, by using free global satellite data, which have characteristics suitable for analysis at a regional scale. In particular the main themes of the applications are water bodies and vegetation phenology. The first application aims at providing procedures for the automatic extraction of water bodies and will lead to the creation and update of an historical archive, which can be analyzed in order to catch the seasonality of water bodies and delineate scenarios of historical flooded areas. The automatic extraction of phenological parameters from satellite data will allow to integrate the existing drought monitoring system with information on vegetation seasonality and to provide further information for the evaluation of food insecurity in the post disaster phase. In the thesis are described the activities carried on for the development of procedures for the automatic processing of free satellite data in order to produce customized layers according to the exigencies in format and distribution of the final users. The main activities, which focused on the development of an automated procedure for the extraction of flooded areas, include the research of an algorithm for the classification of water bodies from satellite data, an important theme in the field of management of the emergencies due to flood events. Two main technologies are generally used: active sensors (radar) and passive sensors (optical data). Advantages for active sensors include the ability to obtain measurements anytime, regardless of the time of day or season, while passive sensors can only be used in the daytime cloud free conditions. Even if with radar technologies is possible to get information on the ground in all weather conditions, it is not possible to use radar data to obtain a continuous archive of flooded areas, because of the lack of a predetermined frequency in the acquisition of the images. For this reason the choice of the dataset went in favor of MODIS (Moderate Resolution Imaging Spectroradiometer), optical data with a daily frequency, a spatial resolution of 250 meters and an historical archive of 10 years. The presence of cloud coverage prevents from the acquisition of the earth surface, and the shadows due to clouds can be wrongly classified as water bodies because of the spectral response very similar to the one of water. After an analysis of the state of the art of the algorithms of automated classification of water bodies in images derived from optical sensors, the author developed an algorithm that allows to classify the data of reflectivity and to temporally composite them in order to obtain flooded areas scenarios for each event. This procedure was tested in the Bangladesh areas, providing encouraging classification accuracies. For the vegetation theme, the main activities performed, here described, include the review of the existing methodologies for phenological studies and the automation of the data flow between inputs and outputs with the use of different global free satellite datasets. In literature, many studies demonstrated the utility of the NDVI (Normalized Difference Vegetation Index) indices for the monitoring of vegetation dynamics, in the study of cultivations, and for the survey of the vegetation water stress. The author developed a procedure for creating layers of phenological parameters which integrates the TIMESAT software, produced by Lars Eklundh and Per Jönsson, for processing NDVI indices derived from different satellite sensors: MODIS (Moderate Resolution Imaging Spectroradiometer), AVHRR (Advanced Very High Resolution Radiometer) AND SPOT (Système Pour l'Observation de la Terre) VEGETATION. The automated procedure starts from data downloading, calls in a batch mode the software and provides customized layers of phenological parameters such as the starting of the season or length of the season and many others

The use of remote sensing data for drought assessment and monitoring in southwest Asia

Drought / Monitoring / Indicators / Assessment / Remote sensing / Asia

Improved leaf area index estimation by considering both temporal and spatial variations

Variations in Leaf Area Index (LAI) can greatly alter output values and patterns of various models that deal with energy flux exchange between the land surface and the atmosphere. Customarily, such models are initiated by LAI estimated from satellite-level Vegetation Indices (VIs) including routinely produced Normalized Difference Vegetation Index (NDVI) products. However, the accuracy from LAI-VI relationships greatly varies due to many factors, including temporal and spatial variations in LAI and a selected VI. In addition, NDVI products derived from various sensors have demonstrated variations in a certain degree on describing temporal and spatial variations in LAI, especially in semi-arid areas. This thesis therefore has three objectives: 1) determine a suitable VI for quantifying LAI temporal variation; 2) improve LAI estimation by considering both temporal and spatial variations in LAI; and 3) evaluate routinely produced NDVI products on monitoring temporal and spatial variations in LAI. The study site was set up in conserved semi-arid mixed grassland in St. Denis, Saskatchewan, Canada. One 600 m - long sampling transect was set up across the rolling typography, and six plots with a size of 40 × 40 m each were randomly designed and each was in a relatively homogenous area. Plant Area Index (PAI, which was validated to obtain LAI), ground hyperspectral reflectance, ground covers (grasses, forbs, standing dead, litter, and bare soil), and soil moisture data were collected over the sampling transect and plots from May through September, 2008. Satellite data used are SPOT 4/5 images and 16-day Moderate Resolution Imaging Spectroradiometer (MODIS) 250m, 1km as well as 10-day SPOT-vegetation (SPOT-VGT) NDVI products from May to October, 2007 and 2008. The results show that NDVI is the most suitable VI for quantifying temporal variation of LAI. LAI estimation is much improved by considering both temporal and spatial variations. Based on the ground reflectance data, the r2 value is increased by 0.05, 0.31, and 0.23 and an averaged relative error is decreased by 1.57, 1.62, and 0.67 in the early, maximum, and late growing season, respectively. MODIS 250m NDVI products are the most useful datasets and MODIS 1km NDVI products are superior to SPOT-VGT 1km composites for monitoring intra-annual spatiotemporal variations in LAI. The proposed LAI estimation approach can be used in other studies to obtain more accurate LAI, and thus this research will be beneficial for grassland modeling