Ground-water conductivity data were obtained from shallow wells in a 12 km2 stream-basin along a 400 m transect, extending from the divide to the stream. The stream, Pringle Branch, is a second-order perennial stream in Hillsborough County, Florida. The shallow stratigraphy consists of 2-3 m of fine sand over a layer of clayey silt and silty clay. Vegetation cover includes grasses on the upper and middle slope, and riparian woodlands on the foot slope and floodplain. Precipitation is about 1.3 m per year. Shallow ground-water conductivity is about 50 uS/cm at the divide. It increases moderately along the mid slope, then increases markedly within the riparian woodlands, reaching a maximum of about 500 uS/cm at 30m from the stream and then decreases to about 150 uS/cm at the stream. The spatial variation of terrain electrical conductivity data collected using electromagnetic methods (EM 31) is similar to the spatial variation of ground-water conductivity.Dry season through wet season monitoring shows that ground-water conductivity in each well varies about 40%, generally following variations in potential evapotranspiration (ETpan). The more than five-fold increase in ground-water conductivity from divide to riparian woodlands is maintained during both dry and wet seasons. The ground-water conductivity in this basin appears to be determined principally by spatial variations in ET and not by temporal variations in ET or interaction with soil minerals. The data suggest that patterns of ground-water conductivity can be used to infer patterns of ET variation within a small basin. A mass transport model constructed to test the hypothesis that evapotranspiration has the dominant effect on ground-water conductivity closely duplicates the observed variation in ground-water conductivity from divide to stream. The model uses two evaporation rates, 0.73 m/y for the grasses and 1.46 m/y for the riparian woodlands, and no contribution from solution of matrix materials.

Potential contamination from hazardous and solid waste landfills stemming from population increase, rapid industrialization, and the proliferation of assembly plants known as the maquiladoras, are of major concern in the U.S.-Mexican border area. Additionally, historical, current, and future stresses on the Hueco Bolson alluvial aquifer in the El Paso/Ciudad Juarez area due to excessive groundwater withdrawal can affect contaminant migration in the area. In the current study, an updated and improved threedimensional numerical groundwater flow and transport model is developed using a current Hueco Bolson groundwater availability model as its basis. The model with contaminant transport is required to access and characterize the extent of vulnerability of the aquifer to potential contamination from landfills in the El Paso/Ciudad Juarez border area. The model developed in this study is very capable of serving as the basis of future studies for water availability, water quality, and contamination assessments in the Hueco Bolson. The implementation of fate and transport modeling and the incorporation of the Visual MODFLOW® pre and post processor, requiring MODFLOW 2000 data conversion, enabled significant enhancements to the numerical modeling and computing capabilities for the Hueco Bolson. The model in the current research was also developed by employing MT3DMS©, ZONEBUDGET, and Visual PEST® for automated calibrations. Simulation results found that the Hueco Bolson released more water from storage than the aquifer was being recharged in response to increased pumping to supply the growing border area population. Hence, significant head drops and high levels of drawdown were observed in the El Paso/Ciudad Juarez area. Predictive simulations were completed representing scenarios of potential contamination from the border area sites. Fate and transport results were most sensitive to hydraulic conductivities, flow velocities, and directions at the sites. Sites that were located within the vicinity of the El Paso Valley and the Rio Grande River, where head differences and permeabilities were significant, exhibited the highest potentials for contaminant migration.

