Characterization of granite matrix porosity and pore-space geometry by in situ and laboratory methods

Most available studies of interconnected matrix porosity of crystalline rocks are based on laboratory investigations; that is, work on samples that have undergone stress relaxation and were affected by drilling and sample preparation. The extrapolation of the results to in situ conditions is therefore associated with considerable uncertainty, and this was the motivation to conduct the ‘in situ Connected Porosity’ experiment at the Grimsel Test Site (Central Swiss Alps). An acrylic resin doped with fluorescent agents was used to impregnate the microporous granitic matrix in situ around an injection borehole, and samples were obtained by overcoring. The 3-D structure of the porespace, represented by microcracks, was studied by U-stage fluorescence microscopy. Petrophysical methods, including the determination of porosity, permeability and P -wave velocity, were also applied. Investigations were conducted both on samples that were impregnated in situ and on non-impregnated samples, so that natural features could be distinguished from artefacts. The investigated deformed granites display complex microcrack populations representing a polyphase deformation at varying conditions. The crack population is dominated by open cleavage cracks in mica and grain boundary cracks. The porosity of non-impregnated samples lies slightly above 1 per cent, which is 2–2.5 times higher than the in situ porosity obtained for impregnated samples. Measurements of seismic velocities (Vp ) on spherical rock samples as a function of confining pressure, spatial direction and water saturation for both non-impregnated and impregnated samples provide further constraints on the distinction between natural and induced crack types. The main conclusions are that (1) an interconnected network of microcracks exists in the whole granitic matrix, irrespective of the distance to ductile and brittle shear zones, and (2) conventional laboratory methods overestimate the matrix porosity. Calculations of contaminant transport through fractured media often rely on matrix diffusion as a retardation mechanism.

Where to find 1.5 million yr old ice for the IPICS "Oldest-Ice" ice core

The recovery of a 1.5 million yr long ice core from Antarctica represents a keystone of our understanding of Quaternary climate, the progression of glaciation over this time period and the role of greenhouse gas cycles in this progression. Here we tackle the question of where such ice may still be found in the Antarctic ice sheet. We can show that such old ice is most likely to exist in the plateau area of the East Antarctic ice sheet (EAIS) without stratigraphic disturbance and should be able to be recovered after careful pre-site selection studies. Based on a simple ice and heat flow model and glaciological observations, we conclude that positions in the vicinity of major domes and saddle position on the East Antarctic Plateau will most likely have such old ice in store and represent the best study areas for dedicated reconnaissance studies in the near future. In contrast to previous ice core drill site selections, however, we strongly suggest significantly reduced ice thickness to avoid bottom melting. For example for the geothermal heat flux and accumulation conditions at Dome C, an ice thickness lower than but close to about 2500 m would be required to find 1.5 Myr old ice (i.e., more than 700 m less than at the current EPICA Dome C drill site). Within this constraint, the resolution of an Oldest-Ice record and the distance of such old ice to the bedrock should be maximized to avoid ice flow disturbances, for example, by finding locations with minimum geothermal heat flux. As the geothermal heat flux is largely unknown for the EAIS, this parameter has to be carefully determined beforehand. In addition, detailed bedrock topography and ice flow history has to be reconstructed for candidates of an Oldest-Ice ice coring site. Finally, we argue strongly for rapid access drilling before any full, deep ice coring activity commences to bring datable samples to the surface and to allow an age check of the oldest ice.

Dating groundwater in the Bohemian Cretaceous Basin: Understanding tracer variations in the subsurface

