1. Introduction
A beach regeneration project requires identifying the fundamental parameters that define the equilibrium profile. As some authors have previously stated [
1], we can define the equilibrium profile as the situation or state that a beach profile reaches in a constant wave situation for a sufficient time. Different authors [
2,
3,
4,
5,
6,
7,
8,
9,
10] have defined the concept of an equilibrium profile. Dean [
11] defined it as “an idealization of the conditions that occur in nature for particular sediment characteristics and stable wave conditions”; later, Dean [
12] added another definition, such as “the balance between constructive and destructive forces”.
Bruun [
13] and, later, Dean [
14] obtained the mathematical formulation that defines the balance profile in beaches, which is of the form:
However, the formulation proposed by Bruun and Dean is only valid for the submerged part of the beach profile, not including the intertidal zone [
15]. This is the reason why some other authors [
16,
17,
18,
19] have developed formulations of the equilibrium profile in two sections—the emerged and submerged profiles—which, when taken together, has been called a bi-parabolic profile (for beaches with a tidal range from 1.80 to 4.92 m). To improve the adjustment of bi-parabolic profiles, Contreras et al. [
20] proposed new design parameters in the Southwest of Spain. However, these authors did not consider the effects of reef platforms on the bi-parabolic equilibrium profile.
In coastal areas where rocky slabs or reefs play a predominant role, from the point of view of coastal morphodynamics, knowledge of reef profiles is important. Muñoz-Perez [
21] was the first researcher who approached, in a theoretical–practical way (using data from Gomez-Pina [
22]), a model that allowed the practical design of supported beaches in a relatively simple way. The studies by Muñoz-Perez et al. [
23] mainly emphasized the relationship between the Dean’s shape parameter with (A
CL) and without (A
SL) a reef:
Gomez-Pina [
1] proposed a bi-parabolic model for beach profiles in tidal seas, based on field data from the Spanish coast, that allows for the quantification of a series of parameters related to the beach profile. They also proposed a model for the treatment with a reef. The formulations and recommendations proposed are applicable to the southern coast of Spain, but the methodology can be applied to any worldwide beach. Gomez-Pina expanded the formulation of Muñoz-Perez, introducing the parameter of "relative wave height" (H/F
L)—or its equivalent term of “relative freeboard” (F
L/H)—one of the classic parameters of greater influence on energy flow transmission in submerged waterproof dykes (such as rocky slabs), as follows:
where B
L is the reef width (m), F
L is the freeboard (m), and H is the wave height (m). These parameters are represented in
Figure 1.
The equilibrium profile extends to a "seaward limit", or "closure depth", beyond which there is no significant material transport and which determines the limit of profile validity.
For reef profiles, Muñoz-Perez [
21] and Gomez-Pina [
1] defined the shape parameters of the emerged parabola. These parameters have been determined in very different coasts, which gives us a good approximation for any other beach. Nevertheless, significant errors have been detected when applying these methods in the Gulf of Cádiz, due to the diversity of tidal range, energy conditions, and sedimentological variations, together with the low number of beaches studied at the time.
Therefore, the objective of this article is determining a new methodology which can be applied to any coast worldwide, by producing formulations which minimize errors for the design parameters of a bi-parabolic profile in reef coastal areas.
2. Methodology
The study area included 71 beaches between the Guadalquivir River and the Strait of Gibraltar (
Figure 2). However, only 25 beaches had a reef flat and were, therefore, selected for analysis. A beach profile spacing of 800 m was decided and, thus, 39 profiles were chosen in total [
24]. From NW to SE, the tidal range decreases from 3.8 to 1.4 m and the wave energy, Hs,
12 (wave height exceeded 12 hours a year), varies from 3.5 to 5 m [
25]. With respect to sediment, as a general rule, D
50 increases from NW (0.25 mm) to SE (0.40 mm).
The maritime climate was obtained using the real (REMRO) and virtual (SIMAR) buoy data sets available at (
http://www.puertos.es). Moreover, cartographic, bathymetric, morphologic, and sediment data were obtained from the “Ecocartographic Study of the Province of Cádiz” [
25] and from other studies carried out previously on the SW coast of Spain [
26,
27,
28,
29].
With the tool “Arctoolbox” in the ArcMap software, the topography, bathymetry, datum, and location of the cross-section profiles were exported to a CAD file. Once the beach profile was obtained, the geometry of the profile was edited in CAD, obtaining the different points in the coordinates (x, h) that define it. The list of points (x, h) were integrated into the database. Sediment data, through ArcMap tables, were also edited and moved into the database.
For the application of our methodology, it was necessary to define the origin of the emerged (O1) and submerged (O2) parables, the emergent (Ae) and submerged (As) shape parameters, the “I12” adjustment parameter (horizontal distance between O1 and O2), the D50 or average diameter of the emerged (De) and submerged (Ds) zones, length of the reef flat (BL), closure depth, and freeboard (FL).
Figure 1 shows the design parameters of the bi-parabolic profiles supported by a reef flat in a tidal sea.
The dimensionless coefficient of the origin of the emerged parabola (CO1) is defined as the relationship between the level of the origin of the emerged parabola (ZO1) and the level of the ZHHWL (Highest High Water Level) of each of the profiles studied. Thus, if the coefficient CO1 is greater than “1”, the origin (O1) is above the HHWL. If the coefficient (CO1) has a value of “1”, the origin of the emerged parabola (O1) coincides with the HHWL. If its value is “0”, the origin of the emerging parabola would be placed at the LLWL (Lowest Low Water Level).
Similarly, the coefficient of the origin of the submerged parabola (CO2) is defined as the relationship between the level of the submerged parabola (ZO2) and the Mean Water Level (ZMWL). Thus, CO2 is “1” when the origin of the submerged parabola is exactly at the Mean level of the sea. Meanwhile, CO2 will be “0” if the origin of the submerged parabola is found on the LLWL or Datum.
The expressions obtained for are fundamental in correctly defining the nourished profile (the subscripts e and s are those referring to the emerged and submerged diameters).
The parameter I
12 indicates the distance between the poles of the emerged and submerged parables. This parameter was initially defined by Gonzalez [
18] and Bernabeu [
30] (see
in Equation (4)). The fitness of this formula has been improved and generalized by adding a correction factor “C”, depending on the zone (see Equation (5)):
The different parameters that define the profile (see
Figure 1) were adjusted to find the best fit to the real profiles. This best fit was achieved by minimizing the mean square error (MSE) between the real profile and the modelled profile.
The design parameters resulting from the present research were compared with those presented by previous authors.
Table 1 summarizes the formulations for the different design parameters, according to each author.
The MSE between the proposed model and the real profile was determined for comparison of the different formulations (see Equation (6)):
where
n is the number of profile points,
ground level at the point “i” of the real profile, and
ground level at the point “i” of the modelled profile.
A new error was designed as the quotient between the MSE calculated for each model and the MSE calculated using Dean’s formula. This result, which was already dimensionless, indicates the improvement of the proposed model with respect to that proposed by Dean. Therefore, a dimensionless error of “1” indicates that the model did not improve anything with respect to that proposed by Dean; while a dimensionless error of “0” indicates that the model fit perfectly to the real profile.
Once the parameters were validated, they can be applied to future beach regeneration projects on the SW coast of Spain.
3. Results and discussion
After analysing 25 beaches (39 profiles) on the SW coast of Spain, three different kind of profiles supported on a reef flat were identified (see
Figure 3):
The Submerged Reef Bi-parabolic Profile (SRBP), a bi-parabolic profile where the reef flat is always submerged under the LLWL or Datum (hL < 0);
The Emerged Reef Bi-parabolic Profile (ERBP), another bi-parabolic profile where the reef flat sometimes emerges over the LLWL (hL > 0); and
The Reef Mono-parabolic Profile, a mono-parabolic profile where the reef flat is sometimes emerged and sometimes submerged, but is always around the LLWL (h
L ≈ 0). If the reef flat were deeply submerged, then a SRBP would appear. Meanwhile, if the reef flat was highly emerged and the length of the reef flat was not too much (B
L < 10·F
L, according to Muñoz-Perez et al. [
23]), then an ERBP would appear.
3.1. Origin (O1) of the Emerged Parabola
After studying the 39 profiles (mono- and bi-parabolic profiles), the coefficients of the origin of the emerged parabolas (CO
1) are shown in
Figure 4 (A). It can be seen that most of the profiles had the origin of the parabola at the HHWL. Thus, the origin of the emerged parabola (O
1) is typically located (63.9%) about the HHWL, and the rest of the profiles are grouped between CO1, ranging from 0.7–0.9 (27.2%). Other positions are negligible. The location of the origin of the emerged parabola at the HHWL confirms what has already been noted by other authors (see
Table 1).
3.2. Origin (O2) of the Submerged Parabola
On the other hand, based on the adjustment of real profiles, the origins of the submerged parabolas (O
2) were mostly located (81.3%) at the MWL (see
Figure 4B). In the rest of the profiles (<20%), the origin was dispersed around the LLWL, with CO
2 coefficients ranging between 0.00 and 0.40. Thus, the identification of the MWL as the origin of the submerged parabola found in
Figure 4B is different to the former statements of other authors, such as Gomez-Pina [
1] and Gonzalez [
18]. This difference may have important repercussions in the nourishment design of the bottom part of the intertidal beach profile.
3.3. Relationship Between Sediment Size and Profile Typology
The mean D
50 of the emerged and submerged samples are presented in
Figure 5.
As can be seen in
Figure 5, the D
50 ranged from 0.22 to 0.39 mm, where submerged samples had a higher proportion of fine particles than emerged samples. Moreover, the maximum difference was observed on the ERBP profile, probably because the location of the reef flat hinders the ability of fine particles to “climb up” to the emerged beach.
3.4. Relationship between Sediment Size and Dean’s Shape Parameter "A"
3.4.1. Emerged Shape Parameter “Ae” for Reef Flat Profiles
As was explained in the methodology, the shape parameters of the emerged parabolas or A
e were adjusted to the different types of profiles by minimizing the MSE.
Figure 6 shows a comparative analysis for bi-parabolic (ERBP and SRBP) and mono-parabolic (RMP) models versus Dean’s, as a function of diameter (D
e). The correlation coefficient was very high and similar for the three typologies (0.978, 0.977, and 0.999, respectively). As can be seen, the A
e of Dean’s formula was always smaller than for any other adjustment to the real topographic data. Moreover, the A
e for SRBP (
Figure 3A) was the largest, while A
e for MRP (
Figure 3C) was the smallest (closer to Dean’s formula); the A
e for ERBP was in the middle (
Figure 3B). As is well-known, the larger the shape parameter A
e, the bigger the slope of the beach. Therefore, the beach slope of the emerged parabola is bigger when the reef flat is submerged, is smaller when the reef flat is emerged, and is smallest when the reef flat is located at about the LLWL. This may be as aerial location of the reef flat (h
L ≥ 0) makes the landward cross-shore transport of the sand by swell waves difficult. Otherwise, the deeper the reef flat (h
L< 0), the lesser the barrier effect.
3.4.2. Submerged Shape Parameter "As" for Reef Flat Rofiles
Similarly, the shape parameters of the submerged parabola (or A
s) were also adjusted to the two types of profiles by minimizing the MSE. Remember that the mono-parabolic type (RMP) does not have a submerged parabola (see
Figure 3C).
Figure 7 shows a comparative analysis for SRBP and ERBP (bi-parabolic profiles) versus Dean’s, as a function of the diameter (D
s). The correlation coefficient was very high and similar, again, for these two typologies (0.983 and 0.999, respectively). As shown, the A–D relationship is identical for both SRBP (
Figure 3A) and ERBP (
Figure 3B). Moreover, there was no difference with Dean’s formula. Thus, the methodology presented herein does not provide any significant improvement for the calculation of the submerged parabola.
For all typologies, the emerged shape parameters A
e (
Figure 6) were higher than the submerged shape parameter A
s (
Figure 7). This means that the emerged parabola has a more vertical profile.
3.5. Parameter "I12"
For the ERBP typology (see
Figure 3B), the parameter I
12 depends on the reef flat length, which will be (see
Figure 1):
With the following formula for X
c,
we obtain
We can simplify the former formulae by using the results presented above (see
Figure 4A,B). The origins of emerged and submerged parabolas are approximately the levels of the HHWL and the MWL, respectively:
Thus, we can obtain the distance between the origins of the emerged and submerged parabolas (I12), using the values of:
BL (reef flat length), which is easily measurable;
h
c (vertical distance from the reef flat to the origin of the submerged parabola O
2) approximately to the MWL, as was shown in
Figure 4B; and
A
e and A
s, which are obtained from
Figure 6 and
Figure 7 by introducing D
e and D
s, respectively (ERBP typology).
For the SRBP typology, the correction coefficient "C" with respect to the parameter I
12B (the theoretical I
12 proposed by Bernabeu [
30] was already mentioned in the methodology) has been analysed, determining a “C” value of 0.965. Thus, substituting and combining Equations (4) and (5), we obtain:
The RMP typology obviously does not have the parameter I
12 (see
Figure 3C), as it is a mono-parabolic profile (i.e., there is no submerged parabola).
3.6. Application of Design Parameters on Real Profiles
The parameters previously determined were validated using the real profiles. These real profiles were modelled to obtain results of quite acceptable adjustments. The equations and parameters used for each of the typologies are shown in
Table 2.
3.7. Comparative Analysis with other Authors
The study was completed with a comparative analysis between the different theories stated by the authors presented in
Table 1 versus the models presented herein, as shown in
Table 2.
The comparative analysis was carried out with respect to the theories proposed by Dean [
14], Muñoz-Perez [
21], Gonzalez [
18], and Gomez-Pina [
1]. The dimensionless MSE, determined as the quotient between the MSE calculated for each model (Equation (6)) and the MSE calculated using Dean’s formula, was estimated for each of the authors and typologies, as summarized in
Table 3.
Thus, in brief, as can be seen in
Table 3, the model that best fit the real profile data was always that proposed in this article. It must be highlighted that the improvement was particularly significant for ERBP; that is, when the reef flat is located above the LLWL.
4. Conclusions
The SW coast of Spain has a wide variety of profiles supported on reef flats. These reef flats interrupt the normal development of the typical profiles in tidal seas. In this paper, three types of profiles were identified, depending on the situation of the reef flat (over, below, or at the LLWL): Emerged Reef Bi-parabolic Profile (ERBP), Submerged Reef Bi-parabolic Profile (SRBP), and Reef Mono-parabolic Profile (RMP).
The emerged and the submerged parts of the profile can be modelled with Dean’s parabolic formulae (he = Ae x2/3 and hs = As x2/3, respectively), where the emerged and submerged shape parameters (Ae and As) are directly related to the sediment diameter (De and Ds). However, its adjustment to real data is not optimal. Thus, based on the analysis of 39 profiles supported on reef flats chosen from 25 beaches located in the SW of Spain, a new model has been presented herein, in order to improve the modelling of the profile.
Among the findings, it should be noted that the origin “O1” of the emerged parabola was confirmed to be located at the HHWL. On the other hand, the origin "O2" of the submerged parabola was placed at the Mean Water Level (MWL), while other authors have established it using the LLWL. Some new and easy-to-apply formulas have been derived, which allow for calculation of the horizontal distance “I12” between “O1” and “O2”.
As for the emerged parabola, the shape parameters found for SRBP, ERBP, and RMP typologies were Ae = 0.360 De0.5, Ae = 0.307 De0.5, and Ae = 0.262 De0.5, respectively, all of which were higher than that in the original Dean’s formula (Ae = 0.214 De0.5). Therefore, the beach slope of the emerged parabola is bigger when the reef flat is submerged, smaller when the reef flat is emerged, and smallest when the reef flat is located at about the LLWL. This may be explained by the aerial location of the reef flat (hL ≥ 0) making it difficult for swell waves to transport sand landward cross-shore. Otherwise, the deeper the reef flat (hL < 0), the lesser the barrier effect.
As for the submerged parabola, there were no appreciable differences between the shape parameters. Moreover, there was no difference from Dean’s formula. Thus, the methodology presented herein does not provide any notorious improvement for the calculation of the submerged parabola.
Finally, it must be highlighted that the improvement is really significant when the adjustment to real profile data of the design parameters proposed in this article is compared to the parameters proposed by other authors; especially when the reef flat is located above the LLWL.