Space-borne study of seasonal, multi-year, and decadal phytoplankton dynamics in the Bay of Biscay

Seasonal and inter-annual variations in phytoplankton community abundance in the Bay of Biscay are studied. Preliminarily processed by the National Aeronautics and Space Administration (NASA) to yield normalized water-leaving radiance and the top-of-the-atmosphere solar radiance, Sea-viewing Wide Field-of-View Sensor (SeaWiFS), Moderate Resolution Imaging Spectroradiometer (MODIS), and Coastal Zone Color Scanner (CZCS) data are further supplied to our dedicated retrieval algorithms to infer the sought for parameters. By applying the National Oceanic and Atmospheric Administration's (NOAA's) Advanced Very High Resolution Radiometer (AVHRR) data, the surface reflection coefficient in the only band in the visible spectrum is derived and employed for analysis. Decadal bridged time series of variations of diatom-dominated phytoplankton and green dinoflagellate Lepidodinium chlorophorum within the shelf zone and the coccolithophore Emiliania huxleyi in the pelagic area of the Bay are documented and analysed in terms of impacts of some biogeochemical and geophysical forcing factors.

The global distribution of phytoplankton size spectrum and size classes from their light-absorption spectra derived from satellite data

The absorption spectra of phytoplankton in the visible domain hold implicit information on the phytoplankton community structure. Here we use this information to retrieve quantitative information on phytoplankton size structure by developing a novel method to compute the exponent of an assumed power-law for their particle-size spectrum. This quantity, in combination with total chlorophyll-a concentration, can be used to estimate the fractional concentration of chlorophyll in any arbitrarily-defined size class of phytoplankton. We further define and derive expressions for two distinct measures of cell size of mixed. populations, namely, the average spherical diameter of a bio-optically equivalent homogeneous population of cells of equal size, and the average equivalent spherical diameter of a population of cells that follow a power-law particle-size distribution. The method relies on measurements of two quantities of a phytoplankton sample: the concentration of chlorophyll-a, which is an operational index of phytoplankton biomass, and the total absorption coefficient of phytoplankton in the red peak of visible spectrum at 676 nm. A sensitivity analysis confirms that the relative errors in the estimates of the exponent of particle size spectra are reasonably low. The exponents of phytoplankton size spectra, estimated for a large set of in situ data from a variety of oceanic environments (similar to 2400 samples), are within a reasonable range; and the estimated fractions of chlorophyll in pico-, nano- and micro-phytoplankton are generally consistent with those obtained by an independent, indirect method based on diagnostic pigments determined using high-performance liquid chromatography. The estimates of cell size for in situ samples dominated by different phytoplankton types (diatoms, prymnesiophytes, Prochlorococcus, other cyanobacteria and green algae) yield nominal sizes consistent with the taxonomic classification. To estimate the same quantities from satellite-derived ocean-colour data, we combine our method with algorithms for obtaining inherent optical properties from remote sensing. The spatial distribution of the size-spectrum exponent and the chlorophyll fractions of pico-, nano- and micro-phytoplankton estimated from satellite remote sensing are in agreement with the current understanding of the biogeography of phytoplankton functional types in the global oceans. This study contributes to our understanding of the distribution and time evolution of phytoplankton size structure in the global oceans.

Surface Ocean CO2 Atlas (SOCAT) gridded data products

As a response to public demand for a well-documented, quality controlled, publically available, global surface ocean carbon dioxide (CO2) data set, the international marine carbon science community developed the Surface Ocean CO2 Atlas (SOCAT). The first SOCAT product is a collection of 6.3 million quality controlled surface CO2 data from the global oceans and coastal seas, spanning four decades (1968–2007). The SOCAT gridded data presented here is the second data product to come from the SOCAT project. Recognizing that some groups may have trouble working with millions of measurements, the SOCAT gridded product was generated to provide a robust, regularly spaced CO2 fugacity (fCO2) product with minimal spatial and temporal interpolation, which should be easier to work with for many applications. Gridded SOCAT is rich with information that has not been fully explored yet (e.g., regional differences in the seasonal cycles), but also contains biases and limitations that the user needs to recognize and address (e.g., local influences on values in some coastal regions).

