USGS Science⌕ Search

SEARCH · USGS Science

Results for “Northeastern Geology”

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

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

At least 847 records · Page 47Linked to original sources

Recovery and reprocessing of legacy geophysical data from the archives of the State Company of Geology and Mining (GEOSURV) of Iraq and Iraq Petroleum Company (IPC)

Aeromagnetic data belonging to the State Company of Geology and Mining of Iraq (GEOSURV) have been recovered from magnetic tapes and early paper maps. In 1974 a national airborne survey was flown by the French firm Compagnie General de Geophysique (CGG). Following the survey the magnetic data were stored on magnetic tapes within an air conditioned archive run by GEOSURV. In 1990, the power supply to the archive was cut resulting in the present‐day poor condition of the tapes. Frontier Processing Company and the U.S. Geological Survey (USGS) have been able to recover over 99 percent of the original digital data from the CGG tapes. Preliminary reprocessing of the data yielded a total magnetic field anomaly map that reveals fine structures not evident in available published maps. Successful restoration of these comprehensive, high quality digital datasets obviates the need to resurvey the entire country, thereby saving considerable time and money. These data were delivered to GEOSURV in a standard format for further analysis and interpretation. A parallel effort by GETECH concentrated on recovering the legacy gravity data from the original field data sheets archived by IPC (Iraq Petroleum Company). These data have been compiled with more recent GEOSURV sponsored surveys thus allowing for the first time a comprehensive digital and unified national gravity database to be constructed with full principal facts. Figure 1 shows the final aeromagnetic and gravity data coverage of Iraq. The only part of Iraq lacking gravity and aeromagnetic data coverage is the mountainous areas of the Kurdish region of northeastern Iraq. Joint interpretation of the magnetic and gravity data will help guide future geophysical investigations by GEOSURV, whose ultimate aim is to discover economical mineral and energy resources.

SEG Technical Program Expanded Abstracts↗

Peak ground motions and site response at Anza and Imperial Valley, California

Power spectra of shear-waves for eighteen earthquakes from the Anza-Imperial Valley region were inverted for source, mid-path Q, site attenuation and site response. The motivation was whether differences in site attenuation (parameterized as t*, r/cQ, where r is distance along ray path near the site, c is shear velocity and Q is the quality factor that parameterizes attenuation) and site response could be correlated with residuals in peak values of velocity or acceleration after removing the affect of distance-dependent attenuation. We decomposed spectra of S-waves from horizontal components of 18 earthquakes from 2010 to 2018 into a common source for each event with ω −2 spectral fall-off at high frequencies and then projected the residuals onto path and site terms following the methodology of Boatwright et al. (Bull Seismol Soc Am 81:1754–1782, 1991). The site terms were constrained to have an amplification at a particular frequency governed by V S30 at two of the sites which had downhole shear-wave logs. The 18 events, 3 < M < 4, had moments between approximately 10 20 and 10 22 dyne-cm, and stress drops between 1 and 100 bars. Average mid-crust attenuation had a Q of 844 reflecting the average path through the crystalline rock of the San Jacinto Mountains. t* for each station corresponded to the geologic environment such that stations on hard rock had low t* (e.g. stations KNW, PFO and RDM) a station in the San Jacinto fault zone (station SND) had a moderate t* of 0.035 s and stations in the Imperial Valley usually had higher t*s . Generally t* correlated with average amplification suggesting that sites characterized by low surface velocities and higher attenuation also have more amplification in the 1–6 Hz band. Residuals of peak values were determined by subtracting the prediction of Boore and Atkinson ( 2008 ). There is a correlation between average amplification and peak velocity, but not peak acceleration. Interestingly, there is less scatter at high values of amplification although there is also less data. Scatter in values of peak velocity and peak acceleration are higher at shorter compared to longer durations. When using a frequency-dependent form for Q , variances are higher, sometimes much higher; the dataset does not support frequency-dependent Q , which is not similar to results from the Imperial Valley and northeastern North America.

California↗

Low streamflow trends at human-impacted and reference basins in the United States

