Next Article in Journal
Flow Velocity Distribution Towards Flowmeter Accuracy: CFD, UDV, and Field Tests
Next Article in Special Issue
A Regional Application of Bayesian Modeling for Coastal Erosion and Sand Nourishment Management
Previous Article in Journal
A New Tool to Estimate Inundation Depths by Spatial Interpolation (RAPIDE): Design, Application and Impact on Quantitative Assessment of Flood Damages
Previous Article in Special Issue
Hydrometeorological Conditions for the Occurrence of Aeolian Processes on the Southern Baltic Coast in Poland
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Field and Numerical Study of Resistance and Resilience on a Sea Breeze Dominated Beach in Yucatan (Mexico)

1
Laboratorio de Ingeniería y Procesos Costeros, Instituto de Ingeniería, Universidad Nacional Autónoma de México, Sisal, Yucatán 97835, Mexico
2
Laboratorio Nacional de Resiliencia Costera, Laboratorios Nacionales CONACYT, Sisal, Yucatán 97835, Mexico
3
Department of Engineering, University of Salento, 73100 Lecce, Italy
4
Department of Engineering, University of Calabria, 87036 Arcavacata di Rende, Italy
5
Department of Earth and Environmental Sciences, University of West Florida, Pensacola, FL 32514, USA
6
Departamento de Ingeniería Civil, Universidad Autónoma del Carmen, Cd. del Carmen, Campeche 24180, Mexico
*
Author to whom correspondence should be addressed.
Water 2018, 10(12), 1806; https://doi.org/10.3390/w10121806
Submission received: 3 November 2018 / Revised: 29 November 2018 / Accepted: 1 December 2018 / Published: 8 December 2018

Abstract

:
The understanding of the beach capability to resist and recover from a disturbance is of paramount importance in coastal engineering. However, few efforts have been devoted to quantifying beach resilience. The present work aims to investigate the shoreline resistance and resilience, associated to a transient disturbance, on a sandy beach. A temporary groin was deployed for 24 h on a micro-tidal sea-breeze dominated beach to induce a shoreline perturbation. Morphological changes were measured by means of beach surveys to estimate the beach perturbation and the further beach recovery after structure removal. An Empirical Orthogonal Function (EOF) analysis of the shoreline position suggests that the first EOF mode describes the spatial-temporal evolution of the shoreline owing to the groin deployment/removal. A new one-line numerical model of beach evolution is calibrated with the field surveys, reproducing both the sediment impoundment and subsequent beach recovery after the structure removal. Thus, a parametric numerical study is conducted to quantify resistance and resilience. Numerical results suggest that beach resistance associated to the presence of a structure decreases with increasing alongshore sediment transport potential, whereas resilience after structure removal is positively correlated with the alongshore diffusivity.

1. Introduction

The stability of an ecosystem depends on both resistance and resilience capability to withstand a given perturbation associated to either natural or anthropogenic disturbances [1,2,3]. The resilience concept has been widely employed in ecological [4] and social [5] sciences and disaster risk reduction [6,7]. However, studies incorporating resilience for coastal engineering applications are scarce [8,9,10] and hence further research is needed [11].
The beach resistance in the coastal vulnerability context can be associated to the amount of change produced by wave events and/or due to the presence of coastal infrastructure which alter the mean pattern of sediment transport in the coastal environment; thus, it is a measure of the beach capability to resist deviations with respect to an equilibrium morphological condition. On the contrary, the beach resilience determines the speed with which the beach morphology features (e.g., dune elevation and shoreline position) return to the pre-disturbed condition [2]. Thus, the knowledge of the beach stability (i.e., resistance and resilience) is fundamental for decision-making regarding mitigation measures against beach erosion. Beach erosion in the northern Yucatan coast is critical at many locations owing to the presence of coastal structures [12,13]. Meyer-Arendt [14] reported that construction of traditional groins, made of timber and rocks, began in the 1950s and increased significantly in the late 1960s with the construction of the port infrastructure. Furthermore, beach erosion has been exacerbated during the past decade owing to the use of impermeable groins and breakwaters. Therefore, structure removal has been considered as a mitigation measure against beach erosion in this region [13]. However, no information regarding the mechanisms controlling the recovery of the shoreline position after structure removal is available.
Temporary groins have been applied in previous studies [15,16,17] to measure alongshore sediment transport and to calibrate sediment transport formulations. More specifically, the beach morphology changes measured in such studies have been considered for the estimation of the K parameter in the CERC equation [18]. However, less efforts have been devoted to investigating beach evolution after coastal structures removal. Recent studies have focused on investigating the morphodynamic responses to seafloor artificial perturbations (e.g., excavated holes and channels) in the nearshore [19,20,21]. Moulton et al. [20] investigated the mechanisms controlling the infill of large excavated holes in the surf zone, finding that downslope gravity-driven bedload transport was important in morphological evolution for bathymetric features with large slopes.
The present work aims to investigate the shoreline resistance and resilience on a sea-breeze dominated beach by means of field observations and numerical modelling. The main findings are that, on sea breeze dominated environments, the: (i) beach resistance to the presence of a groin is negatively correlated with the alongshore sediment transport potential; and (ii) beach resilience after the structure removal is positively correlated with alongshore diffusivity.
The outline of the paper is as follows. Firstly, the study area is presented in Section 2. The experimental setup, numerical model and data analysis are described in the Materials and Methods section (Section 3). Section 4 presents the field observations during the experiment and the numerical model calibration and verification. Then, a discussion on the mechanisms controlling the shoreline resistance and resilience is presented (Section 5). Finally, concluding remarks are given in Section 6.

2. Study Area

The study area is located on a barrier island in the northern Yucatan Peninsula (Figure 1a), at the fishing village of Sisal (see Figure 1b). This coastal region is characterized by a micro-tidal range, intense sea breeze conditions, a mild continental shelf and low energy waves [22]. The field experiment was conducted between the Port and the Sisal Pier (see Figure 1c). Winds, offshore waves and mean sea level have been measured over the past years in order to characterize the main forcing mechanisms affecting the coastal region in this area (Figure 1b,c).
Wind conditions in the study area are dominated by synoptic scale patterns (i.e., Bermuda-Azores, easterly winds, cold-fronts and tropical waves) and local sea breezes [23]. The NE sea breeze winds are present throughout the year but are more frequent and intense in May. On the other hand, Central America Cold Surge events, associated with cold-front passages, are more frequent during winter months. Cold-fronts, usually originated in the Rocky Mountains [24], are characterized by sustained winds (W > 15 m s−1) from the NNW and a high-pressure system (Figure 2a). Therefore, the mean wave climate in the study area is associated to locally generated NE waves owing to sea breeze events, with significant wave height Hs < 1 m (at h = 10 m water depth) (Figure 2b). More energetic NNW swell waves (Hs > 2 m and Tp > 7 s), associated to cold-fronts, occur during winter months (see Figure 2b). Furthermore, the presence of tropical storms is ubiquitous in the area [25,26]. Tidal regime is mixed, predominantly diurnal, with a spring and neap tidal range of 0.8 m and 0.1 m, respectively [27].
Torres-Freyermuth et al. [22] conducted a field experiment to characterize the nearshore circulation during both intense local sea breeze events and synoptic Norte events. They found that during intense sea breeze events the alongshore currents significantly increase inside the surf and swash zones. Therefore, in the present case of a coast subjected to sea breeze conditions, the highly oblique winds (see Figure 2a) play a major role in driving longshore currents and consequently longshore transport [28,29,30].
The beach in the study area is composed of sand with median grain size, D50, equal to 0.3 mm [31]. Furthermore, it presents a shoreline orientation 14° south of the E-W orientation [22], which is altered near the port’s jetty and the Sisal Pier. The nearshore bathymetry is characterized by the presence of a sand bar system (Figure 3), where the outer bar is relatively alongshore uniform and inner bars present a high seasonal variability. Analysis of the shoreline variability, during the intense sea breeze season (May to September, 2015), at two transects bounding the study area suggests a small (<3 m) cross-shore variation at these locations (not shown). Therefore, the temporal groin experiment was conducted between these two transects located at the middle section between the two structures (Figure 3).

