Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Application of GIS Spatial Analysis for the Assessment of Storm Surge Inundation Risks in the Guangdong–Macao–Hong Kong Great Bay Area
Previous Article in Journal
Evolution and Analysis of Water Yield under the Change of Land Use and Climate Change Based on the PLUS-InVEST Model: A Case Study of the Yellow River Basin in Henan Province
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interrelationship between Wall and Beach Erosion in Loc An, Vietnam: Remote Sensing and Numerical Modeling Approaches

1
Faculty of Water Resource Engineering, College of Engineering, Can Tho University, Can Tho 94115, Vietnam
2
Institute of Civil Engineering, Thuyloi University, Ha Noi 11515, Vietnam
3
Faculty of Civil Engineering, Thuyloi University, Ha Noi 11515, Vietnam
4
Institute of Liberal Arts and Sciences, Tohoku University, 41 Kawauchi, Aoba-ku, Sendai 980-8576, Japan
*
Author to whom correspondence should be addressed.
Water 2024, 16(17), 2553; https://doi.org/10.3390/w16172553
Submission received: 21 July 2024 / Revised: 5 September 2024 / Accepted: 7 September 2024 / Published: 9 September 2024
(This article belongs to the Section Oceans and Coastal Zones)

Abstract

:
Beach erosion and coastal protection are complex and interconnected phenomena that have a substantial impact on coastal environments worldwide. Among the various coastal protection measures, seawalls have been widely implemented to mitigate erosion and protect coastal assets. However, the interrelationship between beach erosion and seawalls remains a critical topic for investigation to ensure effective and sustainable coastal management strategies. Seawalls impact the shoreline, particularly through the “end effect”, where the seawall functions similarly to a groin, causing erosion on the downdrift side relative to the direction of wave approach. This study provides a detailed analysis of the interplay between beach erosion and seawall structures in Loc An, Vietnam, employing both remote sensing and numerical approaches. Sentinel-2 images were employed together with an analytical solution to observe the shoreline change at the Loc An sand spit and to determine input values for the numerical model. Based on the shoreline dynamics, a numerical scheme was employed to study the shoreline evolution after the construction of a seawall. Our findings show that the shoreline evolution can be divided into three stages: (1) The first stage corresponds to the elongation of the sand spit without interference from coastal structures. (2) The second stage shows the effect of jetties on the shoreline, as signaled by the buildup of sand updrift of the jetties. (3) The third stage shows the effectiveness of the seawall, where the shoreline reaches its equilibrium condition. The study provides a quick and simple method for estimating shoreline diffusivity (ε) in situations where measured data is scarce.

1. Introduction

Coastal areas are among the most dynamic and vital regions on Earth, serving as the frontier between land and sea [1]. These zones are not only crucial for biodiversity and ecosystem services but also support a significant portion of the global population through tourism, fisheries, and as gateways for trade and transportation. Among these environments, river mouths and inlets hold a special place because of their ecological, economic, and social importance. These areas experience substantial changes in morphological features, driven by a complex interplay of natural factors such as typhoons [2], river flows [3,4], tides [5,6], waves [7,8], and longshore currents [9,10]. The constant reshaping of these landscapes poses challenges for their physical stability and the communities that depend on them.
To mitigate these challenges and stabilize river mouths and inlets, various solutions have been implemented [11,12,13], with the construction of coastal structures like jetties and seawalls being prominent. While the impact of these interventions on shoreline dynamics has been extensively studied globally [14,15,16,17,18], research in Vietnam remains scarce [19]. This gap is primarily due to the lack of measured data, which is crucial for understanding and predicting coastal behavior under the influence of both natural events and human interventions.
With the advancement of remote sensing technology, monitoring and measuring the coastal environment has become easier than ever. This modern technique provides an effective tool for studying the coastal zone, offering coverage at various resolutions and regular temporal frequencies [20]. In recent years, remote sensing has been increasingly used in a wide range of ocean engineering applications, including the monitoring of coastal features such as barrier islands [21], shoreline dynamics [22], changes in mangrove forests [23], and even floating litter on the sea surface [24]. Furthermore, novel approaches have been introduced to enhance coastal applications [25,26,27].
This article aims to explore the effectiveness of combining traditional models (such as shoreline modeling) with modern techniques (like remote sensing) in coastal studies, with a specific focus on the Loc An sand spit at the Loc An estuary in Ba Ria-Vung Tau province (Figure 1). In this study, the data extracted from Sentinel-2 images using remote sensing techniques will be used to calibrate and validate the simple shoreline model. Despite the limited availability of field data, this study combines remote sensing techniques with a simple numerical model to explore the influence of coastal structures on shoreline evolution in regions where direct measurements are scarce. This research adds to the understanding of coastal management and demonstrates the potential of integrating multiple scientific methods to overcome data limitations, offering useful insights for similar studies in other parts of the world.
The study area is influenced by prevailing northeast winds, resulting in a net longshore sediment transport (LST) from northeast to southwest. This has led to the elongation of the Loc An sand spit and a tendency for the river mouth to shift southwestward. The waves, predominantly from the east-northeast direction, with significant heights ranging from 1.51 m to 3.0 m, generate longshore currents that drive sediment transport along the coast in the northeast-to-southwest direction [5]. As the Loc An estuary is located along the East Sea of Vietnam, the semi-diurnal tide also plays a role in the morphological changes of the estuary. Information on the tidal conditions is shown in Figure 2. Here, the water level is measured at the Vung Tau Oceanographic Station, which is located approximately 30 km southwest of the Loc An estuary (see Figure 1).

2. Materials and Methods

2.1. Remote Sensing

