USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth Science”

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

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

At least 1,135 records · Page 63Linked to original sources

Prediction, time variance, and classification of hydraulic response to recharge in two karst aquifers

Many karst aquifers are rapidly filled and depleted and therefore are likely to be susceptible to changes in short-term climate variability. Here we explore methods that could be applied to model site-specific hydraulic responses, with the intent of simulating these responses to different climate scenarios from high-resolution climate models. We compare hydraulic responses (spring flow, groundwater level, stream base flow, and cave drip) at several sites in two karst aquifers: the Edwards aquifer (Texas, USA) and the Madison aquifer (South Dakota, USA). A lumped-parameter model simulates nonlinear soil moisture changes for estimation of recharge, and a time-variant convolution model simulates the aquifer response to this recharge. Model fit to data is 2.4% better for calibration periods than for validation periods according to the Nash–Sutcliffe coefficient of efficiency, which ranges from 0.53 to 0.94 for validation periods. We use metrics that describe the shapes of the impulse-response functions (IRFs) obtained from convolution modeling to make comparisons in the distribution of response times among sites and between aquifers. Time-variant IRFs were applied to 62% of the sites. Principal component analysis (PCA) of metrics describing the shapes of the IRFs indicates three principal components that together account for 84% of the variability in IRF shape: the first is related to IRF skewness and temporal spread and accounts for 51% of the variability; the second and third largely are related to time-variant properties and together account for 33% of the variability. Sites with IRFs that dominantly comprise exponential curves are separated geographically from those dominantly comprising lognormal curves in both aquifers as a result of spatial heterogeneity. The use of multiple IRF metrics in PCA is a novel method to characterize, compare, and classify the way in which different sites and aquifers respond to recharge. As convolution models are developed for additional aquifers, they could contribute to an IRF database and a general classification system for karst aquifers.

South Dakota, Texas↗

Combined impacts of current and future dust deposition and regional warming on Colorado River Basin snow dynamics and hydrology

The Colorado River provides water to 40 million people in seven western states and two countries and to 5.5 million irrigated acres. The river has long been overallocated. Climate models project runoff losses of 5–20% from the basin by mid-21st century due to human-induced climate change. Recent work has shown that decreased snow albedo from anthropogenic dust loading to the CO mountains shortens the duration of snow cover by several weeks relative to conditions prior to western expansion of the US in the mid-1800s, and advances peak runoff at Lees Ferry, Arizona, by an average of 3 weeks. Increases in evapotranspiration from earlier exposure of soils and germination of plants have been estimated to decrease annual runoff by more than 1.0 billion cubic meters, or ~5% of the annual average. This prior work was based on observed dust loadings during 2005–2008; however, 2009 and 2010 saw unprecedented levels of dust loading on snowpacks in the Upper Colorado River Basin (UCRB), being on the order of 5 times the 2005–2008 loading. Building on our prior work, we developed a new snow albedo decay parameterization based on observations in 2009/10 to mimic the radiative forcing of extreme dust deposition. We convolve low, moderate, and extreme dust/snow albedos with both historic climate forcing and two future climate scenarios via a delta method perturbation of historic records. Compared to moderate dust, extreme dust absorbs 2× to 4× the solar radiation, and shifts peak snowmelt an additional 3 weeks earlier to a total of 6 weeks earlier than pre-disturbance. The extreme dust scenario reduces annual flow volume an additional 1% (6% compared to pre-disturbance), a smaller difference than from low to moderate dust scenarios due to melt season shifting into a season of lower evaporative demand. The sensitivity of flow timing to dust radiative forcing of snow albedo is maintained under future climate scenarios, but the sensitivity of flow volume reductions decreases with increased climate forcing. These results have implications for water management and suggest that dust abatement efforts could be an important component of any climate adaptation strategies in the UCRB.

Arizona, Colorado, New Mexico, Utah, Wyoming↗

Large scale snow water status monitoring: Comparison of different snow water products in the upper Colorado basins

