USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrology Journal”

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 343 records · Page 19Linked to original sources
USGS1995

Nitrate transport and transformation processes in unsaturated porous media

A series of experiments was conducted on two contrasting agricultural soils to observe the influence of soil texture, preferential flow, and plants on nitrate transport and denitrification under unsaturated conditions. Calcium nitrate fertilizer was applied to the surface of four large undisturbed soil cores (30 cm diameter by 40 cm height). Two of the cores were a structured clay obtained from central Missouri and two were an unstructured fine sand obtained from central Florida. The cores were irrigated daily and maintained at a matric potential of -20 kPa, representative of soil tension in the rooting zone of irrigated agricultural fields. Volumetric water content (θ), concentration of nitrate-N in the soil solution, and nitrous oxide flux at the surface, 10, 20, and 30 cm were monitored daily. Leaching loss of surface-applied N0 3 − -N was significant in both the sand and the clay. In unplanted sand cores, almost all of the applied nitrate was leached below 30 cm within 10 days. Gaseous N loss owing to denitrification was no greater than 2% of the nitrate-N applied to the unplanted sand cores and, in general, was less than 1 %. Although leaching was somewhat retarded in the clay cores, about 60% of the applied nitrate-N was leached from the unplanted clay soil in 5–6 weeks. Under unsaturated conditions, the clay had little to no tendency to denitrify despite the greater moisture content of the clay and retarded leaching of nitrate in the clay. The planted sand cores had surprisingly large gaseous N loss owing to denitrification, as much as 17% of the nitrate-N. Results from both the clay and sand experiments show that the dynamics of nitrate transport and transformation in unsaturated soils are affected by small, localized variations in the soil moisture content profile, the gaseous diffusion coefficient of the soil, the rate at which the nitrate pulse passes through the soil, the solubility of N 2 O and N 2 and the diffusion of the gasses through the soil solution, and development of a water content profile in the soil. Limited dentrification in the clay soil was due to a limited volume of soil available for infiltration after internal catchment and the development of denitrifying conditions resulting from the presence of an extensive macropore system.

Journal of Hydrology↗
USGS1994

Peak flow rate and recession-curve characteristics of a karst spring in the Inner Bluegrass, central Kentucky

The flow rate at the terminal spring of a 1929 ha karst ground-water catchment has been continuously monitored for 2 years, and 108 identifiable events were analyzed. The peak flow rates followed a beta frequency distribution with parameters α = 0.365 and γ = 1.135. Events were separated into high-flow and low-flow. High-flow events had characteristics attributable to pipe flow. Correlation and stepwise regression were used to develop peak flow rate prediction equations for the combined 108 events and for the 81 low-flow events. The portion of the recession curve identified as pipe flow was a watershed constant and time invariant. The base flow was seasonal, increasing in the winter to approximately 0.071 m 3 s −1 and decreasing in the summer to approximately 0.014 m 3 s −1 .

Journal of Hydrology↗
USGS1994

Simulation of long-term thermal characteristics of three Estonian lakes

A one-dimensional surface energy-balance lake model, coupled to a thermodynamic model of lake ice, is used to simulate variations in the temperature of and evaporation from three Estonian lakes: Karujärv, Viljandi and Kirjaku. The model is driven by daily climate data, derived by cubic-spline interpolation from monthly mean data, and was run for periods of 8 years (Kirjaku) up to 30 years (Viljandi). Simulated surface water temperature is in good agreement with observations: mean differences between simulated and observed temperatures are from −0.8°C to +0.1°C. The simulated duration of snow and ice cover is comparable with observed. However, the model generally underpredicts ice thickness and overpredicts snow depth. Sensitivity analyses suggest that the model results are robust across a wide range (0.1–2.0 m −1) of lake extinction coefficient: surface temperature differs by less than 0.5°C between extreme values of the extinction coefficient. The model results are more sensitive to snow and ice albedos. However, changing the snow (0.2–0.9) and ice (0.15–0.55) albedos within realistic ranges does not improve the simulations of snow depth and ice thickness. The underestimation of ice thickness is correlated with the overestimation of snow cover, since a thick snow layer insulates the ice and limits ice formation. The overestimation of snow cover results from the assumption that all the simulated winter precipitation occurs as snow, a direct consequence of using daily climate data derived by interpolation from mean monthly data.

Karujärv, Kirjaku, Viljandi↗
USGS1994

