Remote sensing, numerical modelling and ground truthing for analysis of lake water quality and temperature

Freshwater accounts for just 2.5% of the earth’s water resources, and its quality and availability are becoming an issue of global concern in the 21st century. Growing human population, over-exploitation of water sources and pressures of global warming mean that both water quantity and quality are affected. In order to effectively manage water quality there is a need for increased monitoring and predictive modelling of freshwater resources. To address these concerns in New Zealand inland waters, an approach which integrates biological and physical sciences is needed. Remote sensing has the potential to allow this integration and vastly increase the temporal and spatial resolution of current monitoring techniques, which typically involve collecting grab-samples. In a complementary way, lake modelling has the potential to enable more effective management of water resources by testing the effectiveness of a range of possible management scenarios prior to implementation. Together, the combination of remote sensing and modelling data allows for improved model initialisation, calibration and validation, which ultimately aid in understanding of complex lake ecosystem processes. This study investigated the use of remote sensing using empirical and semi-analytical algorithms for the retrieval of chlorophyll a (chl a), tripton, suspended minerals (SM), total suspended sediment (SS) and water surface temperature. It demonstrated the use of spatially resolved statistical techniques for comparing satellite estimated and 3-D simulated water quality and temperature. An automated procedure was developed for retrieval of chl a from Landsat Enhanced Thematic Mapper (ETM+) imagery, using 106 satellite images captured from 1999 to 2011. Radiative transfer-based atmospheric correction was applied to images using the Second Simulation of the Satellite in the Solar Spectrum model (6sv). For the estimation of chl a over a time series of images, the use of symbolic regression resulted in a significant improvement in the precision of chl a hindcasts compared with traditional regression equations. Results from this investigation suggest that remote sensing provides a valuable tool to assess temporal and spatial distributions of chl a. Bio-optical models were applied to quantify the physical processes responsible for the relationship between chl a concentrations and subsurface irradiance reflectance used in regression algorithms, allowing the identification of possible sources of error in chl a estimation. While the symbolic regression model was more accurate than traditional empirical models, it was still susceptible to errors in optically complex waters such as Lake Rotorua, due to the effect of variations of SS and CDOM on reflectance. Atmospheric correction of Landsat 7 ETM+ thermal data was carried out for the purpose of retrieval of lake water surface temperature in Rotorua lakes, and Lake Taupo, North Island, New Zealand. Atmospheric correction was repeated using four sources of atmospheric profile data as input to a radiative transfer model, MODerate resolution atmospheric TRANsmission (MODTRAN) v.3.7. The retrieved water temperatures from 14 images between 2007 and 2009 were validated using a high-frequency temperature sensor deployed from a mid-lake monitoring buoy at the water surface of Lake Rotorua. The most accurate temperature estimation for Lake Rotorua was with radiosonde data as an input into MODTRAN, followed by Moderate Resolution Imaging Spectroradiometer (MODIS) Level 2, Atmospheric Infrared Sounder (AIRS) Level 3, and NASA data. Retrieved surface water temperature was used for assessing spatial heterogeneity of surface water temperature simulated with a three-dimensional (3-D) hydrodynamic model (ELCOM) of Lake Rotoehu, located approximately 20 km east of Lake Rotorua. This comparison demonstrated that simulations reproduced the dominant horizontal variations in surface water temperature in the lake. The transport and mixing of a geothermal inflow and basin-scale circulation patterns were inferred from thermal distributions from satellite and model estimations of surface water temperature and a spatially resolved statistical evaluation was used to validate simulations. This study has demonstrated the potential of accurate satellite-based thermal monitoring to validate water surface temperature simulated by 3-D hydrodynamic models. Semi-analytical and empirical algorithms were derived to determine spatial and temporal variations in SS in Lake Ellesmere, South Island, New Zealand, using MODIS band 1. The semi-analytical model and empirical model had a similar level of precision in SS estimation, however, the semi-analytical model has the advantage of being applicable to different satellite sensors, spatial locations, and SS concentration ranges. The estimations of SS concentration (and estimated SM concentration) from the semi-analytical model were used for a spatially resolved validation of simulations of SM derived from ELCOM-CAEDYM. Visual comparisons were compared with spatially-resolved statistical techniques. The spatial statistics derived from the Map Comparison Kit allowed a non-subjective and quantitative method to rank simulation performance on different dates. The visual and statistical comparison between satellite estimated and model simulated SM showed that the model did not perform well in reproducing both basin-scale and fine-scale spatial variation in SM derived from MODIS satellite imagery. Application of the semi-analytical model to estimate SS over the lifetime of the MODIS sensor will greatly extend its spatial and temporal coverage for historical monitoring purposes, and provide a tool to validate SM simulated by 1-D and 3-D models on a daily basis. A bio-optical model was developed to derive chl a, SS concentrations, and coloured dissolved organic matter /detritus absorption at 443 nm, from MODIS Aqua subsurface remote sensing reflectance of Lake Taupo, a large, deep, oligotrophic lake in North Island, New Zealand. The model was optimised using in situ inherent optical properties (IOPs) from the literature. Images were atmospherically corrected using the radiative transfer model 6sv. Application of the bio-optical model using a single chl a-specific absorption spectrum (a*ϕ(λ)) resulted in low correlation between estimated and observed values. Therefore, two different absorption curves were used, based on the seasonal dominance of phytoplankton phyla with differing absorption properties. The application of this model resulted in reasonable agreement between modelled and in situ chl a concentrations. Highest concentrations were observed during winter when Bacillariophytes (diatoms) dominated the phytoplankton assemblage. On 4 and 5 March 2004 an unusually large turbidity current was observed originating from the Tongariro River inflow in the south-east of the lake. In order to resolve fine details of the plume, empirical relationships were developed between MODIS band 1 reflectance (250 m resolution) and SS estimated from MODIS bio-optical features (1 km resolution) were used estimate SS at 250 m resolution. Complex lake circulation patterns were observed including a large clockwise gyre. With the development of this bio-optical model MODIS can potentially be used to remotely sense water quality in near real time, and the relationship developed for B1 SS allows for resolution of fine-scale features such turbidity currents

