Field-scale water transport in unsaturated crystalline rock

Safe disposal of toxic wastes in geologic formations requires minimal water and gas movement in the vicinity of storage areas, Ventilation of repository tunnels or caverns built in solid rock can desaturate the near field up to a distance of meters from the rock surface, even when the surrounding geological formation is saturated and under hydrostatic pressures. A tunnel segment at the Grimsel test site located in the Aare granite of the Bernese Alps (central Switzerland) has been subjected to a resaturation and, subsequently, to a controlled desaturation, Using thermocouple psychrometers (TP) and time domain reflectometry (TDR), the water potentials psi and water contents theta were measured within the unsaturated granodiorite matrix near the tunnel wall at depths between 0 and 160 cm. During the resaturation the water potentials in the first 30 cm from the rock surface changed within weeks from values of less than -1.5 MPa to near saturation. They returned to the negative initial values during desaturation, The dynamics of this saturation-desaturation regime could be monitored very sensitively using the thermocouple psychrometers, The TDR measurements indicated that water contents changed dose to the surface, but at deeper installation depths the observed changes were within the experimental noise. The field-measured data of the desaturation cycle were used to test the predictive capabilities of the hydraulic parameter functions that were derived from the water retention characteristics psi(theta) determined in the laboratory. A depth-invariant saturated hydraulic conductivity k(s) = 3.0 x 10(-11) m s(-1) was estimated from the psi(t) data at all measurement depths, using the one-dimensional, unsaturated water flow and transport model HYDRUS Vogel er al., 1996, For individual measurement depths, the estimated k(s) varied between 9.8 x 10(-12) and 6.1 x 10(-11) The fitted k(s) values fell within the range of previously estimated k(s) for this location and led to a satisfactory description of the data, even though the model did not include transport of water vapor.

Resolving Cl and SO4 Profiles in a Clay-Rich Rock Sequence

The chloride and sulfate concentration profiles in a 260 m thick clay-rich Mesozoic sediment sequence have been analyzed by various methods. Chloride data generally indicate a good consistency between different methods if anion exclusion is accounted for in leaching tests. For sulfate, however, there is an apparent inconsistency between leaching data and those obtained from the other methods, which points to the dissolution of a sulfur-bearing mineral. Traces of diagenetic gypsum seem to be a likely source, but other sulfur minerals cannot be ruled out.

36Cl production rate from K-spallation in the European Alps (Chironico landslide, Switzerland)

The abundant production of in situ cosmogenic 36Cl from potassium renders 36Cl measurements in K-rich rocks or minerals, such as K-feldspars, potentially useful for precisely dating rock surfaces, either in single-nuclide or in multi-nuclide studies, for example combined with 10Be measurements in quartz. However, significant discrepancies in experimentally calibrated 36Cl production rates from spallation of potassium (36PK-sp), referenced to sea-level/high-latitude (SLHL), limit the accuracy of 36Cl dating from K-rich lithologies. We present a new 36Cl calibration using K-feldspars, in which K-spallation is the most dominant 36Cl production pathway (>92% of total 36Cl), thus minimizing uncertainties from the complex multi-pathway 36Cl production systematics. The samples are derived from boulders of an ∼13.4 ka-old landslide in the Swiss Alps (∼820 m, 46.43°N, 8.85°E). We obtain a local 36PK-sp of 306 ± 16 atoms 36Cl (g K)−1 a−1 and an SLHL 36PK-sp of 145.5 ± 7.7 atoms 36Cl (g K)−1 a−1, when scaled with a standard scaling protocol (‘Lm’). Applying this SLHL 36PK-sp to determine 36Cl exposure ages of K-feldspars from 10Be-dated moraine boulders yields excellent agreement, confirming the validity of the new SLHL 36PK-sp for surface exposure studies, involving 36Cl in K-feldspars, in the Alps.

Quaternary uplift rates of the Central Anatolian Plateau, Turkey: insights from cosmogenic isochron-burial nuclide dating of the Kızılırmak River terraces

The Central Anatolian Plateau (CAP) in Turkey is a relatively small plateau (300 × 400 km) with moderate average elevations of ∼1 km situated between the Pontide and Tauride orogenic mountain belts. Kızılırmak, which is the longest river (1355 km) within the borders of Turkey, flows within the CAP and slowly incises into lacustrine and volcaniclastic units before finally reaching the Black Sea. We dated the Cappadocia section of the Kızılırmak terraces in the CAP by using cosmogenic burial and isochron-burial dating methods with 10Be and 26Al as their absolute dating can provide insight into long-term incision rates, uplift and climatic changes. Terraces at 13, 20, 75 and 100 m above the current river indicate an average incision rate of 0.051 ± 0.01 mm/yr (51 ± 1 m/Ma) since ∼1.9 Ma. Using the base of a basalt fill above the modern course of the Kızılırmak, we also calculated 0.05–0.06 mm/yr mean incision and hence rock uplift rate for the last 2 Ma. Although this rate might be underestimated due to normal faulting along the valley sides, it perfectly matches our results obtained from the Kızılırmak terraces. Although up to 5–10 times slower, the Quaternary uplift of the CAP is closely related to the uplift of the northern and southern plateau margins respectively.