The northern section of the Bohemian Cretaceous Basin has been the site of intensive U exploitation with harmful impacts on groundwater quality. The understanding of groundwater flow and age distribution is crucial for the prediction of the future dispersion and impact of the contamination. State of the art tracer methods (3H, 3He, 4He, 85Kr, 39Ar and 14C) were, therefore, used to obtain insights to ageing and mixing processes of groundwater along a north–south flow line in the centre of the two most important aquifers of Cenomanian and middle Turonian age. Dating of groundwater is particularly complex in this area as: (i) groundwater in the Cenomanian aquifer is locally affected by fluxes of geogenic and biogenic gases (e.g. CO2, CH4, He) and by fossil brines in basement rocks rich in Cl and SO4; (ii) a thick unsaturated zone overlays the Turonian aquifer; (iii) a periglacial climate and permafrost conditions prevailed during the Last Glacial Maximum (LGM), and iv) the wells are mostly screened over large depth intervals. Large disagreements in 85Kr and 3H/3He ages indicate that processes other than ageing have affected the tracer data in the Turonian aquifer. Mixing with older waters (>50 a) was confirmed by 39Ar activities. An inverse modelling approach, which included time lags for tracer transport throughout the unsaturated zone and degassing of 3He, was used to estimate the age of groundwater. Best fits between model and field results were obtained for mean residence times varying from modern up to a few hundred years. The presence of modern water in this aquifer is correlated with the occurrence of elevated pollution (e.g. nitrates). An increase of reactive geochemical indicators (e.g. Na) and radiogenic 4He, and a decrease in 14C along the flow direction confirmed groundwater ageing in the deeper confined Cenomanian aquifer. Radiocarbon ages varied from a few hundred years to more than 20 ka. Initial 14C activity for radiocarbon dating was calibrated by means of 39Ar measurements. The 14C age of a sample recharged during the LGM was further confirmed by depleted stable isotope signatures and near freezing point noble gas temperature. Radiogenic 4He accumulated in groundwater with concentrations increasing linearly with 14C ages. This enabled the use of 4He to validate the dating range of 14C and extend it to other parts of this aquifer. In the proximity of faults, 39Ar in excess of modern concentrations and 14C dead CO2 sources, elevated 3He/4He ratios and volcanic activity in Oligocene to Quaternary demonstrate the influence of gas of deeper origin and impeded the application of 4He, 39Ar and 14C for groundwater dating.

Origin and evolution of reactive and noble gases dissolved in matrix pore water

Reactive and noble gases dissolved in matrix pore water of low permeable crystalline bedrock were successfully extracted and characterized for the fist time based on drillcore samples from the Olkiluoto investigation site (SW Finland). Interaction between matrix pore water and fracture groundwater occurs predominately by diffusion. Changes in the chemical and isotopic composition of gases dissolved in fracture groundwater are transmitted and preserved in the pore water. Absolute concentrations, their ratios and the stable carbon isotope signature of hydrocarbon gases dissolved in pore water give valuable indications about the evolution of these gases in the nearby-flowing fracture groundwaters. Inert noble gases dissolved in matrix pore water and their isotopes combined with their in-situ production and accumulation rates deliver information about the residence time of pore water.

The origin of oriented lakes: evidence from the Bolivian Amazon

The presence of hundreds of rectangular and oriented lakes is one of the most striking characteristics of the Llanos de Moxos (LM) landscape in the Bolivian Amazon. Oriented lakes also occur in the Arctic coastal plains of Russia, Alaska and Canada and along the Atlantic Coastal Plain from northeast Florida to southeast New Jersey and along the coast of northeast Brazil. Many different mechanisms have been proposed for their formation. In the LM, Plafker's (1964) tectonic model, in which subsidence results from the propagation of bedrock faults through the foreland sediments, is the most accepted. However, this model has not been verified. Here, we present new results from stratigraphic transects across the borders of three rectangular and oriented lakes in the LM. A paleosol buried under mid-Holocene sediments is used as a stratigraphic marker to assess the vertical displacement of sediments on both sides of the alleged faults. Our results show that there is no vertical displacement and, therefore, that Plafker's model can be ruled out. We suggest that, among all the proposed mechanisms behind lake formation, the combined action of wind and waves is the most likely. The evidence from the LM provides new hints for the formation of oriented lakes worldwide.

Assessing the Environmental Hazard of Using Seawater for Ore Processing at the Lasail Mine Site in the Sultanate of Oman

