# High-resolution, large-scale laboratory measurements of a sandy beach and dynamic cobble berm revetment

Chris E. Blenkinsopp, Paul M. Bayle, Daniel C. Conley, Gerd Masselink, Emily Gulson, Isabel Kelly, Rafael Almar, Ian L. Turner, Tom E. Baldock, Tomas Beuzen, et al.

► To cite this version:

Chris E. Blenkinsopp, Paul M. Bayle, Daniel C. Conley, Gerd Masselink, Emily Gulson, et al. High-resolution, large-scale laboratory measurements of a sandy beach and dynamic cobble berm revetment. Scientific Data, 2021, 10.1038/s41597-021-00805-1. hal-04720546

HAL Id: hal-04720546

https://hal.science/hal-04720546v1

Submitted on 4 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC DATA 110110
0111101
11011110
011101101

<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

DATA DESCRIPTOR

# High-resolution, large-scale laboratory measurements of a sandy beach and dynamic cobble berm revetment

Chris E. Blenkinsopp ID 1✉, Paul M. Bayle<fnref n="1" />, Daniel C. Conley ID 2, Gerd Masselink², Emily Gulson², Isabel Kelly², Rafael Almar ID 3, Ian L. Turner ID 4, Tom E. Baldock 5, Tomas Beuzen 4,6, Robert T. McCall 7, Huub Rijper 7,8,9, Ad Reniers 9, Peter Troch 10, David Gallach-Sanchez ID 10,11, Alan J. Hunter 12, Oscar Bryan 12, Gwyn Hennessey 1, Peter Ganderton², Marion Tissier 9, Matthias Kudella 13 & Stefan Schimmels ID 13

High quality laboratory measurements of nearshore waves and morphology change at, or near prototype-scale are essential to support new understanding of coastal processes and enable the development and validation of predictive models. The DynaRev experiment was completed at the GWK large wave flume over 8 weeks during 2017 to investigate the response of a sandy beach to water level rise and varying wave conditions with and without a dynamic cobble berm revetment, as well as the resilience of the revetment itself. A large array of instrumentation was used throughout the experiment to capture: (1) wave transformation from intermediate water depths to the runup limit at high spatio-temporal resolution, (2) beach profile change including wave-by-wave changes in the swash zone, (3) detailed hydro and morphodynamic measurements around a developing and a translating sandbar.

## Background & Summary

High quality field and numerical investigations are providing new insights into a wide variety of coastal processes and coastal protection solutions<fnref n="1" />,². However, numerical modelling approaches are not yet capable of accurately reproducing all coastal hydro and morphodynamic phenomena, and the difficulties involved in capturing field data in the desired wave, tide and wind conditions mean that controlled laboratory wave flume experiments remain extremely valuable. Large-scale experiments³,⁴ are particularly valuable as they mostly avoid scaling issues, and improvements in the instrumentation and measurement techniques available mean that the quality and resolution of data continues to improve and provide new insights.

The DynaRev experiment was designed to investigate the response of a sand beach and the resilience of a dynamic cobble berm revetment to constant wave forcing and a rising water level at large-scale in a controlled laboratory environment through high spatio-temporal resolution morphology measurements (Fig. 1). A dynamic cobble berm revetment is a nature-based coastal protection approach which consists of a cobble ridge constructed around the high tide runup limit to artificially mimic composite beaches⁵. This commonly occurring beach type

<footnote n="1">¹Centre for Infrastructure, Geotechnics and Water Engineering, Department of Architecture and Civil Engineering, University of Bath, Bath, BA2 7AY, UK. ²School of Biological and Marine Sciences, Plymouth University, Drake Circus, PL4 8AA, Plymouth, UK. ³IRD-LEGOS, UMR5566, 18 Av. Edouard Belin, 31400, Toulouse, France. ⁴Water Research Laboratory, School of Civil and Environmental Engineering, UNSW Sydney, NSW, 2052, Australia. ⁵School of Civil Engineering, University of Queensland, Brisbane, QLD, 4072, Australia. ⁶Department of Statistics, University of British Columbia, Vancouver, Canada. ⁷Department of Marine and Coastal Systems, Deltares, Boussinesqweg 1, 2629 HV, Delft, the Netherlands. ⁸Royal Boskalis Westminster N.V., Rosmolenweg 20, 3356 LK, Papendrecht, the Netherlands. ⁹Faculty of Civil Engineering and Geosciences, Delft University of Technology, Stevinweg 1, 2628 CN, Delft, The Netherlands. ¹⁰Department of Civil Engineering, Ghent University, Technologiepark 60, Ghent, B-9052, Belgium. ¹¹DEME Group, Scheldedijk 30, 2070, Zwijndrecht, Belgium. ¹²Department of Mechanical Engineering, University of Bath, Bath, BA2 7AY, UK. ¹³Forschungszentrum Küste (FZK), Merkurstraße 11, 30419, Hannover, Germany. ✉e-mail: c.blenkinsopp@bath.ac.uk</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** (a) Schematic of flume setup showing primary instrument locations (see Table 1). The yellow shaded area represents the sand volume and the dark grey shaded area is the permanent 1:6 impermeable slope. The black solid and dashed horizontal lines indicate the minimum ($z_{wl}$ = 4.5 m) and maximum ($z_{wl}$ = 4.9 m) water levels. (b) Close up of the dynamic cobble berm revetment geometry after construction corresponding to the grey box in (a). The minimum water level used for revetment testing ($z_{wl}$ = 4.6 m) is shown as a solid horizontal line and the dashed line indicates the maximum water level. The light grey region indicates the constructed dynamic revetment and the dot-dashed line shows the beach profile prior to revetment construction. (c) Photograph of the constructed dynamic revetment on the underlying sand beach. The yellow line indicates the initial line of the revetment crest.</figcaption></figure>

