Emerging trends in the sea state of the Beaufort and Chukchi seas

The sea state of the Beaufort and Chukchi seas is controlled by the wind forcing and the amount of ice-free water available to generate surface waves. Clear trends in the annual duration of the open water season and in the extent of the seasonal sea ice minimum suggest that the sea state should be increasing, independent of changes in the wind forcing. Wave model hindcasts from four selected years spanning recent conditions are consistent with this expectation. In particular, larger waves are more common in years with less summer sea ice and/or a longer open water season, and peak wave periods are generally longer. The increase in wave energy may affect both the coastal zones and the remaining summer ice pack, as well as delay the autumn ice-edge advance. However, trends in the amount of wave energy impinging on the ice-edge are inconclusive, and the associated processes, especially in the autumn period of new ice formation, have yet to be well-described by in situ observations. There is an implicit trend and evidence for increasing wave energy along the coast of northern Alaska, and this coastal signal is corroborated by satellite altimeter estimates of wave energy.

Ocean waves across the Arctic: attenuation due to dissipation dominates over scattering for periods longer than 19 s

The poorly understood attenuation of surface waves in sea ice is generally attributed to the combination of scattering and dissipation. Scattering and dissipation have very different effects on the directional and temporal distribution of wave energy, making it possible to better understand their relative importance by analysis of swell directional spreading and arrival times. Here we compare results of a spectral wave model – using adjustable scattering and dissipation attenuation formulations – with wave measurements far inside the ice pack. In this case, scattering plays a negligible role in the attenuation of long swells. Specifically, scattering-dominated attenuation would produce directional wave spectra much broader than the ones recorded, and swell events arriving later and lasting much longer than observed. Details of the dissipation process remain uncertain. Average dissipation rates are consistent with creep effects but are 12 times those expected for a laminar boundary layer under a smooth solid ice plate.

Airborne remote sensing of ocean wave directional wavenumber spectra in the marginal ice zone

Interactions between surface waves and sea ice are thought to be an important, but poorly understood, physical process in the atmosphere-ice-ocean system. In this work, airborne scanning lidar was used to observe ocean waves propagating into the marginal ice zone (MIZ). These represent the first direct spatial measurements of the surface wave field in the polar MIZ. Data were compared against two attenuation models, one based on viscous dissipation and one based on scattering. Both models were capable of reproducing the measured wave energy. The observed wavenumber dependence of attenuation was found to be consistent with viscous processes, while the spectral spreading of higher wavenumbers suggested a scattering mechanism. Both models reproduced a change in peak direction due to preferential directional filtering. Floe sizes were recorded using co-located visible imagery, and their distribution was found to be consistent with ice breakup by the wave field.

Ray-theoretical modeling of secondary microseism P-waves

Secondary microseism sources are pressure fluctuations close to the ocean surface. They generate acoustic P-waves that propagate in water down to the ocean bottom where they are partly reflected, and partly transmitted into the crust to continue their propagation through the Earth. We present the theory for computing the displacement power spectral density of secondary microseism P-waves recorded by receivers in the far field. In the frequency domain, the P-wave displacement can be modeled as the product of (1) the pressure source, (2) the source site effect that accounts for the constructive interference of multiply reflected P-waves in the ocean, (3) the propagation from the ocean bottom to the stations, (4) the receiver site effect. Secondary microseism P-waves have weak amplitudes, but they can be investigated by beamforming analysis. We validate our approach by analyzing the seismic signals generated by Typhoon Ioke (2006) and recorded by the Southern California Seismic Network. Back projecting the beam onto the ocean surface enables to follow the source motion. The observed beam centroid is in the vicinity of the pressure source derived from the ocean wave model WAVEWATCH IIIR. The pressure source is then used for modeling the beam and a good agreement is obtained between measured and modeled beam amplitude variation over time. This modeling approach can be used to invert P-wave noise data and retrieve the source intensity and lateral extent.

DopSCAT: A mission concept for simultaneous measurements of marine winds and surface currents

A radar scatterometer operates by transmitting a pulse of microwave energy toward the ocean's surface and measuring the normalized (per-unit-surface) radar backscatter coefficient (σ°). The primary application of scatterometry is the measurement of near-surface ocean winds. By combining σ° measurements from different azimuth angles, the 10 m vector wind can be determined through a Geophysical Model Function (GMF), which relates wind and backscatter. This paper proposes a mission concept for the measurement of both oceanic winds and surface currents, which makes full use of earlier C-band radar remote sensing experience. For the determination of ocean currents, in particular, the novel idea of using two chirps of opposite slope is introduced. The fundamental processing steps required to retrieve surface currents are given together with their associated accuracies. A detailed description of the mission proposal and comparisons between real and retrieved surface currents are presented. The proposed ocean Doppler scatterometer can be used to generate global surface ocean current maps with accuracies better than 0.2 m/s at a spatial resolution better than 25 km (i.e., 12.5 km spatial sampling) on a daily basis. These maps will allow gaining some insights on the upper ocean mesoscale dynamics. The work lies at a frontier, given that the present inability to measure ocean currents from space in a consistent and synoptic manner represents one of the greatest weaknesses in ocean remote sensing.