The Lasail mining area (Sultanate of Oman) was contaminated by acid mine drainage during the exploitation and processing of local and imported copper ore and the subsequent deposition of sulphide-bearing waste material into an unsealed tailings dump. In this arid environment, the use of seawater in the initial stages of ore processing caused saline contamination of the fresh groundwater downstream of the tailings dump. After detection of the contamination in the 1980s, different source-controlled remediation activities were conducted including a seepage water collection system and, in 2005, surface sealing of the tailings dump using an HDPE-liner to prevent further infiltration of meteoric water. We have been assessing the benefits of the remediation actions undertaken so far. We present chemical and isotopic (δ18O, δ 2H, 3H) groundwater data from a long-term survey (8–16 years) of the Wadi Suq aquifer along a 28 km profile from the tailings dump to the Gulf of Oman. Over this period, most metal concentrations in the Wadi Suq groundwater decreased below detection limits. In addition, in the first boreholes downstream of the tailings pond, the salinity contamination has decreased by 30 % since 2005. This decrease appears to be related to the surface coverage of the tailings pond, which reduces flushing of the tailings by the sporadic, but commonly heavy, precipitation events. Despite generally low metal concentrations and the decreased salinity, groundwater quality still does not meet the WHO drinking water guidelines in more than 90 % of the Wadi Suq aquifer area. The observations show that under arid conditions, use of seawater for ore processing or any other industrial activity has the potential to contaminate aquifers for decades.

Recharge Areas and Geochemical Evolution of Groundwater in an Alluvial Aquifer System in the Sultanate of Oman

A regional hydrogeochemical model was developed to evaluate the geochemical evolution of different groundwaters in an alluvial aquifer system in the Interior of Oman. In combination with environmental isotopes the model is able to extract qualitative and quantitative information about recharge, groundwater flow paths and hydraulic connections between different aquifers. The main source of water to the alluvial aquifer along the flow paths ofWadi Abyadh andWadi M’uaydin in the piedmont is groundwater from the high-altitude areas of the Jabal Akhdar and local infiltration along the wadi channels. In contrast, the piedmont alluvial aquifer alongWadi Halfayn is primarily replenished by lateral recharge from the ophiolite foothills to the east besides smaller contributions from the Jabal Akhdar and local infiltration. Further down gradient in the Southern Alluvial Plain aquifer a significant source of recharge is direct infiltration of rain and surface runoff, originating from a moisture source that approaches Oman from the south. The model shows that the main geochemical evolution of the alluvial groundwaters occurs along the flow path from the piedmont to the Southern Alluvial Plain, where dedolomitization is responsible for the observed changes in the chemical and carbon isotope composition in these waters.

Helium transfer from water into quartz crystals: A new approach for porewater dating

Several important fundamental and applied problems require a quantification of slow rates of groundwater flow. To resolve these problems helium appears to be a promising tracer. In this contribution we discuss a new approach, which gives the helium inventory in a rock – pore water system by using the relevant mineral record, i.e., without extraction and investigation of the porewater samples. Some U- and Th-poor minerals such as quartz (quartz separates from Permo-Carboniferous Formation, sandstone–shale interlayering, Molasses Basin, Northern Switzerland, hereafter PCF, are used in this study) contain excessive helium having migrated into their internal helium-accessible volume (HAV) from the surrounding porewater [I.N. Tolstikhin, B.E. Lehmann, H.H. Loosli, A. Gautschi, Helium and argon isotopes in rocks, minerals and related groundwaters: a case study in Northern Switzerland, Geochim. Cosmochim. Acta 60 (1996) 1497–1514]. These volumes are estimated by using helium as a nano-size penetrating tool, i.e., by saturation of the minerals with helium under controlled pressure–temperature conditions and subsequent measurements of the helium-saturated concentrations. In the quartz separates HAV/total volume ratios vary from 0.017% to 0.16%; along with the measured initial (unsaturated) He concentration the HAV gives the internal helium pressure, the mean value obtained for 7 samples (25 sample aliquots) is P=0.45F0.15 atm (1 r). The product of helium pressure and solubility (7.35_10_3 cc STP He/cc H2O for the temperature and salinity of PCF aquifers reported in [F.J. Pearson, W. Balderer, H.H. Loosli, B.E. Lehmann, A. Matter, T. Peters, H. Schmassmann, A. Gautschi, Applied Isotope Hydrogeology–A Case Study in Northern Switzerland, Elsevier Amsterdam, 1991, 439 pp.]) is the mineral-derived He concentration in the respective porewater, CPW=0.0035F0.0017 cc He/cc H2O. This value is in full accord with measured He concentrations in PCF aquifers, CPCF, varying from 0.0045 to 0.0016 cc He/cc H2O. This agreement validates the proposed approach and also shows that the mineral–porewater helium–concentration equilibrium has been established. Indeed, estimates of the He-migration rates through our quartz samples show that in ~6000 years the internal pressure should equilibrate with He-concentration in related porewater of PCF, and this time interval is short compared to independent estimates [I.N. Tolstikhin, B.E. Lehmann, H.H. Loosli, A. Gautschi, Helium and argon isotopes in rocks, minerals and related groundwaters: a case study in Northern Switzerland, Geochim. Cosmochim. Acta 60 (1996) 1497–1514]. The helium inventory in the rock–porewater assemblage shows that helium abundance in pore waters is indeed important. In shale samples (with ~15% porosity) porewaters contain more helium than the host minerals altogether. Porewater heliumconcentration profiles, available from the mineral record, along with helium production rates are input parameters allowing model(s) of helium migration through a hydrological structure to be developed. Quite high helium concentrations in PCF porewaters imply slow removal mechanisms, which will be discussed elsewhere.

