Experimental Observation on Beach Evolution Process with Presence of Artificial Submerged Sand Bar and Reef

For observation on the influence mechanism of environmentally and aesthetically friendly artificial submerged sand bars and reefs in a process-based way, a set of experiments was conducted in a 50 m-long flume to reproduce the cross-shore beach morphodynamic process under four irregular wave conditions. The beach behavior is characterized by the scarp (indicating erosion) and the breaker bar (indicating deposition), respectively, and the scarp location can be formulated as a linear equation regarding the natural exponential of the duration time. Overall, main conclusions are: (1) the cross-shore structure of significant wave height and set-up is mainly determined by the artificial reef (AR); (2) the cross-shore distribution of wave skewness, asymmetry, and undertow (indicating shoaling and breaking) is more affected by the artificial submerged sand bar (ASB); (3) the ASB deforms and loses its sand as it attenuates incident waves, which leads to a complex sediment transport pattern; (4) the scarp retreat is related to the beach state, which can be changed by the AR and the ASB; (5) the AR, the ASB, and their combination decrease wave attack on the beach. In conclusion, this study proves positive effects of the AR and the ASB in beach protection through their process-based influences on beach behaviors and beach states for erosive waves.


Introduction
Coastal erosion is a problem at many coastal sites, caused by natural effects as well as human activities [1,2], and the remedial measures fall into soft and hard engineering. Soft engineering is the mechanical placement of material on the beach (beach nourishment) or within the nearshore area (shoreface or bar nourishment). Beach nourishment, introduced in the 1920s [3], is now a technique used world-wide to maintain the eroding beaches along sandy coasts [4][5][6][7]. The U.S. Atlantic coast has confirmed the positive feedback between coastal development and beach nourishment for decades [8], where the mean rate of shoreline change has reversed from −55 cm/year (erosion) prior to 1960 to about +5 cm/year (accretion) after 1960, despite inevitable shoreline retreat over large spatial scales around the world due to the rising sea-level [9]. Shoreface nourishment, introduced in the 1990s, acts as a wave filter such as reef and meanwhile aims at nourishing the nearshore zone with long-term processes [10][11][12]. However, the performance of shoreface nourishment is less effective than beach nourishment [5], because the nourished sand is redistributed by both cross-shore and long-shore transport. Generally, the effect and the longevity of soft nourishment vary greatly from place to place due to multiple factors in the complex coastal process. Beyond those, the long-term ecological influence [13][14][15][16] and the material borrowed site [17][18][19][20] are the current challenges at different levels for various coasts. Hence, soft nourishment is not economically or environmentally suitable for some sites and can even lead to an erosion "hot spot" in a successful nourishment [21], and proper structures are still needed for erosion control.
Meanwhile, hard engineering structures were designed to perform multiple functions among the natural coastal processes and basically comprise groynes, detached breakwaters, and reefs [1]. Detailed and complete investigations on the efficiency and the impact of natural and artificial structures in recent decades were carried out from the perspective of shoreline management [22][23][24], with important experiences and conclusions. First, positive attitudes were expressed towards natural reefs and stream deltas [22]. As an alternative method for dealing with increasing erosion pressures with limited resources, artificial reefs range from rock mounds parallel to the shore to prefabricated units with various types of structures, such as a concrete box, tunnel, dome, leg, or more complex shapes [25,26] with features of submergence and porosity. The deployment objectives in recent research have shifted from improving fisheries as a resource to rehabilitation of the marine ecosystem, improving surfing or diving conditions, recreational fishing, and coastline protection [27]. Second, some structures are most effective when coupled with other measures. For example, low crest structures (LCSs) perform well when groynes acting as artificial headlands are present and can help retain the nourished sand. Since many coasts have already built engineering structures, further design should incorporate and utilize the effect of the existing coastal structures [28][29][30][31]. In addition, modified structures such as submerged artificial headlands (SAH) and detached artificial headlands (DAH) were applied to enhance the water exchange and thus improve the water environment [32]. Third, the key problems, pointed out by Ranasinghe and Turner [33] and associated with but not limited to the use of LCSs, are the identification of environmental and structural parameters governing the mode (i.e., erosion or accretion), the magnitude (e.g., size of salient), and the response of the adjacent shoreline. The morphodynamic mode refers to the fundamental principle and the hydrodynamic and sediment mechanism, which is yet to be fully understood, and the latter emphasizes the longshore current and sediment transport.
Field measurements detailed the phenomena of hydrodynamic and morphological evolution of the nearshore area [34] with structures [35,36], beach nourishment [37][38][39][40], and shoreface nourishment [41][42][43][44][45][46], which also provided the foundation for the establishment of prediction models using numerical simulation [47][48][49][50][51][52]. Basin experimental models were scaled down from the prototype to test the influence of engineering scenarios on both cross-shore and long-shore processes [53][54][55][56]. For the cross-shore responses to extreme surge or wave conditions, large-scale flume experiment projects [57][58][59][60][61] were carried out to elucidate the cross-shore (shoaling, surf, and swash zone) physical processes of the morphodynamic system. Furthermore, the cross-shore hydrodynamic forces and the resulting sediment transport mechanism have been continuously investigated by fixed-bed and movable-bed flume experiments for typical profiles of dune beaches, beaches with bars and reef flat, and barriers. Achievements were made in short wave and infragravity wave dynamics, wave skewness, undertow, boundary-layer dynamics, and the relation with sediment transport [62][63][64][65][66][67] under a range of wave conditions (regular, bichromatic, and random waves). For typical fine beach sand (e.g., 0.1-0.5 mm), using natural sand scaling by geometric ratio in small wave flume may lead to model sand moving as silt, not sand. Hence, a series of light-weight model (LWM) experiments [68][69][70][71][72] were designed using light-weight model sand to cover the modeling of a fine sand beach profile, especially a large cross-shore area in a limited flume length. In spite of all the extensive research, it is still difficult to make a direct estimation of the nearshore influence of multiple measures. For instance, even for an offshore structure or a sand berm, hydrodynamic and morphodynamic changes are delicate for its cross-shore location, not to mention the increasingly complex coastal situation with integrated engineering. Besides, some questions remain to be answered: is there any difference between the actions of soft and hard interventions on the erosion features?; how do they differ from each other in a morphodynamic system?; when a combination of measures is used, which of them can be taken as having the primary impact on the hydrodynamic and sediment processes in shoaling, surfing, and swash zones?
To answer these questions, further research is needed to investigate the influence regime of soft and hard interventions on the morphodynamic mechanism. From these perspectives, we designed a set of movable bed experiments on the basis of the Beidaihe nourishment programs to investigate the cross-shore morphodynamics in this study. Beidaihe in Hebei Province, China is a famous national tourist destination and a demonstration site of beach nourishment, where submerged breakwaters (in the form of prefabricated sandbags and artificial reefs), artificial headlands (made of rippled-rock), submerged sand bars, grouted rubble revetments, and groynes were built. These offshore interventions are important elements in the whole beach restoration design. Their proper adoption may help to reduce the erosion by large waves without further sand filling and yet barely affect the water environment and aesthetic aspects.
The novelties of this process-based study include: (1) two ecological offshore interventions, i.e., artificial submerged sand bar (ASB) and artificial reef (AR), were firstly experimentally investigated for their combination influences on cross-shore hydrodynamic and morphological process; (2) the comparison was made of morphodynamic responses to individual ASB, AR, and the combination of them; (3) quantitative relations of scarp and breaker bar (indicating beach erosion and deposition process, respectively) were established, which illustrated the influence mechanism of the offshore interventions; (4) the changed wave breaking and morphological features were found to be linked with overall morphodynamic state, which implied the role of offshore interventions in beach system.

