A note on Stokes production of turbulence kinetic energy in the oceanic mixed layer: observations in the Baltic Sea

Shear- and convection-driven turbulence coexists with wind-generated surface gravity waves in the upper ocean. The turbulent Reynolds stresses in the oceanic mixed layer can therefore interact with the shear of the wave-generated Stokes drift velocity to extract energy from the surface waves and inject it into turbulence, thus augmenting the mean shear-driven turbulence. Stokes production of turbulence kinetic energy (TKE) is difficult to measure in the field, since it requires simultaneous measurement of the turbulent stress and the Stokes drift profiles in the water column. However, it is readily inferred using second moment closure models of the oceanic mixed layer provided: (1) wave properties are available, along with the usual water mass properties, and radiative and air–sea fluxes needed to drive the mixed layer model and (2) the model skill can be assessed by comparing the model results against the observed dissipation rates of TKE. Comprehensive measurements made during the Reynolds 2002 campaign in the Baltic Sea have made the estimation of Stokes production possible, and in this paper, we report on the effort and the conclusions reached. Measurements of air–sea exchange parameters and water mass properties during the campaign allowed a mixed layer model to be run and the turbulent stress in the water column to be inferred. Simultaneous wave spectrum measurements enabled Stokes drift profile to be deduced and wave breaking to be included in the model run, and the Stokes production of TKE in the water column estimated. Direct measurements of the TKE dissipation rate from an upward traversing microstructure profiler were used to assure that the model could reproduce the turbulent dissipation rate in the water column. The model results indicate that the Stokes production of TKE in the mixed layer is of the same order of magnitude as the shear production and must therefore be included in mixed layer models.     

Observation-based evaluation of surface wave effects on currents and trajectory forecasts

Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.

     

Numerical study of surface water circulation around Sekisei Lagoon, southwest Japan

This paper discusses the variability of surface currents around Sekisei Lagoon using a nested grid ocean circulation model. We developed a triple-nested grid system that consists of a coarse-resolution (1/60° or ∼1.85 km) model off Taiwan, an intermediate-resolution (1/300° or ∼370 m) model around the Yaeyama Islands, and a fine-resolution (1/900° or ∼123 m) model of Sekisei Lagoon. The nested grid system was forced by wind and heat flux calculated from six-hourly atmospheric reanalysis data and integrated over the period from May to July 2003. The coarse-resolution model was driven by lateral boundary conditions calculated from daily ocean reanalysis data to include realistic variation of the Kuroshio and mesoscale eddies with spatial scales of ∼500–700 km in the open ocean. The tidal forcing was included in the intermediate-resolution model by interpolating sea level data obtained from a data-assimilative tidal model. The results were then used to drive the fine-resolution model to simulate the surface water circulation around Sekisei lagoon. Model results show that (1) currents inside the lagoon are mainly driven by tide and wind; (2) there exists a strong southwestward current along the bottom slope in the southeast portion of the lagoon; the current is mainly driven by remote mesoscale eddies and at times intensified by the local wind; (3) the flow relaxation scheme is effective in reducing biases along the open boundaries. The simulated currents were used to examine the retention and dispersion of passive particles in the surface layer. Results show that the surface dispersion in the strong open ocean current region is significantly higher than that inside the lagoon.     

Observation-based evaluation of surface wave effects on currents and trajectory forecasts

Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.     

A numerical study of island wakes in the Southern California Bight

With the existence of eight substantial islands in the Southern California Bight, the oceanic circulation is significantly affected by island wakes. In this paper a high-resolution numerical model (on a 1 km grid), forced by a high-resolution wind (2 km), is used to study the wakes. Island wakes arise due both to currents moving past islands and to wind wakes that force lee currents in response. A comparison between simulations with and without islands shows the surface enstrophy (i.e., area-integrated square of the vertical component of vorticity at the surface) decreases substantially when the islands in the oceanic model are removed, and the enstrophy decrease mainly takes place in the areas around the islands. Three cases of wake formation and evolution are analyzed for the Channel Islands, San Nicolas Island, and Santa Catalina Island. When flows squeeze through gaps between the Channel Islands, current shears arise, and the bottom drag makes a significant contribution to the vorticity generation. Downstream the vorticity rolls up into submesoscale eddies. When the California Current passes San Nicolas Island from the northwest, a relatively strong flow forms over the shelf break on the northeastern coast and gives rise to a locally large bottom stress that generates anticyclonic vorticity, while on the southwestern side, with an adverse flow pushing the main wake current away from the island, positive vorticity has been generated and a cyclonic eddy detaches into the wake. When the northward Southern California Countercurrent passes the irregular shape of Santa Catalina Island, cyclonic eddies form on the southeastern coast of the island, due primarily to lateral stress rather than bottom stress; they remain coherent as they detach and propagate downstream, and thus they are plausible candidates for the submesoscale “spirals on the sea” seen in many satellite images. Finally, the oceanic response to wind wakes is analyzed in a spin-up experiment with a time-invariant wind that exhibits strips of both positive and negative curl in the island lee. Corresponding vorticity strips in the ocean develop through the mechanism of Ekman pumping.     

