Landsat-7 ETM+ radiometric normalization comparison for northern mapping applications

Relative radiometric normalization has long been performed to generate consistency among individual Landsat scenes for production of composites containing multiple scenes. Normalization methods have relied on matching identical and assumed invariant features in both images of an overlapping pair, or on invariant targets that are not necessarily the same features. Problems with overlap normalization methods include sensitivity to outliers in overlap data caused by atmospheric or land cover change between scenes, which can lead to radiometric error propagation across a mosaic caused by a normalized scene becoming a reference for the subsequent scene entered into the mosaic. Solutions to such problems include interactive outlier removal to generate a normalization function using a ‘no change’ data set and methods that are robust against outliers to automatically generate normalization functions with minimal user input. This paper compares two normalization methods that use a robust regression technique called Theil-Sen with an established overlap normalization method. The first method uses Theil-Sen regression to generate a normalization function between overlap regions, while the second uses Theil-Sen to normalize to coarse-resolution composite reflectance data from the SPOT VEGETATION (VGT) sensor. The results of the normalizations were evaluated in two ways: (1) using statistics generated between overlap regions; and (2) separately using coarse-resolution data as a reference. Both overlap normalization methods performed almost identically; however, Theil-Sen was faster and easier to implement than its traditional counterpart due to its insensitivity to outliers and capability for full automation. While overlap and coarse-resolution normalizations each outperformed the other when evaluated against its calibration set, error propagation caused by outliers in overlap samples was avoided in the normalization to coarse-resolution imagery. Advantages offered by normalization to coarse-resolution data using robust regression, including full automation, make this method particularly attractive for generation of large area mosaics containing 100 Landsat scenes or more.     

Estimating tropical forest biomass more accurately by integrating ALOS PALSAR and Landsat-7 ETM+ data

Integration of multisensor data provides the opportunity to explore benefits emanating from different data sources. A fusion between fraction images derived from spectral mixture analysis of Landsat-7 ETM+ and phased array L-band synthetic aperture radar (PALSAR) is introduced. The aim of this fusion is to improve the estimation accuracy of above-ground biomass (AGB) in lowland mixed dipterocarp forest. Spectral mixture analysis was applied to decompose a mixture of spectral components of Landsat-7 ETM+ into vegetation, soil, and shade fractions. These fraction images were integrated with PALSAR data using the discrete wavelet transform (DWT) and Brovey transform. As a comparison, spectral reflectance of Landsat-7 ETM+ was fused directly with PALSAR data. Backscatter of horizontal–horizontal and horizontal–vertical polarizations was also used to estimate AGB. Forest inventory was carried out in 77 randomly distributed plots, the data being used for either model development or validation. A local allometric equation was applied to calculate AGB per plot. Regression models were developed by integrating field measurements of 50 sample plots with remotely sensed data, e.g. fraction images, reflectance of Landsat-7 ETM+, and PALSAR data. The models developed were validated using 27 independent sample plots. The results showed that not all fused images significantly improved the accuracy of AGB estimation. The model based on Brovey transform using the reflectance of Landsat-7ETM+ and PALSAR produced an R² of only 0.03–0.10. By contrast, fusion between PALSAR data and fraction images using Brovey transform improved the accuracy of R² to 0.33–0.46. Further improvement in the accuracy of estimating AGB was observed when DWT was applied to integrate PALSAR with the reflectance of Landsat-7ETM+ (R² = 0.69–0.72) and PALSAR with fraction images (R² = 0.70–0.75).     

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.     

Retrieval of Component Temperatures based on Landsat-7 ETM+ Remote Sensing Data

Land surface temperature plays an important role in drought monitoring and Simulation of surface heat flux.In arid and semi\|arid regions,the Two\|Source Energy Balance model (TSEB) is commonly used to calculate the heat flux of the earth’s surface.Taking the typical irrigated area of the middle reaches of Heihe as the research area,the four Landsat\|7 ETM+ remote sensing images are selected.The soil surface temperature and canopy temperature were retrieved by combining vegetation index with TSEB model.The decomposition algorithm of soil surface temperature and vegetation canopy temperature is mainly discussed.The results showed that soil surface temperature and vegetation canopy temperature had good temporal and spatial consistency.The inversion accuracy of soil surface temperature and vegetation canopy temperature is indirectly verified by surface net radiation and surface heat flux.The calculated values of surface net radiation and surface heat flux correlate well with the observed values,and the correlation coefficient is greater than 0.92.The linear regression analysis of surface net radiation and surface heat flux shows that the fitting accuracy is high.The soil surface temperature and canopy temperature obtained by surface temperature decomposition are feasible for monitoring drought in typical areas and simulating surface heat flux.     

