Remote Sensing of Wave Overtopping on Dynamic Coastal Structures

: The development of coastal regions combined with rising sea levels is leading to an increasing risk of coastal ﬂooding caused by wave overtopping of natural beaches and engineered coastal structures. Previous measurements of wave overtopping have been obtained for static coastal structures using ﬁxed current meters and depth sensors or tanks. These are unsuitable for dynamically stable coastal protection structures however, because the geometry of these structures is expected to evolve under wave action. This study investigates the potential to use elevated 2D laser scanners (Lidar) to remotely sense the ﬂow volumes overtopping the time-varying crest of a porous dynamic cobble berm revetment. Two different analysis methods were used to estimate the wave-by-wave overtopping volumes from measurements of the time-varying free surface elevation with good agreement. The results suggest that the commonly used EurOtop parameterisation can be used to estimate overtopping discharge to an acceptable precision. An advantage of the remote sensing approach reported here is that it enables the spatial distribution of overtopping discharge and inﬁltration rate to be measured. It was found that the overtopping discharge on a porous dynamic revetment decays rapidly landward of the structure crest, and that this has implications for safety and structure design.


Introduction
Continued population growth combined with sea-level rise and climate change-related increases in storm intensity are placing ever-greater stress on the coastal zone, leading to increased reliance on existing coastal protection structures or the installation of new defences. Management of wave overtopping is a key criterion in the design of coastal protection structures, which must be designed to restrict overtopping to tolerable levels during design storm conditions.
Historical design guidance has typically focussed on limiting time-averaged overtopping discharge to minimise coastal flood risk. However, experimental evidence suggests that the peak discharge during individual overtopping events can be several orders of magnitude larger than the mean value and poses substantial hazards to people or vehicles [1].
As a result, recent design guidance also specifies acceptable limits on peak overtopping discharge for different structure types and safety considerations [2].
Wave overtopping of fixed coastal structures is typically measured using overtopping tanks which capture the flow volumes that pass over the structure crest in the laboratory [3][4][5]. Overtopping volumes are then estimated by measuring the water volume or mass collected within the overtopping tank. Field measurement of wave overtopping using tanks has also been undertaken but is limited to a relatively small number of studies [4,6,7] due to the practical challenges of applying this method in field conditions. A recent alternative to this approach is the WireWall system which uses a 3D array of capacitance wires mounted in a frame at the crest of a coastal structure [8]. This approach is designed to overcome some of the limitations of overtopping tanks which include tank capacity limits, challenges in measuring longshore variability and the lack of depth and velocity data for the captured overtopping flows. These approaches are limited to static structures. Here, we explore the use of Lidar to determine overtopping rates on dynamic coastal structures or dynamic morphology.
Dynamic cobble berm revetments are coastal structures that are designed to mimic the gravel barriers which form around the high tide shoreline of composite beaches and are fronted by a typically dissipative sand beach which is exposed at lower tide levels [9,10]. These ridges, whether natural or artificial, consist of pebbles or cobbles and are expected to continuously reshape under wave action, but maintain a coherent structure which provides overtopping protection to the hinterland [11]. Because the geometry of such structures varies naturally in response to wave forcing and water levels, it is not possible to design a fixed crest elevation to provide a defined level of wave overtopping protection. Nonetheless, it is important to understand the overtopping performance of these structures to determine their suitability for a particular location. Measurement of overtopping on a non-fixed structure is challenging however because the cross-shore location and elevation of the structure crest constantly changes during energetic overtopping conditions and as such, the use of fixed devices such as overtopping tanks or WireWall is unsuitable.
In their review of beach overtopping and overwash processes, Matias and Masselink [12] note that only a small number of researchers have presented field measurements of overtopping hydrodynamics on sand and gravel barriers [13][14][15][16][17]. Furthermore, these are typically confined to overtopping velocities and depths for a small number of events due to the challenges of obtaining measurements on a mobile feature during the storm conditions generally associated with overtopping events. It is noted that on natural barriers, the term "overwash" is typically used to describe flows which substantially exceed the crest and deposit sediment onto the backbarrier, while "overtopping" is used to describe events which just reach the crest and tend to lead to crest buildup [18]. Hereafter in this paper, the term "overtopping" will be used to describe flows which exceed the crest of both natural barriers and coastal structures.
The most comprehensive measurements of overtopping hydrodynamics on gravel and sand barriers were obtained during the large-scale BARDEX laboratory experiments [19,20]. These studies used video and an array of ultrasonic altimeters to obtain wave-by-wave measurements of flow velocities, intrusion distances and discharge, taking into account the variation in barrier crest location as the barrier evolved under wave action.
Hofland et al. [21] demonstrated the use of a single 2D Lidar to measure "virtual" wave overtopping volumes and discharge. This was done by integrating the flow depths landward of a chosen "virtual crest" location on a non-overtopped 1:3 planar dike slope in a large-scale laboratory flume. They hypothesised that the peak volume for each virtual overtopping event would represent a reasonable estimate of the volume of water that would have overtopped the virtual crest if the slope were truncated at that location. The results obtained using this method compared well to established parameterisations for wave overtopping on a plane slope. The approach was further developed by Oosterloo et al. [22] for application to oblique wave overtopping. Almeida et al. [17] used a different approach to estimate overtopping discharge at a large gravel barrier at Loe Bar, England using 2D Lidar. They assumed that a timeseries of overtopping discharge could be estimated as the product of the flow depth at the crest (h c ) and the shoreline velocity (u s ) and reported the first field measurements of overtopping discharge by individual waves for 12 events, recording peak values of up to 780 L/s/m.
An understanding of wave overtopping on dynamic coastal structures is essential to enable robust design, however minimal useful data currently exists due to limitations in traditional overtopping measurement techniques to obtain reliable laboratory or field data on a non-fixed structure. This paper will investigate the ability of 2D scanning Lidar to measure the spatial distribution of overtopping discharge on a dynamic cobble berm revetment structure using data from two large-scale laboratory experiments. The measured overtopping volumes are compared to the general wave overtopping formulation presented in [2] and the relationship between the overtopping rate and the deficit in freeboard is investigated.

