SEARCH · Geology Reports

Results for “Input”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

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

At least 1,315 records · Page 73Linked to original sources

USGS2024

A case study and comparison of mid-winter warming and solar driven wet slab avalanche cycles

Wet slab avalanches are poorly understood and often difficult to forecast. Yet, wet slab avalanches can be destructive and may become more common in a changing climate. As the onset of wet avalanches moves earlier in the winter season due to climate change, understanding snowpack and meteorological characteristics of wet slab avalanches will become increasingly important. In this study, we examined two recent late-January wet slab cycles triggered by warming and solar input in the Rocky Mountains of Idaho and Montana, United States. We posed two questions to help us understand these potentially increasingly frequent, mid-winter wet slab cycles: i) what are the weather and snowpack patterns of two mid-winter wet slab cycles and ii) how did these non-rain-on-snow, mid-winter warming events compare to the historical climate normal? In both locations, a short-lived ridge of high pressure strengthened over the region from late January to early February 2024. Temperatures rapidly increased to 6° C at middle and upper-elevation locations and remained above freezing for 48 to 72 hours. Wet slab avalanches (Montana: n=68, Idaho: n=11) released on a layer of facets resting on a crust formed in late December and slab depths averaged 40 to 150 cm. Avalanches primarily occurred on southeast, south, and southwest aspects at middle and upper elevations. Additionally, positive January monthly temperature trends exist in both locations from 1990 to 2020, and these two cycles from 2024 highlight impacts of midwinter warming events with unstable snowpack conditions. These cycles also highlight how a relatively uncommon mid-winter warming event with low radiative input can produce very large destructive wet slab avalanches. Understanding mid-winter warming events and associated wet slab avalanche cycles will help us prepare and forecast for potentially more common scenarios like this in the future.

Idaho, Montana↗

USGS2006

Dealing with uncertainty and sensitivity issues in process-based models of carbon and nitrogen cycles in northern forest ecosystems

Many process-based models on carbon (C) and nitrogen (N) cycles have been developed for northern forest ecosystems. These models are widely used to evaluate the long-term decisions in forest management dealing with effects like particulate pollution, productivity and climate change. Regarding climate change, one of the key questions that have sensitive political implications is whether northern forests will sequester atmospheric C or not. Whilst many process-based models have been tested for accuracy by evaluating or validating against observed data, few have dealt with the complexity of the incorporated procedures to estimate uncertainties associated with model predictions or the sensitivity of these predictions to input factors in a systematic, inter-model comparison fashion. In general, models differ in their underlying attempts to match natural complexities with assumed or imposed model structure and process formulations to estimate model parameters, to gather data and to address issues on scope, scale and natural variations. Uncertainties may originate from model structure, estimation of model parameters, data input, representation of natural variation and scaling exercises. Model structure relates to the mathematical representation of the processes modelled and the type of state variables that a model contains. The modelling of partitioning among above- and below-ground C and N pools and the interdependence among these pools remain a major source of uncertainty in model structure and error propagation. Most soil C models use at least three state variables to represent the different types of soil organic matter (SOM). This approach results in creating three artificial SOM pools, assuming that each one contains C compounds with same turnover rate. In reality, SOM consists of many different types of C compounds with widely different turnover rates. Uncertainty in data and parameter estimates are closely linked. Data uncertainties are associated with high variations in estimating forest biomass, productivity and soil organic matter and may be incomplete for model initialization, calibration, validation and sensitivity analysis of generalized predictor models. The scale at which a model is being used also affects the level of uncertainty, as the errors in the prediction of the C and N dynamics differ from the site to the landscape levels and across climatic regions. If the spatial or temporal scale of a model application is changed, additional uncertainty arises from neglecting natural variability in system variables in time and space. Uncertainty issues are also intimately related to model validation and sensitivity analysis. The estimation of uncertainties is needed to inform decision process, in order to detect the possible corridor of development. Uncertainty in this context is an essential measure of quality for stakeholder and decision makers.

Conference Paper↗

USGS1996

Trends in nutrients

