A Mean Dynamic Topography of the Mediterranean Sea computed from altimetric data, in-situ measurements and a general circulation model

In the Mediterranean Sea, where the mean circulation is largely unknown and characterized by smaller scales and less intensity than in the open ocean, the interpretation of altimetric Sea Level Anomalies (SLA) is rather difficult. In the context of operational systems such as MFS (Mediterranean Forecasting System) or MERCATOR, that assimilate the altimetric information, the estimation of a realistic Mean Dynamic Topography (MDT) consistent with altimetric SLA to be used to reconstruct absolute sea level is a crucial issue. A method is developed here to estimate the required MDT combining oceanic observations as altimetric and in-situ measurements and outputs from an ocean general circulation model (OGCM).In a first step, the average over the 1993–1999 period of dynamic topography outputs from MFS OGCM provides a first guess for the computation of the MDT. Then, in a second step, drifting buoy velocities and altimetric data are combined using a synthetic method to obtain local estimates of the mean geostrophic circulation which are then used to improve the first guess through an inverse technique and map the MDT field (hereafter the Synthetic Mean Dynamic Topography or SMDT) on a 1/8° resolution grid.Many interesting current patterns and cyclonic/anticyclonic structures are visible on the SMDT obtained. The main Mediterranean coastal currents are well marked (as the Algerian Current or the Liguro–Provenço–Catalan Current). East of the Sicily channel, the Atlantic Ionian Stream divides into several main branches crossing the Ionian Sea at various latitudes before joining at 19°E into a unique Mid-Mediterranean Jet. Also, strong signatures of the main Mediterranean eddies are obtained (as for instance the Alboran gyre, the Pelops, Ierapetra, Mersa-Matruh or Shikmona anticyclones and the Cretan, Rhodes or West Cyprius cyclones). Independent in-situ measurements from Sea Campaigns NORBAL in the North Balearic Sea and the North Tyrrhenian Sea and SYMPLEX in the Sicily channel are used to validate locally the SMDT: deduced absolute altimetric dynamic topography compares well with in-situ observations. Finally, the SMDT is used to compute absolute altimetric maps in the Alboran Sea and the Algerian Current. The use of absolute altimetric signal allows to accurately follow the formation and propagation of cyclonic and anticyclonic eddies in both areas.     

Seasonal study of dissolved CH4, CO2 and N2O in a shallow tidal system of the bay of Cádiz (SW Spain)

