Search USGSSearch

SEARCH · Search USGS

Results for “Computer contribution”

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 145 records · Page 8Linked to original sources

Evapotranspiration rate measurements of vegetation typical of ground-water discharge areas in the Basin and Range carbonate-rock aquifer system, Nevada and Utah, September 2005-August 2006

Evapotranspiration was measured at six eddy-correlation sites for a 1-year period between September 1, 2005, and August 31, 2006. Five sites were in phreatophytic shrubland dominated by greasewood, and one site was in a grassland meadow. The measured annual evapotranspiration ranged from 10.02 to 12.77 inches at the shrubland sites and 26.94 inches at the grassland site. Evapotranspiration rates correlated to measured vegetation densities and to satellite-derived vegetation indexes. Evapotranspiration rates were greater at sites with denser vegetation. The primary water source supporting evapotranspiration was water derived from local precipitation at the shrubland sites, and ground water at the grassland site. Measured precipitation, ranging from 6.21 to 11.41 inches, was within 20 percent of the computed long-term annual mean. The amount of ground water consumed by phreatophytes depends primarily on local precipitation and vegetation density. The ground-water contribution to local evapotranspiration ranged from 6 to 38 percent of total evapotranspiration at the shrubland sites, and 70 percent of total evapotranspiration at the grassland site. Average depth to water ranged from 7.2 to 32.4 feet below land surface at the shrubland sites, and 3.9 feet at the grassland site. Water levels declined throughout the growing season and recovered during the non-growing season. Diurnal water-level fluctuations associated with evapotranspiration were evident at some sites but not at others.

Nevada, Utah

Sediment concentrations and loads upstream from and through John Redmond Reservoir, east-central Kansas, 2010–19

Streambank erosion and reservoir sedimentation are primary concerns of resource managers in Kansas and throughout many regions of the United States and negatively affect flood control, water supply, and recreation. The Cottonwood and upper Neosho Rivers drain into John Redmond Reservoir, and since reservoir completion in 1964, there has been substantial conservation-pool sedimentation and storage loss in John Redmond Reservoir, causing storage capacity losses more rapidly than most other Federal reservoirs in Kansas. The U.S. Geological Survey (USGS), in cooperation with the Kansas Water Office, has monitored water quality (temperature, specific conductance, and turbidity) on the Cottonwood River (upstream from the reservoir) and Neosho River (upstream and downstream from the reservoir) since 2007 with additional sites added in 2009. The purpose of this report is to quantify suspended-sediment concentrations, loads, and yields entering and exiting John Redmond Reservoir during January 1, 2010, through December 31, 2019. Three water-quality monitoring sites were upstream from the reservoir (Cottonwood River near Plymouth, Kansas [USGS site 07182250; hereinafter referred to as “Cottonwood”]; Neosho River at Burlingame Road near Emporia, Kans. [USGS site 07179750; hereinafter referred to as “Burlingame”]; and Neosho River at Neosho Rapids, Kans. [USGS site 07182390; hereinafter referred to as “Neosho Rapids”]), and one water-quality monitoring site was downstream from the reservoir (Neosho River at Burlington, Kans. [USGS site 07182510; hereinafter referred to as “Burlington”]). The Neosho Rapids streamgage is downstream from the confluence of the Cottonwood and upper Neosho Rivers and has a contributing drainage area accounting for 91 percent of the total contributing drainage area to John Redmond Reservoir. Continuously measured streamflow, water quality, and discrete water-quality data were used to develop updated regression models to compute suspended-sediment concentrations, loads, and yields upstream and downstream from John Redmond Reservoir in east-central Kansas. Several turbidity sensors were deployed during the analysis period, and there are no established relations between the sensors; therefore, individual models for each sensor were developed. Model statistics for the turbidity and suspended-sediment concentration linear regression models were better (based on the coefficient of determination, root mean square error, and model standard percentage error) than the streamflow and suspended-sediment concentration linear regression models, indicating better model performance. Computed concentrations, loads, and yields do not account for the ungaged 9 percent of the drainage basin downstream from the Neosho Rapids streamgage. Mean daily suspended-sediment loads upstream from the reservoir were largest at Neosho Rapids (2,250 tons), second largest at Cottonwood (2,180 tons), and smallest at Burlingame (624 tons). Streamflow at Burlington was predominately regulated by reservoir releases, and mean daily suspended-sediment loads were smaller (286 tons) than at upstream sites. Among the upstream sites, Cottonwood had the largest mean daily suspended-sediment concentration (179 milligrams per liter [mg/L]), followed by Neosho Rapids (162 mg/L), and Burlingame (108 mg/L). Burlington had the smallest mean daily suspended-sediment concentration of all sites (46 mg/L). Annual reservoir trapping efficiency ranged from 82 to 94 percent, and the largest sediment mass trapped was during 2019 (2,230,000 tons). Reservoir storage decreased an estimated 7,750 acre-feet during 2010 and 2014–19. Using the mean trapping efficiency to estimate suspended-sediment loads during years with missing data (2011–13), the total estimated reservoir storage lost to sedimentation for the analysis period (2010–19) was 8,690 acre-feet, about 17 percent of the remaining storage space reported in 2007. The mean annual sedimentation rate during the analysis period (747 acre-feet per year) was about 85 percent larger than the design sedimentation rate (404 acre-feet per year) originally projected during construction. Different reservoir outflow management strategies, including operating near normal capacity as opposed to higher flood pool levels, could reduce the total reservoir storage lost by 3 percent (about 261 acre-feet), which is equal to 14 percent of the total sediment removed during the dredging operation in 2016. During the study period, about 56 percent of the total suspended-sediment load was transported during streamflows greater than the National Weather Service flood action stage at the upstream sites (0.1–5 percent of the record; Cottonwood mean: 48 percent; Burlingame mean: 40 percent; Neosho Rapids mean: 78 percent). Disproportionately large sediment loads were delivered during short periods of time, and localized efforts of stream erosion protection (streambank stabilization, riparian buffers) were likely to be overwhelmed. Precipitation frequency and intensity are projected to continue to increase in this region; therefore, future sediment reduction strategies that account for extreme episodic events may be beneficial. Changes to reservoir outflow management could also minimize sediment accumulation while still preserving flood control. Continued investigation of sediment reduction measures is necessary for future mitigation with the understanding that sedimentation rate is largely driven by high flows. Results from this study can be used to calibrate sediment models, explore sediment reduction strategies, highlight the importance of continued water-quality monitoring to determine effectiveness and changes in sediment transport, and assess the ability of John Redmond Reservoir to support designated uses into the future.

