Surface CO2 measurements in the English Channel and Southern Bight of North Sea using voluntary observing ships

Ships of opportunity have been used to investigate ocean–atmosphere CO₂ fluxes in the English Channel and Southern Bight of the North Sea. Continuous underway measurements of the fugacity of seawater carbon dioxide (fCO₂^sw), chlorophyll, temperature and salinity have been performed along 26 transects during the spring and autumn periods. The spatial fCO₂^sw distribution along the Channel and Southern Bight is modulated by the photosynthetic activity, temperature changes and water mixing between inputs from the North Atlantic Ocean and riverine discharges. The seasonal variability of fCO₂^sw is assessed and discussed in terms of the biology and temperature effects, these having similar impacts. The variation of fCO₂^sw shows similar interannual patterns, with lower values in spring. The annual average of air–sea CO₂ fluxes places the English Channel as neutral area of CO₂ uptake. The spring and autumn data allow differentiating between distal and proximal continental areas. The Southern Bight shows a tendency towards net CO₂ uptake on the distal continental shelf, whereas the Scheldt and Thames Plumes show a CO₂ source behaviour on the proximal continental shelves.     

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 exchange of mercury vapour over the Gulf of Gda sk and southern Baltic Sea

The total gaseous mercury (TGM) in air over the coastal station at Hel and over the southern Baltic Sea was measured during the summer and winter conditions. Recorded 30-min resolution TGM data showed both higher concentrations and variability during the summer compared to the winter conditions. The summer TGM data ranged from 1.1 to 7.5 ng m⁻³, while the winter data ranged from 0.8 to 4.4 ng m⁻³. The TGM content in air over the southern Baltic Sea indicated that, in general, during the summer conditions, the sea-to-air transport of gaseous mercury dominated, while during the winter season, a tendency of gaseous mercury to sink into the water has been found. The evidences of enhanced water-to-air transfer of mercury vapour were noted, in particular, over the shallow waters of the Gulf of Gda sk under the strong water-to-air temperature gradients. Obtained results indicate that under such conditions, the coastal waters could act as a significant source of mercury vapour that may contribute to the overall budget of atmospheric mercury over the Baltic proper.     

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.     

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.     

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.     

Coupled renewal model of ocean viscous sublayer, thermal skin effect and interfacial gas transfer velocity

A previously developed renewal model included parameterizations for the thermal skin effect and interfacial gas transfer velocity. The more readily available cool skin data were used for an adjustment of the gas transfer parameterization. In this work, the renewal concept is extended to include the velocity difference across the viscous sublayer and to account for the stage of surface wave development. As a result, the empirical coefficients that enter the renewal model have been specified more accurately using laboratory data on the surface wind drift current. In addition, the coefficient linking the cool skin and gas transfer parameterization formulas has been determined from the probability distribution function for renewal events. A comparison of the upgraded renewal model with the thermal skin data collected during the COARE and more recent field programs and with gas transfer data collected during GasEx-01 experiment suggests that the renewal model can be a useful tool for producing a physically based parameterization for the interfacial CO₂ transfer velocity. Model uncertainties associated with the effect of surface films are discussed.     

Dynamics of inorganic nutrient species in the Bohai seawaters

Within the frame of a Sino-German Joint Research Program, two cruises of “R/V Dong Fang Hong 2” were carried out in September–October 1998 and April–May 1999, respectively, to understand the dynamics of nutrients in the Bohai. Nutrient species (NO₃⁻, NO₂⁻, NH₄⁺, PO₄³⁻ and SiO₃²⁻) are determined colorimetrically on board for five anchor and 30 grid stations. In situ incubation experiments are performed to determine planktonic nutrient uptake and benthic exchange flux. Nutrient concentrations display short-term variability and seasonal change in the Bohai, with higher levels in shallow coastal waters than in the Central Bohai. The influence of riverine discharge on nutrient levels can be seen from salinity isopleths, nutrient distribution and species ratios. Near-bottom (nb) waters have similar nutrient concentrations as to the surface waters in the Central Bohai, whereas stratification takes place in the Bohai Strait and North Yellow Sea. In situ incubation experiments provide evidence that the uptake ratio (i.e. N, P) by phytoplankton is proportional to the ratios among nutrient species in ambient waters. Based on the data of this study and previously publications, a preliminary estimate of nutrient budgets via riverine input and atmospheric deposition is established. The results indicate that atmospheric deposition gains importance over rivers in delivering nutrients into the Bohai and sustain the new production, following recent decrease in riverine inflow caused by drought periods in North China and damming practices. A historical review of nutrient data indicates that concentrations of nitrogen increase and phosphorus and silica decrease in the Central Bohai over last 40 years. This potentially has an important influence on the health of ecosystem in Bohai (e.g. food web and community structure), though further study is needed to examine the scenario in more detail.     

Evolving the subspace of the three-dimensional multiscale ocean variability: Massachusetts Bay

