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.     

A global heat and freshwater forcing dataset for ocean models

A global dataset based on the ECMWF Re-Analyses (ERA) is presented that can be used as surface boundary conditions for ocean models with sea-ice components. The definition of these conditions is based on bulk formulae. To study the mean ocean circulation, a mean annual cycle on a daily basis was constructed from ERA for all relevant parameters including wind stress. Continental runoff is considered by using information about the catchment areas of the rivers and about the main drainage basins. The bulk formulae were extended by using sea ice concentration.To estimate meridional heat transports (MHT) and to avoid any drift in ocean model simulations, the heat and fresh water budgets have been closed by applying an inverse procedure to fine-tune the fluxes towards observed transports. To improve the MHTs on the Southern Hemisphere the winds and the short wave radiation at southern higher latitudes should be corrected. Furthermore, tests were performed concerning short wave radiation which was increased in the tropics and decreased in the subsidence zones.The heat and fresh water fluxes are assessed by using a scheme of Macdonald and Wunsch based on hydrographic sections. The net heat fluxes of ERA and of the forcing dataset are consistent with the heat flux divergences and convergences estimated by this scheme except for parts of the South Atlantic and the Indian Ocean sector of the Southern Ocean where none of these datasets is consistent with these estimates. In the subtropical South Indian Ocean the forcing dataset is consistent with these estimates while ERA are not. The flux components of ERA and the forcing dataset were compared to several observational datasets (SRB, SOC, HOAPS, GPCP, and CMAP). For each component, at least one of these datasets (especially HOAPS) supports the effects of the inverse procedure and the bulk formulae almost globally with some regional exceptions: short wave radiation in the tropical oceans and the subtropical North Atlantic, latent heat flux at higher latitudes, and precipitation in the northern North Atlantic.Comparisons to the NCAR/NCEP Re-Analyses (NRA) (versions 1 and 2) and the ECHAM model in place of ERA lead to similar results. In the North Atlantic the net heat fluxes of the model based datasets approach the hydrographic estimate with increasing resolution. Applied to any ocean/sea-ice model and compared to ERA, the forcing dataset would induce only a relative small net sea-surface buoyancy loss.A comparison of the forcing dataset to measurements made using one buoy deployed in the western Pacific warm pool and five buoys deployed in the subduction region of the Northeast Atlantic shows that at the site of the first buoy the net heat fluxes of the forcing dataset are in poorer agreement than those of ERA. At the sites of two subduction buoys both datasets show the same level of agreement within the error bars specified. At the sites of the three remaining subduction buoys the forcing dataset shows a marginal improvement on ERA.     

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.     

Sensitivity of near-surface Tropical Instability Waves to submonthly wind forcing in the tropical Atlantic

This study demonstrates the sensitivity of the near-surface properties in the tropical Atlantic Ocean to the high-frequency of the winds in numerical simulations. At intra-seasonal timescales (2–50 days), two distinct period ranges dominate the variability in the upper ocean: periods between 2 and 20 days, which are essentially wind-forced and periods between 20 and 50 days, due mostly to Tropical Instability Waves (TIWs). Using a numerical model forced by different wind fields, it is shown that the characteristics of the intra-seasonal variability in the ocean surface mixed-layer are strongly dependent on the wind forcing. Submonthly winds are shown to force large variability in the upper ocean that can strikingly decrease the amplitude of the TIWs in the mixed-layer and their imprint on the horizontal distribution of sea surface temperatures. Wind products containing too much energy at submonthly periods thus prevent wind-forced simulations from reproducing a realistic surface signature of TIWs, when compared to satellite observations of sea surface temperature. In addition, submonthly wind variability may be responsible for part of the observed interannual variability of the TIW signature in the temperature. The impact of submonthly winds is strongest in the mixed-layer: beneath the mixed-layer, all simulations show similar characteristics of the TIWs.     

High resolution profiles of vertical particulate organic matter export off Cape Blanc,Mauritania: Degradation processes and ballasting effects

