USGS Science⌕ Search

SEARCH · USGS Science

Results for “Studies in Maine Geology”

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 1,135 records · Page 63Linked to original sources

Seismic instrumentation plan for the Hawaiian Volcano Observatory

The seismic network operated by the U.S. Geological Survey’s Hawaiian Volcano Observatory (HVO) is the main source of authoritative data for reporting earthquakes in the State of Hawaii, including those that occur on the State’s six active volcanoes (Kīlauea, Mauna Loa, Hualālai, Mauna Kea, Haleakalā, Lō‘ihi). Of these volcanoes, Kīlauea and Mauna Loa are considered “very high threat” in a report on the rationale for a National Volcanic Early Warning System (NVEWS) (Ewert and others, 2005). This seismic instrumentation plan assesses the current state of HVO’s seismic network with respect to the State’s active volcanoes and calculates the number of stations that are needed to upgrade the current network to provide a seismic early warning capability for forecasting volcanic activity. Further, the report provides proposed priorities for upgrading the seismic network and a cost assessment for both the installation costs and maintenance costs of the improved network that are required to fully realize the potential of the early warning system. HVO has operated seismometers on the Island of Hawai‘i since 1912. Currently, the seismic network includes more than 70 stations from four different organizations. Generally, the Island of Hawai‘i has most of the seismic stations in the network (and most of the activity), with the density of seismic stations increasing from the northern part of the island to the south-southeast. The strength of the current network, based on theoretical detection and location capabilities, is at the summit of Kīlauea Volcano and Kīlauea’s upper East Rift Zone and Pu‘u ‘Ō‘ō—where few, if any, upgrades need to be made to the seismic network. The network in the region between Kīlauea and Mauna Loa is slightly weaker, as is the summit of Mauna Loa. In general, the rift zones of each volcano are more poorly monitored seismically than the summits and thus require a greater number of stations to achieve a volcanic early warning capability for monitoring seismicity. Priorities for new seismic installations on the volcanoes depend on several factors, including current activity, historical activity, population exposure, and current network quality. On Kīlauea, new installations on the middle East Rift Zone, lower East Rift Zone, and lower Southwest Rift Zone appear to be the highest priorities. On Mauna Loa, improvements to the summit seismic network should be prioritized based on the analysis of the data, followed by the installation of a sparse network on both rift zones. Once installed, the next priority would be to create denser seismic networks on the rift zones, particularly where eruptions could quickly threaten populated areas (middle Northeast Rift Zone, lower Southwest Rift Zone). On Hualālai, analysis of the data indicates that the Northwest Rift Zone is the most important priority, particularly where it runs through the population center of Kalaoa. Hualālai’s South Rift Zone appears to be the lowest priority for additional seismic instrumentation of any rift zone on Kīlauea, Mauna Loa, or Hualālai because of its low historical activity and lack of population exposure. Mauna Kea and Haleakalā have less active historical eruptive activity and thus have more modest proposed upgrades to seismic instrumentation. The installation of new seismic stations is only the first part of building a volcanic early warning capability for seismicity in the State of Hawaii. Additional personnel will likely be required to study the volcanic processes at work under each volcano, analyze the current seismic activity at a level sufficient for early warning, build new tools for monitoring, maintain seismic computing resources, and maintain the new seismic stations.

Hawaii↗

Simulation of ground-water flow to assess geohydrologic factors and their effect on source-water areas for bedrock wells in Connecticut

Generic ground-water-flow simulation models show that geohydrologic factors?fracture types, fracture geometry, and surficial materials?affect the size, shape, and location of source-water areas for bedrock wells. In this study, conducted by the U.S. Geological Survey in cooperation with the Connecticut Department of Public Health, ground-water flow was simulated to bedrock wells in three settings?on hilltops and hillsides with no surficial aquifer, in a narrow valley with a surficial aquifer, and in a broad valley with a surficial aquifer?to show how different combinations of geohydrologic factors in different topographic settings affect the dimensions and locations of source-water areas in Connecticut. Three principal types of fractures are present in bedrock in Connecticut?(1) Layer-parallel fractures, which developed as partings along bedding in sedimentary rock and compositional layering or foliation in metamorphic rock (dips of these fractures can be gentle or steep); (2) unroofing joints, which developed as strain-release fractures parallel to the land surface as overlying rock was removed by erosion through geologic time; and (3) cross fractures and joints, which developed as a result of tectonically generated stresses that produced typically near-vertical or steeply dipping fractures. Fracture geometry is defined primarily by the presence or absence of layering in the rock unit, and, if layered, by the angle of dip in the layering. Where layered rocks dip steeply, layer-parallel fracturing generally is dominant; unroofing joints also are typically well developed. Where layered rocks dip gently, layer-parallel fracturing also is dominant, and connections among these fractures are provided only by the cross fractures. In gently dipping rocks, unroofing joints generally do not form as a separate fracture set; instead, strain release from unroofing has occurred along gently dipping layer-parallel fractures, enhancing their aperture. In nonlayered and variably layered rocks, layer-parallel fracturing is absent or poorly developed; fracturing is dominated by well-developed subhorizontal unroofing joints and steeply dipping, tectonically generated fractures and (or) cooling joints. Cross fractures (or cooling joints) in nonlayered and variably layered rocks have more random orientations than in layered rocks. Overall, nonlayered or variably layered rocks do not have a strongly developed fracture direction. Generic ground-water-flow simulation models showed that fracture geometry and other geohydrologic factors affect the dimensions and locations of source-water areas for bedrock wells. In general, source-water areas to wells reflect the direction of ground-water flow, which mimics the land-surface topography. Source-water areas to wells in a hilltop setting were not affected greatly by simulated fracture zones, except for an extensive vertical fracture zone. Source-water areas to wells in a hillside setting were not affected greatly by simulated fracture zones, except for the combination of a subhorizontal fracture zone and low bedrock vertical hydraulic conductivity, as might be the case where an extensive subhorizontal fracture zone is not connected or is poorly connected to the surface through vertical fractures. Source-water areas to wells in a narrow valley setting reflect complex ground-water-flow paths. The typical flow path originates in the uplands and passes through either till or bedrock into the surficial aquifer, although only a small area of the surficial aquifer actually contributes water to the well. Source-water areas in uplands can include substantial areas on both sides of a river. Source-water areas for wells in this setting are affected mainly by the rate of ground-water recharge and by the degree of anisotropy. Source-water areas to wells in a broad valley setting (bedrock with a low angle of dip) are affected greatly by fracture properties. The effect of a given fracture is to channel the

Scientific Investigations Report↗

Isotopic studies of the late Archean plutonic rocks of the Wind River Range, Wyoming

Isotopic studies of the Rb-Sr and U-Th-Pb systems in whole-rock samples and the U-Pb systematics for zircons document the existence of two late Arehean intrusive events in the Wind River Range. All of the systems examined indicate an age of ∼2,630 ± 20 m.y. for the Louis Lake batholith. Apparent ages for the Bears Ears pluton range from 2,504 ± 40 m.y. to 2,575 ± 50 m.y. The scatter in apparent ages for the Bears Ears pluton does not appear to be primarily the result of disturbance by postintrusive events, but it may be explained by an isotopically inhomogenous magma at the time of intrusion. Data for a few samples indicate that the Wind River Range was affected locally by a postmagmatic hydrothermal event that was approximately Tertiary in age. This event lowered δ 18 O values and disturbed parent-daughter relationships in most of the isotopic systems investigated, but it was recent enough that there is no demonstrable effect in the Pb-Pb system. The Bears Ears pluton has some chemical and petrologic features that are similar to those reported for the granites in the Granite Mountains to the east. These granites are spatially associated with low-temperature uranium deposits of Tertiary age and have been shown to have lost large amounts of uranium during the early to middle Tertiary. U-Pb systematics indicate, however, that the low to moderate uranium contents and highly variable Th/U values noted for the Bears Ears pluton are best interpreted as being primary features. If uranium was lost after magma generation, the loss most likely occurred at the time of intrusion. Such a loss could account for uraniferous Precambrian pegmatites southwest of the main part of the Range. The two intrusive units apparently were derived from different protoliths that were formed during early to middle Archean. Initial isotopic ratios and petrochemistry for the Louis Lake batholith are consistent with an early Archean trondhjemitic to tonalitic source. The protolith for the Bears Ears pluton must have been more evolved and somewhat younger. Inconsistencies as to the degree of evolution of this protolith, as inferred from isotopic and trace-element data, suggest that the protolith may have been subjected to high-grade meta-morphism that caused loss of Rb and U prior to generation of the magma.

Wyoming↗

Forecasting selenium discharges to the San Francisco Bay-Delta Estuary: Ecological effects of a proposed San Luis Drain extension