Origin and fate of surface drift in the oceanic convergence zones of the eastern Pacific

This study investigates the structure and intensity of the surface pathways connecting to and from the central areas of the large-scale convergence regions of the eastern Pacific Ocean. Surface waters are traced with numerical Lagrangian particles transported in the velocity field of three different ocean models with horizontal resolutions that range from ¼° to 1/32°. The connections resulting from the large-scale convergent Ekman dynamics agree qualitatively but are strongly modulated by eddy variability that introduces meridional asymmetry in the amplitude of transport. Lagrangian forward-in-time integrations are used to analyze the fate of particles originating from the central regions of the convergence zones and highlight specific outflows not yet reported for the southeastern Pacific when using the currents at the highest resolutions (1/12° and 1/32°). The meridional scales of these outflows are comparable to the characteristic width of the fine-scale striation of mean currents.

Characterizing, modelling and understanding the climate variability of the deep water formation in the North-Western Mediterranean Sea

Observing, modelling and understanding the climate-scale variability of the deep water formation (DWF) in the North-Western Mediterranean Sea remains today very challenging. In this study, we first characterize the interannual variability of this phenomenon by a thorough reanalysis of observations in order to establish reference time series. These quantitative indicators include 31 observed years for the yearly maximum mixed layer depth over the period 1980–2013 and a detailed multi-indicator description of the period 2007–2013. Then a 1980–2013 hindcast simulation is performed with a fully-coupled regional climate system model including the high-resolution representation of the regional atmosphere, ocean, land-surface and rivers. The simulation reproduces quantitatively well the mean behaviour and the large interannual variability of the DWF phenomenon. The model shows convection deeper than 1000 m in 2/3 of the modelled winters, a mean DWF rate equal to 0.35 Sv with maximum values of 1.7 (resp. 1.6) Sv in 2013 (resp. 2005). Using the model results, the winter-integrated buoyancy loss over the Gulf of Lions is identified as the primary driving factor of the DWF interannual variability and explains, alone, around 50 % of its variance. It is itself explained by the occurrence of few stormy days during winter. At daily scale, the Atlantic ridge weather regime is identified as favourable to strong buoyancy losses and therefore DWF, whereas the positive phase of the North Atlantic oscillation is unfavourable. The driving role of the vertical stratification in autumn, a measure of the water column inhibition to mixing, has also been analyzed. Combining both driving factors allows to explain more than 70 % of the interannual variance of the phenomenon and in particular the occurrence of the five strongest convective years of the model (1981, 1999, 2005, 2009, 2013). The model simulates qualitatively well the trends in the deep waters (warming, saltening, increase in the dense water volume, increase in the bottom water density) despite an underestimation of the salinity and density trends. These deep trends come from a heat and salt accumulation during the 1980s and the 1990s in the surface and intermediate layers of the Gulf of Lions before being transferred stepwise towards the deep layers when very convective years occur in 1999 and later. The salinity increase in the near Atlantic Ocean surface layers seems to be the external forcing that finally leads to these deep trends. In the future, our results may allow to better understand the behaviour of the DWF phenomenon in Mediterranean Sea simulations in hindcast, forecast, reanalysis or future climate change scenario modes. The robustness of the obtained results must be however confirmed in multi-model studies.

Ocean bottom pressure records from the Cascadia array and short surface gravity waves

The ocean bottom pressure records from eight stations of the Cascadia array are used to investigate the properties of short surface gravity waves with frequencies ranging from 0.2 to 5 Hz. It is found that the pressure spectrum at all sites is a well-defined function of the wind speed U10 and frequency f, with only a minor shift of a few dB from one site to another that can be attributed to variations in bottom properties. This observation can be combined with the theoretical prediction that the ocean bottom pressure spectrum is proportional to the surface gravity wave spectrum E(f) squared, times the overlap integral I(f) which is given by the directional wave spectrum at each frequency. This combination, using E(f) estimated from modeled spectra or parametric spectra, yields an overlap integral I(f) that is a function of the local wave age inline image. This function is maximum for f∕fPM = 8 and decreases by 10 dB for f∕fPM = 2 and f∕fPM = 30. This shape of I(f) can be interpreted as a maximum width of the directional wave spectrum at f∕fPM = 8, possibly equivalent to an isotropic directional spectrum, and a narrower directional distribution toward both the dominant low frequencies and the higher capillary-gravity wave frequencies.

The impact of tides and waves on near-surface suspended sediment concentrations in the English Channel

