Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Field Guides”

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 595 records · Page 33Linked to original sources

Semipermanent GPS (SPGPS) as a volcano monitoring tool: Rationale, method, and applications

Semipermanent GPS (SPGPS) is an alternative to conventional campaign or survey-mode GPS (SGPS) and to continuous GPS (CGPS) that offers several advantages for monitoring ground deformation. Unlike CGPS installations, SPGPS stations can be deployed quickly in response to changing volcanic conditions or earthquake activity such as a swarm or aftershock sequence. SPGPS networks can be more focused or more extensive than CGPS installations, because SPGPS equipment can be moved from station to station quickly to increase the total number of stations observed in a given time period. SPGPS networks are less intrusive on the landscape than CGPS installations, which makes it easier to satisfy land-use restrictions in ecologically sensitive areas. SPGPS observations are preferred over SGPS measurements because they provide better precision with only a modest increase in the amount of time, equipment, and personnel required in the field. We describe three applications of the SPGPS method that demonstrate its utility and flexibility. At the Yellowstone caldera, Wyoming, a 9-station SPGPS network serves to densify larger preexisting networks of CGPS and SGPS stations. At the Three Sisters volcanic center, Oregon, a 14-station SPGPS network complements an SGPS network and extends the geographic coverage provided by 3 CGPS stations permitted under wilderness land-use restrictions. In the Basin and Range province in northwest Nevada, a 6-station SPGPS network has been established in response to a prolonged earthquake swarm in an area with only sparse preexisting geodetic coverage. At Three Sisters, the estimated precision of station velocities based on annual ~ 3 month summertime SPGPS occupations from 2009 to 2015 is approximately half that for nearby CGPS stations. Conversely, SPGPS-derived station velocities are about twice as precise as those based on annual ~ 1 week SGPS measurements. After 5 years of SPGPS observations at Three Sisters, the precision of velocity determinations is estimated to be 0.5 mm/yr in longitude, 0.6 mm/yr in latitude, and 0.8 mm/yr in height. We conclude that an optimal approach to monitoring volcano deformation includes complementary CGPS and SPGPS networks, periodic InSAR observations, and measurements from in situ borehole sensors such as tiltmeters or strainmeters. This comprehensive approach provides the spatial and temporal detail necessary to adequately characterize a complex and evolving deformation pattern. Such information is essential to multi-parameter models of magmatic or tectonic processes that can help to guide research efforts, and also to inform hazards assessments and land-use planning decisions.

Journal of Volcanology and Geothermal Research↗

Physics to fish—Understanding the factors that create and sustain native fish habitat in the San Francisco Estuary