Benchmark problems for reactive transport modeling of the generation and attenuation of acid rock drainage

Acid rock drainage (ARD) is a problem of international relevance with substantial environmental and economic implications. Reactive transport modeling has proven a powerful tool for the process-based assessment of metal release and attenuation at ARD sites. Although a variety of models has been used to investigate ARD, a systematic model intercomparison has not been conducted to date. This contribution presents such a model intercomparison involving three synthetic benchmark problems designed to evaluate model results for the most relevant processes at ARD sites. The first benchmark (ARD-B1) focuses on the oxidation of sulfide minerals in an unsaturated tailing impoundment, affected by the ingress of atmospheric oxygen. ARD-B2 extends the first problem to include pH buffering by primary mineral dissolution and secondary mineral precipitation. The third problem (ARD-B3) in addition considers the kinetic and pH-dependent dissolution of silicate minerals under low pH conditions. The set of benchmarks was solved by four reactive transport codes, namely CrunchFlow, Flotran, HP1, and MIN3P. The results comparison focused on spatial profiles of dissolved concentrations, pH and pE, pore gas composition, and mineral assemblages. In addition, results of transient profiles for selected elements and cumulative mass loadings were considered in the intercomparison. Despite substantial differences in model formulations, very good agreement was obtained between the various codes. Residual deviations between the results are analyzed and discussed in terms of their implications for capturing system evolution and long-term mass loading predictions.

(Table 1) K-Ar age of alkaline igneous rocks from the Jan Mayen fracture zone

Bottom morphology of the Jan Mayen transform fracture zone and rock chemistry data show that petrological and chemical specific features of igneous rocks can result from higher permeability of the transform fracture zone and deeper penetration of ocean water into the lithosphere in comparison with rift zones of the Kolbeinsey and Mohn's mid-ocean ridges. Age of alkaline magmatism of the Jan Mayen fracture zone is similar to that of rift zones due to palingenesis of metamorphosed and hydrated mantle and crustal rocks.

Sedimentological, physical, bulk geochemical, elemental and rock magnetic properties of the Middle to Late Pleistocene Rodderberg (Eifel/Germany)

The sediment record from Rodderberg potentially provides a climate and environmental record spanning at least the last ca 130 ka. Results from a low resolution pilot study reveal characteristic fluctuations that can be related to global climate variability as reflected in marine isotope stages and document the potential of this site for continuous and high-resolution investigations of the Middle to Late Pleistocene. Here we document the tentative lithology drilled, and show how the elemental composition can be interpreted with regard to lake level fluctuations, related redox conditions, but also to grain-size distribution and changes in lacustrine productivity. Finally, based on major lithological changes, a preliminary depth/age model is suggested that allows reassessing published luminescence ages from the same site.

Appendix A. Whole-rock major, trace and rare earth element analyses of the Ulubey volcanics, Appendix B. Analytical standards and detection limits for analyses

