Reevaluation of historical ocean heat content variations with time-varying XBT and MBT depth bias corrections

As reported in former studies, temperature observations obtained by expendable bathythermographs (XBTs) and mechanical bathythermographs (MBTs) appear to have positive biases as much as they affect major climate signals. These biases have not been fully taken into account in previous ocean temperature analyses, which have been widely used to detect global warming signals in the oceans. This report proposes a methodology for directly eliminating the biases from the XBT and MBT observations. In the case of XBT observation, assuming that the positive temperature biases mainly originate from greater depths given by conventional XBT fall-rate equations than the truth, a depth bias equation is constructed by fitting depth differences between XBT data and more accurate oceanographic observations to a linear equation of elapsed time. Such depth bias equations are introduced separately for each year and for each probe type. Uncertainty in the gradient of the linear equation is evaluated using a non-parametric test. The typical depth bias is +10 m at 700 m depth on average, which is probably caused by various indeterminable sources of error in the XBT observations as well as a lack of representativeness in the fall-rate equations adopted so far. Depth biases in MBT are fitted to quadratic equations of depth in a similar manner to the XBT method. Correcting the historical XBT and MBT depth biases by these equations allows a historical ocean temperature analysis to be conducted. In comparison with the previous temperature analysis, large differences are found in the present analysis as follows: the duration of large ocean heat content in the 1970s shortens dramatically, and recent ocean cooling becomes insignificant. The result is also in better agreement with tide gauge observations. On leave from the Meteorological Research Institute of the Japan Meteorological Agency.     

Meridional heat transport determined with expandable bathythermographs—Part II: South Atlantic transport

Fourteen temperature sections collected between July 2002 and May 2006 are analyzed to obtain estimates of the meridional heat transport variability of the South Atlantic Ocean. The methodology proposed in Part I is used to calculate the heat transport from temperature data obtained from high-density XBT profiles taken along transects from Cape Town, South Africa to Buenos Aires, Argentina. Salinity is estimated from Argo profiles and CTD casts for each XBT temperature observation using statistical relationships between temperature, latitude, longitude, and salinity computed along constant-depth surfaces. Full-depth temperature/salinity profiles are obtained by extending the profiles to the bottom of the ocean using deep climatological data. The meridional transport is then determined by using the standard geostrophic method, applying NCEP-derived Ekman transports, and requiring that salt flux through the Bering Straits be conserved. The results from the analysis indicate a mean meridional heat transport of 0.54 PW (PW=10¹⁵ W) with a standard deviation of 0.11 PW. The geostrophic component of the heat flux has a marked annual cycle following the variability of the Brazil Malvinas Confluence Front, and the geostrophic annual cycle is 180° out of phase with the annual cycle observed in the Ekman fluxes. As a result, the total heat flux shows significant interannual variability with only a small annual cycle. Uncertainties due to different wind products and locations of the sections are independent of the methodology used.     

Detection of systematic errors in XBT data and their correction

Comparison experiment between XBT of T-7 probe and CTD was conducted at 15 stations in the sea area centered on 29°N, 135°E in December 1985. There were systematic errors in XBT temperature profiles in comparison with CTD temperature profiles. The main cause of errors was attributed to an error in the free-fall speed of the XBT probes which was provided by the XBT maker. A previous equation for depth correction proposed by Heinmilleret al. (1983) could not give effective correction for our data. A new equation between the probe depth and the elapsed time from landing of the probe on the water was obtained by the method of adjusting temperature gradients of XBT profiles to those of CTD profiles. This equation agreed with the theoretical result given by Seaver and Kuleshov (1982) much better than that of Heinmilleret al. (1983). Systematic errors due to a scatter of values of the reference resistance and variation of B-constant of thermistors used in XBT also seemed to exist. After an adjustment using the temperature difference between XBT and CTD in the mixed layer with depths of about 100 m, the standard deviation of temperature difference between XBT and CTD from the surface to the depth of 750 m was 0.14°C.     

