Search NASA⌕ Search

SEARCH · Search NASA

Results for “snowpack dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

Feedbacks Between Surface Deformation and Permafrost Degradation in Ice Wedge Polygons, Arctic Coastal Plain, Alaska

In the past three decades, an abrupt, pan-Arctic acceleration of ice wedge melting has transformed tundra landscapes, spurring the formation of hummock-like features known as high-centered polygons (HCPs). This rapid geomorphic transition profoundly alters regional hydrology and influences surface emissions of CO 2 and CH 4 . In Arctic Alaska, most recent instances of ice wedge degradation have arrested within 15–20 years of inception, stabilizing HCP microtopography. However, feedbacks between ground surface deformation and permafrost stability are incompletely understood, limiting our capacity to predict trajectories of landscape evolution in a still warmer future. In this work, we use field data from a site near Prudhoe Bay, Alaska, to develop a modeling-based framework for assessing the strength of positive (i.e., exacerbating) feedbacks on ice wedge degradation, focusing on the importance of heterogeneity in surface drainage and microtopographic conditions. Our simulations suggest that, when troughs are narrow, positive feedbacks on ice wedge melting (associated with thermokarst pool formation) are relatively weak. Positive feedbacks are markedly stronger beneath wide troughs, such as those that form above older, larger ice wedges. Seasonal thaw abruptly accelerates once a talik begins to form beneath wide and deep thermokarst pools. Once a talik initiates, winter severity and snowpack thickness increase in importance as predictors of thaw intensity in summer. Our results indicate that meter-scale heterogeneity in polygonal microtopography potentially exerts strong, nonlinear controls on thermokarst trajectories. These findings are useful for predicting future thermokarst dynamics and for interpreting the results from coarser-resolution land surface models operating at greater spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗

Geophysical Monitoring Shows that Spatial Heterogeneity in Thermohydrological Dynamics Reshapes a Transitional Permafrost System

Climate change is causing rapid changes of Arctic ecosystems. Yet, data needed to unravel complex subsurface processes are very rare. Using geophysical and in-situ sensing, this study closes an observational gap associated with thermohydrological dynamics in discontinuous permafrost systems. It highlights the impact of vegetation and snow thickness distribution on subsurface thermohydrological properties and processes. Large snow accumulation near tall shrubs insulates the ground and allows for rapid and downward heat flow. Thinner snowpack above graminoid results in surficial freezing and prevents water from infiltrating into the subsurface. Analyzing short term disturbances, we found that lateral flow could be a driving factor in talik formation. Inter-annual measurements show that deep permafrost temperatures increased by about 0.2°C over two years. The results, which suggest that snow-vegetation-subsurface processes are tightly coupled, will be useful for improving predictions of Arctic feedback to climate change, including how subsurface thermohydrology influences CO 2 and CH 4 fluxes.

58 GEOSCIENCES↗

An overview of the Western United States Dynamically Downscaled Dataset (WUS-D3)

Abstract. Predicting future climate change over a region of complex terrain, such as the western United States (US), remains challenging due to the low resolution of global climate models (GCMs). Yet the climate extremes of recent years in this region, such as floods, wildfires, and drought, are likely to intensify further as climate warms, underscoring the need for high-quality and high-resolution predictions. Here, we present an ensemble of dynamically downscaled simulations over the western US from 1980–2100 at 9 km grid spacing, driven by 16 latest-generation GCMs. This dataset is titled the Western US Dynamically Downscaled Dataset (WUS-D3). We describe the challenges of producing WUS-D3, including GCM selection and technical issues, and we evaluate the simulations' realism by comparing historical results to temperature and precipitation observations. The future downscaled climate change signals are shaped in physically credible ways by the regional model's more realistic coastlines and topography. (1) The mean warming signals are heavily influenced by more realistic snowpack. (2) Mean precipitation changes are often consistent with wetting on the windward side of mountain complexes, as warmer, moister air masses are uplifted orographically during precipitation events. (3) There are large fractional precipitation increases on the lee side of mountain complexes, leading to potentially significant changes in water resources and ecology in these arid landscapes. (4) Increases in precipitation extremes are generally larger than in the GCMs, driven by locally intensified atmospheric updrafts tied to sharper, more realistic gradients in topography. (5) Changes in temperature extremes are different from what is expected by a shift in mean temperature and are shaped by local atmospheric dynamics and land surface feedbacks. Because of its high resolution, comprehensiveness, and representation of relevant physical processes, this dataset presents a unique opportunity to evaluate societally relevant future changes in western US climate.

