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.     

The Mercator global ocean operational analysis system: Assessment and validation of an 11-year reanalysis

This paper presents Prototype Système 2 Global (PSY2G), the first Mercator global Ocean General Circulation Model (OGCM) to assimilate along-track sea level anomaly (SLA) satellite data. Based on a coarse resolution ocean model, this system was developed mainly for climatic purposes and will provide, for example, initial oceanic states for coupled ocean-atmosphere seasonal predictions. It has been operational since 3 September 2003 and produces an analysis and a two-week forecast for the global ocean every week. The PSY2G system uses an incremental assimilation scheme based on the Cooper and Haines [Cooper, M., Haines, K., 1996. Data assimilation with water property conservation. J. Geophys. Res., 101, 1059-1077.] lifting–lowering of isopycnals. The SLA increment is obtained using an optimal interpolation method then the correction is partitioned into baroclinic and barotropic contributions. The baroclinic ocean state correction consists of temperature, salinity and geostrophic velocity increments and the barotropic correction is a barotropic velocity increment. A reanalysis (1993–2003) was carried out that enabled the PSY2G system to perform its first operational cycle. All available SLA data sets (TOPEX/Poséïdon, ERS2, Geosat-Follow-On, Jason1 and Envisat) were assimilated for the 1993–2003 period. The major objective of this study is to assess the reanalysis from both an assimilation and a thermodynamic point of view in order to evaluate its realism, especially in the tropical band which is a key region for climatic studies. Although the system is also able to deliver forecasts, we have mainly focused on analysis. These results are useful because they give an a priori estimation of the qualities and capabilities of the operational ocean analysis system that has been implemented. In particular, the reanalysis identifies some regional biases in sea level variability such as near the Antarctic Circumpolar Current, in the eastern Equatorial Pacific and in the Norwegian Sea (generally less than 1 cm) with a small seasonal cycle. This is attributed to changes in mean circulation and vertical stratification caused by the assimilation methodology. But the model's low resolution, inaccurate physical parameterisations (especially for ocean–ice interactions) and surface atmospheric forcing also contribute to the occurrence of the SLA biases. A detailed analysis of the thermohaline structure of the ocean reveals that the isopycnal lifting–lowering tends to diffuse vertically the main thermocline. The impact on temperature is that the surface layer (0–200 m) becomes cooler whereas in deeper waters (from 500 to 1500 m), the ocean becomes slightly warmer. This is particularly true in the tropics, between 30°N and 30°S. However it can be demonstrated that the assimilation improves the variability in both surface currents and sub-surface temperature in the Equatorial Pacific Ocean.     

Interannual variation of the Polar Front in the Japan/East Sea from summertime hydrography and sea level data

The Polar Front in the Japan/East Sea separates the southern warm water region from the northern cold water region. A merged TOPEX/POSEIDON and ERS-1/2 altimeter dataset and upper water temperature data were used to determine the frontal location and to examine the structure of its interannual variability from 1993 to 2001. The identified frontal location, where sea surface height gradient has a maximum about 10–20 cm over the horizontal distance of 100 km, corresponds well to the maximum subsurface horizontal temperature gradient. The front migrates more widely (36°N–41°N) in the western part of the sea than in the eastern part. The interannual migration induces large variability in upper water temperatures and sea surface height in the western region. Responsible physical mechanisms were studied using a reduced-gravity model. Differences between inflow and outflow change the total volume of warm water, and total warm water volume change in the warm water region uniformly pushes the front in the meridional direction across its mean position in the model simulation. Interannual variation of wind stress causes relatively wide migration of the modeled front in the western part.     

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.     

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

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

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₁₀.     

Decadal sea level trends in the Bay of Biscay from tide gauges, GPS and TOPEX