Meridional heat transport determined with expendable bathythermographs—Part I: Error estimates from model and hydrographic data

Heat transports estimated CTD data collected during the World Ocean Circulation Experiment (WOCE) along the January 1993 30°S hydrographic transect (A10) and the output from a numerical model show a mean heat transport of 0.40 and 0.55±0.24 PW (standard deviation), respectively. The model shows a large annual cycle in heat transport (more than 30% of the variance) with a maximum (minimum) heat transport in July (February) of 0.68 (0.41) PW. Using these data, a method is proposed and evaluated to calculate the heat transport from temperature data obtained from a trans-basin section of expendable bathythermographs (XBTs) profiles. In this method, salinity is estimated from Argo profiles and CTD casts for each XBT temperature observation using statistical relationships between temperature, latitude, longitude and salinity computed along constant-depth surfaces. Full-depth temperature/salinity profiles are obtained by extending the profiles to the bottom of the ocean using deep climatological data. The meridional transport is then determined by using the standard geostrophic method, applying NCEP-derived Ekman transports, and requiring that the salt flux through the Bering Straits be conserved. The results indicate that the methods described here can provide heat transport estimates with a maximum uncertainty of ±0.18 PW (1 PW=10¹⁵ W). Most of this uncertainty is due to the climatology used to estimate the deep structure and issues related to not knowing the absolute velocity field and most especially characterizing barotropic motions. Nevertheless, when the methodology is applied to temperatures collected along 30°S (A10) and direct model integrations, the results are very promising. Results from the numerical model suggest that ageostrophic non-Ekman motions can contribute less than 0.05 PW to heat transport estimates in the South Atlantic.     

Using satellite altimetry to correct mean temperature and salinity fields derived from Argo floats in the ocean regions around Australia

We present results from a suite of methods using in situ temperature and salinity data, and satellite altimetric observations to obtain an enhanced set of mean fields of temperature, salinity (down to 2000-m depth) and steric height (0/2000 m) for a time-specific period (1992–2007). Firstly, the improved global sampling resulting from the introduction of the Argo program, enables a representative determination of the large-scale mean oceanic structure. However, shortcomings in the coverage remain. High variability western boundary current eddy fields, continental slope and shelf boundaries may all be below their optimal sampling requirements. We describe a simple method to supplement and improve standard spatial interpolation schemes and apply them to the available data within the waters surrounding Australia (100°E–180°W; 50°S–10°N). This region includes a major current system, the East Australian Current (EAC), complex topography, unique boundary currents such as the Leeuwin Current, and large ENSO related interannual variability in the southwest Pacific. We use satellite altimetry sea level anomalies (SLA) to directly correct sampling errors in in situ derived mean surface steric height and subsurface temperature and salinity fields. The surface correction is projected through the water column (using an empirical model) to modify the mean subsurface temperature and salinity fields. The errors inherent in all these calculations are examined. The spatial distribution of the barotropic–baroclinic balance is obtained for the region and a ‘baroclinic factor’ to convert the altimetry SLA into an equivalent in situ height is determined. The mean fields in the EAC region are compared with independent estimates on repeated XBT sections, a mooring array and full-depth CTD transects.     

Inter-manufacturer Difference and Temperature Dependency of the Fall-Rate of T-5 Expendable Bathythermograph

The fall-rate of the T-5 expendable bathythermograph (XBT) produced by Tsurumi Seiki (TSK) Co., Ltd and that by Sippican Inc., are intercompared by a series of contemporaneous and colocated measurements with conductivity-temperature-depth (CTD) profilers. It is confirmed that the fall-rates of the two manufacturers' T-5 differ by about 5 percent, despite the fact that they had been believed to be identical for many years. The cause of the difference is discussed on the basis of a detailed cross-examination of the two T-5 models. It is found for the first time that the two models are different in several respects. The manufacturer's fall-rate equation is only applicable to the Sippican T-5, for which Boyd and Linzell's (1993) equation seems to be slightly more accurate. Kizu et al.'s (2005) equation gives a clearly less biased depth than the manufacturers' equation for the TSK T-5. It is also found that the fall-rates of both T-5 models are dependent on water temperature, perhaps because of viscosity. The temperature-dependency of the fall-rate of the TSK T-5 is larger than that of the Sippican T-5.     