Kansas

A holistic modelling approach to project the evolution of inlet-interrupted coastlines over the 21st century

Approximately one quarter of the World’s sandy beaches, most of which are interrupted by tidal inlets, are eroding. Understanding the long-term (50-100 year) evolution of inlet-interrupted coasts in a changing climate is therefore of great importance for coastal zone planners and managers. This study therefore focuses on the development and piloting of an innovative model that can simulate the climate-change driven evolution of inlet-interrupted coasts at 50-100 year time scales, while taking into account the contributions from catchment-estuary-coastal systems in a holistic manner. In this new model, the evolution of inlet-interrupted coasts is determined by: (1) computing the variation of total sediment volume exchange between the inlet-estuary system and its adjacent coast, and (2) distributing the computed sediment volume along the inlet-interrupted coast as a spatially and temporally varying quantity. The exchange volume, as computed here, consists of three major components: variation in fluvial sediment supply; basin (or estuarine) infilling due to the sea-level rise-induced increase in accommodation space; and estuarine sediment volume change due to variations in river discharge. To pilot the model, it is here applied to three different catchment-estuary-coastal systems: the Alsea estuary (Oregon, USA), Dyfi estuary (Wales, UK), and Kalutara inlet (Sri Lanka). Results indicate that all three systems will experience sediment deficits by 2100 (i.e. sediment importing estuaries). However, processes and system characteristics governing the total sediment exchange volume, and thus coastline change, vary markedly among the systems due to differences in geomorphic settings and projected climatic conditions. These results underline the importance of accounting for the different governing processes when assessing the future evolution of inlet-interrupted coastlines.

Frontiers in Marine Science

Nitrogen and phosphorus in streams of the Great Miami River Basin, Ohio, 1998–2000

Sources and loads of nitrogen and phosphorus in streams of the Great Miami River Basin were evaluated as part of the National Water-Quality Assessment program. Water samples were collected by the U.S. Geological Survey from October 1998 through September 2000 (water years 1999 and 2000) at five locations in Ohio on a routine schedule and additionally during selected high streamflows. Stillwater River near Union, Great Miami River near Vandalia, and Mad River near Eagle City were selected to represent predominantly agricultural areas upstream from the Dayton metropolitan area. Holes Creek near Kettering is in the Dayton metropolitan area and was selected to represent an urban area in the Great Miami River Basin. Great Miami River at Hamilton is downstream from the Dayton and Hamilton-Middletown metropolitan areas and was selected to represent mixed agricultural and urban land uses of the Great Miami River Basin. Inputs of nitrogen and phosphorus to streams from point and nonpoint sources were estimated for the three agricultural basins and for the Great Miami River Basin as a whole. Nutrient inputs from point sources were computed from the facilities that discharge one-half million gallons or more per day into streams of the Great Miami River Basin. Nonpoint-source inputs estimated in this report are atmospheric deposition and commercial-fertilizer and manure applications. Loads of ammonia, nitrate, total nitrogen, orthophosphate, and total phosphorus from the five sites were computed with the ESTIMATOR program. The computations show nitrate to be the primary component of instream nitrogen loads, and particulate phosphorus to be the primary component of instream phosphorus loads. The Mad River contributed the smallest loads of total nitrogen and total phosphorus to the study area upstream from Dayton, whereas the Upper Great Miami River (upstream from Vandalia) contributed the largest loads of total nitrogen and total phosphorus to the Great Miami River Basin upstream from Dayton. An evaluation of monthly mean loads shows that nutrient loads were highest during winter 1999 and lowest during the drought of summer and autumn 1999. During the 1999 drought, point sources were the primary contributors of nitrogen and phosphorus loads to most of the study area. Nonpoint sources, however, were the primary contributors of nitrogen and phosphorus loads during months of high streamflow. Nonpoint sources were also the primary contributors of nitrogen loads to the Mad River during the 1999 drought, owing to unusually large amounts of ground-water discharge to the stream. The Stillwater River Basin had the highest nutrient yields in the study area during months of high streamflow; however, the Mad River Basin had the highest yields of all nutrients except ammonia during the months of the 1999 drought. The high wet-weather yields in the Stillwater River Basin were caused by agricultural runoff, whereas high yields in the Mad River Basin during drought resulted from the large, sustained contribution of ground water to streamflow throughout the year. In the basins upstream from Dayton, an estimated 19 to 25 percent of the nonpoint source of nitrogen and 4 to 5 percent of the nonpoint source of phosphorus that was deposited or applied to the land was transported into streams.

Ohio

Annual estimates of water and solute export from 42 tributaries to the Yukon River

Annual export of 11 major and trace solutes for the Yukon River is found to be accurately determined based on summing 42 tributary contributions. These findings provide the first published estimates of tributary specific distribution of solutes within the Yukon River basin. First, we show that annual discharge of the Yukon River can be computed by summing calculated annual discharges from 42 tributaries. Annual discharge for the tributaries is calculated from the basin area and average annual precipitation over that area using a previously published regional regression equation. Based on tributary inputs, we estimate an average annual discharge for the Yukon River of 210 km 3  year –1 . This value is within 1% of the average measured annual discharge at the U.S. Geological Survey gaging station near the river terminus at Pilot Station, AK, for water years 2001 through 2005. Next, annual loads for 11 solutes are determined by combining annual discharge with point measurements of solute concentrations in tributary river water. Based on the sum of solutes in tributary water, we find that the Yukon River discharges approximately 33 million metric tons of dissolved solids each year at Pilot Station. Discharged solutes are dominated by cations calcium and magnesium (5.65 × 10 9 and 1.42 × 10 9  g year –1 ) and anions bicarbonate and sulphate (17.3 × 10 9 and 5.40 × 10 9  g year –1 ). These loads compare well with loads calculated independently at the three continuous gaging stations along the Yukon River. These findings show how annual solute yields vary throughout a major subarctic river basin and that accurate estimates of total river export can be determined from calculated tributary contributions.

Hydrological Processes

Simulations of ground-water flow and particle pathline analysis in the zone of contribution of a public-supply well in Modesto, eastern San Joaquin Valley, California

