Search NASA⌕ Search

SEARCH · Search NASA

Results for “Surface runoff”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Using Historical Precipitation, Temperature, and Runoff Observations to Evaluate Evaporation Formulations in Land Surface Models

Key to translating soil moisture memory into subseasonal precipitation and air temperature forecast skill is a realistic treatment of evaporation in the forecast system used - in particular, a realistic treatment of how evaporation responds to variations in soil moisture. The inherent soil moisture-evaporation relationships used in today's land surface models (LSMs), however, arguably reflect little more than guesswork given the lack of evaporation and soil moisture data at the spatial scales represented by regional and global models. Here we present a new approach for evaluating this critical aspect of LSMs. Seasonally averaged precipitation is used as a proxy for seasonally-averaged soil moisture, and seasonally-averaged air temperature is used as a proxy for seasonally-averaged evaporation (e.g., more evaporative cooling leads to cooler temperatures)~ the relationship between historical precipitation and temperature measurements accordingly mimics in certain important ways nature's relationship between soil moisture and evaporation. Additional information on the relationship is gleaned from joint analysis of precipitation and streamflow measurements. An experimental framework that utilizes these ideas to guide the development of an improved soil moisture-evaporation relationship is described and demonstrated.

Koster, Randal D.↗
NASA NTRS2020

Unprecedented Atmospheric Conditions (1948–2019) Drive the 2019 Exceptional Melting Season Over the Greenland Ice Sheet

Understanding the role of atmospheric circulation anomalies on the surface mass balance of the Greenland ice sheet (GrIS) is fundamental for improving estimates of its current and future contributions to sea level rise. Here, we show, using a combination of remote sensing observations, regional climate model outputs, reanalysis data, and artificial neural networks, that unprecedented atmospheric conditions (1948–2019) occurring in the summer of 2019 over Greenland promoted new record or close-to-record values of surface mass balance (SMB), runoff, and snowfall. Specifically, runoff in 2019 ranked second within the 1948–2019 period (after 2012) and first in terms of surface mass balance negative anomaly for the hydrological year 1 September 2018–31 August 2019. The summer of 2019 was characterized by an exceptional persistence of anticyclonic conditions that, in conjunction with low albedo associated with reduced snowfall in summer, enhanced the melt–albedo feedback by promoting the absorption of solar radiation and favored advection of warm, moist air along the western portion of the ice sheet towards the north, where the surface melt has been the highest since 1948. The analysis of the frequency of daily 500 hPa geopotential heights obtained from artificial neural networks shows that the total number of days with the five most frequent atmospheric patterns that characterized the summer of 2019 was 5 standard deviations above the 1981–2010 mean, confirming the exceptional nature of the 2019 season over Greenland.

Greenland ice sheet↗
NASA NTRS2000

Refinements to SSiB with an Emphasis on Snow-Physics: Evaluation and Validation Using GSWP and Valdai Data

Refinements to the snow-physics scheme of SSiB (Simplified Simple Biosphere Model) are described and evaluated. The upgrades include a partial redesign of the conceptual architecture to better simulate the diurnal temperature of the snow surface. For a deep snowpack, there are two separate prognostic temperature snow layers - the top layer responds to diurnal fluctuations in the surface forcing, while the deep layer exhibits a slowly varying response. In addition, the use of a very deep soil temperature and a treatment of snow aging with its influence on snow density is parameterized and evaluated. The upgraded snow scheme produces better timing of snow melt in GSWP-style simulations using ISLSCP Initiative I data for 1987-1988 in the Russian Wheat Belt region. To simulate more realistic runoff in regions with high orographic variability, additional improvements are made to SSiB's soil hydrology. These improvements include an orography-based surface runoff scheme as well as interaction with a water table below SSiB's three soil layers. The addition of these parameterizations further help to simulate more realistic runoff and accompanying prognostic soil moisture fields in the GSWP-style simulations. In intercomparisons of the performance of the new snow-physics SSiB with its earlier versions using an 18-year single-site dataset from Valdai Russia, the version of SSiB described in this paper again produces the earliest onset of snow melt. Soil moisture and deep soil temperatures also compare favorably with observations.

