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.     

Late summer stratification, internal waves, and turbulence in the Yellow Sea

Microstructure profiling measurements at two locations in the Yellow Sea (a deeper central basin and a local shelf break) were analyzed focusing on tidal and internal-wave induced turbulence near the bottom and in the pycnocline. A classical three-layer density structure consisting of weakly stratified surface and bottom boundary layers and a narrow sharp pycnocline is developed by the end of warm season. Turbulence in the surface layer was not influenced by the tidal forcing but by the diurnal cycle of buoyancy flux and wind forcing at the sea surface. The enhanced dissipation and diffusivity generated by the shear stress at the seafloor was found in the water interior at heights 10–15 m above the bottom with a phase shift of ~ 5–6 m/h. No internal waves, turbulence, or mixing were detected in the pycnocline in the central basin, in contrast to the pycnocline near the local shelf break wherein internal waves of various frequencies were observed all the time. The thickness of the surface layer near the local shelf break slightly exceeded that of the bottom layer (20 vs. 18 m). A 5–6 m high vertical displacement of the pycnocline, which emerged during the low tide, was arguably caused by the passage of an internal soliton of elevation. During this episode, the gradient Richardson number decreased below 0.25 due to enhanced vertical shear, leading to local generation of turbulence with dissipation rates exceeding the background level by an order of magnitude.     

Wind-induced mesoscale circulation off the Ebro delta, NW Mediterranean: a numerical study

It is well known that the general circulation on the Catalan continental slope is dominated by a quasi-permanent southwestward geostrophic jet associated to the so-called Catalan front [Millot, C., 1987. Circulation in the western Mediterranean sea. Oceanol Acta 10, 143–149; Font, J., Salat, J., Tintoré, J., 1988. Permanent features of the circulation in the Catalan Sea. Oceanol. Acta 9, 51–57]. On the continental shelf, however, the flow is modified by the action of friction which enhances also other nonlinear interactions. Several authors have hypothesized that the shelf circulation is anticyclonic north of the Ebro delta [Salat, J., Manriquez, M., Cruzado, A., 1978. Hidrografia del golfo de Sant Jordi. Campaña Delta (Abril 1970). Investigación Pesquera 42 (2), 255–272; Ballester, A., Castellví, J., 1980. Estudio hidrográfico y biológico de las plataformas continentales españolas: I. Efecto de los efluentes de una planta de energía nuclear en el Golfo de San Jorge (Febrero 1975–Octubre 1976). Informes Técnicos del Instituto de Investigaciones Pesqueras 76, 70 pp.]. A quasi-3D finite element code based on the shallow-water equations has been used to explore the effect of several mechanisms which might be responsible for such a local circulation pattern, and in particular of wind. The obtained numerical results suggest that the basic anticyclonic structure of the mean flow is controlled by the bathymetry and that the clockwise-rotating mean flow pattern is not a permanent circulation feature. It is seen that the characteristic local wind stress fields—computed through interpolation of the records of a local network of meteo stations—may ‘enhance' or ‘delete' the anticyclonic gyre depending on the sign of their relative vorticity. According to the analysis of a 2-yr record of local wind data, the net contribution of wind events with a duration longer than 24 h is to reinforce the anticyclonic circulation (over 70% of these wind fields supply negative relative vorticity to the study area).     

A numerical model of the long term flow along the Malin-Hebrides shelf

