A simple band ratio technique to quantify the colored dissolved and detrital organic material from ocean color remotely sensed data

Spectral reflectances of the ocean, R, as derived from ocean color remote sensing data at four wavelengths (412, 443, 490, and 555 nm), can be used to form two ratios of spectral reflectance, namely R(412)/R(443), and R(490)/R(555), thereafter denoted R₄₄₃⁴¹² and R₅₅₅⁴⁹⁰. The former is mainly sensitive to the colored dissolved organic material (CDOM), albeit influenced by the algal content as depicted by the chlorophyll concentration, ([Chl]); in contrast, the latter is essentially depending on [Chl], although it is also influenced by CDOM. Therefore the signatures of CDOM and [Chl] which are not truly separable, can nevertheless be identified by considering simultaneously the two ratios. The concomitant variations in these ratios can be established via a bio-optical model developed for Case 1 waters. This model implicitly includes a “mean” relationship between CDOM and [Chl], and thus produces a unique curve relating R₄₄₃⁴¹² to R₅₅₅⁴⁹⁰. Deviations with respect to this mean relationship can be introduced through a factor Φ, with Φ > 1 (excess) or < 1 (deficit), applied to the CDOM-[Chl] ratio. A family of R₄₄₃⁴¹²-R₅₅₅⁴⁹⁰ curves is thus generated, in correspondence with the discrete values given to Φ; this “grid” (or numerically, a 2-D lookup table) allows the Φ-[Chl] couple to be unambiguously derived for any R₄₄₃⁴¹²-R₅₅₅⁴⁹⁰ couple. By applying this straightforward algorithm to actual reflectance ratios derived from ocean color imagery, the relative anomalies in CDOM with respect to its standard (Chl-related) values can be efficiently assessed. Within the global ocean (discarding the coastal zones), the Φ factor is widely varying, between at least ? and 3, and is roughly log-normally distributed around ~ 1 (no anomaly). The spatial distributions of the Φ factor in the whole ocean are strongly featured according to latitude, season, and hydrographic regimes, and these features are regularly reproducible, from year to year (2002-2007). This simple method is also validated against available in situ data, and its results compare favorably, for instance, to those of the GSM (Garver-Siegel-Maritorena) inversion method, in terms of retrieved CDOM concentrations and distribution patterns.     

An assessment of chlorophyll-a algorithms available for SeaWiFS in coastal and open areas of the Bay of Bengal and Arabian Sea

Three ocean colour algorithms, OC4v6, Carder and OC5 were tested for retrieving Chlorophyll-a (Chla) in coastal areas of the Bay of Bengal and open ocean areas of the Arabian Sea. Firstly, the algorithms were run using ~ 80 in situ Remote Sensing Reflectance, (R_rs(λ)) data collected from coastal areas during eight cruises from January 2000 to March 2002 and the output was compared to in situ Chla. Secondly, the algorithms were run with ~ 20 SeaWiFS R_rs(λ) and the results were compared with coincident in situ Chla. In both cases, OC5 exhibited the lowest log₁₀-RMS, bias, had a slope close to 1 and this algorithm appears to be the most accurate for both coastal and open ocean areas. Thirdly the error in the algorithms was regressed against Total Suspended Material (TSM) and Coloured Dissolved Organic Material (CDOM) data to assess the co-variance with these parameters. The OC5 error did not co-vary with TSM and CDOM. OC4v6 tended to over-estimate Chla > 2 mg m⁻³ and the error in OC4v6 co-varied with TSM. OC4v6 was more accurate than the Carder algorithm, which over-estimated Chla at concentrations > 1 mg m⁻³ and under-estimated Chla at values < 0.5 mg m⁻³. The error in Carder Chla also co-varied with TSM. The algorithms were inter-compared using > 5500 SeaWiFS R_rs(λ) data from coastal to offshore transects in the Northern Bay of Bengal. There was good agreement between OC4v6 and OC5 in open ocean waters and in coastal areas up to 2 mg m⁻³. There was a strong divergence between Carder and OC5 in open ocean and coastal waters. OC4v6 and Carder tended to over-estimate Chla in coastal areas by a factor of 2 to 3 when TSM > 25 g m⁻³. We strongly recommend the use of OC5 for coastal and open ocean waters of the Bay of Bengal and Arabian Sea. A Chla time series was generated using OC5 from 2000 to 2003, which showed that concentrations at the mouths of the Ganges reach a maxima (~ 5 mg m⁻³) in October and November and were 0.08 mg m⁻³ further offshore increasing to 0.2 mg m⁻³ during December. Similarly in early spring from February to March, Chla was 0.08 to 0.2 mg m⁻³ on the east coast of the Bay.     

Estimating phytoplankton biomass in coastal waters of Alaska using airborne remote sensing

