SEARCH · USGS Science

Results for “Final Report”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,693 records · Page 94Linked to original sources

USGS1985

Tectonic implications of the Indian Run Formation— A newly recognized sedimentary mélange in the northern Virginia Piedmont

Sedimentary melange in the northeastern part of Fairfax County, Virginia, contains both mesoscopic and mappable fragments of Accotink Schist, Lake Barcroft Metasandstone, metagabbro, and ultramafic rocks as well as smaller fragments of other rock types. This melange was originally mapped as the Sykesville Formation, a major precursory sedimentary melange in northern Virginia and Maryland. The fragments of Accotink Schist and Lake Barcroft Metasandstone within the Sykesville were considered to be rip-ups of these units over which the Sykesville slid when finally emplaced. More recent study has shown that fragments of Accotink and Lake Barcroft are restricted to a certain area of sedimentary melange originally defined as Sykesville, and this part of the melange is now considered to be a separate mappable unit, here named the Indian Run Formation. The Indian Run underlies the sequence Accotink Schist and Lake Barcroft Metasandstone which is here formally named the Annandale Group. The Indian Run is intruded by the Occoquan Granite of Cambrian age, so it is of Cambrian or Late Proterozoic age. The Sykesville Formation (restricted) is a much more extensive unit than the Indian Run Formation and is characterized by its contained olistoliths of the Peters Creek Schist, the unit that tectonically overlies it. The Sykesville and Peters Creek constitute a precursory melange-allochthon pair which is here termed a 'tectonic motif.' The Indian Run-Annandale pair then forms a tectonically lower motif, and the overlying pair, the Yorkshire Formation-Piney Branch Complex, forms a tectonically higher motif. The Chopawamsic Formation and underlying sedimentary melange in the area south of Fairfax County may form a tectonic motif beneath the Indian Run-Annandale tectonic motif. Thus, three and perhaps four repetitions of precursory melange-allochthon pairs occur in northern Virginia. Other percursory melanges and motifs may occur in the Maryland Piedmont to the north. The tectonic setting of the motif formation and assemblage is uncertain at this time. A model involving the obduction of several separate sheets onto the ancestral North American continental margin is appealing in that it involves the closing and destruction of a marginal basin, a relatively simple concept. This model fails, however, to supply a source for the sediment necessary to form the precursory melanges. A trench-slope origin would supply the vast amount of needed sediment by accretion. The precursory melangeallochthon motifs would then be stacked near the base of the trench slope. This model is appealing and requires a rather complicated assemblage of continental, arc, and oceanic fragments, such as those that occur on many modern continental margins.

Virginia↗

USGS1973

Glacial and postglacial geologic history of Isle Royale National Park, Michigan

Isle Royale was overridden by glacial ice during each of the four major glaciations of the Pleistocene Epoch, and each successive glaciation essentially obliterated all direct evidence of preceding glaciations on the island. In the waning phase of the last major glaciation, the Wisconsin Glaciation, the frontal ice margin retreated northward from at least the greater part of the Lake Superior basin, then readvanced into the basin during Valders time, about 11,000 years ago. We can attribute to the Valders ice the final aspect of glaciation on Isle Royale, including both erosional and depositional features. It is impossible to estimate the quantity of glacial debris or other surficial materials that might have been present on Isle Royale prior to the Valders readvance, but the readvancing ice appears to have removed most of what might have been present, as judged by the thin surficial cover on the eastern two-thirds of the island today. During the Valders retreat, a series of lakes formed in the Lake Superior basin in front of the retreating ice margin. The retreating ice opened successively lower outlets, and thus the general trend of lake elevations is downward. Distinct lake stages reflect periods of relative stability during which well- defined shoreline features developed. The ice front forming the north margin of the earlier lakes probably remained south of Isle Royale until about the time of glacial Lake Beaver Bay, when it retreated to a position straddling Isle Royale west of Lake Desor. Abundant deposits of glacial debris were left upon the newly deglaciated west end of the island, and the ice front remained stable long enough to build a complex of ice-margin deposits across the island. Shorelines formed by the glacial lake associated with this ice front are found on the western part of the island about 200 feet above present Lake Superior. Subsequent renewed and complete retreat of the ice margin from Isle Royale was rapid enough that only a minor amount of glacial debris was deposited on the central and eastern parts of the island. When the ice margin reached the north edge of the Lake Superior basin, Lake Minong was formed, and the entire basin was filled for the first time since the Valders advance. Lake Minong marked a relatively stable episode in the history of the basin, and its beaches are among the best developed of the abandoned shoreline features on Isle Royale. Lake Minong beaches and later lower beaches are best developed on the southwest end of Isle Royale, where abundant glacial debris provided easily worked materials for beach construction.

Michigan↗

USGS2017

Geology and mining history of the Southeast Missouri Barite District and the Valles Mines, Washington, Jefferson, and St. Francois Counties, Missouri

