Search USGSSearch

SEARCH · Search USGS

Results for “General Information Product”

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 883 records · Page 49Linked to original sources

Evaluation of water quality and bulk atmospheric deposition in the Guanella Pass area, Clear Creek and Park Counties, Colorado, water year 1995

A study of water quality and bulk atmospheric deposition in the Guanella Pass area during water year 1995 (October 1994 to September 1995) has provided information to characterize conditions in the South Clear Creek and Geneva Creek Basins about 40 miles west of Denver, Colorado, prior to a possible roadreconstruction project. Analytical results of water samples collected from streams, ground water, and lakes and reservoirs indicate that water in the study area generally contained small concentrations of dissolved solids (low-flow stream sample median, 43 milligrams per liter; ground-water median, 104 milligrams per liter) and suspended sediment (high-flow stream sample median, 115 milligrams per liter). Nutrient, trace-element, and organic-carbon concentrations also were small. Acidic pH and large trace-element concentrations were associated with areas affected by natural acid-sulfate weathering or mining. Increases in particulate nutrients, total recoverable trace elements, and suspended sediment were related to seasonal snowmelt and thunderstorms. Benthic macroinvertebrate and periphyton (algae) samples collected from selected sites in August 1995 indicated large variability in characteristics among sites. Macroinvertebrate density ranged from 6.5 to 1,300 organisms per square meter, and periphyton biovolume ranged from 0.023 to 4,700 x 10 6 cubic microns per square centimeter. The discharge of road runoff from any single drainage feature (ditch or culvert) was small compared to streamflows. Analytical results from samples of road runoff indicated a large range in specific conductance, pH in a neutral range, large values of turbidity, and calcium bicarbonate type water. Sodium or magnesium and chloride, however, were predominant in some samples and were probably related to applications of sodium and magnesium chloride to the roads. Although dissolved nutrients and trace elements were detected in samples of road runoff, particulate forms had the largest concentrations. Suspended-sediment concentrations in road runoff were large, and the sediment was primarily fine grained. Large concentrations of dissolved and total recoverable nutrients, total recoverable trace elements, total organic carbon, and suspended sediment in road runoff could cause increased concentrations in streams, lakes, and reservoirs. Ground-water concentrations of chloride, nitrite plus nitrate, and dissolved iron and manganese·could increase if the large concentrations present in road runoff infiltrate shallow aquifers. Bulk atmospheric deposition was collected at three sites with a paved road surface, five sites with a gravel or dirt surface, and five sites with collectors at least 500 feet away from a road. The median bulk -atmospheric-deposition rate for gravel (or dirt) roads was 106 times the median rate at a reference site (located at least 500 feet away from a road) and 39 times the median rate at a typical paved road site. If a stream, lake, or reservoir is located near the road, bulk-atmospheric- deposition products may settle directly onto open water or be washed in by surface runoff. Because the predominant road-runoff trace-element phase is particulate, potential toxic effects on receiving water may be minimal when compared to aquatic-life standards and regulations in the State of Colorado, which are primarily based on dissolved concentrations. U.S. Environmental Protection Agency primary maximum contaminant levels for drinking water were exceeded for dissolved uranium (proposed) in one ground-water sample.

Colorado

Snake River Fall Chinook Salmon research and monitoring

In Chapter 1, we report on development and application of an integrated population model (IPM) for the natural-origin fall Chinook salmon population upstream of Lower Granite Dam. This year’s efforts represent the third update to the model. Initial efforts focused on generating juvenile and adult abundance estimates, with estimates of uncertainty, for informing the life-cycle model and estimating the effects of covariates on key demographic parameters. The goals of this year’s report are to 1) describe the modifications and advances made since the previous report, 2) to annually update and report the abundance estimates and other quantities used in the model, 3) to provide annual estimates of population parameters estimated by the IPM, and 4) to outline the next year’s tasks for advancing and/or applying the model. Since our last report on the life-cycle model, we have made a number of changes including: 1) incorporating jack abundance and age-structure data into the observation model, 2) changing smolt-to-adult survival (SAR) for subyearling and yearling to partial SARs that represent the joint probability surviving and entering the ocean at a given juvenile age, 3) combining age categories for rarely observed ages, 4) using scale data from unmarked fish to estimate age structure, and 5) generating composite life-cycle demographic parameters (cumulative capacity and productivity) from stage-specific parameters. We also generated juvenile abundance estimates, extended the model to include three additional brood years (1992– 2021), and ran the model to forecast returns to Lower Granite in 2022. For posterior medians of life stage-specific parameters, we estimated a mean productivity of 438 natural-origin juvenile recruits per female spawner, a capacity of 1.36 million juveniles, and a mean smolt-to-adult survival (SAR) of 1.2%. We detected strong density-dependent regulation, with juvenile recruits per spawner declining to about 50 juvenile recruits per female spawner at high spawner abundance. Across the entire life cycle, these stage-specific parameters resulted in a median cumulative intrinsic productivity of 1.93 adult female recruits per female spawner and a median equilibrium abundance of 2,851 female spawners (7,842 total spawners). Annual juvenile productivity varied from about 250–1,000 juveniles per spawner but displayed no temporal trends or patterns. For the three most recent brood years added to the model, recruits per spawner were higher than average but well within the range of uncertainty observed over the entire time series. In contrast to juvenile recruitment variability, SAR varied considerably among years and exhibited two periods of high survival (1996–2001 and 2007–2012) when SAR ranged from 2% to 6% and cumulative productivity ranged from 2 to 8 recruits per spawner. Partial SARs revealed that yearling outmigrants contributed substantially to the high SARs in the first high-survival period, but the second period was dominated by subyearlings. Yearlings contributed >30% to SAR in most years prior to 2007, and <30% since 2007. Our two-stage IPM provides a wealth of information about population dynamics affecting two key life-stage transitions (spawner to juvenile, and juvenile to spawner) centered on passage at Lower Granite Dam. By summarizing these stage-specific demographic parameters across the entire life cycle, this information will be useful for informing the recovery status of this threatened population. Whereas previous versions introduced hydrosystem and ocean covariates into the model, this phase of model development focused on solidifying the underlying model structure by introducing the concept of partial SARs and developing composite productivity and capacity as a function of underlying stage-specific parameters. Given this advancement, our next steps are to re-incorporate covariates into the model, specifically to understand how different factors affect partial SARs of subyearling and yearlings. Longer term model developments include:1) incorporating hatchery fish to explicitly estimate their survival as an alternative method for estimating natural-origin age composition, 2) expanding the model’s structure to include the three major spawning aggregates, 3) more explicitly modeling hydrosystem effects including transportation, and 4) using the model to assess retrospective and prospective management actions. In 2022, the U.S. Geological Survey (USGS) focused adult salmon survey efforts in the Snake River on deepwater redd searches and fish collection for parentage-based tagging (PBT) analyses. We use used a boat-mounted underwater video camera to count 99 deepwater redds at 16 of the 29 sites surveyed. Redd depths averaged 4.4 m. In conjunction with the Idaho Power Company, we collected genetic samples from 318 live fall Chinook salmon (Oncorhynchus tshawytscha) and 19 carcasses at 40 unique geographic locations that spanned 91 river kilometers. Eighty fish were collected at three sites (High Range [rkm 332.3], Dug Bar [rkm 315.4], and Three Creek [rkm 384.0]), which accounted for 23% of all collected fish in 2022. Most (333 fish) post-spawned salmon were collected from early to mid-November just after the peak of spawning. A summary of 2021 PBT results produced by the Idaho Power Company can be found in Appendix A.2. Beach seining and PIT tagging of subyearling fall Chinook salmon was conducted in Snake and Salmon rivers to obtain information on population metrics and growth as well as to provide data for ongoing life-cycle modeling. In the Snake River, we collected 7,496 subyearlings, tagged 4,139, and recaptured 502 (12.1%). Using 8-mm tags in 45–49-mm fish allowed us to represent an additional 25% of the juvenile population through PIT tagging beyond just using standard 9- and 12-mm tags. In the Salmon River, we captured 206 natural subyearlings with the majority (52%) of fish being captured at two sites: rkm 20 and 26. We tagged 145 subyearlings and recaptured 9 fish. Many of the subyearlings we tagged in the Snake River were detected passing Lower Granite Dam, but only 4 fish tagged in the Salmon River were detected. In total we detected 484 (11.3%) tagged fish at Lower Granite Dam, and detection rates varied by tag size and passage route. More subyearlings were detected passing via the removable spill weir (RSW) earlier in the season while more fish were detected passing through the juvenile fish bypass system (JBS) earlier in the season while more fish were detected passing via the removable spill weir (RSW) later in the season. In general, fish tagged with 12-mm PIT tags had higher detection rates than fish tagged with smaller tags. Survival to Lower Granite Dam was low and ranged from 0.22 to 0.36. Season-wide, growth of subyearlings was higher in the lower reach than in the upper reach of the Snake River.

Idaho, Oregon, Washington

Geologic and hydrogeologic frameworks of the Biscayne aquifer in central Miami-Dade County, Florida

