USGS ScienceSearch

SEARCH · USGS Science

Results for “Global Change Biology”

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.

552 records · Page 30Linked to original sources

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

Spatio-temporal evolution of distributed volcanic fields, case studies—Sierra Chichinautzin and Michoacán-Guanajuato, México

An analysis of 1,375 volcanoes in the Michoacán-Guanajuato (1,148 volcanoes in a 26,200 square-kilometer area) and Sierra Chichinautzin (227 volcanoes in a 3,500 square-kilometer area) volcanic fields in central Mexico identified patterns in the spatial and temporal distribution of past eruptions. A cluster agglomerative hierarchical method and kernel analysis confirmed that the Michoacán-Guanajuato volcanic field comprises four volcanic fields (Valle de Santiago, Uruapan, Apatzingán, and Pátzcuaro volcanic fields) controlled by different fault systems, indicating that it is not a single volcanic field but rather a group of volcanic fields (a “superfield”), each of which has distinct characteristics. In the Sierra Chichinautzin volcanic field, well-constrained isotopic ages were used to build a model of how the spatial distribution of the eruptions has changed over time. Two new 40Ar/39Ar ages from a locally recognized volcanic feature near the town of El Cantil, herein called El Cantil volcano (1,537±17 kilo-annum [ka]) and the volcanic feature at Cerro el Elefante (herein called El Elefante dome) (1,485±92 ka) belong to the oldest volcanic group identified in the Sierra Chichinautzin volcanic field, confirming the timing of the beginning of monogenetic volcanism in the region. Based on the volcanic groups identified in the Sierra Chichinautzin volcanic field, the youngest volcanism (less than 35 ka) is found only in the central-western sector of the field. Principal component analysis determined the directional trends of feeder dikes only for vents <10 ka in the Sierra Chichinautzin volcanic field. Possible magma migration paths through the crust were identified using seismic data from both volcanic fields using an earthquake catalog from 1973 to 2023, which includes 9,016 earthquakes in the Michoacán-Guanajuato volcanic field and 841 in the Sierra Chichinautzin volcanic field. The spatial distribution of the hypocenters does not highlight any trend that could be associated with superficial movement of magma in the Sierra Chichinautzin volcanic field. In the Michoacán-Guanajuato volcanic field, however, eight seismic swarms since 1997 have been detected. These swarms are interpreted to result from ascending magma. Strengthening monitoring systems and reinforcing mitigation measures to address volcanic hazards and risk are important means of preparing for future eruptions in both regions. Analysis such as those herein provide insights as to where an eruption might occur and may help mitigate volcanic hazards.

Michoacán-Guanajuato, Sierra Chichinautzin

Sackung at Bald Eagle ridge, central Colorado: An updated interpretation of ridge-spreading movement, structures, and mechanisms from 50 years (1975–2025) of U.S. Geological Survey research

Slow gravitational failures of mountain peaks and ridges are poorly understood. Herein, we report on 50 years of studies at a slowly spreading castellate ridge in the Sawatch Range in central Colorado. The orientations of geomorphic-structural features indicate that the fractured Precambrian granitic rock underlying the ridge has extended and spread northwestward toward the formerly glacier-covered Busk Creek valley. Results from surveying, field-based geomorphic-structural mapping using lidar, rock mass quality measurements, a passive seismic survey, and satellite radar provide a major update to research started by U.S. Geological Survey researchers in the 1970s and 1980s. New insights include a recognition that the entire ridge has slowly moved by concurrent sliding along an inferred northwest dipping, compound basal-slip surface (or zone), and through the formation of multiple grabens by normal faulting and flexural toppling along sets of pre-existing fractures that dip moderately (∼45°) to the southeast. We were unable to distinguish the presence of a sudden and strong contrast in seismic velocity across the inferred slip surface. Movement during the 50-year study period has been episodic and gradually decreasing, in correspondence with decreasing cumulative annual precipitation and increasing mean annual air temperatures. The fastest moving area, just upslope from the glacier trimline, had an average horizontal velocity of 3 – 4 mm/yr. Evidence suggests that movement started as a paraglacial response mechanism, but because of the site’s proximity to the Rio Grande Rift, we cannot exclude earthquake shaking as a mechanism for initiation or enhancement of slope movement. An estimate of longer-term horizontal movement from the exposed basal-slip surface at the uphill side of the ridgetop graben is ∼1.1 mm/yr for the 13 – 14 ky post-glacial period. Broad implications of our work are that: (1) long-term measurements (decades or longer) of slope movement can add insights into how sackungen form and evolve through time; (2) the identification of thrust faults and toes in zones of compression near valley bottoms can be crucial for interpreting sackung failure mechanisms, and (3) the use of passive-seismic techniques to identify the depth to a slip surface may not be successful in granitic terrain dominated by planar fractures and subtle changes in rock-mass characteristics.

Colorado

Overview of the Cenozoic geology of the northern Harrat Rahat volcanic field, Kingdom of Saudi Arabia