Executive Summary The Bureau of Reclamation (Reclamation) operates the Central Valley Project (CVP), one of the nation’s largest water projects. Reclamation has an ongoing need to improve the scientific basis for adaptive management of the CVP and, by extension, joint operations with California’s State Water Project. The U.S. Geological Survey (USGS) works cooperatively with the Bureau of Reclamation to provide scientific support for the management of Reclamation’s CVP project. Major habitat restoration efforts and a new water-diversion point are planned to benefit delta smelt ( Hypomesus transpacificus ) and other species of concern while ensuring the reliability of water supply. In addition, various flow actions and management activities have been identified as possible methods to increase populations of delta smelt and salmonid ( Oncorhynchus spp.) runs of concern. The overarching goal of this cooperative project was to provide Reclamation with the scientific information needed to evaluate the efficacy of ongoing and future adaptive management actions and to improve the scientific basis for more flexible CVP operations that would achieve water-supply reliability and fish protection. The research and monitoring described in this report comprises the period 2015–19 and focuses on management issues related to native fish species of concern, especially delta smelt. Conserving the delta smelt population while providing a reliable water supply is a primary management and policy issue in California. Our approach for this cooperative project is based on the “physics to fish” concept, the idea that high-quality habitat is generated and sustained by the interaction between physical processes and the landscape. These interactions create a template for chemical and biological processes that can change across a variety of spatial and temporal scales. Following this concept, this project (hereafter referred to as “the physics to fish project”) included monitoring and studies of water flows, sediments, water quality, and invertebrate and fish dynamics across a range of spatial and temporal scales and in regions relevant to resource managers tasked with managing water supplies and ecosystem health in the San Francisco Estuary. The intent of this approach was to document the habitat conditions, important processes, and interactions among them that create high-quality habitat for native fishes so that the likely effects of future management actions (for example, habitat restoration) can be objectively assessed at the local (site-specific), regional (within subregions of the estuary), and landscape (across the entire estuary and beyond) scales. Hydrodynamically, the upper estuary (landward of Carquinez Strait) is characterized by a fixed volume of tidally exchanged water (for example, tidal prism) that interacts with the existing channel network and bathymetry to create regions with differing hydrodynamics. Our results indicate that careful study of construction or reoperation of existing infrastructure to perform management actions can help (1) improve the accuracy of hydrodynamic models; (2) further understanding of ecological effects; and (3) enhance abilities to predict ecological outcomes. At the local scale, we developed a new concept called the Lagrangian to Eulerian (LE) ratio that can be used as a tool for understanding the importance of various hydrodynamic processes in specific channels or channel networks and for forecasting transport dynamics. Channels with LE ratios<1 in a channel network or in a dead-end slough are hydrodynamically able to develop an exchange zone between two parcels of water that may have different chemical and physical properties. In a dead-end channel, there is a landward region with long residence time (no-exchange zone) and a seaward region with short residence time (high-exchange zone) that are well mixed with seaward waters. At the transition (exchange zone) between the high and no-exchange regions, a gradient will form in water-quality constituents that differ in concentration between the landward and seaward waters. Turbidity affects fish habitat and has declined through time in the San Francisco Estuary. Average turbidity across the Sacramento–San Joaquin Delta (hereafter referred to as “the Delta”) is dependent on annual hydrology. In dry years, the region around Cache Slough (known regionally as the “Cache Slough Complex”) in the northern Delta is generally more turbid than Suisun Bay and the lower Sacramento River. When the Yolo By-Pass (known regionally as “Yolo Bypass”), a large flood bypass that runs parallel to the Sacramento River in the northern Delta, is not flooding and river flows are lower, sediment is usually transported into the Cache Slough Complex because flood tides dominate ebb tides, resulting in transport of suspended sediment from seaward areas of the upper estuary into the Cache Slough Complex. These hydrodynamic conditions also favor the formation of turbidity maximums (TMs) in the Cache Slough Complex. The TMs are areas of higher suspended-sediment concentration, providing higher-turbidity habitat favored by some fishes, including delta smelt, and they can also concentrate other constituents, including phytoplankton and organic carbon that can be important in food webs. Pelagic primary production by phytoplankton is the basis for Delta food webs supporting pelagic fishes such as delta smelt; however, phytoplankton abundance in the Delta has declined during recent decades. We examined how nutrients, hydrodynamics, and other factors affect phytoplankton blooms. Based on our results, we developed three new concepts of phytoplankton bloom formation in the Delta, each associated with a distinct set of hydrologic conditions. First, productivity cascades highlighted how local processes can contribute to phytoplankton blooms observed at the regional scale. Second, we observed phytoplankton blooms in the upper San Francisco Estuary that were associated with transport out of Yolo By-Pass (transport blooms). Third, we also documented a series of phytoplankton blooms that were in the confluence area at the landward edge of Suisun Bay. The conditions leading to creation of confluence phytoplankton blooms are not yet understood, but the confluence region connects the Cache Slough Complex with Suisun Marsh. Therefore, blooms in this area have the potential to spread to large areas of the Delta. At the landscape scale, the distribution of the invasive clams ( Potamocorbula amurensis and Corbicula fluminea , hereafter referred to as “ Corbicula ”) is driven by salinity. At smaller spatial scales, the distribution of either species is sensitive to multiple factors affecting survival and reproduction, complicating efforts to predict distribution and abundance without considering local-scale conditions across the area of interest. In the Cache Slough Complex, the area landward of the exchange zone in regions with LE ratio<1 were characterized by low abundances of Corbicula probably because recruits from seaward areas are not transported past the exchange zone and because there are no landward tributaries with adult Corbicula to provide an upstream source of recruits. Corbicula biomass was highest near or downstream from the exchange zone consistent with Corbicula grazing on phytoplankton produced in the exchange zone or transported from the no-exchange zone. The severity of Corbicula grazing could be reduced by manipulating the hydrodynamic characteristics of waterways; however, the beneficial and harmful effects on the organisms meant to benefit from increased phytoplankton production, including zooplankton and fish species of concern, should be thoroughly examined before manipulating hydrodynamic characteristics. The distribution of fishes at the landscape scale is generally driven by the position of the salinity field in the estuary. The physics to fish project compared distributions of fishes at Ryer Island, a tidal wetland in Suisun Bay and a region of variable salinity, with fish distributions at the Cache Slough Complex, a freshwater region. At Ryer Island, there was an absence of freshwater invasive species and an abundance of native species, such as Sacramento splittail ( Pogonichthys macrolepidotus ), tule perch ( Hysterocarpus traskii ), and Sacramento pikeminnow ( Ptychocheilus grandis ). The native species were almost exclusively captured in wetland and nearshore shallow-water habitat regardless of water-quality conditions. In the Cache Slough Complex, our regional scale objective was to elucidate how hydrodynamic-physical habitat interactions drive fish-community structure. Our studies showed that dendritic channel systems were better able to support native species, while intertidal habitats supported those species best able to exploit the transient character of the habitat. Habitats upstream from the exchange zone were especially important in supporting high numbers of native fishes relative to within or downstream from the exchange zone. Many of the native species were associated with tidal marsh in the no-exchange zone. More pelagic-oriented, mobile species, such as Striped Bass ( Morone saxatilis ), threadfin shad ( Dorosoma petenense ), and Sacramento pikeminnow, were more affected by water-quality conditions, such as turbidity. The physics to fish concept developed in this project provides a framework for designing individual projects and for considering the cumulative effects of multiple projects in a region, using the LE ratio as a guiding metric. The physics to fish concept may also provide a suitable framework for coordinating management actions. Tidal wetlands can function in several ways in the hydrodynamic framework. Relatively small tidal wetlands with short channel networks and with LE ratios>1 are not able to maintain a landward no-exchange zone or an exchange zone. This likely means that any contributions to pelagic food webs would be limited to resources derived from wetland vegetation, which can include dissolved and particulate organic matter (detritus) and populations of consumers that can increase in abundance based on those resources. The fate of the contributed production from these channels depends on the characteristics of the receiving waters seaward of the tidal wetland. If these channels join a large system such as Suisun Bay, then any contribution is likely to be rapidly dispersed in the larger volume; however, the channel junction might provide a focal point for consumers, such as fishes, to congregate and feed on material leaving the wetland on ebb tides before it is dispersed in the larger volume. Fishes might also access these resources by entering the wetland. The physics to fish project has established a foundation and several new concepts for understanding how habitat restoration can benefit native fish populations at the local and regional levels. Many of the ideas regarding habitat restoration and channel modifications outlined in this report could help guide management actions that could improve conditions for native fishes at little or no water cost beyond water already dedicated to other management actions. A complete list of products originating from this work is provided in appendix 1.

California↗

A comparison of macroinvertebrate and habitat methods of data collection in the Little Colorado River Watershed, Arizona 2007

The Arizona Department of Environmental Quality (ADEQ) and the U.S. Environmental Protection Agency (USEPA) Ecological Monitoring and Assessment Program (EMAP), use different field methods for collecting macroinvertebrate samples and habitat data for bioassessment purposes. Arizona’s Biocriteria index was developed using a riffle habitat sampling methodology, whereas the EMAP method employs a multi-habitat sampling protocol. There was a need to demonstrate comparability of these different bioassessment methodologies to allow use of the EMAP multi-habitat protocol for both statewide probabilistic assessments for integration of the EMAP data into the national (305b) assessment and for targeted in-state bioassessments for 303d determinations of standards violations and impaired aquatic life conditions. The purpose of this study was to evaluate whether the two methods yield similar bioassessment results, such that the data could be used interchangeably in water quality assessments. In this Regional EMAP grant funded project, a probabilistic survey of 30 sites in the Little Colorado River basin was conducted in the spring of 2007. Macroinvertebrate and habitat data were collected using both ADEQ and EMAP sampling methods, from adjacent reaches within these stream channels. All analyses indicated that the two macroinvertebrate sampling methods were significantly correlated. ADEQ and EMAP samples were classified into the same scoring categories (meeting, inconclusive, violating the biocriteria standard) 82% of the time. When the ADEQ-IBI was applied to both the ADEQ and EMAP taxa lists, the resulting IBI scores were significantly correlated (r=0.91), even though only 4 of the 7 metrics in the IBI were significantly correlated. The IBI scores from both methods were significantly correlated to the percent of riffle habitat, even though the average percent riffle habitat was only 30% of the stream reach. Multivariate analyses found that the percent riffle was an important attribute for both datasets in classifying IBI scores into assessment categories. Habitat measurements generated from EMAP and ADEQ methods were also significantly correlated; 13 of 16 habitat measures were significantly correlated (p<0.01). The visual-based percentage estimates of percent riffle and pool habitats, vegetative cover and percent canopy cover, and substrate measurements of percent fine substrate and embeddedness were all remarkably similar, given the different field methods used. A multivariate analysis identified substrate and flow conditions, as well as canopy cover as important combinations of habitat attributes affecting both IBI scores. These results indicate that similar habitat measures can be obtained using two different field sampling protocols. In addition, similar combinations of these habitat parameters were important to macroinvertebrate community condition in multivariate analyses of both ADEQ and EMAP datasets. These results indicate the two sampling methods for macroinvertebrates and habitat data were very similar in terms of bioassessment results and stressors. While the bioassessment category was not identical for all sites, overall the assessments were significantly correlated, providing similar bioassessment results for the cold water streams used in this study. The findings of this study indicate that ADEQ can utilize either a riffle-based sampling methodology or a multi-habitat sampling approach in cold water streams as both yield similar results relative to the macroinvertebrate assemblage. These results will allow for use of either macroinvertebrate dataset to determine water quality standards compliance with the ADEQ Indexes of Biological Integrity, for which threshold values were just recently placed into the Arizona Surface Water Quality Standards. While this survey did not include warm water desert streams of Arizona, we would predict that EMAP and ADEQ sampling methodologies would provide similar bioassessment results and would not be significantly different, as we have found that the percent riffle habitat in cold and warm water perennial, wadeable streams is not significantly different. However, a comparison study of sampling methodologies in warm water streams should be conducted to confirm the predicted similarity of bioassessment results. ADEQ will continue to implement a monitoring strategy that includes probabilistic monitoring for a statewide ecological assessment of stream conditions. Conclusions from this study will guide decisions regarding the most appropriate sampling methods for future probabilistic monitoring sample plans.

