The Effect of ENSO on Phytoplankton Composition in the Pacific Ocean

The effect of climate variability on phytoplankton communities was assessed for the tropical and sub-tropical Pacific Ocean between 1998 and 2005 using an established biogeochemical assimilation model. The phytoplankton communities exhibited wide range of responses to climate variability, from radical shifts in the Equatorial Pacific, to changes of only a couple of phytoplankton groups in the North Central Pacific, to no significant changes in the South Pacific. In the Equatorial Pacific, climate variability dominated the variability of phytoplankton. Here, nitrate, chlorophyll and all but one of the 4 phytoplankton types (diatoms, cyanobacteria and coccolithophores) were strongly correlated (p less than 0.01) with the Multivariate El Nino Southern Oscillation Index (MEI). In the North Central Pacific, MEI and chlorophyll were significantly (p<0.01) correlated along with two of the phytoplankton groups (chlorophytes and coccolithophores). Ocean biology in the South Pacific was not significantly correlated with MEI. During La Ni a events, diatoms increased and expanded westward along the cold tongue (correlation with MEI, r=-0.81), while cyanobacteria concentrations decreased significantly (r=0.78). El Nino produced the reverse pattern, with cyanobacteria populations increasing while diatoms plummeted. The diverse response of phytoplankton in the different major basins of the Pacific suggests the different roles climate variability can play in ocean biology

Climate Variability and Phytoplankton in the Pacific Ocean

The effect of climate variability on phytoplankton communities was assessed for the tropical and sub-tropical Pacific Ocean between 1998 and 2005 using an established biogeochemical assimilation model. The phytoplankton communities exhibited wide range of responses to climate variability, from radical shifts in the Equatorial Pacific, to changes of only a couple of phytoplankton groups in the North Central Pacific, to no significant changes in the South Pacific. In the Equatorial Pacific, climate variability dominated the variability of phytoplankton. Here, nitrate, chlorophyll and all but one of the 4 phytoplankton types (diatoms, cyanobacteria and coccolithophores) were strongly correlated (p<0.01) with the Multivariate El Nino Southern Oscillation Index (MEI). In the North Central Pacific, MEI and chlorophyll were significantly (p<0.01) correlated along with two of the phytoplankton groups (chlorophytes and coccolithophores). Ocean biology in the South Pacific was not significantly correlated with MEI. During La Nina events, diatoms increased and expanded westward along the cold tongue (correlation with MEI, r=-0.81), while cyanobacteria concentrations decreased significantly (r=0.78). El Nino produced the reverse pattern, with cyanobacteria populations increasing while diatoms plummeted. The diverse response of phytoplankton in the different major basins of the Pacific suggests the different roles climate variability can play in ocean biology

Interannual Variation in Phytoplankton Primary Production at a Global Scale

We used the NASA Ocean Biogeochemical Model (NOBM) combined with remote sensing data via assimilation to evaluate the contribution of four phytoplankton groups to the total primary production. First, we assessed the contribution of each phytoplankton groups to the total primary production at a global scale for the period 1998-2011. Globally, diatoms contributed the most to the total phytoplankton production ((is)approximately 50%, the equivalent of 20 PgCy1). Coccolithophores and chlorophytes each contributed approximately 20% ((is) approximately 7 PgCy1) of the total primary production and cyanobacteria represented about 10% ((is) approximately 4 PgCy1) of the total primary production. Primary production by diatoms was highest in the high latitudes ((is) greater than 40 deg) and in major upwelling systems (Equatorial Pacific and Benguela system). We then assessed interannual variability of this group-specific primary production over the period 1998-2011. Globally the annual relative contribution of each phytoplankton groups to the total primary production varied by maximum 4% (1-2 PgCy1). We assessed the effects of climate variability on group-specific primary production using global (i.e., Multivariate El Nio Index, MEI) and "regional" climate indices (e.g., Southern Annular Mode (SAM), Pacific Decadal Oscillation (PDO) and North Atlantic Oscillation (NAO)). Most interannual variability occurred in the Equatorial Pacific and was associated with climate variability as indicated by significant correlation (p (is) less than 0.05) between the MEI and the group-specific primary production from all groups except coccolithophores. In the Atlantic, climate variability as indicated by NAO was significantly correlated to the primary production of 2 out of the 4 groups in the North Central Atlantic (diatoms/cyanobacteria) and in the North Atlantic (chlorophytes and coccolithophores). We found that climate variability as indicated by SAM had only a limited effect on group-specific primary production in the Southern Ocean. These results provide a modeling and data assimilation perspective to phytoplankton partitioning of primary production and contribute to our understanding of the dynamics of the carbon cycle in the oceans at a global scale

Forecasting Ocean Chlorophyll in the Equatorial Pacific