The Harrat Rahat volcanic field, located in the west-central part of the Kingdom of Saudi Arabia, is one of the larger Cenozoic harrats among the more than 17 harrats situated upon the Arabia Plate. The map plate contained herein shows, at a scale of 1:100,000, the mapped volcanic geology of northern Harrat Rahat, which consists of the northernmost one-fifth of Harrat Rahat. Northern Harrat Rahat has an area of about 3,340 square kilometers (km 2 ), of which 2,567 km 2 is covered by Harrat Rahat volcanic rocks and deposits, and it encompasses more than 900 exposed vents (that is, craters, cryptodomes, maars, and scoria cones), 289 of which are isolated by younger volcanic deposits and have not been correlated with the 234 volcanic rock units distinguished by geologic mapping. Northern Harrat Rahat is of special interest owing to its proximity to the city of Al Madinah al Munawwarah, which sits within, and is continuing to expand southward over, the north end of the volcanic field. Al Madinah is home to an expanding population, currently at more than 2 million residents, together with the intermittent addition of approximately 3 million pilgrims during Hajj and Umrah (religious visitations). The center of Al Madinah is less than 8 km from lava flows of the only confirmed historically documented eruption, which occurred in 1256 C.E. (654 A.H.). Earlier prehistoric lava flows also encroached into the area of the present-day city limits, as demonstrated by volcanic rocks exposed widely throughout the city in roadcuts, parks, and excavations for new buildings, although no evidence has been found of any other than the 1256 C.E. lava having reached that area during times of human habitation. Geologic mapping was undertaken by the U.S. Geological Survey in collaboration with the Saudi Geological Survey. The features of primary interest within the map area are scoria cones, lava flows, lava domes, craters, and pyroclastic deposits from the Quaternary, which have compositions of basalt, hawaiite, mugearite, benmoreite, and trachyte. The geologic mapping was published by Downs and others (2019) on a single sheet at 1:75,000 scale and two detailed sheets at 1:25,000 scale, accompanied by detailed explanations of the geology. This report presents the geology on a single sheet at 1:100,000 scale and provides condensed geologic explanations for the convenience of readers of this volume. Some minor errors of nomenclature and lava source regions that do not change fundamental interpretations have been corrected herein.

Professional Paper

ShakeAlert® version 3: Expected performance in large earthquakes

The ShakeAlert earthquake early warning (EEW) system partners along with U.S. Geological Survey (USGS) licensed operators deliver EEW alerts to the public and trigger automated systems when a significant earthquake is expected to impact California, Oregon, or Washington. ShakeAlert’s primary goal is to provide usable warning times before the arrival of damaging shaking. EEW is most likely to achieve this goal in large‐magnitude earthquakes. In recent years, ShakeAlert has gone through a series of upgrades to its underlying scientific algorithms aimed at improved performance during large earthquakes. Version 3 of this software recently went live in the production system and includes improvements to all algorithms. The main seismic algorithms that detect an earthquake and characterize its location, magnitude, and fault rupture orientation are faster than older versions. Other key changes include: using real‐time geodetic data to characterize the magnitude growth in large earthquakes; the introduction of an alert pause procedure to compromise between speed near the epicenter and improved accuracy at larger distances; and the inclusion of a nonergodic site‐response model in the ground‐motion predictions. ShakeAlert has achieved its primary goal of usable warning times before strong shaking at some locations in real‐time operations in recent M 6 earthquakes. Using offline tests, we demonstrate usable warning times are possible for many sites with peak shaking values of modified Mercalli intensity (MMI) 7–8 in M 7+ earthquakes and also for many MMI 8–9 sites in M 8+ earthquakes. ShakeAlert partners use a variety of MMI and magnitude thresholds in deciding when to alert their users within bounds set by the USGS. Our study shows that there is room to raise the magnitude thresholds up to about M 5.5 without adversely affecting performance in large earthquakes. The ground‐motion criteria are more complex owing to a significant drop‐off in warning times between the MMI 4 and 5 levels of predicted shaking. However, widely used ShakeAlert products, such as the MMI 3 and 4 contour products, can provide sufficiently long warning times before strong shaking in moderate‐to‐great earthquakes to enable a range of protective actions.

Bulletin of the Seismological Society of America

Water-withdrawal and groundwater-level fluctuations, physical and infiltration characteristics of Ohio River sediment, and isotope-based estimates of water sources relative to a riverbank filtration system in an alluvial aquifer, Louisville, Kentucky, 2000–20