Vertical carbon fluxes between the surface and 2500 m depth were estimated from in situ profiles of particle size distributions and abundances me/asured off Cape Blanc (Mauritania) related to deep ocean sediment traps. Vertical mass fluxes off Cape Blanc were significantly higher than recent global estimates in the open ocean. The aggregates off Cape Blanc contained high amounts of ballast material due to the presence of coccoliths and fine-grained dust from the Sahara desert, leading to a dominance of small and fast-settling aggregates. The largest changes in vertical fluxes were observed in the surface waters (<250 m), and, thus, showing this site to be the most important zone for aggregate formation and degradation. The degradation length scale (L), i.e. the fractional degradation of aggregates per meter settled, was estimated from vertical fluxes derived from the particle size distribution through the water column. This was compared with fractional remineralization rate of aggregates per meter settled derived from direct ship-board measurements of sinking velocity and small-scale O₂ fluxes to aggregates measured by micro-sensors. Microbial respiration by attached bacteria alone could not explain the degradation of organic matter in the upper ocean. Instead, flux feeding from zooplankton organisms was indicated as the dominant degradation process of aggregated carbon in the surface ocean. Below the surface ocean, microbes became more important for the degradation as zooplankton was rare at these depths.     

Heat,volume and chemical fluxes from submarine venting: A synthesis of results from the Rainbow hydrothermal field, 36°N MAR

High-temperature hydrothermal activity occurs in all ocean basins and along ridge crests of all spreading rates. While it has long been recognized that the fluxes associated with such venting are large, precise quantification of their impact on ocean biogeochemistry has proved elusive. Here, we report a comprehensive study of heat, fluid and chemical fluxes from a single submarine hydrothermal field. To achieve this, we have exploited the integrating nature of the non-buoyant plume dispersing above the Rainbow hydrothermal field, a long-lived and tectonically hosted high-temperature vent site on the Mid-Atlantic Ridge. Our calculations yield heat and volume fluxes for high-temperature fluids exiting the seafloor of ～0.5 GW and 450 L s^?1, together with accompanying chemical fluxes, for Fe, Mn and CH₄ of ～10, ～1 and ～1 mol s^?1, respectively. Accompanying fluxes for 25 additional chemical species that are associated with Fe-rich plume particles have also been calculated as they are transported away from the Rainbow vent site before settling to the seabed. High-temperature venting has been found to recur at least once every ～100 km along all slow-spreading ridges investigated to-date, with half of all known sites on the Mid-Atlantic Ridge occurring as long-lived and tectonically hosted systems. If these patterns persist along all slow- and ultraslow-spreading ridges, high-temperature venting of the kind reported here could account for ～50% of the on-axis hydrothermal heat flux along ～30,000 km of the ～55,000 km global ridge crest.     

A note on boundary conditions for salt and freshwater balances

The derivation of surface boundary conditions for salt and freshwater budgets, in the presence of an imbalance of evaporation and precipitation, is revisited in order to: (i) clarify the physical concepts involved and (ii) derive consistent approximations. Proper surface boundary conditions can be obtained by assuming that the air–sea interface is a material surface for the salt continuum. This requires all the mass exchange at the air–sea interface to be fulfilled by the freshwater continuum. The need of diffusive fluxes is crucial for this assumption to hold. If the air–sea interface varies in time (either by seawater advection, an imbalance of evaporation minus precipitation, or a combination of both) there must be a corresponding salt flux, in order for the transport across such moving surface to be nil. Here it is essential to distinguish the flux vector of a property from its transport across a (possibly moving) surface. Only in a motionless ocean surface approximation, which includes the classical rigid lid approximation, the advective and diffusive salt fluxes perpendicular to the surface cancel out exactly.     

Mass, Heat and Salt Balances in the Eastern Barents Sea Obtained by Inversion of Hydrographic Section Data