Empirical airborne remote-sensing relationships were examined to estimate chlorophyll a concentration in the first optical depth (chl_FOD) of coastal waters of Afgonak/Kodiak Islands during July-August 2002. Band-ratio and spectral-curvature models were tested using satellite remote-sensing reflectance (R_rs(λ)) measurements. Additional shipboard and airborne R_rs(λ) data were also analysed to evaluate consistency of proposed chl_FOD-R_rs(λ) relationships. Validation of chlorophyll algorithms was performed using data collected in the northern-part of the Gulf of Alaska and Bering Sea during 1996, 2002, and 2003 cruises. Likewise, oceanographic conditions during the surveys were typified to interpret variability of chl_FOD fields. The SeaWiFS band-ratio algorithm OC2d was the most sensitive R_rs combination (R_rs(509)/R_rs(553)) to detect chl_FOD variability. Conversely, OC2a (R_rs(412)/R_rs(553)) had the lowest performance to derive chl_FOD values. No valid statistical regressions were established for spectral-curvature relationships in the blue spectrum (< 500 nm). Fertile waters (> 5 mg m^− 3) were preferentially located over shallow banks (∼50 m) and at the entrance of the bays. The approach used in this study to derive chl_FOD values could be universal for Alaskan coastal waters. However, chl_FOD-R_rs(λ) relationships must be calibrated locally for a given season.     

Remote chlorophyll-a retrieval in turbid, productive estuaries: Chesapeake Bay case study

Accurate remote assessment of phytoplankton chlorophyll a (chla) concentration is particularly challenging in turbid, productive waters. Recently a conceptual model containing reflectance in three spectral bands in the red and near infra-red range of the spectrum was suggested for retrieving chla concentrations in turbid productive waters; it was calibrated and validated in lakes and reservoirs in Nebraska and Iowa. The objective of this paper is to evaluate the performance of this three band model as well as its special case, the two-band model to estimate chla concentration in Chesapeake Bay, as representative of estuarine Case II waters, and to assess the accuracy of chla retrieval. To evaluate the model performance, dual spectroradiometers were used to measure subsurface spectral radiance reflectance in the visible and near infra-red range of the spectrum. Water samples were collected concurrently and contained widely variable chla (9 to 77.4 mg/m³) and total suspended solids (7-65 mg/L dry wt). Colored dissolved organic matter (CDOM) absorption at 440 nm was 0.20 to 2.50 m^− 1; Secchi disk transparency ranged from 0.28 to 1.5 m. The two- and three-band models were spectrally tuned to select the spectral bands for most accurate chla estimation. Strong linear relationships were established between analytically measured chla and both the three-band model [R^− 1(675)-R^− 1(695)] × R(730) and the two-band model R(720)/R(670), where R(λ) is reflectance at wavelength λ. The three-band model accounted for 81% of variation in chla and allowed estimation of chla with a root mean square error (RMSE) of less than 7.9 mg/m³, whereas the two-band model accounted for 79% of chla variability and RMSE of chla estimation was below 8.4 mg/m³. The three-band model with MERIS spectral bands allows accurate chla estimation with RMSE below 9.1 mg/m³. Two-band model with SeaWiFS bands and MODIS 667 nm and 748 nm bands can estimate chla with RMSE below 11 mg/m³. The findings underlined the rationale behind the conceptual model and demonstrated the robustness of this algorithm for chla retrieval in turbid, productive estuarine waters.     

New models for retrieving and partitioning the colored dissolved organic matter in the global ocean: Implications for remote sensing

Despite the importance of CDOM to upper ocean biogeochemical processes and optics, our current understanding of its spatial and temporal distributions and the factors controlling these distributions is very limited. This eventually prevents an understanding of its relationship to the pool of dissolved organic carbon in coastal and open oceans. This work aims to present a new approach for accurate modeling of absorption spectra of CDOM (a_cdom) and deriving information on its composition in global ocean waters. The modeling approach uses measurements (in situ) of the remote sensing reflectances at two wavelengths (denoted ₄₄₃⁵⁵⁵R_rs) to estimate a_cdom(350) and a_cdom(412), applies them to determine two spectral slopes of an exponential curve fit (S) and a hyperbolic curve fit (γ), derives an appropriate parameter (γ^o) for grading the CDOM compositional changes from a_cdom (350) and γ, and finally employs a_cdom(350), S, and γ^o in a modified exponential model to describe a_cdom(λ) as a function of wavelength. The robustness of this model was rigorously tested on three independent datasets, such as NOMAD in situ data, NOMAD SeaWiFS match-ups data and IOCCG simulated data (all of them contain a_cdom(λ) and R_rs(λ)), which represent a variety of waters within coastal and offshore regions around the world. Accuracy of the retrievals found with the new models was generally excellent, with MRE (mean relative error) and RMSE (root mean square error) of − 5.64-3.55% and 0.203-0.318 for the NOMAD in situ datasets, and − 5.63 to −0.98% and 0.136-0.241 for the NOMAD satellite datasets respectively (for λ₄₁₂ to λ₆₇₀). When used with SeaWiFS images collected over the regional and global waters, the new model showed the highest surface abundances of CDOM within the subpolar gyres and continental shelves dominated by terrestrial inputs (and perhaps local production) of colored dissolved materials, and the lowest surface abundances of CDOM in the central subtropical gyres and the open oceans presumably regulated by photobleaching phenomenon, bacterial activity and local processes. Significant interseasonal and interannual seasonal changes in the terrestrially-derived CDOM distributions were noticed from these new products that closely corresponded with the global mean runoff/river discharge induced by climate change/warming scenarios.     