Selenium discharges to the San Francisco Bay-Delta Estuary (Bay-Delta) could change significantly if federal and state agencies (1) approve an extension of the San Luis Drain to convey agricultural drainage from the western San Joaquin Valley to the North Bay (Suisun Bay, Carquinez Strait, and San Pablo Bay); (2) allow changes in flow patterns of the lower San Joaquin River and Bay-Delta while using an existing portion of the San Luis Drain to convey agricultural drainage to a tributary of the San Joaquin River; or (3) revise selenium criteria for the protection of aquatic life or issue criteria for the protection of wildlife. Understanding the biotransfer of selenium is essential to evaluating effects of selenium on Bay-Delta ecosystems. Confusion about selenium threats to fish and wildlife stem from (1) monitoring programs that do not address specific protocols necessary for an element that bioaccumulates; and (2) failure to consider the full complexity of the processes that result in selenium toxicity. Past studies show that predators are more at risk from selenium contamination than their prey, making it difficult to use traditional methods to predict risk from environmental concentrations alone. This report presents an approach to conceptualize and model the fate and effects of selenium under various load scenarios from the San Joaquin Valley. For each potential load, progressive forecasts show resulting (1) water-column concentration; (2) speciation; (3) transformation to particulate form; (4) particulate concentration; (5) bioaccumulation by invertebrates; (6) trophic transfer to predators; and (7) effects on those predators. Enough is known to establish a first-order understanding of relevant conditions, biological response, and ecological risks should selenium be discharged directly into the North Bay through a conveyance such as a proposed extension of the San Luis Drain. The approach presented here, the Bay-Delta selenium model, determines the mass, fate, and effects of selenium released to the Bay-Delta through use of (1) historical land-use, drainage, alluvial-fill, and runoff databases; (2) existing knowledge concerning biogeochemical reactions and physiological parameters of selenium (e.g., speciation, partitioning between dissolved and particulate forms, and bivalve assimilation efficiency); and (3) site-specific data mainly from 1986 to 1996 for clams and bottom-feeding fish and birds. Selenium load scenarios consider effluents from North Bay oil refineries and discharges of agricultural drainage from the San Joaquin Valley to enable calculation of (a) a composite freshwater endmember selenium concentration at the head of the estuary; and (b) a selenium concentration at a selected seawater location (Carquinez Strait) as a foundation for modeling. Analysis of selenium effects also takes into account the mode of conveyance for agricultural drainage (i.e., the San Luis Drain or San Joaquin River); and flows of the Sacramento River and San Joaquin River on a seasonal or monthly basis. Load scenarios for San Joaquin Valley mirror predictions made since 1955 of a worsening salt (and by inference, selenium) build-up exacerbated by an arid climate and massive irrigation. The reservoir of selenium in the San Joaquin Valley is sufficient to provide loading at an annual rate of approximately 42,500 pounds of selenium to a Bay-Delta disposal point for 63 to 304 years at the lower range of projections presented here, even if influx of selenium from the California Coast Ranges could be curtailed. Disposal of wastewaters on an annual basis outside of the San Joaquin Valley may slow the degradation of valley resources, but drainage alone cannot alleviate the salt and selenium build-up in the San Joaquin Valley, at least within a century. Load scenarios also show the different proportions of selenium loading to the Bay-Delta. Oil refinery loads from 1986 to 1992 ranged from 8.5 to 20 pounds of selenium per day; with treatment and cleanup, loads decreased to 3.0 pounds of selenium per day in 1999. In contrast, San Joaquin Valley agricultural drainage loads disposed of in a San Luis Drain extension could range from 45 to 117 pounds of selenium per day across a set of historical and future conditions. Components of this valley-wide load include five source subareas (i.e., Grassland, Westlands, Tulare, Kern, and Northern) defined by water and drainage management. Loads vary per subarea mainly because of proximity of the subarea to geologic sources of selenium and irrigation history. Loads from the Sacramento River, depending on flow conditions, range from 0.8 to 10 pounds of selenium per day. Loads from the San Joaquin River vary depending on restoration and flow conditions, which are considered. A consistent picture of ecological risk emerges under modeled selenium discharges from a proposed San Luis Drain extension. The threat to the estuary is greatest during low flow seasons and critically dry years. Where selenium undergoes reactions typical of low flow or longer residence time, highly problematic bioaccumulation in prey (food) is forecast. Surf scoter, greater and lesser scaup, and white sturgeon appear to be most at risk because these Bay-Delta predators feed on deposit and filter-feeding bivalves. Recent findings add Sacramento splittail and Dungeness crab to that list. During the low flow season of critically dry years, forecasted selenium concentrations in water, particulate matter, prey (diet), and predator tissue exceed guidelines with a high certainty of producing adverse effects under the most likely load scenario from a proposed San Luis Drain extension. High flows afford some protection under certain conditions in modeled San Joaquin River scenarios. However, meeting a combined goal of releasing a specific load during maximum flows and keeping selenium concentrations in the river below a certain objective to protect against bioaccumulation may not always be attainable. Management of the San Joaquin River on a constant concentration basis also could create problematic bioaccumulation during a wet year, especially during the low flow season, because high flows translate to high loads that are not always offset by seasonal river inflows. Prior to refinery cleanup, selenium contamination was sufficient to threaten reproduction in key species within the Bay-Delta ecosystems and human health advisories were posted based on selenium concentrations in tissues of diving ducks. During this time, selenium concentrations in the Bay-Delta were well below the most stringent recommended water quality criterion [1 microgram per liter (1 µg/L)]. Enhanced biogeochemical transformations to bioavailable particulate selenium and efficient bioaccumulation by bivalves characterized the system. If these biogeochemical conditions continue to prevail and agricultural selenium sources replace or exceed refinery sources, ecological forecasts suggest the risk of adverse effects will be difficult to eliminate under an out-of-valley resolution to the selenium problem. The Bay-Delta selenium model presented here is a systematic approach for conducting forecasts of the ecological effects from selenium on aquatic food webs. It is a new tool that links and models the major processes leading from loads through consumer organisms to predators. It also is a feasible approach for site-specific analysis and could provide a framework for developing new protective selenium foodweb guidelines and predator criteria. Model components that help ensure understanding ecosystems and the basis of environmental protection are (1) contaminant concentrations and speciation in sources, such as particulate material, that most influence bioavailability; (2) bioaccumulation models that calculate concentrations in diet, specifically in bivalves of the Bay-Delta that act as sensitive indicators of selenium contamination; (3) food-web type that determines what animals are threatened and when; and (4) multiple media concentrations (water, particulate material, and tissue of prey and predators) that, in-combination, determine risk or hazard.

California↗

Stratigraphic and geotectonic relationships in northern Vermont and southern Quebec

Stratified rocks of early and middle Paleozoic age form a belt of northeast-trending anticlinoria and synclinoria of middle Paleozoic age, in northern Vermont and adjacent parts of southern Quebec . The foreland margin of this belt, in the Champlain and St. Lawrence valleys to the west, is cut by eastward-dipping thrust faults of middle Paleozoic age and by later northeast-trending normal faults. The Green Mountain anticlinorium, which is the dominant structure of the region, is flanked to the west, on the foreland, by the St. Albans-Hinesburg-Middlebury synclinorium and to the east, in the midst of the folded belt, by the Connecticut Valley-Gaspe synclinorium. The principal thrust faults, notably the Champlain and Philipsburg thrusts, are in the west limb of the St. Albans-Hinesburg-Middlebury synclinorium. East of the Connecticut Valley-Gaspe synclinorium is the Boundary Mountain anticlinorium, in eastern Vermont and adjacent New Hampshire and along the international boundary between Quebec and Maine. Two contrasting intergradational lithic assemblages, the graywacke-shale assemblage and the carbonate-quartzite assemblage, characterize the protolith of the bedded rocks. The graywacke-shale assemblage includes thick sections of lower Paleozoic strata, portions of which lap both gradationally and unconformably westward on the foreland, particularly in Quebec ; it also includes middle Paleozoic strata that offlap eastward away from the foreland. The carbonate-quartzite assemblage laps both unconformably and gradationally eastward over the graywacke-shale assemblage in sections of the middle Paleozoic east of the axis of the Green Mountain anticlinorium. Stratigraphic correlation has become well established in the foreland belt where numerous distinctive and fossiliferous strata, chiefly of the carbonate-quartzite assemblage, have escaped metamorphism. It is also fairly clear in sections in the eastern foreland and western part of the Green Mountain anticlinorium, where the strata of the carbonate-quartzite assemblage extend eastward and interfinger with rocks of the graywacke-shale terrane. Rocks that are entirely of the graywacke-shale assemblage have been correlated in the present study. The stratified rocks west of the axis of the Green Mountain anticlinorium are of Cambrian(?), Cambrian, and Ordovician age; those to the east range in age from Cambrian to Devonian. The geotectonic setting of the region is the once mobile belt of the Appalachian orthogeosyncline, which is at the southeastern margin of the stable continental block, or craton, of North America. The orthogeosyncline was a belt, chiefly of subsidence, that embraced two parallel and adjoining longitudinal zones: the eugeosynclinal zone, which was more mobile, and the miogeosynclinal zone, which was less mobile. Second- and third-generation geosynclines are superimposed not only on the orthogeosyncline but also on adjoining parts of the craton. Uplift, and finally folding, gradually superseded subsidence in the orthogeosyncline. Local uplift, chiefly within the eugeosynclinal zone, provided most of the clastic sediments, principally those of the graywacke-shale assemblage. General stabilization of the western part of the orthogeosyncline at the end of the Ordovician was accompanied by eastward migration of the miogeosynclinal zone. Localized uplift within the eastern part of the orthogeosyncline at this time is marked by unconformities referred to the Taconic disturbance. Folding and uplift after the Early Devonian is shown by angular unconformities referred to the Acadian and Appalachian orogenies. The interpretation of the geotectonic relations of the bedded rocks is aided by critical features of the magmatic activity that began with, accompanied, and followed the diastrophism.

