Geology topics

Michelle A. Walvoord

Publications and source records attributed to Michelle A. Walvoord.

At least 19 recordsLinked to original sources

USGS2026

Inverse modeling for post-wildfire hydrologic process understanding over snow-dominated watersheds

Wildfire-induced changes to vegetation and soil properties can alter snowpack dynamics and hydrologic processes in ways that complicate predictions of post-wildfire water availability. To better understand these effects, we coupled the Parameter ESTimation in C++ (PEST++) software with the Weather Research & Forecasting Hydrologic model (WRF-Hydro) across three burned, snow-dominated watersheds in the western United States. We incorporated three burn-severity-based land cover categories into WRF-Hydro and used an iterative Ensemble Smoother (iES) within PEST++ to estimate model parameters controlling snow accumulation and melt, infiltration, and streamflow generation. This framework supported three complementary objectives: (1) assessing parameter sensitivities under pre- and post-wildfire conditions, (2) comparing pre- and post-wildfire parameter ensembles to identify wildfire-driven shifts in hydrologic processes, and (3) using counterfactual simulations to isolate the hydrologic response to wildfire from interannual climate variability. Our findings indicate that model parameter sensitivity varied by domain and shifted following wildfire. Maximum snow albedo and snow albedo decay were most influential in the higher-elevation East Troublesome domain, while the rain-snow partitioning temperature threshold dominated in the lower-elevation Chippy Creek and Las Conchas domains. Following wildfire, streamflow sensitivity shifted away from snowpack parameters toward soil hydraulic parameters, consistent with wildfire-driven reductions in evapotranspiration and increases in soil moisture. Incorporating wildfire effects led to complex changes in simulated snowpack accumulation and streamflow timing, generally leading to earlier snowmelt freshet initiation and cessation and reduced late season flows. Simulations were more sensitive to interannual climate variations than wildfire land cover changes. These results highlight both the potential and challenges of ensemble-based inverse modeling for understanding post-wildfire hydrologic change in snow-dominated watersheds.

Colorado, Montana, New Mexico↗

USGS2026

As above, so below? A framework for integrating long-term water quantity trends reveals divergent patterns in groundwater and low streamflow across the United States

Climate, land-use, and disturbance drive long-term global trends in groundwater levels and streamflow. At large scales, these trends are typically considered separately, despite the well-established concept that groundwater and surface water comprise a single resource. Joint trend assessment at national scales is challenging because it requires pairing and aggregating data from spatially disparate streamflow and groundwater monitoring sites for which no established framework exists. Here, we evaluate alternative approaches for integrating groundwater and streamflow data to enable joint trend analysis—a critical step toward understanding how water-budget components respond concurrently and interactively to environmental drivers. Mann–Kendall trends were computed for individual groundwater (annual mean depth) and streamflow (annual low of 7 d averages) sites across the U.S over 21- (2000–2020), 31- (1990–2020), and 41-year (1980–2020) periods. Regional Kendall trends were calculated using five spatially contiguous and noncontiguous regional classifications for aggregation based on subsurface (e.g. aquifer, geology) and surface (e.g. watershed, landscape) characteristics. Site-level results revealed contrasting trends, with tendencies toward increasing low flows (wetting) and increasing groundwater depths (drying). Agreement between streamflow and groundwater trends increased with regional aggregation and longer timeframes, though persistent skew toward streamflow wetting and groundwater drying remained. Results varied by region and trend period, with notable consistencies: unified drying in the West/Southwest and wetting in the Upper Midwest. Directional mismatches in long-term trends were prominent in the High Plains and Mississippi Alluvial Plain, whereas near-term mismatches were most evident in the Northwest. Aggregation by hydrologic landscape regions (HLR) yielded the greatest agreement between groundwater and streamflow trends. These findings indicate that coupled responses may represent combined influences of climate, relief, and geology, as captured by HLR, more strongly than geography or geology alone. Integrated water availability assessments may benefit from a multi-characteristic classification framework to treat groundwater and surface water as a unified resource.

