Search USGS⌕ Search

SEARCH · Search USGS

Results for “Internal Report”

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 901 records · Page 50Linked to original sources
USGS2020

Simulation of discharge, water-surface elevations, and water temperatures for the St. Louis River estuary, Minnesota-Wisconsin, 2016–17

The St. Louis River estuary is a large freshwater estuary, next to Duluth, Minnesota, that encompasses the headwaters of Lake Superior. The St. Louis River estuary is one of the most complex and compromised near-shore systems in the upper Great Lakes with a long history of environmental contamination caused by logging, mining, paper mills, and other heavy industrial activities. Presently (2020), a widely available, science-based assessment tool capable of evaluating ecosystem-level responses to remediation and restoration projects has not existed for the estuary. To address this need, the U.S. Geological Survey (USGS) built a predictive, mechanistic, three-dimensional hydrodynamic model for the estuary using the Environmental Fluid Dynamics Code framework. In the current version, the model can simulate continuous discharge, water-surface elevations, water temperature, and flow velocity, although the modular framework allows for future additions of water-quality modeling. The model was calibrated using data collected from April 2016 through November 2016 and validated with data collected from April 2017 through November 2017. The four types of data used to evaluate model performance were water-surface elevations, discharge, water temperature, and flow velocities. Streamflow and temperature boundary condition data included a mixture of USGS streamgage data, Minnesota Department of Natural Resources gage data, and estimates derived from the gage data. The model was able to simulate the water-surface elevations with generally good agreement between the simulated and measured values for both years at the daily time step. Specifically, the model was able to demonstrate excellent agreement with the measured data with Nash-Sutcliffe efficiency coefficients greater than 0.8 for all three locations; however, the model was unable to produce hourly water-surface elevations with such accuracy for 2016–17. Discharge was more dynamic than the water-surface elevations, both for the measured and simulated data. Generally, most of the discharge ranged from −650 to 1,200 cubic meters per second, but the constantly changing flux exiting the estuary into Lake Superior (positive flows) and entering the estuary from Lake Superior (negative flows) occurred throughout the year. Even upstream at the St. Louis River at Oliver, Wisconsin, gage (USGS station 0402403250), the effect of flows into the estuary from Lake Superior did occur, demonstrating the strong effect of the Lake Superior seiche on flows for the estuary. From a performance standpoint, the model was able to simulate discharge with generally good agreement in both years, although the 2017 validation was better than the 2016 calibration period. For the daily Nash-Sutcliffe efficiency coefficients, the simulated values were 0.98, 0.62, 0.49, and 0.71 for the Oliver gage; the Superior Bay entry channel at Superior, Wisc., (USGS station 464226092005600); the Superior Bay Duluth Ship Canal at Duluth, Minn., (USGS station 464646092052900); and total entries (combination of the Superior entry and Duluth entry), respectively. For the hourly evaluation criteria, the model performed poorly, with Nash-Sutcliffe efficiency coefficients less than 0 for the two entries into Lake Superior; therefore, as a predictor of discharge at the hourly scale, the model performed worse than using the measured data average. Similar to discharge, the model was a good predictor of flow velocity at the daily time scale but had difficulty matching the measured data at the hourly scale. For discharge and flow velocity, matching at subdaily time steps for a system as complicated as the St. Louis River estuary is considered difficult because the match is highly sensitive to coordinating the exact measurement location to the simulated value. The final calibration target was water temperature, calibrated for the Oliver gage and the Duluth entry. For calibration purposes, the Duluth entry was the more important water temperature target because the Oliver gage was more of an internal check on the model. The Nash-Sutcliffe efficiency coefficients for the Duluth entry were high; hourly Nash-Sutcliffe efficiency coefficients at the Duluth entry were either at or greater than 0.7 for both years, and daily values were 0.84 and 0.82 for 2016 and 2017, respectively.

Minnesota, Wisconsin↗
USGS2011

Summary of November 2010 meeting to evaluate turbidite data for constraining the recurrence parameters of great Cascadia earthquakes for the update of national seismic hazard maps