Simulation of airflow around rain gauges: Comparison of LES with RANS models

Wind is responsible for systematic errors that affect rain gauge measurements. The authors investigate the use of computational fluid dynamics (CFD) to calculate airflow around rain gauges by applying a high-resolution large eddy simulation (LES) model to determine the flow fields around a measuring system of two rain gauges. The simulated air flow field is characterized by the presence of massive separation which induces the formation and shedding of highly unsteady eddies in the detached shear layers and wakes. Parts of these detached structures occur over the orifice of the rain gauges and may substantially affect the dynamics of the raindrops in this critical region. Non-dissipative LES methods used with fine enough meshes can successfully predict these eddies and their associated fluctuations. The authors compare statistics from LES with steady-state Reynolds averaged Navier–Stokes (RANS) simulations using the k–ε and shear stress transport k–ω turbulence models. They find that both RANS and LES models predict similar mean velocity distributions around the rain gauges. However, they determine the distribution of the resolved turbulent kinetic energy (TKE) to be strongly dependent on the RANS model used. Neither RANS model predictions of TKE are close to those of LES. The authors conclude that the failure of RANS to predict TKE is an important limitation, as TKE is needed to scale the local velocity fluctuations in stochastic models used to calculate the motion of raindrops in the flow field.     

A three-dimensional surface wave–ocean circulation coupled model and its initial testing

A theoretical framework to include the influences of nonbreaking surface waves in ocean general circulation models is established based on Reynolds stresses and fluxes terms derived from surface wave-induced fluctuation. An expression for the wave-induced viscosity and diffusivity as a function of the wave number spectrum is derived for infinite and finite water depths; this derivation allows the coupling of ocean circulation models with a wave number spectrum numerical model. In the case of monochromatic surface wave, the wave-induced viscosity and diffusivity are functions of the Stokes drift. The influence of the wave-induced mixing scheme on global ocean circulation models was tested with the Princeton Ocean Model, indicating significant improvement in upper ocean thermal structure and mixed layer depth compared with mixing obtained by the Mellor–Yamada scheme without the wave influence. For example, the model–observation correlation coefficient of the upper 100-m temperature along 35° N increases from 0.68 without wave influence to 0.93 with wave influence. The wave-induced Reynolds stress can reach up to about 5% of the wind stress in high latitudes, and drive 2–3 Sv transport in the global ocean in the form of mesoscale eddies with diameter of 500–1,000 km. The surface wave-induced mixing is more pronounced in middle and high latitudes during the summer in the Northern Hemisphere and in middle latitudes in the Southern Hemisphere.     

Wave effects in global ocean modeling: parametrizations vs. forcing from a wave model

One of the main challenges of the Copernicus Marine Service is the implementation of coupled ocean/waves systems that accurately estimate the momentum and energy fluxes provided by the atmosphere to the ocean. This study aims to investigate the impact of forcing the Nucleus for European Modelling of the Ocean (NEMO) ocean model with forecasts from the wave model of Météo-France (MFWAM) to improve classical air-sea flux parametrizations, these latter being mostly driven by the 10-m wind. Three wave-related processes, namely, wave-state-dependent stress, Stokes drift-related effects (Stokes-Coriolis force, Stokes drift advection on tracers and on mass), and wave-state-dependent surface turbulence, are examined at a global scale with a horizontal resolution of 0.25°. Three years of sensitivity simulations (2014–2016) show positive feedback on sea surface temperature (SST) and currents when the wave model is used. A significant reduction in SST bias is observed in the tropical Atlantic Ocean. This is mainly due to the more realistic momentum flux provided by the wave model. In mid-latitudes, the most interesting impact occurs during the summer stratification, when the wind is low and the wave model produces a reduction in the turbulence linked with wave breaking. Magnitudes of the large-scale currents in the equatorial region are also improved by 10% compared to observations. In general, it is shown that using the wave model reduces on average the momentum and energy fluxes to the ocean in tropical regions, but increases them in mid-latitudes. These differences are in the order of 10 to 20% compared with the classical parametrizations found in stand-alone ocean models.     

