Global database of surface ocean particulate organic carbon export fluxes diagnosed from the 234Th technique

Through the processes of the biological pump, carbon is exported to the deep ocean in the form of dissolved and particulate organic matter. There are several ways by which downward export fluxes can be estimated. The great attraction of the 234Th technique is that its fundamental operation allows a downward flux rate to be determined from a single water column profile of thorium coupled to an estimate of POC/234Th ratio in sinking matter. We present a database of 723 estimates of organic carbon export from the surface ocean derived from the 234Th technique. Data were collected from tables in papers published between 1985 and 2013 only. We also present sampling dates, publication dates and sampling areas. Most of the open ocean Longhurst provinces are represented by several measurements. However, the Western Pacific, the Atlantic Arctic, South Pacific and the South Indian Ocean are not well represented. There is a variety of integration depths ranging from surface to 220m. Globally the fluxes ranged from -22 to 125 mmol of C/m**2/d. We believe that this database is important for providing new global estimate of the magnitude of the biological carbon pump.

(Table 1) Core top Globorotalia truncatulinoides d18O ratios of Atlantic Ocean samples

We investigate whether the oxygen isotope ratio in the test of Globorotalia truncatulinoides can serve as a proxy for intermediate depth (200-500 m) density. Since intermediate depth horizontal density gradients are associated with the vertical shear of upper ocean flows, this proxy could provide a tool for reconstructing past ocean circulation. The spatial pattern of core top Gr. truncatulinoides d18O in the Atlantic Ocean mimics the upper ocean density gradients associated with the major ocean currents. To better constrain the controls on the calcification depth(s) of Gr. truncatulinoides, we attempt to simulate the surface sediment data set using water column temperature and salinity conditions above the core sites. We predicted foraminiferal d18O for each core site assuming (1) the calcification occurs at a single depth and (2) the initial calcification is at the surface and the subsequent calcification is at 800 m water depth. The predicted d18O best resembles measured d18O of Gr. truncatulinoides when using (1) a single depth calcification at 350 m or (2) a two-depth approximation with 30% surface and 70% 800-m calcification. This result gives us confidence in the ability of d18O in Gr. truncatulinoides to proxy lateral density gradients at the intermediate depths associated with upper ocean flow.

Studies on ocean surface layer responses to atmospheric forcing in the North Indian Ocean

The present study on upper ocean responses to atmospheric forcing (associated with cyclone passage) in North Indian Ocean revealed significant variability between AS and BoB. The analysis of cyclone frequency during 1947 to 2006 exhibited lesser frequency of cyclones in AS than that of BoB. The analysis also revealed significant reduction in cyclone frequency after the year 1976 with substantial reduction during monsoon season. The long term SST data at selected points in AS and BoB could not reveal any relation with reduction in cyclone frequency. However the SLP at same locations exhibited considerable increase during mid 1970’s, which could have contributed to the observed reduction in cyclone frequency after the year 1976.The response in waves during cyclone passage exhibited significant asymmetry on either side of the track in AS and BoB and the response is observed at 100’s of kilometers away from the track. The significant clockwise rotation in wave direction is observed on the right side of the track starting from near the track to far away locations, which existed for a longer duration. However, the anticlockwise rotation in wave direction is observed over a shorter distance on the left side of the track and dissipated immediately.Inertial oscillation is observed in surface current and in the mixed layer temperature associated with cyclone passage, which revealed the role of relative location(s) on either side of the track. The inertial peak closer to the local inertial period indicates maximum transfer of energy during the cyclone passage in both AS and BoB. The absence of strong inertial oscillation even with clockwise rotation in surface current and wind indicates the dominant role of duration of strong wind in generating inertial oscillation.The oceanic response associated with cyclone passage reveal the variable response(s) which depends on cyclone intensity, the proximity to track and cyclone translation speed. It is observed that resonance with wind generates higher response in surface current, wave and SST on the right side of the track and it lasts for a longer duration. The maximum oceanic response is observed at a few kilometers away on right side of the track. However lesser rightward bias in the location of maximum cooling is observed for cyclones with low cyclone translation speed. The response on the left side of the track is less and is limited over a shorter distance and dissipates immediately. It is observed that the ocean response, in general, increases with intensity of cyclones. However the differential cooling produced by the same intensity cyclones in AS and in BoB indicates the dominant role of low cyclone translation speed in oceanic response.The surface cooling exhibited strikingly differential responses between AS and BoB. The TMI-SST and buoy observations exhibited significant cooling for a longer duration in AS compared to that of BoB. The spatial extent of cooling is also much higher in AS than that of BoB. The wide spread cooling associated with cyclone passage in AS indicates the dominant role of thermal structure in oceanic response in AS than that of BoB.

Historical reconstruction of the Atlantic Meridional Overturning Circulation from the ECMWF operational ocean reanalysis

