An idealized model of tidal dynamics in the North Sea: resonance properties and response to large-scale changes

An idealized process-based model is developed to investigate tidal dynamics in the North Sea. The model geometry consists of a sequence of different rectangular compartments of uniform depth, thus, accounting for width and depth variations in a stepwise manner. This schematization allows for a quick and transparent solution procedure. The solution, forced by incoming Kelvin waves at the open boundaries and satisfying the linear shallow water equations on the f plane with bottom friction, is in each compartment written as a superposition of eigenmodes, i.e. Kelvin and Poincaré waves. A collocation method is employed to satisfy boundary and matching conditions. First, the general resonance properties of a strongly simplified geometry with two compartments, representing the Northern North Sea and Southern Bight, are studied. Varying the forcing frequency while neglecting bottom friction reveals Kelvin and Poincaré resonance. These resonances continue to exist (but with lower amplification and a modified spatial structure) when adding the Dover Strait as a third compartment and separating the solutions due to forcing from either the north or the south only. Including bottom friction dampens the peaks. Next, comparison with tide observations along the North Sea coast shows remarkable agreement for both semi-diurnal and diurnal tides. This result is achieved with a more detailed geometry consisting of 12 compartments fitted to the coastline of the North Sea. Further simulations emphasize the importance of Dover Strait and bottom friction. Finally, it is found that a sea level rise of 1 m, uniformly applied to the entire North Sea, amplifies the M2-elevation amplitudes almost everywhere along the coast, with an increase of up to 8 cm in Dover Strait. Bed level changes of ±1 m, uniformly applied to the Southern Bight only, imply weaker changes, with changes in coastal M2-elevation amplitudes below 5 cm.     

海啸波近岸共振响应的数值模拟及分析

应用有限差分方法求解非线性浅水长波方程,建立了海啸波产生和传播的二维数值模型;对太平洋地震引起的夏威夷群岛海啸波进行模拟,并将模拟结果和测站实测值进行比较,验证了模型的正确性。利用快速傅里叶变换对数值模拟结果进行谱分析,得到整个计算区域的能量谱分布,并给出了发生能量聚集的位置及相应的谱峰周期。讨论了海啸波传播到近岸时可能产生的共振响应现象,发现海啸波和近岸的共振响应不仅与近岸复杂地形有关,还与海啸波传播到近岸时的波浪入射方向有关。     

Generation of baroclinic topographic waves by a tropical cyclone impacting a low-latitude continental shelf

Numerical model experiments have been performed to analyze the low-latitude baroclinic continental shelf response to a tropical cyclone. The theory of coastally trapped waves suggests that, provided appropriate slope, latitude, stratification and wind stress, bottom-intensified topographic Rossby waves can be generated by the storm. Based on a scale analysis, the Nicaragua Shelf is chosen to study propagating topographic waves excited by a storm, and a model domain is configured with simplified but similar geometry. The model is forced with wind stress representative of a hurricane translating slowly over the region at 6 km h⁻¹. Scale analysis leads to the assumption that baroclinic Kelvin wave modes have minimal effect on the low-frequency wave motions along the slope, and coastal-trapped waves are restricted to topographic Rossby waves. Analysis of the simulated motions suggests that the shallow part of the continental slope is under the influence of barotropic topographic wave motions and at the deeper part of the slope baroclinic topographic Rossby waves dominate the low-frequency motions. Numerical solutions are in a good agreement with theoretical scale analysis. Characteristics of the simulated baroclinic waves are calculated based on linear theory of bottom-intensified topographic Rossby waves. Simulated waves have periods ranging from 153 to 203 h. The length scale of the waves is from 59 to 87 km. Analysis of energy fluxes for a fixed volume on the slope reveals predominantly along-isobath energy propagation in the direction of the group velocity of a topographic Rossby wave. Another model experiment forced with a faster translating hurricane demonstrates that fast moving tropical cyclones do not excite energetic baroclinic topographic Rossby waves. Instead, robust inertial oscillations are identified over the slope.     

Modelling the influence of spatially varying hydrodynamics on the cross-sectional stability of double inlet systems

The cross-sectional stability of double inlet systems is investigated using an exploratory model that combines Escoffier’s stability concept for the evolution of the inlet’s cross-sectional area with a two-dimensional, depth-averaged (2DH) hydrodynamic model for tidal flow. The model geometry consists of four rectangular compartments, each with a uniform depth, associated with the ocean, tidal inlets and basin. The water motion, forced by an incoming Kelvin wave at the ocean’s open boundary and satisfying the linear shallow water equations on the f -plane with linearised bottom friction, is in each compartment written as a superposition of eigenmodes, i.e. Kelvin and Poincaré waves. A collocation method is employed to satisfy boundary and matching conditions. The analysis of resulting equilibrium configurations is done using flow diagrams.

