Laboratory experiments were conducted to measure the heat flux from seafoam continuously generated in natural seawater. Using a control volume technique, heat flux was calculated from foam and foam-free surfaces as a function of ambient humidity (ranged from 40% to 78%), air–water temperature difference (ranged from –9°C to 0°C), and wind speed (variable up to 3 m s^-1). Water-surface skin temperature was imaged with a calibrated thermal infrared camera, and near-surface temperature profiles in the air, water, and foam were recorded. Net heat flux from foam surfaces increased with increasing wind speed and was shown to be up to four times greater than a foam-free surface. The fraction of the total heat flux due to the latent heat flux was observed for foam to be 0.75, with this value being relatively constant with wind speed. In contrast, for a foam-free surface the fraction of the total heat flux due to the latent heat flux decreased at higher wind speeds. Temperature profiles through foam are linear and have larger gradients, which increased with wind speed, while foam free surfaces show the expected logarithmic profile and show no variation with temperature. The radiometric surface temperatures show that foam is cooler and more variable than a foam-free surface, and bubble-resolving thermal images show that radiometrically transparent bubble caps and burst bubbles reveal warm foam below the cool surface layer, contributing to the enhanced variability.

Flow in rivers and the coastal ocean is controlled by the frictional force exerted on the water by riverbed or seabed roughness. The frictional force is typically characterized by a drag coefficient C_d, which is estimated from bulk measurements and often assumed constant. Here, we demonstrate a relationship between bed roughness and water surface turbulence that can be used to make remote estimates of C_dC_d, and validate this relationship by comparing remotely sensed estimates of C_d to those from in situ measurements. Thus, our results provide an approach for estimating bottom roughness and C_d based entirely on remotely sensed data, including their spatial variability, which can improve modeling of river discharge and morphodynamics in data-poor regions.

This is the first of a 2‐part series concerning remote observation and wave‐by‐wave analysis of the onset of breaking in the surf zone. In the surf zone, breaking waves drive nearshore circulation, suspend sediment, and promote air–sea gas exchange. Nearshore wave model predictions often diverge from in situ measurements near the break point location because common parameterizations do not account for the rapid changes that occur near the onset of breaking. This work presents extensive methodology to combine data from a line‐scanning LIDAR and thermal infrared cameras to detect breaking, classify breaker type, and measure geometric wave parameters on a wave‐by‐wave basis, which can be used to improve breaker parameterizations. Over 2,600 non‐breaking and 1,600 breaking waves are analyzed from data collected at the USACE Field Research Facility in Duck, NC, including 413 spilling and 111 plunging waves for which the onset of breaking was observed. Wave height is estimated using a spatio‐temporal method for wave tracking that preserves the sea surface elevation maximum and overcomes field of view limitations. Methods for estimating instantaneous wave speed are refined by fitting a skewed Gaussian function to each wave profile before tracking the peaks. Wave slope is estimated from a linear fit to the upper 80% of the wave face, which provides a robust metric and strong correlation with geometric wave slope defined relative to mean sea level. Finally, breaking wave face foam coverage is analyzed to assess common model assumptions about roller length for wave energy dissipation parameterizations.

This is the second of a two‐part series concerning remote observation and wave‐by‐wave analysis of the onset of breaking for spilling and plunging waves in the surf zone. Nearshore phase‐averaged and phase‐resolving wave models parameterize and directly simulate wave breaking and require realistic critical values of key wave parameters, such as the depth‐limited breaking index γ, steepness, or phase speed to initialize wave breaking. Using LIDAR line‐scans and infrared imagery, we observe over 1,600 breaking waves at the US Army Corps of Engineers Field Research Facility (FRF) in Duck, NC, and examine these parameters on a wave‐by‐wave basis at the onset of breaking for 413 spilling and 111 plunging waves. We find that γ is maximum near the onset of breaking at values consistent with those previously observed at the FRF, but that γ for plunging waves (0.73 ≤ γ_P ≤ 0.81) is greater than γ for spilling waves (0.63 ≤ γ_S ≤ 0.71). Direct estimates of wave face slope are maximum at the onset of breaking, approximately 22° for spilling and 30° for plunging waves. Using the relationship between γ and wave face slope, we develop a threshold for the onset of breaking that is a linear function of the two parameters. Wave face slope and γ are further used together to quantify whether a spilling‐ or plunging‐type breaker is more likely. We test the Miche steepness limit on our depth‐limited breaking data and find it correctly predicts only 10% of the plunging breakers and none of the spilling breakers in the surf zone.

Quantifying energy dissipation due to wave breaking remains an essential but elusive goal for studying and modeling air‐sea fluxes of heat, gas, and momentum. Previous observations have shown that lifetimes of bubble plumes and surface foam are directly related to the dissipated energy. Specifically, the foam decay time can be used to estimate the timescale of the subsurface bubble plume and the energy dissipated in the breaking process. A mitigating factor is that the foam decay time can be significantly affected by the surfactant concentration. Here we present an experimental investigation of a new technique that exploits the thermal signature of cooling foam to infer wave breaking dynamics. The experiments were conducted in a laboratory wave tank using artificial seawater with and without the addition of a surfactant. We show that the time from the start of the breaking process to the onset of cooling scales with the bubble plume decay time and the dissipated energy, and is not significantly affected by the presence of additional surfactants. We confirm observations from the field of the spatial variability of the temperature of foam generated by an individual breaking event, which has implications for inferring the spatial variability of bubble plume depth.

