Environmental Changes and Cultural Transitions in SW Iberia during the Early-Mid Holocene

: The SW coast of the Iberian Peninsula experiences a lack of palaeoenvironmental and archaeological data. With the aim to ﬁll this gap, we contribute with a new palynological and geochemical dataset obtained from a sediment core drilled in the continental shelf of the Algarve coast. Archaeological data have been correlated with our multi-proxy dataset to understand how human groups adapted to environmental changes during the Early-Mid Holocene, with special focus on the Mesolithic to Neolithic transition. Vegetation trends indicate warm conditions at the onset of the Holocene followed by increased moisture and forest development ca . 10–7 ka BP, after which woodlands are progressively replaced by heaths. Peaks of aridity were identiﬁed at 8.2 and 7. 5 ka BP. Compositional, textural, redox state, and weathering of source area geochemical proxies indicates abrupt palaeoceanographic modiﬁcations and gradual terrestrial changes at 8.2 ka BP, while the 7.5 ka BP event mirrors a decrease in land moisture availability. Mesolithic sites are mainly composed of seasonal camps with direct access to the coast for the exploitation of local resources. This pattern extends into the Early Neolithic, when these sites coexist with seasonal and permanent occupations located in inland areas near rivers. Changes in settlement patterns and dietary habits may be inﬂuenced by changes in coastal environments caused by the sea-level rise and the impact of the 8.2 and 7.5 ka BP climate events.


Introduction
The Early-Mid Holocene period is characterized by several environmental changes, such as the postglacial marine transgression, an increase of temperatures, and the existence of diverse episodes of ice-rafted debris in the North Atlantic Ocean. Littoral areas suffered a rapid transformation resulting from the inundation of the fluvial valleys and the progressive development of diverse coastal features. From a different perspective, the vegetation adapted to the new climatic-edaphic conditions and reacted to abrupt episodes, the so called Bond events, which are associated to dry conditions in the Mediterranean region [1]. However, the Early-Mid Holocene was not only defined by unique environmental changes, but by important cultural transformations within hunter-gatherer groups and the adoption of new subsistence strategies.
The process of neolithization in Southwestern Iberia is controversially disputed and different theories attempt to explain its origin. Some authors interpret it as a maritime pioneer colonisation of depopulated regions by Neolithic groups from the Western Mediterranean [2], while others propose indigenous origins defined by the continuity of local Most of all the Mesolithic sites are found in the Cape São Vicente region. They constitute shell-middens of different sizes (Castelejo, Barranco das Quebradas, Rocha das Gaivotas) and open air sites (Armação Nova, Vale Boi) interpreted as seasonal occupations based on the exploitation of coastal and/or lithic resources [12,13]. Near the Straits of Gibraltar, the only site for this period is considered as a single occupation based on the exploitation of marine resources and hunting (Palmones), although there is some controversy about its chronology [14,15]. The majority of these sites is located close to marine environments or in transitional areas between fluvial and marine influences.
Several Neolithic sites have been discovered in both areas. In the Cape São Vicente region, shell-middens (Alcalar 7, Ribeira de Alcantarilha) and temporary open-air sites (Vale Santo I, Padrão) based on the seasonal exploitation of marine and lithic resources are still present [12]. Moreover, some Mesolithic sites revealed later Neolithic layers (Castelejo, Barranco das Quebradas, Rocha das Gaivotas, Vale Boi) [12,13,16]. Two necropolis caves (Ibn-Ahmmar, Algarão da Goldra) were also assigned to this period. There is one site considered a residential basecamp (Cabranosa), a settlement suggesting an obvious sedentary occupation with hut structures and grave burials (Castelo Belinho), and a possible long-term settlement in a beach environment (Praia do Forte Novo) which is considered a salt production site by some authors [17]. Many of these new Neolithic sites are located in fluvial environments some kilometres inland, if compared with those belonging to the Mesolithic (Alcalar 7, Ribeira de Alcantarilha, Castelo Belinho, Ibne-Ahmmar, Algarão da Goldra). In the southeastern area of the Gulf of Cádiz it is possible to identify varying typologies of archaeological sites with Neolithic occupations: a single human occupation (El Retamar) based on the exploitation of marine resources, hunting, and animal domestication [18,19], a permanent settlement (Campo de Hockey) with a large necropolis [20][21][22], and a necropolis cave (La Dehesilla) [23]. With the exception of the cave situated in a mountainous area, all other sites are located in the marine-fluvial environments of the Cádiz Bay.

