The Origin of Sand and Its Colour on the South-Eastern Coast of Spain: Implications for Erosion Management

: Sand colour can give important information about mineral composition and, consequently, sediment source areas and input systems. Beach appearance, which is mostly linked to sand colour, has a relevant economic function in tourist areas. In this paper, the colour of 66 sand samples, collected along both natural and nourished beaches in the western Mediterranean coast of Spain, were assessed in CIEL*a*b* 1976 colour space. The obtained results showed relevant di ﬀ erences between natural and artiﬁcially nourished beaches. The colour of many nourished beaches generally di ﬀ ers from the native one because the origin of the injected sand is di ﬀ erent. The native sand colour coordinates’ range is: L* (40.16–63.71); a* ( − 1.47–6.40); b* (7.48–18.06). On the contrary, for nourished beaches’ the colour range is: L* (47.66–70.75); a*(0.72–5.16); b* (5.82–18.82). Impacts of beach nourishment on the native sand colour were studied at San Juan beach, the most popular one along the study area. Nourishment works were performed after severe erosion, usually linked to anthropic activities / structures and storm events, but also to increase beach width and hence beneﬁt tourism.


Introduction
Coastal sedimentological studies allow the understanding of both natural processes and anthropogenic interventions, e.g., coastal dynamics and required management actions. Beach sediments characteristics play an essential role in coastal ecology but also in the choice of a pleasant site by beach users [1]. The most important characteristics of the sediment are texture (grain size and sorting), mineralogy, colour and morphology of grains (roundness and sphericity). These characteristics should be preserved during anthropogenic actions since their modifications can cause undesired impacts. First, sediment colour and texture directly influence the perception of beach users [2] and coastal scenery characteristics [3]: beach sediment colour is one out of the 18 natural parameters considered in the probably most commonly used landscape assessment method, i.e., the Coastal Scenic Evaluation [4]. This methodology, applied in >1000 coastal sites around the world [5], including the study area of this paper [6], is based on enquires that revealed that golden or white sand is better valued than dark sand (natural or nourished), to nourish beaches attempting to employ natural-like sand colour, and even to detect future changes in the sand colour or the presence of contaminants.

Study Area
The studied area covers about 200 kilometres of coastline of the western Mediterranean (SE Spain) ( Figure 1). It includes 60 beaches from San Pedro del Pinatar (Murcia Province) to Dénia (Alicante Province, Figure 1, Table 1). The studied coast has very diverse habitats: dunes in the south [30], high cliffs of the Betic Cordillera in the north coast of Alicante [31], as well as beaches located on the Nueva Tabarca Island [23]. From a geological standpoint, this coastline is divided into three different domains: the southern sector mainly belongs to the Bajo Segura basin [32,33], the northern one pertains to the Prebetic zone of the External Zones of the Betic Range [34], while the Island of Tabarca is related to the Internal Zones of the Betic Range [35]. The different geological contexts in which the beaches are included establish broadly three different source areas of sediments, which constitute an intrinsic factor that determines natural sand composition and hence the native beach colour.

