The Zanclean mega ﬂ ood of the Mediterranean – Searching for independent evidence

About six million years ago, the Mediterranean Sea underwent a period of isolation from the ocean and widespread salt deposition known as the Messinian Salinity Crisis (MSC), allegedly leading to a kilometer-scale level drawdown by evaporation. One of the competing scenarios proposed for the termination of this environmental crisis 5.3 million years ago consists of a mega ﬂ ooding event re ﬁ lling the Mediterranean Sea through the Strait of Gibraltar: the Zanclean ﬂ ood. The main evidence supporting this hypothesis is a nearly 390 km long and several hundred meters deep erosion channel extending from the Gulf of Cádiz (Atlantic Ocean) to the Algerian Basin (Western Mediterranean), implying the excavation of ca. 1000 km 3 of Miocene sediment and bedrock. Based on the understanding obtained from Pleistocene onshore mega ﬂ ooding events and using ad-hoc hydrodynamic modeling, here we explore two predictions of the Zanclean outburst ﬂ ood hypothesis: 1) The formation of similar erosion features at sills communicating sub-basins within the Mediterranean Sea, speci ﬁ cally at the Sicily Sill; and 2) the accumulation of the eroded materials as mega ﬂ ood deposits in areas of low ﬂ ow energy. Recent data show a 6-km-wide amphitheater-shaped canyon preserved at the Malta Escarpment that may represent the erosional expression of the Zanclean ﬂ ood after ﬁ lling the western Mediterranean and spilling into the Eastern Basin. Next to that canyon, a ~1600 km 3 accumulation of chaotic, seismically transparent sediment has been found in the Ionian Sea, compatible in age and facies with mega ﬂ ood deposits. Another candidate mega ﬂ ood deposit has been identi ﬁ ed in the Alborán Sea in the form of elongated sedimentary bodies that parallel the ﬂ ooding channel and are seismically characterized by chaotic and discontinuous strati ﬁ ed re ﬂ ections, that we interpret as equivalent to gravel and boulder megabars described in terrestrial mega ﬂ ood settings. Numerical model predictions show that sand deposits found at the Miocene/Pliocene (M/P) boundary in ODP sites 974 and 975 (South Balearic and Tyrrhenian seas) are consistent with suspension transport from the Strait of Gibraltar during a ﬂ ooding event at a peak water discharge of ~10 8 m 3 s


A pre-scientific myth about the origin of the Mediterranean Sea
In his Historia Naturalis (~77 AD), Pliny the Elder reports on a legend popular among the inhabitants of southern Iberia, telling how the Mediterranean Sea was born.According to this myth, the Mediterranean Sea was deserted and cut-off from the Ocean, and it was Hercules who dug an inlet with his sword between Jebel-el-Mina (Africa) and the rock of Gibraltar (Europe).This allowed the ocean to flow into the Mediterranean Basin, where "it was before excluded", thus "changing the face of Nature".All through the history of western culture, a myriad of later classical writers including Galileo Galilei's Dialogues (1632) or Jacinto Verdaguer's poem Atlàntida (1876) were inspired by Pliny's account, elaborating on the scenario of a desiccated Mediterranean Sea flooded by Atlantic waters.None of these accounts linked explicitly the myth to field evidence.Even after Steno and Lyell set the principles of Geology and the Messinian stage was recognized as a pan-Mediterranean evaporitic phase in the late nineteenth century by Mayer-Eymar in 1867 (Selli, 1960), nobody linked it with Pliny's accounts, perhaps because their catastrophic nature was at odds with the principle of gradualism, a scientific pillar deeply rooted in the birth of Geology.
Gradualism was eventually challenged when Bretz (1925) and Pardee (1942) set the floor for a paradigm change in geology recognizing outburst floods of unprecedented magnitude as the main agent for landscape formation in the Scablands (Washington state, NW USA) during the Pleistocene.Bretz and Pardee described flood deposits, giant erosion coulees and giant 'ripples' that are today widely accepted to be the result of the catastrophic emptying of Lake Missoula in Montana, about seventeen thousand years ago (Benito, 2003).

