USGS Science⌕ Search

SEARCH · USGS Science

Results for “Final Report”

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

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

At least 1,585 records · Page 88Linked to original sources

No detection of Batrachochytrium salamandrivorans after non-native salamander putative pet release in Massachusetts, USA

The disease caused by the amphibian chytrid fungus Batrachochytrium salamandrivorans (Bsal) has caused mass mortality of salamanders across Europe (Stegen et al. 2017; amphibians are not known to be susceptible to the disease (Martel et al. 2014). Bsal puts North American amphibian biodiversity at risk (Martel et al. 2014; Gray et al. 2015) despite imports of potentially infected amphibians (USFWS 2016; WAPTR 2018; Grear et al. 2021). Recently, the U.S. has announced the finalization of the 2016 interim rule to the Lacey Act (USFWS 2016) to restrict imports of all species of salamanders from 20 genera to the list of injurious amphibians and has added 16 genera to the injurious amphibian list in a new interim rule in response to updated research on salamander susceptibility to Bsal infection (USFWS 2025). A large-scale survey across the USA failed to detect Bsal in over 10,000 individuals (Waddle et al. 2020), and ongoing efforts by North American colleges and universities are similarly conducting Bsal surveillance without more recent detection (e.g., Hyman and Hill 2024). Concerned amphibian and disease specialists nevertheless are preparing for a Bsal invasion into North America (Grant et al. 2017; Olson et al.2024). Despite the initial survey results, the risk of Bsal incursion into North America remains due to broad host susceptibility of North American amphibian species (Gray et al. 2023). The major driver of Bsal introduction risk in the U.S. is expected to be captive amphibian release (Cunningham et al. 2015; Grear et al. 2021). Bsal primarily affects salamanders but anuran species can also be vectors (e.g., Towe et al. 2021; Monzon et al. 2022; Gray et al. 2023; Olson et al. 2024), increasing the risk of a released amphibian carrying Bsal into the natural environment (Grear et al. 2021). An early survey of pet salamanders in the United States did not detect Bsal (Klocke et al. 2017), but even one animal released with Bsal could have a major effect on local amphibians. Early detection and rapid response can improve host outcomes (Bletz et al. 2024; Gray et al. 2018), but the lack of timely knowledge about unwanted release of non-native pets to natural systems causes inevitable delay. Herein, we describe our experience responding to a report of a released non-native amphibian species, common in the pet trade (Gerson 2012; Maximo et al. 2021). We suggest improvements to early detection and rapid response efforts for amphibians and amphibian disease

Massachusetts↗

Conversion and comparison of the mathematical, three-dimensional, finite-difference, ground-water flow model to the modular, three-dimensional, finite-difference, ground-water flow model for the Tesuque aquifer system in northern New Mexico

The three-dimensional finite-difference groundwater model (using a mathematical groundwater flow code) of the Tesuque aquifer system in northern New Mexico was converted to run using the U.S. Geological Survey 's modular groundwater flow code. Results from the final versions of the predevelopment and 1947 to 2080 transient simulations of the two models are compared. A correlation coefficient of 0.9905 was obtained for the match in block-by-block head-dependent fluxes for predevelopment conditions. There are, however, significant differences in at least two specific cases. In the first case, a difference is associated with the net loss from the Pojoaque River and its tributaries to the aquifer. The net loss by the river is given as 1.134 cu ft/sec using the original groundwater model, which is 38.1% less than the net loss by the river of 1.8319 cu ft/sec computed in this study. In the second case, the large difference is computed for the transient decline in the hydraulic head of a model block near Tesuque Pueblo. The hydraulic-head decline by 2080 is, using the original model, 249 ft, which is 14.7% less than the hydraulic head of 292 ft computed by this study. In general, the differences between the two sets of results are not large enough to lead to different conclusions regarding the behavior of the system at steady state or when pumped. (USGS)

Open-File Report↗

A study of global sand seas

The birth of the idea that led to this publication on "Global Sand Seas" dates back to the late 1920's. At that time I was engaged in a study of the Coconino Sandstone of Arizona's Grand Canyon. Considerable controversy existed then as to whether this sandstone was a subaqueous deposit or was composed of wind-formed dunes. It became apparent that definitive literature was sparse or lacking on types of dunes, global distribution of these types, the mechanics of their development, the precise nature of their internal structure of cross-stratificiation, and the relation of wind systems to these sand forms. Especially lacking were data on criteria that could confidently be used in the recognition of ancient dunes. The common denominator in this publication is eolian sand bodies. Although the book is concerned primarily with desert sand seas, the subject matter is not restricted to deserts; it includes many references to deposits of coastal sand and to sand bodies in humid climates. Nor does the book deal exclusively with dunes, which, according to most definitions, involve mounds or hills. Many references are made to sand sheets, sand stringers, and other types of sand deposits that have no prominent topographic expression. All sand bodies accumulated by the action of wind are discussed. Chapters A-J of this publication are primarily topical. Chapters cover the grain texture, the color, and the structure of modern dunes and other eolian sands. Special treatment is given to the relation of wind data to dune interpretation, the evolution of form in current-deposited sand bodies as determined from experimental studies, and the discriminant analysis technique for differentiating between coastal and inland desert sands. This topical part of the publication also includes an analysis of criteria used in ancient deposits to interpret their eolian genesis and a consideration of economic application of the principles described, including a discussion of potentials and problems associated with eolian hydrocarbon reservoirs. The final chapters present a discussion of the morphology and distribution of dunes as determined largely from Landsat images. Chapter K of the publication is devoted to descriptions of major sand seas based largely on thematic maps derived from Landsat (ERTS) mosaics. Although inclusion herein of the actual mosaics proved to be impractical, the maps derived from them do show the distribution and abundance of various dune types and the relations of these types to certain associated features, such as bedrock, water bodies, and juxtaposed dunes. Furthermore, sand roses included with each of these maps enable the user to draw conclusions on the probable relations of wind strength and direction to dune type in a particular area. Regional studies (chapter K) were a team effort. Analysis of the Landsat (ERTS) mosaics and mapping boundaries of individual dune types were by Carol Breed. Synthesis of the rather voluminous literature and preparation of abstracts covering it was by Camilla MacCauley. Actual preparation of maps was by Franci Lennartz and later by Sarah Andrews. The gathering of data on wind, the calculation of wind roses, and the interpretation of their relations to sand bodies were by Steven Fryberger, assisted by Gary Dean.

Professional Paper↗

Net-infiltration map of the Navajo Sandstone outcrop area in western Washington County, Utah

As populations grow in the arid southwestern United States and desert bedrock aquifers are increasingly targeted for future development, understanding and quantifying the spatial variability of net infiltration and recharge becomes critically important for inventorying groundwater resources and mapping contamination vulnerability. A Geographic Information System (GIS)-based model utilizing readily available soils, topographic, precipitation, and outcrop data has been developed for predicting net infiltration to exposed and soil-covered areas of the Navajo Sandstone outcrop of southwestern Utah. The Navajo Sandstone is an important regional bedrock aquifer. The GIS model determines the net-infiltration percentage of precipitation by using an empirical equation. This relation is derived from least squares linear regression between three surficial parameters (soil coarseness, topographic slope, and downgradient distance from outcrop) and the percentage of estimated net infiltration based on environmental tracer data from excavations and boreholes at Sand Hollow Reservoir in the southeastern part of the study area. Processed GIS raster layers are applied as parameters in the empirical equation for determining net infiltration for soil-covered areas as a percentage of precipitation. This net-infiltration percentage is multiplied by average annual Parameter-elevation Regressions on Independent Slopes Model (PRISM) precipitation data to obtain an infiltration rate for each model cell. Additionally, net infiltration on exposed outcrop areas is set to 10 percent of precipitation on the basis of borehole net-infiltration estimates. Soils and outcrop net-infiltration rates are merged to form a final map. Areas of low, medium, and high potential for ground-water recharge have been identified, and estimates of net infiltration range from 0.1 to 66 millimeters per year (mm/yr). Estimated net-infiltration rates of less than 10 mm/yr are considered low, rates of 10 to 50 mm/yr are considered medium, and rates of more than 50 mm/yr are considered high. A comparison of estimated net-infiltration rates (determined from tritium data) to predicted rates (determined from GIS methods) at 12 sites in Sand Hollow and at Anderson Junction indicates an average difference of about 50 percent. Two of the predicted values were lower, five were higher, and five were within the estimated range. While such uncertainty is relatively small compared with the three order-of-magnitude range in predicted net-infiltration rates, the net-infiltration map is best suited for evaluating relative spatial distribution rather than for precise quantification of recharge to the Navajo aquifer at specific locations. An important potential use for this map is land-use zoning for protecting high net-infiltration parts of the aquifer from potential surface contamination.

Utah↗

Magnitude and frequency of floods for rural streams in Georgia, South Carolina, and North Carolina, 2017—Results

Reliable estimates of the magnitude and frequency of floods are an important part of the framework for hydraulic-structure design and flood-plain management in Georgia, South Carolina, and North Carolina. Annual peak flows measured at U.S. Geological Survey streamgages are used to compute flood‑frequency estimates at those streamgages. However, flood‑frequency estimates also are needed at ungaged stream locations. A process known as regionalization was used to develop regression equations to estimate the magnitude and frequency of floods at ungaged locations. A multistate approach was used to update estimates of the magnitude and frequency of floods in rural, ungaged basins in Georgia, South Carolina, and North Carolina. Annual peak-flow data through September 2017 were analyzed for 965 streamgages with 10 or more years of data on rural streams in Georgia, South Carolina, North Carolina, and adjacent parts of Alabama, Florida, Tennessee, and Virginia. Flood‑frequency estimates of the 50‑, 20‑, 10‑, 4‑, 2‑, 1‑, 0.5‑, and 0.2‑percent annual exceedance probability streamflows, which correspond to flood-recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years, respectively, were computed for the 965 streamgages following national guidelines. As part of the computation of flood‑frequency estimates for the streamgages, an updated value for the regional skew coefficient (0.048) was developed using a Bayesian generalized least squares regression model. The new regional skew has a mean square error or average variance of prediction of 0.092. Additionally, basin characteristics for these stations were computed using a geographical information system. Exploratory analyses on the 965 streamgages confirmed the five hydrologic regions for Georgia, South Carolina, and North Carolina defined in a previous rural flood‑frequency study. From the 965 streamgages, streamgages with 30 or more years of record were used to complete a peak-flow trend analysis. Of the 965 streamgages, 164 streamgages were found to be redundant and were excluded from the regional regression analyses. Data from the remaining 801 streamgages (292 in Georgia, 75 in South Carolina, 303 in North Carolina, 15 in Alabama, 12 in Florida, 39 in Tennessee, and 65 in Virginia) were used in a regional regression analysis relating basin characteristics to flood‑frequency estimates. This analysis, based on generalized least squares regression, was used to develop a set of predictive equations to estimate the 50‑, 20‑, 10‑, 4‑, 2‑, 1‑, 0.5‑, and 0.2‑percent annual exceedance probability streamflows for rural, ungaged basins in Georgia, South Carolina, and North Carolina. The final set of predictive equations are all functions of drainage area and percentage of the drainage basin within each of the five hydrologic regions. Average errors of prediction for these regression equations range from 35.8 to 44.4 percent. Flood‑frequency estimates also were computed for 72 regulated (for example, a streamgage where flow is altered by a dam or weir) streamgages in Georgia, South Carolina, and North Carolina with 20 or more years of post-regulation record using data through water year 2019. The water year is the annual period from October 1 through September 30 and is designated by the year in which the period ends. Of the 72 regulated streamgages, 18 had pre-regulated periods of record that also were analyzed as part of this study. Flow adjustments were applied to historic peaks and large floods from the pre-regulated period, if available, for use in the post-regulation frequency analysis. Estimates of large floods provide valuable information in frequency analysis and, thus, were included in the post-regulation frequency analysis.

Georgia, North Carolina, South Carolina↗

Electrical Resistivity and Seismic Surveys at the Nevada Test Site, Nevada, April 2007

In April 2007, the USGS collected direct-current (DC) electrical resistivity data and shear- (S) and compressional- (P) wave seismic data to provide new detail of previously mapped, overlapping fault splays at two administrative areas in the Nevada Test Site (NTS). In NTS Area 7, we collected two-dimensional DC resistivity data along a transect crossing the Yucca Fault parallel to, and between, two transects along which resistivity data were collected in a previous study in 2006. In addition, we collected three-dimensional DC resistivity data in a grid that overlies part of the 2007 transect. The DC resistivity data show that the fault has a footwall that is more conductive than the hanging wall and an along-strike progression of the fault in a location where overlapping splays are present. Co-located with the northernmost of the two 2006 DC resistivity transects, we acquired S- and P-wave seismic data for both reflection and refraction processing. The S-wave data are corrupted by large amounts of converted (P-wave) energy likely due to the abundance of fractured caliche in the shallow subsurface. The P-wave data show minimal reflected energy, but they show clear refracted first arrivals. We have inverted these first arrival times to determine P-wave seismic velocity models. The seismic model for the transect in Area 7 shows low velocities extending to the base of the model at the location of the Yucca Fault, as well as low velocities at the eastern end of the transect, in the vicinity of the adjacent crater. These new surveys provide further detail about the geometry of the Yucca Fault in this location where it shows two overlapping splays. We collected P- and S-wave seismic data along a transect in the southern part of NTS Area 2, corresponding with the location of a 2006 DC resistivity transect that targeted a set of small faults identified with field mapping. Again, the S-wave data are difficult to interpret. The P-wave data show clear first arrivals that we inverted, yielding a velocity model that shows lateral heterogeneity similar to the 2006 DC resistivity models. Finally, we collected P-wave data along a second transect in Area 2, located north of the first line and in an area of a very minor fault that was targeted by another 2006 DC resistivity survey. The P-wave refraction velocity model shows generally high velocities, with a zone of somewhat lower velocities in the central part of the transect. The position of the low velocity zone corresponds with the location of a minor fault, though it is unclear whether the two are related. Together, these results demonstrate the value of geophysical data for mapping the subsurface extent of faults. The 2007 DC resistivity data complement the 2006 data and provide important new detail of the overlapping fault splays. The seismic data demonstrate the ability of P-wave refraction methods to identify the damage zones at faults, and they show the difficulties associated with S-wave methods in areas with caliche. Combining all of the geophysical data from the Area 7 studies, we are able to develop a coherent interpretation of the relation between the site geology, the fault, and the observations.

Open-File Report↗

Notes on the availability of mid-Atlantic multichannel seismic reflection profiles 11 and 12-A, -B, -C, and -D

Available are two multichannel profiles collected by Digicon Geophysical Corporation in 1975 using a 48-channel streamer (3600 m long). and a 27.9 cu. liter air gun array. They were processed in Denver on the Phoenix "I" by William C. Patterson. The processing included demultiplexing and resampling, geometry and common-depth-point definition, velocity analysis, noise muting, band-pass filtering, time-variant filtering, time-variant deconvolution, and automatic gain control (AGC) scaling, prior to the final profile playout. The release includes parts of two lines off the mid-atlantic continental margin (see map) over the Baltimore Canyon trough. These profiles were collected as part of a regional grid over the offshore Atlantic sedimentary basins, as a part of a continuing program to assess the resource potential using non-proprietary data Line 11 is a cross shelf profile (233 ken long) taken across the shelf, slope and upper rise off Virginia (southern Baltimore Canyon trough). Line 12 (shotpoint 1Ol- 6011, parts A, 8, C and D) is an along-the-shelf profiler and stretches from the Cape Hatteras IPOD line to line 2, at the northern part of the Baltimore Canyon trough), a distance of 591 km. These profiles may be viewed at the U.S. Geological Survey, Atlantic Gulf of Mexico, Office of Marine Geology, Quissett campus, Building B, Woods Hole, Massachusetts 02543, and at the U.S. Geological Survey, Branch of Oil and Gas Resources, Denver Federal Center, Building 25, Denver. Colorado 80225. Copies of the profiles can be purchased from the Nations Geophysical and Solar-Terrestrial Data Center, Environmental Data Service- NOAA, Code D621 Boulder, Colorado 80303.