This report summarizes a meeting of geologists, marine sedimentologists, geophysicists, and seismologists that was held on November 18–19, 2010 at Oregon State University in Corvallis, Oregon. The overall goal of the meeting was to evaluate observations of turbidite deposits to provide constraints on the recurrence time and rupture extent of great Cascadia subduction zone (CSZ) earthquakes for the next update of the U.S. national seismic hazard maps (NSHM). The meeting was convened at Oregon State University because this is the major center for collecting and evaluating turbidite evidence of great Cascadia earthquakes by Chris Goldfinger and his colleagues. We especially wanted the participants to see some of the numerous deep sea cores this group has collected that contain the turbidite deposits. Great earthquakes on the CSZ pose a major tsunami, ground-shaking, and ground-failure hazard to the Pacific Northwest. Figure 1 shows a map of the Pacific Northwest with a model for the rupture zone of a moment magnitude M w 9.0 earthquake on the CSZ and the ground shaking intensity (in ShakeMap format) expected from such an earthquake, based on empirical ground-motion prediction equations. The damaging effects of such an earthquake would occur over a wide swath of the Pacific Northwest and an accompanying tsunami would likely cause devastation along the Pacifc Northwest coast and possibly cause damage and loss of life in other areas of the Pacific. A magnitude 8 earthquake on the CSZ would cause damaging ground shaking and ground failure over a substantial area and could also generate a destructive tsunami. The recent tragic occurrence of the 2011 M w 9.0 Tohoku-Oki, Japan, earthquake highlights the importance of having accurate estimates of the recurrence times and magnitudes of great earthquakes on subduction zones. For the U.S. national seismic hazard maps, estimating the hazard from the Cascadia subduction zone has been based on coastal paleoseismic evidence of great earthquakes over the past 5,000 years. The instrumental catalog of earthquakes is of little use for constraining the hazard of the CSZ, because there are virtually no recorded earthquakes on most of the plate interface of the CSZ. There are no historical accounts in the past 150 years of large earthquakes on most of the CSZ. Until about 20 years ago, some interpreted this lack of recent and historical earthquakes as an indicator that the subduction zone was slipping aseismically and could not produce a great earthquake. The work of Brian Atwater and others, in the late 1980s and the 1990s (Atwater, 1987, 1992; Atwater and others, 1995; Nelson and others, 1996; Clague, 1997; Atwater and Hemphill-Haley, 1997; Atwater and others, 2004) demonstrated that submerged forests, buried soils, tsunami deposits, and liquefaction along and near the coast were compelling evidence of repeated great Cascadia earthquakes over at least the past 5,000 years. Atwater and Hemphill-Haley (1997) concluded from paleoseismic evidence at Willapa Bay, Washington, that great earthquakes ruptured the CSZ with an average recurrence time of about 500 years. The date of the last great CSZ earthquake, January 26, 1700, was established from historical records of the so-called orphan tsunami in Japan that is inferred to have been produced by this earthquake (Satake and others, 1996, 2003; Atwater and others, 2005) and is consistent with tree-ring data from drowned forests in Washington and Oregon. From modeling the observations of the tsunami, Satake and others (2003) estimated a moment magnitude of about 9.0 for this earthquake. Many other paleoseismic sites have been investigated along the Pacific Northwest coast from Vancouver Island to northern California and show evidence of great CSZ earthquakes. Nelson and others (2006) summarized the dates found from these studies and proposed correlations between sites indicating the extent of rupture for individual events. Dating of inferred tsunami deposits in Bradley Lake, Oregon by Kelsey and others (2005), as well as tsunami and subsidence evidence from Six Rivers, Oregon (Kelsey and others, 2002) and Coquille River (Witter and others, 2003), indicates that there were probably M w 8 ruptures in the southern portion of the CSZ in addition to the M w 9 events that rupture the whole length of the CSZ (Nelson and others, 2006). A parallel development over the past 20 years or more is the use of deep-sea turbidite deposits for identifying and dating great Cascadia earthquakes over the past 10,000 years (Adams, 1990; Goldfinger and others, 2003, 2008, in press; Goldfinger, 2011). Turbidites are sediment deposits in the deep ocean from turbidity currents, which are energetic flows of sediment and water along the continental shelf and slope. Adams (1990), using the counts of turbidites in deep-sea cores off the coast of Oregon and Washington collected and analyzed by Griggs (1969) and Griggs and others (1969), proposed that these turbidites were caused by the shaking of great Cascadia earthquakes. Part of his reasoning was that the number (13) of turbidite deposits that occurred since deposition of the Mazama Ash 7,000 years ago gave a recurrence time of about 500 years, consistent with that derived from the coastal submergence data. Adams (1990) also proposed the “confluence test” which evaluates the number of turbidites for submarine channels that form a confluence. He reported that the number of turbidites in the single downstream channel equaled the number in each of the tributary channels. He reasoned that this indicated that the turbidites in each tributary were simultaneously triggered and were, therefore, caused by a common forcing agent. He concluded that shaking from extended ruptures of great Cascadia earthquakes was the most likely cause of these turbidites. Based on the paleoseismic evidence of past great earthquakes, the hazard from the Cascadia subduction zone was included in the 1996 U.S. NSHM (Frankel and others, 1996), which were the basis for seismic provisions in the 2000 International Building Code. These hazard maps used the paleoseismic studies to constrain the recurrence rate of great CSZ earthquakes. Goldfinger and his colleagues have since collected many more deep ocean cores and done extensive analysis on the turbidite deposits that they identified in the cores (Goldfinger and others, 2003, 2008, in press; Goldfinger, 2011). Using their dating of the sediments and correlation of features in the logs of density and magnetic susceptibility between cores, they developed a detailed chronology of great earthquakes along the CSZ for the past 10,000 years (Goldfinger and others, in press). These correlations consist of attempting to match the peaks and valleys in logs of density and magnetic susceptibility between cores separated, in some cases, by hundreds of kilometers. Based on this work, Goldfinger and others (2003, 2008, in press) proposed that the turbidite evidence indicated the occurrence of great earthquakes (M w 8) that only ruptured the southern portion of the CSZ, as well as earthquakes with about M w 9 that ruptured the entire length of the CSZ. For the southernmost portion of the CSZ, Goldfinger and others (in press) proposed a recurrence time of M w 8 or larger earthquakes of about 230 years. This proposed recurrence time was shorter than the 500 year time that was incorporated in one scenario in the NSHM’s. It is important to note that the hazard maps of 1996 and later also included a scenario or set of scenarios with a shorter recurrence time for M w 8 earthquakes, using rupture zones that are distributed along the length of the CSZ (Frankel and others, 1996; Petersen and others, 2008). Originally, this scenario was meant to correspond to the idea that some of the 500-year averaged ruptures seen in the paleoseismic evidence could have been a series of M w 8 earthquakes that occurred over a short period of time (a few decades), rather than M w 9 earthquakes. Figure 2 shows the logic tree for the CSZ used in the 2008 NSHM’s (Petersen and others, 2008). This logic tree includes whole CSZ rupture earthquakes (M w 8.8–9.2) and partial CSZ rupture earthquakes (M w 8.0–8.7). In this latest version of the NSHM’s, the effective recurrence time of earthquakes on the CSZ with moment magnitudes greater than or equal to 8.0 over the various models is about 270 years (Petersen and others, 2008). This recurrence time applies to the entire CSZ, so that the hazard from great earthquakes was approximately equal along the whole zone, although the hazard estimates taper on the northern and southern ends of the CSZ, because of the way rupture zones of M w 8 earthquakes were distributed along the strike of the CSZ. The NSHM will be updated in 2013, as part of the standard update cycle that corresponds to the update cycle of the national model building codes that are based on the seismic hazard maps. A meeting was necessary to assemble a wide group of experts to hear Dr. Goldfinger explain his methodology for dating and correlating the turbidites and for developing the earthquake chronology. The overall goal of the workshop was to evaluate observations of turbidite deposits to provide constraints on the recurrence times and rupture extents of great Cascadia subduction zone earthquakes for the next update of the NSHM. Before the meeting, participants were supplied with the U.S. Geological Survey (USGS) Professional Paper of Goldfinger and others (in press), as well as material from Brian Atwater and Alan Nelson. The agenda of the meeting was developed by Art Frankel, with assistance from Chris Goldfinger, Brian Atwater, Alan Nelson, Mark Petersen, and Craig Weaver. The meeting was hosted by Chris Goldfinger of Oregon State University. We stress that it is difficult to evaluate in a two-day meeting the large amount of work that Goldfinger and his colleagues have done over the past 15 years or more. This meeting is the first step in a process that develops the inputs to the update of the national maps. The conclusions of this workshop will be discussed and possibly modified at the regional Pacific Northwest workshop for the hazard maps to be held in early 2012. Vetting new research results using informed expert opinion is an integral part of updating the national maps and does not reflect on the veracity of these results.

