A numerical model of the long term flow along the Malin-Hebrides shelf

A three dimensional hydrodynamic model of the Malin-Hebrides shelf region is used to investigate the spatial variability of the wind and tidally induced residual flow in the region and the influence of flow from the Irish Sea and along the shelf edge. By this means it is possible to understand the spatial variability in the long term observed flow fields in the region and the range of driving forces producing this flow. The model uses a sigma coordinate grid in the vertical with a finer grid in the near surface and near bed shear layers. The vertical diffusion of momentum in the model is parameterised using an eddy viscosity coefficient which is derived from turbulence energy closure models. Two different turbulence models are used to compute the eddy viscosity, namely a two-equation (itq²−q²ℓ) model which has prognostic equations for both turbulence energy and mixing length and a simpler model in which the mixing length is a specified algebraic function of the water depth.The wind induced response to spatially and temporally constant orthogonal wind stresses, namely westerly and southerly winds of 1 N m⁻², are derived from the model. By using orthogonal winds and assuming linearity, then to first order the response to any wind direction can be derived. Computed flows show a uniform wind driven surface layer of magnitude about 3% of the wind speed and direction 15 ° to the right of the wind, in deep water. Currents at depth particularly in the shelf edge and near coastal region show significant spatial variability which is related to variations in bottom topography and the coastline.Calculations show that tidal residual flows are only significant in the near coastal regions where the tidal current is strong and exhibits spatial variability. Flow into the region from the Irish Sea through the North Channel although having its greatest influence in the near coastal region, does affect currents near the shelf edge region. Again the spatial variability of the flow is influenced by topographic effects.A detailed examination of wind induced current profiles together with turbulence, mixing length and viscosity, at a number of locations in the model from deep ocean to shallow near coastal, shows that both turbulence models yield comparable results, with the mixing length in the two equation model showing a similar dependence to that specified in the simpler turbulence model.Calculations clearly show that flow along the shelf edge area to the west of Ireland and from the Irish Sea entering the region, together with local wind forcing can have a major effect upon currents along the Malin-Hebrides shelf. The flow fields show significant spatial variability in the region, comparable to those deduced from long term tracer measurements. The spatial variability found in the calculations suggests that a very intense measurement programme together with inflow measurements into the area is required to understand the circulation in the region, and provide data sets suitable for a rigorous model validation.     

The wind-forced response on a buoyant coastal current: Observations of the western Gulf of Maine plume

The freshwater plume in the western Gulf of Maine is being studied as part of an interdisciplinary investigation of the physical transport of a toxic alga. A field program was conducted in the springs of 1993 and 1994 to map the spatial and temporal patterns of salinity, currents and algal toxicity. The observations suggest that the plume's cross-shore structure varies markedly as a function of fluctuations in alongshore wind forcing. Consistent with Ekman drift dynamics, upwelling favorable winds spread the plume offshore, at times widening it to over 50 km in offshore extent, while downwelling favorable winds narrow the plume width to as little as 10 km.Using a simple slab model, we find qualitative agreement between the observed variations of plume width and those predicted by Ekman theory for short time scales of integration. Near surface current meters show significant correlations between cross-shore currents and alongshore wind stress, consistent with Ekman theory. Estimates of the terms in the alongshore momentum equation calculated from moored current meter arrays also indicate a dominant Ekman balance within the plume. A significant correlation between alongshore currents and winds suggests that interfacial drag may be important, although inclusion of a Raleigh drag term does not significantly improve the alongshore momentum balance.     

Effects of surfactant on the growth of wind waves — Simultaneous observations with optical and microwave sensors