The U.S. Geological Survey and the Louisville Water Company, from 2000 to 2020, cooperatively investigated physical and hydrologic conditions in the Ohio River alluvial aquifer and adjacent riverbed sediments in the northeast part of Louisville, Kentucky, to increase knowledge of processes affecting water sources to collector wells during operation of a riverbank filtration system (RBF). The RBF, in 2020, included an initial collector well (CW1) completed in 1999 and a tunnel and collector well system (TCWS) composed of four collector wells (CW2–CW5) completed in late 2010 that were connected by a tunnel that routed water to a pumping station and a treatment plant. Daily mean water-withdrawal rates typically ranged from about 10.4 to 20.7 million gallons per day (Mgal/d) from 2005 to 2007, then were cyclic from 2008 to 2010 and ranged from as low as about 9 to 13 Mgal/d in winter-spring months up to about 19 to 20 Mgal/d in summer-autumn months. From 2011 through 2017, peak daily mean water-withdrawal rates from CW2–CW5 fluctuated between winter-spring lows of about 20–30 Mgal/d to summer-autumn peak rates of about 41.6 Mgal/d in 2013 to about 48.2 Mgal/d in 2012. From 2017 to 2020, TCWS peak daily mean water-withdrawal rates decreased from about 45.7 Mgal/d in 2017 to about 37.5 Mgal/d in 2020 and peak daily mean rates from CW1 ranged from about 20.7 Mgal/d in 2020 to about 21.6 Mgal/d in 2018. Periods in 2019–20 with tunnel pumping station daily mean water-withdrawal rates smaller than about 30 Mgal/d had generally stable water levels in CW2 and indicated the sustainable yield of the TCWS collector wells under the prevailing hydrologic conditions. Drawdown during July–August of 2008 and 2009 was associated with CW1, the sole operating collector well. Groundwater-flow directions outside the CW1 cone of depression during July–August of 2008 and 2009 were from areas near the alluvial valley margin toward the Ohio River and in 2009 toward Harrods Creek. Collector wells CW2 through CW5 derived water during 2011–19 surveys from surface-water infiltration, as indicated by drawdown that extended beneath the Ohio River and Harrods Creek. Increases in drawdown toward the alluvial valley margins in July–August of 2013 and 2015, compared with July 2011, indicated less available recharge from those flow directions and decreased yield of groundwater from the aquifer after progressively longer TCWS water withdrawals. Drawdown at CW2–CW5 and at observation wells between them lessened in October 2019 and indicated that onshore groundwater levels can recover when collector well withdrawals decrease. Overlapping cones of depression in 2011–15 potentiometric surfaces centered around CW2, CW3, CW4, and CW5, and drawdown along divides between those wells indicated well interference between collector wells. Well-interference effects increased during 2013 and 2015 survey periods, despite decreased water-withdrawal rates compared with July 2011, and expanded to include CW1 when its use resumed. Groundwater-level fluctuations in observation wells near the Ohio River were consistent with gradual depletion of groundwater yield during tunnel pumping station operation at river elevations of about 420 feet above the National Geodetic Vertical Datum of 1929 (NGVD 29) and enhanced recharge of the aquifer when river elevations were above flood stage of about 430 feet above NGVD 29. Groundwater-level declines in observation wells near the Ohio River were most apparent in summer and some autumn months between 2012 and 2017 during periods when tunnel pumping station mean daily water-withdrawal rates were about 35 Mgal/d or more. Periods of relatively stable groundwater levels during 2019–20 in observation wells between CW2–CW5, when tunnel pumping station water-withdrawal rates were about 30–33 Mgal/d or less, also indicated that the prevailing sustainable yield was about 30 Mgal/d or less. Wells associated with CW3 and CW4 had the largest increases in slopes of best-fit lines relating water-level difference to water-withdrawal rate between 2012 and 2020, an indication of decreased efficiencies of those wells to supply water or of changes in groundwater yield around the well. Increased slopes of best-fit lines relating water-level difference to water-withdrawal rate from 2019 to 2020, when groundwater levels were otherwise relatively stable, also indicated that the specific capacity of one or more TCWS collector wells may have declined during that period. Riverbed surface elevation, roughness, and hardness changes between hydroacoustic surveys indicated active accretion and erosion of bed sediment offshore from the riverbank filtration system in 2002–4. Riverbed hydraulic conductivity computed from 2013 to 2015 in tests at 11 sites ranged from 0.2 to 230 feet per day. Ten temporary piezometer sites with adjacent operating collector wells had drawdowns that ranged from about 4 to 11.7 feet and hammer blows greater than or equal to 49 counts, which indicated a relatively great density of alluvial sediments and the potential for diminished riverbed infiltration rates. Contributions of Ohio River infiltration to collector wells estimated from deuterium and oxygen-18 isotopic ratios in representative Ohio River and groundwater samples were considerably larger in warmer months than cooler months and ranged from 80–90 percent in CW1–CW5 in May 2020 to 100 percent in May 2019 at CW2, CW4, and CW5. Groundwater contributions to water withdrawals by collector wells increased in cooler autumn months and ranged from 40 percent at CW3 and CW5 in November 2020 to 70 percent at CW2, CW4, and CW5 in November–December 2019. November 2019 and November 2020 isotopic compositions of water from the collector wells had enriched deuterium and oxygen-18 isotopic ratios consistent with mixtures of evaporation-modified water from the Ohio River and meteoric-origin groundwater.

Indiana, Kentucky

Benthic habitat map of Olowalu Reef, Maui, Hawaii—Geomorphological structure, biological cover, and geologic zonation determined with spectral, lidar, and acoustic data

The fringing coral reef off Olowalu, Maui, Hawaii, has been identified as a local conservation priority site. In 2007, the National Oceanic and Atmospheric Administration (NOAA) produced a benthic habitat map of the Hawaiian Islands that was used as a foundation for this study. To support place-based management of the reef in the future, the U.S. Geological Survey (USGS) mapped the geologic zone, major and dominant geomorphological structure, biological cover type, and percent of biological cover for 11 square kilometers (km 2 ) of Olowalu Reef at a minimum mapping unit (MMU) of 100 square meters (m 2 ) to create a benthic habitat map. Heads-up digitization was employed on 0.50-meter (m) natural color satellite orthoimagery with ancillary 1-m acoustic backscatter imagery from single-scan sonar (sound navigation and ranging). A 1-m, 4-m, and 8-m digital bathymetric model (DBM) was interpolated from bathymetric lidar (light detection and ranging), and various geomorphometric layers derived from the DBMs were used for habitat interpretation. Still-frame imagery of the seafloor extracted from vessel-towed underwater video transects on Olowalu Reef served as ground validation points ( n =870) during active mapping and accuracy assessment points ( n =216) for thematic accuracy assessment. Thematic accuracy was cross-validated by the Hawai‘i Department of Land and Natural Resources Division of Aquatic Resources. Final thematic accuracy was 88.8 percent for major structure, 85.6 percent for dominant structure, 86.0 percent for major biological cover, and 78.6 percent for type and percent of major biological cover. Reef and hardbottom constituted 52 percent of the total mapped habitat, comprising mostly aggregate reef (31 percent) and pavement (11 percent), with large swaths of spur-and-groove (9 percent). Of this hardbottom, 17 percent was covered with moderate (10 to <50 percent) coral and 27 percent with high coral cover (50 to <90 percent). High (50 to <90 percent) macroalgae cover dominated the continuous sand sheets in offshore bank/shelf zones. The map created in this study supplements the NOAA 2007 map and expands on the observations made by USGS sampling of the reef. The NOAA 2007 map and our map differed in total areal extent by a negligible 6 m 2 and were in general thematic agreement. Our map is intended to serve as a baseline for public access, general research, local-level management, and reef change for future studies.

Hawaii

Phase 1 technical implementation plan for the expansion of the ShakeAlert earthquake early warning system to Alaska