Experimental Setup
Two experiments (DynaRev1 and 2) designed to investigate the performance of a dynamic cobble berm revetment under a rising water level and varying wave conditions were undertaken in the Large Wave Flume (Großer Wellenkanal, GWK) in Hannover, Germany. The GWK flume is 309 m long, 7 m deep and 5 m wide and equipped with a combined piston-flap type wavemaker. All coordinates are given as the distance from the wave paddle rest position (x = 0 m) and elevation above the horizontal flume bed (z = 0 m). DynaRev1 (1DR) is described in detail in [10,23] and was completed during August and September 2017. In summary, a dynamic cobble berm revetment with an initial slope of 1:6 (tanβ = 0.167) was constructed using well-sorted, rounded cobbles (D 50 = 63 mm, D 85 /D 15 = 1.32, porosity, ∅ = 0.41) on a sand beach profile (D 50 = 0.33 mm) which had developed after 20 h of constant wave forcing (H s = 0.8 m, T p = 6.0 s) (see Figure 1a). The water level was then increased in steps of 0.1 m from z WL = 4.6 m to 4.9 m, with 7 h of testing at all water levels except z WL = 4.9 m which lasted 17 h. The experimental testing was broken down into a series of test runs ranging in duration from 20 min to 3 h. At all water levels, the same two-hour long wave paddle signal was generated to produce an identical timeseries of waves at the wave paddle, taking water depth into account. Following the water level rise increments, the short-term response of the revetment was measured at the final water level (z WL = 4.9 m) for a range of different wave conditions (0.00791 < H o /L o < 0.0142) during a further 13 h of testing (see Figure 2a). The final profile of the revetment and underlying sand beach are shown in Figure 1a. It was observed that there was substantial loss of sediment from beneath the revetment during the experiment, an increase in the gradient to tanβ = 0.24, and only a small rise in crest elevation (see Bayle et al., 2020).
DynaRev2 (2DR) was completed during December 2019 and was almost identical to DynaRev1 except that the revetment was constructed using poorly-sorted angular pebbles and cobbles (D 50 = 44 mm, D 85 /D 15 = 3.79) leading to a lower porosity, ∅ = 0.35. The beach and revetment profiles were translated 11 m further from the paddle due to a reduced availability of sand and a slightly different test series order was used at the final water level (Figure 2b). Due to the greater difficulty of shaping the irregular stone used in DynaRev2, the initial revetment profile was more irregular than that for DynaRev1. Note also that the wider grading and greater angularity (leading to increased interlocking) of the stone used during DynaRev2 led to the development of a much more pronounced crest feature and increased gradient (tanβ = 0.25) once waves began to overtop the structure (see x = 272.8 m in Figure 1b).
During both experiments, an array of three SICK LMS511 2D scanning Lidar were installed in the flume roof to measure surf and swash processes at 25 Hz and with an angular resolution of 0.167 • . In this paper, only the most landward Lidar which was just seaward of the revetment toe and covered approximately a 30 m swath was used (x = 255 m, z = 11.8 m (1DR), x = 266 m, z = 11.8 m (2DR)). The analysis presented in this manuscript focusses on the most elevated water level testing at z WL = 4.8 m and 4.9 m when overtopping of the revetment crest occurred. In both experiments, the range of wave conditions tested was relatively limited (Figure 2a,b).
However, the variability of the revetment geometry and crest elevation meant that a large range of relative freeboard (R c /H o ) and deficit in freeboard values (R def ) were captured: where z c is the revetment crest elevation, z SWL is the elevation of the still water level and R 2% is the predicted 2% exceedance wave runup elevation. Experimentally and theoretically, on a plane slope, the overtopping scales with the square of the deficit in freeboard [24]. In this paper, R 2% above the measured MWL at the shoreline is estimated using a runup equation developed specifically for composite beaches using data from field experiments [10,17,25] and the DynaRev experiments reported here: where H m0 is the spectral significant wave height at the revetment toe and tanβ is the revetment gradient.

