Spectral optimization for constituent retrieval in Case 2 waters II: Validation study in the Chesapeake Bay

Coastal waters (Case 2) are generally more optically complex than oceanic waters and contain much higher quantities of colored detrital matter (CDM, a combination of dissolved organic matter and detrital particulates) as well as suspended sediment. Exclusion of CDM in the retrieval can lead to an overestimation of chlorophyll a concentration (C). We present a validation of a Case 2 version of the coupled spectral optimization algorithm (SOA) for simultaneous atmospheric correction and water parameter retrieval using Sea-viewing Wide Field-of-view Sensor (SeaWiFS) satellite ocean color data. Modeling of water constituents uses the Garver, Siegel and Maritorena (GSM) semi-analytic bio-optical model locally tuned for Chesapeake Bay. This includes a parameterization for CDM through its absorption spectrum.SOA-retrieved C and CDM are compared with in situ measurements in Chesapeake Bay. Results are also compared with output from two alternate models 1) the standard algorithm (Std) and 2) the standard atmospheric correction combined with the locally tuned GSM model (StdGSM). The comparisons indicate that the SOA is a viable alternative to both given models in Chesapeake Bay. In contrast, StdGSM appears to require improvement before it can be considered for operational use in these waters. Perhaps the most important result is the high-quality of CDM retrievals with the SOA. They suggest that there is value added using the SOA method in Chesapeake waters, as the Std method does not retrieve CDM. In a companion paper we describe in detail the model implementation, and its accuracy and limitations when applied to the Chesapeake Bay.     

Atmospheric correction of ENVISAT/MERIS data over inland waters: Validation for European lakes

Traditional methods for aerosol retrieval and atmospheric correction of remote sensing data over water surfaces are based on the assumption of zero water reflectance in the near-infrared. Another type of approach which is becoming very popular in atmospheric correction over water is based on the simultaneous retrieval of atmospheric and water parameters through the inversion of coupled atmospheric and bio-optical water models. Both types of approaches may lead to substantial errors over optically-complex water bodies, such as case II waters, in which a wide range of temporal and spatial variations in the concentration of water constituents is expected. This causes the water reflectance in the near-infrared to be non-negligible, and that the water reflectance response under extreme values of the water constituents cannot be described by the assumed bio-optical models. As an alternative to these methods, the SCAPE-M atmospheric processor is proposed in this paper for the automatic atmospheric correction of ENVISAT/MERIS data over inland waters. A-priori assumptions on the water composition and its spectral response are avoided by SCAPE-M by calculating reflectance of close-to-land water pixels through spatial extension of atmospheric parameters derived over neighboring land pixels. This approach is supported by the results obtained from the validation of SCAPE-M over a number of European inland water validation sites which is presented in this work. MERIS-derived aerosol optical thickness, water reflectance and water pigments are compared to in-situ data acquired concurrently to MERIS images in 20 validation match-ups. SCAPE-M has also been compared to specific processors designed for the retrieval of lake water constituents from MERIS data. The performance of SCAPE-M to reproduce ground-based measurements under a range of water types and the ability of MERIS data to monitor chlorophyll-a and phycocyanin pigments using semiempirical algorithms after SCAPE-M processing are discussed. It has been found that SCAPE-M is able to provide high accurate water reflectance over turbid waters, outperforming models based on site-specific bio-optical models, although problems of SCAPE-M to cope with clear waters in some cases have also been identified.     

Coupling and comparing a spatially- and temporally-detailed eutrophication model with an ecosystem network model: An initial application to Chesapeake Bay

Coastal waters are modeled for a variety of purposes including eutrophication remediation and fisheries management. Combining these two approaches provides insights which are not available from either approach independently. Coupling is confounded, however, by differences in model formulations and “currencies.” We present here an initial coupling of a spatially- and temporally-detailed eutrophication model, CE-QUAL-ICM, with a network fisheries model, Ecopath. We list commonalities between the models and present algorithms and software for the exchange of information. The models are applied to the central portion of Chesapeake Bay for a contemporary summer period. After comparison of the representations of Chesapeake Bay by the two models, an illustrative example one-way, off-line, coupling is presented. In an initial examination of a 20% increase in predation on phytoplankton by a small, highly-exploited fish (Atlantic menhaden, Brevoortia tyrannus), computed reduction in phytoplankton biomass is accompanied by increased production due to enhanced nutrient recycling. Minimal impact on the structure of the food web or on biomass of higher-trophic level organisms is computed. The algorithms and software can be adapted to alternate eutrophication models and Ecopath applications and provide the first, necessary, steps for subsequent coupling with the time-variable Ecosim model.     