We illustrate the ability to monitor the status of snow water content over large areas by using a spatially distributed snow accumulation and ablation model that uses data from a weather forecast model in the upper Colorado Basin. The model was forced with precipitation fields from the National Weather Service (NWS) Multi-sensor Precipitation Estimator (MPE) and the Tropical Rainfall Measuring Mission (TRMM) data-sets; remaining meteorological model input data were from NOAA's Global Forecast System (GFS) model output fields. The simulated snow water equivalent (SWE) was compared to SWEs from the Snow Data Assimilation System (SNODAS) and SNOwpack TELemetry system (SNOTEL) over a region of the western US that covers parts of the upper Colorado Basin. We also compared the SWE product estimated from the special sensor microwave imager (SSM/I) and scanning multichannel microwave radiometer (SMMR) to the SNODAS and SNOTEL SWE data-sets. Agreement between the spatial distributions of the simulated SWE with MPE data was high with both SNODAS and SNOTEL. Model-simulated SWE with TRMM precipitation and SWE estimated from the passive microwave imagery were not significantly correlated spatially with either SNODAS or the SNOTEL SWE. Average basin-wide SWE simulated with the MPE and the TRMM data were highly correlated with both SNODAS ( r = 0.94 and r = 0.64; d.f. = 14 – d.f. = degrees of freedom) and SNOTEL ( r = 0.93 and r = 0.68; d.f. = 14). The SWE estimated from the passive microwave imagery was significantly correlated with the SNODAS SWE ( r = 0.55, d.f. = 9, p = 0.05) but was not significantly correlated with the SNOTEL-reported SWE values ( r = 0.45, d.f. = 9, p = 0.05).The results indicate the applicability of the snow energy balance model for monitoring snow water content at regional scales when coupled with meteorological data of acceptable quality. The two snow water contents from the microwave imagery (SMMR and SSM/I) and the Utah Energy Balance forced with the TRMM precipitation data were found to be unreliable sources for mapping SWE in the study area; both data sets lacked discernible variability of snow water content between sites as seen in the SNOTEL and SNODAS SWE data. This study will contribute to better understanding the adequacy of data from weather forecast models, TRMM, and microwave imagery for monitoring status of the snow water content.

Colorado, Utah, Wyoming↗

Hillslope-scale experiment demonstrates role of convergence during two-step saturation

Subsurface flow and storage dynamics at hillslope scale are difficult to ascertain, often in part due to a lack of sufficient high-resolution measurements and an incomplete understanding of boundary conditions, soil properties, and other environmental aspects. A continuous and extreme rainfall experiment on an artificial hillslope at Biosphere 2's Landscape Evolution Observatory (LEO) resulted in saturation excess overland flow and gully erosion in the convergent hillslope area. An array of 496 soil moisture sensors revealed a two-step saturation process. First, the downward movement of the wetting front brought soils to a relatively constant but still unsaturated moisture content. Second, soils were brought to saturated conditions from below in response to rising water tables. Convergent areas responded faster than upslope areas, due to contributions from lateral subsurface flow driven by the topography of the bottom boundary, which is comparable to impermeable bedrock in natural environments. This led to the formation of a groundwater ridge in the convergent area, triggering saturation excess runoff generation. This unique experiment demonstrates, at very high spatial and temporal resolution, the role of convergence on subsurface storage and flow dynamics. The results bring into question the representation of saturation excess overland flow in conceptual rainfall-runoff models and land-surface models, since flow is gravity-driven in many of these models and upper layers cannot become saturated from below. The results also provide a baseline to study the role of the co-evolution of ecological and hydrological processes in determining landscape water dynamics during future experiments in LEO.

Hydrology and Earth System Sciences↗

Hydroclimatic regimes: a distributed water-balance framework for hydrologic assessment, classification, and management

