Multitemporal assessment of the geomorphologic evolution of the Restinga of Marambaia,Rio de Janeiro,Brazil

The Restinga of Marambaia is an emerged sand bar located between the Sepetiba Bay and the South Atlantic Ocean, on the south‐east coast of Brazil. The objective of this study was to observe the geomorphologic evolution of the coastal zone of the Restinga of Marambaia using multitemporal satellite images acquired by multisensors from 1975 to 2004. The images were digitally segmented by a region growth algorithm and submitted to an unsupervised classification procedure (ISOSEG) followed by a raster edit based on visual interpretation. The image time‐series showed a general trend of decrease in the total sand bar area with values varying from 80.61 km² in 1975 to 78.15 km² in 2004. The total area calculation based on the 1975 and 1978 Landsat MSS data was shown to be super‐estimated in relation to the Landsat TM, Landsat ETM+, and CBERS‐2 CCD data. These differences can also be associated to the relatively poorer spatial resolution of the MSS data, nominally 79 m, against the 20 m of the CCD data and 30 m of the TM and ETM+ data. For the estimates of the width in the central portion of the sand bar the variation was from 158 m (1975) to 100 m (2004). The formation of a spit in the northern region of the study area was visually observed. The area of the spit was estimated, with values varying from 0.82 km² (1975) to 0.55 km² (2004).     

Changes in the glaciers of Chandra–Bhaga basin,Himachal Himalaya,India, between 1980 and 2010 measured using remote sensing

This study reports the glacier changes of Chandra–Bhaga basin, northwest Himalaya, India, from 1980 to 2010. Satellite remote-sensing data from the Landsat Multispectral Scanner (MSS) and Thematic Mapper (TM), the Linear Imaging Self Scanning Sensor (LISS) and Advanced Wide Field Sensor (AWiFS) of the Indian Remote Sensing (IRS) series, and the Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) were used to study the changes in glacier parameters such as glacier area, length, snout elevation, and the impact of glacier topographical parameters (glacier slope, aspect, and altitude range) on the glacier changes. It was found that the total glaciated area had shrunk to 368.2 km² in 2010 from 377.6 km² in 1980, a loss of 2.5%. The average position of glacier terminuses retreated by 465.5 ± 169.1 m from 1980 to 2010 with an average rate of 15.5 ± 5.6 m year^?1. The decadal scale analysis showed that the average rate of retreat had increased the most in the recent decade. A moraine-dammed lake located in the study region was found to have expanded in area from (0.65 ± 0.01) km² in 1980 to (1.26 ± 0.03) km² in 2010. Glaciers with steep slope and less altitude range have lost more area than the glaciers having gentle slope and greater altitude range.     

Frontal changes in the Manimahesh and Tal Glaciers in the Ravi basin,Himachal Pradesh,northwestern Himalaya (India), between 1971 and 2013

The Manimahesh and Tal Glaciers are located in the Budhil fifth-order sub-basin of the Ravi, Himachal Himalaya, Northwestern Himalaya (India). These glaciers were analysed using high- (Corona KH-4A) to medium- (Landsat TM/ETM+/OLI, ASTER) spatial resolution satellite data between 1971 and 2013, along with extensive field measurements (2011–2014) of frontal changes. The results show that the Manimahesh and Tal Glaciers retreated by 157 ± 34 m (4 ± 1 m year^–1) and 45 ± 34 m (1 ± 1 m year^–1), respectively, whereas, the total area lost is estimated at 0.21 ± 0.01 km² (0.005 km² year^–1) and 0.010 ± 0.003 km² (0.0002 km² year^–1), respectively, between 1971 and 2013. The rate of retreat is significantly lower than that previously reported. Our field measurements (2011–2014) also suggest a retreating trend and validate the measured glacier changes using remotely sensed temporal data.     

Changes in snow and glacier cover in an arid watershed of the western Kunlun Mountains using multisource remote-sensing data