Neural network estimation of chlorophyll a from MERIS full resolution data for the coastal waters of Galician rias (NW Spain)

In typical Case 2 waters, accurate remote sensing retrieval of chlorophyll a (chla) is still a challenging task. In this study, focusing on the Galician rias (ΝW Spain), algorithms based on neural network (NN) techniques were developed for the retrieval of chla concentration in optically complex waters, using Medium Resolution Imaging Spectrometer (MERIS) data. There is considerable interest in the accurate estimation of chla for the Galician rias, because of the economic and social importance of the extensive culture of mussels, and the high frequency of harmful algal events. Fifteen MERIS full resolution (FR) cloud-free images paired with in situ chla data (for 2002-2004 and 2006-2008) were used for the development and validation of the NN. The scope of NN was established from the clusters obtained using fuzzy c-mean (FCM) clustering techniques applied to the satellite-derived data. Three different NNs were developed: one including the whole data set, and two others using only points belonging to one of the clusters. The input data for these latter two NNs was chosen depending on the quality level, defined on the basis of quality flags given to each data set. The fitting results were fairly good and proved the capability of the tool to predict chla concentrations in the study area. The best prediction was given for the NN trained with high-quality data using the most abundant cluster data set. The performance parameters in the validation set of this NN were R² = 0.86, mean percentage error (MPE) = − 0.14, root mean square error (RMSE) = 0.75 mg m^− 3, and relative RMSE = 66%. The NN developed in this study detected accurately the peaks of chla, in both training and validation sets. The performance of the Case-2-Regional (C2R) algorithm, routinely used for MERIS data, was also tested and compared with our best performing NN and the sea-truthing data. Results showed that this NN outperformed the C2R, giving much higher R² and lower RMSE values.This study showed that the combination of in situ data and NN technology improved the retrieval of chla in Case 2 waters, and could be used to obtain more accurate chla maps. A local-based algorithm for the chla retrieval from an ocean colour sensor with the characteristics of MERIS would be a great support in the quantitative monitoring and study of harmful algal events in the coastal waters of the Rias Baixas. The limitations and possible improvements of the developed chla algorithms are also discussed.     

Validation of SeaWiFS chlorophyll a concentrations in the Southern Ocean: A revisit

Surface chlorophyll a concentrations (C_a, mg m^− 3) in the Southern Ocean estimated from SeaWiFS satellite data have been reported in the literature to be significantly lower than those measured from in situ water samples using fluorometric methods. However, we found that high-resolution (∼ 1 km²/pixel) daily SeaWiFS C_a (C_a^SWF) data (SeaDAS4.8, OC4v4 algorithm) was an accurate measure of in situ C_a during January-February of 1998-2002 if concurrent in situ data measured by HPLC (C_a^HPLC) instead of fluorometric (C_a^Fluor) measurements were used as ground truth. Our analyses indicate that C_a^Fluor is 2.48 ± 2.23 (n = 647) times greater than C_a^HPLC between 0.05 and 1.5 mg m^− 3 and that the percentage overestimation of in situ C_a by fluorometric measurements increases with decreasing concentrations. The ratio of C_a^SWF/C_a^HPLC is 1.12 ± 0.91 (n = 96), whereas the ratio of C_a^SWF/C_a^Fluor is 0.55 ± 0.63 (n = 307). Furthermore, there is no significant bias in C_a^SWF (12% and − 0.07 in linear and log-transformed C_a, respectively) when C_a^HPLC is used as ground truth instead of C_a^Fluor. The high C_a^Fluor/C_a^HPLC ratio may be attributed to the relatively low concentrations of chlorophyll b (C_b/C_a = 0.023 ± 0.034, n = 482) and relatively high concentrations of chlorophyll c (C_c/C_a = 0.25 ± 0.59, n = 482) in the phytoplankton pigment composition when compared to values from other regions. Because more than 90% of the waters in the study area, as well as in the entire Southern Ocean (south of 60° S), have C_a^SWF between 0.05 and 1.5 mg m^− 3, we consider that the SeaWiFS performance of C_a retrieval is satisfactory and for this C_a range there is no need to further develop a “regional” bio-optical algorithm to account for the previous SeaWiFS “underestimation”.     

Application of a generalized additive model (GAM) for estimating chlorophyll-a concentration from MODIS data in the Bohai and Yellow Seas,China

