Intercomparison of ocean colour band-ratio algorithms for chlorophyll concentration in the Subtropical Front east of New Zealand

This study intercompared the performance of eight band-ratio chlorophyll-a algorithms which together can be used to process measurements from the ocean colour satellite sensors CZCS, OCTS, SeaWiFS, MODIS, MERIS, and GLI. The study area included Subtropical, Subtropical Front and Subantarctic waters east of New Zealand, and Case 1 waters of the New Zealand northeast continental shelf. Over 170 co-incident measurements of spectral normalised water-leaving radiance and near-surface concentration of chlorophyll-a were made on nine research voyages between 1998 and 2000. The studentised bootstrap method was used to identify statistically significant bias in algorithm products relative to in situ measurements. The band-ratio algorithms used by CZCS, OCTS and SeaWiFS missions systematically underestimated chlorophyll-a concentration in the offshore regions by between 21% and 45%, but showed no systematic bias in the continental shelf waters. The band-ratio algorithms applicable to the MODIS and MERIS sensors had no clear bias with respect to in situ measurements in offshore waters, but had a positive bias of 20% over the continental shelf. The proposed GLI band-ratio algorithm led to estimates that were negatively biased with respect to in situ measurement offshore (− 30%), and positively biased over the continental shelf (20%). The results were consistent with unusually high values of absorption in the blue part of the spectrum (443-490 nm) compared to the green part (∼ 550 nm) by phytoplankton pigments in the offshore waters, and high chlorophyll-specific absorption over the continental shelf.     

An assessment of the black ocean pixel assumption for MODIS SWIR bands

Recent studies show that an atmospheric correction algorithm using shortwave infrared (SWIR) bands improves satellite-derived ocean color products in turbid coastal waters. In this paper, the black pixel assumption (i.e., zero water-leaving radiance contribution) over the ocean for the Moderate Resolution Imaging Spectroradiometer (MODIS) SWIR bands at 1240, 1640, and 2130 nm is assessed for various coastal ocean regions. The black pixel assumption is found to be generally valid with the MODIS SWIR bands at 1640 and 2130 nm even for extremely turbid waters. For the MODIS 1240 nm band, however, ocean radiance contribution is generally negligible in mildly turbid waters such as regions along the U.S. east coast, while some slight radiance contributions are observed in extremely turbid waters, e.g., some regions along the China east coast, the estuary of the La Plata River. Particularly, in the Hangzhou Bay, the ocean radiance contribution at the SWIR band 1240 nm results in an overcorrection of atmospheric and surface effects, leading to errors of MODIS-derived normalized water-leaving radiance at the blue reaching ~ 0.5 mW cm^− 2 μm^− 1 sr^− 1. In addition, we found that, for non-extremely turbid waters, i.e., the ocean contribution at the near-infrared (NIR) band < ~ 1.0 mW cm^− 2 μm^− 1 sr^− 1, there exists a good relationship in the regional normalized water-leaving radiances between the red and the NIR bands. Thus, for non-extremely turbid waters, such a red-NIR radiance relationship derived regionally can possibly be used for making corrections for the regional NIR ocean contributions without using the SWIR bands, e.g., for atmospheric correction of ocean color products derived from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS).     

Validation of satellite ocean color primary products at optically complex coastal sites: Northern Adriatic Sea, Northern Baltic Proper and Gulf of Finland