Open-File Report↗

Monitoring and modeling to predict Escherichia coli at Presque Isle Beach 2, City of Erie, Erie County, Pennsylvania

The Lake Erie shoreline in Pennsylvania spans nearly 40 miles and is a valuable recreational resource for Erie County. Nearly 7 miles of the Lake Erie shoreline lies within Presque Isle State Park in Erie, Pa. Concentrations of Escherichia coli (E. coli) bacteria at permitted Presque Isle beaches occasionally exceed the single-sample bathing-water standard, resulting in unsafe swimming conditions and closure of the beaches. E. coli concentrations and other water-quality and environmental data collected at Presque Isle Beach 2 during the 2004 and 2005 recreational seasons were used to develop models using tobit regression analyses to predict E. coli concentrations. All variables statistically related to E. coli concentrations were included in the initial regression analyses, and after several iterations, only those explanatory variables that made the models significantly better at predicting E. coli concentrations were included in the final models. Regression models were developed using data from 2004, 2005, and the combined 2-year dataset. Variables in the 2004 model and the combined 2004-2005 model were log10 turbidity, rain weight, wave height (calculated), and wind direction. Variables in the 2005 model were log10 turbidity and wind direction. Explanatory variables not included in the final models were water temperature, streamflow, wind speed, and current speed; model results indicated these variables did not meet significance criteria at the 95-percent confidence level (probabilities were greater than 0.05). The predicted E. coli concentrations produced by the models were used to develop probabilities that concentrations would exceed the single-sample bathing-water standard for E. coli of 235 colonies per 100 milliliters. Analysis of the exceedence probabilities helped determine a threshold probability for each model, chosen such that the correct number of exceedences and nonexceedences was maximized and the number of false positives and false negatives was minimized. Future samples with computed exceedence probabilities higher than the selected threshold probability, as determined by the model, will likely exceed the E. coli standard and a beach advisory or closing may need to be issued; computed exceedence probabilities lower than the threshold probability will likely indicate the standard will not be exceeded. Additional data collected each year can be used to test and possibly improve the model. This study will aid beach managers in more rapidly determining when waters are not safe for recreational use and, subsequently, when to issue beach advisories or closings.