Snow and glaciers in the mountain watersheds of the Tarim River basin in western China provide the primary water resources to cover the needs of downstream oases. Remote sensing provides a practical approach to monitoring the change in snow and glacier cover in those mountain watersheds. This study investigated the change in snow and glacier cover in one such mountain watershed using multisource remote-sensing data, including the Moderate Resolution Imaging Spectroradiometer (MODIS), Landsat (Multispectral Scanner (MSS), Thematic Mapper (TM), and Enhanced Thematic Mapper Plus (ETM+)), Corona, and Google Earth^TM imagery. With 10 years’ daily MODIS snow-cover data from 2002 to 2012, we used two de-cloud methods before calculating daily snow-cover percentage (SCP), annual snow-cover frequency (SCF), and annual minimum snow-cover percentage (AMSCP) for the watershed. Mann–Kendall analysis showed no significant trend in any of those snow-cover characterizations. With a total of 22 Landsat images from 1967 to 2011, we used band ratio and supervised classification methods for snow classification for Landsat TM/ETM+ images and MSS images, respectively. The Landsat snow-cover data were divided into two periods (1976–2002 and 2004–2011). Statistical tests indicated no significant difference in either the variance or mean of SCPs between the two periods. Three glaciers were identified from Landsat images of 1998 and 2011, and their total area increased by 12.6%. In addition, three rock glaciers were also identified on both the Corona image of 1968 and the Google high-resolution image of 2007, and their area increased by 2.5%. Overall, based on multisource remote-sensing data sets, our study found no evidence of significant changes in snow and glacier cover in the watershed.     

Assessment of deforestation in the Lower Amazon floodplain using historical Landsat MSS/TM imagery

The floodplain forests bordering the Amazon River have outstanding ecological, economic, and social importance for the region. However, the original distribution of these forests is not well known, since they have suffered severe degradation since the 16th century. The previously published vegetation map of the Amazon River floodplain (Hess et al., 2003), based on data acquired in 1996, shows enormous difference in vegetation cover classes between the regions upstream and downstream of the city of Manaus. The upper floodplain is mostly covered by forests, while the lower floodplain is predominantly occupied by grasses and shrubs.This study assesses deforestation in the Lower Amazon floodplain over a ~ 30 year period by producing and comparing a historical vegetation map based on MSS/Landsat images acquired in the late 1970s with a recent vegetation map produced from TM/Landsat images obtained in 2008. The maps were generated through the following steps: 1) normalization and mosaicking of images for each decade; 2) application of a linear mixing model transformation to produce vegetation, soil and shade fraction-images; and 3) object-oriented image analysis and classification. For both maps, the following classes were mapped: floodplain forest, non-forest floodplain vegetation, bare soil and open water. The two maps were combined using object-level Boolean operations to identify time transitions among the mapped classes, resulting in a map of the land cover change occurred over ~ 30 years. Ground information collected at 168 ground points was used to build confusion matrices and calculate Kappa indices of agreement. A survey strategy combining field observations and interviews allowed the collection of information about both recent and historical land cover for validation purposes. Kappa values (0.77, 0.75 and 0.75) indicated the good quality of the maps, and the error estimates were used to adjust the estimated deforested area to a value of 3457 km² ± 1062 km² (95% CI) of floodplain deforestation over the ~ 30 years.     

Near‐infrared discrimination of leafless saltcedar in wintertime Landsat TM