Rahimi, Stefan (ORCID:0000000331884462)↗

Temperature and Snow-Mediated Moisture Controls of Summer Photosynthetic Activity in Northern Terrestrial Ecosystems Between 1982 and 2011

Recent warming has stimulated the productivity of boreal and Arctic vegetation by reducing temperature limitations. However, several studies have hypothesized that warming may have also increased moisture limitations because of intensified summer drought severity. Establishing the connections between warming and drought stress has been difficult because soil moisture observations are scarce. Here we use recently developed gridded datasets of moisture variability to investigate the links between warming and changes in available soil moisture and summer vegetation photosynthetic activity at northern latitudes (greater than 45N) based on the Normalized Difference Vegetation Index (NDVI) since 1982. Moisture and temperature exert a significant influence on the inter-annual variability of summer NDVI over about 29% (mean r(sup 2) = 0.29 +/ 0.16) and 43% (mean r(sup 2 = 0.25 +/- 0.12) of the northern vegetated land, respectively. Rapid summer warming since the late 1980s (approximately 0.7deg C) has increased evapotranspiration demand and consequently summer drought severity, but contrary to earlier suggestions it has not changed the dominant climate controls of NDVI over time. Furthermore, changes in snow dynamics (accumulation and melting) appear to be more important than increased evaporative demand in controlling changes in summer soil moisture availability and NDVI in moisture-sensitive regions of the boreal forest. In boreal North America, forest NDVI declines are more consistent with reduced snowpack rather than with temperature-induced increases in evaporative demand as suggested in earlier studies. Moreover, summer NDVI variability over about 28% of the northern vegetated land is not significantly associated with moisture or temperature variability, yet most of this land shows increasing NDVI trends. These results suggest that changes in snow accumulation and melt, together with other possibly non-climatic factors are likely to play a significant role in modulating regional ecosystem responses to the projected warming and increase in evapotranspiration demand during the coming decades.

Snowpack↗

The Impact of Ice Sheet Surface Representation on Surface Mass Balance in the Goddard Earth Observing System

Surface conditions and their impacts on surface mass balance (SMB) play important roles in ice sheet dynamics and ice sheet interactions with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) – an ecosystem of models and reanalyses – represent ice sheet and glacier SMB components, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. We find the successful representation of ice sheet surface mass balance in MERRA-2, and similar systems, is partly due to compensating biases in surface energy budget. In the atmospheric system, subsequent changes to cloud microphysics and in the longwave radiative transfer model are found to have reduced biases in surface net longwave radiation fluxes while exacerbating surface net shortwave radiation biases, producing erroneously high near-surface temperatures and surface melt in summer months. These issues can be exacerbated by poor experiment initialization and the use of two-moment cloud microphysics within the GEOS ocean-atmosphere coupled system. Here, we document the spatial and temporal extent and causes of these biases in the ice sheet surface energy budget across GEOS systems and implement a range of model improvements to mitigate these issues. We examine the combined impact on these changes on regional energy budget and SMB using both free running and replay experiments (to simulate the impact in reanalyses). Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet SMB in GMAO’s future reanalyses and forecasting systems.

Lauren C Andrews↗

Performance Tests of Snow-Related Variables Over the Tibetan Plateau and Himalayas Using a New Version of NASA GEOS-5 Land Surface Model that Includes the Snow Darkening Effect