Coccolithophore surface distributions in the North Atlantic and their modulation of the air-sea flux of CO2 from 10 years of satellite Earth observation data

Coccolithophores are the primary oceanic phytoplankton responsible for the production of calcium carbonate (CaCO3). These climatically important plankton play a key role in the oceanic carbon cycle as a major contributor of carbon to the open ocean carbonate pump (similar to 50 %) and their calcification can affect the atmosphere-to-ocean (air-sea) uptake of carbon dioxide (CO2) through increasing the seawater partial pressure of CO2 (pCO(2)). Here we document variations in the areal extent of surface blooms of the globally important coccolithophore, Emiliania huxleyi, in the North Atlantic over a 10-year period (1998-2007), using Earth observation data from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS). We calculate the annual mean sea surface areal coverage of E. huxleyi in the North Atlantic to be 474 000 +/- 104 000 km(2), which results in a net CaCO3 carbon (CaCO3-C) production of 0.14-1.71 Tg CaCO3-C per year. However, this surface coverage (and, thus, net production) can fluctuate inter-annually by -54/+81% about the mean value and is strongly correlated with the El Nino/Southern Oscillation (ENSO) climate oscillation index (r = 0.75, p < 0.02). Our analysis evaluates the spatial extent over which the E. huxleyi blooms in the North Atlantic can increase the pCO(2) and, thus, decrease the localised air-sea flux of atmospheric CO2. In regions where the blooms are prevalent, the average reduction in the monthly air-sea CO2 flux can reach 55%. The maximum reduction of the monthly air-sea CO2 flux in the time series is 155 %. This work suggests that the high variability, frequency and distribution of these calcifying plankton and their impact on pCO(2) should be considered if we are to fully understand the variability of the North Atlantic air-to-sea flux of CO2. We estimate that these blooms can reduce the annual N. Atlantic net sink atmospheric CO2 by between 3-28 %.

Coccolithophore species as indicators of surface oceanographic conditions in the vicinity of Azores islands

During summer 2008 and spring 2009, surface oceanographic surveys were carried out around three islands of the Azores archipelago (Terceira, Sao Miguel and Santa Maria) to assess the phytoplankton distribution and associated physico-chemical processes. The Azores archipelago is a major feature in the biogeochemical North Atlantic Subtropical Gyre (NAST) province although its influence on the productivity of the surrounding ocean is poorly known. Surface phytoplankton was studied by microscopy and HPLC (High Precision Liquid Chromatography). The mean values for biomass proxy Chlorophyll a (Chla) ranged from 0.04 to 0.55 mu g L-1 (Chla maximum = 0.86 mu g L-1) and coccolithophores were the most abundant group, followed by small flagellates, Cyanobacteria, diatoms and dinoflagellates being the least abundant group. The distribution of phytoplankton and coccolithophore species in particular presented seasonal differences and was consistent with the nearshore influence of warm subtropical waters from the south Azores current and colder subpolar waters from the north. The satellite-derived circulation patterns showed southward cold water intrusions off Terceira and northward warm water intrusions off Santa Maria. The warmer waters signal was confirmed by the subtropical coccolithophore assemblage, being Discosphaera tubifera a constant presence under these conditions. The regions of enhanced biomass, either resulting from northern cooler waters or from island induced processes, were characterized by the presence of Emiliania huxleyi. Diatoms and dinoflagellates indicated coastal and regional processes of nutrient enrichment and areas of physical stability, respectively.

Variation in ocean colour may help predict cod and haddock recruitment

