A global,high-resolution (30-m) inland water body dataset for 2000: first results of a topographic–spectral classification algorithm

The science and management of terrestrial ecosystems require accurate, high-resolution mapping of surface water. We produced a global, 30-m-resolution inland surface water dataset with an automated algorithm using Landsat-based surface reflectance estimates, multispectral water and vegetation indices, terrain metrics, and prior coarse-resolution water masks. The dataset identified 3,650,723 km² of inland water globally – nearly three quarters of which was located in North America (40.65%) and Asia (32.77%), followed by Europe (9.64%), Africa (8.47%), South America (6.91%), and Oceania (1.57%). Boreal forests contained the largest portion of terrestrial surface water (25.03% of the global total), followed by the nominal ‘inland water’ biome (16.36%), tundra (15.67%), and temperate broadleaf and mixed forests (13.91%). Agreement with respect to the Moderate-resolution Imaging Spectroradiometer water mask and Landsat-based national land-cover datasets was very high, with commission errors <4% and omission errors <14% relative to each. Most of these were accounted for in the seasonality of water cover, snow and ice, and clouds – effects which were compounded by differences in image acquisition date relative to reference datasets. The Global Land Cover Facility (GLCF) inland surface water dataset is available for open access at the GLCF website (http://www.landcover.org).     

Integrating global land cover products for improved forest cover characterization: an application in North America

Six widely used coarse-resolution global land cover data-sets – Global Land Cover Characterization (GLCC), Global Land Cover 2000 (GLC2000), GlobCover land cover product (GlobCover), MODIS land cover product (MODIS LC), the University of Maryland land cover product (UMD LC), and the MODIS Vegetation Continuous Fields tree cover layer (MODIS VCF) disagree substantially in their estimates of forest cover. Employing a regression tree model trained on higher-resolution, Landsat-based data, these multisource multiresolution maps were integrated for an improved characterization of forest cover over North America. Evaluated using a withheld test sample, the integrated percent forest cover (IPFC) data-set has a root mean square error of 11.75% – substantially better than the 17.37% of GLCC, 17.61% of GLC2000, 17.96% of GlobCover, 15.23% of MODIS LC, 19.25% of MODIS VCF, and 15.15% of UMD LC, respectively. Although demonstrated for forest, this approach based on integration of multiple products has potential for improved characterization of other land cover types as well.     

Evaluation of global land cover maps for cropland area estimation in the conterminous United States

Global land cover data could provide continuously updated cropland acreage and distribution information, which is essential to a wide range of applications over large geographical regions. Cropland area estimates were evaluated in the conterminous USA from four recent global land cover products: MODIS land cover (MODISLC) at 500-m resolution in 2010, GlobCover at 300-m resolution in 2009, FROM-GLC and FROM-GLC-agg at 30-m resolution based on Landsat imagery circa 2010 against the US Department of Agriculture survey data. Ratio estimators derived from the 30-m resolution Cropland Data Layer were applied to MODIS and GlobCover land cover products, which greatly improved the estimation accuracy of MODISLC by enhancing the correlation and decreasing mean deviation (MDev) and RMSE, but were less effective on GlobCover product. We found that, in the USA, the CDL adjusted MODISLC was more suitable for applications that concern about the aggregated county cropland acreage, while FROM-GLC-agg gave the least deviation from the survey at the state level. Correlation between land cover map estimates and survey estimates is significant, but stronger at the state level than at the county level. In regions where most mismatches happen at the county level, MODIS tends to underestimate, whereas MERIS and Landsat images incline to overestimate. Those uncertainties should be taken into consideration in relevant applications. Excluding interannual and seasonal effects, R² of the FROM-GLC regression model increased from 0.1 to 0.4, and the slope is much closer to one. Our analysis shows that images acquired in growing season are most suitable for Landsat-based cropland mapping in the conterminous USA.     

Production of global land cover data – GLCNMO

Abstract

Global land cover is one of the fundamental contents of Digital Earth. The Global Mapping project coordinated by the International Steering Committee for Global Mapping has produced a 1-km global land cover dataset – Global Land Cover by National Mapping Organizations. It has 20 land cover classes defined using the Land Cover Classification System. Of them, 14 classes were derived using supervised classification. The remaining six were classified independently: urban, tree open, mangrove, wetland, snow/ice, and water. Primary source data of this land cover mapping were eight periods of 16-day composite 7-band 1-km MODIS data of 2003. Training data for supervised classification were collected using Landsat images, MODIS NDVI seasonal change patterns, Google Earth, Virtual Earth, existing regional maps, and expert's comments. The overall accuracy is 76.5% and the overall accuracy with the weight of the mapped area coverage is 81.2%. The data are available from the Global Mapping project website (http://www.iscgm.org/). The MODIS data used, land cover training data, and a list of existing regional maps are also available from the CEReS website. This mapping attempt demonstrates that training/validation data accumulation from different mapping projects must be promoted to support future global land cover mapping.     

Updating Landsat-based forest cover maps with MODIS images using multiscale spectral-spatial-temporal superresolution mapping

With the high deforestation rates of global forest covers during the past decades, there is an ever-increasing need to monitor forest covers at both fine spatial and temporal resolutions. Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat series images have been used commonly for satellite-derived forest cover mapping. However, the spatial resolution of MODIS images and the temporal resolution of Landsat images are too coarse to observe forest cover at both fine spatial and temporal resolutions. In this paper, a novel multiscale spectral-spatial-temporal superresolution mapping (MSSTSRM) approach is proposed to update Landsat-based forest maps by integrating current MODIS images with the previous forest maps generated from Landsat image. Both the 240 m MODIS bands and 480 m MODIS bands were used as inputs of the spectral energy function of the MSSTSRM model. The principle of maximal spatial dependence was used as the spatial energy function to make the updated forest map spatially smooth. The temporal energy function was based on a multiscale spatial-temporal dependence model, and considers the land cover changes between the previous and current time. The novel MSSTSRM model was able to update Landsat-based forest maps more accurately, in terms of both visual and quantitative evaluation, than traditional pixel-based classification and the latest sub-pixel based super-resolution mapping methods The results demonstrate the great efficiency and potential of MSSTSRM for updating fine temporal resolution Landsat-based forest maps using MODIS images.     

Multiscale divergence between Landsat- and lidar-based biomass mapping is related to regional variation in canopy cover and composition

Background

Satellite-based aboveground forest biomass maps commonly form the basis of forest biomass and carbon stock mapping and monitoring, but biomass maps likely vary in performance by region and as a function of spatial scale of aggregation. Assessing such variability is not possible with spatially-sparse vegetation plot networks. In the current study, our objective was to determine whether high-resolution lidar-based and moderate-resolution Landsat-base aboveground live forest biomass maps converged on similar predictions at stand- to landscape-levels (10 s to 100 s ha) and whether such differences depended on biophysical setting. Specifically, we examined deviations between lidar- and Landsat-based biomass mapping methods across scales and ecoregions using a measure of error (normalized root mean square deviation), a measure of the unsystematic deviations, or noise (Pearson correlation coefficient), and two measures related to systematic deviations, or biases (intercept and slope of a regression between the two sets of predictions).

Results

Compared to forest inventory data (0.81-ha aggregate-level), lidar and Landsat-based mean biomass predictions exhibited similar performance, though lidar predictions exhibited less normalized root mean square deviation than Landsat when compared with the reference plot data. Across aggregate-levels, the intercepts and slopes of regression equations describing the relationships between lidar- and Landsat-based biomass predictions stabilized (i.e., little additional change with increasing area of aggregates) at aggregate-levels between 10 and 100 ha, suggesting a consistent relationship between the two maps at landscape-scales. Differences between lidar- and Landsat-based biomass maps varied as a function of forest canopy heterogeneity and composition, with systematic deviations (regression intercepts) increasing with mean canopy cover and hardwood proportion within forests and correlations decreasing with hardwood proportion.

Conclusions

Deviations between lidar- and Landsat-based maps indicated that satellite-based approaches may represent general gradients in forest biomass. Ecoregion impacted deviations between lidar and Landsat biomass maps, highlighting the importance of biophysical setting in determining biomass map performance across aggregate scales. Therefore, regardless of the source of remote sensing (e.g., Landsat vs. lidar), factors affecting the measurement and prediction of forest biomass, such as species composition, need to be taken into account whether one is estimating biomass at the plot, stand, or landscape scale.

     

Comparison of MODIS vegetation continuous field — based forest density maps with IRS-LISS III derived maps

The study compared forest cover maps derived using coarse resolution vegetation continuous fields (MODIS VCF; 500m resolution) with the maps derived from medium resolution (24m; IRS LISS-III) data. The comparison of VCF, per cent tree cover product, for the years 2000 to 2004 with LISS III forest density class maps of 2001 and 2003 was carried out for two sites representing hilly (Uttarakhand) and undulating terrains (Madhya Pradesh). Slicing VCF to corresponding forest crown cover, i.e., 0–10%, 10–40%, 40–70% and >70% produced considerable difference in forest area estimates when compared to original LISS III derived crown cover area. The corresponding value range in VCF for 0–10% of actual forest cover were 0–31% and 0–25% in 2 sites respectively, and the respective limit was consistent at 1–20% when VCF range were sliced with respect to upscaled LISS III at 500m resolution. Similarly, all other class limits were also found through iterative process. These limits were similar, within a site, across five years. Spatial Kappa match between these two data indicated higher match in 40–70% class, and also in undulating site. When compared at same resolution, similar forest area cover estimated with weighted area upscaling gave closest match. The study is useful in knowing the usability and limits of VCF product, and utility of spatial Kappa.     

A Spatio-Temporal Enhancement Method for medium resolution LAI (STEM-LAI)

Satellite remote sensing has been used successfully to map leaf area index (LAI) across landscapes, but advances are still needed to exploit multi-scale data streams for producing LAI at both high spatial and temporal resolution. A multi-scale Spatio-Temporal Enhancement Method for medium resolution LAI (STEM-LAI) has been developed to generate 4-day time-series of Landsat-scale LAI from existing medium resolution LAI products. STEM-LAI has been designed to meet the demands of applications requiring frequent and spatially explicit information, such as effectively resolving rapidly evolving vegetation dynamics at sub-field (30 m) scales. In this study, STEM-LAI is applied to Moderate Resolution Imaging Spectroradiometer (MODIS) based LAI data and utilizes a reference-based regression tree approach for producing MODIS-consistent, but Landsat-based, LAI. The Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) is used to interpolate the downscaled LAI between Landsat acquisition dates, providing a high spatial and temporal resolution improvement over existing LAI products. STARFM predicts high resolution LAI by blending MODIS and Landsat based information from a common acquisition date, with MODIS data from a prediction date. To demonstrate its capacity to reproduce fine-scale spatial features observed in actual Landsat LAI, the STEM-LAI approach is tested over an agricultural region in Nebraska. The implementation of a 250 m resolution LAI product, derived from MODIS 1 km data and using a scale consistent approach based on the Normalized Difference Vegetation Index (NDVI), is found to significantly improve accuracies of spatial pattern prediction, with the coefficient of efficiency (E) ranging from 0.77–0.94 compared to 0.01–0.85 when using 1 km LAI inputs alone. Comparisons against an 11-year record of in-situ measured LAI over maize and soybean highlight the utility of STEM-LAI in reproducing observed LAI dynamics (both characterized by r² = 0.86) over a range of plant development stages. Overall, STEM-LAI represents an effective downscaling and temporal enhancement mechanism that predicts in-situ measured LAI better than estimates derived through linear interpolation between Landsat acquisitions. This is particularly true when the in-situ measurement date is greater than 10 days from the nearest Landsat acquisition, with prediction errors reduced by up to 50%. With a streamlined and completely automated processing interface, STEM-LAI represents a flexible tool for LAI disaggregation in space and time that is adaptable to different land cover types, landscape heterogeneities, and cloud cover conditions.     

