Shorelines are constantly modified by natural processes and human activities. An accurate detection of shoreline changes in relation to coastal works provides information on coastal dynamics and is, therefore, it is very important for future planning.
3.1. Shoreline Hardening Evolution
The Santa Marta Bay beaches have been subjected to continuous intervention, mainly aimed at stopping erosional processes. For example, in the 1970s, the Puertos de Colombia Company conducted beach dredging and filling works. Subsequently, the Central French Laboratory of Hydraulics conducted the “Defending the Santa Marta Beach” study, who’s technical recommendations were not fully implemented. This gradually contributed to the continuous recession of the width of the beach and to losing its equilibrium profile, thus ultimately exposing the foundations of the pedestrian Rodrigo de Bastidas State Park in the northern sector of the study area (Sector 4). In 2000, a mitigation effort was conducted with implementing geotextile containers for straight circular trochoidal sediments, which were planned as containment breakwaters for the particulate material. However, this intervention was not successful due to construction implementation issues and poor filling material selection.
Furthermore, a beach nourishment was conducted from September 26, 2000, to September 4, 2001. This beach nourishment used properly classified quarry sand and concentrated mostly in Sector 4, with an average width of 35.0 m to the north and 17.0 m to the south. However, on more than one occasion, strong waves and breeze completely washed away all the nourishment material, transporting it along the shoreline, especially to the south. A second beach nourishment was conducted in Sectors 2 and 3 as a beach recovery measure. In these sectors, the beach no longer existed because the ocean had totally displaced the sedimentary material and constructions along the shoreline. Consequently, the beach area increased to an average width of 22.0 m. In both cases, the design elevation was controlled at 0.60 m during construction (CIOH, 2004).
Moreover, two rigid breakwaters were built in 2002 (Fig. 2, image from 2003). The first breakwater was placed in the northern area, and the second one was built in the center of the bay, approximately 30 m away from the promenade. These works were supplemented with nourishments using classified external material with similar granulometric characteristics to the sand found at this beach. However, in the same manner as with previous nourishment dredging operations, most of the material was transported south by currents, waves, and wind.
In the same year (2002), construction started for a pier in the southern portion of the Manzanares River mouth with a footbridge perpendicular to the shoreline and a “T” ending driven into piles (Fig. 2, image from 2003). The last construction, and perhaps the most important one, is the service infrastructure built during 2009–2012. This infrastructure contained two breakwaters and one seawall mainly comprising quarry rocks (limestone). These structures aimed at protecting and dissipating wave energy for docking boats (Fig. 2, image from 2011).
The total constructed length of these four structures is 2,693.55 m, whereas the maritime beach in Santa Marta Bay has a linear length of 2,322.32 m. The most outstanding structure is the international marina, with a service infrastructure of 2,201.16 m in length, followed by a 207.73 m pier, and finally the two breakwaters with a length of 147.19 and 137.74 m, respectively (Fig. 2, image from 2019).
Hard structures have been the first management strategy for coastal erosion issues along the Colombian Caribbean coast. In the first months of 2016, there were 1,484 hard structures, with the highest concentrations in tourist cities (Rangel-Buitrago et al., 2018). As a result of these works, Santa Marta Bay evidence irreversible shoreline modifications, with negative visual impacts and the loss of public access to the ocean water area, mainly in Sector 3.
3.2. Shoreline Erosion and Accretion Assessment
3.2.1. Long-Term Shoreline Evolution
The main historical shoreline changes along the Colombian Caribbean coast area reflect both erosional and accretion events, with the latter in highly localized areas (Correa et al. 2005a). The erosion of the Colombian Caribbean coast seems to have accelerated since the 70s and 80s as the size of coastal cities has also grown (Navarrete-Ramírez 2014).
Approximately 50% of this area is suffering erosional processes that cause receding coasts (Rangel-Buitrago et al. 2018). However, erosion rates in the Colombian Caribbean are very varied (Correa and Vernette 2004; Rangel-Buitrago and Posada 2005; Posada and Henao 2008; Manzolli et al. 2020; Villate et al. 2020b). In the Department of Magdalena, where this study is located, maximum rates of − 15 m.y− 1 have been observed (Ciénaga- Tasajera Sector) (Ingeominas 2008).
The results from the Santa Marta Bay shoreline changes from 1985 to 2019 reveal a complex spatial evolution, where the periods of accretion and erosion are interspersed and correlated to the different phases of rigid structure installation. The results from our long-term shoreline change assessments (1985–2019) are presented in Table 3 and Fig. 3. The positive EPR and NSM values represent shoreline progradation, and the negative values indicate retrogradation, associated to erosion processes within the sector. The bay was subdivided into four sections depending on the anthropic interventions developed over time. Here, it was possible to detect significant changes along the shoreline. In fact, both in the different periods assessed and among the four sectors defined, accretion processes prevailed.
Table 3
Computed shoreline change rates using the NSM (m) and EPR (m.y− 1) methods. Total period 1985–2019; period P1 1985–1991; period P2 1991–2005; period P3 2005–2009, subdivided into total area; Sector 1; Sector 2; Sector 3; and Sector 4. (see Fig. 1 for location of each sector).
| | Total Period | | P1 1985–1991 | | P2 1991–2009 | | P3 2009–2019 |
| | NSM | EPR | | NSM | EPR | | NSM | EPR | | NSM | EPR |
Total Area | Maximum | 89.4 | 2.6 | | 45.9 | 8.7 | | 30.3 | 1.6 | | 70.3 | 7.8 |
Minimum | -29.6 | -0.9 | | -18.4 | -3.1 | | -11.2 | -2.6 | | -46.1 | -5.9 |
Median | 10.7 | 0.4 | | 0.1 | 0.0 | | 3.4 | 0.2 | | 0.0 | 0.0 |
Average | 15.1 | 0.5 | | 5.3 | 0.9 | | 6.6 | 0.3 | | 3.9 | 0.4 |
Sector 1 | Maximum | 69.4 | 2.0 | | 45.9 | 8.7 | | 29.0 | 1.5 | | 11.4 | 1.3 |
Minimum | -14.8 | -0.4 | | -6.5 | -1.6 | | -10.1 | -2.6 | | -11.9 | -1.3 |
Median | 5.0 | 0.2 | | 6.5 | 1.1 | | 2.6 | 0.1 | | -0.6 | -0.1 |
Average | 9.4 | 0.3 | | 9.4 | 1.6 | | 6.1 | 0.3 | | -0.3 | 0.0 |
Sector 2 | Maximum | 76.3 | 2.2 | | 42.9 | 7.2 | | 25.7 | 1.4 | | 65.9 | 7.3 |
Minimum | -13.4 | -0.6 | | -9.3 | -1.6 | | 0.5 | 0.0 | | -46.1 | -5.9 |
Median | 41.5 | 1.2 | | 23.2 | 3.9 | | 13.1 | 0.7 | | 0.7 | 0.1 |
Average | 40.5 | 1.2 | | 18.8 | 3.1 | | 11.5 | 0.6 | | 8.9 | 1.0 |
Sector 3 | Maximum | 28.3 | 1.1 | | 7.1 | 1.2 | | 30.3 | 1.6 | | 0 | 0 |
Minimum | -13.3 | -0.6 | | -4.0 | -0.7 | | 0.8 | 0.0 | | 0 | 0 |
Median | 11.3 | 0.5 | | 2.3 | 0.4 | | 11.0 | 0.6 | | 0 | 0 |
Average | 11.3 | 0.5 | | 1.6 | 0.3 | | 13.7 | 0.7 | | 0 | 0 |
Sector 4 | Maximum | 89.4 | 2.6 | | 2.6 | 0.4 | | 29.2 | 1.5 | | 70.3 | 7.8 |
Minimum | -29.6 | -0.9 | | -18.4 | -3.1 | | -11.2 | -0.6 | | -11.9 | -1.3 |
Median | -5.1 | -0.2 | | -10.6 | -1.8 | | 2.2 | 0.1 | | 1.4 | 0.2 |
Average | -2.9 | -0.1 | | -9.8 | -1.6 | | 3.0 | 0.2 | | 4.0 | 0.4 |
In general, Santa Marta Bay experienced shoreline progradation from 1985 to 2019 at an average range of 15.1 m of beach width (NSM) and an average progradation rate of 0.5 m.y− 1, thereby being classified into the C4 class of coastal evolution, i.e., it is an area with stable shoreline variations (Del Río et al. 2013; Natesan et al. 2015). However, the maximum and minimum rates of 2.6 and − 0.9 m.y− 1, respectively, reveal a moderate long-term accretion trend. Further, ap-proximately 70% of the transects along the shoreline revealed a wide beach width, whereas ap-proximately 29% of the transects experienced erosion, and 1% remained stable.
Spatially, long-term changes denote higher progradation values in the central area of the Santa Marta Bay, in Sector 3. This progradation has been mainly influenced by the disposition of the coastal works in relation to coastal dynamics. Conversely, a process of shoreline transgression was identified, mainly in Sector 4 and the external S of Sector 1 (Fig. 3). In general, in Sectors 1, 2, and 3, the end rate point denotes positive mean values (0.3, 1.2, and 0.5 m.y− 1, respectively), and only Sector 4 reports average negative values, with a change rate of − 0.1 m.y− 1 (Table 3 and Fig. 4).Taking the 1985 shoreline as a foundation, we can infer that Santa Marta Bay exhibits stable prograding behavior. However, the bay also evidences different manifestations of coastal movement. For example, Sector 4 reports a net loss of beach area, even with a constant intervention and the installation of rigid structures, whereas Sector 2 has gained area. The other two sectors (1 and 3) have remained more stable, mainly due to engineering works. Sector 2 progradation can be explained through sediment contributions from the Manzanares River, which also influences a part of Sector 1. Nevertheless, in Sector 4, as there is no significant sediment contribution, and due to marine hydrodynamic processes, sediments are distributed to the extremes of the sector adjacent to the rigid structures (Figs. 3 and 4).
The behavior of different sectors regarding accretion, erosion, or stability processes is basically influenced by how anthropic structures interact with local hydrodynamics. This influence has already been discussed in different studies, such as (Di Paola et al., 2020; Manno et al., 2016; Pranzini, 2018; Rangel-Buitrago et al., 2018).
In the past century, the Colombian Caribbean coast has experienced a population increase, with intense development and urbanization, in addition to the drastic increase in the number of tourists on vacation (Rangel-Buitrago et al. 2018). The increase in the use and occupation of these sectors has required a larger number of coastal protective measures to mitigate coastal hazards, such as erosion processes, and the construction of service structures to meet the ever-increasing demands.
In Santa Marta Bay, shoreline change trends have intensified since the port construction. Subsequently, several other works have been incorporated to mitigate erosional processes. Since then, the urban beach adjacent to the port structure has faced continuous erosion processes, up to the intervention from authorities through mitigation measures, such as building hard coastal protection structures throughout the affected area. These efforts were aimed at protecting beaches, stabilizing beach profiles, attenuating wave energy, and sustaining tourism structures.
3.2.2. Short-Term Shoreline Changes
Short-term shoreline movement data reveal significant mobility in the periods analyzed. The different sectors exhibit shoreline changes, both in the sense of beach area loss and expansion. The division of analysis periods is linked to anthropic interventions.
For each period, i.e., 1985–1991, 1991–2009, and 2009–2019, we identified shoreline variability averages for each transect with respect to their distance, observing that variations were closely related to the engineering works conducted in the three anthropic intervention cycles for rigid structures.
The NSM and EPR are expressed for each period assessed in Figs. 5 and 6. Table 3 provides detailed information about the mean, median, maximum, and minimum for each sector and period.
The 1985–1991 Period. In this first period, which represents a 6-year interval, the only rigid structure was the promenade Rodrigo de Bastidas State Park (Sector 4). In 1985, a part of its structure (south of Sector 4) became exposed without the presence of a beach area. The rest of the sector reported an average beach width of 17 m. In the following years (1987 to 1991), the shoreline experienced multiple transgressions, which reduced the sandy area in this sector. The mean EPR for this Sector was − 1.6 m.y− 1 with a mean retrograde of − 9.8 m.
According to the statistical comparative analysis conducted based on the NSM values from 1985 to 1991, the shoreline exhibited stability with the mean NSM values of 5.3 m, an EPR value of 0.9 m.y− 1, and the maximum and minimum rates of 8.7 (Sector 1) and − 3.1 (Sector 4) m.y− 1. Sector 3 evidenced a slight predisposition to progradation with an average gain of 1.6 m. Even when the entire Sector 3 was evaluated (Fig. 6), we could detect small areas in which erosional processes were present. Moreover, in Sectors 1 and 2, shoreline progradation was clearly present, with an average shoreline progress rate of 1.6 and 3.1 m.y− 1, respectively. Here, Sector 2 reports a larger shoreline progradation with an average gain of 18.8 m. These two sectors are directly influenced by the mouth of the Manzanares River, which is the only sediment source in the bay. However, small shoreline transgression segments are observed in the southern border of Sector 1 and in the northern area of Sector 2 (Fig. 6).
Unlike the other three sectors, Sector 4 evidenced an intense erosion from 1985 to 1991, with a maximum shoreline transgression rate of − 3.1 m.y− 1. Here, 42% of the transects in this sector report shoreline dislocation rates exceeding − 2 m.y− 1, which imply that this sector is classified as C3 (moderate erosion) in the coastal evolution for this period.
The 1991–2009 Period (Containment Work Construction Period). The second period is the longest as it spans 18 y. In this period, an average progradation rate of 0.3 m.y− 1 was recorded, with the maximum and minimum rates of 1.6 and − 2.6 m.y− 1, respectively. In general, it is more stable than the first period, evidencing variations among the different sectors assessed. In this period, Sector 1 reduced its progradation rates from 1.6 to 0.3 m.y− 1 because of the installation of some structures in its northern area. For Sector 2, the trend of positive shoreline advancement was sustained, but at a lower average progradation rate (0.6 m.y− 1), despite not presenting any transects with erosion. Sector 3 remained quite stable at an average progradation rate of 0.7 m.y− 1. This shoreline gain rate occurred before anthropic interventions. However, as of the 2009 construction of the International Marina of Santa Marta, this bay sector has not experienced any further changes. In contrast to the first period, Sector 4 reversed the erosional trend and reported the largest number of shoreline gain transects (59%). The average progradation rate was 0.2 m.y− 1. Figures 3 and 4 depict sectors with some erosional processes, mainly in the central portion, and progradation near coastal works. The construction of seawalls provided larger shoreline stability, with a higher tendency for erosion in the central part of the sector.
The 2009–2019 Period (Post Structure Construction Period). Average shoreline progradation rates increased with the completion of the coastal structures, remaining stable at 0.4 m.y− 1 from 2009 to 2019, with the maximum and minimum rates of 7.8 and − 5.9 m.y− 1, respectively. Furthermore, when the different sections are assessed individually, an increase may be observed in the number of erosion process profiles from the previous period. Although several profiles reported negative values, accretion rates remained stable (EPR = 0 m·y− 1) for Sector 1. In this period, accretion rates are expected to increase again at a mean progradation rate of 1.0 m·y− 1, which is higher than the previous period, owing to the installation of the marina structure. This mainly favored an accumulation of sediment next to the S structure (Figs. 5 and 6). The construction of the marina in Sector 3 eliminated this stretch of beach and became significant for the 0.4 m·y− 1 accretion rates of Sector 4. The main accretion profiles were located next to the marina structure.
The main feature detected during these different coastal works construction phases was the maintenance of the beach area, reflected in the overall shoreline stability corresponding to Table 3. The results revealed a high shoreline variation and that the area had been highly vulnerable to coastal erosion before the different interventions.
In retrospect, the erosional or negative variability condition was mitigated with the construction of breakwaters and service infrastructure that generated an accretion condition that remained active until 2019, without ignoring the specific areas wherein a negative variability trend may still be observed. For example, we can mention the area next to the breakwater in Sector N (Sector 4), without affecting the width of the beach, which required 16 y to recover supported by sediment deposits between the protective and service structures erected.
The results obtained through the DSAS analysis reveal similar accretion and erosion patterns, as found in other studies that have also used this approach. A study conducted in Puducherry, India (Misra and Ramakrishnan 2020) assessed the coastal geomorphological impact of a beach restoration project to contain coastal erosion by constructing artificial reefs and beach nourishment. The results reveal the growth of a small stretch of beach, as well as the subsequent coastal stabilization, with approximately 76% of the coast in accretion. Similar to our study, this paper also observed accretion and erosion patterns, and considered this mitigation measure as beneficial for protecting this eroded shoreline.
(Zulfakar et al. 2020) also found similar results with accretion and erosion sectors in Kuala Nerus, Terengganu (Malaysia). In 2016, coastal protective structures were built after several erosion events had taken place. The results denoted that most of the shoreline eroded during preconstruction had been recovered after the structure had been constructed.
3.2.3. Coastal Zone Management
As the coastal population of the Colombian Caribbean increases and coastal erosion becomes more severe, more pressure is exerted on the Colombian government to solve this problem at all levels (Rangel-Buitrago et al. 2018).
As erosion is essentially generated by the conflict among natural processes, retreat of the shoreline, and human activities, and its solution necessarily involves a comprehensive management of the coastal zone. Attempts to stabilize the position of the shoreline through coastal works have proven to be ineffective to control the phenomenon and often generate negative long-term impacts. However, in some extreme cases, this is the fastest and most efficient way to defend public or private assets.
Although it is common for hard engineering structures, such as breakwaters and piers, to actually promote erosion (Phillips and Jones 2006; Rangel-Buitrago et al. 2018), in Santa Marta Bay, this was not the case because the assessments conducted in the last 34-year evidence a clear favoring of dynamics after the works had been completed.
Even considering the importance of structures in reducing erosion rates at a local scale, the impacts related to changes in beach dynamics with the installation of rigid structures must be assessed (Molina et al. 2019) observed erosion areas downstream ports and breakwaters and matching linings/seawalls, because of sediment interruptions caused by these structures. Consequently, rigid interventions have been conducted throughout Santa Marta Bay to contain the initial natural erosion of the northern sector and then reduce the effects from these structures. Commonly, the development of the erosion processes in adjacent areas (driftward) leads to the establishment of new structures (Cooper et al. 2009). This effect seems to be the best explanation for why a total of 1,484 hard structures were identified along the Colombian Caribbean coast in 2016, with the highest concentrations located in tourist cities (Rangel-Buitrago et al. 2018).
Another important factor that must be considered is the changes in the local wind dynamics. Large buildings, port structures, or engineering works, whether onshore or offshore, serve as a wind screen, thus changing current and wave train circulation (Posada and Henao 2008). This section of Santa Marta Bay is the most urbanized in the department and exhibits a large number of alterations resulting in a high shoreline hardening rate. Anthropic structures locally constrain sand circulation and limit contribution areas owing to the occupation of beaches and their surrounding areas, wherein sandbanks may often be found, such as the presence of boardwalks and seawalls.
Unfortunately, from a governance standpoint, coastal erosion management in Colombia fails due to a weak institutional framework coupled with diluted and compromised coastal erosion management regulations (Rangel-Buitrago et al. 2018). Finally, the literature agrees that shoreline intervention must be preceded by protocols that guarantee their effectiveness and longevity, requiring studies on their positive and negative impacts, duly measured both in magnitude and scope.