Characteristics of the spring and fall phytoplankton blooms in spawning areas on the Scotian Shelf, Canada, were estimated from remote sensing data. These blooms, along with anomalies in the North Atlantic Oscillation, were used to explain variation in the recruitment of 4 populations of cod and haddock. We tested the effects of the timing of the bloom using the chlorophyll a (chl a) signal, the maximum amount of chl a, the timing of the diatom bloom, and the maximum relative dominance of diatoms on the recruitment (to Age 1) of cod and haddock on the Scotian Shelf. Models were run separately for the effects of the spring and fall blooms. Only 3 of 10 models tested (0-lag) explained significant (80 to 92%) variation in recruitment. However, the performance of these models was not consistent across populations or species, suggesting that generalities about how spring and fall phytoplankton blooms affect recruitment cannot yet be made. The differences among models suggest that fish larvae are probably adapted locally to food production and thus indirectly to the characteristics of the phytoplankton bloom, which in turn are influenced by regional (meso-scale) oceanographic conditions.

Generalized ocean color inversion model for retrieving marine inherent optical properties

Ocean color measured from satellites provides daily, global estimates of marine inherent optical properties (IOPs). Semi-analytical algorithms (SAAs) provide one mechanism for inverting the color of the water observed by the satellite into IOPs. While numerous SAAs exist, most are similarly constructed and few are appropriately parameterized for all water masses for all seasons. To initiate community-wide discussion of these limitations, NASA organized two workshops that deconstructed SAAs to identify similarities and uniqueness and to progress toward consensus on a unified SAA. This effort resulted in the development of the generalized IOP (GIOP) model software that allows for the construction of different SAAs at runtime by selection from an assortment of model parameterizations. As such, GIOP permits isolation and evaluation of specific modeling assumptions, construction of SAAs, development of regionally tuned SAAs, and execution of ensemble inversion modeling. Working groups associated with the workshops proposed a preliminary default configuration for GIOP (GIOP-DC), with alternative model parameterizations and features defined for subsequent evaluation. In this paper, we: (1) describe the theoretical basis of GIOP; (2) present GIOP-DC and verify its comparable performance to other popular SAAs using both in situ and synthetic data sets; and, (3) quantify the sensitivities of their output to their parameterization. We use the latter to develop a hierarchical sensitivity of SAAs to various model parameterizations, to identify components of SAAs that merit focus in future research, and to provide material for discussion on algorithm uncertainties and future emsemble applications.

Real-time sea surface air temperature retrieval over Southern Ocean using AMSR-E data

The AMSR-E satellite data and in-situ data were applied to retrieve sea surface air temperature (Ta) over the Southern Ocean. The in-situ data were obtained from the 24~(th) -26~(th) Chinese Antarctic Expeditions during 2008-2010. First, Ta was used to analyze the relativity with the bright temperature (Tb) from the twelve channels of AMSR-E, and no high relativity was found between Ta and Tb from any of the channels. The highest relativity was 0.38 (with 23.8 GHz). The dataset for the modeling was obtained by using in-situ data to match up with Tb, and two methods were applied to build the retrieval model. In multi-parameters regression method, the Tbs from 12 channels were used to the model and the region was divided into two parts according to the latitude of 50°S. The retrieval results were compared with the in-situ data. The Root Mean Square Error (RMS) and relativity of high latitude zone were 0.96℃and 0.93, respectively. And those of low latitude zone were 1.29 ℃ and 0.96, respectively. Artificial neural network (ANN) method was applied to retrieve Ta.The RMS and relativity were 1.26 ℃ and 0.98, respectively.

The response of phytoplankton to climate variability associated with the North Atlantic Oscillation

