The time harmonic problem of propagating hydroacoustic waves generated in the ocean by a vertically oscillating ice block in arctic zones is discussed. The generated acoustic modes can result in orbital displacements of fluid parcels sufficiently high that may contribute to deep ocean currents and circulation. This mechanism adds to current efforts for explaining ocean circulation from a snowball earth Neoproterozoic Era to greenhouse earth arctic conditions and raises a challenge as the extent of ice blocks shrinks towards an ice-free sea. Surprisingly, unlike the free-surface setting, here it is found that the higher acoustic modes exhibit a larger contribution.
1. Introduction
The majority of research on ocean circulation neglects the compressibility of water, which is mostly justified. However, it is well known that considering the slight compressibility of water may give rise to compression-type waves [1], known as hydroacoustic waves, or acoustic-gravity waves (under free-surface conditions); see also [2]. It has been shown recently that these waves may contribute to deep ocean currents and circulation [3], as they are continuously generated in the ocean by wind-wave interactions [4, 5], interaction of nearly opposing waves [1, 6–10], and submarine earthquakes [2].
To our knowledge, the vast majority of literature related to this type of waves considers a free-surface boundary condition. As a result, the first mode (often referred to as the zeroth mode) is a surface gravity mode, and its evanescent compartments turn into progressive acoustic modes above a certain cut-off frequency [2]. However, in the case of a rigid boundary replacing the free-surface, propagating and trapped acoustic modes might behave differently as elegantly presented by [11]. In this paper we limit our study to the generation of progressive acoustic modes, by an oscillating block of horizontal ice sheet, and their propagation is confined by horizontal rigid boundaries: the sea floor at the bottom and ice sheets at the top. Obviously, under these conditions and neglecting the elasticity of the boundaries, there can be no gravity surface mode, and all energy is transferred to the acoustic modes. Thus, studying such conditions sheds light on wave energy transfer, as well as water transportation and circulation in confined seas. Here, we present results relevant to arctic zones, from a Neoproterozoic Era [12] to current arctic zones conditions. Oscillations of large ice blocks may be triggered by atmospheric and ocean currents [13], localised wind storms, ice quakes, or other violent geophysical processes.
2. Theoretical Analysis2.1. Assumptions and Governing Equation
The system considered here is two-dimensional Cartesian coordinates (x,z) with the origin in the midpoint of the ice block; the x-axis is horizontally parallel to the ice sheet (and sea floor), and the z-axis is vertically upwards as shown in Figure 1. The ocean has a constant depth h. The ice block, assumed infinitely long with a total width of 2b, is oscillating vertically at the angular frequency ω and a small amplitude ζ0 such that the displacement of the ice block is given by(1)ζx,t=ζ0Hb2-x2expiωt,where t denotes time, H is the Heaviside step function, and i≡-1.
Schematic representation of the flow domain.
The governing equation is the standard two-dimensional compressible wave equation [14, 15]. Applying the appropriate boundary conditions and using standard techniques the following horizontal and vertical fluid velocities, un and wn, induced by the nth acoustic mode can then be derived (see the Appendix):(2)unx,z,t=∓2iζ0μncosμnz+hsinknbexp-iknx-ωt-π/2kncosμnh,x>b,-2iζ0μncosμnz+hcosknbexp-iknx-ωtkncosμnh,x<b,(3)wnx,z,t=∓2ζ0μn2sinμnz+hsinknbexp-iknx-ωt-π/2kn2cosμnh,x>b,-2ζ0μn2cosμnz+hsinknbexp-iknx-ωtkn2cosμnh,x<b,where μn is the eigenvalue and kn is the horizontal wavenumber.
3. Results and Discussion3.1. Displacements
In order to obtain the horizontal and vertical particle displacement components Xn;Zn we integrate (2) with respect to time, so that, for x>b, we can write Xn=Asinknx-ωt and Zn=Bcosknx-ωt, where Xn/A2+Zn/B2=1 is the equation of an ellipse with a horizontal semiaxis A and a vertical semiaxis B given by(4)A=±2ζ0ωμncosμnz+hsinknbkncosμnh;B=∓2ζ0ωμn2sinμnz+hsinknbkn2cosμnh.
The orbital behaviours of the acoustic modes are periodic with depth as shown in Figure 2, with n/2 wavelengths fitting the depth h, in opposition to the free-surface setting in which 2n-1/4 wavelengths fit the depth [3, 16]. For the numerical calculations we considered a speed of sound c=1500m/s, h=4500m, b=104m, ζ0=1m, and ω=5 rad/s. With the current setting, the obtained displacements correspond to horizontal and vertical velocities reaching up to tens of centimetres per second for the first mode, which is comparable to those obtained by an equivalent underwater earthquake with a free surface [2, 3]. However, a distinguishable result here is that the leading mode is not necessarily the first acoustic, but the highest acoustic mode (which can be obtained from (A.9)) given by(5)Nmax=ωhπc,where the special brackets represent the floor function. For the case of n=Nmax, the displacements and velocities may become one order of magnitude larger than those of the first acoustic mode, n=1, and thus, unlike the free-surface problem, all progressive acoustic modes have to be considered.
Semiaxis of water parcels for acoustic modes, from top to bottom: n=1,2,3,4; h=4500m, c=1500m/s, b=104m, ζ0=1m, and ω=5rad/s.
Note that the choice of ζ0=1m far exceeds observed amplitudes of oscillating sea-ice blocks in the late Common Era [17]. Nevertheless this choice allows a proper comparison with the free-surface underwater earthquake problem [2] emphasising the importance of the higher modes and directly reflects on two actual scenarios. The first scenario describes an underwater earthquake in the arctic ocean, in particular in winter if the sea surface is completely frozen. Recall that in the absence of gravity and elasticity, while the mathematical solution in hand is valid for a rigid sea floor and an oscillating ice block at the surface (oscillating ice block problem), it is equally valid for the setting with a rigid frozen surface and an oscillating sea floor block at the bottom (e.g., earthquake in the arctic ocean). The second scenario considers the snowball earth Neoproterozoic Era, where sea-ice thickness might have exceeded a kilometre. Under such conditions, ice quakes and other violent geophysical phenomena could have caused oscillations with amplitudes that are comparable with the chosen ζ0=1m. Further qualitative analysis through the different climate eras is given in the last section.
3.2. Dispersion Relation
The dispersion relation is given in nondimensional form by(6)μnknh2sinhμnh=0.For the incompressible case ωh/c→0, and therefore μn→kn. In this case there is no nontrivial solution, and no waves are formed. For the compressible case, the wavenumber of the acoustic mode is given by(7)knh=ω2h2c2-n2π21/2,n=1,…,Nmax,where Nmax is the maximum number of progressive acoustic modes. Thus, the cut-off frequency for each mode is given by ωc=nπc/h, so that, at a given depth and relatively low frequencies ω<ωc, only evanescent modes can exist, whereas when ω>ωc there should be at least one nonevanescent acoustic mode. For the case of h=4500m the cut-off frequency is ωc=π/3. The dispersion relations of knh versus ωh/c are shown in Figure 3. It is easy to show that the group velocities, defined as cg=dω/dkn, differ from one mode to another though all tend to the sound speed c if the normalised frequency ωh/c is large enough. These can be obtained from the slopes of the dispersion relations presented in Figure 3.
Dispersion of modes; h=4500m and c=1500m/s.
3.3. Energy Aspects
The total wave energy is composed of two parts: kinetic and elastic potential; with the absence of a free surface there is no gravity potential energy. Following similar steps as given by [2] we obtain an expression for the kinetic energy per unit volume:(8)Eκ,n=12ρζ02ω4μnsin2μnh+2μnhsin2knbhc2kn4cos2μnhand the elastic energy is Eϵ,n=Eκ,n/4. Above the cut-off frequency the oscillation of the ice block results in energy “transfer” to the corresponding acoustic mode/modes (Figure 4). However, below the cut-off frequency, the fluid becomes incompressible and no energy transfer occurs. With regard to this point it is worth noting that below the cut-off frequency the stratification of the ocean or the elasticity of the sea floor and ice sheets, which was neglected in the presented analysis, might become important [18]. Considering stratification may give rise to internal gravity waves that are believed to be a major driver of deep ocean mixing [19]. On the other hand, considering the elasticity of the boundaries, it is anticipated that oscillations below the cut-off frequency would give rise to a Scholte-type modes that propagate through both interfaces, with the sea floor and with the ice sheet [20]. Other types of trapped modes are discussed in detail by [11].
Kinetic energy versus frequency; h=4500m, c=1500m/s, b=104m, and ζ0=1m.
It is also notable that the fact that the higher acoustic modes have a larger contribution on water transport compared to the lower modes suggests that the flow field would become temporally and spatially nonuniform, if stratification is considered, as also found in observations [21–25]. Nevertheless, stratification or consequently variations in the speed of sound do not practically have a direct effect on the fluid motion induced by a hydroacoustic wave, whose group velocity remains almost unaltered. Hence, stratification is not considered in the current analysis, though a numeric comparison with the effects of internal gravity waves due to stratification is briefly presented in the following section.
3.4. From Snowball to Greenhouse Earth
The water motion in the ocean plays a prominent role in distributing heat, originated from solar radiation or hydrothermal activities, around the globe. In free-surface open seas, ocean water results in humidification and increasing the temperature of the surrounding air, which in turn forms storms and rain. Without currents and underwater streams localised temperatures become more severe leading to extreme global weather. Hence, water motion in the ocean acts as a global climate regulator.
The flow generated by the motion of the ice block is equivalent to the superposed flow by two line sources located at the edges of the ice block x=±b. Such superposition can be either destructive or constructive, resulting in various peak velocities, as given qualitatively in Figure 5. A minimal flow field is created when the oscillating ice block has a width 2b multiples of a wavelength for each acoustic mode; that is, 2b=NLn, where Ln=2π/kn is the wavelength, and N=1,2,3,…. On the other hand, a maximal flow field is induced when 2b=Ln(N-1/2). Thus, the mutual existence of multiple acoustic modes is expected to create regions of minimal and maximal peak velocities as confirmed by Figure 6. It is worth mentioning here that deep ocean currents induced by internal gravity waves can travel at speeds varying from 0.02 to 0.1 m/s, whereas surface currents (that do not exist in snowball earth) can travel at speeds as high as 2.5 m/s. Thus, from Figure 6, it is clear that internal gravity waves become more important in a greenhouse earth, as opposed to hydroacoustic waves, keeping in mind that the latter, despite the lower induced peak velocities, would still span the entire depth.
Peak velocities induced by four different oscillating ice blocks, b=0.1,1,10, and 100 km, with various amplitudes 0<ζ0<1; h=4500m and c=1500m/s.
Peak velocities from various ice block widths representing timeline from snowball to greenhouse earth, from right to left. The amplitude ζ0=0.0001b; h=4500m and c=1500m/s.
It is believed that during the Neoproterozoic Era the earth was fully covered by snow, with ice blocks that extended over hundreds of metres thickness [12, 26]. It was suggested by [26] that during the Neoproterozoic Era ice could have extended to hundreds or thousands of kilometres, with sea ice that could have covered up to 14×106km2 of the Arctic Ocean during winter [27]. Thus, ice blocks could have been of the scale of tectonic plates [26], with similar vertical displacements, that probably exceeded our choice of ζ0=1m, in particular with violent geophysical events that characterised that period. In this regard, it is also expected that multiple extremely long period oscillations could have formed high amplitude progressive waves causing circulations at velocity magnitudes of centimetres per second that are many orders of magnitudes higher than those from the suggested geothermal mechanism (i.e., micrometers per second) by [12] for that era. Here, it is also notable that, due to the differences in scales and distribution of events, the significance of two mechanisms could be decoupled in the space-time domain.
In the eighties of the last century the mean extent of the ice covering the arctic ocean had shrank to a value of less than 16 million squared kilometres, decreasing at a rate of 3% per decade [13]. Clearly, the ice pack thickness was dramatically reduced, and as the sea ice thins and shrinks further, it is expected that ice block oscillations will take place at relatively smaller areas and have less impact. In addition, recent observations suggest that, under current conditions, ice block oscillation amplitude may not exceed a few millimetres [17], producing a weak pressure signature on fluid parcels residing through the water column. Due to the different nature of waves radiated from ice-sheet movements, which is relatively a new scientific discipline, or even near sea glaciers, the available high quality data goes back to a few years only. It is well known that hundreds of ice quakes occur annually in localised regions, some of which are related to glacial-sea interface and have magnitudes that may exceed 5 [28]. At such discontinuous space-time interaction events, any impact is expected to be localised, and the contribution to water circulation becomes far less efficient, compared to snowball earth conditions. A qualitative demonstration of the change of impact of the ice block oscillations from snowball to greenhouse earth is given in Figure 6 where the x-axis represents the timeline, from right to left.
The Arctic Ocean is relatively poor in marine life, in general, and in plants, in particular. Although Phytoplankton is amongst the few marine plants in the Arctic Ocean, they are essential for sustaining a healthy ocean. Phytoplankton feed on nutrients from currents and water transport. As the greenhouse effect sustains and more areas of ice packs shrink and thin, one can expect a decrease in the generation of acoustic, by the present mechanism. On the other hand, if the sea becomes completely ice free [29], that gives rise to local formation of acoustic by interacting surface gravity waves [7, 9, 10]. The transition between the different mechanisms would possibly disturb global current patterns. Such disturbance is expected to affect nutrients transport and Phytoplankton population and potentially cause severe climate changes [30, 31], which is left for future work.
AppendixDerivation of the Displacements
The governing equation is the standard two-dimensional compressible wave equation [14, 15]:(A.1)φtt-c2∇2φ=0,-h≤z≤0,where φ is the velocity potential.
The upper boundary condition is given by(A.2)φz=ζt,z=0,and the bottom boundary condition is simply(A.3)φz=0,z=-h.Since here we only consider time harmonic problems we can write(A.4)φx,z,t=fx,zexpiωt.Substituting (A.4) into (A.1), (A.2), and (A.3) and making use of (1) give(A.5)∇2f+ω2c2f=0,-h≤z≤0,fz=ζ0Hb2-x2,z=0,fz=0,z=-h.
In the space domain the upper boundary condition is inhomogeneous. To overcome this difficulty we apply a Fourier transformation to the wavenumber domain, where F(w,z) is the transform of f(x,z), defined by(A.6)Fw,z=∫-∞∞exp-iwxfx,zdxand the inverse transform, defined by(A.7)fx,z=12π∫-∞∞expiwxFw,zdx.Defining w=k+iλ, where k and λ are real positive numbers, the eigenvalues can be written as(A.8)w=0,±ks,±kn,±iλn,n=1,2,…,where(A.9)ks=ωc;kn=ks2-μn21/2,μn<ks;λn=μn2-ks21/2,μn>ks.
Using standard techniques as given by [2], in which a Fourier transform is applied over the horizontal axis, we obtain an expression for f(x,z). For the range |x|>b we write(A.10)fx,z=±2ζ0∑n=1Nμncosμnz+hsinknbexp-iknx-π/2kn2cosμnh∓∑n=N+1∞μncosμnz+hsinhλnbexp∓iλnxλn2cosμnh,where the upper signs of ± and ∓ in (A.10) are for x>b and the lower ones are for x<b. On the other hand, for the range x<b, we write(A.11)fx,z=2ζ0∑n=1Nμncosμnz+hcosknbexp-iknxkn2cosμnh-∑n=N+1∞μncosμnz+hcoshλnxexp-λnbλn2cosμnh+cosksz+hkssinksh.
The first terms of right-hand side of (A.10) and (A.11) represent the acoustic modes, whereas the second terms represent the evanescent acoustic modes, which exponentially decay from the source and thus will not be considered in the sequel. The third term of (A.11) represents the surge mode which exists only below the oscillating ice block at x<b.
Ignoring the evanescent and surge terms in (A.10) and (A.11), substituting into the potential equation (A.4), and differentiating with respect to x and z give the horizontal and vertical velocities, respectively.
Competing Interests
The author declares that there are no competing interests.
Acknowledgments
The author thanks Carl Wunsch for raising the question that resulted in this paper and Raffaele Ferrari for useful discussions.
Longuet-HigginsM. S.A theory of the origin of microseisms195024313510.1098/rsta.1950.0012MR0040887YamamotoT.Gravity waves and acoustic waves generated by submarine earthquakes198212758210.1016/0261-7277(82)90016-x2-s2.0-49049131479KadriU.Deep ocean water transport by acoustic-gravity waves2014119117925793010.1002/2014jc0102342-s2.0-84915758565ArdhuinF.LavanantT.ObrebskiM.MariéL.RoyerJ.-Y.d'EuJ.-F.HoweB. M.LukasR.AucanJ.A numerical model for ocean ultra-low frequency noise: wave-generated acoustic-gravity and Rayleigh modes201313443242325910.1121/1.48188402-s2.0-84885415476KedarS.Longuet-HigginsM.WebbF.GrahamN.ClaytonR.JonesC.The origin of deep ocean microseisms in the North Atlantic Ocean2008464209177779310.1098/rspa.2007.0277ZBL1185.860222-s2.0-42949104473FarrellW.MunkW.Booms and busts in the deep2010409KadriU.Wave motion in a heavy compressible fluid: revisited201549505710.1016/j.euromechflu.2014.07.008MR3281921KadriU.Triad resonance between a surface-gravity wave and two high frequency hydro-acoustic waves201655115716110.1016/j.euromechflu.2015.09.008MR34281992-s2.0-84948457039KadriU.AkylasT. R.On resonant triad interactions of acoustic—gravity waves2016788, article R11210.1017/jfm.2015.721MR34820392-s2.0-84951311006KadriU.StiassnieM.Generation of an acoustic-gravity wave by two gravity waves, and their subsequent mutual interaction2013735, article R6910.1017/jfm.2013.539MR31643812-s2.0-84901828422LintonC. M.McIverP.Embedded trapped modes in water waves and acoustics2007451-2162910.1016/j.wavemoti.2007.04.009ZBL1231.760462-s2.0-35648955038AshkenazyY.GildorH.LoschM.TzipermanE.Ocean circulation under globally glaciated snowball earth conditions: steady-state solutions2014441244310.1175/jpo-d-13-086.12-s2.0-84892649271ParkinsonC. L.CavalieriD. J.GloersenP.ZwallyH. J.ComisoJ. C.Arctic sea ice extents, areas, and trends, 1978–199619991049208372085610.1029/1999jc9000822-s2.0-0033567611StoneleyR.The effect of the ocean on Rayleigh waves1926134935610.1111/j.1365-246x.1926.tb05382.xLambH.19326thCambridge, UKCambridge University PressMR1317348JensenF. B.KupermanW. A.PorterM. B.SchmidtH.2011New York, NY, USASpringer10.1007/978-1-4419-8678-8MacAyealD. R.OkalE. A.AsterR. C.BassisJ. N.Seismic and hydroacoustic tremor generated by colliding Icebergs2008113310.1029/2008jf0010052-s2.0-84940590617EyovE.KlarA.KadriU.StiassnieM.Progressive waves in a compressible-ocean with an elastic bottom201350592993910.1016/j.wavemoti.2013.03.003MR30709502-s2.0-84878718430FerrariR.WunschC.Ocean circulation kinetic energy: reservoirs, sources, and sinks200941125328210.1146/annurev.fluid.40.111406.1021392-s2.0-65649100474BromirskiP. D.StephenR. A.Response of the ross ice shelf, antarctica, to ocean gravity-wave forcing2012536016317210.3189/2012AoG60A0582-s2.0-84870516702D'AsaroE. A.MorisonJ. H.Internal waves and mixing in the Arctic Ocean1992392S459S48410.1016/s0198-0149(06)80016-62-s2.0-0027036985HalleC.PinkelR.Internal wave variability in the Beaufort Sea during the winter of 1993/199420031087, article 321010.1029/2000JC0007032-s2.0-0141975330LevineM. D.PaulsonC. A.MorisonJ. H.Internal waves in the arctic ocean: comparison with lower-latitude observations198515680080910.1175/1520-0485(1985)015<0800:iwitao>2.0.co;2PinkelR.Near-inertial wave propagation in the western Arctic200535564566510.1175/JPO2715.12-s2.0-21844480097PlueddemannA. J.Internal wave observations from the Arctic environmental drifting buoy1992978126191263810.1029/92JC010982-s2.0-0027098720HoffmanP. F.KaufmanA. J.HalversonG. P.SchragD. P.A neoproterozoic snowball earth199828153811342134610.1126/science.281.5381.13422-s2.0-0032575722WeissJ.Scaling of fracture and faulting of ice on earth200324218522710.1023/A:10232931173092-s2.0-0038135039JoughinI.Greenland rumbles louder as glaciers accelerate200631157681719172010.1126/science.11244962-s2.0-33645228769OverlandJ. E.WangM.When will the summer Arctic be nearly sea ice free?201340102097210110.1002/grl.503162-s2.0-84879924442SewallJ. O.SloanL. C.Disappearing arctic sea ice reduces available water in the american west2004316L0620910.1029/2003GL019133FrancisJ. A.VavrusS. J.Evidence linking arctic amplification to extreme weather in mid-latitudes2012396L0680110.1029/2012GL051000