Utility of remote sensing data in retrieval of water quality consituents concentrations in coastal water of New Jersey

Three important optical properties used for monitoring coastal water quality are the concentrations of chlorophyll (CHL), color dissolved organic matter (CDOM) and total suspended materials (TSM). Ocean color remote sensing, a technique to collect color data by detection of upward radiance from a distance (Bukata et al.,1995), provides a synoptic view for determining these concentrations from upwelling radiances. In the open ocean (Case-I), it is not difficult to derive empirical algorithms relating the received radiances to surface concentrations of water quality parameters. In coastal waters (Case-Il), there are serious unresolved problems in extracting chlorophyll concentration because of high concentration of suspended particles (Gordon and Morel, 1983). There are three basic approaches to estimate optical water quality parameters from remotely sensed spectral data based on the definitions given by Morel & Gordon (1980): (1) an empirical method, in which statistical relationships between the upward radiance at the sea surface and the quantity of interest are taken into account; (2) a semiempirical method, in which the spectral characteristics of the parameters of interest are known and some modeling of the physics is introduced; and (3) an analytical method, in which radiative transfer models are used to extract the inherent optical properties (lOPs) and a suite of analysis methods can be used to optimally retrieve the water constituents from the remotely sensed upwelling radiance or irradiance reflectance signal. The focus of this research is the modification and application of analytical and statistical algorithms to characterize the physically based surface spectral reflectance for the waters of the Hudson/Raritan Estuary and to retrieve the water constituent concentrations from the NASA Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and LIght Detection And Ranging (LIDAR) signals. The approaches used here are based on the unique capabilities of AVIRIS and LIDAR data which can potentially provide a better understanding of how sunlight interacts with estuarine/inland water, especially when complemented with in situ measurements for analysis of water quality parameters and eutrophication processes. The results of analysis in forms of thematic maps are then input into geographic information system (GIS) of the study site for use by water resource managers and planners for better monitoring and management of water quality condition

Improving Chlorophyll-a Estimation from Sentinel-2 (MSI) in the Barents Sea using Machine Learning

This article addresses methodologies for remote sensing of ocean Chlorophyll-a (Chl-a), with emphasis on the Barents Sea. We aim at improving the monitoring capacity by integrating in situ Chl-a observations and optical remote sensing to locally train machine learning (ML) models. For this purpose, in situ measurements of Chl-a ranging from 0.014–10.81 mg/m 3 , collected for the years 2016–2018, were used to train and validate models. To accurately estimate Chl-a, we propose to use additional information on pigment content within the productive column by matching the depth-integrated Chl-a concentrations with the satellite data. Using the optical images captured by the multispectral imager instrument on Sentinel-2 and the in situ measurements, a new spatial window-based match-up dataset creation method is proposed to increase the number of match-ups and hence improve the training of the ML models. The match-ups are then filtered to eliminate erroneous samples based on the spectral distribution of the remotely sensed reflectance. In addition, we design and implement a neural network model dubbed as the ocean color net (OCN), that has performed better than existing ML-based techniques, including the Gaussian process Regression (GPR), regionally tuned empirical techniques, including the ocean color (OC3) algorithm and the spectral band ratios, as well as the globally trained Case-2 regional/coast colour (C2RCC) processing chain model C2RCC-networks. The proposed OCN model achieved reduced mean absolute error compared to the GPR by 5.2%, C2RCC by 51.7%, OC3 by 22.6%, and spectral band ratios by 29%. Moreover, the proposed spatial window and depth-integrated match-up creation techniques improved the performance of the proposed OCN by 57%, GPR by 41.9%, OC3 by 5.3%, and spectral band ratio method by 24% in terms of RMSE compared to the conventional match-up selection approach

Relating an archive of in situ vertical chlorophyll-a profiles to concurrent remotely sensed surface data

