Search USGSSearch

SEARCH · Search USGS

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 181 records · Page 10Linked to original sources

Spatial aspects of building and population exposure data and their implications for global earthquake exposure modeling

This paper discusses spatial aspects of the global exposure dataset and mapping needs for earthquake risk assessment. We discuss this in the context of development of a Global Exposure Database for the Global Earthquake Model (GED4GEM), which requires compilation of a multi-scale inventory of assets at risk, for example, buildings, populations, and economic exposure. After defining the relevant spatial and geographic scales of interest, different procedures are proposed to disaggregate coarse-resolution data, to map them, and if necessary to infer missing data by using proxies. We discuss the advantages and limitations of these methodologies and detail the potentials of utilizing remote-sensing data. The latter is used especially to homogenize an existing coarser dataset and, where possible, replace it with detailed information extracted from remote sensing using the built-up indicators for different environments. Present research shows that the spatial aspects of earthquake risk computation are tightly connected with the availability of datasets of the resolution necessary for producing sufficiently detailed exposure. The global exposure database designed by the GED4GEM project is able to manage datasets and queries of multiple spatial scales.

Natural Hazards

Cokriging of compositional balances including a dimension reduction and retrieval of original units

Compositional data constitutes a special class of quantitative measurements involving parts of a whole. The sample space has an algebraic-geometric structure different from that of real-valued data. A subcomposition is a subset of all possible parts. When compositional data values include geographical locations, they are also regionalized variables. In the Earth sciences, geochemical analyses are a common form of regionalized compositional data. Ordinarily, there are measurements only at data locations. Geostatistics has proven to be the standard for spatial estimation of regionalized variables but, in general, the compositional character of the geochemical data has been ignored. This paper presents in detail an application of cokriging for the modelling of compositional data using a method that is consistent with the compositional character of the data. The uncertainty is evaluated by a Monte Carlo procedure. The method is illustrated for the contents of arsenic and iron in groundwaters in Bangladesh, which have the peculiarity of being measured in milligrams per litre, units for which the sum of all parts does not add to a constant. Practical results include maps of estimates of the geochemical elements in the original concentration units, as well as measures of uncertainty, such as the probability that the concentration may exceed a given threshold. Results indicate that probabilities of exceedance in previous studies of the same data are too low.

Journal of the Southern African Institute of Minin

Setting up and configuring a total station: Version 1.0: Southeast coast network standard operating procedure 1.2.17

The following standard operating procedure (SOP) outlines the process for setting up and configuring a total station to collect accurate x, y, and z coordinate data. Total stations allow accurate spatial data to be collected and tied to a permanent benchmark. These data can be used to detect small geomorphic changes between site surveys. Many different types of total stations and surveying gear are available, and newer models will be available in the future. This SOP outlines basic steps for using the Trimble S6 robotic total station and the TSC3 data collector. These instructions do not cover detailed care and maintenance of the Trimble S6 robotic total station or the TSC3 data collector.

Southeast Coast Network Standard Operating Procedu

Best practices for distributing and deploying U.S. Geological Survey Shiny applications

A goal of U.S. Geological Survey (USGS) research is to create actionable science for U.S. Department of the Interior managers and partners. Properly managed and released software tools can quickly, accurately, and easily translate scientific results to aid managers’ decision-making. The USGS policies on software development and distribution affect how authors of USGS software products must navigate the review and approval process. Understanding these rules and procedures will save time and effort in making science useful and visible through software development. This document outlines USGS policies on software development as they pertain to Shiny applications, which use R code to create interactive web-based applications. The requirements and procedures described here will also be relevant for other USGS software applications. This document is not intended to be a step-by-step guide on building a Shiny application or to go into detail about any of the referenced software programs. An understanding of the referenced software programs, including R, Git, and hypertext markup language/cascading style sheets (HTML/CSS) is needed to use this document efficiently, though links for more thorough guidance are included throughout.

Techniques and Methods

Construction of a 3-arcsecond digital elevation model for the Gulf of Maine