Evaluations of the lithostratigraphy, lithofacies, paleontology, ichnology, depositional environments, and cyclostratigraphy from 11 test coreholes were linked to geophysical interpretations, and to results of hydraulic slug tests of six test coreholes at the Snapper Creek Well Field (SCWF), to construct geologic and hydrogeologic frameworks for the study area in central Miami-Dade County, Florida. The resulting geologic and hydrogeologic frameworks are consistent with those recently described for the Biscayne aquifer in the nearby Lake Belt area in Miami-Dade County and link the Lake Belt area frameworks with those developed for the SCWF study area. The hydrogeologic framework is characterized by a triple-porosity pore system of (1) matrix porosity (mainly mesoporous interparticle porosity, moldic porosity, and mesoporous to megaporous separate vugs), which under dynamic conditions, produces limited flow; (2) megaporous, touching-vug porosity that commonly forms stratiform groundwater passageways; and (3) conduit porosity, including bedding-plane vugs, decimeter-scale diameter vertical solution pipes, and meter-scale cavernous vugs. The various pore types and associated permeabilities generally have a predictable vertical spatial distribution related to the cyclostratigraphy. The Biscayne aquifer within the study area can be described as two major flow units separated by a single middle semiconfining unit. The upper Biscayne aquifer flow unit is present mainly within the Miami Limestone at the top of the aquifer and has the greatest hydraulic conductivity values, with a mean of 8,200 feet per day. The middle semiconfining unit, mainly within the upper Fort Thompson Formation, comprises continuous to discontinuous zones with (1) matrix porosity; (2) leaky, low permeability layers that may have up to centimeter-scale vuggy porosity with higher vertical permeability than horizontal permeability; and (3) stratiform flow zones composed of fossil moldic porosity, burrow related vugs, or irregular vugs. Flow zones with a mean hydraulic conductivity of 2,600 feet per day are present within the middle semiconfining unit, but none of the flow zones are continuous across the study area. The lower Biscayne aquifer flow unit comprises a group of flow zones in the lower part of the aquifer. These flow zones are present in the lower part of the Fort Thompson Formation and in some cases within the limestone or sandstone or both in the uppermost part of the Pinecrest Sand Member of the Tamiami Formation. The mean hydraulic conductivity of major flow zones within the lower Biscayne aquifer flow unit is 5,900 feet per day, and the mean value for minor flow zones is 2,900 feet per day. A semiconfining unit is present beneath the Biscayne aquifer. The boundary between the two hydrologic units is at the top or near the top of the Pinecrest Sand Member of the Tamiami Formation. The lower semiconfining unit has a hydraulic conductivity of less than 350 feet per day. The most productive zones of groundwater flow within the two Biscayne aquifer flow units have a characteristic pore system dominated by stratiform megaporosity related to selective dissolution of an Ophiomorpha-dominated ichnofabric. In the upper flow unit, decimeter-scale vertical solution pipes that are common in some areas of the SCWF study area contribute to high vertical permeability compared to that in areas without the pipes. Cross-hole flowmeter data collected from the SCWF test coreholes show that the distribution of vuggy porosity, matrix porosity, and permeability within the Biscayne aquifer of the SCWF is highly heterogeneous and anisotropic. Groundwater withdrawals from production well fields in southeastern Florida may be inducing recharge of the Biscayne aquifer from canals near the well fields that are used for water-management functions, such as flood control and well-field pumping. The SCWF was chosen as a location within Miami-Dade County to study the potential for such recharge to the Biscayne aquifer from the C–2 (Snapper Creek) canal that roughly divides the well field in half. Geologic, hydrogeologic, and hydraulic information on the aquifer collected during construction of monitoring wells within the SCWF could be used to evaluate the groundwater flow budget at the well-field scale.

Florida

Depositional sequence stratigraphy and architecture of the cretaceous ferron sandstone: Implications for coal and coalbed methane resources - A field excursion

This Field Excursion will visit outcrops of the fluvial-deltaic Upper Cretaceous (Turonian) Ferron Sandstone Member of the Mancos Shale, known as the Last Chance delta or Upper Ferron Sandstone. This field guide and the field stops will outline the architecture and depositional sequence stratigraphy of the Upper Ferron Sandstone clastic wedge and explore the stratigraphic positions and compositions of major coal zones. The implications of the architecture and stratigraphy of the Ferron fluvial-deltaic complex for coal and coalbed methane resources will be discussed. Early works suggested that the southwesterly derived deltaic deposits of the the upper Ferron Sandstone clastic wedge were a Type-2 third-order depositional sequence, informally called the Ferron Sequence. These works suggested that the Ferron Sequence is separated by a type-2 sequence boundary from the underlying 3rd-order Hyatti Sequence, which has its sediment source from the northwest. Within the 3rd-order depositional sequence, the deltaic events of the Ferron clastic wedge, recognized as parasequence sets, appear to be stacked into progradational, aggradational, and retrogradational patterns reflecting a generally decreasing sediment supply during an overall slow sea-level rise. The architecture of both near-marine facies and non-marine fluvial facies exhibit well defined trends in response to this decrease in available sediment. Recent studies have concluded that, unless coincident with a depositional sequence boundary, regionally extensive coal zones occur at the tops of the parasequence sets within the Ferron clastic wedge. These coal zones consist of coal seams and their laterally equivalent fissile carbonaceous shales, mudstones, and siltstones, paleosols, and flood plain mudstones. Although the compositions of coal zones vary along depositional dip, the presence of these laterally extensive stratigraphic horizons, above parasequence sets, provides a means of correlating and defining the tops of depositional parasequence sets in both near-marine and non-marine parts of fluvial-deltaic depositional sequences. Ongoing field studies, based on this concept of coal zone stratigraphy, and detailed stratigraphic mapping, have documented the existence of at least 12 parasequence sets within the Last Chance delta clastic wedge. These parasequence sets appear to form four high frequency, 4th-order depositional sequences. The dramatic erosional unconformities, associated with these 4th-order sequence boundaries, indicate that there was up to 20-30 m of erosion, signifying locally substantial base-level drops. These base-level drops were accompanied by a basin ward shift in paleo-shorelines by as much as 5-7 km. These 4th-order Upper Ferron Sequences are superimposed on the 3rd-order sea-level rise event and the 3rd-order, sediment supply/accommodation space driven, stratigraphie architecture of the Upper Ferron Sandstone. The fluvial deltaic architecture shows little response to these 4th-order sea-level events. Coal zones generally thicken landward relative to the mean position of the landward pinch-out of the underlying parasequence set, but after some distance landward, they decrease in thickness. Coal zones also generally thin seaward relative to the mean position of the landward pinch-out of the underlying parasequence set. The coal is thickest in the region between this landward pinch-out and the position of maximum zone thickness. Data indicate that the proportion of coal in the coal zone decreases progressively landward from the landward pinch-out. The effects of differential compaction and differences in original pre-peat swamp topography have the effect of adding perturbations to the general trends. These coal zone systematics have major impact on approaches to exploration and production, and the resource accessment of both coal and coalbed methane.

Brigham Young University Geology Studies

Hydrogeological framework, numerical simulation of groundwater flow, and effects of projected water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma

This report describes a study of the hydrology, hydrogeological framework, numerical groundwater-flow models, and results of simulations of the effects of water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma. The purpose of the study was to provide analyses, including estimating equal-proportionate-share (EPS) groundwater-pumping rates and the effects of projected water use and droughts, pertinent to water management of the Beaver-North Canadian River alluvial aquifer for the Oklahoma Water Resources Board. The Beaver-North Canadian River alluvial aquifer consists of unconsolidated sand, gravel, silt, and clay in varying proportions that underlies the Beaver and North Canadian River Valleys for approximately 175 miles (mi) from the Oklahoma Panhandle to the western edge of Oklahoma City in central Oklahoma. The aquifer as delineated for this study varies from 4 to 12 mi wide and is as thick as 308 feet (ft) in the northwest where the aquifer includes the Ogallala Formation. There are two distinct but in most areas hydraulically connected alluvial units that compose the Beaver-North Canadian River alluvial aquifer: a Quaternary-age topographically higher terrace deposit and a topographically lower, younger alluvium along the active river channel that includes active and Quaternary-age alluvium. The Beaver River composes the headwaters of the North Canadian River, which begins at the confluence of the Beaver River and Wolf Creek. The aquifer is divided for water management into two geographic areas: Reach I upstream from Canton Dam and Reach II downstream from Canton Dam. Reach I covers an area of approximately 874 square miles (mi 2 ), and Reach II covers an area of approximately 371 mi 2 . The Beaver-North Canadian River alluvial aquifer crosses several climatic zones, from semiarid in the west to continental subhumid in the east. Mean annual precipitation varies from 23.5 inches (in.) in the western part of this aquifer to 35.7 in. in the east. Surface-water demands were met through numerous temporary and permanent surface-water diversions from the Beaver and North Canadian Rivers during the period of study. During the study period, seven diversions removed a mean annual 2,000 acre-feet (acre-ft) of water from Reach I. There were 14 diversions from Reach II with a mean annual permitted volume of approximately 81,000 acre-ft, including diversion into the Lake Hefner Canal for the Oklahoma City public water supply. During the period of this study, 17 temporary surface-water diversion permits were active in Reach I, with total permitted volumes of 2,000 acre-ft, and 41 diversions were active in Reach II, with total permitted volumes of 38,000 acre-ft. The total water use for each temporary permit was assumed to be taken over the 3-month period allotted to temporary withdrawal permits. The groundwater-use analysis full period of record, 1967&ndash;2011, was divided into two sub-intervals because of varying water use, 1970&ndash;80 and 1981&ndash;2011. Groundwater use in Reach I and Reach II was substantially greater from 1970 to 1980 compared to the rest of the period, and the sub-period 1981&ndash;2011 was used because this period includes recent population growth and modern irrigation methods. The total mean annual groundwater use in Reach I was 15,309 acre-feet per year (acre-ft/yr) during 1967&ndash;2011; 20,724 acre-ft/yr during 1970&ndash;80, and 13,739 acre-ft/yr during 1981&ndash;2011. Total mean annual groundwater use in Reach II was similar but slightly less than in Reach I, with 14,098 acre-ft/yr during 1967&ndash;2011; 19,963 acre-ft/yr during 1970&ndash;80; and 12,285 acre-ft/yr during 1981&ndash;2011. Irrigation composed 72 percent of groundwater use in Reach I and 48 percent of groundwater use in Reach II during the 1967&ndash;2011 period. Public water supply was a much smaller proportion of total groundwater use in Reach I (15 percent) than in Reach II (39 percent). The proportion of groundwater use for power was 10 percent in Reach I and 5.2 percent in Reach II. All other water-use categories in Reach I only composed 2.2 percent of groundwater use in Reach I. In Reach II, industrial, mining, and commercial categories combined accounted for 4.4 percent of groundwater use; recreation, fish, and wildlife groundwater use accounted for 2.3 percent; and nonirrigated agriculture accounted for 1.5 percent of groundwater use. Permian-age bedrock underlies the Beaver-North Canadian River alluvial aquifer. In the east, the Dog Creek Shale, the Duncan Sandstone, and the Blaine and Chickasha Formations, none of which are notable sources of groundwater in the study area, underlie the Beaver-North Canadian River alluvial aquifer. In the northwestern part of Reach I, bedrock is composed of the Rush Springs and Marlow Formations, which are productive aquifers in some areas. The Cloud Chief Formation is not a source of groundwater. One hydrogeological unit was delineated in the Beaver-North Canadian River alluvial aquifer, composed of the terrace deposits and alluvium, with limited flow between this unit and bedrock units. Groundwater in this aquifer generally flows from northwest to southeast and across the aquifer toward the Beaver and North Canadian Rivers. Groundwater recharge from precipitation was estimated for the entire Beaver-North Canadian River alluvial aquifer and then itemized for both reaches by using a soil-water-balance (SWB) model. At two locations in Reach I, a water-table fluctuation method was used to estimate local recharge. Total mean annual groundwater recharge from the soil-water-balance method was estimated to be approximately 136,400 acre-ft in Reach I and 82,400 acre-ft in Reach II; the mean annual recharge for both reaches combined was approximately 218,800 acre-ft. Two sites in Reach I located at observation wells with continuous water-level measurements and nearby streamflow-gaging stations with precipitation gages were used to estimate the percentage of precipitation that becomes groundwater recharge. The Woodward site was located at observation well OW-4 near the Woodward, Okla. (07237500), streamflow-gaging station. Total precipitation and recharge for the Woodward and Seiling sites were calculated for the water year 2013. The Woodward site had a total of 14.18 in. of precipitation and 6.3 in. of recharge was calculated, equaling 44 percent of precipitation. The mean percentage of precipitation that was estimated to become recharge in the SWB model for the period 1980&ndash;2011 at that location was 9.2 percent, although adjacent SWB-model cells were as high as 20 percent of precipitation. The Seiling site had a total of 26.84 in. of precipitation during the water year 2013, and a total of 6.9 in. of recharge was estimated, equaling 25.9 percent of precipitation. At the Seiling site, the mean percentage of precipitation that became recharge in the SWB model for the period 1980&ndash;2011 was 23.0 percent. The principal inflow to the Beaver-North Canadian River alluvial aquifer was estimated to be surface recharge from precipitation, and plant evapotranspiration was estimated to be the greatest discharge, followed by stream and lake base flow, groundwater pumping, and flow to seeps and springs along the eastern margin of the aquifer. Reach I also included inflow from the High Plains aquifer as lateral inflow of groundwater, though this flow was estimated to be a very minor component of the total water budget. Most of the Beaver and North Canadian Rivers were determined to be gaining streamflow from groundwater, but several reaches in Reach I upstream from Wolf Creek were determined to be losing streamflow through infiltration to the aquifer. Aquifer hydrogeologic characteristics were estimated from borehole lithologic logs, well-construction information, and published aquifer tests and during numerical model calibration. The maximum saturated aquifer thickness in Reach I was estimated to be 308 ft, and the mean thickness was estimated to be 36 ft. The maximum saturated thickness in Reach II was estimated to be 86 ft, and the mean thickness was estimated to be 29 ft. Mean hydraulic conductivity of Reach I was estimated to be 70 feet per day (ft/d) with a range of 7&ndash;279 ft/d. Mean hydraulic conductivity in Reach II was estimated to be 92 ft/d with a range of 4&ndash;279 ft/d. Both reach models were calibrated manually by using trial-and-error adjustment of recharge, hydraulic conductivity, specific yield, and conductance of boundary conditions. The Reach I model used 28 head observations during the steady-state period of 1980 and 487 head observations during the transient period of 1981&ndash;2011. The root-mean-square error of head residuals (observed minus simulated head) was 3.86 ft, and 83 percent of head residuals were between -5 and 5 ft. The Reach II model was calibrated to 75 steady-state head observations and 134 head observations during the transient period. The root-mean-square error of head residuals for that reach was 3.58 ft, and similar to Reach I, 85 percent of residuals were between -5 and 5 ft. Several analyses were performed by using the numeric groundwater-flow models as predictive tools, including estimating the EPS pumping rate for both reaches. The EPS is defined by the Oklahoma Water Resources Board as an annual per-acre groundwater-pumping rate that will reduce saturated thickness in half of the aquifer to 5 ft or less over a period of 20 years; additional estimates were made for periods of 40 and 50 years. Other analyses included using models to estimate the effects of groundwater pumping and a prolonged drought on groundwater in storage and streamflow and lake storage of water. The EPS pumping rate was found to be approximately 0.57 acre-feet per acre per year ([acre-ft/acre]/yr) in Reach I and 0.73 (acre-ft/acre)/yr in Reach II for a 20-year period. For a 40-year period, the annual EPS pumping rate was determined to be 0.54 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. For a 50-year period, the EPS pumping rate was determined to be 0.53 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. Groundwater pumping at the 2011 rate for 50 years resulted in a 3.6-percent decrease in the amount of water in groundwater storage in Reach I and a decrease of 2.5 percent in the amount of groundwater in storage in Reach II. A cumulative 32-percent increase in pumping greater than the 2011 rate over a period of 50 years caused a decrease in groundwater storage of 4.0 percent in Reach I and 3.3 percent in Reach II. A hypothetical severe drought was simulated by using aquifer recharge flow rates during the drought year of 2011 for a period of 10 years. All other flows including evapotranspiration and groundwater pumping were set at estimated 2011 rates. The hypothetical drought caused a decrease in water in aquifer storage by about 7 percent in Reach I and 7 percent in Reach II. Another analysis of the effects of hypothetical drought estimated the effects of drought on streamflow and lake storage. The hypothetical drought was simulated by decreasing recharge by 75 percent for a selected 10-year period (1994&ndash;2004) during the 1980&ndash;2011 simulation. In Reach I, the amounts of water stored in Canton Lake and streamflow at the Seiling, Okla., streamflow-gaging station were analyzed. Streamflow at the Seiling station decreased by a mean of 75 percent and was still diminished by 10 percent after 2011. In Reach II, the effect of drought on the streamflow at the Yukon, Okla., streamflow-gaging station was examined. The greatest mean streamflow decrease was approximately 60 percent during the simulated drought, and after 2011, the mean decrease in streamflow was still about 5 percent. Canton Lake storage decreased by as much as 83 percent during the simulated drought and did not recover by 2011.

Oklahoma

Iron-ore resources of the United States including Alaska and Puerto Rico, 1955

The importance of iron ore, the basic raw material of steel, as a fundamental mineral, resource is shown by the fact that about 100 million long tons of steel is used annually in the economy of the United States, as compared with a combined total of about 5 million long tons of copper, lead, zinc, and aluminum. Satisfying this annual demand for steel requires about 110 million tons of iron ore and 70 million tons of scrap iron and steel. The average annual consumption of iron ore in the United States from 1951 to 1955, inclusive, was about 110 million long tons, which is about twice the annual average from 1900 to 1930. Production of iron ore in the United States in this 5-year period averaged approximately 100 million long tons annually, divided by regions as follows (in percent): Lake Superior, 84.1; southeastern, 6.7; western, 6.7; northeastern, 1.4; and central and gulf, 1.1. Mining of iron ore began in the American Colonies about 1619, and for 225 years it was limited to eastern United States where fuel and markets were readily available. Production of iron ore from the Lake Superior region began in 1846; the region became the leading domestic source by 1890, and the Mesabi range in Minnesota has been the world's most productive area since 1896. Proximity of raw materials, water transportation, and markets has resulted in centralization of the country's iron and steel industry in the lower Great Lakes area. Increased imports of iron ore being delivered to eastern United States as well as demands for steel in nearby markets have given impetus to expansion in the steel-making capacity in this area. The four chief iron-ore minerals - hematite, liminite, magnetite, and siderite - are widely distributed but only locally form deposits of sufficient tonnage and grade to be commercially valuable at the present time. The iron content of these minerals, of which hematite is the most important, ranges from 48 percent in siderite to 72 percent in magnetite, but as these minerals are associated with other rock-forming minerals, the iron content of marketable ore has a lower range from 30 to 67 percent. Chemical constituents other than iron also are important in determining the marketability of iron ore. Although some iron ores can be used in the blast furnace as mined, others must first be improved either chemically by reduction of undesirable constituents, or physically by aggregation. Phosphorus and sulfur particularly are common deleterious elements; excessive silica is also undesirable but within certain limits can be controlled by additional flux. Lime and magnesia are beneficial in specified amounts because of their fluxing qualities, and a small amount of alumina improves the fluidity of slag. Manganese is especially desirable as a deoxidizing and desulfurizing agent. Titanium, chromium, and nickel must also be considered in the use of ore containing these elements. The principal iron-ore deposits in the United States have been formed by three processes. Hematite-bearing bedded deposits such as those at Birmingham, Ala., are marine sedimentary rocks which, except for weathering along the outcrop, have remained practically unaltered since deposition. Deposits of the Lake Superior region, also in sedimentary strata, originally had a slightly lower iron content than those at-Birmingham, but ore bodies of hematite and limonite were formed by removal of other constituents in solution after deposition of the beds, with a relative increase of iron content in the material remaining. Limestone adjacent to igneous intrusions has been replaced by magnetite deposits at Cornwall, Pa., and by hematite-magnetite deposits near Cedar City, Utah. Magnetite deposits in New Jersey and in the Adirondack Mountains of New York are generally believed to have been formed by replacement of grains of other minerals in metamorphic rocks. Iron-ore resources are made up of reserves of iron ore, material usable under existing economic and technologic conditions; and potential ore, material likely to become usable under more favorable conditions. The tonnage and grade of material of combined reserves and potential ore in each of the deposits known or believed to contain at least 200,000 long tons of iron-ore resources are tabulated in this report, and numerous sources of additional information are given in a selected bibliography. The total domestic iron-ore resources are estimated at approximately 75,000 million long tons of crude ore. About 10,000 million tons of the resources is reserves of crude ore that will probably yield 5,500 million tons of concentrates and direct-shipping ore. About 65,000 million tons is potential ore and may yield 25,000 million tons of concentrates and some direct-shipping ore.