Shallow ground water in the eastern San Joaquin Valley is affected by high nitrate and uranium concentrations and frequent detections of pesticides and volatile organic compounds (VOC), as a result of ground-water development and intensive agricultural and urban land use. A single public-supply well was selected for intensive study to evaluate the dominant processes affecting the vulnerability of public-supply wells in the Modesto area. A network of 23 monitoring wells was installed, and water and sediment samples were collected within the approximate zone of contribution of the public-supply well, to support a detailed analysis of physical and chemical conditions and processes affecting the water chemistry in the well. A three-dimensional, steady-state local ground-water-flow and transport model was developed to evaluate the age of ground water reaching the well and to evaluate the vulnerability of the well to nonpoint source input of nitrate and uranium. Particle tracking was used to compute pathlines and advective travel times in the ground-water flow model. The simulated ages of particles reaching the public-supply well ranged from 9 to 30,000 years, with a median of 54 years. The age of the ground water contributed to the public-supply well increased with depth below the water table. Measured nitrate concentrations, derived primarily from agricultural fertilizer, were highest (17 milligrams per liter) in shallow ground water and decreased with depth to background concentrations of less than 2 milligrams per liter in the deepest wells. Because the movement of water is predominantly downward as a result of ground-water development, and because geochemical conditions are generally oxic, high nitrate concentrations in shallow ground water are expected to continue moving downward without significant attenuation. Simulated long-term nitrate concentrations indicate that concentrations have peaked and will decrease in the public-supply well during the next 100 years because of the low nitrate concentrations in recharge beneath the urban area and the increasing proportion of urban-derived ground water reaching the well. The apparent lag time between peak input concentrations and peak concentrations in the well is about 20 to 30 years. Measured uranium concentrations were also highest (45 micrograms per liter) in shallow ground water, and decreased with depth to background concentrations of about 0.5 microgram per liter. Naturally-occurring uranium adsorbed to aquifer sediments is mobilized by oxygen-rich, high-alkalinity water. Alkalinity increased in shallow ground water in response to agricultural development. As ground-water pumping increased in the 1940s and 1950s, this alkaline water moved downward through the ground-water flow system, mobilizing the uranium adsorbed to aquifer sediments. Ground water with high alkalinity and high uranium concentrations is expected to continue to move deeper in the system, resulting in increased uranium concentrations with depth in ground water. Because alkalinity (and correspondingly uranium) concentrations were high in shallow ground water beneath both the urban and the agricultural land, long-term uranium concentrations in the public-supply well are expected to increase as the proportion of uranium-affected water contributed to the well increases. Assuming that the alkalinity near the water table remains the same, the simulation of long-term alkalinity in the public-supply well indicates that uranium concentrations in the public-supply well will likely approach the maximum contaminant level; however, the time to reach this level is more than 100 years because of the significant proportion of old, unaffected water at depth that is contributed to the public-supply well.

California

Peak-discharge frequency and potential extreme peak discharge for natural streams in the Brazos River basin, Texas

The 2-, 5-, 10-, 25-, 50-, and 100-year peak discharges were estimated for 186 streamflow-gaging stations with at least 8 years of data for natural streams in and near the Brazos River Basin, Texas. Multiple regression equations were developed to estimate peak-discharge frequency for the 2-, 5-, 10-, 25-, 50-, and 100-year recurrence intervals for each of three hydrologic regions that compose the Brazos River Basin. The equations for each region are a function of significant basin characteristics (explanatory variables). The significant explanatory variables among six that were tested are the contributing drainage area and stream slope for regions 1 and 2 and the contributing drainage area for region 3. For the three sets of equations, the coefficient of determination ranges from 0.59 to 0.93, and the standard error ranges from 0.184 to 0.391 log units. A larger coefficient of determination and a lower standard error generally are associated with the equations for hydrologic regions 2 and 3. Statistics from the regression analysis allow computation of the prediction interval associated with a given significance level for a peak-discharge frequency estimate. The regression equations can be used to estimate peak discharges for sites at, near, or away from sites with streamflow-gaging stations. The potential extreme peak-discharge curves as related to contributing drainage area were estimated for each of the three hydrologic regions from measured extreme peaks of record at 186 sites with streamflow-gaging stations and from measured extreme peaks at 37 sites without streamflow-gaging stations in and near the Brazos River Basin. The potential extreme peak-discharge curves generally are similar for hydrologic regions 1 and 2, and the curve for region 3 consistently is below the curves for regions 1 and 2, which indicates smaller peak discharges.

Texas

Temporal variation and scale in movement-based resource selection functions

A common population characteristic of interest in animal ecology studies pertains to the selection of resources. That is, given the resources available to animals, what do they ultimately choose to use? A variety of statistical approaches have been employed to examine this question and each has advantages and disadvantages with respect to the form of available data and the properties of estimators given model assumptions. A wealth of high resolution telemetry data are now being collected to study animal population movement and space use and these data present both challenges and opportunities for statistical inference. We summarize traditional methods for resource selection and then describe several extensions to deal with measurement uncertainty and an explicit movement process that exists in studies involving high-resolution telemetry data. Our approach uses a correlated random walk movement model to obtain temporally varying use and availability distributions that are employed in a weighted distribution context to estimate selection coefficients. The temporally varying coefficients are then weighted by their contribution to selection and combined to provide inference at the population level. The result is an intuitive and accessible statistical procedure that uses readily available software and is computationally feasible for large datasets. These methods are demonstrated using data collected as part of a large-scale mountain lion monitoring study in Colorado, USA.

Statistical Methodology

Achieving interpretable machine learning by functional decomposition of black-box models into explainable predictor effects

Machine learning (ML) models are often based on complex black-box architectures that are difficult to interpret. This interpretability problem can hinder the use of ML in fields like medicine, ecology, and insurance, and has boosted research in interpretable machine learning (IML). Here, we propose a novel approach for the functional decomposition of black-box predictions, which is a core concept of IML. This approach replaces the prediction function with a surrogate model consisting of simpler subfunctions, providing insights into the direction and strength of the main feature contributions and their interactions. Our method is based on a concept termed “stacked orthogonality”, which ensures that the main effects capture as much functional behavior as possible. To compute the subfunctions, we combine neural additive modeling with an efficient post-hoc orthogonalization procedure. Our method yielded plausible results in an analysis of stream biological condition in the Chesapeake Bay watershed (United States).

Chesapeake Bay watershed

Simulation of runoff and recharge and estimation of constituent loads in runoff, Edwards aquifer recharge zone (outcrop) and catchment area, Bexar County, Texas, 1997-2000