Regional Settings
The study area comprises the southwestern coastal region of the Iberian Peninsula framed by the Atlantic Ocean and located north of the Africa-Iberia plate boundary. The westernmost part is a Meso-Cenozoic sedimentary basin overlying Carboniferous basement, which extends offshore as far as 100 km south [24]. However, most of the coast is dominated by Neogene basins defined by the presence of several estuaries conforming littoral lowlands. Some of them are sheltered by spits, with marshlands extending several kilometers inland [25,26]. These characteristics are mainly due to the morphology of the coast and the prevailing winds from SW that generate a littoral drift towards the East and SE of the Gulf of Cádiz, being responsible for the formation of beach and barrier systems that partially close some estuaries [25,27].
The SW Iberian Peninsula has a Mediterranean climate with oceanic influence that can be defined by two climatic zones described by Köppen-Geiger [28]: a temperate climate with hot summers in oriental and inland areas (Csa), and a temperate climate with warm summers in occidental areas and the littoral (Csb). Considering both zones, the annual average minimum temperature oscillates between 12.5 • C and 15 • C, and the annual average maximum temperature ranges from 22.5 • C to 25 • C [29]. The average total annual precipitation fluctuates between 500 mm and 1000 mm, mostly concentrated during winter season [29] The thermomediterranean belt of the coastal lowlands is defined by different ecosystems and vegetal ecological associations. Woodlands are mainly composed of evergreen (Quercus suber, Q. rotundifolia) and deciduous oaks (Q. faginea, Q. estremadurensis), together with some Mediterranean species, such as Olea europaeae subs. sylvestris or Juniperus oxycedrus, among others. Pre-forest scrubs are dominated by Quercus coccifera, Phillyrea angustifolia, Pistacia lentiscus, Rhamnus sp., and Juniperus sp, while riparian communities are formed by Alnus glutinosa, Salix sp., Fraxinus angustifolia, and Tamarix sp. [30][31][32][33]. Heathlands are well represented by different genera of Ericaceae that grow on acidic low fertility soils, under high humidity levels and strong oceanic conditions [34]. In littoral areas and interdunal valleys the vegetation is mainly composed of Pinus pinaster, P. pinea and P. halepensis. The most representative herbaceous vegetation of the stable dunes are: Polygonum maritimum, Artemisia crithmifolia, Eryngium maritimum, Cyperus capitatus, Anthemis maritima, Silene sp., Juniperus sp. [35]. Freshwater marshlands include species of Cyperaceae and Isoetes spp., Typha sp., Hypericum tomentosum, and Phragmites australis, while saltmarshes are dominated by species of Chenopodiaceae and some taxa such as Apium graveolens, Aster tripolium, Limonium, Juncus maritimus, and others [30,36,37].

