In a slender water body with a large aspect ratio, the angle of wind waves relative to the direction normal to the shoreline may exceed 45°, resulting in the emergence of cuspate forelands and the segmentation of the water body. The BG model was used to predict the segmentation of a rectangular water body by wind waves when the probability of occurrence of the wind direction is given by a circular or elliptic distribution, and the segmentation of a rectangular water body into a circular or elliptic lake was predicted in each case. The segmentation of a shallow water body with a triangular or crescent shape was also predicted together with the prediction of lakeshore changes when a rocky or sandy island exists in a circular lake.
- closed water body
- wind waves
- lakeshore changes
- BG model
- cuspate foreland
In a shallow water body, beach changes may take place owing to wind waves. In a narrow water body with a large aspect ratio, the angle of wind waves relative to the direction normal to the shoreline may exceed 45°, and the shoreline may become unstable because the fetch distance in the direction of the principal axis of the water body is sufficiently large for waves with significant energy to be generated [1, 2]. Therefore, cuspate forelands that develop from both shores of a narrow water body connect with each other, resulting in the segmentation of the water body into smaller rounded lakes [3, 4]. For example, Figure 1 shows the segmentation of a lagoon facing the Chukchi Sea in Russia [3, 4]. In this example, five elliptic lakes can be observed as a result of segmentation, and their axes are parallel to each other. Figure 2 shows an enlarged satellite image of the rectangular area in Figure 1, and in this image, lake segmentation at a primitive stage can be seen with the alternate development of cuspate forelands. Regarding these phenomena, the division and reduction of the fetch distance owing to the formation of a large shoreline protrusion associated with shoreline instability under high-wave-angle conditions and the resulting change in the wave field are key factors. Ashton et al.  predicted that the forelands formed along the shoreline connect with each other, resulting in the segmentation of the water body into smaller rounded lakes. Uda et al.  predicted the three-dimensional (3-D) segmentation of a shallow rectangular water body using the BG model (a model for predicting 3-D beach changes based on Bagnold’s concept) . Uda et al. studied the emergence and mergence of small lakes and their segmentation using the same model , assuming that the wind blew from all directions between 0 and 360° with the same probability of occurrence and intensity, that is, a circular distribution of the probability. The segmentation into elliptic shapes, as shown in Figure 1, was not predicted in their study. It may be accomplished, assuming that the probability of occurrence of the wind direction is given by an elliptic distribution, similarly to the case of oriented lakes . In this study, the segmentation of a rectangular water body was predicted, given a circular or elliptic distribution of the probability of occurrence of the wind direction.
In the coastal area, the segmentation of a shallow water body with a triangular or crescent shape can also be observed. To study the mechanism of segmentation of such a water body, several examples of segmentation together with the development of sand spits along the lakeshore were examined in Lagoa de Mangueira in Brazil, Lake Saroma and Lake Kitaura in Japan. Then, the BG model was used to investigate the segmentation of a shallow water body with a triangular or crescent shape, and 3-D beach changes during the segmentation of a shallow water body into small lakes were predicted.
When wind waves are incident to the lakeshore in a closed water body with a rocky or sandy island, topographic changes may occur on the lee of the island because of the wave-sheltering effect. Since a rocky island is fixed at a location in a closed water body, the wave-sheltering effect of the island is constant with time, and the lakeshore converges to a certain stable form after the wave action for a sufficiently long time. When a sandy island is located in a closed water body, however, the island itself can deform owing to the action of wind waves, resulting in the successive change in wave field. Thus, more complicated lakeshore changes will occur. Here, Lake Balkhash located in Kazakhstan  was selected as an example, and the BG model was used for predicting lakeshore changes when a rocky or sandy island exists in a circular lake.
2. Examples of cuspate forelands in lake
2.1. Lakeshore in Lagoa de Mangueira
Figure 3 shows an example of segmentation and the development of sand spits along the lakeshore of Lagoa de Mangueira (location: 33°09′59”S, 52°49′32”W) in Brazil . The crescent lake is 100 km long and 10 km wide at the center of the lake. Many cuspate forelands have developed along the lakeshore and, in particular, the intervals of the cuspate forelands formed on the west shore become short near the south end of the lake. On the other hand, sand spits with similar shapes and cuspate forelands have developed along the east and west shores, respectively, in the north part of the lake. This is a typical example of segmentation and the development of sand spits in a crescent lake.
2.2. Lakeshore in Lake Saroma
Figure 4 shows a satellite image of an abandoned inlet located at the east end of Lake Saroma in Hokkaido, Japan (Location: 44°08′03”N, 143°59′04″E) . This water body has a triangular shape, and segmentation of a water body can be seen near the east end. When enlarging the rectangular area in the satellite image of Figures 4 and 5 is obtained. For the wind rose in Lake Saroma, the predominant wind direction is WNW, resulting in the eastward development of sand spits. On the south shore is Tofutsu fishing port, as shown in Figure 5. East of this fishing port, the width of the water body gradually decreases, and sand spits A-E develop together with pairs of sand spits A’ and E’. Of these sand spits, sand spits A and A’ are the largest and divide the water body into two. In the vicinity of sand spit A’, wind waves cannot be generated in the presence of the westerly wind, resulting in no development of sand spits. However, east of sand spit A’, wind waves can develop, and the size of the sand spits increases eastward in the order of sand spits B, C, and D.
2.3. Lakeshore in Lake Kitaura
Lake Kitaura located in Ibaraki Prefecture is a shallow lake with an area of 35.2 km2 and 25 km length in the north-south direction, as shown in Figure 6. The formation of cuspate forelands in this lake was discussed in , and here we refer this results. This lake is located in the lowland surrounded by Kashima and Namegata tablelands with elevations of 40 and 30 m on the east and west sides, respectively. Thus, wind waves can be generated without a significant sheltering effect by hills or mountains. Because the direction of the principal axis of Lake Kitaura is N18°W, the predominant wind of NNE blows at an angle of 40.5° clockwise relative to the direction of the principal axis. Because of this oblique wind direction, wind waves are incident at a large incidence angle to the shoreline, resulting in the formation of the protruding shoreline on the west shore. In particular, an enlarged satellite image of two subareas,
2.4. Lake Balkhash
Lake Balkhash has 450 and 200 km lengths in the E-W and S-N directions, respectively (Figure 8). Figure 9 shows an enlarged satellite image of the rectangular area in Figure 8. Island A is located at a location of 46°34′53.99”N and 78°50′17.47″E at the central part of the lake near the east end, and a cuspate foreland of 14 km length extends between island A and the lakeshore. On the shore opposite to island A, a triangular cuspate foreland B is formed with a barrier island. The sand bar extending northwestward from Island A is symmetric with respect to the centerline of the cuspate foreland, and the length of the sand bar is longer than the width of the island. From this, it is inferred that the cuspate foreland extended from the land to Island A by the sand supply from Island A and the land, and connected to Island A.
3. Model for predicting lakeshore changes
For the calculation of the segmentation of a rectangular water body, the BG model employed for the calculation of oriented lakes  was used. Given a local fetch distance
In this calculation, a coordinate system (
When a grid point was located on land and the downslope condition of
When the grid point was located in the lake,
Here, is the net sand transport flux,
In estimating the intensity of sand transport near the berm top and at the depth of closure, the intensity of sand transport was linearly reduced to 0 near the berm height or the depth of closure to prevent sand from being deposited in the zone higher than the berm height and the beach from being eroded in the zone deeper than the depth of closure .
4. Calculation conditions
Lakeshore changes in a rectangular water body with an aspect ratio of 5 owing to wind waves were first predicted when wind blew from all directions between 0 and 360° with the same probability of occurrence and intensity (Case 1) or blew at an angle of 45° relative to the principal axis of the rectangular water body with an elliptic probability of occurrence and intensity (Case 2), as shown in Figure 11 . Then, lakeshore changes in triangle- and crescent-shaped shallow water bodies with a flatbed were predicted in Cases 3 and 4, respectively. In all cases, the water depth of the flatbed, the berm height, and the initial beach slope were set to 3 m, 1 m, and 1/20, respectively. Figure 12 shows the initial topography in each case. Random perturbations with the amplitude △
|Wind velocity||20 m/s|
|Berm height, ||1 m|
|Depth of closure, ||3 m|
|Equilibrium slope, tan||1/20|
|Coefficient of sand transport|
|Calculation cases||Case 1: rectangular water body, circular probability distribution|
Case 2: rectangular water body, elliptic probability distribution
Case 3: segmentation of a triangular water body
Case 4: segmentation of a crescent-shaped water body
Cases 5–8: topographic changes around an island located in a circular lake
|Duration of calculation||106 h (105 steps) in Cases 1–4,|
5 × 105 h (5 × 104 steps) in Cases 5–8
|Boundary conditions||Shoreward and landward ends, |
Right and left boundaries,
In predicting lakeshore changes when a rocky or sandy island is located in a closed water body, four calculations were carried out, as shown in Figure 13. In each case, a circular lake with a radius of 1000 m and a solid bottom of a constant depth of 3 m was set for the calculation domain. In this circular lake, a rocky or sandy island with a radius of 200 m was set at locations deviating from the center of the lake. The foreshore slope of the lakeshore was assumed to be 1/20. In the present study, the incident angle of waves to the mean shoreline exceeds 45° at certain locations of the lakeshore, resulting in shoreline instability. Therefore, a small perturbation with the amplitude Δ
5.1. Segmentation of water body given circular distribution of probability (Case 1)
Figure 14 shows the calculation results for the segmentation of a slender, rectangular water body with a longshore length of 4.5 km, and a width of 0.9 km (aspect ratio = 5), assuming that the probability of occurrence of wind direction was given by a circular distribution . When wind waves were incident to the lakeshore, several cuspate forelands with irregular shapes developed along the shoreline in the beginning. After 2 × 104 steps, the cuspate forelands merged with each other, resulting in a reduction in their number, and sand bars with a hound’s-tooth shape were formed. This development of cuspate forelands well explains the formation of the lakeshore, as shown in Figure 2. After 4 × 104 steps, sand bars extended to the opposite shores, and the water body was about to separate into two lakes, and then the water body had separated into two completely independent lakes. Finally, two completely rounded lakes were formed.
The distributions of the wave height and longshore sand transport alter in response to the wind direction at each time. The formation of cuspate forelands and rounded lakes over time, however, strongly depends on the mean (
The mean sand transport flux after 4 × 104 and 5 × 104 steps in Case 1 can be drawn, as shown in Figure 16 . Intensive sand transport flux occurred along the shoreline of a cuspate foreland at the central part of the water body, enhancing further development of a cuspate foreland. Also, intensive sand transport took place near the right corner of the slender water body because of a large aspect ratio of the water body, which induced the formation of a circular lake.
5.2. Segmentation of water body given elliptic distribution of probability (Case 2)
In Case 2, wind blew from the direction of 45° with respect to the principal axis of the slender lake, that is, the probability of occurrence of the wind direction is given by an elliptic distribution . Uda et al.  predicted the formation of oriented lakes  using the BG model and showed that oriented lakes can develop when the probability of occurrence of the wind direction is given by an elliptic distribution. Here, the segmentation of a rectangular water body was predicted, assuming that the probability of occurrence was given by an elliptic distribution.
Figure 17 shows the predicted results of the lake averaged over 103 steps in Case 2 . Cuspate forelands with an asymmetric form had developed on both shores and inclined rightward (leftward) on lower (upper) shorelines in the beginning. Then, the cuspate forelands had merged to increase their size and moved rightward (leftward) on lower (upper) shorelines. These results are in good agreement with those obtained by Uda et al.  concerning the development of sand spits and cuspate forelands owing to the shoreline instability. Because the principal axis of the wind direction is at an angle of 45° relative to the shoreline, and the effect of wind blowing from the land to the lake can be neglected along lower shoreline, the oblique component of waves incident from the left had a higher probability than that of waves incident from the right. As a result, rightward sand transport predominantly caused the formation of a cuspate foreland with an asymmetric shape along the
5.3. Segmentation of a triangular or crescent-shaped water body (Cases 3 and 4)
Figure 18 shows the results of the calculation of the segmentation of a triangular water body, assuming that the probability of occurrence of wind direction was given by a circular distribution . Although the results are similar to those in , numerical simulation was carried out with changing the size of the lake because of the revision in Eq. (2b). Segmentation rapidly occurred in the vicinity of the right end of the triangular water body, and elliptic lakes were formed in the area between
Similarly, Figure 19 shows the results of the segmentation of a crescent-shaped water body with time. Rapid segmentation occurred in the vicinity of the both ends of the crescent water body in the beginning. In the area between
5.4. Lakeshore changes in circular lake with a rocky or sandy island
Figure 20 shows the lakeshore changes in Case 5 with a rocky island in a circular lake, assuming that the probability of occurrence of wind direction was given by a circular distribution. Under the condition, a wave-shelter zone was primarily formed on the lee of the island against wind waves incident from
Figure 21 shows the same results in Case 6. In this case, the wave-sheltering effect due to the island was strengthened than that in Case 5 because of the proximity of the island to the lakeshore, the cuspate foreland rapidly developed together with the formation of a large cuspate foreland. After 5 × 104 steps, a headland with a circular head was formed. Because the distance between the island and lakeshore decreased, the wave-sheltering effect increased, resulting in the greater cuspate foreland behind the island and erosion on the opposite shore.
The lakeshore changes in Case 7 with a sandy island in a circular lake are shown in Figure 22. When waves were incident to the sandy island, the island deformed by the action of waves incident from
Figure 23 shows the same results in Case 8. The initial circular island significantly deformed owing to the action of wind waves incident from the direction of
The mean sand transport fluxes averaged over 1000 steps between 1.9 × 104 + 1 and 2 × 104 steps in Cases 5 and 7 with the same arrangement of an island are shown in Figure 24. In Case 5 with a rocky island, the intensive sand transport flux occurred on both sides of the island with decreasing the intensity behind the island, whereas in Case 7, strong sand transport flux toward the tips of the sand bar occurred along the shoreline of sand bars. When setting point O at the center of the circular lake, and points
In the case of Lake Balkhash, part of sand comprised of the island was considered to be transported northwestward, forming a long slender sand bar. Such topographic changes can be explained by the mergence of a slender sand bar extended from the island and the sand bar extended in the opposite direction, as shown in the results after 5 × 104 steps in Case 5 and after 2 × 104 steps in Case 7.
Specific geomorphological features associated with shoreline instability under a high-wave-angle condition on the lakeshore, such as the development of sand spits, in several elongated water bodies were investigated, and the segmentation of a water body was numerically predicted using the BG model. It was concluded that a rectangular water body segmented into circular (elliptic) lakes when the probability of occurrence of the wind direction was given by a circular (elliptic) distribution. In each case, the wave-sheltering effect of the cuspate forelands played a primary role. Also, the mergence and segmentation of triangular and crescent-shaped slender water bodies were predicted using the BG model. It was further used for predicting the lakeshore changes when a rocky or sandy island exists in a circular lake. The deformation of a sandy island and mergence of the sandy island to the lakeshore were predicted well.
Ashton A, Murray AB, Arnault O. Formation of coastline features by large–scale instabilities induced by high-angle waves. Nature. 2001; 414:296-300
Ashton A, Murray AB High-angle wave instability and emergent shoreline shapes: 1. Modeling of sand waves, flying spits, and capes. Journal of Geophysics Research. 2006; 111:F04011. DOI: 10.1029/2005JF000422
Zenkovich VP. Processes of Coastal Development. New York: Interscience Publishers; 1967. p. 751
Ashton A, Murray AB, Littlewood R, Lewis DA, Hong P. Fetch limited self-organization of elongate water bodies. Geology. 2009; 37:187-190
Uda T, Serizawa M, Miyahara S. Numerical simulation of three-dimensional segmentation of elongated water body using BG model. Proc. 33rd ICCE, sediment. 65, 2012. pp. 1-11
Serizawa M, Uda T, San-nami T, Furuike K. Three-dimensional model for predicting beach changes based on Bagnold’s concept. Proc. 30th ICCE, 2006; pp. 3155-3167
Uda T, Serizawa M, Miyahara S, San-nami T. Prediction of segmentation and mergence of shallow water bodies by wind waves using BG model. Proc. Coastal Dynamics. 2013; Paper No. 166, pp. 1729-1740
Uda T, Serizawa M, San-nami T, Miyahara S. Prediction of formation of oriented lakes. Proc. 34th ICCE, ASCE. 2014. pp. 1-12
Scheffers AM, Kelletat DH. Lakes of the World with Google Earth: Understanding our Environment. Coastal Research Library, Vol. 16. Switzerland: Springer International Publishing; 2016. p. 293
Miyahara S, Uda T, Serizawa M, San-nami T. Prediction of effects of artificial alteration on segmentation of a slender water body. Coastal Sediments ‘15. 2015; 164:1-14
Wilson BW. Numerical prediction of ocean waves in the North Atlantic for December, 1959. Deut. Hydrogr. Zeit, Jahrgang 18, Heft 3. 1965. pp. 114-130
Goda Y. Revisiting Wilson’s formulas for simplified wind-wave prediction, Journal of Waterway. Port, Coastal and Ocean Engineering. 2003; 129(2):93-95
Serizawa M, Uda T, Miyahara S. Segmentation of water body given probability of occurrence of wind direction by circular or elliptic distribution. Proc. 35th Conf. Coastal Eng., sediment. 6. 2016. pp. 1-15
San-nami T, Uda T, Gibo M, Ishikawa T, Miyahara S, Serizawa M. Change in carbonate beach triggered by construction of a bridge on Irabu Island and its simulation using BG model. Asian and Pacific Coasts 2013, Proc. 7th International Conf. 2013. pp. 24-31
Seppälä M. Wind as a Geomorphic Agent in Cold Climates. New York: Cambridge University Press; 2004. p. 358