Multi-scale modeling of Puget Sound using an unstructured-grid coastal ocean model: from tide flats to estuaries and coastal waters

Water circulation in Puget Sound, a large complex estuary system in the Pacific Northwest coastal ocean of the United States, is governed by multiple spatially and temporally varying forcings from tides, atmosphere (wind, heating/cooling, precipitation/evaporation, pressure), and river inflows. In addition, the hydrodynamic response is affected strongly by geomorphic features, such as fjord-like bathymetry and complex shoreline features, resulting in many distinguishing characteristics in its main and sub-basins. To better understand the details of circulation features in Puget Sound and to assist with proposed nearshore restoration actions for improving water quality and the ecological health of Puget Sound, a high-resolution (around 50 m in estuaries and tide flats) hydrodynamic model for the entire Puget Sound was needed. Here, a three-dimensional circulation model of Puget Sound using an unstructured-grid finite volume coastal ocean model is presented. The model was constructed with sufficient resolution in the nearshore region to address the complex coastline, multi-tidal channels, and tide flats. Model open boundaries were extended to the entrance of the Strait of Juan de Fuca and the northern end of the Strait of Georgia to account for the influences of ocean water intrusion from the Strait of Juan de Fuca and the Fraser River plume from the Strait of Georgia, respectively. Comparisons of model results, observed data, and associated error statistics for tidal elevation, velocity, temperature, and salinity indicate that the model is capable of simulating the general circulation patterns on the scale of a large estuarine system as well as detailed hydrodynamics in the nearshore tide flats. Tidal characteristics, temperature/salinity stratification, mean circulation, and river plumes in estuaries with tide flats are discussed.     

Linking 1D coastal ocean modelling to environmental management: an ensemble approach

The use of a one-dimensional interdisciplinary numerical model of the coastal ocean as a tool contributing to the formulation of ecosystem-based management (EBM) is explored. The focus is on the definition of an experimental design based on ensemble simulations, integrating variability linked to scenarios (characterised by changes in the system forcing) and to the concurrent variation of selected, and poorly constrained, model parameters. The modelling system used was previously specifically designed for the use in “data-rich” areas, so that horizontal dynamics can be resolved by a diagnostic approach and external inputs can be parameterised by nudging schemes properly calibrated. Ensembles determined by changes in the simulated environmental (physical and biogeochemical) dynamics, under joint forcing and parameterisation variations, highlight the uncertainties associated to the application of specific scenarios that are relevant to EBM, providing an assessment of the reliability of the predicted changes. The work has been carried out by implementing the coupled modelling system BFM-POM1D in an area of Gulf of Trieste (northern Adriatic Sea), considered homogeneous from the point of view of hydrological properties, and forcing it by changing climatic (warming) and anthropogenic (reduction of the land-based nutrient input) pressure. Model parameters affected by considerable uncertainties (due to the lack of relevant observations) were varied jointly with the scenarios of change. The resulting large set of ensemble simulations provided a general estimation of the model uncertainties related to the joint variation of pressures and model parameters. The information of the model result variability aimed at conveying efficiently and comprehensibly the information on the uncertainties/reliability of the model results to non-technical EBM planners and stakeholders, in order to have the model-based information effectively contributing to EBM.     

Tidal-driven dynamics and mixing processes in a coastal ocean model with wetting and drying

A three-dimensional sigma coordinate numerical model with wetting and drying (WAD) and a Mellor–Yamada turbulence closure scheme has been used in an idealized island configuration to evaluate how tidally driven dynamics and mixing are affected by inundation processes. Comprehensive sensitivity experiments evaluate the influence of various factors, including tidal amplitudes (from 1- to 9-m range), model grid size (from 2 to 16 km), stratification, wind, rotation, and the impact of WAD on the mixing. The dynamics of the system involves tidally driven basin-scale waves (propagating anticlockwise in the northern hemisphere) and coastally trapped waves propagating around the island in an opposite direction. The evolutions of the surface mixed layer (SML) and the bottom boundary layer (BBL) under different forcing have been studied. With small amplitude tides, wind-driven mixing dominates and the thickness of the SML increases with time, while with large-amplitude tides, tidal mixing dominates and the thickness of the BBL increases with time. The inclusion of WAD in the simulations increases bottom stress and impacts the velocities, the coastal waves, and the mixing. However, the impact of WAD is complex and non-linear. For example, WAD reduces near-coast currents during flood but increases currents during ebb as water drains from the island back to the sea. The impacts of WAD, forcing, and model parameters on the dynamics are summarized by an analysis of the vorticity balance for the different sensitivity experiments.     