The study presents and discusses the application of in situ data from the ocean color component of the Aerosol Robotic Network (AERONET-OC) to assess primary remote sensing products from the Moderate Resolution Imaging Spectroradiometer (MODIS) on the AQUA platform and from the Sea-viewing Wide-Field-of-view Sensor (SeaWiFS) on the OrbView-2 spacecraft. Three AERONET-OC European coastal sites exhibiting different atmospheric and marine optical properties were considered for the study: the Acqua Alta Oceanographic Tower (AAOT) in the northern Adriatic Sea representing Case-1 and Case-2 moderately sediment dominated waters; and, the Gustaf Dalen Lighthouse Tower (GDLT) in the northern Baltic Proper and the Helsinki Lighthouse Tower (HLT) in the Gulf of Finland, both characterized by Case-2 waters dominated by colored dissolved organic matter (CDOM). The analysis of MODIS derived normalized water-leaving radiance at 551 nm, L_WN(551), has shown relatively good results for all sites with uncertainties of the order of 10% and biases ranging from − 1 to − 4%. Larger uncertainty and bias have been observed at 443 nm for the AAOT (i.e., 18 and − 7%, respectively). At the same center wavelength, results for GDLT and HLT have exhibited much larger uncertainties (i.e., 56 and 67%, respectively) and biases (i.e., 18 and 25%, respectively), which undermine the possibility of presently using remote sensing L_WN data at the blue center wavelengths for bio-optical investigations in the Baltic Sea. An evaluation of satellite derived aerosol optical thickness, τ_a, has shown uncertainties and biases of the order of tens of percent increasing with wavelength at all sites. Specifically, MODIS derived τ_a at 869 nm has shown an overestimate of 71% at the AAOT, 101% at GDLT and 91% at HLT, respectively. This result highlights the effects of a limited number of aerosol models for the atmospheric correction process, and might also indicate the need of applying a vicarious calibration factor to the remote sensing data at the 869 nm center wavelength to remove the effects of uncertainties in the atmospheric optical model and the space sensor radiometric calibration. Similar results have been obtained from the analysis of SeaWiFS data. Finally, in view of illustrating the possibility of increasing the accuracy of satellite regional radiometric products, AERONET-OC data have been applied to reduce systematic errors in MODIS and Medium Resolution Imaging Spectrometer (MERIS) L_WN data likely due to the atmospheric correction process. Results relying on MODIS match-ups for the Baltic Sites (i.e., GDLT and HLT) and MERIS matchups for the AAOT, have indicated a substantial reduction of both uncertainty and bias in the blue and red center wavelengths.     

Estimating the altitude of aerosol plumes over the ocean from reflectance ratio measurements in the O2 A-band

A methodology is proposed to infer the altitude of aerosol plumes over the ocean from reflectance ratio measurements in the O₂ absorption A-band (759 to 770 nm). The reflectance ratio is defined as the ratio of the reflectance in a first spectral band, strongly attenuated by O₂ absorption, and the reflectance in a second spectral band, minimally attenuated. For a given surface reflectance, simple relations are established between the reflectance ratio and the altitude of an aerosol layer, as a function of atmospheric conditions and the geometry of observation. The expected accuracy for various aerosol loadings and models is first quantified using an accurate, high spectral resolution, radiative transfer model that fully accounts for interactions between scattering and absorption. The method is developed for POLDER and MERIS, satellite sensors with adequate spectral characteristics. The simulations show that the method is only accurate over dark surfaces when aerosol optical thickness at 765 nm is relatively large (> 0.3). In this case, the expected accuracy is on the order of ± 0.5 km or ± 0.2 km for POLDER or MERIS respectively. More accurate estimates are obtained with MERIS, since in this case the spectral reflectance ratio is more sensitive to aerosol altitude. However, a precise spectral calibration is needed for MERIS. The methodology is applied to MERIS and POLDER imagery acquired over marine surfaces. The estimated aerosol altitude is compared with in situ lidar profiles of backscattering coefficient measured during the AOPEX-2004 experiment for MERIS, or obtained with the space-borne lidar CALIOP for POLDER. The retrieved altitudes agree with lidar measurements in a manner consistent with theory. These comparisons demonstrate the potential of the differential absorption methodology for obtaining information on aerosol altitude over dark surfaces.     

Sizing snow grains using backscattered solar light

In this article, we describe a technique to determine dry snow grain size from optical observations. The method is based on analysis of the snow reflectance in the near-infrared region, in particular, the Medium Resolution Imaging Spectrometer (MERIS) band at 865 nm, which is common to many spaceborne optical sensors, is used. In addition, the algorithm is applied to the Moderate Resolution Imaging Spectroradiometer (MODIS) 1240 nm band. It is found that bands located at 1020 and 1240 nm are the most suitable for snow grain size remote-sensing applications. The developed method is validated using MODIS observations over flat snow deposited on a lake ice in Hokkaido, Japan.     

