Video-Based Nearshore Bathymetric Inversion on a Geologically Constrained Mesotidal Beach during Storm Events

: Although geologically constrained sandy beaches are ubiquitous along wave-exposed coasts, there is still a limited understanding of their morphological response, particularly under storm conditions, which is mainly due to a critical lack of nearshore bathymetry observations. This paper examines the potential to derive bathymetries from video imagery under challenging wave conditions in order to investigate headland control on morphological beach response. For this purpose, a video-based linear depth inversion algorithm is applied to three consecutive weeks of frames collected during daylight hours from a single ﬁxed camera located at La Petite Chambre d’Amour beach (Anglet, SW France). Video-derived bathymetries are compared against in situ topo-bathymetric surveys carried out at the beginning and end of the ﬁeld experiment in order to assess the performance of the bathymetric estimates. The results show that the rates of accretion/erosion within the surf zone are strongly inﬂuenced by the headland, whereas the beach morphological response can be classiﬁed into three main regimes depending on the angle of wave incidence θ p : (1) under deﬂection conﬁguration ( θ p > 0°), the alongshore sediment transport was trapped at the updrift side of the headland, promoting sand accretion. (2) Under shadowed conﬁguration ( θ p < 0°), the interruption of the longshore current drove a deﬁcit of sand supply at the downdrift side of the headland, leading to an overall erosion in the surf zone. (3) Under shore-normal conﬁguration ( θ p = 0°), rip channels developed, and up-state beach transition was observed. A comparison between video-derived bathymetries and surveys shows an overall root mean square error (RMSE) around 0.49 to 0.57 m with a bias ranging between − 0.36 and − 0.29 m. The results show that video-derived bathymetries can provide new insight into the morphological change driven by storm events. The combination of such inferred bathymetry with video-derived surface current data is discussed, showing great potential to address the coupled morphodynamics system under time-varying wave conditions.


