Impact of restriction of the Atlantic-Mediterranean gateway on the Mediterranean Outflow Water and eastern Atlantic circulation during the Messinian

[1] Messinian foraminiferal stable oxygen and carbon isotopes of the Montemayor-1 core (Guadalquivir Basin, SW Spain) have been investigated. This record is exceptional to study the Mediterranean Outflow Water (MOW) impact on the Atlantic meridional overturning circulation (AMOC) and global climate during the Messinian because the core is near the Guadalhorce Corridor, the last Betic gateway to be closed during the early Messinian. Our results allow dating accurately its closure at 6.18 Ma. Constant benthic dO values, high difference between benthic and planktonic dO, and low sedimentation rates before 6.18 Ma indicate the presence of a two-layer water column, with bottom winnowing due to an enhanced Mediterranean outflow current. The enhanced contribution of dense MOW to the North Atlantic Ocean likely fostered the formation of North Atlantic Deep Water (NADW). After 6.18 Ma, benthic dO values parallel that of the global glacioeustatic curve, the difference between benthic and planktonic dO is low, and sedimentation rates considerably increased. This indicates a good vertical mixing of the water column, interruption of the MOW, and a dominant glacioeustatic control on the isotopic signatures. According to the role of MOW in the modern Atlantic thermohaline circulation, the reduction of the MOW after the closure of the Guadalhorce Corridor might have resulted in a decreased NADW formation rate between 6.0 and 5.5 Ma weakening the AMOC and promoting northern hemisphere cooling. After the Gibraltar Strait opening, the restoration of the MOW and related salt export from the Mediterranean could have promoted an enhanced NADW formation.