Runoff-based indicators of terrestrial water availability are appropriate for humid regions, but have tended to limit our basic hydrologic understanding of drylands – the dry-subhumid, semiarid, and arid regions which presently cover nearly half of the global land surface. In response, we introduce an indicator framework that gives equal weight to humid and dryland regions, accounting fully for both vertical (precipitation + evapotranspiration) and horizontal (groundwater + surface-water) components of the hydrologic cycle in any given location – as well as fluxes into and out of landscape storage. We apply the framework to a diverse hydroclimatic region (the conterminous USA) using a distributed water-balance model consisting of 53 400 networked landscape hydrologic units. Our model simulations indicate that about 21% of the conterminous USA either generated no runoff or consumed runoff from upgradient sources on a mean-annual basis during the 20th century. Vertical fluxes exceeded horizontal fluxes across 76% of the conterminous area. Long-term-average total water availability (TWA) during the 20th century, defined here as the total influx to a landscape hydrologic unit from precipitation, groundwater, and surface water, varied spatially by about 400 000-fold, a range of variation ~100 times larger than that for mean-annual runoff across the same area. The framework includes but is not limited to classical, runoff-based approaches to water-resource assessment. It also incorporates and reinterprets the green- and blue-water perspective now gaining international acceptance. Implications of the new framework for several areas of contemporary hydrology are explored, and the data requirements of the approach are discussed in relation to the increasing availability of gridded global climate, land-surface, and hydrologic data sets.

Hydrology and Earth System Sciences↗

Higher-order statistical moments and a procedure that detects potentially anomalous years as two alternative methods describing alterations in continuous environmental data

Statistics of central tendency and dispersion may not capture relevant or desired characteristics of the distribution of continuous phenomena and, thus, they may not adequately describe temporal patterns of change. Here, we present two methodological approaches that can help to identify temporal changes in environmental regimes. First, we use higher-order statistical moments (skewness and kurtosis) to examine potential changes of empirical distributions at decadal extents. Second, we adapt a statistical procedure combining a non-metric multidimensional scaling technique and higher density region plots to detect potentially anomalous years. We illustrate the use of these approaches by examining long-term stream temperature data from minimally and highly human-influenced streams. In particular, we contrast predictions about thermal regime responses to changing climates and human-related water uses. Using these methods, we effectively diagnose years with unusual thermal variability and patterns in variability through time, as well as spatial variability linked to regional and local factors that influence stream temperature. Our findings highlight the complexity of responses of thermal regimes of streams and reveal their differential vulnerability to climate warming and human-related water uses. The two approaches presented here can be applied with a variety of other continuous phenomena to address historical changes, extreme events, and their associated ecological responses.

Hydrology and Earth System Sciences↗

A 2-D process-based model for suspended sediment dynamics: A first step towards ecological modeling

In estuaries suspended sediment concentration (SSC) is one of the most important contributors to turbidity, which influences habitat conditions and ecological functions of the system. Sediment dynamics differs depending on sediment supply and hydrodynamic forcing conditions that vary over space and over time. A robust sediment transport model is a first step in developing a chain of models enabling simulations of contaminants, phytoplankton and habitat conditions. This works aims to determine turbidity levels in the complex-geometry delta of the San Francisco estuary using a process-based approach (Delft3D Flexible Mesh software). Our approach includes a detailed calibration against measured SSC levels, a sensitivity analysis on model parameters and the determination of a yearly sediment budget as well as an assessment of model results in terms of turbidity levels for a single year, water year (WY) 2011. Model results show that our process-based approach is a valuable tool in assessing sediment dynamics and their related ecological parameters over a range of spatial and temporal scales. The model may act as the base model for a chain of ecological models assessing the impact of climate change and management scenarios. Here we present a modeling approach that, with limited data, produces reliable predictions and can be useful for estuaries without a large amount of processes data.

Hydrology and Earth System Sciences↗

Long-term trends in a Dimictic Lake

