Local generation of internal solitary waves in an oceanic pycnocline

Similar documents
Numerical simulations of the local generation of internal solitary waves in the Bay of Biscay

Conditions for local generation of NonLinear Internal Waves in a pycnocline: a numerical study.

Generation of internal solitary waves in a pycnocline by an internal wave beam: anumericalstudy

Nonlinear properties of internal gravity waves Part II : non-homogeneous medium

Reflecting tidal wave beams and local generation of solitary waves in the ocean thermocline

arxiv: v1 [nlin.ps] 20 May 2015

INTERNAL GRAVITY WAVES

Internal Solitary Waves in the Mozambique Channel: Observations. and Interpretation

Laboratory Modeling of Internal Wave Generation in Straits

Journal of Marine Systems

Hepeng Zhang Ben King Bruce Rodenborn

Generation and Propagation of Internal Solitary Waves on the Continental Shelf and Slope

Generation and Propagation of Internal Solitary Waves on the Continental Shelf and Slope

Chapter 14 Internal Solitary Waves System in the Mozambique Channel

The Luzon Strait experiment

Internal Wave Generation and Scattering from Rough Topography

Fission of a weakly nonlinear interfacial solitary wave at a step

Topographic Effects on Stratified Flows

Laboratory experiments on the generation of internal tidal beams over steep slopes

Generation of Internal Tides and Internal Solitary Waves on the Continental Shelf

Continuously stratified nonlinear low-mode internal tides

A laboratory study of low-mode internal tide scattering by finite-amplitude topography

University of Bristol - Explore Bristol Research. Link to publication record in Explore Bristol Research PDF-document.

Shear instabilities in a tilting tube

Goals of this Chapter

Application of an internal tide generation model to baroclinic spring-neap cycles

Non-Reflecting Internal Wave Beam Propagation in the Deep Ocean

Internal Waves. Thierry Dauxois

Generation and Evolution of Internal Waves in Luzon Strait

Tide-Topography Interactions: Asymmetries in Internal Wave Generation due to Surface Trapped Currents

Internal solitary waves induced by flow over a ridge: With applications to the northern South China Sea

The effect of a background shear current on large amplitude internal solitary waves

Using Ocean Acoustics to Improve Large Shallowwater Soliton Simulations

Instability of a coastal jet in a two-layer model ; application to the Ushant front

Ocean self-attraction and loading (SAL) and internal tides dissipation implementation within an unstructured global tide-surge model

MARINE RESEARCH. Journal of. Continuously stratified nonlinear low-mode internal tides. Volume 66, Number 3. by Karl R. Helfrich 1 ABSTRACT

Undular Bores and the Morning Glory

A Statistical Investigation of Internal Wave Propagation in the Northern South China Sea

A process study of tidal mixing over rough topography

Internal Waves in the Vicinity of the Kuroshio Path

On the interaction between internal tides and wind-induced near-inertial currents at the shelf edge

Scattering of Semidiurnal Internal Kelvin Wave at Step Bottom Topography

Nonlinear internal wave penetration via parametric subharmonic instability

Lecture 11: Internal solitary waves in the ocean

Kinematic Effects of Differential Transport on Mixing Efficiency in a Diffusively Stable, Turbulent Flow

INTERNAL SOLITARY WAVES IN THE ATMOSPHERE AND OCEAN

Flows Induced by 1D, 2D and 3D Internal Gravity Wavepackets

Tide-topography interaction?

A Model Study of Internal Tides in Coastal Frontal Zone*

FOURTH ORDER NONLINEAR EVOLUTION EQUATIONS FOR GRAVITY-CAPILLARY WAVES IN THE PRESENCE OF A THIN THERMOCLINE IN DEEP WATER

On the Calculation of Available Potential Energy in Internal Wave Fields

Stability of meridionally-flowing grounded abyssal currents in the ocean

Generation o f internal solitary waves in a pycnocline by an internal wave beam: a numerical study

Developing a nonhydrostatic isopycnalcoordinate

Resonant excitation of trapped coastal waves by free inertia-gravity waves

Refraction and Shoaling of Nonlinear Internal Waves at the Malin Shelf Break*

Internal Wave Generation in Straits

Internal solitary wave transformation over the bottom step: loss of energy

Nonlinear Disintegration of the Internal Tide

Limiting amplitudes of fully nonlinear interfacial tides and solitons

The Evolution of Large-Amplitude Internal Gravity Wavepackets

5. Two-layer Flows in Rotating Channels.

Long-term evolution of strongly nonlinear internal solitary waves in a rotating channel

Internal Wave Modeling in Oceanic General Circulation Models