The U.S. Geological Survey developed a watershed model (Hydrological Simulation Program—FORTRAN) to simulate runoff and recharge and to estimate constituent loads in surface-water runoff in the Edwards aquifer recharge zone (outcrop) and catchment area in Bexar County, Texas. Rainfall and runoff data collected during 1970–98 from four gaged basins in the outcrop and catchment area were used to calibrate and test the model. The calibration parameters were applied in simulations of the four calibration basins and six ungaged basins that compose the study area to obtain runoff and recharge volumes for 4 years, 1997–2000. In 1997, simulated runoff from the study area was 5.62 inches. Simulated recharge in the study area was 7.85 inches (20 percent of rainfall). In 1998, simulated runoff was 11.05 inches; simulated recharge was 10.99 inches (25 percent of rainfall). In 1999, simulated runoff was 0.66 inch; simulated recharge was 3.03 inches (19 percent of rainfall). In 2000, simulated runoff was 5.29 inches; simulated recharge was 7.19 inches (21 percent of rainfall). During 1997– 2000, direct infiltration of rainfall accounted for about 56 percent of the total Edwards aquifer recharge in Bexar County. Streamflow losses contributed about 37 percent of the recharge; flood impoundment contributed 7 percent. The simulated runoff volumes were used with event-mean-concentration data from basins in the study area and from other Bexar County basins to compute constituent loads and yields for various land uses. Annual loads for suspended solids, dissolved solids, dissolved nitrite plus nitrate nitrogen, and total lead were consistently largest from undeveloped land and smallest from commercial land or transportation corridors. Annual loads and yields varied with rainfall, with the maximum loads produced in the wettest year (1998) and the minimum loads produced in the driest year (1999).

Texas

Estimated ground-water discharge by evapotranspiration from Death Valley, California, 1997-2001

The U.S. Geological Survey, in cooperation with the National Park Service and Inyo County, Calif., collected field data from 1997 through 2001 to accurately estimate the amount of annual ground-water discharge by evapotranspiration (ET) from the floor of Death Valley, California. Multispectral satellite-imagery and National Wetlands Inventory data are used to delineate evaporative ground-water discharge areas on the Death Valley floor. These areas are divided into five general units where ground-water discharge from ET is considered to be significant. Based upon similarities in soil type, soil moisture, vegetation type, and vegetation density; the ET units are salt-encrusted playa (21,287 acres), bare-soil playa (75,922 acres), low-density vegetation (6,625 acres), moderate-density vegetation (5,019 acres), and high-density vegetation (1,522 acres). Annual ET was computed for ET units with micrometeorological data which were continuously measured at six instrumented sites. Total ET was determined at sites that were chosen for their soil- and vegetated-surface conditions, which include salt-encrusted playa (extensive salt encrustation) 0.17 feet per year, bare-soil playa (silt and salt encrustation) 0.21 feet per year, pickleweed (pickleweed plants, low-density vegetation) 0.60 feet per year, Eagle Borax (arrowweed plants and salt grass, moderate-density vegetation) 1.99 feet per year, Mesquite Flat (mesquite trees, high-density vegetation) 2.86 feet per year, and Mesquite Flat mixed grasses (mixed meadow grasses, high-density vegetation) 3.90 feet per year. Precipitation, flooding, and ground-water discharge satisfy ET demand in Death Valley. Ground-water discharge is estimated by deducting local precipitation and flooding from cumulative ET estimates. Discharge rates from ET units were not estimated directly because the range of vegetation units far exceeded the five specific vegetation units that were measured. The rate of annual ground-water discharge by ET for each ET unit was determined by fitting the annual ground-water ET for each site with the variability in vegetation density in each ET unit. The ET rate representing the midpoint of each ET unit was used as the representative value. The rate of annual ground-water ET for the playa sites did not require scaling in this manner. Annual ground-water discharge by ET was determined for all five ET units: salt-encrusted playa (0.13 foot), bare-soil playa (0.15 foot), low-density vegetation (1.0 foot), moderate-density vegetation (2.0 feet), and high-density vegetation (3.0 feet), and an area of vegetation or bare soil not contributing to ground-water discharge unclassified (0.0 foot). The total ground-water discharge from ET for the Death Valley floor is about 35,000 acre-feet and was computed by summing the products of the area of each ET unit multiplied by a corresponding ET rate for each unit.

Water-Resources Investigations Report

Streamflow conditions in the Guadalupe River Basin, south-central Texas, water years 1987-2006— An assessment of streamflow gains and losses and relative contribution of major springs to streamflow

The U.S. Geological Survey, in cooperation with the Edwards Aquifer Authority, assessed available streamflow data in the Guadalupe River Basin to determine streamflow gains and losses and the relative contribution of flow from major springs - Comal Springs, San Marcos Springs, and Hueco Springs - to streamflow in reaches of the Guadalupe River and its tributaries. The assessment is based primarily on long-term (1987-2006) and short-term (January 1999, August 1999, August 2000, and August 2006) streamflow conditions. For each analysis period, the ratio of flow from the major springs (measured at the spring source) to the sum of inflows (measured at the source of inflow to the river system) is computed for reaches of the Comal River and San Marcos River that include springflows from major springs, and for Guadalupe River reaches downstream from Canyon Dam. The ratio of springflow to the sum of inflows to the reach is an estimate of the contribution of flows from major springs to streamflow. For 1987-2006, the ratio of springflow from the major springs to the sum of inflows for the most upstream reach that includes inflow from all three major springs, Guadalupe River - above Comal River to Gonzales, is 27 percent. At the lowermost downstream reach, Guadalupe River - Bloomington to the San Antonio River, the percentage of the sum of inflows attributed to springflow is 18 percent. At that lowermost reach, the ratio of Canyon Lake releases to the sum of inflows was 20 percent. For the short-term periods August 2000 and August 2006 (periods of relatively low flow), springflow in the reach Guadalupe River - above Comal River to Gonzales accounted for 77 and 78 percent, respectively, of the sum of inflows in that reach. At the lowermost reach Guadalupe River - Bloomington to San Antonio River, springflow was 52 and 53 percent of the sum of inflows, respectively, during August 2000 and August 2006 (compared with 18 percent during 1987-2006); and during August 2000 and August 2006, the ratios of Canyon Lake releases to the sum of inflows were less than 10 percent (compared with 20 percent during 1987-2006)

Texas

Evaluating external nutrient and suspended-sediment loads to Upper Klamath Lake, Oregon, using surrogate regressions with real-time turbidity and acoustic backscatter data

