Multi-scale modelling of coastal, shelf and global ocean dynamics
Methods for widening the range of resolved scales (i.e.
performing multi-scale simulations) in ocean sciences and
engineering are developing rapidly, now allowing multiscale
ocean dynamics studies. Having recourse to grid
nesting has been and still is a popular method for increasing
marine models’ resolution when and where needed and for
easily allowing the use of different dynamics at different
resolution. However, this is not the only way to achieve this
goal. Various techniques for modifying locally the grid
resolution or dealing with complex-geometry domains are
available. For instance, composite, structured grids and
unstructured meshes offer an almost infinite geometrical
flexibility.
This special issue focuses on multi-scale modelling
of coastal, shelf and global ocean dynamics, including the
development of new methodologies and schemes and their
applications to ocean process studies. Several articles focus
on numerical aspects of unstructured mesh space discretisation.
Danilov (2010) shows that the noise developing on
triangular meshes on which the location of the variables is
inspired by Arakawa’s C-grid is the largest for regimes
close to geostrophic balance. The noise can be reduced by
specific operators but cannot be entirely suppressed,
“making the triangular C-grid a suboptimal choice for
large-scale ocean modelling”. Then, the companion articles
of Blaise et al. (2010) and Comblen et al. (2010) describe
the space and time discretisation of a three-dimensional,
baroclinic, finite element model based on the discontinuous
Galerkin (DG) technique. This is a significant step forward
in the field of finite element ocean modelling, though this
model cannot yet be regarded as suitable for tackling
realistic applications. Ueckermann and Lermusiaux (2010)
also consider DG finite element techniques, focusing on
biological-physical dynamics in regions with complex
bathymetric features. They compare low- to high-order
discretisations, both in time and space, for regimes in which
biology dominates, advection dominates or terms are
balanced. They find that higher-order schemes on relatively
coarse grids generally perform better than low-order
schemes on fine grids. Kleptsova et al. (2010) assess
various advection schemes for z-coordinate, threedimensional
models in which flooding and drying is taken
into account. In this study, the ability to conserve
momentum is regarded as the main criterion for selecting
a suitable method. On the other hand, Massmann (2010) assesses automatic differentiation for obtaining the adjoint
of an unstructured mesh, tidal model of the European
continental shelf.
Two articles deal with grid nesting. Nash and Hartnett
(2010) introduce a flooding and drying method that can be
used in structured, nested grid systems. This can be
regarded as an alternative to flooding and drying techniques
that are being developed for unstructured mesh models (e.g.
Karna et al. 2010). Then, Haley and Lermusiaux (2010)
derive conservative time-dependent structured finite volume
discretisations and implicit two-way embedded
schemes for primitive equations with the intent to resolve
tidal-to-mesoscale processes over large multi-resolution
telescoping domains with complex geometries including
shallow seas with strong tides, steep shelf breaks and deep
ocean interactions. The authors present realistic simulations
with data assimilation in three regions with diverse
dynamics and show that their developments enhance the
predictive capability, leading to better match with ocean
data.
Various multi-scale, realistic simulations are presented.
Using a finite element ice model and a slab ocean as in
Lietaer et al. (2008), Terwisscha van Scheltinga et al.
(2010) model the Canadian Arctic Archipelago, focusing on
the pathways for freshwater and sea-ice transport from the
Arctic Ocean to the Labrador Sea and the Atlantic Ocean.
The unstructured mesh can represent the complex geometry
and narrow straits at high resolution and allows improving
transports of water masses and sea ice. Walters et al. (2010)
have recourse to an unstructured mesh model to study tides
and current in Greater Cook Strait (New Zealand). They
identify the mechanisms causing residual currents. By
means of the unstructured mesh Finite Volume Coastal
Ocean Model (FVCOM), Wang et al. (2010) study the
hydrodynamics of the Bohai Sea. Xu et al. (2010) simulate
coastal and urban inundation due to storm surges along US
East and Gulf Coasts. A sensitivity analysis reveals the
importance of precise topographic data and the need for a
bottom drag coefficient accounting for the presence of
mangroves. Finally, Yang and Khangaonkar (2010) resort to
FVCOM to simulate the three-dimensional circulation of
Puget Sound, a large complex estuary system in the Pacific
Northwest coastal ocean, including variable forcing from
tides, the atmosphere and river inflows. Comparisons of
model estimates with measurements for tidal elevation,
velocity, temperature and salinity are deemed to be
promising, from larger-scale circulation features to nearshore
tide flats.
This special issue suggests that numerical techniques for
multi-scale space discretisation are progressively becoming
mature. One direction for future progress lies in the
improvement of time discretisation methods for the new
generation models, so that they can successfully compete
with finite difference, structured mesh models based on
(almost) constant resolution grids that have been developed
and used over the past 40 years (e.g. Griffies et al. 2009).