Relationship between phytoplankton bloom and wind stress in the sub-polar frontal area of the Japan/East Sea

The seasonal dynamics of phytoplankton blooms in the central Japan/East Sea (JES) show pronounced year-to-year variability based on Sea-viewing Wide Field-of-view Sensor (SeaWiFS; 19972003) and Moderate Resolution Imaging Spectroradiometer (MODIS)/Terra (20002003) observations. Wind seems to strongly influence this variability. To study the relationship between wind and bloom initiation, we analyzed daily, remotely sensed wind stress data (Active Microwave Instrument–wind [AMI–wind], NASA Scatterometer [NSCAT], and Quick Scatterometer [QuickSCAT]: 19972003) and daily chlorophyll concentrations based on ocean color data (SeaWiFS and MODIS). The results agreed well with the hypotheses; in spring, blooms began 615 days after wind stress weakened. Fall blooms started 39 days after wind strengthened. We also simulated seasonal changes using a simple light–nutrient model using two values for the respiration ratio: 10% and 20%. The use of 20% seemed to reproduce the timing of the spring bloom quite well but underestimated the absolute level of chlorophyll concentration. On the other hand, using 10% produced a better estimation of the chlorophyll concentration but failed to match the timing. Neither of the model runs reproduced the timing of the fall bloom well.     

The interaction of tidal advection, diffusion and mussel filtration in a tidal channel

Time series measurements of flow and pigment concentrations (Chl) in the Menai Strait have revealed that the strong residual flow in a tidal channel ( 500 m³ s^− 1) transports phytoplankton from the open sea into the channel where much of it is consumed by suspension feeders, mainly in commercial beds of Mytilus edulis. The progressive depletion of phytoplankton along the channel results in a strong horizontal gradient of plankton and hence Chl. Tidal displacement of this gradient causes large (± 50% of mean) oscillations of Chl in the vicinity of the mussel beds. Vertical mixing by the strong tidal flows is sufficiently vigorous for most of the tidal cycle to ensure that downward diffusion can resupply the near-bed layer although there are indications of some transient depletion around slack water.This paradigm of the interaction of advection, diffusion and filtration determining the distribution of plankton and its supply to mussels has been encapsulated in a series of simple models forced only by boundary values. In the first, a 1-D model of tidal flow in the channel reproduces the principal features of the observed currents including the unusually large spatial change in phase of the currents and the variation of the residual transport with tidal range. The flow field from this physical model is used to drive a second model based on the advection diffusion equation for Chl with a source at the Irish Sea boundary and a sink over the mussel bed. This model illustrates the formation of a strong Chl gradient along the channel and simulates the amplitude and phase of the M₂ oscillations of Chl and the development of the M₄ variation apparent in the observations. This second model has been extended to 2-D over the mussel beds to allow investigation of the effects of water column mixing. The model indicates that only for a short period ( 30 min), close to slack water, is mixing sufficiently reduced to permit the development of a depletion boundary layer and then only within  1 m from the bottom, a result which is consistent with the observations.     

Validation of a hybrid optimal interpolation and Kalman filter scheme for sea surface temperature assimilation

A hybrid data assimilation scheme designed for operational assimilation of satellite sea surface temperatures (SST) into an ocean model has been developed and validated against in-situ observations. The scheme consists of an optimal interpolation (OI) part and a greatly simplified Kalman filter (KF) part.The OI is performed only in the longitudinal and latitudinal directions. A climatological field is used as a background field for the interpolation. It is constructed by fitting daily averages of satellite SST to the annual mean, annual, and semiannual harmonics in a 20 km by 20 km grid. The background error covariance is approximated by a spatially varying two-dimensional exponential covariance model. The parameters of the covariance model are fitted to the deviations of the satellite data from the background field using data from a full year.The simplified KF uses ocean model forecasts as a background field. It is based on the assumption that it is possible to neglect horizontal SST covariances in the filter and that the typical time scale for vertical mixing in the mixed layer is much shorter than the average time between observations. We therefore assume that the error variance in a column of water is evenly spread out throughout the mixed layer. The result of these simplifications is a computationally very efficient KF.A one year validation of the scheme is performed for year 2001 using an operational eddy resolving ocean model covering the North Sea and the Baltic Sea. It is found that assimilation of sea surface temperature data reduces the model root mean square error from 1.13 °C to 0.70 °C. The hybrid scheme is found to reduce the root mean square error slightly more than the simplified KF without OI to 0.66 °C. The inclusion of spatially varying satellite error variances does not improve the performance of the scheme significantly.     