3. Materials and Methods

A description of field observations and data analysis is presented in this section. Furthermore, the numerical model employed in this work is also described.

3.1. Field Experiment

The field experiment was conducted in Spring 2015 to investigate the beach stability owing to the presence/removal of an artificial perturbation in the swash zone. Beach surveys before, during and after the structure deployment allow us to investigate shoreline resistance and resilience. We focused on a short-term sea breeze event due to: (i) the important role that sea-breeze events play in the sediment transport in the study area; (ii) the difficulties for conducting the beach surveys during more energetic wave conditions (storm conditions); and (iii) the labor-intensiveness required for obtaining high- spatial and temporal resolution morphology data for a longer period.

3.1.1. Temporary Groin

The temporary groin built for this study was based on the design proposed by [16]. The groin dimensions are consistent with the typical structures found along the northern Yucatan coast. The structure was made of 0.19 m thick wood-sections of dimensions 2.40 m by 1.20 m, with holes (0.08 m diameter) drilled and covered with a 63 μm sieve cloth in order to avoid a returning offshore flow near the structure [16]. The groin consisted of seven wood sections, installed over a frame made of iron pipes and clamps; 10 m are inside the surf/swash zone within the region of active transport and the remaining 4.4 m lie on the dry beach resulting on a total length of 14.4 m. Furthermore, sandbags were spread out along the base of the structure to avoid bed scouring that could lead to sediment bypassing underneath the structure. The sand bags consisted of polypropylene woven raffia bags (60 by 100 cm) filled at approximately 2/3 of their capacity with sand.
The deployment of each section started from the land toward the sea, allowing a 0.20 m overlap between sections. Each section consisted of three vertical pipes pounded 1.5 m into the sand bed using a hammer and one horizontal pipe, holding the three pipes with scaffold clamps (Figure 4a). The original design from [16] was improved by including two horizontal members, at the down-drift side of the structure (see Figure 4a), perpendicular to the groin and attached with clamps to an additional scaffold frame. This design provided additional resistance to alongshore forces induced by wind, waves and currents. The groin deployment took approximately three hours for a team of 12 people.

3.1.2. Data Collection

Different sensors were deployed to characterize the environmental conditions occurring during the field experiment. Wind data was measured every minute using a weather station located in a tower installed near the Sisal Port (Figure 1c). Offshore wave conditions were recorded at 10 m water depth using an RDI Acoustic Doppler Current Profiler (ADCP) located 11 km offshore (see Figure 1b for instrument location). During the experiment, three breaker lines were observed at the outer and inner bars and the inner-surf/swash zone transition. Moreover, an Acoustic Doppler Velocimeter (ADV) Nortek Vector, located onshore the inner bar (Figure 5a,b) at 0.2 m above the seabed, acquired high-frequency (16 Hz) velocity measurements during 48 h (27 May to 29 May). The instantaneous velocities were measured in the XYZ coordinate system, where velocities are defined such that u, v and w velocities correspond to the x (cross-shore), y (alongshore) and z (vertical) directions, respectively.
A temporary groin was deployed in the inner surf/swash zone on the morning (0800 local time) of 27 May and was removed 24-h later on the morning of 28 May 2015. Beach morphology was surveyed along 15 survey lines, covering the up- and down- drift sides of the temporary structure (Figure 5a). A Leica Differential Global Positioning System (DGPS) was employed using Real Time Kinematics (RTK) for conducting high-resolution topographic surveys. The equipment in RTK mode and Kinematic (phase) moving mode has a horizontal and vertical accuracy of 10 mm and 20 mm, respectively. The alongshore distance between transects varies from 2 to 6 m, with the highest resolution corresponding to those transects located close to the structure (Figure 5a). The DGPS beach surveys were conducted every two hours for 24-h to evaluate the beach resistance owing to the structure presence. Furthermore, measurements continued after the structure removal, with the same two-hour temporal resolution for 10 h and then were resumed with a lower temporal resolution (i.e., 29 May, 3 June), continuing until the beach was fully recovered. A total of 20 beach surveys were conducted. Control lines, located 50 m updrift and downdrift from the structure location (Figure 5a), were surveyed weekly to assess the natural beach variability in this area. It is worth to notice that the survey lines only cover the swash and inner surf zone and hence do not cover the entire surf zone which extends offshore. The field data is available via author’s request following the instruction in the supplementary material section.
The mean sea level was acquired every minute by a tidal gauge located inside the Sisal Port (Figure 1c). Mean sea level in the study area presents a cyclic annual variation [32], showing a minimum in July and a maximum in October. Therefore, when a short-term experiment is carried out (days to weeks), the seasonality of the mean sea level (msl) elevation should be considered. The vertical datum of the surveys corresponds to the MEX97 geoid [33] and the difference between the vertical datum and the mean water elevation during the experiment is approximately 0.2 m. Beach surveys in the present paper are referenced to the msl during the experiment period (May–June 2015). A summary of the field data collected is presented in Table 1.

3.1.3. Field Data Analysis

The shoreline position (z = 0) was extracted from each of the cross-shore DGPS survey lines. An Empirical Orthogonal Function (EOF) analysis of the shoreline position [34,35,36], with respect to the initial shoreline, was performed in order to investigate the dominant modes of variability during the experiment. The EOF analysis allows to represent the coastline data, x(y, t), as a summation of n spatial (en(y)) and temporal functions (cn(t)), where the variance explained decreases with the mode number, as follows:
x ( y , t ) = n = 1 N c n ( t ) e n ( y )
The shoreline variability associated to the presence of the temporal groin can be evaluated by analyzing the resulting spatial and temporal functions.
Wind magnitude and direction, from the meteorological tower, were averaged every 5 min. Wave statistics measured by the ADCP were computed every 60 min based on pressure and velocity measured at 2 Hz. On the other hand, the ADVs velocity information was removed when the correlation value was less than 80% with the aim of identifying potentially inaccurate measurements. The measured velocities were averaged over 512-s intervals to ensure stationary conditions and further estimate the (cross- and along- shore) currents inside the surf zone.