Materials and Methods
Twin sediment cores 19-01 and 19-02 (referenced as GeoB23519-01/-02 in the MARUM core repository in Bremen, Germany) were retrieved in the Algarve continental shelf (37 • 00.656, 008 • 52.247 and 37 • 00.654, 008 • 52.26, respectively) at ca. 65 m water depth and approximately 3 km from the coast using a vibracorer. The 3.64 m-long core 19-01 was sampled for palynological analysis and the neighboring core 19-02 (3.62 m-length) is added to increase data and resolution on sedimentation rates and chronology. In both cores, stratigraphic units were identified by visual description of the sediment.
From bottom to top (Figure 2), both cores are composed of greenish fine sands covered by muddy sediments and capped by coarse silt to silty sand (uppermost 1.5 m). Shells are common as well-preserved specimens and the clastic fraction is composed by bioclasts, siliciclastics, and lithic fragments.
X-ray fluorescence scans (XRF) were measured using an ITRAX core scanner (resolution 2 mm, 20 s exposure time, 30 kV, 55 mA). Carbonate (Ca, Sr), siliciclastic (Si, Al, Ti, Zr, Rb), organic matter (Br, Mo incoherent and coherent scatter [Mo inc, Mo coh]), and redox geochemical (Mn, Fe) proxies were chosen ( Figure 2) and some useful ratios were analyzed (Ca/Si, Si/Al, Zr/Al, Zr/Rb, Mn/Fe, Mo inc/coh) [38]. Anomalous values were removed by comparison to the core photographs to minimize the effect of coarse grains, cracks, and surface irregularities. The remaining data were normalized by count per row to unit by column and smoothed with an arbitrary 11-point running mean filter to reduce the inherent data noise.
Seven articulated bivalves from core 19-01 and six more articulated valves from 19-02 were 14 C-dated at the Beta Analytic Inc. and the University of California-Irvine (USA). Calibration was performed with CALIB 8.2 [39] using the Marine20 calibration dataset [40] and different ∆R calculated from the original datasets of [41][42][43][44][45] by using the online application of [46] (Table 1). After comparison of the sedimentological and geochemical logs, it was possible to correlate both cores by depth. The age model is based on a simple linear interpolation and sedimentation rates were calculated for each dated interval ( Figure 3).
Palynological samples were collected from core 19-01 with an interval of 3-7 cm. All samples were chemically treated with HCl to remove carbonates, KOH to remove humic acids, and sodium polytungstate (SPT: 3Na 2 WO 4 ·9WO 3 ·H 2 O) at 2.0-2.1 cm 3 for densimetric separation. The final residue obtained after the treatment was mounted on slides with the use of glycerol mixed with phenol. Palynomorphs were counted using an optical microscope at 400× and 1000× to a minimum pollen sum of 150 terrestrial pollen grains.
Fossil pollen grains, spores, and non-pollen palynomorphs were identified using published keys [47][48][49][50][51][52] and the modern pollen reference collection of the CSIC in Madrid (Spain). Microcharcoal particles >125 µm were counted alongside the identification of pollen grains and interpreted as indicators of regional fires [53]. Pollen and microcharcoal concentrations (grains/gr and particles/gr of dry sediment, respectively) were estimated by adding one Lycopodium clavatum tablet to each sample [54]. Diagrams were plotted versus age using TiliaIT software (version 2.    Pollen percentages for terrestrial taxa were calculated against the main sum of terrestrial grains, while percentages for aquatics, spores, and Pinus were calculated against the total sum of all pollen and spores. Pinus grains were excluded from the main sum because they are considered to be over-represented in marine sediment cores because of its extensive dispersal ability and buoyancy (Hopkins, 1950). However, morphometric analysis of their pollen grains was done using a measurement of the grain diameter, excluding sacci [55][56][57].

Age-Depth Model
According to the proposed age model, the period between 12 ka and 6 ka BP is placed in the lowermost 1.5 m of the core. From ca. 12 ka BP until ca. 8 ka BP, the sedimentation rate increases from 0.06 mm/yr to 0.69 mm/yr which corresponds to the maximum value for the Holocene in this core. From ca. 8 ka BP to 6 ka BP the depositional rate decreases to 0.27 mm/yr ( Figure 4a).

Geochemistry
The selected elements can be grouped by similar behavior along time allowing the identification of some common periods and events despite the considered parameter.
Zr/Rb ratio has been linked to grain size changes [58] and sediments older than ca. 8.1 ka BP show larger values and a clearer decreasing trendthan younger sediments ( Figure 4b).
Ti, Fe, and Mn run almost parallel during the studied period and they can be linked to weathering of siliciclastic rocks [59]. These elements show reverse trends and peaks as compared to Zr/Rb (Figure 4c), showing large ratios and enhanced increase for the 8.2-6 ka BP period than for the Early Holocene.
Al, Si, Ca, and Sr can be correlated in the long term ( Figure 4d), while they show some differences for the shorter cycles. Al and Si are interpreted as representative of terrestrial input [38]. Comparable Ca and Sr trends represent carbonates [38], and reflect high bioclastic content. All elements are components of the dominant clastic fraction, while siliciclastics depend on the terrestrial supply, carbonate variability can be due to changes on marine productivity (in situ biogenic carbonate) and this could be responsible for the out-of-phase signal for the shorter cycles. All of them show higher values for the 8-12 ka BP period, with larger amplitudes for Si, Ca, and Sr compared to the 6-8 ka BP period, when Al shows wider amplitudes reflecting the clay-richer composition of the sediments ( Figure 2).
It is worth to mention the long-and short-term correlation of Br, Mo inc/coh, and Fe/Mn (Figure 4e). Br and the Mo inc/coh ratio have been used as a proxy of organic matter in marine sediments [60][61][62]. The Fe vs Mn ratios can indicate changes in the redox state [63,64], with Mn-displaced values indicating oxic conditions. All these values show a steady increase, broken by a sudden rise from 8.3 to 8 ka BP. However, they are lower for Br, mirroring the changes in clastic fraction (Al, Si, Ca, Sr).
Ca/Si ratio points to the carbonated vs. siliciclastic sources [65,66], while Si/Al and Zr/Al are used as textural proxies like sorting or grain size [38,67] (Figure 4f). These ratios show no clear trends, as those mentioned before, but they can be divided in three stages: from 12 until 8.8 ka BP, their record show low amplitude cycles that shorten their period in time, Ca/Si decreases while Zr/Al and Si/Al remain, showing a slowly increasing trend; from that time until ca. 7.8 ka BP, the amplitude increases while their period shortens, both Ca/Si and Si/Al increase but Zr/Al shows its average lowest values; from ca. 7.8 ka BP to 6 ka BP, the period of the cycles grows, while their amplitude slightly decreases without reaching values of the previous core section. Ca/Si and Al/Si values show a falling trend while Zr/Al remains nearly stable.
According to these observations, the 12-6 ka BP lapse can be split in three periods ( Figure 4): The sedimentation rate was the highest. There are marked and sudden changes for a decrease in clastic input and an increase in organic components, while the decrease in grain size and increase in weathering proxies is fast but gradual. • 8-6 ka BP. The sedimentation rate decreases from 0.69 to 0.27 mm/yr. Clastics and grain size fall to their minimum average values but amplitude of changes increases for sorting (grains vs. matrix) proxies. Organic matter preservation increases, coincidental to sub-oxic to anoxic conditions, as well as the weathering of source areas.

Geochemical Trends and Events, Sea Level and Palaeoceanographic Changes
Comparison of the evolution of geochemical and sedimentary parameters against regional and global forcings (sea level, climate) ( Figure 6), allows us to correlate them and to identify links among the environment, these forcings, and human populations.
The ca. 10-6 ka BP period is known as the Holocene Thermal Maximum (HTM) [68][69][70][71], a period of warm and humid climate as compared to previous and following periods and only interrupted by a short cooling period known as the 8 ka or the 8.2 ka BP event [72,73] (Figure 6a).
Sea level rose from its previous minimum during the Last Glacial Maximum with maximum rates for the Algarve between 11-8 ka BP (from 9 to 7.4 mm/yr), that decreased between 8-7 ka BP (4.5 mm/yr), and another time from 7 to 6 ka BP (2.0 mm/yr) and fell to 0.7 mm/yr for the 6-5 ka BP period [74]. Coeval to this change, sedimentary records of core 19-01 gradually decrease in grain size (Zr/Rb ratio, Figure 6b) that can be interpreted as the result of deepening (drowning by flooding), increasing distance to the wave base, and the action of coastal currents and sedimentation rates changed accordingly (Figure 6c).
Sea surface temperature decreased in the Gulf of Cádiz [75] and Alborán Sea [76], and this change correlates with the increase in the organic matter content of core 19-01 (Figure 6d).
Weathering proxies are indicative of wet and warm conditions which prevail after 8.2 ka BP (Figure 6e). These conditions are coherent with the alleged characteristics of the HTM and the African Humid Period [77,78]. Stumpf et al. (2011) [78] showed a decrease in the illite/kaolinite ratio and interpreted this as related to a change in dust sources, but it is also coherent with an increase in weathering of the source areas.
The decrease in clastic fraction (Figure 6f) can also be linked to the increase in weathering and the growing water depth. In addition, the continuity of terrestrial influx (Si/Al, Zr/Al, Figure 6g) vs. the gradual increase in organic content and the abrupt and short increase in carbonates around 8.2 ka BP (Ca/Si, Figure 6g) indicate a rise in marine productivity.
However, one of the main features of the geochemical record of core 19-01 is the abrupt change in marine productivity and clastic dynamics and gradual weathering of source areas along with waves/shallow currents. These changes relate to the 8.2 ka BP event. Thornalley et al. (2009) [79] interpret that there was an abrupt switch from a stratified upper ocean to well-mixed waters around 8.4 ka ago caused by changes in the relations between the subpolar and subtropical gyres. Bazzicalupo et al. (2020) [80] showed from records of the Alborán Sea, that present oceanic gyres system developed around 8 ka ago. Also, models results indicate the changes in oceanic gyres around 8.2 ka BP as being responsible for the present day North Atlantic circulation [81]. Thus, it seems possible that the observed abrupt changes in the geochemical record of core 19-01 can be related to those palaeoceanographic changes, while the gradual ones are interpreted to link to terrestrial climate variations. Algarve and Central Portugal [74] and global [85] and Zr/Rb ratio (notice the reversed scale); (c) sea surface temperature (SST) reconstructions for the Gulf of Cádiz (alkenone Uk37 SST derived for core M39-008 [75]) and Mg/Ca SST derived for the Alborán Sea (ALB2 [76]); (d) Ti, Fe, and Mn, weathering proxies, as indicative of land moisture; (e) Al, Si, Sr and Ca, representative of clastic sedimentation, serve as proxies for activity of marine currents; (f) Ca/Si, carbonate vs. siliciclastics, and Zr/Al and Si/Al, proxies of siliciclastic grains vs. matrix, as indicative of land derived vs. bioclastic (marine) inputs; (g) sedimentation rate for core 19-01 for the 12-6 ka BP period.

Vegetation Trends and Evolution of Coastal Environments
During the Early Holocene, the post-glacial marine transgression was progressively drowning the exposed continental shelf of the southwestern Iberian coast with a rapid sea level rise between ca. 13 ka BP and 7-6.5 ka BP, although with different responses depending on the geographic location and topography [27]. Initially, until ca. 10 ka BP, the first stages of the river valleys inundation resulted in pre-estuarine palaeovalleys and transitional areas of fluvial-saltmarsh deposits as the main coastal landscapes [25,[86][87][88][89]. For this period until ca. 9.2 ka BP, the core 19-01 indicates the existence of forests composed of mesophilous trees, and Mediterranean trees and shrubs which defined the onset of the Holocene as a warm phase. Mediterranean and high mountain pines were also present, although in low values, being that the latter (Pinus sylvestris-nigra type) may be a reminiscence of the continental climate of the previous glacial period, as reported in other deposits [90]. These conditions have been also inferred in continental deposits in the Guadiana Estuary and the Medina Lagoon, where the presence of open land indicators and high-mountain pines were recorded, respectively [91,92]. In this regard, it has been found that the hydrological response in some Western Mediterranean areas during the Early Holocene maintained prolonged arid conditions, with shallow lake levels followed by an increase of moisture ca. 10-9 ka BP [93]. This may have been linked to insolation and seasonality changes due to the orbital variability, the presence of the Laurentide and Fennoscandian ice sheets in the North Atlantic Hemisphere, and a series of ice rafted debris (IRD) provoking immediate ocean surface coolings associated with dry conditions in the Mediterranean region [1,94,95]. Some of these events fall into this period ( Figure 7) and two of them appear to have triggered a reaction in the vegetation with a decrease of mesophilous trees and an increase of Asteraceae immediately after 10.3 ka BP, and a retraction of mesophilous and Mediterranean forests with a parallel rise of Mediterranean shrubland and xerophytic taxa after 9.4 ka BP. Diverse records from Greenland show anomalies in the mean δ 18 O curves ca. 9.95 and 9.3 ka BP interpreted as temperature irregularities, but also to changes in moisture sources and/or transport paths [96,97]. The slight delay in the vegetation reacting to these events indicates that changes in terrestrial landscapes were not immediate.
Regarding the evolution of the coast, the period between ca. 10-7 ka BP is also coincident with a phase of intense marine influence due to the marine transgression into the river valleys of the SW Iberian margin, transformed in estuaries by drowning [25,[86][87][88][89]. As a result, salt marshes developed along the littoral of the Gulf of Cádiz and a landward shift of the boundary between marine and fresh waters took place, with a consequent increase of saline environments [25]. High values of foraminiferal linings were recorded between ca. 9.5 to 8 ka BP, probably related to the development of estuaries. From ca. 8 ka BP onwards and abrupt drop on the foraminiferal linings gives way to a progressive increase of Isoetes, characteristic of seasonal fresh marshland environments [98,99]. The increase of wetland taxa seems to parallel the slight rise of riparian communities registered during this period, which may be a consequence of the influence of the river in the hinterland, as already recorded in other continental cores [86,100].
At ca. 10-9 ka BP an abrupt increase in moisture was identified in several lake records of the Western Mediterranean [93], and the progressive expansion and maximum values of oakwoods between ca. 10 to 7 ka BP seem to confirm this trend. Despite slight differences in the chronologies, this rise of temperate and Mediterranean forest taxa is also registered in other marine cores drilled in the Atlantic margin of southwestern Iberia [101][102][103], as well as in some continental deposits of this area [91]. During this period, two peaks of microcharcoals were recorded and interpreted as episodes of increased regional wildfires at 8.8 and 8.4 ka BP (Figure 7). Several forest contractions were also identified at different points and with diverse characteristics but it is noteworthy to highlight two episodes. At 8.2 ka BP there is an increase of xerophytic elements and afterwards an abrupt drop of the forest taxa took place, culminating with low values of mesophilous/Mediterranean trees and a rise of Asteraceae ca. 8.1 ka BP. Greenland ice cores reflect well-defined anomalies in the period between 8.4-8 ka BP with very low values of δ18O around 8.2 ka BP [97]. Although its origins remain unclear, the 8.2 ka BP event has been linked to a fresh water influx into the North Atlantic that would have provoked changes in temperatures and thermohaline circulation and thus, in moisture availability [1,104]. Indeed, the Medina Lagoon record shows a lake-level decrease for the period between ca. 8.5-7.8 ka BP, which has been related to global climatic instability centered on 8.2 ka BP, concluding in a desiccation phase [105]. Events of forest decrease were also recorded in the SW Atlantic margin between ca. 8.6-8 ka BP, as well as in the Alborán Sea between ca. 8.3-8 ka BP [101][102][103]106,107]. In continental deposits, some forest setbacks and expansion of open-ground taxa were recorded in the Guadiana Valley and the Medina Lagoon at ca. 7.8 and 8.2 ka BP, respectively [91,108]; however, several sequences do not mirror any vegetal changes for this period [15,109].
Another important crisis was also identified in the GeoB235-19-01 core at ca. 7.5 ka BP, defined by a visible peak of xerophytes and an abrupt drop of mesophilous trees, riparian woodland, and Mediterranean woodland and shrubland, pointing to a rapid decrease of moisture availability (Figures 6e and 7). Diverse episodes of forest contractions were identified between 7.5 and 7 ka BP in different cores of the Atlantic Iberian Margin and the Alborán Sea, as well as in continental sequences [91,[101][102][103]106,107,109]. Diverse palaeorecords from SE Iberia suggest changes towards increased aridity between 7.8-7.3 ka BP [110].
Regarding the coastal evolution, the core 19-01 recorded the maximum sea level rise rates for the Algarve between 11-7 ka BP, but in most of the estuaries located in this area (Alvor, Alcantarilha, Quarteira, Carcavai, Tinto & Odiel, Guadalete) the maximum flooding of the river valleys was registered at ca. 7-6 ka BP, resulting in the development of open estuaries, inundated channel banks, tidal flats, and even the landward migration of estuarine barriers [25,[86][87][88][89]. The fluctuation between marine vs freshwater markers recorded in the core 19-01 seems to reflect the marine pulses and the diverse evolution of littoral features and edaphic conditions, depending on factors such as the accommodation space of the river valleys. During this millennium, mesophilous and Mediterranean woodland progressively decrease, while Mediterranean shrubland seem to stabilize at the same time that heaths start a stepwise expansion from this point onwards. An important increase of Cichorioideae occur ca. 6 ka BP, correlating the rise of Isoetes and hygro-hydrophytic herbs. The ambiguous significance of Cichorioideae makes it difficult to understand this episode, and its presence may suggest an important development of marsh/wetland communities, or either the development of a semi-open landscape paralleling the expansion of shrublands to the detriment of forests.
After ca. 6 ka BP, rates of sea level rise fell to 0.7 mm/yr in the Algarve [74], while in most of the estuaries of this area rates are of ca. 2.6 mm/yr and oscillations below 1 m defined the last 5000 years [27,111,112]. Some coastal elements like lagoons, spits, and barriers developed during this period [113].

Human Groups in Dynamic Littoral Habitats
The rapid sea level rise occurring from ca. 13 to 7 ka BP was one of the main influential factors affecting the evolution of the coastal morphology and, therefore, human settlements. It is assumed that if the currently known Mesolithic sites were distributed along the coast, other hunter-gatherer occupations may have existed in the exposed continental shelf before its flooding due to marine transgression, which in that case may have been erased by the sea. The absence of any human settlement previous to ca. 9 ka BP in the studied area, with exception of one date in Barranco das Quebradas [13], is a clear example of the consequences of the sea level rise in this region.
From ca. 9 ka BP onwards, all the preserved Mesolithic sites are located at or above 50 m asl, except Castelejo and Palmones (supplementary material), and with direct access to the sea, considering the coastline at that time ( Figure 8). The lack of data in the littoral of the Gulf of Cádiz during this time span is critical, but it is especially dramatic in the central and eastern sector. The fact that most of the preserved sites are located in the cliffy area could be explained by a greater spatial impact of the sea level rise on the coastal plains at that time opened to the sea, since most of the spits and barriers developed after ca. 6 ka BP [27,[111][112][113]. For the period between ca. 9-7.5 ka BP, recent palaeodemographic studies propose a steady growth in population reflected in a concentration of archaeological sites in Central-South Portugal (Muge and Sado estuaries) and across the Algarve, although here in a more dispersed pattern [114]. Most of the sites concentrated in Western Algarve are shellmiddens and with exemption of Vale Boi, where a human tooth dating of Mesolithic but without a funerary context was described [16], all the occupations in this region are seasonal camps focused on the collection of local resources, mainly marine food but also flint material [13]. Some hypotheses correlate these sites with the existence of some (semi)permanent basecamps on which they would depend and that would have similar characteristics to those from the Sado and Tejo estuaries, with elaborated habitational structures, human burials, and broad-spectrum subsistence strategies [6,12,115,116]. However, no (semi)permanent basecamps have been found in the Algarve region yet, and it seems to be a hard task considering the high level of human impact of this area during historic times. Regarding the eastern sector of the Gulf of Cádiz, environmental differences, and the specific cultural trajectories in the "Alborán territory" [8], the nearest area where most of the pre-Neolithic sites are located, prevent us from establishing any comparison between archaeological sites of both zones. Some authors consider that major changes in the social, technological, economic, and in settlement patterns of hunter-gatherer groups did not occur with the transition from the Late Pleistocene to the Holocene, but in relation to the 8.2 ka BP event. An example would be the development of the complex shell-midden basecamps located in inland areas following the course of the Tejo and Sado Rivers, as opposed to the seasonal shell-middens located in the Algarve coast. This change of settlement patterns is explained as the result of a decline in the availability of the marine resources due to the joint action of a drop in the upwelling intensity together with an increased aridity ca. 8.2 ka BP [116]. This would have made the palaeoestuaries of these rivers stable resource-rich brackish environments, unlike the coastal unstable ecosystems adversely affected by this episode [116]. Likewise, results from the Pena d'Agua rock-shelter show important changes in the acquisition of raw material and technological features, which would be explained by modifications in mobility patterns from residential to logistic as an adaptive response to changing environments from 8.2 ka BP onwards [117]. Considering this hypothesis, these changes may have accelerated the confluence of human groups to settle in areas with stable and available resources, resulting in a (semi)sedentary pattern and social complexity [116][117][118]. However, there is no sign of the impact of the 8.2 ka BP event in the SW Atlantic Iberian margin in terms of settlement patterns. Notwithstanding, the bias caused by the lack of information needs more research in the study area to discern whether the seasonality of these sites reflect important regional differences, or if it is related to the existence of some (semi)permanent undiscovered basecamps.
Regarding the diet, not only shell-middens from Western Algarve display a specialization in the exploitation of marine resources, but also the seasonal occupation of Palmones revealed an intense gathering of shellfish [15]. Isotopic analysis of the Vale Boi teeth showed that, although marine resources were a prominent element, the diet of these Mesolithic groups was based on a mix of terrestrial and aquatic resources [12]. However, there may have been important differences between sites, because isotopic data from the Sado shell-middens corroborated frequencies of marine diet ca. 25% [119], while in the Muge area the marine diet is above 50% [120].
In the occidental area of the Algarve, topshells were the most abundant species consumed during the Early Mesolithic, while limpets and gooseneck barnacles were preferred during the Late Mesolithic and Early Neolithic [121]. A reduction in the size of the shellfish collected was also identified in the transition from Mesolithic to Neolithic at Vale Boi, Rocha das Gaivotas, and Armação Nova, and it was considered a consequence of overexploitation [12] or a decrease in the foraging efficiency [121]. However, in some Mediterranean areas such as the Pego palaeolagoon, a reduction of lagoon bivalve size during the Late Mesolithic has been observed, related to a decrease in resource productivity [122]. Despite this there was gently increase in ocean organic productivity linked to the slow overall cooling of the waters during the Mid Holocene, the warmer climate during the Holocene Thermal Maximum caused a weakening of the coastal upwelling, briefly interrupted during the short cooling episodes [79]. This alteration in the coastal palaeoproductivity suggests that the modification in consumption patterns and the size decrease of shellfish in SW Iberian sites may have been a consequence of these environmental alterations rather than being caused by cultural factors.
Commonly, the expansion of the Neolithic in both the western and eastern sectors of the Gulf of Cádiz is considered to have occurred around 7.5 ka BP [9,12]. Corresponding to the period, different sequences reflect an aridity crisis affecting the vegetation of the SW Atlantic Iberian coast and the Alborán Sea [91,[101][102][103]106,107,109] and the geochemistry dataset reflects a decrease of land moisture (Figure 6e). Some Mesolithic sites exhibit layers of Neolithic occupations (Castelejo, Barranco das Quebradas, Rocha das Gaivotas, Vale Boi, Palmones), although with hiatuses of several hundred years. The existence of several seasonal camps (Alcalar 7, Ribeira de Alcantarilha, Vale Santo I, Padrão, El Retamar) suggests the maintenance of previous subsistence strategies until this period. However, Neolithic times also involve new settlement patterns with sites located near rivers, shell deposits thinner than those in the Mesolithic, and caves. Only two sites suggest clear sedentary occupation (Castelo Belinho, Campo de Hockey) and the presence of seasonal camps reveals that the mobility was still key in the provisioning of different resources. Although the aim of this paper is not to discuss the origins of the Neolithic in the study area, the coexistence of previous (Mesolithic) and new (Neolithic) settlement patterns would be coherent with the idea of an integration of the two populations, rather than with the disappearance of the hunter-gatherer way of life [120].
It is noteworthy to mention that except the pre-existing sites and some exceptions such as Praia do Forte Novo, El Retamar, and Campo de Hockey, all the Neolithic sites are located a few kilometres inland. Some changes in dietary habits regarding aquatic resources may be linked to this new settlement pattern, with the increased consumption of cockles and clams [121] typical from transitional brackish environments (rías), in contrast to the rocky-like shellfish preferred during the Mesolithic. Some ideas of what may lie behind this change in settlement patterns are: • Different populations would have a different perception of the same habitat, so they develop preferences for settling in distinct areas.

•
The maximum flooding of some river valleys would have forced human groups to settle in areas less close to the coastline. It should be also considered that archaeological sites located in the mouth of these rivers could have been buried or destroyed after this flooding episode and only those located in inland areas are preserved. Regarding the resource exploitation and subsistence strategies, the introduction of a farming economic system resulted in a major dietary shift towards terrestrial food during this transitional period [123]. In relation to the agricultural practices, the first evidences were dated ca. 7.5-7 ka BP in Eastern Andalusia and Central-South Portugal, while there is not any direct evidence of their origins in the Occidental Andalusia and the Algarve area [124]. Livestock and grazing activities have direct testimony in zooarchaeological material. With respect to palynology, anthropogenic indicators are most likely to be identified in continental and especially archaeological deposits than in marine cores. However, the first manifestations of agricultural activities are only documented ca. 6 ka BP in areas with previous and important Mesolithic occupations [124].
In relation to the preservation of human settlements, the lack of data in the central coast of the Gulf of Cádiz for both Mesolithic and Neolithic periods must be considered. Some archaeological sites and material identified in the surface [3] are ascribed to the Early Neolithic, but the absence of systematic excavations and datings makes it presently impossible to shed light on the evolution of human groups during this time span. However, these settlements display a similar pattern compared with those of the Algarve, located close to the coast but some kilometers inland and near the Tinto, Odiel, Guadalquivir, and Guadalete rivers. An important issue to consider is related to the several high energy wave events that have been identified in this area during this period. In the Guadalquivir estuary, several high energy events were identified between 9.9-9.2 ka BP and 8.2-7.8 ka BP, some of them of unclear origin, some identified as storms, and one of them at 9.1 ka BP of tsunamigenic origins, while ca. 7-6.8 ka BP a palaeotsunami occurred [109,125]. The westernmost Algarve area was also hit by a high energy wave event, only identified in the Alvor estuary ca. 6.4 ka BP [86]. Thus, not only the sea level rise but the hazardous nature of the coastal plain in this region, which seems to be constantly subject to these extreme events as well as affected by subsidence, may hamper the detection of potential archaeological sites.

Conclusions
This paper presents a new palaeonvironmental dataset that contributes to a better understanding of the evolution of coastal landscapes in the littoral of the Gulf of Cádiz. The multi-proxy approach combining palynological, geochemical, and sedimentological data has provided a framework to better contextualize human occupations in SW Iberia during the Early-Mid Holocene (ca. 12-6 ka BP). A review of the state of the art regarding the archaeological research in this area was done considering the known archaeological sites with a reliable chronological control, which revealed an important lack of data. Some of our most important conclusions can be summed up as follows: • The onset of the Holocene is characterized by an increase in temperature, but weak continental conditions seem to remain until ca. 10 ka BP. Between ca. 10-7 ka BP there is an increase of forest values composed of mesophilous and Mediterranean taxa reflecting a rise of moisture, after which an important decrease in the forest cover occurred followed by a stepwise rise of heathlands ca. 7 ka BP.

•
Peaks of aridity indicators were identified at 8.2 and 7. 5 ka BP. The 8.2 ka BP event seems to affect vegetation with some delay, while the 7.5 ka BP event has an immediate impact in different taxa.

•
Changes detected in the geochemical record at ca. 8.2 ka BP seem to lead to abrupt palaeoceanographic modifications, but smooth gradual terrestrial changes. The first were responsible of modifications in the current system that could have affected the coastal productivity. The 7.5 ka BP event mirror a decrease in land moisture availability.

•
Holocene sea level rise shaped the coastal morphology and influenced the settlement patterns of human groups in both Mesolithic and Early Neolithic times, as well as the preservation of archaeological sites. High energy events and the subsidence to which certain areas are subjected have hampered the preservation of any potential remain. • Subsistence strategies based on the aquatic resource exploitation were common to hunter-gatherer groups during the Mesolithic, but also during the early stages of the Neolithic.

•
Only seasonal camps were identified among the Mesolithic sites in the study area, contrary to the (semi)permanent sites in South-Central Portugal. This may be tentatively interpreted as an effect of the lack of data due to the marine transgression or different regional patterns. • Mesolithic sites are located along the coastline with direct access to the sea. Although some of them persist during the Early Neolithic, in this period most of the sites are located near rivers some kilometers inland. Several hypotheses were presented to understand this change in settlement patterns. • Changes in dietary habits and the characteristics of some shellfish species during the Late Mesolithic/Early Neolithic seem to have been related to environmental changes rather than to cultural preferences or human overexploitation.

•
There are no clear evidences of the origins of agriculture in the studied area, which may be due to an archaeological bias or a consequence of a later adoption of this practice.