Atmospheric correction of optical imagery from MODIS and Reanalysis atmospheric products

In this paper we analyze the differences obtained in the atmospheric correction of optical imagery covering bands located in the Visible and Near Infra-Red (VNIR), Short-Wave Infra-Red (SWIR) and Themal-Infrared (TIR) spectral regions when atmospheric profiles extracted from different sources are used. In particular, three sensors were used, Compact High Resolution Imaging Spectrometer (CHRIS), Advanced Spaceborne Thermal Emission and Reflection radiometer (ASTER) and Landsat5 Thematic Mapper (TM), whereas four atmospheric profiles sources were considered: i) local soundings launched near the sensor overpass time, ii) Moderate Resolution Radiometer (MODIS) atmospheric profiles product (MOD07), iii) Atmospheric Correction Parameter Calculator (ACPC) generated by the National Center for Environmental Prediction (NCEP) and iv) Modified Atmospheric Profiles from Reanalysis Information (MAPRI), which includes data from NCEP and National Center of Atmospheric Research (NCAR) Reanalysis project but interpolated to 34 atmospheric levels and resampled to 0.5° × 0.5°. MODIS aerosol product (MOD04) was also used to extract Aerosol Optical Thickness (AOT) values at 550 nm. Analysis was performed for three test dates (12th July 2003, 18th July 2004 and 13th July 2005) over an agricultural area in Spain. Results showed that air temperature vertical profiles were similar for the four sources, whereas dew point temperature profiles showed significant differences at some particular levels. Atmospheric profiles were used as input to MODTRAN4 radiative transfer code in order to compute atmospheric parameters involved in atmospheric correction, with the aim of retrieving surface reflectances in the case of VNIR and SWIR regions, and Land Surface Temperature (LST) in the case of the TIR region. For the VNIR and SWIR region, significant differences depending on the atmospheric profile used were not found, particularly in the Visible region in which the AOT content is the main parameter involved in the atmospheric correction. In the case of TIR, differences depending on the atmospheric profile used were appreciable, since in this case the main parameter involved in the atmospheric correction is the water vapor content, which depends on the vertical profile. In terms of LST retrieval from ASTER data (2004 test case), all profiles provided satisfactory results compared to the ones obtained when using a local sounding, with errors of 0.3 K for ACPC and MAPRI cases and 0.7 K for MOD07. When retrieving LST from TM data (2005 test case), errors for MOD07 and MAPRI were 0.6 and 0.9 K respectively, whereas ACPC provided an error of 2 K. The results presented in this paper show that the different atmospheric profile sources are useful for accurate atmospheric correction when local soundings are not available. In particular, MOD07 product provides atmospheric information at the highest spatial resolution, 5 km, although its use is limited from 2000 to present, whereas MAPRI provides historical information from 1970 to present, but at lower spatial resolution.     

A scalable approach to mapping annual land cover at 250 m using MODIS time series data: A case study in the Dry Chaco ecoregion of South America

Land use and land cover (LULC) maps from remote sensing are vital for monitoring, understanding and predicting the effects of complex human-nature interactions that span local, regional and global scales. We present a method to map annual LULC at a regional spatial scale with source data and processing techniques that permit scaling to broader spatial and temporal scales, while maintaining a consistent classification scheme and accuracy. Using the Dry Chaco ecoregion in Argentina, Bolivia and Paraguay as a test site, we derived a suite of predictor variables from 2001 to 2007 from the MODIS 250 m vegetation index product (MOD13Q1). These variables included: annual statistics of red, near infrared, and enhanced vegetation index (EVI), phenological metrics derived from EVI time series data, and slope and elevation. For reference data, we visually interpreted percent cover of eight classes at locations with high-resolution QuickBird imagery in Google Earth. An adjustable majority cover threshold was used to assign samples to a dominant class. When compared to field data, we found this imagery to have georeferencing error < 5% the length of a MODIS pixel, while most class interpretation error was related to confusion between agriculture and herbaceous vegetation. We used the Random Forests classifier to identify the best sets of predictor variables and percent cover thresholds for discriminating our LULC classes. The best variable set included all predictor variables and a cover threshold of 80%. This optimal Random Forests was used to map LULC for each year between 2001 and 2007, followed by a per-pixel, 3-year temporal filter to remove disallowed LULC transitions. Our sequence of maps had an overall accuracy of 79.3%, producer accuracy from 51.4% (plantation) to 95.8% (woody vegetation), and user accuracy from 58.9% (herbaceous vegetation) to 100.0% (water). We attributed map class confusion to limited spectral information, sub-pixel spectral mixing, georeferencing error and human error in interpreting reference samples. We used our maps to assess woody vegetation change in the Dry Chaco from 2002 to 2006, which was characterized by rapid deforestation related to soybean and planted pasture expansion. This method can be easily applied to other regions or continents to produce spatially and temporally consistent information on annual LULC.     