Arizona↗

Chemical analyses in the World Coal Quality Inventory

The main objective of the World Coal Quality Inventory (WoCQI) was to collect and analyze a global set of samples of mined coal during a time period from about 1995 to 2006 (Finkelman and Lovern, 2001). Coal samples were collected by foreign collaborators and submitted to country specialists in the U.S. Geological Survey (USGS) Energy Program. However, samples from certain countries, such as Afghanistan, India, and Kyrgyzstan, were collected collaboratively in the field with USGS personnel. Samples were subsequently analyzed at two laboratories: the USGS Inorganic Geochemistry Laboratory located in Denver, CO and a commercial laboratory (Geochemical Testing, Inc.) located in Somerset, PA. Thus the dataset, which is in Excel (2003) format and includes 1,580 samples from 57 countries, does not have the inter-laboratory variability that is present in many compilations. Major-, minor-, and trace-element analyses from the USGS laboratory, calculated to a consistent analytical basis (dry, whole-coal) and presented with available sample identification information, are sorted alphabetically by country name. About 70 percent of the samples also have data from the commercial laboratory, which are presented on an as-received analytical basis. The USGS initiated a laboratory review of quality assurance in 2008, covering quality control and methodology used in inorganic chemical analyses of coal, coal power plant ash, water, and sediment samples. This quality control review found that data generated by the USGS Inorganic Geochemistry Laboratory from 1996 through 2006 were characterized by quality practices that did not meet USGS requirements commonly in use at the time. The most serious shortcomings were (1) the adjustment of raw sample data to standards when the instrument values for those standards exceeded acceptable limits or (2) the insufficient use of multiple standards to provide adequate quality assurance. In general, adjustment of raw data to account for instrument 'drift' is an acceptable practice within strictly defined limits. During the denoted period, USGS required that the maximum adjustment of instrument values, guided by calibration standards, was not allowed to exceed 10 percent. However, in some cases, the Inorganic Geochemistry Laboratory released data that were adjusted by more than 10 percent and (or) were not constrained by an adequate number of control standards. Original instrument values no longer exist for about 80 percent of the analyses during this period; therefore, the acceptability of drift corrections for most of the samples analyzed cannot be determined. For these reasons, the WoCQI data from the USGS Inorganic Geochemistry Laboratory should be used with care. For more information, individuals may contact laboratory management at EnergyLabs@usgs.gov with specific questions about particular datasets or analytical attributes. Standard USGS sampling methods were provided and recommended to collaborators, but the analyzed samples may or may not be representative of their locale; for some samples, only limited information is available concerning sample provenance. Single samples cannot represent spatial or temporal variability within a coal area. Geochemical datasets of U.S. coals can be found in the COALQUAL database (Bragg and others, 1997) and the National Coal Quality Inventory (Hatch and others, 2006), as only non-U.S. sample data are presented in the WoCQI. Although the WoCQI does not contain worldwide coverage of coal deposits, it is truly a unique and valuable compilation. The information in the WoCQI should prove useful for identifying possible areas for future global coal research.

Open-File Report↗

Volcano hazards from Mount Rainier, Washington

Mount Rainier—at 4393 meters (14,410 feet) the highest peak in the Cascade Range—is a dormant volcano whose load of glacier ice exceeds that of any other mountain in the conterminous United States. This tremendous mass of rock and ice, in combination with great topographic relief, poses a variety of geologic hazards, both during inevitable future eruptions and during the intervening periods of repose. The volcano’s past behavior is the best guide to possible future hazards. The written history of Mount Rainier encompasses the period since about A.D. 1820, during which time one or two small eruptions, many small debris flows, and several small debris avalanches have occurred. This time interval is far too brief to serve as a basis for estimating the future behavior of a volcano that is several hundreds of thousands of years old. Fortunately, prehistoric deposits record the types, magnitudes, and frequencies of past events, and show which areas were affected by them. At Mount Rainier, as at other Cascade volcanoes, deposits produced since the latest ice age (approximately during the past 10,000 years) are well preserved. Studies of these deposits reveal that we should anticipate potential hazards from some phenomena that only occur during eruptions and from others that may occur without eruptive activity. Tephra falls, pyroclastic flows and pyroclastic surges, ballistic projectiles, and lava flows occur only during eruptions. Debris avalanches, debris flows, and floods commonly accompany eruptions, but can also occur during dormant periods. This report (1) explains the various types of hazardous geologic phenomena that could occur at Mount Rainier, (2) shows areas that are most likely to be affected by the different phenomena, (3) estimates the likelihood that the areas will be affected, and (4) recommends actions that can be taken to protect lives and property. It builds upon and revises a similar document prepared by D.R. Crandell in 1973. Our revision was motivated by the availability of new information about Mount Rainier’s geologic history, by advances in the field of volcanology, and by the need to assess hazards in a more quantitative manner than in Crandell’s pioneering report.

Washington↗

Poroelastic rebound along the Landers 1992 earthquake surface rupture