The growth of wind waves was investigated in a circulating tank over slick and “clean” surfaces. The slick surface was produced through the aging of seawater, and the “clean” surface was obtained by overflowing. Water-surface slopes and microwave backscattering were simultaneously measured with an optical slope gauge and a continuous-wave radar, respectively. Dependencies of the mean-square slope and the radar backscattered power on the wind-friction velocity over these two surfaces were found to be different. A delayed growth of ripples over the slick surface was identified from optical and microwave measurements. The suppression of ripples by surface slicks was most significant at low winds up to the wind-friction velocity of 15 cm s⁻¹. A reduction of −18 dB in the mean-square slope is found at the wind-friction velocity of 11 cm s⁻¹. The exponent of power-law dependence of radar cross section on wind-friction velocity over the slick surface was much larger than that over the “clean” surface. The radar backscattering from the slick surface was reduced by −40 dB at light winds and saturated at high winds with the wind-friction velocities smaller than 11 cm s⁻¹ and larger than 24 cm s⁻¹, respectively.     

The structure of dissipation in the western Irish Sea front

We report on an intensive campaign in the summer of 2006 to observe turbulent energy dissipation in the vicinity of a tidal mixing front which separates well mixed and seasonally stratified regimes in the western Irish Sea. The rate of turbulent dissipation ε was observed on a section across the front by a combination of vertical profiles with the FLY dissipation profiler and horizontal profiles by shear sensors mounted on an AUV (Autosub). Mean flow conditions and stratification were obtained from a bed mounted ADCP and a vertical chain of thermistors on a mooring. During an Autosub mission of 60 h, the vehicle, moving at a speed of ~ 1.2 m s^− 1, completed 10 useable frontal crossings between end points which were allowed to move with the mean flow. The results were combined with parallel measurements of the vertical profile of ε which were made using FLY for periods of up to 13 h at positions along the Autosub track. The two data sets, which show a satisfactory degree of consistency, were combined to elucidate the space–time variation of dissipation in the frontal zone. Using harmonic analysis, the spatial structure of dissipation was separated from the strong time dependent signal at the M₄ tidal frequency to yield a picture of the cross-frontal distribution of energy dissipation. A complementary picture of the frontal velocity field was obtained from a moored ADCP and estimates of the mean velocity derived from the thermal wind using the observed density distribution. which indicated the presence of a strong (0.2 m s^− 1) jet-like flow in the high gradient region of the front. Under neap tidal conditions, mean dissipation varied across the section by 3 orders of magnitude exceeding 10^− 2 W m^− 3 near the seabed in the mixed regime and decreasing to 10^− 5 W m^− 3. in the strongly stratified interior regime. The spatial pattern of dissipation is consistent in general form with the predictions of models of tidal mixing and does not reflect any strong influence by the frontal jet.     

Turbulent mixing and internal tides in Gaoping (Kaoping) Submarine Canyon, Taiwan

Turbulent overturning on scales greater than 10 m is observed near the bottom and in mid-depth layers within the Gaoping (formerly spelled Kaoping) Submarine Canyon (KPSC) in southern Taiwan. Bursts of strong turbulence coexist with bursts of strong sediment concentrations in mid-depth layers. The turbulence kinetic energy dissipation rate in some turbulence bursts exceeds 10^− 4 W kg^− 1, and the eddy diffusivity exceeds 10^− 1 m² s^− 1. Within the canyon, the depth averaged turbulence kinetic energy dissipation rate is ~ 7 × 10^− 6 W kg^− 1, and the depth averaged eddy diffusivity is ~ 10^− 2 m² s^− 1. These are more than two orders of magnitude greater than typical values in the open ocean, and are much larger than those found in the Monterey Canyon where the strong turbulent mixing has also been. The interaction of tidal currents with the complex topography in Gaoping Submarine Canyon is presumably responsible for the observed turbulent overturning via shear instability and the breaking of internal tides and internal waves at critical frequencies. Strong 1st-mode internal tides exist in KPSC. The depth averaged internal tidal energy near the canyon mouth is ~ 0.17 m² s^− 2. The depth integrated internal tidal energy flux at the mouth of the canyon is ~ 14 kW m^− 1, propagating along the axis of the canyon toward the canyon head. The internal tidal energy flux in the canyon is 3–7 times greater than that found in Monterey Canyon, presumably due to the more than 10 times larger barotropic tide in the canyon. Simple energy budget calculations conclude that internal tides alone may provide energy sufficient to explain the turbulent mixing estimated within the canyon. Further experiments are needed in order to quantify the seasonal and geographical distributions of internal tides in Gaoping Submarine Canyon and their effects on the sediment flux in the canyon.     

Westward moving waves or eddies (Storms) on the Subtropical/Azores Front near 32.5°N? Interpretation of the Eulerian currents and temperature records at moorings 155 (35.5°W) and 156 (34.4°W)

Three Argos buoy-years of Lagrangian data in westward-moving cyclonic eddies, or Storms, near 32.5°N, together with hydrographic measurements, have shown that Storms move westward at nearly 3 km day⁻¹. Water in eddies can be trapped and moved westward by advection within the eddy or by phase propagation of the eddy pattern, so we cannot say that the flow field (or Eulerian mean) is 3 km day⁻¹ westward. Two moorings (155 and 156) deployed in the Storm Corridor have provided a further 8 instrument-years of Eulerian data. The temperature and current records confirmed that two Storms a year passed each mooring over the 2-year measurement period. As expected, there is a lag of 1.3 month at mooring 155 (which is 102 km to the west of mooring 156) with respect to conditions at mooring 156. The progressive vector diagrams (PVDs) derived from the current meter records exhibit fairly regular X (east or zonal) and Y (north or meridional) displacement scales that repeat with semi-annual periodicity (SAP). The dominant current signal is the north component of the SAP, which reaches an amplitude of 18 cm s⁻¹ for the upper layer or near surface current record (242-m depth). The geostrophic north component values derived from altimetry were in good agreement with the upper layer current meter measurements. The large north component amplitude was not interpreted as evidence for Rossby Waves but rather due to the passage of nine eddies (eight complete) of alternate sign (cyclonic, anticyclonic) passing the mooring rigs during the 2-year deployment period. The Y scale shows that the near surface characteristic or mean maximum azimuthal speed is about 35 cm s⁻¹ for cyclonic eddies or Storms, and that this value is reduced to 4 cm s⁻¹ at 1400-m depth. The residual or mean Eulerian currents range from 8 cm s⁻¹ for the upper layer currents to 1 cm s⁻¹ for the deeper currents at 1400-m depth and are predominantly westward. Simple theoretical considerations and idealised numerical simulations show that the mean westward Eulerian current depends markedly on whether the eddy centres pass to the north or south of the rigs. This raises the question as to what do we mean by Eulerian residual currents, even for relatively long records (2 years). It is shown that the strong near surface westward current (6 km day⁻¹) measured at mooring 155 is largely due to a westward-moving eddy field with variable centre offsets. The magnitude of the near surface east–west component of flow was estimated as eastward at 2 cm s⁻¹. The north–south component of mean flow was southward at 2 cm s⁻¹. The deeper records gave a weak westward flow of 1 cm s⁻¹ but did not show a significant southward component for the mean Eulerian flow field. 7.4 float-years of Lagrangian ALACE data in the Subtropical Front region near 740 dbar gave mean east–west flows that were <0.5 cm s⁻¹. Overall, it is shown that the eddy structures propagate westward mainly by phase propagation (i.e. a westward-moving pattern with no westward advection for the current meter to measure), though plane Rossby Wave dynamics appeared inappropriate. Theoretical and modeling considerations show that a speed of 3-km day⁻¹ westward is too large a value for the self-advection of eddies due to the beta effect.     

Interannual variability of the shelf-slope front position between 75° and 50° W