An approach for characterizing the distribution of shrubland ecosystem components as continuous fields as part of NLCD

Characterizing and quantifying distributions of shrubland ecosystem components is one of the major challenges for monitoring shrubland vegetation cover change across the United States. A new approach has been developed to quantify shrubland components as fractional products within National Land Cover Database (NLCD). This approach uses remote sensing data and regression tree models to estimate the fractional cover of shrubland ecosystem components. The approach consists of three major steps: field data collection, high resolution estimates of shrubland ecosystem components using WorldView-2 imagery, and coarse resolution estimates of these components across larger areas using Landsat imagery. This research seeks to explore this method to quantify shrubland ecosystem components as continuous fields in regions that contain wide-ranging shrubland ecosystems. Fractional cover of four shrubland ecosystem components, including bare ground, herbaceous, litter, and shrub, as well as shrub heights, were delineated in three ecological regions in Arizona, Florida, and Texas. Results show that estimates for most components have relatively small normalized root mean square errors and significant correlations with validation data in both Arizona and Texas. The distribution patterns of shrub height also show relatively high accuracies in these two areas. The fractional cover estimates of shrubland components, except for litter, are not well represented in the Florida site. The research results suggest that this method provides good potential to effectively characterize shrubland ecosystem conditions over perennial shrubland although it is less effective in transitional shrubland. The fractional cover of shrub components as continuous elements could offer valuable information to quantify biomass and help improve thematic land cover classification in arid and semiarid areas.     