Study Area
The studied area covers about 200 kilometres of coastline of the western Mediterranean (SE Spain) (Figure 1). It includes 60 beaches from San Pedro del Pinatar (Murcia Province) to Dénia (Alicante Province, Figure 1, Table 1). The studied coast has very diverse habitats: dunes in the south [30], high cliffs of the Betic Cordillera in the north coast of Alicante [31], as well as beaches located on the Nueva Tabarca Island [23]. From a geological standpoint, this coastline is divided into three different domains: the southern sector mainly belongs to the Bajo Segura basin [32,33], the northern one pertains to the Prebetic zone of the External Zones of the Betic Range [34], while the Island of Tabarca is related to the Internal Zones of the Betic Range [35]. The different geological contexts in which the beaches are included establish broadly three different source areas of sediments, which constitute an intrinsic factor that determines natural sand composition and hence the native beach colour.  Table 1. Currents in Alicante area, observation period: 2006-2020 (source: www.puertos.es, accessed in January 2020).
Different natural elements reinforce the tourist attraction of the studied coastal scenery, such as cliffs, rocky shores, golden sand, dunes, turquoise water colour, and other coastal landscape features such as islands, lagoons, caves, rock ridges, and irregular headlands. However, most of the studied sites have been affected by infrastructures such as groins, seawalls, houses/buildings, parking, and nourishment works among others. Such human actuations/structures are mostly linked to the tourism demand and strong urban growth in the area and they have been especially frequent during recent years [36]. These anthropic structures were recently observed [6] and took into consideration 26 parameters-18 physical and 8 human-that delineate coastal scenery, according to Coastal Scenic Evaluation (CSE) Methodology [4,37], applied in different countries  Table 1. Currents in Alicante area, observation period: 2006-2020 (source: www.puertos.es, accessed in January 2020). Different natural elements reinforce the tourist attraction of the studied coastal scenery, such as cliffs, rocky shores, golden sand, dunes, turquoise water colour, and other coastal landscape features such as islands, lagoons, caves, rock ridges, and irregular headlands. However, most of the studied sites have been affected by infrastructures such as groins, seawalls, houses/buildings, parking, and nourishment works among others. Such human actuations/structures are mostly linked to the tourism Water 2020, 12, 377 4 of 25 demand and strong urban growth in the area and they have been especially frequent during recent years [36]. These anthropic structures were recently observed [6] and took into consideration 26 parameters-18 physical and 8 human-that delineate coastal scenery, according to Coastal Scenic Evaluation (CSE) Methodology [4,37], applied in different countries [5,8,38,39]. This semi-objective method utilizes weighting parameters and fuzzy logic mathematics to obtain the scenery Decision Value parameter (D). This value allows the coastal areas to be classified into 5 classes, from Class I (extremely attractive/natural sites) to Class V (very unattractive/urbanized sites). Table 1. Location and main characteristics of studied beaches: Map number (Figure 1), beach name, municipality, beach typology and sand condition with their information source.

Sampling Design
A total of 66 sand samples were gathered in [16][17][18][19] February 2019, at sixty Mediterranean beaches (Table 1). Each sample was collected in the middle of the backshore (i.e., in the dry beach) and contained approximately 300 g of sand. First, the surface sand was removed and then, the sample was gathered in the first 20 cm of depth. On the longest beaches such as "La Mata", "San Juan", "Carrer la Mar", and "Ponent", several samples (2-3) were collected (Table 1). In this paper, only sandy beaches were investigated. For this reason, the north-eastern region has some coastal sectors not sampled due to the presence of cliffs or beaches mainly composed of gravel and boulders ( Figure 1). Complementary information about the existence of anthropic actions on the sampled beaches was collected from both press news and coastal management actions by involved municipalities (Table 1).

Sand Colour Measurement and Analysis
The colour of sand samples was determined in the laboratory of the University of Florence (Italy) using the CIEL*a*b* 1976 colour space [25]. It required about 100 g of sand per sample to measure the colour. Moreover, mineralogical composition was quantified under a binocular microscope (Nikon SMZ 1500, Nikon Corp., Japan) from the Earth Sciences laboratory of the Alicante University (Spain). For this analysis a representative sample (about 200 sand grains) was analysed and classified taking into account three main groups of components: quartz, (i.e., leucocratic component), rock fragments, and shell fragments (i.e., melanocratic components), then the values were recalculated to 100% in order to synthesize and examine the results.
Principal Component Analysis (PCA) was carried out in order to establish multivariate relationships between samples as well as connections among them considering the CIEL*a*b* coordinates as input data.
Equations (1)-(3) were applied to quantify the chromatic change in a representative nourished beach (sample No. 35, Table 1). For this, a native sample was required and measured in the CIEL*a*b* colour space with the rest of the 66 samples.
Chroma (C*) = (a* 2 + b* 2 ) 1/2 (2) The distance between colours is given by the Euclidean distance: ∆E*ab = [(∆L*) 2 + (∆a*) 2 + (∆b*) 2 ] 1/2 (3) Sedimentological analysis (particle size distribution) was performed via dry sieving at 1 phi interval in order to have a more detailed database. For this characterization, initially, each sample contained 200 g of sand. The software used for data analysis was "Gradistat" [53] and granulometric parameters were calculated according to Folk and Ward (1957) [54]. Table 1 shows the main beaches under study, and their typology according to the Bathing Area Registration and Evaluation (BARE) classification system [55]. This classification is based on the level of urbanization, accessibility, etc., and establishes four basic types: remote, rural, village, and urban. Remote areas are mainly defined by difficulty of access. They may be adjacent to rural areas and, occasionally, village environments. Usually, remote beaches are not supported by public transport and have very limited (0-5) or no temporary summer housing [55]. Rural areas are located outside the urbanized environment. They are not accessible by public transport and present virtually no facilities. Housing is limited in number, generally ranging from 0 to 10. Rural beaches do not have permanent community focal centres (e.g., schools, shops, cafes and bars) and they have little or no beachfront development [55]. Village areas are located outside the focal urban environment and related to a small and permanent population. Village beaches may be reached by public or private transport [55]. Urban areas have large populations with well-established public services and a wellmarked central business district. In the proximity of urban areas can be found commercial activities such as fishing/boating harbours and marinas. Urban beaches are located within or adjacent to the urban area and are generally freely open to the public [55].