A system-wide description of the seafloor topography is a basic requirement for most coastal oceanographic studies. The necessary detail of the topography obviously varies with application, but for many uses, a nominal resolution of roughly 100 m is sufficient. Creating a digital bathymetric grid with this level of resolution can be a complex procedure due to a multiplicity of data sources, data coverages, datums and interpolation procedures. This report documents the procedures used to construct a 3-arcsecond (approximately 90-meter grid cell size) digital elevation model for the Gulf of Maine (71°30' to 63° W, 39°30' to 46° N). We obtained elevation and bathymetric data from a variety of American and Canadian sources, converted all data to the North American Datum of 1983 for horizontal coordinates and the North American Vertical Datum of 1988 for vertical coordinates, used a combination of automatic and manual techniques for quality control, and interpolated gaps using a surface-fitting routine.

Maine

Measurement of flow under ice covers in North America

A substantial proportion of natural streams in the United States and Canada are affected by ice cover during the winter. To substantiate the currently used procedures for measuring streamflow during the winter, the U.S. Geological Survey (USGS) and Water Survey of Canada (WSC) began independent, coordinated programs for research and development related to the measurement of streamflow under an ice cover. Detailed measurements of vertical velocity profiles under ice covers in field settings were collected by each agency in accordance with standardized guidelines. The data were then compiled into a joint database. This paper presents a description of the two measurement programs, describes the structure and format of the joint database, and provides preliminary summaries of the data. Ongoing research efforts by the USGS and WSC are described briefly to give examples of the use of the joint database.

Journal of Hydraulic Engineering

Guidelines for determining flood flow frequency: Bulletin #17B of the Hydrology Subcommittee

In December 1967, Bulletin No. 15, "A Uniform Technique for Determining Flood Flow Frequencies," was issued by the Hydrology Committee of the Water Resources Council. The report recommended use of the Pearson Type III distribution with log transformation of the data (log-Pearson Type III distribution) as a base method for flood flow frequency studies. As pointed out in that report, further studies were needed covering various aspects of flow frequency determinations. In March 1976, Bulletin 17, "Guidelines for Determining Flood Flow Frequency" was issued by the Water Resources Council. The guide was an extension and update of Bulletin No. 15. It provided a more complete guide for flood flow frequency analysis incorporating currently accepted technical methods with sufficient detail to promote uniform application. It was limited to defining flood potentials in terms of peak discharge and exceedance probability at locations where a systematic record of peak flood flows is available. The recommended set of procedures was selected from those used or described in the literature prior to 1976, based on studies conducted for this purpose at the Center for Research in Water Resources of the University of Texas at Austin (summarized in Appendix 14) and on studies by the Work Group on Flood Flow Frequency. The "Guidelines" were revised and reissued in June 1977 as Bulletin 17A. Bulletin 17B is the latest effort to improve and expand upon the earlier publications. Bulletin 17B provides revised procedures for weighting a station skew value with the results from a generalized skew study, detecting and treating outliers, making two station comparisons, and computing confidence limits about a frequency curve. The Work Group that prepared this revision did not address the suitability of the original distribution or the generalized skew map. Major problems are encountered when developing guides for flood flow frequency determinations. There is no procedure or set of procedures that can be adopted which, when rigidly applied to the available data, will accurately define the flood potential of any given watershed. Statistical analysis alone will not resolve all flood frequency problems. As discussed in subsequent sections of this guide, elements of risk and uncertainty are inherent in any flood frequency analysis. User decisions must be based on properly applied procedures and proper interpretation of results considering risk and uncertainty. Therefore, the judgment of a professional experienced in hydrologic analysis will enhance the usefulness of a flood frequency analysis and promote appropriate application. It is possible to standarize many elements of flood frequency analysis. This guide describes each major element of the process of defining the flood potential at a specific location in terms of peak discharge and exceedance probability. Use is confined to stations where available records are adequate to warrant statistical analysis of the data. Special situations may require other approaches. In those cases where the procedures of this guide are not followed, deviations must be supported by appropriate study and accompanied by a comparison of results using the recommended procedures. As a further means of achieving consistency and improving results, the Work Group recommends that studies be coordinated when more than one analyst is working currently on data for the same location. This recommendation holds particularly when defining exceedance probabilities for rare events, where this guide allows more latitude. Flood records are limited. As more years of record become available at each location, the determination of flood potential may change. Thus, an estimate may be outdated a few years after it is made. Additional flood data alone may be sufficient reason for a fresh assessment of the flood potential. When making a new assessment, the analyst should incorporate in his study a review of earlier estimates. Where differences appear, they should be acknowledged and explained.