The Heterogeneity of the Abiotic and Biotic Components of Techno-Ecosystems: View from Space and from The Earth

The article discusses the use of different methods to assess the spatial structure of the biotic and abiotic component of the techno-ecosystems of cooling ponds of thermal and nuclear power plants in Belarus and Ukraine. Field studies provide point information on the parameters of biotic and abiotic variables, whereas statistical mapping and remote sensing methods makes it possible to generalize visualization and provides a holistic picture of the spatial heterogeneity of conditions and distribution of biota. Spatial changes in cooling ponds are mainly determined by hydrodynamic conditions in the pond and the presence of anthropogenic flow of cooling water. Studies have shown that the distribution of quantitative indicators of plankton in the surface layer in the Khmelnitsky NPP cooling pond was mostly determined by the position of individualstations and sections of the pond surface relative to manmade technical and wind circulation flows. In contrary, thermal conditions had the greatest influence on the characteristics of the plankton distribution in the Lukomskoye Lake - cooling pond of Lukomlskaya TPP

Comparison of methods to estimate lake-surface-water temperature using Landsat 7 ETM+ and MODIS imagery : case study of a large shallow subtropical lake in southern Brazil

Water temperature regulates many processes in lakes; therefore, evaluating it is essential to understand its ecological status and functioning, and to comprehend the impact of climate change. Although few studies assessed the accuracy of individual sensors in estimating lake-surface-water temperature (LSWT), comparative analysis considering different sensors is still needed. This study evaluated the performance of two thermal sensors, MODIS and Landsat 7 ETM+, and used Landsat methods to estimate the SWT of a large subtropical lake. MODIS products MOD11 LST and MOD28 SST were used for comparison. For the Landsat images, the radiative transfer equation (RTE), using NASA’s Atmospheric Correction Parameter Calculator (AtmCorr) parameters, was compared with the single-channel algorithm in different approaches. Our results showed that MOD11 obtained the highest accuracy (RMSE of 1.05 C), and is the recommended product for LSWT studies. For Landsat-derived SWT, AtmCorr obtained the highest accuracy (RMSE of 1.07 C) and is the recommended method for small lakes. Sensitivity analysis showed that Landsat-derived LSWT using the RTE is very sensitive to atmospheric parameters and emissivity. A discussion of the main error sources was conducted. We recommend that similar tests be applied for Landsat imagery on different lakes, further studies on algorithms to correct the cool-skin effect in inland waters, and tests of different emissivity values to verify if it can compensate for this effect, in an effort to improve the accuracy of these estimates