A reconstruction of the Atlantic Meridional Overturning Circulation (MOC) for the period 1959–2006 has been derived from the ECMWF operational ocean reanalysis. The reconstruction shows a wide range of time-variability, including a downward trend. At 26N, both the MOC intensity and changes in its vertical structure are in good agreement with previous estimates based on trans-Atlantic surveys. At 50N, the MOC and strength of the subpolar gyre are correlated at interannual time scales, but show opposite secular trends. Heat transport variability is highly correlated with the MOC but shows a smaller trend due to the warming of the upper ocean, which partially compensates for the weakening of the circulation. Results from sensitivity experiments show that although the time-varying upper boundary forcing provides useful MOC information, the sequential assimilation of ocean data further improves the MOC estimation by increasing both the mean and the time variability.

Isolating the signal of ocean global warming

Identifying the signature of global warming in the world's oceans is challenging because low frequency circulation changes can dominate local temperature changes. The IPCC fourth assessment reported an average ocean heating rate of 0.21 ± 0.04 Wm−2 over the period 1961–2003, with considerable spatial, interannual and inter-decadal variability. We present a new analysis of millions of ocean temperature profiles designed to filter out local dynamical changes to give a more consistent view of the underlying warming. Time series of temperature anomaly for all waters warmer than 14°C show large reductions in interannual to inter-decadal variability and a more spatially uniform upper ocean warming trend (0.12 Wm−2 on average) than previous results. This new measure of ocean warming is also more robust to some sources of error in the ocean observing system. Our new analysis provides a useful addition for evaluation of coupled climate models, to the traditional fixed depth analyses.

Coupled simulations of the Indian monsoon intraseasonal oscillation using a fine-resolution mixed-layer ocean model

The intraseasonal variability (ISV) of the Indian summer monsoon is dominated by a 30–50 day oscillation between “active” and “break” events of enhanced and reduced rainfall over the subcontinent, respectively. These organized convective events form in the equatorial Indian Ocean and propagate north to India. Atmosphere–ocean coupled processes are thought to play a key role the intensity and propagation of these events. A high-resolution, coupled atmosphere–mixed-layer-oceanmodel is assembled: HadKPP. HadKPP comprises the Hadley Centre Atmospheric Model (HadAM3) and the K Profile Parameterization (KPP) mixed-layer ocean model. Following studies that upper-ocean vertical resolution and sub-diurnal coupling frequencies improve the simulation of ISV in SSTs, KPP is run at 1 m vertical resolution near the surface; the atmosphere and ocean are coupled every three hours. HadKPP accurately simulates the 30–50 day ISV in rainfall and SSTs over India and the Bay of Bengal, respectively, but suffers from low ISV on the equator. This is due to the HadAM3 convection scheme producing limited ISV in surface fluxes. HadKPP demonstrates little of the observed northward propagation of intraseasonal events, producing instead a standing oscillation. The lack of equatorial ISV in convection in HadAM3 constrains the ability of KPP to produce equatorial SST anomalies, which further weakens the ISV of convection. It is concluded that while atmosphere–ocean interactions are undoubtedly essential to an accurate simulation of ISV, they are not a panacea for model deficiencies. In regions where the atmospheric forcing is adequate, such as the Bay of Bengal, KPP produces SST anomalies that are comparable to the Tropical Rainfall Measuring Mission Microwave Imager (TMI) SST analyses in both their magnitude and their timing with respect to rainfall anomalies over India. HadKPP also displays a much-improved phase relationship between rainfall and SSTs over a HadAM3 ensemble forced by observed SSTs, when both are compared to observations. Coupling to mixed-layer models such as KPP has the potential to improve operational predictions of ISV, particularly when the persistence time of SST anomalies is shorter than the forecast lead time.

The impact of global warming on the tropical Pacific Ocean and El Niño

The El Niño–Southern Oscillation (ENSO) is a naturally occurring fluctuation that originates in the tropical Pacific region and affects ecosystems, agriculture, freshwater supplies, hurricanes and other severe weather events worldwide. Under the influence of global warming, the mean climate of the Pacific region will probably undergo significant changes. The tropical easterly trade winds are expected to weaken; surface ocean temperatures are expected to warm fastest near the equator and more slowly farther away; the equatorial thermocline that marks the transition between the wind-mixed upper ocean and deeper layers is expected to shoal; and the temperature gradients across the thermocline are expected to become steeper. Year-to-year ENSO variability is controlled by a delicate balance of amplifying and damping feedbacks, and one or more of the physical processes that are responsible for determining the characteristics of ENSO will probably be modified by climate change. Therefore, despite considerable progress in our understanding of the impact of climate change on many of the processes that contribute to El Niño variability, it is not yet possible to say whether ENSO activity will be enhanced or damped, or if the frequency of events will change.

