This thesis presents two interrelated studies that consider nutrient management and seasonal changes in recharge on agricultural lands within the context of source water protection. The research focuses first on the management of the risk to groundwater quality through the implementation of various nutrient management practices and secondly considers the dynamic nature of the transport pathway to the groundwater system associated with seasonal changes in climate and hydrology. The combined results provide insight into several of the key factors influencing the protection of groundwater sources within the agricultural landscape. Field work was completed between 2009 and 2010 on an agricultural field near the City of Woodstock, Ontario. The site is located within a source water protection area; the two-year travel time zone of the Thornton Well Field which represents the primary water supply for the City of Woodstock and which has experienced chronic increases in nitrate concentrations over the last few decades. The wells are completed in glacial overburden consisting of intermingling sand and gravel till aquifers which overly a limestone bedrock aquifer. Agricultural best or beneficial management practices (BMPs) field have been implemented and monitored since 2004. The BMPs were adopted in order to reduce nitrogen losses to the aquifer, and consisted of a reduction in nitrogen fertilizer application rates over a series of agricultural fields located near the well The first study is a one year experiment designed to compare alternative nutrient management practices for corn. Combinations of fertilizer treatments with or without a legume cover crop (red clover) were assessed. The fertilizer treatments studied were: a polymer coated urea (slow-release fertilizer) applied at planting, a conventional urea applied at planting, side-dress treatment of a solution of urea and ammonium nitrate in water containing 28% nitrogen with two different application rates applied in the early summer, and a control. The legume cover crop was incorporated in the soil in the previous fall, and acts as a slow release fertilizer as nitrogen is made available to the following crop as the plants decompose. Treatments were compared based on crop yield, overall economic return, and the potential for nitrate leaching. The potential for nitrate leaching was evaluated with bi-weekly shallow soil core during the growing season, and deep soil cores taken before planting, after harvest and the following spring. The deep cores allowed changes in nitrate storage below the rooting zone to be assessed. The results of this study highlight the importance of timing of fertilizer applications and rate of fertilizer applications. Treatments which provide a delay in the release or application of fertilizer, the polymer-coated urea, the calculator-rate side-dress and the clover cover crop, were found to be advantageous. The polymer-coated urea treatments and side-dress treatments were found to reduce leaching compared to the conventional urea treatment. Treatments with the clover cover crops were not found to reduce crop yields or increase leaching potential, and lower fertilizer costs associated to this practice were found to have a positive economic effect. Plots treated with the high-rate side-dress fertilizer application lost more nitrate to the subsurface compared to the other treatment options, and an economic disadvantage was observed as yields did not compensate for higher fertilizer costs. The study highlights the advantages of the different treatments under study, which may be used to inform policy makers and farmers in the selection of economically and environmentally sustainable nutrient management BMP options. Groundwater monitoring at the site over the years has indentified interesting recharge dynamics, particularly in the vicinity of an ephemeral stream which develops annually during spring and winter melt events in a low lying area of the study site. It was hypothesized that rapid recharge could occur beneath the stream allowing for surface water to quickly reach groundwater, posing a threat to municipal water wells. The current framework of source water protection does not take into account the potential risk posed by this type recharge event. At this field site, rapid infiltration associated with this type of event may pose a risk to drinking water quality due to the proximity of the stream to the pumping wells and the nature of the aquifer. The second study examines rapid groundwater recharge processes beneath the ephemeral stream during the course of a spring melt in 2010. The goals of the study were to quantify recharge at one location beneath the stream and to assess whether temperature variations above the water table can be used as a tracer to reasonably estimate recharge during a short live recharge event. A novel housing for the temperature sensors was designed in order to deploy and position them into gravelly materials within the vadose zone, which reduced the potential for the formation of preferential pathways and permitted the retrieval of the sensors at a later date. Field data were collected during the course of the spring melt period from a network of groundwater monitoring wells and subsurface temperature sensors. Spatial and temporal changes in groundwater geochemistry, hydraulic head and temperature were were used to characterize recharge dynamics at the field site. Recharge beneath a segment of the ephemeral stream was quantified through the numerical analysis of the field data using Hydrus 1-D, a one-dimensional numerical model designed to simulate soil water flow and heat transport in variably saturated porous media. Site specific data were used to create the model domain, provide estimates of physical parameters, and to define initial and time variable boundary conditions. Model parameters were first calibrated by simulating periods where it was expected that soils would be gravity drained with minimal soil water flow, and then further refined by simulating the period when the ephemeral stream was present. A final set of parameters was determined, and the initial gravity drained conditions were re-simulated. The model was able to reproduce field observations under different flow scenarios using the final set of parameters, suggesting that the conceptual model and final model domain representative of the actual field conditions. The successful simulation of the field data sets under the different flow scenarios also increases confidence in the uniqueness of the model results. The model estimated that 0.15 m of recharge occurred beneath the instrumented site during the period between March 9th and March 22nd of 2010 when the ephemeral stream was present. This represents approximately a third of the expected total annual recharge for this location. Regional changes in hydraulic head, groundwater temperature and groundwater chemistry provided additional insight into the dynamic nature of the recharge process during the spring meld period and further illustrated the spatial variability of the aquifers’ response to the stream. The study found that the use of temperature as a tracer provided useful and quantifiable insight into recharge phenomena. The results of this study suggest that high rates of rapid recharge occur beneath the ephemeral stream, and are spatially variable. This type of focused infiltration that occurs during the spring melt may represent a risk to municipal water quality if the infiltrating waters are carrying contaminants.