Model results show that internally generated spatial variations in the water motion are essential for the existence of stable equilibria with two inlets open. In the hydrodynamic model used in the paper, both radiation damping into the ocean and basin depth effects result in these necessary spatial variations. Coriolis effects trigger an asymmetry in the stable equilibrium cross-sectional areas of the inlets. Furthermore, square basin geometries generally correspond to significantly larger equilibrium values of the inlet cross-sections. These model outcomes result from a competition between a destabilising (caused by inlet bottom friction) and a stabilising mechanism (caused by spatially varying local pressure gradients over the inlets).

     

硬化粗质地基冻土层的速度参数与卓越周期的关系

本文给出了对共振现象出现概率的理论分析与实验数据间关系的思考。样品取自位于半空间的硬化粗质冻土地基,包括石质冻层。参考从三个场地14次地震记录的纵、横波速求得的土层相对频率特征,进行了对其共振频率的评估。纵、横波速在前贝加尔和后贝加尔地区的粗质冻土层中随湿度和温度而变化。将这些概括化的速度值与一个共振频率和纵、横波度之间的关系式结合起来,用于评估被调查土层的纵、横波速的概率值。结果显示,通过综合在实验室尺度获得的速度值和现场测量获得的速度值的数据,则利用地震记录波来确定在松散土层内传播速度的概率值是可能的。     

热弹性介质中波传播特征

深层-超深层油气地震勘探涉及高温介质地震波传播问题,热弹介质参数对地震波传播有重要影响.含弛豫时间修正项的Lord-Shulman双曲型耦合热弹波动方程从理论上预测了热弹性介质中存在快纵波、慢纵波(一种准静态慢纵波,简称热波)和横波的传播,两个纵波为热耗散衰减波而横波不受介质热特性的影响.本文结合平面波频散分析和格林函数法数值模拟,详细研究两个热耗散衰减波的频散和衰减特征,着重分析热导率、热膨胀系数及比热的变化对波速和衰减的影响.研究表明热导率作为主要参数决定了波速与衰减的临界变化,热膨胀系数对波速和衰减的幅度有明显影响,比热则兼顾了前两个热弹系数的影响特征.最后,利用热弹性动力学频率域的二阶格林函数进行波场快照数值模拟,展示热弹性介质中纵波、横波和热波的传播行为.     

Response of large-scale coastal basins to wind forcing: influence of topography

Because wind is one of the main forcings in storm surge, we present an idealised process-based model to study the influence of topographic variations on the frequency response of large-scale coastal basins subject to time-periodic wind forcing. Coastal basins are represented by a semi-enclosed rectangular inner region forced by wind. It is connected to an outer region (represented as an infinitely long channel) without wind forcing, which allows waves to freely propagate outward. The model solves the three-dimensional linearised shallow water equations on the f plane, forced by a spatially uniform wind field that has an arbitrary angle with respect to the along-basin direction. Turbulence is represented using a spatially uniform vertical eddy viscosity, combined with a partial slip condition at the bed. The surface elevation amplitudes, and hence the vertical profiles of the velocity, are obtained using the finite element method (FEM), extended to account for the connection to the outer region. The results are then evaluated in terms of the elevation amplitude averaged over the basin’s landward end, as a function of the wind forcing frequency. In general, the results point out that adding topographic elements in the inner region (such as a topographic step, a linearly sloping bed or a parabolic cross-basin profile), causes the resonance peaks to shift in the frequency domain, through their effect on local wave speed. The Coriolis effect causes the resonance peaks associated with cross-basin modes (which without rotation only appear in the response to cross-basin wind) to emerge also in the response to along-basin wind and vice versa.     

A reduced model for nonlinear dispersive waves in a rotating environment

Abstract

The simplest model for geophysical flows is one layer of a constant density fluid with a free surface, where the fluid motions occur on a scale in which the Coriolis force is significant. In the linear shallow water limit, there are non-dispersive Kelvin waves, localized near a boundary or near the equator, and a large family of dispersive waves. We study weakly nonlinear and finite depth corrections to these waves, and derive a reduced system of equations governing the flow. For this system we find approximate solitary Kelvin waves, both for waves traveling along a boundary and along the equator. These waves induce jets perpendicular to their direction of propagation, which may have a role in mixing. We also derive an equivalent reduced system for the evolution of perturbations to a mean geostrophic flow.     

