Research Article - (2018) Volume 21, Issue 2
Keywords: Sea surface irradiance; SeaWiFS; MODIS; MERIS; Generalized linear model
The monitoring of photosynthetically available radiation (PAR) in the oceanic environment is relevant for meteo-oceanographic processes such as heat fluxes within the surface layer [1]. Diurnal sea surface temperature (SST) variability [2] and the mixed-layer deepening/shoaling [1]. Biologically, PAR regulates phytoplankton distribution in the water column considering the availability of light energy at different depths [3]. Estimation of daily PAR (Einstein m-2 d-1) from ocean colour remote sensing (OCRS) is also important for monitoring the oceanic primary productivity (PP) and the subsequent assimilation of carbon by phytoplankton in the photosynthesis process [4]. Accurate estimation of daily PAR from satellite observations is therefore a prerequisite to provide a global coverage of physical and biogeochemical parameters [4,5].
The assessment of PAR algorithms from different sensors and atmospheric conditions requires comprehensive efforts of validation by in situ measurements [6-12]. Satellite estimates of daily PAR has also been evaluated in PP exercises [13-17]. Although PAR applications are considered relevant in meteo-oceanographic studies and primary productivity exercises, orbital data assessments against observational data measurements in the Southwestern South Atlantic are still scarce in the literature [15] especially in Brazilian waters, which makes it difficult to evaluate the spatiotemporal variability of incoming PAR and its dependence on forcing factors.
The present study aims to evaluate OCRS PAR estimates in the southeastern Brazilian continental margin, located at the Southwestern South Atlantic continental margin, with sensors which have a sufficient set of spectral bands, such as Sea-viewing Wide Field-of-view Sensor (SeaWiFS), MODerate Imaging Spectroradiometer onboard on Aqua satellite (MODIS) and MEdium Resolution Imaging Spectrometer (MERIS). All those sensors were operational during the periods of in situ data acquisition. An empirical regional model was adjusted with in situ data collected in summer and winter of 2001 and 2002 in the study region.
In situ PAR data
In situ sampling was conducted during 4 mesoscale cruises in austral summer and winter of 2001-2002 in a region located at the northern portion of the Brazilian Southeastern continental margin delimited between Cape of São Tomé (ST), in Rio de Janeiro state (22°S), and São Sebastião Island (SSI) in São Paulo state (24°S), Southwestern South Atlantic (Figure 1a and 1b). PAR data were acquired with a quantum scalar surface reference sensor QSR-240 (Biospherical Instruments Inc.). A spectral irradiance model as described by [18] was also used to determine the available energy in the sea surface, validated with the in situ measurements. The total irradiance (Watt m-2) was calculated as a function of time (t, hours), geographic position (latitude) according to the in situ station coordinates, cloud cover observations during the cruises, day length (hours) and day of the year (Julian calendar). Some adjustments were applied for seasonal variations as suggested by [19].
Satellite PAR data
Level 1 SeaWiFS daily images were acquired during the austral summer and winter cruises of 2001 and 2002. As MERIS and MODIS onboard satellite Aqua data are only available starting from May and July of 2002, respectively, concomitant daily PAR was acquired only for the winter cruise of 2002. The images for the three sensors were obtained from the OceanColor Web page (http://oceancolor.gsfc.nasa.gov/cms/) supported by the Ocean Biology Processing Group (OBPG) at NASA's Goddard Space Flight Center. The images were processed to Level 2 daily PAR product according this algorithm was first applied to SeaWifs [20,21] and currently is also used on MODIS and MERIS (among other sensors).
OCRS estimation of PAR (Einstein m-2 d-1) is derived from the solar irradiance (Es, mW cm-2 μm-1) integrated in the visible range of the electromagnetic radiation (400-700 nm) over the day length defined by latitude and date of acquisition [21]. The implementation of this algorithm (Equation (1)) depends on the availability of the irradiance at the top of the atmosphere limited by saturation clouds:
(1)
Wher Ed(λ) (mW cm-2 μm-1) is the downward irradiance after the interaction with the atmosphere, Sα refers to the spherical albedo, A is the albedo of clouds and aerosols on cloud-surface path and can be reduced to Sα when in ideal weather conditions. ES is the solar irradiance that should reach the sea surface if A did not exist [22]. At the end of the process, the PAR is obtained in units of mW cm-2 μm-1 and converted to Einstein m-2 d-1 by a factor of 1.193 with a small percentage of error regardless of weather conditions [23]. The atmospheric correction algorithm is described in Ref. [24] based on [25] according to Equation (2):
Lt=[Lr+La+tdv+Lwc+tdv+Lw]tgvtgsfp (2)
Where Lt is the total upwelling radiance observed by the satellite sensor after the interaction with the ocean and atmosphere Lr is the contribution from the Rayleigh molecular scattering La is the aerosol contribution Lwc is the contribution by whitecaps and foam above the sea surface and Lw is the water-leaving radiance. The radiance L is spectral resolved and has units of Wm-2 nm-1 sr-1 . The term is the transmittance from diffuse radiation along the ocean to sensor trajectory tgv and tgs are the transmittance by atmospheric gases in the viewing sensor direction and in the Sun’s direction. Finally fp is a polarization-adjustment factor. Both daily PAR and atmospheric correction algorithms were applied to the L1 satellite images following the NASA OBPG descriptions using the SeaWiFS Data Analysis System (SeaDAS) version 7.3.
Level 2 images were reprojected to the geographic coordinate system Datum World Geodetic System 1984 (WGS84), preserving the nominal spatial resolution of 1 km at nadir for the whole scene using the nearest neighbour algorithm. When two or more images of the same sensor on the same day were available, an average composite was calculated for overlaid pixels. The study area was defined in terms of latitudes 20°S – 26°S and longitudes 40°W – 46°W (Figure 1b).
Statistical analysis
Algorithm performance was assessed based on the root mean square error (RMSE), the average error (bias), the percentage error %error and the coefficient o f determination (r2). The criteria to determine the best relative performance is based on the RMSE, bias, %error and r2, in this sequence, as suggested by [26]. A box plot chart where values below {Q1‐1.5×(Q3‐Q1)} and above {Q3+1.5×(Q3‐Q1)} where Q1 and Q3 are the first and third quartiles, respectively, are identified a priori as outliers. Outliers were checked by visual interpretation of the corresponding satellite image confirming or not a possible pixel flag, and then removed previously to statistical analysis.
Algorithm performance assessment
Overall, PAR derived with SeaWiFS, MODIS and MERIS data slightly overestimates the in situ measurements (Figure 2a-2c). The best performance was obtained for SeaWiFS PAR with RMSE=1.77 Einstein m-2 d-1 , bias=1.63 Einstein m-2 d-1 that corresponds to a %error=3.95% higher than the measured values and r2 ≈ 1.0 (Figure 2(a)). Sequentially, the PAR estimates with MODIS data (Figure 2(b)) were statistically better than with MERIS data (Figure 2(c)). Scatterplots of each satellite PAR against in situ observations reinforced the statistics (Figure 2a-2c).
Main sources of uncertainties are related to cloud cover, time differences between satellite overpass and in situ sampling, and absorption throughout the optical path Frouin et al. [9]. compared SeaWiFS and MODIS PAR values with in situ observations measured at a stationary site off Chesapeake Bay in the North Atlantic. Results were presented as daily, weekly, and monthly uncertainties averaged between 2005 and 2010, also showing an overestimation for SeaWiFS(MODIS) daily PAR RMSE=6.49(6.77) Einstein m-2 d-1, bias = 2.83(1.85) Einstein m-2 d-1 and r2 = 0.87(0.86). Their result corroborates with the performance observed in the present study before the GLM adjustment.
SeaWiFS PAR product has a disadvantage which is the lack of diurnal variability in clouds, due to the use of a single sensor in a nearnoon orbit [27]. For completely clear sky situations, the PAR estimates derived from satellite data are in much better agreement with the measurements. Another input error of about 1 Einstein m-2 d-1 (2 to 3%) is attributed to the accuracy of the irradiance model used in conjunction with the field measurement [9]. There was a small seasonal variation in the ratio of satellite-derived and measured PAR values. [21] analysed SeaWiFS PAR estimates in relation to measurements from two moored buoys, one at the relatively high latitude of British Columbia (Halibut Bank, 49°N) and the other in the Equatorial Pacific (0°N). The authors observed a RMSE equal to 6.2 Einstein m-2 d-1 for both sites and a lower bias in Halibut (0.93) than in the Pacific (2.9 Einstein m-2 d-1).
A wide range of %error (5% to 73%)was observed by Vazyula et al. [11] when comparing MODIS PAR with in situ measurements taken during a scientific cruise from the Baltic to the White Sea from the end of July to the beginning of August 2014. Lalibert et al. [10] evaluated MODIS PAR Level-3 processing at high northern latitudes and obtained a %error between 17% and 20%. In the present study, the %error of satellite daily PAR estimates was ≤ 4.54% which is considered satisfactory for modelling oceanic PP with an accuracy higher than 90% [28]. However, the performance statistics show a systematic tendency of overestimation that may impact some PP models estimates up to 29% [13,14].
Adjustment of a generalized linear model (GLM)
Considering the observed systematic difference, a Generalized Linear Model was adjusted to all three sensors individually, according to Equations (3)-(5):
Sea WiFSPAR*= 0.99 [SeaWiFS PAR]-1.32 (3)
MODIS PAR* = 0.97[MODIS] - 0.37 (4)
MERIS PAR* = 0.32[MERIS PAR] + 24.21 (5)
Where PAR* is the adjusted value for each specific sensor; the first term on right is the slope and the second term is the offset Frouin et al. [9]. Inferred that the PAR values may be corrected (reduced) by a factor of 1.02 and 1.03 Einstein m-2 d-1 for SeaWiFS and MODIS PAR since PAR was computed as the product of the clear sky value. In the present study we found a reduction of 1.32 and 0.37 Einstein m-2 d-1, for SeaWiFS and MODIS PAR, respectively. With the adjusted GLM, SeaWiFS, MODIS and MERIS PAR*RMSE values decreased to 0.68, 0.34, and 1.61 Einstein m-2 d-1 respectively, while bias decreased to close zero values for all sensors (Figure 2d-2f). After GLM adjustment, the best performance was obtained with MODIS PAR* followed by SeaWiFS PAR*and MERIS PAR*. The relatively better performance obtained with MODIS PAR*in relation to SeaWiFS PAR*can be associated to less cloudiness during the winter in response to precipitation variability characterized by a rainy summer in the region [29]. Also, typical lower sea surface temperature values during winter may reduce the evaporation and consequently the formation of clouds during the day resulting in better estimates from OCRS.
Dogliotti et al. [15] observed a bias equal to 10 Einstein m-2 d-1( %error~ 48%) with MODIS PAR in the Argentine shelf (39°S to 55°S) and adjacent region during austral spring, late summer and late winter seasons. According to those authors the PAR estimation had a good r2 in relation with in situ measurements explaining ~70% of variance (N=36) suggesting that a GLM adjustment should be applied.
Accuracy of heliosynchronous satellite daily PAR products is limited by the absence of information about diurnal variability of clouds [9]. Hourly observations can be obtained from geostationary satellites [30] and integrated over time to provide daily values Qi et al. [5]. compared PAR estimates from MODIS and Geostationary Operational Environmental Satellite system (GOES) in the Gulf of Mexico region for the entire year of 2013. According to the authors, MODIS daily PAR tends to be lower than GOES PAR suggesting that when MODIS PAR is used in estimating PP, PP may be underestimated.
Ramon et al. [17] developed a MERIS PAR algorithm based on the NASA OBPG operational algorithm [21] accounting for a diurnal variability of clouds. The daily MERIS PAR estimates were compared with in situ measurements acquired at four sites located at midlatitudes: BOUSSOLE buoy in the Mediterranean Sea (2009-2015), CCE1 and CCE2 moorings off the California coast (2008-2015 and 2011-2015, respectively) and COVE off the East cost of United States (2003-2014). The RMSE error was 8.10 (%error= 23.74%) Einstein m d without statistical correction of diurnal cloud variability and 8.5 (%error= 24.90%) Einstein m-2 d-1 with cloud correction. The bias was reduced from 3.27 Einstein m-2 d-1 (%error= 9.57%) without cloud correction to 2.65 Einstein m-2 d-1 ( %error= 7.77%) with cloud correction. In both cases, r2 values were > 0.80 suggesting that an adjustment of a GLM could practically approach to zero value the systematic overestimation.
Comparisons of SeaWiFS, MODIS and MERIS PAR estimates with in situ measurements show agreement, with more than 90% of accuracy. However, all three satellite estimates were biased, by about 3.95% (SeaWIFS), 4.13% (MODIS) and 4.54% (MERIS). The overestimation may be due to combined effects of cloud cover, atmospheric corrections, satellite’s overpass time and sensor characteristics. A generalized linear model was adjusted regionally correcting the observed systematic differences. The regionally adjusted GLM effectively improved the estimation of PAR derived with SeaWiFS, MODIS and MERIS data. MERIS PAR still requires an additional effort for better tuning which could be worthwhile given that Ocean Land Colour Instrument (OLCI) has been designed with similar specifications with addition of six extra spectral bands. In this way PAR estimates derived from other ocean colour sensors currently in operation such as the Visible Infrared Imaging Radiometer Suite (VIIRS) should also be validated in the study region. The combination of geostationary and polar orbit satellites may allow retrieving information about daily atmospheric/clouds variability decreasing the uncertainties still present in current PAR algorithms and models. Important to note that the comparative analyses has been made using a relatively limited data set, which is insufficient to derive definite conclusions about accuracy. Additional imagery and in situ measurements needs to be analysed, considering different atmospheric conditions and seasons. At present, an effort is being made to bring together a longer time series of data acquired in the study region as part of the ANTARES Latin-American network (www.antares.ws). The assessment of PAR products is part of a larger effort to define the best input algorithms and products to be applied in regional primary productivity models.
These results showed that the GLM was a good choice to correct the systematic deviation present on estimates with the three sensors and that a simple linear model was able to improve the PAR estimate regionally.
DEPROAS project was funded by CNPq/PRONEX-FINEP Brazilian Government. The present study was partially funded by the Inter- American Institute for Global Changes Research (IAI/US-NSF GEO-1128040 grant) and National Council for Scientific and Technological Development (CNPq, 132971/2014-3). The authors would like to thank the Ocean Biology Processing Group (OBPG) for SeaDAS and OCRS data distribution.