Mapping burned areas from Landsat TM/ETM+ data with a two-phase algorithm: Balancing omission and commission errors

Maps of burned area have been obtained from an automatic algorithm applied to a multitemporal series of Landsat TM/ETM+ images in two Mediterranean sites. The proposed algorithm is based on two phases: the first one intends to detect the more severely burned areas and minimize commission errors. The second phase improves burned patches delimitation using a hybrid contextual algorithm based on logistic regression analysis, and tries to minimize omission errors. The algorithm was calibrated using six study sites and it was validated for the whole territory of Portugal (89,000 km²) and for Southern California (70,000 km²). In the validation exercise, 65 TM/ETM+ scenes for Portugal and 35 for California were used, all from the 2003 fire season. A good agreement with the official burned area perimeters was shown, with kappa values close to 0.85 and low omission and commission errors (< 16.5%). The proposed algorithm could be operationally used for historical mapping of burned areas from Landsat images, as well as from future medium resolution sensors, providing they acquire images in two bands of the Short Wave Infrared (1.5-2.2 μm).     

Comparing the spectral signal of burned surfaces between Landsat 7 ETM+ and Landsat 8 OLI sensors

There is a long history of the use of Landsat data in burned land mapping mainly due to certain characteristics of the Landsat imagery including the spatial, spectral, and temporal data resolution, the low cost (Landsat data are now freely available), and the existence of an almost 35-year historical archive (excluding Landsat 1–3). Landsat 8 (Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS)) was launched on 11 February 2013 and it captures data in three new bands along with two additional thermal bands. However, is the spectral signal of burned surfaces in satellite remote-sensing data of Landsat series consistent and robust enough to allow the successful application of the techniques developed so far for Landsat 8? In this article, we compare the spectral signal of burned surfaces between Landsat 7 Enhanced Thematic Mapper Plus (ETM+) and Landsat 8 OLI sensors using five case studies that correspond to five large fire events in different biophysical environments in Greece, for which both Landsat 7 ETM+ and Landsat 8 OLI data were available. From the comparative analysis using histogram data plots of burned (post-fire image) and vegetated (pre-fire image) areas, spectral signature plots and separability indices of certain land-cover types, estimated using the same sampling areas over both satellite images, a general consistency was observed between the two sensors. Slight differences between the sensors were attributed to differences in the acquisition dates and were related to the type of vegetation rather than the sensors used to record the satellite images. Neither sensor provided improved discrimination over the other.     

Application of geographically weighted regression to fill gaps in SLC-off Landsat ETM+ satellite imagery

Landsat 7 enhanced thematic mapper plus (ETM+) satellite imagery is an important data source for many applications. However, the scan line corrector (SLC) failed on 31 May 2003. As a result of the SLC failure, about 22% of the image data is missing in each scene; this is especially pronounced away from nadir. In this article, a local regression method called geographically weighted regression (GWR) is introduced for filling the gaps of the Landsat ETM+ imagery, and it is compared with kriging/cokriging for this purpose. The case studies show that the GWR approach is an effective technique to fill gaps in Landsat ETM+ imagery, although the image restoration is still not perfect. GWR performed marginally better than the complex cokriging method, which too has proven to be an effective method, but is computationally intensive. Although there are visible seam lines at the edges of the filled wide gaps in some bands, the validation results – including RMSE values, error distribution maps, and classification results for the case studies – demonstrate that the DN values estimated by GWR are in fact closer to those of the original image than the corresponding values estimated by kriging/cokriging.     

A comparison of empirical and neural network approaches for estimating corn and soybean leaf area index from Landsat ETM+ imagery