Sea level time series derived from TOPEX altimeter and from tide gauge measurements in the Bay of Biscay (Eastern North Atlantic Ocean) are used to investigate the regional sea level rise. Altimetry sea level anomalies are computed from TOPEX measurements and resampled into fixed along track bins in order to obtain consistent records and to approach as close as possible to the coast. Tidal corrections are critical in the estimation of sea level trends; therefore an additional analysis has been carried out in order to identify and correct for residual tidal signals. The obtained mean sea level rise in the region is + 3.09 ± 0.21 mm/year over the period 1993–2002. Tide gauges have been corrected for vertical land motions by means of collocated GPS measurements. GPS-corrected tide gauges and nearby altimetry trends are found to be statistically consistent. The relationship with the atmospheric forcing is also investigated revealing that 15% of this sea level rise is attributed to the atmospheric pressure effects.     

Picophytoplankton and nanophytoplankton abundance and distribution in the southeastern Beaufort Sea (Mackenzie Shelf and Amundsen Gulf) during Fall 2002

The distribution of picophytoplankton (0.2–2 µm) and nanophytoplankton (2–20 µm) in the Beaufort Sea–Mackenzie Shelf and Amundsen Gulf regions during autumn, 2002 is examined relative to their ambient water mass properties (salinity, temperature and nutrients: nitrate + nitrite, phosphate, and silicate) and to the ratio of variable to maximum fluorescence, Fv/Fm. Total phytoplankton and cell abundances (< 20 µm) were mainly correlated with salinity. Significant differences in picophytoplankton cell numbers were found among waters near the mouth of the Mackenzie River, ice melt waters and the underlying halocline water masses of Pacific origin. Picophytoplankton was the most abundant phytoplankton fraction during the autumnal season, probably reflecting low nitrate concentrations (surface waters average ~ 0.65 µM). The ratio Fv/Fm averaged 0.44, indicating that cells were still physiologically active, even though their concentrations were low (max Chl a = 0.9 mg m^− 3). No significant differences in Fv/Fm were evident in the different water masses, indicating that rate limiting conditions for photosynthesis and growth were uniform across the whole system, which was in a pre-winter stage, and was probably already experiencing light limitation as a result of shortening day lengths.     

Effect of wave boundary layer on sea-to-air dimethylsulfide transfer velocity during typhoon passage

A full-spectral third-generation ocean wind–wave model (Wavewatch-III) implemented in the South China Sea is used to investigate the effects of the wave boundary layer on the drag coefficient and the sea-to-air transfer velocity of dimethylsulfide (DMS) during passage of Typhoon Wukong (September 5–11, 2000) with a maximum sustained wind speed of 38 m s^− 1. The model is driven by the reanalyzed surface winds (1° × 1°, four times daily) from the National Centers for Environmental Prediction. It is found that the wave boundary layer evidently enhances (16.5%) the drag coefficient (in turn increases the momentum flux across the air–sea interface), and reduces (13.1%) the sea-to-air DMS transfer velocity (in turn decreases the sea-to-air DMS flux). This indicates the possibility of important roles of wave boundary layer in atmospheric DMS contents and global climate system.     

1/32° real-time global ocean prediction and value-added over 1/16° resolution

