Search

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.

The NASA PACE mission is a hyper-spectral radiometer planned for launch in the next decade. It is intended to provide new information on ocean biogeochemical constituents by parsing the details of high resolution spectral absorption and scattering. It is the first of its kind for global applications and as such, poses challenges for design and operation. To support pre-launch mission development and assess on-orbit capabilities, the NASA Global Modeling and Assimilation Office has developed a dynamic simulation of global water-leaving radiances, using an ocean model containing multiple ocean phytoplankton groups, particulate detritus, particulate inorganic carbon (PIC), and chromophoric dissolved organic carbon (CDOC) along with optical absorption and scattering processes at 1 nm spectral resolution. The purpose here is to assess the skill of the dynamic model and derived global radiances. Global bias, uncertainty, and correlation are derived using available modern satellite radiances at moderate spectral resolution. Total chlorophyll, PIC, and the absorption coefficient of CDOC (aCDOC), are simultaneously assimilated to improve the fidelity of the optical constituent fields. A 5-year simulation showed statistically significant (P < 0.05) comparisons of chlorophyll (r = 0.869), PIC (r = 0.868), and a CDOC (r =0.890) with satellite data. Additionally, diatoms (r = 0.890), cyanobacteria (r = 0.732), and coccolithophores (r = 0.716) were significantly correlated with in situ data. Global assimilated distributions of optical constituents were coupled with a radiative transfer model (Ocean-Atmosphere Spectral Irradiance Model, OASIM) to estimate normalized water-leaving radiances at 1 nm for the spectral range 250-800 nm. These unassimilated radiances were within 0.074 mW/sq cm/micron/sr of MODIS-Aqua radiances at 412, 443, 488, 531, 547, and 667 nm. This difference represented a bias of 10.4% (model low). A mean correlation of 0.706 (P < 0.05) was found with global distributions of MODIS radiances. These results suggest skill in the global assimilated model and resulting radiances. The reported error characterization suggests that the global dynamical simulation can support some aspects of mission design and analysis. For example, the high spectral resolution of the simulation supports investigations of band selection. The global nature of the radiance representations supports investigations of satellite observing scenarios. Global radiances at bands not available in current and past missions support investigations of mission capability. PACE, ocean color, water-leaving radiances, biogeochemical model, radiative transfer model