Cascadia↗
USGS2008

Mycobacteriosis in striped bass

Mycobacteriosis is a bacterial disease in which striped bass (rockfish) may be disfigured as a result of skin ulcers and internal lesions. The bass may also be skinny or in extremely poor condition due to the chronic nature of this wasting disease. Stripers are a highly prized target species for both recreational anglers and commercial fishermen. As such, the economic impact of diseased and devalued fish could be significant. In addition, some of the mycobacteria that commonly infect fishes can cause infections in people and therefore are a human health concern. The total extent to which the disease is occurring along the Eastern seaboard is unknown but the disease has been reported from stripers taken from North Carolina to New York. During 1998-99, skin ulcers attributed to mycobacterial infection were observed in up to 28% of the striped bass from some Virginia tributaries of the Chesapeake Bay. Data obtained during 2002—2003 from fish harvested in Virginia and Maryland waters indicated that, at least in some areas, over 80% of striped bass may be infected with the mycobacteria that are associated with the disease. Given the persistence over the last 8 years of this mycobacteriosis outbreak, this does not appear to be a short-term problem.

Fact Sheet↗
USGS2010

Public review draft: A method for assessing carbon stocks, carbon sequestration, and greenhouse-gas fluxes in ecosystems of the United States under present conditions and future scenarios

The Energy Independence and Security Act of 2007 (EISA), Section 712, authorizes the U.S. Department of the Interior to develop a methodology and conduct an assessment of the Nation's ecosystems focusing on carbon stocks, carbon sequestration, and emissions of three greenhouse gases (GHGs): carbon dioxide, methane, and nitrous oxide. The major requirements include (1) an assessment of all ecosystems (terrestrial systems, such as forests, croplands, wetlands, shrub and grasslands; and aquatic ecosystems, such as rivers, lakes, and estuaries), (2) an estimation of annual potential capacities of ecosystems to increase carbon sequestration and reduce net GHG emissions in the context of mitigation strategies (including management and restoration activities), and (3) an evaluation of the effects of controlling processes, such as climate change, land use and land cover, and wildlfires. The purpose of this draft methodology for public review is to propose a technical plan to conduct the assessment. Within the methodology, the concepts of ecosystems, carbon pools, and GHG fluxes used for the assessment follow conventional definitions in use by major national and international assessment or inventory efforts. In order to estimate current ecosystem carbon stocks and GHG fluxes and to understand the potential capacity and effects of mitigation strategies, the method will use two time periods for the assessment: 2001 through 2010, which establishes a current ecosystem GHG baseline and will be used to validate the models; and 2011 through 2050, which will be used to assess future potential conditions based on a set of projected scenarios. The scenario framework is constructed using storylines of the Intergovernmental Panel on Climate Change (IPCC) Special Report Emission Scenarios (SRES), along with initial reference land-use and land-cover (LULC) and land-management scenarios. An additional three LULC and land-management mitigation scenarios will be constructed for each storyline to enhance carbon sequestration and reduce GHG fluxes in ecosystems. Input from regional experts and stakeholders will be solicited to construct realistic and meaningful scenarios. The methods for mapping the current LULC and ecosystem disturbances will require the extensive use of both remote-sensing data and in-situ (for example, forest inventory data) to capture and characterize landscape-change events. For future potential LULC and ecosystem disturbances, key drivers such as socioeconomic, policy, and climate assumptions will be used in addition to biophysical data. The product of these analyses will be a series of maps for each future year for each scenario. These annual maps will form the basis for estimating carbon storage and GHG emissions. For terrestrial ecosystems, carbon storage, carbon-sequestration capacities, and GHG emissions under the current and projected future conditions will be assessed using the LULC and ecosystem-disturbance estimates in map format with a spatially explicit biogeochemical ensemble modeling system that incorporates properties of management activities (such as tillage or harvesting) and properties of individual ecosystems (such as elevation, vegetation characteristics, and soil attributes). For aquatic ecosystems, carbon burial in sediments and GHG fluxes are functions of the current and projected future stream flow and sediment transports, and therefore will be assessed using empirical modeling methods. Validation and uncertainty analysis methods described in the methodology will follow established guidelines to assess the quality of the assessment results. The U.S. Environmental Protection Agency's Level II ecoregions map (which delineates 24 ecoregions for the Nation) will be the practical instrument for developing and delivering assessment results. Consequently, the ecoregion will be the reporting unit of the assessment because the mitigation scenarios, assessment results, validation, and uncertainty analysis will be

