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.     

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.     

A novel method for estimating vertical eddy diffusivities using diurnal signals with application to western Long Island Sound

We present an approach that allows the estimation of vertical eddy diffusivity coefficients from buoy measurements made at two or more depths. By measuring the attenuation and phase lag of a scalar signal generated periodically at the surface as it propagates downwards, the vertical eddy diffusivity coefficients can be calculated as K_v = ωΔz²/2ln²(α₂/α₁), where α₂/α₁ is the ratio of the real amplitudes at frequency ω at the two depths separated by Δz = z₂ − z₁; as K_V = ωΔz²/2φ², where φ is the phase lag at the frequency ω; or as K_v = iωΔz²/ln²(U₂/U₁), where U₂/U₁ is the ratio of the complex signal amplitudes at the two depths. The method requires that horizontal fluxes be small at the ω frequency and that the signal-to-noise ratios at the two depths allow the determination of the amplitude and phase of ω.Application of this method to summertime 2004 western Long Island Sound oxygen and temperature buoy measurements at two depths provides a time-series of two-day average vertical eddy diffusivity estimates. Using these eddy diffusivities in conjunction with measured vertical concentration gradients, we obtain a time-series of vertical transport rates for oxygen and heat and estimate mean downward fluxes for June and July as 150–260 mMol m^− 2 day^− 1 and 100–400 W m^− 2 respectively. These estimates are of a similar magnitude to sub-pycnocline O₂ and heat demands of 240 ± 200 mMol m^− 2 day^− 1 and 180 ± 60 W m^− 2 that we infer from simple budgets, implying that vertical transport is significant to both budgets.The eddy coefficients obtained from the independent O₂ and temperature measurements have a 68% correlation, and the O₂ flux estimates show a correlation of 41% to measured rates of change in bottom dissolved oxygen levels. Our results indicate that extended time-series of eddy diffusivity coefficients can be obtained from in situ buoy measurements and the method shows promise as a way to constrain the vertical transport variability in budgets of dissolved materials in estuaries.     

Air–sea CO2 fluxes in the Caribbean Sea from 2002–2004

Air–sea fluxes in the Caribbean Sea are presented based on measurements of partial pressure of CO₂ in surface seawater, pCO_2sw, from an automated system onboard the cruise ship Explorer of the Seas for 2002 through 2004. The pCO_2sw values are used to develop algorithms of pCO_2sw based on sea surface temperature (SST) and position. The algorithms are applied to assimilated SST data and remotely sensed winds on a 1° by 1° grid to estimate the fluxes on weekly timescales in the region. The positive relationship between pCO_2sw and SST is lower than the isochemical trend suggesting counteracting effects from biological processes. The relationship varies systematically with location with a stronger dependence further south. Furthermore, the southern area shows significantly lower pCO_2sw in the fall compared to the spring at the same SST, which is attributed to differences in salinity. The annual algorithms for the entire region show a slight trend between 2002 and 2004 suggesting an increase of pCO_2sw over time. This is in accord with the increasing pCO_2sw due the invasion of anthropogenic CO₂. The annual fluxes of CO₂ yield a net invasion of CO₂ to the ocean that ranges from − 0.04 to − 1.2 mol m^− 2 year^− 1 over the 3 years. There is a seasonal reversal in the direction of the flux with CO₂ entering into the ocean during the winter and an evasion during the summer. Year-to-year differences in flux are primarily caused by temperature anomalies in the late winter and spring period resulting in changes in invasion during these seasons. An analysis of pCO_2sw before and after hurricane Frances (September 4–6, 2004), and wind records during the storm suggest a large local enhancement of the flux but minimal influence on annual fluxes in the region.     

Fluxes of dissolved inorganic carbon in three estuarine systems of the Cantabrian Sea (north of Spain)