Executive Summary The conference report accompanying the fiscal year (FY) 2022 Consolidated Appropriations Act (Public Law 117–103) for the U.S. Department of the Interior and related agencies directed the U.S. Geological Survey (USGS) to “work with the State of Alaska to develop an implementation plan to be completed within two years in order to put ShakeAlert/Earthquake Early Warning in Alaska” (p. 29). Congress included $1 million in the FY 2022 appropriation to conduct this effort. The USGS Earthquake Hazards Program, along with partner organizations, has developed the ShakeAlert earthquake early warning (EEW) system for the West Coast, which currently operates in California, Oregon, and Washington. The purpose of the system and its alert delivery partners is to reduce the impact of earthquakes and save lives and property by delivering ShakeAlert-powered alerts that are transmitted to the public via mass notification technologies, and by providing more detailed data streams to institutional users and commercial service providers to trigger automated, user-specific, protective actions. ShakeAlert was designed in such a way that it could be expanded to other U.S. regions with high earthquake risk, after the build-out of seismic and geodetic networks to support ShakeAlert in a specified region is completed and the necessary funding is secured for long-term operation and maintenance. When an earthquake occurs, seismic waves radiate from the rupturing fault like waves on a pond. It is these waves that people feel as earthquake shaking and that can cause damage to structures. Using networks of ground-motion sensors and sophisticated computer algorithms, ShakeAlert can detect an earthquake seconds after it begins, calculate its location and magnitude, and estimate the resulting intensity of shaking. Early warnings of impending shaking are then sent to people and systems that may experience damaging shaking, allowing them to take appropriate protective actions. Depending on the user’s distance from the earthquake, alerts may be delivered before, during, or after the arrival of strong shaking. There will almost always be a region near the earthquake epicenter where alerts arrive after damaging shaking has begun. The ShakeAlert system updates its ground-motion estimates as an earthquake grows larger. In response to the FY 2022 congressional direction, the USGS worked with the State of Alaska to devise this implementation plan for ShakeAlert expansion to Alaska. The USGS engaged with the Alaska Division of Homeland Security and Emergency Management (DHS&EM) and the Alaska Division of Geological and Geophysical Surveys (DGGS). A cooperative agreement was awarded to the Alaska Earthquake Center (AEC) at the University of Alaska Fairbanks (UAF) for their contributions to the plan and their work coordinating with other networks in Alaska. The USGS engaged with the Alaska Seismic Hazards Safety Commission (ASHSC) throughout the process. The USGS also held a series of Alaska stakeholder engagements. The process of developing the implementation plan was facilitated by contracted staff from Corner Alliance, which is a government consulting firm. This implementation plan describes the details and estimates the costs for a Phase 1 expansion of the ShakeAlert system to Alaska. A geographically limited Phase 1 goal was chosen that covers the highest risk and most populated areas of Alaska. The areas proposed encompass the State’s main population centers and 90 percent of the State’s population. This Phase 1 design is considered very challenging and ambitious from the viewpoint of network operators. The lessons learned if this plan is implemented could be used to consider subsequent phases to expand EEW beyond Phase 1 in Alaska in the future. ShakeAlert is built on the foundation of the sensor networks and data processing infrastructure of the USGS-led Advanced National Seismic System (ANSS). This implementation plan calls for a total of 450 high-quality, real-time EEW-capable ANSS seismic stations in Alaska: 270 new stations, 160 upgraded stations, and 20 existing stations. These seismic station numbers are based on a station spacing of 10 kilometers (km) in urban areas, 20 km in seismic source areas that endanger population centers, and 40 km in other areas. The associated costs also include support for some EEW-capable global navigation satellite system (GNSS) stations, with a focus on improving warnings for large subduction zone earthquakes. For effective EEW, ShakeAlert requires low-latency, high-availability, robust telemetry links to deliver continuous, real-time data from field stations to the data centers. The Alaska data processing hardware infrastructure would follow the general design for fail-safe operation that is used for the ShakeAlert system on the West Coast. The ShakeAlert architecture uses two independent layers: the production layer for earthquake processing and the alert layer to make alerting decisions and serve alerts to users. This implementation plan includes two geographically separated data centers in Alaska, each with two fully independent production and alert layers using the same system design developed for the West Coast. As of March 2024, the ShakeAlert system is at version 3.0.1, with more advanced versions in the development and testing pipeline. ShakeAlert originally used two algorithms to determine the location and magnitude of earthquakes using seismic data. A third algorithm that can calculate very large magnitudes of very large earthquakes with geodetic data was added in March 2024. ShakeAlert publishes several data and alert products to meet the needs of different users. All messages include the location of the earthquake, either as a point or a line, and its magnitude. Ground-shaking estimates are published in two forms, as ground-motion contours and a map grid. Providing adequate warning time for strong shaking (the “target threshold”) requires sending alerts at a threshold lower than that strong shaking level (the “alert threshold”). The thresholds for public alerting in Alaska would be a joint USGS and State decision. To have the greatest benefit, ShakeAlert-powered alerts would be delivered to institutional users and individuals by all practical pathways. The USGS alert layer can support thousands of institutional users and alert redistributors, but the USGS does not have the mission nor the infrastructure and expertise to perform mass notifications to the public or implement automatic actions for end users of the alerts. To meet this need, ShakeAlert recruits private sector “technology enablers” that have the necessary expertise to develop end-user implementations using EEW alerts with the goal of stimulating an EEW industry. Earthquake early warning alerts are useless if people do not know how to respond to them. Although the alert messages include instructions about what to do (drop, cover, and hold on), alerts are more effective if people have been trained in advance. Messages about ShakeAlert’s capabilities, limitations, and benefits could be integrated with existing earthquake education programs, including State-run programs. Therefore, ShakeAlert would coordinate with both public and private partners and stakeholders through various partnerships and agreements to accomplish consistent and ongoing public earthquake hazard education. The estimated capital cost of completing the computing infrastructure and sensor networks for the Phase 1 ShakeAlert expansion to Alaska is approximately $66 million in 2024 dollars. The annual operation and maintenance cost of the completed system is estimated to be $12 million per year in 2024 dollars when fully built out.

Alaska

Permafrost–wildfire interactions: active layer thickness estimates for paired burned and unburned sites in northern high latitudes

