Comparison of MODIS and Landsat TM5 images for mapping tempo–spatial dynamics of Secchi disk depths in Poyang Lake National Nature Reserve,China

Landsat has successfully been applied to map Secchi disk depth of inland water bodies. Operational use for monitoring a dynamic variable like Secchi disk depth is however limited by the 16‐day overpass cycle of the Landsat system and cloud cover. Low spatial resolution Moderate Resolution Imaging Spectroradiometer (MODIS) image captured twice a day could potentially overcome these problems. However, its potential for mapping Secchi disk depth of inland water bodies has so far rarely been explored. This study compared two image sources, MODIS and Landsat Thematic Mapper (TM), for mapping the tempo–spatial dynamics of Secchi disk depth in Poyang Lake National Nature Reserve, China. Secchi disk depths recorded at weekly intervals from April to October in 2004 and 2005 were related to 5 Landsat TM and 22 MODIS images respectively. Two multiple regression models including the blue and red bands of Landsat TM and MODIS respectively explained 83% and 88% of the variance of the natural logarithm of Secchi disk depth. The standard errors of the predictions were 0.20 and 0.37 m for Landsat TM and MODIS‐based models. A high correlation (r = 0.94) between the predicted Secchi disk depth derived from the two models was observed. A discussion of advantages and disadvantages of both sensors leads to the conclusion that MODIS offers the possibility to monitor water transparency more regularly and cheaply in relatively big and frequently cloud covered lakes as is with Poyang Lake.     

An initial validation of Landsat 5 and 7 derived surface water temperature for U.S. lakes,reservoirs, and estuaries

ABSTRACT

The United States Harmful Algal Bloom and Hypoxia Research Control Act of 2014 identified the need for forecasting and monitoring harmful algal blooms (HAB) in lakes, reservoirs, and estuaries across the nation. Temperature is a driver in HAB forecasting models that affects both HAB growth rates and toxin production. Therefore, temperature data derived from the U.S. Geological Survey Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper Plus thermal band products were validated across 35 lakes and reservoirs, and 24 estuaries. In situ data from the Water Quality Portal (WQP) were used for validation. The WQP serves data collected by state, federal, and tribal groups. Discrete in situ temperature data included measurements at 11,910 U.S. lakes and reservoirs from 1980 through 2015. Landsat temperature measurements could include 170,240 lakes and reservoirs once an operational product is achieved. The Landsat-derived temperature mean absolute error was 1.34°C in lake pixels >180 m from land, 4.89°C at the land-water boundary, and 1.11°C in estuaries based on comparison against discrete surface in situ measurements. This is the first study to quantify Landsat resolvable U.S. lakes and reservoirs, and large-scale validation of an operational satellite provisional temperature climate data record algorithm. Due to the high performance of open water pixels, Landsat satellite data may supplement traditional in situ sampling by providing data for most U.S. lakes, reservoirs, and estuaries over consistent seasonal intervals (even with cloud cover) for an extended period of record of more than 35 years.     

Automatic detection of coastal plumes using Landsat TM/ETM+ images

Riverine fresh water outflows create coastal plumes that are distinguished from surrounding sea water by their specific spectral signature. Coastal waters are unique ecosystems, and they are very important in terms of living resources and oceanographic processes. River plumes and coastal turbid waters have important effects on coastal marine ecosystems, and they also influence marine life cycles, sediment distribution, and pollution. Remote sensing and digital image-processing techniques provide an effective tool to detect and monitor these plume zones over large areas. The primary goal of this study was automatic detection and monitoring of coastal plume zones using multispectral Landsat Thematic Mapper/Enhanced Thematic Mapper Plus (TM/ETM+) imagery. For that purpose, the proposed algorithm exploits spectral features of the multispectral images by using feature extraction and decision-making steps. The procedure has two main stages: (1) some pre-processing operations were applied to the images in order to extract the plume core reflectance values with maximum turbidity and offshore water mass reflectance values; (2) a k-means algorithm was applied with initial seed values of reflectance computed from the pre-processing stage to classify coastal plume zones. Spatial pattern and variability of optical characteristics of coastal plume zones were then defined following the results of the classification process. The algorithm was automatically applied in three different regions with three multispectral Landsat images acquired on different dates, and yielded a very high classification accuracy in detecting coastal plume zones.     

Linear 3-D transformations of Landsat 5 TM satellite images for the enhancement of archaeological signatures during the phenological cycle of crops

