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.     

Turbulent mixing and internal tides in Gaoping (Kaoping) Submarine Canyon, Taiwan

Turbulent overturning on scales greater than 10 m is observed near the bottom and in mid-depth layers within the Gaoping (formerly spelled Kaoping) Submarine Canyon (KPSC) in southern Taiwan. Bursts of strong turbulence coexist with bursts of strong sediment concentrations in mid-depth layers. The turbulence kinetic energy dissipation rate in some turbulence bursts exceeds 10^− 4 W kg^− 1, and the eddy diffusivity exceeds 10^− 1 m² s^− 1. Within the canyon, the depth averaged turbulence kinetic energy dissipation rate is ~ 7 × 10^− 6 W kg^− 1, and the depth averaged eddy diffusivity is ~ 10^− 2 m² s^− 1. These are more than two orders of magnitude greater than typical values in the open ocean, and are much larger than those found in the Monterey Canyon where the strong turbulent mixing has also been. The interaction of tidal currents with the complex topography in Gaoping Submarine Canyon is presumably responsible for the observed turbulent overturning via shear instability and the breaking of internal tides and internal waves at critical frequencies. Strong 1st-mode internal tides exist in KPSC. The depth averaged internal tidal energy near the canyon mouth is ~ 0.17 m² s^− 2. The depth integrated internal tidal energy flux at the mouth of the canyon is ~ 14 kW m^− 1, propagating along the axis of the canyon toward the canyon head. The internal tidal energy flux in the canyon is 3–7 times greater than that found in Monterey Canyon, presumably due to the more than 10 times larger barotropic tide in the canyon. Simple energy budget calculations conclude that internal tides alone may provide energy sufficient to explain the turbulent mixing estimated within the canyon. Further experiments are needed in order to quantify the seasonal and geographical distributions of internal tides in Gaoping Submarine Canyon and their effects on the sediment flux in the canyon.     

The structure of dissipation in the western Irish Sea front

We report on an intensive campaign in the summer of 2006 to observe turbulent energy dissipation in the vicinity of a tidal mixing front which separates well mixed and seasonally stratified regimes in the western Irish Sea. The rate of turbulent dissipation ε was observed on a section across the front by a combination of vertical profiles with the FLY dissipation profiler and horizontal profiles by shear sensors mounted on an AUV (Autosub). Mean flow conditions and stratification were obtained from a bed mounted ADCP and a vertical chain of thermistors on a mooring. During an Autosub mission of 60 h, the vehicle, moving at a speed of ~ 1.2 m s^− 1, completed 10 useable frontal crossings between end points which were allowed to move with the mean flow. The results were combined with parallel measurements of the vertical profile of ε which were made using FLY for periods of up to 13 h at positions along the Autosub track. The two data sets, which show a satisfactory degree of consistency, were combined to elucidate the space–time variation of dissipation in the frontal zone. Using harmonic analysis, the spatial structure of dissipation was separated from the strong time dependent signal at the M₄ tidal frequency to yield a picture of the cross-frontal distribution of energy dissipation. A complementary picture of the frontal velocity field was obtained from a moored ADCP and estimates of the mean velocity derived from the thermal wind using the observed density distribution. which indicated the presence of a strong (0.2 m s^− 1) jet-like flow in the high gradient region of the front. Under neap tidal conditions, mean dissipation varied across the section by 3 orders of magnitude exceeding 10^− 2 W m^− 3 near the seabed in the mixed regime and decreasing to 10^− 5 W m^− 3. in the strongly stratified interior regime. The spatial pattern of dissipation is consistent in general form with the predictions of models of tidal mixing and does not reflect any strong influence by the frontal jet.     

Estimation of TKE dissipation rates in dense bottom plumes using a Pulse Coherent Acoustic Doppler Profiler (PC-ADP) — Structure function approach