Launched in 2015 and 2017, respectively, Sentinel 2A and 2B are twin satellites orbiting the poles, established under the European Copernicus initiative [28]. Their purpose is to perform systematic observations of the atmosphere, land, and oceans. These satellites are equipped with a multispectral imager (MSI) that captures imagery in 13 spectral bands ranging from 443 nm to 2190 nm, offering broad coverage with a high spatial resolution. The resolution varies, with most bands in the visible and very near infrared (VNIR) spectrum offering a 10 m resolution, while the short-wave infrared (SWIR) bands have a 60 m resolution. Sentinel 2A and 2B boast a temporal resolution that can reach up to every 5 days, although this varies with latitude, decreasing as one moves further from the equator [29]. This capability enables frequent monitoring and observation of the earth’s land, atmosphere, and oceans.
In this study, Sentinel-2 images covering the period from 2017 to 2023, which depict the sand spit at the Loc An River mouth (Figure 1), were selected for analysis. The details of the Sentinel-2 images are presented in Table 1.
Shoreline detection was performed based on the Normalized Difference Water Index (NDWI) [30] as follows:
N D W I = X g r e e n X n i r X g r e e n + X n i r
where Xgreen and Xnir are the GREEN (B3) and NIR (B8) spectral bands with the central wavelengths of 560 nm and 842 nm, respectively. The extracted shorelines were mapped together to observe the changes in the shoreline. Since the position of the shoreline is influenced by tidal variations, which are particularly important in meso- and macro-tidal regions, tidal correction was applied to the shoreline positions using the method presented in Hoang et al. [31]. Hourly water levels collected at the Vung Tau Oceanographic Station from 2015 to 2023 were utilized. The beach slope was also obtained from the bathymetric map in the topography survey reported by Royal Haskoning DHV for the coastline at Ho Tram Resort, which is 10 km northeast of the Loc An coast [32].

2.2. One-Line Model

The analytical method to analyze the development of shoreline changes along the Loc An sand spit was based on the concept of the one-line model, as depicted in Figure 3. This approach operates under the premise of the one-line theory, which suggests that the beach maintains a consistent profile, implying that all underwater contours remain parallel to each other. Hence, a single line can represent seaward or landward changes of the shoreline, with this line directly aligning with the shoreline itself, as shown in Figure 3. Figure 3 implies that the shoreline change is determined by the net rate of sand entering and leaving the beach segment over a given time down to a certain depth as
Δ y Δ t = 1 D B + D C Δ Q Δ x ± q
where Δyt is the shoreline change; DB is the berm height and DC is the depth of closure, respectively; ΔQt is the net rate of sand entering and leaving the beach segment (Δx); and q is the source (q0) of sink (qs) of sand along the segment (Figure 3).
The morphological changes of the shoreline’s shape, denoted as y(x,t), are governed by the One-Line Contour Equation (OLCE) [33], a formula first introduced by Pelnard-Considére in 1956 [34]. According to Larson et al. (1987) [35], under conditions of a small breaking wave angle and slight shoreline curvature, the OLCE simplifies into the diffusivity equation (DSE). The DSE is expressed as
y t = ε 2 y x 2
In this context, ε represents the diffusivity of the shoreline, determining the rate at which the system reacts to external influences. The value of ε is primarily influenced by the height and angle of the breaking waves, as highlighted in various studies within the literature, including works by Ashton and Murray (2006 a,b) [36,37], Falques (2003, 2005) [38,39], and Walton and Dean (2010) [40].
It should be noted that various analytical solutions have been derived for the DSE as presented in books [41,42]. Therefore, by applying appropriate boundary conditions for shoreline evolution, numerous analytical solutions for shoreline change can be derived from Equation (3), which primarily fall into two categories: infinite shorelines and finite shorelines. Those in the first category were studied and introduced by Pelnard Considère for the case of shoreline evolution downdrift of a groin [34] or shoreline evolution of a river delta with infinite shorelines [35]. In the second category, Tanaka and Nadaoka [17] introduced an analytical solution for shoreline evolution bounded by two groins and Duy et al. [19] presented an analytical solution for the formation and deformation of a river delta with finite beach length.

3. Results of Satellite Image Analysis

3.1. Stage 1: Elongation of the Sand Spit

As seen in Figure 4, there are three stages of shoreline evolution. In the first stage, the shoreline of the Loc An sand spit was in its natural condition, as there were no structures on the shoreline (Figure 4a). The shorelines on 5 December 2015 and 14 January 2016 were also added to Figure 4a to highlight the elongation of the sand spit due to LST. This stage of the Loc An sand spit was detailed in the study by Duc Anh et al. [18]. According to Duc Anh et al. [18], the Loc An sand spit was extended westward about 2500 m from 1988 to 2017 by the longshore sediment transport at a rate of 2 × 105 m3/year from northeast to southwest. This direction of LST coincides with the dominant northeast wave direction.

3.2. Stage 2: Intercept of LST by a Jetty

In the second stage, the jetties were constructed (Figure 4b) to stabilize the river mouth. At this stage, the LST was intercepted by the eastward jetty, leading to the accumulation of sand and the seaward advance of the shoreline adjacent to the eastward side of the jetty. According to satellite images, the eastward jetty was constructed in 2017, and the westward jetty was constructed in 2018. The construction of the jetties has caused a change in the direction of the river flow from westward (1988 to 2017) to southward (2018–present). This directional change in river flow can cause sand to bypass to the east side of the jetty and contribute to the buildup of sediment immediately eastward of the jetties.

3.3. Stage 3: Construction of the Seawall and LST in between Two Boundaries

Finally, the shoreline was bounded by both a jetty and a seawall, as seen in Figure 4c. The seawall is clearly visible in Figure 4d, which was taken during June 2022. It should be noted that a groin was constructed eastward of the seawall using geo-tubes to prevent beach erosion along the Loc An sand spit. However, this groin was destroyed within a few months of its construction. As a result, the seawall was built to stabilize the shoreline (Figure 4d).
This seawall was constructed using a prefabricated non-metallic reinforced concrete solution. The revetment is assembled using A-shaped non-metallic reinforced concrete blocks, each 4 m high, 4.1 m wide, and either 2 m or 1.5 m long. These blocks are connected using male and female grooves to form the framework of the seawall. The M300 concrete wave deflector wall is attached to the A-shaped blocks with M300 concrete square piles, measuring 25 cm × 25 cm and 9 m long. Behind the revetment, the earthfill is raised to an elevation of +3.5 m with continuous ART 11D double-layer geotextile. The foot of the revetment is protected against erosion with 60 cm thick rubble stone, 4.2 m wide, placed on two layers of ART 11D geotextile and combined with M300 concrete-coated gabions, 0.8 m wide and 1 m high, which are durable in seawater environments. The A-shaped concrete blocks are connected using vertical sliding locks with a lock depth greater than 2 cm.
Because the seawall at the Loc An sand spit is a semi-infinite seawall where flanking (erosion occurring beside the seawall) is not possible, this seawall can function similarly to a groin or a jetty [35].
The shoreline changes of the Loc An sand spit from 2015–2023 are presented in Figure 5. As can be seen in Figure 5, the shoreline evolution of Loc An can be divided into three stages. The first stage corresponds to the elongation of the sand spit because of longshore sediment transport with no constraint from coastal structures. This stage includes the shorelines on 17 January 2015, 20 January 2016, and 6 January 2017.
In the second stage, the shoreline was affected by the construction of jetties around 31 August 2017. The influence of the jetties became clear from the buildup of the shorelines on 7 February 2018 and 13 January 2019.
The third stage includes the shorelines on 16 February 2022 and 13 March 2023. In this stage, a seawall had been constructed, hence the shorelines were affected by both the jetties and the seawall. Interestingly, there was no significant change in the shoreline position during this stage, showing that the shoreline had approached its equilibrium condition. This phenomenon agrees well with the analytical solution proposed by Tanaka and Nadaoka (1982) [17].

4. Application of One-Line Model

The three stages of shoreline evolution at the Loc An sand spit can be sketched as in Figure 6. In this section, two analytical solutions for shoreline changes, derived from a well-known theoretical one-line model, were used to examine the sandy beach evolution effected by one boundary (stage 2) and two boundaries (stage 3) in Figure 6.

4.1. Solution for Sandy Beach with One Boundary

Beginning with Pelnard-Considere’s initial work in 1956 [34], which presented an analytical solution for shoreline change on the updrift side of a groin, numerous researchers have since developed analytical tools for a variety of applications. A key study of interest for this research was conducted by Larson et al. (1987) [35]. Their study specifically addressed the open (unbounded) coast hypothesis, significantly simplifying the mathematical challenge. Within the range of analytical solutions offered by Larson et al. (1987) [35], they included a scenario addressing shoreline change in the presence of groins and jetties (Figure 7).
As seen in Figure 7, a groin disrupts the alongshore movement of sand, resulting in a build-up on the side facing the incoming waves (updrift side) and a depletion or erosion on the opposite side (downdrift side). The equation characterizing the part where accumulation occurs is
y ( x , t ) = 2 tan θ ε t π e x x 0 2 / 4 ε t x x 0 2 e r f c x x 0 2 ε t
where θ is the angle of breaking wave crests relative to an axis set parallel to the trend of the shoreline, t is the time, and erfc is the complementary error function.
As shown in Figure 4b, the shoreline change at this stage (stage 2) is markedly influenced by the jetties functioning as a groin. Hence, Equation (4) can be applied to model the shoreline’s evolution during this phase. The shoreline data extracted from Sentinel-2 images at this stage were compared with the modeled shoreline to ascertain the values of ε and θ.
To facilitate a comparison between the shorelines extracted from Sentinel-2 images and the modeled shorelines (Equation (4)), the images were converted into local coordinates as depicted in Figure 8.
The idea for determining ε and θ is to make a comparison between the modeled shoreline and the measured shoreline (extracted from Sentinel-2 images). The modeled shoreline is obtained using Equation (4). However, Equation (4) requires the values of ε and θ for the calculation. Therefore, iteration was performed to determine ε and θ. To evaluate the agreement between the modeled shoreline and the extracted one, the root mean square error (RMSE) was used as in the following equation:
R M S E = y p i y o i 2 n
where ypi is the predicted shoreline position at point i and yoi is the observed shoreline position at point i along the shoreline, and n is the number of observation points along the shoreline.
The iteration process will stop when the minimum RMSE is obtained. The initial values of the parameters required in Equation (4) are shown in Table 2.
At this stage, the shoreline extracted from the Sentinel-2 image dated 27 July 2019 was selected as the measured shoreline. The reason for choosing the shoreline on this day, was that the seawall had not yet been constructed, allowing Equation (4) to apply. It is clear in Figure 9a that the construction of the jetties began around 31 August 2017. Therefore, the shoreline as of 31 August 2017 is considered the initial shoreline. The modeling time can be determined as the duration between 31 August 2017 and 27 July 2019, or t = 1.9 years.
Figure 10 shows the comparison between the measured shoreline and the computed shorelines, whereby the measured shoreline is extracted from the Sentinel-2 image captured on 27 July 2019, and the computed shorelines were obtained using Equation (4). As illustrated in Figure 10, several values of ε were used in the calculations, and a value of ε = 230 m2/day demonstrates good agreement with the measured shoreline. Furthermore, the RMSE values presented in Figure 11 confirm that the best value of ε is 230 m2/day.
At the same time during the iteration to determine the value of ε, various values of θ were also input into the model for calculation. The final values of ε and θ are shown in Table 3.
In order to clarify the suitability of Equation (4) for solving the problem at the Loc An sand spit, the model was applied with another period while still keeping the values of ε = 230 m2/day and θ = 18 degrees. At this stage, the shoreline position on 26 July 2020 was used. Hence, the modeling time is determined as the duration between 31 August 2017 and 26 July 2020, or t = 2.9 years. The validation results are shown in Figure 12. It can be seen that the RMSE in the validation step is 8.9 m which is quite close to the RMSE in the calibration step. This indicates the good performance of the model.

4.2. Solution for Sandy Beach with Two Boundaries

As illustrated in Figure 4c (Stage 3), the shoreline evolution at this stage is influenced by two structures: the jetties and the seawall. Equation (4) was no longer valid and a new analytical solution had to be employed as in the study of [19]. Tanaka and Nadaoka (1982) [17] introduced an analytical solution for shoreline evolution between the two jetties as sketched in Figure 13.
Figure 13 shows that the analytical solution proposed by Tanaka and Nadaoka (1982) [17] agrees well with the case of Loc An shoreline in stage 3, bounded by the jetties and the seawall. Therefore, the analytical solution of Tanaka and Nadaoka (1982) [17] can be used to study the shoreline evolution of the Loc An sand spit in stage 3, and the solution is presented as follows:
y = L . tan θ . n = 1 1 n . 4 ε t π L 2 . exp L 2 4 ε t x L + 1 2 n 2 + x L + 1 2 n . erfc L 4 ε t x L + 1 2 n + L . tan θ . n = 0 1 n . 4 ε t π L 2 . exp L 2 4 ε t x L + 1 2 n 2 x L + 1 2 n . erfc L 4 ε t x L + 1 2 n
where L is the length of the shoreline bounded by the jetties and the seawall, L = 960 m (Figure 8), and x′ is defined as
x = x x 1
The temporal variation of the shoreline upstream of the jetties, denoted as y1 in Figure 6, will be examined by plotting it under two scenarios: one with a single boundary (the jetties) and the other with two boundaries (the jetties and the seawall). Furthermore, we utilized the dimensionless form of the shorelines to establish the criteria for when the shoreline approaches its equilibrium. The dimensionless shoreline position and the dimensionless time are defined as follows [17]:
η 1 = 2 y 1 L tan θ
where η1 represents the dimensionless form of the shoreline, y1 is the position of the shoreline upstream of the jetties, and L is the distance between the jetties and the seawall (as shown in Figure 8).
τ = 4 ε L 2 t
where τ is the dimensionless time.
The dimensionless shoreline positions are plotted in Figure 13 with the normal axis in Figure 14a and the log-log plot in Figure 14b. As can be seen in Figure 14, in the initial stage of the shoreline evolution (τ < 0.5), there is good agreement between the shoreline positions with one boundary and two boundaries. This indicates there is no effect of the second boundary (the seawall in the case of Loc An beach). However, for τ > 0.5, there is a significant gap between the shoreline positions of one boundary and two boundaries, indicating that the longshore sediment transport has reached the second boundary and the shoreline changes were affected by it. It is also observed from Figure 14 that after τ = 2.0, the shoreline position affected by two boundaries will reach the equilibrium condition, as indicated by a horizontal line. This phenomenon is similar to the study of Duy et al. [19].

4.3. Comparison with Field Measured Data

In this section, the shoreline positions in stage 3, modeled by Equation (6), will be compared with the shoreline positions extracted from the satellite images to determine whether the shoreline position at Loc An beach has reached its equilibrium condition.
Equation (6) requires a straight shoreline as the initial condition for the calculation (Figure 13). This condition cannot be met with the Loc An sand spit because the shoreline at Stage 3 had already been influenced by the jetties and was no longer straight. Hence, a numerical approach was employed to model the shoreline evolution with an arbitrary initial shoreline.
In order to apply the numerical model, an initial shoreline is required. This initial shoreline must be established after the construction of the seawall to ensure that the shoreline is affected by both the jetties and the seawall. Based on a photo taken by a contractor, the seawall was constructed around 25 September 2021 (Figure 15b). Additionally, a Sentinel-2 image captured on 29 October 2021 also showed the existence of the seawall, indicated by a straight shoreline (Figure 15b). Therefore, the shoreline position on 29 October 2021 was chosen as the initial shoreline for the numerical model.
The comparison between the modeled shorelines and the measured shorelines is shown in Figure 16. In this figure, the modeled shorelines were obtained from the numerical model by discretizing Equation (6). The measured shorelines were extracted from Sentinel-2 images. It should be noted that the shoreline positions presented in Figure 16 were transformed to local coordinates (Figure 8). The locations of the seawall and the jetties, as well as the direction of prevailing offshore waves, are also presented in Figure 16. Interestingly, the shoreline position appears to be stable after 2 years of the seawall’s construction, especially at the location close to the jetties. This shows the effectiveness of the seawall in stabilizing the shoreline along the Loc An sand spit. At the location down-drift of the seawall, a disagreement between the modeled and measured shorelines can be observed, which is attributed to the buildup of the shoreline. This discrepancy is caused by wave refraction, which was not considered in the numerical model.
The evolution of the shoreline at the jetties (x = 960 m) is plotted in Figure 17 to highlight the equilibrium of the shoreline achieved some time after the construction of the seawall (see Figure 6 for a graphical illustration of y1). In this Figure, the numerical shorelines are represented by a blue line, and the measured shoreline positions are denoted by circles. Although there are fluctuations in the measured shorelines, they follow the trend of the numerical calculations. The numerical results predict that, after 4 years of the seawall’s construction, the shoreline at the Loc An sand spit will stabilize. This underscores the effectiveness of the seawall.
It should be noted that the equilibrium shoreline position, denoted by the dashed red line, is also presented in the figure. Based on the equilibrium shoreline position and the numerical shorelines, it can be concluded that the shoreline bounded by the jetties and the seawall at the Loc An sand spit will be stable by mid-2027. This equilibrium condition of the shoreline at the Loc An sand spit is based on the assumption that sand bypass is neglected. However, it is clear that sand bypass at the eastern jetty and sand transport along the seawall are inevitable. This phenomenon is disregarded in the simple model of this study to quickly analyze basic beach behavior with an acceptable level of accuracy.

5. Discussion

5.1. Tidal Effect to the Shoreline Positions

As shown in Figure 2, the tidal range is approximately 3 m, which can significantly affect the shoreline position as determined from satellite images. Consequently, a detailed discussion of the impact of tides on shoreline shifting would be useful. Figure 18 illustrates the shoreline positions at the jetties (y1) before and after tidal corrections are applied. It is evident that the maximum shoreline shift due to tidal effects at the Loc An river mouth is approximately 20 m (on 12 April 2022).

5.2. Determination of ε Value