Assimilation of RAPID array observations into an ocean model

Experiments assimilating the RAPID dataset of deep temperature and salinity profiles at 26.5°N on the western and eastern Atlantic boundaries into a 1° global NEMO ocean model have been performed. The meridional overturning circulation (MOC) is then assessed against the transports calculated directly from observations. The best initialization found for this short period was obtained by assimilating the EN3 upper-ocean hydrography database prior to 2004, after which different methods of assimilating 5-day average RAPID profiles at the western boundary were tested. The model MOC is strengthened by ∼ 2 Sv giving closer agreement with the RAPID array transports, when the western boundary profiles are assimilated only below 900 m (the approximate depth of the Florida Straits, which are not well resolved) and when the T,S observations are spread meridionally from 10 to 35°N along the deep western boundary. The use of boundary-focused covariances has the largest impact on the assimilation results, otherwise using more conventional Gaussian covariances has a very local impact on the MOC at 26°N with strong adverse impacts on the MOC stream function at higher and lower latitudes. Even using boundary-focused covariances only enables the MOC to be strengthened for ∼ 2 years, after which the increased transport of warm waters leads to a negative feedback on water formation in the subpolar gyre which then reduces the MOC. This negative feedback can be mitigated if EN3 hydrography data continue to be assimilated along with the RAPID array boundary data. Copyright © 2012 Royal Meteorological Society and Crown in the right of Canada.

Atmosphere-ocean coupled processes in the Madden-Julian oscillation

The Madden-Julian oscillation (MJO) is a convectively coupled 30-70 day (intraseasonal) tropical atmospheric mode that drives variations in global weather, but which is poorly simulated in most atmospheric general circulation models. Over the past two decades, field campaigns and modeling experiments have suggested that tropical atmosphere-ocean interactions may sustain or amplify the pattern of enhanced and suppressed atmospheric convection that defines the MJO, and encourage its eastward propagation through the Indian and Pacific Oceans. New observations collected during the past decade have advanced our understand of the ocean response to atmospheric MJO forcing and the resulting intraseasonal sea surface temperature (SST) fluctuations. Numerous modeling studies have revealed a considerable impact of the mean state on MJO ocean-atmosphere coupled processes, as well as the importance of resolving the diurnal cycle of atmosphere--upper-ocean interactions. New diagnostic methods provide insight to atmospheric variability and physical processes associated with the MJO, but offer limited insight on the role of ocean feedbacks. Consequently, uncertainty remains concerning the role of the ocean in MJO theory. Our understanding of how atmosphere-ocean coupled processes affect the MJO can be improved by collecting observations in poorly sampled regions of MJO activity, assessing oceanic and atmospheric drivers of surface fluxes, improving the representation of upper-ocean mixing in coupled-model simulations, designing model experiments that minimize mean-state differences, and developing diagnostic tools to evaluate the nature and role of coupled ocean-atmosphere processes over the MJO cycle.

Balanced Ocean-Data Assimilation near the Equator

The question is addressed whether using unbalanced updates in ocean-data assimilation schemes for seasonal forecasting systems can result in a relatively poor simulation of zonal currents. An assimilation scheme, where temperature observations are used for updating only the density field, is compared to a scheme where updates of density field and zonal velocities are related by geostrophic balance. This is done for an equatorial linear shallow-water model. It is found that equatorial zonal velocities can be detoriated if velocity is not updated in the assimilation procedure. Adding balanced updates to the zonal velocity is shown to be a simple remedy for the shallow-water model. Next, optimal interpolation (OI) schemes with balanced updates of the zonal velocity are implemented in two ocean general circulation models. First tests indicate a beneficial impact on equatorial upper-ocean zonal currents.

Intercomparison and validation of the mixed layer depth fields of global ocean syntheses

Intercomparison and evaluation of the global ocean surface mixed layer depth (MLD) fields estimated from a suite of major ocean syntheses are conducted. Compared with the reference MLDs calculated from individual profiles, MLDs calculated from monthly mean and gridded profiles show negative biases of 10–20 m in early spring related to the re-stratification process of relatively deep mixed layers. Vertical resolution of profiles also influences the MLD estimation. MLDs are underestimated by approximately 5–7 (14–16) m with the vertical resolution of 25 (50) m when the criterion of potential density exceeding the 10-m value by 0.03 kg m−3 is used for the MLD estimation. Using the larger criterion (0.125 kg m−3) generally reduces the underestimations. In addition, positive biases greater than 100 m are found in wintertime subpolar regions when MLD criteria based on temperature are used. Biases of the reanalyses are due to both model errors and errors related to differences between the assimilation methods. The result shows that these errors are partially cancelled out through the ensemble averaging. Moreover, the bias in the ensemble mean field of the reanalyses is smaller than in the observation-only analyses. This is largely attributed to comparably higher resolutions of the reanalyses. The robust reproduction of both the seasonal cycle and interannual variability by the ensemble mean of the reanalyses indicates a great potential of the ensemble mean MLD field for investigating and monitoring upper ocean processes.

