Next Article in Journal
Site Selection of Offshore Solar Farm Deployment in the Aegean Sea, Greece
Previous Article in Journal
Noise of Internal Solitary Waves Measured by Mooring-Mounted Hydrophone Array in the South China Sea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shoaling Wave Shape Estimates from Field Observations and Derived Bedload Sediment Rates

by
Tarandeep S. Kalra
1,*,
Steve E. Suttles
2,
Christopher R. Sherwood
2,
John C. Warner
2,
Alfredo L. Aretxabaleta
2,* and
Gibson R. Leavitt
1,3
1
Aperture Inc., Contractor to the U.S. Geological Survey, Herndon, VA 20171, USA
2
U.S. Geological Survey, Woods Hole, MA 02543, USA
3
Contractor to the U.S. Geological Survey, University of North Carolina, Chapel Hill, NC 27514, USA
*
Authors to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2022, 10(2), 223; https://doi.org/10.3390/jmse10020223
Submission received: 30 December 2021 / Revised: 28 January 2022 / Accepted: 30 January 2022 / Published: 8 February 2022
(This article belongs to the Section Coastal Engineering)

Abstract

:
The shoaling transformation from generally linear deep-water waves to asymmetric shallow-water waves modifies wave shapes and causes near-bed orbital velocities to become asymmetrical, contributing to net sediment transport. In this work, we used two methods to estimate the asymmetric wave shape from data at three sites. The first method converted wave measurements made at the surface to idealized near-bottom wave-orbital velocities using a set of empirical equations: the “parameterized” waveforms. The second method involved direct measurements of velocities and pressure made near the seabed: the “direct” waveforms. Estimates from the two methods were well correlated at all three sites (Pearson’s correlation coefficient greater than 0.85). Both methods were used to drive bedload-transport calculations that accounted for asymmetric waves, and the results were compared with a traditional excess-stress formulation and field estimates of bedload transport derived from ripple migration rates based on sonar imagery. The cumulative bedload transport from the parameterized waveform was 25% greater than the direct waveform, mainly because the parameterized waveform did not account for negative skewness. Calculated transport rates were comparable to rates estimated from ripple migration except during the largest event, when calculated rates were as much as 100 times greater, which occurred during high period waves.

1. Introduction

Waves propagating towards the coast undergo transformations from their sinusoidal wave shape in deep water to non-linear wave shape in shallow water. “Skewed”, as used here, describes wave shapes, bottom-orbital velocities, and sediment transport generated by waves that cannot be adequately represented by linear (Airy) wave theory (e.g., [1]). The fundamental equations for continuity of mass and momentum used as the starting point for most wave theories are linear, but the surface boundary conditions are not. Linear wave theory neglects the non-linear components by assuming that the ratio of wave height/wavelength (wave steepness) is small. As waves shoal and steepen, this assumption no longer holds, and non-linear equations are needed to accurately represent the waves [1]. Outside the surf zone, this shoaling transformation causes the near-bottom wave orbital velocity to differ between the wave crest and trough cycles, referred to as velocity skewness [2]. As waves propagate further into the surf zone, the differential duration in crest and trough half cycles can lead to a differential acceleration in each half cycle, referred to as acceleration asymmetry [3].
The effect of shoaling waves on resulting sediment dynamics has been quantified in various studies. In previous flume experiments [4,5,6], skewed waves led to onshore sediment transport near the bed. Under mild wave conditions, velocity skewness has been found to be the main driver of onshore sediment transport, generating onshore sandbar migration [7,8,9]. The combined effect of different asymmetric wave shapes in the presence of currents acting against the direction of waves and showed that under decreasing skewed-asymmetric wave shapes, net offshore sediment transport increased [10]. Onshore sediment transport has been also related to ripple migration associated with velocity skewness in several field observations [11,12,13]. Other than flume and field data, numerical models [14,15,16] have shown that acceleration asymmetry caused high onshore sediment transport and contributed towards beach recovery. In addition, process-based morphodynamic modeling studies have shown the impact of wave-induced cross-shore sediment transport that caused inlet closure during weak inlet current conditions [17,18,19]. In another study [20], a numerical model in an idealized inlet domain that accounted for asymmetric waves showed that larger waves contributed to the migration and growth of shoals on ebb-tidal deltas. Using a realistic domain [21], it was shown that the skewness of the waves was responsible for the enhancement of sediment transport into the inlet and the evolution of the ebb-tidal delta during storm conditions. Recent modeling [22,23,24,25,26] studies have shown the importance of wave-driven transport under combined action during wave-current interaction.
These studies highlight the importance of estimating the asymmetric waveforms that affect nearshore sediment transport. Other efforts [27,28] emphasized that understanding wave-driven sediment transport is fundamental to the improvement of morphodynamic models. In the presence of waves, morphodynamic models need to account for both suspended and bedload transport. Of the two modes of sediment transport, bedload transport is most affected by near-bottom wave orbital wave motions [29], thus necessitating the need to estimate asymmetric waveforms in modeling frameworks. Asymmetric waveforms can be computed by wave-resolving models [30,31], but these models are computationally expensive and impractical for regional-scale morphodynamic models. One approach to estimating asymmetric waveforms used a statistical method [32] that applied a genetic algorithm to obtain velocity skewness based on field observations made under shoaling and breaking wave conditions. Another approach to quantify asymmetry from shoaling waves is to use the bulk wave statistics and local depth [33,34]. The approach mentioned in [34] would be referred to as Ruessink’s method. A recent study [35] compared the modified version of the approach mentioned in [33] to Ruessink’s method [34] to simulate long-term morphodynamic change. They [35] found that Ruessink’s method provided a better representation of wave and near-bed orbital velocities, required less calibration and, therefore, resulted in more realistic morphodynamic predictions. A modification to Ruessink’s approach was used in a study [36] to accelerate the model calculations to capture the effects of nourishment on cross-shore bedload fluxes at different locations on an idealized planar beach. Ruessink’s parameterization also showed a good skill predicting shoreward evolution of wave skewness and asymmetry under high-energy waves ( H s ~   2   m ) during offshore sandbar migration events [37].
In the present work, we analyzed the efficacy of the widely used Ruessink’s parameterization by taking observations from three sites located outside the surf zone. At these sites, conditions (including depth and wave characteristics) only generate velocity skewness in the waveform shape; acceleration skewness becomes more important in shallower water. The three sites were located at Fire Island, New York, Martha’s Vineyard Coastal Observatory, Massachusetts, and Matanzas Inlet, Florida. The efficacy of Ruessink’s method was determined by parameterized estimates computed from bulk surface-wave statistics with direct observations of skewed waveforms made near the bottom. While the parameterizations depended on bulk wave statistics (significant wave height, representative wave period) and water depth, and were translated to the bottom using linear wave theory, the direct observations were obtained from measurements of near-bottom velocities, as described later in this work. The characteristics of the skewed waveforms from both methods include: crest cycle velocity, ( u c ) , trough cycle velocity, ( u t ) , crest cycle period, ( T c ) , trough cycle period ( T t ) , and accelerative periods of crest and trough cycle, ( T c u ) and ( T t u ) , respectively. Following this comparison, the skewed waveforms obtained from the parameterized and direct methods were used to drive a bedload-transport formulation that accounted for asymmetric waveforms [38]. We compared the bedload transport results with the widely used Meyer–Peter and Müller (MPM) equation [39] and bedform-transport rates from sonar data that tracked the ripple migration on the seafloor. The bedload-transport rate comparisons provide an insight into the implementation of asymmetric bedload formulation in morphodynamic models that would rely on Ruessink’s method to first obtain the asymmetric waveform in a way that is analogous to the “parameterized” calculations described in the present observational work. The paper is organized as follows: The dataset, instrumentation, and analytical methods used to obtain the skewed waveforms from parameterized and direct methods along with the methods used to obtain bedload are described in Section 2. The results are presented in Section 3, followed by a discussion in Section 4. Section 5 presents the conclusions.

2. Materials and Methods

2.1. Field Study Sites

The datasets used in the present work were obtained at three field study sites: Fire Island, New York (FI), Martha’s Vineyard Coastal Observatory, MA (MVCO), and Matanzas Inlet, FL (MI). The FI dataset was part of a field campaign [40,41] to study the along-shore variations in waves and sediment processes. FI data from the site shown in Figure 1a correspond to a mean depth of 16 m located approximately 2 km offshore from the southern shoreline of western Fire Island in New York. The MVCO dataset [42] was collected at a site with an average depth of 14 m offshore located approximately 2 km from the southern shoreline of Martha’s Vineyard (Figure 1b). The Matanzas Inlet dataset was collected [43] to study the effects of cross-shore currents on waves, and measurements were made at a site with an average depth of 9 m located approximately 0.8 km offshore from the eastern shoreline (Figure 1c). Table 1 mentions the filenames containing the datasets used to obtain near bottom and surface wave measurements.

2.2. Parameterized Method: Skewed Waveform Estimation from Surface Wave Data

The parameterized method to compute the near-bottom wave forms uses surface bulk wave parameters. The parameterized method to find skewed waveform characteristics ( u c r e s t , u t r o u g h , T c , T t , T c u , T t u ) is described in Appendix A and requires local values of significant wave height H S , a representative wave period T , depth h , and wave orbital velocity amplitude u w . Surface wave spectra were measured with an acoustic Doppler current profiler (ADCP; Table 1). The ADCP beam orbital velocities measured at bins near the surface were used to calculate wave spectra at the three sites using the WavesMon software. An exception to this methodology was the Matanzas Inlet site where, after the first 31 days of the 80-day deployment, measurements from a near-bed pressure sensor (SeaBird Seagauge) were substituted because the ADCP malfunctioned. In these pressure-sensor data, the mean H S was 7% smaller and the mean T p was 8% larger than those measured with the ADCP for the 31 days when both the sensors worked. The methods described by [44] were used to obtain near-bed estimates of u w and T from the surface wave spectra. Note that for the application of Appendix A, we used a representative wave period ( T ) that was defined as T  = m 0 m 1 , where m 0 is the variance of water surface elevation and m 1 is the first moment of the spectral density of surface elevation. Table 2 provides sensor specifications used from each of the three sites. The characteristics of the parameterized waveform obtained using this method are described later with the use of subscript “( p )”.