The Southeast Missouri Barite District and the Valles Mines are located in Washington, Jefferson, and St. Francois Counties, Missouri, where barite and lead ore are present together in surficial and near-surface deposits. Lead mining in the area began in the early 1700’s and extended into the early 1900’s. Hand mining of lead in the residuum resulted in widespread pits (also called shafts or diggings), and there was some underground mining of lead in bedrock. By the 1860’s barite was recovered from the residuum by hand mining, also resulting in widespread diggings, but generally not underground mines in bedrock. Mechanized open-pit mining of the residuum for barite began in the 1920’s. Barite production slowed by the 1980’s, and there has not been any barite mining since 1998. Mechanized barite mining resulted in large mined areas and tailings ponds containing waste from barite mills. The U.S. Environmental Protection Agency (EPA) has determined that lead is present in surface soils in Washington and Jefferson Counties at concentrations exceeding health-based screening levels. Also, elevated concentrations of barium, arsenic, and cadmium have been identified in surface soils, and lead concentrations exceeding the Federal drinking-water standard of 15 micrograms per liter have been identified in private drinking-water wells. Potential sources of these contaminants are wastes associated with barite mining, wastes associated with lead mining, or unmined natural deposits of barium, lead, and other metals. As a first step in helping EPA determine the source of soil and groundwater contamination, the U.S. Geological Survey (USGS), in cooperation with the EPA, investigated the geology and mining history of the Southeast Missouri Barite District and the Valles Mines. Ore minerals are barite (barium sulfate), galena (lead sulfide), cerussite (lead carbonate), anglesite (lead sulfate), sphalerite (zinc sulfide), smithsonite (zinc carbonate), and chalcopyrite (copper-iron sulfide). The Cambrian Potosi Dolomite is the most important formation for the ore deposits, followed by the Eminence Dolomite. Because galena, sphalerite, and barite are less soluble than dolomite, chemical weathering of the ore-bearing dolomite bedrock resulted in the concentration of ore minerals in the residuum. Most of the barite and lead mining was in the residuum, which averages 10 to 15 feet thick. Lead mining by French explorers may have begun in 1719 along Old Mines Creek at Cabanage de Renaudiere, which was followed shortly by the discovery of lead and the development of lead mines at Mine Renault (also called Forche a Renault Mine), Old Mines, and at other places along the Big River, Mineral Fork, and Forche a Renault Creek. Lead mining began sometime between 1775 and 1780 at Mine a Breton, the name of which was later changed to Potosi. Other mining areas were developed in the early part of the 19th century, including Fourche a Courtois (Palmer Mines), the French Diggings, and the Richwoods Mines. Zinc became a valuable resource after the Civil War, and the Valles Mines was an important supplier of zinc as well as lead, with at least some production up until the 1920’s. Lead mining declined in the early part of the 20th century as mining in the Old Lead Belt, Mine La Motte, and the Tri-State District expanded. The earliest lead mines were diggings in the residuum and were round holes (shafts) about 4 feet in diameter dug with pick and shovel about 15–20 feet deep, with drifts dug a short distance laterally from the bottom of the shafts. This mining process was repeated a short distance away until a large area was covered with pits. Some mining in bedrock began by about 1800, with shafts as deep as 170 feet and as much as several hundred feet of lateral drifts. Smelting of the lead ore to elemental lead was first done using a log furnace, which was inefficient; estimates have been made that only about 50 percent of the lead was recovered, and the remainder was lost to the ashes (slags) and to volatilization. Starting in 1798, ash furnaces were used to smelt the ashes from the log furnaces. These two furnaces were worked in tandem for many years but were gradually replaced by other furnaces, including the Scotch hearth. Estimates of lead recovery as high as 80–90 percent have been made for the Scotch hearth. By the mid-1870’s the air furnace was being used, also with estimated lead recovery as high as 80–90 percent. Zinc furnaces were built when zinc became a valuable commodity, but much of the zinc ore was shipped out of the area, either to a smelter in St. Louis, Missouri, or to other smelters. The total lead and zinc production from the Southeast Missouri Barite District and the Valles Mines is estimated at 180,000 tons of lead and 60,000 tons of zinc. An estimated 97,000 tons of lead and an estimated 120,000 tons of zinc were lost during smelting. The estimated losses do not include losses at the mine site during mining and preparation for smelting, such as the loss of fine-grained galena during hand cleaning or the discarding of zinc ore before its value was known, for which no estimates are available. Hand mining for barite in the residuum was active by at least the 1860’s and peaked from 1905 to the 1930’s when several thousand people were engaged in barite mining. Hand mining (diggings) and cleaning of the ore was done in much the same way as earlier lead mining, with the additional use of a rattle box to further clean the barite. Mechanized open-pit mining of old barite diggings began in 1924 to recover barite left behind by hand mining, and washing plants were used to clean the clay from the barite. Hand mining, however, continued to thrive, and washer plants began to close temporarily in 1931; nearly all of the barite produced before 1937 was by hand mining. By the 1940’s, however, all barite mining was mechanized. Mechanized mining used shovels powered by steam, gasoline, or electricity (and by the 1950’s draglines and front-end loaders) to mine the residuum. The ore was loaded onto rail cars (and by the 1940’s, trucks) for shipment to washer plants. Clay was removed from the barite using a log washer, and a jig was used to concentrate the barite. Overflow from the log washers was waste and went to a mud (tailings) pond. The coarse jig tailings went to tailings piles or were used as railroad ballast and, later, to create roads within the mine pit. Some barite was ground, depending on its final use, and some ground barite was bleached using a hot solution of sulfuric acid to remove impurities such as iron minerals and lead sulfide (galena). An earlier bleaching process used lead-lined tanks. Large quantities of water were required for milling the barite; some was recirculated water and the remainder came from dammed streams or was pumped from wells. Tailings and wastewater were impounded behind dikes that were built across small valleys and were increased in height as necessary using washer waste and any overburden that had been stripped. In some cases, dikes were built across valleys that had already been mined for barite. The total production of barite from the Southeast Missouri Barite District and the Valles Mines is estimated to have been about 13.1 million tons. Most of the barite production was from Washington County. Hand mining and processing of barite was inefficient. Estimates of barite recovery range from less than one-fourth to about one-half because pillars between the shafts in the residuum needed to be left unmined for stability. With mechanized mining, large amounts of barite were lost during the milling process. It has been estimated that about 30 percent of the barite was lost and that about two-thirds of the lost barite was fine-grained and was discharged to the tailings ponds. Some galena was lost to the tailings ponds. A 1972 inventory of tailings ponds by the Missouri Geological Survey identified 67 ponds in the Southeast Missouri Barite District (there are more than this currently documented). Results from samples from four ponds that were drilled were used to estimate that the 67 ponds contained almost 39 million tons (or cubic yards) of tailings averaging about 5 percent barite, for a potential reserve of 1.935 million tons of barite. It is not known how much lead was removed during barite mining, either by hand or mechanized mining and processing, how much lead was recovered, or how much lead went as fines to the tailing ponds or as coarse material to mine roads or was otherwise lost.

Missouri↗

USGS1989

Assessing the relationship between Section 404 and wetland losses: a feasibility study

The primary objective of the Clean Water Act of 1977 (33 U.S.C. 1251) is to restore and maintain the physical, chemical, and biological integrity of the Nation's waters. Section 404 of the Act regulates the discharge of dredged or fill materials into wetlands and represents the primary Federal authority for regulation of wetland alterations. Since its inception, the Section 404 program has been controversial in regard to the extent to which it was intended to provide wetlands regulation. Section 404 requires those who wish to discharge dredged or fill material into waters of the United States, which include many wetlands, to first obtain a Federal permit. The Environmental Protection Agency (EPA) has overall responsibility for administration of the Section 404 program and promulgates guidelines that must be followed in issuing permits. In addition, EPA has the final authority to prohibit specific discharges if the environmental impacts are unacceptable. The U.S. Army Corps of Engineers (Corps) issues Section 404 permits, which can be of two types. Individual Permits are issued following case-by-case reviews of proposed discharges. General Permits, which can be either nationwide or regional in scope, are authorized by the Corps for categories of activities that are similar in nature and that have only minimal individual and cumulative adverse environmental impacts. EPA, the National Marine Fisheries Service (NMFS), the Fish and Wildlife Service (Service), and State natural resource agencies review and comment on permit applications and offer recommendations on appropriate mitigation measures. Although comments from the Service and other natural resource agencies are advisory in nature (EPA's veto authority excepted), they can serve as the basis for modifying, conditioning, or denying a Section 404 permit. In 1986, in a survey conducted by the National Ecology Research Center, Service personnel indicated interest in additional information concerning both wetland trends and the impacts of activities authorized by the Corps under Section 404 (Roelle 1986). Although there is some information concerning wetland losses for certain geographic areas and for the Nation as a whole (Frayer et al. 1983; Tiner 1984), there appears to be little information on how these losses relate to the Section 404 permitting process. The primary objective of this study was to determine the feasibility of estimating wetland losses in relationship to Individual and General Permits issued under Section 404. A secondary objective was to assemble data on acceptance and implementation of specific mitigation recommendations offered by the Service and other natural resource agencies in connection with development activities on wetlands examined. At present, 26 categories of Nationwide permits have been authorized by the Corps. Nationwide permit 26 was of particular interest in this study because it specifically authorizes discharges into wetlands under certain circumstances (see Figure 1 and related text for a more complete discussion of circumstances under which Nationwide Permit 26 is applicable). All subsequent references to Nationwide or General permits pertain to Nationwide Permit 26.

Report↗

USGS2015

Impacts of fire management on aboveground tree carbon stocks in Yosemite and Sequoia & Kings Canyon National Parks