A data and dynamics driven approach to estimate, decompose, organize and analyze the evolving three-dimensional variability of ocean fields is outlined. Variability refers here to the statistics of the differences between ocean states and a reference state. In general, these statistics evolve in time and space. For a first endeavor, the variability subspace defined by the dominant eigendecomposition of a normalized form of the variability covariance is evolved. A multiscale methodology for its initialization and forecast is outlined. It combines data and primitive equation dynamics within a Monte-Carlo approach.The methodology is applied to part of a multidisciplinary experiment that occurred in Massachusetts Bay in late summer and early fall of 1998. For a 4-day time period, the three-dimensional and multivariate properties of the variability standard deviations and dominant eigenvectors are studied. Two variability patterns are discussed in detail. One relates to a displacement of the Gulf of Maine coastal current offshore from Cape Ann, with the creation of adjacent mesoscale recirculation cells. The other relates to a Bay-wide coastal upwelling mode from Barnstable Harbor to Gloucester in response to strong southerly winds. Snapshots and tendencies of physical fields and trajectories of simulated Lagrangian drifters are employed to diagnose and illustrate the use of the dominant variability covariance. The variability subspace is shown to guide the dynamical analysis of the physical fields. For the stratified conditions, it is found that strong wind events can alter the structures of the buoyancy flow and that circulation features are more variable than previously described, on multiple scales. In several locations, the factors estimated to be important include some or all of the atmospheric and surface pressure forcings, and associated Ekman transports and downwelling/upwelling processes, the Coriolis force, the pressure force, inertia and mixing.     

Steady sailing performance of a hybrid-sail assisted bulk carrier

The steady sailing performance of a sail-assisted bulk carrier is investigated utilising towing-tank derived hydrodynamic derivatives and wind tunnel measured aerodynamic properties of the sails and the ship. The aerodynamic characteristics investigated include the ship hull at the fully-loaded draught, the sail–sail interaction effects for two sets of four identical hybrid-sails, and the sail–hull interaction effects for the same two sets of identical sails in the presence of the selected bulk carrier hull-form. This is in addition to lift–drag measurements of single isolated sails of each shape. The form of the two sets of soft sails was rectangular and triangular. This paper is concerned with assessing the benefits of a sail-assisted ship operation, and hence a steady-state rather than complete time-domain integrations of the governing equations are reported. The results of the completed analysis suggest that the benefits of the derived sail generated driving force are greater than the overhead of equipping the ship with a selected system of hybrid-sails. Sail-assisted ships could represent an important contribution to an improving global environment by reducing the demands for a driving force through the propeller.     

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.     

The shoaling of nutrient-enriched subsurface waters as a mechanism to sustain primary productivity off Central Baja California during El Niño winters

Using CalCOFI data for coastal shallow stations (above 100 m depth), higher than expected nitrate concentrations were detected in near-surface high-temperature waters off of Central Baja California during some El Niño winters. Though recent data are not available for Central Baja California, past El Niño data, though limited, showed nitrate concentrations above 16 μM at temperatures above 16 °C, and nitrate concentrations between 1 and 2 μM at 19 °C, while the previously established relationship of temperature and nitrate for California Current waters predicts nitrate depletion above 14 or 15 °C. The anomalous, high temperature–high nitrate enrichment events documented in Central Baja California were detected as shallow as 9 m and as deep as 73 m, were associated with low-oxygen (between 2 and 4 ml/l) and high-salinity (between 33.8 and 34.3 psu) waters, and occurred during the winter months of an El Niño year. Using recent data for San Diego, CA, similar but weaker enrichment events were detected for the El Niño winter of 1997–1998. The periodic shoaling of a subsurface subtropical water mass of high temperature, high salinity, low oxygen and high nutrients during some El Niño winters is proposed to cause periodic enrichment and to maintain productivity during warming events in this area. Enrichment events were not detected off Ensenada, in Northern Baja California, possibly due to the amplification of the onshore flow during El Niño there, or due to the Ensenada front. The proposed mechanism of periodic enrichment of nutrient-depleted surface waters during some El Niño winters by subsurface waters from the California Undercurrent may explain the following: (1) survival of giant kelp forests at their southern limit in Central Baja California documented during past El Niño events in warm waters, (2) the rapid recovery and high carrying capacity of giant kelp documented after the mass disappearance during El Niño 1997–1998, and (3) the increase in the extent of mesotrophic chlorophyll detected in the area during the 1997–1998 and 1982–1983 El Niño events.     

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.     

Quantitative estimation of the influence of surface thermal fronts over chlorophyll concentration at the Patagonian shelf

