Interactions between tides and other frequencies in the Indonesian seas

Interactions of tidal constituents and the transfer of energy from the tidal frequencies to other frequencies are investigated using 3-D tidal simulations for the Indonesian seas, focusing on an area of active internal tides. Semidiurnal tides strongly affect diurnal tides; however, semidiurnal tides are essentially unaffected by diurnal tides. The semidiurnal and diurnal constituents interact with each other through non-linear interference, both destructive and constructive. Semidiurnal tides generate harmonics at nearly the diurnal frequency and higher vertical wavenumbers. In Ombai Strait, these harmonics are out of phase with the diurnal tides and interact destructively with the diurnal tides, effectively negating the diurnal response in some locations. However, this is not a general response, and interactions differ between locations. Energy is also transferred from both semidiurnal and diurnal tides to other frequencies across the spectrum, with more energy originating from semidiurnal tides. These energy transfers are not homogeneous, and the spectral responses differ between the Makassar and Ombai Straits, with the region east of Ombai showing a more active surface response compared to a more intense benthic response in Makassar. In deep water away from topography, velocity spectra generally follow the Garrett–Munk (GM) relation. However, in areas of internal tide generation, spectral density levels exceed GM levels, particularly between 4 and 8 cycles per day (cpd), indicating increased non-linear interactions and energy transfer through resonant interactions. The model indicates strong surface trapping of internal tides, with surface velocity spectra having significantly higher energy between 4 and 8 cpd even 100 km away from the prominent sill generating the internal tides.     

A model study of tidal distributions in the Celtic and Irish Sea regions determined with finite volume and finite element models

An unstructured mesh tidal model of the west coast of Britain, covering the Celtic Sea and Irish Sea is used to compare tidal distributions computed with finite element (FE) and finite volume (FV) models. Both models cover an identical region, use the same mesh, and have topography and tidal boundary forcing from a finite difference model that can reproduce the tides in the region. By this means, solutions from both models can be compared without any bias towards one model or another. Two-dimensional calculations show that for a given friction coefficient, there is more damping in the FV model than the FE model. As bottom friction coefficient is reduced, the two models show comparable changes in tidal distributions. In terms of mesh resolution, calculations show that for the M₂ tide, the mesh is sufficiently fine to yield an accurate solution over the whole domain. However, in terms of higher harmonics of the tide, in particular the M₆ component, its small-scale variability in near-shore regions which is comparable to the mesh of the model, suggests that the mesh resolution is insufficient in the near-coastal regions. Even with a finer mesh in these areas, without detailed bottom topography and a spatial varying friction depending on bed types and bed forms, which is not available, model skill would probably not be improved. In addition in the near-shore region, as shown in the literature, the solution is sensitive to the form of the wetting/drying algorithm used in the model. Calculations with a 3D version of the FV model show that for a given value of k, damping is reduced compared to the 2D version due to the differences in bed stress formulation, with the 3D model yielding an accurate tidal distribution over the region.     

Experiments with a numerical model of coastal currents and tidal mixing fronts

A three-dimensional numerical model is used to simulate the development of disturbances on shelf-sea coastal currents and fronts. The model, which has a free surface, uses a finite difference grid ☐ scheme based on sigma coordinates. It has a semi-implicit scheme for the barotropic flow and a hydrid advection scheme to retain sharp fronts. The results demonstrate that (i) eddy formation follows changes at the inflow of a coastal current, (ii) a simple radiation boundary condition at the outflow produces nearly identical results for different outflow boundary positions, (iii) eddy growth, with matching behaviour of surface and bottom fronts, follows a small displacement on a tidal mixing front and (iv) effects of friction and mixing can significantly alter the behaviour of the front and the relative strength of the cyclonic and anticyclonic eddies formed.     

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.     

A 3-D finite volume model for sediment transport in coastal waters

Ocean Dynamics - A 3-D model has been developed to simulate sediment transport and bed change induced by currents and waves in coastal waters. The currents are calculated with the 3-D...     

Development of a 2-km resolution ocean model covering the coastal seas around Japan for operational application

Ocean Dynamics - In order to expand the coastal ocean monitoring and forecasting system of the Japan Meteorological Agency from the Seto Inland Sea to the entire coastal seas of Japan, a 2-km...     

A numerical study of the string function using a primitive equation ocean model

Abstract

We use results from a primitive-equation ocean numerical model (SCRUM) to test a theoretical 'string function' formulation put forward by Tyler and Käse in another article in this issue. The string function acts as a stream function for the large-scale potential energy flow under the combined beta and topographic effects. The model results verify that large-scale anomalies propagate along the string function contours with a speed correctly given by the cross-string gradient. For anomalies having a scale similar to the Rossby radius, material rates of change in the layer mass following the string velocity are balanced by material rates of change in relative vorticity following the flow velocity. It is shown that large-amplitude anomalies can be generated when wind stress is resonant with the string function configuration.     