A 1/32° global ocean nowcast/forecast system has been developed by the Naval Research Laboratory at the Stennis Space Center. It started running at the Naval Oceanographic Office in near real-time on 1 Nov. 2003 and has been running daily in real-time since 1 Mar. 2005. It became an operational system on 6 March 2006, replacing the existing 1/16° system which ceased operation on 12 March 2006. Both systems use the NRL Layered Ocean Model (NLOM) with assimilation of sea surface height from satellite altimeters and sea surface temperature from multi-channel satellite infrared radiometers. Real-time and archived results are available online at http://www.ocean.nrlssc.navy.mil/global_nlom. The 1/32° system has improvements over the earlier system that can be grouped into two categories: (1) better resolution and representation of dynamical processes and (2) design modifications. The design modifications are the result of accrued knowledge since the development of the earlier 1/16° system. The improved horizontal resolution of the 1/32° system has significant dynamical benefits which increase the ability of the model to accurately nowcast and skillfully forecast. At the finer resolution, current pathways and their transports become more accurate, the sea surface height (SSH) variability increases and becomes more realistic and even the global ocean circulation experiences some changes (including inter-basin exchange). These improvements make the 1/32° system a better dynamical interpolator of assimilated satellite altimeter track data, using a one-day model forecast as the first guess. The result is quantitatively more accurate nowcasts, as is illustrated by several model-data comparisons. Based on comparisons with ocean color imagery in the northwestern Arabian Sea and the Gulf of Oman, the 1/32° system has even demonstrated the ability to map small eddies, 25–75 km in diameter, with 70% reliability and a median eddy center location error of 22.5 km, a surprising and unanticipated result from assimilation of altimeter track data. For all of the eddies (50% small eddies), the reliability was 80% and the median eddy center location error was 29 km. The 1/32° system also exhibits improved forecast skill in relation to the 1/16° system. This is due to (a) a more accurate initial condition for the forecast and (b) better resolution and representation of critical dynamical processes (such as upper ocean – topographic coupling via mesoscale flow instabilities) which allow the model to more accurately evolve these features in time while running in forecast mode (forecast atmospheric forcing for the first 5 days, then gradually reverting toward climatology for the remainder of the 30-day forecast period). At 1/32° resolution, forecast SSH generally compares better with unassimilated observations and the anomaly correlation of the forecast SSH exceeds that from persistence by a larger amount than found in the 1/16° system.     

Atmospheric CO2 measurements and error analysis on seasonal air–sea CO2 fluxes in the Bay of Biscay

Atmospheric molar fraction of CO₂ (xCO₂^atm) measurements obtained on board of ships of opportunity are used to parameterize the seasonal cycle of atmospheric xCO₂ (xCO₂^atm) in three regions of the eastern North Atlantic (Galician and French offshore and Bay of Biscay). Three selection criteria are established to eliminate spurious values and identify xCO₂^atm data representative of atmospheric background values. The filtered data set is fitted to seasonal curve, consisting of an annual trend plus a seasonal cycle. Although the fitted curves are consistent with the seasonal evolution of xCO₂^atm data series from land meteorological stations, only ship-board measurements can report the presence of winter xCO₂^atm minimum on Bay of Biscay. Weekly air–sea CO₂ flux differences (mmol C·m^− 2 day^− 1) produced by the several options of xCO₂^atm usually used (ship-board measurements, data from land meteorological stations and annually averaged values) were calculated in Bay of Biscay throughout 2003. Flux error using fitted seasonal curve relative to on board measurements was minimal, whereas land stations and annual means yielded random (− 0.2 ± 0.3 mmol C·m^− 2·day^− 1) and systematic (− 0.1 ± 0.4 mmol C·m^− 2 day^− 1), respectively. The effect of different available sources of sea level pressure, wind speed and transfer velocity were also evaluated. Wind speed and transfer velocity parameters are found as the most critical choice in the estimate of CO₂ fluxes reaching a flux uncertainty of 7 mmol C·m^− 2·day^− 1 during springtime. The atmospheric pressure shows a notable relative effect during summertime although its influence is quantitatively slight on annual scale (0.3 ± 0.2 mmol C·m^− 2·day^− 1). All results confirms the role of the Bay of Biscay as CO₂ sink for the 2003 with an annual mean CO₂ flux around − 5 ± 5 mmol C m^− 2 day^− 1.     

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.     

Changes in the shelf macrobenthic community over large temporal and spatial scales in the Bohai Sea, China

