Geology ReportsSearch

SEARCH · Geology Reports

Results for “Journal of Analytical and Applied Pyrolysis”

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.

1,689 records · Page 20Linked to original sources

Methods to evaluate and improve the modeling of rupture directivity in assessment of seismic hazard

In recent years, there have been several advancements related to the modelling of near-source effects of earthquake rupture on strong ground shaking, leading to an improved characterization of ground motions and resulting seismic hazard. Some of these modifications have stemmed from physics-based numerical modelling of the earthquake rupture process, using physics-based dynamic rupture simulations. These contributions have led to a better understanding of how fault rupture characteristics, geometry, and the style of faulting can interact with the hypocenter-dependence on the path from source to site that may ultimately guide the development of seismic directivity models. Moving forward, the application of modern techniques can be used to incorporate these source characteristics and near-fault ground motion behavior that contribute to the azimuthally varying effects that result in rupture directivity. One example is the application of machine learning methods to support more automated integration of new predictor variables in model development and open more evaluation opportunities to access residuals. Here, we utilize several techniques to take advantage of the plethora of synthetic data and its ability to supplement preexisting trends observed in data. We showcase two examples of how models can be either developed, expanded upon, or constrained using artificial neural network model (ANNs). We evaluate the performance of the ANN with existing methods, comparing misfit, potential limitations, and ability to continue to improve upon these methods in the future. One approach uses a set of simulations with corresponding synthetic ground motions from the Southern California Earthquake Center (SCEC) CyberShake study to develop a ground motion model adapted to incorporate seismic directivity information using an ANN. This large database (TBs) enables us to train the model to capture magnitude, period, and distance variations and how these parameters relate to amplification from hypocenters located along finite-faults. In some cases, there is reduced misfit from better representing source features that aren’t included in base ground motion models that neglect hypocenter location (e.g. azimuthal variation, source-to-site terms). Another ANN method uses a shallow-layered neural network model to better fit a hypocenter-independent model. This method adjusts the median and aleatory variability to account for the averaged impact of various hypocenter distributions to fit the underlying directivity adjustment model. This method serves as a template to apply to other directivity models, improving computational efficiency and more readily enabling integration in hazard codes.

California

Ground deformation and gravity for volcano monitoring

Introduction When magma accumulates or migrates, it can cause pressurization and related ground deformation. Characterization of surface deformation provides important constraints on the potential for future volcanic activity, especially in combination with seismic activity, gas emissions, and other indicators. A wide variety of techniques and instrument types have been applied to the study of ground deformation at volcanoes (sidebar, p. 2; Dzurisin, 2000, 2003, 2007). Geodetic instruments include continuously recording Global Navigation Satellite System (GNSS; of which the United States’ Global Positioning System is one example) stations (fig. D1), borehole tiltmeters, and interferometric synthetic aperture radar (InSAR) measurements (from satellites, occupied and unoccupied aircraft systems, and ground-based sensors). Additional geodetic measurements like continuous- and survey-mode gravity (fig. D2) can contribute substantially to interpreting these data. Borehole strainmeters (see chapter K , this volume, by Hurwitz and Lowenstern, 2024) also have outstanding utility for monitoring deformation, although because of cost and permitting challenges, we do not include them as part of standard volcano monitoring networks for U.S. volcanoes. Still other techniques like light detection and ranging (lidar), structure from motion, and optical satellite data can be used to derive gross topographic changes, which can be used to map volcanic deposits, infer eruption rates, and gain insights into the source processes associated with eruptive activity (see chapter G , this volume, on tracking surface changes caused by volcanic activity; Orr and others, 2024). Experience has shown that no single geodetic monitoring technique is adequate to detect and track the entire range of ground-motion patterns that occur at volcanoes, primarily because of the temporal and spatial diversity of volcano deformation (fig. D3). Similarly, the magnitude of surface deformation varies widely. Geodetic monitoring strategies should therefore include multiple techniques and instrument types to cover a wide range of spatial and temporal scales. In identifying recommendations for geodetic instrumentation for volcano monitoring networks, we attempted to maximize the diversity of instrument types to measure the full range of deformation signals and minimize their expense and number; thus, we do not include several well-known deformation-monitoring techniques in our recommendations. Extensometers, for example, measure strains over distances of a few meters and have an excellent record of success in detecting changes in preeruptive localized ground motion across existing cracks, including at Mount St. Helens, Washington (Iwatsubo and others, 1992), and Piton de la Fournaise, Réunion Island (Peltier and others, 2006). Despite being relatively inexpensive, extensometers are best used primarily when localized ground displacements (for example, ground cracks) need to be tracked, and are not necessary at all volcanoes. In considering volcano deformation monitoring strategies, two complicating factors are deserving of special attention. First, not all deformation is driven by subsurface magmatic activity—for example, at many large stratovolcanoes (for example, Mount Rainier), flank collapses and landslides are significant geologic hazards (Reid and others, 2001) that may occur even in the absence of magmatic activity. Monitoring the stability of volcanoes is thus another critical application of geodetic monitoring networks to inform hazard assessment. One of the most famous examples of edifice instability is the large flank collapse that initiated the May 18, 1980, eruption of Mount St. Helens. Deformation monitoring had detected a bulge on the north flank of the mountain in April 1980 that was expanding by several meters per day (Lipman and others, 1981). Given that flank collapses can happen at any time during a period of volcanic unrest (or even outside a period of unrest), the capability to assess edifice stability is critical. Second, although volcanoes are commonly treated as idealized structures that erupt from single points, like centralvent stratovolcanoes, many are characterized by long rift zones from which eruptions may originate, and distributed volcanic fields are characterized by broadly spaced vents. For example, linear dikes are common at Kīlauea, Mauna Loa, and between Mount Shasta and Medicine Lake in California. At Kīlauea, one of these linear dikes emerged more than 40 kilometers (km) away from the summit of the volcano during the lower East Rift Zone eruption in 2018. Other volcanic fields, like Lassen volcanic center, California, or the San Francisco Volcanic Field, Arizona, have many small vents spread over a wide area. Although the instrumentation guidelines presented in this chapter remain phrased for central-vent volcanoes, they should be modified as needed in the context of the eruptive characteristics of each individual volcanic system. Spatial analysis of geodetic network coverage could help to ensure adequate instrumentation in areas where volcanism can occur over a broad area as opposed to a central vent. As an example, consider the adjacent volcanoes Mount Shasta and Medicine Lake. If station locations are chosen based only on the distance from the centers of the volcanoes, then any geodetic anomalies between the two volcanoes—an area of potential volcanism as indicated by the presence of volcanic features—may remain undetected by ground-based instrumentation. The spatial analysis is accomplished via a grid of pressure point sources (Mogi, 1958) evenly distributed across the map area, at a depth of 5 km in this example (fig. D4). Each source is inflated until predicted deformations exceed the GNSS white noise uncertainty estimates at one site (Langbein, 2017; Murray and Svarc, 2017). This volume of detectable magma provides a measure of the quality of the coverage (fig. D4). The results indicate that, as of 2022, there is a large area between Mount Shasta and Medicine Lake volcano with existing mapped dikes in which a substantial amount of magma could intrude without being detected geodetically. Applying this style of analysis to individual volcanic systems can provide a guide for designing network geometry given the expected locations of future eruptions.

