Landsat ETM+ images in the estimation of seasonal lake water quality in boreal river basins

Remote Analysis of the Chlorophyll-a Concentration Using Sentinel-2 MSI Images in a Semiarid Environment in Northeastern Brazil

In this paper, the authors use remote-sensing images to monitor the water quality of reservoirs located in the semiarid region of Northeast Brazil. Sentinel-2 MSI TOA Level 1C reflectance images were used to remotely estimate the concentration of chlorophyll-a (chl-a), the main indicator of the trophic state of aquatic environments, in five reservoirs in the state of Ceará, Brazil. A three-spectral band retrieval model was calibrated using 171 water samples, collected from November 2015 through July 2018 in 5 reservoirs. For validation, 71 additional samples, collected from August 2018 through December 2019, were used to ensure a robust accuracy assessment. The TOA Level 1C products performed very well, achieving a relative RMSE of 28% and r2 = 0.80. Data on wind direction and speed, solar radiation and reservoir volume were used to generate a conceptual model to analyze the behavior of chl-a in the surface waters of the Castanhão reservoir. During 2019, the reservoir water quality showed strong variation, with concentration fluctuating from 30 to 95 µg/L We showed that the end of the dry season is marked by strong eutrophic conditions corresponding to very low water inflows into the reservoir. During the rainy season there is a large decrease in the chl-a concentration following the increase of the lake water storage. During the following dry season, satellite data show a progressive improvement of the trophic state controlled by wind intensity that promotes a better mixing of the reservoir waters and inhibiting the development of most phytoplankton.

Modeling and Spatiotemporal Mapping of Water Quality through Remote Sensing Techniques: A Case Study of the Hassan Addakhil Dam

With its high water potential, the Ziz basin is one of the most important basins in Morocco. This paper aims to develop a methodology for spatiotemporal monitoring of the water quality of the Hassan Addakhil dam using remote sensing techniques combined with a modeling approach. Firstly, several models were established for the different water quality parameters (nitrate, dissolved oxygen and chlorophyll a) by combining field and satellite data. In a second step, the calibration and validation of the selected models were performed based on the following statistical parameters: compliance index R2, the root mean square error and p-value. Finally, the satellite data were used to carry out spatiotemporal monitoring of the water quality. The field results show excellent quality for most of the samples. In terms of the modeling approach, the selected models for the three parameters (nitrate, dissolved oxygen and chlorophyll a) have shown a good correlation between the measured and estimated values with compliance index values of 0.62, 0.56 and 0.58 and root mean square error values of 0.16 mg/L, 0.65 mg/L and 0.07 µg/L for nitrate, dissolved oxygen and chlorophyll a, respectively. After the calibration, the validation and the selection of the models, the spatiotemporal variation of water quality was determined thanks to the multitemporal satellite data. The results show that this approach is an effective and valid methodology for the modeling and spatiotemporal mapping of water quality in the reservoir of the Hassan Addakhil dam. It can also provide valuable support for decision-makers in water quality monitoring as it can be applied to other regions with similar conditions.

Multiple Images Improve Lake CDOM Estimation: Building Better Landsat 8 Empirical Algorithms across Southern Canada

Coloured dissolved organic matter (CDOM) is an important water property for lake management. Remote sensing using empirical algorithms has been used to estimate CDOM, with previous studies relying on coordinated field campaigns that coincided with satellite overpass. However, this requirement reduces the maximum possible sample size for model calibration. New satellites and advances in cloud computing platforms offer opportunities to revisit assumptions about methods used for empirical algorithm calibration. Here, we explore the opportunities and limits of using median values of Landsat 8 satellite images across southern Canada to estimate CDOM. We compare models created using an expansive view of satellite image availability with those emphasizing a tight timing between the date of field sampling and the date of satellite overpass. Models trained on median band values from across multiple summer seasons performed better (adjusted R2 = 0.70, N = 233) than models for which imagery was constrained to a 30-day time window (adjusted R2 = 0.45). Model fit improved rapidly when incorporating more images, producing a model at a national scale that performed comparably to others found in more limited spatial extents. This research indicated that dense satellite imagery holds new promise for understanding relationships between in situ CDOM and satellite reflectance data across large areas.

Assessment of Eutrophication and DOC Sources Tracing in the Sea Area around Dajin Island Using CASI and MODIS Images Coupled with CDOM Optical Properties