Using a global ocean biogeochemical model combined with a forecast of physical oceanic and atmospheric variables from the NASA Global Modeling and Assimilation Office, we assess the skill of a chlorophyll concentrations forecast in the Equatorial Pacific for the period 2012-2015 with a focus on the forecast of the onset of the 2015 El Nio event. Using a series of retrospective 9-month hindcasts, we assess the uncertainties of the forecasted chlorophyll by comparing the monthly total chlorophyll concentration from the forecast with the corresponding monthly ocean chlorophyll data from the Suomi-National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite (S-NPP VIIRS) satellite. The forecast was able to reproduce the phasing of the variability in chlorophyll concentration in the Equatorial Pacific, including the beginning of the 2015-2016 El Nio. The anomaly correlation coefficient (ACC) was significant (p less than 0.05) for forecast at 1-month (R=0.33), 8-month (R=0.42) and 9-month (R=0.41) lead times. The root mean square error (RMSE) increased from 0.0399 microgram ch1 L(exp -1) for the 1-month lead forecast to a maximum of 0.0472 microgram ch1 L(exp -1) for the 9-month lead forecast indicating that the forecast of the amplitude of chlorophyll concentration variability was getting worse. Forecasts with a 3-month lead time were on average the closest to the S-NPP VIIRS data (23% or 0.033 microgram ch1 L(exp -1)) while the forecast with a 9-month lead time were the furthest (31% or 0.042 microgram ch1 L(exp -1)). These results indicate the potential for forecasting chlorophyll concentration in this region but also highlights various deficiencies and suggestions for improvements to the current biogeochemical forecasting system. This system provides an initial basis for future applications including the effects of El Nio events on fisheries and other ocean resources given improvements identified in the analysis of these results

Interannual Variation in Phytoplankton Class-specific Primary Production at a Global Scale

Phytoplankton is responsible for over half of the net primary production on earth. The knowledge on the contribution of various phytoplankton groups to the total primary production is still poorly understood. Data from satellite observations suggest that for upwelling regions, photosynthetic rates by microplankton is higher than that of nanoplankton but that when the spatial extent is considered, the production by nanoplankton is comparable or even larger than microplankton. Here, we used the NASA Ocean Biogeochemical Model (NOBM) combined with remote sensing data via assimilation to evaluate the contribution of 4 phytoplankton groups to the total primary production. Globally, diatoms were the group that contributed the most to the total phytoplankton production (approx. 50%) followed by coccolithophores and chlorophytes. Primary production by diatoms was highest in high latitude (>45 deg) and in major upwelling systems (Equatorial Pacific and Benguela system). We assessed the effects of climate variability on the class-specific primary production using global (i.e. Multivariate El Nino Index, MEI) and 'regional' climate indices (e.g. Southern Annular Mode (SAM), Pacific Decadal Oscillation (PDO) and North Atlantic Oscillation (NAO)). Most interannual variability occurred in the Equatorial Pacific and was associated with climate variability. These results provide a modeling and data assimilation perspective to phytoplankton partitioning of primary production and contribute to our understanding of the dynamics of the carbon cycle in the oceans at a global scale

Sensitivity of Simulated Global Ocean Carbon Flux Estimates to Forcing by Reanalysis Products

Reanalysis products from MERRA, NCEP2, NCEP1, and ECMWF were used to force an established ocean biogeochemical model to estimate air-sea carbon fluxes (FCO2) and partial pressure of carbon dioxide (pCO2) in the global oceans. Global air-sea carbon fluxes and pCO2 were relatively insensitive to the choice of forcing reanalysis. All global FCO2 estimates from the model forced by the four different reanalyses were within 20% of in situ estimates (MERRA and NCEP1 were within 7%), and all models exhibited statistically significant positive correlations with in situ estimates across the 12 major oceanographic basins. Global pCO2 estimates were within 1% of in situ estimates with ECMWF being the outlier at 0.6%. Basin correlations were similar to FCO2. There were, however, substantial departures among basin estimates from the different reanalysis forcings. The high latitudes and tropics had the largest ranges in estimated fluxes among the reanalyses. Regional pCO2 differences among the reanalysis forcings were muted relative to the FCO2 results. No individual reanalysis was uniformly better or worse in the major oceanographic basins. The results provide information on the characterization of uncertainty in ocean carbon models due to choice of reanalysis forcing

Global Trends in Ocean Phytoplankton: A New Assessment Using Revised Ocean Colour Data