海潮误差对GRACE时变重力场解算的影响研究

海潮误差是 GRACE 时变重力场反演中重要的误差源,目前发布的海潮模型中主要包含振幅较大的主潮波分量模型,在时变重力场反演中次潮波的影响也是不可忽略的,因此,GRACE 时变重力场反演中的海潮误差主要包括受限于海潮模型误差和次潮波影响.本文利用轨道模拟方法检测了短周期潮波的混频周期以及次潮波对ΔC20, ΔC30的时序特征,并进一步通过轨道模拟结果分析了海潮误差对时变重力场反演的影响,然后通过实测数据解算分析了海潮误差对当前 GRACE 时变重力场解算的影响,研究发现:(1) 利用轨道模拟能够有效地检测短周期潮波的混频周期;(2)时变重力场解算过程中,次潮波的影响大于海潮模型误差的影响;(3)海潮模型误差以及次潮波影响是当前 GRACE 没有达到基准精度的重要因素之一.     

A numerical study of equilibrium states in tidal network morphodynamics

The long-term morphodynamic evolution of tidal networks on tidal flats is investigated using a two-dimensional numerical model. We explore the physical processes related to the development of the morphology and the presence of equilibrium configurations. Tidal networks are simulated over a rectangular domain representing a tidal platform, a different setting compared to estuaries (subject to riverine influence) and lagoons (offshore bars constricting the flow). In the early and middle phases of the tidal network evolution, large sediment patches with rhombus-like shape form and gradually migrate in the flood direction, even though the overall sediment flux is ebb-directed. A cross-section-averaged “equilibrium” state is asymptotically approached after about 500 years. The area and peak discharge of the lower flat cross-sections at year 500 approximately show a 1:1 relationship, which is in agreement with field observations. We also show that model results are consistent with the Q-A relationship (peak discharge Q versus cross-sectional area A), which is obtained under the assumption of a constant Chézy friction.     

A numerical study of the effect of the marginal sea on coastal upwelling in a non-linear inertial model

Inertia theory and the finite element method are used to investigate the effect of marginal seas on coastal upwelling. In contrast to much previous research on wind-driven upwelling, this paper does not consider localized wind effects, but focuses instead on temperature stratification, the slope of the continental shelf, and the background flow field. Finite element method, which is both faster and more robust than finite difference method in solving problems with complex boundary conditions, was developed to solve the partial differential equations that govern coastal upwelling. Our results demonstrate that the environment of the marginal sea plays an important role in coastal upwelling. First, the background flow at the outer boundary is the main driving force of upwelling. As the background flow strengthens, the overall velocity of cross-shelf flow increases and the horizontal scale of the upwelling front widens, and this is accompanied by the movement of the upwelling front further offshore. Second, temperature stratification determines the direction of cross-shelf flows, with strong stratification favoring a narrow and intense upwelling zone. Third, the slope of the continental shelf plays an important role in controlling the intensity of upwelling and the height that upwelling may reach: the steeper the slope, the lower height of the upwelling. An additional phenomenon that should be noted is upwelling separation, which occurs even without a local wind force in the nonlinear model.     

A numerical study of horizontal dispersion in a macro tidal basin

Tidal circulation in Cobscook Bay, a macro tidal basin, is simulated using the three-dimensional, nonlinear, finite element ocean model, QUODDY_dry. Numerical particles are released from various transects in the bay at different tidal phases and tracked for several tidal cycles. Initially, nearby particles in the main tidal channel experience a great deal of spreading and straining, and after a few tidal cycles, they are separated in different parts of the bay. The fundamental mechanism for particle dispersion is the chaotic advection that arises from long tidal excursions passing through many residual eddies. A loosely correlated, inverse relationship between the two dimensionless parameters, ν (the ratio of the residual current to the tidal current) and λ (the ratio of the tidal excursion to the main topographic scale), can be constructed for large values of ν. Several Lagrangian statistical measures are used to quantify and distinguish dispersion regimes in different parts of Cobscook Bay. It is found that the effective Lagrangian dispersion coefficient can be estimated using the product of the magnitude of residual currents and the tidal excursion.     

Far-field effects of tidal energy extraction in the Minas Passage on tidal circulation in the Bay of Fundy and Gulf of Maine using a nested-grid coastal circulation model