Method
For the prototype area Beidaihe of the experiment, the tides are regular diurnal with a multi-annual average tidal range of 0.74 m. Most of the time, beaches are exposed to a low-energy wave climate, i.e., about 91.23% of the recorded significant wave heights (Hs) are less than 0.6 m and only 1.01% of Hs over 0.9 m, and the average peak period (Tp) is 4.5 s. However, the relatively stable beach profiles still suffered summer storm waves attack with wave height over 2.0 m, and some "hot spots" (severe erosion locations) were found. The movable bed experiment ( Figure 1) was conducted in a wave flume (50 m long, 0.8 m wide, and 1.2 m deep) at the Laboratory of Hydraulic and Harbour Engineering, Tongji University to simulate the cross-shore morphological evolution of fine-sand beach profiles under complex hydrodynamic processes in surf and swash zones. As a grain size range of 0.1-0.5 mm was involved, a light-weight model was adopted to reduce the scale effects due to limitations of the experimental facilities. Only half of the flume width was used, causing the width of bed model to be 0.40 m, which enabled the long-shore uniformity of morphological changes. Four model wave conditions (J1-J4) were generated by a piston-type wave generator with active reflection compensation based on the field measurements, as detailed in Table 1. To answer these questions, further research is needed to investigate the influence regime of soft and hard interventions on the morphodynamic mechanism. From these perspectives, we designed a set of movable bed experiments on the basis of the Beidaihe nourishment programs to investigate the cross-shore morphodynamics in this study. Beidaihe in Hebei Province, China is a famous national tourist destination and a demonstration site of beach nourishment, where submerged breakwaters (in the form of prefabricated sandbags and artificial reefs), artificial headlands (made of rippledrock), submerged sand bars, grouted rubble revetments, and groynes were built. These offshore interventions are important elements in the whole beach restoration design. Their proper adoption may help to reduce the erosion by large waves without further sand filling and yet barely affect the water environment and aesthetic aspects.
The novelties of this process-based study include: (1) two ecological offshore interventions, i.e., artificial submerged sand bar (ASB) and artificial reef (AR), were firstly experimentally investigated for their combination influences on cross-shore hydrodynamic and morphological process; (2) the comparison was made of morphodynamic responses to individual ASB, AR, and the combination of them; (3) quantitative relations of scarp and breaker bar (indicating beach erosion and deposition process, respectively) were established, which illustrated the influence mechanism of the offshore interventions; (4) the changed wave breaking and morphological features were found to be linked with overall morphodynamic state, which implied the role of offshore interventions in beach system.

Method
For the prototype area Beidaihe of the experiment, the tides are regular diurnal with a multiannual average tidal range of 0.74 m. Most of the time, beaches are exposed to a low-energy wave climate, i.e., about 91.23% of the recorded significant wave heights (Hs) are less than 0.6 m and only 1.01% of Hs over 0.9 m, and the average peak period (Tp) is 4.5 s. However, the relatively stable beach profiles still suffered summer storm waves attack with wave height over 2.0 m, and some "hot spots" (severe erosion locations) were found. The movable bed experiment ( Figure 1) was conducted in a wave flume (50 m long, 0.8 m wide, and 1.2 m deep) at the Laboratory of Hydraulic and Harbour Engineering, Tongji University to simulate the cross-shore morphological evolution of fine-sand beach profiles under complex hydrodynamic processes in surf and swash zones. As a grain size range of 0.1-0.5 mm was involved, a light-weight model was adopted to reduce the scale effects due to limitations of the experimental facilities. Only half of the flume width was used, causing the width of bed model to be 0.40 m, which enabled the long-shore uniformity of morphological changes. Four model wave conditions (J1-J4) were generated by a piston-type wave generator with active reflection compensation based on the field measurements, as detailed in Table 1.

Experimental Design and Instrumentation
The adoption of light-weight sediment particles with a lower density yet coarser grain size than natural sand [73] means that the sediment grain size is not scaled by the geometric scale yet leads to basically the same incipient motion and settling process. The relaxation on the similitude of the geometric particles indicated the focus on the similarity of the evolution rate, rather than its duration, to the final morphology [74].
In this study, the beach profile model was made of resin sand to reproduce the evolution feature. The density of the resin sand was 1.40-1.45 g/cm 3 with grain size ranging from 0.15 to 0.18 mm, which maintained the similitude of the hydrodynamic and sediment process by satisfying the Froude number, the Shields number, and the Rouse number under the geometric scale factor of 10 (prototype/model). The details of the scaling calculations can be found in [75].
Considering the capacity and the precision of the wave maker under the condition of model scaling, four sets of characteristic wave height and period were selected, which covered normal, moderate, and storm conditions. Four irregular model waves were generated based on these parameters shown in Table 1, and the control point was located at the front side of AR with water depth of 0.47 m. The maximum values of measured profile slopes in the field were adopted, i.e., the beach slope just after the nourishment was 1:10, and the average offshore slope was 1:200. As shown in Figure 2, four cross-shore beach profile types were investigated: (a) a nourished beach with a slope of 1:10 and a mild shore-face part with a slope of 1:200 (profile type named as B-N); (b) a beach with an AR installed 20 m offshore the beach (profile type named as B-AR); (c) a beach with an ASB placed 10 m offshore (profile type named as B-ASB); and (d) a beach with both an AR and an ASB located as in B-AR and B-ASB (profile type named as B-ASB-AR). The movable bed was obstructed by a fixed end offshore with slope of 1:10. A total of 16 experimental tests were conducted with the same duration of 90 min, as listed in Table 1. The model AR was 1.80 m long and 0.30 m high and consisted of polymethyl methacrylate blocks. The ASB was an isosceles trapezoid shape with the crest width of 3.00 m and depth of 0.20 m and slopes of 1:10. The coordinate system for comparison analysis is shown in Figure 2a.
Hydrodynamic measurements were made by 13 pre-calibrated capacitive-type wave gauges (with accuracy of ±2 mm) and three acoustic doppler velocimeters (ADVLab, Nortek, with an accuracy of 0.5% of measurement ±1 mm/s) in a fixed cross-shore arrangement as in Table 2. Among all the tests, the most significant profile change was found in the beach and ASB area during the evolution process captured by videos. The final profiles were obtained from photos, which were calibrated by rulers fixed on the side-wall of the flume. The geometric accuracy was estimated to be 3 mm for video images and 2 mm for photos.