Bulletin

Digital classification of Landsat data for vegetation and land-cover mapping in the Blackfoot River watershed, southeastern Idaho

This paper documents the procedures, results, and final products of a digital analysis of Landsat data used to produce a vegetation and landcover map of the Blackfoot River watershed in southeastern Idaho. Resource classes were identified at two levels of detail: generalized Level I classes (for example, forest land and wetland) and detailed Levels II and III classes (for example, conifer forest, aspen, wet meadow, and riparian hardwoods). Training set statistics were derived using a modified clustering approach. Environmental stratification that separated uplands from lowlands improved discrimination between resource classes having similar spectral signatures. Digital classification was performed using a maximum likelihood algorithm. Classification accuracy was determined on a single-pixel basis from a random sample of 25-pixel blocks. These blocks were transferred to small-scale color-infrared aerial photographs, and the image area corresponding to each pixel was interpreted. Classification accuracy, expressed as percent agreement of digital classification and photo-interpretation results, was 83.0:t 2.1 percent (0.95 probability level) for generalized (Level I) classes and 52.2:t 2.8 percent (0.95 probability level) for detailed (Levels II and III) classes. After the classified images were geometrically corrected, two types of maps were produced of Level I and Levels II and III resource classes: color-coded maps at a 1:250,000 scale, and flatbed-plotter overlays at a 1:24,000 scale. The overlays are more useful because of their larger scale, familiar format to users, and compatibility with other types of topographic and thematic maps of the same scale.

Professional Paper

Procedures for Collecting and Processing Aquatic Invertebrates and Fish for Analysis of Mercury as Part of the National Water-Quality Assessment Program

Mercury studies conducted as part of the U.S. Geological Survey's National Water-Quality Assessment (NAWQA) Program have included nationwide reconnaissance samplings of hundreds of stream sites, as well as detailed, process-oriented research at selected sites. These reconnaissance and detailed studies are intended to provide a better understanding of methylmercury bioaccumulation in stream ecosystems over a range of environmental settings. This publication describes trace-element-clean techniques used for collection and processing of aquatic invertebrates and fish to be analyzed for total mercury, methylmercury, and stable isotopes as part of NAWQA studies.

Open-File Report

New reporting procedures based on long-term method detection levels and some considerations for interpretations of water-quality data provided by the U.S. Geological Survey National Water Quality Laboratory

This report describes the U.S. Geological Survey National Water Quality Laboratory's approach for determining long-term method detection levels and establishing reporting levels, details relevant new reporting conventions, and provides preliminary guidance on interpreting data reported with the new conventions. At the long-term method detection level concentration, the risk of a false positive detection (analyte reported present at the long-term method detection level when not in sample) is no more than 1 percent. However, at the long-term method detection level, the risk of a false negative occurrence (analyte reported not present when present at the long-term method detection level concentration) is up to 50 percent. Because this false negative rate is too high for use as a default 'less than' reporting level, a more reliable laboratory reporting level is set at twice the determined long-term method detection level. For all methods, concentrations measured between the laboratory reporting level and the long-term method detection level will be reported as estimated concentrations. Non-detections will be censored to the laboratory reporting level. Adoption of the new reporting conventions requires a full understanding of how low-concentration data can be used and interpreted and places responsibility for using and presenting final data with the user rather than with the laboratory. Users must consider that (1) new laboratory reporting levels may differ from previously established minimum reporting levels, (2) long-term method detection levels and laboratory reporting levels may change over time, and (3) estimated concentrations are less certain than concentrations reported above the laboratory reporting level. The availability of uncensored but qualified low-concentration data for interpretation and statistical analysis is a substantial benefit to the user. A decision to censor data after they are reported from the laboratory may still be made by the user, if merited, on the basis of the intended use of the data.

Open-File Report

Preliminary report on a study to estimate flood volumes of small rural streams in Ohio: Methods, site selection, and data base