Standard hydrological section data, collected in the eastern Barents Sea in September 1997, have been analyzed using a variational data assimilation technique. This method allows us to obtain temperature, salinity and velocity fields that are consistent with observations and dynamically balanced within the framework of a steady-state model describing large-scale nearly geostrophic circulation. Error bars of the optimized fields are computed by explicit inversion of the Hessian matrix. The optimized velocity field is in agreement with independent velocity observations derived from surface drifter trajectories in the southwestern part of the Barents Sea. Optimized fields provide the following estimates of integral characteristics of the circulation in the region: i) the North Cape current transport is 2.12 ± 0.25 Sv; ii) the Karskie Vorota Strait throughflow is 0.7 ± 0.06 Sv; iii) heat flux with Atlantic water is 4.7 ± 0.16⋅10¹¹ W; iv) salt import from the Atlantic Ocean is 7.41 ± 0.46⋅10³ kg/s. The imbalance of the heat budget in the eastern part of the Barents Sea indicates the presence of statistically insignificant surface heat fluxes which are less than 1 W/m². This revised version was published online in July 2006 with corrections to the Cover Date.     

What processes drive the ocean heat transport?

The ocean contributes to regulating the Earth’s climate through its ability to transport heat from the equator to the poles. In this study we use long simulations of an ocean model to investigate whether the heat transport is carried primarily by wind-driven gyres or whether it is dominated by deep circulations associated with abyssal mixing and high latitude convection. The heat transport is computed as a function of temperature classes. In the Pacific and Indian ocean, the bulk of the heat transport is associated with wind-driven gyres confined to the thermocline. In the Atlantic, the thermocline gyres account for only 40% of the total heat transport. The remaining 60% is associated with a circulation reaching down to cold waters below the thermocline. Using a series of sensitivity experiments, we show that this deep heat transport is primarily set by the strength and patterns of surface winds and only secondarily by diabatic processes at high latitudes in the North Atlantic. Abyssal mixing below 2000 m has hardly any impact on ocean heat transport. A major implication is that the role of the ocean in regulating Earth’s climate strongly depends on how surface winds change across different climates in both hemispheres at low and high latitudes.     

Application of underwater optical data to estimation of primary productivity

We conducted time-series observations of optical fields near the base of the euphotic zone (approximately 40 m) using moored automatic optical sensors at a time-series station in the Western Pacific Subarctic Gyre from March 2005 to July 2006 (with some gaps). We used the ratio of photosynthetically available radiation at the surface (surface PAR) to in situ quantum irradiance (in situ QI) at about 40 m as an index of opacity (surface PAR/in situ QI), which began to increase in the middle of April and peaked between the end of June and the middle of July 2005. This ratio then decreased toward winter. The ratio increased again beginning in January 2006, and large peaks were observed in June and July 2006. As an index of chlorophyll abundance we used the ratio of spectral irradiance at wavelengths of 555 and 443 nm (Ed₅₅₅/Ed₄₄₃) at about 40 m; seasonal variability of this ratio synchronized well with the attenuation coefficient “k” estimated with surface PAR, in situ QI, and BLOOMS depth. We estimated primary productivity (PP) using Ed₅₅₅/Ed₄₄₃ and an empirical equation based on a previous model but improved on the basis of shipboard observations. Estimated PP agreed well with observed PP. Seasonal variability of estimated PP was synchronized with that of organic carbon flux observed by sediment traps from approximately 150, 540, 1000, and 5000 m. This study demonstrates that time-series observations of in situ optical fields could contribute to the estimation of primary productivity and the study of the biological pump in the ocean.     

Ice–ocean interactions during shallow convection under conditions of steady winds: three-dimensional numerical studies