Data Analysis
The water levels and velocity components directly measured by wave gauge and ADV provided the transformation of wave energy and wave shape. Then, the characteristics of the bed profiles were achieved with a local dynamic changing process. The undertow and the total sediment transport were estimated from wave and topography measurements.

Hydrodynamics
The mean water level at each wave gauge indicated the wave set-up and set-down for positive and negative values. Spectral analysis (using Welch's approach) was conducted to calculate the significant wave height using 10 min-long data excluding a start period of 60 s.
Waves contribute to onshore sediment transport by the non-linear wave shape of the incident waves and induce streaming in the wave boundary layer, while offshore transport is related to return flow and phase coupling between wave groups and accompanying (bound) long waves [76]. In nearshore areas where the intra-wave kinematics are crucial, these mechanisms are described by

Data Analysis
The water levels and velocity components directly measured by wave gauge and ADV provided the transformation of wave energy and wave shape. Then, the characteristics of the bed profiles were achieved with a local dynamic changing process. The undertow and the total sediment transport were estimated from wave and topography measurements.

Hydrodynamics
The mean water level at each wave gauge indicated the wave set-up and set-down for positive and negative values. Spectral analysis (using Welch's approach) was conducted to calculate the significant wave height using 10 min-long data excluding a start period of 60 s. Waves contribute to onshore sediment transport by the non-linear wave shape of the incident waves and induce streaming in the wave boundary layer, while offshore transport is related to return flow and phase coupling between wave groups and accompanying (bound) long waves [76]. In nearshore areas where the intra-wave kinematics are crucial, these mechanisms are described by parameterization, as the detailed processes are difficult to isolate from each other. In this study, these parameters were compared for different profile types to determine the influences induced by the AR or the ASB.
In addition to the significant wave height Hs, a representative long wave height H s_long for wave energy in the low-frequency band (0.004-0.05 Hz, 20-250 s) [77] was calculated using the same method as for H s with the frequency domain scaled down to the model scale, corresponding to 0.01-0.15 Hz (6.67-100 s).
where η is the free surface elevation, η is the wave set-up, and − is the time averaging operator; H is the Hilbert transform. A larger value of Sk indicates larger bed shear stresses at the wave crest than at the wave trough, which indicates the tendency of the wave crest to carry the suspended sediment onshore. Abrupt accelerations (Asy) are associated with thinner boundary layers and thus higher bed shear stresses for a given velocity magnitude [81]. The undertow u (offshore-directed depth-averaged current forced by the wave energy dissipation) was estimated as the balance due to a shoreward mass flux including wave-induced (Q w ) and roller-induced (Q r ) mass flux [80]: where h 0 is the water depth at rest, c is the wave speed, ρ r is the density of the fluid in the aerated roller taken as a constant 650 kg/m 3 , A r is the roller area, and T p is the peak period. Equation (4) is used to estimate Q w based on the assumption that the orbital velocity does not vary vertically. The roller area is given by A r = 0.9H 2 , where H is the local wave height [82]. The growth of u is hereinafter referred to as the growth in its absolute value. The raw data of velocity components should be pre-processed before any calculation because of the measurement principles. The velocity data should be excluded when the correlation (COR) of the ADV data is less than 70% or the SNR (Signal Noise Ratio) is less than 20 dB. Then, the spikes were detected using the 3D phase-space thresholding method [83] and were replaced by cubic interpolation in this study. Based on the modified velocity series, the Shields number θ was calculated as follows: where ρ s is the density of the sediment particle, and ρ is the water density, g is the gravitational acceleration, D is the sediment grain size, U w is the bottom orbital velocity of wave, and f w is the wave friction coefficient calculated by the method of Soulsby [84], covering the identification of a laminar, smooth turbulent or rough turbulent regime.
f w = B(R w ) −N for laminar and smooth turbulent regime, where R w = U w A ν is the wave Reynolds number; B = 2, N = 0.16 for the laminar flow and B = 0.0521, N = 0.187 for the smooth turbulent flow; A = U w T/2π is the semi-orbital excursion of wave orbital motion at the edge of the boundary layer, and T p is used for T here. k s is the Nikuradse equivalent sand grain roughness, being 2.5D 50 . The Shields numbers at wave gauges calculated by Soulsby method can be verified by the results calculated from ADV measurement.

Profile Morphology
As the incident waves propagated onshore, the beach face was eroded, showing the vertical cut, usually called the wave-attack scarp. Meanwhile, the eroded sand was transported offshore and formed an underwater bar named a breaker bar. Therefore, the beach profile analysis focused on the development of the scarp and the breaker bar to investigate the influences of different nourishment strategies.
The top and the bottom positions of beach scarp were calculated from minimum and maximum values of the second derivative of the measured profiles [85]. In this study, the larger horizontal coordinate (sx) of the two limits was chosen as the final scarp location for analysis. The height of the breaker bar was defined as the maximum deposition thickness z bar (the maximum elevation relative to the initial bed level) with its location at the horizontal coordinate of x bar . All above location indicators are shown in Figure 3. motion at the edge of the boundary layer, and Tp is used for T here. ks is the Nikuradse equivalent sand grain roughness, being 2.5D50. The Shields numbers at wave gauges calculated by Soulsby method can be verified by the results calculated from ADV measurement.

Profile Morphology
As the incident waves propagated onshore, the beach face was eroded, showing the vertical cut, usually called the wave-attack scarp. Meanwhile, the eroded sand was transported offshore and formed an underwater bar named a breaker bar. Therefore, the beach profile analysis focused on the development of the scarp and the breaker bar to investigate the influences of different nourishment strategies.
The top and the bottom positions of beach scarp were calculated from minimum and maximum values of the second derivative of the measured profiles [85]. In this study, the larger horizontal coordinate (sx) of the two limits was chosen as the final scarp location for analysis. The height of the breaker bar was defined as the maximum deposition thickness zbar (the maximum elevation relative to the initial bed level) with its location at the horizontal coordinate of xbar. All above location indicators are shown in Figure 3.
In addition, some sand ripple appeared on the fore-slope and the crest of the ASB, with a maximum length and height of 5.80 cm and 0.60 cm, respectively, which barely affected the largescale morphological changes on the beach.

Sediment Transport
The total sediment transport rate ̅ was estimated using the sediment mass conservation equation based on the difference between captured beach profiles as in Equation (9): where qs is the instantaneous sediment transport rate, z is the bottom elevation at the cross-shore position x, and p is the sediment bed porosity, assumed to be loosely packed and homogeneous along the beach profile (p = 0.4 used here). The total transport rate at x = 0 is set to be zero. Therefore, the ̅ at a given cross-shore location x during a wave climate duration ∆ can be obtained by the discrete In addition, some sand ripple appeared on the fore-slope and the crest of the ASB, with a maximum length and height of 5.80 cm and 0.60 cm, respectively, which barely affected the large-scale morphological changes on the beach.

Sediment Transport
The total sediment transport rate q s was estimated using the sediment mass conservation equation based on the difference between captured beach profiles as in Equation (9): where q s is the instantaneous sediment transport rate, z is the bottom elevation at the cross-shore position x, and p is the sediment bed porosity, assumed to be loosely packed and homogeneous along the beach profile (p = 0.4 used here). The total transport rate at x = 0 is set to be zero. Therefore, the q s at a given cross-shore location x during a wave climate duration ∆t can be obtained by the discrete version of Equation (10) with the measured bed level changes (∆z b ): where q s > 0 (q s < 0) corresponds to the onshore (offshore) sediment transport.

Results
In this section, the experimental results are presented in terms of different profile types. First, the hydrodynamic patterns along the profiles were analyzed by parameterization of the wave processes. Then, the profile evolution comprising the scarp and the breaker bar was unified by model establishment to enable comparison among the four profile types. Finally, the effect of the ASB was studied on the deformation process, with analysis of the influence of the AR offshore.

Hydrodynamics
The cross-shore changes of hydrodynamic and sediment transport are illustrated in Figure 4. The complex coastal processes were analyzed by parameterization in terms of wave energy (significant wave height H s , long wave height H s_long ), mean water level (η), mass flux (wave skewness Sk, wave asymmetry Asy, and undertow u) and the total sediment transport rate q s .  (f4)) changed in the same trend as those in B-ASB but with smaller absolute values. This was because the incident waves were attenuated by the AR before passing to the ASB, and bottom friction and wave breaking on the ASB further dissipated the wave energy before attack on the beach slope. Generally, the cross-shore sediment transport coupled with wave breaking varied for different profile types. In B-N and B-AR, sediment transport was confined on the beach, where the upper beach was eroded, and sediment transported offshore until settlement at the lower area. Comparing Figure  4 (g1) and (g2), the maximum ̅ was diminished by the AR (by 44%-64%), with the location moving shoreward under larger wave climates. In B-ASB and B-ASB-AR, the profile change could cover nearly the whole horizontal length, because the ASB itself covered about one-third of the whole profile and acted as a sediment source. The offshore sediment transport was generated on the beach and the fore-slope of the ASB, and onshore sediment barely existed in B-N and B-AR tests. Special cases were B-ASB-J4 and B-ASB-AR-J3, where onshore sediment transport dominated the area between the beach and the ASB.

Evolution process
Beach nourishment often causes disturbances, e.g., beach scarps and dunes, that eventually As Figure 4a1 shows, H s in B-N decreased as the incident waves propagated onshore, despite slight local increases at W2, W7, and W9. The wave energy of the low frequency domain (Figure 4b1) decreased along the propagation, and a large incident wave generated large H s-long at each wave gauge. The mean water level η (Figure 4c1) peaked at W11 for all incident waves, while there was another significant growth at W3-W5 in B-N-J4. Hence, in addition to the main shoaling zone on the beach slope, there were two local areas with a relatively strong shoaling effect for J4 waves and one for J1-J3. This was also verified by Sk (Figure 4d1), with peaks at the maximum dissipation. The negative Asy (Figure 4e1) with a large absolute value occurred at W2 and W1, representing intensive breaking. Correspondingly, the absolute value of undertow u ( Figure 4f1) increased with the incident waves. The results of B-N were set as the basis to analyze the relative changes of results in the other profile types as follows.
In tests of B-AR, although the wave energy at W13 was slightly enhanced because of reflection by AR, H s ( Figure 4a2) and H s-long (Figure 4b2) on its leeside were decreased, with maximum reductions of 0.05 m (at W5-W11) and 0.02 m (at W1 and W2), respectively. The reduction in wave energy was heightened with the incident wave energy. The shoaling effect was also reduced compared with B-N, and the peak of η at W11 was preserved only in B-AR-J4 (Figure 4c2). The Sk values at W1 and W2 were reduced under all wave climates (Figure 4d2). In the shoreface part, the Sk value increased at W4-W9 and decreased at W10-W13 in B-AR-J1. Then, as the incident wave grew, the wave gauge points turned to have decreased values relative to the corresponding B-N tests. For the cross-shore pattern of Asy (Figure 4e2), reduction in the absolute value at W1 meant weaker wave breaking relative to B-N. In the leeside of AR, u was decreased greatly in the range of 0.003-0.025 m/s, and u at W13 was increased by 0-0.009 m/s. The influence of the ASB was more complex. In B-ASB-J1 and B-ASB-J2, H s (Figure 4a3) and H s-long (Figure 4b3) changed slightly, with the changes less than 0.003 m and 0.006 m, respectively, which showed a small growth in the crest and the fore-slope of the ASB. For the waves of J3 and J4, H s barely changed at W11 and W13 but increased at W9-W10 (by less than 0.006 m and 0.004 m for J3 and J4, respectively), covering the fore-slope and the fore-part of the crest, and then decreased significantly at W2-W8 (by 13% and 23% for J3 and J4). Meanwhile, the cross-shore pattern of H s-long was similar to that in the B-N test but was different on the beach. In particular, H s and H s-long at W1 were increased for J2 and J3, with the maximum growth of 0.011 m in H s-long in B-ASB-J3. On the basis of η, Sk and Asy (Figure 4c3,d3,e3), wave shoaling, breaking, and bottom friction due to the ASB predominated in the wave transformation of J2-J4. Hence, there were local peaks of Sk and Asy between W6 and W10 covering the ASB, which indicated a stronger energy dissipation area rather than a beach slope. Moreover, u was reduced at W11-W13 and increased around the ASB (Figure 4f3). Beyond that, the ASB decreased u in the lee-side in a wider range for larger incident waves.
In comparison with the results in B-AR and B-ASB, the cross-shore hydrodynamics in B-ASB-AR showed a mixture of both features. The wave energy attenuation in B-ASB-AR was similar to that in B-AR, and the maximum reduction of H s (Figure 4a4) and H s-long (Figure 4b4) ranged between 39% and 50%. The hydrodynamics at W11-W13 depended on the AR rather than the ASB. Although η (Figure 4c4) resembled that in B-AR, Sk, Asy and u (Figure 4d4,e4,f4) changed in the same trend as those in B-ASB but with smaller absolute values. This was because the incident waves were attenuated by the AR before passing to the ASB, and bottom friction and wave breaking on the ASB further dissipated the wave energy before attack on the beach slope.
Generally, the cross-shore sediment transport coupled with wave breaking varied for different profile types. In B-N and B-AR, sediment transport was confined on the beach, where the upper beach was eroded, and sediment transported offshore until settlement at the lower area. Comparing Figure 4g1,g2, the maximum q s was diminished by the AR (by 44-64%), with the location moving shoreward under larger wave climates. In B-ASB and B-ASB-AR, the profile change could cover nearly the whole horizontal length, because the ASB itself covered about one-third of the whole profile and acted as a sediment source. The offshore sediment transport was generated on the beach and the fore-slope of the ASB, and onshore sediment barely existed in B-N and B-AR tests. Special cases were B-ASB-J4 and B-ASB-AR-J3, where onshore sediment transport dominated the area between the beach and the ASB.

Evolution process
Beach nourishment often causes disturbances, e.g., beach scarps and dunes, that eventually equilibrate via sediment transport processes [85]. The beach face evolution (shown in Figure 5) in all the tests started with a typical post-nourishment scarp erosion process summarized by Erikson et al. [86] as follows: continued wave attack on the scarp (in the present study) or dune caused undercutting or notching (removal of scarp or dune foot), which led to the visible tension cracks in the upper part; the resulting overhang eventually turned into sand sliding downward or toppling over, leaving a vertical face. Then, the fallen sand was carried offshore by the undertow accompanied by further beach undercutting in the next cycle. As a result, the eroded sand settled and quickly formed a breaker bar, where wave breaking happened continuously. The breaker bar gradually increased in volume and migrated offshore and thus weakened the onshore wave attack, which impeded the recession of the scarp location. Overall, the AR, the ASB, and their combination significantly reduced the beach change (upper erosion and lower deposition) and also caused different morphodynamic features.  This study focused on the different influences induced by the AR and the ASB, individually and in combination, on the evolution process. Relations were found between the dimension indicators and the duration time. The fitting parameters were assumed to reflect the comprehensive impact of the incident waves and the nourishment method, which are explored in the Discussion section.

Scarp
The scarp erosion is related to wave run-up, which reflects the different influences due to the incident waves and the beach profiles. Based on the beach evolution process (Figure 5), the change of scarp location with time was derived with a time interval of 10 min, as in Figure 6, which shows that the scarp location gradually retreated with time at a decreasing rate. The retreat rate Vsx in each time interval was estimated by calculating the first derivative of the scarp location data sx over time t (Vsx This study focused on the different influences induced by the AR and the ASB, individually and in combination, on the evolution process. Relations were found between the dimension indicators and the duration time. The fitting parameters were assumed to reflect the comprehensive impact of the incident waves and the nourishment method, which are explored in the Discussion section.

Scarp
The scarp erosion is related to wave run-up, which reflects the different influences due to the incident waves and the beach profiles. Based on the beach evolution process (Figure 5), the change of scarp location with time was derived with a time interval of 10 min, as in Figure 6, which shows that the scarp location gradually retreated with time at a decreasing rate. The retreat rate V sx in each time interval was estimated by calculating the first derivative of the scarp location data sx over time t (V sx = d(sx)/dt). It can be inferred that sx increased with t until V sx monotonically reduced to zero, corresponding to the equilibrium state with the maximum sx. As a result, a natural exponential function (exp) was applied to capture this process feature, leading to a scarp retreat model consisting of Equations (11) and (12): where the unit of t is minute. Equation (11) was firstly fitted to acquire parameter T e , which was then used to fit the sx data on Equation (12) to obtain parameters P1 and P2. under present hydrodynamics, and V0 represents the initial retreat rate. Thus, P2 is the estimation from experimental data on the maximum sx in the equilibrium state, and P1 represents the distance before the scarp location reaches P2, i.e., the difference between the minimum and the maximum of sx. The scarp model lays the foundation for quantifying the morphodynamic influence on the beach erosion process of the offshore artificial reef and submerged sand bar.   The fitting results are illustrated in Table 3, and the evaluation shows relatively good agreement. It is noted that P1 should equal to V 0 × T e . Therefore, V 0 is revised by P1/T e in this exponential scarp retreat model. The exponential part exp(−t/T e ) is practically regarded as a non-dimensional indicator for the time-varying feature, and the dimensions of T e and V 0 are time and speed, respectively. When t exceeds T e , exp(−t/T e ) is smaller than 0.36. Parameter T e indicates the time to reach equilibrium under present hydrodynamics, and V 0 represents the initial retreat rate. Thus, P2 is the estimation from experimental data on the maximum sx in the equilibrium state, and P1 represents the distance before the scarp location reaches P2, i.e., the difference between the minimum and the maximum of sx. The scarp model lays the foundation for quantifying the morphodynamic influence on the beach erosion process of the offshore artificial reef and submerged sand bar. The smaller P2 with smaller P1 means the maximum scarp location moving seaward with a reduced distance, and smaller T e can be interpreted as a quicker way to equilibrium. Comparing the parameter values in B-AR with those in B-N under the same wave, it was found that the AR decreased P2 (by 0.23-0.82 m), P1 (by 0.07-0.56 m), and T e (by 3-32 min), except for the increased T e under J4. The same comparison was made between B-ASB and B-N, and P1 and P2 were generally reduced by the ASB but less than the AR, except for the increase in B-ASB-J2. In addition, T e could be increased and decreased by the ASB. The results in B-ASB-AR were similar to those in B-AR, while T e was smaller than that in B-N. Overall, the AR showed a more significant effect in preventing scarp recession than the ASB, leading to smaller values of P2 and P1. However, there was no definite trend in the variation of T e , which is discussed with an explanation referring to the beach state.

Breaker Bar
As Figure 5 shows, the breaker bar migrates offshore with the plunging breaker during beach evolution. For the location of the breaker bar (x bar ), the bar first grows in height (z bar ) until the local elevation is high enough for the occurrence of wave breaking, which eventually turns into a new trough or terrace with the bar location further offshore (decrease in x bar ). In other words, there is a growth limit of bar height z bar for a specific location x bar on the beach, which indicates that the breaker bar moves along a certain space trajectory. Hence, this trajectory can be used to estimate the bar height based on its location.
Eichentopf et al. [87] proposed a linear function fitting, based on data from a large-scale flume experiment, represented as expressed in Equation (13) (R-square = 0.85). Likewise, the linear fitting by Equation (14) was conducted for our experiment, as illustrated in Table 4 and Figure 7.   Figure 7 (a-d). Therefore, the wave climate was not a restriction on this relation. Note that the target location (maximum elevation relative to the initial slope) should always be located on the beach (20.5 m < x < 27.3 m), i.e., the constant background slope, which is assumed to be the restriction of this linear relation.
In the experiment by Eichentopf et al. [87], the origin was set at the intersection of the still water level and the beach slope with the same axis direction that we set ( Figure 2). The beach profile, with a slope of 1:15, consisted of well-sorted sand with a narrow grain size distribution (d50 = 0.25 mm, d10 = 0.154 mm, d90 = 0.372 mm), which covered our target grain size. Hence, we converted the coordinate system as in Eichentopf et al. [87] and generated Equation (15) based on zbar = -0.05xbar + 1.17 to make the comparison.
Referring to Equation (13) (15)) with a depth of 0.064 m (calculated from a slope of 1/10). The ratio of the two water depths basically meets a geometric scale factor of 10. This result also supports the validity of using this light-weight sediment experiment to simulate fine sediment transport.

Artificial Submerged Sand Bar (ASB)
Concerning the ASB area where the water depth started to decrease from the fore-slope to the minimum over the crest, significant variation in hydrodynamic parameters occurred due to shoaling and breaking. The beach protection and nourishment were achieved through sacrificing the ASB, which lost its sand.
In the final profile shown in Figure 4 (h3) and (h4), the shape of the ASB evolved from the symmetrical shape to an asymmetrical shape with a steep fore-slope and mild rear-slope and even a bar-trough over the crest. The breaking zone on the ASB was basically located on the fore-slope and The evaluation parameters show that this linear relation was verified by this experiment. Considering the measurement precision and the horizontal range, the fitting parameters (q1 and q2) obtained in different profile types were regarded as equal. Therefore, the AR, the ASB, and their combination in this study affected the breaker bar location by decreasing the incident waves, thus resulting in a shoreward breaking location, not by changing the space trajectory. Furthermore, the fitting based on all the dimension data was conducted to attain the average results (q1 = −0.049 and q2 = 1.17) in this experiment. Eichentopf et al. [87] also pointed out that the coefficients q1 and q2 may depend on beach slope and sediment properties. Different incident waves resulted in different bar locations, which all accorded with the solid line shown in Figure 7a-d. Therefore, the wave climate was not a restriction on this relation. Note that the target location (maximum elevation relative to the initial slope) should always be located on the beach (20.5 m < x < 27.3 m), i.e., the constant background slope, which is assumed to be the restriction of this linear relation.
In the experiment by Eichentopf et al. [87], the origin was set at the intersection of the still water level and the beach slope with the same axis direction that we set ( Figure 2). The beach profile, with a slope of 1:15, consisted of well-sorted sand with a narrow grain size distribution (d 50 = 0.25 mm, d 10 = 0.154 mm, d 90 = 0.372 mm), which covered our target grain size. Hence, we converted the coordinate system as in Eichentopf et al. [87] and generated Equation (15) based on z bar = −0.05x bar + 1.17 to make the comparison.
Referring to Equation (13) (15)) with a depth of 0.064 m (calculated from a slope of 1/10). The ratio of the two water depths basically meets a geometric scale factor of 10. This result also supports the validity of using this light-weight sediment experiment to simulate fine sediment transport.

