Suction-controlled experiments on Boom clay

The understanding of the thermo-hydro-mechanical behaviour of a clay barrier is needed for the prediction of its final in situ properties after the hydration and thermal transient in a radioactive waste repository.

As part of the CEC 1990–1994 R&D programme on radioactive waste management and storage, the CEA (Fr), CIEMAT (Sp), ENRESA (Sp), SCK · CEN (B), UPC (Sp) and UWCC (UK) have carried out a joint project on unsaturated clay behaviour (Volckaert et al., 1996). The aim of the study is to analyse and model the behaviour of a clay-based engineered barrier during its hydration phase under real repository conditions. The hydro-mechanical and thermo-hydraulic models developed in this project have been coupled to describe stress/strain behaviour, moisture migration and heat transfer. A thermo-hydraulic model has also been coupled to a geochemical code to describe the migration and formation of chemical species.

In this project, suction-controlled experiments have been performed on Boom clay (B), FoCa clay (Fr) and Almeria bentonite (Sp). The aim of these experiments is to test the validity of the interpretive model developed by Alonso and Gens (Alonso et al., 1990), and to build a database of unsaturated clay thermo-hydro-mechanical parameters. Such a database can then be used for validation exercises in which in situ experiments are simulated.

The Boom clay is a moderately swelling clay of Rupellian age. It is studied at the SCK · CEN in Belgium as a potential host rock for a radioactive waste repository. In this paper, suction-controlled experiments carried out on Boom clay by SCK · CEN are described. SCK · CEN has performed experiments to measure the relation between suction, water content and temperature and the relation between suction, stress and deformation. The applied suction-control techniques and experimental setups are detailed. The results of these experiments are discussed in the perspective of the model of Alonso and Gens. The influence of temperature on water uptake was rather small. The measured swelling-collapse behaviour can be explained by the Alonso and Gens model.     

Site effects at a vertical accelerometer array near Paducah, Kentucky

The downhole vertical accelerometer array VSAP near Paducah, KY, consists of three-component accelerometers at the surface, the top of the McNairy Formation (−41 m), and the top of the Paleozoic bedrock (—99 m). The array is at the northern end of the Mississippi Embayment, and it was installed to verify the ground-motion modeling for the site, assuming a significant earthquake in the New Madrid Seismic Zone. Accelerograms from 4.2 and 2.0 m_b earthquakes were used to check aspects of the modeling pertaining to linear behavior of the soil column, and to review the soil column models derived by drilling and geotechnical methods and through the use of high-resolution P- and SH-wave seismic refraction in reflection techniques. Results of the study indicate that for the linear case the soil column models derived by the two techniques are equivalent, and that the most important boundary in the soil column, with respect to amplification of the ground motions, is the interface between the limestone bedrock and soil.     

Mineral inclusions in diamonds from the Sputnik kimberlite pipe, Yakutia

The Sputnik kimberlite pipe is a small “satellite” of the larger Mir pipe in central Yakutia (Sakha), Russia. Study of 38 large diamonds (0.7-4.9 carats) showed that nine contain inclusions of the eclogitic paragenesis, while the remainder contain inclusions of the peridotitic paragenesis, or of uncertain paragenesis. The peridotitic inclusion suite comprises olivine, enstatite, Cr-diopside, chromite, Cr-pyrope garnet (both lherzolitic and harzburgitic), ilmenite, Ni-rich sulfide and a Ti-Cr-Fe-Mg-Sr-K phase of the lindsleyite-mathiasite (LIMA) series. The eclogitic inclusion suite comprises omphacite, garnet, Ni-poor sulfide, phlogopite and rutile. Peridotitic ilmenite inclusions have high Mg, Cr and Ni contents and high Nb/Zr ratios; they may be related to metasomatic ilmenites known from peridotite xenoliths in kimberlite. Eclogitic phlogopite is intergrown with omphacite, coexists with garnet, and has an unusually high TiO₂ content. Comparison with inclusions in diamonds from Mir shows general similarities, but differences in details of trace-element patterns. Large compositional variations among inclusions of one phase (olivine, garnet, chromite) within single diamonds indicate that the chemical environment of diamond crystallisation changed rapidly relative to diamond growth rates in many cases. P-T conditions of formation were calculated from multiphase inclusions and from trace element geothermobarometry of single inclusions. The geotherm at the time of diamond formation was near a 35 mW/m² conductive model; that is indistinguishable from the Paleozoic geotherm derived by studies of xenoliths and concentrate minerals from Mir. A range of Ni temperatures between garnet inclusions in single diamonds from both Mir and Sputnik suggests that many of the diamonds grew during thermal events affecting a relatively narrow depth range of the lithosphere, within the diamond stability field. The minor differences between inclusions in Mir and Sputnik may reflect lateral heterogeneity in the upper mantle.     