During a hydrographic survey in January 2006 the spreading of inflowing saline water was observed in the Arkona Basin (Western Baltic Sea). Two bottom mounted ‘pulse coherent’ acoustic Doppler profilers (PC-ADP) were used to measure the near-bottom current field of the dense plume with a high temporal (1 s) and spatial resolution (5 cm). In order to estimate the dissipation rate of turbulent kinetic energy () a structure function approach was applied to the beam velocity data. Simultaneous measurements with a microstructure shear profiler (MSS) and an acoustic Doppler velocimeter (ADV) supplied independent data for the verification of the structure function method. Additional measurements with standard CTD, near-bottom towed and vessel mounted acoustic Doppler current profilers (ADCP) completed the data set.The estimated dissipation rates from the structure function approach fit well with the values derived from the ADV and the MSS probe. It is shown that the structure function approach is a reliable and easily applicable method to derive estimates of TKE dissipation rates from PC-ADP beam velocities. The observed dissipation rates ranged between 5 · 10^− 6 and 1 · 10^− 8 W kg^− 1 depending on the hydrographic conditions. Inside the plume the dissipation rates exceeded that of the overlaying brackish water by two orders of magnitude. Since the noise level of velocity data in pulse coherent mode is considerably lower than in the Doppler mode the PC-ADP can also be used for estimates in marine environments with low turbulence level. Reynolds stresses estimated from the PC-ADP and the ADV agreed well at the same depth level. TKE production derived from PC-ADP measurements compared reasonably well with the dissipation rate of TKE in a varying environment.     

Experiments on the survival of six brackish macro-invertebrates from the Baltic Sea after dredged spoil coverage and its implications for the field

Physical disturbance by disposal of dredged materials in estuarine and coastal waters may result in burial of benthic fauna. Survival rates depend on a variety of factors including the type and amount of disposed materials and the lifestyle of the organisms. Laboratory burial experiments using six common macrobenthic invertebrates from a brackish habitat of the western Baltic Sea were performed to test the organisms' escape reaction to dredged material disposal. Experimental lab-results were then extrapolated to a field situation with corresponding bottom topography and covering layer thicknesses at experimental field disposal study sites. Resulted survival rates were then verified by comparison with results of an earlier field study at the same disposal sites.Our experimental design in the lab included the disposal of two types of dredged material (i.e. ‘till’ and ‘sand/till mixture’) and two covering layer depths (i.e. 10–20 cm and 14–40 cm). All three bivalves Arctica islandica (Linnaeus), Macoma balthica (Linnaeus), Mya arenaria (Linnaeus) and the polychaete Nephtys hombergii (Savigny) successfully burrowed to the surface of a 32–41 cm deposited sediment layer of till or sand/till mixture and restored contact with the overlying water. These high escape potentials could partly be explained by the heterogeneous texture of the till and sand/till mixture with ‘voids’. The polychaete Bylgides (Harmothoe) sarsi (Malmgren) successfully burrowed through a 16 cm covering layer whereas the polychaete Lagis koreni (Malmgren) showed almost no escaping reaction. No general differences in escape behaviour after burial were detected between our test species from the brackish habitat and those reported in the literature for the same species in marine environments. However, a size-dependence in mobility of motile polychaetes and M. arenaria was apparent within our study. In comparison to a thick coverage, thin covering layers (i.e. 15–16 cm and 20 cm) increased the chance of the organisms (N. hombergii and M. arenaria) to reach the sediment surface after burial. This was not observed for the other test species. While crawling upward to the new sediment surfaces burrowing velocities of up to 8 cm d^− 1 were observed for the bivalves and up to 20 cm d^− 1 for N. hombergii. Between 17 and 79% of the test organisms showed burrowing activity after experimental burial. The survival rate (defined as the ability to regained contact with the sediment surface) ranged from 0 to 33%, depending on species and on burial depth. The organisms reached the sediment surface by burrowing (polychaetes and bivalves) and/or by extending their siphons to the new sediment surface (bivalves). The extrapolation of laboratory survival rates to the two disposal sites was obtained based on the in situ thicknesses of the dredged spoil layers measured by multi-beam echo sounder. This resulted in total average survival rate estimates for the test species of 45 and 43% for the two disposal sites. The results obtained during the laboratory tests and the following extrapolation to the field were verified by the range of results from a previous field study, using grab sampling shortly before and after a disposal event in June 2001. The effect of dredged material disposal on the tested Baltic Sea benthic macrofauna was assessed by extrapolating the verified laboratory results to the field.     

Using a random displacement model to simulate turbulent particle motion in a baroclinic frontal zone: A new implementation scheme and model performance tests