Maps of surface displacement following the 1992 Landers, California, earthquake, generated by interferometric processing of ERS-1 synthetic aperture radar (SAR) images, reveal effects of various postseismic deformation processes along the 1992 surface rupture. The large-scale pattern of the postseismic displacement field includes large lobes, mostly visible on the west side of the fault, comparable in shape with the lobes observed in the coseismic displacement field. This pattern and the steep displacement gradient observed near the Emerson-Camp Rock fault cannot be simply explained by afterslip on deep sections of the 1992 rupture. Models show that horizontal slip occurring on a buried dislocation in a Poisson's material produces a characteristic quadripole pattern in the surface displacement field with several centimeters of vertical motion at distances of 10–20 km from the fault, yet this pattern is not observed in the postseismic interferograms. As previously proposed to explain local strain in the fault step overs [ Peltzer et al ., 1996b], we argue that poroelastic rebound caused by pore fluid flow may also occur over greater distances from the fault, compensating the vertical ground shift produced by fault afterslip. Such a rebound is explained by the gradual change of the crustal rocks' Poisson's ratio value from undrained (coseismic) to drained (postseismic) conditions as pore pressure gradients produced by the earthquake dissipate. Using the Poisson's ratio values of 0.27 and 0.31 for the drained and undrained crustal rocks, respectively, elastic dislocation models show that the combined contributions of afterslip on deep sections of the fault and poroelastic rebound can account for the range change observed in the SAR data and the horizontal displacement measured at Global Positioning System (GPS) sites along a 60-km-long transect across the Emerson fault [ Savage and Svarc , 1997]. Using a detailed surface slip distribution on the Homestead Valley, Kickapoo, and Johnson Valley faults, we modeled the poroelastic rebound in the Homestead Valley pull apart. A Poisson's ratio value of 0.35 for the undrained gouge rocks in the fault zone is required to account for the observed surface uplift in the 3.5 years following the earthquake. This large value implies a seismic velocity ratio V p /V s of 2.1, consistent with the observed low V s values of fault zone guided waves at shallow depth [ Li et al ., 1997]. The SAR data also reveal postseismic creep along shallow patches of the Eureka Peak and Burnt Mountain faults with a characteristic decay time of 0.8 years. Coseismic, dilatant hardening (locking process) followed by post-seismic, pore pressure controlled fault creep provide a plausible mechanism to account for the decay time of the observed slip rate along this section of the fault.

Journal of Geophysical Research B: Solid Earth↗

Development of relations of stream stage to channel geometry and discharge for stream segments simulated with Hydrologic Simulation Program-Fortran (HSPF), Chesapeake Bay Watershed and adjacent parts of Virginia, Maryland, and Delaware

The U.S. Geological Survey (USGS), U.S. Environmental Protection Agency (USEPA), Chesapeake Bay Program (CBP), Interstate Commission for the Potomac River Basin (ICPRB), Maryland Department of the Environment (MDE), Virginia Department of Conservation and Recreation (VADCR), and University of Maryland (UMD) are collaborating to improve the resolution of the Chesapeake Bay Regional Watershed Model (CBRWM). This watershed model uses the Hydrologic Simulation Program-Fortran (HSPF) to simulate the fate and transport of nutrients and sediment throughout the Chesapeake Bay watershed and extended areas of Virginia, Maryland, and Delaware. Information from the CBRWM is used by the CBP and other watershed managers to assess the effectiveness of water-quality improvement efforts as well as guide future management activities. A critical step in the improvement of the CBRWM framework was the development of an HSPF function table (FTABLE) for each represented stream channel. The FTABLE is used to relate stage (water depth) in a particular stream channel to associated channel surface area, channel volume, and discharge (streamflow). The primary tool used to generate an FTABLE for each stream channel is the XSECT program, a computer program that requires nine input variables used to represent channel morphology. These input variables are reach length, upstream and downstream elevation, channel bottom width, channel bankfull width, channel bankfull stage, slope of the floodplain, and Manning's roughness coefficient for the channel and floodplain. For the purpose of this study, the nine input variables were grouped into three categories: channel geometry, Manning's roughness coefficient, and channel and floodplain slope. Values of channel geometry for every stream segment represented in CBRWM were obtained by first developing regional regression models that relate basin drainage area to observed values of bankfull width, bankfull depth, and bottom width at each of the 290 USGS streamflow-gaging stations included in the areal extent of the model. These regression models were developed on the basis of data from stations in four physiographic provinces (Appalachian Plateaus, Valley and Ridge, Piedmont, and Coastal Plain) and were used to predict channel geometry for all 738 stream segments in the modeled area from associated basin drainage area. Manning's roughness coefficient for the channel and floodplain was represented in the XSECT program in two forms. First, all available field-estimated values of roughness were compiled for gaging stations in each physiographic province. The median of field-estimated values of channel and floodplain roughness for each physiographic province was applied to all respective stream segments. The second representation of Manning's roughness coefficient was to allow roughness to vary with channel depth. Roughness was estimated at each gaging station for each 1-foot depth interval. Median values of roughness were calculated for each 1-foot depth interval for all stations in each physiographic province. Channel and floodplain slope were determined for every stream segment in CBRWM using the USGS National Elevation Dataset. Function tables were generated by the XSECT program using values of channel geometry, channel and floodplain roughness, and channel and floodplain slope. The FTABLEs for each of the 290 USGS streamflow-gaging stations were evaluated by comparing observed discharge to the XSECT-derived discharge. Function table stream discharge derived using depth-varying roughness was found to be more representative of and statistically indistinguishable from values of observed stream discharge. Additionally, results of regression analysis showed that XSECT-derived discharge accounted for approximately 90 percent of the variability associated with observed discharge in each of the four physiographic provinces. The results of this study indicate that the methodology developed to generate FTABLEs for every s

Chesapeake Bay Watershed↗

Volcano hazards from Mount Rainier, Washington, revised 1998

Mount Rainier—at 4393 meters (14,410 feet) the highest peak in the Cascade Range—is a dormant volcano whose load of glacier ice exceeds that of any other mountain in the conterminous United States. This tremendous mass of rock and ice, in combination with great topographic relief, poses a variety of geologic hazards, both during inevitable future eruptions and during the intervening periods of repose. The volcano’s past behavior is the best guide to possible future hazards. The written history of Mount Rainier encompasses the period since about A.D. 1820, during which time one or two small eruptions, several small debris avalanches, and many small lahars (debris flows originating on a volcano) have occurred. This time interval is far too brief to serve as a basis for estimating the future behavior of a volcano that is several hundreds of thousands of years old. Fortunately, prehistoric deposits record the types, magnitudes, and frequencies of past events, and show which areas were affected by them. At Mount Rainier, as at other Cascade volcanoes, deposits produced since the latest ice age (approximately during the past 10,000 years) are well preserved. Studies of these deposits reveal that we should anticipate potential hazards from some phenomena that only occur during eruptions and from others that may occur without eruptive activity. Tephra falls, pyroclastic flows and pyroclastic surges, ballistic projectiles, and lava flows occur only during eruptions. Debris avalanches, lahars, and floods commonly accompany eruptions, but can also occur during dormant periods. This report (1) explains the various types of hazardous geologic phenomena that could occur at Mount Rainier, (2) shows areas that are most likely to be affected by the different phenomena, (3) estimates the likelihood that the areas will be affected, and (4) recommends actions that can be taken to protect lives and property. It builds upon and revises a similar document prepared by D.R. Crandell in 1973. Our revision was motivated by the availability of new information about Mount Rainier’s geologic history, by advances in the field of volcanology, and by the need to assess hazards in a more quantitative manner than in Crandell’s pioneering report.

Washington↗

BioData: A national aquatic bioassessment database

