The Congo River estuary is characterised by a deep canyon that connects the river to the deep ocean by cutting through the continental shelf. Its estuary is influenced by high river discharge and micro-tidal conditions, with a large depth and limited vertical mixing. This restricts the supply of oxygen from the surface waters to the more saline bottom waters, leading to hypoxic and anoxic zones. We

Melt ponds are pools of meltwater forming principally on Arctic sea ice during the melt season. The albedo of melt ponds is a key component of the surface energy balance. For this reason, various melt pond schemes have been developed for climate models. These schemes require assumptions on the physical processes governing melt ponds as well as a knowledge of the atmospheric state, which are not perfectly

We report on the development and implementation of a novel physics-based algorithm that includes the effects of subsurface current on buoyant particles in a Lagrangian particle tracking model, specifically the Larval TRANSport Lagrangian (LTRANS) model. In the upper ocean, the direction and the magnitude of horizontal current change with water depth due to local effects such as the Ekman balance and

A 22-year (1994–2015) HYbrid Coordinate Ocean Model (HYCOM) global reanalysis is used to study the generation mechanism of the counter-wind characteristic of the South China Sea Warm Current (SCSWC) in winter seasons. Our analysis clearly demonstrates that the counter-wind SCSWC in winter is driven by wind relaxation. An Extended Empirical Orthogonal Function (EEOF) analysis of the hourly Sea Surface

We investigate the impact of several parameters on the lifecycle of an anticyclonic eddy lying in a topographic depression (a bowl), similar to the Lofoten Vortex and the Rockall Trough eddy cases. We observe that the vortex merger with submesoscale coherent vortices generated at depth allows the eddy to grow in size, and intensify at depth. Wintertime convection is also shown to directly intensify

This work presents simulations of ocean circulation in the southwestern South Atlantic Ocean with emphasis on the São Paulo Plateau region, a unique physiographic feature at the foot of the continental slope offshore southeast Brazil, where the flow is dominated by ocean eddies and dipoles. Surface dynamics south of 23°S presents a heterogeneous circulation pattern intimately related to the presence

The Northern Pacific Ocean swell environment is generated by extratropical cyclones originating at mid-latitude regions that travel thousands of miles, unimpeded, across the Pacific Ocean. Some of these massive swells reach the Hawaiian Islands causing damage to the shoreline and coastal infrastructure. Understanding the cause and relation between extratropical (EX) cyclone trends and associated swells

Coastal flooding models are used to predict the timing and magnitude of inundation during storms, both for real-time forecasting and long-term design. However, there is a need for faster flooding predictions that also represent flow pathways and barriers at the scales of critical infrastructure. This need can be addressed via subgrid corrections, which use information at smaller scales to ‘correct’

We describe a new framework to represent wave effects on currents, including several approximations (Stokes drift, Bernoulli head, and quasi-static pressure) and a parameterization of the vertical mixing due to wave breaking. The framework improves over existing methods not limited by water depth or monochromatic assumptions. The approximations are validated with spectra and current profiles from a

Non-hydrostatic ocean models are widely required in coastal and estuarine studies but are limited in their application due to their very high computational expense. To alleviate this issue, this study developed a highly computationally efficient non-hydrostatic ocean model named MERF v3.0 (the third version of the Marine Environment Research and Forecasting model) in which a semi-implicit and variable

A parameterization of new ice formation in a wave field has been developed. Detailed processes of new sea ice growth with multiple ice types (nilas, grease, and pancake ice) and associated wave attenuation by different sea ice types are presented and implemented into a coupled ROMS and SWAN model. Sea ice formation includes growth and transformation of frazil, grease, pancake, nilas and solid ice.

Bias correction of reanalysis-based wind stress using scatterometer derived equivalent neutral wind has been a common practice in producing the forcing datasets used in recent global ocean model intercomparisons (OMIPs). Here we systematically evaluate the effect of this wind correction procedure on the simulation of the Pacific North Equatorial Countercurrent (NECC) with multiple sets of model experiments

Forecast errors of subsurface temperature and salinity are substantially reduced with an efficient, two-step, multiscale Ensemble Optimal Interpolation (EnOI) system, applied to a near-global eddy-resolving ocean model. A critical element of any data assimilation system is the background error covariance, which for EnOI is typically a static ensemble of anomalies from a long model run. Here, we construct

An improved estimation of the surface currents in the Levantine Basin of the Mediterranean sea is crucial for a wide range of applications, including pollutants transport and nutrients distribution. This estimation remains challenging due to the scarcity or shortcomings of various data types used for this purpose. In this paper, we present an objective validation of a variational assimilation algorithm