Mocko, David M.↗
NASA NTRS1994

Simulating fluxes from heterogeneous land surfaces: Explicit subgrid method employing the biosphere-atmosphere transfer scheme (BATS)

A vectorized version of the biosphere-atmosphere transfer scheme (VBATS) is used to study moisture, energy, and momentum fluxes from heterogeneous land surfaces st the scale of an atmospheric model (AM) grid cells. To incorporate subgrid scale inhomogeneity, VBATS includes two important features: (1) characterization of the land surface (vegetation and soil parameters) at N subgrid points within an AM grid cell and (2) explicit distribution of climate forcing (precipitation, clouds, etc.) over the subgrid. In this study, VBATS is used in stand-alone mode to simulate a single AM grid cell and to evaluate the effects of subgrid scale vegetation and climate specification on the surface fluxes and hydrology. It is found that the partitioning of energy can be affected by up to 30%, runoff by 50%, and surface stress in excess of 60%. Distributing climate forcing over the AM grid cell increases the Bowen ratio, as a result of enhanced sensible heat flux and reduced latent heat flux. The combined effect of heterogeneous vegetation and distribution of climate is found to be dependent on the dominat vegetation class in the AM grid cell. Development of this method is part of a larger program to explore the importance of subgrid scale processes in regional and global climate simulations.

Seth, Anji↗
DOE OSTI2023

Milone eTape Liquid Level Sensor Laboratory Calibration with a Commercial 12.2 cm HS-Flume and 3D Printed 4 cm Micro-HS-Flumes

The NGEE Arctic Rainfall Simulator (NARS) is a variable intensity rainfall simulator (RFS) with a frame design based on the Humphry et al. (2002) RFS and a water delivery system based on the Walnut Gulch (Paige et al., 2004) RFS. The NARS uses an aluminum frame that is fully deconstructable for transportation to field locations and a water system that enables variable rain intensity. Prior to field deployment, H-flume discharge and the corresponding Milone eTape Liquid Level Sensor (eTape) resistance values were measured in the laboratory so that discharge measurements from the NGEE Arctic Rainfall Simulator could be automated. eTape resistance was measured with increasing fluid heights for the full range of the eTape to calculate sensor detection limits and resolution. The eTape was calibrated to both a commercially available 12.2 cm HS-flume and a 3D printed 4 cm Micro-HS-flume. This data package contains two .csv files, one for the eTape calibration and the other for the flume calibrations, and two .stl files to 3D print the 4 cm Micro-HS-Flume design with a 1 cm x 3.6 cm stilling well opening. The .stl files can be opened in most 3D printing or CAD programs/software. Two .kml files of the laboratory location and broader area where testing was conducted are also included in this data package.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Data and scripts associated with the manuscript evaluating the hydrologic responses of the Pacific Northwest watersheds to wildfires (v2)

This data package is associated with the publication “Evaluating Post-fire Watershed Response to Varying Burn Severity and Precipitation Regimes Using Fully-distributed and Integrated Hydrologic Models” submitted to Journal of Hydrology (Li et al. 2025). In this study, we employed the Advanced Terrestrial Simulator (ATS), an integrated watershed model that couples surface flow, subsurface flow, and canopy biophysical processes, to investigate post-fire hydrologic responses in a few selected watersheds with varying burn severity.The data package contains the required input data (meteorological forcing, Leaf Area Index, wildfire burn severities, etc.) to run the model, configuration files, the Jupyter notebooks in Python to pre-process and post-process data, the figures in the manuscript, and the modeling output files. The variables include watershed-averaged evapotranspiration, watershed-averaged surface/subsurface/canopy water content, and river discharge at watershed outlet.The data package contains a file-level metadata that lists and describes all the files contained in the data package (ATS_flmd.csv), a data dictionary file that defines columns headers across all csv files contained in the data package (ATS_dd.csv), a data package level readme file (the current file), and four zipped folders.The ‘data’ folder provides data needed to run the model in .h5, .i2s, .xyz, .shp, and .exo formats. The sub-folders are for each data types. The ‘model’ folder provides input files (.xml format) and essential model outputs. Each sub-folder provides the files from each simulated watershed. The ‘notebooks’ folder provides the Jupyter notebooks (.ipynb format) for pre- and post- processing model files, and for producing the figures in the manuscript. The ‘figures’ folder provides the figures associated with manuscript in .pdf and .png formats.The ‘model’ folder and the ‘data’ folder have been split into 5GB-large pieces using the Linux command ‘split -b 5120m model.zip model.zip.’ and ‘split -b 5120m data.zip data.zip.’, respectively. They can be merged back using the Linux command ‘cat model.zip.* > model.zip’ and ‘cat data.zip.* > data.zip’, respectively.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2008