BioData is a U.S. Geological Survey (USGS) web-enabled database that for the first time provides for the capture, curation, integration, and delivery of bioassessment data collected by local, regional, and national USGS projects. BioData offers field biologists advanced capabilities for entering, editing, and reviewing the macroinvertebrate, algae, fish, and supporting habitat data from rivers and streams. It offers data archival and curation capabilities that protect and maintain data for the long term. BioData provides the Federal, State, and local governments, as well as the scientific community, resource managers, the private sector, and the public with easy access to tens of thousands of samples collected nationwide from thousands of stream and river sites. BioData also provides the USGS with centralized data storage for delivering data to other systems and applications through automated web services. BioData allows users to combine data sets of known quality from different projects in various locations over time. It provides a nationally aggregated database for users to leverage data from many independent projects that, until now, was not feasible at this scale. For example, from 1991 to 2011, the USGS Idaho Water Science Center collected more than 816 bioassessment samples from 63 sites for the National Water Quality Assessment (NAWQA) Program and more than 477 samples from 39 sites for a cooperative USGS and State of Idaho Statewide Water Quality Network (fig. 1). Using BioData, 20 years of samples collected for both of these projects can be combined for analysis. BioData delivers all of the data using current taxonomic nomenclature, thus relieving users of the difficult and time-consuming task of harmonizing taxonomy among samples collected during different time periods. Fish data are reported using the Integrated Taxonomic Information Service (ITIS) Taxonomic Serial Numbers (TSN's). A simple web-data input interface and self-guided, public data-retrieval web site provides access to bioassessment data. BioData currently accepts data collected using two national protocols: (1) NAWQA and (2) U.S. Environmental Protection Agency (USEPA) National Rivers and Streams Assessment (NRSA). Additional collection protocols are planned for future versions.

Fact Sheet↗

Uranium-bearing lignite and its relation to the White River and Arikaree formations in northwestern South Dakota and adjacent states

In northwestern South Dakota and adjacent areas uranium-bearing lignite beds occur at many horizons in the Hell Creek formation of late Cretaceous age and the overlying Ludlow, Tongue River, and Sentinel Butte members of the Fort Union formation of Paleocene age. Uranium analyses of 275 surface and auger samples and about 1,000 core samples show that many of the lignite beds contain 0. 005 to 0. 02 percent uranium with concentrations of 0. 05 to 0.10 percent uranium in the lignite ash. Analytical data indicate that the region contains an aggregate of at least 47,500, 000 tons of lignite with an average grade of slightly more than .0. 008 percent containing 3, 900 tons of uranium. Almost a fifth of the estimated reserves are adapted to strip mining and are in beds averaging about 4 feet in thickness. Uranium concentrations of this magnitude in lignite indicate that these deposits upon the development of proper utilization techniques and processes may be an important future source of uranium. Recent discoveries of ore-grade deposits of autunite-bearing lignite and secondary uranium minerals in carbonaceous sandstone at Cave Hills and Slim Buttes indicate that northwestern South Dakota and adjacent areas may containimportant reserves of uranium-ore. The stratigraphic units containing the uraniferous lignite beds have a combined thickness of about 1, 500 feet and are unconformably overlapped by 300 feet or more of tuffaceous sandstone and bentonitic claystone of the White River and Arikaree formations of Oligocene and Miocene age. The stratigraphically highest lignite beds in the local sequence have the greatest concentration of uranium,, and the uranium content is greatest at the top of thick lignite beds, diminishing progressively downward to a vanishing point in their lower parts. Variations in permeability of the rock overlying the mineralized lignite beds seem to be reflected in the intensity of uranium mineralization. Most of the known uranium-bearing lignite deposits in the region are closely overlain by the White River and Arikaree formations. Field evidence indicates that the uranium in the lignite is independent of the age of the formation in which the lignite occurs and that the uranium has been concentrated by downward and laterally moving ground water from the overlying mildly radioactive, tuffaceous rocks. The White River and Arikaree formations have about 1.2 times more uranium than the average sedimentary rock. The uranium content of spring water from these formations is 30 times as great or greater than that of normal ground water. Field relations suggest that the uranium is of secondary origin and has been introduced subsequent to the accumulation and marked regional warping of the lignite beds and associated rocks. During the transportation of the uranium-bearing volcanic materials from their place of origin to the site of deposition in South Dakota, the uranium is believed to have been held as a finely disseminated constituent in volcanic ash. Subsequent release or displacement of the uranium is thought to have resulted from weathering and chemical breakdown, thus freeing and allowing the transportation of the uranium, probably in ionic form, to the lignite by downward and laterally moving groundwaters 0 Gar bonaceous material in the path of these uranium-bearing waters is believed to have acted as a receptor that extracted the uranium as a result of an ion exchange mechanism or by the formation of organo-metallic compounds Geologic factors that seem most significant in controlling the distribution and concentration of uranium in Dakota lignites are as follows: 1) stratlgraphic proximity of the lignite to the base of the White River formation; 2) permeability of the rocks overlying the lignites; 3) adsorptive properties and porosities of the lignitic constituents; 4) present and past position of the groundwater table; and 5) the amount of uranium in the original White River and Arikaree sediments. Individual maps showing the extent, thickness, and variations in mineral content of the important deposits in the Table Mountain, Cave Hills, Slim Buttes, Lodgepole, and Medicine Pole areas are included. Conditions controlling the concentration of uranium are described and their application as guides to finding additional reserves by the presently held concepts are explained and illustrated.

Montana;North Dakota;South Dakota;Wyoming↗

Reestablishing a foundational species: limitations on post-wildfire sagebrush seedling establishment

Improving post-wildfire restoration of foundational plant species is crucial for conserving imperiled ecosystems. We sought to better understand the initial establishment of sagebrush ( Artemisia sp.), a foundational shrubland species over a vast area of western North America, in the first 1–2 years post-wildfire, a critical time period for population recovery. Field data from 460 sagebrush populations sampled across the Great Basin revealed several patterns. Sagebrush seedlings were uncommon in the first 1–2 years after fire, with none detected in 69% of plots, largely because most fires occurred in areas of low resistance to invasive species and resilience to disturbance (hereafter, R&R). Post-fire aerial seeding of sagebrush dramatically increased seedling occupancy, especially in low R&R areas, which exhibited a 3.4-fold increase in occupancy over similar unseeded locations. However, occupancy models and repeat surveys suggested exceptionally high mortality, as occupancy rates declined by as much as 50% between the first and second years after fire. We found the prevalence of “fertile island” microsites (patches beneath fire-consumed sagebrush) to be the best predictor of seedling occupancy, followed by aerial seeding status, native perennial grass cover, and years since fire. In populations where no sagebrush seeding occurred, seedlings were most likely to occur in locations with a combination of high fertile island microsite cover and close proximity to a remnant sagebrush plant. These important attributes were only present in 13% of post-fire locations, making them rare across the Great Basin. Finally, in the absence of fertile islands and remnant plants, seedling establishment was not observed in any unseeded areas, and rarely in seeded locations. Thus, local extirpation of sagebrush could have important, long-term implications for sagebrush reestablishment following future fires if there are no mature individuals to leave behind fertile islands or serve as remnant individuals. These findings highlight the importance of landscape legacy effects and could help guide where and how big sagebrush restoration is conducted in the future.

