On the seasonal variations of salinity of the tropical Atlantic mixed layer

The physical processes controlling the mixed layer salinity (MLS) seasonal budget in the tropical Atlantic Ocean are investigated using a regional configuration of an ocean general circulation model. The analysis reveals that the MLS cycle is generally weak in comparison of individual physical processes entering in the budget because of strong compensation. In evaporative regions, around the surface salinity maxima, the ocean acts to freshen the mixed layer against the action of evaporation. Poleward of the southern SSS maxima, the freshening is ensured by geostrophic advection, the vertical salinity diffusion and, during winter, a dominant contribution of the convective entrainment. On the equatorward flanks of the SSS maxima, Ekman transport mainly contributes to supply freshwater from ITCZ regions while vertical salinity diffusion adds on the effect of evaporation. All these terms are phase locked through the effect of the wind. Under the seasonal march of the ITCZ and in coastal areas affected by river (7°S:15°N), the upper ocean freshening by precipitations and/or runoff is attenuated by vertical salinity diffusion. In the eastern equatorial regions, seasonal cycle of wind forced surface currents advect freshwaters, which are mixed with subsurface saline water because of the strong vertical turbulent diffusion. In all these regions, the vertical diffusion presents an important contribution to the MLS budget by providing, in general, an upwelling flux of salinity. It is generally due to vertical salinity gradient and mixing due to winds. Furthermore, in the equator where the vertical shear, associated to surface horizontal currents, is developed, the diffusion depends also on the sheared flow stability.

Geochemical synthesis for the Effingen Member in boreholes at Oftringen, Gösgen and Küttigen