This dissertation reports on the methods and results of a three-phased investigation to estimate the annual volume of ephemeral-channel-focused groundwater recharge attributable to urbanization (urban-enhanced groundwater recharge) in the Sierra Vista subwatershed of southeastern Arizona, USA. Results were used to assess a prior estimate. The first research phase focused on establishment of a study area, installation of a distributed network of runoff gages, gaging for stage, and transforming 2008 stage data into time series of volumetric discharge, using the continuous slope-area method. Stage data were collected for water years 2008 - 2011. The second research phase used 2008 distributed runoff data with NWS DOPPLER RADAR data to optimize a rainfall-runoff computational model, with the aim of identifying optimal site-specific distributed hydraulic conductivity values and model-predicted infiltration. The third research phase used the period-of-record runoff stage data to identify study-area ephemeral flow characteristics and to estimate channel-bed infiltration of flow events. Design-storm modeling was used to identify study-area predevelopment ephemeral flow characteristics, given the same storm event. The difference between infiltration volumes calculated for the two cases was attributed to urbanization. Estimated evapotranspiration was abstracted and the final result was equated with study-area-scale urban-enhanced groundwater recharge. These results were scaled up to the Sierra Vista subwatershed: the urban-enhanced contribution to groundwater recharge is estimated to range between 3270 and 3635 cubic decameters (between 2650 and 2945 acre-feet) per year for the period of study. Evapotranspirational losses were developed from estimates made elsewhere in the subwatershed. This, and other sources of uncertainty in the estimates, are discussed and quantified if possible.

This thesis presents two interrelated studies that consider nutrient management and seasonal changes in recharge on agricultural lands within the context of source water protection. The research focuses first on the management of the risk to groundwater quality through the implementation of various nutrient management practices and secondly considers the dynamic nature of the transport pathway to the groundwater system associated with seasonal changes in climate and hydrology. The combined results provide insight into several of the key factors influencing the protection of groundwater sources within the agricultural landscape. Field work was completed between 2009 and 2010 on an agricultural field near the City of Woodstock, Ontario. The site is located within a source water protection area; the two-year travel time zone of the Thornton Well Field which represents the primary water supply for the City of Woodstock and which has experienced chronic increases in nitrate concentrations over the last few decades. The wells are completed in glacial overburden consisting of intermingling sand and gravel till aquifers which overly a limestone bedrock aquifer. Agricultural best or beneficial management practices (BMPs) field have been implemented and monitored since 2004. The BMPs were adopted in order to reduce nitrogen losses to the aquifer, and consisted of a reduction in nitrogen fertilizer application rates over a series of agricultural fields located near the well The first study is a one year experiment designed to compare alternative nutrient management practices for corn. Combinations of fertilizer treatments with or without a legume cover crop (red clover) were assessed. The fertilizer treatments studied were: a polymer coated urea (slow-release fertilizer) applied at planting, a conventional urea applied at planting, side-dress treatment of a solution of urea and ammonium nitrate in water containing 28% nitrogen with two different application rates applied in the early summer, and a control. The legume cover crop was incorporated in the soil in the previous fall, and acts as a slow release fertilizer as nitrogen is made available to the following crop as the plants decompose. Treatments were compared based on crop yield, overall economic return, and the potential for nitrate leaching. The potential for nitrate leaching was evaluated with bi-weekly shallow soil core during the growing season, and deep soil cores taken before planting, after harvest and the following spring. The deep cores allowed changes in nitrate storage below the rooting zone to be assessed. The results of this study highlight the importance of timing of fertilizer applications and rate of fertilizer applications. Treatments which provide a delay in the release or application of fertilizer, the polymer-coated urea, the calculator-rate side-dress and the clover cover crop, were found to be advantageous. The polymer-coated urea treatments and side-dress treatments were found to reduce leaching compared to the conventional urea treatment. Treatments with the clover cover crops were not found to reduce crop yields or increase leaching potential, and lower fertilizer costs associated to this practice were found to have a positive economic effect. Plots treated with the high-rate side-dress fertilizer application lost more nitrate to the subsurface compared to the other treatment options, and an economic disadvantage was observed as yields did not compensate for higher fertilizer costs. The study highlights the advantages of the different treatments under study, which may be used to inform policy makers and farmers in the selection of economically and environmentally sustainable nutrient management BMP options. Groundwater monitoring at the site over the years has indentified interesting recharge dynamics, particularly in the vicinity of an ephemeral stream which develops annually during spring and winter melt events in a low lying area of the study site. It was hypothesized that rapid recharge could occur beneath the stream allowing for surface water to quickly reach groundwater, posing a threat to municipal water wells. The current framework of source water protection does not take into account the potential risk posed by this type recharge event. At this field site, rapid infiltration associated with this type of event may pose a risk to drinking water quality due to the proximity of the stream to the pumping wells and the nature of the aquifer. The second study examines rapid groundwater recharge processes beneath the ephemeral stream during the course of a spring melt in 2010. The goals of the study were to quantify recharge at one location beneath the stream and to assess whether temperature variations above the water table can be used as a tracer to reasonably estimate recharge during a short live recharge event. A novel housing for the temperature sensors was designed in order to deploy and position them into gravelly materials within the vadose zone, which reduced the potential for the formation of preferential pathways and permitted the retrieval of the sensors at a later date. Field data were collected during the course of the spring melt period from a network of groundwater monitoring wells and subsurface temperature sensors. Spatial and temporal changes in groundwater geochemistry, hydraulic head and temperature were were used to characterize recharge dynamics at the field site. Recharge beneath a segment of the ephemeral stream was quantified through the numerical analysis of the field data using Hydrus 1-D, a one-dimensional numerical model designed to simulate soil water flow and heat transport in variably saturated porous media. Site specific data were used to create the model domain, provide estimates of physical parameters, and to define initial and time variable boundary conditions. Model parameters were first calibrated by simulating periods where it was expected that soils would be gravity drained with minimal soil water flow, and then further refined by simulating the period when the ephemeral stream was present. A final set of parameters was determined, and the initial gravity drained conditions were re-simulated. The model was able to reproduce field observations under different flow scenarios using the final set of parameters, suggesting that the conceptual model and final model domain representative of the actual field conditions. The successful simulation of the field data sets under the different flow scenarios also increases confidence in the uniqueness of the model results. The model estimated that 0.15 m of recharge occurred beneath the instrumented site during the period between March 9th and March 22nd of 2010 when the ephemeral stream was present. This represents approximately a third of the expected total annual recharge for this location. Regional changes in hydraulic head, groundwater temperature and groundwater chemistry provided additional insight into the dynamic nature of the recharge process during the spring meld period and further illustrated the spatial variability of the aquifers’ response to the stream. The study found that the use of temperature as a tracer provided useful and quantifiable insight into recharge phenomena. The results of this study suggest that high rates of rapid recharge occur beneath the ephemeral stream, and are spatially variable. This type of focused infiltration that occurs during the spring melt may represent a risk to municipal water quality if the infiltrating waters are carrying contaminants.