Effects of land disposal of municipal sewage sludge on fate of nitrates in soil, streambed sediment, and water quality

This study was undertaken to determine the effects of sewage-sludge disposal at the Lowry sewage-sludge-disposal area, near Denver, Colorado, on ground- and surface-water quality, to determine the fate of nitrates from sludge leachate, and to determine the source areas of leachate and the potential for additional leaching from the disposal area. Sewage-sludge disposal began in 1969. Two methods were used to apply the sludge: burial and plowing. Also, the sludge was applied both in liquid and cake forms. Data in this report represent the chemical composition of soil and streambed sediment from seven soil- and four streambed-sampling sites in 1986, chemical and bacterial composition of ground water from 28 wells from 1981 to 1987, and surface-water runoff from seven water-sampling sites from 1984 to 1987. Ground water samples were obtained from alluvial and bedrock aquifers. Samples of soil, streambed sediment, ground water and surface water were obtained for onsite measurement and chemical analysis. Measurements included determination of nitrogen compounds and major cations and anions, fecal-coliform and -streptococcus bacteria, specific conductance, and pH. Thirteen wells in the alluvial aquifer in Region 3 of the study area contain water that was probably affected by sewage-sludge leachate. The plots of concentration of nitrate with time show seasonal trends and trends caused by precipitation. In addition to yearly fluctuation, there were noticeable increases in ground-water concentrations of nitrate that coincided with increased precipitation. After 3 years of annual ground-water-quality monitoring and 4 years of a quarterly sampling program, it has been determined that leachate from the sewage-sludge-disposal area caused increased nitrite plus nitrate (as nitrogen) concentration in the alluvial ground water at the site. Soil analyses from the disposal area indicate that organic nitrogen was the dominant form of nitrogen in the soil. As a result of investigations at the research site, it has been determined that a potentially large source of contamination exists in the soils of the study area owing to increased concentrations of nitrogen, sodium, calcium, magnesium, sulfate, bicarbonate, and chloride because of sewage disposal. Continued monitoring of surface and ground water for nitrogen and the other ions previously mentioned is required to assess long-term effects of municipal sludge disposal on water quality.

Journal of Hydrology↗
USGS1994

The use of slug tests to describe vertical variations in hydraulic conductivity

Multilevel slug tests provide one means of obtaining estimates of hydraulic conductivity on a scale of relevance for contaminant transport investigations. A numerical model is employed here to assess the potential of multilevel slug tests to provide information about vertical variations in hydraulic conductivity under conditions commonly faced in field settings. The results of the numerical simulations raise several important issues concerning the effectiveness of this technique. If the length of the test interval is of the order of the average layer thickness, considerable error may be introduced into the conductivity estimates owing to the effects of adjoining layers. The influence of adjoining layers is dependent on the aspect ratio (length of test interval/well radius) of the tesy interval and the flow properties of the individual layers. If a low-permeability skin is present at the well, the measured vertical variations will be much less than the actual variations, owing to the influence of the skin conductivity on the parameter estimates. A high-permeability skin can also produce apparent vertical variations that are much less than the actual, owing to water flowing vertically along the conductive skin. In cases where the test interval spans a number of layers, a slug test will yield an approximate thickness-weighted average of the hydraulic conductivities of the intersected layers. In most cases, packer circumvention should not be a major concern when packers of 0.75 m or longer are employed. Results of this study are substantiated by recently reported field tests that demonstrate the importance of well emplacement and development activities for obtaining meaningful estimates from a program of multilevel slug tests.

Journal of Hydrology↗
USGS1994

Comparison of energy-budget evaporation losses from two morphometrically different Florida seepage lakes

