Large to intermediate-scale aquifer heterogeneity in fine-grain dominated alluvial fans (Cenozoic As Pontes Basin,northwestern Spain): insight based on three-dimensional geostatistical reconstruction

Facies reconstructions are used in hydrogeology to improve the interpretation of aquifer permeability distribution. In the absence of sufficient data to define the heterogeneity due to geological processes, uncertainties in the distribution of aquifer hydrofacies and characteristics may appear. Geometric and geostatistical methods are used to understand and model aquifer hydrofacies distribution, providing models to improve comprehension and development of aquifers. However, these models require some input statistical parameters that can be difficult to infer from the study site. A three-dimensional reconstruction of a kilometer scale fine-grain dominated Cenozoic alluvial fan derived from more than 200 continuously cored, closely spaced, and regularly distributed wells is presented. The facies distributions were reconstructed using a genetic stratigraphic subdivision and a deterministic geostatistical algorithm. The reconstruction is only slightly affected by variations in the geostatistical input parameters because of the high-density data set. Analysis of the reconstruction allowed identification in the proximal to medial alluvial fan zones of several laterally extensive sand bodies with relatively higher permeability; these sand bodies were quantified in terms of volume, mean thickness, maximum area, and maximum equivalent diameter. These quantifications provide trends and geological scenarios for input statistical parameters to model aquifer systems in similar alluvial fan depositional settings.     

Surface complexation modeling for predicting solid phase arsenic concentrations in the sediments of the Mississippi River Valley alluvial aquifer, Arkansas, USA

The potential health impact of As in drinking water supply systems in the Mississippi River Valley alluvial aquifer in the state of Arkansas, USA is significant. In this context it is important to understand the occurrence, distribution and mobilization of As in the Mississippi River Valley alluvial aquifer. Application of surface complexation models (SCMs) to predict the sorption behavior of As and hydrous Fe oxides (HFO) in the laboratory has increased in the last decade. However, the application of SCMs to predict the sorption of As in natural sediments has not often been reported, and such applications are greatly constrained by the lack of site-specific model parameters. Attempts have been made to use SCMs considering a component additivity (CA) approach which accounts for relative abundances of pure phases in natural sediments, followed by the addition of SCM parameters individually for each phase. Although few reliable and internally consistent sorption databases related to HFO exist, the use of SCMs using laboratory-derived sorption databases to predict the mobility of As in natural sediments has increased. This study is an attempt to evaluate the ability of the SCMs using the geochemical code PHREEQC to predict solid phase As in the sediments of the Mississippi River Valley alluvial aquifer in Arkansas. The SCM option of the double-layer model (DLM) was simulated using ferrihydrite and goethite as sorbents quantified from chemical extractions, calculated surface-site densities, published surface properties, and published laboratory-derived sorption constants for the sorbents. The model results are satisfactory for shallow wells (10.6 m below ground surface), where the redox condition is relatively oxic or mildly suboxic. However, for the deep alluvial aquifer (21-36.6 m below ground surface) where the redox condition is suboxic to anoxic, the model results are unsatisfactory.     

Sedimentological study of the Chernobyl NPP site to schematise radionuclide migration conditions

The approach, methods and results of a sedimentological study of a near-surface stratum of Late Pleistocene-Holocene deposits in the near-zone (5–10 km radius) of the Chernobyl Nuclear Power Plant are presented. Sedimentological analyses are carried out at three levels of detail: regional-, local- and object-scale. The unsaturated zone and unconfined aquifer at the site are composed of two main genetic types of deposits, aeolian and alluvial, including several dynamic facies. Consideration of lithological properties leads to following ranking of the main genetic sediments facies with respect to radionuclide migration retardation potential: aeolian <alluvial channel <alluvial overbank <alluvial abandoned channel. Based on sedimentological interpretations, the geological environment is schematised into a set of typical geological sections possessing different radionuclide retardation potential.     

Hydrogeochemical characterization of major factors affecting the quality of groundwater in southern Iran,Janah Plain

