SEARCH · Geology Reports

Results for “American Geology”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,513 records · Page 84Linked to original sources

USGS2020

Testing reproducibility of vitrinite and solid bitumen reflectance measurements in North American unconventional source-rock reservoir petroleum systems

An interlaboratory study (ILS) was conducted to test reproducibility of vitrinite and solid bitumen reflectance measurements in six mudrock samples from United States unconventional source-rock reservoir petroleum systems. Samples selected from the Marcellus, Haynesville, Eagle Ford, Barnett, Bakken and Woodford are representative of resource plays currently under exploitation in North America. All samples are from marine depositional environments, are thermally mature (T max >445 °C) and have moderate to high organic matter content (2.9–11.6 wt% TOC). Their organic matter is dominated by solid bitumen, which contains intraparticle nano-porosity. Visual evaluation of organic nano-porosity (pore sizes < 100 nm) via SEM suggests that intraparticle organic nano-pores are most abundant in dry gas maturity samples and less abundant at lower wet gas/condensate and peak oil maturities. Samples were distributed to ILS participants in forty laboratories in the Americas, Europe, Africa and Australia; thirty-seven independent sets of results were received. Mean vitrinite reflectance (VR o ) values from all ILS participants range from 0.90 to 1.83% whereas mean solid bitumen reflectance (BR o ) values range from 0.85 to 2.04% (no outlying values excluded), confirming the thermally mature nature of all six samples. Using multiple statistical approaches to eliminate outlying values, we evaluated reproducibility limit R, the maximum difference between valid mean reflectance results obtained on the same sample by different operators in different laboratories using different instruments. Removal of outlying values where the individual signed multiple of standard deviation was >1.0 produced lowest R values, generally ≤0.5% (absolute reflectance), similar to a prior ILS for similar samples. Other traditional approaches to outlier removal (outside mean ± 1.5*interquartile range and outside F10 to F90 percentile range) also produced similar R values. Standard deviation values < 0.15*(VR o or BR o ) reduce R and should be a requirement of dispersed organic matter reflectance analysis. After outlier removal, R values were 0.1%–0.2% for peak oil thermal maturity, about 0.3% for wet gas/condensate maturity and 0.4%–0.5% for dry gas maturity. That is, these R values represent the uncertainty (in absolute reflectance) that users of vitrinite and solid bitumen reflectance data should assign to any one individual reported mean reflectance value from a similar thermal maturity mudrock sample. R values of this magnitude indicate a need for further standardization of reflectance measurement of dispersed organic matter. Furthermore, these R values quantify realistic interlaboratory measurement dispersion for a difficult but critically important analytical technique necessary for thermal maturity determination in the source-rock reservoirs of unconventional petroleum systems.

Marine and Petroleum Geology↗

USGS2007

Origin of the Lyme Dome and implications for the timing of multiple Alleghanian deformational and intrusive events in southern Connecticut

Geologic mapping, structural analysis, and geochronology in the area of the Lyme dome, southern Connecticut provides constraints on the origin of the rocks in the core of the dome, the absolute timing of the principal deformational and thermal events attributed to Alleghanian orogenesis, and the processes that generated the dome. Detrital zircon geochronology in combination with ages on intrusive rocks brackets the deposition of quartzite in the core of the dome sometime between ca. 925 and 620 Ma. Granite and granodiorite intruded the Neoproteorozic metasedimentary rocks in the core of the dome at ca. 620 to 610 Ma. Four major early Permian events associated with the Alleghanian orogeny affected the rocks in the Lyme dome area. Syn-tectonic migmatization and widespread penetrative deformation (D1, ca. 300 - 290 Ma) included emplacement of alaskite at 290 ?? 4 Ma during regional foliation development and aluminosilicate-orthoclase metamorphic conditions. Rocks of the Avalon terrane may have wedged between Gander cover rocks and Gander basement in the core of the Lyme during D1. Limited structural evidence for diapiric uplift of the Lyme dome indicates that diapirism started late in D1 and was completed by D2 (ca. 290 - 280 Ma) when horizontal WNW contractional stresses dominated over vertical stresses. Second sillimanite metamorphism continued and syn-tectonic D2 granite pegmatite (288 ?? 4 Ma) and the Joshua Rock Granite Gniess (284 ?? 3 Ma) intruded at this time. North-northwest extension during D3 (ca. 280 - 275 Ma) led to granitic pegmatite intrusion along S3 cleavage planes and in extensional zones in boudin necks during hydraulic failure and decompression melting. Intrusion of a Westerly Granite dike at 275 ?? 4 Ma suggests that D3 extension was active, and perhaps concluding, by ca. 275 Ma. Late randomly oriented but gently dipping pegmatite dikes record a final stage of intrusion during D4 (ca. 275 - 260 Ma), and a switch from NNW extension to vertical unloading and exhumation. Monazite and metamorphic zircon rim ages record this event at ca. 259 Ma. The evolution of the Lyme dome involved D1 mylonitization, intrusion, and migmatization during north-directed contraction, limited late D1 diapirism, D2 migmatization during WNW contraction with associated flexural flow and fold interference, D3 NNW horizontal extension and decompression melting, and final D4 vertical extension and rapid exhumation. Late regional uplift, extension, and normal faulting at higher crustal levels may have been caused by diapiric rise of the lower crust, below the structural level of the Lyme dome. The rocks record no evidence of Acadian metamorphism or deformation, suggesting that the Gander zone here was not tectonically juxtaposed with Avalon until the Alleghanian orogeny.

American Journal of Science↗

USGS2000

Late Holocene Peléan-style eruption at Tacaná volcano, Mexico and Guatemala: past, present, and future hazards

Tacaná volcano, located on the border between Mexico and Guatemala, marks the northern extent of the Central American volcanic chain. Composed of three volcanic structures, it is a volcanic complex that has had periodic explosive eruptions for at least the past 40 k.y. The most recent major eruption occurred at the San Antonio volcano, the youngest volcanic edifice forming the complex, about 1950 yr ago. The Peléan style eruption, issued from the southwest part of the dome, and swept a 30° sector with a hot block and ash flow that traveled about 14 km along the Cahoacán ravine. Deposits from this event are well exposed around the town of Mixcun and were therefore given the name of that town, the Mixcun flow deposit. The Mixcun flow deposit is, in the channel facies, a light gray, massive, thick (>10 m), matrix-supported unit with dispersed lithic clasts of gravel to boulder size, divisible in some sections into a variable number of flow units. The overbank facies is represented by a thin (<1 m), massive, matrix-supported unit. In both of these facies the deposit has disseminated charcoal, fumarolic pipes, and juvenile lithics with cooling joints. The Mixcun flow deposit contains clasts of (1) light gray, dense andesite, (2) dark gray, glassy and banded andesite, and (3) minor altered red andesite from the edifice, set in a matrix of sand and silt. The Mixcun flow deposit covers an area of at least 25 km 2 and has a minimum estimated volume of 0.12 km 3 . Basaltic-andesite inclusions (54% SiO 2 ) and various signs of disequilibrium in the mineral assemblage of the two-pyroxene andesitic products (60%–63% SiO 2 ) suggest that magma mixing may have triggered the eruption. Following deposition of the Mixcun flow deposit andesitic to dacitic (62%–64% SiO 2 ) lava flows were extruded and a dacitic dome (64.4% SiO 2 ) at the San Antonio summit formed. Syn-eruptive and posteruptive lahars flooded the main drainages of the Cahoacán and Izapa-Mixcun valleys in the area of the present city of Tapachula (population 250000) and the pre-Hispanic center of Izapa. Three radiocarbon ages date this event between A.D. 25 and 72 (range ±1σ, 38 B.C.–A.D. 216), which correlates with a halt in construction at Izapa (Hato phase of ca. 50 B.C.–A.D. 100), probably due to temporary abandonment of the city caused by lahars. Another similar event would produce extensive damage to the towns (population of about 68,000 people) now built upon the Mixcun flow deposit. The main summit of Tacaná volcano continues to show signs of fumarolic activity; the most recent period of activity in 1985–1986 culminated in a minor phreatic explosion.

Tapachula↗

USGS1993

The Geysers-Clear Lake area, California: Thermal waters, mineralization, volcanism, and geothermal potential

Manifestations of a major thermal anomaly in the Geysers-Clear Lake area of northern California include the late Pliocene to Holocene Clear Lake Volcanics, The Geysers geothermal field, abundant thermal springs, and epithermal mercury and gold mineralization. The epithermal mineralization and thermal springs typically occur along high-angle faults within the broad San Andreas transform fault system that forms the western boundary of the North American plate in this area. The young volcanic rocks overlie Mesozoic marine rocks of the Great Valley sequence which have been thrust above the coeval Franciscan Complex and penecontemporaneously dropped back down along low-angle detachment faults.Many of the waters of the region are nonmeteoric as defined by their isotopic signature. One type of isotopically shifted water emerges from or near Great Valley sequence rocks and is the most chloride rich. It is interpreted to be evolved connate in origin. A second type, evolved meteoric water, has moderate chloride contents, high boron contents, and high B/Cl ratios and is found locally in Franciscan rocks, notably at the Sulphur Bank mercury mine, where it probably results from near-closed-system, repeated boiling of meteoric water in host rocks that also contribute organic components to the water. At the Sulphur Bank mine fracturing of otherwise impermeable Franciscan rocks by faulting has created a localized zone of permeability in which thermal water boils repeatedly with limited venting to the surface. Boron-rich fluids were apparently present at depth in The Geysers geothermal field when intrusion of silicic magma occurred, because the concealed intrusion of felsite is surrounded by a halo of tourmaline-bearing hornfels. The volume of this poorly dated early to middle Quaternary intrusive body probably exceeds the 100 km 3 of erupted Clear Lake Volcanics. Similar intrusions may have occurred in the eastern part of the area at Wilbur Springs and the McLaughlin mine, where gold deposition and evidence of hydrothermal phenomena suggest more magmatic activity than is indicated by the small exposed bodies of early Quaternary basaltic lava. The Clear Lake Volcanics are the present locus of volcanism in the northern Coast Ranges and progressively older volcanic centers occur to the south. Geophysical data suggest that a large silicic magma body may be centered north of The Geysers steam field, providing the heat for the geothermal field.Geothermal power production has peaked at The Geysers and pressure declines indicate significant depletion of the fluid resource. The vapor-dominated field evolved from a preexisting hydrothermal system within fractured, otherwise impermeable Franciscan metamorphic rocks. A deep water table of saline fluid has been postulated to be present under the steam field, but no chloride-rich water has been found at drillable depth. We propose that recently discovered, isotopically shifted steam in the northwest Geysers area indicates the presence not of deep connate water but xather of boiled-down, boron-rich Franciscan evolved meteoric water. This water is likely to be present in limited quantities and will not provide a significant hot water resource for geothermal power production at The Geysers field or from the main Clear Lake volcanic field.

Economic Geology↗

USGS1984

Geologic evolution, sedimentation, and paleoenvironments of the Angola Basin and adjacent Walvis Ridge: Synthesis of results of Deep Sea Drilling Project Leg 75

The section recovered at Site 530 (Holes 53OA and 530B) consists of eight sedimentary units and one basalt unit. The composition of the basalt recovered in Hole 53OA is distinct from typical mid-ocean ridge basalts (MORBs) but is similar to that of Hawaiian tholeiites and basalt from the central part of Walvis Ridge. Throughout most of its history, the southern Angola Basin received large volumes of redeposited material in the form of turbidites and, most recently, debris-flow deposits. Most of this material was derived from Walvis Ridge to the south, but thickness trends of acoustic units suggest that some of the sediment was derived from the African continental margin to the east. The basal sedimentary unit (Albian to Santonian) at Site 530 contains 262 beds of black shale that are interbedded with green and red claystone. Black shale makes up less than 10% of the total section, but in two cores of early Turonian age, black shale beds compose about 50% of the section. The black shales contain up to 19% organic carbon (average of about 5%) that is mainly of autochthonous marine origin but with significant contributions from terrigenous organic matter. The origin of these more- and less-reduced interbedded lithologies with varying amounts and types of organic matter, and variable amounts of pelagic, hemipelagic, and turbiditic sediment is complex and cannot be explained by any one simple process. Many factors affecting the concentration of dissolved oxygen in the bottom waters of the Angola Basin varied throughout the middle Cretaceous to produce bottom-water conditions that fluctuated between mildly oxic and oxygen-deficient, but most of the time bottom-waters and sediment-interstitial waters were sufficiently oxic to permit the accumulation of red oxidized sediment. A relatively complete sedimentary record of the Cretaceous/Tertiary boundary was recovered within a sequence of mudstone and marlstone turbidites in Hole 530A. There is a significant increase in the concentration of iridium above background levels at the boundary. High concentrations of many other elements also occur within the same stratigraphic interval as the iridium anomaly. Furthermore, there is a marked decrease in CaCO3 in the Tertiary strata above the iridium anomaly which suggests that the production of shallow-water carbonate also may have been affected by whatever caused elevated concentrations of iridium and other elements. These observations are consistent with the asteroid-impact theory proposed to explain the worldwide occurrence of an iridium anomaly at the Cretaceous/Tertiary boundary. The Cenozoic history of the Angola Basin was controlled mainly by (1) restriction of bottom-water flow from the south by Walvis Ridge; (2) development of glaciation on Antarctica; (3) opening of circulation passages in the southern oceans; (4) rapid turnover of cold, nutrient-rich waters that resulted in high productivity of diatoms; (5) influx of terrigenous sediment mainly by turbidity currents; and (6) production and preservation of carbonate sediment. The most distinctive Cenozoic event recorded in the section at Site 530 is the beginning of extensive glaciation on Anarctica and concomitant initiation of modern thermohaline bottom-water circulation that is manifested as a middle Eocene to middle Oligocene unconformity or compressed section accompanied by a drastic decrease in accumulation of CaCO3. Diatom abundances in HPC cores from Walvis Ridge (Site 532) and Angola Basin (Hole 53OB) indicate that Benguela upwelling in these areas began in the late Miocene, reached a peak in the late Pliocene to early Pleistocene, and declined thereafter. Short-term variations in sediment composition at Site 532 are manifested as cyclic variations in concentrations of clay, CaCO3, and organic carbon with average periodicities of about 30-60 k.y. The main variability that produced the cycles probably was the influx of terrigenous clastic material which diluted the CaCO3. The sediment at Site 532 also contains several percent organic carbon that is dominantly of marine origin, but with significant terrigenous components. Data from multichannel seismic, gravity, and magnetic surveys were used to define the regional stratigraphic and structural evolution of Walvis Ridge and adjacent Cape and Angola basins. Six structural provinces are recognized, four on Walvis Ridge and two additional provinces that correspond to the Cape and Angola basins. The two eastern structural provinces on Walvis Ridge are underlain by continental crust. The two western structural provinces are underlain by oceanic basement. Two main directions of faults are evident in seismic profiles, one trending N 10° and one trending N 60°. The N 60° trend corresponds to the general orientation of the northern and southern flanks of Walvis Ridge as well as to the dominant direction of fracture zones. During the first phase of separation of Africa from South America (ca. 120-130 m.y. ago), a voluminous mass of volcanics was emplaced simultaneous with the emplacement of basalt in the Parana Basin of Brazil and the Kaokoveld Region of South Africa. This period of volcanism also formed the series of seaward-dipping internal basement reflectors that are characteristic of the two structural provinces of Walvis Ridge. A system of fault blocks developed in the brittle upper part of the newly formed crust. During the second phase of rifting, which ended before late Aptian, more tilted fault blocks were created in the upper brittle stratified continental crust. Magnetic lineations in basement rocks in the Angola and Cape basins in the vicinity of Walvis Ridge are not distinct but suggest that oceanic crust began to be emplaced between 120 and 112 m.y. ago (Barremian to early Aptian). At least part of the oceanic crust of the central plateau of eastern Walvis Ridge (structural province 3) may have been emplaced before any oceanic crust formed in the adjacent basins. A ridge jump occurred during the late Aptian to early Albian in the southern part of the Angola Basin which translated the previously formed oceanic crust and its overlying evaporite deposits on the South American side. Several ridge jumps occurred on both sides of Walvis Ridge during the Late Cretaceous and early Tertiary to produce a 500-km-long segment of mid-ocean ridge.

Initial Reports of the D.S.D.P.↗

USGS2013

Natural-color and color-infrared image mosaics of the Colorado River corridor in Arizona derived from the May 2009 airborne image collection

The Grand Canyon Monitoring and Research Center (GCMRC) of the U.S. Geological Survey (USGS) periodically collects airborne image data for the Colorado River corridor within Arizona (fig. 1) to allow scientists to study the impacts of Glen Canyon Dam water release on the corridor’s natural and cultural resources. These data are collected from just above Glen Canyon Dam (in Lake Powell) down to the entrance of Lake Mead, for a total distance of 450 kilometers (km) and within a 500-meter (m) swath centered on the river’s mainstem and its seven main tributaries (fig. 1). The most recent airborne data collection in 2009 acquired image data in four wavelength bands (blue, green, red, and near infrared) at a spatial resolution of 20 centimeters (cm). The image collection used the latest model of the Leica ADS40 airborne digital sensor (the SH52), which uses a single optic for all four bands and collects and stores band radiance in 12-bits. Davis (2012) reported on the performance of the SH52 sensor and on the processing steps required to produce the nearly flawless four-band image mosaic (sectioned into map tiles) for the river corridor. The final image mosaic has a total of only 3 km of surface defects in addition to some areas of cloud shadow because of persistent inclement weather during data collection. The 2009 four-band image mosaic is perhaps the best image dataset that exists for the entire Arizona part of the Colorado River. Some analyses of these image mosaics do not require the full 12-bit dynamic range or all four bands of the calibrated image database, in which atmospheric scattering (or haze) had not been removed from the four bands. To provide scientists and the general public with image products that are more useful for visual interpretation, the 12-bit image data were converted to 8-bit natural-color and color-infrared images, which also removed atmospheric scattering within each wavelength-band image. The conversion required an evaluation of the histograms of each band’s digital-number population within each map tile throughout the corridor and the determination of the digital numbers corresponding to the lower and upper one percent of the picture-element population within each map tile. Visual examination of the image tiles that were given a 1-percent stretch (whereby the lower 1- percent 12-bit digital number is assigned an 8-bit value of zero and the upper 1-percent 12-bit digital number is assigned an 8-bit value of 255) indicated that this stretch sufficiently removed atmospheric scattering, which provided improved image clarity and true natural colors for all surface materials. The lower and upper 1-percent, 12-bit digital numbers for each wavelength-band image in the image tiles exhibit erratic variations along the river corridor; the variations exhibited similar trends in both the lower and upper 1-percent digital numbers for all four wavelength-band images (figs. 2–5). The erratic variations are attributed to (1) daily variations in atmospheric water-vapor content due to monsoonal storms, (2) variations in channel water color due to variable sediment input from tributaries, and (3) variations in the amount of topographic shadows within each image tile, in which reflectance is dominated by atmospheric scattering. To make the surface colors of the stretched, 8-bit images consistent among adjacent image tiles, it was necessary to average both the lower and upper 1-percent digital values for each wavelength-band image over 20 river miles to subdue the erratic variations. The average lower and upper 1-percent digital numbers for each image tile (figs. 2–5) were used to convert the 12-bit image values to 8-bit values and the resulting 8-bit four-band images were stored as natural-color (red, green, and blue wavelength bands) and color-infrared (near-infrared, red, and green wavelength bands) images in embedded geotiff format, which can be read and used by most geographic information system (GIS) and image-processing software. The tiff world files (tfw) are provided, even though they are generally not needed for most software to read an embedded geotiff image. All image data are projected in the State Plane (SP) map projection using the central Arizona zone (202) and the North American Datum of 1983 (NAD83). The map-tile scheme used to segment the corridor image mosaic followed the standard USGS quarter-quadrangle (QQ) map borders, but the high resolution (20 cm) of the images required further quarter segmentation (QQQ) of the standard QQ tiles, where the image mosaic covered a large fraction of a QQ map tile (segmentation shown in (figure 6), where QQ_1 to QQ_4 shows the number convention used to designate a quarter of a QQ tile). To minimize the size of each image tile, each image or map tile was subset to only include that part of the tile that had image data. In addition, some QQQ image tiles within a QQ tile were combined when adjacent QQQ map tiles were small. Thus, some image tiles consist of combinations of QQQ map tiles, some consist of an entire QQ map tile, and some consist of two adjoining QQ map tiles. The final image tiles number 143, which is a large number of files to list on the Internet for both the natural-color and color-infrared images. Thus, the image tiles were placed in seven file folders based on the one-half-degree geographic boundaries within the study area (fig. 7). The map tiles in each file folder were compressed to minimize folder size for more efficient downloading. The file folders are sequentially referred to as zone 1 through zone 7, proceeding down river (fig. 7). The QQ designations of the image tiles contained within each folder or zone are shown on the index map for each respective zone (figs. 8–14).

Arizona↗

USGS1992

Impact origin of the Avak Structure, Arctic Alaska, and genesis of the Barrow gas fields

Geophysical and subsurface geologic data suggest that the Avak structure, which underlies the Arctic Coastal Plain 12 km southeast of Barrow, Alaska, is a hypervelocity meteorite or comet impact structure. The structure is a roughly circular area of uplifted, chaotically deformed Upper Triassic to Lower Cretaceous sedimentary rocks 8 km in diameter that is bounded by a ring of anastomosing, inwardly dipping, listric normal faults 12 km in diameter. A zone of gently outward-dipping sedimentary country rocks forms a discontinuous ring of "rim anticlines" within the peripheral ring of normal faults. Beyond these anticlines, the sedimentary rocks are almost flat-lying. Basement consists of strongly deformed Ordovician and Silurian argillite. Density and acoustic impedance con rasts between the argillite and the overlying strata produce gravity and seismic-reflection signatures that define a ring of anticlines around the disturbed zone and a structural high surrounded by an annular structural low at its center. In the adjacent Barrow gas fields, the tops of the informally named Neocomian "pebble shale" unit and the gas-producing Lower Jurassic Barrow sand (local usage) lie at average subsea depths of 488 m and 670 m, respectively. In the Avak 1 well, drilled on the central high, the pebble shale and the Barrow sand lie near the surface, documenting more than 500 m of relative uplift at the high. The cores in this well have steep dips (30-90 degrees), mixed breccia with Franklinian argillite clasts 10 and 90 m above basement, quartz grains with shock mosaicism and multiple sets of shock lamellae, oriented concussion fractures in sand-size quartz grains, and shatter cones resembling those found in the peripheral zones of well-documented impact structures. In addition, above-background levels o fractured quartz grains in Barrow sand were found as far as 19 km beyond the rim of the Avak structure. Data concerning the age of the Avak structure are not definitive. If submarine landslide deposits in the upper part of the Aptian and Albian Torok Formation, in the subsurface 200 km to the east, were triggered by the Avak event, then the Avak meteorite struck a submerged marine shelf about 100 + or - 5 Ma. However, the impact features found at Avak (shatter cones, concussion fractures, shock lamellae and shock mosaicism in quartz grains, and widespread cataclasis) characterize the distal zones of meteorite impact structures. Fused rocks, plastic deformation, and shock-metamorphic minerals found in more proximal zones of impact structures are apparently missing. These observations, and the lack of Avak ejecta in cuttings and cores from the Torok Formation and Nanushuk Group (Albian to middle Cenomanian) in surrounding test wells, indicate that the impact event postdated these beds. In this case, the Avak meteorite struck a Late Cretaceous or Tertiary marine shelf or coastal plain between the Cenomanian (ca. 95 Ma), and deposition of the basal beds of the overlying late Pliocene and Quaternary Gubik Formation (ca. 3 Ma).

Alaska↗

USGS1997

A geologic history of the north-central Appalachians, part 2: The Appalachian basin from the Silurian through the Carboniferous

The north-central Appalachians occupy a critical position within the 3000+ km-long Appalachian orogen, lying southwest of the boundary between the central and northern Appalachians (CNAB). The one-billion-year-long history of tectonic activity in the north-central Appalachians includes the assembly and breakup of a late Proterozoic supercontinent, the creation and evolution of the Appalachian orogen during the Paleozoic, and the Mesozoic transformation of the active orogen into a passive margin during Pangea's disassembly. An important part of the Middle and Late Paleozoic history is the formation and growth of the Appalachian basin, an enormous, elongate continental basin lying cratonward of the active Appalachian internides along Laurentia's eastern margin. The Appalachian basin developed out of the Taconic orogeny, the second of the four orogenies that formed and modified the character of the Appalachian orogen. Prior to the Taconic orogeny, the eastern margin of Laurentia consisted of a broad carbonate shelf facing the Theic Ocean. The Taconic orogeny obduced various Theic components, including microcontinents, magmatic arcs, and accretionary prisms, onto the carbonate shelf. The resulting Taconic highlands formed a topographic barrier between Theia and the craton's interior. Because tectonic activity in the internides continued more or less throughout the remainder of the Paleozoic, the Appalachian basin never had a direct connection with Theia. The coarse-grained molasse from the newly uplifted Taconic highland spread westward over most of the basin during the Early Silurian. Carbonate deposition gradually encroached eastward over the basin as the siliciclastic input from the southeast waned. A resurgence in the earliest Late Silurian, perhaps related to orogenesis in Newfoundland and the Maritimes, expanded the clastic wedge somewhat. Before long, however, carbonate deposition once again dominated most of the north-central basin for the remainder of the Silurian and into the Early Devonian. The Early-to-Middle Devonian Acadian orogeny began introducing siliciclastic material into the eastern part of the Appalachian basin, thereby ending the largely paralic environment that persisted from the Late Silurian. A number of deltas formed around local sediment-input centers during the Middle Devonian. Continued uplift and/or orogenesis must have continued into the Late Devonian because a vast amount of terrigenous sediment was introduced into the Appalachian basin to form the vast Upper Devonian Catskill delta. Although the principal tectonic activity appears to have centered in New England, the large volume of sediment input suggests that orogenesis may have extended southward to the north-central Appalachians. However, direct evidence of actual Acadian deformation and metamorphism is lacking at this latitude. The north-central Appalachian basin underwent a significant change during the Carboniferous. Prior to this time, the shoreline shifted laterally through transgressions and regressions largely as a function of sediment input and regional subsidence. In the Carboniferous, vertical movements in the basin became more important. During the Early Carboniferous, much of the Mauch Chunk delta was eroded in the central and western parts of the basin. At the end of the Early Carboniferous, an unconformity truncated progressively older rocks to the north toward New York. New sediment transport patterns and depositional environments within the basin during the Late Carboniferous probably reflect significant topographic and tectonic changes, not only in the hinterland, but in the basin itself, perhaps in conjunction with climatic changes. Extensive interfluvial swamps on a lower delta plain accumulated thick organic deposits which would become coal. This delta plain graded southwestward into tidal and marine environments. The Alleghany orogeny in the Early Permian interrupted deposition in much of the north-central Appalachian basin and profoundly altered its structure.

American Journal of Science↗

USGS2021

Paragenesis of an orogenic gold deposit: New insights on mineralizing processes at the Grass Valley District, California

The Grass Valley orogenic gold district in the Sierra Nevada foothills province, central California, is the largest historical gold producer of the North American Cordillera. Gold mineralization is associated with shallowly dipping north-south veins hosted by the 160 Ma Grass Valley granodiorite to the southwest of the Grass Valley fault and steeply dipping east-west veins in accreted oceanic rocks to the northeast of this major fault. Quartz veins from both vein types show well-preserved primary textural relationships. Using a combination of petrographic and microanalytical techniques, the paragenetic sequence of minerals within the veins and the compositions of ore minerals were determined to constrain the mechanisms of quartz vein formation and gold deposition. The veins are composed of early quartz that formed through cooling of hydrothermal fluids derived from a geopressured reservoir at depth. The early quartz shows growth zoning in optical cathodoluminescence and contains abundant growth bands of primary inclusions. The primary inclusion assemblages and myriads of crosscutting secondary fluid inclusions have been affected by postentrapment modification, suggesting that early quartz formation was postdated by pronounced pressure fluctuations. These pressure fluctuations, presumably involving changes from lithostatic to hydrostatic conditions, may be related to fault failure of the host structure as predicted by the fault-valve model. Fluid flow associated with pressure cycling took place along microfractures and grain boundaries resulting in extensive recrystallization of the early quartz. Deposition of pyrite, arsenopyrite, and first-generation gold from these hydrothermal fluids causing recrystallization of the early quartz occurred as a result of wall-rock sulfidation. The gold forms invisible gold in the compositionally zoned pyrite or micron-sized inclusions within pyrite growth zones. The latest growth zones in euhedral quartz crystals that formed in association with this stage of the paragenesis contain very rare primary fluid inclusions that have not been affected by postentrapment modification. The hydrothermal system transitioned entirely to hydrostatic conditions immediately after formation of the latest quartz, explaining the preservation of the primary fluid inclusions. The formation of minor quartz in open spaces was followed by the deposition of second-generation native gold and telluride minerals that are commonly associated with base metal sulfides. Ore formation at this stage of the paragenesis is attributed to the rapid decompression of hydrothermal fluids escaping from the geopressured part of the crust into the overlying hydrostatic realm. There is no fluid inclusion evidence that this pressure drop resulted in fluid immiscibility of the hydrothermal fluids. Fluid inclusion evidence suggests that the north-south veins formed at a paleodepth of ~8 km, whereas the east-west veins appear to have formed at ~10 to 11 km below surface, confirming previous inferences that the NE-dipping Grass Valley reverse fault accommodated a large displacement. The findings of the study at Grass Valley have significant implications for the model for orogenic gold deposits, as the reconstruction of the paragenetic relationships provides evidence for the occurrence of two discrete events of gold introduction that occurred at different conditions during the evolution of the hydrothermal system.

California↗

USGS2015

Application of U-Th-Pb phosphate geochronology to young orogenic gold deposits: New age constraints on the formation of the Grass Valley gold district, Sierra Foothills province, California

The Grass Valley orogenic gold district in the Sierra Nevada foothills province, central California, the largest historic gold producer of the North American Cordillera, comprises both steeply dipping east-west (E-W) veins located along lithologic contacts in accreted ca. 300 and 200 Ma oceanic rocks and shallowly dipping north-south (N-S) veins hosted by the Grass Valley granodiorite; the latter have yielded about 70 percent of the 13 million ounces of historic lode gold production in the district. The oceanic host rocks were accreted to the western margin of North America between 200 and 170 Ma, metamorphosed to greenschist and amphibolite facies, and uplifted between 175 and 160 Ma. Large-scale magmatism in the Sierra Nevada occurred between 170-140 Ma and 120-80 Ma, with the Grass Valley granodiorite being emplaced during the older episode of magmatism. Uranium-lead isotopic dating of hydrothermal xenotime yielded the first absolute age of 162±5 Ma for the economically more significant N-S veins. The vein-hosted xenotime, as well as associated monazite, are unequivocally of hydrothermal origin as indicated by textural and chemical characteristics, including grain shape, lack of truncated growth banding, lack of a Eu anomaly, and low U and Th concentrations. Furthermore, the crack-seal texture of the veins, with abundant wallrock slivers, suggests their formation as a result of episodic fluid flow possibly related to reoccurring seismic events, rather than a period of fluid exsolution from an evolving magma. The N-S veins are temporally distinct from a younger 153-151 Ma gold event that was previously reported for the E-W veins. Overlapping U-Pb zircon (159.9±2.2 Ma) and 40Ar/39Ar biotite and hornblende (159.7±0.6 to 161.9±1.4 Ma) ages and geothermobarometric calculations indicate that the Grass Valley granodiorite was emplaced at ca. 160 Ma at elevated temperatures (~800°C) within approximately 3 km of the paleosurface and rapidly cooled to the ambient temperature of the surrounding country rocks (<300°C). The age of the granodiorite is indistinguishable from that of the N-S veins, as recorded by the U-Pb age of xenotime in those veins. Consequently, the N-S veins must have formed between 162 and 157 Ma, the maximum permissive age of magma emplacement and the youngest permissive xenotime U-Pb age, respectively, during an E- to ENE-directed compressional regime. The geochemistry of the Grass Valley granodiorite is consistent with it being the product of arc magmatism. It served as a receptive host for mineralization, but it is has no direct genetic relationship to gold mineralization. Initial uplift of the intrusive mass correlates with the initial voluminous fluid flow event and vein formation at depths of no greater than 3 km. The E-W gold-bearing veins hosted within greenschist-facies country rocks adjacent to the intrusion formed during a second hydrothermal event 5-10 million years later than the magmatism and were contemporaneous with a shift to a transtensional deformation denoted by sinistral strike-slip faulting.

California↗

USGS2021

National Park Service Vegetation Mapping Inventory Program: Great Smoky Mountains National Park vegetation mapping project

The National Park Service (NPS) Vegetation Mapping Inventory (VMI) Program is an effort to classify, describe, and map existing vegetation communities in national park units throughout the United States. The NPS VMI Program is managed by the NPS Natural Resource Stewardship and Science Inventory and Monitoring Program and provides baseline vegetation information to natural resource managers, researchers, and ecologists. The U.S. Geological Survey Upper Midwest Environmental Sciences Center, NatureServe, and NPS Great Smoky Mountains National Park (GRSM, also referred to as the “Park”) have completed vegetation classification and mapping of GRSM, including the Foothills Parkway, for the NPS VMI Program. Mappers, ecologists, and botanists collaborated to affirm vegetation types of GRSM and to determine how best to map the vegetation types by using aerial imagery. A vegetation classification developed in 2003 by NatureServe and the NPS served as a foundation to further classify and map the vegetation types of the Park. Data from an additional 10 vegetation plots supported vegetation types either rare or not documented in the 2003 classification. Data from 203 verification sites were collected to test the field key to vegetation types and the application of vegetation types to a sample set of map polygons. Furthermore, data from 972 accuracy assessment (AA) sites were collected (of which 966 were used to test accuracy of the vegetation map layer). This GRSM vegetation mapping project identified 112 vegetation types consisting of 105 association types in the U.S. National Vegetation Classification (USNVC), 2 “park-special” types, 1 “map-special” type, and 4 cultural types in the USNVC. To map the vegetation and land cover of GRSM, 52 map classes were developed. Of these 52 map classes, 46 represent natural (including ruderal) vegetation types, most of which types are recognized in the USNVC. For the remaining 6 of the 52 map classes, 4 represent USNVC cultural types for agricultural and developed areas, and 2 represent non-USNVC types for nonvegetated open water and nonvegetated rock. Features were interpreted from viewing four-band digital aerial imagery using digital onscreen three-dimensional stereoscopic workflow systems in geographic information systems; digital aerial imagery was collected during September 23–October 30, 2015. The interpreted data were digitally and spatially referenced, thus making the spatial-database layers usable in a geographic information system. Polygon units were mapped to either a 0.5- or 0.25- hectare (ha) minimum mapping unit, depending on vegetation type. A geodatabase containing several feature-class layers and tables provides the locations and data of USNVC vegetation types (vegetation map layer), vegetation plots, verification sites, AA sites, project boundary extent, and aerial image centers and flight lines. Covering 210,875 ha, the feature-class layer and related tables for the vegetation map layer provide 34,084 polygons of detailed attribute data when special modifiers are not considered (average polygon size of 6.2 ha) and 36,589 polygons of detailed attribute data when special modifiers are considered (average polygon size of 5.8 ha). Each map polygon is assigned a map-class code and name and, when applicable, are linked to USNVC classification tables within the geodatabase. The vegetation map extent includes the administrative boundary for GRSM and the Foothills Parkway. A summary report, generated from the vegetation map layer, concludes that the 46 map classes representing natural (including ruderal) vegetation types apply to 99.2% of polygons (33,797 polygons; average size of 6.2 ha) and cover 98.6% of the Park (207,971.4 ha). Further broken down, map classes representing natural vegetation types indicate that the Park is 97.7% forest and woodland (205,882.5 ha), 0.6% shrubland (1,174.6 ha), and 0.4% herbaceous (914.3 ha). Map classes representing cultural vegetation types apply to 0.8% of polygons (259 polygons; average size of 4.9 ha) and cover 0.6% of the Park (1,277.4 ha). Map classes representing nonvegetation open and flowing water and unvegetated rock apply to 0.08% of polygons (28 polygons; average size of 58.1 ha) and cover 0.8% of the Park (1,625.9 ha). A thematic AA study was completed of map classes representing the natural (including ruderal) vegetation types of the Park. Initial AA results were discussed with NPS staff from the Park. Following input from NPS staff on how to handle map classes that fell below accuracy standards, adjustments were made to the vegetation map layer. Final results indicate an overall accuracy of 80.64% (kappa index of 79.96% for chance agreements) based on data from 966 of the 972 AA sites. Most individual map-class themes exceed the NPS VMI Program standard of 80% with a 90% confidence interval. The GRSM vegetation mapping project delivers many geospatial and vegetation data products, including an in-depth project report discussing methods and results, which includes map classification and map-class descriptions. This suite of products also includes descriptions and a field key to vegetation types; a database of vegetation plots, verification sites, and AA sites; digital images of field sites; field data sheets; digital aerial imagery; hardcopy and digital maps; a geodatabase of vegetation and land cover (map layer), field sites (vegetation plots, verification sites, and AA sites), aerial imagery index, project boundary, and metadata; and a contingency table listing AA results. Geospatial products are projected in the Universal Transverse Mercator, Zone 17 North, by using the North American Datum of 1983. Information on the NPS VMI Program and completed mapping projects are on the internet at https://www.nps.gov/im/vegetation-inventory.htm.

North Carolina, Tennessee↗

USGS2017

Baseline assessment of groundwater quality in Pike County, Pennsylvania, 2015

The Devonian-age Marcellus Shale and the Ordovician-age Utica Shale, which have the potential for natural gas development, underlie Pike County and neighboring counties in northeastern Pennsylvania. In 2015, the U.S. Geological Survey, in cooperation with the Pike County Conservation District, conducted a study that expanded on a previous more limited 2012 study to assess baseline shallow groundwater quality in bedrock aquifers in Pike County prior to possible extensive shale-gas development. Seventy-nine water wells ranging in depths from 80 to 610 feet were sampled during June through September 2015 to provide data on the presence of methane and other aspects of existing groundwater quality in the various bedrock geologic units throughout the county, including concentrations of inorganic constituents commonly present at low values in shallow, fresh groundwater but elevated in brines associated with fluids extracted from geologic formations during shale-gas development. All groundwater samples collected in 2015 were analyzed for bacteria, dissolved and total major ions, nutrients, selected dissolved and total inorganic trace constituents (including metals and other elements), radon-222, gross alpha- and gross beta-particle activity, dissolved gases (methane, ethane, and propane), and, if sufficient methane was present, the isotopic composition of methane. Additionally, samples from 20 wells distributed throughout the county were analyzed for selected man-made volatile organic compounds, and samples from 13 wells where waters had detectable gross alpha activity were analyzed for radium-226 on the basis of relatively elevated gross alpha-particle activity. Results of the 2015 study show that groundwater quality generally met most drinking-water standards for constituents and properties included in analyses, but groundwater samples from some wells had one or more constituents or properties, including arsenic, iron, manganese, pH, bacteria, sodium, chloride, sulfate, total dissolved solids, and radon-222, that did not meet (commonly termed failed or exceeded) primary or secondary maximum contaminant levels (MCLs) or Health Advisories (HA) for drinking water. Except for iron, dissolved and total concentrations of major ions and most trace constituents generally were similar. Only 1 of 79 well-water samples had any constituent that exceeded a MCL, with an arsenic concentration of about 30 micrograms per liter (µg/L) that was higher than the MCL of 10 µg/L. However, total arsenic concentrations were higher than the HA of 2 µg/L in samples from another 12 of 79 wells (about 15 percent). Secondary maximum contaminant levels (SMCLs) were exceeded most frequently by pH and concentrations of iron and manganese. The pH was outside of the SMCL range of 6.5–8.5 in samples from 24 of 79 wells (30 percent), ranging from 5.5 to 9.2; more samples had pH values less than 6.5 than had pH values greater than 8.5. Total iron concentrations typically were much greater than dissolved iron concentrations, indicating substantial presence of iron in particulate phase, and exceeded the SMCL of 300 µg/L more often [35 of 79 samples (44 percent)] than dissolved iron concentrations [samples from 8 of 79 wells (10 percent)]. Total manganese concentrations exceeded the SMCL of 50 µg/L in samples from 31 of 79 wells (39 percent) and the HA of 300 µg/L in samples from 13 of 79 wells (about 16 percent). A few (1–2) samples had concentrations of sodium, chloride, sulfate, or TDS higher than the SMCLs of 60, 250, 250, and 500 mg/L, respectively. However, dissolved sodium concentrations were higher than the HA of 20 mg/L in samples from 15 of 79 wells (nearly 20 percent). Total coliform bacteria were detected in samples from 25 of 79 wells (32 percent) but Escherichia coli were not detected in any sample. Radon-222 activities ranged from 11 to 5,100 picocuries per liter (pCi/L), with a median of 1,440 pCi/L, and exceeded the proposed and the alternate proposed drinking-water standards of 300 and 4,000 pCi/L, respectively, in samples from 60 of 79 wells (75 percent) and in samples from 2 of 79 wells (3 percent), respectively. Groundwater samples from all wells were analyzed for dissolved methane by one contract laboratory that determined water from 19 of the 79 wells (24 percent) had concentrations of methane greater than the reporting level of 0.010 milligrams per liter (mg/L) with a maximum methane concentration of 2.5 mg/L. Methane concentrations in 18 replicate samples submitted to a second laboratory for dissolved gas and isotopic analysis generally were higher by as much as a factor of 2.7 from those determined by the first laboratory, indicating potential bias related to combined sampling and analytical methods, and therefore, caution needs to be used when comparing methane results determined by different methods. The isotopic composition of methane in 9 of 10 samples with sufficient dissolved methane (about 0.3 mg/L) for isotopic analysis is consistent with values reported for methane of microbial origin produced through carbon dioxide reduction; an isotopic shift in 1 or 2 samples may indicate subsequent methane oxidation. The low concentrations of ethane relative to methane in these samples further indicate that the methane may be of microbial origin. Groundwater samples with relatively elevated methane concentrations (near or greater than 0.3 mg/L) also had chemical compositions that differed in some respects from groundwater with relatively low methane concentrations (less than 0.3 mg/L) by having higher pH (greater than 8) and higher concentrations of sodium, lithium, boron, fluoride, arsenic, and bromide and chloride/bromide ratios indicative of mixing with a small amount of brine of probable natural occurrence. The spatial distribution of groundwater compositions differs by topographic setting and lithology and generally shows that (1) relatively dilute, slightly acidic, oxygenated, calcium-carbonate type waters tend to occur in the uplands underlain by the undivided Poplar Gap and Packerton members of the Catskill Formation in southwestern Pike County; (2) waters of near neutral pH with the highest amounts of hardness (calcium and magnesium) generally occur in areas of intermediate altitudes underlain by other members of the Catskill Formation; and (3) waters with pH values greater than 8, low oxygen concentrations, and the highest arsenic, sodium, lithium, bromide, and methane concentrations can be present in deep wells in uplands but most frequently occur in stream valleys, especially at low altitudes (less than about 1,200 feet above North American Vertical Datum of 1988) where groundwater may be discharging regionally, such as to the Delaware River in northern and eastern Pike County. Thus, the baseline assessment of groundwater quality in Pike County prior to gas-well development shows that shallow (less than about 1,000 feet deep) groundwater generally meets primary drinking-water standards for inorganic constituents but varies spatially, with methane and some constituents present in high concentrations in brine (and connate waters from gas and oil reservoirs) present at low to moderate concentrations in some parts of Pike County.

Pennsylvania↗

USGS2016

Baseline assessment of groundwater quality in Wayne County, Pennsylvania, 2014

The Devonian-age Marcellus Shale and the Ordovician-age Utica Shale, geologic formations which have potential for natural gas development, underlie Wayne County and neighboring counties in northeastern Pennsylvania. In 2014, the U.S. Geological Survey, in cooperation with the Wayne Conservation District, conducted a study to assess baseline shallow groundwater quality in bedrock aquifers in Wayne County prior to potential extensive shale-gas development. The 2014 study expanded on previous, more limited studies that included sampling of groundwater from 2 wells in 2011 and 32 wells in 2013 in Wayne County. Eighty-nine water wells were sampled in summer 2014 to provide data on the presence of methane and other aspects of existing groundwater quality throughout the county, including concentrations of inorganic constituents commonly present at low levels in shallow, fresh groundwater but elevated in brines associated with fluids extracted from geologic formations during shale-gas development. Depths of sampled wells ranged from 85 to 1,300 feet (ft) with a median of 291 ft. All of the groundwater samples collected in 2014 were analyzed for bacteria, major ions, nutrients, selected inorganic trace constituents (including metals and other elements), radon-222, gross alpha- and gross beta-particle activity, selected man-made organic compounds (including volatile organic compounds and glycols), dissolved gases (methane, ethane, and propane), and, if sufficient methane was present, the isotopic composition of methane. Results of the 2014 study show that groundwater quality generally met most drinking-water standards, but some well-water samples had one or more constituents or properties, including arsenic, iron, pH, bacteria, and radon-222, that exceeded primary or secondary maximum contaminant levels (MCLs). Arsenic concentrations were higher than the MCL of 10 micrograms per liter (µg/L) in 4 of 89 samples (4.5 percent) with concentrations as high as 20 µg/L; arsenic concentrations were higher than the Health Advisory level of 2 µg/L in 27 of 89 samples (30 percent). Total iron concentrations exceeded the secondary maximum contaminant level (SMCL) of 300 µg/L in 9 of 89 samples (10 percent). The pH ranged from 5.4 to 9.3 and did not meet the SMCL range of greater than 6.5 to less than 8.5 in 27 samples (30 percent); 22 samples had pH values less than 6.5, and 5 samples had pH values greater than 8.5. Total coliform bacteria were detected in 22 of 89 samples (25 percent); Escherichia coli were detected in only 2 of those 22 samples. Radon-222 activities ranged from 25 to 7,400 picocuries per liter (pCi/L), with a median of 2,120 pCi/L, and exceeded the proposed drinking-water standard of 300 pCi/L in 86 of 89 samples (97 percent); radon-222 activities were higher than the alternative proposed standard of 4,000 pCi/L in 12 of 89 samples (13.5 percent). Water from 8 of the 89 wells (9 percent) had concentrations of methane greater than the reporting level of 0.24 milligrams per liter (mg/L) with the detectable methane concentrations ranging from 0.74 to 9.6 mg/L. Of 16 replicate samples submitted to another laboratory with a lower reporting level of 0.0002 mg/L, 15 samples had detectable methane concentrations that ranged from 0.0011 to 9.7 mg/L. Of these 15 samples, low levels of ethane (0.00032 to 0.0017 mg/L) were detected in 6 of 7 samples with methane concentrations greater than 0.75 mg/L. The isotopic composition of methane in 6 of 8 samples with sufficient dissolved methane (about 1 mg/L) for isotopic analysis is consistent with a predominantly thermogenic methane source (sample carbon isotopic ratio δ13C CH4 values ranging from -56.36 to -45.97 parts per thousand (‰) and hydrogen isotopic ratio δD CH4 values ranging from -233.1 to -141.1 ‰). However, the low levels of ethane relative to methane indicate that the methane may be of microbial origin and subsequently underwent oxidation. Isotopic compositions indicated a possibly mixed thermogenic and microbial source (carbon dioxide reduction process) for the methane in 1 of the 8 samples (δ 13 C CH4 of -63.72 and δD CH4 of -192.3 ‰) and potential oxidation of microbial and (or) thermogenic methane in the remaining sample (δ 13 C CH4 of -46.56 and δD CH4 of -79.7 ‰). Groundwater samples with relatively elevated methane concentrations (near or greater than 1 mg/L) had a chemical composition that differed in some respects (pH, selected major ions, and inorganic trace constituents) from groundwater with relatively low methane concentrations (less than 0.75 mg/L). The seven well-water samples with the highest methane concentrations (from about 1 to 9.6 mg/L) also had among the highest pH values (8.1 to 9.3, respectively) and the highest concentrations of sodium, lithium, boron, fluoride, arsenic, and bromide. Relatively elevated concentrations of some other constituents, such as barium, strontium, and chloride, commonly were present in, but not limited to, those well-water samples with elevated methane. Groundwater samples with the highest methane concentrations had chloride/bromide ratios that indicate mixing with a small amount of brine (0.02 percent or less, by volume) similar in composition to that reported for gas and oil well brines in Pennsylvania. Most other samples with low methane concentrations (less than about 1 mg/L) had chloride/bromide ratios that indicate predominantly man-made sources of chloride, such as road salt, septic systems, and (or) animal waste. Although naturally occurring brines may originate from deeper parts of the aquifer system, the man-made sources are likely to affect shallow groundwater. Geochemical modeling showed that the water chemistry of samples with elevated pH, sodium, lithium, bromide, and alkalinity could result from dissolution of calcite (calcium carbonate) combined with cation exchange and mixing with a small amount of brine. Through cation exchange reactions (which are equivalent to processes in a water softener) calcium ions released by calcite dissolution are exchanged for sodium ions on clay minerals. The spatial distribution of groundwater compositions generally shows that (1) relatively dilute, slightly acidic, oxygenated, calcium-carbonate type waters tend to occur in the uplands along the western border of Wayne County; (2) waters of near neutral pH with the highest amounts of hardness (calcium and magnesium) generally occur in areas of intermediate altitudes; and (3) waters with pH values greater than 8, low oxygen concentrations, and the highest arsenic, sodium, lithium, bromide, and methane concentrations can occur in deep wells in uplands but most frequently occur in stream valleys, especially at low elevations (less than about 1,200 ft above North American Vertical Datum of 1988) where groundwater may be discharging regionally, such as to the Delaware River. Thus, the baseline assessment of groundwater quality in Wayne County prior to gas-well development shows that shallow (less than about 1,000 ft deep) groundwater is generally of good quality, but methane and some constituents present in high concentrations in brine (and produced waters from gas and oil wells) may be present at low to moderate concentrations in some parts of Wayne County.

Pennsylvania↗

USGS2021

Status and understanding of groundwater quality in the northern Sierra Nevada foothills domestic-supply aquifer study units, 2015–17—California GAMA Priority Basin Project

Groundwater quality in the northern Sierra Nevada foothills region of California was investigated as part of California State Water Resources Control Board (SWRCB) Groundwater Ambient Monitoring Assessment Priority Basin Project (GAMA-PBP). The region was divided into two study units: the Yuba-Bear watersheds (YBW) study unit and the American-Cosumnes-Mokelumne watersheds (ACMW) study unit. The GAMA-PBP made a spatially unbiased assessment of aquifer systems used for domestic drinking-water supply in the study region, which are predominantly composed of fractured, hard-rock aquifers of varying lithology. These assessments characterized the quality of raw groundwater to evaluate ambient conditions in the domestic-supply aquifer and not the quality of treated drinking water. The study included three components: (1) a status assessment , which characterized the quality of groundwater resources used for domestic drinking-water supply in the YBW and ACMW study units; (2) an understanding assessment , which evaluated natural and anthropogenic explanatory factors that could potentially affect groundwater quality in the study region; and (3) a comparative assessment between the groundwater resources used for domestic and public drinking-water supply in the study region. The status assessment was based on data collected by the GAMA-PBP from 74 sites in the YBW study unit during 2015–16 and 67 sites in the ACMW study unit from 2016 to 2017. To contextualize water-quality results, concentrations of water-quality constituents in ambient groundwater were compared to regulatory and non-regulatory benchmarks typically used by the State of California and Federal agencies as health-based or aesthetic standards for public drinking water. The status assessment used a grid-based method to estimate proportions of groundwater resources with concentrations approaching or exceeding benchmark thresholds. This method provides spatially unbiased results and allows inter-comparability with similar groundwater-quality assessments. Inorganic constituents with health-based benchmarks were present at high relative concentration (RC), meaning they exceeded the benchmark threshold, in 5.4 and 10 percent of domestic-supply aquifer systems in the YBW and ACMW study units, respectively. Inorganic constituents with aesthetic-based benchmarks were detected at high-RCs in 20 and 28 percent of the YBW and ACMW study units, respectively. The inorganic constituents present at high RC were arsenic, barium, boron, molybdenum, strontium, nitrate, adjusted gross-alpha particle activity, chloride, total dissolved solids, specific conductance, iron, manganese, and hardness. Groundwater samples were tested for presence or absence of three microbial indicators (total coliform, Escherichia coli , and Enterococci ). At least one microbial indicator was present in 26 and 28 percent of the YBW and ACMW study units, respectively. At least one organic constituent was detected in 30 and 42 percent of the YBW and ACMW study units, respectively. Organic constituents were not present at high RC, but tetrachloroethene (PCE), trichloroethene (TCE), and toluene were detected in the YBW study unit at moderate RC (between the benchmark concentration and one-tenth of the benchmark concentration). Methyl tert -butyl ether (MTBE) and chloroform were present at low RC (less than one-tenth of the benchmark concentration) in the YBW and ACMW study units with detection frequencies greater than 10 percent. Perchlorate, a constituent of special interest, was detected in 31 and 41 percent of the YBW and ACMW study units, respectively, at either low or moderate RCs. Relations among select water-quality constituents and potential explanatory factors were evaluated using statistical and graphical approaches. Nitrate, microbial indicators, and perchlorate were all correlated to elevation-dependent variables relating to climate, land use, and recharge condition. Isotopic and dissolved noble-gas tracers indicated these water-quality constituents are associated with recharge conditions associated with irrigation during the summer dry-season, which is common in areas of rural-residential or agricultural land uses. Higher concentrations of iron and manganese were primarily associated with anoxic groundwater in aquifers of metasedimentary lithology. Increased hardness was primarily associated with anoxic groundwater in aquifers of mafic-ultramafic or metavolcanics lithologies at lower elevations in the study region in the Melones fault zone. Chloroform and MTBE were associated with shallow groundwater (wells depths less than 130 m) under oxic and anoxic redox conditions, respectively. The comparative assessment evaluated differences between the aquifer systems used for domestic- and public-supply in study region based on (1) well-construction characteristics, and (2) water quality. Analysis of over 60,000 well-completion reports in the study region showed that although domestic-supply wells span the deepest depth zones in regional aquifers, median depths for public-supply wells were significantly greater than those of domestic-supply wells in both study units. Water-quality data from more than 300 public-supply wells in the study region were assessed using a spatially weighted method for calculation aquifer-scale proportions and compared with the domestic-supply assessment results. Detections of inorganic constituents at high RC and detection frequencies for organic constituents were generally similar between the domestic- and public-supply aquifer systems in both study units, with a few notable exceptions in the ACMW study unit: nitrate was greater for the public- compared to domestic-supply aquifer system and both manganese, hardness, and MTBE were greater in the domestic- compared to public-supply aquifer system. These differences are likely related to contrasting land uses, aquifer lithologies, landscape positions, and depths characterizing domestic- and public-supply wells in the ACMW study unit. Overall, fewer samples from domestic-supply wells in the northern Sierra Nevada foothills exceeded health-based benchmarks compared to aesthetic-based benchmarks for groundwater quality. Exceedences of health-based benchmarks were primarily caused by nitrate and coliform bacteria, which were associated with recharge from diverted surface water used primarily for irrigation. Exceedences of aesthetic-based benchmarks were primarily caused by iron, managanese, and hardness, which were associated with geologic factors. Regional irrigation practices and aquifer lithology can affect groundwater quality in fractured-rock aquifers in the northern Sierra Nevada foothills used for domestic drinking-water supply.

California↗

USGS2011

Thermal profiles for selected river reaches in the Stillaguamish River basin, Washington, August 2011

Datums Horizontal coordinate information is referenced to the North American Datum of 1983 (NAD 83). Abstract Longitudinal profiles of near-streambed temperature were collected for eight river reaches in the Stillaguamish River basin, Washington, during August 2011, to provide information about areas of groundwater discharge to streams. During summer, groundwater discharge can be a source of cold water to streams that regulates warm stream temperatures creating cold-water thermal refugia for native stream biota including salmon and trout. To assess areas of groundwater discharge to streams, temperature was measured using a probe with an internal datalogger towed behind a watercraft moving downstream at ambient stream velocity. The data were referenced to location, concurrently surveyed with a Global Positioning System, during collection of the water temperature data. Data are presented as Microsoft Excel® files consisting of date and time, near-streambed water temperature, and latitude and longitude. Introduction Longitudinal profiles of near-streambed temperatures surveyed at ambient river velocity in a Lagrangian framework provide information about potential areas of groundwater discharge as well as salmonid habitat and thermal refugia (Vaccaro and Maloy, 2006). Longitudinal thermal profiles have previously been surveyed in several rivers in Washington, including the Yakima River and tributaries (Vaccaro and others, 2008) and the Nooksack River (Cox and others, 2005). This report presents eight thermal profiles within the Stillaguamish River basin including parts of the North Fork Stillaguamish River, South Fork Stillaguamish River, Jim Creek, and Pilchuck Creek (fig. 1). This data augments previous investigations of longitudinal temperature variations within the Stillaguamish River and tributaries by thermal infrared radar by the Washington State Department of Ecology (Watershed Sciences, 2002), and may be used as a tool to develop a better understanding of groundwater/surface-water interactions within the Stillaguamish River basin. Purpose and Scope The purpose of this report is to present longitudinal thermal profiles of stream temperature of streams within the Stillaguamish River basin including the North Fork Stillaguamish River, the South Fork Stillaguamish River, Pilchuck Creek, and Jim Creek. This data may be used to determine zones of groundwater discharge and improve understanding of the relation between the groundwater and surface water systems of the Stillaguamish River basin. Description of Study Area The Stillaguamish River basin is in northwestern Washington and is bounded to the east by the Cascade Mountains, to the west by Puget Sound, to the north by the Skagit River basin, and to the south by the Snohomish River basin (fig. 1). The Stillaguamish River basin is characterized by cool, wet winters and warm, dry summers. Mean annual discharge of the North Fork Stillaguamish River (North Fork Stillaguamish River near Arlington, Washington, USGS gaging station 12167000) for water years 1929-2010 is 1,898 ft 3 /s and mean annual discharge of the South Fork Stillaguamish River (South Fork Stillaguamish River near Granite Falls, Washington gaging station 12161000) for water years 1929-1980 is 1,071 ft 3 /s. Jim Creek is a tributary of the South Fork Stillaguamish River and Pilchuck Creek is a tributary of the mainstem Stillaguamish River. Thermal Profile Survey Continuous water temperature and Global Positioning System (GPS) data were collected at 3-second intervals while drifting downstream at ambient stream velocity in a Lagrangian framework following the method of Vaccaro and Maloy (2006) for Pilchuck Creek between river mile (RM) 0.0 and 3.7 (table 1); the North Fork Stillaguamish River between RM 0.0 and 34.2 (tables 2-5); South Fork Stillaguamish River between RM 17.7 and 33.4 (tables 6-7); and Jim Creek between RM 0.0 and 7.0 (table 8). Profiling at ambient stream velocity in a Lagrangian framework tracks a parcel of water as it moves downstream during the day; departures from the diurnal heating cycle may be due to groundwater input, surface-water inflows, or riparian shading. Continuous temperature was measured using a Solinst® Levelogger LT temperature probe verified by a National Institute of Standards and Technology (NIST) certified thermistor and position data was measured using a Garmin® GPSmap® 60Csx for the eight surveys during August 15-26, 2011. The temperature probe was towed behind a watercraft following the stream thalweg and dragged along the streambed except when in-stream obstacles prevented probe movement downstream. The location of each temperature measurement was determined by relating the time stamp of the GPS data to the temperature data. If a GPS location was not recorded at the same time as a temperature measurement, the location of the temperature measurement was determined by linear interpolation of the two GPS known locations that bracket the time of the temperature measurement. A 0.5-mi gap exists between the beginning of the North Fork Stillaguamish datasets collected on August 18 (table 4) and August 22 (table 5) because of inadequate equilibration of the temperature probe to ambient stream temperature during the initial 0.5 mi of the August 22 survey. Distribution of Information An Excel file of tables 1-8 that include the thermal-profile data for each longitudinal thermal profile is available at http://pubs.usgs.gov/ds/654/ds654_tables.xls. Table 1. Temperature and Global Positioning System location data for the Pilchuck Creek (RM 0.0-3.7), August 15, 2011. Table 2. Temperature and Global Positioning System location data for the North Fork Stillaguamish River (RM 30.0-34.2), August 16, 2011. Table 3. Temperature and Global Positioning System location data for the North Fork Stillaguamish River (RM 17.6-30.0), August 17, 2011. Table 4. Temperature and Global Positioning System location data for the North Fork Stillaguamish River (RM 9.5-17.6), August 18, 2011. Table 5. Temperature and Global Positioning System location data for the North Fork Stillaguamish River (RM 0.0-9.0), August 22, 2011. Table 6. Temperature and Global Positioning System location data for the South Fork Stillaguamish River (RM 25.9-33.4), August 24, 2011. Table 7. Temperature and Global Positioning System location data for the South Fork Stillaguamish River (RM 17.7-25.9), August 26, 2011. Table 8. Temperature and Global Positioning System location data for Jim Creek (RM 0.0-7.0), August 25, 2011. References Cited Cox, S.E., Simonds, F.W., Doremus, L., Huffman, R.L., and Defawe, R.M., 2005, Ground water/surface water interactions and quality of discharging ground water in streams of the lower Nooksack River Basin, Whatcom County, Washington: U.S. Geological Survey Scientific Investigations Report 2005-5255, 46 p Vaccaro, J.J., Keys, M.E., Julich, R.J., and Welch, W.B., 2008, Thermal profiles for selected river reaches in the Yakima River basin, Washington: U.S. Geological Survey Data Series 342 (Available at http://pubs.usgs.gov/ds/342/). Vaccaro, J.J., and Maloy, K.J., 2006, A thermal profile method to identify potential ground-water discharge areas and preferred salmonid habitats for long river reaches: U.S. Geological Survey Scientific Investigations Report 2006-5136, 16 p. Watershed Sciences, LLC, 2002, Aerial surveys in the Stillaguamish and Skagit River Basins-Thermal infrared and color videography: Corvallis, Oreg., Water Sciences, for Washington Department of Ecology, 28 p.

Data Series↗

USGS2001

The geology of a part of Acadia and the nature of the Acadian orogeny across Central and Eastern Maine

The zone of Acadian collision between the Medial New England and Composite Avalon terranes is well preserved in Maine. A transect from northwest (Rome) to southeast (Camden) crosses the eastern part of Medial New England comprising the Central Maine basin, Liberty-Orrington thrust sheet, and Fredericton trough, and the western part of Composite Avalon, including the Graham Lake, Clarry Hill, and Clam Cove thrust sheets. U-Pb geochronology of events before, during, and after the Acadian orogeny helps elucidate the nature and distribution of tectonostratigraphic belts in this zone and the timing of some Acadian events in the Northern Appalachians. The Central Maine basin consists of sedimentary and volcanic rocks of Middle Ordovician (∼470 to ∼460 Ma) age overlain with probable conformity by latest Ordovician(?) through earliest Devonian marine rift and flysch sedimentary rocks; these are intruded by weakly to undeformed plutonic rocks of Early and Middle Devonian age (∼399–378 Ma). The Fredericton trough consists of Early Silurian gray pelite and sandstone to earliest Late Silurian calcareous turbidite, deformed and variably metamorphosed prior to the emplacement of Late Silurian (∼422 Ma) and Early to Late Devonian (∼418 to ∼368 Ma) plutons. The Liberty-Orrington thrust sheet consists of Cambrian(?)-Ordovician (>∼474 to ∼469 Ma and younger) clastic sedimentary and volcanic rocks intruded by highly deformed Late Silurian (∼424 to ∼422 Ma) and Devonian (∼418 to ∼389 Ma) plutons, possibly metamorphosed in Late Silurian time (prior to ∼417 Ma), and metamorphosed to amphibolite facies in Early to Middle Devonian time (∼400 to ∼381 Ma). The Graham Lake thrust sheet contains possible Precambrian rocks, Cambrian sedimentary rocks with a volcanic unit dated at ∼503 Ma, and Ordovician rocks with possible Caradocian Old World fossils, metamorphosed and deformed in Silurian time and intruded by mildly to undeformed Late Silurian (∼421 Ma) and Late Devonian (∼371 to ∼368 Ma) plutons. The Clarry Hill thrust sheet consists of poorly studied, highly metamorphosed Cambrian (?) rocks. The Clam Cove thrust sheet contains highly deformed Precambrian limestone, shale, sandstone, and conglomerate, metamorphosed to epidote amphibolite facies and intruded by a mildly deformed pluton dated at ∼421 Ma. Metamorphism, deformation, and voluminous intrusive igneous activity of Silurian age are common to both the most southeastern parts of Medial New England and the thrust sheets of Composite Avalon. In contrast to Medial New England, the thrust sheets of Composite Avalon show only modest effects of Devonian deformation and metamorphism. Regional stratigraphic relations, paleontologic findings, and U-Pb geochronology suggest that the Graham Lake, Clarry Hill, and Clam Cove thrust sheets are far-traveled allochthons that were widely separated from Medial New England in the Silurian. One hundred nine (109) new U-Pb analyses of zircon, monazite, and sphene from 25 samples of metamorphosed, stratified, and intrusive igneous rocks are used to decipher the history of events along the transect. In our view, many of the structures within Medial New England were formed during the closing of the back-arc ocean during the waning stages of the Taconian orogeny. These features include the precursors to the upright folds found in Silurian sedimentary rocks, the west-facing Liberty-Orrington thrust sheet, and the Silurian plutons of Medial New England. Acadian tectonic features, including flysch sedimentation, igneous activity, deformation with nappe emplacement, and metamorphism record the progressive loading of Medial New England by a stack of thrust nappes emplaced in latest Silurian to Middle Devonian time. The Acadian orogeny is a prolonged event, lasting from earliest Late Silurian to the Late Devonian, whose evolution involved: (1) convergence between Medial New England and Composite Avalon along an east-dipping subduction zone from earliest Late Silurian to Early Devonian time; (2) collision and concurrent delamination of lithospheric mantle beneath Medial New England in Early Devonian time resulting in deformation, high-grade metamorphism, and intrusive igneous activity in the most eastern part of Medial New England and the western parts of Composite Avalon; (3) Early to Middle Devonian northwest-migrating penetrative deformations of the Acadian Main stage, including northwest-directed thrusting and recumbent folding followed by tightening of folds possibly produced in the waning stages of the Taconian orogeny and forming folded wave trains of isoclinal folds (∼419 to ∼404 Ma), (4) asymmetrical folds produced by east-west shortening (∼399 to ∼380 Ma); and (5) final westward emplacement of Composite Avalon thrust sheets onto Medial New England in Late Devonian time (∼380 to ∼371 Ma). The present boundary between rocks of Medial New England and Composite Avalon at the surface is the Sennebec Pond fault, a high-angle fault that cuts the Graham Lake and Clarry Hill thrusts, and is intruded by the Mt. Waldo pluton (∼371), one of several Late Devonian plutons (∼371 to ∼367 Ma) that mark the end of the Acadian orogeny in coastal Maine. The thrust sheets of Composite Avalon carried rocks of peri-Gondwanan affinity and with possible Old World fauna many kilometers to the west over Medial New England, thus potentially hiding the original suture at depth under the Gulf of Maine.

Maine↗

USGS1995

Geologic framework of a transect of the central Brooks Range: Regional relations and an alternative to the Endicott Mountains allochthon

This paper evaluates the geologic framework and tectonic development of the central Brooks Range based on a transect through the range and Arctic foothills. A geologic cross section constructed through the transect is confirmed by comparing the retrodeformed section with the regional distribution of lithofacies in the central Brooks Range. Stratigraphic relations in the retrodeformed section are further explained by comparing them to similar stratigraphic relations in the Ikpikpuk-Umiat basin under the Arctic coastal plain. The structural framework of the central Brooks Range and Arctic foothills consists of fold nappes, thrust faults, and detached folds that sole in decollements and late-stage high-angle faults. In the central Brooks Range, shortening is by north-directed thrust faulting and folding of mostly Paleozoic rocks, and transport of any individual thrust sheet relative to underlying rocks is less than 30 km. In the middle of the range, imbricate blocks of lower Paleozoic basement are exposed in the core of the Doonerak anticline, and thrust sheets of stratigraphically higher Paleozoic rocks that overlie basement are exposed in the limbs of the anticline. In the northeast part of the anticline, the Amawk thrust emplaces Silurian and Upper Devonian rocks on a succession of Lower Mississippian an stratigraphically higher rocks that have been detached from the underlying basement along the Blarney Creek thrust. The Slatepile fault system, a system of high-angle faults in the north limb of the Doonerak anticline, drops the core and part of the north limb of the anticline down, giving the impression that the succession of Lower Mississippian and stratigraphically higher rocks that lie on basement south of the system high-angle faults extends under the Upper Devonian rocks that extensively crop out north of the high-angle faults. In the Arctic foothills, the mostly Paleozoic rocks of the north-central Brooks Range extend under Lower Cretaceous rocks of the North Slope foreland basin, and blind thrusts that sole in the Paleozoic rocks ramp up into the Lower Cretaceous and stratigraph cally higher rocks. Also in the Arctic foothills, a thrust sheet that contains the Arctic foothills assemblage overlies rocks of the north-central Brooks Range and Lower Cretaceous rocks of the North Slope foreland basin. Thrust transport of the Arctic foothills assemblage more than 40 km from south of the Doonerak anticline took place during the Early Cretaceous, but thrusting that deformed rocks of the North Slope foreland basin took place during the early Tertiary, with the vertical uplift of the Doonerak anticline being a late-formed feature. Conclusions based on the retrodeformed cross section contrast significantly with previous work in which the Upper Devonian and stratigraphically higher rocks north of the Doonerak anticline are considered part of the Endicott Mountains allochthon, a regional allochthon that extends the breadth of the Brooks Range. In these models, Upper Devonian and younger rocks in the north-central Brooks Range have been thrust-transported 90 or 200 km from south of the Doonerak anticline, and emplacement of the allochthon could reflect as much as 885 km of tectonic shortening. The Lower Mississippian and stratigraphically higher rocks together with the underlying basement in the northeast part of the anticline are considered to be in a window in the Endicott Mountains allochthon and to extend northward beneath allochthonous Upper Devonian and stratigraphically higher rocks in the north-central Brooks Range. Lithofacies patterns in rocks in the central Brooks Range are consistent with the retrodeformed cross section and imply plausible Upper Devonian and Carboniferous depositional systems. Thick Upper Devonian and Lower Mississippian(?) clastic prisms were deposited in basins north of the Doonerak anticline. Mississippian carbonate rocks that overlie these clastic prisms were deposited in differentially subsiding shelf environments that included rocks in the Doonerak anticline. Restored across the Blarney Creek thrust, the Mississippian shelf carbonate rocks that presently lie north of the Doonerak anticline are those that were deposited on basement in the anticline. A carbonate ramp at the south edge of these shelf deposits extends east-southeast across the central Brooks Range and in th retrodeformed section lies south of the Doonerak anticline where Upper and Middle(?) Devonian shaly rocks thicken to the south. Unrestored, the ramp would extend across the Doonerak anticline. Restored Late Devonian and Carboniferous lithofacies patterns in the central Brooks Range also are plausible from a regional perspective and have implications for exploration of basins under the Arctic coastal plain. The Late Devonian to Early Mississippian(?) basins in the north-central Brooks Range are part of a system of Early(?) Devonian to Early Mississippian(?) clastic basins that extend the length of the Brooks Range and include basins under the Arctic coastal plain. These basins are a template for depositional patterns in overlying rocks. Marine shelves between these basins where Mississippian strata unconformably lie on basement, such as in the northeastern Brooks Range and the Doonerak anticline, have depositional histories that are in contrast to areas that overlie the basi s. The resulting stratigraphic framework, together with the structural framework in the basement rocks that controlled the basins, has had a profound effect on the structural style of the fold belt, the salient effect being folds and thrust faults that are not orthogonal to the direction of structural transport. Stratigraphic relations exposed in the fold belt, especially the distribution of potential source rocks, likely model little-explored basins that underlie the North Slope foreland basin.

Alaska↗

USGS2020

Sediment lithology and borehole erosion testing, American and Sacramento Rivers, California

Executive Summary A field investigation intended to measure the potential for erosion of sediments beside the American and Sacramento Rivers near Sacramento, California, is described. The study featured two primary components: (1) drilling and soil sampling to reveal lithology, down to depths matching the local river thalweg, where possible, and (2) borehole erosion tests (BETs) as described by Briaud and others (2017) at many of the same locations. The latter test involved drilling a vertical hole, measuring its diameter profile, inserting a hollow drilling rod to almost the bottom of the hole, and pumping fluid through the drilling rod at a known discharge for a chosen time interval. The hole was then resurveyed to establish an erosion rate (change in borehole radius divided by duration of flow event) as a function of depth, and the test was repeated. This test was performed with water as the erosive fluid at 12 locations, with 1 test repeated with drilling mud. Lithology holes were drilled at these same locations and an additional five locations. Drilling operations took place on river left and river right on the American River and river left (left bank, when looking downstream) on the Sacramento River. The drilling to acquire sediment samples and reveal lithology involved the use of a mobile drilling rig equipped with a 6-inch (in.) auger, a 140-pound pneumatic hammer to drive split spoon and Calmod samplers, and a piston to push Shelby tube samplers to obtain samples of clayey material. Blow count (hammer blows per 6-in. sampler advance) was recorded while sampling, and the process was logged using standard U.S. Army Corps of Engineers (USACE), Sacramento District procedures. Sediment samples were identified and described in the field per ASTM D2488 and then delivered to a USACE laboratory and to Texas A&M University for additional laboratory analysis. The BETs were performed with the same drilling rig that performed the drilling for definition of lithology. In most instances, tests were limited to regions above the water table, to avoid slumping of the borehole and heaving sands pushing into the hole. Most of the tests featured sediments that were primarily silty sand or sandy silt. The testing procedure involved comparing borehole profiles before and after passing an assumed constant discharge through a drilling rod to the bottom of the drilled hole. Discharge and water losses were logged during the testing procedure, and water losses into the walls of the drilled hole were typically less than 5 percent of the introduced volume. For the tests performed with water, the coefficient of variation of the discharge ranged from 4.5 to 28 percent, with a mean of 13 percent, but the mean discharge appeared to be reasonably steady over the typical test duration of 10–30 minutes. It was thus assumed that discharge was constant and water losses during the tests were neglected. Coefficients of variation of the discharge for the three tests performed with drilling mud were much higher (20–50 percent), but erosion rates were much smaller. Resolution of the borehole caliper-reported diameter was 0.1 in. and several of the tests lasted for 10 minutes. With boreholes measured twice, before and after each test, and averaged, these numbers correspond to an apparent erosion rate (radius change divided by test duration) of 0.3 inches per hour (in/hr), which is a theoretical lower bound on what could be measured with this approach and equipment. In practice, 0.5 in/hr appears to be a more realistic lower bound on the detectable erosion rate, based on inspection of computed changes and erosion rates. Three flow speeds (5, 8, and 12 feet per second; ft/s) were targeted for the tests. Because of equipment limitations, it was not possible in the field to reach an average of 12 ft/s throughout any given borehole, although much higher flow speeds were reached locally in some cases. Most tests featured at least two different flow rates, and the borehole was typically surveyed at least twice for each condition, to allow averaging to reduce the influence of random diameter measurement errors. Errors arising from out-of-round boreholes appeared to be uncommon. Briaud and others (2017) recommend stepped increases in the flow rate during a borehole test. This approach was taken during initial testing but proved to be problematic. The drilled hole would be enlarged by the first (smaller) discharge, and then it would be difficult to reach the desired higher flow speed because of the larger annulus between the drilled hole and the drilling rod that supplied the water for testing. This was largely solved by starting with a high discharge and, in many cases, maintaining it for subsequent tests with the average flow speed decreasing as the hole enlarged. Several different measures of erosion rate were computed and investigated by comparison to lithological profiles. The vertically averaged erosion rate for each hole was computed, but this result does not reveal vertical variability of erodibility; and the mean flow speed within the hole is not a good representation of the speed when attempting to determine a relationship between erosion rate and flow speed. Instead, for each 6-inch layer within the hole, vertically averaged erosion rates and local flow speeds were computed and plotted. Where possible, the soil type for each layer was identified. For later laboratory analysis, project protocol dictated collection of Shelby tube samples whenever clay was encountered. Plots of erosion rate versus flow speed displayed scatter that indicate that several other factors influence the erosion potential of the soil. Blow count was not a good predictor variable; it is better correlated with soil type than erodibility. Soils were classified as sand, silt, or clay, depending on which soil type dominated within a sample. In general, those classified as sand and silt did not reveal clear patterns allowing erosion rate to be computed directly from flow speed, but the test results define the range and bounds on the erosion rate. Results for clay were slightly clearer with the erosion rate increasing with flow speed, once a threshold had been reached. In this case, the erosion rate appeared to change near a speed of 7 ft/s; above this threshold, erosion rates jumped from less than 2 in/hr to greater than 3 in/hr. Even for soils with similar classifications, large differences in erodibility were observed between sites and in different layers within an individual hole. One potential means of dealing with this problem would be to perform more tests at each site to allow establishment of relationships between flow speed and erodibility for individual layers within a borehole. The maximum number of tests performed at a site in this study was four, but in some cases, results are available for only one or two flow events. Comparison of data to a set of Erosion Function Apparatus tests that provide better resolution of the vertical variation in the erosion rate versus flow speed relationship would allow further investigation of this idea. It was hypothesized that drilling mud could expand the utility of the test in soft sands by reducing the likelihood of slumping that would be interpreted as erosion. The one test that was performed with drilling mud indicated that it greatly reduced the erosion rate of the soils encountered. It yielded very different results from the test performed at the same site with water. Erosion rate is often expressed as a function of shear stress applied to a soil. In order to compute shear stress on the walls of the drilled hole, one must assume a form for the relationship between flow speed and shear stress and select a friction factor that is often estimated empirically from head loss, observed water-surface profiles, surface roughness, or other data not available in this report. One methodology for computing shear stress from flow speed is discussed in this report, but the test results have been presented in terms of erosion rate versus flow speed to avoid assuming values that are not verifiable via the field data collected in this study. Erosion rate was computed from directly measured values (sequential borehole profiles) and flow speed was computed directly from measured quantities (discharge and borehole geometry). The BET has seen limited application, primarily in clayey soils, whereas most of the soils encountered in this study were primarily sand or silt. The objective of the BET is to determine the erodibility of in situ soil below the ground or riverbed surface. The BET is simple in principle and has the advantage of revealing erodibility of in situ sediments below the ground or riverbed surface; it appears to be very useful in clayey soils, based on previously published work, but is more difficult to apply in sandy soils where slumping and water losses within the hole during testing are more likely to occur. The BET did reveal a large variation in the results both laterally and vertically, even for the same soil-type classification. It is thus recommended that the results be applied considering these spatial variations rather than attempting to universally assign an erosion-rate relationship to a particular soil type. Results have been provided showing the results by site and by sediment classification (sand, silt, and clay), to allow either approach. Where possible, it is important to rely on site-specific results because the erosion-rate relationship for a given soil type varied by site. Data collected during this project have been made publicly available online via the U.S. Geological Survey (USGS) Sciencebase database. The measured borehole profiles, discharge, lithology log sheets, and photos are available in the data release that accompanies this report (see Work and Livsey (2019) in the “Selected References” section for the appropriate link).

California↗