Mixing and advection of a cold water cascade over the Wyville Thomson Ridge

The Wyville Thomson Ridge forms part of the barrier to the meridional circulation across which cold Nordic Sea and Arctic water must traverse to reach the Atlantic Ocean. Overflow rates across the ridge are variable (but can be dramatic at times), and may provide a subtle indicator of significant change in the circulation in response to climate change. In spring 2003, a series of CTD sections were conducted during a large overflow event in which Norwegian Sea Deep Water (NSDW) cascaded down the southern side of the ridge into the Rockall Trough at a rate of between 1 and 2 Sv. The NSDW was partially mixed with overlying North Atlantic Water (NAW), and comprised about 1/3rd of the cascading water. The components of NAW and NSDW in the overflow were sufficiently large that there must have been a significant divergence of the inflow through the Faroe-Shetland Channel, and of the outflow through the Faroe Bank Channel.As the plume descended, its temperature near the sea bed warmed by over 3 °C in about a day. Although the slope was quite steep (0.03), the mean speed of the current (typically 0.36 m s⁻¹) was too slow for significant entrainment of NAW to occur (the bulk Richardson number was of order 5). However, very large overturns (up to 50 m) were evident in some CTD profiles, and it is demonstrated from Thorpe scale estimates that the warming of the bottom waters was due to mixing within the plume. It is likely that some of the NSDW had mixed with NAW before it crossed the ridge. The overflow was trapped in a gully, which caused it to descend to great depth (1700 m) at a faster rate, and with less modification due to entrainment, than other overflows in the North Atlantic. The water that flowed into the northern part of the Rockall Trough had a temperature profile that ranged from about 3 to 8 °C. Water with a temperature of >6 °C probably escaped into the Iceland Basin, between the banks that line the north-western part of the Trough. Colder water (< 6 °C) must have travelled down the eastern side of the Rockall Bank, and may have had a volume flux of up to 1.5 Sv.     

The western boundary current in the Bransfield Strait,Antarctica

The Bransfield Strait west of the Antarctic Peninsula has been considered as a highly productive region for all trophic levels from primary production, to zooplankton aggregations, especially krill, to birds and mammals. The western boundary current, referred to as the Bransfield Current, plays an important role in determining the transport and retention of biota in the Bransfield Strait. Following the study of surface current characteristics in the strait using 39 tracks of mixed-layer drifters deployed between 1988 and 1990, a high-resolution transect of temperature, salinity and current measurements crossing the Bransfield Current was conducted between 13 and 14 March 2004, for understanding its horizontal and vertical structure and dynamics. The results from current, temperature and salinity measurements using a vessel mounted narrow band acoustic doppler current profiler and conductivity–temperature–depth (CTD) sensors revealed the magnitude of this current of approximately 50 cm/s within a horizontal distance of 15 km associated with a narrow and deep density front 4–6 km wide and 500 m deep. The comparison between the direct current measurements and the geostrophic current estimates from the density field implies that the Bransfield Current is geostrophically balanced. The mechanism forming this current is explored with Sverdrup dynamics. Results indicate that the negative wind stress curl and β-effect lead to a southwestward transport in the Bransfield Strait. When this transport is restricted by land and shelves, a narrow western boundary current is formed.     

Cross-slope zonation of erosion and deposition in the Faeroe-Shetland Channel,North Atlantic Ocean