Bulletin

Water Quality Conditions in Upper Klamath and Agency Lakes, Oregon, 2006

The U.S. Geological Survey Upper Klamath Lake water quality monitoring program gathered information from multiparameter continuous water quality monitors, physical water samples, dissolved oxygen production and consumption experiments, and meteorological stations during the June-October 2006 field season. The 2006 study area included Agency Lake and all of Upper Klamath Lake. Seasonal patterns in water quality were similar to those observed in 2005, the first year of the monitoring program, and were closely related to bloom dynamics of the cyanobacterium (blue-green alga) Aphanizomenon flos-aquae (AFA) in the two lakes. High dissolved oxygen and pH conditions in both lakes before the bloom declined in July, which coincided with seasonal high temperatures and resulted in seasonal lows in dissolved oxygen and decreased pH. Dissolved oxygen and pH in Upper Klamath and Agency Lakes increased again after the bloom recovered. Seasonal low dissolved oxygen and decreased pH coincided with seasonal highs in ammonia and orthophosphate concentrations. Seasonal maximum daily average temperatures were higher and minimum dissolved oxygen concentrations were lower in 2006 than in 2005. Conditions potentially harmful to fish were influenced by seasonal patterns in bloom dynamics and bathymetry. Potentially harmful low dissolved oxygen and high un-ionized ammonia concentrations occurred mostly at the deepest sites in the Upper Klamath Lake during late July, coincident with a bloom decline. Potentially harmful pH conditions occurred mostly at sites outside the deepest parts of the lake in July and September, coincident with a heavy bloom. Instances of possible gas bubble formation, inferred from dissolved oxygen data, were estimated to occur frequently in shallow areas of Upper Klamath and Agency Lakes simultaneously with potentially harmful pH conditions. Comparison of the data from monitors in nearshore areas and monitors near the surface of the water column in the open waters of Upper Klamath Lake revealed few differences in water quality dynamics. Median daily temperatures were higher in nearshore areas, and dissolved oxygen concentrations were periodically higher as well during periods of high AFA bloom. Differences between the two areas in water quality conditions potentially harmful to fish were not statistically significant (p < 0.05). Chlorophyll a concentrations varied temporally and spatially throughout Upper Klamath Lake. Chlorophyll a concentrations indicated an algal bloom in late June and early July that was followed by an algae bloom decline in late July and early August and a subsequent recovery in mid-August. Sites in the deepest part of the lake, where some of the highest chlorophyll a concentrations were observed, were the same sites where the lowest dissolved oxygen concentrations and the highest un-ionized ammonia concentrations were recorded during the bloom decline, indicating cell senescence. Total phosphorus concentrations limited the initial algal bloom in late June and early July. The rate of net dissolved oxygen production (that is, production in excess of community respiration) and consumption (due to community respiration) in the lake water column as measured in light and dark bottles, respectively, ranged from 2.79 to -2.14 milligrams of oxygen per liter per hour. Net production rate generally correlated positively with chlorophyll a concentration, except episodically at a few sites where high chlorophyll a concentrations resulted in self-shading that inhibited photosynthesis. The depth of photic zone was inversely correlated with chlorophyll a concentration. Calculations of a 24-hour change in dissolved oxygen concentration indicated that oxygen-consuming processes predominated at the deep trench sites and oxygen-producing processes predominated at the shallow sites. In addition, calculations of the 24-hour change in dissolved oxygen indicate that oxygen-consuming processes in the water column di

Scientific Investigations Report

Preliminary hydrogeologic appraisal of selected hydrothermal systems in northern and central Nevada

Several hydrothermal systems in northern and central Nevada were explored in a hydrogeologic reconnaissance. The systems studied comprise those at Stillwater and Soda Lakes-Upsal Hogback in the Carson Desert, Gerlach, Fly Ranch-Granite Range, and Double Hot Springs in the Black Rock Desert, Brady's Hot Springs, Leach Hot Springs in Grass Valley, Buffalo Valley Hot Springs, and Sulphur Hot Springs in Ruby Valley. The investigation focused on (1) delineating of areas of high heat flow associated with rising thermal ground water, (2) determining the nature of the discharge parts of the hydrothermal systems, (3) estimating heat discharge from the systems, (4) estimating water discharge from the systems, (5) obtaining rough estimates of, conductive heat flow outside areas of hydrothermal discharge, and (6) evaluating several investigative techniques that would yield the required information quickly and at relatively low cost. The most useful techniques were shallow test drilling to obtain geologic, hydraulic, and thermal data and hydrogeologic mapping of the discharge areas. The systems studied are in the north-central part of the Basin and Range province. Exposed volcanic rocks of latest Tertiary and Quaternary age are chiefly basaltic. Basaltic terranes are generally regarded as less favorable for geothermal resources than terranes that contain large volumes of young volcanic mocks of felsic to intermediate composition. Most of the known hydrothermal systems are associated with Basin and Range faults which are caused by crustal extension across the province. An area of high heat flow centered at Battle Mountain and possibly other areas of high heat flow may be related to crustal heat sources. However, some of the hydrothermal systems studied appear to be related to deep circulation of meteoric water in areas of 'normal' regional heat flow rather than to shallow-crustal heat sources. Discharge temperatures of thermal springs in the region range from slightly above mean annual air temperature (8?-12?C at most places) to boiling or slightly hotter. Geochemical data indicate that, in the major systems, subsurface temperatures at which thermal waters equilibrate with reservoir rocks range from 150? to more than 200?C. These data also indicate that the major systems are of the hot-water type rather than the vapor-dominated type. Depths of thermal-water circulation probably range from 2 to 6 kilometres in areas of 'normal' regional heat flow (~2 heatflow units) and from 1 to 3 kilometres in areas of high heat flow (~3-4 heat-flow units) such as near Battle Mountain. Most of the heat is discharged from the hydrothermal systems studied by (1) conduction through near-surface materials heated as a consequence of thermal-water convection, (2) convection as springflow, and (3) convection as steam discharge from spring pools, vents, fumaroles, and cracks. The mate of heat discharge by radiation from warm ground and by convection as lateral ground-water outflow is believed to be small in most systems and is not estimated. Estimates of net heat discharge from the systems studied range from about 0.8 x 106 calories per second at Buffalo Valley Hot Springs to about 14 x 106 calories per second at Stillwater. These estimates represent the approximate magnitude of the excess heat discharge from the thermal areas that results from the upward convection of hot water from deep sources. Water discharges from the hydrothermal systems by springflow, evapotranspiration, steam discharge, and lateral ground-water outflow. Estimated discharges range from about 0.2 x 106cubic metres per year from the Buffalo Valley Hot Springs system to about 3 x 106 cubic metres per year from the Stillwater system. In most of the hydrothermal systems studied and, by inference, in other similar systems in northern and central Nevada, the scale for potential .commercial development for production of electricity or for other uses may be constrained by the

Open-File Report

Geohydrology and ground-water resources of Philadelphia, Pennsylvania

The aquifers underlying the 134.6-square-mile city of Philadelphia are divided by the Fall Line into the unconsolidated aquifers (chiefly sand and gravel) of the Coastal Plain and the consolidated-rock aquifers (chiefly schist of the Wissahickon Formation) of the Piedmont. Ground water is present under confined and unconfined conditions. The principal units of the confined-aquifer system are the lower and middle sands of the Potomac-Raritan-Magothy aquifer system. The lower sand unit is the most productive aquifer in Philadelphia. The median yield of wells screened in the lower sand unit is 275 gal/min (gallons per minute), and yields of some wells are as high as 1,350 gal/min. The median specific capacity is 16 (gal/min)/ft (gallons per minute per foot of drawdown). The principal units of the unconsolidated unconfined-aquifer system are the upper sand unit of the Potomac-Raritan-Magothy aquifer system and the informally named Trenton gravel. The median yield of wells tapping these two undifferentiated units is 90 gal/min, and yields of some wells are as high as 1,370 gal/min. The median specific capacity is 12 (gal/min)/ft. The consolidated unconfined-aquifer system consists mainly of the Wissahickon Formation. The median yield of nondomestic wells that tap the Wissahickon Formation is 45 gal/min, and yields are as high as 350 gal/min. The median specific capacity is 0.5 (gal/min)/ft. Urbanization has considerably modified the hydrologic cycle in Philadelphia. Impervious surfaces have reduced recharge areas and evapotranspiration and have increased direct runoff. Leakage from the water-distribution system, which is supplied from the Delaware and Schuylkill Rivers, was about 60 to 72 Mgal/d (million gallons per day) in 1980. Groundwater infiltration to sewers is estimated to be as much as 135 Mgal/d when the water table is high. The potentiometric surface of the lower sand unit has been lowered substantially by pumping. By 1954, cones of depression were more than 50 ft (feet) below sea level at the U.S. Naval Base and more than 70 ft below sea level along the Delaware River northeast of the naval base. As a result of withdrawals, declining heads in the lower sand unit caused water to flow downward from the overlying unconsolidated deposits and the water table to decline below sea level along the Delaware River. Beginning in the mid1960's, ground-water withdrawals from the lower sand unit decreased, and, by 1979, water levels had risen 25 ft at the U.S. Naval Base and 45 ft farther north along the Delaware River. As of 1985, water levels in the lower sand unit were controlled largely by pumping in nearby parts of New Jersey. Urbanization also has caused substantial degradation of the quality of ground water in Philadelphia. By 1945, the quality of water in the unconfined aquifer system began to deteriorate as contaminants present at the land surface migrated down- ward. Withdrawal of water from the deeper confined aquifers caused a head decline that resulted in downward movement of contaminated water from the overlying unconfined aquifer system. Consequently, water in the confined aquifers deteriorated progressively in chemical quality so it resembles water in the unconfined aquifer system. The concentration of dissolved solids in water samples collected during 1979-80 ranged from 90 to 4,480 mg/L (milligrams per liter). The average concentration of 778 mg/L was 45 percent higher than that of samples collected during 1945-58. Water from the unconfined unconsolidated aquifers generally had the highest dissolved-solids concentration. The concentration of dissolved iron in water samples collected during 197980 ranged from 0 to 220 mg/L and exceeded 0.30 mg/L in 71 percent of the samples. The average concentration of 17 mg/L was nearly 30 percent higher than that of samples collected during 1945-58. Many wells have been abandoned because of elevated iron concentrations. The concentration of dissolved manganese in water