In optically complex waters, it is important to evaluate the accuracy of the standard satellite chlorophyll-a (chl-a) concentration algorithms, and to develop accurate algorithms for monitoring the dynamics of chl-a concentration. In this study, the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite remote-sensing reflectance and concurrent in situ measured chl-a (2010–2013) were used to evaluate the standard OC3M algorithm (ocean chlorophyll-a three-band algorithm for MODIS) and Graver–Siegel–Maritorena model version 1 (GSM01) algorithm for estimating chl-a concentration in the Bohai and Yellow Seas (BYS). The results showed that the chl-a algorithms of OC3M and GSM01 with global default parameters presented poor performance in the BYS (the mean absolute percentage difference (MAPD) and coefficient of determination (R²) of OC3M are 222.27% and 0.25, respectively; the MAPD and R² of GSM01 are 118.08% and 0.07, respectively). A novel statistical algorithm based on the generalized additive model (GAM) was developed, with the aim of improving the satellite-derived chl-a accuracy. The GAM algorithm was established using the in situ measured chl-a concentration as the output variable, and the MODIS above water remote-sensing reflectance (visible bands at 412, 443, 469, 488, 531, 547, 555, 645, 667, and 678 nm) and bathymetry (water depth) as input variables. The MAPD and R² calculated between the GAM and the in situ chl-a concentration are 39.96% and 0.67, respectively. The results suggest that the GAM algorithm can yield a superior performance in deriving chl-a concentrations relative to the standard OC3M and GSM01 algorithms in the BYS.     

Evaluation of SeaWiFS chlorophyll algorithms in the Southwestern Atlantic and Southern Oceans

Bio-optical measurements of spectral upwelling radiance and surface chlorophyll-a concentration have been conducted during 15 cruises between 1995 and 2004. The bio-optical data were divided into two sub-sets: the Southwestern Atlantic Ocean (SwAO), comprising a variety of biogeochemical provinces, from the oligotrophic waters in the South Atlantic gyre to the coastal waters influenced by La Plata River and Patos Lagoon discharge, and the Southern Ocean (SO) data set, comprising sampling stations south of the mean position of the Polar Front, with most stations being located in the vicinity of the Antarctic Peninsula. We derived regional chlorophyll algorithms for both regions and comparisons were made with the NASA's OC4v4 (operational algorithm) and OC2v4. For the Southwestern Atlantic region, the NASA OC4v4 algorithm presented a reasonable performance (r²=0.87, rmse-L=0.475, N=136) as compared to the revised algorithm for SwAO data (r²=0.89, rmse-L=0.426, N=136). A few stations under strong river plume influence were not considered in the analyses. These were detected by a higher reflectance at 670 nm, at low in situ chlorophyll concentration (<2 mg m⁻³). These results show that empirical algorithms applied to in-situ radiance data have a limited ability to extract accurate chlorophyll estimates below a 30% uncertainty level. For Southern Ocean stations, a 2-band linear-type model was generated (r²=0.64, rmse-L=0.347, N=77), which significantly improved the bias (6.4%) as compared to NASA's OC4v4 algorithm (bias=−21.7%). An evaluation of some published high-latitude algorithms on our data set has shown a better performance by taxon-specific models, even from distant regions. A validation experiment of the normalized spectral water-leaving radiances and chlorophyll-a SeaWiFS products was also conducted using the FURG-SwAO/SO data set, through a match-up exercise. Despite the relatively low number of pairs of radiometric measurements, SeaWiFS estimations compare well with in situ data (0.77<r²<0.98, N=21), although the satellite estimate show a marked bias (−35.6%) in the blue band nL_w (412). Regarding the chlorophyll-a concentration, an overall agreement was observed (r²=0.77, rmse-L=0.66, N=28), with a mean absolute percentage difference of 66%, which is above the goal generally accepted of 35% for satellite ocean color chlorophyll estimates. For the studied Southern Ocean area (mainly the Bransfield Strait), NASA's OC4v4 algorithm systematically underestimates chlorophyll above 0.2 mg m⁻³, as previously demonstrated by other researchers.     

Detecting the red tide algal blooms from satellite ocean color observations in optically complex Northeast-Asia Coastal waters