Introduction
[2] At present, the Mediterranean connects with the Atlantic by the Strait of Gibraltar.The water mass exchange throughout the Strait of Gibraltar is characterized by an antiestuarine circulation pattern (Figure 1a) [Wüst, 1961].This anti-estuarine circulation pattern was definitively established after the opening of the Strait of Gibraltar [Nelson, 1990], although it has been suggested that the movements of water masses reversed to estuarine-type circulation during the early Pleistocene [Huang and Stanley, 1972].Lower salinity surface waters from the North Atlantic flow superficially eastward into the Mediterranean (Figure 1a).On the other hand, strong evaporation and production of dense, saline intermediate and deep waters in the eastern Mediterranean forces a high-velocity density driven bottom current westward as MOW through the Strait of Gibraltar.This mass of water is mainly fed by the Levantine Intermediate Water (LIW) [Bryden and Stommel, 1984], formed by convection in the Eastern Mediterranean [Marshall and Schott, 1999;Hernández-Molina et al., 2011], and in lesser extent by the Western Mediterranean Deep Water (WMDW), which is formed in the Gulf of Lion during cold and windy winters [MEDOC Group, 1970;Bryden and Stommel, 1984;Lacombe et al., 1985].Today, the LIW contributes to 2/3 of the MOW while WMDW only represents 1/3 of the MOW [Millot, 1999].Therefore, the LIW might be more important in controlling the MOW than WMDW.The MOW is characterized by higher d 13 C and d 18 O than the Atlantic waters [Vergnaud-Grazzini, 1983;Sierro et al., 2005].
[3] The outflow of the MOW into the Eastern Atlantic has a significant effect on the Atlantic oceanic circulation as well as on the global climate.The dense MOW mixes with the North Atlantic intermediate waters forming a high salinity tongue (Figure 1b) that contributes to the momentum of the Atlantic meridional overturning circulation (AMOC) [Reid, 1979;Rahmstorf, 1998;Bigg and Wadley, 2001;Bigg et al., 2003].In turn, the AMOC is the driving force for the Atlantic Ocean circulation, and even the global thermohaline circulation [Bethoux et al., 1999].Global thermohaline circulation affects the global radiation budget and global carbon cycling and can thus produce major climate changes [Brown et al., 1989;Bigg et al., 2003;Murphy et al., 2009].Therefore, a reduction or interruption of the MOW could have a critical impact both on the AMOC and on the global circulation, as well as on the global climate.Without the contribution of the saline MOW, the formation of dense water, which triggers the AMOC, would have most likely not taken place steadily in the North Atlantic [Rahmstorf, 1998;Bethoux et al., 1999].
[5] All these studies have analyzed the history and impact of the MOW on the Atlantic circulation after the end of Messinian salinity crisis (MSC), when the Atlantic-Mediterranean connections were reestablished through the Strait of Gibraltar.However, little is known about the impact of the MOW during the Messinian, when the connections were through the Betic and Rifian Corridors [van der Laan et al., 2012] or after the cessation of the MOW due to the closure of these corridors.Keigwin et al. [1987] questioned that the MSC had any effect on the deep circulation in the North Atlantic.On the contrary, Zhang and Scott [1996] reported the presence of the MOW at the northeastern Atlantic Ocean, at least reaching 50 N of latitude, during the Messinian.Moreover, Pb and Nd isotopic studies also pointed out to the influence of the MOW in the NE Atlantic during the Messinian [Abouchami et al., 1999].Apart from these works deciphering the influence of the MOW in the distant Atlantic, no study has focused on the areas close to the Atlantic-Mediterranean connections in the Betics.
[7] The Guadalquivir Basin, located in the south of the Iberian Peninsula, represents the Atlantic side of the Betic Corridors that extended through southern Spain during the early late Miocene [Benson et al., 1991;Martín et al., 2001Martín et al., , 2009;;Braga et al., 2002].The last active Betic gateway was the Guadalhorce Corridor, which controlled the Messinian pre-evaporitic circulation in the western Mediterranean, and allowed the MOW to enter the Atlantic Ocean [Martín et al., 2001].The Guadalhorce Corridor was a NW-SE trending strait with an estimated maximum width of 5 km and maximum water depth of 120 m [Martín et al., 2001].The corridor was filled by sediment displaying huge unidirectional sedimentary structures indicating Mediterranean waters flowing out into the Atlantic at estimated current velocities of about 1.0-1.5 m/s [Martín et al., 2001].This water mass circulation is consistent with the siphon model of Benson et al. [1991] stating that prior to the closure of the Betic Corridors, the water exchange between the Mediterranean and the Atlantic during the Messinian was characterized by Atlantic inflow through the Rifian Corridors and MOW through the Guadalhorce Corridor [Benson et al., 1991;Martín et al., 2001].After the closure of the Guadalhorce Corridor in the early Messinian [Martín et al., 2001], MOW was interrupted and circulation was restricted to the Rifian Corridors [Esteban et al., 1996].Later the closure of the Rifian Corridors in the late Messinian [Krijgsman et al., 1999a] caused the isolation of the Mediterranean Sea and, consequently, a hydrographical deficit that triggered the onset of the MSC with deposition of extensive evaporites in the central and deeper parts of the Mediterranean [Hsü et al., 1973[Hsü et al., , 1977]].
[8] In this study, we examine the impact of the MOW during the Messinian close to the Betic corridors.We analyze foraminiferal stable O and C isotope composition in the Montemayor-1 core (SW Spain) (Figure 2).The core site is located close to the Guadalhorce Corridor, the last Betic corridor to be closed [Martín et al., 2001[Martín et al., , 2009]], and shows a continuous Messinian record accurately dated by magnetobiostratigraphic methods [Larrasoaña et al., 2008].The Guadalquivir Basin was well connected with the Atlantic Ocean during the MSC, so there was neither desiccation nor evaporite deposition.Furthermore, according to the siphon model of Benson et al. [1991], Mediterranean outflow took place only throughout the Betic Corridors.Therefore, the location of the core is exceptional to study the effect of the MOW on the eastern Atlantic Ocean circulation during the Messinian, as well as its possible impact on global climate changes.
[9] The main aims of this study are: 1) to assess the impact of the MOW on the AMOC during the Messinian; and 2) to precisely date the closure of the Guadalhorce Corridor.
[11] The Guadalquivir foreland basin was formed as a consequence of the Betic Cordillera compressional overthrusting during the early middle Miocene [Sanz de Galdeano and Vera, 1992;Riaza and Martínez del Olmo, 1996;Sanz de Galdeano and Rodríguez-Fernández, 1996;Martín et al., 2009;Braga et al., 2010].During the Serravallian, the Atlantic-Mediterranean connection started to be restricted in the northeastern edge of the Guadalquivir Basin, in the Prebetic Domain of the Betic Cordillera [Aguirre et al., 2007;Martín et al., 2009;Braga et al., 2010].The progressive tectonic uplifting of the Betic mountain chain led to a progressive closure of this seaway, originating the so-called North Betic Strait during the latest middle Miocene-earliest late Miocene (topmost Serravallian-earliest Tortonian) [Aguirre et al., 2007;Martín et al., 2009;Braga et al., 2010].The final closure of the North Betic Strait took place during the early Tortonian [Sierro et al., 1996;Martín et al., 2009;Braga et al., 2010] and the Guadalquivir Basin was established as a wide, marine embayment only opened to the Atlantic Ocean [Martín et al., 2009].
[12] After the cessation of the North Betic Strait, other Betic gateways connected the Atlantic and the Mediterranean through the Guadalquivir Basin.They were progressively closed during the late Miocene.In the late Tortonian, the Dehesas de Guadix Corridor and the Granada Basin were the main Atlantic-Mediterranean connections [Esteban et al., 1996;Braga et al., 2003;Betzler et al., 2006;Martín et al., 2009].After their closure, the Guadalhorce Corridor was the only connection during the earliest Messinian [Martín et al., 2001].This last Betic Corridor became closed in the early Messinian [Martín et al., 2001].Since its closure, the Rifian Corridors were the unique Atlantic-Mediterranean gateways [Esteban et al., 1996].