Similar to the Stokes drift in periodic gravity waves, we introduce a second-order Stokes effect on the particle drift in isolated linear surface disturbances (pulses) in shallow water. For a linear disturbance with the shape and length scale of a KdV solitary wave, the model results agree surprisingly well with the observed drift in moderately steep solitary waves in the laboratory, as well as with

Numerical wave models are used for a wide range of applications, from the global ocean to coastal scales. Here we report on significant improvements compared to the previous hindcast detailed in Part 2 of the present study by Rascle and Ardhuin (2013). This result was obtained by updating forcing fields, adjusting the spectral discretization and retuning wind wave growth and swell dissipation parameters

The three-dimensional transport pathways, the time and space scales of vertical transport, and the dispersion characteristics (single-, pair- and multi-particle statistics) of submesoscale currents at an upper-ocean front are investigated using material points (tracer particles) that advect with the local fluid velocity. Coherent vortex filaments and eddies, generated and sustained during the evolution

The formation of Dense Shelf Water (DSW) and Antarctic Bottom Water (AABW) in the Southern Ocean is an essential part of the thermohaline circulation, and understanding this phenomenon is crucial for studying the global climate. AABW is formed as DSW flows down the continental slope and mixes with the surrounding waters. However, DSW formation and its descent remains a poorly resolved issue in many

This study considers the isopycnal eddy transport of mass and passive tracers in eddy-resolving double-gyre quasigeostrophic oceanic circulation. Here we focus on advective transport, whereas a companion paper focuses on eddy-induced diffusive tracer transport. To work towards parameterising eddy tracer transport we quantify the eddy tracer flux using a transport tensor with eddies defined using a

Analysis of the mean (wave-averaged) momentum balance is a common approach used to explain the physical forcing driving wave set-up and mean currents in the nearshore zone. Traditionally this approach has been applied to phase-averaged models but has more recently been applied to phase-resolving models using post-processing, whereby model output is used to calculate each of the momentum terms. While

Various strategies are proposed for enforcing stress-free boundary conditions on deformed domains using a weak-form-based discretization in a Cartesian frame of reference. Due to the irregularity of the computational domain and the particular type of boundary condition, a coupling between the Cartesian velocity components is introduced. As such, a different computational kernel is required for the

Waves in large lakes put coastal communities and vessels under threat, and accurate wave predictions are needed for early warnings. While physics-based numerical wave models such as WAVEWATCH III (WW3) are useful to provide spatial information to supplement in situ observations, they require intensive computational resources. An attractive alternative is machine learning (ML) methods, which can potentially

The implicit equal-weights variational particle smoother (IEWVPS) is a combination of the particle filter (PF) and weak-constraint 4-dimensional variational (4D-Var) method, that inherits the merits of both. The IEWVPS avoids the filter degeneracy of particle filters through an implicit equal-weights scheme and reduces the root mean square deviations (RMSDs) by introducing the 4D-Var method. This method

Wave breaking in the ocean affects the height of extreme waves, energy dissipation, and interaction between the atmosphere and upper ocean. Numerical modelling is a critical step in understanding the physics of wave breaking and offers insight that is hard to gain from field data or experiments. High-fidelity numerical modelling of three-dimensional breaking waves is extremely challenging. Conventional

This study provides an interpretation of isopycnal eddy transport for mass and passive tracers in double-gyre eddy-resolving oceanic circulation. This paper focuses on a transport/diffusion tensor representation of the eddy tracer flux, and a companion paper will focus on advective eddy-induced tracer and mass transports. We use a spatial filter to separate the large and small scales, which leads to

The efficacy of an ocean observing, analysis, and forecasting system for the Mid-Atlantic Bight and the Gulf of Maine is explored using the concept of array modes. The analysis-forecast system is based on a triply nested configuration of the Regional Ocean Modeling System (ROMS) in conjunction with 4-dimensional variational (4D-Var) data assimilation. The array modes identify the degrees of freedom

Storm surge and coastal flooding predictions can require high resolution of critical flow pathways and barriers, typically with simulations using grids/meshes with millions of cells/elements to represent a coastal region. However, the cost of this resolution can slow forecasts during a storm. To add resolution when and where it is needed, previous studies have used adaptive mesh methods, which update

Wave bottom reflection induced by a changing bathymetry has a major impact on the near-shore environment, where harbors and other infrastructure are based. However, this reflection is only partially accounted for in large scale numerical wave prediction models. The present work implements the analytical bottom reflection source term derived in Yevnin and Toledo (2018) to WAVEWATCH III, a numerical

The deep eastern boundary current (DEBC) in the southeast Atlantic Ocean plays an important role in the global meridional overturning circulation, by carrying a comparable fraction of the North Atlantic Deep Water transport toward the Southern Ocean to that of the deep western boundary current. At the same time, the southeast Atlantic Ocean is constantly influenced by energetic Agulhas rings that are

