Two gravel pit lakes in central Texas were examined over the course of two years with upgradient and downgradient piezometer installations. Groundwater and lake water were sampled bimonthly for nutrients, water levels, and groundwater chemistry, and in addition, rain and lake gauges and mini-piezometers were installed, depth surveys were conducted, and a simple 2D flow model was constructed. The project goal was evaluating and examining flow dynamics and chemical effects as groundwater flows to surface water and back to groundwater with the intent to understand the effects that gravel pit lake systems have on connected shallow groundwater. Both lake systems were shown to be flow-through systems that influence the water quality by decreasing the dissolved nutrients in the groundwater in their vicinity while oxygenating the water and altering the pH. However, the lakes are also prone to high levels of evaporation, meaning that minor improvements to water quality come at the cost of decreasing the quantity of water in storage within the aquifer. Similar groundwater and mine lake systems may show comparable tendencies, providing new information for water managers, regulators, and stake-holders about the potential roles that non-remediated gravel pit lakes may play in local ecosystems and aquifer dynamics.

To investigate the influence of flow rate increment on the solute transport parameter of immobile zones in a karst system, a dye tracer test was conducted in the Downtown Salado Spring Complex (DSSC) comprising three springs: Big Boiling, Anderson, and Doc Benedict springs. The Multiflow two-region nonequilibrium model (2RNE) was used to simulate the breakthrough curve (BTC) of the springs, and changes in the solute transport parameters in response to flow rate increment were observed. The simulation result showed that the 2RNE model was capable of reproducing the BTC of all the DSSC springs, with an R-squared value greater than 0.9 in all flow rate increment scenarios. The research demonstrates that a positive correlation will exist between the flow rate and solute transport parameter of the immobile zones if the tracer transport to the spring is truly influenced by immobile zones. In contrast, a negative correlation will exist between the flow rate and mass transfer coefficient if the immobile zone has less influence. Overall, the research provides insights into contaminant movement in karst by documenting how tracers are retained in the immobile fluid zone.

Removal of CBOD(5) and nitrogen from septic tank effluent was evaluated in four horizontal subsurface flow (HSSF) wetlands. An intermittently loaded cell was compared to a continuously loaded control cell, with both treatments receiving the same weekly volume. The intermittent cell was rapidly drained and "rested" for 24-hr, then refilled in steps, twice weekly. Two media with different particle sizes but similar porosities were also compared. The two media, light weight expanded shale and gravel, were both continuously loaded. As hypothesized, the wetland cell that was intermittently loaded had higher dissolved oxygen, greater ammonia removal, and greater nitrate production than the continuously loaded cells. Areal NH(3)-N removal for the intermittently loaded cell was 0.90 g m(-2) d(-1) compared to 0.47 g m(-2) d(-1) for the control. Ammonia removal was also higher in continuously loaded gravel cells than in cells with expanded shale. Ammonia-N removal was an order of magnitude lower in a similar SSF wetland that had been in operation for 3 years. However, CBOD(5), total suspended solids, and total nitrogen did not vary substantially among the treatments.

The Naval Weapons Industrial Reserve Plant (NWIRP) in McGregor, Texas began manufacturing explosives in 1980 and several hazardous chemicals were discovered in lakes and streams surrounding the site in 1998. Contaminants traveled to local lakes and streams much faster than initially predicted. This research estimated contaminant travel times and identified locations where monitoring wells should be installed to yield the greatest reductions in uncertainties in travel-time predictions. To this end, groundwater and particle-tracking models for NWIRP site were built to predict hydraulic heads and contaminant travel times. Next, parameter (hydraulic conductivities) uncertainties, parameter identifiabilities, observation (hydraulic heads) worth, and predictive (contaminant travel times) uncertainties were quantified. Parameter uncertainties were reduced by up to 92%; a total of 19 of 158 parameters were at least moderately identifiable; travel-time uncertainties were reduced up to 92%. Additionally, travel-time predictions and post-calibration parameter distributions were generated using the null-space Monte Carlo (NSMC) technique. NSMC predicted that conservative tracers exited the flow system within a year, which matches with field data. Finally, an observations-worth analysis found that additional 11 more measurements would reduce travel-time uncertainties by factors from 1.04 to 4.3 over existing data if monitoring wells were installed at the suggested locations.

This study investigates the hydraulic connection of surface karst features within the Northern segment of the Edwards Balcones Fault Zone Aquifer, using a combination of 2D and pseudo-3D Electrical Resistivity Tomography (ERT) at an outcrop near Salado, Texas. The study site features several surface fractures whose hydrological functions are not well understood. Nine ERT profiles and two pseudo-3D models were used to evaluate the connection between surface fractures and subsurface karst conduits. Karst features at the study site were physically evaluated using characteristics such as morphology, which resulted in the identification of three surface fractures (F1, F2, and F3). The ERT results showed several high-resistivity anomalies interpreted as a poorly fractured zone and low-resistivity water-filled conduits within the Edwards Formation. Furthermore, the result reveals that slow hydraulic connectivity exists in F1 and F2; however, F3 presents a low-resistivity zone that extends vertically into the subsurface, which suggests that F3 may serve as a potential recharge feature to the Edwards Aquifer. These findings are corroborated by a water percolation test, as water penetrated more at F3 compared to F1 and F2. This study showed that the combined application of 2D and pseudo-3D ERT can successfully delineate potential recharge pathways in an exposed karst system, thereby constituting a supportive approach providing critical insight into recharge and the vulnerability of karst aquifers to contamination.