The development of the Messinian Salinity Crisis theory
Only in the decade of the 1970's did the hypothesis of a flood ending a Mediterranean desiccation enter the scientific literature, when DSDP drilling confirmed large deposits of evaporites in the deepest regions of the Mediterranean Sea (Hsü et al., 1973b;Ryan et al., 1973).The question as to how did this Messinian Salinity Crisis (MSC) end emerged from the barely 15 cm M/P transition in ODP 974B (Tyrrhenian Sea; Fig. 1), which meant an "abrupt palaeoceanographic change created by the re-establishment of an open-marine environment in the earliest Pliocene" (Iaccarino and Bossio, 1999).The term flood or reflooding of the Mediterranean became a common place in the scientific literature studying the MSC (Hsü et al., 1973a, b;Comas et al., 1996), although it referred vaguely to a geologically-rapid transition from evaporitic to open marine conditions.Sedimentological analysis of the drilling cores (grain-size and microfauna) cannot resolve gradual changes shorter than a few millimeters in the core (i.e., a few thousand years in duration).For example, based on DSDP sites 121-139, K. J. Hsü et al. (1973a) speculated about a refill of the Mediterranean lasting "probably less than 1000 years".Onshore evidence for a rapid Mediterranean Sea level rise at the end of the MSC came from the abundance of Giberttype deltas in the Gulf of Lyons and the Tyrrhenian Sea (Clauzon, 1978(Clauzon, , 1982;;Breda et al., 2009) and from the abrupt transit from paleosol to foraminifera-rich, deep marine laminites in the M/P contact in Pissouri (Orszag-Sperber et al., 2000).These observations suggested that the sea level rise was rapid in sedimentological terms and that the Strait of Gibraltar was the scenario for the reconnection between the Atlantic Ocean and the Mediterranean Sea, reestablishing the normal marine conditions ca. 5.3 million years ago.
The main evidence for a restriction from the Ocean and a Mediterranean-wide environmental change is the widespread presence of gypsum deposits along the Mediterranean coast, and a 1 to 3 kmthick layer of halite in the abyssal areas of the Mediterranean, generally interpreted as the result of a severe restriction of the connection to the Atlantic Ocean.This halite layer is outcropping today in Calabria and it is widely present in the subsurface of the Caltanissetta Basin in Sicily, where it is commercially exploited in salt mines, and it is conspicuous in the deep Mediterranean basins from marine seismic reflection data (Lofi, 2010(Lofi, , 2018)).Its mass had been estimated at > 3 10 18 kg (Blanc, 2000;Ryan, 2009), around 10% of the salt contained in the global ocean, but a recent study based on a dense compilation of seismic prospection surveys, estimates its volume at only 821 to 927 10 6 km 3 of (Haq et al., 2020), ca.4% of today's salt dissolved in the ocean.
However, salt precipitation does not necessarily imply a desiccation.Natural salt pans do not require a level drawdown by desiccation in order to accumulate salt; a persistent but restricted connection that allows no outflow of brine but only inflow of normal saline water can lead to salt precipitation with no drawdown (e.g., today's Garabogazköl in Turkmenistan).The main evidence actually supporting a kilometric drawdown for the Mediterranean during the late Neogene is the Mediterranean-wide Messinian Erosion Surface (MES) (Lofi, 2010(Lofi, , 2018) ) recorded in the seismic stratigraphy of the basin margins.
The Zanclean flood hypothesis is thus entangled with a wider discussion about the MSC: whether there was a kilometric-scale drawdown after the precipitation of the marginal primary gypsum.The Zanclean (5.33 to 3.60 million year ago) is the first stage of the Pliocene and its beginning coincides with the termination of the Messinian (last stage of the Miocene).A megaflood-like termination of the MSC requires a significant drawdown of the Mediterranean level.The occurrence of a significant drawdown of the Mediterranean during the Messinian (be it a full or a partial desiccation) is currently under strong debate and far from a consensus.Let us start summarizing the main reasons that question the kilometric drawdown scenario.A geochemical reexamination of the aforementioned drillings has questioned the shallowwater formation of the deep evaporite cores (Hardie and Lowenstein, 2004).Some studies suggest that the evaporitic deposition occurred under no substantial sea level drawdown but in a deep-basin under deep-water in a deep-basin condition, i.e., in a Mediterranean Sea filled to a level close to today's (e.g., Lugli et al., 2015;Roveri et al., 2001).This model is also supported by the homogeneity of isotopic signatures between distant synchronous gypsum deposits, suggesting a connection between the Mediterranean sub-basins and hence a high water level (García-Veigas et al., 2018).
Questioning a significant lowering of the level of the Mediterranean requires a mechanism to explain the MES alternative to subaerial exposure.The most developed such alternative mechanism so far is the hyperpycnal submarine cascading erosion (Roveri et al., 2014b), according to which the erosion of the margins resulted from sea water currents from the platforms towards the abyssal areas as a result of high density related to evaporation.As for the abrupt M/P transition, the alternative explanation to the Zanclean flood is a density-driven outpumping (Marzocchi et al., 2016) of the Mediterranean Sea.In this model, the Mediterranean level raises due to fresh water inputs from the Paratethys during the Lago-Mare stage, reconnecting and pumping it out to the denser Atlantic Ocean.Salt deposition was followed by a pan-Mediterranean phase of sedimentation in a fresh to brackish water environment, represented by the so-called "Lago-Mare" sedimentary facies, which includes microfossils and ostracods originating from the Paratethys realm (present Carpathian and Black Sea areas) (Krijgsman et al., 2010).This deposition preceded the reestablishment of normal open marine conditions (Zanclean stage), and its presence at nearly the present sea level in outcrops in the margins around the Mediterranean Sea suggests that those fresh/brackish waters were covering a full and connected Mediterranean.However, the nature and environment of this last stage of the MSC are to date far from well understood since it also includes marine fish (Carnevale et al., 2006(Carnevale et al., , 2017)).
Studies that argue for a kilometer-scale drawdown (not necessarily implying a full desiccation) place it temporally either following (Escutia and Maldonado, 1992;Ryan, 2008;Estrada et al., 2011;Urgeles et al., 2011) or preceding (Bache et al., 2009(Bache et al., , 2015) ) the precipitation of the 1 km of halite in the abyssal regions.Within the drawdown hypothesis, the erosional surfaces in the margins of the sea are attributed to subaerial processes such as river incision or wave erosion.Further evidence supporting is based on the relative abundance of volcanic samples dating to the MSC around the Mediterranean, which according to modeling results can be explained by the decompression of the crust due to the removal of the water column (Sternai et al., 2017).
As for the deep geodynamic causes for these changes in oceanic connectivity at the Strait of Gibraltar, computer simulations suggest that desiccation of the basin, as well as its later reflooding at the earliest Pliocene, may have been caused both by the vertical motions related to deep processes in the lithospheric mantle, namely the sink and detachment of a lithospheric slab underneath the Gibraltar Arc (Garcia-Castellanos and Villaseñor, 2011;Jiménez-Munt et al., 2019).The combination of these processes with erosion by the inflowing water at the gateway may explain the long initial stage of salt precipitation (gypsum and halite) and an extremely rapid basin refill.
The aim of the present article is to summarize recent and new stratigraphic evidence and modeling results that test the megaflood hypothesis (Garcia-Castellanos et al., 2009).The geophysical data we examine target the erosional and depositional impact of a large water flow entering the Mediterranean through the Alborán Sea and then entering the Eastern Mediterranean over the Malta Escarpment.Because such a megaflooding event would require a significant sea level drawdown in the Mediterranean, we will discuss the implications of the findings for the evolution of the MSC.

First models and quantifications of the Zanclean flood
The quantification of a geologically-rapid refill after the MSC was addressed specifically only in this century, first based on the geomorphology of the Strait of Gibraltar area and then based on geophysical and modeling techniques (Blanc, 2002(Blanc, , 2006;;Garcia-Castellanos et al., 2009).According to Blanc (2002), an unprecedented stream shaped the Strait of Gibraltar after a capture of the Atlantic Ocean took place by retrogressive erosion of a river draining into the desiccated Alborán Sea (a process also known as river piracy).Based on assumptions on the size of the gateway during the Late Messinian and on collapse structures near the straits, that study concluded that the refill was cataclysmic, following an "instant breakthrough of the Atlantic waters".
Meanwhile, cores from drilling funded by the Africa-Europe tunnel project (Esteras et al., 2000) showed a > 200 m-deep erosional trough carved in the hard flysch units (Late Cretaceous to Neogene in age) outcropping around the strait.A similar eastward-oriented Pliocene incision is observed further to the east, in the Mediterranean side of the strait (Alborán Sea; Fig. 2).Both features had previously been interpreted as subaerial (fluvial) erosion by a forerunner stream during the Messinian sea level drawdown (Blanc, 2002;Esteras et al., 2000).However, the U-shape of the Messinian erosion surface (MES) in the Alboran Sea, the lack of a significant catchment, and the fact that this erosional channel crosses the natural drainage divide between the Atlantic and the Mediterranean (80 km in the Gulf of Cádiz and > 200 km in the Alborán Basin) argue against subaerial erosion during the MSC drawdown (Garcia-Castellanos et al., 2009).Numerical modeling of this erosion channel assuming that it was excavated by the Atlantic waters overtopping the Strait of Gibraltar, led to estimations of the flood discharge of up to 100 Sv (1 Sverdrup = 10 6 m 3 s −1 ), with a duration of 2 years or less (Garcia-Castellanos et al., 2009) (see Section 5).The initial stages of the flood, according to that model, may have lasted up to a few thousand years before the inflow of water becomes significant relative to the Mediterranean water balance.By the time the flood overcomes lake evaporation and significantly refills the Mediterranean basin, erosion rates are predicted higher than a few mm yr −1 , preventing any tectonic or eustatic change from closing the seaway again and making the flooding process irreversible.This argues against the possibility of multiple flooding events (Rouchy and Martin, 1992) or against calm periods intercalated during the flood.It conflicts with the interpretation of a ravinement surface in the Gulf of Lions as the result of wave erosion during the flood (Bache et al., 2009;Estrada et al., 2011) and suggests that this might have been formed during climatically-controlled changes of lake level during the drawdown phase, similar to the wave-cut surfaces in the Malta Escarpment that will be shown below (Micallef et al., 2019).The accelerating flooding rates due to inlet erosion (Garcia-Castellanos et al., 2009) precludes the stagnation of the flooding process once the water discharge has become substantial in comparison with the water budget of the Mediterranean and has significantly risen its level.
In summary, while observations remain unclear on whether the crisis involved an important desiccation of the Mediterranean Sea (necessary condition for a catastrophic flood), models do provide some ideas on how to independently test this scenario.In this article, we summarize recent results from seismic stratigraphy and hydrological modeling assessing two predictions implied by the Zanclean flood thesis: 1) that at a few thousand cubic kilometers of rock eroded from the sills should have deposited elsewhere in the Mediterranean Sea; and 2) that flood erosion and sedimentation patterns similar to those described in Gibraltar should be found near the Sicily Sill, in the area where the flooded western Mediterranean should overspill into the eastern basin.