Pennsylvania↗

Final project memorandum: Identifying conservation objectives for the Gulf Coast habitats of the black skimmer and gull-billed tern

Many shorebirds and nearshore waterbirds are of conservation concern across the Gulf of Mexico due to stressors such as human disturbance, predation, and habitat loss and degradation. Conservation and protection of these birds is important for the functioning of healthy ecosystems and for maintaining biodiversity in North America. Consequently, resource managers along the Gulf need decision-aiding tools that can help to answer important conservation questions for different species (e.g., how much area should be targeted by management actions to meet a species’ needs). To address this need, project researchers developed statistical models that could help identify habitat conservation objectives and actions for bird species taking into account different Gulf coast conservation scenarios that might occur in response to sea-level rise. The project focused specifically on the Black Skimmer ( Rynchops niger ) and Gull-billed Tern ( Gelochelidon nilotica ), two species designated as U.S. Fish and Wildlife Service Species of Conservation Concern and Gulf Coast Joint Venture Priority Species. These two birds are also representative of a variety of other beach and barrier-island nesting birds whose nesting habitats are threatened by sea-level rise (e.g., Least Tern, Snowy and Wilson’s Plover). The statistical models linked each bird’s abundance to habitat characteristics that could be influenced by different management actions. This information could be used to identify conservation objectives under different conservation scenarios.