Open-File Report↗
USGS2023

Survey evaluation of Florida’s freshwater fisheries long-term monitoring program

Natural resource monitoring programs benefit from routine evaluation. Here, Florida’s statewide Freshwater Fisheries Long-Term Monitoring (LTM) program is used to show how stakeholder surveys can be integral to this process. In 2022, an online questionnaire was sent to internal stakeholders, i.e., state agency personnel who collect, enter, or use freshwater fisheries data for fisheries and habitat management purposes. The survey’s primary objective was to evaluate the program at its 15-year mark; secondary objectives were to compare results with a similar survey conducted at the 4-year mark, compare results among respondents based on experience and functional role, and develop recommendations for strategic initiatives to further improve the program. The survey consisted of 43 questions across six sections of program evaluation: demographics; field sampling; data entry, summary, and reporting; management decision support; programmatic views; and additional input. Respondents generally had positive views of the LTM program, but the survey revealed differences among respondents with different functional roles (e.g., fisheries researchers and managers viewed the decisional value, priority, and sample sizes of LTM data more favorably than habitat managers) while highlighting high-priority future initiatives (e.g., database development). Our results demonstrate the utility of stakeholder surveys as an important step in evaluating monitoring programs.

Florida↗
USGS2002

Challenges in observational seismology

Earthquake seismology became a quantitative scientific discipline after instruments were developed to record seismic waves in the late 19th century ( Dewey and Byerly, 1969 ; Chapter 1 by Agnew). Earthquake seismology is essentially based on field observations. The great progress made in the past several decades was primarily due to increasingly plentiful and high-quality data that are readily distributed. Our ability to collect, process, and analyze earthquake data has been accelerated by advances in electronics, communications, computers, and software (see Chapter 85 edited by Snoke and Garcia-Fernandez). Instrumental observation of earthquakes has been carried out for a little over 100 years by seismic stations and networks of various sizes, from local to global scales (see Chapter 87 edited by Lahr and van Eck). The observed data have been used, for example, (1) to compute the source parameters of earthquakes, (2) to determine the physical properties of the Earth's interior, (3) to test the theory of plate tectonics , (4) to map active faults, (5) to infer the nature of damaging ground shaking, and (6) to carry out seismic hazard analysis. Construction of a satisfactory theory of the earthquake process has not yet been achieved within the context of physical laws. Good progress, however, has been made in building a physical foundation of the earthquake source process, partly as a result of research directed toward earthquake prediction. This chapter is intended for a general audience. Technical details are not given, but relevant references and chapters in this Handbook are referred to. The first part of this chapter presents a brief overview of the observational aspects of earthquake seismology, concentrating on instrumental observations of seismic waves generated by earthquakes (i.e., seismic monitoring), and readers are referred to Chapter 49 by Musson and Cecic for noninstrumental observations. A few key developments and practices are summarized by taking a general view, since many national and regional developments have been chronicled in national and institutional reports (see Chapter 79 edited by Kisslinger). In the latter part of this chapter, the nature of seismic monitoring and some challenges in observational seismology are discussed from a personal perspective. Comments of a technical or philosophical nature are given in the Notes at the end of the chapter, and they are referenced by superscript numbers in the text.

International Geophysics↗
USGS1989

Mineral occurrence and drill-hole location map of the International Falls 1 degree by 2 degrees quadrangle, Minnesota and Ontario

A project to assess the mineral potential of the U.S. part of the International Fall 1° x 2° quadrangle was undertaken by the U.S. Geological Survey (USGS) in the early 1984, as part of the Conterminous United States Mineral Assessment Program (CUSMAP). Because glacial drift covers most of the basement rocks south of the U.S.-Canada border and mineral production in this part of Minnesota has been virtually nonexistent, data from the adjacent part of Ontario has been used to understand the geologic framework and to evaluate mineral-deposit models relevant to the assessment of the U.S. part of the International Falls quadrangle. This map and accompanying tables show the locations of mineral occurrences, prospects and mines in both Canada and the United States, and summarize their lithologic association and the characteristic metals found in them. Locations of drill holes within the United States and completed before 1987 are also included on the map.

Minnesota, Ontario↗
USGS2010

After a century: Revised Paleogene coal stratigraphy, correlation, and deposition, Powder River Basin, Wyoming and Montana