A comparison between vertical motions measured by ADCP and inferred from temperature data

Exact Solitary Wave Solutions in Shallow Water

An Atlas of Oceanic Internal Solitary Waves (February 2004) by Global Ocean Associates Prepared for Office of Naval Research Code 322 PO

Sediment resuspension from shoaling internal solitary waves. Daniel Bourgault, Institut des sciences de la mer de Rimouski, Rimouski, QC, Canada

Scattering of Internal Gravity Waves at Finite Topography

Experimental Investigations of Nonlinear Internal Wave Generation and Evolution

Erratic internal waves at SIO Pier. data and wavelet analysis courtesy of E. Terrill, SIO

Combined E ect of Rotation and Topography on Shoaling. Oceanic Internal Solitary Waves. Roger Grimshaw. Chuncheng Guo. Karl Helfrich.

Energy Budget of Nonlinear Internal Waves near Dongsha

196 7 atmospheric oscillations:

Effect of the North Equatorial Counter Current on the generation and propagation of internal solitary waves off the Amazon shelf (SAR observations)

Internal boundary layers in the ocean circulation

Lecture 12: Transcritical flow over an obstacle

Understanding inertial instability on the f-plane with complete Coriolis force

Donald Slinn, Murray D. Levine

New England Shelf and the Gulf of Maine

Figure 1. adiabatically. The change in pressure experienced by the parcel is. dp = -ρ o gξ

Statistical Prediction of Ocean Circulation and Trajectories

Island Wakes in Shallow Water

Modelling the morning glory of the Gulf of Carpentaria

Impact of stratification and climatic perturbations to stratification on barotropic tides

Periodic variability of currents induced by topographically trapped waves in the coastal zone in the Gulf of Finland*

Internal inertio-gravity waves in the laboratory: Mechanisms, properties, and impacts

Numerical simulations of the interaction of internal waves with a shelf break

Analysis of the 500 mb height fields and waves: testing Rossby wave theory

Paper C Numerical studies of nonhydrostatic pressure effects on the wind forced boundary layer dynamics

What Is a Soliton? by Peter S. Lomdahl. Solitons in Biology

Exponential asymptotics theory for stripe solitons in two-dimensional periodic potentials

Tsunamis and ocean waves

2013 Annual Report for Project on Isopycnal Transport and Mixing of Tracers by Submesoscale Flows Formed at Wind-Driven Ocean Fronts

Generation of undular bores and solitary wave trains in fully nonlinear shallow water theory

Experiments on capillary-gravity waves of solitary type on deep water

Internal gravity waves and topography: emission, propagation and reflection.

British Colombia Knight Inlet Strait of Georgia Strait of Juan de Fuca

Internal Solitary Wave-Induced Bed Failure in a Stratified Water Column

Impact of Typhoons on the Western Pacific Ocean DRI: Numerical Modeling of Ocean Mixed Layer Turbulence and Entrainment at High Winds

Transcription:

Abstract Local generation of internal solitary waves in an oceanic pycnocline Nicolas Grisouard 1,2 and Chantal Staquet 1 1 Laboratoire des Ecoulements Géophysiques et Industriels, Grenoble, France 2 Courant Institute of Mathematical Sciences, New-York, USA grisouard@cims.nyu.edu Oceanic observations have provided evidence of the generation of internal solitary waves due to an internal tidal beam impinging on a seasonal pycnocline from below, a process referred to as local generation. We shall present the first direct numerical simulations of such a generation process with a fully nonlinear non-hydrostatic model (MITgcm) for an idealised twodimensional configuration. The forcing mechanism is a simple resonance so that, for a mode-n wave to be generated in the pycnocline, the horizontal phase speed of the beam has to match the phase velocity of that mode. This mode next evolves into an internal solitary wave. We shall provide examples of internal solitary waves trapped in the pycnocline as first, second and third modes. This work has been done in collaboration with Theo Gerkema (NIOZ, NL). We shall also present two-dimensional numerical simulations of situations that are relevant for the Bay of Biscay in summer. 1 Introduction Oceanic observations from the Bay of Biscay, Portugal seas, Mozambique Channel or Mascarene Ridge have provided evidence of the generation of internal solitary waves due to an internal tidal beam impinging on a seasonal pycnocline from below, a process referred to as local generation (as opposed to the more direct generation over topography). This process was proposed by New and Pingree (1990) to explain observations of internal solitary waves (ISWs) in the central Bay of Biscay. They argued that these ISWs are not directly related to topography (the continental slope), but are generated by an internal tidal beam hitting the seasonal pycnocline from below. Measurements of the beam location were consistent with ray paths computation based upon the local stratification and upon the bathymetry of the Bay of Biscay (Pingree and New, 1991). Further evidence of local generation of ISWs was presented by New and Pingree (1992). Since the beginning of the 2000s, new remote sensing techniques namely SAR (New and da Silva, 2002; da Silva et al., 2002) and ocean colour imaging (Azevedo et al., 2006) have confirmed the occurrence of this mechanism in spring and summer, when the seasonal pycnocline is present. Few attempts have so far been made to explain quantitatively the underlying physical mechanism. The effect on a pycnocline of an incident plane internal gravity wave was studied experimentally and theoretically by Delisi and Orlanski (1975) for the first time. Following the work of New & Pingree in the 90 s, the theoretical model of Delisi and Orlanski (1975) was considered again by Thorpe (1998). This work was further extended by Akylas et al. (2007), in the limit of long waves and when an internal wave beam (instead of a simple plane wave of infinite horizontal extent) impinges on the pycnocline. Akylas et al. (2007) addressed both the linear development of the interfacial waves and proposed a model equation for the far-field evolution involving weakly nonlinear, weakly nonhydrostatic and radiative effects which admits soliton-like solutions. Their result agreed with the finding of Gerkema (2001), namely that the initial phase is essentially linear, nonlinear (and nonhydrostatic) effects becoming important in 7 th Int. Symp. on Stratified Flows, Rome, Italy, August 22-26, 2011 1

