<p>The analysis of geologic, hydrologic and hydrogeologic data interpreted to give the characteristics of the Cape Flats aquifer showed the quality of groundwater from the aquifer is suitable for development as a water resource. The conceptual model of the Cape Flats sand shows an unconfined sandy aquifer, grading into semi-confined conditions in some places where thick lenses of clay and peat exists. Recharge rates through the saturated zone of the Cape Flats aquifer have been determined by water table fluctuation (WTF), rainfall-recharge relationship, soil water balance and chloride mass balance methods (CMB). Recharge rates using the WTF vary considerably between wet and dry years and between locations, with a range of 17.3% to 47.5%. Values obtained from empirical rainfall-recharge equation (method 2) agree with those of the WTF. Recharge estimates from the water balance model are comparatively lower but are within the range calculated using empirical method 2 (i.e. 87 &ndash / 194 mm or 4 &ndash / 21% of MAP). These recharge rates also agree with estimates from the series of other methods applied to sites located in the north-western coast of Western Cape and are comparable to recharge rates obtained elsewhere in the world.</p>

Sediment delivery is a major problem in the Green River, Kentucky, home of 71 of the state’s 103 known mussel species and 151 fish species. The river also provides water for many of its surrounding counties. This research focuses on how suspended sediment loads, grain size, and sediment concentration during runoff events are related to watershed characteristics. The research characterized suspended sediment loads, grain size, and sediment concentration during runoff events and how they were related to watershed characteristics such as hydro-climatic regime, watershed size, geology and soils, topography and landuse conditions and land cover conditions. The study focused on Brush Creek and Pitman Creek watersheds in the Upper Green River Basin. This research can help in the planning and development of effective environmental strategies by screening out mitigation measures that would not be effective for implementation to minimize sediment load and suspended sediment concentration in the Green River, thereby improving the water quality of the river. Water quality was monitored using data sondes positioned at selected sites in the two watersheds. Water samples were collected during turbidity thresholds of 100 NTU and analyzed for suspended sediment concentrations. Regression models between ‘discharge and stage’ and also between ‘average turbidity and suspended sediment concentration’ were formulated and load estimates were made and compared. Four sets of samples were collected, two at Brush Creek on 11 April (Brush Creek’s event 1) and 3 May (Brush Creek’s event 2) and the other two at Pitman Creek on the 12 February (Pitman Creek’s event 1) and 3 March (Pitman Creek’s event 2) all in the year 2008. The suspended sediment samples collected for all four events were well graded but had relatively more silt than clay and sand. This could be due to the fact that more time and energy was needed to break the bonds in clay minerals or particles and also to the fact that more energy was also needed to transport sand compared to silt. Brush Creek watershed’s particles had smaller grain sizes than Pitman Creek watershed’s particles. All four events showed clockwise hysteresis indicating that most of the sediments from both watersheds during the events were derived from the bed and banks of the channel or area adjacent to the channel. The 11 April event (Brush Creek’s event 1) produced an estimated load of 1.1 x 105 kg and a sediment yield of 5.3 x 102 kg/km2. The 3 May event (Brush Creek’s event 2) produced an estimated load of 3.8 x 104 kg and a sediment yield of 1.8 x 102 kg/km2. Brush Creek watershed’s estimated load for the period compared was 4.9 x 105 kg and a sediment yield of 2.3 x 103 kg/km2 (53 kg/km2/day). The 12 February event (Pitman Creek’s event 1) produced an estimated load of 2.9 x 105 kg and a sediment yield of 8.4 x 102 kg/km2. The 3 March event (Pitman Creek’s event 2) produced an estimated load of 5.7 x 105 kg and a sediment yield of 1.6 x 103 kg/km2. Pitman Creek watershed’s estimated load for the period compared was 1.1 x 106 kg and a sediment yield of 3.1 x 103 kg/km2 (71 kg/km2/day). Pitman Creek watershed’s higher number of stream network per unit area, its high elevation and relief, its high percentage of erodible soil per unit area, its lesser area of protection of erodible soil by its vegetation compared to Brush Creek watershed’s are responsible for its higher sediment load and yield.