Coupled retrieval of aerosol optical thickness, columnar water vapor and surface reflectance maps from ENVISAT/MERIS data over land

An algorithm for the derivation of atmospheric parameters and surface reflectance data from MEdium Resolution Imaging Specrometer Instrument (MERIS) on board ENVIronmental SATellite (ENVISAT) images has been developed. Geo-rectified aerosol optical thickness (AOT), columnar water vapor (CWV) and spectral surface reflectance maps are generated from MERIS Level-1b data over land. The algorithm has been implemented so that AOT, CWV and reflectance products are provided on an operational manner, making no use of ancillary parameters apart from those attached to MERIS products. For this reason, it has been named Self-Contained Atmospheric Parameters Estimation from MERIS data (SCAPE-M). The fundamental basis of the algorithm and applicable error figures are presented in the first part of this paper. In particular, errors of ± 0.03, ± 4% and ± 8% have been estimated for AOT, CWV and surface reflectance retrievals, respectively, by means of a sensitivity analysis based on a synthetic data set simulated under a usual MERIS scene configuration over land targets. The assumption of a fixed aerosol model, the coarse spatial resolution of the AOT product and the neglection of surface reflectance directional effects were also identified as limitations of SCAPE-M. Validation results are detailed in the second part of the paper. Comparison of SCAPE-M AOT retrievals with data from AErosol RObotic NETwork (AERONET) stations showed an average Root Mean Square Error (RMSE) of 0.05, and an average correlation coefficient R² of about 0.7-0.8. R² values grew up to more than 0.9 in the case of CWV after comparison with the same stations. A good correlation is also found with the MERIS Level-2 ESA CWV product. Retrieved surface reflectance maps have been successfully compared with reflectance data derived from the Compact High Resolution Imaging Spectrometer (CHRIS) on board the PRoject for On-Board Autonomy (PROBA) in the first place. Reflectance retrievals have also been compared with reflectance data derived from MERIS images by the Bremen AErosol Retrieval (BAER) method. A good correlation in the red and near-infrared bands was found, although a considerably higher proportion of pixels was successfully processed by SCAPE-M.     

Modeling spectral reflectance of optically complex waters using bio-optical measurements from Tokyo Bay

This study presents an approach for optimally parameterizing a reflectance model. A parameterization scheme is realized based on a comprehensive bio-optical data set, including subsurface downwelling and upwelling irradiance spectra, absorption spectra of particle and dissolved substances, as well as chlorophyll and total suspended matter concentrations at 45 stations near Tokyo Bay between 1982 and 1984. The irradiance reflectance model is implemented with three-component inherent optical property submodels.In this parameterization scheme, an unsupervised classification was applied in the hyper-spectral space of reflectance, leading to three spectrally distinct optical water types. The reflectance model was parameterized for the entire data set, and then parameterized for each of the water types. The three sets of type-specific model parameters, which define corresponding IOP submodels, are believed to accommodate differences in the optical properties of the in-water constituents. The parameterized reflectance model was evaluated by both reconstructing measured reflectance spectra and solving for the nonlinear inverse problem to retrieve in-water constituent concentrations. The model accuracy was significantly improved in the forward direction for classified waters over that of non-classified waters, but no significant improvement was achieved in the retrieval accuracy (inverse direction). A larger data set with greater resolution of constituent inherent optical properties would likely improve the modeling results.     

Retrieval of oceanic chlorophyll concentration with relevance vector machines

