Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Detailed procedure”

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 289 records · Page 16Linked to original sources
USGS2026

Wavelet Inversion for SliP (WISP): Open-source earthquake slip modeling software

Models of the spatiotemporal evolution of earthquake slip, termed finite-fault models, are a critical component of rapid earthquake and tsunami response, earthquake forecasting, seismic ground-motion estimates, and studies of earthquake kinematics. Here, we detail a newly released finite-fault modeling software, Wavelet Inversion for SliP (WISP), in use at the U.S. Geological Survey’s National Earthquake Information Center (NEIC) and available to the public. WISP version 1.1.0 allows inversion of teleseismic body and surface waves, as well as local strong-motion, static and dynamic Global Navigation Satellite System, and satellite imagery (e.g., Interferometric Synthetic Aperture Radar) observations on single or multiple planar fault segments. The software is used in NEIC rapid response of earthquakes M w ≥ 7, generally resulting in a published model within the first few hours after the event origin time. The rupture location and dimensions are then used as inputs to downstream products to estimate earthquake shaking, predict loss, and model the likelihood of secondary hazards, namely landslides and liquefaction. WISP is also used in research studies to evaluate the characteristics of complex ruptures including multifault ruptures and earthquake doublets, among others. The WISP version 1.1.0 software release is composed of Python-wrapped FORTRAN code to accomplish the inversion procedure. A simple command line interface facilitates ease of use even for those with only a cursory knowledge of Python scripting. WISP version 1.1.0 includes a Jupyter Notebook tutorial demonstrating use of the software for modeling the 2015 M w 8.3 Illapel, Chile, earthquake. In parallel with the tutorial, we demonstrate the typical usage of the WISP software using the M w 8.3 Illapel earthquake example here.

Seismological Research Letters↗
USGS2023

Mapping abandoned uranium mine features using Worldview-3 imagery in portions of Karnes, Atascosa and Live Oak Counties, Texas

Worldview-3 (WV3) 16-band multispectral data were used to map exposed bedrock and mine waste piles associated with legacy open-pit mining of sandstone-hosted roll-front uranium deposits along the South Texas Coastal Plain. We used the “spectral hourglass” approach to extract spectral endmembers representative of these features from the image. This approach first requires calibrating the imagery to reflectance, then masking for vegetation, followed by spatial and spectral data reduction using a principal component analysis-based procedure that reduces noise and identifies homogeneous targets which are “pure” enough to be considered spectral endmembers. In this case, we used a single WV3 image which covered an ~11.5 km by ~19.5 km area of Karnes, Atascosa and Live Oak Counties, underlain by mined rocks from the Jackson Group and Catahoula Formation. Up to 58 spectral endmembers were identified using a further multi-dimensional class segregation method and were used as inputs for spectral angle mapper (SAM) classification. SAM classification resulted in the identification of at least 117 mine- and mine waste-related features, most of which were previously unknown. Class similarity was further evaluated, and the dominant minerals in each class were identified by comparison to spectral libraries and measured samples of actual Jackson Group uranium host rocks. Redundant classes were eliminated, and SAM was run a second time using a reduced set of 23 endmembers, which were found to map these same features as effectively as using the full 58 set of endmembers, but with significantly reduced noise and spectral outliers. Our classification results were validated by evaluating detailed scale mapping of three known mine sites (Esse-Spoonamore, Wright-McCrady and Garbysch-Thane) with published ground truth information about the vegetation cover, extent of erosion and exposure of waste pile materials and/or geologic information about host lithology and mineralization. Despite successful demonstration of the utility of WV3 data for inventorying mine features, additional landscape features such as bare agricultural fields and oil and gas drill pads were also identified. The elimination of such features will require combining the spectral classification maps presented in this study with high-quality topographic data. Also, the spectral endmembers identified during the course of this study could be useful for larger-scale mapping efforts using additional well-calibrated WV3 images beyond the coverage of our initial study area.

Texas↗
USGS2002

Preliminary assessment of using tree-tissue analysis and passive-diffusion samplers to evaluate trichloroethene contamination of ground water at Site SS-34N, McChord Air Force Base, Washington, 2001

Two low-cost innovative sampling procedures for characterizing trichloroethene (TCE) contamination in ground water were evaluated for use at McChord Air Force Base (AFB) by the U.S. Geological Survey, in cooperation with the U.S. Air Force McChord Air Force Base Installation Restoration Program, in 2001. Previous attempts to characterize the source of ground-water contamination in the heterogeneous glacial outwash aquifer at McChord site SS-34N using soil-gas surveys, direct-push exploration, and more than a dozen ground-water monitoring wells have had limited success. The procedures assessed in this study involved analysis of tree-tissue samples to map underlying ground-water contamination and deploying passive-diffusion samplers to measure TCE concentrations in existing monitoring wells. These procedures have been used successfully at other U.S. Department of Defense sites and have resulted in cost avoidance and accelerated site characterization. Despite the presence of TCE in ground water at site SS-34N, TCE was not detected in any of the 20 trees sampled at the site during either early spring or late summer sampling. The reason the tree tissue procedure was not successful at the McChord AFB site SS-34N may have been due to an inability of tree roots to extract moisture from a water table 30 feet below the land surface, or that concentrations of TCE in ground water were not large enough to be detectable in the tree tissue at the sampling point. Passive-diffusion samplers were placed near the top, middle, and bottom of screened intervals in three monitoring wells and TCE was observed in all samplers. Concentrations of TCE from the passive-diffusion samplers were generally similar to concentrations found in samples collected in the same wells using conventional pumping methods. In contrast to conventional pumping methods, the collection of ground-water samples using the passive-diffusion samples did not generate waste purge water that would require hazardous-waste disposal. In addition, the results from the passive-diffusion samples may show that TCE concentrations are stratified across some screened intervals. The overall results of the limited test of passive-diffusion samplers at site SS-34N were similar to more detailed tests conducted at other contaminated sites across the country and indicate that further evaluation of the use of passive-diffusion samplers at McChord site SS-34N is warranted.