The one-dimensional hydrodynamic ice model, DYRESM-WQ-I, was modified to simulate ice cover and thermal structure of dimictic Lake Mendota, Wisconsin, USA, over a continuous 104-year period (1911–2014). The model results were then used to examine the drivers of changes in ice cover and water temperature, focusing on the responses to shifts in air temperature, wind speed, and water clarity at multiyear timescales. Observations of the drivers include a change in the trend of warming air temperatures from 0.081 °C per decade before 1981 to 0.334 °C per decade thereafter, as well as a shift in mean wind speed from 4.44 m s −1 before 1994 to 3.74 m s −1 thereafter. Observations show that Lake Mendota has experienced significant changes in ice cover: later ice-on date(9.0 days later per century), earlier ice-off date (12.3 days per century), decreasing ice cover duration (21.3 days per century), while model simulations indicate a change in maximum ice thickness (12.7 cm decrease per century). Model simulations also show changes in the lake thermal regime of earlier stratification onset (12.3 days per century), later fall turnover (14.6 days per century), longer stratification duration (26.8 days per century), and decreasing summer hypolimnetic temperatures (−1.4 °C per century). Correlation analysis of lake variables and driving variables revealed ice cover variables, stratification onset, epilimnetic temperature, and hypolimnetic temperature were most closely correlated with air temperature, whereas freeze-over water temperature, hypolimnetic heating, and fall turnover date were more closely correlated with wind speed. Each lake variable (i.e., ice-on and ice-off dates, ice cover duration, maximum ice thickness, freeze-over water temperature, stratification onset, fall turnover date, stratification duration, epilimnion temperature, hypolimnion temperature, and hypolimnetic heating) was averaged for the three periods (1911–1980, 1981–1993, and 1994–2014) delineated by abrupt changes in air temperature and wind speed. Average summer hypolimnetic temperature and fall turnover date exhibit significant differences between the third period and the first two periods. Changes in ice cover (ice-on and ice-off dates, ice cover duration, and maximum ice thickness) exhibit an abrupt change after 1994, which was related in part to the warm El Niño winter of 1997–1998. Under-ice water temperature, freeze-over water temperature, hypolimnetic temperature, fall turnover date, and stratification duration demonstrate a significant difference in the third period (1994–2014), when air temperature was warmest and wind speeds decreased rather abruptly. The trends in ice cover and water temperature demonstrate responses to both long-term and abrupt changes in meteorological conditions that can be complemented with numerical modeling to better understand how these variables will respond in a future climate.

Wisconsin↗

Ordinary kriging as a tool to estimate historical daily streamflow records

Efficient and responsible management of water resources relies on accurate streamflow records. However, many watersheds are ungaged, limiting the ability to assess and understand local hydrology. Several tools have been developed to alleviate this data scarcity, but few provide continuous daily streamflow records at individual streamgages within an entire region. Building on the history of hydrologic mapping, ordinary kriging was extended to predict daily streamflow time series on a regional basis. Pooling parameters to estimate a single, time-invariant characterization of spatial semivariance structure is shown to produce accurate reproduction of streamflow. This approach is contrasted with a time-varying series of variograms, representing the temporal evolution and behavior of the spatial semivariance structure. Furthermore, the ordinary kriging approach is shown to produce more accurate time series than more common, single-index hydrologic transfers. A comparison between topological kriging and ordinary kriging is less definitive, showing the ordinary kriging approach to be significantly inferior in terms of Nash–Sutcliffe model efficiencies while maintaining significantly superior performance measured by root mean squared errors. Given the similarity of performance and the computational efficiency of ordinary kriging, it is concluded that ordinary kriging is useful for first-order approximation of daily streamflow time series in ungaged watersheds.

Hydrology and Earth System Sciences↗

Parameter regionalization of a monthly water balance model for the conterminous United States

A parameter regionalization scheme to transfer parameter values from gaged to ungaged areas for a monthly water balance model (MWBM) was developed and tested for the conterminous United States (CONUS). The Fourier Amplitude Sensitivity Test, a global-sensitivity algorithm, was implemented on a MWBM to generate parameter sensitivities on a set of 109 951 hydrologic response units (HRUs) across the CONUS. The HRUs were grouped into 110 calibration regions based on similar parameter sensitivities. Subsequently, measured runoff from 1575 streamgages within the calibration regions were used to calibrate the MWBM parameters to produce parameter sets for each calibration region. Measured and simulated runoff at the 1575 streamgages showed good correspondence for the majority of the CONUS, with a median computed Nash–Sutcliffe efficiency coefficient of 0.76 over all streamgages. These methods maximize the use of available runoff information, resulting in a calibrated CONUS-wide application of the MWBM suitable for providing estimates of water availability at the HRU resolution for both gaged and ungaged areas of the CONUS.

Hydrology and Earth System Sciences↗

An index of floodplain surface complexity