The Janah alluvial aquifer is located in southern Iran with an arid climate. The type of groundwater in this aquifer is dominantly of sodium chloride and total dissolved solid of groundwater samples range from 1.63 to 335 g/L which confirms that groundwater quality has been severely degraded by salinization. Hydrogeochemical and isotopic investigations were conducted to identify the source of salinity. Total dissolved solids and major ion concentrations were measured at 51 selected sampling sites including springs, wells and surface waters. In addition stable isotopic composition (oxygen-18 and deuterium) was measured in 6 sampling points.The study indicates that the sources of salinity of the Janah aquifer include dissolution of salt diapir and evaporite rocks, a geothermal spring and intrusion of the river water which function individually or together in different parts of the aquifer. Based on the hydrogeochemical and geological studies conceptual flow models were prepared for different parts of the aquifer which illustrate how each source of salinity deteriorates the quality of the alluvial aquifer. We proposed few remediation methods including construction of cemented channel and sealed basins to improve groundwater quality. These methods would prevent infiltration of low quality water into the alluvial aquifer.     

Assessment of naturally occurring arsenic contamination in the groundwater of Sarkisla Plain (Sivas/Turkey)

High arsenic levels in groundwater of the aquifers, belonging to the Pliocene terrestrial layers and Quaternary alluvial sediments, have become a significant problem for the inhabitants living in Sarkisla (Turkey). The main objective of this study was to determine the origin and arsenic contamination mechanisms of the Sarkisla drinking water aquifer systems. The highest arsenic concentrations were found in Pliocene layers and alluvial sediments with concentrations ranging from 2.1 to 155 mg/kg. These rocks are the main aquifers in the study area, and most of the drinking groundwater demand is met by these aquifers. Groundwater from the Pliocene aquifer is mainly Ca-HCO₃ and Ca-SO₄ water type with high EC values reaching up to 3,270 μS/cm, which is due to the sulfate dissolution in some parts of the alluvial aquifer. Stable isotope values showed that the groundwater was of meteoric origin. Tritium values for the groundwater were between 8.31 and 14.06 TU, representing a fast circulation in the aquifer. Arsenic concentrations in the aquifers were between 0.5 and 345 μg/L. The highest arsenic concentrations detected in the Pliocene aquifer system reached up to 345 μg/L with an average value of 60.38 μg/L. The arsenic concentrations of the wells were high, while the springs had lower arsenic concentrations. These springs are located in the upper parts of the study area where the rocks are less weathered. The hydrogeochemical properties demonstrated that the water–rock interaction processes in sulfide-bearing rocks were responsible for the remarkably high groundwater arsenic contamination in the study area. In the study area, the arsenic levels determined in groundwater exceeded the levels recommended by the WHO. Therefore, it is suggested that this water should not be used for drinking purposes and new water sources should be investigated.     

Approaches to surface complexation modeling of Uranium(VI) adsorption on aquifer sediments