Water-Resources Investigations Report↗
USGS2021

Oyster model inventory: Identifying critical data and modeling approaches to support restoration of oyster reefs in coastal U.S. Gulf of Mexico waters

Executive Summary Along the coast of the U.S. Gulf of Mexico, the eastern oyster ( Crassostrea virginica ) plays important ecological and economic roles. Commercial landings from this region account for more than 50 percent of all U.S. landings; these oyster reefs also provide varied ecosystem services, including nursery habitat for many fish and macroinvertebrate species, shoreline protection, and water-quality maintenance. Declining trends in both total oyster production and functional reef area across this region have spurred investment in restoration of oyster resources, with specific calls for restoration projects to develop a network of reefs and identify broodstock and sanctuary reef restoration sites. Decision making related to restoration and establishment of a network of oyster reefs in the Gulf of Mexico requires information on both the environment and the effects of the environment on the oyster life cycle (including larval movement, survival, oyster recruitment, reproduction, growth, and mortality). Here, we examined the current state of data and model development in this region with the goal of providing an overview of oyster modeling approaches and an inventory of available data and existing oyster models. This report is meant to provide an overview to managers for understanding existing efforts and identify a path forward to most efficiently inform oyster resource management and restoration planning in moving from a single reef management approach to a reef network management approach. Numerous models related to some aspect of the oyster life cycle have been built, calibrated, and validated for various Gulf of Mexico estuaries over the last few decades (over 30 models identified). These models, which could inform site restoration, can be classified into four approaches: (1) oyster Habitat Suitability Index (HSI) models; (2) larval transport models; (3) on-reef oyster models that may include oyster growth, mortality and reproduction, and substrate persistence; and (4) coupled larval transport on-reef metapopulation models that simulate the entire oyster life cycle. The data requirements, model complexity and assumptions, and transferability vary by approach. Specifically, some approaches may offer greater accessibility, flexibility, and transferability spatially or temporally, with minimal data input, but only provide broad information to support site selection. In contrast, other approaches may require significant site-specific data for their construction and validation but may provide more accurate and location-specific data to support site selection for broodstock reefs. Regardless of modeling approach used, data on environmental drivers, such as salinity, water temperature, or water flow impacting oyster metabolism and movement, are required at appropriate spatial and temporal scales. While numerous data collection platforms, environmental models, and research products exist within Gulf of Mexico estuaries to provide important environmental data to use as drivers in the oyster models, significant variability in temporal and spatial coverage of the data, and variation in the availability of future condition models, exists across estuaries. This variation influences the spatial and temporal scales at which oyster models may be developed and impacts the calibration and validation of the oyster models within a given estuary, affecting its potential ability to address specific management or restoration questions. While multiple modeling approaches exist for informing site selection of broodstock or sanctuary oyster reefs, the development, calibration, and validation of a single modeling platform presents the most efficient, transferable, and useful tool for managers across the Gulf of Mexico. The development of a single modeling platform would involve using standardized input variables, governing equations, and assumptions for the modeled oyster processes and outputs, and for standardized calibration and validation procedures that could be applied within each estuary. The differences among estuary applications would require substituting only estuary-specific environmental data, and calibrating and validating the modeling approach with local oyster data. Two modeling approaches likely to be useful include (1) development of a general geospatial HSI modeling framework that could be applied consistently across estuaries and (2) a mechanistic coupled larval transport on-reef metapopulation model requiring only estuarine specific calibration and hydrodynamic models. Both approaches benefit from existing work across multiple Gulf of Mexico estuaries and could provide valuable support for oyster restoration, but may differ in their ability to address specific questions related to oyster restoration. HSI models specifically guide restoration practitioners in determining suitable habitat based on available data. The HSI approach, while currently more widely used and accessible, requires more development of larval suitability and larval input and output components in order to inform reef connectivity. A metapopulation approach considering the full oyster life cycle that simulates both on-reef oyster growth, mortality, reproduction, substrate persistence, and larval transport (ideally with larval growth and mortality) would provide the greatest detail and level of understanding but requires significant up-front investment. The larval oyster model and on-reef oyster model are usually developed independently for systems, although the two approaches can be coupled to represent the entire oyster life cycle in order to characterize and assess a reef metapopulation. This approach may be less accessible and much more data-intensive, however, and it requires some expertise to run and apply to inform oyster resource management. Ultimately, the development of single modeling platforms for each of these approaches would provide flexible tools applicable across all Gulf of Mexico oyster supporting estuaries. By using a single platform for model development, testing, calibrating and validating, and evaluation of modeled future scenarios, oyster restoration scientists and managers would not only be able to examine different scenario outcomes within a single estuary, but could also have comparable modeled results to evaluate potential outcomes, across estuaries and regions, that are not confounded by varying modeled data inputs, governing equations, assumptions, or user judgement.

Alabama, Florida, Louisiana, Mississippi, Texas↗
USGS2016

Automatic delineation of seacliff limits using lidar-derived high-resolution DEMs in southern California

Seacliff erosion is a serious hazard with implications for coastal management and is often estimated using successive hand-digitized cliff tops or bases (toe) to assess cliff retreat. Even if efforts are made to standardize manual digitizing and eliminate subjectivity, the delineation of cliffs is time-consuming and depends on the analyst's interpretation. An automatic procedure is proposed to extract cliff edges from high-resolution lidar-derived bare-earth digital elevation models, generalized coastal shoreline vectors, and approximate measurements of distance between the shoreline and the cliff top. The method generates orthogonal transects and profiles with a minimum spacing equal to the digital elevation model resolution. The method also extracts the xyz coordinates for each profile for the cliff top and toe, as well as second major inflections along the profile. Over 75% of the automated cliff top points and 78% of the toe automated points are within 95% confidence interval of the hand-digitized top and toe lines, and over 79% of the digitized top points and 84% of the digitized toe points are within the 95% confidence interval of the automated top and toe lines along a stretch of coast in Del Mar, California. Outlier errors were caused by either the failure to remove all vegetation from the bare-earth digital elevation model or errors of interpretation. The automatic method was further applied between Point Conception and Los Angeles Harbor, California. This automatic method is repeatable, takes advantage of detailed topographic information within high-resolution digital elevation models, and is more efficient than hand-digitizing.