210Pb dating of the Miaoergou ice core from the eastern Tien Shan, China

In 2005, two ice cores with lengths of 58.7 and 57.6 m respectively to bedrock were recovered from the Miaoergou flat-topped glacier (43 degrees 03 ' 19 '' N, 94 degrees 19 ' 21 '' E; 4512 m a.s.l.), eastern Tien Shan. Pb-210 dating of one of the ice cores (57.6 m) was performed, and an age of AD 1851 +/- 6 at a depth of 35.2 m w.e. was determined. For the period AD 1851-2005, a mean annual net accumulation of 229 +/- 7 mm w.e. a(-1) was calculated. At the nearby oasis city of Hami (similar to 80 km from the Miaoergou flat-topped glacier) the annual precipitation rate is 38 mm w.e. a(-1), hence glacial meltwater is a major water supply for local residents. The surface activity concentration of Pb-210(ex) was found to be similar to 400 mBq kg(-1), which is higher than observed at other continental sites such as Belukha, Russia, and Tsambagarav, Mongolia, which have surface activity concentrations of 280 mBq kg(-1). The Pb-210 dating agrees well with the chronological sequence deduced from the annual-layer counting resulting from the seasonalities of delta O-18 and trace metals for the period AD 1953-2005, and beta-activity horizons resulting from atmospheric nuclear testing during the period AD 1962-63. We conclude that Pb-210 analysis is a suitable method for obtaining a continuous dating of the Miaoergou ice core for similar to 160 years, which can also be applied to other ice cores recovered from the mountains of western China.

Continental degassing of 4He by surficial discharge of deep groundwater

Radiogenic He is produced by the decay of uranium and thorium in the Earth’s mantle and crust. From here, it is degassed to the atmosphere and eventually escapes to space. Assuming that all of the 4He produced is degassed, about 70% of the total He degassed from Earth comes from the continental crust. However, the outgoing flux of crustal He has not been directly measured at the Earth’s surface and the migration pathways are poorly understood. Here we present measurements of helium isotopes and the long-lived cosmogenic radio-isotope Kr in the deep, continental-scale Guarani aquifer in Brazil and show that crustal He reaches the atmosphere primarily by the surficial discharge of deep groundwater. We estimate that He in Guarani groundwater discharge accounts for about 20% of the assumed global flux from continental crust, and that other large aquifers may account for about 33%. Old groundwater ages suggest that He in the Guarani aquifer accumulates over half- to one-millionyear timescales. We conclude that He degassing from the continents is regulated by groundwater discharge, rather than episodic tectonic events, and suggest that the assumed steady state between crustal production and degassing of He, and its resulting atmospheric residence time, should be re-examined.