Lidar Measurements of Wave Overtopping
The 2D scanning Lidar used in the experiments reported here enable high spatiotemporal resolution measurements of both the time-varying bed elevation (when the bed is exposed) and the swash surface elevation (when the bed is submerged). The Lidar dataset was converted from polar to cartesian coordinates and interpolated onto a horizontal grid (∆x = 0.1 m). These data were then post-processed using the method of Almeida et al. [26] to generate separate timeseries of: (1) continuous beach profile measurements z(x, t) updated on a swash-by-swash timescale from which the time-varying crest location can be obtained (see Figure 1); and (2) swash hydrodynamics including swash/overtopping flow depths h(x, t) and horizontal shoreline position X s (t) (see Figure 3a). (d) Timeseries of overtopping volume using the peak volume (V v ) and sum of positive discharge (V uh ) methods. The red markers indicate the overtopping volume estimates V v,max = 0.087 m 3 (red diamond) and V uh,max = 0.096 m 3 (red circle) for this event. Note that V v is calculated relative to the bed elevation preceding the overtopping event and so does not return to zero in this case because the overtopping leads to minor bed accretion.
Wave overtopping events were identified by applying a zero-upcrossing method to the shoreline timeseries with the horizontal crest location x c (t) as the datum. Only events where the peak water depth at the crest exceeded 0.05 m were considered as valid overtopping events in the analysis presented here. As such the overtopping percentage values presented in Figure 2e,f are smaller than those shown by Bayle et al. [10] who did not apply an overtopping depth threshold.
Due to the constantly evolving nature of the revetment geometry, all experimental runs were broken down into 20 min "segments". During each segment, the cross-shore location of the crest (x c ) was assumed to be constant and was determined from the Lidar-derived bed profile data at the mid-point of the segment (e.g., after 10 min). The crest was taken as the location were the seaward slope of the revetment transitioned to either a flat structure top (1DR, see Figure 1a) or a peaked crest which sloped landward in the positive x-direction (2DR, see Figure 1b). In general, the crest also corresponded with the highest point on the revetment profile. While the cross-shore location of the crest and revetment gradient (tanβ) was assumed to be fixed within each segment, the crest elevation (z c ) was updated on a wave-by-wave basis using the Lidar-derived bed profile data. The 20 min segment length was chosen to minimise the effect of crest movement (|∆x c | < 0.5 m; |∆z c | < 0.05 m) and any variability in the revetment slope, but is shorter than would typically be used to obtain stable overtopping statistics on a fixed structure. Over all 20 min segments, an average of 205 swash events were measured, with 26.7 (13%) of these events overtopping the structure crest. A sensitivity analysis using 40 and 60-min segments was undertaken which showed comparable results to those presented below.
Wave overtopping volumes were estimated on a wave-by-wave basis using two different approaches. The "peak volume" method is comparable to the "virtual" overtopping volume approach [21,22]. Timeseries of overtopping flow volume were produced by integrating the timeseries of flow depths landward of the crest position x c (t): where h c represents the depth at the crest, h c+1 the depth at the next landward position and h s the depth immediately seaward of the time varying shoreline position X s (t), thus h s+1 = 0 m. For every identified overtopping event, the overtopping volume V v,max is assumed to be the peak value of the overtopping volume timeseries V v (t) (see Figure 3). Note that [21] did not actually measure overtopping of a structure crest, but analyzed peak flow volumes landward of a fixed point on a non-overtopped planar slope. For the experimental setup investigated here, where waves overtopped the crest of a porous structure, there are a number of potential sources of error in this approach: 1.
An air void fraction α = 0 is assumed, however for violent overtopping flows on an irregular stone surface this is unlikely to be the case and will lead to an overestimate of overtopping volume. Indeed, very rapid overtopping flows were occasionally observed to "jump" as they passed over the revetment crest during 2DR leading to a short-lived air cavity beneath the overtopping flow. This effect led to anomalously large maxima in the depth timeseries during initial overtopping which were removed by discounting volume peaks during the first 0.25 s of these events.