Montemayor-1 Core
[13] The studied material is the Montemayor-1 core, a 260 m long core that has been drilled in the northwestern margin of the lower Guadalquivir Basin close to Moguer (SW Spain) (Figures 2 and 3).This core includes marine sediments that can be divided into four lithostratigraphic units [see a detail description in Pérez-Asensio et al., 2012] (Figure 3): the Niebla Formation (Tortonian), the Arcillas de Gibraleón Formation (latest Tortonian-Messinian), the Arenas de Huelva Formation (early Pliocene), and the Arenas de Bonares Formation (late Pliocene-Pleistocene).
[14] In this study, we analyzed an interval of 70 m, from 240 to 170 m (from 6.67 Ma to 5.7 Ma according to the age model.See below), including Messinian sediments from the Arcillas de Gibraleón Formation (Figure 3).In this interval, a total of 132 samples were collected with a sampling interval of 0.5 m.Samples were wet sieved over a 63 mm mesh and dried out in an oven at 40 C. A representative split was dry-sieved over a 125 mm mesh to estimate the planktonic/benthic ratio (P/B ratio henceforth), calculated as [P/(P + B)] as a proxy for relative sea level change.

Stable Isotope Analyses
[15] Stable isotope analyses (d 18 O and d 13 C) were performed on 10 individuals of Cibicidoides pachydermus for benthic foraminifera and 20 individuals of Globigerina bulloides for planktonic foraminifera separated from the size fraction >125 mm.Foraminiferal shells were ultrasonically cleaned, and washed with distilled water prior to the analyses.The isotopic analyses were performed on a Finnigan MAT 251 mass spectrometer connected to a Kiel I (prototype) Figure 3. Age model, estimation of the sedimentation rate (estimated in cm/kyr) and magnetobiostratigraphic framework for the Montemayor-1 core.Magnetostratigraphy is based on the ATNTS2004 [Lourens et al., 2004].Type 1 samples have higher quality than types 2 and 3 samples.CrRM is the characteristic remanent magnetization.Biostratigraphy has been established using the planktonic foraminiferal events (PF events) of Sierro et al. [1993]  preparation device for carbonates at the Leibniz-Laboratory for Radiometric Dating and Isotope Research, Kiel, Germany.Results are given in d-notation in per mil, and are reported on the Vienna Pee Dee belemnite (VPDB) scale.The VPDB scale is defined by a certain value of the National Bureau of Standards (NBS) carbonate standard NBS-19.On the basis of the international and lab-internal standard material, the analytical reproducibility is <AE0.05‰for d 13 C, and <AE0.07‰ for d 18 O.

Spectral Analyses
[16] Spectral analysis was performed in order to identify the nature and significance of periodic changes in the benthic d 18 O record.The analysis was carried out in the time domain using the software PAST [Hammer et al., 2001] and the REDFIT procedure of Schulz and Mudelsee [2002].This procedure allows assessing the spectral analysis with unevenly spaced samples.Spectral peaks over the 95% confidence interval (CI) were considered significant.