Net primary productivity of forest stands in New Hampshire estimated from Landsat and MODIS satellite data

Background

A simulation model that relies on satellite observations of vegetation cover from the Landsat 7 sensor and from the Moderate Resolution Imaging Spectroradiometer (MODIS) was used to estimate net primary productivity (NPP) of forest stands at the Bartlett Experiment Forest (BEF) in the White Mountains of New Hampshire.

Results

Net primary production (NPP) predicted from the NASA-CASA model using 30-meter resolution Landsat inputs showed variations related to both vegetation cover type and elevational effects on mean air temperatures. Overall, the highest predicted NPP from the NASA-CASA model was for deciduous forest cover at low to mid-elevation locations over the landscape. Comparison of the model-predicted annual NPP to the plot-estimated values showed a significant correlation of R² = 0.5. Stepwise addition of 30-meter resolution elevation data values explained no more than 20% of the residual variation in measured NPP patterns at BEF. Both the Landsat 7 and the 250-meter resolution MODIS derived mean annual NPP predictions for the BEF plot locations were within ± 2.5% of the mean of plot estimates for annual NPP.

Conclusion

Although MODIS imagery cannot capture the spatial details of NPP across the network of closely spaced plot locations as well as Landsat, the MODIS satellite data as inputs to the NASA-CASA model does accurately predict the average annual productivity of a site like the BEF.     