Evaporation was computed by the energy-budget method for two north Florida lakes with similar surface areas but different depths, for the period May 1989 to December 1990. Lake Barco, in north-central Florida, is shallow, with an average depth of 3 m; Lake Five-O, in the Florida panhandle, is considerably deeper, with an average depth of 9.5 m. As a result, the thermal regime and seasonal evaporation rates of the lakes are different. Evaporation from the shallower lake was higher than that from the deeper lake in the winter and spring. In the late summer and autumn, however, the situation is reversed. Evaporation from the shallow lake is directly related to the amount of incoming shortwave radiation because of its limited ability to store energy. The lag in evaporation at the deeper lake is a function of the greater amount of heat that it seasonally stores and releases. The difference in annual evaporation between Lake Barco (151 cm year −1 ) and Lake Five-O (128 cm year −1 ) is related to differences in regional climatic conditions between the two sites. Additionally, higher than normal evaporation rates at the two lakes are probably related to drought conditions experienced in north Florida during 1990, which resulted in higher temperatures and more incoming radiation. Monthly evaporation at Lake Barco could usually be estimated within 10% of the energy-budget evaporation using a constant pan coefficient. This lake may be representative of other shallow lakes that do not store considerable heat. Monthly evaporation at Lake Five-O, however, could not be estimated accurately by using an annual pan coefficient because of the large seasonal influence of change in stored heat. Monthly mass-transfer evaporation compared well with energy-budget evaporation at Lake Barco, but did not compare well at Lake Five-O. These errors may also be associated with changes in heat storage. Thus, the thermal regime of the lake must be considered to estimate accurately the seasonal evaporation rates from a deep lake.

Journal of Hydrology↗
USGS1994

Relationships between atmospheric circulation and snowpack in the Gunnison River basin, Colorado

In this study, winter mean 700 mbar height anomalies over the eastern North Pacific Ocean and the western USA are related to variability in snowpack accumulations measured on or about 1 April at 21 snowcourse stations within and near the Gunnison River basin in Colorado. Results indicate that lower than normal snowpack accumulations are primarily associated with positive 700 mbar height anomalies (anomalous anticyclonic circulation) over the western USA. Moist air from the Pacific Ocean is moved to the north of the western USA along the western margin of the anomalous anticyclonic circulation. In contrast, higher than normal snowpack accumulations are associated with negative 700 mbar height anomalies (anomalous cyclonic circulation) over the western USA and over most of the eastern North Pacific Ocean. The anomalous cyclonic circulation over the western USA enhances the movement of moisture from the Pacific Ocean into the southern and central parts of the West. Results also indicate that variability in winter mean 700 mbar height anomalies can explain over 50% of the variability in snowpack accumulations in the Gunnison River basin. The significant linear relationships between 700 mbar height anomalies and snowpack accumulations in the Gunnison River basin can be used in conjunction with general circulation model simulations of 700 mbar height anomalies for future climatic conditions to estimate future snowpack accumulations in the Gunnison River basin.

Journal of Hydrology↗
USGS1994

Model for dolomite formation in northwest Florida

Petrographic methods are used to examine the nature of the Floridan post-Oligocene unconformity and diagenetic responses in the associated carbonate units. The study addresses the depositional environment and the character of sediments associated with the unconformity, types and timing of diagenetic changes, and mode of dolomitization for the Suwannee Limestone dolomites and basal Hawthorn Group dolosilts. The purpose of the study is to determine if the original texture of geologic materials is preserved more frequently in rocks deposited in a supratidal environment in the uppermost Floridan carbonates, to examine possible sources of magnesium for the dolomitization that has occurred in the Floridan aquifer, and to determine the chemical environment of dolomite formation. On the basis of lithology, the tri-county study area can be divided into two regions; one in which the uppermost Floridan carbonate is calcareous or locally silicified, and one in which the aquifer system is capped with dolomite. Lateral interfingering between the two regions is apparent, and it is inferred that the boundary between the limestone and dolomite represents a geochemical boundary between a mixed marine and freshwater environment and a subaerially exposed freshwater environment. Where dolomite has replaced the uppermost Suwannee Limestone, the dolomite is overlain by a thin unit of dolosilt. The predominant source of magnesium for the dolomitization that occurred in the study area is probably seawater.

Journal of Hydrology↗
USGS1994

The travel-time ellipse: An approximate zone of transport

A zone of transport for a well is defined as the area in the horizontal plane bounded by a contour of equal ground-water travel time to the well. For short distances and ground-water travel times near a well, the potentiometric surface may be simulated analytically as that for a fully penetrating well in a uniform flow field. The zone of transport for this configuration is nearly elliptical. A simple method is derived to calculate a travel-time ellipse that approximates the zone of transport for a well in a uniform flow field. The travel-time ellipse was nearly congruent with the exact solution for the theoretical zone of transport for ground-water travel times of at least 10 years and for aquifer property values appropriate for southeastern Minnesota. For distances and travel times approaching infinity, however, the ellipse becomes slightly wider at its midpoint and narrower near its upgradient boundary than the theoretical zone of transport. The travel-time ellipse also may be used to simulate the plume area surrounding an injection well. However, the travel-time ellipse is an approximation that does not account for the effect of dispersion in enlarging the true area of an injection plume or zone of transport; hence, caution is advised in the use and interpretation of this simple construction.