The importance of including directional and spectral light in simulations of ocean radiative transfer was investigated using a coupled biogeochemical-circulation-radiative model of the global oceans. The effort focused on phytoplankton abundances, nutrient concentrations and vertically-integrated net primary production. The importance was approached by sequentially removing directional (i.e., direct vs. diffuse) and spectral irradiance and comparing results of the above variables to a fully directionally and spectrally-resolved model. In each case the total irradiance was kept constant; it was only the pathways and spectral nature that were changed. Assuming all irradiance was diffuse had negligible effect on global ocean primary production. Global nitrate and total chlorophyll concentrations declined by about 20% each. The largest changes occurred in the tropics and sub-tropics rather than the high latitudes, where most of the irradiance is already diffuse. Disregarding spectral irradiance had effects that depended upon the choice of attenuation wavelength. The wavelength closest to the spectrally-resolved model, 500 nm, produced lower nitrate (19%) and chlorophyll (8%) and higher primary production (2%) than the spectral model. Phytoplankton relative abundances were very sensitive to the choice of non-spectral wavelength transmittance. The combined effects of neglecting both directional and spectral irradiance exacerbated the differences, despite using attenuation at 500 nm. Global nitrate decreased 33% and chlorophyll decreased 24%. Changes in phytoplankton community structure were considerable, representing a change from chlorophytes to cyanobacteria and coccolithophores. This suggested a shift in community function, from light-limitation to nutrient limitation: lower demands for nutrients from cyanobacteria and coccolithophores favored them over the more nutrient-demanding chlorophytes. Although diatoms have the highest nutrient demands in the model, their relative abundances were generally unaffected because they only prosper in nutrient-rich regions, such as the high latitudes and upwelling regions, which showed the fewest effects from the changes in radiative simulations. The results showed that including directional and spectral irradiance when simulating the ocean light field can be important for ocean biology, but the magnitude varies with variables and regions. The quantitative results are intended to assist ocean modelers when considering improved irradiance representations relative to other processes or variables associated with the issues of interest.

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.

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.

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. 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 were the group that contributed the most to the total phytoplankton production (50, the equivalent of 20 PgC y-1. Coccolithophores and chlorophytes each contributed to 20 (7 PgC y-1 of the total primary production and cyanobacteria represented about 10 (4 PgC y(sub-1) of the total primary production. Primary production by diatoms was highest in high latitude (45) 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 PgC y-1. We assessed the effects of climate variability on the class-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 0.05) between the MEI and the class-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 (diatomscyanobacteria) and in the North Atlantic (chlorophytes and coccolithophores). We found that climate variability as indicated by SAM had only a limited effect on the class-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.

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 PgC·y1). Coccolithophores and chlorophytes each contributed approximately 20% ((is) approximately 7 PgC·y1) of the total primary production and cyanobacteria represented about 10% ((is) approximately 4 PgC·y1) 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 PgC·y1). We assessed the effects of climate variability on group-specific primary production using global (i.e., Multivariate El Niño 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.

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.

Results from twin control simulations of the preindustrial CO2 gas exchange (natural flux of CO2) between the ocean and the atmosphere are presented here using the NASA-GISS climate model, in which the same atmospheric component (modelE2) is coupled to two different ocean models, the Russell ocean model and HYCOM. Both incarnations of the GISS climate model are also coupled to the same ocean biogeochemistry module (NOBM) which estimates prognostic distributions for biotic and abiotic fields that influence the air-sea flux of CO2. Model intercomparison is carried out at equilibrium conditions and model differences are contrasted with biases from present day climatologies. Although the models agree on the spatial patterns of the air-sea flux of CO2, they disagree on the strength of the North Atlantic and Southern Ocean sinks mainly because of kinematic (winds) and chemistry (pCO2) differences rather than thermodynamic (SST) ones. Biology/chemistry dissimilarities in the models stem from the different parameterizations of advective and diffusive processes, such as overturning, mixing and horizontal tracer advection and to a lesser degree from parameterizations of biogeochemical processes such as gravitational settling and sinking. The global meridional overturning circulation illustrates much of the different behavior of the biological pump in the two models, together with differences in mixed layer depth which are responsible for different SST, DIC and nutrient distributions in the two models and consequently different atmospheric feedbacks (in the wind, net heat and freshwater fluxes into the ocean).

Forward GEOS-5 AGCM simulations of CO2, with transport constrained by analyzed meteorology for 2009-2010, are examined. The CO2 distributions are evaluated using AIRS upper tropospheric CO2 and ACOS-GOSAT total column CO2 observations. Different combinations of surface C02 fluxes are used to generate ensembles of runs that span some uncertainty in surface emissions and uptake. The fluxes are specified in GEOS-5 from different inventories (fossil and biofuel), different data-constrained estimates of land biological emissions, and different data-constrained ocean-biology estimates. One set of fluxes is based on the established "Transcom" database and others are constructed using contemporary satellite observations to constrain land and ocean process models. Likewise, different approximations to sub-grid transport are employed, to construct an ensemble of CO2 distributions related to transport variability. This work is part of NASA's "Carbon Monitoring System Flux Pilot Project,"

A new empirical approach is developed for ocean color remote sensing. Called the Empirical Satellite Radiance-In situ Data (ESRID) algorithm, the approach uses relationships between satellite water-leaving radiances and in situ data after full processing, i.e., at Level-3, to improve estimates of surface variables while relaxing requirements on post-launch radiometric re-calibration. The approach is evaluated using SeaWiFS chlorophyll, which is the longest time series of the most widely used ocean color geophysical product. The results suggest that ESRID 1) drastically reduces the bias of ocean chlorophyll, most impressively in coastal regions, 2) modestly improves the uncertainty, and 3) reduces the sensitivity of global annual median chlorophyll to changes in radiometric re-calibration. Simulated calibration errors of 1% or less produce small changes in global median chlorophyll (less than 2.7%). In contrast, the standard NASA algorithm set is highly sensitive to radiometric calibration: similar 1% calibration errors produce changes in global median chlorophyll up to nearly 25%. We show that 0.1% radiometric calibration error (about 1% in water-leaving radiance) is needed to prevent radiometric calibration errors from changing global annual median chlorophyll more than the maximum interannual variability observed in the SeaWiFS 9-year record (+/- 3%), using the standard method. This is much more stringent than the goal for SeaWiFS of 5% uncertainty for water leaving radiance. The results suggest ocean color programs might consider less emphasis of expensive efforts to improve post-launch radiometric re-calibration in favor of increased efforts to characterize in situ observations of ocean surface geophysical products. Although the results here are focused on chlorophyll, in principle the approach described by ESRID can be applied to any surface variable potentially observable by visible remote sensing.