Kelvin wave generation by a semi-infinite barrier

A closed form solution is given for the problem of simple harmonic plane waves incident obliquely on a semi-infinite barrier. The amplitude of the Kelvin wave component is considered and it is found that increased amplitudes of the Kelvin waves can be expected for certain wave numbers, depending upon the angle of incidence.     

Propagation of Kelvin waves along irregular coastlines in finite-difference models

In this paper, we examine the behavior of internal Kelvin waves on an f-plane in finite-difference models using the Arakawa C-grid. The dependence of Kelvin wave phase speed on offshore grid resolution and propagation direction relative to the numerical grid is illustrated by numerical experiments for three different geometries: (1) Kelvin wave propagating along a straight coastline; (2) Kelvin wave propagating at a 45° angle to the numerical grid along a stairstep coastline with stairstep size equal to the grid spacing; (3) Kelvin wave propagating at a 45° angle to the numerical grid along a coarse resolution stairstep coastline with stairstep size greater than the grid spacing. It can be shown theoretically that the phase speed of a Kelvin wave propagating along a straight coastline on an Arakawa C-grid is equal to the analytical inviscid wave speed and is not dependent on offshore grid resolution. However, we found that finite-difference models considerably underestimate the Kelvin wave phase speed when the wave is propagating at an angle to the grid and the grid spacing is comparable with the Rossby deformation radius. In this case, the phase speed converges toward the correct value only as grid spacing decreases well below the Rossby radius. A grid spacing of one-fifth the Rossby radius was required to produce results for the stairstep boundary case comparable with the straight coast case. This effect does not appear to depend on the resolution of the coastline, but rather on the direction of wave propagation relative to the grid. This behavior is important for modeling internal Kelvin waves in realistic geometries where the Rossby radius is often comparable with the grid spacing, and the waves propagate along irregular coastlines.©1998 Published by Elsevier Science Limited. All rights reserved     

Tsunami–tide interactions: A Cook Inlet case study

First, we investigated some aspects of tsunami–tide interactions based on idealized numerical experiments. Theoretically, by changing total ocean depth, tidal elevations influence the speed and magnitude of tsunami waves in shallow regions with dominating tidal signals. We tested this assumption by employing a simple 1-D model that describes propagation of tidal waves in a channel with gradually increasing depth and the interaction of the tidal waves with tsunamis generated at the channel's open boundary. Important conclusions from these studies are that computed elevations by simulating the tsunami and the tide together differ significantly from linear superposing of the sea surface heights obtained when simulating the tide and the tsunami separately, and that maximum tsunami–tide interaction depends on tidal amplitude and phase. The major cause of this tsunami–tide interaction is tidally induced ocean depth that changes the conditions of tsunami propagation, amplification, and dissipation. Interactions occur by means of momentum advection, bottom friction, and variable water flux due to changing total depth and velocity. We found the major cause of tsunami–tide interactions to be changing depth. Secondly, we investigate tsunami–tide interactions in Cook Inlet, Alaska, employing a high-resolution 2-D numerical model. Cook Inlet has high tides and a history of strong tsunamis and is a potential candidate for tsunami impacts in the future. In agreement with previous findings, we find that the impacts of tsunamis depend on basin bathymetries and coastline configurations, and they can, in particular, depend on tsunami–tide interactions. In regions with strong tides and tsunamis, these interactions can result in either intensification or damping of cumulative tsunami and tide impacts, depending on mean basin depth, which is regulated by tides. Thus, it is not possible to predict the effect of tsunami–tide interaction in regions with strong tides without making preliminary investigations of the area. One approach to reduce uncertainties in tsunami impact in regions with high tides is to simulate tsunamis together with tidal forcing.     

Modelling the influence of spatially varying hydrodynamics on the cross-sectional stability of double inlet systems