In this communication, we evaluate the performance of the relevance vector machine (RVM) for the estimation of biophysical parameters from remote sensing data. For illustration purposes, we focus on the estimation of chlorophyll-a concentrations from remote sensing reflectance just above the ocean surface. A variety of bio-optical algorithms have been developed to relate measurements of ocean radiance to in situ concentrations of phytoplankton pigments, and ultimately most of these algorithms demonstrate the potential of quantifying chlorophyll-a concentrations accurately from multispectral satellite ocean color data. Both satellite-derived data and in situ measurements are subject to high levels of uncertainty. In this context, robust and stable non-linear regression methods that provide inverse models are desirable.Lately, the use of the support vector regression (SVR) has produced good results in inversion problems, improving state-of-the-art neural networks. However, the SVR has some deficiencies, which could be theoretically alleviated by the RVM. In this paper, performance of the RVM is evaluated in terms of accuracy and bias of the estimations, sparseness of the solutions, robustness to low number of training samples, and computational burden. In addition, some theoretical issues are discussed, such as the sensitivity to training parameters setting, kernel selection, and confidence intervals on the predictions.Results suggest that RVMs offer an excellent trade-off between accuracy and sparsity of the solution, and become less sensitive to the selection of the free parameters. A novel formulation of the RVM that incorporates prior knowledge of the problem is presented and successfully tested, providing better results than standard RVM formulations, SVRs, neural networks, and classical bio-optical models for SeaWIFS, such as Morel, CalCOFI and OC2/OC4 models.     

Simultaneous retrieval of oceanic and atmospheric parameters for ocean color imagery by spectral optimization: a validation

We report application and validation of a spectral optimization algorithm for processing SeaWiFS data in Case 1 waters. The algorithm couples a simplified aerosol model with a sophisticated water-reflectance model to simultaneously retrieve both atmospheric and ocean parameters. Two of the retrieved ocean properties—the absorption coefficient of colored detrital material and the chlorophyll a concentration—are validated by comparison with “surface” truth obtained with airborne and space-borne sensors. We show that employing a more complete water reflectance model significantly improves the decoupling between the oceanic and atmospheric optical signals. Methodologies for applying the algorithm to Case 2 waters and for delineating terrestrial vs. marine chromophoric dissolved organic matter (CDOM) are suggested.     

Assessment of water quality in Lake Garda (Italy) using Hyperion

For testing the integration of the remote sensing related technologies into the water quality monitoring programs of Lake Garda (the largest Italian lake), the spatial and spectral resolutions of Hyperion and the capability of physics-based approaches were considered highly suitable. Hyperion data were acquired on 22nd July 2003 and water quality was assessed (i) defining a bio-optical model, (ii) converting the Hyperion at-sensor radiances into subsurface irradiance reflectances, and (iii) adopting a bio-optical model inversion technique. The bio-optical model was parameterised using specific inherent optical properties of the lake and light field variables derived from a radiative transfer numerical model. A MODTRAN-based atmospheric correction code, complemented with an air/water interface correction was used to convert Hyperion at-sensor radiances into subsurface irradiance reflectance values. These reflectance values were comparable to in situ reflectance spectra measured during the Hyperion overpass, except at longer wavelengths (beyond 700 nm), where reflectance values were contaminated by severe atmospheric adjacency effects. Chlorophyll-a and tripton concentrations were retrieved by inverting two Hyperion bands selected using a sensitivity analysis applied to the bio-optical model. The sensitivity analysis indicated that the assessment of coloured dissolved organic matter was not achievable in this study due to the limited coloured dissolved organic matter concentration range of the lake, resulting in reflectance differences below the environmental measurement noise of Hyperion. The chlorophyll-a and tripton image-products were compared to in situ data collected during the Hyperion overpass, both by traditional sampling techniques (8 points) and by continuous flow-through systems (32 km). For chlorophyll-a the correlation coefficient between in situ point stations and Hyperion-inferred concentrations was 0.77 (data range from 1.30 to 2.16 mg m^− 3). The Hyperion-derived chlorophyll-a concentrations also match most of the flow-through transect data. For tripton, the validation was constrained by variable re-suspension phenomena. The correlation coefficient between in situ point stations and Hyperion-derived concentrations increased from 0.48 to 0.75 (data range from 0.95 to 2.13 g m^− 3) if the sampling data from the re-suspension zone was avoided. The comparison of Hyperion-derived tripton concentrations and flow-through transect data exhibited a similar mismatch. The results of this research suggest further studies to address compatibilities of validation methods for water body features with a high rate of change, and to reduce the contamination by atmospheric adjacency effects on Hyperion data at longer wavelengths in Alpine environment. The transferability of the presented method to other sensors and the ability to assess water quality independent from in situ water quality data, suggest that management relevant applications for Lake Garda (and other subalpine lakes) could be supported by remote sensing.     

