Streams draining coniferous plantations contain higher loadings of hydrogen ion, aluminium, sulphate and nitrate, in comparison with streams in adjacent grasslands. Almost all of this ion-load is transported to streams via subsurface water-pathways. An incontrovertible, physical characterization of these pathways within a natural, layered hillslope, has yet to be presented. This research has sought to provide such a characterization for two hillslopes - one afforested with conifers, the other an improved grassland. Much of the uncertainty associated with the identification of soil-water-pathways stems from an inadequate characterization of the errors imposed by the use of each measurement technique. This research has, therefore, compared the predictions of a number of quasi-independent field and analytical techniques, to attempt to lessen the impact of measurement error upon the observed response of the true hydrological system. The impact of conifers upon the detailed water-pathways and lumped catchment response was monitored to educe any changes in the hydrological response which could account for the increased loading of acidic solutes within forest streams. The results of the analysis, indicated that the pathways of water through hillslopes could be predicted from the response of hydrological properties averaged over control volumes of soil-pores. The accuracy of these solutions was proven by the concordance of the response of all of the properties contained within the Darcy-Richards equation. The marked horizon development within the ferric podzol soil of the instrumented forest hillslope, in particular the presence of an indurated B horizon, deflects most percolation laterally within the 0/A and A/E horizons. This pathway was indicted by the results of techniques which included numerical and approximative calculations, discontinuities between the state-dependent hydraulic conductivity of each soil horizon, and the generation of steep, vertical potential gradients in layered porous media. The instrumented grassland hillslope was ploughed 11 years prior to instrumentation. This greatly increased the conductivity of the controlling B horizon, allowing almost all flow to percolate to depth. During winter-storms, the forest hillslope generated flows smaller than those within the grassland hillslope, concomitant with the 29 percent difference in the rainfall-runoff behaviour of the catchment areas. This increased loss of runoff within the afforested areas, may result from the high losses of wetted-canopy-evaporation (39 percent of gross-precipitation) from the Sitka spruce (Picea sitchensis. Bong. Carr.) trees. Individual conifer trees growing on the steep, ferric podzol hillslope appeared to enhance the lateral deflection of flow within the O/A and A/E horizons, probably as a result of their platy root systems, and the high rates of precipitation input to soil at the stem-base. The enhancement of both lateral near-surface flow and below-canopy ion concentrations could, therefore, generate the chemical signatures characteristic of streams draining coniferous forests.

Debris flows are a known hazard in southern California where growing numbers of people are moving into the urban-wildland interface, threatening lives and property. A common location to see a debris flow head scarp is the upper one-third to one-half of an unburned slope at or near the head of a first-order catchment, particularly in areas of relatively shallow soils overlying bedrock. Unburned, relatively steep slopes with gently rounded shoulders and thin soil over bedrock in southern California were investigated to determine if there is a position on these types of slopes where near-surface water levels and the associated pore pressures are relatively and consistently higher during and after rainfall events than the rest of the slope, resulting in an area of preferential shallow slope failure and debris flow initiation. It was hypothesized that this position, if it exists, would be on the upper one-third to one-half of the slope near a change from a shallower slope to a steeper slope (the slope shoulder). It was further hypothesized that elevated subsurface pore pressures at this location would contribute to it being an area of preferential shallow slope failure. The near-surface water levels at two field sites in southern California were monitored for three field seasons. In the laboratory, a meso-scale simulator was constructed and used to replicate field conditions using an adjustable artificial slope and simulated rainfall. The field research showed that areas of higher water levels can exist on the upper one-third to one-half of hillslopes meeting the designated criteria. The laboratory simulations showed elevated water levels in the same general area as the field data. Laboratory simulations also suggested that this is an area of preferential shallow slope failure. The angle of the slope influenced how long a slope took to fail and how much water was needed to do so, with gentler slopes requiring more time and approximately double the amount of water than steeper slopes.

No description available.