The Gulf Stream, rings and North Atlantic eddy structures from remote sensing (Altimeter and SeaWiFS)

Seasonal SeaWiFS chlorophyll a concentrations cycles and annual changes of altimeter Sea Level Anomaly are derived for the subtropical North Atlantic near  35°N and along a Gulf Stream axis. Spatial structure of SeaWiFS, is defined in terms of deviations from a local seasonal cycle and examined in relation to altimeter eddy structure. In the subtropical region near 35°N, SeaWiFS structure is evident during the spring bloom period with a scale of  430 km, or about twice the eddy scale. A Gulf Stream axis has been selected as a region where the Sea Level Anomaly variance is a maximum. Eddy propagation speeds and scales are examined. Cold-core (cyclonic) rings correspond to areas of high SeaWiFS chlorophyll a. Warm-core (anticyclonic) rings relate to areas of low chlorophyll concentration. Both SeaWiFS structure and eddy structure have a spatial scale of  450 km or twice the ring scale along the Gulf Stream axis. SeaWiFS chlorophyll anomalies and Altimeter Sea Level Anomaly structure have an overall negative correlation coefficient of r = − 0.34. Swirl currents between eddies redistribute surface chlorophyll concentrations and can spatially bias maximum and minimum concentration levels off eddy centre.     

On the key role of nutrient data to constrain a coupled physical–biogeochemical assimilative model of the North Atlantic Ocean

A sequential assimilative system has been implemented into a coupled physical–biogeochemical model (CPBM) of the North Atlantic basin at eddy-permitting resolution (1/4°), with the long-term goal of estimating the basin scale patterns of the oceanic primary production and their seasonal variability. The assimilation system, which is based on the SEEK filter [Brasseur, P., Verron, J., 2006. The SEEK filter method for data assimilation in oceanography: a synthesis. Ocean Dynamics. doi: 10.1007/s10236-006-0080-3], has been adapted to this CPBM in order to control the physical and biogeochemical components of the coupled model separately or in combination. The assimilated data are the satellite Topex/Poseidon and ERS altimetric data, the AVHRR Sea Surface Temperature observations, and the Levitus climatology for salinity, temperature and nitrate.In the present study, different assimilation experiments are conducted to assess the relative usefulness of the assimilated data to improve the representation of the primary production by the CPBM. Consistently with the results obtained by Berline et al. [Berline, L., Brankart, J-M., Brasseur, P., Ourmières, Y., Verron, J., 2007. Improving the physics of a coupled physical–biogeochemical model of the North Atlantic through data assimilation: impact on the ecosystem. J. Mar. Syst. 64 (1–4), 153–172] with a comparable assimilative model, it is shown that the assimilation of physical data alone can improve the representation of the mixed layer depth, but the impact on the ecosystem is rather weak. In some situations, the physical data assimilation can even worsen the ecosystem response for areas where the prior nutrient distribution is significantly incorrect. However, these experiments also show that the combined assimilation of physical and nutrient data has a positive impact on the phytoplankton patterns by comparison with SeaWiFS ocean colour data, demonstrating the good complementarity between SST, altimetry and in situ nutrient data. These results suggest that more intensive in situ measurements of biogeochemical nutrients are urgently needed at basin scale to initiate a permanent monitoring of oceanic ecosystems.     

A dense water outflow from the Ross Sea, Antarctica: Mixing and the contribution of tides