Regional and seasonal variability of chlorophyll-a in Chesapeake Bay as observed by SeaWiFS and MODIS-Aqua

Concentrations of the phytoplankton pigment chlorophyll-a (C_a) provide indicators of nutrient over-enrichment that has negatively affected Chesapeake Bay, U.S.A. C_a time-series from the National Aeronautics and Space Administration (NASA) Sea-viewing Wide Field-of-view Sensor (SeaWiFS) and Moderate Resolution Imaging Spectroradiometer aboard the Aqua spacecraft (MODIS-Aqua) provide observations on temporal and spatial scales that far exceed current field and aircraft sampling strategies. These sensors provide consistent, frequent, and high density data to potentially complement ongoing Bay monitoring activities. We used the in situ Water Quality Monitoring Data set of the Chesapeake Bay Program to evaluate decade-long time-series of SeaWiFS and MODIS-Aqua C_a retrievals in the Bay. The accuracy of the retrievals generally degraded with increasing latitude as the optical complexity increases northward. C_a derived using empirical (“band ratio”) algorithms overestimated in situ measurements by 10-50 and 40-100% for SeaWiFS and MODIS-Aqua, respectively, but with limited variability. C_a derived using spectral-matching algorithms showed less bias for both sensors, but with significant variability and sensitivity to radiometric errors. Regionally-tuned empirical algorithms performed best throughout the Bay, offering a combination of reasonable accuracy and high spatial coverage. The radiometric spectral resolution used as input to the algorithms strongly influenced the quality of C_a retrievals from both sensors. These results establish a baseline quantification of algorithm and sensor performance in a variable and stressed ecosystem against which novel approaches might be compared.     

Comparison of three SeaWiFS atmospheric correction algorithms for turbid waters using AERONET-OC measurements

The use of satellites to monitor the color of the ocean requires effective removal of the atmospheric signal. This can be performed by extrapolating the aerosol optical properties in the visible from the near-infrared (NIR) spectral region assuming that the seawater is totally absorbant in this latter part of the spectrum. However, the non-negligible water-leaving radiance in the NIR which is characteristic of turbid waters may lead to an overestimate of the atmospheric radiance in the whole visible spectrum with increasing severity at shorter wavelengths. This may result in significant errors, if not complete failure, of various algorithms for the retrieval of chlorophyll-a concentration, inherent optical properties and biogeochemical parameters of surface waters.This paper presents results of an inter-comparison study of three methods that compensate for NIR water-leaving radiances and that are based on very different hypothesis: 1) the standard SeaWiFS algorithm (Stumpf et al., 2003; Bailey et al., 2010) based on a bio-optical model and an iterative process; 2) the algorithm developed by Ruddick et al. (2000) based on the spatial homogeneity of the NIR ratios of the aerosol and water-leaving radiances; and 3) the algorithm of Kuchinke et al. (2009) based on a fully coupled atmosphere-ocean spectral optimization inversion. They are compared using normalized water-leaving radiance nL_w in the visible. The reference source for comparison is ground-based measurements from three AERONET-Ocean Color sites, one in the Adriatic Sea and two in the East Coast of USA.Based on the matchup exercise, the best overall estimates of the nL_w are obtained with the latest SeaWiFS standard algorithm version with relative error varying from 14.97% to 35.27% for λ = 490 nm and λ = 670 nm respectively. The least accurate estimates are given by the algorithm of Ruddick, the relative errors being between 16.36% and 42.92% for λ = 490 nm and λ = 412 nm, respectively. The algorithm of Kuchinke appears to be the most accurate algorithm at 412 nm (30.02%), 510 (15.54%) and 670 nm (32.32%) using its default optimization and bio-optical model coefficient settings.Similar conclusions are obtained for the aerosol optical properties (aerosol optical thickness τ(865) and the Ångström exponent, α(510, 865)). Those parameters are retrieved more accurately with the SeaWiFS standard algorithm (relative error of 33% and 54.15% for τ(865) and α(510, 865)).A detailed analysis of the hypotheses of the methods is given for explaining the differences between the algorithms. The determination of the aerosol parameters is critical for the algorithm of Ruddick et al. (2000) while the bio-optical model is critical for the algorithm of Stumpf et al. (2003) utilized in the standard SeaWiFS atmospheric correction and both aerosol and bio-optical model for the coupled atmospheric-ocean algorithm of Kuchinke. The Kuchinke algorithm presents model aerosol-size distributions that differ from real aerosol-size distribution pertaining to the measurements. In conclusion, the results show that for the given atmospheric and oceanic conditions of this study, the SeaWiFS atmospheric correction algorithm is most appropriate for estimating the marine and aerosol parameters in the given turbid waters regions.     