Floodplain surface topography is an important component of floodplain ecosystems. It is the primary physical template upon which ecosystem processes are acted out, and complexity in this template can contribute to the high biodiversity and productivity of floodplain ecosystems. There has been a limited appreciation of floodplain surface complexity because of the traditional focus on temporal variability in floodplains as well as limitations to quantifying spatial complexity. An index of floodplain surface complexity (FSC) is developed in this paper and applied to eight floodplains from different geographic settings. The index is based on two key indicators of complexity, variability in surface geometry (VSG) and the spatial organisation of surface conditions (SPO), and was determined at three sampling scales. FSC, VSG, and SPO varied between the eight floodplains and these differences depended upon sampling scale. Relationships between these measures of spatial complexity and seven geomorphological and hydrological drivers were investigated. There was a significant decline in all complexity measures with increasing floodplain width, which was explained by either a power, logarithmic, or exponential function. There was an initial rapid decline in surface complexity as floodplain width increased from 1.5 to 5 km, followed by little change in floodplains wider than 10 km. VSG also increased significantly with increasing sediment yield. No significant relationships were determined between any of the four hydrological variables and floodplain surface complexity.

Hydrology and Earth System Sciences↗

Influence of groundwater on distribution of dwarf wedgemussels (Alasmidonta heterodon) in the upper reaches of the Delaware River, northeastern USA

The remaining populations of the endangered dwarf wedgemussel (DWM) ( Alasmidonta heterodon ) in the upper Delaware River, northeastern USA, were hypothesized to be located in areas of greater-than-normal groundwater discharge to the river. We combined physical (seepage meters, monitoring wells and piezometers), thermal (fiber-optic distributed temperature sensing, infrared, vertical bed-temperature profiling), and geophysical (electromagnetic-induction) methods at several spatial scales to characterize known DWM habitat and explore this hypothesis. Numerous springs were observed using visible and infrared imaging along the river banks at all three known DWM-populated areas, but not in adjacent areas where DWM were absent. Vertical and lateral groundwater gradients were toward the river along all three DWM-populated reaches, with median upward gradients 3 to 9 times larger than in adjacent reaches. Point-scale seepage-meter measurements indicated that upward seepage across the riverbed was faster and more consistently upward at DWM-populated areas. Discrete and areally distributed riverbed-temperature measurements indicated numerous cold areas of groundwater discharge during warm summer months; all were within areas populated by DWM. Electromagnetic-induction measurements, which may indicate riverbed geology, showed patterning but little correlation between bulk streambed electromagnetic conductivity and areal distribution of DWM. In spite of complexity introduced by hyporheic exchange, multiple lines of research provide strong evidence that DWM are located within or directly downstream of areas of substantial focused groundwater discharge to the river. Broad scale thermal-reconnaissance methods (e.g., infrared) may be useful in locating and protecting other currently unknown mussel populations.

Delaware River↗

Towards simplification of hydrologic modeling: Identification of dominant processes

The Precipitation–Runoff Modeling System (PRMS), a distributed-parameter hydrologic model, has been applied to the conterminous US (CONUS). Parameter sensitivity analysis was used to identify: (1) the sensitive input parameters and (2) particular model output variables that could be associated with the dominant hydrologic process(es). Sensitivity values of 35 PRMS calibration parameters were computed using the Fourier amplitude sensitivity test procedure on 110 000 independent hydrologically based spatial modeling units covering the CONUS and then summarized to process (snowmelt, surface runoff, infiltration, soil moisture, evapotranspiration, interflow, baseflow, and runoff) and model performance statistic (mean, coefficient of variation, and autoregressive lag 1). Identified parameters and processes provide insight into model performance at the location of each unit and allow the modeler to identify the most dominant process on the basis of which processes are associated with the most sensitive parameters. The results of this study indicate that: (1) the choice of performance statistic and output variables has a strong influence on parameter sensitivity, (2) the apparent model complexity to the modeler can be reduced by focusing on those processes that are associated with sensitive parameters and disregarding those that are not, (3) different processes require different numbers of parameters for simulation, and (4) some sensitive parameters influence only one hydrologic process, while others may influence many

Hydrology and Earth System Sciences↗

On the probability distribution of daily streamflow in the United States