Using lidar and effective LAI data to evaluate IKONOS and Landsat 7 ETM+ vegetation cover estimates in a ponderosa pine forest

Structural and functional analyses of ecosystems benefit when high accuracy vegetation coverages can be derived over large areas. In this study, we utilize IKONOS, Landsat 7 ETM+, and airborne scanning light detection and ranging (lidar) to quantify coniferous forest and understory grass coverages in a ponderosa pine (Pinus ponderosa) dominated ecosystem in the Black Hills of South Dakota. Linear spectral mixture analyses of IKONOS and ETM+ data were used to isolate spectral endmembers (bare soil, understory grass, and tree/shade) and calculate their subpixel fractional coverages. We then compared these endmember cover estimates to similar cover estimates derived from lidar data and field measures. The IKONOS-derived tree/shade fraction was significantly correlated with the field-measured canopy effective leaf area index (LAI_e) (r²=0.55, p<0.001) and with the lidar-derived estimate of tree occurrence (r²=0.79, p<0.001). The enhanced vegetation index (EVI) calculated from IKONOS imagery showed a negative correlation with the field measured tree canopy effective LAI and lidar tree cover response (r²=0.30, r=−0.55 and r²=0.41, r=−0.64, respectively; p<0.001) and further analyses indicate a strong linear relationship between EVI and the IKONOS-derived grass fraction (r²=0.99, p<0.001). We also found that using EVI resulted in better agreement with the subpixel vegetation fractions in this ecosystem than using normalized difference of vegetation index (NDVI). Coarsening the IKONOS data to 30 m resolution imagery revealed a stronger relationship with lidar tree measures (r²=0.77, p<0.001) than at 4 m resolution (r²=0.58, p<0.001). Unmixed tree/shade fractions derived from 30 m resolution ETM+ imagery also showed a significant correlation with the lidar data (r²=0.66, p<0.001). These results demonstrate the power of using high resolution lidar data to validate spectral unmixing results of satellite imagery, and indicate that IKONOS data and Landsat 7 ETM+ data both can serve to make the important distinction between tree/shade coverage and exposed understory grass coverage during peak summertime greenness in a ponderosa pine forest ecosystem.     

Modelling and mapping some soil surface properties of Central Kelkit Basin in Turkey by using Landsat-7 ETM+ images

Using remote sensing and geographical information system (GIS) tools, we investigated relationships among Landsat-7 Enhanced Thematic Mapper Plus (Landsat-7 ETM+) satellite bands and some soil surface variables of Central Kelkit Basin in Turkey. We collected 164 geo-referenced surface soil samples (0–20 cm) from field studies in 2008 and analysed the data to determine soil variables such as calcium carbonate (CaCO₃), pH, electrical conductivity (EC), organic matter (OM), nitrogen (N), phosphorus (P), exchangeable potassium (EK), and texture (clay, silt, and sand). Utilizing geographic references and determined soil variables, we established a point database in GIS to overlay and extract corresponding digital number (DN) values of Landsat-7 ETM+ bands. We studied possible relationships between the soil variables and the bands' DN values by using correlation (Pearson) analysis. We found that four soil variables including pH, OM, CaCO₃, and N have significant correlations with band 5 (short wave infrared) DN values. We conducted linear regression analysis to develop models for these four variables. Residual statistics and plots showed that the developed models were robust. After running the models in GIS, we mapped pH, CaCO₃, OM, and N contents of the basin.     

Hyperion, IKONOS, ALI, and ETM+ sensors in the study of African rainforests