Gulf of Mexico↗

2023 Coastal master plan: ICM-wetlands – Submerged aquatic vegetation (SAV) updates

Submerged aquatic vegetation (SAV) provides critical structural habitat for valuable nekton and wildlife species across coastal ecosystems and can buffer the negative effects of land loss. Landscape change and restoration efforts across coastal Louisiana can impact the occurrence, coverage, and species assemblages of SAV, and changes to these foundational species can have cascading impacts across food webs. To support the 2023 Coastal Master Plan efforts, a unique SAV model was developed to assess coverage and occurrence of SAV across aquatic waterbodies in response to environmental variables evaluated. This effort created a spatial model describing the probability of presence of SAV across the study area in response to changing conditions over the modeled time period. To develop the initial coverage data layer, we used remotely sensed Normalized Difference Vegetation Index (NDVI) and modified Normalized Difference Water Index (mNDWI) data from 2015-2018 to identify areas containing variable vegetation and water spectral reflectance. Key environmental variables evaluated included total suspended sediments (TSS), salinity, and physical exposure. Seasonal estimates for TSS and salinity were used, as research indicates that seasonal environmental variability is a significant driver for SAV establishment. Seasonal salinity was derived from Coast-wide Reference Monitoring Station (CRMS) data, and seasonal TSS was estimated from hyperspectral imagery. Estimates of physical exposure have previously been provided by calculating fetch (the distance across water over which waves can propagate), but this proved to be too computationally intensive to be feasible, and we found distance to land to be a reasonable proxy for exposure. To represent geographic conditions and historical factors influences on SAV establishment and occurrence (e.g., variables too numerous and complex to model) we developed a basin variable that served as a proxy for complex historical, or prior, conditions, determined by the forested, fresh, intermediate, brackish, or saline (FFIBS) score. The final model included spring TSS, spring salinity, distance to land, and the basin prior. The model performed well for the area evaluated, correctly classifying SAV (as present or absent) 89% of the time (Kappa = 580). SAV probability of presence responded as expected to change in these environmental variables, with likelihood of occurrence decreasing in response to increasing spring TSS, spring salinity, and distance to land. However, the model was more accurate at predicting absence (true negative = 0.940) than predicting presence (true positive = 0.626), suggesting that the scale of the model may limit the ability to predict presence. Moreover, the simplicity of the model limited the accuracy in highly dynamic environments, for example near the outflow of diversions or areas of significant changes in salinity or TSS. Through incorporating underwater communities like SAV, this master plan provides a holistic view of coastal change and restoration. To create healthy ecological structure and function in wetland habitats, the submergent communities must be considered alongside the emergent habitats. As the benefits of SAV are increasingly recognized, both here in Louisiana and beyond, SAV restoration and the use of SAV communities in assessing and improving ecological condition are becoming more common.