In this section, the value of ε is calculated based on measured data to discuss the value of ε determined by the iteration shown in the previous sections. According to Rosati et al. (2002) [43], the value of ε is calculated using the following equation:
ε = K H 2 C g b 8 ρ ρ s ρ 1 1 n 1 D B + D C
Here, K represents the dimensionless coefficient for the sediment transport rate formula, H is the wave height, and Cg is the group velocity. The subscript “b” denotes the quantity at the wave-breaking point. ρ and ρs represent the density of sea water and sand, ρ = 1025 kg/m3 and ρs = 2650 kg/m3, respectively, n denotes the sediment porosity (n = 0.4), DB is the berm height, and DC is the depth of closure.
According to Del Valle et al. (1993) [44], K varies with the median grain size (D50) of the sediment and can be calculated using the empirical relationship
K = 1.4 e 2.5 D 50
According to a survey map from Duc Anh (2021) [5], the median grain size of sand in the study area (Loc An sand spit) ranges from 0.56 to 0.85 mm (Figure 19). By taking the interval of 0.05 mm between 0.56 mm and 0.85 mm, the values of K can be determined as in Table 4.
The breaking wave height (Hb) was calculated based on the measured waves from 18 September 2016 to 9 October 2016 at the location as shown by the wave buoy in Figure 15. The automatic method for calculating breaking waves proposed by Sana [45] was employed and this resulted in Hb = 0.54 m.
The group velocity Cgb is calculated based on the following equation:
C g b = g H b κ
where the value of the breaker index, κ, is often taken as 0.78 [46].
The depth of closure, DC, was determined from the bathymetric map in the topography survey reported by Royal Haskoning DHV for the coastline at Ho Tram resort which is 10 km North-East of the Loc An coast. According to the report of Royal Haskoning DHV, the value of DC is 6.5 m [32].
From the value of DC, the berm height (DB) can be determined as follows:
D B = 0.32 D C = 2.1   ( m )
By substituting all the values into Equation (10), we obtain the values of ε as in Table 5.
Figure 20 shows the comparison between the values of ε calculated by two methods. The blue bars represent the variation of ε calculated by Equation (5) based on measured data, while the red line shows the value of ε determined by the simple method (see Figure 11). The result from the simple method falls within the upper and lower bounds of the ε value range. Hence, the simple method can determine the value of ε in cases of data scarcity.

6. Conclusions

Shoreline evolution of a sand spit affected by coastal structures at Loc An, Ba Ria–Vung Tau province, Vietnam was investigated by combining remote sensing techniques and a numerical model. In this approach, a numerical method was applied based on a simple model that can be quickly used to determine the diffusion coefficient (ε) for LST, which accounts for wave energy and sediment grain size on the beach. Although the analytical solutions are based on several approximations, they have been proven effective in numerous study areas worldwide, particularly where field measurement data are limited. The results of this study can support the sustainable development of the coastal area using placed structures as erosion control measure. The main findings of this study can be summarized as follows:
-
The evolution of the Loc An sand spit was divided into three stages in which the first stage showed the elongation of the sand spit because of LST without human intervention. The second stage showed the shoreline buildup upstream of the jetties. In the third stage, the shoreline at the Loc An sand spit was affected by both the jetties and the seawall and the shoreline approached equilibrium conditions.
-
The value of shoreline diffusivity and the breaking wave angle at the Loc An coastline were determined as 230 m2/day and 18 degrees, respectively.
-
Although erosion occurred after the construction of the seawall, this erosion was limited and the shoreline still progressed towards equilibrium conditions.
-
A comparison between the measured and calculated data of shoreline evolution at the jetty was conducted. Although there were some fluctuations in the measured data, they generally followed the trend of the calculated shoreline evolution. This indicates the reliability of the simple shoreline model used in this study.
-
There is a good agreement between the simple method proposed in this study and the measured data for calculating ε. Hence, the method can determine the value of ε for data-scarce localities.
-
If the sand bypass at the eastern jetty and sediment transport along the seawall are neglected, the shoreline at the Loc An sand spit is expected to reach its equilibrium condition by mid-2027.

Author Contributions