A regional algorithm for separating light absorption by chlorophyll-a and coloured detrital matter in the Black Sea,using 480–560 nm bands from ocean colour scanners

The aim of this study is to modify the regional algorithm for Moderate Resolution Imaging Spectroradiometer (MODIS) and Medium-spectral Resolution Imaging Spectrometer (MERIS) bands using newly available data of seasonal and spatial variability of light absorption by all optically active components in the Black Sea, and to obtain a merged product based on data retrieved from all the colour scanners that have operated since September 1997. Comparison of chlorophyll-a concentration (chl-a) simulated by the standard National Aeronautics and Space Administration (NASA) algorithm with in situ chl-a measurements showed that the NASA algorithm provided incorrect chl-a assessment of Black Sea shelf and deep-sea waters during spring?summer. Originally the standard NASA algorithm could be applied if there was a high correlation between light absorption by phytoplankton (a_ph) and that by coloured dissolved and suspended organic matter (a_CDM), which is not the case in the Black Sea. Consequently, development of the correct regional chl-a algorithm requires splitting of light absorption into a_ph and a_CDM. This issue has been resolved by the proposed regional algorithm developed for the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) using remote-sensing reflectance in three (as minimum) spectral bands from 480 to 560 nm. Operation of the SeaWiFS and MERIS colour scanners ceased in December 2010 and April 2012, respectively, while the MODIS scanner is still working on the Terra and Aqua satellites. In this research, level 2 data (products of standard atmosphere correction at three bands filtered by masks/flags) of SeaWiFS, MODIS (on Terra and Aqua satellites), and MERIS scanners were retrieved for their mission lifetime. The regional algorithm was validated independently for each scanner, based on the adequacy of the algorithm-derived chl-a and a_CDM to in situ-measured data for the same day. The results suggest a satisfactory accuracy of the modified regional algorithm.     

AERONET-based surface reflectance validation network (ASRVN) data evaluation: Case study for railroad valley calibration site

The AERONET-based Surface Reflectance Validation Network (ASRVN) is an operational processing system developed for validation of satellite derived surface reflectance products at regional and global scales. The ASRVN receives 50 × 50 km² subsets of MODIS data centered at AERONET sites along with AERONET aerosol and water vapor data, and performs an atmospheric correction. The ASRVN produces surface bidirectional reflectance factor (BRF), albedo, parameters of the Ross-Thick Li-Sparse (RTLS) BRF model, as well as Hemispherical-Directional Reflectance Factor (HDRF), which is required for comparison with the ground-based measurements. This paper presents a comparison of ASRVN HDRF with the ground-based HDRF measurements collected during 2001-2008 over a bright calibration Railroad Valley, Nevada site as part of the MODIS land validation program. The ground measurements were conducted by the Remote Sensing Group (RSG) at the University of Arizona using an ASD spectrometer. The study reveals a good agreement between ASRVN and RSG HDRF for both MODIS Terra and Aqua with rmse ~ 0.01-0.025 in the 500 m MODIS land bands B1-B7. Obtained rmse is below uncertainties due to the spatial and seasonal variability of the bright calibration 1 km² area. While two MODIS instruments have a similar rmse in the visible bands, MODIS Aqua has a better agreement (lower rmse) with the ground data than MODIS Terra at wavelengths 0.87-2.1 μm. An independent overall good agreement of two MODIS instruments with the ground data indicates that the relative calibration of MODIS Terra and Aqua at medium-to-bright reflectance levels for the stated time period is significantly better than uncertainties of the ASRVN and ground data.     

Validation of MODIS derived aerosol optical depth over the Yangtze River Delta in China