物流监控中的蜂窝无线定位

利用现有的蜂窝网络,通过测量车栽移动终端的位置特征参数TOA和TDOA,结合扇区信息进行数据融合,构造一种基于移动台位置的动态定位算法,提出一种改进型数据融合模型.详细分析基于TDOA的Chan算法及特点,对Chan算法计算过程的中间阶段数据进行数据融合,通过定义可信度函数,从而实现一种蜂窝网络定位新方法.本文的定位方法为物流企业业务的监控和优化提供有力的技术支持.     

Evaluation of shortwave infrared atmospheric correction for ocean color remote sensing of Chesapeake Bay

The NASA Moderate Resolution Imaging Spectroradiometer onboard the Aqua platform (MODIS-Aqua) provides a viable data stream for operational water quality monitoring of Chesapeake Bay. Marine geophysical products from MODIS-Aqua depend on the efficacy of the atmospheric correction process, which can be problematic in coastal environments. The operational atmospheric correction algorithm for MODIS-Aqua requires an assumption of negligible near-infrared water-leaving radiance, nL_w(NIR). This assumption progressively degrades with increasing turbidity and, as such, methods exist to account for non-negligible nL_w(NIR) within the atmospheric correction process or to use alternate radiometric bands where the assumption is satisfied, such as those positioned within shortwave infrared (SWIR) region of the spectrum. We evaluated a decade-long time-series of nL_w(λ) from MODIS-Aqua in Chesapeake Bay derived using NIR and SWIR bands for atmospheric correction. Low signal-to-noise ratios (SNR) for the SWIR bands of MODIS-Aqua added noise errors to the derived radiances, which produced broad, flat frequency distributions of nL_w(λ) relative to those produced using the NIR bands. The SWIR approach produced an increased number of negative nL_w(λ) and decreased sample size relative to the NIR approach. Revised vicarious calibration and regional tuning of the scheme to switch between the NIR and SWIR approaches may improve retrievals in Chesapeake Bay, however, poor SNR values for the MODIS-Aqua SWIR bands remain the primary deficiency of the SWIR-based atmospheric correction approach.     

The determination of dust cloud altitudes from a satellite using hyperspectral measurements in the gaseous absorption band

An approach to derive dust layer optical thickness and top height using top-of-atmosphere (TOA) reflectance in the oxygen A-band is introduced. The algorithm is similar to that developed by the authors for the case of water clouds. It is based on the fitting of spectral TOA reflectance measurements in the narrow band around 760 nm using results of the exact radiative transfer calculations for a given dust layer model. The accuracy of the technique with respect to the uncertainty in a priori assumption of the dust single-scattering albedo is discussed. The algorithm is applied to satellite hyperspectral measurements over the Atlantic Ocean.     

Top-of-atmosphere flux retrievals from CERES using artificial neural networks

