SEARCH · USGS Science

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.

SourceResearch connectionFrom yearThrough year

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

At least 955 records · Page 53Linked to original sources

USGS1987

Hydrology of area 59, northern Great Plains and Rocky Mountain coal provinces, Colorado and Wyoming

Hydrologic information and analysis aid in decisions to lease federally owned coal and to prepare necessary Environmental Assessments and Impact Study reports. This need has become even more critical with the enactment of Public Law 95-87, the "Surface Mining Control and Reclamation Act of 1977." This act requires an appropriate regulatory agency to issue permits, based on the review of permit-application data to assess hydrologic impacts. This report, which partially fulfills this requirement, is one in a series of nationwide coal province reports that present information thematically, through the use of a brief text and accompanying maps, graphs, charts, or other illustrations for single hydrologic topics. The report broadly characterizes the hydrology of Area 59 in north-central Colorado and southeastern Wyoming. The report area, located within the South Platte River basin, covers a 16,000-square-mile area of diverse geology, topography, and climate. This diversity results in contrasting hydrologic characteristics. The South Platte River, the major stream in the area, and most of its tributaries originate in granitic mountains and flow into and through the sedimentary rocks of the Great Plains. Altitudes range from less than 5,000 feet to more than 14,000 feet above sea level. Precipitation in the mountains may exceed 40 inches annually, much of it during the winter, and produces deep snowpacks. Snowmelt during the spring and summer produces most streamflow. Transmountain diversion of water from the streams on the western slope of the mountains also adds to the streamflow. Precipitation in the plains is as little as 10 inches annually. Streams that originate in the plains are ephemeral. Streamflow quality is best in the mountains, where dissolved-solids concentrations are generally small. Concentrations increase in the plains as streams flow through sedimentary basins, and as urbanization and irrigation increase. The quality of some mountain streams is affected by drainage from previous metalmining areas, as indicated by greater trace-element concentrations and smaller pH values. However, the large trace-element concentrations decrease rapidly downstream from the metal-mining areas. Because the climate is semiarid in most of the area, the soils are not adequately leached; therefore, flows in ephemeral streams usually have larger concentrations of dissolved solids than flows in perennial streams. Ground water is available throughout the area; yields range from less than 0.1 gallons per minute in the fractured granite aquifer in the mountains to more than 2,000 gallons per minute in the alluvial aquifer of the South Platte River valley. Major bedrock aquifers in order of decreasing age are the Laramie-Fox Hills, Arapahoe, Denver, and Dawson; these aquifers are used for municipal, domestic, and livestock supplies. Alluvial aquifers supply the high-yield irrigation wells. The best quality ground water is found at the center of the major bedrock aquifers, where dissolved-solids concentrations are less than 200 milligrams per liter. The poorest-quality water is usually found near the edges of these aquifers. Water in the coal-bearing Laramie and Denver Formations is locally affected by coal deposits, causing dissolved-solids concentrations to be relatively large. Only one coal mine is now operating in Area 59, the Coors Energy Company surface coal mine, which produced 100,000 short tons of subbituminous coal from the Upper Cretaceous Laramie Formation in 1982. Past coal-mining operations removed more than 130 million tons of coal and lignite from Area 59,99 percent of which came from underground mines. The largest coal production was in Weld and Boulder Counties, Colorado. Hydrologic problems related to surface mining are erosion, sedimentation, decline in water levels, disruption of aquifers, and degradation of water quality. Because the semiarid mine areas have very little runoff, and the major streams have large buffer and dilution capacities, the effects of mining on surface water is minimal. However, effects on ground water may be much more severe and long-lasting.

Colorado, Wyoming↗

USGS2012

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↗