The overall goal of this study was to develop a comprehensive tool for assessing the effectiveness of selected BMPs on both hydrology and water quality and to demonstrate the applicability of the system by considering 1) temporally and spatially changing land use management practice in an agricultural watershed and 2) interaction between surface and ground water over the entire system. A user interface and Dynamic Agricultural Non-point Source Assessment Tool (DANSAT) were developed to achieve this goal. DANSAT is the only distributed-parameter, physically-base, continuous-simulation, and multi-soil layer model for simulating impacts of agricultural BMPs on hydrology and water qulality in small agricultural watersheds. DANSAT was applied to QNB plot (18m à 27m) and two agricultural watersheds in Virginia, including Owl Run watershed (1140 ha) and QN2 in the Nomini Creek watershed (216 ha), to evaluate the model components and its performance in predicting runoff, sediment yield, and pesticide load. DANSAT performed well in predicting total runoff and temporal variations in surface runoff for both field-scale and watershed-scale applications. Total percent errors between the measured and predicted results were less than 10% except for one case (39.8% within a subwatershed of Owl Run watershed), while the daily Nash-Sutcliffe model efficiencies were greater than 0.5 in all applications. Predicted total sediment yields were within ±35% of observed values in all applications. However, the performance of DANSAT in predicting temporal trend and spatial distribution of sediment loads was acceptable only within Owl Run watershed, where high correlations between flow rates and sediment loads exist. The predicted total pesticide loads were within ±100% of observed values. DANSAT failed to simulate the temporal occurrence of pesticide loads with a 0.42 daily Nash-Sutcliffe efficiency value. The Dual-Simulation (DS) was developed within the linked ground water approach to resolve problems encountered due to the existence of different temporal scales between DANSAT and the existing ground water models such as MODFLOW and MT3D. The linked approach performed better in predicting the seasonal trend of total runoff compared to the integrated approach by showing an increase in monthly Nash-Sutcliffe efficiency value from 0.53 to 0.60. Surface and subsurface output variables were sensitive to the changes in spatially distributed soil parameters such as total porosity and field capacity. A maximum grid size of 100 m was recommended to be appropriate for representing spatial distribution of topographic, land use, and soil characteristics based on accuracy analysis during the GIS manipulation processes. Larger time-step based on predefined acceptable maximum grid size, decreased the computational time dramatically. Overall sensitivity to different grid sizes and time-steps was smallest for hydrology components followed by sediment and pesticide components. Dynamic crop rotation was considered by DANSAT, and the model successfully simulated the impacts of temporal and spatial changes in crop rotations on hydrology and water quality for both surface and subsurface areas. DANSAT could prove to be a useful tool for non-point source pollution managers to assess the relative effectiveness of temporally and spatially changing BMPs on both surface and ground water quantity and quality. / Ph. D.

The need to consider groundwater and surface water as a single resource has fostered the interest of the scientific community on the interactions between surface water and groundwater. The region below and alongside rivers where surface hydrology and subsurface hydrology concur is the hyporheic zone. This is the region where water exchange determines many biogeochemical and ecological processes of great impact on the functioning of rivers. However, the complex processes taking place in the hyporheic zone require a multidisciplinary approach. The combination of innovative point and distributed techniques originally developed in separated disciplines is of great advantage for the indirect identification of water exchange in the hyporheic zone. Distributed techniques using temperature as a tracer such as fiber-optic distributed temperature sensing can identify the different components of groundwater-surface water interactions based on their spatial and temporal thermal patterns at the sediment-water interface. In particular, groundwater, interflow discharge and local hyporheic exchange flows can be differentiated based on the distinct size, duration and sign of the temperature anomalies. The scale range and resolution of fiber-optic distributed temperature sensing are well complemented by geophysics providing subsurface structures with a similar resolution and scale. Thus, the use of fiber-optic distributed temperature sensing to trace flux patterns supported by the exploration of subsurface structures with geophysics enables spatial and temporal investigation of groundwater-surface water interactions with an unprecedented level of accuracy and resolution. In contrast to the aforementioned methods that can be used for pattern identification at the interface, other methods such as point techniques are required to quantify hyporheic exchange fluxes. In the present PhD thesis, point methods based on hydraulic gradients and thermal profiles are used to quantify hyporheic exchange flows. However, both methods are one-dimensional methods and assume that only vertical flow occurs while the reality is much more complex. The study evaluates the accuracy of the available methods and the factors that impact their reliability. The applied methods allow not only to quantify hyporheic exchange flows but they are also the basis for an interpretation of the sediment layering in the hyporheic zone. For upscaling of the previous results three-dimensional modelling of flow and heat transport in the hyporheic zone combines pattern identification and quantification of fluxes into a single framework. Modelling can evaluate the influence of factors governing groundwater-surface water interactions as well as assess the impact of multiple aspects of model design and calibration of high impact on the reliability of the simulations. But more importantly, this modelling approach enables accurate estimation of water exchange at any location of the domain with unparalleled resolution. Despite the challenges in 3D modelling of the hyporheic zone and in the integration of point and distributed data in models, the benefits should encourage the hyporheic community to adopt an integrative approach comprising from the measurement to the upscaling of hyporheic processes.