Water Supply Paper

Assessing native and introduced fish predation on migrating juvenile salmon in Priest Rapids and Wanapum Reservoirs, Columbia River, Washington, 2009--11

Hydroelectric development on the mainstem Columbia River has created a series of impoundments that promote the production of native and non-native piscivores. Reducing the effects of fish predation on migrating juvenile salmonids has been a major component of mitigating the effects of hydroelectric development in the Columbia River basin. Extensive research examining juvenile salmon predation has been conducted in the lower Columbia River. Fewer studies of predation have been done in the Columbia River upstream of its confluence with the Snake River; the most comprehensive predation study being from the early 1990s. The Public Utility District No. 2 of Grant County, Washington initiated a northern pikeminnow removal program in 1995 in an attempt to reduce predation on juvenile salmonids. However, there has been no assessment of the relative predation within the Priest Rapids Project since the removal program began. Further, there is concern about the effects of piscivores other than northern pikeminnow ( Ptychocheilus oregonensis ), such as channel catfish ( Ictalurus punctatus ), smallmouth bass ( Micropterus dolomieu ), and walleye ( Sander vitreus , formerly Stizostedion vitreum ). The Public Utility District No. 2 of Grant County, Washington and the Priest Rapids Coordinating Committee requested that the U.S. Geological Survey, in collaboration with the Washington Department of Fish and Wildlife, assist them in evaluating the effects of native and introduced predatory fish on migrating juvenile salmon. From 2009 to 2010, we conducted sampling in the 103 kilometers (64 river miles) of the Columbia River from the tailrace of Rock Island Dam downstream to the tailrace of Priest Rapids Dam. To assess predation, we used electrofishing to collect northern pikeminnow, smallmouth bass, and walleye to analyze their diets during 2009 and 2010. In 2009, we used methods to allow comparisons to a previous study conducted in 1993. During 2009, we also used an alternate sampling strategy using habitat data and geographic information system software to select sites and allocate samples. In 2010, we used the data collected during 2009 to further refine our sampling design, with the intent of using the data collected during 2010 to formulate a design strategy for implementation during 2011. Based on the results of 2011, we would then propose a strategy for future studies. However, during 2011, our efforts were redirected to specifically address factors that may be affecting steelhead trout survival in the Priest Rapids Reservoir, Columbia River. We used the catch and diet data collected in 2009 and 2010 to estimate relative abundance, consumption, and predation indices for northern pikeminnow and smallmouth bass. Despite extensive sampling in the study area in 2009 and 2010, very few channel catfish and walleye were captured. The mean total lengths of northern pikeminnow were much lower than those observed in 1993; suggesting that efforts to remove northern pikeminnow in the study area may be shifting the population towards smaller fish. The northern pikeminnow predation index values were lower in 2009 than in the 1993 study. The reduced predation levels observed may be due to the prevalence of smaller pikeminnow in our catches than in catches reported in 1993. Predation by smallmouth bass was lower in 2009 than in 2010, and generally was greater than predation for northern pikeminnow. Predation for northern pikeminnow was concentrated in the tailrace areas of Priest Rapids, Wanapum, and Rock Island Dams; predation for smallmouth bass was concentrated in the forebay and mid-reservoir sections of the study area. Our results indicate areas where control measures for smallmouth bass could be concentrated to reduce predation in the Priest Rapids Project.

Washington

Southern Salish Sea Habitat Map Series: Admiralty Inlet

In 2010 the Environmental Protection Agency, Region 10 initiated the Puget Sound Scientific Studies and Technical Investigations Assistance Program, designed to support research in support of implementing the Puget Sound Action Agenda. The Action Agenda was created in response to Puget Sound having been designated as one of 28 estuaries of national significance under section 320 of the U.S. Clean Water Act, and its overall goal is to restore the Puget Sound Estuary's environment by 2020. The Southern Salish Sea Mapping Project was funded by the Assistance Program request for proposals process, which also supports a large number of coastal-zone- and ocean-management issues. The issues include the recommendations of the Marine Protected Areas Work Group to the Washington State Legislature (Van Cleve and others, 2009), which endorses a Puget Sound and coast-wide marine conservation needs assessment, gap analysis of existing Marine Protected Areas (MPA) and recommendations for action. This publication is the first of four U.S. Geological Survey Scientific Investigation Maps that make up the Southern Salish Sea Mapping Project. The remaining three map blocks to be published in the future, located south of Admiralty Inlet, are shown in figure 1. Puget Sound is a deep, fjord-type estuary covering an area of 2,330 km 2 in the Pacific Northwest region of the United States (fig. 1). It is connected to the ocean by the Strait of Juan de Fuca, a turbulent passage approximately 160 km in length and 22 km wide at its west end, expanding to over 40 km wide at its east end (Thomson, 1994). During the Pleistocene, the area was occupied several times by lobes of continental ice, resulting in a complex basin-fill of glacial and interglacial deposits that are locally as thick as 1100 m (Johnson and others, 2001). The last glaciation, called the Fraser glaciation, began after 28,800&plusmn;740 14 C yr B.P. when ice started a slow expansion (Clague, 1981). At peak advance the westward Juan de Fuca lobe reached the edge of the continental shelf through the Juan de Fuca Strait shortly before 14,460&plusmn;200 14 C yr B.P. (Herzer and Bornhold, 1982). The southward Puget lobe advanced to its terminal position in Puget Sound by around 14,150 14 C yr B.P. (Porter and Swanson, 1998). Ice retreated from its maximum to northern Whidbey Island by 13,650&plusmn;350 14 C yr B.P. (Dethier and others, 1995). Retreating glaciers resulted in a thick sequence of ice-contact, glacial-marine sediment, and early post-glacial sediments (Linden and Schurrer, 1988). These deposits have experienced the effects of a marine transgression followed by regression, resulting in a sea-level several tens of meters lower than the present day (Linden and Schurrer, 1988). A second transgression brought sea level to about the present level by around 5,470&plusmn;120 14 C yr B.P. (Clague and others, 1982) establishing the present oceanographic and geologic environment Puget Sound is separated into four interconnected basins; Whidbey, Central (Main), Hood Canal, and South (Thomson, 1994). The Whidbey, Central, and Hood Canal basins are the three main branches of the Puget Sound estuary and are separated from the Strait of Juan de Fuca by a double sill at Admiralty Inlet. The Admiralty Inlet map area includes the Inlet and a portion of the Whidbey Basin (fig. 1). The shallower South Basin is separated by a sill at Tacoma Narrows and is highly branched with numerous finger inlets. Flow within Puget Sound is dominated by tidal currents of as much as 1 m/s at Admiralty Inlet, reducing to approximately 0.5 m/s in the Central Basin (Lavelle and others, 1988). The lack of silt and clay-sized sediments in the Admiralty Inlet map area is likely a result of the strong currents (see Ground-Truth Studies for the Admiralty Inlet Map Area, sheet 3). The subtidal component of flow reaches approximately 0.1 m/s and is driven by density gradients arising from the contrast in salty ocean water at the entrance and freshwater inputs from stream flow (Lavelle and others, 1988). The total freshwater input to Puget Sound is approximately 3.4 x 10 6 m 3 /day, primarily from the Skagit River (Cannon, 1983). The subtidal circulation mostly consists of a two-layered flow in the basins with fresher water exiting at the surface and saltier water entering at depth (Ebbesmeyer and Cannon, 2001). In general, surface waters flow north and deeper waters flow south; variations arise from wind effects that can drive a surface current in the same direction as the wind, and a baroclinic response in the lower layer to about 100-m depth (Matsuura and Cannon, 1997). Oceanographic properties are influenced by temporal forcing parameters such as reduced stream flow during the 2000-01 drought that increased surface salinity and decreased differences between surface and bottom waters (Newton and others, 2003). On offshore seismic-reflection profiles, Pleistocene strata (excluding latest Pleistocene glacial and post-glacial deposits) form a distinct seismic unit, bounded below by pre-Tertiary or Tertiary basement and above by typically flat-lying latest Pleistocene to Holocene deposits that fill in erosional or depositional relief (Johnson and others, 2001). Cores from central Puget Sound have accumulation rates that range from 85 to 1200 mg/cm 2 /yr, or 0.12 to 2.4 cm/yr; the highest accumulation rates are near the southern end of central Puget Sound (Carpenter and others, 1985). Carpenter and others (1985) un-weighted arithmetic mean of accumulation rates for central Puget Sound deeper stations is 480&plusmn;340 (&plusmn; one standard deviation) mg/cm 2 /yr. Lavelle and others (1985) also found rates as high as 1200 mg/cm 2 /yr over the past approximately 70 years in cores in the Central Basin off of and north and south of Elliott Bay. Puget Sound basin rates are comparable to rates in midshelf silt deposits on the Washington coast north of the Columbia River (Nittrouer and others, 1979). The deep subtidal (in other words, below SCUBA depths) habitats of Puget Sound are relatively poorly known. A few subtidal surveys exist for several habitat types from the 1960s and 1970s (reviewed in Dethier, 1990), using grab and box core data. The Dethier (1990) review divides habitat up into Coast and Marine Ecological Classification Standard (CMECS) substrate, water column energy, and depth zones but does not attempt to map these habitats, rather it is an inventory of habitats found in the area and the flora and fauna associated with each habitat. The approach of the Southern Salish Sea Mapping project is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data (the undersea equivalent of satellite remote-sensing data in terrestrial mapping), acoustic backscatter, seafloor video, seafloor photography, and bottom-sediment sampling data. This approach is based in part on methods presented and data collection and product needs identified at the Washington State Seafloor Mapping Workshop (Washington State Seafloor Mapping Workshop Steering Committee, 2008), attended by coastal and marine managers and scientists. The map products display seafloor geomorphology and substrate, and identify potential marine benthic habitats. It is emphasized that the more interpretive habitat and geology maps rely on the integration of multiple, new high-resolution datasets and that mapping at small scales would not be possible without such data. Oceanographic current and wave data is not included in this analysis, however, the accompanying geographic information system (GIS) data set is designed and intended to be combined with oceanographic and biologic data sets assembled by others in the future and some of the GIS data has already been incorporated in the unpublished Nature Conservancy Benthic Habitats of Puget Sound database. This publication includes four map sheets, explanatory text, and a descriptive pamphlet. Each map sheet is published as a portable document format (PDF) file. ESRI ArcGIS compatible geotiffs (for example, bathymetry) and shapefiles (for example video observation points) will be available for download in the data catalog associated with this publication (Cochrane, 2015). An ArcGIS Project File with the symbology used to generate the map sheets is also provided. For those who do not own the full suite of ESRI GIS and mapping software, the data can be read using ESRI ArcReader, a free viewer that is available at http://www.esri.com/software/arcgis/arcreader/index.html .