Louisiana↗

Detrital zircon geochronology of pre- and syncollisional strata, Acadian orogen, Maine Appalachians

The Central Maine Basin is the largest expanse of deep-marine, Upper Ordovician to Devonian metasedimentary rocks in the New England Appalachians, and is a key to the tectonics of the Acadian Orogeny. Detrital zircon ages are reported from two groups of strata: (1) the Quimby, Rangeley, Perry Mountain and Smalls Falls Formations, which were derived from inboard, northwesterly sources and are supposedly older; and (2) the Madrid, Carrabassett and Littleton Formations, which were derived from outboard, easterly sources and are supposedly younger. Deep-water deposition prevailed throughout, with the provenance shift inferred to mark the onset of foredeep deposition and orogeny. The detrital zircon age distribution of a composite of the inboard-derived units shows maxima at 988 and 429 Ma; a composite from the outboard-derived units shows maxima at 1324, 1141, 957, 628, and 437 Ma. The inboard-derived units have a greater proportion of zircons between 450 and 400 Ma. Three samples from the inboard-derived group have youngest age maxima that are significantly younger than the nominal depositional ages. The outboard-derived group does not share this problem. These results are consistent with the hypothesised provenance shift, but they signal potential problems with the established stratigraphy, structure, and (or) regional mapping. Shallow-marine deposits of the Silurian to Devonian Ripogenus Formation, from northwest of the Central Maine Basin, yielded detrital zircons featuring a single age maximum at 441 Ma. These zircons were likely derived from a nearby magmatic arc now concealed by younger strata. Detrital zircons from the Tarratine Formation, part of the Acadian foreland-basin succession in this strike belt, shows age maxima at 1615, 980 and 429 Ma. These results are consistent with three episodes of zircon recycling beginning with the deposition of inboard-derived strata of the Central Maine Basin, which were shed from post-Taconic highlands located to the northwest. Next, southeasterly parts of this succession were deformed in the Acadian orogeny, shedding detritus towards the northwest into what remained of the basin. Finally, by Pragian time, all strata in the Central Maine Basin had been deformed and detritus from this new source accumulated as the Tarratine Formation in a new incarnation of the foreland basin. Silurian-Devonian strata from the Central Maine Basin have similar detrital zircon age distributions to coeval rocks from the Arctic Alaska and Farewell terranes of Alaska and the Northwestern terrane of Svalbard. We suggest that these strata were derived from different segments of the 6500-km-long Appalachian-Caledonide orogen.

Maine, New Hampshire↗

The hydrogeology of the Tully Valley, Onondaga County, New York: an overview of research, 1992-2012

