Environmental tracers as indicators of karst conduits in groundwater in South Dakota,USA

Environmental tracers sampled from the carbonate Madison aquifer on the eastern flank of the Black Hills, South Dakota, USA indicated the approximate locations of four major karst conduits. Contamination issues are a major concern because these conduits are characterized by direct connections to sinking streams, high groundwater velocities, and proximity to public water supplies. Objectives of the study were to estimate approximate conduit locations and assess possible anthropogenic influences associated with conduits. Anomalies of young groundwater based on chlorofluorocarbons (CFCs), tritium, and electrical conductivity (EC) indicated fast moving, focused flow and thus the likely presence of conduits. δ¹⁸O was useful for determining sources of recharge for each conduit, and nitrate was a useful tracer for assessing flow paths for anthropogenic influences. Two of the four conduits terminate at or near a large spring complex. CFC apparent ages ranged from 15 years near conduits to >50 years in other areas. Nitrate-N concentrations >0.4 mg/L in groundwater were associated with each of the four conduits compared with concentrations ranging from <0.1 to 0.4 mg/L in other areas. These higher nitrate-N concentrations probably do not result from sinking streams but rather from other areas of infiltration.     

Water balance in the playa-lakes of an arid environment,Monegros, NE Spain

The playa-lakes of the Monegros desert in north-east Spain are saline wetlands in an arid environment, a rare phenomenon in Europe. These extremely valuable habitats are threatened by changes associated with agricultural expansion and incorporation of new irrigated areas. An understanding of the present hydrologic regime will enable changes to be identified, particularly those brought about by flooding and pollution caused by irrigation surplus. This study sets out to show the results of applying a daily water balance in three selected playa-lakes. The balance was in two parts and consisted of: (1) the average balance for all the endorheic basin using the BALAN_11 program, and (2) the water balance in some playa-lakes, applying discharge flows obtained from the previous balance. The resulting volumes of water were converted to water depths and contrasted with reference volumes taken from field and Landsat images. The model was calibrated by applying various hypotheses of function which enabled the results to be adjusted. The proposed balance is an acceptable reproduction of field water measurements during this period, and underlines the consistency of the conceptual model. The methodology used is appropriate for understanding the playa-lakes function and for monitoring them for conservation purposes.     

In situ estimation of the effective chemical diffusion coefficient of a rock matrix in a fractured aquifer

An in situ method of estimating the effective diffusion coefficient for a chemical constituent that diffuses into the primary porosity of a rock is developed by abruptly changing the concentration of the dissolved constituent in a borehole in contact with the rock matrix and monitoring the time-varying concentration. The experiment was conducted in a borehole completed in mudstone on the campus of the University of the Free State in Bloemfontein, South Africa. Numerous tracer tests were conducted at this site, which left a residual concentration of sodium chloride in boreholes that diffused into the rock matrix over a period of years. Fresh water was introduced into a borehole in contact with the mudstone, and the time-varying increase of chloride was observed by monitoring the electrical conductivity (EC) at various depths in the borehole. Estimates of the effective diffusion coefficient were obtained by interpreting measurements of EC over 34 d. The effective diffusion coefficient at a depth of 36 m was approximately 7.8×10^?6 m²/d, but was sensitive to the assumed matrix porosity. The formation factor and mass flux for the mudstone were also estimated from the experiment.     

Uncertainty as the overlap of alternate conditional distributions

An important task in modern geostatistics is the assessment and quantification of resource and reserve uncertainty. This uncertainty is valuable support information for many management decisions. Uncertainty at specific locations and uncertainty in the global resource is of interest. There are many different methods to build models of uncertainty, including Kriging, Cokriging, and Inverse Distance. Each method leads to different results. A method is proposed to combine local uncertainties predicted by different models to obtain a combined measure of uncertainty that combines good features of each alternative. The new estimator is the overlap of alternate conditional distributions.     

Numerical simulation of migration–accumulation of oil resources

Numerical simulation of oil migration and accumulation is to describe the history of oil migration and accumulation in basin evolution. It is of great value to the evaluation of oil resources and to the determination of the location and amount of oil deposits. This thesis discusses the characteristics of petroleum geology and permeation fluid mechanics. For the three-dimensional problems of Dongying hollow of Shengli Petroleum Oil Field, it puts forward a new model and a kind of modified method of upwind finite difference fractional steps implicit interactive scheme. For the famous hydraulic experiment of secondary migration–accumulation, the numerical simulation test has been done, and both the computational and experimental results are basically identical. For the actual problem of Dongying hollow, the numerical simulation test and the actual conditions are basically coincident. Thus, the well-known problem has been solved.     

A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow

We consider conjunctive surface-subsurface flow modeling, where surface water flow is described by the shallow water equations and ground water flow by Richards’ equation for the vadose zone. Coupling between the models is based on the continuity of flux and water pressure. Numerical approximation of the coupled model using the framework of discontinuous Galerkin (DG) methods is formulated. In the subsurface, the local discontinuous Galerkin (LDG) method is used to approximate ground water velocity and hydraulic head; a DG method is also used to approximate surface water velocity and elevation. This approach allows for a weak coupling of the models and the use of different approximating spaces and/or meshes within each regime. A simplified LDG method based on continuous approximations to water head is also described. Numerical results that investigate physical and numerical aspects of surface–subsurface flow modeling are presented. This work was supported by National Science Foundation grant DMS-0411413.     

A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.     

Groundwater resources in the Jabal Al Hass region,northwest Syria: an assessment of past use and future potential

In many cases, the development of groundwater resources to boost agricultural production in dry areas has led to a continuous decline in groundwater levels; this has called into question the sustainability of such exploitation. In developing countries, limited budgets and scarce hydrological data often do not allow groundwater resources to be assessed through groundwater modeling. A case study is presented of a low-cost water-balance approach to groundwater resource assessments in a 1,550 km² semi-arid region in northwestern Syria. The past development of irrigated agriculture and its effect on the groundwater system were studied by analysis of Landsat images and long-term groundwater level changes, respectively. All components of the groundwater balance were determined. Groundwater recharge was estimated using the chloride mass balance method. Over the past three decades, groundwater levels have declined, on average, 23 m, coinciding with a two-fold increase in the groundwater-irrigated area. Groundwater resources are currently depleted by a value that lies between 9.5×10⁶ and 118×10⁶ m³ year^?1, which is larger than can be compensated for by a future decrease in natural discharge or changes in boundary conditions. However, groundwater resources are likely to be sufficient to supply domestic and livestock needs in the area.     

Fractional crystallization models and B–Be–Li systematics at Mt Somma-Vesuvius volcano (Southern Italy)

A compilation of B–Be–Li data on rocks that cover the entire eruptive history of Somma-Vesuvius is presented and interpreted in the light of evolution models for the Somma-Vesuvius rocks. Using major and trace element data, fractional crystalllization models are presented for different geochemical units. These data were used to constrain the source mineralogy of the Somma-Vesuvius rocks (ol-opx-cpx-gar-amp of 0.4-0.3-0.1-0.1-0.1), the amount of sediment added (5–10%) and the melt fraction from batch partial melting computations (0.05–0.1). From the B–Li data it is inferred that the main process responsible for the B isotopic signature is sediment recycling. However, the B–Li data show a major variation in Li abundances respect to B which is explained with Li dehydration before the fluid enriched the mantle wedge that produced the arc magmas. The Somma-Vesuvius B isotope composition is intermediate between that of the Campi Flegrei and the broad field of the Eolian Island arc. A low Be isotopes in the recent volcanic rocks can be explained as: (a) the top 1–22 m of the incoming sediment is accreted, (b) large amounts of sediment erosion, (c) a slow rate of subduction which have provoked a long magmatic history for the Vesuvius magma, (d) the sediment component takes several Myr longer than the subducting plate to reach the magma source region beneath Italy.