The geomorphology and morphometry of the banded terrain in Hellas basin, Mars

Hellas basin is a large impact basin situated in the southern highlands of Mars. The north-western part of the basin has the lowest elevation (-7.5 km) on the planet and contains a possibly unique terrain type, which we informally call “banded terrain”. The banded terrain is made up of smooth-looking banded deposits that display signs of viscous behavior and a paucity of superimposed impact craters. In this study, we use newly acquired high spatial resolution images from the High Resolution Imaging Science Experiment (HiRISE) in addition to existing datasets to characterize the geomorphology, the morphometry and the architecture of the banded terrain. The banded terrain is generally confined to the NW edge of the Alpheus Colles plateau. The individual bands are ~3–15 km-long, ~0.3 km-wide and are separated by narrow inter-band depressions, which are ~65 m-wide and ~10 m-deep. The bands display several morphologies that vary from linear to concentric forms. Morphometric analysis reveals that the slopes along a given linear or lobate band ranges from 0.5° to 15° (average~6°), whereas the concentric bands are located on flatter terrain (average slope~2–3°). Crater-size frequency analysis yields an Amazonian-Hesperian boundary crater retention age for the terrain (~3 Gyr), which together, with the presence of very few degraded craters, either implies a recent emplacement, resurfacing, or intense erosion. The apparent sensitivity to local topography and preference for concentrating in localized depressions is compatible with deformation as a viscous fluid. In addition, the bands display clear signs of degradation and slumping at their margins along with a suite of other features that include fractured mounds, polygonal cracks at variable size-scales, and knobby/hummocky textures. Together, these features suggest an ice-rich composition for at least the upper layers of the terrain, which is currently being heavily modified through loss of ice and intense weathering, possibly by wind.

Potential desiccation cracks on Mars: A synthesis from modeling, analogue-field studies, and global observations

Potential desiccation polygons (PDPs) are polygonal surface patterns that are a common feature in Noachian-to-Hesperian-aged phyllosilicate- and chloride-bearing terrains and have been observed with size scales that range from cm-wide (by current rovers) to 10s of meters-wide. The global distribution of PDPs shows that they share certain traits in terms of morphology and geologic setting that can aid identification and distinction from fracturing patterns caused by other processes. They are mostly associated with sedimentary deposits that display spectral evidence for the presence of Fe/Mg smectites, Al-rich smectites or less commonly kaolinites, carbonates, and sulfates. In addition, PDPs may indicate paleolacustrine environments, which are of high interest for planetary exploration, and their presence implies that the fractured units are rich in smectite minerals that may have been deposited in a standing body of water. A collective synthesis with new data, particularly from the HiRISE camera suggests that desiccation cracks may be more common on the surface of Mars than previously thought. A review of terrestrial research on desiccation processes with emphasis on the theoretical background, field studies, and modeling constraints is presented here as well and shown to be consistent with and relevant to certain polygonal patterns on Mars.

The role of climate-driven chemical weathering on soil production

Climate plays an important role in controlling rates of weathering and weathered regolith production. Regolith production functions, however, seldom take climate parameters into account. Based on a climate-dependent weathered regolith production model, at low denudation rates, relative regolith thicknesses are less sensitive to changes in precipitation rates, while at high denudation rates, small changes in climatic parameters can result in complete stripping of hillslopes. This pattern is compounded by the long residence times and system response times associated with low denudation rates, and vice versa. As others have shown, the transition between regolith-mantled and bedrock slopes is dependent on the ratio of denudation to production. Here, we further suggest that this is itself a function of precipitation rate and temperature. We suggest that climatic parameters can be easily incorporated into existing soil production models and that such additions improve the predictive power of soil production models. (C) 2013 Elsevier B.V. All rights reserved.

Drilling into debris-rich basal ice at the bottom of the NEEM (Greenland) borehole