Collisional and post-collisional volcanic rocks in the Ulubey (Ordu) area at the western edge of the Eastern Pontide Tertiary Volcanic Province (EPTVP) in NE Turkey are divided into four suites; Middle Eocene (49.4-44.6 Ma) aged Andesite-Trachyandesite (AT), Trachyandesite-Trachydacite-Rhyolite (TTR), Trachydacite-Dacite (TD) suites, and Middle Miocene (15.1 Ma) aged Trachybasalt (TB) suite. Local stratigraphy in the Ulubey area starts with shallow marine environment sediments of the Paleocene-Eocene time and then continues extensively with sub-aerial andesitic to rhyolitic and rare basaltic volcanism during Eocene and Miocene time, respectively. Petrographically, the volcanic rocks are composed primarily of andesites/trachyandesites, with minor trachydacites/rhyolites, basalts/trachybasalts and pyroclastics, and show porphyric, hyalo-microlitic porphyric and rarely glomeroporphyric, intersertal, intergranular, fluidal and sieve textures. The Ulubey (Ordu) volcanic rocks indicate magma evolution from tholeiitic-alkaline to calc-alkaline with medium-K contents. Primitive mantle normalized trace element and chondrite normalized rare earth element (REE) patterns show that the volcanic rocks have moderate light rare earth element (LREE)/heavy rare earth element (HREE) ratios relative to E-Type MORB and depletion in Nb, Ta and Ti. High Th/Yb ratios indicate parental magma(s) derived from an enriched source formed by mixing of slab and asthenospheric melts previously modified by fluids and sediments from a subduction zone. All of the volcanic rocks share similar incompatible element ratios (e.g., La/Sm, Zr/Nb, La/Nb) and chondrite-normalized REE patterns, indicating that the basic to acidic rocks originated from the same source. The volcanic rocks were produced by the slab dehydration-induced melting of an existing metasomatized mantle source, and the fluids from the slab dehydration introduced significant large ion lithophile element (LILE) and LREE to the source, masking its inherent HFSE-enriched characteristics. The initial 87Sr/86Sr (0.7044-0.7050) and eNd (-0.3 to +3.4) ratios of the volcanics suggest that they originated from an enriched lithospheric mantle source with low Sm/Nd ratios. Integration of the geochemical, petrological and isotopical with regional and local geological data suggest that the Tertiary volcanic rocks from the Ulubey (Ordu) area were derived from an enriched mantle, which had been previously metasomatized by fluids derived from subducted slab during Eocene to Miocene in collisional and post-collisional extension-related geodynamic setting following Late Mesozoic continental collision between the Eurasian plate and the Tauride-Anatolide platform.

Determination of sulfite and manganese in sediments from the K/T boundary

The solid-state-physics technique of electron spin resonance (ESR) has been employed in an exploratory study of marine limestones and impact-related deposits from Cretaceous-Tertiary (KT) boundary sites including Spain (Sopelana and Caravaca), New Jersey (Bass River), the U.S. Atlantic continental margin (Blake Nose, ODP Leg 171B/1049/A), and several locations in Belize and southern Mexico within -600 km of the Chicxulub crater. The ESR spectra of SO3(1-) (a radiation-induced point defect involving a sulfite ion substitutional for CO3(2-) which has trapped a positive charge) and Mn(2+) in calcite were singled out for analysis because they are unambiguously interpretable and relatively easy to record. ESR signal strengths of calcite-related SO3(1-) and Mn(2+) have been studied as functions of stratigraphic position in whole-rock samples across the KT boundary at Sopelana, Caravaca, and Blake Nose. At all three of these sites, anomalies in SO3(1-) and/or Mn(2+) intensities are noted at the KT boundary relative to the corresponding background levels in the rocks above and below. At Caravaca, the SO3(1-) background itself is found to be lower by a factor of 2.7 in the first 30,000 years of the Tertiary relative to its steady-state value in the last 15,000 years of the Cretaceous, indicating either an abrupt and quasi-permanent change in ocean chemistry (or temperature) or extinction of the marine biota primarily responsible for fixing sulfite in the late Cretaceous limestones. An exponential decrease in the Mn(2+) concentration per unit mass calcite, [Mn(2+)], as the KT boundary at Caravaca is approached from below (1/e characteristic length =1.4 cm) is interpreted as a result of post-impact leaching of the seafloor. Absolute ESR quantitative analyses of proximal impact deposits from Belize and southern Mexico group naturally into three distinct fields in a twodimensional [SO3(1-)]-versus-[Mn(2+)] scatter plot. These fields contain (I) limestone ejecta clasts, (II) accretionary lapilli, and (III) a variety of SO3(1-) -depleted/Mn(2+) enriched impact deposits. Data for the investigated non-impact-related Cretaceous and Tertiary marine limestones (Spain and Blake Nose) fall outside of these three fields. With reference to thes enon-impact deposits, fields I, II, and III can be respectively characterized as Mn(2+) -depleted, SO3(1-) -enhanced, and SO3(1-) -depleted. It is proposed that (1) field I represents calcites from the Yucatin Platform, and that the Mn(2+) -depleted signature can be used as an indicator of primary Chicxulub ejecta in deep marine environments and (2) field II represents calcites that include a component formed in the vapor plume, either from condensation in the presence of CO2/SO3(1-) -rich vapors, or reactions between CaO and CO2/SO3 rich vapors, and that this SO3(1-) -enhanced signature can be used as an indicator of impact vapor plume deposits. Given these two propositions, the ESR data for the Blake Nose deposits are ascribed to the presence of basal coarse calcitic Chicxulub ejecta clasts, while the finer components that are increasingly represented toward the top are interpreted to contain high- SO3(1-) calcite from the vapor plume. The apparently-undisturbed Bass River deposit may contain even higher concentrations of vapor-plume calcite. None of the three components included in field III appear to be represented at distal, deep marine KT-boundary sites; this field may include several types of impact-related deposits of diverse origins and diagenetic histories.