ROCK OUTCROPS IN THE CANADIAN SHIELD: AN INVESTIGATION OF CONTAMINANT TRANSPORT FROM SURFACE SOURCES IN FRACTURED ROCK AQUIFERS

Flow, recharge and transport dynamics in fractured rock aquifers with low lying rock outcrops is a largely unexplored area of study in hydrogeology. The purpose of this thesis is to examine these topics in an agricultural area in Eastern Ontario. The study consists of a regional scale groundwater quality study, an infiltration experiment that considers bacteria transport from the ground surface to a well, and a numerical modelling study that tests the parameters that affect surface infiltration of a tracer from a rock outcrop to a deeper horizontal fracture. In the water quality study, approximately 65% of the samples contained total coliform, 16% contained E. coli, and 1% contained nitrate-N at greater than 5 mg/L. Occurrence of E. coli increased when considering seasonality, where wells were drilled on rock outcrops, and for shallow well intervals. Nitrate-N did not occur above the Guidelines for Canadian Drinking Water Quality (Health Canada, 2012) of 10 mg/L. Rapid arrival times were observed in the infiltration study for both the microspheres (30 minutes) and a dye tracer (45 minutes) in a well approximately 6.0 m in horizontal and 2.8 m in vertical distance from the tracer source. Transport velocities were approximately 38.9 m/day for the dye tracer and 115.2 m/day for the colloidal tracer. Results of the model runs indicate that overburden can provide an effective protective layer to transport in fractures, that high groundwater velocities occur in larger fracture apertures and higher gradients dilute tracer concentrations, and that lower groundwater velocities occur with smaller fracture apertures and lower gradients result in elevated tracer concentrations. Lower rainfall rates, larger fracture apertures, early tracer time, larger gradients, and lower water levels maintained unsaturated conditions for longer time periods such that tracer transport was delayed until saturated conditions were attained. The overall heterogeneity of this aquifer environment creates a source water protection conundrum where the water quality is generally good, while transport can occur very quickly in proximity to rock outcrops and in areas with limited overburden.

Surface wave and storm surge prediction during Hurricane Sandy

Hurricane Sandy was the largest storm on historical record in the Atlantic Ocean basin with extensive coastal damage caused by large waves and high storm surge. The primary objectives of this thesis are to compare and evaluate three different spatially-varying surface wind fields of Hurricane Sandy to investigate the impact of the differences between the complex wind fields on predictions of the sea surface evolution, and to evaluate the impact of the storm on the hydrodynamics in Great South Bay (GSB) and the discharge of ocean water into the back-barrier bay from overwash over Fire Island. Three different spatially-varying surface wind fields were evaluated and compared to wind observations, including the parametric Holland (1980) model (H80), the parametric Generalized Asymmetric Holland Model (GAHM), and results from the WeatherFlow Regional Atmospheric Modelling System (WRAMS). The winds were used to drive the coupled Delft3D-SWAN hydrodynamic and ocean wave models on a regional grid. The results indicate that the WRAMS wind field produces wave model predictions in the best agreement with significant wave height observations, followed by the GAHM and H80 wind fields and that a regional atmospheric wind model is best for hindcasting hurricane waves and water levels when detailed observations are available, while a parametric vortex model is best for forecasting hurricane sea surface conditions. Using a series of four connected Delft3D-SWAN grids to achieve finer resolution over Fire Island and GSB, a higher resolution WRAMS was used to predict waves and storm surge. The results indicate that strong local winds have the largest influence on water level fluctuations in GSB. Three numerical solutions were conducted with varying extents of barrier island overwash. The simulations allowing for minor and major overwash indicated good agreement with observations in the east end of GSB and suggest that island overwash provided a significant contribution of ocean water to GSB during the storm. Limiting the overwash in the numerical model directly impacts the total discharge into GSB from the ocean through existing inlets. The results of this study indicate that barrier island overwash had a significant impact on the water levels in eastern GSB.