3.2. Numerical Model

3.2.1. Model Description

Numerical simulations have been conducted by means of a newly proposed morphodynamic model, named General Shoreline beach 1.0 (GSb), belonging to the one-line model typology [37]. This typology assumes that the beach cross-shore profile remains unchanged [38,39], thereby allowing beach change to be described uniquely in terms of the shoreline position. The peculiarity of the GSb model consists of simulating shoreline evolution based on a longshore transport formula/procedure suitable at any coastal mound: sand, gravel, cobbles, shingle and rock beaches [40,41,42,43,44]. The GSb model presents one calibration coefficient solely, KGSb, which does not depend on the grain size diameter and depends on the alongshore gradient in breaking wave height [45]. The proposed general formula/procedure considers an energy flux approach combined with an empirical/statistical relationship between the wave-induced forcing and the number of moving units. GSb model allows to determine short-term (daily base) or long-term (years base) shoreline change for arbitrary combinations and configurations of structures (groins, jetties, detached breakwaters and seawalls) and beach fills that can be represented on a modelled reach of coast. A demo version of the numerical model can be downloaded by following the instructions in the supplementary material section.

3.2.2. Data Analysis

The results from the numerical simulations have been adopted to investigate the shoreline beach resistance and resilience at the considered stretch of coast. The resistance index, RS and resilience index, RL, proposed by [2], to investigate ecological stability, were adapted for the present study. RS is then defined as,
R S ( t 0 ) = 1 2 | Δ S 0 | ( l + | Δ S 0 | )
where ∆S0 represents the cross-shore distance, at t0 = 720 h, between the perturbed shoreline and the unperturbed shoreline, in vicinity of the groin and l is the perturbation length which equals to the groin length. RS ranges between 0 and 1; minimal resistance (largest effect) to the beach perturbation corresponds to smaller values of RS. Similarly, the resilience index RL(t) is defined as,
R L ( t ) = 2 | Δ S 0 | ( | Δ S | + | Δ S 0 | ) 1
where ∆S represents the cross-shore distance, at time t, between the perturbed shoreline and the unperturbed shoreline with RL = 1 corresponding to a fully recovered shoreline (i.e., return to the pre-disturbed condition). Shoreline beach resistance and resilience are illustrated in Figure 6.

4. Results

4.1. Field Observations

4.1.1. Forcing Conditions

Wind velocity time series shows a diurnal variability associated to sea-breeze events (Figure 7a). A maximum wind speed of 16 m s−1, corresponding to the sea breeze event peak, was recorded during the experiment in the afternoon (1630 local time) of 27 May (Figure 7a). Offshore wave conditions, measured at 10 m water depth (Figure 7b), are highly correlated with local winds (Figure 7a). The wave height increased from Hs = 0.3 m measured at 10:00 to Hs = 1.0 m measured at 1800, with mean wave direction approaching from the NE (Figure 7b). The temporary structure was deployed during neap tides (light gray shade in Figure 7c) in order to decrease the influence of the tide on the effective length of the groin and hence restricting the swash zone width [16].
Sea breeze events increased the wave energy inside the surf zone (Figure 8a). The alongshore current velocity Vy becomes negligible during the land breeze (0300 h), whereas it reaches Vy > 0.3 m s−1 near the sea breeze peak (1800 h) (Figure 8b). Thus, westward currents dominated the surf zone hydrodynamics during the measured period (red solid-line in Figure 8b). On the other hand, the mean cross-shore current Vx (blue solid-line in Figure 8b) shows negligible current intensity at this location/elevation during the 48 h period.

4.1.2. Observed Morphology Evolution

Beach morphology changes are evaluated by analyzing the high spatial and temporal resolution RTK DGPS survey data (Table 1 and Figure 9). The measured beach survey before the structure deployment (0730 local time) shows the alongshore uniformity of the beach contour lines (Figure 9a). However, the deployment of a temporary impermeable structure (at x = 0 m) induces a shoreline perturbation owing to the sediment impoundment at the east side of the structure and erosion at the downdrift (west) side (a time-lapse video of the experiment is included as supplementary material). The beach survey conducted 24 h after the structure deployment (0842 of the 28 May 2015) shows significant changes for h > −0.6 m (Figure 9b). These observations are consistent with the westward alongshore current associated to sea breeze events (Figure 8b). The shoreline contour reaches a maximum advance of 6 m and maximum retreat of less than −4 m at the updrift and downdrift side, respectively. The calculated sediment volume impoundment, at the updrift side of the structure, reached 70 m3 in 24 h (2 to 6 m3/m), whereas the volume loss at the downdrift for the same period is less than 40 m3. Differences between up- and down- drift volume changes might be ascribed to the limited alongshore spatial coverage of the topographic measurements.
Alongshore uniformity for the submerged area (h < 0 m) was observed 24 h after the structure removal (0825 of the 29 May 2015). However, a clear perturbation in the subaerial beach profile was still present (Figure 9c). The disturbance smooths out during the following days, returning to the pre-disturbed condition (i.e., straight and parallel contours) 144 h after the structure removal (3 June 2015) (Figure 9d).
The EOF analysis of the shoreline position allowed us to identify the spatial and temporal shoreline evolution before and after the structure removal. The first mode of variability, which represents more than 91% of the total variance, has been associated solely to the groin influence. The spatial function of the first mode (e1(y) in Figure 10a) describes the shoreline perturbation with positive/negative values associated to accretion/erosion at the updrift/downdrift side of the temporal groin. Furthermore, the temporal evolution (c1(t) in Figure 10b) shows an increase from zero to a maximum value in the first 24 h (before structure removal). After the structure removal, c1(t) shows slower decrease reaching a value close to zero at the end of the measured period. The second and third modes (not shown) represent 7.5% of the total variance and show spatial and temporal functions not associated to the groin presence. The beach recovery time can be calculated from the decay time of the perturbation described by the first temporal function (c1(t)) which is nearly 144 h.

4.2. Numerical Modelling