To test a hypothesis that leafless riparian canopies enable accurate multi‐spectral discrimination of saltcedar (Tamarix ramosissima Ledeb.) from other native species, winter Landsat TM5 data (16 November 2005) were analysed for a reach of the Arkansas River in Colorado, USA. Supporting spectroscopic analysis confirmed that saltcedar could not easily be discriminated from other riparian vegetation using TM5 data when in‐leaf, but bare branches could be easily distinguished due to much lower reflectance than other riparian cover. Use of TM Band 4 (B4) allowed differentiation of wintertime saltcedar into four qualitative density classes judged from high‐resolution low‐oblique aerial photography: high (76%–100%), medium (51%–75%), low (16%–50%), and none (0%–15%). Spectral overlap was removed from the B4 saltcedar classification using TM Band 5 (B5) thresholds to eliminate low‐reflectant wet areas and higher‐reflectant multi‐year darkened weed canopies. The accuracy of a classification algorithm that used B5 thresholds followed by a B4 density slice was judged against high‐resolution aerial photography as providing 98% discrimination of saltcedar cover from other riparian cover and about 90% discrimination of the qualitative density classes. Applying this method to the 2835 km² riparian corridor study area, 1298 km² (45.78%) was identified as containing saltcedar, with over 43% having medium or greater density.     

Analysis of land use/cover changes and urban expansion of Nairobi city using remote sensing and GIS

We used three Landsat images together with socio‐economic data in a post‐classification analysis to map the spatial dynamics of land use/cover changes and identify the urbanization process in Nairobi city. Land use/cover statistics, extracted from Landsat Multi‐spectral Scanner (MSS), Thematic Mapper (TM) and Enhanced Thematic Mapper plus (ETM+) images for 1976, 1988 and 2000 respectively, revealed that the built‐up area has expanded by about 47?km². The road network has influenced the spatial patterns and structure of urban development, so that the expansion of the built‐up areas has assumed an accretive as well as linear growth along the major roads. The urban expansion has been accompanied by loss of forests and urban sprawl. Integration of demographic and socio‐economic data with land use/cover change revealed that economic growth and proximity to transportation routes have been the major factors promoting urban expansion. Topography, geology and soils were also analysed as possible factors influencing expansion. The integration of remote sensing and Geographical Information System (GIS) was found to be effective in monitoring land use/cover changes and providing valuable information necessary for planning and research. A better understanding of the spatial and temporal dynamics of the city's growth, provided by this study, forms a basis for better planning and effective spatial organization of urban activities for future development of Nairobi city.     

Subpixel measurement of tropical forest cover using AVHRR data

Abstract

Tropical forest assessment using data from the Advanced Very High Resolution Radiometer (AVHRR) may lead to inaccurate estimates of forest cover in regions of small subpixel forest or non-forest patches and in regions where the pattern of clearance is particularly convoluted. Test sites typifying these two patterns were chosen in Ghana and Rondonia, respectively. To capture the subpixel proportions of forest cover, a linear mixture model was applied to two AVHRR test images over the test sites. The model produced image outputs in which pixel intensities indicated the proporton of forest cover per km². For comparison, supervised maximum likelihood classifications were also performed. The outputs were assessed against classified Landsat TM scenes, converted to proportions maps and coregistered to the AVHRR images. An empirical method was applied for determining the critical forest cover per km² needed for an AVHRR pixel to be classified as forest. The critical values exceeded 50 per cent, indicating a tendency for AVHRR classification to underestimate forest cover. This was confirmed by comparing estimates of total forest cover obtained from the AVHRR and TM classifications. In the case of Ghana, a more accurate estimate of forest cover was obtained from the AVHRR mixture model than from the classification. Both mixture model outputs were found to be well correlated with those from Landsat TM. Further work should test the robustness of the approach adopted here when applied to much larger areas.     

Classification accuracy and trend assessments of land cover- land use changes from principal components of land satellite images