The goal of this research was to compare narrowband hyperspectral Hyperion data with broadband hyperspatial IKONOS data and advanced multispectral Advanced Land Imager (ALI) and Landsat-7 Enhanced Thematic Mapper Plus (ETM+) data through modeling and classifying complex rainforest vegetation. For this purpose, Hyperion, ALI, IKONOS, and ETM+ data were acquired for southern Cameroon, a region considered to be a representative area for tropical moist evergreen and semi-deciduous forests. Field data, collected in near-real time to coincide with satellite sensor overpass, were used to (1) quantify and model the biomass of tree, shrub, and weed species; and (2) characterize forest land use/land cover (LULC) classes.The study established that even the most advanced broadband sensors (i.e., ETM+, IKONOS, and ALI) had serious limitations in modeling biomass and in classifying forest LULC classes. The broadband models explained only 13-60% of the variability in biomass across primary forests, secondary forests, and fallows. The overall accuracies were between 42% and 51% for classifying nine complex rainforest LULC classes using the broadband data of these sensors. Within individual vegetation types (e.g., primary or secondary forest), the overall accuracies increased slightly, but followed a similar trend. Among the broadband sensors, ALI sensor performed better than the IKONOS and ETM+ sensors.When compared to the three broadband sensors, Hyperion narrowband data produced (1) models that explained 36-83% more of the variability in rainforest biomass, and (2) LULC classifications with 45-52% higher overall accuracies. Twenty-three Hyperion narrowbands that were most sensitive in modeling forest biomass and in classifying forest LULC classes were identified and discussed.     

Change detection of coral reef habitat using Landsat-5 TM,Landsat 7 ETM+ and Landsat 8 OLI data in the Red Sea (Hurghada,Egypt)

Owing to continuing touristic developments in Hurghada, Egypt, several coral reef habitats have suffered major deterioration between 1987 and 2013, either by being bleached or totally lost. Such alterations in coral reef habitats have been well observed in their varying distributions using change detection analysis applied to a Landsat 5 image representing 1987, a Landsat 7 image representing 2000, and a Landsat 8 image representing 2013. Different processing techniques were carried out over the three images, including but not limited to rectification, masking, water column correction, classification, and change detection statistics. The supervised classifications performed over the three scenes show five significant marine-related classes, namely coral, sand subtidal, sand intertidal, macro-algae, and seagrass, in different degrees of abundance. The change detection statistics obtained from the classified scenes of 1987 and 2000 reveal a significant increase in the macro-algae and seagrass classes (93 and 47%, respectively). However, major decreases of 41, 40, and 37% are observed in the sand intertidal, coral, and sand subtidal classes, respectively. On the other hand, the change detection statistics obtained from the classified scenes of 2000 and 2013 revealed increases in sand subtidal and macro-algae classes by 14 and 19%, respectively, while major decreases of 49%, 46% and 74% are observed in the sand intertidal, coral, and seagrass classes, respectively.     

A comparison of gap-filling approaches for Landsat-7 satellite data

The purpose of this study is to assess the relative performance of four different gap-filling approaches across a range of land-surface conditions, including both homogeneous and heterogeneous areas as well as in scenes with abrupt changes in landscape elements. The techniques considered in this study include: (1) Kriging and co-Kriging; (2) geostatistical neighbourhood similar pixel interpolator (GNSPI); (3) a weighted linear regression (WLR) algorithm; and (4) the direct sampling (DS) method. To examine the impact of image availability and the influence of temporal distance on the selection of input training data (i.e. time separating the training data from the gap-filled target image), input images acquired within the same season (temporally close) as well as in different seasons (temporally far) to the target image were examined, as was the case of using information only within the target image itself. Root mean square error (RMSE), mean spectral angle (MSA), and coefficient of determination (R²) were used as the evaluation metrics to assess the prediction results. In addition, the overall accuracy (OA) and kappa coefficient (kappa) were used to assess a land-cover classification based on the gap-filled images. Results show that all of the gap-filling approaches provide satisfactory results for the homogeneous case, with R² > 0.93 for bands 1 and 2 in all cases and R² > 0.80 for bands 3 and 4 in most cases. For the heterogeneous example, GNSPI performs the best, with R² > 0.85 for all tested cases. WLR and GNSPI exhibit equivalent accuracy when a temporally close input image is used (i.e. WLR and GNSPI both have an R² equal to 0.89 for band 1). For the case of abrupt changes in scene elements or in the absence of ancillary data, the DS approach outperforms the other tested methods.     