In 1981, the U.S. Geological Survey, in cooperation with the Ohio Department of Transportation and the Federal Highway Administration, began a 7-year flood-volume study of small rural basins in Ohio. This report summarizes the methods of study and describes reconnaissance and site-selection procedures, locations and characteristics of the stations, instrumentation, and methods of collecting and storing data. The first phase of this study involved an intensive field reconnaissance of about 7,000 sites, of which 32 basins were selected for detailed analysis. Drainage areas for the basins varied from 0.13 to 6.45 square miles, and main-channel slopes ranged from 7.6 to 276 feet per mile. Five years of 5-minute rainfall-runoff data will be colledted for each study site. These data will be used to calibrate and verify a rainfall-runoff model for each basin. The calibrated model will be used in conjunction with 80 years of National Weather Service 5-minute precipitation data to synthesize a representative 80-year streamflow record at each site. A Log-Pearson Type III frequency distribution will be applied to each record to define the magnitudes and frequencies of flood volumes at each site. These data will be used to develop regionalized multiple regression models for estimating flood-volume magnitudes and frequencies at small rural ungaged sites in Ohio. The report also summarizes rainfall-runoff data collected from July 1981 through September 1983, but does not interpret the data. An average of eleven event periods per site were monitored where maximum 5-minute rainfall intensities varied from 0.02 to .067 inches and maximum peak discharges varied from 1 to 1,130 cubic feet per second.

Ohio

Levels at streamflow gaging stations

This manual establishes the surveying procedures for setting gages at a streamflow gaging station to datum and for checking them periodically for errors caused by vertical movement of the gage-supporting structures. The surveying terms and concepts used are explained; and the details of testing, adjusting, and operating the instruments are outlined. Notekeeping, adjusting level circuits, checking gages, summarizing results, locating the nearest National Geodetic Vertical Datum of 1929 bench mark, and relating the gage datum to the national datum are described.

Open-File Report

Organic and isotopic indicators for sorting of sedimentary organic matter along a marginal submarine canyon

Submarine canyons are incised features of many continental margins that can have significant influence on the hydrodynamic distribution of sediments and organic matter eroded and deposited from the continents. Baltimore Canyon, on the mid-Atlantic margin of the United States, contains a complex set of sedimentary processes that simultaneously create unique benthic habitats and control the deposition of organic matter. Along the canyon axis, loci of net erosion, net deposition, and intense winnowing each host diverse faunal assemblages and varying mixtures of sedimentary organic matter derived both from production in the overlying water column and from mobilized sediments. Bioavailable components of this deposited organic matter sustain benthic communities, while recalcitrant components can contribute to long-term carbon burial in the deep sea. However, commonly employed bulk geochemical analyses provide little information about the relative bioavailability or depositional history of sedimentary organic matter. Here we employ a range of organic and isotopic analyses to explore in more detail how canyon-specific sediment dynamics determine the sorting of organic matter from shelf to open ocean. In combination with bulk geochemical characteristics, we subjected surface sediments from water depths of ∼200–1200 m in Baltimore Canyon to a sequential extraction procedure, isolating nonpolar and polar lipid classes, an acid-soluble fraction, and an acid-insoluble fraction. Each class was analyzed for carbon and nitrogen quantities and stable isotope ratios, and radiocarbon content where possible, along with compound-specific carbon and nitrogen isotope analysis of individual amino acids in the acid-hydrolysed fraction. We find different organic matter sources and depositional history recorded in the properties of younger, bioavailable organic matter components (polar lipids, amino acids) in comparison to the older, more recalcitrant components (nonpolar lipids, acid-insoluble fraction). These differences in source and bioavailability of organic matter vary along the canyon, correlating with grain size and erosion/deposition dynamics, and may help shape the benthic faunal assemblages. Additionally, our results suggest that determining the relative concentrations of acid-soluble and acid-insoluble organic matter may provide an easily accessible method to improve our understanding of the nutritional quality of sediments for benthic fauna than more commonly used bulk carbon or nitrogen concentrations.

Delaware, New Jersey

Continuous fields of land cover for the conterminous United States using Landsat data: First results from the Web-Enabled Landsat Data (WELD) project