The snow darkening effect, i.e. the reduction of snow albedo, is caused by absorption of solar radiation by absorbing aerosols (dust, black carbon, and organic carbon) deposited on the snow surface. This process is probably important over Himalayan and Tibetan glaciers due to the transport of highly polluted Atmospheric Brown Cloud (ABC) from the Indo-Gangetic Plain (IGP). This effect has been incorporated into the NASA Goddard Earth Observing System model, version 5 (GEOS-5) atmospheric transport model. The Catchment land surface model (LSM) used in GEOS-5 considers 3 snow layers. Code was developed to track the mass concentration of aerosols in the three layers, taking into account such processes as the flushing of the compounds as liquid water percolates through the snowpack. In GEOS-5, aerosol emissions, transports, and depositions are well simulated in the Goddard Chemistry Aerosol Radiation and Transport (GO CART) module; we recently made the connection between GOCART and the GEOS-5 system fitted with the revised LSM. Preliminary simulations were performed with this new system in "replay" mode (i.e., with atmospheric dynamics guided by reanalysis) at 2x2.5 degree horizontal resolution, covering the period 1 November 2005 - 31 December 2009; we consider the final three years of simulation here. The three simulations used the following variants of the LSM: (1) the original Catchment LSM with a fixed fresh snowfall density of 150 kg m-3 ; (2) the LSM fitted with the new snow albedo code, used here without aerosol deposition but with changes in density formulation and melting water effect on snow specific surface area, (3) the LSM fitted with the new snow albedo code as same as (2) but with fixed aerosol deposition rates (computed from GOCART values averaged over the Tibetan Plateau domain [Ion.: 60-120E; lat.: 20-50N] during March-May 2008) applied to all grid points at every time step. For (2) and (3), the same setting on the fresh snowfall density as in (1) was used.

Yasunari, Tppei J.↗

Snow Reflectance from LANDSAT-4 Thematic Mapper

In California 75% of the agricultural water supply comes from the melting Sierra Nevada snowpack. Basin-wide albedo measurements from the LANDSAT-4 Thematic Mapper could be used to better forecast the timing of the spring runoff, because these data can be combined with solar radiation calculations to estimate the net radiation budget. The TM is better-suited for this purpose than the MSS because of its large dynamic range. Saturation still occurs in bands 1-4, but is severe only in band 1. Differentiation of snow optical grain size is possible with TM band 4 through a moderately clear atmosphere. The TM band 5 can discriminate clouds from snow, and the combination of bands 2 and 5 appears best for snow mapping.

Dozier, J.↗

Regional Impacts of Ice Sheet Surface Representation Improvements in an Earth System Model

Ice sheet surface conditions and their subsequent impact on surface mass balance play an important role in ice sheet dynamics and ice sheet interaction with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) ecosystem of models and reanalyses model ice sheet and glacier surface mass balance with a moderately complex snow and ice module, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. With this configuration, Greenland near surface air temperatures and ice sheet melt are well represented in MERRA-2, GMAO’s current atmospheric reanalysis; however, the similarity to observations in some regions is due, in part, to compensating biases in surface energy budget. Recent improvements to GEOS’ cloud microphysics and longwave radiation parameterizations reduced biases in surface net longwave radiation and exacerbated surface net shortwave radiation biases, effectively increasing near surface temperatures and melting, particularly during the summer months. Here, we demonstrate the Earth system response to an improved representation of the ice sheet surface focused on minimizing the surface energy biases and increasing spatial and temporal representativeness. New changes to the surface albedo parameterization, including the introduction of a prognostic aerosol darkening scheme, help mitigate shortwave biases, while modifications to surface roughness characteristics improve near surface sensible heat fluxes. Focusing on the Arctic and Greenland, we examine the combined impact on these changes on regional weather conditions and ice sheet surface mass balance using a high resolution atmospheric only experiment, and we use a coupled model configuration to assess the impact on oceanic conditions and subseasonal-to-seasonal prediction. Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet mass balance in GMAO’s future reanalyses and forecasting systems.

Lauren C Andrews↗

Regional Impacts of Ice Sheet Surface Representation Improvements in an Earth System Model

Ice sheet surface conditions and their subsequent impact on surface mass balance play an important role in ice sheet dynamics and ice sheet interaction with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) ecosystem of models and reanalyses model ice sheet and glacier surface mass balance with a moderately complex snow and ice module, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. With this configuration, Greenland near surface air temperatures and ice sheet melt are well represented in MERRA-2, GMAO’s current atmospheric reanalysis; however, the similarity to observations in some regions is due, in part, to compensating biases in surface energy budget. Recent improvements to GEOS’ cloud microphysics and longwave radiation parameterizations reduced biases in surface net longwave radiation and exacerbated surface net shortwave radiation biases, effectively increasing near surface temperatures and melting, particularly during the summer months. Here, we demonstrate the Earth system response to an improved representation of the ice sheet surface focused on minimizing the surface energy biases and increasing spatial and temporal representativeness. New changes to the surface albedo parameterization, including the introduction of a prognostic aerosol darkening scheme, help mitigate shortwave biases, while modifications to surface roughness characteristics improve near surface sensible heat fluxes. Focusing on the Arctic and Greenland, we examine the combined impact on these changes on regional weather conditions and ice sheet surface mass balance using a high resolution atmospheric only experiment, and we use a coupled model configuration to assess the impact on oceanic conditions and subseasonal-to-seasonal prediction. Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet mass balance in GMAO’s future reanalyses and forecasting systems.