The diffusive and in situ fluxes of dissolved inorganic carbon (DIC) and total alkalinity (TA) have been measured and an estimation has been made of the water–atmosphere fluxes of CO₂ in three estuarine systems of the Cantabrian Sea during the spring of 1998. Each of these systems undergoes a different anthropogenic influence. The diffusive fluxes of dissolved inorganic carbon and total alkalinity obtained present values ranging between 0.54–2.65 and 0.0–2.4 mmol m⁻² day⁻¹, respectively. These ranges are in agreement with those of other coastal systems. The in situ fluxes are high and extremely variable (35–284 mmol TA m⁻² day⁻¹, 43–554 mmol DIC m⁻² day⁻¹ and 22–261 mmol dissolved oxygen (DO) m⁻² day⁻¹), because the systems studied are very heterogeneous. The values of the ratio of the in situ fluxes of TA and DIC show on average that the rate of dissolution of CaCO₃ is 0.37 times that of organic carbon oxidation. Equally, the interval of variation of the relationship between the benthic fluxes of inorganic carbon and oxygen (F_DIC/F_DO) is very wide (0.3–13.9), which demonstrates the different contributions made by the processes of aerobic and anaerobic degradation of the organic matter, as well as by the dissolution–precipitation of CaCO₃. The water–atmosphere fluxes of CO₂ present a clear dependence on the salinity. The brackish water of these systems (salinity<20), where maximum fluxes of 989 mmol m⁻² day⁻¹ have been estimated, act as a source of CO₂ to the atmosphere. The more saline zones of the estuary (salinity>30) act as a sink of CO₂, with fluxes between −5 and −10 mmol m⁻² day⁻¹.     

Evaluating transfer velocity-wind speed relationship using a long-term series of direct eddy correlation CO2 flux measurements

Long-term observations of the marine atmospheric boundary layer were performed by an eddy correlation system, which was set-up on a platform in the Baltic Sea. In this experiment the three-dimensional wind vector and the turbulent fluxes of momentum, sensible and latent heat and CO₂ were measured for one and a half years. Simultaneously the CO₂ partial pressure p_CO2 in surface water was measured by a submersible autonomous moored instrument for CO₂ at the platform in 7-m depth. The high-resolution eddy correlation measurements of the atmospheric CO₂ flux F_CO2, together with the measurements of the CO₂ partial pressure differences between air and sea Δp_CO2 led to a long-term data set which provided the possibility to investigate the parameterization of the CO₂ transfer velocity k as a function of 10-m wind speed u in a statistical manner. From half-hour mean CO₂ fluxes and CO₂ partial pressure differences, k was calculated using k = F_CO2 / (K₀Δp_CO2), with K₀ the CO₂ solubility. The half-hour mean data points, used for the determination of the k–u parameterization, show large scatter. However, assuming a linear, quadratic dependency the analysis yields: k₆₆₀ = 0.365u² + 0.46u (k at 20 °C and salinity 35 psu) with a correlation coefficient of r² = 0.81. The large scatter indicates that the kinetics of the air–sea CO₂ transfer velocity is not only a function of the wind speed alone, but might also be controlled by other environmental parameters and mechanisms, such as sea state and surface coverage with surfactants.     

An approach to estimation of near-surface turbulence and CO2 transfer velocity from remote sensing data

The air–sea CO₂ exchange is primarily determined by the boundary-layer processes in the near-surface layer of the ocean since it is a water-side limited gas. As a consequence, the interfacial component of the CO₂ transfer velocity can be linked to parameters of turbulence in the near-surface layer of the ocean. The development of remote sensing techniques provides a possibility to quantify the dissipation of the turbulent kinetic energy in the near-surface layer of the ocean and the air–sea CO₂ transfer velocity on a global scale. In this work, the dissipation rate of the turbulent kinetic energy in the near-surface layer of the ocean and its patchiness has been linked to the air–sea CO₂ transfer velocity with a boundary-layer type model. Field observations of upper ocean turbulence, laboratory studies, and the direct CO₂ flux measurements are used to validate the model. The model is then forced with the TOPEX POSEIDON wind speed and significant wave height to demonstrate its applicability for estimating the distribution of the near-surface turbulence dissipation rate and gas transfer velocity for an extended (decadal) time period. A future version of this remote sensing algorithm will incorporate directional wind/wave data being available from QUIKSCAT, a now-cast wave model, and satellite heat fluxes. The inclusion of microwave imagery from the Special Sensor Microwave Imager (SSM/I) and the Synthetic Aperture Radar (SAR) will provide additional information on the fractional whitecap coverage and sea surface turbulence patchiness.     

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.     