Rain-generated lenses of fresher water at the ocean surface affect satellite remote sensing of salinity, mixed-layer dynamics, and air-sea exchange of heat, momentum, and gases. Understanding how rain and wind generate turbulence at the ocean surface is important in modeling the generation and evolution of these fresh lenses. This paper discusses the use of the active controlled flux technique (ACFT) to determine relative levels of turbulence in the top centimeter of the ocean surface in the presence of rain. ACFT measurements were made during the 2016 second Salinity Processes in the Upper-ocean Regional Study (SPURS-2) in the eastern equatorial Pacific Ocean. The data show that at wind speeds below 4 m s^-1, the turbulence dissipation rate at the ocean surface (as parameterized by the water-side surface renewal time constant) is correlated with the instantaneous rain rate. However, at higher wind speeds, the wind stress dominates turbulence production and rain is not a significant source of turbulence. There is also evidence that internal waves can be a significant source of turbulence at the ocean surface under non-raining conditions when a diurnal warm layer is present.

Ship-based X-band radar observations of rain were collected with high spatial resolution during the 2016 and 2017 Salinity Processes in the Upper-ocean Regional Study 2 (SPURS-2) field experiments in the eastern tropical Pacific Ocean. These observations were collected with a repurposed marine radar that is not typically used for weather monitoring. The radar images captured during SPURS-2 show the spatial extent and variable intensity of rain at a horizontal resolution of 180 m within 30 km of the ship. When analyzed alongside collocated measurements of oceanic and atmospheric properties collected during SPURS-2, the radar-derived rain maps enable a clearer understanding of the impact of spatially and temporally varying freshwater fluxes on ocean salinity. Ocean surface freshening, measured by ship gauges, is found to be affected by local rain accumulation, and also by prior rain accumulation in surrounding locations that was measured by radar. In one example, the X-band marine radar measured rain directly ahead of the ship’s path. The ship then sampled a near-surface freshening signature within the time period expected based on the ship speed, ship heading, and rain area measured by the radar. ­

During the second Salinity Processes in the Upper-ocean Regional Study (SPURS-2) field experiments in 2016 and 2017 in the eastern tropical Pacific Ocean, the surface salinity profiler (SSP) measured temperature and salinity profiles in the upper 1.1 m of the ocean. The SSP captured the response of the ocean surface to 35 rain events, providing insight into the generation and evolution of rain-formed fresh layers. This paper describes the measurements made with the SSP during SPURS-2 and quantifies the fresh layers in terms of their vertical salinity gradients between 0.05 m and 1.1 m, ΔS_{1.1 - 0.05 m}. For the 35 rain events sampled with the SSP in 2016 and 2017, the maximum value of ΔS_{1.1 - 0.05 m} is well correlated with the accumulated rainfall. The maximum value of ΔS_{1.1 - 0.05 m} is shown to be linearly proportional to the maximum rain rate and inversely proportional to the wind speed. This wind speed-dependent relationship shows a high degree of scatter, reflecting that the vertical salinity gradient formed during any individual rain event depends on the complex interaction between the local ocean dynamics and the highly variable forcing from rain and wind.

Sea Surface Temperature (SST) modifies the turbulent mixing, drag, and pressure gradients within the marine atmospheric boundary layer that accelerate near‐surface flow from cool to warm SST and decelerate the flow from warm to cool SST. This phenomenon is well documented on scales of 100–1,000 km (the oceanic mesoscale); however, the nature of this air–sea coupling at scales on the order of 1–10 km (the submesoscale) remains unknown. The Advanced Spaceborne Thermal Emission and Reflection Radiometer can be used to study submesoscale phenomena because the high‐resolution infrared and near‐infrared images can used to estimate both SST and wind speed. Observations of dramatic temperature and wind gradients along the Gulf Stream landward edge are used to examine the surface wind response to submesoscale fronts in SST. Our analysis indicates that SST‐induced wind speed perturbations are observed at the scales of order 1–10 km, significantly smaller than previously suggested.

Along-track interferometric synthetic aperture radar measurements of a small fiberglass-hull boat were made in various wind and wave conditions and for different measurement geometries and boat speeds. The data collected show three different cases: 1) the boat signature is visible both in the backscattered power and interferometric phase images; 2) the boat signature is visible only in the interferometric phase image, and 3) the boat signature is not visible in either image. From a preliminary analysis of the data, we conclude that the angle between radar look direction and the nominal velocity vector of the boat significantly affects boat detection. The worst cases for detection are when those vectors are collinear or oppositely directed. The best detection cases appear to be for the case, when boat velocity vector and radar look direction are orthogonal or when the boat is stationary.

Tides and river slope are fundamental characteristics of estuaries, but they are usually undersampled due to deficiencies in the spatial coverage of water level measurements. This study aims to address this issue by investigating the use of airborne lidar measurements to study tidal statistics and river slope in the Columbia River estuary. Eight plane transects over a 12-h period yield at least eight independent measurements of water level at 2.5-km increments over a 65-km stretch of the estuary. These data are fit to a sinusoidal curve and the results are compared to seven in situ gauges. In situ– and lidar-based tide curves agree to within a root-mean-square error of 0.21 m, and the lidar-based river slope estimate of 1.8 × 10^−5 agrees well with the in situ–based estimate of 1.4 × 10^−5 (4 mm km^−1 difference). Lidar-based amplitude and phase estimates are within 10% and 8°, respectively, of their in situ counterparts throughout most of the estuary. Error analysis suggests that increased measurement accuracy and more transects are required to reduce the errors in estimates of tidal amplitude and phase. However, the results validate the use of airborne remote sensing to measure tides and suggest this approach can be used to systematically study water levels at a spatial density not possible with in situ gauges.

Highlights

We measured infrared emissivity of seawater and sea foam in a laboratory experiment.

We developed a method to estimate emissivity for incidence angles up to 85°.

Foam emissivity is higher than water for all wavelengths and angles > 65°.

The difference between foam and water emissivity increases with incidence angle.

This paper shows results of simulations of single point and complex target signatures of along track interferometric SAR radar. The SAR images of the moving target with motion in along track direction show the typical smearing and displacement of the target signature in case of squint-looking radar. The interferometric image shows the rainbow-like signature for moving target. The behaviors of targets signatures in a good agreement with the obtained experimental results.