The Effingen Member is a low-permeability rock unit of Oxfordian age (ca. 160 Ma) that occurs across northern Switzerland. It comprises sandy calcareous marls and (argillaceous) limestones. This report describes the hydrogeochemistry, mineralogy and supporting physical properties of the Effingen Member in three boreholes in the Jura-Südfuss area: Oftringen, Gösgen and Küttigen, where it is 220–240 m thick. The top of the Effingen Member is at 420, 66 and 32 m depths at the three sites. Core materials are available from Oftringen and Gösgen, whereas information from Küttigen is limited to cuttings, in-situ hydrogeological testing and geophysical logging. Hydrogeological boundaries of the Effingen Member vary between locations. Ground-water flows were identified during drilling at the top (Geissberg Member), but not at the base, of the Effingen Member at Oftringen, at the base (Hauptrogenstein Formation) of the Effingen Member at Gösgen, and in a limestone layer (Gerstenhübel unit) within the Effingen Member at Küttigen. The marls and limestones of the Effingen Member have carbonate contents of 46–91 wt.-% and clay-mineral contents of 5–37 wt.-%. Pyrite contents are up to 1.6 wt.-%, but no sulphate minerals were detected by routine analyses. Clay minerals are predominantly mixed-layer illite-smectite, illite and kaolinite, with sporadic traces of chlorite and smectite. Veins filled with calcite ± celestite occur through the Effingen Member at Oftringen but not at Gösgen or Küttigen. They formed at 50–70 ºC from externally derived fluids, probably of Miocene age. Water contents are 0.7–4.2 wt.-%, corresponding to a water-loss porosity range of 1.9–10.8 vol.-%. Specific surface areas, measured by the BET method, are 2–30 m2/g, correlating with clay-mineral contents. Water activity has been measured and yielded surprisingly low values down to 0.8. These cannot be explained by pore-water salinity alone and include other effects, such as changes in the fabric due to stress release or partial saturation. Observed variations in measurements are not fully understood. Cation exchange capacity (CEC) and exchangeable cation populations have been studied by the Ni-en method. CEC, derived from the consumption of the index cation Ni, is 9–99 meq/kgrock at a solid:liquid ratio of 1, correlating with the clay-mineral content. Cation concentrations in Ni-en extract solutions are in the order Na+≥Ca2+>Mg2+>K+>Sr2+. However, the analytical results from the Ni-en extractions have additional contributions from cations originating from pore water and from mineral dissolution reactions that occurred during extraction, and it was not possible to reliably quantify these contributions. Therefore, in-situ cation populations and selectivity coefficients could not be derived. A suite of methods have been used for characterising the chemical compositions of pore waters in the Effingen Member. Advective displacement was used on one sample from each Oftringen and Gösgen and is the only method that produces results that approach complete hydrochemical compositions. Aqueous extraction was used on core samples from these two boreholes and gives data only for Cl- and, in some cases, Br-. Out-diffusion was used on core samples from Oftringen and similarly gives data for Cl- and Br- only. For both aqueous extraction and out-diffusion, reaction of the experimental water with rock affected concentrations of cations, SO42 and alkalinity in experimental solutions. Another method, centrifugation, failed to extract pore water. Stable isotope ratios (δ18O and δ2H) of pore waters in core samples from Oftringen were analysed by the diffusive exchange method and helium contents of pore water in Oftringen samples were extracted for mass spectrometric analysis by quantitative outgassing of preserved core samples. Several lines of evidence indicate that drillcore samples might not have been fully saturated when opened and subsampled in the laboratory. These include comparisons of water-loss porosities with physical porosities, water-activity measurements, and high contents of dissolved gas as inferred from ground-water samples. There is no clear proof of partial saturation and it is unclear whether this might represent in-situ conditions or is due to exsolution of gas due to the pressure release since drilling. Partial saturation would have no impact on the recalculation of pore-water compositions from aqueous extraction experiments using water-loss porosity data. The largest uncertainty in the pore-water Cl- concentrations recalculated from aqueous extraction and out-diffusion experiments is the magnitude of the anion-accessible fraction of water-loss porosity. General experience of clay-mineral rich formations suggests that the anion-accessible porosity fraction is very often about 0.5 and generally in a range of 0.3 to 0.6 and tends to be inversely correlated with clay-mineral contents. Comparisons of the Cl- concentration in pore water obtained by advective displacement with that recalculated from aqueous extraction of an adjacent core sample suggests a fraction of 0.27 for an Oftringen sample, whereas the same procedure for a Gösgen sample suggests a value of 0.64. The former value for anion-accessible porosity fraction is presumed to be unrepresentative given the local mineralogical heterogeneity at that depth. Through-diffusion experiments with HTO and 36Cl- suggest that the anion-accessible porosity fraction in the Effingen Member at Oftringen and Gösgen is around 0.5. This value is proposed as a typical average for rocks of the Effingen Member, bearing in mind that it varies on a local scale in response to the heterogeneity of lithology and pore-space architecture. The substantial uncertainties associated with the approaches to estimating anion-accessible porosity propagate into the calculated values of in-situ pore-water Cl- concentrations. On the basis of aqueous extraction experiments, and using an anion-accessible porosity fraction of 0.5, Cl- concentrations in the Effingen Member at Oftringen reach a maximum of about 14 g/L in the centre. Cl- decreases upwards and downwards from that, forming a curved depth profile. Cl- contents in the Effingen Member at Gösgen increase with depth from about 3.5 g/L to about 14 g/L at the base of the cored profile (which corresponds to the centre of the formation). Out-diffusion experiments were carried out on four samples from Oftringen, distributed through the Effingen Member. Recalculated Cl- concentrations are similar to those from aqueous extraction for 3 out of the 4 samples, and somewhat lower for one sample. Concentrations of other components, i.e. Na+, K+, Ca2+, Mg2+, Sr2+, SO42- and HCO3- cannot be obtained from the aqueous extraction and out-diffusion experimental data because of mineral dissolution and cation exchange reactions during the experiments. Pore-water pH also is not constrained by those extraction experiments. The only experimental approach to obtain complete pore-water compositions for samples from Oftringen and Gösgen is advective displacement of pore water. The sample from Oftringen used for this experiment is from 445 m depth in the upper part of the Effingen Member and gave eluate with 16.5 g/L Cl- whereas aqueous extraction from a nearby sample indicated about 9 g/L Cl-. The sample from Gösgen used for advective displacement is from 123 m depth in the centre of the Effingen Member sequence and gave eluate with about 9 g/L Cl- whereas aqueous extraction gave 11.5 g/L Cl-. In both cases the pore waters have Na-(Ca)-Cl compositions and SO42- concentrations of about 1.1 g/L. The Gösgen sample has a Br/Cl ratio similar to that of sea water, whereas this ratio is lower for the Oftringen sample. Taking account of uncertainties in the applied experimental approaches, it is reasonable to place an upper limit of ca. 20 g/L on Cl- concentration for pore water in the Effingen Member in this area. There are major discrepancies between pore-water SO42- concentrations inferred from aqueous extraction or out-diffusion experiments and those obtained from advective displacement in both the Oftringen and Gösgen cases. A general conclusion is that all or at least part of the discrepancies are attributable to perturbation of the sulphur system and enhancement of SO42- by sulphate mineral dissolution and possibly minor pyrite oxidation during aqueous extraction and out-diffusion. Therefore, data for SO42- calculated from those pore-water sampling methods are considered not to be representative of in-situ conditions. A reference pore-water composition was defined for the Effingen Member in the Jura Südfuss area. It represents the probable upper limits of Cl- contents and corresponding anion and cation concentrations that are reasonably constrained by experimental data. Except for Cl- and possibly Na+ concentrations, this composition is poorly constrained especially with respect to SO42- and Ca2+ concentrations, and pH and alkalinity. Stable isotope compositions, δ18O and δ2H, of pore waters in the Effingen Member at Oftringen plot to the right of the meteoric water line, suggesting that 18O has been enriched by water-rock exchange, which indicates that the pore waters have a long residence time. A long residence time of pore water is supported by the level of dissolved 4He that has accumulated in pore water of the Effingen Member at Oftringen. This is comparable with, or slightly higher than, the amounts of 4He in the Opalinus Clay at Benken. Ground waters were sampled from flowing zones intersected by boreholes at the three locations. The general interpretation is that pore waters and ground-water solutes may have similar origins in Mesozoic and Cenozoic brackish-marine formations waters, but ground-water solutes have been diluted rather more than pore waters by ingress of Tertiary and Quaternary meteoric waters. The available hydrochemical data for pore waters from the Effingen Member at these three locations in the Jura-Südfuss area suggest that the geochemical system evolved slowly over geological periods of time, in which diffusion was an important mechanism of solute transport. The irregularity of Cl- and δ18O profiles and spatial variability of advective ground-water flows in the Malm-Dogger system suggests that palaeohydrogeological and hydrochemical responses to changing tectonic and surface environmental conditions were complex.