Vermont, Quebec↗

Geology and undiscovered resource assessment of the potash-bearing, Middle Devonian (Givetian), Prairie Evaporite, Elk Point Basin, Canada and United States

The U.S. Geological Survey (USGS) assessed undiscovered potash resources in the Elk Point Basin in Canada and the United States as part of a global mineral resource assessment. The Elk Point Basin is a large, Middle Devonian (Givetian) intracratonic evaporite basin covering approximately 1,200,000 square kilometers (km 2 ) and filled mainly with marine evaporite and minor clastic sedimentary rocks that contain stratabound potash-bearing salt. The potash-bearing salt is concentrated in four stratigraphic members (Patience Lake, Belle Plaine, White Bear, and Esterhazy) in the upper 100 meters (m) of the Prairie Evaporite and are separated by beds of halite (NaCl) that contain lesser—presently non-economic—amounts of sylvite (KCl) and carnallite (KMgCl 3 ·6H 2 O). The principal ore-bearing salt contains mainly sylvite. Four permissive tracts were defined that permit the presence of undiscovered stratabound potash (both sylvite- and carnallite-bearing salt) using geological criteria. Permissive tracts are defined by the spatial extent of each stratigraphic member that is at least 1 m thick, are less than 3 kilometers (km) from the surface, contain at least 4 percent equivalent potassium oxide (K 2 O), and contain the currently known resources. The permissive tracts include known potash deposits and potash occurrences as wells or mines not in production and show where undiscovered potash resources may be present. Well data are used to define the extent, thickness, average K 2 O equivalent grades, and volumes of each member. Data were supplied by the Saskatchewan Geological Survey or were obtained from published National Instrument (NI) 43-101 technical reports and other published reports, such as annual 10-K reports or news releases. The Elk Point Basin is the world’s largest source of potash, producing 23.0 million metric tons (Mt) of potassium chloride (KCl) (the equivalent of about 14.4 Mt of K 2 O) in 2018. In terms of global importance, the Elk Point Basin may contain 40 to greater than 50 percent of the world’s potash resources. Since 1962, potash companies have mined more than 1.5 trillion metric tons of ore containing 605 Mt of KCl (the equivalent of about 380 Mt of K 2 O). The total value of the ore produced through 2018 is on the order of $70 trillion (CAD). Potash is currently produced from eight conventional and three underground solution mines at depths ranging from 900 m to nearly 1,800 m. Estimates of the amount of potash in the Elk Point Basin vary considerably and the data and methods used in those estimations are not well documented. Known potash resources are approximately 99 billion metric tons (Bt) of ore containing 22 Bt of K 2 O equivalent. As a result of new mine openings and increased production capacity at existing mines, the total production capacity of mines in the Elk Point Basin has increased significantly (to about 32.8 Mt of KCl or 22.8 Mt of K 2 O equivalent per year). Additional production capacity of about 31 Mt of KCl (or 17 Mt of K 2 O equivalent) per year could be realized over the next decade if several current (as of 2019) exploration and development projects reach production status. Stratabound potash-bearing salt of the Prairie Evaporite presently underlies a total area of about 188,000 km 2 and has a total volume of about 2,690 cubic kilometers (km 3 ). Post-depositional solution processes considerably modified the mineralogy and presence of the potash-bearing salt. These changes had a profound effect on the volume and grade of potash resources that remained in the Prairie Evaporite and are a major consideration of exploration and mining operations as well as in this assessment of undiscovered potash resources. This USGS assessment includes the locations and possible amounts of undiscovered potash resources in the Prairie Evaporite. Volumes for each stratigraphic member were computed using member thicknesses and areal extent modified by actual, estimated geologic loss owing to salt dissolution and extraction ratios, as well as estimated distribution of carnallite and sylvite. Both sylvite- and carnallite-bearing salts were assessed for potash in this study. The assessment uses modern published grade and tonnage data. The amount of undiscovered potash is estimated by using Monte Carlo simulations to combine volume estimates of the potash-bearing members with probability distributions for average grade and bulk density. Mean potash grades (expressed as percentage of K 2 O equivalent) calculated using drill core analyses are 17.76 for the Patience Lake Member, 15.98 for the Belle Plaine Member, 10.66 for the White Bear Member, and 15.30 for the Esterhazy Member. Geologic losses reported as extraction ratios during mining may range from 27.5 to 41.6 percent and are dependent on mining method and local geologic conditions. The assessment determined that mean estimated undiscovered K 2 O equivalent resources for the Patience Lake, Belle Plaine, White Bear, and Esterhazy Members are 340, 220, 34, and 190 Bt, respectively, and estimated a total mean of 790 Bt for the entire Prairie Evaporite above a depth of 3 km. The total mineralized rock tonnage is estimated to be about 5,000 Bt. Most of the assessed potash is located within Saskatchewan with lesser amounts in Alberta and Manitoba as well as Montana and North Dakota within the United States. Although carnallite is mined for potash in Europe, it has historically been avoided in mining plans for potash-producing companies in Saskatchewan because of mining, processing, and grade considerations. Carnallite-rich salt is locally present in concentrations and volumes that could be a significant resource of magnesium chloride (MgCl 2 ) obtained as a byproduct of processing the carnallite for potash. Previously estimated reserves (not NI 43-101 compliant) of mineralized material from 1955 to 2019 are 695 Mt at 22.1 percent MgCl 2 . The total amount of K 2 O equivalent as carnallite was estimated during this USGS assessment to be about 120 Bt (or 180 Bt KCl). With uncertainties in defining the areal extent of carnallite in each of the potash-bearing members, the amount of MgCl 2 as carnallite in the Elk Point Basin could be approximately 180 Bt.

Elk Point basin↗

Geology, water-quality, hydrology, and geomechanics of the Cuyama Valley groundwater basin, California, 2008--12