Testing a MODIS Global Disturbance Index across North America

Large-scale ecosystem disturbances (LSEDs) have major impacts on the global carbon cycle as large pulses of CO₂ and other trace gases from terrestrial biomass loss are emitted to the atmosphere during disturbance events. The high temporal and spatial variability of the atmospheric emissions combined with the lack of a proven methodology to monitor LSEDs at the global scale make the timing, location and extent of vegetation disturbance a significant uncertainty in understanding the global carbon cycle. The MODIS Global Disturbance Index (MGDI) algorithm is designed for large-scale, regular, disturbance mapping using Aqua/Moderate Resolution Imaging Spectroradiometer (MODIS) Land Surface Temperature (LST) and Aqua/MODIS Enhanced Vegetation Index (EVI) data. The MGDI uses annual maximum composite LST data to detect fundamental changes in land-surface energy partitioning, while avoiding the high natural variability associated with tracking LST at daily, weekly, or seasonal time frames. Here we apply the full Aqua/MODIS dataset through 2006 to the improved MGDI algorithm across the woody ecosystems of North America and test the algorithm by comparison with confirmed, historical wildfire events and the windfall areas of documented major hurricanes. The MGDI accurately detects the location and extent of wildfire throughout North America and detects high and moderate severity impacts in the windfall area of major hurricanes. We also find detections associated with clear-cut logging and land-clearing on the forest-agricultural interface. The MGDI indicates that 1.5% (195,580 km²) of the woody ecosystems within North America was disturbed in 2005 and 0.5% (67,451 km²) was disturbed in 2006. The interannual variability is supported by wildfire detections and official burned area statistics.     

Evaluation of ASTER and MODIS land surface temperature and emissivity products using long-term surface longwave radiation observations at SURFRAD sites

Land surface temperature (LST) and emissivity are key parameters in estimating the land surface radiation budget, a major controlling factor of global climate and environmental change. In this study, Terra Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) and Aqua MODerate resolution Imaging Spectroradiometer (MODIS) Collection 5 LST and emissivity products are evaluated using long-term ground-based longwave radiation observations collected at six Surface Radiation Budget Network (SURFRAD) sites from 2000 to 2007. LSTs at a spatial resolution of 90 m from 197 ASTER images during 2000-2007 are directly compared to ground observations at the six SURFRAD sites. For nighttime data, ASTER LST has an average bias of 0.1 °C and the average bias is 0.3 °C during daytime. Aqua MODIS LST at 1 km resolution during nighttime retrieved from a split-window algorithm is evaluated from 2002 to 2007. MODIS LST has an average bias of − 0.2 °C. LST heterogeneity (defined as the Standard Deviation, STD, of ASTER LSTs in 1 × 1 km² region, 11 × 11 pixel in total) and instrument calibration error of pyrgeometer are key factors impacting the ASTER and MODIS LST evaluation using ground-based radiation measurements. The heterogeneity of nighttime ASTER LST is 1.2 °C, which accounts for 71% of the STD of the comparison, while the heterogeneity of the daytime LST is 2.4 °C, which accounts for 60% of the STD. Collection 5 broadband emissivity is 0.01 larger than that of MODIS Collection 4 products and ASTER emissivity. It is essential to filter out the abnormal low values of ASTER daily emissivity data in summer time before its application.     

Air temperature estimation with MSG-SEVIRI data: Calibration and validation of the TVX algorithm for the Iberian Peninsula