California↗
USGS1977

Population ecology of the mallard: IV. A review of duck hunting regulations, activity, and success, with special reference to the mallard

This, the fourth in a series of reports on the mallard, ( Anas platyrhynchos ), deals at length with the harvest of mallards by waterfowl hunters. Long-term summaries of duck hunting regulations (1948- 1974), Migratory Bird Hunting Stamp sales (1934-1974), Hunter Questionnaire (1952-1974), Duck Wing Collection (1960-1974), and Hunter Performance (1965-1972) Survey data for the United States are presented and discussed. Similar data from Canada are also summarized. Mallard harvest figures for 1961-1974 are presented by Mallard Harvest Area, of which 100 are defined for the United States and 14 for Canada, as well as by State or Province and flyway. During the 23-year period beginning in 1952, an average of 1.6 million adult and 0.2 million junior waterfowl hunters accumulated almost 12.3 million hunter-days of recreation and a harvest of 11.2 million ducks each year. Hunter reports indicate that mallards made up about 43% (5.5 million annually) of the ducks taken before 1960, when mallard regulations were less restrictive; the Duck Wing Survey indicates that mallards have made up 33% of the harvest (3.6 million annually) since 1960. The age and sex compositions and the chronological distribution of the mallard harvest are examined in detail. Among the patterns noted are peak harvests during the first few days of the season in many States, alternately increasing and decreasing annual age ratios, and sex ratios that suggest differential migration of adult drakes and hunter selectivity for males. It is estimated that almost 19% of the ducks shot down are not retrieved. Relationships between duck hunting regulations and hunter behavior are examined briefly. Hunter compliance with mallard bag limits, hunter selectivity of mallards by sex, and, to a lesser extent, the size of the unretrieved kill are all sensitive to the particular bag limit regulations in effect. Survey data are also examined for relationships between harvest and various hunting regulations: starting time and day of the week for opening day, opening date, season length, split seasons, daily shooting hours, and daily bag limits. Tables are presented relating changes in duck and mallard harvests to season length and bag limit, and examples of what effects changes in other regulations have on harvest are also given. The evaluation of bag limit regulations, one of the most important tools used in managing harvest, is carried a step further with the development of a procedure for calculating expected hunter success under a wide variety of bag limit regulations. This method appears promising for evaluating point-limit as well as fixed-limit regulations. In addition, it may provide useful measurements of the degree of hunter selectivity induced by various types of bag limit regulations, an increasingly important aspect of harvest management. Finally, this study clearly demonstrates that the effects of a particular regulation can differ dramatically from area to area, so it is usually necessary to evaluate each proposal on a State-by-State basis.

Resource Publication↗
USGS2015

Nutrient attenuation in rivers and streams, Puget Sound Basin, Washington

Nutrients such as nitrogen and phosphorus are important for aquatic ecosystem health. Excessive amounts of nutrients, however, can make aquatic ecosystems harmful for biota because enhanced growth and decay cycles of aquatic algae can reduce dissolved oxygen in the water. In Puget Sound marine waters, low dissolved oxygen concentrations are observed in a number of marine nearshore areas, and nutrients have been identified as a major stressor to the local ecosystem. Delivery of nutrients from major rivers in the Puget Sound Basin to the marine environment can be large. Therefore, it is important to identify factors related to how nutrients are retained (attenuated) within streams and rivers in the Puget Sound Basin. Physical, chemical, and biological factors related to nutrient attenuation were identified through a review of related scientific literature. Numerous empirical modeling approaches for estimating nutrient attenuation in streams and rivers also were compiled, and a subset of these models was applied to the Puget Sound Basin. In particular, models based on the physical characteristics of a river reach (RivR-N model) and on the physical and biological features of a river reach (v f model) were used and compared for the 17 major rivers draining to the Puget Sound. Data on the relative amount of instream attenuation (the fraction of nutrient input removed per kilometer of stream reach) showed some common and general themes. Firstly, headwater reaches throughout the Puget Sound Basin tend to be better than the main stems of the major rivers at attenuating nitrate and orthophosphorus (ortho-P). Secondly, rivers are more efficient at attenuating nitrate than ortho-P, probably because of the close relation between phosphorus and suspended sediment, which was not captured fully in the models. Thirdly, when comparing the RivR-N and v f models for nitrate, physical characteristics of the channel may be more effective predictors of relative nitrate attenuation for main stem reaches, whereas biological factors may be more effective predictors in headwater reaches. These results are explained in terms of four primary factors of attenuation: sinuosity, channel slope, specific discharge, and uptake velocity (v f ) of the reach. A simple scoring procedure based on these four factors showed that reaches where attenuation scores were high had higher relative attenuation of nutrients from the RivR-N and v f models. This attenuation "scorecard" can be used to quickly assess the potential for a given reach to attenuate nutrients. Seasonal relative attenuation at three case studies was greater in summer months (July through September) and much lower and almost constant from January through June. An analysis of relative attenuation across a range of nutrient concentrations showed that, at some point, relative instream attenuation is minimized. For nitrate, relative attenuation reached a minimum value greater than 3 milligrams of nitrogen per liter (mg N/L) during low flow and 1 mg N/L during high flow. For orthophosphate, minimum relative attenuation was observed at about 0.1 milligram of phosphorus per liter (mg P/L) for both low- and high-flow conditions. Generally, the temporal dynamics of nutrient attenuation are dependent on the travel time through a given reach, the proportion of flow in contact with the sediment, and the amount of biological activity. Improved understanding of nutrient attenuation in Puget Sound Basin will benefit from the compilation of more detailed data for specific discharge, channel slope, and channel sinuosity in Puget Sound streams and rivers. Additionally, field studies examining upstream-downstream changes in nutrient load and field-based measurements of v f are needed. From a management perspective, preservation and improvement of instream nutrient attenuation should focus on increasing the travel time through a reach and contact time of water sediment (reactive) surfaces and lowering nutrient concentrations (and loads) to avoid saturation of instream attenuation and increase attenuation efficiency. These goals can be reached by maintaining and restoring channel-flood plain connectivity, maintaining and restoring healthy riparian zones along streams, managing point and nonpoint nutrient loads to streams and rivers, and restoring channel features that promote attenuation such as the addition of woody debris and maintaining pool-riffle morphologies. Many of these management approaches are already being undertaken during projects aimed to restore quality salmon habitat. Therefore, there is a dual benefit to these projects that also may lead to enhanced potential for nitrogen and phosphorus attenuation.