Mixing in the surface boundary layer of a tropical freshwater reservoir

A combined observational-modeling study was conducted to investigate turbulence mixing, and the relation to surface forcing, in the surface boundary layer (SBL) of a tropical, high-altitude, freshwater reservoir. A suite of vertical profiles of temperature microstructure, collected at three different stations of one-day duration each, provided estimates of dissipation rates of turbulence kinetic energy, , and temperature variance, χ. Numerical simulations of and χ, using state-of-the-art, public domain, two-equation turbulence closure models, compared favorably with the observations and reproduced the dynamics of daytime wind mixing as well as the vertical and temporal turbulence structure during nighttime convective conditions.Two independent estimates of vertical eddy diffusivities in the stably stratified (daytime) SBL, computed from the microstructure measurements, agreed closely, and the near surface heat and buoyancy fluxes, computed from the diffusivities, were similar to those computed independently from surface meteorology. Model generated eddy diffusivities agreed closely with the observed values, except those generated by K profile parameterization (KPP) model simulations. The good agreement provides confidence that nutrient fluxes in the SBL may be accurately computed from the models when forced with regularly measured surface meteorological parameters. The consequences are important for estimation of daily primary productivity rates in the euphotic zone and the ability to predict algal blooms such as those observed in the present reservoir.     

Simple Lagrangian formulation of bubbly flow in a turbulent boundary layer (bubbly boundary layer flow)

For the theoretical consideration of a system for reducing skin friction, a mathematical model was derived to represent, in a two-phase field, the effect on skin friction of the injection of micro air bubbles into the turbulent boundary layer of a liquid stream. Based on the Lagrangian method, the equation of motion governing a single bubble was derived. The random motion of bubbles in a field initially devoid of bubbles was then traced in three dimensions to estimate void fraction distributions across sections of the flow channel, and to determine local bubble behavior. The liquid phase was modeled on the principle of mixing length. Assuming that the force exerted on the liquid phase was equal to the fluid drag generated by bubble slip, an equation was derived to express the reduction in turbulent shear stress. Corroborating experimental data were obtained from tests using a cavitation tunnel equipped with a slit in the ceiling from which bubbly water was injected. The measurement data provided qualitative substantiation of the trend shown by the calculated results with regard to the skin friction ratio between cases with and without bubble injection as function of the distance downstream from the point of bubble injection.List of symbols B law of wall constant - C _f local coefficient of skin friction - C _f0 local coefficient of skin friction in the absence of bubbles - d _b bubble diameter [m] - g acceleration of gravity [m/s²] - k ₁ k₄ proportional coefficient - k _L turbulent energy of the liquid phase [m²/s²] - L representative length [m] - l _b mean free path of a bubble [m] - m _A added mass of a single bubble [kg] - m _b mass of a single bubble [kg] - N _x ,N _y ,N _z force perpendicular to the wall or ceiling exerted on a bubble adhering to that wall or ceiling [N] - P absolute pressure [Pa] - Q _G rate of air supply [/min] - q _L ⁽ⁱ⁾ turbulent velocity at the ith time increment [m/s] - R> _ex Reynolds number defined by Eq. 32 - T _*L integral time scale of the liquid phase [s] - U velocity of the main stream [m/s] - ,¯v,¯w time-averaged velocity components [m/s] - u,v,w turbulent velocity components [m/s] - û _L ,v_L root mean square values of liquid phase turbulence components in thex- and y-directions [m/s] - V volume of a single bubble [m³] - X,Y,Z components of bubble displacement [m] - x _s ,y _s ,z _s coordinate of a random point on a sphere of unit diameter centered at the coordinate origin - root mean square of bubble displacement in they-direction in reference to the turbulent liquid phase velocity [m] - local void fraction - _m mean void fraction in a turbulent region - regular random number - R _v increment of the horizontal component of the force acting on a single bubble, defined by Eq. 22 [N] - t time increment [s] - ₁ reduction of turbulent stress [N/m²] - _L rate of liquid energy dissipation [m²/s³] - _m coefficient defined by Eq. 30 - law of wall constant in the turbulent region in absence of bubbles - ₁ law of wall constant in the turbulent region in presence of bubbles     

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.     