During 2004, 10 samplings were performed in order to measure dissolved methane (CH₄), carbon dioxide (CO₂) and nitrous oxide (N₂O) in the surface waters of Río San Pedro, a tidal creek in the salt marsh area of the Bay of Cádiz (SW Spain). The inner partvs of the creek is affected by the inputs coming from an intensive fish farm and the drainage of an extensive salt marsh area.Dissolved CH₄, CO₂ and N₂O concentrations ranged from 11 to 88 nM, 36 to 108 μM and 14 to 50 nM, respectively. Surface waters were in all cases oversaturated with respect to the atmosphere, reaching values of up to 5000% for CH₄, 1240% for CO₂ and 840% for N₂O. Dissolved CH₄, CO₂ and N₂O showed a significant tidal and seasonal variability. Over a tidal cycle, concentrations were always highest during low tide, which points to the influence of the inputs from the fish farm effluent and the drainage of the adjacent salt marsh area, as well as in situ production within the system. Dissolved CH₄, CO₂ and N₂O seasonal patterns were similar and showed maximum concentrations in summer conditions. Using four different parameterizations to calculate the gas transfer coefficients [Liss, P.S. and Merlivat, L., 1986. Air-sea exchange rates: introduction and synthesis. In P. Buat-Ménard (Ed.), The Role of Air-Sea Exchanges in Geochemical Cycling. Reidel, Dordrecht, The Netherlands, p. 113–127.; Clark, J.F., Schlosser, P., Simpson, H.J., Stute, M., Wanninkhof, R., and Ho, D.T., 1995. Relationship between gas transfer velocities and wind speeds in the tidal Hudson River determined by the dual tracer technique. In: B. Jähne and E. Monahan (Eds.), Air-Water Gas Transfer: AEON Verlag and Studio, Hanau, Germany, pp. 785–800.; Carini, S., Weston, N., Hopkinson, G., Tucker, J., Giblin, A. and Vallino, J., 1996. Gas exchanges rates in the Parker River estuary, Massachusetts. Biol. Bull., 191: 333–334.; Kremer, J.N., Reischauer, A. and D'Avanzo, C., 2003. Estuary-specific variation in the air-water gas exchange coefficient for oxygen. Estuaries, 26: 829–836.], the averaged air–water fluxes of CH₄, CO₂ and N₂O from the creek to the atmosphere ranged between 34 and 150 μmol CH₄ m^− 2 day^− 1, 73 and 177 mmol CO₂ m^− 2 day^− 1 and 24 and 62 μmol N₂O m⁻² day⁻¹, respectively.     

Parameterisation of clastic sediments including benthic structures

The sediment transport processes in the south-western Baltic Sea are predicted by means of a numerical model in the project DYNAS. There are two sediment parameters that influence the results of modelling remarkably: critical shear stress velocity and bottom roughness. This paper presents the way how to parameterise these factors and extrapolate them into the investigation area. The critical shear stress velocity is parameterised basing on grain size data, combining approximations after Hjulström [Hjulström, F., 1935: Studies in the morphological activity of rivers as illustrated by the river Fyris. Geological Institution of University of Uppsala: Bulletin (25): 221–528.], Shields [Shields, A., 1936: Anwendung der Ähnlichkeits-Mechanik und der Turbulenzforschung auf die Geschiebebewegung. Mitteilungen der Preussischen Versuchsanstalt für Wasserbau und Schiffahrt (26): 26 pp.] and Bohling [Bohling, B., 2003: Untersuchungen zur Mobilität natürlicher und anthropogener Sedimente in der Mecklenburger Bucht. unpublished doctoral thesis, Mathematisch-Naturwissenschaftliche Fakultät, Ernst-Moritz-Arndt-Universität Greifswald/Germany, 156 pp.]. The roughness length, in the case of absence of macro zoo-benthos and their structures, is parameterised basing on grain size too employing Soulsby [Soulsby, R.L., 1997: Dynamics of Marine Sands: a Manual for Practical Applications. London, Thomas Telford Publications. 249 pp.], Nielsen [Nielsen, P., 1983: Analytical determination of nearshore wave height variation due to refraction shoaling and friction. Coastal Engineering 7, 233–251.] and Yalin [Yalin, M.S., 1977: Mechanics of Sediment Transport. Pergamon Press, New York. 298 pp.]. No equivalent simple parameterisations for biologically caused bed roughness exist. Here, findings of Friedrichs [Friedrichs, M., 2004: Flow-induced effects of macro zoo-benthic structures on the near-bed sediment transport. Dissertation, Universität Rostock, 80 S.] and estimations by the DYNAS biologists group were combined in order to derive roughness lengths from abundance measurements of four previously selected key species which represent the originators of the dominating benthic structures at the sea floor in the south-western Baltic Sea. Critical shear stress velocity and bed roughness are known at few sample sites only. They were extrapolated into the larger investigation area using a proxy-target concept. The mean near bottom milieu (bathymetry, median grain size, salinity, oxygen) which was derived using results from numerical modelling serves as the proxy. Since the milieu parameters are measured at the sampling sites for which the target parameters have been determined, a combined hierarchical and supervised classification was employed to transfer the local knowledge into the unknown investigation area.     

Recent sedimentary study of the shelf of the Basque country

The Northern Iberian margin of the Spanish Basque country (provinces of Gipuzkoa and Viscaia) is characterized by a narrow continental platform, which receives inputs of riverine particulate matter from the numerous riverine systems located within the Basque country. This particulate matter is subsequently deposited within the Bay of Biscay, and Gouf de Capbreton [Frouin, R., Fiuza, A.F.G., Ambar, I., Boyd, T.J., 1990. Observations of a poleward surface current off the coasts of Portugal and Spain during winter. Journal of Geophysical Research 95 (C1), 679–691]. The main goal of this study is to establish a map of the surface sediment distribution of the Basque continental shelf and more specifically to map the muddy patch located at the eastern side of that continental shelf.Three oceanographic cruises were conducted in 2003 and 2004. From these campaigns 340 surface samples, 12 short cores and 3 gravity cores have been collected over the mid and outer shelf from depths ranging between 50 m and 150 m deep. 3 seismic profiles were obtained across the shelf mud patch using a Sparker device.Sediment grain-size analyses were performed by the classical physical method of sieving and use of settling columns. The POC (Particular Organic Carbon) amounts in sediment and water samples were determined using the Strickland and Parsons' method [Strickland, J.D.H., Parsons, T.R., 1972. Determination of particulate carbon. In : A practical handbook of seawater analysis. Fisheries ResearchBoard of Canada, Ottawa, pp. 207−211] as adapted by Etcheber [Etcheber, H., 1981. Comparaison des diverses méthodes d'évaluation des teneurs en matières en suspension et en carbone organique particulaire des eaux marines du plateau continental aquitain. Journal de Recherche Océanographique VI (2), 37−42]. Radioisotopic measurements (²¹⁰Pb_exc) were made using a semi-planar germanium detector coupled to a multichannel analyser. Radiographical analysis was performed with an X-ray equipment (SCOPIX®) coupled with a radioscopy instrumentation and processing unit.Firstly, a detailed sedimentological map of this shelf has been produced and secondly, geophysical surveys have precisely mapped the geometry of the main mud patch on the continental shelf. In the mud patch itself the rates of sedimentation are between 0.13 and 0.50 cm yr^− 1. The maximum rate of sedimentation is located in the central mud patch, whereas the minimum rate of deposition occurs close to the rocky outcrops. These results seem to be in agreement with the estimation of the total thickness of the mud patch revealed by seismic profiles. The central part corresponds to the maximum thickness of 7 m.Interpretations of the associated oceanic current forcing factors (current direction, wave fetch and wind directional modes) relating to the identified sediment depositional zones are also undertaken.     