Thermal infrared (IR) imagery is used to quantify the high spatial and temporal variability of dissipation due to wave breaking in the surf zone. The foam produced in an actively breaking crest, or wave roller, has a distinct signature in IR imagery. A retrieval algorithm is developed to detect breaking waves and extract wave roller length using measurements taken during the Surf Zone Optics 2010 experiment at Duck, NC. The remotely derived roller length and an in situ estimate of wave slope are used to estimate dissipation due to wave breaking by means of the wave-resolving model by Duncan (1981). The wave energy dissipation rate estimates show a pattern of increased breaking during low tide over a sand bar, consistent with in situ turbulent kinetic energy dissipation rate estimates from fixed and drifting instruments over the bar. When integrated over the surf zone width, these dissipation rate estimates account for 40–69% of the incoming wave energy flux. The Duncan (1981) estimates agree with those from a dissipation parameterization by Janssen and Battjes (2007), a wave energy dissipation model commonly applied within nearshore circulation models.

Knowledge of the intensity and spatial-temporal distribution of rainfall over the ocean is critical in understanding the global hydrological cycle. However, rain has proven difficult to measure over the ocean due to problems associated with platform motion and flow distortion combined with the spatial and temporal variability of rainfall itself. Underwater acoustical rain gauges avoid these issues by using the loud and distinctive underwater sound generated by raindrops on the ocean surface to detect and quantify rainfall. Here, the physics and operation of and results from an instrument that uses underwater ambient sound to measure rainfall rate and wind speed are presented. Passive Aquatic Listener (PAL) instruments were mounted on a buoy deployed at Ocean Station P and on 13 Argo profilers that were deployed as part of the US National Aeronautics and Space Administration-sponsored Salinity Processes in the Upper-ocean Regional Study (SPURS) field experiment in the North Atlantic Ocean. The PALs provide near-continuous measurements of rain rate and wind speed during the two-year period over the SPURS study region defined by the Argo profilers. Comparisons of PAL data with rain and wind measured by other techniques, including direct in situ observations and satellite measurements, show good agreement for both rain rate and wind speed.

Under conditions with a large solar flux and low wind speed, a stably stratified warm layer forms at the ocean surface. Evaporation can then lead to an increase in salinity in the warm layer. A large temperature gradient will decrease density enough to counter the density increase caused by the salinity increase, forming a stable positive salinity anomaly at the surface. If these positive salinity anomalies are large in terms of the change in salinity from surface to the base of the gradient, if their areal coverage is a significant fraction of the satellite footprint, and if they persist long enough to be in the satellite field of view, they could be relevant for calibration and validation of L-band microwave salinity measurements. A towed, surface-following profiler was deployed from the N/O Thalassa during the Subtropical Atlantic Surface Salinity Experiment (STRASSE). The profiler measured temperature and conductivity in the surface ocean at depths of 10, 50, and 100 cm. The measurements show that positive salinity anomalies are common at the ocean surface for wind speeds less than 4 m s^-1 when the average daily insolation is >300 W m^-2 and the sea-to-air latent heat flux is greater than zero. A semiempirical model predicts the observed dependence of measured anomalies on environmental conditions. However, the model results and the field data suggest that these ocean surface salinity anomalies are not large enough in terms of the salinity difference to significantly affect microwave radiometric measurements of salinity.

The infrared signature of breaking waves at large incidence angles was investigated using laboratory experiments and a radiometric model. Infrared imagery of the water surface at incidence angles greater than 70° shows multipath reflections for both spilling and plunging waves generated using a programmable wave maker. For the spilling breakers, the multipath signature was initially distinct from the breaking wave front roller signature but then merged to create a single large bright distributed target. For the plunging breakers, the roller and multipath signatures overlapped from the inception of breaking. The radiance of the multipath reflection was higher than the surrounding water for simulated cold sky conditions and lower for a simulated warm sky. A specular double-reflection model successfully predicted the presence of multipath reflection but the magnitude was sensitive to small uncertainties in geometry, wave slope, and input temperatures. The results show that multipath reflection from breaking waves is characteristic of large incidence angle infrared measurements and increases the area and magnitude of the infrared signature of breaking waves compared to the background.

We investigate the relationship between turbulence statistics and coherent structures (CS) in an unstratified reach of the Snohomish River estuary using in situ velocity measurements and surface infrared (IR) imaging. Sequential IR images are used to estimate surface flow characteristics via a particle-image-velocimetry (PIV) technique, and are conditionally sampled to delineate the surface statistics of bottom-generated CS, or boils. In the water column, we find that turbulent kinetic energy (TKE) production exceeds dissipation near the bed but is less than dissipation in the midwater column and that TKE flux divergence closes a significant portion of the measured imbalance. The surface boundary leads to divergence in upwelling CS, and leads to the redistribution of vertical TKE to the horizontal. Very near the surface, statistical anisotropy is observed at length scales larger than the depth H (3%u20135 m), while boil-scale motions of O(1)m are nearly isotropic and exhibit a –5/3 turbulent cascade to smaller scales. Conditional sampling suggests that TKE dissipation in boils is approximately 2 times greater on average than dissipation in ambient flow. Similarly, surface boils are marked by significantly greater velocity variance, upwelling, divergence, and TKE flux divergence than ambient flow regions. Coherent structures and their surface manifestation, therefore, play an important role in the vertical transport of TKE and the water column distribution of dissipation, and are an important component of the TKE budget.

Recently, researchers at the Applied Physics Laboratory at the University of Washington collected a unique dataset by suspending two cameras, one infrared and one electro-optical, from a balloon. This apparatus was then used to image objects drifting on the surface of Lake Washington. The authors took that data and built a processing stream to track the movements of those drifting surface objects.