A three dimensional hydrodynamic model of the Malin-Hebrides shelf region is used to investigate the spatial variability of the wind and tidally induced residual flow in the region and the influence of flow from the Irish Sea and along the shelf edge. By this means it is possible to understand the spatial variability in the long term observed flow fields in the region and the range of driving forces producing this flow. The model uses a sigma coordinate grid in the vertical with a finer grid in the near surface and near bed shear layers. The vertical diffusion of momentum in the model is parameterised using an eddy viscosity coefficient which is derived from turbulence energy closure models. Two different turbulence models are used to compute the eddy viscosity, namely a two-equation (itq²−q²ℓ) model which has prognostic equations for both turbulence energy and mixing length and a simpler model in which the mixing length is a specified algebraic function of the water depth.The wind induced response to spatially and temporally constant orthogonal wind stresses, namely westerly and southerly winds of 1 N m⁻², are derived from the model. By using orthogonal winds and assuming linearity, then to first order the response to any wind direction can be derived. Computed flows show a uniform wind driven surface layer of magnitude about 3% of the wind speed and direction 15 ° to the right of the wind, in deep water. Currents at depth particularly in the shelf edge and near coastal region show significant spatial variability which is related to variations in bottom topography and the coastline.Calculations show that tidal residual flows are only significant in the near coastal regions where the tidal current is strong and exhibits spatial variability. Flow into the region from the Irish Sea through the North Channel although having its greatest influence in the near coastal region, does affect currents near the shelf edge region. Again the spatial variability of the flow is influenced by topographic effects.A detailed examination of wind induced current profiles together with turbulence, mixing length and viscosity, at a number of locations in the model from deep ocean to shallow near coastal, shows that both turbulence models yield comparable results, with the mixing length in the two equation model showing a similar dependence to that specified in the simpler turbulence model.Calculations clearly show that flow along the shelf edge area to the west of Ireland and from the Irish Sea entering the region, together with local wind forcing can have a major effect upon currents along the Malin-Hebrides shelf. The flow fields show significant spatial variability in the region, comparable to those deduced from long term tracer measurements. The spatial variability found in the calculations suggests that a very intense measurement programme together with inflow measurements into the area is required to understand the circulation in the region, and provide data sets suitable for a rigorous model validation.     

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.     

Large eddy simulations for quasi-2D turbulence in shallow flows: A comparison between different subgrid scale models

In this study, the performance of the horizontal large eddy simulation module, developed at the University of Leuven (HLES-KULeuven module) is assessed. A comparison between different subgrid scale models has been carried out. The study is concerned with the non-rotating and unstratified flows. The results of the simulation for an oscillatory backward facing (BFS) flow are presented in case of an expanding flume based on a one-length scale approach and a two-length scale approach. Three subgrid scale (SGS) models have been tested: Smagorinsky SGS model (Smagorinsky, J., (1963). General circulation experiments with the primitive equations, I. the basic experiments. Monthly Weather Review, 91(3), 99–164), Uittenbogaard SGS model (Uittenbogaard, R.E., and van Vossen, B., (2004). Subgrid-scale model for quasi-2D turbulence in shallow water. Shallow Flows. Jirka and Uijttewaal (Eds.), Taylor & Francis Group, London, ISBN 90 5809 700 5) and a proposed two-length scale approach. The first two models are considered to be a one-length scale models. A simulation without a subgrid scale model for the horizontal mixing has also been conducted. In all simulations, a quadratic friction model parameterizes the dissipation produced by the 3D-subdepth scale turbulence. The two-length scale concept uses a newly mixing length formulation for the quasi-2D turbulence and doesn't depend on the filter width in contrast to the one-length scale approach, in which the mixing length is function of the filter width. The outputs of the HLES-KULeuven module have been compared with the experimental data taken from Stelling, G.S., and Wang, L.X., (1984). Experiments and computations on separating flow in an expanding flume. Dept. Civil Engineering, Delft University of Technology, Report 2–84.). The two-length scale approach has been validated with experimental data from SERC Flood Channel Facility at HR Wallingford. In general, there is a qualitative agreement with the experimental data. It has also been found that the two-length scale approach produces more elongated and less isotropic vortex than the one-length scale models.     

Primary productivity of the Palmer Long Term Ecological Research Area and the Southern Ocean