A three-dimensional hydrodynamic ocean model coupled to a thermohydrodynamic model for young sea ice is applied to study shallow haline convection in the central Greenland Sea, with an emphasis on sub-mesoscale ice–ocean interactions. Two types of young sea ice are distinguished; i.e., frazil and pancake ice, both acting different on surface heat, salt, and momentum fluxes. Two scenarios are considered: (a) continued frazil-ice production during steady winds, and (b) the same scenario but with the intermittent formation of pancake ice during a short intervening period of low winds. Brine release due to new-ice production creates shallow convection in both cases. Under conditions of continued frazil-ice production, ice streaks develop at the sea surface, finally becoming oriented roughly parallel to the wind. These streaks are the result of convective plumes that induce organized patterns of convergent and divergent surface currents. Frazil-ice is herded into convergence zones where it becomes as thick as 6 m within 24 h. The studies suggest a strong relationship between streak spacing and the penetration depth of convection, given by an aspect ratio in the range of 2–3. After pancake ice has been formed, however, the organized ice streaks vanish, developing into complex patterns of pancake ice. This finding is in agreement with recent field observations in the Greenland Sea Is-Odden ice tongue. With the existence of pancake ice, moreover, the surface-averaged buoyancy flux decreases and is determined from the integral of local sub-mesoscale ice–ocean interactions.     

Recent changes in the North Atlantic circulation simulated with eddy-permitting and eddy-resolving ocean models

Three eddy-permitting (1/4°) versions and one eddy-resolving (1/12°) version of the OCCAM ocean model are used to simulate the World Ocean circulation since 1985. The first eddy-permitting simulation has been used extensively in previous studies, and provides a point of reference. A second, improved, eddy-permitting simulation is forced in the same manner as the eddy-resolving simulation, with a dataset based on a blend of NCEP re-analysis and satellite data. The third eddy-permitting simulation is forced with a different dataset, based on the ERA-40 re-analysis data. Inter-comparison of these simulations in the North Atlantic clarifies the relative importance of resolution and choice of forcing dataset, for simulating the mean state and recent variability of the basin-scale circulation in that region. Differences between the first and second eddy-permitting simulations additionally reveal an erroneous influence of sea ice on surface salinity, dense water formation, and the meridional overturning circulation. Simulations are further evaluated in terms of long-term mean ocean heat transport at selected latitudes (for which hydrographic estimates are available) and sea surface temperature errors (relative to observations). By these criteria, closest agreement with observations is obtained for the eddy-resolving simulation. In this simulation, there is also a weak decadal variation in mid-latitudes, with heat transport strongest, by around 0.2 PW, in the mid-1990s. In two of the eddy-permitting simulations, by contrast, heat transport weakens through the study period, by up to 0.4 PW in mid-latitudes. The most notable changes of heat transport in all simulations are linked to a weakening of the subpolar gyre, rather than changes in the meridional overturning circulation. It is concluded that recent changes in the structure of mid-latitude heat transport in the North Atlantic are more accurately represented if eddies are explicitly resolved.     

Radiolarian fluxes from the southern Bay of Bengal: sediment trap results

A study of radiolarian fluxes collected during 1991–93 from time-series sediment traps deployed at 1071 and 3010 m water depth in the southern Bay of Bengal (SBBT) yielded 40 species/groups of radiolarians. Among the order Polycystina, the species of sub-order Spumellaria were by far the most abundant (∼95%) followed by sub-order Nassellaria (5%). This is contrary to reports from the Atlantic and Pacific Oceans and is attributed to the prevailing hyposaline condition resulting from the monsoonal rainfall. Higher radiolarian fluxes occurred during March–May, when moderate salinity and a high sea surface temperature (SST) regime prevailed at the trap site. R-mode cluster analysis of the radiolarian flux data revealed three assemblages represented by the cooler (A) and warmer (C) surface dwelling fauna (0–50 m) dominated by spumellarians, and a deeper dwelling (B) sub-surface fauna (50–100 m) associated with deep dwelling (>100 m) nassellarian species. Spongaster tetras tetras, a surface water radiolarian species, exhibited its preference for high SST and moderate salinity conditions during the pre-monsoon season (March–May). Radiolarian fluxes responded to seasonal changes in SST and salinity variations due to the monsoonal precipitation, and the freshwater runoff from the Indian rivers causing a hyposaline condition in the Bay of Bengal. Results imply that the radiolarian assemblages in the down core data may reveal the monsoonal history in the geological past.     