consists of a lower foreshore of sand and a backshore ridge constructed of gravel or cobbles that stabilises the upper beach and provides overtopping protection. Dynamic revetment structures contrast with static coastal defence structures as they are specifically designed to reshape under wave attack. In addition to the morphology data, high-resolution measurements of nearshore hydrodynamic processes were also collected.

DynaRev took place over a 2-month period from August to September 2017 in the 309 m long Large Wave Flume (Großer Wellenkanal, GWK), Hannover, Germany. A total of 141.6 hours of testing under wave action was completed. This testing comprised two "phases", with each phase being split into a series of "runs" varying from 20 minutes to 3 hours in duration. The beach profile was only reset between the two phases, thus all runs had a different antecedent morphology corresponding to the beach profile at the end of the preceding run.

**Phase SB - Unmodified sand beach response to a rising water level:** Starting with a plane 1:15 sand slope, the evolution of the beach profile was measured under constant wave forcing ($H_s$ = 0.8 m, $T_p$ = 6.0 s) for 20 hours. The mean water level in the flume was then raised from an initial elevation $z_{wl}$ = 4.5 m by a total of 0.4 m in incremental steps of 0.1 m (38 hours of water level rise testing). Following the completion of the water level rise increments, the short-term response of the beach was measured at the final water level ($z_{wl}$ = 4.9 m) for a range of different wave conditions expected to produce both erosion and accretion.

**Phase DR - Dynamic cobble berm revetment response to a rising water level:** Again starting with a manually reshaped 1:15 plane slope, a sand beach was measured as it evolved under the same constant wave conditions as used in Phase SB for 20 hours to provide a natural beach profile on which to construct the dynamic revetment. Following this, the same water level increments and test durations as for Phase SB were applied. Prior to the first water level increment, a cobble revetment was installed at the location of the sand beach berm and was designed such that its crest height was at the elevation of the $R_{2\%}$ runup level measured during Phase SB for the second water level increment to ensure significant overtopping as the water level was increased. The sand foreshore and dynamic revetment were then allowed to reshape under constant wave conditions over the remaining water level increments, with the test durations at each water level mirroring those in Phase SB (38 hours of water level rise testing). Finally, higher energy storm waves were used at the end of the final water level increment to investigate revetment resilience to higher energy conditions.

The availability to researchers of large-scale measurements of nearshore hydro and morphodynamics at the spatio-temporal resolution achieved during DynaRev is very limited. Potential uses for the datasets obtained during the DynaRev test program are wide-ranging and include: the assessment of dynamic cobble berm revetment performance⁶, the investigation of nearshore processes such as the formation and dynamics of nearshore sandbars⁷, the response of sandy coasts to a rising sea level⁸, morphology change in the swash zone⁹, wave-by-wave sediment transport rates¹⁰, air entrainment in breaking waves¹¹ and the development of numerical models¹².
<table><thead><tr><th>Abbrev.</th><th>Instrument</th><th>Purpose (measurement units)</th><th>x (m)</th><th>z (m)</th></tr></thead><tbody><tr><td>WG1</td><td>Wave gauge</td><td></td><td>50</td><td>—</td></tr><tr><td>WG2</td><td>Wave gauge</td><td rowspan="3">Array 1: Water surface elevation in the deep flume section, η (m)</td><td>51.9</td><td>—</td></tr><tr><td>WG3</td><td>Wave gauge</td><td>55.2</td><td>—</td></tr><tr><td>WG4</td><td>Wave gauge</td><td>60</td><td>—</td></tr><tr><td>WG5</td><td>Wave gauge</td><td></td><td>160</td><td>—</td></tr><tr><td>WG6</td><td>Wave gauge</td><td rowspan="3">Array 2: Water surface elevation in the deep flume section, η (m)</td><td>161.9</td><td>—</td></tr><tr><td>WG7</td><td>Wave gauge</td><td>165.2</td><td>—</td></tr><tr><td>WG8</td><td>Wave gauge</td><td>170</td><td>—</td></tr><tr><td>ADV1</td><td>Nortek Vector</td><td>Flow velocity, u, v, w (ms<sup>-1</sup>) – shoaling waves</td><td>180</td><td>2.5</td></tr><tr><td>ADV2</td><td>Nortek Vector</td><td>Flow velocity, u, v, w (ms<sup>-1</sup>) – surf zone</td><td>235</td><td>3.67</td></tr><tr><td>ADV3</td><td>Nortek Vector</td><td></td><td>242</td><td>4.22</td></tr><tr><td>WGADV1</td><td>Wave gauge</td><td>Water surface elevation at ADV1 location, η (m)</td><td>180</td><td>2.5</td></tr><tr><td>PTADV2</td><td>Pressure transducer</td><td>Pressure at ADV2 location, P (kPa)</td><td>235</td><td>3.67</td></tr><tr><td>PTADV3</td><td>Pressure transducer</td><td>Pressure at ADV3 location, P (kPa)</td><td>242</td><td>4.22</td></tr><tr><td>PT3</td><td>Pressure transducer</td><td>Pressure between the surf zone/bar processes instrument rigs, P (kPa)</td><td>231.7</td><td>4.13</td></tr><tr><td>LID1</td><td>SICK LMS511 2D Lidar</td><td rowspan="2">High spatio-temporal resolution water surface profile, η (m) – surf zone</td><td>230.04</td><td>11.76</td></tr><tr><td>LID2</td><td>SICK LMS511 2D Lidar</td><td>242.02</td><td>11.85</td></tr><tr><td>LID3</td><td>SICK LMS511 2D Lidar</td><td>Swash surface profile, η (m), Beach/revetment profile, z (m)</td><td>254.99</td><td>11.82</td></tr><tr><td>CAM</td><td>Vivotek IB9381-HT high resolution camera</td><td>Swash zone imagery</td><td>Adjustable (276–280 m)</td><td>11.8</td></tr><tr><td>MB</td><td>Reson 7125 Multibeam</td><td>Bubble cloud, Bathymetry, x,z (dB)</td><td>Adjustable</td><td>Adjustable</td></tr><tr><td>FARO</td><td>FARO Focus 3D (Lidar)</td><td>3D topography (m)</td><td>Adjustable</td><td>Adjustable</td></tr><tr><td>RFID</td><td>Instrumented cobbles</td><td>Cobble movement</td><td>97 cobbles placed at 3 depths along the revetment centreline</td><td></td></tr></tbody></table>