Includes bibliographical references (leaves 56-67).Knowledge of the vertical distribution of phytoplankton in the upper ocean is essential for accurate estimates of primary production. Satellite remote sensing has given scientists an unprecedented view of near-surface chlorophyll distribution and other surface conditions, including sea surface temperature and wind data, from regional to global scales but little information on the dynamics below the surface. As a result estimates of global production tend to use regional profile averages but these methods oversimplify the smaller scale dynamics, particularly in coastal regions where productivity is highly variable on time scales of weeks. A pilot study by computer science honours students in 2006 showed the viability of using a Dynamic Bayesian Network (DBN) in predicting a representative profile per pixel of a satellite map based on a database of time series satellite surface data. In this study, 5813 in situ profiles were obtained from the highly dynamic upwelling region around the southwestern coastline of southern Africa. The samples were collected between 1988 and 2006 between the coast and the continental slope. The region was divided into three sub-regions according to biophysical processes: the west Coast; the west Agulhas Bank; and the east Agulhas Bank. Of the 5813 profiles, 5557 were included in the sub-regions. Two consecutive processes were then applied to the profile database. First, the profiles were clustered using a k-means clustering program which produced 16 representative clusters

Remote Sensing of Water Quality Parameters over Lake Balaton by Using Sentinel-3 OLCI

Source at https://doi.org/10.3390/w10101428.The Ocean and Land Color Instrument (OLCI) onboard Sentinel 3A satellite was launched in February 2016. Level 2 (L2) products have been available for the public since July 2017. OLCI provides the possibility to monitor aquatic environments on 300 m spatial resolution on 9 spectral bands, which allows to retrieve detailed information about the water quality of various type of waters. It has only been a short time since L2 data became accessible, therefore validation of these products from different aquatic environments are required. In this work we study the possibility to use S3 OLCI L2 products to monitor an optically highly complex shallow lake. We test S3 OLCI-derived Chlorophyll-a (Chl-a), Colored Dissolved Organic Matter (CDOM) and Total Suspended Matter (TSM) for complex waters against in situ measurements over Lake Balaton in 2017. In addition, we tested the machine learning Gaussian process regression model, trained locally as a potential candidate to retrieve water quality parameters. We applied the automatic model selection algorithm to select the combination and number of spectral bands for the given water quality parameter to train the Gaussian Process Regression model. Lake Balaton represents different types of aquatic environments (eutrophic, mesotrophic and oligotrophic), hence being able to establish a model to monitor water quality by using S3 OLCI products might allow the generalization of the methodology

Optical remote sensing of water quality parameters retrieval in the Barents Sea

This thesis addresses various aspects of monitoring water quality indicators (WQIs) using optical remote sensing technologies. The dynamic nature of aquatic systems necessitate frequent monitoring at high spatial resolution. Machine learning (ML)-based algorithms are becoming increasingly common for these applications. ML algorithms are required to be trained by a significant amount of training data, and their accuracy depends on the performance of the atmospheric correction (AC) algorithm being used for correcting atmospheric effects. AC over open oceanic waters generally performs reasonably well; however, limitations still exist over inland and coastal waters. AC becomes more challenging in the high north waters, such as the Barents Sea, due to the unique in-water optical properties at high latitudes, long ray pathways, as well as the scattering of light from neighboring sea ice into the sensors’ field of view adjacent to ice-infested waters. To address these challenges, we evaluated the performances of state-of-the-art AC algorithms applied to the high-resolution satellite sensors Landsat-8 Operational Land Imager (OLI) and Sentinel-2 Multispectral Instrument (MSI), both for high-north (Paper II) and for global inland and coastal waters (Paper III). Using atmospherically corrected remote sensing reflectance (Rrs ) products, estimated after applying the top performing AC algorithm, we present a new bandpass adjustment (BA) method for spectral harmonization of Rrs products from OLI and MSI. This harmonization will enable an increased number of ocean color (OC) observations and, hence, a larger amount of training data. The BA model is based on neural networks (NNs), which perform a pixel-by-pixel transformation of MSI-derived Rrs to that of OLI equivalent for their common bands. In addition, to accurately retrieve concentrations of Chlorophyll-a (Chl-a) and Color Dissolved Organic Matter (CDOM) from remotely sensed data, we propose in the thesis (Paper 1) an NN-based WQI retrieval model dubbed Ocean Color Net (OCN). Our results indicate that Rrs retrieved via the Acolite Dark Spectrum Fitting (DSF) method is in best agreement with in-situ Rrs observations in the Barents Sea compared to the other methods. The median absolute percentage difference (MAPD) in the blue-green bands ranges from 9% to 25%. In the case of inland and coastal waters (globally), we found that OC-SMART is the top performer, with MAPD Rrs products for varying optical regimes than previously presented methods. Additionally, to improve the analysis of remote sensing spectral data, we introduce a new spatial window-based match-up data set creation method which increases the training data set and allows for better tuning of regression models. Based on comparisons with in-water measured Chl-a profiles in the Barents Sea, our analysis indicates that the MSI-derived Rrs products are more sensitive to the depth-integrated Chl-a contents than near-surface Chl-a values (Paper I). In the case of inland and coastal waters, our study shows that using combined OLI and BA MSI-derived Rrs match-ups results in considerable improvement in the retrieval of WQIs (Paper III). The obtained results for the datasets used in this thesis illustrates that the proposed OCN algorithm shows better performance in retrieving WQIs than other semi-empirical algorithms such as the band ratio-based algorithm, the ML-based Gaussian Process Regression (GPR), as well as the globally trained Case-2 Regional/Coast Colour (C2RCC) processing chain model C2RCC-networks, and OC-SMART. The work in this thesis contributes to ongoing research in developing new methods for merging data products from multiple OC missions for increased coverage and the number of optical observations. The developed algorithms are validated in various environmental and aquatic conditions and have the potential to contribute to accurate and consistent retrievals of in-water constituents from high-resolution satellite sensors

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

