# Generation of Hydroacoustic Waves by an Oscillating Ice Block in Arctic Zones

Authors: Usama Kadri
Journal: Advances in Acoustics and Vibration
 (2016)

---

## Abstract

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.


---

## Body

## 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 Analysis

### 2.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-x2exp⁡iωt,where t denotes time, H is the Heaviside step function, and i≡-1.Figure 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+hsin⁡knbexp⁡-iknx-ωt-π/2kncos⁡μnh,x>b,-2iζ0μncos⁡μnz+hcos⁡knbexp-iknx-ωtkncos⁡μnh,x<b,(3)wnx,z,t=∓2ζ0μn2sin⁡μnz+hsin⁡knbexp⁡-iknx-ωt-π/2kn2cos⁡μnh,x>b,-2ζ0μn2cos⁡μnz+hsin⁡knbexp-iknx-ωtkn2cos⁡μnh,x<b,where μn is the eigenvalue and kn is the horizontal wavenumber.


## 2.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-x2exp⁡iωt,where t denotes time, H is the Heaviside step function, and i≡-1.Figure 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+hsin⁡knbexp⁡-iknx-ωt-π/2kncos⁡μnh,x>b,-2iζ0μncos⁡μnz+hcos⁡knbexp-iknx-ωtkncos⁡μnh,x<b,(3)wnx,z,t=∓2ζ0μn2sin⁡μnz+hsin⁡knbexp⁡-iknx-ωt-π/2kn2cos⁡μnh,x>b,-2ζ0μn2cos⁡μnz+hsin⁡knbexp-iknx-ωtkn2cos⁡μnh,x<b,where μn is the eigenvalue and kn is the horizontal wavenumber.


## 3. Results and Discussion

### 3.1. Displacements

In order to obtain the horizontal and vertical particle displacement componentsXn;Zn we integrate (2) with respect to time, so that, for x>b, we can write Xn=Asin⁡knx-ωt and Zn=Bcos⁡knx-ω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+hsin⁡knbkncos⁡μnh;B=∓2ζ0ωμn2sin⁡μnz+hsin⁡knbkn2cos⁡μnh.The orbital behaviours of the acoustic modes are periodic with depth as shown in Figure2, 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.Figure 2
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.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μnsin⁡2μnh+2μnhsin2⁡knbhc2kn4cos2⁡μ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].Figure 4
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 blockx=±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.Figure 5
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.Figure 6
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.


## 3.1. Displacements

In order to obtain the horizontal and vertical particle displacement componentsXn;Zn we integrate (2) with respect to time, so that, for x>b, we can write Xn=Asin⁡knx-ωt and Zn=Bcos⁡knx-ω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+hsin⁡knbkncos⁡μnh;B=∓2ζ0ωμn2sin⁡μnz+hsin⁡knbkn2cos⁡μnh.The orbital behaviours of the acoustic modes are periodic with depth as shown in Figure2, 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.Figure 2
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.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μnsin⁡2μnh+2μnhsin2⁡knbhc2kn4cos2⁡μ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].Figure 4
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 blockx=±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.Figure 5
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.Figure 6
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.



---

*Source: 8076108-2016-07-28.xml*