## Surf Zone Instrumentation

Rigs were reset to maintain constant instrument elevations above the bed at the start of every test, thus all elevations are presented in cm relative to the local bed and given the notation h.

<table><caption>**Table 1.** Summary of the measurement instruments deployed during the experiment including: Instrument type, measurement purpose, measurement units and primary instrument locations noting that some instruments were moved during the experiment as described in the manuscript.</caption><thead><tr><th>Abbrev.</th><th>Instrument</th><th>Purpose (measurement units)</th><th>x (m)</th><th>h (cm)</th></tr></thead><tbody><tr><td>PT1</td><td>Pressure transducer</td><td>Pressure, P (kPa)</td><td rowspan="4">226.5</td><td>45</td></tr><tr><td>OBS1</td><td>Optical backscatter sensor</td><td rowspan="2">Suspended sediment concentration, C (kg/m<sup>3</sup>)</td><td>10</td></tr><tr><td>OBS2</td><td>Optical backscatter sensor</td><td>5</td></tr><tr><td>RPR1</td><td>Ripple Profiler</td><td>Bed profile, z (m)</td><td>76</td></tr><tr><td>EM1</td><td>Valeport Electromagnetic Current Meter</td><td rowspan="2">Flow velocity, u,v (ms<sup>-1</sup>)</td><td rowspan="2">233.5</td><td>5</td></tr><tr><td>EM2</td><td>Valeport Electromagnetic Current Meter</td><td>10</td></tr><tr><td>PT2</td><td>Pressure transducer</td><td>Pressure, P (kPa)</td><td rowspan="4">233.5</td><td>45</td></tr><tr><td>OBS3</td><td>Optical backscatter sensor</td><td rowspan="2">Suspended sediment concentration, C (kg/m<sup>3</sup>)</td><td>10</td></tr><tr><td>OBS4</td><td>Optical backscatter sensor</td><td>5</td></tr><tr><td>RPR2</td><td>Ripple Profiler</td><td>Bed profile, z (m)</td><td>75</td></tr><tr><td>EM3</td><td>Valeport Electromagnetic Current Meter</td><td rowspan="2">Flow velocity, u,v (ms<sup>-1</sup>)</td><td rowspan="2"></td><td>11</td></tr><tr><td>EM4</td><td>Valeport Electromagnetic Current Meter</td><td>5.5</td></tr></tbody></table>

## Methods

In this section, the experimental facility and test program are described, followed by the details of the instrumentation.

**Experimental setup and morphology.** The GWK large wave flume is 309 m long, 7 m deep and 5 m wide with a combined piston-flap type wavemaker. A schematic of the experimental setup is shown in Fig. 1. All coordinates are given as the distance from the wave paddle rest position (x = 0 m), elevation above the horizontal flume bed (z = 0 m) and across-flume distance from the centreline (y = 0 m). The flume was filled with fresh water from the Mitteland canal which runs adjacent to the facility.

A large suite of instruments was deployed during the experiment and is detailed below. All instruments were logged by PCs connected to a local area network with a shared timeserver to ensure time-synchronisation. Table 1 lists all instruments and their locations within the flume, and the primary instrument positions are shown in
Fig. 1 (noting that some instruments were moved in response to water level increases and/or evolving beach morphology).

Both phases of the experiment used an initially planar sand beach with a gradient of 1:15 which was placed on top of a permanent 1:6 asphalt slope, with a minimum sand depth of 3.1 m beneath the active part of the profile (seaward of the maximum runup limit, x = 278 m). The beach was constructed using 1660 m³ of medium-coarse quartz sand (D₅₀ = 330 µm, D₉₀ = 650 µm and D₁₀ = 200 µm) from the GWK facility's material store. The sand had a density of 2650 kg/m³ and dry bulk density of 1680 kg/m³ giving a porosity of 0.37. A 25 m long layer of sand with a thickness of 0.5 m was installed in front of the slope in order to provide an additional supply of sediment. The toe of this layer was located at x = 161 m, the toe of the beach slope at x = 188.5 m and the top of the slope was at x = 283 m, z = 6.8 m (Fig. 1a).