the second phase, when the perturbation of the pycnocline propagates away, steepens and may break up into ISWs. Gerkema (2001) solved the linear problem analytically in terms of vertical modes. All parameters were representative of the oceanic configuration of the Bay of Biscay and were fixed, except for the density jump across the interface, which was varied. Gerkema (2001) concluded that the density jump across the pycnocline has to be moderate for the displacement amplitude of the interface to be maximum. If the density jump is stronger, the beam is reflected at the interface; if weaker, the beam is transmitted across the discontinuity as an evanescent wave. All these studies are limited in that they address at best weakly nonlinear and weakly nonhydrostatic effects and treat the pycnocline as an interface, i.e. a layer of infinitesimal thickness (with the exception of Maugé and Gerkema (2008) who used stratification profiles obtained from measurements in the Bay of Biscay). In the present paper, we present numerical results obtained with the MIT general circulation model; we are thus able to relax all these constraints at once: we allow a finite thickness of the pycnocline as well as fully nonlinear and nonhydrostatic effects. Our aim is to investigate when ISWs can be obtained in this more general situation, what the influence is of the finite thickness of the pycnocline upon the wave structure and whether the criteria derived by Delisi and Orlanski (1975), Gerkema (2001) and Akylas et al. (2007) for maximum pycnocline displacement still hold. As a second step, we consider a configuration guided by measurements in the Bay of Biscay to investigate the occurence of ISWs in a realistic context. This work has been published in two recent papers, which we summarize below. 2 Generation of internal solitary waves in a pycnocline by an internal wave beam: a numerical study In Grisouard et al. (2011), we addressed the local generation of ISWs, when an internal wave beam impinges on a pycnocline. As just mentioned, nonlinear, nonhydrostatic numerical experiments were conducted for this purpose, the vertical density profile being continuous and made up of a lower layer with constant value N 0 of the Brunt-Väisälä frequency, a pycnocline of finite thickness and a thin homogeneous upper layer. We showed that ISWs can be generated and that the finite thickness of the pycnocline allows different vertical modes to be excited, depending on the parameters. We ran numerical experiments showing that mode-1, -2 or -3 internal solitary waves can develop (Figure 1). We next proposed two different approaches to predict the conditions of occurrence of a mode-n pycnocline wave, which depend fundamentally upon the finite thickness of the pycnocline. In the far-field approach, that is when harmonics of the forcing frequency have appeared, we showed that the observed mode-n weakly nonlinear pycnocline waves have the same horizontal phase speed as the incident wave beam. We demonstrated this result heuristically by using a modal decomposition based upon the actual density profile and a frequency of value higher than N 0. For a mode-1 wave, this result extends to a pycnocline of finite thickness the linear criterion found for Akylas et al. (2007). This result allowed us to design numerical experiments so that a mode-1, -2 or -3 pycnocline wave develops, the density jump across the pycnocline increasing with the mode number (or equivalently, the horizontal wavelength of the wave beam decreasing with the mode number). One might wonder why the phase speed of a beam should be matched with the phase speed of trapped waves, that are physically some distance away from the beam. We observed in Section 3.2 that the evolution towards ISWs actually takes one or two wavelengths only, a much shorter distance than for weakly nonlinear waves described e.g. by a KdV model. This is very 7 th Int. Symp. on Stratified Flows, Rome, Italy, August 22-26, 2011 2