The shelf-slope front (SSF) is a continuous shelf-break front running from the Tail of the Grand Banks to Cape Hatteras, North Carolina, separating colder and less-saline continental shelf waters from warmer and more saline slope waters. Time series containing mean monthly SSF positions were produced along each of 26 longitude lines between 75° and 50°W by workers located at Bedford Institute of Oceanography by digitizing individual frontal charts and computing mean monthly latitudinal positions over a 29-year (1973–2001) period. After removing seasonal variability at each longitude, interannual variability (IAV) of the SSF position at each longitude was computed as the annual mean of all monthly SSF position anomalies for each year over the 29-year period. Despite some missing data, a longitude-time plot reveals alternating bands of offshore (late-1970s, late-1980s, late-1990s) and onshore (early-1980s, early-1990s, early-2000s) annual mean SSF anomaly values, exhibiting a period of approximately 10 years. Annual mean SSF anomaly amplitudes are largest in the east, with maxima of O (± 100 km) located east of 60° W for years when data are available. Empirical orthogonal function (EOF) modes 1–4 (accounting for > 90% of the variance) form a set of basis functions that describe the SSF anomaly data and allow reconstruction of the entire data set since missing data are relatively few (14%). A complex empirical orthogonal function (CEOF) analysis using the “reconstructed” data reveals a wavelength scale of approximately 20° of longitude, a distance nearly equal to the entire study domain, along with steady, westward phase propagation of SSF anomalies over approximately the same distance. Speed calculations for the westward-propagating features yield a value of approximately 1.2 to 2.4 cm s^− 1 (1 to 2 km d^− 1), with annual mean SSF anomalies thus requiring about 4 years to propagate from the Tail of the Grand Banks in the east to Cape Hatteras, North Carolina, in the west. This propagation speed and the timing of the SSF positional anomalies at the Tail of the Grand Banks for the 29-year study period are in agreement with speeds computed for the propagation of quasi-decadal salinity anomalies through the Labrador Sea and the time of their arrival at the Tail of the Grand Banks. The small westward SSF anomaly propagation speed is an order of magnitude smaller than the associated currents, in agreement with a highly damped flow-through system originating from both Davis Strait and the West Greenland Current as discussed by other workers. Observations from both southern and northern portions of the study domain, within both continental shelf and slope waters, show that interannual changes in the volume of shelf water along with shelf water bulk properties exhibit a strong relationship with IAV of the SSF position over long time periods.     

Krill behaviour as recorded by acoustic doppler current profilers in the Gullmarsfjord

Vertical distribution of sound scattering layers were observed using bottom deployed acoustic doppler current profilers (ADCP) during early spring of 1996 and autumn of 1997 in the Gullmarsfjord on the Swedish west coast. Variations in relative backscatter were interpreted in relation to horizontal water velocities, oxygen saturation as well as differences in the light, salinity and temperature regimes. Net catches revealed that much of the backscatter below 20-m depth was associated with the presence of krill, principally Meganyctiphanes norvegica.Horizontal currents seemed to influence the migration and distribution of krill, which showed weak vertical migration patterns with low abundance during periods of strong intermediate in- and outflows, while during periods with weaker currents, a more regular diel migration occurred. Horizontal water velocities >5 cm s⁻¹ seemed to have the potential to decrease the peak in the backscatter profile. Mean vertical migration rates of krill was 1 cm s⁻¹, while maximum vertical migration rates were estimated to be 2.5–3 cm s⁻¹. The range of the vertical migration was different in 1997 due to severe oxygen deficiency in the bottom water, which prevented the krill from descending >80 m. The commencement of vertical migration correlated closely to the seasonal light conditions. The descent was immediately triggered by sunrise, while ascent occurred with a delay of about 1 h at sunset.     

Suspended sediment and erosion dynamics in Kugmallit Bay and Beaufort Sea during ice-free conditions