Vegetation Continuous Field (VCF) layers of 30 m percent tree cover, bare ground, other vegetation and probability of water were derived for the conterminous United States (CONUS) using Landsat 7 Enhanced Thematic Mapper Plus (ETM+) data sets from the Web-Enabled Landsat Data (WELD) project. Turnkey approaches to land cover characterization were enabled due to the systematic WELD Landsat processing, including conversion of digital numbers to calibrated top of atmosphere reflectance and brightness temperature, cloud masking, reprojection into a continental map projection and temporal compositing. Annual, seasonal and monthly WELD composites for 2008 were used as spectral inputs to a bagged regression and classification tree procedure using a large training data set derived from very high spatial resolution imagery and available ancillary data. The results illustrate the ability to perform Landsat land cover characterizations at continental scales that are internally consistent while retaining local spatial and thematic detail.

Remote Sensing Letters

Oil, biological communities and contingency planning

The Oil Pollution Act of 1990 mandates the inclusion of a fish and wildlife response plan in the National Contingency Plan (NCP) and the creation of Area Committees that must develop an Area Contingency Plan (ACP). Area Contingency Plans must include a detailed annex containing a Fish and Wildlife and Sensitive Environments Plan. Tank vessels, offshore facilities, and certain onshore facilities must have response plans consistent with the requirements of the NCP and the ACP. New regulations to supersede the Type A and B procedures of the Natural Resources Damage Assessment Regulations are being developed for oil spills. Currently, four assessment methods have been proposed: (1) Type A, (2) comprehensive (Type B), (3) intermediate (between types A and B), and (4) compensation tables. The Oil Spill Liability Trust Fund is approaching its ceiling of $1 billion, but only $50 million has been appropriated. Effective biological contingency planning requires extensive knowledge of (1) the environmental fate of petroleum, (2) the effects of petroleum on organisms, (3) the existing biological resources, and (4) the establishment of a system of biological priorities. The characteristics and fate of petroleum and the biological effects of petroleum are reviewed. Assessment of biological resources includes plant and animal distributions, important habitat, endangered or threatened species, and economic considerations. The establishment by Area Committees of priorities for environmental protection, injury assessment, and restoration will promote efficient spill response. Three special issues are discussed: (1) improving our ability to restore natural resources, (2) the potential role of biological diversity in spill response planning, and (3) planning for animal rehabilitation.

Book chapter

Flood-frequency analyses from paleoflood investigations for Spring, Rapid, Boxelder, and Elk Creeks, Black Hills, western South Dakota

