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

Toward the estimation of river discharge variations using MODIS data in ungauged basins

This study investigates the capability of the Moderate resolution Imaging Spectroradiometer (MODIS) to estimate river discharge, even for ungauged sites. Because of its frequent revisits (as little as every 3 h) and adequate spatial resolution (250 m), MODIS bands 1 and 2 have significant potential for mapping the extent of flooded areas and estimating river discharge even for medium-sized basins. Specifically, the different behaviour of water and land in the Near Infrared (NIR) portion of the electromagnetic spectrum is exploited by computing the ratio (C/M) of the MODIS channel 2 reflectance values between two pixels located within (M) and outside (C), but close to, the river. The values of C/M increase with the presence of water and, hence, with discharge. Moreover, in order to reduce the noise effects due to atmospheric contribution, an exponential smoothing filter is applied, thus obtaining C/M⁎. Time series of hourly mean flow velocity and discharge between 2005 and 2011 measured at four gauging stations located along the Po river (Northern Italy) are employed for testing the capability of C/M⁎ to estimate discharge/flow velocity. Specifically, the meanders and urban areas are considered the best locations for the position of the pixels M and C, respectively. Considering the optimal pixels, the agreement between C/M⁎ and discharge/flow velocity is fairly good with values in the range of 0.65–0.77. Additionally, the application to ungauged sites is tested by deriving a unique regional relationship between C/M⁎ and flow velocity valid for the whole Po river and providing only a slight deterioration of the performance. Finally, the sensitivity of the results to the selection of the C and M pixels is investigated by randomly changing their location. Also in this case, the agreement with in situ observations of velocity is fairly satisfactory (r ~ 0.6). The obtained results demonstrate the capability of MODIS to monitor discharge (and flow velocity). Therefore, its application for a larger number of sites worldwide will be the object of future studies

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

Assessing the role of EO in biodiversity monitoring: options for integrating in-situ observations with EO within the context of the EBONE concept

The European Biodiversity Observation Network (EBONE) is a European contribution on terrestrial monitoring to GEO BON, the Group on Earth Observations Biodiversity Observation Network. EBONE’s aims are to develop a system of biodiversity observation at regional, national and European levels by assessing existing approaches in terms of their validity and applicability starting in Europe, then expanding to regions in Africa. The objective of EBONE is to deliver: 1. A sound scientific basis for the production of statistical estimates of stock and change of key indicators; 2. The development of a system for estimating past changes and forecasting and testing policy options and management strategies for threatened ecosystems and species; 3. A proposal for a cost-effective biodiversity monitoring system. There is a consensus that Earth Observation (EO) has a role to play in monitoring biodiversity. With its capacity to observe detailed spatial patterns and variability across large areas at regular intervals, our instinct suggests that EO could deliver the type of spatial and temporal coverage that is beyond reach with in-situ efforts. Furthermore, when considering the emerging networks of in-situ observations, the prospect of enhancing the quality of the information whilst reducing cost through integration is compelling. This report gives a realistic assessment of the role of EO in biodiversity monitoring and the options for integrating in-situ observations with EO within the context of the EBONE concept (cfr. EBONE-ID1.4). The assessment is mainly based on a set of targeted pilot studies. Building on this assessment, the report then presents a series of recommendations on the best options for using EO in an effective, consistent and sustainable biodiversity monitoring scheme. The issues that we faced were many: 1. Integration can be interpreted in different ways. One possible interpretation is: the combined use of independent data sets to deliver a different but improved data set; another is: the use of one data set to complement another dataset. 2. The targeted improvement will vary with stakeholder group: some will seek for more efficiency, others for more reliable estimates (accuracy and/or precision); others for more detail in space and/or time or more of everything. 3. Integration requires a link between the datasets (EO and in-situ). The strength of the link between reflected electromagnetic radiation and the habitats and their biodiversity observed in-situ is function of many variables, for example: the spatial scale of the observations; timing of the observations; the adopted nomenclature for classification; the complexity of the landscape in terms of composition, spatial structure and the physical environment; the habitat and land cover types under consideration. 4. The type of the EO data available varies (function of e.g. budget, size and location of region, cloudiness, national and/or international investment in airborne campaigns or space technology) which determines its capability to deliver the required output. EO and in-situ could be combined in different ways, depending on the type of integration we wanted to achieve and the targeted improvement. We aimed for an improvement in accuracy (i.e. the reduction in error of our indicator estimate calculated for an environmental zone). Furthermore, EO would also provide the spatial patterns for correlated in-situ data. EBONE in its initial development, focused on three main indicators covering: (i) the extent and change of habitats of European interest in the context of a general habitat assessment; (ii) abundance and distribution of selected species (birds, butterflies and plants); and (iii) fragmentation of natural and semi-natural areas. For habitat extent, we decided that it did not matter how in-situ was integrated with EO as long as we could demonstrate that acceptable accuracies could be achieved and the precision could consistently be improved. The nomenclature used to map habitats in-situ was the General Habitat Classification. We considered the following options where the EO and in-situ play different roles: using in-situ samples to re-calibrate a habitat map independently derived from EO; improving the accuracy of in-situ sampled habitat statistics, by post-stratification with correlated EO data; and using in-situ samples to train the classification of EO data into habitat types where the EO data delivers full coverage or a larger number of samples. For some of the above cases we also considered the impact that the sampling strategy employed to deliver the samples would have on the accuracy and precision achieved. Restricted access to European wide species data prevented work on the indicator ‘abundance and distribution of species’. With respect to the indicator ‘fragmentation’, we investigated ways of delivering EO derived measures of habitat patterns that are meaningful to sampled in-situ observations

