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.     

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.     

Modeling gas dissolution in deepwater oil/gas spills

Gases in deepwater oil/gas spills can lose considerable amounts of the gas phase due to dissolution in water. Gas dissolution has a significant impact on the behavior of the oil/gas jet/plume because of its impact on the buoyancy. A method is presented in this paper for computing gas dissolution that covers a broad range of water depth, from shallow water where gases behave as ideal ones under low pressure to deepwater where gases behave as non-ideal ones under high pressures. The method presented also accounts for the spherical and non-spherical shapes of gas bubbles. The gas dissolution computations are validated by comparing the computed results with observed data from previously conducted laboratory experiments. The gas dissolution computation module is then integrated with a model for underwater oil/gas jets/plumes by Yapa and Zheng [J. Hydraul. Res. 35 (5) (1997) 673]. Scenario simulations are presented to show the impacts of gas dissolution on the behavior of jets/plumes. These scenarios show the impact of dissolution on the behavior of the jet/plume. The comparison of results using ideal gas conditions and non-ideal gas conditions is also shown.     

The combined effect of rain and wind on air–water gas exchange: A feasibility study

A series of experiments were conducted at University of Delaware's Air–Sea Interaction Laboratory to examine the combined effects of rain and wind on air–water gas exchange. During this study, ASIL WRX I, a combination of 3 rain rates and 4 wind speeds were used, for a total of 12 different environmental conditions. The SF₆ evasion method was used to determine the bulk gas transfer velocities, and airside profiles of wind and CO₂ were used to estimate flux–profiles of momentum and carbon dioxide. In addition to measurements of fluxes with and without rain in a wind–wave boundary layer, measurements of wave properties were also obtained. Rain is shown to alter the wind profile in the flume, and dampen surface waves. Also, SF₆ evasion indicates that with the present experimental setup, for most of the experimental conditions, rain and wind combine linearly to influence air–water gas exchange. Flux–profile relationships for marine atmospheric boundary layers, which were performed to scale up to field measurements, were explored by a comparison between SF₆-derived bulk fluxes and airside CO₂ profile measurements.     

Lengthscales of motions that control air–water gas transfer in grid-stirred turbulence

The relationship between the gas transfer velocity and turbulent lengthscales is investigated experimentally in a grid-stirred turbulent flow. The horizontal velocity field at the water surface is measured using particle image velocimetry (PIV). The gas transfer velocity for oxygen is obtained through reaeration experiments. In addition, the gas transfer process by surface-renewal eddies is visualized using laser-induced fluorescence (LIF) technique, in which carbon dioxide is used as the tracer gas. The definition of the Taylor microscale holds that the root-mean-square (RMS) of the surface divergence is expressed by the square root of the turbulent kinetic energy divided by the Taylor microscale. Experimentally obtained data support this scaling. They show the gas transfer velocity to be proportional to the square root of the RMS of the surface divergence. These experimental results imply that the Taylor microscale is an important parameter for gas transfer velocity at the air–water interface. These relations indicate that a nondimensional gas transfer velocity is proportional to the − 1/4 power of a turbulent-macroscale Reynolds number, which is similar to a small-eddy model, assuming that turbulent eddies with the Kolmogorov scale control the gas transfer process. However, this Reynolds number dependence does not necessarily mean the superiority of turbulent eddies with the Kolmogorov scale in the gas transfer. The LIF visualizations in horizontal and vertical planes close to the air–water interface indicate that the horizontal CO₂-concentration field has a fine spatial pattern, which resembles that of the surface divergence field, and that surface-renewal motions observed in the vertical plane have a larger lengthscale than the Kolmogorov scale. We infer from both PIV and LIF results that the Taylor microscale is an important lengthscale for air–water gas transfer.     

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.     

A conceptual model of the strongly tidal Columbia River plume

The Columbia River plume is typical of large-scale, high discharge, mid-latitude plumes. In the absence of strong upwelling winds, freshwater from the river executes a rightward turn and forms an anticyclonic bulge before moving north along the Washington coast. In addition to the above dynamics, however, the river plume outflow is subject to large tides, which modify the structure of the plume in the region near the river mouth. Observations based on data acquired during a summer 2005 cruise indicate that the plume consists of four distinct water masses; source water at the lift-off point, and the tidal, re-circulating and far-field plumes. In contrast to most plume models that describe the discharge of low-salinity estuary water into ambient high-salinity coastal water, we describe the Columbia plume as the superposition of these four plume types.We focus primarily on a conceptual summary of the dynamics and mutual interaction of the tidal and re-circulating plumes. The new tidal plume flows over top of the re-circulating plume and is typically bounded by strong fronts. Soon after the end of ebb tide, it covers roughly 50–100% of the re-circulating plume surface area. The fronts may penetrate well below the re-circulating plume water and eventually spawn internal waves that mix the re-circulating plume further. The re-circulating plume persists throughout the tidal cycle and corresponds to a freshwater volume equivalent to 3–4 days of river discharge. Finally, the plume water masses are distinguished from one another in term of surface chlorophyll concentration, suggesting that the above classification may also describe different biological growth regimes. The low-salinity re-circulating plume serves as an extension of the estuary into the coastal ocean, or an “estuary at sea”, because residence times during periods of high river flow are greater than those in the estuary.     