Analysis of three years of boundary layer observations over the Gulf of Mexico and its shores

Boundary layer observations were made over the Gulf of Mexico over a 3-year period in order to develop and test methods for estimating surface fluxes and boundary layer wind fields. In addition to routinely available buoy and CMAN surface data, six 915 MHz radar wind profilers (RWPs) and RASS profilers were mounted on oil platforms and on the shore. Estimates of surface momentum, sensible heat, and latent heat fluxes have been made from the surface observations using the COARE software. Simulations by the National Weather Service's Eta meteorological model are compared with the observations of surface fluxes and wind profiles. The boundary layer is found to be unstable over 90% of the time, and latent heat fluxes are about five to ten times larger than sensible heat fluxes, as usually found over tropical oceans. Eta model simulations of surface fluxes are within about ±50% of COARE estimates of the fluxes based on surface observations. Most of the time, COARE-derived fluxes at 11 sites are within a factor of two of each other at any given hour. In multi-day case studies, COARE calculations are found to agree with Eta model simulations of these fluxes and parameters within a factor of two most of the time. Eta model simulations of wind speeds in the boundary layer tend to exceed the RWP observations by 1–2 m s⁻¹ near shore and by 2–6 m s⁻¹ at distances of 100–200 km offshore.     

Summer and winter air–sea CO2 fluxes in the Southern Ocean

The seasonal variability of the carbon dioxide (CO₂) system in the Southern Ocean, south of 50°S, is analysed from observations obtained in January and August 2000 during OISO cruises conducted in the Indian Antarctic sector. In the seasonal ice zone, SIZ (south of 58°S), surface ocean CO₂ concentrations are well below equilibrium during austral summer. During this season, when sea-ice is not obstructing gas exchange at the air–sea interface, the oceanic CO₂ sink ranges from −2 to −4 mmol/m²/d in the SIZ. In the permanent open ocean zone, POOZ (50–58°S), surface oceanic fugacity fCO₂ increases from summer to winter. The seasonal fCO₂ variations (from 10 to 30 μatm) are relatively low compared to seasonal amplitudes observed in the subtropics or the subantarctic zones. However, these variations in the POOZ are large enough to cross the atmospheric level from summer to winter. Therefore, this region is neither a permanent CO₂ sink nor a permanent CO₂ source. In the POOZ, air–sea CO₂ fluxes calculated from observations are about −1.1 mmol/m²/d in January (a small sink) and 2.5 mmol/m²/d in August (a source). These estimates obtained for only two periods of the year need to be extrapolated on a monthly scale in order to calculate an integrated air–sea CO₂ flux on an annual basis. For doing this, we use a biogeochemical model that creates annual cycles for nitrate, inorganic carbon, total alkalinity and fCO₂. The changing pattern of ocean CO₂ summer sink and winter source is well reproduced by the model. It is controlled mainly by the balance between summer primary production and winter deep vertical mixing. In the POOZ, the annual air–sea CO₂ flux is about −0.5 mol/m²/yr, which is small compared to previous estimates based on oceanic observations but comparable to the small CO₂ sink deduced from atmospheric inverse methods. For reducing the uncertainties attached to the global ocean CO₂ sink south of the Polar Front the regional results presented here should be synthetized with historical and new observations, especially during winter, in other sectors of the Southern Ocean.     

Seasonal and interannual variability in particle fluxes of carbon,nitrogen and silicon from time series of sediment traps at Ocean Station P, 1982–1993: relationship to changes in subarctic primary productivity