Hydrodynamic loads on ice-class propellers during propeller-ice interaction

Hydrodynamic loads on a propeller blocked with simulated ice were studied using a cavitation tunnel. Comparative predictions were made using a panel method. The propeller was a model of the Canadian Coast Guard's R-class icebreake propeller, and the ice block was simulated using a solid blockage. Experimental results show the open water performance of the propeller, its performance behind a blockage, and the effects of cavitation in these conditions, as well as the loading on the simulated ice block. Panel method predictions were made of the time series propeller performance in the blocked flow. Cavitation during propellerice interaction resulted in a reduction of mean suction load on the ice block. Block load measurements indicated an increase in the oscillation about the mean value of the loads, with a variation in the phase of the loading with respect to blade position as compared with the non-cavitating results. Comparisons of panel method results with the measured block loads support the reliability of the dynamic measurements.List of symbols D propeller diameter - F block drag load - K _T thrust coefficient,T/(n ² D ⁴) - K _B block load coefficient,F/(n ² D ⁴) - K _Q torque coefficient,Q/(n ² D ⁵) - Q propeller torque - T propeller thrust - n propeller rotational speed - J propeller advance coefficientV _A/(nD) - P _A ambient pressure at propeller - P _ATM atmospheric pressure - P _V vapour pressure of water - V _A propeller advance speed - dissolved gas content - _s saturated dissolved gas content at atmospheric pressure - _o open water propeller efficiency - cavitation number, (P _A –P _V )/(0.5(nD)²) - density of water     

Surf-riding and oscillations of a ship in quartering waves

The behavior of a ship encountering large regular waves from astern at low frequency is the object of investigation, with a parallel study of surf-riding and periodic motion paterns. First, the theoretical analysis of surf-riding is extended from purely following to quartering seas. Steady-state continuation is used to identify all possible surf-riding states for one wavelength. Examination of stability indicates the existence of stable and unstable states and predicts a new type of oscillatory surf-riding. Global analysis is also applied to determine the areas of state space which lead to surf-riding for a given ship and wave conditions. In the case of overtaking waves, the large rudder-yaw-surge oscillations of the vessel are examined, showing the mechanism and conditions responsible for loss of controllability at certain vessel headings.List of symbols c wave celerity (m/s) - C(p) roll damping moment (Ntm) - g acceleration of gravity (m/s²) - GM metacentric height (m) - H wave height (m) - I _x ,I _z roll and yaw ship moments of inertia (kg m²) - k wave number (m^–1) - K _H ,K _W ,K _R hull reaction, wave, rudder, and propeller - K _p forces in the roll direction (Ntm) - m ship mass (kg) - n propeller rate of rotation (rpm) - N _H ,N _W ,N _R hull reaction, wave, rudder, and propeller - N _P moments in the yaw direction (Ntm) - p roll angular velocity (rad/s) - r rate-of-turn (rad/s) - R(,x) restoring moment (Ntm) - Res(u) ship resistance (Nt) - t time (s) - u surge velocity (m/s) - U vessel speed (m/s) - v sway velocity (m/s) - W ship weight (Nt) - x longitudinal position of the ship measured from the wave system (m) - x _G ,z _G longitudinal and vertical center of gravity (m) - x _S longitudinal position of a ship section (S), in the ship-fixed system (m) - X _H ,X _W ,X _R hull reaction, wave, rudder, and propeller - X _P forces in the surge direction (Nt) - y transverse position of the ship, measured from the wave system (m) - Y _H ,Y _W ,Y _R hull reaction, wave, rudder, and propeller - Y _p forces in the sway direction (Nt) - z _Y vertical position of the point of action of the lateral reaction force during turn (m) - z _W vertical position of the point of action of the lateral wave force (m) Greek symbols angle of drift (rad) - rudder angle (rad) - wavelength (m) - position of the ship in the earth-fixed system (m) - water density (kg/m³) - angle of heel (rad) - heading angle (rad) - _e frequency of encounter (rad/s) Hydrodynamic coefficients K roll added mass - N _v ,N _r yaw acceleration coefficients - N _v N _r N _rr N _rrv ,N _vvr yaw velocity coefficients K. Spyrou: Ship behavior in quartering waves - X _u surge acceleration coefficient - X _u X _vr surge velocity coefficients - Y _v ,Y _r sway acceleration coefficients - Y _v ,Y _r ,Y _vv ,Y _rr ,Y _vr sway velocity coefficients European Union-nominated Fellow of the Science and Technology Agency of Japan, Visiting Researcher, National Research Institute of Fisheries Engineering of Japan     

