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.

Interaction of tidal flow with a complex topography and bathymetry including headlands, islands, coral reefs and shoals create a rich submesoscale field of tidal jets, vortices, unsteady wakes, lee eddies and free shear layers, all of which impact marine ecology. A unique and detailed view of the submesoscale variability in a part of the Great Barrier Reef lagoon, Australia, that includes a number of small islands was obtained by using a “stereo” pair of 2-m-resolution visible-band images that were acquired just 54 s apart by the WorldView-3 satellite. Near-surface current and vorticity were extracted at a 50-m-resolution from those data using a cross-correlation technique and an optical-flow method, each yielding a similar result. The satellite-derived data are used to test the ability of the second-generation Louvain-la-Neuve ice-ocean model (SLIM), an unstructured-mesh, finite element model for geophysical and environmental flows, to reproduce the details of the currents in the region. The model succeeds in simulating the large-scale (> 1 km) current patterns, such as the main current and the width and magnitude of the jets developing in the gaps between the islands. Moreover, the order of magnitude of the vorticity and the occurrence of some vortices downstream of the islands are correctly reproduced. The smaller scales (< 500 m) are resolved by the model, although various discrepancies with the data are observed. The smallest scales (< 50 m) are unresolved by both the model- and image-derived velocity fields. This study shows that high-resolution models are able to a significant degree to simulate accurately the currents close to a rugged coast. Very-high-resolution satellite oceanography stereo images offer a new way to obtain snapshots of currents near a complex topography that has reefs, islands and shoals, and is a potential resource that could be more widely used to assess the predictive ability of coastal circulation models.