Lauren C. Andrews↗

Diverging Trends in Rain-on-Snow Over High Mountain Asia

Rain-on-snow (ROS) over snow-dominated regions such as High Mountain Asia (HMA) modulates snowmelt and runoff and is key contributor in influencing water availability and hazards (e.g., floods and landslides). We studied the trends in ROS in HMA over the past two decades from 2001 to 2018 using the land surface model Noah-MP driven by an ensemble precipitation data set. Our results show that changes in precipitation phase and rainfall are altering ROS. Because of the strong physical heterogeneity and atmospheric dynamics of HMA, ROS characteristics and trends are region-dependent and ROS occurs predominantly over the Indus, Ganges-Brahmaputra, and northwestern basins. In the Indus, ROS representing ∼5% of the annual precipitation and ∼20% of the annual snowmelt, has an increasing trend. This is contrary to the Ganges-Brahmaputra characterized by decreasing ROS trends, where it represents ∼11% of the annual precipitation and ∼60% of the annual snowmelt. In the northwestern basins, ROS has bidirectional trends due to elevation patterns and trends in rainfall, and it constitutes ∼5 to ∼10% of the annual precipitation. Increasing trends in ROS over Indus contribute to reducing the snowpack in late summer, with concerns of reduced water availability and increased groundwater exploitation. Similarly, because of its high amount and contribution to snowmelt, the decreasing ROS trends in the Ganges-Brahmaputra will have consequences of decreased recharge from the headwaters and exacerbated use of groundwater unless increasing trends in rainfall compensate for the decreasing snowmelt. These results provide new insights on ROS-driven changes in the hydrological cycle over HMA.

rain-on-snow↗

Snow reflectance from Landsat-4 Thematic Mapper

In California 75 percent of the agricultural water supply comes from the melting Sierra Nevada snowpack. Basin-wide spectral albedo measurements from the Landsat-4 Thematic Mapper (TM) could be used to better forecast the timing of the spring runoff, because these data can be combined with solar radiation calculations to estimate the net radiation balance. The TM is better-suited for this purpose than the Multispectral Scanner because of its larger dynamic range. Saturation still occurs in bands 1-4, but is severe only in TM1 (0.45-0.52 micron). Snow reflectance in TM2 (0.43-0.61 micron) is typical of the visible wavelength region, where reflectance is almost insensitive to crystal size but sensitive to contamination. TM4 (0.78-0.90 micron) allows estimation of effective optical grain size and thereby spectral extension throughout the near-infrared. TM5 (1.57-1.78 microns) can discriminate clouds from snow.

Dozier, J.↗

Model and remote-sensing-guided experimental design and hypothesis generation for monitoring snow-soil–plant interactions

In this study, we develop a machine-learning (ML)-enabled strategy for selecting hillslope-scale ecohydrological monitoring sites within snow-dominated mountainous watersheds, with a particular focus on snow-soil–plant interactions. Data layers rely on spatial data layers from both remote sensing and hydrological model simulations. Specifically, a Landsat-based foresummer drought sensitivity index is used to define the dependency of the annual peak plant productivity on the Palmer drought severity index in the early growing season. Hydrological simulations provide the spatiotemporal dynamics of near-surface soil moisture and snow depth. In this framework, a regression analysis identifies the key hydrological variables relevant to the spatial heterogeneity of drought sensitivity. We then apply unsupervised clustering to these key variables, using the Gaussian mixture model, to group hillslopes into several zones that have divergent relationships regarding soil moisture, snow dynamics, and drought sensitivity. Using the datasets collected in the East River Watershed (Crested Butte, Colorado, United States), results show that drought sensitivity is significantly correlated with model-derived soil moisture and snow-free timing over space and time. The relationship is, however, non-linear, such that the correlation decreases above a threshold elevation and in a heavy snow year due to large snowpacks, lateral flow, and soil storage limitations. Clustering is then able to define the zones that have high or low sensitivity to drought, as well as the mid-elevation regions where sensitivity is associated with the topographic aspect and net potential radiation. In addition, the algorithm identifies the most representative hillslopes with road/trail access within each zone for installing monitoring sites. Our method also aims to significantly increase the use of ML and model-simulation results to guide critical zone and watershed monitoring activities.

