Ground-Penetrating Radar Study of Progradational Units in Holocene Coastal Plains: Carchuna Beach (SE Spain)

: The internal geometries of progradational units on the coastal plain in southeastern Spain (Carchuna Beach, Granada) since the maximum Holocene marine transgression (ca. 7000 cal BP) have been analysed using ground-penetrating radar (GPR). Radargram analysis shows the presence of sedimentary structures consistent with the progradation of coastal plains, interpreted as laminations associated with the foreset and topset. Laterally, ﬁve middle scale progradational units (H2, H3, H4, H5 and H6) are identiﬁed. Each of these units can be grouped into two major units. The general trend towards progradation is interrupted at the boundary of the units. Minor scale units (prograding wedges), deposited at topographically lower areas, are detected in GPR proﬁles. The radargrams also show channel erosional surfaces and sigmoidal landward dipping reﬂections interpreted as washover channel deposits, with the erosion and transference of sediments between the beach and the back-ridge domain. Furthermore, the saturation of coastline sediments with saltwater rapidly attenuates the GPR signal. Thus, the GPR is also a good tool to detect saltwater intrusions in shallow aquifers.


Introduction
Coastal accretion and erosion processes are closely related to climatic conditions, to sea level oscillations [1][2][3] and to anthropogenic changes [4]. In recent decades, the increase in population in coastal areas has been considerable. This entails, synchronously, the development of urban projects, the construction of infrastructures and the exploitation of natural resources, without taking into account, in many cases, coastal environmental impacts. Knowledge of current and recent coastal dynamics is a valuable tool in predicting natural risks and hazards in these coastal areas, and therefore, for sustainable development. In this regard, examining the Holocene coastal environments provides valuable information [4].
Studies conducted in the southern region of the Iberian Peninsula have identified the maximum of the Holocene marine transgression as occurring approximately 7000 years ago [5][6][7]. Since then, sea level oscillations have been minimal, remaining near the current level [6,8,9]. In this region, in Holocene highstand shallow marine deposits, Hernández Molina et al. [10] have described infralittoral prograding wedge (IPW) stages, internally characterised by a sedimentary prograding wedge with downlap in the lower boundary and toplap in the upper boundary, as well as internal downlap surfaces. Lobo et al. [11] studied the internal structures of these offshore deposits and suggested that there are three hierarchical patterns. They described two major units, which may correspond to those defined by Hernández-Molina et al. [10], and four intermediate and 12 minor units, which the authors correlated with progradational units. In high-resolution seismic profiles of infralittoral zones, Fernández-Salas et al. [12] observed a more complex distribution of IPWs due to the lateral deposition of individual prograding wedges related to longshore transport. These bodies result from the accumulation of small sedimentary units that would reflect minor sea-level fluctuations.
Conversely, in adjacent coastal plains, prograding bodies have also been described for this time interval. Thus, in zones with an adequate sediment input, two phases of coastal progradation occurred, separated by an erosional gap that dates back to approximately 2400 years ago [13,14]. These bodies are correlated with the two major units described by Lobo et al. [11] in offshore zones. They described three progradational units in the first unit (H1, H2 and H3, sensu Goy et al. [9] and three others in the second unit (H4, H5 and H6 sensu Goy et al. [9]). The bases of these units (H2, H3, H4, H5 and H6, respectively) have been dated to 5400, 4200, 2700, 1900 and 500 years before present (BP), respectively, in studies conducted by Zazo et al. [6] on the Mediterranean coast.
Carchuna Beach, in the province of Granada (SE Spain, Figure 1), has been one of the outcrops on which studies of these progradational units of the coastal plain have focused [9,[12][13][14][15][16]. This beach is 4 km long and is bordered to the west by Cape Sacratif and to the east by the tip of Cerro Gordo ( Figure 1).
Geosciences 2020, 10, x FOR PEER REVIEW 2 of 11 suggested that there are three hierarchical patterns. They described two major units, which may correspond to those defined by Hernández-Molina et al. [10], and four intermediate and 12 minor units, which the authors correlated with progradational units. In high-resolution seismic profiles of infralittoral zones, Fernández-Salas et al. [12] observed a more complex distribution of IPWs due to the lateral deposition of individual prograding wedges related to longshore transport. These bodies result from the accumulation of small sedimentary units that would reflect minor sea-level fluctuations.
Conversely, in adjacent coastal plains, prograding bodies have also been described for this time interval. Thus, in zones with an adequate sediment input, two phases of coastal progradation occurred, separated by an erosional gap that dates back to approximately 2400 years ago [13,14]. These bodies are correlated with the two major units described by Lobo et al. [11] in offshore zones. They described three progradational units in the first unit (H1, H2 and H3, sensu Goy et al. [9] and three others in the second unit (H4, H5 and H6 sensu Goy et al. [9]). The bases of these units (H2, H3, H4, H5 and H6, respectively) have been dated to 5400, 4200, 2700, 1900 and 500 years before present (BP), respectively, in studies conducted by Zazo et al. [6] on the Mediterranean coast.
Carchuna Beach, in the province of Granada (SE Spain, Figure 1), has been one of the outcrops on which studies of these progradational units of the coastal plain have focused [9,[12][13][14][15][16]. This beach is 4 km long and is bordered to the west by Cape Sacratif and to the east by the tip of Cerro Gordo ( Figure 1).  Human activity does not allow us to observe the prograding coastal plain units (b). Aerial photograph of the studied area in 1946. The Holocene prograding coastal plain units (H2, H3, H4, H5 and H6) described [12] have been superimposed (c).
Geologically, two sets can be differentiated in the sector: an underlying metamorphic basement and an overlying sedimentary cover, which is the current coastal plain. The metamorphic basement (1 in Figure 1) forms part of the Alpujárride Complex of the Internal Zones of the Betic Cordillera and consists of Palaeozoic schists, Permian-Triassic phyllites and Triassic marbles [17]. In the detritic sediments of the overlying sedimentary cover, three units are differentiated: on the one hand, the mappable conglomerate deposits in the north zone and at the bottom of the gullies (2 in Figure 1) are interpreted as alluvial deposits of Pleistocene-Middle Holocene age [17]; on the other hand, to the south, in the coastal zone, sand and gravel facies of Holocene age are associated with progradational units of the coastal plain [12][13][14] (3 in Figure 1). Based on photointerpretation and 14 C dating [13,14] Goy et al. [9] suggest that these structures were generated after the maximum Holocene transgression and that up to five of the six progradational units described in the SE of the Iberian Peninsula are separated from each other by small erosion episodes (H2, H3, H4, H5 and H6 in Figure 1).
One of the characteristics of this beach is that it is not associated with the mouths of rivers with significant sedimentation rates. In fact, the sedimentary contribution is related to regional coastal drift towards the west-southwest from the Adra River [14]. Another noteworthy characteristic of this beach is the barrier that Cape Sacratif forms, which determines the eastward formation of the aforementioned progradational units [13,14].
Despite previous studies on this beach, the internal morphology of these bodies has not been described due to the poor quality of the outcrops, which are partly covered by more recent alluvial deposits (4 in Figure 1a) and by the current greenhouses ( Figure 1b). Thus, this study aims to characterise the internal structure of these coastal deposits and to analyse in detail the morphologies of the sedimentary bodies in the boundaries of these middle scale units (Figure 1c). This study will provide new information on events during the Holocene age defined in the literature in SE Spain.
Our study uses ground-penetrating radar (GPR). This geophysical prospecting method has already been successfully used to characterise the internal structures of sedimentary bodies, particularly the internal structure of coastal units [18][19][20][21][22][23][24][25]. Therefore, we propose that GPR is an appropriate and useful tool in this case.

Material and Methods
GPR was used as a geophysical survey method. It is a non-destructive method based on the study of high-frequency electromagnetic wave propagation in the subsurface. Penetration depth and resolution depend on the electromagnetic properties of the soil materials and on the frequency of the electromagnetic signal emitted by the antenna used: high-frequency antennas provide high resolution and low penetration depth, whereas low-frequency antennas provide low resolution and high penetration depth [26,27].
The ground-penetrating radar used in this study was the RAMAC/GPR system, model Pro-Ex, manufactured by MALA GEOSCIENCE (http://www.malags.com/home). GPR profiles were acquired using 100-MHz and 250-MHz shielded antennas (Figures 2-4) in these progradational units, next to the current beach ( Figure 1c).
The GPR data acquisition software Ground Vision (http://www.ingeodav.fcen.uba.ar/Curso/ prospeccion%20geofisica/Martingarcia/Georradar/GroundVision%20Manual.pdf) was used in this study because it allows us to set, for each recording, the measurement parameters, including the antenna frequency, time window, sampling frequency or the zero time of the profile. Thus, sampling frequencies of 1089 MHz and 2540 MHz were used for the 100-MHz and 250-MHz antennas, respectively, with a 426.9-189.8 ŋs time window (512 samples per profile) and with a 0.03 m profile interval.
In total, 67 profiles were recorded, totalling 4500 m. In this study, seven profiles were selected to provide an overview of the distribution of sedimentary bodies in this coastal plain. The position of each unit is shown in Figure 1. The radargrams, and their respective locations, are shown in  prograding wedges characterised by downlap in the lower boundary and toplap in the upper boundary resulting from small sea-level falls (approximately 0.3-0.5 m) and mild storm activity (some erosion tracks are identified). Locally, convex reflections, reaching inclinations from 0° to 6°, have been interpreted as berm ridge facies (bf in Figure 2) sensu Montesa et al. [24]. On the top of this facies, the landward terminations grade downlap to onlap. Toward the sea, the bf facies are limited by an erosive surface ( Figure 2).  The signal recorded in the field was processed using the software Reflexw [Sandmeier, K. J. REFLEXW Version 7.0]. First, time zero (first arrival of a radar wave) was adjusted using the static correction filter "move startime" (9 ŋs for the 100-MHz antenna and 15 ŋs for the 250-MHz antenna). In profiles perpendicular to the coastline, with an elevation difference greater than one metre, the topography was corrected using the options "3D-topography" and "topography migration". In profiles parallel to the coastline, with minimal grade difference, these corrections were unnecessary. The filter "subtract DC shift" was used to remove residual voltage, and the filter "subtract mean" (dewow) was used to remove the direct wave. The option "automatic gain control" was selected to amplify the signal. Lastly, the filter "background removal" was used to remove background noise.
Geosciences 2020, 10, x FOR PEER REVIEW 7 of 11 been interpreted as washover channel deposits, with the erosion and transference of sediments between the beach and the back-barrier domain [24,30,33].  The velocity of electromagnetic waves in these sediments was calculated using two methods. The first was an experimental method because the depth of the water table and the time required for the wave to arrive and be reflected (in the radargram) was measured in some of the wells near the coastline (Figure 1b). This velocity was confirmed in the data processing stage, by the second method, using the GPR processing software Reflexw, (specifically the option velocity adaptation, by fitting Geosciences 2020, 10, 277 6 of 11 diffraction hyperbola). Both methods indicate a value of 0.11 m/ŋs, which is similar to those assessed in beach-ridge deposits in other sectors [21,28,29].

Conclusions
The GPR profiles recorded with the 250-MHz shielded antenna provide a good balance between penetration depth and resolution for these sandy coastal plain sediments. This method differentiates the southeastward progradation of sedimentary bodies with planar low-angle cross-bedding separated by both foreset and topset reactivation surfaces. These sedimentary structures are typical of prograding beach face deposits in a regressive context. The profile perpendicular to the coastline shows parallel laminations dipping towards the sea associated with the current foreshore. Conversely, the water table level is identified at 3.5 m depth (varying slightly, depending on the sector), which matches the reflection present in the profiles and data measured directly in wells located in the area. Furthermore, the presence of saltwater in sediments near the coastline rapidly

Results and Discussion
The profiles recorded using the 250-MHz antenna more clearly show the internal structures of the sedimentary bodies than those recorded using the 100-MHz antenna. The penetration of radar electromagnetic waves can be quite variable, depending mostly on the electrical conductivity of the subsoil. High conductivity increases signal attenuation. Therefore, the freshwater-saltwater interface is inferred from the absence of reflected radar energy from the high-conductivity zones. This behavior allows for using the GPR for determining moisture content and its variation, flow processes and water quality [30]. Accordingly, profile 25, parallel to the coastline (Figure 2), has a reflector (from 3.2 to 3.6 m in depth) under which the signal weakens and is then lost, which is associated with the water table. The behaviour of the electromagnetic signal differs between profiles perpendicular to the coastline (see profile 2 in Figure 2): near the sea (in the extreme SW of profile 2), already in the saltwater zone, the signal fails to penetrate deeply. This shows the geometry of the freshwater-saltwater interface and is, therefore, highly useful for deducing saltwater intrusions in the continental zone. Additionally, in this profile, the erosive surface that separates the current beach and the Holocene progradational unit is detected.
In general, washed sands and gravels are the predominant lithologies in the sector, with little or no proportion in the pelitic matrix, which allows for electromagnetic wave transmission with low attenuation. Laterally, the presence of mud, deposited in enclosed and low-energy environments, can justify the signal attenuation [24].
Profile 25 is recorded in Unit H2 sensu Goy et al. [9]. The GPR data show the two stacked planar sets, approximately 1.5-2 m in power each, separated by reactivation surfaces between layers with different slopes. Figure 2 shows the interval between 500 and 600 m of profile 25 in which these structures are clearly visible. In this profile, a very gentle slope is detected on the stoss side, with higher slopes in the slipface zone. Inside these structures, the planar cross-bedding has a high angle in the lower body and a somewhat low angle in the upper position, migrating towards the southeast in all cases (toward the current coast). These sedimentary structures are likely equivalent to those described in the sedimentological exposure of beach ridge deposits [1,2,18]. These prograding beach facies (pf in Figures 2-4) have also been described in studies using GPR [19,24,25]. Towards the top, separated by an erosion surface, subhorizontal laminations likely correspond to the topset zone. Hyperbola associated with the roots of the current groundcover and with surface bioturbation are frequently found in the most superficial section of the profile.
These bodies of unit H2 are interpreted as southeastward-migrating coastal progradations, associated with a regressive period, equivalent to that described in other sectors of the coast in SE Spain [9] or in adjacent infralittoral zones [11,12]. Minor progradational units, which are deposited in topographically lower areas, are identified in the distal section of unit H2, from 660 to 760 m in profile 25: Figure 2 highlights the different vertical positions of two of these units, which highlights the decrease in accommodation space (at metre 700 in profile 25; Figure 2). These structures are prograding wedges characterised by downlap in the lower boundary and toplap in the upper boundary resulting from small sea-level falls (approximately 0.3-0.5 m) and mild storm activity (some erosion tracks are identified). Locally, convex reflections, reaching inclinations from 0 • to 6 • , have been interpreted as berm ridge facies (bf in Figure 2) sensu Montesa et al. [24]. On the top of this facies, the landward terminations grade downlap to onlap. Toward the sea, the bf facies are limited by an erosive surface ( Figure 2).
The last section of the profile (from 750 to 790 m) shows levels with high dips and high amplitudes intersecting the previous record, which could be interpreted as erosion tracks, equivalent to those described by Baumann et al. [31] in coastal environments. Bodies with onlap morphologies are deposited on these surfaces, which could correspond to a small event of rising stage deposits associated with the base of unit H3. Subsequently, from metre 785 of profile 25, the electromagnetic signal is considerably attenuated. Therefore, the information provided is very scarce.
Profile 2, perpendicular to the previous profile, shows no cross-bedding, reflecting a virtually null dip. Towards the top and separated from an erosion surface, parallel laminations deeping seaward appear. Theses facies are interpreted as comprising sedimentary structures associated with the current foreshore ( Figure 2). Units H3, H4, H5 and H6 have morphologies equivalent to those described for unit H2. Figure 3 shows profile 27, which was recorded in the distal section of unit H3. As in the previous case, this profile highlights minor units with prograding wedges, which show downlap in the lower boundary and toplap in the upper boundary. These bodies could be associated, also in this case, with small regressive pulses at the end of unit H3 (from 150-250 m in Figure 3), with small sea-level falls, as shown by the stepwise decrease in elevation of these bodies on the radargram (Figure 3). This event could be correlated with the gap at the boundary of units H3 and H4 detected in other sectors of SE Spain [6,9], which dates to 2700 years BP according to those authors. In turn, this event could be correlated with one of the six major Holocene changes in vegetation identified by Jalut et al. [32] in the Mediterranean area, from southeastern France to southeastern Spain. These authors consider that this event corresponds to an aridification phase that occurred around 2850-1730 cal BP. This arid episode was a regional response to more global climatic changes and determined the changes in the vegetation cover [32].
From metre 270 of profile 27, an erosional surface is detected that shows onlap, which could correspond to rising stage deposits in the base of unit H4. Berm ridges facies are interpreted (br in Figure 3). The top of the profile shows subhorizontal laminations, separated from unit H3-H4 by an erosional surface. This upper set corresponds to the topset zone. Hyperbola associated with the roots of the current groundcover are also frequently found.
Profile 28 records the first metres of unit H4 (Figure 3), including erosion tracks that are decametric in lateral length and metric in depth and which could correspond to large erosion surfaces associated with the base of channels, equivalent to those described in other coastal sectors [3,33,34], which could be related to increased strong intensity. These events could be correlated with those described in south Spain, both on the Mediterranean coast [9] and on the Atlantic coast [35] at approximately 2400 years BP. Beyond metre 100 in this profile, the electromagnetic signal lacks quality, showing considerable attenuation and disturbance that could be associated with an increase in the pelitic fraction, as described in similar sedimentary environments in other sectors [36].
The boundary between units H4 and H5 was drawn at metre 200 of profile 29 (Figure 4), showing a relatively high body, with erosion tracks, on which there are onlap deposits that are most likely rising stage deposits. This H4/H5 boundary has been dated to 1900 years BP in other regions of the SE of Spain [9]. In the same profile, from metre 300, H5 progradational facies occur (pf in Figure 4) and show stacking of two planar sets separated by reactivation surfaces between layers with different slopes that are very similar to those of unit H2 and which could be interpreted as sedimentary structures typical of prograding beach ridges in a regressive context. These deposits could be associated with a stage of highstand deposits.
Morphologies similar to those described in the boundary between units H4 and H5 are also observed in the boundary between unit H5 and H6 in profile 32 from 150 to 170 m: unit H6 starts at approximately 160 m and has an internal structure equivalent to the previous structures.
The morphology of its base shows erosion and stacking of bodies migrating from a relatively high area.
Sigmoidal landward dipping reflections have been detected in the H6 units (wf in Figure 4). In other regions, equivalent radar facies have been associated with the active migration of washover fans in the back-ridge environment [24,25,30,37], in relation to storms or tsunami events. Laterally, sub-horizontal continuous reflections (landward dipping) appear (wcf in Figure 4). This facies has been interpreted as washover channel deposits, with the erosion and transference of sediments between the beach and the back-barrier domain [24,30,33].

Conclusions
The GPR profiles recorded with the 250-MHz shielded antenna provide a good balance between penetration depth and resolution for these sandy coastal plain sediments. This method differentiates the southeastward progradation of sedimentary bodies with planar low-angle cross-bedding separated by both foreset and topset reactivation surfaces. These sedimentary structures are typical of prograding beach face deposits in a regressive context. The profile perpendicular to the coastline shows parallel laminations dipping towards the sea associated with the current foreshore. Conversely, the water table level is identified at 3.5 m depth (varying slightly, depending on the sector), which matches the reflection present in the profiles and data measured directly in wells located in the area. Furthermore, the presence of saltwater in sediments near the coastline rapidly attenuates the signal, thereby preventing its penetration in the subsoil. Thus, the GPR is also a good tool to detect both water tables and saltwater intrusions in shallow aquifers.
In the GPR data, we identified a zone with strong erosional processes that can be correlated with an event 2400 years BP that marks the boundary between the two phases of progradation previously described by other authors both in the coastal plain and in infralittoral prograding wedges.
The GPR method can also be used to differentiate mid-scale units (H2, H3, H4 and H5). The distal part of each of the progradational units described in the coastal plain shows small regressive cycles represented by prograding wedges (minor scale units), deposited at lower elevations. These were likely the result of small sea-level falls that may be correlated with the periods of increased aridity described by other authors on the coast of SE Spain.
Highly advanced erosional surfaces, detected towards the base of the mid-scale units, could be associated with increases in storm intensity, combined with small sea level rises. Related to high-wave energy and the associated storm, washover deposit appears. These deposits are characterised by gently landward dipping, sub-parallel stratification.