Abstract

Turbulence characteristics in the coastal ocean bottom boundary layer are measured using a submersible Particle Image Velocimetry (PIV) system with a sample area of 20 × 20 cm². Measurements are performed in the New York Bight at elevations ranging from 10 cm to about 1.4 m above the seafloor. Recorded data for each elevation consists of 130 s of image pairs recorded at 1 Hz. After processing, the data at each elevation consist of 130 instantaneous spatial velocity distributions within the sample area. The vertical distribution of mean velocity indicates the presence of large-scale shear even at the highest measurement station. The flow also undergoes variations at timescales longer than the present data series.

Spatial spectra of the energy and dissipation are calculated from individual vector maps. The data extend well beyond the peak in the dissipation spectrum and demonstrate that the turbulence is clearly anisotropic even in the dissipation range. The vector maps are also patched together to generate extended velocity distributions using the Taylor hypothesis. Spectra calculated from the extended data cover about three decades in wavenumber space. For the overlapping range the extended spectra show small differences from those determined using the instantaneous distributions. Use of the Taylor hypothesis causes “contamination” of the extended spectra with surface waves. Nevetheless, the results still indicate that the turbulence is also anisotropic at low wavenumbers (energy containing eddies). The vertical component of velocity fluctuations at energy containing scales is significantly damped as the bottom is approached, while the horizontal component maintains a similar energy level at all elevations.

Different methods of estimating the turbulent energy dissipation are compared. Several of these methods are possible only with 2D data, such as that provided by PIV, including a “direct” method, which is based on measured components of the deformation tensor. Estimates based on assumptions of isotropy are typically larger than those based on the direct method (using available velocity gradients and least number of assumptions), but the differences vary from 30% to 100%.

Abstract

Six sets of particle image velocimetry (PIV) data from the bottom boundary layer of the coastal ocean are examined. The data represent periods when the mean currents are higher, of the same order, and much weaker than the wave-induced motions. The Reynolds numbers based on the Taylor microscale (Re _λ ) are 300–440 for the high, 68–83 for the moderate, and 14–37 for the weak mean currents. The moderate–weak turbulence levels are typical of the calm weather conditions at the LEO-15 site because of the low velocities and limited range of length scales. The energy spectra display substantial anisotropy at moderate to high wavenumbers and have large bumps at the transition from the inertial to the dissipation range. These bumps have been observed in previous laboratory and atmospheric studies and have been attributed to a bottleneck effect. Spatial bandpass-filtered vorticity distributions demonstrate that this anisotropy is associated with formation of small-scale, horizontal vortical layers. Methods for estimating the dissipation rates are compared, including direct estimates based on all of the gradients available from 2D data, estimates based on gradients of one velocity component, and those obtained from curve fitting to the energy spectrum. The estimates based on vertical gradients of horizontal velocity are higher and show better agreement with the direct results than do those based on horizontal gradients of vertical velocity. Because of the anisotropy and low turbulence levels, a −5/3 line-fit to the energy spectrum leads to mixed results and is especially inadequate at moderate to weak turbulence levels. The 2D velocity and vorticity distributions reveal that the flow in the boundary layer at moderate speeds consists of periods of “gusts” dominated by large vortical structures separated by periods of more quiescent flows. The frequency of these gusts increases with Re _λ , and they disappear when the currents are weak. Conditional sampling of the data based on vorticity magnitude shows that the anisotropy at small scales persists regardless of vorticity and that most of the variability associated with the gusts occurs at the low-wave-number ends of the spectra. The dissipation rates, being associated with small-scale structures, do not vary substantially with vorticity magnitude. In stark contrast, almost all the contributions to the Reynolds shear stresses, estimated using structure functions, are made by the high- and intermediate-vorticity-magnitude events. During low vorticity periods the shear stresses are essentially zero. Thus, in times with weak mean flow but with wave orbital motion, the Reynolds stresses are very low. Conditional sampling based on phase in the wave orbital cycle does not show any significant trends.

Abstract

Measurements of temperature, velocity, and microscale velocity shear were made from the research submarine F. A. Forel in the near-surface mixed layer of Lake Geneva under conditions of moderate winds of 6–8 m s⁻¹ and of net heating at the water surface. The submarine carried arrays of thermistors and a turbulence package, including airfoil shear probes. The rate of dissipation of turbulent kinetic energy per unit mass, estimated from the variance of the shear, is found to be lognormally distributed and to vary with depth roughly in accordance with the law of the wall at the measurement depths, 15–20 times the significant wave height. Measurements revealed large-scale structures, coherent over the 2.38-m vertical extent sampled by a vertical array of thermistors, consisting of filaments tilted in the wind direction. They are typically about 1.5 m wide, decreasing in width in the upward direction, and are horizontally separated by about 25 m in the downwind direction. Originating in the upper thermocline, they are characterized in the mixed layer by their relatively low temperature and low rates of dissipation of turbulent kinetic energy and by an upward vertical velocity of a few centimeters per second.