A linear 3-D transformation that can be used for the enhancement of crop marks related to buried archaeological features is developed and presented in this article. The methodology is based on three steps: (1) recalculation of the ground narrowband spectroradiometric measurements to the multispectral Landsat 5 Thematic Mapper (TM) sensor, based on the relative response filter of the sensor, (2) application of principal component analysis (PCA) transformation in order to determine the initial axes used for the orthogonal transformation, followed by (3) a 3-D rotation of the PCA axes. The linear coefficients of the transformation were retrieved and adjusted to different phenological stages of the crops. The transformation was successfully evaluated using both in situ measurements and Landsat 5 TM images in two different archaeological case studies. The proposed transformation tends to enhance archaeological signatures better than other established vegetation indices or algorithms, while the methodology can be expanded to any other multispectral satellite images using only the visible and very near-infrared part of the spectrum.     

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.     

A 20-year Landsat water clarity census of Minnesota's 10,000 lakes

A 20-year comprehensive water clarity database assembled from Landsat imagery, primarily Thematic Mapper and Enhanced Thematic Mapper Plus, for Minnesota lakes larger than 8 ha in surface area contains data on more than 10,500 lakes at five-year intervals over the period 1985–2005. The reliability of the data was evaluated by examining the precision of repeated measurements on individual lakes within short time periods using data from adjacent overlapping Landsat paths and by comparing water clarity computed from Landsat data to field-collected Secchi depth data. The agreement between satellite data and field measurements of Secchi depth within Landsat paths was strong (average R² of 0.83 and range 0.71–0.96). Relationships between late-summer Landsat and field-measured Secchi depth for the combined statewide data similarly were strong (r² of 0.77–0.80 for individual time periods and r² = 0.78 for the entire database). Lake clarity has strong geographic patterns in Minnesota; lakes in the south and southwest have low clarity, and lakes in the north and northeast tend to have the highest clarity. This pattern is evident at both the individual lake and the ecoregion level. Mean water clarity in the Northern Lakes and Forest and North Central Hardwood Forest ecoregions in central and northern Minnesota remained stable from 1985 to 2005 while decreasing water clarity trends were detected in the Western Corn Belt Plains and Northern Glaciated Plains ecoregions in southern Minnesota, where agriculture is the predominant land use. Mean water clarity at the statewide level also remained stable with an average around 2.25 m from 1985 to 2005. This assessment demonstrates that satellite imagery can provide an accurate method for obtaining comprehensive spatial and temporal coverage of key water quality characteristics that can be used to detect trends at different geographic scales.     

A procedure to determine the particulate content of shallow water from Thematic Mapper data

A procedure for the determination of chlorophyll and suspended sediment concentration in shallow water from Landsat Thematic Mapper data has been developed, based on numerical computation.     

Effects of sediments and coloured dissolved organic matter on remote sensing of chlorophyll-a using Landsat TM/ETM+ over turbid waters

Remote sensing of chlorophyll-a is challenging in water containing inorganic suspended sediments (i.e. non-volatile suspended solids, NVSS) and coloured dissolved organic matter (CDOM). The effects of NVSS and CDOM on empirical remote-sensing estimates of chlorophyll-a in inland waters have not been determined on a broad spatial and temporal scale. This study evaluated these effects using a long-term (1989–2012) data set that included chlorophyll-a, NVSS, and CDOM from 39 reservoirs across Missouri (USA). Model comparisons indicated that the machine-learning algorithm BRT (boosted regression trees, validation Nash–Sutcliffe coefficient = 0.350) was better than linear regression (validation Nash–Sutcliffe coefficient = 0.214) for chlorophyll-a estimate using Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) imagery. Only a small proportion of BRT model residuals could be explained by sediments or CDOM, and the observed trends in BRT residuals were different from the theoretical effects expected from NVSS and CDOM. Our results also indicated a small systematic bias by the BRT model, but it was not likely caused by NVSS or CDOM.     

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.     

Improving 30 m global land-cover map FROM-GLC with time series MODIS and auxiliary data sets: a segmentation-based approach