The stratigraphy, correlation, mapping, and depositional history of coal-bearing strata in the Paleogene Fort Union and Wasatch Formations in the Powder River Basin were mainly based on measurement and description of outcrops during the early 20th century. Subsequently, the quality and quantity of data improved with (1) exploration and development of oil, gas, and coal during the middle 20th century and (2) the onset of coalbed methane (CBM) development during the late 20th and early 21st centuries that resulted in the drilling of more than 26,000 closely spaced wells with accompanying geophysical logs. The closeness of the data control points, which average 0.5 mi (805 m) apart, made for better accuracy in the subsurface delineation and correlation of coal beds that greatly facilitated the construction of regional stratigraphic cross sections and the assessment of resources. The drillhole data show that coal beds previously mapped as merged coal zones, such as the Wyodak coal zone in the Wyoming part of the Powder River Basin, gradually thinned into several discontinuous beds and sequentially split into as many as 7 hierarchical orders westward and northward. The thinning and splitting of coal beds in these directions were accompanied by as much as a ten-fold increase in the thicknesses of sandstone-dominated intervals within the Wyodak coal zone. This probably resulted from thrust loading by the eastern front of the Bighorn uplift accompanied by vertical displacement along lineaments that caused subsidence of the western axial part of the Powder River Basin during Laramide deformation in Late Cretaceous and early Tertiary time. Accommodation space was thereby created for synsedimentary alluvial infilling that controlled thickening, thinning, splitting, pinching out, and areal distribution of coal beds. Equally important was differential subsidence between this main accommodation space and adjoining areas, which influenced the overlapping, for example, of the Dietz coal zone in Montana, over the Wyodak coal zone in Wyoming. Correlation in a circular track of the Wyodak coal zone in the southern part of the basin also demonstrates overlapping with lower coal zones. Recognition of this stratigraphic relationship has led to revision of the correlations and nomenclature of coal beds because of inconsistency within these zones as well as those below and above them, which have long been subjects of controversy. Also, it significantly changes the traditional coal bed-to-bed correlations, and estimates of coal and coalbed methane resources of these coal zones due to thinning and pinching out of beds. More notably, thickness isopach, orientation, and distribution of the merged Wyodak coal bodies in the south-southeast part of the basin suggest that differential movement of lineament zones active during the Cretaceous was not a major influence on coal accumulation during the Paleocene. Improved knowledge of alluvial depositional environments as influenced by external and internal paleotectonic conditions within the Powder River Basin permits more accurate correlation, mapping, and resource estimation of the Fort Union and Wasatch coal beds. The result is a better understanding of the sedimentology of the basin infill deposits in relation to peat bog accumulation.

Montana, Wyoming↗
USGS2007

The health hazards of volcanic ash--A guide for the public

This document has been prepared by the International Volcanic Health Hazard Network (IVHHN), Cities and Volcanoes Commission, GNS Science and the United States Geological Survey (USGS) to promote the safety of those who experience volcanic ashfall. This guide explains the potential health effects of volcanic ash and gives details on how to protect yourself and your family in the event of a volcanic ashfall.

Report↗
USGS1979

Application of a new Raman microprobe spectrometer to nondestructive analysis of sulfate and other ions in individual phases in fluid inclusions in minerals

Rosasco et al . (1975), reported the first successful application of laser-excited Raman spectroscopy for the identification and nondestructive partial analysis of individual solid, liquid, and gaseous phases in selected fluid inclusions. We report here the results of the application of a new instrument, based on back-scattering, that eliminates many of the previous stringent sample limitations and hence greatly expands the range of applicability of Raman spectroscopy to fluid inclusions. Fluid inclusions in many porphyry copper deposits contain 5–10 μm ‘daughter’ crystals thought to be anhydrite but too small for identification by the previous Raman technique. Using the new instrument, we have verified that such daughter crystals in quartz from Bingham, Utah, are anhydrite. They may form by leakage of hydrogen causing internal autooxidation of sulfide ion. Daughter crystals were also examined in apatite (Durango, Mexico) and emerald (Muzo, Colombia). Valid analyses of sulfur species in solution in small fluid inclusions from ore deposits would be valuable, but are generally impossible by conventional methods. We present a calibration procedure for analyses for SO 4 2− in such inclusions from Bingham, Utah (12,000 ± 4000 ppm ) and Creede, Colo. (probably < 500 ppm). A fetid Brazilian quartz, originally thought to contain liquid H 2 S, is shown to contain only HS − in major amounts.

Geochimica et Cosmochimica Acta↗
USGS2020

Rethinking a groundwater flow system using a multiple-tracer geochemical approach: A case study in Moab-Spanish Valley, Utah

The Glen Canyon Group Aquifer (GCGA) is the sole source of public water supply for the city of Moab, Utah, a domestic and international tourist destination. Population and tourism growth are likely to target the GCGA for future water resources, but our analysis indicates that additional withdrawals would likely be sourced from groundwater storage and not be sustained by recharge. A quantitative estimate of groundwater discharge from the GCGA is problematic because the downgradient aquifer boundary is the Colorado River, and groundwater discharge to the river is very small compared to the river flow. A water budget based on a conceptual model of GCGA discharging into an adjacent alluvial Valley-Fill Aquifer (VFA) was reported by Sumsion (1971) and numerous subsequent studies have repeated and utilized this water budget. The GCGA contains stable isotopes, tritium, 3 He/ 4 He ratios, dissolved solids, and sulfate concentrations that contrast with the VFA, indicating it is instead recharged by local streams rather than from the GCGA. Water-budget calculations, based on: (1) measured spring discharge and streamflow gains, (2) horizontal gradients in VFA groundwater age, and (3) GCGA outcrop area vadose-zone pore waters are all less than previously thought. Using a lumped parameter model and 14 C groundwater ages, we estimate recharge to the deeper GCGA (DGCGA) to be 4.2 ± 2.3 × 10 6 m 3 /yr, which is approximately equal to the measured discharge from wells and springs.

Utah↗
USGS2010

A method for assessing carbon stocks, carbon sequestration, and greenhouse-gas fluxes in ecosystems of the United States under present conditions and future scenarios