Age Model
[17] The age model of the Montemayor-1 core was established using a combination of paleomagnetism, biostratigraphy, and stable oxygen isotope stratigraphy (Figure 3).The reversed chron C3r is almost continuously recorded since a discontinuity is detected close to the boundary between chrons C3r and C3n, which prevent us from using this magnetostratigraphic reversal boundary as a chronological datum in the upper part of the core [Pérez-Asensio et al., 2012].
[18] To complete the age model above the normal magnetic chron C3An.1n, we have used the glacial stage TG 22 as a tie point, which was astronomically calibrated at 5.79 Ma [Krijgsman et al., 2004].We identified this glacial stage by means of stable oxygen isotope stratigraphy (Figure 4).The benthic oxygen isotope record shows two distinctly pronounced "paired" glacial peak stages that we identify as the glacial stages TG 20 and TG 22 (Figure 4) according to the nomenclature established by Shackleton et al. [1995].These two stages are easily identifiable because they are the most pronounced "paired" glacial peaks at the end of a progressively increasing trend in the d 18 O record along the chron C3r (Figure 4).Both the isotopic trend and the two "paired" glacial maxima are evident and have been observed in other cores at global scale including the Pacific Ocean [Shackleton et al., 1995], the Atlantic Ocean [Hodell et al., 2001;Vidal et al., 2002], as well as in sediments from the Rifian Corridors [Hodell et al., 1994].
[19] Using the TG 22 as a tie point in the chronological framework, the estimated age for the TG 20 based on the reconstructed sedimentation rate (14.8 cm/kyr) is 5.75 Ma (Figures 3 and 5).This age estimation is coincident with the age estimated by Krijgsman et al. [2004] using astronomical tuning.
[20] Further to the TG 20 and TG 22, we have identified the rest of the glacial stages following the nomenclature of Shackleton et al. [1995].These authors found that the benthic O isotope record was controlled by 41-kyr cycles related to orbital obliquity.Hodell et al. [1994] also showed that the benthic isotope record from the Salé Briqueterie core at the Rifian Corridors reflects obliquity induced changes.
The benthic oxygen isotope record from the Montemayor-1 core is also mainly controlled by 41-kyr cycles related to orbital obliquity (Figure 6).
[21] Since the TG glacial stages of Shackleton et al. [1995] were related to obliquity, we use the same methodology to identify the rest of the TG stages in the chron C3r.Further, in order to confirm the reliability of the identification of TG 20 and TG 22 in the Montemayor-1 core, we plotted the benthic isotope record versus age (Ma), including the TG 22 datum (5.79 Ma), and counted obliquity cycles backward (Figure 5).Comparing the benthic oxygen isotope record from the Montemayor-1 core with its obliquity component (Figure 5) we identified up to stage TG 32 for chron C3r.However, Hodell et al. [1994] and van der Laan et al.
[2005] recognized one extra obliquity cycle, up to TG 34 stage.The discrepancy in just one obliquity cycle (1 glacial stage and 1 interglacial stage) can be due to the fact that Hodell et al. [1994] identified glacial stages versus depth instead of time.This could result in counting cycles of different periodicities.On the other hand, van der Laan et al.
[2005], considered the glacial stage TG 30, which is a precession-related signal, at the same level of the rest of the obliquity-related cycles.