FROM-GLC (Fine Resolution Observation and Monitoring of Global Land Cover) is the first 30 m resolution global land-cover map produced using Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) data. Due to the lack of temporal features as inputs in producing FROM-GLC, considerable confusion exists among land-cover types (e.g. agriculture lands, grasslands, shrublands, and bareland). The Moderate Resolution Imaging Spectrometer (MODIS) provides high-temporal frequency information on surface cover. Other auxiliary bioclimatic, digital elevation model (DEM), and world maps on soil-water conditions are possible sources for improving the accuracy of FROM-GLC. In this article, a segmentation-based approach was applied to Landsat imagery to down-scale coarser-resolution MODIS data (250 m) and other 1 km resolution auxiliary data to the segment scale based on TM data. Two classifiers (support vector machine (SVM) and random forest (RF)) and two different strategies for use of training samples (global and regional samples based on a spatial temporal selection criterion) were performed. Results show that RF based on the global use of training samples achieves an overall classification accuracy of 67.08% when assessed by test samples collected independently. This is better than the 64.89% achieved by FROM-GLC based on the same set of test samples. Accuracies for vegetation cover types are most substantially improved.     

Production of Landsat ETM+ reference imagery of burned areas within Southern African savannahs: comparison of methods and application to MODIS

Accurate production of regional burned area maps are necessary to reduce uncertainty in emission estimates from African savannah fires. Numerous methods have been developed that map burned and unburned surfaces. These methods are typically applied to coarse spatial resolution (1 km) data to produce regional estimates of the area burned, while higher spatial resolution (<30 m) data are used to assess their accuracy with little regard to the accuracy of the higher spatial resolution reference data. In this study we aimed to investigate whether Landsat Enhanced Thematic Mapper (ETM+)‐derived reference imagery can be more accurately produced using such spectrally informed methods. The efficacy of several spectral index methods to discriminate between burned and unburned surfaces over a series of spatial scales (ground, IKONOS, Landsat ETM+ and data from the MOderate Resolution Imaging Spectrometer, MODIS) were evaluated. The optimal Landsat ETM+ reference image of burned area was achieved using a charcoal fraction map derived by linear spectral unmixing (k = 1.00, a = 99.5%), where pixels were defined as burnt if the charcoal fraction per pixel exceeded 50%. Comparison of coincident Landsat ETM+ and IKONOS burned area maps of a neighbouring region in Mongu (Zambia) indicated that the charcoal fraction map method overestimated the area burned by 1.6%. This method was, however, unstable, with the optimal fixed threshold occurring at >65% at the MODIS scale, presumably because of the decrease in signal‐to‐noise ratio as compared to the Landsat scale. At the MODIS scale the Mid‐Infrared Bispectral Index (MIRBI) using a fixed threshold of >1.75 was determined to be the optimal regional burned area mapping index (slope = 0.99, r ² = 0.95, SE = 61.40, y = Landsat burned area, x = MODIS burned area). Application of MIRBI to the entire MODIS temporal series measured the burned area as 10 267 km² during the 2001 fire season. The char fraction map and the MIRBI methodologies, which both produced reasonable burned area maps within southern African savannah environments, should also be evaluated in woodland and forested environments.     

Improvements in land surface temperature retrieval based on atmospheric water vapour content and atmospheric temperature

Three methods are currently used to retrieve land surface temperatures (LSTs) from thermal infrared data supplied by the Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) sensors: the radiative transfer equation, mono-window, and generalized single-channel algorithms. Most retrieval results obtained using these three methods have an average error of more than 1 K. But if the regional mean atmospheric water vapour content and temperature are supplied by in situ radiosounding observations, the mono-window algorithm is able to provide better results, with a mean error of 0.5 K. However, there are no in situ radiosounding data for most regions. This article provides an improved method to retrieve LST from Landsat TM and ETM+ data using atmospheric water vapour content and atmospheric temperature, which can be obtained from remote-sensing data. The atmospheric water vapour content at the pixel scale was first calculated from Moderate Resolution Imaging Spectroradiometer (MODIS) data. The emissivities of various land covers and uses were then defined by Landsat TM or ETM+ data. In addition, the temperature–vegetation index method was applied to map area-wide instantaneous near-surface air temperatures. The parameters of mean atmospheric water vapour content and temperature and land surface emissivity were finally inputted to the mono-window algorithm to improve the LST retrieval precision. Our results indicate that this improved mono-window algorithm gave a significantly better retrieval of the estimated LST than that using the standard mono-window algorithm, not only in dry and elevated mountain regions but also in humid regions, as shown by the bias, standard deviation (σ), and root mean square deviation (RMSD). In Madoi County, the improved mono-window algorithm validated against the LST values measured in situ produced a bias and RMSD of –0.63 K and 0.91 K, respectively, compared with the mono-window algorithm’s bias and RMSD of –1.08 K and 1.27 K. Validated against the radiance-based method, the improved algorithm shows bias and RMSD values of –1.08 K and 1.27 K, respectively, compared with the initial algorithm’s bias and RMSD –1.65 K and 1.75 K. Additionally, the improved mono-window algorithm also appeared to be more accurate than the mono-window algorithm, with lower error values when validated against in situ measurement and the radiance-based method in the validation area in Zhangye City, Gansu Province, China. Remarkable LST accuracy improvements are shown by the improved mono-window algorithm, with better agreement not only with the in situ measurements but also with the simulated LSTs in the two validation areas, indicating the soundness and suitability of this method.     