Flood-frequency analyses for the Black Hills area are important because of severe flooding of June 9-10, 1972, that was caused by a large mesoscale convective system and caused at least 238 deaths. Many 1972 peak flows are high outliers (by factors of 10 or more) in observed records that date to the early 1900s. An efficient means of reducing uncertainties for flood recurrence is to augment gaged records by using paleohydrologic techniques to determine ages and magnitudes of prior large floods (paleofloods). This report summarizes results of paleoflood investigations for Spring Creek, Rapid Creek (two reaches), Boxelder Creek (two subreaches), and Elk Creek. Stratigraphic records and resulting long-term flood chronologies, locally extending more than 2,000 years, were combined with observed and adjusted peak-flow values (gaged records) and historical flood information to derive flood-frequency estimates for the six study reaches. Results indicate that (1) floods as large as and even substantially larger than 1972 have affected most of the study reaches, and (2) incorporation of the paleohydrologic information substantially reduced uncertainties in estimating flood recurrence. Canyons within outcrops of Paleozoic rocks along the eastern flanks of the Black Hills provided excellent environments for (1) deposition and preservation of stratigraphic sequences of late-Holocene flood deposits, primarily in protected slack-water settings flanking the streams; and (2) hydraulic analyses for determination of associated flow magnitudes. The bedrock canyons ensure long-term stability of channel and valley geometry, thereby increasing confidence in hydraulic computations of ancient floods from modern channel geometry. Stratigraphic records of flood sequences, in combination with deposit dating by radiocarbon, optically stimulated luminescence, and cesium-137, provided paleoflood chronologies for 29 individual study sites. Flow magnitudes were estimated from elevations of flood deposits in conjunction with hydraulic calculations based on modern channel and valley geometry. Reach-scale paleoflood chronologies were interpreted for each study reach, which generally entailed correlation of flood evidence among multiple sites, chiefly based on relative position within stratigraphic sequences, unique textural characteristics, or results of age dating and flow estimation. The FLDFRQ3 and PeakfqSA analytical models (assuming log-Pearson Type III frequency distributions) were used for flood-frequency analyses for as many as four scenarios: (1) analysis of gaged records only; (2) gaged records with historical information; (3) all available data including gaged records, historical flows, paleofloods, and perception thresholds; and (4) the same as the third scenario, but ?top fitting? the distribution using only the largest 50 percent of gaged peak flows. The PeakfqSA model is most consistent with procedures adopted by most Federal agencies for flood-frequency analysis and thus was (1) used for comparisons among results for study reaches, and (2) considered by the authors as most appropriate for general applications of estimating low-probability flood recurrence. The detailed paleoflood investigations indicated that in the last 2,000 years all study reaches have had multiple large floods substantially larger than in gaged records. For Spring Creek, stratigraphic records preserved a chronology of at least five paleofloods in approximately (~) 1,000 years approaching or exceeding the 1972 flow of 21,800 cubic feet per second (ft 3 /s). The largest was ~700 years ago with a flow range of 29,300-58,600 ft 3 /s, which reflects the uncertainty regarding flood-magnitude estimates that was incorporated in the flood-frequency analyses. In the lower reach of Rapid Creek (downstream from Pactola Dam), two paleofloods in ~1,000 years exceeded the 1972 flow of 31,200 ft 3 /s. Those occurred ~440 and 1,000 years ago, with flows of 128,000-256,000 and 64,000-128,000 ft 3 /s, respectively. Five smaller paleofloods of 9,500-19,000 ft 3 /s occurred between ~200 and 400 years ago. In the upper reach of Rapid Creek (above Pactola Reservoir), the largest recorded floods are substantially smaller than for lower Rapid Creek and all other study reaches. Paleofloods of ~12,900 and 12,000 ft 3 /s occurred ~1,000 and 1,500 years ago. One additional paleoflood (~800 years ago) was similar in magnitude to the largest gaged flow of 2,460 ft 3 /s Boxelder Creek was treated as having two subreaches because of two tributaries that affect peak flows. During the last ~1,000 years, paleofloods of ~39,000-78,000 ft 3 /s and 40,000-80,000 ft 3 /s in the upstream subreach have exceeded the 1972 peak flow of 30,800 ft 3 /s. One other paleoflood was similar to the second largest gaged flow (16,400 ft 3 /s in 1907). For the downstream subreach, paleofloods of 61,300-123,000 ft 3 /s and 52,500-105,000 ft 3 /s in the last ~1,000 years have substantially exceeded the 1972 flood (50,500 ft 3 /s). Four additional paleofloods had flows between 14,200 and 33,800 ft 3 /s. The 1972 flow on Elk Creek (10,400 ft 3 /s) has been substantially exceeded at least five times in the last 1,900 years. The largest paleoflood (41,500-124,000 ft 3 /s) was ~900 years ago. Three other paleofloods between 37,500 and 120,000 ft 3 /s occurred between 1,100 and 1,800 years ago. A fifth paleoflood of 25,500-76,500 ft 3 /s was ~750 years ago. Considering analyses for all available data (PeakfqSA model) for all six study reaches, the 95-percent confidence intervals about the low-probability quantile estimates (100-, 200-, and 500-year recurrence intervals) were reduced by at least 78 percent relative to those for the gaged records only. In some cases, 95-percent uncertainty intervals were reduced by 99 percent or more. For all study reaches except the two Boxelder Creek subreaches, quantile estimates for these long-term analyses were larger than for the short-term analyses. The 1972 flow for the Spring Creek study reach (21,800 ft 3 /s) corresponds with a recurrence interval of ~400 years. Recurrence intervals are ~500 years for the 1972 flood magnitudes along the lower Rapid Creek reach and the upstream subreach of Boxelder Creek. For the downstream subreach of Boxelder Creek, the large 1972 flood magnitude (50,500 ft 3 /s) exceeds the 500-year quantile estimate by about 35 percent. The recurrence interval of ~100 years for 1972 flooding along the Elk Creek study reach is small relative to other study reaches along the eastern margin of the Black Hills. All of the paleofloods plot within the bounds of a national envelope curve, indicating that the national curve represents exceedingly rare floods for the Black Hills area. Elk Creek, lower Rapid Creek, and the downstream subreach of Boxelder Creek all have paleofloods that plot above a regional envelope curve; in the case of Elk Creek, by a factor of nearly two. The Black Hills paleofloods represent some of the largest known floods, relative to drainage area, for the United States. Many of the other largest known United States floods are in areas with physiographic and climatologic conditions broadly similar to the Black Hills-semiarid and rugged landscapes that intercept and focus heavy precipitation from convective storm systems. The 1972 precipitation and runoff patterns, previous analyses of peak-flow records, and the paleoflood investigations of this study support a hypothesis of distinct differences in flood generation within the central Black Hills study area. The eastern Black Hills are susceptible to intense orographic lifting associated with convective storm systems and also have high relief, thin soils, and narrow and steep canyons-factors favoring generation of exceptionally heavy rain-producing thunderstorms and promoting runoff and rapid concentration of flow into stream channels. In contrast, storm potential is smaller in and near the Limestone Plateau area, and storm runoff is further reduced by substantial infiltration into the limestone, gentle topography, and extensive floodplain storage. Results of the paleoflood investigations are directly applicable only to the specific study reaches and in the case of Rapid Creek, only to pre-regulation conditions. Thus, approaches for broader applications were developed from inferences of overall flood-generation processes, and appropriate domains for application of results were described. Example applications were provided by estimating flood quantiles for selected streamgages, which also allowed direct comparison with results of at-site flood-frequency analyses from a previous study. Several broad issues and uncertainties were examined, including potential biases associated with stratigraphic records that inherently are not always complete, uncertainties regarding statistical approaches, and the unknown applicability of paleoflood records to future watershed conditions. The results of the paleoflood investigations, however, provide much better physically based information on low-probability floods than has been available previously, substantially improving estimates of the magnitude and frequency of large floods in these basins and reducing associated uncertainty.