Time-series analysis of remote-sensed chlorophyll and environmental factors in the Santa Monica–San Pedro Basin off Southern California

The time-series of remote-sensed surface chlorophyll concentration measured by SeaWiFS radiometer from September 1997 to December 2001 and the relevant hydrological and meteorological factors (remote-sensed sea surface temperature, atmospheric precipitation, air temperature and wind stress) in Santa Monica Bay and adjacent waters off southern California were analyzed using wavelet and cross-correlation statistical methods. All parameters exhibited evident seasonal patterns of variation. Wavelet analysis revealed salient long-term variations most evident in air temperature during El Niño 1997–1998 and in wind stress during La Niña 1998–1999. Short-period (<100 days) variations of remote-sensed chlorophyll biomass were mostly typical to spring seasons. Chlorophyll biomass was significantly correlated with air temperature and wind stress: an increase of chlorophyll biomass followed with 5–6-day time lag an increase of wind stress accompanied by a simultaneous decrease of air temperature. The mechanism of these variations was an intensification of phytoplankton growth resulting from the mixing of water column by wind stress and entrainment of nutrients into the euphotic layer.     

Air–sea gas transfer velocity estimates from the Jason-1 and TOPEX altimeters: Prospects for a long-term global time series

Estimation of global and regional air–sea fluxes of climatically important gases is a key goal of current climate research programs. Gas transfer velocities needed to compute these fluxes can be estimated by combining altimeter-derived mean square slope with an empirical relation between transfer velocity and mean square slope derived from field measurements of gas fluxes and small-scale wave spectra [Frew, N.M., Bock, E.J., Schimpf, U., Hara, T., Hauβecker, H., Edson, J.B., McGillis, W.R., Nelson, R.K., McKenna, S.P., Uz, B.M., Jähne, B., 2004. Air–sea gas transfer: Its dependence on wind stress, small-scale roughness and surface films, J. Geophys. Res., 109, C08S17, doi: 10.1029/2003JC002131.]. We previously reported initial results from a dual-frequency (Ku- and C-band) altimeter algorithm [Glover, D.M., Frew, N.M., McCue, S.J., Bock, E.J., 2002. A Multi-year Time Series of Global Gas Transfer Velocity from the TOPEX Dual Frequency, Normalized Radar Backscatter Algorithm, In: Gas Transfer at Water Surfaces, editors: Donelan, M., Drennan, W., Saltzman, E., and Wanninkhof, R., Geophysical Monograph 127, American Geophysical Union, Washington, DC, 325–331.] for estimating the air–sea gas transfer velocity (k) from the mean square slope of short wind waves (40–100 rad/m) and derived a 6-year time series of global transfer velocities based on TOPEX observations. Since the launch of the follow-on altimeter Jason-1 in December 2001 and commencement of the TOPEX/Jason-1 Tandem Mission, we have extended this time series to 12 years, with improvements to the model parameters used in our algorithm and using the latest corrected data releases. The prospect of deriving multi-year and interdecadal time series of gas transfer velocity from TOPEX, Jason-1 and follow-on altimeter missions depends on precise intercalibration of the normalized backscatter. During the Tandem Mission collinear phase, both satellites followed identical orbits with a mere 73-s time separation. The resulting collocated, near-coincident normalized radar backscatter (σ°) data from both altimeters present a unique opportunity to intercalibrate the two instruments, compare derived fields of transfer velocity and estimate the precision of the algorithm. Initial results suggest that the monthly gas transfer velocity fields generated from the two altimeters are very similar. Comparison of along-track Ku-band and C-band σ° during the collinear phase indicates that observed discrepancies are due primarily to small offsets between TOPEX and Jason-1 σ°. The Jason-1 k values have an apparent bias of + 4% relative to TOPEX, while the precision estimated from the two observation sets is 5–7% and scales with k. The resultant long-term, global, mean k is 16 cm/h.     