Conceptualization, H.T., N.T.V. and D.V.D.; methodology, H.T., N.T.V. and D.V.D.; software, H.T., D.V.D. and N.Q.D.A.; validation, H.T., D.V.D. and N.Q.D.A.; formal analysis, H.T., N.T.V., D.V.D. and N.Q.D.A.; investigation, H.T. and D.V.D.; resources, H.T., N.Q.D.A. and D.V.D.; data curation, H.T., D.V.D. and N.T.V.; writing—original draft preparation, D.V.D.; writing—review and editing, H.T., N.T.V. and D.V.D.; visualization, H.T., D.V.D. and N.T.V.; supervision, H.T. and N.T.V.; funding acquisition, H.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Datasets used in this study are contained within the article. Further data requests should be addressed to the corresponding author or [email protected].

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Crossland, C.; Kremer, H. Coastal zones: Ecosystems under pressure. In Proceedings of the Oceans and Coasts at RIO+ 10-Global Conference, Paris, France, 3–7 December 2001; UNESCO: Paris, France, 2001. [Google Scholar]
  2. Almar, R.; Marchesiello, P.; Almeida, L.P.; Thuan, D.H.; Tanaka, H.; Viet, N.T. Shoreline response to a sequence of typhoon and monsoon events. Water 2017, 9, 364. [Google Scholar] [CrossRef]
  3. Viet, N.T.; Tanaka, H. Study on the effect of morphology change on salinity distribution in the Dinh An estuary, lower Mekong River of Vietnam. J. Coast. Res. 2007, SI 50, 268–272. [Google Scholar] [CrossRef]
  4. Duy, D.V.; Ty, T.V.; Than, C.T.N.; Thang, C.N.; Hong, H.T.C.; Viet, N.T.; Tanaka, H. Coastal Erosion Caused by River Mouth Migration on a Cuspate Delta: An Example from Thanh Hoa, Vietnam. Hydrology 2023, 10, 189. [Google Scholar] [CrossRef]
  5. Duc Anh, N.Q. Comprehensive Study on the Sand Spit Elongation at Tidal Inlets in Vietnam. Ph.D. Dissertation, Tohoku University, Sendai, Japan, 2021; p. 184. [Google Scholar]
  6. Quang, D.N.; Ngan, V.H.; Tam, H.S.; Viet, N.T.; Tinh, N.X.; Tanaka, H. Long-Term Shoreline Evolution Using DSAS Technique: A Case Study of Quang Nam Province, Vietnam. J. Mar. Sci. Eng. 2021, 9, 1124. [Google Scholar] [CrossRef]
  7. Tinh, N.X.; Thanh, T.M.; Tanaka, H.; Viet, N.T.; Mitobe, Y.; Dien, D.C. Numerical investigation of the effect of seasonal variations of depth-of-closure on shoreline evolution. Int. J. Sediment Res. 2021, 36, 1–16. [Google Scholar]
  8. Thanh, T.M.; Tanaka, H.; Mitobe, Y.; Viet, N.T.; Almar, R. Seasonal variation of morphology and sediment movement on nha trang coast, vietnam. J. Coast. Res. 2018, 81, 22–31. [Google Scholar]
  9. Duy, D.V.; Tanaka, H.; Mitobe, Y.; Anh, N.Q.; Viet, N.T. Sand spit elongation and sediment balance at cua lo inlet in central Vietnam. J. Coast. Res. 2018, 81, 32–39. [Google Scholar]
  10. Duy, D.V.; Ty, T.V.; Thanh, T.N.; Minh, H.V.T.; De, C.V.; Duong, V.H.T.; Dan, T.C.; Viet, N.T.; Tanaka, H. Sand Spit Morphology at an Inlet on Phu Quoc Island, Vietnam. Water 2023, 15, 1941. [Google Scholar] [CrossRef]
  11. António, M.H.; Fernandes, E.H.; Muelbert, J.H. Impact of Jetty Configuration Changes on the Hydrodynamics of the Subtropical Patos Lagoon Estuary, Brazil. Water 2020, 12, 3197. [Google Scholar] [CrossRef]
  12. Liang, T.-Y.; Chang, C.-H.; Hsiao, S.-C.; Huang, W.-P.; Chang, T.-Y.; Guo, W.-D.; Liu, C.-H.; Ho, J.-Y.; Chen, W.-B. On-site investigations of coastal erosion and accretion for the northeast of Taiwan. J. Mar. Sci. Eng. 2022, 10, 282. [Google Scholar] [CrossRef]
  13. Siemes, R.W.A.; Duong, T.M.; Willemsen, P.W.J.M.; Borsje, B.W.; Hulscher, S.J.M.H. Morphological Response of a Highly Engineered Estuary to Altering Channel Depth and Restoring Wetlands. J. Mar. Sci. Eng. 2023, 11, 2150. [Google Scholar] [CrossRef]
  14. Tanaka, H.; Shuto, N. Field investigation at a mouth of small river. In Coastal Engineering 1992; American Society of Civil Engineers: Reston, VA, USA, 1993; pp. 2486–2499. [Google Scholar]
  15. Tanaka, H.; Lee, H.-S. Influence of jetty construction on morphology and wave set-up at a river mouth. Coast. Eng. J. 2003, 45, 659–683. [Google Scholar] [CrossRef]
  16. Lawson, S.K.; Tanaka, H.; Udo, K.; Hiep, N.T.; Tinh, N.X. Morphodynamics and Evolution of Estuarine Sandspits along the Bight of Benin Coast, West Africa. Water 2021, 13, 2977. [Google Scholar] [CrossRef]
  17. Tanaka, N.; Nadaoka, K. Development and Application of a Numerical Model for the Prediction of Shoreline Changes; Tech. Note 436; Port and Harbour Research Institute, Ministry of Transport: Yokosuka, Japan, 1982; p. 40. (In Japanese)
  18. Duc Anh, N.Q.; Duy, D.V.; Tanaka, H.; Tinh, N.X.; Viet, N.T. Elongation of sand spit at Loc An river mouth, Southern Vietnam. J. Jpn. Soc. Civ. Eng. Ser. B3 (Ocean. Eng.) 2018, 74, I_695–I_700. [Google Scholar]
  19. Duy, D.V.; Tanaka, H.; Larson, M.; Viet, N.T. A Theory for Estuarine Delta Formation with Finite Beach Length under Sediment Supplied from the River. J. Mar. Sci. Eng. 2022, 10, 947. [Google Scholar] [CrossRef]
  20. Mishra, D.R.; Gould, R.W. Preface: Remote Sensing in Coastal Environments. Remote Sens. 2016, 8, 665. [Google Scholar] [CrossRef]
  21. Waqas, M.; Nazeer, M.; Shahzad, M.I.; Zia, I. Spatial and temporal variability of open-ocean barrier islands along the Indus Delta region. Remote Sens. 2019, 11, 437. [Google Scholar] [CrossRef]
  22. Zanutta, A.; Lambertini, A.; Vittuari, L. UAV Photogrammetry and Ground Surveys as a Mapping Tool for Quickly Monitoring Shoreline and Beach Changes. J. Mar. Sci. Eng. 2020, 8, 52. [Google Scholar] [CrossRef]
  23. Pham, T.D.; Le, N.N.; Ha, N.T.; Nguyen, L.V.; Xia, J.; Yokoya, N.; To, T.T.; Trinh, H.X.; Kieu, L.Q.; Takeuchi, W. Estimating mangrove above-ground biomass using extreme gradient boosting decision trees algorithm with fused sentinel-2 and ALOS-2 PALSAR-2 data in can Gio biosphere reserve, Vietnam. Remote Sens. 2020, 12, 777. [Google Scholar] [CrossRef]
  24. Papageorgiou, D.; Topouzelis, K.; Suaria, G.; Aliani, S.; Corradi, P. Sentinel-2 Detection of Floating Marine Litter Targets with Partial Spectral Unmixing and Spectral Comparison with Other Floating Materials (Plastic Litter Project 2021). Remote Sens. 2022, 14, 5997. [Google Scholar] [CrossRef]
  25. Abdelhady, H.U.; Troy, C.D.; Habib, A.; Manish, R. A Simple, Fully Automated Shoreline Detection Algorithm for High-Resolution Multi-Spectral Imagery. Remote Sens. 2022, 14, 557. [Google Scholar] [CrossRef]
  26. Vos, K.; Splinter, K.D.; Palomar-Vázquez, J.; Pardo-Pascual, J.E.; Almonacid-Caballer, J.; Cabezas-Rabadán, C.; Kras, E.C.; Luijendijk, A.P.; Calkoen, F.; Almeida, L.P.; et al. Benchmarking satellite-derived shoreline mapping algorithms. Commun. Earth Environ. 2023, 4, 345. [Google Scholar] [CrossRef]
  27. Zollini, S.; Dominici, D.; Alicandro, M.; Cuevas-González, M.; Angelats, E.; Ribas, F.; Simarro, G. New Methodology for Shoreline Extraction Using Optical and Radar (SAR) Satellite Imagery. J. Mar. Sci. Eng. 2023, 11, 627. [Google Scholar] [CrossRef]
  28. Baba, M.W.; Thoumyre, G.; Bergsma, E.W.J.; Daly, C.J.; Almar, R. Deriving Large-Scale Coastal Bathymetry from Sentinel-2 Images Using an HIGH-Performance Cluster: A Case Study Covering North Africa’s Coastal Zone. Sensors 2021, 21, 7006. [Google Scholar] [CrossRef]
  29. Bergsma, E.W.; Almar, R. Coastal coverage of ESA’Sentinel 2 mission. Adv. Space Res. 2020, 65, 2636–2644. [Google Scholar] [CrossRef]
  30. McFeeters, S.K. The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features. Int. J. Remote Sens. 1996, 17, 1425–1432. [Google Scholar] [CrossRef]
  31. Hoang, V.C.; Tanaka, H.; Mitobe, Y. A method for correcting tidal effect on shoreline position extracted from an image with unknown capture time. Geosciences 2017, 7, 62. [Google Scholar] [CrossRef]
  32. Royal_Haskoning_DHV. Coastal Erosion Protection—Ho Tram Strip. In Coastal Erosion Protection Concept Study; Ho Tram Project Company Ltd.: Ba Ria-Vung Tau, Vietnam, 2018; p. 72. [Google Scholar]
  33. Ciccaglione, M.C.; Buccino, M.; Calabrese, M. Beaches in a Semi-Insulated Compartment: Engineering Tools from the Diffusion Theory. Estuar. Coast. Shelf Sci. 2014, 301, 108726. [Google Scholar] [CrossRef]
  34. Pelnard-Considère, R. Essai de theorie de l’evolution des formes de rivage en plages de sable et de galets. Journées de L’hydraulique 1957, 4, 289–298. [Google Scholar]
  35. Larson, M.; Hanson, H.; Kraus, N.C. Analytical Solutions of the One-Line Model of Shoreline Change; Coastal Engineering Research Center: Vicksburg, MS, USA, 1987. [Google Scholar]
  36. Ashton, A.D.; Murray, A.B. High-angle wave instability and emergent shoreline shapes: 1. Modeling of sand waves, flying spits, and capes. J. Geophys. Res. Earth Surf. 2006, 111, F04011. [Google Scholar] [CrossRef]
  37. Ashton, A.D.; Murray, A.B. High-angle wave instability and emergent shoreline shapes: 2. Wave climate analysis and comparisons to nature. J. Geophys. Res. Earth Surf. 2006, 111, F04012. [Google Scholar] [CrossRef]
  38. Falqués, A. On the diffusivity in coastline dynamics. Geophys. Res. Lett. 2003, 30, 2119. [Google Scholar] [CrossRef]
  39. Falqués, A.; Calvete, D. Large-scale dynamics of sandy coastlines: Diffusivity and instability. J. Geophys. Res. Ocean. 2005, 110, C03007. [Google Scholar] [CrossRef]
  40. Walton, T.L., Jr.; Dean, R.G. Longshore sediment transport via littoral drift rose. Ocean. Eng. 2010, 37, 228–235. [Google Scholar] [CrossRef]
  41. Carslaw, H.; Jaeger, J. Conduction of Heat in Solids; Clarendon Press: Oxford, UK, 1959; p. 510. [Google Scholar]
  42. Crank, P. The Mathematics of Diffusion, 2nd ed.; Oxford University Press: Oxford, UK, 1975; p. 421. [Google Scholar]
  43. Rosati, J.D.; Walton, T.; Bodge, K. Longshore sediment transport. In Coastal Engineering Manual; U.S. Army Corps of Engineers: Washington, DC, USA, 2002; p. 119. [Google Scholar]
  44. del Valle, R.; Medina, R.; Losada, M.A. Dependence of coefficient K on grain size. J. Waterw. Port Coast. Ocean. Eng. 1993, 119, 568–574. [Google Scholar] [CrossRef]
  45. Sana, A. Teaching fundamental concepts of coastal engineering using excel spreadsheet. Comput. Appl. Eng. Educ. 2017, 25, 304–310. [Google Scholar] [CrossRef]
  46. Weggel, R.J. Maximum breaker height. J. Waterw. Harb. Coast. Eng. Div. 1972, 98, 529–548. [Google Scholar] [CrossRef]
