An Alternative Solution to Erosion Problems at Punta Bete-Punta Maroma, , : Conciliating Tourism and Nature Author(s): Itxaso Odériz, Edgar Mendoza, Cristina Leo, Gabriel Santoyo, Rodolfo Silva, Rubí Martínez, Ernesto Grey, and Raúl López Source: Journal of Coastal Research, 71(sp1):75-85. Published By: Coastal Education and Research Foundation DOI: http://dx.doi.org/10.2112/SI71-009.1 URL: http://www.bioone.org/doi/full/10.2112/SI71-009.1

BioOne (www.bioone.org) is a nonprofit, online aggregation of core research in the biological, ecological, and environmental sciences. BioOne provides a sustainable online platform for over 170 journals and books published by nonprofit societies, associations, museums, institutions, and presses. Your use of this PDF, the BioOne Web site, and all posted and associated content indicates your acceptance of BioOne’s Terms of Use, available at www.bioone.org/page/terms_of_use. Usage of BioOne content is strictly limited to personal, educational, and non-commercial use. Commercial inquiries or rights and permissions requests should be directed to the individual publisher as copyright holder.

BioOne sees sustainable scholarly publishing as an inherently collaborative enterprise connecting authors, nonprofit publishers, academic institutions, research libraries, and research funders in the common goal of maximizing access to critical research. Journal of Coastal Research SI 71 75-85 Coconut Creek, Florida 2014

An Alternative Solution to Erosion Problems at Punta Bete-Punta Maroma, Quintana Roo, Mexico: Conciliating Tourism and Nature

Itxaso Odériz†*, Edgar Mendoza†, Cristina Leo§, Gabriel Santoyo‡, Rodolfo Silva†, Rubí Martínez†, Ernesto Grey††, and Raúl López†† †Instituto de Ingeniería, Universidad Nacional § Huaribe-OHL Desarrollos ‡ Ecoturismo Tres Rios †† Tecnoceano Autónoma de México, Mexico City, México , México México Cancún, México www.cerf-jcr.org

ABSTRACT Odériz, I.; Mendoza, E.; Leo, C.; Santoyo, G.; Silva, R.; Martínez, R.; Grey E., and López, R., 2014. An alternative solution to erosion problems at Punta Bete-Punta Maroma , Quintana Roo, Mexico : Conciliating tourism and nature. In: Silva, R., and Strusińska-Correia, A. (eds.), Coastal Erosion and Management along Developing Coasts: Selected Cases. Journal of Coastal Research, Special Issue, No. 71, pp. 75–85. Coconut Creek (Florida), ISSN 0749-0208.

www.JCRonline.org This paper presents an erosion problem found between Punta Bete and Punta Maroma, on the , Mexico. The case is worth studying as its features are quite complex: ecologically speaking, the area contains part of the Mesoamerican barrier reef, the second largest in the world, an important area of seagrass and large mangrove fields; on the other hand regarding economic aspects, the main activity is tourism, which is of vital importance to the region and nationally. These two opposing interests need to be conciliated in the adoption of a solution to problem of erosion on the coast as both are being degraded at the present. The chronic erosion reported in recent years is due to natural causes, such as the slenderizing of the coral barriers and to anthropogenic causes such as local constructions on the shoreline that interrupt the longitudinal sediment transport. This work describes the characterization of the zone and analyses of marine climate and beach sediment. The hydrodynamics have been numerically modelled with WAPO and COCO models and the morphological response of profiles with the XBeach model. Three possible solutions are proposed and numerically tested, each of which present strengths and weakness. The three alternatives are discussed in order to help decision making as some form of intervention is urgently required.

ADDITIONAL INDEX WORDS: Beach erosion, Riviera Maya, touristic developments.

INTRODUCTION over 70% of the gross domestic product of the State of Quintana Beaches and dunes provide storm protection and damage Roo (Fay, 2013). reduction against climatologic events, for both infrastructure and

natural systems beyond the beach, (Kobayashi et al., 2009). For example, the damage found on the coast around Galveston, Texas, was found to be less where the dune was higher (Doran et al., 2009). This natural protection has also recreational, biological and economic benefits. These are some of the many reasons for preventing, or stopping, coastal erosion. This paper presents the case of the erosion of the beach between Punta Bete and Punta Maroma (henceforth referred to as PBPM). The area is located on the Mexican Caribbean, forming part of the Riviera Maya, in the state of Quintana Roo (Figure 1). The study area hosts a wide variety of activities with different characteristics as well as an ecologically invaluable environment. It is also an important tourist destination, so there are many conflicts between economic, social, ecological and political interests that increase the complexity of any possible solution to beach erosion. In the ecological part of the equation,

the area contains a section of the Mesoamerican barrier reef, the Figure 1. Location of the study area, Punta Bete-Punta Maroma. second largest in the world and an important area of seagrass;