Forest biomass on Sierra Nevada landscapes constitutes one of the largest carbon stocks in California, and its stability is tightly linked to the factors driving fire regimes. Research suggests that fire suppression, logging, climate change, and present management practices in Sierra Nevada forests have altered historic patterns of landscape carbon storage, and over a century of fire suppression and the resulting accumulation in surface fuels have been implicated in contributing to recent increases in high severity, stand-replacing fires. For over 30 years, fire management at Yosemite (YOSE) and Sequoia & Kings Canyon (SEKI) national parks has led the nation in restoring fire to park landscapes; however, the impacts on the stability and magnitude of carbon stocks have not been thoroughly examined. The purpose of this study is to quantify relationships between recent fire patterns and aboveground tree carbon stocks in YOSE and SEKI. Our approach focuses on evaluating fire effects on 1) amounts of aboveground tree carbon on the landscape, and 2) rates of carbon accumulation by individual trees. In 2010, we compiled a database of existing plot data for our analyses. In 2011, our field crews acquired vegetation data and collected tree growth cores from 105 plots. In 2012, we completed an interpretive component and began data analyses. In 2013, processing of tree cores began. In 2014, final processing of tree cores, data analyses, and manuscript preparation was conducted. The work for this project was facilitated through an interagency agreement between the National Park Service and the U.S. Geological Survey, and through a Cooperative Ecosystems Studies Unit (CESU) agreement with the University of Washington. In order to accurately quantify landscape-level carbon stocks, our analyses accounted for major sources of measurement errors, propagating those errors as we scaled plot-based carbon density estimates up to landscape-level totals. Using Monte Carlo simulation methods, we found that vegetation type mapping error was the largest source of uncertainty, while measurement uncertainties contributed by tree diameter measurements and tree diameter–biomass allometry equations were relatively minor. For some forest types, we found differences in aboveground tree carbon densities between burned and unburned areas. For example, mean carbon density in burned red fir forests was estimated to be ~29% lower versus unburned areas. Alternative measures of fire history, such as time since fire and number of times burned, were poorly related to carbon densities. Within YOSE, we evaluated the stability of landscape carbon pools by quantifying carbon stocks in areas of varying degrees of departure from historic fire return intervals. Of the ~25 Tg of total aboveground tree carbon in YOSE, ~10 Tg is contained within relatively stable areas (the next fire is unlikely to be high severity and stand-replacing), ~10 Tg occurs in areas deemed moderately stable, and the remaining ~5 Tg within relatively unstable areas. We compared our landscape carbon estimates in YOSE to remotely-sensed carbon estimates from the NASA–CASA project and found that the two methods roughly agree. Our analysis and comparisons suggest, however, that fire severity should be integrated into future carbon mapping efforts. We illustrate this with an example using the 2013 Rim Fire, which we estimate burned an area containing over 5 Tg of aboveground tree carbon, but likely left a large fraction of that carbon on the landscape if one accounts for fire severity.

California↗

USGS2006

A field guide to the central, creeping section of the San Andreas fault and the San Andreas Fault Observatory at Depth

This field trip is along the central section of the San Andreas fault and consists of eight stops that illustrate surface evidence of faulting, in general, and features associated with active fault creep, in particular. Fault creep is slippage along a fault that occurs either in association with small-magnitude earthquakes or without any associated large-magnitude earthquakes. Another aspect of the trip is to highlight where there are multiple fault traces along this section of the San Andreas fault zone in order to gain a better understanding of plate-boundary processes. The first stop is along the Calaveras fault, part of the San Andreas fault system, at a location where evidence of active fault creep is abundant and readily accessible. The stops that follow are along the San Andreas fault and at convenient locations to present and discuss rock types juxtaposed across the fault that have been transported tens to hundreds of kilometers by right-lateral motion along the San Andreas fault. Stops 6 and 7 are examples of recent studies of different aspects of the fault: drilling into the fault at the depth of repeating magnitude (M) 2 earthquakes with the San Andreas Fault Observatory at Depth (SAFOD) and the geological, geophysical, and seismological study of M 6 earthquakes near the town of Parkfield. Along with the eight official stops on this field trip are 12 “rolling stops”—sites of geologic interest that add to the understanding of features and processes in the creeping section of the fault. Many of the rolling stops are located where stopping is difficult to dangerous; some of these sites are not appropriate for large vehicles (buses) or groups; some sites are not appropriate for people at all. We include photographs of or from many of these sites to add to the reader's experience without adding too many stops or hazards to the trip. An extensive set of literature is available for those interested in the San Andreas fault or in the creeping section, in particular. For more scientifically oriented overviews of the fault, see Wallace ( 1990 ) and Irwin ( 1990 ); for a more generalized overview with abundant, colorful illustrations, see Collier ( 1999 ). Although the presence of small sections of the San Andreas fault was known before the great 1906 San Francisco earthquake, it was only after that event and subsequent geologic investigations reported in Lawson ( 1908 ) that showed the fault as a long structure, extending all the way from east of Los Angeles into northern California. Prentice ( 1999 ) described the importance of the 1908 “Lawson report” and how it pivotally influenced the understanding of the San Andreas. Hill ( 1981 ) presented a wonderful introduction to the evolution of thought on the San Andreas. Geologic maps and maps of the most recently active fault trace in the creeping section, or large parts of it, include those by Brown ( 1970 ), Dibblee ( 1971 , 1980 ), and Wagner et al. ( 2002 ); detailed geologic maps are discussed at various stops in this guide. Various aspects of the creeping section of the San Andreas fault have been the focus of many geologic field trips in the past few decades. Guidebooks for some of those trips include those by Gribi ( 1963a , 1963b ), Brabb et al. ( 1966 ), Rogers ( 1969 ), Bucknam and Haller ( 1989 ), Harden et al. ( 2001 ), and Stoffer ( 2005 ). The creeping section of the San Andreas fault zone lies between areas that experienced large-displacement surface breakage during great earthquakes in 1857 and 1906 (Fig. 1 inset). Burford and Harsh ( 1980 ) divided the creeping section into three segments: (1) a northwest section where the creep rate increases to the southeast in step-like increments, (2) a central section where the creep rate is relatively constant at a maximum value of ∼30 mm/yr (∼1.2 in/yr), and (3) a southeast section where the creep rate decreases to the southeast (Fig. 2 ). The rate of slip along the creeping section of the fault zone has been measured using creepmeters, alignment arrays, and laser distance-measuring devices. The aperture of measurements over which these measurements are made ranges from 10 m (∼33 ft) (creepmeters) to 100 m (∼330 ft) (alignment arrays) to kilometers and tens of kilometers (laser measuring devices). Creepmeter and alignment-array measurements are here termed “near-fault” measurements; laser measurements over distances of 1–2 km (∼0.6–1.2 mi) are termed “intermediate-scale” measurements; laser measurements over tens of kilometers (miles) are termed “broadscale” measurements. Comparisons among near-fault, intermediate-scale, and broadscale measurements and geologic maps show that the northwest part of the creeping section of the fault is composed of two narrow zones of active deformation, one along the San Andreas fault and one along the Calaveras-Paicines fault, whereas the central and southeast sections are both composed of a single relatively narrow zone of deformation. The southeast section is transitional to a locked zone southeast of Cholame; a locked fault is one that slips only in association with a moderate to large earthquake. Throughout the creeping section of the San Andreas fault zone, broadscale measurements generally indicate more deformation than near-fault and intermediate-scale measurements, which are in reasonably close agreement except at Monarch Peak (Mustang Ridge), near the center of the creeping section and our Stop 5 ((Figs. 1 ) and 2 ). Features that we see on this trip include offset street curbs, closed depressions (sag ponds), fault scarps (steep slopes formed by movement along a fault), a split and displaced tree, offset fence lines, fresh fractures, and offset road lines (Fig. 3 is a sketch showing some of the landforms that represent deformation by an active fault). We also see evidence of long-term maturity of the San Andreas fault, as indicated by fault features and displaced rock types (Fig. 4 ). Finally, we will visit sites of ongoing research into the processes associated with earthquakes and their effects. Discussions include drilling into the San Andreas fault at the SAFOD drill site and the 2004 Parkfield earthquake and its effects and implications.

California↗

USGS2017