The Clouds and the Earth's Radiant Energy System (CERES) instruments on the Terra spacecraft provide accurate shortwave (SW), longwave (LW) and window (WN) region top-of-atmosphere (TOA) radiance measurements from which TOA radiative flux values are obtained by applying Angular Distribution Models (ADMs). These models are developed empirically as functions of the surface and cloud properties provided by coincident high-resolution imager measurements over CERES field-of-view. However, approximately 5.6% of the CERES/Terra footprints lack sufficient imager information for a reliable scene identification. To avoid any systematic biases in regional mean radiative fluxes, it is important to provide TOA fluxes for these footprints. For this purpose, we apply a feedforward error-backpropagation Artificial Neural Network (ANN) technique to reproduce CERES/Terra ADMs relying only on CERES measurements. All-sky ANN-based angular distribution models are developed for 10 surface types separately for shortwave, longwave and window TOA flux retrievals. To optimize the ANN performance, we use a partially connected first hidden neuron layer and compact training sets with reduced data noise. We demonstrate the performance of the ANN-based ADMs by comparing TOA fluxes inferred from ANN and CERES anisotropic factors. The global annual average bias in ANN-derived fluxes relative to CERES is less than 0.5% for all ANN scene types. The maximum bias occurs over sea ice and permanent snow surfaces. For all surface types, instantaneous ANN-derived TOA fluxes are self-consistent in viewing zenith angle to within 9% for shortwave, 3.5% and 3% longwave daytime and nighttime, respectively.     

LANDSAT test of diffuse reflectance models for aquatic suspended solids measurement

LANDSAT radiance data were used to test mathematical models relating diffuse reflectance to aquatic suspended solids concentration. Digital CCT data for LANDSAT passes over the Bay of Fundy, Nova Scotia were analyzed on a General Electric Co. Image 100 multispectral analysis system. Three data sets were studied separately and together in all combinations with and without solar angle correction. Statistical analysis and chromaticity analysis show that a nonlinear relationship between LANDSAT radiance and suspended solids concentration is better at curve-fitting than a linear relationship. In particular, the quasi-single-scattering diffuse reflectance model developed by Gordon and coworkers is corroborated. The Gordon model applied to 33 points of MSS 5 data combined from three dates produced r = 0.98.     

Web-enabled Landsat Data (WELD): Landsat ETM+ composited mosaics of the conterminous United States

Since January 2008, the U.S. Department of Interior / U.S. Geological Survey have been providing free terrain-corrected (Level 1T) Landsat Enhanced Thematic Mapper Plus (ETM+) data via the Internet, currently for acquisitions with less than 40% cloud cover. With this rich dataset, temporally composited, mosaics of the conterminous United States (CONUS) were generated on a monthly, seasonal, and annual basis using 6521 ETM+ acquisitions from December 2007 to November 2008. The composited mosaics are designed to provide consistent Landsat data that can be used to derive land cover and geo-physical and bio-physical products for detailed regional assessments of land-cover dynamics and to study Earth system functioning. The data layers in the composited mosaics are defined at 30 m and include top of atmosphere (TOA) reflectance, TOA brightness temperature, TOA normalized difference vegetation index (NDVI), the date each composited pixel was acquired on, per-band radiometric saturation status, cloud mask values, and the number of acquisitions considered in the compositing period. Reduced spatial resolution browse imagery, and top of atmosphere 30 m reflectance time series extracted from the monthly composites, capture the expected land surface phenological change, and illustrate the potential of the composited mosaic data for terrestrial monitoring at high spatial resolution. The composited mosaics are available in 501 tiles of 5000 × 5000 30 m pixels in the Albers equal area projection and are downloadable at http://landsat.usgs.gov/WELD.php. The research described in this paper demonstrates the potential of Landsat data processing to provide a consistent, long-term, large-area, data record.     

Revisiting empirical ocean-colour algorithms for remote estimation of chlorophyll-a content on a global scale

Remote-sensing data can be useful for investigating the bio-optical properties of the ocean. Among these bio-optical properties, chlorophyll-a content is of great importance. The standard NASA empirical ocean-colour (OC) algorithms are used widely to estimate global chlorophyll-a content. Despite their simplicity and effectiveness, these regression-based models have two shortcomings that we investigate here: (1) the general form of the models is a fourth-order polynomial that results in multicollinearity, and (2) the models have the same parameters for all ocean regions (i.e. they use global approaches). To resolve the first issue, we use partial least squares (PLS), which allows for an orthogonal transformation such that the covariance between the transformed independent variables and the dependent variable is maximized. To investigate the second issue, we use geographically weighted regression (GWR) to reveal the spatial variation of estimated parameters, demonstrating how the global model underperforms in some locations. GWR results show that model coefficients vary substantially between eastern and western portions of the same ocean basin. By including sea-surface temperature (SST) as an additional independent variable in the PLS model, we also develop a new approach that provides additional explanatory power and makes the global estimation of chlorophyll-a content more valid.     