2.3. Direct Method: Skewed Waveform Estimation from Near-Bed Data

The direct method uses near-bed observations to determine a representative skewed waveform. Waveforms were observed directly from near-bed velocity measurements made with acoustic Doppler velocimeters (ADVs; Table 3). The measured time series of wave velocity directly obtained from ADV recordings are in two components, u r a w and v r a w close to the seafloor in bursts (ref. Table 3 for site-specific burst length) at an elevation close to the seafloor but above the wave boundary layer. The process is described here for one component, that is, u r a w for the sake of brevity; a similar analysis was performed for v r a w .   u r a w was decomposed into a slowly varying mean tidal component u , an incident-wave band u ˜ i n c i d e n t , an infragravity wave band u ˜ i n f r a g r a v i t y and noise associated with turbulence and instrument error u .
u r a w = u + u ˜ i n c i d e n t + u ˜ i n f r a g r a v i t y + u
Velocities from each burst time series were first detrended to remove u . Incident-wave velocities u ˜ i n c i d e n t with periods of 4–20 s were isolated using a Butterworth bandpass filter. The lower limit of 4 s removed higher-frequency motions that were assumed to not affect the calculations of bedload at the bottom, thus reducing the noise in the signal. The upper limit of 20 s excluded waves in the infragravity band. The infragravity band was removed so that only the sea-swell frequency band remained, to match the frequency range in Ruessink’s data [34] and to apply an asymmetric bedload-transport equation that was formulated with the exclusion of the infragravity band. Figure 2 compares the raw and incident wave spectrum from a storm event in the Fire Island dataset that corresponds with an H S of 3.7 m associated with a corresponding T P of 11.6 s.
A principal components analysis (PCA) was applied to the u ˜ i n c i d e n t and v ˜ i n c i d e n t wave signals to determine the axis of dominant wave motions. PCA determined the dominant wave direction and the corresponding velocity signal (referred to as u ˜ ) for each burst. Representative waveforms were found within burst time series by (a) identifying individual waves using zero-up crossings, (b) determining characteristics ( u c r e s t ,   u t r o u g h , T c r e s t , T t r o u g h , T c u , and T t u ) of the individual waveforms, (c) calculating parameters that describe waveform according to the formula of [3], and computing a weighted average of those parameters. Zero up-crossings times were used to find short time series of velocity u ( t ) with period T i for each wave in the time series, where subscript “ i ” refers to the individual waveforms in a wave burst time series. From these individual wave motions, the crest and trough velocities u c , i ,   u t , i and crest and trough periods T c , i and T t , i were determined. In the next step, three parameters corresponding to the individual waveforms were calculated that included the amplitude of wave orbital velocity as defined by:
u w , i = 0.5   ( u c ,   i u t ,   i )  
and their velocity skewness ( R i ) according to [45]:
R i = u c ,   i   ( u c ,   i u t ,   i )  
The inputs of non-dimensional u w , i , T i and R i were used to calculate the third parameter b 1 , i a fitting parameter that results in a free-stream velocity description [3] by finding the maximum and minimum positions of velocity in the individual waveforms. The equations to calculate b 1 , i are provided in the Matlab script “abreu_fit” provided by [46]. The weighted means of individual waveform u w , i , T i , R i , and b 1 , i were used in the Matlab script “abreu_fstream” [46] to obtain the representative waveform for each burst time series (Figure 2b as an example). Weighting was based on the product of the squared wave-orbital amplitude and wave period, as:
w i = u w , i 2 T i  
The characteristics of the parameterized waveform obtained using this method are described later with the use of subscript “ d ”.

2.4. Calculation of Bedload for Skewed Waveforms

We derived three estimates of wave-driven cross-shore bedload transport using the Matanzas Inlet dataset because that site had the highest proportion of skewed waveforms (Section 3.4). Both the parameterized and direct skewed waveforms were used to drive a bedload transport formulation that accounted for skewed waves [38], referred to as asymmetric bedload formulation. Results were also calculated using a traditional excess-stress bedload formulation, the Meyer–Peter and Müller (MPM) equation [39]. Appendix B and Appendix C detail the equations used for the asymmetric and MPM formulations. The grain size used in the calculations was d 50 = 0.23   mm [43]. In contrast to the asymmetric bedload formulations, which represent waveforms with wave asymmetry indices, the MPM formulation (Equation (A43)) uses the entire time series of instantaneous wave motions ( u ˜ i n c i d e n t ) in the bedload calculations.

2.5. Sonar Imagery for Bedform Transport Rate from Ripple Migration

Measurements of bedform transport, a proxy for bedload transport, were determined from sonar imagery that tracked ripple migration on the seafloor. The analysis was performed for the shallowest site (Matanzas Inlet) where good images of active ripple migration were available. Images of the seabed were obtained every 30 min from a 2.25 MHz Imagenex model 881A fan-beam rotating imaging sonar mounted 0.45 m (nominal) above the bed. The raw images of backscatter intensity from the sonar were corrected for slant range and rotated into geographical coordinates (positive x = east). This resulted in images with a horizontal resolution of 0.01 m and covering a circular region with a 5 m radius (Figure 3a). In these images, the light and dark pixel values corresponded to acoustic backscatter intensity (light indicates high backscatter that revealed the presence of ripples). Manual inspection was used to find a region with ripples that were well defined for most of the time series. This corresponded to a 1-m × 1-m square (Figure 3a) orientated in the onshore direction located 1–2 m from the center of the image. The ripple wavelength ( λ r ) was determined in the sub-patch region using the two-dimensional fast Fourier transform method described in these studies [47,48,49,50].
We used a single cross-shore transect (shown in Figure 3a) located at 0.05 m distance from the southern edge of the square to form a time-stacked image (Figure 3b). Manual trials showed that the quality of the time stacked image remained unchanged for ripple tracking purposes by choosing any other current transect within the square region. The cross-shore position of pixel intensity was determined using edge detection [51,52]. The cross-shore migration rate ( d r ) was determined as the change in cross-shore position at half-hour intervals. The resulting volumetric bedform transport rates ( q b ,   m 3 / m ) estimated from the sonar data were computed [11] as:
q b = ( 1 ε ) f r η r d r
assuming a bed porosity factor ε = 0.4 , dimensionless ripple shape factor f r = 0.5 , and a constant ripple height to ripple wavelength ratio η r λ r = 0.16 [11].

3. Results

3.1. Surface Wave Statistics

MVCO had the least-energetic waves of the three sites (Figure 4). Fire Island was the deepest site and had the highest proportion of high H S . Matanzas Inlet had the highest proportion of waves with high T P . Energetic waves (high H S and T P ) and shallower water depths would be expected to generate larger wave skewness.

3.2. Comparison of Representative Wave Period and Orbital Velocity from Surface Wave Data and Near-Bed Wave Data

Peak wave orbital velocity (uw,d and uw,p) and representative wave periods (Tw,d and Tw,p) are compared in Figure 5 for the three sites, based on both (parameterized and direct) methods described in Section 2.2 and Section 2.3 where subscript “w” refers to the representative waveforms obtained from the two methods (parameterized and direct). The correlation (Pearson’s correlation coefficient, r) between uw,d and uw,p was above 90% at all three sites and the best linear fit (identified by the slope value) was obtained at the shallowest site (Matanzas Inlet; Figure 5a,c,e). At Fire Island and Matanzas Inlet, uw,p was underpredicted compared to uw,d as significant wave heights (Hs) increased (Figure 5a,e). The correlation between Tw,d and Tw,p was found to be above 90% at all three sites (Figure 5b,d,f). Tw,p was consistently underpredicted compared to Tw,d at all three sites for almost all wave heights. This pattern can be attributed to the damping of high-frequency motions in the near-bottom measurements, leading to relatively higher values of Tw,d. The linear fits to Tw,d and Tw,p were better in the shallower depths at MVCO (Figure 5d) and Matanzas Inlet (Figure 5f), compared to those at Fire Island (Figure 5b).

3.3. Skewed Waveform Characteristics: Parameterized and Direct Methods

We used the equations described in Appendix A and methods described in Section 2.2 and Section 2.3 to determine the parameterized and direct representative skewed waveform characteristics which included u c r e s t , u t r o u g h ,   T c , T t , T c u , and T t u . These characteristics were used to calculate the wave-driven bedload calculations (Appendix B). The comparison of u c r e s t , d   and   u c r e s t , p   (Figure 6a) and u t r o u g h , d   and   u t r o u g h , p (Figure 6b) showed a r > 0.9 and that linear fit improved for the Matanzas Inlet site compared to the other two sites. The correlations between T c , d   and   T c , p   ,   T t , d   and   T t , p   ,   T c u , d   and   T c d , p     and   T t u , d   and   T t u , p   remain high ( r > 0.85 ) (Figure 6c–f). The pattern of linear fit for half periods was consistent with the results of the representative wave period (i.e., Figure 5b,d,f).
The E1 error was highest for the shallowest site of Matanzas Inlet for all waveform characteristics (Table 4). E2 error that accounted for scaling the error with the RMS values showed that the Matanzas Inlet site had the least error corresponding to the velocity characteristics of the waveform ( u w ,   u c r e s t , u t r o u g h ). On the contrary, E2 error was highest for period characteristics ( T ,     T c , T t , T c u , and T t u ) at the Matanzas Inlet site.