Scientific Investigations Report

Hidden system identification: Basin modeling as a tool for examining sedimentary geothermal resource potential

Three-dimensional (3D) geologic and temperature models have been developed for the onshore U.S. Gulf Coast. The results from these models identify areas of moderate- to high-temperature (90°-150°C and >150°C; respectively) geothermal resources at depths <6 km. This modeling study addresses the fundamental challenge of predicting where opportune temperature and lithology coincide. Unlike traditional geothermal systems with surface expressions of hydrothermal circulation (e.g., hot springs, fumaroles, sinter), sedimentary geothermal systems (SGS) are generally hidden. Historically, simplified efforts to predict subsurface temperatures in sedimentary basins have focused on linear temperature extrapolation that does not consider the variable thermal properties of different lithologies or lithologic changes with depth (e.g., compaction, lithification). Therefore, the need to understand basin architecture and predict temperatures in 3D within SGS is paramount to identifying geothermal resources and determining economic feasibility. Basin modeling software has long been used to characterize the subsurface conditions of sedimentary basins, including temperature, in the pursuit of finding hydrocarbons. This tool can also be adapted to evaluate the potential of geothermal resources in a sedimentary basin by predicting the confluence of desirable temperatures and reservoir lithologies. In this work, PetroMod basin modeling software was used to create a regional geologic model of the onshore U.S. Gulf Coast, covering over 500,000 km 2 calibrated to temperature data from wells. Inputs include structural surfaces from commercial databases, lithology information derived from published literature, and corrected bottom-hole temperatures (BHT) from over 6,000 wells. The resulting 3D geologic model can be used to predict temperatures throughout the basin. Maps were exported showing the depth, depositional unit, and reservoir lithology at which temperatures of 90°C and 150°C were reached, revealing over 400,000 km 2 of moderate- to high-temperature resources at depths <6 km. These maps function as a first-order screening tool to identify areas where low-, moderate-, or high-grade resource potential may exist, based on temperature and if optimal reservoir lithologies or depositional units of interest are present. Depending on the success criteria of a project, the same maps can be exported for any isotherm or incorporate other 1407 Gardner and Birdwell subsurface properties. The methodology employed in this work can be applied in any sedimentary basin with available subsurface data. Further calibration incorporating other data, including pressure and porosity, can expand the utility of basin modeling for geothermal evaluations. Basin modeling is a powerful but underutilized tool for identifying prospective geothermal resources in sedimentary basins.

Conference Paper

The value of electronic tagging and tracking studies for understanding fish–hypoxia interactions

Objective Hypoxia (i.e., low dissolved oxygen availability) is a natural phenomenon but can also be induced by human activities (e.g., nutrient enrichment from runoff). Given that dissolved oxygen is essential for aquatic life, periods of hypoxia tend to have negative consequences (e.g., sublethal disturbances, mortality) for most fishes. Extensive laboratory research has documented hypoxia thresholds and physiological and behavioral outcomes for a variety of freshwater and marine fishes and there is also an extensive body of fieldwork assessing population-level responses (e.g., survival, distribution). Yet, studying how individual fish respond to hypoxia in the wild has proved challenging; electronic tagging and tracking tools (e.g., biotelemetry, biologging) have made it easier to study individual responses to hypoxia and complement other tools like hydroacoustics that tend to focus on population-level responses. Methods We review what has been learned from contemporary studies that employ electronic tagging and tracking tools to understand how fish respond to hypoxia in the wild. Topics explored include identifying and validating hypoxia thresholds, habitat compression, connectivity, mortality, physiological and bioenergetic consequences, and extreme weather conditions. We also consider what we have learned about evaluating various management plans for hypoxia and reflect on how electronic tags have also been used in aquaculture systems and in hybrid studies that combine laboratory and field research. We highlight fishes in Lake Erie as a research narrative to demonstrate how electronic tagging and tracking have markedly improved our understanding of a longstanding hypoxia issue. Results Our synthesis revealed that electronic tagging and tracking have provided critical information on how fish respond to hypoxia in the field, revealing complex trade-offs and compensatory mechanisms as well as cryptic hypoxia-induced mortality. Beyond just illuminating space use and mortality, tags equipped with various sensors are revealing how fish deal with hypoxia in real time in terms of physiology, bioenergetics, and behavior. Conclusions As electronic tagging and tracking methods experience further innovation and are increasingly applied to understand the effects of hypoxia on fish, we expect more unanticipated findings about the effects of hypoxia on fish in all aquatic ecosystems, which will strengthen our ability to manage and mitigate hypoxia. Combining tools and approaches (e.g., lab and field) is perhaps the best way to generate comprehensive understanding.

Transactions of the American Fisheries Society

Evaluation of best management practices at an edge-of-field site in the Eagle Creek watershed, Ohio, 2012–20