Modeling the flow regime near the source in underwater gas releases

Recent progress in calculating gas bubble sizes in a plume, based on phenomenological approaches using the release conditions is a significant improvement to make the gas plume models self-reliant. Such calculations require details of conditions Near the Source of Plume (NSP); (i.e. the plume/jet velocity and radius near the source), which inspired the present work. Determining NSP conditions for gas plumes are far more complex than that for oil plumes due to the substantial density difference between gas and water. To calculate NSP conditions, modeling the early stage of the plume is important. A novel method of modeling the early stage of an underwater gas release is presented here. Major impact of the present work is to define the correct NSP conditions for underwater gas releases, which is not possible with available methods as those techniques are not based on the physics of flow region near the source of the plume/jet. We introduce super Gaussian profiles to model the density and velocity variations of the early stages of plume, coupled with the laws of fluid mechanics to define profile parameters. This new approach, models the velocity profile variation from near uniform, across the section at the release point to Gaussian some distance away. The comparisons show that experimental data agrees well with the computations.     

水下气体释放源附近的流型建模（英文）

Recent progress in calculating gas bubble sizes in a plume, based on phenomenological approaches using the release conditions is a significant improvement to make the gas plume models self-reliant. Such calculations require details of conditions Near the Source of Plume(NSP);(i.e. the plume/jet velocity and radius near the source), which inspired the present work. Determining NSP conditions for gas plumes are far more complex than that for oil plumes due to the substantial density difference between gas and water. To calculate NSP conditions, modeling the early stage of the plume is important. A novel method of modeling the early stage of an underwater gas release is presented here. Major impact of the present work is to define the correct NSP conditions for underwater gas releases, which is not possible with available methods as those techniques are not based on the physics of flow region near the source of the plume/jet. We introduce super Gaussian profiles to model the density and velocity variations of the early stages of plume, coupled with the laws of fluid mechanics to define profile parameters. This new approach, models the velocity profile variation from near uniform, across the section at the release point to Gaussian some distance away. The comparisons show that experimental data agrees well with the computations.     

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.     

A study on bubble plume behavior in stratified water

Plumes of air and carbon-dioxide (CO₂) bubbles in stratified water were studied experimentally and numerically. It is important to understand the plume behavior of droplets or bubbles in the ocean in marine environmental engineering. In sequestration of CO₂ in the ocean, liquid CO₂ is injected in the form of droplets, and thermal stratification in enclosed seas is possibly destroyed by a bubble generator. This study focuses on the relationship between intrusion depth and stratification intensity, gas flow rate, and bubble size. A desktop-sized tank was used to achieve no-background-flow conditions for salt stratification. The results from the air-bubble experiments indicated that the larger the bubble size, the smaller the intrusion depth. We also observed the behavior of CO₂ bubbles. They are different from air in that they dissolve in water and are reduced in volume. Our numerical simulation method for two-phase flow was validated by comparisons with the experiments. Received: August 13, 2001 / Accepted: October 9, 2001     

Inorganic carbon fluxes through the boundaries of the Greenland Sea Basin based on in situ observations and water transport estimates