2.
Infiltration into the porous revetment structure will mean that the estimated peak swash volume is smaller than the total amount of water that overtopped the crest. However, from an engineering perspective, the volumes presented here represent the flow volume actually present landward of the structure crest that could lead to safety concerns or inundation risk.
The second "sum of positive discharge" method used to estimate wave overtopping volume is based on [17] which estimated instantaneous overtopping discharge on a gravel barrier as: where h c is updated at every timestep and us is the shoreline velocity calculated as the first derivative of the shoreline position timeseries: By summing the flow volume passing the crest since the start of each overtopping event, a timeseries of overtopping flow volume can be estimated: and the overtopping volume for each identified overtopping event is estimated as: where T + is the duration of positive shoreline velocity values (u s > 0 m/s) for each overtopping event.
This methodology is similar to that used by previous researchers [27] to estimate wave overtopping of fixed structures in the laboratory where velocimeters are used to measure flow velocity over the structure crest. For the experimental setup investigated here, a number of potential sources of error exist:

1.
As with the peak volume method, an air void fraction α = 0 is assumed, however for violent overtopping flows on an irregular stone surface this is unlikely to be the case. The presence of air voids will increase h c and lead to an overestimate of overtopping volume. Note that the flow "jumps" described above occur landward of the structure crest and so the large depths associated with this effect were not measured at the crest location.

2.
The flow depth at the structure crest was sometimes observed to be noisy during energetic overtopping events due to splashing, potentially leading to overestimates of the depth and hence overtopping volume. The depth data was despiked to remove short-lived depth maxima in order to minimise this effect.

3.
The shoreline velocity u s is assumed to be representative of the depth-averaged flow velocity at the barrier crest u c . Schüttrumpf and Oumeraci [27] demonstrated that bore front velocity compared well with flow velocity measurements (micro-propeller) for waves overtopping an impermeable sea dike in the laboratory. In the current experiment, the irregular and porous nature of the structure may be expected to lead to enhanced deceleration of the swash front after overtopping. As a result, u s is expected to provide a reasonable estimate of u c during the initial stages of overtopping but underestimate u c as flow reversal is approached.
For both methods, overtopping discharge q v ("peak volume") or q uh ("sum of positive discharge") during a 20 min segment (t = 1200 s) is defined as: where V max is V v,max for the peak volume method and V uh,max for the sum of positive discharge method.