This paper analyses a methodology to identify sub-series of waves parameters and wind speed able to explain the overall variability of the input dataset. To this end, the K-means clustering technique is applied to a 40-year long time series of hindcast data off the Genoa coastline (NW Italy). K-means aims to group the data in a reduced number of clusters, represented by as many modes of variability

Flash rips and surf eddies are transient horizontal structures of the order of 10 to 100 m, which can be generated in the surfzone in the absence of bathymetric irregularities. They are traditionally evaluated in a depth averaged setting which involves intrinsic horizontal shear instabilities and the direct generation of vorticity by short-crested waves. In this article, we revisit the processes of

In recent years, high-resolution global ocean circulation models have begun to include astronomical tidal forcing to concurrently simulate ocean circulation, barotropic and internal tides. The inclusion of tides is of great significance for understanding three-dimensional (3-D) ocean mixing and developing internal tide prediction. However, due to uncertainty in topography and damping terms, improper

Internal tides generated from multiple sites create complicated interference patterns, often making it difficult to identify the origin, direction of propagation, and dissipation of individual waves. To identify these constituent waves, we apply a directional Fourier filter (DFF) method to analyse multimodal, broad-wavenumber, internal wave fields found over varying topography. We apply the method

The D’Urville Sea, East Antarctica, is a major source of Dense Shelf Water (DSW), a precursor of Antarctic Bottom Water (AABW). AABW is a key water mass involved in the worldwide ocean circulation and long-term climate variability. The properties of AABW in global climate models suffer from several biases, making climate projections uncertain. These models are potentially omitting or misrepresenting

Remotely sensed optical data are fundamental to be integrated into biogeochemical models, since the key role of the optical properties on lower trophic dynamics. In this paper, it is shown how ocean optics satellite products are used to constrain the optical algorithm adopted in the operational biogeochemical model system MedBFM that produces analyses, forecasts and reanalyses for the Mediterranean

The influence of ocean surface currents on the global wind-wave field is revisited. State-of-the-art numerical spectral wave model simulations with and without surface currents taken from an eddy resolving global ocean reanalysis were compared. As a global average, simulations forced with currents display significantly better agreement with altimeter derived wave heights. The bias and root mean square

The Red Sea Project (RSP) is based on a coastal lagoon with over 90 pristine islands. The project intends to transform the Red Sea coast into a world-class tourist destination. To better understand the regional dynamics and water exchange scenarios in the lagoon, a high-resolution numerical model is implemented. The general and tidal circulation dynamics are then investigated with a particular focus

It is typical for low-resolution simulations of the ocean to miss not only small- but also large-scale patterns of the flow dynamics compared with their high-resolution analogues. It is usually attributed to the inability of coarse-grid models to properly reproduce the effects of the unresolved small-scale dynamics on the resolved large scales. In part, the reason for that is that coarse-grid models

Two- and three-dimensional models are proposed for ocean-wave attenuation due to scattering by ice floes in the marginal ice zone, in which the attenuation rate depends on the horizontal size of the individual floes. The scattering models are shown to reproduce the behaviour of wave attenuation over short wave periods. However, it is shown that scattering alone cannot explain the observed asymptotic

In order to estimate the spatially varying bottom friction coefficients (BFCs) in the Bohai, Yellow and East China Seas (BYECS), multi-mission satellite observations from TOPEX/Poseidon, Jason-1 and Jason-2 are assimilated into a two-dimensional tidal model with an adjoint method. This work aims to: (1) study the characteristics of spatially varying BFCs in BYECS; and (2) find out the possible sources

This paper examines factors contributing to the remarkable longevity of coherent vortices in the subtropical westward flows. Baroclinic vortices embedded in large-scale vertical shears generate Rossby waves which form an opposite sign eddy associated with inertial Taylor columns on the beta-plane. The combination of the vortex and lee Rossby wave can be viewed as a hetonic dipole that induces meridional

Oscillating turbulent bottom boundary layers (BBLs) occur in lakes and coastal oceans. At the mesoscale level, their kinematics are usually characterized by assuming either laminar or steady turbulent flow, and applying analytical solutions or semi-empirical correlations; e.g. log-law, Stokes’ second problem, inertial dissipation method (IDM), Batchelor fitting to temperature microstructure method

In this article, Part 1 of a two-part series, we run and evaluate the skill of a regional physical–biogeochemical stochastic ocean model based on NEMO. The domain covers the Bay of Biscay at 1/36° resolution, as a case study for open-ocean and coastal shelf dynamics. We generate model ensembles based on assumptions about errors in the atmospheric forcing, the ocean model parameterizations and in the

We describe a numerical treatment for a two-layer coupled model developed for the investigation of submarine landslides and resulting tsunami generation over irregular bathymetry. The landslide model is formulated in a Cartesian coordinate system oriented with the still water level in order to facilitate coupling between water and ground motions. Motions in the upper water layer are simulated using