Washington↗
USGS2026

Creating 3D point clouds, digital elevation models, and orthomosaics from historical aerial imagery through structure from motion aided photogrammetry

Detailed and accurate historical data are crucial to characterizing landscape change, whether it is caused by natural phenomena or of anthropogenic origin. Aerial imagery has been periodically collected by government and nongovernment agencies, for a broad range of purposes, since the early 20th century. The acquisition of aerial images involves using standardized equipment with large-format, carefully calibrated cameras and lens assemblies. Historical archives of aerial images preserve a record of the landscape and are invaluable sources of information that predate satellite and lidar. Structure from motion (SfM) aided photogrammetry provides an efficient means for creating digital three-dimensional data from the vast archives of historical aerial imagery. This report presents techniques, methods, best practices, and a standard workflow for creating point clouds, digital elevation models, and orthomosaic digital scans of overlapping (stereographic) vertical aerial photographs using photogrammetry software. A list of critical issues relevant to SfM-aided photogrammetry of historical aerial imagery is provided, along with potential contingencies and solutions. The workflow was tested in many geographic study areas and with historical imagery datasets of different scales and sources. Three case studies are included. By integrating known SfM methods, these procedures enable the consistent and optimal extraction of historical quantitative data from archives like the U.S. Geological Survey (USGS) Aerial Photo Single Frame archive in the USGS EarthExplorer data portal.

Techniques and Methods↗
USGS2009

Hydrodynamic modeling of tsunamis from the Currituck landslide

Tsunami generation from the Currituck landslide offshore North Carolina and propagation of waves toward the U.S. coastline are modeled based on recent geotechnical analysis of slide movement. A long and intermediate wave modeling package (COULWAVE) based on the non-linear Boussinesq equations are used to simulate the tsunami. This model includes procedures to incorporate bottom friction, wave breaking, and overland flow during runup. Potential tsunamis generated from the Currituck landslide are analyzed using four approaches: (1) tsunami wave history is calculated from several different scenarios indicated by geotechnical stability and mobility analyses; (2) a sensitivity analysis is conducted to determine the effects of both landslide failure duration during generation and bottom friction along the continental shelf during propagation; (3) wave history is calculated over a regional area to determine the propagation of energy oblique to the slide axis; and (4) a high-resolution 1D model is developed to accurately model wave breaking and the combined influence of nonlinearity and dispersion during nearshore propagation and runup. The primary source parameter that affects tsunami severity for this case study is landslide volume, with failure duration having a secondary influence. Bottom friction during propagation across the continental shelf has a strong influence on the attenuation of the tsunami during propagation. The high-resolution 1D model also indicates that the tsunami undergoes nonlinear fission prior to wave breaking, generating independent, short-period waves. Wave breaking occurs approximately 40-50??km offshore where a tsunami bore is formed that persists during runup. These analyses illustrate the complex nature of landslide tsunamis, necessitating the use of detailed landslide stability/mobility models and higher-order hydrodynamic models to determine their hazard.

Marine Geology↗
USGS1992

Decomposition techniques

Sample decomposition is a fundamental and integral step in the procedure of geochemical analysis. It is often the limiting factor to sample throughput, especially with the recent application of the fast and modern multi-element measurement instrumentation. The complexity of geological materials makes it necessary to choose the sample decomposition technique that is compatible with the specific objective of the analysis. When selecting a decomposition technique, consideration should be given to the chemical and mineralogical characteristics of the sample, elements to be determined, precision and accuracy requirements, sample throughput, technical capability of personnel, and time constraints. This paper addresses these concerns and discusses the attributes and limitations of many techniques of sample decomposition along with examples of their application to geochemical analysis. The chemical properties of reagents as to their function as decomposition agents are also reviewed. The section on acid dissolution techniques addresses the various inorganic acids that are used individually or in combination in both open and closed systems. Fluxes used in sample fusion are discussed. The promising microwave-oven technology and the emerging field of automation are also examined. A section on applications highlights the use of decomposition techniques for the determination of Au, platinum group elements (PGEs), Hg, U, hydride-forming elements, rare earth elements (REEs), and multi-elements in geological materials. Partial dissolution techniques used for geochemical exploration which have been treated in detail elsewhere are not discussed here; nor are fire-assaying for noble metals and decomposition techniques for X-ray fluorescence or nuclear methods be discussed.

Journal of Geochemical Exploration↗
USGS1961

Time of travel of water in the Potomac River, Cumberland to Washington