Magnitude of flood flows for selected annual exceedance probabilities for streams in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Transportation, determined the magnitude of flood flows at selected annual exceedance probabilities (AEPs) at streamgages in Massachusetts and from these data developed equations for estimating flood flows at ungaged locations in the State. Flood magnitudes were determined for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent AEPs at 220 streamgages, 125 of which are in Massachusetts and 95 are in the adjacent States of Connecticut, New Hampshire, New York, Rhode Island, and Vermont. AEP flood flows were computed for streamgages using the expected moments algorithm weighted with a recently computed regional skewness coefficient for New England. Regional regression equations were developed to estimate the magnitude of floods for selected AEP flows at ungaged sites from 199 selected streamgages and for 60 potential explanatory basin characteristics. AEP flows for 21 of the 125 streamgages in Massachusetts were not used in the final regional regression analysis, primarily because of regulation or redundancy. The final regression equations used generalized least squares methods to account for streamgage record length and correlation. Drainage area, mean basin elevation, and basin storage explained 86 to 93 percent of the variance in flood magnitude from the 50- to 0.2-percent AEPs, respectively. The estimates of AEP flows at streamgages can be improved by using a weighted estimate that is based on the magnitude of the flood and associated uncertainty from the at-site analysis and the regional regression equations. Weighting procedures for estimating AEP flows at an ungaged site on a gaged stream also are provided that improve estimates of flood flows at the ungaged site when hydrologic characteristics do not abruptly change. Urbanization expressed as the percentage of imperviousness provided some explanatory power in the regional regression; however, it was not statistically significant at the 95-percent confidence level for any of the AEPs examined. The effect of urbanization on flood flows indicates a complex interaction with other basin characteristics. Another complicating factor is the assumption of stationarity, that is, the assumption that annual peak flows exhibit no significant trend over time. The results of the analysis show that stationarity does not prevail at all of the streamgages. About 27 percent of streamgages in Massachusetts and about 42 percent of streamgages in adjacent States with 20 or more years of systematic record used in the study show a significant positive trend at the 95-percent confidence level. The remaining streamgages had both positive and negative trends, but the trends were not statistically significant. Trends were shown to vary over time. In particular, during the past decade (2004–2013), peak flows were persistently above normal, which may give the impression of positive trends. Only continued monitoring will provide the information needed to determine whether recent increases in annual peak flows are a normal oscillation or a true trend. The analysis used 37 years of additional data obtained since the last comprehensive study of flood flows in Massachusetts. In addition, new methods for computing flood flows at streamgages and regionalization improved estimates of flood magnitudes at gaged and ungaged locations and better defined the uncertainty of the estimates of AEP floods.

Massachusetts↗

USGS2016

Methods for estimating selected spring and fall low-flow frequency statistics for ungaged stream sites in Iowa, based on data through June 2014

A statewide study was led to develop regression equations for estimating three selected spring and three selected fall low-flow frequency statistics for ungaged stream sites in Iowa. The estimation equations developed for the six low-flow frequency statistics include spring (April through June) 1-, 7-, and 30-day mean low flows for a recurrence interval of 10 years and fall (October through December) 1-, 7-, and 30-day mean low flows for a recurrence interval of 10 years. Estimates of the three selected spring statistics are provided for 241 U.S. Geological Survey continuous-record streamgages, and estimates of the three selected fall statistics are provided for 238 of these streamgages, using data through June 2014. Because only 9 years of fall streamflow record were available, three streamgages included in the development of the spring regression equations were not included in the development of the fall regression equations. Because of regulation, diversion, or urbanization, 30 of the 241 streamgages were not included in the development of the regression equations. The study area includes Iowa and adjacent areas within 50 miles of the Iowa border. Because trend analyses indicated statistically significant positive trends when considering the period of record for most of the streamgages, the longest, most recent period of record without a significant trend was determined for each streamgage for use in the study. Geographic information system software was used to measure 63 selected basin characteristics for each of the 211streamgages used to develop the regional regression equations. The study area was divided into three low-flow regions that were defined in a previous study for the development of regional regression equations. Because several streamgages included in the development of regional regression equations have estimates of zero flow calculated from observed streamflow for selected spring and fall low-flow frequency statistics, the final equations for the three low-flow regions were developed using two types of regression analyses—left-censored and generalized-least-squares regression analyses. A total of 211 streamgages were included in the development of nine spring regression equations—three equations for each of the three low-flow regions. A total of 208 streamgages were included in the development of nine fall regression equations—three equations for each of the three low-flow regions. A censoring threshold was used to develop 15 left-censored regression equations to estimate the three fall low-flow frequency statistics for each of the three low-flow regions and to estimate the three spring low-flow frequency statistics for the southern and northwest regions. For the northeast region, generalized-least-squares regression was used to develop three equations to estimate the three spring low-flow frequency statistics. For the northeast region, average standard errors of prediction range from 32.4 to 48.4 percent for the spring equations and average standard errors of estimate range from 56.4 to 73.8 percent for the fall equations. For the northwest region, average standard errors of estimate range from 58.9 to 62.1 percent for the spring equations and from 83.2 to 109.4 percent for the fall equations. For the southern region, average standard errors of estimate range from 43.2 to 64.0 percent for the spring equations and from 78.1 to 78.7 percent for the fall equations. The regression equations are applicable only to stream sites in Iowa with low flows not substantially affected by regulation, diversion, or urbanization and with basin characteristics within the range of those used to develop the equations. The regression equations will be implemented within the U.S. Geological Survey StreamStats Web-based geographic information system application. StreamStats allows users to click on any ungaged stream site and compute estimates of the six selected spring and fall low-flow statistics; in addition, 90-percent prediction intervals and the measured basin characteristics for the ungaged site are provided. StreamStats also allows users to click on any Iowa streamgage to obtain computed estimates for the six selected spring and fall low-flow statistics.

Iowa↗

USGS2017

Genomics of Arctic cod

The Arctic cod (Boreogadus saida) is an abundant marine fish that plays a vital role in the marine food web. To better understand the population genetic structure and the role of natural selection acting on the maternally-inherited mitochondrial genome (mitogenome), a molecule often associated with adaptations to temperature, we analyzed genetic data collected from 11 biparentally-inherited nuclear microsatellite DNA loci and nucleotide sequence data from from the mitochondrial DNA (mtDNA) cytochrome b (cytb) gene and, for a subset of individuals, the entire mitogenome. In addition, due to potential of species misidentification with morphologically similar Polar cod (Arctogadus glacialis), we used ddRAD-Seq data to determine the level of divergence between species and identify species-specific markers. Based on the findings presented here, Arctic cod across the Pacific Arctic (Bering, Chukchi, and Beaufort Seas) comprise a single panmictic population with high genetic diversity compared to other gadids. High genetic diversity was indicated across all 13 protein-coding genes in the mitogenome. In addition, we found moderate levels of genetic diversity in the nuclear microsatellite loci, with highest diversity found in the Chukchi Sea. Our analyses of markers from both marker classes (nuclear microsatellite fragment data and mtDNA cytb sequence data) failed to uncover a signal of microgeographic genetic structure within Arctic cod across the three regions, within the Alaskan Beaufort Sea, or between near-shore or offshore habitats. Further, data from a subset of mitogenomes revealed no genetic differentiation between Bering, Chukchi, and Beaufort seas populations for Arctic cod, Saffron cod (Eleginus gracilis), or Walleye pollock (Gadus chalcogrammus). However, we uncovered significant differences in the distribution of microsatellite alleles between the southern Chukchi and central and eastern Beaufort Sea samples of Arctic cod. Finally, using ddRAD-Seq data, we identified species-specific markers and in conjunction with mitogenome data, identified an Arctic cod x Polar cod hybrid in western Canadian Beaufort Sea. Overall, the lack of genetic structure among Arctic cod within the Bering, Chukchi and Beaufort seas of Alaska is concordant with the absence of geographic barriers to dispersal and typical among marine fishes. Arctic cod may exhibit a genetic pattern of isolation-by-distance, whereby populations in closer geographic proximity are more genetically similar than more distant populations. As this signal is only found between our two fartherest localities, data from populations elsewhere in the species’ global range are needed to determine if this is a general characteristic. Further, tests for selection suggested a limited role for natural selection acting on the mitochondrial genome of Arctic cod, but do not exclude the possibility of selection on genes involved in nuclear-mitogenome interactions. Unlike previous genetic assessment of Arctic cod sampled from the Chukchi Sea, the high levels of genetic diversity found in Arctic cod assayed in this study, across regions, suggests that the species in the Beaufort and Chukchi seas does not suffer from low levels of genetic variation, at least at neutral genetic markers. The large census size of Arctic cod may allow this species to retain high levels of genetic diversity. In addition, we discovered the presence of hybridization between Arctic and Polar cod (although low in frequency). Hybridization is expected to occur when environmental changes modify species distributions that result in contact between species that were previously separated. In such cases, hybridization may be an evolutionary mechanism that promotes an increase in genetic diversity that may provide species occupying changing environments with locally-adapted genotypes and, therefore, phenotypes. Natural selection can only act on the standing genetic variation present within a population. Therefore, given its higher levels of genetic diversity in combination with a large population size, Arctic cod may be resilient to current and future environmental change, as high genetic diversity is expected to increase opportunities for positive selection to act on genetic variants beneficial in different environments, regardless of the source of that genetic variation.