USGS2015

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±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±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±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±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±340 (± 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↗

USGS2013

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’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–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–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–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↗

USGS2006

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↗

USGS2007

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↗

USGS2020

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↗

USGS1978

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° to 118° F (-48.3° to 47.8° 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↗

USGS1977

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° to 118°F (-48.3° to 47.8°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↗

USGS2021

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↗

USGS2010

Changes in streamflow and the flux of nutrients in the Mississippi-Atchafalaya River Basin, USA, 1980-2007

Nutrients and freshwater delivered by the Mississippi and Atchafalaya Rivers drive algal production in the northern Gulf of Mexico, which eventually results in the widespread occurrence of hypoxic bottom waters along the Louisiana and Texas coast. Researchers have demonstrated a relation between the extent of the hypoxic zone and the magnitude of streamflow, nutrient fluxes, and nutrient concentrations in the Mississippi River, with springtime streamflows and fluxes being the most predictive. In 1999 the U.S. Geological Survey (USGS) estimated the flux of nitrogen, phosphorus, and silica at selected sites in the Mississippi Basin and to the Gulf of Mexico for 1980-1996. These flux estimates provided the baseline information used by the Mississippi River/Gulf of Mexico Watershed Nutrient Task Force to develop an Action Plan for reducing hypoxia in the northern Gulf of Mexico. The primary goal of the Action Plan was to achieve a reduction in the size (areal extent) of the hypoxic zone from an average of approximately 14,000 square kilometers in 1996-2000 to a 5-year moving average of less than 5,000 square kilometers by 2015. Improved statistical models and adjusted maximum likelihood estimation using USGS Load Estimator (LOADEST) software were used to estimate annual and seasonal nutrient fluxes for 1980-2007 at selected sites on the Mississippi River and its tributaries. These data provide a means to evaluate the influence of natural and anthropogenic effects on delivery of water and nutrients to the Gulf of Mexico; to define subbasins that are the most important contributors of nutrients to the gulf; and to investigate the relations among streamflow, nutrient fluxes, and the size and duration of the Gulf of Mexico hypoxic zone. A comparative analysis between the baseline period of 1980-1996 and 5-year moving averages thereafter indicate that the average annual streamflow and fluxes of total nitrogen, nitrate, orthophosphate, and silica to the Gulf of Mexico have decreased. However, the flux of total phosphorus between the baseline period and subsequent 5-year periods has increased. The average spring (April, May, and June) streamflow and fluxes of silica, total nitrogen, nitrate, and orthophosphate to the Gulf of Mexico also decreased, whereas the spring flux of total phosphorus has increased. Similar changes in streamflow and nutrient flux were observed at many sites Buxtonwithin the basin. The inputs of water, total nitrogen, and total phosphorus from the major subbasins of the Mississippi-Atchafalaya River Basin as a percentage of the to-the-gulf totals have increased from the Ohio River Basin, decreased from the Missouri River Basin, and remained relatively unchanged from the Upper Mississippi, Red, and Arkansas River Basins. Changes in streamflow and nutrient fluxes are related, but short-term variations in sources of streamflow and nutrients complicate the interpretation of factors that affect nutrient delivery to the Gulf of Mexico. Parametric time-series models are used to try and separate natural variability in nutrient flux from changes due to other causes. Results indicate that the decrease in annual nutrient fluxes that has occurred between the 1980-1996 baseline period and more recent years can be largely attributed to natural causes (climate and streamflow) and not management actions or other human controlled activities in the Mississippi-Atchafalaya River Basin. The downward trends in total nitrogen, nitrate, ammonium, and orthophosphate that were detected at either the Mississippi River near St. Francisville, La., or the Atchafalaya River at Melville, La., occurred prior to 1995. In spite of the general decrease in nutrient flux, the average size of the Gulf of Mexico hypoxic zone has increased between 1997 and 2007. The reasons for this are not clear but could be due to the type or nature of nutrient delivery. Whereas the annual flux of total nitrogen to the Gulf of Mexico has decreased, the proporti

Scientific Investigations Report↗

USGS2018

Santa Barbara and Foothill groundwater basins Geohydrology and optimal water resources management—Developed using density dependent solute transport and optimization models

Groundwater has been a part of the city of Santa Barbara’s water-supply portfolio since the 1800s; however, since the 1960s, the majority of the city’s water has come from local surface water, and the remainder has come from groundwater, State Water Project, recycled water, increased water conservation, and as needed, seawater desalination. Although groundwater from the Santa Barbara and Foothill groundwater basins only accounts for a small percentage of the long-term supply, it is an important source of supplemental water during times of surface-water shortages. During the late 1980s and early 1990s, production wells extracted additional groundwater to compensate for drought related water-delivery shortfalls from other sources; in response, water levels declined substantially in the Santa Barbara and Foothill groundwater basins (below sea level in the Santa Barbara groundwater basin). In coastal basins that have groundwater extraction near shore, seawater intrusion is often a problem. Seawater intrusion in the Santa Barbara groundwater basin is thought to be more limited than in other coastal basins because of an offshore fault that acts as a partial barrier to groundwater flow. During the late 1980s and early 1990s, seawater intrusion was observed in the Santa Barbara groundwater basin, as indicated by increased chloride concentrations at several monitoring wells that ranged from 200 ft to 1,300 ft from the ocean and as close as 2,900 ft to the nearest pumping well. This demonstrated that seawater can intrude into the Santa Barbara groundwater basin when groundwater levels fall below sea level near the coast. The city of Santa Barbara is interested in developing a better understanding of the sustainability of its groundwater supplies. In 2014, California adopted historic legislation to manage its groundwater: the Sustainable Groundwater Management Act (SGMA). The SGMA requires the development and implementation of “Groundwater Sustainability Plans” in 127 priority groundwater basins; although Santa Barbara was not a designated priority basin, the city is taking steps to achieve sustainability. Sustainability was defined in the SGMA in terms of avoiding undesirable results: significant and unreasonable groundwater-level declines, reduction in groundwater storage, seawater intrusion, water-quality degradation, land subsidence, and surface-water depletion. In this project, a cooperative study between the U.S. Geological Survey (USGS) and the city of Santa Barbara, sustainable yield is defined as the volume of groundwater that can be pumped from storage without causing water-level drawdowns and the associated increases in seawater intrusion (as indicated by increases in measured chloride concentrations) at selected wells. In order to estimate the sustainability of Santa Barbara’s groundwater basins, a three-dimensional density-dependent groundwater-flow and solute-transport model (the Santa Barbara Flow and Transport Model, or SBFTM) was developed on the basis of an existing groundwater-flow model. To simulate seawater intrusion to the Santa Barbara Basin under various management strategies, the SBFTM uses the USGS code SEAWAT to simulate salinity transport and variable-density flow. The completed SBFTM was coupled with a management optimization tool, in this case a multi-objective evolutionary algorithm, to determine optimal pumping strategies that maximize the sustainable yield and at the same time satisfy user-defined drawdown and chloride-concentration constraints. As part of this study, a three-dimensional hydrogeologic framework model was developed to quantify the extent and hydrogeologic characteristics of the Santa Barbara and Foothill groundwater basins and to help define the discretization and hydraulic properties used in the SBFTM. The development of the hydrogeologic framework model required the collection and reconciliation of geologic and geophysical data from existing maps, reports, and databases, along with geologic and hydrologic data from recently drilled wells. These data were integrated into a three-dimensional hydrogeologic framework model that defines the stratigraphy and geometry of the aquifer zones and the major geologic structures in the basin. The hydrogeologic framework model also quantifies the variation in sediment grain size within each aquifer zone as the percentage of coarse-grained sediment. Previous studies indicated that there are two principal water-producing zones in the Santa Barbara groundwater basin, the upper and lower producing zones; an additional thin, productive zone was identified as part of this study. This “middle producing zone” is not as areally extensive as the upper and lower producing zones and only exists in the coastal part of Storage Unit I. These producing zones are bounded at depth by less productive shallow, middle, and deep zones. Two versions of the SBFTM were constructed: an initial-condition model and a modern transient model. The initial-condition model is a long-term transient model that simulates flow and solute-transport conditions during a period with limited anthropogenic influences preceeding the modern transient model. The simulation-transient model simulates flow and transport conditions from 1929 through 2013; however, because of data availability, the focus of the model calibration was 1972–2013. The SBFTM was calibrated to measured groundwater levels and drawdown, as well as measured chloride concentrations and change in concentrations, using a combination of automated and trial-and-error parameter-estimation techniques. A sensitivity analysis indicated that, in general, the SBFTM was most sensitive to recharge- and pumping-distribution parameters, specifically those controlling the amount of small-catchment recharge and the distribution of water extraction by hydrogeologic layer for production wells. The model was also sensitive to parameters controlling stream-recharge rates, horizontal and vertical hydraulic conductivity, and porosity. From 1929 to 1971, most of the water entering the area represented by the SBFTM was from creek and small-catchment recharge, and the majority of water leaving the SBFTM area was from pumping, discharge to creeks, and drains. In addition, about 37 percent of the total pumpage came from a net reduction in groundwater storage. From 1972 to 2013, the amount of water entering and leaving the SBFTM was fairly similar as that from 1929 to 1971, except the reduction in pumpage added about 17,000 acre-ft of water to storage. During this later period, there were also times of storage loss. For example, during July 1990, a month when approximately 705 acre-ft of groundwater was pumped in the study area, the pumpage was much greater than all sources of recharge combined, and about 382 acre-ft of water was removed from groundwater storage. Simulated hydraulic heads replicated the observed data to an acceptable matching of the measured water-level, flow direction, and vertical gradients. Simulated hydrographs for selected wells were in good agreement with the measured data, with an average residual of -2.7 ft and a standard deviation of 14.5 ft, indicating that the simulated heads, on average, underestimated the observed water levels. An examination of the model fit indicated that most of the discrepancies were lower simulated heads at wells proximal to production well sites. The simulated chloride concentrations reasonably matched the rising limbs of the measured breakthrough curves in terms of timing and magnitude; however, the simulation overestimated the chloride concentrations on the falling limbs. The overestimation of low chloride concentrations was attributed to the model overestimating the advance of the chloride front during periods of heavy pumping and underestimating the retreat of the chloride front during periods of low pumping. These simulation errors would result in a conservative response by local water managers to seawater intrusion. The SBFTM was used to develop a collection of predictive simulations optimized to produce pumping schedules that maximize yield, subject to a set of constraints and competing objectives. The simulations were grouped as scenarios that differed in their time horizon, initial conditions for groundwater levels and chloride concentrations, as well as precipitation, which was incorporated into the model through simulated recharge. Overall, five scenarios were developed in a multi-objective framework to obtain optimal pumping rates for all of the wells managed by the city, while minimizing excessive drawdown and seawater intrusion. For the current study, complexities in the simulation model and the optimization formulation required additional considerations. Incorporating the solute-transport equations to simulate chloride transport added a highly nonlinear process that is solved iteratively in each time step of the groundwater-flow model. These nonlinearities, coupled with the highly refined grid in the current model, creates challenges for many traditional optimization methods. Therefore, an optimization method was needed that could address nonlinear relationships as well as a very large problem size. Lastly, the optimization problem was reformulated to include multiple objectives without requiring convergence to a single solution. This approach, guided by the city’s objectives, allowed the maximum extraction of information from the complex simulation. Borg, a multi-objective evolutionary algorithm, was chosen as the optimization algorithm for this study for several reasons: (1) it is very computationally efficient; (2) it can run in parallel; (3) it requires little user input; and (4) it can solve for multiple competing objectives. The first three points allow the algorithm to proceed toward the optimal solutions at the fastest possible rate. The fourth point is advantageous for large, complex optimization problems because it is difficult to formulate the optimization problem in a way that produces only one optimal solution. The problem formulation consisted of four competing objectives and a constraint set in accordance with the main concerns of the city. The objectives were maximizing total pumpage, minimizing seawater intrusion, minimizing total drawdown in production wells, and minimizing the maximum drawdown. The constraints were pump capacity, meeting drinking-water standards for chloride, maintaining a specified minimum flowrate to a groundwater treatment plant, and maintaining minimum water levels in pumping wells. The decision variables either were quarterly pumpage by well or total pumpage by basin. Five optimization scenarios were developed that allow the decision makers to evaluate a range of optimal solutions for a variety of water levels and chloride concentrations as well as potential future climatic conditions. Three scenarios (1, 2, and 5) were multi-objective optimization formulations that allowed for variations in management preferences and climatic conditions. The other two scenarios (3 and 4) were designed to examine the optimization results to answer specific questions. Scenario 1 described the best-case sustainable yield assuming a “full” basin (that is, high initial water levels) and typical climate conditions for 10 years. Scenario 2 also started with a “full” basin; however, this was followed by a 10-year drought. Scenario 3 determined if an “empty” basin (that is, low initial water levels) would recover to full conditions (1998 conditions) given climate assumptions and optimal pumping schedules from scenarios 1 and 2. Scenario 4 was designed to produce decision rules that can be used by water managers to help choose an optimal pumping schedule based on measured water-level or chloride data. Scenario 5 identified future pumping schedules based on short-term climate variations during a 2-year management horizon. The results from scenarios 1 and 2 described the differences in maximum pumpage in the basin under typical and dry long-term climate projections, respectively. The scenario 1 results indicated the maximum 10-year pumpage of the basin was about 31,300 acre-ft under typical conditions and controlling simulated seawater intrusion and drawdowns. For scenario 2, less recharge over the 10-year dry climate produced a maximum pumpage estimate of 30,000 acre-ft to control seawater intrusion and drawdowns. The larger pumpage for scenario 1 resulted in more seawater intrusion, but less total drawdown, compared to that of scenario 2. Results for scenarios 3 and 4 showed the basin’s response to management actions combined with climate projections. Both scenarios used the optimal pumping schedules and the 10-year climates from scenarios 1 and 2. The scenario 3 results showed that under minimal pumping, the basin did not fully recover to 1998 water levels within 10 years under either climate scenario. The relatively larger recharge from the typical climate resulted in less drawdown at coastal monitoring wells after the 10-year recovery period than that from the dry climate. The location of the seawater intrusion front was not appreciably different between the scenarios, however. Scenario 4 used the optimal results from scenarios 1 and 2 to produce decision-rule curves that illustrated the pumpage for each basin, given measured levels of chloride concentration or drawdown. This allowed the use of additional measurements at monitoring wells to assess future management decisions on the basis of the sensitivity of observations of drawdown and seawater intrusion to various pumping rates. Scenario 5 allowed managers to investigate the effects of short-term climate variations on optimal pumping schedules. Three specific 2-year simulations were optimized: typical-to-dry (scenario 5A), dry-to-typical (scenario 5B), and dry-to-dry (scenario 5C). The most noteable result from scenario 5 was the overall reduction in optimal pumpage for most schedules in scenario 5C, when the climate is simulated as dry-to-dry. There are also many optimal pumping schedules that produced an overall increase in waterlevels over the two-year simulation period, regardless of climatic condition. Similar to scenario 2, the scenario 5C results represents conservative yield estimates under a minimal-precipitation climatic condition.

California↗

USGS1998

Preliminary lithostratigraphy, interpreted geophysical logs and hydrogeologic characteristics of the 98th Street core hole, Albuquerque, New Mexico

Core samples, cuttings, and numerous geophysical logs obtained from the 1560 ft (475.5 m) core hole drilled at 98th Street on the west side of Albuquerque provide key stratigraphic and hydraulicproperty information for the upper clastic sediments of the Santa Fe Group, which form the principal aquifer in the region. The core hole and an adjacent water-level monitoring well were drilled cooperatively by the U.S. Geological Survey (USGS) and the City of Albuquerque and investigated in collaboration with the New Mexico Bureau of Mines and Mineral Resources and the New Mexico Office of the State Engineer to improve understanding of aquifer characteristics and controls on ground-water availability and quality. The 751.5 ft (229 m) of core samples recovered from the core hole are the only undisturbed samples of nonlithified sediments of the upper part of the Santa Fe Group that have been collected in this area. These samples have allowed us, for the first time, to directly observe and characterize the lithic and sedimentologic features of this part of the section, and to correlate the detailed geologic features with geophysical-log characteristics, magnetic susceptibility measurements, hydraulic variables, and trace-element geochemistry. The adjacent well was designed to be an areally representative ground-water level and water-quality monitoring well for the Santa Fe Group aquifer. This report chiefly addresses the lithologic, stratigraphic, and hydrogeologic features determined from the 98th Street core hole; other reports address related characteristics. Previous geologic studies predicted the stratigraphy at the site to be, from the land surface downward: 1) Quaternary alluvial and eolian valley-border sediments; 2) fluvial sand and gravel of the upper unit of the Santa Fe Group (Ceja Member of the Santa Fe Formation of Kelly, 1978; equivalent to the Sierra Ledrones Formation of Machette (1978a); 3) downward-fining basin-floor silty clay deposits and 4) fluvial sandy and silty facies of the middle unit of the Santa Fe Group (the Middle Red Member of Bryan and McCann, 1937, and Lambert, 1968). New geologic interpretations indicate that the drill site is in a fault block bounded by east-dipping normal faults and the oblique Atrisco-Rincon fault zone. Core-hole sampling recovered 760.6 ft (231.8 m) of core, in core segments 2.1-2.375 in. (5.3-6 cm) in diameter, and 0.2-10 ft (6.1 cm-3 m) long. The core hole was cased with centered 3-in. PVC casing, and is available for geophysical logging. The monitoring-well hole contains four piezometers at depths of 1544 ft (470.6 m), 1112 ft (338.9 m), 749 ft (228.3 m), and 458 ft (139.6 m). Sediments in the core are loose to weakly cemented gravel, sand, silt, and clay, and lithified sandstone. Laboratory analyses of particle-size distributions of 28 channel samples show that most silty sand samples are uniformly graded and poorly sorted; medium sand samples are moderately sorted. Six principal sediment types are used to describe the core; these sediment types are repeated in various combinations throughout the core and are used to define 22 lithologic units in the cored interval. The six principal sediment types contain sequences of beds having similar modal grain size and sedimentary structure, and are listed in decreasing abundance: 1) Silty fine sand, poorly sorted, containing a coarse silt matrix. Geophysical logs show highly variable baselines with deflections that are related to clay beds and sequences of silt, clay, and sorted fine sand. Density values of 2.12-2.25 g/cc and porosity values of 30-35 percent are typical. 2) Medium sand, moderately to poorly sorted. Geophysical logs show baselines of low variability with deflections that are related to clay beds and sequences of silt, clay, and sorted fine sand. Density values of 2.05-2.20 g/cc and porosity values of 30-35 percent are typical. 3) Clayey sandy silt, poorly sorted, locally microlaminated clay and silt, generally nonplastic. Geophysical logs show highly variable baselines with deflections that are related to sequences of clay and fine sand. Density values of 2.1-2.2 g/cc and porosity values of 30-40 percent are typical. 4) Silt and clay, characteristically red to reddish brown and medium to high plasticity, massive to indistinctly microlaminated. Geophysical logs show variable baselines with broad, high-amplitude compound spikes that are related to sequences of silt and fine sand. Density values of 2.12-2.25 g/cc and porosity values of >45 percent are typical. 5) Sand and gravel, poorly sorted. Geophysical logs show variable baselines with deflections that are related to sequences of silty and sorted fine sand. 6) Sandstone, fine-to-medium grained, poorly sorted, cemented chiefly by calcite, which fills the original pore space. Geophysical logs show density values >2.25 g/cc and porosity values <30 percent. The 22 lithologic units are correlated with recognized basin-floor fluvial lithofacies (Hawley, 1996), which include sand and gravel (lithofacies I), sand with lenses of pebbly sand, silt, and silty clay (lithofacies II), and interbedded sand, silt, and silty clay (modified lithofacies III, IV, IX). The sediments in the core hole are correlated with three informal lithostratigraphic units. The top unit, 0-19 ft (0-5.8 m) depth, consists of Quaternary eolian sand and valley-border alluvium. Coarsegrained deposits in the 19-97 ft (5.8-29.6 m) interval are correlated with the upper unit of the Santa Fe Group. The fine-grained section in the 97-787 ft (29.6-239.9 m) interval is correlated tentatively with the middle unit of the Santa Fe Group. This section contains thick sequences of laminated red and olivebrown clay and silt overbank deposits (441-787 ft) in the distinctive Atrisco member of Connell and others (1998). The Atrisco is correlated with fine-grained zones in numerous wells throughout the central Albuquerque metropolitan area, and is recognized as a zone that separates the upper Santa Fe aquifer from underlying middle Santa Fe deposits. The lower section of the middle unit of the Santa Fe, 787-1500 ft (239.9-457.2 m) depth, includes an upper sequence of moderately sorted channel-fill medium sand, and a lower sequence of sand, silt, and clay overbank deposits. The age of the cored interval is not known precisely. The upper Santa Fe gravel is related regionally to a through-flowing river system that was established in the Rio Grande rift valleys in Early Pliocene time, >4.5 MA. The middle Santa Fe unit is dated tentatively by correlation with a fossiliferous section, in which sandy beds that directly underlie the upper Santa Fe are Late Miocene (Hemphellian), 4.6- 8.9 MA. Further, the middle Santa Fe unit, with dominantly normal magnetic polarity, may have been deposited during closely spaced normal magnetic chrons 5.9-8.3 Ma. Four hydrostratigraphic units summarize the hydrogeologic framework for the 98th Street site: 1) Quaternary valley-border deposits, 2) upper Santa Fe sand and gravel deposits, 3) middle Santa Fe overbank deposits, and 4) middle Santa Fe channel-sand deposits. Empirical values of horizontal hydraulic conductivity estimated from core samples reveal a previously unknown contrast in hydraulic conductivity in the lowest two hydrostratigraphic units. Correlations among numerous wells show that the distinctively fine-grained Atrisco member, with estimated hydraulic conductivities (K) of <0.02-17 ft/day, is a laterally extensive barrier to vertical ground-water flow. The underlying unit that contains moderately sorted medium sand is a potential aquifer production zone that should be investigated further. Laboratory determination of vertical hydraulic conductivity values for fine-grained core samples range from 10 -2 to 10 -7 ft/day; recompacted sandy samples have K values of 1 to 10 -2 ft/day. Results of tests conducted with increasing effective stress show that K values of all samples decrease with decreasing porosity. Comparison of K values from laboratory, empirical, and calculated geophysical values shows discrepancies of 1-3 orders of magnitude (ft/day), indicating that additional analyses of core samples and geophysical data are necessary for future characterization of the Santa Fe Group aquifer.

New Mexico↗

USGS1985

Effects of flood control alternatives on fish and wildlife resources of the Malheur-Harney lakes basin

Malheur Lake is the largest freshwater marsh in the western contiguous United States and is one of the main management units of the Malheur National Wildlife Refuge in southeastern Oregon. The marsh provides excellent waterfowl production habitat as well as vital migration habitats for birds in the Pacific flyway. Water shortages have typically been a problem in this semiarid area; however, record snowfalls and cool summers have recently caused Malheur Lake to rise to its highest level in recorded history. This has resulted in the loss of approximately 57,000 acres of important wildlife habitat as well as extensive flooding of local ranches, roads, and railroad lines. Because of the importance of the Refuge, any water management plan for the Malheur-Harney Lakes Basin needs to consider the impact of management alternatives on the hydrology of Malheur Lake. The facilitated modeling workshop described in this report was conducted January 14-18, 1985, under the joint sponsorship of the Portland Ecological Services Field Office and the Malheur National Wildlife Refuge, Region 1, U.S. Fish and Wildlife Service (FWS). The Portland Field Office is responsible for FWS reporting requirements on Federal water resource projects while the Refuge staff has management responsibility for much of the land affected by high water levels in the Malheur-Harney Lakes Basin. The primary objective of the workshop was to begin gathering and analyzing information concerning potential fish and wildlife impacts, needs, and opportunities associated with proposed U.S. Army Corps of Engineers (COE) flood control alternatives for Malheur Lake. The workshop was structured around the formulation of a computer model that would simulate the hydrologic effects of the various alternatives and any concommitant changes in vegetation communities and wildlife use patterns. The simulation model is composed of three connected submodels. The Hydrology submodel calculates changes in lake volume, elevation, and surface area, as well as changes in water quality, that result from the proposed water management projects (upstream storage, upstream diversions, drainage canals) and the no action alternative. The Vegetation submodel determines associated changes in the areal extent of wetland and upland vegetation communities. Finally, the Wildlife submodel calculates indices of abundance or habitat suitability for colonial nesting birds (great egret, double-crested cormorant, white-faced ibis), greater sandhill crane, diving ducks, tundra swan, dabbling ducks, and Canada goose based on hydrologic and vegetation conditions. The model represents the Malheur-Harney Lakes Basin, but provides water quantity and quality indicators associated with additional flows that might occur in the Malheur River Basin. Several management scenarios, representing various flood control alternatives and assumptions concerning future runoff, were run to analyze model behavior. Scenario results are not intended as an analysis of all potential management actions or assumptions concerning future runoff. Rather, they demonstrate the type of analysis that could be conducted if the model was sufficiently refined and tested. Early in a model development project, the process of building the model is usually of greater benefit than the model itself. The model building process stimulates interaction among agencies, assists in integrating existing information, and helps identify research needs. These benefits usually accrue even in the absence of real predictive power in the resulting model. This workshop initiated interaction among the primary State and Federal resource and development agencies in a nonadversarial forum. The exchange of information and expertise among agencies provided the FWS with the best information currently available for use in the Planning Aid Letter it will develop at the Reconnaissance state of the COE study. If the COE subsequently initiates a Feasability Study, this information will be refined further and will aid the FWS in preparing its Coordination Act Report on any flood control alternative proposed by the COE. The model building and testing process also helped identify model limitations and more general information needs that should be evaluated for further study prior to preparation of an FWS Coordination Act Report. Major needs associated with the Hydrology submodel include a more detailed representation of hydrologic units (separately consider Harney Lake, Mud Lake, and Malheur Lake or the three hydrological units within Malheur Lake, rather than a combined lake system) and explicitly representation of groundwater storage and discharge in water budget calculations. A better representation of the hydrological units will require more detailed topographic data for the basin, capacity-elevation and elevation-surface area curves for each unit, and better water flow data between the units. Additional water quality parameters and constraints on proposed canal operation due to conditions in the Malheur River might also be added. Key Vegetation submodel needs include fine-tuning existing vegetation relationships in the model and adding relationships to address the influence of historical conditions on vegetation development, effects of very rapid changes in lake level, effects of wildlife populations (e.g., carp, muskrat), responses of vegetation to habitat management actions (e.g, haying, grazing, burning), and better representation of sago pondweed dynamics. A complementary geographic information system might also be developed for spatial analyses. Major needs that should be evaluated for the Wildlife submodel include addition of other wildlife species that have important effects on habitat on the Refuge (e.g., carp, muskrat) and consideration of additional life-cycle requisites and controlling variable for species presently in the model. Some of these limitations could perhaps be overcome if historical data on habitat conditions were developed to use with historical data on wildlife populations.

Oregon↗

USGS1989

Geology of the Arabian Peninsula; shield area of western Saudi Arabia

Western Arabia lies within the low-latitude desert of north Africa and the Middle East, the core being the Arabian segment of the African Shield. The core of complex basement rocks accounts for about 670,000 km 2 , or one-third of the Arabian Peninsula. Reconnaissance mapping of these crystalline rocks, together with bordering sedimentary rocks and volcanic flows, begun in 1950, resulted during the next 13 years in a series of geologic and geographic maps without extensive texts. The maps served as general guides for development of natural resources, including water supplies, ore deposits, and building materials. An intensive exploration program that began in 1963 and involved numerous geologists has vastly increased geologic information. Rainfall in Arabia is meager and episodic, and vegetation is sparse except in isolated copses on the crest of the Hejaz Range. Comparison of flora with similar species in the Sudan, where records of rainfall have long been kept, allows evaluation of mean annual precipitation. Wandering bedouin following fodder created a delicate balance between population and water supply-now disturbed by wells drilled in alluvium and lava fields. A trapezoidal region of Precambrian crystalline rocks lies along the northeast flank of the Red Sea, with two long prongs extending northwest and southeast for a total of 1,800 km. These basement rocks of the Arabian Shield are well exposed on the uplands, scarp mountains, and coastal pediments where the Phanerozoic cover rocks have been stripped as a result of Paleozoic epeirogeny and Tertiary ramping. The shield outcrops are divided into three tectonic provinces by N. 45°W.- trending shear zones of the Najd fault system of latest Proterozoic and possibly earliest Paleozoic time. The southwestern province, the 'Asir1 upland, was sharply uplifted and tilted to the northeast during the Neogene. The northwestern province, consisting of the Ash Shifa'- Hisma upland as well as Jabal Shammar farther east, similarly was uplifted and tilted. These two provinces are separated by the flat-lying median N ajd province, which is chiefly bounded by the principal Najd faults. The outcrops of the shield rocks are of the Late Proterozoic Eonupper Riphean to Vendian or Infracambrian epochs, including the Ediacarian System. The most reliable isotopic ages range from about 900 to 560 m.y., but some Middle Proterozoic rocks may be present in the easternmost shield. The rocks are divided into six lithostratigraphic sequences, two plutonic suites, and an ophiolitic suite. The mafic and ultramafic volcanic and plutonic rocks of the ophiolitic suite everywhere were emplaced tectonically and are probably of different ages in different places. Some ophiolite occurs as obducted blocks, but most is highly deformed and altered to serpentinite in fault zones that mostly define sutures between different tectonic blocks or terranes within the shield. Three of the lithostratigraphic sequences consist of mafic to silicic volcanic rocks and volcanic-derived clastic rocks which, with their subvolcanic plutonic rocks of a dioritic suite, probably formed in oceanic island arcs during convergent plate tectonism. These rocks make up the primary, or first-formed, crust of the shield. Chemical analyses show that the primary shield rocks, regardless of age, are principally calc-alkalic with some associated tholeiitic varieties. Most of the layered rocks are andesitic, but they range from basalt to dacite and in places contain intercalated pillow basalt, marble, chert, and carbonaceous or graphitic schist. Most of the plutonic rocks of the dioritic suite are dioritic, but they range from gabbro to trondhjemite and rarely contain potassium feldspar. The sequences and an associated dioritic suite become younger toward the eastern shield, that is, the primary crust of the shield youngs toward the east. Two western sequences consist of the Jiddah (Samran) and BaishBahah Groups and range in radiometric age from about 900 to 800 m.y.; the eastern sequence consists of the Halaban (Hulayfah) Group and ranges from 800 to about 700 m.y. During subsequent orogeny, most of the rocks were intensely deformed and mostly metamorphosed to upper greenschist facies, but rising in places to the almandine-amphibolite facies. Two other lithostratigraphic sequences with an associated plutonic granitic suite are the products of two mountain-building episodes during which the primary crust was greatly thickened and converted into craton. The two sequences, including largely the Ablah (Al Ays) and Murdama (Shammar) Groups, consist of abundant sedimentary rocks, commonly arkosic, that are the erosional products of the orogenic mountains. They are several thousand meters thick. Less abundant calc-alkalic to alkalic volcanic rocks, commonly dacitic and rhyolitic, are intercalated with the sedimentary rocks. The plutonic rocks of the granitic suite in association with both sequences have syntectonic and posttectonic phases, are products of the orogenies, and are the principal new ingredients making up the craton. Gneiss domes were a significant part of these cratonization orogenies. In association with orogenic crustal heating, some of the low-density, more silicic tonalitic and trondhjemitic rocks of the primary crust rose as gneiss domes. Partial melting in the middle or lower crust below the gneiss domes produced large volumes of granitic magma that intruded the gneiss domes as granodioritic batholiths. The Ablah Group and the older part of the granitic suite are about 775 to 740 Ma old and are associated with the Ablah orogeny and early cratonization in the western and earlier formed half of the shield. The Murdama (Shammar) Group and the younger part of the granitic suite are about 660 to 580 Ma old and are associated with the culminant orogeny and late cratonization that was shieldwide. The granitic suite during both orogenies consists of early, syntectonic granodiorite batholiths associated with the gneiss domes and late, posttectonic monzogranite plutons. Only during the culminant orogeny, late magmatic evolution produced syenogranite and alkali-feldspar granite commonly in circular and ring-structured plutons and with associated explosive volcanic deposits (Shammar Group); final products, some of which have economic potential, were peralkalic and peraluminous. The late plutonism of the culminant orogeny was distinctly bimodal in that subordinate gabbroic rocks are associated with the granites. Various building blocks or terranes of the andesitic and dioritic primary crust were collisionally agglomerated during the Ablah orogeny, early cratonization, whereas the entire shield as currently exposed was further collisionally accreted and compressionally consolidated during the culminant orogeny, final cratonization. Thousands of kilometers of oceanic crust had to be subducted in about 300 m.y. to form the large primary crust of the Arabian Shield. The inevitable collisional events during consumption of such a large volume of oceanic crust invariably led to numerous collisional orogenies that collectively encompass the widely known Pan African tectonic episode. The youngest lithostratigraphic sequence, the Jubaylah Group, is essentially postcratonic, although it is the end product of the collisional culminant orogeny. Final east-west compression of the entire shield from about 580 to 560 m.y. caused the craton to fracture along the large northwest-trending, left-lateral faults and elsewhere along lesser, northeast-trending, right-lateral, conjugate faults of the N ajd fault system. Erosional products of this more localized deformation were the sedimentary rocks of the Jubaylah Group, which also includes intercalated andesitic to basaltic volcanic rocks of a mafic alkalic compositional trend. The collisional edge of an old continental plate (or tectonic fragments thereof), suspected on the eastern edge of the Arabian Shield, has not been shown with certainty to be exposed. Presumably, widespread contamination from such an old continental crust affects U/Pb, Sm/Nd, Rb/Sr, and common lead ratios in the young plutonic rocks of the easternmost shield. One mass of anorthosite near Jabal Khida' on the central eastern edge of the shield may be a fragment of this old continental plate in that associated granodiorite may be as old as 1,600 to 1,800 Ma. Epeirogenic uplift, erosion, and cooling of the uppermost shield during Early and Middle Cambrian time is indicated by an average fission track age of 510±52 m.y., on sphene from diorite (hornblende K-Ar age of 615±12 m.y.) in the southwestern part of the shield. The hiatus was followed by extensive deposition of the Cambro-Ordovician Saq Sandstone in the north and northeast and the Wajid Sandstone in the southeast and south of the shield. The Cambrian Siq Sandstone had already been deposited in the northern part. During the middle and late Paleozoic, broad epeirogeny caused further erosion of the shield until marine transgression deposited the Upper Permian Khuff Formation at least in the eastern part of the shield. In the southwestern shield, the nonmarine Upper Triassic Khums Sandstone was deposited variably on Wajid or Precambrian rocks and is overlain by limestone of the middle Upper Jurassic Amran Formation. Except for shallow marine sandstone of problematic Cretaceous age deposited on the Amran Formation in the south.western shield and on Precambrian rocks in the northwestern shield, the younger beds on the shield are Paleocene and younger, with the possibility that the lowermost are upper Maestrichtian. The early Tertiary beds contain vertebrate fossils of coastal marine or estuarine environment 250 km east of the Red Sea in the central shield. Marginal marine sediments were deposited in a western tongue of the latest Tethys Sea as late as Eocene on the western shield and at least as far south as Jiddah. The great harrats of flood basalt erupted on th~ western shield during late Oligocene and early Miocene at the same time a 2,000-kmlong continental rift valley developed along the future Red Sea axis. Within this rift valley, Baid freshwater tuffaceous lakebeds were deposited between mafic and silicic volcanoes. During late early Miocene time, the Red Sea opened at a rate of 4.4 cm/yr in a firststage movement while continental dikes and swarms of oceanic tholeiitic dikes, gabbro, and granophyre plutonic rocks were intruded into the rift sedimentary and volcanic rocks at the newly formed continental margin. The continental margin was deformed and greatly extended at this time. About 14 or 15 m.y., as the first-stage spreading stopped, the Red Sea Escarpment rose; its erosion caused deposition of coarse conglomerate of the Bathan Formation. About 3,000 m of evaporite was deposited on the young Red Sea oceanic crust during the late Miocene desiccation crisis. A second stage of sea-floor spreading about 4-5 m.y. produced the Red Sea axial trough, consisting of oceanic crust, as well as renewed uplift and tilting of the three tectonic provinces in response to compression from counterclockwise rotation against the Dead Sea Rift. This late movement caused widespread major stream capture, especially along the wadis that formerly drained southwesterly or northwesterly, the channels turning westward through narrow gorges to the coastal plain and the Red Sea.

Professional Paper↗

USGS2003

Hydrology and water quality of Elkhead Creek and Elkhead Reservoir near Craig, Colorado, July 1995–September 2001

The U.S. Geological Survey, in cooperation with the Colorado River Water Conservation District, collected and analyzed baseline streamflow and water-quality information for Elkhead Creek and water-quality and trophic-state information for Elkhead Reservoir from July 1995 through September 2001. In the study area, Elkhead Creek is a meandering, alluvial stream dominated by snowmelt in mountainous headwaters that produces most of the annual discharge volume and discharge peaks during late spring and early summer. During most of water year 1996 (a typical year), daily mean discharge at station 09246400 (downstream from the reservoir) was similar to daily mean discharge at station 09246200 (upstream from the reservoir). Flow-duration curves for stations 09246200 and 09246400 were nearly identical, except for discharges less than about 10 cubic feet per second. Specific conductance generally had an inverse relation to discharge in Elkhead Creek. During late fall and winter when discharge was small and derived mostly from ground water, specific conductance was high, whereas during spring and early summer, when discharge was large and derived mostly from snowmelt, specific conductance was low. Water temperatures in Elkhead Creek were smallest during winter, about 0.0 degrees Celsius ( o C), and largest during summer, about 20–25 o C. Concentrations of major ions, nutrients, trace elements, organic carbon, and suspended sediment in Elkhead Creek indicated no substantial within-year variability and no substantial differences in variability from one year to the next. A seasonal pattern in the concentration data was evident for most constituents. The seasonal concentration pattern for most of the dissolved constituents followed the seasonal pattern of specific conductance, whereas some nutrients, some trace elements, and suspended sediment followed the seasonal pattern of discharge. Statistical differences between station 09246200 (upstream from the reservoir) and station 09246400 (downstream from the reservoir) were indicated for specific conductance, dissolved calcium, magnesium, sodium, and sulfate, acid-neutralizing capacity, and dissolved solids. Trend analysis indicated upward temporal trends for pH, dissolved ammonia plus organic nitrogen, total nitrogen, and total phosphorus at station 09246200; upward temporal trends for dissolved and total ammonia plus organic nitrogen, total nitrogen, and total phosphorus were indicated at station 09246400. No downward trends were indicated for any constituents. Annual loads for dissolved constituents during water years 1996–2001 were consistently larger at station 09246400 than at station 09246200, except for silica and sulfate. Mean monthly loads for dissolved constituents followed the seasonal pattern of discharge, indicating that most of the annual loads were transported during March–June. Annual dissolved nutrient loads at stations 09246400 and 09246200 were not substantially different, except for total phosphorus and total nitrogen loads, which were smaller at the downstream station than at the upstream station, most likely due to biological uptake and settling in the reservoir. Mean annual suspended-sediment load during water years 1996–2001 was about 87-percent smaller at the downstream station than at the upstream station. Temperature in Elkhead Reservoir varied seasonally, from about 0 o C during winter when ice develops on the reservoir to about 20 o C during summer. Specific conductance varied from minimums of 138 to 169 microsiemens per centimeter at 25 o C (µS/cm) during snowmelt inflow to maximums of 424 to 610 µS/cm during early spring low flow (April). Median pH in the reservoir ranged from 7.2 to 8.0 at all sites near the surface. Median dissolved oxygen ranged from 7.1 to 7.2 milligrams per liter (mg/L) in near-surface samples and from 4.8 to 5.6 mg/L in near-bottom samples. During reservoir stratification, specific conductance generally was largest in the epilimnion, resulting from warm and relatively concentrated water from Elkhead Creek that was routed through the reservoir in the relatively warm epilimnion. The pH in the epilimnion generally increased from May to September, probably a result of algal productivity. In the hypolimnion, pH decreased slightly with depth in the July and September, probably a result of biomass decay processes and a lack of circulation during stratification. Concentrations of nutrients in both near-surface and near-bottom samples from Elkhead Reservoir were highest during snowmelt inflow (April–May). Total phosphorus concentrations in near-surface samples generally were largest during runoff, whereas total phosphorus concentrations in near-bottom samples generally were largest during July or September. Concentrations of nitrite plus nitrate in near-surface samples were substantially depleted by biological uptake during July, September, and October, compared to near-bottom samples. Variations in concentration of chlorophyll- a in near-surface samples were large during the growing season with peak seasonal concentrations during runoff or late summer and fall. Trophic state for Elkhead reservoir ranged from oligotrophic to eutrophic.

Colorado↗

USGS1960

Geology and mineral deposits of the St. Regis-Superior area, Mineral County, Montana

The St. Regis-Superior area occupies about 300 square miles in northwestern Montana and includes parts of the Squaw Peak Range and Coeur d'Alerie Mountains of the northern Rocky Mountains physiographic province. Nearly 50,000 feet of metasedimentary rocks of the Precambrian Belt series, chiefly varieties of quartzite and argillite, underlies most of the area. The Belt series is informally subdivided with reference to the top of the Wallace formation into lower and upper parts. In this area, the lower part of the Belt series is divided into the Prichard, Burke and Revett, St. Regis, and Wallace formations, in order of decreasing age, and the upper part of the Belt series or the Missoula group is divided »into the Spruce, Lupine, Sloway, and Bouchard formations, and an unnamed feldspathic quartzite at Rock Rabbit Ridge, also from oldest to youngest. Formations in the lower part of the Belt series are correlated with formations of the same names in the Coeur d'Alene district, and formations in the upper part of the Belt series are tentatively correlated in part with formations of the Missoula group in the vicinity of Missoula, Mont. Paleozoic quartzite, shale, limestone, and dolomite crop out in several localities in the southeastern part of the area. The limestone unit contains fragments of a single species of Glossopleura of early Middle Cambrian age which, together with lithologic similarities, has been used to correlate at least the quartzite, shale, and limestone part of this Paleozoic sequence with the Flathead sandstone, Gordon shale, and Damnation limestone sequence known elsewhere in northwestern Montana. Several small diabasic dikes and sills are present in the area, generally associated with northwestward-trending faults. The major faults generally trend northwestward and are considered to be part of the Lewis and Clark structural line. The Osburn fault, the major element of the Lewis and Clark line through the Coeur d'Alene district and western Mineral County, has been traced to the east edge of the St. Regis- Superior area, thus extending the mapped length of the fault to about 100 miles. Evidence indicates that this major fault has diminished in intensity in this area and that most of the stress has been relieved along the Boyd Mountain fault that apparently splits from the Osburn fault a few miles west of St. Regis. Stratigraphic and structural evidence indicates a strike-slip right-lateral movement of about 3 miles along the Osburn fault. Horizontal stratigraphic separation along the Boyd Mountain fault indicates a right-lateral movement of about 13 miles. Low-grade regional metamorphism of the sedimentary rocks in the area has caused the recrystallization of quartz grains and the formation of sericite. Argillite and quartzite have been converted to phyllite and foliated quartzite by dynamic metamorphism in a large area north of the central part of the Osburn fault. Some of the shear zones contain a large amount of introduced carbonate minerals. From 1901 through 1953 this area has produced 8,086,827 pounds of zinc, 7,932,958 pounds of lead, 2,053,715 pounds of copper, 584,168 fine ounces of silver, and 588 fine ounces of gold. The lead, sine, and silver have come chiefly from veins in highly foliated rocks near the Osburn fault zone. The attitudes of these veins are controlled in large part by the cleavage. The principal ore minerals are galena, sphalerite, tetrahedrite, and boulangerite, and the gangue minerals are quartz, carbonate minerals, and barite. Most of the copper has come from the Amador vein where chalcopyrite and bornite are the principal ore minerals, and the gangue minerals are pyrite, quartz, and carbonate minerals. The Amador vein occurs in a belt of copper deposits that extends westward into the Coeur d'Alene district. These copper deposits may be genetically associated with diorite dikes and sills lying within the same belt. Fluorspar has been found in three closely spaced prospects along a northward- trending zone of brecciation and small-scale folding in Dry Creek valley. Incomplete production records show that 781 tons of fluorspar has been shipped from 2 of these prospects.

Montana↗

USGS2010

Discriminating silt-and-clay from suspended-sand in rivers using side-looking acoustic profilers

The ability to accurately monitor suspended-sediment flux in rivers is needed to support many types of studies, because the sediment that typically travels in suspension affects geomorphology and aquatic habitat in a variety of ways (e.g. bank and floodplain deposition, bar morphology, light penetration and primary productivity, tidal wetland deposition in the context of sea-level rise, sediment-associated contaminants, reservoir sedimentation and potential erosion during dam removal, among others). In addition, human-induced changes to the landscape have resulted in substantially altered suspended-sediment loads (Syvitski et al., 2005). Thus, accurate monitoring of suspended-sediment flux is necessary for informed resource management of rivers. Because of this need, a variety of techniques have been developed and applied for suspendedsediment monitoring. The traditional approach in the United States, which was developed and has been used extensively by the U.S. Geological Survey (USGS), is to collect an isokinetic, velocity-weighted sample from a river cross-section, analyze the sample in the laboratory, and use water-discharge records to compute a record of suspended-sediment flux (Guy, 1969, Guy, 1970, Edwards and Glysson, 1999, Porterfield, 1972). The labor and expense associated with this traditional approach is substantial such that the number of USGS gages reporting daily records of suspended-sediment flux decreased from 364 in 1981 to 120 in 2003 (Osterkamp et al., 2004). Also, the traditional sampling approach is limited with respect to the temporal resolution that can be achieved, thus requiring the use of approximate relations between suspended-sediment concentration and water discharge to fill gaps between samples. To address these limitations, several indirect or "surrogate" measures have been investigated (see e.g. Gray and Gartner, 2009) most notably optical backscatter (i.e. turbidity), laser-diffraction, and acoustic backscatter. These indirect techniques rely on measurements of ancillary properties that correlate with suspended-sediment concentration and particle size and thus require the collection of traditional samples for calibration. Through in situ deployments, these methods can provide the high temporal resolution that cannot be achieved through traditional sampling. Here we focus on the evaluation of acoustic profiling techniques (e.g. acoustic-Doppler sideways-looking profilers, or ADPs). One major advantage of acoustic profiling is the ability to concurrently measure water velocity (using Doppler-shift methods) and suspended-sediment concentration such that suspended-sediment flux can be directly computed using data from a single instrument. Acoustic-Doppler profilers have become popular for measuring water velocity and discharge in rivers, through both moving-boat operations and from fixed deployments such as bank-mounted sideways-looking instruments (Hirsch and Costa, 2004, Muste et al., 2007). The method presented herein is most suited to sideways-looking applications as a complement to the "index velocity" technique, whereby an index velocity from a sideways-looking instrument is related to the cross-section average velocity (determined from moving-boat discharge measurements) as a means for developing a continuous water-discharge record (Ruhl and Simpson, 2005). Topping et al. (2007) presented a method for discriminating silt-and-clay from suspended sand, using single frequency ADPs. This method takes advantage of the relations among acoustic backscatter, sediment-induced acoustic attenuation, suspended-sediment concentration (SSC), and particle size distribution (PSD). Backscatter is the amount of sound scattered back and received at the transducer while sediment-induced attenuation is the amount of sound scattered in other directions and absorbed by the sediment particles. Both of these parameters can be measured with an ADP, and their different dependencies on SSC and PSD allow for the discrimination of suspended silt-and-clay from suspended sand. Topping et al. (2007) describe application of the method at several sites along the Colorado River in Grand Canyon, and herein we present an example application of the technique for the Gunnison River, CO. However, the methods general applicability in rivers has yet to be evaluated due to a lack of concurrent acoustic and sediment data at a range of sites. To this end, the objective of the analysis presented herein is to evaluate the potential general applicability of the method, drawing from the extensive USGS database on SSC and PSD. We refer to it as "potential" general applicability because it relies on the theory underlying the previous empirical results. Use of the theoretical relations is necessary due to the lack of concurrent ADP and SSC/PSD data, but also serves the additional purpose of providing further justification of the empirical calibrations developed for the Colorado and Gunnison Rivers.

Conference Paper↗