There is growing recognition that rigorous skill assessment is required to understand the ability of ocean biological models to represent ocean processes and distributions. Statistical analysis of model results with observations represents the most quantitative form of skill assessment, and this principle serves as well for data assimilation models. However, skill assessment for data assimilation requires special consideration. This is because there are three sets of information in the free-run model, data, and the assimilation model, which uses Data assimilation information from both the flee-run model and the data. Intercom parison of results among the three sets of information is important and useful for assessment, but is not conclusive since the three information sets are intertwined. An independent data set is necessary for an objective determination. Other useful measures of ocean biological data assimilation assessment include responses of unassimilated variables to the data assimilation, performance outside the prescribed region/time of interest, forecasting, and trend analysis. Examples of each approach from the literature are provided. A comprehensive list of ocean biological data assimilation and their applications of skill assessment, in both ecosystem/biogeochemical and fisheries efforts, is summarized.

Although modem ocean color sensors, such as MODIS and SeaWiFS are often considered global missions, in reality it takes many days, even months, to sample the ocean surface enough to provide complete global coverage. The irregular temporal sampling of ocean color sensors can produce biases in monthly and annual mean chlorophyll estimates. We quantified the biases due to sampling using data assimilation to create a "truth field", which we then sub-sampled using the observational patterns of MODIS and SeaWiFS. Monthly and annual mean chlorophyll estimates from these sub-sampled, incomplete daily fields were constructed and compared to monthly and annual means from the complete daily fields of the assimilation model, at a spatial resolution of 1.25deg longitude by 0.67deg latitude. The results showed that global annual mean biases were positive, reaching nearly 8% (MODIS) and >5% (SeaWiFS). For perspective the maximum interannual variability in the SeaWiFS chlorophyll record was about 3%. Annual mean sampling biases were low (<3%) in the midlatitudes (between -40deg and 40deg). Low interannual variability in the global annual mean sampling biases suggested that global scale trend analyses were valid. High latitude biases were much higher than the global annual means, up to 20% as a basin annual mean, and over 80% in some months. This was the result of the high solar zenith angle exclusion in the processing algorithms. Only data where the solar angle is <75deg are permitted, in contrast to the assimilation which samples regularly over the entire area and month. High solar zenith angles do not facilitate phytoplankton photosynthesis and consequently low chlorophyll concentrations occurring here are missed by the data sets. Ocean color sensors selectively sample in locations and times of favorable phytoplankton growth, producing overestimates of chlorophyll. The biases derived from lack of sampling in the high latitudes varied monthly, leading to artifacts in the apparent seasonal cycle from ocean color sensors. A false secondary peak in chlorophyll occurred in May-August, which resulted from lack of sampling in the Antarctic.

