Skin layer recovery of free-surface wakes: Relationship to surface renewal and dependence on heat flux and background turbulence

The thermal signatures of free-surface wakes observed in the open ocean show that the recovery of the cool skin layer is related to the degree of surface mixing and to ambient environmental conditions. Wakes produced by two surface-piercing cables of O(10−2 m) in diameter are analyzed using infrared imagery. Under low-wind-speed conditions when the swell and surface current were aligned, the wakes exhibited distinctive patchlike features of O(1 m) in diameter that were generated by the passage of individual waves. The time t* required by the skin layer to recover from these disturbances is compared to the surface-renewal timescale τ used in heat and gas flux models. At low wind speeds, t* is comparable to τ, but at moderate wind speeds the agreement is poor. The spatial and temporal variations in the skin temperature of these wakes are related to a wave Reynolds number used to characterize the strength of the disturbance due to the waves. The recovery process is characterized in terms of the restoring internal energy flux Jr which is proportional to both the initial thickness and the thermal recovery rate of the skin layer and was found to be directly related to the strength of the surface disruption. Comparison of the wake results with laboratory and other field measurements of breaking waves implies that Jr is also a strong function of the net heat flux and background turbulence, which relate directly to the existing environmental conditions such as wind stress and sea state. Our results demonstrate that Jr may vary by several orders of magnitude, depending on the environmental conditions

A Note on the Phillips Spectral Framework for Ocean Whitecaps

There has been a recent upsurge in interest in quantifying kinematic, dynamic, and energetic properties of wave breaking in the open ocean, especially in severe sea states. The underpinning observational and modeling framework is provided by the seminal paper of O. M. Phillips. In this note, a fundamental issue contributing to the scatter in results between investigators is highlighted. This issue relates to the choice of the independent variable used in the expression for the spectral density of the mean breaking crest length per unit area. This note investigates the consequences of the different choices of independent variable presently used by various investigators for validating Phillips model predictions for the spectral density of the breaking crest length per unit area and the associated spectral breaking strength coefficient. These spectral measures have a central role in inferring the associated turbulent kinetic energy dissipation rate and the momentum flux to the upper ocean from breaking wave observations

Microscale wave breaking and air-water gas transfer

Laboratory results showing that the air-water gas transfer velocity k is correlated with mean square wave slope have been cited as evidence that a wave-related mechanism regulates k at low to moderate wind speeds [Jähne et al., 1987; Bock et al., 1999]. Csanady [1990] has modeled the effect of microscale wave breaking on air-water gas transfer with the result that k is proportional to the fractional surface area covered by surface renewal generated during the breaking process. In this report we investigate the role of microscale wave breaking in gas transfer by determining the correlation between k and AB, the fractional area coverage of microscale breaking waves. Simultaneous, colocated infrared (IR) and wave slope imagery is used to verify that AB detected using IR techniques corresponds to the fraction of surface area covered by surface renewal in the wakes of microscale breaking waves. Using measurements of k and AB made at the University of Washington wind-wave tank at wind speeds from 4.6 to 10.7 m s−1, we show that k is linearly correlated with AB, regardless of the presence of surfactants. This result is consistent with Csanady's [1990] model and implies that microscale wave breaking is likely a fundamental physical mechanism contributing to gas transfer

Sea breeze forcing of estuary turbulence and air-water CO2 exchange

The sea breeze is often a dominant meteorological feature at the coastline, but little is known about its estuarine impacts. Measurements at an anchored catamaran and meteorological stations along the Hudson River and New York Bay estuarine system are used to illustrate some basic characteristics and impacts of the feature. The sea breeze propagates inland, arriving in phase with peak solar forcing at seaward stations, but several hours later at up-estuary stations. Passage of the sea breeze front raises the water-to-air CO2 flux by 1–2 orders of magnitude, and drives turbulence comparable to spring tide levels in the upper meter of the water column, where most primary productivity occurs in this highly turbid system. Modeling and observational studies often use remotely-measured winds to compute air-water fluxes (e.g., momentum, CO2), and this leads to a factor of two flux error on sea breeze days during the study

Wave breaking in developing and mature seas

In response to the growing need for robust validation data for Phillips (1985) breaking wave spectral framework, we contribute new field results observed from R/P FLIP for the breaking crest length distributions, Λ, during two different wind-wave conditions, and breaking strength during one wind-wave condition. The first experiment in Santa Barbara Channel had developing seas and the second experiment in the central Pacific Ocean south of Hawaii had mature seas. These are among the first experiments to use dissipation rate measurements probing up into the breaking crest together with simultaneous measurements of breaking crest length distributions. We directly measured the effective breaking strength parameter to be inline image in mature seas with wave age, inline image, of 40–47. We also found that the velocity scale of the breaking dissipation rate peak decreases with increasing wave age. Further, the breaking crest length spectrum falls off slower than the inline image behavior predicted by Phillips (1985). The integrated dissipation rate was consistently higher for mature seas compared to developing seas due to higher energy and momentum fluxes from the wind

A Multisensor Comparison of Ocean Wave Frequency Spectra from a Research Vessel during the Southern Ocean Gas Exchange Experiment

Obtaining accurate measurements of wave statistics from research vessels remains a challenge due to the platform motion. One principal correction is the removal of ship heave and Doppler effects from point measurements. Here, open-ocean wave measurements were collected using a laser altimeter, a Doppler radar microwave sensor, a radar-based system, and inertial measurement units. Multiple instruments were deployed to capture the low- and high-frequency sea surface displacements. Doppler and motion correction algorithms were applied to obtain a full 1D (0.035–1.3 ± 0.2 Hz) wave spectrum. The radar-based system combined with the laser altimeter provided the optimal low- and high-frequency combination, producing a frequency spectrum in the range from 0.035 to 1.2 Hz for cruising speeds ≤3 m s−1 with a spectral rolloff of f−4 Hz and noise floor of −20/−30 dB. While on station, the significant wave height estimates were comparable within 10%–15% among instrumentation. Discrepancies in the total energy and in the spectral shape between instruments arise when the ship is in motion. These differences can be quantified using the spectral behavior of the measurements, accounting for aliasing and Doppler corrections. The inertial sensors provided information on the amplitude of the ship’s modulation transfer function, which was estimated to be ~1.3 ± 0.2 while on station and increased while underway [2.1 at ship-over-ground (SOG) speed; 4.3 m s−1]. The correction scheme presented here is adequate for measurements collected at cruising speeds of 3 m s−1 or less. At speeds greater than 5 m s−1, the motion and Doppler corrections are not sufficient to correct the observed spectral degradation

Variations in Ocean Surface Temperature due to Near-Surface Flow: Straining the Cool Skin Layer

The aqueous thermal boundary layer near to the ocean surface, or skin layer, has thickness O(1 mm) and plays an important role in controlling the exchange of heat between the atmosphere and the ocean. Theoretical arguments and experimental measurements are used to investigate the dynamics of the skin layer under the influence of an upwelling flow, which is imposed in addition to free convection below a cooled water surface. Previous theories of straining flow in the skin layer are considered and a simple extension of a surface straining model is posed to describe the combination of turbulence and an upwelling flow. An additional theory is also proposed, conceptually based on the buoyancy-driven instability of a laminar straining flow cooled from above. In all three theories considered two distinct regimes are observed for different values of the Péclet number, which characterizes the ratio of advection to diffusion within the skin layer. For large Péclet numbers, the upwelling flow dominates and increases the free surface temperature, or skin temperature, to follow the scaling expected for a laminar straining flow. For small Péclet numbers, it is shown that any flow that is steady or varies over long time scales produces only a small change in skin temperature by direct straining of the skin layer. Experimental measurements demonstrate that a strong upwelling flow increases the skin temperature and suggest that the mean change in skin temperature with Péclet number is consistent with the theoretical trends for large Péclet number flow. However, all of the models considered consistently underpredict the measured skin temperature, both with and without an upwelling flow, possibly a result of surfactant effects not included in the models

Riverine skin temperature response to subsurface processes in low wind speeds