Ichthyoplankton in the neritic and coastal zone of Antarctica and Subantarctic islands: A review

Since the article published by Loeb et al. [Loeb, V.J., Kellermann, A., Koubbi, P., North, A.W., White, M., 1993. Antarctic larval fish assemblages: a review. Bull. Mar. Sci. 53(2), 416–449.] about Antarctic ichthyoplankton, many surveys were carried out in different sectors of the Southern Ocean focusing on different aspects of the ecology of fish larvae. Some of these researches were conducted in the Subantarctic Kerguelen Islands and others on the continental shelf off Terre Adélie and Georges V land. Oceanographic and geographic features influence fish larvae ecology such as island mass effects, gyres, canyons. Antarctic fishes show also temporal segregation of spawning which induces temporal succession of early stage larvae. This avoids competition and probably the predation on early stages for species having few recruits. In that case, we have to understand how these larvae can deal with the match–mismatch with their preys and how they find sufficient food to survive. But our knowledge on Antarctic fish larvae is still insufficient as we do not know larvae for quite a lot of species and because of the difficulty to sample during winter.     

Multivariate assimilation in MERCATOR project: New statistical parameters from forecast error estimation

The new operational prototype of Mercator (french Global Ocean Data Assimilation Experiment contribution) is composed of a North Atlantic primitive equation ocean model OPA (Ocean Parallel Algorithm between 20°S and 70°N, [Madec, G., P. Delecluse, M. Imbard and C. Lévy (1998). OPA8.1 ocean general circulation model reference manuel. Notes du pôle de modélisation IPSL. n°11: 91p]) and of a multivariate and multidata assimilation scheme [De Mey, P. and M. Benkiran (2002). “A multivariate reduced-order optimal interpolation method and its application in Mediterranean basin-scale circulation.” Ocean Forecasting : Conceptual basis and application, Pinardi, N., Springer Verlag.] This system has already given some significant improvements from previous Mercator configurations (M. Benkiran, personal communication). However some biases on ocean state still remain in the tropics where the reduced-order optimal interpolation scheme is suspected to be ill-parameted in the model forecast error. Indeed the guess error covariance matrix is decomposed into an error variance value and a spatio-temporal correlation function which are assumed to have some “good” properties (spatial homogeneity of the correlation function, constant ratio between signal and error variance). This study shows how we can use ensemble methods to validate these assumptions. We can see that the correlation function can reach negative values locally, mostly in regions of high variability contradictory with the homogeneous hypothesis. The reduced space used in the operational configuration is based on the signal seasonal Empirical Orthogonal Functions (EOFs). An empirical relationship between signal and error variance has been set and the correlation function is the same on every dimension of the reduced space. By projection of the estimated guess error variance onto the reduced space, we find a repartition of this quantity quite different to what was set in the system. The error statistics is found to be inhomogeneous compared to hypothesis made in the assimilation scheme. These two new parameters tested separately in the assimilation scheme gives significant improvements of the forecast and analysis results. This is particularly obvious in the tropics. But relationship between signal and error statistics (as assumed in the optimal interpolation) is found to be complex.     

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.     

