SEARCH · USGS Science
Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
High-resolution airborne radiometric surveys are covering more ground than ever to provide insights into unseen geology, mineral resource potential, and possible health hazards.
Contraction along the Yakutat–North America plate boundary drives 4.6–9.0 mm/year Holocene rock uplift rates along Earth's fastest slipping (≥49 mm/year) ocean–continent transform fault, the Fairweather Fault. Between Icy Point and Lituya Bay, the near-vertical Fairweather fault focuses rock uplift and rapid right-lateral slip by accommodating both vertical and fault-parallel strain during oblique-slip and separate, predominantly strike-slip ruptures. Unusually high uplift rates, indicated by radiocarbon and luminescence dating, result from a 10-km-wide, asymmetric, positive flower structure along a 20°, ∼30-km-long restraining double bend in the Fairweather fault. The principal reverse fault in the flower structure, the offshore, blind Icy Point–Lituya Bay fault, ruptures no more than every 460–1040 years evidenced by uplifted Holocene shorelines. Maximum 3–5 m coseismic uplifts imply 3.1–10 m dip slip per event and earthquake magnitudes of M w 7.0–7.5. The Yakutat block collides obliquely into North America, and our model entails oblique slip on the Fairweather fault with and without corupture on the reverse fault. Oblique slip is evident by vertically offset (>25 m) fluvial and marine terraces and by the primary Fairweather fault strand that strikes >20° to the west of plate-boundary motion.
n the Gulf of Alaska, a series of marine heat waves depleted Pacific cod ( Gadus macrocephalus ) biomass to the lowest abundance ever recorded and led to the fishery’s closure in 2020. Although the fishery has been productive for decades, this collapse may have historical precedents. Traditional knowledge holders refer to cod as ‘the fish that stop’, and there is a suggested period of decline in the 1930s. Here we conduct a catch reconstruction of the early commercial fishery (1864–1950), confirming a rapid catch decline in the 1920s and 1930s. Next, we evaluate evidence for possible drivers. We document changes to demand and technology that contributed to declining catch. However, we also find both qualitative and quantitative evidence of depletion, suggesting catch declines were not driven entirely by social factors. Overfishing may have contributed to localized catch declines as evidenced by declining catch rates in heavily fished localities. We also find evidence for climate as a driver of regional decline, with the period of catch decline characterized by up to 2°C higher temperatures as compared to the earlier period of high fisheries production. Our analysis underscores the importance of understanding long-term drivers of fisheries productivity and the value of linking fisheries and climate histories.
Introduction The bedrock geology of the Bellows Falls 7.5- x 15-minute quadrangle, Vermont and New Hampshire, consists of polydeformed Ordovician to Devonian metasedimentary, metavolcanic, and metaplutonic rocks of the Connecticut Valley trough, Bronson Hill anticlinorium (or Bronson Hill terrane), and the Central Maine terrane. Previous work in this area includes a 1:62,500-scale published map and text (1946), State geologic maps of New Hampshire (1997) and Vermont (2011), and various maps and reports presented largely as parts of field trip guidebooks. A provisional open-file map of the geology of the Vermont part of the Bellows Falls 7.5- x 15-minute quadrangle, completed in 1997, is incorporated and revised on this map based on additional fieldwork. This study recognizes three major structural levels from west to east and from lowest to highest: (1) autochthonous rocks of the Connecticut Valley trough; (2) allochthonous rocks of the New Hampshire sequence and Bronson Hill arc in the Monroe thrust sheet (or nappe), including the Skitchewaug nappe; and (3) allochthonous rocks of the Fall Mountain thrust sheet or nappe.
Coseismic slip on the Patton Bay splay fault system during the 1964 M w 9.2 Great Alaska Earthquake contributed to local tsunami generation and vertically uplifted shorelines as much as 11 m on Montague Island in Prince William Sound (PWS). Sudden uplift of 3.7–4.3 m caused coastal lagoons along the island's northwestern coast to gradually drain. The resulting change in depositional environment from marine lagoon to freshwater muskeg created a sharp, laterally continuous stratigraphic contact between silt and overlying peat. Here, we characterize the geomorphology, sedimentology, and diatom ecology across the 1964 earthquake contact and three similar prehistoric contacts within the stratigraphy of the Hidden Lagoons locality. We find that the contacts signal instances of abrupt coastal uplift that, within error, overlap the timing of independently constrained megathrust earthquakes in PWS—1964 Common Era, 760–870 yr BP, 2500–2700 yr BP, and 4120–4500 yr BP. Changes in fossil diatom assemblages across the inferred prehistoric earthquake contacts reflect ecological shifts consistent with repeated draining of a lagoon system caused by >3 m of coseismic uplift. Our observations provide evidence for four instances of combined megathrust-splay fault ruptures that have occurred in the past ∼4,200 years in PWS. The possibility that 1964-style combined megathrust-splay fault ruptures may have repeated in the past warrants their consideration in future seismic and tsunami hazards assessments.
The depth level at which porphyry Cu–forming magmas fractionated and exsolved mineralizing fluids is actively debated. In the classic model, extensive magma fractionation occurs in large, upper crustal magma chambers, and concomitant fluid exsolution leads to forceful expulsion of residual magmas in the form of porphyry dikes, stocks, and breccia pipes, which subsequently serve as pathways for the mineralizing fluids. In contrast, some recent studies highlighting the role of deep crustal magma fractionation in the production of fertile magmas essentially deny the existence of upper crustal magma chambers at the time of mineralization. To address this, we conducted a detailed thermobarometric investigation of 13 intermediate to felsic, porphyritic intrusive rocks related to porphyry-skarn Cu mineralization at Santa Rita and Hanover-Fierro, New Mexico, United States, representing two premineralization magmas (61–60 Ma), seven synmineralization magmas (60–58 Ma), and four late- to postmineralization magmas (58–57 Ma). For each sample, the pressure of last magma crystallization before final magma ascent to the current exposure level was reconstructed based on Al-in-hornblende barometry of small hornblende inclusions trapped within quartz phenocrysts and through titanium-in-quartz (TitaniQ) thermobarometry of the host quartz phenocrysts themselves. Since quartz is one of the last crystallizing magmatic minerals, and no significant phenocryst growth could have occurred in small dikes and stocks after final magma emplacement, quartz phenocrysts and their contained hornblende inclusions record the depth of last magma crystallization before final magma ascent. When present, hornblende phenocrysts and hornblende inclusions within other major phenocrysts were also analyzed. Both quartz and hornblende barometers return consistent average pressures of 3.2 ± 0.4 kbar for the entire suite of pre- to postmineralization magmas, corresponding to depths of 11 to 14 km. The synmineralization magmas return even more consistent average pressures of 3.1 ± 0.2 kbar, corresponding to a depth of 12 ± 1 km. The volume of the mineralizing porphyry dikes and stocks at the emplacement level is far too small to have provided all the fluids and metals required to form the observed ore deposits. Therefore, the majority of the ore-forming fluids must have originated from the magmas that crystallized at 12 ± 1 km depth. The ore deposits, conversely, formed at ~5-km paleodepth. This implies that most of the mineralizing fluids traveled an average vertical distance of ~7 km from their magmatic source to the eventual site of ore precipitation. The relatively unaltered nature and low veining degree of deeper parts of mineralized porphyry dikes and stocks suggest that the fluid transport through these intrusive bodies occurred mostly at near-solidus conditions by means of fluid percolation along grain boundaries. In summary, our results suggest that (1) a large, upper crustal pluton exists ~7 km beneath the Santa Rita and Hanover-Fierro deposits; (2) abundant phenocryst crystallization occurred at this depth level; and (3) this pluton was the main source for the exsolution of ore-forming fluids. However, the investigated rocks have elevated whole-rock Sr/Y ratios, indicating magma fractionation at deep crustal levels. As a result, our preferred model is a combination of the two end-member models introduced above, with most magma fractionation having occurred in the deep crust and with residual, intermediate to felsic melts having ascended and accumulated at 11 to 14 km paleodepth, where they continued to crystallize with comparatively little crystal-liquid separation, before some of these magmas ascended further to shallow levels and quenched to porphyries.
The surficial geologic map of the Owlshead Mountains 30' x 60' quadrangle depicts the distribution and characteristics of surficial-deposit materials and neotectonic deformation for an area of approximately 5,000 square kilometers (km 2 ) located in the western Basin and Range Province of eastern California. The map represents a new compilation of the surficial geology that encompasses deposits within the late Pliocene to Quaternary. The map is based primarily on new mapping conducted between 2001 and 2009. Map compilation was supported by field observations distributed across the map area, combined with reference to several published and unpublished mapping sources that mostly emphasized neotectonic deformation. The surficial-deposit units included in the map follow a classification scheme that systematically denotes depositional process, relative age, and any secondary sedimentologic or morphologic characteristics. Identification, correlation, and age estimation of map units are based primarily on the relative degree of development of certain time-dependent characteristics such as surface morphology, including local dissection and surface preservation, surface clast modification, and degree of soil development; these characteristics are implicitly incorporated into unit designations. The map represents a detailed and regionally uniform synthesis of the late Neogene geology for this large area that provides a framework applicable to many interpretative studies, such as regional patterns of deposition and dissection; surface drainage development and evolution; and the distribution, style, and timing of neotectonic deformation.
There is increasing awareness of the global diversity of tidal forested wetlands (TFWs) and their significance in the provision of ecosystem services. These ecosystems, including mangrove forests, tidal freshwater forested wetlands, supratidal forests and transitional forests together span multiple climatic zones, geomorphic settings, and inundation and salinity regimes. We utilise case studies across five continents to demonstrate the state of knowledge among TFWs. Intertidal mangroves are the best-defined of the TFWs thanks to decades of research on their geomorphology, hydrology and ecology across their broad distribution. Non-mangrove forest settings, however, demonstrate more diverse hydrological, biochemical and vegetation conditions. In many cases, non-mangrove forests are situated at upper intertidal or supratidal elevations, where surface waters and groundwater are subject to interactions between tides freshwater inputs. Salinity datasets show variations ranging from tidal freshwater forested wetlands and ‘low-salinity mangroves’ to mesohaline or marine salinities, often with high temporal variability. While the floristic composition of non-mangrove forests vary among biogeographic regions, locally dominant TFW species are commonly distributed beyond the tidal niche into non-tidal wetland and upland forests. This presents challenges for traditional remote sensing approaches to ecosystem mapping, which are mostly lacking for non-mangrove forests. Geomorphic approaches and developments in machine learning offer opportunities to address this.
The Chehalis basin is located between the Cascade arc and the Coast Range in southwest Washington State. It consists of sedimentary and volcanic rocks deposited throughout the Cenozoic and is underlain by the Siletzia terrane, a thick, large igneous province accreted in the Eocene. Here, we constrain evolution of the Chehalis basin depocenter using isochore maps derived from isostatic gravity anomalies, borehole data, and stratigraphy for several time periods: the Eocene (47.6–36.8 Ma), the Oligocene (36.8–20 Ma), and the Neogene (20–0 Ma). Our results suggest that local subsidence is driven by deformation on faults that bound and intersect the basin. We see northward depocenter migration, interpreted as a shift in deformation from the northwest striking Cedar Creek fault to the west striking Doty fault. We interpret these data in terms of the long-term north-south shortening and clockwise rotation of the Cascadia forearc. During the Eocene, the Cedar Creek fault was preferentially aligned with north-south shortening, but became less active when clockwise rotation brought the Doty fault into east-west alignment with the stress field. Sediment accumulation rates decreased from 196 m/Myr in the Paleogene (∼40–20 Ma) to 27 m/Myr in the Neogene (∼20–0 Ma) as Miocene uplift of the Coast Range limited accommodation space and shifted the depositional environment from marine to fluvial. Our results are consistent with the geologic evolution of the Portland and Tualatin basins to the south and reveal a depositional and structural history uniquely shaped by clockwise rotation of the Cascadia forearc.
Using a geology-based assessment methodology, the U.S. Geological Survey estimated undiscovered, technically recoverable mean conventional resources of 12 million barrels of oil and 184 billion cubic feet of gas in the Buda Limestone of Texas.
The Outer California Borderland (OCB) is an active transform plate boundary offshore Southern California, where the relationship between faulting and submarine mass transport deposits (MTDs) remains poorly understood. Onshore paleoseismic data provide high-resolution earthquake records, whereas marine geophysical data capture longer-term histories. Offshore fault systems pose hazards to infrastructure and dense coastal populations, particularly when linked to submarine landslides. We present new high-resolution geophysical data set (cruise SR2303), including bathymetric and CHIRP sub-bottom data integrated with legacy seismic reflection data and chronostratigraphic constraints from ODP Site 1012 to examine Quaternary MTD recurrence and tectonic controls in the Cortes Basin, OCB. Bathymetry shows deformational features, including slide scarps and previously unmapped fault segments with evidence of Holocene activity. CHIRP profiles reveal 10 stacked MTDs in the East Cortes Basin and 8 in the West Cortes Basin, spanning ∼752 ka with an average recurrence of ∼83.6 ± 1 ka. Acoustic imaging shows 7 MTD intervals coinciding with fault offset increments and fault growth suggesting earthquake-triggered mass wasting. A strong association between MTD occurrences and sea-level extremes also supports glacio-eustatic contribution to slope failure. Stratigraphic correlations suggest quasi-synchronous MTDs across the eastern and western areas, likely triggered by larger eathquakes in the Quaternary. Although the identified MTDs occur relatively far from the Southern California coast, they still pose a potential tsunamigenic hazard requiring further assessment. Moreover, if linked to earthquakes along major strike-slip faults, for example, the Ferrelo fault, the MTDs may provide valuable proxies to constrain rupture scenarios and fault connectivity within the understudied OCB.
The Southern California Borderland (SCB) is a topographically complex region on the active continental margin that hosts varied hardground habitats, including ferromanganese (FeMn) crusts and phosphorites, marine minerals being considered for resource extraction. The SCB is influenced by seasonal upwelling and terrestrial inputs, and has a well-defined oxygen-minimum zone (OMZ). We analyzed megafaunal community composition, density, and diversity across 41 ROV video transects at 10 SCB hard substrate sites spanning depth (378–2765 m), temperature (1.79–7.96 °C), and oxygen (3.91–105.67 μM) gradients, with varying rock types (FeMn crusts, phosphorites, other rocks). We counted 32,426 individuals representing 146 unique taxa over a total area of 21,935 m 2 . Echinodermata contributed 48 % of the total abundance; Cnidaria 24 %; Porifera 13 %; Annelida 6 %; Arthropoda 6 %; Chordata 2 %; Mollusca, Foraminifera, and Hemichordata <1 %. Megafauna communities showed clear heterogeneity, with density, diversity, and community composition varying among transects within sites and across sites. Rock type recorded no influence on megafaunal density, but habitats that included FeMn crusts had higher diversity and distinct taxonomic composition. Location, depth, and oxygen were the most important variables explaining variance among megafaunal communities, with distinct communities formed at deeper depths and within the OMZ. This study provides baseline information about the megafauna inhabiting SCB hardgrounds, particularly those associated with minerals considered for their resource potential. It extends existing documentation of such communities to greater depths, reveals increased representation of sponges and corals below 1000 m, and highlights the variability within and among deep-sea megafaunal communities, offering context for informed policy decisions.
Using a geology-based assessment methodology, the U.S. Geological Survey (USGS) estimated undiscovered, technically recoverable mean resources of 16 million barrels of oil and 348 billion cubic feet of gas in conventional reservoirs of the Norphlet Formation in the U.S. Gulf Coast region.
The Stibnite-Yellow Pine district of central Idaho was mined from the early 1900s until the 1990s, extracting gold, antimony, tungsten, and mercury from veins and disseminated and replacement ores in mountainous terrain along the headwaters of the Salmon River. Mining during the two World Wars supplied critical antimony and tungsten to the war efforts. Recent exploration has delineated mineral resources of over 187 metric tons Au, 274 metric tons Ag, and 93,000 metric tons Sb. Mineralization is hosted in Cretaceous Idaho batholith granitic rocks and a sequence of Neoproterozoic to Paleozoic metasedimentary strata of carbonate and siliciclastic compositions. Historical studies outlined some of the complex paragenesis but debated the absolute age of mineralization. New petrographic and geochronologic work documents a sequence of five hydrothermal events in the Stibnite-Yellow Pine district. Event 1 is related to Cretaceous magmatic and hydrothermal activity and includes events ranging in age from 86 to 75 Ma, including sparse quartz-molybdenite veins dated at 86 Ma. Disseminated gold mineralization of event 2 is associated with sericitic alteration and sulfidation of igneous biotite and replacement of plagioclase by potassium feldspar, largely in granodiorite. Gold is present in zoned arsenian pyrite in both disseminated ores and in crosscutting carbonate-quartz veins containing pyrite and arsenopyrite. The large Yellow Pine deposit, localized at a dilatant bend in the Meadow Creek fault, hosts such disseminated and vein gold. Event 2 is interpreted as the major gold-forming event; 40 Ar/ 39 Ar ages of sericite and potassium feldspar alteration range, respectively, from 70 to 59 and 66 to 56 Ma. The long span is interpreted to reflect the age of gold mineralization and local overprinting by event 3. A narrower range from 66 to 61 Ma is interpreted to date the peak of gold mineralization and alteration. Event 3, tungsten mineralization with scheelite, is texturally later than event 2 gold and localized along the Meadow Creek structure. Event 3 scheelite has been dated by isotope dilution-thermal ionization mass spectrometry (ID-TIMS) and laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) U-Pb methods at 57 Ma. Event 4, best developed in the West End area, includes gold-silver–bearing quartz-carbonate-pyrite veins and breccias with epithermal textures and potassium feldspar alteration envelopes. Adularia from event 4 yields 40 Ar/ 39 Ar plateau ages of 52 to 51 Ma. Event 5 antimony and mercury mineralization consists of stibnite veins and breccia cements at the Yellow Pine and Hangar Flat deposits as well as cinnabar veins and replacements at the peripheral Fern and Hermes deposits; it is constrained by an LA-ICP-MS U-Pb date on scheelite (ca. 47 Ma) intergrown with stibnite. Minor propylitic and argillic alteration is evident in 47 Ma igneous dikes, which do not contain economic mineralization. The Au-Sb-W ores in the Stibnite-Yellow Pine mining district formed over an extended time period from about 70 to 45 Ma in multiple pulses that were localized along the Meadow Creek fault zone. Each event corresponds to episodes of magmatism and/or hydrothermal activity in the region. Insignificant event 1 skarn and molybdenum mineralization is similar in age to the Thompson Creek porphyry molybdenum deposit in central Idaho. Event 2 gold mineralization occurred during a magmatic gap in central Idaho but was synchronous with magmatism in the Bitterroot lobe further north; event 2 is similar in age to orogenic gold-arsenic mineralization at the Beartrack mine in eastern Idaho. Event 3 scheelite mineralization coincides with tungsten mineralization at the Quartz Creek deposit, late magmatism in the Bitterroot lobe, and rapid exhumation of the Atlanta lobe of the Idaho batholith. Event 4 gold mineralization is coincident with the onset of regional Challis magmatism and extension. Event 5 antimony and mercury mineralization is time-equivalent to epithermal gold mineralization in the nearby Thunder Mountain volcanic field and the peak of Challis magmatism.
Using a geology-based assessment methodology, the U.S. Geological Survey estimated undiscovered, technically recoverable mean conventional resources of 13.3 billion barrels of oil and 111 trillion cubic feet of gas in the central Arabian Peninsula.
Population in western Sarpy County, Nebraska, has steadily increased over the last several decades and has led to increased groundwater use for domestic purposes. To meet the increase in demand, the Papio-Missouri River Natural Resources District is seeking to use all available sources of groundwater in western Sarpy County. Additionally, elevated groundwater nitrate plus nitrite as nitrogen concentrations were detected, indicating the need to better understand the groundwater quality of the area. Although the general geology of the area is understood, the area does not have detailed information on the extent of the various aquifers, particularly the Dakota aquifer. To characterize these aquifers, the Papio-Missouri River Natural Resources District invested in airborne electromagnetic surveys of the area to better understand the subsurface geology. Although these surveys improved understanding of the groundwater systems in the area, the Papio-Missouri River Natural Resources District wanted to integrate the subsurface information with available water-quality and groundwater-level data. In response, the U.S. Geological Survey, in cooperation with the Papio-Missouri River Natural Resources District, the Nebraska Natural Resources Commission, and the Nebraska Department of Natural Resources, assembled geologic, hydrogeologic and nitrate plus nitrite as nitrogen information for the selected area into a three-dimensional visualization computer software package called GeoScene3D. The completed GeoScene3D project was assembled to provide a visualization of the groundwater systems and associated water-quality results in Sarpy County and to provide the Papio-Missouri River Natural Resources District managers with information that can be used to make more informed groundwater resource-planning decisions in the future. This report details the development of a three-dimensional model created within GeoScene3D to visualize the subsurface, particularly the Dakota Sandstone in western Sarpy County.
The bedrock geology of the 7.5-minute Eagle Lake quadrangle, Essex County, New York, consists of deformed and metamorphosed Mesoproterozoic gneisses of the Adirondack Highlands unconformably overlain by weakly deformed lower Paleozoic sedimentary rocks of the Champlain Valley. The Mesoproterozoic rocks occur on the eastern edge of the Adirondack Highlands and represent an extension of the Grenville Province of Laurentia. Granulite facies Mesoproterozoic paragneiss, marble, and amphibolite hosted the emplacement of an anorthosite-mangerite-charnockite-granite (AMCG) suite, now exposed mostly as orthogneiss, at approximately 1.18–1.15 giga-annum (Ga, billion years before present). The earliest of four phases of deformation (D1) predated AMCG magmatism and is characterized by gneissosity, rarely preserved F1 isoclinal folds, and migmatite in the paragneiss host rocks. A sample of hornblende quartz syenite from the AMCG suite, collected from an abandoned railroad cut on Old Furnace Road, yielded a U-Pb zircon age of 1,149±10 million years before present. D2 deformation produced a composite penetrative gneissosity, migmatite, and isoclinal F2 folds. Towards the end of D2, felsic magmatism (including the regionally extensive Lyon Mountain Granite Gneiss, abbreviated “LMG”) spread by penetrative migration as semiconcordant alkali feldspar granite sheets subparallel to S2 into the previously deformed lithologies. The LMG crystallized at approximately 1.15 to 1.14 Ga and displays synkinematic F2 folds thus constraining the time of D2 deformation. Exhumation of the Marcy anorthosite began during D3 along a mylonitic extensional detachment, as a type of core complex. Protracted D3 produced F3 folds exhibited in regional domes and basins, such as the Hammondville antiform, reactivation of the S2 foliation, partial melting, metamorphism, metasomatism, iron ore remobilization, and intrusion of magnetite-bearing pegmatite both as layer-parallel sills and crosscutting dikes. D4 created NE- and NW-trending boudinage, local high-grade ductile shear zones, and crosscutting granitic pegmatite dikes. Kilometer (km)-scale lineaments readily observed in lidar data are Ediacaran mafic dikes and Phanerozoic brittle faults. Lower Paleozoic rocks are part of the Early Cambrian to Late Ordovician great American carbonate bank on the ancient margin of Laurentia. The Potsdam Sandstone preserves the Cambrian stratigraphy in outliers above the Great Unconformity. The Paleozoic rocks are weakly folded and block faulted. Parts of the quadrangle are covered by undifferentiated glacial deposits, but much of the quadrangle contains only a variably thick, veneer of unmapped glacial till over significant areas of exposed bedrock. The map also shows waste rock piles and locations of historical mining operations. This study was undertaken to improve our understanding of the bedrock geology in the Adirondack Highlands, establish a modern framework for 1:24,000-scale bedrock geologic mapping in the Adirondack Mountains, and provide a modern context for historical mines. This Scientific Investigations Map of the Eagle Lake 7.5-minute quadrangle consists of a map sheet, an explanatory pamphlet, and a geographic information system database that includes bedrock geologic units, faults, outcrops, and structural geologic information. The map sheet includes a bedrock geologic map, a correlation of map units, a description of map units, an explanation of map symbols, and two cross sections. The explanatory pamphlet includes a discussion of the geology.
Household well water can be degraded by contaminants from the land's surface, but private well owners lack means to protect the source water from neighboring disturbances. Rural residents of southwestern Wisconsin, USA, rely on private well water, and the combination of land use and fractured carbonate bedrock makes groundwater vulnerable to contamination. To identify the extent, sources, and risk factors of private well contamination, randomly selected wells sampled during two-day periods in fall (n = 301) and spring (n = 529) were analyzed for nitrate and indicator bacteria, and a subset (n = 138) was sampled across four seasonal events for analysis of pathogens and microbial source tracking markers by quantitative polymerase chain reaction. Risk factors representing land use, hydrology, geology, and well construction were analyzed for associations with contamination in multivariable models. The importance of risk factors varied by contaminant, illustrating the multifaceted nature of rural groundwater quality. Nitrate contamination was associated with agricultural land use, and wells with casings that extended below a shale aquitard accessed less contaminated water than those drawing water from above it. Human fecal microbes were detected in 64 wells (46%), and rainfall was the key risk factor for contamination, indicating that wastewater from septic systems was available to contaminate wells when transport conditions were favorable. Manure microbes from cattle/ruminants and pigs were detected in 33 and 13 wells, respectively, and concentrations increased with the hectarage of cultivated land near wells. Pathogen genes for viruses, bacteria, and protozoa were detected in 66 wells (48%), including more detections of zoonotic than human-specific pathogens, and human Bacteroides , an indicator of wastewater, was an equivocal predictor of pathogen presence in private wells. Characterizing important elements of the setting, like geology, and identifying sources and risk factors for contaminants can inform landscape-level policies to protect groundwater quality.