3.4. Wave Skewness

Wave-orbital velocity skewness drives bedload transport and can be used as a diagnostic tool to indicate transport potential. Non-dimensional wave-orbital velocity skewness (hereafter referred to as skewness, S u ) was calculated at the bottom using the ADV-based direct method as:
S u = u ˜ 3 ( t ) ( u ˜ r m s ) 3
where u ˜ is wave-orbital velocity, angle brackets indicate time-averaging over the burst, and subscript rms indicates root mean square over the burst. S u tends to increase with Ursell number ( U r ) above a threshold of about U r = 0.04 (Figure 7). This is consistent with the observations of [34] who mentioned that U r 0.04 have near-zero skewness. Most of the observations at the two deeper sites (93% at FI and 75% at MVCO) fell below this threshold. In contrast, only 25% of the MI data had a U r 0.04 .
For U r > 0.04 , about 22% of the waveforms at MI were associated with negative skewness (similar to the results in [34]) (Figure 7). We examined the directional wave spectrum, D s p e c = f ( ω , θ ) of bottom skewness, where ω is the frequency in Hz and θ is the direction in degrees to determine whether negative skewness was related to spectral characteristics. The events associated with positive and negatively skewed waves ( S u ) in the directional wave spectrum were binned separately. We first averaged D s p e c along all directions to get mean frequency spectra ( D s p e c ω , Figure 8). The spectra associated with positive S u (green shaded region) frequencies are shifted toward lower frequencies compared to spectra associated with negative S u (red shaded region), which are shifted towards higher frequencies at all three sites (Figure 8a–c). Positively skewed waves also had narrower frequency distributions. The differences in mean frequency and standard deviation about the mean for spectra with positive versus negative skewness were significant with ( p < 0.01 ). We also found that the statistical skewness of the spectral distributions was higher for positively skewed waves (significantly different with p < 0.001 ). This pattern confirmed a shift of the spectral distribution towards lower frequencies associated with positively skewed waves.
Next, we averaged D s p e c across all frequencies to get spectra ( D s p e c θ ) varying over direction and binned based on the sign of skewness (Figure 9a–c). Positively skewed waves (green lines) were associated with a smaller directional spread of energy while negatively skewed waves (red) corresponded to a larger directional spread of energy (Figure 9a–c). This was confirmed by the kurtosis of the direction distributions, which was significantly different ( p < 0.001 ) with positively skewed waves exhibiting more peakiness (leptokurtic), and negatively skewed waves exhibiting broader directional spread, especially at MVCO. In summary, data from all three sites (Figure 8 and Figure 9) show that spectra associated with positively skewed waves were distributed over a lower peak frequency with a narrower directional spread. Conversely, spectra for negatively skewed waves were distributed over higher frequencies and with a larger directional spread of energy at all three sites. A similar pattern was noticed in data [12] from an 18-h storm event. They concluded that the sign of the skewed waves originated from wave-wave interactions, and the spectral shift towards higher frequencies originated from mixed sea and swell conditions, whereas unimodal swell conditions with lower frequencies led to positive skewness, consistent with Figure 8.

3.5. Estimating Wave-Driven Cross-Shore Bedload Using Skewed Waveform Parameters

We compared our three estimates of wave-driven cross-shore bedload transport (Section 2.4) using the Matanzas Inlet dataset because that site had the highest proportion of skewed waveforms (Section 3.4). The first two estimates were calculated using a slightly modified version of the SANTOSS formula ([38]; Appendix B). This formula requires estimates of wave asymmetry and calculates transport separately for each part of the wave cycle. We derived asymmetry values from surface waves (our “parameterized” method using the [34] approach; Appendix A) and from near-bottom measurements (our “direct” method; Appendix B) and compared it with a traditional excess stress-based bedload formulation (i.e., the MPM method detailed in Appendix C). Figure 10a,b compare instantaneous ( q b ) and time-integrated net ( 0 t q b   d t ) bedload transport rates where positive values in Figure 10a indicated bedload along the wave aligned axis and negative values are opposite to that. The bedload values reach local peaks during similar time intervals between all three methods and are associated with high wave energy events (refer to Figure 4) that result in larger skewness (Figure 10a). The parameterized method’s q b exceeded the direct one during 51% of time instances because it always predicted a positive skewness and as a result, the parameterized q b showed an increase of 25% in the ( 0 t q b   d t ) (Figure 10b). In both the direct and MPM method, there is a change in the sign of q b during certain periods although the magnitude of negative q b is much smaller over the entire time series for both methods. For all the events where q b > 0 , the ( 0 t q b   d t ) with MPM method was 700% smaller when compared to the direct skewed waveform based ( 0 t q b   d t ) .
Analysis of the time-stacked ripple images identified seven periods that were continuous, which allowed tracking of the ripple front. Manual interpretation (red lines, Figure 11a) was used to detect the starting and ending point of each continuous period. The ripple wavelength in these periods ranged from 0.065 to 45 cm. Using the distance traversed by ripple wave fronts over the duration of each segment, we found ripple migration rates ranging from 0.3 cm/day to 0.48 cm/day. The corresponding significant wave height ( H s ) and representative wave period ( T ) at the bottom ranged between 0.4 to 1.5 m and 5.8 to 14.4 s respectively (Figure 11b,c). The Ursell ( U r ) number during these periods ranged from 0.024 to 0.23 (Figure 11d). The relatively low values of these parameters (Figure 11b–d) for the periods when ripples could be tracked suggest the sonar was only capable of providing bedform transport estimates during weak events. During other times, ripples were not clearly visible in the images.
A comparison of instantaneous q b from sonar imagery and empirical methods (Figure 11e; note logarithmic y-axis) shows that there was substantial scatter within events, and that the four methods often overlapped. The q b   estimates are typically within 2 orders of magnitude with each other. The exception is the strongest event, between 06–09 March, when the sonar q b   estimates were lower than q b   calculated from both the asymmetric waveform methods and the MPM method. There are periods (02–03 February, 06–07 February, and 03–04 April) when the parameterized method results in zero q b (not visible on the logarithmic axis) because the calculated stress did not exceed the threshold critical stress in crest or trough cycles.
The RMSE of the cumulative bedload calculated from the asymmetric bedload method is lower than MPM for two events (1 and 7) but in all other events, it is greater (Table 5). For the asymmetric bedload, events 1 and 2 are significantly (positively) correlated ( r 1 ) with the sonar data explaining between 56–62% of variance. Events 4 and 5 also have a positive correlation but are not significant. For MPM, event 2 is the only one for which the explaining variance is above 50%.

4. Discussion

4.1. Skewness Parameterization for Operational Morphodynamic Models

In observational datasets, the near-bottom velocity skewness of the wave shape is often calculated with Equation (6) [53] and u 3 , can be used to estimate a net sediment transport rate driven by velocity-skewed oscillatory flows [4]. At the three study sites, parameterized S u is always positive while the direct skewness changes sign due to the factor of u 3 (also seen in Section 3.4) (Figure 12). Observations [12] showed the dependence of bedform migration direction on the sign of skewness where positively skewed waves resulted in an onshore migration of ripples and negatively skewed waves resulted in the offshore migration of ripples. However, in this study, the parameterized method that depends on Ruessink’s approach always predicts positively skewed waveform shapes, indicating onshore wave-driven sediment transport. Therefore, the application of this methodology in models may overpredict onshore wave-driven transport and require one to calibrate velocity skewness based on field data and/or modify the current and wave contributions [35].

4.2. Bedload Comparison from Various Methods

The large differences between the cumulative bedload (Figure 10b) from asymmetric bedload formulations and the MPM method may be attributed to several factors, all of which enhance transport in the direction of wave propagation. The asymmetric bedload formulation was designed to account for differential crest and trough cycle stresses and for processes associated with free surface wave propagation. These processes included the effect of modifying crest and trough periods due to surface waves that led to a longer crest period and shorter trough period experienced by sediments at the bed ([54], Equations (38) and (39) in Appendix B of this paper). Another important process is the increase in crest cycle stress and a decrease in trough cycle stress due to surface wave propagation (Equation (40)), also referred to as “wave streaming” or “boundary layer streaming” [38,55,56]. Using a high-resolution CFD model for waves and sediment that included a free surface [57] showed that progressive wave streaming produced a 60% increase in onshore sediment transport under nonbreaking wave conditions (only velocity skewed). Another mechanism that affects the asymmetric bedload calculation under surface waves is the adjustment of phase lag parameters under crest and trough to account for the differential horizontal sediment advection [38,58]. This mechanism leads to a decrease in effective settling velocity of sediments during crest (decreasing phase lag during crest) and an increase during trough (increasing phase lag during trough). However, the biggest difference between the parameterized method for wave shape used in conjunction with the non-linear transport formulation is the lack of negative skewness.
There was only limited sonar imagery suitable for extracting bedform transport rates for comparison with bedload calculations. The low pixel intensities and loss of bottom reflectors during storms due to plane-bed conditions or suspended material compromised most of the record. We explored the sensitivity of our migration-rate estimates to the starting point in the first two events (02–03 February and 06–07 February) in the time stack (Figure 11e). The resulting variations in migration rates varied as much as a factor of three, which is small compared to the order-of-magnitude changes in our bedform-transport estimates.
The choice of parameters (bed porosity factor, ripple shape factor, and ripple height to ripple wavelength mentioned in Section 2.5) would modify the results but not impact the order of magnitude changes in the sonar-based bedform transport. In future experiments, the ripple imagery could be improved by averaging over several sweeps to obtain clearer images of the ripples. This could reduce uncertainty in crest tracking and might also reduce the effort of manual interpretation and identification of the start and end points of migrating ripples in the time-stacked image. Despite the shortcomings, the dataset provides seven instances for comparison during weak to mild wave conditions and shows that transport rates from the asymmetric bedload formulations fall within an order of magnitude of rates from sonar data. This study is the first to compare and show that the asymmetric bedload formulations fall in the range of sonar data for all events (except event 5, likely due to the larger period waves). Overall, these results highlight the challenge in assessing the relative performance of bedload formulations and comparing it with field observations.