A carbon budget for the exchange of total dissolved inorganic carbon C_T between the Greenland Sea and the surrounding seas has been constructed for winter and summer situations. An extensive data set of C_T collected over the years 1994–1997 within the European Sub-polar Ocean Programmes (ESOP1 and ESOP2) are used for the budget calculation. Based on these data, mean values of C_T in eight different boxes representing the inflow and outflow of water through the boundaries of the Greenland Sea Basin are estimated. The obtained values are then combined with simulated water transports taken from the ESOP2 version of the Miami Isopycnic Coordinate Ocean Model (MICOM). The fluxes of inorganic carbon are presented for three layers; a surface mixed layer, an intermediate layer and a deep layer, and the imbalance in the fluxes are attributed to air–sea exchange, biological fixation of inorganic carbon, and sedimentation. The main influx of carbon is found in the surface and the deep layers in the Fram Strait, and in the surface waters of direct Atlantic origin, whereas the main outflux is found in the surface layer over the Jan Mayen Fracture Zone and the Knipovich Ridge, transporting carbon into the Atlantic Ocean via the Denmark Strait and towards the Arctic Ocean via the Norwegian Sea, respectively. The flux calculation indicates that there is a net transport of carbon out of the Greenland Sea during wintertime. In the absence of biological activity, this imbalance is attributed to air sea exchange, and requires an oceanic uptake of CO₂ of 0.024±0.006 Gt C yr⁻¹. The flux calculations from the summer period are complicated by biological fixation of inorganic carbon, and show that data on organic carbon is required in order to estimate the air–sea exchange in the area.     

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.     

Air–sea gas exchange at extreme wind speeds measured by autonomous oceanographic floats

Measurements of the air–sea fluxes of N₂ and O₂ were made in winds of 15–57 m s^− 1 beneath Hurricane Frances using two types of air-deployed neutrally buoyant and profiling underwater floats. Two “Lagrangian floats” measured O₂ and total gas tension (GT) in pre-storm and post-storm profiles and in the actively turbulent mixed layer during the storm. A single “EM-APEX float” profiled continuously from 30 to 200 m before, during and after the storm. All floats measured temperature and salinity. N₂ concentrations were computed from GT and O₂ after correcting for instrumental effects. Gas fluxes were computed by three methods. First, a one-dimensional mixed layer budget diagnosed the changes in mixed layer concentrations given the pre-storm profile and a time varying mixed layer depth. This model was calibrated using temperature and salinity data. The difference between the predicted mixed layer concentrations of O₂ and N₂ and those measured was attributed to air–sea gas fluxes F_BO and F_BN. Second, the covariance flux F_CO(z) = wO₂′(z) was computed, where w is the vertical motion of the water-following Lagrangian floats, O₂′ is a high-pass filtered O₂ concentration and (z) is an average over covariance pairs as a function of depth. The profile F_CO(z) was extrapolated to the surface to yield the surface O₂ flux F_CO(0). Third, a deficit of O₂ was found in the upper few meters of the ocean at the height of the storm. A flux F_SO, moving O₂ out of the ocean, was calculated by dividing this deficit by the residence time of the water in this layer, inferred from the Lagrangian floats. The three methods gave generally consistent results. At the highest winds, gas transfer is dominated by bubbles created by surface wave breaking, injected into the ocean by large-scale turbulent eddies and dissolving near 10-m depth. This conclusion is supported by observations of fluxes into the ocean despite its supersaturation; by the molar flux ratio F_BO/F_BN, which is closer to that of air rather than that appropriate for Schmidt number scaling; by O₂ increases at about 10-m depth along the water trajectories accompanied by a reduction in void fraction as measured by conductivity; and from the profile of F_CO(z), which peaks near 10 m instead of at the surface.At the highest winds O₂ and N₂ are injected into the ocean by bubbles dissolving at depth. This, plus entrainment of gas-rich water from below, supersaturates the mixed layer causing gas to flux out of the near-surface ocean. A net influx of gas results from the balance of these two competing processes. At lower speeds, the total gas fluxes, F_BO, F_BN and F_CO(0), are out of the ocean and downgradient.     

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.     

Measurement of air–water CO2 transfer at four coastal sites using a chamber method

We measured the air–water CO₂ flux in four coastal regions (two coral reefs, one estuary, and one coastal brackish lake) using a chamber method, which has the highest spatial resolution of the methods available for measuring coastal air–water gas flux. Some of the measurements were considerably higher than expected from reported wind-dependent relationships. The average k₆₀₀ values for Shiraho Reef, Fukido Reef, Fukido River, and Lake Nakaumi were 1.5 ± 0.6, 3.2 ± 0.3, 0.69 ± 0.26, and 2.2 ± 0.9 (mean ± S.D.) times larger than the wind-dependent relationships. Results were compared with current-dependent relationships and vertical turbulence intensity (VTI). VTI is an index of water-surface stirring and is calculated from near-surface vertical velocity. Although some measurements from the reefs and river closely matched those expected from wind-dependent relationships, others were considerably higher. All data were correlated with VTI and were qualitatively explained by bottom macro-roughness enhancement. In Lake Nakaumi, results tended to differ from the wind-dependent relationships, and the difference between the measured and expected gas-transfer velocity was correlated with biological DO changes and/or the intensity of density stratification. We found these factors to have important effects on coastal gas flux. In addition, the chamber method was an effective tool for evaluating coastal gas flux.     

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

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