California, Idaho, Nevada, Oregon, Utah↗

Identification of metrics to monitor salt marsh integrity on National Wildlife Refuges in relation to conservation and management objectives

Executive Summary Most salt marshes in the US have been degraded by human activities, and threats from physical alterations, surrounding land-use, species invasions, and global climate change persist. Salt marshes are unique and highly productive ecosystems with high intrinsic value to wildlife, and many National Wildlife Refuges (NWRs) have been established in coastal areas to protect large tracts of salt marsh and wetland-dependent species. Various management practices are employed routinely on coastal NWRs to restore and enhance marsh integrity and ensure ecosystem sustainability. Prioritizing NWR salt marshes for application of management actions and choosing among multiple management options requires scientifically-based methods for assessing marsh condition. Monitoring is integral to structured decision-making (SDM), a formal process for decomposing a decision into its essential elements. Within a natural resource context, SDM involves identifying management objectives, alternative management actions, and expected management outcomes. The core of SDM is a set of criteria for measuring system performance and evaluating management responses. Therefore, use of SDM to frame natural resource decisions leads to logical selection of monitoring attributes that are linked explicitly to management needs. We used SDM to guide selection of variables for monitoring the ecological integrity of salt marshes within the National Wildlife Refuge System (NWRS). Our objectives were to identify indicators of salt marsh integrity that are effective across large geographic regions, responsive to a wide range of threats, and feasible to implement within funding and staffing constraints of the NWRS. In April, 2008, we engaged interdisciplinary experts in a week-long rapid prototyping SDM workshop to define the essential elements of salt marsh management decisions on refuges throughout the northeastern, southwestern, and northwestern US, corresponding to respective Regions 5, 2, and 1 of the US Fish and Wildlife Service (FWS). Through this process we identified measurable attributes for monitoring salt marsh ecosystems that are integrated into conservation practice and target management objectives. The following salt marsh attributes were identified through the SDM process either for describing state condition to determine management needs or for evaluating the achievement of management objectives: historical condition and geomorphic setting; ditch density; surrounding land use; ratio of open water area to vegetation area; rate of pesticide application; environmental contaminant concentration; change in marsh surface elevation relative to sea level rise; tidal range and groundwater level; surface topography; salinity; and species composition and abundance of vegetation, invasive species, invertebrates, nekton, and breeding and wintering birds. The identified attributes were too broadly defined to serve as operational monitoring variables. Therefore, we tested specific metrics for quantifying most of these attributes in summers of 2008 and 2009. The first four attributes in the above list can be characterized by office-based analysis of existing GIS data layers. The remaining attributes require field-based methods for assessment. We were forced to exclude a small number of attributes from field tests due to inconsistent data (pesticide application rate, environmental contaminant concentrations) or requirements that exceeded the scope of this project (change in marsh surface elevation; surface topography; benthic invertebrates; wintering birds). We evaluated potential metrics for evaluating all remaining field attributes. In partnership with NWRS biologists, we tested rapid versus intensive metrics for monitoring field attributes (tidal range and groundwater level; marsh surface elevation; salinity; and species composition and abundance of vegetation, invasive species, nekton, and breeding birds) at coastal refuges throughout FWS Region 5. Seven refuges participated in metric testing in 2008: Rachel Carson (ME), Parker River (MA), Wertheim (NY), E. B. Forsythe (NJ), Bombay Hook (DE), Prime Hook (DE), and Eastern Shore of Virginia Complex (VA). These seven and two additional refuges participated in metric testing in 2009: Rhode Island Complex (RI) and Stewart B. McKinney (CT). We based all field metrics on existing protocols for salt marsh assessment. Sampling locations were determined randomly within delineated marsh study units (MSUs) at each refuge. Detailed field methods are provided in appendices to this report. Measurements for individual metrics were averaged across samples within MSUs during each year of sampling. Each year, correlation or regression analysis was conducted on average measurements across MSUs within each attribute set to identify redundant metrics. Statistical redundancy between a pair of metrics within an attribute set (i.e., correlation or regression slopes with p-values < 0.05) was considered justification for eliminating one of the pair from the regional set of monitoring metrics. Decisions regarding metric elimination versus retention were based on feasibility of monitoring, considering such factors as sampling time, resources required, and potential for regional standardization in implementation. The result of these tests is a reduced suite of monitoring metrics that targets NWRS management decisions and is practicable for implementing on a regional scale. Based on these tests, we recommend the following list of metrics for monitoring integrity of NWRS salt marshes (marsh attribute category is in parentheses): (historical condition and geomorphic setting) position of marsh in the landscape, marsh shape, degree of fill and/or fragmentation, degree of tidal flushing, amount of aquatic edge; (ditch density) ranking of ditch density from none to severe; (surrounding land use) relative proportion of agricultural land in a 150-m buffer around the marsh, relative proportion of natural land in a 150-m buffer around the marsh, relative proportion of natural land in a 1-km buffer around the marsh; (ratio of open water area to vegetation area) ratio of open water to emergent herbaceous wetlands within the marsh; (marsh surface elevation) elevation referenced to NAVD88 in a representative area of the marsh; (tidal range and groundwater level) percent of time the marsh surface is flooded during deployment of a continuous water-level monitor at a representative marsh location, mean depth of surface flooding as measured by a continuous water-level monitor at a representative location; (salinity) salinity measured in surface water; (vegetation community) vegetation species richness using the point-intercept method in 100-m diameter survey plots, percent cover of various marsh community types within 100-m diameter survey plots; (invasive species abundance) percent cover of invasive plant species measured using the point-intercept method in 100-m diameter survey plots; (nekton community) nekton density, nekton species richness, length of Fundulus heteroclitus ; (breeding bird community) abundance of Willets counted per point during standard call-broadcast surveys, summed abundance of tidal marsh obligate species (Clapper Rail, Willet, Saltmarsh Sparrow, Seaside Sparrow) counted per point during standard call-broadcast surveys. Metrics describing the historical condition, geomorphic setting, and broad landscape features can be assessed using existing GIS databases. Our results support use of rapid methods to assess the majority of field metrics; only those used to describe the nekton community must be measured using intensive methods (throw traps or ditch nets, dependant on habitat configuration). Implementation of these metrics for quantitative assessment of NWRS salt marsh integrity in FWS Region 5 requires developing sampling designs for each refuge. Additionally, it is important to determine how the monitoring information will be used within a management context. SDM should be used to complete the analysis of salt marsh management decisions. The next steps would involve 1) prioritizing and weighting the management objectives; 2) predicting responses to individual management actions in terms of objectives and metrics; 3) using multiattribute utility theory to convert all measurable attributes to a common utility scale; 4) determining the total management benefit of each action by summing utilities across objectives; and 5) maximizing the total management benefits within cost constraints for each refuge. This process would allow the optimum management decisions for NWRS salt marshes to be selected and implemented based directly on monitoring data and current understanding of marsh responses to management actions. Monitoring the outcome of management actions would then allow new monitoring data to be incorporated into subsequent decisions.