5. Conclusions

The use of skewed waveforms in sediment-transport calculations represents an advance in computing bedload transport that accounts for wave-driven processes in morphodynamic models. The present work compares parameterizations of skewed waveforms from bulk surface-wave statistics with waveforms obtained from direct measurements near the seafloor, using observations from three different sites under relatively calm conditions ( u W < 0.76 m/s). Overall, the correlations for all the waveform characteristics (velocity and period in each half cycle) were significant and remained above 0.85. The two methods produced skewed waveform characteristics that matched well at the three sites. The cumulative bedload calculated using the parameterized method was about 25% greater than that calculated using the direct method. These results demonstrate the efficacy of using the parameterized method to obtain representative waveforms for use in asymmetric bedload formulations. The direct measurements revealed instances of negative skewness at all three sites, but the equations used for the parameterized method cannot generate negative skewness. Episodes of negative skewness drive offshore transport, which explains why the parameterized approach predicts a larger onshore wave-driven bedload than the direct method. While this work shows the efficacy of Ruessink’s method in estimating asymmetric waveform characteristics, that approach relies on local measurements and does not account for shoaling history [29]. Other studies [59,60] have included non-local and offshore wave parameters to estimate wave asymmetry. They [60] studied the effects of offshore wave steepness, offshore spectral bandwidth, and beach slope and showed a 50% better prediction of asymmetric parameters nearing breaking regimes. Incorporating the effects of negative skewness in parameterizations and testing the effect of adding non-local wave conditions to estimate asymmetric waveforms could potentially improve the predictions of bedload. In addition, more observations at shallower depths where breaking waves lead to higher acceleration asymmetry could improve the empiricisms involved in the asymmetric bedload formulation, which is currently limited to non-breaking waves.

Author Contributions

Conceptualization, C.R.S. and J.C.W.; methodology, C.R.S., T.S.K. and S.E.S.; software, T.S.K., S.E.S., C.R.S. and G.R.L.; validation, T.S.K. and S.E.S.; Analysis, T.S.K., S.E.S., A.L.A. and C.R.S.; writing—original draft preparation, T.S.K. and S.E.S.; writing—review and editing, T.S.K., A.L.A., S.E.S., C.R.S. and J.C.W.; supervision, J.C.W. All authors have read and agreed to the published version of the manuscript.

Funding

USGS Coastal and Marine Hazards and Resources Program.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Table 1 specifies the filenames for accessing the data used in the above analysis. The datasets have been published on https://stellwagen.er.usgs.gov/, which is a Trusted Digital Repository for distributing time series data. These files can be downloaded from the following weblinks: 1. Fire Island [40,41]: https://stellwagen.er.usgs.gov/FIREISLAND14-a.html. 2. MVCO [42]: https://stellwagen.er.usgs.gov/mvco_15-a.html. 3. Matanzas Inlet [43]: https://stellwagen.er.usgs.gov/matanzas_18-a.html. All weblinks were last accessed on 29 January 2022.

Acknowledgments

We thank Kai Parker at the U.S. Geological Survey (USGS), St. Petersburg, Florida for providing his feedback to improve the quality of the paper. We thank Christie Hegermiller, Tom Hsu and Yashar Rafati for valuable discussions throughout the length of the project. We would like to acknowledge our colleague Jin-Si Over for her contribution in creating Figure 1 of the manuscript. We extend our gratitude to the team of scientists and engineers at USGS, Woods Hole Coastal, and Marine Science Center for gathering, processing, and publishing the observational datasets used in this work. The work by Tarandeep Kalra and Gibson Leavitt was done under contract to the USGS.

Conflicts of Interest

Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

Appendix A. Determining Skewed Waveform from Bulk Wave Statistics

The inputs required to determine the skewed waveform parameters from the wave spectrum include Hs (significant wave height), T (representative wave period), u w (amplitude of wave orbital velocity), and h (water depth). Using these inputs, one can calculate Ursell number ( U r ) [61] as:
U r = 3 a w k 4 ( k h ) 3  
where a w = 0.5 H s and k is the wave number computed from the dispersion relationship [39] that depends on T and h.
The third moments of skewness S u and asymmetry A u from 34,000 time series of natural near-bed orbital motion were computed according to [34]. They combined S u and A u into a non-dimensional skewed parameter B and found the phase angle ψ between skewness parameters depended on U r . B and ψ as a function of U r can be expressed as:
B = p 1 + p 2 p 1 1 + exp ( p 3 ln ( U r ) p 4 )
ψ = π 2 + ( π 2 ) t a n h ( p 5 U r p 6 )
where p 1 = 0.0 , p 2 = 0.85 , p 3 = 0.471 , p 4 = 0.297 , p 5 = 0.815 and p 6 = 0.672 .
S u and A u can be computed from B and ψ as follows:
S u = B c o s ( ψ )
A u = B s i n ( ψ )
The analytical expressions [3] for time series of velocity and acceleration by curve fitting using parameters r and .
u ( t ) = u w ( 1 r 2 ) [ s i n ( ω t ) + r s i n ( ) 1 + 1 r 2 ] [ 1 r c o s ( ω t + ) ]
a ( t ) = u w ω ( 1 r 2 ) [ c o s ( ω t ) r c o s ( ) r 2 1 + 1 r 2 s i n ( ) s i n ( ω t + ) ] [ 1 r c o s ( ω t + ) ] 2
where u w is the amplitude of wave orbital velocity and ω is the angular frequency expressed as:
ω = 2 π T b o t
r is an index of nonlinearity, is a waveform parameter whose values determined the profile of the waveform (sawtooth shaped, 1st-order cnoidal wave etc.) and t ranges from 0 t 2 π . To obtain r and , we utilized the skewed parameters, B and ψ obtained above from [34], written as:
r = 2 b R u e s s i n k b R u e s s i n k 2 + 1
where is defined as:
= φ π 2
where the parameter b R u e s s i n k is computed from B as:
b R u e s s i n k = 2 B B 2 + 9
A detailed discussion [3] on the implications of the choice of parameters r and to derive various nonlinear wave profiles. While the velocity in crest and trough along with skewness parameters can be obtained from analytical expressions [3], a more consistent way from a model implementation perspective was to use the formulations mentioned in [46]. From here onwards, the next set of equations used in step 1 are described in [46] that replace r by b (different from b R u e s s i n k ) and utilize with its sign reversed to get the time series for velocity and acceleration and are written as:
u ( t ) = f ( b , ) sin ( t ) b s i n ( ) 1 + b 2 2 b c o s ( t )
a ( t ) = f ( t , ) ( 1 + b 2 ) c o s ( t ) 2 c o s ( ) + 2 b 2 s i n s i n ( t ) [ 1 + b 2 2 b c o s ( t ) ] 2
where b can be written as:
b = r ( 1 + P )
where parameter P depends on r as follows:
P = 1 r 2
By knowing the maximum and minimum values of velocity, an analytic expression for velocity skewness ( R ) can be obtained from b and and is expressed as:
R = 1 2 ( 1 + b s i n ( ) )  
Similarly, acceleration skewness ( β ) is obtained from r and by using the following equations:
β = 0.5 + ( β r _ 0 0.5 )   s i n ( 0.5 π | | ) F 0 ) s i n ( 0.5 π F 0 )
where β r _ 0 is written as:
β r _ 0 = 0.5 ( 1 + r ) ,   i f   r 0.5  
β r _ 0 = 4 r ( 1 + r ) 4 r ( 1 + r ) + 1 , i f   r > 0.5
and where F 0 is written as:
F 0 = 1 0.27 ( 2 r ) ( 2.1 ) , i f   r 0.5
  F 0 = 0.59 + 0.14 ( 2 r ) ( 6.2 ) , i f   r 0.5
After finding velocity skewness ( R ) and acceleration asymmetry ( β ) parameters, the time of zero up-crossing and down-crossing of the waveform can be calculated from parameter c as:
t z u = s i n 1 ( c )
t z d = 2 c o s 1 ( c ) + t z u
where c = b s i n ( ) . The dimensional wave period in crest ( T c ) can be found from t z d and t z u (dimensionalized by dividing with angular frequency ω ).
T c = ( t z d t z u ) ω
T t = T b o t T c
The dimensional wave period in acceleration during crest and trough cycle is calculated based on t z u and t z u (Figure 2b) as:
T c u = ( t m t z u ) ω
T t u = ( t m t z d ) ω
where t m corresponds to the maximum or minimum phase in the expression defining time series of acceleration and can be written as:
t m = s i n 1 [ 4 c ( b 2 c 2 ) + ( 1 + b 2 ) ( 1 + b 2 2 c 2 ) ( 1 + b 2 ) 2 4 c 2 ]
Next, the non-dimensional peak crest and trough velocity are computed as:
u ^ c n o n d i m = ( 1 + c )
u ^ t n o n d i m = ( u ^ c n o n d i m 2 )
The crest and trough velocity are dimensionalized using the amplitude of wave orbital velocity ( u ^ ) and described as:
u c r e s t = u ^ c n o n d i m u w
u t r o u g h = u ^ t n o n d i m u w
By the end of the first step, one obtains a skewed waveform with the following parameters: velocity skewness ( R ) and acceleration asymmetry ( β ), accelerative/total wave period ( T c u ,   T c ,   T t u   and   T t ), and peak velocity ( u c r e s t   and   u t r o u g h ) in each half cycle.