inland there are broad mangrove fields; all of which make up a complex and delicate ecosystem. Economically, the zone has PBPM is within Mexico´s most vulnerable areas to hurricanes; great touristic potential, according to INEGI (Mexican National 47 of these storms have been registered here from 1948 to 2007 Institute for Statistics and Geography) the touristic affluence in (Silva et al., 2012). Recently, the occurrence of hurricanes e.g., Quintana Roo is the highest in Mexico and tourism contributes Gilbert in 1988 and Wilma in 2005 (www.nhc.noaa.gov/ for ______more details), caused severe damage to infrastructure, economic DOI: 10.2112/SI71-009.1 received 12 February 2014; accepted in activities and social wellbeing. revision 24 August 2014. In recent years the area has suffered severe beach erosion *Corresponding author: [email protected] which has increased the vulnerability of the infrastructure there. © Coastal Education & Research Foundation 2014

76 Odériz et al. ______

As a consequence, the tourist industry may decline and, in addition to the local economic repercussions, the attractive, positive image of Quitana Roo´s beaches could affect tourism at national level. Any erosion in PBPM is considered to have a negative effect, so adequate coastal management is required and for this the features of the area need to be characterized in order to find a solution which can conciliate the social and economic interests with the environmental elements that made this site so important a touristic destination. This work, then, focuses on the characterization of the study area as the means to acquiring sensibility in order to propose alternative solutions within the aforementioned constraints, numerically test them, choose the optimal solution and issue coastal management recommendations at PMPB.

CHARACTERIZATION OF THE STUDY AREA Figure 2. Marine elements and environment at PMPB. Littoral Cell: Land and Marine Environments The littoral cell is physically delimited by the Bete (PM) and Maroma (PM) tips, which are two natural headlands, Figure 1. Present Sedimentology Characteristics Similar to many other littoral cells in the Mexican Caribbean, The small bay formed between the two headlands shows the PMPB coastal area is largely covered by vegetation (mainly almost no sediment interchange with neighbouring beaches, mangrove) as can be seen in Figure 1. The area is located over which lets us classify the bay as a closed littoral cell. The sand an aquifer which discharges fresh water to the sea. There are few from the beach face is moved offshore in the form of cross-shore surface rivers in the area, but close to the central point between transport during storms; most of this material is usually the two headlands there are three small (Figure 2). deposited in the shallower part of the study area, but if the storm The marine zone can be divided in two areas: a) the bay, and is intense enough, the sediment is carried offshore from near the b) the reef, as can be seen in Figure 2. Within these, ten types of cliff. In both cases, once the sand has been taken away from environments can be distinguished: a shallow sandy area, beach, it is almost impossible that it get back, as the mean wave seagrass, a landward , a seaward coral reef, a seaward climate is of quite low energy. transition zone, a sand channel, patches of Gorgonacea, a deep sandy area, a pronounced cliff and the continental shelf. In the reef section the reef crest is well defined and is approximately 4 km long, varying in distance from the shore from 250 to 800 m. Up to a point, it is the existence of the reef which is responsible for the formation of PM. In the 1990s a reduction in the coral reef mass was documented (De’ath et al., 2012; Gadner et al., 2003; Hughes et al., 2003). This was due to disease, blanking and the growth of epibiont organisms (Lessios et al., 2001; Toller, 2001) As a result PM receeded. This observation coincides with the knowledge we have, that the crest height of a reef controls the amount of wave energy dissipation and that the distance between the reef-crest and the shore determines the vulnerability of the beach to morphological changes, (Ruiz de Alegria-Arzaburu et al., 2013). In this same Figure 3. Location and definition of the sand profile samples, study they found that a reef-crest degradation of 1 m can result coordinates UTM 16N. in an increase in incoming wave energy of up to 10%.

The bay section is the southern part of the study area and has shallow sandy areas with some small patches of seagrass. It is At the North, the beach is defined by mangrove, rocky and characterized by the absence of a coral reef, as coral is unable to sandy areas and at the South the beach is sandy. The sediment of grow here because of the superficial and subterranean fresh both areas is mainly of biogenic origin, formed by oolite, water discharges (Bacabes del Mar, 2013). At the southern limit calcareous shells, and coralline fragments, (Carranza et al., of the study zone, PB, an incipient patch of coral reef can be 1996). These sources generate new sediment at very low rates, found due to the absence of fresh water discharge spots. This so the system has sediment deficit and thus the resilience of the reef gives shelter to the southern part of the study area and is, beach is very low. To get an idea of the stability of the beach, probably, responsible for the existence of Punta Bete. seven sand profile samples (drying beach, swash zone and surf Geological evidence of the effect of the coral reef in the stability zone) were taken along the study area in the locations shown in and regression of PM are the beach ridges found there that can Figure 3. be seen in the left panel of Figure 2.

Journal of Coastal Research, Special Issue No. 71, 2014 An Alternative Solution To Erosion Problems At Punta Bete-Punta Maroma, Quintana Roo, Mexico 77 ______

To get an idea of the stability of the beach, seven sand profile samples (drying beach, swash zone and surf zone) were taken along the study area in the locations shown in Figure 3. Table 1 summarizes the main parameters of the sand samples, i.e. d50, fine material portion (%f), roundness (R), sphericity (Sp), shape factor (Sh) and fall velocity (w). The mean measured density of all the samples was 2.2 g/cm3. In Table 1 shows the parameters of the samples: D stands for drying beach, Sw is swash zone and Su is surf zone.

Table 1. Sediment relevant parameters.

Site Sample d50 %f R Sp Sh w (mm) (m/s) D 0.4 0.33 0.87 0.82 0.6 0.04 MK Sw 0.2 0.42 0.87 0.83 0.7 0.01 Figure 4. Recent infrastructure constructed along the PMPB beach. Su 0.2 1.73 0.86 0.82 0.7 0.01 D 0.4 0.14 0.88 0.84 0.7 0.03 TR2 Sw 0.4 0.42 0.87 0.81 0.6 0.04 Although a full topographic survey was not available, ten Su 0.4 3.59 0.78 0.86 0.6 0.03 dune profiles from the southern part of the study area (Figure 4) D 0.4 0.16 0.88 0.84 0.7 0.03 were provided by Consultores en Gestión, Política y TR1 Sw 0.6 0.42 0.86 0.82 0.6 0.05 Planificación Ambiental. These profiles were recorded between Su 0.8 0.85 0.87 0.82 0.6 0.07 April and May 2013 and are shown in Figure 5. D 0.3 0.53 0.88 0.85 0.7 0.02 K2 Sw 0.2 1.21 0.88 0.87 0.7 0.02 Su 0.1 4.52 0.87 0.87 0.7 0.01 KN Sw 0.3 0.08 0.88 0.86 0.7 0.03 KV Sw 0.2 0.39 0.86 0.82 0.6 0.02 D 0.3 0.31 0.88 0.85 0.7 0.03 PM Sw 0.3 0.24 0.88 0.85 0.7 0.03 Su 0.2 2.25 0.86 0.82 0.6 0.02

The results of the sediment analysis show that at points K2 and MK the grain size is larger in the drying beach, which indicates profile instability and that the sand is being lost from these areas. The largest sediment sizes all along the study area were found in section TR1, probably because it is close to the river discharges and some terrigenous sediments can get there. This section, according to the grain size distribution along the profile seems to be reasonably stable. In all the other sample sections it was found that the grain size distribution along the beach profile decreases away from the shoreline, which indicates that these profiles are more stable than those at K2 and MK.

Shoreline and Dune Profiles Unfortunately there is almost no historical data available for the study area; only some satellite imagery and the records of human actions in the last decade. A summary of the perturbations which took place in recent years is presented in Figure 4. The construction of groins and breakwaters and the placing of geotextile bags is visible in several sections of the Figure 5. Dune profiles from the south of the study area. coast; most of these structures were installed using touristic and recreational philosophy; disregarding the physical processes involved. In Figure 4 it can also be seen that the main The beach in the areas located through the profiles P1, P2 and consequence of the infrastructure selected is the interruption of P10 has been preserved and native vegetation up to where the the longshore sediment transport. This mechanism is highly mangrove begins can still be found. However, the rest of the undesirable as the natural capacity of the beach to distribute the profiles were modified by the construction of a golf course and sand and keep a sound state disappears. other infrastructure (P3, P5, P6, P7, P8 and P9) and only the foredune is still in a sound state. Profile P4 is in a relatively

Journal of Coastal Research, Special Issue No. 71, 2014 78 Odériz et al. ______good state of conservation, as only the top of the dune has been The methodology used, in order to define a feasible solution, flattened. comprised: the characterization of local hydrodynamics, the The evolution of the shoreline in recent years (2006 to 2010) numerical modelling of wave and current fields and the was assessed by comparing satellite photographs, Figure 6 numerical modelling of the beach response to storms (using shows detailed images for 3 zones. In A, PM, the erosion is available beach profiles) followed by the numerical testing of evident, as mentioned earlier this is mainly due to the reduction three possible alternative solutions and the selection of that with of coral reef protection. The shoreline retreated here almost 10 the best performance. m. The shoreline in the central section, B, has remained more or less in the same position, with only a small variation; some parts Marine Climate show erosion, others show accretion. Section C, PB is by far the Wave and wind data from 1948 to 2010 were taken from the most damaged as it has suffered from the naturally driven Atlas of Wave Climate by Silva et al., (2007). This data are the erosion (including the effects of Hurricane Dean in 2007) and result of a reanalysis based on the hybrid wave model WAM- also the effect of the badly planned infrastructure that blocks the HURAC, WAMDI-Group (1988) and Silva et al., (2002). longshore sediment transport. From 2006 to 2010 an average of Annual and seasonal statistical data were used to define the 29 m of beach was lost. mean and extreme regimes; the results are shown in Figure 7.

Figure 7. Mean period-significant wave height joint probability (a) and wave height and direction rose (b).