Environmental Research: Water↗

USGS2026

Ice sheet dynamics drive pronounced changes in the subsurface freshwater-saltwater interface

Saltwater is migrating into freshwater aquifers globally with water quality and biogeochemical implications, yet saltwater intrusion in glaciated regions is sparsely investigated. Field observations suggest that groundwater head in glaciated systems is influenced by ice sheet forcings and provides evidence that seawater infiltrated into offshore aquifers during past deglaciation events. To understand links between ice sheet dynamics, groundwater head, and saltwater intrusion, we use numerical models to explore the effects of deglaciation on nearshore head and salinity distributions. We find that ice sheet thinning diminishes groundwater head, and the resulting shift in subsurface pressure gradients drives rapid landward movement of the subsurface freshwater-saltwater interface up to 4.0 km or 1.3 m per m ice sheet loss. Results highlight an overlooked saltwater intrusion mechanism that aligns with field observations and affects glaciated coastlines undergoing ice sheet retreat, underscoring the need to consider this mechanism in studies of contemporary coastal water quality.

Geophysical Research Letters↗

USGS2026

A review and synthesis of post-wildfire shifts in hydrologic processes and streamflow generation mechanisms

Critical water supply watersheds in the western United States (WUS) are impacted by wildfires, with potential negative effects on water quality and quantity. Scientific understanding is currently insufficient to deliver estimates of wildfire consequences for water quantity that are regionally accurate. Regional variability in the directionality and magnitude of post-wildfire shifts in streamflow generation fuels uncertainty in estimates of wildfire effects on water supply. In this work we provide a narrative review of wildfire effects on hydrologic processes and the resulting changes in streamflow generation mechanisms with a focus on the WUS, incorporating other global regions when pertinent. A conceptual model summary of wildfire effects on streamflow generation emphasizes: (1) precipitation seasonality, (2) synchrony of precipitation and potential evapotranspiration, (3) net shifts in interception, evaporation, and transpiration relative to total annual precipitation, (4) vegetation changes, including compensatory uptake and type conversion, (5) degree of overlap in rainfall rates and infiltration, (6) fire extent and severity, (7) burn scar positioning (e.g. in headwaters or proximal to watershed outlet), (8) scale-dependent groundwater leakage, (9) near-surface water storage reduction, and (10) soil to groundwater connectivity. Ongoing gaps and challenges include separating the influences of precipitation variability, water withdrawals, and post-fire land management; compound and overlapping disturbances; and lack of pre-fire data. Notable future opportunities include: harnessing ever-improving gridded and remotely sensed precipitation and fire-effects data; linking geophysical, isotopic tracer, and geochemical signatures to diagnose hydrologic changes; leveraging physically based and data-driven model advancements; and analyzing streamflow generation recovery trajectories across diverse watersheds.

western United States↗

USGS2025

Permafrost–wildfire interactions: active layer thickness estimates for paired burned and unburned sites in northern high latitudes