The paper evaluated the accuracy of classifying Land Cover-Land Use (LCLU) types and assessed the trends of their changes from Principal Components (PC) of Land satellite (Landsat) images. The accuracy of the image classification of LCLU was evaluated using the confusion matrices and assessed with cross-referencing of samples of LCLU types interpreted and classified from System Pour l’Observation de la Terre (SPOT) images and topographical map. LCLU changes were detected, quantified, and statistically analysed. The interpretation error of the composite image of Landsat Enhanced Thematic Mapper Plus (Landsat ETM⁺) (2006) was high compared with that from the PC image of Landsat ETM⁺ (2006). From 1986 – 2006 the area covered by settlements increased by 0.8% (230,380.00 km²), agricultural land decreased by 7.5% (1009.40 km²), vegetation cover decreased by 0.9% (114.00 km²) while waterbody increased by 0.2% (25.91 km²). Also, from 1986 – 2006 the average annual rates of change in the area of settlements was 6.7%. Agricultural land and bare land showed fluctuations of change rates from 6.7% and 5.0% annually in 1986 and 2006 respectively. The quantitative evidences of LCLU changes revealed the growth of settlements. The conversions of land from agriculture to urban land represent the most significant land cover changes. The rate of change was as high as 4.8% for settlements while agricultural lands were converted at 5.0% per year. The Principal Component Analysis (PCA) of the Landsat images and supervised classification method used made it possible to classify and determine the area of LCLU classes from the set of Landsat images without prior depiction and delimitation of individual LCLU type. It permitted the measurement of area of each LCLU class at a high accuracy level and kept the level of error relatively constant. The PCA analysis in this study affirms the previous research findings. Future research works should focus on the use of remotely sensed images with high temporal and spatial resolutions such as Quick Bird and SPOT 6 to develop effective and accurate LCLU change mapping and monitoring at the local scale.

The PCA technique has been used quite widely to study changes in land cover and land use in many ‘developed’ countries but much still needs to be done in developing and undeveloped countries where land cover and land use change is poorly mapped and knowledge of such changes is very important for planning development of the country.     

Long‐term monitoring of land cover changes based on Landsat imagery to improve hydrological modelling in West Africa

The spatial and temporal variability of land cover changes is a fundamental parameter to integrate when modelling water resources in order to reproduce the relations between rainfall and surface flow more precisely. This is particularly important in West Africa, where the land cover has been changing for more than 40 years under the combined impact of climatic effects and human activities. In this study, we evaluated the potential of Landsat imagery to monitor the vegetation cover in the upper Niger watershed (120 000 km²) using archive images from MSS, TM and ETM+ sensors covering three periods of time around 1975, 1985, and 2000. Because of the heterogeneity of the acquisition dates, the spatial and spectral resolution of the images, and the scale of analysis, we chose a simple system of classification. Pretreatments were applied to reduce variations between the images. Vegetation indices (NDVI) were then calculated and subsequently thresholded using the same land‐cover classification system. The thresholds were then optimized by automated recursive calculations of confusion matrices and control parcels. Our results revealed that although the accuracy was not perfect, it was nevertheless possible to estimate changes using an unconventional spatio‐temporal scale. The resulting changes were characterized by a moderate trend to deforestation with a corresponding increase in bare soils, soils with sparse vegetation, and shrublands. The spatial layers produced were then combined with a soil map to incorporate changes in surface conditions in the hydrological modelling of the Niger River.     

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).     

Mapping population density distribution at multiple scales in Zhejiang Province using Landsat Thematic Mapper and census data

Population density is usually calculated from the census data, but it is dynamic over time and updating population data is often challenging because it is time-consuming and costly. Another problem is that population data for public use are often too coarse, such as at the county scale in China. Previous research on population estimation mainly focused on megacities due to their importance in socio-economic conditions, but has not paid much attention to the township or village scale because of the sparse population density and less importance in economic conditions. In reality, population density in townships and villages plays an important role in land-use/cover change and environmental conditions. It is an urgent task to timely update population density at the township and cell-size scales. Therefore, this article aims to develop an approach to estimate population density at the township scale and at a cell size of 1 km by 1 km through downscaling the population density from county to township and then to cell size. We estimated population density using Landsat Thematic Mapper (TM) and census data in Zhejiang Province, China. Landsat TM images in 2010 were used to map impervious surface area (ISA) distribution using a hybrid approach, in which a decision tree classifier was used to extract ISA data and cluster analysis was used to further modify the ISA results. A population density estimation model was developed at the county scale, and this model was then transferred to the township scale. The population density was finally redistributed to cell-size scale based on the assumption that population only occupied the sites having ISA. This research indicates that most townships have residuals within ±50 persons/km² with a root mean squared error (RMSE) of 71.56 persons/km², and a relative RMSE of 27.6%. The spatial patterns of population density distribution at the 1 km² cell size are much improved compared to the township and county scales. This research indicates the importance of using the ISA for population density estimation, where ISA can be accurately extracted from remotely sensed data.     