In Figure 7 it can be seen that the dominant waves comes from East and SE directions with significant wave height (Hs) lower than 1.5 m and a mean period (Tm) of 4-6 s. In winter, due to cold fronts, extreme events with directions NNE, NE and NEE, with Hs of 2-3 m and Tm of 6-8 s, occur frequently. In summer there are sporadic extreme events with higher energy: Hs of 4-13 m and Tm of 6-12 s with dominant incoming directions from E and ESE. Hurricanes have a strong negative impact on the stability of Figure 6. Shoreline evolution, 2006-2010. the beach, causing severe erosion. Three hurricanes had substantial effects: Gilbert in 1988 with a minimum central pressure of 925 mb and maximum sustained winds of 130 km/h After the analysis of the sediments, shoreline evolution and when it was close to PBPM; Wilma, 2005, passed close to the the available profiles, it can be said that the present dynamics of zone with a central pressure of 930 mb and wind velocities of the PMPB beach system are of low resilience as the beach and 120 km/h, in 2007 Dean was a weaker storm: 960 mb central dune elements have been degraded by both natural and pressure and 75 km/h wind velocity (www.noaa.gov). anthropogenic actions. In addition, as the headlands that govern Regarding wind climate, the dominant directions are E and the morphology are retreating, the shoreline in between is also NEE, with a medium velocity of 5 m/s. In winter, the main eroding, a process worsened by the infrastructure that interrupts directions are ENE and E with velocities of 10 m/s due to cold the longshore sediment transport. Therefore, taking into account fronts. The highest wind velocities occur in summer, reaching that the tourism industry here is very important for the local and 50 m/s. national economy, some form of engineering solution must be Astronomic tide time series recorded at , found. The main goals of any conservation action should be: to approximately 17 km north from the study area (see Figure 1) ensure the continuity of tourist activities, to offer security to were taken as valid for the study area; the tide regime in PBPM visitors, to recover the resilience of the system and to preserve is microtidal and semidiurnal with a mean range of 20 cm. the valuable environmental elements described earlier. The storm surge was obtained from synthetic data computed by Duran (2010) with the numerical model MATO (Posada et al., 2007), which solves the 2D Non Linear Shallow Water

Journal of Coastal Research, Special Issue No. 71, 2014 An Alternative Solution To Erosion Problems At Punta Bete-Punta Maroma, Quintana Roo, Mexico 79 ______

Equations (NLSWE). The estimated sea levels at the shoreline cell, and another with high resolution (6 × 6 m), that only as a function of wind direction and velocity are shown in Figure covered the southern part of the domain. This division was made 8. as the southern area is exposed to larger waves (very few reefs exist here), so this is the area with worst erosion. Figure 9 shows an example of the results obtained with WAPO on the low resolution grid for Hs=2 m, T=8 s, incident direction NEE and meteorological tide of 2m. It can be seen, as was mentioned before, that the southern shoreline receives higher waves, up to 2.5 m high, while the area around PM is naturally protected by the reef barrier. This scenario corresponds to severe storm conditions but the general pattern is similar in all the modelled scenarios. These results are in agreement with previous works regarding the role of the coral reef as a natural defense (i.e. Figure 8. Sea levels at PMPB as a function of wind direction and Alvarez-Filip, 2011; Gourlay and Colleter, 2005; Hearn, 1999; velocity. Lowe et al., 2005). Figure 10 shows an example of the results of wave

propagation and wave induced currents on the high resolution Modelling of Wave Propagation and Currents grid for Hs=2 m, T=7 s, incident direction NEE and null To represent the hydrodynamic conditions in the study area meteorological tide which correspond to a mild storm; again the two models were used: WAPO, Silva et al., (2005), a numerical southern part receives waves of around 2.5 m and the wave wave propagation model based on the Modified Mild Slope induced currents reveal a pattern close to PB that clearly explain Equation; and COCO, Mendoza et al., (2007), which computes why the sediments are being lost to the south. wave induced currents solving the NLSWE. A total of 80 combinations of significant wave heights (1, 2 and 4 m), mean wave periods (7, 8, 10 and 12 s), wave directions (NEE, E, ESE and SE) and meteorological tides (0, 1 and 2 m) were modelled. With all these conditions, mean regime, mild storm, severe storm and hurricane scenarios were covered.

Figure 10. Propagation of waves (a) and currents pattern (b) with high resolution mesh for Hs=2 m; Tp=7 s, dir=NEE, tide=0 m.

Response of the Beach Profile The numerical model used to analyze the beach profile behavior under storm conditions was X-Beach, Roelvink et al., (2009). The only topographic data available was the modelled Figure 9. Propagation of waves, mesh with low resolution for H=2 m; profiles presented in Figure 5. The main focus of this numerical T=8 s, dir=NEE, tide=2 m. exercise was to assess the outcome of the “do nothing” solution as well as verify the vulnerability of the actual state of the dune- The hydrodynamics were estimated in two separate grids, one beach system to mild (T1 and T2) and severe (T3, T4 and T5) of low resolution (15 × 15 m), which covered the whole littoral storms and hurricanes (T6 and T7). Table 2 summarizes the