Monitoring and predicting global ocean biogeochemistry and marine ecosystems is one of the biggest challenges for the coming decade. In operational systems, biogeochemical (BGC) models are forced – or coupled – with physical ocean models fields that are generally constrained by data assimilation of temperature, salinity and sea level anomalies. Yet, while physical data assimilation substantially improves

In this Part 2 article of a two-part series, observations based on satellite missions were used to evaluate the empirical consistency of model ensembles generated via stochastic modelling of ocean physics and biogeochemistry. A high-resolution Bay of Biscay configuration was used as a case study to explore the model error subspace in both the open and coastal ocean. In Part 1 of this work, three experiments

Two methods for setting tides in regional numerical models focussed on replication of internal waves are discussed in this paper. A traditional procedure prescribes the tidal forcing at the boundaries of the model domain. As it follows from our analysis, this method works well in some cases but does not allow to control identical generation conditions in all points of the model grid. Tidal ellipses

Numerical storm surge models are essential to forecasting coastal flood hazard and informing the design of coastal defences. However, such models rely on a variety of inputs, many of which carry uncertainty. An awareness and understanding of the sensitivity of model outputs with respect to those uncertain inputs is therefore essential when interpreting model results. Here, we use an unstructured-mesh

We address ocean modeling capability that has grown exponentially while ocean observation growth has not maintained pace, a situation leading to seemingly degraded forecast skill when model resolution is increased. Skill in predicting ocean instabilities such as mesoscale eddies requires satellite and in situ observations continually correcting numerical model conditions. Observations constrain positions

It has been known for some time that the ocean basins are populated by what is known as “zonal jets”, “deep zonal jets”, or “striations”. Since the oceanic flow is, at least weakly, chaotic, it is not known whether the positions of the jets are “deterministic”, that is, entirely determined by external parameters. A number of theories have been proposed to explain them, some of them predicting zonal

We describe a new 4-D variational assimilation system, called CYCLOCIM, to estimate the climatological seasonal cycle of the residual mean ocean circulation. CYCLOCIM assimilates monthly mean potential temperature and salinity data from the World Ocean Atlas, and CFC-11, CFC-12 and natural radiocarbon measurements for the deep ocean from the Global Data Analysis Project, Version 2. CYCLOCIM’s control

Basin-scale wind curl anomalies in the North Pacific associated with El Niño increase westward Kuroshio transport in the Luzon Strait, whereas those associated with La Niña reduce it. However, the Kuroshio current in the Luzon Strait is further affected by mesoscale eddies and local wind stress. Numerical experiments are conducted to demonstrate that the anomalous low-level circulation over the Philippine

With increased computing power, the horizontal grid-spacing of regional ocean models is decreasing to the point where they can directly simulate lee waves. Although oceanic lee waves can be inherently nonhydrostatic, such as in the abyssal ocean or in the Gulf Stream, regional ocean models are frequently run in hydrostatic mode to avoid the computational expense of solving the nonhydrostatic pressure

Abstract Understanding the trajectory of oil droplets in crossflow jets is important to estimate the pathways of hydrocarbons and to plan countermeasures. We report experimental results of an oil jet with release velocity around 1.5 m/s in a crossflow of 0.3 m/s. The hydrodynamics of the jets obtained with the Large Eddy Simulation (LES) were used to predict the migration of the oil droplets. Two Lagrangian

Abstract Coastal hazard is rarely driven by only one source, as exemplified by the compound flooding from Hurricane Harvey in Galveston Bay in 2017. A 3D creek-to-ocean model is developed to explicitly resolve, without grid nesting, the marine (combination of atmospheric forcing and tides), fluvial and pluvial extremes for this extreme event. We first thoroughly assess the model skills using all available

Abstract Resonant triad interaction is an important mechanism via which the energy of internal tides (ITs) is dissipated. In this study, based on a two-dimensional high-resolution non-hydrostatic model, resonant triad interaction over mid-ocean ridges and corresponding energetics are investigated. Impacts of topographic criticality (the ratio of topographic slope to internal wave slope) on resonant

Abstract Short waves are of key importance for nearshore dynamics, particularly under storms, where they contribute to extreme water levels and drive large morphological changes. Therefore, it is crucial to model accurately the propagation and dissipation of storm waves in the nearshore area. In this paper, field observations collected in contrasted environments and conditions are combined with predictions

Abstract In this work, sea ice concentrations in the regions of Baffin Bay, Hudson Bay, the Labrador Sea and the Gulf of St. Lawrence were estimated using a standalone regional configuration of Los Alamos Sea ice Model, CICE. A sensitivity analysis of the mixed layer depth criteria that affects the predicted sea ice formation/ablation has been carried out to tune the model. Sea ice concentration from