After the first water level rise of Phase DR, a dynamic cobble berm revetment was constructed on the modified sand beach profile. The revetment was composed of 9.375 m³ (15 tonnes) of well sorted rounded granite cobbles with characteristics Dmax = 90 mm, Dmin = 50 mm, D50 = 63 mm, D85/D15 = 1.32, bulk density = 1600 kg/m³ and porosity = 0.41. The toe of the revetment was located at x = 256.8 m, z = 4.77 m, with a 1:6 slope leading to the crest at x = 260.7 m, z = 5.42 m. The overall height and width of the constructed revetment was 0.65 m and 7.3 m respectively. The revetment slope was selected based on guidance for recharge of shingle beaches¹³ and the crest elevation was designed to be at the elevation of the R₂% runup level for the second water level increment measured during Phase SB using the Lidar.

The top of the revetment extended horizontally from the crest until it intersected with the sand beach at x = 264.1 m, z = 5.42 m. Note that due to the slope of the modified sand profile approaching that of the designed revetment at the installation location, it was necessary to dig out 7.2 m³ of sand to enable the designed cobble volume to be placed (see Fig. 1).

**Test program.** The experiment was divided into two phases corresponding to sand beach (Phase SB) and dynamic revetment (Phase DR) testing. Within each phase, the profile was monitored as it evolved under wave forcing and increasing water level. Testing within each phase was undertaken at 5 different water levels (0.1 m increments), and at each water level the experiment was divided into “runs” of increasing duration as the rate of morphological change reduced (133 runs in total). An overview of the test program is provided in Table 2 and the details of all runs are listed in the dataset associated with this paper. The initial case for both phases was a 1:15 planar sand beach with a water level zwl = 4.5 m and as previously noted the beach profile was only reset between the two phases, thus all runs had a different antecedent morphology corresponding to the beach profile the end of the preceding run.

*Phase SB - Unmodified sand beach response.* Starting with an initially planar slope and a water level zwl = 4.5 m, the beach was first allowed to evolve naturally under constant wave forcing (Hs = 0.8 m,Tp = 6.0 s). The mean water level in the flume was raised by a total of 0.4 m in steps of 0.1 m. Measurements were undertaken for a period of 20 and 17 hours for the first (zwl = 4.5 m) and final (zwl = 4.9 m) water levels, and for 7 hours at the intermediate levels. This testing was divided into 63 runs with durations ranging from 20 minutes to 3 hours. Run names for this phase are given as SB < WL increment>_<<Run No.>, where water level (WL) increments are numbered 0 for the initial water level of 4.5 m to 4 for zwl = 4.9 m and run numbering is started from 1 for each WL increment.

Following the completion of the WL increments, “resilience testing” was completed to investigate the short-term response of the beach to a range of different wave conditions (“tests”) expected to produce both erosion and accretion. This testing was undertaken at the highest water level (zwl = 4.9 m). Each test was divided into 3 to 7 runs with durations ranging from 20 to 60 minutes. These runs were labelled SBE for erosive cases and SBA for cases expected to cause accretion, numbered according to test number and then run number, e.g. SBE1_3 for erosive test 1, run 3.

*Phase DR - Dynamic cobble berm revetment response.* Initially, a 1:15 planar sand beach was allowed to reshape naturally under constant wave conditions (Hs = 0.8 m, Tp = 6.0 s) for 20 hours, repeating the first WL increment of Phase SB (zWL = 4.5 m) to provide a natural beach profile on which to construct the dynamic cobble berm revetment. The cobble revetment was installed at the location of the sand beach berm according to the configuration given above. The revetment was designed such that it would be overtopped significantly as the water-level rose. The sand foreshore and dynamic revetment were then reshaped by waves (constant conditions; Hs = 0.8 m, Tp = 6.0 s) for the remaining water level increments, with the test durations at each water level mirroring those in Phase SB. Run names for this phase are given as DR < WL increment>_<<Run No.>, where WL increments and run numbers follow those for Phase SB.

After completion of the WL increments, “resilience testing” of the revetment under varying wave conditions was undertaken at the highest water level, zwl = 4.9 m. Each test was divided into 2 to 4 runs with durations ranging from 20 to 60 minutes. These runs were labelled DRE for erosive cases and DRR for cases expected to allow the revetment to recover, and numbered as per the Phase SB resilience tests.

Finally, to investigate the effect of recharging the revetment, 2.5 m³ of additional cobbles, corresponding to a 0.2 m thick layer were placed on the front face of the revetment. Following this recharge, the response of the revetment to a range of different high energy, erosive wave cases was measured. These runs were labelled DRN and numbered using the same notation as the resilience tests.

**Wave conditions.** Wave paddle steering signals were generated according to the JONSWAP spectrum (using a peak enhancement coefficient of 3.3) specified using significant wave height, Hs and peak wave period Tp. For Phases SB and DR constant wave forcing was applied, Hs = 0.8 m andTp = 6 s. This wave condition was chosen to
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Example wave measurements. (a) Timestack of water depth measured by the Lidar throughout the surf and swash zones. (b) Timeseries of water surface elevation at *x* = 225 m as indicated by the vertical dashed line in (a). (c) Measured free-surface profile through the surf and swash zone at the time indicated by the horizontal solid line in (a). Note that the measurements capture the splash-up generated by a breaking wave at *x* = 235.5 m.</figcaption></figure>