Primitives as building blocks for constructing land cover maps

Land cover maps play an integral role in environmental management. However, countries and institutes encounter many challenges with producing timely, efficient, and temporally harmonized updates to their land cover maps. To address these issues we present a modular Regional Land Cover Monitoring System (RLCMS) architecture that is easily customized to create land cover products using primitive map layers. Primitive map layers are a suite of biophysical and end member maps, with land cover primitives representing the raw information needed to make decisions in a dichotomous key for land cover classification. We present best practices to create and assemble primitives from optical satellite using computing technologies, decision tree logic and Monte Carlo simulations to integrate their uncertainties. The concept is presented in the context of a regional land cover map based on a shared regional typology with 18 land cover classes agreed on by stakeholders from Cambodia, Laos PDR, Myanmar, Thailand, and Vietnam. We created annual map and uncertainty layers for the period 2000–2017. We found an overall accuracy of 94% when taking uncertainties into account. RLCMS produces consistent time series products using free long term historical Landsat and MODIS data. The customizable architecture can include a variety of sensors and machine learning algorithms to create primitives and the best suited smoothing can be applied on a primitive level. The system is transferable to all regions around the globe because of its use of publicly available global data (Landsat and MODIS) and easily adaptable architecture that allows for the incorporation of a customizable assembly logic to map different land cover typologies based on the user's landscape monitoring objectives     