Executive Summary Suspended-sediment and total phosphorus loads were computed for two sites in the Upper Klamath Basin on the Wood and Williamson Rivers, the two main tributaries to Upper Klamath Lake. High temporal resolution turbidity and acoustic backscatter data were used to develop surrogate regression models to compute instantaneous concentrations and loads on these rivers. Regression models for the Williamson River site showed strong correlations of turbidity with total phosphorus and suspended-sediment concentrations (adjusted coefficients of determination [Adj R 2 ]=0.73 and 0.95, respectively). Regression models for the Wood River site had relatively poor, although statistically significant, relations of turbidity with total phosphorus, and turbidity and acoustic backscatter with suspended sediment concentration, with high prediction uncertainty. Total phosphorus loads for the partial 2014 water year (excluding October and November 2013) were 39 and 28 metric tons for the Williamson and Wood Rivers, respectively. These values are within the low range of phosphorus loads computed for these rivers from prior studies using water-quality data collected by the Klamath Tribes. The 2014 partial year total phosphorus loads on the Williamson and Wood Rivers are assumed to be biased low because of the absence of data from the first 2 months of water year 2014, and the drought conditions that were prevalent during that water year. Therefore, total phosphorus and suspended-sediment loads in this report should be considered as representative of a low-water year for the two study sites. Comparing loads from the Williamson and Wood River monitoring sites for November 2013–September 2014 shows that the Williamson and Sprague Rivers combined, as measured at the Williamson River site, contributed substantially more suspended sediment to Upper Klamath Lake than the Wood River, with 4,360 and 1,450 metric tons measured, respectively. Surrogate techniques have proven useful at the two study sites, particularly in using turbidity to compute suspended-sediment concentrations in the Williamson River. This proof-of-concept effort for computing total phosphorus concentrations using turbidity at the Williamson and Wood River sites also has shown that with additional samples over a wide range of flow regimes, high-temporal-resolution total phosphorus loads can be estimated on a daily, monthly, and annual basis, along with uncertainties for total phosphorus and suspended-sediment concentrations computed using regression models. Sediment-corrected backscatter at the Wood River has potential for estimating suspended-sediment loads from the Wood River Valley as well, with additional analysis of the variable streamflow measured at that site. Suspended-sediment and total phosphorus loads with a high level of temporal resolution will be useful to water managers, restoration practitioners, and scientists in the Upper Klamath Basin working toward the common goal of decreasing nutrient and sediment loads in Upper Klamath Lake.

Oregon

PESTools – A Python toolkit for processing PEST-related information

PESTools is an open-source Python package for processing and visualizing information associated with the parameter estimation software PEST and PEST++. While PEST output can be reformatted for post- processing in spreadsheets or other menu-driven software packages, that approach can be error-prone and time-consuming. Managing information from highly parameterized models with thousands of parameters and observations presents additional challenges. PESTools consists of a set of Python object classes to facilitate efficient processing and visualization of PEST-related information. Processing and visualization of observation residuals, objective function contributions, parameter and observation sensitivities, parameter correlation and identifiability, and other common PEST outputs have been implemented. PESTools is integrated with the pyemu software package for linear-based computer model uncertainty analyses, allowing for efficient computations using the Jacobian Matrix without any external utilities or files. The use of dataframe objects (pandas Python package) facilitates rapid subsetting and querying of large datasets, as well as the incorporation of ancillary information such as observation locations, times, measurement types, and other associated information. PESTools’ object methods can be easily scripted with concise code, or alternatively, the use of IPython notebooks allows for live interaction with the information. PESTools is designed to streamline workflows and provide deeper insight into model behavior, enhance troubleshooting, and improve transparency in the calibration process.

Conference Paper

Simulating transport of volatile organic compounds in the unsaturated zone using the computer model R-UNSAT

Subsurface spills of gasoline and other petroleum products are a common environmental problem throughout the industrialized world. The U.S. Environmental Protection Agency has estimated that 40 percent of the more than 200,000 retail service stations in the United States have had accidental releases of petroleum hydrocarbons to the subsurface (U.S. Environmental Protection Agency, 1991). Restoration of a contaminated aquifer to regulatory standards is a technically difficult problem even when best engineering strategies are applied. Natural attenuation, a remediation strategy that relies on intrinsic physical, chemical, and biological processes to decrease contaminant concentrations, is gaining widespread acceptance in aquifer restoration efforts (Tremblay and others, 1995). The potential for successful remediation by natural attenuation depends on the fate of the organic constituents of the spilled product, which may include additives such as methyl tert-butyl ether (MTBE). These compounds can dissolve in ground water, adsorb to subsurface sediments, volatilize and diffuse through the unsaturated zone, or undergo chemical and biological reactions (fig. 1). Volatilization and biodegradation near the water table are two processes that can contribute significantly to the natural attenuation of volatile organic compounds (VOCs) in shallow ground water (McAllister and Chiang, 1994). To date, quantitative information on the rates at which these processes occur has been limited. R-UNSAT, a computer model designed for quantifying rates of volatilization and biodegradation of organic compounds near the water table, was developed and documented by the U.S. Geological Survey (USGS) and is now available to the public. R-UNSAT also can be applied, however, to other unsaturated-zone transport problems that involve gas diffusion, such as radon migration, and the deposition of compounds from the atmosphere to shallow ground water. This fact sheet describes the transport model and demonstrates its capabilities through applications to point- and nonpoint-source contamination.

Fact Sheet

Spatial and temporal variation of stream chemistry associated with contrasting geology and land-use patterns in the Chesapeake Bay watershed—Summary of results from Smith Creek, Virginia; Upper Chester River, Maryland; Conewago Creek, Pennsylvania; and Difficult Run, Virginia, 2010–2013