Eddy energy sources and mesoscale eddies in the Sea of Okhotsk

Based on eddy-permitting ocean circulation model outputs, the mesoscale variability is studied in the Sea of Okhotsk. We confirmed that the simulated circulation reproduces the main features of the general circulation in the Sea of Okhotsk. In particular, it reproduced a complex structure of the East-Sakhalin current and the pronounced seasonal variability of this current. We established that the maximum of mean kinetic energy was associated with the East-Sakhalin Current. In order to uncover causes and mechanisms of the mesoscale variability, we studied the budget of eddy kinetic energy (EKE) in the Sea of Okhotsk. Spatial distribution of the EKE showed that intensive mesoscale variability occurs along the western boundary of the Sea of Okhotsk, where the East-Sakhalin Current extends. We revealed a pronounced seasonal variability of EKE with its maximum intensity in winter and its minimum intensity in summer. Analysis of EKE sources and rates of energy conversion revealed a leading role of time-varying (turbulent) wind stress in the generation of mesoscale variability along the western boundary of the Sea of Okhotsk in winter and spring. We established that a contribution of baroclinic instability predominates over that of barotropic instability in the generation of mesoscale variability along the western boundary of the Sea of Okhotsk. To demonstrate the mechanism of baroclinic instability, the simulated circulation was considered along the western boundary of the Sea of Okhotsk from January to April 2005. In April, the mesoscale anticyclonic eddies are observed along the western boundary of the Sea of Okhotsk. The role of the sea ice cover in the intensification of the mesoscale variability in the Sea of Okhotsk was discussed.     

Coupled assimilation for an intermediated coupled ENSO prediction model

The value of coupled assimilation is discussed using an intermediate coupled model in which the wind stress is the only atmospheric state which is slavery to model sea surface temperature (SST). In the coupled assimilation analysis, based on the coupled wind–ocean state covariance calculated from the coupled state ensemble, the ocean state is adjusted by assimilating wind data using the ensemble Kalman filter. As revealed by a series of assimilation experiments using simulated observations, the coupled assimilation of wind observations yields better results than the assimilation of SST observations. Specifically, the coupled assimilation of wind observations can help to improve the accuracy of the surface and subsurface currents because the correlation between the wind and ocean currents is stronger than that between SST and ocean currents in the equatorial Pacific. Thus, the coupled assimilation of wind data can decrease the initial condition errors in the surface/subsurface currents that can significantly contribute to SST forecast errors. The value of the coupled assimilation of wind observations is further demonstrated by comparing the prediction skills of three 12-year (1997–2008) hindcast experiments initialized by the ocean-only assimilation scheme that assimilates SST observations, the coupled assimilation scheme that assimilates wind observations, and a nudging scheme that nudges the observed wind stress data, respectively. The prediction skills of two assimilation schemes are significantly better than those of the nudging scheme. The prediction skills of assimilating wind observations are better than assimilating SST observations. Assimilating wind observations for the 2007/2008 La Niña event triggers better predictions, while assimilating SST observations fails to provide an early warning for that event.     

Numerical reconstruction of trajectory of small-size surface drifter in the Mediterranean sea

In this study, we addressed the effects of wind-induced drift on Lagrangian trajectories of surface sea objects using high-resolution ocean forecast and atmospheric data. Application of stochastic Leeway model for prediction of trajectories drift was considered for the numerical reconstruction of the Elba accident that occurred during the period 21.06.2009–22.06.2009: a person on an inflatable raft was lost in the vicinity of the Elba Island coast; from the initial position, the person on a raft drifted southwards in the open sea and later he was found on a partially deflated raft during rescue operation. For geophysical forcing, we used high-resolution currents from the Mediterranean Forecasting System and atmospheric wind from the European Centre for Medium-Range Weather Forecasts. To investigate the effect of wind on trajectory behavior, numerical simulations were performed using different categories of drifter-like particles similar to a person on an inflatable raft. An algorithm of spatial clustering was used to differentiate the most probable search areas with a high density of particles. Our results showed that the simulation scenarios using particles with characteristics of draft-limited sea drifters provided better prediction of an observed trajectory in terms of the probability density of particles.     