Comparison of methods for LUCC monitoring over 50 years from aerial photographs and satellite images in a Sahelian catchment

Land use/cover change (LUCC) is a major indicator of the impact of climate change and human activity, particularly in the Sahel, where the land cover has changed greatly over the past 50 years. Aerial and satellite sensors have been taking images of the Earth's surface for several decades. These data have been widely used to monitor LUCC, but many questions remain concerning what type of pre-processing should be carried out on image resolutions and which methods are most appropriate for successfully mapping patterns and dynamics in both croplands and natural vegetation. This study considers these methodological questions. It uses multi-source imagery from 1952 to 2003 (aerial photographs, Corona, Landsat Multispectral Scanner (MSS), Landsat Thematic Mapper (TM) and Satellite Pour l'Observation de la Terre (SPOT) 5 images) and pursues two objectives: (i) to implement and compare a number of processing chains on the basis of multi-sensor data, in order (ii) to accurately track and quantify LUCC in a 100 km² Sahelian catchment over 50 years. The heterogeneity of the spatial and spectral resolution of the images led us to compare post-classification methods aimed at producing coherent diachronic maps based on a common land-cover nomenclature. Three main approaches were tested: pixel-based classification, vector grid-based on-screen interpretation and object-oriented classification. Within the automated approaches, we also examined the influence of spectral synthesis and spatial homogenization of the data through the use of composite bands (principal component analysis (PCA) and indices) and by resampling images at a common resolution. Classification accuracy was estimated by computing confusion matrices, by analysing overall change in the relative areas of land use/cover types and by studying the geographical coherence of the changes. These analyses indicate that on-screen interpretation is the most suitable approach for providing coherent, valid results from the multi-source images available over the study period. However, satisfactory classifications are obtained with the pixel-based and object-oriented approaches. The results also show significant sensitivity, depending on the method considered, to the combinations of bands used and to resampling. Lastly, the 50-year trends in LUCC point out a large increase in croplands and erosional surfaces with sparse vegetation and a drastic reduction in woody covers.     

Mapping land-use and land-cover change along Bolivia's Corredor Bioceánico with CBERS and the Landsat series: 1975–2008

This study uses a combination of Landsat series data (Multispectral Scanner or MSS, Thematic Mapper or TM and Enhanced Thematic Mapper or ETM+) to map land-use and land-cover change (LULCC) from 1975 to 2001. It extends the land change record to 2008 using Chinese–Brazil Earth Resources Satellite (CBERS)-2 and CBERS-2B data on a multi-scene level. It also establishes a methodology to correct for systematic distortion inherent in CBERS imagery without the loss of information present in Landsat 7 ETM+ imagery post-2003. Image analysis focuses on a 63 000 km² strip of land along a main highway and railroad in southeastern Bolivia named the Corredor Bioceánico. This strip of land is one of the most important agriculturally driven deforestation hotspots in Latin America. It is also located in one of the most poorly understood forest biomes in the world in terms of LULCC – Southern Hemisphere seasonally dry tropical forests – which have very high conservation values globally. Over the 33-year study period, approximately 12 000 km² of forest was lost among the three sub-regions – which is an area nearly the size of Connecticut. Evidence suggests that agriculture-driven deforestation is pushing into sensitive areas threatening globally important ecosystems such as those in the Chaco, Chiquitano and Pantanal as well as noteworthy protected areas. The results also show that imagery of CBERS-2 and CBERS-2B can help to fill the imagery gap created by Landsat ETM's Scan Line Corrector (SLC) failure in 2003. They can help to extend the land change record forward in time.     