Understanding the controls on water movement on forested uplands is critical in predicting the potential effects of disturbance on the sustainability of water resources. I examined the controls on vertical and lateral water movement on forested uplands on a range of landforms (coarse textured outwash, fine textured moraine) and time periods (individual events, during snowmelt, through the growing season, annually, and long-term) at the Utikuma Region Study Area (URSA) on the sub-humid Boreal Plains of Alberta, Canada. To quantify vertical and lateral water movement, hydrometric and tracer measurements were made under natural and experimental conditions at plot and hillslope scales. Vertical flow and unsaturated zone storage dominated hydrologic response to snowmelt and rainfall at the plot and hillslope scales. Plot-scale snowmelt infiltration was greater than near-surface runoff, and when runoff occurred it was limited to south-facing outwash hillslopes underlain by concrete frost. Rainfall simulation studies showed that even under the extreme conditions tested, vertical flow and storage dominated the hydrologic response. Soils at field capacity and precipitation inputs of 15-20 mm or greater at high intensities were required to generate lateral flow via the transmissivity feedback mechanism. The threshold soil moisture and precipitation conditions are such that lateral flow will occur infrequently under natural conditions. Seasonal vertical water movement under natural conditions was greater on outwash than moraine uplands. The maximum downward vertical movement occurred in response to snowmelt, with little subsequent movement over the growing season. Recharge following snowmelt was similar for outwash and moraine sites and was followed by declining water tables through the growing season. Tracer estimates of long-term root zone drainage were low, while estimates of recharge for the moraine were high, raising questions about the appropriateness of this method for these sites. These results emphasize the dominance of vertical relative to lateral water flow on Boreal Plain uplands. Detailed understanding of the controls on water movement can be used to predict the potential effects of disturbance on hydrology and water resources. / Ecology