Assimilation of satellite ocean color data is a relatively new phenomenon in ocean sciences. However, with routine observations from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS), launched in late 1997, and now with new data from the Moderate Resolution Imaging Spectroradometer (MODIS) Aqua, there is increasing interest in ocean color data assimilation. Here SeaWiFS chlorophyll data were assimilated with an established thre-dimentional global ocean model. The assimilation improved estimates of hlorophyll and primary production relative to a free-run (no assimilation) model. This represents the first attempt at ocean color data assimilation using NASA satellites in a global model. The results suggest the potential of assimilation of satellite ocean chlorophyll data for improving models.

Results from the climatologically forced coupled ice/ocean/biogeochemical model that covers the Arctic and North Atlantic Oceans are presented and compared to the chlorophyll fields of satellite-derived ocean color measurements. Biogeochemical processes in the model are determined from the interactions among four phytoplankton functional groups (diatoms, chlorophytes, cyanobacteria and coccolithophores) and four nutrients (nitrate, ammonium, silicate and dissolved iron). The model simulates the general large-scale pattern in April, May and June, when compared to both satellite-derived and in situ observations. The subpolar North Atlantic was cool in the 1980s and warm in the latter 1990s, corresponding to the CZCS and SeaWiFS satellite observing periods, respectively. The oceanographic conditions during these periods resemble the typical subpolar upper ocean response to the NAO+ and NAO-phases, respectively. Thus, we use the atmospheric forcing composites from the two NAO phases to simulate the variability of the mid-ocean bloom during the satellite observing periods. The model results show that when the subpolar North Atlantic is cool, the NAO+ case, more nutrients are available in early spring than when the North Atlantic is warm, the NAO-case. However, the NAO+ simulation produces a later bloom than the NAO-simulation. This difference in the bloom times is also identified in SeaWiFS and CZCS satellite measurements. In the model results, we can trace the difference to the early diatom bloom due to a warmer upper ocean. The higher nutrient abundance in the NAO+ case did not provide larger total production than in the NAO- case, instead the two cases had a comparable area averaged amplitude. This leads us to conclude that in the subpolar North Atlantic, the timing of the spring phytoplankton bloom depends on surface temperature and the magnitude of the bloom is not significantly impacted by the nutrient abundance.

For both terrestrial and ocean carbon cycle science objectives, a hyperspectral geostationary sensor should enable the development of new remote sensing measurements for important but as yet unobservable variables, and with the overall goal of linking both terrestrial and ocean carbon cycle processes to climate variability. For terrestrial research, accurate estimates of carbon, water and energy (CWE) exchange between the terrestrial biosphere and atmosphere are needed to identify the geographical locations of carbon sources/sinks and to improve regional climate models and global climate change assessments. It is an enormous challenge to estimate CWE exchange from the infrequent temporal coverage provided by most polar-orbiting satellites, and without benefit of spectral indices that capture vegetation responses to stress conditions that down-regulate photosynthesis. Physiological status can be better assessed with spectral indices based on continuous, narrow (5 nm) bands, as can seasonal and annual terrestrial productivity. For coastal and ocean constituents, narrow-band observations in the ultraviolet and visible are essential to investigate the variability, dynamics and biogeochemical cycles of the world s coastal and open ocean regions, which will in turn help in measuring ocean productivity and predicting the variability of Ocean carbon uptake and its role in climate change scenarios. The GSFC Carbon Team has been pursuing a geostationary hyperspectral instrument, which would revolutionize our knowledge of biological processes on land, in the ocean, and along the coast by providing multiple, diurnal coverage. Preliminary studies in Goddard's Instrument Synthesis and Analysis Laboratory (ISAL) indicate that we can meet many of our science requirements: full spectral coverage (360-1000 nm); narrow bandwidths (5-10 nm); adequate ground resolution (100-200 m); and continental-scale coverage 4-6 times per day; all the while achieving a signal to noise ratio of between 500 and 1000 to 1. However, an innovative and bold focal plane design and a large mirror (1.8 meter diameter) would be required. The development of our science requirements and the results of the initial design study will be presented as well as our most recent technological developments.