Flood erosion features
Apart from the ubiquitous MES discordance in the seismic sedimentary record of the Mediterranean Sea, two erosional features have so far been described that seem incompatible with subaerial exposure during the MSC drawdown or with processes other than an unprecedented discharge of water.One is at the seaway between the Atlantic Ocean and the Mediterranean Sea, and the other at the gateway between the western and eastern Mediterranean basins.

Alborán Sea erosive channel
A remarkable erosive channel has been described in the deepest areas of the Messinian Erosion Surface (MES) of the Alborán Sea, particularly in the vicinity of the Strait of Gibraltar.Fig. 2 shows a significantly updated version of the depth to the MES in this region, based on a more dense and extended seismic reflection database relative to Estrada et al. (2011).The erosion channel is in lateral continuation with the erosional trough incised into the flysch units in the Strait of Gibraltar (Blanc, 2002;Esteras et al., 2000) and extends along 390 km with a W-E trend from the Gulf of Cádiz (Atlantic Ocean) (Garcia-Castellanos et al., 2009), crossing the deepest parts of the Alborán Sea (Figs. 2 & 3) until the Algerian Basin.This erosion channel geometry varies from 2 km width and > 200 m depth in the Atlantic side (were it dissects the basement) to > 15 km width and > 500 m depth in the Alborán Sea (were it erodes Miocene sediment).While Blanc (2002) linked the Atlantic section of this feature to the reflooding of the Mediterranean Sea, its Alborán continuation was linked to subaerial erosion by rivers during the MSC drawdown (Blanc, 2002;Campillo et al., 1992).Later investigations based on a wider geophysical database and numerical modeling lead to the alternative interpretation of an outburst flood erosion due to the breaching of the Straight of Gibraltar (Garcia-Castellanos et al., 2009;Estrada et al., 2011), after which the subsequent flood-retrogressive erosion propagated westwards, 80 km into the Atlantic domain.This hypothesis is consistent with the channel's U-shaped section, which contrasts with the V-shaped valley associated to fluvial erosion (Fig. 3).It also explains why this erosional feature has continuity westwards from the Strait of Gibraltar into the Gulf of Cádiz, crossing the natural divide formed by the Betic-Rif orogen.
The morphological setting of the Alborán Basin is inherited from Late Neogene E-W extension coeval with African-Eurasian convergence and with volcanism (Duggen et al., 2003).As a result, the relief of the Alborán Basin was already abundant in volcanic seamounts and structural highs (Fig. 2).Adapting to this inherited landscape, the channel splits in two W-E branches only to merge again further to the E (Fig. 2).They are separated by a seamount (the Vizconde de Eza High, VdE in Fig. 2) that may have been a topographic obstacle during the flood.The channel incises vertically up to 600 m into the upper Miocene deposits, eroding MSC deposits that do not include the halite evaporites (Mobile Unit) found in the deeper basins of the Mediterranean.This is probably due to the shallower depth of this basin during that time (e.g., Booth-Rea et al., 2018).The channel has been subsequently deformed during the Plio-Quaternary by mud diapirism in the Western Alboran Basin and the compressional tectonism as a consequence of the ca.4-5 mm yr −1 of the ongoing convergence rate between Africa and Eurasia (Estrada et al., 2018).

