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.

Nickel isotope fractionation during laterite Ni ore smelting and refining: implications for tracing the sources of Ni in smelter-affected soils

Nickel isotope ratios were measured in ores, fly ash, slags and FeNi samples from two metallurgical plants located in the Goiás State, Brazil (Barro Alto, Niquelândia). This allowed investigating the mass-dependent fractionation of Ni isotopes during the Ni-laterite ore smelting and refining. Feeding material exhibits a large range of δ60Ni values (from 0.02 ± 0.10 ‰ to 0.20 ± 0.05 ‰, n=7), explained by the diversity of Ni-bearing phases, and the average of δ60Nifeeding materials was found equal to 0.08 ± 0.08‰ (2SD, n=7). Both δ60Ni values of fly ash (δ60Ni = 0.07 ± 0.07‰, n=10) and final FeNi produced (0.05 ± 0.02 ‰, n=2) were not significantly different from the feeding materials ones. These values are consistent with the very high production yield of the factories. However, smelting slags present the heaviest δ60Ni values of all the smelter samples, with δ60Ni ranging from 0.11 ± 0.05 ‰ to 0.27 ± 0.05 ‰ (n=8). Soils were also collected near and far from the Niquelândia metallurgical plant, to evaluate the potential of Ni isotopes for tracing the natural vs anthropogenic Ni in soils. The Ni isotopic composition of the non-impacted topsoils developed on ultramafic rocks ranges from -0.26 ± 0.09 ‰ to -0.04 ± 0.05 ‰ (n=20). On the contrary, the Ni isotopic composition of the non-ultramafic topsoils, collected close to the plant, exhibit a large variation of δ60Ni, ranging from -0.19 ± 0.13 ‰ up to 0.10 ± 0.05 ‰ (n=4). This slight but significant enrichment in heavy isotopes highlight the potential impact of smelting activity in the surrounding area, as well as the potential of Ni isotopes for discerning anthropogenic samples (heavier δ60Ni values) from natural ones (lighter δ60Ni values). However, given the global range of published δ60Ni values (from -1.03 to 2.5 ‰) and more particularly those associated to natural weathering of ultramafic rocks (from -0.61 to 0.32‰), the use of Ni isotopes for tracing environmental contamination from smelters will remain challenging.