To assess the water resources of the Cuyama Valley groundwater basin in Santa Barbara County, California, a series of cooperative studies were undertaken by the U.S. Geological Survey and the Santa Barbara County Water Agency. Between 2008 and 2012, geologic, water-quality, hydrologic and geomechanical data were collected from selected sites throughout the Cuyama Valley groundwater basin. Geologic data were collected from three multiple-well groundwater monitoring sites and included lithologic descriptions of the drill cuttings, borehole geophysical logs, temperature logs, as well as bulk density and sonic velocity measurements of whole-core samples. Generalized lithologic characterization from the monitoring sites indicated the water-bearing units in the subsurface consist of unconsolidated to partly consolidated sand, gravel, silt, clay, and occasional cobbles within alluvial fan and stream deposits. Analysis of geophysical logs indicated alternating layers of finer- and coarser-grained material that range from less than 1 foot to more than 20 feet thick. On the basis of the geologic data collected, the principal water-bearing units beneath the monitoring-well sites were found to be composed of younger alluvium of Holocene age, older alluvium of Pleistocene age, and the Tertiary-Quaternary Morales Formation. At all three sites, the contact between the recent fill and younger alluvium is approximately 20 feet below land surface. Water-quality samples were collected from 12 monitoring wells, 27 domestic and supply wells, 2 springs, and 4 surface-water sites and were analyzed for a variety of constituents that differed by site, but, in general, included trace elements; nutrients; dissolved organic carbon; major and minor ions; silica; total dissolved solids; alkalinity; total arsenic and iron; arsenic, chromium, and iron species; and isotopic tracers, including the stable isotopes of hydrogen and oxygen, activities of tritium, and carbon-14 abundance. Of the 39 wells sampled, concentrations of total dissolved solids and sulfate from 38 and 37 well samples, respectively, were greater than the U.S. Environmental Protection Agency’s secondary maximum contaminant levels. Concentrations greater than the maximum contaminant levels for nitrate were observed in five wells and were observed for arsenic in four wells. Differences in the stable-isotopic values of hydrogen and oxygen among groundwater samples indicated that water does not move freely between different formations or between different zones within the Cuyama Valley. Variations in isotopic composition indicated that recharge is derived from several different sources. The age of the groundwater, expressed as time since recharge, was between 600 and 38,000 years before present. Detectable concentrations of tritium indicated that younger water, recharged since the early 1950s, is present in parts of the groundwater basin. Hydrologic data were collected from 12 monitoring wells, 56 domestic and supply wells, 3 surface-water sites, and 4 rainfall-gaging stations. Rainfall in the valley averaged about 8 inches annually, whereas the mountains to the south received between 12 and 19 inches. Stream discharge records showed seasonal variability in surface-water flows ranging from no-flow to over 1,500 cubic feet per second. During periods when inflow to the valley exceeds outflow, there is potential recharge from stream losses to the groundwater system Water-level records included manual quarterly depth-to-water measurements collected from 68 wells, time-series data collected from 20 of those wells, and historic water levels from 16 wells. Hydrographs of the manual measurements showed declining water levels in 16 wells, mostly in the South-Main zone, and rising water levels in 14 wells, mostly in the Southern Ventucopa Uplands. Time-series hydrographs showed daily, seasonal, and longer-term effects associated with local pumping. Water-level data from the multiple-well monitoring sites indicated seasonal fluctuations as great as 80 feet and water-level differences between aquifers as great as 40 feet during peak pumping season. Hydrographs from the multiple-well groundwater monitoring sites showed vertical hydraulic gradients were upward during the winter months and downward during the irrigation season. Historic hydrographs showed water-level declines in the Southern-Main, Western Basin, Caliente Northern-Main, and Southern Sierra Madre zone ranging from 1 to 7 feet per year. Hydrographs of wells in the Southern Ventucopa Uplands zone showed several years with marked increases in water levels that corresponded to increased precipitation in the Cuyama Valley. Investigation of hydraulic properties included hydraulic conductivity and transmissivity estimated from aquifer tests performed on 63 wells. Estimates of horizontal hydraulic conductivity ranged from about 1.5 to 28 feet per day and decreased with depth. The median estimated hydraulic conductivity for the older alluvium was about five times that estimated for the Morales Formation. Estimates of transmissivity ranged from 560 to 163,400 gallons per day per foot and decreased with depth. The median estimated transmissivity for the younger alluvium was about three times that estimated for the older alluvium. Geomechanical analysis included land-surface elevation changes at five continuously operating global positioning systems (GPS) and land-subsidence detection at five interferometric synthetic aperture radar (InSAR) reference points. Analysis of data collected from continuously operating GPS stations showed the mountains to the south and west moved upward about 1 millimeter (mm) annually, whereas the station in the center of the Southern-Main zone moved downward more than 7 mm annually, indicating subsidence. It is likely that this subsidence is inelastic (permanent) deformation and indicates reduced storage capacity in the aquifer sediments. Analysis of InSAR data showed local and regional changes that appeared to be dependent, in part, on the time span of the interferogram, seasonal variations in pumping, and tectonic uplift. Long-term InSAR time series showed a total maximum detected subsidence rate of approximately 12 mm per year at one location and approximately 8 mm per year at a second location, while short-term InSAR time series showed maximum subsidence of about 15 mm at one location and localized maximum uplift of about 10 mm at another location.

California↗

Peak-flow frequency estimates based on data through water year 2001 for selected streamflow-gaging stations in South Dakota

Numerous users, including the South Dakota Department of Transportation, have continuing needs for peak-flow information for the design of highway infrastructure and many other purposes. This report documents results from a cooperative study between the South Dakota Department of Transportation and the U.S. Geological Survey to provide an update of peak-flow frequency estimates for South Dakota. Estimates of peak-flow magnitudes for 2-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence intervals are reported for 272 streamflow-gaging stations, which include most gaging stations in South Dakota with 10 or more years of systematic peak-flow records through water year 2001. Recommended procedures described in Bulletin 17B were used as primary guidelines for developing peak-flow frequency estimates. The computer program PEAKFQ developed by the U.S. Geological Survey was used to run the frequency analyses. Flood frequencies for all stations were initially analyzed by using standard Bulletin 17B default procedures for fitting the log-Pearson III distribution. The resulting preliminary frequency curves were then plotted on a log-probability scale, and fits of the curves with systematic data were evaluated. In many cases, results of the default Bulletin 17B analyses were determined to be satisfactory. In other cases, however, the results could be improved by using various alternative procedures for frequency analysis. Alternative procedures for some stations included adjustments to skew coefficients or use of user-defined low-outlier criteria. Peak-flow records for many gaging stations are strongly influenced by low- or zero-flow values. This situation often results in a frequency curve that plots substantially above the systematic record data points at the upper end of the frequency curve. Adjustments to low-outlier criteria reduced the influence of very small peak flows and generally focused the analyses on the upper parts of the frequency curves (10- to 500-year recurrence intervals). The most common alternative procedures involved several different methods to extend systematic records, which was done primarily to address biases resulting from nonrepresentative climatic conditions during several specific periods of record and to reduce inconsistencies among multiple gaging stations along common stream channels with different periods of record. In some cases, records for proximal stations could be combined directly. In other cases, the two-station comparison procedure recommended in Bulletin 17B was used to adjust the mean and standard deviation of the logs of the systematic data for a target station on the basis of correlation with concurrent records from a nearby long-term index station. In some other cases, a 'mixed-station procedure' was used to adjust the log-distributional parameters for a target station, on the basis of correlation with one or more index stations, for the purpose of fitting the log-Pearson III distribution. Historical adjustment procedures were applied to peak-flow frequency analyses for 17 South Dakota gaging stations. A historical adjustment period extending back to 1881 (121 years) was used for 12 gaging stations in the James and Big Sioux River Basins, and various other adjustment periods were used for additional stations. Large peak flows that occurred in 1969 and 1997 accounted for 13 of the 17 historical adjustments. Other years for which historical peak flows were used include 1957, 1962, 1992, and 2001. A regional mixed-population analysis was developed to address complications associated with many high outliers for the Black Hills region. This analysis included definition of two populations of flood events. The population of flood events that composes the main body of peak flows for a given station is considered the 'ordinary-peaks population,' and the population of unusually large peak flows that plot substantially above the main body of peak flows on log-probability scale is co

South Dakota↗

Primary deposition and early diagenetic effects on the high saturation accumulation of gas hydrate in a silt dominated reservoir in the Gulf of Mexico

On continental margins, high saturation gas hydrate systems (>60% pore volume) are common in canyon and channel environments within the gas hydrate stability zone, where reservoirs are dominated by coarse-grained, high porosity sand deposits. Recent studies, including the results presented here, suggest that rapidly deposited, silt-dominated channel-levee environments can also host high saturation gas hydrate accumulations. Here we present several sedimentological data sets, including sediment composition, biostratigraphic age from calcareous nannofossils, grain size, total organic carbon (TOC), C/N elemental ratio, δ 13 C-TOC , CaCO 3 , total sulfur (TS), and δ 34 S-TS from sediments collected with pressure cores from a gas hydrate rich, turbidite channel-levee system in the Gulf of Mexico during the 2017 UT-GOM2-1 Hydrate Pressure Coring Expedition. Our results indicate the reservoir is composed of three main lithofacies, which have distinct sediment grain size distributions (type A-silty clay to clayey silt, type B-clayey silt, and type C-sandy silt to silty sand) that are characteristic of variable turbidity current energy regimes within a Pleistocene (< 0.91 Ma) channel-levee environment. We document that the TOC in the sediments of the reservoir is terrestrial in origin and contained within the fine fraction of each lithofacies, while the CaCO 3 fraction is composed of primarily reworked grains, including Cretaceous calcareous nannofossils, and part of the detrital load. The lack of biogenic grains within the finest grained sediment intervals throughout the reservoir suggests interevent hemipelagic sediments are not preserved, resulting in a reservoir sequence of silt dominated, stacked turbidites. We observe two zones of enhanced TS at the top and bottom of the reservoir that correspond with enriched bulk sediment δ 34 S, indicating stalled or slowly advancing paleo-sulfate-methane transition zone (SMTZ) positions likely driven by relative decreases in sedimentation rate. Despite these two diagenetic zones, the low abundance of diagenetic precipitates throughout the reservoir allowed the primary porosity to remain largely intact, thus better preserving primary porosity for subsequent pore-filling gas hydrate. In canyon, channel, and levee environments, early diagenesis may be regulated via sedimentation rates, where high rates result in rapid progression through the SMTZ and minimal diagenetic mineralization and low rates result in the stalling of the SMTZ, enhancing diagenetic mineralization. Here, we observed some enhanced pyritization to implicate potential sedimentation rate changes, but not enough to consume primary porosity, resulting in a high saturation gas hydrate reservoir. These results emphasize the important implications of sedimentary processes, sedimentation rates, and early diagenesis on the distribution of gas hydrate in marine sediments along continental margins.

Alabama, Louisiana, Mississippi↗

Sources, instream transport, and trends of nitrogen, phosphorus, and sediment in the lower Tennessee River basin, 1980-96