A recent revision of the NASA global ocean colour record shows changes in global ocean chlorophyll trends. This new 18-year time series now includes three global satellite sensors, the Sea-viewing Wide Field of view Sensor (SeaWiFS), Moderate Resolution Imaging Spectroradiometer (MODIS-Aqua), and Visible Infrared Imaging Radiometer Suite (VIIRS). The major changes are radiometric drift correction, a new algorithm for chlorophyll, and a new sensor VIIRS. The new satellite data record shows no significant trend in global annual median chlorophyll from 1998 to 2015, in contrast to a statistically significant negative trend from 1998 to 2012 in the previous version. When revised satellite data are assimilated into a global ocean biogeochemical model, no trend is observed in global annual median chlorophyll. This is consistent with previous findings for the 1998-2012 time period using the previous processing version and only two sensors (SeaWiFS and MODIS). Detecting trends in ocean chlorophyll with satellites is sensitive to data processing options and radiometric drift correction. The assimilation of these data, however, reduces sensitivity to algorithms and radiometry, as well as the addition of a new sensor. This suggests the assimilation model has skill in detecting trends in global ocean colour. Using the assimilation model, spatial distributions of significant trends for the 18-year record (1998-2015) show recent decadal changes. Most notable are the North and Equatorial Indian Oceans basins, which exhibit a striking decline in chlorophyll. It is exemplified by declines in diatoms and chlorophytes, which in the model are large and intermediate size phytoplankton. This decline is partially compensated by significant increases in cyanobacteria, which represent very small phytoplankton. This suggests the beginning of a shift in phytoplankton composition in these tropical and subtropical Indian basins

The Effects of Chlorophyll Assimilation on Carbon Fluxes in a Global Biogeochemical Model

In this paper, we investigated whether the assimilation of remotely-sensed chlorophyll data can improve the estimates of air-sea carbon dioxide fluxes (FCO2). Using a global, established biogeochemical model (NASA Ocean Biogeochemical Model, NOBM) for the period 2003-2010, we found that the global FCO2 values produced in the free-run and after assimilation were within -0.6 mol C m(sup -2) y(sup -1) of the observations. The effect of satellite chlorophyll assimilation was assessed in 12 major oceanographic regions. The region with the highest bias was the North Atlantic. Here the model underestimated the fluxes by 1.4 mol C m(sup -2) y(sup -1) whereas all the other regions were within 1 mol C m(sup -2) y(sup -1) of the data. The FCO2 values were not strongly impacted by the assimilation, and the uncertainty in FCO2 was not decreased, despite the decrease in the uncertainty in chlorophyll concentration. Chlorophyll concentrations were within approximately 25% of the database in 7 out of the 12 regions, and the assimilation improved the chlorophyll concentration in the regions with the highest bias by 10-20%. These results suggest that the assimilation of chlorophyll data does not considerably improve FCO2 estimates and that other components of the carbon cycle play a role that could further improve our FCO2 estimates

Technical Report Series on Global Modeling and Data Assimilation

MERRA products were used to force an established ocean biogeochemical model to estimate surface carbon inventories and fluxes in the global oceans. The results were compared to public archives of in situ carbon data and estimates. The model exhibited skill for ocean dissolved inorganic carbon (DIC), partial pressure of ocean CO2 (pCO2) and air-sea fluxes (FCO2). The MERRA-forced model produced global mean differences of 0.02% (approximately 0.3 microns) for DIC, -0.3% (about -1.2 (micro) atm; model lower) for pCO2, and -2.3% (-0.003 mol C/sq m/y) for FCO2 compared to in situ estimates. Basin-scale distributions were significantly correlated with observations for all three variables (r=0.97, 0.76, and 0.73, P<0.05, respectively for DIC, pCO2, and FCO2). All major oceanographic basins were represented as sources to the atmosphere or sinks in agreement with in situ estimates. However, there were substantial basin-scale and local departures

Assessing the Skill of Chlorophyll Forecasts: Latest Development and Challenges Ahead Using the Case of the Equatorial Pacific

Using a global ocean biogeochemical model combined with a forecast of physical oceanic and atmospheric variables from the NASA Global Modeling and Assimilation Office, we assess the skill of a chlorophyll concentrations forecast in the Equatorial Pacific for the period 2012-2015 with a focus on the forecast of the onset of the 2015 El Nino event. Using a series of retrospective 9-month hindcasts, we assess the uncertainties of the forecasted chlorophyll by comparing the monthly total chlorophyll concentration from the forecast with the corresponding monthly ocean chlorophyll data from the Suomi-National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite (S-NPP VIIRS) satellite. The forecast was able to reproduce the phasing of the variability in chlorophyll concentration in the Equatorial Pacific, including the beginning of the 2015-2016 El Nino. The anomaly correlation coefficient (ACC) was significant (p less than 0.05) for forecast at 1-month (R=0.33), 8-month (R=0.42) and 9-month (R=0.41) lead times. The root mean square error (RMSE) increased from 0.0399 microgram chl L(exp -1) for the 1-month lead forecast to a maximum of 0.0472 microgram chl L(exp -1) for the 9-month lead forecast indicating that the forecast of the amplitude of chlorophyll concentration variability was getting worse. Forecasts with a 3-month lead time were on average the closest to the S-NPP VIIRS data (23% or 0.033 microgram chl L(exp -1)) while the forecast with a 9-month lead time were the furthest (31% or 0.042 microgram chl L(exp -1)). These results indicate the potential for forecasting chlorophyll concentration in this region but also highlights various deficiencies and suggestions for improvements to the current biogeochemical forecasting system. This system provides an initial basis for future applications including the effects of El Nino events on fisheries and other ocean resources given improvements identified in the analysis of these results