A coupled oscillator model of shelf and ocean tides

The resonances of tides in the coupled open ocean and shelf are modeled by a mechanical analogue consisting of a damped driven larger mass and spring (the open-ocean) connected to a damped smaller mass and spring (the shelf). When both masses are near resonance, the addition of even a very small mass can significantly affect the oscillations of the larger mass. The influence of the shelf is largest if the shelf is resonant with weak friction. In particular, an increase of friction on a near-resonant shelf can, perhaps surprisingly, lead to an increase in ocean tides. On the other hand, a shelf with large friction has little effect on ocean tides. Comparison of the model predictions with results from numerical models of tides during the ice ages, when lower sea levels led to a much reduced areal extent of shelves, suggests that the predicted larger tidal dissipation then is related to the ocean basins being close to resonance. New numerical simulations with a forward global tide model are used to test expectations from the mechanical analogue. Setting friction to unrealistically large values in Hudson Strait yields larger North Atlantic _M2$M_2$ amplitudes, very similar to those seen in a simulation with the Hudson Strait blocked off. Thus, as anticipated, a shelf with very large friction is nearly equivalent in its effect on the open ocean to the removal of the shelf altogether. Setting friction in shallow waters throughout the globe to unrealistically large values yields even larger open ocean tidal amplitudes, similar to those found in simulations of ice-age tides. It thus appears that larger modeled tides during the ice ages can be a consequence of enhanced friction in shallower water on the shelf in glacial times as well as a reduced shelf area then. Single oscillator and coupled oscillator models for global tides show that the maximum extractable power for human use is a fraction of the present dissipation rate, which is itself a fraction of global human power consumption.     

Regional and global effects of southern ocean constraints in a global model

Global ocean circulation models do not usually take high-latitude processes into account in an adequate form due to a limited model domain or insufficient resolution. Without the processes in key areas contributing to the lower part of the global thermohaline circulation, the characteristics and flow of deep and bottom waters often remain unrealistic in these models. In this study, various sections of the Bremerhaven Regional Ice Ocean Simulation model results are combined with a global inverse model by using temperature, salinity, and velocity constraints for the Hamburg Large Scale Geostrophic ocean general circulation model. The differences between the global model with and without additional constraints from the regional model demonstrate that the Weddell Sea circulation exerts a significant influence on the course of the Antarctic Circumpolar Current with consequences for Southern Ocean water mass characteristics and the spreading of deep and bottom waters in the South Atlantic. The influence of the Ross Sea is found to be less important in terms of global influences. However, regional changes in the Pacific sector of the Southern Ocean are found to be of Ross Sea origin. The additional constraints change the hydrographic conditions of the global model in the vicinity of the Antarctic Circumpolar Current in such a way that transport values, e.g., in Drake Passage no longer need to be prescribed to obtain observed transports. These changes not only improve the path and transport of the Antarctic Circumpolar Current but affect the meso- and large-scale circulation. With a higher (lower) mean Drake Passage transport, the mean Weddell Gyre transport is lower (higher). Furthermore, an increase (decrease) in the Antarctic Circumpolar Current leads to a decrease (increase) of the circum-Australian flow, i.e., a decrease (increase) of the Indonesian Throughflow.     

The equatorial dynamics of a shallow,homogeneous ocean

Abstract