Observed impact of mesoscale circulation on oceanic response to Typhoon Man-Yi (2007)

The oceanic response to a typhoon, where mesoscale ocean circulations co-exist, was investigated by analyzing the independent observations of profiling floats data at three different locations, satellite altimetry data near the eye of Typhoon Man-Yi (2007) before and after its passage, and synthetic aperture radar data taken during the typhoon’s passage. In spite of the nearly symmetric wind pattern around the eye, the distribution of mesoscale eddies had a major impact on the surface currents and mixed layer (ML) depths. As a result, the entrainment of the water below the ML into the ML was affected by the mesoscale circulation and became asymmetric, which accounted for most of the changes observed in the temperature profiles. Changes in the isotherms were driven primarily by the westward propagation of the mesoscale pattern rather than by the typhoon-induced shoaling. The typhoon-induced shoaling could have played a significant role in the generation of high-frequency (e.g., near-inertial) oscillations and/or sub-mesoscale structures. Although a similar or even greater energy flux was observed at the surface, the entrainment within the anticyclonic circulation was weaker than that within the cyclonic circulation and at the edge of the anticyclonic circulation because of the thick pre-existing ML. A strong ocean response to Typhoon Man-Yi (2007) within a cyclonic circulation or at the edge of an anticyclonic circulation, rather than within an anticyclonic eddy, has implications for the role of mesoscale ocean circulations in better understanding and forecasting the typhoon intensity.     

Numerical simulation of the circulation within the Perth Submarine Canyon,Western Australia

Surface and sub-surface currents along the ocean boundary of Western Australia were simulated using Regional Ocean Modelling System (ROMS) to examine the circulation within the Perth Canyon. Two major current systems influenced the circulation within the canyon: (1) The Leeuwin current interacted weakly with the canyon as the majority of the canyon was below the depth of the Leeuwin current and (2) Leeuwin undercurrent interacted strongly with the canyon, forming eddies within the canyon at depths of 400–800 m. The results indicated that within the canyon, the current patterns changed continuously although there were some repeated patterns. Recurrent eddies produced regions where upwelling or downwelling dominated during the model runs. Deep upwelling was stronger within the canyon than elsewhere on the shelf, but vertical transport in the upper ocean was strong everywhere when wind forcing was applied. Upwelling alone appeared to be insufficient to transport nutrients to the euphotic zone because the canyon rims were deep. Increased upwelling, combined with entrapment within eddies and strong upwelling-favourable winds, which could assist mixing, may account for the high productivity attributed to the canyon. The Leeuwin current is otherwise a strong barrier to the upwelling of nutrients.     

An assessment of the impact of surface currents on wave modeling in the Southern Ocean

This paper presents an assessment of the impact of the ocean circulation on modeled wave fields in the Southern Ocean, where a systematic positive bias of the modeled wave height against altimetry data has been reported. The inclusion of ocean currents in the wave model considerably reduces the positive bias of the simulated wave height for high southern latitudes. The decrease of wave energy in the presence of currents is almost exclusively related to the reduction of the relative wind, caused by an overall co-flowing current field associated with the Antarctic Circumpolar Current. Improvements of the model results are also found for the peak period and the mean period against a long-term moored buoy. At the mooring location, the effect of currents is greater for larger and longer waves, suggesting remotely generated swells are more influenced by the currents than local waves. However, an additional qualitative analysis using high-resolution currents in a finer grid nested to the global coarser grid shows that typical resolution of global hydrodynamic reanalysis is not sufficient to resolve mesoscale eddies, and as a consequence, the simulation of mesoscale wave patterns can be compromised. The results are also discussed in terms of the accuracy of forcing fields.     

Numerical simulation and analysis of the mean coastal circulation off California