The Mackenzie River is the largest river on the North American side of the Arctic and its huge freshwater and sediment load impacts the Canadian Beaufort Shelf. Huge quantities of sediment and associated organic carbon are transported in the Mackenzie plume into the interior of the Arctic Ocean mainly during the freshet (May to September). Changing climate scenarios portend increased coastal erosion and resuspension that lead to altered river-shelf-slope particle budgets. We measured sedimentation rates, suspended particulate matter (SPM), particle size and settling rates during ice-free conditions in Kugmallit Bay (3–5 m depth). Additionally, measurements of erosion rate, critical shear stress, particle size distribution and resuspension threshold of bottom sediments were examined at four regionally contrasting sites (33–523 m depth) on the Canadian Beaufort Shelf using a new method for assessing sediment erosion. Wind induced resuspension was evidenced by a strong relationship between SPM and wind speed in Kugmallit Bay. Deployment of sediment traps showed decreasing sedimentation rates at sites along an inshore–offshore transect ranging from 5400 to 3700 g m^− 2 day^− 1. Particle settling rates and size distributions measured using a Perspex settling chamber showed strong relationships between equivalent spherical diameter (ESD) and particle settling rates (r² = 0.91). Mean settling rates were 0.72 cm s^− 1 with corresponding ESD values of 0.9 mm. Undisturbed sediment cores were exposed to shear stress in an attempt to compare differences in sediment stability across the shelf during September to October 2003. Shear was generated by vertically oscillating a perforated disc at controlled frequencies corresponding to calibrated shear velocity using a piston grid erosion device. Critical (Type I) erosion thresholds (u) varied between 1.1 and 1.3 cm s^− 1 with no obvious differences in location. Sediments at the deepest site Amundsen Gulf displayed the highest erosion rates (22–54 g m^− 2 min^− 1) with resuspended particle sizes ranging from 100 to 930 µm for all sites. There was no indication of biotic influence on sediment stability, although our cores did not display a fluff layer of unconsolidated sediment. Concurrent studies in the delta and shelf region suggest the importance of a nepheloid layer which transports suspended particles to the slope. Continuous cycles of resuspension, deposition, and horizontal advection may intensify with reduction of sea ice in this region. Our measurements coupled with studies of circulation and cross-shelf exchange allow parameterization and modeling of particle dynamics and carbon fluxes under various climate change scenarios.     

Estimation of the surface and mid-depth currents from Argo floats in the Pacific and error analysis

With the large deployment, the Array for Real-time Geostrophic Oceanography program has great potential for measuring the ocean currents both on the surface and at mid-depth. However the positioning error of fixes in a trajectory varies from 150 m to 1000 m, and thus created difficulty for accurate estimations of the surface and mid-depth currents. Also the reliability of the estimated surface and mid-depth currents requires accurate error estimations.In this study a new sequential method of Argo float surface trajectory tracking and extrapolating is proposed based on Kalman Filter (KF), under the presumption that a surface trajectory of Argo float is dominated by a constant current plus inertial oscillation. This trajectory tracking and extrapolating method is able to reduce the positioning uncertainties of Argo surface trajectories and provides error estimations. When this method was applied to extrapolate the positions when float resurfacing and descending, the estimation error of the mid-depth currents can be reduced. Utilizing this method in the Pacific, surface and mid-depth currents were estimated from surface trajectories of Argo floats from 2001 to 2004, along with their detailed error estimations. The average error for surface currents is about 4.4 cm s^− 1 which is equivalent to the accuracy order (5 cm s^− 1) of the Surface Velocity Program drifters. The estimation error of the mid-depth currents at 1000 db is reduced to about 0.21 cm s^− 1 without considering the effect of vertical shear.This study shows that the surface trajectory from Argo float provides a new means to measure surface circulations in the global ocean at real time, and that the estimated mid-depth current could be one of the important sources to improve the understanding for ocean dynamic.     

Circulation on the Armorican shelf (Bay of Biscay) in autumn