Sublimation of water ice mixed with silicates and tholins: Evolution of surface texture and reflectance spectra, with implications for comets

The surfaces of many objects in the Solar System comprise substantial quantities of water ice sometimes mixed with minerals and/or organic molecules. The sublimation of the ice changes the structural and optical properties of these objects. We present laboratory data on the evolution of the structure and the visible and near-infrared spectral reflectance of icy surface analogues of cometary ices, made of water ice, complex organic matter (tholins) and silicates, as they undergo sublimation under low temperature (<-70°C) and pressure (10-⁵mbar) conditions inside the SCITEAS simulation chamber. As the water ice sublimated, we observed in situ the formation of a porous sublimation lag deposit, or sublimation mantle, at the top of the ice. This mantle is a network of filaments made of the non-volatile particles. Organics or phyllosilicates grains, able to interact via stronger inter-particulate forces than olivine grains, can form a foam-like structure having internal cohesiveness, holding olivine grains together. As this mantle builds-up, the band depths of the sub-surface water ice are attenuated until complete extinction under only few millimeters of mantle. Optically thick sublimation mantles are mainly featureless in the near infrared. The absorption bands of the minerals present in the mantle are weak, or even totally absent if minerals are mixed with organics which largely dominate the VIS–NIR reflectance spectrum. During sublimation, ejections of large fragments of mantle, triggered by the gas flow, expose ice particles to the surface. The contrast of brightness between mantled and ice-exposed areas depends on the wavelength range and the dust/ice ratio considered. We describe how the chemical nature of the non-volatiles, the size of their particles, the way they are mixed with the ice and the dust/ice mass ratio influence the texture, activity and spectro-photometric properties of the sublimation mantles. These data provide useful references for interpreting remote-sensing observations of comets and also icy satellites or trans-neptunian objects.