In 1997, the U.S. Geological Survey (USGS) began an assessment of the lower Tennessee River Basin as part of the National Water-Quality Assessment Program. Existing nutrient and sediment data from 1980 to 1996 were compiled, screened, and interpreted to estimate watershed inputs from nutrient sources, provide a general description of the distribution and transport of nutrients and sediments in surface water, and evaluate trends in nutrient and sediment concentrations in the lower Tennessee (LTEN) River Basin. Nitrogen inputs from major sources varied widely among tributary basins in the LTEN River Basin. Point source wastewater discharges contributed between 0 and 0.61 tons per square mile per year [(tons/mi 2 )/yr]. Of the nonpoint sources of nitrogen for which inputs were estimated (atmospheric deposition, nitrogen fixation, fertilizer application, and livestock waste) livestock waste contributed the largest input in about two-thirds (7 out of 11) of the tributary basins, and fertilizer application contributed the largest input in the remaining 4 basins. Nitrogen input from fertilizer application was the most variable spatially among the nonpoint sources of nitrogen, ranging from 1.5 to 23 (tons/mi 2 )/yr. Atmospheric deposition estimates varied the least from basin to basin, ranging from 1.6 to 2.0 (tons/mi 2 )/yr. Estimates of nitrogen input from livestock waste ranged between 2.0 to 13 (tons/mi 2 )/yr. The percentage of the input from each of these nonpoint sources that entered the surface-water system is not known. Wastewater discharge contributed between 0 and 0.14 (ton/mi 2 )/yr of phosphorus to tributary basins. Livestock waste contributed most of the input in 8 out of the 11 basins, and fertilizer application contributed the most in the remaining 3 basins. Estimates of phosphorus input for fertilizer application ranged from 0.35 to 5.1 (tons/mi 2 )/yr and from 0.62 to 4.3 (tons/mi 2 )/yr from livestock waste. Reservoirs on the main stem of the Tennessee River and on the Duck and Elk Rivers affect nutrient transport because hydrodynamic conditions in the reservoirs promote assimilation by aquatic plants and deposition of particulate matter. Observed decreases in total nitrite plus nitrate and dissolved-orthophosphorus concentrations in reservoirs or at sites downstream of reservoirs during summer months were probably related to seasonality of plant growth. Nutrient and sediment data used to estimate annual instream loads and yields were compiled from various water-quality monitoring programs and represent the best available data in the LTEN River Basin, but these data have several characteristics that limit accuracy of load estimates. Many of the monitoring programs were not designed with the objective of annual load estimation, and data representing storm transport are, therefore, sparse; sampling and analytical methods varied through time and among the monitoring programs, hampering spatial and temporal comparisons. The load estimates computed from these data are useful for evaluating broad spatial patterns of instream load, and comparisons of instream load to inputs, but may not be sufficiently accurate for local-scale evaluations of water quality. Estimates of the mean annual instream load of total nitrogen entering (Chattanooga, Tenn.) and leaving (Paducah, Ky.) the LTEN River Basin were 29,000 and 60,000 tons per year (tons/yr), respectively. These estimates represent a gain of 31,000 tons/yr, on average, across the area (18,930 mi 2 ) between these inlet and outlet sites. The sum of the mean annual instream load from gaged tributaries to the main stem within the study unit was 14,000 tons/yr; however, this number cannot be directly compared with the gain between the inlet and outlet sites because (1) the gaged area represents only 30 percent of the total area and (2) the period of record at many tributary sites did not correspond with the period of record at the inlet or outlet sites. Estimates of mean annual instream load of total phosphorus at the inlet and outlet sites of the LTEN River Basin were 1,300 and 5,000 tons/yr, respectively, representing a gain of 3,700 tons/yr, on average, across the study unit. The sum of the gaged tributary load, representing only 28 percent of the area contributing to the main stem, was 4,300 tons/yr. Although this number cannot be closely compared with the gain throughout the study unit, for the same reasons given for total nitrogen, a general comparison suggests that the main stem of the Tennessee River and the tributary embayments along the main stem function as a sink for total phosphorus, removing a substantial amount from the water column through deposition or assimilation. The estimates of inputs can be compared and correlated with yields (area-normalized instream loads); significant correlations between estimates of inputs and yields might be useful as predictive tools for instream water quality where monitoring data are not available. Yields of nitrogen correlated moderately well with inputs from nonpoint sources, based on 1992 estimates. Nitrogen yield was highest [3.5 (tons/mi 2 )/yr] for Town Creek, for which the balance of nonpoint-source inputs to agricultural lands (fertilizer application plus nitrogen fixation plus livestock waste minus harvest) was also the highest [15 (tons/mi 2 )/yr]. Nitrogen yield was low [1.0 (tons/mi 2 )/yr] for the Buffalo River, for which the balance of agricultural nonpoint-source input was correspondingly low [3.2 (tons/mi 2 )/yr, the second lowest]. Correlation of wastewater discharge with yield was poor, and contrasted with the significant correlation between wastewater discharge and median nitrogen concentration during low streamflow. The poor correlation between wastewater discharge and annual yield was expected, however, as wastewater discharge is a small fraction compared with annual yield. In contrast with nitrogen, phosphorus yield did not correlate well with any estimated inputs or land-use types for the tributary basins. Phosphorus yield was highest [1.1 and 0.93 (tons/mi 2 )/yr] at two sites along the Duck River and at Elk River near Prospect [0.89 (ton/mi 2 )/yr]; however, estimates of inputs at these sites were in the middle of their respective ranges. The influence of the outcrop of phosphatic limestone formations of the brown-phosphate districts in the lower Duck and lower Elk River Basins might be responsible for the poor correlation between estimated inputs and yields of phosphorus. The outcrop pattern of these phosphatic limestones are an important factor to consider as regional boundaries are established for attainable, region-specific water-quality criteria for total phosphorus. Estimates of sediment input from cropland soil erosion in 1992 ranged from 51 to 540 (tons/mi 2 )/yr among the major hydrologic units in the LTEN River Basin. Information was not available to estimate this input for individual tributaries. Sediment yield estimates ranged from 65 to 263 (tons/mi 2 )/yr for the three tributary monitoring basins for which instream data were available, and from 17 to 26 (tons/mi 2 )/yr for the Tennessee River at South Pittsburg and at Pickwick Landing Dam, respectively. Lower sediment yields for the main stem sites compared with the tributary sites is probably due to sediment deposition in the main stem of the Tennessee River and tributary embayments along the main stem. Most of the significant trends in nutrient concentrations from about 1985 to about 1995 were decreasing trends, except for total nitrite plus nitrate, which increased at one site on the Elk River. The spatial distribution of decreasing trends of total nitrogen and total ammonia corresponds with the spatial variation among basins in wastewater loading rate. The time period of observed trends corresponds to the period of improvements in municipal treatment, thus decreases in wastewater effluent concentrations of nitrogen might be responsible for the decreasing trend in instream concentrations at these sites. Concentrations of total phosphorus did not decrease during this period at these sites, as might have been expected considering the reductions in wastewater input of phosphorus during this period.

Alabama, Georgia, Kentucky, Mississippi, Tennessee↗

Detrital zircon geochronology of Cretaceous and Paleogene strata across the south-central Alaskan convergent margin

