Resolving diffusion in clay minerals at different time scales: Combination of experimental and modeling approaches

Since no single experimental or modeling technique provides data that allow a description of transport processes in clays and clay minerals at all relevant scales, several complementary approaches have to be combined to understand and explain the interplay between transport relevant phenomena. In this paper molecular dynamics simulations (MD) were used to investigate the mobility of water in the interlayer of montmorillonite (Mt), and to estimate the influence of mineral surfaces and interlayer ions on the water diffusion. Random Walk (RW) simulations based on a simplified representation of pore space in Mt were used to estimate and understand the effect of the arrangement of Mt particles on the meso- to macroscopic diffusivity of water. These theoretical calculations were complemented with quasielastic neutron scattering (QENS) measurements of aqueous diffusion in Mt with two pseudo-layers of water performed at four significantly different energy resolutions (i.e. observation times). The size of the interlayer and the size of Mt particles are two characteristic dimensions which determine the time dependent behavior of water diffusion in Mt. MD simulations show that at very short time scales water dynamics has the characteristic features of an oscillatory motion in the cage formed by neighbors in the first coordination shell. At longer time scales, the interaction of water with the surface determines the water dynamics, and the effect of confinement on the overall water mobility within the interlayer becomes evident. At time scales corresponding to an average water displacement equivalent to the average size of Mt particles, the effects of tortuosity are observed in the meso- to macroscopic pore scale simulations. Consistent with the picture obtained in the simulations, the QENS data can be described using a (local) 3D diffusion at short observation times, whereas at sufficiently long observation times a 2D diffusive motion is clearly observed. The effects of tortuosity measured in macroscopic tracer diffusion experiments are in qualitative agreement with RW simulations. By using experimental data to calibrate molecular and mesoscopic theoretical models, a consistent description of water mobility in clay minerals from the molecular to the macroscopic scale can be achieved. In turn, simulations help in choosing optimal conditions for the experimental measurements and the data interpretation. (C) 2014 Elsevier B.V. All rights reserved.

Environmental control on the occurrence of high-coercivity magnetic minerals and formation of iron sulfides in a 640 ka sediment sequence from Lake Ohrid (Balkans)

The bulk magnetic mineral record from Lake Ohrid, spanning the past 637 kyr, reflects large-scale shifts in hydrological conditions, and, superimposed, a strong signal of environmental conditions on glacial–interglacial and millennial timescales. A shift in the formation of early diagenetic ferrimagnetic iron sulfides to siderites is observed around 320 ka. This change is probably associated with variable availability of sulfide in the pore water. We propose that sulfate concentrations were significantly higher before  ∼  320 ka, due to either a higher sulfate flux or lower dilution of lake sulfate due to a smaller water volume. Diagenetic iron minerals appear more abundant during glacials, which are generally characterized by higher Fe / Ca ratios in the sediments. While in the lower part of the core the ferrimagnetic sulfide signal overprints the primary detrital magnetic signal, the upper part of the core is dominated by variable proportions of high- to low-coercivity iron oxides. Glacial sediments are characterized by high concentration of high-coercivity magnetic minerals (hematite, goethite), which relate to enhanced erosion of soils that had formed during preceding interglacials. Superimposed on the glacial–interglacial behavior are millennial-scale oscillations in the magnetic mineral composition that parallel variations in summer insolation. Like the processes on glacial–interglacial timescales, low summer insolation and a retreat in vegetation resulted in enhanced erosion of soil material. Our study highlights that rock-magnetic studies, in concert with geochemical and sedimentological investigations, provide a multi-level contribution to environmental reconstructions, since the magnetic properties can mirror both environmental conditions on land and intra-lake processes.

(Table T1) Geochemical composition of tochilinite and related minerals of ODP Hole 173-1068A