Typology of Beaches
The selected study case is an excellent example of coastal development under high touristic pressure, with many coastal sites showing very urbanized conditions (34 studied beaches classified as Urban in Table 2) and the rest of the cases are located in urbanized environments (13 Village beaches) or they show the presence of houses (14 Rural beaches). Only a very few coastal sites (5) can be classified as natural beaches, i.e., beaches with few anthropic interventions. They were classified as Remote. As a consequence of the urban and touristic pressure that affects the study area, some beaches have been modified by anthropogenic activities, but others have preserved their original state, i.e., only receive sand through natural processes. For this reason, a new classification was made according to sand condition: "Nourished" and "Native" (Table 1).
Nourished beaches are common in the study area [41, [49][50][51][56][57][58], causing anthropogenic modifications in the original appearance of the coast [59]. Construction of dams on land that retain sediments, massive construction of buildings near the coast, creating new harbours and breakwaters, beach nourishment, as well as sand movements from beach to beach are some examples. These coastal modifications that have been carried out since approximately the middle of the 20th century [60], have changed the native sand colour among other sediment characteristics. Figures A1-A5 show the current state (February 2019) of the sampled sites. Sand colour observed in these images is a result of multiple factors acting simultaneously: the camera used, hour, luminosity of the moment (e.g., sunny or shaded), beach surface characteristics, etc. All these extrinsic determinants modify the perception of the real intrinsic sand colour, which is mainly a consequence of the mineralogical composition. For these reasons it is necessary to apply an objective methodology to quantify the chromatic appearance of samples, i.e., in a laboratory with a quasi-uniform colour space method (CIEL*a*b*).

Sand Colour Analysis
An interesting diversity of colours emerged along the study area ( Figure 2, Table A1) mainly due to two factors: the different origin of sediments (intrinsic factor) and the anthropic modifications that have been carried out since approximately the middle of the 20th century (extrinsic factor).

Sand Colour Analysis
An interesting diversity of colours emerged along the study area ( Figure 2, Table A1) mainly due to two factors: the different origin of sediments (intrinsic factor) and the anthropic modifications that have been carried out since approximately the middle of the 20th century (extrinsic factor). Regarding the colour parameters (the L*, a* and b* values) expressed in Figure 3a,b, the following considerations can be argued: These results highlight the relevant underlying influence of the geological context on the final appearance of the coastal deposit. This geological frame defines the natural source area of beach grains. In this sense, the above first group includes the southern beaches. Their geological context is mainly related to the Bajo Segura basin. This source area is mostly constituted of Neogene-Quaternary sediments mainly related to continental and coastal deposits [61]. Also, this area is characterized by unaltered coastal stretches. The northern beaches related to the External Zones of the Betic range are grouped in the second group. They are the most modified beaches and only a few of them are natural (Figure 3a). The dark colour of the northernmost beaches (included in the third group) is related to the contributions of different rivers and streams (e.g., Girona River and "Barranc de L'Alberca"). However, the darkest sands correspond to those included in the fourth group. These samples are included in the geological context of the Tabarca Island (Figure 3a). In this island the outcrops are mainly materials from Internal Zones of the Betic range. The low Lightness of these natural beaches is justified by the high content of dark rock fragments that comes from the black dolomites and diabases, which outcrop along the island [23,52].
Two samples differ significantly from the general trend defined above. On the one hand, sample No. 62 differs from the rest of the samples located in Tabarca Island due to the influence of the yellowish calcarenites (Miocene age) that outcrop at this point of the island [23,34]. On the other hand, sample number 51 shows the highest Lightness value (i.e., 70.75, Table 3, Figure 3a) because this sand is from an inland quarry.

The a* and *b Plane
The parameters a* and b* (Table A1) are positively correlated among them (R 2 = 0.89; p-value <0.05; Figure 3b and Table A2). The found values were more uniform than the Lightness ones. On the one hand, the parameter a* ranges from −1. Neogene-Quaternary sediments mainly related to continental and coastal deposits [61]. Also, this area is characterized by unaltered coastal stretches. The northern beaches related to the External Zones of the Betic range are grouped in the second group. They are the most modified beaches and only a few of them are natural (Figure 3a). The dark colour of the northernmost beaches (included in the third group) is related to the contributions of different rivers and streams (e.g., Girona River and "Barranc de L'Alberca"). However, the darkest sands correspond to those included in the fourth group. These samples are included in the geological context of the Tabarca Island ( Figure  3a). In this island the outcrops are mainly materials from Internal Zones of the Betic range. The low Lightness of these natural beaches is justified by the high content of dark rock fragments that comes from the black dolomites and diabases, which outcrop along the island [23,52].
Two samples differ significantly from the general trend defined above. On the one hand, sample No. 62 differs from the rest of the samples located in Tabarca Island due to the influence of the yellowish calcarenites (Miocene age) that outcrop at this point of the island [23,34]. On the other hand, sample number 51 shows the highest Lightness value (i.e., 70.75, Table 3, Figure 3a) because this sand is from an inland quarry.

The a* and *b Plane
The parameters a* and b* (Table A1) are positively correlated among them (R 2 = 0.89; p-value <0.05; Figure 3b and Table A2). The found values were more uniform than the Lightness ones. On the one hand, the parameter a* ranges from −1.    Table A1). Colour (yellow, green and blue) refers to the geological domains of the study area (respectively, The Bajo Segura basin, The Betic External Zones, and the Tabarca Island which belongs to the Betic Internal Zones).