The field data were employed to calibrate the GSb model. Thus, a 14.4 m long groin (10 m wet plus 4.4 m dry) was positioned at the center of the domain. The alongshore model computational domain was assumed equal to 41 m. Model grid cell resolution, DX, has been set equal to 1 m with a total number of cells, NX, equal to 41, whereas the model experiment has been simulated adopting a calculation time step, DT, equal to 0.005 h. For a direct comparison with field measurements, the recording time step of the output files has been set to 1 h. The median grain size, D50, has been set to 0.3 mm and the closure depth, h*, 0.8 m. Hourly wave conditions recorded at the ADCP located 11 km offshore the beach have been adopted as input to the GSb model (Table 2).
Table 3 shows the values of the RMSE (Root Mean Square Error) from the comparison of the observed and calculated (for different KGSb values) shoreline positions, yi,obs and yi,GSb, respectively, at 23 h after the groin deployment; the RMSE is defined as,
R M S E = i = 1 N ( y i , G S b y i , o b s ) 2 N
where N is the number of transects along the considered shoreline. A KGSb = 0.01 value determined the lowest RMSE, showing a better agreement between the calculated and the observed shoreline positions.
Numerical simulations were conducted for two different lateral boundary conditions: pinned or moving lateral boundary. Figure 11 shows the comparison of the simulated shoreline positions obtained after 23 h for the two different lateral boundaries; as expected, the shoreline shapes differ in vicinity of the lateral boundaries but overlap in vicinity of the groin. If a moving lateral boundary condition is selected, the boundary will move a specified distance over a certain time period. GSb lateral boundaries have been selected as moving boundaries.
The first beach survey was assumed as the initial shoreline position in the numerical model. In particular, at the down-drift (x = 0 m) and up-drift (x = 41 m) boundaries of the computational domain, the observed specific distances from the first beach survey, equal to −1.4 m and 1.8 m, respectively, over a period of 24 h, have been assumed. On the other hand, for the post structure removal condition the computation duration has been extended to 168 h taking as initial condition the survey at t = 24 h.
Figure 12c–g show that the calibrated numerical model satisfactorily predicts the downdrift shoreline evolution, whereas the model is not capable of fully reproducing the shoreline advance, at t = 21 h, in the updrift side (Figure 12h).
The model capability to predict the shoreline resilience after the structure removal was investigated (verification). The assumed initial shoreline position corresponds to the beach survey performed immediately before the structure removal and the numerical model is run without the structure using the daily mean conditions as measured for the following seven days. The numerical model calculated the drastic change occurring during the first 24 h after the groin removal (Figure 13a–e). Furthermore, it calculated the beach recovery occurring after approximately 7 days (Figure 13f). Therefore, within the framework of the field data gained in the investigation (sea-breeze conditions), the model can be considered as a reliable tool to conduct a numerical study on beach resistance and resilience for the adopted study area.

5. Discussion: Shoreline Stability

Once the GSb model has been calibrated and verified, it has been used to investigate the beach resistance and resilience phenomena. The main limitation in the field data is related to the limited alongshore spatial coverage (40 m) of the topographic measurements. Therefore, for the numerical setup, we considered a 10 m long groin located in the middle of a 200 m long shoreline computational domain. For the parametric study we employed different cases (Table 4) encompassing constant (low-energy) wave conditions during a 720 h period. The latter allows us to assess the resistance and resilience sensitivity to different forcing conditions on a longer time-scale basis. For all cases, the shoreline position at the end of the simulation is shown in Figure 14a. A positive correlation between the shoreline distance change and the sediment transport was observed (Figure 14b) and hence R S ( t 0 ) decreases as the sediment transport rate Q [40,41] increases (Figure 14c). Values of Δ S ( t 0 ) were estimated for wave conditions in Table 4.
The numerical model was also used to simulate the beach recovery phenomena (i.e., resilience) after the structure removal. Thus, the numerical model is initialized with the shoreline from Test 6A at t0 = 720 h (e.g., yellow line in Figure 14a), subjected to different wave forcing conditions (Table 5), without the presence of the structure. The shoreline recovery is significant during the first 24-h for all cases (Figure 15a), consistent with the field observations (Figure 9); afterward, it continues at a lower rate. The numerical results were used to compute the temporal evolution of R L using Equation (2) for each case (Figure 14b). Contrary to beach resistance, numerical results suggest that the beach resilience is not controlled by the longshore transport potential and it depends on the alongshore diffusivity G given by [40,41],
G = 2 μ ( h * + B ) H 1 / 50 , b 5 / 2 cos 2 θ b
where μ is assumed equal to 0.15 m1/2 s−1, h * is the closure depth, B is the berm elevation, H 1 / 50 , b is the value of H 1 / 50 at breaking [40,41], θ b is the wave angle breaking with respect to the mean rectilinear trend of the shoreline. The diffusivity is associated to the longshore spreading of a shoreline perturbation owing to its departure from equilibrium for the existing forcing. The numerical simulations show that, for a given value of H 1 / 50 , the beach resilience increases as the value of G increases (Figure 15b and Table 5). Therefore, alongshore diffusivity plays an important role on beach resilience in the study area.

6. Conclusions

A field and numerical study of shoreline resistance and resilience was conducted on a sea-breeze dominated sandy beach. The following conclusions were found:
(1)
Analysis of high- spatial and temporal resolution field observations showed a high growth rate of the perturbation due to a groin disturbance and a lower decaying rate after the disturbance removal.
(2)
A new shoreline evolution model was calibrated with field observation and was able to reproduce both the perturbation growth and decay observed in the field.
(3)
A parametric numerical study suggests that shoreline resistance decreases with wave conditions enhancing alongshore sediment transport; whereas the resilience increases as a function of the alongshore diffusivity.

Supplementary Materials

A demo version of the GSb numerical model, for Mac and Windows systems, has been made available for the scientific community and can be downloaded at: www.scacr.eu and field measurements are available at http://ocse.mx/en/experimento/beach-resilience-to-coastal-structures-brics upon request. A time-lapse video of the beach perturbation due to the temporary groin disturbance is included as supplementary material.

Author Contributions

Conceptualization, G.M., A.T.-F., P.A.T.; Writing-Original Draft Preparation, G.M., A.T.-F., G.R.T., A.F.; Writing-Review & Editing, G.M., A.T.-F., G.R.T., A.F., L.L.; Field work, G.M., A.T.-F., J.L., L.P.-A.; Field data Analysis, G.M., Structure design, L.P.-A., J.L., P.A.T., Numerical modeling, G.R.T., A.F., L.L.; Numerical data Analysis, G.R.T., A.F., L.L.; Funding Acquisition, G.M., A.T.-F.

Funding

This research was funded by CONACYT through the Cátedras-CONACyT (Project 1146), Investigación Científica Básica (Project 284819) and the Laboratorio Nacional de Resiliencia Costera (Project LN 293354). Additional financial support was provided by PAPIIT DGAPA UNAM (IN101218) and Instituto de Ingeniería UNAM.

Acknowledgments