Onondaga Creek begins approximately 15 miles south of Syracuse, New York, and flows north through the Onondaga Indian Nation, then through Syracuse, and finally into Onondaga Lake in central New York. Tully Valley is in the upper part of the Onondaga Creek watershed between U.S. Route 20 and the Valley Heads end moraine near Tully, N.Y. Tully Valley has a history of several unusual hydrogeologic phenomena that affected past land use and the water quality of Onondaga Creek; the phenomena are still present and continue to affect the area today (2014). These phenomena include mud volcanoes or mudboils, landslides, and land-surface subsidence; all are considered to be naturally occurring but may also have been influenced by human activity. The U.S. Geological Survey (USGS), in cooperation with the U.S. Environmental Protection Agency and the Onondaga Lake Partnership, began a study of the Tully Valley mudboils beginning in October 1991 in hopes of understanding (1) what drives mudboil activity in order to remediate mudboil influence on the water quality of Onondaga Creek, and (2) land-surface subsidence issues that have caused a road bridge to collapse, a major pipeline to be rerouted, and threatened nearby homes. Two years into this study, the 1993 Tully Valley landslide occurred just over 1 mile northwest of the mudboils. This earth slump-mud flow was the largest landslide in New York in more than 70 years (Fickies, 1993); this event provided additional insight into the geology and hydrology of the valley. As the study of the Tully Valley mudboils progressed, other unusual hydrogeologic phenomena were found within the Tully Valley and provided the opportunity to perform short-term, small-scale studies, some of which became graduate student theses—Burgmeier (1998), Curran (1999), Morales-Muniz (2000), Baldauf (2003), Epp (2005), Hackett, (2007), Tamulonis (2010), and Sinclair (2013). The unusual geology and hydrology of the Tully Valley, having been investigated for more than two decades, provides the basis for this report.

New York↗

Dacitic ash-flow sheet near Superior and Globe, Arizona

Remnants of a dacitic ash-flow sheet near Globe, Miama, and Superia, Arizona cover about 100 square miles; before erosion the area covered by the sheet was at least 400 square miles and perhaps as much as 1,500 square miles. Its maximum thickness is about 2,000 feet, its average thickness is about 500 feet, and its original volume was at least 40 cubic miles. It was erupted on an eroded surface with considerable relief. The main part of the deposit was thought by early workers to be a lava flow. Even after the distinctive character of welded tuffs and related rocks was discovered, the nature and origin of this deposit remained dubious because textures did not correspond to those in other welded tuff bodies. Yet a lava flow as silicic as this dacite would be viscous instead of spreading out as an extensive sheet. The purpose of this investigation has been to study the deposit, resolve the inconsistencies, and deduce its origin and history. Five stratigraphic zones are distinguished according to differences in the groundmass. From bottom to top the zones are basal tuff, vitrophyre, brown zone, gray zone, and white zone. The three upper zones are distinguished by colors on fresh surfaces, for each weathers to a similar shade of light reddish brown. Nonwelded basal tuff grades upward into the vitrophyre, which is a highly welded tuff. The brown and gray zones consist of highly welded tuff with a lithoidal groundmass. Degree of welding decreases progressively upward through the gray and the white zones, and the upper white zone is nonwelded. Textures are clearly outlined in the lower part of the brown zone, but upward they become more diffuse because of increasing devitrification. In the white zone, original textures are essentially obliterated, and the groundmass consists of spherulites and microcrystalline intergrowths. The chief groundmass minerals are cristobalite and sanidine, with lesser quartz and plagioclase. Phenocrysts comprise about 40 percent of the rock, and their relative proportions are fairly uniform. Almost three-fourths of the phenocrysts are plagioclase, one-tenth quartz, one-tenth biotite, and the remainder sanidine, magnetite, and hornblende, with accessory sphene, zircon, and appetite. Pumice fragments are nearly equidimensional near the top of the sheet, and downward they become progressively more flattened until they finally disappear. The zones and the pumice fragment flattening ration (ratio of length to height) provide means for recognizing several faults within the sheet. Twelve new chemical analyses are nearly uniform in composition. If named according to chemical composition, the rock would be a quartz latite, but when named according to phenocrysts, it is a dacite. From the field occurrence and the interpretation of relict textures, it is concluded that the deposit is an ash-flow sheet containing large amounts of welded tuff, and that it was emplaced by a type of nuee ardente instead of a lava flow or air-fall shower. The nature of zoning and trend of flattening ratios indicate a series of eruptions in rapid enough succession for the sheet to form a single cooling unit. Except in the lower part of the sheet, original textures were obscured by devitrification and crystallization during cooling. Nearly uniform mineralogy and chemistry suggest a single magnetic source. A nearly circular area, about 3? miles in diameter, of altered dacite and earlier volcanic rocks, bounded by intricately faulted and brecciated older rocks, may be the site of a caldera that represents the source of the eruptions.

Open-File Report↗