We use hydrographic, current, and microstructure measurements, and tide-forced ocean models, to estimate benthic and interfacial mixing impacting the evolution of a bottom-trapped outflow of dense shelf water from the Drygalski Trough in the northwestern Ross Sea. During summer 2003 an energetic outflow was observed from the outer shelf ( 500 m isobath) to the  1600 m isobath on the continental slope. Outflow thickness was as great as  200 m, and mean speeds were  0.6 m s^− 1 relative to background currents exceeding  1 m s^− 1 that were primarily tidal in origin. No outflow was detected on the slope in winter 2004, although a thin layer of dense shelf water was present on the outer shelf. When the outflow was well-developed, the estimated benthic stress was of order one Pascal and the bulk Froude number over the upper slope exceeded one. Diapycnal scalar diffusivity (K_z) values in the transition region at the top of the outflow, estimated from Thorpe-scale analysis of potential density and measurements of microscale temperature gradient from sensors attached to the CTD rosette, were of order 10^− 3−10^− 2 m² s^− 1. For two cases where the upper outflow boundary was particularly sharply defined, entrainment rate w_e was estimated from K_z and bulk outflow parameters to be  10^− 3 m s^− 1 ( 100 m day^− 1). A tide-forced, three-dimensional primitive equation ocean model with Mellor-Yamada level 2.5 turbulence closure scheme for diapycnal mixing yields results consistent with a significant tidal role in mixing associated with benthic stress and shear within the stratified ocean interior.     

Assessment of sea surface temperature observational networks in the Baltic Sea and North Sea

The satellite and in situ Sea Surface Temperature (SST) observational networks in the Baltic Sea and North Sea are evaluated based on the quality of the gridded SST products generated from the networks. A multi-indicator approach is applied in the assessment. It includes evaluation of data quality, effective data coverage, field reconstruction error and model nowcast error. The results show that the best available full-coverage SST product is generated by assimilating the SST observations to obtain a yearly mean model bias of 0.07 °C and RMSE of 0.64 °C. The effective data coverage rate is 31% by using AVHRR (Advanced Very High Resolution Radiometer) data from NOAA (National Ocean and Atmosphere Administration) satellites 12, 14 and 16. The data redundancy increases rapidly with the number of infrared sensors. Using either NOAA satellite 12 or all 3 satellites makes a small difference with regard to derived effective coverage and the ocean model nowcast error. The influence of using the in situ SST observations in the SST field reconstruction is negligibly small. Instead, the major role of in situ SST observations is in calibrating the satellite observations. To study the relative importance of data quality and data coverage, an assessment is done for two satellite products: one product is based entirely on NOAA 12 data and has larger coverage but lower quality. The other product is a subset of the SAF products (derived from NOAA 14 and 16) and has lower coverage but higher quality. Based on current monitoring, modelling and assimilation technology, the results suggest that the data quality is an important factor in further improving the quality of the gridded SST products. Recommendations are made for possible further improvements of the existing SST observational networks.     

Improving the physics of a coupled physical–biogeochemical model of the North Atlantic through data assimilation: Impact on the ecosystem

Several studies on coupled physical–biogeochemical models have shown that major deficiencies in the biogeochemical fields arise from the deficiencies in the physical flow fields. This paper examines the improvement of the physics through data assimilation, and the subsequent impact on the ecosystem response in a coupled model of the North Atlantic. Sea surface temperature and sea surface height data are assimilated with a sequential method based on the SEEK filter adapted to the coupling needs. The model domain covers the Atlantic from 20°S to 70°N at eddy-permitting resolution. The biogeochemical model is a NPZD-DOM model based on the P3ZD formulation. The results of an annual assimilated simulation are compared with an annual free simulation.With assimilation, the representation of the mixed layer depth is significantly improved in mid latitudes, even though the mixed layer depth is generally overestimated compared to the observations. The representation of the mean and variance of the currents is also significantly improved.The nutrient input in the euphotic zone is used to assess the data assimilation impact on the ecosystem. Data assimilation results in a 50% reduction of the input due to vertical mixing in mid-latitudes, and in a four- to six-fold increase of the advective fluxes in mid-latitudes and subtropics. Averaged zonally, the net impact is a threefold increase for the subtropical gyre, and a moderate (20–30%) decrease at mid and high latitudes.Surface chlorophyll concentration increases along the subtropical gyre borders, but little changes are detected at mid and high latitudes. An increase of the primary production appears along the Gulf Stream path, but it represents only 12% on average for mid and high latitudes. In the subtropical gyre centre, primary production is augmented but stays underestimated (20% of observations). These experiments show the benefits of physical data assimilation in coupled physical–biogeochemical applications.     