A Landsat surface reflectance dataset for North America, 1990-2000

The Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center has processed and released 2100 Landsat Thematic Mapper and Enhanced Thematic Mapper Plus surface reflectance scenes, providing 30-m resolution wall-to-wall reflectance coverage for North America for epochs centered on 1990 and 2000. This dataset can support decadal assessments of environmental and land-cover change, production of reflectance-based biophysical products, and applications that merge reflectance data from multiple sensors [e.g., the Advanced Spaceborne Thermal Emission and Reflection Radiometer, Multiangle Imaging Spectroradiometer, Moderate Resolution Imaging Spectroradiometer (MODIS)]. The raw imagery was obtained from the orthorectified Landsat GeoCover dataset, purchased by NASA from the Earth Satellite Corporation. Through the LEDAPS project, these data were calibrated, converted to top-of-atmosphere reflectance, and then atmospherically corrected using the MODIS/6S methodology. Initial comparisons with ground-based optical thickness measurements and simultaneously acquired MODIS imagery indicate comparable uncertainty in Landsat surface reflectance compared to the standard MODIS reflectance product (the greater of 0.5% absolute reflectance or 5% of the recorded reflectance value). The rapid automated nature of the processing stream also paves the way for routine high-level products from future Landsat sensors.     

Landsat 8 vs. Landsat 5: A comparison based on urban and peri-urban land cover mapping

An image dataset from the Landsat OLI spaceborne sensor is compared with the Landsat TM in order to evaluate the excellence of the new imagery in urban landcover classification. Widely known pixel-based and object-based image analysis methods have been implemented in this work like Maximum Likelihood, Support Vector Machine, k-Nearest Neighbor, Feature Analyst and Sub-pixel. Classification results from Landsat OLI provide more accurate results comparing to the Landsat TM. Object-based classifications produced a more uniform result, but suffer from the absorption of small rare classes into large homogenous areas, as a consequence of the segmentation, merging and the spatial parameters in the spatial resolution (30 m) of Landsat images. Based exclusively on the overall accuracy reports, the SVM pixel-based classification from Landsat 8 proved to be the most accurate for the purpose of mapping urban land cover, using medium spatial resolution imagery.     

Estimating fractional green vegetation cover of Mongolian grasslands using digital camera images and MODIS satellite vegetation indices

ABSTRACT

Fractional green vegetation cover (FVC) is a useful indicator for monitoring grassland status. Satellite imagery with coarse spatial but high temporal resolutions has been preferred to monitor seasonal and inter-annual FVC dynamics in wide geographic area such as Mongolian steppe. However, the coarse spatial resolution can cause a certain uncertainty in the satellite-based FVC estimation, which calls attention to develop a robust statistical test for the relationship between field FVC and satellite-derived vegetation indices. In the arid and semi-arid Mongolian steppe, nadir pointing digital camera images (DCI) were collected and used to produce a FVC dataset to support the evaluation of satellite-based FVC retrievals. An optimal DCI processing method was determined with respect to three color spaces (RGB, HIS, L*a*b*) and six green pixel classification algorithms, from which a country-wide dataset of DCI-FVC was produced and used for evaluating the accuracy of satellite-based FVC estimates from MODIS vegetation indices. We applied three empirical and three semi-empirical MODIS-FVC retrieval models. DCI data were collected from 96 sites across the Mongolian steppe from 2012 to 2014. The histogram algorithm using the hue (H) value of the HIS color space was the optimal DCI method (r² = 0.94, percent root-mean-square-error (RMSE) = 7.1%). For MODIS-FVC retrievals, semi-empirical Baret model was the best-performing model with the highest r² (0.69) and the lowest RMSE (49.7%), while the lowest MB (+1.1%) was found for the regression model with normalized difference vegetation index (NDVI). The high RMSE (>50% or so) is an issue requiring further enhancement of satellite-based FVC retrievals accounting for key plant and soil parameters relevant to the Mongolian steppe and for scale mismatch between sampling and MODIS data.     