A major objective of the Palmer Long Term Ecological Research (Palmer LTER) project is to obtain a comprehensive understanding of the various components of the Antarctic marine ecosystem. Phytoplankton production plays a key role in this so-called high nutrient, low chlorophyll environment, and factors that regulate production include those that control cell growth (light, temperature, and nutrients) and those that control cell accumulation rate and hence population growth (water column stability, grazing, and sinking). Sea ice mediates several of these factors and frequently conditions the water column for a spring bloom which is characterized by a pulse of production restricted in both time and space. This study models the spatial and temporal variability of primary production within the Palmer LTER area west of the Antarctic Peninsula and discusses this production in the context of historical data for the Southern Ocean. Primary production for the Southern Ocean and the Palmer LTER area have been computed using both light-pigment production models [Smith, R.C., Bidigare, R.R., Prézelin, B.B., Baker, K.S., Brooks, J.M., 1987. Optical characterization of primary productivity across a coastal front. Mar. Biol. (96), 575–591; Bidigare, R.R., Smith, R.C., Baker, K.S., Marra, J., 1987. Oceanic primary production estimates from measurements of spectral irradiance and pigment concentrations. Global Biogeochem. Cycles (1), 171–186; Morel, A., Berthon, J.F., 1989. Surface pigments, algal biomass profiles and potential production of the euphotic layer—relationships reinvestigated in view of remote-sensing applications. Limnol. Oceanogr. (34), 1545–1562] and an ice edge production model [Nelson, D.M., Smith, W.O., 1986. Phytoplankton bloom dynamics of the western Ross Sea ice edge: II. Mesoscale cycling of nitrogen and silicon. Deep-Sea Res. (33), 1389–1412; Wilson, D.L., Smith, W.O., Nelson, D.M., 1986. Phytoplankton bloom dynamics of the Western Ross Sea ice edge: I. primary productivity and species-specific production. Deep-Sea Res., 33, 1375–1387; Smith, W.O., Nelson, D.M., 1986. Importance of ice edge phytoplankton production in the Southern Ocean. BioScience (36), 251–257]. Chlorophyll concentrations, total photosynthetically available radiation (PAR) and sea ice concentrations were derived from satellite data. These same parameters, in addition to hydrodynamic conditions, have also been determined from shipboard and Palmer Station observations during the LTER program. Model results are compared, sensitivity studies evaluated, and productivity of the Palmer LTER region is discussed in terms of its space time distribution, seasonal and interannual variability, and overall contribution to the marine ecology of the Southern Ocean.     

Benthic foraminifera in the surface sediments of the Beaufort Shelf and slope, Beaufort Sea, Canada: Applications and implications for past sea-ice conditions

This paper presents new data on distribution patterns of modern benthic foraminifera and other microfossils from the Canadian Arctic, specifically the Beaufort Shelf and slope. The material was collected in June to August of 2004 and is the first of its kind in this area to be collected since 1970. We examined the smaller sizes (45–63µm) as well as > 63µm and discovered that many species had been severely underrepresented in previous studies. Deep sea forms, that had been overlooked previously, were common on the shelf; two species (Elphidiella arctica and Ammotium cassis) appeared in preliminary results to be indicators of methane seepage; and it was possible to make determinations of sea-ice coverage using a combination of foraminifera and tintinnids (planktic ciliates). Our data indicated the presence of many of the same species as previous studies from this area, but improved techniques of sample processing greatly increased the number of specimens and species found (particularly the small deep sea arctic species Buliminella hensoni and Bolivina arctica) which provide much more reliable data for paleoceanographic determinations. One of the primary objectives for this work was to provide baseline data to help determine paleo-ice cover; these data cover a broad range of conditions on the Beaufort Shelf that make it possible to achieve this objective as well as improving what it is known about the assemblages on this shelf as compared to other arctic shelf areas, such as the Siberian Shelf).     

Nutrient and organic matter patterns across the Mackenzie River, estuary and shelf during the seasonal recession of sea-ice