Boundary currents and internal waves determine cross-slope zonation of erosion and deposition in the Faeroe-Shetland Channel. Currents were measured at 8 and 34–50 m above the bottom at three mooring sites (502, 595 and 708 m depth) for 14 days. The structure of the water column was evaluated from CTD sections, and included nepheloid layers and particulate matter concentrations. Indicators for recent deposition in the sediment (organic carbon, phytopigments, ²¹⁰Pb) were measured at eight stations across the slope. Strong near-bottom currents at the upper slope sustain down-slope particle transport in a benthic nepheloid layer, which is eroded under the influence of critically reflecting M₂ internal tidal waves at 350–550 m, where the major pycnocline meets the sloping bottom. Beam attenuation profiles confirmed the presence of intermediate nepheloid layers intruding into the Channel along the major pycnocline, and elevated concentrations of particulate matter and chlorophyll-a were measured at this depth. Near-bottom currents decreased with depth, thus allowing particle deposition down the slope. Inventories of excess ²¹⁰Pb activity in the sediment deeper than 600 m were higher than what was expected on the basis of atmospheric input of ²¹⁰Pb and production in the water column, thus indicating additional lateral inputs. Simple calculations showed that off-slope input of particles from areas shallower than 600 m may be responsible for the enhanced deposition at greater depths.     

Corrections to ocean surface forcing in the California Current System using 4D variational data assimilation

The option for surface forcing correction, recently developed in the 4D-variational (4DVAR) data assimilation systems of the Regional Ocean Model System (ROMS), is presented. Assimilation of remotely-sensed (satellite sea surface height anomaly and sea surface temperature) and in situ (from mechanical and expendable bathythermographs, Argo floats and CTD profiles) oceanic observations has been applied in a realistic, high resolution configuration of the California Current System (CCS) to sequentially correct model initial conditions and surface forcing, using the Incremental Strong constraint version of ROMS-4DVAR (ROMS-IS4DVAR). Results from both twin and real data experiments are presented where it is demonstrated that ROMS-IS4DVAR always reduces the difference between the model and the observations that are assimilated. However, without corrections to the surface forcing, the assimilation of surface data can degrade the temperature structure at depth. When using surface forcing adjustment in ROMS-IS4DVAR the system does not degrade the temperature structure at depth, because differences between the model and surface observations can be reduced through corrections to surface forcing rather than to temperature at depth. However, corrections to surface forcing can generate abnormal spatial and temporal variability in the structure of the wind stress or surface heat flux fields if not properly constrained. This behavior can be partially controlled via the choice of decorrelation length scales that are assumed for the forcing errors. Abnormal forcing corrections may also arise due to the effects of model error which are not accounted for in IS4DVAR. In particular, data assimilation tends to weaken the alongshore wind stress in an attempt to reduce the rate of coastal upwelling, which seems to be too strong due to other sources of error. However, corrections to wind stress and surface heat flux improve systematically the ocean state analyses. Trends in the correction of surface heat fluxes indicate that, given the ocean model used and its potential limitations, the heat flux data from the Coupled Ocean–Atmosphere Mesoscale Prediction System (COAMPS) used to impose surface conditions in the model are generally too low except in spring-summer, in the upwelling region, where they are too high. Comparisons with independent data provide confidence in the resulting forecast ocean circulation on timescales ～14 days, with less than 1.5 °C, 0.3 psu, and 9 cm RMS error in temperature, salinity and sea surface height anomaly, respectively, compared to observations.     

Estimating the compensation irradiance in the ocean: The importance of accounting for non-photosynthetic uptake of inorganic carbon