Remote sensing of phytoplankton in the Arctic Ocean : development, tuning and evaluation of new algorithms

Thèse en cotutelle : Université Laval, Québec, Canada, Philosophiæ doctor (Ph. D.) et Wuhan University, Wuhan, Chine.Au cours des dernières décennies, l'augmentation de la production primaire (PP) dans l'océan Arctique (AO) a en partie été associée à une augmentation de la biomasse phytoplanctonique, comme l'ont montré des études de télédétection. La concentration en chlorophylle a (Chl), un indicateur de la biomasse phytoplanctonique, est un facteur clé qui peut biaiser les estimations de la PP quand elle comporte des erreurs de mesure. En d'autres mots, une estimation précise de la Chl est cruciale pour améliorer notre connaissance de l'écosystème marin et de sa réponse au changement climatique en cours. Cependant, la télédétection de la couleur de l'océan dans l'océan Arctique présente plusieurs défis. Tout d'abord, il est bien connu que l'échec des algorithmes standards de la couleur de l'océan dans l'AO est dû à l'interférence des matières colorées et détritiques (CDM) dans le spectre visible, mais comment et dans quelle mesure cela va biaiser l'estimation de la Chl reste inconnu. En outre, la Chl étant un facteur clé utilisé pour estimer la PP, la propagation des erreurs des estimations de la Chl aux estimations de la PP doit être évaluée. Le dernier mais le plus important est qu'un algorithme robuste avec une incertitude raisonnable, en particulier pour les eaux côtières complexes et productives, n'est pas encore disponible. Pour résoudre ces problèmes, dans cette étude, nous avons d'abord compilé une grande base de données bio-optiques in situ dans l'Arctique, à partir de laquelle nous avons évalué de manière approfondie les algorithmes de couleur de l'océan actuellement disponibles du point de vue des impacts des CDM. Nous avons constaté que plus le niveau de CDM par rapport à la Chl dans la colonne d'eau était élevé, plus il biaisait les estimations de la Chl. L'analyse de sensibilité des estimations de la PP sur la Chl a montré que l'erreur des estimations de la Chl était amplifiée de 7% lorsqu'elle était passée dans l'estimation du PP en utilisant un modèle de PP résolu spectralement et verticalement. En outre, pour obtenir de meilleurs résultats, nous avons optimisé un algorithme semi-analytique (Garver-Siegel-Maritorena, GSM) pour l'AO en ajoutant la bande spectrale de 620 nm qui est moins affectée par le CDM et le signal ici est généralement élevé pour les eaux riches en CDM, devenant anisi important pour le GSM afin d'obtenir des estimates précises de la Chl. Notre algorithme ajusté, GSMA, n'a amélioré la précision que de 8% pour l'AO, mais l'amélioration pour les eaux côtières a atteint 93%. Enfin, étant donné que les algorithmes qui n'exploitent pour la plupart que les parties bleue et verte du spectre visible sont problématiques pour les eaux très absorbantes/obscures, nous avons introduit un modèle d'émission de fluorescence pour tenir compte des propriétés bio-optiques du phytoplancton dans la partie rouge du spectre visible. En se couplant avec le GSMA, le nouvel algorithme à spectre complet, FGSM, a encore amélioré la précision des estimations de la Chl de ~44% pour les eaux eutrophes. À l'avenir, des couplages sont nécessaires à des fins de validation en ce qui concerne l'application satellitaire. De plus, de nouvelles approches pouvant être appliquées pour détecter le maximum de chlorophylle sous la surface (SCM), les efflorescences en bordure de glace et/ou sous la glace, les types fonctionnels de phytoplancton (PFT), sont attendues.In the recent decades, the raise of primary production (PP) in the Arctic Ocean (AO) is mainly driven by the increase of phytoplankton biomass as multiple remote sensing studies have suggested. Chlorophyll a concentration (Chl), a proxy of phytoplankton biomass, is a key factor that biases PP estimates. In terms of bottom-up control, accurate Chl estimate is crucial to improve our knowledge of marine ecosystem and its response to ongoing climate change. However, there are several challenges of ocean color remote sensing in the Arctic Ocean. Firstly, it is well known that the failure of standard ocean color algorithms in the AO is due to the interference of colored and detrital material (CDM) in the visible spectrum, but how and to what extend it will bias the estimation of Chl remains unknown. Besides, Chl as a key factor used to estimate PP, error propagation from Chl estimates to PP estimates needs to be assessed. The last but most important is that a robust algorithm with reasonable uncertainty, especially for the complex and productive coastal waters, is not yet available. To address these problems, in this study, we first compiled a large Arctic in situ bio-optical database, based on which we thoroughly evaluated presently available ocean color algorithms from a perspective of the impacts of CDM. We found that the higher the level of CDM relative to Chl in the water column, the larger it would bias Chl estimates. Sensitivity analysis of PP estimates on Chl showed that the error of Chl estimates was amplified within 7% when passed into the estimation of PP using a spectrally- and vertically-resolved PP model. Besides, to obtain better results, we tuned GSM for the AO by adding 620 waveband which is less affected by CDM and the signal here is generally high for CDM-rich waters thus become important for GSM to retrieve accurate Chl estimates. Our tuned algorithm, GSMA, merely improved the accuracy by 8% for the AO, but the improvement for coastal waters reached up to 93%. Finally, given that algorithms that only exploits visible spectrum are problematic for highly-absorbing/dark waters, we introduced the fluorescence emission model to account for the bio-optical properties of phytoplankton in the near infrared spectrum. By coupling with GSMA, the novel full-spectrally algorithm, FGSM, further improved the accuracy of Chl estimates by ~44% for eutrophic waters. In the future, matchups are needed for validation purposes with respect to satellite application. Moreover, new approaches that can be applied to detect subsurface chlorophyll maximum (SCM), ice-edge and/or under-ice blooms, phytoplankton functional types (PFT) and so on are expected