Reduction of skin friction by microbubbles and its relation with near-wall bubble concentration in a channel

Determination of the flow structure near the wall is essential for a clear insight into the phenomenon of skin friction reduction by microbubbles in a turbulent boundary layer. An important parameter, is the bubble concentration or void fraction in the wall region in drag-reducing conditions. The purpose of this paper is to show drag-reducing effects due to microbubbles in a water channel and, more importantly, to show the dependence of the drag-reduction values on the near-wall void fraction. A two-dimensional channel with an aspect ratio of 10 was specially built for this purpose with provisions for air injection through porous plates. Skin friction was directly measured by a miniature floating element transducer with a 5-mm circular sensing disk mounted flush on the top wall 67 channel-heights downstream of the injector. The wall friction in the presence of air bubbles was found to be reduced under the same bulk velocity when compared with the value without air. Detailed void fraction profiles across the channel were obtained by a sampling probe and a fiber-optic probe. Better collapse of the drag reduction data, independent of different profile shapes, was found when plotted against the near-wall void fraction than against a cross-sectional mean void fraction. While this dependence reconfirms that the phenomena are essentially inner-region dependent, the lack of influence of the bubble distribution patterns away from the wall implies lack of outer region influence.     

Seasonal variation of riverine nutrient inputs in the northern Bay of Biscay (France), and patterns of marine phytoplankton response

Seasonal variations in nutrient inputs are described for the main rivers (Loire and Vilaine) flowing into the northern Bay of Biscay. The river plumes are high in N/P ratio in late winter and spring, but not in the inner plume during the summer. Conservative behavior results in most nutrients entering the estuary and eventually reaching the coastal zone. Temporal and spatial aspects of phytoplankton growth and nutrient uptake in the northern Bay of Biscay distinguish the central area of salinity 34 from the plume area. The first diatom bloom appears offshore in late winter, at the edge of the river plumes, taking advantage of haline stratification and anticyclonic “weather windows.” In spring, when the central area of the northern shelf is phosphorus-limited, small cells predominate in the phytoplankton community and compete with bacteria for both mineral and organic phosphorus. At that period, river plumes are less extensive than in winter, but local nutrient enrichment at the river mouth allows diatom growth. In summer, phytoplankton become nitrogen-limited in the river plumes; the central area of the shelf is occupied by small forms of phytoplankton, which are located on the thermocline and use predominantly regenerated nutrients.     

Asymmetry of Columbia River tidal plume fronts

Columbia River tidal plume dynamics can be explained in terms of two asymmetries related to plume-front depth and internal wave generation. These asymmetries may be an important factor contributing to the observed greater primary productivity and phytoplankton standing crop on the Washington shelf. The tidal plume (the most recent ebb outflow from the estuary) is initially supercritical with respect to the frontal internal Froude number F_R on strong ebbs. It is separated from the rotating plume bulge by a front, whose properties are very different under upwelling vs. downwelling conditions. Under summer upwelling conditions, tidal plume fronts are sharp and narrow (< 20–50 m wide) on their upwind or northern side and mark a transition from supercritical to subcritical flow for up to 12 h after high water. Such sharp fronts are a source of turbulent mixing, despite the strong stratification. Because the tidal plume may overlie newly upwelled waters, these fronts can mix nutrients into the plume. Symmetry would suggest that there should be a sharp front south of the estuary mouth under summer downwelling conditions. Instead, the downwelling tidal plume front is usually diffuse on its upstream side. Mixing is weaker, and the water masses immediately below are low in nutrients. There is also an upwelling–downwelling asymmetry in internal wave generation. During upwelling and weak wind conditions, plume fronts often generate trains of non-linear internal waves as they transition from a supercritical to a subcritical state. Under downwelling conditions, internal wave release is less common and the waves are less energetic. Furthermore, regardless of wind conditions, solition formation almost always begins on the south side of the plume so that the front “unzips” from south to north. This distinction is important, because these internal waves contribute to vertical mixing in the plume bulge and transport low-salinity water across the tidal plume into the plume bulge.F_R and plume depth are key parameters in distinguishing the upwelling and downwelling situations, and these two asymmetries can be explained in terms of potential vorticity conservation. The divergence of the tidal outflow after it leaves the estuary embeds relative vorticity in the emerging tidal plume water mass. This vorticity controls the transition of the tidal plume front to a subcritical state and consequently the timing and location of internal wave generation by plume fronts.