Advancing regional satellite-based assessment of phytoplankton size structure in a subtropical Bight

This study presents an algorithm for globally retrieving chlorophyll a (Chl-a) concentrations of phytoplankton functional types (PFTs) from multi-sensor merged ocean color (OC) products or Sentinel-3A (S3) Ocean and Land Color Instrument (OLCI) data from the GlobColour archive in the frame of the Copernicus Marine Environmental Monitoring Service (CMEMS). The retrieved PFTs include diatoms, haptophytes, dinoflagellates, green algae and prokaryotic phytoplankton. A previously proposed method to retrieve various phytoplankton pigments, based on empirical orthogonal functions (EOF), is investigated and adapted to retrieve Chl-a concentrations of multiple PFTs using extensive global data sets of in situ pigment measurements and matchups with satellite OC products. The performance of the EOF-based approach is assessed and cross-validated statistically. The retrieved PFTs are compared with those derived from diagnostic pigment analysis (DPA) based on in situ pigment measurements. Results show that the approach predicts well Chl-a concentrations of most of the mentioned PFTs. The performance of the approach is, however, less accurate for prokaryotes, possibly due to their general low variability and small concentration range resulting in a weak signal which is extracted from the reflectance data and corresponding EOF modes. As a demonstration of the approach utilization, the EOF-based fitted models based on satellite reflectance products at nine bands are applied to the monthly GlobColour merged products. Climatological characteristics of the PFTs are also evaluated based on ten years of merged products (2002−2012) through inter-comparisons with other existing satellite derived products on phytoplankton composition including phytoplankton size class (PSC), SynSenPFT, OC-PFT and PHYSAT. Inter-comparisons indicate that most PFTs retrieved by our study agree well with previous corresponding PFT/PSC products, except that prokaryotes show higher Chl-a concentration in low latitudes. PFT dominance derived from our products is in general well consistent with the PHYSAT product. A preliminary experiment of the retrieval algorithm using eleven OLCI bands is applied to monthly OLCI products, showing comparable PFT distributions with those from the merged products, though the matchup data for OLCI are limited both in number and coverage. This study is to ultimately deliver satellite global PFT products for long-term continuous observation, which will be updated timely with upcoming OC data, for a comprehensive understanding of the variability of phytoplankton composition structure at a global or regional scale.

Chlorophyll-a (chl-a) is an important parameter of water quality and its concentration can be directly retrieved from satellite observations. The Ocean and Land Color Instrument (OLCI), a new-generation water-color sensor onboard Sentinel-3A and Sentinel-3B, is an excellent tool for marine environmental monitoring. In this study, we introduce a new machine learning model, Light Gradient Boosting Machine (LightGBM), for estimating time-series chl-a concentration in Fujian’s coastal waters using multitemporal OLCI data and in situ data. We applied the Case 2 Regional CoastColour (C2RCC) processor to obtain OLCI band reflectance and constructed four spectral indices based on OLCI feature bands as supplementary input features. We also used root-mean-square error (RMSE), mean absolute error (MAE), median absolute percentage error (MAPE), and R2 as performance indicators. The results indicate that the addition of spectral indices can easily improve the prediction accuracy of the model, and normalized fluorescence height index (NFHI) has the best performance, with an RMSE of 0.38 µg/L, MAE of 0.22 µg/L, MAPE of 28.33%, and R2 of 0.785. Moreover, we used the well-known band ratio and three-band methods for chl-a estimation validation, and another two OLCI chl-a products were adopted for comparison (OC4Me chl-a and Inverse Modelling Technique (IMT) Neural Net chl-a). The results confirmed that the LightGBM model outperforms the traditional methods and OLCI chl-a products. This study provides an effective remote sensing technique for coastal chl-a concentration estimation and promotes the advantage of OLCI data in ocean color remote sensing.