Analyse der Wasserfarbe von Seen mithilfe räumlich hoch und mittel auflösender Satelliten

Remote sensing techniques can assist traditional lake monitoring approaches by supplying spatial information on optically active lake ecology indicators, i.e. chlorophyll-a (CHL), total suspended matter (TSM), coloured dissolved organic matter (CDOM), and, especially in optically shallow waters, water depth and substrate composition. The present thesis provides an overview on the current research status concerning lake remote sensing and the benefit of time series analyses for lake ecology. To investigate the suitability of Sentinel-2 and Landsat 8 for lake monitoring and their combination with other sensors this thesis focused on two study areas with highly different optical characteristics, i.e. the oligotrophic Lake Starnberg (southern Germany) and the mesotrophic-eutrophic Lake Kummerow (northern Germany). Using the bio-optical model WASI-2D, Sentinel-2A turned out to be suited for retrieving low TSM and CDOM values. The high spatial resolution enabled the differentiation between bare ground and areas covered by submerged aquatic vegetation. Water depth estimations performed well until half Secchi disk depth. Cross-sensor comparisons demonstrated high correlation of CHL among timely acquired, spatially high and medium resolved sensors. Evaluations with in situ data showed that most of the sensor-in situ match-ups were within an uncertainty range of in situ measurements. Analysing a 9-year MERIS time series with FUB/WeW revealed unprecedented information on temporal trends and seasonal behaviour of CHL, TSM and CDOM at the study area Lake Kummerow. Combining CHL, retrieved with the Modular Inversion and Processing System, from different satellite sensors (MODIS, Landsat 7/ 8, Sentinel-2A) enabled detailed observations of phytoplankton development. Such combinations are a step forward to future lake analyses which may integrate remote sensing data, in situ measurements and environmental modelling.Fernerkundungstechniken können das Seemonitoring mit räumlichen Informationen über optisch aktive Indikatoren der Gewässerökologie liefern, z.B. Chlorophyll-a (CHL), suspendierte Schwebstoffe (TSM), Gelbstoffe (CDOM) und insbesondere in optisch flachen Gewässern, Wassertiefe und Substratbedeckung. Die vorliegende Arbeit gibt einen Überblick über den aktuellen Forschungsstand zur Seefernerkundung und den Nutzen von Zeitreihenanalysen für die Seeökologie. Um die Eignung von Sentinel-2 und Landsat 8 für ein Seenmonitoring und deren Kombination mit anderen Sensoren zu untersuchen, konzentrierte sich diese Arbeit auf zwei Untersuchungsgebiete mit sehr unterschiedlichen optischen Eigenschaften: den oligotrophen Starnberger See (Süddeutschland) und den mesotroph-eutrophen Kummerower See (Norddeutschland). Mit dem bio-optischen Modell WASI-2D erwies sich Sentinel-2A als geeignet, um niedrige TSM- und CDOM-Werte zu bestimmen. Die hohe räumliche Auflösung ermöglichte eine Unterscheidung zwischen unbewachsenem und mit Makrophyten bewachsenem Untergrund. Die Wassertiefenbestimmung verlief bis zur halben Sichttiefe gut. Sensorübergreifende Vergleiche zeigten eine hohe Korrelation von CHL zwischen zeitnah erfassten, räumlich mittel und hoch aufgelösten Sensoren. Auswertungen mit in-situ-Daten zeigten, dass die meisten Sensor-in-situ-Match-ups innerhalb eines Unsicherheitsbereichs von in-situ-Messungen lagen. Die Analyse einer 9-jährigen MERIS-Zeitreihe mit FUB/WeW ergab neue Informationen über zeitliche Trends und saisonales Verhalten von CHL, TSM und CDOM im Untersuchungsgebiet Kummerow See. Die Kombination von CHL aus verschiedenen Satellitensensoren (MODIS, Landsat 7/ 8, Sentinel-2A) mit dem Modular Inversion and Processing System ermöglichte detaillierte Beobachtungen der Phytoplanktonentwicklung. Solche Kombinationen sind ein Schritt für zukünftigen Gewässeranalysen, die Fernerkundungsdaten, in-situ-Messungen und Umweltmodellierung integrieren sollten