Figure 1: Generation of a mode-1 ISW by an internal wave beam. (a): Magnification of the isopycnals around the center line of the pycnocline, with a train of 3 deformations explicitly indicated with arrows (the leading deformation has almost passed by the rear deformation of the preceeding train). (b): Horizontal velocity field allowing to locate the internal wave beam; values range from 7.5 mm.s 1 (dark blue) to 7.5 mm.s 1 (red). likely due to the large amplitude of the forcing which allows the development of nonlinear structures in the vicinity of the beam impact, that propagate with the same speed as the forcing phase speed. The near-field approach consists of deriving simple geometrical conditions ensuring that a mode-n wave develops at the beam impact (the model holds for n 2 only). In this approach, the dynamics is linear, a simple plane wave is considered in place of a wave beam and a piecewise three-layer stratification is assumed. We showed that a mode-n wave is forced if a simple relation is satisfied which involves the thickness of the pycnocline, the angle of the refracted internal wave in the pycnocline and the horizontal wavelength of this wave. This model was originally designed as a qualitative argument but in spite of its great simplicity, it also gave a good quantitative agreement with our numerical data for n 2. It explains why a beam with a short wavelength generates high mode waves or, equivalently, why such high mode waves are favoured by a strongly stratified pycnocline. When related to the conclusion of the far-field approach, this model provides a simple analytical expression for the phase speed of a mode-n ISW. We discussed the more general situation in which the wave beam has an intermediate phase speed between those of modes n and n + 1. With the help of an example, we conjectured that both modes could be selected in this case and that the lower the mode-number is, the larger its bandwidth is. The computation of the amplitude of those modes would be required to answer precisely this problem. This requires to solve the nonlinear nonhydrostatic Boussinesq equations forced by the wave beam, a rather involved task. 7 th Int. Symp. on Stratified Flows, Rome, Italy, August 22-26, 2011 3

3 Numerical simulations of the local generation of internal solitary waves in the Bay of Biscay In Grisouard and Staquet (2010), we presented two-dimensional numerical experiments using parameters representative of the Bay of Biscay when local generation is observed. Two first experiments addressed the problem of an internal tide beam impinging on a pycnocline initially at rest. The first experiment E1, using a mid to late summer stratification that is typical of the Bay of Biscay when local generation is observed and a realistic internal tide beam did not lead to a convincing generation of ISWs, disagreeing with the observations. The second experiment E2 used a slightly reduced density jump across the seasonal pycnocline similar to that observed in late spring. In this case mode-1 ISWs were generated. The results of E1 and E2 agreed with the predictions of the far-field condition derived in GSG10. Nonetheless, the distance between two trains of ISWs in E2 displays significant discrepancies with the observations. In the Bay of Biscay, this distance is about 50 km, which is also the dominant wavelength of the internal tide in the upper 100 m (referred to as the upper tide), regardless of the characteristics of the beam. In E2 by contrast, the ISW trains were separated by 25 to 30 km, close enough to the horizontal wavelength of the forcing beam for it to generate mode-1 ISWs. A possible explanation for this discrepancy relies upon the interaction of the beam with the upper tide generated at the shelf break and propagating throughout the domain. A third numerical experiment E3 was thus designed, forced by the internal tide generated by a barotropic flow on the continental slope. The stratification was the same as in E1 and as in the Bay of Biscay when local generation occurs. This experiment showed that mode-1 ISWs are generated, each train being now separated by about 55 km, in agreement with observations. To investigate the role of the upper tide in the generation process of these ISWs, a fourth and last experiment E4 was designed, which forcing included the upper tide but not the beam. No ISWs appeared but the interfacial displacement could be superposed on the carrier wave of the ISWs observed in the third experiment. To sum up, E1 and E2 showed that an internal tide beam impinging on a seasonal pycnocline can generate ISWs but that the mechanism is very sensitive to the amplitude of the density jump across the seasonal pycnocline, in agreement with GSG10 predictions, as well as the predictions of previous works (Gerkema, 2001; Akylas et al., 2007). E2 also showed that mode-1 ISWs are generated but with quantitatively different features than those of the ISWs observed in the Bay of Biscay. On the other hand, the upper tide alone does not evolve into trains of ISWs. Only when an internal tide beam interacts with the upper tide is the formation of ISWs observed, with physical characteristics close to the observations. The propagation speeds and distance between trains of ISWs are now controlled by the upper tide but the existence of these ISWs is most likely due to an energy input from the inernal tide beam in the pycnocline. Conjectures on the mechanism of how energy is transferred to the pycnocline in the presence of the upper tide can be proposed. For instance, the upper tide can be modelled as an oscillating horizontal shear flow. As this oscillating motion is synchronous with those induced in the pycnocline by the beam, the phase shift between the beam and the upper tide may be such that the phase speed of mode-1 ISW would be reduced by Doppler effect at the location of the impact of the beam and hence brought closer to the beam phase speed. If it is brought close enough, the selection condition discussed in Grisouard and Staquet (2010) would predict that mode-1 ISWs are generated. Further analytical work is required to address the validity of this conjecture. We are quite aware that several effects have been neglected. First, our numerical experiments are two-dimensional. Also, Gerkema et al. (2008) showed that the internal tide beam 7 th Int. Symp. on Stratified Flows, Rome, Italy, August 22-26, 2011 4