Distribution of void fraction in bubbly flow through a horizontal channel: Bubbly boundary layer flow, 2nd report

A method of enveloping the hull with a sheet of microbubbles is discussed. It forms part of a study on means of reducing the skin friction acting on a ship's hull. In this report, a bubble traveling through a horizontal channel is regarded as a diffusive particle. Based on this assumption, an equation based on flow flux balance is derived for determining the void fraction in approximation. The equation thus derived is used for calculation, and the calculation results are compared with reported experimental data. The equation is further manipulated to make it compatible with a mixing length model that takes into account the presence of bubbles in the liquid stream. Among the factors contained in the equation thus derived, those affected by the presence of bubbles are the change of mixing length and the difference in the ratio of skin friction between cases with and without bubbles. These factors can be calculated using the mean void fraction in the boundary layer determined by the rate of air supply into the flow field. It is suggested that the ratio between boundary layer thickness and bubble diameter could constitute a significant parameter to replace the scale effect in estimating values applicable to actual ships from corresponding data obtained in model experiments.List of symbols a ₁ proportionality constant indicating directionality of turbulence - B law-of-the-wall constant - C _f local skin-friction coefficient in the presence of bubbles - C _f0 local skin-friction coefficient in the absence of bubbles - d _b bubble diameter (m) - g acceleration of gravity (m/s²) - j _g flow flux of gas phase accountable to buoyancy (m/s) - j _t flow flux of gas phase accountable to turbulence (m/s) - k ₄ constant relating reduction of liquid shear stress by bubble presence to decrease of force imparted to bubble by its displacement due to turbulence - l _b mixing length of gas phase (m) - l _m mixing length of liquid phase (m) - l _mb diminution of liquid phase mixing length by bubble presence (m) - Q _G rate of air supply to liquid stream (l/min) - q _/g velocity of bubble rise (m/s) - 2R height of horizontal channel (m) - T _* integral time scale (s) - U _m mean stream velocity in channel (m/s) - U friction velocity in channel (m/s) - V volume of a bubble (m³) - u, ¯ v time-averaged stream velocities inx- andy-directions, respectively (m/s) - u, v turbulent velocity components inx- andy-directions, respectively (m/s) - v root mean square of turbulence component in they-direction (m/s) - root mean square of bubble displacement iny-direction with reference to turbulent liquid phase velocity (m) - y displacement from ceiling (m) - local void fraction - _m mean void fraction in boundary layer - _m constant relating local void fraction to law-of-the-wall constant - _t reduction of turbulent stress (N/m²) - law-of-the-wall constant in turbulent liquid region in absence of bubbles - ₁ law-of-the-wall constant in turbulent liquid region in presence of bubbles - ₂ law-of-the-wall constant in gas phase - _m constant indicating representative turbulence scale (m) - viscosity (Pa × s) - v kinematic viscosity (m²/s) - density (kg/m³) Suffixes G gas - L liquid - 0 absence of bubbles     

Carbonate chemistry dynamics and carbon dioxide fluxes across the atmosphere–ice–water interfaces in the Arctic Ocean: Pacific sector of the Arctic