MODIS derived aerosol optical depths (AODs) at 550 nm are compared with sunphotometer CE318 measurements at 7 sites located at Yangtze River Delta (YRD) in China from July to October, 2007. The evaluation result indicates that MODIS AODs (Collection 5, C005) are in good agreement with those from CE318 in dense vegetation regions, but show more differences in those regions with complex underlying surface (such as at lake water and urban surface sites). Reasons for these differences are discussed after removing cases with significant errors caused by validation scheme. The final validation result shows that MODIS AODs are in good agreement with CE318 with a correlation coefficient of 0.85 and RMS of 0.15. 90% of MODIS cases fall in the range of Δτ = ± 0.05 ± 0.20τ, indicating MODIS aerosol retrieval algorithm, aerosol models and surface reflectance estimate are generally suitably reasonable for aerosol retrieval in YRD. However, MODIS AODs show a systemic errors with fitted line of y = 0.75x + 0.13, indicating underestimation of AOD when aerosol loadings are high. Aerosol models and surface reflectance estimations are dominant sources of MODIS aerosol retrieval errors.     

Comparison of MERIS,MODIS, SeaWiFS-derived particulate organic carbon,and in situ measurements in the South China Sea

Medium Resolution Imaging Spectrometer (MERIS), Moderate Resolution Imaging Spectroradiometer (MODIS), and Sea-viewing Wide Field-of-view Sensor (SeaWiFS) particulate organic carbon (POC) concentration products for the South China Sea (SCS) were compared with in situ data collected from October 2007 to December 2013. Spectral remote-sensing reflectance (R_rs,λ) was also measured to help understand POC algorithm performance. A strict comparison of the satellite-derived POC and in situ measurements showed that MERIS, MODIS, and SeaWiFS underestimated in situ values by 29.1, 11.7, and 31.5%, respectively. Similar results were obtained with a relaxed matching criterion. Through analysis of the causes of product uncertainty, the results suggested that satellite retrieval of R_rs,λ and the global POC algorithm both have an impact on inversion accuracy. However, the formulation of the POC algorithm seems to be more critical. When a regional algorithm was developed to obtain satellite-derived POC, both the strict and relaxed comparison results showed significant improvement, but for coastal waters, both algorithms had larger errors. Other factors affecting the comparison are also discussed.     

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.     

Assessment of SeaWiFS,MODIS, and MERIS ocean colour products in the South China Sea

The Sea-viewing Wide Field-of-view Sensor (SeaWiFS), Moderate Resolution Imaging Spectroradiometer (MODIS), and Medium Resolution Imaging Spectrometer (MERIS) remote-sensing radiometric and chlorophyll-a (chl-a) concentration products for the South China Sea (SCS) from October 2003 to May 2010 were assessed using in situ data. A strict spatiotemporal match-up method was used to minimize the temporal variability effects of atmosphere and seawater around the measurement site. A comparison of the remote-sensing reflectance (R_rs(λ)) of the three sensors with in situ values from the open waters of the SCS showed that the mean absolute percentage difference varied from 13% to 55% in the 412–560 nm spectral range. Generally, the MERIS radiometric products exhibited higher typical uncertainties and bias than the SeaWiFS and MODIS products. The R_rs(443) to R_rs(555/551/560) band ratios of the satellite data were in good agreement with in situ observations for these sensors. The SeaWiFS, MODIS, and MERIS chl-a products overestimated in situ values by 74%, 42%, and 120%, respectively. MODIS retrieval accuracy was better than those of the other sensors, with MERIS performing the worst. When the match-up criteria were relaxed, the assessment results degraded systematically. Therefore, strict spatiotemporal match-up is recommended to minimize the possible influences of small-scale variation in geophysical properties around the measurement site. Coastal and open-sea areas in the SCS should be assessed separately because their biooptical properties are different and the results suggest different atmospheric correction problems.     

Comparison and evaluation of Medium Resolution Imaging Spectrometer leaf area index products across a range of land use