Groundwater can be a critical issue to be considered in civil engineering, mining engineering and interdisciplinary fields. Karstic structures and aquifers enclosing groundwater are potential risks in case they are not studied in detail. Enclosed groundwater can result in floods at pit bottom or can cause instabilities of permanent pit slopes. This study is about analyses of flooding possibility at the pit bottom and possible instabilities of pit slopes in the Elbistan-&Ccedil / &ouml / llolar open cast coal mine due to the presence of a karstic aquifer under the lignite formation. Thickness and permeability of the bottom clay formation under the lignite bed are necessary critical parameters for investigating a possible water rush from a confined aquifer in limestone formation underneath the bottom clay. These parameters were changed, and water flow quantities towards the pit bottom were determined by finite element models. Critical values of these parameters were investigated considering the lack of accurate site investigation information regarding the thickness and permeability of bottom clay. Possible strength loss, fracturing, and thus permeability increase in bottom clay due to a confined aquifer were studied. In flooding and slope stability analyses Phase2 software based on finite element method is used. Results of analyses showed that as reported thickness of bottom clay is around 120 m at the pit bottom and permeability values are in orders of magnitudes of 10-8 m/s, no serious flooding problems are expected to occur unless the thickness of bottom clay layer drops down to around 20 m, and the permeability of this layer reaches an order of magnitude of 10-5 m/s. Mechanical effects of confined aquifer on slopes and bottom clay displacements were investigated, and thus fracturing and failure possibilities of bottom clay and permanent slope were assessed. Slope and pit bottom displacements increased to meter levels for less than 60 m bottom clay thicknesses. Whereas 50-60 m bottom clay thickness can be critical for cracking, 20 m bottom clay thickness was found to be critical for water rush to the pit bottom. With reported bottom clay thickness of 120 m and with 25o slope angle permanent slope factor of safety was found to be 1.2, and this value was not effected unless clay layer thickness drops below 70 m levels. Higher than 32o overall slope angle there will be a risk of slope failure for permanent and production slopes, reflected by safety factors less than one, in the stability analyses.

Density-dependent groundwater flow in the Floridan aquifer system (FAS) depends on chloride concentrations and fluid temperature. Previous studies addressing the role of chloride concentration and temperatures on groundwater flow in the FAS have relied on observation data or simplified two-dimensional numerical models. A three-dimensional hydrologic analysis of FAS in peninsular Florida was performed using a modified version of SUTRA (SUTRA-MS) capable of simulating multi-species solute and heat transport. SUTRA-MS was developed during this investigation and is capable of reproducing results for several problems with known solutions.The model was developed using available geometric and hydraulic parameter data and calibrated using hydraulic head, chloride concentrations, and temperatures representative of conditions prior to significant groundwater pumpage from the FAS. The calibrated model is capable of reproducing observed pressures and temperatures but in general ov er-simulates chloride concentrations. The inability of the model to simulate observed chloride concentrations suggests chloride concentrations in the FAS are not in equilibrium with current sea level. Previous hydrologic studies of the FAS have attributed anomalous chloride concentrations to incomplete flushing of relict seawater that entered the aquifer during previous sea-level highstands.Three hypothetical, sinusoidal sea-level changes occurring over 100,000-years were used to evaluate how the aquifer responds to sea-level fluctuations. Model results indicate pressure equilibrates most rapidly and is followed by temperatures and then chloride concentrations. Confining unit thicknesses directly affect response times of pressure, temperature, and chloride concentrations in the FAS.Simulation of the system with ("geothermal case") and without ("isothermal case") the geothermal component reveals that the inflow of seawater from the Florida Straits would be similar without the heat f low but the distribution would differ significantly. The addition of heat flow also reduces the asymmetry of the circulation. Simulations evaluating aquifer responses to sea-level fluctuations and the thermal component indicate that the complicated three-dimensional setting of the FAS is a key component of the groundwater flow system and steady state conditions may not exist for relatively thick coastal aquifers that have experienced multiple sea-level cycles.

Knowledge of fracture orientation and density is significant for reservoir and aquifer characterization. In this study, field experiments are designed to estimate fracture parameters in situ from seismic and GPR (radar) data. The seismic experiment estimates parameters of orientation, density, and filling material. The GPR experiment estimates channel flow geometry and aperture. In the seismic study, lines of 2D data are acquired in a vertically fractured limestone at three different azimuths to look for differences in seismic velocities. A sledgehammer, vertical source and a multicomponent, Vibroseis source are used with multicomponent receivers. Acquisition parameters of frequency, receiver spacing and source-to-receiver offset are varied. The entire suite of seismic body waves and Rayleigh waves is analyzed to characterize the subsurface. Alford rotations are used to determine fracture orientation and demonstrate good results when geophone orientation is taken into account. Results indicate that seismic anisotropy is caused by regional faulting. Average fracture density of less than 5% and water table depth estimates are consistent with field observations. Groundwater flow direction has been observed by others to cross the fault trend and is subparallel to a secondary fracture set. In this study, seismic anisotropy appears unrelated to this secondary fracture set. Vp/Vs and Poisson's ratio values indicate a dolomite lithology. Sledgehammer and Vibroseis data provide consistent results. In the GPR experiment, reflection profiles are acquired through common-offset profiling perpendicular to the dominant flow direction. High frequency waves are used to delineate fluid flow paths through a subhorizontal fracture and observe tracer channeling. Channeling of flow is expected to control solute transport. Changes in radar signal are quantitatively associated with changes in fracture filling material from an innovative method using correlation coefficients. Mapping these changes throughout the survey area reveals the geometry of the flow path of each injected liquid. The tracer is found to be concentrated in the center of the survey area where fracture apertures are large. This demonstrates that spatial variations in concentration are controlled by fluid channel geometry.

No description available.

In subhumid to arid climates throughout the world, recharge to groundwater in urban areas is often found to be higher than pre-urbanization rates, despite an increased percentage of impermeable surfaces. Groundwater recharge in the city of Lethbridge is substantially higher than recharge rates prior to urbanization, resulting in the formation of perched water table conditions. High perched water table conditions, typically at depths between one and 2.5 metres, have created problems for the City and University of Lethbridge, including the increased occurrence of slope failures along nearby coulees. This study estimates of the volume of excess water available for groundwater recharge through the practices of urban turfgrass irrigation, and water storage. Between May and September, 1990-1996 irrigation was applied far above evapotranspiration demands, resulting in large volumes of water available for groundwater recharge in the Varsity Village subdivision of the City of Lethbridge. The relationship between the amount of water applied and the development of perched water table systems was strong enough that equations between inputs and water table depth could be derived, and used to predict water table elevation. / xiii, 190 leaves : ill., maps ; 28 cm.

In the numerical simulation of fluid flow and solute transport in porous media, finite element method (FEM) has long been utilized and has been proven to be efficient. In this work, an alternative approach called random walk particle tracking (RWPT) method is proposed. In this method, a finite number of particles represent the distribution of a solute mass. Each particle carries a certain fraction of the total mass and moves in the porous media according to the velocity field. The proposed RWPT model is established on a scientific software platform OpenGeoSys (OGS), which is an open source initiative for numerical simulation of thermo-hydro-mechanical-chemical (THMC) processes in porous media. The flow equation is solved using finite element method in OGS. The obtained hydraulic heads are numerically differentiated to obtain the velocity field. The particle tracking method does not solve the transport equation directly but deals with it in a physically stochastic manner by using the velocity field. Parallel computing concept is included in the model implementation to promote computational efficiency. Several benchmarks are developed for the particle tracking method in OGS to simulate solute transport in porous media and pore space. The simulation results are compared to analytical solutions and other numerical methods to test the presented method. The particle tracking method can accommodate Darcy flow as it is the main consideration in groundwater flow. Furthermore, other flow processes such as Forchheimer flow or Richards flow can be combined with as well. Two applications indicate the capability of the method to handle theoretical real-world problems. This method can be applied as a tool to elicit and discern the detailed structure of evolving contaminant plumes. / Bei der numerischen Simulation von Strömung und Stofftransport in porösen Medien hat die Nutzung der Finite-Elemente-Methode (FEM) eine lange Tradition und wird sich als effizient erweisen. In dieser Arbeit wird ein alternativer Ansatz, die random walk particle tracking (RWPT) Methode vorgeschlagen. Bei diesem Verfahren stellt eine endliche Anzahl von Partikeln die Verteilung eines gelösten Stoffes dar. Jedes Teilchen trägt einen bestimmten Bruchteil der Gesamtmasse und bewegt sich in den porösen Medien gemäß des Geschwindigkeitsfeldes. Das vorgeschlagene RWPT Modell basiert auf der wissenschaftlichen Softwareplattform OpenGeoSys (OGS), die eine Open-Source-Initiative für die numerische Simulation thermo-hydro-mechanisch-chemischen (THMC) in porösen Medien darstellt. Die Strömungsgleichung wird in OGS mit der Finite-Elemente-Methode gelöst. Der Grundwasserstand wird numerisch berechnet, um das Geschwindigkeitsfeld zu erhalten. Die Partikel-Tracking-Methode löst die Transportgleichung nicht direkt, sondern befasst sich mit ihr in einer physikalisch stochastische Weise unter Nutzung des Geschwindigkeitsfeldes. Zur Berücksichtigung der Recheneffizienz ist ein Parallel Computing-Konzept in der Modell-Implementierung enthalten. Zur Simulation des Stofftransports in porösen Medien und im Porenraum wurden mehrere Benchmarks für die Partikel-Tracking-Methode in OGS entwickelt. Die Simulationsergebnisse werden mit analytischen Lösungen und andere numerische Methoden verglichen, um die Aussagefähigkeit des vorgestellten Verfahrens zu bestätigen. Mit der Partikel-Tracking-Methode kann die Darcy-Strömung gelöst werden, die das wichtigste Kriterium in der Grundwasserströmung ist. Außerdem bewältigt die Methode auch andere Strömungsprozesse, wie die Forchheimer-Strömung und die Richards-Strömung. Zwei Anwendungen zeigen die Leistungsfähigkeit der Methode bei der prinzipiellen Handhabung von Problemen der realen Welt. Die Methode kann als ein Instrument zur Aufdeckung Erkennung der detaillierte Struktur von sich entwickelnden Schadstofffahnenangewendet werden.

<p>The analysis of geologic, hydrologic and hydrogeologic data interpreted to give the characteristics of the Cape Flats aquifer showed the quality of groundwater from the aquifer is suitable for development as a water resource. The conceptual model of the Cape Flats sand shows an unconfined sandy aquifer, grading into semi-confined conditions in some places where thick lenses of clay and peat exists. Recharge rates through the saturated zone of the Cape Flats aquifer have been determined by water table fluctuation (WTF), rainfall-recharge relationship, soil water balance and chloride mass balance methods (CMB). Recharge rates using the WTF vary considerably between wet and dry years and between locations, with a range of 17.3% to 47.5%. Values obtained from empirical rainfall-recharge equation (method 2) agree with those of the WTF. Recharge estimates from the water balance model are comparatively lower but are within the range calculated using empirical method 2 (i.e. 87 &ndash / 194 mm or 4 &ndash / 21% of MAP). These recharge rates also agree with estimates from the series of other methods applied to sites located in the north-western coast of Western Cape and are comparable to recharge rates obtained elsewhere in the world.</p>

A study was undertaken to better understand the hydrogeologic framework of the Delaware County, Indiana. Arc View GIS 3-D and Spatial Analysts along with VISUAL MODFLOW are used to study the flow patterns by developing a 3-D model of major aquifers in the county, both glacial and bedrock. The GIS platform facilitates the time consuming task of preparing data input and output. In addition, major recharge zones are also identified in GIS using soil and slope data.The bedrock of Delaware County is composed of Silurian bedrock, which is overlain by glacial drift. The drift is mostly till that is interbedded with eight sand and gravel layers which are horizontal and discontinuous. Both, bedrock and sand & gravel glacial aquifers are good sources of groundwater in the county. More than 3000 water wells are located in the area. It is very interesting to note that 50% of the water-wells in the area are in confined sand and gravel and the remaining wells are in carbonaceous bedrock. The bedrock contains numerous preglacial karsts valley systems which are areas of high transmissivity solution features, also adding further interest to modeling this region.The hydraulic conductivity of sand and gravel is assumed to be 0.0015 m/s whereas for bedrock 0.00025 m/s is used. The streambeds in the area are also the source of groundwater discharge. The hydraulic conductivity of streambed is equal to 0.00028 m/s. Recharge to the area occurs as net recharge, which varies spatially depending on the nature of soil.Several simplifying assumptions were made for the conceptualization and simulation of flow in the basins. The main assumptions are as follows: 1) groundwater is in steady-state, 2) pumping does not significantly affect the level of hydraulic heads; therefore only high capacity pumping wells are simulated, 3) Net recharge from precipitation varies spatially, 4) flow in the bedrock aquifer occurs in the uppermost 40 m and is horizontal 5) vertical flow is assumed to be controlled by intervening sand and gravel units. Ground-water flow in the basins was conceptualized as a three-dimensional flow system. The model boundaries selected to represent natural hydrologic boundaries include (1) river leakage boundaries along major rivers; (2) a constant head boundary around the aquifer boundaries; (3) a general head boundary along the major streams; (4) drain boundaries along major drains; and 5) Evapotranspiration losses.The simulated region is an area of 398 square miles. The model consists of over 3600 cells and employed a regular grid spacing of 6o x 60. A variable grid was designed to provide additional detail in areas of special interest and thus allowed these areas to be simulated more accurately in the model.The model was calibrated using a manual trial-and-error adjustment of parameters. Hydraulic conductivity values, and streambed conductance were adjusted during successive simulations until the flow pattern matches the regional flow direction as computed from the water-well derived (static water level) potentiometric surface. The computed potentiometric surface is an adequate or reasonable match on a regional scale, with the general trend of SE-NW. It is observed that the model is extremely sensitive to changes in horizontal hydraulic conductivity and recharge in the form of precipitation. The model is least sensitive to streambed vertical hydraulic conductivity.The water budget for the calibrated model represents the distribution of groundwater inflow and outflow during calibration. The data indicate that 16.5% of the inflow to the modeled groundwater system is flow across model boundaries and 83% from effective recharge from precipitation, and the rest from streams and rivers. Outflow consists of 2% pumpage, 5% seepage to streams and drains, and 93% is flow across boundaries.It is expected that this study will be beneficial to improve the understanding of groundwater in Delaware County, including both vertical and horizontal flow and interaction of flow between surface and groundwater. Also, the results of the modeling study can be used as a predictive tool for long-term management and monitoring of water resources in the region. / Department of Geology