Over the past 20 years, the Bohai Sea has been subjected to a considerable human impact through over-fishing and pollution. Together with the influence of the Yellow River cut-off, the ecosystem experienced a dramatic change. In order to integrate available information to detect any change in macrobenthic community structure and diversity over space and time, data collected during the 1980s and the 1990s from 3 regions of the Bohai Sea (Laizhou Bay, 16 stations, 37–38°N, 119–120.5°E; central Bohai Sea, 25 stations, 38–39°N, 119–121°E; eastern Bohai Bay, 12 stations, 38–39°N, 118.5–119°E) were reanalyzed in a comparative way by means of a variety of statistical techniques. A considerable change in community structure between the 1980s and the 1990s and over the geographical regions at both the species and family level were revealed. After 10 years, there was a considerable increase in abundance of small polychaetes, bivalves and crustaceans but decreased number of echinoderms. Once abundant in Laizhou Bay in the 1980s, a large echinoderm Echinocardium cordatum and a small mussel Musculista senhousia almost disappeared from the surveying area in the 1990s. Coupled with the increased abundance was the increased species richness in general whereas evenness was getting lower in central Bohai Sea and Bohai Bay but increased in Laizhou Bay. K-dominance plot showed the same trend as evenness J′. After 10 years, the macrobenthic diversity in the Bohai Sea as a whole was slightly reduced and a diversity ranking of central Bohai Sea > Laizhou Bay > eastern Bohai Bay over space was also suggested. Sediment granulometry and organic content were the two major agents behind the observed changes.     

Remotely-sensed chlorophyll a observations of the northern Red Sea indicate seasonal variability and influence of coastal reefs

The biological dynamics of the open northern Red Sea (21.5°–27.5° N, 33.5°–40° E) have not been studied extensively, due in part to both the inaccessibility of this desert region and political considerations. Remotely-sensed chlorophyll a data therefore provide a framework to investigate the primary patterns of biological activity in this oceanic basin. Monthly chlorophyll a data from the 8-year Sea-viewing Wide Field-of-View sensor (SeaWiFS) mission, and data from the Moderate Resolution Imaging Spectroradiometer (MODIS), were analyzed with the Goddard Earth Sciences Data and Information Services Center (GES DISC) online data analysis system “Giovanni”. The data indicate that despite the normal low chlorophyll a concentrations (0.1–0.2 mg m^− 3) in these oligotrophic waters, there is a characteristic seasonal bloom in March–April in the northernmost open Red Sea (24° to 27.5° N) concurrent with minimum sea surface temperature. The location of the highest chlorophyll concentrations is consistent with a linear box model [Eshel, G., and Naik, N.H., 1997. Climatological coastal jet collision, intermediate water formation, and the general circulation of the Red Sea. J. Phys. Oceanogr. 27(7), 1233–1257.] of Red Sea circulation. Two years in the data set exhibited a different seasonal cycle consisting of a relatively weak northern spring bloom and elevated chlorophyll concentrations to the south (21.5° to 24° N).The data also indicate that large coral reef complexes may be sources of either nutrients or chlorophyll-rich detritus and sediment, enhancing chlorophyll a concentration in waters adjacent to the reefs.     

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.     

Protist entrapment in newly formed sea ice in the Coastal Arctic Ocean

Protist abundance and taxonomic composition were determined in four development stages of newly formed sea ice (new ice, nilas, young ice and thin first-year ice) and in the underlying surface waters of the Canadian Beaufort Sea from 30 September to 19 November 2003. Pico- and nanoalgae were counted by flow cytometry whereas photosynthetic and heterotrophic protists ≥ 4 µm were identified and counted by inverted microscopy. Protists were always present in sea ice and surface water samples throughout the study period. The most abundant protists in sea ice and surface waters were cells < 4 µm. They were less abundant in sea ice (418–3051 × 10³ cells L^− 1) than in surface waters (1393–5373 × 10³ cells L^− 1). In contrast, larger protists (≥ 4 µm) were more abundant in sea ice (59–821 × 10³ cells L^− 1) than in surface waters (22–256 × 10³ cells L^− 1). These results suggest a selective incorporation of larger cells into sea ice. The ≥ 4 µm protist assemblage was composed of a total number of 73 taxa, including 12 centric diatom species, 7 pennate diatoms, 11 dinoflagellates and 16 flagellates. The taxonomic composition in the early stage of ice formation (i.e., new ice) was very similar to that observed in surface waters and was composed of a mixed population of nanoflagellates (Prasinophyceae and Prymnesiophyceae), diatoms (mainly Chaetoceros species) and dinoflagellates. In older stages of sea ice (i.e., young ice and thin first-year ice), the taxonomic composition became markedly different from that of the surface waters. These older ice samples contained relatively fewer Prasinophyceae and more unidentified nanoflagellates than the younger ice. Diatom resting spores and dinoflagellate cysts were generally more abundant in sea ice than in surface waters. However, further studies are needed to determine the importance of this winter survival strategy in Arctic sea ice. This study clearly shows the selective incorporation of large cells (≥ 4 µm) in newly formed sea ice and the change in the taxonomic composition of protists between sea ice and surface waters as the fall season progresses.     