Environmental and legal aspects of karst areas

 Environmental impacts on karst settings are common as they are more sensitive than those of other rock terrains. Regulatory procedures that are effective in other rock terrains are not necessarily applicable to karst settings. Development and exploitation by man that affect the karst hydrology regime can trigger catastrophic events and result in numerous legal actions where the effects of changes go beyond property boundaries. A great variety of regulations and examples of litigation exist for karst areas. Received: 25 October 1994 · Accepted: 17 October 1995     

Ore lead and sulphur isotope pattern from the low-temperature veins of the Catalonian Coastal Ranges (NE Spain)

The lead isotopic composition of galenas from low-temperature veins of the Catalonian Coastal Ranges (NE Spain) displays a small but distinct range of values: 18.318 ≤²⁰⁶Pb/²⁰⁴Pb ≤ 18.678; 15.678 ≤²⁰⁷Pb/²⁰⁴Pb ≤ 15.767 and 38.534 ≤ ²⁰⁸Pb/²⁰⁴Pb ≤ 39.189. The data define a linear trend in the ²⁰⁸Pb/²⁰⁴Pb versus ²⁰⁶Pb/²⁰⁴Pb diagram that can be interpreted in terms of mixing of lead from different reservoirs. The lead appears to be derived from a solution that preferentially leached granites of Late-Hercynian age. Sulphur isotopic compositions of the galenas range from −8.5 to + 7.3‰ but variations in a single deposit are small (1–2‰), pointing to different sources of sulphur. In a δ³⁴S versus ²⁰⁶Pb/²⁰⁴Pb plot no clear correlation trends appear. Lead isotope data from deposits of similar characteristics situated within the Hercynian orogenic belt as the Linares-La Carolina district (S Spain) and the Inglesiente-Sulcis area (Sardinia) show slight differences indicating a source area of similar lead isotopic composition. In the CCR and Linares-La Carolina districts, the main source of lead is related to the Hercynian granites of each area, although small contributions from metasediments cannot be ruled out. The contribution of lead from metasedimentary rocks is particularly seen in the Sardinian (Hercynian vein-type) ores. Received: 8 August 1996 / Accepted: 5 December 1996     

Harmonic Analysis of Sedimentary Cyclic Sequences in Kansas, Midcontinent, USA

Several stratigraphic sequences in the Upper Carboniferous (Pennsylvanian) in Kansas (Midcontinent, USA) were analyzed quantitatively for periodic repetitions. The sequences were coded by lithologic type into strings of datasets. The strings then were analyzed by an adaptation of a one-dimensional Fourier transform analysis and examined for evidence of periodicity. The method was tested using different states in coding to determine the robustness of the method and data. The most persistent response is in multiples of 8–10 ft (2.5–3.0 m) and probably is dependent on the depositional thickness of the original lithologic units. Other cyclicities occurred in multiples of the basic frequency of 8–10 with persistent ones at 22 and 30 feet (6.5–9.0 m) and large ones at 80 and 160 feet (25–50 m). These levels of thickness relate well to the basic cyclothem and megacyclothem as measured on outcrop. We propose that this approach is a suitable one for analyzing cyclic events in the stratigraphic record.     

Trough structures in the Western syenite of Kûngnât, S Greenland: mineralogy and mechanism of formation

In the Western layered syenite stock of the Kûngnât complex, close to the intrusion wall, there are two groups of regularly spaced, horizontally aligned, stacks of infilled troughs. These features resemble, to some extent, the trough layering within the UZa trough zone of the gabbroic Skaergaard¹ intrusion. Such features have never before been described in syenitic rocks. The lower group comprises six stacks of infilled troughs; the upper group comprises seven. The two groups are vertically separated by about 10 m of sparsely layered syenite. The infilled troughs in the upper group are stacked vertically; those in the lower group gradually shift laterally and appear to converge over an antiformal structure. It is not possible to measure the plunge of the troughs due to outcrop geometry. Troughs are up to about 80 cm deep. The material filling the troughs is more melanocratic than the syenite that surrounds them. Individual troughs exhibit normal grading. The syenite within the troughs contains large amphibole and biotite oikocrysts. Faintly layered syenite occurs laterally between stacks. Minerals in the host syenite are more ferroan than those in the troughs, this is thought to represent re-equilibration with residual melt. Feldspar textures, irregular pyroxene and apatite zoning, and the presence of irregular shaped amphibole patches within pyroxene grains all imply alteration by a pervasive deuteric fluid in both troughs and host syenite. We suggest the following hypothesis for the formation of the troughs. A mafic crystal sheath which grows on the walls of the magma chamber becomes detached and flows downwards. The descending slurry forms a series of regularly spaced streams which erode trough shaped channels into the crystal pile at the base of the chamber. Uneroded crystal material separates the channels. Crystals are deposited from these streams. Successive streams flow over each other to generate stacks of troughs.     

Erosion rates of alpine bedrock summit surfaces deduced from in situ Be and Al

We have measured the concentration of in situ produced cosmogenic ¹⁰Be and ²⁶Al from bare bedrock surfaces on summit flats in four western U.S. mountain ranges. The maximum mean bare-bedrock erosion rate from these alpine environments is 7.6 ± 3.9 m My⁻¹. Individual measurements vary between 2 and 19 m My⁻¹. These erosion rates are similar to previous cosmogenic radionuclide (CRN) erosion rates measured in other environments, except for those from extremely arid regions. This indicates that bare bedrock is not weathered into transportable material more rapidly in alpine environments than in other environments, even though frost weathering should be intense in these areas. Our CRN-deduced point measurements of bedrock erosion are slower than typical basin-averaged denudation rates ( 50 m My⁻¹). If our measured CRN erosion rates are accurate indicators of the rate at which summit flats are lowered by erosion, then relief in the mountain ranges examined here is probably increasing.

We develop a model of outcrop erosion to investigate the magnitude of errors associated with applying the steady-state erosion model to episodically eroding outcrops. Our simulations show that interpreting measurements with the steady-state erosion model can yield erosion rates which are either greater or less than the actual long-term mean erosion rate. While errors resulting from episodic erosion are potentially greater than both measurement and production rate errors for single samples, the mean value of many steady-state erosion rate measurements provides a much better estimate of the long-term erosion rate.     

Artificial periodic irregularities in the auroral ionosphere

Artificial periodic irregularities (API) are produced in the ionospheric plasma by a powerful standing electromagnetic wave reflected off the F region. The resulting electron-density irregularities can scatter other high-frequency waves if the Bragg scattering condition is met. Such measurements have been performed at mid-latitudes for two decades and have been developed into a useful ionospheric diagnostic technique. We report here the first measurements from a high-latitude station, using the EISCAT heating facility near Troms0, Norway. Both F-region and lower-altitude ionospheric echoes have been obtained, but the bulk of the data has been in the E and D regions with echoes extending down to 52-km altitude. Examples of API are shown, mainly from the D region, together with simultaneous VHP incoherent-scatter-radar (ISR) data. Vertical velocities derived from the rate of phase change during the irregularity decay are shown and compared with velocities derived from the ISR. Some of the API-derived velocities in the 75–115-km height range appear consistent with vertical neutral winds as shown by their magnitudes and by evidence of gravity waves, while other data in the 50–70-km range show an unrealistically large bias. For a comparison with ISR data it has proved difficult to get good quality data sets overlapping in height and time. The initial comparisons show some agreement, but discrepancies of several metres per second do not yet allow us to conclude that the two techniques are measuring the same quantity. The irregularity decay time-constants between about 53 and 70 km are compared with the results of an advanced ion-chemistry model, and height profiles of recorded signal power are compared with model estimates in the same altitude range. The calculated amplitude shows good agreement with the data in that the maximum occurs at about the same height as that of the measured amplitude. The calculated time-constant agrees very well with the data below 60 km but is larger above 60 km by a factor of up to 2 at 64 km. The comparisons with the model are considered to be a good basis for more refined comparisons.