Nitrate (NO⁻3) export can vary widely among forested watersheds with similar nitrogen loading, geology, and vegetation, which suggests the importance of understanding differing internal retention mechanisms. Transport should be studied at the hillslope scale because the hillslope is the smallest unit with spatial and temporal resolution to reflect many relevant NO⁻3 retention and transport (flow-generation) processes, and headwater forested watersheds are largely comprised of sections of hillslopes. I conducted two experiments to elucidate subsurface flow dynamics and NO⁻3 transport and retention mechanisms on a constructed experimental hillslope model. In the first experiment, I tested whether decadal pedogenetic changes in soil properties in the experimental hillslope used by Hewlett and Hibbert (1963) would lead to changes in recession flow. I repeated (twice) their seminal experiment, whose results led to the development of the Variable Source Area paradigm, by also saturating, covering, and allowing the experimental hillslope to drain until it no longer yielded water. In the historical experiment there was fast drainage for 1.5 d, followed by slow drainage for ~140 d, which led the authors to conclude that recession flow in unsaturated soil could sustain baseflow throughout droughts. This long, slow drainage period was not reproduced in my experiments. Shapes of the drainage curves in my experiments were similar to the historical curve, but slow drainage was truncated, ending after 17 and 12 d, due likely to a leak in the boundary conditions, rather than to pedogenetic changes since the historical experiment. Leakage to bedrock, analogous to the leak in the hillslope model, is a commonly observed phenomenon and this study highlights how that can reduce drainage duration and the contribution of moisture from soils to support baseflow. In the second experiment, I tested whether movement of NO⁻3, which is considered a mobile ion, would be delayed relative to movement of water through a hillslope. I added concentrated pulses of ¹⁵NO⁻3 and a conservative tracer (²H₂O) on the same experimental hillslope, which was devegetated and irrigated at hydrologic steady state. Retention of the ¹⁵NO⁻3 tracer was high in the soil surface (0–10 cm) layer directly where the tracer was added. The portion of the ¹⁵NO⁻3 tracer that passed through this surface layer was further retained/removed in deeper soil. The reduction in the peaks in δ¹⁵N breakthrough was an order of magnitude larger than in δ₂H breakthrough at the outlet 5 m downslope of the tracer addition. The peaks in δ¹⁵N were also delayed relative to the peaks in δ₂H by 1, 6, 9 and 18.5 d for slope distances of 0, 2, 4, and 5 m, respectively, from tracer addition to the outlet. The excess mass of ¹⁵NO⁻3 recovered at the outlet was less than 3% of the original tracer mass injected. Nitrification and denitrification were estimated to be roughly 1:1 and were large fluxes relative to lateral transport into and out of the riparian zone. This tracer experiment shows that bedrock leakage, coupled with multiple retention/removal mechanisms can significantly delay export of added NO⁻3 with implications of additional NO⁻3 sink strength at the watershed scale. / Ph. D. / Nitrate (NO₃⁻) export can vary widely among forested watersheds with similar nitrogen loading, geology, and vegetation, which suggests the importance of understanding differing internal process mechanisms. I conducted two experiments to illustrate how water and NO₃⁻ moved on a constructed hillslope model. In the first experiment, I quantified differences in soil properties in the hillslope model used by Hewlett and Hibbert (1963). Then I repeated (twice) the seminal drainage experiment described in Hewlett and Hibbert (1963). The same hillslope (21.8°; 40%) was wetted up, covered, and allowed to drain until water stopped exiting at the outlet. In the historical experiment there was fast drainage for 1.5 d, followed by slow drainage for ~140 d, which led the authors to hypothesize that slow drainage in surface soil could continually contribute water to streams even during droughts. This long, slow drainage period was not reproduced in my experiments. Drainage was similar at the beginning of drainage between my experiments and the historical experiment, but in my experiment the slow drainage ended earlier (after 17 and 12 d) due likely to a leak in the constructed hillslope model, rather than to significant changes that occurred in the soil itself since the original experiment. This leak in the hillslope model is similar to leakage to bedrock, which is commonly observed in natural hillslopes. In the second experiment, I tested whether NO₃⁻ and water would move through a hillslope at the same rate. I added concentrated pulses of NO₃⁻ (as ¹⁵NO₃⁻ and water (as ²H₂O) on the same devegetated experimental hillslope. Retention of the ¹⁵NO₃⁻ tracer was high in the surface (0–10 cm) where the tracer was added, with little change in the immediately surrounding soil, despite high rates of water input immediately after tracer addition and throughout the experiment. The portion of the ¹⁵NO₃⁻ tracer that passed through the surface layer was further processed by microbes in deeper soil as it traveled downslope. This body of work shows that bedrock leakage, coupled with multiple retention mechanisms throughout the soil profile, can significantly delay export of added NO₃⁻ at the watershed scale.

For this dissertation I studied groundwater and surface water interactions in the Kenai Lowlands, Alaska. In particular, I examine two important aspects of groundwater and surface water interactions: 1) Groundwater's influence on surface-water temperature; and 2) Groundwater's role in forming hydrologic flow paths that can connect uplands to streamside wetlands and streams. Chapter 2 investigates the controls on stream temperature in salmon-bearing headwater streams in two common hydrogeologic settings: 1) drainage-ways, which are low-gradient streams that flow through broad valleys; and 2) discharge-slopes, which are high gradient streams that flow through narrow valleys. The results from chapter 2 showed significant differences in stream temperatures between the two hydrogeologic settings. Observed stream temperatures were higher in drainage-way sites than in discharge-slope sites, and showed strong correlations as a continuous function with the calculated topographic metric flow-weighted slope. Additionally, modeling results indicated that the potential for groundwater discharge to moderate stream temperature is not equal between the two hydrogeologic settings, with groundwater having a greater moderating effect on stream temperature at the low gradient drainage-way sites. Chapter 3 examines the influence of groundwater on juvenile coho salmon winter habitat along the Anchor River. Two backwater habitats were selected from the larger set of 25 coho overwintering sites from a previous study for an in-depth hydrologic analysis. The results from chapter 3 showed that the type of groundwater discharge (i.e., focused versus diffuse groundwater discharge) can play an important role in determining habitat suitability in these backwater sites. During winter, focused discharge from a local groundwater seep maintained higher surface-water temperatures and higher concentrations of dissolved oxygen compared to the site with more diffuse groundwater discharge. Chapter 4 investigates the linkages along hydrologic flow paths among alder (Alnus spp.) stands, streamside wetlands, and headwater streams. Chapter 4 tested four related hypotheses: 1) groundwater nitrate concentrations are greater along flow paths with alder compared to flow paths without alder; 2) on hillslopes with alder, groundwater nitrate concentrations are highest when alder stands are located near the streamside wetlands at the base of the hillslope; 3) primary production of streamside wetland vegetation is N limited and wetlands are less N limited when alder stands are located nearby along flow paths; and 4) stream reaches at the base of flow paths with alder have higher nitrate concentrations than reaches at the base of flow paths without alder. The results from chapter 4 showed that groundwater nitrate concentrations were highest along flow paths with alder, however no difference was observed between flow paths with alder located near versus alder located further from streamside wetlands. Vegetation had a greater response to N fertilization in streamside wetlands that were connected to flow paths without alder and less when alder stands were near. Finally, higher nitrate concentrations were measured in streams at the base of flow paths with alder. The combined results of this dissertation showed that, in the Kenai Lowlands, groundwater and surface water interactions have a direct influence on the local ecology and that a fundamental understanding of the hydrology can aid in the successful management and protection of this unique and important ecosystem.