Results
[22] The benthic oxygen isotope record shows stable values around 1‰ before 6.35 Ma, and then it decreases reaching a minimum of À0.93‰ at 6.18 Ma (Figure 7).After 6.18 Ma, it exhibits fluctuations with a trend toward heavier values reaching the maximum values at 5.79 Ma (TG 22), and then it decreases toward lighter values (Figure 7).The benthic stable O isotope curves of the Montemayor-1 core and site 1085 [Vidal et al., 2002] reveal different trends for the time interval before 6.18 Ma (Figure 5).On the contrary, both curves show a parallel trend with similar fluctuations after 6.18 Ma.In this interval, heavy values of d 18 O from the Montemayor-1 core can be easily matched with glacial peaks from site 1085 (Figure 5).
[23] The planktonic oxygen isotope record exhibits a fluctuating trend with values ≤0‰ before 6.18 Ma and values >0‰ after this age (Figure 7).Both the benthic and planktonic O isotopic curves follow different trends before 6.18 Ma.After 6.18 Ma, the planktonic O isotopic curve parallels that of the benthic one except for some low values around 5.79 Ma (Figure 7).The similar trend of both records is confirmed by a statistically significant positive correlation (r = 0.415).The Dd 18 O benthic-planktonic (Dd 18 O b-p ) curve reflects these differences before and after 6.18 Ma between the benthic and planktonic O isotope values, being higher before 6.18 Ma (Figure 7).The difference reaches average values close to 0 at around 6.2 Ma (Figure 7).
[24] The planktonic carbon isotope record exhibits significant fluctuations with average values around À0.8‰ before 6.18 Ma (Figure 7).After 6.18 Ma, the d 13 C decreases reaching its lowest values from 6.05 to 5.85 Ma.Then, it drastically increases from average values of À1‰ to 0‰ at 5.85 Ma, and remains with relatively high average values around 0 until 5.77 Ma.At this age, another carbon shift took place recovering average values of À1‰ (Figure 7).The benthic carbon isotope record shows fluctuations with average values around 0.4‰ before 6.18 Ma.After 6.18 Ma, the d 13 C decreases reaching its lowest values from 6 to 5.9 Ma.Finally, it increases and remains with relatively high values around 0.4‰ from 5.9 to 5.7 Ma (Figure 7).