We acknowledge the field support provided by students and researchers at the Laboratorio de Ingeniería y Procesos Costeros at UNAM, especially from Gonzalo Uriel Martín Ruiz, David A. Gracia, Elena Ojeda, Tonatiuh Mendoza, José Carlos Pintado-Patiño, José Alberto Zamora, Daniel Toxtega, Martín Ezquivelzeta, Miguel Ángel Valencia, Paola Espadas, Luis Ángel Gallegos, Daniel Pastrana, Jesús Aragón, Pedro Cabañas, Alejandro Astorga, Enna López, Rafael Meza, Wilmer Rey and Marcos García. Special thanks to Elena Ojeda and David Gracia for providing the time-lapse video of the experiment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pimm, S.L. The complexity and the stability of ecosystems. Nature 1984, 307, 321–326. [Google Scholar] [CrossRef]
  2. Orwin, K.H.; Wardle, D.A. New indices for quantifying the resistance and resilience of soil biota to exogenous disturbances. Soil Biol. Biochem. 2004, 36, 1907–1912. [Google Scholar] [CrossRef]
  3. Folke, C. Resilience: The emergence of a perspective for social-ecological systems analyses. Glob. Environ. Chang. 2006, 16, 253–267. [Google Scholar] [CrossRef]
  4. Holling, C.S. Resilience and stability of ecological systems. Ann. Rev. Ecol. Syst. 1973, 4, 1–23. [Google Scholar] [CrossRef]
  5. Garmezy, N.; Masten, A.S.; Tellegen, A. The study of stress and competence in children: A building block for developmental psychopathology. Child Dev. 1984, 55, 97–111. [Google Scholar] [CrossRef] [PubMed]
  6. Alexander, D.E. Resilience against earthquakes: Some practical suggestions for planners and managers. J. Seism. Earthq. Eng. 2012, 13, 109–115. [Google Scholar]
  7. Alexander, D.E. Resilience and disaster risk reduction: An etymological journey. Nat. Hazards Earth Syst. Sci. 2013, 13, 2707–2716. [Google Scholar] [CrossRef]
  8. Tomasicchio, G.R.; Sanchez Arcilla, A.; D’Alessandro, F.; Ilic, S.; James, M.; Fortes, C.J.E.M.; Sancho, F.; Schüttrumpf, H. Large-scale flume experiments on dune erosion processes. J. Hydraul. Res. 2011, 49, 20–30. [Google Scholar] [CrossRef]
  9. D’Alessandro, F.; Tomasicchio, G.R. Wave-dune interaction and beach resilience in large-scale physical model tests. Coast. Eng. 2016, 116, 15–25. [Google Scholar] [CrossRef]
  10. Smith, E.R.; D’Alessandro, F.; Tomasicchio, G.R.; Gailani, J.Z. Nearshore placement of a sand dredged mound. Coast. Eng. 2017, 126, 1–10. [Google Scholar] [CrossRef]
  11. Martinez, L.; Taramelli, A.; Silva, R. Resistance and Resilience: Facing the Multidimensional Challenges in Coastal Areas. J. Coast. Res. 2017, 77 (Suppl. 1), 1–6. [Google Scholar] [CrossRef] [Green Version]
  12. Lira-Pantoja, M.; Torres-Freyermuth, A.; Appendini, C.; Fernández, D.; Salles, P.; Mendoza, E.; López, J.; Pedrozo-Acuña, A. Chronic beach erosion induced by coastal structures in Chelem, Yucatan. Coast. Eng. Proc. 2012, 1, 125. [Google Scholar] [CrossRef]
  13. Medellín, G.; Mariño-Tapia, I.; Euán-Ávila, J. The influence of a seawall on postnourishment evolution in a sea-breeze-dominated microtidal beach. J. Coast. Res. 2015, 31, 1449–1458. [Google Scholar] [CrossRef]
  14. Meyer-Arendt, K.J. Recreational development and shoreline modification along the north coast of Yucatán, Mexico. Tour. Geogr. 2001, 3, 87–104. [Google Scholar] [CrossRef]
  15. Bodge, K.R.; Dean, R.G. Short-Term Impoundment of Longshore Sediment Transport; TR CERC-87-7; WES: Vicksburg, MS, USA, 1987. [Google Scholar]
  16. Wang, P.; Kraus, N.C. Longshore sediment transport rate measured by short-term impoundment. J. Waterw. Port Coast. Ocean Eng. 1999, 125, 118–126. [Google Scholar] [CrossRef]
  17. Martin-Grandes, I.; Hughes, J.; Simmonds, D.J.; Chadwick, A.J.; Reeve, D.E. Novel methodology for one line model calibration using impoundment on mixed beach. Coast. Dyn. 2009. [Google Scholar] [CrossRef]
  18. Rosati, J.D.; Walton, T.L.; Bodge, K. Longshore Sediment Transport, in Coastal Engineering Manual, Part III; King, D.B., Ed.; Technical Report EM 1110-2-1100; U.S. Army Corps of Engineers: Washington, DC, USA, 2002. [Google Scholar]
  19. Elgar, S.; Raubenheimer, B.; Thomson, J.; Moulton, M. Resonance in an evolving hole in the swash zone. J. Waterw. Port Coast. Ocean Eng. 2012, 138, 299–302. [Google Scholar] [CrossRef]
  20. Moulton, M.; Elgar, S.; Raubenheimer, B. A surf zone morphological diffusivity estimated from the evolution of excavated holes. Geophys. Res. Lett. 2014, 41, 4628–4636. [Google Scholar] [CrossRef]
  21. Moulton, M.; Elgar, S.; Raubenheimer, B.; Warner, J.C.; Kumar, N. Rip currents and alongshore flows in single channels dredged in the surf zone. J. Geophys. Res. Oceans 2017, 122, 3799–3816. [Google Scholar] [CrossRef]
  22. Torres-Freyermuth, A.; Puleo, J.A.; DiCosmo, N.; Allende-Arandia, M.E.; Chardón-Maldonado, P.; López-González, J.; Figueroa, B.; Ruiz de Alegría-Arzaburú, A.; Figlus, J.; Roberts, T.; et al. Nearshore hydrodynamics on a sea breeze dominated beach during intense wind events. Cont. Shelf Res. 2017, 151, 40–52. [Google Scholar] [CrossRef]
  23. Figueroa-Espinoza, B.; Salles, P.; Zavala, J. On the Wind Power Potential in the northwest of the Yucatan Peninsula in Mexico. Atmosfera 2014, 27, 77–89. [Google Scholar] [CrossRef]
  24. Pérez, E.P.; Magaña, V.; Caetano, E.; Kusunoki, S. Cold surge activity over the Gulf of Mexico in a warmer climate. Front. Earth Sci. 2014, 2, 19. [Google Scholar] [CrossRef]
  25. Appendini, C.M.; Pedrozo-Acuña, A.; Meza-Padilla, R.; Torres-Freyermuth, A.; Cerezo-Mota, R.; López-González, J.; Ruiz-Salcines, P. On the role of climate change on wind waves generated by tropical cyclones in the Gulf of Mexico. Coast. Eng. J. 2017, 59. [Google Scholar] [CrossRef]
  26. Martell-Dubois, R.; Silva-Casarin, R.; Mendoza-Baldwin, E.G.; Muñoz-Pérez, J.J.; Cerdeira-Estrada, S.; Escalante-Mancera, E.; Laiz, I. Spectral bimodality of waves produced by hurricanes in the Caribbean coastal zone off Mexico. Cienc. Mar. 2018, 44, 33–48. [Google Scholar] [CrossRef] [Green Version]
  27. Cuevas-Jiménez, A.; Euán-Ávila, J. Morphodynamics of carbonate beaches in the Yucatán Peninsula. Cienc. Mar. 2009, 35, 307–320. [Google Scholar] [CrossRef]
  28. Hubertz, J.M. Observations of local wind effects on longshore currents. Coast. Eng. 1986, 10, 275–288. [Google Scholar] [CrossRef]
  29. King, D.B. Wave and Beach Processes Modeling for Sabine Pass to Galveston Bay, Texas, Shoreline Erosion Feasibility Study; Technical Report ERDC/CHL TR-07-6; USACE Engineer Research and Development Center: Vicksburg, MS, USA, 2007; 150p. [Google Scholar]
  30. Burnette, C.; Dally, W.R. The Longshore Transport Enigma and Analysis of a 10-Year Record of Wind-Driven Nearshore Currents. J. Coast. Res. 2018, 34, 26–41. [Google Scholar] [CrossRef]
  31. Wellmann, N. Analysis of Near-Shore Sediment Samples from Sisal Beach (Mexico) Comparing Effects of Sea Breeze and el Norte Events. Ph.D. Thesis, Hochschule für Technik, Wirtschaft und Kultur Leipzig, Leipzig, Germany, 2014. [Google Scholar]
  32. Zavala-Hidalgo, J.; Morey, S.L.; O’Brien, J.O. Seasonal circulation on the western shelf of the Gulf of Mexico using a high-resolution numerical model. J. Geophys. Res. 2003, 108, 3389. [Google Scholar] [CrossRef]
  33. MEXICO97. Available online: https://www.ngs.noaa.gov/GEOID/MEXICO97/ (accessed on 20 August 2018).
  34. Medellín, G.; Medina, R.; Falqués, A.; González, M. Coastline sand waves on a low-energy beach at “El Puntal” spit, Spain. Mar. Geol. 2008, 250, 143–156. [Google Scholar] [CrossRef]
  35. Miller, J.K.; Dean, R.G. Shoreline Variability via Empirical Orthogonal Function Analysis: Part I Temporal and Spatial Characteristics. Coast. Eng. 2006, 54, 111–131. [Google Scholar] [CrossRef]
  36. Muñoz-Pérez, J.J.; Medina, R.; Tejedor, B. Evolution of longshore beach contour lines determined by the E.O.F. method. Sci. Mar. 2001, 65, 393–402. [Google Scholar] [CrossRef]
  37. Frey, A.E.; Connell, K.J.; Hanson, H.; Larson, M.; Thomas, R.C.; Munger, S.; Zundel, A. GenCade Version 1 Model Theory and User’s Guide; Technical Report ERDC/CHL TR-12-25; U.S. Army Engineer Research and Development Center: Vicksburg, MS, USA, 2012. [Google Scholar]
  38. Bruun, P. Coast Erosion and the Development of Beach Profiles; Technical Memorandum N. 44; Department of the Army, Beach Erosion Board, Corps of Engineers: Washington, DC, USA, 1954.
  39. Dean, R.G. Equilibrium beach profiles: Characteristics and applications. J. Coast. Res. 1990, 71, 53–84. [Google Scholar]
  40. Tomasicchio, G.R.; Lamberti, A.; Guiducci, F. Stone movement on a reshaped profile. In Proceedings of the 24th International Conference on Coastal Engineering, Kobe, Japan, 23–28 October 1994; Volume 2, pp. 1625–1640. [Google Scholar]
  41. Lamberti, A.; Tomasicchio, G.R. Stone mobility and longshore transport at reshaping breakwaters. Coast. Eng. 1997, 29, 263–289. [Google Scholar] [CrossRef]
  42. Tomasicchio, G.R.; Archetti, R.; D’Alessandro, F.; Sloth, P. Long-shore transport at berm breakwaters and gravel beaches. Coast. Struct. 2007, 65–76. [Google Scholar] [CrossRef]
  43. Tomasicchio, G.R.; D’Alessandro, F.; Barbaro, G.; Malara, G. General longshore transport model. Coast. Eng. 2013, 71, 28–36. [Google Scholar] [CrossRef]
  44. Tomasicchio, G.R.; D’Alessandro, F.; Barbaro, G.; Musci, E.; De Giosa, T.M. Longshore transport at shingle beaches: An independent verification of the general model. Coast. Eng. 2015, 104, 69–75. [Google Scholar] [CrossRef]
  45. Ozasa, H.; Brampton, A.H. Mathematical modeling of beaches backed by seawalls. Coast. Eng. 1980, 4, 47–64. [Google Scholar] [CrossRef]