Tochilinite (approximately FeS(Mg,Fe)(OH)2) is locally abundant in Hole 1068A serpentinites from Cores 173-1068A-21R and 22R. It occurs in veins, as fillings in void space, and in intergrowths with serpentine and andradite. An apparently related mineral, but with Ca and Al largely replacing Mg, occurs in association with, and possibly as a replacement of, pyrrhotite in serpentinite breccias from the bottom of Core 173-1068A-20R. The transition from Mg-Fe-rich brucite tochilinites to Ca- and S-rich carbonate tochilinites is consistent with increasing sulfur and oxygen activity upsection. Tochilinite has been reported at other sites on the Iberia Abyssal Plain and is abundant to the point of being a rock-forming mineral in several samples from Site 1068. Rather than being a mineralogical curiosity, tochilinite appears to be common and a major sink for sulfur in the upper serpentinites of the Iberia Abyssal Plain.

Zircons, isotope ratios and element analyses from tonalite and oxide gabbro of the mid-ocean ridge from ODP Hole 176-735B

A limiting factor in the accuracy and precision of U/Pb zircon dates is accurate correction for initial disequilibrium in the 238U and 235U decay chains. The longest-lived-and therefore most abundant-intermediate daughter product in the 235U isotopic decay chain is 231Pa (T1/2 = 32.71 ka), and the partitioning behavior of Pa in zircon is not well constrained. Here we report high-precision thermal ionization mass spectrometry (TIMS) U-Pb zircon data from two samples from Ocean Drilling Program (ODP) Hole 735B, which show evidence for incorporation of excess 231Pa during zircon crystallization. The most precise analyses from the two samples have consistent Th-corrected 206Pb/238U dates with weighted means of 11.9325 ± 0.0039 Ma (n = 9) and 11.920 ± 0.011 Ma (n = 4), but distinctly older 207Pb/235U dates that vary from 12.330 ± 0.048 Ma to 12.140 ± 0.044 Ma and 12.03 ± 0.24 to 12.40 ± 0.27 Ma, respectively. If the excess 207Pb is due to variable initial excess 231Pa, calculated initial (231Pa)/(235U) activity ratios for the two samples range from 5.6 ± 1.0 to 9.6 ± 1.1 and 3.5 ± 5.2 to 11.4 ± 5.8. The data from the more precisely dated sample yields estimated DPazircon/DUzircon from 2.2-3.8 and 5.6-9.6, assuming (231Pa)/(235U) of the melt equal to the global average of recently erupted mid-ocean ridge basaltic glasses or secular equilibrium, respectively. High precision ID-TIMS analyses from nine additional samples from Hole 735B and nearby Hole 1105A suggest similar partitioning. The lower range of DPazircon/DUzircon is consistent with ion microprobe measurements of 231Pa in zircons from Holocene and Pleistocene rhyolitic eruptions (Schmitt (2007; doi:10.2138/am.2007.2449) and Schmitt (2011; doi:10.1146/annurev-earth-040610-133330)). The data suggest that 231Pa is preferentially incorporated during zircon crystallization over a range of magmatic compositions, and excess initial 231Pa may be more common in zircons than acknowledged. The degree of initial disequilibrium in the 235U decay chain suggested by the data from this study, and other recent high precision datasets, leads to resolvable discordance in high precision dates of Cenozoic to Mesozoic zircons. Minor discordance in zircons of this age may therefore reflect initial excess 231Pa and does not require either inheritance or Pb loss.

Geochemical analyses of secondary minerals from ODP Leg 115 basalts (Tables 1+2)

