NONLINEAR SOLUTIONS OF THE BQUSSINESQ EQUATION AND COMPARISONS WITH FIELD OBSERVATIONS1

ABSTRACT: Two models based on variations of the one dimensional Boussinesq equation are presented. The models are formulated by treating each of the two nonlinearities appearing in the equation separately, and the resulting equations are solved by implicit finite difference techniques. The models have been used to predict drawdowns created from large scale pumping at an open coal mine pit dewatering project in Poland. The predictions were found to be in good agreement with data collected over a period of seven years.     

OPTIMAL GROUND WATER MANAGEMENT IN TWO-AQUIFER SYSTEMS1

ABSTRACT; This paper presents a numerical model for the prediction of optimal ground water withdrawal from a two-aquifer system by observing a set of constraints determined by the ecological conditions of the ground water basin. The aquifer system consists of an upper unconfined and a lower confined aquifer with a leaky stratum between them. It is assumed that water is withdrawn from the confined aquifer only, but the unconfined aquifer will also be affected due to the leakiness of the layer separating the upper and lower aquifers. Simulation and linear programming are employed for developing a computer model for the optimal management of such systems, with the objectives of determining withdrawal rates for predetermined ground water levels.     

INFLUENCES OF INCREASED SAND DELIVERY ON THE MORPHOLOGY OF SAND AND GRAVEL CHANNELS1

ABSTRACT A flume study was conducted to examine (1)changes in the particle-size distribution of sediments in riffles due to the proportion of sand in transport and the total rate of bedload transport at the time the riffle is deposited and (2) the effect of high sand transport rates on the stability of gravel riffles. The median particle size of sediment deposited in the riffle was larger than that of the sediment in transport. Small but significant (a = 0.05) decreases in the median particle size of riffle sediments resulted as the sand-to-gravel ratio. Increased concentrations of sand in transport caused previously stable gravel riffles to undergo scour. These results, in combination with information from other studies, suggest that an alluvial channel with pool-riffle sequences and with sand and gravel beds may respond to an increased delivery of sand by reducing form roughness. Form roughness can be reduced by degrading riffles and filling pools. Subsequent responses may be increases in width-to-depth ratio and slope.     

Development of a Model for Simulation of Solute Transport in a Stream–Aquifer System

A model for simulation of solute transport in a dynamic stream–aquifer system is developed by integrating four existing sub-models. Interface packages are created to link the sub-models. The developed model is successfully used to simulate chloride transport in the stream–aquifer system of the Arkansas River and the Equus Beds Aquifer in Kansas and demonstrates that chloride concentration in the aquifer decreases in the vicinity of the simulated channel over time.     

Surfactant dissolution and mobilization of LNAPL contaminants in aquifers

Improper disposal, accidental spills and leaks of non-aqueousphase liquids (NAPL) such as gasoline, fuel oil and creosote result in long-term persistent sources of groundwater pollution.Column and 2-D tanks experiments were conducted to evaluate the use of surfactant-enhanced recovery of light non-aqueous phase liquids (LNAPL) in groundwater aquifers. These experiments focused on the use of surfactants to promote dissolution and mobilization in addition to evaluating the increase of aqueous phase permeability as residual NAPL is recovered. Further experiments are presented that show the innovative use of surfactants during primary pumping to recover free product canpotentially increase the amount of free product recovered, canpotentially reduce the amount of residual NAPL remaining afterprimary pumping and performs better than the use of surfactantsto mobilize trapped residual NAPL.     

USING GEOLOGY TO IMPROVE FLOOD HAZARD MANAGEMENT ON ALLUVIAL FANS ‐ AN EXAMPLE FROM LAUGHLIN,NEVADA1

A study of the piedmont of the Newberry Mountains near Laughlin, Nevada, demonstrates that geologic information can improve the scientific basis of flood‐hazard management on alluvial fans in desert areas. Comparison of geologic information against flood insurance rate maps (FIRMs) reveals flaws in conventional methods for flood hazard delineation in this setting. Geologic evidence indicates that large parts of the Newberry piedmont have been isolated from significant flooding for at least the past 10,000 years. This contrasts with existing FIRMs that include large tracts of nonflood prone land in the 100‐year and 500‐year flood hazard zones and exclude areas of indisputably flood prone land from the regulatory flood plain. From the basis of the geology, flood hazards on at least one‐third of the piedmont are mischaracterized on the regulatory maps. The formal incorporation of geologic data into flood hazard studies on desert piedmonts could significantly reduce this type of discrepancy and substantially reduce the scope, hence cost, of more elaborate engineering studies and hazard mitigation strategies. The results of this study affirm the value of new Federal Emergency Management Agency (FEMA) recommendations for characterizing alluvial fan flood hazards and support an argument for mandating geological studies in the regulatory process.     

