Geology ReportsSearch

SEARCH · Geology Reports

Results for “Issue Paper”

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

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

At least 577 records · Page 32Linked to original sources

Applications of optical sensors for high-frequency water-quality monitoring and research

The recent commercial availability of in-situ optical sensors, together with new techniques for data collection and analysis, provides the opportunity to monitor a wide range of water-quality constituents over time scales during which environmental conditions actually change. Traditional approaches for data collection (daily to monthly discrete samples) are often limited by high sample collection, processing, and analytical costs, difficult site access, and logistical challenges, particularly for long-term sampling at a large number of sites. Optical sensors that continuously measure constituents in the environment by absorbance or fluorescence properties (Figure 1) have had a long history of use in oceanography for measuring highly resolved concentrations and fluxes of organic matter, nutrients, and algal material. However, much of the work using commercially-available optical sensors in rivers and streams has taken place in only the last few years. Figure 1. [NOT SHOWN] Optical sensor technology is now sufficiently developed to warrant broader application for research and monitoring in coastal and freshwater systems, and the United States Geological Survey (a U.S. science agency) is now using these sensors in a variety of research and monitoring programs to better understand water quality in-situ and in real-time. Examples are numerous and range from the applications of nitrate sensors for calculating loads to estuaries susceptible to hypoxia (Pellerin et al., 2014) to the use of fluorometers to estimate methymercury fluxes (Bergamaschi et al., 2011) and disinfection byproduct formation (Carpenter et al., 2013). Transmitting these data in real-time provides information that can be used for early trend detection, help identify monitoring gaps critical for water management, and provide science-based decision support across a range of issues related to water quality, freshwater ecosystems, and human health. Despite the value of these sensors, collecting data that meet high-quality standards requires investment in and adherence to tested and established methods and protocols for sensor operation and data management (Pellerin et al., 2013). For example, optical sensor measurements can be strongly influenced by a variety of matrix effects, including water temperature, inner filtering from highly colored water, and scattering of light by suspended particles (Downing et al., 2012). Characterizing and correcting sensors for these effects – as well as the continued development of common methodologies and protocols for sensor use – will be critical to ensuring comparable measurements across sites and over time. In addition, collaborative efforts such as the Nutrient Sensor Challenge (www.nutrients-challenge.org) will continue to accelerate the development, production and use of affordable, reliable and accurate sensors for a range of environments. REFERENCES Bergamaschi .B.A., Fleck J.A., Downing B.D., Boss E., Pellerin B.A., Ganju N.K., Schoellhamer D.H., Byington A.A., Heim W.A., Stephenson M., Fujii R. (2011), Methyl mercury dynamics in a tidal wetland quantified using in situ optical measurements. Limnology and Oceanography, 56(4): 1355-1371. Carpenter K.D., Kraus T.E.C., Goldman J.H., Saraceno J., Downing B.D., Bergamaschi B.A., McGhee G., Triplett T. (2013), Sources and Characteristics of Organic Matter in the Clackamas River, Oregon, Related to the Formation of Disinfection By-products in Treated Drinking Water: U.S. Geological Survey Scientific Investigations Report 2013–5001, 78 p. Downing .B.D., Pellerin B.A., Bergamaschi B.A., Saraceno J., Kraus T.E.K. (2012), Seeing the light: The effects of particles, temperature and inner filtering on in situ CDOM fluorescence in rivers and streams. Limnology and Oceanography: Methods, 10: 767-775. Pellerin B.A., Bergamaschi B.A., Downing B.D., Saraceno J., Garrett J.D., Olsen L.D. (2013), Optical Techniques for the Determination of Nitrate in En

Conference Paper

Tin

Tin (Sn) is one of the first metals to be used by humans. Almost without exception, tin is used as an alloy. Because of its hardening effect on copper, tin was used in bronze implements as early as 3500 B.C. The major uses of tin today are for cans and containers, construction materials, transportation materials, and solder. The predominant ore mineral of tin, by far, is cassiterite (SnO 2 ). In 2015, the world’s total estimated mine production of tin was 289,000 metric tons of contained tin. Total world reserves at the end of 2016 were estimated to be 4,700,000 metric tons. China held about 24 percent of the world’s tin reserves and accounted for 38 percent of the world’s 2015 production of tin. The proportion of scrap used in tin production is between 10 and 25 percent. Unlike many metals, tin recycling is relatively efficient, and the fraction of tin in discarded products that get recycled is greater than 50 percent. Only about 20 percent of the world’s identified tin resources occur as primary hydrothermal hard-rock veins, or lodes. These lodes contain predominantly high-temperature minerals and almost invariably occur in close association with silicic, peraluminous granites. About 80 percent of the world’s identified tin resources occur as unconsolidated secondary or placer deposits in riverbeds and valleys or on the sea floor. The largest concentration of both onshore and offshore placers is in the extensive tin belt of Southeast Asia, which stretches from China in the north, through Thailand, Burma (also referred to as Myanmar), and Malaysia, to the islands of Indonesia in the south. Furthermore, tin placers are almost always found closely allied to the granites from which they originate. Other countries with significant tin resources are Australia, Bolivia, and Brazil. Most hydrothermal tin deposits belong to what can be thought of as a superclass of porphyry-greisen deposits. The hydrothermal tin deposits are all characterized by a close spatial, temporal, and genetic association with highly differentiated, peraluminous porphyritic granite intrusions. The intrusions form pegmatites; disseminated ore; parallel or subparallel, greisen-bordered sheeted veins that either cross-cut the intrusion or are peripheral to it; skarns; and (or) limestone replacements that contain different amounts of cassiterite, molybdenite, and wolframite. The tectonic settings of tin-bearing granites are relatively well understood and of limited variety. Tin and tungsten deposits and their associated igneous rocks are found mainly in continental settings. Historically, prospecting for tin has been carried out by the time-honored methods of panning, drilling, trenching, and assaying. Geophysical and geochemical surveys have been employed to cover large areas more rapidly, isolating areas of possible tin deposits so that drilling can be more effective and less costly. Elemental concentrations and relationships of the lithophile elements, especially barium, lithium, niobium, potassium, rubidium, and zirconium, are the most reliable chemical indicators of ore-forming processes and tin-bearing potential. The average human diet includes an intake of about 10 milligrams per day of tin. Ingestion of tin in significantly greater amounts than 10 milligrams per day may lead to a stomach ache, anemia, and liver and kidney problems. Exposure to some organo-tin compounds can interfere with brain and nervous system function and, in severe cases, can cause death. Extended inhalation of tin oxide—an issue mainly for those people who work in the tin industry—results in a higher potential to develop stannosis, which is a mild disease of the lungs caused by the inhalation of tin-bearing dust. Inorganic tin is poorly absorbed by the body, and no evidence exists for the carcinogenicity of metallic tin and tin compounds in humans. Most placer tin deposits are mined by open pit and (or) dredging methods. Mining of alluvial placers in modern streambeds and riverbeds is likely to increase the amount of sediment delivered downstream. This, combined with potential diversion of rivers and streams, may negatively affect downstream ecosystems. Many of the placer deposits located in Burma, Indonesia, Malaysia, and Thailand are located offshore. Most offshore placer tin deposits are mined by dredging methods, which have the potential to negatively affect benthic, midwater, and pelagic ecosystems. In a congressionally mandated U.S. Department of Defense study of strategic minerals published in 2013, tin has the greatest shortfall amount (insufficient supply to meet demand) at \$416 million; this amount is more than twice that of antimony ($182 million), which is the strategic mineral with the next largest shortfall amount (U.S. Department of Defense, 2013). The United States imported 75 percent of its tin supply in 2015. During the period 2012–15, these imports were from, in descending order of amount imported, Peru, Indonesia, Malaysia, and Bolivia. A promising advancement concerning research into the origin of tin deposits is the recent development of a reliable method of analyzing tin isotopes in cassiterite. Although the mechanism of transport and deposition of tin is fairly well understood, the means by which tin is incorporated into the parent magma at the points of magma generation and ascent needs further investigation. Tin metallogenic provinces worldwide are well known. Consequently, any undiscovered tin deposits will likely be spatially close to known deposits or extensions of the same.

Professional Paper

Meeting summary - Coastal meteorology and oceanography: Report of the third prospectus development team of the U.S. Weather Research Program to NOAA and NSF

U.S. Weather Research Program (USWRP) prospectus development teams (PDTs) are small groups of scientists that are convened by the USWRP lead scientist on a one-time basis to discuss critical issues and to provide advice related to future directions of the program. PDTs are a principal source of information for the Science Advisory Committee, which is a standing committee charged with the duty of making recommendations to the Program Office based upon overall program objectives. PDT-1 focused on theoretical issues, and PDT-2 on observational issues; PDT-3 is the first of several to focus on more specialized topics. PDT-3 was convened to identify forecasting problems related to U.S. coastal weather and oceanic conditions, and to suggest likely solution strategies. There were several overriding themes that emerged from the discussion. First, the lack of data in and over critical regions of the ocean, particularly in the atmospheric boundary layer, and the upper-ocean mixed layer were identified as major impediments to coastal weather prediction. Strategies for data collection and dissemination, as well as new instrument implementation, were discussed. Second, fundamental knowledge of air-sea fluxes and boundary layer structure in situations where there is significant mesoscale variability in the atmosphere and ocean is needed. Companion field studies and numerical prediction experiments were discussed. Third, research prognostic models suggest that future operational forecast models pertaining to coastal weather will be high resolution and site specific, and will properly treat effects of local coastal geography, orography, and ocean state. The view was expressed that the exploration of coupled air-sea models of the coastal zone would be a particularly fruitful area of research. PDT-3 felt that forecasts of land-impacting tropical cyclones, Great Lakes-affected weather, and coastal cyclogenesis, in particular, would benefit from such coordinated modeling and field efforts. Fourth, forecasting for Arctic coastal zones is limited by our understanding of how sea ice forms. The importance of understanding air-sea fluxes and boundary layers in the presence of ice formation was discussed. Finally, coastal flash flood forecasting via hydrologic models is limited by the present accuracy of measured and predicted precipitation and storm surge events. Strategies for better ways to improve the latter were discussed.