A two-dimensional numerical model is applied to a coastal ocean wherein alongshore elevation and density gradients, normally calculated by a three-dimensional model, are instead supplied by climatologically averaged data for the California Current System between 25 and 40°N. Surface wind stress is also obtained from climatological data. Both surface and bottom boundary layers are resolved in the model calculations; a second moment turbulence closure submodel supplies vertical diffusivities. Near steady state solutions are possible when surface buoyancy flux is imposed at the surface.Model results are as follows: Southward wind stress produces a broad equatorward current with an embedded coastal jet in accordance with previous studies. Positive wind stress curl reduces the jet current and produces a poleward undercurrent which then surfaces as the curl is increased. The jet currents are reduced and poleward flow increases as bottom steepness increases; to a lesser extent, inclusion of the beta effect has a similar effect. The existence of near bottom, poleward or equatorward flow is explained rather simply in terms of the bottom stress resulting from the alongshore balance of surface wind stress and vertically integrated pressure gradient, the latter involving the alongshore surface elevation and density gradient. A further finding is that the upwelling circulation associated with wind stress is confined to the top 200 to 300 m of the ocean along the California coast.     

Maintenance of the mean kinetic energy in the global ocean by the barotropic and baroclinic energy routes: the roles of JEBAR and Ekman dynamics

In order to determine the maintenance mechanisms of the currents of the global ocean, this study investigates the budget of the annual mean kinetic energy (KE) in a high-resolution (0.1° × 0.1°) semi-global ocean simulation. The analysis is based on a separation of the mean KE using the barotropic (i.e., depth-averaged) and baroclinic (the residual) components of velocity. The barotropic and baroclinic KEs dominate in higher and lower latitudes, respectively, with their global average being comparable to each other. The working rates of wind forcing on the barotropic and baroclinic circulations in the global ocean are 243 and 747 gigawatts, respectively. This study presents at least three new results for the budget of the barotropic KE. Firstly, an energy diagram is rederived to show that the work of the barotropic component of the horizontal pressure gradient (HPG) is connected to the work related to the joint effect of baroclinicity and bottom relief (JEBAR), and then to the budget of potential energy (PE). Secondly, the model analysis shows that the globally averaged work of the barotropic HPG (which is connected to the work related to JEBAR and then to the budget of the PE) is nearly zero. This indicates that the wind- and buoyancy-induced barotropic circulations in the global ocean are of the same strength with opposite sign. Thirdly, it is found that the work of the wind forcing on the barotropic component of the simulated Antarctic Circumpolar Current (ACC) is canceled by the combined effect, in equal measure, of the work of the barotropic HPG and the work of dissipative processes for mean KE. This result makes a significant contribution to the discussion on the depth-integrated momentum balance of the ACC. The barotropic KE is dissipated by the effects of bottom frictional stress, lateral frictional stress, and the Reynolds stress, of which more than half is attributed to an unexpectedly large contribution from biharmonic horizontal friction. Future studies should pay more attention to the role of biharmonic friction used in high-resolution numerical models.     

On the meridional circulation and balance of momentum in the Southern Ocean of POP

 The circulation of the Southern Ocean is studied in the eddy-resolving model POP (Parallel Ocean Program) by an analysis of zonally integrated balances. The TEM formalism (Transformed Eulerian Mean) is extended to include topography and continental boundaries, thus deviations from a zonally integrated state involve transient and standing eddies. The meridional circulation is presented in terms of the Eulerian, eddy-induced, and residual streamfunctions. It is shown that the splitting of the meridional circulation into Ekman and geostrophic transports and the component induced by subgrid and Reynolds stresses is identical to a particular form of the zonally integrated balance of zonal momentum. In this balance, the eddy-induced streamfunctions represent the interfacial form stresses by transient and standing eddies and the residual streamfunction represents the acceleration of the zonal current by density fluxes in a zonally integrated frame. The latter acceleration term is directly related to the surface flux of density and interior fluxes due to the resolved and unresolved eddies. The eddy-induced circulation is extremely vigorous in POP. In the upper ocean a shallow circulation, reversed in comparison to the Deacon cell and mainly due to standing eddies, appears to the north of Drake Passage latitudes, and in the Drake Passage belt of latitudes a deep-reaching cell is induced by transient eddies. In the resulting residual circulation the Deacon cell is largely cancelled and the residual advection of the zonal mean potential density is balanced by diapycnal eddy and subgrid fluxes which are strong in the upper few hundred meters but small in the ocean interior. The balance of zonal momentum is consistent with other eddy-resolving models; a new aspect is the clear identification of density effects in the zonally integrated balance. We show that the wind stress and the stress induced by the residual circulation drive the eastward current, whereas both eddy species result in a braking. Finally, we extend the Johnson–Bryden model of zonal transport to incorporate all relevant terms from the zonal momentum balance. It is shown that wind stress and induction by the residual circulation carry an eastward transport while bottom form stress and the stress induced by standing eddies yield westward components of transport. Received: 26 June 2001 / Accepted: 2 November 2001     