Numerous ecological problems of continental shelf ecosystems require a refined knowledge of the evolution of suspended sediment concentrations (SSC). The present investigation focuses on the spatial and temporal variabilities of near-surface SSC in coastal waters of the English Channel (western Europe) by exploiting numerical predictions from the Regional Ocean Modeling System ROMS. Extending previous investigations of ROMS performances in the Channel, this analysis refines, with increased spatial and temporal resolutions, the characterization of near-surface SSC patterns revealing areas where concentrations are highly correlated with evolutions of tides and waves. Significant tidal modulations of near-surface concentrations are thus found in the eastern English Channel and the French Dover Strait while a pronounced influence of waves is exhibited in the Channel Islands Gulf. Coastal waters present furthermore strong SSC temporal variations, particularly noticeable during storm events of autumn and winter, with maximum near-surface concentrations exceeding 40 mg l−1 and increase by a factor from 10 to 18 in comparison with time-averaged concentrations. This temporal variability strongly depends on the granulometric distribution of suspended sediments characterized by local bi-modal contributions of silts and sands off coastal irregularities of the Isle of Wight, the Cotentin Peninsula and the southern Dover Strait.

Reconstructability of 3-Dimensional Upper Ocean Circulation from SWOT Sea Surface Height Measurements

Utilizing the framework of effective surface quasi-geostrophic (eSQG) theory, we explored the potential of reconstructing the 3D upper ocean circulation structures, including the balanced vertical velocity (w) field, from high-resolution sea surface height (SSH) data of the planned SWOT satellite mission. Specifically, we utilized the 1/30°, submesoscale-resolving, OFES model output and subjected it through the SWOT simulator that generates the along-swath SSH data with expected measurement errors. Focusing on the Kuroshio Extension region in the North Pacific where regional Rossby numbers range from 0.22 to 0.32, we found that the eSQG dynamics constitutes an effective framework for reconstructing the 3D upper ocean circulation field. Using the modeled SSH data as input, the eSQG-reconstructed relative vorticity (ζ) and w fields are found to reach a correlation of 0.7–0.9 and 0.6–0.7, respectively, in the 1,000m upper ocean when compared to the original model output. Degradation due to the SWOT sampling and measurement errors in the input SSH data for the ζ and w reconstructions is found to be moderate, 5–25% for the 3D ζ field and 15-35% for the 3D w field. There exists a tendency for this degradation ratio to decrease in regions where the regional eddy variability (or Rossby number) increases.

On Growth Rate of Wind Waves: Impact of Short-Scale Breaking Modulations

The wave generation model based on the rapid distortion concept significantly underestimates empirical values of the wave growth rate. As suggested before, inclusion of the aerodynamic roughness modulations effect on the amplitude of the slope-correlated surface pressure could potentially reconcile this model approach with observations. This study explores the role of short-scale breaking modulations to amplify the growth rate of modulating longer waves. As developed, airflow separations from modulated breaking waves result in strong modulations of the turbulent stress in the inner region of the modulating waves. In turn, this leads to amplifying the slope-correlated surface pressure anomalies. As evaluated, such a mechanism can be very efficient for enhancing the wind-wave growth rate by a factor of 2-3.

Inertial dissipation method applied to derive turbulent fluxed over the ocean during the Surface of the Ocean, Fluxes and Interactions with the Atmosphere Atlantic Stratocumulus Transition Experiment (SOFIA/ASTEX) and Structure des Echanges Mer-Atmosphere, Proprietes des Heterogeneites Oceaniques: Recherche Experimentale (SEMAPHORE) experiments with low to moderate wind speed

The transfer coefficients for momentum and heat have been determined for 10 m neutral wind speeds (U-10n) between 0 and 12 m/s using data from the Surface of the Ocean, Fluxes and Interactions with the Atmosphere (SOFIA) and Structure des Echanges Mer-Atmosphere, Proprietes des Heterogeneites Oceaniques: Recherche Experimentale (SEMAPHORE) experiments. The inertial dissipation method was applied to wind and pseudo virtual temperature spectra from a sonic anemometer, mounted on a platform (ship) which was moving through the turbulence held. Under unstable conditions the assumptions concerning the turbulent kinetic energy (TKE) budget appeared incorrect. Using a bulk estimate for the stability parameter, Z/L (where Z is the height and L is the Obukhov length), this resulted in anomalously low drag coefficients compared to neutral conditions. Determining Z/L iteratively, a low rate of convergence was achieved. It was concluded that the divergence of the turbulent transport of TKE was not negligible under unstable conditions. By minimizing the dependence of the calculated neutral drag coefficient on stability, this term was estimated at about -0.65Z/L. The resulting turbulent fluxes were then in close agreement with other studies at moderate wind speed. The drag and exchange coefficients for low wind speeds were found to be C-en x 10(3) = 2.79U(10n)(-1) + 0.66 (U-10n < 5.2 m/s), C-en x 10(3) = C-hn x 10(3) = 1.2 (U-10n greater than or equal to 5.2 m/s), and C-dn x 10(3) = 11.7U(10n)(-2) + 0.668 (U-10n < 5.5 m/s), which imply a rapid increase of the coefficient values as the wind decreased within the smooth flow regime. The frozen turbulence hypothesis and the assumptions of isotropy and an inertial subrange were found to remain valid at these low wind speeds for these shipboard measurements. Incorporation of a free convection parameterization had little effect.