For both terrestrial and ocean carbon cycle science objectives, high resolution (less than l0 nm) imaging spectrometers capable of acquiring multiple regional to global scale observations per day should enable the development of new remote sensing measurements for important but as yet unobservable variables, with the overall goal of linking both terrestrial and ocean carbon cycle processes to climate variability. For terrestrial research, accurate estimates of carbon, water and energy (CWE) exchange between the terrestrial biosphere and atmosphere a needed to id- the geographical locations and temporal dynamics of carbon sources/sinks and to improve regional climate models and climate change assessments. It is an enormous challenge to estimate CWE exchange from the infrequent temporal coverage and sparse spectral information provided by most single polar-orbiting, earth-looking satellite. The available satellite observations lack a sufficient number of well-placed narrow bands from which to derive spectral indices that capture vegetation responses to stress conditions associated with down-regulation of photosynthesis. Physiological status can best be assessed with spectral indices based on continuous, narrow bands in the visible/near infrared spectra, as can seasonal and annual terrestrial productivity. For coastal and ocean constituents, narrow-band observations in the ultraviolet and visible are essential to investigate the variability, dynamics and biogeochemical cycles of the world's coastal and open ocean regions, which will in turn help in measuring ocean productivity and predicting the variability of ocean carbon uptake and its role in climate change.

Oceanic phytoplankton may play an important role in the cycling of carbon on the Earth, through the uptake of carbon dioxide in the process of photosynthesis. Although they are ubiquitous in the global oceans, their abundances and dynamics are difficult to estimate, primarily due to the vast spatial extent of the oceans and the short time scales over which their abundances can change. Consequently, the effects of oceanic phytoplankton on biogeochemical cycling, climate change, and fisheries are not well known. In response to the potential importance of phytoplankton in the global carbon cycle and the lack of comprehensive data, NASA and the international community have established high priority satellite missions designed to acquire and produce high quality ocean color data (Table 6.1). Ten of the missions are routine global observational missions: the Ocean Color and Temperature Sensor (OCTS), the Polarization and Directionality of the Earth's Reflectances sensor (POLDER), Sea-viewing Wide Field-of-view Sensor (SeaWiFS), Moderate Resolution Imaging Spectrometer-AM (MODIS-AM), Medium Resolution Imaging Spectrometer (MERIS), Global Imager (GLI), MODIS-PM, Super-GLI (S-GLI), and the Visible/Infrared Imager and Radiometer Suite (VIIRS) on the NPOESS Preparatory Project (NPP) and the National Polar-orbiting Operational Environmental Satellite System (NPOESS). In addition, there are several other missions capable of providing ocean color data on smaller scales. Most of these missions contain the spectral band complement considered necessary to derive oceanic chlorophyll concentrations and other related parameters. Many contain additional bands that can provide important ancillary information about the optical and biological state of the oceans.

Satellite observations of global ocean chlorophyll from SeaWiFS were recently reprocessed to incorporate calibration and algorithm improvements. Here, comparisons are made between the newly reprocessed SeaWiFS Level-3 chlorophyll product and the previous version using in situ measurements. The results show that the newly reprocessed SeaWiFS data matches up better with the surface measurements than the previous version did. Globally, the slope of the match-ups improves to 0.85 from 0.78 in log-log scale. A significant trend that contributed to this improvement was the overall decrease in SeaWiFS chlorophyll levels less than 1.0 mg m(exp -3). Regional analyses reveal that the matchups improve in every oceanic basin, except the Antarctic. However, SeaWiFS continues to exhibit poor correspondence with in situ data in the North Atlantic where the match-ups have a slope of 0.54. Also, an examination of monthly images for May 1999 revealed that the number and magnitude of high-value chlorophyll pixels had increased in the high-latitude open ocean of the South Pacific.