Relationship between thermally forced surface wind and sea surface temperature gradient

An important part of the influence of the oceans on the atmosphere is through direct radiation, sensible heat flux and release of latent heat of evaporation, whereby all of these processes are directly related to the surface temperature of the oceans. A main effect of the atmosphere on the oceans is through momentum exchange at the air-ocean interface, and this process is directly related to the surface wind stress. The sea surface temperature (SST) and the surface wind stress are the two important components in the air-ocean system. If SST is given, a thermally forced boundary layer atmospheric circulation can be simulated. On the other hand, if the surface wind stress is given, the wind-driven ocean waves and ocean currents can be computed.The relationship between SST and surface wind is a coupling of the atmosphere and the oceans. It changes a one-way effect (ocean mechanically driven by atmosphere, or atmosphere thermally forced by oceans) into two-way air-sea interactions. Through this coupling the SST distribution, being an output from an ocean model, leads to the thermally forced surface winds, which feeds back into the ocean model as an additional forcing.Based on Kuo's planetary boundary layer model a linear algebraic equation is established to link the SST gradient with the thermally forced surface wind. The surface wind blows across the isotherms from cold to warm region with some deflection angle to the right (left) in the Northern (Southern) Hemisphere. Results from this study show that the atmospheric stratification reduces both the speed and the deflection angle of the thermally forced wind, however, the Coriolis' effect increases the wind speed in stable atmosphere (Ri>10^–4) and increases the deflection angle.     

Benchmarking the mesoscale variability in global ocean eddy-permitting numerical systems

The role of data assimilation procedures on representing ocean mesoscale variability is assessed by applying eddy statistics to a state-of-the-art global ocean reanalysis (C-GLORS), a free global ocean simulation (performed with the NEMO system) and an observation-based dataset (ARMOR3D) used as an independent benchmark. Numerical results are computed on a 1/4 ^° horizontal grid (ORCA025) and share the same resolution with ARMOR3D dataset. This “eddy-permitting” resolution is sufficient to allow ocean eddies to form. Further to assessing the eddy statistics from three different datasets, a global three-dimensional eddy detection system is implemented in order to bypass the need of regional-dependent definition of thresholds, typical of commonly adopted eddy detection algorithms. It thus provides full three-dimensional eddy statistics segmenting vertical profiles from local rotational velocities. This criterion is crucial for discerning real eddies from transient surface noise that inevitably affects any two-dimensional algorithm. Data assimilation enhances and corrects mesoscale variability on a wide range of features that cannot be well reproduced otherwise. The free simulation fairly reproduces eddies emerging from western boundary currents and deep baroclinic instabilities, while underestimates shallower vortexes that populate the full basin. The ocean reanalysis recovers most of the missing turbulence, shown by satellite products , that is not generated by the model itself and consistently projects surface variability deep into the water column. The comparison with the statistically reconstructed vertical profiles from ARMOR3D show that ocean data assimilation is able to embed variability into the model dynamics, constraining eddies with in situ and altimetry observation and generating them consistently with local environment.     

Floating marine debris surface drift: Convergence and accumulation toward the South Pacific subtropical gyre

Whatever its origin is, a floating particle at the sea surface is advected by ocean currents. Surface currents could be derived from in situ observations or combined with satellite data. For a better resolution in time and space, we use satellite-derived sea-surface height and wind stress fields with a 1/3° grid from 1993 to 2001 to determine the surface circulation of the South Pacific Ocean. Surface currents are then used to compute the Lagrangian trajectories of floating debris. Results show an accumulation of the debris in the eastern-centre region of the South Pacific subtropical gyre ([120°W; 80°W]–[20°S; 40°S]), resulting from a three-step process: in the first two years, mostly forced by Ekman drift, the debris drift towards the tropical convergence zone (30°S). Then they are advected eastward mostly forced by geostrophic currents. They finally reach the eastern-centre region of the South Pacific subtropical gyre from where they could not escape.