This report introduces a graphical procedure for estimating the time required for water to travel down the Potomac River in the reach extending from Cumberland, Md., to Washington, D.C. The time of travel varies with the flow of the river; so the stage of the river at the lower end of the reach--the gaging station on the Potomac River near Washington, D.C.--is used as an index of flow. To develop the procedure, the reach between Cumberland and Washington was divided into five subreaches, delineated by six gaging stations. The average of the mean velocities of the river at adjacent gaging stations was used as the mean velocity in .the intervening subreach, and a unit mass of water was assumed to travel at a rate equal to the mean velocity of the river. A statistical analysis of possible variations in travel time between Cumberland and Washington indicated that the shortest travel time corresponding to a given stage near Washington would be about 80 percent of the most probable travel time. The report includes a flow-duration curve and a flow-frequency chart for use in estimating discharge at the gaging station near Washington and subsequently the travel time of Potomac River water without knowledge of stage. The flow-duration curve shows the percentage of time during which specified discharges were equaled or exceeded in the past, and it can be used to predict future flow in connection with long-range planning. The flow-frequency chart shows the time distribution of flow by months and can be used to make a more nearly accurate estimate of discharge in any given month than could be made from the flow-duration curve. The method used to develop the time-of-travel charts is described in sufficient detail to make it usable as a guide for similar studies on other rivers, where the velocity of flow is relatively unaffected by dams and pools in the reach being studied.

Circular↗
USGS2005

Analysis and mapping of post-fire hydrologic hazards for the 2002 Hayman, Coal Seam, and Missionary Ridge wildfires, Colorado

Wildfires caused extreme changes in the hydrologic, hydraulic, and geomorphologic characteristics of many Colorado drainage basins in the summer of 2002. Detailed assessments were made of the short-term effects of three wildfires on burned and adjacent unburned parts of drainage basins. These were the Hayman, Coal Seam, and Missionary Ridge wildfires. Longer term runoff characteristics that reflect post-fire drainage basin recovery expected to develop over a period of several years also were analyzed for two affected stream reaches: the South Platte River between Deckers and Trumbull, and Mitchell Creek in Glenwood Springs. The 10-, 50-, 100-, and 500-year flood-plain boundaries and water-surface profiles were computed in a detailed hydraulic study of the Deckers-to-Trumbull reach. The Hayman wildfire burned approximately 138,000 acres (216 square miles) in granitic terrain near Denver, and the predominant potential hazard in this area is flooding by sediment-laden water along the large tributaries to and the main stem of the South Platte River. The Coal Seam wildfire burned approximately 12,200 acres (19.1 square miles) near Glenwood Springs, and the Missionary Ridge wildfire burned approximately 70,500 acres (110 square miles) near Durango, both in areas underlain by marine shales where the predominant potential hazard is debris-flow inundation of low-lying areas. Hydrographs and peak discharges for pre-burn and post-burn scenarios were computed for each drainage basin and tributary subbasin by using rainfall-runoff models because streamflow data for most tributary subbasins were not available. An objective rainfall-runoff model calibration method based on nonlinear regression and referred to as the ?objective calibration method? was developed and applied to rainfall-runoff models for three burned areas. The HEC-1 rainfall-runoff model was used to simulate the pre-burn rainfall-runoff processes in response to the 100-year storm, and HEC-HMS was used for runoff hydrograph generation. Post-burn rainfall-runoff parameters were determined by adjusting the runoff-curve numbers on the basis of a weighting procedure derived from the U.S. Soil Conservation Service (now the National Resources Conservation Service) equation for precipitation excess and the effect of burn severity. This weighting procedure was determined to be more appropriate than simple area weighting because of the potentially marked effect of even small burned areas on the runoff hydrograph in individual drainage basins. Computed water-peak discharges from HEC-HMS models were increased volumetrically to account for increased sediment concentrations that are expected as a result of accelerated erosion after burning. Peak discharge estimates for potential floods in the South Platte River were increased by a factor that assumed a volumetric sediment concentration (Cv) of 20 percent. Flood hydrographs for the South Platte River and Mitchell Creek were routed down main-stem channels using watershed-routing algorithms included in the HEC-HMS rainfall-runoff model. In areas subject to debris flows in the Coal Seam and Missionary Ridge burned areas, debris-flow discharges were simulated by 100-year rainfall events, and the inflow hydrographs at tributary mouths were simulated by using the objective calibration method. Sediment concentrations (Cv) used in debris-flow simulations were varied through the event, and were initial Cv 20 percent, mean Cv approximately 31 percent, maximum Cv 48 percent, Cv 43 percent at the time of the water hydrograph peak, and Cv 20 percent for the duration of the event. The FLO-2D flood- and debris-flow routing model was used to delineate the area of unconfined debris-flow inundation on selected alluvial fan and valley floor areas. A method was developed to objectively determine the post-fire recovery period for the Hayman and Coal Seam burned areas using runoff-curve numbers (RCN) for all drainage basins for a 50-year period. A

Scientific Investigations Report↗
USGS2023

Interlaboratory comparison of testing hydraulic, elastic, and failure properties in compression: Lessons learned

Many geoscientific problems require us to exploit synergies of experimental and numerical approaches, which in turn lead to questions regarding the significance of experimental details for validation of numerical codes. We report results of an interlaboratory comparison regarding experimental determination of mechanical and hydraulic properties of samples from five rock types, three sandstone varieties with porosities ranging from 5% to 20%, a marble, and a granite. The objective of this study was to build confidence in the participating laboratories’ testing approaches and to establish tractable standards for several physical properties of rocks. We addressed the issue of sample-to-sample variability by investigating the variability of basic physical properties of samples of a particular rock type and by performing repeat tests. Compressive strength of the different rock types spans an order of magnitude and shows close agreement between the laboratories. However, differences among stress–strain relations indicate that the external measurement of axial displacement and the determination of system stiffness require special attention, apparently more so than the external load measurement. Furthermore, post-failure behavior seems to exhibit some machine-dependence. The different methods used for the determination of hydraulic permeability, covering six orders of magnitude for the sample suite, yield differences in absolute values and pressure dependence for some rocks but not for others. The origin of the differences in permeability, in no case exceeding an order of magnitude, correlate with the compressive strength and potentially reflect a convolution of end plug–sample interaction, sample-to-sample variability, heterogeneity on sample scale, and/or anisotropy, the last two aspects are notably not accounted for by the applied evaluation procedures. Our study provides an extensive data set apt for “benchmarking” considerations, be it regarding new laboratory equipment or numerical modeling approaches.