As the northern high-latitude permafrost zone experiences accelerated warming, permafrost has become vulnerable to widespread thaw. Simultaneously, wildfire activity across northern boreal forest and Arctic/subarctic tundra regions impacts permafrost stability through the combustion of insulating organic matter, vegetation, and post-fire changes in albedo. Efforts to synthesis the impacts of wildfire on permafrost are limited and are typically reliant on antecedent pre-fire conditions. To address this, we created the FireALT dataset by soliciting data contributions that included thaw depth measurements, site conditions, and fire event details with paired measurements at environmentally comparable burned and unburned sites. The solicitation resulted in 52 466 thaw depth measurements from 18 contributors across North America and Russia. Because thaw depths were taken at various times throughout the thawing season, we also estimated end-of-season active layer thickness (ALT) for each measurement using a modified version of the Stefan equation. Here, we describe our methods for collecting and quality-checking the data, estimating ALT, the data structure, strengths and limitations, and future research opportunities. The final dataset includes 48 669 ALT estimates with 32 attributes across 9446 plots and 157 burned–unburned pairs spanning Canada, Russia, and the United States. The data span fire events from 1900 to 2022 with measurements collected from 2001 to 2023. The time since fire ranges from 0 to 114 years. The FireALT dataset addresses a key challenge: the ability to assess impacts of wildfire on ALT when measurements are taken at various times throughout the thaw season depending on the time of field campaigns (typically June through August) by estimating ALT at the end-of-season maximum. This dataset can be used to address understudied research areas, particularly algorithm development, calibration, and validation for evolving process-based models as well as extrapolating across space and time, which could elucidate permafrost–wildfire interactions under accelerated warming across the high-northern-latitude permafrost zone. The FireALT dataset is available through the Arctic Data Center ( https://doi.org/10.18739/A2RN3092P , Talucci et al., 2024).

Arctic

Wind River subbasin restoration: Annual report of U.S. Geological Survey activities January 2024 through December 2024

We sampled juvenile wild Oncorhynchus mykiss (Steelhead Trout) in headwater streams of the Wind River, Washington, to characterize population attributes and investigate life-history metrics, particularly migratory patterns, and early life-stage survival. Additionally, this sampling provides an index of abundance of non-native Brook Trout Salvelinus confluentus . We used passive integrated transponder (PIT) tagging and a series of instream PIT-tag interrogation systems (PTISs) to track juveniles and adults. The Wind River subbasin is considered a wild Steelhead refuge by Washington Department of Fish and Wildlife (WDFW). No hatchery Steelhead Trout have been released in the Wind River subbasin since 1997, and hatchery adults are estimated at less than one percent of spawners in most years. Over 25 years of Steelhead Trout status and trend monitoring and research in the subbasin is contributing to understanding of population response to numerous restoration actions in the subbasin, including removal of Hemlock Dam from Trout Creek in 2009, which had an inadequate adult ladder and contributed to increased water temperatures, to near lethal levels for salmonids. Data from our study, and companion work by WDFW, are contributing to the Columbia Basin Fish and Wildlife Program (2009) Research, Monitoring, and Evaluation (RM&E) Strategy of Fish Population Status Monitoring. Specifically, this work addresses the sub-strategies of 1) Assessing the Status and Trends of Diversity of Natural Origin Fish Populations and Uncertainties Research regarding differing life histories of a wild Steelhead Trout population, 2) Assessing the Status and Trend of Adult Natural Origin Fish Populations, and 3) Monitoring and Evaluating the Effectiveness of Tributary Habitat Actions Relative to Environmental, Physical, or Biological Performance Objectives. During summer and fall 2024, we PIT-tagged 2,075 Steelhead Trout parr (age-0 and age-1), in the Trout Creek and upper Wind River watersheds. Age-0 and age-1 parr were at high abundance due to strong spawner numbers in 2023 and 2024 (Estimate of 932 in 2023; and a moderate number of spawners in 2024 though an estimate is not available at this time; 24-year median = 619; Charlie Cochran, WDFW Fish Biologist, personal commun., 2025). An additional 100 parr (age-0 and age-2 or older) were tagged to provide fish for estimating detection efficiencies at PTISs. Steelhead Trout parr were recaptured and detected through repeat headwater sampling, smolt trapping, instream PTISs and Columbia River mainstem PIT-tag detection. We maintained and upgraded six instream PTISs to detect PIT-tagged Steelhead Trout parr, smolts, and adults, providing data for population assessments, and life-cycle research. Detection data from PIT-tagged adult Steelhead Trout at our four primary PTISs allow assessment of adult escapement to tributary watersheds within the Wind River subbasin. Most adult Steelhead detected were tagged by WDFW at an adult trap in the lower Wind River and some as smolts at screw traps. Adult Steelhead Trout detection efficiency estimates at our primary PTIS in Trout Creek have been greater than 99 percent during eight of the past ten years and have exceeded 97% at our primary PTIS in the Wind River during nine of the past ten years. Adult escapement estimates to tributary watersheds are helping evaluate the efficacy of the 2009 removal of Hemlock Dam from rkm 2.0 of Trout Creek. Detections at the instream PTISs have shown trends of age-0 and age-1 Steelhead Trout parr emigration from natal areas during summer and fall, in addition to the expected movement of parr and smolts in spring. Our data suggest that most fish from a cohort that migrate downstream will do so at age-1 for additional rearing in the subbasin. We have estimated that from 15% to 56% of parr tagged as age-0 fish in headwater areas make downstream migrations at age-1 for additional rearing (median = 29%). We have estimated that from 1% to 27% of Steelhead Trout parr, tagged as age-1 fish, make downstream migrations during fall (median = 3.5%): this is especially pronounced in the upper Wind River portion of the watershed (median 10%). These findings raise questions about parr rearing habitat and whether migrations are density- or habitat-quality driven, and answers to such questions are long-term goals of this study as active and passive habitat restoration occurs. Repeat sampling at sites in the subbasin within and between years has enabled assessment of juvenile Steelhead Trout growth patterns. Growth rates (relative change in weight) of age-0 PIT-tagged parr during summer have been similar across the subbasin, though slightly lower in the Trout Creek watershed (Trout Creek sites, 0.57 - 1.25%/day, Wind River sites = 0.94 – 2.48%/day). Summer growth rates have been lower for age-1 parr in the Trout Creek watershed than the upper Wind River watershed (Trout Creek sites, = 0.03 - 0.14%/day, Wind River sites = 0.05 – 0.19%/day). Yearly relative growth was similar across the subbasin for both age-0 and age-1 tagged parr. During 2024 range of mean daily relative change in weight from the previous year for age-0 tagged fish was: Trout Creek sites = 0.96 – 1.41%/day, Wind River sites = 0.80 – 1.30 %/day; and for age1 tagged fish was: Trout Creek sites = 0.37 - 0.50 %/day, Wind River sites = no age-1 tagged recaptures. Non-native Brook Trout are present in the subbasin, chiefly the Trout Creek watershed. Mean percent-of-catch that is Brook Trout, at four sample sites in Trout Creek, has declined from the period 1998 – 2003 to the period 2011 – 2024 (range of decrease = 19 - 82%). Percent-of-catch (of 64 samples, 52 < 10%; 3 > 20%) and number of Brook Trout at the Trout Creek sites from 2011 through 2024, though variable, have generally declined. Habitat restoration evaluation efforts are critical to ensure efficient use of money and resources but require long-term monitoring of population metrics. Detailed viable salmonid population and life-history data, such as that provided by PIT-tagging and instream PTIS networks inform fisheries policy and management and enable assessment of long-term effects of habitat restoration actions such as the removal of Hemlock Dam on Trout Creek and proposed major instream habitat restoration in the upper Wind River.

Washington

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

Preserving and increasing water resources—Natural infrastructure in dryland streams in Baja California Sur, Mexico

The Los Planes watershed of Baja California Sur, Mexico, and its underlying aquifer are experiencing groundwater decline owing to low average annual rainfall (28.1 centimeters per year) and rising water demand from population growth and agricultural activities. This decline in water availability can lead to desertification—a process that changes arable land to desert by degrading soil and vegetation—and can pose serious challenges to livelihoods that depend on the land. To address these issues, a ranch in the Los Planes watershed has installed many natural infrastructures in dryland streams (NIDS) in channels for soil and water conservation. In 2022, the U.S. Geological Survey (USGS) began working with regional researchers and land managers to investigate the effects of NIDS on natural biological, geochemical, and physical processes and determine the efficacy of NIDS for water augmentation in the Los Planes watershed. The USGS also worked with local academic institutions and nonprofit organizations to create public educational opportunities focused on the area’s hydrogeology. These and other collaborative efforts with the U.S. Water Partnership and Innovaciones Alumbra aim at enhancing water resources in the Baja California Sur region and promoting water security and safeguarding community well-being. Resume La cuenca de Los Planes, ubicada en Baja California Sur, México, y su acuífero subyacente, están sufriendo una disminución de las aguas subterráneas debido a la baja precipitación media anual (28.1 centímetros por año) y la alta demanda de agua por parte de una población creciente y la actividad agrícola. Esta disminución de la disponibilidad de agua puede conducir a la desertificación—un proceso que por medio de la degradación del suelo y la vegetación convierte a la tierra cultivable en desierto—representando un serio desafío para los medios de vida de las personas. Para abordar estos problemas, un rancho en la cuenca de Los Planes ha instalado numerosas obras de Infraestructura Natural en Arroyos de Tierras Áridas (INATS) para conservación del suelo y del agua. En 2022, el Servicio Geológico de los Estados Unidos (USGS, por sus siglas en inglés) comenzó a trabajar con investigadores regionales y gestores de tierras para estudiar los efectos de INATS en los procesos biológicos, geoquímicos y físicos, y determinar su eficacia en el aumento de los recursos hídricos en la cuenca de Los Planes. El USGS se ha asociado con instituciones académicas y organizaciones locales sin fines de lucro para crear oportunidades educativas públicas centradas en la hidrogeología de la zona. Estos y otros esfuerzos colaborativos con la Asociación del Agua de Estados Unidos (U.S. Water Partnership) e Innovaciones Alumbra, tienen como objetivo mejorar el uso de los recursos hídricos en la región de Baja California Sur, promover la seguridad hídrica y proteger el bienestar de la comunidad.

Baja California Sur

Eruptive history of northern Harrat Rahat—Volume, timing, and composition of volcanism over the past 1.2 million years

Harrat Rahat, one of several large, basalt-dominated volcanic fields in the western part of the Kingdom of Saudi Arabia, is a prime example of continental, intraplate volcanism. Excellent exposure makes this an outstanding site to investigate changing volcanic flux and composition through time. We present 93 40 Ar/ 39 Ar ages and 6 36 Cl surface-exposure ages for volcanic deposits throughout northern Harrat Rahat that, integrated with a new geologic map, define 12 eruptive stages. Exposed volcanic deposits in the study area erupted less than 1.2 million years ago (Ma), and 214 of 234 identified eruptions occurred less than 570 thousand years ago (ka). Two eruptions were in the Holocene, including a historically described basaltic eruption in 1256 C.E. and a trachyte eruption newly recognized as Holocene (4.2±5.2 ka). An estimated approximately 82 cubic kilometers (km 3 ; dense rock equivalent) of volcanic products can be documented as having erupted since 1.2 Ma, though this is a lower limit because of concealment of deposits older than 570 ka. Over the last 570 thousand years (k.y.), the average eruption rate was 0.14 cubic kilometers per thousand years (km 3 /k.y.), but volcanism was episodic with periods alternating between low (0.04–0.06 km 3 /k.y.) and high (0.1–0.3 km 3 /k.y.) effusion rates. Before 180 ka, eruptions vented from the volcanic field’s dominant eastern vent axis and from a subsidiary, diffuse, western vent axis. After 180 ka, volcanism focused along the eastern vent axis, and the composition of volcanism varied systematically along its length from basalt dominated in the north to trachyte dominated in the south. We hypothesize that these compositional variations younger than 180 k.y. reflect the growth of a mafic intrusive complex beneath the southern part of the vent axis, which led to the development of evolved magmas. Lastly, these new age data allow for a reassessment of the volcanic recurrence interval at northern Harrat Rahat. Based on available data, volcanism in northern Harrat Rahat over the last 180 k.y. is poorly described using a Poisson distribution with a single recurrence interval. Instead, data for northern Harrat Rahat are better described using a mixed exponential distribution that is applicable for volcanic systems characterized by two different eruptive states, where one state with a longer recurrence interval corresponding to periods of low eruption frequency and one state with a shorter recurrence interval corresponding to periods of high eruption frequency. The preferred model for northern Harrat Rahat over the last 180 k.y. uses a long recurrence interval of 4.0 k.y. and a short recurrence interval of 0.22 k.y.