Three drifters drogued at 65 m were launched on a transect on the Armorican shelf of the Bay of Biscay for 4 years. The experiments were conducted in autumn. They revealed a north-westward, poleward current over the 100 m isobath and a very weak eastward current over depths comprised between 120 and 150 m. A model was used to assess the role of residual tidal currents and wind-induced circulation. The results show that the former are quite weak and the latter do not explain the average velocity of over 10 cm s^− 1. It is thought that this current is mainly driven by the density gradient induced by the breakdown of stratification. Hydrological data and satellite images from the period are discussed, in the light of this hypothesis.     

Estimates of vertical mixing during a Lagrangian experiment off the Galician coast

An adjoint 1-D model was used to determine vertical diffusivity coefficients from temperature profiles collected within a filament escaping from the Galician coast following an upwelling event. The optimisation scheme ended with relatively high diffusivity values within the thermocline (9×10⁻⁵ m² s⁻¹). Such high values are relevant for biogeochemical exchanges between surface and deep waters in stratified areas.The optimised values were several orders of magnitude higher than the bulk of diffusivity measurements recorded with a free-falling device; however, the optimisation solution was consistent with the arithmetic mean of the measurements in the thermocline (7.7×10⁻⁵ m² s⁻¹), giving more weight to the few largest values. Below the thermocline, the data assimilation method failed because of the three-dimensional nature of the advective field of the upwelling system. Ignoring this advective forcing in the model led to estimates that were two orders of magnitude too high.The results suggest that turbulent mixing is a random process where a few intense events determine the average mixing that drives the long-term evolution of the water column structure. This statistical property is very important when one wants to use instantaneous diffusivity measurements for modelling purposes.     

Parameterization of air–sea gas fluxes at extreme wind speeds

