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.     

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.     

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.     

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.     

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.     

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.     

Phytoplankton patch patterns: Seascape anatomy in a turbulent ocean

Marine phytoplankton experience competition, predation, infection and aggregation occurring across distances of micrometres to centimetres. However, the consequences of these interactions influence global processes, such as climate and fisheries productivity. There is a long-standing default assumption that these global processes cannot be traced to plankton distributions and interactions below a few metres because of the homogenising effect of turbulence [Hutchinson, G.E., 1961. The paradox of the plankton. Am. Nat. 95, 137–146.; Siegel, D.A., 1998. Resource competition in a discrete environment: Why are plankton distributions paradoxical? Limnol. Oceanogr. 43, 1133–1146.]. We show that, in active turbulence, phytoplankton patches, on the order of 10 cm, have repeatable asymmetry and regular spacing over distances of centimetres to tens of metres. The regularity and hierarchical nature of the patches in mixed ocean water means that phytoplankton are distributed in a dynamic, but definite seascape topography, where groups of patches coalesce between intermittent turbulent eddies. These patches may link large scale processes and microscale interactions, acting as fundamental components of marine ecosystems that influence grazing efficiency, taxonomic diversity, and the initiation of aggregation and subsequent carbon flux.     

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.     

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.     

Comparative measurements on the flow structure of a marine propeller wake between an open free surface and closed surface flows

In order to investigate the effects of a free surface on the wake behind a rotating propeller, experiments were carried out in a circulating water channel for two cases: one with an open free surface and one with a closed free surface covered by a rigid plate. Four hundred instantaneous velocity fields were measured using a two-frame particle image velocimetry (PIV) technique at four different blade phases. These were ensemble-averaged to investigate the time-averaged flow structure in the near-wake region. For a surface ship, the flow behind the propeller is influenced by the hull wake and the free surface. The phase-averaged mean velocity fields show the potential wake and the viscous wake formed by the boundary layers developed on the blade surfaces. The interaction between the bilge vortices and the incoming flow along the ship’s hull deforms the wake structure. Tip vortices are generated periodically, and the slipstream contraction occurs in the near-wake region. The free surface was found to affect the axial velocity component and vortex structure behind the propeller. As the flow goes downstream, the tip and trailing vortices dissipate due to turbulent diffusion and active mixing with adjacent vortices.     

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.     

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.     

Northern Adriatic Sea surface circulation and temperature/pigment fields in September and October 1997

Satellite images of surface chlorophyll-a concentration measured by the sea-viewing wide field-of-view sensor (SeaWiFS) and of sea surface temperature derived from advanced very high resolution radiometer (AVHRR) measurements, combined with in-situ drifter measurements of surface currents, and ancillary wind, Po River discharge and surface salinity data, are used to describe the surface dynamics in the northern Adriatic during the period September–October 1997.The satellite observations revealed very complex mesoscale dynamics, with time scales of a day or two and length scales of about 10 km, including the meandering and instability of basin-scale currents (e.g., the western coastal layer), jets/filaments and eddies. In addition, the two typical patterns of the Po River plume are observed and qualitatively explained in terms of wind forcing. A basin-wide double gyre pattern spreads the rich runoff water across most of the northern Adriatic from mid-September to early October, following Bora wind events and under stratified sea conditions. In contrast, in late October the Po plume is confined to the coast due to weaker winds and de-stratified conditions. This variability in the Po River plume extension is also confirmed by in-situ salinity measurements.     

Study of the seasonal cycle of the biogeochemical processes in the Ligurian Sea using a 1D interdisciplinary model

A one-dimensional coupled physical–biogeochemical model has been built to study the pelagic food web of the Ligurian Sea (NW Mediterranean Sea). The physical model is the turbulent closure model (version 1D) developed at the GeoHydrodynamics and Environmental Laboratory (GHER) of the University of Liège. The ecosystem model contains 19 state variables describing the carbon and nitrogen cycles of the pelagic food web. Phytoplankton and zooplankton are both divided in three size-based compartments and the model includes an explicit representation of the microbial loop including bacteria, dissolved organic matter, nano-, and microzooplankton. The internal carbon/nitrogen ratio is assumed variable for phytoplankton and detritus, and constant for zooplankton and bacteria. Silicate is considered as a potential limiting nutrient of phytoplankton's growth. The aggregation model described by Kriest and Evans in (Proc. Ind. Acad. Sci., Earth Planet. Sci. 109 (4) (2000) 453) is used to evaluate the sinking rate of particulate detritus. The model is forced at the air–sea interface by meteorological data coming from the “Côte d'Azur” Meteorological Buoy. The dynamics of atmospheric fluxes in the Mediterranean Sea (DYFAMED) time-series data obtained during the year 2000 are used to calibrate and validate the biological model. The comparison of model results within in situ DYFAMED data shows that although some processes are not represented by the model, such as horizontal and vertical advections, model results are overall in agreement with observations and differences observed can be explained with environmental conditions.     

Hydrographic cruises off northwest Africa: the Canary Current and the Cape Ghir region