Leaf area index (LAI) is a commonly required parameter when modelling land surface fluxes. Satellite based imagers, such as the 300 m full resolution (FR) Medium Spectral Resolution Imaging Spectrometer (MERIS), offer the potential for timely LAI mapping. The availability of multiple MERIS LAI algorithms prompts the need for an evaluation of their performance, especially over a range of land use conditions. Four current methods for deriving LAI from MERIS FR data were compared to estimates from in-situ measurements over a 3 km × 3 km region near Ottawa, Canada. The LAI of deciduous dominant forest stands and corn, soybean and pasture fields was measured in-situ using digital hemispherical photography and processed using the CANEYE software. MERIS LAI estimates were derived using the MERIS Top of Atmosphere (TOA) algorithm, MERIS Top of Canopy (TOC) algorithm, the Canada Centre for Remote Sensing (CCRS) Empirical algorithm and the University of Toronto (UofT) GLOBCARBON algorithm. Results show that TOA and TOC LAI estimates were nearly identical (R² > 0.98) with underestimation of LAI when it is larger than 4 and overestimation when smaller than 2 over the study region. The UofT and CCRS LAI estimates had root mean square errors over 1.4 units with large (∼ 25%) relative residuals over forests and consistent underestimates over corn fields. Both algorithms were correlated (R² > 0.8) possibly due to their use of the same spectral bands derived vegetation index for retrieving LAI. LAI time series from TOA, TOC and CCRS algorithms showed smooth growth trajectories however similar errors were found when the values were compared with the in-situ LAI. In summary, none of the MERIS LAI algorithms currently meet performance requirements from the Global Climate Observing System.     

Estimation of cyanobacterial pigments in a freshwater lake using OCM satellite data

Cyanobacteria represent a major harmful algal group in fresh to brackish water environments. Lac des Allemands, a freshwater lake of 49 km² southwest of New Orleans, Louisiana on the upper end of the Barataria Estuary, provides a natural laboratory for remote characterization of cyanobacterial blooms because of their seasonal occurrence. The Oceansat-1 satellite Ocean Colour Monitor (OCM) provides measurements similar to SeaWiFS but with higher spatial resolution, and this work is the first attempt to use OCM measurements to quantify cyanobacterial pigments. The satellite signal was first vicariously calibrated using SeaWiFS as a reference, and then corrected to remove the atmospheric effects using a customized atmospheric correction procedure. Then, empirical inversion algorithms were developed to convert the OCM remote sensing reflectance (R_rs) at bands 4 and 5 (centered at 510.6 and 556.4 nm, respectively) to concentrations of phycocyanin (PC), the primary cyanobacterial pigment. A holistic approach was used to minimize the influence of other optically active constituents on the PC algorithm. Similarly, empirical algorithms to estimate chlorophyll a (Chl a) concentrations were developed using OCM bands 5 and 6 (centered at 556.4 and 669 nm, respectively). The best PC algorithm (R² = 0.7450, p < 0.0001, n = 72) yielded a root mean square error (RMSE) of 36.92 μg/L with a relative RMSE of 10.27% (PC from 2.75 to 363.50 μg/L, n = 48). The best algorithm for Chl a (R² = 0.7510, p < 0.0001, n = 72) produced an RMSE of 31.19 μg/L with a relative RMSE of 16.56% (Chl a from 9.46 to 212.76 μg/L, n = 48). While more field data are required to further validate the long-term performance of these algorithms, currently they represent the best protocol for establishing a long time-series of cyanobacterial blooms in the Lac des Allemands using OCM data.     

Algorithm for retrieval of aerosol optical properties over the ocean from the Geostationary Ocean Color Imager