Journal of Coastal Research, Special Issue No. 71, 2014 80 Odériz et al. ______modelled scenarios, where Mt stands for meteorological tide and (1) Only a few of the original coastal dunes are conserved and t for simulation time. in a sound state, which means that resilience is decreasing as the sand supply to the beaches during storms is limited or Table 2. Hydrodynamic parameters used in modeling Xbeach. stopped. (2) Disregarding the causes, the reduction of the reef Storm Hs (m) Tp (s) Dir Mt (m) t (h) dimensions is increasing the wave energy that arrives T1 2 8 - 0 7 onshore and thus the coastline is moving landward, looking T2 2 8 ESE/E 0.8 7 for a new stable state. T3 2 8 ENE 1.3 7 (3) The combination of density, sphericity and shape factor of T4 2 8 ESEE 1.6 7 the sediment results in easily transported sand; this, and the T5 4 10 ESE/E 0.8 3 T6 4 10 ENE 1.3 3 lack of sediment sources, means that the littoral cell has a T7 4 10 ESEE 1.6 3 sediment deficit. (4) The numerical results and the satellite imagery show the Figure 11 shows the results of X-Beach modeling, in which it process of permanent erosion in the system. can be seen that the profiles with a sound dune (profiles 1, 2, 4 (5) Several structures, aligned perpendicular to the coast, i.e., and 10) can endure mild and severe storms with less damage and groins, have been constructed along the beach interrupting having the possibility of a natural recovery. In contrast, the the longshore sediment transport and thus worsening the profiles with the most damaged dunes are rapidly eroded. An sediment deficit in some areas. important feature revealed by these results is that an important (6) The stability of the beaches is often threatened by storms part of the sand taken by the storms is deposited in the nearshore and hurricanes. area and only a small part is lost offshore; this means that it is All the above factors lead us to describe the PBPM system as available for artificial nourishment. having low vulnerability, in a state of permanent erosion with a low possibility of auto-regeneration and in urgent need of intervention in order to preserve the natural value of the area and

its social and economic interests.

SHORELINE PROTECTION ALTERNATIVES The solutions suggested assume that regenerating and protecting the beach in certain segments can be the means to recovering the protection function of the dune-beach system in the entire littoral cell and thus preserving its attractiveness and, hence, its economic viability. It is desirable that the solution be a combination of rigid structures, intended to control wave energy, and soft actions (i.e. artificial nourishment), to control sediment transport. Submerged structures were chosen as the rigid elements as they have low visual impact, produce stilted wave breaking and reduce wave transmission, which helps water renewal in the protected area Three alternatives were proposed and numerically tested. For each alternative the same 80 wave conditions used in the characterization were modelled and the profile response was evaluated with X-beach model in the ten profiles shown in Figure 5. The description of each alternative and the numerical findings are presented below. It is important to remark here that all the alternatives include beach nourishment using sand taken from the nearshore area of the system. Alternative A has four medium to large submerged structures placed in critical positions, to recover the protective function of the reef barriers. A sketch of the proposed location of the structures is shown in Figure 12. All the barriers are crested at the mean sea level. Structures A1 and A2 of alternative A are intended to extend the length of the existing reef barriers, thus reducing the wave energy that reaches the shoreline. At the same time these Figure 11. Response of the profiles to storms. structures may help to retain part of the sediment that travels offshore and cannot return the beach. In turn, structures A3 and As a result of the characterization of the study area, the main A4 are also wave control structures which would reduce the problems detected in the study area are wave energy that is reaching PM as a result of the reef slenderizing so that the new stable position of the coastline can

Journal of Coastal Research, Special Issue No. 71, 2014 An Alternative Solution To Erosion Problems At Punta Bete-Punta Maroma, Quintana Roo, Mexico 81 ______be found as soon as possible thereby reducing the dry sand loss. behind each barrier. This array of structures will not The UTM coordinates of the tips of the structures, their lengths immediately benefit the whole littoral cell; only the most and crest heights are shown in Table 3. damaged sections. Anyway, the distribution of the sand from the nourishment will strengthen the beach in the long term as long as no severe storm or hurricane event occurs. The coordinates of the tips of the structures and their main dimensions are shown in Table 4.

Figure 12. Location of the submerged structures in alternative A.

Table 3. Location and dimensions of the structures in alternative A

Figure 13. Location of the submerged structures in alternative B. Structure East (m) North (m) Length (m)