Washington

Hydrology and water quality of Shell Lake, Washburn County, Wisconsin, with special emphasis on the effects of diversion and changes in water level on the water quality of a shallow terminal lake

Shell Lake is a relatively shallow terminal lake (tributaries but no outlets) in northwestern Wisconsin that has experienced approximately 10 feet (ft) of water-level fluctuation over more than 70 years of record and extensive flooding of nearshore areas starting in the early 2000s. The City of Shell Lake (City) received a permit from the Wisconsin Department of Natural Resources in 2002 to divert water from the lake to a nearby river in order to lower water levels and reduce flooding. Previous studies suggested that water-level fluctuations were driven by long-term cycles in precipitation, evaporation, and runoff, although questions about the lake&rsquo;s connection with the groundwater system remained. The permit required that the City evaluate assumptions about lake/groundwater interactions made in previous studies and evaluate the effects of the water diversion on water levels in Shell Lake and other nearby lakes. Therefore, a cooperative study between the City and U.S. Geological Survey (USGS) was initiated to improve the understanding of the hydrogeology of the area and evaluate potential effects of the diversion on water levels in Shell Lake, the surrounding groundwater system, and nearby lakes. Concerns over deteriorating water quality in the lake, possibly associated with changes in water level, prompted an additional cooperative project between the City and the USGS to evaluate efeffects of changes in nutrient loading associated with changes in water levels on the water quality of Shell Lake. Numerical models were used to evaluate how the hydrology and water quality responded to diversion of water from the lake and historical changes in the watershed. The groundwater-flow model MODFLOW was used to simulate groundwater movement in the area around Shell Lake, including groundwater/surface-water interactions. Simulated results from the MODFLOW model indicate that groundwater flows generally northward in the area around Shell Lake, with flow locally converging toward the lake. Total groundwater inflow to Shell Lake is small (approximately 5 percent of the water budget) compared with water entering the lake from precipitation (83 percent) and surface-water runoff (13 percent). The MODFLOW model also was used to simulate average annual hydrologic conditions from 1949 to 2009, including effects of the removal of 3 billion gallons of water during 2003&ndash;5. The maximum decline in simulated average annual water levels for Shell Lake due to the diversion alone was 3.3 ft at the end of the diversion process in 2005. Model simulations also indicate that although water level continued to decline through 2009 in response to local weather patterns (local drought), the effects of the diversion decreased after the diversion ceased; that is, after 4 years of recovery (2006&ndash;9), drawdown attributable to the diversion alone decreased by about 0.6 ft because of increased groundwater inflow and decreased lake-water outflow to groundwater caused by the artificially lower lake level. A delayed response in drawdown of less than 0.5 ft was transmitted through the groundwater-flow system to upgradient lakes. This relatively small effect on upgradient lakes is attributed in part to extensive layers of shallow clay that limit lake/groundwater interaction in the area. Data collected in the lake indicated that Shell Lake is polymictic (characterized by frequent deep mixing) and that its productivity is limited by the amount of phosphorus in the lake. The lake was typically classified as oligotrophic-mesotrophic in June, mesotrophic in July, and mesotrophic-eutrophic in August. In polymictic lakes like Shell Lake, phosphorus released from the sediments is not trapped near the bottom of the lake but is intermittently released to the shallow water, resulting in deteriorating water quality as summer progresses. Because the productivity of Shell Lake is limited by phosphorus, the sources of phosphorus to the lake were quantified, and the response in water quality to changes in phosphorus inputs were evaluated by means of eutrophication models. During 2009, the total input of phosphorus to Shell Lake was 1,730 pounds (lb), of which 1,320 lb came from external sources (76 percent) and 414 lb came from internal loading from sediments in the lake (24 percent). The largest external source was from surface-water runoff, which delivered about 52 percent of the total phosphorus load compared with about 13 percent of the water input. The second largest source was from precipitation (wetfall and dryfall), which delivered 19 percent of the load compared to about 83 percent of the water input. Contributions from septic systems and groundwater accounted for about 3 and 2 percent, respectively. Increased runoff raises water levels in the lake but does not necessarily increase phosphorus loading because phosphorus concentrations in the tributaries decline during increased flow, possibly because of shorter retention times in upstream wetlands. Phosphorus loading to the lake in 2009 represented what occurred after a series of dry years; therefore, this information was combined with data from 2011, a wet year, to estimate phosphorus loading during a range of hydrologic conditions by estimating loading from each component of the phosphorus budget for each year from 1949 to 2011. Comparisons of historical water-quality records with historical water levels and applications of a hydrodynamic model (Dynamic Lake Model, DLM) and empirical eutrophication models were used to understand how changes in water level and the coinciding changes in phosphorus loading affect the water quality of Shell Lake. DLM simulations indicate that large changes in water level (approximately 10 ft) affect the persistence of stratification in the lake. During periods with low water levels, the lake is a well-mixed, polymictic system, with water quality degrading slightly as summer progresses. During periods with high water levels, the lake is more stratified, and phosphorus from internal loading is trapped in the hypolimnion and released later in summer, which results in more extreme seasonality in water quality and better clarity in early summer. Results of eutrophication model simulations using a range in external phosphorus inputs illustrate how water quality in Shell Lake (phosphorus and chlorophyll a concentrations and Secchi depths) responds to changes in external phosphorus loading. Results indicate that a 50-percent reduction in external loading from that measured in 2009 would be required to change phosphorus concentrations from 0.018 milligram per liter (mg/L) (measured in 2009) to 0.012 mg/L (estimated for the mid-1800s from analysis of diatoms in sediment cores). Such reductions in phosphorus loading cannot be accomplished by targeting septic systems or internal loading alone because septic systems contribute only about 3 percent of the phosphorus input to the lake, and internal loading from the sediments of Shell Lake contributes only about 25 percent of phosphorus input. Complete elimination of phosphorus from septic systems and internal loading would decrease the phosphorus concentrations in the lake by 0.003&ndash;0.004 mg/L. Therefore, reducing phosphorus concentration in the lake more than by 0.004 mg/L requires decreasing phosphorus loading from surface-water contributions, primarily runoff to the lake. Reconstructed changes in water quality from 1860 to 2010, based on changes in the diatom communities archived in the sediments and eutrophication model simulations, suggest that anthropogenic changes in the watershed (sawmill construction in 1881; the establishment of the village of Shell Lake; and land-use changes in the 1920s, including increased agriculture) had a much larger effect on water quality than the natural changes associated with fluctuations in water level. Although the effects of natural changes in water level on water quality appear to be small, changes in water level do have a modest effect on water quality, primarily manifested as small improvements during higher water levels. Fluctuations in water level, however, have a larger effect on the seasonality of water-quality patterns, with better water quality, especially increased Secchi depths, in early summer during years with high water levels.

Wisconsin

Petroleum system modeling capabilities for use in oil and gas resource assessments

Summary: Petroleum resource assessments are among the most highly visible and frequently cited scientific products of the U.S. Geological Survey. The assessments integrate diverse and extensive information on the geologic, geochemical, and petroleum production histories of provinces and regions of the United States and the World. Petroleum systems modeling incorporates these geoscience data in ways that strengthen the assessment process and results are presented visually and numerically. The purpose of this report is to outline the requirements, advantages, and limitations of one-dimensional (1-D), two-dimensional (2-D), and three-dimensional (3-D) petroleum systems modeling that can be applied to the assessment of oil and gas resources. Primary focus is on the application of the Integrated Exploration Systems (IES) PetroMod? software because of familiarity with that program as well as the emphasis by the USGS Energy Program on standardizing to one modeling application. The Western Canada Sedimentary Basin (WCSB) is used to demonstrate the use of the PetroMod? software. Petroleum systems modeling quantitatively extends the 'total petroleum systems' (TPS) concept (Magoon and Dow, 1994; Magoon and Schmoker, 2000) that is employed in USGS resource assessments. Modeling allows integration of state-of-the-art analysis techniques, and provides the means to test and refine understanding of oil and gas generation, migration, and accumulation. Results of modeling are presented visually, numerically, and statistically, which enhances interpretation of the processes that affect TPSs through time. Modeling also provides a framework for the input and processing of many kinds of data essential in resource assessment, including (1) petroleum system elements such as reservoir, seal, and source rock intervals; (2) timing of depositional, hiatus, and erosional events and their influences on petroleum systems; (3) incorporation of vertical and lateral distribution and lithologies of strata that compose the petroleum systems; and (4) calculations of pressure-volume-temperature (PVT) histories. As digital data on petroleum systems continue to expand, the models can integrate these data into USGS resource assessments by building and displaying, through time, areas of petroleum generation, migration pathways, accumulations, and relative contributions of source rocks to the hydrocarbon components. IES PetroMod? 1-D, 2-D, and 3-D models are integrated such that each uses the same variables for petroleum systems modeling. 1-D burial history models are point locations, mainly wells. Maps and cross-sections model geologic information in two dimensions and can incorporate direct input of 2-D seismic data and interpretations using various formats. Both 1-D and 2-D models use data essential for assessments and, following data compilation, they can be completed in hours and retested in minutes. Such models should be built early in the geologic assessment process, inasmuch as they incorporate the petroleum system elements of reservoir, source, and seal rock intervals with associated lithologies and depositional and erosional ages. The models can be used to delineate the petroleum systems. A number of 1-D and 2-D models can be constructed across a geologic province and used by the assessment geologists as a 3-D framework of processes that control petroleum generation, migration, and accumulation. The primary limitation of these models is that they only represent generation, migration, and accumulation in two dimensions. 3-D models are generally built at reservoir to basin scales. They provide a much more detailed and realistic representation of petroleum systems than 1-D or 2-D models because they portray more fully the temporal and physical relations among (1) burial history; (2) lithologies and associated changes through burial in porosity, permeability, and compaction; (3) hydrodynamic effects; and (4) other parameters that influence petroleum gen

Open-File Report

Assessment and management of dead-wood habitat

The Bureau of Land Management (BLM) is in the process of revising its resource management plans for six districts in western and southern Oregon as the result of the settlement of a lawsuit brought by the American Forest Resource Council. A range of management alternatives is being considered and evaluated including at least one that will minimize reserves on O&C lands. In order to develop the bases for evaluating management alternatives, the agency needs to derive a reasonable range of objectives for key issues and resources. Dead-wood habitat for wildlife has been identified as a key resource for which decision-making tools and techniques need to be refined and clarified. Under the Northwest Forest Plan, reserves were to play an important role in providing habitat for species associated with dead wood (U.S. Department of Agriculture Forest Service and U.S. Department of the Interior Bureau of Land Management, 1994). Thus, the BLM needs to: 1) address the question of how dead wood will be provided if reserves are not included as a management strategy in the revised Resource Management Plan, and 2) be able to evaluate the effects of alternative land management approaches. Dead wood has become an increasingly important conservation issue in managed forests, as awareness of its function in providing wildlife habitat and in basic ecological processes has dramatically increased over the last several decades (Laudenslayer et al., 2002). A major concern of forest managers is providing dead wood habitat for terrestrial wildlife. Wildlife in Pacific Northwest forests have evolved with disturbances that create large amounts of dead wood; so, it is not surprising that many species are closely associated with standing (snags) or down, dead wood. In general, the occurrence or abundance of one-quarter to one-third of forest-dwelling vertebrate wildlife species, is strongly associated with availability of suitable dead-wood habitat (Bunnell et al., 1999; Rose et al., 2001). In Oregon and Washington, approximately 150 species of wildlife are reported to use dead wood in forests (O’Neil et al., 2001). Forty-seven sensitive and special-status species are associated with dead wood (Appendix A). These are key species for management consideration because concern over small or declining populations is often related to loss of suitable dead-wood habitat (Marshall et al., 1996). Primary excavators (woodpeckers) also are often the focus of dead-wood management, because they perform keystone functions in forest ecosystems by creating cavities for secondary cavity-nesters (Martin and Eadie, 1999; Aubry and Raley, 2002). A diverse guild of secondary cavity-users (including swallows, bluebirds, several species of ducks and owls, ash-throated flycatcher, flying squirrel, bats, and many other species) is unable to excavate dead wood, and therefore relies on cavities created by woodpeckers for nesting sites. Suitable nest cavities are essential for reproduction, and their availability limits population size (Newton, 1994). Thus, populations of secondary cavity-nesters are tightly linked to the habitat requirements of primary excavators. Although managers often focus on decaying wood as habitat for wildlife, the integral role dead wood plays in ecological processes is an equally important consideration for management. Rose et al. (2001) provide a thorough review of the ecological functions of dead wood in Pacific Northwest forests, briefly summarized here. Decaying wood functions in: soil development and productivity, nutrient cycling, nitrogen fixation, and carbon storage. From ridge tops, to headwater streams, to estuaries and coastal marine ecosystems, decaying wood is fundamental to diverse terrestrial and aquatic food webs. Wildlife species that use dead wood for cover or feeding are linked to these ecosystem processes through a broad array of functional roles, including facilitation of decay and trophic interactions with other organisms (Marcot, 2002; Marcot, 2003). For example, by puncturing bark and fragmenting sapwood, woodpeckers create sites favorable for wood-decaying organisms (Farris et al., 2004), which in turn create habitat for other species and facilitate nutrient cycling. Small mammals that use down wood for cover function in the dispersal of plant seeds and fungal spores (Carey et al., 1999). Resident cavitynesting birds may regulate insect populations by preying on overwintering arthropods (Jackson, 1979; Kroll and Fleet, 1979). These examples illustrate how dead wood not only directly provides habitat for a large number of wildlife species, but also forms the foundation of functional webs that critically influence forest ecosystems (Marcot, 2002; Marcot, 2003). The important and far-reaching implications of management of decaying wood highlight the need for conservation of dead-wood resources in managed forests. Consideration of the key ecological functions of species associated with dead wood can help guide management of dead wood in a framework consistent with the paradigm of ecosystem management (Marcot and Vander Heyden, 2001; Marcot, 2002.) As more information is revealed about the ecological and habitat values of decaying wood, concern has increased over a reduction in the current amounts of dead wood relative to historic levels (Ohmann and Waddell, 2002). Past management practices have tended to severely reduce amounts of dead wood throughout all stages of forest development (Hansen et al., 1991). The large amounts of legacy wood that characterize young post-disturbance forests are not realized in managed stands, because most of the wood volume is removed at harvest for economic and safety reasons. Mid-rotation thinning is used to “salvage” some mortality that might otherwise occur due to suppression, so fewer snags are recruited in mid-seral stages. Harvest rotations of 80 years or less truncate tree size in managed stands, and thus limit the production of large-diameter wood. As a consequence of these practices, dead wood has been reduced by as much as 90% after two rotations of managed Douglas-fir (Rose et al., 2001). Large legacy deadwood is becoming a scarce, critical habitat that will take decades to centuries to replace. Furthermore, management continues to have important direct and indirect effects on the amount and distribution of dead wood in forests. Current guidelines for managing dead wood may be inadequate to maintain habitat for all associated species because they largely focus on a single use of dead wood (nesting habitat) by a small suite of species (cavity-nesting birds), and may under represent the sizes and amounts of dead wood used by many wildlife species (Rose et al., 2001, Wilhere, 2003).

Open-File Report

Energy and minerals division tight oil and gas committee: Activities and commodity report for 2019-2020

In 2019, total daily tight oil and gas production increased in the United States month over month, with annualized growth of 14% for oil and 12% for gas. Those gains leveled off in the first quarter of 2020 due to aggressive price competition and increases in international production. Then came the pandemic with a substantially larger dose of economic turmoil, driving down demand due in part to shelter in place orders and safety concerns around travel. Between March and May, tight oil and gas production dropped by nearly 2 million bpd and almost 5 Bcf/day before beginning to recover. Production has continued to increase for the most part through the second half of 2020, but drilling remains subdued throughout most of the U.S. and uncertainty around long term demand along with the current price environment and general state of the economy has contributed to layoffs throughout the industry. Some shale-gas production has declined recently, but a few areas have seen expansion due to construction of LNG facilities along the East Coast of the U.S. (e.g., the Haynesville Formation). Current U.S. shale-gas production is still higher now than in 2019, with daily production of almost 71 Bcf as of October 2020 driven in large part by increased production from the Marcellus Shale in the Appalachian Basin and shales within the Permian Basin. Shale liquids production is down by around a million bpd to approximately 7.1 million (September 2020; U.S. EIA) from pre-pandemic production levels at the end of 2019 and beginning of 2020. Tight oil production remains dominated by plays in the Permian Basin as well as the Bakken and Eagle Ford Formations. On the development and production front, new enhanced oil recovery approaches for tight shale reservoirs are being more widely implemented. Natural gas or CO2 injection is currently being utilized in the Bakken Formation, Eagle Ford Formation, Anadarko Basin, and the Permian Basin to optimize injection sequences and boost recovery. Refracturing of existing wells to reduce drilling costs, improve production, and prolong well productive life has also begun to occur more widely in developed plays. International interest in exploiting hydrocarbons from unconventional reservoirs continues to develop, with active exploration projects on most continents. Europe remains relatively underexplored as compared to North America, although a total of 141 exploration and appraisal wells with a possible shale-gas exploration component have been spudded, including horizontal legs from vertical wells. Shale exploration has made a breakthrough in China with shale gas output in 2019 of 10 billion cubic meters (35.3 Bcf), 60% of which was produced from Sinopec’s Fuling Shale Gas field. Lacustrine shale oil exploration has also been successful in the Sichuan and Ordos Basins in central China, Junggar and Tarim Basins in northwest China, and Songliao Basin in north China, and Bohai Bay Basins in northeast China as of 2018. South America’s potential as an unconventional shale gas and oil province is mainly in Argentina and Brazil, where the production from Neuquen Basin’s tight shale of the Vaca Muerta Formation has been steadily increasing since 2016, but only 4% of the shale resource has been developed thus far. According to International Energy Agency’s report in 2013, Brazil holds the 9th largest unconventional gas reserves. Brazil has shale oil and gas potential in the Parana, Solimoes and Amazon Basins and is actively producing from the oil shale unit of the Irati Formation. In 2019, the Brazil energy ministry launched REATE 2020 to boost onshore investments that include the expectation of drilling an experimental unconventional well in the northeast region. For this inaugural report, the new AAPG EMD Tight Oil and Gas Committee (TO&G; formerly the Shale Gas & Liquids and Tight Gas Sands committees) has developed new commodity report requirements for contributors. This includes shorter annual reports focused on new developments, play concepts, along with the typical updates on production and new drilling in the play areas they cover. We are also asking contributors to collect background geologic and production related information into a document that summarizes important features of the plays they cover that will be stored on the TO&G webpage along with our commodity reports. TO&G is currently working to expand the number of contributors to cover more play areas and replace committee and advisory board members that have recently stepped down. Changes to committee leadership occurred in October as recent chairs transition to EMD elected positions.

Activities and Commodities Report

Summary appraisals of the nation's ground-water resources – Souris-Red-Rainy region

