skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, February 13 until 2:00 AM ET on Friday, February 14 due to maintenance. We apologize for the inconvenience.


Title: Active Layer Groundwater Flow: The Interrelated Effects of Stratigraphy, Thaw, and Topography
Abstract

The external drivers and internal controls of groundwater flow in the thawed “active layer” above permafrost are poorly constrained because they are dynamic and spatially variable. Understanding these controls is critical because groundwater can supply solutes such as dissolved organic matter to surface water bodies. We calculated steady‐state three‐dimensional suprapermafrost groundwater flow through the active layer using measurements of aquifer geometry, saturated thickness, and hydraulic properties collected from two major landscape types over time within a first‐order Arctic watershed. The depth position and thickness of the saturated zone is the dominant control of groundwater flow variability between sites and during different times of year. The effect of water table depth on groundwater flow dwarfs the effect of thaw depth. In landscapes with low land‐surface slopes (2–4%), a combination of higher water tables and thicker, permeable peat deposits cause relatively constant groundwater flows between the early and late thawed seasons. Landscapes with larger land‐surface slopes (4–10%) have both deeper water tables and thinner peat deposits; here the commonly observed permeability decrease with depth is more pronounced than in flatter areas, and groundwater flows decrease significantly between early and late summer as the water table drops. Groundwater flows are also affected by microtopographic features that retain groundwater that could otherwise be released as the active layer deepens. The dominant sources of groundwater, and thus dissolved organic matter, are likely wet, flatter regions with thick organic layers. This finding informs fluid flow and solute transport dynamics for the present and future Arctic.

 
more » « less
Award ID(s):
1637459 1936769 1753731
PAR ID:
10455211
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Volume:
55
Issue:
8
ISSN:
0043-1397
Page Range / eLocation ID:
p. 6555-6576
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    A small imbalance in plant productivity and decomposition accounts for the carbon (C) accumulation capacity of peatlands. As climate changes, the continuity of peatland net C storage relies on rising primary production to offset increasing ecosystem respiration (ER) along with the persistence of older C in waterlogged peat. A lowering in the water table position in peatlands often increases decomposition rates, but concurrent plant community shifts can interactively alter ER and plant productivity responses. The combined effects of water table variation and plant communities on older peat C loss are unknown. We used a full‐factorial 1‐m3mesocosm array with vascular plant functional group manipulations (Unmanipulated Control, Sedge only, and Ericaceous only) and water table depth (natural and lowered) treatments to test the effects of plants and water depth on CO2fluxes, decomposition, and older C loss. We used Δ14C and δ13C of ecosystem CO2respiration, bulk peat, plants, and porewater dissolved inorganic C to construct mixing models partitioning ER among potential sources. We found that the lowered water table treatments were respiring C fixed before the bomb spike (1955) from deep waterlogged peat. Lowered water table Sedge treatments had the oldest dissolved inorganic14C signature and the highest proportional peat contribution to ER. Decomposition assays corroborated sustained high rates of decomposition with lowered water tables down to 40 cm below the peat surface. Heterotrophic respiration exceeded plant respiration at the height of the growing season in lowered water table treatments. Rates of gross primary production were only impacted by vegetation, whereas ER was affected by vegetation and water table depth treatments. The decoupling of respiration and primary production with lowered water tables combined with older C losses suggests that climate and land‐use‐induced changes in peatland hydrology can increase the vulnerability of peatland C stores.

     
    more » « less
  2. Abstract

    Predicting rainfall‐induced landslide motion is challenging because shallow groundwater flow is extremely sensitive to the preexisting moisture content in the ground. Here, we use groundwater hydrology theory and numerical modeling combined with five years of field monitoring to illustrate how unsaturated groundwater flow processes modulate the seasonal pore water pressure rise and therefore the onset of motion for slow‐moving landslides. The onset of landslide motion at Oak Ridge earthflow in California’s Diablo Range occurs after an abrupt water table rise to near the landslide surface 52–129 days after seasonal rainfall commences. Model results and theory suggest that this abrupt rise occurs from the advection of a nearly saturated wetting front, which marks the leading edge of the integrated downward flux of seasonal rainfall, to the water table. Prior to this abrupt rise, we observe little measured pore water pressure response within the landslide due to rainfall. However, once the wetting front reaches the water table, we observe nearly instantaneous pore water pressure transmission within the landslide body that is accompanied by landslide acceleration. We cast the timescale to reach a critical pore water pressure threshold using a simple mass balance model that considers variable moisture storage with depth and explains the onset of seasonal landslide motion with a rainfall intensity‐duration threshold. Our model shows that the seasonal response time of slow‐moving landslides is controlled by the dry season vadose zone depth rather than the total landslide thickness.

     
    more » « less
  3. Abstract

    Emission of CO2from tropical peatlands is an important component of the global carbon budget. Over days to months, these fluxes are largely controlled by water table depth. However, the diurnal cycle is less well understood, in part, because most measurements have been collected daily at midday. We used an automated chamber system to make hourly measurements of peat surface CO2emissions from chambers root‐cut to 30 cm. We then used these data to disentangle the relationship between temperature, water table and heterotrophic respiration (Rhet). We made two central observations. First, we found strong diurnal cycles in CO2flux and near‐surface peat temperature (<10 cm depth), both peaking at midday. The magnitude of diurnal oscillations was strongly influenced by shading and water table depth, highlighting the limitations of relying on daytime measurements and/or a single correction factor to remove daytime bias in flux measurements. Second, we found mean daily Rhethad a strong linear relationship to the depth of the water table, and under flooded conditions, Rhetwas small and constant. We used this relationship between Rhetand water table depth to estimate carbon export from both Rhetand dissolved organic carbon over the course of a year based on water table records. Rhetdominates annual carbon export, demonstrating the potential for peatland drainage to increase regional CO2emissions. Finally, we discuss an apparent incompatibility between hourly and daily average observations of CO2flux, water table and temperature: water table and daily average flux data suggest that CO2is produced across the entire unsaturated peat profile, whereas temperature and hourly flux data appear to suggest that CO2fluxes are controlled by very near surface peat. We explore how temperature‐, moisture‐ and gas transport‐related mechanisms could cause mean CO2emissions to increase linearly with water table depth and also have a large diurnal cycle.

     
    more » « less
  4. Abstract

    Groundwater flow direction within the critical zone of headwater catchments is often assumed to mimic land surface topographic gradients. However, groundwater hydraulic gradients are also influenced by subsurface permeability contrasts, which can result in variability in flow direction and magnitude. In this study, we investigated the relationship between shallow groundwater flow direction, surface topography, and the subsurface topography of low permeability units in a headwater catchment at the Hubbard Brook Experimental Forest (HBEF), NH. We continuously monitored shallow groundwater levels in the solum throughout several seasons in a well network (20 wells of 0.18–1.1 m depth) within the upper hillslopes of Watershed 3 of the HBEF. Water levels were also monitored in four deeper wells, screened from 2.4 to 6.9 m depth within glacial drift of the C horizon. We conducted slug tests across the well network to determine the saturated hydraulic conductivity (Ksat) of the materials surrounding each well. Results showed that under higher water table regimes, groundwater flow direction mimics surface topography, but under lower water table regimes, flow direction can deviate as much as 56 degrees from surface topography. Under these lower water table conditions, groundwater flow direction instead followed the topography of the top of the C horizon. The interquartile range ofKsatwithin the C horizon was two orders of magnitude lower than within the solum. Overall, our results suggest that the land surface topography and the top of the C horizon acted as end members defining the upper and lower bounds of flow direction variability. This suggests that temporal dynamics of groundwater flow direction should be considered when calculating hydrologic fluxes in critical zone and runoff generation studies of headwater catchments that are underlain by glacial drift.

     
    more » « less
  5. Abstract

    Mountain block systems are critical to water resources and have been heavily studied and modeled in recent decades. However, due to lack of field data, there is little consistency in how models represent the mountain block subsurface. While there is a large body of research on subsurface heterogeneity, few studies have evaluated the effect that common conceptual choices modelers make in mountainous systems have on simulated hydrology. Here we simulate the hydrology of a semi‐idealized headwater catchment using six common conceptual models of the mountain block subsurface. These scenarios include multiple representations of hydraulic conductivity decaying with depth, changes in soil depth with topography, and anisotropy. We evaluate flow paths, discharge, and water tables to quantify the impact of subsurface conceptualization on hydrologic behavior in three dimensions. Our results show that adding higher conductivity layers in the shallow subsurface concentrates flow paths near the surface and increases average saturated flow path velocities. Increasing heterogeneity by adding additional layers or introducing anisotropy increases the variance in the relationship between the age and length of saturated flow paths. Discharge behavior is most sensitive to heterogeneity in the shallow subsurface layers. Water tables are less sensitive to layering than they are to the overall conductivity in the domain. Anisotropy restricts flow path depths and controls discharge from storage but has little effect on governing runoff. Differences in the response of discharge, water table depth, and residence time distribution to subsurface representation highlight the need to consider model applications when determining the level of complexity that is needed.

     
    more » « less