Hydrologic analysis of the High Plains aquifer system in Box Butte County, Nebraska /

Mode of access: Internet.

A Comparative Analysis of Two Slug Test Methods in Puget Lowland Glacio-Fluvial Sediments near Coupeville, WA

Two different slug test field methods are conducted in wells completed in a Puget Lowland aquifer and are examined for systematic error resulting from water column displacement techniques. Slug tests using the standard slug rod and the pneumatic method were repeated on the same wells and hydraulic conductivity estimates were calculated according to Bouwer & Rice and Hvorslev before using a non-parametric statistical test for analysis. Practical considerations of performing the tests in real life settings are also considered in the method comparison. Statistical analysis indicates that the slug rod method results in up to 90% larger hydraulic conductivity values than the pneumatic method, with at least a 95% certainty that the error is method related. This confirms the existence of a slug-rod bias in a real world scenario which has previously been demonstrated by others in synthetic aquifers. In addition to more accurate values, the pneumatic method requires less field labor, less decontamination, and provides the ability to control the magnitudes of the initial displacement, making it the superior slug test procedure.

Numerical investigation of seawater intrusion at Gooburrum, Bundaberg, Queensland, Australia

Seawater intrusion in coastal agricultural areas due to groundwater abstraction is a major environmental problem along the northeastern coast of Australia. Management options are being explored using numerical modelling, however, questions remain concerning the appropriate level of sophistication in models, choice of seaward boundary conditions, and how to accommodate heterogeneity and data uncertainty. The choice of seaward boundary condition is important since it affects the amount of salt transported into the aquifers and forms the focus of the present study. The impact of this boundary condition is illustrated for the seawater-intrusion problem in the Gooburrum aquifers, which occur within Tertiary sedimentary strata. A two-dimensional variable-density groundwater and solute-transport model was constructed using the computer code 2DFEMFAT (Cheng et al. 1998). The code was tested against an experiment for a steady-state freshwater-saltwater interface and against the Elder (Elder 1967) free-convection problem. Numerical simulations show that the imposition of the commonly-used equivalent hydrostatic freshwater heads, combined with a constant salt concentration at the seaward boundary, results in overestimated seawater intrusion in the lower Gooburrum aquifer. Since the imposition of this boundary condition allows water flow across the boundary, which subsequently takes salt into the aquifer, a careful check is essential to estimate whether too much mass of salt is introduced.

Two-dimensional approximation for tide-induced watertable fluctuations in a sloping sandy beach

The prediction of watertable fluctuations in a coastal aquifer is important for coastal management. However, most previous approaches have based on the one-dimensional Boussinesq equation, neglecting variations in the coastline and beach slope. In this paper, a closed-form analytical solution for a two-dimensional unconfined coastal aquifer bounded by a rhythmic coastline is derived. In the new model, the effect of beach slope is also included, a feature that has not been considered in previous two-dimensional approximations. Three small parameters, the shallow water parameter (epsilon), the amplitude parameter (a) and coastline parameter (beta) are used in the perturbation approximation. The numerical results demonstrate the significant influence of both the coastline shape and beach slopes on tide-driven coastal groundwater fluctuations. (c) 2004 Elsevier Ltd. All rights reserved.

Similarity solution of axisymmetric flow in porous media

Applications of the axisymmetric Boussinesq equation to groundwater hydrology and reservoir engineering have long been recognised. An archetypal example is invasion by drilling fluid into a permeable bed where there is initially no such fluid present, a circumstance of some importance in the oil industry. It is well known that the governing Boussinesq model can be reduced to a nonlinear ordinary differential equation using a similarity variable, a transformation that is valid for a certain time-dependent flux at the origin. Here, a new analytical approximation is obtained for this case. The new solution,, which has a simple form, is demonstrated to be highly accurate. (c) 2005 Elsevier Ltd. All rights reserved.

Application of a coupled ground-surface water flow model to simulate periodic groundwater flow influenced by a sloping boundary, capillarity and vertical flows

Theoretical developments as well as field and laboratory data have shown the influence of the capillary fringe on water table fluctuations to increase with the fluctuation frequency. The numerical solution of a full, partially saturated flow equation can be computationally expensive. In this paper, the influence of the capillary fringe on water table fluctuations is simplified through its parameterisation into the storage coefficient of a fully-saturated groundwater flow model using the complex effective porosity concept [Nielsen, P., Perrochet, P., 2000. Water table dynamics under capillary fringes: experiments and modelling. Advances in Water Resources 23 (1), 503-515; Nielsen, P., Perrochet, P., 2000. ERRATA: water table dynamics under capillary fringes: experiments and modelling (Advances in Water Resources 23 (2000) 503-515). Advances in Water Resources 23, 907-908]. The model is applied to sand flume observations of periodic water table fluctuations induced by simple harmonic forcing across a sloping boundary, analogous to many beach groundwater systems. While not providing information on the moisture distribution within the aquifer, this approach can reasonably predict the water table fluctuations in response to periodic forcing across a sloping boundary. Furthermore, he coupled ground-surface water model accurately predicts the extent of the seepage face formed at the sloping boundary. (C) 2005 Elsevier Ltd. All rights reserved.

Regional-scale, fully coupled modelling of stream-aquifer interaction in a tropical catchment

The planning and management of water resources in the Pioneer Valley, north-eastern Australia requires a tool for assessing the impact of groundwater and stream abstractions on water supply reliabilities and environmental flows in Sandy Creek (the main surface water system studied). Consequently, a fully coupled stream-aquifer model has been constructed using the code MODHMS, calibrated to near-stream observations of watertable behaviour and multiple components of gauged stream flow. This model has been tested using other methods of estimation, including stream depletion analysis and radon isotope tracer sampling. The coarseness of spatial discretisation, which is required for practical reasons of computational efficiency, limits the model's capacity to simulate small-scale processes (e.g., near-stream groundwater pumping, bank storage effects), and alternative approaches are required to complement the model's range of applicability. Model predictions of groundwater influx to Sandy Creek are compared with baseflow estimates from three different hydrograph separation techniques, which were found to be unable to reflect the dynamics of Sandy Creek stream-aquifer interactions. The model was also used to infer changes in the water balance of the system caused by historical land use change. This led to constraints on the recharge distribution which can be implemented to improve model calibration performance. (c) 2006 Elsevier B.V. All rights reserved.

Characterisation of sea-water intrusion in the Pioneer Valley, Australia using hydrochemistry and three-dimensional numerical modelling

Sea-water intrusion is actively contaminating fresh groundwater reserves in the coastal aquifers of the Pioneer Valley,north-eastern Australia. A three-dimensional sea-water intrusion model has been developed using the MODHMS code to explore regional-scale processes and to aid assessment of management strategies for the system. A sea-water intrusion potential map, produced through analyses of the hydrochemistry, hydrology and hydrogeology, offsets model limitations by providing an alternative appraisal of susceptibility. Sea-water intrusion in the Pioneer Valley is not in equilibrium, and a potential exists for further landward shifts in the extent of saline groundwater. The model required consideration of tidal over-height (the additional hydraulic head at the coast produced by the action of tides), with over-height values in the range 0.5-0.9 m giving improved water-table predictions. The effect of the initial water-table condition dominated the sensitivity of the model to changes in the coastal hydraulic boundary condition. Several salination processes are probably occurring in the Pioneer Valley, rather than just simple landward sea-water advancement from modern sources of marine salts. The method of vertical discretisation (i.e. model-layer subdivision) was shown to introduce some errors in the prediction of watertable behaviour.

Swash-aquifer interaction in the vicinity of the water table exit point on a sandy beach

The coupling of sandy beach aquifers with the swash zone in the vicinity of the water table exit point is investigated through simultaneous measurements of the instantaneous shoreline (swash front) location, pore pressures and the water table exit point. The field observations reveal new insights into swash-aquifer coupling not previously gleaned from measurements of pore pressure only. In particular, for the case where the exit point is seaward of the observation point, the pore pressure response is correlated with the distance between the exit point and the shoreline in that when the distance is large the rate of pressure drop is fast and when the distance is small the rate decreases. The observations expose limitations in a simple model describing exit point dynamics which is based only on the force balance on a particle of water at the sand surface and neglects subsurface pressures. A new modified form of the model is shown to significantly improve the model-data comparison through a parameterization of the effects of capillarity into the aquifer storage coefficient. The model enables sufficiently accurate predictions of the exit point to determine when the swash uprush propagates over a saturated or a partially saturated sand surface, potentially an important factor in the morphological evolution of the beach face. Observations of the shoreward propagation of the swash-induced pore pressure waves ahead of the runup limit shows that the magnitude of the pressure fluctuation decays exponentially and that there is a linear increase in time lags, behavior similar to that of tidally induced water table waves. The location of the exit point and the intermittency of wave runup events is also shown to be significant in terms of the shore-normal energy distribution. Seaward of the mean exit point location, peak energies are small because of the saturated sand surface within the seepage face acting as a "rigid lid'' and limiting pressure fluctuations. Landward of the mean exit point the peak energies grow before decreasing landward of the maximum shoreline position.

Tidal impacts on riparian salinities near estuaries

Groundwater-dependent riparian biota is known to be sensitive to changes in soil and groundwater salinity in estuarine systems. The groundwater flow and salinity behaviour in a phreatic aquifer adjoining a partially penetrating, tidal. estuary is investigated through two-dimensional numerical experiments for a lateral cross-section, which explore the influence of factors, such as aquifer and soil materials, tidal amptitudes, and regional groundwater hydraulic gradients. The density contrast between estuarine water and the fresh groundwater drives saltwater penetration of the aquifer even in the case of a marked groundwater hydraulic gradient towards the estuary. We show that tidal fluctuations in estuaries can significantly affect the groundwater salinity distribution in adjacent density-stratified phreatic aquifers. This has consequences for the expected distribution of salinity-sensitive biota in the hyporheic zone as well as vegetation and fauna dependent on water in the riparian soil and aquifer. The shape of the dense saltwater wedge propagating into the adjacent groundwater system is also modified by the estuarine tidal signal, although this effect appears to have only minor influence on the maximum distance penetrated into the aquifer (i.e., location of the 'toe' of the wedge). Tide-induced changes to riparian groundwater salinity are advection-driven, as evidenced by the modified time-averaged groundwater flow dynamics. (c) 2006 Elsevier B.V. All rights reserved.

Zinc deposits and related mineralization of the Burketown mineral field, including the world-class Century Deposit, Northern Australia: Fluid inclusion and stable isotope evidence for basin fluid sources

The stratiform Century Zn-Pb deposit and the discordant Zn-Pb lode deposits of the Burketown mineral field, northern Australia, host ore and gangue minerals with primary fluid inclusions that have not been affected by the Isan orogeny, thus providing a unique opportunity to investigate the nature of the ore-forming brines. All of the deposits are hosted in shales and siltstones belonging to the Isa superbasin and comprise sphalerite, pyrite, carbonate, quartz, galena, minor chalcopyrite, and minor illite. According to Pb model ages, the main ore stage of mineralization at Century formed at I575 Ma, some 20 m.y. after deposition of the host shale sequence. Microthermometry on undeformed, primary fluid inclusions hosted in porous sphalerite shows that the Zn at Century was transported to the deposit by a homogeneous, Ca2+- and Na+-bearing brine with a salinity of 21.6 wt percent NaCl equiv. delta D-fluid of the fluid inclusion water ranges from -89 to -83 per mil, consistent with a basinal brine that evolved from meteoric water. Fluid inclusion homogenization temperatures range between 74 degrees and 125 degrees C, which are lower than the 120 degrees to 160 degrees C range calculated from vitrinite reflectance and illite crystallinity data from the deposit. This discrepancy indicates that mineralization likely formed at 50 to 85 Mpa, corresponding to a depth of 1,900 to 3,100 m. Transgressive galena-sphalerite veins that cut stratiform mineralization at Century and breccia-filled quartz-dolomite-sphalerite-galena veins in the discordant Zn-Pb lodes have Pb model ages between 1575 and 1485 Ma. Raman spectroscopy and microthermometry reveal that the primary fluid inclusions in these veins contain Ca2+, Na+. but they have lower salinities between 23 and 10 wt percent NaCl equiv and higher delta D-fluid values ranging from -89 to -61 per mil than fluid inclusions in porous sphalerite from Century. Fluid inclusion water from sphalerite in one of the lode deposits has delta O-18(fluid) values of 1.6 and 2.4 per mil, indistinguishable from delta O-18(fluid) values between -0.3 to +7.4 per mil calculated from the isotopic composition of coexisting quartz, dolomite, and illite. The trend toward lower salinities and higher delta D-fluid values relative to the earlier mineralizing fluids is attributed to mixing between the fluid that formed Century and a seawater-derived fluid from a different source. Based on seismic data from the Lawn Hill platform and paragenetic and geochemical results from the Leichhardt River fault trough to the south, diagenetic aquifers in the Underlying Calvert superbasin appear to have been the most likely sources for the fluids that formed Century and the discordant lode deposits. Paragenetically late sphalerite and calcite cut sphalerite, quartz, and dolomite in the lode deposits and contain Na+-dominated fluid inclusions with much lower salinities than their older counterparts. The isotopic composition of calcite also indicates delta O-18(fluid) from 3.3 to 10.7 per mil, which is larger than the range obtained from synmineralization minerals, supporting the idea that a unique fluid source was involved. The absolute timing of this event is unclear, but a plethora of Pb model, K-Ar, and Ar-40/Ar-39 ages between 1440 and 1300 Ma indicate that a significant volume of fluid was mobilized at this time. The deposition of the Roper superbasin from ca. 1492 +/- 4 Ma suggests that these late veins formed from fluids that may have been derived from aquifers in overlying sediments of the Roper superbasin. Clear, buck, and drusy quartz in veins unrelated to any form of Pb-Zn mineralization record the last major fluid event in the Burketown mineral field and form distinct outcrops and ridges in the district. Fluid inclusions in these veins indicate formation from a low-salinity, 300 degrees +/- 80 degrees C fluid. Temperatures approaching 300 degrees C recorded in organic matter adjacent to faults and at sequence boundaries correspond to K-Ar ages spanning 1300 to 1100 Ma, which coincides with regional hydrothermal activity in the northern Lawn Hill platform and the emplacement of the Lakeview Dolerite at the time of assemblage of the Rodinia supercontinent.

Stochastic modelling and simulations for solute transport in porous media

A stochastic model for solute transport in aquifers is studied based on the concepts of stochastic velocity and stochastic diffusivity. By applying finite difference techniques to the spatial variables of the stochastic governing equation, a system of stiff stochastic ordinary differential equations is obtained. Both the semi-implicit Euler method and the balanced implicit method are used for solving this stochastic system. Based on the Karhunen-Loeve expansion, stochastic processes in time and space are calculated by means of a spatial correlation matrix. Four types of spatial correlation matrices are presented based on the hydraulic properties of physical parameters. Simulations with two types of correlation matrices are presented.

Indagine modellistica sulla tenuta idraulica delle difese arginali deputate alla protezione idraulica del territorio

La prima parte di questo lavoro di tesi tratta dell’interazione tra un bacino di laminazione e il sottostante acquifero: è in fase di progetto, infatti, la costruzione di una cassa di espansione sul torrente Baganza, a monte della città di Parma. L’obiettivo di tale intervento è di ridurre il rischio di esondazione immagazzinando temporaneamente, in un serbatoio artificiale, la parte più pericolosa del volume di piena che verrebbe rilasciata successivamente con portate che possono essere agevolmente contenute nel tratto cittadino del torrente. L’acquifero è stato preliminarmente indagato e monitorato permettendone la caratterizzazione litostratigrafica. La stratigrafia si può riassumere in una sequenza di strati ghiaioso-sabbiosi con successione di lenti d’argilla più o meno spesse e continue, distinguendo due acquiferi differenti (uno freatico ed uno confinato). Nel presente studio si fa riferimento al solo acquifero superficiale che è stato modellato numericamente, alle differenze finite, per mezzo del software MODFLOW_2005. L'obiettivo del presente lavoro è di rappresentare il sistema acquifero nelle condizioni attuali (in assenza di alcuna opera) e di progetto. La calibrazione è stata condotta in condizioni stazionarie utilizzando i livelli piezometrici raccolti nei punti d’osservazione durante la primavera del 2013. I valori di conducibilità idraulica sono stati stimati per mezzo di un approccio geostatistico Bayesiano. Il codice utilizzato per la stima è il bgaPEST, un software gratuito per la soluzione di problemi inversi fortemente parametrizzati, sviluppato sulla base dei protocolli del software PEST. La metodologia inversa stima il campo di conducibilità idraulica combinando osservazioni sullo stato del sistema (livelli piezometrici nel caso in esame) e informazioni a-priori sulla struttura dei parametri incogniti. La procedura inversa richiede il calcolo della sensitività di ciascuna osservazione a ciascuno dei parametri stimati; questa è stata valutata in maniera efficiente facendo ricorso ad una formulazione agli stati aggiunti del codice in avanti MODFLOW_2005_Adjoint. I risultati della metodologia sono coerenti con la natura alluvionale dell'acquifero indagato e con le informazioni raccolte nei punti di osservazione. Il modello calibrato può quindi essere utilizzato come supporto alla progettazione e gestione dell’opera di laminazione. La seconda parte di questa tesi tratta l'analisi delle sollecitazioni indotte dai percorsi di flusso preferenziali causati da fenomeni di piping all’interno dei rilevati arginali. Tali percorsi preferenziali possono essere dovuti alla presenza di gallerie scavate da animali selvatici. Questo studio è stato ispirato dal crollo del rilevato arginale del Fiume Secchia (Modena), che si è verificato in gennaio 2014 a seguito di un evento alluvionale, durante il quale il livello dell'acqua non ha mai raggiunto la sommità arginale. La commissione scientifica, la cui relazione finale fornisce i dati utilizzati per questo studio, ha attribuito, con molta probabilità, il crollo del rilevato alla presenza di tane di animali. Con lo scopo di analizzare il comportamento del rilevato in condizioni integre e in condizioni modificate dall'esistenza di un tunnel che attraversa il manufatto arginale, è stato realizzato un modello numerico 3D dell’argine mediante i noti software Femwater e Feflow. I modelli descrivono le infiltrazioni all'interno del rilevato considerando il terreno in entrambe le porzioni sature ed insature, adottando la tecnica agli elementi finiti. La tana è stata rappresentata da elementi con elevata permeabilità e porosità, i cui valori sono stati modificati al fine di valutare le diverse influenze sui flussi e sui contenuti idrici. Per valutare se le situazioni analizzate presentino o meno il verificarsi del fenomeno di erosione, sono stati calcolati i valori del fattore di sicurezza. Questo è stato valutato in differenti modi, tra cui quello recentemente proposto da Richards e Reddy (2014), che si riferisce al criterio di energia cinetica critica. In ultima analisi è stato utilizzato il modello di Bonelli (2007) per calcolare il tempo di erosione ed il tempo rimanente al collasso del rilevato.

Theory and application of inverse problems in groundwater: numerical, laboratory and field studies.

Aquifers are a vital water resource whose quality characteristics must be safeguarded or, if damaged, restored. The extent and complexity of aquifer contamination is related to characteristics of the porous medium, the influence of boundary conditions, and the biological, chemical and physical processes. After the nineties, the efforts of the scientists have been increased exponentially in order to find an efficient way for estimating the hydraulic parameters of the aquifers, and thus, recover the contaminant source position and its release history. To simplify and understand the influence of these various factors on aquifer phenomena, it is common for researchers to use numerical and controlled experiments. This work presents some of these methods, applying and comparing them on data collected during laboratory, field and numerical tests. The work is structured in four parts which present the results and the conclusions of the specific objectives.