Plant foliage density expressed as leaf area index (LAI) is used in many ecological, meteorological, and agronomic models, and as a means of quantifying crop spatial variability for precision farming. LAI retrieval using spectral vegetation indices (SVI) from optical remotely sensed data usually requires site-specific calibration values from the surface or the use of within-scene image information without surface calibrations to invert radiative transfer models. An evaluation of LAI retrieval methods was conducted using (1) empirical methods employing the normalized difference vegetation index (NDVI) and a new SVI that uses green wavelength reflectance, (2) a scaled NDVI approach that uses no calibration measurements, and (3) a hybrid approach that uses a neural network (NN) and a radiative transfer model without site-specific calibration measurements. While research has shown that under a variety of conditions NDVI is not optimal for LAI retrieval, its continued use for remote sensing applications and in analysis seeking to develop improved parameter retrieval algorithms based on NDVI suggests its value as a “benchmark” or standard against which other methods can be compared. Landsat-7 ETM+ data for July 1 and July 8 from the Soil Moisture EXperiment 2002 (SMEX02) field campaign in the Walnut Creek watershed south of Ames, IA, were used for the analysis. Sun photometer data collected from a site within the watershed were used to atmospherically correct the imagery to surface reflectance. LAI validation measurements of corn and soybeans were collected close to the dates of the Landsat-7 overpasses. Comparable results were obtained with the empirical SVI methods and the scaled SVI method within each date. The hybrid method, although promising, did not account for as much of the variability as the SVI methods. Higher atmospheric optical depths for July 8 leading to surface reflectance errors are believed to have resulted in overall poorer performance for this date. Use of SVIs employing green wavelengths, improved method for the definition of image minimum and maximum clusters used by the scaled NDVI method, and further development of a soil reflectance index used by the hybrid NN approach are warranted. More importantly, the results demonstrate that reasonable LAI estimates are possible using optical remote sensing methods without in situ, site-specific calibration measurements.     

Performance of atmospheric and topographic correction methods on Landsat imagery in mountain areas

An effective removal of atmospheric and topographic effects on remote-sensing imagery is an essential preprocessing step for mapping land cover accurately in mountain areas. Various techniques that remove these effects have been proposed and consist of specific combinations of an atmospheric and a topographic correction (TC) method. However, it is possible to generate a wide range of new combined correction methods by applying alternative combinations of atmospheric and TC methods. At present, a systematic overview of the statistical performance and data input requirement of preprocessing techniques is missing. In order to assess the individual and combined impacts of atmospheric and TC methods, 15 permutations of two atmospheric and/or four TC methods were evaluated statistically and compared to the uncorrected imagery. Furthermore, results of the integrated ATCOR3 method were included. This evaluation was performed in a study area in the Romanian Carpathian mountains. Results showed that the combination of a transmittance-based atmospheric correction (AC), which corrects the effects of Rayleigh scattering and water-vapour absorption, and a pixel-based C or Minnaert TC, which account for diffuse sky irradiance, reduced the image distortions most efficiently. Overall results indicated that TC had a larger impact than AC and there was a trade-off between the statistical performance of preprocessing techniques and their data requirement. However, the normalized difference vegetation index analysis indicated that atmospheric methods resulted in a larger impact on the spectral information in bands 3 and 4.     

Assessment of fire severity and species diversity in the southern Appalachians using Landsat TM and ETM+ imagery

Relatively little is known about the disturbance ecology of large wildfires in the southern Appalachians. The occurrence of a 4000-ha wildfire in the Linville Gorge Wilderness area in western North Carolina has provided a rare opportunity to study a large fire with a range of severities. The objectives of this study were to 1) assess the potential for using multi-temporal Landsat imagery to map fire severity in the southern Appalachians, 2) examine the influences of topography and forest community type on the spatial pattern of fire severity; and 3) examine the relationship between predicted fire severity and changes in species richness. A non-linear regression equation predicted a field-based composite burn index (CBI) as a function of change in the Normalized Burn Ratio (dNBR) with an R² of 0.71. Fire severity was highest on drier landforms located on upper hillslopes, ridges, and on southwest aspects, and was higher in pine communities than in other forest types. Predicted CBI was positively correlated with changes in species richness and with the post-fire cover of pine seedlings (Pinus virginiana, P. rigida, and P. pungens), suggesting that burn severity maps can be used to predict community-level fire effects across large landscapes. Despite the relatively large size of this fire for the southern Appalachians, severity was strongly linked to topographic variability and pre-fire vegetation, and spatial variation in fire severity was correlated with changes in species richness. Thus, the Linville Gorge fire appears to have generally reinforced the ecological constraints imposed by underlying environmental gradients.     

A consistency analysis of surface reflectance and leaf area index retrieval from overlapping clear-sky Landsat ETM+ imagery