Projected 21st century decrease in marine productivity: a multi-model analysis

Changes in marine net primary productivity (PP) and export of particulate organic carbon (EP) are projected over the 21st century with four global coupled carbon cycle-climate models. These include representations of marine ecosystems and the carbon cycle of different structure and complexity. All four models show a decrease in global mean PP and EP between 2 and 20% by 2100 relative to preindustrial conditions, for the SRES A2 emission scenario. Two different regimes for productivity changes are consistently identified in all models. The first chain of mechanisms is dominant in the low- and mid-latitude ocean and in the North Atlantic: reduced input of macro-nutrients into the euphotic zone related to enhanced stratification, reduced mixed layer depth, and slowed circulation causes a decrease in macro-nutrient concentrations and in PP and EP. The second regime is projected for parts of the Southern Ocean: an alleviation of light and/or temperature limitation leads to an increase in PP and EP as productivity is fueled by a sustained nutrient input. A region of disagreement among the models is the Arctic, where three models project an increase in PP while one model projects a decrease. Projected changes in seasonal and interannual variability are modest in most regions. Regional model skill metrics are proposed to generate multi-model mean fields that show an improved skill in representing observation-based estimates compared to a simple multi-model average. Model results are compared to recent productivity projections with three different algorithms, usually applied to infer net primary production from satellite observations.

Solar forcing of North Atlantic surface temperature and salinity over the past millennium

There were several centennial-scale fluctuations in the climate and oceanography of the North Atlantic region over the past 1,000 years, including a period of relative cooling from about AD 1450 to 1850 known as the Little Ice Age1. These variations may be linked to changes in solar irradiance, amplified through feedbacks including the Atlantic meridional overturning circulation2. Changes in the return limb of the Atlantic meridional overturning circulation are reflected in water properties at the base of the mixed layer south of Iceland. Here we reconstruct thermocline temperature and salinity in this region from AD 818 to 1780 using paired δ18O and Mg/Ca ratio measurements of foraminifer shells from a subdecadally resolved marine sediment core. The reconstructed centennial-scale variations in hydrography correlate with variability in total solar irradiance. We find a similar correlation in a simulation of climate over the past 1,000 years. We infer that the hydrographic changes probably reflect variability in the strength of the subpolar gyre associated with changes in atmospheric circulation. Specifically, in the simulation, low solar irradiance promotes the development of frequent and persistent atmospheric blocking events, in which a quasi-stationary high-pressure system in the eastern North Atlantic modifies the flow of the westerly winds. We conclude that this process could have contributed to the consecutive cold winters documented in Europe during the Little Ice Age.

Understanding the performance of the FLake model over two African Great Lakes

The ability of the one-dimensional lake model FLake to represent the mixolimnion temperatures for tropical conditions was tested for three locations in East Africa: Lake Kivu and Lake Tanganyika's northern and southern basins. Meteorological observations from surrounding automatic weather stations were corrected and used to drive FLake, whereas a comprehensive set of water temperature profiles served to evaluate the model at each site. Careful forcing data correction and model configuration made it possible to reproduce the observed mixed layer seasonality at Lake Kivu and Lake Tanganyika (northern and southern basins), with correct representation of both the mixed layer depth and water temperatures. At Lake Kivu, mixolimnion temperatures predicted by FLake were found to be sensitive both to minimal variations in the external parameters and to small changes in the meteorological driving data, in particular wind velocity. In each case, small modifications may lead to a regime switch, from the correctly represented seasonal mixed layer deepening to either completely mixed or permanently stratified conditions from similar to 10 m downwards. In contrast, model temperatures were found to be robust close to the surface, with acceptable predictions of near-surface water temperatures even when the seasonal mixing regime is not reproduced. FLake can thus be a suitable tool to parameterise tropical lake water surface temperatures within atmospheric prediction models. Finally, FLake was used to attribute the seasonal mixing cycle at Lake Kivu to variations in the near-surface meteorological conditions. It was found that the annual mixing down to 60m during the main dry season is primarily due to enhanced lake evaporation and secondarily to the decreased incoming long wave radiation, both causing a significant heat loss from the lake surface and associated mixolimnion cooling.