The compensation irradiance, the irradiance at which net photosynthesis is zero over a 24-h period, was estimated at station ALOHA (22°45′N, 158°W) from analysis of ¹⁴C uptake rates measured from 8 January 1989 to 13 June 1990 at depths ranging from 5 to 175 m. The estimates were made on the basis of linear regressions of the difference between light bottle and dark bottle ¹⁴C uptake in the light-limited region of the euphotic zone and determination of the depth at which the difference between the uptake rates was zero. About half of the non-photosynthetic ¹⁴C uptake at the compensation irradiance could be attributed to chemolithoautotrophy; the remainder was presumably due to anaplerotic processes. Deriving the compensation irradiance by extrapolating dawn-to-dawn light-bottle uptake above the compensation irradiance to zero resulted in underestimation of the compensation irradiance by a factor of 2. We estimated the compensation irradiance at station ALOHA to be 0.054 mol-photons m⁻² d⁻¹, about 0.11% of surface 400–700 nm radiation and 1% of surface 475-nm (blue) light.     

Seasonal and inter-annual biogeochemical variations in the Porcupine Abyssal Plain 2003–2005 associated with winter mixing and surface circulation

We present a 3-year multidisciplinary biogeochemical data set taken in situ at the Porcupine Abyssal Plain (PAP) time-series observatory in the Northeast Atlantic (49°N, 16.5°W; water depth ∼4850 m) for the period 2003 to 2005. The high-resolution year-round autonomous measurements include temperature, salinity, chlorophyll-a (derived from in situ chlorophyll-fluorescence) and inorganic nitrate, all at a nominal depth of 30 m on an Eulerian observatory mooring. This study compares these in situ time-series data with satellite chlorophyll-a data, regional data from a ship of opportunity, mixed-layer depth measurements from profiling Argo floats and lateral advection estimates from altimetry. This combined and substantial data set is used to analyse seasonal and inter-annual variability in hydrography and nitrate concentrations in relation to convective mixing and lateral advection. The PAP observatory site is in the inter-gyre region of the North Atlantic where convective mixing ranges from 25 m in the summer to over 400 m in winter when nutrients are supplied to the surface. Small inter-annual changes in the winter mixed layer can result in large changes in nitrate supply and productivity. However the decrease in maximum winter nitrate over the three-year period, from 10 to 7 mmol m⁻³, cannot be fully explained by convective mixing. Trajectories leading to the PAP site, computed from altimetry-derived geostrophic velocities, confirm that lateral advection cannot be ignored at this site and may be an important process along with convective mixing. Over the three years, there is an associated decrease in new production calculated from nitrate assimilation from 85.4 to 40.3±4.3 gCm⁻² a⁻¹. This confirms year-to-year variability in primary production seen in model estimates for the region. The continuous in situ dataset also shows inter-annual variation in the timing of the spring bloom due to variations in heat flux; the 2005 bloom occurred earlier than in 2004.     

Total kinetic energy in four global eddying ocean circulation models and over 5000 current meter records

We compare the total kinetic energy (TKE) in four global eddying ocean circulation simulations with a global dataset of over 5000, quality controlled, moored current meter records. At individual mooring sites, there was considerable scatter between models and observations that was greater than estimated statistical uncertainty. Averaging over all current meter records in various depth ranges, all four models had mean TKE within a factor of two of observations above 3500 m, and within a factor of three below 3500 m. With the exception of observations between 20 and 100 m, the models tended to straddle the observations. However, individual models had clear biases. The free running (no data assimilation) model biases were largest below 2000 m. Idealized simulations revealed that the parameterized bottom boundary layer tidal currents were not likely the source of the problem, but that reducing quadratic bottom drag coefficient may improve the fit with deep observations. Data assimilation clearly improved the model-observation comparison, especially below 2000 m, despite assimilated data existing mostly above this depth and only south of 47 °N. Different diagnostics revealed different aspects of the comparison, though in general the models appeared to be in an eddying-regime with TKE that compared reasonably well with observations.     

Long-term hydrographic changes at 52 and 66°W in the North Atlantic Subtropical Gyre & Caribbean