A trace metal (Pb, Cd, Zn, Cu) balance for surface waters in the eastern Gotland Basin, Baltic Sea

Trace metal inventories and transboundary fluxes were estimated for the surface waters of the eastern Gotland Basin, and a mass balance for this water body was derived.The study area was delineated by a box which bordered vertically by the sea surface and the halocline and horizontally by the coast of Gotland and the coast of Latvia/Lithuania. For the calculation of the trace metal inventories in the box, a high-resolution sampling was carried out during different seasons. Additional sediment trap studies were performed to calculate the vertical particulate trace metal fluxes through the halocline.The following transboundary fluxes have been estimated: (a) atmospheric input of trace metals; (b) lateral transport from adjacent bodies of surface water; (c) diapycnal fluxes of dissolved trace metals at the halocline; (d) sinking of trace metals, associated with particles, through the halocline, (e) vertical advection through the halocline.The lateral transport into and out of the surface box is important for the metals (Cd, Cu, Zn) with “nutrient-like” behaviour. For particle reactive elements like lead, vertical sedimentation and lateral transport out of the box as much as atmospheric input and lateral input into the box are in the same order of magnitude. Turbulent diapycnal mixing plays a minor role compared to the advective fluxes.The estimated range for the residence times for the region under investigation are  0.5–1.3 years for Pb;  2–11 years for Zn,  7–36 years for Cd and 42–89 years for Cu. This demonstrates that the system reacts very fast for particle reactive elements like Pb, while for copper sedimentation processes are not the preferential sink and can be neglected.     

Numerical simulation of the barotropic tides in the Tunisian Shelf and the Strait of Sicily

We have investigated the barotropic tides in the Tunisian shelf and the Strait of Sicily using the Regional Ocean Modelling System (ROMS) with very high-resolution. Model performance was evaluated with respect to tide gauge, satellite data, and current meter measurements. The model fields faithfully reproduced the major feature of the barotropic tidal currents and agreed well with existing tidal elevation and phase observations. General features for the various semidiurnal constituents are nearly similar to each other with maximum amplitude in the Gulf of Gabes. The larger tidal currents occur over the continental shelves. In the Adventure Bank, the current is essentially of diurnal type whereas in the Gulf of Gabes it is of semidiurnal type.Tidal energy lost, which is primarily due to bottom stress dissipation, is predominantly in the Gulf of Gabes ( 61%), the Strait of Sicily, and the Strait of Messina. The forcing function for internal tides shows for both M₂ and K₁ constituents, significant spatial variability in the Strait of Sicily. This suggests that some internal tides will be generated in these regions and could thus explain the observed strong diurnal internal waves in the Adventure Bank.     

Application of a SEEK filter to a 1D biogeochemical model of the Ligurian Sea: Twin experiments and real in-situ data assimilation