Air temperature can be estimated from remote sensing by combining information in thermal infrared and optical wavelengths. The empirical TVX algorithm is based on an estimated linear relationship between observed Land Surface Temperature (LST) and a Spectral Vegetation Index (NDVI). Air temperature is assumed to be equal to the LST corresponding to the effective full vegetation cover, and is found by extrapolating the line to a maximum value of NDVI_max. The algorithm has been tested and reported in the literature previously. However, the effect of vegetation types and climates and the potential variation in NDVI of the effective full cover has not been subject for investigation. The present study proposes a novel methodology to estimate NDVI_max that uses observed air temperature to calibrate the NDVI_max for each vegetation type. To assess the validity of this methodology, we have compared the accuracy of estimates using the new NDVI_max and the previous NDVI_max that have been proposed in literature with MSG-SEVIRI images in Spain during the year 2005. In addition, a spatio-temporal assessment of residuals has been performed to evaluate the accuracy of retrievals in terms of daily and seasonal variation, land cover, landscape heterogeneity and topography. Results showed that the new calibrated NDVI_max perform well, with a Mean Absolute Error ranging between 2.8 °C and 4 °C. In addition, vegetation-specific NDVI_max improve the accuracy compared with a unique NDVI_max.     

Land cover classification with coarse spatial resolution data to derive continuous and discrete maps for complex regions

Regularly updated land cover information at continental or national scales is a requirement for various land management applications as well as biogeochemical and climate modeling exercises. However, monitoring or updating of map products with sufficient spatial detail is currently not widely practiced due to inadequate time-series coverage for most regions of the Earth. Classifications of coarser spatial resolution data can be automatically generated on an annual or finer time scale. However, discrete land cover classifications of such data cannot sufficiently quantify land surface heterogeneity or change. This study presents a methodology for continuous and discrete land cover mapping using moderate spatial resolution time series data sets. The method automatically selects sample data from higher spatial resolution maps and generates multiple decision trees. The leaves of decision trees are interpreted considering the sample distribution of all classes yielding class membership maps, which can be used as estimates for the diversity of classes in a coarse resolution cell. Results are demonstrated for the heterogeneous, small-patch landscape of Germany and the bio-climatically varying landscape of South Africa. Results have overall classification accuracies of 80%. A sensitivity analysis of individual modules of the classification process indicates the importance of appropriately chosen features, sample data balanced among classes, and an appropriate method to combine individual classifications. The comparison of classification results over several years not only indicates the method's consistency, but also its potential to detect land cover changes.     

Object-based analysis and change detection of major wetland cover types and their classification uncertainty during the low water period at Poyang Lake, China

Productive wetland systems at land-water interfaces that provide unique ecosystem services are challenging to study because of water dynamics, complex surface cover and constrained field access. We applied object-based image analysis and supervised classification to four 32-m Beijing-1 microsatellite images to examine broad-scale surface cover composition and its change during November 2007-March 2008 low water season at Poyang Lake, the largest freshwater lake-wetland system in China (> 4000 km²). We proposed a novel method for semi-automated selection of training objects in this heterogeneous landscape using extreme values of spectral indices (SIs) estimated from satellite data. Dynamics of the major wetland cover types (Water, Mudflat, Vegetation and Sand) were investigated both as transitions among primary classes based on maximum membership value, and as changes in memberships to all classes even under no change in a primary class. Fuzzy classification accuracy was evaluated as match frequencies between classification outcome and a) the best reference candidate class (MAX function) and b) any acceptable reference class (RIGHT function). MAX-based accuracy was relatively high for Vegetation (≥ 90%), Water (≥ 82%), Mudflat (≥ 76%) and the smallest-area Sand (≥ 75%) in all scenes; these scores improved with the RIGHT function to 87-100%. Classification uncertainty assessed as the proportion of fuzzy object area within a class at a given fuzzy threshold value was the highest for all classes in November 2007, and consistently higher for Mudflat than for other classes in all scenes. Vegetation was the dominant class in all scenes, occupying 41.2-49.3% of the study area. Object memberships to Vegetation mostly declined from November 2007 to February 2008 and increased substantially only in February-March 2008, possibly reflecting growing season conditions and grazing. Spatial extent of Water both declined and increased during the study period, reflecting precipitation and hydrological events. The “fuzziest” Mudflat class was involved in major detected transitions among classes and declined in classification accuracy by March 2008, representing a key target for finer-scale research. Future work should introduce Vegetation sub-classes reflecting differences in phenology and alternative methods to discriminate Mudflat from other classes. Results can be used to guide field sampling and top-down landscape analyses in this wetland.     