In July–September 1997 two hydrographic lines were done in the western N. Atlantic along longitudes of 52 and 66°W as part of the WOCE one-time hydrographic survey of the oceans. Each of these two lines approximately repeated earlier ones done during the International Geophysical Year(s) (IGY) and the mid-1980s. Because of this repeated sampling, long-term hydrographic changes in the water masses can be examined. In this report, we focus on temperature and salinity changes within the subtropical gyre mainly between latitudes of 20 and 35°N and compare our results to those presented by Bryden et al. (1996), who examined changes along a zonal line at 24°N, most recently occupied in 1992. Since this most recent 24°N section in 1992, substantial changes have occurred in the western part of the subtropical gyre at the depths of the Labrador Sea Water (LSW). In particular, we see clear evidence for colder, fresher Labrador Sea Water throughout the gyre on our two recent sections that was not yet present in 1992 at similar longitudes along 24°N. At shallower depths inhabited by waters that are an admixture of Mediterranean (MW) and Antarctic Intermediate Waters (AAIW), our recent survey shows an increase in salinity, which can only be attributed to changes in water masses on potential temperature or neutral density surfaces. Furthermore, waters above the MW/AAIW layer and into the deeper part of the main pycnocline have continued to become saltier and warmer throughout the 40-year period spanned by our sections. These latter changes have been dominantly due to a vertical sinking of density surfaces as T/S changes in density surfaces are small, but depths of individual T/S horizons have increased with time. The net change since the IGY shows a mean temperature increase between 800 and 2500 m depth at a rate of 0.57°C/century with a corresponding steric sea level rise of 1 mm/yr, and a net downward heave with small values near the top and bottom, and a maximum rate of −2.7 m/yr at 1800 m depth. Changes in the deep Caribbean indicate a warming since the IGY due to temperature increases of the inflowing source waters in the subtropical gyre at 1800m depth, but no significant change in the deep salinity.     

Deep Caribbean Sea warming

Data collected from hydrographic stations occupied within the Venezuelan and Columbian basins of the Caribbean Sea from 1922 through 2003 are analyzed to study the decadal variability of deep temperature in the region. The analysis focuses on waters below the 1815-m sill depth of the Anegada–Jungfern Passage. Relatively dense waters (compared to those in the deep Caribbean) from the North Atlantic spill over this sill to ventilate the deep Caribbean Sea. Deep warming at a rate of over 0.01 °C decade^–1 below this sill depth appears to have commenced in the 1970s after a period of relatively constant deep Caribbean Sea temperatures extending at least as far back as the 1920s. Conductivity–temperature–depth station data from World Ocean Circulation Experiment Section A22 along 66°W taken in 1997 and again in 2003 provide an especially precise, albeit geographically limited, estimate of this warming over that 6-year period. They also suggest a small (0.001 PSS-78, about the size of expected measurement biases) deep freshening. The warming is about 10 times larger than the size of geothermal heating in the region, and is of the same magnitude as the average global upper-ocean heat uptake over a recent 50-year period. Together with the freshening, the warming contributes about 0.012 m decade^–1 of sea level rise in portions of the Caribbean Sea with bottom depths around 5000 m.     

Hydrodynamic controls on cold-water coral growth and carbonate-mound development at the SW and SE Rockall Trough Margin,NE Atlantic Ocean

Long-term (⩽1-year) records obtained by seabed observatories (BOBO) and repeated (24-h) CTD casts show the presence of a highly energetic environment in and around two cold-water carbonate-mound provinces, on the Southwest and Southeast Rockall Trough (SW and SE RT) margin. Carbonate mounds, covered with a thriving coral cover, are embedded mainly in the Eastern North Atlantic Water (ENAW) and are observed in a confined bathymetric zone between 600 and 1000 m water depth. Cold-water corals seem to be restricted in their growth by temperature and food availability. The presence of living corals on top of the carbonate mounds appears linked to the presence of internal waves and tidal currents in the water column, and consequently carbonate mound structures are shaped by the local hydrodynamic regime. Mound clusters have an elongated shape perpendicular to the regional contours and corresponding to the direction of the highest current speeds. On the SW RT margin temperature, salinity and current speed reflect a diurnal tidal pattern, causing maximum temperature variations at 900 m depth of more than 3 °C. Current speeds up to 45 cm s⁻¹ occur, and a residual current of 10 cm s⁻¹ is directed along the slope to the southwest. At the SE RT margin the temperature of the bottom water fluctuates more than 1 °C with a semi-diurnal tidal cyclicity. Amplitudes of average and peak current speeds here are comparable with those measured on the southwest margin, but the residual current in this area is directed to the northeast. Tidal currents and internal waves at both margins force the formation of intermediate and bottom nepheloid layers and bring fresh food particles with increased velocity to the mounds. The distribution of corals in both mound areas is considered directly related to the presence of enhanced turbidity. An increase in temperature can be directly related to an increase in the amount of particles in the water column. Current velocity increases when a transition occurs from cold to warm waters. High current velocities prevent local sedimentation but provide sufficient food particles to the corals, so that the corals thrive at the mound summits.     