Inter-comparison of remote sensing platforms for height estimation of mango and avocado tree crowns

To support the adoption of precision agricultural practices in horticultural tree crops, prior research has investigated the relationship between crop vigour (height, canopy density, health) as measured by remote sensing technologies, to fruit quality, yield and pruning requirements. However, few studies have compared the accuracy of different remote sensing technologies for the estimation of tree height. In this study, we evaluated the accuracy, flexibility, aerial coverage and limitations of five techniques to measure the height of two types of horticultural tree crops, mango and avocado trees. Canopy height estimates from Terrestrial Laser Scanning (TLS) were used as a reference dataset against height estimates from Airborne Laser Scanning (ALS) data, WorldView-3 (WV-3) stereo imagery, Unmanned Aerial Vehicle (UAV) based RGB and multi-spectral imagery, and field measurements. Overall, imagery obtained from the UAV platform were found to provide tree height measurement comparable to that from the TLS (R² = 0.89, RMSE = 0.19 m and rRMSE = 5.37 % for mango trees; R² = 0.81, RMSE = 0.42 m and rRMSE = 4.75 % for avocado trees), although coverage area is limited to 1–10 km² due to battery life and line-of-sight flight regulations. The ALS data also achieved reasonable accuracy for both mango and avocado trees (R² = 0.67, RMSE = 0.24 m and rRMSE = 7.39 % for mango trees; R² = 0.63, RMSE = 0.43 m and rRMSE = 5.04 % for avocado trees), providing both optimal point density and flight altitude, and therefore offers an effective platform for large areas (10 km²–100 km²). However, cost and availability of ALS data is a consideration. WV-3 stereo imagery produced the lowest accuracies for both tree crops (R² = 0.50, RMSE = 0.84 m and rRMSE = 32.64 % for mango trees; R² = 0.45, RMSE = 0.74 m and rRMSE = 8.51 % for avocado trees) when compared to other remote sensing platforms, but may still present a viable option due to cost and commercial availability when large area coverage is required. This research provides industries and growers with valuable information on how to select the most appropriate approach and the optimal parameters for each remote sensing platform to assess canopy height for mango and avocado trees.     

Forest cover classification using Landsat ETM+ data and time series MODIS NDVI data

Forest cover plays a key role in climate change by influencing the carbon stocks, the hydrological cycle and the energy balance. Forest cover information can be determined from fine-resolution data, such as Landsat Enhanced Thematic Mapper Plus (ETM+). However, forest cover classification with fine-resolution data usually uses only one temporal data because successive data acquirement is difficult. It may achieve mis-classification result without involving vegetation growth information, because different vegetation types may have the similar spectral features in the fine-resolution data. To overcome these issues, a forest cover classification method using Landsat ETM+ data appending with time series Moderate-resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) data was proposed. The objective was to investigate the potential of temporal features extracted from coarse-resolution time series vegetation index data on improving the forest cover classification accuracy using fine-resolution remote sensing data. This method firstly fused Landsat ETM+ NDVI and MODIS NDVI data to obtain time series fine-resolution NDVI data, and then the temporal features were extracted from the fused NDVI data. Finally, temporal features combined with Landsat ETM+ spectral data was used to improve forest cover classification accuracy using supervised classifier. The study in North China region confirmed that time series NDVI features had significant effects on improving forest cover classification accuracy of fine resolution remote sensing data. The NDVI features extracted from time series fused NDVI data could improve the overall classification accuracy approximately 5% from 88.99% to 93.88% compared to only using single Landsat ETM+ data.     

