Assessment of Remote Sensing Techniques Applicability for Beach Morphology Mapping: A Case Study of Hvar Island, Central Adriatic, Croatia

: This study investigates the quality and accuracy of remote sensing data in beach surveys based on three different data sources covering a 10-year period (2011–2021). Orthophotos from State Geodetic Administration Geoportal and satellite imagery from Google Earth were compared with orthophotos generated from UAV using ArcGIS Pro and Drone2Map. The beach area and length of 20 beaches on the island of Hvar were measured using each data source from different years. The average deviation for beach area ( − 2.3 to 5.6%) and length ( − 1 to 2.7%) was determined (without outliers). This study conﬁrms that linear feature measurement is more accurate than polygon-based measurement. Hence, smaller beach areas were associated with higher errors. Furthermore, it was observed that morphological complexity of the beach may also affect the measurement accuracy. This work showed that different remote sensing sources could be used for relatively accurate beach surveys, as there is no statistically signiﬁcant difference between the calculated errors. However, special care should always be addressed to the deﬁnition of errors.


Introduction
Coastal areas, and beaches in particular, are dynamic geomorphological features that are constantly changing. Natural processes, such as waves, tides, rainfalls, relative sea-level changes, combined with anthropogenic activities, lead to changes in the morphological characteristics of beaches. Some of them change rapidly, within a few hours, others slowly, over several decades [1]. Beaches make up 40% of the world's coastline, 70% of which are subject to erosion [1]. More precision is recently provided for sandy beaches revealing that 24% of them are subject to erosion, and 28% are prograding, while 48% are relatively stable [2]. Growing anthropic pressures combined with climate change and sea-level change, e.g., as in References [3][4][5][6][7], makes beaches today highly vulnerable features.
There are many different methods for coastal research, e.g., for coastal monitoring, vulnerability, or risk assessment [8][9][10][11]. However, all research studies need detailed, precise, and high-resolution data. In the past, beach surveys used simple measuring instruments that combined rope, strike, clinometer, and maps. For a time, lower accuracy GNSS was used for mapping and profiling beaches. Today, remote sensing techniques (satellite-and UAV-based) are most commonly used and generally combined with GNSS field measurements. In the last decade, remote sensing data have become readily available and have better spatial and temporal resolution, so the number of papers using a remote sensing methods constantly increase ( Figure 1). This study provides a brief overview of remote sensing techniques in geomorphology, particularly coastal geomorphology, i.e., beach surveys. This study aims to investigate the quality and accuracy of satellite and aerial photo data in beach surveys. This is effectuated through the systematic measurements of the size and length of 20 gravel pocket beaches on Hvar Island from different satellite images and aerial photographs (orthophotos) versus recent UAV-based measurements through the definition of its average errors. Measurements were effectuated on sources covering a 10-year period with the aim to provide an averaged beach size and length database representing the reference point for further decadal analyses of morphological changes of beaches.

Remote Sensing Techniques in Geomorphological Investigations
An important task in geomorphology is to document and analyze landform change. Direct observations can be done for a short period of time and limited areas, while different remote sensing techniques are available for longer time scales. Aerial photography and satellite images today provide an important means of monitoring landforms. This is particularly applicable to features which are in constant change. Likewise, unconsolidated sediments form a beach body that could be subject to frequent geomorphic activity, so significant changes in form may occur.
Remote sensing is applicable to many different types of geomorphological studies: fluvial, e.g., as in References [12,13], karst, e.g., as in References [14,15], glacial, e.g., as in Reference [16], and particularly in coastal geomorphology, e.g., as in References [17,18]. Figure 1 shows a number of Web of sciences core collection (WoSCC) research articles in last 20 years for topics of "Remote sensing" and topic of "Remote sensing" + "coastal geomorphology". The figure shows trends in research on the topics. It can be easily showed that, although the number of articles of "Remote sensing" + "coastal geomorphology" is significantly smaller that "Remote sensing" articles, the trends are quite similar.

Satellite Images
The development of satellite imagery is accompanied by its increase in resolution and quality. Consequently, data from the early days could only be used at small scales, while today's imagery reaches a resolution of less than 0.5 m. Moreover, satellite images have a limitation in small scale measurements considering that they are taken from higher altitudes [19].
Satellite imagery has tremendous advantages, for example, in the study of vast or difficult-to-access areas, such as aeolian landforms [20] or glacial processes [21]. Even though rapid development, increasing global coverage, and many free or relatively cheap data allow easier and more frequent investigations, they could not replace traditional geomorphological fieldwork. Typically, the highest resolution satellite is not precise enough to determine geomorphological changes at small scales, requiring more advanced techniques or fieldwork at local scales [22].
Satellite imagery provides detailed multispectral features and allows advanced geomorphological analyses. The most common research is related to geomorphological mapping, i.e., distribution of geomorphological features and recognition of processes [23,24]. However, other processes, such as anthropogenic impacts that reflect morphological changes, are also often studied. One of such widespread influences relates, e.g., to studies of land cover changes [25][26][27][28][29]. In addition, multispectral features are often used to study hazards, e.g., floods [30], or forest fires, e.g., as in Reference [31]. Furthermore, reforestation, which usually occurs after deagrarization, thus, could have an impact on reducing gullying processes influencing beach erosion [32,33].
Today, thanks to the precise and high spatial resolution of satellite images, digital elevation models (DEM) and digital surface models (DSM) can be extracted in high resolution, which is revealed to be important in geomorphological analyses. Some studies resulted in relatively accurate models from e.g., WorldView 2-3 [34,35].

Unmanned Aerial Vehicle (UAV)
At the beginning of the 21st century, an unmanned aerial vehicle (UAV), also called a Remotely Piloted Aircraft System (RPAS) or drone, became one of the most commonly used methods in geosciences because of its affordability, ease of use, high spatial resolution, and image quality. The UAV is usually used on a local scale. For larger areas, it is better to combine satellite data or available data from other airborne systems, such as aircraft. The acquisition of images by UAVs can be automatic, meaning that the mission is planned in advance by the software, or it can be carried out directly by the pilot. The collected georeferenced images are usually processed using photogrammetric techniques, resulting in orthophotos, DEM, DSM, or point clouds. Low altitude UAV data collection provides very high spatial resolution from 0.01 m to 0.1 m.
UAV provides a fast and efficient survey that is suitable for geomorphological surveys. It can provide a multitemporal comparison for the studies of short-term landform change. Its disadvantage relates to spatial and temporal limitations. Depending on the drone's battery, it lasts about 20-30 min on average; on the other hand, the distance is limited to~2 km on average due to the signal quality. It provides repeated surveys with high-resolution products that can monitor a wide range of geomorphic processes [36]. Consequently, one of the fastest-growing applications of UAVs is right in geomorphology [37]. Compared to satellite imagery, UAV offers even higher quality and resolution, making it applicable at a smaller scale; however, it requires fieldwork. It has been used for geomorphological mapping [38], karst surveys [39], analyses of river flows [40], and many others. Often, authors combine satellite imagery, UAVs, or other remote sensing techniques to obtain more accurate and higher quality data.

Investigations of Beach Morphology
Beach surveys and techniques can be direct or indirect [41]. Direct methods and techniques refer to in situ measurements using GNSS or topographic measurements [32,33,42] and indirect referring to remote sensing data from satellites and aerial photographs [43,44]. Prior to remote sensing capabilities, a standard technique in beach surveys related to terrain profiling and evolved through echo sounders [45]. Today, a huge database of satellite imagery of relatively high quality, covering a longer period of time and a large spatial coverage, has influenced a rapidly growing number of beach surveys ( Figure 1). The availability and simplicity of UAVs have also contributed to the growing number of beach surveys. Table 1 shows some recent papers that applied remote sensing techniques to beach surveys. Using different data sources, they obtain models with very high resolution from 2 m to less than 0.01 m. It can be seen that UAV models provide slightly better resolution (<0.1). Satellite images, especially those from the last century, have low resolution up to 80 m, while recent ones have much better, even high, resolution. Moreover, recent investigations [46] obtained the resolution of 0.5 m from automated extraction of coastlines (WorldView satellite). According to the listed data, remote sensing techniques have proven to be time-saving and useful methods for various beach surveys. Since a beach is a very dynamic form, it responds to changes caused by waves, storms, tides, and to anthropogenic influences. The erosion or deposition of beach material is monitored and studied through remote sensing techniques all around the world. Satellite remote sensing allows long-term studies of beach morphology changes, while UAVs are better suited for short-term changes, such as changes inferred form an intense event. In general, the majority of beach surveys using UAVs focus on multitemporal studies [41].
Satellite and multispectral imagery have been used for various studies related to beaches: their development [47][48][49], morphodynamics [50][51][52], vulnerability [53], or, e.g., beach erosion, on a larger scale [54][55][56][57][58]. In recent years, the use of UAV has become one of the most commonly used methods for monitoring [17,59] or short-term morphological alteration [43,60]. UAV is also combined with archival maps and photographs to detect beach evolution [61]. Data collected by UAV are processed with photogrammetric techniques to obtain DEM, DSM, orthophoto, and point cloud. In addition, the data are also used for multispectral analysis, depending on the specification of the UAV camera.
The number of beach surveys along the Croatian Eastern Adriatic Coast has increased in the last decade. Accordingly, a research methodology is being developed from surveys being conducted by direct in situ measurements in the field, using GNSS receivers, handheld cameras, and usually in combination with other methods (e.g., repeat photography) [32,33,42,62], towards the use of photogrammetric SfM methods by means of a handheld camera that provide creation of high-resolution models [63,64]. Recent studies use UAV for data collection, whose images are also processed with photogrammetric techniques and result in very high-resolution models [61,65]. UAV products should significantly improve coastal and beach surveys in Croatia [10].

Study Area
Hvar is the longest and fourth-largest Croatian island located in Central Dalmatia along the Eastern Adriatic Coast ( Figure 2). It has a particular elongated shape and is oriented in an east-west direction. The island's geology is relatively simple, consisting mostly of limestone and dolomite from the Cretaceous period, Paleogene flysch, and sporadically Quaternary deposits. At the central part of the island, the highest altitude reaches 628 m, and it gradually decreases to the east and west. The research area includes 20 pocket beaches. Most of them are composed of gravel or pebbles. To obtain representative statistical data, we selected beaches along the entire island coast ranging from~150 m 2 to~1800 m 2 (8 on the northern side and 12 on the southern side of the island).
The island has a Mediterranean climate-Csa, with hot and sunny summers and mild, rainy winters. It is the sunniest part of Croatia, often called the sunny island, and also one of the most important touristic destinations in Dalmatia. Thus, the beaches are a very important natural resource for tourism. Today, they are also under increasing anthropogenic pressure. These activities reinforce natural changes, such as erosion, so that beach replenishment is very often required [33,66].
The southern side of the island is affected by sirocco (south-east quadrant) winds that have the greatest impact on the coast, forming larger waves influencing short-term changes in beach morphology. The northern side is dominated by bora (north-east quadrant) wind. The current warm period led to more frequent extreme events, such as increasing precipitation extremes [67]. Such short-term rainfall events may cause important morphological changes of beaches [32].
Tides in the Adriatic are of semidiurnal type. In the Northern Adriatic, their amplitudes are higher than in the rest of the Mediterranean [4,68], which, thus, could influence beach morphology [17,59]. In the Central Adriatic, they are much lower, with an average amplitude of 25 cm [29,[69][70][71], so they do not have an important effect on beach morphology. The influence of average tides on beach area measurements is tested on Lučišće beach ( Figure 2, beach no. 6), a beach without anthropogenic influence. The beach was surveyed with a GNSS receiver and recorded with a UAV on high and low tide on the same day. The collected data were processed and analyzed. First estimation of area difference between low and high tide was around 6%.

Data Acquisition
To investigate the accuracy of aerial orthophotos and satellite imagery as a base for further studies of decadal changes of beach morphology, we effectuated repeated measurements of beach area and length of 20 defined beaches from different sources and years.
For data collection, we used UAV: DJI Phantom 4 Pro v.2.0 (Da-Jiang Innovations, Shenzhen, China). The drone has an RGB camera FC6310 with a resolution of 20 megapixels. The mission was not planned but was coordinated by the pilot on-site. The average flight altitude was 20-30 m, and the overlap between images was 70-80%. Flight duration averaged 15 min. During the fieldwork, about 3500 images were collected at 20 different study locations. We also used Trimble GNSS GeoxH 6000 (Trimble Navigation Limited, Sunnyvale, CA, USA) to collect ground control points (GCP), 3-8 at each beach. In the moment of collecting GCP's positions (x,y,z), GNSS receiver was connected to Croatian Positioning System (CROPOS) which obtained higher horizontal and vertical accuracy (<0,1 m). Surveying was done in accordance with Croatian laws, including registration of drones, reservation of the portal AMC-Airspace Management Cell and permits from the State Geodetic Administration. The fieldwork took place from 5 to 11 November 2020 (first part), and from 16 to 22 May 2021 (second part).
The orthophotos were created using the software Drone2Map, which is based on structure-motion algorithms (SfM). SfM is an automatic photogrammetric technique for generating orthophotos, DEM, DSM, and point clouds from overlapping images [72,73]. Even though drones provide GNSS metadata to the images, ground control points (GCP, collected by Trimble GNSS receiver) were used for better accuracy of the model. Drone2Map combines imagery from UAVs and imported GCPs that are manually merged with the imagery. Tie-points were automatically determined on all images using SfM algorithm. Photogrammetric data processing resulted in high-resolution models.
This study used historical satellite imagery from the desktop application Google Earth Pro (GE) (Google LLC, Mountain View, CA, USA). The oldest imagery covering the island of Hvar is from 1985; however, since this was almost 40 years ago, the resolution was very low, so it could not be used.
Furthermore, due to the elongated shape of the island, satellite imagery had not always covered the entire island at the same time (e.g., the western and eastern part of the island was not recoded by imagery in the same year). Consequently, we analyzed all available satellite images at Google Earth Pro, which covered the island in the same year and had acceptable resolution (Table 2). Thus, we selected 5 generations of satellite imagery: 2013, 2016, 2018, 2019, and 2020. We also used orthophotos provided by the State Geodetic Administration from Geoportal (further in the text-SGA). Orthophotos have been produced from aerial photogra-

Measurements and Accuracy Assessment
This study measured beach areas and lengths from three different data sources. Beach is defined as a non-vegetated sediment body (Figures 3 and 4). If vegetation was too high and masked the view, the beach boundary was approximated using the method of generalization. Some of the beaches were modified by anthropogenic activities, so the boundary was determined following constructions (walls or buildings). All beach measurements were done in the same way, manually in the scale 1:500-1:1000. We measured each beach area and lengths min 3 times and then calculated the average values.
From the orthophotos provided by the State Geodetic Administration (SGA), beaches were measured using the Area measurement tool in the Geoportal Web and connecting the Web Mapping Service with the software GIS (https://geoportal.dgu.hr/services/dof/wms, accessed on 15 October 2021). No important differences between the two have been observed; consequently, they are not further discussed.
From satellite images, beaches were measured in the application Google Earth Pro using the ruler tool-polygon-which has been shown to be relatively accurate for measurements at scales smaller than 1:30,000 [74]. In addition, finally, orthophotos generated by a UAV were analyzed in GIS using ArcGIS Pro 2.7.3. software (ESRI Company, Redlands, CA, USA). The beach areas and lengths were done using the area measurement tool (Figure 3). Two sites (beaches Grabovac and Zečja) were not measured due to low image resolution (high vegetation or shade).
The average values of beach areas and lengths were calculated for each beach, based on all available sources for each year, and separately according to the sources (Geoportal and satellite images).
All measured values were examined with Shapiro-Wilk test to check if the beach areas and lengths values are normally distributed. The accuracy of the values was calculated using several statistical methods: ∆A i -area error, which is used for RMSE calculation; and RMSE (Root Mean Square Error), which presents standard deviation of the residuals and % ERROR -percentage error.
A SGA/GE = beach area measured on Geoportal/Google Earth, A U AV = beach area measured on U AV.
UAV, combined with GNSS, has been shown in many works to be an efficient tool for high-resolution results [17,41,75]. The spatial resolution of UAV orthophotos is the most accurate and has ten times better resolution then SGA or GE (Table 2), so it was chosen as a reference value. All measured values from GE or SGA were compared to UAV. Finally, we used Kruskal-Wallis non-parametric method to test the statistical significance ( Figure 5).

Results
Measurements from three different data sources (Geoportal (SGA), Google Earth (GE), and UAV) from different years yielded 198 beach area and 200 beach length values for twenty locations, which are shown in the Figure 6.

Beach Area
The average beach areas, obtained from Geoportal (SGA) and Google Earth (GE), are less than the UAV, with the exception of 2021GE (Table 3). In general, the beach areas measured on the Geoportal orthophoto have slightly smaller values than those measured from Google Earth (22.7 m 2 and 5.7 m 2 less than the UAV, respectively). The Geoportal (SGA) and Google Earth (GE) strongly correlate with UAV values, (r 2 ) ranged from 0.981 to 0.997.
In addition, the average values from Geoportal and Google Earth were compared separately with the UAV values, and a very small difference was found (SGA r 2 = 0.971; GE r 2 = 0.991).
The average values of the beach area are presented chronologically (Figure 7). It can be seen that the deviation from the reference value also changes chronologically, as shown by the trend line. In the first part of the decade, the deviations are larger, while, in the recent period, they are lower. In addition, comparing different sources (satellite-orthophoto), it was observed that values of satellite images (GE) are slightly closer to the average value than the SGA orthophoto.
The mean area of each beach measured from different sources overlap (Figure 8). Data from Table 3 and Figure 8 show that smaller beaches revealed slightly larger differences.   As the Shapiro-Wilk test showed, the deviation from the normal distribution (p = 0.03) nonparametric tests were further chosen.
The accuracy error was tested first using the RMSE ( Table 4). The GE2021 satellite imagery revealed the lowest error, while the SGA2011 orthophoto showed the highest. In general, Geoportal has an average error 28.08 m 2 higher than Google Earth. When the RMSE is ordered chronologically, the trend points downwards, which means a reduction in error over the years (approaching the present) (the newer the source the smaller the error). In addition, as can be seen in Figure 8, there is a difference between the measured areas quality between small and large beaches. The calculated average RMSE for smaller beaches (with an area of less than 650 m 2 ) is about 20 m 2 higher than the RMSE for larger beaches (Table 4). The measurement error was also tested with the percentage error. It ranges from −4.1 to 8.1% for all measurements, and, when the outliers are removed, the result changes from −2.3% to 5.6%. The boxplot diagram (Figure 9) shows slightly larger errors for the Geoportal orthophoto measurements with errors from −6.1% to 8.5%. The 2017 measurements SGA have the largest error from −9.3% to 11.9%. The satellite image measurements show much better results. The average error ranges from −3.8% to 7.9%. The 2021 images provide the lowest error, −2.7% to 6.1%, with respect to the reference value. The calculated percentage errors were tested with the non-parametric Kruskal-Wallis test, which showed that there is no statistically significant difference (p = 0.572) between the calculated errors from the listed sources.

Beach Length
The average values of beach length show relatively similar results between the sources ( Table 5). The smallest length value was measured on Jagodna beach, 17.47 m, while the longest relates to Lučišće (Brusje) beach, 127.41 m. Their average length value (UAV) is 65.24 m. SGA and GE average measurements have very similar values; they differ by only 0.24 m. The lowest value is measured on the SGA2017, while the highest value relates to SGA2014. The average value of Geoportal differs by 0.89 m from the UAV, while the value of Google Earth differs by 0.65 m. Each source was compared to UAV, the reference value which revealed very strong correlation (r 2 ≥ 0.997). All measured values overlap ( Figure 10). Only some larger beaches show very small differences. These differences are tested with the RMSE ( Table 6). The lowest RMSE was calculated for the SGA2014, and the highest error for the SGA2011. The average error is slightly lower for the Google Earth values than for the Geoportal values, with a very small difference of 0.1 m. However, difference in RMSE was observed according to beach size (length), 2.03, for large beaches, and 1.24 for small beaches. The percentage error showed that the average deviation of SGA and GE from the UAV value was −2.35 to 2.18. When the outliers are removed, the deviation is smaller, from −1.01 to 2.68%. The boxplot diagram ( Figure 11) shows that the average deviation of beach lengths obtained from Google Earth measurements is smaller than those from Geoportal. However, the difference is very small (< ±1%). The Kruskal-Wallis test showed that there is no statistically significant deviation (p = 0.322).

Discussion
This work has demonstrated that there are no significant statistical differences between values obtained from different sources. All measured beach areas and lengths have showed relatively similar values. However, there is a difference in the error range. In this study, the models generated by UAV resulted in measurements of very high resolution (0.02-0.05 m), while satellite imagery and available orthophotos proved to be slightly less precise (compared to UAV), with an average deviation of −4 to 8% for beach area, and −2 to 2% for beach length. Excluding the outliers, the deviation reaches −2.3 to 5.6% for beach area, and −1.01 to 2.68% for beach length, while error (RMSE) reaches 7.2% for beach area, and 2.5% for beach length. Both, percentage error and RMSE, showed similar range of deviation. As the measurements cover in total a 10-year period, natural and anthropogenic processes [40] may also affect the obtained results.
Depending on the size of the beach, differences in measurement errors were observed. Larger beaches revealed lower errors in beach areas but higher variation in beach length measurements, while it is just the opposite in the case of small beaches. This may be related to the beach boundary definition, which is often more challenging in the case of small beaches, as approximations revealed to be more often. Even if digitizing errors are low in the case of small beaches, they represent a high percentage error.
Beach length measurement errors were revealed to be very low. In the case of length measurement, large beaches have somewhat larger errors. This may be related to more complex morphology (more curved shoreline) ( Figure 12) compared to small, simple beaches.
Trend analysis of the calculated error over the 10-year period showed that the error value decreases chronologically, which may be related to the quality of satellite imagery and resolution increase through the years.
Most studies that use aerial images in geomorphology focus on linear features (e.g., shoreline change, changes of watercourses), while only a small number of them examine features as a whole. Solazzo et al. [75] studied aeolian forms and tested the accuracy of dune measurements from two different aerial images. They selected 4 different dunes with large surface areas (4000-8000 m 2 ) and calculated a deviation of 4%, which is quite similar to the deviation of beach area measurements. Swanson et al. [76] calculated uncertainty in aerial photography measurement for river width and areas. They used digitized and georeferenced aerial images from 1985 to 2008 and found out that river width error was 2-6% for the recent data, while those from 1985 had a higher error. River area measurement error was 3-12%. The error decreases chronologically, so recent photos are most accurate. In addition, measurements on a smaller polygon, e.g., river islands, resulted in higher error then those on the larger areas. Those results could be related to Hvar Island beaches, where the calculated error is higher on small beach areas and also decreases chronologically.
Apart from geomorphological landform measurement, accuracy could also be tested using, e.g., cadastral data. Lopes and Nogueira [74] compared Google Earth satellite imagery in terms of point, line, and polygon measurements with the accurate cadastral data. They found that the values measured on linear features differ for 0.44%, while polygon measurements differ for 3.54%. The calculated errors are close to errors obtained here in the case of beach area measurement on the island of Hvar. In general, higher errors are to be expected when measuring geomorphological features than other objects due to the problem of boundary definition; consequently, some approximation is required.
Many circumstances can affect the accuracy of measurement. The accuracy of beach length measurement is very important because a low error could lead to significant discrepancies in the final results [56]. Therefore, to decrease the measurement error to a minimum, the number of measurements can be increased. Here, beach area and length measurements were done at least 3 times for each source and year. The complexity of beach morphology could affect the detection of the coastline in orthophotos [43]. For sandy beaches, the definition of coastline is usually easier, due to the simple recognition of the wet-dry boundary [17], which leads to higher measurement accuracy. Complex beaches with different sediment sizes can affect the lower measurement accuracy. In addition, vegetation, shade, or even deposited banks of seagrass species (Posidonia) can affect the difficulty of shoreline detection. In this case, two beaches from different sources could not be measured in our study.

Conclusions
This study investigated the accuracy of remote sensing techniques in beach surveys. Several available remote sensing sources were found to be sufficiently accurate for beach measurements. A measured beach area and length values from the SGA and GE sources correlate very strongly with the UAV (r 2 ≥ 0.97). Their calculated error for beach area (7.2%) and length (2.5%) was found to be in a different range. Although most of the beaches are small, according to UAV measurements, with the surface area ranging from 175.9 to 1828.0 m 2 and a length from 17.47 to 127.41 m, the calculated error is higher in the case of beaches of smaller surface areas and in the case of beaches of larger lengths.
As it is not possible to accurately distinguish the short term morphological changes from most of the data sources available, particularly from older ones, they are part of the obtained errors. The acquired geodatabase on Hvar Island beaches is crucial for further investigations of beach evolution and determination of beach vulnerability, as well as for future beach management, as Hvar Island is one of the most important touristic destination along the eastern Adriatic coast.
The results have shown that different remote sensing sources could be used for an accurate geomorphological survey of beaches, considering the calculated error. Thus, remote sensing data, which is free and readily available, may allow a spatial and temporal study of beach morphology over decades. In addition, it could be globally applied for similar investigations, not only in coastal geomorphology, but also in other scientific disciplines.