Numerical artifacts can limit accurate simulation of turbulent particle motion when Lagrangian particle-tracking models are implemented in hydrodynamic models with stratified conditions like fronts. Yet, modeling of individual particle motion in frontal regions is critical for understanding sediment dynamics as well as the transport and retention of planktonic organisms. The objective of this research was to develop a numerical technique to accurately simulate turbulent particle motions in a particle-tracking model embedded within a hydrodynamic model of a frontal zone. A new interpolation scheme, the ‘water column profile’ scheme, was developed and used to implement a random displacement model for turbulent particle motions. A new interpolation scheme was necessary because linear interpolation schemes caused artificial aggregation of particles where abrupt changes in vertical diffusivity occurred. The new ‘water column profile’ scheme was used to fit a continuous function (a tension spline) to a smoothed profile of vertical diffusivities at the x–y particle location. The new implementation scheme was checked for artifacts and compared with a standard random walk model using (1) Well Mixed Condition tests, and (2) dye-release experiments. The Well Mixed Condition tests confirmed that the use of the ‘water column profile’ interpolation scheme for implementing the random displacement model significantly reduced numerical artifacts. In dye-release experiments, high concentrations of Eulerian tracer and Lagrangian particles were released at the same location up-estuary of the salt front and tracked for 4 days. After small differences in initial dispersal rates, tracer and particle distributions remained highly correlated (r = 0.84 to 0.99) when a random displacement model was implemented in the particle-tracking model. In contrast, correlation coefficients were substantially lower (r = 0.07 to 0.58) when a random walk model was implemented. In general, model performance tests indicated that the ‘water column interpolation’ scheme was an effective technique for implementing a random displacement model within a hydrodynamic model, and both could be used to accurately simulate diffusion in a highly baroclinic frontal region. The new implementation scheme has the potential to be a useful tool for investigating the influence of hydrodynamic variability on the transport of sediment particles and planktonic organisms in frontal zones.     

Analysis of the hydrophysical structure of the Sea of Azov in the period of the bottom anoxia development

The hydrophysical and hydrochemical structure of the Sea of Azov, with developed bottom anoxia, was studied during the RV “Akvanavt” cruise from July 31 to August 03, 2001. The anoxic zone with a thickness from 0.5 to 4 m above the bottom was found in all deep regions of the Sea. Concentrations of hydrochemical parameters were similar to the pronounced anoxic conditions (about 90 mmol m^− 3 of hydrogen sulfide, 17 mmol m^− 3 of ammonia, 6 mmol m^− 3 of phosphate, 7 mmol m^− 3 of total manganese). The hydrophysical structure was characterized by the uniform distribution of temperature in the upper 6–7 m mixed layer (UML). Below this a thin (0.4–0.8 m) thermocline layer was observed, just above the anoxic waters. Formation of this phenomenon was connected with that summer weather conditions. Intensive rains led to increased influx of river waters in June. That resulted in large input of allochtonous organic matter (OM) and inorganic nutrients; the latter were consumed on the additional autochthonous organic matter production. In July the weather was characterized by a significant rise in the daily averaged air temperature and large oscillations of temperature during the day. In this period a wind of constant direction was absent, but wind bursts were observed. The completed analyses showed that the formation of such a structure could be connected with the following factors: (i) positive growth trends of the daily averaged temperature and the daily oscillations of temperature, (ii) presence of wind bursts. The joint action of these factors resulted in the formation of the UML. The amplitude of wind bursts determined the depth of UML, and the value of trend determined the value of the temperature change in the thermocline. An initial presence of bottom halocline (caused by the Black Sea water influx to the bottom of the Sea of Azov) prevented the heating of the bottom layer and therefore led to an increase of vertical gradient of temperature in the thermocline. The spatial distribution of the turbulent exchange coefficient confirmed the existence of a “stagnation” area located above the anoxia zone, which is also, apparently, the reason for its occurrence.     

Dynamics of suprabenthos-zooplankton communities around the Balearic Islands (western Mediterranean): Influence of environmental variables and effects on the biological cycle of Aristeus antennatus