Artificial Submerged Sand Bar (ASB)
Concerning the ASB area where the water depth started to decrease from the fore-slope to the minimum over the crest, significant variation in hydrodynamic parameters occurred due to shoaling and breaking. The beach protection and nourishment were achieved through sacrificing the ASB, which lost its sand.
In the final profile shown in Figure 4h3,h4, the shape of the ASB evolved from the symmetrical shape to an asymmetrical shape with a steep fore-slope and mild rear-slope and even a bar-trough over the crest. The breaking zone on the ASB was basically located on the fore-slope and the seaward part of the crest, where the evolution was captured by video (as shown in Figure 8). The local net sediment transport rate at every 10 min q s was estimated to investigate the effect of the changed profile on the sediment behavior. eroded sediment was deposited near x = 12.0 m and formed a one-peak shape on the crest, which enhanced the wave breaking and shoaling effect. Then, the peak was steadily eroded by continued wave breaking and meanwhile also received the sediment from the middle crest, which led to the aforementioned dynamic feature from x = 10.5 m to the crest with a large magnitude ̅ (Figure 8 (g)). However, in B-ASB-J4 (Figure 8 (d)), the wave breaking area covered most of the fore-slope and the crest, where a large amount of sediment was transported offshore during the initial 30 min ( Figure  8 (h)), and several small bar-trough shapes developed. A two-peak shape (a large one accompanied with a small one) finally formed after 50 min (Figure 8 (d)), when a repetitive dynamic state developed after a small duration of overall onshore sediment transport (Figure 8 (h)).
In contrast with B-ASB, the morphological changes were significantly reduced in B-ASB-AR (Figure 8 (i)-(l)) under the same wave climate. The original shape of the ASB was basically retained under J1-J4, and local sediment transport was directed offshore during the initial responses ( Figure  8 (n)-(p)), except for B-ASB-AR-J1 with a little sediment transport around the ASB (Figure 8 (m)). The sediment transport entered a stably changing state earlier than in the corresponding test of B-ASB, with onshore-directed sediment transport of a small magnitude in the final stage. In general, the AR greatly limited the wave attenuation by the ASB and prohibited the onshore sediment from the ASB.