The roles of nitrogen (N) and phosphorus (P) as key nutrients determining the trophic status of water bodies are examined, and evidence reviewed for trends in concentrations of N and P species which occur in freshwaters, primarily in northern temperate environments. Data are reported for water bodies undergoing eutrophication and acidification, especially water bodies receiving increased nitrogen inputs through the atmospheric deposition of nitrogen oxides (NO x ). Nutrient loading on groundwaters and surface freshwaters is assessed with respect to causes and rates of change, relative rates of change for N and P, and implications of change for the future management of lakes, rivers and groundwaters. In particular, the nature and emphasis of studies for N species and P fractions in lakes versus rivers and groundwaters are contrasted. This review paper primarily focuses on results from North America and Europe, particularly for the UK where a wide range of data sets exists. Few nutrient loading data have been published on water bodies in less developed countries; however, some of the available data are presented to provide a global perspective. In general, N and P concentrations have increased dramatically (>20 times background concentrations) in many areas and causes vary considerably, ranging from urbanization to changes in agricultural practices.

Hydrological Processes↗

USGS1997

Disturbance of biological soil crusts: Impacts on potential wind erodibility of sandy desert soils in southeastern Utah

Friction threshold velocities (FTVs) were determined for biological soil crusts in different stages of recovery. Particles on the surface of crusts that had been relatively undisturbed for at least 20 years were found to have significantly higher FTVs than those that had been disturbed 5, 10 or 1 years previously (376, 87, and 46 cm sec -1 , respectively). FTV's for crust breakage was also much higher for undisturbed crusts when compared to the previously disturbed crusts (573, 148, and 88 cm sec -1 , respectively). All crusted surfaces were more stable than bare sand, which had an FTV of 16 cm sec -1 . Disturbance treatments were then applied to the three crustal classes. Disturbance significantly reduced the FTVs of all classes by 73±92 per cent. Comparing crustal FTVs with mean and high monthly wind speeds found in this region, it was observed that only crusts that had been undisturbed for approximately 20 years or more were able to protect soil surfaces from wind gusts expected on the average of once a month. Other crustal classes, as well as all disturbance treatments, had FTVs lower or equal to that of commonly occurring winds in this region. Because most of the crustal biomass occurs in the top 0-3 mm of soils, even slight soil loss can negatively influence stability and nutrient inputs to this ecosystem.

Utah↗

USGS2006

Unsaturated zone flow processes

Water flow in the unsaturated zone is greatly influenced by unsaturated hydrostatics (water content, energy, pressure, and retention) and by unsaturated hydrodynamics (diffuse flow and preferential flow). Important multiphase processes include the transport of gases, nonaqueous liquids, and solid particles. Numerous means are available for determination of unsaturated conditions and properties, both measurement (of moisture state, water retention, and dynamic characteristics) and through various formulas and models that are mostly empirical in nature, but in some cases incorporating insight into unsaturated‐zone physical processes. Applications to practical problems include models and techniques relating to distributions of water and energy, fluxes at the land surface, inputs, outputs, and fluxes within the unsaturated zone, all of which are frequently complicated by heterogeneity and preferential flow. Further scientific advance requires new measurement techniques and theoretical constructs that more adequately represent the important physical processes within practical modeling schemes.

Book chapter↗

USGS2001

Erosion and sediment delivery following removal of forest roads

Erosion control treatments were applied to abandoned logging roads in California, with the goal of reducing road-related sediment input to streams and restoring natural hydrologic patterns on the landscape. Treatment of stream crossings involved excavating culverts and associated road fill and reshaping streambanks. A variety of techniques were applied to road benches, which included decompacting the road surface, placing unstable road fill in more stable locations, and re-establishing natural surface drainage patterns. Following treatment and a 12-year recurrence-interval storm, some road reaches and excavated stream crossings showed evidence of mass movement failures, gullying, bank erosion and channel incision. Post-treatment erosion from excavated stream crossings was related to two variables: a surrogate for stream power (drainage area × channel gradient) and the volume of fill excavated from the channel. Post-treatment erosion on road reaches was related to four explanatory variables: method of treatment, hillslope position (upper, mid-slope or lower), date of treatment, and an interaction term (hillslope position × method of treatment). Sediment delivery from treated roads in upper, middle and lower hillslope positions was 10, 135 and 550 m 3 of sediment per kilometre of treated roads, respectively. In contrast, inventories of almost 500 km of forest roads in adjacent catchments indicate that untreated roads produced 1500 to 4700 m 3 of sediment per kilometre of road length. Erosion from 300 km of treated roads contributed less than 2 per cent of the total sediment load of Redwood Creek during the period 1978 to 1998. Although road removal treatments do not completely eliminate erosion associated with forest roads, they do substantially reduce sediment yields from abandoned logging roads.