As the northern high-latitude permafrost zone experiences accelerated warming, permafrost has become vulnerable to widespread thaw. Simultaneously, wildfire activity across northern boreal forest and Arctic/subarctic tundra regions impacts permafrost stability through the combustion of insulating organic matter, vegetation, and post-fire changes in albedo. Efforts to synthesis the impacts of wildfire on permafrost are limited and are typically reliant on antecedent pre-fire conditions. To address this, we created the FireALT dataset by soliciting data contributions that included thaw depth measurements, site conditions, and fire event details with paired measurements at environmentally comparable burned and unburned sites. The solicitation resulted in 52 466 thaw depth measurements from 18 contributors across North America and Russia. Because thaw depths were taken at various times throughout the thawing season, we also estimated end-of-season active layer thickness (ALT) for each measurement using a modified version of the Stefan equation. Here, we describe our methods for collecting and quality-checking the data, estimating ALT, the data structure, strengths and limitations, and future research opportunities. The final dataset includes 48 669 ALT estimates with 32 attributes across 9446 plots and 157 burned–unburned pairs spanning Canada, Russia, and the United States. The data span fire events from 1900 to 2022 with measurements collected from 2001 to 2023. The time since fire ranges from 0 to 114 years. The FireALT dataset addresses a key challenge: the ability to assess impacts of wildfire on ALT when measurements are taken at various times throughout the thaw season depending on the time of field campaigns (typically June through August) by estimating ALT at the end-of-season maximum. This dataset can be used to address understudied research areas, particularly algorithm development, calibration, and validation for evolving process-based models as well as extrapolating across space and time, which could elucidate permafrost–wildfire interactions under accelerated warming across the high-northern-latitude permafrost zone. The FireALT dataset is available through the Arctic Data Center ( https://doi.org/10.18739/A2RN3092P , Talucci et al., 2024).

Arctic↗

USGS2025

Disparate groundwater responses to wildfire

Post-wildfire investigations of groundwater response reveal a range of outcomes, varying from substantial increases to notable decreases in recharge and baseflow, with some studies indicating negligible or short-lived effects. This review assesses these varied responses within five critical categories: climate, vegetation, hydrogeology, fire characteristics, and the cryosphere, examining both short-term (within 2 years) and intermediate (2–10 years post-fire) effects. Despite considerable variability, some consistent patterns emerge. For instance, in hydroclimatic settings where water input and evaporative demand cycles are out of sync, post-wildfire groundwater responses tend to be positive (i.e., increased flux or storage), whereas under low fire severity conditions or in vegetation types that quickly recover, groundwater responses tend to be negative (i.e., decreased flux or storage). We synthesize relevant findings into a compendium of testable hypotheses aimed at explaining the spatiotemporal variability in observed post-wildfire groundwater responses. A recurring theme is the critical influence of the pre-wildfire groundwater regime on expected response and recovery. We identify opportunities for specific improvements in post-wildfire monitoring and modeling that would further advance capabilities to predict groundwater response. A key area for further research is understanding how wildfire effects on snow dynamics and other cryospheric processes translate to changes in groundwater.

WIREs Water↗

USGS2025

A review of post-wildfire adaptations of surface-water-quality models: Synthesis, gaps, and opportunities

As wildfires increasingly affect water-supply watersheds, the demand for models to predict water-quality responses is increasing. This work reviews and synthesizes existing post-wildfire applications of water-quality models in the context of geographic and ecohydrological distribution, hydrologic and water-quality response process representation, model parameterization, model and input data scales, model calibration data availability, as well as calibration and performance evaluation approaches. Emphasis is placed on models that simulate water-quality output, rather than sediment and erosional response as the primary focus. Here, identified gaps and opportunities to advance the post-wildfire application of water-quality models include: 1. applying models in under-represented geographic and ecohydrologic regions, 2. simulating multiple streamflow generation mechanisms, including groundwater, with an emphasis on shifting dominant flow pathways as the landscape recovers following wildfire, 3. adding studies that include the simulation of metals, 4. incorporating more biogeochemical and in-stream processes to model applications, 5. applying finer spatial and temporal resolution of precipitation data input as well as finer spatial resolution hydrologic response units, 6. implementing fully distributed grid or element models or finer resolution response units to capture burn severity heterogeneity, 7. collecting enhanced water-quality data for model calibration and validation, 8. conducting model-intercomparison studies, and 9. developing model parameter value guidance in post-wildfire applications. These identified gaps and opportunities may assist users in deciding on key processes and approaches to consider in modeling post-wildfire water-quality conditions.

Science of the Total Environment↗

USGS2024

Opportunities and challenges for precipitation forcing data in post-wildfire hydrologic modeling applications

The frequency and extent of wildfires have increased in recent decades with immediate and cascading effects on water availability in many regions of the world. Precipitation is used as primary input to hydrologic models and is a critical driver of post-wildfire hydrologic hazards including debris flows, flash floods, water-quality effects, and reservoir sedimentation. These models are valuable tools for understanding the hydrologic response to wildfire but require accurate precipitation data at suitable spatial and temporal resolutions. Wildfires often occur in data-sparse, headwater catchments in complex terrain, and post-wildfire hydrologic effects are particularly sensitive to high-intensity, short-duration precipitation events, which are highly variable and difficult to measure or estimate. Therefore, the assessment and prediction of wildfire-induced changes to watershed hydrology, including the associated effects on ecosystems and communities, are complicated by uncertainty in precipitation data. When direct measurements of precipitation are not available, datasets of indirect measurements or estimates are often used. Choosing the most appropriate precipitation dataset can be difficult as different datasets have unique trade-offs in terms of spatial and temporal accuracy, resolution, and completeness. Here, we outline the challenges and opportunities associated with different precipitation datasets as they apply to post-wildfire hydrologic models and modeling objectives. We highlight the need for expanded precipitation gage deployment in wildfire-prone areas and discuss potential opportunities for future research and the integration of precipitation data from disparate sources into a common hydrologic modeling framework.

WIREs Water↗

USGS2024

Thermo-hydrologic processes governing supra-permafrost talik dynamics in discontinuous permafrost near Umiujaq (Québec, Canada)

Widespread supra-permafrost talik formation is currently recognized as a critical mechanism that could accelerate permafrost thaw in the Arctic (e.g., Connon et al. 2018; Farquharson et al. 2022). However, the trajectory of permafrost dynamics following talik formation may prove difficult to predict. Physically-based cryohydrogeologic models provide a powerful tool for understanding processes and factors controlling talik dynamics and, ultimately, how permafrost will respond to climate change. Such models are typically used to represent multiple non-linear processes relevant for groundwater systems in cold regions, such as coupled heat and groundwater movement, including freeze-thaw dynamics and the effects on the surface energy balance and the subsurface thermal and hydraulic properties (Lamontagne-Hallé et al. 2020). Though cryohydrogeologic modeling advances have been made in simulating talik dynamics, few applications have been tested against robust long-term hydrometeorological and subsurface observations.

Quebec↗

USGS2024

A history of cryohydrogeology modeling and recent advancements through the integration of solute transport

Groundwater flow systems and permafrost are interrelated because permafrost thaw enhances permeability, while groundwater flow can advect heat and accelerate permafrost thaw (McKenzie et al. 2021). Given amplified climate change in cold regions, there is renewed interest in ‘cryohydrogeology’, the study of groundwater in cold regions. Many data-driven studies have shown that permafrost thaw is leading to activated aquifers and increased baseflow across the pan-Arctic region (e.g. Walvoord and Striegl 2007, Evans et al. 2020). Empirical evidence of a subsurface ‘replumbing’ (Walvoord and Kurylyk 2016) in permafrost regions raises questions about the fate of sequestered contaminants in the North (Langer et al. 2023). We will discuss the history of and emerging opportunities in cryohydrogeological modeling, with a focus on recent contaminant transport modeling.

Conference Paper↗

USGS2023

Wind-modulated groundwater discharge along a microtidal Arctic coastline

Groundwater discharge transports dissolved constituents to the ocean, affecting coastal carbon budgets and water quality. However, the magnitude and mechanisms of groundwater exchange along rapidly transitioning Arctic coastlines are largely unknown due to limited observations. Here, using first-of-its-kind coastal Arctic groundwater timeseries data, we evaluate the magnitude and drivers of groundwater discharge to Alaska's Beaufort Sea coast. Darcy flux calculations reveal temporally variable groundwater fluxes, ranging from −6.5 cm d −1 (recharge) to 14.1 cm d −1 (discharge), with fluctuations in groundwater discharge or aquifer recharge over diurnal and multiday timescales during the open-water season. The average flux during the monitoring period of 4.9 cm d −1 is in line with previous estimates, but the maximum discharge exceeds previous estimates by over an order-of-magnitude. While the diurnal fluctuations are small due to the microtidal conditions, multiday variability is large and drives sustained periods of aquifer recharge and groundwater discharge. Results show that wind-driven lagoon water level changes are the dominant mechanism of fluctuations in land–sea hydraulic head gradients and, in turn, groundwater discharge. Given the microtidal conditions, low topographic relief, and limited rainfall along the Beaufort Sea coast, we identify wind as an important forcing mechanism of coastal groundwater discharge and aquifer recharge with implications for nearshore biogeochemistry. This study provides insights into groundwater flux dynamics along this coastline over time and highlights an oft overlooked discharge and circulation mechanism with implications towards refining solute export estimates to coastal Arctic waters.

Alaska↗

USGS2023

Integrated water resources trend assessments: State of the science, challenges, and opportunities for advancement

Water is vital to human life and healthy ecosystems. Here we outline the current state of national-scale water resources trend assessments, identify key gaps, and suggest advancements to better address critical issues related to changes in water resources that may threaten human development or the environment. Questions like, “Do we have less suitable drinking water now than we had 20 years ago?” or “Are flood events more common now than they were in the past?” prompted improvements in data, trend estimation methods, and modeling frameworks to track changes in, and better understand how land use and climate influence four water resources domains: surface and groundwater quantity and quality. However, continued advancement in trend assessments to better address issues related to changes in water availability is needed. Areas of need include more timely and efficient delivery of water resources trend results and improved capacity to estimate trends at unmonitored locations. Additional integration pieces include increased understanding of groundwater–surface water interactions, incorporation of both quantity and quality trends into water availability estimates, and the refinement of trend metrics to account for the competing needs of society and ecological integrity. Coupled with improved driver attribution studies, these components will better inform current and future water resources management.

Journal of the American Water Resource Association↗

USGS2023

Modeling post-wildfire hydrologic response: Review and future directions for applications of physically based distributed simulation

Wildfire is a growing concern as climate shifts. The hydrologic effects of wildfire, which include elevated hazards and changes in water quantity and quality, are increasingly assessed using numerical models. Post-wildfire application of physically based distributed models provides unique insight into the underlying processes that affect water resources after wildfire. This work reviews and synthesizes post-wildfire applications of physically based distributed models by examining the scales and geographic/ecohydrologic distribution of model applications, hydrologic response process representation, model parameterization, and model performance metrics. Highlighted gaps and opportunities for advancing physically based distributed hydrologic response modeling after wildfire include the following: (a) applying models in under-represented geographic (S. America, Africa, Asia) and ecohydrologic regions (arid or dry subhumid climates), (b) incorporating all four major streamflow generation mechanisms (infiltration excess, saturation excess, subsurface storm flow, and groundwater flow), (c) representing integrated vadose zone and saturated zone processes to better capture subsurface streamflow generation, (d) building new remotely sensed model parameterization methods for precipitation interception, infiltration, and overland flow that account for burn severity and recovery, (e) incorporating distributed state variables (e.g., soil moisture, groundwater levels) in model performance assessment, (f) designing model intercomparison studies, including field datasets specifically for post-wildfire model development and validation, (g) linking mechanistic vegetation regrowth models with hydrologic models to improve simulation of process shifts as ecosystems recover, and (h) creating a new community modeling framework to integrate modeling advances across the wildfire science community.

Earth's Future↗

USGS2023

Long-term, high-resolution permafrost monitoring reveals coupled energy balance and hydrogeologic controls on talik dynamics near Umiujaq (Nunavik, Québec, Canada)

Rising temperatures in the Arctic and subarctic are driving the rapid thaw of permafrost by reducing permafrost cooling, increasing active layer thickness, and promoting talik formation. In this study, the cyrohydrogeology of a permafrost mound located within the discontinuous permafrost zone near Umiujaq (Nunavik, Québec, Canada) is characterized through the analysis of a dataset covering more than two decades of monitoring. This dataset captures a high degree of interannual variability in air temperature and ground thermal conditions, as well as the formation and closure of a supra-permafrost talik. Data indicate that variable saturation and advective heat transport directly contribute to the expansion and contraction of the talik. Data further indicate the presence of two distinct thermo-hydrologic settings resulting from differences in surface conditions, as well as subsurface thermal and flow regimes. The first, found at the top of the mound feature, is characterized by very low moisture contents (< 0.05 m 3 /m 3 ), while the second, found at the side of the mound feature, shows higher annual moisture contents that strongly influence the dynamics of heat and groundwater flow. The data were synthesized into a detailed conceptual model of the cyrohydrogeological dynamics that highlights the important role of hydrogeological characterization and long-term datasets in understanding the effects of groundwater flow on seasonal frost and permafrost dynamics. Specifically, the results presented here show that in the absence of long-term datasets, longer-period transient phenomena such as talik opening and closure may be misrepresented as uni-directional feedback loops, as opposed to highly-dynamic temporary phenomena.

Nunavik, Quebec↗

USGS2023

Wildfire-induced shifts in groundwater discharge to streams identified with paired air and stream water temperature analyses

Within the western United States, increasingly severe and frequent wildfires may alter the magnitude, timing, and quality of water exported from burned areas by streams. Post-fire hydrologic studies often focus on peak stream flow responses to shifts in runoff generation or on annual streamflow yield response to changes in evapotranspiration following fire. However, the magnitude and duration of wildfire effects on groundwater recharge, changes in subsurface routing, and consequences for stream low flows sourced predominately by baseflow are poorly understood. Here, we demonstrate an approach using the amplitude and phase of paired annual air and stream water temperature signals to broadly identify changes in watershed subsurface flow contributions after fire. Watersheds were classified using pre-fire temperature data, as having air-coupled (i.e., reduced apparent groundwater signature), deep groundwater, or shallow groundwater stream temperature signals. Changes in pre- and post-fire paired air and stream water temperature metrics were compared for locations (n = 17) spanning a large range of physiographic and climatic conditions across the western United States. Pre- and post-fire comparisons were computed by quantile using bootstrapped confidence intervals (ci = 95), as well as in aggregate using Kruskal-Wallis and post-hoc Dunn tests. Statistical comparisons of pre- and post-fire temperature metrics suggest that overall, watersheds classified as having minimal groundwater influence are the most likely to experience fire-induced subsurface hydrologic change. More specifically, watersheds classified as having air-coupled or shallow groundwater signals experienced increases in the magnitude of groundwater discharge, with more stable annual thermal regimes post-fire that are less-coupled to ambient air temperature. These findings form the basis of a conceptual framework for watershed resistance to subsurface hydrologic change following fire that can be broadly applied as a first approximation for water management, impacts on aquatic habitat, and post-wildfire response planning.

California, Idaho, Oregon, Washington↗

USGS2023

Sentinel responses of Arctic freshwater systems to climate: linkages, evidence, and a roadmap for future research

While the sentinel nature of freshwater systems is now well recognized, widespread integration of freshwater processes and patterns into our understanding of broader climate-driven Arctic terrestrial ecosystem change has been slow. We review the current understanding across Arctic freshwater systems of key sentinel responses to climate, which are attributes of these systems with demonstrated and sensitive responses to climate forcing. These include ice regimes, temperature and thermal structure, river baseflow, lake area and water level, permafrost-derived dissolved ions and nutrients, carbon mobilization (dissolved organic carbon, greenhouse gases, and radiocarbon), dissolved oxygen concentrations, lake trophic state, various aquatic organisms and their traits, and invasive species. For each sentinel, our objectives are to clarify linkages to climate, describe key insights already gained, and provide suggestions for future research based on current knowledge gaps. We suggest that tracking key responses in Arctic freshwater systems will expand understanding of the breadth and depth of climate-driven Arctic ecosystem changes, provide early indicators of looming, broader changes across the landscape, and improve protection of freshwater biodiversity and resources.

Arctic Science↗

USGS2022

Sea-level rise and warming mediate coastal groundwater discharge in the Arctic

Groundwater discharge is an important mechanism through which fresh water and associated solutes are delivered to the ocean. Permafrost environments have traditionally been considered hydrogeologically inactive, yet with accelerated climate change and permafrost thaw, groundwater flow paths are activating and opening subsurface connections to the coastal zone. While warming has the potential to increase land-sea connectivity, sea-level change has the potential to alter land-sea hydraulic gradients and enhance coastal permafrost thaw, resulting in a complex interplay that will govern future groundwater discharge dynamics along Arctic coastlines. Here, we use a recently developed permafrost hydrological model that simulates variable-density groundwater flow and salinity-dependent freeze-thaw to investigate the impacts of sea-level change and land and ocean warming on the magnitude, spatial distribution, and salinity of coastal groundwater discharge. Results project both an increase and decrease in discharge with climate change depending on the rate of warming and sea-level change. Under high warming and low sea-level rise scenarios, results show up to a 58% increase in coastal groundwater discharge by 2100 due to the formation of a supra-permafrost aquifer that enhances freshwater delivery to the coastal zone. With higher rates of sea-level rise, the increase in discharge due to warming is reduced to 21% as sea-level rise decreased land-sea hydraulic gradients. Under lower warming scenarios for which supra-permafrost groundwater flow was not established, discharge decreased by up to 26% between 1980 and 2100 for high sea-level rise scenarios and increased only 8% under low sea-level rise scenarios. Thus, regions with higher warming rates and lower rates of sea-level change (e.g. northern Nunavut, Canada) will experience a greater increase in discharge than regions with lower warming rates and higher rates of sea-level change. The magnitude, location and salinity of discharge have important implications for ecosystem function, water quality, and carbon dynamics in coastal zones.

Environmental Research Letters↗

USGS2022

Heterogeneous patterns of aged organic carbon export driven by hydrologic flow paths, soil texture, fire, and thaw in discontinuous permafrost headwaters

Climate change is thawing and potentially mobilizing vast quantities of organic carbon (OC) previously stored for millennia in permafrost soils of northern circumpolar landscapes. Climate-driven increases in fire and thermokarst may play a key role in OC mobilization by thawing permafrost and promoting transport of OC. Yet, the extent of OC mobilization and mechanisms controlling terrestrial-aquatic transfer are unclear. We demonstrate that hydrologic transport of soil dissolved OC (DOC) from the active layer and thawing permafrost to headwater streams is extremely heterogeneous and regulated by the interactions of soils, seasonal thaw, fire, and thermokarst. Repeated sampling of streams in eight headwater catchments of interior Alaska showed that the mean age of DOC for each stream ranges widely from modern to ∼2,000 years B.P. Together, an endmember mixing model and nonlinear, generalized additive models demonstrated that Δ 14 C-DOC signature (and mean age) increased from spring to fall, and was proportional to hydrologic contributions from a solute-rich water source, related to presumed deeper flow paths found predominantly in silty catchments. This relationship was correlated with and mediated by catchment properties. Mean DOC ages were older in catchments with >50% burned area, indicating that fire is also an important explanatory variable. These observations underscore the high heterogeneity in aged C export and difficulty of extrapolating estimates of permafrost-derived DOC export from watersheds to larger scales. Our results provide the foundation for developing a conceptual model of permafrost DOC export necessary for advancing understanding and prediction of land-water C exchange in changing boreal landscapes.

Alaska↗