Spatial variability of the barotropic M2 constituent tidal current over the Rías Baixas Galician shelf (NW Spain)

Using current meters and systematically repeated vessel mounted ADCP data legs, we describe the role of the semidiurnal and diurnal constituents in the tidal currents and investigate the mesoscale variability of the M₂ tidal currents field along three cross-shelf legs over the Galician shelf. We found that the estimated values of the semi-major axes of the main semidiurnal tidal constituents are closely related to those obtained in previous research and that the tidal currents are predominantly semidiurnal. Amplitudes of the M₂ and S₂ semi-major axes constituents are comparable, and important fortnightly variability should thus be expected in the amplitude of the semidiurnal tidal currents. Vertical profiles of the semi-major axes of the main tidal constituents reveal that, in the absence of stratification, and far from the bottom frictional effects, the tidal currents are mostly barotropic. Over the slope, the tidal ellipses tilt from a north–south orientation toward the northeast–southwest direction and, as we go further into the shelf, their amplitude is increased, a feature that is reasonably explained by the transition from the dynamics of a Poincaré wave offshore to the dynamics of a northward alongshore propagating Kelvin wave over the shelf. The Ría de Vigo exerts a notable influence over the tidal currents, acting as an extension of the shelf and introducing a meridional variation in the orientation of the ellipses from the alongshore direction in the southern straight coast toward the northeast–southwest direction at the latitude of the Rías. Finally, the spatial fitting methods we tried have proved to be successful in VMADCP data detiding in this region.     

Statistical comparison of trace metal concentrations in wild mussels (Mytilus galloprovincialis) in selected sites of Galicia and Gulf of Biscay (Spain)

Concentrations of five trace metals (Cd, Hg, Pb, Cu and Zn) in wild mussels (Mytilus galloprovincialis) from the Galician and Gulf of Biscay areas in Spain were monitored from 2000 to 2004. A pool of mussel soft tissue was prepared using 50 or more individuals representing the available size range (35–60 mm) present at the sampling points (A Coruña, Pontevedra and Vigo–Galicia, NW of Spain—and Avilés, Bilbao and Santander–Gulf of Biscay). Metals were analysed by atomic absorption spectrometry, AAS, (flame–AAS: Cu and Zn; graphite furnace–AAS: Cd and Pb; flow injection-cold vapour–AAS: total Hg). Quality of the chemical analyses was assessed by participation in periodic QUASIMEME intercalibration exercises carried out from 2000 to 2004, while the present study was developed . Univariate statistical studies, Anova (Analysis of the Variance) and Principal Component Analysis (autoscaled data, Varimax rotation) were carried out. Differences between the two areas as well as certain temporal trends were found and, in general, the Gulf of Biscay samples showed higher metallic contents, mainly close to Avilés, where iron and steel factories thrown residues for years into a nearby coastal area.     

Application of new parameterizations of gas transfer velocity and their impact on regional and global marine CO2 budgets

One of the dominant sources of uncertainty in the calculation of air–sea flux of carbon dioxide on a global scale originates from the various parameterizations of the gas transfer velocity, k, that are in use. Whilst it is undisputed that most of these parameterizations have shortcomings and neglect processes which influence air–sea gas exchange and do not scale with wind speed alone, there is no general agreement about their relative accuracy.The most widely used parameterizations are based on non-linear functions of wind speed and, to a lesser extent, on sea surface temperature and salinity. Processes such as surface film damping and whitecapping are known to have an effect on air–sea exchange. More recently published parameterizations use friction velocity, sea surface roughness, and significant wave height. These new parameters can account to some extent for processes such as film damping and whitecapping and could potentially explain the spread of wind-speed based transfer velocities published in the literature.We combine some of the principles of two recently published k parameterizations [Glover, D.M., Frew, N.M., McCue, S.J. and Bock, E.J., 2002. A multiyear time series of global gas transfer velocity from the TOPEX dual frequency, normalized radar backscatter algorithm. In: Donelan, M.A., Drennan, W.M., Saltzman, E.S., and Wanninkhof, R. (Eds.), Gas Transfer at Water Surfaces, Geophys. Monograph 127. AGU,Washington, DC, 325–331; Woolf, D.K., 2005. Parameterization of gas transfer velocities and sea-state dependent wave breaking. Tellus, 57B: 87–94] to calculate k as the sum of a linear function of total mean square slope of the sea surface and a wave breaking parameter. This separates contributions from direct and bubble-mediated gas transfer as suggested by Woolf [Woolf, D.K., 2005. Parameterization of gas transfer velocities and sea-state dependent wave breaking. Tellus, 57B: 87–94] and allows us to quantify contributions from these two processes independently.We then apply our parameterization to a monthly TOPEX altimeter gridded 1.5° × 1.5° data set and compare our results to transfer velocities calculated using the popular wind-based k parameterizations by Wanninkhof [Wanninkhof, R., 1992. Relationship between wind speed and gas exchange over the ocean. J. Geophys. Res., 97: 7373–7382.] and Wanninkhof and McGillis [Wanninkhof, R. and McGillis, W., 1999. A cubic relationship between air−sea CO2 exchange and wind speed. Geophys. Res. Lett., 26(13): 1889–1892]. We show that despite good agreement of the globally averaged transfer velocities, global and regional fluxes differ by up to 100%. These discrepancies are a result of different spatio-temporal distributions of the processes involved in the parameterizations of k, indicating the importance of wave field parameters and a need for further validation.     