An ocean-atmosphere radiative model (OARM) evaluates irradiance availability and quality in the water column to support phytoplankton growth and drive ocean thermodynamics. An atmospheric component incorporates spectral and directional effects of clear and cloudy skies as a function of atmospheric optical constituents, and spectral reflectance across the air-sea interface. An oceanic component evaluates the propagation of spectral and directional irradiance through the water column as a function of water, five phytoplankton groups, and chromophoric dissolved organic matter. It tracks the direct and diffuse streams from the atmospheric component, and a third stream, upwelling diffuse irradiance. The atmospheric component of OARM was compared to data sources at the ocean surface with a coefficient of determination (r2) of 0.97 and a root mean square of 12.1%.

We are midway into our 5th consecutive year of nearly continuous, high quality ocean color observations from space. The Ocean Color and Temperature Scanner/Polarization and Directionality of the Earth's Reflectances (OCTS/POLDER: Nov. 1996 - Jun. 1997), the Sea-viewing Wide Field-of-view Sensor (SeaWiFS: Sep. 1997 - present), and now the Moderate Resolution Imaging Spectrometer (MODIS: Sep. 2000 - present) have and are providing unprecedented views of chlorophyll dynamics on global scales. Global synoptic views of ocean chlorophyll were once a fantasy for ocean color scientists. It took nearly the entire 8-year lifetime of limited Coastal Zone Color Scanner (CZCS) observations to compile seasonal climatologies. Now SeaWIFS produces comparably complete fields in about 8 days. For the first time, scientists may observe spatial and temporal variability never before seen in a synoptic context. Even more exciting, we are beginning to plausibly ask questions of interannual variability. We stand at the beginning of long-time time series of ocean color, from which we may begin to ask questions of interdecadal variability and climate change. These are the scientific questions being addressed by users of the 18-year Advanced Very High Resolution Radiometer time series with respect to terrestrial processes and ocean temperatures. The nearly 5-year time series of ocean color observations now being constructed, with possibilities of continued observations, can put us at comparable standing with our terrestrial and physical oceanographic colleagues, and enable us to understand how ocean biological processes contribute to, and are affected by global climate change.

Oceanic phytoplankton may play an important role in the cycling of carbon on the Earth, through the uptake of carbon dioxide in the process of photosynthesis. Although they are ubiquitous in the global oceans, their abundances and dynamics are difficult to estimate, primarily due to the vast spatial extent of the oceans and the short time scales over which their abundances can change. Consequently, the effects of oceanic phytoplankton on biogeochemical cycling, climate change, and fisheries are not well known. In response to the potential importance of phytoplankton in the global carbon cycle and the lack of comprehensive data, the National Aeronautics and Space Administration (NASA) and the international community have established high priority satellite missions designed to acquire and produce high quality ocean color data. Seven of the missions are routine global observational missions: the Ocean Color and Temperature Sensor (OCTS), the Polarization and Directionality of the Earth's Reflectances sensor (POLDER), Sea-viewing Wide Field-of-view Sensor (SeaWiFS), Moderate Resolution Imaging Spectrometer-AM (MODIS-AM), Medium Resolution Imaging Spectrometer (MERIS), Global Imager (GLI), and MODIS-PM. In addition, there are several other missions capable of providing ocean color data on smaller scales. Most of these missions contain the spectral band complement considered necessary to derive oceanic pigment concentrations (i.e., phytoplankton abundance) and other related parameters. Many contain additional bands that can provide important ancillary information about the optical and biological state of the oceans. Any individual ocean color mission is limited in ocean coverage due to sun glint and clouds. For example, one of the first proposed missions, the SeaWiFS, can provide about 45% coverage of the global ocean in four days and only about 15% in one day.