An aerosol retrieval algorithm for the first Geostationary Ocean Color Imager (GOCI) to be launched in March 2010 onboard the Communication, Ocean, and Meteorological Satellite (COMS) is presented. The algorithm retrieves aerosol optical depth (AOD), fine-mode fraction (FMF), and aerosol type in 500 m × 500 m resolution. All the products are retrieved over clear water which is defined by surface reflectance ratio between 640 nm and 860 nm (SRR) less or equal to 2.5, while only AOD is retrieved over turbid water (SRR > 2.5) due to high surface reflectance. To develop optimized algorithm for the target area of GOCI, optical properties of aerosol are analyzed from extensive observation of AERONET sunphotometers to generate lookup table. Surface reflectance of turbid water is determined from 30-day composite of Rayleigh- and gas corrected reflectance. By applying the present algorithm to MODIS top-of-the atmosphere reflectance, three different aerosol cases dominated by anthropogenic aerosol contains black carbon (BC), dust, and non-absorbing aerosol are analyzed to test the algorithm. The algorithm retrieves AOD, and size information together with aerosol type which are consistent with results inferred by RGB image in a qualitative way. The comparison of the retrieved AOD with those of MODIS collection 5 and AERONET sunphotometer observations shows reliable results. Especially, the application of turbid water algorithm significantly increases the accuracy in retrieving AOD at Anmyon station. The sensitivity study between MODIS and GOCI instruments in terms of relative sensitivity and scattering angle shows promising applicability of the present algorithm to future GOCI measurements.     

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.     

Assessing the potential of SeaWiFS and MODIS for estimating chlorophyll concentration in turbid productive waters using red and near-infrared bands

Bio-optical algorithms for remote estimation of chlorophyll-a concentration (Chl) in case-1 waters exploit the upwelling radiation in the blue and green spectral regions. In turbid productive waters other constituents, that vary independently of Chl, absorb and scatter light in these spectral regions. As a consequence, the accurate estimation of Chl in turbid productive waters has so far not been feasible from satellite sensors. The main purpose of this study was to evaluate the extent to which near-infrared (NIR) to red reflectance ratios could be applied to the Sea Wide Field-of-View Sensor (SeaWiFS) and the Moderate Imaging Spectrometer (MODIS) to estimate Chl in productive turbid waters. To achieve this objective, remote-sensing reflectance spectra and relevant water constituents were collected in 251 stations over lakes and reservoirs with a wide variability in optical parameters (i.e. 4 ≤ Chl ≤ 240 mg m^− 3; 18 ≤ Secchi disk depth ≤ 308 cm). SeaWiFS and MODIS NIR and red reflectances were simulated by using the in-situ hyperspectral data. The proposed algorithms predicted Chl with a relative random uncertainty of approximately 28% (average bias between − 1% and − 4%). The effects of reflectance uncertainties on the predicted Chl were also analyzed. It was found that, for realistic ranges of R_rs uncertainties, Chl could be estimated with a precision better than 40% and an accuracy better than ± 35%. These findings imply that, provided that an atmospheric correction scheme specific for the red-NIR spectral region is available, the extensive database of SeaWiFS and MODIS images could be used to quantitatively monitor Chl in turbid productive waters.     

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.     

Vicarious calibration of GOCI for the SeaDAS ocean color retrieval

A key on-orbit calibration step for satellite remote sensing of ocean color is the vicarious calibration. This establishes the final gains for each spectral band on the sensor that minimize bias in the retrieved ocean color signal. The vicarious calibration is specific to the instrument and the atmospheric correction algorithm. The vicarious calibration gains for the Geostationary Ocean Color Imager (GOCI) are presented here, which were derived to optimize the performance of NASA’s standard atmospheric correction algorithm as implemented in the l2gen code and distributed through the SeaDAS open-source software package. Following NASA’s protocols, the near-infrared (NIR) bands were calibrated first, and the visible bands were then calibrated relative to this fixed NIR calibration. The gain for the 745-nm NIR band was derived using a fixed aerosol model, which was chosen based on the Angstrom Coefficients derived from MODIS on Aqua (MODISA). For the vicarious gains of the visible bands, two sources for the target water-leaving radiances were tested: matchups from MODISA and climatological data from SeaWiFS. A validation analysis using AERONET-OC data shows an improvement in sensor performance when compared with results using the current vicarious gains and results using no vicarious calibration. Good agreement was found in vicarious gains derived using both concurrent MODISA and climatological SeaWiFS as vicarious calibration data sources. These results support the use of a concurrent sensor for the vicarious calibration when in situ data are not available and demonstrate that using climatology from a well-calibrated sensor like SeaWiFS for the vicarious calibration is a valid alternative when it is not possible to use a concurrent sensor or in situ data. We recommend using the gains derived from concurrent GOCI matchups with MODISA for GOCI processing in SeaDAS/l2gen.     

