Processes and evolution of the Pleistocene coastal sedimentary succession of Es Codolar (Southern Eivissa, Balearic Islands, Western Mediterranean): insights from soft-sediment deformation structures

We analyze the evolution of the undeformed Middle to Late Pleistocene deposits of Es Codolar (Southern Eivissa, Western Mediterranean). The outcrop records a succession characterized by the alternation of aeolian, colluvial and alluvial fan deposits and palaeosols that result in a complex stratigraphic architecture. In this area, aeolian beds, colluvial deposits and palaeosols are exposed along sea-cliffs for almost 500 m, allowing detailed descriptions both of the general sedimentological and geomorphological features of the Middle to Late Pleistocene deposits. Several different types of soft-sediment deformation structures are described (load-casts structures, injection structures, water-scape structures, rizoconcretions), which will help us in the understanding of the climatic evolution and the syn- post-depositional processes. In this way, main processes triggering the formation of these structures seem to be sea level changes together with a wetter environment during warmer climatic episodes.


Introduction
In the western Mediterranean, Pleistocene climatic changes resulted in several wet and dry periods (Martrat et al. 2004) combined with sea level fluctuations (Dorale et al. 2010). The past climate change record of Eivissa is recorded in the alternation of aeolianites and palaeosols, representing cold and warmer periods, respectively (Zazo 1999;Del Valle et al. 2015;Del Valle 2016).
Pleistocene sedimentary deposits of Eivissa (Balearic Islands), outcropping patchily along the coasts, provide an excellent record of glacial-interglacial cycles and eustatic fluctuations. The main representative deposits are the aeolianites, also present at other locations in the Balearic Islands (Clemmensen et al. 2001;Nielsen et al. 2004;Fornós et al. 2009Fornós et al. , 2012aPomar 2016) and western Mediterranean (Andreucci et al. 2010(Andreucci et al. , 2014Pappalardo et al. 2013;El-Asmar 1994). The quick lithification of the aeolian deposits after subaerial exposure allows preserving a high-resolution sedimentological record (Fornós et al. 2012a, b). Sedimentary successions containing Pleistocene coastal aeolianites at Eivissa use to display wide sorts of deposits characterized by colluvial facies and palaeosols and some thin shallow marine facies interbedded. In this sense, the variety of deposits and sometimes the rugged geomorphology of the coasts show a complex architecture Del Valle et al. 2016a, b;Del Valle 2016).
The main goal of this article is to present the sedimentary and climatic evolution of the Middle to Late Pleistocene deposits of Es Codolar (Southern Eivissa, Western Mediterranean). The outcrop records a succession characterized by the alternation of aeolian, colluvial and alluvial fan deposits and palaeosols that result in a tangled stratigraphic This article is part of a Topical Collection in Environmental Earth Sciences on Earth Surface Processes and Environment in a Changing World: Sustainability, Climate Change and Society, guest edited by Alberto Gomes, Horácio García, Alejandro Gomez, Helder I. Chaminé.
* Laura del Valle laura.del-valle@uib.eu; DeLa14061@studmail.ubbcluj.ro 1 architecture. In this area, aeolian beds, colluvial deposits and palaeosols are exposed along sea-cliffs for almost 500 m, allowing detailed descriptions both of the general sedimentological and geomorphological features of the Middle to Late Pleistocene deposits.
The second goal of this study is to describe for the first time the presence of soft-sediment deformation structures at Es Codolar succession and its implications for the evolution of the sedimentary succession. For this reason, identified soft-sediment deformation structures are described and deformation mechanisms are proposed.
The term soft-sediment deformation refers to several processes which change and disrupt the sedimentary structures of unconsolidated sediments. Deformation usually occurs rapidly, close to the surface, during or shortly after the deposition, and before significant diagenesis (Owen et al. 2011). The soft-sediment deformation structures are widely described in the literature (e.g., Lowe 1975Lowe , 1976Allen 1982;Mills 1983;Owen 1987Owen , 1996Maltman 1994). Deformations features in aeolian deposits have long been recognized in the geologic record (e.g., Allen 1982;Moretti 2000;Olabode 2014). They affect mostly multilayer deposits, especially sand-marl sequences, being concentrated along bedding planes or single beds, but they can as well spread through the whole deposit (Horváth et al. 2005). Most of them are related to a substantial decrease in shear resistance in water-saturated and unconsolidated sediments. They have been attributed to pore pressure changes most probably caused by fluid escape during fluidization and liquefaction (Lowe 1975;Allen 1982;Owen 1987;Rossetti 1999). Liquefaction and fluidization are the main processes which allow a temporary change from solid to liquid-like behavior in silts and sands (Allen 1982). Many natural processes may induce liquefaction and fluidization such as uneven loading, overloading, wave-originated cyclical and/or impulsive stresses, sudden changes in earthquakes and water table (Allen 1982;Owen 1987Owen , 1996Molina et al. 1998;Rossetti 1999;Moretti 2000).