Environmental Earth Sciences↗
USGS2009

Simulation of Groundwater-Level and Salinity Changes in the Eastern Shore, Virginia

Groundwater-level and salinity changes have been simulated with a groundwater model developed and calibrated for the Eastern Shore of Virginia. The Eastern Shore is the southern part of the Delmarva Peninsula that is occupied by Accomack and Northampton Counties in Virginia. Groundwater is the sole source of freshwater to the Eastern Shore, and demands for water have been increasing from domestic, industrial, agricultural, and public-supply sectors of the economy. Thus, it is important that the groundwater supply be protected from overextraction and seawater intrusion. The best way for water managers to use all of the information available is usually to compile this information into a numerical model that can simulate the response of the system to current and future stresses. A detailed description of the geology, hydrogeology, and historical groundwater extractions was compiled and entered into the numerical model. The hydrogeologic framework is composed of a surficial aquifer under unconfined conditions, a set of three aquifers and associated overlying confining units under confined conditions (the upper, middle, and lower Yorktown-Eastover Formation), and an underlying confining unit (the St. Marys Formation). An estimate of the location and depths of two major paleochannels was also included in the framework of the model. Total withdrawals from industrial, commercial, public-supply, and some agricultural wells were compiled from the period 1900 through 2003. Reported pumpage from these sources increased dramatically during the 1960s and 70s, up to currently about 4 million gallons per day. Domestic withdrawals were estimated on the basis of population census districts and were assigned spatially to the model on the assumption that domestic users are located close to roads. A numerical model was created using the U.S. Geological Survey (USGS) code SEAWAT to simulate both water levels and concentrations of chloride (representing salinity). The model was calibrated using 605 predevelopment and transient water-level observations that are associated predominantly with 20 observation nests of wells sited across the study area. Sampling for groundwater chemistry at these sites revealed that chloride has not increased significantly in the last 20 years. Environmental tracers in the samples also indicated that the water in the surficial aquifer is typically years to decades old, whereas water in the confined aquifers is typically centuries to millennia old. The calibration procedure yielded distributions of hydraulic conductivity and storage coefficients of the aquifers and confining units that are based on 21 pilot points, but vary smoothly across the study area. The estimated values are consistent with other measurements of these properties measured previously on cores and during hydraulic tests at various well fields. Simulations performed with the model demonstrated that the calibrated model can reproduce the observed historical water levels fairly well (R2 = 0.93). The chloride concentrations were also simulated, but a match with chloride concentrations was more difficult to achieve (R2 = 0.16) because of the lack of sufficient data and the unknown exact behavior of the entire transition zone in the millennia leading up to the present day. Future pumping scenarios were simulated through 2050, with pumping set to either 2003 rates or total permitted withdrawal rates. Water levels in 2050 are predicted to be lower than current levels by a few feet where stresses are currently heaviest but potentially by tens of feet if total permitted withdrawals are extracted at current low-stressed sites. Simulations of chloride concentrations through 2050 revealed some potential for seawater intrusion in the areas of Cape Charles, Chincoteague, east of the town of Exmore, and east of the town of Accomac, but precise estimates of concentration increases are highly uncertain. Simulation results were also used to estimate that the down

Scientific Investigations Report↗
USGS2007

An initial-abstraction, constant-loss model for unit hydrograph modeling for applicable watersheds in Texas

Estimation of representative hydrographs from design storms, which are known as design hydrographs, provides for cost-effective, riskmitigated design of drainage structures such as bridges, culverts, roadways, and other infrastructure. During 2001?07, the U.S. Geological Survey (USGS), in cooperation with the Texas Department of Transportation, investigated runoff hydrographs, design storms, unit hydrographs,and watershed-loss models to enhance design hydrograph estimation in Texas. Design hydrographs ideally should mimic the general volume, peak, and shape of observed runoff hydrographs. Design hydrographs commonly are estimated in part by unit hydrographs. A unit hydrograph is defined as the runoff hydrograph that results from a unit pulse of excess rainfall uniformly distributed over the watershed at a constant rate for a specific duration. A time-distributed, watershed-loss model is required for modeling by unit hydrographs. This report develops a specific time-distributed, watershed-loss model known as an initial-abstraction, constant-loss model. For this watershed-loss model, a watershed is conceptualized to have the capacity to store or abstract an absolute depth of rainfall at and near the beginning of a storm. Depths of total rainfall less than this initial abstraction do not produce runoff. The watershed also is conceptualized to have the capacity to remove rainfall at a constant rate (loss) after the initial abstraction is satisfied. Additional rainfall inputs after the initial abstraction is satisfied contribute to runoff if the rainfall rate (intensity) is larger than the constant loss. The initial abstraction, constant-loss model thus is a two-parameter model. The initial-abstraction, constant-loss model is investigated through detailed computational and statistical analysis of observed rainfall and runoff data for 92 USGS streamflow-gaging stations (watersheds) in Texas with contributing drainage areas from 0.26 to 166 square miles. The analysis is limited to a previously described, watershed-specific, gamma distribution model of the unit hydrograph. In particular, the initial-abstraction, constant-loss model is tuned to the gamma distribution model of the unit hydrograph. A complex computational analysis of observed rainfall and runoff for the 92 watersheds was done to determine, by storm, optimal values of initial abstraction and constant loss. Optimal parameter values for a given storm were defined as those values that produced a modeled runoff hydrograph with volume equal to the observed runoff hydrograph and also minimized the residual sum of squares of the two hydrographs. Subsequently, the means of the optimal parameters were computed on a watershed-specific basis. These means for each watershed are considered the most representative, are tabulated, and are used in further statistical analyses. Statistical analyses of watershed-specific, initial abstraction and constant loss include documentation of the distribution of each parameter using the generalized lambda distribution. The analyses show that watershed development has substantial influence on initial abstraction and limited influence on constant loss. The means and medians of the 92 watershed-specific parameters are tabulated with respect to watershed development; although they have considerable uncertainty, these parameters can be used for parameter prediction for ungaged watersheds. The statistical analyses of watershed-specific, initial abstraction and constant loss also include development of predictive procedures for estimation of each parameter for ungaged watersheds. Both regression equations and regression trees for estimation of initial abstraction and constant loss are provided. The watershed characteristics included in the regression analyses are (1) main-channel length, (2) a binary factor representing watershed development, (3) a binary factor representing watersheds with an abundance of rocky and thin-soiled terrain, and (4) curve numb