he Energy Independence and Security Act of 2007 (EISA), Section 712, mandates the U.S. Department of the Interior to develop a methodology and conduct an assessment of the Nation’s ecosystems, focusing on carbon stocks, carbon sequestration, and emissions of three greenhouse gases (GHGs): carbon dioxide, methane, and nitrous oxide. The major requirements include (1) an assessment of all ecosystems (terrestrial systems, such as forests, croplands, wetlands, grasslands/shrublands; and aquatic ecosystems, such as rivers, lakes, and estuaries); (2) an estimate of the annual potential capacities of ecosystems to increase carbon sequestration and reduce net GHG emissions in the context of mitigation strategies (including management and restoration activities); and (3) an evaluation of the effects of controlling processes, such as climate change, land-use and land-cover change, and disturbances such as wildfires. The concepts of ecosystems, carbon pools, and GHG fluxes follow conventional definitions in use by major national and international assessment or inventory efforts. In order to estimate current ecosystem carbon stocks and GHG fluxes and to understand the potential capacity and effects of mitigation strategies, the method will use two time periods for the assessment: 2001 through 2010, which establishes a current ecosystem carbon and GHG baseline and will be used to validate the models; and 2011 through 2050, which will be used to assess potential capacities based on a set of scenarios. The scenario framework will be constructed using storylines of the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emission Scenarios (SRES), along with both reference and enhanced land-use and land-cover (LULC) and land-management parameters. Additional LULC and land-management mitigation scenarios will be constructed for each storyline to increase carbon sequestration and reduce GHG fluxes in ecosystems. Input from regional experts and stakeholders will be solicited to construct these scenarios. The methods for mapping the current LULC and ecosystem disturbances will require the extensive use of both remote-sensing data and field-survey data (for example, forest inventories) to capture and characterize landscape-changing events. For potential LULC changes and ecosystem disturbances, key drivers such as socioeconomic and climate changes will be used in addition to the biophysical data. The result of these analyses will be a series of maps for each future year for each scenario. These annual maps will form the basis for estimating carbon storage and GHG emissions. For terrestrial ecosystems, carbon storage, carbon-sequestration capacities, and GHG emissions under the present conditions and future scenarios will be assessed using the LULC-change and ecosystem-disturbance estimates in map format with a spatially explicit biogeochemical ensemble modeling system that incorporates properties of management activities (such as tillage or harvesting) and properties of individual ecosystems (such as energy exchange, vegetation characteristics, hydrological cycling, and soil attributes). For aquatic ecosystems, carbon burial in sediments and fluxes of GHG are functions of the present and future potential stream flow and sediment transport and will be assessed using empirical hydrological modeling methods. Validation and uncertainty analysis methods described in the methodology will follow established guidelines to assess the quality of the assessment results. The U.S. Environmental Protection Agency’s Level II ecoregions map will be the practical instrument for developing and delivering assessment results. Consequently, the ecoregion (there are 22 modified ecoregions) will be the reporting unit of the assessment because the scenarios, assessment results, validation, and uncertainty analysis will be produced at that scale. The implementation of these methods will require collaborations among various Federal agencies, State agencies, nongovernmental organizations, and the science community. Using the method described in this document, the assessment can be completed in approximately 3 to 4 years. The primary deliverables will be assessment reports containing tables, charts, and maps that will present the estimated GHG parameters annually for 2001 through 2050 by ecosystem, pool, and scenario. The results will permit the evaluation of a range of policies, mitigation options, and research topics, such as the demographic, LULC-change, or climate-change effects on carbon stocks, carbon sequestration, and GHG fluxes in ecosystems.

Scientific Investigations Report↗
USGS1960

Gravity survey in part of the Snake River Plain, Idaho — A preliminary report

During the early summer of 1959, a total of 1,187 gravity stations were occupied on the western part of the Snake River plain in Idaho. An area of 2,000 square miles extending from Glenns Ferry, Idaho, to Caldwell, Idaho, was covered with a station density of one station per two square miles. An additional 1,200 square miles of surrounding area, mainly from Caldwell, Idaho, to the Oregon-Idaho state line, was covered with a density of one station per seven square miles. The mean reproducibility of the observed gravities of these stations was 0.05 milligal, with a maximum discrepancy of 0.2 milligal. Gravity data were reduced to simple Bouguer values using a combined free-air and Bouguer correction of 0.06 milligal per foot. The only anomalies found with closure in excess of 10 milligals are two elongated highs, orientated northwest-southeast, with the northwestern high offset to the northeast by 10 miles. The smaller of these highs extends from Meridian, Idaho, to Nyssa, Oregon, and the larger extends from Swan Falls, Idaho, to Glenns Ferry, Idaho. The maximum value recorded is a simple Bouguer value of -66.5 milligals with respect to the International Ellipsoid. Gradients on the sides of these highs are largest on the northeast sides, reaching six milligals per mile in places. Graticule interpretations of a profile across the southeastern high using a density contrast of 0.3 gm per cubic centimeter indicate an accumulation of lava reaching a thickness of at least 28,000 feet. The Snake River investigation was made for the purpose of searching out, defining, and interpreting gravity anomalies present on the western part of the Snake River lava plain in Idaho. In particular, it was desired to further define gradients associated with the gravity high shown by the regional work of Bonini and Lavin (1957). It was not planned to cover any specific area, but rather to let the observed anomalies determine the course of the field work. The study was undertaken as part of a project on Volcanism and Crustal Deformation, supervised by L. C. Pakiser of the U.S. Geological Survey. Professor Rodgers of the Geophysics Department of the Colorado School of Mines acted as an advisor.

Idaho↗
USGS1984

A workshop model simulating fate and effect of drilling muds and cuttings on benthic communities