It is shown that the linear equatorial dynamics of a shallow ocean is characterized by two boundary layers of width γ^? L and γL (γ is the Ekman number of the flow, assumed small, and L is a horizontal dimension of the basin). In the γ^? layer stress in the bottom Ekman layer is comparable to that in the surface Ekman layer. In the γ layer vertical friction is important throughout the depth of the ocean. Should the Rossby number ? be so large as to invalidate a linear theory (? > γ^5/3), then inertial effects become important at a distance ?^2/5 L from the equator. The role played in the circulation of the basin by the non-linear equatorial current first studied by Charney (1960) is shown to be similar to that of the γ layer of the linear theory. Though lateral friction is unimportant in a linear model of the flow, shear layers at the equator are found to be a necessary feature of non-linear flow.     

Multi-scale modeling for dynamics of structure-soil-structure interactions

This paper presents two cases of multi-scale modeling in the context of dynamic structure soil-structure interaction. The first case concerns the behaviour of reinforced soils. It is shown that such system may involve both shear and bending effect at the leading order, which corresponds to a second gradient material. The second case addresses the seismic response of soils in presence of a densely urbanized city. It appears that the effect of resonance of the whole buildings in interaction actually modify the seismic response. In both cases the theoretical approach is completed by a validation through analogous samples tested on a shaking table.     

Diurnal-period,wind-forced ocean variability on the inner shelf off Concepción,Chile

We characterize the response of diurnal-period ocean current variability to the sea breeze using measurements of current velocity taken off the mouth of the Itata River and wind stress collected at Hualpen Point (central Chile) in spring of 2007 and summer of 2006 and 2008. During these three periods, the winds are predominately towards the northeast, following the coastal topography, with the highest variability found in the near-diurnal and synoptic frequency bands. The sea breeze amplitude is intermittent in time and is associated with synoptic-scale variability on the order of three to 15 days, so that the diurnal-period winds (and currents) are enhanced when the alongshore wind (i.e. upwelling-favorable) is strong. The water current variability in the near-diurnal band is significant, explaining up to 40% (spring 2007) of the total current variance in the first 15 m depth.     

Multi-scale support vector algorithms for hot spot detection and modelling

The algorithmic approach to data modelling has developed rapidly these last years, in particular methods based on data mining and machine learning have been used in a growing number of applications. These methods follow a data-driven methodology, aiming at providing the best possible generalization and predictive abilities instead of concentrating on the properties of the data model. One of the most successful groups of such methods is known as Support Vector algorithms. Following the fruitful developments in applying Support Vector algorithms to spatial data, this paper introduces a new extension of the traditional support vector regression (SVR) algorithm. This extension allows for the simultaneous modelling of environmental data at several spatial scales. The joint influence of environmental processes presenting different patterns at different scales is here learned automatically from data, providing the optimum mixture of short and large-scale models. The method is adaptive to the spatial scale of the data. With this advantage, it can provide efficient means to model local anomalies that may typically arise in situations at an early phase of an environmental emergency. However, the proposed approach still requires some prior knowledge on the possible existence of such short-scale patterns. This is a possible limitation of the method for its implementation in early warning systems. The purpose of this paper is to present the multi-scale SVR model and to illustrate its use with an application to the mapping of Cs¹³⁷ activity given the measurements taken in the region of Briansk following the Chernobyl accident.     

Stable and unstable barotropic shelf waves in a coastal current

Abstract

We consider the linearized stability of a barotropic coastal current flowing parallel to a straight coastline over a continental shelf and slope whose depth varies monotonically with distance from the coast. Some necessary conditions for stability and various semi-circle theorems are reviewed for general current profiles and bottom topography. A criterion for topography to be a destabilizing influence is derived. Some general results for stable waves are also described. Analytic solutions are obtained for a piece-wise linear current profile and the exponential depth profile (Buchwald and Adams, 1968). Dispersion diagrams are obtained for a monotonic current profile, where it is shown that the effect of topography is destabilizing, and for a triangular current profile. The dispersion diagrams generally contain a finite number (usually one or two) of unstable waves, and a set of stable waves, which may be infinite in number. The results are applied to some specific coastal regimes.     