Journal of Hydrology↗
USGS1994

The effect of a confining unit on the geochemical evolution of ground water in the Upper Floridan aquifer system

In west-central Florida, sections of the Upper Floridan aquifer system range in character from confined to leaky to unconfined. The confining unit is the Hawthorn Formation, a clay-rich sequence. The presence or absence of the Hawthorn Formation affects the geochemical evolution of the ground water in the Upper Floridan aquifer system. Mass-balance and mass-transfer models suggest that, in unconfined areas, the geochemical reactions are dolomite dissolution, ion exchange (Mg for Na, K), sulfate reduction, calcite dissolution, and CO 2 exchange. In the areas in which the Hawthorn Formation is leaky, the evolution of the ground water is accounted for by ion exchange, sulfate reduction, calcite dissolution, and CO 2 exchange. In the confined areas, no ion exchange and only limited sulfate reduction occur, and the chemical character of the ground water is consistent with dolomite and gypsum dissolution, calcite precipitation, and CO 2 ingassing. The Hawthorn Formation acts both as a physical barrier to the transport of CO 2 and organic matter and as a source of ion-exchange sites, but the carbonate-mineral reactions are largely unaffected by the extent of confinement of the Upper Floridan aquifer.

Journal of Hydrology↗
USGS1996

Transverse dispersion of contaminants in fractured permeable formations

Our interest in understanding some of the mechanisms involved in the mineralization of the Great Bend Prairie aquifer of Kansas by salt water originating from Permian bedrock formations, which are fractured sandstones, has prompted this basic study. The fractured permeable formatiodis represented by a simplified conceptual model incorporating two sets of oblique and parallel fractures embedded in permeable blocks. The domain is initially divided into a completely freshwater zone overlying a completely saline water zone. However, the sharp interface originally existing between the saline and fresh water is subject to dispersion because of mixing in fracture intersections and between the fracture flow and the permeable block flow. Simulations based on the use of an appropriate numerical model developed in the present study have helped us to characterize dispersion of the sharp interface and creation of the transition zone. Relationships of transverse and longitudinal dispersion in the domain are also determined.

Journal of Hydrology↗
USGS1996

Analytical solutions to non-Fickian subsurface dispersion in uniform groundwater flow

Analytical solutions are obtained by the Fourier transform technique for the one-, two-, and three-dimensional transport of a conservative solute injected instantaneously in a uniform groundwater flow. These solutions account for dispersive non-linearity caused by the heterogeneity of the hydraulic properties of aquifer systems and can be used as building blocks to construct solutions by convolution (principle of superposition) for source conditions other than slug injection. The dispersivity is assumed to vary parabolically with time and is thus constant for the entire system at any given time. Two approaches for estimating time-dependent dispersion parameters are developed for two-dimensional plumes. They both require minimal field tracer test data and, therefore, represent useful tools for assessing real-world aquifer contamination sites. The first approach requires mapped plume-area measurements at two specific times after the tracer injection. The second approach requires concentration-versus-time data from two sampling wells through which the plume passes. Detailed examples and comparisons with other procedures show that the methods presented herein are sufficiently accurate and easier to use than other available methods.

Journal of Hydrology↗
USGS2015

Sediment yields from small, steep coastal watersheds of California

Global inventories of sediment discharge to the ocean highlight the importance of small, steep watersheds (i.e., those having drainage areas less than 100,000 km2 and over 1000 m of relief) that collectively provide a dominant flux of sediment. The smallest of these coastal watersheds (e.g., those that have drainage areas less than 1000 km2) can represent a large portion of the drainage areas of active margin coasts, such as California’s coast, but remain almost universally unmonitored. Here we report on the suspended-sediment discharge of several small coastal watersheds (10-56 km2) of the Santa Ynez Mountains, California, that were found to have ephemeral discharge and suspended-sediment concentrations ranging between 1 and over 200,000 mgL-1. Sediment concentrations were weakly correlated with discharge (r2 = 0.10–0.25), and all types of hysteresis patterns were observed during high flows (clockwise, counterclockwise, no hysteresis, and complex). Sediment discharge varied strongly with time and was measurably elevated in one watershed following a wildfire. Although sediment yields varied by over 100-fold across the watersheds (e.g., 15 – 2100 tkm-2 yr -1during the relatively wet 2005 water year), the majority of sediment discharge (65-80%) occurred during only 1% of the time for all watersheds. Furthermore, sampling of dozens of high flow events provides evidence that sediment yields were generally related to peak discharge yields, although these relationships were not consistent across the watersheds. These results suggest that small watersheds of active margins can provide large fluxes of sediment to the coast, but that the rates and timing of this sediment discharge is more irregular in time – and thus more difficult to characterize – than the better monitored and studied watersheds that are 1000-100,000 km2.