A major challenge to cold regions hydrology and northern water resources management lies in predicting runoff dynamically in the context of warming-induced changes to the rates and patterns of ground thaw and drainage. Meeting this challenge requires new knowledge of the mechanisms and rates of ground thaw and their implications to water drainage and storage patterns and processes. The study carries out to evaluate the concept of energy-based runoff in the perspective of ground heat flux, soil thaw and liquid moisture content, tortuosity of snow-free area, preferential flow and discharge of the hillslope. Based on field measurements, coupled energy and water flow is simulated in the Area of Interest (AOI) with a half-hour time interval by the distributed hydrological model, GEOtop. In the field, the saturated hydraulic conductivity varies exponentially between the superficial organic layer and the underlying mineral layer. In the simulation, the parameters of the soil physical properties are input by fourteen uneven layers below the ground surface. Starting from the initially frozen state, the process of soil thaw is simulated with dynamic variables such as soil liquid moisture and ice content, hydraulic conductivity, thermal conductivity and heat capacity. The simulated frost table depths are validated by 44-point measurements and the simulation of point soil temperature is also compared to data measured in an excavated soil pit. As a result, the frost table topography is dominated by both the snow-free pattern and the energy fluxes on the ground surface. The rate and magnitude of runoff derived from snow drift and the ice content of frozen soil is greatly influenced by the frost table topography. According to the simulation, the frost table depth is closely regressed with the ground surface temperature by a power function. As soil thawing progresses, ground heat flux reduces gradually and the rate of soil thaw becomes small when the frost table descends. Along with the snow-free area expanding, the average soil moisture of the AOI increases prior to that time when the average frost table is less than 25 cm deep. The snow-free patches expand heterogeneously in the AOI, which causes the spatial and temporal variation of hydraulic conductivity due to the non-uniform frost table depth. According to the simulation, the transit time of the flow through the AOI decreases to the shortest span on May 13 with the average frost table of 10 cm. Before this date, the time lag between snowmelt percolation and slope runoff is about 8-10 hours; while after this date, the time lag is no more than 5 hours. The pattern of the preferential flow in the AOI highly depends on the frost table topography. When the snow-free patches are widely scattered and the average frost table is between 0 and 10 cm, the preferential flow paths are inhibited. With soil thaw progresses, the preferential flow paths are prominent with the largest single contributing area occurring when the average frost table is between 10 cm to 15 cm. When the average frost table reaches 25 cm, the importance of preferential flow is not apparent, and matrix flow prevails.

In this study, I examined the effect of watershed urbanization on the invasion of alien woody species in riparian forests. This study was conducted in three major steps: 1) estimating the degree of watershed urbanization using impervious surface maps extracted from remote sensing images; 2) examining the effect of urbanization on hydrologic regime; and 3) investigating a relationship between watershed urbanization and ecosystem invasibility of a riparian forest. I studied twelve riparian forests along urban-rural gradients in Austin, Texas. Hydrologic regimes were quantified by transfer function (TF) models using four-year daily rainfall-streamflow data in two study periods (10/1988-09/1992 and 10/2004-09/2008) between which Austin had experienced rapid urbanization. For each study period, an impervious surface map was generated from Landsat TM image by a support vector machine (SVM) with pairwise coupling. SVM more accurately estimated impervious surface than other subpixel mapping methods. Ecosystem invasibilities were assessed by relative alien cover (RAC) of riparian woody species communities. The results showed that the effects of urbanization differ by hydrogeologic conditions. Of the study watersheds, seven located in a hillslope region experienced the diminishing peakflows between the two study periods, which are contrary to current urban hydrologic model. I attributed the decreased peakflows to land grading that transformed a hillslope into a stair-stepped landscape. In the rest of the watersheds, peakflow diminished between the two study periods perhaps due to the decrease in stormwater infiltration and groundwater pumpage that lowered groundwater level. In both types of watersheds, streamflow rising during a storm event more quickly receded as watershed became more urbanized. This study found a positive relationship between RAC and watershed impervious surface percentage. RAC was also significantly related to flow recession and canopy gap percentages, both of which are indicators of hydrologic disturbance. These results suggest that urbanization facilitated the invasion of alien species in riparian forests by intensifying hydrologic disturbance. The effects of urbanization on ecosystems are complex and vary by local hydrologeologic conditions. These results imply that protection of urban ecosystems should be based on a comprehensive and large-scale management plan.