Discussion
6.1.Identification of the MOW in the Northeastern Atlantic and Age of the Closure of the Guadalhorce Corridor 6.1.1.MOW Presence Before 6.18 Ma [25] The benthic O isotopic values from the Montemayor-1 core depart from the fluctuating global trend, based on Vidal et al. [2002], from the base of the studied interval to 6.18 Ma (Figure 5).The presence of a dense, i.e., highly saline, bottom water mass affecting the study area might account for the observed deviation toward heavier values in the O isotope record.
[26] The Montemayor-1 core site is located adjacent to the Guadalhorce Corridor, the last Atlantic-Mediterranean gateway providing a connection between the Mediterranean and Atlantic oceans through the Betic Cordillera prior to the onset of the MSC (Figure 8a).Based on the vicinity of the core site to this gateway it appears most likely that the benthic habitats of the study area were bathed by the highly saline MOW prior to 6.18 Ma (Figure 8a).
[27] In the Montemayor-1 core, benthic O values remain approximately constant at around 1‰ before 6.35 Ma.This suggests a more or less constant MOW flux as indicated by the scarcity of reactivation surfaces in the sedimentary structures of the Guadalhorce Corridor [Martín et al., 2001].The benthic stable O isotope values of 1‰ reflect the density of the MOW that is the product of mixing between the exported Mediterranean water and the ambient Atlantic water.This value is in the range of recent benthic stable O isotope values of the upper core of the MOW, which ranges between 1 and 2‰ [Rogerson et al., 2011].Before 6.35 Ma the exchange through the Guadalhorce Corridor is expected to be maximal because of the increased density between the MOW and Atlantic surface waters (Dd 18 O b-p ) and the high sea level indicated by the P/B ratios (Figure 7).The benthic C isotope record could partially reflect MOW activity.Using the modern analogue, high values of benthic d 13 C are associated to MOW because of the low residence time of this water mass [Vergnaud-Grazzini, 1983;Schönfeld and Zahn, 2000;Raddatz et al., 2011;Rogerson et al., 2011].Thus, relatively high benthic C isotopic values around 6.4-6.5 Ma might be the result of increased MOW presence (Figure 7).
[28] The presence of the MOW before 6.18 Ma is also supported by the paleobathymetry of the study area during this period.At the present day, the MOW entering the Atlantic flows along the western Iberian continental slope between 400 and 1,500 m water depth [Schönfeld and Zahn, 2000;Llave et al., 2006].In the study area, the presence of R-mode Anomalinoides flinti assemblage and Planulina ariminensis before 6.18 Ma (below 227.5 m core depth in Pérez-Asensio et al., 2012) indicates that the current flowed along the middle and upper slope [Pérez-Asensio et al., 2012].This is within the depth range of MOW flow at the present day.
[29] A comparison between benthic and planktonic O isotopic records from the Montemayor-1 core offers another indication of the MOW presence before 6.18 Ma.The decoupling and high difference between the benthic and planktonic O isotopic signals before 6.18 Ma (Figure 7) are indicative of a two-layer water column with the presence of MOW at the seafloor.Moreover, the lowest sedimentation rates are estimated before 6.18 Ma (Figure 3).Thus, winnowing by the MOW might most likely account for this result.6.1.2.Response of the MOW Plume to the Restriction of the Guadalhorce Corridor (6.35 to 6.18 Ma) [30] Between 6.35 and 6.18 Ma, the benthic O isotopic record underwent a significant negative excursion of about 2‰, from 1‰ to a minimum of À0.93‰ (Figures 5 and 7).This noteworthy isotopic shift suggests a MOW product with a lighter O isotopic signature.The ambient Atlantic water has a lower density than the MOW [Price and O'Neill-Baringer, 1994].Hence, the local minimum in the benthic d 18 O from 6.35 to 6.18 Ma could be explained by a higher proportion of the lighter ambient Atlantic water in the final MOW product as occurs at present day outside the upper core of the MOW plume [Rogerson et al., 2011].In addition, this fact could reflect a gradual reduction in the exchange through the Guadalhorce Corridor as it is indicated by the reduction of density between the MOW and the Atlantic surface waters (Dd 18 O b-p ) and relatively low sea level (Figure 7).The reduced exchange from 6.35 to 6.18 Ma is also reflected by the decrease in the benthic C isotopes (Figure 7).Similarly, the diminished outflow of Mediterranean waters in the Rifian Corridors is indicated by a negative excursion in the planktonic C isotopes at 6.0 Ma [van der Laan et al., 2012].
[31] Alternatively, the local minimum in the benthic d 18 O might reflect the isotopic signature of less dense Atlantic waters.Increased salinity of Mediterranean source waters due to the restriction of the Guadalhorce Corridor could produce a denser MOW that mixes faster with ambient Atlantic waters reducing its density and consequently shoaling on the slope [Rogerson et al., 2012].Therefore, the location of the Montemayor-1 core would be beneath the MOW between 6.35 and 6.18 Ma.Our data from the Montemayor-1 core do not allow us to rule out any of these two alternative explanations.6.1.3.Cessation of the MOW at 6.18 Ma [32] After 6.18 Ma, the benthic d 18 O parallels that of the global benthic O record, with main glacial stages easily distinguishable in both records (Figure 5).Therefore, benthic O isotopic values at the study area were primarily controlled by global glacioeustatic fluctuations.Most likely, the striking change in the d 18 O record at 6.18 Ma can be linked with the cessation of the dense MOW influence in the study area (Figure 8b).
[33] A comparison between benthic and planktonic O isotopic records from the Montemayor-1 core also supports the interruption of the MOW after 6.18 Ma.Benthic and planktonic O isotopic values covary after 6.18 Ma indicating enhanced vertical mixing and absence of MOW.This is also supported by low Dd 18 O b-p and positive statistical correlation between benthic and planktonic O isotopic values (r = 0.415) (Figure 7).The covariation and the decreasing difference between planktonic and benthic O isotopes might have also been favored by the shallowing-upward trend that reduces the difference between benthic and planktonic isotopic values.
[34] Moreover, sedimentation rate could be affected by the interruption of the MOW at 6.18 Ma.After this date, sedimentation rates drastically increases from 6.3 to 14.8 cm/kyr (Figure 3) likely due to the interruption of the MOW and the resulting progradation of depositional systems along the axis of the Guadalquivir Basin [Sierro et al., 1996] (Figure 8b).
[35] The closure of this corridor at 6.18 Ma caused the final interruption of MOW.The observed age for the end of this Betic gateway is consistent with planktonic foraminifera recorded in the sediments filling the Guadalhorce Corridor that indicate that the strait was open at least from 7.2 to 6.3 Ma [Martín et al., 2001].Further, it is also coincident with the first land mammal exchange between Africa and the Iberian Peninsula prior to the MSC taking place at 6.1 Ma [Garcés et al., 1998].Therefore, according to our data, the closure of the Guadalhorce Corridor (6.18 Ma) predates by about 220-kyr the onset of the MSC, dated at 5.96 AE 0.02 Ma by Krijgsman et al. [1999b].

Impact of the MOW on the Eastern North Atlantic Ocean Circulation
[36] The MOW is essential to maintain the meridional overturning circulation (AMOC) because it increases the North Atlantic density gradient and can restart the AMOC after its collapse [Rogerson et al., 2012].It appears likely that any fluctuation of the MOW can alter the AMOC.During the Messinian period, the increase of the NADW presence in the South Atlantic Ocean between 6.6 and 6.0 Ma [Billups, 2002] suggests an enhanced NADW formation due to the suspected salinity increase of the MOW before the closure of the Atlantic-Mediterranean gateways.During this interval, MOW affected our study area before the closure of the Guadalhorce Corridor at 6.18 Ma (Figures 5 and 7).It has to be mentioned however, that the increase in NADW flow in the South Atlantic Ocean could have also been caused by the shoaling of the Central American Seaway (Panamanian Seaway) during the middlelate Miocene [Billups, 2002;Nisancioglu et al., 2003;Butzin et al., 2011].
[37] The MOW is predominantly fed by the LIW [Bryden and Stommel, 1984].During the Quaternary, cold stadials have been associated with production of a denser LIW [Cacho et al., 2000] that enhanced the current activity and depth of settling of the MOW [Schönfeld and Zahn, 2000;Rogerson et al., 2005].Similarly, during the late Miocene global cooling [Zachos et al., 2001;Murphy et al., 2009] a denser LIW could have enhanced the MOW via increasing its buoyancy loss, which is consistent with the observed increase of the NADW formation and AMOC intensification during this time interval.
[38] The final closure of the Guadalhorce Corridor might have produced a dramatic reduction of the MOW, which then could only flow through the Rifian Corridors until their final closure at around 6.0 Ma [Krijgsman et al., 1999a].This led to a diminished NADW production as is recorded between 6.0 and 5.5 Ma in the South Atlantic [Billups, 2002].Concomitantly, the cessation of the MOW increased the freshwater input in the North Atlantic reducing or interrupting the NADW formation [Rahmstorf, 1998].This might have had a critical impact on the AMOC, global thermohaline circulation and global climate [Bethoux et al., 1999].Specifically, the reduction or interruption of the NADW by MOW cessation would reduce the AMOC leading to northern hemisphere cooling [Clark et al., 2002].This is supported by development of northern hemisphere ice sheets in Greenland margin [Fronval and Jansen, 1996;Thiede et al., 1998].
[39] Finally, with the restoration of the Atlantic-Mediterranean connections through the Strait of Gibraltar, the MOW enhanced the NADW formation as it is indicated by the rise in NADW flow around the Miocene-Pliocene boundary in the South Atlantic and the western equatorial Atlantic Ocean [King et al., 1997;Billups, 2002].

Conclusions
[40] The Messinian stable isotope records from the Montemayor-1 core located in the lower Guadalquivir Basin (SW Spain) allow the accurate dating of the closure of the Guadalhorce Corridor, the last Betic Corridor connecting the Atlantic and the Mediterranean before the MSC, at 6.18 Ma.
[41] Before the closure of the Guadalhorce Corridor, the export of highly saline intermediate waters as [42] Acknowledgments.We acknowledge the comments provided by reviewers Michael Rogerson and Marit-Solveig Seidenkrantz which substantially improved the quality of an earlier version of the manuscript.This paper is part of the research projects CGL2010-20857 and CGL2009-11539/BTE of the Ministerio de Ciencia e Innovación of Spain and the Research Group RNM-190 of the Junta de Andalucía.J.N.P.A. has been funded by a research F.P.U.scholarship (ref.AP2007-00345) provided by the Ministerio de Educación of Spain.We are grateful to N. Andersen (Leibniz-Laboratory, Kiel, Germany) for stable isotope analysis.We also thank L. Vidal who provided the stable isotope data of benthic foraminifera from the ODP Site 1085 (SE Atlantic).

Figure 1 .
Figure 1.(a) Present salinity and circulation patterns at the Strait of Gibraltar and in the Mediterranean (based on Wüst [1961]).The dense and saline Mediterranean intermediate and deep waters form a bottom outflow, while lower salinity surface waters from the Atlantic enters the Mediterranean representing an anti-estuarine circulation pattern.(b) Present-day general circulation pattern at the eastern North Atlantic and Mediterranean Sea (based on Hernández-Molina et al. [2011] and Pinardi and Masetti [2000]).Mediterranean Outflow Water (MOW), North Atlantic Deep Water (NADW), Levantine Intermediate Water (LIW), Western Mediterranean Deep Water (WMDW) and WMDW formation are indicated.

Figure 2 .
Figure 2. Geological setting of the lower Guadalquivir Basin.Asterisk indicates the location of the Montemayor-1 core.

Figure 4 .
Figure 4. Benthic (Cibicidoides pachydermus) stable oxygen isotope record versus core depth (gray line) and three-point running average of the benthic oxygen isotope record (black line) of the Montemayor-1 core.Glacial stages TG 20 and TG 22 defined by Shackleton et al. [1995] in benthic oxygen isotope record of the ODP site 846 are indicated.

Figure 5 .
Figure 5. Benthic (Cibicidoides pachydermus) stable oxygen isotope record versus time (gray lines) and three-point running average of the benthic oxygen isotope record (black lines) from the Site 1085 [Vidal et al., 2002] and Montemayor-1 core.The obliquity component of the benthic isotope record is extracted by sinusoidal curve fitting with a period of 41 kyr.The onset of the Messinian salinity crisis (MSC) is indicated (vertical gray bar in the lower diagram).The vertical black line at 6.18 Ma marks the end of the Atlantic-Mediterranean connection through the Guadalhorce Corridor.Glacial stages from TG 20 to TG 32 defined by Shackleton et al. [1995] and from C3An.d 18 O.2 to C3An.d 18 O.8 according to the nomenclature used in Hodell et al. [1994] are indicated.The two vertical black solid lines mark glacial stages TG 20 and 22 from the Montemayor-1 core and Site 1085.Vertical gray dashed lines mark glacial stages and the obliquity component of the benthic oxygen isotope record of the Montemayor-1 core.

Figure 6 .
Figure 6.Power spectrum of the benthic (Cibicidoides pachydermus) stable oxygen isotope record of the Montemayor-1 core with main orbital periodicities indicated in kiloyears.The 95% confidence interval (CI) is indicated.

Figure 7 .
Figure 7. Benthic (Cibicidoides pachydermus) and planktonic (Globigerina bulloides) d 18 O and d 13 C records in ‰ VPDB, difference between d 18 O and d 13 C values of benthic and planktonic foraminifera (Dd 18 O b-p and Dd 13 C b-p ), as well as the planktonic-benthic ratio (P/B ratio) versus time (gray lines) and three-point running average (black lines) of the Montemayor-1 core.The vertical dashed line at 6.18 Ma marks the end of the Atlantic-Mediterranean connection through the Guadalhorce Corridor.

Figure 8 .
Figure 8. Paleogeographic and paleoceanographic evolution of the lower Guadalquivir Basin during the Messinian (based on Martín et al., 2009): (a) situation before 6.18 Ma, when the Mediterranean Outflow Water (MOW) reaches the studied core (asterisk) and (b) situation after 6.18 Ma, when the MOW was interrupted.The black thick arrow marks the progradation of the main depositional systems along the axis of the Guadalquivir Basin.
MOW contributed to increased NADW formation in the North Atlantic Ocean.During this period, constant benthic d 18 O values that depart from the global glacioeustatic trend, high Dd 18 O b-p , and low sedimentation rates indicate a strong water stratification and bottom water winnowing due to the MOW flow.After the closure of the corridor, benthic d 18 O values, which parallel that of the global glacioeustatic curve, low Dd 18 O b- p , and high sedimentation rates suggest an improved vertical mixing of the water column and interruption of MOW.Changes in the stable O isotope composition of the subsurface North Atlantic water masses during this time interval are primarily controlled by glacioeustatic processes.The cessation of the MOW might have resulted in decreased NADW formation between 6.0 and 5.5 Ma weakening the AMOC and promoting northern hemisphere cooling.Then, the restoration of the MOW due to the opening of the Strait of Gibraltar around the Miocene-Pliocene boundary would have favored an enhanced NADW formation.
and first occurrence of Globorotalia puncticulata.Position of glacial stage TG 22 is indicated.Question marks show uncertainties in the chronology and sedimentation rate.