A broad-perspective analysis of the ground-water resources and present and possible future water development and management in the Souris-Red-Rainy Region is presented. The region includes the basins of the Souris River within Montana and North Dakota; the Red River of the North in South Dakota, North Dakota, and Minnesota; and the Rainy River within Minnesota. The region includes 59,645 square miles, mostly in North Dakota and Minnesota. The terrain is relatively flat, but ranges in altitude from 2,541 to 750 feet. Annual average precipitation ranges from 14 inches in the west to 28 inches in the east and about 75 percent of it is rain. The mean annual snowfall ranges from 32 inches in the west to 64 inches in the east. Temperatures range from -55&deg; to 118&deg; F (-48.3&deg; to 47.8&deg; C). Irrigation is needed at least part of the time to assure crop production, particularly in the western part of the region. Sand and gravel deposits in the drift form the most important freshwater aquifers. Other aquifers are found in at least parts of the region in the Precambrian, Paleozoic, Cretaceous, and Tertiary rocks. The potentiometric surface in the bedrock aquifers generally decreases in altitude toward the Red River of the North, indicating that the general direction of ground-water movement is toward the river. Ground water with less than 3,000 milligrams per liter dissolved solids is available throughout the region. Ground water with less than 1,000 milligrams per liter occurs in most of the region east of the Red River of the North and in most of the shallow aquifers west of the river. The total volume of water available from storage having less than 3,000 milligrams per liter dissolved solids is estimated to be 5x10 8 acre-feet. In addition to the fresh and slightly saline water, the region has abundant highly mineralized water that can be considered as a resource. Yields of wells in individual bedrock aquifers are generally less than 100 gallons per minute but locally yields may be as much as 500 gallons per minute and more. Yields in drift aquifers are frequently less than 100 gallons per minute but range from 5 to 1,000 gallons per minute. In a few places outwash yields more than 1,000 gallons per minute. Ground water is the sole or a primary source of water supply in much of the region, including supplies for irrigation, domestic and livestock, municipal, and industrial needs. Reportedly, the potential irrigation development is 1,550,000 acres, as compared with 50,200 acres in 1975. Both ground- and surface-water supplies would be required to meet these demands. Rural domestic and livestock water supplies are derived almost entirely from ground-water sources. Smaller communities and towns generally rely on ground water, and the cities and industries use ground water, surface water, or both. The municipalities using surface water generally depend upon reservoir storage. Water quality rather than quantity is the greater water-supply problem for many communities in the region. Increased demands on both ground-water and surface-water supplies likely will be made in the future. Storage of surface water in the ground-water reservoirs during times of surplus for withdrawal during times of scarcity would aid in meeting these demands. The surplus (flood) water is of better chemical quality than underlying ground water in parts of the western half of the region. Freshwater could be stored in saline- or freshwater aquifers, and pumped out later, as needed. Thus, the ground-water reservoirs have a definite present and potential role in water management. To understand the hydrologic system for management purposes there is a need to determine more adequately the geologic and hydrologic characteristics of existing aquifers and the location of new aquifers. Also, as pumping and other stresses on any part of the hydrologic system affect other parts of the system, monitoring programs ideally should be started and maintained to detect changes and determine effects of the stresses. Many alternatives are available for managing water in the region. Some of these are operational and others are undergoing research. Adequate hydrologic information is needed to aid in solving problems of water supply, use, and pollution.

Minnesota, Montana, North Dakota, South Dakota

Summary appraisals of the Nation's ground-water resources; Souris-Red-Rainy region

A broad-perspective analysis of the ground-water resources and present and possible future water development and management in the Souris-Red-Rainy Region is presented. The region includes the basins of the Souris River within Montana and North Dakota; the Red River of the North in South Dakota, North Dakota, and Minnesota; and the Rainy River within Minnesota. The region includes 59,645 square miles, mostly in North Dakota and Minnesota. The terrain is relatively flat, but ranges in altitude from 2,541 to 750 feet. Annual average precipitation ranges from 14 inches in the west to 28 inches in the east and about 75 percent of it is rain. The mean annual snow fall ranges from 32 inches in the west to 64 inches in the east. Temperatures range from -55&deg; to 118&deg;F (-48.3&deg; to 47.8&deg;C). Irrigation is needed at least part of the time to assure crop production, particularly in the western part of the region. Sand and gravel deposits in the drift form the most important fresh-water aquifers. Other aquifers are found in at least parts of the region in the Precambrian, Paleozoic, Cretaceous, and Tertiary rocks. The potentiometric surface in the bedrock generally decreases in altitude toward the Red River of the North, indicating that the general direction of ground-water movement is toward the river. Ground-water with less than 3,000 milligrams per liter dissolved solids is available throughout the region. Ground water with less than 1,000 milligrams per liter occurs in most of the region east of the Red River of the North and in most of the shallow aquifers west of the river. The total volume of water available from storage having less than 3,000 milligrams per liter dissolved solids is estimated to be 5 x 10 8 acre-feet. In addition to the fresh and slightly saline water, the region has abundant highly mineralized water that can be considered as a resource. Yields of wells in individual bedrock aquifers are generally less than 100 gallons per minute but locally yields may be as much as 500 gallons per minute and more. Yields in drift aquifers are frequently less than 100 gallons per minute but range from 5 to 1,000 gallons per minute. In a few places outwash yields more than 1,000 gallons per minute. Ground water is the sole or a primary source of water supply in much of the region, including supplies for irrigation, domestic and livestock, municipal, and industrial needs. Reportedly, the potential irrigation development is 1,550,000 acres, as compared with 50,200 acres in 1975. Both ground- and surface-water supplies would be required to meet these demands. Rural domestic and livestock water supplies are derived almost entirely from ground-water sources. Smaller communities and towns generally rely on ground water, and the cities and industries use ground water, surface water, or both. The municipalities using surface water generally depend upon reservoir storage. Water quality rather than quantity is the greater water-supply problem for many communities in the region. Increased demands on both ground-water and surface-water supplies likely will be made in the future. Storage of surface water in the ground-water reservoirs during times of surplus for withdrawal during times of scarcity would aid in meeting these demands. The surplus (flood) water is of better chemical quality than underlying ground water in parts of the western half of the region. Fresh water could be stored in saline- or fresh-water aquifers, and pumped out later, as needed. Thus, the ground-water reservoirs have a definite present and potential role in water management. To understand the hydrologic system for management purposes there is a need to determine more adequately the geologic and hydrologic characteristics of existing aquifers and the location of new aquifers. Also, as pumping and other stresses on any part of the hydrologic system affect other parts of the system, monitoring programs ideally should be started and maintained to detect changes and determine effects of the stresses. Many alternatives are available for managing water in the region. Some of these are operational and others are undergoing research. Adequate hydrologic information is needed to aid in solving problems of water supply, use and pollution.

Minnesota, Montana, North Dakota, South Dakota

Assessment of contaminant trends in plumes and wells and monitoring network optimization at the Badger Army Ammunition Plant, Sauk County, Wisconsin

Soil and groundwater at the Badger Army Ammunition Plant (BAAP), Sauk County, Wisconsin, were affected by several contaminants as a result of production and waste disposal practices common during its operation from 1942 to 1975. Three distinct plumes of contaminated groundwater originate on BAAP property and extend off-site, as identified by previous studies. Routine sampling of groundwater quality from a network of monitoring wells and off-site private wells has been performed since 1990, although the number of wells monitored and the monitoring frequency have varied as the approved monitoring plan was modified. During the period of monitoring from 1990 to 2018, numerous site investigations and remedial actions were conducted to address the sources of contamination, contaminated soils, and groundwater. Concentrations of contaminants reportedly decreased between 2000 and 2012 within all three plumes. Five or six contaminants of concern (COCs) were identified for each of the three plumes. An independent assessment of the contaminant plumes and of the monitoring network was conducted using groundwater-quality data collected from more than 600 wells between 2000 and 2018. In a study conducted by the U.S. Geological Survey (USGS), in cooperation with the Army Environmental Command, a consistent data aggregation and interpolation scheme was applied to derive the likely maximum groundwater plume extents in four 3-year time periods between 2000 and 2018. The plume extent was defined by the Enforcement Standard for each COC and represents the maximum concentration observed in each 3-year time period. The plume boundary analysis shows that the spatial extent of groundwater contamination decreased for most COCs during the study period. Some plume boundaries are not well delineated by the existing monitoring network, particularly the downgradient edge of the Propellant Burning Ground plume. Maps identify the plume boundary in each time period, the sampling well network used to delineate the plume, and wells that were sampled in the 2010–12 period but not sampled in the 2015–18 period. A series of statistical analyses using the Monitoring and Remediation Optimization System, version 3.0, program were applied to the available COC concentration data for two distinct periods, 2000 to 2012 and 2013 to 2018, with the break between periods coinciding with changes to the monitoring network in 2013. Trends in the concentration of COCs in individual wells varied, although generally more wells had decreasing than had increasing concentrations for most COCs in both time periods. The exceptions were ethyl ether in the 2004–12 period and 2,6-dinitrotoluene in the 2013–18 period, for which more wells had an increasing trend. Spatial moment analysis of concentration data from the well network was used to assess the stability of each plume for the COCs. During the 2000–12 period, most of the contaminant plumes for which data were sufficient to complete the analysis were either decreasing or stable in mass and size. The exceptions were carbon tetrachloride (associated solely with the Propellant Burning Ground plume) and 2,4-dinitrotoluene and 2,6-dinitrotoluene (in the Deterrent Burning Ground plume), which showed an increasing trend in mass. No COCs showed an increasing trend in plume mass in the 2013–18 period. Some wells with increasing trends in concentration or with concentrations greater than the enforcement standard are near the tail of a plume, where increased monitoring may be of value to better define future plume boundaries. A spatial optimization analysis covering the 2013–18 period identified six wells that provided information redundant to that from other wells. A temporal optimization analysis identified optimal sampling frequencies for 125 wells. Remedial actions directed at the Propellant Burning Ground plume coincided with a general decrease in plume mass and size, although in specific areas and depths, the plume size for specific contaminants may still be increasing.

Wisconsin