Appendix B. Calculation of Wave-Driven Bedload

Calculation of wave-driven bedload due to asymmetric waveforms was done using the Santos formulations [38]. After obtaining the skewed waveform, calculations are performed to get non-dimensionalized shear stress (Shields parameter) for the full wave cycle (time-averaged absolute Shields parameter) and then the half wave cycle. For the full wave cycle, the ripple dimensions (height and length) [62] are determined from the maximum mobility number ( φ ^ ) that is defined as:
φ ^ = 1.27   u w 2 ( s 1 ) g d 50
where u w is the amplitude of orbital velocity and the denominator term s = ρ s / ρ , where ρ s and ρ are the densities of sediment and water, respectively, g is the gravitational constant, and d 50 is median grain size. The denominator term ( s 1 ) g d 50 appears repeatedly in the following text and is only defined here for the sake of brevity. Above a φ ^ of 240, the ripple regime transitions into a sheet flow regime. Once the bedform is determined and ripple dimensions are obtained, the time-averaged Shields parameter along with the current and wave friction factor for the full wave cycle are computed iteratively. The iterations need to be performed because the time-averaged Shields parameter depends on the friction factor calculation which in turn depends on the time-averaged Shields parameter. Note that for the current friction factor and near-bottom current velocity ( u d ) calculation, the edge of the wave boundary layer needs to be determined [63]. Since only the wave-driven bedload is considered in this work, u d and corresponding current friction factors are equal to zero.
Next, the effect of surface waves that modify the Lagrangian motion at the sediment bed by increasing the crest period and shortening the trough wave period is accounted for. The change in wave period is dependent on the wave propagation velocity ( c w ) and the horizontal grain displacement during each half cycle ( u ^ T π   ) and can be parameterized as:
Δ T c = ( ζ u ^ π c w 2 ζ u ^   )
T c _ n e w   = T c + Δ T c
Δ T t = ( ζ u ^ π c w + 2 ζ u ^   )
T t _ n e w   = T t Δ T t
where ζ = 0.55 [64], Δ T c and Δ T t are the change in crest and trough period that is added and subtracted in crest and trough wave periods to obtain the new half-cycle periods ( T c _ n e w ,   T t _ n e w ) . Based on the new half-cycle periods, the accelerative components are also modified proportionally to get T c u _ n e w   and T t u _ n e w   and can be written as:
T c u _ n e w   = ( T c u   T c _ n e w T c   )
T t u _ n e w   = ( T t u   T t _ n e w T t   )
After the wave period modification from the above two steps, the magnitude of crest and trough Shields parameter is computed. Each half cycle Shields parameter is then modified by the effects of progressive surface waves that account for the vertical gradient of horizontal momentum caused by the vertical orbital motion of surface waves. This phenomenon is parameterized through the contribution of a wave-averaged Reynolds stress term ( τ w R e ). The effect of τ w R e is to increase the crest cycle stress (direction of wave propagation) and decrease the trough cycle stress. τ w R e is expressed as:
τ w R e = ρ f w d   α w ( u ^ ) ( 3 ) 2 c w
where α w = 0.424, f w d is the full-cycle wave-current friction factor and c w is the wave speed.
Following the calculation of Shields parameter in each half cycle, the next step is to obtain a phase lag parameter that determines the distribution of sediment load between each half cycle. The phase lag parameter is comprised of two terms that account for the ratio of sediment stirring height ( l b e d f o r m 2 ( T i T i u ) w s i ) and effects of progressive surface waves that induce horizontal sediment flux gradients ( 1 ε u ^ i c w ) [58] in each half cycle. The phase lag parameter in each half cycle can be written as:
P i = α ( 1 ε u ^ i c w ) ( l b e d f o r m 2 ( T i T i u ) w s i )  
where α and ε are set to 8.2 and 1.7, u ^ i is the peak half cycle velocity, c w is the wave, l b e d f o r m corresponds to ripple height or sheet flow thickness, T i is the duration in a particular half cycle while T i u is the duration of accelerating flow in that cycle and w s i is the sediment settling velocity. w s i is calculated using the settling velocity for still water [39] and modified due to the vertical orbital velocity in each half-cycle caused by progressive surface waves [38,58].
If the phase lag parameter in each half cycle is less than one, then all the sediment entrained in each half cycle is transported within that half cycle. In that case, the bedload transport mechanism is only controlled by the Shields parameter calculated in that half cycle. If the phase lag parameter is greater than one, then the sediment entrained in one half cycle can be transported to the other half cycle. The resulting volumetric bedload transport per unit width ( q b ) is given by:
q b = ( | θ c |   T c   ( Ω c c   + T c   2 T c   Ω t c ) θ c | θ c | + | θ t |   T t ( Ω t t   + T t   2 T t u   Ω c t   ) θ t | θ t | T ) ( s 1 ) g d 50 3  
where s = ρ s / ρ , ρ s and ρ are the densities of sediment and water respectively, d 50 is median grain size, θ is the non-dimensional bed stress, (Shields parameter), θ c and θ t are the Shields parameter in the crest and trough, respectively and T c   and T t   indicate crest and trough cycle wave periods, respectively. T c u   and T t u refer to the wave period of accelerating flow within crest and trough cycles. Ω c c is the sediment bedload that is entrained during the wave crest and transported during the crest period, Ω t c   is the sediment bedload that is entrained during the wave crest and transported through the trough period, Ω t t   is the sediment bedload that is entrained during the wave trough and transported during the trough period, Ω c t is the sediment bedload that is entrained during the wave trough and transported during the crest period.

Appendix C. Calculation of Meyer-Peter and Müller (MPM) Bedload

In MPM, bedload transport depends on excessive shear stress, above a threshold for motion, acting on a uniform bed of grains, and the bedload is defined as:
q b = {   γ ( θ θ c r ) 1.5 ( s 1 ) g d 50 3 ,   i f ,   θ > θ c r   θ = 0.0 ,   i f ,   θ   θ c r   }
where γ = 8 is a constant of proportionality (typically ranging between 8 and 12) and θ c r is the critical Shields parameter that can be calculated for a given grain size [39]. The wave skin friction Shields parameter θ is given by
θ = 1 2 f w u ˜ | u ˜ | g ( s 1 ) d 50
where u ˜ is the dominant wave signal (Section 2.2) and f w is the wave friction factor [65] and is written as:
f w = exp [ 5.5 ( 2.5 d 50 ω b r u ˜ b r ) 0.2 6.3 ]
where u ˜ b r is the representative near-bottom wave orbital velocity for each wave burst is calculated [44] and is written as:
u b r = 2 i S u v d f  
where, S u v is the Welch’s power spectral density for the near-bed horizontal velocity components,   d f corresponds to the frequency bins in the spectrum and subscript i corresponds to the individual wave bursts and ω b r is the angular frequency ( ω ) obtained from the spectral amplitudes of wave velocities ( u i ) and is expressed as:
ω b r =   ω i u i   u i
Note that the velocity signal undergoes a filtration process to only include incident wave band as described in Section 2.2.