Air-sea gas fluxes are generally defined in terms of the air/water concentration difference of the gas and the gas transfer velocity, k_L. Because it is difficult to measure k_L in the ocean, it is often parameterized using more easily measured physical properties. Surface divergence theory suggests that infrared (IR) images of the water surface, which contain information concerning the movement of water very near the air-water interface, might be used to estimate k_L. Therefore, a series of experiments testing whether IR imagery could provide a convenient means for estimating the surface divergence applicable to air-sea exchange were conducted in a synthetic jet array tank embedded in a wind tunnel. Gas transfer velocities were measured as a function of wind stress and mechanically generated turbulence; laser-induced fluorescence was used to measure the concentration of carbon dioxide in the top 300 µm of the water surface; IR imagery was used to measure the spatial and temporal distribution of the aqueous skin temperature; and particle image velocimetry was used to measure turbulence at a depth of 1 cm below the air-water interface. It is shown that an estimate of the surface divergence for both wind-shear driven turbulence and mechanically generated turbulence can be derived from the surface skin temperature. The estimates derived from the IR images are compared to velocity field divergences measured by the PIV and to independent estimates of the divergence made using the laser-induced fluorescence data. Divergence is shown to scale with k_L values measured using gaseous tracers as predicted by conceptual models for both wind-driven and mechanically generated turbulence.

Estuarine fronts are well known to influence transport of waterborne constituents such as phytoplankton and sediment, yet due to their ephemeral nature, capturing the physical driving mechanisms and their influence on stratification and mixing is difficult. We investigate a repetitive estuarine frontal feature in the Snohomish River Estuary that results from complex bathymetric shoal/channel interactions. In particular, we highlight a trapping mechanism by which mid-density water trapped over intertidal mudflats converges with dense water in the main channel forming a sharp front. The frontal density interface is maintained via convergent transverse circulation driven by the competition of lateral baroclinic and centrifugal forcing. The frontal presence and propagation give rise to spatial and temporal variations in stratification and vertical mixing. Importantly, this front leads to enhanced stratification and suppressed vertical mixing at the end of the large flood tide, in contrast to what is found in many estuarine systems. The observed mechanism fits within the broader context of frontogenesis mechanisms in which varying bathymetry drives lateral convergence and baroclinic forcing. We expect similar trapping-generated fronts may occur in a wide variety of estuaries with shoal/channel morphology and/or braided channels and will similarly influence stratification, mixing, and transport.

Thermal infrared (IR)-based particle image velocimetry (PIV) is used to measure the evolution of velocity, turbulent kinetic energy (TKE), and the TKE dissipation rate at the water surface in the tidally influenced Snohomish River. Patterns of temperature variability in the IR imagery arise from disruption of the cool-skin layer and are used to estimate the 2-D velocity field. Comparisons of IR-based PIV mean velocity made with a colocated acoustic velocimeter demonstrate high correlation. IR-based PIV provides detailed measurements of previously inaccessible surface velocities and turbulence statistics.

An ongoing limitation of common regression-based infrared (IR) satellite sea surface temperature (SST) algorithms has been the lack of sufficient in situ skin temperature measurements for derivation of the algorithm coefficients. Since IR brightness temperatures respond to the skin temperature, use of the more numerous subsurface observations to tune the algorithms introduces uncertainty into the resulting SST products. Coincident in situ skin and subsurface SST measurements from three years of cruises are used to derive parallel skin and subsurface multichannel SST (MCSST)-type regression algorithms to determine the extent to which improved accuracy can be obtained using the skin measurements. Through use of only coincident measurements, the advantage offered by the greater volume of available subsurface observations is eliminated. Surprisingly, we find no accuracy improvement using skin SST algorithms relative to algorithms derived from the research-grade ship-borne subsurface temperature measurements used in our analysis. However, better accuracy was found relative to algorithms derived from subsurface observations whose accuracy was degraded to that of buoys. The results are robust with regard to satellite resolution, collocation criteria, geographical regions, and time of day.

The accuracy differences are found to be generally consistent with the effects of: (1) increased measurement uncertainty of radiometric measurements relative to research-grade subsurface observations, and (2) differences in spatial variability between the skin SST and temperature-at-depth. The subsurface algorithms are regenerated after degrading the subsurface measurements by adding increasing levels of Gaussian white noise to determine the amplitude of the additional variability required to ensure equal accuracy between the skin and subsurface products. The required supplemental noise ranges between 0.10 and 0.17 K for all data combined and generally decreases with tighter collocation windows and higher-resolution satellite observations. Variogram analysis and filtering of the in situ measurements suggest that differences in measurement uncertainty between the infrared radiometers and the subsurface sensors can explain 0.07–0.10 K of the required noise, while differences in spatial variability with depth can account for up to 0.07–0.10 K of the residual noise. A key consequence is that spatial averages of the skin temperature over satellite footprints of 2 km or more, while potentially biased in the mean, may exhibit less variance relative to point samples of the subsurface temperature than to the actual radiometric skin temperature.

Images of river surface features that reflect the bathymetry and flow in the river have been obtained using remote sensing at microwave, visible, and infrared frequencies. The experiments were conducted at Jetty Island near the mouth of the Snohomish River at Everett, Washington, where complex tidal flow occurs over a varied bathymetry, which was measured as part of these experiments. An X band (9.36 GHz) Doppler radar was operated from the river bank and produced images of normalized radar cross sections and radial surface velocities every 20 min over many tidal cycles. The visible and infrared instruments were flown in an airplane. All of these techniques showed surface evidence of frontal features, flow over a sill, and flow conditioned by a deep hole. These features were modeled numerically, and the model results correspond well to the remote observations. In situ measurements made near the hole showed that changes in measured velocities correlate well with the occurrence of the features in the images. In addition to tidal phase, the occurrence of these features in the imagery depends on tidal range. The surface roughness observed in the imagery appears to be generated by the bathymetry and flow themselves rather than by the modulation of wind waves.

