Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Long-Period Variations of the Eccentricity Vector Valid also for Near Circular Orbits around a Non-Spherical Body

This paper studies the long period variations of the eccentricity vector of the orbit of an artificial satellite, under the influence of the gravity field of a central body. We use modified orbital elements which are non-singular at zero eccentricity. We expand the long periodic part of the corresponding Lagrange equations as power series of the eccentricity. The coefficients characterizing the differential system depend on the zonal coefficients of the geopotential, and on initial semi-major axis, inclination, and eccentricity. The differential equations for the components of the eccentricity vector are then integrated analytically, with a definition of the period of the perigee based on the notion of “free eccentricity”, and which is also valid for circular orbits. The analytical solution is compared to a numerical integration. This study is a generalization of (Cook, Planet. Space Sci., 14, 1966): first, the coefficients involved in the differential equations depend on all zonal coefficients (and not only on the very first ones); second, our method applies to nearly circular orbits as well as to not too eccentric orbits. Except for the critical inclination, our solution is valid for all kinds of long period motions of the perigee, i.e., circulations or librations around an equilibrium point.     

Investigating the mechanisms leading to the deglaciation of past continental northern hemisphere ice sheets with the CLIMBER–GREMLINS coupled model

A coupling procedure between a climate model of intermediate complexity (CLIMBER-2.3) and a 3-dimensional thermo-mechanical ice-sheet model (GREMLINS) has been elaborated. The resulting coupled model describes the evolution of atmosphere, ocean, biosphere, cryosphere and their mutual interactions. It is used to perform several simulations of the Last Deglaciation period to identify the physical mechanisms at the origin of the deglaciation process. Our baseline experiment, forced by insolation and atmospheric CO₂, produces almost complete deglaciation of past northern hemisphere continental ice sheets, although ice remains over the Cordilleran region at the end of the simulation and also in Alaska and Eastern Siberia. Results clearly demonstrate that, in this study, the melting of the North American ice sheet is critically dependent on the deglaciation of Fennoscandia through processes involving switches of the thermohaline circulation from a glacial mode to a modern one and associated warming of the northern hemisphere. A set of sensitivity experiments has been carried out to test the relative importance of both forcing factors and internal processes in the deglaciation mechanism. It appears that the deglaciation is primarily driven by insolation. However, the atmospheric CO₂ modulates the timing of the melting of the Fennoscandian ice sheet, and results relative to Laurentide illustrate the existence of threshold CO₂ values, that can be translated in terms of critical temperature, below which the deglaciation is impeded. Finally, we show that the beginning of the deglaciation process of the Laurentide ice sheet may be influenced by the time at which the shift of the thermohaline circulation from one mode to the other occurs.     

Examples of clustering around high redshift radio-galaxies

Published in Astrofizika, Vol. 38, No. 4, pp. 677–680, October–December, 1995.     

A 20 cm2 CCD mosaic camera for a dark matter search part I : Mechanics,optics and cryogeny

A 20cm ² CCD mosaic camera has been especially built to search for dark galactic halo objects by the gravitational microlensing effect. The sensitive area is made of 16 edge-buttable CCDs developped by Thomson-CTS, with 23×23 µm ² pixels. The 35 kg camera housing and mechanical equipment is presented. The associated electronics and data acquisition system are described in a separate paper. The camera resides at the focal plane of a 40 cm, f/10, Ferson reflector. The instrument has been in operation since December 1991 at the La Silla Observatory (ESO).     

WIMPs and solar evolution code

The Saclay solar evolution code is used to check the effect of WIMPs on solar evolution. In this paper we study the effects of various types of Cosmion-matter interactions, give constraints on the crosssections compatible with the measured neutrino rate of 2 SNU on chlorine, and relate these constraints to ongoing dark matter detection experiments.Unité associée au CNRS UA 280, F-75251 Paris Cedex 05, France.On leave from LPC, Collège de France.     

Geometry and kinematics of recent deformation in the Mondy–Tunka area (south-westernmost Baikal rift zone, Mongolia–Siberia)

We used satellite imagery and field data to investigate the south‐westernmost Baikal rift zone. We focus our study in the Mondy and Ikhe Ukhgun valleys, site of an Mw = 6.9 seismic event in 1950. Surface deformations are observed along the E–W‐trending Mondy strike‐slip fault and along the Ikhe Ukhgun thrust. The Mondy fault system is 80 km long and is composed of four segments 10–15 km long. These segments are characterized by subvertical planes with left‐lateral movements. The Ikhe Ukhgun thrust is 20 km long, dips 40° to the south and shows reverse movement with a left‐lateral component. These observations are consistent with the present‐day regional NNE–SSW compression and with the focal mechanism of the 1950 Mondy earthquake that was recently re‐evaluated. These features, like those observed in the Tunka basin, demonstrate a recent change of regional strain regime from transtension to transpression that we place before the Late Pleistocene.     

Wavelet analysis of inter‐annual variability in the runoff regimes of glacial and nival stream catchments,Bow Lake,Alberta

Continuous wavelet analyses of hourly time series of air temperature, stream discharge, and precipitation are used to compare the seasonal and inter‐annual variability in hydrological regimes of the two principal streams feeding Bow Lake, Banff National Park, Alberta: the glacial stream draining the Wapta Icefields, and the snowmelt‐fed Bow River. The goal is to understand how water sources and flow routing differ between the two catchments. Wavelet spectra and cross‐wavelet spectra were determined for air temperature and discharge from the two streams for summers (June–September) 1997–2000, and for rainfall and discharge for the summers of 1999 and 2000. The diurnal signal of the glacial runoff was orders of magnitude higher in 1998 than in other years, indicating that significant ice exposure and the development of channelized glacial drainage occurred as a result of the 1997–98 El Niño conditions. Early retreat of the snowpack in 1997 and 1998 led to a significant summer‐long input of melt runoff from a small area of ice cover in the Bow River catchment; but such inputs were not apparent in 1999 and 2000, when snow cover was more extensive. Rainfall had a stronger influence on runoff and followed quicker flow paths in the Bow River catchment than in the glacial catchment. Snowpack thickness and catchment size were the primary controls on the phase relationship between temperature and discharge at diurnal time scales. Wavelet analysis is a fast and effective means to characterize runoff, temperature, and precipitation regimes and their interrelationships and inter‐annual variability. The technique is effective at identifying inter‐annual and seasonal changes in the relative contributions of different water sources to runoff, and changes in the time required for routing of diurnal meltwater pulses through a catchment. However, it is less effective at identifying changes/differences in the type of the flow routing (e.g. overland flow versus through flow) between or within catchments. Copyright © 2003 John Wiley & Sons, Ltd.     

Groundwater flow and contaminant transport modelling at an air weapons range

Numerical modelling was done at the Cold Lake Air Weapons Range, Canada, to test whether the dissolved RDX and nitrate detected in groundwater come from the same sources, and to predict whether contamination poses a threat to the surface water receptors near the site. Military live fire training activities may indeed pose a risk of contamination to groundwater resources, however field investigations on military bases are quite recent, and little information is available on the long-term behaviour of munition residues related contaminants. Very limited information was available about the contaminant source zones, which were assigned based on our knowledge of current training activities. The RDX plume was well represented with the model, but the heterogeneous distribution of nitrate concentrations was more difficult to reproduce. It was nonetheless determined that both contaminants originate from the same areas. According to the model, both contaminants should reach the nearby river, but concentrations in the river should remain very low if the source zone concentration does not change. Finally, the model allowed the recommendation of a new location for the main bombing target, which would offer added protection to the river and the lake into which it flows.     

Effectiveness of sulfate-reducing passive bioreactors for treating highly contaminated acid mine drainage: I. Effect of hydraulic retention time

Sulfate-reducing passive bioreactors have proved to be an effective technology for the treatment of acid mine drainage (AMD) contaminated waters over relatively short periods of time (1–5 a). However, long-term efficiency can be limited by several factors including problems related to the hydraulic properties of the reactive mixture. In this study, the effect of two hydraulic retention times (HRTs) of 7.3 d and 10 d on the performance of passive bioreactors was evaluated over an 11-month period for the treatment of a highly contaminated AMD. Evolution of the porosity and hydraulic conductivity of the reactive mixture was also evaluated during the 15-month operation of two bioreactors. Results indicated that bioreactors were effective at both HRTs for increasing the pH and alkalinity of contaminated water and for SO₄ and metal removal (60–82% for Fe and up to 99.9% for Cd, Ni and Zn). Although the quality of treated effluent was significantly improved with the 10 d HRT compared to the 7.3 d HRT, results showed that the higher HRT reduced the porosity and the permeability of the reactive mixture which might lead to hydraulic related problems and, eventually, to limited efficiency in long-term operation compared to a shorter HRT. The choice of HRT for a passive bioreactor must therefore consider both the desired quality of treated effluent and the potential for deterioration of hydraulic properties in the reactive mixture.