Mapping rice field anopheline breeding habitats in Mali,West Africa,using Landsat ETM+ sensor data

The aim of this study was to determine whether remotely sensed data could be used to identify rice-related malaria vector breeding habitats in an irrigated rice growing area near Niono, Mali. Early stages of rice growth show peak larval production, but Landsat sensor data are often obstructed by clouds during the early part of the cropping cycle (rainy season). In this study, we examined whether a classification based on two Landsat Enhanced Thematic Mapper (ETM)+ scenes acquired in the middle of the season and at harvesting times could be used to map different land uses and rice planted at different times (cohorts), and to infer which rice growth stages were present earlier in the season. We performed a maximum likelihood supervised classification and evaluated the robustness of the classifications with the transformed divergence separability index, the kappa coefficient and confusion matrices. Rice was distinguished from other land uses with 98% accuracy and rice cohorts were discriminated with 84% accuracy (three classes) or 94% (two classes). Our study showed that optical remote sensing can reliably identify potential malaria mosquito breeding habitats from space. In the future, these ‘crop landscape maps' could be used to investigate the relationship between cultivation practices and malaria transmission.     

Mapping rice field anopheline breeding habitats in Mali, West Africa, using Landsat ETM+ sensor data

The aim of this study was to determine whether remotely sensed data could be used to identify rice-related malaria vector breeding habitats in an irrigated rice growing area near Niono, Mali. Early stages of rice growth show peak larval production, but Landsat sensor data are often obstructed by clouds during the early part of the cropping cycle (rainy season). In this study, we examined whether a classification based on two Landsat Enhanced Thematic Mapper (ETM)+ scenes acquired in the middle of the season and at harvesting times could be used to map different land uses and rice planted at different times (cohorts), and to infer which rice growth stages were present earlier in the season. We performed a maximum likelihood supervised classification and evaluated the robustness of the classifications with the transformed divergence separability index, the kappa coefficient and confusion matrices. Rice was distinguished from other land uses with 98% accuracy and rice cohorts were discriminated with 84% accuracy (three classes) or 94% (two classes). Our study showed that optical remote sensing can reliably identify potential malaria mosquito breeding habitats from space. In the future, these 'crop landscape maps' could be used to investigate the relationship between cultivation practices and malaria transmission.     

Modelling relationships between NDVI and precipitation during vegetative growth cycles

This paper describes how NOVI responds to daily precipitation for vegetation under water-limited growth conditions (e.g., in arid or semi-arid environments). An analytical model was developed to address timing to peak NOVI response and total duration of NOVI response to a rainfall event with the passage of time. The model also describes how effectively the vegetation uses precipitated water. A case study in the Sandhills area of Nebraska, with NOVI calculated from Landsat-MSS images, showed that the model can accurately portray NOVI time series, and then can be used to derive some useful information about relations of NOVI and precipitation. For the vegetation in the test area, we found that the peak vegetative response to a precipitation event shifts from about 14 days after the rainfall event at the beginning of a growing season to 25 days at the late stage of the season, and then shifts back to 12 days at the end of the season. The duration of response, similar to the peak response in behavior, is 22, 61, and 25 days, respectively.     

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.     

Estimating aboveground biomass using Landsat 7 ETM+ data across a managed landscape in northern Wisconsin, USA