Professional Paper

Chemostratigraphy of the Cretaceous Hue Shale in Arctic Alaska: Exploring paleoceanographic controls on trace element enrichment, organic matter accumulation, and source-rock evolution

We document chemostratigraphy in an outcrop of late Albian to early Campanian (∼103–82 Ma) marine source rocks to better understand paleoenvironmental controls on trace element (TE) enrichment and organic matter accumulation in the distal Colville foreland basin of Arctic Alaska and how those drivers are linked to arc volcanism and successions of Cretaceous oceanographic and climatic biogeochemical events. This unique, 113-m-thick section of Cretaceous Hue Shale deposited during a series of previously undocumented Arctic Cretaceous oceanic anoxic events (Lease et al., 2024) is the only known exposure of thermally immature (0.48–0.52% R o , random vitrinite reflectance) Hue Shale in Arctic Alaska. Strata comprise mainly clay-rich mudstone with elevated total organic carbon (TOC) and hydrogen index values reaching 26.3 wt% (mean = 7.5 wt%) and 689 mg hydrocarbon (HC)/g TOC (mean = 385 mg HC/g TOC), respectively. Maceral composition consists predominantly of fluorescent amorphous organic matter, with abundant brightly fluorescent alginite, including Tasmanites , acritarchs, and Leiosphaeridia . Discrete layers of volcanic ash (preserved as bentonite) are present throughout the section and provide quantitative age control based on U–Pb dates. Chemostratigraphic trends are investigated to advance our understanding of local oceanographic conditions and controls on stratigraphic and temporal heterogeneity of Brookian source rocks. Concurrent sedimentary enrichment in Mo, U, V, Pb, and Cu across the Albian–Cenomanian boundary of the exposed basal gamma-ray zone, may reflect anoxic to euxinic benthic redox conditions favoring organic matter accumulation and preservation. Fluctuating degrees of anoxia-euxinia are inferred throughout the overlying Hue Shale succession, reflected by varying patterns of TE enrichment and TE–TOC covariation. Elevated C org /P molar ratios (>250) across most of the section, with several values exceeding 690, signify that enhanced biological productivity is sustained throughout deposition. Enhanced productivity, recorded by both C org /P and excess Ba, also parallels increases in source rock richness (elevated TOC and S 2 values) during the late Albian–early Cenomanian and late Cenomanian – Turonian. Enhanced productivity and variations in oceanic circulation/stratification likely both drove changes in benthic redox conditions that favored organic carbon accumulation and preservation. Increased continental arc volcanism (e.g., Okhotsk–Chukotka volcanic belt) and High Arctic Large Igneous Province magmatic eruptions throughout the Cretaceous, inferred to have influenced nutrient cycling and local aqueous nutrient availability, also have been invoked as potential drivers of organic carbon burial and source-rock development across the sedimentary sequence. Results presented here document the organic-rich and oil-prone source-rock quality of the Hue Shale in the distal part of the Colville foreland basin and bolster the potential for a Cretaceous petroleum system beneath the eastern North Slope.

Alaska

Methodology and technical input for the 2025 U.S. List of Critical Minerals—Assessing the potential effects of mineral commodity supply chain disruptions on the U.S. economy

The Secretary of the Interior, acting through the Director of the U.S. Geological Survey, is tasked by section 7002 (“Mineral Security”) of title VII (“Critical Minerals”) of the Energy Act of 2020 (Public Law 116–260, December 27, 2020, 116th Congress) with reviewing and revising the methodology used to evaluate mineral commodity supply risk and the U.S. List of Critical Minerals (LCM) no less than every 3 years. Following two previous LCM assessments, this analysis represents the latest technical input for evaluating each mineral commodity’s supply risk and determining their recommended status on the LCM. We evaluated mineral commodity supply risk using two criteria: (1) an economic effects assessment that quantified the potential effects of various trade disruption scenarios on the U.S. economy, and (2) an examination of whether the mineral commodity’s U.S. supply chain relied on a sole domestic producer that represented a single point of failure. For the first criterion, postdisruption equilibrium quantities and prices for each mineral commodity were calculated based on their price elasticities of supply and demand and the availability of excess production capacity for each yearlong foreign trade disruption scenario. Subsequently, a nonlinear optimization routine was used with detailed economic input-output tables to estimate the potential economic effects on the U.S. economy of over 1,200 scenarios for 84 mineral commodities. After accounting for the probability of each scenario’s occurrence, the overall results are presented in terms of changes in U.S. gross domestic product (GDP) by individual industry and the economy overall. The results, which ranged from a net decrease in U.S. GDP of nearly $4.5 billion to a net increase of $33 million, largely reflect U.S. import dependency and world production concentration. Using the Jenks natural breaks optimization method, a statistical classification technique, we categorized the mineral commodities into several classes based on this overall risk quantification. Mineral commodities with annualized probability-weighted net decreases in U.S. GDP greater than $2 million were recommended for inclusion on the LCM. If a mineral commodity did not meet the threshold for inclusion on the LCM under the first criterion, its domestic supply chain was examined under the second criterion, which recommended a mineral commodity for inclusion on the LCM if there was only a single domestic producer. Ultimately, the two criteria resulted in the recommendation of the addition of six mineral commodities (in descending risk order, potash, silicon, copper, silver, rhenium, and lead) to and the removal of two mineral commodities (arsenic and tellurium) from the LCM. By using an economic effects assessment, the results of this analysis provide a prioritization that can also be compared directly against other risk analyses and the cost of various risk mitigation strategies.