Abstract

Results from three separate velocity profilers operated nearly simultaneously in the northwest Atlantic in 1975 are used to form a composite shear spectrum over vertical wavelengths from 100 m down to a few centimeters. This exercise constitutes an intercomparison of the three different measurement techniques and reveals a shear spectrum which is approximately fiat at a WKB-scaled level from k = 0.01 cpm through k ₀ ≈ 0.1 cpm, then falls as k ⁻¹ to a buoyancy wavenumber k ₀ = (N ³/ε)^1/2 determined by the local average Väisälä frequency N and the volume-averaged dissipation rate ε. Various consequences of the observed shear spectral shape are explored.

Abstract

Maps of monthly self-calibrating Palmer Drought Severity Index (SC-PDSI) have been calculated for the period of 1901–2002 for Europe (35°–70°N, 10°W–60°E) with a spatial resolution of 0.5° × 0.5°. The recently introduced SC-PDSI is a convenient means of describing the spatial and temporal variability of moisture availability and is based on the more common Palmer Drought Severity Index. The SC-PDSI improves upon the PDSI by maintaining consistent behavior of the index over diverse climatological regions. This makes spatial comparisons of SC-PDSI values on continental scales more meaningful.

Over the region as a whole, the mid-1940s to early 1950s stand out as a persistent and exceptionally dry period, whereas the mid-1910s and late 1970s to early 1980s were very wet. The driest and wettest summers on record, in terms of the amplitude of the index averaged over Europe, were 1947 and 1915, respectively, while the years 1921 and 1981 saw over 11% and over 7% of Europe suffering from extreme dry or wet conditions, respectively.

Trends in summer moisture availability over Europe for the 1901–2002 period fail to be statistically significant, both in terms of spatial means of the drought index and in the area affected by drought. Moreover, evidence for widespread and unusual drying in European regions over the last few decades is not supported by the current work.

Abstract

Concurrent measurements of the rate of dissipation of turbulent kinetic energy and the void fraction and size distribution of near-surface bubbles are described. Relatively high dissipation rates and void fractions are found in bubble bands produced by Langmuir circulation. The mean dissipation rates observed in the bands are close to those at which the dynamics of algae is significantly affected. The data are used to test basic assumptions underpinning models of subsurface bubbles and associated air–sea gas transfer. A simple model is used to examine the qualitative effect of Langmuir circulation on the vertical diffusion of bubbles and the representation of Langmuir circulation in models of gas transfer. The circulation is particularly effective in vertical bubble transfer when bubbles are injected by breaking waves to depths at which they are carried downward by the circulation against their tendency to rise. The estimated value of the ratio r of the eddy diffusivity of particles (resembling bubbles) K _p to the eddy viscosity K _z depends on depth z and on the form selected for K _z . The effects of nonoverlapping or superimposed Langmuir cells of different size may be very different. Multiple nonoverlapping cells of similar scales with K _z independent of depth can result in concentration profiles that resemble those of a law-of-the-wall K _z . It is demonstrated that model prediction of bubble distributions and of gas transfer (which is related to bubble submergence time) is sensitive to K _z and to the size distribution of Langmuir circulation cells.

Abstract

The climatically sensitive zone of the Arctic Ocean lies squarely within the domain of the North Atlantic oscillation (NAO), one of the most robust recurrent modes of atmospheric behavior. However, the specific response of the Arctic to annual and longer-period changes in the NAO is not well understood. Here that response is investigated using a wide range of datasets, but concentrating on the winter season when the forcing is maximal and on the postwar period, which includes the most comprehensive instrumental record. This period also contains the largest recorded low-frequency change in NAO activity—from its most persistent and extreme low index phase in the 1960s to its most persistent and extreme high index phase in the late 1980s/early 1990s. This long-period shift between contrasting NAO extrema was accompanied, among other changes, by an intensifying storm track through the Nordic Seas, a radical increase in the atmospheric moisture flux convergence and winter precipitation in this sector, an increase in the amount and temperature of the Atlantic water inflow to the Arctic Ocean via both inflow branches (Barents Sea Throughflow and West Spitsbergen Current), a decrease in the late-winter extent of sea ice throughout the European subarctic, and (temporarily at least) an increase in the annual volume flux of ice from the Fram Strait.

Abstract