South Dakota

Comparison of groundwater-model construction methods, representations of glacial geology, model designs, and groundwater-model flow simulations within Elkhart County, Indiana

Automated data-processing methods allow hydrologists to efficiently incorporate digital well-record datasets into the construction of hydrostratigraphic frameworks for groundwater-flow models. The method selected to construct the hydrostratigraphic framework can affect the extent of geologic heterogeneity that can be included in the model. The detail generated from a hydrostratigraphic framework can affect groundwater simulation results. The effects of detail on model accuracy, groundwater-flow simulations, and particle-tracking simulations are described in this study. This report compares differences in hydrostratigraphic frameworks and results of groundwater models using (1) a method that incorporates more hydrologic judgment at the expense of using limited lithologic data and (2) a method that is more automated and uses all available lithologic data. The study additionally evaluates the effect of model discretization and inclusion of more (or less) geologic detail on simulation results. Two methods were used to create hydrostratigraphic frameworks of glacial deposits in the St. Joseph River Basin. One method, referred to as the subjective method, manually identifies stratigraphic boundaries using a sample of well logs from State databases and uses two-dimensional kriging to create three model layers of the study area. Indicator kriging is used to define aquifer extent in each layer. The second method, referred to as the objective method, uses three-dimensional kriging to automatically create a detailed heterogeneous model of the study area using all wells logs from the State database. The objective method increases detail in the vertical by greatly increasing the number of computer groundwater model layers from 3 to 30. In Elkhart County, Indiana, a previously published model represents the product of the subjective method, and a newly calibrated model of the same area represents the product of the objective method. An automated calibration procedure was used with the objective model (derived from the objective method) for Elkhart County. The two most-sensitive parameters for the Elkhart County objective model are horizontal hydraulic conductivity of the sand and the combined sand and gravel/gravel deposits. Vertical hydraulic conductivity of the fine-grained and intermediate-sized deposits could not be estimated, possibly indicating major flow paths are along a continuously connected series of sand and gravel deposits and not through a confining layer. The statistics measuring model calibration accuracy for the objective model were slightly better than statistics for the subjective model (model derived from the subjective method) of Elkhart County, but the hydraulic conductivities and flow rates for the two models were different. The mean absolute errors between simulated and measured groundwater levels are 2.04 and 2.16 feet for the objective and subjective models, respectively. Simulated seepage losses from and groundwater discharges to measured stream reaches in the objective model were evenly balanced in terms of over and under simulations of measured values; the subjective model tended to overpredict measured groundwater discharge to streams. The overprediction may be related to the 58 percent greater total inflow and outflow through the subjective model. The greater flow rate through the subjective model results from higher horizontal hydraulic conductivities in the subjective model than in the objective model. Horizontal hydraulic conductivity ranged from 23.9 to 111 feet per day in the objective model and generally ranged from 170 to 370 feet per day in the subjective model. The improvement in calibration statistics for the objective model relative to the subjective model may be from increased detail in how the objective model represents the distribution of fine- and coarse-grained deposits. The improvement also could be associated with the difference in methods used to represent the continuity of the confining unit. The effect of differences in horizontal hydraulic conductivity distributions between the two models for Elkhart County is evident in the groundwater-flow paths simulated by the objective and subjective models. At a withdrawal well location, the flow lines produced by the objective model indicate a wider contributing area than that for the subjective model. The discontinuous confining unit represented in the objective model provided the opportunity for groundwater flow to split into an upper and lower path. The split in flow simulated by the objective model at one location was independently supported by bromide concentrations in groundwater; the subjective model did not duplicate the split in flow.