In this study, the consistency of systematic retrievals of surface reflectance and leaf area index was assessed using overlap regions in adjacent Landsat Enhanced Thematic Mapper-Plus (ETM+) scenes. Adjacent scenes were acquired within 7-25 days apart to minimize variations in the land surface reflectance between acquisition dates. Each Landsat ETM+ scene was independently geo-referenced and atmospherically corrected using a variety of standard approaches. Leaf area index (LAI) models were then applied to the surface reflectance data and the difference in LAI between overlapping scenes was evaluated. The results from this analysis show that systematic LAI retrieval from Landsat ETM+ imagery using a baseline atmospheric correction approach that assumes a constant aerosol optical depth equal to 0.06 is consistent to within ±0.61 LAI units. The average absolute difference in LAI retrieval over all 10 image pairs was 26% for a mean LAI of 2.05 and the maximum absolute difference over any one pair was 61% for a mean LAI of 1.13. When no atmospheric correction was performed on the data, the consistency in LAI retrieval was improved by 1%. When a scene-based dense, dark vegetation atmospheric correction algorithm was used, the LAI retrieval differences increased to 28% for a mean LAI of 2.32. This implies that a scene-based atmospheric correction procedure may improve the absolute accuracy of LAI retrieval without having a major impact on retrieval consistency. Such consistency trials provide insight into the current limits concerning surface reflectance and LAI retrieval from fine spatial resolution remote sensing imagery with respect to the variability in clear-sky atmospheric conditions.     

Monitoring on-orbit calibration stability of the Terra MODIS and Landsat 7 ETM+ sensors using pseudo-invariant test sites

The ability to detect and quantify changes in the Earth's environment depends on sensors that can provide calibrated, consistent measurements of the Earth's surface features through time. A critical step in this process is to put image data from different sensors onto a common radiometric scale. This work focuses on monitoring the long-term on-orbit calibration stability of the Terra Moderate Resolution Imaging Spectroradiometer (MODIS) and the Landsat 7 (L7) Enhanced Thematic Mapper Plus (ETM+) sensors using the Committee on Earth Observation Satellites (CEOS) reference standard pseudo-invariant test sites (Libya 4, Mauritania 1/2, Algeria 3, Libya 1, and Algeria 5). These sites have been frequently used as radiometric targets because of their relatively stable surface conditions temporally. This study was performed using all cloud-free calibrated images from the Terra MODIS and the L7 ETM+ sensors, acquired from launch to December 2008. Homogeneous regions of interest (ROI) were selected in the calibrated images and the mean target statistics were derived from sensor measurements in terms of top-of-atmosphere (TOA) reflectance. For each band pair, a set of fitted coefficients (slope and offset) is provided to monitor the long-term stability over very stable pseudo-invariant test sites. The average percent differences in intercept from the long-term trends obtained from the ETM + TOA reflectance estimates relative to the MODIS for all the CEOS reference standard test sites range from 2.5% to 15%. This gives an estimate of the collective differences due to the Relative Spectral Response (RSR) characteristics of each sensor, bi-directional reflectance distribution function (BRDF), spectral signature of the ground target, and atmospheric composition. The lifetime TOA reflectance trends from both sensors over 10 years are extremely stable, changing by no more than 0.4% per year in its TOA reflectance over the CEOS reference standard test sites.     

Comparison of data gap-filling methods for Landsat ETM+ SLC-off imagery for monitoring forest degradation in a semi-deciduous tropical forest in Mexico