Aboveground biomass (AGB; Mg/ha) is defined in this study as a biomass of growing stock trees greater than 2.5 cm in diameter at breast height (dbh) for stands >5 years and all trees taller than 1.3 m for stands <5 years. Although AGB is an important variable for evaluating ecosystem function and structure across the landscape, such estimates are difficult to generate without high-resolution satellite data. This study bridges the application of remote sensing techniques with various forest management practices in Chequamegon National Forest (CNF), Wisconsin, USA by producing a high-resolution stand age map and a spatially explicit AGB map. We coupled AGB values, calculated from field measurements of tree dbh, with various vegetation indices derived from Landsat 7 ETM+ data through multiple regression analyses to produce an initial biomass map. The initial biomass map was overlaid with a land-cover map to generate a stand age map. Biomass threshold values for each age category (e.g., young, intermediate, and mature) were determined through field observations and frequency analysis of initial biomass estimates by major cover types. We found that AGB estimates for hardwood forests were strongly related to stand age and near-infrared reflectance (r²=0.95) while the AGB for pine forests was strongly related to the corrected normalized difference vegetation index (NDVIc; r²=0.86). Separating hardwoods from pine forests improved the AGB estimates in the area substantially, compared to overall regression (r²=0.82). Our AGB results are comparable to previously reported values in the area. The total amount of AGB in the study area for 2001 was estimated as 3.3 million metric tons (dry weight), 76.5% of which was in hardwood and mixed hardwood/pine forests. AGB ranged from 1 to 358 Mg/ha with an average of 70 and a standard deviation of 54 Mg/ha. The AGB class with the highest percentage (16.1%) was between 81 and 100 Mg/ha. Forests with biomass values >200 Mg/ha accounted for less than 3% of the study area and were usually associated with mature hardwood forests. Estimated AGB was validated using independent field measurements (R²=0.67, p<0.001). The AGB and age maps can be used as baseline information for future landscape level studies such as quantifying the regional carbon budget, accumulating fuel, or monitoring management practices.     

Evaluation of gap-filling methods for Landsat 7 ETM+ SLC-off image for LULC classification in a heterogeneous landscape of West Africa

ABSTRACT

The Landsat mission which has existed over five decades has remained at the forefront of providing consistent moderate spatial and temporal resolution optical images of the earth. The failure of the scan line corrector (SLC) on board the Landsat 7 Enhanced Thematic Mapper Plus (ETM+) in May 2003 has permanently resulted in data gaps on each Landsat 7 scene. Due to the obvious negative impacts on the image usability, a number of methods have been developed to fill the no-data areas in the image. This study assessed the performance of four Landsat 7 ETM+ SLC-off gap-filling methods in a highly heterogeneous landscape of West Africa for two different seasons (dry and rainy). The methods considered are: (1) Weighted Linear Regression (WLR) integrated with Laplacian Prior Regularization Method (LPRM), (2) Localised Linear Histogram Matching (LLHM), (3) Neighbourhood Similar Pixel Interpolator (NSPI) and (4) Geostatistical Neighbourhood Similar Pixel Interpolator (GNSPI). All the images used were Landsat 7 ETM+ SLC-off images, temporally close and from the same season for each set of time step. Visual comparison, mean, and standard deviations of the histograms of all bands of only the filled areas were used to assess the results. Additionally, overall accuracy (OA), kappa coefficient (κ), and balanced accuracy (BA) per class were used to evaluate a land use/cover (LULC) classification based on the gap-filled images. Visually, all the four methods were able to completely fill the gaps in the Landsat 7 ETM+ SLC-off image. They all look similar and spatially continuous with no anomalies or artefacts on them. The histograms from each band for only the filled areas for all the four methods also gave similar means and standard deviations in most cases. All the four gap-filling methods provided satisfactory results (OA >96% and κ> 0.937 in all methods for images in the dry season and OA >93% and κ> 0.877 for the image in the rainy season) in the land cover classification considering the complexity of the study area. But the GNSPI was superiority in all cases with the highest OA of 97.1% and κ of 0.947 in the dry season and OA of 94.6% and κ of 0.899 in the rainy season. This implies that the GNSPI is more robust in gap filling of Landsat 7 ETM+ SLC-off images than the other three methods in a heterogeneous landscape of West Africa regardless of the season. This study suggests that gap filling of Landsat 7 ETM+ SLC-off images will help to increase the number of Landsat images needed to build time-series data for a data-scarce region such as West Africa.     

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.     

Evaluation of dynamic linkages between evapotranspiration and land-use/land-cover changes with Landsat TM and ETM+ data