The cross-sectional stability of double inlet systems is investigated using an exploratory model that combines Escoffier’s stability concept for the evolution of the inlet’s cross-sectional area with a two-dimensional, depth-averaged (2DH) hydrodynamic model for tidal flow. The model geometry consists of four rectangular compartments, each with a uniform depth, associated with the ocean, tidal inlets and basin. The water motion, forced by an incoming Kelvin wave at the ocean’s open boundary and satisfying the linear shallow water equations on the f -plane with linearised bottom friction, is in each compartment written as a superposition of eigenmodes, i.e. Kelvin and Poincaré waves. A collocation method is employed to satisfy boundary and matching conditions. The analysis of resulting equilibrium configurations is done using flow diagrams. Model results show that internally generated spatial variations in the water motion are essential for the existence of stable equilibria with two inlets open. In the hydrodynamic model used in the paper, both radiation damping into the ocean and basin depth effects result in these necessary spatial variations. Coriolis effects trigger an asymmetry in the stable equilibrium cross-sectional areas of the inlets. Furthermore, square basin geometries generally correspond to significantly larger equilibrium values of the inlet cross-sections. These model outcomes result from a competition between a destabilising (caused by inlet bottom friction) and a stabilising mechanism (caused by spatially varying local pressure gradients over the inlets).     

Seismoelectric numerical simulation in 2D vertical transverse isotropic poroelastic medium

Seismoelectric coupling in an electric isotropic and elastic anisotropic medium is developed using a primary–secondary formulation. The anisotropy is of vertical transverse isotropic type and concerns only the poroelastic parameters. Based on our finite difference time domain algorithm, we solve the seismoelectric response to an explosive source. The seismic wavefields are computed as the primary field. The electric field is then obtained as a secondary field by solving the Poisson equation for the electric potential. To test our numerical algorithm, we compared our seismoelectric numerical results with analytical results obtained from Pride's equation. The comparison shows that the numerical solution gives a good approximation to the analytical solution. We then simulate the seismoelectric wavefields in different models. Simulated results show that four types of seismic waves are generated in anisotropic poroelastic medium. These are the fast and slow longitudinal waves and two separable transverse waves. All of these seismic waves generate coseismic electric fields in a homogenous anisotropic poroelastic medium. The tortuosity has an effect on the propagation of the slow longitudinal wave. The snapshot of the slow longitudinal wave has an oval shape when the tortuosity is anisotropic, whereas it has a circular shape when the tortuosity is isotropic. In terms of the Thomsen parameters, the radiation anisotropy of the fast longitudinal wave is more sensitive to the value of ε, while the radiation anisotropy of the transverse wave is more sensitive to the value of δ.     

简谐平面弹性波在楔形域中传播

本文利用自动满足楔形域任一面自由边界条件的纵、横驻波系,满足另一面上的初始应力条件,求得了初始简谐平面P波入射时楔形域中波传播的解析解.     

Inertia-gravity waves in the troposphere and lower stratosphere associated with a jet stream exit region

Radar measurements at Aberystwyth (52.4°N, 4.1°W) of winds at tropospheric and lower stratospheric heights are shown for 12–13 March 1994 in a region of highly curved flow, downstream of the jet maximum. The perturbations of horizontal velocity have comparable amplitudes in the troposphere and lower stratosphere with downward and upward phase propagation, respectively, in these two height regions. The sense of rotation with increasing height in hodographs of horizontal perturbation velocity derived for hourly intervals show downwards propagation of energy in the troposphere and upward propagation in the lower stratosphere with vertical wavelengths of 1.7 to 2.3 km. The results indicate inertia-gravity waves propagating in a direction similar to that of the jet stream but at smaller velocities. Some of the features observed contrast with those of previous observations of inertia-gravity waves propagating transverse to the jet stream. The interpretation of the hodographs to derive wave parameters has taken account of the vertical shear of the background wind transverse to the direction of wave propagation.     

On the generation and propagation of internal solitary waves in the southern Andaman Sea: A numerical study

Numerous internal solitary waves(ISWs) have been observed in the southern Andaman Sea. In this study, the two-dimensional Massachusetts Institute of Technology general circulation model is applied to investigate the dynamics of ISWs and explore the effects of the bottom topography and tidal forcing on the generation and propagation of ISWs in the southern Andaman Sea. The results show that the large-amplitude depression ISWs are mainly generated via the oscillating tidal flow over the sill of the Great Channel, and the generation of ISWs is subject to the lee wave regime. The Dreadnought Bank cannot generate ISWs itself; however, it can enhance the amplitudes of eastward-propagating ISWs generated from sill A, owing to constructive interference of internal tide generation between the sill of the Great Channel and the Dreadnought Bank. The eastward-propagating ISWs generated by the eastern shallow sill near the continental shelf can propagate to the shelf, where they evolve into elevation waves because of the shallow water. Sensitivity runs show that both the semidiurnal and diurnal tides over the sill of the Great Channel can generate ISWs in this area. However, the ISWs generated by diurnal tides are much weaker than those generated by semidiurnal tides. Mixed tidal forcing has no significant effect on the generation of ISWs.     