Determining the rate of forest conversion in Mato Grosso,Brazil, using Landsat MSS and AVHRR data

Abstract

AVHRR-LAC thermal data and Landsat MSS and TM spectral data were used to estimate the rate of forest clearing in Mato Grosso, Brazil, between 1981 and 1984. The Brazilian state was stratified into forest and non-forest. A list sampling procedure was used in the forest stratum to select Landsat MSS scenes for processing based on estimates of fire activity in the scenes. Fire activity in 1984 was estimated using AVHRR-LAC thermal data. Slate-wide estimates of forest conversion indicate that between 1981 and 1984, 353966 ha ±77 000 ha (0·4 percent of the state area) were converted per year. No evidence of reforestation was found in this digital sample. The relationship between forest clearing rate (based on MSS-TM analysis)and fire activity (estimated using AVHRR data)was noisy (R²= 0·41). The results suggest that AVHRR data may be put to better use as a stratification tool rather than as a subsidiary variable in list sampling.     

Application of multitemporal Landsat data to monitor land cover changes in the Eastern Nile Delta region,Egypt

Due to the progressive increase in development of desert land in Egypt, the demand for efficient and accurate land cover change information is increasing. In this study, we apply the methodology of post‐classification change detection to map and monitor land cover change patterns related to agricultural development and urban expansion in the desert fringes of the Eastern Nile Delta region. Using a hybrid classification approach, we employ multitemporal Landsat Thematic Mapper (TM)/Enhanced Thematic Mapper Plus (ETM+) images from 1984, 1990 and 2003 to produce three land cover/land‐use maps. Post‐classification comparison of these maps was used to obtain ‘from–to’ statistics and change detection maps. The change detection results show that agricultural development increased by 14% through the study period. The average annual rate of land reclamation during 1990–2003 (4511 ha a^?1) was comparable to that during 1984–1990 (4644 ha a^?1), reflecting a systematic national plan for desert reclamation that went into effect. We find that the increase in urbanization (by ca 21 300 ha) during 1990–2003 was predominantly due to encroachment into traditionally cultivated land at the fringes of urban centres. Our results accurately quantify the land cover changes and delineate their spatial patterns, demonstrating the utility of Landsat data in analysing landscape dynamics over time. Such information is critical for making efficient and sustainable policies for resource management.     

Assessing the potential of VEGETATION sensor data for mapping snow and ice cover: A Normalized Difference Snow and Ice Index

The VEGETATION (VGT) sensor in SPOT 4 has four spectral bands that are equivalent to Landsat Thematic Mapper (TM) bands (blue, red, near-infrared and mid-infrared spectral bands) and provides daily images of the global land surface at a 1-km spatial resolution. We propose a new index for identifying and mapping of snow/ice cover, namely the Normalized Difference Snow/Ice Index (NDSII), which uses reflectance values of red and mid-infrared spectral bands of Landsat TM and VGT. For Landsat TM data, NDSII is calculated as NDSII_TM=(TM3-TM5)/(TM3+TM5); for VGT data, NDSII is calculated as NDSII_VGT=(B2-MIR)/(B2+MIR). As a case study we used a Landsat TM image that covers the eastern part of the Qilian mountain range in the Qinghai-Xizang (Tibetan) plateau of China. NDSIITM gave similar estimates of the area and spatial distribution of snow/ice cover to the Normalized Difference Snow Index (NDSI=(TM2-TM5)/(TM2+TM5)) which has been proposed by Hall et al. The results indicated that the VGT sensor might have the potential for operational monitoring and mapping of snow/ice cover from regional to global scales, when using NDSII_VGT.     

Estimating sub-pixel components of a semi-arid woodland