and swash zone has been successfully demonstrated¹⁷. Typically, Lidar requires bubbles to be present on the water surface to ensure that the incident laser light is scattered sufficiently to obtain a valid detection. During the experiment described here, it was found that the instruments performed better than during previous field deployments¹⁷⁻¹⁹, with valid return signals even when levels of aeration were very low or in some cases, non-existent. It is thought that this was due to the presence of fine sediment in the water column which caused light to be scattered from the water surface. Example wave data obtained using the Lidar array is shown in Fig. 2.

**Morphology measurements.** The emergent and submerged beach profile, between $x = 183$ m and $x = 270$ m was measured at the end of each run using a mechanical roller attached to the overhead trolley which ran along the centre of the flume. Figure 3a shows an example profile measurement.

A Reson SeaBat 7125 multibeam echo-sounder was deployed to obtain pilot measurements of the bubble clouds generated by wave breaking¹¹ and non-intrusive, regular measurements of the submerged beach profile. The echo-sounder was mounted on a vertical arm fixed to the overhead trolley of the mechanical profiler. The receiver was oriented in the vertical plane and aligned centrally along the length of the flume. A range of different cross-shore locations, depths and angles were tested to optimise data collection leading to a primary deployment position of $x = 223.71$, $z = 3.8$ m and an angle of 30 above the horizontal. The instrument has a 128° opening angle 0.54 beam divergence angle, operates at a frequency of 400 kHz and measurements in units of dB were collected at 1 ping per second. Note that the shallow depths and presence of bubble clouds during wave sequences make regular detection of the changing bed difficult using conventional processing methods, however new algorithms which make use of the double acoustic reflection from the water surface to the bed and back to the receiver are being developed and will be reported in future works. Due to the pilot nature of this deployment, the multiple instrument positions and orientations used, the size of the dataset and the large quantity of noisy data, the multi-beam dataset is not provided in the downloadable dataset.

Wave-by-wave measurements of the changing beach face profile were obtained using the landward-most Lidar located at $x = 255$ m. Lidar detects the uppermost surface at each scan position within the swash zone – either swash surface (when submerged) or the emergent bed (between swash events). By separating the “swash” and “bed” signals within the Lidar dataset using a variance-based approach²⁰ (see Fig. 3b) it is possible to obtain the beach profile landward of the swash rundown position between every swash event (Fig. 3c). The quoted error range for the Lidar is ±6 mm, however testing has demonstrated that for a stationary sand or cobble bed, this range is reduced to approximately ±0.95 mm.

Measurements of the entire three-dimensional bathymetry were obtained at irregular intervals when the flume was drained using a FARO Focus 3D terrestrial laser scanner. A total of 11 surveys of this type were completed throughout the duration of the experiment.

**Surf Zone/Sandbar measurements.** Two measurement rigs were installed immediately landward and seaward of the predicted sandbar location and each housed an array of instrumentation designed to measure hydrodynamics, sediment transport and morphological change during bar formation and migration. The main instrument mounting bars for these rigs were located at $x = 226.5$ and 233.5 m. Each of the measurement rigs was fixed to the walls on a mechanism such that they could be lifted and lowered manually to the bed after each run to ensure that all instruments remained a constant height above the evolving bed (see Table 1).

Each rig was equipped with the following instruments which were sampled at 8 Hz: 2 optical backscatter sensors (OBS) mounted at 5 and 10 cm from the bed, two electromagnetic current meters (EMCM) at elevations of 5 and 10 cm above the bed and a pressure transducer (PT) mounted 45 cm above the bed. The error ranges of the EMCMs and PTs are approximately ±0.015 ms⁻¹ and ±0.6 Pa respectively. Finally, a ripple profile scanner (RPS) was mounted 75 cm above the bed to obtain local bed profile measurements along a 0.9 m transect. The RPS on each rig was sampled alternately for one minute to avoid crosstalk between instruments.
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Example morphology data. (a) An example beach profile as measured by the mechanical profiler (black) and the swash zone profile obtained from the Lidar data (blue). (b) Separation of bed (black dots) and swash data at $x=253.8$ m (blue), $x=255.3$ m (red) and $x=256.8$ m (orange) for an example section of data. The mean bed elevation between each swash event is shown in white. (c) Bed elevation change relative to the initial profile in the swash zone at the wave-by-wave timescale. (d) Beach profile data showing the evolution of the sand beach and dynamic revetment modified from Bayle *et al.*<sup>6</sup>. The revetment surface is marked with a thicker line.</figcaption></figure>

In addition to the two rigs, two Nortek ADVs were located at $x=235$ and 242 m, maintained at a height 15 cm above the bed and sampled at 25 Hz. Each ADV was co-located with a pressure transducer and an additional standalone pressure transducer was installed at $x=231.7$ m, $z=4.13$ m. The error range for the ADVs for the velocities measured is approximately $\pm0.01$ ms<sup>-1</sup>.

Note that the two surf zone rigs described here were present for the entirety of Phase SB and the first 20 hours of the Phase DR testing. The instruments and scaffold rigs were removed during installation of the dynamic cobbles revetment to avoid the risk of damage due to impact from stray cobbles from the revetment. Example post-processed data from the seaward surf zone rig is presented in Fig. 4.