Despite widespread and ongoing implementation of conservation practices throughout the Chesapeake Bay watershed, water quality continues to be degraded by excess sediment and nutrient inputs. While the Chesapeake Bay Program has developed and maintains a large-scale and long-term monitoring network to detect improvements in water quality throughout the watershed, fewer resources have been allocated for monitoring smaller watersheds, even though water-quality improvements that may result from the implementation of conservation practices are likely to be first detected at smaller watershed scales. In 2010, the U.S. Geological Survey partnered with the U.S. Environmental Protection Agency and the U.S. Department of Agriculture to initiate water-quality monitoring in four selected small watersheds that were targeted for increased implementation of conservation practices. Smith Creek watershed is an agricultural watershed in the Shenandoah Valley of Virginia that is dominated by cattle and poultry production, and the Upper Chester River watershed is an agricultural watershed on the Eastern Shore of Maryland that is dominated by row-cropping activities. The Conewago Creek watershed is an agricultural watershed in southeastern Pennsylvania that is characterized by mixed agricultural activities. The fourth watershed, Difficult Run, is a suburban watershed in northern Virginia that is dominated by medium density residential development. The objective of this study was to investigate spatial and temporal variations in water chemistry and suspended sediment in these four relatively small watersheds that represent a range of land-use patterns and underlying geology to (1) characterize current water-quality conditions in these watersheds, and (2) identify the dominant sources, sinks, and transport processes in each watershed. The general study design involved two components. The first included intensive routine water-quality monitoring at an existing streamgage within each study area (including continuous water-quality monitoring as well as discrete water-quality sampling) to develop a detailed understanding of the temporal and hydrologic variability in stream chemistry and sediment transport in each watershed. The second component involved extensive water-quality monitoring at various sites throughout each watershed to develop a detailed understanding of spatial patterns. Both components were used to improve understanding of sources and transport processes affecting stream chemistry, including nutrients and suspended sediments, and their implications for detecting long-term trends related to best management practices. This report summarizes the results of monitoring that was performed from April 2010 through September 2013. Individual Small Watershed Summaries Summaries for each of the four small watersheds are presented below. Each watershed has a more descriptive and detailed section in the report, but these summaries may be particularly useful for some watershed managers and stakeholders desiring slightly less technical detail. Smith Creek Smith Creek is a 105.39-mi 2 watershed within the Shenandoah Valley that drains to the North Fork Shenandoah River. The long-term Smith Creek base-flow index is 72.3 percent, indicating that on average, approximately 72 percent of Smith Creek flow was base flow, which suggests that Smith Creek streamflow is dominated by groundwater discharge rather than stormwater runoff. A series of cluster and principal components analyses demonstrated that the majority of the variability in Smith Creek water quality could be attributed to hydrologic and seasonal variability. Statistically significant positive correlations with flow were observed for turbidity, suspended sediments, total nitrogen, ammonium, orthophosphate, iron, total phosphorus, and the ratio of calcium to magnesium. Statistically significant inverse correlations with flow were observed for specific conductance, magnesium, δ 15 N of nitrate, pH, bicarbonate, calcium, and δ 18 O of nitrate. Of particular note, flow and nitrate were not statistically significantly correlated, likely because of the relatively complex concentration-discharge relationship observed in continuous and discrete datasets. Statistically significant seasonal patterns were observed for numerous water-quality constituents: water temperature, turbidity, orthophosphate, total phosphorus, suspended-sediment concentration, and silica were higher during the warm season, but pH, dissolved oxygen, and sulfate were higher during the cool season. Surrogate regression models were developed to compute sediment and nutrient loads in Smith Creek using the continuous water-quality monitors. The mean Smith Creek in-stream sediment load was approximately 6,900 tons per year, with nearly 90 percent of the sediment load over the 3-year study period contributed during the eight largest storm events during that period. The Smith Creek total phosphorus load was approximately 21,000 pounds of phosphorus per year, with the majority of the load contributed during stormflow periods, although a substantial phosphorus load still occurs during base-flow conditions. The Smith Creek total nitrogen load was approximately 400,000 pounds per year, with total nitrogen accumulation less dominated by stormflow contributions (as was the case for sediment and total phosphorus) and strongly affected by base-flow export of nitrogen from the basin. Extensive water-quality monitoring throughout the Smith Creek watershed revealed how the complex geology and hydrology interacted to result in variable water chemistry. During relatively dry and low base-flow periods, much of the discharge in Smith Creek was contributed by a single dominant spring—Lacey Spring. During wetter base-flow periods, the flows in Smith Creek were largely generated by a mixture of headwater springs and forested mountain tributaries with very different geochemical composition. The headwater springs generally issued from limestone bedrock and were characterized as having relatively high nitrate, specific conductance, calcium, and magnesium, as well as relatively low concentrations of phosphorus, ammonium, iron, and manganese. The undeveloped, high-gradient, forested mountain sites were generally characterized by low ionic strength waters with low nutrient concentrations. Nitrate isotope data from the limestone springs generally were consistent with manure-derived nitrogen sources (such as cattle and poultry), although the possibility of other mixed sources cannot be excluded. Nitrate isotope data from the undeveloped, high-gradient forested mountain sites were more consistent with nitrogen from undisturbed soils, atmospheric deposition, or nitrogen fixation. Regardless of the nitrogen source, oxygen isotope data indicate that the nitrate was largely a result of nitrification. Land-use data indicate that manure sources of nitrogen dominated watershed nitrogen inputs. Phosphorus sources were less well studied. The presence of a single point-source discharge near the town of New Market contributed the majority of the phosphorus to Smith Creek under base-flow conditions, but nonpoint sources of phosphorus dominated the loading to Smith Creek during stormflow periods. Implementation of conservation practices increased in the Smith Creek watershed during the study period, and even though a broad range of practice types was implemented, the most common practices included stream fencing (for cattle exclusion), the development of nutrient management plans, conservation crop rotation, and the planting of cover crops. While the implementation of these conservation practices is encouraging, results indicate small increases in nitrate concentrations at the streamgage over the last 29 years, concurrent with small decreases in nitrate fluxes. It will likely be years before the cumulative effect of these practices can be detected in the Smith Creek water quality, and the magnitude of the effect of these conservation practices detected in Smith Creek will depend largely on whether nutrient loading (of manure and commercial fertilizer) is reduced over time. Upper Chester River The Upper Chester River watershed includes the 36-square-mile (mi 2 ) watershed area around several nontidal tributaries that drain into the tidal Chester River. The streamgage is on Chesterville Branch, the largest nontidal tributary (approximately 6.12 mi 2 ) and is the site for continuous water-quality monitoring for this project. The base-flow index at Chesterville Branch is about 72 percent and indicates that, as in most of the Coastal Plain, groundwater is the greatest contributor to streamflow. As such, more than 90 percent of the nitrogen in the stream is in the form of nitrate from groundwater. Continuous and discrete data collected at Chesterville Branch show the effects of streamflow and season on water quality. Significantly positive correlations with flow were observed for ammonium, dissolved and total phosphorus, sediment, and turbidity as runoff carried these constituents from the land surface into Chesterville Branch. Other constituents that increased significantly with flow include potassium, sulfate, iron, and manganese, which are likely contributed from near-stream areas and ponds with high organic-matter content. Total nitrogen, pH, and specific conductance, along with chemical constituents associated with groundwater inputs including nitrate, calcium, ratio of calcium to magnesium, silica, bicarbonate, and sodium, were negatively correlated with flow because concentrations of these constituents were diluted by runoff. Seasonal differences in water chemistry, which are most likely related to increased biologic effects on the uptake and release of chemicals in the stream and near-stream areas, also were observed. Water temperature, orthophosphate, δ 15 N of nitrate, bicarbonate, sodium, and the ratio of sodium to chloride were higher during the warm season, and dissolved oxygen, total nitrogen, nitrate, magnesium, sulfate, and manganese were higher during the cool season. Surrogate-regression models developed by using continuous water-quality data showed that the annual sediment load for the 2013 water year was about 2,600 tons, with more than 90 percent of this sediment contributed during two storms. The total phosphorus load in 2013 was about 13,000 pounds with more than 90 percent contributed during the same two storms as sediment. The load of total nitrogen, 140,000 pounds, accumulated steadily throughout the 2013 water year as nitrate in groundwater continuously discharged into the stream. The same two large storms that contributed 90 percent of the suspended-sediment and total phosphorus load only contributed about 20 percent of the annual total nitrogen load. Extensive water-quality monitoring of stream base flow throughout the Upper Chester River watershed identified how differences in land use and hydrogeology affected water chemistry. In parts of the watershed with well-drained soil and thick sandy aquifer sediments, concentrations of nitrate and other chemicals associated with fertilizer and lime application increased in streams as agricultural land use increased. More than 90 percent of the nitrogen in streams from these areas was in the form of nitrate, and concentrations ranged from about 5 milligrams per liter (mg/L) to 8 mg/L as nitrogen in the two largest tributaries. Stream nitrate concentrations were about 1 mg/L as nitrogen where soils were more poorly drained, the surficial aquifer sediments were thinner, and forests and wetlands were more widespread than agriculture. Nitrate isotope data were consistent with inorganic fertilizers ± atmospheric deposition and N 2 fixation as sources of nitrogen, and with nitrification as the dominant nitrate-forming process. Nitrate reduction was indicated by elevated δ 15 N and δ 18 O values in some samples from streams draining watersheds with poorly drained soils. An analysis of land-use data and SPARROW modeling input data attributed almost 90 percent of the nitrogen sources in the Upper Chester River watershed to inorganic fertilizer and fixation of atmospheric nitrogen by legumes, which is in agreement with the isotopic characteristics of nitrate in this watershed. Local sources of manure are limited in this area. Total phosphorus concentrations during base flow ranged from below detection to about 0.2 mg/L. Stream phosphorus concentrations during base flow were generally lower than those measured during storms because most phosphorus transport likely occurs as phosphorus attached to sediment particles during runoff. Because manure is not widely used in this area, the major source of phosphorus is likely fertilizer. The implementation of conservation practices in the Upper Chester River watershed increased substantially during the study period, with a total implementation of 1,194 U.S. Department of Agriculture-compliant practices. The most frequently used practices were oriented towards nutrient and sediment control, including cover crops, nutrient management planning, conservation crop rotation, conservation tillage, and irrigation management. The current Chesapeake Bay model for this area predicts that implementation of best management practices should result in a 13-percent decrease in overall delivery of nitrogen to the Upper Chester River. Because most nitrogen travels through the groundwater system for years to decades before being discharged to streams, the time period of monitoring was not sufficient to see the effects of these practices on water quality. The magnitude of the effect that may eventually be detected will depend on the degree to which nitrate leaching into the groundwater system is reduced over time. Loadings of phosphorus and sediment are primarily transported during large runoff events and are difficult to control and analyze for trends because of their timing and episodic nature. Conewago Creek Conewago Creek has two primary monitoring locations—one near the middle of the 47-mi 2 watershed and the other near the outlet just upstream of the Susquehanna River. The base-flow index was 47.3 percent for 2012–2013, indicating that on average, approximately 53 percent of the streamflow in Conewago Creek exited the watershed as surface flow, which suggests that the stormwater runoff was somewhat greater than groundwater discharge (base flow). A series of cluster and principal components analyses demonstrated that the majority of the variability in the Conewago Creek water quality could be attributed to hydrologic and seasonal variability. Statistically significant positive correlations with flow were observed at both monitoring sites for ammonium, total phosphorus, orthophosphate, iron, and manganese; additionally, at the upstream monitoring station, total nitrogen demonstrated a statistically significant positive correlation with flow. Statistically significant inverse correlations with flow were observed at both sites for water temperature, specific conductance (at the downstream site only), sulfate, chloride, calcium, and magnesium. Statistically significant seasonal patterns were observed for several water-quality constituents. Water temperature, phosphorus (upstream site only), and orthophosphate were higher during the warm season, and nitrate and total nitrogen (upstream site only) were higher during the cool season. Surrogate regression models were developed to compute sediment and nutrient load in Conewago Creek by using the continuous water-quality monitors and water-quality samples. Conewago Creek sediment load was approximately 9,900 tons in 2012 and approximately 18,900 tons in 2013, with nearly 80 percent of the sediment load in 2013 contributed by the three largest storm events. Annual total nitrogen loads could not be estimated due to poor model performance. The addition of continued monitoring or a continuously recording nitrate sensor could improve estimates of total nitrogen loads. During 2012 and 2013, phosphorus loads in Conewago Creek were approximately 50,000 pounds in each year. Combining data from one high-flow synoptic sampling with the data from routine sampling revealed how the geology and hydrology interact to result in variable water chemistry throughout the Conewago Creek watershed. The areas above the upstream gage in the headwaters are generally underlain by forested non-carbonate bedrock and are characterized by relatively low nitrate, specific conductance, calcium, and magnesium, as well as relatively low concentrations of phosphorus, ammonium, iron, and manganese. The more developed, agricultural areas below the upstream site were generally characterized by higher ionic strength waters with higher nutrient and metal concentrations. An analysis of land-use data and SPAtially Referenced Regressions On Watershed (SPARROW) modeling data indicates that manure sources of nitrogen dominate the input of nitrogen to the watershed. Implementation of conservation practices increased in the Conewago Creek watershed during the study period, and while a broad range of practice types were implemented, the most common practices included residue and tillage management, cover crops, nutrient management, terracing, and stream fencing (for animal exclusion or bank restoration). While the implementation of these conservation practices is encouraging, the cumulative effects of these practices probably will not be detected in Conewago Creek water quality for several years. The magnitude of the effects of these conservation practices on water quality in Conewago Creek will depend largely on the extent to which nutrient loading (septic, manure, and commercial fertilizer) and sediment-producing activities are reduced over time. Difficult Run The Difficult Run watershed is a 57.82-mi 2 watershed that drains to the Potomac River. The long-term Difficult Run base-flow index (from 1936 to 2010) was 57.9, indicating that approximately 58 percent of streamflow exited the watershed as base flow and 42 percent as stormflow; however, with continued development and urbanization of the watershed, the base-flow index has decreased to 50 percent during the last 20 years. This base-flow index was less than those of the other watersheds evaluated in this study, likely because the Difficult Run watershed largely is underlain by crystalline piedmont metamorphic rocks and has a greater proportion of impervious urban land cover. A series of cluster and principal components analyses indicated that most of the variability in Difficult Run water quality could be attributed to hydrologic variability and seasonality. Statistically significant positive correlations with flow were observed for turbidity, dissolved oxygen, suspended sediments, ammonium, orthophosphate, iron, and total phosphorus. Statistically significant inverse correlations with flow were observed for water temperature, pH, specific conductance, bicarbonate, calcium, magnesium, nitrate, δ 15 N of nitrate, and silica. Statistically significant seasonal patterns were observed for numerous water-quality constituents: water temperature, ammonium, orthophosphate, and δ 15 N of nitrate were higher during the warm season, and dissolved oxygen, nitrate, and manganese were higher during the cool season. Surrogate regression models were developed to compute sediment and nutrient loading rates. The Difficult Run sediment load was approximately 8,000 tons per year, with greater than 95 percent of the sediment load in the 2013 water year contributed by the seven largest storm events. The total phosphorus load in Difficult Run was approximately 14,000 pounds of phosphorus per year, with the majority of the load contributed during stormflow periods. The total nitrogen load in Difficult Run is estimated to have been approximately 140,000 pounds per year, with total nitrogen accumulation less dominated by stormflow contributions than that of phosphorus and strongly affected by base-flow export of nitrogen from the basin. Extensive water-quality monitoring throughout the Difficult Run watershed revealed relatively uniform generation of flow per unit of watershed area, as well as spatial variation in water quality that is strongly related to land-use activities. Elevated nitrate concentrations were observed in a subset of monitoring sites that are inversely correlated with population density and positively correlated to the septic system density within each subwatershed. The majority of the elevated nitrate concentrations for these sites are hypothesized to be caused by nitrate leaching from septic systems, more so than homeowner fertilizer usage among these subwatersheds that have lower population densities than other parts of the watershed. Nitrate isotope data, temporal patterns in the water-quality data, mass-balance computations, and a separate land-use analysis all generally indicate that leachate from septic systems was the likely source of the elevated nitrate. Another group of water-quality sites have relatively low nitrogen concentrations, are located in areas that are served by city sewer lines, and have experienced stream restoration activities. A final group of sites drained the areas with the highest imperviousness and had strongly elevated specific conductance, chloride, and sodium, which were likely caused by a combination of road salting and other anthropogenic sources draining these urbanized areas in the watershed. A fourth group of sites represents a mixture of water sources and had water quality similar to that at the Difficult Run streamgage. Analysis of the nitrate isotope data generally indicates a broad range of composition indicative of mixed natural and anthropogenic nitrogen sources. Implementation of conservation practices increased in the Difficult Run watershed during the study period, and while a broad range of practice types was implemented, the most common practices included stream restoration. While the implementation of these conservation practices is encouraging, the cumulative effect of these practices probably will not be detected in Difficult Run water quality for several years.