In 2010, the U.S. Geological Survey worked in partnership with the Great Lakes Restoration Initiative and the Natural Resources Conservation Service to identify farm fields in priority watersheds—watersheds critically important to Great Lake health. Intensive best management practices were implemented in five States to test the efficacy of best management practices, also referred to as agricultural conservation practices, to reduce sediment and nutrient runoff on agricultural fields. These fields were chosen as representatives of the priority watersheds because their farming practices and geographical conditions were common among farms in those watersheds. A farm site located near Findlay, Ohio, in the Eagle Creek watershed used two conservation practices, cover crop (Natural Resources Conservation Service practice standard 340) and variable rate technology (VRT; Natural Resources Conservation Service practice standard 590), which represent common practices in the Maumee priority watershed. This study monitored surface runoff and subsurface tile runoff at the Eagle Creek watershed site. The effects of cover crop and VRT nutrient application on sediment and nutrient runoff were assessed from October 2012 to September 2020. Cover crops were applied in fall of 2016, 2017, and 2018, coinciding with VRT applications. The following parameters were analyzed as part of the study: runoff discharge; peak discharge; and concentrations and loads of total phosphorus, particulate phosphorus, dissolved reactive phosphorus, and suspended sediment. High peak runoff events disproportionally affected nutrient concentrations and loads in runoff. A threshold of the 85th percentile of the peak discharge was chosen to assess the effect of cover crop and VRT nutrient application on nutrient mitigation at each gage. Runoff events below the 85th percentile were considered small, and runoff events over the 85th percentile were considered large. Concentrations from collected water-quality samples were used for the analysis; estimated concentrations were excluded from the statistical evaluation. All surface runoff parameters during small runoff events, except surface runoff volume, were significantly lower in the period with cover crops and VRT implementation. Contrarily, subsurface runoff parameters during small events were not statistically different between the periods with and without cover crops and VRT implementation. Surface and subsurface parameters during large events were also not statistically different between the periods with and without cover crops and VRT implementation. These results indicate cover crops and VRT may improve water quality during small runoff events, but additional best management practices that mitigate large runoff events may lead to greater water-quality improvement given the contribution of large events to overall losses.

Ohio

Constraining landslide frequency across the United States to inform county-level risk reduction

Informative landslide hazard estimates are needed to support landslide mitigation strategies to reduce landslide risk across the United States. Whereas existing national-scale landslide susceptibility products assess where landslides are likely to occur, they do not address how often , which is a critical element of landslide hazard and risk assessments. In particular, the U.S. Federal Emergency Management Agency's National Risk Index (NRI) requires landslide frequency estimates to inform expected annual loss estimates. We present county-level landslide frequency (landslides per area per year) estimates for the 50 US states. We applied Bayesian negative binomial regression to estimate both the expected (average) reported landslide frequency and full distribution of annual landslide counts for each county. We compared a suite of models that used combinations of landslide-susceptible area, probability of potentially triggering earthquakes, frequency of potentially triggering precipitation, and ecological region as predictors. We trained our models with landslide inventory data from counties with the most comprehensive records available nationwide and used zero-inflated negative binomial distributions as an incompleteness model to correct for temporal reporting gaps. We selected a preferred frequency model to inform the NRI based on information criteria and physically plausible parameter estimates. The model showed that average annual reported landslide frequencies vary by 5 orders of magnitude across US counties, ranging from 0.002 (0.00015–0.05) landslides 1000 km −2 yr −1 in Kusilvak Census Area, Alaska, to 29 (19–46) landslides 1000 km −2 yr −1 in Lake County, California, reflecting the country's strong variations in landslide susceptibility, earthquake probability, and other factors for which ecological region serves as a proxy. Counties with estimated frequencies in the top 20 % of all counties are predominately along the West Coast of the continental United States, in mountainous regions of the Pacific Northwest and Intermountain West, in locally steep or earthquake-prone regions of the Midwest and Southeast, along the Appalachians, in southern and southeastern Alaska, and on some Hawaiian islands. By examining the number of landslides predicted in 99th percentile years for each county, we identified that 26 % of US counties likely have potential for widespread landsliding with more than 10 landslides 1000 km −2 yr −1 , even when such large events have not been reported in the training data for that county. Overall, our results better represent the range of possible landslide frequencies and spatial variations than previous national-scale estimates reported in the NRI, and our approach can inform other risk-reduction and loss-mitigation efforts across the United States and globally.

Natural Hazards and Earth System Sciences

The 2023 US 50-State National Seismic Hazard Model: Overview and implications

The US National Seismic Hazard Model (NSHM) was updated in 2023 for all 50 states using new science on seismicity, fault ruptures, ground motions, and probabilistic techniques to produce a standard of practice for public policy and other engineering applications (defined for return periods greater than ∼475 or less than ∼10,000 years). Changes in 2023 time-independent seismic hazard (both increases and decreases compared to previous NSHMs) are substantial because the new model considers more data and updated earthquake rupture forecasts and ground-motion components. In developing the 2023 model, we tried to apply best available or applicable science based on advice of co-authors, more than 50 reviewers, and hundreds of hazard scientists and end-users, who attended public workshops and provided technical inputs. The hazard assessment incorporates new catalogs, declustering algorithms, gridded seismicity models, magnitude-scaling equations, fault-based structural and deformation models, multi-fault earthquake rupture forecast models, semi-empirical and simulation-based ground-motion models, and site amplification models conditioned on shear-wave velocities of the upper 30 m of soil and deeper sedimentary basin structures. Seismic hazard calculations yield hazard curves at hundreds of thousands of sites, ground-motion maps, uniform-hazard response spectra, and disaggregations developed for pseudo-spectral accelerations at 21 oscillator periods and two peak parameters, Modified Mercalli Intensity, and 8 site classes required by building codes and other public policy applications. Tests show the new model is consistent with past ShakeMap intensity observations. Sensitivity and uncertainty assessments ensure resulting ground motions are compatible with known hazard information and highlight the range and causes of variability in ground motions. We produce several impact products including building seismic design criteria, intensity maps, planning scenarios, and engineering risk assessments showing the potential physical and social impacts. These applications provide a basis for assessing, planning, and mitigating the effects of future earthquakes.