Reprocessing the MODIS Leaf Area Index products for land surface and climate modelling

Land surface and climate modelling requires continuous and consistent Leaf Area Index (LAI). High spatiotemporal resolution and long-time record data are more in demand nowadays and will continue to be in the future. MODIS LAI products meet these requirements to some degree. However, due to the presence of cloud and seasonal snow cover, the instrument problems and the uncertainties of retrieval algorithm, the current MODIS LAI products are spatially and temporally discontinuous and inconsistent, which limits their application in land surface and climate modelling. To improve the MODIS LAI products on a global scale, we considered the characteristics of the MODIS LAI data and made the best use of quality control (QC) information, and developed an integrated two-step method to derive the improved MODIS LAI products effectively and efficiently on a global scale. First, we used the modified temporal spatial filter (mTSF) method taking advantage of background values and QC information at each pixel to do a simple data assimilation for relatively low quality data. Then we applied the post processing-TIMESAT (A software package to analyze time-series of satellite sensor data) Savitzky-Golay (SG) filter to get the final result. We implemented the method to 10 years of the MODIS Collection 5 LAI data. In comparison with the LAI reference maps and the MODIS LAI data, our results showed that the improved MODIS LAI data are closer to the LAI reference maps in magnitude and also more continuous and consistent in both time-series and spatial domains. In addition, simple statistics were used to evaluate the differences between the MODIS LAI and the improved MODIS LAI.     

Thematic accuracy of the National Land Cover Database (NLCD) 2001 land cover for Alaska

The National Land Cover Database (NLCD) 2001 Alaska land cover classification is the first 30-m resolution land cover product available covering the entire state of Alaska. The accuracy assessment of the NLCD 2001 Alaska land cover classification employed a geographically stratified three-stage sampling design to select the reference sample of pixels. Reference land cover class labels were determined via fixed wing aircraft, as the high resolution imagery used for determining the reference land cover classification in the conterminous U.S. was not available for most of Alaska. Overall thematic accuracy for the Alaska NLCD was 76.2% (s.e. 2.8%) at Level II (12 classes evaluated) and 83.9% (s.e. 2.1%) at Level I (6 classes evaluated) when agreement was defined as a match between the map class and either the primary or alternate reference class label. When agreement was defined as a match between the map class and primary reference label only, overall accuracy was 59.4% at Level II and 69.3% at Level I. The majority of classification errors occurred at Level I of the classification hierarchy (i.e., misclassifications were generally to a different Level I class, not to a Level II class within the same Level I class). Classification accuracy was higher for more abundant land cover classes and for pixels located in the interior of homogeneous land cover patches.     

Urban growth modeling of Kathmandu metropolitan region, Nepal

The complexity of urban system requires integrated tools and techniques to understand the spatial process of urban development and project the future scenarios. This research aims to simulate urban growth patterns in Kathmandu metropolitan region in Nepal. The region, surrounded by complex mountainous terrain, has very limited land resources for new developments. As similar to many cities of the developing world, it has been facing rapid population growth and daunting environmental problems. Three time series land use maps in a fine-scale (30 m resolution), derived from satellite remote sensing, for the last three decades of the 20th century were used to clarify the spatial process of urbanization. Based on the historical experiences of the land use transitions, we adopted weight of evidence method integrated in cellular automata framework for predicting the future spatial patterns of urban growth. We extrapolated urban development patterns to 2010 and 2020 under the current scenario across the metropolitan region. Depending on local characteristics and land cover transition rates, this model produced noticeable spatial pattern of changes in the region. Based on the extrapolated spatial patterns, the urban development in the Kathmandu valley will continue through both in-filling in existing urban areas and outward rapid expansion toward the east and south directions. Overall development will be greatly affected by the existing urban space, transportation network, and topographic complexity.