A New Fall-Rate Equation for T-5 Expendable Bathythermograph (XBT) by TSK

The accuracy of the manufacturer’s fall-rate equation for the T-5 Model of expendable bathythermograph (XBT) has been investigated based on about 300 collocated pairs of XBT-CTD (Conductivity-Temperature-Depth profiler) measurements in various climatological regions. We found that the equation systematically overestimates depth by about 5% for the T-5 produced by Tsurumi Seiki, Co. Ltd. (TSK), but almost no bias is associated with the T-5 produced by Sippican, Inc., in USA. The cause of this difference is not clear, because the two manufacturers’ T-5 probes are reported to have identical shape and weight in water. We propose a new fall-rate equation for the TSK T-5: z(t) = 6.54071t - 0.0018691t ², where z(t) is depth in meters at time, t, in seconds.     

Fractionation during remineralization of organic matter in the ocean

Remineralization ratios (–O₂:P, C_org.:P, N:P) in the ocean are estimated from ocean tracer data using a new approach, which takes into account the effects of local exchange across neutral surfaces. This approach is applied to temperature, salinity, phosphate, nitrate, dissolved oxygen, alkalinity, and dissolved inorganic carbon data from the low- and mid-latitude Pacific, Indian, and South Atlantic Oceans. The consideration of local exchange effects tends to reduce the –O₂:P and C_org.:P remineralization estimates above 1500 m compared to earlier estimates. Below 1500 m, exchange effects can be neglected (except in the South Atlantic) and earlier estimates appear robust. In the deep South Atlantic, the consideration of these effects leads to increased –O₂:P and C_org.:P remineralization ratio estimates, bringing them more in line with the robust deep ocean estimates. For reasonable, open ocean mixing coefficient values and several choices for phosphate remineralization rate profiles, –O₂:P (C_org.:P) remineralization ratios in the ocean increase from about 140 (100) at 750 m depth to about 170 (130) at 1500 m and remain so deeper down. Such an increase down through the upper ocean thermocline implies significant fractionation during remineralization of organic matter—nutrients are released higher in the water column than inorganic carbon. These results also argue for a –O₂:P (C_org.:P) uptake ratio in new production of about 140–150 (100–110). N:P remineralization ratios decrease from about 15 at 750 m to about 12 at 1500–2000 m. This may reflect a “true” N:P remineralization (and uptake) ratio of about 16, modified by denitrification.These results imply that applications of derived, quasi-conservative tracers, based on the assumption of constant remineralization ratios, may be subject to significant error for depths less than 1500 m. In addition, present Ocean General Circulation Models of the natural carbon cycle in the ocean–atmosphere system assume remineralization to occur without fractionation but have problems simulating observed, pre-industrial levels of atmospheric pCO₂, given observed ocean inventories of alkalinity and dissolved inorganic carbon. Implementation of uptake and (depth-dependent) remineralization ratios estimated here would likely reduce this problem considerably. Furthermore, calculations with a simple global carbon cycle model show that fractionation in the modern ocean, as estimated in the present work, has reduced atmospheric pCO₂ by more than 20 ppm below the level it would have had without fractionation.     

