Rasmus Fensholt

MODIS leaf area index products: from validation to algorithm improvement

Global products of vegetation green Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation absorbed by vegetation (FPAR) are being operationally produced from Terra and Aqua Moderate Resolution Imaging Spectroradiometers (MODIS) at 1-km resolution and eight-day frequency. This paper summarizes the experience of several collaborating investigators on validation of MODIS LAI products and demonstrates the close connection between product validation and algorithm refinement activities. The validation of moderate resolution LAI products includes three steps: 1) field sampling representative of LAI spatial distribution and dynamic range within each major land cover type at the validation site; 2) development of a transfer function between field LAI measurements and high resolution satellite data to generate a reference LAI map over an extended area; and 3) comparison of MODIS LAI with aggregated reference LAI map at patch (multipixel) scale in view of geo-location and pixel shift uncertainties. The MODIS LAI validation experiences, summarized here, suggest three key factors that influence the accuracy of LAI retrievals: 1) uncertainties in input land cover data, 2) uncertainties in input surface reflectances, and 3) uncertainties from the model used to build the look-up tables accompanying the algorithm. This strategy of validation efforts guiding algorithm refinements has led to progressively more accurate LAI products from the MODIS sensors aboard NASAs Terra and Aqua platforms

Evaluation of fraction of absorbed photosynthetically active radiation products for different canopy radiation transfer regimes: Methodology and results using Joint Research Center products derived from SeaWiFS against ground-based estimations

[1] This paper discusses the quality and the accuracy of the Joint Research Center (JRC) fraction of absorbed photosynthetically active radiation (FAPAR) products generated from an analysis of Sea-viewing Wide Field-of-view Sensor (SeaWiFS) data. The FAPAR value acts as an indicator of the presence and state of the vegetation and it can be estimated from remote sensing measurements using a physically based approach. The quality of the SeaWiFS FAPAR products assessed in this paper capitalizes on the availability of a 6-year FAPAR time series over the full globe. This evaluation exercise is performed in two phases involving, first, an analysis of the verisimilitude of the FAPAR products under documented environmental conditions and, second, a direct comparison of the FAPAR values with ground-based estimations where and when the latter are available. This second phase is conducted following a careful analysis of problems arising for performing such a comparison. This results in the grouping of available field information into broad categories representing different radiative transfer regimes. This strategy greatly helps the interpretation of the results since it recognizes the various levels of difficulty and sources of uncertainty associated with the radiative sampling of different types of vegetation canopies.

Evaluation of MODIS and NOAA AVHRR vegetation indices with in situ measurements in a semi‐arid environment

Much effort has been made in recent years to improve the spectral and spatial resolution of satellite sensors to develop improved vegetation indices reflecting surface conditions. In this study satellite vegetation indices from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Advanced Very High Resolution Radiometer (AVHRR) are evaluated against two years of in situ measurements of vegetation indices in Senegal. The in situ measurements are obtained using four masts equipped with self‐registrating multispectral radiometers designed for the same wavelengths as the satellite sensor channels. In situ measurements of the MODIS Normalized Difference Vegetation Index (NDVI) and AVHRR NDVI are equally sensitive to vegetation; however, the MODIS NDVI is consistently higher than the AVHRR NDVI. The MODIS Enhanced Vegetation Index (EVI) proved more sensitive to dense vegetation than both AVHRR NDVI and MODIS NDVI. EVI and NDVI based on the MODIS 16‐day constrained view angle maximum value composite (CV‐MVC) product captured the seasonal dynamics of the field observations satisfactorily but a standard 16‐day MVC product estimated from the daily MODIS surface reflectance data without view angle constraints yielded higher correlations between the satellite indices and field measurements (R 2 values ranging from 0.74 to 0.98). The standard MVC regressions furthermore approach a 1 : 1 line with in situ measured values compared to the CV‐MVC regressions. The 16‐day MVC AVHRR data did not satisfactorily reflect the variation in the in situ data. Seasonal variation in the in situ measurements is captured reasonably with R 2 values of 0.75 in 2001 and 0.64 in 2002, but the dynamic range of the AVHRR satellite data is very low—about a third to a half of the values from in situ measurements. Consequently the in situ vegetation indices were emulated much better by the MODIS indices than by the AVHRR NDVI.