Downscaling MODIS-derived water maps with high-precision topographic data in a shallow lake

ABSTRACT

Remotely sensed imagery is the most efficient and widely used data source to monitor the water area changes. However, a trade-off always exists between temporal resolution and spatial resolution for satellite images. Taking the southern Dongting Lake as an example, this study was conducted to develop a method of downscaling the Moderate Resolution Imaging Spectroradiometer (MODIS)-derived coarse spatial resolution water maps in shallow lakes with high-precision digital elevation model. The main principle of the method is to identify and adjust the horizontal location errors of the waterlines extracted from coarse-resolution data by analysing and modifying the elevation leaps using finer-scale topography information. Moving average filter was used to smooth the errors of waterlines caused by the geometric inaccuracies and classification uncertainties of the coarse data. The optimal local window size of the moving average filter was selected automatically using an exponential decay function model and a curvature algorithm for each pixel in the waterlines. In reference to Landsat Thematic Mapper data, the accuracy of the downscaling result is distinctly higher than that of the original MODIS normalized difference water index-derived water maps. The presented method is proved to be an effective tool for acquiring water maps of shallow lake with high spatio-temporal resolution using coarse- or moderate-resolution satellite imagery and high-precision topographic data.     

Obtaining crop-specific time profiles of NDVI: the use of unmixing approaches for serving the continuity between SPOT-VGT and PROBA-V time series

The study examined the potential of two unmixing approaches for deriving crop-specific normalized difference vegetation index (NDVI) profiles so that upon availability of Project for On-Board Autonomy – Vegetation (PROBA-V) imagery in winter 2013, this new data set can be combined with existing Satellite Pour l’Observation de la Terre – VEGETATION (SPOT-VGT) data despite the differences in spatial resolution (300 m of PROBA-V versus 1 km of SPOT-VGT). To study the problem, two data sets were analysed: (1) a set of 10 temporal NDVI images, with 300 and 1000 m spatial resolution, from the state of São Paulo (Brazil) synthesized from 30 m Landsat Thematic Mapper (TM) images, and (2) a corresponding set of 10 observed Moderate Resolution Imaging Spectroradiometer (MODIS) images (250 m spatial resolution). To mimic the influence of noise on the retrieval accuracy, different sensor/atmospheric noise levels were applied to the first data set. For the unmixing analysis, a high-resolution land-cover (LC) map was used. The LC map was derived beforehand using a different set of Landsat TM images. The map distinguishes nine classes, with four different sugarcane stages, two agricultural sub-classes, plus forest, pasture, and urban/water. Unmixing aiming at the retrieval of crop-specific NDVI profiles was done at administrative level. For the synthesized data set it was demonstrated that the ‘true’ NDVI temporal profiles of different land-cover classes (from 30 m TM data) can generally be retrieved with high accuracy. The two simulated sensors (PROBA-V and SPOT-VGT) and the two unmixing algorithms gave similar results. Analysing the MODIS data set, we also found a good correspondence between the modelled NDVI profiles (both approaches) and the (true) Landsat temporal endmembers.     

Depth‐variant spectral characteristics of submersed aquatic vegetation detected by Landsat 7 ETM+