Remote Sensing Monitoring of Land Surface Temperature (LST)

This book is a collection of recent developments, methodologies, calibration and validation techniques, and applications of thermal remote sensing data and derived products from UAV-based, aerial, and satellite remote sensing. A set of 15 papers written by a total of 70 authors was selected for this book. The published papers cover a wide range of topics, which can be classified in five groups: algorithms, calibration and validation techniques, improvements in long-term consistency in satellite LST, downscaling of LST, and LST applications and land surface emissivity research

Continuous detection of small-scale changes in Scots pine dominated stands using dense Sentinel-2 time series

Climate change and severe extreme events, i.e., changes in precipitation and higher drought frequency, have a large impact on forests. In Poland, particularly Norway spruce and Scots pine forest stands are exposed to disturbances and have, thus experienced changes in recent years. Considering that Scots pine stands cover approximately 58% of forests in Poland, mapping these areas with an early and timely detection of forest cover changes is important, e.g., for forest management decisions. A cost-efficient way of monitoring forest changes is the use of remote sensing data from the Sentinel-2 satellites. They monitor the Earth’s surface with a high temporal (2–3 days), spatial (10–20 m), and spectral resolution, and thus, enable effective monitoring of vegetation. In this study, we used the dense time series of Sentinel-2 data from the years 2015–2019, (49 images in total), to detect changes in coniferous forest stands dominated by Scots pine. The simple approach was developed to analyze the spectral trajectories of all pixels, which were previously assigned to the probable forest change mask between 2015 and 2019. The spectral trajectories were calculated using the selected Sentinel-2 bands (visible red, red-edge 1–3, near-infrared 1, and short-wave infrared 1–2) and selected vegetation indices (Normalized Difference Moisture Index, Tasseled Cap Wetness, Moisture Stress Index, and Normalized Burn Ratio). Based on these, we calculated the breakpoints to determine when the forest change occurred. Then, a map of forest changes was created, based on the breakpoint dates. An accuracy assessment was performed for each detected date class using 861 points for 46 classes (45 dates and one class representing no changes detected). The results of our study showed that the short-wave infrared 1 band was the most useful for discriminating Scots pine forest stand changes, with the best overall accuracy of 75%. The evaluated vegetation indices underperformed single bands in detecting forest change dates. The presented approach is straightforward and might be useful in operational forest monitoring

Interannual Variations and Trends in Global Land Surface Phenology Derived from Enhanced Vegetation Index During 1982-2010

Land swiace phenology is widely retrieved from satellite observations at regional and global scales, and its long-term record has been demonstmted to be a valuable tool for reconstructing past climate variations, monitoring the dynamics of terrestrial ecosystems in response to climate impacts, and predicting biological responses to future climate scenarios. This srudy detected global land surface phenology from the advanced very high resolution radiometer (AVHRR) and the Moderate Resolution Imaging Spectroradiometer (MODIS) data from 1982 to 2010. Based on daily enhanced vegetation index at a spatial resolution of 0.05 degrees, we simulated the seasonal vegetative trajectory for each individual pixel using piecewise logistic models, which was then used to detect the onset of greenness increase (OGI) and the length of vegetation growing season (GSL). Further, both overall interannual variations and pixel-based trends were examIned across Koeppen's climate regions for the periods of 1982-1999 and 2000-2010, respectively. The results show that OGI and OSL varied considerably during 1982-2010 across the globe. Generally, the interarmual variation could be more than a month in precipitation-controlled tropical and dry climates while it was mainly less than 15 days in temperature-controlled temperate, cold, and polar climates. OGI, overall, shifted early, and GSL was prolonged from 1982 to 2010 in most climate regions in North America and Asia while the consistently significant trends only occurred in cold climate and polar climate in North America. The overall trends in Europe were generally insignificant. Over South America, late OGI was consistent (particularly from 1982 to 1999) while either positive or negative OSL trends in a climate region were mostly reversed between the periods of 1982-1999 and 2000-2010. In the Northern Hemisphere of Africa, OGI trends were mostly insignificant, but prolonged GSL was evident over individual climate regions during the last 3 decades. OGI mainly showed late trends in the Southern Hemisphere of Africa while GSL was reversed from reduced GSL trends (1982-1999) to prolonged trends (2000-2010). In Australia, GSL exhibited considerable interannual variation, but the consistent trend lacked presence in most regions. Finally, the proportion of pixels with significant trends was less than I% in most of climate regions although it could be as large as 10%