An extended time series of particle fluxes at 3800 m was recorded using automated sediment traps moored at Ocean Station Papa (OSP, 50°N, 145°W) in the northeast Pacific Ocean for more than a decade (1982–1993). Time-series observations at 200 and 1000 m, and short-term measurements using surface-tethered free-drifting sediment traps also were made intermittently. We present data for fluxes of total mass (dry weight), particulate organic carbon (POC), particulate organic nitrogen (PON), biogenic Si (BSi), and particulate inorganic carbon (PIC) in calcium carbonate. Mean monthly fluxes at 3800 m showed distinct seasonality with an annual minimum during winter months (December–March), and maximum during summer and fall (April–November). Fluxes of total mass, POC, PIC and BSi showed 4-, 10-, 7- and 5-fold increases between extreme months, respectively. Mean monthly fluxes of PIC often showed two plateaus, one in May–August dominated by <63 μm particles and one in October–November, which was mainly >63 μm particles. Dominant components of the mass flux throughout the year were CaCO₃ and opal in equal amounts. The mean annual fluxes at 3800 m were 32±9 g dry weight g m⁻² yr⁻¹, 1.1±0.5 g POC m⁻² yr⁻¹, 0.15±0.07 g PON m⁻² yr⁻¹, 5.9±2.0 g BSi m⁻² yr⁻¹ and 1.7±0.6 g PIC m⁻² yr⁻¹. These biogenic fluxes clearly decreased with depth, and increased during “warm” years (1983 and 1987) of the El Niño, Southern Oscillation cycle (ENSO). Enhancement of annual mass flux rates to 3800 m was 49% in 1983 and 36% in 1987 above the decadal average, and was especially rich in biogenic Si. Biological events allowed estimates of sinking rates of detritus that range from 175 to 300 m d⁻¹, and demonstrate that, during periods of high productivity, particles sink quickly to deep ocean with less loss of organic components. Average POC flux into the deep ocean approximated the “canonical” 1% of the surface primary production.     

Role of the Agulhas Current in Indian Ocean circulation and associated heat and freshwater fluxes

A reduced estimate of Agulhas Current transport provides the motivation to examine the sensitivity of Indian Ocean circulation and meridional heat transport to the strength of the western boundary current. The new transport estimate is 70 Sv, much smaller than the previous value of 85 Sv. Consideration of three case studies for a large, medium and small Agulhas Current transport demonstrate that the divergence of heat transport over the Indian Ocean north of 32°S has a sensitivity of 0.08 PW per 10 Sv of Agulhas transport, and freshwater convergence has a sensitivity of 0.03×10⁹ kg s⁻¹ per 10 Sv of transport. Moreover, a smaller Agulhas Current leads to a better silica balance and a smaller meridional overturning circulation for the Indian Ocean. The mean Agulhas Current transport estimated from time-series current meter measurements is used to constrain the geostrophic transport in the western boundary region in order to re-evaluate the circulation, heat and freshwater transports across 32°S. The Indonesian Throughflow is taken to be 12 Sv at an average temperature of 18°C. The constrained circulation exhibits a vertical–meridional circulation with a net northward flow below 2000 dbar of 10.1 Sv. The heat transport divergence is estimated to be 0.66 PW, the freshwater convergence to be 0.54×10⁹ kg s⁻¹, and the silica convergence to be 335 kmol s⁻¹. Meridional transports are separated into barotropic, baroclinic and horizontal components, with each component conserving mass. The barotropic component is strongly dependent on the estimated size of the Indonesian Throughflow. Surprisingly, the baroclinic component depends principally on the large-scale density distribution and is nearly invariant to the size of the overturning circulation. The horizontal heat and freshwater flux components are strongly influenced by the size of the Agulhas Current because it is warmer and saltier than the mid-ocean. The horizontal fluxes of heat and salt penetrate down to 1500 m depth, suggesting that warm and salty Red Sea Water may be involved in converting the intermediate and upper deep waters which enter the Indian Ocean from the Southern Ocean into warmer and saltier waters before they exit in the Agulhas Current.     

Diapycnal nutrient fluxes on the northern boundary of Cape Ghir upwelling region