Satellite observations of global ocean chlorophyll span over two decades. However, incompatibilities between processing algorithms prevent us from quantifying natural variability. We applied a comprehensive reanalysis to the Coastal Zone Color Scanner (CZCS) archive, called the NOAA-NASA CZCS Reanalysis (NCR) Effort. NCR consisted of 1) algorithm improvement (AI), where CZCS processing algorithms were improved using modernized atmospheric correction and bio-optical algorithms, and 2) blending, where in situ data were incorporated into the CZCS AI to minimize residual errors. The results indicated major improvement over the previously available CZCS archive. Global spatial and seasonal patterns of NCR chlorophyll indicated remarkable correspondence with modern sensors, suggesting compatibility. The NCR permits quantitative analyses of interannual and interdecadal trends in global ocean chlorophyll.

The global ocean chlorophyll archive produced by the Coastal Zone Color Scanner (CZCS) was revised using compatible algorithms with the Sea-viewing Wide Field-of-view Sensor (SeaWIFS), and both were blended with in situ data. This methodology permitted a quantitative comparison of decadal changes in global ocean chlorophyll from the CZCS (1979-1986) and SeaWiFS (Sep. 1997-Dec. 2000) records. Global seasonal means of ocean chlorophyll decreased over the two observational segments, by 8% in winter to 16% in autumn. Chlorophyll in the high latitudes was responsible for most of the decadal change. Conversely, chlorophyll concentrations in the low latitudes increased. The differences and similarities of the two data records provide evidence of how the Earth's climate may be changing and how ocean biota respond. Furthermore, the results have implications for the ocean carbon cycle.

A coupled ocean general circulation, biogeochemical, and radiative model was constructed to evaluate and understand the nature of seasonal variability of chlorophyll and nutrients in the global oceans. Biogeochemical processes in the model are determined from the influences of circulation and turbulence dynamics, irradiance availability. and the interactions among three functional phytoplankton groups (diatoms. chlorophytes, and picoplankton) and three nutrients (nitrate, ammonium, and silicate). Basin scale (greater than 1000 km) model chlorophyll results are in overall agreement with CZCS pigments in many global regions. Seasonal variability observed in the CZCS is also represented in the model. Synoptic scale (100-1000 km) comparisons of imagery are generally in conformance although occasional departures are apparent. Model nitrate distributions agree with in situ data, including seasonal dynamics, except for the equatorial Atlantic. The overall agreement of the model with satellite and in situ data sources indicates that the model dynamics offer a reasonably realistic simulation of phytoplankton and nutrient dynamics on synoptic scales. This is especially true given that initial conditions are homogenous chlorophyll fields. The success of the model in producing a reasonable representation of chlorophyll and nutrient distributions and seasonal variability in the global oceans is attributed to the application of a generalized, processes-driven approach as opposed to regional parameterization and the existence of multiple phytoplankton groups with different physiological and physical properties. These factors enable the model to simultaneously represent many aspects of the great diversity of physical, biological, chemical, and radiative environments encountered in the global oceans.

A coupled ocean general circulation, biogeochemical, and radiative model was constructed to evaluate and understand the nature of seasonal variability of chlorophyll and nutrients in the global oceans. Biogeochemical processes in the model were determined from the influences of circulation and turbulence dynamics, irradiance availability, and the interactions among three functional phytoplankton groups (diatoms, chlorophytes, and picoplankton) and three nutrients (nitrate, ammonium, and silicate). Basin scale (>1000 km) model chlorophyll seasonal distributions were statistically positively correlated with CZCS chlorophyll in 10 of 12 major oceanographic regions, and with SeaWiFS in all 12. Notable disparities in magnitudes occurred, however, in the tropical Pacific, the spring/summer bloom in the Antarctic, autumn in the northern high latitudes, and during the southwest monsoon in the North Indian Ocean. Synoptic scale (100-1000 km) comparisons of satellite and in situ data exhibited broad agreement, although occasional departures were apparent. Model nitrate distributions agreed with in situ data, including seasonal dynamics, except for the equatorial Atlantic. The overall agreement of the model with satellite and in situ data sources indicated that the model dynamics offer a reasonably realistic simulation of phytoplankton and nutrient dynamics on basin and synoptic scales.