Wave transformation across a macrotidal shore platform under low to moderate energy conditions

We investigate how waves are transformed across a shore platform as this is a central question in rock coast geomorphology. We present results from deployment of three pressure transducers over four days, across a sloping, wide (~200 m) cliff‐backed shore platform in a macrotidal setting, in South Wales, United Kingdom. Cross‐shore variations in wave heights were evident under the predominantly low to moderate (significant wave height < 1.4 m) energy conditions measured. At the outer transducer 50 m from the seaward edge of the platform (163 m from the cliff) high tide water depths were 8+ m meaning that waves crossed the shore platform without breaking. At the mid‐platform position water depth was 5 m. Water depth at the inner transducer (6 m from the cliff platform junction) at high tide was 1.4 m. This shallow water depth forced wave breaking, thereby limiting wave heights on the inner platform. Maximum wave height at the middle and inner transducers were 2.41 and 2.39 m, respectively, and significant wave height 1.35 m and 1.34 m, respectively. Inner platform high tide wave heights were generally larger where energy was up to 335% greater than near the seaward edge where waves were smaller. Infragravity energy was less than 13% of the total energy spectra with energy in the swell, wind and capillary frequencies accounting for 87% of the total energy. Wave transformation is thus spatially variable and is strongly modulated by platform elevation and the tidal range. While shore platforms in microtidal environments have been shown to be highly dissipative, in this macro‐tidal setting up to 90% of the offshore wave energy reached the landward cliff at high tide, so that the shore platform cliff is much more reflective. Copyright © 2017 John Wiley & Sons, Ltd.     

Fluid transport by solitary waves along growing faults: A field example from the South Eugene Island Basin, Gulf of Mexico

The Red Fault system is one of the main growth faults found in the South Eugene Island Basin, a salt withdrawal minibasin located offshore Louisiana, in the Gulf of Mexico. This fault system corresponds to a lateral boundary between fluid overpressured compartments. In addition, there is a set of observations indicating that the Red Fault system exhibits rapid episodic migration of fluids. This fault represents an example of preferential pathway for the upward episodic migration of overpressured hydrocarbons from deep, heavily pressured, compartments on time scales of years. The migrations of fluids into active growing faults could take the form of propagating surges (solitary waves) that propagate upward along the fault planes in a wave-like manner at km/yr. Solitary waves represent a very efficient mechanism for the upward transport of fluids along growth faults in sedimentary basins generating its own permeability. In addition, this mechanism is compatible with the fact that the fault plane is observed to sustain a static pore fluid pressure difference between its two sides. The propagation of solitary waves in active growth faults appears as a fundamental mechanism to understand the nature of upward fast migration of fluids along active growth faults in compartimentalized sedimentary basins.     

地震波在热弹性介质中的衰减性质及其相应Q值的计算方法

该文得出热弹性介质中波的基本形式有准弹性纵波、横波和具有有限传播速度的温度波,并指出横波在传播过程中同其在完全弹性介质中一样是不衰减的,而准弹性纵波和温度波是相互伴随的、衰减的。其次,给出了具有两种不同传播速度的准弹性纵波的基本形式。最后,将地震P波看作由初始位移扰动的准弹性纵波,就固定频率和固定波长两种情况分别给出了相应Qr值的计算公式和近似公式。通常人们是利用多台的地震波资料在空间域中计算Qr值,而本文得到的时间域中Qr值的近似公式,为利用单台地震P波资料计算Qr值提供了一种方法。     

The wave field patterns of the propagation of longitudinal and transverse elastic waves in grain-oriented rocks

In the work under consideration, on the basis of data analysis on the velocities of quasi-longitudinal Vp and quasi-transverse Vs waves—measured for rock samples at high pressures and temperatures, and analogous velocities, calculated on the basis of the grain orientation distribution function, reconstructed from the neutron-diffraction textural experiment—the indicative inconsistencies between the experimental and model characteristics were inferred. The theoretical analysis of the wave field patterns of the propagation of longitudinal and transverse elastic waves in the anisotropic media is carried out. It is established that, in the general case, in the anisotropic inhomogeneous media the velocities of Vp and Vs propagation, measured experimentally and obtained from the modeling, cannot coincide due to the existence of the physical coupling between the vibrations of two types: transverse and longitudinal vibrations.