Daily streamflows are often represented by flow duration curves (FDCs), which illustrate the frequency with which flows are equaled or exceeded. FDCs have had broad applications across both operational and research hydrology for decades; however, modeling FDCs has proven elusive. Daily streamflow is a complex time series with flow values ranging over many orders of magnitude. The identification of a probability distribution that can approximate daily streamflow would improve understanding of the behavior of daily flows and the ability to estimate FDCs at ungaged river locations. Comparisons of modeled and empirical FDCs at nearly 400 unregulated, perennial streams illustrate that the four-parameter kappa distribution provides a very good representation of daily streamflow across the majority of physiographic regions in the conterminous United States (US). Further, for some regions of the US, the three-parameter generalized Pareto and lognormal distributions also provide a good approximation to FDCs. Similar results are found for the period of record FDCs, representing the long-term hydrologic regime at a site, and median annual FDCs, representing the behavior of flows in a typical year.

Hydrology and Earth System Sciences↗

The importance of parameterization when simulating the hydrologic response of vegetative land-cover change

Computer models of hydrologic systems are frequently used to investigate the hydrologic response of land-cover change. If the modeling results are used to inform resource-management decisions, then providing robust estimates of uncertainty in the simulated response is an important consideration. Here we examine the importance of parameterization, a necessarily subjective process, on uncertainty estimates of the simulated hydrologic response of land-cover change. Specifically, we applied the soil water assessment tool (SWAT) model to a 1.4 km 2 watershed in southern Texas to investigate the simulated hydrologic response of brush management (the mechanical removal of woody plants), a discrete land-cover change. The watershed was instrumented before and after brush-management activities were undertaken, and estimates of precipitation, streamflow, and evapotranspiration (ET) are available; these data were used to condition and verify the model. The role of parameterization in brush-management simulation was evaluated by constructing two models, one with 12 adjustable parameters (reduced parameterization) and one with 1305 adjustable parameters (full parameterization). Both models were subjected to global sensitivity analysis as well as Monte Carlo and generalized likelihood uncertainty estimation (GLUE) conditioning to identify important model inputs and to estimate uncertainty in several quantities of interest related to brush management. Many realizations from both parameterizations were identified as behavioral in that they reproduce daily mean streamflow acceptably well according to Nash–Sutcliffe model efficiency coefficient, percent bias, and coefficient of determination. However, the total volumetric ET difference resulting from simulated brush management remains highly uncertain after conditioning to daily mean streamflow, indicating that streamflow data alone are not sufficient to inform the model inputs that influence the simulated outcomes of brush management the most. Additionally, the reduced-parameterization model grossly underestimates uncertainty in the total volumetric ET difference compared to the full-parameterization model; total volumetric ET difference is a primary metric for evaluating the outcomes of brush management. The failure of the reduced-parameterization model to provide robust uncertainty estimates demonstrates the importance of parameterization when attempting to quantify uncertainty in land-cover change simulations.

Hydrology and Earth System Sciences↗

Streamflow characteristics from modelled runoff time series: Importance of calibration criteria selection

Ecologically relevant streamflow characteristics (SFCs) of ungauged catchments are often estimated from simulated runoff of hydrologic models that were originally calibrated on gauged catchments. However, SFC estimates of the gauged donor catchments and subsequently the ungauged catchments can be substantially uncertain when models are calibrated using traditional approaches based on optimization of statistical performance metrics (e.g., Nash–Sutcliffe model efficiency). An improved calibration strategy for gauged catchments is therefore crucial to help reduce the uncertainties of estimated SFCs for ungauged catchments. The aim of this study was to improve SFC estimates from modeled runoff time series in gauged catchments by explicitly including one or several SFCs in the calibration process. Different types of objective functions were defined consisting of the Nash–Sutcliffe model efficiency, single SFCs, or combinations thereof. We calibrated a bucket-type runoff model (HBV – Hydrologiska Byråns Vattenavdelning – model) for 25 catchments in the Tennessee River basin and evaluated the proposed calibration approach on 13 ecologically relevant SFCs representing major flow regime components and different flow conditions. While the model generally tended to underestimate the tested SFCs related to mean and high-flow conditions, SFCs related to low flow were generally overestimated. The highest estimation accuracies were achieved by a SFC-specific model calibration. Estimates of SFCs not included in the calibration process were of similar quality when comparing a multi-SFC calibration approach to a traditional model efficiency calibration. For practical applications, this implies that SFCs should preferably be estimated from targeted runoff model calibration, and modeled estimates need to be carefully interpreted.