Complexity embedded in coastal management leads to numerous questions as to how inherent spatial and temporal linkages among evapotranspiration (ET), depth to groundwater and land-use/land-cover change (LUCC) could affect the dynamics among these seemingly unrelated events. This article aims to address such unique dynamics in the nexus of physical geography and ecohydrology. To understand such dynamic linkages, a case study was carried out in a fast growing coastal region – the southern Laizhou Bay in Shandong Province, China – by identifying the coastal LUCC at the decadal scale in association with the variations of ET with the aid of Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM+) data. In such a coastal landscape evolutionary assessment, findings show that the major patterns of land use and land cover (LULC) in the study area are farmland, saline-alkali land, developed land, salt land and beach land. Over a 20-year time frame, declining groundwater trends were observed, while ET increased gradually with changing LULC. By using the surface energy balance algorithm for land (SEBAL) with Landsat TM/ETM+ images and additional environmental data, the concomitant response of ET variations due to LUCC becomes lucid among three significantly correlated pairs including fractional vegetation cover (FVC), land surface temperature (LST) and soil heat flux. The dynamic linkages between ET and LULC were finally confirmed with such a pair-wise analysis.     

Comparing MODIS and ETM+ data for regional and global land classification

Nearly simultaneous reflectance data sets from the Landsat 7 Enhanced Thematic Mapper Plus (ETM+), at 30-m resolution, and the Terra satellite instrument MODIS, at 500-m resolution, are compared for their ability to map fractional coverage of surface types over large areas. Lower spatial resolution MODIS classification results are generally comparable those of ETM+, with discrepancies for some regions with mixed surface types. Analysis of laboratory and field spectra suggests an ambiguity, the “brightness ambiguity”, which can prevent accurate area estimation of pixels having two or more surface types. This ambiguity, plus general mathematical inversion issues, can account for the discrepancy. Thus, occasional high-resolution measurements, as from Landsat 7, are necessary to refine estimations of large area surface types from MODIS and similar lower spatial resolution instruments.     

Modelling forest stand volume and tree density using Landsat ETM+ data

Estimation of stand volume and tree density in a large area using remotely sensed data has considerable significance for sustainable management of natural resources. In this paper, we explore likely relationships between forest stand characteristics and Landsat Enhanced Thematic Mapper Plus (ETM+) reflectance values. We used multivariate regression technique to predict stand volume and tree density. The result showed that a linear combination of greenness and difference vegetation index (DVI) were better predictors of stand volume (adjusted R² = 43%; root mean square error (RMSE) = 97.4 m³ ha^?1) than other ETM+ bands and vegetation indices. In addition, the regression model with ETM4 (near infrared band) and ETM5 (first shortwave band) as independent variables was a better predictor of tree density (adjusted R² = 73.4%; RMSE = 170.13 ha^?1) than other combinations of ETM⁺ bands and vegetation indices. Results obtained from this study demonstrate the significant relationship between forest stand characteristics and ETM+ reflectance values and the utility of transformed bands in modelling stand volume and tree density. Based on the results of this study, we conclude that ETM+ data are useful to estimate forest volume and density and to gain insights into its structural characteristics in our study area. Forest managers could use ETM+ data for gaining insights into stand characteristics and generating maps required for developing forest management plans and identifying locations within stands that require treatments and other interventions.     

Identifying ring geomorphic features in Lebanon using ASTER and Landsat 7 ETM+ images

Satellite images enable us to identify several geomorphologic features on the Earth's surface that could not be easily recognized on ground surface or by using conventional methods. This is mainly attributed to the optical advantages of remote sensing techniques. Thus, suspicious geomorphologic signatures can be observed on the terrain surface. These are mostly geologic‐controlled. Linear aspects are given most attention in many geological studies to reflect subsurface, hidden structures. In this study, using ENVI‐4.3 and ERDAS Imagine‐9.1 software to analyse ASTER and Landsat 7 ETM+ images of Lebanon, we exposed a miscellany of geomorphic ring structures on different geologic formations of the Lebanese terrain. The diameter of the recognized ring structures ranges from a few hundreds of metres to several kilometres. Preliminary field surveys were carried out on some of these structures in order to identify their origin. All the recognized structures are circular and semi‐rounded in their geomorphology; however, they are attributed either to subsurface crust processes or to the impact of falling meteoritic from outer space. The former were indicated through rock deformations and existence of intrusive bodies. Meteorite impact was induced from the glassy and metallic materials prevailing within the ring structure area, as well as from the interference of some rings.