Earthquake Spectra

Ground-water quality near an inactive landfill and sludge-spreading area, Tallahassee, Florida

Ground-water quality of a 120-acre area encompassing an inactive landfill southwest of Tallahassee, Florida, was assessed by installing and sampling 21 monitoring wells. Sludge from a municipal sewage-treatment plant was also applied to the land surface near the landfill site. Water samples were collected from the 21 monitoring wells from June 1987 through September 1990 and analyzed for concentrations of major ions, nitrogen species, trace metals, field measurements, and more than 140 organic constituents. Water quality in the Upper Floridan aquifer at the study site was compared to background water quality in other areas of Leon County. Ground-water quality of the study site was evaluated by the location of wells sampled in relation to the landfill and the sludge-spreading areas. Results of nonparametric statistical analysis indicated that potassium and nitrite plus nitrate as nitrogen concentrations in water samples from wells completed in the Upper Floridan aquifer at the study site were significantly different at the 5 percent probability level from those in samples from background wells completed in this aquifer in Leon County. Median potassium concentrations were 0.4 mg/L (milligrams per liter) in water from wells at the study site and 0.7 mg/L in water from background wells in Leon County. Median nitrite plus nitrate concentrations were 6.48 mg/L in water samples from wells at the study site and 0.51 mg/L in water samples from background wells. Graphical comparison of the water quality characteristics in water from six categories of wells; upgradient, landfill, adjacent to the landfill, downgradient onsite, downgradient offsite, and background Leon County, indicated that sodium, bicarbonate, sulfate, iron, manganese, and dissolved solids concentrations and specific conductance values were highest within the landfill. Upgradient and downgradient well categories include wells in areas affected by sludge spreading. Localized reducing conditions were indicated in a downgradient area by median iron concentrations of 930 and 30 ug/L (micrograms per liter) in samples from wells only 5 feet apart. Nitrite plus nitrate concentrations were lowest in water from wells in the landfill compared to water from the other categories. Concentrations of most trace elements and organic compounds were mostly below detection limits, although State maximum contaminant levels of 1 ug/L for benzene and vinyl chloride and 3 ug/L for tetrachloroethene were exceeded in 14 samples. Concentrations of sodium, potassium, chloride, sulfate, nitrite plus nitrate, total organic and ammonia nitrogen, and ammonia nitrogen in ground water were statistically analyzed in areas affected by sludge-spreading and other areas. Only nitrite plus nitrate and chloride concentrations were significantly different between the two areas. Median nitrite plus nitrate and chloride concentrations in sludge-spreading areas were 6.9 mg/L and 2.9 mg/L, respectively, compared to 1.1 mg/L and 1.8 mg/L in other areas.

Florida

The spatially adaptable filter for error reduction (SAFER) process: Remote sensing-based LANDFIRE disturbance mapping updates

LANDFIRE (LF) has been producing periodic spatially explicit vegetation change maps (i.e., LF disturbance products) across the entire United States since 1999 at a 30 m spatial resolution. These disturbance products include data products produced by various fire programs, field-mapped vegetation and fuel treatment activity (i.e., events) submissions from various agencies, and disturbances detected by the U.S. Geological Survey Earth Resources Observation and Science (EROS)-based Remote Sensing of Landscape Change (RSLC) process. The RSLC process applies a bi-temporal change detection algorithm to Landsat satellite-based seasonal composites to generate the interim disturbances that are subsequently reviewed by analysts to reduce omission and commission errors before ingestion them into LF’s disturbance products. The latency of the disturbance product is contingent on timely data availability and analyst review. This work describes the development and integration of the Spatially Adaptable Filter for Error Reduction (SAFER) process and other error and latency reduction improvements to the RSLC process. SAFER is a random forest-based supervised classifier and uses predictor variables that are derived from multiple years of pre- and post-disturbance Landsat band observations. Predictor variables include reflectance, indices, and spatial contextual information. Spatial contextual information that is unique to each contiguous disturbance region is parameterized as Z scores using differential observations of the disturbed regions with its undisturbed neighbors. The SAFER process was prototyped for inclusion in the RSLC process over five regions within the conterminous United States (CONUS) and regional model performance, evaluated using 2016 data. Results show that the inclusion of the SAFER process increased the accuracies of the interim disturbance detections and thus has potential to reduce the time needed for analyst review. LF does not track the time taken by each analyst for each tile, and hence, the relative effort saved was parameterized as the percentage of 30 m pixels that are correctly classified in the SAFER outputs to the total number of pixels that are incorrectly classified in the interim disturbance and are presented. The SAFER prototype outputs showed that the relative analysts’ effort saved could be over 95%. The regional model performance evaluation showed that SAFER’s performance depended on the nature of disturbances and availability of cloud-free images relative to the time of disturbances. The accuracy estimates for CONUS were inferred by comparing the 2017 SAFER outputs to the 2017 analyst-reviewed data. As expected, the SAFER outputs had higher accuracies compared to the interim disturbances, and CONUS-wide relative effort saved was over 92%. The regional variation in the accuracies and effort saved are discussed in relation to the vegetation and disturbance type in each region. SAFER is now operationally integrated into the RSLC process, and LANDFIRE is well poised for annual updates, contingent on the availability of data.

Fire

Groundwater-level elevations in the bedrock aquifers of the Denver Basin aquifer system, Elbert County, Colorado, 2015–23