Comparison of Overtopping Measurement Methods
Hofland et al. [21] and Oosterloo et al. [22] demonstrated that the peak volume method could be used to obtain valid overtopping measurements on relatively smooth, impermeable dike structures, however it has not been tested on a constantly evolving, porous structure such as a dynamic cobble berm revetment. Due to the limitations of fixed instrumentation for measuring overtopping on a dynamic revetment it was not possible to validate the Lidar-based methodologies outlined in Section 2.2 against traditional overtopping measurement techniques. Consequently, the following section compares results obtained using the two different methodologies. Figure 4 compares the discharge measured using the peak volume and sum of positive discharge methods for all 20 min segments during both 1DR and 2DR. For 1DR, there is reasonably good agreement, with all values within a factor of two of each other. There is a tendency for q v < q uh for higher values of discharge (q v > 2.5 l/s/m). For 2DR, q v > q uh for lower values of discharge (q v < 2.5 l/s/m), possibly due to higher levels of aeration and splashing in the flow caused by overtopping events passing over the prominent crest. Figure 5 compares the estimated overtopping volume for all individual overtopping events measured during both experiments (3006 events (1DR); 2665 events (2DR)). While the data are scattered, there is a reasonable level of agreement, with 82% and 74% of values within a factor of two agreement for 1DR and 2DR respectively.
The results presented above indicate that while there is not total agreement between the two methods, the majority of values are within a factor of two agreement which is considered an acceptable result given the complexity of overtopping flows on a rough, porous and constantly evolving structure and the inherent uncertainties in the methodologies. For the remainder of the manuscript, the peak volume method is used to estimate overtopping volume or discharge as this method is simpler and is likely to be more reliable because it has been previously validated on impermeable structures and is less sensitive to spikes in the depth or velocity timeseries.

Overtopping Discharge at the Revetment Crest
The timeseries presented in Figure 2c,d demonstrates the variability of the crest elevation throughout the two experiments. As noted in [10], the rounded nature of the cobbles used in 1DR limits crest build-up and so there is only a 0.22 m increase in crest elevation during the test runs analysed here (Figure 2c). By contrast for 2DR, the poorly sorted angular stone allowed a peaked crest to form once the structure was regularly overtopped (z WL ≥ 4.8 m) and this feature adjusted rapidly to changes in water level and wave conditions. It is evident that the crest elevation increases during the first few hours after the changes in water level to z WL = 4.8 m (0 h) and z WL = 4.9 m (7 h) leading to a gradual reduction in the freeboard deficit which can be observed in Figure 2e-h. The observed reduction in freeboard deficit with time leads to an associated decrease in the percentage of wave overtopping (Figure 2e,f) and overtopping discharge (Figure 2g,h), both of which approximately stabilise between 12 and 24 h at mean values of 8.6% (σ = 1.5%) and 2.03 l/s/m (σ = 0.62 l/s/m) respectively. Note that predicted runup elevation for 2DR is consistently larger than for 1DR despite using the same wave conditions at the wave paddle due to slightly larger depths seaward of the revetment and hence greater wave heights measured at the revetment toe.
EurOtop [2] provides a widely used general formulation for the prediction of wave overtopping on coastal structures: where H m0 is the spectral significant wave height at the structure toe, α is the revetment gradient, ξ m−1,0 is the surf-similarity parameter calculated using the spectral wave period T m-0,1 and γ b , γ f , γ β , γ v are influence factors to account for a berm, roughness, wave obliquity and a wall at the end of the slope respectively. Figure 6 shows the overtopping discharge normalised using the approach presented by EurOtop [2] as a function of the normalised freeboard deficit R*: While there is substantial scatter (common for large-scale overtopping measurements), there is evidence of a linear relationship between R* 2 and the normalised overtopping discharge as previously noted by Ibrahim and Baldock [24]. Figure 7 presents a scatter plot comparing the measured discharge rates with those predicted by the EurOtop equation where all reduction factors are taken as unity except for the roughness factor γ f . Little previous research has been done to determine the roughness factor for cobble slopes and no design values are provided in [2], however recent numerical work [28] suggests values of γ f between 0.62 and 0.75 for a dynamic cobble berm revetment. Their work suggested that overtopping increases as the porosity of the revetment reduces due to sand infilling and γ f ∝ 1/p v where p v is the relative pore volume of the revetment. It was found that the values of γ f that best fitted the present data were γ f = 0.50 and γ f = 0.70 for 1DR and 2DR respectively. It is suggested that the lower value for 1DR accounts for the higher porosity of the revetment (∅ 1DR = 0.41; ∅ 2DR = 0.35). Using these values for γ f enables predictions of overtopping discharge to within a factor of three which is the approximate accuracy quoted for Equation (10) given in [2]. Furthermore, in agreement with results from other types of structures, the largest deviations tend to be at smaller discharges, particularly for 2DR where overtopping of small discharges is consistently underestimated. Figure 8 shows the exceedance probabilities for individual wave overtopping discharge measured between 12 and 24 h of testing when z WL = 4.9 m and the crest elevation and cross-shore location were relatively constant (R c,1DR = 0.54 m (σ = 0.026 m) and R c,2DR = 0.61 m (σ = 0.022 m)). Note that the relative freeboard R c /H m0 is almost identical for the two datasets at 1.26 and 1.27 for 1DR and 2DR respectively. For both experiments, the measured data agrees well with a 2-parameter Weibull distribution using the shape and scale factors calculated using the guidance in EurOtop [2] for probabilities of exceedance larger than approximately 10%, however the probability of larger discharges is smaller than that suggestedin [2]. This was also observed by Oosterloo et al. [22] and in this case may be due to greater relative underestimation of the largest overtopping volumes caused by greater infiltration over the longer overtopping duration.