**Swash zone measurements.** The swash zone was monitored by a high definition IP camera (Vivotek IB9381-HT) which was used in RGB mode, the frame rate was 10 fps with a resolution of 2560 x 1920 px. The camera was mounted in the flume roof at $z=11.8$ m landward of the runup limit, facing the wave paddle. The cross-shore position of the camera varied with the water level in the range $x=267$ m to 280 m. A series of ground control points (GCPs) were positioned within the camera field of view to enable generation of rectified timestack images. The position of these GCPs was surveyed using the FARO Focus 3D terrestrial laser scanner.

The timestack images of swash flow are complimented by the data from the most landward Lidar which monitored flow depths and bed elevations within the swash zone. Separation of the “bed” and “swash” using variance criteria<sup>20</sup> as described above enables not only extraction of wave-by-wave bed elevations, but also estimates of the shoreline timeseries and depth-averaged flow velocity<sup>21</sup> and capture of the bore collapse process<sup>19</sup>. Example swash zone measurements are presented in Fig. 5.

**Instrumented cobbles.** The movement of individual cobbles within the dynamic revetment was monitored using an RFID tracking system similar to that previously used in field experiments<sup>22</sup>. The RFID system consists of three components: Passive Integrated Transponder (PIT) tags, the module reader and the antenna.
<figure><img src="image_6.png" /></figure>

**Fig. 4** Timeseries data from surf zone rig 1, *x* = 226.5 m. (a) Water depth derived from pressure transducer data, (b) cross-shore flow velocity measured 5 cm (blue) and 10 cm (red) above the bed using EMCMs, and (c) suspended sediment concentrations 5 cm (blue) and 10 cm (red) above the bed measured using OBS.

<figure><img src="image_7.png" /></figure>

**Fig. 5** Example swash data. (a) Video timestack extracted from the high definition video. (b) Timestack of water depth extracted from the Lidar data with the timeseries of shoreline position added in red.

Texas Instruments TRPGR30ATGA PIT tags with a unique identification number and a detection range of 0.6 m were installed in 97 cobbles. The tags were placed inside 5 mm diameter holes drilled into the short axis of the cobbles and sealed using epoxy glue. Following PIT installation, the cobbles were washed, dried and painted in 3 different colours: 20 cobbles were painted pink and placed on the bottom layer of the revetment (at the sand interface) during its construction; 30 cobbles were painted orange and placed 20 cm above the bottom of the revetment (mid layer); 47 cobbles were painted green and placed at the toe and on the top layer of the revetment. All cobbles were placed along the centre line of the revetment in groups of 3 cobbles at 0.4 m cross-shore intervals. An additional 7 cobbles were initially placed at the revetment toe. Finally, the crest line of the revetment was painted yellow to enable modification of the crest by waves to be easily observed (Fig. 1c). Further details of the instrumented cobble placement are provided by Bayle *et al.*⁶ and the 'DynaRev_RFID.xlsx' spreadsheet provided in the dataset associated with this paper details the initial cobble positions and locations in each RFID survey.

The RFID reader used here was a Texas Instrument Series 2000 RI-STU-251B which transmits a radio frequency of 130.2 kHz and was connected to a logging computer via an RS232 serial connection. A 120 dB beeper was used to provide an audible beep when a PIT was detected. A Texas Instrument Ri-ANT-G02E antenna was connected to the module reader. The antenna measured 20 cm by 20 cm and was attached to a telescopic pole (up to 5 m long) to allow cobble detection from the side of the flume, avoiding the need for the operator to walk on, and potentially damage the revetment. Instrumented cobble surveys were completed at the end of each water level increment and day of testing during Phase DR by passing the antenna over the revetment surface in a systematic manner. The identification number and cross-shore position of each detected cobble was recorded for each survey.

## Data Records

The data detailed in this paper is available for download from https://doi.org/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 https://doi.org/10.5281/zenodo.3860196²⁴.
<table><caption>Table 3. Data files associated with the DynaRev experiment available from https://doi.org/10.5281/zenodo.3889796²³.</caption><thead><tr><th>Filename</th><th>Data description</th><th>Instruments (ref. Table 1)</th></tr></thead><tbody><tr><td>DynaRev_TestProgram.xlsx</td><td>Complete list of test cases</td><td>—</td></tr><tr><td>DynaRev_Profiles.mat</td><td>Beach profiles measured after each run (x, z)</td><td>Mechanical profiler</td></tr><tr><td>DynaRev_Paddle_Files.zip</td><td>Wave paddle driver files in ascii format</td><td>Wave paddle</td></tr><tr><td>DynaRev_DAQ.mat</td><td>Timeseries data collected by the central data acquisition system:<ul><li>Wave gauges - surface elevation, η (m)</li><li>ADVs - flow velocity, u, v, w (ms<sup>-1</sup>)</li><li>PTs - pressure, P (kPa)</li><li>Paddle stroke (m)</li></ul></td><td>WG1 to 8, WGADV1<br>ADV1, ADV2, ADV3<br>PTADV2, PTADV3<br>Measured wave paddle stroke</td></tr><tr><td>DynaRev_SurfZone.mat</td><td>Timeseries data from the surf zone rigs:<ul><li>PTs - pressure, P (kPa)</li><li>EMCMs - flow velocity, u, v (ms<sup>-1</sup>)</li><li>OBS - sediment concentration, C (gL<sup>-1</sup>)</li></ul></td><td>PT1, PT2, PT3<br>OBS1 to OBS 4<br>EM1 to EM4</td></tr><tr><td>DynaRev_Lidar_<Phase>&lt;WL increment>&lt;Run No. &gt;.mat</td><td>Timeseries x, z data from the combined Lidar array in.mat format. The data for each run is stored in a separate file, e.g. "DynaRev_Lidar_SB1-5. mat" contains the data for Phase SB, WL 1 (z<sub>wl</sub> = 4.6 m), Run 1.</td><td>LID1, LID2, LID3</td></tr><tr><td>DynaRev_TimeStack.mat</td><td>Image timestack of swash zone</td><td>CAM</td></tr><tr><td>DynaRev_RFID.xlsx</td><td>Table containing instrumented cobble positions</td><td>RFID</td></tr><tr><td>DynaRev_3Dscans.zip</td><td>Point cloud data (x,y,z (m)) from 11 3D Lidar scans of the morphology in ".xyz" format</td><td>FARO</td></tr><tr><td>DynaRev_Lib</td><td>Scripts for post-processing raw instrument data</td><td></td></tr></tbody></table>

