Browsing by Author "Pollyea, Ryan M."
Now showing 1 - 15 of 15
Results Per Page
Sort Options
- Advancing the Global Land Grant Institution: Creating a Virtual Environment to Re-envision Extension and Advance GSS-related Research, Education, and CollaborationHall, Ralph P.; Polys, Nicholas F.; Sforza, Peter M.; Eubank, Stephen D.; Lewis, Bryan L.; Krometis, Leigh-Anne H.; Pollyea, Ryan M.; Schoenholtz, Stephen H.; Sridhar, Venkataramana; Crowder, Van; Lipsey, John; Christie, Maria Elisa; Glasson, George E.; Scherer, Hannah H.; Davis, A. Jack; Dunay, Robert J.; King, Nathan T.; Muelenaer, Andre A.; Muelenaer, Penelope; Rist, Cassidy; Wenzel, Sophie (Virginia Tech, 2017-05-15)The vision for this project has emerged from several years of research, teaching, and service in Africa and holds the potential to internationalize education at Virginia Tech and in our partner institutions in Malawi. The vision is simple, to develop a state-of-the-art, data rich, virtual decision-support and learning environment that enables local-, regional-, and national-level actors in developed and developing regions to make decisions that improve resilience and sustainability. Achieving these objectives will require a system that can combine biogeophysical and sociocultural data in a way that enables actors to understand and leverage these data to enhance decision-making at various levels. The project will begin by focusing on water, agricultural, and health systems in Malawi, and can be expanded over time to include any sector or system in any country. The core ideas are inherently scalable...
- Advective Heat Transport and the Salt Chimney Effect: A Numerical AnalysisCanova, David P.; Fischer, Mark P.; Jayne, Richard S.; Pollyea, Ryan M. (Hindawi, 2018-07-09)We conducted numerical simulations of coupled fluid and heat transport in an offshore, buried salt diapir environment to determine the effects of advective heat transport and its relation to the so-called “salt chimney effect.” Model sets were designed to investigate (1) salt geometry, (2) depth-dependent permeability, (3) geologic heterogeneity, and (4) the relative influence of each of these factors. Results show that decreasing the dip of the diapir induces advective heat transfer up the side of the diapir, elevating temperatures in the basin. Depth-dependent permeability causes upwelling of warm waters in the basin, which we show to be more sensitive to basal heat flux than brine concentration. In these model scenarios, heat is advected up the side of the diapir in a narrower zone of upward-flowing warm water, while cool waters away from the diapir flank circulate deeper into the basin. The resulting fluid circulation pattern causes increased discharge at the diapir margin and fluid flow downward, above the crest of the diapir. Geologic heterogeneity decreases the overall effects of advective heat transfer. The presence of low permeability sealing horizons reduces the vertical extent of convection cells, and fluid flow is dominantly up the diapir flank. The combined effects of depth-dependent permeability coupled with geologic heterogeneity simulate several geologic phenomena that are reported in the literature. In this model scenario, conductive heat transfer dominates in the basal units, whereas advection of heat begins to affect the middle layers of the model and dominates the upper units. Convection cells split by sealing layers develop within the upper units. From our highly simplified models, we can predict that advective heat transport (i.e., thermal convection) likely dominates in the early phases of diapirism when sediments have not undergone significant compaction and retain high porosity and permeability. As the salt structures mature into more complex geometries, advection will diminish due to the increase in dip of the salt-sediment interface and the increased hydraulic heterogeneity due to complex stratigraphic architecture.
- Analyzing a 10-Year Cave Drip Record in James Cave, Virginia: Implications for Storage and Recharge in Shallow Appalachian Karst SystemsGroce-Wright, Nigel C. (Virginia Tech, 2021-07-16)Karst aquifers, characterized by soluble rocks such as limestone and dolostone, provide drinking water to 20-25% of the world's population and are thus critical global water sources. However, recent work suggests that rapid alteration of karst aquifers due to the impact of climate change on precipitation patterns may affect recharge to these aquifers. Much of the research on recharge in karst aquifers has relied on using patterns of spring discharge to quantify recharge. Spring outlets allow for continuous monitoring of discharge from karst aquifers, making them easily accessible monitoring sites. However, because springs can integrate multiple flow paths, it is difficult to rely on spring discharge patterns to get information on where and how karst aquifers are receiving recharge. Monitoring closer to the source of recharge through the measurement of cave drips allows for a more accurate analysis of recharge timing and mechanisms. In this study, I conducted recession analyses on cave drip hydrographs from a 10-year record (2008-2018) of three drip monitoring stations within James Cave (Pulaski Co., VA) to: 1) examine differences in hydrologic characteristics of the epikarst (the zone of soil and weathered bedrock above a karst aquifer); 2) quantify the storage volume of the epikarst and 3) investigate seasonal, and annual trends in recharge. Results of recession analysis show heterogeneity in epikarst hydrologic characteristics, reflected by calculations of the recession coefficient, , and storage volume. Calculations of the recession coefficient show subtle differences between the three drip sites, suggestive of spatial heterogeneity in permeability and storage in the overlying epikarst. The storage volume calculations show that during the recharge season (winter- spring), up to 95% of recharge through the unsaturated zone to the cave occurs through rapid pathways (i.e., fractures), and 5% through diffuse pathways (i.e., pores). However, during the recession period (spring-summer), when evapotranspiration is active, recharge through cave drips decreases and occurs predominantly through diffuse flow. Combined, these results underscore the importance of both spatial and temporal characterization of drip rates and other recharge inputs into karst aquifer systems.
- High density oilfield wastewater disposal causes deeper, stronger, and more persistent earthquakesPollyea, Ryan M.; Chapman, Martin C.; Jayne, Richard S.; Wu, Hao (Nature, 2019)Oilfield wastewater disposal causes fluid pressure transients that induce earthquakes. Here we show that, in addition to pressure transients related to pumping, there are pressure transients caused by density differences between the wastewater and host rock fluids. In northern Oklahoma, this effect caused earthquakes to migrate downward at ~0.5 km per year during a period of high-rate injections. Following substantial injection rate reductions, the downward earthquake migration rate slowed to ~0.1 km per year. Our model of this scenario shows that the density-driven pressure front migrates downward at comparable rates. This effect may locally increase fluid pressure below injection wells for 10+ years after substantial injection rate reductions. We also show that in north-central Oklahoma the relative proportion of high-magnitude earthquakes increases at 8+ km depth. Thus, our study implies that, following injection rate reductions, the frequency of high-magnitude earthquakes may decay more slowly than the overall earthquake rate.
- Hydrothermal Transport in the Panama Basin and in Brothers Volcano using Heat Flow, Scientific Deep Sea Drilling and Mathematical ModelsKolandaivelu, Kannikha Parameswari (Virginia Tech, 2019-02-15)Two-thirds of submarine volcanism in the Earth's ocean basins is manifested along mid-ocean ridges and the remaining one-third is revealed along intraoceanic arcs and seamounts. Hydrothermal systems and the circulation patterns associated with these volcanic settings remove heat from the solid Earth into the deep ocean. Hydrothermal circulation continues to remove and redistribute heat in the crust as it ages. The heat and mass fluxes added to the deep ocean influence mixing in the abyssal ocean thereby affecting global thermohaline circulation. In addition to removing heat, hydrothermal processes extract chemical components from the oceanic and carry it to the surface of the ocean floor, while also removing certain elements from seawater. The resulting geochemical cycling has ramifications on the localized mineral deposits and also the biota that utilize these chemical fluxes as nutrients. In this dissertation, I analyze observed conductive heat flow measurements in the Panama Basin and borehole thermal measurements in Brothers Volcano and use mathematical models to estimate advective heat and mass fluxes, and crustal permeability. In the first manuscript, I use a well-mixed aquifer model to explain the heat transport in a sediment pond in the inactive part of the Ecuador Fracture Zone. This model yields mass fluxes and permeabilities similar to estimates at young upper oceanic crust suggesting vigorous convection beneath the sediment layer. In the second manuscript, I analyze the conductive heat flow measurements made in oceanic between 1.5 and 5.7 Ma on the southern flank of the Costa Rica Rift. These data show a mean conductive heat deficit of 70%, and this deficit is explained by various hydrothermal advective transport mechanisms, including outcrop to outcrop circulation, transport through faults, and redistribution of heat by flow of hydrothermal fluids in the basement. In the third manuscript, I analyze the borehole temperature logs for two sites representative of recharge and discharge areas of hydrothermal systems in the Brothers Volcano. I develop upflow and downflow models for fluids in the borehole and formation resulting in estimated of flow rates and permeabilities. All three independent research works are connected by the common thread of utilizing relatively simple mathematical concepts to get new insights into hydrothermal processes in oceanic crust.
- Identification of Recharge Source Areas in a Fractured Crystalline-rock Aquifer in Ploemeur, FranceHumm, Cathleen Hana (Virginia Tech, 2021-06-17)Characterizing and preserving available groundwater resources within crystalline rocks is pertinent to understanding and predicting resources for ecosystems worldwide. Crystalline-rock aquifers, with favorable structure and climate, can be pumped year-round to meet local domestic demand. The Ploemeur hydrogeologic site, near the southern coast of Brittany, France, is characterized by a structurally complex fractured mica-schist and granite confined aquifer system. A contact zone, which acts as the main localized flow path through the aquifer, separates the two crystalline units, and a sub-vertical permeable fault zone cross-cuts the crystalline bedrock and contact zone. Using field observations, recharge estimates, and a calibrated three-dimensional numerical multi-zone MODFLOW 6 model, we present preferential flow paths of recharge infiltrating the complex geology of the Ploemeur hydrogeological site during pumping conditions. Using MODPATH to track groundwater and recharge path lines, we determine that water extracted from the aquifer originates from higher elevation areas west of the pumping site. Particle tracking analyses indicate that precipitation simulated over the pumping zone takes a minimum of two years to reach the pumping wells and travels up to 100 m in distance. Analyses of the water budget of the aquifer system using Zonebudget show that storage contributes significantly to the productivity of the system. Based on these analyses, we determine that recharge mechanisms such as piston flow and preferential flow play important roles in the Ploemeur hydrogeologic site. Though the Ploemeur site is unique in its composition and geometry, the methods used to characterize and monitor the aquifer can be applied to fractured crystalline-rock aquifers globally. Fractured crystalline-rock aquifers make up 10% of the region's freshwater sources, thus understanding their flow mechanisms contributes greatly to the management of freshwater resources.
- Implications of Permeability Uncertainty During Three-phase CO2 Flow in a Basalt Fracture NetworkGierzynski, Alec Owen (Virginia Tech, 2016-12-15)Recent studies suggest that continental flood basalts may be suitable for geologic carbon sequestration due to fluid-rock reactions that mineralize injected CO₂ on relatively short time-scales. Flood basalts also possess a permeability structure favorable for injection, with alternating high-permeability (flow margin) and low-permeability (flow interior) layers. However, little information exists on the behavior of CO₂ as it leaks through fractures characteristic of the flow interior, particularly at conditions near the critical point for CO₂. In this study, a two-dimensional 5 × 5 m model of a fracture network is built based on high-resolution LiDAR scans of a Columbia River Basalt flow interior taken near Starbuck, WA. Three-phase CO₂ flow is simulated using TOUGH3 (beta) with equation of state ECO2M for 10 years simulation time. Initial conditions comprise a hydrostatic pressure profile corresponding to 750-755 m below ground surface and a constant temperature of 32° C. Under these conditions, the critical point for CO₂ occurs 1.5 meters above the bottom of the domain. Matrix permeability is assumed to be constant, based on literature values for the Columbia River Basalt. Fracture permeability is assigned based on a lognormal distribution of random values with mean and standard deviation based on measured fracture aperture values and in situ permeability values from literature. In order to account for fracture permeability uncertainty, CO₂ leakage is simulated in 50 equally probable realizations of the same fracture network with spatially random permeability constrained by the lognormal permeability distribution. Results suggest that fracture permeability uncertainty has some effect on the distribution of CO₂ within the fractures, but network geometry is the primary control in determining flow paths. Fracture permeability uncertainty has a larger influence on fluid pressure, and can affect the location of the critical point within ~1.5 m. Uncertainty in fluid pressure was found to be highest along major flow paths below channel constrictions, indicating permeability at a few key points can have a large influence on fluid pressure distribution.
- Implications of permeability uncertainty within engineered geologic fluid systemsJayne Jr, Richard Scott (Virginia Tech, 2019-10-07)Carbon-capture and sequestration (CCS) in geologic reservoirs is one strategy for reducing anthropogenic CO2 emissions from large-scale point source emitters. Recent developments have shown that basalt reservoirs are highly effective for permanent mineral trapping on the basis of CO2-water-rock interactions, which result in the formation of carbonate minerals. However, the injection of super-critical CO2 into the subsurface causes a disturbance in the pressure, temperature, and chemical systems within the target reservoir. How the ambient conditions change in response to a CO2 injection ultimately affects the transport and fate of the injected CO2. Understanding the behavior and transport of CO2 within a geologic reservoir is a difficult problem that is only exacerbated by heterogeneities within the reservoir; for example, permeability can be highly heterogeneous and exhibits significant control on the movement of CO2. This work is focused on constraining the permeability uncertainty within a flood basalt reservoir, specifically the Columbia River Basalt Group (CRBG). In order to do so, this dissertation is a culmination of four projects: (1) a geostatistical analysis resulting in a spatial correlation model of regional scale permeability within the CRBG, (2) a Monte Carlo-type modeling studying investigating the effects that permeability uncertainty has on the injectivity and storativity of the CRBG as a storage reservoir, (3) a modeling study utilizing 1-, 2-, and 3-D numerical models to investigate how the thermal signature of the CO2-water system evolves during a CO2 injection, and (4) a Monte Carlo-type modeling study focused on the integrity of the CRBG as a CO2 storage reservoir through a probabilistic assessment of static threshold criteria.
- A Numerical and Statistical Analysis of the Fractured Rock Aquifer System in Ploemeur, France to Quantify Local and Regional RechargeLaw, Stacey E. (Virginia Tech, 2019-08-14)Groundwater recharge is an essential metric for understanding and protecting groundwater resources. Quantifying this parameter remains extremely challenging due to the uncertainties associated with the extent to which the vadose zone affects groundwater movement and the highly heterogeneous nature of the aquifer systems being monitored. The difficulty surrounding recharge quantification is compounded when considering a fractured rock aquifer system, where classification and modeling is complicated by highly complex structural geology. However, the ability to distinguish the character and geometry of fractured rock aquifers is indispensable for quantifying recharge to evaluate sustainable yields, as well as for implementing protective measures to manage these systems. The primary intention of this study is to assess the hydrogeologic properties that have led the unique recharge signals within the fractured crystalline-rock aquifer system near Ploemeur, France. Infiltration and groundwater movement are characterized via time-series hydraulic head and precipitation data collected at daily, monthly, yearly, and at decadal intervals. In spite of the nearly one million cubic meters of groundwater extraction, measured drawdowns are marginal, suggesting that local and regional recharge plays a significant role in moderating water-level declines and raising questions as to the origins of the substantial inflow required to sustain this complex system. A roughly two-month lag has been observed between seasonal water level and monthly precipitation at Ploemeur, which has previously been attributed solely to slow vertical migration of water through the low-permeability micaschist layer to the fractured contact zone and interconnected fault. However, results from this study suggest that a significant portion of the observed lag can be attributed to vadose-zone processes, particularly the thickness of the vadose zone. This investigation also reveals a recharge signal that continues throughout the calendar year, departing from the traditional simplified concept that recharge quantity is essentially equivalent to the value of evapotranspiration subtracted from infiltration.
- Numerical Investigations of Geologic CO2 Sequestration Using Physics-Based and Machine Learning Modeling StrategiesWu, Hao (Virginia Tech, 2020-08-06)Carbon capture and sequestration (CCS) is an engineering-based approach for mitigating excess anthropogenic CO2 emissions. Deep brine aquifers and basalt reservoirs have shown outstanding performance in CO2 storage based on their global widespread distribution and large storage capacity. Capillary trapping and mineral trapping are the two dominant mechanisms controlling the distribution, migration, and transportation of CO2 in deep brine aquifers and basalt reservoirs. Understanding the behavior of CO2 in a storage reservoir under realistic conditions is important for risk management and storage efficiency improvement. As a result, numerical simulations have been implemented to understand the relationship between fluid properties and multi-phase fluid dynamics. However, the physics-based simulations that focus on the uncertainties of fluid flow dynamics are complicated and computationally expensive. Machine learning method provides immense potential for improving computational efficiency for subsurface simulations, particularly in the context of parametric sensitivity. This work focuses on parametric uncertainty associated with multi-phase fluid dynamics that govern geologic CO2 storage. The effects of this uncertainty are interrogated through ensemble simulation methods that implement both physics-based and machine learning modeling strategies. This dissertation is a culmination of three projects: (1) a parametric analysis of capillary pressure variability effects on CO2 migration, (2) a reactive transport simulation in a basalt fracture system investigating the effects of carbon mineralization on CO2 migration, and (3) a parametric analysis based on machine learning methods of simultaneous effects of capillary pressure and relative permeability on CO2 migration.
- Pore-scale Interfacial and Transport Phenomena in Hydrocarbon ReservoirsFang, Chao (Virginia Tech, 2019-06-10)Exploring unconventional hydrocarbon reservoirs and enhancing the recovery of hydrocarbon from conventional reservoirs are necessary for meeting the society's ever-increasing energy demand and requires a thorough understanding of the multiphase interfacial and transport phenomena in these reservoirs. This dissertation performs pore-scale studies of interfacial thermodynamics and multiphase hydrodynamics in shale reservoirs and conventional oil-brine-rock (OBR) systems. In shale gas reservoirs, the imbibition of water through surface hydration into gas-filled mica pores was found to follow the diffusive scaling law, but with an effective diffusivity much larger than the self-diffusivity of water molecules. The invasion of gas into water-filled pores with width down to 2nm occurs at a critical invasion pressure similar to that predicted by the classical capillary theories if effects of disjoining pressure and diffusiveness of water-gas interfaces are considered. The invasion of oil droplets into water-filled pores can face a free energy barrier if the pressure difference along pore is small. The computed free energy profiles are quantitatively captured by continuum theories if capillary and disjoining pressure effects are considered. Small droplets can invade a pore through thermal activation even if an energy barrier exists for its invasion. In conventional oil reservoirs, low-salinity waterflooding is an enhanced oil recovery method that relies on the modification of thin brine films in OBR systems by salinity change. A systematic study of the structure, disjoining pressure, and dynamic properties of these thin brine films was performed. As brine films are squeezed down to sub-nanometer scale, the structure of water-rock and water-oil interfaces changes marginally, but that of the electrical double layers in the films changes greatly. The disjoining pressure in the film and its response to salinity change follow the trend predicted by the DLVO theory, although the hydration and double layer forces are not simple additive as commonly assumed. A notable slip between the brine film and the oil phase can occur. The role of thin liquid films in multiphase transport in hydrocarbon reservoirs revealed here helps lay foundation for manipulating and leveraging these films to enhance hydrocarbon production and to minimize environmental damage during such extraction.
- PVTX and Raman Spectral Properties of Fluids at Elevated Pressures and TemperaturesSublett, David Matthew Jr. (Virginia Tech, 2020-01-08)Fluids are associated with a wide range of physical and chemical processes in the Earth, including transporting and concentrating important ore elements such as Cu, Au, Zn, and Pb. Significant amounts of fluid may be generated as a result of dehydration or decarbonation reactions, and the volatile content of a magma is directly linked to the explosivity of eruptions. In most cases, small amounts of the fluids involved in the formation or alteration of rocks are trapped within minerals in the form of fluid inclusions. These fluid inclusions may be studied to understand the composition and pressure and temperature of the original fluid involved in the geologic process of interest, however, an understanding of the composition of the fluid as well as how the fluid behaves under changing pressure and temperature conditions is essential to reconstruct the fluid evolution path based on data obtained from fluid inclusions. Several analytical techniques are involved in the study of fluids, including fluid inclusion microthermometry and Raman spectroscopy. Microthermometry is the heating/cooling of fluid inclusions to observe and record temperatures of phase changes which, in turn, are used to determine properties such as salinity (based on the freezing point depression of liquid), or density based on the temperature at which all phases within the fluid inclusion homogenize to a single phase. Raman spectroscopy is a non-destructive analytical technique that measures the vibrational frequency of molecules in a given material. The Raman spectral properties of fluids act as a "fingerprint" of the chemical species within the fluid and serve to identify both the presence of chemical species, such as H2O, N2, CO2, and CH4, and the density of the fluid. Microthermometric and Raman spectroscopic experiments involving synthetic fluid systems are necessary to elucidate the pressure-volume-temperature-composition (PVTX) and Raman spectral behavior of the fluid systems, which then aids in the study and characterization of natural fluids. In chapter 1, the partitioning of NaCl and KCl between coexisting immiscible fluid phases during boiling is experimentally determined at temperatures and pressures relevant to magmatic-hydrothermal systems using synthetic fluid inclusions. The partitioning behavior is then combined with literature data to calculate the Na/K ratio of the original silicate melt phase in a magma body before the exsolution of a fluid phase. In chapter 2, we explore the Raman spectral behavior of N2, CO2, and CH4 in pure, single-component systems from PT conditions corresponding to the liquid-vapor curve to elevated temperatures and pressures, and relate the changes in the spectral behavior to changes in the bonding environment of the molecules through intermolecular attraction and repulsion. In chapter 3, the observations and relationships determined for pure fluids and described in chapter 2 are used to explore the Raman spectral properties of N2, CO2, and CH4 in the N2-CO2-CH4 ternary system and the manner in which the spectral behavior of each component in the system varies with changing temperature, pressure, molar volume, and fugacity.
- Regional-Scale Impacts of Fluid Composition and Geologic Structure for Injection-Induced Seismicity in the Southern U.S. MidcontinentKonzen, Graydon Leo (Virginia Tech, 2020-06-24)Over the last decade, an increase in earthquake occurrence in Oklahoma and Kansas has been linked to oilfield wastewater injection disposal, particularly into the regionally underpressured Arbuckle Group. The Arbuckle is hydraulically connected to Precambrian basement through an extensive fracture system, which transmits pressure perturbations from wastewater injections to seismogenic depths. Previous studies have convincingly attributed induced seismicity to pore pressure diffusion and solid elastic stressing, both resulting from fluid waste injection. Recent work adds to the physical understanding of injection-induced seismicity by demonstrating that the density differential between injection fluids and formation brines may also drive fluid pressure into the seismogenic basement. In this thesis, variable density groundwater flow is modeled in a numerical simulation comprising parts of the Anadarko Basin, the Anadarko Shelf, the Cherokee Platform, and the Nemaha Fault Zone as well as injection data from 2006-2018. Results show buoyancy forces interacting with regional stratigraphic dip to force density-driven pressure transients into the deep Anadarko Basin, aligning with previously unexplained earthquakes in that region.
- The Role of Chemistry in Fracture Pattern Development and Opportunities to Advance Interpretations of Geological MaterialsLaubach, S. E.; Lander, R. H.; Criscenti, L. J.; Anovitz, L. M.; Urai, J. L.; Pollyea, Ryan M.; Hooker, J. N.; Narr, W.; Evans, M. A.; Kerisit, S. N.; Olson, J. E.; Dewers, T.; Fisher, D.; Bodnar, Robert J.; Evans, B.; Dove, Patricia M.; Bonnell, L. M.; Marder, M. P.; Pyrak-Nolte, L. (AGU, 2019-09-13)Fracture pattern development has been a challenging area of research in the Earth sciences for more than 100 years. Much has been learned about the spatial and temporal complexity inherent to these systems, but severe challenges remain. Future advances will require new approaches. Chemical processes play a larger role in opening‐mode fracture pattern development than has hitherto been appreciated. This review examines relationships between mechanical and geochemical processes that influence the fracture patterns recorded in natural settings. For fractures formed in diagenetic settings (~50 to 200 °C), we review evidence of chemical reactions in fractures and show how a chemical perspective helps solve problems in fracture analysis. We also outline impediments to subsurface pattern measurement and interpretation, assess implications of discoveries in fracture history reconstruction for process‐based models, review models of fracture cementation and chemically assisted fracture growth, and discuss promising paths for future work. To accurately predict the mechanical and fluid flow properties of fracture systems, a processes‐based approach is needed. Progress is possible using observational, experimental, and modeling approaches that view fracture patterns and properties as the result of coupled mechanical and chemical processes. A critical area is reconstructing patterns through time. Such data sets are essential for developing and testing predictive models. Other topics that need work include models of crystal growth and dissolution rates under geological conditions, cement mechanical effects, and subcritical crack propagation. Advances in machine learning and 3‐D imaging present opportunities for a mechanistic understanding of fracture formation and development, enabling prediction of spatial and temporal complexity over geologic timescales. Geophysical research with a chemical perspective is needed to correctly identify and interpret fractures from geophysical measurements during site characterization and monitoring of subsurface engineering activities.
- Temporal Dynamics of Groundwater Flow Direction in a Glaciated, Headwater CatchmentBenton, Joshua Robert (Virginia Tech, 2020-05-12)Shallow groundwater flow in the critical zone of steep headwater mountain catchments is often assumed to mimic surface topography. However, groundwater flow is influenced by other variables, such as the elevation of the water table and subsurface hydraulic conductivity, which can result in temporal variations in both magnitude and direction of flow. In this study, I investigated the temporal variability of groundwater flow in the soil zone (solum) within the critical zone of a headwater catchment at the Hubbard Brook Experimental Forest in North Woodstock, NH. Groundwater levels were continuously monitored throughout several seasons (March 2019 to Jan 2020) in a network of wells comprising three hillslope transects within the upper hillslopes of the catchment. Five clusters of three wells per cluster were screened from 0.18 – 1.1 m depth at the base of the solum. Water levels were also monitored in five deeper wells, screened from 2.4 - 6.9 m depth within glacial sediments of the C horizon. I conducted 47 slug tests across the well network to determine hydraulic properties of the aquifer materials surrounding each well. In addition, our team conducted a large-scale auger investigation mapping soil horizon depths and thicknesses. Results show that the magnitude of hydraulic gradients and subsurface hydrologic fluxes varied at each site with respect to changing water-table elevation, having a maximum range of 0.12 m/m and 9.19 x 10-6 m/s, respectively. The direction of groundwater flow had an arithmetic mean deviating from surface topography by 2-10 degrees, and a total range that deviated from surface topography by as much as 51 degrees. During lower water table regimes, groundwater flow direction deviated from the ground surface, but under higher water table regimes, in response to recharge events, flow direction mimicked surface topography. Within most of the well clusters, there is an observable connection between the slope direction of the top of the C horizon and the direction of groundwater flow during lower water table regimes. Slug test results show the interquartile range of saturated hydraulic conductivity (Ksat) within the C horizon (1.5×10-7 to 9.8×10-7 m/s) is two orders of magnitude lower than the interquartile range of Ksat values within the solum (2.9×10-5 to 5.2×10-5 m/s). Thus, the C horizon is on average a confining unit relative to the solum that may constrict groundwater flow below the solum. Additionally, results from the larger scale auger investigation suggest that deviations in subsurface topography of the C horizon may be generalizable at the larger hillslope scale. Overall, these results indicate that 1) shallow groundwater flow direction and magnitude within this headwater catchment are dynamic and can deviate from surface topography, and 2) the subsurface topography of the C horizon can influence groundwater flow direction. These results imply that temporal dynamics of groundwater flow direction and magnitude should be considered when characterizing subsurface flow in critical zone studies. Additionally, knowledge of subsurface topography of confining units may provide constraints on the temporal variability of groundwater flow direction.