Mineralogical Characterization
Mineralogical composition of the studied sands is very diverse, showing variations in the type of minerals and also in the relative proportion in which they are represented. The colour of each grain is directly related to its mineralogical composition and, therefore, the global colour of the deposit is the result of the presence of certain grains and their abundance.
Leucocratic components (mainly transparent grains with angular morphology) were identified, and they correspond mainly to quartz (SiO 2 ). They are very well represented in all samples in variable proportions (from 8.08% to 72.26%, Table 3). Melanocratic components correspond mainly to both rock fragments (basically calcarenites and dark carbonate rocks) and other minerals such as pyroxenes, amphiboles, and biotites. In general terms, melanocratic components are more abundant than leucocratic ones, especially at Tabarca Island ( Figure 4). Their abundance varies along the study area from 27.74% to 91.92% (Table 3). Finally, bioclasts (i.e., small fragments of shells, sea urchins, calcareous algae, corals, etc.) were also counted and their proportion in the study area is very low, almost insignificant. Due to the dark aspect of bioclasts and their low abundance, they were included in the melanocratic fraction.
Concerning the abundance of leucocratic and melanocratic components, 14 samples contained less than 50% of leucocratic components, while the rest (52) contained ≥50% of melanocratic components (Figure 4). Table 3 shows median, mean, as well as minimum and maximum values of the colour data (CIEL*a*b*) and mineralogical composition, i.e., percentage of leucocratic and melanocratic components.