Maryland, Pennsylvania, Virginia

Evaluation of development options for Alaska North Slope viscous and heavy oil

Current estimates of discovered viscous and heavy oil in Alaska’s North Slope are 12 billion barrels of oil-in-place and 12–18 billion barrels of oil-in-place, respectively (see Appendix 1 for conversion to SI units). Since the early 1990s to the end of 2010, cumulative viscous oil production has amounted to 150 million barrels, and there has been no commercial production of heavy oil. During the last three decades, the industry has been challenged to develop technologies to commercially produce these untapped oil resources in this Arctic environment. In this paper, the general locations and geologic properties of the viscous oil-bearing West Sak/Schrader Bluff and heavy oil-bearing Ugnu stratigraphic intervals are described first. The geologic variability within these deposits and the evolution of technology have forced an incremental development approach, requiring costly field testing at the pilot scale of innovative extraction techniques. Although viscous oil is currently produced, its development is not mature, and firms appear to be still spending large sums on new approaches to improve recovery. The analysis specifies a representative viscous oil project and then applies a “real options” framework using simulation to determine whether the risked expected project value is sufficient to fund required expenditures on extraction process research and field testing. Computations show available field test funds to be highly sensitive to the operator’s hurdle rate of return as well as the range in magnitude of potential State revenues. The contribution of the paper is solving this problem using an approach where the extreme low return and high scenarios need only be specified, and where the uncertainties are modeled with beta distributions based on historical data or expert opinion.