Biogeochemical Characterisation of an Alum Dosed Stream: Implications for Phosphate Cycling in Lake Rotoehu

Lakes are highly valued freshwater ecosystems which provide many goods and services upon which humans are reliant. Eutrophication of waterways, driven by the over-enrichment of nutrients such as phosphorus and nitrogen, is a threat to the future quality of water resources globally. Lake restoration methods are increasingly being employed to improve eutrophic waterways, via numerous catchment derived and in-lake approaches. However, spatio-temporal heterogeneity in physicochemical and biogeochemical conditions within lakes may restrict the efficacy of restoration methods. This thesis explores in-stream alum dosing as a lake remediation strategy for the purpose of dissolved phosphorus inactivation, and its physicochemical and biogeochemical interaction and fate within Lake Rotoehu, New Zealand. Laboratory and field measurements were undertaken at Lake Rotoehu, on New Zealand’s North Island. The shallow, polymictic lake resides in an active volcanic area with geothermal inputs, and is subject to several management issues including elevated nutrient concentrations, invasive macrophytes, and frequent cyanobacteria blooms. The geothermal Waitangi Springs, which discharge into Lake Rotoehu, contribute ~69 % of the lake’s total ionic content, and are responsible for enhanced concentrations of biologically limiting nutrients including phosphorus, silicon, nitrate, ammonium, and iron. In an effort to curb in-lake phosphorus levels, in-stream alum-dosing has been employed in order to floc out dissolved reactive phosphorus (DRP), through chemical adsorption and sedimentation. The results presented here, culminate from an investigation of the physicochemical and geochemical dynamics across the mixing zone from the Waitangi Springs geothermal stream outlet across Te Wairoa Bay to the main lake body. A combination of approaches was used: two field experiments with fixed location and transect measurements, laboratory analysis and geochemical speciation modelling with PHREEQC. Results show sharp changes in physicochemical water properties across the mixing zone within the bay: pH, O2 and dissolved reactive phosphorus values increased with distance from the stream outlet, whereas major ion concentrations, temperature and conductivity values decreased. Initial in-stream phosphorus stripping through alum dosing is effective in reducing the DRP load by ~50 % of background concentration. As alum is introduced to the stream water (~pH 6) it precipitates to form amorphous aluminium hydroxide (Al(OH)3(am)) and adsorbs phosphorus via hydroxyl ligand exchange. However, elevated levels of iron in amorphous hydrous ferric hydroxides Fe(OH)3(am) are also likely to be contributing to natural phosphorus binding capacity. Sediment core data also indicated that settled Al(OH)3(am) floc and Fe(OH)3(am) particulates were primarily concentrated within the inner portion of Te Wairoa Bay near the Waitangi Springs outlet. Surface water physicochemical and geochemical concentrations were spatially resolved and indicated distinct mixing boundaries (pH, DO, temp) and patchiness (Al, Fe, DRP) within Te Wairoa Bay associated with the locations and configuration of dense submerged macrophytes (Ceratophyllum demersum). Geochemical speciation modelling also indicated that the primary dissolved Al species was Al(OH)4- under the observed daytime conditions, and that Fe was primarily in colloidal form Fe(OH)3(am), which was confirmed through diffusive gradients in thin films (DGT) measurements. A diel sampling experiment also confirmed that as alum-dosed water enters the lake, daytime biogeochemical conditions driven by C. demersum alter physicochemical water properties from ~pH 6 to ~pH 9, and DO (supersaturation) via increased photosynthesis. Diel-fluxes in geochemical solubility (Al, Fe, and DRP) responded to physicochemical shifts (pH, O2,) and demonstrate that C. demersum have the capacity to influence Al solubility and DRP availability within Te Wairoa Bay. This work highlights the complexity of biogeochemical processes within aquatic freshwater ecosystems. Moreover, the results emphasise the need to account for the significant spatial and temporal heterogeneity of physicochemical parameters in the development of effective lake remediation strategies

Reports of planetary geology and geophysics program, 1988

This is a compilation of abstracts of reports from Principal Investigators of NASA's Planetary Geology and Geophysics Program, Office of Space Science and Applications. The purpose is to document in summary form research work conducted in this program during 1988. Each report reflects significant accomplishments within the area of the author's funded grant or contract