The importance of ocean temperature to global biogeochemistry

Variations in the mean temperature of the ocean, on time scales from millennial to millions of years, in the past and projected for the future, are large enough to impact the geochemistry of the carbon, oxygen, and methane geochemical systems. In each system, the time scale of the temperature perturbation is key. On time frames of 1-100 ky, atmospheric CO₂ is controlled by the ocean. CO₂ temperature-dependent solubility and greenhouse forcing combine to create an amplifying feedback with ocean temperature; the CaCO₃ cycle increases this effect somewhat on time scales longer than ∼5-10 ky. The CO₂/T feedback can be seen in the climate record from Vostok, and a model including the temperature feedback predicts that 10% of the fossil fuel CO₂ will reside in the atmosphere for longer than 100 ky. Timing is important for oxygen, as well; the atmosphere controls the ocean on short time scales, but ocean anoxia controls atmospheric pO₂ on million-year time scales and longer. Warming the ocean to Cretaceous temperatures might eventually increase pO₂ by approximately 25%, in the absence of other perturbations. The response of methane clathrate to climate change in the coming century will probably be small, but on longer time scales of 1-10 ky, there may be a positive feedback with ocean temperature, amplifying the long-term climate impact of anthropogenic CO₂ release.     

Interactions between barotropic tides and mesoscale processes in deep ocean and shelf regions

Ocean Dynamics - The interactions between barotropic tides and mesoscale processes were studied using the results of a numerical model in which tidal forcing was turned on and off. The research...     

Perturbation of regional ocean tides due to coastal dikes

The tidal regime modeling system for ocean tides in the seas bordering the Korean Peninsula is designed to cover an area that is broad in scope and size, yet provide a high degree of resolution in coastal development areas, including the Saemangeum area in the eastern Yellow Sea and the Ariake Sea in Japan, where serious environmental problems have occurred after the completion of interior tidal dikes. With this simulation system, we have estimated the changes in tidal regime due to barriers at Saemangeum and Isahaya Bay in the Ariake Sea. Some results in terms of perturbations in tidal elevations due to the construction of coastal dikes are presented and discussed.     

Interactions between wind and tidally induced currents in coastal and shelf basins

This paper addresses the impact of atmospheric variability on ocean circulation in tidal and non-tidal basins. The data are generated by an unstructured-grid numerical model resolving the dynamics in the coastal area, as well as in the straits connecting the North Sea and Baltic Sea. The model response to atmospheric forcing in different frequency intervals is quantified. The results demonstrate that the effects of the two mechanical drivers, tides and wind, are not additive, yet non-linear interactions play an important role. There is a tendency for tidally and wind-driven circulations to be coupled, in particular in the coastal areas and straits. High-frequency atmospheric variability tends to amplify the mean circulation and modify the exchange between the North and the Baltic Sea. The ocean response to different frequency ranges in the wind forcing is area-selective depending on specific local dynamics. The work done by wind on the oceanic circulation depends strongly upon whether the regional circulation is tidally or predominantly wind-driven. It has been demonstrated that the atmospheric variability affects the spring-neap variability very strongly.     

Unstructured,anisotropic mesh generation for the Northwestern European continental shelf,the continental slope and the neighbouring ocean

A new mesh refinement strategy for generating high quality unstructured meshes of the Northwestern European continental shelf, the continental slope and the neighbouring ocean is presented. Our objective is to demonstrate the ability of anisotropic unstructured meshes to adequately address the challenge of simulating the hydrodynamics occurring in these three regions within a unique mesh. The refinement criteria blend several hydrodynamic considerations as the tidal wave propagation on the continental shelf and the hydrostatic consistency condition in steep areas. Several meshes illustrate both the validity and the efficiency of the refinement strategy. The selection of the refinement parameters is discussed. Finally, an attempt is made to take into account tidal ellipses, providing another cause for anisotropy in the mesh.