4.1. Characteristics of Shoreline Change
The long-term variation of the shoreline location was examined.
Figure 5 shows the temporal variation of the alongshore-averaged shoreline location along a 7.2-km stretch of the marine drive region from 1986 to 2006. Over that period, the shoreline locations in March and September indicate a retreating trend of 1.1 m/year and 0.8 m/year, respectively. On average, the beach width has decreased approximately 20 m in the last two decades. The overall eroding trend is considered to be related to the construction of Kanazawa Port and the development of the Tedori River basin [
2,
3,
4]. It is quite important to examine the possible causes of accelerated erosion, but such an examination is beyond the scope of the present study. The shoreline retreat was most accelerated during the periods 1988 to 1991 and 2001 to 2004.
The features of temporal variation are generally similar in spring and autumn surveys. Quantitatively, the year-to-year variation is more significant in September, and the anomaly from the long-term trend reached 10 m in some years. The seasonal variations from September to March and from April to September are plotted in
Figure 6. Note that the shoreline retreat during autumn and winter (September to March) is slightly greater than the shoreline advance during spring and summer (April to September), and this difference has increased recently. This is related to the aforementioned difference in the rate of shoreline retreat in March and September. Moreover, the magnitude of shoreline change has been becoming larger recently.
Figure 7 demonstrates the temporal variation of mean significant wave heights observed at Kanazawa Port in contrasting seasons (autumn and winter: from the beginning of October to the end of March; spring and summer: from the beginning of April to the end of September). The wave height in the autumn and winter indicates a weakly increasing trend. This is considered to be related to the recent increase in shoreline variation during winter.
The variation of shoreline location was not uniform over the stretch and period; several kinds of variation with different spatial and temporal scales coexist.
Figure 8 represents the alongshore comparison of the rate of shoreline variation in the first (1986 to 1995) and second (1996 to 2005) halves of the study period. Note that a three-point moving-average was conducted in space to smooth out the local variation. The rate of variation was computed by linear regression. In general, a significant erosional trend is observed, except for the northern end of the study area where an accretionary trend is observed related to the influence of the construction of nearby Taki Port. In the first half of the study period, areas with shoreline advance and retreat appear alternatively. In the second half, the erosional trend became quite dominant. The locations where strong retreat of the shoreline was observed in the second half of the study period roughly correspond to the areas of shoreline advance in the first half. Shoreline retreat over the whole period was more significant in the northern part of the study area. The rate of shoreline retreat was not uniform during the study period. In some locations, it varied significantly between the first and second halves of the study period.
On the basis of the cross-shore survey record, the variation of shoreline location and sediment volume per unit width was computed for each line over the period 1998 to 2010. Although the shoreline analysis indicated a strong retreating trend in the range from 1.5 to 3.0 m/year, the results of volumetric change indicated quite different features; the change in sediment volume is generally small, and in fact, even a weakly increasing tendency was observed. This result implies that the trends in the variation of seabed-level change differ between nearshore and offshore regions. To clarify the difference related to the cross-shore locations, temporal variations of depth contour lines were examined for each line. As shown in
Figure 9, the cross-shore seabed slope increases for lines that are located more to the south, which is consistent with the previous field observations [
1,
6]. This feature is more pronounced in areas where the water depth exceeds 7 m. On lines H40 and H60, located in the central and north parts of the study area, the depth contour lines offshore are slightly advanced with time, whereas the contour lines near the shoreline indicate a retreating trend. As a result, the slope of the seabed is tending to decrease. Because the data coverage of cross-shore survey data used in this study (approximately 10 years) is not sufficiently long and the magnitude of change is not very different from the order of sounding accuracy, further efforts are needed to continuously monitor bed-level changes and to clarify the possible cause of these peculiar changes.
4.2. Characteristics of Systematic Sandbar Migration
Large-scale multiple sandbars have evolved in the study area. In
Figure 9, it is clearly seen that interannual movements of these bars are significant, and the migration of the bars has been repeated several times.
Figure 10 presents typical examples of observed interannual variations of cross-shore profiles on line H40 from November 2000 to November 2007. The bars migrate in a net offshore direction, and the transition of an individual bar consists of generation, seaward migration and decay. When the bar evolves most significantly, the bar height can reach 3 to 4 m. Similar types of temporal variation are observed on other lines. This kind of periodic movement of multiple sandbars is similar to the net offshore migration of bars reported at various locations around the world including the Netherlands [
13,
14,
17], New Zealand [
15,
17] and Japan [
16,
17,
18].
The relationship between the crest and trough elevations is plotted in
Figure 11. During the evolution of the bar (roughly corresponding to a crest in elevation of
Zc > −4 m), the water depth at the crest increases nearly proportionally to the depth at the trough. The slope of a linear regression line based on the least square method is in the range from 1.8 to 1.9. These values are slightly larger than those determined for the west coast of the United States (= 1.69) [
26] and for various Japanese coasts (= 1.47) [
5]. In the area where
Zc < −4 m, the trough elevation slightly increases with an increase of crest depth on the southern lines (H01 and H03), whereas on the northern lines (H40 and H60), a nearly constant value is taken around the closure depth. As a result, the bar height initially increases with depth at the crest (
i.e., seaward distance) as shown in
Figure 12. It reaches its maximal value around
Zc = −4.2 m. The maximum bar height observed is as large as 4.8 m. Further offshore, the bar height decays progressively. This is consistent with a previous report for various Japanese coasts [
5]. The temporal variation of the crest location and bar height is shown in
Figure 13 for each line. In this figure, the centers of the circles indicate the locations of bar crests, and the areas of the circles are proportional to the bar heights. These figures clearly indicate that the bars have repeated systematic interannual migrations on all of the lines. The periodic cycles consist of generation, seaward migration and offshore decay. Bars were first generated at a distance of 100 to 200 m from the shore and then migrated seaward. The bar height increased with cross-shore distance until it reached 400 to 700 m. When the offshore bar evolved most significantly, the height became as large as 3 to 4 m. Then, the bar migrated further seaward and diminished at an offshore distance of around 1000 m. After the outermost bar diminished, the next bar replaced the former one. The velocity of offshore migrations was approximately 100 m/year within 400 m of the coast. In the area lying approximately 400 to 800 m from the shore, it was accelerated to 150 to 200 m/year. Further offshore, the migration velocity was not more than 100 m/year. The life span of an individual bar from generation to offshore decay is approximately eight years. The return period, which is defined as the interval between each cycle, is estimated to be three to four years. During the study period, three or four cycles of periodic movement were recognized.
An alongshore comparison of bar migration was conducted next to examine the plan-view shape and migration characteristics of the bar system.
Figure 14 displays examples of the comparisons of cross-shore profiles on the survey lines observed during 2003. In the figure, triple- or four-barred profiles are observed on each transect. From the observation of aerial photographs, the individual bars are generally considered to be continuous in the alongshore direction. To facilitate the example of an alongshore continuous bar configuration, estimated crest lines are plotted by a broken line and open circles in
Figure 14. The seaward distance of the bar crests from the shoreline decreases as the location of the line moves southward (
i.e., the seabed slope becomes larger). As a result, the plan shape of an individual bar is slightly curved with respect to the shoreline. In
Figure 14, the alongshore variability can be clearly seen. For example, the outermost bar on line H01 is substantially decayed on line H03. Further north on lines H40 and H60, the bar is considered to have disappeared already. On the contrary, the third bar (from the offshore) on line H60 has not evolved yet on H40, H03 and H01. The interannual movements of the longshore bar, which corresponds to the crest line in
Figure 14, were tracked and compared on the four lines. The comparison shows that the seaward migration of the bar progressed from north (H60) to south (H01). In summary, the longshore bar is slightly closer to the shoreline as the alongshore location moves southward. The individual bar migrates net seaward, keeping the slightly curved plan shape with respect to the shoreline, and finally decays from the north to the south.
The cyclic movements of the crest lines are plotted in
Figure 15. The centers of the circles indicate the bar crest, and the areas of the circle are proportional to the bar height. To distinguish a series of periodic sequences of generation, migration and decay, crest markers that are considered to belong to the same generation of bar (
i.e., generated in nearly the same period) are indicated by the same color in the figure. For 1998, the bar indicated in yellow evolved as the outermost bar. Behind it, the next (purple) and the third (blue) bars also evolved. As time goes on, the outermost bar (yellow) migrates seaward and finally disappears from the north to the south. The cyclic movement indicates common features among the four lines. After one cycle, a configuration similar to that of 1998 reappeared in 2001. Another cycle was observed during 2001 and 2005, when the transition seems quite similar to the previous one. The bar behavior in the third cycle is different in some ways to the first and second cycles. Around 2007, the outermost bar (blue) disappeared in most of the area, and then, the next bar (red) started accelerating seaward migration. The migration is significant on the three lines H03, H40 and H60. In contrast, the old outermost bar (blue) remains on H01, and therefore the movement of the next generation of the bar (red) on H01 is much smaller. It is deduced that such an alongshore difference in bar migration resulted in the disconnection of the longshore bar between H01 and H03. From 2008 to 2009, the new bar (red) was considered to have realigned with the old bar (blue). It is assumed that the bars located more landward (light blue and red) have also realigned at the same time. The disconnection and realignment of a bar system has been reported at several NOM sites as ‘bar switching’ (e.g., [
15]).
4.3. Principal Modes of Sandbar Migration
Next, EOF analysis was applied to cross-shore profiles on the four lines. The first empirical eigenfunctions explain over 99% of the total variation. The second and third eigenfunctions explain 70% to 80% of the remaining depth variation. The variance explained by higher modes is much smaller. Therefore, the focus is mainly placed on the top three modes of variation. Furthermore, variation of cross-shore profiles appears most clearly in the EOF results up to X = 1300 m, so even though the cross-shore range of EOF analysis extends to X = 2400 m, plots are limited to the shorter distance.
The temporal and spatial eigenfunctions for the first mode of long-term variation are shown in
Figure 16. The spatial function of the first mode (
e1) expresses the mean profile of the seabed during the study period and is generally similar among the four lines. Detailed inspection indicates that, on line H01, the profile of
e1 has local maxima that correspond to the preferential locations of bar crests. The temporal coefficient of the first mode (
C1) is shown and compared to the cumulative change of sediment volume per unit width (Δ
V) with respect to the value at 1998. It is noted that the time coefficients are plotted after being multiplied by −1 (−
C1) in the figure for comparison. Note that on all four lines,
C1 is very closely related to the volumetric change induced by the imbalance of alongshore sediment transport or cross-shore exchange of sediment with the backshore and dunes or offshore regions.
The combination of the second and third modes corresponds to the cyclic behavior of the sandbar system.
Figure 17 describes the temporal and spatial eigenfunctions for the second and third modes. Both of the temporal functions,
C2 and
C3, show periodic variations. The phase difference between these two curves is approximately a quarter of the variation period. Similar features of phase shift are seen for the spatial eigenfunctions; the spatial functions
e2 and
e3 have features such that when one curve becomes a local maximum or minimum, the other crosses the horizontal axis. This implies that the relative difference in space between
e2 and
e3 is approximately a quarter wavelength. By synthesis of the characteristics of temporal and spatial functions, the combination of the second and third modes expresses the periodic seaward migration of bars.
Figure 18 shows the summation of these two modes. The reconstructed data clearly reproduce the cross-shore periodic movement of the outer bar with an approximately four-year cycle. Note also that each bar first appears near the shoreline and then migrates seaward. The size of the bar increases with offshore migration. The net offshore movement continues throughout the bar cycle at this site. After the crest location exceeds 700 to 800 m, the bar begins to degenerate and finally disappears. The EOF results indicate that the life span of individual bars is around eight years on these lines. Although these two modes of bar migration are accompanied by short-term changes in sediment volume among the sub-areas in the profile, it can be deduced that these two modes do not contribute to the long-term variation of sediment volume over a whole transect because there are no long-term trends in the variations of
C2 and
C3.
In the EOF analysis results, the magnitudes of the temporal functions of the second and third modes (C2 and C3) are generally larger than those for higher (fourth to sixth) modes (C4 to C6). Concerning the spatial functions, the magnitudes of the second to sixth modes (e2 to e6) are comparable in the bar migration zone where large-scale bed-level variation occurs. Because the variation in each mode is expressed by the product of the temporal and spatial functions, the second and third modes are dominant compared with higher modes in the bar migration zone. Near the shoreline, however, the magnitudes of e4 to e6 become larger than e2 and e3 in relatively many cases. Accordingly, the influences of the fourth to sixth modes are not negligible over the local area near the shoreline. Inspection of the spatial and temporal functions of higher (fourth to sixth) modes revealed that the long-term retreating trend of the shoreline is expressed by the fifth mode for lines H01, H40 and H60 and the sixth mode for line H03. It was also seen that the short-term (year-to-year) variation is substantial in the fourth to sixth modes. These higher modes should be taken into account when the long-term trend and year-to-year variation in the shoreline location need to be reproduced.
Finally, the relationship between bar migration and shoreline movement is examined. In general, the location of the shoreline is affected by various factors with different time scales. The annual survey record reflects the combination of the long-term (decadal) trend of erosion, with middle-term (several years) cyclic variations related to interannual bar migration and short-term variations (seasonal variation and short-term changes induced by storms). In the EOF analysis, these factors are decomposed into different eigenmodes. The middle-term shoreline movement related to NOM behavior appears in the second and third modes. Corresponding to the progress of NOM cycles, the cyclic variation of bed level is observed around
X = 0 in
Figure 18a to
Figure 18d.
Figure 19 compares shoreline location changes in the original survey with reconstructed EOF results for the first three (
Z1-3 =
C1e1 +
C2e2 +
C3e3) and five (
Z1-5 =
C1e1 +
… +
C5e5) modes for line H40. The reconstructed variation expressed by
Z1-5 adequately reproduced the long-term retreating trend of the shoreline. Among the overall variation,
Z1-3 extracts oscillatory variations with amplitudes of approximately 8 m without containing the long-term trend. Hereafter, this kind of cyclic variation accompanied by the NOM phenomena is examined to clarify the relation between bar configurations and shoreline movement.
Figure 20 and
Figure 21 demonstrate typical examples of cross-shore profiles corresponding to times when the reconstructed EOF results up to the third mode (
Z1-3) indicate the retreat (
Figure 20) and advance (
Figure 21) of the shoreline. The shoreline tends to retreat when the outermost bar is decaying and a new bar is being generated at a seaward distance of some 200 m (
i.e., when a triple bar configuration is newly established). This feature is common in lines H03, H40 and H60. In contrast, the shoreline experiences an advance when the outer bar is most evolved. More specifically, on lines H40 and H60, shoreline advances are observed when the old outermost bar remained, but with substantial decay, and the next bar evolved dominantly as the middle bar in a triple bar configuration. On lines H01 and H03, the shoreline experiences an advance when the new outermost bar has evolved significantly in the double bar configuration (the old outermost bar has almost completely diminished). On lines H40 and H60, the shoreline location is neutral when the old outermost bar has completely decayed, and the profile has a double bar configuration.
As a typical example, consider the variation on line H40. The spatial functions in
Figure 17c show that the second mode on line H40 corresponds to a triple bar configuration with a new bar generated at around
X = 200 m, whereas the third mode indicates the double bar configuration (the third bar has not evolved yet). As was mentioned previously, the temporal coefficients of these modes have a phase shift of approximately a quarter of a period. Accordingly, when the second mode is dominant (
i.e., a large
C2 magnitude), the influence of the third mode decreases
(i.e., corresponding to a small
C3 magnitude). Therefore, the bar system is close to a double bar configuration when the third mode is dominant
(e.g., 11/2003 in
Figure 10), while the transition to a well-developed triple bar system is observed when the second mode is dominant instead (e.g., 10/2004). The cross-shore distribution of
e2 indicated that the third bar evolves at around
X = 200 m when
C2 is maximum. At the same time, over the landward area corresponding to
X < 200 m, where
e2 is negative, the bed level decreased, and correspondingly, the shoreline location retreated
(e.g., 10/2004 in
Figure 10). Therefore, the generation of a new bar and resulting transition into a triple bar configuration in a bar system may induce a temporal retreat of the shoreline, whereas the evolution of the outer bar accompanies an advance in the shoreline
(e.g., at 10/2006 when
C2 had a negative peak). Note here that this kind of middle-term shoreline change over a time scale of years is cyclic and does not contribute to long-term changes in shoreline location within a decadal time scale.
These results imply that on sandy beaches similar to those on the Chirihama Coast where large-scale bar systems have evolved, the temporal variation of shoreline locations with a time scale of several years may appear with the systematic transition to a multiple bar system configuration. Accordingly, it is very difficult on such sandy beaches to estimate the variation of sediment volume over the transect from only the change in the cross-shore shoreline location with an assumption of parallel translation in the cross-shore profile. Care should be taken on this point from an engineering point of view.