Amphitheater-headed canyon in the Malta Escarpment
The Noto Canyon is a ca.27 km long, 15 km wide, and 1.5 km deep geomorphic feature in the northern Malta Escarpment (SE Sicily margin; Fig. 4) (Micallef et al., 2018).It includes a canyon with a Ushaped cross-section carved in Mesozoic limestone.The head of this canyon has an amphitheater-shaped head, 6 km in width and 700 m high.This is very unusual for a submarine canyon and conspicuously different from the other canyons along the escarpment (Micallef et al., 2019).The strike of the Noto Canyon changes by 90°at a distance of 8 km upslope from its mouth.The canyon walls have a slope of 70°.Upslope of the Noto Canyon head, the seafloor has a gentler gradient of 7°and hosts a dendritic network of small tributary canyons.At the base of the Noto Canyon head, a sediment layer 100-200 m thick has accumulated at the bottom of the canyon (Fig. 4b).
The Noto canyon is located next to a set of massive, seismically chaotic to transparent deposits identified through multichannel seismic reflection data in the western Ionian Basin and interpreted as a sedimentary body of up to 860 m in thickness and 1600 km 3 in volume (see next section).Based on the inference that this sedimentary unit is compatible with a megaflood deposit, the Noto Canyon has been interpreted as having formed via excavation by a large subaerial flow of water caused by the overflow of the Western Mediterranean and its spilling into the Eastern Mediterranean basin during Zanclean flood (Micallef et al., 2018).However, this erosional feature is not directly dated and no stratigraphic correlation exists between its excavation and the Ionian stratigraphy and drillings.
Further to the south along the same escarpment, a series of concave slope breaks is observed at a consistent depth of 2400-2500 m below present sea level and with a total length of about 100 km (Fig. 5).These are interpreted as palaeoshorelines and shore platforms excavated by wave erosion related to a lowstand of the Mediterranean Sea below its present sea level (Micallef et al., 2019).An isostatic restoration of those palaeoshorelines and shore platforms on the northern Malta Escarpment suggest an evaporative drawdown of 1800-2000 m below present sea level.While these do not directly support the occurrence of an abrupt outburst flood refilling the eastern Mediterranean, the shorelines do support a large hydraulic head available for the refilling waters pouring into the eastern basin with a flow energy comparable to, if not larger than, that at the Strait of Gibraltar.The bottom of the Noto Canyon (2300-2600 m) is close to the depth of those shorelines and to the base of the megaflood deposit, allowing us to speculate that the canyon may have been excavated after a long period when the level of the Eastern Mediterranean remained at that level.After isostatic restoration by removing the weight of the Plio-Quaternary sediment and the water column (Micallef et al., 2019), the original depth of the base of the canyon the during the MSC would range between 1750 and 2050 m.