Remote sensing of cyanobacteria-dominant algal blooms and water quality parameters in Zeekoevlei, a small hypertrophic lake, using MERIS

Eutrophication and cyanobacterial algal blooms present an increasing threat to the health of freshwater ecosystems and to humans who use these resources for drinking and recreation. Remote sensing is being used increasingly as a tool for monitoring these phenomena in inland and near-coastal waters. This study uses the Medium Resolution Imaging Spectrometer (MERIS) to view Zeekoevlei, a small hypertrophic freshwater lake situated on the Cape Flats in Cape Town, South Africa, dominated by Microcystis cyanobacteria. The lake's small size, highly turbid water, and covariant water constituents present a challenging case for both algorithm development and atmospheric correction. The objectives of the study are to assess the optical properties of the lake, to evaluate various atmospheric correction procedures, and to compare the performance of empirical and semi-analytical algorithms in hypertrophic water. In situ water quality parameter and radiometric measurements were made simultaneous to MERIS overpasses. Upwelling radiance measurements at depth 0.66 m were corrected for instrument self-shading and processed to water-leaving reflectance using downwelling irradiance measurements and estimates of the vertical attenuation coefficient for upward radiance, K_u, generated from a simple bio-optical model estimating the total absorption, a(λ), and backscattering coefficients, b_b(λ). The normalised water-leaving reflectance was used for assessing the accuracy of image-based Dark Object Subtraction and 6S Radiative Transfer Code atmospheric correction procedures applied to MERIS. Empirical algorithms for estimating chlorophyll a (Chl a), Total Suspended Solids (TSS), Secchi Disk depth (z_SD) and absorption by CDOM (a_CDOM) were derived from simultaneously collected in situ and MERIS measurements. The empirical algorithms gave high correlation coefficient values, although they have a limited ability to separate between signals from covariant water constituents. The MERIS Neural Network algorithms utilised in the standard Level 2 Case 2 waters product and Eutrophic Lakes processor were also used to derive water constituent concentrations. However, these failed to produce reasonable comparisons with in situ measurements owing to the failure of atmospheric correction and divergence between the optical properties and ranges used to train the algorithms and those of Zeekoevlei. Maps produced using the empirical algorithms effectively show the spatial and temporal variability of the water quality parameters during April 2008. On the basis of the results it is argued that MERIS is the current optimal sensor for frequent change detection applications in inland waters. This study also demonstrates the considerable potential value for simple TOA algorithms for hypertrophic systems. It is recommended that regional algorithm development be prioritized in southern Africa and that remote sensing be integrated into future operational water quality monitoring systems.     

Fields of non-linear regression models for atmospheric correction of satellite ocean-color imagery

Remote sensing of ocean color from space, a problem that consists of retrieving spectral marine reflectance from spectral top-of-atmosphere reflectance, is considered as a collection of similar inverse problems continuously indexed by the angular variables influencing the observation process. A general solution is proposed in the form of a field of non-linear regression models over the set T of permitted values for the angular variables, i.e., as a map from T to some function space. Each value of the field is a regression model that performs a direct mapping from the top-of-atmosphere reflectance to the marine reflectance. Since the spectral components of the field take values in the same variable vector space, the retrievals in individual spectral bands are not independent, i.e., the solution is not just a juxtaposition of independent models for each spectral band. A scheme based on ridge functions is developed to approximate this solution to an arbitrary accuracy, and is applied to the retrieval of marine reflectance in Case 1 waters, for which optical properties are only governed by biogenic content. The statistical models are evaluated on synthetic data as well as actual data originating from the SeaWiFS instrument, taking into account noise in the data. Theoretical performance is good in terms of accuracy, robustness, and generalization capabilities, suggesting that the function field methodology might improve atmospheric correction in the presence of absorbing aerosols and provide more accurate estimates of marine reflectance in productive waters. When applied to SeaWiFS imagery acquired off California, the function field methodology gives generally higher estimates of marine reflectance than the standard SeaDAS algorithm, but the values are more realistic.