Base flow recession curve during a dry period is a distinct hydrologic signature of a watershed. The base flow recession analysis for both streamflow and spring flow has been extensively studied in the literature. Studies have shown that the recession behaviors during the early stage and the late stage are different in many watersheds. However, research on the transition from early stage to late stage is limited and the hydrologic control on the transition is not completely understood. In this dissertation, a novel cumulative regression analysis method is developed to identify the transition flow objectively for individual recession events in the well-studied Panola Mountain Research Watershed in Georgia, USA. The streamflow at the watershed outlet is identified when the streamflow at the perennial stream head approaches zero, i.e., flowing streams contract to perennial streams. The identified transition flows are then compared with observed flows when the flowing stream contracts to the perennial stream head. As evidenced by a correlation coefficient of 0.90, these two characteristics of streamflow are found to be highly correlated, suggesting a fundamental linkage between the transition of base flow recession from early to late stages and the drying up of ephemeral streams. At the early stage, the contraction of ephemeral streams mostly controls the recession behavior. At the late stage, perennial streams dominate the flowing streams and groundwater hydraulics governs the recession behavior. The ephemeral stream densities vary from arid regions to humid regions. Therefore, the characteristics of transition flow across the climate gradients are also tested in 40 watersheds. It is found that climate, which is represented by climate aridity index, is the dominant controlling factor on transition flows from early to late recession stages. Transition flows and long-term average base flows are highly correlated with a correlation coefficient of 0.82. Long-term average base flow and the transition flow of recession are base flow characteristics at two temporal scales, i.e., the long-term scale and the event scale during a recession period. This is a signature of the co-evolution of climate, vegetation, soil, and topography at the watershed scale. The characteristics of early and late recession are applied for quantifying human impacts on streamflow in agricultural watersheds with extensive groundwater pumping for irrigation. A recession model is developed to incorporate the impacts of human activities (such as groundwater pumping) and climate variability (such as evapotranspiration) on base flow recession. Groundwater pumping is estimated based on the change of observed base flow recession in watersheds in the High Plains Aquifer. The estimated groundwater pumping rate is found consistent compared with the observed data of groundwater uses for irrigation. Besides streamflow recession analysis, this dissertation also presents a novel spring recession model for Silver Springs in Florida by incorporating groundwater head, spring pool altitude, and net recharge into the existing Torricelli model. The results show that the effective springshed area has continuously declined since 1988. The net recharge has declined since the 1970s with a significant drop in 2002. Subsequent to 2002, the net recharge increased modestly but not to the levels prior to the 1990s. The decreases in effective springshed area and net recharge caused by changes in hydroclimatic conditions including rainfall and temperature, along with groundwater withdrawals, contribute to the declined spring flow.