Climatic changes in the Northern Hemisphere have led to remarkable environmental changes in the Arctic Ocean, which is surrounded by permafrost. These changes include significant shrinking of sea-ice cover in summer, increased time between sea-ice break-up and freeze-up, and Arctic surface water freshening and warming associated with melting sea-ice, thawing permafrost, and increased runoff. These changes are commonly attributed to the greenhouse effect resulting from increased atmospheric carbon dioxide (CO₂) concentration and other non-CO₂ radiatively active gases (methane, nitrous oxide). The greenhouse effect should be most pronounced in the Arctic where the largest air CO₂ concentrations and winter–summer variations in the world for a clean background environment were detected. However, the air–land–shelf interaction in the Arctic has a substantial impact on the composition of the overlying atmosphere; as the permafrost thaws, a significant amount of old terrestrial carbon becomes available for biogeochemical cycling and oxidation to CO₂. The Arctic Ocean's role in determining regional CO₂ balance has been ignored, because of its small size (only  4% of the world ocean area) and because its continuous sea-ice cover is considered to impede gaseous exchange with the atmosphere so efficiently that no global climate models include CO₂ exchange over sea-ice. In this paper we show that: (1) the Arctic shelf seas (the Laptev and East-Siberian seas) may become a strong source of atmospheric CO₂ because of oxidation of bio-available eroded terrestrial carbon and river transport; (2) the Chukchi Sea shelf exhibits the strong uptake of atmospheric CO₂; (3) the sea-ice melt ponds and open brine channels form an important spring/summer air CO₂ sink that also must be included in any Arctic regional CO₂ budget. Both the direction and amount of CO₂ transfer between air and sea during open water season may be different from transfer during freezing and thawing, or during winter when CO₂ accumulates beneath Arctic sea-ice; (4) direct measurements beneath the sea ice gave two initial results. First, a drastic pCO₂ decrease from 410 μatm to 288 μatm, which was recorded in February–March beneath the fast ice near Barrow using the SAMI-CO₂ sensor, may reflect increased photosynthetic activity beneath sea-ice just after polar sunrise. Second, new measurements made in summer 2005 beneath the sea ice in the Central Basin show relatively high values of pCO₂ ranging between 425 μatm and 475 μatm, values, which are larger than the mean atmospheric value in the Arctic in summertime. The sources of those high values are supposed to be: high rates of bacterial respiration, import of the Upper Halocline Water (UHW) from the Chukchi Sea (CS) where values of pCO₂ range between 400 and 600 μatm, a contribution from the Lena river plume, or any combination of these sources.     

Gas transfer velocities of CO2 and CH4 in a tropical reservoir and its river downstream

We have measured simultaneously the methane (CH₄) and carbon dioxide (CO₂) surface concentrations and water–air fluxes by floating chambers (FC) in the Petit-Saut Reservoir (French Guiana) and its tidal river (Sinnamary River) downstream of the dam, during the two field experiments in wet (May 2003) and dry season (December 2003). The eddy covariance (EC) technique was also used for CO₂ fluxes on the lake. The comparison of fluxes obtained by FC and EC showed little discrepancies mainly due to differences in measurements durations which resulted in different average wind speeds. When comparing the gas transfer velocity (k₆₀₀) for a given wind speed, both methods gave similar results. On the lake and excluding rainy events, we obtained an exponential relationship between k₆₀₀ and U₁₀, with a significant intercept at 1.7 cm h^− 1, probably due to thermal effects. Gas transfer velocity was also positively related to rainfall rates reaching 26.5 cm h⁻¹ for a rainfall rate of 36 mm h^− 1. During a 24-h experiment in dry season, rainfall accounted for as much as 25% of the k₆₀₀. In the river downstream of the dam, k₆₀₀ values were 3 to 4 times higher than on the lake, and followed a linear relationship with U₁₀.     

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

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

CO2 partial pressure in Northeast Atlantic and adjacent shelf waters: Processes and seasonal variability

CO₂ partial pressure in surface water was measured in the Northeast Atlantic and in the Hebride Shelf/North Sea area during a cruise with R.V. Poseidon in June 1991. A mean p_CO2 of 303 μatm was found in the Atlantic between 50°N and 60°N. For an atmospheric CO₂ content of 357.5 ppm(v) this corresponds to a partial pressure difference of −55 μatm. This supports the view that the subarctic Atlantic is a significant sink within the CO₂ cycle between the ocean and the atmosphere. A comparison of our measurements with other data reveals that the p_CO2 distribution changes significantly during May/June. This explained by seasonal warming, CO₂ exchange with the atmosphere and biomass production. The contribution by each of these processes to the seasonal variations is calculated. It was found that during a plankton bloom the production of biomass is the dominating factor and may lower seawater p_CO2 by almost 100 μatm. The shelf areas are charactrized by strong p_CO2 gradients which are explained by water exchange with the Atlantic, temperature effects and biomass production.     

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.     

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.