Mixture models were applied to Landsat-MSS and -TM data in a semi-arid woodland in central western New South Wales, Australia to extract information on soil, herbage and tree cover. There was a significant correlation (r² = 0-71) between estimated and ground data for tree cover using the TM data, with a mean error of ± 143 per cent, and a mean error of ±11-2 per cent for the bare soil cover estimates. Only general trends were observed using the MSS data. The main areas of confusion were between senesced herbage and soil, and between green grass and the green leaves of trees.     

Landsat mapping of annual inundation (1979–2006) of the Macquarie Marshes in semi-arid Australia

Measuring inundation over long timeframes is essential for understanding the responses of large floodplain wetlands on regulated rivers, such as the internationally Ramsar listed Macquarie Marshes (2000 km²) in central-eastern Australia. We used near-spring Landsat images (Multispectral Scanner (MSS) and Thematic Mapper (TM) imagery) over 28 years (1979–2006) and classified for inundation, integrating water and vegetation response using Iterative Self-Organizing Data Analysis (ISODATA) clustering. A spatially explicit inundation index showed that zones inundated with high frequency were mostly in the northern region. Change detection of inundation indices over three consecutive water management periods (period 1 (1979–1987), period 2 (1988–1996) and period 3 (1997–2006)) showed that zones inundated with high frequency across the Macquarie Marshes contracted, equating to the loss of three or more spring floods from each 9-year period, despite no corresponding change in annual catchment or local rainfall. Landsat represents the only effective available long-term information for analysing long-term changes in inundation patterns of floodplain wetlands.     

Lake-area mapping in the Tibetan Plateau: an evaluation of data and methods

Lake area derived from remote-sensing data is a primary data source, because changes in lake number and area are sensitive indicators of climate change. These indicators are especially useful when the climate change is not convoluted with a signal from direct anthropogenic activities. The data used for lake-area mapping is important, to avoid introducing unnecessary uncertainty into long-term trends of lake-area estimates. The methods for identifying waterbodies from satellite data are closely linked to the quality and efficiency of surface-water differentiation. However, few studies have comprehensively considered the factors affecting the selection of data and methods for mapping lake area in the Tibetan Plateau (TP), nor of evaluating their consequences. This study tests the dominant data sets (Landsat and Moderate Resolution Imaging Spectroradiometer (MODIS) data) and the methods for automated waterbody mapping on 14 large lakes (>500 km²) distributed across different climate zones of the TP. Seasonal changes in lake area and data availability from Landsat imagery are evaluated. Data obtained in October is optimal because in this month the lake area is relatively stable. The data window can be extended to September and November if insufficient data is available in October. Grouping data into three-year bins decreases the effects of year-to-year seasonal variability and provides a long-term trend that is suitable for time series analysis. The Landsat data (Multispectral Scanner, MSS; Thematic Mapper, TM; Enhanced Thematic Mapper Plus, ETM+; and Operational Land Imager, OLI) and MODIS data (MOD09A1) showed good performance for lake-area mapping. The Otsu method is used to determine the optimal threshold for distinguishing water from non-water features. Several water extraction indices, namely NDWI_McFeeters, NDWI_Xu, and AWEI_non-shadow, yielded high overall classification accuracy (92%), kappa coefficient (0.83), and user’s accuracy (~90%) for lake-water classification using Landsat data. The MODIS data using NDWI_McFeeters and NDWI_Xu showed consistent lake area (r² = 0.99) compared with Landsat data on the corresponding date with root mean square error (RMSE) values of 86.87 and 103.33 km² and mean absolute error (MAE) values of 25.7 and 29.04 km², respectively. The MODIS data is suitable for great lake mapping, which is the case for the large lakes in the TP. Although automated water extraction indices exhibited high accuracy in separating water from non-water, visual examination and manual editing are still necessary. Combined with recent Chinese high-resolution satellites, these remotely sensed imageries will provide a wealth of data for studies of lake dynamics and long-term lake evolution in the TP.