The sea area around Dajin Island in the Pearl River Estuary is the second-largest habitat in China for the Indo-Pacific humpback dolphin (Sousa Chinensis). However, the rapid economic development of this area brings potential threats to the aquatic ecology around Dajin Island. Real-time monitoring and evaluation of the ecological health of the sea area are urgent. In this study, band ratio and single-band inversion algorithms were performed to obtain Chlorophyll-a (Chl-a) and Suspended Sediment Concentration (SSC), relying on both Compact Airborne Spectrographic Imager (CASI) and Moderate resolution Imaging Spectrometer (MODIS) images. The CASI/Chl-a with high spatial resolution was adopted to assess the eutrophication status, while the dissolved organic carbon (DOC) concentration and chromophoric dissolved organic matter (CDOM) optical properties were used to derive the material composition and sources. The results suggest that the study area is under a low to medium eutrophication state with evenly distributed low Chl-a concentration. However, higher Chl-a is observed in the outer estuary with MODIS/Chl-a. The relatively high DOC concentration, especially in the north, where aquaculture is practiced, and near the estuary’s main axis, i.e., east Dajin Island, indicates that the eutrophication state might be underestimated using satellite chlorophyll alone. CDOM optical properties indicated that terrestrial materials are the DOC’s primary material sources, but the DOC derived from fishery aquaculture cannot be ignored. The low Chl-a concentration is likely due to the turbulent hydrodynamic regime caused by jet flow and reciprocating flow in this marine area. Comprehensive observation, including the assessment of different technological platforms, is suggested for the aquatic environment.

Pattern of Turbidity Change in the Middle Reaches of the Yarlung Zangbo River, Southern Tibetan Plateau, from 2007 to 2017

Turbidity is an important indicator of riverine conditions, especially in a fragile environment such as the Tibetan Plateau. Remote sensing, with the advantages of large-scale observations, has been widely applied to monitor turbidity change in lakes and rivers; however, few studies have focused on turbidity change of rivers on the Tibetan Plateau. We investigated the pattern of turbidity change in the middle reaches of the Yarlung Zangbo River, southern Tibetan Plateau, based on multispectral satellite imagery and in situ measurements. We developed empirical models from in situ measured water leaving reflectance and turbidity, and applied the best performed s-curve models on satellite imagery from Sentinel-2, Landsat 8, and Landsat 5 to derive turbidity change in 2007–2017. Our results revealed an overall decreasing spatial trend from the upper to lower streams. Seasonal variations were observed with high turbidity from July to September and low turbidity from October to May. Annual turbidity showed a temporally slightly declining trend from 2007 to 2017. The pattern of turbidity change is affected by the confluence of tributaries and the changes in precipitation and vegetation along the river. These findings provide important insights into the responses of riverine turbidity to climate and environmental changes on the Tibetan Plateau.

Spatial Variability and Detection Levels for Chlorophyll-a Estimates in High Latitude Lakes Using Landsat Imagery

Monitoring lakes in high-latitude areas can provide a better understanding of freshwater systems sensitivity and accrete knowledge on climate change impacts. Phytoplankton are sensitive to various conditions: warmer temperatures, earlier ice-melt and changing nutrient sources. While satellite imagery can monitor phytoplankton biomass using chlorophyll a (Chl) as a proxy over large areas, detection of Chl in small lakes is hindered by the low spatial resolution of conventional ocean color satellites. The short time-series of the newest generation of space-borne sensors (e.g., Sentinel-2) is a bottleneck for assessing long-term trends. Although previous studies have evaluated the use of high-resolution sensors for assessing lakes’ Chl, it is still unclear how the spatial and temporal variability of Chl concentration affect the performance of satellite estimates. We discuss the suitability of Landsat (LT) 30 m resolution imagery to assess lakes’ Chl concentrations under varying trophic conditions, across extensive high-latitude areas in Finland. We use in situ data obtained from field campaigns in 19 lakes and generate remote sensing estimates of Chl, taking advantage of the long-time span of the LT-5 and LT-7 archives, from 1984 to 2017. Our results show that linear models based on LT data can explain approximately 50% of the Chl interannual variability. However, we demonstrate that the accuracy of the estimates is dependent on the lake’s trophic state, with models performing in average twice as better in lakes with higher Chl concentration (>20 µg/L) in comparison with less eutrophic lakes. Finally, we demonstrate that linear models based on LT data can achieve high accuracy (R2 = 0.9; p-value < 0.05) in determining lakes’ mean Chl concentration, allowing the mapping of the trophic state of lakes across large regions. Given the long time-series and high spatial resolution, LT-based estimates of Chl provide a tool for assessing the impacts of environmental change.