Statistical Analysis
In order to find correlations among the different colour and mineralogical parameters, the Pearson correlation statistic coefficient was obtained (Table A2). On the one hand, the highest correlation was found among a* and b* parameters (R 2 = 0.89; p-value = 2.2 × 10 −16 ). On the other hand, L* parameter shows a good correlation with the leucocratic component (R 2 = 0.62; p-value = 2.187 × 10 −8 ). Quartz grains normally present light colours and their high proportion in samples determines sand clarity. In this way, the content of quartz grains is inversely proportional to the a* and b* colour parameters as was previously concluded for the Belgian coast [9].
Principal Component Analysis (PCA) is a multivariate-statistical technique used for reducing the dimensionality of data. The resulting analysis is shown in a two-dimensional plot composed of

Statistical Analysis
In order to find correlations among the different colour and mineralogical parameters, the Pearson correlation statistic coefficient was obtained (Table A2). On the one hand, the highest correlation was found among a* and b* parameters (R 2 = 0.89; p-value = 2.2 × 10 −16 ). On the other hand, L* parameter shows a good correlation with the leucocratic component (R 2 = 0.62; p-value = 2.187 × 10 −8 ).
Quartz grains normally present light colours and their high proportion in samples determines sand clarity. In this way, the content of quartz grains is inversely proportional to the a* and b* colour parameters as was previously concluded for the Belgian coast [9].
Principal Component Analysis (PCA) is a multivariate-statistical technique used for reducing the dimensionality of data. The resulting analysis is shown in a two-dimensional plot composed of two axes or components, PC1 and PC2 ( Figure 5). PC1 explained 63.57% of the total variance and PC2 32.83%. The cumulative proportion of the second principal component was 96.39% (the data were previously normalized). In the graph, the three colour variables (L*, a*, and b*) are represented with arrows and each one of the 66 investigated sites is located in the graph according to its relation with the variables, a technique applied previously in another sand colour study with important implications for beach nourishment and coastal management [26]. The resulting distribution of the samples in the graph is due to both their geographical location mainly related to the geological context and anthropic actions carried out on the deposit (i.e., nourishment works). Four groups of beaches were identified in the graph of Figure 5 according to their geographical location ("neighbouring beaches"): (i) the southernmost beaches (n = 7); (ii) the northernmost beaches (n = 4); (iii) beaches of Santa Pola (n = 3); and (iv) beaches of Tabarca Island (n = 5).
Moreover, very interesting groupings can be observed according to sand condition consulted in the bibliography (Table 1). A group of natural samples (+10) from different municipalities is shown towards the arrows of the colour variables (a* and b*, Figure 5). On the left side of the graph ( Figure 5) there are generally those beaches identified as natural or without alterations in colour of sand. The sand composition of the northern beaches was also considered as native (Table 1) where sediment supplies mainly arrive from the river and nearby streams. The sand of Tabarca beaches is also native and its source area is the island itself and seabed. In the central sector of Figure 5 there are beaches located in the central-northern part of the study area. Some of them correspond to modified beaches by nourishment works and emplacement structures, e.g., harbours (No. 38,39,41,46,48,53,57).
Water 2020, 12, 377 11 of 26 ( Figure 5) there are generally those beaches identified as natural or without alterations in colour of sand. The sand composition of the northern beaches was also considered as native (Table 1) where sediment supplies mainly arrive from the river and nearby streams. The sand of Tabarca beaches is also native and its source area is the island itself and seabed. In the central sector of Figure 5 there are beaches located in the central-northern part of the study area. Some of them correspond to modified beaches by nourishment works and emplacement structures, e.g., harbours (No. 38, 39, 41, 46, 48, 53, 57).  Table A1 for numerical information of each variable. Green colour: native samples; red colour: altered samples (information obtained from bibliography data, Table 1).

Grain Size Parameters
The sediment studied from the 66 samples was classified mainly as coarse sand (78.8%), very coarse sand (12.6%), and medium sand (6.5%). Beaches in the southernmost part of the study area  Table A1 for numerical information of each variable. Green colour: native samples; red colour: altered samples (information obtained from bibliography data, Table 1).

