Hydraulic Characteristics of Rough Fractures in Linear Flow under Normal and Shear Load

Summary A hydro-mechanical testing system, which is capable of measuring both the flow rates and the normal and shear displacement of a rock fracture, was built to investigate the hydraulic behaviour of rough tension fractures. Laboratory hydraulic tests in linear flow were conducted on rough rock fractures, artificially created using a splitter under various normal and shear loading. Prior to the tests, aperture distributions were determined by measuring the topography of upper and lower fracture surfaces using a laser profilometer. Experimental variograms of the initial aperture distributions were classified into four groups of geostatistical model, though the overall experimental variograms could be well fitted to the exponential model. The permeability of the rough rock fractures decayed exponentially with respect to the normal stress increase up to 5 MPa. Hydraulic behaviours during monotonic shear loading were significantly affected by the dilation occurring until the shear stress reached the peak strength. With the further dilation, the permeability of the rough fracture specimens increased more. However, beyond shear displacement of about 7 to 8 mm, permeability gradually reached a maximum threshold value. The combined effects of both asperity degradation and gouge production, which prohibited the subsequent enlargement of mean fracture aperture, mainly caused this phenomenon. Permeability changes during cyclic shear loading showed somewhat irregular variations, especially after the first shear loading cycle, due to the complex interaction from asperity degradations and production of gouge materials. The relation between hydraulic and mechanical apertures was analyzed to investigate the valid range of mechanical apertures to be applied to the cubic law. Received June 12, 2001; accepted February 26, 2002 Published online September 2, 2002     

Poroelastic analysis of cover-collapse sinkhole formation by piezometric surface drawdown

Where the water table is above the soil-rock contact in karst regions, cover-collapse sinkholes in the soil and soft sediment above the rock commonly occur as a result of drawdown of the piezometric surface in the karst aquifer. Transient stresses and pore pressures around soil voids at the soil-rock contact can cause hydraulic fracturing of the soil near the wall of the void. After the first such fracture, successive sloughing of soil propagates the soil void rapidly to the surface, resulting in a cover-collapse sinkhole. Sinkhole formation by this mechanism should be strongly a function of rate and magnitude of piezometric surface drawdown, permeability and tensile strength of the soil, and the size, depth, and geometry of the initial soil void. Large soil voids and those with walls that are partly planar or of low curvature are most susceptible to hydraulic fracture and the resulting progression to sinkhole formation.     

Trace element mobility during sub-seafloor alteration of basaltic glass from Ocean Drilling Program site 953 (off Gran Canaria)

Trace element concentrations of altered basaltic glass shards (layer silicates) and zeolites in volcaniclastic sediments drilled in the volcanic apron northeast of Gran Canaria during Ocean Drilling Program (ODP) leg 157 document variable element mobilities during low-temperature alteration processes in a marine environment. Clay minerals (saponite, montmorillonite, smectite) replacing volcanic glass particles are enriched in transition metals and rare earth elements (REE). The degree of retention of REE within the alteration products of the basaltic glass is correlated with the field strength of the cations. The high field-strength elements are preferentially retained or enriched in the alteration products by sorption through clay minerals. Most trace elements are enriched in a boundary layer close to the interface mineral-altered glass. This boundary layer has a key function for the physico-chemical conditions of the subsequent alteration process by providing a large reactive surface and by lowering the fluid permeability. The release of most elements is buffered by incorporation into secondary precipitates (sodium-rich zeolites, phillipsite, Fe- and Mn-oxides) as shown by calculated distribution coefficients between altered glasses and authigenic minerals. Chemical fluxes change from an open to a closed system behavior during prograde low-temperature alteration of volcaniclastic sediments with no significant trace metal flux from the sediment to the water column.     

Origin of iodine in volcanic fluids: I results from the Central American Volcanic Arc

The largest reservoir of crustal iodine is found in marine sediments, where it is closely associated with organic material. This presence, together with the existence of a long-lived, cosmogenic radioisotope ¹²⁹I (t_1/2 = 15.7 Ma), make this isotopic system well suited for the study of sediment recycling in subduction zones. Reported here are the results of ¹²⁹I/I ratios in volcanic fluids, collected during a comprehensive study of fluids and gases in the Central American Volcanic Arc. ¹²⁹I/I ratios, together with I, Br, and Cl concentrations, were determined in 79 samples from four geothermal centers and a number of crater lakes, fumaroles, hot springs, and surface waters in Costa Rica, Nicaragua, and El Salvador. Geothermal and volcanic fluids were found to have iodine concentrations substantially higher than values in seawater or meteoric waters. ¹²⁹I/I ratios in most of the geothermal fluids are below the preanthropogenic input ratio of 1500 × 10⁻¹⁵, demonstrating that recent anthropogenic additions are largely absent from the volcanic systems. The majority of the ¹²⁹I/I ratios are between 500 and 800 × 10⁻¹⁵. These ratios indicate minimum iodine ages between 25 and 15 Ma, in good agreement with the age of subducted sediments in this region. In all four geothermal systems, however, a few samples were found with iodine ages older than 40 Ma—that is, considerably below the expected age range for subducted sediments from the Cocos Plate. These samples probably reflect the presence of iodine derived from sediments in older accreted oceanic terraines. The iodine ages indicate that the magmatic end member for the volcanic fluids originates in the deeper parts of the subducted sediment column, with small additions from older iodine mobilized from the overlying crust. The high concentrations of iodine in geothermal fluids, combined with the observed iodine ages, demonstrate that remobilization in the main volcanic zone (and probably also in the forearc area) is an important part in the overall marine cycle of iodine and similar elements.     

The major-ion composition of silurian seawater

One-hundred fluid inclusions in Silurian marine halite were analyzed in order to determine the major-ion composition of Silurian seawater. The samples analyzed were from three formations in the Late Silurian Michigan Basin, the A-1, A-2, and B Evaporites of the Salina Group, and one formation in the Early Silurian Canning Basin (Australia), the Mallowa Salt of the Carribuddy Group. The results indicate that the major-ion composition of Silurian seawater was not the same as present-day seawater. The Silurian ocean had lower concentrations of Mg²⁺, Na⁺, and SO₄²⁻, and much higher concentrations of Ca²⁺ relative to the ocean’s present-day composition. Furthermore, Silurian seawater had Ca²⁺ in excess of SO₄²⁻. Evaporation of Silurian seawater of the composition determined in this study produces KCl-type potash minerals that lack the MgSO₄-type late stage salts formed during the evaporation of present-day seawater. The relatively low Na⁺ concentrations in Silurian seawater support the hypothesis that oscillations in the major-ion composition of the oceans are primarily controlled by changes in the flux of mid-ocean ridge brine and riverine inputs and not global or basin-scale, seawater-driven dolomitization. The Mg²⁺/Ca²⁺ ratio of Silurian seawater was ∼1.4, and the K⁺/Ca²⁺ ratio was ∼0.3, both of which differ from the present-day counterparts of 5 and 1, respectively. Seawaters with Mg²⁺/Ca²⁺ <2 facilitate the precipitation of low-magnesian calcite (mol % Mg < 4) marine ooids and submarine carbonate cements whereas seawaters with Mg²⁺/Ca²⁺ >2 (e.g., modern seawater) facilitate the precipitation of aragonite and high-magnesian calcite. Therefore, the early Paleozoic calcite seas were likely due to the low Mg²⁺/Ca²⁺ ratio of seawater, not the pCO₂ of the Silurian atmosphere.     

Physical and chemical properties of the waters of saline lakes and their importance for deep-water renewal: Lake Issyk-Kul, Kyrgyzstan

The relationships between electrical conductivity, temperature, salinity, and density are studied for brackish Lake Issyk-Kul. These studies are based on a newly determined major ion composition, which for the open lake shows a mean absolute salinity of 6.06 g kg⁻¹. The conductivity-temperature relationship of the lake water was determined experimentally showing that the lake water is about 1.25 times less conductive than seawater diluted to the same absolute salinity as that of the lake water. Based on these results, an algorithm is presented to calculate salinity from in-situ conductivity measurements. Applied to the field data, this shows small but important vertical salinity variations in the lake with a salinity maximum at 200 m and a freshening of the surface water with increasing proximity to the shores. The algorithm we adopt to calculate density agrees well with earlier measurements and shows that at 20°C and 1 atm Lake Issyk-Kul water is about 530 g m⁻³ denser than seawater at the same salinity. The temperature of maximum density at 1 atm is about 0.15°C lower than that for seawater diluted to the same salinity. Despite its small variations, salinity plays an important role, together with temperature changes, in the static stability and in the production of deep-water in this lake. Changes in salinity may have had important consequences on the mixing regime and the fate of inflowing river water over geological time. Uncharged silicic acid is negligible for the stability of the water column except near an ∼15 m thick nepheloid layer observed at the bottom of the deep basin.     

Lithosphere structure of European sedimentary basins from regional three-dimensional gravity modelling

A three-dimensional (3D) density model, approximated by two regional layers—the sedimentary cover and the crystalline crust (offshore, a sea-water layer was added), has been constructed in 1° averaging for the whole European continent. The crustal model is based on simplified velocity model represented by structure maps for main seismic horizons—the “seismic” basement and the Moho boundary. Laterally varying average density is assumed inside the model layers. Residual gravity anomalies, obtained by subtraction of the crustal gravity effect from the observed field, characterize the density heterogeneities in the upper mantle. Mantle anomalies are shown to correlate with the upper mantle velocity inhomogeneities revealed from seismic tomography data and geothermal data. Considering the type of mantle anomaly, specific features of the evolution and type of isostatic compensation, the sedimentary basins in Europe may be related into some groups: deep sedimentary basins located in the East European Platform and its northern and eastern margins (Peri-Caspian, Dnieper–Donets, Barents Sea Basins, Fore–Ural Trough) with no significant mantle anomalies; basins located on the activated thin crust of Variscan Western Europe and Mediterranean area with negative mantle anomalies of −150 to −200×10⁻⁵ ms⁻² amplitude and the basins associated with suture zones at the western and southern margins of the East European Platform (Polish Trough, South Caspian Basin) characterized by positive mantle anomalies of 50–150×10⁻⁵ ms⁻² magnitude. An analysis of the main features of the lithosphere structure of the basins in Europe and type of the compensation has been carried out.