Over the last few decades, the coastal regions throughout the world have experienced incidences of algal blooms, which are harmful or otherwise toxic because of their potential threat to humans as well as marine organisms, owing to accelerated eutrophication from human activities and certain oceanic processes. Previous studies have found that correct identification of these blooms remains a great challenge with the standard bio-optical algorithms applied to satellite ocean color data in optically complex coastal waters containing high concentrations of the interfered dissolved organic and particulate inorganic materials. Here a new method called the red tide index (RI) is presented which is capable of identifying potential areas of harmful algal blooms (HABs) from SeaWiFS ocean color measurements representing the typical Case-2 water environments off the Korean and Chinese coasts. The RI method employs the water-leaving radiances (L_w), collected from in-situ radiometric measurements of three SeaWiFS bands centered at 443 nm, 510 nm and 555 nm, to achieve derivation of indices that are then related to absorbing characteristics of harmful algae (i.e., L_w at 443 nm) from which a best fit with a cubic polynomial function with correlation coefficient of R² = 0.91 is obtained providing indices of higher ranges for HABs and lower and slightly reduced ranges for turbid and non-bloom waters. Similar indices derived from the use of remote sensing reflectance (R_rs), normalized water-leaving radiance (nL_w) and combination of both are found rather inadequate to characterize the variability of the encountered bloom. In order to quantify the HABs in terms of chlorophyll (Chl), an empirical relationship is established between the RI and in-situ Chl in surface waters from about 0.4-71 mg m^− 3, which yields a Red tide index Chlorophyll Algorithm (RCA) based on an exponential function with correlation coefficient R² = 0.92. The established methods were extensively tested and compared with the performances of standard Ocean Chlorophyll 4 (OC4) algorithm and Local Chlorophyll Algorithm (LCA) using SeaWiFS images collected from typical red tide waters of Korean South Sea (KSS), East China Sea (ECS), Yellow Sea (YS) and Bohai Sea (BS) during 1999-2002. The standard spectral ratio algorithms, the OC4 and LCA, yielded large errors in Chl retrievals for coastal areas, besides providing false information about the encountered HABs in KSS, ECS, YS and BS waters. On the contrary, the RI coupled with the standard spectral ratios yielded comprehensive information about various ranges of algal blooms, while RCA Chl showing a good agreement with in-situ data led to enhanced understanding of the spatial and temporal variability of the recent HAB occurrences in high scattering and absorbing waters off the Korean and Chinese coasts.     

Estimating chlorophyll a concentrations from remote-sensing reflectance in optically shallow waters

A multi-spectral classification and quantification technique is developed for estimating chlorophyll a concentrations, Chl, in shallow oceanic waters where light reflected by the bottom can contribute significantly to the above-water remote-sensing reflectance spectra, R_rs(λ). Classification criteria for determining bottom reflectance contributions for shipboard R_rs(λ) data from the west Florida shelf and Bahamian waters (1998-2001; n = 451) were established using the relationship between R_rs(412)/R_rs(670) and the spectral curvature about 555 nm, [R_rs(412) ? R_rs(670)]/R_rs(555)². Chlorophyll concentrations for data classified as “optically deep” and “optically shallow” were derived separately using best-fit cubic polynomial functions developed from the band-ratios R_rs(490)/R_rs(555) and R_rs(412)/R_rs(670), respectively. Concentrations for transitional data were calculated from weighted averages of the two derived values. The root-mean-square error (RMSE_log10) calculated for the entire data set using the new technique was 14% lower than the lowest error derived using the best individual band-ratio. The standard blue-to-green, band-ratio algorithm yields a 26% higher RMSE_log10 than that calculated using the new method. This study demonstrates the potential of quantifying chlorophyll a concentrations more accurately from multi-spectral satellite ocean color data in oceanic regions containing optically shallow waters.     

Assessing structural effects on PRI for stress detection in conifer forests

The retrieval of indicators of vegetation stress from remote sensing imagery is an important issue for the accurate assessment of forest decline. The Photochemical Reflectance Index (PRI) has been demonstrated as a physiological index sensitive to the epoxidation state of the xanthophyll cycle pigments and to photosynthetic efficiency, serving as a proxy for short-term changes in photosynthetic activity, stress condition, and pigment absorption, but highly affected by illumination conditions, viewing geometry and canopy structure. In this study, a diurnal airborne campaign was conducted over Pinus sylvestris and Pinus nigra forest areas with the Airborne Hyperspectral Scanner (AHS) to evaluate the effects of canopy structure on PRI when used as an indicator of stress in a conifer forest. The AHS airborne sensor was flown at two times (8:00 GMT and 12:00 GMT) over forest areas under varying field-measured stress levels, acquiring 2 m spatial resolution imagery in 80 spectral bands in the 0.43-12.5 μm spectral range. Five formulations of PRI (based on R₅₃₁ as a xanthophyll-sensitive spectral band) were calculated using different reference wavelengths, such as PRI₅₇₀ (reference band R_REF = R₅₇₀), and the PRI modifications PRI_m1 (R_REF = R₅₁₂), PRI_m2 (R_REF = R₆₀₀), PRI_m3 (R_REF = R₆₇₀), and PRI_m4 (R_REF = R₅₇₀, R₆₇₀), along with other structural indices such as NDVI, SR, OSAVI, MSAVI and MTVI2. In addition, thermal bands were used for the retrieval of the land surface temperature. A radiative transfer modeling method was conducted using the LIBERTY and INFORM models to assess the structural effects on the PRI formulations proposed, studying the sensitivity of PRI_m indices to detect stress levels while minimizing the effects caused by the conifer architecture. The PRI indices were related to stomatal conductance, xanthophyll epoxidation state (EPS) and crown temperature. The modeling analysis showed that the coefficient of variation (CV) for PRI was 50%, whereas the CV for PRIm1 (band R512 as a reference) was only 20%. Simulation and experimental results demonstrated that PRI_m1 (R_REF = R₅₁₂) was less sensitive than PRI (R_REF = R₅₇₀) to changes in Leaf Area Index (LAI) and tree densities. PRI₅₁₂ was demonstrated to be sensitive to EPS at both leaf (r² = 0.59) and canopy level (r² = 0.40), yielding superior performance than PRI₅₇₀ (r² = 0.21) at the canopy level. In addition, PRI₅₁₂ was significantly related to water stress indicators such as stomatal conductance (Gs; r² = 0.45) and water potential (Ψ; r² = 0.48), yielding better results than PRI₅₇₀ (Gs, r² = 0.21; Ψ, r² = 0.21) due to the structural effects found on the PRI₅₇₀ index at the canopy level.     