Spatial Distribution of Overtopping Discharge
An advantage of using remote sensing methods to estimate wave overtopping is that the location of measurement points is not fixed, making it possible to investigate the spatial distribution of overtopping flows. This is particularly relevant on a porous structure where the mass of flow passing over the structure is not conserved due to high levels of infiltration (unlike on impermeable sea dikes, e.g., [29]).
The analysis below focusses on the 1DR experiment where the structure crest is relatively flat and so there are reduced complications caused by flow acceleration on the landward slope of the revetment crest in 2DR. Figure 9a presents the cross-shore distribution of q v relative to the structure crest for all 20 min segments during 1DR where the entire structure was not overtopped such that flow reached the sand behind the revetment (see Figure 1a). In Figure 9a, the values of q v at each cross-shore position (∆x = 0.1 m) are calculated using Equation (9) by evaluating the maximum overtopping volume (V v,max ) captured landward of each measurement location for every overtopping event. Due to infiltration, the maximum overtopping volumes and hence q v decays rapidly landward of the revetment crest.
The normalised discharge distributions in Figure 9b indicate a level of self-similarity and can be described using power law decay, with the best-fit power law decay function found to be: where q c is the overtopping discharge at the crest, x c is the horizontal location of the revetment crest and R x,max − x c represents the incursion distance, defined here as the horizontal distance between the revetment crest and the maximum recorded shoreline position R x,max during each 20 min segment. While the landward sloping crest of the 2DR revetment will influence the spatial distribution of overtopping discharge, Figure 10 indicates that the distributions are relatively similar, if slightly more variable than for 1DR and Equation (12) can be reasonably applied to the 2DR data. Figure 9. (a) Cross-shore distribution of overtopping discharge q v relative to the revetment crest location during 1DR where the data from each 20 min segment is represented using a different colour. Note that 20 min segments where the entire structure was overtopped such that overtopping events reached the sand beach landward of the revetment (x = 264.1 m) were excluded from this analysis. (b) Overtopping discharge normalised by the overtopping discharge at the crest as a function of distance from the revetment crest normalised by the measured incursion distance for each segment. The grey lines indicate the data for each segment and the solid black line shows the best fit power law and exponential decay functions respectively. Figure 10. Overtopping discharge normalised by the overtopping discharge at the crest during 2DR as a function of incursion distance for each segment. The grey lines indicate the data for each segment and the solid black line shows Equation (12). Note that 20 min segments where the entire structure was overtopped such that overtopping events reached the sand beach landward of the revetment (x = 275.1 m) were excluded from this analysis.