Ages of detrital zircons are reported from ten samples of Lower Cretaceous to Paleogene metasandstones and sandstones from the Chugach Mountains, Talkeetna Mountains, and western Alaska Range of south-central Alaska. Zircon ages are also reported from three igneous clasts from two conglomerates. The results bear on the regional geology, stratigraphy, tectonics, and mineral resource potential of the southern Alaska convergent margin. Chugach Mountains - The first detrital zircon data are reported here from the two main components of the Chugach accretionary complex - the inboard McHugh Complex and the outboard Valdez Group. Detrital zircons from sandstone and two conglomerate clasts of diorite were dated from the McHugh Complex near Anchorage. This now stands as the youngest known part of the McHugh Complex, with an inferred Turonian (Late Cretaceous) depositional age no older than 91-93 Ma. The zircon population has probability density peaks at 93 and 104 Ma and a smattering of Early Cretaceous and Jurassic grains, with nothing older than 191 Ma. The two diorite clasts yielded Jurassic U-Pb zircon ages of 179 and 181 Ma. Together, these findings suggest a Mesozoic arc as primary zircon source, the closest and most likely candidate being the Wrangellia composite terrane. The detrital zircon sample from the Valdez Group contains zircons as young as 69 and 77 Ma, consistent with the previously assigned Maastrichtian to Campanian (Late Cretaceous) depositional age. The zircon population has peaks at 78, 91, 148, and 163 Ma, minor peaks at 129, 177, 330, and 352 Ma, and no concordant zircons older than Devonian. A granite clast from a Valdez Group conglomerate yielded a Triassic U-Pb zircon age of 221 Ma. Like the McHugh Complex, the Valdez Group appears to have been derived almost entirely from Mesozoic arc sources, but a few Precambrian zircons are also present. Talkeetna Mountains - Detrital zircons ages were obtained from southernmost metasedimentary rocks of the Talkeetna Mountains (schist of Hatcher Pass) and, immediately to the south, the northernmost sedimentary sequence of the Matanuska forearc basin (Arkose Ridge Formation). Detrital zircons from the Paleogene Arkose Ridge Formation are as young as 61 and 70 Ma; the population is dominated by a single Late Cretaceous peak at 76 Ma; the oldest zircon is 181 Ma. Sedimentological evidence clearly shows that the conglomeratic Arkose Ridge Formation was derived from the Talkeetna Mountains; our detrital zircon data support this inference. Zircons dated at ca. 90 Ma in the Arkose Ridge sample suggest that buried or unmapped plutons of this age may exist in the Talkeetnas. This is a particularly interesting age as it corresponds to the age of the supergiant Pebble gold-molybdenum-copper porphyry prospect near Iliamna and suggests a new area of prospectivity for Pebble-type deposits. The schist of Hatcher Pass, which was previously assigned a Jurassic depositional age, yielded surprisingly young Late Cretaceous detrital zircons, the youngest at 75 Ma. The probability density curve has four Cretaceous peaks from 76 to 102 Ma, a pair of Late Jurassic peaks at 155 and 166 Ma, three Early Jurassic to Late Triassic peaks at 186, 197, and 213 Ma, minor Carboniferous peaks at 303 and 346 Ma, and a minor Paleoproterozoic peak at 1828 Ma. The schist of Hatcher Pass was largely derived from Mesozoic arc sources, most likely the Wrangellia composite terrane, with some contribution from one or more older, inboard sources, probably including the Yukon-Tanana terrane. We postulate that the schist of Hatcher Pass represents metamorphosed rocks of the Valdez Group that were subducted and then exhumed along the Chugach terrane's 'backstop' during Paleogene transtension. Western Alaska Range - Six detrital zircon samples were collected from a little studied belt of turbidites in Tyonek quadrangle on strike with the Kahiltna assemblage of the central Alaska Range.

Alaska↗

Aeromagnetic map and interpretation of geophysical data from the Condrey Mountain Roadless Area, Siskiyou County, California

The Condrey Mountain Roadless Area lies within the Klamath Mountains geologic province, a province composed of four eastward-dipping imbricate thrust slices or belts consisting predominantly of marine arc-related volcanic and sedimentary rocks, but also including significant amounts of ultramafic and other ophiolitic rocks (Irwin, 1981). From west to east the four Klamath Mountains thrust slices are called the western Jurassic belt, the western Paleozoic and Triassic belt, the central metamorphic belt, and the eastern Klamath belt. The Condrey Mountain Road less Area is located within a structural window in the western Paleozoic and Triassic belt, a window in which the structural dome comprised of the Condrey Mountain Schist is exposed (Coleman and others, 1983). North, east and south of the roadless area, the Condrey Mountain Schist is separated from the overlying western Paleozoic and Triassic belt by a low-angle regional thrust, whereas to the west the two units are in high-angle fault contact (Coleman and others, 1983). The Condrey Mountain schist consists mainly of sedimentary and volcanic rocks metamorphosed to greenschist facies. Metasedimentary rocks are exposed over most of the eastern and central parts of the structural window. Metavolcanic rocks occupy the western part of the window but also occur in small exposures within the metasedimentary rocks and along the eastern and southern margins of the window. Tabular bodies of metaserpentinite, the largest of which crops out near White Mountain, are contained within the metasedimentary rocks. Densities of hand samples from 11 sites scattered throughout the Condrey Mountain window average 2.66&plusmn;0.05 g/cm 3 . The four samples of metavolcanic rocks yielded a higher average density than the seven samples of metasedimentary rocks (2.71&plusmn;0.04 g/cm 3 versus 2.63&plusmn;0.03 g/cm 3 ). Along the western edge of the Condrey Mountain Road less Area, numerous narrow north-trending zones of mineralized schist extend from near Copper Butte on the south, northward to Elliot Creek (Coleman and others, 1983). These zones contain abundant pyrite mineralization associated with pyrrhotite, chalcopyrite, qalena, and sphalerite. The Blue Ledge mine area, 2 mi north-northeast of Copper Butte, is located in a zone of mineralized schist. The western Paleozoic and Triassic belt that nearly surrounds the Condrey Mountain Schist is a melange of sedimentary, volcanic, and ultramafic rocks metamorphosed to amphibolite facies (Coleman and others, 1983). Only two samples of the metamorphic melange were collected near the Condrcy Mountain Road less Area, but extensive sampling of this unit southwest of the roadless area yielded an average sample density of 2.86&plusmn;0.15 g/cm 3 (112 samples) (Jachens and others, 1983).

California↗

Application of health-based screening levels to ground-water quality data in a state-scale pilot effort

A state-scale pilot effort was conducted to evaluate a Health-Based Screening Level (HBSL) approach developed for communicating findings from the U.S. Geological Survey (USGS) National Water-Quality Assessment Program in a human-health context. Many aquifers sampled by USGS are used as drinking-water sources, and water-quality conditions historically have been assessed by comparing measured contaminant concentrations to established drinking-water standards and guidelines. Because drinking-water standards and guidelines do not exist for many analyzed contaminants, HBSL values were developed collaboratively by the USGS, U.S. Environmental Protection Agency (USEPA), New Jersey Department of Environmental Protection, and Oregon Health & Science University, using USEPA toxicity values and USEPA Office of Water methodologies. The main objective of this report is to demonstrate the use of HBSL approach as a tool for communicating water-quality data in a human-health context by conducting a retrospective analysis of ground-water quality data from New Jersey. Another important objective is to provide guidance on the use and interpretation of HBSL values and other human-health benchmarks in the analyses of water-quality data in a human-health context. Ground-water samples collected during 1996-98 from 30 public-supply, 82 domestic, and 108 monitoring wells were analyzed for 97 pesticides and 85 volatile organic compounds (VOCs). The occurrence of individual pesticides and VOCs was evaluated in a human-health context by calculating Benchmark Quotients (BQs), defined as ratios of measured concentrations of regulated compounds (that is, compounds with Federal or state drinking-water standards) to Maximum Contaminant Level (MCL) values and ratios of measured concentrations of unregulated compounds to HBSL values. Contaminants were identified as being of potential human-health concern if maximum detected concentrations were within a factor of 10 of the associated MCL or HBSL (that is, maximum BQ value (BQmax) greater than or equal to 0.1) in any well type (public supply, domestic, monitoring). Most (57 of 77) pesticides and VOCs with human-health benchmarks were detected at concentrations well below these levels (BQmax less than 0.1) for all three well types; however, BQmax values ranged from 0.1 to 3,000 for 6 pesticides and 14 VOCs. Of these 20 contaminants, one pesticide (dieldrin) and three VOCs (1,2-dibromoethane, tetrachloroethylene, and trichloroethylene) both (1) were measured at concentrations that met or exceeded MCL or HBSL values, and (2) were detected in more than 10 percent of samples collected from raw ground water used as sources of drinking water (public-supply and (or) domestic wells) and, therefore, are particularly relevant to human health. The occurrence of multiple pesticides and VOCs in individual wells also was evaluated in a human-health context because at least 53 different contaminants were detected in each of the three well types. To assess the relative human-health importance of the occurrence of multiple contaminants in different wells, the BQ values for all contaminants in a given well were summed. The median ratio of the maximum BQ to the sum of all BQ values for each well ranged from 0.83 to 0.93 for all well types, indicating that the maximum BQ makes up the majority of the sum for most wells. Maximum and summed BQ values were statistically greater for individual public-supply wells than for individual domestic and monitoring wells. The HBSL approach is an effective tool for placing water-quality data in a human-health context. For 79 of the 182 compounds analyzed in this study, no USEPA drinking-water standards or guidelines exist, but new HBSL values were calculated for 39 of these 79 compounds. The new HBSL values increased the number of detected pesticides and VOCs with human-health benchmarks from 65 to 77 (of 97 detected compounds), thereby expanding the basis for interpreting contaminant-occu

Scientific Investigations Report↗

Strata-bound Fe-Co-Cu-Au-Bi-Y-REE deposits of the Idaho Cobalt Belt: Multistage hydrothermal mineralization in a magmatic-related iron oxide copper-gold system

Mineralogical and geochemical studies of strata-bound Fe-Co-Cu-Au-Bi-Y-rare-earth element (REE) deposits of the Idaho cobalt belt in east-central Idaho provide evidence of multistage epigenetic mineralization by magmatic-hydrothermal processes in an iron oxide copper-gold (IOCG) system. Deposits of the Idaho cobalt belt comprise three types: (1) strata-bound sulfide lenses in the Blackbird district, which are cobaltite and, less commonly, chalcopyrite rich with locally abundant gold, native bismuth, bismuthinite, xenotime, allanite, monazite, and the Be-rich silicate gadolinite-(Y), with sparse uraninite, stannite, and Bi tellurides, in a gangue of quartz, chlorite, biotite, muscovite, garnet, tourmaline, chloritoid, and/or siderite, with locally abundant fluorapatite or magnetite; (2) discordant tourmalinized breccias in the Blackbird district that in places have concentrations of cobaltite, chalcopyrite, gold, and xenotime; and (3) strata-bound magnetite-rich lenses in the Iron Creek area, which contain cobaltiferous pyrite and locally sparse chalcopyrite or xenotime. Most sulfide-rich deposits in the Blackbird district are enclosed by strata-bound lenses composed mainly of Cl-rich Fe biotite; some deposits have quartz-rich envelopes. Whole-rock analyses of 48 Co- and/or Cu-rich samples show high concentrations of Au (up to 26.8 ppm), Bi (up to 9.16 wt %), Y (up to 0.83 wt %), ∑REEs (up to 2.56 wt %), Ni (up to 6,780 ppm), and Be (up to 1,180 ppm), with locally elevated U (up to 124 ppm) and Sn (up to 133 ppm); Zn and Pb contents are uniformly low (≤821 and ≤61 ppm, respectively). Varimax factor analysis of bulk compositions of these samples reveals geochemically distinct element groupings that reflect statistical associations of monazite, allanite, and xenotime; biotite and gold; detrital minerals; chalcopyrite and sparse stannite; quartz; and cobaltite with sparse selenides and tellurides. Significantly, Cu is statistically separate from Co and As, consistent with the general lack of abundant chalcopyrite in cobaltite-rich samples. Paragenetic relations determined by scanning electron microscopy indicate that the earliest Y-REE-Be mineralization preceded deposition of Co, Cu, Au, and Bi. Allanite, xenotime, and gadolinite-(Y) commonly occur as intergrowths with and inclusions in cobaltite; the opposite texture is rare. Monazite, in contrast, forms a poikiloblastic matrix to cobaltite and thin rims on allanite and xenotime, reflecting a later metamorphic paragenesis. Allanite and xenotime also show evidence of late dissolution and reprecipitation, forming discordant rims on older anhedral allanite and xenotime and separate euhedral crystals of each mineral. Textural data suggest extensive deformation of the deposits by folding and shearing, and by pervasive recrystallization, all during Cretaceous metamorphism. Sensitive high resolution ion microprobe U-Pb geochronology by Aleinikoff et al. (2012) supports these paragenetic interpretations, documenting contemporaneous Mesoproterozoic growth of early xenotime and crystallization of megacrystic A-type granite on the northern border of the district. These ages are used together with mineralogical and geochemical data from the present study to support an epigenetic, IOCG model for Fe-Co-Cu-Au-Bi-Y-REE deposits of the Idaho cobalt belt. A sulfide facies variant of IOCG deposits is proposed for the Blackbird district, in which reducing hydrothermal conditions favored deposition of sulfide minerals over iron oxides. This new model includes Mesoproterozoic vein mineralization and related Fe-Cl metasomatism that formed the biotite-rich lenses, a predominantly felsic magmatic source for metals in the deposits, given their local abundance of Y, REEs, and Be, and a major sedimentary component in the hydrothermal fluids based on independent sulfur isotope and boron isotope data for sulfides and ore-related tourmaline, respectively.

Economic Geology↗

Potential sources of analytical bias and error in selected trace element data-quality analyses

Potential sources of analytical bias and error associated with laboratory analyses for selected trace elements where concentrations were greater in filtered samples than in paired unfiltered samples were evaluated by U.S. Geological Survey (USGS) Water Quality Specialists in collaboration with the USGS National Water Quality Laboratory (NWQL) and the Branch of Quality Systems (BQS). Causes for trace-element concentrations in filtered samples to exceed those in associated unfiltered samples have been attributed to variability in analytical measurements, analytical bias, sample contamination either in the field or laboratory, and (or) sample-matrix chemistry. These issues have not only been attributed to data generated by the USGS NWQL but have been observed in data generated by other laboratories. This study continues the evaluation of potential analytical bias and error resulting from matrix chemistry and instrument variability by evaluating the performance of seven selected trace elements in paired filtered and unfiltered surface-water and groundwater samples collected from 23 sampling sites of varying chemistries from six States, matrix spike recoveries, and standard reference materials. Filtered and unfiltered samples have been routinely analyzed on separate inductively coupled plasma-mass spectrometry instruments. Unfiltered samples are treated with hydrochloric acid (HCl) during an in-bottle digestion procedure; filtered samples are not routinely treated with HCl as part of the laboratory analytical procedure. To evaluate the influence of HCl on different sample matrices, an aliquot of the filtered samples was treated with HCl. The addition of HCl did little to differentiate the analytical results between filtered samples treated with HCl from those samples left untreated; however, there was a small, but noticeable, decrease in the number of instances where a particular trace-element concentration was greater in a filtered sample than in the associated unfiltered sample for all trace elements except selenium. Accounting for the small dilution effect (2 percent) from the addition of HCl, as required for the in-bottle digestion procedure for unfiltered samples, may be one step toward decreasing the number of instances where trace-element concentrations are greater in filtered samples than in paired unfiltered samples. The laboratory analyses of arsenic, cadmium, lead, and zinc did not appear to be influenced by instrument biases. These trace elements showed similar results on both instruments used to analyze filtered and unfiltered samples. The results for aluminum and molybdenum tended to be higher on the instrument designated to analyze unfiltered samples; the results for selenium tended to be lower. The matrices used to prepare calibration standards were different for the two instruments. The instrument designated for the analysis of unfiltered samples was calibrated using standards prepared in a nitric:hydrochloric acid (HNO 3 :HCl) matrix. The instrument designated for the analysis of filtered samples was calibrated using standards prepared in a matrix acidified only with HNO 3 . Matrix chemistry may have influenced the responses of aluminum, molybdenum, and selenium on the two instruments. The best analytical practice is to calibrate instruments using calibration standards prepared in matrices that reasonably match those of the samples being analyzed. Filtered and unfiltered samples were spiked over a range of trace-element concentrations from less than 1 to 58 times ambient concentrations. The greater the magnitude of the trace-element spike concentration relative to the ambient concentration, the greater the likelihood spike recoveries will be within data control guidelines (80–120 percent). Greater variability in spike recoveries occurred when trace elements were spiked at concentrations less than 10 times the ambient concentration. Spike recoveries that were considerably lower than 90 percent often were associated with spiked concentrations substantially lower than what was present in the ambient sample. Because the main purpose of spiking natural water samples with known quantities of a particular analyte is to assess possible matrix effects on analytical results, the results of this study stress the importance of spiking samples at concentrations that are reasonably close to what is expected but sufficiently high to exceed analytical variability. Generally, differences in spike recovery results between paired filtered and unfiltered samples were minimal when samples were analyzed on the same instrument. Analytical results for trace-element concentrations in ambient filtered and unfiltered samples greater than 10 and 40 μg/L, respectively, were within the data-quality objective for precision of ±25 percent. Ambient trace-element concentrations in filtered samples greater than the long-term method detection limits but less than 10 μg/L failed to meet the data-quality objective for precision for at least one trace element in about 54 percent of the samples. Similarly, trace-element concentrations in unfiltered samples greater than the long-term method detection limits but less than 40 μg/L failed to meet this data-quality objective for at least one trace-element analysis in about 58 percent of the samples. Although, aluminum and zinc were particularly problematic, limited re-analyses of filtered and unfiltered samples appeared to improve otherwise failed analytical precision. The evaluation of analytical bias using standard reference materials indicate a slight low bias for results for arsenic, cadmium, selenium, and zinc. Aluminum and molybdenum show signs of high bias. There was no observed bias, as determined using the standard reference materials, during the analysis of lead.

Scientific Investigations Report↗

Depletion of 13C in Cretaceous marine organic matter: Source, diagenetic, or environmental sigal?

Geochemical studies of Cretaceous strata rich in organic carbon (OC) from Deep Sea Drilling Project (DSDP) sites and several land sections reveal several consistent relationships among amount of OC, hydrocarbon generating potential of kerogen (measured by pyrolysis as the hydrogen index, HI), and the isotopic composition of the OC. First, there is a positive correlation between HI and OC in strata that contain more than about 1% OC. Second, percent OC and HI often are negatively correlated with carbon isotopic composition (σ 13 C) of kerogen. The relationship between HI and OC indicates that as the amount of organic matter increases, this organic matter tends to be more lipid rich reflecting the marine source of the organic matter. Cretaceous samples that contain predominantly marine organic matter tend to be isotopically lighter than those that contain predominantly terrestrial organic matter. Average σ 13 C values for organic matter from most Cretaceous sites are between −26 and −28‰, and values heavier than about −25‰ occur at very few sites. Most of the σ 13 C values of Miocene to Holocene OC-rich strata and modern marine plankton are between −16 to −23‰. Values of σ 13 C of modern terrestrial organic matter are mostly between −23 and −33‰. The depletion of terrestial OC in 13 C relative to marine planktonic OC is the basis for numerous statements in the literature that isotopically light Cretaceous organic matter is of terrestrial origin, even though other organic geochemical and(or) optical indicators show that the organic matter is mainly of marine origin. A difference of about 5‰ in σ 13 C between modern and Cretaceous OC-rich marine strata suggests either that Cretaceous marine planktonic organic matter had the same isotopic signature as modern marine plankton and that signature has been changed by diagenesis, or that OC derived from Cretaceous marine plankton was isotopically lighter by about 5‰ relative to modern plankton OC. Diagenesis does not produce a significant shift in σ 13 C in Miocene to Holocene sediments, and therefore probably did not produce the isotopically light Cretaceous OC. This means that Cretaceous marine plankton must have had σ 13 C values that were about 5‰ lighter than modern marine plankton, and at least several per mil lighter than Cretaceous terrestrial vegetation. The reason for these lighter values, however, is not obvious. It has been proposed that concentrations of CO 2 were higher during the middle Cretaceous, and this more available CO 2 may be responsible for the lighter σ 13 C values of Cretaceous marine organic matter.

Marine Geology↗

Paragenesis and chemistry of multistage tourmaline formation in the sullivan Pb-Zn-Ag deposit, British Columbia

Detailed petrographic study, scanning electron microscope imaging, and electron microprobe analyses of tourmalines from the Sullivan Pb-Zn-Ag massive sulfide deposit (British Columbia, Canada) document multiple paragenetic stages and large compositional variations. The tourmalines mainly belong to two common solid-solution series: dravite-schorl and dravite-uvite. Ca- and Fe-rich feruvite and alkali-deficient tourmalines are present locally. Products of tourmaline-forming stages include (from oldest to youngest): (1) rare Fe-rich dravite-schorl within black tourmalinite clasts in footwall fragmental rocks; (2) widespread Mg-rich, very fine grained, felted dravite in the footwall (the main type of tourmaline in the footwall tourmalinite pipe); (3) recrystallized, Fe-rich dravite-schorl (locally Ca-Fe feruvite) in the tourmalinite pipe, which preferentially occurs near postore gabbroic intrusions; (4) Mg-rich dravite or uvite associated with chlorite-pyrrhotite and chlorite-albite-pyrite-altered rocks in the shallow footwall and hanging wall; (5) discrete Mg-rich tourmaline grains associated with chlorite and discordant Mg-rich tourmaline rims which occur on disseminated Fe-rich schorl in the bedded Pb-Zn-Ag ores. The timing of rare Fe-rich schorl in the bedded ores is uncertain, but it most likely occurred during or between stages 2 and 3. The different paragenetic stages and their respective tourmaline compositions are interpreted in terms of a multistage evolution involving contributions from: (1) variable mixtures of synsedimentary, Fe-rich hydrothermal fluids and entrained seawater; (2) postore, Fe-rich, gabbro-related hydrothermal fluids; and (3) postore metamorphic reactions. Early synsedimentary, Fe-rich hydrothermal fluids which contained little or no entrained seawater formed Fe-rich black tourmalinite clasts locally in the footwall. The major type of tourmaline in the footwall tourmalinite pipe is Mg rich, recording seawater entrainment under high water/rock conditions, rather than control by the chemical composition of the original host sediments. Rare Fe-rich schorl within the bedded Pb-Zn-Ag ores is believed to have formed on the sea floor by reaction of an Fe-rich brine pool with detrital aluminous sediments. Postore emplacement of gabbro sills and local dikes in the footwall produced Fe-rich hydrothermal fluids, which were responsible for formation of minor Fe-rich dravite-schorl which overprinted earlier dravite. Postore, but synsedimentary, hydrothermal alteration involving entrained seawater was responsible for deposition of dravite and uvite in the hanging wall and for dravite in the brown tourmalinites of the shallow footwall. Mg-rich dravite-uvite associated with chlorite and in discordant rims on schorl in the bedded ores formed by sulfide-silicate reactions during greenschist facies regional metamorphism.

Economic Geology↗

A comparison of synthetic flowpaths derived from light detection and ranging topobathymetric data and National Hydrography Dataset High Resolution Flowlines

Bathymetric and topobathymetric light detection and ranging (lidar) digital elevation models created for the Delaware River were provided to the National Geospatial Program and used to evaluate synthetic flowpath extraction from bathymetric/topobathymetric lidar survey data as a data source for improving the density, distribution, and connectivity of the National Hydrography Dataset High Resolution Flowline Network. As the surface-water component of The National Map, the National Hydrography Dataset maintains the Nation’s drainage network flow information and geometries for surface-water features used in hydrologic, hydraulic, and other science and engineering disciplines. The regional lidar survey for the Delaware River between Hancock, New York, and Trenton, New Jersey, was collected for the U.S. Geological Survey using the Experimental Advanced Airborne Research Lidar sensor system and processed by the Coastal National Elevation Database Applications Program. Using 1 percent of the maximum flow accumulation value for the surveyed Delaware River corridor as the flow accumulation threshold for grid cells at 1-, 5-, and 10-meter resolution created 223 to 283 kilometers of synthetic flowpaths potentially representing the river channel thalweg, which is the deepest point in a riverbed cross-section. There was potential for improving the High Resolution National Hydrography Dataset (HR NHD) Flowline network in places where the Delaware River channel, depicted as an Artificial Path in the HR NHD, is offset from the extracted synthetic river flowpath which sometimes appeared better positioned than the Artificial Path to represent the river thalweg. For the same area, using 0.05 percent of the maximum flow accumulation at the 1-, 5-, and 10-meter resolutions extracted 744 to 1,317 kilometers of synthetic flowpaths, with extracted synthetic flowpaths representing the main river channel and additional synthetic flowpaths representing tributaries or streams adjacent to the main channel. Overlaying these results with the HR NHDFlowline Network indicates that some of the additional synthetic flowpaths are connected to or extend HR NHD stream/river feature types. Some disconnected or isolated synthetic flowpaths not included in stream/river feature types were validated in orthoimagery and U.S. Topo Maps and provide examples of how extracted synthetic flowpaths could be used to delineate new stream/river features. Other additional extracted synthetic flowpaths depict linear features such as canals, tree lines, roads, or linear topographic depressions. For some river reaches where obstructions to flow or where low-relief topographic or bathymetric surfaces alter the flow direction, the software tool used to develop the flow direction grid did not calculate a primary flowpath for the river channel. Based on the results of this analysis, site conditions for the Delaware River corridor did not affect the quality of lidar bathymetric survey data. However, depending on the resolution of the lidar bathymetric digital elevation models (BDEMs), site conditions do have different effects on results for extracted synthetic flowpaths. We found that synthetic flowpaths extracted from 1-meter resolution lidar DEMs had more varied flow directions around in-channel landforms that obstructed flow than synthetic flowpaths extracted from 5- or 10-meter resolution lidar DEMs. As a result the 1-meter resolution DEM created some isolated or discontinuous synthetic flowpath segments where the 5- and 10-meter DEMs developed more continuous flowpaths. In this case the river bed upstream from the in-channel obstruction is shallower than the river bed downstream. Under these conditions the 1-meter resolution DEM provided synthetic flowpaths delineating a potential river thalweg. In this same area, the software solution modified (virtually raised) the river bed in the 5- and 10-meter resolution DEMs and flattened the bathymetric surface to create a continuous downstream flow direction, which caused trellis-patterned synthetic flowpaths to form. Under different site conditions and converse to the above development of synthetic flowpaths at different resolutions, at an abandoned river flood plain (terrace) with low relief that is adjacent to the river channel, the flow direction grid for the 1-meter resolution DEM developed continuous synthetic flowpath corresponding to a HR NHD Flowline network stream/river feature that connected to the main river channel but the larger resolution DEMs created isolated or disconnected synthetic flowpaths. A project to continue an evaluation of benefits of or issues caused by extracting synthetic flowpaths to enhance the HR NHD could include a study to assess the potential for merging surface-water flowpaths extracted from lidar topobathymetry and 3D Elevation Program digital elevation models. The merged DEM approach to synthetic flowpath extraction could extend the HR NHDFlowline network and enhance flow accumulations that might develop better flow direction grids in low-relief areas. Because of the confined lateral extent of the Delaware River, the lidar DEMs were not used to create catchments or watersheds; however, the merged DEM approach could also be tested as a resource for enhancing HR NHD catchments and watersheds. This lidar DEM synthetic flowpath extraction project supports the National Geospatial Program efforts to collect and produce high-quality lidar data to provide 3-dimensional representations of natural feature and aligns with the National Spatial Data Infrastructure to improve utilization of geospatial data. The results also can be useful for understanding strategies that can help maintain quality data in the HR NHD programs. KEYWORDS: bathymetric, digital elevation model, extracted synthetic flowpath, lidar, High Resolution National Hydrography Dataset, topobathymetric

New Jersey↗