## Technical Validation

All data was collected using well-established coastal field and/or laboratory techniques using commercially available instrumentation. Post-processing was undertaken to remove outliers and convert spatial data to the x, y, z coordinate system defined above.

The profiler system provides the beach profile data directly in the local coordinate system (x, z). A visual check was completed directly after each profile to ensure no obvious measurement errors. Where errors were detected, the profile was repeated. The elevation data was interpolated onto a 0.025 m cross-shore grid.

The output from each Lidar provides the distance to the nearest target for every angle within each 2D scan at 25 Hz. This data was converted to local Cartesian coordinates (x, z) based on the position and orientation of each Lidar within the flume and interpolated onto a 0.1 m cross-shore grid. Outliers were only obtained where an object or person was positioned within the Lidar scan and these were removed manually. The exact location and orientation of the Lidar array was confirmed through comparison with the mechanical beach profiler data when no waves were running (see Fig. 3a). A RMSE smaller than 0.014 m was obtained.

Data from the wave gauges, ADVs, PTADV1 and PTADV2 (see Table 1) were sampled by the central GWK data acquisition system at 25 Hz. All wave gauges were calibrated at regular intervals throughout the experiment using a standard procedure. For each calibration, the water level was lowered from 5 m to 0.5 m in increments of 0.3 m and the voltage from all wave gauges at each water level was recorded for 180 s to create a calibration function relating water level to voltage. Wave gauge data was provided by the GWK system as a timeseries of water surface elevation in metres relative to the mean water level. ADV data was provided as u, v, w velocities (ms⁻¹) and the pressure data were corrected for atmospheric pressure and provided in kPa.

In the surf zone, PTs were sampled at 8 Hz, corrected for atmospheric pressure and provided in kPa. EMCM data was sampled directly as u, v velocities at 8 Hz, no further post-processing was undertaken. The time-varying free surface elevations obtained from the Lidar data were compared with point measurements from pressure transducers PT1, PT2 and PT3 and wave gauge WGADV1 (see Table 1). For all runs the signals matched closely with zero lag.

All optical backscatter sensors were calibrated after the experiment to provide sediment concentration (gL⁻¹) by applying the method of Betteridge et al.²⁵ using sand from DynaRev in a specially constructed sediment tower at the University of Plymouth.

Camera timesteps were processed by extracting a line of pixels along the flume centreline and rectified using surveyed ground control points within the camera field of view.

## Code availability

All code provided in DynaRev_Lib is written in MATLAB (R2019b). This folder contains the scripts used to process the raw data in order to obtain the post-processed data provided within the repository.

The 3D Lidar point clouds described in Table 3 are provided in ".xyz" format which can be opened using the open source CloudCompare software package. The filename for each scan includes the date collected and the run after which the scan was completed, e.g. 20170918_DR2_7. xyz was completed after Run DR2-7 on 18th September, 2017. A table providing the timings and notes about each scan is included within the DynaRev_3Dscans data record.

Received: 11 September 2020; Accepted: 14 December 2020; Published: 20 January 2021
## References

1. Almeida, L. P., Masselink, G., McCall, R. & Russell, P. E. Storm overwash of a gravel barrier: Field measurements and XBeach-G modelling. *Coast. Eng.* **120**, 22–35 (2017).

2. Briganti, R., Dodd, N., Incelli, G. & Kikkert, G. Numerical modelling of the flow and bed evolution of a single bore-driven swash event on a coarse sand beach. *Coast. Eng.* **142**, 62–76 (2018).

3. Masselink, G. et al. Large-scale Barrier Dynamics Experiment II (BARDEX II): Experimental design, instrumentation, test program, and data set. *Coast. Eng.* **113**, 3–18 (2016).

4. Eichentopf, S., van der Zanden, J., Cáceres, I., Baldock, T. E. & Alsina, J. M. Influence of storm sequencing on breaker bar and shoreline evolution in large-scale experiments. *Coast. Eng.* **157**, 103659 (2020).

5. Allan, J. C. & Komar, P. Environmentally compatible berm and artificial dune for shore protection. *Shore and Beach* **721**, 9–16 (2004).

6. Bayle, P. M. et al. Performance of a dynamic cobble berm revetment for coastal protection. *Coast. Eng.* **159**, 103712 (2020).

7. Price, T. D., Ruessink, B. G. & Castelle, B. Morphological coupling in multiple sandbar systems - a review. *Earth Surf. Dyn.* **2**, 309–321 (2014).