Implications for Dynamic Cobble Berm Revetment Design
The ability to measure the spatial variation of overtopping discharge landward of the crest is a significant advantage of the Lidar-based method used here. The results presented in Figure 10a indicate that discharge q v decays rapidly landward of the crest, leading to a substantial reduction in the risk to people safety or inundation with distance from the crest. For example, if the current large-scale experiments are considered to be at prototype scale, Figure 4 shows that the time-averaged overtopping discharge at the crest, q v exceeded the recommended tolerable limit for an aware pedestrian of 0.1 l/s/m [30] during 88% of all 20 min segments during 1DR. However, this value was not exceeded during any segment only 2 m landward of the crest.
The decay of overtopping discharge landward of a dynamic cobble berm revetment crest is described well using the power law relationship given in Equation (12) for the current experiment. Further it was demonstrated that overtopping discharge at the crest can be predicted using the EurOtop equation to within approximately a factor of 3 if an appropriate value of γ f is used (Figure 7). This performance is comparable to that for static coastal structures. Thus, if the incursion distance (R x,max − x c ) can be reasonably estimated, it is possible to predict the spatial distribution of overtopping discharge. Figure 11 indicates that the incursion distance beyond the revetment crest increases with freeboard deficit as might be expected, but the data from the present study is very scattered. This scatter and difference between results for 1DR and 2DR is likely due to the variable backbarrier geometry. For 2DR, the landward sloping gradient behind the crest means that any events which overtop the crest will run down the backbarrier and this leads to a minimum measured incursion distance of 1.5 m. For 1DR, the barrier crest is approximately flat and as such there is no physically defined minimum incursion distance and the relationship between incursion and runup deficit is stronger. While scattered, the results indicate that it may be possible to approximate (R x,max − x c ) as a function of R def and use this value along with Equations (10) and (12) to estimate the spatial distribution of overtopping discharge. The best fit linear regression for the current data suggests R x,max − x c ≈ 11.4R def , however it is suggested that additional experimental data is required to better understand this relationship, in particular the effect of barrier geometry. The ability to predict the spatial distribution of wave overtopping discharge would aid the dynamic cobble berm design process as it would enable the specification of a design crest elevation and structure width that would maintain tolerable overtopping discharge at a particular location. Such a method would also enable more robust safety assessments of existing revetment structures and naturally occurring composite beach ridges. It is noted that the dynamic nature of dynamic cobble berm revetments means that the structure geometry would not be constant and so a conservative design process would be required.

Infiltration of Overtopping Volume
The decay of overtopping discharge landward of the revetment crest observed in Figures 9 and 10 is primarily a result of volume loss through infiltration into the porous revetment. This infiltration of the overtopping flow is an important characteristic contributing to the stability of coastal barriers and dynamic revetments. However, if infiltration rates are too great this could lead to scour below the revetment and subsequent sub-surface erosion of the barrier.
The ability of the Lidar to capture flow volumes enables the time decay of volume after overtopping to be measured in order to estimate the infiltration rate. Figure 3d shows the typical decay of the overtopped flow volume landward of the crest due to infiltration into the bed between t = 4.3 s and t = 5.7 s and indicates a gradually reducing rate of volume loss with time as the depth of the overtopped flow volume on the revetment crest reduces. Figure 12 presents the infiltration flow rate as a function of the mean depth of the overtopping flow, where infiltration rate S i (t) is calculated during the period after the overtopping event has reached its landward limit and the overtopped flow volume is decaying (e.g., t = 4.3 s to t = 5.7 s in Figure 3d) as: where l is the cross-shore distance landward of the revetment crest covered by time-varying overtopping flow. Figure 11. Incursion distance relative to the cross-shore crest location as a function of the freeboard deficit. Black crosses represent data for 1DR and the 2DR data are shown as grey circles. Note that 20 min segments where the entire structure was overwashed such that overtopping events reached the sand beach landward of the revetment were excluded from this analysis. The analysis focusses on overtopping events with a single peak in the volume timeseries that occurred in the 3-h period between 18 and 21 h for both experiments and discounts any overtopping events where multiple overriding waves overtopped the revetment, because this greatly influences the cross-barrier decay of the flow volume. Furthermore, only events with an intrusion distance larger than 0.5 m but smaller than the structure crest width were considered. Figure 12 indicates that instantaneous infiltration rates averaged over the overtopping volume are two to three orders of magnitude larger than previously observed on sand beaches [31], with median values of −0.07 ms −1 (−0.11 ms −1 ) and rates of up to −0.87 ms −1 (−1.37 ms −1 ) measured during experiment 1DR (2DR). While the data is scattered, there is a clear trend for higher infiltration rates with larger flow depths when the pressure head at the bed is relatively large. Infiltration rates tend to be larger for 2DR, primarily because of larger measured flow depths. It may be expected that for a constant value of mean flow depth h, infiltration rates would be higher for 1DR where the well-sorted, rounded cobbles for 1DR would lead to a more porous structure and higher infiltration rates (∅ 1DR = 0.41; ∅ 2DR = 0.35) . This is not obvious in the results presented here, though the scatter is substantial. Due to the intermittent nature of wave overtopping, it is assumed that the bed is unsaturated during initial overtopping. Packwood [32] presented a model to estimate the influence of an unsaturated porous bed on wave runup on sandy beaches which has been successfully employed in the XBeach-G model to calculate infiltration rate S i . If we assume hydrostatic pressure variation with depth: where K is the hydraulic conductivity, h is the flow depth and δ w f is the time-varying thickness of the wetting front (the depth of the infiltrating flow volume below the revetment surface). Note that due to the presence of cobbles, the wetting front thickness will increase at a greater rate than the decrease in overlying overtopping volume depth. Equation (14) includes a dependence of S i on the overtopping volume flow depth h which is observed in the current data in Figure 12. However, Equation (14) assumes that the cobble bed is unsaturated and as such is ventilated beneath the wetting front, leading to the result that S i → K as the water depth h → 0. Values of hydraulic conductivity for cobble beds are not well documented, but K is expected to be approximately 0.1 to 1 m/s [33,34]. The current data does not appear to agree with this result as observed in Figure 12 which indicates that S i → 0 for h → 0. Several possible reasons for this exist: (1) Equation (13) is based on conservation of mass and assumes that volume is only lost through infiltration into the cobbles, as such it must be applied over the entire overtopping volume. This means that the infiltration rates presented are averaged over cross-shore distances of the order of metres and local rates may be different. (2) The cross-shore resolution of the Lidar data (0.1 m) is not sufficient to resolve individual cobbles, consequently the accuracy of the depth estimates is limited by the irregularity of the cobble surface and the effect of this uncertainty will be largest at small values of h. Furthermore, the values of depth presented are also averaged over the entire time-varying overtopping volume. (3) Due to the relatively low gradient of the underlying sand bed, infiltrating water may not always drain seaward between consecutive overtopping events meaning that the bed is not always ventilated. (4) It is possible that due to the small thickness of the cobble layer overlying the sand bed (generally < 0.5 m) and the low gradient of the sand, the wetting front for larger events reaches the sand bed during the early stages of overtopping, meaning that the bed is not ventilated as assumed by Equation (14).