Antipodal acoustic thermometry: 1960, 2004

On 21 March 1960, sounds from three 300-lb depth charges deployed at 5.5-min intervals off Perth, Australia were recorded by the SOFAR station at Bermuda. The recorded travel time of these signals, about 13,375 s, is a historical measure of the ocean temperature averaged across several ocean basins. The 1960 travel time measurement has about 3-s precision. High-resolution global ocean state estimates for 2004 from the “Estimating the Circulation and Climate of the Ocean, Phase II” (ECCO2) project were combined with ray tracing to determine the paths followed by the acoustic signals. The acoustic paths are refracted geodesics that are slightly deflected by either small-scale topographic features in the Southern Ocean or the coast of Brazil. The refractive influences of intense, small-scale oceanographic features, such as Agulhas Rings or eddies in the Antarctic Circumpolar Current, greatly reduce the necessary topographic deflection and cause the acoustic paths to meander in time. The ECCO2 ocean state estimates, which are constrained by model dynamics and available data, were used to compute present-day travel times. Measured and computed arrival coda were in good agreement. Based on recent estimates of warming of the upper ocean, the travel-time change over the past half-century was nominally expected to be about −9 s, but little difference between measured (1960) and computed (2004) travel times was found. Taking into account uncertainties in the 1960 measurements, the 2004 ocean state estimates, and other approximations, the ocean temperature averaged along the sound channel axis over the antipodal paths has warmed at a rate less than about 4.6 m °C yr⁻¹ (95% confidence). At this time, the estimated uncertainties are comparable in size to the expected warming signal, however.     

Behaviors of dissolved and particulate Co,Ni, Cu,Zn, Cd and Pb during a mesoscale Fe-enrichment experiment (SEEDS II) in the western North Pacific

During mesoscale Fe enrichment (SEEDS II) in the western North Pacific ocean, we investigated dissolved and particulate Co, Ni, Cu, Zn, Cd and Pb in seawater from both field observation and shipboard bottle incubation of a natural phytoplankton assemblage with Fe addition. Before the Fe enrichment, strong correlations between dissolved trace metals (Ni, Zn and Cd) and PO₄³⁻, and between particulate trace metals (Ni, Zn and Cd) and chlorophyll-a were obtained, suggesting that biogeochemical cycles mainly control the distributions of Ni, Zn and Cd in the study area. Average concentrations of dissolved Co, Ni, Cu, Zn, Cd and Pb in the surface mixed layer (0–20 m) were 70 pM, 4.9, 2.1, 1.6, 0.48 nM and 52 pM, respectively, and those for the particulate species were 1.7 pM, 0.052, 0.094, 0.46, 0.037 nM and 5.2 pM, respectively. After Fe enrichment, chlorophyll-a increased 3 fold (up to 3 μg L⁻¹) during developing phases of the bloom (<12 days). Mesozooplankton biomass also increased. Particulate Co, Ni, Cu and Cd inside the patch hinted at an increase in the concentrations, but there were no analytically significant differences between concentrations inside and outside the patch. The bottle incubation with Fe addition (1 nM) showed an increase in chlorophyll-a (8.9 μg L⁻¹) and raised the particulate fraction up to 3–45% for all the metals, accompanying changes in Si/P, Zn/P and Cd/P. These results suggest that Fe addition lead to changes in biogeochemical cycling of trace metals. The comparison between the mesoscale Fe enrichment and the bottle incubation experiment suggests that although Fe was a limiting factor for the growth of phytoplankton, the enhanced biomass of mesozooplankton also limited the growth of phytoplankton and the transformation of trace metal speciation during the mesoscale Fe enrichment. Sediment trap data and the elemental ratios taken up by phytoplankton suggest that export loss was another reason that no detectable change in the concentrations of particulate trace metals was observed during the mesoscale Fe enrichment.