The high-resolution sea surface temperature (SST) is very useful for its coverage near to coast to study the fine scale processes and its impact on multi-scale atmospheric process. The work presents the validation of Multiscale Ultra High Resolution (MUR) SST products with three coastal buoy and three open ocean buoy observations in the Bay of Bengal (BoB). The comparisons show the correlation coefficient (R) more than 0.80 with coastal buoys and root mean square error (RMSE) less than 0. $1^{o}$C. In the open ocean, the value of R is higher and RMSE is lesser compared to the coastal observation. This is first of kind validation of the SST data with coastal buoys. The results provided the usefulness of the MUR SST in the BoB region. The front detection method is applied on December SST. A strong SST gradient of value 0. 1 $\circ$c. k$m^{-1}$ is observed in the head and westem BoB. We conclude that the MUR SST is good in quality in the coastal and open ocean to study the small scale feature in the BoB.

Abstract. As offshore wind farm development expands, accurate wind resource forecasting over the ocean is needed. One important yet relatively unexplored aspect of offshore wind resource assessment is the role of sea surface temperature (SST). Models are generally forced with reanalysis data sets, which employ daily SST products. Compared with observations, significant variations in SSTs that occur on finer timescales are often not captured. Consequently, shorter-lived events such as sea breezes and low-level jets (among others), which are influenced by SSTs, may not be correctly represented in model results. The use of hourly SST products may improve the forecasting of these events. In this study, we examine the sensitivity of model output from the Weather Research and Forecasting model (WRF) 4.2.1 to different SST products. We first evaluate three different data sets: the Multiscale Ultrahigh Resolution (MUR25) SST analysis, a daily, 0.25∘ × 0.25∘ resolution product; the Operational Sea Surface Temperature and Ice Analysis (OSTIA), a daily, 0.054∘ × 0.054∘ resolution product; and SSTs from the Geostationary Operational Environmental Satellite 16 (GOES-16), an hourly, 0.02∘ × 0.02∘ resolution product. GOES-16 is not processed at the same level as OSTIA and MUR25; therefore, the product requires gap-filling using an interpolation method to create a complete map with no missing data points. OSTIA and GOES-16 SSTs validate markedly better against buoy observations than MUR25, so these two products are selected for use with model simulations, while MUR25 is at this point removed from consideration. We run the model for June and July of 2020 and find that for this time period, in the Mid-Atlantic, although OSTIA SSTs overall validate better against in situ observations taken via a buoy array in the area, the two products result in comparable hub-height (140 m) wind characterization performance on monthly timescales. Additionally, during hours-long flagged events (< 30 h each) that show statistically significant wind speed deviations between the two simulations, both simulations once again demonstrate similar validation performance (differences in bias, earth mover's distance, correlation, and root mean square error on the order of 10−1 or less), with GOES-16 winds validating nominally better than OSTIA winds. With a more refined GOES-16 product, which has been not only gap-filled but also assimilated with in situ SST measurements in the region, it is likely that hub-height winds characterized by GOES-16-informed simulations would definitively validate better than those informed by OSTIA SSTs.

<strong class="journal-contentHeaderColor">Abstract.</strong> As offshore wind farm development expands, accurate wind resource forecasting over the ocean is needed. One important yet relatively unexplored aspect of offshore wind resource assessment is the role of sea surface temperature (SST). Models are generally forced with reanalysis data sets, which employ daily SST products. Compared with observations, significant variations in SSTs that occur on finer timescales are often not captured. Consequently, shorter-lived events such as sea breezes and low-level jets (among others), which are influenced by SSTs, may not be correctly represented in model results. The use of hourly SST products may improve the forecasting of these events. In this study, we examine the sensitivity of model output from the Weather Research and Forecasting model (WRF) 4.2.1 to different SST products. We first evaluate three different data sets: the Multiscale Ultrahigh Resolution (MUR25) SST analysis, a daily, 0.25<span class="inline-formula">^∘</span> <span class="inline-formula">×</span> 0.25<span class="inline-formula">^∘</span> resolution product; the Operational Sea Surface Temperature and Ice Analysis (OSTIA), a daily, 0.054<span class="inline-formula">^∘</span> <span class="inline-formula">×</span> 0.054<span class="inline-formula">^∘</span> resolution product; and SSTs from the Geostationary Operational Environmental Satellite 16 (GOES-16), an hourly, 0.02<span class="inline-formula">^∘</span> <span class="inline-formula">×</span> 0.02<span class="inline-formula">^∘</span> resolution product. GOES-16 is not processed at the same level as OSTIA and MUR25; therefore, the product requires gap-filling using an interpolation method to create a complete map with no missing data points. OSTIA and GOES-16 SSTs validate markedly better against buoy observations than MUR25, so these two products are selected for use with model simulations, while MUR25 is at this point removed from consideration. We run the model for June and July of 2020 and find that for this time period, in the Mid-Atlantic, although OSTIA SSTs overall validate better against in situ observations taken via a buoy array in the area, the two products result in comparable hub-height (140 <span class="inline-formula">m</span>) wind characterization performance on monthly timescales. Additionally, during hours-long flagged events (<span class="inline-formula">&lt;</span> 30 <span class="inline-formula">h</span> each) that show statistically significant wind speed deviations between the two simulations, both simulations once again demonstrate similar validation performance (differences in bias, earth mover's distance, correlation, and root mean square error on the order of 10<span class="inline-formula">^−1</span> or less), with GOES-16 winds validating nominally better than OSTIA winds. With a more refined GOES-16 product, which has been not only gap-filled but also assimilated with in situ SST measurements in the region, it is likely that hub-height winds characterized by GOES-16-informed simulations would definitively validate better than those informed by OSTIA SSTs.

Satellite remote sensing of turbidity in Lake Xingkai using eight years of OLCI observations.

The proxy for phytoplankton biomass, Chlorophyll a (Chl a), is an important variable to assess the health and state of the oceans which are under increasing anthropogenic pressures. Prior to the operational use of satellite ocean-colour Chl a to monitor the oceans, rigorous assessments of algorithm performance are necessary to select the most suitable products. Due to their inaccessibility, the oligotrophic open-ocean gyres are under-sampled and therefore under-represented in global in situ data sets. The Atlantic Meridional Transect (AMT) campaigns fill the sampling gap in Atlantic oligotrophic waters. In-water underway spectrophotometric data were collected on three AMT field campaigns in 2016, 2017 and 2018 to assess the performance of Sentinel-3A (S3-A) and Sentinel-3B (S3-B) Ocean and Land Colour Instrument (OLCI) products. Three Chl a algorithms for OLCI were compared: Processing baseline (pb) 2, which uses the ocean colour 4 band ratio algorithm (OC4Me); pb 3 (OL_L2M.003.00) which uses OC4Me and a colour index (CI); and POLYMER v4.8 which models atmosphere and water reflectance and retrieves Chl a as a part of its spectral matching inversion. The POLYMER Chl a for S-3A OLCI performed best. The S-3A OLCI pb 2 tended to under-estimate Chl a especially at low concentrations, while the updated OL_L2M.003.00 provided significant improvements at low concentrations. OLCI data were also compared to MODIS-Aqua (R2018 processing) and Suomi-NPP VIIRS standard products. MODIS-Aqua exhibited good performance similar to OLCI POLYMER whereas Suomi-NPP VIIRS exhibited a slight under-estimate at higher Chl a values. The reasons for the differences were that S-3A OLCI pb 2 Rrs were over-estimated at blue bands which caused the under-estimate in Chl a. There were also some artefacts in the Rrs spectral shape of VIIRS which caused Chl a to be under-estimated at values >0.1 mg m-3. In addition, using in situ Rrs to compute Chl a with OC4Me we found a bias of 25% for these waters, related to the implementation of the OC4ME algorithm for S-3A OLCI. By comparison, the updated OLCI processor OL_L2M.003.00 significantly improved the Chl a retrievals at lower concentrations corresponding to the AMT measurements. S-3A and S-3B OLCI Chl a products were also compared during the Sentinel-3 mission tandem phase (the period when S-3A and S-3B were flying 30 sec apart along the same orbit). Both S-3A and S-3B OLCI pb 2 under-estimated Chl a especially at low values and the trend was greater for S-3A compared to S-3B. The performance of OLCI was improved by using either OL_L2M.003.00 or POLYMER Chl a. Analysis of coincident satellite images for S-3A OLCI, MODIS-Aqua and VIIRS as composites and over large areas illustrated that OLCI POLYMER gave the highest Chl a concentrations and percentage (%) coverage over the north and south Atlantic gyres, and OLCI pb 2 produced the lowest Chl a and % coverage.

Abstract. The distribution and variation of phytoplankton size class (PSC) are key to understanding ocean biogeochemical processes and ecosystems. Remote sensing of the PSC in the East China Sea (ECS) remains a challenge, although many algorithms have been developed to estimate PSC. Here based on a local dataset from the ECS, a regional model was tuned to estimate the PSC from the spectral features of normalized phytoplankton absorption (aph) using a principal component analysis approach. Before applying the refined PSC model to MODIS (Moderate Resolution Imaging Spectroradiometer) data, reconstructing satellite remote sensing reflectance (Rrs) at 412 and 443 nm was critical through modeling them from Rrs between 469 and 555 nm using multiple regression analysis. Satellite-derived PSC results compared well with those derived from pigment composition, which demonstrated the potential of satellite ocean color data to estimate PSC distributions in the ECS from space. Application of the refined PSC model to the reconstructed MODIS data from 2003 to 2016 yielded the seasonal distributions of the PSC in the ECS, suggesting that the PSC distributions were heterogeneous in both temporal and spatial scales. Micro-phytoplankton were dominant in coastal waters throughout the year, especially in the Changjiang estuary. For the middle shelf region, the seasonal shifts from the dominance of micro- and nano-phytoplankton in the winter and spring to the dominance of nano- and pico-phytoplankton in the summer and autumn were observed. Pico-phytoplankton were especially dominant in the Kuroshio region in the spring, summer, and autumn. The seasonal variations of the PSC in the ECS were probably affected by a combination of the water column stability, upwelling, sea surface temperature, and the Kuroshio. Additionally, human activity and riverine discharge might also influence the PSC distribution in the ECS, especially in the coastal region.

High quality independent ground measurements that are traceable to metrology standards, with a full uncertainty budget, are required for validation over the lifetime of ocean-colour satellite missions. In this paper, we used radiometric Fiducial Reference Measurements (FRM) collected during four Atlantic Meridional Transect (AMT) field campaigns from 2016 to 2019 to assess the performance of radiometric products from the Ocean and Land Colour Instrument (OLCI) aboard Sentinel-3A (S-3A) and 3B (S‐3B), the Moderate Resolution Imaging Spectroradiometer instrument aboard Aqua (MODIS-Aqua), and the Visible Infrared Imaging Radiometer Suite instrument aboard Suomi NPP and NOAA-20 (Suomi-VIIRS and NOAA-20 VIIRS). The AMT provides one of the few sampling platforms that make high-quality in situ radiometric measurements in oligotrophic, low chlorophyll-a oceanic waters for ocean colour satellite validation. In situ data were acquired and processed following established FRM protocols, calibrated to metrology standards, referenced to inter-comparison exercises and with a full uncertainty budget. From these we selected an uncertainty threshold, which we used as part of a matchup procedure that takes into account the temporal and spatial variability of both the in situ and satellite data. Three atmospheric correction models were compared for S-3A and S‐3B OLCI radiometric products; the standard OLCI IPF-OL-2, POLYMER and NASA SeaDAS l2gen. Based on the round-robin comparison, POLYMER provided the best performance in the retrieval of water-leaving radiances. The analysis showed that Suomi-VIIRS and MODIS-Aqua performed better than NOAA-20 VIIRS, and comparably with S‐3B OLCI standard products. The S-3A OLCI standard product outperformed the NASA products. The S-3A OLCI and S‐3B OLCI instruments were also compared during their tandem phase, which showed that S‐3B OLCI radiances were systematically higher than S-3A OLCI across the spectrum.

Towards a novel approach for Sentinel-3 synergistic OLCI/SLSTR cloud and cloud shadow detection based on stereo cloud-top height estimation

&amp;lt;p&amp;gt;The first Copernicus Sentinel-3 satellite, Sentinel-3A, was launched in early 2016, and its twin Sentinel-3B in April 2018. The Sentinel-3 constellation is now fully operational with Sentinel-3B satellite flying in the same orbit plan with a phase difference of 140&amp;amp;#176;. This constellation provides a unique consistent, long-term collection of marine and land data for operational analysis, forecasting and environmental and climate monitoring. The marine centre is part of the Sentinel-3 Payload Data Ground Segment, located at EUMETSAT. This centre together with the existing EUMETSAT facilities provides a routine centralised service for operational meteorology, oceanography, and other Sentinel-3 marine users as part of the European Commission's Copernicus programme. The EUMETSAT marine centre delivers operational Sea Surface Temperature, Ocean Colour and Sea Surface Topography data products based on the measurements from the Sea and Land Surface Temperature Radiometer (SLSTR), Ocean and Land Colour Instrument (OLCI) and Synthetic Aperture Radar Altimeter (SRAL), all aboard Sentinel-3 satellites. All products have been developed together with ESA and industry partners and EUMETSAT is responsible for the production, distribution, performance and future evolution of Level-2 marine products. We will give an overview of the scientific characteristics and algorithms of all marine Level-2 products, as well as instrument calibration and product validation results based on on-going Sentinel-3 Cal/Val activities. Information will be also provided about the current status of the product dissemination and the future evolutions that are envisaged. Also, we will provide information how to access Sentinel-3 data from EUMETSAT and where to look for further information.&amp;lt;/p&amp;gt;

Sea Surface Temperature (SST) products have been available from a number of operational and experimental satellites for more than thirty years. These SST products have been evaluated in the global as well as regional areas. In China, the meteorological satellites, Fengyun (FY) series and oceanic satellites, Haiyang (HY) series, have been designed with infrared or microwave channels for SST observations. The FY-2 series are geostationary-orbiting satellites. The Visible and Infrared Spin Scan Radiometer (VISSR-2) onboard FY-2C, 2D, 2E has infrared split window channels for the observation of SST. The brightness temperature data of the FY-2D are analyzed and compared with simultaneous brightness temperature data from the Imager onboard the Multi-functional Transport Satellite (MTSAT-1R). The results show poor calibration of FY-2D VISSR-2 infrared channels, which are not capable of retrieving valid SST products. The Chinese Ocean Color and Temperature Scanner (COCTS) onboard HY-1A, 1B also has infrared split window channels. The SST products are operational delivered. The intercomparisons of COCTS SST products with AVHRR and MODIS SST products are carried out. The results show negative bias exists for COCTS SST products. The SST products derived from the Visible and Infrared Radiometer (VIRR) onboard FY-3A are under investigation.

Surface oceanic fronts are regions characterized by high biological activity. Here, Sea Surface Temperature (SST) fronts are analyzed for the period 2003–2019 using the Multi-scale Ultra-high Resolution (MUR) SST product in northern Patagonia, a coastal region with high environmental variability through river discharges and coastal upwelling events. SST gradient magnitudes were maximum off Chiloé Island in summer and fall, coherent with the highest frontal probability in the coastal oceanic area, which would correspond to the formation of a coastal upwelling front in the meridional direction. Increased gradient magnitudes in the Inner Sea of Chiloé (ISC) were found primarily in spring and summer. The frontal probability analysis revealed the highest occurrences were confined to the northern area (north of Desertores Islands) and around the southern border of Boca del Guafo. An Empirical Orthogonal Function analysis was performed to clarify the dominant modes of variability in SST gradient magnitudes. The meridional coastal fronts explained the dominant mode (78% of the variance) off Chiloé Island, which dominates in summer, whereas the SST fronts inside the ISC (second mode; 15.8%) were found to dominate in spring and early summer (October–January). Future efforts are suggested focusing on high frontal probability areas to study the vertical structure and variability of the coastal fronts in the ISC and its adjacent coastal ocean.

A joint effort between the Copernicus Climate Change Service (C3S) and the Group for High Resolution Sea Surface Temperature (GHRSST) has been dedicated to an intercomparison study of eight global gap-free sea surface temperature (SST) products to assess their accurate representation of the SST relevant to climate analysis. In general, all SST products show consistent spatial patterns and temporal variability during the overlapping time period (2003–18). The main differences between each product are located in the western boundary current and Antarctic Circumpolar Current regions. Linear trends display consistent SST spatial patterns among all products and exhibit a strong warming trend from 2012 to 2018 with the Pacific Ocean basin as the main contributor. The SST discrepancy between all SST products is very small compared to the significant warming trend. Spatial power spectral density shows that the interpolation into 1° spatial resolution has negligible impacts on our results. The global mean SST time series reveals larger differences among all SST products during the early period of the satellite era (1982–2002) when there were fewer observations, indicating that the observation frequency is the main constraint of the SST climatology. The maturity matrix scores, which present the maturity of each product in terms of documentation, storage, and dissemination but not the scientific quality, demonstrate that ESA-CCI and OSTIA SST are well documented for users’ convenience. Improvements could be made for MGDSST and BoM SST. Finally, we have recommended that these SST products can be used for fundamental climate applications and climate studies (e.g., El Niño).

An underwater volcanic eruption off the Vava’u island group in Tonga on 7 August 2019 resulted in the creation of floating pumice on the ocean’s surface extending over an area of 150 km2. The pumice’s far-reaching effects from its origin in the Tonga region to Fiji and the methods of automatic detection using satellite imagery are described, making it possible to track the westward drift of the pumice raft over 43 days. Level 2 Moderate Resolution Imaging Spectroradiometer (MODIS), Visible Infrared Imaging Radiometer Suite (VIIRS), Sentinel-3 Ocean and Land Color Instrument (OLCI), and Sentinel-3 Sea and Land Surface Temperature Radiometer (SLSTR) imagery of sea surface temperature, chlorophyll-a concentration, quasi-surface (i.e., Rayleigh-corrected) reflectance, and remote sensing reflectance were used to distinguish consolidated and fragmented rafts as well as discolored and mesotrophic waters. The rafts were detected by a 1 to 3.5 °C enhancement in the MODIS-derived “sea surface temperature” due to the emissivity difference of the raft material. Large plumes of discolored waters, characterized by higher satellite reflectance/backscattering of particles in the blue than surrounding waters (and corresponding to either submersed pumice or associated white minerals), were associated with the rafts. The discolored waters had relatively lower chlorophyll-a concentration, but this was artificial, resulting from the higher blue/red reflectance ratio caused by the reflective pumice particles. Mesotrophic waters were scarce in the region of the pumice rafts, presumably due to the absence of phytoplanktonic response to a silicium-rich pumice environment in these tropical oligotrophic environments. As beach accumulations around Pacific islands surrounded by coral shoals are a recurrent phenomenon that finds its origin far east in the ocean along the Tongan trench, monitoring the events from space, as demonstrated for the 7 August 2019 eruption, might help mitigate their potential economic impacts.