An assessment of photosynthetic light use efficiency from space: Modeling the atmospheric and directional impacts on PRI reflectance

Estimation of photosynthetic light use efficiency (ε) from satellite observations is an important component of climate change research. The photochemical reflectance index, a narrow waveband index based on the reflectance at 531 and 570 nm, allows sampling of the photosynthetic activity of leaves; upscaling of these measurements to landscape and global scales, however, remains challenging. Only a few studies have used spaceborne observations of PRI so far, and research has largely focused on the MODIS sensor. Its daily global coverage and the capacity to detect a narrow reflectance band at 531 nm make it the best available choice for sensing ε from space. Previous results however, have identified a number of key issues with MODIS-based observations of PRI. First, the differences between the footprint of eddy covariance (EC) measurements and the MODIS footprint, which is determined by the sensor's observation geometry make a direct comparison between both data sources challenging and second, the PRI reflectance bands are affected by atmospheric scattering effects confounding the existing physiological signal. In this study we introduce a new approach for upscaling EC based ε measurements to MODIS. First, EC-measured ε values were “translated” into a tower-level optical PRI signal using AMSPEC, an automated multi-angular, tower-based spectroradiometer instrument. AMSPEC enabled us to adjust tower-measured PRI values to the individual viewing geometry of each MODIS overpass. Second, MODIS data were atmospherically corrected using a Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm, which uses a time series approach and an image-based rather than pixel-based processing for simultaneous retrievals of atmospheric aerosol and surface bidirectional reflectance (BRDF). Using this approach, we found a strong relationship between tower-based and spaceborne reflectance measurements (r² = 0.74, p < 0.01) throughout the vegetation period of 2006. Swath (non-gridded) observations yielded stronger correlations than gridded data (r² = 0.58, p < 0.01) both of which included forward and backscatter observations. Spaceborne PRI values were strongly related to canopy shadow fractions and varied with different levels of ε. We conclude that MAIAC-corrected MODIS observations were able to track the site-level physiological changes from space throughout the observation period.     

Spectral variations in the near-infrared ocean reflectance

The optical properties of natural waters beyond the visible range, in the near-infrared (NIR, 700-900 nm), have received little attention because they are often assumed to be mostly determined by the large absorption coefficient of pure water, and because of methodological difficulties. It is now growingly admitted that the NIR represents a potential optical source of unambiguous information about suspended sediments in turbid waters, thence the need for better understanding the NIR optical behaviour of such waters. It has recently been proposed (Ruddick et al., Limnology and Oceanography. 51, 1167-1179, 2006) that the variability in the shape of the surface ocean reflectance spectrum in the NIR is negligible in turbid waters. In the present study, we show, based on both in situ and remote sensing data, that the shape of the ocean reflectance spectrum in the NIR does vary in turbid to extremely turbid waters. Space-borne ocean reflectance data were collected using 3 different sensors (SeaWiFS, MODIS/Aqua and MERIS) over the Amazon, Mackenzie and Rio de la Plata turbid river plumes during extremely clear atmospheric conditions so that reliable removal of gas and aerosol effects on reflectance could be achieved. In situ NIR reflectance data were collected in different European estuaries where extremely turbid waters were found. In both data sets, a flattening of the NIR reflectance spectrum with increasing turbidity was observed. The ratio of reflectances at 765 nm and 865 nm, for instance, varied from ca. 2 down to 1 in our in situ data set, while a constant value of 1.61 had been proposed based on theory in a previous study. Radiative transfer calculations were performed using a range of realistic values for the seawater inherent optical properties, to determine the possible causes of variations in the shape of the NIR reflectance spectrum. Based on these simulations, we found that the most significant one was the gradual increase in the contribution of suspended sediments to the color of surface waters, which often leads to the flattening of the reflectance spectrum. Changes in the scattering and absorption properties of particles also contributed to variations in the shape of the NIR surface ocean reflectance spectrum. The impact of such variations on the interpretation of ocean color data is discussed.