Inter-annual variability in the timing of phytoplankton spring bloom and phytoplankton community structure in the central North Atlantic Ocean was quantified using ocean color data and continuous plankton recorder (CPR) data. This variability was related to the North Atlantic Oscillation using correlation analysis and multivariate auto-regression models. The initiation of the spring bloom derived from CPR phytoplankton color index data is similar to that derived from satellite chlorophyll, and exhibits a nominal correlation with the sea surface temperature (SST) and the North Atlantic Oscillation (NAO). The extrapolated spring bloom timing suggested later initiation of blooms in the mid-1980s and earlier initiation of blooms in the 1990s. The climatological phytoplankton community structure in the central North Atlantic is dominated by diatoms, except for a shift in community composition favoring dinoflagellates in August. The ratio of diatoms to total phytoplankton abundance and the ratio of dinoflagellates to total phytoplankton abundance are both closely correlated with the NAO and SST. The extended time series of phytoplankton community structure between 1985 and 2009, deduced from the time series of SST and NAO over the same interval, showed a decadal shift away from diatoms towards dinoflagellates. The linkages between the NAO, and changes in stratification and phytoplankton processes occur over a larger scale than previously observed.

Joint distributions of waves and rain

The transfer of gases between the atmosphere and ocean is affected by a number of processes, of which wave action and rainfall are two of potential significance. Efforts have been made to quantify separately their contributions; however such assessments neglect the interaction of these phenomena. Here we look at the correlation statistics of waves and rain to note which regions display a strong association between rainfall and the local sea state. The conditional probability of rain varies from ~0.5% to ~15%, with most of the equatorial belt (which contains the ITCZ) showing a greater likelihood of rain at the lowest sea states. In contrast the occurrence of rain is independent of wave height in the Southern Ocean. The 1997/98 El Niño enhances the frequency of rain in some Pacific regions, with this change showing some association with wave conditions.

Transfer Across the Air-Sea Interface

The efficiency of transfer of gases and particles across the air-sea interface is controlled by several physical, biological and chemical processes in the atmosphere and water which are described here (including waves, large- and small-scale turbulence, bubbles, sea spray, rain and surface films). For a deeper understanding of relevant transport mechanisms, several models have been developed, ranging from conceptual models to numerical models. Most frequently the transfer is described by various functional dependencies of the wind speed, but more detailed descriptions need additional information. The study of gas transfer mechanisms uses a variety of experimental methods ranging from laboratory studies to carbon budgets, mass balance methods, micrometeorological techniques and thermographic techniques. Different methods resolve the transfer at different scales of time and space; this is important to take into account when comparing different results. Air-sea transfer is relevant in a wide range of applications, for example, local and regional fluxes, global models, remote sensing and computations of global inventories. The sensitivity of global models to the description of transfer velocity is limited; it is however likely that the formulations are more important when the resolution increases and other processes in models are improved. For global flux estimates using inventories or remote sensing products the accuracy of the transfer formulation as well as the accuracy of the wind field is crucial.

The Ocean Colour Climate Change Initiative: III. A round-robin comparison on in-water bio-optical algorithms

Satellite-derived remote-sensing reflectance (Rrs) can be used for mapping biogeochemically relevant variables, such as the chlorophyll concentration and the Inherent Optical Properties (IOPs) of the water, at global scale for use in climate-change studies. Prior to generating such products, suitable algorithms have to be selected that are appropriate for the purpose. Algorithm selection needs to account for both qualitative and quantitative requirements. In this paper we develop an objective methodology designed to rank the quantitative performance of a suite of bio-optical models. The objective classification is applied using the NASA bio-Optical Marine Algorithm Dataset (NOMAD). Using in situRrs as input to the models, the performance of eleven semi-analytical models, as well as five empirical chlorophyll algorithms and an empirical diffuse attenuation coefficient algorithm, is ranked for spectrally-resolved IOPs, chlorophyll concentration and the diffuse attenuation coefficient at 489 nm. The sensitivity of the objective classification and the uncertainty in the ranking are tested using a Monte-Carlo approach (bootstrapping). Results indicate that the performance of the semi-analytical models varies depending on the product and wavelength of interest. For chlorophyll retrieval, empirical algorithms perform better than semi-analytical models, in general. The performance of these empirical models reflects either their immunity to scale errors or instrument noise in Rrs data, or simply that the data used for model parameterisation were not independent of NOMAD. Nonetheless, uncertainty in the classification suggests that the performance of some semi-analytical algorithms at retrieving chlorophyll is comparable with the empirical algorithms. For phytoplankton absorption at 443 nm, some semi-analytical models also perform with similar accuracy to an empirical model. We discuss the potential biases, limitations and uncertainty in the approach, as well as additional qualitative considerations for algorithm selection for climate-change studies. Our classification has the potential to be routinely implemented, such that the performance of emerging algorithms can be compared with existing algorithms as they become available. In the long-term, such an approach will further aid algorithm development for ocean-colour studies.