Discussion
The experimental results reveal the mechanisms that influence the morphological evolution by the offshore AR and ASB. In this section, the similarities and the differences among the four profile types are discussed with some relationships and comparison on the basis of previous works.

Beach Erosion
The beach profile behavior was characterized by a scarp and a breaker bar indicating erosion and deposition, respectively. The scarp retreat was caused by wave run-up, whose process was described by an exponential model with parameters depending on the incident wave energy because For the B-ASB profile type, the ASB changed very slightly under J1 with a very small q s (Figure 8a,e), which was directed offshore during the first 30 min and then stayed stable at nearly zero. The sediment response to J2 (Figure 8f) started roughly from x = 11.0 m, showed a repetitive dynamic feature that the direction of q s alternated between offshore and onshore with the same magnitude until reaching zero from 40 min to 50 min, and then went in another alternation. During the initial 50 min in B-ASB-J3 (Figure 8c), the incident waves broke at the middle part of the crest, leading to erosion near x = 13.0 m with offshore sediment transport (Figure 8g). As a result, the eroded sediment was deposited near x = 12.0 m and formed a one-peak shape on the crest, which enhanced the wave breaking and shoaling effect. Then, the peak was steadily eroded by continued wave breaking and meanwhile also received the sediment from the middle crest, which led to the aforementioned dynamic feature from x = 10.5 m to the crest with a large magnitude q s (Figure 8g). However, in B-ASB-J4 (Figure 8d), the wave breaking area covered most of the fore-slope and the crest, where a large amount of sediment was transported offshore during the initial 30 min (Figure 8h), and several small bar-trough shapes developed. A two-peak shape (a large one accompanied with a small one) finally formed after 50 min (Figure 8d), when a repetitive dynamic state developed after a small duration of overall onshore sediment transport (Figure 8h).
In contrast with B-ASB, the morphological changes were significantly reduced in B-ASB-AR (Figure 8i-l) under the same wave climate. The original shape of the ASB was basically retained under J1-J4, and local sediment transport was directed offshore during the initial responses (Figure 8n-p), except for B-ASB-AR-J1 with a little sediment transport around the ASB (Figure 8m). The sediment transport entered a stably changing state earlier than in the corresponding test of B-ASB, with onshore-directed sediment transport of a small magnitude in the final stage. In general, the AR greatly limited the wave attenuation by the ASB and prohibited the onshore sediment from the ASB.