Diffusion-driven transport in clayrock formations

Clay mineral-rich sedimentary formations are currently under investigation to evaluate their potential use as host formations for installation of deep underground disposal facilities for radioactive waste (e.g. Boom Clay (BE), Opalinus Clay (CH), Callovo-Oxfordian argillite (FR)). The ultimate safety of the corresponding repository concepts depends largely on the capacity of the host formation to limit the flux towards the biosphere of radionuclides (RN) contained in the waste to acceptably low levels. Data for diffusion-driven transfer in these formations shows extreme differences in the measured or modelled behaviour for various radionuclides, e. g. between halogen RN (Cl-36, I-129) and actinides (U-238,U-235, Np-237, Th-232, etc.), which result from major differences between RN of the effects on transport of two phenomena: diffusion and sorption. This paper describes recent research aimed at improving understanding of these two phenomena, focusing on the results of studies carried out during the EC Funmig IP on clayrocks from the above three formations and from the Boda formation (HU). Project results regarding phenomena governing water, cation and anion distribution and mobility in the pore volumes influenced by the negatively-charged surfaces of clay minerals show a convergence of the modelling results for behaviour at the molecular scale and descriptions based on electrical double layer models. Transport models exist which couple ion distribution relative to the clay-solution interface and differentiated diffusive characteristics. These codes are able to reproduce the main trends in behaviour observed experimentally, e.g. D-e(anion) < D-e(HTO) < D-e(cation) and D-e(anion) variations as a function of ionic strength and material density. These trends are also well-explained by models of transport through ideal porous matrices made up of a charged surface material. Experimental validation of these models is good as regards monovalent alkaline cations, in progress for divalent electrostatically-interacting cations (e.g. Sr2+) and still relatively poor for 'strongly sorbing', high K-d cations. Funmig results have clarified understanding of how clayrock mineral composition, and the corresponding organisation of mineral grain assemblages and their associated porosity, can affect mobile solute (anions, HTO) diffusion at different scales (mm to geological formation). In particular, advances made in the capacity to map clayrock mineral grain-porosity organisation at high resolution provide additional elements for understanding diffusion anisotropy and for relating diffusion characteristics measured at different scales. On the other hand, the results of studies focusing on evaluating the potential effects of heterogeneity on mobile species diffusion at the formation scale tend to show that there is a minimal effect when compared to a homogeneous property model. Finally, the results of a natural tracer-based study carried out on the Opalinus Clay formation increase confidence in the use of diffusion parameters measured on laboratory scale samples for predicting diffusion over geological time-space scales. Much effort was placed on improving understanding of coupled sorption-diffusion phenomena for sorbing cations in clayrocks. Results regarding sorption equilibrium in dispersed and compacted materials for weakly to moderately sorbing cations (Sr2+, Cs+, Co2+) tend to show that the same sorption model probably holds in both systems. It was not possible to demonstrate this for highly sorbing elements such as Eu(III) because of the extremely long times needed to reach equilibrium conditions, but there does not seem to be any clear reason why such elements should not have similar behaviour. Diffusion experiments carried out with Sr2+, Cs+ and Eu(III) on all of the clayrocks gave mixed results and tend to show that coupled diffusion-sorption migration is much more complex than expected, leading generally to greater mobility than that predicted by coupling a batch-determined K-d and Ficks law based on the diffusion behaviour of HTO. If the K-d measured on equivalent dispersed systems holds as was shown to be the case for Sr, Cs (and probably Co) for Opalinus Clay, these results indicate that these cations have a D-e value higher than HTO (up to a factor of 10 for Cs+). Results are as yet very limited for very moderate to strongly sorbing species (e.g. Co(II), Eu(III), Cu(II)) because of their very slow transfer characteristics. (C) 2011 Elsevier Ltd. All rights reserved.

Provenance discrimination of Lower Cretaceous synrift sandstones (eastern Iberian Chain Spain) Constraints from detrital modes heavy minerals and geochemistry