Monitoring changes in phytoplankton abundance and composition in the Northwest Atlantic: a comparison of results obtained by Continuous Plankton Recorder sampling and colour satellite imagery

Phytoplankton abundance in the NW Atlantic was measured by continuous plankton recorder (CPR) sampling along tracks between Iceland and the western Scotian Shelf from 1998 to 2006, when sea-surface chlorophyll (SSChl) measurements were also being made by ocean colour satellite imagery using the SeaWiFS sensor. Seasonal and inter-annual changes in phytoplankton abundance were examined using data collected by both techniques, averaged over each of four shelf regions and four deep ocean regions. CPR sampling had gaps (missing months) in all regions and in the four deep ocean regions satellite observations were too sparse between November and February to be of use. Average seasonal cycles of SSChl were similar to those of total diatom abundance in seven regions, to those of the phytoplankton colour index in six regions, but were not similar to those of total dinoflagellate abundance anywhere. Large inter-annual changes in spring bloom dynamics were captured by both samplers in shelf regions. Changes in annual (or 8 months) averages of SSChl did not generally follow those of the CPR indices within regions and multi-year averages of SSChl, and the three CPR indices were generally higher in shelf than in deep ocean regions. Remote sensing and CPR sampling provide complementary ways of monitoring phytoplankton in the ocean: the former has superior temporal and spatial coverage and temporal resolution, and the latter provides better taxonomic information.

From silk to satellite: half a century of ocean colour anomalies in the Northeast Atlantic

Changes in phytoplankton dynamics influence marine biogeochemical cycles, climate processes, and food webs, with substantial social and economic consequences. Large-scale estimation of phytoplankton biomass was possible via ocean colour measurements from two remote sensing satellites – the Coastal Zone Color Scanner (CZCS, 1979-1986) and the Sea-viewing Wide Field-of-view Sensor (SeaWiFS, 1998-2010). Due to the large gap between the two satellite eras and differences in sensor characteristics, comparison of the absolute values retrieved from the two instruments remains challenging. Using a unique in situ ocean colour dataset that spans more than half a century, the two satellite-derived chlorophyll-a (Chl-a) eras are linked to assess concurrent changes in phytoplankton variability and bloom timing over the Northeast Atlantic Ocean and North Sea. Results from this unique re-analysis reflect a clear increasing pattern of Chl-a, a merging of the two seasonal phytoplankton blooms producing a longer growing season and higher seasonal biomass, since the mid-1980s. The broader climate plays a key role in Chl-a variability as the ocean colour anomalies parallel the oscillations of the Northern Hemisphere Temperature (NHT) since 1948.

Identifying four phytoplankton functional types from space: An ecological approach

Deriving maps of phytoplankton taxa based on remote sensing data using bio-optical properties of phytoplankton alone is challenging. A more holistic approach was developed using artificial neural networks, incorporating ecological and geographical knowledge together with ocean color, bio-optical characteristics, and remotely sensed physical parameters. Results show that the combined remote sensing approach could discriminate four major phytoplankton functional types (diatoms, dinoflagellates, coccolithophores, and silicoflagellates) with an accuracy of more than 70%. Models indicate that the most important information for phytoplankton functional type discrimination is spatio-temporal information and sea surface temperature. This approach can supply data for large-scale maps of predicted phytoplankton functional types, and an example is shown.