8. Ranasinghe, R., Callaghan, D. & Stive, M. Estimating coastal recession due to sea level rise: Beyond the Bruun rule. *Clim. Change.* **110**, 561–574 (2012).

9. Chardón-Maldonado, P., Pintado-Patiño, J. C. & Puleo, J. A. Advances in swash-zone research: Small-scale hydrodynamic and sediment transport processes. *Coast Eng.* **115**, 8–25 (2016).

10. Blenkinsopp, C. E., Turner, I. L., Masselink, G. & Russell, P. E. Swash zone sediment fluxes: Field observations. *Coast. Eng.* **58**, 28–44 (2011).

11. Bryan, O., Bayle, P. M., Blenkinsopp, C. E. & Hunter, A. J. Breaking wave imaging using lidar and sonar. *IEEE J. Oceanic Eng.* **45**, 1–11 (2019).

12. McCall, R., Rijper, R.H. & Blenkinsopp, C.E. Towards the development of a morphological model for composite sand-gravel beaches. In *Proc. 9th International Conference on Coastal Sediments* (Eds. Wang, P., Rosati, J.D. & Vallee, M.) 1889–1900 (World Scientific, 2019).

13. Powell, K.A. Dissimilar Sediments - Model tests of replenished beaches using widely graded sediments. Technical Report SR350 (HR Wallingford, 1993).

14. Dean, R.G. Heuristic model of sand transport in the surf zone, In *Proc. 1st Australian Conference on Coastal Engineering* 208-214 (Institution of Engineers, Australia, 1973).

15. Hattori, M. & Kawamata, R. Onshore-Offshore transport and beach change. *Coast. Eng.* **72**, 1175–1193 (1980).

16. Dalrymple, R.A. & Thompson W. Study of equilibrium beach profiles. In *Proc. 15th International Conference on Coastal Engineering* 1277–1296 (ASCE, 1976).

17. Martins, K. et al. High resolution monitoring of wave transformation in the surf zone using a LiDAR scanner array. *Coast. Eng.* **128**, 37–43 (2017).

18. Puleo, J. A. et al. A Comprehensive Field Study of Swash-zone Processes, Part 1: Experimental Design With Examples of Hydrodynamic and Sediment Transport Measurements. *J. Waterw. Port Coast. 140*, 14–28 (2013).

19. Bergsma, E. W. J., Blenkinsopp, C. E., Martins, K., Almar, R. & Almeida, L. P. Bore collapse and wave run-up on a sandy beach. *Cont. Shelf Res.* **174**, 132–139 (2019).

20. Almeida, L. P., Masselink, G., Russell, P. E. & Davidson, M. A. Observations of gravel beach dynamics during high energy wave conditions using a laser scanner. *Geomorph.* **228**, 15–27 (2015).

21. Blenkinsopp, C. E., Turner, I. L., Masselink, G. & Russell, P. E. Validation of volume continuity method for estimation of cross-shore swash flow velocity. *Coast. Eng.* **57**, 953–958 (2010).

22. Allan, J. C., Hart, R. & Tranquili, J. V. The use of passive integrated transponder (pit) tags to trace cobble transport in a mixed sand-and-gravel beach on the high-energy Oregon coast, USA. *Mar. Geol.* **232**, 63–86 (2006).

23. Schimmels, S., Blenkinsopp, C. DynaRev - Dynamic Coastal Protection: Resilience of Dynamic Revetments Under Sea Level Rise. *Zenodo* https://doi.org/10.5281/zenodo.3889796 (2020).

24. Schimmels, S., Blenkinsopp, C. DynaRev - Dynamic Coastal Protection: Resilience of Dynamic Revetments Under Sea Level Rise. *Zenodo* https://doi.org/10.5281/zenodo.3860196 (2020).

25. Betteridge, K. F. E., Thorne, P. D. & Cooke, R. D. Calibrating multi-frequency acoustic backscatter systems for studying near-bed suspended sediment transport processes. *Cont. Shelf. Res.* **28**, 227–235 (2007).

## Acknowledgements

This project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No 654110, HYDRALAB+. P. Bayle and G. Hennessey are supported by a PhD scholarship through the EPSRC CDT in Water Informatics: Science & Engineering (WISE). A. Hunter and O. Bryan were supported by a University of Bath Alumni Fund grant. T. Baldock acknowledges support from the Australian Research Council through Discovery grant DP140101302. R. McCall and H. Rijper acknowledge funding provided by the "Hydro- and morphodynamics during extreme events" Deltares Strategic Research Program. The authors would like to thank Aaron Barrett and Jak McCarroll for their help with decommissioning the experiment and all staff at the GWK facility for their assistance before, during and after the DynaRev experiment.

## Author contributions

C.E.B. conceived the project with input from D.C., G.M., I.L.T., T.E.B., R.M. and A.R. P.M.B. and C.E.B. did the primary experiment design and managed the experiment logistics, further experiment design input was provided by D.C., G.M., I.L.T., T.E.B. and A.H. P.M.B. led the experimental team with all authors providing assistance at various stages during the experiment including: setting up equipment, logging data, undertaking manual investigation and decommissioning. S.S. and M.K. provided the primary in-house planning and experiment assistance.

## Competing interests

The authors declare no competing interests.

## Additional information

Correspondence and requests for materials should be addressed to C.E.B.

Reprints and permissions information is available at www.nature.com/reprints.
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver http://creativecommons.org/publicdomain/zero/1.0/
applies to the metadata files associated with this article.

© The Author(s) 2021, corrected publication 2021