Evolution of the open-sea eddy ALGERS'98 in the Algerian Basin with Lagrangian trajectories and remote sensing observations

The westward evolution of an open-sea anticyclonic eddy along the western Algerian Basin is shown, for the first time, by means of 15 buoy trajectories and remote sensing observations. For 3 months, the buoy trajectories described several anticyclonic loops in periods of 4–21 days. The eddy's movement, translation, and rotation were separated with a kinematic model, resulting in a mean translation speed of 2 km/day, which fits the self-propulsion speed predicted on theoretical models for isolated eddies on a beta plane. Fluctuations in translation speed were associated with advection of the mean flow and topographic interactions. Both mechanisms changed the eddy's horizontal shape from circular to elliptical, inducing fluctuations in its swirl velocity and solid-body rotation. The initial stage of the eddy is an isolated asymmetric dipole, comprised by a small cyclone and a large anticyclone, the latter generated from a frontal instability, which under the Coriolis term acquires anticyclonic relative vorticity. During its first days of life, the anticyclonic eddy was shallow Ro=0.9 and small (diameter less than 50 km). Later on, it reached a diameter of 150 km and a vertical structure of 3 km (Ro=0.1). A retrospective analysis with infrared images shows that the eddy's generation took place at about 3–4°E. Then, the eddy completed a counterclockwise circuit never before reported in other studies and ended up at the entrance of the Algerian Basin, where the interaction with the topography and the coastal instability induced its decay. The eddy's life span was 10 months. Computations of the heating rate following clusters of buoy trajectories show fluctuations throughout the eddy's journey, induced by advection and a seasonal warming.     

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.     

On the distribution of dissolved lead in the Loire estuary and the North Biscay continental shelf, France

The dissolved lead was studied in the whole salinity gradient of the system composed of the Loire estuary and the North Biscay continental shelf. About 130 samples were collected in winter 2001 and spring 2002 during Nutrigas and Gasprod campaigns (Programme PNEC-Golfe de Gascogne, RV Thalassa) and metal measurements were conducted on board by Potentiometric Stripping Analysis. In the Loire estuary, levels of dissolved lead ranged from 0.15 to 0.24 nM and from 0.04 to 0.26 nM in winter and spring, respectively. Compared to the concentrations reported in 1987 and 1990 (0.4–1.7 nM; Boutier, B., Chiffoleau, J.F., Auger, D., Truquet, I., 1993. Influence of the Loire river on dissolved lead and cadmium concentrations in coastal waters of Brittany. Estuar. Coast. Shelf S., 36:133–143, Estuarine, Coastal and Shelf Science 36, 133–143) our study indicated much lower values. The fall in concentration in the estuary could be attributed to the stopping of activity of Octel, a big manufacturer of tetra alkyl lead. Discharge in dissolved metal to the continental shelf by the Loire river was assessed as 7.5 and 1.9 kg day^− 1 for winter and spring, respectively. On the continental shelf, levels of dissolved lead varied within 0.06 and 0.27 nM in winter (0.15 ± 0.06 nM, sd = 1.96, n = 49), whereas concentrations measured in spring were in the range 0.06–0.17 nM (0.09 ± 0.03 nM, sd = 1.96, n = 60). This difference in metal concentration was related to the amounts of rainfall that have fallen over the continental shelf: estimations of inputs by this way (74 and 32 kg day^− 1 in winter and spring, respectively) appeared to be significantly higher than inputs from the Loire river (7.5 and 1.9 kg day^− 1 in winter and spring, respectively). The distributions of dissolved metal in the surface waters highlighted the role of suspended particular matter (SPM) for a rapid “trapping” of lead near the mouth of the estuary. The vertical distributions showed, in the stratified area, a biological transfer of lead between winter and spring from surface waters to the halocline.     