Uranium(VI) adsorption onto aquifer sediments was studied in batch experiments as a function of pH and U(VI) and dissolved carbonate concentrations in artificial groundwater solutions. The sediments were collected from an alluvial aquifer at a location upgradient of contamination from a former uranium mill operation at Naturita, Colorado (USA). The ranges of aqueous chemical conditions used in the U(VI) adsorption experiments (pH 6.9 to 7.9; U(VI) concentration 2.5 · 10⁻⁸ to 1 · 10⁻⁵ M; partial pressure of carbon dioxide gas 0.05 to 6.8%) were based on the spatial variation in chemical conditions observed in 1999-2000 in the Naturita alluvial aquifer. The major minerals in the sediments were quartz, feldspars, and calcite, with minor amounts of magnetite and clay minerals. Quartz grains commonly exhibited coatings that were greater than 10 nm in thickness and composed of an illite-smectite clay with occluded ferrihydrite and goethite nanoparticles. Chemical extractions of quartz grains removed from the sediments were used to estimate the masses of iron and aluminum present in the coatings. Various surface complexation modeling approaches were compared in terms of the ability to describe the U(VI) experimental data and the data requirements for model application to the sediments. Published models for U(VI) adsorption on reference minerals were applied to predict U(VI) adsorption based on assumptions about the sediment surface composition and physical properties (e.g., surface area and electrical double layer). Predictions from these models were highly variable, with results overpredicting or underpredicting the experimental data, depending on the assumptions used to apply the model. Although the models for reference minerals are supported by detailed experimental studies (and in ideal cases, surface spectroscopy), the results suggest that errors are caused in applying the models directly to the sediments by uncertain knowledge of: 1) the proportion and types of surface functional groups available for adsorption in the surface coatings; 2) the electric field at the mineral-water interface; and 3) surface reactions of major ions in the aqueous phase, such as Ca²⁺, Mg²⁺, HCO₃⁻, SO₄²⁻, H₄SiO₄, and organic acids. In contrast, a semi-empirical surface complexation modeling approach can be used to describe the U(VI) experimental data more precisely as a function of aqueous chemical conditions. This approach is useful as a tool to describe the variation in U(VI) retardation as a function of chemical conditions in field-scale reactive transport simulations, and the approach can be used at other field sites. However, the semi-empirical approach is limited by the site-specific nature of the model parameters.     

Modelling of the construction of the Rhône alluvial plain since 15 000 years BP

A numerical model presented here develops a three-dimensional image of alluvial media on an elementary scale significant for groundwater flow modelling. The model was tested on the alluvial plain of the Rhône River (France), on a scale of several kilometres and, from geomorphological observations and dating, reproduced the construction of this alluvial plain from ≈15 000 years BP to the present. The history of the alluvial plain during the Late Glacial and Holocene periods is summarized. Through most of this time, the River Rhône has maintained a braided pattern, with the exception of two incising phases with a meander pattern. The model does not use any physically based equations or water representation. The main processes governing the construction of the plain are modelled by simple rules chosen according to geometrical or empirical laws taken from the literature or as modelling assumptions. Using multi-agent concepts of distribution and interaction of elementary entities, these sedimentary rules are applied to 'sediment' entities or to conceptual 'erosion' entities that simulate local deposition and erosion of sediments. The sedimentation model reproduces the various climatic periods during which the sediments were deposited by simulating genetic periods and associated modelled processes. For each period, the model was constrained by quantitative field data such as altitude of ancient channels and deposits or thickness of sediments. The general geometry of the alluvial deposits was satisfactorily reproduced. During the simulation, characteristic large-scale features emerge despite the use of local rules. The model results are discussed with reference to other approaches, such as geostatistical or Boolean models, and the applicability of the model to other less documented alluvial plains is outlined.     

Application of tritium in precipitation and in groundwater of the Kouris catchment (Cyprus) for description of the regional groundwater flow

The Kouris catchment is located in the south of the Troodos massif in Cyprus. It constitutes one of the biggest catchments of the island with important freshwater resources. Geologically, the catchment includes an ophiolitic complex outcropping in the north which is overlaid by sedimentary rocks in the south. The hydrology is driven by a Mediterranean climate, a mountainous topography, and a complex distribution of the hydrogeological properties resulting from the complex geology.To improve the understanding of groundwater hydrology of the Kouris catchment, 176 groundwater and precipitation samples were collected and their ³H contents were analyzed. The three-dimensional ³H transport in the groundwater was simulated by the PMPATH code. For numerical modelling, a regional input function of ³H in precipitation was constructed from a linear regression between data for Cyprus and for neighboring meteorological stations. The calculated residence times for the groundwaters in the sedimentary aquifer and Pillow Lavas were greater than 48 a and were considerably greater than those of the ophiolitic complex (14–30 a). The calibrated aquifer porosities were in a range of 0.05–0.06. The PMPATH model was applied for delineation of spring catchments that were represented by quite narrow zones of lengths up to 5 km.Another contribution resulting from the ³H analysis was a better understanding of the river–aquifer interactions. In most of the southern part, the lithified sediments received only negligible amounts of water from the rivers, while the alluvial aquifer contained mostly water infiltrated from rivers. The largest springs in the southern part, associated with the alluvial aquifer, also discharged water identical to that in the rivers.     