Conclusions
The results presented in this paper indicate that 2D scanning Lidar can be used to capture wave overtopping of porous dynamic coastal structures with evolving geometry. Two different analysis methods were used to estimate wave-by-wave overtopping volumes and agreed to within a factor of two for the majority of events, with better agreement for time-averaged discharge. This result gives confidence that Lidar is a valid tool for the measurement of wave overtopping, however further validation against overtopping tanks on a fixed, porous structure would be beneficial.
The commonly used EurOtop parameterisation for estimating overtopping discharge was found to provide estimates of time-averaged wave overtopping to within a factor of three, which is comparable to the expected accuracy for typical fixed coastal protection structures. Use of an appropriate friction factor is essential for robust application of the EurOtop equation. For the current dataset the friction factor was found here to be smaller for round, well-sorted compared to angular, poorly-sorted cobbles. It is hypothesised that this difference is due to the increased porosity of the rounded cobble revetment leading to higher infiltration rates. To investigate this further, the time-variation of the overtopped flow volume was used to obtain the first estimates of infiltration rate into cobble beaches. These data indicate infiltration rates 2-3 orders of magnitude larger than on sand beaches, but there was no evidence of higher infiltration rates for the rounded cobble revetment.
The ability to capture the time-variation of the entire overtopping volume, rather than just flow volume or velocity and depth measurements at one location, enables the cross-shore distribution of overtopping discharge to be measured. It is demonstrated that the overtopping discharge decays rapidly landward of the crest due to the porous nature of the structure and can be described using a power law function. This is valuable for the design of porous structures where the rapid decay of overtopping discharge means that inundation and safety concerns can potentially be reduced by using an appropriate structure width and setting back access routes and infrastructure.

Data Availability Statement:
The data from the DynaRev1 experiment analysed within this manuscript is available for download from DOI 10.5281/zenodo.3889796. Additional metadata is provided within each *.mat file detailing how the data from each instrument is stored. Note also that all raw, unprocessed data is available at DOI 10.5281/zenodo.3855650.