Whitecap coverage in coastal environment for steady and unsteady wave field conditions

Breaking waves represent a “key” parameter for many applications involved with a large number of environmental phenomena. In particular, it is well recognized that the whitecap cover induced by breaking waves allows substantial enhancement of heat, momentum, gas and particle transfer at the air–sea interface. A large number of studies were conducted during the last decades on the variation of the whitecap fraction, commonly noted W. The results presented in this paper deal with the evolution of the whitecap coverage in coastal zone. In such areas, the wave field is often unsteady with an important variety of sea state developments. The present analysis is based on an extensive series of data obtained during an experimental campaign which took place on the Mediterranean coast in 2001. The results allow observation of the influence of the sea state conditions of the wave field on the whitecap coverage. In addition, this paper confirms the occurrence of a peak in the variations of the whitecap fraction with the wave age for coastal areas as suggested by Lafon et al. [Lafon, C., Piazzola, J., Forget, P., Le Calvé, O. and Despiau, S., 2004. Analysis of the variations of the whitecap fraction as measured in a coastal zone. Boundary-Layer Meteorol., 111: 339-360.]. A wave age dependent model for the whitecap fraction is then proposed, which takes into account both the wind and the wave influence, and hence, is characteristics of the different sea state conditions.     

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.     

Cross-shelf circulation and phytoplankton distribution at the summertime New England shelfbreak front

We investigate aspects of the secondary (cross-shelf) circulation at the Middle Atlantic Bight shelfbreak front using high-resolution data collected on the New England Shelf in August 2002. The alongshore shelfbreak jet coincides with the front at the seaward edge of the cold pool (remnant winter shelf water) and there is a suggestion of a cross-stream convergence centered at the jet core. Despite indications of convergence we found no evidence of a surface subduction on the seaward side of the front. At depth 70 m near the shelfbreak there was a patch of chlorophyll, located within a local temperature–salinity maximum which, though significantly below the euphotic zone, appeared to be photo-acclimated and viable. The chlorophyll feature could be the result of a local subduction by a larger scale eddy circulation seaward of the front.Dye tracer experiments directly observed the convergence at the foot of the shelfbreak front and subsequent upwelling of bottom boundary layer water along the shoreward side of the shelfbreak front. But, we found no evidence that this upwelling influenced productivity at the front. Further, since there was no cross-shelf maximum in subsurface chlorophyll at the front, we conclude that this productivity is in general, sustained by a broadly distributed local vertical nutrient flux from an underlying nutrient reservoir.     

Interactions of wind and density driven currents in North Sea ROFIs—a model study

Three different versions of a baroclinic three-dimensional circulation model of the North Sea are used to obtain information on the wind and density interactions in the North Sea ROFIs (Regions Of Freshwater Influence): the standard version with fully prognostic treatment of salinity and temperature is compared to a barotropic model run on the same grid on the one hand and to an also fully prognostic model run on a four times coarser grid on the other hand. In order to gain knowledge on the wind and density interactions, two opposing wind directions are chosen for investigation, namely a time of strong north wind, 21st–28th April 1982, and a time of strong southwest wind, 22nd–24th May 1982. In the April case the effect of the salinity gradients on the border of the ROFIs of Rhine, Weser, Ems and Elbe, i.e. along the continental shore, is shown to lead to a clear enhancement of the mean surface currents. In May this result is partly disguised by the additional effect of the thermocline in the deeper parts of the North Sea, i.e. in the classical shelf sea regime region. Nevertheless, the same pattern of enhanced mean surface currents along the coast is detected and is of the same order of magnitude as in the April case. It is thus concluded that although the circulation in the North Sea is reversed by the wind, the density induced component of the general circulation is modified only slightly.