Spectral variations along depth profiles were compared using two subsets of a Landsat 7 Enhanced Thematic Mapper (ETM+) scene to test the difference between submersed aquatic vegetation (SAV) and non‐vegetated bare substrate in their depth‐induced spectral variation. Field‐surveyed water depth and SAV cover along transects were overlaid with the satellite image of Lake Pontchartrain, LA, USA. Digital numbers on the survey transects for each band and for band ratios were correlated with depth and vegetation cover. Band 1/band 3 correlated well with depth in both SAV and bare substrates, indicating that this ratio least reflects the effect of SAV. The ratio of bands 2 and 1 correlated best with vegetation cover within the shallow estuarine waters. Correlations between depth and the ratio of band 2/band3 showed contrasting results between the two substrate types (SAV and bare), suggesting that the depth‐induced variations in the band ratio can be used as indicators of SAV.     

Drought monitoring of the wetland in the Tumen River Basin between 1991 and 2016 using Landsat TM/ETM+

Wetland areas are known as ‘the kidneys of the Earth’ because they provide important functions towards stabilizing the environment, long-term protection of water sources, effectively minimizing sediment loss, purifying surface water from industrial and agricultural pollutants, and enhancing aquifer recharge. The condition of water supply in wetlands directly affects the growth of wetland plants and local biodiversity. Therefore, drought monitoring is vital in wetlands. In this study, Vegetation Temperature Condition Index (VTCI) derived from normalized difference vegetation index (NDVI) and land surface temperature (LST) is used to observe the drought status of the wetland in the cross-border (China and North Korea) Tumen River Basin from 1991 to 2016. For this purpose, the Landsat Thematic Mapper/Enhanced Thematic Mapper Plus (TM/ETM+) data for six periods were used for the analysis. Soil moisture maps acquired from the China Meteorological Administration Land Data Assimilation System Version 1.0 (CLDAS-V1.0) were then introduced for validating the reliability of the drought monitoring method. The results showed that most areas with a normal moisture level (decreased 25.8%) began experiencing slight drought (increased 29.7%). The coefficient of determination (R²) between VTCI and soil moisture showed values of 0.69, 0.32, and 0.2 for 0–5 cm, 0–10 cm, and 10–20 cm thicknesses, respectively. Although climate change probably contributes to the formation of drought by decreasing precipitation (50 mm decrease in Chinese section) and increasing temperature (0.5°C increase in North Korean section), human activities such as surges in daily water consumption appear as the main threats that leading to droughts in this wetland.     

Deriving land surface temperature from Landsat 5 and 7 during SMEX02/SMACEX

A sequence of five high-resolution satellite-based land surface temperature (T_s) images over a watershed area in Iowa were analyzed. As a part of the SMEX02 field experiment, these land surface temperature images were extracted from Landsat 5 Thematic Mapper (TM) and Landsat 7 Enhanced Thematic Mapper (ETM) thermal bands. The radiative transfer model MODTRAN 4.1 was used with atmospheric profile data to atmospherically correct the Landsat data. NDVI derived from Landsat visible and near-infrared bands was used to estimate fractional vegetation cover, which in turn was used to estimate emissivity for Landsat thermal bands. The estimated brightness temperature was compared with concurrent tower based measurements. The mean absolute difference (MAD) between the satellite-based brightness temperature estimates and the tower based brightness temperature was 0.98 °C for Landsat 7 and 1.47 °C for Landsat 5, respectively. Based on these images, the land surface temperature spatial variation and its change with scale are addressed. The scaling properties of the surface temperature are important as they have significant implications for changes in land surface flux estimation between higher-resolution Landsat and regional to global sensors such as MODIS.     

Leaf area index mapping in northern Canada

Leaf area index (LAI) is an important structural vegetation parameter that is commonly derived from remotely sensed data. It has been used as a reliable indicator for vegetation's cover, status, health and productivity. In the past two decades, various Canada-wide LAI maps have been generated by the Canada Centre for Remote Sensing (CCRS). These products have been produced using a variety of very coarse satellite data such as those from SPOT VGT and NOAA AVHRR satellite data. However, in these LAI products, the mapping of the Canadian northern vegetation has not been performed with field LAI measurements due in large part to scarce in situ measurements over northern biomes. The coarse resolution maps have been extensively used in Canada, but finer resolution LAI maps are needed over the northern Canadian ecozones, in particular for studying caribou habitats and feeding grounds.

