Next Article in Journal
Rapid Recruitment of Symbiotic Algae into Developing Scleractinian Coral Tissues
Next Article in Special Issue
Combining Numerical Simulations and Normalized Scalar Product Strategy: A New Tool for Predicting Beach Inundation
Previous Article in Journal
Experimental Study on the Distribution of Wave-Induced Excess Pore Pressure in a Sandy Seabed around a Mat Foundation
Previous Article in Special Issue
Downscaling Future Longshore Sediment Transport in South Eastern Australia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of Free Long Wave Generation on the Shoaling of Forced Infragravity Waves

1
Hidromod ltd., Porto Salvo, Lisbon, 2740-278 Lisbon, Portugal
2
Centre for Marine and Environmental Research, University of Algarve, 8005-139 Faro, Portugal
3
School of Civil Engineering, University of Queensland, St Lucia QLD 4072, Australia
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2019, 7(9), 305; https://doi.org/10.3390/jmse7090305
Submission received: 9 August 2019 / Revised: 28 August 2019 / Accepted: 30 August 2019 / Published: 4 September 2019
(This article belongs to the Special Issue Dynamics of the Coastal Zone)

Abstract

:
Different conceptual models for forced infragravity (long) waves exist in the literature, which suggest different models for the behavior of shoaling forced waves and the possible radiation of free long waves in that process. These are discussed in terms of existing literature. A simple numerical model is built to evaluate the wave shape in space and time during shoaling of forced waves with concurrent radiation of free long waves to ensure mass continuity. The same qualitative results were found when performing simulations with the COULWAVE model using the radiation stress term in the momentum equation to force the generation and propagation of bound waves. Both model results indicate a strong frequency dependence in the shoaling rate and on the lag of the total long wave with respect to the forcing, consistent with observations in the literature and more complex evolution models. In this approach, a lag of the long wave is only observed in the time domain, not in the space domain. In addition the COULWAVE is used to investigate dissipation rates of incident free and forced long waves inside the surf zone. The results also show a strong frequency dependence, as previously suggested in the literature, which can contribute to the total rate of decay of the incident forced wave after short wave breaking.

1. Introduction

The importance of forced long waves, or bound waves, to the oceans dynamics is well established in the literature, especially in shallow waters. However, in nature, some processes related to their generation, propagation, shoaling and dissipation are yet not fully understood, with eventual conflicting concepts to explain the same processes. A plausible partial explanation is that those concepts emerge from different interpretations of the bound wave theory first introduced by Longuet-Higgins and Stewart [1]. For steady groups, two solutions have been proposed for forced waves by wave groups, which differ in the definition of the water level about which the forced waves oscillate. These indicate a forced wave oscillating either around or solely below the mean water level for the short waves, respectively (Figure 1). The difference between the two interpretations is the constant of integration, which can be taken as the still water level or mean water level under the group, which for a steady condition may be arbitrarily chosen giving the long wave response as [1]
η ¯ = S x x ρ ( g h c g 2 ) + c o n s t a n t
and which has been widely adopted to estimate the forced wave amplitude from the magnitude of the radiation stress ( S x x ) derived from the short wave amplitude or wave envelope. Here, ρ is the constant water density, c g is the group velocity, g is the acceleration due to gravity and h is the local water depth.
This solution emerges from the mass and momentum balance equations including the radiation stress forcing function as
η ¯ t + ( h u ¯ ) x = 0
u ¯ t + g η ¯ x = 1 ρ h S x x x .
where η ¯ and u ¯ are the long wave surface displacement and horizontal velocity, respectively. The solution domain is from the bed to the free surface, and initial conditions are discussed further below.
Equation (1) is the forced part of solution to the non-homogenous differential equations, which allows free waves (the homogenous solution) to be added to the forced wave term [3]. The additional free waves may become important in a transient solution, but these have most often been ignored in analysis and modelling of bound wave propagation by treating this as a quasi-steady process. Transient solutions to similar problems dates back to Proudman [4], who proposed dynamical theories of water motions induced by a moving pressure field. Whitham [5] first addressed this issue with respect to the bound wave problem. Later, Molin [6] and Mei and Benmoussa [7] extended the investigation to sloping bottoms. The difference in interpretation of the forced wave as illustrated in Figure 1 is important during shoaling. With a bound wave with both positive and negative components then shoaling (an increase in amplitude in shallower water) does not necessarily require radiation of free waves to ensure continuity of volume, but with a bound wave of pure depression, then an increase in amplitude of that depression requires radiation of free waves to ensure continuity. Radiation of those free waves then changes the shape of the total infragravity wave during shoaling, and may also influence the commonly observed lag between incident bound waves and the incoming short wave envelope.
Nielsen and Baldock [8] presented simple analytical transient solutions to the response of the surface elevation to both non-resonant (finite water depth conditions) and resonant (shallow water conditions) radiation stress forcing. This solution has the same initial conditions as those proposed for startup of a submarine landslide that may trigger tsunami [9,10], i.e., a flat water surface and no flow velocity. In summary, for a 1-D scenario and assuming an abrupt onset of a non-resonant positive forcing with a Gaussian shape, three waves are generated: a purely negative forced wave, propagating phase locked with the forcing, and two free waves with positive amplitude and opposite direction of propagation. While the shape of all three waves are the same as the forcing, the free wave amplitudes depend on the ratio between the wave group velocity and the shallow water wave celerity as
( A f r e e + , A f r e e ) = ( A f o r c e d 2 [ 1 + c g g h ] , A f o r c e d 2 [ 1 c g g h ] )
where A f o r c e d ,   A f r e e + and A f r e e are the amplitudes of the forced and the two free waves respectively.
At the resonant condition (short waves in shallow water), the free long waves arising from the initial condition cannot radiate away, forming a N-shaped wave, with a positive hump leading the forced wave. The total long wave surface grows linearly in time, but with the shape of the derivative of the forcing, rather than the inverse shape as in the steady solution. Changes in the shape of the bound wave from a symmetrical depression to an asymmetrical depression plus surge leading wave has been reported to occur during the shoaling processes, where again purely forced [11] and a mix of forced and free waves [8] interpretations have been provided.
During short wave breaking, the assumption that bound waves are released as free waves is commonly adopted [12,13,14]. However, strong infragravity wave dissipation has also been observed, suggesting that bound waves are not released during the breaking process, but remain phase locked (forced) in the surf zone decaying with the short wave forcing [15,16,17]. Breaking of shorter infragravity waves close to the shoreline [18], friction [19] and energy transfer [20] and dissipation due to the existence of turbulence in the surf zone [13] have been proposed or suggested as additional or alternate reasons for dissipation of forced wave energy.
Even though the generation of free waves due to the unsteady nature of forced bound waves have been discussed in different studies [7,8,21], it is not clearly addressed in others, where the interpretation of the results relies on the quasi-steady solution [11,12,13,22]. To date, generation and radiation of free waves during bound wave shoaling has not been documented in observations, but the shape of the long wave and its phase (lag) relative to the forcing is observed to change. In part this may be due to errors in wave generation, the predominance of shallow water conditions, and the likely small magnitude of the free waves. Nevertheless, it is important to comprehend the possible effects of such free wave generation in the interpretation of propagating and shoaling infragravity waves. The analysis is based on numerical results and previous published laboratory data.
The paper is organised as follows: In Section 2, aspects related to the shape change of the bound wave during shoaling, the possible consequences of the generation of free waves on the observed shoaling rate and lag of the long wave are discussed. A simple semi-analytical model is constructed to evaluate the model concepts of Nagase and Mizuguchi [21] and Nielsen and Baldock [8] and to quantify how the group frequency influences the apparent shoaling rate and lag of the long waves. In Section 3, the possibility of dissipation of long waves by surf zone turbulence is investigated using the COULWAVE Boussinesq wave model, and particularly if the frequency dependence suggested by Battjes [13] occurs in the model results. Final remarks are presented in Section 4.