A wind-wave flume is used to determine the extent to which the thermal boundary layer (TBL) at a wind-forced air-water interface is completely renewed from below. We measure skin temperature, T_skin , radiometrically, temperature immediately below the TBL, T_subskin , using a temperature profiler, and net heat flux using the gradient flux technique. The T_skin probability density function, p(T_skin), and surface renewal time scale were measured using passive and active infrared imaging techniques, respectively. We find that the mean percentile rank of T_subskin in p(T_skin) is 99.90, implying that complete surface renewal occurs. This result suggests an alternative to radiometric measurement of T_skin through the simple combination of an infrared camera and an in situ temperature sensor. Comparison of the temperature difference across the TBL to the expected cooling implies that a significant portion of events only partially renew the TBL. This result should impact efforts to improve air-sea transfer models.

A Fourier-based method is presented to process video observations of water waves and calculate the speed distribution of breaking crest lengths. The method has increased efficiency and robust statistics compared with conventional algorithms that assemble distributions from tracking individual crests in the time domain. The method is tested using field observations (video images of whitecaps) of fetch-limited breaking waves during case studies with low (6.7 m s^-1), moderate (8.5 m s^-1), and high (12.6 m s^-1) wind speeds. The method gives distributions consistent with conventional algorithms, including breaking rates that are consistent with direct observations. Results are applied to obtain remote estimates of the energy dissipation associated with wave breaking.

Energy dissipation by breaking water waves is quantified indirectly using remote observations (digital video recordings) and directly using in situ observations (acoustic Doppler velocity profiles). The analysis is the first validation using field data to test the Duncan-Phillips formulation relating energy dissipation to the spectral distribution of whitecap speeds and lengths. Energy dissipation estimates are in agreement over two orders of magnitude, and demonstrate a promising method for routine observation of wave breaking dynamics. Breaking statistics are partitioned into contributions from waves at the peak of the wave-height spectrum and waves at higher frequencies in the spectrum. Peak waves are found to be only 10% of the total breaking rate, however peak waves contribute up to 75% of the total dissipation rate. In addition, breaking statistics are found to depend on the peak wave steepness and the energy input by the wind.

Surface disruptions by boils during strong tidal flows over a rocky sill were observed in thermal infrared imagery collected at the Snohomish River estuary in Washington State. Locations of boil disruptions and boil diameters at the surface were quantified and are used to test an idealized model of vertical boil propagation. The model is developed as a two-dimensional approximation of a three-dimensional vortex loop, and boil vorticity is derived from the flow shear over the sill. Predictions of boil disruption locations were determined from the modeled vertical velocity, the sill depth, and the over-sill velocity. Predictions by the vertical velocity model agree well with measured locations (rms difference 3.0 m) and improve by using measured velocity and shear (rms difference 1.8 m). In comparison, a boil-surfacing model derived from laboratory turbulent mixed-layer wakes agrees with the measurements only when stratification is insignificant.

Microwave modulation by swell waves and its relation to microbreaking waves were investigated in an ocean experiment. Simultaneous collocated microwave and infrared (IR) measurements of wind waves and swell on the ocean were made. The normalized radar cross section sigma₀ and the skin temperature T_skin were both modulated by the swell, but with differing phases. In general, sigma₀ maxima occurred on the front face, whereas T_skin maxima occurred on the rear face of the swell. Infrared imagery has shown that swell-induced microbreaking occurred at or near the swell crest and that the resulting warm wakes occurred on the rear face of the wave. When tilt and range modulations are taken into account, the location of microbreaking also accounts for the maximum of sigma₀ occurring on the front face of the swell. Thus, microbreaking waves generated near the crest of low-amplitude swell can produce microwave and IR signatures with the observed phase. The relationship between microwave and IR signals was further emphasized by comparing microwave Doppler spectra with simultaneous IR and visible images of the sea surface from the same location. When small and microscale breaking waves were present, Doppler spectra exhibited characteristics that are similar to those from whitecaps, having peaks with large Doppler offsets and polarization ratios near unity. When no microbreakers were present, Doppler offsets and polarization ratios were much smaller in accordance with a composite surface scattering theory.

Extensive comparisons are made of the infrared sea surface skin temperature T_skin measured by the Calibrated Infrared In situ Measurement System (CIRIMS) and the Marine-Atmospheric Emitted Radiance Interferometer (M-AERI). Data were collected from four separate deployments on the NOAA research vessel (R/V) Ronald H. Brown and the U.S. Coast Guard (USCG) Polar Sea over a wide range of latitudes and environmental conditions. The deployment time totaled roughly 6 months over a 4-yr period and resulted in over 7000 comparison values.

The mean offset between the two instruments showed that CIRIMS consistently measured a lower temperature than the M-AERI, but by less than 0.10°C. This mean offset was found to be dependent upon sky condition, wind speed, and ship roll, which implies the offset is likely due to uncertainty in the emissivity. The CIRIMS T_skin was recomputed using two alterative emissivity values, one based on emissivity measured by the M-AERI and the other based on a wind-speed-dependent model. In both cases, the recomputation of the CIRIMS T_skin significantly reduced the mean offset. The overall standard deviation between the M-AERI and CIRIMS T_skin was 0.16°C, did not significantly depend on environmental conditions, and was within the expected values of instrument and comparison uncertainties. These comparisons demonstrate the success of CIRIMS in achieving good agreement with the M-AERI over a wide range of conditions. The results also highlight the importance of the sea surface emissivity when measuring the ocean surface skin temperature.

The design and performance of a shipboard-integrated system for underway skin and bulk temperature is presented. The system consists of the Calibrated Infrared In situ Measurement System (CIRIMS) and through-hull temperature sensors. The CIRIMS is an autonomous shipboard radiometer system that measures the sea surface skin temperature T_skin for validation of satellite-derived sea surface temperature products. General design considerations for shipboard radiometer systems are discussed and the philosophy behind the CIRIMS design is presented. Unique features of the design include a constant temperature housing to stabilize instrument drift, a two-point dynamic calibration procedure, separate sky- and sea-viewing radiometers for simultaneous measurements, and the ability to use an infrared transparent window for environmental protection. Laboratory testing and field deployments are used to establish an estimated error budget, which includes instrumentation and environmental uncertainties.

The combination of this testing and field comparison to the Marine-Atmosphere Emitted Radiance Interferometer (M-AERI) and Infrared SST Autonomous Radiometer (ISAR) instruments indicates that the CIRIMS meets the design goal of ±0.10°C accuracy. Temperature and pressure sensors were installed in custom-designed through-hull ports on the NOAA research vessel (R/V) Ronald H. Brown and the University of Washington R/V Thomas G. Thompson to complement the CIRIMS measurements. The ports allow sensors to be installed while the ship is in water and can accommodate a variety of sensors. The combined system provides the ability to measure near-surface temperature profiles from the skin to a depth of 5 m while underway.

Infrared (IR) imagery of microbreaking waves in the ocean and laboratory showed modulation of breaking by swell and paddle-generated waves, respectively. Skin temperature also was modulated by the long waves, with the maxima occurring on the rear face of the long waves in both the laboratory and the field. The IR imagery from the ocean and laboratory showed that long-wave-induced microbreaking occurred at or near the long wave crest, and the resulting warm wakes occurred on the rear face. Thus, microbreaking waves generated near the crest of low-amplitude long waves can produce modulation with the maxima on the rear face. This mechanism was shown to be responsible for modulation of the measured in the laboratory and also likely contributed to the modulation observed in the field.

In aerial surveys conducted during the Tropical Ocean–Global Atmosphere Coupled Ocean-Atmosphere Response Experiment and the low-wind component of the Coupled Boundary Layer Air-Sea Transfer (CBLAST-Low) oceanographic field programs, sea surface temperature (SST) variability at relatively short spatial scales (O(50 m) to O(1 km)) was observed to increase with decreasing wind speed. A unique set of coincident surface and subsurface oceanic temperature measurements from CBLAST-Low is used to investigate the subsurface expression of this spatially organized SST variability, and the SST variability is linked to internal waves. The data are used to test two previously hypothesized mechanisms for SST signatures of oceanic internal waves: a modulation of the cool-skin effect and a modulation of vertical mixing within the diurnal warm layer.

Under conditions of weak winds and strong insolation (which favor formation of a diurnal warm layer), the data reveal a link between the spatially periodic SST fluctuations and subsurface temperature and velocity fluctuations associated with oceanic internal waves, suggesting that some mechanism involving the diurnal warm layer is responsible for the observed signal. Internal-wave signals in skin temperature very closely resemble temperature signals measured at a depth of about 20 cm, indicating that the observed internal-wave SST signal is not a result of modulation of the cool-skin effect. Numerical experiments using a one-dimensional upper ocean model support the notion that internal-wave heaving of the warm-layer base can produce alternating bands of relatively warm and cool SST through the combined effects of surface heating and modulation of wind-driven vertical shear.

This paper describes a unique validation data set acquired from a marine intensive observing period (IOP) conducted on board the NOAA Ship Ronald H. Brown (RHB) during the 2004 Aerosol and Ocean Science Expedition (AEROSE) in the tropical North Atlantic Ocean from 29 February to 26 March 2004. The radiometric and in situ data complement includes marine observations of the Saharan air layer (SAL), including two significant Saharan dust outbreaks over the Atlantic Ocean. Because the impact of tropospheric dust aerosols on satellite infrared (IR) radiometric observations has not yet been fully characterized, the AEROSE data are particularly valuable for IR sensor validation. Shipboard radiometric data germane to satellite validation include observations from a Marine Atmospheric Emitted Radiance Interferometer (M-AERI), a Calibrated Infrared In situ Measurement System (CIRIMS), and Microtops handheld sunphotometers. Among other things, these data provide, for the first time, coincident IR spectra of the dry, dusty SAL from both the uplooking M-AERI and the downlooking Atmospheric Infrared Sounder (AIRS) on board the Aqua satellite. In situ data collected throughout the cruise include Vaisala RS80/90 radiosondes, launched 3-hourly to include Aqua overpass times. The Aqua matchup profiles provide data for validation of AIRS in the presence of high dust loading, along with temperature and water vapor profile retrievals of the SAL. The frequency of sonde launches also enables validation of coincident uplooking M-AERI boundary layer profile retrievals. Preliminary analyses of the AEROSE data are presented here. Focused AEROSE validation studies are the subjects of separate papers.

The oceanic near-surface temperature profile must be accurately characterized to enable precise determination of air–sea heat exchange and satellite retrievals of sea surface temperature. An improved solar transmission parameterization is integrated into existing models for the oceanic warm layer and cool skin within the Tropical Ocean and Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE) bulk flux model to improve the accuracy of predictions of the temperature profile and corresponding heat flux components. Application of the revised bulk flux model to data from 12 diverse cruises demonstrates that the improved parameterization results in significant changes to the predicted cool-skin effect and latent heat fluxes at low wind speeds with high solar radiation due to reduced absorption of solar radiation just below the surface.

Daytime skin-layer cooling is predicted to increase by 0.03 K on average but by more than 0.25 K for winds below 1 m s^-1 and surface irradiance exceeding 900 W m². Predicted changes to the warm-layer correction were smaller but exceeded 0.1 K below 1 m s^-1. Average latent and sensible heat fluxes changed by 1 W m^-2, but the latent flux decreased by 5 W m^-2 near winds of 0.5 m s^-1 and surface irradiance of 950 W m^-2. Comparison with direct observations of skin-layer cooling demonstrated, in particular, that use of the improved solar transmission model resulted in the reduction of previous systematic overestimates of diurnal skin-layer warming. Similar results can be achieved using a simplified treatment of solar absorption with an appropriate estimate of the fraction of incident solar radiation absorbed within the skin layer.