Suspended material, nutrients and organic matter in Mackenzie River water were tracked along a 300 km transect from Inuvik (Northwest Territories, Canada), across the estuarine salinity gradient in Kugmallit Bay, to offshore marine stations on the adjacent Mackenzie Shelf. All particulates measured (SPM, POC, PN, PP) declined by 87–95% across the salinity gradient and levels were generally below conservative mixing. Organic carbon content of suspended material decreased from 3.1% in the river to 1.7% in shelf surface waters while particulate C:N concurrently decreased from 17.1 to 8.6. Nitrate and silicate concentrations declined by more than 90% across the salinity gradient, with nitrate concentrations often below the conservative mixing line. Phosphate concentrations increased from 0.03 μmol/L in the river to 0.27 μmol/L over shelf waters, thereby changing the inorganic nutrient regime downstream from P to N limitation. Dissolved organic carbon decreased conservatively offshore while dissolved organic N and P persisted at high levels in the Mackenzie plume relative to river water, increasing 2.7 and 25.3 times respectively. A deep chlorophyll-a maximum was observed at two offshore stations and showed increases in most nutrients, particulates and organic matter relative to the rest of the water column. During river passage through the Mackenzie estuary, particulate matter, dissolved organic carbon and inorganic nutrients showed sedimentation, dilution and biological uptake patterns common to other arctic and non-arctic estuaries. Alternatively, inorganic content of particles increased offshore and dissolved organic N and P increased substantially over the shelf, reaching concentrations among the highest reported for the Arctic Ocean. These observations are consistent with the presence of a remnant ice-constrained (‘stamukhi’) lake from the freshet period and a slow flushing river plume constrained by sea-ice in close proximity to shore. Nutrient limitation in surface shelf waters during the ARDEX cruise contributed to the striking deep chlorophyll-a maximum at 21 m where phytoplankton communities congregated at the margin of nutrient-rich deep ocean waters.     

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.     

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.     

Turbulence in coastal fronts near the mouths of Block Island and Long Island Sounds

Measurements of turbulence were performed in four frontal locations near the mouths of Block Island Sound (BIS) and Long Island Sound (LIS). These measurements extend from the offshore front associated with BIS and Mid-Atlantic Bight Shelf water, to the onshore fronts near the Montauk Point (MK) headland, and the Connecticut River plume front. The latter feature is closely associated with the major fresh water input to LIS. Turbulent kinetic energy (TKE) dissipation rate, ε, was obtained using shear probes mounted on an autonomous underwater vehicle. Offshore, the BIS estuarine outflow front showed, during spring season and ebb tide, maximum TKE dissipation rate, ε, estimates of order 10^− 5 W/kg, with background values of order 10^− 6 to 10^− 9 W/kg. Edwards et al. [Edwards, C.A., Fake, T.A., and Bogden, P.S., 2004a. Spring–summer frontogenesis at the mouth of Block Island Sound: 1. A numerical investigation into tidal and buoyancy-forced motion. Journal of Geophysical Research 109 (C12021), doi:10.1029/2003JC002132.] model this front as the boundary of a tidally driven, baroclinically adjusted BIS flow around the MK headland eddy. At the entrance to BIS, near MK, two additional fronts are observed, one of which was over sand waves. For the headland site front east of MK, without sand waves, during ebb tide, ε estimates of 10^− 5 to 10^− 6 W/kg were observed. The model shows that this front is at the northern end of an anti-cyclonic headland eddy, and within a region of strong tidal mixing. For the headland site front further northeast over sand waves, maximum ε estimates were of order 10^− 4 W/kg within a background of order 10^− 7–10^− 6 W/kg. From the model, this front is at the northeastern edge of the anti-cyclonic headland eddy and within the tidal mixing zone. For the Connecticut River plume front, a surface trapped plume, during ebb tide, maximum ε estimates of 10^− 5 W/kg were obtained, within a background of 10^− 6 to 10^− 8 W/kg. Of all four fronts, the river plume front has the largest finescale mean-square shear, S² ~ 0.15 s^− 2. All of the frontal locations had local values of the buoyancy Reynolds number indicating strong isotropic turbulence at the dissipation scales. Local values of the Froude number indicated shear instability in all of the fronts.     

On the parameterization of interleaving and turbulent mixing using CTD data from the Azores Frontal Zone