2. Bound Wave Shoaling

During bound wave shoaling two processes have often been observed: changes in shape, with the development of a leading positive surge and a lag of the long wave in time relative to the short wave envelope; and a shoaling rate dependence on the slope and frequency of the bound wave [12,13,18,23,24]. From the static quasi equilibrium approach, the positive surge ahead of the short wave group can be interpreted as the stronger response of surface elevation to the radiation stress gradients in shallow water and an asymmetry in the short wave envelope as the wave group shoals. The main argument presented by Baldock [11] when analysing single transient wave groups was the very close match of the spatial gradients of the wave envelope (radiation stress) with the spatial gradients of the long, as in Longuet-Higgins and Stewart [1]. A similar observation from measured data is presented in Figure 2, but for an individual wave group from a random wave case from Baldock and Huntley [25]. The lag of the long wave relative to the short wave envelope is also very small in this case, but appears larger in the time domain than in the space domain. It is important to note that these data were obtained with 2nd order generation, which minimises the generation of free (error) waves at the wavemaker (see Orzaghova et al., [26] for other examples).
The relationship between the wave envelope and the bound wave is commonly analysed based on time series from a point measurement. This is despite the fact that because the propagating wave group is not steady and not symmetric during shoaling the wave shape measured in the time domain is not the same as the instantaneous shape observed in the space domain, see direct comparisons by Baldock [11] and Figure 2 below. In the left panel in Figure 2 the long wave is in antiphase with the forcing and both are asymmetric. In the right panel, no larger leading wave is observed, even though this is present in the left panel, which means that wave only is generated or amplified very close to the shore, in shallow water. With measurements in the time domain, lags between the forcing and bound wave are often reported [22,27,28]. That occurs even when no lag is observed in the spatial domain (Figure 2).
In a time series analysis, the bound wave lag is often associated with the intensity of energy transfer from the forcing to the long wave, where a greater lag means more energy transfer. For the same beach slope, higher frequency bound waves appear to lag more than lower frequency bound waves, resulting in higher shoaling rates at higher frequencies. This difference in shoaling rates is clearly evident in data from Baldock et al. [24] and Janssen et al. [12]. Battjes et al. [13] showed that this shoaling rate is determined by the wave group length relative to the bottom slope, a normalised beach slope, varying from approximately ~ h 1 4 (free wave shoaling rate) for relatively lower frequencies to ~ h 5 2 (steady solution shoaling rate) for relatively higher frequencies. It is noted that the shoaling at higher frequencies corresponds to the amplitude growth proposed for the quasi-steady solution of Longuet-Higgins and Stewart [1,29], despite no energy transfer or lag being proposed at that time.
In the literature, two terms in the energy balance equation are relevant for modelling the energy exchanges between short and long waves. To the authors’ understanding there are still some difficulties in the interpretation of their role in the energy transfer from short to bound waves. The problem has been introduced by Longuet-Higgins and Stewart [1,29], followed by Whitham [5], who attempted to clarify the issue addressed by the previous authors. In the full energy balance introduced by Philips [30] and Schaffer [3], and using their notation for clarity, the conservation of total energy is given by
E t + W x = 0
where E is the total energy density and W is the total energy flux. The short wave energy equation is
t { E s 1 2 ρ ( h + η ) U s 2 } + { U c E s + W s 1 2 ρ ( h + η ) U U s 2 } + S x x U c x U s S x x x = 0
where subscript s and c denote the short and long waves (current) respectively. The equation for the free slow varying current or long wave is
E c t + W c x + U S x x x = 0
The two terms related to energy transfer between short wave (via S x x ) and long waves, with U as the bound long wave velocity (equivalent to u ¯ in (2) and (3) are
S xx x   U S x x x ,
S x x U x
While Battjes et al. [13] have considered ( 6 ) as the main driver of energy transfer, Henderson et al. [31], Guedes et al. [32] and others considered (7) as the principal mechanism. Schaffer [3] (see page 564), taking in consideration free and forced long waves, suggested (6) as being the interaction between S x x and the forced wave current, and (7) as the interaction between S x x and the free current (free long waves). When considering the short wave and free current (tide) interaction (7) was the only term considered by Longuet-Higgins and Stewart [29], with the forced bound wave being a product of that interaction.
Janssen et al. [12], following the work of Bowers [33] and Van Leeuwen [34], presented analytical and numerical solutions for the bound wave amplitude and phase shift induced by changes in depth. Battjes et al. [13], assuming quasi-steady waves, derived a phase-average rate of energy transfer based on the phase lag between the wave envelope and the forced wave ( Δ ϕ - the deviation from the equilibrium solution, π), the amplitudes of radiation stress ( S ^ ) and the forced long wave velocity ( U ^ )
R U S x x x 1 2 κ U ^ ( f ) S ^ ( f ) sin Δ ϕ
where κ = 2 π f / c g is the forced wave number, at individual frequencies, f , and (ˆ) denotes real amplitudes. Good agreement between this model and laboratory data was found.
Generally, the discussion regarding bound wave shoaling requires the assumption of steady wave conditions, even though the possibility of free wave generation during the shoaling process is recognized in some studies. With a dynamic process, Nagase and Mizuguchi [21] suggested that the observed smaller growth rate compared to the steady solution of Longuet-Higgins and Stewart [1] of the bound wave is a consequence of the superposition of forced and free waves generated due to the transient behaviour of the bound (forced) wave on the slope, which is similar to the suggestion of Nielsen and Baldock [8]. Based on this assumption, even if the shoaling of the forced wave part is independent of the wave group frequency and follows the equilibrium solution, due to the superposition of forced and free waves, an apparent distinct rate of energy transfer (as defined by the rate of change of the total wave amplitude) would exist for wave groups with different length and the same group velocity, c g . The apparent lag of the negative pulse would also depend on the group frequency.
Using the conceptual models of Nagase and Mizuguchi and Nielsen and Baldock, it is straightforward to calculate the evolution of the total long wave expected from superposition of propagating forced wave and free long waves, as illustrated in Figure 3. Consider a Gaussian-shaped forced wave with amplitude A and length (Lscale) with its center location at xi
η = A e ( x i x L s c a l e ) ,
propagating with constant form over a slope, and with an increasing amplitude proportional to x ~ h 5 / 2 (a shoaling rate proportional to the equilibrium solution), where x is the horizontal distance. At each new time step, the gain in amplitude by the bound wave (an increasing negative depression) is added to the forward free wave as a positive elevation with same shape as (9), which is shoaling also as a free wave ( x ~ h 1 / 4 ), to ensure continuity (1). The Gaussian depression is traveling at slower velocity, at the wave group velocity cg, which is defined in terms of the short wave period T and wave number k,
c g = g T 4 π tanh k h ( 1 + 2 k h sinh 2 k h ) ,
while the positive free wave travels faster, following the shallow water wave velocity ( g h ). Therefore, the two sets of waves slowly separate. This conceptual model is designed by adopting the definition of the forced wave proposed by Mei [2] (right panel of Figure 1). The definition by Longuet-Higgins and Stewart [1] would not necessarily require a positive free wave to balance the increasing magnitude of the forced wave depression. The results show that for groups with the same mean primary wave frequency (i.e., the same group velocity), the gain (change in the crest to trough height) and phase (lag of the long wave trough relative to the forcing) of the total infragravity wave varies with the group frequency.
Thus, even though the forced bound wave shoaling is the same as the equilibrium solution, an apparent frequency dependence appears in the shoaling rate of the total long wave. This is illustrated in Figure 4, which shows how frequency influences the total infragravity wave amplitude at a fixed position on the virtual slope. The results are qualitatively similar to the numerical results in Madsen et al. [23], their Figure 20, obtained from a Boussinesq wave model with improved dispersion in deeper water.
In this model approach, the only factor affecting the apparent changes in the shoaling rate of the total infragravity wave is the behaviour of the free wave relative to the wave group (forced wave, Equation (9). Over the same distance, the free waves generated by shorter wave groups will tend to propagate further apart from the bound wave than the free waves generated by longer wave groups, since the difference in phase speed between the forced wave and the free long wave is greater. In addition, for a longer group, the free wave needs to travel a longer distance to get ahead of the bound wave. Note that due to the difference between the free wave speed g h and group velocity, c g , a lag is expected for the total long wave (Figure 3a) and the magnitude varies according to the length of the forcing. Interestingly, the lags are also only observed in the time domain, not in the space domain. This is consistent with the observations of Baldock [11] and also with the data in Figure 2.
Based on these lags, hypothetical rates of energy transfer are calculated using Equation (8). As illustrated in Figure 5, the apparent rate of energy transfer also seems to increase with increasing frequency of the forced wave. This is the same behaviour as that predicted by the energy transfer model of Janssen et al. [12]. The apparent shoaling rate for the long wave groups are closer to ~ h 1 / 4 , i.e., similar to free long waves, while short groups are closer to ~ h 5 / 2 , qualitatively matching the different shoaling behavior in the mild-slope and the steep-slope regimes described by Battjes et al. [13]. The results presented for this conceptual model are quite simple, but they indicate that adding free waves into the bound wave shoaling process may change the results and therefore the interpretation of energy transfers, and phase lags of the forced wave.
In order to verify this hypothesis a set of simulations using the COULWAVE model [35] which adopts the fully nonlinear Boussinesq equations, is presented. This model is widely used for wave modelling and has been well-verified in previous studies, including surf beat studies [36]. The model solves the conservative form of weakly dispersive and fully nonlinear Boussinesq equations in 1-D, following [37]
H t + H U α x + D c = 0
H U α t + H U α 2 + g H η x + g H D x + U α D c R b x = 0
where H = η + h is the total water depth, η is the surface elevation, h water depth, U α velocity at water depth z α . D c and D x are the higher order terms that include bottom turbulence and dispersive properties and R b x is the wave breaking dissipation term. The simulations are performed introducing the radiation stress forcing term into the momentum Equation (12), using the right hand term in Equation (1), as the driving mechanism of bound wave generation, with radiation of free waves (4) as required to ensure continuity. Following Nielsen and Baldock [8], seven cases of Gaussian-shaped wave group envelopes, similar to Equation (8), with the same primary mean wave period of T = 6   s (and the same group velocity, Equation (9)), but different group lengths, varying from 200 m to 345 m , are tested. The domain starts with a long section of horizontal bottom ( h = 20   m ) connected to a slope ( β = 0.025 ) followed by a long plateau ( h = 5   m ), similar to Figure 3b. The length of the plateau is sufficient for the free waves, generated during the shoaling process, to propagate completely away from the forced waves so that the different waves are separated in space and time. The initial condition is flat water, with the forced wave at the initial time-step balanced by the two free waves as proposed by Nielsen and Baldock [8]. The length of the initial horizontal domain is also sufficient so the free waves propagate completely away from the forcing region, such that when the forcing zone propagates to the start of the sloping beach, only the forced wave remains as a pure depression, representing the conditions proposed in the model of Mei [2] (Figure 1b).
Similar to the results from the conceptual model, Figure 5, these specific examples, Figure 6, show that in the shoaling zone the shorter infragravity waves have larger apparent amplitudes than the longer infragravity waves. Also, the shorter the long wave, the greater is the apparent lag with respect to the forcing. However, after the forcing reaches the constant depth plateau, the shoaling stops and the free waves propagate away ahead of the forcing region, leaving behind only forced waves, which all have with same amplitude. The differences in amplitude of the leading free waves at around t ≈ 800 s on the plateau, Figure 6b, are due to the different lengths of the forcing relative to the length of the slope where the changes are occurring, with shorter groups resulting in a smaller free wave amplitude.
The numerical calculations for these same specific examples also show that for all cases no lag occurs in the spatial domain between the forcing and the bound wave (see Figure 7) but that a lag is observed in the time domain (Figure 8), similar to that observed in the laboratory data (Figure 3). The difference between observations in the spatial domain and the time domain is important, since most measurements of the bound wave lag occur from data obtained in the time-domain. In this model, the forced wave part of the total wave is by definition in antiphase with the forcing, so it does not lag the forcing. However, since the system is unsteady and evolving as the wave passes a fixed location, the total solution has an apparent lag in the time domain. Overall, the results suggests that when considering the approach of Mei [2] (Figure 1b) for the bound wave problem via radiation stress forcing, both the shape and shoaling rates could be affected by the generation of positive leading free waves required to balance the shoaling forced wave. Further investigation is needed in order to verify if free waves generation occurs in laboratory or natural conditions. This is difficult, since long spatial domains are required to observe separation of free and forced waves with very similar celerity. Indeed, in shallow water the waves do not separate at all.

3. Bound Wave Dissipation after Short Wave Breaking

One important question yet to be fully answered is what happens to the bound wave after short wave breaking? Different datasets indicate bound wave release [12,13], while others indicate dissipation of the forced wave in the surf zone [17,24,36,38]. Baldock [15] suggested that the relationship between short and infragravity wave with water depth at the breakpoint is important to determine whether the bound wave is released or if it remains forced, decaying with the dissipation of the short waves. During the breaking process, three different conditions are possible: (i) at the breaking point both the short and forced waves are not in shallow water, c g < g h and the length of the bound wave is not much longer than h; (ii) the bound wave is much longer than h, but the short waves are not, hence the bound wave number differs from the wave number of a free wave with same frequency; (iii) both short and bound wave are in shallow water at the break point. Bound wave release is considered to occur when both short and bound waves are in shallow waters, therefore being strictly valid only for the third condition. While the first possibility is more likely for deep-water wave breaking, the second condition can occur at the edge of the surf zone, especially during storm conditions [15].
Baldock and O’Hare [39] compared the behaviour of free non-breaking waves and forced infragravity waves inside the surf zone. The laboratory experiment was designed so the effects of short wave breaking on forced and free waves could be evaluated. With one set of cases based on two free waves, a short breaking wave and a smaller and longer non-breaking wave; and a second set of cases based on fully modulated bichromatic groups. The results have shown that while small or nearly no dissipation occurred for the free non-breaking waves, strong dissipation occurred for short bound waves starting at the break point. The authors then observed that the short wave breaking was occurring when the waves were not in shallow water, justifying the forced wave decay in the surf zone. The numerical investigation presented in Moura and Baldock [36] showed, for the same cases, that no significant breakpoint excursion occurs at the group frequencies, therefore there was no breakpoint forcing at these frequencies. This removes the possibility of interference between the bound wave and breakpoint force long waves, which has been previously suggested as an explanation for the infragravity wave amplitude decay in the surf zone [3].
In order to test whether the breaking depth affects the behaviour of the bound wave, the simulations using the COULWAVE model for two fully-modulated bichromatic cases presented by Moura and Baldock [36] are extended to a larger range of wave conditions, where the amplitude of the fully-modulated wave groups are progressively reduced, hence shifting the break point to shallower water. The COULWAVE model has two options to simulate wave breaking, the eddy viscosity model described in Kennedy et al. [40] and an adaptation to this model, named the transport-based model [41], which is the one selected for this study. Both models use the same eddy viscosity dissipation term and a triggering mechanism (local time derivative of the surface elevation exceeding a threshold) to initiate wave breaking. For further details of the modelling see Moura [42]. For both cases, strong bound wave dissipation, starting at the edge of the surf zone, have been observed in the data [24] and reproduced numerically [36]. The numerical results in Figure 9 indicate that independently of the water depth at the short wave breakpoint the bound wave is dissipated from the breakpoint up to the swash zone. The small outgoing waves corroborate that the bound wave dissipation occurs for all tested cases. This means that when short wave breaking occurs in shallow water, and hence bound wave release is expected, strong dissipation is still occurring. This implies that the energy dissipation is due to turbulence generated by the short waves. This supports the suggestion of Battjes et al. [13], who proposed that infragravity waves that are relatively short compared to the surf zone width are more easily dissipated than longer waves whose its lengths are only a fraction of the surf zone. Hence, dissipation is also strongly frequency dependent.
To extend this analysis further, using the COULWAVE model, and following the experiment layout in Baldock and O’Hare [39], the effects of a surf zone on the dissipation of non-breaking free long waves is numerically investigated by two sets of simulations, with and without a surf zone. The first set of simulations with a surf zone comprise of a short breaking monochromatic wave (w1) with an amplitude of 5 cm and frequency of 0.8 Hz paired with longer free monochromatic waves (w2) with an amplitude of 0.05 cm and frequencies varying from 0.09 to 0.5 Hz, to represent free long waves inside a surf zone. The second set of simulations comprised of the same monochromatic free long waves (w2) but with no short waves and hence no surf zone. The model results for the long waves (w2) with and without a surf zone are then compared in terms of total, incident and outgoing wave amplitude.
The standing wave total amplitude patterns generated by w2 in the absence or presence of a surf zone are significantly different (Figure 10). While for both scenario a reduction of the nodal structures occurs towards the higher frequencies, the process is intensified for the cases with a surf zone present. Thus, dissipation of free long waves by eddy turbulence inside a surf zone is certainly frequency dependent based on the COULWAVE model, and it is likely this explains the decay of the incident forced wave also, as proposed by Battjes et al. [13], with some additional reduction in amplitude since the forcing is also being reduced by short wave breaking.

4. Discussion

There are two interpretations of the forced wave in the literature, as indicated in Figure 1. The model proposed by Mei [2] appears to require the radiation of free waves if the forced wave is to change amplitude by shoaling, which is suggested by Nagase and Mizuguchi [21] and Nielsen and Baldock [8]. However, radiation of such free waves during shoaling has yet to be observed with any degree of certainty. The model of Longuet-Higgins and Stewart [1] does not appear to require such radiation of free waves to ensure continuity, since the surface elevation has bot negative and positive components, both of which can grow. This study has adopted the conceptual model approach of Mizuguchi and Nagase and Nielsen and Baldock and constructed a simple analytical model to investigate the shoaling behavior of the total long wave if forced and radiating free waves co-exist. The results clearly show an apparent growth of the long wave at rates less than that proposed by Longuet-Higgins and Stewart (Figure 5), which has been observed in many data sets, and theoretically analysed by alternate means by Janssen et al. [12] and Battjes et al. [13] The growth rate and lag of the long wave behind the forcing are frequency dependent, in agreement with the mild and steep slope shoaling regimes proposed by Battjes et al. However, this analysis shows that the lag only occurs in the time domain, which is how most data is collected, but not in the spatial domain. The same qualitative results were found when performing simulations with the COULWAVE model using the radiation stress term in the momentum equation to force the generation and propagation of bound waves. This is again supported by the limited data sets available, see Figure 3. With this model approach, the forced waves (not the total wave) follow the equilibrium solution of Longuet-Higgins and Stewart, and the shoaling of the forced wave part is not frequency dependent. We note that radiation of free waves during shoaling has yet to be verified in observations, and further work with RANS models is ongoing to determine if this is the case.
We have also used the COULWAVE to extend experimental data sets that investigated the dissipation of the incident forced and free long waves in the surf zone. The results clearly suggest that the dissipation rate is dependent on the long wave frequency or length, for both forced and free waves, in agreement with Battjes et al. [13] irrespective of the depth at which short wave break. This suggests that the reduction in forced wave amplitude inside the surf zone is a combination of the reduced forcing, as originally proposed by Longuet-Higgins and Stewart [1] and further analysed by Baldock [15], and turbulent dissipation as suggested by Battjes et al [13].

Author Contributions

T.M. and T.E.B. conceived and designed the modelling; T.M. performed the modelling and analysed the results; T.M. and T.E.B. wrote the paper.

Funding

Theo Moura gratefully acknowledges a scholarship from the Brazilian National Council for Scientific and Technological Development (CNPq). This work was supported by Australian Research Council grants LP100100375 and DP13101122.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Longuet-Higgins, M.S.; Stewart, R.W. Radiation stress and mass transport in gravity waves, with application to ‘surf beats’. J. Fluid Mech. 1962, 13, 481–504. [Google Scholar] [CrossRef]
  2. Mei, C.C. The Applied Dynamics of Ocean Surface Waves; World Scientific: Singapore, 1989; p. 740. [Google Scholar]
  3. Schäffer, H.A. Infragravity waves induced by short-wave groups. J. Fluid Mech. 1993, 247, 551. [Google Scholar] [CrossRef]
  4. Proudman, J. The Effects on the Sea of Changes in Atmospheric Pressure. Geophys. J. Int. 1929, 2, 197–209. [Google Scholar] [CrossRef]
  5. Whitham, G.B. Mass, momentum and energy flux in water waves. J. Fluid Mech. 1962, 12, 135–147. [Google Scholar] [CrossRef]
  6. Molin, B. On the generation of long-period second-order free-waves due to changes in the bottom profile. Ship Res. Inst. Pap. 1982, 1–23. [Google Scholar]
  7. Mei, C.C.; Benmoussa, C. Long waves induced by short-wave groups over an uneven bottom. J. Fluid Mech. 1984, 139, 219. [Google Scholar] [CrossRef]
  8. Nielsen, P.; Baldock, T.E. И-Shaped surf beat understood in terms of transient forced long waves. Coast. Eng. 2010, 57, 71–73. [Google Scholar] [CrossRef]
  9. Tinti, S.; Bortolucci, E.; Chiavettieri, C. Tsunami Excitation by Submarine Slides in Shallow-water Approximation. Pure Appl. Geophys. 2001, 158, 759–797. [Google Scholar] [CrossRef]
  10. Didenkulova, I.; Nikolkina, I.; Pelinovsky, E.; Zahibo, N. Tsunami waves generated by submarine landslides of variable volume: Analytical solutions for a basin of variable depth. Nat. Hazards Earth Syst. Sci. 2010, 10, 2407–2419. [Google Scholar] [CrossRef]
  11. Baldock, T.E. Long wave generation by the shoaling and breaking of transient wave groups on a beach. Proc. R. Soc. A Math. Phys. Eng. Sci. 2006, 462, 1853–1876. [Google Scholar] [CrossRef]
  12. Van Dongeren, A.R.; Janssen, T.T.; Battjes, J.A. Long waves induced by short-wave groups over a sloping bottom. J. Geophys. Res. Space Phys. 2003, 108, 108. [Google Scholar]
  13. Battjes, J.A.; Bakkenes, H.J.; Janssen, T.T.; Van Dongeren, A.R. Shoaling of subharmonic gravity waves. J. Geophys. Res. Space Phys. 2004, 109, 109. [Google Scholar] [CrossRef]
  14. Dong, G.; Ma, X.; Perlin, M.; Ma, Y.; Yu, B.; Wang, G. Experimental study of long wave generation on sloping bottoms. Coast. Eng. 2009, 56, 82–89. [Google Scholar] [CrossRef]
  15. Baldock, T.E. Dissipation of incident forced long waves in the surf zone—Implications for the concept of “bound” wave release at short wave breaking. Coast. Eng. 2012, 60, 276–285. [Google Scholar] [CrossRef]
  16. Padilla, E.M.; Alsina, J.M. Transfer and dissipation of energy during wave group propagation on a gentle beach slope. J. Geophys. Res. Oceans 2017, 122, 6773–6794. [Google Scholar] [CrossRef] [Green Version]
  17. Padilla, E.M.; Alsina, J.M. Long Wave Generation Induced by Differences in the Wave-Group Structure. J. Geophys. Res. Oceans 2018, 123, 8921–8940. [Google Scholar] [CrossRef]
  18. Van Dongeren, A.R.J.A.; Battjes, J.; Janssen, T.; Van Noorloos, J.; Steenhauer, K.; Steenbergen, G.; Reniers, A.J.H.M. Shoaling and shoreline dissipation of low-frequency waves. J. Geophys. Res.-Oceans 2007, 112, C02011. [Google Scholar] [CrossRef]
  19. Bowen, A.J.; Henderson, S.M. Observations of surf beat forcing and dissipation. J. Geophys. Res. Space Phys. 2002, 107, 107. [Google Scholar]
  20. Raubenheimer, B.; Thomson, J.; Elgar, S.; Herbers, T.H.C.; Guza, R.T. Tidal modulation of infragravity waves via nonlinear energy losses in the surfzone. Geophys. Res. Lett. 2006, 33. [Google Scholar] [CrossRef] [Green Version]
  21. Nagase, S.; Mizuguchi, M. Laboratory experiment on long wave generation by time-varying breakpoint. Coast. Eng. 1997, 1307–1320. [Google Scholar]
  22. Masselink, G. Group bound long waves as a source of infragravity energy in the surf zone. Cont. Shelf Res. 1995, 15, 1525–1547. [Google Scholar] [CrossRef]
  23. Madsen, P.A.; Sørensen, O.; Schäffer, H. Surf zone dynamics simulated by a Boussinesq type model. Part II: Surf beat and swash oscillations for wave groups and irregular waves. Coast. Eng. 1997, 32, 289–319. [Google Scholar] [CrossRef]
  24. Baldock, T.E.; Huntley, D.A.; Bird, P.A.D.; O’hare, T.; Bullock, G.N. Breakpoint generated surf beat induced by bichromatic wave groups. Coast. Eng. 2000, 39, 213–242. [Google Scholar] [CrossRef]
  25. Baldock, T.E.; Huntley, D.A. Long-wave forcing by the breaking of random gravity waves on a beach. Proc. R. Soc. Lond. A. Phys. Eng. Sci. 2002, 458, 2177–2201. [Google Scholar] [CrossRef]
  26. Orszaghova, J.; Taylor, P.H.; Borthwick, A.G.; Raby, A.C. Importance of second-order wave generation for focused wave group run-up and overtopping. Coast. Eng. 2014, 94, 63–79. [Google Scholar] [CrossRef] [Green Version]
  27. Mansard, E.; Barthel, V. Shoaling properties of bounded long waves. Coast. Eng. 1984, 1985, 798–814. [Google Scholar]
  28. List, J.H. A model for the generation of two dimensional surf beat. J. Geophys. Res. 1992, 97, 5623–5635. [Google Scholar] [CrossRef]
  29. Longuet-Higgins, M.S.; Stewart, R. Changes in the form of short gravity waves on long waves and tidal currents. J. Fluid Mech. 1960, 8, 565–583. [Google Scholar] [CrossRef]
  30. Phillips, O.M. The Dynamics of the Upper Ocean; CUP Archive, 1966. [Google Scholar]
  31. Henderson, S.M.; Guza, R.T.; Elgar, S.; Herbers, T.H.C. Refraction of Surface Gravity Waves by Shear Waves. J. Phys. Oceanogr. 2006, 36, 629–635. [Google Scholar] [CrossRef] [Green Version]
  32. Guedes, R.M.; Bryan, K.R.; Coco, G.; Bryan, R.; Coco, G. Observations of wave energy fluxes and swash motions on a low-sloping, dissipative beach. J. Geophys. Res. Oceans. 2013, 118, 3651–3669. [Google Scholar] [CrossRef]
  33. Bowers, E.C. Low Frequency Waves in Intermediate Water Depths. In Proceedings of the 23rd International Conference on Coastal Engineering, Venice, Italy, 4–9 October 1992; pp. 832–845. [Google Scholar]
  34. Van Leeuwen, P.J. Low Frequency Wave Generation Due to Breaking Wind Waves. Ph.D. Thesis, Delft University of Technol, Delft, The Netherlands, 1992. [Google Scholar]
  35. Lynett, P.J. Nearshore wave modeling with high-order Boussinesq-type equations. J. Waterw. Port Coast. Ocean Eng. 2006, 132, 348–357. [Google Scholar] [CrossRef]
  36. Moura, T.; Baldock, T.E. New evidence of breakpoint forced long waves: Laboratory, numerical, and field observations. J. Geophys. Res. Oceans. 2018, 123, 2716–2730. [Google Scholar] [CrossRef]
  37. Kim, D.-H.; Lynett, P.J.; Socolofsky, S.A. A depth-integrated model for weakly dispersive, turbulent, and rotational fluid flows. Ocean Model. 2009, 27, 198–214. [Google Scholar] [CrossRef]
  38. Contardo, S.; Symonds, G. Infragravity response to variable wave forcing in the nearshore. J. Geophys. Res. Oceans 2013, 118, 7095–7106. [Google Scholar] [CrossRef]
  39. Smith, J.M.; Baldock, T.E.; O’Hare, T.J. Energy transfer and dissipation during surf beat conditions. Coast. Eng. 2005, 4, 1212–1224. [Google Scholar]
  40. Kennedy, A.B.; Chen, Q.; Kirby, J.T.; Dalrymple, R.A. Boussinesq Modeling of Wave Transformation, Breaking, and Runup. I: 1D. J. Waterw. Port Coast. Ocean Eng. 2000, 126, 39–47. [Google Scholar] [CrossRef]
  41. Løvholt, F.; Lynett, P.; Pedersen, G. Simulating run-up on steep slopes with operational Boussinesq models; capabilities, spurious effects and instabilities. Nonlinear Process. Geophys. 2013, 20, 379–395. [Google Scholar] [CrossRef] [Green Version]
  42. Moura, T. Infragravity Wave Forcing in the Surf and Swash Zone; School of Civil Engineering, the University of Queensland: Brisbane, Australia, 2016. [Google Scholar]
Figure 1. Two different graphical interpretations of the bound wave forced by a short wave groups (black line) by Longuet-Higgins and Stewart [1] with positive and negative elevations (red line), and Mei [2] with purely negative elevation (green line).
Figure 1. Two different graphical interpretations of the bound wave forced by a short wave groups (black line) by Longuet-Higgins and Stewart [1] with positive and negative elevations (red line), and Mei [2] with purely negative elevation (green line).
Jmse 07 00305 g001
Figure 2. (a) Spatial visualization of short wave envelope (dashed black line), total long wave surface elevation (grey dashed), long wave with setup removed (full black line) and mean setup (full grey line), over bottom topography shown (solid black line). (b) Time evolution of short wave envelope (dashed line) and total long wave (full line) at x = −3.75 m (location A in left panel). Random wave case J6033A from Baldock and Huntley [25].
Figure 2. (a) Spatial visualization of short wave envelope (dashed black line), total long wave surface elevation (grey dashed), long wave with setup removed (full black line) and mean setup (full grey line), over bottom topography shown (solid black line). (b) Time evolution of short wave envelope (dashed line) and total long wave (full line) at x = −3.75 m (location A in left panel). Random wave case J6033A from Baldock and Huntley [25].
Jmse 07 00305 g002
Figure 3. (a) Schematic representation of the surface elevation at a fixed position for a superposition of forced and free waves. Dot-dashed line - bound wave propagating with group velocity and amplitude A f o r c e d . Dashed line - free wave generated during the shoaling process propagating with with respective amplitude A f r e e + . Full line - total signal with amplitude A t o t a l . (b) Spatial representation of bound wave shoaling. Black line, superposition of incident bound wave (IBW) and an incident forward free wave (IFW) generated to balance changes in the forced solution during shoaling. Dotted-dashed line is the representation of the radiation stress forcing.
Figure 3. (a) Schematic representation of the surface elevation at a fixed position for a superposition of forced and free waves. Dot-dashed line - bound wave propagating with group velocity and amplitude A f o r c e d . Dashed line - free wave generated during the shoaling process propagating with with respective amplitude A f r e e + . Full line - total signal with amplitude A t o t a l . (b) Spatial representation of bound wave shoaling. Black line, superposition of incident bound wave (IBW) and an incident forward free wave (IFW) generated to balance changes in the forced solution during shoaling. Dotted-dashed line is the representation of the radiation stress forcing.
Jmse 07 00305 g003
Figure 4. Variation in total infragravity wave amplitude at a fixed position on the virtual slope for wave groups of different frequency using the conceptual model of bound and free wave superposition illustrated in Figure 3. Wave groups have the same group velocity but varying frequency.
Figure 4. Variation in total infragravity wave amplitude at a fixed position on the virtual slope for wave groups of different frequency using the conceptual model of bound and free wave superposition illustrated in Figure 3. Wave groups have the same group velocity but varying frequency.
Jmse 07 00305 g004
Figure 5. (a) Amplitude variation of total infragravity wave during shoaling for different group frequencies, grey dashed line ( ~ h 5 / 2 ), grey dash-dotted line ( ~ h 1 / 4 ). (b) Variation in theoretical rates of energy transfer calculated based on the observed lags using Equation (7) for different group frequencies.
Figure 5. (a) Amplitude variation of total infragravity wave during shoaling for different group frequencies, grey dashed line ( ~ h 5 / 2 ), grey dash-dotted line ( ~ h 1 / 4 ). (b) Variation in theoretical rates of energy transfer calculated based on the observed lags using Equation (7) for different group frequencies.
Jmse 07 00305 g005
Figure 6. Waves generated by Gaussian-shaped wave group envelopes with different forcing length but with the same propagation speed. Each line represents a different group length, from shortest 200 m (blue line), to the longest 345 m (dark red). Shoaling zone (a) and horizontal plateau (b).
Figure 6. Waves generated by Gaussian-shaped wave group envelopes with different forcing length but with the same propagation speed. Each line represents a different group length, from shortest 200 m (blue line), to the longest 345 m (dark red). Shoaling zone (a) and horizontal plateau (b).
Jmse 07 00305 g006
Figure 7. Spatial visualization of a short (a) and long (b) Gaussian-shape radiation stress (black line) and total long wave surface elevation (red) at different time instants and location over the bottom profile (c). The vertical dashed lines indicate the zero lag between forcing and bound wave.
Figure 7. Spatial visualization of a short (a) and long (b) Gaussian-shape radiation stress (black line) and total long wave surface elevation (red) at different time instants and location over the bottom profile (c). The vertical dashed lines indicate the zero lag between forcing and bound wave.
Jmse 07 00305 g007
Figure 8. Time series of surface displacement (at x = −100 m in Figure 7c), normalized by its absolute minimum value (red) and Gaussian-shape radiation stress normalized by its maximum value (black), for the longer forcing case shown in Figure 7b.
Figure 8. Time series of surface displacement (at x = −100 m in Figure 7c), normalized by its absolute minimum value (red) and Gaussian-shape radiation stress normalized by its maximum value (black), for the longer forcing case shown in Figure 7b.
Jmse 07 00305 g008
Figure 9. Incident (full lines) and outgoing (dash lines) infragravity surface elevation at f g for bichromatic cases 1055 (a), 1060 (b) from Baldock et al. [24]. Colours are for different fully modulated wave groups with short wave amplitude 0.05 m (blue), 0.035 m (red), 0.025 m (green) 0.015 m (black).
Figure 9. Incident (full lines) and outgoing (dash lines) infragravity surface elevation at f g for bichromatic cases 1055 (a), 1060 (b) from Baldock et al. [24]. Colours are for different fully modulated wave groups with short wave amplitude 0.05 m (blue), 0.035 m (red), 0.025 m (green) 0.015 m (black).
Jmse 07 00305 g009
Figure 10. Cross-shore nodal pattern for long waves of different frequencies with (red) and without (black) a surf zone. Frequencies vary from f = 0.195 (a), 0.265 (b) and 0.44 (c) Hz.
Figure 10. Cross-shore nodal pattern for long waves of different frequencies with (red) and without (black) a surf zone. Frequencies vary from f = 0.195 (a), 0.265 (b) and 0.44 (c) Hz.
Jmse 07 00305 g010

Share and Cite

MDPI and ACS Style

Moura, T.; Baldock, T.E. The Influence of Free Long Wave Generation on the Shoaling of Forced Infragravity Waves. J. Mar. Sci. Eng. 2019, 7, 305. https://doi.org/10.3390/jmse7090305

AMA Style

Moura T, Baldock TE. The Influence of Free Long Wave Generation on the Shoaling of Forced Infragravity Waves. Journal of Marine Science and Engineering. 2019; 7(9):305. https://doi.org/10.3390/jmse7090305

Chicago/Turabian Style

Moura, Theo, and Tom E. Baldock. 2019. "The Influence of Free Long Wave Generation on the Shoaling of Forced Infragravity Waves" Journal of Marine Science and Engineering 7, no. 9: 305. https://doi.org/10.3390/jmse7090305

APA Style

Moura, T., & Baldock, T. E. (2019). The Influence of Free Long Wave Generation on the Shoaling of Forced Infragravity Waves. Journal of Marine Science and Engineering, 7(9), 305. https://doi.org/10.3390/jmse7090305

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop