Search

The relationship between aquifer features and short term spring responses to recharge events of both spring discharge and physiochemical parameters of the discharged water was studied using a process-based flow and transport mode and in the respective software, a pipe network model was coupled to MODFLOW. The application of the new modeling tool demonstrated that variations of pyhsicochemical parameters depend to a large extent on the intensity and duration of recharge events.

A tracer test in a carbonate aquifer is analyzed using the method of moments and two analytical advection-dispersion models (ADMs) as well as a numerical model. The numerical model is a coupled continuum-pipe flow and transport model that accounts for two different flow components in karstified carbonate aquifers, i.e., rapid and often turbulent conduit flow and Darcian flow in the fissured porous rock. All techniques employed provide reasonable fits to the tracer breakthrough curve (TBC) measured at a spring. The resulting parameter estimates are compared to investigate how each conceptual model of flow and transport processes that forms the basis of the analyses affects the interpretation of the tracer test. Numerical modeling results suggest that the method of moments and the analytical ADMs tend to overestimate the conduit volume because part of the water discharged at the spring is wrongly attributed to the conduit system if flow in the fissured porous rock is ignored. In addition, numerical modeling suggests that mixing of the two flow components accounts for part of the dispersion apparent in the measured TBC, while the remaining part can be attributed to Taylor dispersion. These processes, however, cannot reasonably explain the tail of the TBC. Instead, retention in immobile-fluid regions as included in a nonequilibrium ADM provides a possible explanation., Introduction Ground water flow in carbonate aquifers is often focused in solution conduits that comprise only a small percentage of total aquifer porosity. The majority of aquifer storage, however, is [...]

Spatial variations in hydraulic conductivity (K) can significantly affect the transport of contaminants in ground water. Conventional field methods, however, rarely provide a description of these variations at the level of detail necessary for reliable transport predictions and effective remediation designs. A direct-push (DP) method, hydrostratigraphic profiling, has been developed to characterize the spatial variability of both electrical conductivity (EC) and hydraulic conductivity in unconsolidated formations in a cost-effective manner. This method couples a dual-rod approach for performing slug tests in DP equipment with high-resolution EC logging. The method was evaluated at an extensively studied site in the Kansas River floodplain. A series of profiles was performed on a surface grid, resulting in a detailed depiction of the three-dimensional distribution of EC and K. Good agreement was found between K estimates obtained from this approach and those obtained using other methods. The results of the field evaluation indicate that DP hydrostratigraphic profiling is a promising method for obtaining detailed information about spatial variations in subsurface properties without the need for permanent wells., Introduction Spatial variations in hydraulic conductivity (K) can have a significant impact on the transport of contaminants in ground water. An understanding of the three-dimensional distribution of K at a [...]

A new method was developed and implemented into an Excel Visual Basic for Applications ( VBAs) algorithm utilizing trigonometry laws in an innovative way to overlap recession segments of time series and create master recession curves ( MRCs). Based on a trigonometry approach, the algorithm horizontally translates succeeding recession segments of time series, placing their vertex, that is, the highest recorded value of each recession segment, directly onto the appropriate connection line defined by measurement points of a preceding recession segment. The new method and algorithm continues the development of methods and algorithms for the generation of MRC, where the first published method was based on a multiple linear/nonlinear regression model approach (Posavec et al. 2006). The newly developed trigonometry-based method was tested on real case study examples and compared with the previously published multiple linear/nonlinear regression model-based method. The results show that in some cases, that is, for some time series, the trigonometry-based method creates narrower overlaps of the recession segments, resulting in higher coefficients of determination R2, while in other cases the multiple linear/nonlinear regression model-based method remains superior. The Excel VBA algorithm for modeling MRC using the trigonometry approach is implemented into a spreadsheet tool ( MRCTools v3.0 written by and available from Kristijan Posavec, Zagreb, Croatia) containing the previously published VBA algorithms for MRC generation and separation. All algorithms within the MRCTools v3.0 are open access and available free of charge, supporting the idea of running science on available, open, and free of charge software. [ABSTRACT FROM AUTHOR]

Karst aquifers represent dual flow systems consisting of a highly conductive conduit system embedded in a less permeable rock matrix. Hybrid models iteratively coupling both flow systems generally consume much time, especially because of the nonlinearity of turbulent conduit flow. To reduce calculation times compared to those of existing approaches, a new iterative equation solver for the conduit system is developed based on an approximated Newton-Raphson expression and a Gauß-Seidel or successive over-relaxation scheme with a single iteration step at the innermost level. It is implemented and tested in the research code CAVE but should be easily adaptable to similar models such as the Conduit Flow Process for MODFLOW-2005. It substantially reduces the computational effort as demonstrated by steady-state benchmark scenarios as well as by transient karst genesis simulations. Water balance errors are found to be acceptable in most of the test cases. However, the performance and accuracy may deteriorate under unfavorable conditions such as sudden, strong changes of the flow field at some stages of the karst genesis simulations. [ABSTRACT FROM AUTHOR]

A new MODFLOW package (Nonlinear Flow Process; NLFP) simulating nonlinear flow following the Forchheimer equation was developed and implemented in MODLFOW-2005. The method is based on an iterative modification of the conductance calculated and used by MODFLOW to obtain an effective Forchheimer conductance. The package is compatible with the different layer types, boundary conditions, and solvers as well as the wetting capability of MODFLOW. The correct implementation is demonstrated using four different benchmark scenarios for which analytical solutions are available. A scenario considering transient flow in a more realistic setting and a larger model domain with a higher number of cells demonstrates that NLFP performs well under more complex conditions, although it converges moderately slower than the standard MODFLOW depending on the nonlinearity of flow. Thus, this new tool opens a field of opportunities to groundwater flow simulation with MODFLOW, especially for core sample simulation or vuggy karstified aquifers as well as for nonlinear flow in the vicinity of pumping wells. [ABSTRACT FROM AUTHOR]

As a result of rock dissolution processes, karst aquifers exhibit highly conductive features such as caves and conduits. Within these structures, groundwater flow can become turbulent and therefore be described by nonlinear gradient functions. Some numerical groundwater flow models explicitly account for pipe hydraulics by coupling the continuum model with a pipe network that represents the conduit system. In contrast, the Conduit Flow Process Mode 2 (CFPM2) for MODFLOW-2005 approximates turbulent flow by reducing the hydraulic conductivity within the existing linear head gradient of the MODFLOW continuum model. This approach reduces the practical as well as numerical efforts for simulating turbulence. The original formulation was for large pore aquifers where the onset of turbulence is at low Reynolds numbers (1 to 100) and not for conduits or pipes. In addition, the existing code requires multiple time steps for convergence due to iterative adjustment of the hydraulic conductivity. Modifications to the existing CFPM2 were made by implementing a generalized power function with a user-defined exponent. This allows for matching turbulence in porous media or pipes and eliminates the time steps required for iterative adjustment of hydraulic conductivity. The modified CFPM2 successfully replicated simple benchmark test problems. [ABSTRACT FROM AUTHOR]

A method to estimate reactive transport parameters as well as geometric conduit parameters from a multitracer test in a karst aquifer is provided. For this purpose, a calibration strategy was developed applying the two-region nonequilibrium model CXTFIT. The ambiguity of the model calibration was reduced by first calibrating the model with respect to conservative tracer breakthrough and later transferring conservative transport parameters to the reactive model calibration. The reactive transport parameters were only allowed to be within a defined sensible range to get reasonable calibration values. This calibration strategy was applied to breakthrough curves obtained from a large-scale multitracer test, which was performed in a karst aquifer of the Swabian Alb, Germany. The multitracer test was conducted by the simultaneous injection of uranine, sulforhodamine G, and tinopal CBS-X. The model succeeds to represent the tracer breakthrough curves (TBCs) of uranine and sulforhodamine G and verifies that tracer-rock interactions preferably occur in the immobile fluid region, although the fraction of this region amounts to only 3.5% of the total water. However, the model failed to account for the long tailing observed in the TBC of tinopal CBS-X. Sensitivity analyses reveal that model results for the conservative tracer transport are most sensitive to average velocity and volume fraction of the mobile fluid region, while dispersion and mass transfer coefficients are least influential. Consequently, reactive tracer calibration allows the determination of sorption sites in the mobile and immobile fluid region at small retardation coefficients. [ABSTRACT FROM AUTHOR]

Ground water in karst terrains is highly vulnerable to contamination due to the rapid transport of contaminants through the highly conductive conduit system. For contamination risk assessment purposes, information about hydraulic and geometric characteristics of the conduits and their hydraulic interaction with the fissured porous rock is an important prerequisite. The relationship between aquifer characteristics and short-term responses to recharge events of both spring discharge and physicochemical parameters of the discharged water was examined using a process-based flow and transport model. In the respective software, a pipe-network model, representing fast conduit flow, is coupled to MODFLOW, which simulates flow in the fissured porous rock. This hybrid flow model was extended to include modules simulating heat and reactive solute transport in conduits. The application of this modeling tool demonstrates that variations of physicochemical parameters, such as solute concentration and water temperature, depend to a large extent on the intensity and duration of recharge events and provide information about the structure and geometry of the conduit system as well as about the interaction between conduits and fissured porous rock. Moreover, the responses of solute concentration and temperature of spring discharge appear to reflect different processes, thus complementing each other in the aquifer characterization. [ABSTRACT FROM AUTHOR]