The Singular Evolutive Extended Kalman (SEEK) filter has been implemented to assimilate in-situ data in a 1D coupled physical-ecosystem model of the Ligurian Sea. The biogeochemical model describes the partly decoupled nitrogen and carbon cycles of the pelagic food web. The GHER hydrodynamic model (1D version) is used to represent the physical forcings. The data assimilation scheme (SEEK filter) parameterizes the error statistics by means of a set of empirical orthogonal functions (EOFs). Twin experiments are first performed with the aim to choose the suitable experimental protocol (observation and estimation vectors, number of EOFs, frequency of the assimilation,…) and to assess the SEEK filter performances. This protocol is then applied to perform real data assimilation experiments using the DYFAMED data base. By assimilating phytoplankton observations, the method has allowed to improve not only the representation of the phytoplankton community, but also of other variables such as zooplankton and bacteria that evolve with model dynamics and that are not corrected by the data assimilation scheme. The validation of the assimilation method and the improvement of model results are studied by means of suitable error measurements.     

1/32° real-time global ocean prediction and value-added over 1/16° resolution

A 1/32° global ocean nowcast/forecast system has been developed by the Naval Research Laboratory at the Stennis Space Center. It started running at the Naval Oceanographic Office in near real-time on 1 Nov. 2003 and has been running daily in real-time since 1 Mar. 2005. It became an operational system on 6 March 2006, replacing the existing 1/16° system which ceased operation on 12 March 2006. Both systems use the NRL Layered Ocean Model (NLOM) with assimilation of sea surface height from satellite altimeters and sea surface temperature from multi-channel satellite infrared radiometers. Real-time and archived results are available online at http://www.ocean.nrlssc.navy.mil/global_nlom. The 1/32° system has improvements over the earlier system that can be grouped into two categories: (1) better resolution and representation of dynamical processes and (2) design modifications. The design modifications are the result of accrued knowledge since the development of the earlier 1/16° system. The improved horizontal resolution of the 1/32° system has significant dynamical benefits which increase the ability of the model to accurately nowcast and skillfully forecast. At the finer resolution, current pathways and their transports become more accurate, the sea surface height (SSH) variability increases and becomes more realistic and even the global ocean circulation experiences some changes (including inter-basin exchange). These improvements make the 1/32° system a better dynamical interpolator of assimilated satellite altimeter track data, using a one-day model forecast as the first guess. The result is quantitatively more accurate nowcasts, as is illustrated by several model-data comparisons. Based on comparisons with ocean color imagery in the northwestern Arabian Sea and the Gulf of Oman, the 1/32° system has even demonstrated the ability to map small eddies, 25–75 km in diameter, with 70% reliability and a median eddy center location error of 22.5 km, a surprising and unanticipated result from assimilation of altimeter track data. For all of the eddies (50% small eddies), the reliability was 80% and the median eddy center location error was 29 km. The 1/32° system also exhibits improved forecast skill in relation to the 1/16° system. This is due to (a) a more accurate initial condition for the forecast and (b) better resolution and representation of critical dynamical processes (such as upper ocean – topographic coupling via mesoscale flow instabilities) which allow the model to more accurately evolve these features in time while running in forecast mode (forecast atmospheric forcing for the first 5 days, then gradually reverting toward climatology for the remainder of the 30-day forecast period). At 1/32° resolution, forecast SSH generally compares better with unassimilated observations and the anomaly correlation of the forecast SSH exceeds that from persistence by a larger amount than found in the 1/16° system.     

Microscale variability in the distributions of large fluorescent particles observed in situ with a planar laser imaging fluorometer

It has been known for decades that particle-size and biomass spectra show regular patterns in the ocean, and that these patterns often show systematic variations with other properties such as total biomass, nutrient concentration, season, and distance (both vertical and horizontal). The recent finding of the ubiquitous nature of layers of phytoplankton < 1 m thick prompted us to explore the fine- and microscale vertical variations of size- and fluorescence-abundance spectra in the ocean. Using a two-dimensional planar laser imaging system mounted on a free-falling platform, we quantified the properties of large fluorescent particles ( 20 μm–2 cm) through the water column, obtaining images every 10–30 cm. These images showed systematic relationships of the spectral properties to total chlorophyll: increased proportions of the smallest particles at high chlorophyll concentrations, and a lengthening of the spectral size range at high total chlorophyll concentrations (more large particles at high chlorophyll concentrations). Further, we observed significant variations of the spectral properties over scales of 1 m and less, and recorded the frequent occurrence of unusual layers of large particles. Our new instrument, which is sensitive to thin layers of enhanced phytoplankton biomass, shows the planktonic community to be highly structured vertically on scales of 1–2 m, particularly within the DCM.     

Circulation, hydrographic and nutrient characteristics of the Cilician Basin, Northeastern Mediterranean Sea

The water mass, circulation and chemical properties of the Cilician Basin, the northeastern Levantine Sea, are described on the basis of three hydrographic cruises performed during May 1997 (spring), July 1998 (summer) and October 2003 (autumn). The hydrographic data reveal the presence of Levantine Surface Water (LSW) and Modified Atlantic Water (MAW) within the upper 90 m layer, Levantine Intermediate Water (LIW) between 90 and 250 m, and Transitional Mediterranean Water (TMW) further below. The temporal variability of the circulation system is manifested by a change in shape, size and intensity of eddies as well as the pathways of the Lattakia Basin coastal current system. The nutrient concentrations varied between nitrate + nitrite = 0.16–0.31 μM, phosphate = 0.02–0.03 μM and silicate = 0.95–1.2 μM for the surface layer during sampling periods. Dissolved nutrient concentrations in the Transitional Mediterranean Water were: 2.1–5.3 μM for NO₃ + NO₂, 0.10–0.21 μM for PO₄ and 5.7–10 μM for Si. The molar ratios of nitrate to phosphate in the water column range between 5 and 20 in the surface layer and reach up to a value of 45 at the top of the nutricline at the depths of 29.05 kg/m³ isopycnal surface for most of the year. Below the nutricline the N / P ratios retain the values around 24–28.     

The Mercator global ocean operational analysis system: Assessment and validation of an 11-year reanalysis

This paper presents Prototype Système 2 Global (PSY2G), the first Mercator global Ocean General Circulation Model (OGCM) to assimilate along-track sea level anomaly (SLA) satellite data. Based on a coarse resolution ocean model, this system was developed mainly for climatic purposes and will provide, for example, initial oceanic states for coupled ocean-atmosphere seasonal predictions. It has been operational since 3 September 2003 and produces an analysis and a two-week forecast for the global ocean every week. The PSY2G system uses an incremental assimilation scheme based on the Cooper and Haines [Cooper, M., Haines, K., 1996. Data assimilation with water property conservation. J. Geophys. Res., 101, 1059-1077.] lifting–lowering of isopycnals. The SLA increment is obtained using an optimal interpolation method then the correction is partitioned into baroclinic and barotropic contributions. The baroclinic ocean state correction consists of temperature, salinity and geostrophic velocity increments and the barotropic correction is a barotropic velocity increment. A reanalysis (1993–2003) was carried out that enabled the PSY2G system to perform its first operational cycle. All available SLA data sets (TOPEX/Poséïdon, ERS2, Geosat-Follow-On, Jason1 and Envisat) were assimilated for the 1993–2003 period. The major objective of this study is to assess the reanalysis from both an assimilation and a thermodynamic point of view in order to evaluate its realism, especially in the tropical band which is a key region for climatic studies. Although the system is also able to deliver forecasts, we have mainly focused on analysis. These results are useful because they give an a priori estimation of the qualities and capabilities of the operational ocean analysis system that has been implemented. In particular, the reanalysis identifies some regional biases in sea level variability such as near the Antarctic Circumpolar Current, in the eastern Equatorial Pacific and in the Norwegian Sea (generally less than 1 cm) with a small seasonal cycle. This is attributed to changes in mean circulation and vertical stratification caused by the assimilation methodology. But the model's low resolution, inaccurate physical parameterisations (especially for ocean–ice interactions) and surface atmospheric forcing also contribute to the occurrence of the SLA biases. A detailed analysis of the thermohaline structure of the ocean reveals that the isopycnal lifting–lowering tends to diffuse vertically the main thermocline. The impact on temperature is that the surface layer (0–200 m) becomes cooler whereas in deeper waters (from 500 to 1500 m), the ocean becomes slightly warmer. This is particularly true in the tropics, between 30°N and 30°S. However it can be demonstrated that the assimilation improves the variability in both surface currents and sub-surface temperature in the Equatorial Pacific Ocean.     