Figure 2: Snapshot at 16 periods of the barotropic tide after the start of the simulation. (a): Contours of a few isopycnals around the pycnocline. A mode-2 ISW is framed. (b): Horizontal velocity field in the domain, the colour scale ranging from -12 (blue) to 12 cm/s (red). in the Bay of Biscay is very sensitive to non-traditional effects. Including these effects would modify the initialisation of experiments E3 and E4. Also, the relative amplitudes of the beam and of the upper tide might be different than those set in E3. Performing three-dimensional realistic non-hydrostatic (and nonlinear) numerical experiments of the Bay of Biscay, still a computational challenge nowadays, appears to be necessary to assess the respective roles of the internal tide beam and the upper tide in the generation of ISWs and in the control of their spatial structure. References Akylas, T. R., Grimshaw, R. H. J., Clarke, S. R., and Tabaei, A. (2007). Reflecting tidal wave beams and local generation of solitary waves in the ocean thermocline. J. Fluid Mech., 593:297 313. Azevedo, A., da Silva, J. C. B., and New, A. L. (2006). On the generation and propagation of internal solitary waves in the southern Bay of Biscay. Deep-Sea Res. Pt. I, 53:927 941. da Silva, J. C. B., New, A. L., Srokosz, M. A., and Smyth, T. J. (2002). On the observability of internal tidal waves in remotely-sensed ocean colour data. Geophys. Res. Lett., 29(12):120000 1. Delisi, D. P. and Orlanski, I. (1975). On the role of density jumps in the reflexion and breaking of internal gravity waves. J. Fluid Mech., 69:445 464. Gerkema, T. (2001). Internal and interfacial tides: Beam scattering and local generation of solitary waves. J. Mar. Res., 59:227 255. 7 th Int. Symp. on Stratified Flows, Rome, Italy, August 22-26, 2011 5

Gerkema, T., Zimmerman, J., Maas, L., and van Haren, H. (2008). Geophysical and astrophysical fluid dynamics beyond the traditional approximation. Rev. Geophys., 46:RG2004, doi:10.1029/2006rg000220. Grisouard, N. and Staquet, C. (2010). Numerical simulations of the local generation of internal solitary waves in the Bay of Biscay. Nonlin. Processes Geophys, 17:575 584. Grisouard, N., Staquet, C., and Gerkema, T. (2011). Generation of internal solitary waves in a pycnocline by an internal wave beam: a numerical study. J. Fluid Mech., 676:491 513. Maugé, R. and Gerkema, T. (2008). Generation of weakly nonlinear nonhydrostatic internal tides over large topography: a multi-modal approach. Nonlinear Proc. Geoph., 15:233 244. New, A. L. and da Silva, J. C. B. (2002). Remote-sensing evidence for the local generation of internal soliton packets in the central Bay of Biscay. Deep-Sea Res. Pt. I, 49:915 934. New, A. L. and Pingree, R. D. (1990). Large-amplitude internal soliton packets in the central Bay of Biscay. Deep-Sea Res., 37(3):513 524. New, A. L. and Pingree, R. D. (1992). Local generation of internal soliton packets in the central Bay of Biscay. Deep-Sea Res., 39(9):1521 1534. Pingree, R. D. and New, A. L. (1991). Abyssal Penetration and Bottom Reflection of Internal Tidal Energy in the Bay of Biscay. J. Phys. Oceanogr., 21:28 39. Thorpe, S. A. (1998). Nonlinear Reflection of Internal Waves at a Density Discontinuity at the Base of the Mixed Layer. J. Phys. Oceanogr., 28:1853 1860. 7 th Int. Symp. on Stratified Flows, Rome, Italy, August 22-26, 2011 6