Scientific Investigations Report↗
USGS1999

Age, double porosity, and simple reaction modifications for the MOC3D ground-water transport model

This report documents modifications for the MOC3D ground-water transport model to simulate (a) ground-water age transport; (b) double-porosity exchange; and (c) simple but flexible retardation, decay, and zero-order growth reactions. These modifications are incorporated in MOC3D version 3.0. MOC3D simulates the transport of a single solute using the method-of-characteristics numerical procedure. The age of ground water, that is the time since recharge to the saturated zone, can be simulated using the transport model with an additional source term of unit strength, corresponding to the rate of aging. The output concentrations of the model are in this case the ages at all locations in the model. Double porosity generally refers to a separate immobile-water phase within the aquifer that does not contribute to ground-water flow but can affect solute transport through diffusive exchange. The solute mass exchange rate between the flowing water in the aquifer and the immobile-water phase is the product of the concentration difference between the two phases and a linear exchange coefficient. Conceptually, double porosity can approximate the effects of dead-end pores in a granular porous media, or matrix diffusion in a fractured-rock aquifer. Options are provided for decay and zero-order growth reactions within the immobile-water phase. The simple reaction terms here extend the original model, which included decay and retardation. With these extensions, (a) the retardation factor can vary spatially within each model layer, (b) the decay rate coefficient can vary spatially within each model layer and can be different for the dissolved and sorbed phases, and (c) a zero-order growth reaction is added that can vary spatially and can be different in the dissolved and sorbed phases. The decay and growth reaction terms also can change in time to account for changing geochemical conditions during transport. The report includes a description of the theoretical basis of the model, a detailed description of input requirements and output options, and the results of model testing and evaluation. The model tests illustrate use of these modifications and demonstrate that accurate solutions can be obtained for these simple cases. Two test cases have no dispersion, illustrating the suitability of this method-of-characteristics model for simulation of advection-dominated transport in ground water.

Water-Resources Investigations Report↗
USGS1996

Developing a temporal database of urban development for the Baltimore/Washington region