Detecting Long Time Changes in Benthic Macroalgal Cover Using Landsat Image Archive

Coastal macroalgae worldwide provide multiple ecological functions and support vital ecosystem services. Thereby, it is important to monitor changes in the extent of benthic macroalgal cover. However, as in situ sampling is costly and time-consuming, areal estimates of macroalgal species cover are often based only on a limited number of samples. This low sampling effort likely yields very biased estimates, as macroalgal communities are often characterized by large spatial variability at multiple spatial scales. Moreover, ecological time series are often short-term, making it impossible to assess changes in algal communities over decades and relate this to different human pressures and/or climate change. The Landsat series satellites have operated for 40 years. In the current study, we tested if the Landsat sensors could be used for mapping the cover of shallow water benthic macroalgae. This study was carried out at two sites in the West Estonian Archipelago, in the northeastern Baltic Sea. Our results show that the Landsat imagery accurately reflected both spatial and temporal variability in benthic algal cover. To conclude, the current methodology can be used to improve the existing assessments of areal macroalgal cover, or to estimate the cover values, in areas and times lacking ecological observations.

Optical Water Type Guided Approach to Estimate Optical Water Quality Parameters

Currently, water monitoring programs are mainly based on in situ measurements; however, this approach is time-consuming, expensive, and may not reflect the status of the whole water body. The availability of Multispectral Imager (MSI) and Ocean and Land Colour Instrument (OLCI) free data with high spectral, spatial, and temporal resolution has increased the potential of adding remote sensing techniques into monitoring programs, leading to improvement of the quality of monitoring water. This study introduced an optical water type guided approach for boreal regions inland and coastal waters to estimate optical water quality parameters, such as the concentration of chlorophyll-a (Chl-a) and total suspended matter (TSM), the absorption coefficient of coloured dissolved organic matter at a wavelength of 442 nm (aCDOM(442)), and the Secchi disk depth, from hyperspectral, OLCI, and MSI reflectance data. This study was based on data from 51 Estonian and Finnish lakes and from the Baltic Sea coastal area, which altogether were used in 415 in situ measurement stations and covered a wide range of optical water quality parameters (Chl-a: 0.5–215.2 mg·m−3; TSM: 0.6–46.0 mg·L−1; aCDOM(442): 0.4–43.7 m−1; and Secchi disk depth: 0.2–12.2 m). For retrieving optical water quality parameters from reflectance spectra, we tested 132 empirical algorithms. The study results describe the best algorithm for each optical water type for each spectral range and for each optical water quality parameter. The correlation was high, from 0.87 up to 0.93, between the in situ measured optical water quality parameters and the parameters predicted by the optical water type guided approach.

Validation and Comparison of Water Quality Products in Baltic Lakes Using Sentinel-2 MSI and Sentinel-3 OLCI Data

Inland waters, including lakes, are one of the key points of the carbon cycle. Using remote sensing data in lake monitoring has advantages in both temporal and spatial coverage over traditional in-situ methods that are time consuming and expensive. In this study, we compared two sensors on different Copernicus satellites: Multispectral Instrument (MSI) on Sentinel-2 and Ocean and Land Color Instrument (OLCI) on Sentinel-3 to validate several processors and methods to derive water quality products with best performing atmospheric correction processor applied. For validation we used in-situ data from 49 sampling points across four different lakes, collected during 2018. Level-2 optical water quality products, such as chlorophyll-a and the total suspended matter concentrations, water transparency, and the absorption coefficient of the colored dissolved organic matter were compared against in-situ data. Along with the water quality products, the optical water types were obtained, because in lakes one-method-to-all approach is not working well due to the optical complexity of the inland waters. The dynamics of the optical water types of the two sensors were generally in agreement. In most cases, the band ratio algorithms for both sensors with optical water type guidance gave the best results. The best algorithms to obtain the Level-2 water quality products were different for MSI and OLCI. MSI always outperformed OLCI, with R² 0.84–0.97 for different water quality products. Deriving the water quality parameters with optical water type classification should be the first step in estimating the ecological status of the lakes with remote sensing.

Large-Scale Retrieval of Coloured Dissolved Organic Matter in Northern Lakes Using Sentinel-2 Data