54 ENVIRONMENTAL SCIENCES↗

Unprecedented Shifts in Hydrology Are Emerging Across California's Critical Basins: An Evaluation From 0.5 to 3.5°C

With advances in climate models and downscaling techniques, stakeholders anticipate high-resolution analysis to inform regional to local changes in water management. Here, we produce hydrologic projections from an ensemble of Earth System Models (ESMs) that were selected and downscaled to support California's 5th Climate Assessment. An ensemble of 19 ESMs was downscaled to a 3-km resolution across California using a statistical-dynamical downscaling approach and subsequently run through two calibrated hydrology models. Although California has been extensively studied in the context of climate change, we provide the first evaluation of the warming thresholds at which hydroclimate metrics demonstrate statistically significant shifts. We show that present-day to near-term warming levels in Klamath and Northern Sierra Nevada basins, which serve as a critical source of water for California, show statistically significant decreases in snowfall and peak snowpack and associated decreases in summer snowmelt and runoff. More generally, shifts in these hydroclimate metrics occur for intermediate elevation basins ranging from 1,315 to 1,455 m (4,314.3–4,773.6 ft), while the warming level of emergence is delayed for lower and higher elevation basins. We also find that several basins already demonstrate significant increases in 5- to 100-year runoff intensities, primarily due to the increasing influence of precipitation falling as rain. Hydroclimate metrics with trends that demonstrate near-term warming levels of emergence are reflected in reanalysis-based observations, suggesting California is entering a fundamentally different hydroclimate regime. While this will likely stress California's water management, the research provided can support when to implement adaptation efforts.

Bass, B. [Univ. of California, Los Angeles, CA (Un↗

Ecohydrological controls on root and microbial respiration in the East River watershed of Colorado

The main objective of this project was to conduct exploratory work to quantify how snow and rain water inputs influence the CO2 coming from the soil surface (soil CO2 flux), and its plant and microbial sources, in the East River watershed, near Crested Butte, Colorado. Knowledge gained from this effort laid the groundwork for a more comprehensive (ongoing) follow-on grant that is using a combination of experiments, field observations, machine learning and simulation modeling to fully disentangle these relationships. New field measurements were made at four locations along an elevational transect on Snodgrass Mountain at Rocky Mountain Biological Laboratory. These sites were chosen to differ in total snowpack and water table depth, and included both deciduous (aspen) and evergreen (spruce/fir) forest types. We used automated measurements of soil CO2 concentrations to quantify the total soil CO2 flux, and the vertical CO2 production within the soil profile at each site. Isotope (radiocarbon, 14C) measurements determined how much of the CO2 emitted from the soil surface came from plant respiration (root metabolism) versus microbial respiration (decomposition of soil organic matter) sources. Supporting data on plant phenology and microbial dynamics provided context for the observed variation in respiration sources. This work was motivated by our overarching hypothesis that quantifying belowground plant and microbial processes separately, and how they are influenced by snow and rain inputs, is necessary for understanding and predicting how the East River watershed ecosystems will respond to future environmental change.

54 ENVIRONMENTAL SCIENCES↗

Climate Dynamics Preceding Summer Forest Fires in California and the Extreme Case of 2018

Recent record-breaking wildfire seasons in California prompt an investigation into the climate patterns that typically precede anomalous summer burned forest area. Using burned-area data from the U.S. Forest Service’s Monitoring Trends in Burn Severity (MTBS) product and climate data from the fifth major global reanalysis produced by the European Centre for Medium-Range Weather Forecasts (ERA5) over 1984–2018, relationships between the interannual variability of antecedent climate anomalies and July California burned area are spatially and temporally characterized. Lag correlations show that antecedent high vapor pressure deficit (VPD), high temperatures, frequent extreme high temperature days, low precipitation, high subsidence, high geopotential height, low soil moisture, and low snowpack and snowmelt anomalies all correlate significantly with July California burned area as far back as the January before the fire season. Seasonal regression maps indicate that a global midlatitude atmospheric wave train in late winter is associated with anomalous July California burned area. July 2018, a year with especially high burned area, was to some extent consistent with the general patterns revealed by the regressions: low winter precipitation and high spring VPD preceded the extreme burned area. However, geopotential height anomaly patterns were distinct from those in the regressions. Extreme July heat likely contributed to the extent of the fires ignited that month, even though extreme July temperatures do not historically significantly correlate with July burned area. While the 2018 antecedent climate conditions were typical of a high-burned-area year, they were not extreme, demonstrating the likely limits of statistical prediction of extreme fire seasons and the need for individual case studies of extreme years.

54 ENVIRONMENTAL SCIENCES↗

Advances in Land Surface Models and Indicators for Drought Monitoring and Prediction

Millions of people across the globe are affected by droughts every year, and recent droughts have highlighted the considerable agricultural impacts and economic costs of these events. Monitoring the state of droughts depends on integrating multiple indicators that each capture particular aspects of hydrologic impact and various types and phases of drought. As the capabilities of land-surface models and remote sensing have improved, important physical processes such as dynamic, interactive vegetation phenology, groundwater, and snow pack evolution now support a range of drought indicators that better reflect coupled water, energy and carbon cycle processes. In this work, we discuss these advances, including newer classes of indicators that can be applied to improve the characterization of drought onset, severity and duration. We utilize a new model-based drought reconstruction to illustrate the role of dynamic phenology and groundwater in drought assessment. Further, through case studies on flash droughts, snow droughts, and drought recovery, we illustrate the potential advantages of advanced model physics and observational capabilities, especially from remote sensing, in characterizing droughts.

Drought↗

Observational Evidence of EHP Effects on the Melting of Snowpack over the Tibetan Plateau

Observational evidences are presented showing that the Indo-Gangetic Plain (IGP) regions, bounded by the high altitude Himalayan mountains, are subject to heavy loading of absorbing aerosols, i.e., black carbon and dust, which can lead to widespread enhancement warming over the Tibetan Plateau and accelerated snowmelt in the western Tibetan Plateau (WTP) and Himalayas. The two pre-monsoon seasons of 2004 and 2005 were strikingly contrasting in terms of the aerosol loading over IGP. The warming of the TP in 2004 relative to 2005 was widespread, covering most of the WTP and Himalayas. This warming is closely linked to patterns of the snow melt. Consistent with the Elevated Heat Pump hypothesis, we find that increased loading of absorbing aerosols over IGP in the pre-monsoon season is associated with increased heating of the upper troposphere by dynamical feedback induced by aerosol heating, and enhances the rate of snowmelt over Himalayas and the WTP in April-May. Composite analysis with more contrasting years also shows that the heating of the troposphere by elevated dust and black carbon aerosols in the boreal spring can lead to widespread enhanced land-atmosphere warming, and accelerated snow melt in the Himalayas and Tibetan Plateau.

Kim, Maeng-Ki↗

Understanding the Hydrogeochemical Response of a Mountainous Watershed Using Integrated Surface‐Subsurface Flow and Reactive Transport Modeling

Abstract Climate change and other disturbances significantly impact hydrogeochemical exports from mountainous headwater catchments such as the Upper Colorado River Basin. Developing a mechanistic understanding of how the physical and chemical processes interact in time and space in an integrated manner is key to quantifying the future impacts of such disturbances. The hydrogeochemical response of a mountainous catchment in the 2010–2019 period is evaluated quantitatively using a high‐resolution model that simulates integrated hydrology, and transport and reactions for selected solutes and minerals. The model assumes that pyrite is present only at depth while calcite is distributed uniformly, and captures the observed C‐Q reasonably well. Distinct C‐Q dynamics are observed in an average (WY16), a wet (WY17), and a dry (WY18) water year. The model also quantifies the water fraction from surface, shallow and deep groundwater compartments using tracers, and suggests greater groundwater contributions to peak stream discharge in the dry WY18. Results demonstrate that calcium concentrations do not change significantly from year to year, while sulfate shows significant temporal variability. Pyrite dissolution is affected by the changing hydrological drivers where it is enhanced in the dry WY18; calcite dissolution supplements calcium dilution under high flow conditions. The model simulates the reaction hotspots controlled by hydrological conditions, and the spatially‐resolved results show that higher soil saturation and less snowpack occur earlier on the south‐facing side than on the north‐facing side. This is a first‐of‐its‐kind demonstration of a model that integrates hydrologic processes, including evapotranspiration, and reactive transport to enable a predictive understanding of hydrogeochemical exports.

54 ENVIRONMENTAL SCIENCES↗