References

  1. Svendsen, I.A. Introduction to Nearshore Hydrodynamics; World Scientific: Singapore, 2006; Volume 24, p. 774. [Google Scholar]
  2. Elfrink, B.; Hanes, D.M.; Ruessink, B.G. Parameterization and simulation of near bed orbital velocities under irregular waves in shallow water. Coast. Eng. 2006, 53, 915–927. [Google Scholar] [CrossRef]
  3. Abreu, T.; Silva, P.A.; Sancho, F.; Temperville, A. Analytical approximate wave form for non-linear waves. Coast. Eng. 2010, 57, 656–667. [Google Scholar] [CrossRef]
  4. Ribberink, J.S.; Al-Salem, A.A. Sediment transport in oscillatory boundary layers in cases of rippled beds and sheet flow. J. Geophys. Res. Ocean. 1994, 99, 12707–12727. [Google Scholar] [CrossRef]
  5. Watanabe, A.; Sato, S. A sheet-flow transport rate formula for asymmetric, forward-leaning waves and currents. Coast. Eng. 2004, 4, 1703–1714. [Google Scholar] [CrossRef]
  6. Silva, P.A.; Temperville, A.; Seabra Santos, F. Sand transport under combined current and wave conditions: A semi-unsteady, practical model. Coast. Eng. 2006, 53, 897–913. [Google Scholar] [CrossRef] [Green Version]
  7. Hsu, T.J.; Elgar, S.; Guza, R.T. Wave-induced sediment transport and onshore sandbar migration. Coast. Eng. 2006, 53, 817–824. [Google Scholar] [CrossRef] [Green Version]
  8. Ruessink, B.G.; Kuriyama, Y.; Reniers, A.J.H.M.; Roelvink, J.A.; Walstra, J.A. Modeling cross-shore sandbar behavior on the timescales of weeks. J. Geophys. Res. Earth Surf. 2007, 112. [Google Scholar] [CrossRef] [Green Version]
  9. Fernández-Mora, A.; Calvete, D.; Falqués, A.; De Swart, H.E. Onshore sandbar migration in the surf zone: New insights into the wave-induced sediment transport mechanisms. Geophys. Res. Lett. 2015, 42, 2869–2877. [Google Scholar] [CrossRef] [Green Version]
  10. Dong, L.P.; Sato, S.; Liu, H. A sheetflow sediment transport model for skewed-non-linear waves combined with strong opposite currents. Coast. Eng. 2013, 71, 87–101. [Google Scholar] [CrossRef]
  11. Traykovski, P.; Hay, A.E.; Irish, J.D.; Lynch, J.F. Geometry, migration, and evolution of wave orbital ripples at LEO. J. Geophys. Res. Ocean. 1999, 104, 1505–1524. [Google Scholar] [CrossRef]
  12. Crawford, A.M.; Hay, A.E. Linear transition ripple migration and wave orbital velocity skewness: Observations. J. Geophys. Res. 2001, 106, 113–128. [Google Scholar] [CrossRef]
  13. Hurther, D.; Thorne, P.D. Suspension and near bed load sediment transport processes above a migrating, sand rippled bed under shoaling waves. J. Geophys. Res. 2011, 116, C07001. [Google Scholar] [CrossRef] [Green Version]
  14. Drake, G.; Calantoni, J. Discrete particle model for sheet flow sediment transport in the nearshore. J. Geophys. Res. 2001, 106, 19859–19868. [Google Scholar] [CrossRef]
  15. Hoefel, F.; Elgar, S. Wave-induced sediment transport and sandbar migration. Science 2003, 299, 1885–1887. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Nielsen, P. Sheet flow sediment transport under waves with acceleration skewness and boundary layer streaming. Coast. Eng. 2006, 53, 749–758. [Google Scholar] [CrossRef]
  17. Ranasinghe, R.; Pattiaratchi, C.; Masselink, G. Morphodynamic model to simulate the seasonal closure of tidal inlets. Coast. Eng. 1999, 37, 1–36. [Google Scholar] [CrossRef]
  18. van der Zanden, J.; van der A, D.A.; Hurther, D.; Cáceres, I.; O’Donoghue, T.; Hulscher, S.J.M.H.; Ribberink, J.S. Bedload and suspended load contributions to breaker bar morphodynamics. Coast. Eng. 2017, 129, 74–92. [Google Scholar] [CrossRef]
  19. Mieras, R.S.; Puleo, J.A.; Anderson, D.; Hsu, T.-J.; Cox, D.T.; Calantoni, J. Relative contributions of bed load and suspended load to sediment transport under skewed-asymmetric waves on a sandbar crest. J. Geophys. Res. Oceans 2019, 124, 1294–1321. [Google Scholar] [CrossRef]
  20. Ridderinkhof, W.; de Swart, H.E.; van der Vegt, M.; Hoekstra, P. Modeling the growth and migration of sandy shoals on ebb-tidal deltas. J. Geophys. Res. Earth Surf. 2016, 121, 1351–1372. [Google Scholar] [CrossRef] [Green Version]
  21. Reniers, A.J.H.M.; De Wit, F.P.; Tissier, M.F.S.; Pearson, S.G.; Brakenhoff, L.B.; Vegt, V.D.; Mol, J.; Van Prooien, B.C. Wave-skewness and current-related ebb-tidal sediment transport: Observations and modeling coastal sediments. In Proceedings of the 9th International Conference on Coastal Sediments, St. Petersburg, FL, USA, 27–31 May 2019; pp. 2018–2028. [Google Scholar] [CrossRef]
  22. Lenstra, K.J.H.; Pluis, S.R.P.M.; Ridderinkhof, W.; Ruessink, G.; Van der Vegt, M. Cyclic channel-shoal dynamics at the Ameland inlet: The impact on waves, tides, and sediment transport. Ocean. Dyn. 2019, 69, 409–425. [Google Scholar] [CrossRef] [Green Version]
  23. De Wit, F.; Tissier, M.; Reniers, A. Characterizing wave shape evolution on an ebb-tidal shoal. J. Mar. Sci. Eng. 2019, 7, 367. [Google Scholar] [CrossRef] [Green Version]
  24. Gazi, A.H.; Purkayastha, S.; Afzal, M.S. The equilibrium scour depth around a pier under the action of collinear waves and current. J. Mar. Sci. Eng. 2020, 8, 36. [Google Scholar] [CrossRef] [Green Version]
  25. Gazi, A.H.; Saud, M.S. A new mathematical model to calculate the equilibrium scour depth around a pier. Acta Geophys. 2020, 68, 181–187. [Google Scholar] [CrossRef]
  26. Zhang, G.; Chen, Y.; Cheng, W.; Zhang, H.; Gong, W. Wave effects on sediment transport and entrapment in a channel-shoal estuary: The Pearl River estuary in the dry winter season. J. Geophys. Res. Ocean. 2021, 126. [Google Scholar] [CrossRef]
  27. Roelvink, D.; Reniers, A.; Van Dongeren, A.; Van Thiel de Vries, J.; McCall, R.; Lescinski, J. Modelling storm impacts on beaches, dunes and barrier islands. Coast. Eng. 2009, 56, 1133–1152. [Google Scholar] [CrossRef]
  28. Warner, J.C.; Armstrong, B.; He, R.; Zambon, J.B. Development of a coupled ocean-atmosphere-wave-sediment transport (COAWST) modeling system. Ocean. Model. 2010, 35, 230–244. [Google Scholar] [CrossRef] [Green Version]
  29. Sherwood, C.R.; Dongeren, A.V.; Doyle, J.; Hegermiller, C.A.; Hsu, T.J.; Kalra, T.S.; Olabarrieta, M.; Penko, A.M.; Rafati, Y.; Roelvink, D.; et al. Modeling the morphodynamics of coastal responses to extreme events: What shape are we in? Annu. Rev. Mar. Sci. 2022, 14, 457–492. [Google Scholar] [CrossRef]
  30. Zijlema, M.; Stelling, G.; Smit, P. SWASH: An operational public domain code for simulating wave fields and rapidly varied flows in coastal waters. Coast. Eng. 2011, 58, 992–1012. [Google Scholar] [CrossRef]
  31. Malej, M.; Mith, J.M.; Salgado-Dominguez, G. Introduction to Phase-Resolving Wave Modeling with FUNWAVE; ERDC/CHL CHETN-I-87; US Army Corps of Engineers: Washington, DC, USA, 2015. [Google Scholar]
  32. Doering, J.C.; Elfrink, B.; Daniel, M.; Hanes, D.M.; Ruessink, G. Parameterization of velocity skewness under waves and its effect on cross-shore sediment transport. Coast. Eng. 2000, 276, 1383–1397. [Google Scholar] [CrossRef]
  33. Isobe, M.; Horikawa, K. Study on water particle velocities of shoaling and breaking waves. Coast. Eng. 1982, 25, 109–123. [Google Scholar] [CrossRef]
  34. Ruessink, B.G.; Ramaekers, G.; Van Rijn, L.C. On the parameterization of the free-stream non-linear wave orbital motion in nearshore morphodynamic models. Coast. Eng. 2012, 65, 56–63. [Google Scholar] [CrossRef]
  35. Albernaz, M.B.; Ruessink, G.; Bert Jagers, H.R.A.; Kleinhans, M.G. Effects of wave orbital velocity parameterization on nearshore sediment transport and decadal morphodynamics. J. Mar. Sci. Eng. 2019, 7, 188. [Google Scholar] [CrossRef] [Green Version]
  36. Chen, W.L.; Dodd, N. An idealised study for the evolution of a shoreface nourishment. Cont. Shelf Res. 2019, 178, 15–26. [Google Scholar] [CrossRef]
  37. Rafati, Y.; Hsu, T.-J.; Elgar, S.; Raubenheimer, B.; Quataert, E.; van Dongeren, A. Modeling the hydrodynamics and morphodynamics of sandbar migration events. Coast. Eng. 2021, 166, 103885. [Google Scholar] [CrossRef]
  38. 38. van der, A.D.A.; Ribberink, J.S.; van der Werf, J.J.; O’Donoghue, T.; Buijsrogge, R.H.; Kranenburg, W.M. Practical sand transport formula for non-breaking waves and currents. Coast. Eng. 2013, 76, 26–42. [Google Scholar] [CrossRef] [Green Version]
  39. Soulsby, R.L. Dynamics of Marine Sands; Thomas Telford Publications: London, UK, 1997. [Google Scholar]
  40. Armstrong, B.N.; Warner, J.C.; List, J.H.; Martini, M.A.; Montgomery, E.T.; Traykovski, P.; Voulgaris, G. Coastal Change Processes Project Data Report for Oceanographic Observations near Fire Island, New York, February through May 2014: U.S. Geological Survey Open-File Report; USGS: Reston, VA, USA, 2015; Volume 1033. [CrossRef]
  41. Armstrong, B.N.; Warner, J.C.; List, J.H.; Martini, M.A.; Montgomery, E.T. Oceanographic Measurements—Fire Island, NY, Nearshore, 2014: U.S. Geological Survey Data Release; USGS: Reston, VA, USA, 2015. [CrossRef]
  42. Martini, M.A.; Montgomery, E.T.; Sherwood, C.R. Oceanographic and Water-Quality Measurements Collected South of Martha’s Vineyard, MA, November–December 2015: U.S. Geological Survey Data Release; USGS: Reston, VA, USA, 2017. [CrossRef]
  43. Suttles, S.E.; Warner, J.C.; Montgomery, E.T.; Martini, M.A. Oceanographic and Water Quality Measurements in the Nearshore Zone at Matanzas Inlet, Florida, January–April, 2018: U.S. Geological Survey Data Release; USGS: Reston, VA, USA, 2018. [CrossRef]
  44. Wiberg, P.L.; Sherwood, C.R. Calculating wave-generated bottom orbital velocities from surface-wave parameters. Comput. Geosci. 2008, 34, 1243–1262. [Google Scholar] [CrossRef]
  45. Silva, P.A.; Abreu, T.; Dominic, V.D.A.A.; Sancho, F.; Ruessink, B.G.; Van der Werf, J.J.; Ribberink, J.S. Sediment transport in nonlinear skewed oscillatory flows: Transkew experiments. J. Hydraul. Res. 2011, 49, 72–80. [Google Scholar] [CrossRef] [Green Version]
  46. Malarkey, J.; Davies, A.G. Free-stream velocity descriptions under waves with skewness and asymmetry. Coast. Eng. 2012, 68, 78–95. [Google Scholar] [CrossRef]
  47. Voulgaris, G.; Morin, J.P. A long-term real time sea bed morphology evolution system in the South Atlantic Bight. In Proceedings of the IEEE/OES 9th Working Conference on Current Measurement Technology, Charleston, SC, USA, 17–19 March 2008; pp. 71–79. [Google Scholar] [CrossRef]
  48. Cheel, R.A.; Hay, A.E. Cross-ripple patterns and wave directional spectra. J. Geophys. Res. 2008, 113, C10009. [Google Scholar] [CrossRef] [Green Version]
  49. Englert, C.M. Development of a Spectral 2-D Fast Fourier Transform Analysis for Sand Ripple Morphology Interpretation. Honors Thesis, Paper 583. Colby College, Waterville, ME, USA, 2010. [Google Scholar]
  50. Nelson, T.R.; Voulgaris, G. Temporal and spatial evolution of wave-induced ripple geometry: Regular versus irregular ripples. J. Geophys. Res. Ocean. 2014, 119, 664–688. [Google Scholar] [CrossRef]
  51. Canny, J. A computational approach to edge detection. IEEE Trans. Pattern Anal. Mach. Intell. 1986, 8, 679–698. [Google Scholar] [CrossRef]
  52. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  53. Elgar, S. Relationships involving third moments and bispectra of a harmonic process. IEEE Trans. Acoust. Speech Signal Process. 1987, 35, 1725–1726. [Google Scholar] [CrossRef]
  54. Ribberink, J.S.; Van der, A.D.A.; Buijsrogge, R.H. Santoss Transport Model: A New Formula for Sand Transport under Waves and Currents; Report number: SANTOSS_UT_IR3; University of Twente: Enschede, The Netherlands; University of Aberdeen: Aberdeen, UK, 2010. [Google Scholar]
  55. Holmedal, L.E.; Myrhaug, D.; Håvard, R. The sea bed boundary layer under random waves plus current. Cont. Shelf Res. 2003, 27, 717–750. [Google Scholar] [CrossRef]
  56. Afzal, M.S.; Holmedal, L.E.; Myrhaug, M. Sediment transport in combined wave–current seabed boundary layers due to streaming. J. Hydraul. Eng. 2021, 147. [Google Scholar] [CrossRef]
  57. Kim, Y.; Cheng, Z.; Hsu, T.-J.; Chauchat, J. A numerical study of sheet flow under monochromatic nonbreaking waves using a free surface resolving Eulerian two-phase flow model. J. Geophys. Res. Ocean. 2018, 123, 4693–4719. [Google Scholar] [CrossRef]
  58. Kranenburg, W.M.; Ribberink, J.S.; Schretlen, J.J.L.M.; Uittenbogaard, R.E. Sand transport beneath waves: The role of progressive wave streaming and other free surface effects. J. Geophys. Res. Earth Surf. 2013, 118, 122–139. [Google Scholar] [CrossRef] [Green Version]
  59. Cienfuegos, R.; Barthelemy, E.; Bonneton, P.; Gondran, X. Nonlinear surf zone wave properties as estimated from boussinesq modelling: Random waves and complex bathymetries. Coast. Eng. 2006, 5, 360–371. [Google Scholar] [CrossRef] [Green Version]
  60. Rocha, M.V.L.; Michallet, H.; Silva, P.A. Improving the parameterization of wave nonlinearities—The importance of wave steepness, spectral bandwidth and beach slope. Coast. Eng. 2017, 121, 77–89. [Google Scholar] [CrossRef]
  61. Doering, J.C.; Bowen, A.J. Parametrization of orbital velocity asymmetries of shoaling and breaking waves using bispectral analysis. Coast. Eng. 1995, 26, 15–33. [Google Scholar] [CrossRef]
  62. O’Donoghue, T.; Doucette, J.S.; Van der Werf, J.J.; Ribberink, J.S. The dimensions of sand ripples in full-scale oscillatory flows. Coast. Eng. 2006, 53, 997–1012. [Google Scholar] [CrossRef]
  63. Kalra, T.S.; Sherwood, C.R.; Warner, J.C.; Rafati, Y.; Hsu, T.J. Investigating bedload transport under asymmetrical waves using a coupled ocean-wave model. In Proceedings of the 9th International Conference on Coastal Sediments, St. Petersburg, FL, USA, 27–31 May 2019; World Scientific: Singapore; pp. 591–604. [Google Scholar] [CrossRef]
  64. Schretlen, J.J.L.M. Ph.D. Thesis, University of Twente, Enschede, The Netherlands, 2010.
  65. Nielsen, P. Coastal Bottom Boundary Layers and Sediment Transport; World Scientific: Singapore, 1992. [Google Scholar]