Earth observation of vegetation status in the Sahelian and Sudanian West Africa - comparison of TERRA MODIS and NOAA AVHRR satellite data

The first year of Moderate Resolution Imaging Spectroradiometer (MODIS) data are compared with National Oceanic and Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVHRR) data for derivation of biophysical variables in Senegal, West Africa. The dynamic range of the two MODIS vegetation indices (VIs)—the continuity vegetation index (CVI) and the enhanced vegetation index (EVI)—is generally much larger than for the NOAA AVHRR normalized difference vegetation index (NDVI) data, indicating the importance of the change in near-infrared wavelength configuration from the NOAA AVHRR sensor to the MODIS sensor. Senegal is characterized by a pronounced gradient in the vegetation density covering a range of agro-climatic zones from arid to humid and it is found that the MODIS CVI values saturate for high VI values while the EVI demonstrates improved sensitivity for high biomass. Compared to NOAA AVHRR the MODIS VIs generally correlate better to the MODIS fraction of absorbed photosynthetically active radiation (fAPAR) absorbed by vegetation canopies and the leaf area index (LAI; the one-sided green leaf area per unit ground area). CVI is found to correlate better to both fAPAR and LAI than is the case for EVI because of the larger dynamic range of the CVI data. This suggests that the problem of background contamination on VIs from soil is not as severe in Senegal as has been found in other semi-arid African areas.

Evaluation of AVHRR PAL and GIMMS 10‐day composite NDVI time series products using SPOT‐4 vegetation data for the African continent

Global 8 km resolution AVHRR (advanced very high resolution radiometer) NDVI (normalized difference vegetation index) 10‐day composite data sets have been used for numerous local to global scale vegetation time series studies during recent years. AVHRR Pathfinder (PAL) NDVI was available from 1981 until 2001, and the new AVHRR GIMMS NDVI was available from 1981 to the present time. A number of aspects potentially introduce noise in the NDVI data set due to the AVHRR sensor design and data processing. NDVI from SPOT‐4 VGT data is considered an improvement over AVHRR, and for this reason it is important to examine how and if the differences in sensor design and processing influence continental scale NDVI composite products. In this study, the quality of these AVHRR NDVI time series are evaluated by the continental scale 1 km resolution SPOT‐4 vegetation (VGT) 10‐day composite (S10) NDVI data. Three years of AVHRR PAL (1998–2000) and seven years of GIMMS (1998–2004) have been compared to 8 km resampled SPOT‐4 VGT (1998–2004) data. The dynamic range of SPOT‐4 VGT NDVI tends to be higher than the AVHRR PAL NDVI, whereas there is an exact match between AVHRR GIMMS NDVI and SPOT‐4 VGT NDVI. Ortho‐regression analysis on annually integrated values of AVHRR PAL/GIMMS and SPOT‐4 VGT on a continental scale reveals high correlations amongst the AVHRR and the SPOT data set, with lowest RMSE (root mean square error) on the GIMMS/SPOT‐4 VGT compared to the PAL/SPOT‐4 VGT. Analyses on decade data likewise show that a linear relation exists between Spot‐4 VGT NDVI and the two AVHRR composite products; GIMMS explaining most of the Spot‐4 VGT NDVI variance compared to PAL. These results show that the AVHRR GIMMS NDVI is more consistent with Spot‐4 VGT NDVI compared to AVHRR PAL versus Spot‐4 VGT NDVI (in terms of RMSE and dynamic range) and can therefore be considered the more accurate long time AVHRR data record. Analyses performed on monthly maximum composites and decade composite data, however, reveal intra‐annual variations in the correlation between SPOT‐4 VGT and the two AVHRR data sets, which are attributed to different cloud masking algorithms. The SPOT‐4 VGT cloud‐screening algorithm is insufficient, thereby suppressing the rainy season NDVI.

Assimilation of SMOS‐derived soil moisture in a fully integrated hydrological and soil‐vegetation‐atmosphere transfer model in Western Denmark

Real surface soil moisture retrieved from the Soil Moisture and Ocean Salinity (SMOS) satellite is downscaled and assimilated in a fully integrated hydrological and soil-vegetation-atmosphere transfer (MIKE SHE SW-ET) model using a bias aware ensemble transform Kalman filter (Bias-ETKF). Satellite-derived soil moisture assimilation in a catchment scale model is typically restricted by two challenges: (1) passive microwave is too coarse for direct assimilation and (2) the data tend to be biased. The solution proposed in this study is to disaggregate the SMOS bias using a higher resolution land cover classification map that was derived from Landsat thermal images. Using known correlations between SMOS bias and vegetation type, the assimilation filter is adapted to calculate biases online, using an initial bias estimate. Real SMOS-derived soil moisture is assimilated in a precalibrated catchment model in Denmark. The objective is to determine if any additional gains can be achieved by SMOS surface soil moisture assimilation beyond the optimized model. A series of assimilation experiments were designed to (1) determine how effectively soil moisture corrections propagate downward in the soil column, (2) compare the efficacy of in situ versus SMOS assimilation, and (3) determine how soil moisture assimilation affects fluxes and discharge in the catchment. We find that assimilation of SMOS improved R2 soil moisture correlations in the upper 5 cm compared to a network of 30 in situ sensors for most land cover classes. Assimilation also brought modest gains in R2 at 25 cm depth but slightly degraded the correlation at 50 cm depth. Assimilation overcorrected discharge peaks.

Mapping and Evaluation of NDVI Trends from Synthetic Time Series Obtained by Blending Landsat and MODIS Data around a Coalfield on the Loess Plateau

The increasingly intensive and extensive coal mining activities on the Loess Plateau pose a threat to the fragile local ecosystems. Quantifying the effects of coal mining activities on environmental conditions is of great interest for restoring and managing the local ecosystems and resources. This paper generates dense NDVI (Normalized Difference Vegetation Index) time series between 2000 and 2011 at a spatial resolution of 30 m by blending Landsat and MODIS (Moderate Resolution Imaging Spectroradiometer) data using the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and further evaluates its capability for mapping vegetation trends around a typical coalfield on the Loss Plateau. Synthetic NDVI images were generated using (1) STARFM-generated NIR (near infrared) and red band reflectance data (scheme 1) and (2) Landsat and MODIS NDVI images directly as inputs for STARFM (scheme 2). By comparing the synthetic NDVI images with the corresponding Landsat NDVI, we found that scheme 2 consistently generated better results (0.70 < R2 < 0.76) than scheme 1 (0.56 < R2 < 0.70) in this study area. Trend analysis was then performed with the synthetic dense NDVI time series and the annual maximum NDVI (NDVImax) time series. The accuracy of these trends was evaluated by comparing to those from the corresponding MODIS time series, and it was concluded that both the trends from synthetic/MODIS NDVI dense time series and synthetic/MODIS NDVImax time series (2000–2011) were highly consistent. Compared to trends from MODIS time series, trends from synthetic time series are better able to capture fine scale vegetation changes. STARFM-generated synthetic NDVI time series could be used to quantify the effects of mining activities on vegetation, but the test areas should be selected with caution, as the trends derived from synthetic and MODIS time series may be significantly different in some areas.

Climate Contributions to Vegetation Variations in Central Asian Drylands: Pre- and Post-USSR Collapse

Central Asia comprises a large fraction of the worlds drylands, known to be vulnerable to climate change. We analyzed the inter-annual trends and the impact of climate variability in the vegetation greenness for Central Asia from 1982 to 2011 using GIMMS3g normalized difference vegetation index (NDVI) data. In our study, most areas showed an increasing trend during 1982-1991, but experienced a significantly decreasing trend for 1992-2011. Vegetation changes were closely coupled to climate variables (precipitation and temperature) during 1982-1991 and 1992-2011, but the response trajectories differed between these two periods. The warming trend in Central Asia initially enhanced the vegetation greenness before 1991, but the continued warming trend subsequently became a suppressant of further gains in greenness afterwards. Precipitation expanded its influence on larger vegetated areas in 1992-2011 when compared to 1982-1991. Moreover, the time-lag response of plants to rainfall tended to increase after 1992 compared to the pre-1992 period, indicating that plants might have experienced functional transformations to adapt the climate change during the study period. The impact of climate on vegetation was significantly different for the different sub-regions before and after 1992, coinciding with the collapse of the Union of Soviet Socialist Republics (USSR). It was suggested that these spatio-temporal patterns in greenness change and their relationship with climate change for some regions could be explained by the changes in the socio-economic structure resulted from the USSR collapse in late 1991. Our results clearly illustrate the combined influence of climatic/anthropogenic contributions on vegetation growth in Central Asian drylands. Due to the USSR collapse, this region represents a unique case study of the vegetation response to climate changes under different climatic and socio-economic conditions.

Comparisons of compositing period length for vegetation index data from polar-orbiting and geostationary satellites for the cloud-prone region of West Africa

Land surface data from MODIS and AVHRR have been extensively used for vegetation monitoring. In cloud-prone areas like West Africa the use of Normalized Difference Vegetation Index (NDVI) data for vegetation monitoring is hampered by persistent cloud cover especially during the rainy season. The new geostationary satellite Meteosat Second Generation (SEVIRI MSG) is the first geostationary satellite suited for vegetation monitoring allowing NDVI to be derived with a 15-minute temporal resolution. For West Africa, MODIS (combined TERRA and AQUA) produce above 85 percent cloud-free pixels in the scene during the entire rainy season using 16-day composite periods. SEVIRI MSG data produces � 98 percent cloud-free pixels during the entire season using a 3-day composite period. Therefore, there is a much higher probability for producing high quality cloud free data using SEVIRI MSG data for a short time composite period compared to Polar Orbiting Environmental Satellite (POES) data, which is expected to substantially improve various applications of satellite based natural resource management, including vegetation monitoring, in West Africa.

Mapping dynamics of deforestation and forest degradation in tropical forests using radar satellite data

Mapping anthropogenic forest disturbances has largely been focused on distinct delineations of events of deforestation using optical satellite images. In the tropics, frequent cloud cover and the challenge of quantifying forest degradation remain problematic. In this study, we detect processes of deforestation, forest degradation and successional dynamics, using long-wavelength radar (L-band from ALOS PALSAR) backscatter. We present a detection algorithm that allows for repeated disturbances on the same land, and identifies areas with slow- and fast-recovering changes in backscatter in close spatial and temporal proximity. In the study area in Madre de Dios, Peru, 2.3% of land was found to be disturbed over three years, with a false positive rate of 0.3% of area. A low, but significant, detection rate of degradation from sparse and small-scale selective logging was achieved. Disturbances were most common along the tri-national Interoceanic Highway, as well as in mining areas and areas under no land use allocation. A continuous spatial gradient of disturbance was observed, highlighting artefacts arising from imposing discrete boundaries on deforestation events. The magnitude of initial radar backscatter, and backscatter decrease, suggested that large-scale deforestation was likely in areas with initially low biomass, either naturally or since already under anthropogenic use. Further, backscatter increases following disturbance suggested that radar can be used to characterize successional disturbance dynamics, such as biomass accumulation in lands post-abandonment. The presented radar-based detection algorithm is spatially and temporally scalable, and can support monitoring degradation and deforestation in tropical rainforests with the use of products from ALOS-2 and the future SAOCOM and BIOMASS missions.