Connecticut, Delaware, Maine, Massachusetts, New J↗

Guidance for parameterizing post-fire hydrologic models with in situ infiltration measurements

Wildfire can alter soil-hydraulic properties, often resulting in an increased prevalence of infiltration-excess overland flow and greater potential for debris-flow hazards. Mini disk tension infiltrometers (MDIs) can be used to estimate soil hydraulic properties, such as field-saturated hydraulic conductivity ( K fs ) and wetting front potential ( H f ), and their spatial variability following wildfire. However, the small (point-scale) footprint of MDI measurements makes it challenging to use these data to parameterize hydrologic models at the hillslope and watershed scales where hydrologic hazards, such as debris flows, initiate. Here, we designed numerical experiments to estimate spatially constant or watershed-scale effective hydrologic parameters (EHPs) that approximate the response of spatially variable hydrologic parameters with distributions derived from MDI measurements at five sites in the southwestern United States. We found that it is possible to define EHPs for both K fs and H f based on the MDI measurements that lead to reasonable approximations of run-off hydrographs at the outlets of small watersheds (<1 km 2 ). We found that watershed EHPs are functions of rainfall characteristics, although they are most sensitive to rainfall intensity and relatively less sensitive to the temporal distribution of rainfall. EHPs are lower than the arithmetic mean of the MDI measurements and are better approximated by the median or geometric mean of the MDI measurements, particularly for storms with recurrence intervals of approximately 1 year or less that commonly initiate post-fire debris flows. This work demonstrated that using the proposed upscaling method to estimate watershed-scale EHPs, as opposed to approximating EHPs based on the arithmetic mean of the MDI measurements, improved the ability of a hydrologic model to identify storms that are likely to produce debris flows. Results improved our ability to link point-scale MDI measurements and watershed-scale EHPs in post-fire settings and helped guide our ability to use MDI data to parameterize post-fire hydrologic models.

Arizona, New Mexico↗

Lead isotope compositions as guides to early gold mineralization: The North Amethyst vein system, Creede district, Colorado

The North Amethyst vein system, which is hosted by approximately 27 Ma Carpenter Ridge Tuff and approximately 26 Ma Nelson Mountain Tuff, has two mineral associations separated by brecciation and sedimentation in the veins. The early association consists of quartz, rhodonite, hematite, magnetite, electrum (Au (sub 0.3-0.5) Ag (sub 0.7-0.5)) , and Mn carbonate, Au-Ag sulfide, Ag sulfosalt, and base metal sulfide minerals. The later mineral association cuts the Mn- and Au-bearing assemblages and consists of quartz, calcite, sericite, chlorite, hematite, adularia, fluorite, base metal sulfides, and Ag-bearing tetrahedrite.Our detailed studies show that the Pb isotope compositions of paragenetically early galenas associated with Au-rich mineralization in the North Amethyst vein system are relatively unradiogenic ( 206 Pb/ 204 Pb: 18.826-18.881, 207 Pb/ 204 Pb: 15.588-15.602, and 208 Pb/ 204 Pb: 37.790-37.926) compared to Pb isotope compositions of galenas formed later at about 25 Ma during Ag and base metal mineralization ( 206 Pb/ 204 Pb: 19.041-19.115, 207 Pb/ 204 Pb: 15.627- 15.672, and 208 Pb/ 204 Pb: 37.829-38.057). New Pb isotope data for the central and southern parts of the Creede district, which are located 5 to 7 km south of the North Amethyst area, agree with the results of a regional study by Doe et al. (1979) that included five samples from the main part of the Creede district. Galenas and adularia from the central and southern Creede district and galenas from the Bondholder district, 7 km north of the North Amethyst area, are similar to the Pb isotope compositions of galenas formed later in the North Amethyst area during Ag and base metal mineralization. Galenas from the Alpha-Corsair vein, which was mined for Ag prior to 1910, are isotopically similar to galenas associated with North Amethyst Au-stage mineralization. This isotopic similarity suggests that unexplored segments of the Alpha-Corsair structure may have the same mineralogy as the North Amethyst Au stage; thus, the Alpha-Corsair structure has the potential for high gold contents.Pb isotope compositions from the late stage of the North Amethyst vein system and from the Bondholder and central and southern Creede mining districts are more radiogenic than the host volcanic rocks of the central cluster of the San Juan volcanic field. Our Pb isotope results indicate that early Au mineralization of the North Amethyst area may represent the product of an older and relatively local hydrothermal system distinct from that of the younger base metal and Ag mineralization found throughout the region. Fluids that deposited Au minerals may have derived their Pb isotope composition by a greater degree of interaction with shallow, relatively less radiogenic volcanic wall rocks. The younger, base metal and Ag-rich mineralization that overprints the Au mineralization in the North Amethyst area clearly has a more radiogenic isotopic signature, which implies that the later mineralization derived a greater component of its Pb from Proterozoic source rocks, or sediments derived from them.Paragenetically early sulfide-rich vein assemblages have the least radiogenic galenas and generally also have the highest Au contents. Thus, identification of paragenetically early vein assemblages with relatively unradiogenic Pb isotope compositions similar to those of the North Amethyst area provides an additional exploration tool for Au in the central San Juan Mountains area.

Colorado↗

Methodological considerations of terrestrial laser scanning for vegetation monitoring in the sagebrush steppe

Terrestrial laser scanning (TLS) provides fast collection of high-definition structural information, making it a valuable field instrument to many monitoring applications. A weakness of TLS collections, especially in vegetation, is the occurrence of unsampled regions in point clouds where the sensor’s line-of-sight is blocked by intervening material. This problem, referred to as occlusion, may be mitigated by scanning target areas from several positions, increasing the chance that any given area will fall within the scanner’s line-of-sight from at least one position. Because TLS collections are often employed in remote regions where the scope of sampling is limited by logistical factors such as time and battery power, it is important to design field protocols which maximize efficiency and support increased quantity and quality of the data collected. This study informs researchers and practitioners seeking to optimize TLS sampling methods for vegetation monitoring in dryland ecosystems through three analyses. First, we quantify the 2D extent of occluded regions based on the range from single scan positions. Second, we measure the efficacy of additional scan positions on the reduction of 2D occluded regions (area) using progressive configurations of scan positions in 1 ha plots. Third, we test the reproducibility of 3D sampling yielded by a 5-scan/ha sampling methodology using redundant sets of scans. Analyses were performed using measurements at analysis scales of 5 to 50 cm across the 1-ha plots, and we considered plots in grass and shrub-dominated communities separately. In grass-dominated plots, a center-scan configuration and 5 cm pixel size sampled at least 90% of the area up to 18 m away from the scanner. In shrub-dominated plots, sampling at least 90% of the area was only achieved within a distance of 12 m. We found that 3 and 5 scans/ha are needed to sample at least ~ 70% of the total area (1 ha) in the grass and shrub-dominated plots, respectively, using 5 cm pixels to measure sampling presence-absence. The reproducibility of 3D sampling provided by a 5 position scan layout across 1-ha plots was 50% (shrub) and 70% (grass) using a 5-cm voxel size, whereas at the 50-cm voxel scale, reproducibility of sampling was nearly 100% for all plot types. Future studies applying TLS in similar dryland environments for vegetation monitoring may use our results as a guide to efficiently achieve sampling coverage and reproducibility in datasets.