A number of methods to overcome the 2003 failure of the Landsat 7 Enhanced Thematic Mapper (ETM+) scan-line corrector (SLC) are compared in this article in a forest-monitoring application in the Yucatan Peninsula, Mexico. The objective of this comparison is to determine the best approach to accomplish SLC-off image gap-filling for the particular landscape in this region, and thereby provide continuity in the Landsat data sensor archive for forest-monitoring purposes. Four methods were tested: (1) local linear histogram matching (LLHM); (2) neighbourhood similar pixel interpolator (NSPI); (3) geostatistical neighbourhood similar pixel interpolator (GNSPI); and (4) weighted linear regression (WLR). All methods generated reasonable SLC-off gap-filling data that were visually consistent and could be employed in subsequent digital image analysis. Overall accuracy, kappa coefficients (κ), and quantity and allocation disagreement indices were used to evaluate unsupervised Iterative Self-Organizing Data Analysis (ISODATA) land-cover classification maps. In addition, Pearson correlation coefficients (r) and root mean squares of the error (RMSEs) were employed for estimates agreement with fractional land cover. The best results visually (overall accuracy > 85%, κ < 9%, quantity disagreement index < 5.5%, and allocation disagreement index < 12.5%) and statistically (r > 0.84 and RMSE < 7%) were obtained from the GNSPI method. These results suggest that the GNSPI method is suitable for routine use in reconstructing the imagery stack of Landsat ETM+ SLC-off gap-filled data for use in forest-monitoring applications in this type of heterogeneous landscape.     

Landsat ETM+ imaging for mineral potential mapping: application to Avaj area,Qazvin, Iran

The Avaj area is located in the northwestern part of Iran, on the boundary between the Orumieh-Dokhtar and Sanandaj-Sirjan geological zones. It corresponds to two different geological sub-zones:Abe-Garm, to the north, which is related to the volcanic belt of Orumieh-Dokhtar; and the Razan sub-zone, to the south, which is related to the Sanandaj-Sirjan zone. Lithological aspects of different rock units indicate the presence of four distinct volcanic groups of andesite, basalt, tuff, and dacite. The present research utilizes various Landsat Enhanced Thematic Mapper Plus (ETM+) image-processing techniques, including false colour composite images, colour composite ratio images, and standard principal component (PC) image analysis on six bands and colour composite selective PC images, filtering and supervised classification. Nonetheless, colour composite selective images may be the most reliable method for enhancement of areas with hydrothermal alteration. All the techniques used clearly show the pervasive alteration of kaolinite-argillite in the area. These alterations are mainly related to the andesite, and tuffs have the same trend as fault directions.     

Comparison of two types of forest disturbance using multitemporal Landsat TM/ETM+ imagery and field vegetation data

Forest disturbances influence many landscape processes, including changes in microclimate, hydrology, and soil erosion. We analyzed the spectral response and temporal progress of two types of disturbances of spruce forest (bark beetle outbreak and clear-cuts) in the central part of Šumava Mountains at the border between the Czech Republic and Germany, Central Europe. The bark beetle (Ips typographus [L.]) outbreak in this region in the last 20 years resulted in regional-scale spruce forest decay. Clear-cutting was done here to prevent further bark-beetle propagation in the buffer zones.The aim of the study is to identify the differences in spectral response between the two types of forest disturbances and their temporal dynamics. General trends were analyzed throughout the study area, with sampled disturbance areas selected to assess the relationship between field vegetation data and their spectral response. Thirteen Landsat TM/ETM+ scenes from 1985 to 2007 were used for the assessment. The following spectral indices were estimated: NDMI, Tasseled Cap (Brightness, Greenness, Wetness), DI, and DI′. The DI′, Wetness, and Brightness indices show the highest sensitivity to forest disturbance for both disturbance types (clear-cuts and bark beetle outbreak). The multitemporal analysis distinguished three different stages of development. The highest spectral differences between the clear-cuts and the bark beetle disturbances were found in the period between 1996 and 2004 with increased levels of forest disturbance (repeated measures ANOVA, Scheffé post hoc test; p ≤ 0.05). Clear-cut disturbance resulted in significantly higher spectral differences from the original forest and occurred as a more discrete event in comparison to bark beetle outbreak.     

Comparison of burned area estimates derived from SPOT-VEGETATION and Landsat ETM+ data in Africa: Influence of spatial pattern and vegetation type