California↗
USGS2019

Probability of streamflow permanence model (PROSPER): A spatially continuous model of annual streamflow permanence throughout the Pacific Northwest

The U.S. Geological Survey (USGS) has developed the PRObability of Streamflow PERmanence (PROSPER) model, a GIS raster-based empirical model that provides streamflow permanence probabilities (probabilistic predictions) of a stream channel having year-round flow for any unregulated and minimally-impaired stream channel in the Pacific Northwest region, U.S. The model provides annual predictions for 2004-2016 at a 30-m spatial resolution based on monthly or annually updated values of climatic conditions and static physiographic variables associated with the upstream basin. Predictions correspond to any pixel on the channel network consistent with the medium resolution National Hydrography Dataset channel network stream grid. Total annual precipitation and percent forest cover were consistently the most important predictor variables among global and most subregional models, which had error rates between 17 and 22%. Probabilities were converted to wet and dry streamflow permanence classes with an associated confidence. Wet and dry classifications were used to derive descriptors that characterize the statistical and spatial distribution of streamflow permanence in three focal basins. Predicted dry channel segments account for 52 to 92% of the stream network across the three focal basins; streamflow permanence decreased during climatically drier years. Predictions are publicly available through the USGS StreamStats platform. Results demonstrate the utility of the PROSPER model as a tool for identifying areas that may be resilient or sensitive to drought conditions, allowing for management efforts that target protecting critical reaches. Importantly, PROSPER’s successful predictive performance can be improved with new datasets of streamflow permanence underscoring the importance of field observations.

Idaho, Montana, Nevada, Oregon, Utah, Washington↗
USGS2020

Effects of John Martin Reservoir on water quality and quantity: Assessment by chemical, isotopic, and mass-balance methods

Water quality and quantity can be influenced by transit through and storage in reservoirs. Assessing such effects can be challenging, however, because of mixing and residence times, and inter-annual net storage and release from both the reservoir itself and surrounding porosity. Here, different methodologies were used to assess the effect of John Martin Reservoir (JMR), located on the Arkansas River, on water volumes and the problematic constituents salinity (total dissolved solids, TDS), selenium (Se), and uranium (U). Methodologies addressed short-term (16 months) and long-term (31 years) effects depending upon data availability. Evaporation was assessed by using isotopes of water to determine 12% short-term evaporation, and by pan evaporation and changes in storage to determine 11% long-term evaporation. Salinity, Se, and U mass balance were assessed by using chloride (Cl−) as an index by which to measure short-term gains or losses between inflows and outflows in the short term. Chloride gain from ungaged inflows skewed those results to overestimate retention. Continuous monitoring of discharge and specific conductance for inflows and outflows, along with discrete sampling for dissolved constituents were used to compute long-term, load-based mass balance. Mild gains of TDS (34,000 ± 15,000 Mg/yr) and U (0.1 ± 0.5 Mg/yr) in JMR were detected. Although the additions are small relative to uncertainty, they indicate little to no retention of TDS and U and likely additions from ungaged inflows. In contrast, an average of 0.6 ± 0.2 Mg/yr or 23% of gaged inflow Se was removed in JMR. The study illustrates the benefit of long-term records for assessing the influence of reservoirs for which net storage and release keep them from approaching steady-state conditions.

Colorado↗
USGS2021

Changing climate drives future streamflow declines and challenges in meeting water demand across the southwestern United States