We present a continent-scale exploration of trends in annual 7-day low streamflows at 2482 U.S. Geological Survey streamgages across the conterminous United States over the past 100, 75, and 50 years (1916–2015, 1941–2015 and 1966–2015). We used basin characteristics to identify subsets of study basins representative of reference basins with streamflow relatively free from human effects (n = 259), and predominantly agricultural basins (n = 78), regulated basins (n = 220), and urban basins (n = 121). Trend significance was computed using the Mann-Kendall test considering short- and long-term persistence. Lag-one autocorrelation tests of detrended 7-day low streamflows for all gage classes show that time-series independence is not an appropriate assumption for annual low streamflow data at many basins. Among all study gages, upward trends (wetter conditions) in 7-day low streamflows outnumbered downward trends (drier conditions) approximately 2–1 for the 75- and 100-year trend periods—50-year trends indicated roughly equal numbers of increases and decreases. Increases in 7-day low streamflow were consistently observed for all time periods throughout much of the northeastern quadrant of the conterminous U.S. including western New England and the Mid-Atlantic, the southeastern Great Lakes basin, northern Ohio River basin, and the Upper Mississippi River and eastern Missouri River basins. Decreases in 7-day low streamflow were consistently observed for all time periods at many gages in the southeastern U.S. and in the northwestern U.S. in much of Idaho and northwestern Washington. Overall, we observed greater percentages of statistically significant trends at gages with human-induced influences than at reference gages. Low-flow trends at agricultural gages were regionally consistent with trends at reference gages. Regulated basins had many statistically significant upward trends for all three time periods tested, which may be attributed in part to substantial increases in dam-related storage prior to 1970. Urban gages had the greatest percentage of significant decreases in 7-day low flows compared to all other gage classes even though most urban gages saw upward trends in mean annual flows. Urban gages also had the greatest percentage of significant increases in low flows second only to regulated gages, highlighting that urban development can increase or decrease low streamflows depending on the basin-specific development.

Conterminous United States↗

Regional magnetic anomalies, crustal strength, and the location of the northern Cordilleran fold-and-thrust belt

The northern Cordilleran fold-and-thrust belt in Canada and Alaska is at the boundary between the broad continental margin mobile belt and the stable North American craton. The fold-and-thrust belt is marked by several significant changes in geometry: cratonward extensions in the central Yukon Territory and northeastern Alaska are separated by marginward re-entrants. These geometric features of the Cordilleran mobile belt are controlled by relations between lithospheric strength and compressional tectonic forces developed along the continental margin. Regional magnetic anomalies indicate deep thermal and compositional characteristics that contribute to variations in crustal strength. Our detailed analysis of one such anomaly, the North Slope deep magnetic high, helps to explain the geometry of the fold-and-thrust front in northern Alaska. This large magnetic anomaly is inferred to reflect voluminous mafic magmatism in an old (Devonian?) extensional domain. The presence of massive amounts of malic material in the lower crust implies geochemical depletion of the underlying upper mantle, which serves to strengthen the lithosphere against thermal erosion by upper mantle convection. We infer that deep-source magnetic highs are an important indicator of strong lower crust and upper mantle. This stronger lithosphere forms buttresses that play an important role in the structural development of the northern Cordilleran fold-and-thrust belt. ?? 2007 The Geological Society of America.

Geology↗

Preliminary metallogenic belt and mineral deposit maps for northeast Asia

The metallogenic belts and locations of major mineral deposits of Northeast Asia are portrayed on Sheets 1-4. Sheet 1 portrays the location of significant lode deposits and placer districts at a scale of 1:7,500,000. Sheets 2-4 portray the metallogenic belts of the region in a series of 12 time-slices from the Archean through the Quaternary at a scale of 1:15,000,000. For all four map sheets, a generalized geodynamics base map, derived from a more detailed map by Parfenov and others (2003), is used as an underlay for the metallogenic belt maps. This geodynamics map underlay permits depicts the major host geologic units and structures that host metallogenic belts. Four tables are included in this report. A hierarchial ranking of mineral deposit models is listed in Table 1. And summary features of lode deposits, placer districts, and metallogenic belts are described in Tables 2, 3, and 4, respectively. The metallogenic belts for Northeast Asia are synthesized, compiled, described, and interpreted with the use of modern concepts of plate tectonics, analysis of terranes and overlap assemblages, and synthesis of mineral deposit models. The data supporting the compilation are: (1) comprehensive descriptions of mineral deposits; (2) compilation and synthesis of a regional geodynamics map the region at 5 million scale with detailed explanations and cited references; and (3) compilation and synthesis of metallogenic belt maps at 15 million scale with detailed explanations and cited references. These studies are part of a major international collaborative study of the Mineral Resources, Metallogenesis, and Tectonics of Northeast Asia that is being conducted from 1997 through 2002 by geologists from earth science agencies and universities in Russia, Mongolia, Northeastern China, South Korea, Japan, and the USA. Companion studies and previous publications are: (1) a detailed geodynamics map of Northeast Asia (Parfenov and 2003); (2) a compilation of major mineral deposit models (Rodionov and Nokleberg, 2000; Rodionov and others, 2000; Obolenskiy and others, 2003); and (3) a database on significant metalliferous and selected nonmetalliferous lode deposits, and selected placer districts (Ariunbileg and others, 2003).

