The discontinuous Galerkin (DG) finite element method is well suited for the modelling, with a relatively small number of elements, of three-dimensional flows exhibiting strong velocity or density gradients. Its performance can be highly enhanced by having recourse to r-adaptivity. Here, a vertical adaptive mesh method is developed for DG finite elements. This method, originally designed for finite difference schemes, is based on the vertical diffusion of the mesh nodes, with the diffusivity controlled by the density jumps at the mesh element interfaces.
The mesh vertical movement is determined by means of a conservative arbitrary Lagrangian–Eulerian (ALE) formulation. Though conservativity is naturally achieved, tracer consistency is obtained by a suitable construction of the mesh vertical velocity field, which is defined in such a way that it is fully compatible with the tracer and continuity equations at a discrete level.
The vertically adaptive mesh approach is implemented in the three-dimensional version of the geophysical and environmental flow Second-generation Louvain-la-Neuve Ice-ocean Model (SLIM 3D; www.climate.be/slim). Idealised benchmarks, aimed at simulating the oscillations of a sharp thermocline, are dealt with. Then, the relevance of the vertical adaptivity technique is assessed by simulating thermocline oscillations of Lake Tanganyika. The results are compared to measured vertical profiles of temperature, showing similar stratification and outcropping events. ... Read more

The Columbia River (CR) estuary is characterized by high river discharge and strong tides that generate high velocity flows and sharp density gradients. Its dynamics strongly affects the coastal ocean circulation. Tidal straining in turn modulates the stratification in the estuary. Simulating the hydrodynamics of the CR estuary and plume therefore requires a multi-scale model as both shelf and estuarine circulations are coupled. Such a model has to keep numerical dissipation as low as possible in order to correctly represent the plume propagation and the salinity intrusion in the estuary. Here, we show that the 3D baroclinic discontinuous Galerkin finite element model SLIM 3D is able to reproduce the main features of the CR estuary-to-ocean continuum. We introduce new vertical discretization and mode splitting that allow us to model a region characterized by complex bathymetry and sharp density and velocity gradients. Our model takes into account the major forcings, i.e. tides, surface wind stress and river discharge, on a single multi-scale grid. The simulation period covers the end of spring-early summer of 2006, a period of high river flow and strong changes in the wind regime. SLIM 3D is validated with in-situ data on the shelf and at multiple locations in the estuary and compared with an operational implementation of SELFE. The model skill in the estuary and on the shelf indicate that SLIM 3D is able to reproduce the key processes driving the river plume dynamics, such as the occurrence of bidirectional plumes or reversals of the inner shelf coastal currents. ... Read more

The containment of genetically modified (GM) pollen is an issue of significant concern for many countries. For crops that are bee-pollinated, model predictions of outcrossing rates depend on the movement hypothesis used for the pollinators. Previous work studying pollen spread by honeybees, the most important pollinator worldwide, was based on the assumption that honeybee movement can be well approximated by Brownian motion. A number of recent studies, however, suggest that pollinating insects such as bees perform Lévy flights in their search for food. Such flight patterns yield much larger rates of spread, and so the Brownian motion assumption might significantly underestimate the risk associated with GM pollen outcrossing in conventional crops. In this work, we propose a mechanistic model for pollen dispersal in which the bees perform truncated Lévy flights. This assumption leads to a fractional-order diffusion model for pollen that can be tuned to model motion ranging from pure Brownian to pure Lévy. We parametrize our new model by taking the same pollen dispersal dataset used in Brownian motion modelling studies. By numerically solving the model equations, we show that the isolation distances required to keep outcrossing levels below a certain threshold are substantially increased by comparison with the original predictions, suggesting that isolation distances may need to be much larger than originally thought. ... Read more

In the context of the emergence of extraterrestrial oceanography, we adapted an existing oceanographic model, SLIM (www.climate.be/slim), to the conditions of Titan, a moon of Saturn. The tidal response of the largest southern lake at Titan’s surface, namely Ontario Lacus, is simulated. SLIM solves the 2D, depth-averaged shallow water equations on an unstructured mesh using the discontinuous Galerkin finite element method, which allows for high spatial resolution wherever needed. The impact of the wind forcing, the bathymetry, and the bottom friction is also discussed. The predicted maximum tidal range is about 0.56 m in the southern part of the lake, which is more than twice as large as the previous estimates (see Tokano, Ocean Dyn 60:(4) 803–817
doi:10.1007/s10236-010-0285-3 2010). The patterns and magnitude of the current are also markedly different from those of previous studies: the tidal motion is not aligned with the major axis of the lake and the speed is larger
nearshore. Indeed, the main tidal component rotates clockwise in an exact period of one Titan day and the tidal currents can reach 0.046 ms−1 close to the shores depending on the geometry and the bathymetry. Except for these
specific nearshore regions, the current speed is less than 0.02 ms−1. Circular patterns can be observed offshore, their rotational direction and size varying along the day. ... Read more