Basalts recovered along the Reunion hotspot track on Ocean Drilling Program (ODP) Leg 115 range in age from 34 Ma at Site 706 to 64 Ma at Site 707. They have undergone various degrees of secondary alteration. Within single holes the amount of alteration can vary from a few percent to near complete replacement of phenocrysts and groundmass by secondary minerals. Olivine appears to be the most susceptible to alteration and in some sections it is the only mineral altered. In other sections, olivine, pyroxene and plagioclase phenocrysts, and groundmass have been completely replaced by secondary minerals. Clays are the predominant form of secondary mineralization. In addition to replacing olivine, pyroxene, glass, and groundmass, clays have filled veins, vesicles, and voids. Minor amounts of calcite, zeolites, and K-feldspar were also detected. The clays that filled vesicles and veins often show color zonations of dark, opaque bands near the edges that grade into tan or green transparent regions in the centers of the veins. The electron microprobe was used to obtain chemical analyses of these veins as well as to characterize isolated clays that replaced specific minerals and filled voids and vesicles.

Sediment geochemistry of ODP Leg 127 sites

The relative effects of paleoceanographic and paleogeographic variations, sediment lithology, and diagenetic processes on the final preserved chemistry of Japan Sea sediments are evaluated by investigating the rare earth element (REE), major element, and trace element concentrations in 59 squeeze-cake whole-round and 27 physical-property sample residues from Sites 794, 795, and 797, cored during ODP Leg 127. The most important variation in sedimentary chemical composition is the increase in SiO2 concentration through the Pliocene diatomaceous sequences, which dilutes most other major and trace element components by various degrees. This biogenic input is largest at Site 794 (Yamato Basin), moderately developed at Site 797 (Yamato Basin), and of only minor importance at Site 795 (Japan Basin), potentially reflecting basinal contrasts in productivity with the Yamato Basin recording greater biogenic input than the Japan Basin and with the easternmost sequence of Site 794 lying beneath the most productive waters. There are few systematic changes in solid-phase chemistry resulting from the opal-A/opal-CT or opal-CT/quartz silica phase transformations. Most major and trace element concentrations are controlled by the aluminosilicate fraction of the sediment, although the effects of diagenetic silica phases and manganese carbonates are of localized importance. REE total abundances (Sum REE) in the Japan Sea are strongly dependent upon the paleoceanographic position of a given site with respect to terrigenous and biogenic sources. REE concentrations at Site 794 overall correspond well to aluminosilicate chemical indices and are strongly diluted by SiO2 within the upper Miocene-Pliocene diatomaceous sequence. Eu/Eu* values at Site 794 reach a maximum through the diatomaceous interval as well, most likely suggesting an association of Eu/Eu* with the siliceous component, or reflecting slight incorporation of a detrital feldspar phase. SumREE at Site 795 also is affiliated strongly with aluminosilicate phases and yet is diluted only slightly by siliceous input. At Site 797, SumREE is not as clearly associated with the aluminosilicate fraction, is correlated moderately to siliceous input, and may be sporadically influenced by detrital heavy minerals originating from the nearby rifted continental fragment composing the Yamato Rise. Ce/Ce* profiles at all three sites increase essentially monotonically with depth and record progressive diagenetic LREE fractionation. The observed Ce/Ce* increases are not responding to changes in the paleoceanographic oxygenation state of the overlying water, as there is no independent evidence to suggest the proper oceanographic conditions. Ce/Ce* correlates slightly better with depth than with age at the two Yamato Basin sites. The downhole increase in Ce/Ce* at Sites 794 and 797 is a passive response to the diagenetic transfer of LREE (except Ce) from sediment to interstitial water. At Site 795, the overall lack of correlation between Ce/Ce* and Lan/Ybn suggests that other processes mask the diagenetic behavior of all LREEs. First-order calculations of the Ce budget in Japan Sea waters and sediment indicate that ~20% of the excess Ce adsorbed by settling particles is recycled within the water column and that an additional ~38% is recycled at or near the seafloor. Thus, because the remaining excess Ce is only ~10% of the total Ce, there is not a large source of Ce to the deeply buried sediment, further suggesting that the downhole increase in Ce/Ce* is a passive response to diagenetic behavior of the other LREEs. The REE chemistry of Japan Sea sediment therefore predicts successive downhole addition of LREEs to deeply buried interstitial waters.