Figure 1. Study sites: (a) Fire Island, New York, (b) Martha’s Vineyard, Massachusetts, (c) Matanzas Inlet, Florida.
Figure 1. Study sites: (a) Fire Island, New York, (b) Martha’s Vineyard, Massachusetts, (c) Matanzas Inlet, Florida.
Jmse 10 00223 g001
Figure 2. (a) Raw ( u r a w ) and filtered incident wave band ( u ˜ i n c i d e n t ) spectra for a representative burst time series during the 14 February 2014, storm event in the Fire Island dataset, (b) individual waveforms (dashed lines) and the representative waveform (solid line) for the burst time series.
Figure 2. (a) Raw ( u r a w ) and filtered incident wave band ( u ˜ i n c i d e n t ) spectra for a representative burst time series during the 14 February 2014, storm event in the Fire Island dataset, (b) individual waveforms (dashed lines) and the representative waveform (solid line) for the burst time series.
Jmse 10 00223 g002aJmse 10 00223 g002b
Figure 3. (a) Image from the rotary fan-beam corresponding to the first image captured on 3 February 2018 at Matanzas Inlet, Florida. The transect chosen for creating the time stacked image is indicated in bold yellow line and the square sub-patch is the region chosen for ripple wavelength calculation. (b) Time stacked image to calculate the ripple migration rate from the transect.
Figure 3. (a) Image from the rotary fan-beam corresponding to the first image captured on 3 February 2018 at Matanzas Inlet, Florida. The transect chosen for creating the time stacked image is indicated in bold yellow line and the square sub-patch is the region chosen for ripple wavelength calculation. (b) Time stacked image to calculate the ripple migration rate from the transect.
Jmse 10 00223 g003aJmse 10 00223 g003b
Figure 4. Datasets of significant wave height and peak wave period estimated from the upward-looking acoustic Doppler current profile respectively from (a) and (b) Fire Island, (c) and (d) MVCO, (e) and (f) Matanzas Inlet.
Figure 4. Datasets of significant wave height and peak wave period estimated from the upward-looking acoustic Doppler current profile respectively from (a) and (b) Fire Island, (c) and (d) MVCO, (e) and (f) Matanzas Inlet.
Jmse 10 00223 g004
Figure 5. Comparison of representative wave orbital velocity, uw (m/s), and representative wave period, T (s) from ADV based near-bottom (direct) and surface wave spectrum (parameterized) measurements at the Fire Island site (a,b), MVCO site (c,d), and Matanzas Inlet site (e,f), respectively. The color of the circles indicates significant wave height in m. The dashed and solid black lines correspond to an identity line and a linear fit to the points, respectively.
Figure 5. Comparison of representative wave orbital velocity, uw (m/s), and representative wave period, T (s) from ADV based near-bottom (direct) and surface wave spectrum (parameterized) measurements at the Fire Island site (a,b), MVCO site (c,d), and Matanzas Inlet site (e,f), respectively. The color of the circles indicates significant wave height in m. The dashed and solid black lines correspond to an identity line and a linear fit to the points, respectively.
Jmse 10 00223 g005
Figure 6. Scatter plot showing a comparison of ADV based direct measurements and surface wave spectrum based parameterized measurements (a) peak crest cycle velocity, u c r e s t (m/s); (b) peak trough cycle velocity, u t r o u g h   (m/s); (c) crest cycle period, T c (s); (d) trough cycle period T t (s); (e) crest cycle accelerative period ( T c u ) , (s); and (f) trough cycle accelerative period, T t u (s). The dashed black line shows an identity line and solid-colored lines correspond to site-specific linear fit to data points.
Figure 6. Scatter plot showing a comparison of ADV based direct measurements and surface wave spectrum based parameterized measurements (a) peak crest cycle velocity, u c r e s t (m/s); (b) peak trough cycle velocity, u t r o u g h   (m/s); (c) crest cycle period, T c (s); (d) trough cycle period T t (s); (e) crest cycle accelerative period ( T c u ) , (s); and (f) trough cycle accelerative period, T t u (s). The dashed black line shows an identity line and solid-colored lines correspond to site-specific linear fit to data points.
Jmse 10 00223 g006
Figure 7. Velocity skewness from ADV data based on representative waveforms obtained from the direct method (Section 2.3) at all the sites as a function of Ursell number. The black dashed line indicates a threshold value of U r = 0.04 .
Figure 7. Velocity skewness from ADV data based on representative waveforms obtained from the direct method (Section 2.3) at all the sites as a function of Ursell number. The black dashed line indicates a threshold value of U r = 0.04 .
Jmse 10 00223 g007
Figure 8. Directional wave spectrum averaged over all directions ( D s p e c ω ) for the three sites (a) Fire Island site, (b) MVCO site, and (c) Matanzas Inlet site. Mean values are indicated by bold lines and shaded regions highlight the standard deviation of wave spectrum. Green bold lines and shaded region bin the events associated with positive wave-orbital velocity skewness. Red bolds line and shaded region bin the events associated with negative skewness.
Figure 8. Directional wave spectrum averaged over all directions ( D s p e c ω ) for the three sites (a) Fire Island site, (b) MVCO site, and (c) Matanzas Inlet site. Mean values are indicated by bold lines and shaded regions highlight the standard deviation of wave spectrum. Green bold lines and shaded region bin the events associated with positive wave-orbital velocity skewness. Red bolds line and shaded region bin the events associated with negative skewness.
Jmse 10 00223 g008
Figure 9. Directional wave spectrum averaged over all frequencies ( D s p e c θ ) for the three sites (a) Fire Island site, (b) MVCO site, and (c) Matanzas Inlet site. Mean values are indicated by bold lines and shaded regions highlight the standard deviation of wave spectrum. Green bold lines and shaded region bin the events associated with positive wave-orbital velocity skewness. Red bold lines and shaded region bin the events associated with negative skewness.
Figure 9. Directional wave spectrum averaged over all frequencies ( D s p e c θ ) for the three sites (a) Fire Island site, (b) MVCO site, and (c) Matanzas Inlet site. Mean values are indicated by bold lines and shaded regions highlight the standard deviation of wave spectrum. Green bold lines and shaded region bin the events associated with positive wave-orbital velocity skewness. Red bold lines and shaded region bin the events associated with negative skewness.
Jmse 10 00223 g009
Figure 10. Comparison of (a) instantaneous, ( q b   ) and (b) cumulative ( 0 t q b   d t ) wave-driven bedload transport using parameterized and direct asymmetric waveform and MPM methods. Positive rates in Figure 10a indicated bedload along the wave-aligned axis and negative values are opposite to that.
Figure 10. Comparison of (a) instantaneous, ( q b   ) and (b) cumulative ( 0 t q b   d t ) wave-driven bedload transport using parameterized and direct asymmetric waveform and MPM methods. Positive rates in Figure 10a indicated bedload along the wave-aligned axis and negative values are opposite to that.
Jmse 10 00223 g010
Figure 11. (a) Time stacked sonar imagery at the Matanzas Inlet site with periods of ripple tracking shown in red, (b) significant wave height, (c) representative wave period with red markings corresponding to the periods of ripple movement, (d) Ursell number with red markings corresponding to the periods of ripple movement, and (e) instantaneous bedload on a log scale in mm2/s comparing the bedload with parameterized (black) and direct (blue) waveforms using the asymmetric bedload formulation (Appendix B), MPM (green) method (Appendix C) and sonar imagery (brown). Event numbers are marked corresponding to the seven ripple migration events tracked in the sonar data (Figure 10e).
Figure 11. (a) Time stacked sonar imagery at the Matanzas Inlet site with periods of ripple tracking shown in red, (b) significant wave height, (c) representative wave period with red markings corresponding to the periods of ripple movement, (d) Ursell number with red markings corresponding to the periods of ripple movement, and (e) instantaneous bedload on a log scale in mm2/s comparing the bedload with parameterized (black) and direct (blue) waveforms using the asymmetric bedload formulation (Appendix B), MPM (green) method (Appendix C) and sonar imagery (brown). Event numbers are marked corresponding to the seven ripple migration events tracked in the sonar data (Figure 10e).
Jmse 10 00223 g011
Figure 12. Near-bottom skewness velocity from direct waveforms from the three sites calculated using Equation (6) compared with parameterized skewness using Equation (A4). Colors indicate significant wave height in m.
Figure 12. Near-bottom skewness velocity from direct waveforms from the three sites calculated using Equation (6) compared with parameterized skewness using Equation (A4). Colors indicate significant wave height in m.
Jmse 10 00223 g012
Table 1. Filenames containing the datasets used to obtain near bottom and surface wave measurements.
Table 1. Filenames containing the datasets used to obtain near bottom and surface wave measurements.
Wave ParameterFire Island, NY
(40.6193, −73.1840)
MVCO, MA
(41.3337, −70.5659)
Matanzas Inlet, FL
(29.7114, −81.2186)
ADV burst9917advb-cal.nc10577vecb-cal.nc11109vecb-a.nc
ADV burst statistics9917advs-cal.nc10577vecs-a.nc11109vecs-a.nc
Workhorse9921whp-cal.nc10571whVp-cal.nc11101whVp-cal.nc and
Seabird Seaguage---
Sonar--11107hffan_raw.cdf
ADV burst9917advb-cal.nc10577vecb-cal.nc11109vecb-a.nc
Table 2. ADCP sensor specifications for the three sites. The sensor height for the velocity measurements changed during the deployment period but the median elevation remained around the nominal elevation (latitude and longitude in decimal degrees).
Table 2. ADCP sensor specifications for the three sites. The sensor height for the velocity measurements changed during the deployment period but the median elevation remained around the nominal elevation (latitude and longitude in decimal degrees).
Site LocationFire Island, NY
(40.6193, −73.1840)
MVCO, MA
(41.3337, −70.5659)
Matanzas Inlet, FL
(29.7114, −81.2186)
Collection period7 February–5 May 2014 12 November 15–December 2015 24 January–13 April 2018
Sensor typeRD Instruments ADCPTRDI VTRDI V
Instrument frequency600 kHz1000 kHz1000 kHz
Burst sampling rate2 Hz2 Hz2 Hz
Burst sampling length1024 s every hour1024 s every hour2048 s every hour
Initial instrument elevation above bottom2.1 m2.4 m2.4 m
Table 3. ADV sensor specifications for the three sites. The sensor height for the velocity measurements changed during the deployment period (distance to the measurement volume).
Table 3. ADV sensor specifications for the three sites. The sensor height for the velocity measurements changed during the deployment period (distance to the measurement volume).
Site LocationFire Island, NYMVCO, MAMatanzas Inlet, FL
Collection period7 February–3 May 2014 17 November–8 December 2015 24 January–13 April 2018
Sensor typeSontek ADVNortek Vector ADVNortek Vector ADV
Instrument acoustic frequency
Burst sampling rate
5000 kHz
8 Hz
6000 kHz
16 Hz
6000 kHz
16 Hz
Burst sampling length1050 s every hour1875 s every hour2048 s every hour
Measurement location (height above bed)0–35 cm64 cm20–40 cm
Pressure sensor height1.14–1.5 m1.69 m1.97 m
Table 4. Errors between direct and parameterized waveform characteristics at all three sites where E1 = R M S E i and E2 = 100 ( R M S E i ) / R M S i d where subscript “i” is for an individual waveform characteristic and subscript “d” refers to the direct method.
Table 4. Errors between direct and parameterized waveform characteristics at all three sites where E1 = R M S E i and E2 = 100 ( R M S E i ) / R M S i d where subscript “i” is for an individual waveform characteristic and subscript “d” refers to the direct method.
Wave ParameterFire Island, NYMVCO, MAMatanzas Inlet, FL
E1E2E1E2E1E2
uw (m/s)0.0419.630.0420.620.0515.89
T (s)0.919.500.788.441.3413.88
ucrest (m/s)0.0419.020.0420.170.0414.76
utrough (m/s)0.0420.280.0421.130.0517.48
Tc (s)0.469.640.408.710.6914.61
Tt (s)0.449.250.388.230.6513.28
Tcu (s)0.3513.890.3212.300.3815.94
Ttu (s)0.219.090.177.960.3414.17
Table 5. Event-based cumulative error (RMSE) and correlation coefficient ( r 1 ) between the calculated transports (both asymmetric direct and MPM) and the sonar estimates. Significant correlations (p < 0.05) are indicated in bold.
Table 5. Event-based cumulative error (RMSE) and correlation coefficient ( r 1 ) between the calculated transports (both asymmetric direct and MPM) and the sonar estimates. Significant correlations (p < 0.05) are indicated in bold.
Event PeriodAsymmetric (Direct)MPM
RMSE r 1 RMSE r 1
10.040.790.11−0.72
20.340.750.060.76
30.42−0.020.22−0.58
41.520.240.580.13
557.90.269.490.36
60.22−0.50.110.2
70.2−0.60.26−0.45
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kalra, T.S.; Suttles, S.E.; Sherwood, C.R.; Warner, J.C.; Aretxabaleta, A.L.; Leavitt, G.R. Shoaling Wave Shape Estimates from Field Observations and Derived Bedload Sediment Rates. J. Mar. Sci. Eng. 2022, 10, 223. https://doi.org/10.3390/jmse10020223

AMA Style

Kalra TS, Suttles SE, Sherwood CR, Warner JC, Aretxabaleta AL, Leavitt GR. Shoaling Wave Shape Estimates from Field Observations and Derived Bedload Sediment Rates. Journal of Marine Science and Engineering. 2022; 10(2):223. https://doi.org/10.3390/jmse10020223

Chicago/Turabian Style

Kalra, Tarandeep S., Steve E. Suttles, Christopher R. Sherwood, John C. Warner, Alfredo L. Aretxabaleta, and Gibson R. Leavitt. 2022. "Shoaling Wave Shape Estimates from Field Observations and Derived Bedload Sediment Rates" Journal of Marine Science and Engineering 10, no. 2: 223. https://doi.org/10.3390/jmse10020223

APA Style

Kalra, T. S., Suttles, S. E., Sherwood, C. R., Warner, J. C., Aretxabaleta, A. L., & Leavitt, G. R. (2022). Shoaling Wave Shape Estimates from Field Observations and Derived Bedload Sediment Rates. Journal of Marine Science and Engineering, 10(2), 223. https://doi.org/10.3390/jmse10020223

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