Owing to the significant societal value of inland water resources, there is a need for cost-effective monitoring of water quality on large scales. We tested the suitability of the recently launched Sentinel-2A to monitor a key water quality parameter, coloured dissolved organic matter (CDOM), in various types of lakes in northern Sweden. Values of a(420)CDOM (CDOM absorption at 420 nm wavelength) were obtained by analyzing water samples from 46 lakes in five districts across Sweden within an area of approximately 800 km2. We evaluated the relationships between a(420)CDOM and band ratios derived from Sentinel-2A Level-1C and Level-2A products. The band ratios B2/B3 (460 nm/560 nm) and B3/B5 (560 nm/705 nm) showed poor relationships with a(420)CDOM in Level-1C and 2A data both before and after the removal of outliers. However, there was a slightly stronger power relationship between the atmospherically-corrected B3/B4 ratio and a(420)CDOM (R2 = 0.28, n = 46), and this relationship was further improved (R2 = 0.65, n = 41) by removing observations affected by light haze and cirrus clouds. This study covered a wide range of lakes in different landscape settings and demonstrates the broad applicability of a(420)CDOM retrieval algorithms based on the B3/B4 ratio derived from Sentinel-2A.

Time Delay Evaluation on the Water-Leaving Irradiance Retrieved from Empirical Models and Satellite Imagery

Temporal delays and spatial randomness between ground-based data and satellite overpass involve important deviations between the empirical model output and real data; these are factors poorly considered in the model calibration. The inorganic matter-generated turbidity in Lake Chapala (Mexico) was taken as a study case to expose the influence of such factors. Ground-based data from this study and historical records were used as references. We take advantage of the at-surface reflectance from Landsat-8, sun-glint corrections, a reduced NIR-band range, and null organic matter incidence in these wavelengths to diminish the physical phenomena-related radiometric artifacts; leaving the spatio-temporal relationships as the principal factor inducing the model uncertainty. Non-linear correlations were assessed to calibrate the best empirical model; none of them presented a strong relationship (<73%), including that based on hourly delays. This last model had the best predictability only for the summer-fall season, explaining 71% of the turbidity variation in 2016, and 59% in 2017, with RMSEs < 24%. The instantaneous turbidity maps depicted the hydrodynamic complexity of the lake, highlighting a strong component of spatial randomness associated with the temporal delays. Reasonably, robust empirical models will be developed if several dates and sampling-sites are synchronized with more satellite overpasses.

Monitoring the Water Quality of Small Water Bodies Using High-Resolution Remote Sensing Data

Remotely sensed data can reinforce the abilities of water resources researchers and decision-makers to monitor water quality more effectively. In the past few decades, remote sensing techniques have been widely used to measure qualitative water quality parameters. However, the use of moderate resolution sensors may not meet the requirements for monitoring small water bodies. Water quality in a small dam was assessed using high-resolution satellite data from RapidEye and in situ measurements collected a few days apart. The satellite carries a five-band multispectral optical imager with a ground sampling distance of 5 m at its nadir and a swath width of 80 km. Several different algorithms were evaluated using Pearson correlation coefficients for electrical conductivity (EC), total dissolved soils (TDS), water transparency, water turbidity, depth, suspended particular matter (SPM), and chlorophyll-a. The results indicate strong correlation between the investigated parameters and RapidEye reflectance, especially in the red and red-edge portion with highest correlation between red-edge band and water turbidity (r2 = 0.92). Two of the investigated indices showed good correlation in almost all of the water quality parameters with correlation higher than 0.80. The findings of this study emphasize the use of both high-resolution remote sensing imagery and red-edge portion of the electromagnetic spectrum for monitoring several water quality parameters in small water areas.

Toward Long-Term Aquatic Science Products from Heritage Landsat Missions

This paper aims at generating a long-term consistent record of Landsat-derived remote sensing reflectance (Rrs) products, which are central for producing downstream aquatic science products (e.g., concentrations of total suspended solids). The products are derived from Landsat-5 and Landsat-7 observations leading to Landsat-8 era to enable retrospective analyses of inland and nearshore coastal waters. In doing so, the data processing was built into the SeaWiFS Data Analysis System (SeaDAS) followed by vicariously calibrating Landsat-7 and -5 data using reference in situ measurements and near-concurrent ocean color products, respectively. The derived Rrs products are then validated using (a) matchups using the Aerosol Robotic Network (AERONET) data measured by in situ radiometers, i.e., AERONET-OC, and (b) ocean color products at select sites in North America. Following the vicarious calibration adjustments, it is found that the overall biases in Rrs products are significantly reduced. The root-mean-square errors (RMSE), however, indicate noticeable uncertainties due to random and systematic noise. Long-term (since 1984) seasonal Rrs composites over 12 coastal and inland systems are further evaluated to explore the utility of Landsat archive processed via SeaDAS. With all the qualitative and quantitative assessments, it is concluded that with careful algorithm developments, it is possible to discern natural variability in historic water quality conditions using heritage Landsat missions. This requires the changes in Rrs exceed maximum expected uncertainties, i.e., 0.0015 [1/sr], estimated from mean RMSEs associated with the matchups and intercomparison analyses. It is also anticipated that Landsat-5 products will be less susceptible to uncertainties in turbid waters with Rrs(660) > 0.004 [1/sr], which is equivalent of ~1.2% reflectance. Overall, end-users may utilize heritage Rrs products with “fitness-for-purpose” concept in mind, i.e., products could be valuable for one application but may not be viable for another. Further research should be dedicated to enhancing atmospheric correction to account for non-negligible near-infrared reflectance in CDOM-rich and extremely turbid waters.

Hybrid forward-selection method-based water-quality estimation via combining Landsat TM, ETM+, and OLI/TIRS images and ancillary environmental data

A simple approach to enable water-management agencies employing free data to create a single set of water quality predictive equations with satisfactory accuracy is proposed. Multiple regression-derived equations based on surface reflectance, band ratios, and environmental factors as predictor variables for concentrations of Total Suspended Solids (TSS) and Total Nitrogen (TN) were derived using a hybrid forward-selection method that considers both p-value and Variance Inflation Factor (VIF) in the forward-selection process. Landsat TM, ETM+, and OLI/TIRS images were jointly utilized with environmental factors, such as wind speed and water surface temperature, to derive the single set of equations. Through splitting data into calibration and validation groups, the coefficients of determination are 0.73 for TSS calibration and 0.70 for TSS validation, respectively. The coefficients of determination for TN calibration and validation are 0.64 and 0.37, respectively. Among all chosen predictor variables, ratio of reflectance of visible red (Band 3 for Landsat TM and ETM+, or Band 4 for Landsat OLI/TIRS) to visible blue (Band 1 for Landsat TM and ETM+, or Band 2 for Landsat OLI/TIRS) has a strong influence on the predictive power for TSS retrieval. Environmental factors including wind speed, remote sensing-derived water surface temperature, and time difference (in days) between the image acquisition and water sampling were found to be important in water-quality quantity estimation. The hybrid forward-selection method consistently yielded higher validation accuracy than that of the conventional forward-selection approach.

Satellite-Derived Bathymetry for Improving Canadian Hydrographic Service Charts

Approximately 1000 Canadian Hydrographic Service (CHS) charts cover Canada’s oceans and navigable waters. Many charts use information collected with techniques that predate the more advanced technologies available to Hydrographic Offices (HOs) today. Furthermore, gaps in survey data, particularly in the Canadian Arctic where only 6% of waters are surveyed to modern standards, are also problematic. Through a Canadian Space Agency (CSA) Government Related Initiatives Program (GRIP) project, CHS is exploring remote sensing techniques to assist with the improvement of Canadian navigational charts. Projects exploring optical/Synthetic Aperture Radar (SAR) shoreline extraction and change detection, as well as optical Satellite-Derived Bathymetry (SDB), are currently underway. This paper focuses on SDB extracted from high-resolution optical imagery, highlighting current results as well as the challenges and opportunities CHS will encounter when implementing SDB within its operational chart production process. SDB is of particular interest to CHS due to its ability to supplement depths derived from traditional hydrographic surveys. This is of great importance in shallow and/or remote Canadian waters where achieving wide-area depth coverage through traditional surveys is costly, time-consuming and a safety risk to survey operators. With an accuracy of around 1 m, SDB could be used by CHS to fill gaps in survey data and to provide valuable information in dynamic areas.

Remote Sensing for Irrigation of Horticultural Crops

This paper reviews the literature on applications of remote sensing for monitoring soil- and crop- water status for irrigation purposes. The review is organized into two main sections: (1) sensors and platforms applied to irrigation studies and (2) remote sensing approaches for precision irrigation to estimate crop water status, evapotranspiration, infrared thermography, soil and crop characteristics methods. Recent literature reports several remote sensing (RS) approaches to monitor crop water status in the cultivated environment. Establishing the right amount of water to supply for different irrigation strategies (maximization of yield or water use efficiency (WUE)) for a large number of crops is a problem that remains unresolved. For each crop, it will be necessary to create a stronger connection between crop-water status and crop yield.

Impact of the spatial resolution of satellite remote sensing sensors in the quantification of total suspended sediment concentration: A case study in turbid waters of Northern Western Australia

The impact of anthropogenic activities on coastal waters is a cause of concern because such activities add to the total suspended sediment (TSS) budget of the coastal waters, which have negative impacts on the coastal ecosystem. Satellite remote sensing provides a powerful tool in monitoring TSS concentration at high spatiotemporal resolution, but coastal managers should be mindful that the satellite-derived TSS concentrations are dependent on the satellite sensor’s radiometric properties, atmospheric correction approaches, the spatial resolution and the limitations of specific TSS algorithms. In this study, we investigated the impact of different spatial resolutions of satellite sensor on the quantification of TSS concentration in coastal waters of northern Western Australia. We quantified the TSS product derived from MODerate resolution Imaging Spectroradiometer (MODIS)-Aqua, Landsat-8 Operational Land Image (OLI), and WorldView-2 (WV2) at native spatial resolutions of 250 m, 30 m and 2 m respectively and coarser spatial resolution (resampled up to 5 km) to quantify the impact of spatial resolution on the derived TSS product in different turbidity conditions. The results from the study show that in the waters of high turbidity and high spatial variability, the high spatial resolution WV2 sensor reported TSS concentration as high as 160 mg L^-1 while the low spatial resolution MODIS-Aqua reported a maximum TSS concentration of 23.6 mg L^-1. Degrading the spatial resolution of each satellite sensor for highly spatially variable turbid waters led to variability in the TSS concentrations of 114.46%, 304.68% and 38.2% for WV2, Landsat-8 OLI and MODIS-Aqua respectively. The implications of this work are particularly relevant in the situation of compliance monitoring where operations may be required to restrict TSS concentrations to a pre-defined limit.

Assessment of nutrient distributions in Lake Champlain using satellite remote sensing

The introduction of nutrients to lakes causing eutrophic conditions is a major problem around the world. Proper monitoring and modeling are important to effectively manage eutrophication in lake waters. The goal is to develop remote sensing models for nutrients, total phosphorus and total nitrogen, in Lake Champlain. The remote sensing models were created using multivariate linear regression with the unique band combinations of Landsat Enhanced Thematic Mapper Plus (ETM+) imagery based on the empirical relationship with the field observations. The resulting models successfully showed nutrient distributions in the most eutrophic part of Lake Champlain, Missisquoi Bay, with reasonable adjusted coefficient of determination values (R(2)=0.81 and 0.75 for total phosphorus and total nitrogen, respectively). The results show the feasibility and the utility of satellite imagery to detect spatial distributions of lake water quality constituents, which can be used to better understand nutrient distributions in Lake Champlain. This approach can be applicable to other lakes experiencing eutrophication assisting decision making when implementing Best Management Practices and other mitigation techniques to lakes.

Hindcasting water clarity from Landsat satellite images of unmonitored shallow lakes in the Waikato region, New Zealand

Cost-effective monitoring is necessary for all investigations of lake ecosystem responses to perturbations and long-term change. Satellite imagery offers the opportunity to extend low-cost monitoring and to examine spatial and temporal variability in water clarity data. We have developed automated procedures using Landsat imagery to estimate total suspended sediments (TSS), turbidity (TURB) in nephlometric turbidity units (NTU) and Secchi disc transparency (SDT) in 34 shallow lakes in the Waikato region, New Zealand, over a 10-year time span. Fifty-three Landsat 7 Enhanced Thematic Mapper Plus images captured between January 2000 and March 2009 were used for the analysis, six of which were captured within 24 h of physical in situ measurements for each of 10 shallow lakes. This gave 32-36 usable data points for the regressions between surface reflectance signatures and in situ measurements, which yielded r (2) values ranging from 0.67 to 0.94 for the three water clarity variables. Using these regressions, a series of Arc Macro Language scripts were developed to automate image preparation and water clarity analysis. Minimum and maximum in situ measurements corresponding to the six images were 2 and 344 mg/L for TSS, 75 and 275 NTU for TURB, and 0.05 and 3.04 m for SDT. Remotely sensed water clarity estimates showed good agreement with temporal patterns and trends in monitored lakes and we have extended water clarity datasets to previously unmonitored lakes. High spatial variability of TSS and water clarity within some lakes was apparent, highlighting the importance of localised inputs and processes affecting lake clarity. Moreover, remote sensing can give a whole lake view of water quality, which is very difficult to achieve by in situ point measurements.