Conference Paper

Water-quality characteristics and contaminants in the rural karst-dominated Spring Mill Lake watershed, southern Indiana

The Spring Mill Lake watershed is located in the Mitchell Plateau, a karst area that developed on Mississippian carbonates in southern Indiana. Spring Mill Lake is a reservoir built in the late 1930s and is located in Spring Mill State Park. Within the park, groundwater from subsurface conduits issues as natural springs and then flows in surface streams to the lake. From 1998 to 2002, surface and subsurface hydrology and water quality were investigated to determine the types and sources of potential contaminants entering the lake. Water samples collected during base flow and a February 2000 storm event were analyzed for selected cations, anions, trace elements, selected U.S. Environmental Protection Agency (EPA) primary and secondary drinkingwater contaminants, nitrogen isotopes, suspended solids, Escherichia coli, and pesticides. All of the water samples met the EPA drinking-water standards for inorganic constituents, except those collected at five sites in August 1999 during a drought. Nitrate nitrogen (NO 3 -N) concentrations were highest during base-flow conditions and displayed a dilutional trend during peak-flow periods. The NO 3 -N concentrations in water samples collected during the 2001 spring fertilizer applications tended to increase from early to late spring. All of the ?? 15 N values were low, which is indicative of either an inorganic source or soil organic matter. Storm discharge contained increased concentrations of total suspended solids; thus, storms are responsible for most of the sediment accumulation in the lake. E. coli levels in 24% of the samples analyzed contained a most probable number (MPN) greater than 235/100 mL, which is the maximum acceptable level set for recreational waters in Indiana. E. coli does appear to be a potential health risk, particularly at Rubble spring. The sources of E. coli found at this spring may include barnyard runoff from a horse barn or wastes from a wastewater treatment facility. The pesticides atrazine, metolachlor, acetochlor, and simazine were detected during the spring of 2001. Atrazine, metolachlor, acetochlor, and simazine are used to suppress weeds during corn and soybean production. Additional sources of atrazine and simazine may result from application to right-of-ways, orchards, and managed forest areas. ?? 2006 Geological Society of America.

Special Paper of the Geological Society of America

The Maryland Coastal Plain Aquifer Information System: A GIS-based tool for assessing groundwater resources

Groundwater is the source of drinking water for ∼1.4 million people in the Coastal Plain Province of Maryland (USA). In addition, groundwater is essential for commercial, industrial, and agricultural uses. Approximately 0.757 × 109 L d –1 (200 million gallons/d) were withdrawn in 2010. As a result of decades of withdrawals from the coastal plain confined aquifers, groundwater levels have declined by as much as 70 m (230 ft) from estimated prepumping levels. Other issues posing challenges to long-term groundwater sustainability include degraded water quality from both man-made and natural sources, reduced stream base flow, land subsidence, and changing recharge patterns (drought) caused by climate change. In Maryland, groundwater supply is managed primarily by the Maryland Department of the Environment, which seeks to balance reasonable use of the resource with long-term sustainability. The chief goal of groundwater management in Maryland is to ensure safe and adequate supplies for all current and future users through the implementation of appropriate usage, planning, and conservation policies. To assist in that effort, the geographic information system (GIS)–based Maryland Coastal Plain Aquifer Information System was developed as a tool to help water managers access and visualize groundwater data for use in the evaluation of groundwater allocation and use permits. The system, contained within an ESRI ArcMap desktop environment, includes both interpreted and basic data for 16 aquifers and 14 confining units. Data map layers include aquifer and ­confining unit layer surfaces, aquifer extents, borehole information, hydraulic properties, time-series groundwater-level data, well records, and geophysical and lithologic logs. The aquifer and confining unit layer surfaces were generated specifically for the GIS system. The system also contains select groundwater-quality data and map layers that quantify groundwater and surface-water withdrawals. The aquifer information system can serve as a pre- and postprocessing environment for groundwater-flow models for use in water-supply planning, development, and management. The system also can be expanded to include features that evaluate constraints to groundwater development, such as insufficient available drawdown, degraded groundwater quality, insufficient aquifer yields, and well-field interference. Ultimately, the aquifer information system is intended to function as an interactive Web-based utility that provides a broad array of information related to groundwater resources in Maryland’s coastal plain to a wide-ranging audience, including well drillers, consultants, academia, and the general public.

Special Paper of the Geological Society of America

Are the Columbia River Basalts, Columbia Plateau, Idaho, Oregon, and Washington, USA, a viable geothermal target? A preliminary analysis

The successful development of a geothermal electric power generation facility relies on (1) the identification of sufficiently high temperatures at an economically viable depth and (2) the existence of or potential to create and maintain a permeable zone (permeability >10 -14 m 2 ) of sufficient size to allow efficient long-term extraction of heat from the reservoir host rock. If both occur at depth under the Columbia Plateau, development of geothermal resources there has the potential to expand both the magnitude and spatial extent of geothermal energy production. However, a number of scientific and technical issues must be resolved in order to evaluate the likelihood that the Columbia River Basalts, or deeper geologic units under the Columbia Plateau, are viable geothermal targets. Recent research has demonstrated that heat flow beneath the Columbia Plateau Regional Aquifer System may be higher than previously measured in relatively shallow (<600 m depth) wells, indicating that sufficient temperatures for electricity generation occur at depths 5 km. The remaining consideration is evaluating the likelihood that naturally high permeability exists, or that it is possible to replicate the high average permeability (approximately 10 -14 to 10 -12 m 2 ) characteristic of natural hydrothermal reservoirs. From a hydraulic perspective, Columbia River Basalts are typically divided into dense, impermeable flow interiors and interflow zones comprising the top of one flow, the bottom of the overlying flow, and any sedimentary interbed. Interflow zones are highly variable in texture but, at depths <600 m, some of them form highly permeable regional aquifers with connectivity over many tens of kilometers. Below depths of ~600 m, permeability reduction occurs in many interflow zones, caused by the formation of low-temperature hydrothermal alteration minerals (corresponding to temperatures above ~35 °C). However, some high permeability (>10-14 m 2 ) interflows are documented at depths up to ~1,400 m. If the elevated permeability in these zones persists to greater depths, they may provide natural permeability of sufficient magnitude to allow their exploitation as conventional geothermal reservoirs. Alternatively, if the permeability in these interflow zones is less than 10-14 m 2 at depth, it may be possible to use hydraulic and thermal stimulation to enhance the permeability of both the interflow zones and the natural jointing within the low-permeability interior portions of individual basalt flows in order to develop Enhanced/Engineered Geothermal System (EGS) reservoirs. The key challenge for an improved Columbia Plateau geothermal assessment is acquiring and interpreting comprehensive field data that can provide quantitative constraints on the recovery of heat from the Columbia River Basalts at depths greater than those currently tested by deep boreholes.

Idaho, Oregon, Washington

Mapping ground water in three dimensions: An analysis of airborne geophysical surveys of the Upper San Pedro River Basin, Cochise County, southeastern Arizona

This report summarizes the results of two airborne geophysical surveys conducted in the upper San Pedro Valley of southeastern Arizona in 1997 and 1999. The combined surveys cover about 1,000 square kilometers and extend from the Huachuca Mountains on the west to the Mule Mountains and Tombstone Hills on the east and from north of the Babocomari River to near the Mexican border on the south. The surveys included the acquisition of high-resolution magnetic data, which were used to map depth to the crystalline basement rocks underlying the sediments filling the basin. The magnetic inversion results show a complex basement morphology, with sediment thickness in the center of the valley ranging from ~237 meters beneath the city of Sierra Vista to ~1,500 meters beneath Huachuca City and the Palominas area near the Mexican border. The surveys also included acquisition of 60-channel time-domain electromagnetic (EM) data. Extensive quality analyses of these data, including inversion to conductivity vs. depth (conductivity-depth-transform or CDT) profiles and comparisons with electrical well logs, show that the electrical conductor mapped represents the subsurface water-bearing sediments throughout most of the basin. In a few places (notably the mouth of Huachuca Canyon), the reported water table lies above where the electrical conductor places it. These exceptions appear to be due to a combination of outdated water-table information, significant horizontal displacement between the wells and the CDT profiles, and a subtle calibration issue with the CDT algorithm apparent only in areas of highly resistive (very dry) overburden. These occasional disparities appear in less than 5 percent of the surveyed area. Observations show, however, that wells drilled in the thick unsaturated zone along the Huachuca Mountain front eventually intersect water, at which point the water rapidly rises high into the unsaturated zone within the wellbore. This rising of water in a wellbore implies some sort of confinement below the thick unsaturated zone, a confinement that is not identified in the available literature. Occasional disparities notwithstanding, maps of the electrical conductor derived from the airborne EM system provide a synoptic view of the presence of water underlying the upper San Pedro Valley, including its three-dimensional distribution. The EM data even show faults previously only inferred from geologic mapping. The magnetic and electromagnetic data together appear to show the thickness of the sediments, the water in the saturated sediments down to a maximum of about 400 meters depth, and even places where the main ground-water body is not in direct contact with the San Pedro River. However, the geophysical data cannot reveal anything directly about hydraulic conductivity or ground-water flow. Estimating these characteristics requires new hydraulic modeling based in part on this report. One concern to reviewers of this report is the effect that clays may have on the electrical conductor mapped with the airborne geophysical system. Although the water in the basin is unusually conductive, averaging 338 microsiemens per centimeter, reasoning cited below suggests that the contribution of clays to the overall conductivity would be relatively small. Basic principles of sedimentary geology suggest that silts and clays should dominate the center of the basin, while sands and gravels would tend to dominate the margins. Although clay content may increase the amplitude of the observed electrical conductors somewhat, it will not affect the depths to the conductor derived from depth inversions. Further, fine-grained sediments generally have higher porosity and tend to lie toward a basin center, a fact in general agreement with the observed geophysical data.