Earth Surface Processes and Landforms↗

USGS2014

CO 2 and CH 4 emissions from streams in a lake-rich landscape: Patterns, controls, and regional significance

Aquatic ecosystems are important components of landscape carbon budgets. In lake-rich landscapes, both lakes and streams may be important sources of carbon gases (CO 2 and CH 4 ) to the atmosphere, but the processes that control gas concentrations and emissions in these interconnected landscapes have not been adequately addressed. We use multiple data sets that vary in their spatial and temporal extent during 2001–2012 to investigate the carbon gas source strength of streams in a lake-rich landscape and to determine the contribution of lakes, metabolism, and groundwater to stream CO 2 and CH 4 . We show that streams emit roughly the same mass of CO 2 (23.4 Gg C yr −1 ; 0.49 mol CO 2 m −2 d −1 ) as lakes at a regional scale (27 Gg C yr −1 ) and that stream CH 4 emissions (189 Mg C yr −1 ; 8.46 mmol CH 4 m −2 d −1 ) are an important component of the regional greenhouse gas balance. Gas transfer velocity variability (range = 0.34 to 13.5 m d −1 ) contributed to the variability of gas flux in this landscape. Groundwater inputs and in-stream metabolism control stream gas supersaturation at the landscape scale, while carbon cycling in lakes and deep groundwaters does not control downstream gas emissions. Our results indicate the need to consider connectivity of all aquatic ecosystems (lakes, streams, wetlands, and groundwater) in lake-rich landscapes and their connections with the terrestrial environment in order to understand the full nature of the carbon cycle.

Global Biogeochemical Cycles↗

USGS2014

Influences of glacial melt and permafrost thaw on the age of dissolved organic carbon in the Yukon River basin

Responses of near-surface permafrost and glacial ice to climate change are of particular significance for understanding long-term effects on global carbon cycling and carbon export by high-latitude northern rivers. Here we report Δ 14 C-dissolved organic carbon (DOC) values and dissolved organic matter optical data for the Yukon River, 15 tributaries of the Yukon River, glacial meltwater, and groundwater and soil water end-member sources draining to the Yukon River, with the goal of assessing mobilization of aged DOC within the watershed. Ancient DOC was associated with glacial meltwater and groundwater sources. In contrast, DOC from watersheds dominated by peat soils and underlain by permafrost was typically enriched in Δ 14 C indicating that degradation of ancient carbon stores is currently not occurring at large enough scales to quantitatively influence bulk DOC exports from those landscapes. On an annual basis, DOC exported was predominantly modern during the spring period throughout the Yukon River basin and became older through summer-fall and winter periods, suggesting that contributions of older DOC from soils, glacial meltwaters, and groundwater are significant during these months. Our data indicate that rapidly receding glaciers and increasing groundwater inputs will likely result in greater contributions of older DOC in the Yukon River and its tributaries in coming decades.

Alaska↗

USGS2014

Differential response of carbon fluxes to climate in three peatland ecosystems that vary in the presence and stability of permafrost

Changes in vegetation and soil properties following permafrost degradation and thermokarst development in peatlands may cause changes in net carbon storage. To better understand these dynamics, we established three sites in Alaska that vary in permafrost regime, including a black spruce peat plateau forest with stable permafrost, an internal collapse scar bog formed as a result of thermokarst, and a rich fen without permafrost. Measurements include year-round eddy covariance estimates of carbon dioxide (CO 2 ), water, and energy fluxes, associated environmental variables, and methane (CH 4 ) fluxes at the collapse scar bog. The ecosystems all acted as net sinks of CO 2 in 2011 and 2012, when air temperature and precipitation remained near long-term means. In 2013, under a late snowmelt and late leaf out followed by a hot, dry summer, the permafrost forest and collapse scar bog were sources of CO 2 . In this same year, CO 2 uptake in the fen increased, largely because summer inundation from groundwater inputs suppressed ecosystem respiration. CO 2 exchange in the permafrost forest and collapse scar bog was sensitive to warm air temperatures, with 0.5 g C m −2 lost each day when maximum air temperature was very warm (≥29°C). The bog lost 4981 ± 300 mg CH 4 m −2 between April and September 2013, indicating that this ecosystem acted as a significant source of both CO 2 and CH 4 to the atmosphere in 2013. These results suggest that boreal peatland responses to warming and drying, both of which are expected to occur in a changing climate, will depend on permafrost regime.