Infrared techniques have been shown to be uniquely capable of detecting and quantifying microscale breaking waves at an air–water interface. Here we extend current capabilities by developing image processing algorithms to measure the crest lengths and velocities of microbreaking waves in a laboratory wind–wave tank. The measurements are used to compute the distribution of crest lengths as a function of speed, ∆(c), introduced by Phillips [1] as a formulation for the distribution of breaking waves. Two methods to determine the crest velocity by applying a particle imaging velocimetry (PIV) algorithm to the infrared imagery are developed and compared to a method based on tracking the centroid of the crest. The crest-PIV method is based on estimation of the velocity of crests identified using a temperature threshold. The image-PIV method is based on a velocity threshold applied to a surface velocity map obtained by using the PIV algorithm over the entire image. Both methods are used to compute the surface turnover rate, which is compared to the frequency of breaking. The methods developed demonstrate the potential for infrared imaging techniques to measure the geometric and kinematic properties of microbreaking waves and are relevant to air–sea flux studies.

Airborne measurements of ocean skin temperature T_s are presented from the Coupled Boundary Layers, Air-Sea Transfer in Low Winds (CBLAST-Low) Pilot Experiment in August 2001 off Martha's Vineyard, MA. We used an infrared (IR) camera with a spatial resolution of 1 m or less and temperature resolution of roughly 0.02°C. Using subframe sampling of the IR imagery, we achieve lower noise and higher spatial resolution than reported by previous investigators using IR radiometers. Fine-scale maps of Ts exhibit horizontal variability over spatial scales ranging from O(10 km) down to O(1 m) that are related to atmospheric and subsurface phenomena under low to moderate wind conditions. Based on supporting measurements of wind and waves, we identify coherent ramp-like structures in T_s with stratification breakdown and meandering streaky features with internal waves. Regional maps of T_s show the standard deviation for the region is ±1.04°C, while the meridional and zonal variability is 0.23°C x km^-1 and 0.27°C x km^-1, respectively. This temperature variability results in meridional and zonal scalar heat flux variability of 7.0 W x m^-2 x km^-1 and 7.6 W x m^-2 x km^-1, respectively. Our results demonstrate the potential for airborne IR imagery accompanied by high-quality ocean data to identify T_s features produced by subsurface circulation.

The role of microscale wave breaking in controlling the air-water transfer of heat and gas is investigated in a laboratory wind-wave tank. The local heat transfer velocity, k_H , is measured using an active infrared technique and the tank-averaged gas transfer velocity, k_G , is measured using conservative mass balances. Simultaneous, colocated infrared and wave slope imagery show that wave-related areas of thermal boundary layer disruption and renewal are the turbulent wakes of microscale breaking waves, or microbreakers. The fractional area coverage of microbreakers, A_B , is found to be 0.10.4 in the wind speed range 4.2–9.3 m s^-1 for cleaned and surfactant-influenced surfaces, and k_H and k_G are correlated with A_B . The correlation of k_H with A_B is independent of fetch and the presence of surfactants, while that for k_G with A_B depends on surfactants. Additionally, A_B is correlated with the mean square wave slope, S², which has shown promise as a correlate for k_G in previous studies. The ratio of k_H measured inside and outside the microbreaker wakes is 3.4, demonstrating that at these wind speeds, up to 75% of the transfer is the direct result of microbreaking. These results provide quantitative evidence that microbreaking is the dominant mechanism contributing to air-water heat and gas transfer at low to moderate wind speeds.

Detailed understanding of the hydrodynamic mechanisms controlling the air-sea exchange of heat and gas requires a method for rapid measurement of the associated transfer velocities. The active controlled flux technique (ACFT), where the temperature decay of a small patch of water heated by an infrared laser is tracked using an infrared imager, has been proposed as a method for making these fast noninvasive measurements of the heat and gas transfer velocities. Here, we report on ACFT measurements of the transfer velocity of heat, k_H, made in the ocean during the Fluxes, Air-sea Interactions and Remote Sensing (FAIRS) experiment (September/October, 2000) and GasEx-01 (January/February, 2001). The results for k_H from both FAIRS and GasEx-01 compare favorably when plotted versus wind speed. However, when scaled to a Schmidt number of 660, the measured k_H values were found to be a factor of two larger than gas transfer velocities measured during GasEx-01. The ACFT-derived k_H values were combined with direct measurements of the bulk-skin oceanic temperature difference to calculate net air-sea heat fluxes during both experiments. Comparison of these values with heat fluxes determined by direct measurements of the latent, sensible, and radiative heat fluxes showed that the ACFT measurements are a factor of seven larger than the direct measurements. One possible theory explaining both the overprediction of the gas transfer velocities and the scale factor between the measured and calculated net heat fluxes is that air-sea exchange is best described by surface penetration rather than surface renewal.

A comparison study of the experimental and theoretical transfer velocities of heat and gas transfer at a wavy air-water interface is undertaken using an active infrared technique and two gas tracers. Applying the surface renewal model formalism [ Danckwerts, 1951], we find that the experimentally evaluated heat transfer velocity is roughly a factor of 2 higher than the transfer velocity of a gas with a low solubility in water when both are referenced to Sc = 600. Potential origins of such a discrepancy are investigated and we propose the use of the random eddy model [ Harriott, 1962 ] to explain our results. The model is an extension of surface renewal to include the eddy approach distance as a new parameter. Numerical simulations of the random eddy model have been performed using a timescale evaluated from the Active Controlled Flux Technique (ACFT) and the characteristics of heat as well as the two gases used in the experiments (He and SF₆). The simulation results show that the transfer velocities of two species, referenced to the same Schmidt number, are different and that their ratio depends on the average value of the approach distance and its distribution. The model as implemented in the present work also predicts changes in the Schmidt number exponent when the hydrodynamics conditions are varied.