Discussion
The experimental results reveal the mechanisms that influence the morphological evolution by the offshore AR and ASB. In this section, the similarities and the differences among the four profile types are discussed with some relationships and comparison on the basis of previous works.

Beach Erosion
The beach profile behavior was characterized by a scarp and a breaker bar indicating erosion and deposition, respectively. The scarp retreat was caused by wave run-up, whose process was described by an exponential model with parameters depending on the incident wave energy because of the same sediment and initial beach slope. Wave breaking occurred on the breaker bar after its formation, and it then sheltered the onshore area and affected the wave run-up. The key to figuring out the indistinct morphodynamic role of offshore interventions was to link the beach behavior and their impact on the hydrodynamics in the surf and the swash zones.
The wave gauge at W5 measured the incident waves for the beach area under offshore influence in this experiment. Therefore, correlation analysis was conducted between the scarp model parameters and the hydrodynamic factors calculated using measurement at W5, which may shed light on the role of the AR and ASB in the beach erosion regime. Apart from the significant wave height H s , long wave height H s-long, and mean water level η, other factors used to estimate the wave run-up included deep water wave length L s , wave steepness H s /L s , the products of wave height and length H s × L s and the Irribarren number ξ (with values between 0.64 and 0.92 calculated by Equation (16)). For each profile type, the Pearson's linear correlation coefficients were calculated separately for P1, P2, and T e in Table 5.
Parameters P1 and P2 of B-N were highly correlated with the hydrodynamic factors of H s , H s-long , η, L s , and the Irribarren number ξ, with coefficients over 0.90 and significance levels above 0.10, which is in line with the widely-used runup model of Stockdon et al. [88]. The coefficients in B-AR show that not only were these correlations enhanced, but also correlations were induced with the wave steepness and the product of the wave height and length. Nevertheless, the ASB basically maintained the correlations in B-N yet removed the influence of η. The coefficients in B-ASB-AR were similar to those in B-ASB, which revealed that the final influence of the combination may be determined by the ASB. In addition, there was no significant correlation confirmed between T e and any of the hydrodynamic factors. Therefore, it is speculated that the hydrodynamics alone cannot lead to a definite T e .  Because of the definite linear space trajectory of the breaker bar in this experiment, as proposed in Section 3.2.3, the effects of the AR and the ASB on the breaker bar were investigated by comparing the dimensions during the migration progress for different profiles under the same wave climate. The breaker bar formed where the wave broke on the beach and moved with the break point. It can be concluded that the wave energy was decreased by offshore interventions, leading to onshore movement of the break point and the breaker bar location with a smaller bar height along the trajectory. Figure 9 shows that the AR and the ASB resulted in a smaller breaker bar, while the effect of the ASB was less than that of the AR. The dimensions of the breaker bar in B-ASB-AR were basically at the same level as in B-AR except for the smaller height and the shoreward location under J1.