Figure 1. Location map showing (a) the Yucatan peninsula at the SE of the Gulf of Mexico, (b) a section of the Yucatan north coast showing the barrier island and wetlands and (c) the study area location, existing monitoring systems (ADCP: Acoustic Doppler Velocimeter Profile; ADV: Acoustic Doppler Velocimeter) and coastal structures.
Figure 1. Location map showing (a) the Yucatan peninsula at the SE of the Gulf of Mexico, (b) a section of the Yucatan north coast showing the barrier island and wetlands and (c) the study area location, existing monitoring systems (ADCP: Acoustic Doppler Velocimeter Profile; ADV: Acoustic Doppler Velocimeter) and coastal structures.
Water 10 01806 g001
Figure 2. (a) Wind and (b) wave roses in the study area. Wind data are taken from the Sisal weather station MeteoSisal (www.weather underground.com 1 January 2009–17 June 2016) and the wave data were collected from an ADCP located at 10 m water depth in front of Sisal (10 December 2013–20 April 2016). WS and Hs stand for wind speed and significant wave height, respectively.
Figure 2. (a) Wind and (b) wave roses in the study area. Wind data are taken from the Sisal weather station MeteoSisal (www.weather underground.com 1 January 2009–17 June 2016) and the wave data were collected from an ADCP located at 10 m water depth in front of Sisal (10 December 2013–20 April 2016). WS and Hs stand for wind speed and significant wave height, respectively.
Water 10 01806 g002
Figure 3. Bathymetry measured two days before the experiment (25 May 2015), showing the study area (■), the ADV (•) and the updrift/downdrift control lines (black solid lines) locations. The study site is located between the jetty of the Sisal port entrance channel (left-hand side, gray line) and the Sisal pier (right-hand side gray feature). The color bar indicates the elevation with respect to the mean sea level.
Figure 3. Bathymetry measured two days before the experiment (25 May 2015), showing the study area (■), the ADV (•) and the updrift/downdrift control lines (black solid lines) locations. The study site is located between the jetty of the Sisal port entrance channel (left-hand side, gray line) and the Sisal pier (right-hand side gray feature). The color bar indicates the elevation with respect to the mean sea level.
Water 10 01806 g003
Figure 4. (a) Temporary groin design made of wood-sheets, lined with sand bags and a scaffold frame of iron pipes. (b) Picture of the shoreline perturbation 12-h after the groin deployment.
Figure 4. (a) Temporary groin design made of wood-sheets, lined with sand bags and a scaffold frame of iron pipes. (b) Picture of the shoreline perturbation 12-h after the groin deployment.
Water 10 01806 g004
Figure 5. (a) Plan view of the survey lines and the structure location, showing the breaker lines and ADV location. (b) Beach profile of the middle transect showing the ADV location with respect to the structure.
Figure 5. (a) Plan view of the survey lines and the structure location, showing the breaker lines and ADV location. (b) Beach profile of the middle transect showing the ADV location with respect to the structure.
Water 10 01806 g005
Figure 6. Definition sketch of shoreline evolution between points A and B for the cases of: (a) beach resistance owing to the groin (original shoreline: blue line; disturbed shoreline: red dashed line) and (b) beach resilience after the structure removal (perturbed shoreline: red-dashed line; shoreline position at time ti after structure removal: yellow line).
Figure 6. Definition sketch of shoreline evolution between points A and B for the cases of: (a) beach resistance owing to the groin (original shoreline: blue line; disturbed shoreline: red dashed line) and (b) beach resilience after the structure removal (perturbed shoreline: red-dashed line; shoreline position at time ti after structure removal: yellow line).
Water 10 01806 g006
Figure 7. Measured (a) wind speed and wind direction (b) offshore significant wave height and wave direction at 10 m water depth and (c) mean sea level variation η ¯ from the Sisal gauge before, during (light gray shade) and after (dark gray shade) the groin deployment.
Figure 7. Measured (a) wind speed and wind direction (b) offshore significant wave height and wave direction at 10 m water depth and (c) mean sea level variation η ¯ from the Sisal gauge before, during (light gray shade) and after (dark gray shade) the groin deployment.
Water 10 01806 g007
Figure 8. Surf zone conditions during and 24 h after the structure disturbance. (a) Significant wave height inside the surf zone and (b) cross- (blue line) and along- (red line) shore currents. Negative values in panel (b) indicate offshore/west-ward currents, while positive values indicate onshore/east-ward currents.
Figure 8. Surf zone conditions during and 24 h after the structure disturbance. (a) Significant wave height inside the surf zone and (b) cross- (blue line) and along- (red line) shore currents. Negative values in panel (b) indicate offshore/west-ward currents, while positive values indicate onshore/east-ward currents.
Water 10 01806 g008
Figure 9. Beach survey (a) right before groin deployment, (b) right before groin removal (24 h later), (c) 24 h after groin removal, (d) 144 h after groin removal. The title of each subplot indicates the weekday and local time of each survey (W = Wednesday, T = Thursday, F = Friday).
Figure 9. Beach survey (a) right before groin deployment, (b) right before groin removal (24 h later), (c) 24 h after groin removal, (d) 144 h after groin removal. The title of each subplot indicates the weekday and local time of each survey (W = Wednesday, T = Thursday, F = Friday).
Water 10 01806 g009
Figure 10. Empirical orthogonal function analysis of the shoreline data during the experiment (27 May–3 June, 2015). Only the (a) spatial function and (b) the temporal evolution of the first mode are shown.
Figure 10. Empirical orthogonal function analysis of the shoreline data during the experiment (27 May–3 June, 2015). Only the (a) spatial function and (b) the temporal evolution of the first mode are shown.
Water 10 01806 g010
Figure 11. Comparison of the simulated shoreline positions obtained after 23 h for the two different lateral boundaries.
Figure 11. Comparison of the simulated shoreline positions obtained after 23 h for the two different lateral boundaries.
Water 10 01806 g011
Figure 12. Measured and calculated shoreline position (a) 2 h, (b) 4 h, (c) 6 h, (d) 8 h, (e) 11 h, (f) 14 h, (g) 16 h and (h) 21 h after the structure deployment. Information on the brackets indicates the corresponding weekday (W = Wednesday, T = Thursday) and local time of field data and GSb simulation.
Figure 12. Measured and calculated shoreline position (a) 2 h, (b) 4 h, (c) 6 h, (d) 8 h, (e) 11 h, (f) 14 h, (g) 16 h and (h) 21 h after the structure deployment. Information on the brackets indicates the corresponding weekday (W = Wednesday, T = Thursday) and local time of field data and GSb simulation.
Water 10 01806 g012
Figure 13. Measured and calculated shoreline position (a) 2 h, (b) 6 h, (c) 8 h, (d) 10 h, (e) 26 h and (f) 1 week after the structure removal. Information on the brackets indicates the corresponding weekday (T = Thursday, F = Friday) and local time of field data and GSb simulation.
Figure 13. Measured and calculated shoreline position (a) 2 h, (b) 6 h, (c) 8 h, (d) 10 h, (e) 26 h and (f) 1 week after the structure removal. Information on the brackets indicates the corresponding weekday (T = Thursday, F = Friday) and local time of field data and GSb simulation.
Water 10 01806 g013
Figure 14. (a) Shoreline position, owing to the 10 m long groin presence, after 720 h (b) the corresponding shoreline distance increases for different wave conditions listed in Table 4 and (c) resistance index temporal evolution for different wave conditions listed in Table 4.
Figure 14. (a) Shoreline position, owing to the 10 m long groin presence, after 720 h (b) the corresponding shoreline distance increases for different wave conditions listed in Table 4 and (c) resistance index temporal evolution for different wave conditions listed in Table 4.
Water 10 01806 g014
Figure 15. Beach resilience after structure removal: (a) shoreline change with respect to time after the structure removal and (b) resilience index temporal evolution.
Figure 15. Beach resilience after structure removal: (a) shoreline change with respect to time after the structure removal and (b) resilience index temporal evolution.
Water 10 01806 g015
Table 1. Measured data, sensor employed, sampling frequency and measured period.
Table 1. Measured data, sensor employed, sampling frequency and measured period.
VariableSensorManufacturerSampling FrequencyStart DateEnd Date
MorphologyGPSLeica Geosystems (Gallen, Switzerland)2–24 h27 May 20153 June 2015
WavesADCPRD Instruments (Poway, CA, USA)30 min27 May 20153 June 2015
Surf zone currentsADV VectorNortek (Rud, Norway)16 Hz27 May 201529 May 2015
Sea levelRadar level sensorOTT HydroMet (Kempten Germany)1 min27 May 20154 June 2015
WindsWeather stationDavis Vantage Pro 2 (Hayward, CA, USA)1 min27 May 20154 June 2015
Table 2. Observed wave conditions from 27 May 2015 to 3 June 2015; M = month, D = day, H = hour, Hs = significant wave height, Tp = peak wave period, θ = wave direction.
Table 2. Observed wave conditions from 27 May 2015 to 3 June 2015; M = month, D = day, H = hour, Hs = significant wave height, Tp = peak wave period, θ = wave direction.
MDHHs (m)Tp (s)θ (°) MDHHs (m)Tp (s)θ (°) MDHHs (m)Tp (s)θ (°)
52770.543.741 529170.472.571 6130.573.336
52780.513.651 529180.492.967 6140.533.437
52790.53.826 529190.55330 6150.53.156
527100.483.638 529200.613.345 6160.44.529
527110.53.436 529210.663.716 6170.44.134
527120.473.226 529220.713.441 6180.334.336
527130.472.777 529230.73.665 6190.324.313
527140.42351 53000.63.362 61100.284.312
527150.383.819 53010.563.536 61110.273.7358
527160.394.143 53020.593.943 61120.283.862
527170.462.724 53030.713.352 61130.313.541
527180.621.432 53040.753.445 61140.313.348
527190.683.528 53050.763.583 61150.243.814
527200.743.622 53060.723.48 61160.244.122
527210.83.726 53070.693.363 61170.224.2358
527220.853.824 53080.641.413 61180.194.26
527231.054.232 53090.63.671 61190.32.5260
52800.844.519 530100.534.227 61200.24.111
52810.833.855 530110.494.214 61210.184.210
52820.93.470 530120.522.945 61220.262.364
52830.853.330 530130.453.733 61230.422.722
52840.783.832 530140.463.815 6200.463.29
52850.743.947 530150.451.653 6210.413.223
52860.664.328 530160.413.29 6220.475.2357
52870.571.432 530170.444.319 6230.413.3352
52880.484.620 530180.462.717 6240.35337
52890.473.446 530190.523.256 6250.313.82
528100.434.525 530200.58314 6260.273.432
528110.394.231 530210.613.544 6270.391.411
528120.363.825 530220.673.739 6280.232.722
528130.383.612 530230.753.832 6290.213.246
528140.373.925 53100.763.916 62100.163.328
528150.384.118 53110.713.721 62110.12366
528160.452.943 53120.783.839 62120.132.6146
528170.533.155 53130.723.532 62130.2247
528180.633.284 53140.673.728 62140.212.158
528190.633.551 53150.63.748 62150.192.5239
528200.633.420 53160.543.445 62160.152.176
528210.833.737 53170.533.9360 62170.272.180
528220.774.155 53180.562.870 62180.291.941
528230.833.454 53190.514.232 62190.292.352
52900.693.453 531100.493.632 62200.292.464
52910.783.634 531110.521.561 62210.342.462
52920.73.919 531120.434.528 62220.362.354
52930.813.856 531130.414.131 62230.362.640
52940.823.938 531140.393.151 6300.322.769
52950.774.130 531150.354.837 6310.382.655
52960.774.154 531160.313.941 6320.322.782
52970.764.129 531170.331.9352 6330.272.7353
52980.63.943 531180.432.755 6340.292.651
52990.613.738 531190.54317 6350.262.762
529100.573.645 531200.583.449 6360.27252
529110.543.717 531210.573.453 6370.269.4282
529120.523.228 531220.663.426 6380.242.953
529130.533.557 531230.653.525 6390.271.821
529140.453.630 6100.73.825 63100.243.4222
529150.393.349 6110.693.753
529160.441.525 6120.673.749
Table 3. Calculated and observed shoreline positions at 14 transects along the shoreline and relative values of RMSE (Root Mean Square Error).
Table 3. Calculated and observed shoreline positions at 14 transects along the shoreline and relative values of RMSE (Root Mean Square Error).
Longshore Distance from Groin (m)−20−14−10−8−6−4−22468101420RMSE
Cross-shore Distance from Baseline (m)Initial Shoreline11.311.611.511.811.912.012.812.713.113.513.613.514.214.9
Field data at 23 h9.99.39.19.19.29.39.518.118.518.317.817.517.116.9
GSb, KGSb = 0.00510.19.89.59.49.510.011.215.716.216.416.516.416.416.63.825
GSb, KGSb = 0.0110.19.89.49.49.410.211.415.816.316.416.416.416.416.63.817
GSb, KGSb = 0.0510.19.49.39.89.811.612.517.016.716.516.416.416.416.64.127
GSb, KGSb = 0.110.18.910.211.211.212.813.318.116.816.816.616.516.516.65.269
Table 4. GSb simulated cases to investigate beach resistance to a 10 m groin, for waves of Tp = 3.5 s, after 720 h (t0) of simulation.
Table 4. GSb simulated cases to investigate beach resistance to a 10 m groin, for waves of Tp = 3.5 s, after 720 h (t0) of simulation.
CaseH1/50 (m)θ (°)ΔS(t0) (m)RS(t0)Q × 10−5 (m3/s)
Test 1A0.232151.330.7651.52
Test 2A0.232302.280.6292.45
Test 3A0.232452.620.5852.46
Test 4A0.310152.470.6044.03
Test 5A0.310303.750.4546.47
Test 6A0.310455.110.3246.50
Table 5. GSb simulated cases to investigate beach resilience associated to the groin removal. The initial condition for all simulations corresponds to the shoreline position for Test 6A at t = 720 h. θo is the deep-water wave angle.
Table 5. GSb simulated cases to investigate beach resilience associated to the groin removal. The initial condition for all simulations corresponds to the shoreline position for Test 6A at t = 720 h. θo is the deep-water wave angle.
CaseH1/3 (m)H1/50 (m) H 1 / 50 , b   ( m ) θo (°) θ b (°) G × 10−2 (m2/s)RL(ti)
Test 1B0.1500.2320.34153.740.880.865
Test 2B0.1500.2320.33306.760.770.815
Test 3B0.1500.2320.31458.710.670.689
Test 4B0.2000.3100.43153.891.600.945
Test 5B0.2000.3100.42307.401.430.920
Test 6B0.2000.3100.394510.291.140.844