Oil and gas exploration and production at marine sites has generated concern over potential environmental impacts resulting from the discharge of spent drilling muds and cuttings. This concern has led to a broad array of publicly and privately sponsored research. This report described a cooperative modeling effort designed to focus information resulting from this research through construction of explicit equations that simulate the potential impacts of discharge drilling fluids (muds) and cuttings on marine communities. The model is the result of collaboration among more than 30 scientists. The principal cooperating organizations were the E.S. Environmental Protection Agency, the U.S. Minerals Management Service, the Offshore Operators Committee, and the Alaska Oil and Gas Association. The overall simulation model can be conceptualized as three connected submodels: Discharge and Plume Fate, Sediment Redistribution, and Benthic Community Effects. On each day of simulation, these submodels are executed in sequence, with flows of information between submodels. The Benthic Community Effects submodel can be further divided into sections that calculate mortality due to burial, mortality due to toxicity, mortality due to resuspension disturbance, and growth of the community. The model represents a series of seven discrete 1-m 2 plots at specified distances along a transect in one direction away from a discharge point. It consists of coupled difference equations for which parameter values can easily be set to evaluate different conditions or to examine the sensitivity of output to various assumptions. Sets of parameter values were developed to represent four general cases or scenarios: (1) a shallow (5 m), cold environment with ice cover during a substantial fraction of the year, such as might be encountered in the Beaufort Sea, Alaska; (2) a shallow (20 m), temperate environment, such as might be encountered in the Gulf of Mexico; (3) a deeper (80 m), temperate environment, such as might be encountered in the Gulf of Mexico; and (4) a very deep (1,000 m) environment, such as might be encountered on the Atlantic slope. The focus of the modeling effort was on the connection of a reasonable representation of physical fate to the biological responses of populations, rather than on highly detailed representations of individual processes. For example, the calculations of physical fate are not as detailed as those in the recently published model of Brandsma et al. (1983). The value of the model described herein is in the broad scope of processes that are explicitly represented and linked together. The model cannot be considered to produce reliable predictions of the quantitative impacts of discharged drilling fluids and cuttings on biological populations at a particular site. Limitations of the model in predicting integrated fate and effects can be traced to three general areas: level of refinement of the algorithms used in the model; lack of understanding of the processes determining fate and effects; and parameter and data values. Despite the limitations, several qualitative conclusions concerning both potential impacts and the importance of various remaining data gaps can be drawn from the modeling effort. These include: (1) Simple, unequivocal conclusions about fate and effects across geographical regions and drilling operations are difficult, if not misleading, due to the large amount of variability in characteristics of discharged materials (e.g., oil content and toxicity), discharge conditions (e.g., duration of drilling operations), physical environments (e.g., water depth, current direction, and sediment disturbance regimes), and biological communities (e.g., intrinsic growth rates). Different combinations of these characteristics can result in substantial differences in simulated environmental fate and biological effects. For examples, simulated recovery in some high-energy environments occurs within months after the cessation of discharge operations, even at heavily impacted sites, whereas simulated recover in some low-energy environments takes years at heavily impacted sites. <2) Considerable difficulties remain in the reliable extrapolation of results from laboratory toxicity experiments to predictions of population effects in the field. (3) The volume of material discharged and duration of operations in the production drilling operations simulated by the model are sufficient to produce substantial simulated biological impacts at some plots, both in terms of differences from a control plot during the period of discharge operations, and in terms of the recovery period following the perturbations. Evaluation of the significance of potential effects involves the following factors: • Definition of a specific spatial and temporal reference frame (e.g., What is the natural variation? Is 1 year to be considered a "long" or "short" time? Is 50 m to be considered a "large" or "trivial" distance? • Consideration of rare or unique resources and particularly sensitive biotic assemblages. • Consideration of the potential for long term, cumulative effects. Some of these aspects are clearly beyond the scope of this modeling efforts (e.g., the model does not simulate the long term fate of resuspended material). The model does, however, contain an internal "reference frame" by comparison to simulated behavior at a control plot. The model, in general, simulates substantial "natural" variation at the reference or control plots, both over time, due to sediment disturbance events in medium to high energy environments, and over space, due to geographically varying conditions, such as water depth and current regime.

Report↗
USGS2000

Evaluation of planar halogenated and polycyclic aromatic hydrocarbons in estuarine sediments using ethoxyresorufin-O-deethylase induction of H4IIE cells

Polycyclic aromatic hydrocarbons (PAHs) and planar halogenated hydrocarbons (PHHs), including polychlorinated dibenzo- p -dioxins (PCDDs), dibenzofurans (PCDFs), and biphenyls (PCBs) were determined in fractionated sediment extracts from the Hudson-Raritan estuary and Newark Bay, New Jersey, USA, as part of a comprehensive risk assessment. Contributions of PCDDs/PCDFs, PCBs, and PAHs to the total toxic equivalents (TEQs) were measured using an H4IIE bioassay and calculated from instrumentally determined concentrations using international toxic equivalency factors. The H4IIE TEQs of whole and fractionated extracts were compared to calculated TEQs to investigate the applicability of the bioassay approach for evaluating 7-ethoxyresorufin- O -deethylease induction by PHHs and PAHs present together in complex mixtures. Although 2,3,7,8-tetrachloro-dibenzo- p -dioxin contributed from 41 to 79% of the calculated TEQs from PCDDs/PCDFs and planar PCBs in all sediments sampled, the PAH-containing fractions accounted for >80% of the total TEQs determined either instrumentally or by bioassay. Calculated TEQs from PAHs, based on reported toxic equivalency factors for only seven PAHs, were severalfold greater than the bioassay-derived TEQs of PAH-only fractions of the sediment extracts. Significant correlations were observed between bioassay and instrumentally determined toxic equivalents in the more purified fractions but not in fractions only purified by size-exclusion or argentate chromatographies alone.

Environmental Toxicology and Chemistry↗
USGS2022

Great expectations for earthquake early warnings on the United States West Coast

In October 2019, California became the first state in the United States to fully activate a public earthquake early warning system—ShakeAlert®—managed by the U.S. Geological Survey. The system was subsequently rolled out in March 2021 in Oregon and May 2021 in Washington. Earthquake early warning (EEW) systems can provide seconds of notice to people and technological systems that shaking is imminent, but their effectiveness depends on recipients’ expectations and actions as well as technical performance. To better understand these dependencies, we surveyed representative samples of adults in California (N = 1219), Oregon (N = 1020), and Washington (N = 1037) in February 2021. Most respondents had experienced earthquakes, but few had lived through violent shaking; most had not followed protective action guidance to Drop, Cover, and Hold On (DCHO) in earthquakes; and most reported no personal or social harm from prior earthquakes. Nevertheless, expectations and perceived usefulness of EEW were high, and higher still for those who expected alerts to be accurate and easy to use, expressed tolerance of missed and erroneous warnings, and expected to be affected by a damaging earthquake in their lifetime. Results suggest opportunities to better align public preferences and expectations with ShakeAlert operations. For example, some respondents preferred lower alerting thresholds than those proposed by government and scientists. Moreover, reported tolerance of warning errors was widespread, but respondents wanted explanations quickly, suggesting a need to further develop post-alert messaging. Findings from this study should be informative for future research on the co-evolution of experiences and expectations with EEW systems.

California, Oregon, Washington↗
USGS2021

SSA task force on diversity, equity, and inclusion: Toward a changing, inclusive future in earthquake science

In the United States, a wide variety of studies show that the geoscience community does not reflect the broader societal makeup (e.g., Velasco and Jaurrieta de Velasco, 2010; Dutt, 2020; Howley, 2020). In fact, only about 10% of all Science, Technology, Engineering, and Mathematics (STEM) Ph.D. degrees are awarded to people of color, although they represent more than a third of the population (Dutt, 2020). These numbers have changed little over the past 40 yr (e.g., Bernard and Cooperdock, 2018; Dutt, 2020). Recent events in the United States have again raised awareness of this discrepancy in many fields; similar divergences may be present throughout the world. This discrepancy represents a tremendous loss of talent and contributes to ongoing bias and racism. As the premier international seismological professional society, Seismological Society of America (SSA) is committed, as stated in its ethics policy, to “freedom and transparency in research and education, which should be conducted in a supportive, inclusive, and respectful environment, free from any discrimination, harassment or bullying.” SSA convened the SSA Diversity, Equity, and Inclusion (DEI) Task Force to identify specific and actionable recommendations that the Society can take to start to address the significant issues surrounding the systematic underrepresentation in seismology and related geosciences by marginalized groups on both short- and long- term time horizons. This report, written by the Task Force, outlines actionable items that SSA can undertake as well as a framework by which progress can be measured to help guide the Society and its members forward as a community to make earthquake science more just, equitable, diverse, and inclusive. We recognize that this document cannot address the multiple individual, cultural, and historical barriers that some members may face (e.g., Núñez et al., 2020), nor can it reflect all possibilities or all recommended practices (e.g., Ali et al., 2021). What this document does represent, however, is the beginning of what we hope to be a significant cultural change in how we conduct earthquake science as an organization, not the end of this process as the Task Force fulfills its charter. As a note, some of the recommendations that the Task Force has outlined here have already been implemented within SSA or are actively being considered. For completeness, we still include them here because this highlights the commitment to change that SSA leadership has embraced.

Seismological Research Letters↗
USGS2000

Factors affecting nutrient trends in major rivers of the Chesapeake Bay Watershed

Trends in nutrient loads and flow-adjusted concentrations in the major rivers entering Chesapeake Bay were computed on the basis of water-quality data collected between 1985 and 1998 at 29 monitoring stations in the Susquehanna, Potomac, James, Rappahannock, York, Patuxent, and Choptank River Basins. Two computer models-the Chesapeake Bay Watershed Model (WSM) and the U.S. Geological Survey's 'Spatially Referenced Regressions on Watershed attributes' (SPARROW) Model-were used to help explain the major factors affecting the trends. Results from WSM simulations provided information on temporal changes in contributions from major nutrient sources, and results from SPARROW model simulations provided spatial detail on the distribution of nutrient yields in these basins. Additional data on nutrient sources, basin characteristics, implementation of management practices, and ground-water inputs to surface water were analyzed to help explain the trends. The major factors affecting the trends were changes in nutrient sources and natural variations in streamflow. The dominant source of nitrogen and phosphorus from 1985 to 1998 in six of the seven tributary basins to Chesapeake Bay was determined to be agriculture. Because of the predominance of agricultural inputs, changes in agricultural nutrient sources such as manure and fertilizer, combined with decreases in agricultural acreage and implementation of best management practices (BMPs), had the greatest impact on the trends in flow-adjusted nutrient concentrations. Urban acreage and population, however, were noted to be increasing throughout the Chesapeake Bay Watershed, and as a result, delivered loads of nutrients from urban areas increased during the study period. Overall, agricultural nutrient management, in combination with load decreases from point sources due to facility upgrades and the phosphate detergent ban, led to downward trends in flow-adjusted nutrient concentrations atmany of the monitoring stations in the watershed. The loads of nutrients, however, were not reduced significantly at most of the monitoring stations. This is due primarily to higher streamflow in the latter years of the monitoring period, which led to higher loading in those years. Results of this study indicate a need for more detailed information on BMP effectiveness under a full range of hydrologic conditions and in different areas of the watershed; an internally consistent fertilizer data set; greater consideration of the effects of watershed processes on nutrient transport; a refinement of current modeling efforts; and an expansion of the non-tidal monitoring network in the Chesapeake Bay Watershed.

Delaware, Maryland, New York, Pennsylvania, Virgin↗