In this study, a new LAI algorithm was developed with particular emphasis over northern Canada using a much finer resolution of remotely sensed data and in situ measurements collected over a wide range of northern arctic vegetation. A statistical relationship was developed between the in situ LAI measurements collected over vegetation plots in northern Canada and their corresponding pixel spectral information from Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) data. Furthermore, all Landsat TM and ETM+ data have been pre-normalized to NOAA AVHRR and SPOT VGT data from the growing season of 2005 to reduce any seasonal or temporal variations. Various spectral vegetation indices developed from the Landsat TM and ETM?+?data were analysed in this study. The reduced simple ratio index (RSR) was found to be the most robust and an accurate estimator of LAI for northern arctic vegetation. An exponential relationship developed using the Theil–Sen regression technique showed an R ² of 0.51 between field LAI measurement and the RSR. The developed statistical relationship was applied to a pre-existing Landsat TM 250 m resolution mosaic for northern Canada to produce the final LAI map for northern Canada ecological zones. Furthermore, the 250 m resolution LAI estimates, per ecological zone, were almost generally lower than those of the CCRS Canada-wide VGT LAI maps for the same ecozones. Validation of the map with LAI field data from the 2008 season, not used in the derivation of the algorithm, shows strong agreement between the in situ LAI measurement values and the map-estimated LAI values.     

Mapping inland lake water quality across the Lower Peninsula of Michigan using Landsat TM imagery

The number, size, and distribution of inland freshwater lakes present a challenge for traditional water-quality assessment due to the time, cost, and logistical constraints of field sampling and laboratory analyses. To overcome this challenge, Landsat imagery has been used as an effective tool to assess basic water-quality indicators, such as Secchi depth (SD), over a large region or to map more advanced lake attributes, such as cyanobacteria, for a single waterbody. The overarching objective of this research application was to evaluate Landsat Thematic Mapper (TM) for mapping nine water-quality metrics over a large region and to identify hot spots of potential risk. The second objective was to evaluate the addition of landscape pattern metrics to test potential improvements in mapping lake attributes and to understand drivers of lake water quality in this region. Field-level in situ water-quality measurements were collected across diverse lakes (n = 42) within the Lower Peninsula of Michigan. A multicriteria statistical approach was executed to map lake water quality that considered variable importance, model complexity, and uncertainty. Overall, band ratio radiance models performed well (R² = 0.65–0.81) for mapping SD, chlorophyll-a, green biovolume, total phosphorus (TP), and total nitrogen (TN) with weaker (R² = 0.37) ability to map total suspended solids (TSS) and cyanobacteria levels. In this application, Landsat TM and pattern metrics showed poor ability to accurately map non-purgable organic carbon (NPOC) and diatom biovolume, likely due to a combination of gaps in temporal overpass and field sampling and lack of signal sensitivity within broad spectral channels of Landsat TM. The composition and configuration of croplands, urban, and wetland patches across the landscape were found to be moderate predictors of lake water quality that can complement lake remote-sensing data. Of the 4071 lakes, over 4 ha in the Lower Peninsula, approximately two-thirds, were identified as mesotrophic (n = 2715). This application highlights how an operational tool might support lake decision-making or assessment protocols to identify hot spots of potential risk.     

Assessment of SeaWiFS,MODIS, and MERIS ocean colour products in the South China Sea

The Sea-viewing Wide Field-of-view Sensor (SeaWiFS), Moderate Resolution Imaging Spectroradiometer (MODIS), and Medium Resolution Imaging Spectrometer (MERIS) remote-sensing radiometric and chlorophyll-a (chl-a) concentration products for the South China Sea (SCS) from October 2003 to May 2010 were assessed using in situ data. A strict spatiotemporal match-up method was used to minimize the temporal variability effects of atmosphere and seawater around the measurement site. A comparison of the remote-sensing reflectance (R_rs(λ)) of the three sensors with in situ values from the open waters of the SCS showed that the mean absolute percentage difference varied from 13% to 55% in the 412–560 nm spectral range. Generally, the MERIS radiometric products exhibited higher typical uncertainties and bias than the SeaWiFS and MODIS products. The R_rs(443) to R_rs(555/551/560) band ratios of the satellite data were in good agreement with in situ observations for these sensors. The SeaWiFS, MODIS, and MERIS chl-a products overestimated in situ values by 74%, 42%, and 120%, respectively. MODIS retrieval accuracy was better than those of the other sensors, with MERIS performing the worst. When the match-up criteria were relaxed, the assessment results degraded systematically. Therefore, strict spatiotemporal match-up is recommended to minimize the possible influences of small-scale variation in geophysical properties around the measurement site. Coastal and open-sea areas in the SCS should be assessed separately because their biooptical properties are different and the results suggest different atmospheric correction problems.