CTD-data obtained in the Azores Frontal Zone using a towed undulating vehicle are analyzed to study the relationship between characteristics of intrusions and mean parameters of the thermohaline field. A self-similar dependence between intrusion intensity and hydrological parameters is obtained. The most well-founded interpretation of the empirical dependence is as follows: (a) the main source supporting intrusive layering is the salt finger convection; (b) the abrupt decrease of intrusion intensity with the reduction of geostrophic Richardson number obtained from the analysis is explained by the beginning of turbulence when salt fingers do not work any longer, so the “driving force” for intrusive motion disappears. These results are consistent with the conclusions of the paper [Kuzmina N.P., Rodionov V.B., 1992. About the influence of baroclinicity upon generation of the thermohaline intrusions in the oceanic frontal zones. Izvestiya Akad. Nauk SSSR, Atmosperic and Oceanic Physics 28 (10–11), 1077–1086]. These conclusions imply that there are three main mechanisms of intrusive layering at oceanic fronts, namely the 2D baroclinic instability of geostrophic flow, the vertical shear instability and the thermohaline instability where the driving source of intrusive motion is double diffusive convection. The baroclinic and thermohaline instabilities can generate intrusions of large vertical scale, while vertical shear instability usually gives rise to thin turbulent layers. Turbulence in these thin layers can prevent salt finger convection and thus destroy the energy source of the intrusive motion conditioned by thermoclinicity. Therefore, the baroclinicity plays two parts in the processes of the intrusive layering: (1) it prevents double-diffusion interleaving by means of turbulence, and (2) it generates intrusions due to the 2D baroclinic instability of geostrophic current. Using features of thermohaline interleaving as a specific tracer of turbulent mixing, we have estimated turbulent mixing coefficient as k_tRi^−0.8 (Ri>1), where Ri is the geostrophic Richardson number. Application of the proposed approach to other frontal zones is discussed.     

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.     

Sea ice and the onshore–offshore gradient in pre-winter zooplankton assemblages in southeastern Beaufort Sea

Zooplankton communities were studied in southeastern Beaufort Sea (Arctic Ocean) in September–October 2002. Cluster analysis and non-metric multidimensional scaling revealed three distinct mesozooplankton assemblages. A neritic assemblage occurred on the Mackenzie Shelf and in Franklin Bay, while distinct off-shelf assemblages prevailed in the Cape Bathurst Polynya and on the Beaufort Slope respectively. Over 95% of the mesozooplankton was comprised of eight copepod taxa. Pseudocalanus spp. contributed predominantly to the discrimination of the three assemblages and was the only significant indicator of the Shelf assemblage. Oithona similis, Oncaea borealis, Metridia longa and Calanus hyperboreus were indicators of the Polynya assemblage. Cyclopina sp. and Microcalanus pygmaeus were indicative of the overall off-shelf community (Polynya and Slope assemblages). The importance of omnivores and carnivores increased from the shelf to the polynya and the slope. Station depth and duration of reduced ice conditions during summer (< 50% ice concentration) underpinned the distribution of the assemblages (r² = 0.71 and 0.45 respectively). The abundance of Pseudocalanus spp. was independent of depth and increased with the duration of reduced ice conditions (r_s = 0.438). The abundance of Cyclopina sp., M. pygmaeus and other indicators of the offshore assemblages followed the opposite trend (r_s = − 0.467 and − 0.5 respectively). Under continued climate warming, a reduction of the ice cover will affect the biogeography of mesozooplankton on and around the Mackenzie Shelf, to the potential advantage of Pseudocalanus spp. and other calanoid herbivores.     

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.     

Comparison between a reanalyzed product by 3-dimensional variational assimilation technique and observations in the Ulleung Basin of the East/Japan Sea

Reanalyzed products from a MOM3-based East Sea Regional Ocean Model with a 3-dimentional variational data assimilation module (DA-ESROM), have been compared with the observed hydrographic and current datasets in the Ulleung Basin (UB) of the East/Japan Sea (EJS). Satellite-borne sea surface temperature and sea surface height data, and in-situ temperature profiles have been assimilated into the DA-ESROM. The performance of the DA-ESROM appears to be efficient enough to be used in an operational ocean forecast system.Comparing with the results from Mitchell et al. [Mitchell, D. A., Watts, D. R., Wimbush, M., Teague, W.J., Tracey, K. L., Book, J. W., Chang, K.-I., Suk, M.-S., Yoon, J.-H., 2005a. Upper circulation patterns in the Ulleung Basin. Deep-Sea Res. II, 52, 1617-1638.], the DA-ESROM fairly well simulates the high variability of the Ulleung Warm Eddy and Dok Cold Eddy as well as the branching of the Tsushima Warm Current in the UB. The overall root-mean-square error between 100 m temperature field reproduced by the DA-ESROM and the observed 100-dbar temperature field is 2.1 °C, and the spatially averaged grid-to-grid correlation between the two temperature fields is high with a mean value of 0.79 for the inter-comparison period.The DA-ESROM reproduces the development of strong southward North Korean Cold Current (NKCC) in summer consistent with the observational results, which is thought to be an improvement of the previous numerical models in the EJS. The reanalyzed products show that the NKCC is about 35 km wide, and flows southward along the Korean coast from spring to summer with maximum monthly mean volume transport of about 0.8 Sv in August–September.     