Tennessee↗

Characterization and evaluation of controls on post-fire streamflow response across western U.S. watersheds

This research investigates the impact of wildfires on watershed flow regimes, specifically focusing on evaluation of fire events within specified hydroclimatic regions in the western United States, and evaluating the impact of climate and geophysical variables on response. Eighty-two watersheds were identified with at least 10 years of continuous pre-fire daily streamflow records and 5 years of continuous post-fire daily flow records. Percent change in annual runoff ratio, low flows, high flows, peak flows, number of zero flow days, baseflow index, and Richards–Baker flashiness index were calculated for each watershed using pre- and post-fire periods. Independent variables were identified for each watershed and fire event, including topographic, vegetation, climate, burn severity, percent area burned, and soils data. Results show that low flows, high flows, and peak flows increase in the first 2 years following a wildfire and decrease over time. Relative response was used to scale response variables with the respective percent area of watershed burned in order to compare regional differences in watershed response. To account for variability in precipitation events, runoff ratio was used to compare runoff directly to PRISM precipitation estimates. To account for regional differences in climate patterns, watersheds were divided into nine regions, or clusters, through k-means clustering using climate data, and regression models were produced for watersheds grouped by total area burned. Watersheds in Cluster 9 (eastern California, western Nevada, Oregon) demonstrate a small negative response to observed flow regimes after fire. Cluster 8 watersheds (coastal California) display the greatest flow responses, typically within the first year following wildfire. Most other watersheds show a positive mean relative response. In addition, simple regression models show low correlation between percent watershed burned and streamflow response, implying that other watershed factors strongly influence response. Spearman correlation identified NDVI, aridity index, percent of a watershed's precipitation that falls as rain, and slope as being positively correlated with post-fire streamflow response. This metric also suggested a negative correlation between response and the soil erodibility factor, watershed area, and percent low burn severity. Regression models identified only moderate burn severity and watershed area as being consistently positively/negatively correlated, respectively, with response. The random forest model identified only slope and percent area burned as significant watershed parameters controlling response. Results will help inform post-fire runoff management decisions by helping to identify expected changes to flow regimes, as well as facilitate parameterization for model application in burned watersheds.

Hydrology and Earth System Sciences↗

Technical note: False low turbidity readings from optical probes during high suspended-sediment concentrations

Turbidity, a measure of water clarity, is monitored for a variety of purposes including (1) to help determine whether water is safe to drink, (2) to establish background conditions of lakes and rivers and detect pollution caused by construction projects and stormwater discharge, (3) to study sediment transport in rivers and erosion in catchments, (4) to manage siltation of water reservoirs, and (5) to establish connections with aquatic biological properties, such as primary production and predator–prey interactions. Turbidity is typically measured with an optical probe that detects light scattered from particles in the water. Probes have defined upper limits of the range of turbidity that they can measure. The general assumption is that when turbidity exceeds this upper limit, the values of turbidity will be constant, i.e., the probe is pegged ; however, this assumption is not necessarily valid. In rivers with limited variation in the physical properties of the suspended sediment, at lower suspended-sediment concentrations, an increase in suspended-sediment concentration will cause a linear increase in turbidity. When the suspended-sediment concentration in these rivers is high, turbidity levels can exceed the upper measurement limit of an optical probe and record a constant pegged value. However, at extremely high suspended-sediment concentrations, optical turbidity probes do not necessarily stay pegged at a constant value. Data from the Colorado River in Grand Canyon, Arizona, USA, and a laboratory experiment both demonstrate that when turbidity exceeds instrument-pegged conditions, increasing suspended-sediment concentration (and thus increasing turbidity) may cause optical probes to record decreasing false turbidity values that appear to be within the valid measurement range of the probe. Therefore, under high-turbidity conditions, other surrogate measurements of turbidity (e.g., acoustic-attenuation measurements or suspended-sediment samples) are necessary to correct these low false turbidity measurements and accurately measure turbidity.

Hydrology and Earth System Sciences↗