Beach State
The AR and the ASB changed the cross-shore hydrodynamics and the sediment transport, which reduced the rate and the magnitude of morphodynamic evolution on the beach. Wright and Short [89] proposed that the hydrodynamic processes and the relative contributions of different mechanisms to sediment transport and morphologic change differ dramatically as functions of the beach state (complete assemblages of depositional forms and coupled hydrodynamic process signatures in the surf zone and beaches). From this perspective, the offshore interventions may cause a transformation in the beach state, at least for a period of time, which was detailed in the section based on wave breaking and morphological features in the experiment. Since the flume experiments were under the normal model waves, analysis focusing on the influence of offshore interventions was limited to cross-shore processes regardless of the long-shore effect.
in Section 3.2.3, the effects of the AR and the ASB on the breaker bar were investigated by comparing the dimensions during the migration progress for different profiles under the same wave climate. The breaker bar formed where the wave broke on the beach and moved with the break point. It can be concluded that the wave energy was decreased by offshore interventions, leading to onshore movement of the break point and the breaker bar location with a smaller bar height along the trajectory. Figure 9 shows that the AR and the ASB resulted in a smaller breaker bar, while the effect of the ASB was less than that of the AR. The dimensions of the breaker bar in B-ASB-AR were basically at the same level as in B-AR except for the smaller height and the shoreward location under J1.