Open-File Report↗

Stable-isotope and solute-chemistry approaches to flow characterization in a forested tropical watershed, Luquillo Mountains, Puerto Rico

The prospect of changing climate has led to uncertainty about the resilience of forested mountain watersheds in the tropics. In watersheds where frequent, high rainfall provides ample runoff, we often lack understanding of how the system will respond under conditions of decreased rainfall or drought. Factors that govern water supply, such as recharge rates and groundwater storage capacity, may be poorly quantified. This paper describes 8-year data sets of water stable isotope composition ( δ 2 H and δ 18 O) of precipitation (4 sites) and a stream (1 site), and four contemporaneous stream sample sets of solute chemistry and isotopes, used to investigate watershed response to precipitation inputs in the 1780-ha Río Mameyes basin in the Luquillo Mountains of northeastern Puerto Rico. Extreme δ 2 H and δ 18 O values from low-pressure storm systems and the deuterium excess ( d -excess) were useful tracers of watershed response in this tropical system. A hydrograph separation experiment performed in June 2011 yielded different but complementary information from stable isotope and solute chemistry data. The hydrograph separation results indicated that 36% of the storm rain that reached the soil surface left the watershed in a very short time as runoff. Weathering-derived solutes indicated near-stream groundwater was displaced into the stream at the beginning of the event, followed by significant dilution. The more biologically active solutes exhibited a net flushing behavior. The d- excess analysis suggested that streamflow typically has a recent rainfall component (∼25%) with transit time less than the sampling resolution of 7 days, and a more well-mixed groundwater component (∼75%). The contemporaneous stream sample sets showed an overall increase in dissolved solute concentrations with decreasing elevation that may be related to groundwater inputs, different geology, and slope position. A considerable amount of water from rain events runs off as quickflow and bypasses subsurface watershed flowpaths, and better understanding of shallow hillslope and deeper groundwater processes in the watershed will require sub-weekly data and detailed transit time modeling. A combined isotopic and solute chemistry approach can guide further studies to a more comprehensive model of the hydrology, and inform decisions for managing water supply with future changes in climate and land use.

Applied Geochemistry↗

Groundwater flow conceptualization of the Pahute Mesa–Oasis Valley Groundwater Basin, Nevada—A synthesis of geologic, hydrologic, hydraulic-property, and tritium data

This report provides a groundwater-flow conceptualization that integrates geologic, hydrologic, hydraulic-property, and radionuclide data in the Pahute Mesa–Oasis Valley (PMOV) groundwater basin, southern Nevada. Groundwater flow in the PMOV basin is of interest because 82 underground nuclear tests were detonated, most near or below the water table. A potentiometric map and nine sets of hydrostratigraphic and hydrologic cross sections supplement the conceptualization. Potentiometric contours indicate that groundwater in the PMOV basin generally flows south-southwest and discharges at Oasis Valley. Groundwater encounters an alternating sequence of low- and high-transmissivity rocks, referred to as dams and pools, respectively, as it moves from east to west across eastern Pahute Mesa. Flow from all Pahute Mesa nuclear tests is to Oasis Valley and is well-constrained by water-level data. Flow converges along a corridor of high transmissivity between Pahute Mesa and Oasis Valley. The location of the lateral PMOV basin boundary is well defined, and this boundary, with a few minor exceptions, represents a no-flow boundary. Some boundary uncertainty exists in the northeastern part of the basin, but potential flow-rate estimates across the northeastern boundary resulting from this uncertainty are small relative to the basin groundwater budget. Recharge in the PMOV basin is derived from episodic pulses of modern water and the diffuse percolation of old water (greater than 1,000 years). Episodic recharge is a minor recharge component observed as a rise in groundwater levels that occurs 3 months to 1 year following a wet winter. Minor amounts of episodic recharge through an unsaturated zone in excess of 1,000 feet (ft) requires preferential flow through faults and fractures. The dominant recharge component is slow, steady, diffuse percolation of old water through the unsaturated zone. A large component of old water recharging the groundwater system is consistent with observations of isotopically light deuterium and oxygen 18 compositions in water from wells on Pahute Mesa and central Oasis Valley. About half the recharge in the PMOV basin is derived from the eastern Pahute Mesa area. The remaining recharge is derived primarily from other highland areas including Timber Mountain, Belted and Kawich Ranges, and Black Mountain. The PMOV groundwater system is nearly steady state, where recharge is balanced by the 5,900 acre-feet per year of natural discharge at Oasis Valley. This assumption is reasonable because the basin is dominated by steady-state conditions, where long-term changes in groundwater storage are minimal. Total groundwater withdrawals from 1963 to 2018 have amounted to less than 10 percent of annual groundwater discharge and less than 0.2 percent of the basin’s groundwater storage. Therefore, present-day (2020) conditions are considered representative of predevelopment (pre-1950) conditions in nearly all areas of the basin. The lower PMOV basin boundary is defined at 4,000 ft below the water table to encompass all underground nuclear tests and tritium plumes. This boundary defines the lower boundary of radionuclide migration. However, nearly all flow and tritium transport occur in the upper 1,600 ft of the saturated zone because a transmissivity-with-depth relation indicates that greater than 90 percent of the transmissivity contributing to groundwater flow occurs within 1,600 ft of the water table. Rocks at deeper depths have low transmissivity because argillic and mineralized alterations plug the fractures. Volcanic rocks form the primary aquifers and confining units in the PMOV basin. Volcanic hydrogeologic units (HGUs) and hydrostratigraphic units (HSUs) have transmissivity distributions that span up to eight orders of magnitude with considerable overlap between distributions. Despite the large overlap between units, mean transmissivities of aquifers are one-to-two orders of magnitude greater than the confining units. However, all volcanic-rock HGUs and HSUs are composite units, meaning that they can function spatially as either an aquifer or confining unit .

Nevada↗

Stable isotope study of fluid inclusions in fluorite from Idaho: Implications for continental climates during the Eocene

Isotopic studies of fluid inclusions from meteoric water-dominated epithermal ore deposits offer a unique opportunity to study paleoclimates because the fluids can provide direct samples of ancient waters. The oxygen and hydrogen isotope compositions of meteoric waters vary because of changes in climatic variables such as mean annual temperature of precipitation, relative humidity, origin and history of air masses, and the isotope composition of the oceans. Inclusion fluids found in fluorite (CaF 2 ) are especially useful because their host is devoid of oxygen or hydrogen, thus precluding postentrapment isotope exchange. Fluorite-hosted fluid inclusions from the Eocene (51-50 Ma) epithermal deposits of the Bayhorse mining district, northeastern Idaho, have low salinities, most less than 0.6 equivalent wt% NaCl, and low to moderate homogenization temperatures (98 to 146 °C), indicating meteoric origins for the fluids. Oxygen and hydrogen isotope data on inclusion fluids are almost identical to those of modern meteoric waters in the area. The equivalence of the isotope composition of the Eocene inclusion fluids and modern meteoric waters indicates that the Eocene climatic conditions were similar to those today. This conclusion supports the climate modeling of Sloan and Barron, who suggested that the climates of continental interiors do not reflect the magnitude of warming preserved by the deep-ocean paleoclimate record during the Eocene.

Geology↗

Comparative mineral chemistry and textures of SAFOD fault gouge and damage-zone rocks

Creep in the San Andreas Fault Observatory at Depth (SAFOD) drillhole is localized to two foliated gouges, the central deforming zone (CDZ) and southwest deforming zone (SDZ). The gouges consist of porphyroclasts of serpentinite and sedimentary rock dispersed in a foliated matrix of Mg-smectite clays that formed as a result of shearing-enhanced reactions between the serpentinite and quartzofeldspathic rocks. The CDZ takes up most of the creep and exhibits differences in mineralogy and texture from the SDZ that are attributable to its higher shearing rate. In addition, a ∼0.2-m-wide sector of the CDZ at its northeastern margin (NE-CDZ) is identical to the SDZ and may represent a gradient in creep rate across the CDZ. The SDZ and NE-CDZ have lower clay contents and larger porphyroclasts than most of the CDZ, and they contain veinlets and strain fringes of calcite in the gouge matrix not seen elsewhere in the CDZ. Matrix clays in the SDZ and NE-CDZ are saponite and corrensite, whereas the rest of the CDZ lacks corrensite. Saponite is younger than corrensite, reflecting clay crystallization under declining temperatures, and clays in the more actively deforming portions of the CDZ have better equilibrated to the lower-temperature conditions.

California↗

Chemical constituent concentrations in stream water, streambed sediment, and soils of Fort Belvoir, Virginia—A characterization of ambient conditions in 2019

Introduction The U.S. Army Fort Belvoir (FTBL) installation is on the banks of the Potomac River in Fairfax County, northeastern Virginia. The installation was founded by the U.S. Army during World War I. It has been home to a variety of military organizations over the course of its more than 100-year history and currently houses more than 145 mission partners. The installation consists of two noncontiguous units, the Main Post, and a smaller area to the northwest, Fort Belvoir North Area (FTNA). FTBL encompasses 8.91 square miles. There is concern that activities on FTBL, including a long history of training, operations, and maintenance, may have resulted in contamination of stream water, streambed sediment, and (or) soils. Of particular concern is the U.S. Environmental Protection Agency (EPA) Target Analyte List (TAL). TAL refers to “the list of inorganic compounds/elements designated for analysis as contained in the version of the EPA Contract Laboratory Program Statement of Work for Inorganics Analysis, Multi-Media, Multi-Concentration in effect as of the date on which the laboratory is performing the analysis” ( https://www.nj.gov/dep/srp/guidance/tcl_tal/ ). Because of the potential for TAL contamination at FTBL, the U.S. Geological Survey (USGS), in cooperation with U.S. Army Fort Belvoir, conducted a survey of FTBL’s stream water, streambed sediment, and soils during calendar year 2019. The terminology “ambient concentrations” is used in this report to represent the concentrations of the TAL and other constituents at the time of sampling. This is in contrast to “background concentrations,” a term that “refers to areas in which the concentrations of chemicals have not been elevated by site activities”. Although some of the samples collected for this project may represent “background concentrations,” there is no assurance that they do, so all data collected are described as having “ambient concentrations.” The purpose of the study was to obtain environmental data to characterize ambient concentrations of EPA TAL constituents in stream water, streambed sediment, and soils in FTBL, Virginia. This report describes methods and results of sampling stream water, streambed sediment, and soils during 2019. The purpose of this report is four-fold: (1) to describe the field sampling methods used to collect stream water, streambed sediment, and soils; (2) to describe the laboratory methods used to analyze the samples; (3) to report summaries of the field and laboratory results; and (4) to report the quality assurance and quality control results.

Virginia↗

Selected metals in sediments and streams in the Oklahoma Part of the Tri-State Mining District, 2000-2006

The abandoned Tri-State mining district includes 1,188 square miles in northeastern Oklahoma, southeastern Kansas, and southwestern Missouri. The most productive part of the Tri-State mining district was the 40-square mile part in Oklahoma, commonly referred to as 'the Picher mining district' in north-central Ottawa County, Oklahoma. The Oklahoma part of the Tri-State mining district was a primary producing area of lead and zinc in the United States during the first half of the 20th century. Sulfide minerals of cadmium, iron, lead, and zinc that remained in flooded underground mine workings and in mine tailings on the land surface oxidized and dissolved with time, forming a variety of oxide, hydroxide, and hydroxycarbonate metallic minerals on the land surface and in streams that drain the district. Metals in water and sediments in streams draining the mining district can potentially impair the habitat and health of many forms of aquatic and terrestrial life. Lakebed, streambed and floodplain sediments and/or stream water were sampled at 30 sites in the Oklahoma part of the Tri-State mining district by the U.S. Geological Survey and the Oklahoma Department of Environmental Quality from 2000 to 2006 in cooperation with the U.S. Environmental Protection Agency, and the Quapaw and Seneca-Cayuga Tribes of Oklahoma. Aluminum and iron concentrations of several thousand milligrams per kilogram were measured in sediments collected from the upstream end of Grand Lake O' the Cherokees. Manganese and zinc concentrations in those sediments were several hundred milligrams per kilogram. Lead and cadmium concentrations in those sediments were about 10 percent and 0.1 percent of zinc concentrations, respectively. Sediment cores collected in a transect across the floodplain of Tar Creek near Miami, Oklahoma, in 2004 had similar or greater concentrations of those metals than sediment cores collected at the upstream end of Grand Lake O' the Cherokees. The greatest concentrations of cadmium, iron, lead, and zinc were detected in sediments beneath an intermittent tributary to Tar Creek, a slough which drains mined areas near Commerce, Oklahoma. In surface water, aluminum and iron concentrations were greatest in the Neosho River, perhaps a result of runoff from areas underlain by shales. The greatest aqueous concentrations of cadmium, lead, manganese, and zinc were measured in water from Tar Creek, the primary small stream draining the Picher mining district with the largest proportion of mined area. Water from the Spring River had greater zinc concentrations than water from the Neosho River, perhaps as a result of a greater proportion of mined area in the Spring River Basin. Dissolved metals concentrations were generally much less than total metals concentrations, except for manganese and zinc at sites on Tar Creek, where seepage of ground water from the mine workings, saturated mine tailings, and/or metalliferous streambed sediments may be sources of these dissolved metals. Iron and lead concentrations generally decreased with increasing streamflow in upstream reaches of Tar Creek, indicating dilution of metals-rich ground water by runoff. Farther downstream in Tar Creek, and in the Neosho and Spring Rivers, metals concentrations tended to increase with increasing streamflow, indicating that most metals in these parts of these streams were associated with runoff and re-suspension of metals precipitated as oxide, hydroxide, and hydroxycarbonate minerals on land surface and streambeds. Estimated total aluminum, cadmium, iron, manganese, and zinc loads generally were greatest in water from the Neosho and Spring Rivers, primarily because of comparatively large streamflows in those rivers. Slight increases in metal loads in the downstream directions on those rivers indicated contributions of metals from inflows of small tributaries such as Tar Creek and from runoff.

Kansas, Missouri, Oklahoma↗

Computer program for simulation of variable recharge with the U. S. Geological Survey modular finite-difference ground-water flow model (MODFLOW)

The Variable-Recharge Package is a computerized method designed for use with the U.S. Geological Survey three-dimensional finitedifference ground-water flow model (MODFLOW-88) to simulate areal recharge to an aquifer. It is suitable for simulations of aquifers in which the relation between ground-water levels and land surface can affect the amount and distribution of recharge. The method is based on the premise that recharge to an aquifer cannot occur where the water level is at or above land surface. Consequently, recharge will vary spatially in simulations in which the Variable- Recharge Package is applied, if the water levels are sufficiently high. The input data required by the program for each model cell that can potentially receive recharge includes the average land-surface elevation and a quantity termed ?water available for recharge,? which is equal to precipitation minus evapotranspiration. The Variable-Recharge Package also can be used to simulate recharge to a valley-fill aquifer in which the valley fill and the adjoining uplands are explicitly simulated. Valley-fill aquifers, which are the most common type of aquifer in the glaciated northeastern United States, receive much of their recharge from upland sources as channeled and(or) unchanneled surface runoff and as lateral ground-water flow. Surface runoff in the uplands is generated in the model when the applied water available for recharge is rejected because simulated water levels are at or above land surface. The surface runoff can be distributed to other parts of the model by (1) applying the amount of the surface runoff that flows to upland streams (channeled runoff) to explicitly simulated streams that flow onto the valley floor, and(or) (2) applying the amount that flows downslope toward the valley- fill aquifer (unchanneled runoff) to specified model cells, typically those near the valley wall. An example model of an idealized valley- fill aquifer is presented to demonstrate application of the method and the type of information that can be derived from its use. Documentation of the Variable-Recharge Package is provided in the appendixes and includes listings of model code and of program variables. Comment statements in the program listings provide a narrative of the code. Input-data instructions and printed model output for the package are included.

Open-File Report↗

Changes in phosphorus and suspended solids loading in the Fox River, northeastern Wisconsin, 1989–2021

The entire Lower Fox River and inner bay of Green Bay, in northeastern Wisconsin, have been listed as impaired by the Wisconsin Department of Natural Resources (WDNR) for low dissolved oxygen and degraded habitat, with total phosphorus (TP) and total suspended solids (TSS) concentrations listed as the likely causes of these impairments. To restore the Fox River and Green Bay, total maximum daily loads (TMDLs) were developed for TP and TSS, and actions were taken throughout the Fox River Basin to improve water quality. In this study, we estimated concentrations and loads of TP, dissolved phosphorus (DP), and TSS at the Lake Winnebago outlet, De Pere, and the mouth of the Fox River from water year (WY) 1989 to WY 2021; described changes in concentrations and loads through time during this period; and compared the concentrations and loads for the most recent 5-year period (WYs 2017–21) with the WDNR criteria for TP impairment and the TMDL loading goals. TP, DP, TSS, and total suspended sediment concentration data were obtained from NEW Water (the brand of the Green Bay Metropolitan Sewerage District), the WDNR, and the U.S. Geological Survey and combined into one dataset. All the TSS and total suspended sediment data were used together with no adjustment factor and are referred to as simply “TSS.” During WYs 1989–2021, mean annual TP concentrations increased from 0.089 milligram per liter (mg/L) at the Lake Winnebago outlet to 0.128 mg/L at the mouth of the Fox River, and concentrations decreased at all three sites from WY 1989 to WY 2021. The most recent (WYs 2017–21) median May–October TP concentrations were just less than the 0.1-mg/L WDNR criterion for TP impairment at the two upstream sites (Lake Winnebago outlet and De Pere) but were slightly greater than the criterion for impairment at the mouth of the Fox River. Mean annual DP concentrations increased from 0.024 mg/L at the Lake Winnebago outlet to 0.036 mg/L at the mouth of the Fox River. DP concentrations increased from WY 1989 to WY 2021 at the Lake Winnebago outlet but not at the other sites. Mean annual TSS concentrations increased from 13.5 mg/L at the Lake Winnebago outlet to 23.9 mg/L at the mouth of the Fox River and have decreased at all three sites from WY 1989 to WY 2021. The recent median May–October TSS concentrations were less than the 20-mg/L WDNR criterion for impairment at all three sites. Streamflow and TP, DP, and TSS loads increased from the Lake Winnebago outlet to the mouth of the Fox River (TP loads increased from 360,000 to 557,000 kilograms per year [kg/yr], DP loads increased from 114,000 to 162,000 kg/yr, and TSS loads increased from 60,400 metric tons per year [t/yr] to 122,600 t/yr). At the Lake Winnebago outlet, DP concentrations and TP and DP loads increased from WY 1989 to WY 2021 because of an increase in DP concentrations in Lake Winnebago resulting from the lake becoming nitrogen limited as a result of biological processes not consuming the DP in the lake and an increase in streamflow leaving the lake. Although TP and TSS concentrations decreased at De Pere and the mouth of the Fox River, there was little change in the loading because of an increase in flow. Flow-normalized TP and TSS loads at De Pere and the mouth of the Fox River decreased possibly because of implementation of agricultural conservation management practices, reductions in point-source discharges in its drainage basin, and deposition of sediment and phosphorus in recently dredged areas of the Lower Fox River. Additional studies are needed to determine the relative importance of each of these actions and whether the decrease in concentrations and flow-normalized loads will continue to be observed in the Fox River.

Wisconsin↗

Late Cenozoic tephrochronology of the Mount Diablo area within the evolving plate-tectonic boundary zone of northern California

We present a tephrochronologic/chronostratigraphic database for the Mount Diablo area and greater San Francisco Bay region that provides a spatial and temporal framework for geologic studies in the region, including stratigraphy, paleogeography, tectonics, quantification of earth surface processes, recurrence of natural hazards, and climate change. We identified and correlated 34 tephra layers within this region using the chemical composition of their volcanic glasses, stratigraphic sequence, and isotopic and other dating techniques. Tephra layers range in age from ca. 65 ka to ca. 29 Ma, as determined by direct radiometric techniques or by correlation to sites where they have been dated. The tephra layers are of Quaternary or Neogene age except for two that are of Oligocene age. We correlated the tephra layers among numerous sites throughout northern California. Source areas of the tephra layers are the Snake River–Yellowstone hotspot trend of northern Nevada, southern Idaho, and western Wyoming; the Nevadaplano caldera complex of central Nevada; the Jemez Mountains–Valles Caldera in northwestern New Mexico; the Southern Nevada volcanic field and related source areas in eastern California and west-central Nevada; the Quien Sabe–Sonoma volcanic centers of the California Coast Ranges; and the young Cascade Range volcanic centers of northeastern California and Oregon.

GSA Memoirs↗

An index of biological integrity for northern Mid-Atlantic Slope drainages

An index of biological integrity (IBI) was developed for streams in the Hudson, Delaware, and Susquehanna River drainages in the northeastern United States based on fish assemblage data from the Mohawk River drainage of New York. The original IBI, developed for streams in the U.S. Midwest, was modified to reflect the assemblage composition and structure present in Mid-Atlantic Slope drainages. We replaced several of the Midwestern IBI metrics and criteria scores because fishes common to the Midwest are absent from or poorly represented in the Northeast and because stream fish assemblages in the Northeast are less rich than those in the Midwest. For all replacement metrics we followed the ecology-based rationale used in the development of each of the metrics of the Midwestern IBI so that the basic theoretical underpinnings of the IBI remained unchanged. The validity of this modified IBI is demonstrated by examining the quality of streams in the Hudson, Delaware, and lower Susquehanna River basins. The relationships between the IBI and other indicators of environmental quality are examined using data on assemblages of fish and benthic macroinvertebrates and on chemical and physical stream characteristics obtained during 1993-2000 by the U.S. Geological Survey's National Water Quality Assessment Program in these three river basins. A principal components analysis (PCA) of chemical and physical variables from 27 sites resulted in an environmental quality gradient as the primary PCA axis (eigenvalue, 0.41 ). Principal components analysis site scores were significantly correlated with such benthic macroinvertebrate metrics as the percentage of Ephemeroptera, Plecoptera, and Trichoptera taxa (Spearman R = -0.66, P < 0.001). Index of biological integrity scores for sites in these three river basins were significantly correlated with this environmental quality gradient (Spearman R = -0.78, P = 0.0001). The northern Mid-Atlantic Slope IBI appears to be sensitive to environmental degradation in all three of the river basins addressed in this study. Adjustment of metric scoring criteria may be warranted, depending on composition of fish species in streams in the study area and on the relative effort used in the collection of fish assemblage data.

Transactions of the American Fisheries Society↗

Influence of groundwater on distribution of dwarf wedgemussels (Alasmidonta heterodon) in the upper reaches of the Delaware River, northeastern USA

The remaining populations of the endangered dwarf wedgemussel (DWM) ( Alasmidonta heterodon ) in the upper Delaware River, northeastern USA, were hypothesized to be located in areas of greater-than-normal groundwater discharge to the river. We combined physical (seepage meters, monitoring wells and piezometers), thermal (fiber-optic distributed temperature sensing, infrared, vertical bed-temperature profiling), and geophysical (electromagnetic-induction) methods at several spatial scales to characterize known DWM habitat and explore this hypothesis. Numerous springs were observed using visible and infrared imaging along the river banks at all three known DWM-populated areas, but not in adjacent areas where DWM were absent. Vertical and lateral groundwater gradients were toward the river along all three DWM-populated reaches, with median upward gradients 3 to 9 times larger than in adjacent reaches. Point-scale seepage-meter measurements indicated that upward seepage across the riverbed was faster and more consistently upward at DWM-populated areas. Discrete and areally distributed riverbed-temperature measurements indicated numerous cold areas of groundwater discharge during warm summer months; all were within areas populated by DWM. Electromagnetic-induction measurements, which may indicate riverbed geology, showed patterning but little correlation between bulk streambed electromagnetic conductivity and areal distribution of DWM. In spite of complexity introduced by hyporheic exchange, multiple lines of research provide strong evidence that DWM are located within or directly downstream of areas of substantial focused groundwater discharge to the river. Broad scale thermal-reconnaissance methods (e.g., infrared) may be useful in locating and protecting other currently unknown mussel populations.

Delaware River↗

Small plate tectonics in the northeastern Pacific

Lithospheric plate motions in the northeastern Pacific were complicated at about 2.5 m.y. B.P. by the movement along a major northeast-trending fault cutting Cascadia Basin. An estimate of the slip rate along this fault gives critical information on the relative motions of four geometrically interdependent blocks. The fault is presently inactive. Seventy km of slip along this fault during 2 m.y. or less gives an average slip rate of about 3-5 cm/yr or greater, and resulting plate motions suggest a significantly greater rate of net subduction along the continental margin off Oregon than off Washington and Vancouver Island. Subduction rate off Oregon is less sensitive to slip rate along this fault than is subduction off Washington.

British Columbia, Oregon, Washington↗

Maps of elevation of top of Pierre Shale and surficial deposit thickness with hydraulic properties from borehole geophysics and aquifers tests within and near Ellsworth Air Force Base, South Dakota, 2020–21

The U.S. Geological Survey, in cooperation with the U.S. Air Force Civil Engineer Center, collected borehole geophysical data and completed simple aquifer tests to estimate the thickness and hydraulic properties of surficial deposits. The purpose of data collection was to create generalized contour maps of Pierre Shale elevation and surficial deposit thickness within and near Ellsworth Air Force Base (study area). Natural gamma and electromagnetic induction data were collected to refine or determine surficial deposit thickness at selected wells. Additionally, data from previous geophysical studies and driller logs were compiled and combined with results from natural gamma and electromagnetic induction data to provide a more spatially complete image of the subsurface. Borehole nuclear magnetic resonance (bNMR) data were collected to estimate hydraulic conductivity and water content of surficial deposits overlying Pierre Shale. Simple aquifer tests using water slugs (slug tests) were completed to estimate hydraulic conductivity of surficial deposits, and results were compared to hydraulic conductivity estimates from bNMR data. All data used to construct maps and estimate hydraulic properties are provided in an accompanying U.S. Geological Survey data release (available at https://doi.org/10.5066/P9FLR79F ). Generalized contour maps were constructed using results from 26 borehole geophysical logs, 35 geophysical transects from previous studies, and 304 wells with driller logs. Pierre Shale elevation generally followed land-surface topography, sloping from high elevations in the north to lower elevations in the south. Topographic highs of Pierre Shale, where present, could act as groundwater divides that potentially affect groundwater flow direction. Surficial deposit thickness varied spatially and ranged from 0 to 86 feet. Surficial deposits generally were thickest in higher elevation areas near ephemeral streams in the northern part of the study area. Hydraulic conductivity estimated from bNMR results using two analytical methods ranged from 0.1 to 2,314 feet per day, whereas hydraulic conductivity estimated from slug tests ranged from 0.001 to 193 feet per day. Hydraulic conductivity estimates from slug tests were plotted with surficial deposit thickness contours instead of bNMR estimates because bNMR estimates were determined to overestimate hydraulic conductivity. Hydraulic conductivity values generally were greater in the southwestern part of study area than the northeastern part.

South Dakota↗