Multi-scale characterization of a complex karst and alluvial aquifer system in southern Germany using a combination of different tracer methods

Water suppliers face major challenges such as climate change and population growth. To prepare for the future, detailed knowledge of water resources is needed. In southern Germany, the state water supplier Zweckverband Landeswasserversorgung provides 3 million people with drinking water obtained from a complex karst and alluvial aquifer system and the river Danube. In this study, a combination of different tracing techniques was used with the goal of a multi-scale characterization of the aquifer system and to gain additional knowledge about groundwater flow toward the extraction wells in the Danube Valley. For the small-scale characterization, selected groundwater monitoring wells were examined using single-borehole dilution tests. With these tests, a wide range of flow behavior could be documented, including fast outflow within just a few hours in wells with good connection to the aquifer, but also durations of many weeks in low-permeability formations. Vertical flow, caused by multiple flow horizons or uprising groundwater, was detected in 40% of the tested wells. A regional multi-tracer test with three injections was used to investigate the aquifer on a large scale. For the highly karstified connection between a swallow hole and a spring group, high flow velocities of around 80 m/h could be documented. Exceptionally delayed arrivals, 250 and 307 days after the injection, respectively showing maximum velocities of 0.44 and 0.39 m/h, were observed in an area where low-permeability sediments overlay the karst conduits. With the chosen methods, a distinct heterogeneity caused by the geological setting could be documented on both scales.

     

Olocenic alluvial aquifer of the River Cornia coastal plain (southern Tuscany, Italy): database design for groundwater management

 Hydrogeological research is in progress, utilizing GIS methods, with the principal aim of modelling the Olocenic alluvial aquifer of the River Cornia coastal plain (southern Tuscany, Italy), which has been exploited for drinking water, irrigation, and industrial uses. A consequence of exploitation has been the appearance of wide seawater intrusion. The alluvial aquifer has recently been subjected to new well fields for the supply of drinking water, with an increase of total average discharge of about 4×10⁶ m³/year. This paper presents results obtained from updating and integrating basic knowledge and structuring the database. The hydrogeological study allowed the recognition of the extension of areas that are characterized by a hydraulic head under the sea level, the progressive salinization of the aquifer, and the increase of water deficit in the aquifer which is produced by a progressive extraction of water superior to the natural recharge. In addition, benefits and disadvantages resulting from the location of new well fields in a hydrogeologically favourable zone, and the boundary conditions for much of the area studied have been defined. The GIS was used as support for making and updating the tabular and spatial database with the aim of integrating the local and regional hydrogeological knowledge. This study will permit the realization of a numerical simulation of the groundwater flow of the aquifer aimed at correcting the management of water resources, by means of the GIS-modelling integration. Received: 23 June 1998 · Accepted: 16 November 1998     

Subaqueous artesian springs and associated spring pits in a Himalayan pond

Subaqueous, bowl-shaped depressions found in a Himalayan pond formed in an abandoned river channel in the Lingti Valley (Spiti, NW India) are spring pits (Quirke 1930 ). The occurrence of the spring pits is restricted to the western end of the pond, where coarse-grained, highly permeable alluvial fan material continues below lacustrine mud deposits. The spring pits formed by active vertical discharge of ground water from an underlying artesian alluvial fan aquifer, confined by the overlying fine-grained lacustrine sediments. The aquifer is continuously recharged by down-slope ground-water flow in the alluvial fan. These small artesian springs are comparable with much larger artesian springs described in the literature and a similar mechanism of formation is proposed. Some similarities in their shapes and mechanisms of formation may indicate that spring pits represent small, nearshore examples of pockmarks. The differences of these features, formed by persistent fluidization from short-lived seismic liquefaction processes, are discussed and the utility of the structures for palaeo-environmental and palaeo-ground-water interpretation is evaluated.     

Hydrogeologic characteristics of the alluvial tuff aquifer of northern Sahand Mountain slopes, Tabriz, Iran

The Tabriz area is a densely populated area of northwestern Iran (more than 1.5 million in population) with a large proportion of its drinking, domestic, industrial and agricultural water supplied from groundwater resources. The average rate of drinking and industrial water use in the city of Tabriz is about 3.45 m³ s^–1. The Plio-Pleistocene unconfined alluvial tuff aquifer (about 1,275 km²), the most important aquifer in the area, has been known for many years as a reliable resource. The greatest estimated thickness of the alluvial tuff lies in the Saidabad area, with 350 m thickness. There are 994 deep and 284 shallow active pumping wells and 83 qanats operate in the alluvial tuff aquifer. The total water withdrawal from all these artificial discharge points has been measured at 72, 3.8 and 17 million m3/year, respectively. Analytical and numerical methods have been applied to the constant rate pumping test data from the Saidabad wellfield (eight pumping and three observation wells). The values of electrical conductivity in the groundwater of alluvial tuff aquifer range from 203 to 960 μS cm^–1 and bicarbonate type water dominates.     

Dealing with the Spatial Synthetic Heterogeneity of Aquifers in the North China Plain: A Case Study of Luancheng County in Hebei Province

The complexity of alluvial-pluvial fan depositional systems makes the detailed characterization of their heterogeneity difficult, yet such a detailed characterization is commonly needed for construction of reliable groundwater models. Traditional models mainly focus on using a single aquifer property to qualitatively or semi-quantitatively characterize the heterogeneity of aquifer, so that they are unable to quantitatively reflect the synthetic heterogeneity of all aquifer properties. In this paper, we propose the heterogeneity synthetic index (HSI) for quantitative characterization of synthetic heterogeneity of an aquifer. The proposed calculation process involves four steps: (1) estimation of the hydraulic conductivity of a sediment sample using the cloud-Markov model, (2) establishment of the sedimentary microfacies distribution model through the Markov chain, (3) characterization of the distribution model of hydrogeological parameters using the improved sequential simulation method according to the “facies-controlled modeling” technique, and (4) application of the entropy weight method to calculate the weight coefficient of the above aquifer properties. The HSI of an aquifer is calculated by superposition of these models according to the corresponding weight coefficient. This approach was applied to the Luancheng aquifer deposit in the southeast Hutuo River alluvial-pluvial fan in the North China Plain (NCP). The results have demonstrated that aquifer 3 which was formed in the middle Pleistocene has the strongest heterogeneity, with an HSI of 0.25–0.75. Aquifer 4 formed in the early Pleistocene shows an intermediate heterogeneity, with the HSI ranging 0.35–0.75. The weakest heterogeneity was found in aquifers 1 and 2 formed in the Holocene and late Pleistocene, with HSI values of 0.40–0.75 and 0.40–0.80, respectively. The heterogeneity of all the four aquifers is relatively strong in the radial direction of the Huai River alluvial-pluvial fan due to the abrupt change of microfacies. In contrast, in the radial direction of the Hutuo River alluvial-pluvial fan, the microfacies change mildly, and the continuity of hydrogeological parameters is better, which has resulted in weaker heterogeneity of the four aquifers in this direction. Findings suggest that the sedimentary environment has significant effects on the aquifer heterogeneity. Considering that there are many aquifer properties, HSI can quantitatively characterize the synthetic heterogeneity of the aquifer and describe the influence of each aquifer property on the synthetic heterogeneity of the aquifer according to its weight coefficient. Thus the HSI approach can be successfully used to deal with the spatial heterogeneity of aquifer and provide a foundation for studies on contaminant transport.     

Geological controls and effects of floodplain asymmetry on river–groundwater interactions in the southeastern Coastal Plain, USA

Channel sediment and alluvial aquifer hydraulic properties exert a major control on river–groundwater interactions. Channels and floodplains are often asymmetrical, resulting in differences in sediment hydraulic properties across the river. Floodplain asymmetry is common along Coastal Plain rivers in South Carolina and North Carolina, USA. The Tar River, North Carolina, has an asymmetrical valley. The study objective was to characterize the effects of floodplain asymmetry and geological controls on river–groundwater interactions. Floodplain and river channel sediments adjacent to the river were characterized with split spoon cores and hand auger samples along a 22-km reach. Hydrogeology was characterized with 38 piezometers and water level recorders in and adjacent to the river. Ground penetrating radar was used to define the shallow stratigraphy. Channel sediments were significantly different between the north and south sides of the river. Hydraulic conductivity and groundwater inputs were greater on the side of the river (north) that contained more permeable fluvial deposits. Groundwater chemistry (δ¹⁸O, specific conductance) data also suggested greater exchange between surface water and groundwater on the north side of the river channel. A conceptual hydrogeological model illustrates that groundwater movement and contaminant transport to the river differs across the channel due to asymmetrical geology.     

Daily scale modelling of aquifer–river connectivity in the urban alluvial aquifer in Langreo, Spain

Evaluating aquifer–river interactions is naturally complex, particularly within urban settings. This is largely due to the difficulties involved in quantifying most elements of the water balance. The ability of numerical models to deal with several dynamic variables simultaneously makes them valuable tools to address this kind of problem. An applied, modeling-based approach to investigate the spatial and temporal variations of aquifer–river connectivity within a shallow urban aquifer is presented. Model development is based on comprehensive field campaigns in Langreo, Spain. Two calibration runs (for summer and winter conditions) were carried out in order to evaluate the spatial distribution of recharge rates. The model suggests that baseflows are largely negligible in comparison with total streamflows. This is mostly attributed to the abrupt nature of the catchment, which prevents the existence of sufficiently large alluvial systems to a great extent. Modelling results also show that aquifer–river connectivity at the study site is constrained by urban pumping as well as by seasonal fluctuations.     

Investigation of the hydrogeochemical processes in an alluvial channel aquifer located in a typical Karoo Basin of Southern Africa

An investigation was conducted to assess the hydrogeochemical processes of an alluvial channel aquifer located in a typical Karoo Basin of Southern Africa. The investigation was aimed at identifying and describing the groundwater chemistry evolution and its contribution to the overall groundwater quality. X-ray fluorescent spectrometry (XRF) and X-ray diffractometry (XRD) analyses were performed on geological samples to identify and quantify the major element oxides and minerals. The study utilises the conventional Piper diagram, bivariate plots and PHREEQC hydrogeochemical model to analyse groundwater chemistry data obtained during the wet (February and May) and dry seasons (August and December) of 2011. The XRF and XRD results show that the channel deposits are dominated by SiO₂ element oxides and quartz minerals, thus elevated concentrations of silicon (Si⁴⁺) were found in the groundwater. Dolomite and calcite minerals were also detected in the unconsolidated aquifer sediments. The detailed study of the alluvial aquifer system has shown that dissolution of dolomite and calcite minerals and ion exchange are the dominant hydrogeochemical processes influencing the groundwater quality. The groundwater evolves from Ca²⁺–Mg²⁺–HCO₃ ^? recharge water that goes through ion exchange with Na⁺ in the clay-silt sediment to give a Na⁺–HCO₃ ^? water type. The groundwater is supersaturated with respect to quartz, dolomite and calcite minerals. The study shows the potential usefulness of simple bivariate plots as a complimentary tool to the conventional methods for analyzing groundwater hydrogeochemical processes.     

Alluvial aquifer recharge enhanced by a natural dam: feasibility assessment based on multidisciplinary characterization (Khuzestan,Southwest Iran)

In the Kushkak Valley (Khuzestan, Southwest Iran) an anticlinal structure has partially impounded an ephemeral stream. This natural impounded area has been chosen for an artificial recharge site due to its current rate of recharge, capability to store water and favorable situation for the construction of man-made barriers to stream flow. The aquifer to be recharged is the Kushkak unconfined aquifer which consists of medium to coarse-grained alluvial deposits that overlie consolidated conglomerate rock. In this semi-arid area with infrequent relatively heavy falls of rain, alluvial aquifer recharge can be an important process that sustains shallow, over-exploited groundwater bodies. In this investigation a multidisciplinary approach including: hydrometerological studies, and a detailed hydrogeochemical survey, have been carried out. Other essential prerequisite parameters for the scheme were also taken into account to determine the suitability of this location for groundwater artificial recharge. The assessment has brought out that (1) the proposed reservoir will conserve a major part of the water being lost, (2) annual runoff of about 0.27 MCM can be injected into the aquifer through recharge from impounded water, (3) hydrochemical data from surface water and from the Kushkak aquifer water demonstrates that dilution and change in compositional trend in the groundwater proximal to the impounded alluvial bed areas would be expected based on the infiltration capacity of this site, and (4) cost–benefit ratio of the project is 1:2 and it is assumed to recover the investment within six years.     

Changes in groundwater chemistry due to metallurgical activities in an alluvial aquifer in the Moa area (Cuba)

The aim of this work is to evaluate the changes in groundwater chemistry in an alluvial aquifer in the Moa area. Surface and ground water, metallurgical waste and various geological material samples were collected in order to evaluate groundwater composition. The results show that the alluvial aquifer is polluted with SO₄^2-, Mg²⁺and heavy metals. According to its major components in the alluvial aquifer, two types of groundwater are identified: magnesi-bicarbonated and sulphate-magnesic. Maximum SO₄^2– and Mg²⁺ contents are more than 1000 mg/L, and are four times higher than the acceptable levels for human consumption of water. The high values of Cr(VI), Ni(II), Mn(II) Fe_(total), SO₄^2– and Mg²⁺ in alluvial aquifers are due to polluted recharge from metallurgical waste from the tailing dam. This recharge is favoured by the preferential flow due to desiccation cracks in metallurgical waste. Geochemical modelling showed that potentially toxic heavy metals might exist largely in the forms of MSO₄^2– and M²⁺ in pore water of SAL metallurgical waste. All samples were supersaturated in goethite and hematite. Results from batch testing indicate that the heavy metals have two origins: natural, due to the existing ultramaphic rocks and laterites, and anthropogenic, by metallurgical waste rich in sulphate and (oxy)hydroxide minerals. These results highlight the need to locate and evaluate a new water source to supply the population of the city of Moa.     

Evaluation of a pumping test of the Snake River Plain aquifer using axial-flow numerical modeling

The Snake River Plain aquifer in southeast Idaho is hosted in a thick sequence of layered basalts and interbedded sediments. The degree to which the layering impedes vertical flow has not been well understood, yet is a feature that may exert a substantial control on the movement of contaminants. An axial-flow numerical model, RADFLOW, was calibrated to pumping test data collected by a straddle-packer system deployed at 23 depth intervals in four observation wells to evaluate conceptual models and estimate properties of the Snake River Plain aquifer at the Idaho National Engineering and Environmental Laboratory. A delayed water-table response observed in intervals beneath a sediment interbed was best reproduced with a three-layer simulation. The results demonstrate the hydraulic significance of this interbed as a semi-confining layer. Vertical hydraulic conductivity of the sediment interbed was estimated to be about three orders of magnitude less than vertical hydraulic conductivity of the lower basalt and upper basalt units. The numerical model was capable of representing aquifer conceptual models that could not be represented with any single analytical technique. The model proved to be a useful tool for evaluating alternative conceptual models and estimating aquifer properties in this application. Electronic Publication