Improving assimilation of SeaWiFS data by the application of bias correction with a local SEIK filter

Ocean-biogeochemical models show typically significant errors in the representation of chlorophyll concentrations. The model state can be improved by the assimilation of satellite chlorophyll data with algorithms based on the Kalman filter. However, these algorithms do usually not account for the possibility that the model prediction contains systematic errors in the form of model bias. Accounting explicitly for model biases can improve the assimilation performance. To study the effect of bias estimation on the estimation of surface chlorophyll concentrations, chlorophyll data from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) are assimilated on a daily basis into the NASA Ocean Biogeochemical Model (NOBM). The assimilation is performed by the ensemble-based SEIK filter combined with an online bias correction scheme. The SEIK filter is simplified here by the use of a static error covariance matrix. The performance of the filter algorithm is assessed by comparison with independent in situ data over the 7-year period 1998–2004. The bias correction results in significant improvements of the surface chlorophyll concentrations compared to the assimilation without bias estimation. With bias estimation, the daily surface chlorophyll estimates from the assimilation show about 3.3% lower error than SeaWiFS data. In contrast, the error in the global surface chlorophyll estimate without bias estimation is 10.9% larger than the error of SeaWiFS data.     

Observations of the Fawn Trough Current over the Kerguelen Plateau from instrumented elephant seals

Due to its great meridional extent and relatively shallow depths, the Kerguelen Plateau constitutes a major barrier to the eastward flowing Antarctic Circumpolar Current in the Indian sector of the Southern Ocean. While most of the Antarctic Circumpolar Current transport is deflected north of the Kerguelen Islands, the remainder ( 50 Sv, 1 Sv = 10⁶ m³ s^− 1) must pass south of the islands, most probably through the Fawn and Princess Elizabeth Troughs. However, the paucity of finely resolved quasi-synoptic hydrographic data in this remote and infrequently sampled area has limited the progress in our knowledge of the regional circulation. Since 2004, a new approach using elephant seals from the Kerguelen Islands as autonomous oceanographic profilers has provided new information on the hydrography over the Kerguelen Plateau, covering the entire Antarctic Zone between the Polar Front and Antarctica, with a mean along-track resolution of about 25 km. These finely resolved bio-logged data revealed details of a strong northeastward current found across the Fawn Trough (sill depth: 2600 m; 56°S, 78°E). This so-called Fawn Trough Current transports cold Antarctic waters found mostly south of the Elan Bank, between the Ice Limit (58°S) and the Antarctic Divergence (64°S) in the eastern Enderby Basin, toward the Australian–Antarctic Basin. Our analysis also demonstrates that the Deep Western Boundary Current, which carries cold Antarctic water along the eastern flank of the southern Kerguelen Plateau collides with Fawn Trough Current at the outlet of the Fawn Trough sill. In other words, the Fawn Trough constitutes a veritable bottleneck, channelling the quasi-totality of the Antarctic Circumpolar flow found south of the Polar Front. Thanks to the unprecedented fine resolution of seal-borne data, a branch of flow centered at the Winter Water isotherm of 1 °C is also revealed along the northern escarpment of the Elan Bank, and then along the southern edge of Heard Island. Further analysis of different supplementary data reveals a complex circulation pattern in the entire Enderby Basin, with several distinctive branches of flow being strongly controlled by prominent topographic features such as the Southwest Indian Ridge, Conrad Rise, Elan Bank, and Kerguelen Plateau. This newly emerged frontal structure refines considerably previous large-scale circulation schematics of the area.     

Multigrid state vector for data assimilation in a two-way nested model of the Ligurian Sea

A system of two nested models composed by a coarse resolution model of the Mediterranean Sea, an intermediate resolution model of the Provençal Basin and a high resolution model of the Ligurian Sea is coupled with a Kalman-filter based assimilation method. The state vector for the data assimilation is composed by the temperature, salinity and elevation of the three models. The forecast error is estimated by an ensemble run of 200 members by perturbing initial condition and atmospheric forcings. The 50 dominant empirical orthogonal functions (EOF) are taken as the error covariance of the model forecast. This error covariance is assumed to be constant in time. Sea surface temperature (SST) and sea surface height (SSH) are assimilated in this system.     

Dissolved organic matter in shelf waters off the Ría de Vigo (NW Iberian upwelling system)

Net in situ production and export of dissolved organic carbon (DOC) and nitrogen (DON) have been studied in shelf waters off the Ría de Vigo (NW Spain), as part of a comprehensive hydrographic survey carried out from September 1994 to September 1995 with a fortnight periodicity. DOC and DON correlated well (r=+0.78), the slope of the regression line being 12.0±0.7 mol-C mol-N⁻¹, about twice the Redfieldian slope of particulate organic matter, 6.5±0.2 mol-C mol-N⁻¹ (r=+0.95). Labile DOC and DON accumulated in the upper 50 m during the upwelling season (March–September), mainly after prolonged periods of wind relaxation, when horizontal flows were reduced. This labile material represented 50% and 35% of the total (dissolved+particulate) organic carbon and nitrogen susceptible of microbial utilisation, which assert the key contribution of dissolved organic matter (DOM) to the export of new primary production in the NW Iberian upwelling system. This surface excess in shelf waters appeared to be formed into the highly productive Ría de Vigo (a large coastal indentation) at net rates of 4.4 μM-C d⁻¹ and 1.3 μM-C d⁻¹ in the inner and outer segments of the embayment respectively, and subsequently exported to the shelf. Once in the shelf, simple dilution with the inert DOM pool of recently upwelled Eastern North Atlantic Central Water (ENACW) occurred. Eventually, the DOM excess produced during the upwelling season is exported to the adjacent open ocean waters by the coastal circulation. Conversely, during the unproductive downwelling season (October–February), the lowest DOC and DON levels were recorded and export was prevented by the characteristic downwelling front associated to the seasonal poleward slope current.     

The biogeochemical cycling of methane in Ria de Vigo, NW Spain: Sediment processing and sea–air exchange

Methane (CH₄) concentrations were measured in the water column, in sediment porewaters, and in atmospheric air, in the Ría de Vigo, NW Spain, during both the onset (April 2003) and at the end of (September 2004) seasonal upwelling. In addition, CH₄ concentration and stable isotopic signatures (δ¹³CH₄) were measured in porewaters, and sediment methanogenesis and aerobic oxidation of CH₄ were determined in sediment incubations. Surface water column CH₄ (2 m depth) was in the range 3–180 nmol l^− 1 (110–8500% saturation) and followed a generally landward increase but with localised maxima in both the inner and middle Ría. These maxima were consistent with CH₄ inputs from underlying porewaters in which CH₄ concentrations were up to 3 orders of magnitude higher (maximum 350 μmol l^− 1). Surface water CH₄ concentrations were approximately three times higher in September than in April, consistent with a significant benthic CH₄ flux driven by enhanced sediment methanogenesis following the summer productivity maximum. CH₄ and δ¹³CH₄ in sediment porewaters and in incubated sediment slurries (20 °C) revealed significant sediment CH₄ oxidation, with an apparent isotopic fractionation factor (r_c) of  1.004. Using turbulent diffusion models of air–sea exchange we estimate an annual emission of atmospheric CH₄ from the Ría de Vigo of 18–44 × 10⁶ g (1.1–2.7 × 10⁶ mol). This estimate is approximately 1–2 orders of magnitude lower than a previous estimate based on a bubble transport model.