Application of Sample Selection Model in Estimating Response Functions for Nitrate Percolation

The Seymour aquifer region of Texas has been identified as containing elevated levels of nitrate in ground water. Various state and federal agencies are currently studying policy options for the region by gathering more site-specific information. However, because of lack of sufficient information, cause and effect relationships between water quality and agricultural practices have not been well established for the region. Some recently available biophysical simulation models have impressive capabilities in generating large amounts of data on environmental pollution resulting from agricultural production practices. In this study, the data generated by a biophysical simulation model were used to estimate the nitrate percolation response functions for the Seymour aquifer region. Interestingly, nitrate percolation values obtained from simulation models often comprise acensoredsample because the non-zero percolation values are only observed under certain climatic events and input levels. It has been shown in the econometric literature that the use of Ordinary Least Squares (OLS) on censored sample data produces biased and inconsistent parameter estimates. Thus, a sample selection model was used in this study to estimate the response functions for nitrate percolation. The study provides some insight into the relationship between nitrate percolation and agricultural production practices. In particular, the study demonstrates the potential of selected design standards in minimizing agricultural nonpoint-source (NPS) pollution for the study area.     

深厚黄土覆盖层上土石坝地震响应特性分析

采用非线性有限单元法和笔者曾提出的动孔压试验曲线法,对某深厚黄土覆盖层上土石坝进行了有效应力法地震响应分析,重点分析大坝的绝对加速度、动位移和动应力等动力响应及动孔隙水压力分布情况。分析结果表明,在现有设计条件下,由于坝基软弱黄土覆盖层较厚,大坝在7度地震作用下地震响应不强烈,但坝基黄土覆盖层会出现液化情况,需采取相应的抗液化工程措施。     

Simulating Nitrate Leaching Profiles in a Highly Permeable Vadose Zone

An approach is developed to simulate leaching of a dissolved chemical constituent in the vadose zone of an aquifer. Specifically, nitrate loading at the water table for different water table depths, for a range of aquifer permeability values, and for different cases of heterogeneity of the aquifer, are considered. Models from the literature are first used to derive soil–water characteristic curves (water retention and hydraulic conductivity) from a grain size distribution curve for unsaturated conditions. Given infiltration from the surface, the initial conditions for the chemical concentration, and the water content profile, leaching of the chemical in the vadose zone is simulated as a function of both time and depth. The methodology is illustrated for a permeable aquifer. Simulations are undertaken using a finite element code for saturated and unsaturated flow. Different scenarios are simulated depending on the heterogeneity of the aquifer and the depth of the water table. Modeling results show that in the example case studied, nitrate concentration loading at the water table does not depend strongly on the position of the water table, but rather on the material properties of the aquifer. The contribution of this endeavor resides in the methodology which allows a prediction of nitrate leaching using only the grain size property of the aquifer. It allows practitioners to obtain a first assessment of leaching with limited data.     

A method for quick assessment of CO2 storage capacity in closed and semi-closed saline formations

Saline aquifers of high permeability bounded by overlying/underlying seals may be surrounded laterally by low-permeability zones, possibly caused by natural heterogeneity and/or faulting. Carbon dioxide (CO₂) injection into and storage in such “closed” systems with impervious seals, or “semi-closed” systems with non-ideal (low permeability) seals, is different from that in “open” systems, from which the displaced brine can easily escape laterally. In closed or semi-closed systems, the pressure buildup caused by continuous industrial-scale CO₂ injection may have a limiting effect on CO₂ storage capacity, because geomechanical damage caused by overpressure needs to be avoided. In this research, a simple analytical method was developed for the quick assessment of the CO₂ storage capacity in such closed and semi-closed systems. This quick-assessment method is based on the fact that native brine (of an equivalent volume) displaced by the cumulative injected CO₂ occupies additional pore volume within the storage formation and the seals, provided by pore and brine compressibility in response to pressure buildup. With non-ideal seals, brine may also leak through the seals into overlying/underlying formations. The quick-assessment method calculates these brine displacement contributions in response to an estimated average pressure buildup in the storage reservoir. The CO₂ storage capacity and the transient domain-averaged pressure buildup estimated through the quick-assessment method were compared with the “true” values obtained using detailed numerical simulations of CO₂ and brine transport in a two-dimensional radial system. The good agreement indicates that the proposed method can produce reasonable approximations for storage–formation–seal systems of various geometric and hydrogeological properties.