In this study we estimate diffusive nutrient fluxes in the northern region of Cape Ghir upwelling system (Northwest Africa) during autumn 2010. The contribution of two co-existing vertical mixing processes (turbulence and salt fingers) is estimated through micro- and fine-structure scale observations. The boundary between coastal upwelling and open ocean waters becomes apparent when nitrate is used as a tracer. Below the mixed layer (56.15±15.56 m), the water column is favorable to the occurrence of a salt finger regime. Vertical eddy diffusivity for salt (K_s) at the reference layer (57.86±8.51 m, CI 95%) was 3×10⁻⁵ (±1.89×10⁻⁹, CI 95%) m² s⁻¹. Average diapycnal fluxes indicate that there was a deficit in phosphate supply to the surface layer (6.61×10⁻⁴ mmol m⁻² d⁻¹), while these fluxes were 0.09 and 0.03 mmol m⁻² d⁻¹ for nitrate and silicate, respectively. There is a need to conduct more studies to obtain accurate estimations of vertical eddy diffusivity and nutrient supply in complex transitional zones, like Cape Ghir. This will provide us with information about salt and nutrients exchange in onshore–offshore zones.     

Lunar cycles and seasonal variations in deposition fluxes of planktic foraminiferal shell carbonate to the deep South Atlantic (central Walvis Ridge)

Several authors have argued that lunar reproductive cycling controls the shell fluxes of planktic foraminifera, one of the major carbonate-producing groups in the global pelagic ocean. A time-series sediment trap at 2700 m depth on the central Walvis Ridge below the South Atlantic central gyre demonstrate for the first time that shell deposition fluxes of Hastigerina pelagica are synchronous with lunar periodicity. Spectral analysis of the 6-month time-series with 8-day resolution showed a strong 30-day cyclicity in the flux maxima of H. pelagica arriving at the ocean floor on average 12.5 days after each full moon. Given a shell settling velocity of about 400 m day⁻¹, which implies about 7 days for settling, this coincides with the pronounced endogenous reproduction rhythm of 5±2 days after full moon as originally observed in laboratory-cultured isolates from off Bermuda in the North Atlantic. By contrast, no endogenous or exogenous lunar periodicity was observed in the deposition flux or size distribution of any of the 27 other shell species from austral winter (August 2000) to austral summer (February 2001). Instead, the deposition fluxes of shell species, the bulk carbonate and the total mass were dominated by a seasonal maximum during austral spring, without any periodicity in the 16–90-day domain of this study. Since H. pelagica exhibits low fluxes with a low burial efficiency, and continuous (re)production is shown by the deposition fluxes of other species, lunar reproductive cycling appears not to affect pelagic carbonate productivity and deep ocean sedimentation fluxes.     

Particulate organic carbon in the global ocean derived from SeaWiFS ocean color

Satellite remote sensing offers new means of quantifying particulate organic carbon, POC, concentration over large oceanic areas. From SeaWiFS ocean color, we derived 10-year data of POC concentration in the surface waters of the global ocean. The 10-year time series of the global and basin scale average surface POC concentration do not display any significant long-term trends. The annual mean surface POC concentration and its seasonal amplitude are highest in the North Atlantic and lowest in the South Pacific, when compared to other ocean basins. POC anomalies in the North Atlantic, North Pacific, and global concentrations seem to be inversely correlated with El Niño index, but longer time series are needed to confirm this relationship. Quantitative estimates of POC reservoir in the oceanic surface layer depend on the choice of what should represent this layer. Global average POC biomass is 1.34 g m^?2 if integrated over one optical depth, 3.62 g m^?2 if integrated over mixed layer depth, and up to 6.41 g m^?2 if integrated over 200-m layer depth (when assumed POC concentration below MLD is 20 mg m^?3). The global estimate of total POC reservoir in the surface 200-m layer of the ocean is 228.61×10¹³ g. We expect that future estimates of POC reservoir may be even larger, when more precise calculations account for deep-water organic-matter maxima in oligotrophic regions, and POC biomass located just below the seasonal mixed layer in spring and summer in the temperate regions.