Air–sea flux measurements of O₂ and N₂ obtained during Hurricane Frances in September 2004 [D'Asaro, E. A. and McNeil, C. L., 2006. Measurements of air–sea gas exchange at extreme wind speeds. Journal Marine Systems, this edition.] using air-deployed neutrally buoyant floats reveal the first evidence of a new regime of air–sea gas transfer occurring at wind speeds in excess of 35 m s^− 1. In this regime, plumes of bubbles 1 mm and smaller in size are transported down from near the surface of the ocean to greater depths by vertical turbulent currents with speeds up to 20−30 cm s^− 1. These bubble plumes mostly dissolve before reaching a depth of approximately 20 m as a result of hydrostatic compression. Injection of air into the ocean by this mechanism results in the invasion of gases in proportion to their tropospheric molar gas ratios, and further supersaturation of less soluble gases. A new formulation for air–sea fluxes of weakly soluble gases as a function of wind speed is proposed to extend existing formulations [Woolf, D.K, 1997. Bubbles and their role in gas exchange. In: Liss, P.S., and Duce, R.A., (Eds.), The Sea Surface and Global Change. Cambridge University Press, Cambridge, UK, pp. 173–205.] to span the entire natural range of wind speeds over the open ocean, which includes hurricanes. The new formulation has separate contributions to air–sea gas flux from: 1) non-supersaturating near-surface equilibration processes, which include direct transfer associated with the air–sea interface and ventilation associated with surface wave breaking; 2) partial dissolution of bubbles smaller than 1 mm that mix into the ocean via turbulence; and 3) complete dissolution of bubbles of up to 1 mm in size via subduction of bubble plumes. The model can be simplified by combining “surface equilibration” terms that allow exchange of gases into and out of the ocean, and “gas injection” terms that only allow gas to enter the ocean. The model was tested against the Hurricane Frances data set. Although all the model parameters cannot be determined uniquely, some features are clear. The fluxes due to the surface equilibration terms, estimated both from data and from model inversions, increase rapidly at high wind speed but are still far below those predicted using the cubic parameterization of Wanninkhof and McGillis [Wannikhof, R. and McGillis, W.R., 1999. A cubic relationship between air–sea CO₂ exchange and wind speed. Geophysical Research Letters, 26:1889–1892.] at high wind speed. The fluxes due to gas injection terms increase with wind speed even more rapidly, causing bubble injection to dominate at the highest wind speeds.     

A dense water outflow from the Ross Sea, Antarctica: Mixing and the contribution of tides

We use hydrographic, current, and microstructure measurements, and tide-forced ocean models, to estimate benthic and interfacial mixing impacting the evolution of a bottom-trapped outflow of dense shelf water from the Drygalski Trough in the northwestern Ross Sea. During summer 2003 an energetic outflow was observed from the outer shelf ( 500 m isobath) to the  1600 m isobath on the continental slope. Outflow thickness was as great as  200 m, and mean speeds were  0.6 m s^− 1 relative to background currents exceeding  1 m s^− 1 that were primarily tidal in origin. No outflow was detected on the slope in winter 2004, although a thin layer of dense shelf water was present on the outer shelf. When the outflow was well-developed, the estimated benthic stress was of order one Pascal and the bulk Froude number over the upper slope exceeded one. Diapycnal scalar diffusivity (K_z) values in the transition region at the top of the outflow, estimated from Thorpe-scale analysis of potential density and measurements of microscale temperature gradient from sensors attached to the CTD rosette, were of order 10^− 3−10^− 2 m² s^− 1. For two cases where the upper outflow boundary was particularly sharply defined, entrainment rate w_e was estimated from K_z and bulk outflow parameters to be  10^− 3 m s^− 1 ( 100 m day^− 1). A tide-forced, three-dimensional primitive equation ocean model with Mellor-Yamada level 2.5 turbulence closure scheme for diapycnal mixing yields results consistent with a significant tidal role in mixing associated with benthic stress and shear within the stratified ocean interior.     

Modelling physical processes of haline stratification in ROFIs with application to the Rhine outflow region

A physical and numerical study is made of the processes governing the stratification and circulation in ROFIs (Regions of Freshwater Influence) where there is an important impact of wind and tides. Observations in the Rhine ROFI showed that the salinity field consists of a mean and a tidally oscillating part. The physical processes are first analysed using the analytical solutions from a one-dimensional two-layer model. A justification is given for the neglection of non-linear advective terms in the equations of momentum and salinity. The dimensionless forms of the solutions can be expressed in terms of a series of dimensionless numbers. It is shown in particular that stratification and cross-shore circulation largely depend on the balance between rotation and turbulent diffusion, which depends in turn on parameters such as the Ekman number, the bottom friction coefficient, the eddy viscosity ratio and the depth of the layer interface. Surface winds either enhance or destroy stratification depending on the wind angle. The response to wind forcing is discussed using classical Ekman theory. To verify the analytical theory numerical tests are performed with a point model including an advanced turbulence closure scheme. Differences arise due to the non-linear interaction between turbulence on the one hand and current shear and stratification on the other hand. It is shown in particular that the amplitude of the tidal forcing and the off-shore horizontal salinity gradient strongly affect the semi-diurnal and semi-monthly variation of stratification. The effect of the wind is found to be in good agreement with the analysis of the two-layer model. Finally, the numerical model is compared with existing observational data in the Rhine ROFI for October 1990.     

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.     

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.     

Interactions of wind and density driven currents in North Sea ROFIs—a model study

Three different versions of a baroclinic three-dimensional circulation model of the North Sea are used to obtain information on the wind and density interactions in the North Sea ROFIs (Regions Of Freshwater Influence): the standard version with fully prognostic treatment of salinity and temperature is compared to a barotropic model run on the same grid on the one hand and to an also fully prognostic model run on a four times coarser grid on the other hand. In order to gain knowledge on the wind and density interactions, two opposing wind directions are chosen for investigation, namely a time of strong north wind, 21st–28th April 1982, and a time of strong southwest wind, 22nd–24th May 1982. In the April case the effect of the salinity gradients on the border of the ROFIs of Rhine, Weser, Ems and Elbe, i.e. along the continental shore, is shown to lead to a clear enhancement of the mean surface currents. In May this result is partly disguised by the additional effect of the thermocline in the deeper parts of the North Sea, i.e. in the classical shelf sea regime region. Nevertheless, the same pattern of enhanced mean surface currents along the coast is detected and is of the same order of magnitude as in the April case. It is thus concluded that although the circulation in the North Sea is reversed by the wind, the density induced component of the general circulation is modified only slightly.     

Dissolved organic matter in shelf waters off the Ría de Vigo (NW Iberian upwelling system)

Net in situ production and export of dissolved organic carbon (DOC) and nitrogen (DON) have been studied in shelf waters off the Ría de Vigo (NW Spain), as part of a comprehensive hydrographic survey carried out from September 1994 to September 1995 with a fortnight periodicity. DOC and DON correlated well (r=+0.78), the slope of the regression line being 12.0±0.7 mol-C mol-N⁻¹, about twice the Redfieldian slope of particulate organic matter, 6.5±0.2 mol-C mol-N⁻¹ (r=+0.95). Labile DOC and DON accumulated in the upper 50 m during the upwelling season (March–September), mainly after prolonged periods of wind relaxation, when horizontal flows were reduced. This labile material represented 50% and 35% of the total (dissolved+particulate) organic carbon and nitrogen susceptible of microbial utilisation, which assert the key contribution of dissolved organic matter (DOM) to the export of new primary production in the NW Iberian upwelling system. This surface excess in shelf waters appeared to be formed into the highly productive Ría de Vigo (a large coastal indentation) at net rates of 4.4 μM-C d⁻¹ and 1.3 μM-C d⁻¹ in the inner and outer segments of the embayment respectively, and subsequently exported to the shelf. Once in the shelf, simple dilution with the inert DOM pool of recently upwelled Eastern North Atlantic Central Water (ENACW) occurred. Eventually, the DOM excess produced during the upwelling season is exported to the adjacent open ocean waters by the coastal circulation. Conversely, during the unproductive downwelling season (October–February), the lowest DOC and DON levels were recorded and export was prevented by the characteristic downwelling front associated to the seasonal poleward slope current.     

Variability of the western Galician upwelling system (NW Spain) during an intensively sampled annual cycle. An EOF analysis approach

The key features of the western Galician shelf hydrography and dynamics are analyzed on a solid statistical and experimental basis. The results allowed us to gather together information dispersed in previous oceanographic works of the region. Empirical orthogonal functions analysis and a canonical correlation analysis were applied to a high-resolution dataset collected from 47 surveys done on a weekly frequency from May 2001 to May 2002. The main results of these analyses are summarized bellow. Salinity, temperature and the meridional component of the residual current are correlated with the relevant local forcings (the meridional coastal wind component and the continental run-off) and with a remote forcing (the meridional temperature gradient at latitude 37°N). About 80% of the salinity and temperature total variability over the shelf, and 37% of the residual meridional current total variability are explained by two EOFs for each variable. Up to 22% of the temperature total variability and 14% of the residual meridional current total variability is devoted to the set up of cross-shore gradients of the thermohaline properties caused by the wind-induced Ekman transport. Up to 11% and 10%, respectively, is related to the variability of the meridional temperature gradient at the Western Iberian Winter Front. About 30% of the temperature total variability can be explained by the development and erosion of the seasonal thermocline and by the seasonal variability of the thermohaline properties of the central waters. This thermocline presented unexpected low salinity values due to the trapping during spring and summer of the high continental inputs from the River Miño recorded in 2001. The low salinity plumes can be traced on the Galician shelf during almost all the annual cycle; they tend to be extended throughout the entire water column under downwelling conditions and concentrate in the surface layer when upwelling favourable winds blow. Our evidences point to the meridional temperature gradient acting as an important controlling factor of the central waters thermohaline properties and in the development and decay of the Iberian Poleward Current.