OCS Study↗

USGS2010

Death Valley regional groundwater flow system, Nevada and California: Hydrogeologic framework and transient groundwater flow model

A numerical three-dimensional (3D) transient groundwater flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the groundwater flow system and previous less extensive groundwater flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect groundwater flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley regional groundwater flow system (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the groundwater flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural groundwater discharge occurring through evapotranspiration (ET) and spring flow; the history of groundwater pumping from 1913 through 1998; groundwater recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were provided by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional groundwater flow system. Groundwater flow in the Death Valley region is composed of several interconnected, complex groundwater flow systems. Groundwater flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional groundwater flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and ET are the dominant natural groundwater discharge processes. Groundwater also is withdrawn for agricultural, commercial, and domestic uses. Groundwater flow in the DVRFS was simulated using MODFLOW-2000, the U.S. Geological Survey 3D finitedifference modular groundwater flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 meters (m) on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins simulated by constant-head boundaries. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient model generally fit observed heads reasonably well (residuals with absolute values less than 10 meters) with two exceptions: in most areas of nearly flat hydraulic gradient the fit is considered moderate (residuals with absolute values of 10 to 20 meters), and in areas of steep hydraulic gradient along the Eleana Range and western part of Yucca Flat, southern part of the Owlshead Mountains, southern part of the Bullfrog Hills, and the north-northwestern part of the model domain (residuals with absolute values greater than 20 meters). Groundwater discharge residuals are fairly random, with as many areas where simulated flows are less than observed flows as areas where simulated flows are greater. The highest unweighted groundwater discharge residuals occur at Death Valley, Sarcobatus Flat (northeastern area), Tecopa, and early observations at Manse Spring in Pahrump Valley. High weighted-discharge residuals were computed in Indian Springs Valley and parts of Death Valley. Most of these inaccuracies in head and discharge can be attributed to insufficient representation of the hydrogeology in the HFM and(or) discharge estimates, misrepresentation of water levels, and(or) model error associated with grid-cell size. The model represents the large and complex groundwater flow system of the Death Valley region at a greater degree of refinement and accuracy than has been possible previously. The representation of detail provided by the 3D digital hydrogeologic framework model and the numerical groundwater flow model enabled greater spatial accuracy in every model parameter. The lithostratigraphy and structural effects of the hydrogeologic framework; recharge estimates from simulated net infiltration; discharge estimates from ET, spring flow, and pumping; and boundary inflow and outflow estimates all were reevaluated, some additional data were collected, and accuracy was improved. Uncertainty in the results of the flow model simulations can be reduced by improving on the quality, interpretation, and representation of the water-level and discharge observations used to calibrate the model and improving on the representation of the HGU geometries, the spatial variability of HGU material properties, the flow model physical framework, and the hydrologic conditions.

California, Nevada↗

USGS2016

Late Holocene volcanism at Medicine Lake Volcano, northern California Cascades

Late Holocene volcanism at Medicine Lake volcano in the southern Cascades arc exhibited widespread and compositionally diverse magmatism ranging from basalt to rhyolite. Nine well-characterized eruptions have taken place at this very large rear-arc volcano since 5,200 years ago, an eruptive frequency greater than nearly all other Cascade volcanoes. The lavas are widely distributed, scattered over an area of ~300 km 2 across the >2,000-km 2 volcano. The eruptions are radiocarbon dated and the ages are also constrained by paleomagnetic data that provide strong evidence that the volcanic activity occurred in three distinct episodes at ~1 ka, ~3 ka, and ~5 ka. The ~1-ka final episode produced a variety of compositions including west- and north-flank mafic flows interspersed in time with fissure rhyolites erupted tangential to the volcano’s central caldera, including the youngest and most spectacular lava flow at the volcano, the ~950-yr-old compositionally zoned Glass Mountain flow. At ~3 ka, a north-flank basalt eruption was followed by an andesite eruption 27 km farther south that contains quenched basalt inclusions. The ~5-ka episode produced two caldera-focused dacitic eruptions. Quenched magmatic inclusions record evidence of intrusions that did not independently reach the surface. The inclusions are present in five andesitic, dacitic, and rhyolitic host lavas, and were erupted in each of the three episodes. Compositional and mineralogic evidence from mafic lavas and inclusions indicate that both tholeiitic (dry) and calcalkaline (wet) parental magmas were present. Petrologic evidence records the operation of complex, multi-stage processes including fractional crystallization, crustal assimilation, and magma mixing. Experimental evidence suggests that magmas were stored at 3 to 6 km depth prior to eruption, and that both wet and dry parental magmas were involved in generating the more silicic magmas. The broad distribution of eruptive events and the relative accessibility and good exposure of lavas, combined with physical and petrologic evidence for multiple and varied mafic inputs, has created an unusual opportunity to understand the workings of this large magmatic system. A combined total of more than 25 intrusive and extrusive events are indicated for late Holocene time. Plutonic inclusions, some with ages as young as Holocene, were also brought to the surface in five of the eruptions. All eruptions took place along northwest- to northeast-trending alignments of vents, reflecting the overall east-west extensional tectonic environment. The interaction of tectonism and volcanism is a dominant influence at this subduction-related volcano, located where the west edge of the extensional Basin and Range Province impinges on the Cascades arc. Ongoing subsidence focused at the central caldera has been documented along with geophysical evidence for a small magma body. This evidence, combined with the frequency of eruptive and intrusive activity in late Holocene time, an active geothermal system, and intermittent long-period seismic events indicate that the volcano is likely to erupt again.

California↗

USGS2001

Large carbon isotope fractionation associated with oxidation of methyl halides by methylotrophic bacteria

The largest biological fractionations of stable carbon isotopes observed in nature occur during production of methane by methanogenic archaea. These fractionations result in substantial (as much as ≈70‰) shifts in δ 13 C relative to the initial substrate. We now report that a stable carbon isotopic fractionation of comparable magnitude (up to 70‰) occurs during oxidation of methyl halides by methylotrophic bacteria. We have demonstrated biological fractionation with whole cells of three methylotrophs (strain IMB-1, strain CC495, and strain MB2) and, to a lesser extent, with the purified cobalamin-dependent methyltransferase enzyme obtained from strain CC495. Thus, the genetic similarities recently reported between methylotrophs, and methanogens with respect to their pathways for C 1 -unit metabolism are also reflected in the carbon isotopic fractionations achieved by these organisms. We found that only part of the observed fractionation of carbon isotopes could be accounted for by the activity of the corrinoid methyltransferase enzyme, suggesting fractionation by enzymes further along the degradation pathway. These observations are of potential biogeochemical significance in the application of stable carbon isotope ratios to constrain the tropospheric budgets for the ozone-depleting halocarbons, methyl bromide and methyl chloride. Methyl bromide (MeBr) and methyl chloride (MeCl) are, respectively, the most abundant volatile bromo- and chlorocarbons in the troposphere and are major contributors to stratospheric ozone destruction (1). Both compounds have natural and human-influenced sources and a predominant sink by reaction with OH in the troposphere (2–4). MeBr also has a bacterial soil sink (5) that represents about 20% of the estimated total removal from the troposphere, and it is likely that a soil sink of similar magnitude exists for MeCl (6). Hence, if an isotopic fractionation is associated with the soil sink, it will influence the isotopic compositions of MeBr and MeCl in the lower atmosphere (7). The δ 13 C value of industrially produced MeBr ranges between −43.5‰ and −66.4‰ (7), but δ 13 C values of tropospheric MeBr and natural sources are not yet known. The δ 13 C of atmospheric MeCl has been measured from −22‰ to −45‰ (8, 9). If carbon isotope ratios are to be used to constrain the budgets of these methyl halides, it is essential to determine the extent of carbon isotope fractionation that occurs during biological degradation of these compounds. Methylotrophic bacteria use C 1 compounds, which are simple organic molecules that contain no carbon–carbon bonds. Strains IMB-1, CC495, and MB2 are as-yet-unnamed facultative methylotrophs isolated from agricultural soil, woodland leaf litter, and coastal seawater, respectively (10–13), environments where methyl halides are produced. They are members of the α subgroup of the Proteobacteria . On the basis of 16S rRNA gene sequences, strains IMB-1 and CC495 show some phylogenetic alignment with the genus Rhizobium (10, 11) and are very closely related to the new genus Pseudoaminobacte r (I. McDonald, personal communication). Strain MB2 aligns within the Ruegeria clade [J. K. Schaefer, K. D. Goodwin, I. R. McDonald, J. C. Murrell and R.S.O., unpublished work]. All of these aerobic bacteria are methylotrophs in that they can grow by using MeBr or MeCl as their sole carbon source, but they do not metabolize methane. They oxidize MeBr, MeCl, and methyl iodide (MeI) to CO 2 . Soil bacteria are known to consume MeBr at the ambient tropospheric mixing ratio of around 10 parts per trillion by volume (5). Preliminary experiments with strain IMB-1 indicate that it can oxidize MeBr at these mixing ratios** and is therefore likely to be characteristic of bacteria associated with MeBr uptake by soils. We examined δ 13 C of MeCl, MeBr, and MeI during oxidation by whole-cell suspensions of IMB-1 and CC495 and also the change in δ 13 C values of the three methyl halides during oxidation by the marine strain MB2. In addition, we measured the fractionation of carbon isotopes during formation of methane thiol (MeSH) from MeCl by the purified cobalamin-dependent enzyme, halomethane:bisulphide/halide ion methyltransferase (11) from CC495, to determine whether this initial step in MeCl degradation could account for the observed fractionation by whole cells. Finally, we determined the fractionation associated with the degradation of MeBr during field studies with agricultural soil by monitoring MeBr concentration and δ 13 C of MeBr in the headspace of flux chambers under fumigation conditions.

Proceedings of the National Academy of Sciences of↗

USGS1980

Handbook of acute toxicity of chemicals to fish and aquatic invertebrates : summaries of toxicity tests conducted at Columbia National Fisheries Research Laboratory, 1965-78

Acute toxicity is a major subject of research at Columbia National Fisheries Research Laboratory for evaluating the impact of toxic chemicals on fishery resources. The Laboratory has played a leading role in developing research technology for toxicity testing and data interpretation. In 1965-78, more than 400 chemicals were tested against a variety of invertebrates and fish species representative of both cold- and warm-water climates. The use of acute toxicity tests for assessing the potential hazard of chemical contaminants to aquatic organisms is well documented (Boyd 1957; Henderson et al. 1960; Sanders and Cope 1966; Macek and McAllister 1970). Static acute toxicity tests provide rapid and (within limits) reproducible concentration-response curves for estimating toxic effects of chemicals on aquatic organisms. These tests provide a database for determining relative toxicity of a large number of chemicals to a variety of species and for estimating acute effects of chemical spills on natural aquatic systems; they also assist in determining priority and design of additional toxicity studies. Acute toxicity tests usually provide estimates of the exposure concentration causing 50% mortality (LC50) to test organisms during a specified period of time. For certain invertebrates, the effective concentration is based on immobilization, or some other identifiable endpoint, rather than on lethality. The application of the LC50 has gained acceptance among toxicologists and is generally the most highly rated test for assessing potential adverse effects of chemical contaminants to aquatic life (Brungs and Mount 1978; American Institute for Biological Sciences 1978 a ). The literature contains numerous papers dealing with the acute toxicity of chemicals to freshwater organisms. However, there is a tremendous need for a concise compendium of toxicity data covering a large variety of chemicals and test species. This Handbook is a compilation of a large volume of acute toxicity data from the Columbia Laboratory and its field laboratories. It presents definitive acute toxicity data on 271 chemicals tested against a variety of freshwater invertebrates and fishes. The chemicals represent all major groups of pesticides, as well as numerous industrial chemicals. This compilation should serve as a useful database for the many agencies and organizations dealing with research and management programs concerned with the impact of chemicals on aquatic resources. The Columbia Laboratory has played a major role in developing currently used standard methodology for static acute toxicity testing. The use of standardized methodology greatly reduces variation in results. The data presented here have been carefully scrutinized to eliminate tests that failed to follow acceptable procedures. Handling of test organisms and procedures for static toxicity tests followed those described by Lennon and Walker (1964) and Macek and McAllister (1970), and conform well with those recommended by Brauhn and Schoettger (1975) and the Committee on Methods for Toxicity Tests with Aquatic Organisms (1975). The species of fish and invertebrates that were tested are listed in phylogenetic order in Tables 1 and 2. Fish were obtained from Federal and State hatcheries as either eggs or fry. Original stocks of invertebrates were collected and cultured from wild populations with no known source of contamination; these populations were replenished regularly. The invertebrates were cultured in the Laboratory by methods similar to those described by Sanders and Cope (1966). Test chemicals usually consisted of technical or analytical grade samples of known purity. Formulations of the chemicals were also tested when available. When purity of test chemicals was known, all calculated concentrations were based on percent active ingredients. Stock solutions were prepared immediately before each test, with commercial grade acetone as the carrier solvent. Occasionally, ethanol or dimethyl-formamide was substituted. Solvent concentrations did not exceed 0.5 mL/L in final dilution water. Test water (dilution water) was reconstituted from deionized water of at least 10 6 ohms resistivity by the addition of appropriate reagent grade chemicals (Marking 1969). Water was buffered to maintain a pH of 7.2 to 7.5, an alkalinity of 30 to 35 mg/L, and a hardness of 40 to 50 mg/L as CaCO 3 . Test water was mixed thoroughly and aerated before transfer into test chambers. Fish were acclimated to dilution water by gradually changing the water in acclimated tanks from 100% well water to 100% reconstituted water over a 1- to 3-day period at the desired testing temperature. Invertebrates were acclimated from well water to dilution water over a 4- to 6-h period. Toxicity tests were conducted under static conditions without aeration, and the organisms were not fed during acclimation or testing. Temperature of test solutions was maintained within ± 1°C of that required for a given test. Toxicity tests with fish were conducted in 18.9-liter (5-gal) wide-mouthed jars containing 15 liters of test solution. Fingerling fish weighing 0.2 to 1.5 g were tested at each concentration. Caution was taken not to exceed 0.8 g of test organisms per liter of solution. Duplicate test chambers were used to accommodate larger fish. Test chambers varied in size for invertebrates, depending on the species used; volume of test solution ranged from 0.25 to 4 liters. At least 10 organisms were exposed to each concentration for all definitive tests. At least six concentrations were used per toxicity test. The tests began upon initial exposure to the toxicant and continued for 96 h. Immobilization tests with invertebrates were conducted for only 48 h. The number of dead or affected organisms in each test chamber were recorded and the dead organisms were removed every 24 h; general observations on the condition of test organisms were also recorded at these times. Toxicity data were analyzed by a statistical method described by Litchfield and Wilcoxon (1949) to determine LC50 (theoretical estimate of the concentration lethal to 50% of the test animals) and 95% confidence intervals. This method is recommended by the American Public Health Association (1971) and by Sprague (1969) for determining median lethal concentrations. The procedure is easily modified for computing a single LC50 when replicate tests are performed.

Resource Publication↗

USGS2011

Two-dimensional streamflow simulations of the Jordan River, Midvale and West Jordan, Utah

The Jordan River in Midvale and West Jordan, Utah, flows adjacent to two U.S. Environmental Protection Agency Superfund sites: Midvale Slag and Sharon Steel. At both sites, geotechnical caps extend to the east bank of the river. The final remediation tasks for these sites included the replacement of a historic sheet-pile dam and the stabilization of the river banks adjacent to the Superfund sites. To assist with these tasks, two hydraulic modeling codes contained in the U.S. Geological Survey (USGS) Multi-Dimensional Surface-Water Modeling System (MD_SWMS), System for Transport and River Modeling (SToRM) and Flow and Sediment Transport and Morphological Evolution of Channels (FaSTMECH), were used to provide predicted water-surface elevations, velocities, and boundary shear-stress values throughout the study reach of the Jordan River. A SToRM model of a 0.7 mile subreach containing the sheet-pile dam was used to compare water-surface elevations and velocities associated with the sheet-pile dam and a proposed replacement structure. Maps showing water-surface elevation and velocity differences computed from simulations of the historic sheet-pile dam and the proposed replacement structure topographies for streamflows of 500 and 1,000 cubic feet per second (ft 3 /s) were created. These difference maps indicated that the velocities associated with the proposed replacement structure topographies were less than or equal to those associated with the historic sheet-pile dam. Similarly, water-surface elevations associated with the proposed replacement structure topographies were all either greater than or equal to water-surface elevations associated with the sheet-pile dam. A FaSTMECH model was developed for the 2.5-mile study reach to aid engineers in bank stabilization designs. Predicted water-surface elevations, velocities and shear-stress values were mapped on an aerial photograph of the study reach to place these parameters in a spatial context. Profile plots of predicted cross-stream average water-surface elevations and cross-stream maximum and average velocities showed how these parameters change along the study reach for two simulated discharges of 1,040 ft 3 /s and 2,790 ft 3 /s. The profile plots for the simulated streamflow of 1,040 ft 3 /s show that the highest velocities are associated with the constructed sheet-pile replacement structure. Results for the simulated streamflow of 2,790 ft 3 /s indicate that the geometry of the 7800 South Bridge causes more backwater and higher velocities than the constructed sheet-pile replacement structure.

Utah↗

USGS2020

Reconnaissance map of the Cenozoic geology in the Carlin basin area, Elko and Eureka counties, Nevada

The middle Miocene Carlin sedimentary basin encompasses a large area between the Adobe Range to the east, the Piñon Range to the south, the southern Independence Mountains and Marys Mountain to the west, and Swales Mountain to the north. The town of Carlin is in the southern part of the basin. The geologic map includes detailed to more reconnaissance mapping of Cenozoic units in the main part of the basin, including different facies of the middle Miocene Humboldt Formation. The mapping was part of a broader study of the Miocene and younger paleogeographic evolution of the region. Earlier work obtained numerous 40 Ar/ 39 Ar and tephra correlation dates on sedimentary and volcanic units in the basin. The basin connected to the east into the Elko sedimentary basin and to the northwest and northeast into similar smaller basins between present-day mountain ranges. Early sediments, largely fluvial, began to accumulate in the lowlands between the surrounding ranges at about 16.5 Ma. The sediments were derived from Paleozoic sedimentary and middle Tertiary volcanic units in the nearby highlands, and flow patterns indicate a general southward flow towards present-day Pine Valley. The eruption of the Palisade Canyon–Marys Mountain rhyolite flows at the southwest end of the basin at 15.3 Ma blocked the southward flow, and a lake began to form in the basin. As the lake grew in extent, sedimentary units around the fringes of the lake included a mixture of inflowing fluvial sediments mixed with the pyroclastic-fall, ash-rich sediments deposited in the lake. The lake margin expanded, and stratigraphic sections record the progressive transition from fluvial to mixed fluvial and lacustrine, and finally to entirely lacustrine. The volcanic rock dam was breached at about 14.7 Ma, the lake drained, and fluvial sediments blanketed the entire basin for an unknown period of time after that. Sedimentation progressively buried existing highlands and bridged gaps between adjacent basins. For example, the Carlin and Elko basins connected across the southern Adobe Range. Normal faulting produced numerous, mostly north- to north-northeast-striking faults that cut the sedimentary units and surrounding highlands largely after sedimentation ceased. The largest fault formed in the eastern third of the basin and tilted all of the sedimentary units in the western two-thirds of the basin, as well as the eastern part of Marys Mountain, to the east. Some offset took place during sedimentation. Many other normal faults of smaller extent and offset cut the sedimentary units. The integration of streams draining the Elko and Carlin basins began after about 9.8 Ma. The streams, which together comprised the early stages of the Humboldt River, flowed regionally southwestward beyond the Carlin basin. As many as thirteen, downward-stepping strath terraces in the Carlin basin record the progressive downcutting into and removal of the middle Miocene sediments. Gravel deposits form a thin veneer on some of the higher terraces. Clasts in those gravel deposits, as well as the overall terrace pattern, indicate southward drainage towards the Humboldt River. The erosion gradually re-exposed the flanks of the surrounding highlands. A brief pause in downcutting allowed the formation of a small lake in the Hemphillian (late Miocene), represented by lacustrine units northwest of Carlin. The sedimentary rocks of the Carlin basin conceal a large segment of the world-class, late Eocene Carlin gold trend, which extends from the southern Independence Mountains south into the Piñon Range. Sedimentation largely buried the Gold Quarry and Mike gold deposits in the northwestern part of the basin. Later faulting and erosion re-exposed the Gold Quarry deposit, but the Mike deposit remains buried. The basin’s sedimentary units conceal potential Paleozoic host rocks, and the sedimentary facies and post-sedimentation faults shown on the map may help guide interpretations of geophysical and other exploration data in the Carlin basin. The current map publication was supported by the USGS National Cooperative Geologic Mapping Program under STATEMAP award number G19AC00383.

Nevada↗

USGS1979

Seismic refraction profile, Kingdom of Saudi Arabia: field operations, instrumentation, and initial results

In February 1978 a seismic deep-refraction profile was recorded by the USGS along a 1000-km line across the Arabian Shield in western Saudi Arabia. The line begins in Paleozoic and Mesozoic cover rocks near Riyadh on the Arabian Platform, leads southwesterly across three major Precambrian tectonic provinces, traverses Cenozoic rocks of the coastal plain near Jizan (Tihamat Asir), and terminates at the outer edge of the Farasan Bank in the southern Red Sea. More than 500 surveyed recording sites were occupied, including 19 in the Farasan Islands. Six shot points were used--five on land, with charges placed mostly below water table in drill holes, and one at sea, with charges placed on the sea floor and fired from a ship. The total charge consumed was slightly in excess of 61 metric tons in 21 discrete firings. Seismic energy was recorded by means of a set of 100 newly developed portable seismic stations. Each station consists of a standard 2-Hz vertical geophone coupled to a self-contained analog recording instrument equipped with a magnetic-tape cassette. The stations were deployed in groups of 20 by five observer teams, each generally consisting of two scientist-technicians and a surveyor-guide. On the day prior to deployment, the instruments were calibrated and programmed for automatic operation by means of a specially designed device called a hand-held tester. At each of ten pre-selected recording time windows on a designated firing day, the instruments were programmed to turn on, stabilize, record internal calibration signals, record the seismic signals at three levels of amplification, and then deactivate. After the final window in the firing sequence, all instruments were retrieved and their data tapes removed for processing. A specially designed, field tape- dubbing system was utilized at shot point camps to organize and edit data recorded on the cassette tapes. The main functions of this system are to concatenate all data from each shot on any given day onto a single shot tape, and to provide hard copy for monitoring recorder performance so that any problems can be corrected prior to the next deployment. Composite digital record sections were produced from the dubbed tapes for each shot point by a portable processing and plotting system. The heart of this system is a DEC PDP 11VO3 computer, which controls a cassette playback unit identical to those used in the recorders and dubbers, a set of discriminators, a time-code translator, a digitizer, and a digital plotter. The system was used to maintain various informational data sets and to produce tabulations and listings of various sorts during the field operations, in addition to its main task of producing digital record sections. Two master clocks, both set to time signals broadcast by the British Broadcasting Corporation, provided absolute time for the recording operations. One was located on the ship and the other was stationed at a base camp on the mainland. The land-based master clock was used to set three additional master clocks located at the other active shot points a few days in advance of each firing, and these clocks were then used to set the internal clocks in the portable seismic stations via the hand-held tester. A master clock signal was also linked to the firing system at each shot point for determination of the absolute shot instant. It is possible to construct a generalized crustal model from examination of the six shot point composite record sections obtained in the field. Such a model rests upon a number of simplifying assumptions and will almost certainly be modified at a later stage of interpretation. The main assumptions are that the crust consists of two homogeneous isotropic layers having no velocity inversion,, that the Mohorovicic discontinuity is sharp, and that effects of surface inhomogeneities and elevation changes can be ignored. The main characteristics of the tentative model are the following: (1) The thickness of th

Open-File Report↗

USGS2016

A spatially explicit suspended-sediment load model for western Oregon

We calibrated the watershed model SPARROW (Spatially Referenced Regressions on Watershed attributes) to give estimates of suspended-sediment loads for western Oregon and parts of northwestern California. Estimates of suspended-sediment loads were derived from a nonlinear least squares regression that related explanatory variables representing landscape and transport conditions to measured suspended-sediment loads at 68 measurement stations. The model gives estimates of model coefficients and their uncertainty within a spatial framework defined by the National Hydrography Dataset Plus hydrologic network. The resulting model explained 64 percent of the variability in suspended-sediment yield and had a root mean squared error value of 0.737. The predictor variables selected for the final model were (1) generalized lithologic province, (2) mean annual precipitation, and (3) burned area (by recent wildfire). Other landscape characteristics also were considered, but they were not significant predictors of sediment transport, were strongly correlated with another predictor variable, or were not as significant as the predictors selected for the final model. The northern Oregon coastal drainages had the highest predicted suspended sediment yields (median yield 475 kilograms per hectare per year) and the Klamath River Basin had the lowest (median yield 53 kilograms per hectare per year). Quaternary deposits were, on average, the largest contributor to incremental suspended-sediment yield even though this lithologic province only makes up 17 percent of the modeling domain. Coast Range sedimentary rocks and Coast Range volcanic rocks had high suspended-sediment yields whereas, in addition to the Klamath terrane, the Western Cascade and High Cascade lithologic provinces had low suspended-sediment yields. Precipitation and the area affected by recent wildfire both positively correlated with suspended-sediment load. Suspended-sediment transport rates predicted by this SPARROW model are less than historical (1956–73) and long‑term (thousands of years) geological rates. This difference likely results, in part, from biases in the data underlying the SPARROW model, probably resulting in predicted suspended-sediment estimates that underestimate actual transport rates. However, the differences also likely owe to natural and human-caused variation in suspended-sediment yields as they respond to changes in climate, vegetation, fire frequency, and land use. In particular, decreases in mean annual suspended-sediment yields within the Umpqua River Basin since 1956–73 may owe to less intense forest harvest, passage of the Oregon Forest Practices Act of 1971, and increased emphasis in habitat protection in recent decades. Such sensitivity may have implications for the spatial and temporal distributions of aquatic and riparian habitats. Knowledge of the regionally important patterns and factors in suspended-sediment sources and transport could support broad-scale, water-quality management objectives and priorities. Because of biases and limitations of this model, however, these results are most applicable for general comparisons and for broad areas such as large watersheds. For example, despite having similar area, precipitation, and land-use, the Umpqua River Basin generates 68 percent more suspended sediment than the Rogue River Basin, chiefly because of the large area of Coast Range sedimentary province in the Umpqua River Basin. By contrast, the Rogue River Basin contains a much larger area of Klamath terrane rocks, which produce significantly less suspended load, although recent fire disturbance (in 2002) has apparently elevated suspended sediment yields in the tributary Illinois River watershed. Fine-scaled analysis, however, will require more intensive, locally focused measurements.

Scientific Investigations Report↗

USGS2014

Porphyry copper assessment of East and Southeast Asia: Philippines, Taiwan (Republic of China), Republic of Korea (South Korea), and Japan

The U.S. Geological Survey collaborated with member countries of the Coordinating Committee for Geoscience Programmes in East and Southeast Asia (CCOP) on an assessment of the porphyry copper resources of East and Southeast Asia as part of a global mineral resource assessment. The assessment covers the Philippines in Southeast Asia, and the Republic of Korea (South Korea), Taiwan (Province of China), and Japan in East Asia. The Philippines host world class porphyry copper deposits, such as the Tampakan and Atlas deposits. No porphyry copper deposits have been discovered in the Republic of Korea (South Korea), Taiwan (Province of China), or Japan. Thirteen geographic areas were delineated as tracts that are permissive for porphyry copper deposits in the assessed area. Individual tracts range from about 3,000 to 100,000 square kilometers in area. Permissive tracts are delineated on the basis of mapped distributions of igneous rocks of specific age ranges that define subduction-related magmatic arcs or magmatic belts that might contain porphyry copper deposits. Most of these magmatic arcs are subduction related, although some porphyry deposits and prospects are present in back-arc or poorly understood tectonic settings. Maps at various scales were used in the compilation; however, the final tract boundaries are intended for use at a scale of 1:1,000,000. Numbers of undiscovered deposits were estimated at different levels of confidence for 10 permissive tracts in the Philippines including one area that extends to eastern Taiwan (Republic of China); permissive tracts in South Korea and Japan are discussed qualitatively. Estimates of numbers of undiscovered deposits were combined with grade and tonnage models using Monte Carlo simulation to estimate amounts of undiscovered resources. Grades and tonnages of known porphyry copper deposits in the study area were compared with global grade and tonnage models to determine the appropriate model for simulation of undiscovered resources. Most of the known deposits are best described as copper-gold subtypes of porphyry copper deposits. For some permissive tracts, a general porphyry copper-gold-molybdenum model was used. Thirty-eight porphyry copper deposits are known in the Philippines; the mean number of undiscovered deposits was estimated to be 28. Mean (arithmetic) resources that could be associated with the undiscovered deposits are 90 million metric tons of copper and 5,800 metric tons of gold, as well as byproduct molybdenum and silver. Additional resources that could be discovered in extensions to known deposits were not evaluated. Assessment results, presented in tables and graphs, indicate expected amounts of total contained metal and mineralized rock in undiscovered deposits at different quantile levels, as well as the arithmetic mean for each tract. The Philippines have a long history of porphyry exploration cycles and mine development, interrupted at times by political and social unrest, environmental concerns, and natural disasters. Changes in mining laws within the region and the recent high price of gold on the world market have prompted renewed interest in porphyry copper deposits in the region. South Korea and Japan have been thoroughly explored for many types of mineral deposits. Available data suggest that the permissive rocks in South Korea typically are too deeply eroded to preserve porphyry copper deposits. Porphyry copper systems may be present in Japan, but are likely to lie at depths greater than the 1 kilometer from the surface protocol adopted for this study. Descriptions of the geologic basis for delineating each tract, the data used, the geologic criteria and rationale for the assessment, and results of the assessment are included in appendixes along with the description of a geographic information system (GIS) that includes tract boundaries, known porphyry copper deposits and significant prospects, and assessment results.

Scientific Investigations Report↗