The Bay of Fundy in eastern Canada has the highest tides in the world. Harnessing the tidal energy in the region has long been considered. In this study, the effects of tidal in-stream energy extraction in the Minas Passage on the three-dimensional (3D) tidal circulation in the Bay of Fundy (BoF) and the Gulf of Maine (GoM) are examined using a nested-grid coastal ocean circulation model based on the Princeton Ocean Model (POM). The nested-grid model consists of a coarse-resolution (~4.5 km) parent sub-model for the GoM and a high-resolution (~1.5 km) child sub-model for the BoF. The tidal in-stream energy extraction in the model is parameterized in terms of nonlinear Rayleigh friction in the momentum equation. A suite of numerical experiments are conducted to determine the ranges of extractable tidal in-stream energy and resulting effects on the 3D tidal circulation over the Bay of Fundy and the Gulf of Maine (BoF-GoM) in terms of the Rayleigh friction coefficients. The 3D model results suggest that the maximum energy extraction in the Minas Passage increases tidal elevations and tidal currents throughout the GoM and reduces tidal elevations and circulation in the upper BoF, especially in the Minas Basin. The far-field effect of tidal energy extraction in the Passage on the 3D tidal circulation in the BoF-GoM is examined in two cases of harnessing tidal in-stream energy from (a) the entire water column and (b) the lower water column within 20 m above the bottom in the Passage. The 3D model results demonstrate that tidal in-stream energy extraction from the lower water column has less impact on the tidal elevations and circulation in the BoF-GoM than the energy extraction from the whole water column in the Minas Passage.     

Maintenance of the mean kinetic energy in the global ocean by the barotropic and baroclinic energy routes: the roles of JEBAR and Ekman dynamics

In order to determine the maintenance mechanisms of the currents of the global ocean, this study investigates the budget of the annual mean kinetic energy (KE) in a high-resolution (0.1° × 0.1°) semi-global ocean simulation. The analysis is based on a separation of the mean KE using the barotropic (i.e., depth-averaged) and baroclinic (the residual) components of velocity. The barotropic and baroclinic KEs dominate in higher and lower latitudes, respectively, with their global average being comparable to each other. The working rates of wind forcing on the barotropic and baroclinic circulations in the global ocean are 243 and 747 gigawatts, respectively. This study presents at least three new results for the budget of the barotropic KE. Firstly, an energy diagram is rederived to show that the work of the barotropic component of the horizontal pressure gradient (HPG) is connected to the work related to the joint effect of baroclinicity and bottom relief (JEBAR), and then to the budget of potential energy (PE). Secondly, the model analysis shows that the globally averaged work of the barotropic HPG (which is connected to the work related to JEBAR and then to the budget of the PE) is nearly zero. This indicates that the wind- and buoyancy-induced barotropic circulations in the global ocean are of the same strength with opposite sign. Thirdly, it is found that the work of the wind forcing on the barotropic component of the simulated Antarctic Circumpolar Current (ACC) is canceled by the combined effect, in equal measure, of the work of the barotropic HPG and the work of dissipative processes for mean KE. This result makes a significant contribution to the discussion on the depth-integrated momentum balance of the ACC. The barotropic KE is dissipated by the effects of bottom frictional stress, lateral frictional stress, and the Reynolds stress, of which more than half is attributed to an unexpectedly large contribution from biharmonic horizontal friction. Future studies should pay more attention to the role of biharmonic friction used in high-resolution numerical models.     

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.     

Analytical and numerical analysis of tides and salinities in estuaries; part I: tidal wave propagation in convergent estuaries

Analytical solutions of the momentum and energy equations for tidal flow are studied. Analytical solutions are well known for prismatic channels but are less well known for converging channels. As most estuaries have a planform with converging channels, the attention in this paper is fully focused on converging tidal channels. It will be shown that the tidal range along converging channels can be described by relatively simple expressions solving the energy and momentum equations (new approaches). The semi-analytical solution of the energy equation includes quadratic (nonlinear) bottom friction. The analytical solution of the continuity and momentum equations is only possible for linearized bottom friction. The linearized analytical solution is presented for sinusoidal tidal waves with and without reflection in strongly convergent (funnel type) channels. Using these approaches, simple and powerful tools (spreadsheet models) for tidal analysis of amplified and damped tidal wave propagation in converging estuaries have been developed. The analytical solutions are compared with the results of numerical solutions and with measured data of the Western Scheldt Estuary in the Netherlands, the Hooghly Estuary in India and the Delaware Estuary in the USA. The analytical solutions show surprisingly good agreement with measured tidal ranges in these large-scale tidal systems. Convergence is found to be dominant in long and deep-converging channels resulting in an amplified tidal range, whereas bottom friction is generally dominant in shallow converging channels resulting in a damped tidal range. Reflection in closed-end channels is important in the most landward 1/3 length of the total channel length. In strongly convergent channels with a single forward propagating tidal wave, there is a phase lead of the horizontal and vertical tide close to 90^o, mimicking a standing wave system (apparent standing wave).