A Potential Beach Monitoring Based on Integrated Methods

: This study focuses on the analysis of sandy beaches by integrating sedimentological, geomorphological, and geophysical investigations. The beach represents an extremely variable environment where different natural processes act simultaneously with human activities, leading to the gathering of different methodologies of the Earth Sciences to study its evolution in space and time. The aim of this research is to propose a potential procedure for monitoring the morpho-sedimentary processes of sandy beaches by analyzing the textural and compositional characteristics of the sands and quantifying the volumes involved in the coastal dynamics. The study area includes two Apulian sandy beaches (Torre Guaceto and Le Dune beach) that are representative of the coastal dynamics of a large sector of the central/northern Mediterranean Sea involving the southern Adriatic Sea and the northern Ionian Sea. Sedimentological and ecological investigations allowed to describe the textural and compositional characteristics of the beach sands by interpreting their sand provenance and the physical/biological interactions within the beach. The topographic surveys carried out with a Terrestrial Laser Scanner and an Optical Total Station, aimed to quantify the variations of sediment volume over time, whereas the Delft3d software was applied to analyze the effects of the dominant wave motion on the sedimentary dynamics. Lastly, the geophysical techniques which included Sub Bottom Proﬁler procedures, Ground Penetrating Radar investigation, and resistivity models enabled us to calculate the sand sediment thickness above the bedrock.


Introduction
Coastal zones can be defined as complex natural ecosystems where hydrodynamic, sedimentary, morphological, and biological conditions and human disturbance interact at different spaces and time scales [1][2][3].Coasts are often seen as fundamental resources to be "exploited", especially for touristic and economic purposes.However, social interests, economic investments, and the protection of natural ecosystems must meet the requirements of Integrated Coastal Zone Management (ICZM), which considers the fragility of coastal ecosystems and landscapes.
As reported from the latest United Nations reports, about 37% of the world's population lives within 100 km of the coast.This area is settled by a wide range of activities and infrastructures that have significantly changed the natural response of littoral sectors to extreme events and sea level rise.Indeed, the scientific community agrees on the increasing exposure of coastal areas to the hazard of erosion due to climate change [4,5] intensifying coast vulnerabilities [6][7][8][9][10], and the attention of decision-makers to find sustainable solutions for their protection, disaster risk reduction, and building community resilience [11][12][13][14].
The most threatened regions are Europe with 86% and Asia with 69% of their coastal ecosystems at risk [15].In Europe, most of the impact zones (15,100 km) are actively retreating, some of them despite coastal protection works (2900 km), and another 4700 km have become artificially stabilized [16].Only the magnitude and nature of erosion change from place to place.
Several European research projects (Eurosion, Conscience, Micore, Peseta, Coastgap) focus on sustainable management along Mediterranean coastal areas.These littoral sectors are mainly wave-dominated and the effect of tides on sedimentation is negligible (microtidal regime).Mediterranean sandy beaches are characterized by terrigenous sediments coming from delta deposits and/or from the erosion of cliffs or headlands occurring in the same coastal site.However, a significant percentage of sandy beaches are constituted by carbonate bioclast components that derive from shells or fragments of the organisms that populate the proximal marine environment of the Mediterranean area [17][18][19].Bioclasts become of primary importance in a geological context characterized by the absence of rivers and deltas.
Apulian sandy beaches (Southern Italy) represent a significant Mediterranean example of beach sand composition variability.The composition is connected to the typical coastal geomorphological and sedimentological features in karst areas, often lacking important fluvial courses capable of transporting large quantities of sediments into the sea [20].At the regional level, many studies analyze the features of Apulian sandy beaches (Figure 1) also showing the processes concerning their susceptibility to erosion.From a granulometric and mineralogical point of view, beach sediments show significantly different features between the Ionian and Adriatic side [21].Sediments are essentially bioclastic on the Ionian side as the littoral sector of Porto Cesareo (Figure 1 [20,22]) and terrigenous on the Adriatic Seaside [23][24][25][26] as it is detected along Torre Mileto, Rosa Marina, Pilone, Torre Canne, and Alimini coastal sectors (Figure 1 [17,21,27]).These beaches are mainly composed of fine-to coarse-grained sands characterized by carbonates (either lithoclasts or bioclasts), quartz, and other minerals such as pyroxene, amphibole, and feldspar [18].In terms of evolutionary tendency, stable conditions are recorded along the Ionian coast [28][29][30], whereas retreating tendencies are detected along the Adriatic littoral stretch [17,18,25,[31][32][33][34][35][36][37].
Evaluating the evolutionary tendency of a beach (accretion/erosion) requires various monitoring techniques, applied to the emerged and submerged sectors, in order to quantify the volume involved in coastal dynamics.In the specific literature, primary studies aimed at the description of sand movements along beaches [38], the recognition of sand textural parameters as a tool to evaluate the health state of the beach [39], the concept of equilibrium profile [40], and the role of dissipative and reflective characteristics [41] for a morphodynamic beach classification [42,43].Later, the study of sandy beaches experienced significant progress.Table 1 represents an overview of the different studies carried out on the beach environment.The list of authors provides an example of the large number of techniques for beach investigation.

Scientific Sector
Author Methodology  For this reason, this work focuses on the mechanism of erosion, transport, and sedimentation of two Apulian sandy littoral sectors through a multidisciplinary approach.We suggest a reliable guideline for monitoring sandy beaches aiming to:
a description of the textural and compositional features of beach sands; 2.
an overall sketch scheme of beach dynamics; 4.
a measurement of the beach sediment thickness above the bedrock and therefore of the sand amount virtually involved in the beach sedimentary dynamics.
The research proposes a methodological approach that could be useful for scientific communities or private/public coastal management.After a brief description of the study areas, the main part of the research focuses on testing the procedures for beach investigation.In the final paragraphs, a suggestion regarding data management is also provided by using the result of this research as an example of data interpretation.

Study Area
The study area includes two Apulian sandy littoral stretches (Figure 2a m and 100-110 m, and it is predominantly covered by bioclastic sands.Within the inner portion of the shelf, Posidonia oceanica meadows form large patches which are replaced at depth by coralligenous platform deposits [95].Le Dune can be considered representative of the Ionian sandy beaches that have experienced a remarkable urbanization of the coastal areas in the last decades.Generally, wave directions are between 150° (15.91%) and 300° (8.47%), and seas from the south (24.48%) and southeast (15.91%) have the highest frequencies [96].A rip current phenomenon is often observed during storm events [22].Torre Guaceto (Figure 2b) is a pocket beach that stretches for about 1 km.The embayment is delimited by two rocky headlands: the south-eastern Lower Pleistocene shallow-water skeletal calcarenites (Punta Penna Grossa outcrop) and an anthropogenic structure of the 20th century in the north-westward [93].The area is characterized by high coastal dunes and a rocky and sandy sea bottom.The beach is affected by a slight human impact, and it can be representative of sandy beaches lying on the northern coastal sector of Brindisi in terms of wind and wave weather, sand composition, and granulometric features.Beach sediments range from coarse to fine sands mainly composed of carbonate grains (bioclasts or lithoclasts) and siliciclastic minerals.Wave directions are mainly between 330 • (16.79%) and 120 • (15.36%).It is also observed that the highest frequency of appearance (22.15%) and the highest significant wave heights (Hs between 3 m and 4 m) derive from the 0 • direction [94].
On the contrary, Le Dune beach (Figure 2c) is located along the Ionian Seaside and stretches for about 800 m.The small embayment has a slight Z-shaped form with an increased sediment accumulation on its western limit.Beach sediments range from very coarse to medium-fine sands and they are mainly made up of bioclast fragments.The pocket beach is bounded landward by aeolian dunes and by promontories where Cretaceous limestone crops out.At deeper depths, the shelf extends between the isobaths at −5 m and 100-110 m, and it is predominantly covered by bioclastic sands.Within the inner portion of the shelf, Posidonia oceanica meadows form large patches which are replaced at depth by coralligenous platform deposits [95].Le Dune can be considered representative of the Ionian sandy beaches that have experienced a remarkable urbanization of the coastal areas in the last decades.Generally, wave directions are between 150 • (15.91%) and 300 • (8.47%), and seas from the south (24.48%) and southeast (15.91%) have the highest frequencies [96].A rip current phenomenon is often observed during storm events [22].

Materials and Method
Since this work is based upon a multidisciplinary approach, some of the main techniques used in the field of sedimentology, petrography, biogeology-paleontology, geomorphology, and applied geophysics were used for beach investigation.However, according to the objectives to be achieved, some procedures were constructed within the current research in order to obtain a detailed outline of beach dynamics and sediment thickness quantification.
3.1.Beach Sand Analysis 3.1.1.Sampling Sampling procedures were carried out every six months along the backshore, foreshore, and shoreface sectors (example in Figure 3).Ground control points were used as collection points to compare the texture and the composition of the sand and to record the sampling depth variations over time.Samplings were carried out from the base of the foredune until 6 m depth of the shoreface (the local storm-wave base in the Adriatic and the Ionian Sea).Along the shoreline, the collecting points were spaced at around 100-200 m from each other, whereas shoreface samples were collected at each meter depth through diving techniques along two or three cross-shore transects.Around 300 gr of sand was collected between 0 and 2 cm depth down from the water/sediment interface by following the standard sampling procedure for marine sediments proposed by [97].
southernmost sector and a developed dune environment on the northwestern sector.

Materials and Method
Since this work is based upon a multidisciplinary approach, some of the main techniques used in the field of sedimentology, petrography, biogeology-paleontology, geomorphology, and applied geophysics were used for beach investigation.However, according to the objectives to be achieved, some procedures were constructed within the current research in order to obtain a detailed outline of beach dynamics and sediment thickness quantification.

Sampling
Sampling procedures were carried out every six months along the backshore, foreshore, and shoreface sectors (example in Figure 3).Ground control points were used as collection points to compare the texture and the composition of the sand and to record the sampling depth variations over time.Samplings were carried out from the base of the foredune until 6 m depth of the shoreface (the local storm-wave base in the Adriatic and the Ionian Sea).Along the shoreline, the collecting points were spaced at around 100-200 m from each other, whereas shoreface samples were collected at each meter depth through diving techniques along two or three cross-shore transects.Around 300 gr of sand was collected between 0 and 2 cm depth down from the water/sediment interface by following the standard sampling procedure for marine sediments proposed by [97].

Grain-Size Analysis
The grain-size analyses were carried out by using the standard procedures provided by the American Society for Testing and Materials and the British Standard.For the sieving, a set of ASTM sieves with meshes of ½ phi from 2 mm to the minimum granulometric fraction (<0.125 mm) was used.The grain fractions with diameters less than 0.062 mm were excluded from the analysis because they were represented by less than 1-2%.In the

Grain-Size Analysis
The grain-size analyses were carried out by using the standard procedures provided by the American Society for Testing and Materials and the British Standard.For the sieving, a set of ASTM sieves with meshes of 1  2 phi from 2 mm to the minimum granulometric fraction (<0.125 mm) was used.The grain fractions with diameters less than 0.062 mm were excluded from the analysis because they were represented by less than 1-2%.In the laboratory, samples were dried in the oven at a temperature of 80 • for 24 h and each individual sample was quartered and set in a sieve column.The sand sediments from 2.0 mm to 0.125 mm were sieved with the vibrating screen for a duration of 20 min.Subsequently, each held fraction was weighed and the results were processed with a specific application for Microsoft Excel (Gradistat© v8, Figure 4), which yield distribution cumulative curves, histograms, and statistically evaluate the main textural parameters: mean size (Mz), sorting (σ), skewness (S k ), and Kurtosis (K G ) (Figure 4).Lastly, a seasonal comparison of the granulometric parameters was carried out in order to evaluate their temporal changes.
to 0.125 mm were sieved with the vibrating screen for a duration of 20 min.Subsequently, each held fraction was weighed and the results were processed with a specific application for Microsoft Excel (Gradistat© v8, Figure 4), which yield distribution cumulative curves, histograms, and statistically evaluate the main textural parameters: mean size (Mz), sorting (σ), skewness (Sk), and Kurtosis (KG) (Figure 4).Lastly, a seasonal comparison of the granulometric parameters was carried out in order to evaluate their temporal changes.

Compositional Analysis and Bioclasts Quantification/Recognition
The most frequent size class of the statistical distribution was investigated through a binocular optical microscope (Figure 5a).The percentages of the main constituents of sands were evaluated to obtain quantitative and qualitative information in terms of petrographic composition.Since the investigation was carried out on Apulian beach sands, which generally include carbonates (lithoclasts or bioclasts) and siliciclastic minerals, the sands were analyzed by considering the diagram proposed by [98,99].This method provides the recognition of three main classes: CE carbonate extrabasinal (lithoclasts), CI carbonate intrabasinal (bioclasts), and NCE non-carbonate extrabasinal grains.In particular, the sediment fraction was spread on a rigid base divided into five fields (Figure 5b).Each field was analyzed in detail under the microscope (Figure 5c) and for each field, the particles belonging to the NCE, CE, and CI class were counted.Once the counting operation for each field was completed, the results obtained were summed and the percentage of each component was calculated and inserted in the final classification diagram (Figure 5d).

Compositional Analysis and Bioclasts Quantification/Recognition
The most frequent size class of the statistical distribution was investigated through a binocular optical microscope (Figure 5a).The percentages of the main constituents of sands were evaluated to obtain quantitative and qualitative information in terms of petrographic composition.Since the investigation was carried out on Apulian beach sands, which generally include carbonates (lithoclasts or bioclasts) and siliciclastic minerals, the sands were analyzed by considering the diagram proposed by [98,99].This method provides the recognition of three main classes: CE carbonate extrabasinal (lithoclasts), CI carbonate intrabasinal (bioclasts), and NCE non-carbonate extrabasinal grains.In particular, the sediment fraction was spread on a rigid base divided into five fields (Figure 5b).Each field was analyzed in detail under the microscope (Figure 5c) and for each field, the particles belonging to the NCE, CE, and CI class were counted.Once the counting operation for each field was completed, the results obtained were summed and the percentage of each component was calculated and inserted in the final classification diagram (Figure 5d).Furthermore, the whole shells and fragmented particles were isolated and analyzed by mean of a binocular optical microscope in more detail (example in Figure 6).In this context, a set of tweezers was used to manually separate the fragmented and unrecognizable fraction from the rest of the sample.The percentage was defined among three Furthermore, the whole shells and fragmented particles were isolated and analyzed by mean of a binocular optical microscope in more detail (example in Figure 6).In this context, a set of tweezers was used to manually separate the fragmented and unrecognizable fraction from the rest of the sample.The percentage was defined among three classes: 0-30%, 30-60%, 60-90%.The whole shells were analyzed and separated on the base of their Phylum to provide a first classification.The class and genus of shells were also evaluated in the case of foraminifers.

Wave Simulations and Hydrodynamic Model
Software Delft3D was used to analyze the correlation between sedimentological parameters and wave processes.As shown in Figure 7, GRID, FLOW, and WAVE modules were applied to calculate the hydrodynamic flow and to simulate the wave motion.In particular, coupled simulations were performed by means of Delft3D-WAVE with Delft3D-FLOW to incorporate wave-current interaction.Delft3DFLOW was used to assess the hydrodynamic flow, and Delft3D-WAVE was used to simulate the wave propagation,  The main bioclast classification was carried out on the first sampling season.Successively, the bioclast evaluation was conducted by analyzing any differences in the Phylum already identified.
It is essential to specify that the recognition of the nature of sands was relatively simpler than expected.It recorded a maximum of 3% of grains with uncertain nature mainly related to the distinction between lithoclasts and bioclasts without a recognizable internal structure.The sandy particles with uncertain origin were not considered, while fossils were included in the class of lithoclasts as they were considered part of older deposits.

Wave Simulations and Hydrodynamic Model
Software Delft3D was used to analyze the correlation between sedimentological parameters and wave processes.As shown in Figure 7, GRID, FLOW, and WAVE modules were applied to calculate the hydrodynamic flow and to simulate the wave motion.In particular, coupled simulations were performed by means of Delft3D-WAVE with Delft3D-FLOW to incorporate wave-current interaction.Delft3DFLOW was used to assess the hydrodynamic flow, and Delft3D-WAVE was used to simulate the wave propagation, based on the SWAN (Simulating WAves Nearshore) model.The model was based on two main datasets: bathymetry and climatic information (waves and wind data).The bathymetric data were provided by the EMODnet website (https://portal.emodnet-bathymetry.eu/, accessed on 1 November 2022).The wave data were acquired from SIT Puglia website (http:// cartografia.sit.puglia.it/download/PRC/Relazione%20Generale_allegato711.pdf,accessed on 30 October 2022) and from a wave buoy located in the Ionian Sea (Datawell Directional

The Terrestrial Laser Scanner (TLS)
The topography of the emerged beach sector was investigated using a Terrestria ser Scanner (TLS) that allows constructing 3D elevation models and quantifying the m phological changes over time.In particular, the instrument acquires the spatial coo nates of many points by measuring the distance between the TLS and the object of s in order to build point clouds and successive three-dimensional elevation models.repeated survey campaigns were carried out during the summer of 2018 and the sum of 2019 by using a Riegl VZ-400 laser scanner, which is characterized by a theoretical ra of 400 m.The instrumental accuracy is ±0.003 m for 50 m (example in Figure 8a).
Due to the length of both beaches, several scans were carried out from four stat with a distance of about 250 m from each other.The point clouds were processed thro different stages: scans registration; multi-scan adjustment and georeferencing; 3D p cloud cleaning; triangulation (mesh) and Digital Terrestrial Model DTM creation w size cell of 50 cm.The Riscan Pro and CloudCompare software were used for point c processing (Figure 8b), filtering, and rasterization, while the elevation correction and comparison of the seasonal DTMs were carried out with ArcMap © 10.1 (Figure 8c).
As the sampling procedure, the TLS acquisition was carried out every six month detect the main morphological changes in the emerged sectors of the beach.The topography of the emerged beach sector was investigated using a Terrestrial Laser Scanner (TLS) that allows constructing 3D elevation models and quantifying the morphological changes over time.In particular, the instrument acquires the spatial coordinates of many points by measuring the distance between the TLS and the object of study in order to build point clouds and successive three-dimensional elevation models.Two repeated survey campaigns were carried out during the summer of 2018 and the summer of 2019 by using a Riegl VZ-400 laser scanner, which is characterized by a theoretical range of 400 m.The instrumental accuracy is ±0.003 m for 50 m (example in Figure 8a).
Due to the length of both beaches, several scans were carried out from four stations with a distance of about 250 m from each other.The point clouds were processed through different stages: scans registration; multi-scan adjustment and georeferencing; 3D point cloud cleaning; triangulation (mesh) and Digital Terrestrial Model DTM creation with a size cell of 50 cm.The Riscan Pro and CloudCompare software were used for point cloud processing (Figure 8b), filtering, and rasterization, while the elevation correction and the comparison of the seasonal DTMs were carried out with ArcMap © 10.1 (Figure 8c).
As the sampling procedure, the TLS acquisition was carried out every six months to detect the main morphological changes in the emerged sectors of the beach.

The Optical Total Station (OTS)
With regards to the shoreface investigation, two profile surveys were performed using the Optical Total Station (OTS) Leica TS15 to collect seasonal bathymetric information.The OTS is composed of an electromagnetic distance measure instrument and electronic theodolite that is also integrated with a computer storage system (Figure 9).The instrument allows the measuring of horizontal and vertical angles as well as the sloping distance between the object and the instrument.

The Optical Total Station (OTS)
With regards to the shoreface investigation, two profile surveys were performed using the Optical Total Station (OTS) Leica TS15 to collect seasonal bathymetric information.The OTS is composed of an electromagnetic distance measure instrument and electronic theodolite that is also integrated with a computer storage system (Figure 9).The instrument allows the measuring of horizontal and vertical angles as well as the sloping distance between the object and the instrument.
During the seasonal field surveys, the measurements were performed along two cross-shore transects in georeferenced points following the direction of the sampling profiles.About 20 bathymetric points were measured for each transect from the shoreline to a 4 m depth on both beaches.During the seasonal field surveys, the measurements were performed along two crossshore transects in georeferenced points following the direction of the sampling profiles.About 20 bathymetric points were measured for each transect from the shoreline to a 4 m depth on both beaches.

Ground Penetrating Radar (GPR)
The GPR investigation was applied for evaluating the sediment thicknesses in the emerged beach sector.The choice of this technique depends on the dynamic nature of this sub-environment that constantly changes with the amount of available material for transport and sedimentation.Generally, these variations are recorded within the sand deposits and particularly in the dune environment, whose stratigraphic horizons can be explored with the GPR technique.The instrument operation is based on the input of high frequency electromagnetic pulses deriving from an antenna, which is placed in contact with the surface to be investigated.The depth of the investigation depends on the frequency used and therefore on the type of antenna.As the signal frequency increases, there is a growth in data resolution but a decreasing in the investigation depth.

Ground Penetrating Radar (GPR)
The GPR investigation was applied for evaluating the sediment thicknesses in the emerged beach sector.The choice of this technique depends on the dynamic nature of this sub-environment that constantly changes with the amount of available material for transport and sedimentation.Generally, these variations are recorded within the sand deposits and particularly in the dune environment, whose stratigraphic horizons can be explored with the GPR technique.The instrument operation is based on the input of high frequency electromagnetic pulses deriving from an antenna, which is placed in contact with the surface to be investigated.The depth of the investigation depends on the frequency used and therefore on the type of antenna.As the signal frequency increases, there is a growth in data resolution but a decreasing in the investigation depth.
A survey campaign through a 200 Mhz antenna was carried out at Torre Guaceto and Le Dune beach (Figure 10a).The choice of this frequency depended on some tests carried out previously with 400 and 40 Mhz in order to detect the most significant radar stratigraphy section (Figure 10c).The data acquisition was performed along a cross-shore transect from the swash zone to the dune environment (Figure 10b).
With regards to the data processing, the Reflexw software was used to analyze and elaborate the 2D radar stratigraphic sections.A survey campaign through a 200 Mhz antenna was carried out at Torre Guaceto and Le Dune beach (Figure 10a).The choice of this frequency depended on some tests carried out previously with 400 and 40 Mhz in order to detect the most significant radar stratigraphy section (Figure 10c).The data acquisition was performed along a cross-shore transect from the swash zone to the dune environment (Figure 10b).
With regards to the data processing, the Reflexw software was used to analyze and elaborate the 2D radar stratigraphic sections.

Sub Bottom Profiler (SBP)
As the GPR investigation, the sub bottom profiler was carried out to quantify the sand thicknesses within the submerged sector of the beach.This technique is widely applied in the field of marine geology.The SBP exploits the elastic properties of the ground to reconstruct the stratigraphic succession of deposits occurring below the seabed.Each surface that marks a lithological transition or any acoustic impedance such as the water/sediment passage represents an elastic discontinuity capable of reflecting part of the seismic energy.The reflected signal is received by a transducer and sent to the visualization program by creating a seismic section.The signal penetration and reflection depend on the frequency and the physical properties of the sediment.

Sub Bottom Profiler (SBP)
As the GPR investigation, the sub bottom profiler was carried out to quantify the sand thicknesses within the submerged sector of the beach.This technique is widely applied in the field of marine geology.The SBP exploits the elastic properties of the ground to reconstruct the stratigraphic succession of deposits occurring below the seabed.Each surface that marks a lithological transition or any acoustic impedance such as the water/sediment passage represents an elastic discontinuity capable of reflecting part of the seismic energy.The reflected signal is received by a transducer and sent to the visualization program by creating a seismic section.The signal penetration and reflection depend on the frequency and the physical properties of the sediment.
In this study case, one marine field survey was carried out for each study area by applying low signal frequencies (85-115 kHz).The investigation was organized by following the same cross-shore transect of the georadar investigation from 6-2 m depth (example in Figure 11).
Moreover, a boat capable of reaching shallow waters was utilized for the investigation and the instrumentation was characterized by an Innomar SBP SES 2000 Compact system.The SBP transducer was installed on one side of the boat (Figure 11a), while the hardware system and the data management laptop were positioned inside the boat.Moreover, the whole system was connected to a GPS RTK, and interfaced with the instrument management software to acquire the exact data location.Regarding data acquisition, SES WIN software was utilized to obtain seismic sections in real-time.The program produces a *.SES file which displays the width of the acoustic signal reflected at different depths; it also contains data relating to the longitude, latitude, depth of the seabed, and acquisition parameters.All these data were transferred to Delph Seismic software to analyze the 2dimensional stratigraphic cross sections, define the bedrock location, and obtain strata thicknesses.In this study case, one marine field survey was carried out for each study area by applying low signal frequencies (85-115 kHz).The investigation was organized by following the same cross-shore transect of the georadar investigation from 6-2 m depth (example in Figure 11).

Geoelectric Investigation (ERT)
A geoelectric investigation was applied to analyze the entire beach profile and to gain data about the foreshore subenvironment.The application of a geoelectric investigation provides information relating to the subsoil electrical resistivity distribution.
As shown in Figure 12, the ERT was carried out across the same GPR cross-shore transect.The investigation profile was extended until 3 m depth in the shoreface to reach the SBP profile (Figure 12c).The measurements were carried out by adopting an IRIS Instruments Syscal Pro Georesistive Meter, two multi-electrode cables at 24-channels with 3 m interelectrode spacing for a total length of 141 m.Proceeding from NNE to SSW, the Moreover, a boat capable of reaching shallow waters was utilized for the investigation and the instrumentation was characterized by an Innomar SBP SES 2000 Compact system.The SBP transducer was installed on one side of the boat (Figure 11a), while the hardware system and the data management laptop were positioned inside the boat.Moreover, the whole system was connected to a GPS RTK, and interfaced with the instrument management software to acquire the exact data location.Regarding data acquisition, SES WIN software was utilized to obtain seismic sections in real-time.The program produces a *.SES file which displays the width of the acoustic signal reflected at different depths; it also contains data relating to the longitude, latitude, depth of the seabed, and acquisition parameters.All these data were transferred to Delph Seismic software to analyze the 2-dimensional stratigraphic cross sections, define the bedrock location, and obtain strata thicknesses.

Geoelectric Investigation (ERT)
A geoelectric investigation was applied to analyze the entire beach profile and to gain data about the foreshore subenvironment.The application of a geoelectric investigation provides information relating to the subsoil electrical resistivity distribution.
As shown in Figure 12, the ERT was carried out across the same GPR cross-shore transect.The investigation profile was extended until 3 m depth in the shoreface to reach the SBP profile (Figure 12c).The measurements were carried out by adopting an IRIS Instruments Syscal Pro Georesistive Meter, two multi-electrode cables at 24-channels with 3 m interelectrode spacing for a total length of 141 m.Proceeding from NNE to SSW, the first cable (steel electrodes 1-24) was placed on the subsoil of the emerged beach sector, while the second cable (graphite electrodes 25-48) was placed on the seabed (Figure 12a).The data were acquired both in Wenner-Schlumberger and dipole-dipole configurations.The former provides high vertical resolution while the latter high lateral resolution-although it is characterized by a lower signal-to-noise ratio than other devices [100].For this reason, reciprocal measures were also carried out to facilitate the data quality control necessary to obtain reliable and well-resolved images.After the quality control (carried out by setting a 10% threshold for repeatability errors, and standard deviation obtained from the stacks), the experimental data were inverted using the RES2DINV code (Geotomo Software, [101]).The resulting model (Figure 12d) was obtained by inverting the dataset relating to the measurements in dipole-dipole configuration, including the GPS-RTK data that were acquired along the same transept for topography correction.

Data Comparison
The first step in beach dynamics interpretation is comparing some of the results deriving from different methodologies applied in this research.Figure 13 represents an ex- Lastly, two excavations were also carried out using shovels along the GPR/ERT transects.This technique was fundamental to calibrate and verify with a direct measurement the exact location of the bedrock extrapolated by the radar stratigraphy sections and the resistivity models.

Data Comparison
The first step in beach dynamics interpretation is comparing some of the results deriving from different methodologies applied in this research.Figure 13 represents an example of a combination of Delft3D models, DTM comparison, sand classification, and mean size distribution.
This type of comparison enables to highlight the processes concerning the sediment transport and the evolutionary tendency of the beach explained by the connection among the simulation of the wave motion, erosion and accretion phases, sand composition, and granulometric features.For instance, at Torre Guaceto beach, the analysis of the emerged beach underlines the presence of a southernmost sector more exposed to the wave motion than the northern sector during storm events.The northern beach sector is characterized by an erosional tendency, especially at the dune base as reported from the DTM investigation (Figure 13a).Moreover, a significant composition variability is detected along the foreshore samples.The sand of Torre Guaceto beach is classified as "Hybrid intrabasinal sand", but the northern-central sector of the emerged beach is characterized by a higher content of bioclast than the southern sector, which instead is richer in siliciclastic content (Figure 13b).In addition, according to the mean size data, an increasing trend of the grain diameter is recorded from the samples collected close to the outermost promontories towards the central part of the beach (Figure 13c), suggesting a bioclast longshore transport in this direction.Indeed, looking at the main wave direction, the more porous and lighter bioclastic sediments are removed from the southern sector and transported to the central part of the beach.Indirectly, the transport of the bioclast component increases the percentages of the siliciclastic minerals in the southern sector, explaining the sudden longshore compositional variation of the emerged sector of Torre Guaceto beach.Comparing this result with the DTM, we found that the sector with the higher bioclast component (northern central sector) represents the more stable part of the beach in terms of evolutionary tendency, whereas the southernmost sector is eroding.

Beach Dynamics
Another example of data integration can be performed by using the results shown in Figure 14.The main sand characteristics (texture and petrography) coupled with the study of the meteorological events and the interpretation of the sediment transport allows to reconstruct the sedimentary balance and propose an overall sketch scheme of beach dynamics (Figure 14d).In the semi-closed coastal system of Le Dune beach, considering the southerly seas as the main wave direction, the scheme shows how the sediment inputs mainly derive from local sources and the sediment outputs involve small amounts of sediments.The gain sediment includes material coming from the erosion of the dunes, rocky shoreface, and lateral headlands as well as the bioclastic sediment transported from offshore.Therefore, the sedimentary dynamics are controlled by a predominant accretion of sand diffusion by aeolian processes especially in the westernmost part of the coast; a cross-shore sediment transport; a nearshore rip circulation characterized by a longshore sediment transport converging in the middle of the embayment; a weak lateral sediment interchange with adjacent littoral sectors and an offshore sand dispersion during storm events.The current morphological configuration of the sea bottom significantly influences the cross-shore sediment transport.The rocky headlands located in the shallow water and the rocky sea bottom make the beach an almost completely closed bay, which is also characterized by the absence of river basins and the presence of strong rip currents that carry large amounts of sediments towards the offshore during storm events.

Sediment Thickness and Beach Profile
One of the main applications deriving from the methodologies tested in this research is the combination of geophysics techniques.This operation is fundamental to reconstruct the entire beach profile from the dune environment to the offshore and to extrapolate the bedrock location.In particular, the merging of the radar section, resistivity model, and sub-bottom profile allows to quantify the sediment thickness within the entire beach environment and to record the location of the sub-environments.For instance, in the case of Torre Guaceto, the beach profile (Figure 15) extends for about 700 m.It represents the connection of the GPR, resistivity model, and SBP navigation transect.The bedrock location ranges between +1 m and −2 m within the emerged beach sector and from −2 m to −10 m in the submerged sector.
Figure 15.Example of beach profile with bedrock location deriving from the merging of GPR + ERT + SBP results.The orange dot line in the map depicts the investigation gap.Within the resistivity model and the GPR profile, the distance axis has been exaggerated for a higher result resolution.As ERT and GPR measurements were performed along the same transect of the emerged beach sector, the figures were overlapped to compare the same bedrock location extrapolated from both models.The yellow area between the ground surface/sea bottom and the bedrock represents the beach sand sediments.

Main Findings
The main outcomes highlight and confirm the literature data about the sand composition variability and the morpho-sedimentary dynamics occurring along the Adriatic and Ionian littoral.The findings represent the features of the Mediterranean coastal types.
The results underline the need to improve methodologies for pocket beach investigation which are difficult to preserve due to their low geological records.They are affected by seasonal changes strictly correlated with the main sediment source.The sediment transport mainly occurs between the dune environment and the offshore characterized by the presence of Posidonia oceanica meadows.In this context, carbonate factory protection and dune restoration measures are crucial to preserve this geological environment where the anthropic pressure is increasingly growing, especially for tourist activities.
The analysis of Torre Guaceto beach provides a significant example of improvement studies in ancient research.Environments at the same depth could have different compositions due to the main sediment selection caused by the current beach dynamics system.Indeed, the qualitative and quantitative changes of the bioclast content along the crossshore result in problematic stratigraphic correlations between adjacent logs in ancient successions located in similar settings.Moreover, this case study points out the presence of distinct dynamics due to the different exposure of the beach sector to wave motion providing information about eroding and prograding phases within the same environment.
The case study of Le Dune beach underlines the existence of many pocket beaches in the world that are only preserved thanks to the bioclast provenance source and the dune environment.In this case, the growth of urbanization, especially along the southern coastal sector, could significantly impact on the amount of sediment available.Although this beach is part of a marine protected area, the lack of collaboration between local decision-makers and the scientific community about the knowledge of pocket beach dynamics has led to the destruction of the dune environment with a consequent establishment of a The most relevant reflector resolution (bedrock) occurs below the sediments characterizing the dune environment, the sandy upper shoreface, and beneath the sand bar located between the lower shoreface and the offshore.Moreover, the bedrock visibly emerges in the rocky lower shoreface.
Considering the geology of the study area, the bedrock is certainly associated with the Pleistocene biocalcarenite (packstone and grainstone texture).

Main Findings
The main outcomes highlight and confirm the literature data about the sand composition variability and the morpho-sedimentary dynamics occurring along the Adriatic and Ionian littoral.The findings represent the features of the Mediterranean coastal types.
The results underline the need to improve methodologies for pocket beach investigation which are difficult to preserve due to their low geological records.They are affected by seasonal changes strictly correlated with the main sediment source.The sediment transport mainly occurs between the dune environment and the offshore characterized by the presence of Posidonia oceanica meadows.In this context, carbonate factory protection and dune restoration measures are crucial to preserve this geological environment where the anthropic pressure is increasingly growing, especially for tourist activities.
The analysis of Torre Guaceto beach provides a significant example of improvement studies in ancient research.Environments at the same depth could have different compositions due to the main sediment selection caused by the current beach dynamics system.Indeed, the qualitative and quantitative changes of the bioclast content along the crossshore result in problematic stratigraphic correlations between adjacent logs in ancient successions located in similar settings.Moreover, this case study points out the presence of distinct dynamics due to the different exposure of the beach sector to wave motion providing information about eroding and prograding phases within the same environment.
The case study of Le Dune beach underlines the existence of many pocket beaches in the world that are only preserved thanks to the bioclast provenance source and the dune environment.In this case, the growth of urbanization, especially along the southern coastal sector, could significantly impact on the amount of sediment available.Although this beach is part of a marine protected area, the lack of collaboration between local decision-makers and the scientific community about the knowledge of pocket beach dynamics has led to the destruction of the dune environment with a consequent establishment of a building area.In this respect, higher restrictive measures should be considered for the protection of these types of environments.
Torre Guaceto and Le Dune represent two different pocket beaches in terms of length, sedimentology, composition, dynamics, and anthropic impact.For this reason, some considerations raised during the investigation.Geophysics provided more detailed results of the submerged sector at Torre Guaceto than at Le Dune beach.The rocky seabed in the latter affected the navigation along the shoreface transect requiring a major support from electrical tomography.The closed beach system, the morphology of the seabed, and the urbanization of the dune environment influenced the progress of the investigations.It would have been interesting to measure the sediment thicknesses along different crossshore transects, but the southern sector was affected by a significant urbanization.
These results are strictly correlated with a two-year beach monitoring that is only able to evaluate sudden and small-scale beach variations.However, by seasonally applying the proposed methodologies over a longer period of time, it would be possible to recognize the real evolution of our beaches.The study of the evolutionary tendency could provide a significant impact on public decision-makers and current legislation on unnecessary protection and defense works.

Final Remarks
The procedure application highlights the need of comparing data from different sectors of the Earth Sciences to enrich the amount of missing information about complex and dynamic systems such as the beach environment.The multidisciplinary approach enables to collect a wide range of data that can be connected and interpreted with each other.The data combination provides a detailed outline of the beach dynamics in qualitative and quantitative terms.
From a sedimentological point of view, the research allows expanding the knowledge about large-scale sedimentary processes by quantifying the volumes involved in the erosion, transport, and deposition dynamics and the potential preservation of specific sub-environments during the sea-level rise or still-stand periods.Furthermore, this study shows how lateral variability in sand composition (not only in grain size) is a widespread feature in present-day beaches.The lateral composition is mainly influenced by very local processes and the sudden variations can be easily misinterpreted in the geological record.
Moreover, the study of sand composition and in particular the analysis of the bioclast component represents a rising technique in coastal erosion investigation.Indeed, measuring the bioclast percentage in the beach sand allows to hypothesize strategies for provenance marine environment safeguarding and monitoring by avoiding unnecessary or too-restrictive protection measures.Lastly, the beach sediment analysis provides significant progress for compatibility studies and nourishment interventions in terms of textural and compositional features.This research tries to respond to the increasingly pressing need to overcome the conflicts between naturalistic requirements and the use of beaches for social and economic purposes, which represent one of the main land use planning issues of coastal regions.

Conclusions
This research described the processes of erosion, transport, and sedimentation of two Apulian sandy beaches through a multidisciplinary approach.Table 2 shows the methodological proposal for sandy beach analysis deriving from the multidisciplinary approach tested in this research.This procedure includes topographic, geophysical, sedimentological, and compositional analyses for monitoring and safeguarding these transitional environments that can be useful for research studies or industrial applications.As shown in Table 2, the procedure for investigating sandy beaches was based upon four main types of analyses: topographic surveys through TLS and OTS; a hydrodynamic model using Delft3D software; a series of geophysics techniques by the use of GPR, ERT and SBP; and a beach sand analysis through sedimentological and compositional applications.Each technique covered a specific beach subenvironment and provided distinct results.Therefore, the procedure could be invalid without the combination of the proposed techniques as the beach dynamics interpretation could be incomplete.
Regarding the geophysical investigation, it was essential to confirm the indirect measures with direct observations such as excavations.The resistivity studies were also correlated to the GPR and SBP surveys for a more detailed outline of the lithology and stratigraphy of the subsoil.Lastly, for the GPR application, the use of 400 and 200 Mhz antennas produced more satisfactory resolutions of the bedrock location than the 40 Mhz antenna.
It is important to emphasize that all the equipment, instrument frequencies, and resolutions were set on the aim of this research and the geological context.
Investigating sandy beaches is not undemanding due to the dynamism of this type of environment.For this reason, as shown in Table 2, the use of GPS supported most of the methodologies in order to respect the exact location of the seasonal investigation.Data comparison should be provided by taking into account the same season of the field measure in order to avoid interpretation issues.
Another suggestion that could improve this type of investigation is the use of the "Sediment transport module" in the application of Delft3D software.Indeed, the Delft3D-SED coupled with the others enables to analyze the effects of dredging on the environment, the sedimentation and resuspension of sediment, and the sand transport.Moreover, a primary phase of data collection on the marine climate of the study area is highly recommended to interpret the beach dynamics.

Figure 1 .
Figure 1.(a) Geographical location of a few sandy beaches (yellow circle) characterizing the Apulian littoral stretch.The white arrows indicate the general alongshore current influencing the sedimen transport on the Adriatic and the Ionian Sea; (b-d) Macroscopic features variability of the Apulian beach sands in terms of color and grain size (modified from[20]).

Figure 1 .
Figure 1.(a) Geographical location of a few sandy beaches (yellow circle) characterizing the Apulian littoral stretch.The white arrows indicate the general alongshore current influencing the sediment transport on the Adriatic and the Ionian Sea; (b-d) Macroscopic features variability of the Apulian beach sands in terms of color and grain size (modified from[20]).
): Torre Guaceto (40,71 lat.N; 17,77 long.E) and Le Dune beach (40,27 lat.N; 17,87 long.E) occurring in the Adriatic and Ionian sectors of the Salento Peninsula, respectively (Brindisi and Lecce provinces).Both beaches are in two of the most important Marine Protected Areas (MPA) which include 15 different habitats of the typical Mediterranean submerged populations and the presence of Posidonia oceanica meadows.They can be considered endmembers of the coastal dynamics of a large sector of the central/northern Mediterranean Sea involving the southern Adriatic Sea and northern Ionian Sea with negligible anthropogenic impact.

Figure 2 .
Figure 2. Coastal features of Torre Guaceto and Le Dune beach.(a) Geographical location of the study areas (orange star); (b) satellite image of Torre Guaceto beach showing a landward well developed dune belt; (c) satellite image of Le Dune beach showing a large, urbanized area on the southernmost sector and a developed dune environment on the northwestern sector.

Figure 3 .
Figure 3. (a) Map showing the sampling collecting points.D = shoreline sample label, BD = foredune base sample, TA-1 = Sample collected along transect A at 1 m depth; (b) Ideal submerged beach profile (not in scale) with the indication of the sampling depths; (c) collection of shoreface samples; (d) small net for marine samples; (e) whiteboard for sample marine marks.

Figure 3 .
Figure 3. (a) Map showing the sampling collecting points.D = shoreline sample label, BD = foredune base sample, TA-1 = Sample collected along transect A at 1 m depth; (b) Ideal submerged beach profile (not in scale) with the indication of the sampling depths; (c) collection of shoreface samples; (d) small net for marine samples; (e) whiteboard for sample marine marks.

25 Figure 5 .
Figure 5. Example of hybrid sand composition analysis; (a) sand analysis at the microscope; (b) sample spread on a rigid base with the field locations; (c) field for counting method; (d) example of counting for sand classification, green: bioclasts, yellow: carbonate lithoclasts; red: quartz.

Figure 5 .
Figure 5. Example of hybrid sand composition analysis; (a) sand analysis at the microscope; (b) sample spread on a rigid base with the field locations; (c) field for counting method; (d) example of counting for sand classification, green: bioclasts, yellow: carbonate lithoclasts; red: quartz.

Figure 6 .
Figure 6.Example of bioclast evaluation procedure; (a) sample bioclast separation with a set of tweezers; (b) isolated bioclasts in test tubes; (c) bioclast classification.

Figure 7 .
Figure 7. Flow chart of the processing steps followed in the application of Delft3D software.

Figure 7 .
Figure 7. Flow chart of the processing steps followed in the application of Delft3D software.

Figure 8 .
Figure 8. Emerged beach topography investigation.(a) TLS Riegl VZ-400 instrument during a survey campaign; (b) Example of a point cloud during a scan registration with Riscan Pro software: the red point represents the instrument location during the measurement acquisition; (c) Digital elevation model (DEM) of the emerged sector in summer 2018 processed with ArcMap © 10.1.

Figure 8 .
Figure 8. Emerged beach topography investigation.(a) TLS Riegl VZ-400 instrument during a survey campaign; (b) Example of a point cloud during a scan registration with Riscan Pro software: the red point represents the instrument location during the measurement acquisition; (c) Digital elevation model (DEM) of the emerged sector in summer 2018 processed with ArcMap © 10.1.

Figure 9 .
Figure 9. Example of OTS field surveys during winter 2018/19.(a) OTS instrument with the reflective system (prism) for bathymetric point measurement; (b) map showing the bathymetric transect location; (c) bathymetric profiles extrapolated from the measurement carried out in winter 2018/19.

9 .
Example of OTS field surveys during winter 2018/19.(a) OTS instrument with the reflective system (prism) for bathymetric point measurement; (b) map showing the bathymetric transect location; (c) bathymetric profiles extrapolated from the measurement carried out in winter 2018/19.

Figure 10 .
Figure 10.GPR procedure example.(a) 200 MHz antenna during a GPR field survey; (b) map showing the location of the measurement transect; (c) example of a radar stratigraphy section deriving from the measurement path A1-A2 by the use of 400 MHz antenna.

Figure 10 .
Figure 10.GPR procedure example.(a) 200 MHz antenna during a GPR field survey; (b) map showing the location of the measurement transect; (c) example of a radar stratigraphy section deriving from the measurement path A1-A2 by the use of 400 MHz antenna.

Figure 11 .
Figure 11.SBP field survey and data processing.(a) SBP transducer during the data acquisition; (b) hardware system and laptops for data monitoring; (c) map showing the navigation path; (d) 2-dimensional stratigraphic cross sections data processing.

Figure 11 .
Figure 11.SBP field survey and data processing.(a) SBP transducer during the data acquisition; (b) hardware system and laptops for data monitoring; (c) map showing the navigation path; (d) 2-dimensional stratigraphic cross sections data processing.

Figure 12 .
Figure 12.Example of ERT field survey.(a) Graphite electrodes positioning; (b) instrument for sending current into the ground; (c) map showing the ERT transect location; (d) resistivity model.

Figure 12 .
Figure 12.Example of ERT field survey.(a) Graphite electrodes positioning; (b) instrument for sending current into the ground; (c) map showing the ERT transect location; (d) resistivity model.

Figure 13 .
Figure 13.Example of emerged beach sector analysis.(a) Comparison between DTM investigation and Delft3D Model; (b) bioclast analysis and (c) mean size investigation for longshore sediment distribution.

Figure 13 .
Figure 13.Example of emerged beach sector analysis.(a) Comparison between DTM investigation and Delft3D Model; (b) bioclast analysis and (c) mean size investigation for longshore sediment distribution.

Figure 14 .
Figure 14.Gathering the main results for beach dynamics investigation.(a) Petrographical features; (b) sand texture of the bioclastic beach sand; (c) wave height distribution; (d) beach dynamics sketch scheme (modified from [102]).

Figure 14 .
Figure 14.Gathering the main results for beach dynamics investigation.(a) Petrographical features; (b) sand texture of the bioclastic beach sand; (c) wave height distribution; (d) beach dynamics sketch scheme (modified from [102]).

Figure 15 .
Figure15.Example of beach profile with bedrock location deriving from the merging of GPR + ERT + SBP results.The orange dot line in the map depicts the investigation gap.Within the resistivity model and the GPR profile, the distance axis has been exaggerated for a higher result resolution.As ERT and GPR measurements were performed along the same transect of the emerged beach sector, the figures were overlapped to compare the same bedrock location extrapolated from both models.The yellow area between the ground surface/sea bottom and the bedrock represents the beach sand sediments.

Table 1 .
Techniques for beach monitoring applied in the different sectors of the Earth Sciences.

Table 2 .
Methodological proposal for monitoring sandy beaches.Each color represents the methodological steps for beach monitoring.