Introduction
Morphological changes of the beach profile are spread over a wide range of spatiotemporal scales that vary from a few hours (storm) to several weeks (changes in swell regimes) at the scale from ripples to sandbars [1]. Such morphological changes are driven by the interactions between flow (waves and currents), sediment transport and the evolving morphology that feeds back into the hydrodynamics, forming the coupled morphodynamic system. Decisions on coastal zone management are mainly based on understanding sediment budgets, which include the subaerial and subtidal zones [2]. Critical to beach response and sediment budget is the nearshore bathymetry, which is the primary source of morphological variability, but it is challenging and expensive to survey in situ. In situ bathymetric data are traditionally obtained using an echo sounder mounted behind a jet ski [3] or a small vessel [4] equipped with a Real-Time Kinematic Differential Global Positioning System sensor (RTK-DGPS) that can retrieve depths within a couple of centimeters of accuracy under calm wave conditions [5]. At best, such a monitoring approach can be performed twice a month on beaches with common periods of low wave conditions and access to the sea for small vessels [6]. Consequently, it is of great interest to find more accessible alternatives to obtain bathymetric data with high spatial and temporal resolution.
A wide range of remote sensing approaches have been proposed to infer subtidal bathymetry by exploiting different signatures associated with water depth. Among these groups of methods are the one based on light penetration in the water column and the subsequent reflection of the seabed (e.g., multi-and hyperspectral depth inversion methods [7][8][9][10]) and the one based on sea surface characteristics. The latter group can broadly be subdivided into methods that rely on wave dissipation patterns caused by depthinduced wave breaking [11][12][13][14][15] and changes in local wave celerity over a varying depth profile [16][17][18][19][20][21]. Much effort has been dedicated to derive bathymetry maps from video imagery through the linear dispersion relationship for free surface waves [18,19,[22][23][24][25][26][27][28][29]. In particular, shore-based video systems have the capability to collect large volumes of data with high temporal and spatial resolution over long periods with the advantage of being relatively low cost [30,31].
The cBathy algorithm, developed by Holman et al. [19], is a spectral depth inversion method that is nowadays the most popular algorithm to obtain two-dimensional bathymetries from video stations [6,25,[32][33][34][35][36][37][38][39][40][41][42][43][44]. cBathy is based on the linear wave dispersion relationship, and therefore, its validity is inherently bounded to the increasing degree of wave non-linearity (finite amplitude effects) as waves approach the shore, leading to larger propagation speeds for higher waves [45,46]. Although there are several studies that have proposed the correction for the finite-amplitude dispersion in shallow water (e.g., [47][48][49]), these non-linear depth inversion approaches require additional information on wave height evolution, which is difficult to acquire remotely. In any case, the cBathy algorithm still provides a good alternative to estimate robust bathymetry maps as it incorporates data assimilation through a Kalman filter to reduce instantaneous depth errors related to the loss of optical wave signature.
The quality of bathymetric estimates derived from cBathy has previously been compared with conventional bathymetric surveys obtained from echo sounders and acoustic altimeters for a number of different beaches (summarized in Table 1). However, cBathy has rarely been tested in beach environments under storm events with waves larger than 2 m and T p >10 s due to the difficulty of conducting in situ surveys in the presence of large waves [38]. Moreover, most video-based depth inversion studies have focused on open beaches but not on geologically constrained beaches, which show more complex and less understood morphodynamic processes [50][51][52][53][54]. Geologically constrained beaches are characterized by significant geological controls due to headlands, reefs, platforms, rock outcrops and islets, which determine beach boundaries, beach morphology and morphodynamics [55,56], making them behave differently from unconstrained beaches [57][58][59][60][61][62][63].
In this study, we explore the potential of the cBathy algorithm to derive nearly continuous (every 30 min during daylight hours) bathymetry maps during a 3-week field experiment in order to examine the morphological evolution of a geologically constrained beach under a wide range of energetic wave and tide conditions. Field measurements are used to assess the capability of cBathy to reproduce underwater morphological features before and after a storm event. Overall, this work aims to extend the limited knowledge we have of cBathy performance under high-energy waves and to explore morphological changes and their driving factors at a geologically constrained beach. The combination of such video-derived nearshore bathymetry with other video-derived products to better understand the morphodynamic system is further discussed.

Study Site and Field Experiment
La Petite Chambre d'Amour, hereafter referred to as LPCA beach, is the last and southernmost beach of Anglet located on the Basque coast in the southwest of France (Figure 1). LPCA beach is bounded by a 500 m long rocky headland in the South (Saint-Martin headland) and characterized by the presence of a 20 m wide submerged rocky reef located approximately 150 m from the headland. According to Wright and Short [64] classification, Anglet is a high-energy intermediate beach with a bimodal sediment distribution with a steep beach face (tan β ≈ 1/10 and D 50 ≈ 2 mm) and a more gentle sloped surf zone (tan β ≈ 1/50 with D 50 ≈ 0.3-0.4 mm) [65] with the presence of a seasonally modulated single to double-bar system [66]. An intensive field campaign was conducted at LPCA beach from 3 to 26 October 2018 with the objective to investigate the wave-induced circulation and the morphological response of the beach under energetic wave conditions. The experiment involved the deployment and data collection of a large array of instruments including current profilers, surf-zone drifters, topo-bathymetric surveys and image acquisition from a fixed videocamera station and a camera equipped Unmmaned Aerial Vehicle (UAV).
Wave and tide data were continuously retrieved (every 30 and 60 min, respectively) from a 6 km offshore wave buoy and from the Bayonne-Boucau tide gauge, which is located within the mouth of the Adour river, 4 km north of the study site. Figure 2 shows the offshore wave and tide conditions during the field experiment. LPCA beach was exposed to Atlantic W-NW long-period swells, where significant wave heights H s ranged from 0.4 to 4.2 m and peak wave periods T p ranged from 5 to 16 s (Figure 2a). Offshore wave conditions were relatively energetic (average H s of 1.5 m and average T p of 12 s) with several high-energy wave events exceeding 2 m. During the experiment, the tidal range varied from 1.1 to 4.4 m for neap tides and spring tides, respectively.
An important parameter governing beach morphological changes is the wave energy flux P (also referred to as "wave power" [67,68]) which can be approximated [69,70] by: where ρ is the seawater density (1025 kg/m 3 ) and g is the gravity acceleration (9.81 m/s 2 ). The longshore component of the energy flux P y can be obtained as: where θ p is the offshore peak wave incidence (angle of wave incidence relative to the shore normal). P y represents the portion of wave power available to drive longshore currents in the surf zone [65,71]. Therefore, P y can be used as a proxy of sediment transport along the coast [67]. The sign of P y depends on θ p ; high positive (negative) values of P y correspond to high-energy swell that arrives obliquely to the shore to the right (left) of the shore normal, whereas near-zero values of P y indicate shore-normally incident waves. As the LPCA beach is geologically constrained by the headland in the south, three main regimes of wave-induced circulation are identified based on the angle of wave incidence [65,72,73]: (1) a headland-directed longshore current that evolves into a rip current along the headland under deflection configuration (θ p > 0°and P y > 0 kW/m), (2) an onshore-directed current and a weak oscillating eddy induced under a shadowed configuration (θ p < 0°and P y < 0 kW/m), and (3) cross-shore motions with rip cells' circulation under a shore-normal configuration (θ p ≈ 0°and P y ≈ 0 kW/m). Figure 2b shows the filtered values of the offshore peak wave incidence (−13 < θ p < 22°) and the longshore component of wave power (−7 < P y < 21 kW/m). θ p and P y time series were smoothed using a 12 h window moving average.

Topo-Bathymetric Surveys
During October 2018, several surveys were conducted with the purpose to quantify intertidal and subtidal morphological variations throughout the experiment (see Figure 2). Two mono-beam bathymetric surveys with 50 m spaced cross-shore transects were carried out at the beginning (5/Oct/2018) and end (26/Oct/2018) of the field experiment. All bathymetries were obtained through Real-Time Kinematic Global Navigation Satellite System (RTK-GNSS) surveyed from a boat. The upper and intertidal beach topography were surveyed using a Post Processing Kinematic Differential Global Navigation Satellite System (PPK-GNSS) carried in on foot on 15 October 2018. In addition, a high-resolution topographic survey employing photogrammetry [74] obtained by a UAV was performed at spring low tide on 24 October 2018 to measure the headland topography as well as the upper shore-face topography. The vertical uncertainty (95% CI) from bathymetric and topographic surveys is estimated to be less than 0.20 and 0.07 m, respectively [65].
It was decided not to merge the bathymetry data collected on 5 October with the topography collected on 15 October, as a high-energy wave event occurred between the two dates (H s up to 4 m on 7 October 2018; see Figure 2), potentially inducing changes in the seabed morphology. On the other hand, the topographic and bathymetric data collected on 24 and 26 October 2018, respectively, were combined and gridded using a natural neighbor interpolation approach accounting for anisotropy to produce a single bathy-topo map. Thus, only the bathymetry data (5/Oct/2018), as well as the topo-bathymetric ensemble (24,26/Oct/2018) representing the morphological conditions during the beginning and end of the experiment, respectively, will be used as ground truth to further validate cBathy estimations.

Video Data
Images capturing LPCA beach were collected from a single camera installed on top of the Biarritz lighthouse (70 m above mean sea level) located near the tip of Saint-Martin headland ( Figure 1). Then, 1 Hz sampled images recorded at 1624 × 1234 px were continuously collected during daylight hours on 5,8-15,18-30/Oct/2018 ( Figure 2). Images containing fog, sun glare and raindrops were visually identified and manually excluded. The remaining images were stabilized [75] and rectified [76] into local coordinates using 16 ground control points (GCPs) spatially distributed both on land and water.

cBathy Algorithm
The cBathy algorithm (Holman et al. [19]; https://github.com/Coastal-Imaging-Research-Network/cBathy-toolbox [44], accessed on 20 January 2020) is based on the linear wave dispersion relationship to estimate depth and thus requires a time series of images with the presence of waves in intermediate or shallow water depths. Images are typically sampled at 1 or 2 Hz over a 17 min collection to yield a burst or stack of images. Data runs are typically collected hourly but can be reduced to half-hourly for field experiments or places where conditions, such as tide or morphology, change rapidly. The image spatial resolution should be adequate to resolve the anticipated dominant wavelengths. For the case of LPCA beach, the wavelength for a typical 10 s wave in 1 m depth is ∼31 m. This means that for a cross-shore spacing (∆x p ) of 5 m, six points will be provided to resolve the wavelength, which is fine to keep the computation time reasonable. As alongshore scales are longer, they can be well resolved with a larger sample spacing (e.g., ∆y p = 10 m). Once the stack of images is recorded, the analysis is carried out sequentially at a series of user-selected analysis points (x m , y m ) for which a group of neighboring pixels is used within a user-specified range (x m ± L x , y m ± L y ) or tile ( Figure 3). The cBathy algorithm consists of three processing steps. The objective of step 1 is to estimate a wavenumber k at the point (x m , y m ) for a user-defined set of candidate frequencies f b . These frequencies f b are typically defined between 0.056 and 0.250 Hz and correspond to wave periods ranging from 4 to 18 s. Each pixel within the analysis tile (x m ± L x , y m ± L y ) is treated as an individual pixel intensity time series and converted to the frequency domain through a Fourier transformation. A cross-spectral matrix is then computed between all possible pixel pairs inside the tile for each of the desired frequency bands f b . A selection of the most coherent frequencies (commonly four; cBathy default value) are identified and extracted through spatial empirical orthogonal function (EOF) to determine the dominant spatial phase of the wave. The corresponding wavenumbers k are derived by fitting the wave patterns of the observed spatial phase to a forward modeled wave train. A skill value is provided to indicate the percentage of the variance explained by the fit, where a skill score of 1 corresponds to a perfect match. Finally, for each f b − k pair, a depthh(x m , y m ) can be estimated using the linear dispersion relationship: where ω(2π f ) is the radial frequency, f is the linear frequency, k is the wavenumber, h is the water depth and g is the gravitational acceleration. Since cBathy estimates depth, not bathymetry, tidal elevations must be subtracted from estimates to yield depths relative to a fixed tidal datum.
In step 2, a single depthĥ(x m , y m ) is computed by combining the multiple f b − k pairs estimated from step 1 through a weighted non-linear least squares fit so that the best correspondence to the linear dispersion relationship is found.
Step 2 depths, similar to those of step 1, must be tide-corrected to yield bathymetry data referenced with respect to a tidal datum.
Step 3 uses a Kalman filter [77] to smooth and average the datum-referenced depths in order to provide a robust and reliable depth estimate, h(x m , y m ). As the Kalman filter depends on both the current and the previous cBathy-estimated bathymetry, this step only works if more than one video-derived bathymetry is available. The purpose of the Kalman filter is to provide an estimate of beach morphology over certain areas even when in situ conditions are not favorable by filling gaps in coverage and objectively averaging new estimates with prior estimates.
cBathy Settings The image time series were rotated and rearranged into stacks of ∼17 min (1024 s) duration, yielding one image stack every 30 min. All images within a stack were rectified by setting the reference level (tidal elevation) according to the time of the middle image of the stack. This was completed to avoid rectifcation inaccuracies (i.e., horizontal pixel shifting) induced by large tidal level changes, as described by Bergsma et al. [36]. The cBathy algorithm was then applied to a total of 330 image stacks. The parameters used for cBathy are shown in Table 2. Time series of water depth were estimated on a 25 × 10 m analysis grid (alongshore × cross-shore spacing) and further referenced to NGF-IGN69 (Global French Levelling) to produce bathymetry maps. For future reference, the zeroelevation contour relative to NGF-IGN69 is located 0.43 m below MSL. The Kalman filter was then consecutively applied over time to yield a stable running average depth by automatically weighting better estimates from prior information. In addition to the Kalman filter, a moving average was applied with a temporal window of two (daylight) days over the whole bathymetry time series in order to remove a spurious depth-diurnal signal, which was probably related to a rectification inaccuracy. It should be noted that for the 5th of October 2018 (the same day as the first bathymetry survey), only one image stack was available, meaning that the Kalman filter could not be applied during that day.

Topo-Bathymetric Surveys Comparison
As shown in Figure 4, a comparison between the October 5 and 26 surveys indicates the formation of an alongshore uniform sandbar centered at x = 300 m with a net deposition of sand (around 1 m of accretion) at the respectively cross-shore and longshore position: x = 300 m, y = 400 m. Moreover, there is an overall erosion at x = 400-600 m as well as beach accretion around x = 100 m. The development of a shallow rip channel is also observed at x = 250 m, y = 400-500 m. Based solely on these survey observations, it is difficult to determine under what specific hydrodynamic conditions the morphology responded in relation to the time between wave events and the pre-existing geology. Hence, it highlights the importance of estimating bathymetry more frequently using remote sensing techniques, the results of which are presented in the following subsections.  However, the remaining shallow water regions, including the sandbar, are fairly reproduced in shape but overestimated in depth (inaccuracies around 1 m), as shown in Figure 5d and the cross-shore and alongshore profile transects in Figure 5e,f. Figure 6 shows the estimated Kalman-filtered cBathy bathymetry corresponding to the topo-bathymetry survey conducted at the end of the field experiment (24,26/Oct/2018). The standard deviation image (Figure 6a) depicts an oblique rip channel within the surf zone range of y = 400-500 m, which is consistent with the surveyed topo-bathy map (Figure 6c). This rip channel is resolved by cBathy estimations but underestimated in depth. cBathy shows an overall RMSE and bias of 0.49 and −0.29 m, respectively, with larger errors near the steep shoreline and reef location, as shown in Figure 6d and the cross-shore and alongshore profile transects in Figure 6e,f. cBathy performance is hypothesized to be lowered by wave breaking due to low tide and relatively large waves (H s = 0.8 m; T p = 11 s; θ p = 14°; tide = −1.5 m).

cBathy Error Assessment
Following the approach proposed by Bouvier et al. [42], a quality assessment was computed for each of the non-filtered cBathy-derived bathymetriesĥ(x m , y m ) by counting the number of points for which the cBathy algorithm returns a physical value (skill higher than 0.5 with associated depth errors lower than 1 m). The quality assessment (Qual) of video-derived bathymetries was assessed in percentage terms in relation to the total number of grid points inside the camera viewfield delimited by values below the zero-elevation contour line. In other words, the quality assessment counts how many non-NaN depth estimates are retrieved from cBathy within the sampled area. Figure 7 shows the cBathy quality assessment for the 330 computed bathymetries according to their corresponding tidal elevation and offshore significant wave heights (H s ). Tidal elevation has an influence on the quality of depth inversion with appreciable errors associated with low tide stage. During low tide, waves break over the reef and sandbar affecting the readability of the wave signature and resulting depth estimates. In addition, non-linear effects affect the estimation of the bathymetry in the outer and inner surf zone as the depth range decreases at low tide. Surprisingly, large waves (H s > 2 m) do not appear to significantly reduce cBathy estimations in the sampled area as opposed to the results of Bouvier et al. [42]. Nevertheless, it is important to note that this quality assessment is based only on the total number of individual depth outputs within the grid, whose performance at the end is given by the same cBathy algorithm. Figure 8 shows the LPCA beach morphological up-state transition from a Low-Tide Terrace (LTT) to a Transverse Bar and Rip (TBR) beach state according to Wright and Short [64] classification. Table 3 shows the field conditions corresponding to    Figure 2). Despite the storm, the sandy bed morphology remained reasonably uniform alongshore, showing a subdued terrace bar still attached to the shore. Figure 8b,g,l shows the beach morphology following the second energetic wave event (14/Oct/2018; H s ≈ 3 m) that occurred under a shadowed configuration. Overall, the large waves eroded the surf zone, yielding little more than a pronounced terrace bar. On 22 October 2018, shorenormally incident waves led to a state transition in between LTT and TBR (e.g., [78]) with the formation of rip channels that did not incise the outer edge bar but increased the trough continuity (see Figure 8c,h,m).

cBathy Video-Derived Morphological Evolution
After 23 October 2018, the angle of wave incidence augmented under a deflection configuration, inducing a longshore current towards the headland. This period was characterized by persistent accretion, as the sediment supply from the longshore drift was probably trapped by the headland. During the high-angle wave event, the rip channel orientation became more oblique, and the successive infilling of the trough caused part of the sandbar to weld to the shore (see Figure 8d,i,n). The latter description coincides with an erosive Transverse Bar and Rip (eTBR) state, as proposed by Price and Ruessink [71] who related the offshore dominant oblique angle of wave incidence with a discontinuous trough and the presence of oblique rip channels.
By the end of the experiment, LPCA beach was characterized by a TBR state with the presence of two deep skewed rip channels (Figure 8e,j,o). Interestingly, the relatively high-energy wave event of 28 October 2018 (H s = 3 m; θ p = 12°; P y = 7 kW/m) resulted in accretion rather than erosion; an overall sand deposition along the bar and the region around the reef is observed in accordance with survey measurements from Figure 4c.

cBathy Video-Derived Profile Response
Different cross-shore transects were selected (Figure 9a-c) to investigate the temporal evolution of the seabed over the 3-week field experiment. Wave and tide conditions are presented in order to associate morphological evolution to hydrodynamic conditions at LPCA beach (Figure 9d). The 2-day smoothed cross-shore profile time-evolution shows progressive sand accumulation along the bar (Figure 9f,g) and the reef (Figure 9e) linked by the increasingly obliquely incident wave conditions after 23 October. Sand deposition is also observed during the energetic wave event of 28 October 2018, which is consistent with the sand accretion observed for the same region in the topo-bathymetric survey conducted at the end of the experiment (Figure 4c). The temporal evolution of the smoothed cross-shore profile transect centered at y = 700 m (Figure 9g) illustrates the temporary dissection of the low tide terrace leading to a rip channel formation on 20-22/Oct/2018, which coincides with the transition to shore-normal incident waves. Overall, the cBathy-derived cross-shore transects (Figure 9h Similarly, Figure 10 shows the LPCA beach 2-day time-averaged vertical profile response over time for three different alongshore transects depicted in Figure 10a-c. A systematic change in seabed elevation is evident after 23 October 2018 with more sand accumulation close to the reef (y ≈ 300 m; Figure 10f,g,i,j). This decrease in depth appears to be related to the transition of the angle of wave incidence from a shadowed configuration to a deflection configuration (see θ p and P y in Figure 10d). The temporal evolution of the alongshore transect corresponding to the shallowest subtidal region (Figure 10e,h) shows an overall bias with respect to the ground truth surveys that might be related to non-linear effects resulting in depth overestimation.

cBathy Performance and Sources of Errors
Video-derived bathymetries were obtained under challenging wave conditions. Nevertheless, the comparison between surveys and cBathy estimates revealed typical errors (RMSE = 0.49 to 0.57 m and bias = −0.36 to −0.29 m) similar to previous works (see Table 1). The very upper part of the beach face was related with less accurate results, as analysis tiles may have contained partly wet/dry pixels at the sea-land interface, mixing unuseful subaerial and acceptable subaqueous signals [54]. The systematic shallow water bias (overestimation of true depth) can also be explained by finite amplitude effects that increase wave celerity as waves shoal and break, reducing the validity of the linear dispersion relationship [38,46].
The current cBathy version does not include the effects of currents and Doppler shifting in the dispersion relation (Equation (3)). For the case of LPCA beach, persistent rip currents (e.g., near the headland and within the rip channels) can potentially interact with incident short-period waves, shortening their wavelength (k increases) and thus underestimating local depth. However, with cBathy typical settings (Table 2), short-period waves and rapid cross-shore depth changes over the sample domain cannot be resolved unless a denser pixel spacing is defined and more (higher) frequencies are accounted for analysis, thus requiring more computational effort [19].
Several studies have reported a significant variation in cBathy performance with increasing wave height [19,36,38,39,42] since breaking waves are spread over a wider surf zone, obscuring the optical wave signal and leading to incorrect estimates of f and k from imagery. Moreover, during storm conditions, the Kalman filter is not designed to account for the systematically biased errors that occur on the timescale as real bathymetric change [38], thus resulting in greater faith deterioration of the previous video-derived bathymetry estimate under energetic conditions compared to calm conditions [39]. Despite this issue, the quality assessment computed using the approach of Bouvier et al. [42] indicates the capability of cBathy to return depth estimates for more than 70% of the domain even under waves larger than 2 m.
As previously suggested by Brodie et al. [38], we acknowledge the fact that the cBathy algorithm should preferably not be used during storm events, as video-derived bathymetries may not be quantitatively reliable. Nevertheless, we believe that cBathy is still capable of providing qualitative insight into the morphological evolution of the seabed in between storms despite challenging wave conditions when beach measurements are often intermittent. This is supported based on the timex images computed throughout the field experiment, in which cBathy was able to reproduce the morphological features fairly well and match them qualitatively with the surveys.

LPCA Beach Morphological Response to Changes in Wave Direction
The 3-week video-derived bathymetry evolution showed that the morphological variability at LPCA beach was essentially influenced by the wave direction regime and the headland control. In contrast to expectations, the high-energy wave events did not necessarily correspond to abrupt erosional state transitions, as evidenced by the wave event of 28 October 2018. Remarkably, a previous study at LPCA beach [66] has even reported steadily onshore migration of the inner bar despite the extreme storms (H s > 5 m) during the outstanding winter of 2013/2014. The latter highlights the complexity of beach morphodynamics in such environments that are geologically controlled, as outlined by different authors [50][51][52][53]. In the case of LPCA beach, it has been demonstrated that the submerged reef and the headland play an important role in determining the wave-induced circulation patterns within the surf zone depending on offshore wave obliquity [65,72,73,79,80]. Obliquely incident waves approaching from the right side of the headland (θ p > 0°) induce a strong longshore current that transports sediment along the surf zone (P y > 0 kW/m) toward the headland that eventually end up being retained by this geological boundary. Although the sediment is assumed to remain accumulated on the updrift side of the headland, strong wave-induced deflection rips can lead to sand bypassing toward the adjacent beach [80]. On the other hand, when obliquely incident waves approach from the left side of the headland (θ p < 0°), the lee of the headland remains protected from these incident waves, but the longshore sediment transport is inhibited. This deficit of sediment supply is the cause of downdrift erosion within the shadowed zone (P y < 0 kW/m). Direct exposition of shore-normally incident waves (θ p ≈ 0°and P y ≈ 0 kW/m) during low tide are associated with rip current activity and rip cell circulations [65,73], resulting in rip channel formation and morphological discontinuities in the surf zone, as shown in Figure 11. In addition to the longshore drift blocked by the headland, it is hypothesized that onshore sediment transport may have also contributed to the formation of the sandbar throughout the experiment. The field experiment was characterized by long period swell and moderately energetic constructive waves with low steepness that could promote slow sandbar building. This would also explain the overall erosion in the cross-shore distance between 400 and 600 m by the end of the experiment, as evidenced in Figure 4.

Perspectives and Future Challenges
The classic beach state classification based on the Dean's parameter [64] and the Relative Tidal Range (RTR) parameter [81] assumes that the cross-shore and alongshore sediment exchange is unconstrained by geology [53]. Thus, traditional models are generally not directly applicable in such complex settings where geological control has a fundamental influence on their morphodynamics [52]. The use of improved morphodynamic models in conjunction with remotely sensed data could help to better understand beach evolution during storm events and assist coastal managers in effective management of geologically controlled beach environments.
Nowadays, the increasing proliferation of low-cost and flexible new platforms such as UAVs [82,83], swift cameras, CoastSnap [84] and online-streaming webcams [31,85,86] offer an attractive option to collect image products and derive measurements of the nearshore. In addition to image-processing techniques that allow estimating the nearshore bathymetry (e.g., cBathy [19,43]), recent works have opened up the possibility to remotely estimate 2D surface currents by tracking the drifting foam, left after the passage of breaking waves, from video imagery (e.g., optical flow-based algorithms [73,87,88]). Future work could be to apply both remote sensing techniques to provide fresh insight into the coupled (hydrodynamic/morphology) morphodynamic system, which is key to the validation of process-based morphodynamic models. So far, processed-based morphodynamic models have been validated based only on surveyed morphological changes and at best with concurrent local wave and flow punctual measurements [89]. However, sandy beach morphodynamics is driven by the interaction of the evolving seabed and rip-cell systems [64]. Measuring in situ the spatial structure of rip current systems requires the deployment of a large number of drifters, which is only possible during a limited duration and under low-to moderate-energy wave conditions (e.g., [90]). In contrast, new optical flow-based algorithms recently showed skill to infer surface currents in the surf zone and to address rip cell circulation [73,88]. Figure 11 shows an example of a 17-min time-averaged image combining both video-based techniques. It highlights the potential to address surf zone morphodynamic processes, i.e., the coupled morphology and hydrodynamics system, at unprecedented frequency and spatial coverage.

Conclusions
Bathymetry maps were estimated every 30 min during daylight hours for three consecutive weeks using a video-based linear depth inversion algorithm in order to examine LPCA beach morphological response under a wide range of wave and tide conditions. The smoothed bathymetry time evolution showed that LPCA beach morphology was controlled by the lateral boundary effects of the headland and was strongly dependent on the incident wave direction. The assessment was carried out by comparing concurrent video-derived bathymetries with in situ topo-bathymetric measurements collected during the beginning and end of the field experiment, as well as with timex images computed throughout the entire experiment. Comparisons between surveyed and video-derived bathymetric estimates showed a similar performance with respect to previous studies with an overall RMSE = 0.49 to 0.57 m and bias = −0.36 to −0.29 m. Overall, the morphological features were in approximately the right places but presented a significant vertical offset at the shallowest parts (reef and beach shore face) where depth was consistently overestimated mainly due to non-linear effects. The results suggest that video-derived bathymetries can provide qualitative information on beach morphology evolution in between storm events when environmental conditions are challenging, but caution should be exercised as morphological features may be quantitatively misrepresented. In conclusion, this study highlights the necessity of further modeling and combining existing video-based techniques to improve our understanding of beach morphodynamics near headlands, particularly under extreme wave conditions.  Acknowledgments: LPCA beach is a monitoring site labeled by the Service National d'Observation (SNO) Dynalit (https://www.dynalit.fr, accessed on 6 June 2022) with additional funding from the Observatoire de la Côte Nouvelle-Aquitaine (OCNA) to ensure the maintenance of the video station. The installation of the video station at the top of Biarritz lighthose was allowed by le Service des Phares et Balises. Storage resources for this study was provided by the computing facilities MCIA (Mésocentre de Calcul Intensif Aquitain) of Univ. Bordeaux and Univ. Pau et Pays de l'Adour.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: