Geology Reports⌕ Search

SEARCH · Geology Reports

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,099 records · Page 61Linked to original sources

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↗

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↗

Groundwater and surface-water interaction within the upper Smith River Watershed, Montana 2006-2010

The 125-mile long Smith River, a tributary of the Missouri River, is highly valued as an agricultural resource and for its many recreational uses. During a drought starting in about 1999, streamflow was insufficient to meet all of the irrigation demands, much less maintain streamflow needed for boating and viable fish habitat. In 2006, the U.S. Geological Survey, in cooperation with the Meagher County Conservation District, initiated a multi-year hydrologic investigation of the Smith River watershed. This investigation was designed to increase understanding of the water resources of the upper Smith River watershed and develop a detailed description of groundwater and surface-water interactions. A combination of methods, including miscellaneous and continuous groundwater-level, stream-stage, water-temperature, and streamflow monitoring was used to assess the hydrologic system and the spatial and temporal variability of groundwater and surface-water interactions. Collectively, data are in agreement and show: (1) the hydraulic connectedness of groundwater and surface water, (2) the presence of both losing and gaining stream reaches, (3) dynamic changes in direction and magnitude of water flow between the stream and groundwater with time, (4) the effects of local flood irrigation on groundwater levels and gradients in the watershed, and (5) evidence and timing of irrigation return flows to area streams. Groundwater flow within the alluvium and older (Tertiary) basin-fill sediments generally followed land-surface topography from the uplands to the axis of alluvial valleys of the Smith River and its tributaries. Groundwater levels were typically highest in the monitoring wells located within and adjacent to streams in late spring or early summer, likely affected by recharge from snowmelt and local precipitation, leakage from losing streams and canals, and recharge from local flood irrigation. The effects of flood irrigation resulted in increased hydraulic gradients (increased groundwater levels relative to stream stage) or even reversed gradient direction at several monitoring sites coincident with the onset of nearby flood irrigation. Groundwater-level declines in mid-summer were due to groundwater withdrawals and reduced recharge from decreased precipitation, increased evapotranspiration, and reduced leakage in some area streams during periods of low flow. Groundwater levels typically rebounded in late summer, a result of decreased evapotranspiration, decreased groundwater use for irrigation, increased flow in losing streams, and the onset of late-season flood irrigation at some sites. The effect of groundwater and surface-water interactions is most apparent along the North and South Forks of the Smith River where the magnitude of streamflow losses and gains can be greater than the magnitude of flow within the stream. Net gains consistently occurred over the lower 15 miles of the South Fork Smith River. A monitoring site near the mouth of the South Fork Smith River gained (flow from the groundwater to the stream) during all seasons, with head gradients towards the stream. Two upstream sites on the South Fork Smith River exhibited variable conditions that ranged from gaining during the spring, losing (flowing from the stream to the groundwater) during most of the summer as groundwater levels declined, and then approached or returned to gaining conditions in late summer. Parts of the South Fork Smith River became dry during periods of losing conditions, thus classifying this tributary as intermittent. The North Fork Smith River is highly managed at times through reservoir releases. The North Fork Smith River was perennial throughout the study period although irrigation diversions removed a large percentage of streamflow at times and losing conditions persisted along a lower reach. The lowermost reach of the North Fork Smith River near its mouth transitioned from a losing reach to a gaining reach throughout the study period. Groundwater and surface-water interactions occur downstream from the confluence of the North and South Fork Smith Rivers, but are less discernible compared to the overall magnitude of the main-stem streamflow. The Smith River was perennial throughout the study. Monitoring sites along the Smith River generally displayed small head gradients between the stream and the groundwater, while one site consistently showed strongly gaining conditions. Synoptic streamflow measurements during periods of limited irrigation diversion in 2007 and 2008 consistently showed gains over the upper 41.4 river miles of the main stem Smith River where net gains ranged from 13.0 to 28.9 cubic feet per second. Continuous streamflow data indicated net groundwater discharge and small-scale tributary inflow contributions of around 25 cubic feet per second along the upper 10-mile reach of the Smith River for most of the 2010 record. A period of intense irrigation withdrawal during the last two weeks in May was followed by a period (early June 2010 to mid-July 2010) with the largest net increase (an average of 71.1 cubic feet per second) in streamflow along this reach of the Smith River. This observation is likely due to increased groundwater discharge to the Smith River resulting from irrigation return flow. By late July, the apparent effects of return flows receded, and the net increase in streamflow returned to about 25 cubic feet per second. Two-dimensional heat and solute transport VS2DH models representing selected stream cross sections were used to constrain the hydraulic properties of the Quaternary alluvium and estimate temporal water-flux values through model boundaries. Hydraulic conductivity of the Quaternary alluvium of the modeled sections ranged from 3x10-6 to 4x10-5 feet per second. The models showed reasonable approximations of the streambed and shallow aquifer environment, and the dynamic changes in water flux between the stream and the groundwater through different model boundaries.

Montana↗

Stable isotope study of water-rock interaction and ore formation, Bayhorse base and precious metal district, Idaho

The Bayhorse base and precious metal district is situated east of the Idaho batholith in south-central Idaho. The ores occur near the Nevada Mountain granitic stock as veins cutting the lower Paleozoic Ramshorn Slate and the Garden Creek Phyllite, and as fillings around breccia fragments within the Bayhorse Dolomite. The veins are dominated by siderite and tetrahedrite, with lesser quartz and galena, whereas the breccia ores dominantly comprise only quartz and galena. Mineralization and intrusive activity were contemporaneous during Cretaceous time. Fluid inclusion and stable isotope data indicate that mineralization formed from hot (ca. 375 degrees -225 degrees C), CO 2 -rich (< or =8.3+ or -1.4 mole %) brines (5-20 wt % NaCl equiv) at confining pressures between 1.1 and 1.7 kbars. Fluid cooling and the resulting CO 2 effervescence were the most important causes of ore deposition.Whole-rock delta 18 O and delta D values from the Garden Creek Phyllite define an isotopically depleted zone (60 km 2 ) around the Nevada Mountain stock and are the result of high-temperature interactions with ancient meteoric waters at water/rock ratios (mass fractions) ranging from 0.002 to 0.09. The delta 18 O (sub H 2 O) and delta D (sub H 2 O) values of the ore fluids ranged from 11.2 to 3.9 per rail and from -55 to -146 per mil, respectively. Comparison of the ore fluid delta 18 O (sub H 2 O) and delta D (sub H 2 O) values with hypothetical waters equilibrated with the Garden Creek Phyllite intricates that the hydrothermal fluids must have also interacted with the basal dolomite of Bayhorse Creek, which underlies the phyllite. Early delta 13 C (sub CO 2 ) values calculated in equilibrium with siderite ranged from -9.1 to -7.4 per mil, whereas later delta 13 C (sub CO 2 ) values obtained from quartz-hosted fluid inclusion extracts range from -4.0 to -1.8 per mil. The delta 13 C (sub CO 2 ) values for the hydrothermal fluids thus also record a transition from early water/rock interactions that were dominated by the Garden Creek Phyllite (organic matter delta 13 C = -23.6 to -20.9ppm) to later interactions that were influenced significantly by the basal dolomite of Bayhorse Creek (delta 13 C = -2.3 to -2.2ppm). The delta 34 S values of main-stage hydrothermal sulfides range from 3.0 to 28.3 per mil and display a distinct mode at approximately 10 per mil. The range of values may be interpreted as either a heterogeneous sedimentary source or mixed sedimentary-magmatic sources. On the basis of mass balance considerations, it is likely that all of the fluid salinity and the dissolved metals were derived from the Garden Creek Phyllite. The paragenetic sequence developed in the ores appears to be related to: (1) processes at the site of deposition, most notably cooling and CO 2 effervescence, and (2) water-rock interactions within the Garden Creek Phyllite in the 60-km 2 isotopically altered area surrounding the Nevada Mountain stock.

Economic Geology↗

Jurassic tectonics of northeastern Nevada and northwestern Utah from the perspective of barometric studies

Jurassic tectonism in the northeastern Great Basin produced varied structures, many closely associated with widespread magmatism at ca. 155–165 Ma and with local metamorphism. Many of the plutons are of suitable mineralogy for Al-in-hornblende barometry, providing the potential for depth data. We have studied conditions of metamorphism in the Pilot Range and barometry for six Jurassic plutons across the northeastern Great Basin. All barometry results are in harmony with pressures estimated from stratigraphic data, requiring little or no tectonic thickening. On the basis of structural styles and barometric data, we divide the northeastern Great Basin into three Jurassic tectonic provinces. An eastern extensional province, largely in western Utah, is characterized by Paleozoic strata that were thrust faulted and then intruded by shallow plutons shortly after or during normal and strike-slip faulting. Extension was probably a short-lived event associated with magmatism, but its west trend indicates a total reorientation of stress at this time, perhaps within transtensional strike-slip zones. A central province of modest, and possibly locally extreme, Jurassic shortening in eastern Nevada is characterized by metamorphosed Paleozoic rocks and by thrusts and kilometer-scale southeast-vergent folds. Upper amphibolite facies, but low pressure (3–4 kbar) metamorphism is present near Jurassic plutons in the Pilot Range and Ruby Mountains, probably indicating metamorphism induced by heat from magmas. In contrast, metamorphism in other ranges, which is known only to be pre–Late Cretaceous, indicates thickening of 10–20 km. This thickening may have entirely postdated the Jurassic. A western province in north-central Nevada is characterized by preserved Jurassic volcanic rocks and shallow plutons, indicating that little erosion, and probably surface uplift, occurred during the late Mesozoic. Folds and thrust faults indicate minor Jurassic shortening but many structures are undated. The low-pressure upper-crustal conditions for demonstrably Jurassic events suggest that higher-pressure metamorphism recorded in the central province is younger (Cretaceous) in age. We suggest that Jurassic structures were caused by distributed minor crustal shortening, manifested mainly as small-scale thrust faults. Local thermal highs created by plutonism produced metamorphic zones in relatively shallow crust. Shortening in the east was manifested by zones of strike-slip, within which plutons were emplaced in tensile niches. Lack of a deep foreland basin and lack of evidence for massive erosion argue against high-relief mountain belts caused by significant crustal shortening. Paleozoic rocks metamorphosed at pressures far in excess of stratigraphic burial are restricted to narrow lenses exhumed during Late Cretaceous and Tertiary extension and are bordered by rocks that always have been part of the shallow crust. The abundant shallow-crustal rocks preserved across the region indicate that a conventional hypothesis of large-scale, regional crustal thickening causing many kilometers of surface uplift and consequent erosion is unlikely to have taken place in the Mesozoic.

Nevada, Utah↗

Status of Water Levels and Selected Water-Quality Conditions in the Sparta-Memphis Aquifer in Arkansas and the Status of Water Levels in the Sparta Aquifer in Louisiana, Spring 2005

The U.S. Geological Survey in cooperation with the Arkansas Natural Resources Commission, the Arkansas Geological Commission, and the Louisiana Department of Transportation and Development has monitored water levels in the Sparta Sand of Claiborne Group and Memphis Sand of Claiborne Group since the 1920's. Ground-water withdrawals have increased while water levels have declined since monitoring was initiated. This report has been produced to describe ground-water levels in the aquifers in the Sparta Sand and Memphis Sand and provide information for the management of this valuable resource. The 2005 potentiometric-surface map of the aquifers in the Sparta Sand and Memphis Sand was constructed using water-level data collected in 333 wells in Arkansas and 120 wells in Louisiana during the spring of 2005. The highest water-level altitude measured in Arkansas was 327 feet above National Geodetic Vertical Datum of 1929 located in Grant County in the outcrop at the western boundary of the study area; the lowest water-level altitude was 189 feet below National Geodetic Vertical Datum of 1929 in Union County. The highest water-level altitude measured in Louisiana was 246 feet above National Geodetic Vertical Datum of 1929 located in Bossier Parish in the outcrop area near the western boundary of the study area; the lowest water-level altitude was 226 feet below National Geodetic Vertical Datum of 1929 in central Ouachita Parish. Three large depressions centered in Columbia, Jefferson, and Union Counties in Arkansas are the result of large withdrawals for industrial and public supplies. In Louisiana, three major pumping centers are in Ouachita, Jackson, and Lincoln Parishes. Water withdrawals from these major pumping centers primarily is used for industrial and public-supply purposes. Withdrawals from Ouachita and Lincoln Parishes and Union County, Arkansas, primarily for industrial purposes, have caused the resulting cones of depression to coalesce so that the -40 foot potentiometric contour encircles the three pumping centers. Seven smaller depressions are evident on the 2005 Sparta-Memphis potentiometric-surface map located in Webster and Winn Parishes, Louisiana, and Calhoun, Cleveland, western Columbia, Desha, and Lafayette Counties, Arkansas. The depression in Calhoun County initially was shown in the 1996-1997 potentiometric surface. The depression in Desha County initially was shown in the 1999 potentiometric surface. The depressions in Webster and Winn Parishes were shown as early as 1975. The depressions in Cleveland, western Columbia, and Lafayette Counties initially were shown in the 2003 potentiometric surface. A map of differences in water-level measurements between 2001 and 2005 was constructed using the difference between water-level measurements from 294 wells in Arkansas and 29 wells in Louisiana. The difference in water levels between 2001 and 2005 ranged from -30.1 to 44.6 feet. The largest rise of 44.6 feet in water level measured was in Union County in Arkansas. The largest decline of 30.1 feet in water level measured was in Columbia County in Arkansas. Areas with a general rise in water levels in Arkansas are shown in Arkansas, Columbia, Craighead, Jefferson, Prairie, and the western half of Union Counties. The area around west-central Union County had rises as much as 44.6 feet, with seven wells showing a rise of 20 feet or greater, which is an annual rise of 5 feet or greater. Areas in Arkansas with a general decline in water level are shown in western Bradley, eastern Calhoun, Cleveland, Cross, Desha, Drew, Lafayette, Lee, Lincoln, Lonoke, Poinsett, and the eastern half of Union Counties. In Louisiana, the water-level difference map showed a general rise in water levels in northern Claiborne, northern Webster, and northwestern Union Parishes mainly because of a decrease in industrial withdrawals in southern Arkansas, particularly Union County. Another rise in water level was indicated in western

Scientific Investigations Report↗