Indiana

Distribution of arsenic, selenium, and other trace elements in high pyrite Appalachian coals: evidence for multiple episodes of pyrite formation

Pennsylvanian coals in the Appalachian Basin host pyrite that is locally enriched in potentially toxic trace elements such as As, Se, Hg, Pb, and Ni. A comparison of pyrite-rich coals from northwestern Alabama, eastern Kentucky, and West Virginia reveals differences in concentrations and mode of occurrence of trace elements in pyrite. Pyrite occurs as framboids, dendrites, or in massive crystalline form in cell lumens or crosscutting veins. Metal concentrations in pyrite vary over all scales, from microscopic to mine to regional, because trace elements are inhomogeneously distributed in the different morphological forms of pyrite, and in the multiple generations of sulfide mineral precipitates. Early diagenetic framboidal pyrite is usually depleted in As, Se, and Hg, and enriched in Pb and Ni, compared to other pyrite forms. In dendritic pyrite, maps of As distribution show a chemical gradient from As-rich centers to As-poor distal branches, whereas Se concentrations are highest at the distal edges of the branches. Massive crystalline pyrite that fills veins is composed of several generations of sulfide minerals. Pyrite in late-stage veins commonly exhibits As-rich growth zones, indicating a probable epigenetic hydrothermal origin. Selenium is concentrated at the distal edges of veins. A positive correlation of As and Se in pyrite veins from Kentucky coals, and of As and Hg in pyrite-filled veins from Alabama coals, suggests coprecipitation of these elements from the same fluid. In the Kentucky coal samples (n = 18), As and Se contents in pyrite-filled veins average 4200 ppm and 200 ppm, respectively. In Alabama coal samples, As in pyrite-filled veins averages 2700 ppm (n = 34), whereas As in pyrite-filled cellular structures averages 6470 ppm (n = 35). In these same Alabama samples, Se averages 80 ppm in pyrite-filled veins, but was below the detection limit in cell structures. In samples of West Virginia massive pyrite, As averages 1700 ppm, and Se averages 270 ppm (n = 24). The highest concentration of Hg (≤ 102 ppm) is in Alabama pyrite veins. Improved detailed descriptions of sulfide morphology, sulfide mineral paragenesis, and trace-element concentration and distribution allow more informed predictions of: (1) the relative rate of release of trace elements during weathering of pyrite in coals, and (2) the relative effectiveness of various coal-cleaning procedures of removing pyrite. For example, trace element-rich pyrite has been shown to be more soluble than stoichiometric pyrite, and fragile fine-grained pyrite forms such as dendrites and framboids are more susceptible to dissolution and disaggregation but less amenable to removal during coal cleaning.

Alabama;Kentucky;West Virginia