Both surface and subsurface processes modulate the surface thermal skin and as such the skin temperature may serve as an indicator for coastal, estuarine, and alluvial processes. Infrared (IR) imagery offers the unique tool to survey such systems, allowing not only to assess temperature variability of the thermal boundary layer, but also to derive surface flow fields through digital particle image velocimetry, optical flow techniques, or spectral methods. In this study, IR time-series imagery taken from a boat moored in the Hudson River estuary is used to determine surface flow, turbulent kinetic energy dissipation rate, and characteristic temperature and velocity length scales. These are linked to subsurface measurements provided by in situ instruments. Under the low wind conditions and weak stratification, surface currents and dissipation rate are found to reflect subsurface mean flow (r^2 = 0.89) and turbulence (r^2 = 0.75). For relatively low dissipation rates, better correlations are obtained by computing dissipation rates directly from wavenumber spectra rather than when having to assume the validity of the Taylor hypothesis. Furthermore, the subsurface dissipation rate scales with the surface length scales (L) and mean flow (U) using ε ∝ U^3/L (r^2 = 0.9). The surface length scale derived from the thermal fields is found to have a strong linear relationship (r^2 = 0.88) to water depth (D) with (D/L) ∼ 13. Such a relation may prove useful for remote bathymetric surveys when no waves are present

The Response of Ocean Skin Temperature to Rain: Observations and Implications for Parameterization of Rain-Induced Fluxes

Rainfall alters the physical and chemical properties of the surface ocean, and its effect on ocean skin temperature and surface heat fluxes is poorly represented in many air-sea interaction models. We present radiometric observations of ocean skin temperature, near-surface (5 cm) temperature from a towed thermistor, and bulk atmospheric and oceanic variables, for 69 rain events observed over the course of 4 months in the Indian Ocean as part of the DYNAMO project. We test a state-of-the-art prognostic model developed by Bellenger et al. (2017, https://doi.org/10.1002/2016JC012429) to predict ocean skin temperature in the presence of rain, and demonstrate a physically motivated modification to the model that improves its performance with increasing rain rate. We characterize the vertical skin-bulk temperature gradient induced by rain and find that it levels off at high rain rates, suggestive of a transition in skin-layer physics that has been previously hypothesized in the literature. We also quantify the small bias that will be present in turbulent sensible heat fluxes parameterized from ocean temperature measurements made at typical “bulk” depths during a rain event. Finally, a wind threshold is observed above which the surface ocean remains well-mixed during a rain event; however, the skin temperature is observed to decrease at all wind speeds in the presence of rain. Keywords: Precipitation, Ocean Skin Temperature, Sea Surface Temperature, Rain Layers, Rain Heat Flux, Flux Parameterizatio

Optical measurements of small deeply penetrating bubble populations generated by breaking waves in the Southern Ocean

Bubble size distributions ranging from 0.5 to 125 μm radius were measured optically during high winds of 13 m s−1 and large-scale wave breaking as part of the Southern Ocean Gas Exchange Experiment. Very small bubbles with radii less than 60 µm were measured at 6–9 m depth using optical measurements of the near-forward volume scattering function and critical scattering angle for bubbles (∼80°). The bubble size distributions generally followed a power law distribution with mean slope values ranging from 3.6 to 4.6. The steeper slopes measured here were consistent with what would be expected near the base of the bubble plume. Bubbles, likely stabilized with organic coatings, were present for time periods on the order of 10–100 s at depths of 6–9 m. Here, relatively young seas, with an inverse wave age of approximately 0.88 and shorter characteristic wave scales, produced lower bubble concentrations, shallower bubble penetration depths, and steep bubble size distribution slopes. Conversely, older seas, with an inverse wave age of 0.70 and longer characteristic wave scales, produced relatively higher bubble concentrations penetrating to 15 m depth, larger bubble sizes, and shallower bubble size distribution slopes. When extrapolated to 4 m depth using a previously published bubble size distribution, our estimates suggest that the deeply penetrating small bubbles measured in the Southern Ocean supplied ∼36% of the total void fraction and likely contributed to the transfer and supersaturation of low-solubility gases