Hydrogeologic investigation, framework, and conceptual flow model of the Antlers aquifer, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statute §82–1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. Every 20 years, the Oklahoma Water Resources Board is required to update the hydrologic investigation on which the maximum annual yield determinations were based. The maximum annual yield allocated per acre of land is used to set the equal-proportionate share pumping rate. The maximum annual yield of 5,913,600 acre-feet per year and equal-proportionate-share of 2.1 acre-feet per acre per year currently (2025) in place for the Antlers aquifer were issued by the Oklahoma Water Resources Board on February 14, 1995. Because more than 20 years have elapsed since the 1995 final order for the Antlers aquifer was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, completed an in-depth hydrologic study that included a hydrogeologic framework and conceptual groundwater-flow model for the 1980–2022 study period. The results of an analysis of land use, long-term climate patterns, streamflow and base-flow patterns, historical groundwater use, as well as groundwater-level fluctuations across the Antlers aquifer are described. In addition, groundwater quality was analyzed for total dissolved solids concentrations and major ions for the Antlers aquifer. An updated hydrogeologic framework was developed that included refining the aquifer boundary in Oklahoma, the creation of new potentiometric surface and saturated thickness of fresh groundwater maps, one multiple-well aquifer test, slug tests, and an analysis of lithologic logs across the aquifer. A conceptual groundwater flow model and water budget were developed by incorporating estimates of recharge from precipitation, saturated-zone evapotranspiration, streambed seepage, lateral groundwater flows, vertical leakage, and withdrawals from groundwater wells.

Oklahoma, Texas↗

Cretaceous sedimentation and tectonism in the southeastern Kaiparowits region, Utah

Upper Cretaceous strata in the southeastern Kaiparowits region of south-central Utah consist of approximately 3,500 feet of interfingering sandstone, mudstone, shale, and coal in the Dakota Formation (oldest), Tropic Shale, Straight Cliffs Formation, and Wahweap Formation (youngest). The formations consist of several depositional facies that can be recognized by characteristic lithologies bedding structures, and fossils; these are the alluvial plain, deltaic plain, lagoonal-paludal, barrier sandstone, and offshore marine facies. The distribution of facies clearly defines the paleogeography of the region during several cycles of marine transgression and regression. The nonmarine beds were deposited on a broad alluvial coastal plain that was bordered on the west and southwest by highlands and on the east and northeast by the Western Interior seaway. The marine beds were deposited whenever the seaway advanced into or across the region. The Dakota Formation and the lower part of the Tropic Shale were deposited in nonmarine and marine environments, while the shoreline advanced generally westward across the region. The middle and upper part of the Tropic Shale and the Tibbet Canyon and Smoky Hollow Members of the Straight Cliffs Formation were deposited in marine and nonmarine environments when the seaway had reached its greatest areal extent and began a gradual northeastward withdrawal. An unconformity at the top of the Smoky Hollow represents a period of erosion and possibly nondeposition before deposition of the John Henry Member of the Straight Cliffs. The John Henry Member grades from nonmarine in the southwest to predominantly marine in the northeast, and was deposited during two relatively minor cycles of transgression and regression. The Drip Tank Member at the top of the Straight Cliffs Formation is a widespread sandstone unit deposited mainly in fluvial environments. Some of the beds in the northeastern part of the region were probably deposited in marine waters during the final incursion of the seaway into the Kaiparowits region. The overlying Wahweap Formation was deposited in nonmarine environments. Slight but continued tectonism during Late Cretaceous time is indicated by lateral changes of facies and thickness variations that coincide at least partly with present structures. These criteria indicate that Laramide tectonism consisted of two phases. An early phase that lasted from about late Albian to late Campanian time included regional subsidence, basin downwarping, and movement on local folds and faults. A later phase that lasted from late Campanian to about late Paleocene time included regional uplift, monoclinal flexing, and probable new faulting, as well as continued basin downwarping and movement on local folds and probably on the older faults. The principal economic resource in the Kaiparowits region is bituminous or subbituminous coal in the john Henry Member. Because basin downwarping and movement on local folds occurred during deposition, the thicker and more continuous coal beds are in the ancestral synclines and tile deeper part of the structural basin. Presently indicated resources total 7.3 billion tons, but considerably larger quantities are probably present in the unexplored parts of the region. Several potential resources include ground water, titaniferous sandstone, and possibly oil and gas.

Utah↗

The quartz crystal deposits of western Arkansas

Extensive geological investigations of domestic deposits of quartz crystals were undertaken by the Geological Survey of the United States Department of the Interior in the fall of 1942 in response to the urgent demand for crystals for war purposes. The work was continued throughout 1943 and took survey geologists into some 20 states in which over 1,400 deposits or reports of deposits were examined. Of the numerous highly varied occurrences of quartz examined, only the deposits in western Arkansas, in the Piedmont and Blue Ridge provinces of Virginia, North Carolina and Georgia and the placer quartz of Mokelumne Hill, California, were of sufficient promise to warrant detailed work. Production in 1943 from these areas totaled slightly less than 4 tons of oscillator quartz, over 85 per cent of which came from Arkansas.The Arkansas deposits occur throughout the thick, deformed Paleozoic shales, sandstones, and cherts exposed along the central belt of the Ouachita Mountains. Steeply dipping fractures closely related to the major folds control the deposition of most quartz.Clear quartz is confined largely to the terminal parts of primary crystals, which have developed without disturbance or interference, are commonly elongate parallel to their C axes and are bounded by relatively simple forms. Deposition of silica during and subsequent to the fracturing of crystals resulted in the formation of complex crystals which are characterized by extensive optical twinning and lineage structures, and commonly are bounded by aggregates of the simpler crystal forms. The principal defects in all types of crystals are twinning, smokiness, cavities, solid inclusions and fractures.The Arkansas quartz deposits include veins, sheeted zones and stock-works. They are largely cavity fillings, apparently deposited by rising, attenuated, hydrothermal solutions, at relatively low temperatures and pressures. Minerals associated with the quartz, which constitutes 90 per cent or more of the cavity fillings, include dickite and carbonaceous material, calcite, adularia, and chlorite. The constituents of the cavity fillings could have been derived principally from magmatic sources or more probably from underlying rocks, with small additions from the rocks enclosing the cavities. The complex vein fabrics apparently resulting from intermittent regional deformation during the deposition of the quartz, and certain structural relations of the deposits, indicate that they were formed in the final stages of the Ouachita orogeny, probably in mid-Pennsylvanian time.

Arkansas↗

Methods for determining magnitude and frequency of floods in California, based on data through water year 2006

Methods for estimating the magnitude and frequency of floods in California that are not substantially affected by regulation or diversions have been updated. Annual peak-flow data through water year 2006 were analyzed for 771 streamflow-gaging stations (streamgages) in California having 10 or more years of data. Flood-frequency estimates were computed for the streamgages by using the expected moments algorithm to fit a Pearson Type III distribution to logarithms of annual peak flows for each streamgage. Low-outlier and historic information were incorporated into the flood-frequency analysis, and a generalized Grubbs-Beck test was used to detect multiple potentially influential low outliers. Special methods for fitting the distribution were developed for streamgages in the desert region in southeastern California. Additionally, basin characteristics for the streamgages were computed by using a geographical information system. Regional regression analysis, using generalized least squares regression, was used to develop a set of equations for estimating flows with 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities for ungaged basins in California that are outside of the southeastern desert region. Flood-frequency estimates and basin characteristics for 630 streamgages were combined to form the final database used in the regional regression analysis. Five hydrologic regions were developed for the area of California outside of the desert region. The final regional regression equations are functions of drainage area and mean annual precipitation for four of the five regions. In one region, the Sierra Nevada region, the final equations are functions of drainage area, mean basin elevation, and mean annual precipitation. Average standard errors of prediction for the regression equations in all five regions range from 42.7 to 161.9 percent. For the desert region of California, an analysis of 33 streamgages was used to develop regional estimates of all three parameters (mean, standard deviation, and skew) of the log-Pearson Type III distribution. The regional estimates were then used to develop a set of equations for estimating flows with 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities for ungaged basins. The final regional regression equations are functions of drainage area. Average standard errors of prediction for these regression equations range from 214.2 to 856.2 percent. Annual peak-flow data through water year 2006 were analyzed for eight streamgages in California having 10 or more years of data considered to be affected by urbanization. Flood-frequency estimates were computed for the urban streamgages by fitting a Pearson Type III distribution to logarithms of annual peak flows for each streamgage. Regression analysis could not be used to develop flood-frequency estimation equations for urban streams because of the limited number of sites. Flood-frequency estimates for the eight urban sites were graphically compared to flood-frequency estimates for 630 non-urban sites. The regression equations developed from this study will be incorporated into the U.S. Geological Survey (USGS) StreamStats program. The StreamStats program is a Web-based application that provides streamflow statistics and basin characteristics for USGS streamgages and ungaged sites of interest. StreamStats can also compute basin characteristics and provide estimates of streamflow statistics for ungaged sites when users select the location of a site along any stream in California.

California↗

Review of literature on the endangered masked bobwhite

The masked bobwhite (Colinus virginianus ridgwayi) once inhabited restricted areas in southern Arizona and middle Sonora, Mexico. It probably never was a widespread and abundant bird. Ornithologists discovered this race in 1884, presumably during its final decline in Arizona due to overgrazing and a series of droughts. . It was gone from Arizona by 1900. Apparently the bird was not seriously threatened in Sonora until the 1940's when the cattle industry increased there. Only small and scattered populations remain in Sonora today. Although behavior of the masked bobwhite is similar to that of the closely-related eastern bobwhite (C. v. virginianus), the desert variety prefers a mesquite-grassland habitat at elevations of from 1,000 to 4,000 feet, has a later and shorter nesting season, and has more striking sexual plumage. The male's black head and 'robin's red' breast readily identify the bird. Average annual precipitation in the heart of its Sonoran habitat is 13.5 inches, with 75 percent of that occurring during the 3-month period of July through September. Foods consist of small weed and grass seeds, supplemented with invertebrates and green vegetational growth during the summer rainy season. Early attempts at reintroduction of masked bobwhites into Arizona failed. Experiments on propagation and reestablishment and a life history study are currently being carried out by the Bureau of Sport Fisheries and Wildlife in cooperation with the Arizona Fish and Game Department.

Resource Publication↗