Water users in Elbert County, Colorado, rely on groundwater from bedrock aquifers in the Denver Basin aquifer system (upper Dawson, lower Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers) for approximately half of their water uses. Withdrawals from the bedrocks aquifers have increased to meet the water use needs of expanding regional population growth and development. The U.S. Geological Survey, in cooperation with the Elbert County Board of County Commissioners, began a study in 2015 to monitor groundwater levels within Elbert County. The primary purpose of this report is to present a summary of groundwater levels measured during the study period (2015–23) and present results from statistical analyses of changes in groundwater-level elevations through time. Discrete groundwater levels were measured at 36 wells within Elbert County. Seven of those wells contained equipment to make and record continuous groundwater-level measurements at hourly intervals. All aquifers, except the lower Dawson aquifer, had only declining groundwater-level elevations in discrete measurements for wells with statistically significant trends. Of the eight statistically significant trends in the lower Dawson aquifer, two wells indicated increasing groundwater-level elevation from discrete measurements. The groundwater-level elevation trend medians in the upper Dawson, lower Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers were −0.23, −0.66, −0.64, −0.39, and −0.63 feet per year, respectively, for discrete groundwater-level elevation measurements. Trends in continuous groundwater-level elevations were in agreement with statistically significant trends in discrete groundwater-level elevations for all wells. The groundwater-level elevation trend medians in this study, compared to the overall trends in a 2015−2018 study, both indicated declining groundwater-level elevations except in the upper Dawson aquifer, where the trend direction was opposite, a positive trend from 2015 to 2018 and a negative trend (declining groundwater elevations) from 2015 to 2023. The change in trends within the upper Dawson aquifer may be affected by differences in the study period and the trend analysis applied. Trends during the 2015–23 study period were compared to departures from the median 2015 groundwater-level elevation for each site in each aquifer. In general, the departures from the 2015 median supported trends observed at each site and correlated spatially with greater departures near the western border of Elbert County. Additionally, 30-year precipitation data showing wet and dry periods were overlaid with the departure from the 2015 median to assess groundwater-level patterns in wells in the five aquifers. Departures from the 2015 median groundwater-level elevations appeared greatest during the dry period between 2020 and 2023. Potentiometric-surface maps of the upper and lower Dawson aquifers created from static April 2023 groundwater elevations indicated groundwater-flow direction is generally from the south to the north. Results of this study could be used to guide additional groundwater monitoring in Elbert County and could aid in long-term planning of water resources.

Colorado

Environmental setting and water-quality issues in the lower Tennessee River basin

The goals of the National Water-Quality Assessment Program are to describe current water-quality conditions for a large part of the Nation's water resources, identify water-quality changes over time, and identify the primary natural and human factors that affect water quality. The lower Tennessee River Basin is one of 59 river basins selected for study. The water-quality assessment of the lower Tennessee River Basin study unit began in 1997. The lower Tennessee River Basin study unit encompasses an area of about 19,500 square miles and extends from Chattanooga, Tennessee, to Paducah, Kentucky. The study unit had a population of about 1.5 million people in 1995. The study unit was subdivided into subunits with relatively homogeneous geology and physiography. Subdivision of the study unit creates a framework to assess the effects of natural and cultural settings on water quality. Nine subunits were delineated in the study unit; their boundaries generally coincide with level III and level IV ecoregion boundaries. The nine subunits are the Coastal Plain, Transition, Western Highland Rim, Outer Nashville Basin, Inner Nashville Basin, Eastern Highland Rim, Plateau Escarpment and Valleys, Cumberland Plateau, and Valley and Ridge.The lower Tennessee River Basin consists of predominantly forest (51 percent) and agricultural land (40 percent). Activities related to agricultural land use, therefore, are the primary cultural factors likely to have a widespread effect on surface- and ground-water quality in the study unit. Inputs of total nitrogen and phosphorus from agricultural activities in 1992 were about 161,000 and 37,900 tons, respectively. About 3.7 million pounds (active ingredient) of pesticides was applied to crops in the lower Tennessee River Basin in 1992. State water-quality agencies identified nutrient enrichment and pathogens as water-quality issues affecting both surface and ground water in the lower Tennessee River Basin. Water-quality data collected by State and Federal agencies between 1980 and 1996 were summarized to characterize surface- and ground-water quality of the subunits with respect to these issues. Median concentrations of nitrogen species generally were less than 1 milligram per liter in surface and ground water in all subunits, and were highest throughout the subunits that had the largest percentages of agricultural land use. Median phosphorus concentrations also were less than 1 milligram per liter in all subunits. Phosphatic limestones present in two subunits had a larger effect on phosphorus concentrations in surface and ground water than did the amount of agricultural land use in these subunits. Median counts of fecal coliform were higher in surface water than in ground water in all subunits. The highest median counts in surface water were in the Valley and Ridge (7,500 colonies per 100 milliliters) and the Outer Nashville Basin subunits (5,000 colonies per 100 milliliters). Highest median counts in ground water were in the Inner and Outer Nashville Basin subunit. Natural setting likely has an important effect with respect to fecal contamination of surface and ground water in the lower Tennessee River Basin.

Alabama, Georgia, Kentucky, Mississippi, Tennessee

Recharge rates and aquifer hydraulic characteristics for selected drainage basins in middle and east Tennessee

Quantitative information concerning aquifer hydrologic and hydraulic characteristics is needed to manage the development of ground-water resources. These characteristics are poorly defined for the bedrock aquifers in Middle and East Tennessee where demand for water is increasing. This report presents estimates of recharge rate, storage coefficient, diffusivity, and transmissivity for representative drainage basins in Middle and East Tennessee, as determined from analyses of stream-aquifer interactions. The drainage basins have been grouped according to the underlying major aquifer, then statistical descriptions applied to each group, in order to define area1 distribution of these characteristics. Aquifer recharge rates are estimated for representative low, average, and high flow years for 63 drainage basins using hydrograph analysis techniques. Net annual recharge during average flow years for all basins ranges from 4.1 to 16.8 in/yr (inches per year), with a mean value of 7.3 in. In general, recharge rates are highest for basins underlain by the Blue Ridge aquifer (mean value11.7 in/yr) and lowest for basins underlain by the Central Basin aquifer (mean value 5.6 in/yr). Mean recharge values for the Cumberland Plateau, Highland Rim, and Valley and Ridge aquifers are 6.5, 7.4, and 6.6 in/yr, respectively. Gravity drainage characterizes ground-water flow in most surficial bedrock aquifer in Tennessee. Accordingly, a gravity yield analysis, which compares concurrent water-level and streamflow hydrographs, was used to estimate aquifer storage coefficient for nine study basins. The basin estimates range from 0.002 to 0.140; however, most estimates are within a narrow range of values, from 0.01 to 0.025. Accordingly, storage coefficient is estimated to be 0.01 for all aquifers in Middle and East Tennessee, with the exception of the aquifer in the inner part of the Central Basin, for which storage coefficient is estimated to be 0.002. Estimates of aquifer hydraulic diffusivity are derived from estimates of the streamflow recession index and drainage density for 75 drainage basins; values range from 3,300 to 130,000 ft^2/d (feet squared per day). Basin-specific and site-specific estimates of transmissivity are computed from estimates of hydraulic diffusivity and specific-capacity test data, respectively. Basin-specific, or areal, estimates of transmissivity range from 22 to 1,300 ft^2/d, with a mean of 240 ft^2/d In general, areal transmissivity is highest for basins underlain by the Cumberland Plateau aquifer (mean value 480 ft^2/d) and lowest for basins underlain by the Central Basin aquifer (mean value 79 ft^2/d). Mean transmissivity values for the Highland Rim, Valley and Ridge, and Blue Ridge aquifer are 320,140, and 120 ft^2/d respectively. Site-specific estimates of transmissivity, computed from specific-capacity data from 118 test wells in Middle and East Tennessee range from 2 to 93,000 ft^2/d with a mean of 2,600 ft^2/d Mean transmissivity values for the Cumberland Plateau, Highland Rim, Central Basin, Valley and Ridge, and Blue Ridge aquifers are 2,800,1,200, 7,800, 390, and 65Oft Id, respectively.