An Overview of Approaches and Challenges for Retrieving Marine Inherent Optical Properties from Ocean Color Remote Sensing

Ocean color measured from satellites provides daily global, synoptic views of spectral water-leaving reflectances that can be used to generate estimates of marine inherent optical properties (IOPs). These reflectances, namely the ratio of spectral upwelled radiances to spectral downwelled irradiances, describe the light exiting a water mass that defines its color. IOPs are the spectral absorption and scattering characteristics of ocean water and its dissolved and particulate constituents. Because of their dependence on the concentration and composition of marine constituents, IOPs can be used to describe the contents of the upper ocean mixed layer. This information is critical to further our scientific understanding of biogeochemical oceanic processes, such as organic carbon production and export, phytoplankton dynamics, and responses to climatic disturbances. Given their importance, the international ocean color community has invested significant effort in improving the quality of satellite-derived IOP products, both regionally and globally. Recognizing the current influx of data products into the community and the need to improve current algorithms in anticipation of new satellite instruments (e.g., the global, hyperspectral spectroradiometer of the NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission), we present a synopsis of the current state of the art in the retrieval of these core optical properties. Contemporary approaches for obtaining IOPs from satellite ocean color are reviewed and, for clarity, separated based their inversion methodology or the type of IOPs sought. Summaries of known uncertainties associated with each approach are provided, as well as common performance metrics used to evaluate them. We discuss current knowledge gaps and make recommendations for future investment for upcoming missions whose instrument characteristics diverge sufficiently from heritage and existing sensors to warrant reassessing current approaches