Idaho↗

Geology, coastal geomorphology, and soils of eastern Santa Cruz Island (Limuw), Channel Islands National Park, California, USA

This one-day field trip explores northeastern Santa Cruz Island ( Limuw , in native Chumash), a part of Channel Islands National Park, USA. The geomorphology of eastern Santa Cruz Island has been controlled largely by active tectonics and sea-level fluctuations. The bedrock is Miocene volcanic rock overlain by Miocene shale and siltstone. The island has experienced Quaternary uplift, perhaps due to movement along an offshore thrust fault. Smaller faults are exposed in sea cliffs and have displaced Miocene rocks. Superimposed upon island uplift, there have been Quaternary sea-level fluctuations from interglacial-glacial climate changes. Interglacial high-sea stands are recorded as marine terraces. The last major interglacial period, ~120,000 years ago, left only small remnants of marine terraces. Most evidence of this high-sea stand was eroded away in the Holocene. However, a prominent marine terrace is preserved at 75–120 m above sea level. Some fossil mollusks from the deposits of this terrace, probably reworked, have given ages as old as Pliocene, but most yield ages of 2.6–2.0 Ma. The age and elevation of this terrace indicate a very low rate of tectonic uplift, similar to nearby Anacapa Island. A low uplift rate explains the absence or scarcity of younger terraces, including that of the last interglacial period. Low stands of sea (glacial periods) exposed the insular shelf, rich in carbonate skeletal sand. During glacial periods, these sands were entrained by the wind, deposited as dunes on marine terraces, and cemented into eolianite. Clay-rich Vertisols with silt mantles have developed on eolianites and terraces of the island, partly from in situ weathering, but also from inputs of Mojave Desert dust during Santa Ana wind events. This guide includes stops at Scorpion Anchorage, Cavern Point, and Potato Harbor. It provides insights into the bedrock, coastal geomorphology, fossiliferous marine terraces, eolianite, Vertisols, and the three formations on eastern Santa Cruz Island: the Santa Cruz Island Volcanics, the Monterey Formation, and the Potato Harbor Formation.

California↗

An international code comparison study on coupled thermal, hydrologic and geomechanical processes of natural gas hydrate-bearing sediments

Geologic reservoirs containing gas hydrate occur beneath permafrost environments and within marine continental slope sediments, representing a potentially vast natural gas source. Numerical simulators provide scientists and engineers with tools for understanding how production efficiency depends on the numerous, interdependent (coupled) processes associated with potential production strategies for these gas hydrate reservoirs. Confidence in the modeling and forecasting abilities of these gas hydrate reservoir simulators (GHRSs) grows with successful comparisons against laboratory and field test results, but such results are rare, particularly in natural settings. The hydrate community recognized another approach to building confidence in the GHRS: comparing simulation results between independently developed and executed computer codes on structured problems specifically tailored to the interdependent processes relevant for gas hydrate-bearing systems. The United States Department of Energy, National Energy Technology Laboratory (DOE/NETL), sponsored the first international gas hydrate code comparison study, IGHCCS1, in the early 2000s. IGHCCS1 focused on coupled thermal and hydrologic processes associated with producing gas hydrates from geologic reservoirs via depressurization and thermal stimulation. Subsequently, GHRSs have advanced to model more complex production technologies and incorporate geomechanical processes into the existing framework of coupled thermal and hydrologic modeling. This paper contributes to the validation of these recent GHRS developments by providing results from a second GHRS code comparison study, IGHCCS2, also sponsored by DOE/NETL. IGHCCS2 includes participants from an international collection of universities, research institutes, industry, national laboratories, and national geologic surveys. Study participants developed a series of five benchmark problems principally involving gas hydrate processes with geomechanical components. The five problems range from simple geometries with analytical solutions to a representation of the world’s first offshore production test of methane hydrates, which was conducted with the depressurization method off the coast of Japan. To identify strengths and limitations in the various GHRSs, study participants submitted solutions for the benchmark problems and discussed differing results via teleconferences. The GHRSs evolved over the course of IGHCCS2 as researchers modified their simulators to reflect new insights, lessons learned, and suggested performance enhancements. The five benchmark problems, final sample solutions, and lessons learned that are presented here document the study outcomes and serve as a reference guide for developing and testing gas hydrate reservoir simulators.

Journal of Marine and Petroleum Geology↗

Geophysical investigations at Hidden Dam, Raymond, California — Flow simulations

Numerical flow modeling and analysis of observation-well data at Hidden Dam are carried out to supplement recent geophysical field investigations at the site (Minsley and others, 2010). This work also is complementary to earlier seepage-related studies at Hidden Dam documented by Cedergren (1980a, b). Known seepage areas on the northwest right abutment area of the downstream side of the dam was documented by Cedergren (1980a, b). Subsequent to the 1980 seepage study, a drainage blanket with a sub-drain system was installed to mitigate downstream seepage. Flow net analysis provided by Cedergren (1980a, b) suggests that the primary seepage mechanism involves flow through the dam foundation due to normal reservoir pool elevations, which results in upflow that intersects the ground surface in several areas on the downstream side of the dam. In addition to the reservoir pool elevations and downstream surface topography, flow is also controlled by the existing foundation geology as well as the presence or absence of a horizontal drain in the downstream portion of the dam. The current modeling study is aimed at quantifying how variability in dam and foundation hydrologic properties influences seepage as a function of reservoir stage. Flow modeling is implemented using the COMSOL Multiphysics software package, which solves the partially saturated flow equations in a two-dimensional (2D) cross-section of Hidden Dam that also incorporates true downstream topography. Use of the COMSOL software package provides a more quantitative approach than the flow net analysis by Cedergren (1980a, b), and allows for rapid evaluation of the influence of various parameters such as reservoir level, dam structure and geometry, and hydrogeologic properties of the dam and foundation materials. Historical observation-well data are used to help validate the flow simulations by comparing observed and predicted water levels for a range of reservoir elevations. The flow models are guided by, and discussed in the context of, the geophysical work (Minsley and others, 2010) where appropriate.

California↗