Hydrodynamic models for the fate of megaflooding sediment
The erosional channel that runs from the Gulf of Cádiz to the Algerian Basin (Fig. 2) implies the excavation of at least 1000 km 3 of Miocene sediment from the Alboran Basin and bedrock (Oligocene flysch) from the Strait of Gibraltar.The submarine threshold at today's Camarinal Sill is thought to have been incised by at least 200 m during the flood (Garcia-Castellanos et al., 2009).This number is key in constraining the peak velocity and discharge of the Zanclean flood (Fig. 6; Garcia-Castellanos et al., 2009;Periáñez et al., 2019).To show this, we calculate the feedback between the flow of water (calculated using Manning's equation) and the enlargement of the inlet by erosion (assumed proportional to basal shear stress).using the same algorithm as Garcia-Castellanos et al. (2009).We show two simulations of the Mediterranean refill using two values of sill rock erodability (Fig. 6).Although erodability k e is a still a poorly understood property of rocks and its quantification is limited by orders of magnitude of uncertainty (Garcia-Castellanos and O'Connor, 2018), the erodability value used for our reference model is constrained by the total incision of the flood, estimated at > 200 m in the present Camarinal Sill (Esteras et al., 2000), The results show that this depth of erosion is reached only if k e > 2 10 −4 m yr −1 Pa -1.5 , which in turn implies a peak water discharge close to 100 Sv (1 Sverdrup = 10 6 m 3 s −1 ).Smaller values of erodability fail to reproduce the observed depth of the erosional trough across Camarinal and the Strait of Gibraltar.
Nevertheless, the fate of these eroded materials has remained unknown, despite the fact that their present location might be an independent source of evidence to validate the megaflood hypothesis.A first attempt to predict the distribution of these flood deposits has been recently performed (Periáñez et al., 2019) using a Lagrangian 2D hydrodynamic model to simulate the transport of material eroded from the Strait of Gibraltar.The transport of rock particles of variable size, either in suspension or as bed load is calculated using standard sediment transport equations.Water circulation during the flood has been estimated from a hydrodynamic model of the whole Mediterranean Sea (Abril and Periáñez, 2016) based on the assumption that the inlet was enlarged by the flow itself (Garcia-Castellanos et al., 2009).The hydrodynamic model solves the 2D depth-averaged fluid dynamic equations incorporating wetting/drying processes.The model finds the topographic lows in the exposed seafloor along which water flows (Periáñez and Abril, 2015).The numerical solution of the equations yield values for the water velocity over the model domain and modifies water depth preserving the volume of water.Sea-bed stresses over the domain due to the water current are calculated for the suspended sediment transport.The incision model used for calculating the enlargement of the Gibraltar inlet is based on an energy balance for the energy dissipated by the stream of water (when gravitational potential energy decreases) and the energy required for eroding the unit volume of substrate.This explains the erosion of the former Gibraltar isthmus and  the Zanclean Channel (both the branch running through the Alborán Sea and the recessive incision channel extending to the Camarinal Sill in the Gulf of Cádiz).The approach uses as initial conditions the water level drop inferred from the Messinian shoreline reported by Urgeles et al. (2011), ca.1200 m, and adopts the present-day bathymetry of the Mediterranean as a proxy for the Messinian hypsometry.Sediment transport is described in a Lagrangian framework both for suspended load and bed load.The transport model for suspended load includes the processes of advection (transport with the water current), diffusion (mixing due to water turbulence), and settling of particles.Settling velocity depends on the particle size.Particles are deposited once they reach the bottom if the bed stress is lower than a critical deposition stress.This critical stress depends on the particle size as well.A detailed description of the model and its parameters (particle sizes, settling velocities, critical stresses) may be seen in Periáñez et al. (2019).Bed load transport occurs when a critical velocity is exceeded that depends on the local water depth and sediment size following (Bilgili et al., 2003).Particles are assumed to travel at one-sixth of the depth averaged current (Periáñez et al., 2019).
We show the areas of preferential sediment deposition in the Mediterranean Sea depending on grain size (Fig. 7) for sediments eroded from the Strait of Gibraltar and precisely at the stage of peak discharge, when the sea level has risen to −170 m.Suspended load is deposited in areas sheltered from the jet of incoming water by the local topography, areas where water current abruptly decreases due to a sudden increase in water depth, and the central zone of eddies, where water current magnitudes also decrease, as well as in areas were the flow changes from more to less restricted.Thus, main deposits are predicted at both sides of the Strait of Gibraltar (i.e., both sides of the Zanclean channel), the Alborán Sea central area and its northern and southern shores during the peak discharge.Once a particle leaves the Alborán Sea, it follows two main routes, one along the African coast and the second south of the Balearic Islands and Sardinia.Particles, as they are transported with the water jet, fall according to their settling velocity and finally they are deposited on the seabed in areas of lower flow energy.The distance of these deposits to the Strait of Gibraltar increases as particle size decreases.As a consequence, very extensive deposits of mainly silt and fine sand are predicted between the Balearic Islands and Sardinia due to energy loss and deposition along a northern current jet (Fig. 7c-e).Only a small fraction of sediment, particularly the finer particles, reaches the eastern Mediterranean through the Sicily Sill.Long deposits may also be seen along the Algerian coast, tracing the route of the stronger southern jet (Fig. 7a).
In contrast, sediment particles transported as bed load are deposited following the main streamlines.Deposits of very coarse sediment eroded from the Strait of Gibraltar and transported as bed load are more localized in space than sediment transported in suspension (Fig. 7b; Periáñez et al., 2019).In interpreting these model results, one must bear in mind that they are valid for the adopted bathymetry, which we take as the present one due to the absence of reconstructions of the Messinian paleo-bathymetry.Therefore, the relevance for the Zanclean flood is in the qualitative relationships between flow and deposition, rather than the exact location predicted for the deposits, particularly in areas of significant tectonic deformation such as the Tyrrhenian and the Alborán seas.

Identification of possible flood deposits in seismic reflection and borehole data
While the volume of eroded bedrock and sediment in Gibraltar and Alborán was estimated at 500 km 3 (Garcia-Castellanos et al., 2009), more recent data in Figs. 2 and 8 suggest this volume could be double that amount, since the erosional trough is now known to extend towards the east until the Algerian basin (Fig. 8).1000 km 3 of sediment evenly distributed over the extent of the western Mediterranean evaporites (4 10 5 km 2 ) would imply a sediment layer of only 2 m, undetectable in the available seismic stratigraphy.However, our model predictions suggest a strong localization of the flood deposits in areas of low energy (Fig. 7), and therefore larger size deposits are expected to be detectable in the seismic record and apparent in the core sedimentary

Possible megabar deposits in the Alborán Sea
Seismic reflection profiles show the presence of an elongated sedimentary body (Fig. 8C to D) associated to the erosion channel shown in Figs. 2 and 3.These deposits have contrasting seismic character with the underlying Miocene and overlying Pliocene deposits.Seismically they are characterized by chaotic, hyperbolic and discontinuous stratified facies.They form an elongated body up to 163 m thick, 35 km long and 7 km wide, located next to the main erosion channel and oriented parallel to its main trend (Fig. 8A).The Upper Miocene and Plio-Quaternary seismic stratigraphy of the Alborán Sea (Estrada et al., 2011;Juan et al., 2016), based on correlation with commercial, DSDP and ODP sites, indicate that this body is sandwiched between the MSC erosive surface (MES) and the Pliocene contourites deposited once the oceanic circulation was restored.The presence of these deposits occurs in the eastern Alborán basin, behind (in the lee side of) a pre-MSC volcanic edifice.
Their elongated plan-view shape, seismic facies, stratigraphic position and their orientation parallel to the flooding channel, are compatible with megabar flood deposits described for Pleistocene outburst flood scenarios such as the Altay floods (Carling et al., 2009;Herget, 2005).We propose that these deposits are related to the formation to the Zanclean erosion channel and represent possible megaflood deposits that were deposited because of the presence of a volcanic obstacle.This edifice may have acted as a topographic obstacle to the impinging flood flow, allowing the accumulation of sediments near the peak of the Zanclean flood, in the protected area at the lee side of the volcano.Similar sedimentary bars of larger grain size (boulder) have also been reported downstream the Snake River, spillway of the Bonneville flood in the Late Pleistocene (O'Connor, 1993), in areas where the flooding waters decreased in flow energy, for instance as a result of entering a less restricted topography.
Chaotic deposits have been reported as filling the bottom of the erosion channel in a seismic profile a few kilometers to the west of the section shown in Fig. 7b (Periáñez et al., 2019), and interpreted as the deposition of flood-erosion products during the late stages of the flood.These are not visible in the Conrad 823 seismic line, perhaps due to the smaller thickness that is undetectable in the seismic record, which in turn was the result of a more constrained flow imposed by the seamounts in this area.However, it must be noted that ODP site 977 (Figs. 1 and 7) drilled the lowermost part of the southern and northern erosion channels, respectively, reaching partially cemented rounded gravels and sands that contained a mixed fauna of Miocene and Pliocene, interpreted as possibly related to the Zanclean reflooding (Comas et al., 1996).ODP site 978 also reported pebble-sized clasts in the equivalent stratigraphic position.Model results in Fig. 7b predict that gravel-to-pebble-sized clasts can be dragged for a few hundreds of km at peak discharge, thus suggesting that the pebbles recovered from the drilling cores could come from the numerous volcanic sea mounts located further to the west, in the center of the Alborán Basin.

Possible megaflood deposits in the Ionian Sea
Multichannel seismic reflection data in the western Ionian Basin exhibit an acoustically chaotic to transparent body of up to 860 m in thickness and 1430-1620 km 3 in volume, covering an area of 11,000 km 2 buried at depth between 3.4 and 4.8 km below sea level (Fig. 9; (Micallef et al., 2018)).Towards the west, this body terminates abruptly against the Malta Escarpment.Although this unit has not been drilled, its stratigraphic position is placed with confidence between the Messinian Upper Unit (UU) and the overlying 200-400 m thick Plio-Pleistocene sediment layer characterized by well-laminated, continuous reflectors.Therefore, its age is at or near the Miocene/Pliocene boundary (5.33 Ma).
The scale of this feature cannot be readily explained by collapse of the eastern margin of the Malta Escarpment.A large submarine mass movement has been discarded as a possible explanation for the origin of Unit 2 due to the lack of nearby scars of sufficient volume and the hard (limestone) lithology of the escarpment (Micallef et al., 2018).A formation by folded Messinian upper gypsum is also discarded based on the seismic velocities measured (2.3-2.6 km/ s), which are significantly lower than the 3.1 km/ s associated with gypsum deposits.Finally, the lack of local sources of sediment (the canyons along the Malta Escarpment combined add up to only 1100 km 3 in volume), suggests that this unit may be a megaflood deposit.The seismic texture and the isopachs of the body indicate that this sedimentary unit represents highenergy deposits sourced in the W, beyond the Malta Escarpment.Its location next to the Noto amphitheater-shaped canyon (SE Sicily margin; Fig. 9) suggests a NW provenance of the flow, likely through the Pelagian Platform.Thus, the interpretation by Micallef et al. (2018) is that the chaotic body records deposition from a megaflood flow as the water refill of the western Mediterranean overtopped the Sicily Sill.

Significance for the Messinian salinity crisis
The Zanclean flood hypothesis is closely entangled to the controversy about the desiccation of the Mediterranean Sea during the Messinian salinity crisis.The occurrence and extent of this desiccation by the end of the Neogene is far from a consensus among the scientific community (e.g., Roveri et al., 2014a).The main evidence questioning a kilometric-scale drawdown right before the Pliocene is the ubiquitous presence of Lago-Mare facies (fresh to brackish water fauna that also include marine fish; Carnevale et al., 2017) during the last stage of the MSC, at elevations close to normal sea level (Orszag-Sperber, 2006).The Lago-Mare is characterized by mollusks, ostracods and/or dinoflagellate cysts of Paratethyan-type (originating from the Black Sea region) indicating a brackish sedimentary environment before the return to normal open-marine environment across the Miocene-Pliocene boundary.The widespread presence of these sediments from southern Spain to the Levant Basin suggests that, right before the Pliocene, the Mediterranean was a single connected body of water with its level above the Sicily Sill (today at −430 m) and that its salinity was lower than the open ocean.To explain the abrupt end of the MSC without a Correlation with DSDP drilling site 374 indicates that Unit 1 corresponds to Pliocene-Quaternary pelagics, contourites and turbidites; Unit 3 correlates with MSC evaporites including a highly reflective and mostly deformed upper unit made of gypsum and marls (3a, green), and a typically reflector-less unit with a basin-fill geometry made of halite and also known as mobile unit (3b, yellow); a lower, highly reflective, discontinuous unit (3c, purple) has never been drilled yet but is attributed to lower gypsum or clastics.Modified from (Micallef et al., 2018).(For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)significant drawdown that can trigger a flood, Marzocchi et al. (2016) proposed a Mediterranean-wide hydrological pump model driven by the strong salinity contrast between the Mediterranean, the Atlantic, and the Paratethys.This can explain the co-existence of dense Atlantic water influx from Gibraltar and lighter brackish water influx from the Paratethys during the Lago-Mare stage prior to the end of the MSC, all under normal sea-level conditions.Supporting this view, the ages of the Neogene marine seaways onshore the Alborán Sea suggest that the Gibraltar Strait (Krijgsman et al., 2018) may have never been fully closed (Manzi et al., 2005), questioning the idea of a completely isolated Mediterranean Sea during the MSC.In order to keep a kilometric drawdown for a period longer than a few thousands of years, a full disconnection of the marine gateways is required (Garcia-Castellanos and Villaseñor, 2011), because otherwise the inflow of water would trigger a fast erosion and reestablish the level of the Mediterranean Sea.Finally, models excluding a significant sea level drawdown explain the erosion features of the Mediterranean margins as the result of hyperpycnal density currents related to storage of high-density water due to evaporation during the MSC (Roveri et al., 2014b).A possible explanation for the intriguing Lago-Mare facies is that it post-dates the flooding.Based on outcrops near Málaga (southern Spain), Do Couto et al. (2014) suggest that the Lago-Mare sediments were deposited after the Atlantic flooding, which would then have taken place at 5.46 Ma instead of the Messinian/Zanclean boundary (5.33 Ma).
To add to these conflicting observations, recent industrial seismicreflection data-sets in the deep Ionian Sea (Sirt Basin, Libya) and in the deep Levant Basin (offshore Syria) show the existence of large fluvial sedimentary systems formed after the halite precipitation (Bowman, 2012;Madof et al., 2019).In the Levant, these fluvial systems are longer than 100 km and interpreted as fluvio-deltaic lobes resting directly atop the Messinian halite and seem coeval with the erosion of deformed salt layers (Bowman, 2012;Kartveit et al., 2019;Madof et al., 2019).These are interpreted as the result of vast paleo-river basins draining areas of Africa and the Middle East region during the lowering of the Mediterranean Sea level.In contrast, seismic stratigraphy in the same area suggests that salt promontories formed by salt flow in the sea floor were truncated by dissolution under a deep, stratified water column, conflicting with the possibility of a seafloor exposure (Gvirtzman et al., 2017).In the absence of drill holes and age determinations, the presumed fluvial lobes are ambiguously correlated to the Lago-Mare phase, but could they alternatively predate it and predate the megaflood?The seismic stratigraphic record to date fails to answer this question so far and points to the need for scientific drilling.
The evidence shown in this chapter supports an abrupt, catastrophic refill of the Mediterranean in the form of an outburst flood caused by the overtopping of the Atlantic through a Neogene isthmus in Gibraltar (Garcia-Castellanos et al., 2009).The same mechanism has been described for tens of Pleistocene overtopping lakes that triggered outburst floods (Garcia-Castellanos and O'Connor, 2018).Such scenario requires a significantly low level in the Mediterranean Sea (hundreds of meters at least) prior to flooding, to be consistent with the erosion and flood deposits reported in the Alborán area and the Malta Escarpment.But it does not require a complete desiccation of neither of the Mediterranean sub basins.

Comparison to other megafloods
The Missoula floods in the Scablands region (Washington, USA) are the first studied and best-known outburst floods that have.Based on information as diverse as the elevation of ice-rafted boulders, megaripple marks, the size of mobilized blocks and 0D step-backwater and 3D hydrodynamic modeling, several tens of events of up to around 10 7 m 3 s −1 in discharge (O'Connor and Baker, 1992;Benito and OʼConnor, 2003;Alho et al., 2010;Larsen and Lamb, 2016) have been described in this region dating between 18.5 and 15 thousand years ago (late Pleistocene).Erosional and deposition features, including km-scale regressive erosion cliffs as the Dry Falls or the Potholes Coulee are attributed to a series of outburst events involving a water flow deeper than 100 m at speeds up to 15 m/ s.Such horseshoe-or amphitheatershaped cliffs, frequent smoking guns for megaflooding, are similar to the Noto Canyon in the Malta Escarpment (Fig. 4).Such features are found widespread across the state of Washington and attributed to the Missoula floods.The Pothole coulee is ca. 2 km wide and therefore narrower than the Noto Canyon (6 km), a priori consistent with the 10 times higher peak discharge derived for the Zanclean flood (Garcia-Castellanos et al., 2009;Abril and Periáñez, 2016).Similar amphitheater-headed valleys have been linked to megaflooding events in Idaho (Lamb et al., 2008).Because of the current impossibility to date the formation of the Noto Canyon, the only indication linking its origin to the Zanclean flood is its morphology, its location next to the possible megaflood deposits in Fig. 9, and the absence of other explanations for its formation (Micallef et al., 2018).
The Altay flood along the Katun River (Altay Republic, Russia) is another case scenario where discharge rates comparable to the Missoula Floods (~10 Sv) have been estimated.Here, outburst-flood sedimentary features along mountain valleys are well exposed due to tributaries of the main flooding channel that were backflooded during the events.An abundant depositional feature is megabars, equivalent to fluvial bars but two orders of magnitude larger and thought to form near the time of peak discharge of outburst floods.The concept in Fig. 10 is based on a gravel megabar exposure cut by a tributary river after the flooding.It shows three cycles of flood deposition separated by slack water deposits, interpreted as an indication for repeated flooding.Accretion results from intense fallout of suspended gravels and is then rearranged by unsteady traction transport during upper stage surges.As the flood discharge wanes, cross-beds may form due ripple and dune migration during the lower stage (Carling et al., 2009).In the analogue case of the Zanclean flood, one single flooding event is envisaged.If, as proposed, the sedimentary feature shown in Fig. 8 corresponds to an equivalent megabar formed during the transition from high to low energy flow in the Alborán Sea, then an upward-fining sequence of mostly suspended sediment would be expected within that unit.However, the potential megabar shown in Fig. 8 has not been drilled so far, and the existing seismic surveys do not have enough resolution to detect grain-size changes or features like those in Fig. 10.Therefore, while we are unable to further check the consistency between the unit identified in Alborán Sea and the megabars in the Altay, the morphological distribution of the sediments (elongated, parallel and next to the main channel) make it an important candidate site for future drilling to find supporting lithological evidence for the Zanclean flood deposits.

Future search for further evidence
A fundamental peculiarity of the Zanclean flood setting relative to Pleistocene outburst floods is that, whereas its water source (the global ocean) is virtually infinite in volume, the sink (the desiccated Mediterranean Sea) is not.During the flood, the western shoreline of the western Mediterranean should raise and shift towards the Strait of Gibraltar, possibly hindering the spatial concentration of deposits.Nevertheless, because of the large size of the Mediterranean basin, the ever-increasing flooding discharge becomes significant to the water budget only when the flood discharge becomes comparable to hydrological deficit (presently 75,000 m 3 s −1 ).Therefore, future geophysical exploration of the western margin of the Algerian basin should address the presence of buried Zanclean flood deposits postulated by the hydrodynamic models.
The numerical modeling results in Fig. 7 provide an explanation for the 15 cm M/P transition found in ODP site 974 in the Tyrrhenian Sea (Fig. 1 for location) (Iaccarino and Bossio, 1999;Pierre et al., 2006) and the 2-3 cm of > 0.3 mm of terrigenous clasts (medium sand) recovered just above the M/P boundary in site 975 at the South Balearic Basin (Iaccarino and Bossio, 1999;Zahn et al., 1999) consisting of "foraminifera and calcareous nannofossils that are reworked from the early marine Messinian, while the ostracod fauna reflects the original paleoenvironment".The inflow of Atlantic waters at unprecedented discharge rates close to 10 8 m 3 s −1 (100 Sv) are proven capable of transporting medium sand and foraminifera at such distances (Fig. 7).In using these models for searching flood erosional and depositional features in the Alborán sea, one must consider the uncertainty induced by the use of today's bathymetry as a proxy for the paleogeography of the region.During the Messinian, most of the Tyrrhenian Sea, for instance, was still not formed by the retreat of the Calabrian Arc.In the Alboran Sea, the paleo bathymetry at the time of flooding may have been significantly different from today, since the lithosphere was younger, thinner, hotter and therefore more buoyant (Booth-Rea et al., 2018), and the compressional tectonics of the ongoing convergence between Africa and Eurasia has raised most of the Alborán Ridge during the Plioquaternary.Since the Messinian, Africa has approached Eurasia by 20 to 100 km (in the west and the east of the Mediterranean Sea, respectively), and the bathymetry has been modified by a layer of Plio-Quaternary deposits of, typically, a few hundreds of meters.The lack of good reconstructions of the bathymetric and geographic changes of the Mediterranean Sea at the time of the flooding limits a more detailed prediction of the locations ideal for future megaflood deposits exploration.
The present review shows recent independent assessments of the Zanclean megaflood hypothesis, and highlights the need for further future research in order to test it.Seismic studies may validate it if clinoform deposits and found in the westernmost Algerian basin and dated to the end of the MSC, for example.But a more direct test may arrive from drilling the sedimentary deposits described in the Alborán basin and in the Ionian Sea.In this regard, the outcome of three drilling proposals (DEMISE, DREAM, and IMMAGE) that are currently undergoing a second round of review by IODP, will be crucial.If approved and drilled, they have the potential to decipher many of the aforementioned open scientific questions regarding the end of the Messinian salinity crisis.

Conclusions
We have summarized a set of observations that provide additional evidence to the previously identified erosional channel at the Strait of Gibraltar, supporting a megaflood of the Mediterranean basin near the limit between the Miocene and the Pliocene.Based on them, we can draw the following conclusions: 1.If the > 200 m erosion at the Strait of Gibraltar was caused by the refilling of the Mediterranean Sea, numerical modeling shows that this erosion implies a flood of an unprecedented discharge, possibly above 100 Sv (10 8 m 3 s −1 ). 2. The amphitheater-headed Noto Canyon, off SE Sicily, is compatible with an overspill the Zanclean flood from the western to the eastern Mediterranean across the Sicily sill.The morphology of the canyon and the Malta Escarpment slope suggest a Mediterranean Sea level at the time of the flood −1750 to −1900 m lower than today.3. A 1600 km 3 sedimentary body to the east of the Noto Canyon, laying atop the Messinian salt layer and below the Plio-Quaternary muds is compatible with a megaflood deposit delivered to the eastern basin and bringuing to an end the drawdown phase.4. A 163 m thick, elongated sediment body of 35 × 7 km, oriented parallel to the main flood channel in the Alborán Sea and located in the lee side of a seamount, is compatible with megabars formed along the largest documented continental outburst floods.

Declaration of Competing Interest
None.

Fig. 1 .
Fig. 1. a) Topographic map of the Ionian Sea and western Mediterranean showing the location of data shown in this article.Red arrows indicate the average motion of Africa relative to Eurasia since 5.3 million years ago.Numbers in red indicate IODP/DSDP drillings mentioned in the text.b) Cartoon (not to scale) showing the timing and the flood erosional and depositional features expected in the eastern and western basins of the Mediterranean as a result of a large water input through the Gibraltar gateway, showing 5 stages: 0) Initial level before significant Atlantic inflow occurred; 1) time of maximum erosion rate at the Strait of Gibraltar; 2) the western basin level reaches the Sicily Sill.3) eastern basin level reaches the Sicily Sill; and 4) the Mediterranean is filled to the normal oceanic level.The mobile unit (Messinian halite), Lago-Mare, and flood deposits discussed in this study are depicted in yellow, green, and red, respectively.(For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)

Fig. 2 .
Fig. 2. Map showing the depth to the base of Pliocene, including the Messinian Erosion surface (MES) and the erosion channel (red line) resulting from the Zanclean megaflood and the main flood paths (orange line).The MES has a polygenetic origin and represents the base of Pliocene (purple line in the seismic profiles.Contours in milliseconds, colour scale bar in seconds.Updated from Estrada et al. (2011).(For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)

Fig. 3 .
Fig. 3. Composite seismic profile showing the Messinian erosion surface (MES; purple line) in the eastern side of the Strait of Gibraltar.Location in Fig. 2.This unconformity is interpreted as the erosion channel excavated into Miocene sediments at the latest MSC or earliest Pliocene.Depth in TWTT (double sound travel time in seconds).Detailed discussion and other seismic profiles in Estrada et al. (2011).(For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)

Fig. 4 .
Fig. 4. Possible gateway for the Zanclean flood into the eastern Mediterranean.(a) Map of south-eastern Sicily showing a seafloor slope gradient map of the Noto Canyon.Isobaths denoted by dotted lines.Location in Fig. 1.(b) Seismic reflection profile CIR-07 showing the Noto Canyon head and the Messinian erosional surface upslope.Location in (a).The sediment resting on the acoustic basement in the ENE (bottom of the canyon) is up to 200 m thick.Modified from Micallef et al. (2018).

Fig. 5 .
Fig. 5. (a) Multichannel seismic reflection profile across the Malta Escarpment intersecting the break of slope (dotted black line) and associated terrace, with the concave break of slope at the base.(Facies B = sub-parallel, convex upward reflectors; facies C = transparent facies).(b) Histogram of the depth of breaks of slope in the northern Malta Escarpment.The concave breaks are concentrated at depths of 2400-2500 mbsl.After isostatic restoration, these breaks have been interpreted as the result of wave erosion during the MSC lowstand of the Mediterranean level at −1800 to −2000 m below present sea level.Modified from Micallef et al. (2019).

Fig. 6 .
Fig.6.Modeling results of the evolution of the Zanclean flood based on the same methodology as inGarcia-Castellanos et al. (2009).These two simulations assume different rock erodability k e for the Gibraltar Sill.The refill of the Mediterranean starts with a 10 m deep inflow river of oceanic water.The reference model (bold lines) erodes this sill down to −200 m below the global sea level, whereas imposing half the erodability (thin lines) leads to a sill eroded by only 155 m, less than observed in geophysical data.a) Water velocity and inflow channel width at the Strait of Gibraltar.b) Water discharge and incision rate.c) Inlet depth (black lines) and level of the Atlantic, the Western Mediterranean, and Eastern Mediterranean (red lines).The Sicily Sill is flooded after 1.9 yr in the reference model (3.6 yr in the less erodible model).Most of the water flow is concentrated within less than two years in the reference model, implying sea level rise rates of up to 10 m/day in the Mediterranean.Circled numbers and vertical dashed lines show the five stages defined in Fig.1(for the reference model only).(For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)

Fig. 7 .
Fig. 7. Sediment transport model results corresponding to the peak discharge of the flood.Top left: logarithm of the water velocity (m/s) in the western Mediterranean at the peak stage of flooding.Bottom left: predicted distribution of bed-load deposits fed to the flow at the Strait of Gibraltar.Right: suspension particles sedimented per unit area (normalized to the maximum value).A total of 20 k particles are released for each grainsize of sand.Areas in white mean no grain is deposited in that location.Note the correlation between areas of deposition and areas of slow velocity.The results point to a strong concentration of most of the deposits in < 1% of the area of the Western Mediterranean.Modified from Periáñez et al. (2019).
These overtopping waters would plunge into the partially desiccated eastern Mediterranean basin, possibly at a level of −1750 to −1900 m as suggested by the wave erosion marks preserved in the escarpment, the base of the Noto Canyon, and the base of the flood deposits (Figs. 4, 5 & 9).

Fig. 8 .
Fig. 8. Seismic profiles showing possible megaflood deposits in the lee side of a volcanic edifice in the eastern Alborán Basin.a) Isobath map of base of the Pliocene (MES) and the main flood paths (orange lines).Location in Fig. 1.The shaded area maps the proposed flood megabar deposit.Contours every 150 milliseconds of double seismic travel time (TWTT), colour scale bar in seconds.b) Seismic cross-section Conrad 826 showing the two erosion channels separated by the volcanic building.Inset shows the location in a bathymetric map.c-d-e) Seismic profiles showing possible megaflood deposits in the lee side of a volcanic edifice in the eastern Alborán Basin.Legend: blue line represents base of Quaternary; The purple line marks the base of Pliocene (top Messinian); the green line is the base of hypothetical flooding deposits and the yellow line corresponds to Messinian salt deposition in the deep Algero-Balear Basin.Vertical scale in seconds (TWTT) and horizontal scale in meters.(For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)

Fig. 10 .
Fig. 10.Sedimentological studies in Pleistocene continental megafloods provide a model for future assessment of the megaflood origin of the possible Zanclean deposits in Fig.9.One frequent feature observed along the spillway of outburst floods consists of lateral megabars off-and parallel to the main flooding channel.This sedimentary facies model summarizes the depositional surfaces across a megabar along the Katun River in the Altay Mnts, were an ice-dam lake burst during the Pleistocene reached a discharge of ca.one million m 3 s −1 , with layers tilting towards the mountain side of the valley.Flood flow is directed towards the reader.The section has been exposed by incision along a tributary valley that had been backflooded during the events.Three cycles of flood deposition are visible separated by slack water silt deposits (black).For each cycle, fine gravel is deposited onto the flooded, rapidly-accreting bar surface, largely from intense fallout from suspension and is then rearranged by unsteady (surges) traction transport during upper stage flow.As flood flow wanes, cross-beds form due to bedload sheet and minor ripple and dune migration across predominantly at lower stage.The stratigraphy of this and other megaflood deposits provide a framework to interpret the outcome of future scientific drilling of the much older and buried Zanclean flood deposits.Modified fromCarling et al. (2009).