Dynamics of suprabenthos and zooplankton were analyzed in two areas located in the NW (off Sóller harbour) and S (off Cabrera Archipelago) of Mallorca (Balearic Islands, western Mediterranean) at depths ranging between 135–780 m. Four stations situated respectively at 150 m (shelf-slope break), and at bathyal depths of 350, 650 and 750 m were sampled at bi-monthly intervals during six cruises performed between August 2003 and June 2004. Suprabenthos showed maximum biomass in both areas from late spring to summer (April to August), while minimum biomass was found in autumn (September–November). Though variable, temporal dynamics of zooplankton showed peaks of biomass in late winter and summer (February and June), while minimals occurred in autumn (August–September) and, at bathyal depths, in April. Suprabenthos (abundance; MDS analyses) showed a sample aggregation as a function of depth (3 groups corresponding to the shelf-slope break, upper slope — over 350 m; and the middle, deeper part of the slope — over 650–750 m), without any separation of hauls by season. By contrast, zooplankton samples were separated by season and not by depth. There was evidence of three seasonal groups corresponding to summer (June 2004–August 2003), autumn–winter (September and November 2003, February 2004), and spring (April 2004), being especially well established off Sóller. In general, suprabenthos was significantly correlated with the sediment variables (e.g. total organic matter content (% OM), potential REDOX), whereas zooplankton was almost exclusively dependent on Chl a at the surface, which suggests two different food sources for suprabenthos and zooplankton. The increase of suprabenthos abundance in April–June was paralleled by a sharp increase (ca. 2.8 times) in the %OM on sediment during the same period, coupled ca. 1–2 months of delay with the peak of surface Chl a recorded in February–March (from satellite imagery data). Suprabenthos biomass was also correlated with salinity close to the bottom, suggesting a link between suprabenthos abundance and changes in the oceanographic condition of water masses close to the bottom. It is suggested that a higher suprabenthos biomass recorded off Sóller in comparison to that off Cabrera in June could, in turn, be related to a seasonal inflow of Levantine Intermediate Water (LIW) in April–June in this area at mid bathyal depths (350–650 m). This trend would be based on: 1) it was evident only at mid-slope depths between 350–750 m, coinciding with the LIW distribution, and 2) it was not recorded among zooplankton (collected throughout the water column). The possible effect of the fluctuations of suprabenthos and zooplankton on higher trophic levels has been explored studying the diet and food consumption rates of the red shrimp Aristeus antennatus, as indicator species by its dominance in bathyal communities. A. antennatus increased its food consumption from February to April–June 2004 off Sóller, which in the case of large (CL > 40 mm) specimens was found in both areas. In addition, there was a shift of diet from winter to spring–early summer. In this last period, A. antennatus preyed upon euphausiids and mesopelagic decapods and fish, while benthos (e.g. polychaetes and bivalves) decreased in the diet. This indicates an increase in the food consumption and probably in the caloric content of the diet in pre-spawning females in April–June 2004, which is synchronized with the period when gonad development begins in A. antennatus females (May–June). Anyway, macrozooplankton, and not suprabenthos, is crucial as a high energetic food source in the coupling between food intake and reproduction in the red shrimp.     

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.     

Swimming for survival: A role of phytoplankton motility in a stratified turbulent environment

We investigate a role for vertical migration in stratified coastal water, where the swimming speed is generally significantly less than the typical turbulent fluctuations in a tidally-mixed bottom layer. In our modelling approach we use a k- turbulence model to describe the physical forcing, a Lagrangian random walk model to describe the vertical displacement of individual cells in response to turbulence and due to cell motility, and a phytoplankton growth model to direct the swimming behaviour of the phytoplankton according to their light and nutrient requirements. The model results show how the cells form a stable subsurface chlorophyll maximum (SCM) at the base of the thermocline where episodic tidal turbulence causes erosion of part of the SCM biomass into the bottom mixed layer (BML). We then focus on the question of whether an ability to swim (weakly, compared to typical bottom layer turbulent intensities) provides any advantage by allowing return to the SCM. Our results show that tidal turbulence in the BML helps both motile and neutrally-buoyant cells by periodically pushing them into the base of the thermocline. Motile cells then have the advantage that they can swim further into the thermocline towards higher light which also reduces the likelihood of being re-mixed back into the BML.     

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.     

Mixing over the steep side of the Cycladic Plateau in the Aegean Sea

Intensive microstructure sampling over the southern slope of the Cycladic Plateau found very weak mixing in the pycnocline, centered on a thin minimum of diapycnal diffusivity with Kρ=1.5×¹⁰−6 m² s^− 1. Below the pycnocline, Kρ increased exponentially in the bottom 200 m, reaching 1 × 10^− 4 m² s^− 1 a few meters above the bottom. Near-bottom mixing was most intense where the bottom slope equaled the characteristic slope of the semi-diurnal internal tide. This suggests internal wave scattering and/or generation at the bottom, a conclusion supported by near-bottom dissipation rates increasing following rising winds and with intensifying internal waves. Several pinnacles on the slope were local mixing hotspots. Signatures included a vertical line of strong mixing in a pinnacle's wake, an hydraulic jump or lee wave over a downstream side of the summit, and a ‘beam’ sloping upward at the near-inertial characteristic slope. Because dissipation rate averages were dominated by strong turbulence, ?/νN² > 100, the effect on Kρ of alternate mixing efficiencies proposed for this range of turbulent intensity is explored.     

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.     