Arizona

Assessing streamflow sensitivity to variations in glacier mass balance

The mountains ringing the Gulf of Alaska (GOA) receive upwards of 4–8 m yr −1 of precipitation (Simpson et al.2005; Weingartner et al. 2005; O’Neel 2012), much of which runs off into productive coastal waters. The alpine landscape is heavily glacierized, and storage and turnover of water by glaciers substantially influences the regional surface water balance (Neal et al. 2010). In turn, the land-to-ocean flux of freshwater impacts the biogeochemistry, physical oceanography, freshwater and marine ecology of the downstream components of the GOA ecosystem (e.g., Royer et al. 2001; Hood and Scott 2008). In this way, the links between terrestrial and ocean ecosystems along the GOA have widespread impacts on regional socioeconomic issues including water and hydropower resources, fish populations, and sea level change (Dorava and Milner 2000; Royer and Grosch 2006; Cherry et al. 2010; Gardner et al. 2013). Moreover, predicting future changes in physical, chemical and biological processes in near-shore ecosystems along the GOA hinges, in part, on developing a robust understanding of water storage and transfer by glaciers through streams to the ocean. Glacierized basins (i.e. presently ice covered as opposed to glaciated, or historically ice covered) are very efficient producers of runoff, yielding 2–10 times greater runoff than similarly sized, non-glacierized basins (Mayo 1984). The unique energy balance that characterizes these basins (Jansson et al. 2003; Hock 2005) results in substantial alterations to streamflow, even when fractional ice coverage is very small (Stahl and Moore 2006). Consistent and precise treatment of glacier runoff is key to accurate assessment of hydrologic, ecological and socioeconomic impacts, but previously used definitions for glacier runoff are variable. They include: 1) meltwater produced as negative annual mass balance (e.g., Fountain and Tangborn 1985); 2) storage changes in the monthly water budget, where solid precipitation is balanced by melt and evaporation (Huss 2011, concept #2); 3) meltwater derived from melting ice only (irrespective of melting snow or mass balance) (Nolin et al. 2010; Huss 2011, concept #1); 4) all meltwater derived from the glacier surface (Cogley et al. 2011, meltwater runoff); 5) total runoff from the glacier surface (meltwater runoff plus rain on the glacier) (Neal et al. 2010). Total glacier runoff (Definitions 4 and 5 above) includes a contribution from annual mass balance, i.e. the sum of accumulation and ablation through a mass balance year (Definition 1), or what has historically been referred to as the “net” balance (Cogley et al. 2011). Indeed, annual balance has been shown to be an important driver of streamflow trends in glacierized basins, with periods of persistent negative annual balance resulting in statistically significant increases in streamflow (e.g., Pellicciotti et al. 2010). However, in maritime climates, anomalies in glacier runoff can be disconnected from annual balance because of the high variability in winter precipitation. For example, positive anomalies in winter accumulation can result in elevated levels of glacier runoff in times of positive annual mass balance (Thayyen and Gergan 2010). Quantifying the impacts of changing glacier geometries (annual balance) on glacier runoff is essential for predicting future changes in streamflow in glacierized basins. However, determining the role that this component plays in total glacier runoff (Definition 5) requires consistent measurements of seasonal (or shorter period) mass balances, measurements of precipitation at multiple locations within a basin, and streamflow measurements in close proximity to a glacier’s terminus. Practical and logistical challenges associated with assembling such data sets typically preclude such partitioning. As a result, most analyses of the relationship between annual mass balance and streamflow rely on some component of model output to compute glacier runoff (e.g. Huss et al. 2008; Kaser et al. 2010). Ultimately, developing an understanding of how total glacier runoff will change in the future is critical for predicting downstream ecological impacts associated with changes in riverine fluxes of water, sediment, and solutes (e.g., metals and nutrients) to near-shore coastal ecosystems. The purpose of this paper is to evaluate relationships among seasonal and annual glacier mass balances, glacier runoff and streamflow in two glacierized basins in different climate settings. We use long-term glacier mass balance and streamflow datasets from the United States Geological Survey (USGS) Alaska Benchmark Glacier Program to compare and contrast glacier-streamflow interactions in a maritime climate (Wolverine Glacier) with those in a continental climate (Gulkana Glacier). Our overall goal is to improve our understanding of how glacier mass balance processes impact streamflow, ultimately improving our conceptual understanding of the future evolution of glacier runoff in continental and maritime climates.

Alaska

Hydrology of the Price River basin, Utah, with emphasis on selected coal-field areas

Data obtained during a hydrologic study of the Price River basin, Utah, are used to describe seasonal variations off low of springs, relation between ground water and surface water, hydraulic properties of the ground-water reservoir, ground-water recharge and discharge, flood characteristics of streams, mineralogic composition and depositional rates of sediments, nutrient and inorganic loading in streams and Scofield Reservoir, and water budgets for selected basins. Additional study and monitoring are needed to detect possible hydrologic changes caused by coal mining. Much of the ground-water discharge from the Star Point Sandstone in the Mesaverde Group in the Wasatch Plateau occurs along faults. In the Book Cliffs, where faulting is less extensive, most of the ground-water discharge is from the Flagstaff Limestone. The Flagstaff Limestone is greatly diffusive, has a small storage coefficient, and contains water which is perched. Springs issuing from the Star Point Sandstone in the Mud Creek drainage (Wasatch Plateau) had recession indexes greater than 365 days per log cycle. Springs issuing at higher altitudes from the Colton Formation and the Flagstaff Limestone in the Soldier Creek area (Book Cliffs) have great seasonal variability, with recession indexes ranging from 24 to 115 days per log cycle. Estimated transmissivities in the Soldier Creek area ranged from 0.003 foot squared per day in the lower part of the Castlegate Sandstone to 0.07 foot squared per day in the Price River Formation. Seepage from the Star Point Sandstone is the major contributor to base flow of the stream in Eccles Canyon (Wasatch Plateau). Gains of as much as 230 gallons per minute occurred near a fault zone which crosses Eccles Canyon at the junction with South Fork Canyon. The potentiometric surface of water in the Blackhawk Formation in the Wasatch Plateau (Mud Creek drainage) and the Book Cliffs (Soldier Creek area) generally is above the coal zones, and dewatering will be necessary in most places during mining. The total reported mine discharge during 1980 was about 2,800 acre-feet. During the 1980 water year, ground-water discharge from the Flagstaff Limestone contributed about 4,400 acre-feet, or 41 percent of the combined streamflow of Coal, Soldier, and Dugout Creeks. Discharge from the Blackhawk and other formations contributed about 2,200 acrefeet, or 21 percent of the streamflow. The recharge to a 4,020-acre outcrop of the Flagstaff Limestone in the Soldier Creek area was determined from measurements of the discharge of springs and seepage gains in streams. The estimated recharge was 2 inches per year, or 12 percent of the April 1 snowpack, or 9 percent of the precipitation during 1980. Benthic-diversity indexes at streams in the Book Cliffs and Wasatch Plateau are related to the size of the bed material. The Shannon-Weiner diversity index for five sites in the Mud Creek drainage ranged from 1.57 to 4.'20 during four sampling periods in 1979-80. In the Soldier Creek drainage, the diversity index for four sampling sites ranged from 0.8 to 2.3 during three sampling periods in 1979-80. Sediment deposition in Scofield Reservoir during 1943-79 was estimated to be 3,000 acre-feet but most of the coal deposition in the reservoir occurred prior to about 1950. An anaerobic zone develops in the reservoir during the summer, but water in the reservoir completely mixes during late summer or early fall.

Utah

Connectivity of Mojave Desert tortoise populations—Management implications for maintaining a viable recovery network

Executive Summary The historic distribution of Mojave desert tortoises ( Gopherus agassizii ) was relatively continuous across the range, and the importance of tortoise habitat outside of designated tortoise conservation areas (TCAs) to recovery has long been recognized for its contributions to supporting gene flow between TCAs and to minimizing impacts and edge effects within TCAs. However, connectivity of Mojave desert tortoise populations has become a concern because of recent and proposed development of large tracts of desert tortoise habitat that cross, fragment, and surround designated conservation areas. This paper summarizes the underlying concepts and importance of connectivity for Mojave desert tortoise populations by reviewing current information on connectivity and providing information to managers for maintaining or enhancing desert tortoise population connectivity as they consider future proposals for development and management actions. Maintaining an ecological network for the Mojave desert tortoise, with a system of core habitats (TCAs) connected by linkages, is necessary to support demographically viable populations and long-term gene flow within and between TCAs. There are four points for wildlife and land-management agencies to consider when making decisions that could affect connectivity of Mojave desert tortoise populations (for example, in updating actions in resource management plans or amendments that could help maintain or restore functional connectivity in light of the latest information): Management of all desert tortoise habitat for persistence and connectivity . Desert tortoise populations continue to decline within most TCAs, and it is unlikely that trends are better in populations outside protected areas. Fragmentation exacerbates negative population trends by breaking large continuous populations into smaller isolated populations. Connectivity within large populations can enhance resilience to localized disturbances due to rescue by neighboring individuals. In contrast, smaller fragmented populations are resistant to rescue by their isolation and thus could suffer irreversible declines to extirpation from a variety of threats and stochastic events. Enhanced threat reduction to reverse declines within TCAs and to maintain occupied habitat in the surrounding matrix would help reduce the variability in population growth rates and improve the resilience of protected populations even while implementing efforts to improve connectivity. Each TCA has unique strengths and weaknesses regarding its ability to support minimum sustainable populations based on areal extent and its ability to support population increases based on landscape connection with adjacent populations. Considering how proposed projects (inside or outside of TCAs) affect connectivity and the ability of TCAs to support at least 5,000 adult tortoises (the numerical goal for each TCA) could help managers to maintain the resilience of TCAs to population declines. The same project, in an alternative location, could have very different impacts on local and regional populations. For example, within the habitat matrix surrounding TCAs, narrowly delineated corridors may not allow for natural population dynamics if they do not accommodate overlapping home ranges along most of their widths so that tortoises reside, grow, find mates, and produce offspring that can replace older tortoises. In addition, most habitat outside TCAs may receive more surface disturbance than habitat within TCAs. Therefore, managing the entire remaining matrix of desert tortoise habitat for permeability may be better than delineating fixed corridors. These concepts apply, especially given uncertainty about long-term condition of habitat, within and outside of TCAs under a changing climate. Ultimately, questions such as “ What are the critical linkages that need to be protected ?” could be better framed as “ How can we manage the remaining habitat matrix in ways that sustain ecological processes and habitat suitability for special status species ?” Land-management decisions made in the context of the latter question may be more conducive to maintenance of a functional ecological network. Limitations on landscape-level disturbance across habitat managed for the desert tortoise Clearly delineating habitat linkages and differentiating them from non-delineated areas by the uses that are permitted or prohibited within them by specific management guidelines can help achieve functional connectivity. Such guidelines would be most effective if they considered and accounted for all surface disturbances (for example, temporary disturbances such as fiberoptic lines or off-highway vehicle routes, right-of-ways, utility-scale solar development, urbanization) to the extent possible. A weighted framework that varies with the permanence or severity of the disturbance, and can be additive to quantify cumulative effects, could be useful (Xiong, 2020). For example, minor roads can alter tortoise movements independently of other features (Peaden and others, 2017; Hromada and others, 2020), but if the isolated dirt road is accompanied by a powerline that encourages raven predation (Xiong, 2020), then the two features together may be additive. Ignoring minor or temporary disturbance on the landscape could result in a cumulatively large impact that is not explicitly acknowledged (Goble, 2009); therefore, understanding and quantifying all surface disturbance on a given landscape is prudent. In California, the Bureau of Land Management established 0.1–1.0 percent caps on new surface-disturbance for TCAs and mapped linkages that address the issues described in number 1 of this list. Nevada, Utah, and Arizona currently do not have surface-disturbance limits. Limits comparable to those in the Desert Renewable Energy Conservation Plan (DRECP) would be 0.5 percent within TCAs and 1 percent within the linkages modeled by Averill-Murray and others (2013). Limits in some areas of California within the Desert Renewable Energy Conservation Plan, such as Ivanpah Valley, are more restrictive, at 0.1 percent. Continuity across the state line in Nevada could be achieved with comparable limits in the adjacent portion of Ivanpah Valley, as well as the Greater Trout Canyon Translocation Area and the Stump Springs Regional Augmentation Site. These more restrictive limits would help protect remaining habitat in the major interstate connectivity pathway through Ivanpah Valley and focal areas of population augmentation that provide additional population connectivity along the western flank of the Spring Mountains. In a recent study that analyzed 13 years of desert tortoise monitoring data, nearly all desert tortoise observations were at sites in which 5 percent or less of the surrounding landscape within 1 kilometer was disturbed (Carter and others, 2020a). To help maintain tortoise habitability and permeability across all other non-conservation-designated tortoise habitat, all surface disturbance could be limited to less than 5-percent development per square kilometer because the 5-percent threshold for development is the point at which tortoise occupation drops precipitously (Carter and others, 2020a). However, although individual desert tortoises were observed at development levels up to 5 percent, we do not know the fitness or reproductive characteristics of these individuals. This level of development also may not allow for long-term persistence of healthy populations that are of adequate size needed for demographic or functional connectivity; therefore, a conservative interpretation suggests that, ideally, development could be lower. Lower development levels would be particularly useful in areas within the upper 5th percentile of connectivity values modeled by Gray and others (2019). Reducing ancillary threats in places where connectivity is restricted to narrow strips of habitat, for example, narrow mountain passes or vegetated strips between solar development, could enhance the functionality of these vulnerable linkages. In such areas, maintaining multiple, redundant linkages could further enhance overall connectivity. Minimization of mortality from roads and maximization of passage under roads . Roads pose a significant threat to the long-term persistence of local tortoise populations, and roads of high traffic volume lead to severe population declines, which ultimately fragments populations farther away from the roads. Three points (a.–c.) pertain to reducing direct mortality of tortoises on the many paved roads that cross desert tortoise habitat and to maintaining a minimal level of permeability across these roads: Tortoise-exclusion fencing tied into culverts, underpasses, overpasses, or other passages below roads in desert tortoise habitat, would limit vehicular mortality of tortoises and provide opportunities for movement across the roads. Installation of shade structures on the habitat side of fences installed in areas with narrow population-depletion zones would limit overheating of tortoises that may pace the fence. Passages below highways could be maintained or retrofitted to ensure safe tortoise access, for example, by filling eroded drop-offs or modifying erosion-control features such as rip-rap to make them safer and more passable for tortoises. Wildlife management agencies could work with transportation departments to develop construction standards that are consistent with hydrologic/erosion management goals, while also incorporating a design and materials consistent with tortoise survival and passage and make the standards widely available. The process would be most effective if the status of passages was regularly monitored and built into management plans. Healthy tortoise populations along fenced highways could be supported by ensuring that land inside tortoise-exclusion fences is not so degraded that it leads to degradation of tortoise habitat outside the exclusion areas. For example, severe invasive plant infestations inside a highway exclusion could cause an increase of invasive plants outside the exclusion area and degrade habitat; therefore, invasive plants inside road rights of way could be mown or treated with herbicide to limit their spread into adjacent tortoise habitat and minimize the risk of these plants carrying wildfires into adjacent habitat. Adaptation of management based on new information . Future research will continue to build upon and refine models related to desert tortoise population connectivity and develop new ones. New models could consider landscape levels of development and be constructed such that they share common foundations to support future synthesis efforts. If model development was undertaken in partnership with entities that are responsible for management of desert tortoise habitat, it would facilitate incorporation of current and future modeling results into their land management decisions. There are specific topics that may be clarified with further evaluation: The effects of climate change on desert tortoise habitat, distribution, and population connectivity; The effects of large-scale fires, especially within repeatedly burned habitat, on desert tortoise distribution and population connectivity; The ability of solar energy facilities or similar developments to support tortoise movement and presence by leaving washes intact; leaving native vegetation intact whenever possible, or if not possible, mowing the site, allowing vegetation to re-sprout, and managing weeds; and allowing tortoises to occupy the sites; and The design and frequency of underpasses necessary to maintain functional demographic and genetic connectivity across linear features, like highways.

Arizona, California, Nevada

City-scale geothermal energy everywhere to support renewable resilience – A transcontinental cooperation

Cities have important and varying incentives to transform their energy sector to all-electric with low carbon emissions. However, they often encounter a number of impediments when attempting to implement such a change. For example, while urban areas have the highest energy demand-density, cities often lack the space for installing additional energy generation and/or long-duration energy storage systems. Cities also have existing environmental issues from energy sources (e.g., pollution from dust, waste heat or noise) that make residents sensitive to energy infrastructure development. Utilizing power from conventional sources, such as natural gas, biomass and hydropower, which usually are distanced from the urban areas, also make cities more vulnerable to supply disruptions. One promising de-carbonizing energy option for cities focuses on their heating and cooling needs, which constitutes around 1/3 of U.S. and 1/2 of European energy consumption (including industrial processes like drying, pasteurization, etc.; Jadun and others, 2017; EU Commission 2022). If heating and cooling loads can be met by geothermal direct-use technologies, then the need for new electric sources can be greatly lessened. Despite the proven efficacy of geothermal energy as a city/community-scale heating and cooling resource, it is currently only a niche resource in the heating and cooling sector, though has significant potential for future growth. Historically, emphasis has been placed on geothermal electricity generation potential that requires higher temperature (greater than 90 °C) resources at drillable depths, but potentially viable areas are geographically limited and typically well removed from urban centers. Key drivers for investments were represented by greater political interest in renewable electricity production, higher revenues and less effort in distributing the produced energy via grids. In contrast, low-temperature (less than 90 °C) geothermal resources can be used directly for heating and cooling almost everywhere and are cost-effective in urban/suburban settings. In addition, the increased prominence of renewable electricity sources, such as wind and solar onto city-scale electric grids, has led to new urgency around questions of energy storage. Underground thermal energy storage (UTES), wherein surplus or waste heat is stored underground for later use, could present a long-duration energy storage solution. From October 2022 through September 2024, a transcontinental consortium consisting of geological surveys, geoscience organizations, industry representatives and universities aims to develop an understanding of the global potential for city-scale geothermal, proposing guidelines to aid in promoting the economic utilization of low temperature geothermal resources. Efforts will focus on providing city managers and other decision makers with the information needed to evaluate and implement suitable city/community-scale geothermal technologies. Funded by the U.S. Geological Survey’s John Wesley Powell Center for Analysis and Synthesis, this interdisciplinary consortium will showcase tools, datasets, and scientific recommendations to accelerate the broader understanding and adoption of renewable energy systems that access geothermal resources. The collaborative research activities include standardization of nomenclatures, resource description and characterization strategies globally. The results from these activities will be combined with a preliminary climate-driven, city-based energy needs related analysis to perform energy supply/demand matching analysis. The identification of city-specific applications that would benefit from the geothermal technologies provides the basis to up-scale city-specific determinations to regional and national assessments of resource estimates. The city-scale geothermal energy research initiative will ultimately provide the synergies and management analysis that can address benefits, environmental impacts, regulatory frameworks, sustainability, and suitability in retrofitted buildings or new as well as existing heating networks.

Conference Paper

Forest conditions in the San Francisco Mountains Forest Reserve, Arizona

On April 12, 1902, President Roosevelt issued a proclamation "for the purpose of consolidating into one reserve the lands heretofore embraced in the San Francisco Mountains Forest Reserves and of including therein the other adjacent lands within Francisco Mountains Forest Reserve, and is described by metes and bound as follows: ''Beginning at the northwest corner of township twenty-two (22) north, range one (1) east, Gila and Salt River meridian, Arizona; thence southerly along the said meridian, allowing for the proper offset on the fifth (5th) standard parallel north, to the southwest corner of Township nineteen (19) south, range one (1) east; thence easterly along the surveyed and unsurveyed township line to the point for the northwest corner of township eighteen (18) north, range four (4) east; thence southerly along the unsurveyed range line to its intersection with the fourth (4th) standard' parallel north; thence easterly along said parallel to the point for the northwest corner of township sixteen (16) north, range five (5) east; thence southerly to the point for the southwest corner of said township; then easterly to the point for the northwest corne1&middot; of township fifteen (15) north, range six (6) east; thence southerly to the point for the southwest corner of section eighteen (18), said township; thence easterly along the unsurveyed section line to the point for the northwest corner of section nineteen (19), township fifteen (15) north, range seven (7) east; thence southerly to the southwest corner of said section; thence easterly along the unsurveyed section lines to the southwest corner of section nineteen (19), township fifteen (15) north, range nine (9) east; thence northerly to the northwest corner of said section; thence easterly along the section line to the southeast corner of section thirteen (13), said township; thence northerly to the northeast corner of section twelve (12), said township; thence easterly along the section lines to the southeast corner of section one (1), township fifteen (15) north, range eleven (11) east; thence northerly along the range line to its intersection with the fourth (4th) standard parallel north; thence westerly along said parallel to the southeast corner of township seventeen (17) north, range eleven (11) east; thence northerly along the surveyed and unsurveyed range line to the point for the northeast corner of township eighteen (18) north, range eleven (11) east; thence westerly to the southeast corner of township nineteen (19) north, range ten (10) east; thence northerly along the range line to its intersection with the fifth (5th) standard parallel north; thence westerly along said parallel to the point for the southeast corner of township twenty-one (21) north, range nine (9) east; thence northerly along the unsurveyed range line, allowing for the proper offset on the sixth (6th) standard parallel north, to the point for the northeast cornet&middot; of township twenty-five (25) north, range nine (9) east; thence westerly along the surveyed and unsurveyed township line to the point for the northwest corner of township twenty-five (25) north, range three (3) east; thence southerly along the surveyed and unsurveyed range line, allowing for the proper offset on the sixth (6th) standard parallel north, to the northeast corner of township twenty-two (22) north, range two (2) east; thence westerly along the township line to the northwest corner of township twenty-two (22) north, range one (1) east, to the place of beginning."

Arizona

Comparing empirical sediment transport modeling approaches in Michigan rivers

Excess or limited fluvial sediment transport can contribute to and exacerbate many environmental issues including nutrient loading, aquatic habitat degradation, flooding, channel navigation dredging, dam operation, and stream degradation or aggradation. However, fluvial sediment transport is difficult and expensive to comprehensively characterize because it can vary substantially both temporally and spatially. Having better estimates of fluvial sediment transport is important for understanding and solving these environmental issues when it is not possible to collect fluvial sediment samples. Different modeling approaches can be used to help estimate suspended sediment when sampling data are limited or unavailable. This study compared dimensionless sediment rating curves (DSRCs) developed in Pagosa Springs Colorado, Minnesota, and Michigan to determine if these DSRCs were suitable to make predictions of suspended sediment for Michigan rivers. Approximately 3,000 suspended sediment samples collected in or near Michigan from the mid-1960s through August 2022 were used to develop two DSRC models. The DSRCs developed in Michigan include a pooled DSRC model which uses nonlinear least squares regression, and a mixed-effects DSRC model which uses a mixed-effects modeling approach. In general, there was not a noticeable improvement in the performance of the Michigan mixed-effects DSRC model over the Michigan pooled DSRC model. The two Michigan DSRCs were evaluated against DSRCs developed for Pagosa Springs and Minnesota. The results showed DSRC models developed from Minnesota and Michigan were similar to each other. In contrast, the Pagosa Springs DSRC predicts higher suspended-sediment concentration (SSC) at low flows and increases at a higher rate due to having a greater exponent. The Pagosa Springs DSRC produces higher SSC predictions that do not approximate the observed data well at most of the Michigan sites in the study. The results suggest that the Pagosa Springs DSRC was not suitable to make predictions of suspended sediment for Michigan rivers. The similarity of the DSRC equations developed for Minnesota and Michigan compared to the Pagosa Springs DSRC equation suggest that there may be regional patterns of SSC in the upper Midwest rivers that differ from those in other areas of the country like Pagosa Springs. A regionally applicable model could be developed and strengthened by combining data from additional midwestern states. Since the Michigan DSRCs goodness-of-fit metrics were comparable to the site-specific simple linear regressions (SLRs) and outperformed them in the aggregate goodness-of-fit metrics, the Michigan DSRCs are suitable to make predictions of suspended sediment in Michigan rivers with limited data. However, the availability of the DSRCs from this study should not diminish the value of collecting physical samples and exploring alternative modeling approaches because of the uncertainty associated with using DSRCs. Approximately 3,000 suspended sediment samples collected in or near Michigan from the mid-1960s through August 2022 were used to develop two DSRC models. The DSRCs developed in Michigan include a pooled DSRC model which uses nonlinear least squares regression, and a mixed-effects DSRC model which uses a mixed-effects modeling approach. In general, there was not a noticeable improvement in the performance of the Michigan mixed-effects DSRC model over the Michigan pooled DSRC model. The two Michigan DSRCs were evaluated against DSRCs developed for Pagosa Springs and Minnesota. The results showed DSRC models developed from Minnesota and Michigan were similar to each other. In contrast, the Pagosa Springs DSRC predicts higher suspended-sediment concentration (SSC) at low flows and increases at a higher rate due to having a greater exponent. The Pagosa Springs DSRC produces higher SSC predictions that do not approximate the observed data well at most of the Michigan sites in the study. The results suggest that the Pagosa Springs DSRC was not suitable to make predictions of suspended sediment for Michigan rivers. The similarity of the DSRC equations developed for Minnesota and Michigan compared to the Pagosa Springs DSRC equation suggest that there may be regional patterns of SSC in the upper Midwest rivers that differ from those in other areas of the country like Pagosa Springs. A regionally applicable model could be developed and strengthened by combining data from additional midwestern states. Since the Michigan DSRCs goodness-of-fit metrics were comparable to the site-specific simple linear regressions (SLRs) and outperformed them in the aggregate goodness-of-fit metrics, the Michigan DSRCs are suitable to make predictions of suspended sediment in Michigan rivers with limited data. However, the availability of the DSRCs from this study should not diminish the value of collecting physical samples and exploring alternative modeling approaches because of the uncertainty associated with using DSRCs.

Michigan

Volcanic processes and geology of Augustine Volcano, Alaska

Augustine Island (volcano) in lower Cook Inlet, Alaska, has erupted repeatedly in late-Holocene and historical times. Eruptions typically beget high-energy volcanic processes. Most notable are bouldery debris avalanches containing immense angular clasts shed from summit domes. Coarse deposits of these avalanches form much of Augustine's lower flanks. A new geologic map at 1:25,000 scale depicts these deposits, these processes. We correlate deposits by tephra layers calibrated by many radiocarbon dates.Augustine Volcano began erupting on the flank of a small island of Jurassic clastic-sedimentary rock before the late Wisconsin glaciation (late Pleistocene). The oldest known effusions ranged from olivine basalt explosively propelled by steam, to highly explosive magmatic eruptions of dacite or rhyodacite shed as pumice flows. Late Wisconsin piedmont glaciers issuing from the mountainous western mainland surrounded the island while dacitic eruptive debris swept down the south volcano flank.Evidence is scant for eruptions between the late Wisconsin and about 2,200 yr B.P. On a few south-flank inliers, thick stratigraphically low pumiceous pyroclastic-flow and fall deposits probably represent this period from which we have no radiocarbon dates on Augustine Island. Eruptions between about 5,350 and 2,200 yr B.P. we know with certainty by distal tephras. On Shuyak Island 100 km southeast of Augustine, two distal fall ashes of Augustinian chemical provenance (microprobe analysis of glass) date respectively between about 5,330 and 5,020 yr B.P. and between about 3,620 and 3,360 yr B.P. An Augustine ash along Kamishak Creek 70 km southwest of Augustine dates between about 3,850 and 3,660 yr B.P. A probably Augustinian ash lying within peat near Homer dates to about 2,275 yr B.P.From before 2,200 yr B.P. to the present, Augustine eruptive products abundantly mantle the island. During this period, numerous coarse debris avalanches swept beyond Augustine's coast, most recently in A.D. 1883. The decapitated summit after the 1883 eruption, replaced by andesite domes of six eruptions since, shows a general process: collapse of steep summit domes, then the summit regrown by later dome eruptions. The island's stratigraphy is based on six or seven coarse-pumice tephra "marker beds." In upward succession they are layers G (2,100 yr B.P.), I (1,700 yr B.P.), H (1,400 yr B.P.), C (1,200-1,000 yr B.P.), M (750 yr B.P.), and B (390 yr B.P.).A coarse, hummocky debris-avalanche deposit older than about 2,100 yr B.P.-or perhaps a stack of three of them-lies along the east coast, the oldest exposed such bouldery diamicts on Augustine Island. Two large debris avalanches swept east and southeast into the sea between about 2,100 and 1,800 yr B.P. A large debris avalanche shed east and east-northeast into the sea between 1,700 and 14,00 yr B.P.Between about 1,400 and 1,100 yr B.P. debris avalanches swept into the sea on the volcano's south, southwest, and north-northwest. Pumiceous pyroclastic fans spread to the southeast and southwest, lithic pyroclastic flows and lahars (?) to the south and southeast. Pyroclastic flows, pyroclastic surges, and lahars swept down the west and south flanks between about 1,000 and 750 yr B.P.A debris avalanche swept into the sea on the west, and a small one on the south-southeast, between about 750 and 400 yr B.P. Large lithic pyroclastic flows shed to the southeast; smaller ones descended existing swales on the southwest and south.Between about 400 yr B.P. and historical time (late 1770s), three debris avalanches swept into the sea on the west-northwest, north-northwest, and north flanks. One of them (West Island) was large and fast: most of it rode to sea far beyond a former sea cliff, and its surface includes geomorphic evidence of having initiating a tsunami. Augustine's only conspicuous lava flow erupted on the north flank.During this prehistoric period numerous domes grew at the volcano's summit, remnants of which form the east and south sides of the present summit-dome complex. Three domes grew below the summit area on the upper south and northwest flanks. In between large eruptions that deposited coarse pumiceous fall beds, many smaller eruptions emplaced beds of sand-sized ash on the volcano flanks.During the past 750 years, beach and back-beach eolian dunes accreted at the southwest coast, forming a ribbed coastwise topography. Lesser dunes grew at the backs of beaches in coves on other flanks.An eruption in 1883 shed a debris avalanche swiftly into the sea on the north-northeast, followed by pyroclastic flows and surges. Eruptions in 1935 and 1963-64 grew summit domes that spilled over the southwest and south flanks and shed coarse rubbly lithic pyroclastic flows down those flanks. Eruptions and 1976 and 1986 grew domes that draped down the north flank and shed voluminous pyroclastic flows to the northeast through north-northwest flanks, when smaller pyroclastic flows and (or) lahars swept down other flanks. A small dome-building eruption in January-March 2006 after this report was all but complete we treat only fleetingly. The largest debris avalanches sweep into the sea at Augustine's coast at speeds inferred between 60 and 80 m/s. Augustine is capable of initiating damaging tsunami to lower Cook Inlet, but geologic evidence for them on the mainland is sporadic and sparse.

Alaska

Groundwater availability of the Mississippi embayment

Groundwater is an important resource for agricultural and municipal uses in the Mississippi embayment. Arkansas ranks first in the Nation for rice and third for cotton production, with both crops dependent on groundwater as a major source of irrigation requirements. Multiple municipalities rely on the groundwater resources to provide water for industrial and public use, which includes the city of Memphis, Tennessee. The demand for the groundwater resource has resulted in groundwater availability issues in the Mississippi embayment including: (1) declining groundwater levels of 50 feet or more in the Mississippi River Valley alluvial aquifer in parts of eastern Arkansas from agricultural pumping, (2) declining groundwater levels of over 360 feet over the last 90 years in the confined middle Claiborne aquifer in southern Arkansas and northern Louisiana from municipal pumping, and (3) litigation between the State of Mississippi and a Memphis water utility over water rights in the middle Claiborne aquifer. To provide information to stakeholders addressing the groundwater-availability issues, the U.S. Geological Survey Groundwater Resources Program supported a detailed assessment of groundwater availability through the Mississippi Embayment Regional Aquifer Study (MERAS). This assessment included (1) an evaluation of how these resources have changed over time through the use of groundwater budgets, (2) development of a numerical modeling tool to assess system responses to stresses from future human uses and climate trends, and (3) application of statistical tools to evaluate the importance of individual observations within a groundwater-monitoring network. An estimated 12 million acre-feet per year (11 billion gallons per day) of groundwater was pumped in 2005 from aquifers in the Mississippi embayment. Irrigation constitutes the largest groundwater use, accounting for approximately 10 million acre-feet per year (9 billion gallons per day) in 2000 from the Mississippi River Valley alluvial aquifer in Arkansas, Louisiana, Mississippi, and Missouri, and to a lesser extent in Illinois, Kentucky, and Tennessee. Predevelopment groundwater flow is represented in the MERAS model as a steady-state stress period, assumed to be prior to 1870. The simulated groundwater-flow budget indicates the largest predevelopment inflow to the system is net recharge to the alluvial aquifer. This inflow is balanced by outflow to gaining streams. Overall, water enters as net recharge to the alluvial aquifer or through outcrop areas of the various hydrogeologic units. Away from the outcrop areas, groundwater flow in the deeper formations is primarily upward into overlying units, ultimately discharging to streams through the alluvial aquifer. Total net recharge and discharge (sum of inflows or outflows) for the model ranged from about 0.66 million acre-feet per year during predevelopment to 20.16 million acre-feet per year by the end of the simulation (final simulated irrigation period in summer of 2006). This change in the model budget reflects increases in withdrawals compared to predevelopment conditions. Cumulative storage within aquifers simulated in the MERAS model indicates overall depletion of 140 million acre-feet (equivalent to 2.8 feet of water covering the entire study area). Postdevelopment inflow to the system is still through net recharge to the alluvial aquifer and the outcrop areas of the several hydrogeologic units, however, the flow between each unit is no longer upward to the alluvial aquifer. Groundwater flow during the summer of 2006 was primarily downward to offset demand from pumping. Early in the model simulation (1870-1920s), the primary components of the water budget were simulated as outflow from stream leakage and inflow from net recharge. As pumpage increased through time, water that would otherwise flow to streams reversed, and net stream leakage became an inflow to the system. The largest reversals began in the mid-1980s, but indications of the reversal began in the early 1960s with a trend in loss of streamflow leakage coupled with the first consistent inflow from storage. While groundwater pumped out of the alluvial aquifer was derived primarily from storage, pumpage out of the middle Claiborne aquifer was derived primarily from other aquifers (up to 15 percent from the alluvial aquifer), followed by flow from storage and net recharge. The potential consequences of climate change have been identified as a major concern facing the sustainability of the Nation's groundwater resources. To address this concern, two climate simulations were developed through the use of the MERAS model by extending the simulation period by 30 years to the year 2038 using extrapolated precipitation based on frequency analysis of historic climate cycles. There is little difference between the dry and wet scenarios in terms of percent water-level change. Both scenarios resulted in 14.6 to 13.9 percent of the area containing more than 100 feet of decline, 14.5 to 13.8 percent containing between 75 and 100 feet of decline, and 15.8 to 15.7 percent containing 51 to 75 feet of decline in the alluvial aquifer. The middle Claiborne aquifer water-level changes also were similar between the two scenarios. These scenarios indicate that even with a 25-percent increase in precipitation from that of the dry scenario, there is little difference in the resultant water levels. This is in large part because of the magnitude of differences between changes in net recharge and changes in pumping. When compared to the volume of water pumped out of the system, the effect of this change in net recharge is negligible. The groundwater-level monitoring network used to construct the 2007 middle Claiborne aquifer potentiometric surface was used as an example case to demonstrate statistical technique and to evaluate the importance of individual groundwater-level observations. To calculate the importance of each water-level observation to a prediction, predictions were specified as water-level altitudes near the end of the dry scenario simulation. These predictions were located near the center of cones of depression. Many of the observations that have a high importance are in close proximity to stressed areas of the aquifer.

Mississippi

Vegetation changes over 12 years in ungrazed and grazed Conservation Reserve Program Grasslands in the central and southern plains

The Conservation Reserve Program (CRP) established under the 1985 Food Security Act has the fundamental objectives of jointly providing economic support to segments of the agricultural community and conservation of natural resources (Osborn, 1997; Heard and others, 2000). Although soil loss on highly erodable lands was the principal natural resource conservation issue addressed in the 1985 CRP, improving water quality and wildlife habitat both became important considerations as the program evolved (Farmer and others, 1988). For example, Best and others (1997) found that production of young birds on CRP fields in the Midwest was ≥15 times the production on row-crop fields because of improved habitat. The increasing importance of wildlife habitat is reflected in continuing refinement of the Environmental Benefits Index (EBI) used by the U.S. Department of Agriculture (USDA) to quantify the potential benefits of enrolling lands in CRP (Osborn, 1997; Ribaudo and others, 2001). The refinements reflect input furnished by federal, state, and non-government organizations seeking greater wildlife habitat quality on CRP lands (Roseberry and David, 1994; Hughes and others, 1995; Millenbah and others, 1996; Patterson and Best, 1996; Rodgers, 1999; Allen and others, 2001). Refinement in the EBI has changed the types of grasses planted on newly enrolled land. In early CRP signups (1 through 11), 71% of new grassland acres were planted to introduced grasses and legumes [Conservation Practice (CP) 1] while 29% of the acres were planted to native grasses (CP2) (Osborn and others, 1992). By the 27th signup in July 2004, over 34.8 million acres (14 million ha) were enrolled in the CRP. More than 73% of these lands were planted to various mixtures of introduced (CP1) or native (CP2) grasses for a minimal contract period of 10 years (USDA, 2004). Continuation of grass plantings under the 2002 Farm Bill may result in CRP lands furnishing grass dominated cover for 20 or more consecutive years. The species of grass established in seeded grasslands can have a major influence on the potential quality of wildlife habitat where vegetation is maintained over a multi-year period. Different species of grass may have comparable abilities to alleviate soil erosion but furnish dissimilar qualities of wildlife habitat (fig. 1). For example, smooth brome, an introduced cool-season grass (grass species and scientific names are presented in table 1), is highly valued for its erosion control and forage attributes (Casler and Carlson, 1995). Switchgrass, a native warm-season grass, also is valued for its soil and water conservation qualities (Moser and Vogel, 1995) but provides greater benefits for some species of wildlife (Clubine, 1995). The quality of nesting and winter cover for ring-necked pheasants (Phasianus colchicus) furnished by smooth brome on northeastern Colorado CRP lands is inferior to that provided by the taller, more robust switchgrass (Allen, 1994). Characteristics of the agricultural landscape surrounding individual CRP fields also play a role in the wildlife habitat potential of CRP plantings (Weber and others, 2002; Nusser and others, 2004). of native and seeded grasslands change in response to the presence (and absence) of physical disturbances such as fire, grazing, tillage, and haying (Hobbs and Huenneke, 1992; Millenbah and others, 1996; Allen and others, 2001; Renfrew and Ribic, 2001; Swengel and Swengel, 2001). The perpetuation of diversity in species composition and vegetation structure following disturbance sustains desirable habitat for a variety of grassland-dependent wildlife (Hall and Willig, 1994; Barnes and others, 1995; King and Savidge, 1995; Granfors and others, 1996; Herkert and others, 1996; Kurzejeski, 1996; Patterson and Best, 1996; Klute and others, 1997). Undisturbed grasslands have lower grass and forb species diversity, greater amounts of dead plant material, decreased as well as seasonally delayed productivity, and diminished structural diversity of vegetation (Peet and others, 1975; Rice and Parenti, 1978; Butler and Briske, 1988; Campa and Winterstein, 1992). Recommendations for the timing of disturbance to increase grass and forb species diversity range from 3 to 8 years following establishment of seeded grasslands in the northern Great Plains and Midwest (Duebbert and others, 1981; Higgens, 1987; Millenbah and others, 1996). The management interval, however, is affected by climatic conditions, soils, grass species, and management history of the individual stand. We quantified changes in vegetation structure and species composition across the typical 10-year contract period in undisturbed southern and central Great Plains CRP fields (fig. 2) planted to introduced and native grasses. In addition, we compared changes in vegetation in fields grazed during the emergency release of 1996 by comparing conditions prior to grazing and two and four years post grazing relative to changes in similar fields that were not grazed. Documentation of long-term changes in vegetation structure and composition for fields planted to common grass seed mixtures across a wide range of environmental conditions provides information to improve long-term wildlife habitat potential, guide program administration, and define management practices that yield economic benefits to operators while still meeting wildlife and conservation objectives. Emergency grazing provisions of the CRP are controversial. Although grazing can alter vegetative characteristics and reduce habitat quality in the short-term (Temple and others, 1999), periodic disturbance may be necessary to maintain habitat quality, and more information is needed assessing long-term effects of emergency grazing on vegetative structure and species composition.

Conference Paper

Numerical-simulation and conjunctive-management models of the Hunt-Annaquatucket-Pettaquamscutt stream-aquifer system, Rhode Island

Numerical-simulation and optimization techniques were used to evaluate alternatives for the conjunctive management of ground- and surface-water resources of the Hunt­-Annaquatucket­-Pettaquamscutt stream-aquifer system in central Rhode Island. Ground-water withdrawals from the Hunt-Annaquatucket-Pettaquamscutt aquifer exceeded 8 million gallons per day during months of peak water use during 1993­98, and additional withdrawals have been proposed to meet growing demands from within and outside of the system boundary. The system is defined by the Hunt­-Annaquatucket­-Pettaquamscutt aquifer, which is composed of glacial stratified deposits, and the network of rivers, brooks, and ponds that overlie and are in hydraulic connection with the aquifer. Nearly all of the water withdrawn, however, is derived from depletions of flow in the rivers, brooks, and ponds that overlie the aquifer. Streamflow depletions are of concern to environmental agencies because of the adverse effects that reductions in streamflow can have on aquatic and riparian ecosystems. A conjunctive-management model of the stream-aquifer system was developed to simultaneously address the water-demand and streamflow-depletion issues. The objective of the model was to maximize total ground-water withdrawal from the aquifer during July, August, and September. These three months are generally the time of year when water-supply demands are largest and streamflows are simultaneously lowest. Total withdrawal from the aquifer was limited by a set of constraints specified in the model. These constraints were (1) maximum rates of streamflow depletion in the Hunt, Annaquatucket, and Pettaquamscutt Rivers; (2) minimum monthly water demands of each of three water-supply systems that withdraw water from the aquifer; and (3) minimum and maximum withdrawal rates at each supply well. The conjunctive-management model was formulated mathematically as a linear program. The model was solved by a response-matrix technique that incorporates the results of transient, numerical simulation of the stream-aquifer system into the constraint set of the linear program. The basis of the technique was the assumption that streamflow-depletion rates in each river were a linear function of ground-water-withdrawal rates at each well. This assumption was shown to be valid for the conditions evaluated in this study, primarily because of the very high transmissivity of the aquifer near many of the wells pumped for water supply. A transient, numerical model of the system was developed to simulate an average annual cycle of monthly withdrawal and hydrologic conditions representative of the 56-year period 1941­96. The transient model was used to generate characteristic streamflow-depletion responses in each river to simulated withdrawals at each well; these characteristic responses, or response coefficients, were then incorporated directly into the streamflow-depletion constraints of the linear program. Four sets of applications of the conjunctive-management model were made to determine whether total ground-water withdrawal from the aquifer during July, August, and September could be increased over the current total withdrawal for alternative definitions of the maximum rates of streamflow depletion allowed in the Hunt, Annaquatucket, and Pettaquamscutt Rivers. Current conditions were defined as the average monthly withdrawal rates at each supply well, water demands of each of the three water-supply systems, and estimated streamflow-depletion rates during the 6-year period 1993­98. Total withdrawal from all wells in the system from July through September during 1993­98 was 506.5 million gallons. Estimated streamflow-depletion rates for 1993­98 were calculated by use of the transient model, with the 1993­98 average monthly withdrawal rates specified at each supply well. Streamflow-depletion rates calculated for July, August, and September averaged 25 percent of the model-calculated pre-withdrawal streamflow rates for the Hunt River, 19 percent for the Annaquatucket River, and 7 percent for the Pettaquamscutt River. The first set of applications of the model were made with the current estimated rates of streamflow depletion in the Hunt, Annaquatucket, and Pettaquamscutt Rivers. Results of these applications indicated that total withdrawal from the aquifer during July, August, and September could be increased from about 8 to 18 percent (from 546.0 to 596.3 million gallons) over the current total withdrawal. The increased withdrawal would require modifications to the current annual withdrawal schedule of each supply well and, for the 18-percent increase, a modified network of supply wells that would include two new wells in the Annaquatucket River Basin. A second set of model applications then was made to determine if current estimated rates of streamflow depletion in the Hunt River could be reduced without increasing current estimated rates of streamflow depletion in the Annaquatucket or Pettaquamscutt Rivers. Decreases in the current rates of streamflow depletion in the Hunt River would result in increased streamflow in the river during these three months. Results showed that current rates of streamflow depletion in the Hunt River during July, August, and September could be decreased from 5 to 15 percent, depending on whether the existing or modified well network was used. Subsequent model applications indicated that substantial increases in total ground-water withdrawal from the aquifer are possible, but would require increased rates of streamflow depletion in the Annaquatucket and Pettaquamscutt Rivers. Maximum increases in the July through September withdrawal from the aquifer of about 39 to 50 percent (from 705.1 to 760.3 million gallons) over the current total withdrawal were calculated when streamflow-depletion rates in the Annaquatucket and Pettaquamscutt Rivers were allowed to increase from current estimated rates to a maximum of 25 percent of the model-calculated pre-withdrawal streamflow for each river during July, August, and September. Alternatively, it was shown that current estimated rates of streamflow depletion in the Hunt River during July, August, and September could be reduced by as much as 35 percent for the maximum allowed increases in streamflow depletion in the Annaquatucket and Pettaquamscutt Rivers; maximum increased withdrawal from the aquifer, however, would range from 8 to 18 percent over the current total withdrawal for the 35-percent reduction in streamflow-depletion rates in the Hunt River. Results of the different applications of the model demonstrate the usefulness of coupling numerical-simulation and optimization techniques for regional-scale evaluation of water-resource management alternatives. The results of the evaluation must be viewed, however, within the limitations of the quality of data available for the Hunt­-Annaquatucket­-Pettaquamscutt stream-aquifer system and representation of the system by a simulation model. An additional limitation of the analysis was the use of an average annual cycle of monthly withdrawal and hydrologic conditions. Ground-water withdrawal strategies may need to be modified to meet streamflow-depletion constraints during extreme hydrologic events, such as droughts. Contributing areas and sources of water to the supply wells also were delineated by use of a steady-state model of the stream-aquifer system. The model was developed to simulate long-term-average ground-water flow and ground-water/ surface-water interactions in the system during the 56-year period 1941­96. Sources of water to the wells consisted of precipitation and wastewater recharge to the aquifer, streamflow leakage from natural stream-channel losses, streamflow leakage caused by induced infiltration, and lateral ground-water inflow from till and bedrock upland areas.

Rhode Island

Forecasting selenium discharges to the San Francisco Bay-Delta Estuary: Ecological effects of a proposed San Luis Drain extension

Selenium discharges to the San Francisco Bay-Delta Estuary (Bay-Delta) could change significantly if federal and state agencies (1) approve an extension of the San Luis Drain to convey agricultural drainage from the western San Joaquin Valley to the North Bay (Suisun Bay, Carquinez Strait, and San Pablo Bay); (2) allow changes in flow patterns of the lower San Joaquin River and Bay-Delta while using an existing portion of the San Luis Drain to convey agricultural drainage to a tributary of the San Joaquin River; or (3) revise selenium criteria for the protection of aquatic life or issue criteria for the protection of wildlife. Understanding the biotransfer of selenium is essential to evaluating effects of selenium on Bay-Delta ecosystems. Confusion about selenium threats to fish and wildlife stem from (1) monitoring programs that do not address specific protocols necessary for an element that bioaccumulates; and (2) failure to consider the full complexity of the processes that result in selenium toxicity. Past studies show that predators are more at risk from selenium contamination than their prey, making it difficult to use traditional methods to predict risk from environmental concentrations alone. This report presents an approach to conceptualize and model the fate and effects of selenium under various load scenarios from the San Joaquin Valley. For each potential load, progressive forecasts show resulting (1) water-column concentration; (2) speciation; (3) transformation to particulate form; (4) particulate concentration; (5) bioaccumulation by invertebrates; (6) trophic transfer to predators; and (7) effects on those predators. Enough is known to establish a first-order understanding of relevant conditions, biological response, and ecological risks should selenium be discharged directly into the North Bay through a conveyance such as a proposed extension of the San Luis Drain. The approach presented here, the Bay-Delta selenium model, determines the mass, fate, and effects of selenium released to the Bay-Delta through use of (1) historical land-use, drainage, alluvial-fill, and runoff databases; (2) existing knowledge concerning biogeochemical reactions and physiological parameters of selenium (e.g., speciation, partitioning between dissolved and particulate forms, and bivalve assimilation efficiency); and (3) site-specific data mainly from 1986 to 1996 for clams and bottom-feeding fish and birds. Selenium load scenarios consider effluents from North Bay oil refineries and discharges of agricultural drainage from the San Joaquin Valley to enable calculation of (a) a composite freshwater endmember selenium concentration at the head of the estuary; and (b) a selenium concentration at a selected seawater location (Carquinez Strait) as a foundation for modeling. Analysis of selenium effects also takes into account the mode of conveyance for agricultural drainage (i.e., the San Luis Drain or San Joaquin River); and flows of the Sacramento River and San Joaquin River on a seasonal or monthly basis. Load scenarios for San Joaquin Valley mirror predictions made since 1955 of a worsening salt (and by inference, selenium) build-up exacerbated by an arid climate and massive irrigation. The reservoir of selenium in the San Joaquin Valley is sufficient to provide loading at an annual rate of approximately 42,500 pounds of selenium to a Bay-Delta disposal point for 63 to 304 years at the lower range of projections presented here, even if influx of selenium from the California Coast Ranges could be curtailed. Disposal of wastewaters on an annual basis outside of the San Joaquin Valley may slow the degradation of valley resources, but drainage alone cannot alleviate the salt and selenium build-up in the San Joaquin Valley, at least within a century. Load scenarios also show the different proportions of selenium loading to the Bay-Delta. Oil refinery loads from 1986 to 1992 ranged from 8.5 to 20 pounds of selenium per day; with treatment and cleanup, loads decreased to 3.0 pounds of selenium per day in 1999. In contrast, San Joaquin Valley agricultural drainage loads disposed of in a San Luis Drain extension could range from 45 to 117 pounds of selenium per day across a set of historical and future conditions. Components of this valley-wide load include five source subareas (i.e., Grassland, Westlands, Tulare, Kern, and Northern) defined by water and drainage management. Loads vary per subarea mainly because of proximity of the subarea to geologic sources of selenium and irrigation history. Loads from the Sacramento River, depending on flow conditions, range from 0.8 to 10 pounds of selenium per day. Loads from the San Joaquin River vary depending on restoration and flow conditions, which are considered. A consistent picture of ecological risk emerges under modeled selenium discharges from a proposed San Luis Drain extension. The threat to the estuary is greatest during low flow seasons and critically dry years. Where selenium undergoes reactions typical of low flow or longer residence time, highly problematic bioaccumulation in prey (food) is forecast. Surf scoter, greater and lesser scaup, and white sturgeon appear to be most at risk because these Bay-Delta predators feed on deposit and filter-feeding bivalves. Recent findings add Sacramento splittail and Dungeness crab to that list. During the low flow season of critically dry years, forecasted selenium concentrations in water, particulate matter, prey (diet), and predator tissue exceed guidelines with a high certainty of producing adverse effects under the most likely load scenario from a proposed San Luis Drain extension. High flows afford some protection under certain conditions in modeled San Joaquin River scenarios. However, meeting a combined goal of releasing a specific load during maximum flows and keeping selenium concentrations in the river below a certain objective to protect against bioaccumulation may not always be attainable. Management of the San Joaquin River on a constant concentration basis also could create problematic bioaccumulation during a wet year, especially during the low flow season, because high flows translate to high loads that are not always offset by seasonal river inflows. Prior to refinery cleanup, selenium contamination was sufficient to threaten reproduction in key species within the Bay-Delta ecosystems and human health advisories were posted based on selenium concentrations in tissues of diving ducks. During this time, selenium concentrations in the Bay-Delta were well below the most stringent recommended water quality criterion [1 microgram per liter (1 µg/L)]. Enhanced biogeochemical transformations to bioavailable particulate selenium and efficient bioaccumulation by bivalves characterized the system. If these biogeochemical conditions continue to prevail and agricultural selenium sources replace or exceed refinery sources, ecological forecasts suggest the risk of adverse effects will be difficult to eliminate under an out-of-valley resolution to the selenium problem. The Bay-Delta selenium model presented here is a systematic approach for conducting forecasts of the ecological effects from selenium on aquatic food webs. It is a new tool that links and models the major processes leading from loads through consumer organisms to predators. It also is a feasible approach for site-specific analysis and could provide a framework for developing new protective selenium foodweb guidelines and predator criteria. Model components that help ensure understanding ecosystems and the basis of environmental protection are (1) contaminant concentrations and speciation in sources, such as particulate material, that most influence bioavailability; (2) bioaccumulation models that calculate concentrations in diet, specifically in bivalves of the Bay-Delta that act as sensitive indicators of selenium contamination; (3) food-web type that determines what animals are threatened and when; and (4) multiple media concentrations (water, particulate material, and tissue of prey and predators) that, in-combination, determine risk or hazard.

California