The Gulf Stream, rings and North Atlantic eddy structures from remote sensing (Altimeter and SeaWiFS)

Seasonal SeaWiFS chlorophyll a concentrations cycles and annual changes of altimeter Sea Level Anomaly are derived for the subtropical North Atlantic near  35°N and along a Gulf Stream axis. Spatial structure of SeaWiFS, is defined in terms of deviations from a local seasonal cycle and examined in relation to altimeter eddy structure. In the subtropical region near 35°N, SeaWiFS structure is evident during the spring bloom period with a scale of  430 km, or about twice the eddy scale. A Gulf Stream axis has been selected as a region where the Sea Level Anomaly variance is a maximum. Eddy propagation speeds and scales are examined. Cold-core (cyclonic) rings correspond to areas of high SeaWiFS chlorophyll a. Warm-core (anticyclonic) rings relate to areas of low chlorophyll concentration. Both SeaWiFS structure and eddy structure have a spatial scale of  450 km or twice the ring scale along the Gulf Stream axis. SeaWiFS chlorophyll anomalies and Altimeter Sea Level Anomaly structure have an overall negative correlation coefficient of r = − 0.34. Swirl currents between eddies redistribute surface chlorophyll concentrations and can spatially bias maximum and minimum concentration levels off eddy centre.     

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

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

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.     

Appendicularian ecophysiology I: Food concentration dependent clearance rate, assimilation efficiency, growth and reproduction of Oikopleura dioica

Three aspects of the appendicularian O. dioica's ecophysiology were measured here: 1) morphological parameters over a wide range of appendicularian sizes, including mature animals in order to document the morphological characteristics inducing reproduction; 2) clearance rate and assimilation efficiency using feeding incubations with different algal concentrations and 3) the effect of food concentration on growth, mortality and reproduction.The relationship between the body carbon weight and the clearance rate follows a power function, with an exponent of 0.91 (± 0.07). The rate of particles retention increases with the food concentration following a Michaelis–Menten relationship (half-saturation constant = 151 ± 22 µg C l^− 1, maximum clearance rate = 12 ± 1 µg C µg C^− 1 d^− 1). The carbon assimilation efficiency decreases with the increasing food concentration. As a result, appendicularian growth which is limited in concentrations lower than 50 µg C l^− 1 is saturated above 100 µg C l^− 1.In immature animals the gonad represents less than 30% of the body volume whereas in mature individuals, its volume varies between 50% and 87% (mean 63%) suggesting that gonad/total volume ratio can be used as indicator of the maturation stages. The gonad weight in mature animals represents 70.3 (± 4.6)% of the total body carbon weight. Two major maturity stages can explain the changes in energy allocation: i) the somatic growth, when less energy is invested in gonad growth when compared to the rest of the body and ii) the maturation phase where most of the assimilated matter is invested in gonad maturation. This process is rapid, lasting only few hours. For this reason we measured completely mature organisms that are generally not measured during the experimental work with appendicularians. In food-limited conditions, the gonad maturation process starts with smaller individuals and ends with smaller reproductive animals having the same gonad to total volume ratio than in unlimited food conditions. The results obtained in this study were used to model the life cycle of O. dioica (see Lombard, F., Sciandra, A. and Gorsky, G., 2009-this volume. Appendicularian ecophysiology. II. Modeling nutrition, metabolism, growth and reproduction of the appendicularian Oikopleura dioica.).