Monitoring turbidity in Tampa Bay using MODIS/Aqua 250-m imagery

We developed an approach to map turbidity in estuaries using a time series (May 2003 to April 2006) of 250-m resolution images from the Moderate Resolution Imaging Spectroradiometer (MODIS) on NASA's Aqua satellite, using Tampa Bay as a case study. Cross-calibration of the MODIS 250-m data (originally designed for land use) with the well-calibrated MODIS 1-km ocean data showed that the pre-launch radiometric calibration of the 250-m bands was adequate. A simple single scattering atmospheric correction provided reliable retrievals of remote sensing reflectance at 645 nm (0.002 < R_rs(645) < 0.015 sr^− 1, median bias = − 7%, slope = 0.95, intercept = 0.00, r² = 0.97, n = 15). A more rigorous approach, using a multiple scattering atmospheric correction of the cross-calibrated at-sensor radiances, retrieved similar R_rs(645). R_rs(645) estimates, after stringent data quality control, showed a close correlation with in situ turbidity (turbidity = 1203.9 × R_rs(645)^1.087, 0.9 < turbidity < 8.0 NTU, r² = 0.73, n = 43). MODIS turbidity imagery derived using the developed approach showed that turbidity in Hillsborough Bay (HB) was consistently higher than that in other sub-regions except in August and September, when higher concentrations of colored dissolved organic matter seem to have caused underestimates of turbidity. In comparison, turbidity in Middle Tampa Bay (MTB) was generally lowest among the Bay throughout the year. Both Old Tampa Bay (OTB) and Low Tampa Bay (LTB) showed marked seasonal variations with higher turbidity in LTB during the dry season and in OTB during the wet season, respectively. This seasonality is linked to wind-driven bottom resuspension events in lower portion of the Bay and river inputs of sediments in the upper portion of the Bay. The Bay also experiences significant interannual variation in turbidity, which was attributed primarily to changes in wind forcing. Compared with the once-per-month, non-synoptic in situ surveys, synoptic and frequent sampling facilitated by satellite remote sensing provides improved assessments of turbidity patterns and thus a valuable tool for operational monitoring of water quality of estuarine and coastal waters such as in Tampa Bay.     

Can a satellite-derived estimate of the fraction of PAR absorbed by chlorophyll (FAPARchl) improve predictions of light-use efficiency and ecosystem photosynthesis for a boreal aspen forest?

We used daily MODerate resolution Imaging Spectroradiometer (MODIS) imagery obtained over a five-year period to analyze the seasonal and inter-annual variability of the fraction of absorbed photosynthetically active radiation (FAPAR) and photosynthetic light use efficiency (LUE) for the Southern Old Aspen (SOA) flux tower site located near the southern limit of the boreal forest in Saskatchewan, Canada. To obtain the spectral characteristics of a standardized land area to compare with tower measurements, we scaled up the nominal 500 m MODIS products to a 2.5 km × 2.5 km area (5 × 5 MODIS 500 m grid cells). We then used the scaled-up MODIS products in a coupled canopy-leaf radiative transfer model, PROSAIL-2, to estimate the fraction of absorbed photosynthetically active radiation (APAR) by the part of the canopy dominated by chlorophyll (FAPAR_chl) versus that by the whole canopy (FAPAR_canopy). Using the additional information provided by flux tower-based measurements of gross ecosystem production (GEP) and incident PAR, we determined 90-minute averages for APAR and LUE (slope of GEP:APAR) for both the physiologically active foliage (APAR_chl, LUE_chl) and for the entire canopy (APAR_canopy, LUE_canopy).The flux tower measurements of GEP were strongly related to the MODIS-derived estimates of APAR_chl (r² = 0.78) but only weakly related to APAR_canopy (r² = 0.33). Gross LUE between 2001 and 2005 for LUE_chl was 0.0241 µmol C µmol^− 1 PPFD whereas LUE_canopy was 36% lower. Time series of the 5-year normalized difference vegetation index (NDVI) were used to estimate the average length of the core growing season as days of year 152-259. Inter-annual variability in the core growing season LUE_chl (µmol C µmol^− 1 PPFD) ranged from 0.0225 in 2003 to 0.0310 in 2004. The five-year time series of LUE_chl corresponded well with both the seasonal phase and amplitude of LUE from the tower measurements but this was not the case for LUE_canopy. We conclude that LUE_chl derived from MODIS observations could provide a more physiologically realistic parameter than the more commonly used LUE_canopy as an input to large-scale photosynthesis models.     

HYDROPT: A fast and flexible method to retrieve chlorophyll-a from multispectral satellite observations of optically complex coastal waters