Beach State
The AR and the ASB changed the cross-shore hydrodynamics and the sediment transport, which reduced the rate and the magnitude of morphodynamic evolution on the beach. Wright and Short [89] proposed that the hydrodynamic processes and the relative contributions of different mechanisms to sediment transport and morphologic change differ dramatically as functions of the beach state (complete assemblages of depositional forms and coupled hydrodynamic process signatures in the surf zone and beaches). From this perspective, the offshore interventions may cause a transformation in the beach state, at least for a period of time, which was detailed in the section based on wave breaking and morphological features in the experiment. Since the flume experiments were under the normal model waves, analysis focusing on the influence of offshore interventions was limited to cross-shore processes regardless of the long-shore effect.
The beach states are classified as two extremes (i.e., fully dissipative and highly reflective states) and four intermediate states, i.e., longshore bar-trough (BT), rhythmic bar and beach (RB), transverse bar and rip (BR), and ridge-runnel or low tide terrace (LT). The surf-scaling parameter , Equation (17) [90], was originally used to distinguish the morphodynamic states.
where is the breaker amplitude estimated by half the breaker wave height Hb, is the incident wave radian frequency, g is the acceleration of gravity, and is the beach or surf zone slope gradient. Hb is calculated using Equation (18) [91] with the initial beach slope taken into account. The beach states are classified as two extremes (i.e., fully dissipative and highly reflective states) and four intermediate states, i.e., longshore bar-trough (BT), rhythmic bar and beach (RB), transverse bar and rip (BR), and ridge-runnel or low tide terrace (LT). The surf-scaling parameter ε, Equation (17) [90], was originally used to distinguish the morphodynamic states.
where α b is the breaker amplitude estimated by half the breaker wave height H b , ω is the incident wave radian frequency, g is the acceleration of gravity, and β is the beach or surf zone slope gradient. H b is calculated using Equation (18) [91] with the initial beach slope taken into account.
H 0 is the equivalent deep-water wave height considering refraction and diffraction. It was concluded that [89]: complete reflection can be expected when ε < 1; as long as ε < 2.0 − 2.5, strong reflection will continue to permit strong standing wave motion, surging breakers and resonance; when ε > 2.5, waves begin to plunge, dissipating energy; and spilling breakers occur when ε > 20. Moreover, the incident wave first broke at the breaker bar, and a second breaking would happen when the onshore beach area was eroded with enough depth under continuous wave attack. Hence, measurements at W5 and W1 were adopted to estimate the initial wave breaking for the outer and the inner surf zones, with the results shown in Table 6. It was found that plunging breakers dominated with ε between 2.5 and 20 except for the inner surf zone in B-N-J4, B-ASB-J3 and -J4, and B-ASB-AR-J3, where surging breakers arose.
The final experimental beach profiles in different tests are shown in Figure 10, where the slope gradients of the shoreward and the seaward breaker bars were marked as s1 and s2. Comparing the morphological features in different profile types, a bar-trough shape formed in B-N-J1, while small terraces without obvious breaking troughs were captured in the other three types of profiles. On the basis of the plunging breaker identified in both the outer and the inner surf zones (Table 6), the beach state in B-N-J1 was identified as BT, and LT was presented in B-AR, B-ASB, and B-ASB-AR under J1. Similarly, the beach states were confirmed as LT by the same criteria in B-ASB-J2, B-ASB-AR-J2, B-AR-J4, and B-ASB-AR-J4. Under the waves of J2 and J3, the beach profiles in B-N and B-AR both developed into a bar-trough shape, which slowly migrated offshore during the process ( Figure 5), with s1 being 0.05, which was linked with the beach state of the RB. Meanwhile, a long milder slope formed in B-ASB-J3 and B-ASB-AR-J3 with slope s1 under 0.05, and it was thus considered as a BR combined with wave breakers changing from plunging to surging ones, which was also the case in B-N-J4 and B-ASB-J4.  Comparing the identification in the B-N and the profile types with offshore interventions under the same wave climate, the AR and the ASB could induce a transformation in beach states towards a reflective state for some time. The ASB individually changed the beach state under J1-J3, while the AR individually changed the beach state under J1 and J4. The beach state in B-ASB-AR shows the result was closer to reflective state between the B-AR and the B-ASB.
Furthermore, the parameters of the scarp model in different tests were compared under the frame of the beach state, as in Table 7. It is noticeable that Te changed in the same trend as P1 and P2 in the same beach state, i.e., smaller Te corresponded to smaller P1 and P2, and larger Te came with larger P1 and P2 under the same type of beach state. Parameter Te can be seen as a representative of the hydrodynamic and the morphodynamic patterns, which had a real impact on the scarp retreat process by altering the time duration to the equilibrium state. It also indicated that comparison of the Te values should be confined within the same beach state, which accounted for the decreased P1 and P2 not necessarily corresponding to reduced Te.
Based on this assumption, the offshore interventions caused the beach state under J1 to transfer from BT to LT with the better protection effect of the AR and the combination compared to the ASB. For tests under J2 and J3, the beach state in B-AR preserved that in B-N, while the ASB and the combination led to more reflective beach states despite the possible increase in P1 and P2 (B-ASB-J2). Comparing the identification in the B-N and the profile types with offshore interventions under the same wave climate, the AR and the ASB could induce a transformation in beach states towards a reflective state for some time. The ASB individually changed the beach state under J1-J3, while the AR individually changed the beach state under J1 and J4. The beach state in B-ASB-AR shows the result was closer to reflective state between the B-AR and the B-ASB.
Furthermore, the parameters of the scarp model in different tests were compared under the frame of the beach state, as in Table 7. It is noticeable that T e changed in the same trend as P1 and P2 in the same beach state, i.e., smaller T e corresponded to smaller P1 and P2, and larger T e came with larger P1 and P2 under the same type of beach state. Parameter T e can be seen as a representative of the hydrodynamic and the morphodynamic patterns, which had a real impact on the scarp retreat process by altering the time duration to the equilibrium state. It also indicated that comparison of the T e values should be confined within the same beach state, which accounted for the decreased P1 and P2 not necessarily corresponding to reduced T e .  Based on this assumption, the offshore interventions caused the beach state under J1 to transfer from BT to LT with the better protection effect of the AR and the combination compared to the ASB. For tests under J2 and J3, the beach state in B-AR preserved that in B-N, while the ASB and the combination led to more reflective beach states despite the possible increase in P1 and P2 (B-ASB-J2). Conditions that decreased P1 and P2 with increased T e were found in cases where the beach state changed, such as B-ASB-J3 and B-AR-J2. All three parameters were decreased by the offshore interventions in the other cases. Overall, the offshore interventions prevented scarp retreat by moving the final position seaward and reducing the time duration or, if not, then by turning into milder beach states toward the reflective state.

Conclusions
This study experimentally investigated the cross-shore morphodynamics of a beach integrated with offshore interventions of an AR, an ASB, and the combination of both. For the beach without offshore interventions, the strong shoaling effect appeared in the offshore area for large incident waves. Wave attacks on the beach caused a scarp, and the fallen sand was carried offshore by undertow, then settled and quickly formed the breaker bar, where the most intensive breaking happened with the largest wave asymmetry and undertow and offshore sediment transport. The incident wave bores, after breaking continuously, attacked the scarp, leading to further beach undercutting, and meanwhile the breaker bar gradually increased in volume and migrated offshore. Therefore, the beach behavior was characterized by the scarp for erosion and by the breaker bar for deposition, respectively, while the scarp retreat process could be formulated as linearly related to the natural exponential of the duration time, and the breaker bar migrated with the bar height linearly depending on its location.
When the offshore interventions were included, the wave energy was dissipated in advance before reaching the beach. For the area in the leeside of the AR or the ASB, wave energy, shoaling, and breaking were reduced, which was more significant for larger waves. The AR caused the sediment transport in offshore-directed and was confined on the beach area with the maximum diminished by 44%-64%. The ASB not only attenuated the incident wave but also acted as a sediment source. The ASB crest was the main breaking area, leading to local erosion and strong undertow, which carried the sand offshore, depositing it on the fore-slope. As a result, the ASB evolved from a symmetrical shape to an asymmetrical shape with a steeper fore-slope, milder rear-slope, and a bar-trough over the crest under large waves. The morphodynamic responses of the combination of the AR and the ASB showed a mixture of the features of the AR and the ASB. The cross-shore distributions of the wave energy and the wave set-up or set-down were mainly controlled by the AR, while wave skewness, asymmetry, and undertow embodied the effect of the ASB.
The AR and the ASB decreased the wave attack on the beach, which resulted in smaller beach profile changes and the break point moving onshore. As a result, the scarp location moved seaward, and the breaker bar appeared at the shoreward location with a higher bar height. The ASB was less efficient than the AR during the process, where the beach state transformation might have been promoted further toward a reflective state under the same wave climate. Their combined effect was basically determined by the AR and the beach state towards a reflective element between the individual AR and ASB.