Study area
The Eivissa island is situated in the south-western Mediterranean Sea and it is considered the third largest (571 km 2 ) and most Westerly island of the Balearic Archipelago (Fig. 1A). Geological structure of the island consists of a series of mainly Middle Triassic to Middle Miocene carbonate deposits thrust sheets trending NE-SW (Rangheard 1971). Muschelkalk limestone and dolomite, as well as Keuper marls and clays made up the Triassic materials of the island, whereas Jurassic limestone and dolomites and Cretaceous to Middle Miocene limestone overlie the Triassic materials and comprise the major reliefs of the island (Fig. 1B). Quaternary fluvial, alluvial and colluvial sediments fill the central basins, whereas Pleistocene coastal successions characterized by beaches to littoral aeolian and colluvial deposits crop out irregularly along the cliffed coasts (Del Valle et al. 2016a). Balearic Islands have been considered tectonically stable since the Pliocene with negligible or very little tectonic movements (Fornós et al. 2012a, b;Just et al. 2011;Sàbat et al. 2011).
The outcrop of Es Codolar (38°51′01″N-1°22′03″E) is extending continuously along 500 m of coastline and approximately 300 m inland. The Pleistocene succession is Fig. 1 A Location of the studied area; B geological map of Eivissa modified by Diaz de Neira et al. (1997); C detail of geological map of the studied area modified from Gil et al. (1997) and D topographic map and stratigraphic logs from the studied area bounded at the base by an unconformity resting on folded upper Jurassic (Malm) marly limestones (Fig. 1C). The succession shows a complex alternation of colluvial-alluvial deposits built up by heterometric matrix to clastsupported breccias, aeolianites and red palaeosols (Del Valle 2016). Es Codolar deposits are capping the steep slopes facing northwest of Puig des Falcó (144 m) and are showing a vertical cliff seaward shaped by wave erosion (Fig. 1D).

Methods
We have conducted sedimentological analyses to identify the main facies and stratigraphic units composing the sequence. Additionally, we also have performed OSL analysis on aeolianites samples to unravel the ages of these deposits and the palaeoenvironmental context that took place during the deposition.

Facies analysis
Facies analysis has been performed acquiring additional six lithostratigraphic logs (Fig. 1D). Logs were analyzed and correlated regarding major unconformities, bounding surfaces and the lateral continuity and architecture of the deposits. Units were defined in terms of brusque/abrupt facies changes.
Information on cross-bed dip direction for paleocurrent analyses was collected. Rock samples were taken for grain size, carbonate content and mineralogical composition. Carbonate content was derived by hydrochloric acid etching by means of Bascomb calcimeter method. All collected samples showed a high degree of cementation, so they were cut and polished to perform grain size analysis. Images on polished faces were taken with a 4× binocular magnifying and Motic Image 2.0 software to make observations on grain composition and size. Sample images were analyzed by means of ImageJ open source software. For the remaining analyses, samples were crushed. The resulting powder was used to identify the color by comparing dry sample powder and Munsell charts with artificial light. Mineralogical content was analyzed through X-ray powder diffraction with the SIEMENS D5000 diffractometer using Cu Kαradiation. Diffraction spectra were recorded from 3° up to 65°2θ with lapses of 0.03° every 3 s at 25 °C and compiled in digital archives for later analysis. Using X-Powder v.2010.01.09 software, a semi-quantitative analysis was carried out comparing the results with the data base Difdata (Downs and Hall-Wallace 2003).

Soft-sediment deformation structures (SSD)
Detailed photographs of the different sections displaying deformations were taken to describe the main features. The shape and geometry of the deformations were described by means of digital image reconstructions combined with the field observations (Chan and Bruhn 2014). Classifications of soft-sediment deformation structures can be used to establish morphologic and genetic systematizations. Thus, in our study area the following characteristics were identified: limited deformations among stratigraphic horizons, lateral continuity of SSD structures at appreciable distances and a confinement between non-deformed strata and its lithological association with psammitic-pelitic deposits. Softsediment deformation structures identified in this study were described and classified according the works of Rossetti (1999), Moretti (2000), Horváth et al. (2005)

Optically stimulated luminescence (OSL) dating
OSL dating relies on the properties of some mineral grains such as quartz, but not only, to store energy resulted from exposure to the environmental radiation field during their burial within deposits, and its later release it in the form of light upon stimulation with light under controlled laboratory conditions. OSL has become a key chronological tool of Quaternary siliciclastic successions and recent contributions have highlighted the suitability of this method for aeolian sediments such as carbonate aeolianites (Murray and Clemmensen 2000;Nielsen et al. 2004;Fornós et al. 2009;Anechitei-Deacu et al. 2018) and their possible correlation with the marine isotopic stages in order to obtain climatic information on the formation of these deposits and their geomorphological processes. The aeolianites collected at Es Codolar, contain between 1 and 8% of siliciclastic grains, primarily quartz. Three rock samples (~ 70 × 70 × 50 cm; 6 kg) were collected for luminescence dating from aeolian deposits at Es Codolar. Aeolianite blocks were extracted from the stratigraphic layers considered to be representative of the succession. Blocks were detached in shadowed low light conditions and wrapped in lightproof material, labeled and documented and transported to the laboratory. The inner part of each sample was extracted in the laboratory and subsequently treated with HCl (30%) and H 2 O 2 (10% and 30%) for carbonate and organic matter removal. The remaining material was sieved to obtain material < 63 µm and different fractions ranging from 63 to 250 µm. The mineral fraction was subsequently etched with 40% HF for 60 min.
The equivalent dose measurements were undertaken using a Risø TL/OSL-DA-20 luminescence reader (Thomsen et al. 2006). The 90 Sr-90 Y beta source was calibrated using gamma-irradiated calibration quartz supplied by Risø National Laboratory (Hansen et al. 2015). For these measurements, the modified single aliquot regenerative protocol (Murray and Wintle 2003) was applied, with preheat of 10 s at 220 °C, a cut heat to 180 °C and an elevated temperature bleach of 40 s at 280 °C at the end of each cycle. Radionuclide concentrations were derived by high-resolution gamma spectrometry and were transformed to dose rates using published conversion factors (Adamiec and Aitken 1998). The cosmic ray contribution was evaluate using published formulae (Prescott and Hutton 1994). Water content estimation was based on the difference in weight of the material from the inner part of the blocks, before and after oven drying, assuming a relative error of 25% for the time-averaged water content.

Stratigraphic succession and facies analyses
Textural and compositional analyses have enabled to differentiate four major sedimentary deposits which represent aeolian and colluvial facies and two types of palaeosols (Fig. 2).

Aeolian facies
Aeolian facies are characterized by very pale brown color (HUE 10 YR 8/2), well-sorted, fine to medium-grained sandstones, with large-scale trough cross stratification. The strata are 0.5 m-3 m thick and are moderately disrupted by root casts (0.50 a 5 cm width and 0.10-1 m in length). Most of them are remains of moulds and calcified root casts. The facies composition is mainly medium-grained carbonate marine sand with very little terrigenous material (i.e., quartz grains). As stated to the diffraction analysis, calcite is the main mineral (80%) and there is a minor presence of clay minerals. These deposits record the wind-transported marine carbonate sand trapping in front of a steep inland cliff. The aeolian levels present an active migration inland towards ESE-SE direction. High-angle cross-stratified layers have a direction of N45°/25°/SE.

Colluvial facies
Colluvial facies are characterized by massive reddish (HUE 5 YR 6/8) to pale brown (HUE 10 YR 8/3) silty matrixsupported breccias, composed of angular clasts forming millimeter to centimeter thick layers disrupted by calcrete levels. Large clasts are aligned downslope (270°-330°NW), although the clasts display modal fabric orientation. Clasts come from the erosion of the basement composed of Jurassic marly limestones. Conforming to the diffraction analysis, the main composition is characterized by carbonates (72%). These deposits record episodic mass movements on slope surfaces characterized mainly by debris-flows transport.

Palaeosols
Two different types of palaeosols have been observed. The first one, around 20-40 cm thick is mainly composed of reddish silts (HUE 5 YR 6/8) with the presence of iron bands. According to the diffraction analysis, quartz is the main mineral (20%) and there is a minor presence of clay minerals. The second one is characterized by massive silt and silty-sandy deposits with nodular forms and lots of bioturbation evidence with some root casts. Its color is very pale brown (HUE 10 YR 8/4) and layers are 0.30-1.5 m thick, rich in carbonate (71%). These palaeosols may represent periods of warmer environments characterized by very little erosion rates on slopes and the development of vegetation.

Stratigraphic units
Es Codolar outcrop is composed of four units. Unit 1 (U1) consists of three colluvial levels interbedded with palaeosols. This unit is bounded at the base by an unconformity resting on the Jurassic basement. Thin calcretes are capping the colluvial-palaeosols sequence at the upper part.
In general, it is possible to observe three repeated sequences at U2, U3 and U4. All of them are made up from the base to the top by calcretes, palaeosol and/or colluvial deposits, aeolian deposits and palaeosol. The lower palaeosol of each unit is filling up the rizoconcretion moulds of the immediately upper sandy level (aeolian facies) (Fig. 2).

Soft-sediment deformation structures
Eleven different types of well-exposed soft-sediment deformation structures were distinguished in the Es Codolar. On characterizing the deformation structures, we considered lithology, size and geometrical characteristics (Fig. 3). Fifteen deformed levels have been observed .

Contorted structures
Contorted structure is used as a general term for all types of features with different degrees of crumpling or complicated folding of the laminae within sand beds (Brenchely and Newall 1977). It was possible to define several types of contorted stratification: small-scale slump, large-scale slump, small-scale load structures, large-scale load structures, deformed lamination, recumbent folds, convolute folds, ball and pillow, concave up path with consolidation lamination. In the study area, three deformation structures were observed: load structures and recumbent folds.

Load structures
Load structures are laterally continuous undulations at an interface between layers that varied in density with relatively denser sediment above and less dense below. Deformation is driven by the gravitational forces associated with the reverse density system and it is dominated by vertical displacements (Owen et al. 2011). Internal laminae are well-recognizable and, typically simulate the external circular form of the load structures (Moretti and Sabato 2007).
These forms are located in between carbonated sandstone and silts. The undulations are between 30 and 50 cm width, Fig. 3 Classification of the soft-sediment deformation structures described at Es Codolar succession. Different driving force systems, deformation and trigger mechanisms for each kind of soft-sediment deformation structure are presented and 20 cm and 30 cm high. They display a regular concave and convex morphology (wavy laminae) along the outcrop (Fig. 4). This deformation is present in all units of the outcrop within the sandstone levels, silty layers and calcretes levels. These structures extend continuously for 200 m along the coastline. The presence of deformation breccia (Fb), flexi caves (fc), convex structures with vertical axial plane characteristic by 20 cm in height and 13 cm in width is associated with these structures (Figs. 5 and 6).

Recumbent folds
These structures consist of cross-sets with overturned internal stratification, which forms recumbent folds with either horizontal or slightly inclined axial planes (Rossetti 1999). These structures were developed in cross-bedded calcretes or loamy sands. Layers under the recumbent folds show sigmoid and concave forms. Recumbent folds at the uppermost part have parabolic morphology and there is no structure on the fold limbs. Recumbent folds are 15-20 cm in height and 30 cm wide. Fold axes are sub-horizontal and generally oblique to the dip of the not deformed foresets strata. The folds have mainly rounded hinge zones. These structures can be observed at the level C2 of the Unit U1 (Fig. 5).

Injection structures
The injection structures are composed of root cast moulds with a size of 0.50-5 cm width and 15-25 cm height in sandstone levels. They are filled by reddish silts from the underlying palaeosol. The extrusive flows are low marked (ef). The base and upper contact are irregular (Figs. 5 and  6). The analysis of the mineralogical composition carried out both on the lower palaeosol and on the materials that fill the voids-rizoconcretions moulds-of the aeolianites coincide. Conversely, there is no coincidence with the palaeosol located overlying the sandstone layer. These structures occur on two levels; in D1 (Unit U2) and level D2 (Unit U3). Nevertheless, we observed other injection structures within several levels on Cl3 (Unit U1), P4 (Unit U2) and P5 (Unit U3).

Water-escape structures
The term water-escape structure-hydroplastic mixing layer-is here used in a general sense to include all layers of sediment which have been internally deformed by hydroplastic flowage but which have not moved as a whole relative to surrounding sediments. Internally, liquefaction layers and pockets may show severely deformed primary structures ( Fig. 7), nearly complete homogenization, or water-escape structures formed by the redistribution of grains during water-escape and re-sedimentation (Lowe 1975;Trincardi et al. 2004). They are formed when sediment is transported by the expulsion of pore water or gas (Lowe 1975;Owen et al. 2011).
The resulting structures show a range of 15-20 cm in height and 3-25 cm in width. Deformation is limited at the Fig. 7 A Water-escape structure developed in a complete deformed sandstone primary stratified bed; extrusive flows (ef), distorted lamination (dl), folded dykes (dk). B Empty volcanoes/bell injection deformed structure developed in a complete deformed sandstone bed.
Note the high presence of bioturbation with abundant root casts (0.1-1.5 cm in width and 0.10 to 15-20 cm weight). Lithological characteristics, sandstone (Sd), loamy sands (Sc) and marl (Cs) top by irregular palaeosol that shows extrusive flows. They are mainly represented by vertical conduits. The vertical development (perpendicular to the primary lamination) is very variable from linear to very irregular. The upper part of the conduits is often represented by conical morphology (wider downward)-like volcanoes-that was probably located at the initial water sediment interface and they end upward with "vent" level. In Fig. 7A, water-escape structures in the sandstone beds forming diapirs or sand volcanoes can be observed. On the other hand, in Fig. 7B, the same structure but completely empty can be observed. These structures are observed in the sandstone levels D2 (Unit U3) and D3 (Unit U4), respectively.

Flames and pillars
Flame structures (Figs. 5 and 8) developed at the interface between sandstone and silts are associated with large-scale load structures (KoçTaşgin et al. 2011). They consist of short, irregular and curved upward protuberances located along sharp boundaries formed between two deformed layers. Flame structures result from diapiric intrusion of fine sediments (Mills, 1983); they formed due to a gravitationally unstable density gradient.
Pillars (Figs. 5 and 8) consist of discrete, straight to slightly sinuous, steeply inclined to vertical paths up to a few centimeters length (usually < 0.10 m) and less than 0.03 m in diameter, which cut sharply through the deformed beds (Rossetti 1999).

Dyke and sill dyke
This structure consists of masses of irregular, convolute laminated or massive sandstone that intruded sharply upwards into the overlying strata to form large (up to 0.65 m long and 0.29 m wide) elongated shapes. Although dykes are morphology similar to pillars, the latter term is reserved to describe discrete, straight to sinuous intrusions, while dykes are larger and show more irregular shapes (Rossetti 1999). The dykes in the study area vary from vertical to steeply dipping and their basal edges attached to underlying sediment sources (Rossetti 1999;Moretti and Sabato 2007).

Ice/sand-wedges
Ice/sand-wedges form by infilling of thermal constriction cracks in permafrost with water or sublimated ice crystal.
A key feature of ice-wedge is the occurrence of vertical to steeply dipping laminae generally 0.5-5.0 mm thick (Mackay 1986 in: Murton and Ballantyne 2017;Wolfe et al. 2018). Cross-cutting relationships among laminae are common. Each lamina represents an ice vein formed by the filling of a single thermal contraction crack with ice. Thermal contraction cracks may fill with snow, ice, mineral particles or organic material, producing a variety of vein and wedge-shaped structures characterized as icewedges, sand-wedges, soil wedges and composite (ice and sediment) wedges (Murton and Ballantyne 2017). These structures have been observed at the levels P4 and P6; characterized by regular size (1-3 cm height and 10-14 cm width) (Fig. 8).

Collapse structures
The slope collapse in most cases occurs when an area has a steep slope and in the sedimentary succession some lower levels are saturated. The horizontality strata are affected, showing the movement of the upper materials due to the plastic behavior of the lower levels, marked by a rupture of the original slope-rugged cut slump. The lower levels are expected to transport the upper ones without affecting the original succession of layers. At the delay of movement derived from the compressive effect of the displacement front, fluids extrude forming cones or cut dikes (Fig. 9). Distinct deformations of this type initiate the movement as result of differences in the hydrostatic gradient (Strachan 2002). Normally, the occurrence of slumps requires the presence of a steep slope (KoçTaşgin et al. 2011).
Others collapse structures are formed when sediment can no longer be supported by its substrate. This may occur in association with karstic solution or substrate weakening (DeMille et al. 1964). The resulting structures can be very diverse (Fig. 10), although they can be overlapped with load structures (Owen et al. 2011).

Cryogenic brecciation
These forms consist of numerous irregular rocks, isolated blocks, or sometimes oriented blocks (Figs. 5, 6 and 10). The blocks have an irregular size generally with 1-5 cm short axis and a 5-15 cm long axis. They are produced by cryotoplastic weathering (Lanfranco et al. 2014).

Others associated structures
Depositional structures may be disrupted and distinctive new structures may form as a result of chemical, physical, (e.g., crystal or concretion) and biological processes (e.g., plants roots, footprints) (Fig. 11).

OSL ages
Three OSL ages have been obtained in the Pleistocene succession of Es Codolar. The measured parameters, relevant for age calculation, are presented in Table 1. The samples were collected in aeolian facies beds in order to establish which the main aeolian events that took place here (Fig. 2). The first aeolian level dated corresponds to U2 and shows an age of 317 ± 20 ka, coinciding with the cold stage MIS 10. The next dated sample was from U3 which gave an age of 227 ± 15 ka, coinciding with the cold stage MIS 8. The last one dated was U4; this coincides with the end of the penultimate glacial stage MIS 6 with an age of 124 ± 9 ka (Figs. 12).

Fig. 10
Picture and line draw of the collapse structure (Cc) associated with the large-scale load structures presence in the outcrops of Es Codolar (Lc); injection structure (Is); slope collapse (rcs), bioturbation (bt), folded breccia (Fb) and dyke (dk); cryogenic brecciation (Cbt) Fig. 11 A Detail picture of the pisolits level. B Calcretes of Es Codolar outcrops. C Roots casts and root casts molds (5-15 cm in width and 50 cm-1 m in height), the x represents two levels with highly small root casts (0.1-2.5 cm in width and 1-10 cm in height approximately) each moment, that is, for example, the presence of aeolian sediments during a more arid moment contributes to the formation of the silty-sandy soil. Main source of the silty material is arid zones of Northern Africa affected by winds which are able to carry fine particles towards the Balearic Islands and form the majority of the soils (Fornós et al. 1997;Fiol et al. 2005;Muhs 2007).
According to the OSL dating from the aeolian levels, the deposition of coastal dunes of units U2, U3 and U4 took place during the MIS 10, MIS 8 and MIS 6 respectively. This leads us to infer that the colluvial deposits and palaeosols of units U1-U4 were formed during the interglacial stages MIS 11, MIS 9, MIS 7 and MIS 5, respectively.
The study of soft-sediment deformation structures helps us to understand the sedimentary and climatic evolution of Es Codolar sequence. In this sequence, the presence in loamy sandstones beds of ice-sand-wedges at levels P4 and P6 (Fig. 2), located just above the aeolian levels D2 and D3, confirms the occurrence of several periods of cold temperatures just after the deposition, during the glacial stage, of D2 and D3 aeolian levels. But these periods of cold temperatures occurred even during the interglacial stages, as shown by the presence of cryogenic breccias at levels P2 and P3 (located below the D1 aeolian level), deposited during the interglacial stage MIS 11.
On the other hand, in Es Codolar sequence, the softsediment deformation structures are considerably larger and frequent into the lower levels than into the upper levels, pointing to an increasing load as responsible for the development of these structures. In this way, structures formed after the deposition of the majority of the sedimentary levels. Moreover, because the liquefaction structures formed in water-saturated sediments, it is more probable than these structures were formed during the interglacial stages, when sea level is higher, the rainfall increase and therefore, there is a greater presence of vegetation. In Fig. 1, a coastal lagoon is observed near the outcrop, indicating the significant proximity of the water table. Thus, in Es Codolar, during the interglacial stages the sediments were saturated due to both a wetter environment and a higher water table, being the maximum high during MIS 5, when sea level was nearly 3 m higher than today (Polyak et al. 2018). Because of these two reasons (1st-maximum sea level during MIS 5, and 2nd-deformation structures larger and more frequent at the lower levels), we infer that the soft-sediment deformation structures were formed at the end of MIS 6 to MIS 5, after the deposition of aeolian unit U4.
Based on the discussion above, Fig. 12 shows the evolution of the Pleistocene sedimentary sequence of Es Codolar. After the formation of the colluvial beds at the base of the sequence, the first aeolian level was deposited during a relative sea level drop in MIS 10 (Fig. 13). Bioturbation structures developed in this level, with dominant palaeo-wind currents from SSW to NNE (Fig. 12A). During MIS 9, a thin sequence (1 m) of colluvial and palaeosol sediments was deposited on top of the previous aeolian level. A second aeolian level was formed during MIS 8 under a predominant wind from S to N. Bioturbation structures also developed into this level (Fig. 12c). During MIS 7, characterized by a higher water level and humidity than in MIS 6 ( Fig. 13), a thin sandy palaeosol was formed, with some intercalation of colluvial beds (Fig. 12d). The upper aeolian level was emplaced during MIS 6 and has the same characteristics as the previous aeolian beds. At the end of this period soft sedimentary deformation structures formed. These structures are larger and frequent in lower levels because of the increasing load as we go lower in the stratigraphic series. The sequence of Es Codolar ends with the deposition of a thin sequence of sandy and loamy palaeosols and calcretes corresponding to MIS 5. This combined with an intensification in the pluviometry and the rise of the water table, produced fluidization structures of these palaeosol into the upper aeolian level.

Conclusions
Results of sedimentary analyses and OSL dating allowed distinguishing the effects of climatic fluctuations during the sedimentation of these Pleistocene deposits. Due to the strong climatic signal registered in Es Codolar sequence and the absence of significant tectonic movement, changes in the environmental agents (i.e., arid-wetter periods and sea level changes) are the most plausible processes controlling the described sediment deformations. The common factor in the formation of soft-sediment structures is that all the origin sediments went through liquefaction-fluidization, which appears to be mostly triggered by a combination of sea level rise and increasing humidity during interglacial periods. It should not be forgotten, however, that wave erosion at the hillslope could also contribute to instability of the outcrop. On the other hand, during glacial periods, the encroaching of new dune fields on Es Codolar added up thick sedimentary bodies to the outcrop package that could have increased the overloading processes. The colder conditions during these periods could have boosted the cryogenic-related structures, which help us in the understanding of the climatic evolution and the syn-and post-depositional processes. The main processes triggering the formation of these structures seem to be sea level changes together with a wetter environment during warmer climatic episodes.
Funding Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.

Fig. 13
Sea level reconstructions for the last 500 Ka based on continuous records of O 18 benthic by Silva et al. (2009); sea level reconstructions for the last 500 Ka based on continuous records, compiled from the Red Sea by Rohling et al. (2009Rohling et al. ( , 2012; O 18 benthonic registers present in the deep waters of Waelbroeck et al., (2002); coral reefs register for Camoin and Webster (2015) and representation of the periods of aeolian sedimentation on the island of Ibiza based on the absolute OSL dating performed here, shaded in pale brown. The black line symbols represent the OSL age ± error