We present a generic innovative algorithm for remote sensing of coastal waters that can deal with a large range of concentrations of chlorophyll-a, SPM and CDOM and their inherent optical properties. The algorithm is based on the exact solutions of the HYDROLIGHT numerical radiative transfer model to support retrieval in optically complex waters with varying sensor wide swath viewing geometry. The algorithm estimates the concentrations by minimizing the difference between observed and modeled reflectance spectra. The use of a look-up table and polynomial interpolation greatly reduces computation time, allowing operational and near-real time processing of large sets of satellite imagery. Because the remote sensing reflectance was tabulated as a function of in-water light absorption and scattering, rather than actual constituents concentrations, the algorithm can be applied with any definition of the specific inherent optical properties of CHL, SPM and CDOM. A statistical measure for the goodness-of-fit and the formal standard errors in the fitted concentrations are provided, thus producing error maps with each thematic chlorophyll image, often lacking in most applications of innovative algorithms. The performance of the algorithm is demonstrated for multispectral observations of the North Sea, a shallow coastal sea with large concentration gradients in SPM (due to resuspension) and CDOM (from riverine influx). The standard errors of estimated chlorophyll-a concentrations ranged between 0.5 and 3 (mg m^− 3) for mean concentrations between 2 and 20 (mg m^− 3), quite acceptable results for these optically complex waters.     

Impact of sea ice on the retrieval of water-leaving reflectance, chlorophyll a concentration and inherent optical properties from satellite ocean color data

Two physical phenomena by which satellite remotely sensed ocean color data are contaminated by sea ice at high latitudes are described through simulations and observations: (1) the adjacency effect that occurs along sea ice margins and (2) the sub-pixel contamination by a small amount of sea ice within an ocean pixel. The signal at the top of the atmosphere (TOA) was simulated using the 6S radiative transfer code that allows modeling of the adjacency effect for various types of sea ice surrounding an open water area. In situ sea ice reflectance spectra used in the simulations were measured prior to and during the melt period as part of the 2004 Canadian Arctic Shelf Exchange Study (CASES). For sub-pixel contamination, the TOA signal was simulated for various surface reflectances obtained by linear mixture of both sea ice and water-leaving reflectances (ρ_w). The standard atmospheric correction algorithm was then applied to the simulated TOA spectra to retrieve ρ_w spectra from which chlorophyll a concentrations (CHL) and inherent optical properties (IOPs) were derived. The adjacency effect was associated with large errors (> 0.002) in the retrieval of ρ_w as far as 24 km from an ice edge in the blue part of the spectrum (443 nm). Therefore, for moderate to high CHL (> 0.5 mg m^− 3), any pixel located within a distance of ∼ 10-20 km from the ice edge were unreliable. It was also found necessary to consider the adjacency effect when the total absorption coefficient (a_t) was to be retrieved using a semi-analytical algorithm. a_t(443) was underestimated by more than 35% at a distance of 20 km from an ice edge for CHL > 0.5 mg m^− 3. The effect on the retrieval of the particle backscattering coefficient (b_bp) was important only for clear waters (CHL ∼ 0.05 mg m^− 3). In contrast, sub-pixel contamination by a small amount of sea ice produced systematic underestimation of ρ_w in the blue because of incorrect interpretation of enhanced reflectance in the near infrared that is attributed to higher concentrations of atmospheric aerosols. In general, sub-pixel contamination was found to result in overestimations of CHL and a_t, and underestimations of b_bp. A simple method was proposed to flag pixels contaminated by adjacency effect.     

Hyperspectral remote sensing of cyanobacterial pigments as indicators for cell populations and toxins in eutrophic lakes

The growth of mass populations of toxin-producing cyanobacteria is a serious concern for the ecological status of inland waterbodies and for human and animal health. In this study we examined the performance of four semi-analytical algorithms for the retrieval of chlorophyll a (Chl a) and phycocyanin (C-PC) from data acquired by the Compact Airborne Spectrographic Imager-2 (CASI-2) and the Airborne Imaging Spectrometer for Applications (AISA) Eagle sensor. The retrieval accuracies of the semi-analytical models were compared to those returned by optimally calibrated empirical band-ratio algorithms. The best-performing algorithm for the retrieval of Chl a was an empirical band-ratio model based on a quadratic function of the ratio of reflectance at 710 and 670 nm (R² = 0.832; RMSE = 29.8%). However, this model only provided a marginally better retrieval than the best semi-analytical algorithm. The best-performing model for the retrieval of C-PC was a semi-analytical nested band-ratio model (R² = 0.984; RMSE = 3.98 mg m⁻³). The concentrations of C-PC retrieved using the semi-analytical model were correlated with cyanobacterial cell numbers (R² = 0.380) and the particulate and total (particulate plus dissolved) pools of microcystins (R² = 0.858 and 0.896 respectively). Importantly, both the empirical and semi-analytical algorithms were able to retrieve the concentration of C-PC at cyanobacterial cell concentrations below current warning thresholds for cyanobacteria in waterbodies. This demonstrates the potential of remote sensing to contribute to early-warning detection and monitoring of cyanobacterial blooms for human health protection at regional and global scales.     

A simple semi-analytical model for remote estimation of chlorophyll-a in turbid waters: Validation

Accurate assessment of phytoplankton chlorophyll-a (chla) concentrations in turbid waters by means of remote sensing is challenging due to the optical complexity of case 2 waters. We have applied a recently developed model of the form [R_rs^? 1(λ₁) ? R_rs^? 1(λ₂)] × R_rs(λ₃) where R_rs(λ_i) is the remote-sensing reflectance at the wavelength λ_i, for the estimation of chla concentrations in turbid waters. The objectives of this paper are (a) to validate the three-band model as well as its special case, the two-band model R_rs^? 1(λ₁) × R_rs(λ₃), using datasets collected over a considerable range of optical properties, trophic status, and geographical locations in turbid lakes, reservoirs, estuaries, and coastal waters, and (b) to evaluate the extent to which the three-band model could be applied to the Medium Resolution Imaging Spectrometer (MERIS) and two-band model could be applied to the Moderate Resolution Imaging Spectroradiometer (MODIS) to estimate chla in turbid waters.The three-band model was calibrated and validated using three MERIS spectral bands (660–670 nm, 703.75–713.75 nm, and 750?757.5 nm), and the 2-band model was tested using two MODIS spectral bands (λ₁ = 662–672, λ₃ = 743–753 nm). We assessed the accuracy of chla prediction in four independent datasets without re-parameterization (adjustment of the coefficients) after initial calibration elsewhere. Although the validation data set contained widely variable chla (1.2 to 236 mg m^? 3), Secchi disk depth (0.18 to 4.1 m), and turbidity (1.3 to 78 NTU), chla predicted by the three-band algorithm was strongly correlated with observed chla (r² > 0.96), with a precision of 32% and average bias across data sets of ? 4.9% to 11%. Chla predicted by the two-band algorithm was also closely correlated with observed chla (r² > 0.92); however, the precision declined to 57%, and average bias across the data sets was 18% to 50.3%. These findings imply that, provided that an atmospheric correction scheme for the red and NIR bands is available, the extensive database of MERIS and MODIS imagery could be used for quantitative monitoring of chla in turbid waters.     

Bio-optical properties in waters influenced by the Mississippi River during low flow conditions

Spatial and temporal patterns of bio-optical properties were studied in the Northern Gulf of Mexico during cruises in April and October of 2000, a year during which the discharge volume from the Mississippi River was unusually low. Highly variable surface Chl a concentrations (0.1 to 17 mg m⁻³) and colored dissolved organic matter (CDOM) absorption (0.07 to 1.1 m⁻¹ at 412 nm) were observed in the study region that generally decreased with increasing salinity waters, being highest nearshore and decreasing at offshore stations. The optical properties of absorption, scattering, and diffuse attenuation coefficients reflected these distributions with phytoplankton particles and CDOM contributing to most of the spatial, vertical, and seasonal variability. The diffuse attenuation coefficient K_d(λ) and spectral remote sensing reflectance R_rs(λ) were linear functions of absorption and backscattering coefficients a(λ) and b_b(λ) through the downward average cosine μ_d and the ratio of variables f/Q at the SeaWiFS wavebands for waters with widely varying bio-optical conditions. Although various R_rs(λ) ratio combinations showed high correlation with surface Chl a concentrations and CDOM absorption at 412 nm, power law equations derived using the R_rs(490)/R_rs(555) and R_rs(510)/R_rs(555) ratios provided the best retrievals of Chl a concentrations and CDOM absorption from SeaWiFS reflectance data.     

Mapping lake CDOM by satellite remote sensing

Given the importance of coloured dissolved organic matter (CDOM) for the structure and function of lake ecosystems, a method to estimate the amount of CDOM in lake waters over large geographic areas would be highly desirable. Advanced Land Imager (ALI) images were acquired in southern Finland (in 2002) and southern Sweden (in 2003) together with in situ measurements of bio-optical properties of 34 lakes (39 measuring stations). Based on this dataset, a band-ratio type algorithm was developed using ALI band 2 and band 3 for estimating CDOM content (absorption of filtrated water at 420 nm) in lakes. Correlation between in situ measured CDOM and the remote sensing estimate of CDOM was high, r²=0.73. The CDOM retrieval algorithm obtained on the basis of two images and in situ data was validated on a third ALI image (eastern Finland, 2002) that was available in the ALI image archive. In situ water-colour monitoring data from 22 lakes (27 measuring stations) in the third image were available in a database of the Finnish Environment Administration. The water-colour data were converted to CDOM absorption values, which were then compared to the results from a third ALI image. The correlation between remotely estimated and in situ CDOM values in the algorithm validation image was high, r²=0.83. These results support the conclusion that CDOM content in lakes over a wide range of concentrations (a_CDOM(420) between 0.68 and 11.13 m⁻¹) can be mapped using Advanced Land Imager data.