Eighteen-year (1985–2002) mean monthly SST Pathfinder data with 9 km spatial resolution have been used to estimate surface gradients by finite differences. Then the seasonal climatological means have been calculated from the intensity of these gradients, and surface thermal fronts present in the Patagonian Continental Shelf (PCS) have been located. Moreover, 6 years (1998–2003) of SeaWiFS data with approximately 4 km spatial resolution have been used to estimate monthly composite images of surface chlorophyll concentration, after which seasonal climatological means distributions have been generated. Both seasonal distributions have been analyzed together and by combining the knowledge of oceanographic processes and phytoplankton responses to light and nutrient availability, regions where the presence of a thermal front affects photosynthetic activity have been identified. Subjective criteria have been applied to define eighteen areas where phytoplankton biomass is influenced by the presence of a thermal front. In these areas, the surface chlorophyll (spatial mean and total), its relationship with the surface chlorophyll of the whole region, and the seasonal evolution of this relationship have been calculated. All frontal areas cover less than 15% of the total surface, but they contribute with over 23% of the phytoplankton annual mean biomass. Considered as a group, during summer they show high chlorophyll values very similar to those in spring. During the cold period, when the water column is vertically mixed in practically the whole of PCS, the influence of physical fronts over the biological production is minimum. The frontal zone image remains clearly defined during summer, when approximately 85% of the area will have a determined mean chlorophyll concentration, while the other 15% has a 2.45 times larger value. While three pattern trends have been identified in the frontal areas, only two of them condition the pattern of the group, due to their horizontal extension.     

Physical–biological coupling in spore dispersal of kelp forest macroalgae

The physical–biological linkages controlling the dispersal of spores produced by macroalgae that reside in kelp forests are complicated and laced with feedbacks. Here we discuss the fundamental elements of these interactions. Biological considerations include spore swimming and sinking speeds, their periods of viability in the plankton, and the height of spore release above the seafloor, which together determine the durations over which spores can be swept by horizontal currents before they contact the seafloor. Morphologies and material properties of canopy forming kelps may also influence the drag exerted on passing waters by the kelps, the plants' ability to persist in the face of rapid flows, and thereby the degree to which impinging currents are redirected around, or slowed within, kelp forests. Macroalgal life histories, and the size of spore sources as controlled by the dimensions of kelp forests and the density and fecundity of individuals within them, influence effective dispersal distances as well. Physical considerations encompass the mean speed, direction, and timescales of variability of currents relative to spore suspension times, the interaction of surface gravity waves with currents in producing turbulence in the benthic boundary layer, wind-driven surface mixing, water stratification, and shoreline bathymetry and substratum roughness, all of which can affect the interplay of vertical and horizontal transport of macroalgal spores. Intricate within-forest processes may induce attenuation of current speeds and consequent reductions in seabed shear, along with simultaneous production of small-scale turbulence in kelp wakes. Slower mean currents and smaller eddy scales in turn may attenuate vertical mixing within forests, thus extending spore suspension times. Further complexities likely arise due to changes in the relative rates of horizontal and vertical dispersion, modifications to the overall profiles of vertical mixing, and the creation of fine-scale secondary flows around kelp individuals and substratum features. Under conditions of more rapid currents, submergence of the surface canopy and the establishment of skimming flows at the canopy–fluid interface may introduce additional coherent flow structures that alter rates of fluid exchange to and from the forest. Many of these coupled physical–biological processes are just beginning to be examined in a rigorous fashion in kelp forests, but their potential importance is clear.     

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.     

Patterns in diversity and community structure of epipelagic copepods from the Brazil–Malvinas Confluence area, south-western Atlantic

Diversity and structure of epipelagic copepod communities were investigated using 70 zooplankton samples collected from the top 50 m of the Brazil–Malvinas Confluence area between 35° and 55°S. Biogeographic patterns were investigated using multivariate analysis. Biodiversity patterns were examined using different univariate indexes. Representatives of 35 species of copepods from 23 genera and 13 families were found. Two zones were delimited based on their copepod composition, one in the subtropical waters of the Brazil current, and the other in the subantarctic waters of the Malvinas current. Both environments displayed contrasting patterns of biodiversity depending upon which element is measured (regional species richness as well as mean point species richness were significantly higher in the subtropical group of stations, whereas taxonomic distinctness was significantly higher in the subantarctic community). Based on these contrasting results, we suggest the use of both kind of indexes when defining priority areas for conservation.     

Coupled transport-reaction pathways and distribution patterns between siliciclastic-carbonate sediments at the Ria de Vigo

This paper examines the linkages between the space-distribution of grain sizes and the relative percentage of the amount of mineral species that result from the mixing process of siliciclastic and carbonate sediments at the Ria de Vigo (NW of Spain).The space-distribution of minerals was initially determined, starting from a detailed mineralogical study based on XRD-Rietveld analysis of the superficial sediments. Correlations between the maps obtained for grain sizes, average fractions of either siliciclastic or carbonates, as well as for individual-minerals, were further stabilised. From this analysis, spatially organized patterns were found between carbonates and several minerals involved in the siliciclastic fraction. In particular, a coupled behaviour is observed between plagioclases and carbonates, in terms of their relative percentage amounts and the grain size distribution.In order to explain these results a conceptual model is proposed, based on the interplay between chemical processes at the seawater–sediment interface and hydrodynamical factors. This model suggests the existence of chemical control mechanisms that, by selective processes of dissolution-crystallization, constrain the mixed environment's long-term evolution, inducing the formation of self-organized sedimentary patterns.     

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.     

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

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