An algorithm for burned area mapping in Africa based on classification trees was developed using SPOT-VEGETATION (VGT) imagery. The derived 1 km spatial resolution burned area maps were compared with 30 m spatial resolution maps obtained with 13 Landsat ETM+ scenes, through linear regression analysis. The procedure quantifies the bias in burned area estimation present in the low spatial resolution burned area map. Good correspondence was observed for seven sites, with values of the coefficient of determination (R²) ranging from 0.787 to 0.983. Poorer agreement was observed in four sites (R² values between 0.257 and 0.417), and intermediate values of R² (0.670 and 0.613) were obtained for two sites. The observed variation in the level of agreement between the Landsat and VGT estimates of area burned results from differences in the spatial pattern and size distribution of burns in the different fire regimes encompassed by our analysis. Small and fragmented burned areas result in large underestimation at 1 km spatial resolution. When large and compact burned areas dominate the landscape, VGT estimates of burned area are accurate, although in certain situations there is some overestimation. Accuracy of VGT burned area estimates also depends on vegetation type. Results showed that in forest ecosystems VGT maps underestimate substantially the amount of burned area. The most accurate estimates were obtained for woodlands and grasslands. An overall linear regression fitted with the data from the 13 comparison sites revealed that there is a strong relationship between VGT and Landsat estimates of burned area, with a value of R² of 0.754 and a slope of 0.803. Our findings indicate that burned area mapping based on 1 km spatial resolution VGT data provides adequate regional information.     

The relationship between land surface temperature and NDVI with remote sensing: application to Shanghai Landsat 7 ETM+ data

The relationship between land surface temperature (LST) and Normalized Difference Vegetation Index (NDVI) associated with urban land‐use type and land‐use pattern is discussed in the City of Shanghai, China using data collected by the Enhanced Thematic Mapper Plus (ETM+) and aerial photographic remote sensing system. There is an apparent correlation between LST and NDVI from the visual interpretation of LST and NDVI contrasts. Mean LST and NDVI values associated with different land‐use types are significantly different. Multiple comparisons of mean LST and NDVI values associated with pairings of each land‐use type are also shown to be significantly different. The result of a regressive analysis shows an inverse correlation relationship between LST and NDVI within all land‐use polygons, the same to each land‐use type, but correlation coefficients associated with land‐use types are different. An analysis on the relationship between LST, NDVI and Shannon Diversity Index (SHDI) shows a positive correlation between LST and SHDI and a negative correlation between NDVI and SHDI. According to the above results, LST, SHDI and NDVI can be considered to be three basic indices to study the urban ecological environment and to contribute to further validation of the applicability of relatively low cost, moderate spatial resolution satellite imagery in evaluating environmental impacts of urban land function zoning, then to examine the impact of urban land‐use on the urban environment in Shanghai City. This provides an effective tool in evaluating the environmental influences of zoning in urban ecosystems with remote sensing and geographical information systems.     

The use of Meteosat and GMS imagery to detect burned areas in tropical environments

This paper describes a methodology of using data acquired by the European Meteosat and the Japanese Geostationary Meteorological Satellite (GMS) geostationary satellites to detect burned areas in different tropical environments. The methodology is based on a multiple threshold approach applied to the thermal radiance and to a spectral index specific for burned surfaces. The Simple Index for Burned Areas (SIBA), also developed in this study, makes use of the information contained in the visible and thermal InfraRed (IR) band available on the geostationary satellites, whose main advantages are the high temporal resolution and the minimal level of pre-processing required. The results obtained with Meteosat data have been evaluated comparing them with NOAA-Advanced Very High Resolution Radiometer (AVHRR) data acquired over the Central Africa forest-savannah areas. For GMS imagery, AVHRR data acquired over the woodland-savannah areas of Northern Territory in Australia have been used. Despite the very low spatial and spectral resolution of the data, accuracy assessment showed at a regional and continental scale the resulting burned area maps could be a valuable source of information for the monitoring of the fire activity and for the assessment of fire impact on tropospheric chemistry.     

Mapping land cover change in a reindeer herding area of the Russian Arctic using Landsat TM and ETM+ imagery and indigenous knowledge

Much of Russia north of the treeline is grazed by reindeer, and this grazing has materially altered the vegetation cover in many places. Monitoring vegetation change in these remote but ecologically sensitive regions is an important task for which satellite remote sensing is well suited. Further difficulties are imposed by the highly dynamic nature of arctic phenology, and by the difficulty of obtaining accurate official data on land cover in arctic Russia even where such data exist. We have approached the problem in a novel fashion by combining a conventional multispectral analysis of satellite imagery with data on current and historical land use gathered by the techniques of social anthropology, using a study site in the Nenets Autonomous Okrug (NAO). A Landsat-7 ETM+ image from the year 2000 was used to generate a current land cover classification. A Landsat-5 TM image was used to generate a land-cover classification for 1988, taking due account of phenological differences and between the two dates. A cautious comparison of these two classifications, again taking account of possible effects of phenological differences, shows that much of the study area has already undergone a notable transformation to grass-dominated tundra, almost certainly as a result of heavy grazing by reindeer. The grazing pattern is quite heterogeneous, and may have reached unsustainable levels in some areas. Finally, we suggest that this situation is unlikely to be unique to our study area and may well be widespread throughout the Eurasian tundra zone, particularly in the west.     