Alaska

Morphology and paleohydrology of intracrater alluvial fans north of Hellas Basin, Mars

Alluvial fans and sinuous ridges are both important records of the history of fluvial activity on Mars, and they often occur together. We present observations of alluvial fans, many of which exhibit inverted relief, in five craters in the region north of Hellas basin. The observed fans ranged in size from ~10 to 820 km 2 . We identified three primary fan surface morphology classes (chute, degraded, and Inverted) as well as many instances where the morphology transitions from proximal chutes (or, rarely, a cratered degraded surface) to distal ridges corresponding to increasing thermal inertia. Clear superposition relationships at contacts between adjacent fans are rarely observed, suggesting interfingered deposits and concurrent fan development across the region. Localized factors appear to influence fan development as there is no systematic trend in the azimuth range of fan location, size of fan or catchment, as well as the degree of crater filling. Water and sediment availability may be controlled by lithology differences and weather patterns. Many of the fans had a mismatch between catchment and fan volume, corresponding to significant amounts of erosion perhaps due to windblown stripping of fine sediment. However, several notable fans exhibited volumes greater than their corresponding catchments. This may reflect uncertainty in the accuracy of the estimated paleosurface, or it may indicate sediment contributions to the fan from outside the mapped catchment. Ridges, inferred to be the resistant remnants of fluvially transported deposits, were used to estimate flow magnitude in fan construction with computed discharges of 60–400 m 3 /s and corresponding supply rate runoff values ~1–20 mm/h. Acknowledging that width-derived discharge values may overestimate flow conditions due to the likelihood of amalgamated channel deposits, this quantification provides important climate constraints. The upper range of runoff values and discharge rates are quite high, and would require either intense rain storms to generate immediate runoff, or longer-term snow accumulation and subsequent melt-runoff, potentially enhanced by rain-on-snow events. Minimum continuous formation time scales of less than a century are computed, but are incompatible with fan morphology (e.g., superposition relationships, embedded craters) and mechanisms to sustain flows. More realistic lower-limit fan construction times, accounting for modeled precipitation rates from the literature, are tens to hundreds of thousands of years. Fans were active in multiple events spanning the Hesperian to Amazonian periods, requiring transient climate conditions to support the fan aggradation.

Icarus