Open-File Report

Precipitation-based flood-inundation maps for the East Fork Little Blue River and tributaries at Lee’s Summit, Missouri, 2024

The U.S. Geological Survey, in cooperation with the City of Lee’s Summit, Missouri, assessed flooding of the East Fork Little Blue River and tributaries for varying precipitation magnitudes and durations, varying antecedent runoff conditions, and projected climate-change conditions. The precipitation scenarios were used to develop a library of flood-inundation maps for a 2.95-mile reach of the East Fork Little Blue River and tributaries within the city. A two-dimensional U.S. Army Corps of Engineers Hydrologic Engineering Center–River Analysis System (HEC–RAS; ver. 6.5) rain-on-grid model was calibrated to selected runoff events representing a range of antecedent runoff conditions and hydrologic responses. Lowest adjacent grades for structures within the nearby study area were incorporated into the terrain, and depth grids and water-surface elevation grids were developed for the study area. Simulated velocities at selected bridge locations were also developed from the model. The model was calibrated using water-surface elevation data collected from water-level loggers (pressure transducers) and streamflow measurements and water-surface elevation measurements made at a reference point during runoff events. The calibrated HEC–RAS model was used to simulate streamflows from design rainfall events of 15-minute to 24-hour durations and ranging from a 100- to 0.1-percent annual exceedance probability (1-year to 1,000-year recurrence intervals). Flood-inundation maps were produced for depths at a reference location of 3 to 16 feet, or a depth exceeding the 0.1-percent annual exceedance probability interval precipitation. The results of each precipitation duration-frequency value were represented by a 1-foot-increment inundation map based on the generated peak streamflow from that rainfall event and the corresponding water-surface elevation at the East Fork Little Blue River reference location. Within the HEC–RAS model, 240 scenarios were developed from the design rainfall events with each of 3 antecedent conditions. Additional scenarios were created to simulate the effects of projected precipitation scenarios on the 100-year recurrence interval, 24-hour storm and the 100-year recurrence interval, 6-hour storm. All simulation results were assigned to a flood-inundation map condition based on the generated peak flow and corresponding water-surface elevation at the East Fork Little Blue River reference location. The flood-inundation maps are shown on a web mapping application made available to the public through the City of Lee’s Summit (hyperlink will be added when available). The flood-inundation maps are tied to real-time precipitation data obtained from the Automated Surface Observing System weather station at the Lee’s Summit Municipal Airport, accessible at https://mesonet.agron.iastate.edu/request/download.phtml?network=MO_ASOS . The availability of these maps, along with information regarding observed rainfall, could help provide emergency management personnel and residents with information that is critical for flood-response activities, such as evacuations and road closures, and for postflood recovery efforts.

Missouri

Three-dimensional temperature maps of the Williston Basin, USA: Implications for deep hot sedimentary and enhanced geothermal resources

As part of U.S. Geological Survey's (USGS) efforts to identify and assess geothermal energy resources of the US, a three-dimensional (3D) geologic and thermal model has been constructed for the Williston Basin, USA. The geologic model consists of all sedimentary units above the Proterozoic and Archean crystalline rock (called basement herein), with a total sedimentary thickness of up to 5 km near the basin center. Twenty-nine geologic units were mapped from interpreted formation tops from 16,465 wells. A 3D temperature model was constructed to a depth of 7 km by constructing a 3D heat flow model for the sedimentary units, followed by estimating underlying temperature using a one-dimensional (1D) analytic solution for heat flow within the underlying crystalline basement. Using the sedimentary basin model, heat flow was simulated in 3D and was calibrated using three temperature datasets: 1) 24 high-confidence static temperature logs (equilibrium thermal profiles), 2) more than15,000 drill stem test (DST) measurements from >7,000 wells, and 3) more than 45,000 bottomhole temperature (BHT) measurements from >14,000 wells. The DST and BHT datasets provide broad spatial coverage, but are lower confidence, primarily because measurements were made prior to attaining thermal equilibrium. DST and BHT measurements were binned regionally to develop representative thermal profiles that generally agree with these lower quality data (hereafter called pseudowell temperature profiles). Layer properties (primarily thermal conductivity and compaction curves) were set to best estimate values, then the heat flow model was calibrated to fit pseudowell and static temperature logs primarily by adjusting basal heat flow to approximate the overall temperature profile. Minor adjustments to thermal conductivity allowed adjusting changes in slope at lithologic contacts. Resulting maps include 3D temperature and basal (bottom of sedimentary units) heat flow estimates, which are used as input for the temperature model of the basement. The crystalline basement temperature model uses an analytic 1D solution to the heat flow equation that requires estimates of heat flow and temperature at the upper boundary (i.e., the sediment/basement contact), radiogenic heat production within the crystalline basement, and reference thermal conductivity (i.e., uncorrected for temperature). Two regions of high heat flow are identified: 1) in western North Dakota along the North American Central Plains Conductivity Anomaly and 2) in eastern Montana near the Poplar dome. Within the sedimentary column in the center of the basin of the basin, an area of approximately 100,000 km2 is predicted to have moderate- to high-temperature geothermal resources (>90 °C) under the thickest sequences of sediments. Where thick insulation and high heat flow coincide, electric-grade resources can be less than 4 km deep. Assuming a maximum feasible drilling depth of 7 km, temperatures are predicted to be as high as 175 °C. The geologic model may be used to identify strata at sufficient temperatures that may have natural permeability or that may have conditions that favor development of enhanced/engineered geothermal systems resources.

Montana, North Dakota, South Dakota

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