Petrography, geochemical whole-rock composition, and chemical analyses of tourmaline were performed in order to determine the source areas of Lower Cretaceous Mora, El Castellar, and uppermost Camarillas Formation sandstones from the Iberian Chain, Spain. Sandstones were deposited in intraplate subbasins, which are bound by plutonic and volcanic rocks of Permian, Triassic, and Jurassic age, Paleozoic metamorphic rocks, and Triassic sedimentary rocks. Modal analyses together with petrographic and cathodoluminescence observations allowed us to define three quartz-feldspathic petrofacies and recognize diagenetic processes that modified the original framework composition. Results from average restored petrofacies are: Mora petrofacies = P/F >1 and Q(r)70 F(r)22 R(r)9; El Castellar petrofacies = P/F >1 and Q(r)57 F(r)25 R(r)18; and Camarillas petrofacies = P/F ∼ zero and Q(r)64 F(r)28 R(r)7 (P—plagioclase; F—feldspar; Q—quartz; R—rock fragments; r—restored composition). Trace-element and rare earth element abundances of whole-rock analyses discriminate well between the three petrofacies based on: (1) the Rb concentration, which is indicative of the K content and reflects the amount of K-feldspar modal abundance, and (2) the relative modal abundance of heavy minerals (tourmaline, zircon, titanite, and apatite), which is reproduced by the elements hosted in the observed heavy mineral assemblage (i.e., B and Li for tourmaline; Zr, Hf, and Ta for zircon; Ti, Ta, Nb, and their rare earth elements for titanite; and P, Y, and their rare earth elements for apatite). Tourmaline chemical composition for the three petrofacies ranges from Fe-tourmaline of granitic to Mg-tourmaline of metamorphic origin. The three defined petrofacies suggest a mixed provenance from plutonic and metamorphic source rocks. However, a progressively major influence of granitic source rocks was detected from the lowermost Mora petrofacies toward the uppermost Camarillas petrofacies. This provenance trend is consistent with the uplift and erosion of the Iberian Massif, which coincided with the development of the latest Berriasian synrift regional unconformity and affected all of the Iberian intraplate basins. The uplifting stage of Iberian Massif pluton caused a significant dilution of Paleozoic metamorphic source areas, which were dominant during the sedimentation of the lowermost Mora and El Castellar petrofacies. The association of petrographic data with whole-rock geochemical compositions and tourmaline chemical analysis has proved to be useful for determining source area characteristics, their predominance, and the evolution of source rock types during the deposition of quartz-feldspathic sandstones in intraplate basins. This approach ensures that provenance interpretation is consistent with the geological context.

Interaction of corroding iron with bentonite in the ABM1 experiment at Äspö, Sweden: A microscopic approach

Bentonite and iron metals are common materials proposed for use in deep-seated geological repositories for radioactive waste. The inevitable corrosion of iron leads to interaction processes with the clay which may affect the sealing properties of the bentonite backfill. The objective of the present study was to improve our understanding of this process by studying the interface between iron and compacted bentonite in a geological repository-type setting. Samples of MX-80 bentonite samples which had been exposed to an iron source and elevated temperatures (up to 115ºC) for 2.5 y in an in situ experiment (termed ABM1) at the Äspö Hard Rock Laboratory, Sweden, were investigated by microscopic means, including scanning electron microscopy, μ-Raman spectroscopy, spatially resolved X-ray diffraction, and X-ray fluorescence. The corrosion process led to the formation of a ~100 mm thick corrosion layer containing siderite, magnetite, some goethite, and lepidocrocite mixed with the montmorillonitic clay. Most of the corroded Fe occurred within a 10 mm-thick clay layer adjacent to the corrosion layer. An average corrosion depth of the steel of 22–35 μm and an average Fe2+ diffusivity of 1–26×10–13 m2/s were estimated based on the properties of the Fe-enriched clay layer. In that layer, the corrosion-derived Fe occurred predominantly in the clay matrix. The nature of this Fe could not be identified. No indications of clay transformation or newly formed clay phases were found. A slight enrichment of Mg close to the Fe–clay contact was observed. The formation of anhydrite and gypsum, and the dissolution of some SiO2 resulting from the temperature gradient in the in situ test, were also identified. © 2014, Clay Minerals Society. All right reserved.