The Ocean Reanalysis Intercomparison project (ORA-IP)

Uncertainty in ocean analysis methods and deficiencies in the observing system are major obstacles for the reliable reconstruction of the past ocean climate. The variety of existing ocean reanalyses is exploited in a multi-reanalysis ensemble to improve the ocean state estimation and to gauge uncertainty levels. The ensemble-based analysis of signal-to-noise ratio allows the identification of ocean characteristics for which the estimation is robust (such as tropical mixed-layer-depth, upper ocean heat content), and where large uncertainty exists (deep ocean, Southern Ocean, sea ice thickness, salinity), providing guidance for future enhancement of the observing and data assimilation systems.

Seasonality of submesoscale flows in the ocean surface boundary layer

A signature of submesoscale flows in the upper ocean is skewness in the distribution of relative vorticity. Expected to result for high Rossby-number flows, such skewness has implications for mixing, dissipation and stratification within the upper ocean. An array of moorings deployed in the Northeast Atlantic for one year as part of the OSMOSIS experiment reveals that relative vorticity is positively skewed during winter even though the scale of the Rossby number is less than 0.5. Furthermore, this skewness is reduced to zero during spring and autumn. There is also evidence of modest seasonal variations in the gradient Rossby number. The proposed mechanism by which relative vorticity is skewed is that the ratio of lateral to vertical buoyancy gradients, as summarized by the inverse gradient Richardson number, restricts its range during winter but less so at other times of the year. These results support recent observations and model simulations suggesting the upper ocean is host to a seasonal cycle in submesoscale turbulence.

Mean and Variability of the Tropical Atlantic Ocean in the CCSM4

This study analyzes important aspects of the tropical Atlantic Ocean from simulations of the fourth version of the Community Climate System Model (CCSM4): the mean sea surface temperature (SST) and wind stress, the Atlantic warm pools, the principal modes of SST variability, and the heat budget in the Benguela region. The main goal was to assess the similarities and differences between the CCSM4 simulations and observations. The results indicate that the tropical Atlantic overall is realistic in CCSM4. However, there are still significant biases in the CCSM4 Atlantic SSTs, with a colder tropical North Atlantic and a hotter tropical South Atlantic, that are related to biases in the wind stress. These are also reflected in the Atlantic warm pools in April and September, with its volume greater than in observations in April and smaller than in observations in September. The variability of SSTs in the tropical Atlantic is well represented in CCSM4. However, in the equatorial and tropical South Atlantic regions, CCSM4 has two distinct modes of variability, in contrast to observed behavior. A model heat budget analysis of the Benguela region indicates that the variability of the upper-ocean temperature is dominated by vertical advection, followed by meridional advection.

Database of diazotrophs in global ocean: abundances, biomass and nitrogen fixation rates

[EN] Marine N2 fixing microorganisms, termed diazotrophs, are a key functional group in marine pelagic ecosystems. The biological fixation of dinitrogen (N2) to bioavailable nitrogen provides an important new source of nitrogen for pelagic marine ecosystems 5 and influences primary productivity and organic matter export to the deep ocean. As one of a series of efforts to collect biomass and rates specific to different phytoplankton functional groups, we have constructed a database on diazotrophic organisms in the global pelagic upper ocean by compiling about 12 000 direct field measurements of cyanobacterial diazotroph abundances (based on microscopic cell counts or qPCR 10 assays targeting the nifH genes) and N2 fixation rates. Biomass conversion factors are estimated based on cell sizes to convert  abundance data to diazotrophic biomass. The database is limited spatially, lacking large regions of the ocean especially in the Indian Ocean. The data are approximately log-normal distributed, and large variances exist in most sub-databases with non-zero values differing 5 to 8 orders of magnitude. 15 Lower mean N2 fixation rate was found in the North Atlantic Ocean than the Pacific Ocean. Reporting the geometric mean and the range of one geometric standard error below and above the geometric mean, the pelagic N2 fixation rate in the global ocean is estimated to be 62 (53–73) TgNyr−1 and the pelagic diazotrophic biomass in the global ocean is estimated to be 4.7 (2.3–9.6) TgC from cell counts and to 89 (40–20 200) TgC from nifH-based abundances. Uncertainties related to biomass conversion factors can change the estimate of geometric mean pelagic diazotrophic biomass in the global ocean by about ±70 %. This evolving database can be used to study spatial and temporal distributions and variations of marine N2 fixation, to validate geochemical estimates and to parameterize and validate biogeochemical models. The database is 25 stored in PANGAEA (http://doi.pangaea.de/10.1594/PANGAEA.774851).