Evolution of neodymium isotopic signature of seawater during the Late Cretaceous: Implications for intermediate and deep circulation

Neodymium isotopic compositions (εNd) have been largely used for the last fifty years as a tracer of past ocean circulation, and more intensively during the last decade to investigate ocean circulation during the Cretaceous period. Despite a growing set of data, circulation patterns still remain unclear during this period. In particular, the identification of the deep-water masses and their spatial extension within the different oceanic basins are poorly constrained. In this study we present new deep-water εNd data inferred from the Nd isotope composition of fish remains and Fe-Mn oxyhydroxide coatings on foraminifera tests, along with new εNd data of residual (partly detrital) fraction recovered from DSDP sites 152 (Nicaraguan Rise), 258 (Naturaliste Plateau), 323 (Bellinghausen Abyssal Plain), and ODP sites 690 (Maud Rise) and 700 (East Georgia Basin, South Atlantic). The presence of abundant authigenic minerals in the sediments at sites 152 and 690 detected by XRD analyses may explain both middle rare earth element enrichments in the spectra of the residual fraction and the evolution of residual fraction εNd that mirror that of the bottom waters at the two sites. The results point towards a close correspondence between the bottom water εNd values of sites 258 and 700 from the late Turonian to the Santonian. Since the deep-water Nd isotope values at these two sites are also similar to those at other proto-Indian sites, we propose the existence of a common intermediate to deep-water water mass as early as the mid-Cretaceous. The water mass would have extended from the central part of the South Atlantic to the eastern part of proto-Indian ocean sites, beyond the Kerguelen Plateau. Furthermore, data from south and north of the Rio Grande Rise-Walvis Ridge complex (sites 700 and 530) are indistinguishable from the Turonian to Campanian, suggesting a common water mass since the Turonian at least. This view is supported by a reconstruction of the Rio Grande Rise-Walvis Ridge complex during the Turonian, highlighting the likely existence of a deep breach between the Rio Grande Rise and the proto-Walvis Ridge at that time. Thus deep-water circulation may have been possible between the different austral basins as early as the Turonian, despite the presence of potential oceanic barriers. Comparison of new seawater and residue εNd data on Nicaraguan Rise suggest a westward circulation of intermediate waters through the Caribbean Seaway during the Maastrichtian and Paleocene from the North Atlantic to the Pacific. This westward circulation reduced the Pacific water influence in the Atlantic, and was likely responsible for more uniform, less radiogenic εNd values in the North Atlantic after 80 Ma. Additionally, our data document an increasing trend observed in several oceanic basins during the Maastrichtian and the Paleocene, which is more pronounced in the North Pacific. Although the origin of this increase still remains unclear, it might be explained by an increase in the contribution of radiogenic material to upper ocean waters in the northern Pacific. By sinking to depth, these waters may have redistributed to some extent more radiogenic signatures to other ocean basins through deep-water exchanges.

Characterizing, modelling and understanding the climate variability of the deep water formation in the North-Western Mediterranean Sea

Observing, modelling and understanding the climate-scale variability of the deep water formation (DWF) in the North-Western Mediterranean Sea remains today very challenging. In this study, we first characterize the interannual variability of this phenomenon by a thorough reanalysis of observations in order to establish reference time series. These quantitative indicators include 31 observed years for the yearly maximum mixed layer depth over the period 1980–2013 and a detailed multi-indicator description of the period 2007–2013. Then a 1980–2013 hindcast simulation is performed with a fully-coupled regional climate system model including the high-resolution representation of the regional atmosphere, ocean, land-surface and rivers. The simulation reproduces quantitatively well the mean behaviour and the large interannual variability of the DWF phenomenon. The model shows convection deeper than 1000 m in 2/3 of the modelled winters, a mean DWF rate equal to 0.35 Sv with maximum values of 1.7 (resp. 1.6) Sv in 2013 (resp. 2005). Using the model results, the winter-integrated buoyancy loss over the Gulf of Lions is identified as the primary driving factor of the DWF interannual variability and explains, alone, around 50 % of its variance. It is itself explained by the occurrence of few stormy days during winter. At daily scale, the Atlantic ridge weather regime is identified as favourable to strong buoyancy losses and therefore DWF, whereas the positive phase of the North Atlantic oscillation is unfavourable. The driving role of the vertical stratification in autumn, a measure of the water column inhibition to mixing, has also been analyzed. Combining both driving factors allows to explain more than 70 % of the interannual variance of the phenomenon and in particular the occurrence of the five strongest convective years of the model (1981, 1999, 2005, 2009, 2013). The model simulates qualitatively well the trends in the deep waters (warming, saltening, increase in the dense water volume, increase in the bottom water density) despite an underestimation of the salinity and density trends. These deep trends come from a heat and salt accumulation during the 1980s and the 1990s in the surface and intermediate layers of the Gulf of Lions before being transferred stepwise towards the deep layers when very convective years occur in 1999 and later. The salinity increase in the near Atlantic Ocean surface layers seems to be the external forcing that finally leads to these deep trends. In the future, our results may allow to better understand the behaviour of the DWF phenomenon in Mediterranean Sea simulations in hindcast, forecast, reanalysis or future climate change scenario modes. The robustness of the obtained results must be however confirmed in multi-model studies.