Subsurface solutional pathways make limestone terrains sensitive to changes in soil properties that regulate flows to the epikarst. This study examines biogeomorphic factors responsible for changed water movements and erosion in fluviokarst slopes deforested 200 years ago along the Kentucky River, Kentucky. In this project, infiltration and water content data from forest and fescue grass soil profiles were analyzed within a detailed overview of system factors regulating hillslope hydrology. Results show that grass has growth and rooting characteristics that tend to create a larger volume of lateral water movement in upper soil layers than occurs under forests. This sets up the current emergent pattern of erosion in which water perches at grass slope bases and overwhelms pre-existing epikarst drainage. Tree roots are able to cause solution at multiple discrete points of entry into fractures and bedding planes, increasing storage capacity and releasing sediment over time. Grass roots do not enter bedrock, and their rooting depth limits diffuse vertical preferential flow in root channels to above one meter. In the areas dense clay soils, flow under grass is conducted sideways either through the regolith or at the bedrock surface. Rapid flow along rock faces in hillslope benches likely moves fines via subsurface routes from the hillslope shoulders, causing the exposure of flat outcrops under grass. Lower growing season evapotranspiration also promotes higher grass summer flow volumes. Gullying occurs at sensitive points where cutters pass from the uphill grassed area into the forest, or where flow across the bedrock surface crosses grass/forest boundaries oriented vertical to the slope. At these locations, loss of the protective grass root mat, coupled with instigation of tree root preferential flow in saturated soils, causes soil pipes to develop. Fluviokarst land management decisions should be based on site-specific slope, soil depth, and epkarst drainage conditions, since zones sensitive to erosion are formed by spatial and temporal conjunctions of a large number of lithologic, karst, soil, climate, and vegetation factors. This study shows that it is the composite of differing influences created by forest and grass that make forests critical for soil retention in high-energy limestone terrains.

A major challenge to cold regions hydrology and northern water resources management lies in predicting runoff dynamically in the context of warming-induced changes to the rates and patterns of ground thaw and drainage. Meeting this challenge requires new knowledge of the mechanisms and rates of ground thaw and their implications to water drainage and storage patterns and processes. The study carries out to evaluate the concept of energy-based runoff in the perspective of ground heat flux, soil thaw and liquid moisture content, tortuosity of snow-free area, preferential flow and discharge of the hillslope. Based on field measurements, coupled energy and water flow is simulated in the Area of Interest (AOI) with a half-hour time interval by the distributed hydrological model, GEOtop. In the field, the saturated hydraulic conductivity varies exponentially between the superficial organic layer and the underlying mineral layer. In the simulation, the parameters of the soil physical properties are input by fourteen uneven layers below the ground surface. Starting from the initially frozen state, the process of soil thaw is simulated with dynamic variables such as soil liquid moisture and ice content, hydraulic conductivity, thermal conductivity and heat capacity. The simulated frost table depths are validated by 44-point measurements and the simulation of point soil temperature is also compared to data measured in an excavated soil pit. As a result, the frost table topography is dominated by both the snow-free pattern and the energy fluxes on the ground surface. The rate and magnitude of runoff derived from snow drift and the ice content of frozen soil is greatly influenced by the frost table topography. According to the simulation, the frost table depth is closely regressed with the ground surface temperature by a power function. As soil thawing progresses, ground heat flux reduces gradually and the rate of soil thaw becomes small when the frost table descends. Along with the snow-free area expanding, the average soil moisture of the AOI increases prior to that time when the average frost table is less than 25 cm deep. The snow-free patches expand heterogeneously in the AOI, which causes the spatial and temporal variation of hydraulic conductivity due to the non-uniform frost table depth. According to the simulation, the transit time of the flow through the AOI decreases to the shortest span on May 13 with the average frost table of 10 cm. Before this date, the time lag between snowmelt percolation and slope runoff is about 8-10 hours; while after this date, the time lag is no more than 5 hours. The pattern of the preferential flow in the AOI highly depends on the frost table topography. When the snow-free patches are widely scattered and the average frost table is between 0 and 10 cm, the preferential flow paths are inhibited. With soil thaw progresses, the preferential flow paths are prominent with the largest single contributing area occurring when the average frost table is between 10 cm to 15 cm. When the average frost table reaches 25 cm, the importance of preferential flow is not apparent, and matrix flow prevails.