Figure 1. The location of the study area: Loc An sand spit, Ba Ria-Vung Tau province.
Figure 1. The location of the study area: Loc An sand spit, Ba Ria-Vung Tau province.
Water 16 02553 g001
Figure 2. Hourly water level collected at the Vung Tau Oceanographic station.
Figure 2. Hourly water level collected at the Vung Tau Oceanographic station.
Water 16 02553 g002
Figure 3. The principle of the one-line model for shoreline change, as re-drawn based on the sketch mentioned in Duy et al. (2022) [19].
Figure 3. The principle of the one-line model for shoreline change, as re-drawn based on the sketch mentioned in Duy et al. (2022) [19].
Water 16 02553 g003
Figure 4. The appearances of jetties and seawalls along the Loc An sand spit.
Figure 4. The appearances of jetties and seawalls along the Loc An sand spit.
Water 16 02553 g004
Figure 5. Shoreline evolution from 2015 to 2023 of the Loc An sand spit.
Figure 5. Shoreline evolution from 2015 to 2023 of the Loc An sand spit.
Water 16 02553 g005
Figure 6. The three stages of shoreline evolution at the Loc An sand spit.
Figure 6. The three stages of shoreline evolution at the Loc An sand spit.
Water 16 02553 g006
Figure 7. Shoreline evolution on the updrift side of a groin, as re-drawn based on the sketch mentioned in Larson et al. (1987) [35].
Figure 7. Shoreline evolution on the updrift side of a groin, as re-drawn based on the sketch mentioned in Larson et al. (1987) [35].
Water 16 02553 g007
Figure 8. Transforming the images from the UTM system to the local coordinate system.
Figure 8. Transforming the images from the UTM system to the local coordinate system.
Water 16 02553 g008
Figure 9. The construction time of the jetties. The red-dashed circle highlights the jetty located at its center.
Figure 9. The construction time of the jetties. The red-dashed circle highlights the jetty located at its center.
Water 16 02553 g009
Figure 10. Comparison between the measured shoreline and the modeled shorelines.
Figure 10. Comparison between the measured shoreline and the modeled shorelines.
Water 16 02553 g010
Figure 11. RMSE between the modeled and measured shorelines.
Figure 11. RMSE between the modeled and measured shorelines.
Water 16 02553 g011
Figure 12. Validation of the analytical model.
Figure 12. Validation of the analytical model.
Water 16 02553 g012
Figure 13. Shoreline evolution between two jetties, as re-drawn based on the sketch mentioned in Tanaka and Nadaoka (1982) [17].
Figure 13. Shoreline evolution between two jetties, as re-drawn based on the sketch mentioned in Tanaka and Nadaoka (1982) [17].
Water 16 02553 g013
Figure 14. Shoreline evolution at the jetties in case of (a) one boundary and (b) two boundaries.
Figure 14. Shoreline evolution at the jetties in case of (a) one boundary and (b) two boundaries.
Water 16 02553 g014
Figure 15. Construction of the seawall.
Figure 15. Construction of the seawall.
Water 16 02553 g015
Figure 16. Comparison between modeled shorelines and measured shorelines.
Figure 16. Comparison between modeled shorelines and measured shorelines.
Water 16 02553 g016
Figure 17. Evolution of the shoreline at the jetties (x = 960 m).
Figure 17. Evolution of the shoreline at the jetties (x = 960 m).
Water 16 02553 g017
Figure 18. Shoreline positions extracted at the jetties (y1) before and after making tidal correction.
Figure 18. Shoreline positions extracted at the jetties (y1) before and after making tidal correction.
Water 16 02553 g018
Figure 19. Distribution of sediment grain size along the coast of Ba Ria–Vung Tau province, as re-produce based on the map mentioned in Duc Anh (2021) [5].
Figure 19. Distribution of sediment grain size along the coast of Ba Ria–Vung Tau province, as re-produce based on the map mentioned in Duc Anh (2021) [5].
Water 16 02553 g019
Figure 20. Comparison of ε values between two methods.
Figure 20. Comparison of ε values between two methods.
Water 16 02553 g020
Table 1. Information about Sentinel-2 images.
Table 1. Information about Sentinel-2 images.
No.DateSensorResolution (m)Data Source
105 December 2015MSI10Sentinel 2–1LC
225 December 2015MSI10Sentinel 2–1LC
314 January 2016MSI10Sentinel 2–1LC
424 March 2016MSI10Sentinel 2–1LC
526 March 2016MSI10Sentinel 2–1LC
608 January 2017MSI10Sentinel 2–1LC
718 April 2017MSI10Sentinel 2–1LC
831 August 2017MSI10Sentinel 2–1LC
905 September 2017MSI10Sentinel 2–1LC
1010 October 2017MSI10Sentinel 2–1LC
1106 February 2018MSI10Sentinel 2–1LC
1225 October 2018MSI10Sentinel 2–1LC
1304 November 2018MSI10Sentinel 2–1LC
1413 January 2019MSI10Sentinel 2–1LC
1527 July 2019MSI10Sentinel 2–1LC
1615 October 2019MSI10Sentinel 2–1LC
1726 July 2020MSI10Sentinel 2–1LC
1904 September 2020MSI10Sentinel 2–1LC
2029 October 2021MSI10Sentinel 2–1LC
2116 February 2022MSI10Sentinel 2–1LC
2228 November 2022MSI10Sentinel 2–1LC
2313 March 2023MSI10Sentinel 2–1LC
2427 April 2023MSI10Sentinel 2–1LC
2506 July 2023MSI10Sentinel 2–1LC
2627 January 2024MSI10Sentinel 2–1LC
Table 2. Initial values of variables used for modeling the shoreline using Equation (4).
Table 2. Initial values of variables used for modeling the shoreline using Equation (4).
Diffusion coefficient, ε (m2/day)Unknown
Breaking wave angle, θ (degrees)Unknown
Modeling time, t (years)1.9
Table 3. Best fit values of ε and θ.
Table 3. Best fit values of ε and θ.
Diffusion coefficient, ε (m2/day)230
Breaking wave angle, θ (degrees)18
Modeling time, t (years)1.9
Table 4. Variation of K based on the values of D50.
Table 4. Variation of K based on the values of D50.
D50 (mm)K
0.560.35
0.600.31
0.650.28
0.700.24
0.750.21
0.800.19
0.850.17
Table 5. Variation of ε based on the values of D50.
Table 5. Variation of ε based on the values of D50.
D50 (mm)Kε (m2/day)
0.560.35347
0.600.31314
0.650.28277
0.700.24245
0.750.21216
0.800.19191
0.850.17168
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Van Duy, D.; Duc Anh, N.Q.; Viet, N.T.; Tanaka, H. Interrelationship between Wall and Beach Erosion in Loc An, Vietnam: Remote Sensing and Numerical Modeling Approaches. Water 2024, 16, 2553. https://doi.org/10.3390/w16172553

AMA Style

Van Duy D, Duc Anh NQ, Viet NT, Tanaka H. Interrelationship between Wall and Beach Erosion in Loc An, Vietnam: Remote Sensing and Numerical Modeling Approaches. Water. 2024; 16(17):2553. https://doi.org/10.3390/w16172553

Chicago/Turabian Style

Van Duy, Dinh, Nguyen Quang Duc Anh, Nguyen Trung Viet, and Hitoshi Tanaka. 2024. "Interrelationship between Wall and Beach Erosion in Loc An, Vietnam: Remote Sensing and Numerical Modeling Approaches" Water 16, no. 17: 2553. https://doi.org/10.3390/w16172553

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop