Tordera River Delta system build up ( NE Iberian Peninsula ) : sedimentary sequences and offshore correlation

The Tordera River Delta is made up of coarse sand bodies which were formed during the last sea level rise. Marine and land seismic profiles and wells show three depositional Holocene sequences corresponding to a) the transgressive base level, b) the progradant intermediate unit and c) the progradant unit which onlaps the Plio-Pleistocene sequence seawards and the granitic basement on the present delta plain and northern coast. The morphology of the sedimentary bodies, isobath maps and marine geophysical data enabled us to gain fresh insights into the evolution of the delta and its prodelta and into the various relict sand bodies on the inner shelf. The evolution of the coastline can be followed during the prograding process of the lower units from the successive coastal and internal lagoons. These coastal lagoons were conditioned by river changes and basement morphology. Once the sea level maximum was attained, the prograding units covered these morphologies, resulting in the present delta configuration.


INTRODUCTION
The geometry and evolution of a delta are the consequence of the fluvial regime, the sediment grain size, the tectonics of the area and the oceanic system of the receiving basin (Liu et al. 2000, Pellegrini 2015).In addition to the classification of deltas as a function of the dominant processes in the delta front area (Wright and Coleman 1972, Galloway 1975), two important factors should be considered to account for delta dynamics: sediment supply and grain size.These factors influence the gradient and channel pattern of the delta plain, the type of shoreline and the response to wave energy and tidal regime, in addition to the resedimentation processes on the subaqueous delta slope (Macpherson et al. 1987, Postma 1990, Orton and Reading 1993).Consequently, the deltas formed under the environmental conditions of a sediment supply and a depositional gradient seawards yield a different depositional architecture in terms of grain size.Coarse sand deltas are characterized by a delta plain with two types of channel patterns: very high bed load channels (gravels) and bed load channels (coarse sand).When the very high bed load channels reach the coast, Gilbert-Type foresets are formed if the basin is sufficiently deep (Orton and Reading 1993).This is the case of the coarse-grained fan deltas, which have a mainly alluvial, lacustrine or glacial origin (Dunne and Hepton 1984, Corner et al. 1990, Martini 1990) or are subjected to significant tectonic influence (Bardaji et al. 1990).The bed load channels in the outflow tend to be concentrated in a network of bifurcating multiple channels separated by numerous bars (Orton and Reading 1993).This is the case on the Losna and Antnsjo sandy braid deltas (Bogen 1983).Along the Mediterranean coast, the coarse sand deltas described consist of gravel pleistocene-Holocene fan deltas influenced by recent tectonic and wave activity (Massari and Parea 1990, Dabrio et al. 1991, Pellegrini et al. 2015).Only two examples of sandy braid deltas have been cited: Goynuk Delta in the gulf of Antalya (Ayranci et al. 2010) and the Tordera delta, analysed here.
At regional level, the stratigraphic and sequential Holocene evolution of small delta systems of the Rivers Fluvià, Besós, Llobregat and Foix has been studied by several authors (Checa et al. 1988, Díaz and Ercilla 1993, Gámez et al. 2009, Velasco et al 2012).These deltas are similar in structure, development and sedimentary processes, with the exception of the coarsegrained Tordera and Besós deltas, both controlled by glacio-eustatic sea level changes, fault activity, and the action of waves and drift currents (Serra et al. 1985, Riba and Colombo 2009, Velasco et al. 2012).The Holocene sediments of the Besós delta are only represented by the upper units owing to the presence of a Plio-Pleistocene high which prevented the sedimentation of the transgressive lower units.
In the absence of studies about coarse sand deltas in the Mediterranean, the Tordera delta is a fine example of this delta type in microtidal conditions dominated by fluvial processes and littoral drift.The stratigraphic reconstruction of this specific whole delta system also offers new insights into its build-up and its chronology.Moreover, the permeability characteristics of coarse sand deltas make them a good geological model as possible oil reservoir.Therefore, given the scarcity of this type of deltaic bodies, the Tordera River is a unique case for studying this kind of formation in the sedimentary record.
Data integration of delta and submerged prodelta geophysical profiles and the detailed sequential study of on-land drills allow us to establish the correlation between units and sequences of the two domains.A chronology of the sedimentary sequences is proposed on the basis of these findings, and geographical distribution of the sedimentary units in addition to a paleogeographic reconstruction are also presented.The model distinguishes three sedimentary sequences: the lower, middle and present Holocene unconformably overlying the Plio-Pleistocene shelf materials and the granitic basement.

GENERAL SETTING
The River Tordera flows into the sea between Malgrat de Mar and Blanes in the NE of the Iberian Peninsula.The upper and middle sectors of the River Tordera follow the trench between the Littoral and Pre-Littoral Cordilleras before veering to the SE and then crossing the coastal range (Serra 1976).The Tordera has a basin area of 879 km 2 and a length of 59 km, and its delta aerial surface is 4.2 km 2 .The river has a typical Mediterranean regime with 7.2 m 3 s -1 mean water discharge, long dry summers and sporadic floods and spring time (Martín-Vide and Llasat 2000, Condal and Ullastres 1996) as registered at Can Serra of up to 1280 m 3 s -1 (ACA 2002) and also at Fogars de Tordera of up to 1836 m 3 s -1 (Liquete 2008).The delta is mainly composed of coarse sand derived from granitoid and metamorphic rocks of the pre-littoral mountains.Because of the coarse texture of the sand, the Tordera delta is characterized by a high permeability that facilitates water infiltration into the delta plain, by the presence of a well-developed coarse sand beach mouth bar (most of the time between floods), and by a high gradient delta front (up to 23°).Sediment supply is carried mainly during high regimes as bed load (>75%) and only less than 20% as suspended load (Rovira et al. 2005).The total contribution of sediment to the sea has been calculated by a number of authors using different methods.Considering the littoral drift necessary to feed the relict sedimentary bodies on the continental shelf, the sediment carried can attain 200000 m 3 a -1 (Copeiro del Villar 1982, Serra et al. 2007).According to Rovira et al. (2005) the river bed load sediment transport amounted to 21500 m 3 a -1 .
The Tordera fluvial-deltaic system resulted from the pulsations of the Holocene sea level rise (between 11 and 6.6 Ka) that led to successive deltaic sequences culminating in the current Tordera delta.These sequences range from a depth of 60 m to the present delta front, which currently progrades over the basal transgressive sequences.Geophysical surveys (Teixidó 2000, Falgàs 2007) and deep wells for the Tordera desalination plant (Serra et al. 2007) have considerably improved our understanding of the structure and evolution of the Tordera delta system.The morpho-dynamic and structural study of the submerged continental shelf has been undertaken by several authors (e.g.Serra 1976, Díaz and Maldonado 1989, Duran et al. 2013, 2014), who have contributed to the knowledge of sedimentary processes and the eustatic evolution that affected the submerged continental shelf linked to the present delta formation.
The River Tordera gives rise to the Maresme sedimentary coastal cell (Fig. 1), which extends to the River Besos, 50 km to the SW.The sediments that reach the mouth of the Tordera are redistributed mostly by the dominant wave climate to the southwest or descend the steep delta talus.The impact of other dynamic factors, such as geostrophic currents and the action of the tides, are negligible (Font et al. 1987, Sorribas et al. 1993, Grobas 2003).In the case of high energy events, part of these sediments may be transported to the internal shelf and Blanes Canyon head zone (Zúñiga et al. 2009, Pedrosa-Pàmies et al. 2013).The effect of the dense water currents (DSW cascading) is also low on the prodelta and surrounding shelf (Lastras and Euroleon 2011).
The wave climate is dominated by a NE to E regime, which shaped the coast and the existing relict sand bodies on the shelf (Díaz andMaldonado 1989, Serra et al. 1990).This regime accounts for the significant accumulation of relict sands during the last low sea level along the Maresme shelf to the SW of the River Tordera.Towards the Blanes Canyon head, only two main sand bars are distinguished (Fig. 1).A transverse sediment transport component must be taken into account given the high wave reflectivity of this coast.

METHODS
Sequential stratigraphy was used to better understand the evolution of the Tordera River Delta and the sedimentary bodies on the continental shelf, following the approach developed by Catuneanu et al. (2009).The findings offer new insights into the facies distribution, the geometry of the strata, and the major surfaces, establishing a chronological order of the basin infill.To this end, we integrated our data with other diverse data sources in order to identify depositional trends and sedimentological patterns of this environment.
The following types of data were used in this study: land wells, marine and land seismic profiles, bathymetry (single and multibeam) and side scan sonar.Side scan sonar is a specialized system for investigating the seafloor morphology: it transmits sound energy and analyses the return signal (echo) that has bounced off the seafloor.The transmitted energy is shaped into a fan that sweeps the seafloor from directly under the towfish to either side.The strength of the return echo is continuously recorded, creating an image from the ocean botton (background structures) that is comparable to aerial photography on land (Able et al. 1987).

Delta plain
The internal structure and the stratigraphy of the Tordera delta were analysed using two high-resolution seismic profiles (Fig. 1, Teixidó 2000, Falgàs 2007) and a network of 50 land wells (Fig. 2) (ACA 2002, Sufi-ACA 2008, Perfibesa-ACA 2001-2002).Seismic data were acquired using a 48-channel digital seismograph, 40-Hz geophones and 5-m shot and receiver spacing (Texidó 2000, Falgàs 2007).Out a total of 77 land wells, 50 were selected for this study, considering depth and number of identified units.The land wells show alternations of fine sediments, sands, and gravels, and are described in the results section.New textural and sedimentological arguments such as fine sediment content were considered in the stratrigraphic interpretation of the wells.This interpretation sheds light on the paleogeographic Holocene evolution of the delta, enabling us to obtain a 3D view of the structural and environmental distribution of the delta.

Prodelta and shelf
The geomorphology of the Tordera prodelta and internal shelf was determined by integrating bathymetric data from the earlier projects ESPACE and Beachmed-  Lastras and Euroleon (2011).In red and orange, respectively, situation of seismic profiles P (Serra et al. 2007) and PL (Díaz and Maldonado 1989).In green and yellow the aerial seismic profiles from Teixidó (2000) and Falgàs (2007) respectively.
GESA-UB and present bathymetric works from UB-CEAB using a SIMRAD EM3000D, a SeaBeam 1185 echo sounder with a frequency of 80-125 khz and a Garmin 76CSx GPS map (Fig. 1).

RESULTS
The sedimentary units are described in this section with respect to their present situation onshore and offshore.The methods for land and sea domains are well data and seismic profiles, respectively.

Delta plain
Three stratigraphic sequences were differentiated from base to top: Sd1, Sd2 and Sd3.Sequence Sd1 lies unconformably over a pre-Holocene substrate or over the granite basement.This basement is especially significant since it plays a major role in the transgressive units and contributes to our understanding of the evolution of the present delta.There are two structural highs of the basement on the delta front and northeast delta plain sector (Perea 2006), as shown on seismic profiles from Teixidó (2000) and Falgàs (2007), well data units and basement distribution.These highs are limited by faults and exert a paleogeografic control over the deposition and distribution of the lower holocene units.
The sequences are coarsening upwards and have two sedimentary units for Sd1 and three for Sd2 and Sd3 (Fig. 3).The limits between the sequences were deduced from textural and geophysical changes using sedimentological data such as the net contact between the upper coarse detrital sequence material and clay or lime at the base of the upper sequence.Although these contact surfaces must be regarded as heterochronous (transgressive or prograding), they also serve as a reference level for correlation purposes.
In the frontal area of the delta, sequence Sd1 is made up of two units, Sd1ll and Sd1sup from bottom to top (Fig. 3).Sd1ll consists of grey silt in the southern hemidelta and silt and grey clay in the northern one.In the internal area of the delta plain, the unit changes laterally to sandier facies (Fig. 2, S-21 without unit Sd1ll and ACA 4).
In the southern hemidelta, Sd1sup shows alternating coarse detrital levels, some of them with bivalve fragments, levels of grey silty clays and levels of silt near the present delta front.Towards the internal delta, this unit consists of silt, sand and silty sand (well S-21, ACA 4, Fig. 2).In the northern hemidelta, which is made up of coarse detrital materials, a lateral facies change can be deduced: the wells closer to the mouth have coarser lithologies (gravel and pebbles P1-S1 well, Fig. 2), whereas those further away show alternating levels of sandy gravel with or without interstitial clay (P6-S5 well, Fig. 2).Clay or silt levels are absent.Granitic basement was reached at depths of between 70 and more than 150 m in the surrounding wells in the front of the delta plain (wells S-21 and P6-S5, respectively).
The Sd2 sequence (Fig. 3) is composed of three units: Sd2ll, Sd2p and Sd2sup from bottom to top.Sd2ll is located in the front part of the delta and consists of fine to very fine sand and grey silt.Towards the delta plain, Sd2ll changes sideways to sandier lithologies, as at well S-21.Sd2p consists of medium sands, gravel, clay and brown silt.This unit is found throughout the delta plain, except in the northeastern part because of the presence of a structural high of the basement at 20-24 m (well B-3-A, Fig. 2).Sd2sup consists of coarse sand, gravel  and silt, displaying coarsening-upwards sequences.Sd-2sup is the first unit present in the entire Holocene delta (Fig. 2).From the correlation between nearby wells, it can be deduced that these materials were deposited in a fluvial environment, as evidenced by their channel and overflow facies and geometry (Fig. 4).Sequence Sd3 (Fig. 3) comprises three units from bottom to top: Sd3ll, Sd3p and Sd3sup.Sd3ll occurs in two sectors of the delta: at the front of the delta plain (subunit Sd3ll1) and on the northern delta plain (subunit Sd3ll2).In the southern hemidelta, Sd3ll1 consists of grey-black clays with organic matter and vegetal remains.In the northern hemidelta, this unit is more variable with alternating sandy grey clay, silt and bioclasts, bivalves, brachiopods and gastropods (well P11-Supr.).Subunit Sd3ll2 is made up of dark-grey silt with abundant organic matter, but with no marine bioclasts (well B-3-A).Subunits Sd3ll1 and Sd3ll2 change to coarse sandy lithologies towards the delta plain.Unit Sd3p is coarser than unit Sd3ll1, with levels of fine and medium sands, gravel, silt and clay deposits.The Sd3sup unit is mainly detrital in the upper delta plain, displaying a lateral change to fine and well-classified sands towards the internal plain.

Prodelta and shelf
The analysis of the high-resolution seismic profiles (Fig. 1 for location) allows us to differentiate four seis-mic sequences: a lower sequence (PH) and the upper sequences Sd1, Sd2 and Sd3, each one formed by internal units.These sequences are defined by their erosive (top of PH) or unconformable limits (Sd1-Sd2 and Sd2-Sd3) and internal reflection configuration (Fig. 5).Basal sequence PH is characterized by sigmoidal morphology reflectors seaward and onlaps the surface towards the coast in contact with the granitic basement.Unit PH is eroded in its upper part, showing an irregular topographic surface.The limits between upper sequences (Sd1-Sd2 and Sd2-Sd3) are unconformable (Fig. 5).Internally, Sd1 displays a basal onlap and downlap relationship and aggrading to sigmoidal reflections (Fig. 6).It is interpreted as an accretionary system tract (TST) based on its internal configuration and reflection terminations.
The lower boundary of Sd2 corresponds to a discontinuity defined by an onlap surface towards the coast that evolves seawards to a downlap surface.Sd2 comprises three internal units (units 2 to 4, from older to younger) that differ in development and distribution (Fig. 5).Unit 2 is formed by clinoforms with oblique progradational and aggrading configurations (Fig. 6).Unit 3 is a wedge-shaped seismic unit that sits basinward downlapping onto unit 1, and it is interpreted as a falling stage system tract (FSST).Unit 4 shows aggrading and onlap reflectors from the coast to downlap on the Sd1 sequence seawards, where the morphology is a sigmoid-oblique.This prograding unit is interpreted as a low-stand system tract (LST).Sd3 consisting of only unit 5, displays a wedge-shaped geometry internally defined by seawards downlapping clinoforms, suggesting that it was formed in a late stage of sea level rise (HST), resulting in progradation and aggradation of coastal sand bodies.This unit corresponds to the top lap delta unit.
These seismic sequences (see Fig. 5 for system tract order) constitute the progressive transition from an HST (corresponding to unit 2), an FSST (unit 3) and an LST (unit 4) to the present HST (unit 5).It should be noted that, in profile P-2, Sd1 corresponds to a transgressive level (TST-1).Transgressive levels are not present in sequences Sd2 and Sd3 because of one of the following factors: a) the erosion of earlier transgressive levels, b) a thickness below the limit of seismic resolution, and c) a sea level rise that does not exceed the sediment supply, resulting in a continuous progradation.The base of the TST (Sd1) and the lower limits of the HST are used to make the correlation with the sedimentary sequences described on the onshore wells, and to make the correlation with the different sand bodies on the offshore (Fig. 7).

Units, evolution and chronology
As regards the delta systems in the NW Mediterranean, it may be affirmed that the eustatic and climatic processes dating from the postglacial transgressive level are common to all deltas.Some parallels may be drawn from the deltas of the Rivers Llobregat, Besós, Tordera and Rhone, (Gámez et al. 2009, Velasco et al. 2012, Boyer et al. 2006, in Fig. 8).
The boundaries between PH-Sd1, Sd1-Sd2, Sd2-Sd3 and Sd3ll-Sd3p are equivalent to the transgressive surface and stabilization depths defined at 70, 60, 38 and 22 m by the above authors, with slight differences in the chronology of units and levels.
In order to determine the differences and similarities between the aforementioned delta examples, the morphoclimatic, sedimentological and dynamic factors of the respective fluvial basins should be considered (Table 1).
This table highlights two prominent features: a) the volume of the river discharge (average and maximum), which accounts for the transport capacity and valley excavation-aggradation rate of the river; and b) the lithology of the fluvial basin, which determines the typology of the river sediment.This is the case of the Besós River Delta, where the transgressive surface at the delta mouth is only 30 m deep, whereas it is 60-70 m in the Llobregat, Tordera and Rhone deltas.This significant difference is attributed to the low river discharge and to the low capacity of excavation, which gives rise to the presence of a Plio-Pleistocene delta basement and is absent in the above examples (Table 1).The second feature accounts for the river sediment characteristics such as fine content and bed load transport, resulting in the prodelta sedimentation and in the presence of a delta aquitard.This is the case of the Llobregat, Besòs and Rhone deltas, with the exception of the Tordera coarse sand delta, where there is no formal aquitard.The second characteristic also accounts for the differences in unit chronology and depth because of the different aquitard structure and the subsidence rate as a function of this structure.
The correlation between land and marine sequences was carried out using the nearest wells on the coast (wells of Sufi-ACA 2008, Perfibesa-ACA 2001-2002) and the seismic data offshore.The boundaries between sedimentary sequences of the Tordera delta and the inner shelf are represented in a graphical model (Fig. 9).The model sections show the evolution of the Tordera Fig. 8. -Comparison of other deltas in the NW Mediterranean with the Tordera delta.Four surfaces with a regional extension are observed: Mfs (in fuchsia and yellow), the top of transgressives units (in green), the top of Sd2 sequence (in orange) and the transgressive surface.
delta system and some aspects of the location of the sedimentary units: e.g.pelitic units Sd3ll1, Sd2ll and Sd1ll observed on the present delta front.The sea level rise did not exceed the height of the granite basement, with the result that the Sd3ll1, Sd2ll and Sd1ll pelitic levels were subsequently covered by the upper detrital units (Sd3p, Sd2p, Sd1 sup) sedimented in the subsequent prograding and aggradational deltaic environments.Deltaic sequences Sd1, Sd2 and Sd3 show continuity seawards with the present relict sedimentary bodies.Four sectors may be differentiated on the basis of the morphological characteristics: A) the Tordera prodelta and Malgrat platform sector, B) the rocky shore littoral prism from Lloret de Mar to Tossa, C) the inner shelf between the Tordera prodelta and the Blanes Canyon, and D) the Blanes Canyon head (Fig. 10).
The analysis of bathymetric, Side Scan Sonar and high-resolution seismic profiles revealed the submerged morphological features, structure and distribution of the sedimentary units in the prodelta and shelf as far as the Blanes Canyon.The inner continental shelf of The Malgrat sector (Fig. 10) at the base of the steep delta front displays alobeted morphology body, A1, which descends gently to a depth of 60 m.A2 corresponds to the present delta front and the high gradient slope characterized by the presence of gravitational sand avalanches and traces of grain flows, as shown in side scan sonar profiles (Serra andMontori 2003, Parra 2006).
Two detrital bodies, C2 and C1, extend to the E-NE between the steep east delta front and the head area of the Blanes Canyon, descending respectively to depths of 35 and 52 m, as shown in profile P1 (Fig. 7).The front of these two bodies is similar to the delta front and the Malgrat ridge bar or spit (Serra et al. 2007), attaining a gradient of up to 23° and a drop of 7 and 20 m, respectively.This sector also displays outcrops of a pre-Holocene substrate that could have conditioned the morphology and the size of these detrital sand bodies (Serra 1976).The sedimentary bodies B2 and B1 (Fig. 10) correspond to the present and to the relict high gradient littoral prisms, respectively.These bodies consist of detrital input from the streams of the locality (Tossa -Blanes).
The morphology of the transgressive clastic units is consistent with the marine dynamic regime that existed during the last eustatic rise, and this regime is similar to the present one (Serra 1976).Relict forms, the internal structure and the morphology of the sea floor (Fig. 10 and Fig. 6) offer new insights into the evolution of the fluvial-deltaic River Tordera system during the Holocene transgression.During the eustatic rise, the resulting bed forms migrated upwards to lower depths.The first Holocene delta lobe, A1, corresponds to sequence Sd1.The C1 and C2 sand bodies are characterized by an external steep seaward slope and an internal concave morphology towards the mainland (Figs 1 and 10) and can be correlated with Sd2 and Sd3, respectively.These features bear a striking resemblance to those described in the spit at the mouth of the River Tordera (bar of Malgrat) and in the relict bar system on the Maresme shelf (Serra et al. 1990, Serra and Sorribas 1993, Díaz and Maldonado 1989).The bars in C1 and C2 must have corresponded to a dynamic situation similar to the present one.However, transgressive conditions caused a retreat of the Tordera accreted riverbed from the Blanes Canyon head to its present position.The palaeochannels at the base of unit Sd1 to the NE and to the SW (Figs 6 and 7) indicate that the former courses of the River Tordera flowed in both directions.It is therefore possible to establish a relationship between the River Tordera and the Blanes Canyon during the last low stage of the sea level and the subsequent eustatic rise.At that time, the sedimentary contribution of the River Tordera constituted the main source of sediments to the Blanes Canyon and to the relict inner shelf sand bodies.
The B1 relict and the present littoral prisms of B2 can also be correlated with Sd3 and Sd3sup, respectively.
As regards sediment dynamics, the topography of the NE Blanes coast acts as a physical barrier against drift currents from the north.La Planassa High also constituted a barrier against drift currents from the NE during the last sea level low and subsequent rise (Figs 1and 10), which is evidenced by the lack of detrital sediment from the NE (Serra-Raventós and Got 1976).This argument is based on the high content of present and relict carbonatic facies (maerl) on the inner shelf between the Blanes Canyon and the coast (Serra 1976, ITGE 1989) and on their growth requirements (Bosence and Wilson 2003, Blake and Maggs 2003, Wilson et al. 2004).This finding contrasts sharply with the sand facies of the SW Tordera coast and the Maresme shelf, which is fed by River Tordera sediment (Serra 1976, Pedrosa-Pàmies et al. 2013).
The sedimentary and seismostratigraphic units obtained from the Tordera delta wells data and marine seismic profiles, respectively, are interpreted to be formed during the Holocene after correlation with the regional eustatic curves (Aloïsi 1986, Verdaguer 1983).The ages assigned to the sequence limits are represented in Figure 11.The erosive surface underlying sequence Sd1 is regionally traceable throughout the continental shelf because of its high erosive signature (PH, Fig. 8).This surface has been described by earlier works (Got 1973, Gámez et al. 2009, Ercilla et al. 2010), and has been interpreted in this work as the base of the transgressive accretionary system tract SD1 (Fig. 6).Chronologically, earlier works established the PH dates as 14.5 ky.In this work the base of the Sd1 unit (Fig. 8) has been assigned to 14 ky, and the top of Sd1 corresponding to 11.5 ky.Both ages are assigned after correlation with well P6-S5 and the eustatic curve (Figs 11 and 2).
Sd3ll is composed of pelitic subunits, Sd3ll1 and Sd3ll2.The age of Sd3ll1 ranges between 8900 and 7500 years BP, during which the sea level rose between -33 and -17 m (P6-S5 well, Fig. 2).This unit consists of grey-blue clays and silts with traces of marine organisms and organic matter, and becomes thinner and coarser towards the delta plain.
Subunit Sd3ll2 is located between -13 and -7.5 m (B-3-A well, Fig. 2) and ranges in age between 7600 and 6800 YBP.It was coeval with Sd3ll1 for less than 300 years.It consists of clay and silt with abundant organic matter but without marine bioclasts.The limit between Sd3p and Sd3ll at 7800 YBP is equivalent to the maximum flooding surface at the Rhone and Llobregat deltas defined by Boyer et al. (2006) and Gámez et al. (2009).Above this level, Sd3p and Sd-3sup correspond to the top Tordera delta units.In accordance with the eustatic curve, the sea level rose up to 9 m above the Sd3ll2 unit.The absence of marine bioclasts suggests that Sd3ll2, Sd3p and Sd3sup are fluvio-deltaic units.Sd3ll1 and Sd3ll2 correspond to brackish and fresh coastal lagoon systems, respectively.These two lagoons were coeval during the delta plain progradation but were separated by the aforementioned granite structural high.When the brackish lagoon was flooded owing to the continuous rise in sea level, the sea was obstructed first by the granite high and subsequently by the prograding units, Sd3p and Sd3sup, which resulted in the downlap cover of the two lagoons.

CONCLUSION
The correlation between the sedimentary and seismic data from the delta and the inner shelf of the River Tordera system and their chronology enable us to define three depositional sequences of the Holocene (Sd1, Sd2 and Sd3).These sequences are made up of several subunits that are characterized by their coarsening upwards, progradational and transgressive nature, and are environmentally conditioned by an unconformable pre-Holocene (Plio-Pleistocene) substrate or crystalline basement.
The current rich and relict carbonate productivity (maerl) and the lack of detritic sediments on the inner and middle shelf between Lloret, La Planassa High and the Blanes Canyon head provides strong evidence of the absence of sedimentary drift transport from the NE.The chronology of sequences Sd1, Sd2 and Sd3 is established after its correlation with the local sea level eustatic curve.The construction of the present delta began with unit SD1 at 14000 YBP.

Fig. 2 .
Fig. 2. -Situation Map of the 50 wells and piezometers on the Tordera Delta plain under study.The colours indicate the different types of well classified as a function of the number of stratigraphic sequences and units present.

Fig. 3 .
Fig. 3. -Standard lithological reference section and depths (comprising all units) of the Tordera Delta Holocene sequences, as recorded at the P6-S5 well.The map in Figure 2 shows the three well types as a function of the sequence content: pink contains all the sequences (as in P6-S5 well), orange does not contain Sd1-Sd2ll (as in ACA-4 well), and green only contains Sd3 and Sd2sup sequences (as in B-3-A well).

Fig. 7 .
Fig. 7. -A, image of slope gradients (modified from UB-Pridesa in Alfaro 2011) and situation of different seismic profiles: P1-P4 and PL.B, interpretation High-resolution seismic profile (see Fig. 5) showing three sequences Sd1, Sd2 and Sd3, PH and basement of the Tordera prodelta and inner shelf.A large paleochannel on the erosive base of Sd1 can be distinguished on P4.

Fig. 11
Fig. 11.-Standard reference section of the Tordera Delta Holocene sequences recorded at P6-S5 well.The location of the P6-S5 well can be seen in Figure 2.