Phytoplankton patch patterns: Seascape anatomy in a turbulent ocean

Marine phytoplankton experience competition, predation, infection and aggregation occurring across distances of micrometres to centimetres. However, the consequences of these interactions influence global processes, such as climate and fisheries productivity. There is a long-standing default assumption that these global processes cannot be traced to plankton distributions and interactions below a few metres because of the homogenising effect of turbulence [Hutchinson, G.E., 1961. The paradox of the plankton. Am. Nat. 95, 137–146.; Siegel, D.A., 1998. Resource competition in a discrete environment: Why are plankton distributions paradoxical? Limnol. Oceanogr. 43, 1133–1146.]. We show that, in active turbulence, phytoplankton patches, on the order of 10 cm, have repeatable asymmetry and regular spacing over distances of centimetres to tens of metres. The regularity and hierarchical nature of the patches in mixed ocean water means that phytoplankton are distributed in a dynamic, but definite seascape topography, where groups of patches coalesce between intermittent turbulent eddies. These patches may link large scale processes and microscale interactions, acting as fundamental components of marine ecosystems that influence grazing efficiency, taxonomic diversity, and the initiation of aggregation and subsequent carbon flux.     

Suspended sediment and erosion dynamics in Kugmallit Bay and Beaufort Sea during ice-free conditions

The Mackenzie River is the largest river on the North American side of the Arctic and its huge freshwater and sediment load impacts the Canadian Beaufort Shelf. Huge quantities of sediment and associated organic carbon are transported in the Mackenzie plume into the interior of the Arctic Ocean mainly during the freshet (May to September). Changing climate scenarios portend increased coastal erosion and resuspension that lead to altered river-shelf-slope particle budgets. We measured sedimentation rates, suspended particulate matter (SPM), particle size and settling rates during ice-free conditions in Kugmallit Bay (3–5 m depth). Additionally, measurements of erosion rate, critical shear stress, particle size distribution and resuspension threshold of bottom sediments were examined at four regionally contrasting sites (33–523 m depth) on the Canadian Beaufort Shelf using a new method for assessing sediment erosion. Wind induced resuspension was evidenced by a strong relationship between SPM and wind speed in Kugmallit Bay. Deployment of sediment traps showed decreasing sedimentation rates at sites along an inshore–offshore transect ranging from 5400 to 3700 g m^− 2 day^− 1. Particle settling rates and size distributions measured using a Perspex settling chamber showed strong relationships between equivalent spherical diameter (ESD) and particle settling rates (r² = 0.91). Mean settling rates were 0.72 cm s^− 1 with corresponding ESD values of 0.9 mm. Undisturbed sediment cores were exposed to shear stress in an attempt to compare differences in sediment stability across the shelf during September to October 2003. Shear was generated by vertically oscillating a perforated disc at controlled frequencies corresponding to calibrated shear velocity using a piston grid erosion device. Critical (Type I) erosion thresholds (u) varied between 1.1 and 1.3 cm s^− 1 with no obvious differences in location. Sediments at the deepest site Amundsen Gulf displayed the highest erosion rates (22–54 g m^− 2 min^− 1) with resuspended particle sizes ranging from 100 to 930 µm for all sites. There was no indication of biotic influence on sediment stability, although our cores did not display a fluff layer of unconsolidated sediment. Concurrent studies in the delta and shelf region suggest the importance of a nepheloid layer which transports suspended particles to the slope. Continuous cycles of resuspension, deposition, and horizontal advection may intensify with reduction of sea ice in this region. Our measurements coupled with studies of circulation and cross-shelf exchange allow parameterization and modeling of particle dynamics and carbon fluxes under various climate change scenarios.