3.1. Long-Term Shoreline Response
The daily variation of beach width
for all beachlines are contoured over four years from June 2015 to July 2019 in
Figure 6a, where the numbers along the
y-axis represent the beachlines as marked in
Figure 3.
Figure 6a shows that the beach width actively varied not only in time but also in space along the coastline, although the variation pattern is not so regular. From January 2018,
in the G1 and G2 areas rapidly increased owing to beach nourishment (23,500 m
3). Therefore, the patterns in
after nourishment should be distinguished from those before it.
For the southern part of the beach (beachlines 1–5), rapidly decreased (i.e., the beach was eroded) from November 2015 and this continued until February 2017. Erosion also occurred along the mid-part of the beach (beachlines 10–15) starting from January 2016. In the northern part of the beach (beachlines higher than 15), increased from March 2016, indicating that the sediments eroded from the middle region moved to the north during this period. These changes in beach width resulted from seasonal fluctuation of total sediment transport due to effective wave action.
To separate the influence of longshore transportation from overall change implying seasonal recurrence, we differentiated
with respect to
and
.
Figure 6b,c display the distributions of
and
, respectively. To remove noise caused by daily variation and obtain longer-period changes,
for each beachline was averaged (monthly) before differentiation.
The pattern of
in
Figure 6b confirms longshore sediment movement. For example, large values of
gradually moved from beachline 10 in September 2016 to beachline 15 in February 2017, which shows that sediment moved to the north from the middle of the beach during this period.
It should be noted that the pattern of
in
Figure 6c shows clearer seasonal repetition, specifically in winter, regardless of nourishment. For example,
increased in the southern part of the beach (beachlines 1–10) in January of each year from 2016 to 2019, while it decreased in the northern part (beachlines 11–20) during almost the same winter period each year. This indicates that sediments at this beach might be transported from north to south along the coast during winter.
Figure 6d shows integrated beach widths along the coast;
are plotted to show the trend of gross variation of beach width. Although the patterns in sectors G1 and G6 are slightly different, the total beach area gradually decreased until the end of 2017 but rapidly increased early in 2018 after beach nourishment. Compared with total beach width,
and
in
Figure 7b,c do not include the beach fill effect and thus show time variation patterns of shoreline changes more clearly.
The integral of the rate
along the coast,
, is also presented in
Figure 6e. This quantity is related to the time variation of LST,
, as described by one-line theory in Equation (2). It should be noted that there are times when
values changed significantly within a period of several days. These periods of rapid change in
are likely related to periods of storm wave attacks, especially S99, which lasted more than 24 h (hereafter S99_24hr+). In
Figure 6, periods under the effect of S99_24hr+ are indicated by vertical gray lines in all panels. In addition, each storm event for S99_24hr+ is marked in
Figure 6e by a circle with a diameter representing the storm intensity. The storm intensity was calculated as the product of (maximum significant wave height during the storm event)
2 and the storm duration. The color of the circle represents the most frequent wave direction during each storm event with respect to the color bar on the left side.
After every storm event when S99 lasted more than 24 h,
increased significantly. Thus, it can be seen that severe wave conditions above a certain threshold can cause changes to the shoreline by increasing
; this change can last for months. This result indicates that LST caused by S99_24hr+ events resulted in conspicuous shoreline changes that lasted longer than changes made by cross-shore sediment transport. In the absence of S99 lasting more than 24 h,
decreased before beach nourishment in early 2018. In
Figure 6f,
is estimated for the south (G1), middle (G4), and north (G6) parts of the beach. After every S99 event, a positive peak in
was observed in G1 first, with peaks in G4 and G6 occurring later in the sequence. This indicates that sediments were transported to the north along the beach after severe wave events regardless of nourishment. Especially in January 2017, when the storm power intensity had sharp peaks,
in the G1 sector subsequently increased. Afterward, gross
and sectoral
in sectors G4 and G6 also increased. This indicates that sediment transport actively occurred in sector G1 and that sediments moved first to G4 and then toward G6.
In addition, weekly and monthly cumulative wave powers
are displayed in
Figure 6g. The cumulative
gives supplemental information on the behavior of LST. Especially in 2018, after beach nourishment, LST significantly increased, although the wave conditions were not severe even in the winter season. In addition, LST had a sharp peak in January 2019 after a slight increase in cumulative
(
Figure 6e). It is likely that sediment transport was accelerated even by small wave actions when the beach was filled with excessive sand after nourishment.
3.2. Short-Term Change in the Winter Season
Figure 7,
Figure 8,
Figure 9 and
Figure 10 show the time variation of wave parameters
estimated by Eqation (3) and beach width during four winter periods between 2015 and 2019, respectively, with an indication of high-wave-event periods. The first panel shows the daily mean of significant wave height
along with the daily mean of wave power
. The second panel shows the daily mode of
along with the daily mode of wave direction (
). Negative or positive values in the second panel imply that the waves propagate from the left or right side of the shore, respectively. In the third panel, the daily integrated beach width along the coast
is plotted with a dashed line to show the trend of gross variation of the daily beach width along with
, as it is related to
(solid line). Since it is an integration of the beach width for all beachlines,
is equivalent to the ‘gross beach width’. Similarly, ‘sectoral beach width’ can be defined as the integration of beach width,
, for the beachlines in each sector as
. In the last panel,
(dashed line) and
(solid line) in sectors G1 and G6 are presented separately. Here, the zero line indicates the average over the observation period and thus the positive and negative values of
and
indicate higher or lower values over the average, respectively. In every panel, columns in gray, blue, and red represent storm periods of S90, S95, and S99, respectively. Since periods over
also qualify as periods over
and
, the columns of S99 include the columns of S95 and S90. Similarly, the columns of S95 include those of S90.
The change of beach width represented by
and
during high-wave events is more conspicuous than that during periods of lower waves. During the winter between 2015 and 2016, as shown in
Figure 7, a total of four S99 events were observed. The first two S99 events were characterized not only by the wave intensity but also by their duration, as they lasted longer than five days. The third and fourth S99 events lasted for two and three days, respectively. In all S99 events, the total beach width
decreased considerably (dashed line in
Figure 7c), which indicates that the beach was rapidly eroded owing to offshore movement of beach sediments by the storm waves, as has been commonly observed (e.g., [
30]). It should be noted that, following the decrease in
, the estimated LST,
, increased; then, the eroded beach width began to recover as the decrease in
, the increase of
, and the recovery of
occurred sequentially for S99 events. The recovery of beach width after strong LST was likely related to
(
Figure 7b) rather than
(
Figure 7a), confirming the role of wave direction. As shown in
Figure 7d, during the winter after the two S99 events,
increased while
decreased owing to the negative wave direction. This indicates that the sediments were eroded from sector G6 and transferred to sector G1 by strong LST. It also shows that the beach width gradually recovered after severe erosions while the sediments were generally moved alongshore, not offshore, by the obliquely incident waves. The beach width eroded by the first S99 event in November did not recover its pre-storm condition during the rest of the winter period. It is likely that the first attack of strong S99 that lasted more than five days resulted in a non-recovery condition. Moreover, the time gap between the first and second S99 events was relatively short (only 10 days), which prevented the beach from recovering. In addition, three short S95 storms developed between the first and second S99 events, which also reduced the beach’s resilience for recovery.
In the winter between 2016 and 2017 (
Figure 8), there were more frequent storm events but with lower intensity and shorter duration than during the previous winter season. Owing to frequent storm attacks, the total beach width gradually decreased as the magnitude of
became lower at the end of the winter compared with that at the beginning (dashed line in
Figure 8c). This indicates loss of sediments by dominant offshore sediment transport under harsh wave conditions over the winter period. In addition,
variation patterns at the two end sectors of G1 and G6 showed clear differences for this winter;
in sector G1 increased, whereas
in sector G6 decreased (dashed lines in
Figure 8d). This shows that sediments were moved from sector G6 to sector G1 by LST. The opposite trend by sector confirmed the role of LST caused by the action of obliquely incident waves in the winter season, as in the previous winter.
Figure 8c also shows that
increased following the increase in
during storm events, confirming that the beach width increased as LST increased, as observed in the previous year. Therefore, although the gross beach width generally decreased owing to stronger offshore sediment transport, there were times when the beach recovered according to the increase in LST.
This specific pattern was also observed in the winter between 2017 and 2018 (
Figure 9), when the beach width increased following the increase of LST during severe storms. However, in general, the beach width increased as
at the beginning of the season was lower than that at the end. This increment of beach width can also be seen clearly in
Figure 9d, where
gradually increased over the winter period for both G1 and G6, indicating that LST during storm events did not contribute to overall shoreline variation. This likely reflects a long low-wave energy period between December 2017 and February 2018, during which no S99 events were observed. Therefore, the shoreline erosion by storm waves at the beginning of the winter was recovered during the time of weaker wave conditions, when obliquely incident waves were dominant.
In the winter between 2018 and 2019 (
Figure 10), the wave conditions were milder compared with those of previous winters. There were no S99 storm events observed and even the number of S90 and S95 were also significantly smaller compared with previous winter seasons. In addition, the level of
was significantly higher compared with that in previous seasons. This high
was due to beach nourishment implemented in early 2018, which provided a significant amount of sediment. The magnitude of wave energy was lower, but the beach width was significantly higher than during other winter seasons (see different scale on the
y-axis of
Figure 10 compared with those of
Figure 7,
Figure 8 and
Figure 9). Despite having the mildest wave conditions, the variation pattern (where LST increments follow beach width reduction owing to CST during storm events) was repeated in this season. Following increments in LST, the eroded beach width gradually began to recover. Moreover, clearly different patterns between
in two sectors were also detected. The increased and decreased beach widths in sectors G1 and G6 resulted from obliquely incident waves during this winter season.
3.3. Event-Scale Analysis
We investigated the patterns of observed shoreline variation in terms of wave conditions during specific events including both severe and rapid variations in cases of extreme storm conditions, and the shoreline response under milder wave conditions.
When extreme storm conditions (S99) lasted more than five days, the beach was rapidly eroded by CST. For example,
Figure 11 shows that, when storm waves attacked the beach on 6 November 2015,
(dashed line) rapidly decreased, indicating erosion. However, the LST estimated by
(solid line) did not start to increase until 7 November 2015, although
peaked on 6 November 2015. This indicates that the beach was first eroded owing to CST but LST responded later to obliquely incident waves. A similar pattern occurred on 22 November 2015, when another S99 storm event started. While
decreased, corresponding to a storm event, LST only started to increase on 23 November 2015, one day after the storm beginning (22 November 2015). The combined effect of CST and LST during extreme storm periods resulted in a characteristic recovery pattern following the storm. In general, the shoreline eroded by CST was recovered in most sectors of the beach. However, in some sectors the shoreline recovery was slower; this occurred when sediment was lost by LST. For example,
in most sectors gradually increased owing to recovery; however, the recovery was not clearly observed in sector G1 as the sediments in this end sector were likely lost by LST.
Figure 12 (December 2016) shows a similar pattern;
decreased when an S99 event occurred and
increased following erosion. Specifically, the beach width in sector G1 sharply increased after 16 December 2016 following the end of an S99 event, while
at other sectors did not show clear changes. This indicates that sediments moved to sector G1 owing to LST caused by the storm events.
Under milder wave conditions, rapid and significant changes in beach width did not occur. Instead, LST caused by obliquely incident waves carried sediments along the beach toward the end sectors. In
Figure 13, neither
or
show significant changes from 15 November until 7 December 2017. However, the beach width in sector G1 gradually increased, while it decreased in sectors G5 and G6. In the corresponding period,
remained less than 2 m but the wave propagating direction,
, actively varied from 0 to +90°, indicating that waves were approaching from the left of the beach (
Figure 3). Therefore, sediments were consistently moved alongshore toward sector G1 by LST, regardless of low-wave energy corresponding to negative values of
.
3.4. Shoreline Response after Beach Nourishment
One of the most significant changes in beach width during the experimental period was caused by beach nourishment implemented in early 2018. As shown in
Figure 6a, the beach width sharply increased in sectors G1–G4 from January 2018, showing that nourished sediments were placed in sectors near the southern end of the beach. As a result, after the nourishment, the total beach width,
, increased from about 30 m
2 in early 2018 to about 150 m
2 in early 2019 (
Figure 6d). When a large amount of sediment was placed on the beach face by nourishment, the shoreline rapidly responded to reach equilibrium in both cross-shore and longshore directions. In
Figure 14, the wave and shoreline parameters are plotted for a period of about one month in June 2018, when the nourishment had been completed. Starting from 8 June 2018,
and the magnitude of
increased for about 10 days until 18 June 2018. However, the increment of wave power was not significant compared with those observed in storm periods such as S90–S99, and even the maximum wave height was less than 2 m. However, the total beach width,
, rapidly decreased from June 8 to 16 June 2018, and the LST estimated from
increased over the same period (
Figure 14c). The significant changes in the beach width and LST regardless of the comparably low wave energy in this short period likely occurred because the beach reached its equilibrium status after nourishment. For example, the beach width in sector G1 sharply decreased while
in sectors G2–G5 did not change significantly. Considering that the nourishing sediments were placed at the southern end of the beach, the rapid change in sector G1 can be explained by the equilibrium process.
In
Figure 15a, the gross beach width is compared between sectors for each winter from 2015. It is clearly seen that the beach width in sector G1 immediately increased immediately after nourishment in the winter between 2017 and 2018. After that, beach width increased with increasing sector number, as the increment of beach width in sector G6 was at a minimum in the winter between 2018 and 2019. This pattern is also observed in
Figure 15b where the beach width distribution in all sectors is contoured. The red dashed arrow in the figure shows that the sediments that were placed in the lower sectors gradually moved to higher sectors as soon as nourishment was implemented, which can be observed in the aerial photographs taken from Google Earth at four different times steps during the observation period (
Figure 15c).
The impact of the beach fill in which sediments were placed at lower sectors only in early 2018 was not immediately reflected in higher sectors such as G5 and G6. Therefore, beach widths in sectors G5 or G6 still displayed an erosional trend in the winter of 2017–2018, which continued from the previous winters. Before the beach fill of early 2018, all sectors showed the same decreasing trend until the winter of 2016–2017. Therefore, if there was no beach fill, the erosion trend would have continued in all beach sectors, although the rate of erosion might have diminished as waves became milder in subsequent winters.