Mapping of land degradation from space: a comparative study of Landsat ETM+ and ASTER data

The purpose of this study is to compare the role of spectral and spatial resolutions in mapping land degradation from space‐borne imagery using Landsat ETM+ and ASTER data as examples. Land degradation in the form of salinization and waterlogging in Tongyu County, western Jilin Province of northeast China was mapped from an ETM+ image of 22 June 2002 and an ASTER image recorded on 24 June 2001 using supervised classification, together with several other land covers. It was found that the mapping accuracy was achieved at 56.8% and higher for moderately degraded (e.g. salinized) farmland, and over 80% for severely degraded land (e.g. barren) from both ASTER and ETM+ data. The spatial resolution of the ASTER data exerts only a negligible effect on the mapping accuracy. The 30 m ETM+ outperforms the ASTER image of both 15 m and 30 m resolution in consistently generating a higher overall accuracy as well as a higher user's accuracy for barren land. The inferiority of ASTER data is attributed to the highly repetitive spectral content of its six shortwave infrared bands. It is concluded that the spectral resolution of an image is not as important as the information content of individual bands in accurately mapping land covers automatically.     

A water index for rapid mapping of shoreline changes of five East African Rift Valley lakes: an empirical analysis using Landsat TM and ETM+ data

Because of the complicated shorelines, inaccessibility and vast spread of some lakes, information on changing shorelines is difficult to acquire. A new water index (WI) has been applied to quantify changes in five saline and non‐saline Rift Valley lakes in Kenya using Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM+) data. The WI is based on a logical combination of the Tasseled Cap Wetness (TCW) index and the Normalized Difference Water Index (NDWI). Using regression analysis with estimated shoreline coordinates, the WI detected the shorelines with an accuracy of 98.4%, which was 22.3% higher than the TCW, and 43.2% more accurate than the NDWI. Change detection was derived using image differencing followed by density slicing and unsupervised classification. The saline lakes (Bogoria, Nakuru and Elementaita) changed more with respect to the ratio of the change in the original surface areas than the non‐saline lakes (Baringo and Naivasha).     

Comparison of regression and geostatistical methods for mapping Leaf Area Index (LAI) with Landsat ETM+ data over a boreal forest

This study compared aspatial and spatial methods of using remote sensing and field data to predict maximum growing season leaf area index (LAI) maps in a boreal forest in Manitoba, Canada. The methods tested were orthogonal regression analysis (reduced major axis, RMA) and two geostatistical techniques: kriging with an external drift (KED) and sequential Gaussian conditional simulation (SGCS). Deterministic methods such as RMA and KED provide a single predicted map with either aspatial (e.g., standard error, in regression techniques) or limited spatial (e.g., KED variance) assessments of errors, respectively. In contrast, SGCS takes a probabilistic approach, where simulated values are conditional on the sample values and preserve the sample statistics. In this application, canonical indices were used to maximize the ability of Landsat ETM+ spectral data to account for LAI variability measured in the field through a spatially nested sampling design. As expected based on theory, SGCS did the best job preserving the distribution of measured LAI values. In terms of spatial pattern, SGCS preserved the anisotropy observed in semivariograms of measured LAI, while KED reduced anisotropy and lowered global variance (i.e., lower sill), also consistent with theory. The conditional variance of multiple SGCS realizations provided a useful visual and quantitative measure of spatial uncertainty. For applications requiring spatial prediction methods, we concluded KED is more useful if local accuracy is important, but SGCS is better for indicating global pattern. Predicting LAI from satellite data using geostatistical methods requires a distribution and density of primary, reference LAI measurements that are impractical to obtain. For regional NPP modeling with coarse resolution inputs, the aspatial RMA regression method is the most practical option.     

Finer resolution observation and monitoring of global land cover: first mapping results with Landsat TM and ETM+ data

We have produced the first 30 m resolution global land-cover maps using Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) data. We have classified over 6600 scenes of Landsat TM data after 2006, and over 2300 scenes of Landsat TM and ETM+ data before 2006, all selected from the green season. These images cover most of the world's land surface except Antarctica and Greenland. Most of these images came from the United States Geological Survey in level L1T (orthorectified). Four classifiers that were freely available were employed, including the conventional maximum likelihood classifier (MLC), J4.8 decision tree classifier, Random Forest (RF) classifier and support vector machine (SVM) classifier. A total of 91,433 training samples were collected by traversing each scene and finding the most representative and homogeneous samples. A total of 38,664 test samples were collected at preset, fixed locations based on a globally systematic unaligned sampling strategy. Two software tools, Global Analyst and Global Mapper developed by extending the functionality of Google Earth, were used in developing the training and test sample databases by referencing the Moderate Resolution Imaging Spectroradiometer enhanced vegetation index (MODIS EVI) time series for 2010 and high resolution images from Google Earth. A unique land-cover classification system was developed that can be crosswalked to the existing United Nations Food and Agriculture Organization (FAO) land-cover classification system as well as the International Geosphere-Biosphere Programme (IGBP) system. Using the four classification algorithms, we obtained the initial set of global land-cover maps. The SVM produced the highest overall classification accuracy (OCA) of 64.9% assessed with our test samples, with RF (59.8%), J4.8 (57.9%), and MLC (53.9%) ranked from the second to the fourth. We also estimated the OCAs using a subset of our test samples (8629) each of which represented a homogeneous area greater than 500 m?×?500 m. Using this subset, we found the OCA for the SVM to be 71.5%. As a consistent source for estimating the coverage of global land-cover types in the world, estimation from the test samples shows that only 6.90% of the world is planted for agricultural production. The total area of cropland is 11.51% if unplanted croplands are included. The forests, grasslands, and shrublands cover 28.35%, 13.37%, and 11.49% of the world, respectively. The impervious surface covers only 0.66% of the world. Inland waterbodies, barren lands, and snow and ice cover 3.56%, 16.51%, and 12.81% of the world, respectively.     

Characterizing changes in cropping patterns using sequential Landsat imagery: an adaptive threshold approach and application to Phoenix,Arizona

Since the 1970s, the Phoenix Active Management Area has experienced rapid urbanization, mostly through land conversions from agricultural lands to urban land use. Rapid urban expansion and population growth have placed unprecedented pressure on agricultural production in this region. Agricultural intensification, in particular double cropping, has been observed globally as an important response to the growing pressure on land. However, the intensification has a number of negative impacts on water quality, biodiversity, and biogeochemical cycles. Thus, quantifying the spatial pattern of cropping intensity is important for natural resource management. In this study, we developed an adaptive threshold approach to map cropping intensity using time series Landsat data and examined the spatiotemporal patterns of cropping intensity in the Phoenix Active Management Area from 1995 to 2010 at 5-year intervals. To map cropping intensity accurately, the adaptive threshold algorithm was designed specifically to address several issues caused by the complex cropping patterns in the study area. The adaptive threshold method has abilities to (1) distinguish true crop cycles from multiple false phenological peaks, (2) minimize errors caused by data noise and missing data, (3) identify alfalfa and interyear crops and to distinguish alfalfa from double crops, and (4) adapt to temporal profiles with different numbers of observations. The adaptive threshold algorithm is effective in characterizing cropping intensity with overall accuracies exceeding 97%. Results show that there is a dramatic decline in the area of total croplands (46.1%), single crops (46.3%), and double crops (43.4%) during the study period. There was a small conversion (1.9%) from single to double crop from 1995 to 2000, whereas a reverse conversion (1.3%) was observed from 2005 to 2010. Updated and accurate information on the spatial distribution of cropping intensity provide important implications on effective and sustainable cropping practices. In addition, joint investigation on cropping patterns and irrigation water use can shed light on future agricultural water demand, which is of paramount importance in this rapidly expanding arid region.