The difference in the fugacities of CO₂ across the diffusive sublayer at the ocean surface is the driving force behind the air-sea flux of CO₂. Bulk seawater fugacity is normally measured several meters below the surface, while the fugacity at the water surface, assumed to be in equilibrium with the atmosphere, is measured several meters above the surface. Implied in these measurements is that the fugacity values are the same as those across the diffusive boundary layer. However, temperature gradients exist at the interface due to molecular transfer processes, resulting in a cool surface temperature, known as the skin effect. A warm layer from solar radiation can also result in a heterogeneous temperature profile within the upper few meters of the ocean. Here we describe measurements carried out during a 14-day study in the equatorial Pacific Ocean (GasEx-2001) aimed at estimating the gradients of CO₂ near the surface and resulting flux anomalies. The fugacity measurements were corrected for temperature effects using data from the ship's thermosalinograph, a high-resolution profiler (SkinDeEP), an infrared radiometer (CIRIMS), and several point measurements at different depths on various platforms. Results from SkinDeEP show that the largest cool skin and warm layer biases occur at low winds, with maximum biases of ~4% and +4%, respectively. Time series ship data show an average CO₂ flux cool skin retardation of about 2%. Ship and drifter data show significant CO₂ flux enhancement due to the warm layer, with maximums occurring in the afternoon. Temperature measurements were compared to predictions based on available cool skin parameterizations to predict the skin-bulk temperature difference, along with a warm layer model.

Coherent structures generated beneath laboratory wind waves were investigated using particle image velocimetry. An algorithm was developed to detect these structures and to determine their size, vorticity, and kinetic energy. As the wind speed increased from 4.5 to 11.0 m s^-1, the maximum vorticity of the coherent structures increased by 40%, their average size increased by 20%, their frequency of occurrence increased 400%, and the fraction of the water surface renewed by coherent structures increased from 0.12 to 0.33. Distributions of the total kinetic energy of the coherent structures as a function of size showed that the most energetic eddies occurred in the size range 0.8–1.6 cm in diameter. The near-surface flow could be divided into areas with one of two distinct characteristics: energetic regions occupied by coherent structures and quiescent regions largely devoid of coherent structures. A surface renewal model for air-water exchange was used to calculate the local transfer velocity in both types of regions. The model predicted that the gas transfer velocities in the energetic regions were 2.8 times larger than in the quiescent regions and that 60% of the total air-water gas flux occurred across the energetic regions at all wind speeds. In addition, the rate of turbulent kinetic energy dissipation was ~2.5 times higher in the energetic regions compared to the quiescent regions at all wind speeds.

The second calibration and intercomparison of infrared radiometers (Miami2001) was held at the University of Miami's Rosenstiel School of Marine and Atmospheric Science (RSMAS) during a workshop held from May to June 2001. The radiometers targeted in these two campaigns (laboratory-based and at-sea measurements) are those used to validate the skin sea surface temperatures and land surface temperatures derived from the measurements of imaging radiometers on earth observation satellites. These satellite instruments include those on currently operational satellites and others that will be launched within two years following the workshop. The experimental campaigns were completed in one week and included laboratory measurements using blackbody calibration targets characterized by the National Institute of Standards and Technology (NIST), and an intercomparison of the radiometers on a short cruise on board the R/V F. G. Walton Smith in Gulf Stream waters off the eastern coast of Florida. This paper reports on the results obtained from the shipborne measurements.

Seven radiometers were mounted alongside each other on the R/V Walton Smith for an intercomparison under seagoing conditions. The ship results confirm that all radiometers are suitable for the validation of land surface temperature, and the majority are able to provide high quality data for the more difficult validation of satellite-derived sea surface temperature, contributing less than 0.1 K to the error budget of the validation. The measurements provided by two prototype instruments developed for ship-of-opportunity use confirmed their potential to provide regular reliable data for satellite-derived SST validation. Four high quality radiometers showed agreements within 0.05 K confirming that these instruments are suitable for detailed studies of the dynamics of air–sea interaction at the ocean surface as well as providing high quality validation data. The data analysis confirms the importance of including an accurate correction for reflected sky radiance when using infrared radiometers to measure SST. The results presented here also show the value of regular intercomparisons of ground-based instruments that are to be used for the validation of satellite-derived data products—products that will be an essential component of future assessments of climate change and variability.

We report the results from a laboratory investigation in which microscale breaking waves were detected using an infrared (IR) imager and two-dimensional (2-D) velocity fields were simultaneously measured using particle image velocimetry (PIV). In addition, the local heat transfer velocity was measured using the controlled flux technique. To the best of our knowledge these are the first measurements of the instantaneous 2-D velocity fields generated beneath microscale breaking waves. Careful measurements of the water surface profile enabled us to make accurate estimates of the near-surface velocities using PIV. Previous experiments have shown that behind the leading edge of a microscale breaker the cool skin layer is disrupted creating a thermal signature in the IR image [Jessup et al., J. Geophys. Res. 102, 23145 (1997)]. The simultaneously sampled IR images and PIV data enabled us to show that these disruptions or wakes are typically produced by a series of vortices that form behind the leading edge of the breaker. When the vortices are first formed they are very strong and coherent but as time passes, and they move from the crest region to the back face of the wave, they become weaker and less coherent. The near-surface vorticity was correlated with both the fractional area coverage of microscale breaking waves and the local heat transfer velocity. The strong correlations provide convincing evidence that the wakes produced by microscale breaking waves are regions of high near-surface vorticity that are in turn responsible for enhancing air–water heat transfer rates.

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 A_B , the fractional area coverage of microscale breaking waves. Simultaneous, colocated infrared (IR) and wave slope imagery is used to verify that A_B 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 A_B 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 A_B, 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.