Society and the environment in the arid southwestern United States depend on reliable water availability, yet current water use outpaces supply. Water demand is projected to grow in the future and climate change is expected to reduce supply. To adapt, water managers need robust estimates of future regional water supply to support management decisions. To address this need, we estimate future streamflow in seven water resource regions in the southwestern U.S. using a new SPAtially Referenced Regressions On Watershed attributes (SPARROW) streamflow model. We present streamflow projections corresponding to input data from seven climate models and two greenhouse gas Representative Concentration Pathways (RCP4.5 and 8.5) for three, thirty-year intervals centered on the 2030s, 2050s, and 2080s, and for a historical thirty year interval centered on the 1990s. Across water resource regions, about half of the RCP4.5 models (51%) and two thirds of the RCP8.5 models (67%) indicate decreases in streamflow in the 2080s relative to the historical period. Models project maximum decreases in streamflow of 36–80% in all water resource regions for all periods and RCPs relative to historical streamflow, and maximum streamflow decreases of up to 20–45% in the 2080s at sites along the Colorado River used for measuring compliance with interstate and international water agreements. Headwaters are projected to experience the greatest declines, with substantial downstream implications. Among these estimates, the streamflows from models forced with RCP8.5 tend to be lower than those forced with RCP4.5. Not all climate models, times, and RCPs project widespread streamflow declines. The most ubiquitous streamflow increases are projected to occur in the 2030s under RCP4.5. Later time periods and enhanced greenhouse gas forcings indicate smaller regions of streamflow increase and lower accumulated streamflows, suggesting that limiting or reducing greenhouse gas concentrations could support future water availability. Although some possible streamflow increases are promising, the modest and spatially limited increases in streamflow projected for later time periods are still unlikely to be sufficient to meet the projected water demand. These results inform the likelihood of future water agreement compliance, and support developing strategies to balance water supply and demand.

Arizona, California, Colorado Idaho, Nevada, New M↗
USGS2022

Predictions and drivers of sub-reach-scale annual streamflow permanence for the upper Missouri River basin: 1989-2018

The presence of year-round surface water in streams (i.e., streamflow permanence) is an important factor for identifying aquatic habitat availability, determining the regulatory status of streams, managing land use change, allocating water resources, and designing scientific studies. However, accurate, high resolution, and dynamic prediction of streamflow permanence that accounts for year-to-year variability at a regional extent is a major gap in modeling capability. Herein, we expand and adapt the U.S. Geological Survey (USGS) PRObability of Streamflow PERmanence (PROSPER) model from its original implementation in the Pacific Northwest (PROSPER PNW ) to the upper Missouri River basin (PROSPER UM ), a geographical region that includes mountain and prairie ecosystems of the northern United States. PROSPER UM is an empirical model used to estimate the probability that a stream channel has year-round flow in response to climatic conditions (monthly and annual) and static physiographic predictor variables of the upstream basin. The structure and approach of PROSPER UM are generally consistent with the PROSPER PNW model but include improved spatial resolution (10 m) and a longer modeling period. Average model accuracy was 81 %. Drainage area, upstream proportion as wetlands, and upstream proportion as developed land cover were the most important predictor variables. The PROSPER UM model identifies decreases in streamflow permanence during climatically drier years, although there is variability in the magnitude across basins highlighting geographically varying sensitivity to drought. Variability in the response of perennial streams to drought conditions among basins in the study area was also observed.

Montana, North Dakota, South Dakota, Wyoming↗
USGS2025

A generalized design flood estimation framework under stationary and non-stationary scenarios

With the ever-ongoing debate over the death of stationarity in flood time series, developing methods for flood frequency analysis that can address both stationarity and non-stationarity simultaneously has become increasingly important for design flood estimation. Existing non-stationary flood frequency analysis (NFFQ) methods are either limited to providing time-varying conditional flood quantile estimates or lack closed-form expression to estimate the design flood over the planning period. We propose a novel framework, MM-NFFQ, that introduces marginal moments (MM) estimation techniques to provide a closed-form expression to estimate design flood under non-stationarity. We demonstrate MM-NFFQ using the LP3 distribution for estimating conditional moments, but in principle, it can work for any 3-parameter distribution. We first show the proposed MM-NFFQ collapses to stationary flood frequency analysis analytically using synthetic data and then demonstrate the MM-NFFQ approach for two basins exhibiting non-stationarity in their flood time series. We further extend the analysis to selected 40 basins across CONUS and find that arid basins exhibit higher deviation from stationarity. Thus, the proposed MM-NFFQ framework can estimate traditional flood frequency curves for both, stationary and non-stationary flood processes, and can also be utilized to analyze the changes in conditional moments and marginal moments over different planning horizons.

Journal of Hydrology X↗