Results are presented from a set of experiments designed to investigate factors that may influence proxy-based reconstructions of large-scale temperature patterns in past centuries. The factors investigated include 1) the method used to assimilate proxy data into a climate reconstruction, 2) the proxy data network used, 3) the target season, and 4) the spatial domain of the reconstruction. Estimates of hemispheric-mean temperature are formed through spatial averaging of reconstructed temperature patterns that are based on either the local calibration of proxy and instrumental data or a more elaborate multivariate climate field reconstruction approach. The experiments compare results based on the global multiproxy dataset used by Mann and coworkers, with results obtained using the extratropical Northern Hemisphere (NH) maximum latewood tree-ring density set used by Briffa and coworkers. Mean temperature reconstructions are compared for the full NH (Tropics and extratropics, land and ocean) and extratropical continents only, withvarying target seasons (cold-season half year, warm-season half year, and annual mean). The comparisons demonstrate dependence of reconstructions on seasonal, spatial, and methodological considerations, emphasizing the primary importance of the target region and seasonal window of the reconstruction. The comparisons support the generally robust nature of several previously published estimates of NH mean temperature changes in past centuries and suggest that further improvements in reconstructive skill are most likely to arise from an emphasis on the quality, rather than quantity, of available proxy data.

Abstract

Seven sets of 2D particle image velocimetry data obtained in the bottom boundary layer of the coastal ocean along the South Carolina and Georgia coast [at the South Atlantic Bight Synoptic Offshore Observational Network (SABSOON) site] are examined, covering the accelerating and decelerating phases of a single tidal cycle at several heights above the seabed. Additional datasets from a previous deployment are also included in the analysis. The mean velocity profiles are logarithmic, and the vertical distribution of Reynolds stresses normalized by the square of the free stream velocity collapse well for data obtained at the same elevation but at different phases of the tidal cycle. The magnitudes of 〈u′u′〉, 〈w′w′〉, and −〈u′w′〉 decrease with height above bottom in the 25–160-cm elevation range and are consistent with the magnitudes and trends observed in laboratory turbulent boundary layers. If a constant stress layer exists, it is located below 25-cm elevation. Two methods for estimating dissipation rate are compared. The first, a direct estimate, is based on the measured in-plane instantaneous velocity gradients. The second method is based on fitting the resolved part of the dissipation spectrum to the universal dissipation spectrum available in Gargett et al. Being undervalued, the direct estimates are a factor of 2–2.5 smaller than the spectrum-based estimates. Taylor microscale Reynolds numbers for the present analysis range from 24 to 665. Anisotropy is present at all resolved scales. At the transition between inertial and dissipation range the longitudinal spectra exhibit a flatter than −5/3 slope and form spectral bumps. Second-order statistics of the velocity gradients show a tendency toward isotropy with increasing Reynolds number. Dissipation exceeds production at all measurement heights, but the difference varies with elevation. Close to the bottom, the production is 40%–70% of the dissipation, but it decreases to 10%–30% for elevations greater than 80 cm.

Abstract

The rate of dissipation of turbulent kinetic energy has been measured with airfoil probes mounted on an autonomous vehicle, Autosub, on constant-depth legs at 2–10 m below the surface in winds up to 14 m s⁻¹. The observations are mostly in an area limited by fetch to 26 km where the pycnocline depth is about 20 m. At the operational depths of 1.55–15.9 times the significant wave height H _s , and in steady winds of about 11.6 m s⁻¹ when the wave age is 11.7–17.2, dissipation is found to be lognormally distributed with a law-of-the-wall variation with depth and friction velocity. Breaking waves, leaving clouds of bubbles in the water, are detected ahead of the Autosub by a forward-pointing sidescan sonar, and the dissipation is measured when the clouds are subsequently reached. Bands of bubbles resulting from the presence of Langmuir circulation are identified by a semiobjective method that seeks continuity of band structure recognized by both forward- and sideways-pointing sidescan sonars. The times at which bands are crossed are determined and are used to relate dissipation rates and other measured parameters to the location of Langmuir bands. Shear-induced “temperature ramps” are identified with large horizontal temperature gradients. The turbulence measurements are consequently related to breaking waves, the bubble clouds, Langmuir circulation, and temperature ramps, and therefore to the principal processes of mixing in the near-surface layer of the ocean, all of which are found to have associated patterns of turbulent dissipation rates. A large proportion of the highest values of dissipation rate occur within bubble clouds. Dissipation is enhanced in the convergence region of Langmuir circulation at depths to about 10 m, and on the colder, bubble containing, side of temperature ramps associated with water advected downward from near the surface. Near the sea surface, turbulence is dominated by the breaking waves; below a depth of about 6H _s the local vertical mixing in stronger Langmuir circulation cells exceeds that produced on average by the shear-induced eddies that form temperature ramps.