Alaska↗

USGS2014

Uncertainty in the fate of soil organic carbon: A comparison of three conceptually different soil decomposition models

Conventional Q10 soil organic matter decomposition models and more complex microbial models are available for making projections of future soil carbon dynamics. However, it is unclear (1) how well the conceptually different approaches can simulate observed decomposition and (2) to what extent the trajectories of long-term simulations differ when using the different approaches. In this study, we compared three structurally different soil carbon (C) decomposition models (one Q10 and two microbial models of different complexity), each with a one- and two-horizon version. The models were calibrated and validated using 4 years of measurements of heterotrophic soil CO 2 efflux from trenched plots in a Dahurian larch ( Larix gmelinii Rupr.) plantation. All models reproduced the observed heterotrophic component of soil CO 2 efflux, but the trajectories of soil carbon dynamics differed substantially in 100 year simulations with and without warming and increased litterfall input, with microbial models that produced better agreement with observed changes in soil organic C in long-term warming experiments. Our results also suggest that both constant and varying carbon use efficiency are plausible when modeling future decomposition dynamics and that the use of a short-term (e.g., a few years) period of measurement is insufficient to adequately constrain model parameters that represent long-term responses of microbial thermal adaption. These results highlight the need to reframe the representation of decomposition models and to constrain parameters with long-term observations and multiple data streams. We urge caution in interpreting future soil carbon responses derived from existing decomposition models because both conceptual and parameter uncertainties are substantial.

Journal of Geophysical Research G: Biogeosciences↗

USGS2014

Lattice Boltzmann methods applied to large-scale three-dimensional virtual cores constructed from digital optical borehole images of the karst carbonate Biscayne aquifer in southeastern Florida

Digital optical borehole images at approximately 2 mm vertical resolution and borehole caliper data were used to create three-dimensional renderings of the distribution of (1) matrix porosity and (2) vuggy megaporosity for the karst carbonate Biscayne aquifer in southeastern Florida. The renderings based on the borehole data were used as input into Lattice Boltzmann methods to obtain intrinsic permeability estimates for this extremely transmissive aquifer, where traditional aquifer test methods may fail due to very small drawdowns and non-Darcian flow that can reduce apparent hydraulic conductivity. Variogram analysis of the borehole data suggests a nearly isotropic rock structure at lag lengths up to the nominal borehole diameter. A strong correlation between the diameter of the borehole and the presence of vuggy megaporosity in the data set led to a bias in the variogram where the computed horizontal spatial autocorrelation is strong at lag distances greater than the nominal borehole size. Lattice Boltzmann simulation of flow across a 0.4 × 0.4 × 17 m (2.72 m 3 volume) parallel-walled column of rendered matrix and vuggy megaporosity indicates a high hydraulic conductivity of 53 m s −1 . This value is similar to previous Lattice Boltzmann calculations of hydraulic conductivity in smaller limestone samples of the Biscayne aquifer. The development of simulation methods that reproduce dual-porosity systems with higher resolution and fidelity and that consider flow through horizontally longer renderings could provide improved estimates of the hydraulic conductivity and help to address questions about the importance of scale.

Water Resources Research↗

USGS2014

Accuracy of travel time distribution (TTD) models as affected by TTD complexity, observation errors, and model and tracer selection

Analytical models of the travel time distribution (TTD) from a source area to a sample location are often used to estimate groundwater ages and solute concentration trends. The accuracies of these models are not well known for geologically complex aquifers. In this study, synthetic datasets were used to quantify the accuracy of four analytical TTD models as affected by TTD complexity, observation errors, model selection, and tracer selection. Synthetic TTDs and tracer data were generated from existing numerical models with complex hydrofacies distributions for one public-supply well and 14 monitoring wells in the Central Valley, California. Analytical TTD models were calibrated to synthetic tracer data, and prediction errors were determined for estimates of TTDs and conservative tracer (NO 3 −) concentrations. Analytical models included a new, scale-dependent dispersivity model (SDM) for two-dimensional transport from the watertable to a well, and three other established analytical models. The relative influence of the error sources (TTD complexity, observation error, model selection, and tracer selection) depended on the type of prediction. Geological complexity gave rise to complex TTDs in monitoring wells that strongly affected errors of the estimated TTDs. However, prediction errors for NO3− and median age depended more on tracer concentration errors. The SDM tended to give the most accurate estimates of the vertical velocity and other predictions, although TTD model selection had minor effects overall. Adding tracers improved predictions if the new tracers had different input histories. Studies using TTD models should focus on the factors that most strongly affect the desired predictions.

California↗

USGS2016

Toward more realistic projections of soil carbon dynamics by Earth system models

Soil carbon (C) is a critical component of Earth system models (ESMs), and its diverse representations are a major source of the large spread across models in the terrestrial C sink from the third to fifth assessment reports of the Intergovernmental Panel on Climate Change (IPCC). Improving soil C projections is of a high priority for Earth system modeling in the future IPCC and other assessments. To achieve this goal, we suggest that (1) model structures should reflect real-world processes, (2) parameters should be calibrated to match model outputs with observations, and (3) external forcing variables should accurately prescribe the environmental conditions that soils experience. First, most soil C cycle models simulate C input from litter production and C release through decomposition. The latter process has traditionally been represented by first-order decay functions, regulated primarily by temperature, moisture, litter quality, and soil texture. While this formulation well captures macroscopic soil organic C (SOC) dynamics, better understanding is needed of their underlying mechanisms as related to microbial processes, depth-dependent environmental controls, and other processes that strongly affect soil C dynamics. Second, incomplete use of observations in model parameterization is a major cause of bias in soil C projections from ESMs. Optimal parameter calibration with both pool- and flux-based data sets through data assimilation is among the highest priorities for near-term research to reduce biases among ESMs. Third, external variables are represented inconsistently among ESMs, leading to differences in modeled soil C dynamics. We recommend the implementation of traceability analyses to identify how external variables and model parameterizations influence SOC dynamics in different ESMs. Overall, projections of the terrestrial C sink can be substantially improved when reliable data sets are available to select the most representative model structure, constrain parameters, and prescribe forcing fields.

Global Biogeochemical Cycles↗

USGS2015

Abundance, stable isotopic composition, and export fluxes of DOC, POC, and DIC from the Lower Mississippi River during 2006–2008

Sources, abundance, isotopic compositions, and export fluxes of dissolved inorganic carbon (DIC), dissolved and colloidal organic carbon (DOC and COC), and particulate organic carbon (POC), and their response to hydrologic regimes were examined through monthly sampling from the Lower Mississippi River during 2006–2008. DIC was the most abundant carbon species, followed by POC and DOC. Concentration and δ 13 C of DIC decreased with increasing river discharge, while those of DOC remained fairly stable. COC comprised 61 ± 3% of the bulk DOC with similar δ 13 C abundances but higher percentages of hydrophobic organic acids than DOC, suggesting its aromatic and diagenetically younger status. POC showed peak concentrations during medium flooding events and at the rising limb of large flooding events. While δ 13 C-POC increased, δ 15 N of particulate nitrogen decreased with increasing discharge. Overall, the differences in δ 13 C between DOC or DIC and POC show an inverse correlation with river discharge. The higher input of soil organic matter and respired CO 2 during wet seasons was likely the main driver for the convergence of δ 13 C between DIC and DOC or POC, whereas enhanced in situ primary production and respiration during dry seasons might be responsible for their isotopic divergence. Carbon export fluxes from the Mississippi River were estimated to be 13.6 Tg C yr −1 for DIC, 1.88 Tg C yr −1 for DOC, and 2.30 Tg C yr −1 for POC during 2006–2008. The discharge-normalized DIC yield decreased during wet seasons, while those of POC and DOC increased and remained constant, respectively, implying variable responses in carbon export to the increasing discharge.

Louisiana↗

USGS2015

Reactive transport modeling of geochemical controls on secondary water quality impacts at a crude oil spill site near Bemidji, MN

Anaerobic biodegradation of organic amendments and contaminants in aquifers can trigger secondary water quality impacts that impair groundwater resources. Reactive transport models help elucidate how diverse geochemical reactions control the spatiotemporal evolution of these impacts. Using extensive monitoring data from a crude oil spill site near Bemidji, Minnesota (USA), we implemented a comprehensive model that simulates secondary plumes of depleted dissolved O 2 and elevated concentrations of Mn 2+ , Fe 2+ , CH 4 , and Ca 2+ over a two-dimensional cross section for 30 years following the spill. The model produces observed changes by representing multiple oil constituents and coupled carbonate and hydroxide chemistry. The model includes reactions with carbonates and Fe and Mn mineral phases, outgassing of CH 4 and CO 2 gas phases, and sorption of Fe, Mn, and H + . Model results demonstrate that most of the carbon loss from the oil (70%) occurs through direct outgassing from the oil source zone, greatly limiting the amount of CH 4 cycled down-gradient. The vast majority of reduced Fe is strongly attenuated on sediments, with most (91%) in the sorbed form in the model. Ferrous carbonates constitute a small fraction of the reduced Fe in simulations, but may be important for furthering the reduction of ferric oxides. The combined effect of concomitant redox reactions, sorption, and dissolved CO 2 inputs from source-zone degradation successfully reproduced observed pH. The model demonstrates that secondary water quality impacts may depend strongly on organic carbon properties, and impacts may decrease due to sorption and direct outgassing from the source zone.

Minnesota↗

USGS2015

Beyond annual streamflow reconstructions for the Upper Colorado River Basin: a paleo-water-balance approach

In this paper, we present a methodology to use annual tree-ring chronologies and a monthly water balance model to generate annual reconstructions of water balance variables (e.g., potential evapotrans- piration ( PET ), actual evapotranspiration ( AET ), snow water equivalent ( SWE ), soil moisture storage ( SMS ), and runoff ( R )). The method involves resampling monthly temperature and precipitation from the instrumental record directed by variability indicated by the paleoclimate record. The generated time series of monthly temperature and precipitation are subsequently used as inputs to a monthly water balance model. The methodology is applied to the Upper Colorado River Basin, and results indicate that the methodology reliably simulates water-year runoff, maximum snow water equivalent, and seasonal soil moisture storage for the instrumental period. As a final application, the methodology is used to produce time series of PET , AET , SWE , SMS , and R for the 1404–1905 period for the Upper Colorado River Basin.

Water Resources Research↗

USGS2015

Accelerating advances in continental domain hydrologic modeling

In the past, hydrologic modeling of surface water resources has mainly focused on simulating the hydrologic cycle at local to regional catchment modeling domains. There now exists a level of maturity among the catchment, global water security, and land surface modeling communities such that these communities are converging toward continental domain hydrologic models. This commentary, written from a catchment hydrology community perspective, provides a review of progress in each community toward this achievement, identifies common challenges the communities face, and details immediate and specific areas in which these communities can mutually benefit one another from the convergence of their research perspectives. Those include: (1) creating new incentives and infrastructure to report and share model inputs, outputs, and parameters in data services and open access, machine-independent formats for model replication or reanalysis; (2) ensuring that hydrologic models have: sufficient complexity to represent the dominant physical processes and adequate representation of anthropogenic impacts on the terrestrial water cycle, a process-based approach to model parameter estimation, and appropriate parameterizations to represent large-scale fluxes and scaling behavior; (3) maintaining a balance between model complexity and data availability as well as uncertainties; and (4) quantifying and communicating significant advancements toward these modeling goals.

Water Resources Research↗

USGS2015

Hydrologic implications of GRACE satellite data in the Colorado River Basin

Use of GRACE (Gravity Recovery and Climate Experiment) satellites for assessing global water resources is rapidly expanding. Here we advance application of GRACE satellites by reconstructing long-term total water storage (TWS) changes from ground-based monitoring and modeling data. We applied the approach to the Colorado River Basin which has experienced multiyear intense droughts at decadal intervals. Estimated TWS declined by 94 km 3 during 1986–1990 and by 102 km 3 during 1998–2004, similar to the TWS depletion recorded by GRACE (47 km 3 ) during 2010–2013. Our analysis indicates that TWS depletion is dominated by reductions in surface reservoir and soil moisture storage in the upper Colorado basin with additional reductions in groundwater storage in the lower basin. Groundwater storage changes are controlled mostly by natural responses to wet and dry cycles and irrigation pumping outside of Colorado River delivery zones based on ground-based water level and gravity data. Water storage changes are controlled primarily by variable water inputs in response to wet and dry cycles rather than increasing water use. Surface reservoir storage buffers supply variability with current reservoir storage representing ∼2.5 years of available water use. This study can be used as a template showing how to extend short-term GRACE TWS records and using all available data on storage components of TWS to interpret GRACE data, especially within the context of droughts.

Water Resources Research↗