Estimating spatial variations of total evaporation using multispectral sensors within the uMngeni catchment,South Africa

Total evaporation is of importance in assessing and managing long-term water use, especially in water-limited environments. Therefore, there is need to account for water utilisation by different land uses for well-informed water resources management and future planning. This study investigated the feasibility of using multispectral Landsat 8 and moderate resolution imaging spectroradiometer (MODIS) remote sensing data to estimate total evaporation within the uMngeni catchment in South Africa, using surface energy balance system. The results indicated that Landsat 8 at 30 m resolution has a better spatial representation of total evaporation, when compared to the 1000 m MODIS. Specifically, Landsat 8 yielded significantly different mean total evaporation estimates for all land cover types (one-way ANOVA; F_4.964?=?87.011, p < 0.05), whereas MODIS failed to differentiate (one-way ANOVA; F_2.853?=?0.125, p = 0.998) mean total evaporation estimates for the different land cover types across the catchment. The findings of this study underscore the utility of the Landsat 8 spatial resolution and land cover characteristics in deriving accurate and reliable spatial variations of total evaporation at a catchment scale.     

Comparison of the Landsat Surface Reflectance Climate Data Record (CDR) and manually atmospherically corrected data in a semi-arid European study area

This study contributes to the quality assessment of atmospherically corrected Landsat surface reflectance data that are routinely generated by the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS). This dataset, named Landsat Surface Reflectance Climate Data Record (Landsat CDR), is available at global scale and offers unprecedented opportunities to land monitoring and management services that require atmospherically corrected Earth observation (EO) data. Our assessment is based on the comparison of the Landsat CDR data against a set of Landsat and DEIMOS-1 images processed to a high degree of accuracy using an industry-standard atmospheric correction algorithm (ATCOR-2). The software package has been used for many years and its correction procedures can be considered consolidated and well-established. The dataset of Landsat and DEIMOS-1 images was acquired over a semi-arid agricultural area located in Lower Austria and was independently corrected by using a manual fine-tuning of ATCOR-2 parameters to reach the highest possible accuracy. Results show a very good correspondence of the surface reflectance in each of the six reflective spectral channels as well as for the NDVI (Normalized Difference Vegetation Index). An additional comparison against a NDVI time series from MODIS revealed also a good correspondence. Coefficients of determination (R²) between the two multi-year and multi-seasonal Landsat/DEIMOS datasets range between 0.91 (blue band) and 0.98 (nIR, SWIR-1 and SWIR-2). The results obtained for our semi-arid test site in Austria confirm previous findings and suggest that automatic atmospheric procedures, such as the one implemented by LEDAPS are accurate enough to be used in land monitoring services that require consistent multi-temporal surface reflectance data.     

Spatial analysis of human-induced vegetation productivity decline over eastern Africa using a decade (2001–2011) of medium resolution MODIS time-series data

Climate variation and land transformations related to exploitative land uses are among the main drivers of vegetation productivity decline and ongoing land degradation in East Africa. We combined analysis of vegetation trends and cumulative rain use efficiency differences (CRD), calculated from 250-m MODIS NDVI time-series data, to map vegetation productivity loss over eastern Africa between 2001 and 2011. The CRD index values were furthermore used to discern areas of particular severe vegetation productivity loss over the observation period. Monthly 25-km Tropical Rainfall Measuring Mission (TRMM) data metrics were used to mask areas of rainfall declines not related to human-induced land productivity loss. To provide insights on the productivity decline, we linked the MODIS-based vegetation productivity map to land transformation processes using very high resolution (VHR) imagery in Google Earth (GE) and a Landsat-based land-cover change map. In total, 3.8 million ha experienced significant vegetation loss over the monitoring period. An overall agreement of 68% was found between the rainfall-corrected MODIS productivity decline map and all reference pixels discernable from GE and the Landsat map. The CRD index showed a good potential to discern areas with ‘severe’ vegetation productivity losses under high land-use intensities.