Share and Cite

MDPI and ACS Style

Medellín, G.; Torres-Freyermuth, A.; Tomasicchio, G.R.; Francone, A.; Tereszkiewicz, P.A.; Lusito, L.; Palemón-Arcos, L.; López, J. Field and Numerical Study of Resistance and Resilience on a Sea Breeze Dominated Beach in Yucatan (Mexico). Water 2018, 10, 1806. https://doi.org/10.3390/w10121806

AMA Style

Medellín G, Torres-Freyermuth A, Tomasicchio GR, Francone A, Tereszkiewicz PA, Lusito L, Palemón-Arcos L, López J. Field and Numerical Study of Resistance and Resilience on a Sea Breeze Dominated Beach in Yucatan (Mexico). Water. 2018; 10(12):1806. https://doi.org/10.3390/w10121806

Chicago/Turabian Style

Medellín, Gabriela, Alec Torres-Freyermuth, Giuseppe Roberto Tomasicchio, Antonio Francone, Peter A. Tereszkiewicz, Letizia Lusito, Leonardo Palemón-Arcos, and José López. 2018. "Field and Numerical Study of Resistance and Resilience on a Sea Breeze Dominated Beach in Yucatan (Mexico)" Water 10, no. 12: 1806. https://doi.org/10.3390/w10121806

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