Ongoing Development of NASA's Global Land Data Assimilation System

NASA's Global Land Data Assimilation System (GLDAS) produces global fields of land surface states (e.g., soil moisture and temperature) and fluxes (e.g., latent heat flux and runoff) by driving offline land surface models with observation-based inputs, using the Land Information System (LIS) software. Since production began in 2001, GLDAS has supported more than 100 scientific investigations and applications. Some examples are GEWEX and NASA Energy and Water Cycle Study (NEWS) global water and energy budget analyses, interpretations of hydrologic data derived from the Gravity Recovery and Climate Experiment (GRACE) satellite mission, and forecast model initiation studies at NOAA and NASA. At the same time, the GLDAS team has continued improve results through the development of new modeling and data assimilation techniques. Here we describe several recent and ongoing innovations. These include global implementation of a runoff routing procedure, GRACE data assimilation, advanced snow cover assimilation, and irrigation modeling.

Rodell, Matthew↗
DOE OSTI2026

Runoff evaluation in an Earth System Land Model for permafrost regions in Alaska

Modeling of hydrological runoff is essential for accurately capturing spatiotemporal feedbacks within the land–atmosphere system, particularly in sensitive regions such as permafrost landscapes. However, substantial uncertainties persist in the terrestrial runoff parameterization schemes used in Earth system and land surface models. This is particularly true in permafrost regions, where landscape heterogeneity is high and reliable observational data are scarce. In this study, we evaluate the performance of runoff parameterization schemes in the Energy Exascale Earth System Model (E3SM) land model (ELM). Our proposed framework leverages simulation results from the Advanced Terrestrial Simulator (ATS), which is a physics-based integrated surface/subsurface hydrologic model that has been successfully evaluated previously in Arctic tundra regions. We used ATS to simulate runoff from 22 representative hillslopes in the Sagavanirktok River basin, located on the North Slope of Alaska, then compared the output with ELM's parameterized representation of total runoff. Results show that (1) ELM's total runoff was the same order of magnitude as the ATS simulations, and both models were similarly variable over time; (2) minor adjustments to coefficients in ELM's runoff parameterization improved the match between the ATS simulation and ELM's parameterized representation of annual and seasonal total runoff; (3) overall, runoff responses in ATS and ELM are more similar in flat hillslope environments compared to steep hillslopes; and (4) shallower active layer thicknesses and higher precipitation simulations resulted in lower correlations between the two models due to greater total runoff. By incorporating the optimized runoff coefficients from the Sagavanirktok River basin into ELM, the simulated total runoff better matched the streamflow observations at a small watershed located on the Seward Peninsula of Alaska. Our findings revealed important insights into the effectiveness of runoff parameterizations in land surface models and pathways for improving runoff coefficients in typical Arctic regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Efficacy of underground aggregate infiltration beds under a permeable pavement system

Urban areas often face significant flooding challenges due to impermeable surfaces that exacerbate stormwater runoff, making effective stormwater management essential. This study demonstrated that combining pervious concrete pavement systems with underground aggregate storage beds and rain gardens can significantly enhance water infiltration and drainage efficiency, providing valuable insights for designing low-impact development strategies in flood-prone regions.

Haselbach, Liv [Lamar Univ., Beaumont, TX (United ↗
DOE OSTI2025

Integrating Contaminant Source Indicators, Water Quality Measures, and Ecotoxicity to Characterize Contaminant Mixtures and Per- and Polyfluoroalkyl Substance (PFAS) Variability in an Urban Watershed

Thousands of chemical contaminants threaten watersheds but are time and cost prohibitive to monitor. Identifying their sources, transport, and ecological risk is limited in heterogeneous urban watersheds. We present an integrative watershed approach using source-specific indicator compounds, common water quality measures, and ecotoxicity assays to examine the distribution of contaminant mixtures in an urbanized watershed. Indicator compound concentrations were temporally and spatially distributed for treated/untreated sewage (sucralose, artificial sweetener), road runoff (diphenyl-guanidine [DPG] and 6PPD-quinone [6PPD-Q], automobile tire additives), and lawncare runoff (aminomethanephosphonic acid (AMPA), major degradant of the herbicide glyphosate). Sucralose was predominately sourced from treated wastewater; measurable concentrations in tributaries indicated raw sewage inputs. DPG and 6PPD-Q concentrations correlated to road density during base flow and were elevated during stormflow. AMPA was measurable spring through fall, especially where lawns were dense. When specific sources dominated flow, water quality measures correlated with wastewater (sulfate, potassium, chloride, and sodium) and road runoff (chromium and lead) indicators. The limited behavioral toxicity observed in exposed zebrafish (Danio rerio) (18%) was not well explained by source-indicators. PFAS concentrations were highly variable spatially but not well explained by our source-specific indicator compounds. Here, more costly compound-specific monitoring may be necessary when multiple sources exist or when unexpected toxicity trends occur.

computer simulations↗
NASA NTRS1994

Statistical analysis of mesoscale rainfall: Dependence of a random cascade generator on large-scale forcing

Under the theory of independent and identically distributed random cascades, the probability distribution of the cascade generator determines the spatial and the ensemble properties of spatial rainfall. Three sets of radar-derived rainfall data in space and time are analyzed to estimate the probability distribution of the generator. A detailed comparison between instantaneous scans of spatial rainfall and simulated cascades using the scaling properties of the marginal moments is carried out. This comparison highlights important similarities and differences between the data and the random cascade theory. Differences are quantified and measured for the three datasets. Evidence is presented to show that the scaling properties of the rainfall can be captured to the first order by a random cascade with a single parameter. The dependence of this parameter on forcing by the large-scale meteorological conditions, as measured by the large-scale spatial average rain rate, is investigated for these three datasets. The data show that this dependence can be captured by a one-to-one function. Since the large-scale average rain rate can be diagnosed from the large-scale dynamics, this relationship demonstrates an important linkage between the large-scale atmospheric dynamics and the statistical cascade theory of mesoscale rainfall. Potential application of this research to parameterization of runoff from the land surface and regional flood frequency analysis is briefly discussed, and open problems for further research are presented.

Over, Thomas, M.↗
NASA NTRS2022

Assessing Bare-Ice Albedo Simulated By MAR Over the Greenland Ice Sheet (2000–2021) and Implications for Meltwater Production Estimates

Surface mass loss from the Greenland ice sheet (GrIS) has accelerated over the past decades, mainly due to enhanced surface melting and liquid water runoff in response to atmospheric warming. A large portion of runoff from the GrIS originates from exposure of the darker bare ice in the ablation zone when the overlying snow melts, where surface albedo plays a critical role in modulating the energy available for melting. In this regard, it is imperative to understand the processes governing albedo variability to accurately project future mass loss from the GrIS. Bare-ice albedo is spatially and temporally variable and contingent on non-linear feedbacks and the presence of light-absorbing constituents. An assessment of models aiming at simulating albedo variability and associated impacts on meltwater production is crucial for improving our understanding of the processes governing these feedbacks and, in turn, surface mass loss from Greenland. Here, we report the results of a comparison of the bare-ice extent and albedo simulated by the regional climate model Modèle Atmosphérique Régional (MAR) with satellite imagery from the Moderate Resolution Imaging Spectroradiometer (MODIS) for the GrIS below 70∘ N. Our findings suggest that MAR overestimates bare-ice albedo by 22.8 % on average in this area during the 2000–2021 period with respect to the estimates obtained from MODIS. Using an energy balance model to parameterize meltwater production, we find this bare-ice albedo bias can lead to an underestimation of total meltwater production from the bare-ice zone below 70∘ N of 42.8 % during the summers of 2000–2021.

albedo↗
DOE OSTI2024

Modeled estimates of permafrost hydrology and related fields for pan-Arctic region over the period 1980-2100

This dataset contains gridded estimates of total runoff (mm/day); subsurface runoff (mm/day); snow water equivalent (SWE, mm/day); vertical upward water flux of combined surface evaporation, evapotranspiration, and snow sublimation (mm/day); and soil temperature (C) across the pan-Arctic basin from simulations of the Permafrost Water Balance Model (PWBM v4). These data were generated to allow researchers to advance understanding of the changes unfolding in the Arctic's terrestrial water cycle under the influences of climate warming, hydrological cycle intensification, and permafrost thaw. The dataset will enable researchers to address important questions spanning multiple disciplines. For example, how will rising river flows exported to the coast impact critical zone ecosystems, near-shore ice dynamics, and ocean biogeochemical processes. Spatial resolution is 25x25 km with N=39926 grid cells spanning the domain. Three forcing datasets were used to drive simulations at an implicit daily time step over the period 1980-2100. Daily air temperature, precipitation, and wind speed data from the WATCH Forcing Data methodology applied to ERA5 reanalysis (W5E5) were used for meteorological forcings over the period 1980-2019. Simulations for the period 1980-2100 used meteorological data from two global climate models; IPSL-CM6A-LR and MPI-ESM1-2-HR. The PWBM hydrological outputs are organized by grid and month in gzipped ASCII text files within each archived tar file. Soil temperature is organized by grid cell and day of year from DOY 205 to DOY 260. No special software is required to work with these data. The package contains spatial domain and river basin information files. These data are most useful for analyses of the spatial and temporal variations in terrestrial hydrology and soil thermal dynamics. A manuscript describing the model simulations and associated analysis has been accepted for publication in The Cryosphere (Rawlins and Karmalkar, 2024).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

An Update on the Hydrological Land Surface Data and Services at NASA GES DISC

"The NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) is one of twelve NASA Earth Observing System data centers that document, process, archive and distribute data from Earth science missions and related projects. The GES DISC hosts many hydrological land surface data products including North American Land Data Assimilation System (NLDAS), Global Land Data Assimilation System (GLDAS), the Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System (FLDAS), and assimilated-Gravity Recovery and Climate Experiment (GRACE). The Land Surface Model (LSM) products contain model output, including heat fluxes, rain, snow, soil temperature, soil moisture and runoff; and observational forcing data, including surface pressure, temperature, precipitation, downward shortwave and longwave radiation, humidity, and wind. The temporal resolution of the hydrology data at the GES DISC ranges from hourly to monthly and the spatial resolution 0.1 degree to 1.0 degree. The variability of Earth’s water and energy cycles can be characterized with these spatially and temporally consistent, and quality controlled hydrological land surface data. These data could support ongoing modeling activities, and help improve our understanding of land-surface-atmosphere interactions and their impact on climate. Over the years, new data is added, and older data is reprocessed and updated. This presentation provides a summary table of the hydrological data and discusses recent updates to the data and data services at the GES DISC. New data products include (1) GRACE Global Version 3.0 drought indicators and (2) FLDAS monthly global and Central Asia daily data. Reprocessed and updated products include (1) NLDAS Version 2.0 and (2) GRACE United States Version 4.0 drought indicators. Lastly, this presentation will discuss the opportunities for accessing some of our hydrology data from the cloud. As the hydrology land surface data is moved to the cloud, users can expect faster service for longer time series, notably the hydrology data rods time series service."

Ashley Heath↗
NASA NTRSDate not supplied

Relating Land Cover Change to Runoff Distribution Using NASA Earth Observations in Riley County, Kansas

In recent years, Riley County, Kansas experienced high levels of flooding potentially due to changes in land use and impervious surface. To analyze the increase in floods and runoff patterns over time, this study contrasted two methods of curve number (CN) calculation, which estimated runoff from rainfall events, from 2006-2020. The conventional calculation method used a look-up table and tracked land use/land cover change. This method resulted in static tables from various inputs such as soil type, surface imperviousness, slope, and crop phenology to estimate runoff versus infiltration. These tables allowed for phenology specific curve numbers and accounted for various farming techniques but required many inputs and lacked flexibility to include seasonal variations in vegetation. In contrast, the dynamic method is a more modern and simple calculation employing normalized difference vegetation index (NDVI) data from Landsat 5, Landsat 7, and Landsat 8 satellites compiled over the rainy season each year to calculate curve numbers based off seasonal vegetation. This dynamic method allowed for more precise analysis of runoff variability because NDVI is a more accurate, real-time measurement for land cover and it could be calculated with greater temporal detail. The results will allow decision makers to make more informed decisions on resiliency strategies to address future flooding.

Trista Brophy↗
DOE OSTI2023

High‐Latitude Stratospheric Aerosol Injection to Preserve the Arctic

Abstract Stratospheric aerosol injection (SAI) has been shown in climate models to reduce some impacts of global warming in the Arctic, including the loss of sea ice, permafrost thaw, and reduction of Greenland Ice Sheet (GrIS) mass; SAI at high latitudes could preferentially target these impacts. In this study, we use the Community Earth System Model to simulate two Arctic‐focused SAI strategies, which inject at 60°N latitude each spring with injection rates adjusted to either maintain September Arctic sea ice at 2030 levels (“Arctic Low”) or restore it to 2010 levels (“Arctic High”). Both simulations maintain or restore September sea ice to within 10% of their respective targets, reduce permafrost thaw, and increase GrIS surface mass balance by reducing runoff. Arctic High reduces these impacts more effectively than a globally focused SAI strategy that injects similar quantities of SO 2 at lower latitudes. However, Arctic‐focused SAI is not merely a “reset button” for the Arctic climate, but brings about a novel climate state, including changes to the seasonal cycles of Northern Hemisphere temperature and sea ice and less high‐latitude carbon uptake relative to SSP2‐4.5. Additionally, while Arctic‐focused SAI produces the most cooling near the pole, its effects are not confined to the Arctic, including detectable cooling throughout most of the northern hemisphere for both simulations, increased mid‐latitude sulfur deposition, and a southward shift of the location of the Intertropical Convergence Zone. For these reasons, it would be incorrect to consider Arctic‐focused SAI as “local” geoengineering, even when compared to a globally focused strategy.

Lee, Walker Raymond↗
DOE OSTI2023

Detecting Permafrost Active Layer Thickness Change From Nonlinear Baseflow Recession

We report permafrost underlies about one fifth of the global land area and affects ground stability, freshwater runoff, soil chemistry, and surface-atmosphere gas exchange. The depth of thawed ground overlying permafrost (active layer thickness) has broadly increased across the Arctic in recent decades, coincident with a period of increased streamflow, especially the lowest flows (baseflow). Mechanistic links between active layer thickness and baseflow have recently been explored using linear reservoir theory, but most watersheds behave as nonlinear reservoirs. We derive theoretical nonlinear relationships between long-term average saturated soil thickness $\bar{η}$(proxy for active layer thickness) and long-term average baseflow. When applied to 38 years of daily streamflow data for the Kuparuk River basin on the North Slope of Alaska, the theory predicts $\bar{η}$ increased 0.17 ± 0.22[2σ] cm a -1 between 1983 and 2020 (6.4 ± cm total). The rate of increase nearly doubled to 0.29 ± 0.31 cm a -1 between 1990 and 2020, during which time local field measurements from Circumpolar Active Layer Monitoring sites indicate the active layer increased 0.31 ± 0.22 cm a -1 . The predicted rate of increase more than doubled again between 2002 and 2020, outpacing a near doubling of observed active layer thickening, consistent with trends in terrestrial water storage inferred from Gravity Recovery and Climate Experiment satellite gravimetry and Modern-Era Retrospective Analysis for Research and Applications climate reanalysis. Overall, hydrologic change is accelerating in the Kuparuk River basin, and we provide a theoretical framework for estimating basin-scale changes in active layer water storage from streamflow measurements.

54 ENVIRONMENTAL SCIENCES↗