The use of MRS in the determination of hydraulic transmissivity: The case of alluvial aquifers

Magnetic Resonance Sounding (MRS) is nowadays accepted as a new geophysical method that can be used for a reliable determination of the ground water content distribution in the top 150 m. A great effort has also been made in MRS development to deduce the hydraulic transmissivity, based on empiric relationships of the permeability with a factor F which is calculated with NMR parameters measured at laboratory scale. To use this relationship under field conditions a calibration coefficient C_T = T_pt / F has to be previously established, which demands the knowledge of the transmissivity T_pt evaluated in the pumping test. The transmissivity can then be calculated at any other site of the same aquifer using the relation T_mrs = C_TF. The C_T values reported suggest a certain relationship with the lithology, but with a great dispersion and contradictory results. MRS surveys carried out in alluvial aquifers in Spain have shown that the value of C_T evaluated at one site may not be valid at another place of the same aquifer, because of the great heterogeneity of this kind of geological environment. The demand of a pumping test at each site where a MRS is measured invalidates the method actually used for MRS transmissivity evaluation. More than 50 MRS have been used to propose a new methodology. The aquifers visited cover a great range of transmissivities (from 2 × 10^− 6 to 9 × 10^− 3 m²/s). The MRS signal amplitude varies between 20 and 1400 nV, the signal/noise ratio is in the range from 0.6 to 42, and the value of the decay time constant varies from 200 to 800 ms. It has been demonstrated that when the transmissivity increases, the value of F decreases, and C_T increases, except for certain groups of MRS taken at the same aquifer or part of one aquifer, for which F increases with T_pt, keeping C_T constant. A function C_T(F) of the type C_T = mF^− n has been obtained that allows the transmissivity evaluation without the need of T_pt. Considering that both values of transmissivity, T_pt and T_mrs, are subjected to deviations due to the experimental errors as well as due to evaluation errors, the prediction achieved by the proposed equation is rather good. To perform a better evaluation of the values of the coefficients m and n it is necessary to have a greater number of MR soundings of good quality and with a trustworthy inversion at locations where a really comparable and good performed pumping test is available, covering a sufficient range of transmissivities. Though the data we have used do not always fulfil these conditions, the result is promising. Once a trustable function is available, the forecast of the transmissivity using MRS will not need the existence of any pumping test in the area. The general extension of this methodology demands the availability of MRS taken at all kinds of geological and hydrogeological environments, which is impossible without the existence of a universal MRS data base.     

Hydrodynamics of floodplain wetlands in a chalk catchment: The River Lambourn, UK

Variations in floodplain channel water levels and valley floor groundwater levels (measured in piezometers and boreholes) are examined at selected points along the course of the River Lambourn, a chalk river in southern England. A local alluvial gravel aquifer in the valley bottom is associated with numerous small wetlands that extend over much of the river's perennial profile. Variations in hydraulic gradient between local borehole levels and/or floodplain channel water levels are described for three sites in the seasonal section of the channel at Bockhampton, East Garston and West Shefford. The results indicate that observed groundwater levels are closely associated with flows from discrete springs at the margins of the channel and floodplain. However, as the floodplain widens and the alluvial gravel aquifer increases in size, the gravel aquifer accounts for a substantial down-valley component of groundwater flow with a diffuse vertical water flux. In the lower catchment, the exchange of flows between the gravel aquifer and the river enables some attenuation of floodplain water-table variability, providing a stable hydrological regime for valley-bottom wetlands. Catchment controls upon the local, valley-bottom, wetland regime are demonstrated with the application of a simple groundwater model developed using MODFLOW. The model is used to simulate groundwater discharge to the river in the upper and lower catchment, in addition to the water level regime at selected points in the valley bottom in the lower catchment. The results demonstrate the importance of taking catchment-scale water flow into account when managing isolated wetlands in a permeable catchment.     

A master equation for reactive solute transport in porous media

The mean value of a density of a cloud of points described by a generalized Liouville equation associated with a convection dispersion equation governing adsorbing solute transport yields a joint concentration probability density. The general technique can be applied for either linear or nonlinear adsorption; here the application is restricted to linear adsorption in one-dimensional transport. The equation generated for the joint concentration probability density is in the general form of a Fokker-Planck equation, but with a suitable coordinate transformation, it is possible to represent it as a diffusion equation with variable coefficients.     

Lithologic control of bedrock meander dimensions in the appalachian valley and ridge province: A comment on a comment

Abrahams' comment relates meander length to channel cross-section shape and recurrence interval so that meander length can be both directly and inversely proportional to rock resistance. This reply notes that either meander length is directly proportional to rock resistance or it is not; one cannot have it both ways. Many Appalachian Valley and Ridge bedrock meanders are shown to be the same size as alluvial meanders, and appear to be somewhat underfit. A hypothesis is proposed where modest discharge increases may have accelerated bedrock meander cutting although present streams remain capable of slowly cutting the meanders.     

Flow and transport through two-dimensional isotropic media of binary conductivity distribution. Part 1: NUMERICAL methodology and semi-analytical solutions

Flow and transport take place in a heterogeneous medium made up from inclusions of conductivity K submerged in a matrix of conductivity K ₀. We consider two-dimensional isotropic media, with circular inclusions of uniform radii, that are placed at random and without overlap in the matrix. The system is completely characterized by the conductivity contrast =K/K ₀ and by the volume fraction n. The flow is uniform in the mean, of velocity U=const. The derivation of the velocity field is achieved by a numerical method of high accuracy, based on analytical elements. Approximate analytical solutions are derived by a few methods: composite elements, effective medium, dilute systems and first-order approximation in logconductivity variance. The latter was employed by Rubin (1995), while the dilute system approximation was used by Eames and Bush (1999) and Dagan and Lessoff (2001). Transport is solved in a Lagrangean framework, with trajectories determined numerically from the velocity field, by particle tracking. Results for the velocity variance and for the longitudinal macrodispersivity, for a few values of and n, are presented in Part 2.     

Flow and transport through two-dimensional isotropic media of binary conductivity distribution. Part 2: NUMERICAL simulations and comparison with theoretical results

In the present part the results of numerical simulations of flow and transport in media made up from circular inclusions of conductivity K that are submerged in a matrix of conductivity K ₀, subjected to uniform mean velocity, are presented. This is achieved for a few values of =K/K ₀ (0.01, 0.1 and 10) and of the volume fraction n (0.05, 0.1 and 0.2). The numerical simulations (NS) are compared with the analytical approximate models presented in Part 1: the composite elements (CEA), the effective medium (EMA), the dilute system (DSA) and the first-order in the logconductivity variance (FOA). The comparison is made for the longitudinal velocity variance and for the longitudinal macrodispersivity. This is carried out for n<0.2, for which the theoretical and simulation models represent the same structure of random and independent inclusions distribution. The main result is that transport is quite accurately modeled by the EMA and CEA for low , for which _L is large, whereas in the case of =10, the EMA matches the NS for n<0.1. The first-order approximation is quite far apart from the NS for the values of examined . This material is based upon work supported by the National Science Foundation under Grant No. 0218914. Authors also wish to thank the Center of Computational Research, University at Buffalo for assistance in running numerical simulations.     

Origin of the high variability of water mineral content in the bedrock aquifers of Southern Madagascar

This study assesses the causes of the high spatial variability of the mineral content of groundwater in crystalline bedrock of Southern Madagascar. Although many kilometres from the coast and at a mean altitude of 400 m a.s.l, wells drilled in this area produce water with electrical conductivities in the range of 300–30,000 μS cm⁻¹ with a high spatial variability. Chemical and isotopic data are used to identify the processes involved in the groundwater mineralization. It is shown that the chemical composition of the groundwater in this region has its origin in (i) normal silicate and carbonate weathering reactions and (ii) input of marine salts, probably via rainfall recharge, modified by evapo-concentrative processes probably including precipitation and re-dissolution of secondary evaporites in the unsaturated zone. To obtain a better understanding of the spatial salinity distribution, well parameters such as yields, weathered zone thickness, weathered materials and morphological positions (upper slope, mid-slope, lower slope or valley bottom) are scrutinized.

A correlation was found between high salinity and low flow, shallow groundwater environments (flat hill tops, valley bottoms, weakly developed and clayey weathered zones) and between low salinity and high flow environments (granular, well-developed weathered zones and situation on valley slopes).     

Comment on “Luminescence chronology of the Sankosh river terraces in the Assam–Buthan foothills of the Himalayas: Implications to climate and tectonics” by Quaternary Geochronology 72, 101364

In a study of which the main objective was to assess the impact of climate change and tectonics on the formation of river terraces along the Sankosh River, the eastern foreland of the Himalayas, the authors obtained geochrono-logical data using luminescence technique. Four strath river terraces (T4–T1) were distinguished within the valley bottom, and alluvial sediments from three terraces (T4, T2 and T1) were dated to the age range from 143 to 14 ka. The alluvial mantels of river terraces were then linked to the scheme of glacial-interglacial cycle. The paper, however, suffers from few inconsistent and missing information, and the assessment on the landscape evolution of the river valley is incomprehensive. The authors relied on feldspars and therefore the IRSL method was used, but “OSL ages” are discussed at the end. They state that (i) tectonics creates space for sediment accommodation and (ii) the luminescence ages flank deposition and incision phases. Despite the fact that the statements are questionable, the ways these happen are not elucidated within the paper. Due to the lack of geochronological data for T3 terrace, its formation is very enigmatic. It is even more mysterious in the light of the data for the other terraces, but the authors made no attempt to explain this riddle. However, that certain inability of reconstructing the history of T3 would shed a shadow on the robustness of the ages obtained.     

Environmental reconstructions of a late devensian terrace sequence. Some preliminary findings

Facies analysis of Severn Main terrace sediments at Eardington, Shropshire is instructive in determining aspects of the sedimentary environment at a ‘proximal’ location in the terrace. Evaluation of the structural and directional properties of the gravels indicates that the terrace is a composite feature comprising two units. The lower unit is notable for the widespread occurrence of large sandy facies, a near absence of massive gravel facies, and the presence of large, locally derived, lithoclasts in excess of 1 m in diameter. Facies associations observed indicate the existence of in-channel bars and large channels with bedforms at the dune-plane bed transition. The directional properties of the unit indicate flows parallel to the main valley axis and although the unit is unlike previously reported examples it is interpreted as being the product of a proximal low sinuosity environment. The upper unit is composed mainly of multistorey, coarse, massive gravel units with limited interbedded sand facies. Fabric analysis indicates progradation of the unit from the west. The facies sequence is similar to those encompassed by facies models proposed for proximal alluvial fan deposits and therefore, on structural and directional grounds, the upper unit is interpreted as an alluvial fan deposit which prograded into the main Severn valley from the adjacent Mor Brook tributary. The implications of the composite nature of the terrace at this point are considered both in terms of terrace correlation and palaeohyrological estimation.     

EFFECTS OF RIVER TRAINING STRUCTURES ON THE SEDIMENT BALANCE OF A SAND - BEARING RIVER

1. INTROOrCTIOXThe upper Tisza in Hungary is a sand bed alluvial river. presenting several river training problems. In order to solve these problems first the laws of free river--bed evolution must be clarified thenthe effects of existing river training s…