A1 498240 2286320 498916 2286980 945 Table 4. Location and dimensions of the structures in alternative B A2 500571 2288831 501279 2289586 1035 A3 504486 2292678 Structure East (m) North (m) Lenght (m) 504651 2293080 435 B1 498916 2286980 A4 504932 2293664 498191 2286440 900 505145 2294198 575 B2 498091 2287433 498054 2287267 170 B3 498208 2287767 The numerical results for mean conditions show that the 498160 2287588 185 structures have little effect on wave propagation and wave B4 499175 2289239 induced currents, but as the wave conditions become more 499238 2289316 100 energetic (mild storm, severe storm and hurricane) the structures B5 499 920 2289598 170 stop high waves reaching the coast. Some currents are found in the protected area, which help to distribute the material from the Alternative C has two groups of small submerged structures nourishment all along the littoral cell. as shown in Figure 14. In this alternative the idea of replicating In alternative B the existing structures placed to protect PM the function of the reefs has been abandoned for the sake of are eliminated, believing that the position and shape of the reducing overall costs. The effects of barriers C1 to C9 in headland can be left to reach its stable state freely, while the controlling wave energy are purely local and their efficiency is protection and artificial stabilization can be focused in the more higher for the mean regime and mild storms than for more damaged sections of the cell. The reduced cost of this alternative energetic situations. On the contrary, they are better working as is an additional reason for adopting this alternative. As can be sediment traps when intense storms occur. Figure 15 shows an seen in Figure 13, where the location of the structures is shown, example of the wave propagation and the numerically produced structure B1 is similar to structure A1 in location, dimensions wave induced current fields. and performance. On the other hand, structures B2 to B5 are The protected area offered by alternative C is smaller than in small submerged barriers intended to control wave energy but the other options, as is its impact across the whole littoral cell. also help to control sediment transport. However, this alternative has the advantage that it can be The numerical results of this alternative predict a local effect constructed in less time and with smaller equipment. of structures B2 to B5 and the possible formation semi-tombolos

Journal of Coastal Research, Special Issue No. 71, 2014 82 Odériz et al. ______

The response of the beach profiles to the small barriers of alternative C was estimated with the X-Beach model for all the scenarios presented in Table 2. Figure 16 shows the results only for those profiles in front of which a barrier was placed, i.e. profiles P2, P3 P5 P7 and P8 of Figure 5. The numerical results show that in general with barriers placed, the erosion caused by the storms is less than without them. In agreement with the wave propagation results, the submerged structures lose protection efficiency as the intensity of the storm increases and, comparing the results in Figure 16 with those presented in Figure 5, for severe storms and hurricanes, the erosion is similar with and without the barriers. A great difference in the performance of alternative C is regarding sediment transport, that is, as the structures are closer to the shore, the possibility of semi-tombolo formation is greater. This means that the longshore sediment transport will be impeded and may be blocked, so the artificial nourishment may not distribute sediment along the entire shoreline, nevertheless it will stay in the most damaged areas and will not easily be removed.

Table 5. Location and dimensions of the structures in alternative C Figure 14. Location of the submerged structures in alternative C. Structure East (m) North (m) Length (m) C1 497933 2286445 497925 2286245 200 C2 497948 2286823 497940 2286623 200 C3 498004 2287196 497970 2286999 200 C4 498096 2287562 498033 2287372 200 C5 498215 2287921 498152 2287731 200 C6 499175 2289239 499238 2289316 100 C7 499397 2289427 499496 2289439 100 C8 499920 2289598 500042 2289684 150 C9 500216 2289876 500295 2289967 120

The retreat of the shoreline is less in profiles 2 and 8 with the structures of option C; on the other hand in profiles 5, 3 and 7 there are variations depending on the hydrodynamic conditions. With respect to the dune height, in profiles P3, P5 and P8 the effect of structure is not evident due to the wall at the edge of the grid, and in profile P2 (storm 1, storm 2), P3 (storm 2) and P7 (storms 4 and 6) the height is greater with the structure. Figure 15. Wave propagation (a) and wave induced currents (b), There is a reduction in the percentage of area eroded. The with alternative C barriers for Hs=2 m; Tp=8 s, dir=NEE, tide=0 m. shoreline retreat is very evident in P2 (storms 4, 5, 6 and 7) and P8 (storms 2, 3 y 4). Modelling shows that structures work, but dune regeneration would be recommendable to obtain a resilient On the other hand, the maintenance costs needed to keep the profile and provide more protection for hydrodynamic barriers and the artificial nourishment working properly may be conditions with longer return periods. important in the event of a severe storm, perhaps even equaling Regardless the solution selected, as said before, artificial the original investment in the case of a strong hurricane. This nourishment is recommended as well as the rigid barriers alternative seems to be most environmental friendly as it will not (preferably after constructing the barriers). To select the require any construction in the marine area nor close to the coral appropriate sand bank for this nourishment several factors were reefs. considered (CEM, 2006):

Journal of Coastal Research, Special Issue No. 71, 2014 An Alternative Solution To Erosion Problems At Punta Bete-Punta Maroma, Quintana Roo, Mexico 83 ______

Grain-size distribution: The mean diameter, d50, of the fill Particle shape: The sphericity is better higher as this has sand should be larger than the original sand in order to ensure greater fall velocities, higher roundness enhances beach stability beach stability and reduce overfill. The optimum overfill ratio is and angular grains are more difficult to be suspended. 1.00 to 1.05. Colour: As important touristic activities are developed in the area, the colour is important for the visual impact of the finished project. The sand found within the littoral cell in the nearshore area has the following characteristics: d50 of 0.542 mm, 1.5% of fine material, roundness of 0.867, sphericity of 0.785 and a shape factor of 0.656; which considering the above mentioned recommendations, can be used for the artificial nourishment suggested.

DISCUSSION AND CONCLUSIONS Conciliating the natural features and ecological services of the coastal area with medium or high density touristic infrastructure is difficult, and often seen as impossible. This is not the case in places like the coast of Quintana Roo, where it is precisely the natural richness of the area which has attracted touristic development. What was not properly considered during the planning and construction of the infrastructure is that Caribbean coasts are highly vulnerable because of their sand characteristics and availability, and that the intensity and energy of hydrodynamic conditions here causes them to be very destructive. As a result, there are many tourist developments which have fallen into a vicious circle in which the infrastructure degrades the natural elements (vegetation, dunes and beaches) leaving the infrastructure itself without natural protection, then the extreme meteorological events damage the infrastructure; as a response, more poorly planned infrastructure is constructed (e.g., groins and other barriers perpendicular to the shoreline) which again degrades the natural elements. This process occurs within the frame of larger natural cycles: atmospheric (el niño, la niña), marine (sea level rise) and ecological (reef slenderizing); which put the littoral systems in critical and chronic situations. Now, political, economic and social aspects make the “do nothing” or the “remove all” scenarios unviable, so a well- planned engineering project must be carried out in order to let all those involved continue to prosper. The present paper suggests a solution combining rigid and soft measures as the means to recover the stability of all the elements in the littoral cell Punta Bete-Punta Maroma. The soft measure consists of artificial nourishment with sand taken from the nearshore area. This sand, as has been numerically demonstrated, was taken from the dry beach and placed underwater by storms so, from this point of view, the soft measure suggested can be viewed as a recovery of the original state of the beach-dune system. It is straight forward to think that this action alone may have a short life as the recurrent storms that occur in the Caribbean will take the sand back to the nearshore area. This is the main motivation for the addition of the rigid part of the proposed solution. Three alternatives for the rigid barriers were considered. Alternative A aims to reproduce the original condition of the littoral cell, in which the wave energy was controlled by the Figure 16. Propagation of waves (a) and current patterns (b) with high coral reefs, in order to stabilize the PB and PM headlands. Given resolution mesh for H=2 m; T=7 s, dir=NE-E, tide=0 m. that the position of the headlands determines the position of the shoreline of the whole littoral cell, a stable state can be reached

Journal of Coastal Research, Special Issue No. 71, 2014 84 Odériz et al. ______in mid to long term (in the short term the stability will depend Bacabes del Mar, 2013. Caracterización Biológica Marina del only on the soft measure). The construction of the large Sistema Ambiental Regional Punta Maroma – Punta Bete. structures in this alternative will prevent the sediment going Technical report. 53p. offshore but not to being transported underwater by mild storms. Carranza-Edwards, A.; Rosales-Hoz, L., and Santiago-Pérez, S., Alternative A presents greatest efficiency in controlling the 1996. A reconnaissance study of carbonates in Mexican wave energy of severe storms and hurricanes; this, the length of beach sands. Sedimentary Geology, 101(3-4), 261-268. the barriers and the depth in which they will be placed (5-7 m), CEM, 2006. Coastal Engineering Manual. U.S. Army Corps of make this the most expensive option. Engineers (USACE), Department of Army, Washington, Alternative B combines the reproduction of the reef action DC. with local protection for the most damaged beach sections. The De'ath, G.; Fabricius, K.E.; Swetman, H., and Puotinen, M., barriers suggested in this alternative will only stabilize the 2012. The 27-year decline of coral cover on the Great headland at PB, and a group of small structures will control the Barrier Reef and its causes. PNAS, 109(44), 17995-17999. waves and work as sediment transport traps. The effects of this Doran, K.S.; Plant N. G.; Stockdon H.F.; Sallenger A. H., and alternative will not affect the whole littoral cell; mainly the Serafin K. A., 2009. Hurricane Ike: Observations and southern and central parts. The small barriers, as they are placed Analysis of Coastal Change. U.S. Geological Survey Open- close to the shore, will stop the sediment in the beach being File Report, 35 p. deposited underwater by storms. Alternative B was the most Duran G., 2010. Análisis del peligro por marea de tormenta en efficient for mild and severe storms. el Golfo de México. Mexico City, Mexico: National Alternative C focuses only on protecting the beach sections University of Mexico, Master´s thesis, 276 p. (in spanish). with worst erosion problems. This alternative has 9 short Fay D., 2013. Tourists as colonizers in Quintana Roo, Mexico. barriers, placed very close to the shore (2 m depth). The main The Canadian Geographer, 32, 691-699. objective of these structures is to dissipate wave energy and to Gardner, T.A.; Cote, I.M.; Gill, J.A.; Grant, A., and Watkinson, produce protected areas where sand settlement is favoured and A.R., 2003. Long-term region wide declines in Caribbean semi-tombolos are formed. This works well with the artificial corals. Science, 301, 958-960. nourishment as the barriers stop the sand moving to the Gourlay, M.R., and Colleter, G., 2005. Wave-generated flow on nearshore, although longshore sediment transport is limited, coral reefs an analysis for two-dimensional horizontal reef- meaning that the benefits will be local and not affect all of the tops with steep faces. Coastal Engineering, 52, 353-387. littoral cell. This alternative shows better efficiency for a mean Hughes, T.P., 2003. Climate change, human impacts, and the regime and mild storms and it is possible that the barriers and resilience of coral reefs. Science, 301, 117-127. the soft measure may need maintenance and/or reconstruction Hearn, C.J., 1999. Wave-breaking hydrodynamics within coral after a strong hurricane. The costs of this maintenance are not reef systems and the effect of changing relative sea level. high and the low initial investment of this alterative is cheaper Journal of Geophysical Research, C104, 30007-30019. than A or B. Kobayashi, N.; Buck, M.; Payo, A., and Johnson, B.D., 2009. It is clear that there is no one, unique, solution and that the Berm and dune erosion during a storm. Journal of economic aspect is not the only one to be taken into account in Waterway, Port, Coastal and Ocean Engineering, 135(1), making the decision. Instead, a coastal management plan is 1-10. needed to determine which dune-beach performance and which Lessios, H.A.; Garrido, M.J., and Kessing, B.D., 2001. of the engineering measures are desirable and compatible with Demographic history of diadema antillarum, a keystone future development and tourism. herbivore on caribbean reefs. Proceedings of the Royal

Society B: Biological Sciences, 268(1483), 2347-2353. ACKNOWLEDGMENTS Lowe, R.J.; Falter, J.L.; Bandet, M.D.; Pawlak, G.; Atkinson, This publication is one of the results of the Latin American M.J.; Monismith, S.G., and Koseff, J.R., 2005. Spectral Regional Network global collaborative project ‘‘EXCEED – wave dissipation over a barrier reef. Journal of Excellence Center for Development Cooperation – Sustainable Geophysical Research, 110, 16. Water Management in Developing Countries’’ consisting of 35 Mendoza, E.; Silva R.; Baquerizo A., and Losada M.A., 2007. universities and research centres from 18 countries on 4 Quad tree based numerical model for coastal continents. The authors would like to acknowledge the support hydrodynamics simulation. In: Ocean Engineering of the German Academic Exchange Service DAAD, the Centro Research Advances. Nova Publishers, p. 243-266. de Tecnologia e Geociências da Universidade Federal de Posada, G.; Silva, R.; and Medina, R., 2008. Modelo numérico Pernambuco, the Fundação de Amparo a Ciência e Tecnologia tridimensional para transporte de un contaminante do Estado de Pernambuco-FACEPE and the Instituto de conservativo. Ingeniería Hidráulica en México. XXIII(1), Ingeniería of the Universidad Nacional Autónoma de México for 5-19. their participation in this EXCEED project. Roelvink, D.; Reniers, A.; van Dongeren, A.; van Thiel de

Vries, J.; McCall, R., and Lescinski, J., 2009. Modelling LITERATURE CITED storm impacts on beaches, dunes and barrier islands. Alvarez-Filip, L.; Cote, I.M.; Gill, J.A.; Watkinson, A.R., and Coastal Engineering, 56(11-12), 1133-1152. Dulvy, N.K., 2011. Region-wide temporal and spatial Ruiz de Alegria-Arzaburu, A.; Mariño-Tapia, I.; Enriquez, C.; variation in Caribbean reef architecture: is coral cover the Silva, R., and González-Leija, M., 2013. The role of whole story?. Global Change Biology, 17, 2470-2477.

Journal of Coastal Research, Special Issue No. 71, 2014 An Alternative Solution To Erosion Problems At Punta Bete-Punta Maroma, Quintana Roo, Mexico 85 ______

fringing coral reefs on beach morphodynamics. Silva, R.; Ruiz, G.; Mariño-Tapia, I.; Posada, G.; Mendoza, E., Geomorphology, 198, 69-83. and Escalante, E., 2012. Manmade vulnerability of the Silva, R.; Govaere G.; Salles P.; Bautista G., and Diaz G., 2002. Cancun beach system: The case of hurricane Wilma. Clean Oceanographic vulnerability to hurricanes on the Mexican - Soil, Air, Water, 40 (9), 911-919. coast. Proceedings of 6th International Workshop on Wave Toller, W.W.; Rowan, R., and Knowlton, N., 2001. Hind Casting and Forecasting (Melbourne), pp.137-151. Repopulation of zooxanthellae in the caribbean corals Silva, R.; Borthwick, A.G.L., and Eatock-Taylor, R., 2005. montastraea annularis and M. faveolata following Numerical implementation of the harmonic modified mild- experimental and disease-associated bleaching. Biological slope equation. Coastal Engineering, 55(5), 391-407. Bulletin, 201(3), 360-373. Silva, R.; Ruiz, G.; Posada, G.; Pérez, D.; Rivillas, G.; Espinal WAMDI-Group, 1988. The WAM Model – a third generation J., and Mendoza, E., 2007. Atlas de clima marítimo de la ocean wave prediction model. Journal of Physical vertiente Atlántica Mexicana. UNAM, Mexico (in spanish). Oceanography, 18, 1775-1810.

Journal of Coastal Research, Special Issue No. 71, 2014