The U.S. Geological Survey (USGS), the University of Maryland Baltimore County (UMBC), and the U.S. Bureau of the Census are working together as a multiagency, multidisciplinary team in developing a temporal database that documents the growth of the Baltimore-Washington metropolitan region. This database consists of urban development, principal transportation, shoreline, and population density change. The urban development theme, considered a primary data layer in the study of urban land transformation resulting from human impact on the land, is the focus of this paper. The Baltimore-Washington Spatial Dynamics and Human Impacts Study builds on earlier research efforts that mapped urban land use change for the San Francisco Bay area (Acevedo and Bell, 1994; Bell and others, 1995; Kirtland and others, 1994). In developing a temporal database (Acevedo and others, in press), the team participants hope to provide data that can be used to study patterns of urban growth; assess ecological, environmental, and climatic impacts of urban change; and model and predict future urbanization patterns and impacts (Clarke and others, 1996). Both the San Francisco and Baltimore-Washington regions were selected because of the rapid urban growth and resulting impacts on their ecosystems. The Chesapeake Bay region in particular has undergone extensive environmental agitation due to the hydrologic problems that have arisen from the increase in impermeable surfaces and structures, that is buildings and pavement that physically cover the soil. Because of the inability of water to percolate into the ground, little purification occurs by filtration. Water runs over paved surfaces and quickly washes high levels of toxins directly into the water system. Toxins like gasoline, oil, and fertilizer have dramatically affected the local streams, rivers, and the bay. ABSTRACT The U.S. Geological Survey, the University of Maryland Baltimore County, and the U.S. Bureau of the Census are developing a temporal database to study urban development in the Baltimore-Washington region. The primary data layer, the extent of urban or built-up areas, was compiled using a geographic information system and historical maps, remotely sensed data, digital land use data, and census information from a variety of sources. Urban land use change has been documented by the Baltimore-Washington Spatial Dynamics & Human Impact Study Team for the last 200 years. The methods, definitions, and collection criteria used to define urban or built-up areas were developed by a multi-disciplinary team that also ensures consistency in collection techniques and documentation methods for subsequent application in other regions. Animation techniques were used to visualize the database and to document the evolution of the region's urban landscape. The database is an important tool to urban and regional planners, ecologists, and global change researchers for measuring trends in urban sprawl, analyzing patterns of water pollution, understanding the impacts of development on ecosystems, and developing predictive modeling techniques to better forecast areas of urban growth. This paper describes the techniques used to map the extent of urban areas for Phase I and does not discuss Phase II in detail because the work is still in progress. In this study, urban development is defined as areas of intensive use, with much of the land covered by structures. The built-up areas are characterized by the existence of a systematic street pattern, and the relative concentration of buildings and associated intensive use areas, such as parking lots. Using this definition, urban development does not refer to political boundaries and may include incorporated or unincorporated areas as well as military reservations. ABSTRACT The U.S. Geological Survey, the University of Maryland Baltimore County, and the U.S. Bureau of the Census are developing a temporal database to study urban development in the Baltimore-Washington region. The primary data layer, the extent of urban or built-up areas, was compiled using a geographic information system and historical maps, remotely sensed data, digital land use data, and census information from a variety of sources. Urban land use change has been documented by the Baltimore-Washington Spatial Dynamics & Human Impact Study Team for the last 200 years. The methods, definitions, and collection criteria used to define urban or built-up areas were developed by a multi-disciplinary team that also ensures consistency in collection techniques and documentation methods for subsequent application in other regions. Animation techniques were used to visualize the database and to document the evolution of the region's urban landscape. The database is an important tool to urban and regional planners, ecologists, and global change researchers for measuring trends in urban sprawl, analyzing patterns of water pollution, understanding the impacts of development on ecosystems, and developing predictive modeling techniques to better forecast areas of urban growth. To build the urban component of the temporal database, a multidisciplinary team was assembled and a phased approach initiated. Expanding on procedures developed for the San Francisco Regional Study (Bell and others, 1995), the team developed data definitions, a classification scheme, compilation criteria, mapping specifications, guidelines for source materials, and metadata specifications to support development of a logically consistent dataset. Extensive documentation procedures were established to ensure consistency in data collection, and for subsequent application to other regions. Phase II was the implementation of the regional mapping effort. ABSTRACT The U.S. Geological Survey, the University of Maryland Baltimore County, and the U.S. Bureau of the Census are developing a temporal database to study urban development in the Baltimore-Washington region. The primary data layer, the extent of urban or built-up areas, was compiled using a geographic information system and historical maps, remotely sensed data, digital land use data, and census information from a variety of sources. Urban land use change has been documented by the Baltimore-Washington Spatial Dynamics & Human Impact Study Team for the last 200 years. The methods, definitions, and collection criteria used to define urban or built-up areas were developed by a multi-disciplinary team that also ensures consistency in collection techniques and documentation methods for subsequent application in other regions. Animation techniques were used to visualize the database and to document the evolution of the region's urban landscape. The database is an important tool to urban and regional planners, ecologists, and global change researchers for measuring trends in urban sprawl, analyzing patterns of water pollution, understanding the impacts of development on ecosystems, and developing predictive modeling techniques to better forecast areas of urban growth. The study area for Phase I consisted of an approximate area of 15- by 15-minute segment centered around the city of Baltimore (fig. 1). Phase I was used as a prototype for the technique development and integration that the multiagency collaborative effort would require. The regional study, Phase II, encompassed a 2-degree square centered on Washington, D.C. With more than 7 million people spread across 39 counties, the Baltimore-Washington region is one the Nation's fastest growing metropolitan areas. The two cities are rapidly merging into one.

Maryland↗
USGS2017

Analysis of artificially matured shales with confocal laser scanning raman microscopy: Applications to organic matter characterization

Raman spectroscopy has been suggested as a method for characterizing the thermal maturity of rocks. The literature contains many empirical correlations between thermal maturity proxies, such as vitrinite reflectance (V Ro ) and pyrolysis-T max , with spectral metrics such as Raman peak-widths, peak-center positions, peak-areas and all manner of differences and ratios of these parameters. However, while these correlations may be convincing for small data sets from limited sample series, broader application of these metrics to disparate and heterogeneous samples proves difficult and there remains no consensus. In this extended abstract, Raman spectroscopy is introduced and the history of Raman analysis of carbonaceous material is briefly outlined, highlighting some of the latent difficulties and potential sources of bias. We suggest the organization of a community working group to establish terminology, guidelines, procedures and standards necessary for the successful development of this technique to characterize organic matter in an accessible, unbiased, and reproducible manner. For the present multi-phase study, immature shale samples from the Bakken and Duvernay formations were subjected to hydrous pyrolysis for 72 hours at temperatures from 280°C to 360°C. Rock residues were mounted and polished for analysis via confocal laser-scanning Raman microscopy and reflectance. The maturation series from the Bakken was randomized for the Phase-I single-blind study to be presented at this conference. For the Phase-II study, solid bitumen reflectance (B Ro ) values for the Duvernay series will be known. Multiple hyperspectral maps were collected from each Bakken sample, with each map consisting of a single diffraction-limited spot-size spectrum per 1 µm 2 in rectangular areas several hundred micrometers on a side. Initial attempts at using basic spectral metrics on small numbers of hand-selected spectra to sort the blind series produced inconclusive results: any number of possible correlations could be found. In an improved approach, the statistics of the full spectral datasets were leveraged to: 1) objectively identify organic carbon types (OCTs) in a given map based on Raman and fluorescence spectral characteristics, 2) identify those OCTs in other maps from the same sample and determine if the heterogeneity of the sample has been adequately characterized, and 3) identify the same OCTs in maps from other samples in the maturation series. In ongoing work, our goals are to: 1) use these analyses of the blind series to develop a hypothesis for a correlation to maturation, 2) test the hypothesis by applying the same analyses to the known Duvernay series (in Phase-II), 3) if necessary refine the hypothesis based on observations from the Duvernay analysis, and 4) finally reveal the true order of the Bakken series to verify if the hypothesized correlation accurately predicts the maturity order of the samples. In this document, we share progress to date. The analysis of one area of interest is detailed showing the differentiation of two OCTs based on Raman and fluorescence spectral features, including the use of 2-factor histograms, Principle Components Analysis (PCA), and Nonlinear Iterative Peak Fitting (NIPF).

Conference Paper↗