Grain Size Parameters
The sediment studied from the 66 samples was classified mainly as coarse sand (78.8%), very coarse sand (12.6%), and medium sand (6.5%). Beaches in the southernmost part of the study area were composed essentially of fine sand, while from the central part to the north, different types of beaches can be observed with rocky shores, gravel, pebbles, with fine, medium and coarse sand. Table A1 shows some of the grain size parameters of each sample.

Colour Modification due to Human Activities
As a consequence of the urban and touristic pressure, artificial beaches were created and deep modifications in the native coastal deposits were carried out. The sand used to nourish the beaches came from different sources: (i) submerged sandbanks (e.g., San Juan beach, No. [34][35], sometimes near to the mouth of harbours; (ii) from other beaches (e.g., El Puerto beach, No. 4); and even (iii) from a quarry inland (Fustera Beach, No. 51).
At the beginning of the 1990s one of the largest beach modifications in the province of Alicante was carried out in San Juan beach ( Figure 6). Thanks to a sample of sand from the late 1980s it was possible to compare the native colour of the sand with the current one of the nourished beach, whose sand grains come from the littoral bottoms in front of the Serra Gelada Mountain (Benidorm, Alicante) [57]. Table 4 shows the colorimetric parameters of two samples from different years (1989/2019) and gathered in the same coastal sector of the beach. Difference of colour values between these two samples was estimated [62]. Equations (1) and (2) show the mathematical expression for quantifying the chromatic change in terms of Hue differences (∆h) and Chroma differences (∆C*). The parameter Total Colour Difference (∆E*, Equation (3)) is also used. The above calculations show that the colour of the beach after nourishment works can change even when sediments come from near submerged areas. The colour change is registered in all the chromatic parameters L*, a*, and b* (Table 4), and differences of −0.4 • (h), −6 (C*), and 7.15 (∆E*) are evidenced by different equations. All these values reflect a significant aesthetical change of the coastal deposit after the nourishment. On the one hand, the negative values of ∆h and ∆C* mean that the new appearance of the altered San Juan beach is duller and more brown than the native aspect. Therefore, the new beach colour is more distant from the ideal golden or light sand appreciated by beach users [6]. On the other hand, the high ∆E* value indicates that the chromatic change is clearly perceptible by local beach users, this occurs when the Total Colour Difference parameter exceeds three units (∆E* > 3) [63]. Another similar example is the Varadero beach (Cuba) where beach nourishment works performed several times in the past decades with nearshore sands changed the original beach colour that was the second lightest one out of 100 in Cuba, now occupying only the fiftieth position in rank [26].
San Juan beach is a clear example showing how the construction on the coast imposes a barrier to wind dynamics that can even cause the disappearance of beaches with the consequent requirement of nourishment works ( Figure 6). In addition, river channelling and sediment retention in upstream reservoirs, which contribute to the erosion of beaches [64,65], have played a significant role in the reduction of the sediment input to this representative beach of south-eastern Spain. Historical bathymetric measures in reservoirs of the study area confirm this issue (Centre for Hydrographic Studies). Indeed, San Juan beach has been the subject of geological outreach activities to promote the Geoheritage of the area, which has also very didactic geological and sedimentary outcrops [66]. In the didactic activities undertaken [67], the sedimentological analyses of the sand (colour, sorting and presence of bioclasts) has been used to understand and reflect on the differences between the natural and the nourishment sediments of the beach, and how important it is to attempt to archive equivalent sands to the original in beach nourishment works. The study case of "Cala Fustera" (No. 51, Table A1, Figures 7 and A4) constitutes another example of significant coastal alteration after nourishment works. In this case, the previous sand deposit was replaced in 2017 by artificial quarry sand. Figure 7 shows visually differences in the colour of the coast of the Fustera Beach (now composed of quarry sand) with respect to the native colour in this coastal area, which can be seen at the Pinets beach. This original deposit is mainly composed of rocky shore with Posidonia oceanica banquettes (that are the characteristics that Fustera The study case of "Cala Fustera" (No. 51, Table A1, Figures 7 and A4) constitutes another example of significant coastal alteration after nourishment works. In this case, the previous sand deposit was replaced in 2017 by artificial quarry sand. Figure 7 shows visually differences in the colour of the coast of the Fustera Beach (now composed of quarry sand) with respect to the native colour in this coastal area, which can be seen at the Pinets beach. This original deposit is mainly composed of rocky shore with Posidonia oceanica banquettes (that are the characteristics that Fustera Beach should have). Unfortunately, no sand sample from the native Fustera Beach deposit was available. However, colour differences have been assessed assuming that the original aspect of the beach was similar to the sand colour of adjacent beaches. In this case, the chromatic change parameters take values of ∆h = 0.01, ∆C* = 8.13 and ∆E = 20.55, indicating that the new beach appearance is much whiter and bright than the original aspect.  A third example of sand colour modification due to management actions is the case of El Puerto beach (sample No. 4) located in Pilar de la Horadada, Alicante. In this case, sand from the El Puerto beach (a site next to a harbour from where sand is usually dredged, Figure 8) was transferred to other beaches such as "El Mojón" or "Mil Palmeras". Figure 8 shows the contact between the old and new coastal deposit. The colour difference of the added sand can be justified because the southern beaches are grouped at the top of Figure 5 with important differences in a* and b* parameters. A third example of sand colour modification due to management actions is the case of El Puerto beach (sample No. 4) located in Pilar de la Horadada, Alicante. In this case, sand from the El Puerto beach (a site next to a harbour from where sand is usually dredged, Figure 8) was transferred to other beaches such as "El Mojón" or "Mil Palmeras". Figure 8 shows the contact between the old and new coastal deposit. The colour difference of the added sand can be justified because the southern beaches are grouped at the top of Figure 5 with important differences in a* and b* parameters.

Conclusions
This research is a new contribution to the diverse geological heritage of the southeast of the Iberian Peninsula. Sand colour of natural beaches is influenced by the geological context in which they are located: External and Internal zones of the Betic Cordillera, and the Bajo Segura basin. These three main source areas of sediments determine natural sand composition and thus, the native beach colour. Because of this, the resulting colour of each beach is the consequence of their mineralogical composition linked to the source and hydrographic basin. In this sense, natural beaches related to the geological context of the Bajo Segura basin are clearer than the rest due to the high content in quartz grains, especially the southernmost ones. On the contrary, the darkest beaches are located in Tabarca Island, whose rocky outcrops mainly belong to the Internal zones of the Betic Cordillera. Therefore, although geology is the basis of the native beach colour, human factors can modify sand characteristics, as in the case of several studied beaches, especially those related to the geological context of the External Zones of the Betic Cordillera.
To measure the colour in sandy sediments it is adequate to apply CIEL*a*b* methodology to objectively quantify the chromatic appearance of each sample. In the study area, generally nourished/altered beaches differ from the native ones. The alteration of the sand colour by anthropic actions transformed the natural colour increasing green and blue tones (a* and b* colour components) and Lightness (positive values in L*). Despite this, a big group of natural beaches was identified with the highest a* and b* values (towards red and yellow tones). Colour and mineralogical parameters are related among them, such as the Lightness and the content of leucocratic components. Many beaches with native sand present a high content of melanocratic components (>60%), while most of the nourished beaches normally show less than 60% of such components. Principal Component Analysis shows good groupings associated with the mineralogical composition (related to the geographical location) and anthropic vs natural beaches. The position of the samples in the PCA graph is very useful to compare the sand colour from different beaches and also to interpret the limited information about nourishment works carried out in the past.
The 66 investigated samples differ in colour parameters and the beaches that receive more tourists have been nourished and tend to be clearer (e.g., "Ponent", "El Postiguet", "San Juan"). Results show that a few nourished beaches, e.g., "Playa del Cura", "Llevant", "Arenal Bol", "El Postiguet", and "Ponent", have a similar colour to the nearby native ones. The rest of the nourished beaches have lower values of a* and b* colour parameters, moving towards green and blue tones,

Conclusions
This research is a new contribution to the diverse geological heritage of the southeast of the Iberian Peninsula. Sand colour of natural beaches is influenced by the geological context in which they are located: External and Internal zones of the Betic Cordillera, and the Bajo Segura basin. These three main source areas of sediments determine natural sand composition and thus, the native beach colour. Because of this, the resulting colour of each beach is the consequence of their mineralogical composition linked to the source and hydrographic basin. In this sense, natural beaches related to the geological context of the Bajo Segura basin are clearer than the rest due to the high content in quartz grains, especially the southernmost ones. On the contrary, the darkest beaches are located in Tabarca Island, whose rocky outcrops mainly belong to the Internal zones of the Betic Cordillera. Therefore, although geology is the basis of the native beach colour, human factors can modify sand characteristics, as in the case of several studied beaches, especially those related to the geological context of the External Zones of the Betic Cordillera.
To measure the colour in sandy sediments it is adequate to apply CIEL*a*b* methodology to objectively quantify the chromatic appearance of each sample. In the study area, generally nourished/altered beaches differ from the native ones. The alteration of the sand colour by anthropic actions transformed the natural colour increasing green and blue tones (a* and b* colour components) and Lightness (positive values in L*). Despite this, a big group of natural beaches was identified with the highest a* and b* values (towards red and yellow tones). Colour and mineralogical parameters are related among them, such as the Lightness and the content of leucocratic components. Many beaches with native sand present a high content of melanocratic components (>60%), while most of the nourished beaches normally show less than 60% of such components. Principal Component Analysis shows good groupings associated with the mineralogical composition (related to the geographical location) and anthropic vs natural beaches. The position of the samples in the PCA graph is very useful to compare the sand colour from different beaches and also to interpret the limited information about nourishment works carried out in the past.
The 66 investigated samples differ in colour parameters and the beaches that receive more tourists have been nourished and tend to be clearer (e.g., "Ponent", "El Postiguet", "San Juan"). Results show that a few nourished beaches, e.g., "Playa del Cura", "Llevant", "Arenal Bol", "El Postiguet", and "Ponent", have a similar colour to the nearby native ones. The rest of the nourished beaches have lower values of a* and b* colour parameters, moving towards green and blue tones, e.g., "Punta Prima", "Marineta Cassiana", "Playa Lisa", "Bernabeu II", "Calas del Cuartel", etc. Despite this, native clear samples generally presented higher values in the a* and b* parameters, this giving a more golden visual appearance, e.g., Ferri Beach, Cervera Beach, "Los Locos", "Les Ortigues", "La Mata" or "Bon-nou" among others. Today, the beach colour is not taken into consideration in some nourishment works, at least in the study area where "El Puerto" in the south or Fustera Beach in the north are clear examples.
This study showed with analytical data a post-nourishment change in sand colour, a common consequence also produced in other coastal sites such as Italy, Belgium or Cuba. Such changes may not be appreciated by national and international tourists as well as by old local beachgoers. The CIEL*a*b* colour space used to characterise sediment colour and Lightness is a useful tool to avoid significant changes in sand colour that can affect coastal scenery, beachgoers preferences, and the biological processes of some species. Native sediment colour is a natural characteristic of the coastal landscape and for this reason it should be maintained, at least in rural beaches and natural parks. It is essential to preserve the most natural coastal areas such as the dunes of Guardamar or the Island of Nueva Tabarca, considered by the Spanish Geological Survey as places of geological interest, mainly due to their geomorphological, sedimentological and stratigraphic particularities. In urban beaches, where all is artificial, a change of sand colour can be accepted in order to improve the beach attractiveness. Regarding future research prospects on the colour of beach sand, it would be interesting to analyse the sand properties against time in order to detect possible changes in the coastal environment.
Sand is a limited and increasingly valuable resource. Its use and relocation generates strong environmental, economic, and social impacts, especially on the coast. For this reason, it is preferable to opt for dune restoration and the elimination of unnecessary structures/constructions (on the coast or inland) that promote beach erosion instead of proposing nourishment projects that involve large economic investments and that may affect different coastal and marine habitats as well as associated species. Finally, the method used in this paper is easy to apply and can be very useful to carry out appropriate erosion management actions by means of beach nourishment works in the future. Consequently, the current conclusions are of particular interest to future actions to preserve the natural environments and original characteristics of native deposits.