Parameterization of air–sea gas fluxes at extreme wind speeds

Air–sea flux measurements of O₂ and N₂ obtained during Hurricane Frances in September 2004 [D'Asaro, E. A. and McNeil, C. L., 2006. Measurements of air–sea gas exchange at extreme wind speeds. Journal Marine Systems, this edition.] using air-deployed neutrally buoyant floats reveal the first evidence of a new regime of air–sea gas transfer occurring at wind speeds in excess of 35 m s^− 1. In this regime, plumes of bubbles 1 mm and smaller in size are transported down from near the surface of the ocean to greater depths by vertical turbulent currents with speeds up to 20−30 cm s^− 1. These bubble plumes mostly dissolve before reaching a depth of approximately 20 m as a result of hydrostatic compression. Injection of air into the ocean by this mechanism results in the invasion of gases in proportion to their tropospheric molar gas ratios, and further supersaturation of less soluble gases. A new formulation for air–sea fluxes of weakly soluble gases as a function of wind speed is proposed to extend existing formulations [Woolf, D.K, 1997. Bubbles and their role in gas exchange. In: Liss, P.S., and Duce, R.A., (Eds.), The Sea Surface and Global Change. Cambridge University Press, Cambridge, UK, pp. 173–205.] to span the entire natural range of wind speeds over the open ocean, which includes hurricanes. The new formulation has separate contributions to air–sea gas flux from: 1) non-supersaturating near-surface equilibration processes, which include direct transfer associated with the air–sea interface and ventilation associated with surface wave breaking; 2) partial dissolution of bubbles smaller than 1 mm that mix into the ocean via turbulence; and 3) complete dissolution of bubbles of up to 1 mm in size via subduction of bubble plumes. The model can be simplified by combining “surface equilibration” terms that allow exchange of gases into and out of the ocean, and “gas injection” terms that only allow gas to enter the ocean. The model was tested against the Hurricane Frances data set. Although all the model parameters cannot be determined uniquely, some features are clear. The fluxes due to the surface equilibration terms, estimated both from data and from model inversions, increase rapidly at high wind speed but are still far below those predicted using the cubic parameterization of Wanninkhof and McGillis [Wannikhof, R. and McGillis, W.R., 1999. A cubic relationship between air–sea CO₂ exchange and wind speed. Geophysical Research Letters, 26:1889–1892.] at high wind speed. The fluxes due to gas injection terms increase with wind speed even more rapidly, causing bubble injection to dominate at the highest wind speeds.     

Effects of small-scale turbulence on the growth of two diatoms of different size in a phosphorus-limited medium

The effect of turbulence on the nutrient flux towards osmotrophic cells is predicted to be size dependent. This should translate into growth. We experimentally followed and modelled the growth of two marine diatoms of different size (Thalassiosira pseudonana, 6 μm in diameter and Coscinodiscus sp., ca. 109 μm in diameter) under still water and turbulent conditions, using a shaker table. Experiments were done with phosphorus-limited cultures and lasted for ca. 5 days. Turbulence enhanced the growth of Coscinodiscus sp. in agreement with theory but not the growth of T. pseudonana, which was actually slightly lower under turbulence. At the end of the experiments there were about 1.7 times as many Coscinodiscus sp. cells in the turbulent treatment than in the still treatment, while for T. pseudonana almost the same cell concentration was found in both conditions. In addition, the Coscinodiscus sp. cells growing under still conditions presented a higher specific alkaline phosphatase activity than those growing in turbulence which indicates a higher need for phosphorus in the still cultures. A simple dynamic model, based on Michaelis–Menten nutrient uptake kinetics, needed nearly no optimisation other than using observed initial conditions of phosphate and cell concentrations. The model showed how an increased nutrient flux towards the cells translates non-linearly into cell growth, most likely by affecting the half-saturation constant (K_M). However, since Coscinodiscus sp. experienced significant mortality and cells partially settled to the bottom of the containers, unequivocal support for the size-dependent effect of turbulence on nutrient uptake will require further experiments and more sophisticated modelling. The mechanisms to connect an increased nutrient flux towards cells with population growth and whether this process is size dependent are important in parameterizing the effects of turbulence on marine plankton in coupled physical–biological models.     

Influence of the heavy fuel spill from the Prestige tanker wreckage in the overlying seawater column levels of copper, nickel and vanadium (NE Atlantic ocean)

The water column above the Prestige wreckage was sampled during two consecutive campaigns: Prestinaut (December 2002) two weeks after the tanker sunk and HidroPrestige0303 (March 2003) one month after the sealing of the main fuel leaks. Samples of the original cargo fuel and the emulsified fuel in the surface of the ocean were also collected. Analysis of the fuel indicated the release of 135 kg of Cu, 1700 kg of Ni and 5300 kg of V from the original fuel to the water column, remaining 35 kg of Cu, 3100 kg of Ni and 13,800 kg of V in the emulsified fuel. The metal partitioning between the water column and the emulsioned floating fuel, Cu > Ni ~ V, are in accordance with the stability index for the metal–nitrogen bond in metalloporphyrins. This release had an impact on dissolved trace metal concentrations in the water column. An increase on dissolved copper (2.8–4.7 nM) and nickel (2.2–8.0 nM) with respect to natural values (1–3 nM for Cu and 1.6–5 nM for Ni) was observed. Values for vanadium (28–35 nM) were in the range of pristine North Atlantic waters (30–36 nM). This contamination was especially observed in the upper water column (0–50 m), associated with the mixing of seawater with the fuel moving upwards, and in deep waters, where the residence time of fuel is higher. Future research in this field should focus on the environmental variables and the processes that control the release of contaminants from fuels for a better assessment of the contamination in oil-spill events.     

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.     

Characteristics of bubble plumes, bubble-plume bubbles and waves from wind-steepened wave breaking

Observations of breaking waves, associated bubble plumes and bubble-plume size distributions were used to explore the coupled evolution of wave-breaking, wave properties and bubble-plume characteristics. Experiments were made in a large, freshwater, wind-wave channel with mechanical wind-steepened waves and a wind speed of 13 m s^− 1. Bubble plumes exhibited a wide range of bubble distributions, physical extent and dynamics. A classification scheme was developed based on plume extent and “optical density” which is the ability of a plume to optically obscure the image of the background until maximum penetration of the plume. Plumes were classified as either dense (obscure) or diffuse (no-obscure). For each class, the plume bubble population size distribution, Φ(r,t), where r is the bubble radius and t the time, was determined. Dense plumes have a large radius peak in Φ and thus are enhanced in large bubbles. Diffuse plumes are well-described by a weakly size decreasing Φ(r,t) for r < 1000 μm and a more strongly size decreasing Φ(r,t) for r > 1000 μm.The bubble-plume formation rate, P, for each class, wave-breaking rate and wave characteristics were measured with respect to fetch. Wave-breaking rate and intensity are strongly fetch-dependent. In general, the trends in P and wave breaking are similar, reaching a maximum at the fetch of maximum wave breaking. The ratio of P for dense to diffuse plumes is even more sensitive to the occurrence of the most intense wave breaking, where dense plume formation is the greatest.Using P and the bubble size population distributions for each plume class, the global bubble-plume, injection size distribution, Ψ_i(r), was calculated. The volume injection rate for the study area was 640 cm³ s^− 1 divided approximately equally between bubbles smaller and larger than r  1700 μm.     

Recovery of North-East Atlantic temperature fields from profiling floats: Determination of the optimal float number from sampling and instrumental error analysis

Argo is an international project that is deploying an array of temperature and salinity profiling floats over the global ocean. Here we use the error formulation derived from Optimal Statistical Interpolation to estimate statistical errors associated with the recovery of the temperature field in the North-East Atlantic ocean. Results indicate that with the present distribution of floats (119 in the considered domain), scales of wavelength larger than 500 km can be recovered with a relative uncertainty (rms error relative to the standard deviation of the field) of about 7% at 50 m, 8% at 200 m and 10% at 1000 m. This corresponds to mean absolute errors of 0.111 °C at 50 m, 0.104 °C at 200 m and 0.073 °C at 1000 m.The splitting of total errors into instrumental and sampling contributions reveals that, in the present scenario, errors are more due to the small number of floats than to instrumental errors, especially at upper levels. For scales larger than 500 km this will hold true until 200–250 floats are deployed (less than 200 for deep levels). In such a simulated scenario, the number of observations and the technology become approximately equally limiting factors for the accuracy of the temperature field mapping, with total relative errors of less than 2% at upper levels and about 3% at 1000 m.