After the NEEM (Greenland) deep ice-core drilling was declared terminated with respect to developing stratigraphic climate reconstructions, efforts were turned toward collecting basal ice-sheet debris and, if possible, drilling into the bedrock itself. In 2010, several meters of banded debris-rich ice were obtained under normal ice-drilling operations with the NEEM version of the Hans Tausen (HT) drill, but further penetration was obstructed by a rock in the path of the drill head at 2537.36 m. During short campaigns in 2011 and 2012, attempts were made to penetrate further using various reinforced ice cutters mounted on the HT drill head, tailored to cut through rock. These had some success in penetrating coarse material, but produced severely damaged cutters. Additionally a 51 mm diameter diamond cutting tipped rock drill was adapted to fit the NEEM drill. With this device, several additional meters of core containing subglacial sediments, rocks and rock fragments were collected. With these tools 1.39 m of additional material were obtained during the 2011 field season, and 7.1 m during 2012. Subglacial water refreezing into the newly formed borehole hindered further penetration, and the bedrock interface was not reached before final closure of the NEEM Camp.

Ce(III) and Ce(IV) (re)distribution and fractionation in a laterite profile from Madagascar: Insights from in situ XANES spectroscopy at the Ce LIII-edge

The distribution of trivalent and tetravalent cerium, Ce(III) and Ce(IV) respectively, in a lateritic profile from Madagascar, has been characterized by X-ray-absorption near-edge structure (XANES) spectroscopy at the Ce LIII-edge on the LUCIA beamline (SOLEIL synchrotron, France). XANES spectra were acquired on bulk-rock samples as well as on specific lateritic minerals or polymineral zones (in-situ measurements) of the tonalite bedrock and the three overlying weathered horizons (C-, B- and A-horizons). Geochemically, the bedrock, and the A- and C-horizons show similar rare earth element content (REE = 363–405 mg/kg). They also display the same positive Ce-anomaly (CeCN/Ce∗ = 1.12–1.45), which is therefore likely to be inherited from the bedrock. In the B-horizon, the higher REE content (REE = 2194 mg/kg) and the larger Ce-anomaly (CeCN/Ce∗ = 4.26) are consistent with an accumulation zone caused by the evaporation of groundwater during the dry season. There is a good agreement between the Ce(III)/Cetotal ratio (XCe(III)) deduced from the positive Ce-anomaly (bulk-rock geochemical data) and that derived from XANES spectroscopy on the same bulk-rock samples (BR-XCe(III)-XANES) in the bedrock, and the C- and B-horizons. In the A-horizon, XANES measurements on bulk rock and minerals revealed a higher BR-XCe(III)-XANES (up to 100%) compared to the XCe(III) deduced from geochemical data (XCe(III) = 79%). The preservation of a positive Ce-anomaly in the A-horizon suggests that the Ce mobilization and redistribution during weathering occurred with no significant Ce fractionation from other trivalent REE. Remarkably, the only investigated sample where cerianite is observed belongs to the B-horizon. Within this horizon, Ce oxidation state varies depending on the microstructural position (porosity, cracks, clay-rich groundmass). The highest Ce(IV) concentrations are measured in cerianite (and aluminophosphates) localized in pores at the vicinity of Mn-rich domains (XCe(III)-XANES = 30–51%). Therefore, Ce fractionation from other REE is attributed to a Ce oxidation and precipitation potentially assisted by oxyhydroxide scavenging. In the C-horizon, Ce(III) and Ce(IV) are mainly distributed in REE-minerals of the rhabdophane group found in pores and cracks. The similarity between the Ce(III) proportion of rhabdophane grains (XCe(III)-XANES = 74–89%) with that of the bedrock (BR-XCe(III)-XANES = 79%) suggests no significant fractionation of Ce(III) and Ce(IV) between solution and mineral during the successive stages of primary REE-mineral alteration, transport in solution and secondary precipitation in the incipient stages of weathering. Overall, our novel spectroscopic approach shows that Ce is not necessarily oxidized nor fractionated from other REE during weathering in lateritic conditions. This implies that like Ce(III), Ce(IV) can be mobilized in aqueous fluids during weathering, possibly thanks to complexation with organic molecules, and can precipitate together with Ce(III) in secondary REE-bearing minerals. The corollary is that (paleo)redox reconstructions in soils and/or sediments based on Ce-anomaly in weathered rocks or minerals must be interpreted with caution.