We present hydrographic data for several sections located along the African coastline and off Cape Ghir, carried out at times of weak surface winds (October 1995 and September 1997). The main sections are near the continental slope, at mean distances between 100 and 150 km from the coastline. North of Cape Ghir (31°N) the geostrophic transport (referenced to 650 m) of North Atlantic Central Water through these sections is 3.7 and 2.0 Sv for 1995 and 1997, respectively. This confirms that a major fraction of the water transport by the Canary Current flows east, into the continental slope off northwest Africa, at latitudes above Cape Ghir. Most of this flow continues south past Cape Ghir, along the coast and probably through the eastern passages of the Canary Archipelago. A significant fraction, however, may escape through surface Ekman transport (0.3–0.5 Sv during the early fall season) and by offshore flow at Cape Ghir (1.1 Sv in September 1997, referenced to 650 m). Despite the weak winds the Cape Ghir filament was clearly visible, characterized by localized coastal upwelling associated to a cyclonic shallow structure and cold (and fresh) waters stretching offshore as a very shallow feature (50–100 m deep). The satellite images show that the surface temperature field is highly variable, in rapid response to the surface winds, always with a core region of relatively cold water and commonly with one or two associated eddies. Our results support the existence of two recirculation cells in the area: a horizontal one that connects the interior eastern boundary currents with the coastal region and a vertical one related to both wind-induced and filament upwelling. The data also show a salinity subsurface maximum at the root of the filament, linked to water inflow from northern latitudes, and a subsurface anticyclonic eddy over the Agadir canyon, likely related to the poleward slope undercurrent.     

Numerical experiments on the prediction of sea surface temperature anomalies in the Gulf of Mexico

The conservation of thermal energy equation applied to the mixed layer of the ocean, has been used to predict the sea surface temperature anomalies (SSTA) and the month-to-month changes in the Gulf of Mexico. The model includes the horizontal transport of heat by mean ocean currents and by turbulent eddies, as well as the heating by short and long wave radiation, evaporation and sensible heat given off to the atmosphere. A comparative study is carried out on the relative importance of the heating and transport terms. An objective verification of the skill of the predictions is presented for each season and for the whole period from March 1986 to February 1987. The predictions using only the heating terms have some skill over the control predictions (persistence and return to normal). The skill is substantially increased when the horizontal transport of heat by turbulent mixing is included in the model. The incorporation in the model of the Ekman wind drift current anomalies computed from the anomalous surface geostrophic wind improves appreciably the skill of the predictions in winter and fall.The mixed layer depth computed using the Kraus and Turner theory with dissipation, shows that the depths in summer and fall are shallower than in spring and winter. The effect of the shallow mixed layer depth in the model becomes apparent in summer and fall, improving the skill of the predictions in these seasons, with respect to the skill obtained using a constant mixed layer depth of 60 m.The incorporation in the model of the cooling in the mixed layer by turbulent entrainment of colder water from the thermocline, does not improve in an appreciable way the average skill of the predictions.     

Spatial distribution of cadmium over a cyclonic eddy in the southern East China Sea

Eddy systems are a unique ecosystem, usually having high biological masses and primary production in the sea. In this study, both particulate and dissolved phases of cadmium in the water column of 15 stations over a cyclonic eddy in the southern East China Sea were determined to obtain their spatial distributions. This allows us to confirm that cyclonic eddy systems play the role of a Cd pump in the sea.Results showed that particulate Cd (PCd) and dissolved Cd (DCd) concentrations in water varied greatly, ranging over two orders and one order of magnitude, respectively. Large spatial variability was found not only for PCd but also for DCd in the upper water, apparently due to the effects of the cyclonic eddy system on the Cd distributions over the study area. DCd accounts for about 99% of the total Cd. For the surface water, DCd concentration at the eddy center was about five times the average of the water surrounding the eddy center. The depth distributions of DCd exhibited a typical surface depletion and a subsequent increase with depth; however, the PCd distribution showed the opposite, i.e. a surface maximum and a subsequent decline with depth. In general, the DCd maxima were found at depths of 600–1000 m, agreeing well with the literature. It reflected the internal biogeochemical cycling of Cd in the water column, which was driven by the utilization of Cd by plankton in the euphotic zone and by the regeneration of Cd at depth. In addition, a remarkably high DCd concentration existed in near-bottom water around the cyclonic eddy center.The horizontal distributions of both PCd and DCd in the upper water shared a common feature with elevated concentrations centering around the eddy center and a decline in concentration with distance from the eddy center. This shows that the cyclonic eddy could bring up the Cd-rich deep water to the surface water around the eddy center and could then expand toward eddy's vicinity via advection–diffusion. It is thus likely that it produces ample biological masses over the eddy system. Hence, this work can demonstrate that the ascending nutrient-rich water driven by the cyclonic eddies can serve as an important source not only for many nutrients but also for Cd in the sea.     

PIV velocity field measurements of flow around a KRISO 3600TEU container ship model

The main purpose of this investigation was to demonstrate a useful application of the particle image velocimetry (PIV) method to analyze the complex flow characteristics around a ship. For a sample illustration, the KRISO 3600TEU container ship model was chosen. The flow structure in the stern and near-wake region of the model has been investigated experimentally in a circulating water channel. Instantaneous velocity fields measured by the PIV velocity field measurement technique have been ensemble-averaged to give details of flow structures such as the spatial distributions of the local mean velocity, vorticity, and turbulent kinetic energy. The free-stream velocity was fixed at U _o = 0.6m/s, and the corresponding Reynolds number based on the length between perpendiculars was about 9.0 × 10⁵. The container ship model shows a complicated three-dimensional flow structure in the stern and near-wake regions. The PIV results clearly revealed the formation of large-scale bilge vortices in the stern region and their effect on the flow in the near-wake. The results shown here provide valuable information for hull form design and the validation of viscous ship flow codes and of turbulence models.