The Mediterranean Sea became disconnected from the world's oceans and mostly desiccated by evaporation about 5.6 million years ago during the Messinian salinity crisis. The Atlantic waters found a way through the present Gibraltar Strait and rapidly refilled the Mediterranean 5.33 million years ago in an event known as the Zanclean flood. The nature, abruptness and evolution of this flood remain poorly constrained. Borehole and seismic data show incisions over 250 m deep on both sides of the Gibraltar Strait that have previously been attributed to fluvial erosion during the desiccation. Here we show the continuity of this 200-km-long channel across the strait and explain its morphology as the result of erosion by the flooding waters, adopting an incision model validated in mountain rivers. This model in turn allows us to estimate the duration of the flood. Although the available data are limited, our findings suggest that the feedback between water flow and incision in the early stages of flooding imply discharges of about 10(8) m(3) s(-1) (three orders of magnitude larger than the present Amazon River) and incision rates above 0.4 m per day. Although the flood started at low water discharges that may have lasted for up to several thousand years, our results suggest that 90 per cent of the water was transferred in a short period ranging from a few months to two years. This extremely abrupt flood may have involved peak rates of sea level rise in the Mediterranean of more than ten metres per day.

-Quantified balanced and restored crustal cross-sections across the NW Zagros Mountains are presented in this work integrating geological and geophysical local and global datasets. The balanced crustal cross-section reproduces the surficial folding and thrusting of the thick cover succession, including the near top of the Sarvak Formation (∼ 90 Ma) that forms the top of the restored crustal cross-section. The base of the Arabian crust in the balanced cross-section is constrained by recently published seismic receiver function results showing a deepening of the Moho from 42 ± 2 km in the undeformed foreland basin to 56 ± 2 km beneath the High Zagros. The internal parts of the deformed crustal cross-section are constrained by new seismic tomographic sections imaging a ∼ 50 • NE-dipping sharp contact between the Arabian and Iranian crusts. These surfaces bound an area of 10 800 km 2 that should be kept constant during the Zagros orogeny. The Arabian crustal cross-section is restored using six different tectonosedimentary domains according to their sedimentary facies and palaeobathymetries, and assuming Airy isostasy and area conservation. While the two southwestern domains were directly determined from well-constrained surface data, the reconstruction of the distal domains to the NE was made using the recent margin model of Wrobel-Daveau et al. (2010) and fitting the total area calculated in the balanced cross-section. The Arabian continental-oceanic boundary, at the time corresponding to the near top of the Sarvak Formation, is located 169 km to the NE of the trace of the Main Recent Fault. Shortening is estimated at ∼ 180 km for the cover rocks and ∼ 149 km for the Arabian basement, including all compressional events from Late Cretaceous to Recent time, with an average shortening rate of ∼ 2 mm yr −1 for the last 90 Ma.

Three‐dimensional modeling that integrates fluvial sediment transport, crustal‐scale tectonic deformation, and lithospheric flexural subsidence is carried out to simulate the landscape and drainage evolution of the Ebro sedimentary basin (NE Iberia). The Ebro Basin underwent a long period of closed intramountain drainage as a result of tectonic topography generation at the Pyrenees, the Iberian Range, and the Catalan Coastal Range. In the late Oligocene, the Catalan Coastal Range underwent extension leading to the formation of the Valencia Trough (NW Mediterranean), but the Ebro Basin remained closed for nearly 15 Myr more before the Ebro River cut through the remnants of the topographic barrier. This drainage opening caused widespread basin incision that shaped spectacular outcrops of the syntectonic and posttectonic infill. Here we investigate the processes controlling these major drainage changes. The modeling results, constrained by a large data set on the tectonic and transport evolution of the area, predict a closed phase characterized by a large lake in the central eastern Ebro Basin. Dry climatic conditions probably lowered the lake level and contributed, together with rift flank uplift, to prolong this endorheic basin stage. The age and amount of reworked sediment after the opening are consistent with an onset of basin incision between 13 and 8.5 Ma as a result of lake capture by escarpment erosion and lake level rise associated with sediment accumulation and wetter climatic conditions. Sea level changes in the Mediterranean had no major impact in the large‐scale drainage evolution of the Ebro Basin.

The Messinian salinity crisis (5.96 to 5.33 million years ago) was caused by reduced water inflow from the Atlantic Ocean to the Mediterranean Sea resulting in widespread salt precipitation and a decrease in Mediterranean sea level of about 1.5 kilometres due to evaporation. The reduced connectivity between the Atlantic and the Mediterranean at the time of the salinity crisis is thought to have resulted from tectonic uplift of the Gibraltar arc seaway and global sea-level changes, both of which control the inflow of water required to compensate for the hydrological deficit of the Mediterranean. However, the different timescales on which tectonic uplift and changes in sea level occur are difficult to reconcile with the long duration of the shallow connection between the Mediterranean and the Atlantic needed to explain the large amount of salt precipitated. Here we use numerical modelling to show that seaway erosion caused by the Atlantic inflow could sustain such a shallow connection between the Atlantic and the Mediterranean by counteracting tectonic uplift. The erosion and uplift rates required are consistent with previous mountain erosion studies, with the present altitude of marine sediments in the Gibraltar arc and with geodynamic models suggesting a lithospheric slab tear underneath the region. The moderate Mediterranean sea-level drawdown during the early stages of the Messinian salinity crisis can be explained by an uplift of a few millimetres per year counteracted by similar rates of erosion due to Atlantic inflow. Our findings suggest that the competition between uplift and erosion can result in harmonic coupling between erosion and the Mediterranean sea level, providing an alternative mechanism for the cyclicity observed in early salt precipitation deposits and calling into question previous ideas regarding the timing of the events that occurred during the Messinian salinity crisis.

[1] We present an interactive 3-D computer program (LitMod3D) developed to perform combined geophysical-petrological modeling of the lithosphere and sublithospheric upper mantle. In contrast to other available modeling software, LitMod3D is built within an internally consistent thermodynamicgeophysical framework, where all relevant properties are functions of temperature, pressure, and composition. By simultaneously solving the heat transfer, thermodynamic, rheological, geopotential, and isostasy (local and flexural) equations, the program outputs temperature, pressure, surface heat flow, density (bulk and single phase), seismic wave velocities, geoid and gravity anomalies, elevation, and lithospheric strength for any given model. These outputs can be used to obtain thermal and compositional models of the lithosphere and sublithospheric upper mantle that simultaneously fit all available geophysical and petrological observables. We illustrate some of the advantages and limitations of LitMod3D using synthetic models and comparing our predictions with those from other modeling methods. In particular, we show that (1) temperature at midlithosphere depths may be overestimated by as much as 200 K when compositional heterogeneities in the mantle and T-P effects are not considered in lithospheric models and (2) the neglect of mantle phase transformations on gravity-based models in thin-crust settings can result in a significant overestimation and underestimation of the derived crustal thickness and its internal density distribution, respectively.