A coupled general ocean circulation, biogeochemical, and radiative model was constructed to evaluate and understand the nature of seasonal variability of chlorophyll and nutrients in the global oceans. The model is driven by climatological meteorological conditions, cloud cover, and sea surface temperature. Biogeochemical processes in the model are determined from the influences of circulation and turbulence dynamics, irradiance availability, and the interactions among three functional phytoplankton groups (diatoms, chorophytes, and picoplankton) and three nutrient groups (nitrate, ammonium, and silicate). Phytoplankton groups are initialized as homogeneous fields horizontally and vertically, and allowed to distribute themselves according to the prevailing conditions. Basin-scale model chlorophyll results are in very good agreement with CZCS pigments in virtually every global region. Seasonal variability observed in the CZCS is also well represented in the model. Synoptic scale (100-1000 km) comparisons of imagery are also in good conformance, although occasional departures are apparent. Agreement of nitrate distributions with in situ data is even better, including seasonal dynamics, except for the equatorial Atlantic. The good agreement of the model with satellite and in situ data sources indicates that the model dynamics realistically simulate phytoplankton and nutrient dynamics on synoptic scales. This is especially true given that initial conditions are homogenous chlorophyll fields. The success of the model in producing a reasonable representation of chlorophyll and nutrient distributions and seasonal variability in the global oceans is attributed to the application of a generalized, processes-driven approach as opposed to regional parameterization, and the existence of multiple phytoplankton groups with different physiological and physical properties. These factors enable the model to simultaneously represent the great diversity of physical, biological, chemical, and radiative environments encountered in the global oceans.

The historical archives of in situ (National Oceanographic Data Center) and satellite (Coastal Zone Color Scanner) chlorophyll data were combined using the blended analysis method of Reynolds [1988] in an attempt to construct an improved climatological seasonal representation of global chlorophyll distributions. The results of the blended analysis differed dramatically from the CZCS representation: global chlorophyll estimates increased 8-35% in the blended analysis depending upon season. Regional differences were even larger, up to 140% in the equatorial Indian Ocean in summer (during the southwest monsoon). Tropical Pacific chlorophyll values increased 25-41%. The results suggested that the CZCS generally underestimates chlorophyll. Regional and seasonal differences in the blended analysis were sufficiently large as to produce a different representation of global chlorophyll distributions than otherwise inferred from CZCS data alone. Analyses of primary production and biogeochemical cycles may be substantially impacted by these results.

Ocean coverages of SeaWiFS and MODIS were assessed for three seasons by considering monthly mean values of surface winds speeds and cloud cover. Mean and maximum coverages combined SeaWiFS and MODIS by considering combined coverages for ten-degree increments of the MODIS orbital mean anomaly. From this analysis the mean and maximum combined coverages for SeaWiFS and MODIS were determined for one and four-day periods for spring, summer, and winter seasons. Loss of coverage due to Sun glint and cloud cover were identified for both the individual and combined cases. Our analyses indicate that MODIS will enhance ocean coverage for all three seasons examined. ne combined SeaWiFS/MODIS show an increase of coverage of 42.2% to 48.7% over SeaWiFS alone for the three seasons studied; the increase in maximum one day coverage ranges from 47.5% to 52.0%. The increase in four-day coverage for the combined case ranged from 31.0% to 35.8% for mean coverage and 33.1 % to 39.2% for maximum coverage. We computed meridional distributions of coverages by binning the data into five-degree latitude bands. Our analysis shows a strong seasonal dependence of coverage. In general the meridional analysis indicates that increase in coverages for SeaWiFS/MODIS over SeaWiFS alone are greatest near the solar declination.