Tennessee

Ground-water flow in the surficial aquifer system and potential movement of contaminants from selected waste-disposal sites at Naval Station Mayport, Florida

Ground-water flow through the surficial aquifer system at Naval Station Mayport near Jacksonville, Florida, was simulated with a two-layer finite-difference model as part of an investigation conducted by the U.S. Geological Survey. The model was calibrated to 229 water-level measurements from 181 wells during three synoptic surveys (July 17, 1995; July 31, 1996; and October 24, 1996). A quantifiable understanding of ground-water flow through the surficial aquifer was needed to evaluate remedial-action alternatives under consideration by the Naval Station Mayport to control the possible movement of contaminants from sites on the station. Multi-well aquifer tests, single-well tests, and slug tests were conducted to estimate the hydraulic properties of the surficial aquifer system, which was divided into three geohydrologic units�an S-zone and an I-zone separated by a marsh-muck confining unit. The recharge rate was estimated to range from 4 to 15 inches per year (95 percent confidence limits), based on a chloride-ratio method. Most of the simulations following model calibration were based on a recharge rate of 8 inches per year to unirrigated pervious areas. The advective displacement of saline pore water during the last 200 years was simulated using a particle-tracking routine, MODPATH, applied to calibrated steady-state and transient models of the Mayport peninsula. The surficial aquifer system at Naval Station Mayport has been modified greatly by natural and anthropogenic forces so that the freshwater flow system is expanding and saltwater is being flushed from the system. A new MODFLOW package (VAR1) was written to simulate the temporal variation of hydraulic properties caused by construction activities at Naval Station Mayport. The transiently simulated saltwater distribution after 200 years of displacement described the chloride distribution in the I-zone (determined from measurements made during 1993 and 1996) better than the steady-state simulation. The advective movement of contaminants from selected sites within the solid waste management units to discharge points was simulated using MODPATH. Most of the particles were discharged to the nearest surface-water feature after traveling less than 1,000 feet in the ground-water system. Most areas within 1,000 feet of a surface-water feature or storm sewer had traveltimes of less than 50 years, based on an effective porosity of 40 percent. Contributing areas, traveltimes, and pathlines were identified for 224 wells at Naval Station Mayport under steady-state and transient conditions by back-tracking a particle from the midpoint of the wetted screen of each well. Traveltimes to contributing areas that ranged between 15 and 50 years, estimated by the steady-state model, differed most from the transient traveltime estimates. Estimates of traveltimes and pathlines based on steady-state model results typically were 10 to 20 years more and about twice as long as corresponding estimates from the transient model. The models differed because the steady-state model simulated 1996 conditions when Naval Station Mayport had more impervious surfaces than at any earlier time. The expansion of the impervious surfaces increased the average distance between contributing areas and observation wells.

Florida

Introduction and methods of analysis for peak streamflow trends and their relation to changes in climate in Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin

Flood-frequency analysis, also called peak-flow frequency or flood-flow frequency analysis, is essential to water resources management applications including critical structure design and floodplain mapping. Federal guidelines for doing flood-frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C is that for drainage basins without major hydrologic alterations, statistical properties of the distribution of annual peak streamflows (peak flows) are stationary; that is, the mean, variance, and skew are constant. The stationarity assumption has been widely accepted within the flood-frequency community; however, a better understanding of long-term climatic persistence and concerns about potential climate change and land-use change has caused a reexamination of the stationarity assumption. Flood-frequency analyses that do not incorporate observed trends and abrupt changes may result in a poor representation of the true flood risk. Bulletin 17C does not offer guidance for incorporating nonstationarities when estimating floods, and it describes a need for studies that incorporate changing climate or basin characteristics. In response to this need and a history of concern regarding nonstationarity peak flows in the region, this study was done by the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, to assess potential nonstationarity in peak flows in the north-central United States. This chapter summarizes the methods used to detect hydroclimatic changes in peak-flow data in the study region. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flow. Four periods were selected for analysis of peak flow, daily streamflow, and climate data. The periods are (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. The climate data consist of monthly time series estimates of temperature, precipitation, potential evapotranspiration, actual evapotranspiration, snowfall, soil moisture storage, snow water equivalent, and runoff on a 3.1-mile by 3.1-mile grid for the conterminous United States. Statistical and graphical analyses were used to investigate potential changes in hydrology and climate. The starting point for these analyses was the initial data analysis of peak flow described in Bulletin 17C, which includes plotting the peak flow and checking for autocorrelation, monotonic trends, and changes points. Analyses were added to examine additional features in the data. To examine potential causal drivers of changes, the climate data were analyzed graphically and statistically. Results are provided in a U.S. Geological Survey data release. The study limitations are documented for users of the results.

Iowa, Illinois, Michigan, Missouri, Minnesota, Mon

Refining the earthquake history of south-central Alaska through lake records

The Alaska–Aleutian subduction zone (AASZ) is one of the world's most seismically active plate boundaries and the source of the 1964 M w 9.2 Great Alaska earthquake–the second largest instrumentally recorded earthquake in the world. Understanding the nature and frequency of such earthquakes is necessary for seismic and tsunami hazard assessment, but instrumental and historical records that span less than 150 years are too short to allow a statistically reliable analysis of earthquake recurrence times. This calls for studies of evidence of past earthquakes, extending the earthquake catalog further back in time. Subduction-zone paleoseismology in south-central Alaska is predominantly based on coastal evidence of land-level changes and tsunamis generated by megathrust earthquakes and preserved in the geological record. A complementary approach is lacustrine paleoseismology, which is still a relatively young discipline in Alaska. However, globally, lake basins are well-established high-resolution and continuous recorders of paleoseismic activity along subduction zones, relying on the identification of underwater landslide deposits and turbidites generated by seismic shaking. As a result, lake basins not only register ground shaking from megathrust earthquakes, but also from intraslab and crustal earthquakes, which are typically not accompanied by significant land-level changes. In this review paper, we combine coastal and lacustrine paleoseismology approaches to refine the south-central Alaskan earthquake history by comparing the paleoseismic records from two lakes (i.e., Eklutna Lake, located in the Chugach Mountain Range, and Skilak Lake, situated on the Kenai Peninsula) with the coastal and crustal earthquake catalog in Alaska. The resulting age ranges of all known megathrust earthquakes involving the Alaskan megathrust between the Kodiak and Prince William Sound (PWS) sections are more precise and accurate for the last 1.3 kyrs BP than the previously published age ranges from coastal records. As a result, this study supports the following key conclusions: (1) The 1964 CE earthquake was an exceptionally strong and unique event in the last 2000 years, rupturing the PWS, Kenai, Barren Islands, and Kodiak sections simultaneously. (2) The high-resolution and seasonal markings of the varved lake records now disentangle for the first time closely timed earthquakes, which was not possible based on the coastal evidence alone. (3) No persistent megathrust rupture boundaries exist. So, the possibility of a full rupture of the entire eastern AASZ, from PWS to Semidi cannot be excluded. (4) The rupture pattern in the eastern AASZ reveals superimposed cycles of multi-asperity ruptures (1964 earthquake) and clustered complementary partial ruptures, or rupture cascades. (5) The PWS section hosts the largest asperity in the eastern AASZ. (6) The shaking record of megathrust earthquakes indicates a time-dependent (quasiperiodic) behavior for the study area, but the observation of complementary clusters means that the hazards will not drop to zero but instead may even increase for a neighboring section. (7) The time-independent behavior of intraplate earthquakes implies that the intraslab hazard did not decrease following the 2016 and 2018 earthquakes. This study utilizes an integrated approach for subduction zone paleoseismology as a solution for unraveling recurrence and rupture patterns in Alaska, which can be applied worldwide.

Alaska

Comparison of two approaches for determining ground-water discharge and pumpage in the lower Arkansas River Basin, Colorado, 1997-98

In March 1994, the Colorado Division of Water Resources (CDWR) adopted “Rules Governing the Measurement of Tributary Ground Water Diversions Located in the Arkansas River Basin” (Office of the State Engineer, 1994); these initial rules were amended in February 1996 (Office of the State Engineer, 1996). The amended rules require users of wells that divert tributary ground water to annually report the water pumped monthly by each well. The rules allow a well owner to report the pumpage measured by a totalizing flowmeter (TFM) or pumpage determined from electrical power data and a power conversion coefficient (PCC) (Hurr and Litke, 1989). Opinions by representatives of the State of Kansas, presented before the Special Master hearing a court case [State of Kansas v. State of Colorado, No. 105 Original (1996)] concerning post-Compact well pumping, stated that the PCC approach does not provide the same level of accuracy and reliability as a TFM when used to determine pumpage. In 1997, the U.S. Geological Survey (USGS), in cooperation with the CDWR, began a 2-year study to compare ground-water pumpage estimates made using the TFM and the PCC approaches. The study area was along the Arkansas River between Pueblo, Colorado, and the Colorado-Kansas State line (fig. 1). The two approaches for estimating ground-water discharge and pumpage were compared for more than 100 wells completed in the alluvial aquifer of the Arkansas River Basin. The TFM approach uses an inline flowmeter to directly measure instantaneous discharge and the total volume of water pumped at a well. The PCC approach uses electrical power consumption records and a power conversion coefficient to estimate the pumpage at ground-water wells. This executive summary describes the results of the comparison of the two approaches. Specifically, (1) the differences in instantaneous discharge measured with three portable flowmeters and measured with an inline TFM are evaluated, and the statistical differences in paired instantaneous discharge between the two approaches are determined; (2) short- and long-term variations in the PCC’s are presented; (3) differences in pumpage between the two approaches are evaluated, and the statistical differences in pumpage between the two approaches are determined; (4) potential sources of discrepancy between pumpage estimates are discussed; and (5) differences in total network pumpage using the two approaches are presented. During the irrigation seasons of 1997 and 1998, instantaneous discharge and electrical power demand were measured at randomly selected wells to determine PCC’s. At more than 100 wells, the PCC’s determined during the 1998 season were applied to total electrical power consumption data that was recorded between the initial and final readings at each network well site in 1998 to estimate total ground-water pumpage. At each site, an inline TFM was installed in a full-flowing, acceptable test section of pipe on the discharge side of the pump where the measurement of discharge was made. Measurements of instantaneous ground-water discharge also were made using three different types of portable flowmeters. The average velocity multiplied by the cross-sectional area of the discharge pipe was used to compute the discharge in gallons per minute. Whenever possible, discharge measurements were made at each network site using all three types of portable flowmeters.

Colorado

Timing and geometry of the Chemehuevi Formation reveal a late Pleistocene sediment pulse into the Lower Colorado River

The Chemehuevi Formation is a distinctive 50−150-m-thick wedge-shaped Pleistocene sedimentary unit deposited by the Colorado River. It lines the perimeters of the river’s floodplains and bedrock canyons for more than 600 km between the mouth of the Grand Canyon and the delta region in the Gulf of California. The formation is composed of a basal tan to light-yellowish-brown and pale-orange mud-dominated facies overlain and interbedded by a light-yellow-brown sand-dominated facies. The unit is one of two extensively exposed aggradational packages in the Lower Colorado River corridor, in addition to a series of other smaller alluvial terrace deposits. The Chemehuevi Formation appears to represent the response of a fully integrated Colorado River system to a significant perturbation, in contrast to the Bullhead Alluvium, which is likely a unique result of Pliocene river integration. The aggradation of the Chemehuevi Formation in the Lower Colorado River corridor may be similarly due to a unique event in the Colorado River system, or it may instead be a well-preserved sedimentary sequence recording typical behavior of the Colorado River below the Grand Canyon in the late Pleistocene. As such, multiple causal mechanisms have been proposed, but no study to date has conclusively explained the Chemehuevi Formation. To help resolve its timing, duration, and origin, we applied post-infrared infrared stimulated luminescence, carbonate U-Th series, and zircon sensitive high-resolution ion microprobe U-Th series geochronology to determine the ages of key exposures of the unit over a wide spatial area. These new data demonstrate that the Chemehuevi Formation was deposited ca. 110−90 ka. The depositional ages collectively overlap, suggesting that deposition occurred rapidly relative to the resolution of the geochronometers. The new depositional timing coincides with a shift from glacial to interglacial conditions after the marine isotope stage 5-6 transition. This observation is consistent with a climate-induced sediment pulse as a causal mechanism, yet correlations with similar deposits in the Colorado River headwaters or in neighboring catchments appear elusive. Potentially, climate transitions between glacial and interglacial periods induced a sediment pulse from hillslopes of the Colorado River system that resulted in the Chemehuevi Formation. An alternative or additional explanation is that the Chemehuevi Formation represents release of lava dam−impounded sediment in the Grand Canyon. The surface geometry of the Chemehuevi Formation projects upstream to the approximate location of lava dams, and the largest possible lava dam impoundment (the Upper Prospect dam) is comparable in volume to the formation. The lava dam hypothesis appears to be a possible explanation for the Chemehuevi Formation. However, tying deposition to a specific lava dam or series of lava dams remains challenging due to discrepancies in timing and volume. The combined effects of a series of lava dams may have led to the Chemehuevi Formation, as the last Pleistocene lava dam eruption coincides with the onset of deposition. Alternatively, the formation may result from the combined effects of both regional climate transitions and the lava dams that created a transient reservoir to compound a climate transition−driven sediment pulse. The geochronologic data presented here do not allow us to distinguish between the lava dam or climate transition hypotheses but will need to be reconciled with any future proposed depositional model.

Arizona, California, Nevada

Spatial distribution of API gravity and gas/oil ratios for petroleum accumulations in Upper Cretaceous strata of the San Miguel, Olmos, and Escondido Formations of the south Texas Maverick Basin—Implications for petroleum migration and charge history

The Maverick Basin of south Texas is currently undergoing active exploration and production of gas and oil from tight sandstone reservoirs. The most productive tight sandstones in the basin are in the Upper Cretaceous San Miguel, Olmos, and Escondido Formations. These units are second only to the Eagle Ford Shale in terms of cumulative production volumes. The structural history of the Maverick Basin, from rifting to subsidence to exhumation, has had a profound effect on the characteristics of these reservoirs and the petroleum resources contained therein. This U.S. Geological Survey review of the production history of these strata reflects a recent shift from conventional production to horizontal drilling (unconventional) that exploits low permeability reservoirs in previously overlooked areas of existing oil and gas fields in southern Texas, typically outside of established field boundaries. To investigate the physical properties of the Maverick Basin hydrocarbon accumulations, this case study compiled American Petroleum Institute (API) gravity measurements and calculated cumulative gas/oil ratios (GOR) for thousands of producing wells from the San Miguel, Olmos, and Escondido Formations. Maps were generated from the compiled well production data to show the spatial heterogeneity of API gravity and GOR values for the three formations within the Maverick Basin and immediately outside the basin to the northeast. Within the Maverick Basin, the spatial patterns of API gravity values indicate lighter oils downdip towards the southern basin edge. GOR values indicative of wet and dry gases within the basin are seen interspersed, with values that correspond to black and heavy oils. Differences in the spatial patterns of the petroleum properties within the Maverick Basin are interpreted as effects of Eocene basin inversion caused by Laramide orogenic deformation, and the resulting reservoir exhumation of basin strata. East of the Maverick Basin, spatial distributions of API gravity and GOR values show progressively heavier oils updip to the northwest, grading to dry gases downdip to the southeast, which correlates to the oil and gas windows of the underlying Eagle Ford Shale. Correlation of API gravity and GOR values from the San Miguel, Olmos, and Escondido Formations with thermal maturity data from the Eagle Ford Shale suggests that the Eagle Ford Shale is the petroleum source, and that petroleum migration was approximately vertical for areas to the east of the Maverick Basin. The discontinuity of API gravity and GOR properties within the Maverick Basin implies a complex petroleum charge history, possibly involving the remigration of petroleum and the addition of petroleum from other source intervals in Mexico, to the southwest. Depressurization of exhumed, overpressured reservoirs of the San Miguel, Olmos, and Escondido Formations can explain the intermittent occurrence of gas production throughout the southern Maverick Basin by exsolution of gas from formation brines and the resulting dry gas flushing of hydrocarbon-charged reservoirs. The introduction of dry gas through flushing can, in turn, explain why the patterns of API gravity and GOR values are so dissimilar in the Maverick Basin. This process has implications for possible future production of unconventional resources from undiscovered tight-gas reservoirs in strata of the San Miguel, Olmos, and Escondido Formations, and a different approach to petroleum exploration may be needed in the Maverick Basin relative to exploration techniques applied in other basins within the northern Gulf of Mexico.

Texas