USGS ScienceSearch

SEARCH · USGS Science

Results for “Applied Computing and Geosciences.”

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,656 records · Page 46Linked to original sources

High-precision earthquake catalog for Minto Flats fault zone, central Alaska, reveals complex and conjugate faulting

The Minto Flats fault zone (MFFZ) in central Alaska is a left‐lateral strike‐slip fault system situated between the continental‐scale right‐lateral Denali and Kaltag‐Tintina faults. The MFFZ has the potential to generate magnitude 7 earthquakes, and it hosted a magnitude 6 earthquake in 1995. It has also produced exotic events, such as very‐low‐frequency earthquakes and nucleation signals. We use network‐matched filtering and relative earthquake relocation techniques to derive a detailed catalog of earthquake locations for the MFFZ. The catalog spans from August 2014 to December 2019, a time period including 13 temporary seismic stations in the region. Our results provide the most complete catalog for the MFFZ and include deeper events, clusters of shallow seismicity, and a complex and segmented fault structure not observed in the original regional catalog. We document right‐lateral strike‐slip faulting, conjugate to the main northeast‐striking left‐lateral faults of the MFFZ. Below Nenana basin, the relocated seismicity reveals northwest‐dipping left‐lateral faults, supporting the inference that deep crustal active faulting is associated with recent basin deformation.

Alaska

Geothermal district energy systems coupled with seasonal underground thermal energy storage: A U.S. techno-economic screening by climate and geology

In the United States, cooling-dominated commercial building loads can cause geothermal heat pump-based district energy systems to accumulate a long-term subsurface thermal imbalance, motivating the incorporation of seasonal underground thermal energy storage. We developed a transferable workflow to evaluate geothermal district systems that pair ground heat exchangers with seasonal underground thermal energy storage. Using standardized hourly loads for seven commercial buildings and a uniform cost framework, we simulated ten U.S. cities with a physics-based ground heat exchanger model, subsurface storage simulations, and economic assessment to isolate the roles of climate and hydrogeology. In cooling-dominated cities, underground thermal energy storage supplied the majority of annual cooling, cutting electricity use and summer peaks substantially while achieving levelized costs comparable to or below conventional chiller-boiler plants. In cooler climates, the storage share shrunk, required borefield size and costs rose, and levelized cost of energy increased nearly linearly with declining underground thermal energy storage fraction, indicating storage fraction as the primary economic lever. Sensitivity analysis showed capital risk dominated by borefield drilling and surface heating, ventilation, and air-conditioning and piping, with underground thermal energy storage costs secondary. This workflow provides a transparent foundation for site-specific design and screening of next-generation geothermal district energy systems.

Arizona, Illinois, Maryland, Michigan, Mississippi

Hydrogeologic conceptualization of the Arroyos subarea and Northeast Air Force Research Laboratory subarea of the Air Force Research Laboratory, Edwards Air Force Base, Antelope Valley, California, 2017–21

Research and testing within the Air Force Research Laboratory at Edwards Air Force Base in the southwestern Mojave Desert has been conducted since the 1950s. This testing requires large quantities of groundwater for cooling and operations involving hazardous constituents that have entered the groundwater. To better understand the potential movement of contaminants through the groundwater system, the U.S. Geological Survey entered into a cooperative agreement with the U.S. Air Force Civil Engineer Center with the purpose of (1) gaining a more detailed understanding of the hydraulic connection between the unconsolidated basin-fill sediment and the underlying weathered and crystalline bedrock and (2) investigating faults and their potential effect on groundwater flow. This study was specific to the Arroyos subarea and Northeast Air Force Research Laboratory subarea of the Air Force Research Laboratory. The refined hydrogeologic characterization was derived by evaluating surface geologic maps, a regional gravity model, new subsurface geophysical surveys, lithology data from boreholes, and groundwater-level data. Refinement of the hydrogeologic characterization of the study area involved developing a map of basin-fill sediment thickness, estimating the geometry and depth to bedrock, determining the presence of weathered bedrock material, and identifying the location of faults and subsurface structures. The thickness of basin-fill sediment was derived from a regional gravity model. Basin-fill sediment ranges from less than 20 feet (ft) thick to about 150 ft thick and increases substantially toward a deep structural basin in the northwestern part of the study area. The top of bedrock estimated from the new surface geophysical surveys is generally deeper than the top of bedrock from the regional gravity model. These differences in interpreted bedrock tops indicate that there may be a zone of weathered bedrock in the subsurface, overlying crystalline bedrock. Additional subsurface data evaluating the geometry and thickness of weathered bedrock are necessary for a comprehensive evaluation of weathered bedrock in the study area. Previously unmapped extensions of the Leuhman Fault, Spring Fault Zone, and Boron Gate Fault were identified in the new surface geophysical surveys. Evaluation of groundwater flow between bedrock and basin-fill sediment, and the effects of faults on groundwater flow, was done using groundwater-level data from 2020. The groundwater table was in crystalline bedrock throughout most of the Arroyos and Northeast Air Force Research Laboratory subareas, with a few exceptions where it was in basin-fill sediment and weathered bedrock. Basin-fill sediment, and likely weathered bedrock, was partially saturated in parts of the study area where basin-fill sediment was as thick as 150 ft. Weathered bedrock material may be present near the Spring Fault Zone, and if so, was likely saturated in 2020. The groundwater-level data near and within Spring Fault Zone indicate that the Spring Fault Zone acts as a vertical-along-strike conduit for, and as an across-strike impediment to, groundwater flow. In addition, steep groundwater-level gradients in the southeastern part of the study area indicate that a subsurface structure affecting groundwater flow could be related to the Leuhman Fault. Sufficient groundwater-level data were not available to conclusively determine if the Arroyos, Rich, and Boron Gate Faults, and their inferred extensions from the geophysical surveys, affect the movement or direction of groundwater flow.

California

Using gridded seismicity to forecast the long-term spatial distribution of earthquakes for the 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

Gridded (or background) seismicity models are a critical component of probabilistic seismic hazard assessments, accounting for off‐fault and smaller‐magnitude earthquakes. They are typically developed by declustering and spatially smoothing an earthquake catalog to estimate a long‐term seismicity rate that can be used to forecast future earthquakes. Here, we present new gridded seismicity models for use in the 2025 National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The previous PRVI NSHM was released in 2003, and our new models incorporate updates to both data and methodology. We utilize an updated earthquake catalog based on improved Puerto Rico Seismic Network data with newly characterized completeness epochs. The catalog is divided into crustal, subduction interface, and intraslab seismicity using new methods and Slab2 subduction zone geometries. To forecast the long‐term spatial distribution of earthquakes, we use an updated methodology developed for the 2023 U.S. 50‐state NSHM, considering three declustering methods and two spatial smoothing methods based on 2D Gaussian kernels. To adapt it for the complex seismotectonics of the region, we also adopt probabilistic methods to account for events with unknown depths and uncertainties in tectonic classification, and develop a new method for spatial scaling to counteract the effects of spatial variability in network coverage while maintaining the use of smaller events. Finally, we test the performance of these spatial models in forecasting the location of M w ≥ 5earthquakes in the region. Our updated methodology improves the representation of epistemic uncertainty relative to the 2003 model, and our results demonstrate the effectiveness of the new measures we have introduced to address heterogeneities in network detection and systematically evaluate forecast performance.

Puerto Rico, U.S. Virgin Islands

Agricultural tile drains increase the susceptibility of streams to longer and more intense streamflow droughts

Streamflow droughts are receiving increased attention worldwide due to their impact on the environment and economy. One region of concern is the Midwestern United States, whose agricultural productivity depends on subsurface pipes known as tile drains to improve trafficability and soil conditions for crop growth. Tile drains accomplish this by rapidly transporting surplus soil moisture and shallow groundwater from fields, resulting in reduced watershed storage. However, no work has previously examined the connection between tile drainage and streamflow drought. Here, we pose the question: does the extent of watershed-level tile drainage lead to an increased susceptibly and magnitude of streamflow droughts? To answer this, we use daily streamflow data for 122 watersheds throughout the Midwestern United States to quantify streamflow drought duration, frequency, and intensity. Using spatial multiple regression models, we find that agricultural tile drainage generates statistically significant ( p < 0.05) increases in streamflow drought duration and intensity while significantly reducing drought frequency. The magnitude of the effect of tile drainage on streamflow drought characteristics is similar to that of water table depth and precipitation seasonality, both of which are known to influence streamflow droughts. Furthermore, projected changes in regional precipitation characteristics will likely drive the installation of additional tile drainage. We find that for each 10% increase in tile-drained watershed area, streamflow drought duration and intensity increase by 0.03 d and 12%, respectively, while frequency decreases by 0.10 events/year. Such increases in tile drainage may lead to more severe streamflow droughts and have a detrimental effect on the socio-environmental usage of streams throughout the Midwest.

Illinois, Indiana, Iowa, Michigan, Minnesota, Ohio

Quantitative risk of earthquake disruption to global copper and rhenium supply

Earthquakes have the potential to substantially affect mining operations, potentially leading to supply chain disruptions and adversely affecting the global economy. This study explores the quantification of earthquake risk to copper and rhenium commodity supply by examining the spatial concentration of high earthquake hazard areas and the commodity-specific mining, smelting, and refining operations across the globe. Because many of the largest facilities are concentrated geographically near the highly seismic regions of South America, East Asia, and the Pacific, there is a potential for cascading effects on the entire supply chain. The analysis indicates that the expected annual disruption of global production is 0.3–1.1 percent for copper mines, 1.8–4.0 percent for smelters, and 1.5–3.3 percent for refineries. Expected annual disruption of global rhenium production capacity is 0.32–1.32 percent. The research highlights that the potential lost revenue from earthquake disruptions is from $315 million to $1.29 billion for copper mining, $1.92 billion to $4.33 billion for copper smelting, $2.06 billion to $4.52 billion for copper refining, and $337,000 to $1.40 million for rhenium production capacity.

Open-File Report

Tapwater exposures, residential risk, and mitigation in a PFAS-impacted-groundwater community

Tapwater (TW) safety and sustainability are priorities in the United States. Per/polyfluoroalkyl substance(s) (PFAS) contamination is a growing public-health concern due to prolific use, widespread TW exposures, and mounting human-health concerns. Historically-rural, actively-urbanizing communities that rely on surficial-aquifer private wells incur elevated risks of unrecognized TW exposures, including PFAS, due to limited private-well monitoring and contaminant-source proliferation in urbanizing landscapes. Here, a broad-analytical-scope TW-assessment was conducted in a hydrologically-vulnerable, Mississippi River alluvial-island community, where PFAS contamination of the shallow-alluvial drinking-water aquifer has been documented, but more comprehensive contaminant characterization to inform decision-making is currently lacking. In 2021, we analyzed 510 organics, 34 inorganics, and 3 microbial groups in 11 residential and community locations to assess (1) TW risks beyond recognized PFAS issues, (2) day-to-day and year-to-year risk variability, and (3) suitability of the underlying sandstone aquifer as an alternative source to mitigate TW-PFAS exposures. Seventy-six organics and 25 inorganics were detected. Potential human-health risks of detected TW exposures were explored based on cumulative benchmark-based toxicity quotients ( ∑ TQ ). Elevated risks ( ∑ TQ ≥ 1) from organic and inorganic contaminants were observed in all alluvial-aquifer-sourced synoptic samples but not in sandstone-aquifer-sourced samples. Repeated sampling at 3 sites over 52–55 h indicated limited variability in risk over the short-term. Comparable PFAS-specific ∑ TQ for spatial-synoptic, short-term (3 days) temporal, and long-term (3 years quarterly) temporal samples indicated that synoptic results provided useful insight into the risks of TW-PFAS exposures at French Island over the long-term. No PFAS detections in sandstone-aquifer-sourced samples over a 3 year period indicated no PFAS-associated risk and supported the sandstone aquifer as an alternative drinking-water source to mitigate community TW-PFAS exposures. This study illustrated the importance of expanded contaminant monitoring of private-well TW, beyond known concerns (in this case, PFAS), to reduce the risks of a range of unrecognized contaminant exposures.

Wisconsin

Cascadia Subduction Zone science: Call for the next generation community seismic velocity model

The Cascadia subduction zone (CSZ) hosts major seismic and tsunami hazards, yet key questions persist about the relationship between margin structure, fluid distribution, episodic tremor and slip, shallow megathrust behavior, shaking and tsunamigenesis, and the resulting hazard estimates. Addressing these problems requires an empirically grounded, three‐dimensional seismic velocity model to illuminate subsurface structure and properties and to provide a basis for geophysical studies such as earthquake simulations and ground‐motion estimation. In May 2024, the National Science Foundation‐funded Cascadia Region Earthquake Science Center (CRESCENT) community velocity model (CVM) working group, with U.S. Geological Survey and regional partners, convened a workshop to identify priorities for such a model. Participants emphasized the features necessary for addressing key science questions, including implementing findability, accessibility, interoperability, and reusability (FAIR) access, capturing along‐strike and along‐dip structural heterogeneity, resolving shallow offshore–onshore structure, constraining elastic properties and quantifying their uncertainties for numerical wave propagation simulations, their validation benchmarks, and supporting associated accurate earthquake ground‐motion simulations and hazard assessments. This article describes the priorities defined in the workshop, and a description of how, guided by these needs, CRESCENT plans to develop multiple generations of a CVM to advance CSZ science and improve seismic and tsunami hazard modeling across the Pacific Northwest. The CVM will span the CSZ from the surface to ∼100 km depth, offshore and east of the Cascades into Idaho (∼132°–110° W) and the southern and northern tectonic regime transitions (∼36°–52° N) to capture the entire tectonic system as well as its surroundings.

Cascadia Subduction Zone

Geologic map of MTM −15032 and −20032 quadrangles, western Ladon basin, Mars

Our mapping efforts focused on the geomorphic terrains and geologic units contained within Mars Transverse Mercator (MTM) −15032 and −20032 quadrangles. These two quadrangles are located along the west side of Ladon basin and span lat −12.5° N. to −22.5° N. and long 325° E. to 330° E. The western part of Ladon basin and its bounding basin ring structures to the west preserved features that help to understand the long history of drainage across the Margaritifer Terra region of Mars. Our geologic map provides new insight into the extent of aqueous activity across western Ladon basin and the adjoining highlands, which includes the northern part of Ladon Valles, Arda Valles, numerous small valleys in the western highlands, and light-toned layered deposits associated with these fluvial features. The Mars Odyssey Thermal Emission Imaging System (THEMIS) infrared (IR) daytime mosaic (100 meters per pixel) was used as the primary base map. We constructed the geologic map of western Ladon basin at 1:1,000,000 scale. We identified 20 geologic units in the map area, which we divided into the following groups: crater units; volcanic units; chaotic units; basin fill units; crater, valley, and channel units; and plateau and highlands units. We observed the following sequence of events: (1) formation of Ladon and Holden impact basins in the middle Noachian, producing a mountainous unit; (2) landscape degradation and infilling during the Late Noachian to Late Hesperian forming terra units; (3) formation of Ladon Valles in the Late Noachian to Early Hesperian by catastrophic flooding, producing channel units; (4) accumulation of sediments in the Late Noachian to Early Amazonian, forming basin fill units; (5) formation of a smooth crater fill unit during the Early to Late Hesperian; (6) deposition during the Late Hesperian to Early Amazonian, creating light-toned layered units within Ladon Valles, Ladon basin, and other smaller valley networks along the western uplands, as well as formation of a light-toned unit along some crater floors; (7) formation a of chaotic unit in the Early to Middle Amazonian, formation of alluvial fans in the Early Amazonian, and eruption of a volcanic unit in the Middle Amazonian; and (8) formation of craters throughout the geologic history of the map region.

Scientific Investigations Map

Satellite interferometry landslide detection and preliminary tsunamigenic plausibility assessment in Prince William Sound, southcentral Alaska

Regional mapping of actively deforming landslides, including measurements of landslide velocity, is integral for hazard assessments in paraglacial environments. These inventories are also critical for describing the potential impacts that the warming effects of climate change have on slope instability in mountainous and cryospheric terrain. The objective of this study is to identify slow-moving landslides in the Prince William Sound region, southcentral Alaska, United States, which has had rapid deglaciation since the mid-1800s, and assess their tsunamigenic plausibility. We use an automated time series persistent scatterer interferometric synthetic aperture radar processing method with 7 years of Sentinel-1 data (2016–22) to identify 43 slow-moving slopes with average velocities ranging from approximately 0.2 to 21 millimeters per year. Landslide presence is confirmed using aerial imagery and previous landslide inventory records. We assess the tsunamigenic plausibility of the landslides using empirically derived estimates of landslide mobility based on modeled landslide volumes. Of the identified landslides, our preliminary analysis suggests that 11 have tsunamigenic potential if they were to fail rapidly and catastrophically. Although our estimate of tsunamigenic plausibility is preliminary and can be refined with additional observations and analyses, it can be used to prioritize ongoing and future hazard assessment, surveillance, and research efforts.

Alaska

Integrated stratigraphic and geochemical analysis of organic-rich intervals of the Lewis Shale in the eastern Washakie Basin, Wyoming

Geological studies in the Cretaceous Western Interior Seaway (KWIS) in North America reveal highly variable sedimentological conditions on its western shore caused by rapidly changing sea level and detrital input during the seaway’s closure. Here we examine a 601-foot (183 meters) continuous core through the lower part of the Maastrichtian Lewis Shale in the eastern Washakie Basin, Wyoming, through integrating stratigraphic and geochemical analyses to better understand organic matter deposition and preservation during the final marine transgression within the seaway. The core penetrates eight organic-rich flooding intervals (F1-F8) and a regional condensed section, informally known as the Asquith marker. The lower portions of the core record sediment input from a southern source, likely the Sierra Madre/Park Ranges, while the upper part records sediment input from a northern source, likely the Granite Mountains. This provenance transition is supported by regional stratigraphic analysis and changes in bulk mineral and trace metal composition. The Asquith marker and early flooding surfaces are enriched in oil-prone, marine organic matter deposited under dysoxic to euxinic conditions, whereas younger flooding surfaces show increased terrigenous input and poorer preservation conditions. The Asquith marker is identified as a prime oil-prone source rock. Although younger flooding surfaces also exhibit favorable source-rock properties, their generative potential is reduced due to increased clastic and terrigenous organic matter deposition caused by regression. These results emphasize the importance of sediment source variability, organic matter preservation, and changing redox conditions to provide insight into the sediment provenance and petroleum potential of the Lewis Shale in the context of the final transgression within the KWIS.

Colorado, Wyoming

Solute transport and modeling of water quality in a small stream

An injection of chloride, sodium, and stable strontium was made at a constant rate for 3 hours into Uvas Creek, Santa Clara County, Calif., to determine the mass transport processes in a small stream. Five observation points were selected within a 610-metre reach of the stream below the injection site. Water samples were collected at the observation points during and immediately after the injection. A mathematical model of the stream was obtained by solving analytically and optimally the one-dimensional mass transport equation of the solutes in the stream. Comparison of field results with a simplified mathematical model indicates the dominance of convection in the behavior of sodium and chloride. The concentration of chloride and sodium can be closely simulated by the model. However, strontium cannot be well represented by the simplified model, which contains a first-order decay-type sink.

California

Global performance of remote sensing-based and reanalysis-driven models to estimate open water evaporation

Evaporation plays an essential role in the water cycle, influencing local and regional climates while directly impacting water availability in lakes. However, directly measuring evaporation over water bodies remains challenging due to the high costs of installing and maintaining the required in situ instrumentation. Although several remote sensing algorithms have been providing evaporation estimates, the lack of a global validation hinders our understanding of their relative uncertainties and performances across different regions. Here, we analyze the performance of a suite of models that leverage satellite data and meteorological reanalysis to estimate evaporation over lakes worldwide. We compare 3 remote sensing-based models, 1 reanalysis-driven model and 1 ensemble approach, using in situ observations from 27 lakes representing a diverse range of geographic and climatic regions. Our results demonstrate that, overall, the ensemble outperformed any individual model in terms of accuracy, with a RMSE and a bias of 1.3 and 0.3 mm day −1 , respectively. These findings highlight the benefits of using an ensemble approach to estimate open water evaporation with satellite-based models at the global scale, leveraging the unique strengths of each model. For the individual models, differences in the representation of heat storage changes and advection effects led to lower values of RMSE and bias, depending on the location and depth of the lakes. This study sets the path for future improvement of open water evaporation algorithms globally, while remote sensing techniques are proven satisfactory to monitoring of water loss in lakes globally, an essential step toward effective large-scale water resources management.

Water Resources Research

The Sand AbrasioN Device for Aeolian Research (SANDAR): A new experimental device for investigating how wind transport affects sand on Earth and Mars

On Earth and Mars, aeolian transport causes sand grains to become abraded, resulting in mineralogic and textural changes. Understanding how sands evolve, or mature, with transport via experimental studies is important for understanding the origins, geologic history, and cycling of sediments, as well as dust production. Previous experimental works have used a variety of methods to simulate aeolian transport in the laboratory, but practical limitations and similitude concerns have limited such research. Here, we present and validate the Sand AbrasioN Device for Aeolian Research (SANDAR), a modified air mill that uses pressurized air to circulate sand around a small abrasion chamber, simulating the effects of aeolian transport. This device is re-circulating to simulate long-distance transport, and it allows for repeated analyses of well-constrained sediment samples, revealing their evolution over time. It is compatible with the grain sizes (74–500 μm) and grain impact velocities (∼0.6–3.7 m/s) typically expected for natural aeolian environments, and is also adaptable for diverse applications simulating different wind conditions. We show that the SANDAR achieves similitude of kinetic energy with respect to saltating sand on both Earth and Mars. SEM and optical microscope imaging reveal that the SANDAR produces microtextures on the surfaces of sand grains similar to those found with natural aeolian transport, demonstrating that it effectively simulates the mechanical effects of aeolian processes. Thus, the SANDAR is a valid tool for use in experimental research to improve our understanding of sedimentary processes across the Solar System.

Aeolian Research

Hydrogeologic conceptual model of groundwater occurrence and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River is a major concern in the southwestern United States where the river provides water to about 40 million people for municipal and industrial use and is used to irrigate about 5.5 million acres of land. Much of the salinity in the Colorado River Basin is derived from natural interactions of surface water and groundwater with various geologic materials (rocks, soils, and alluvial deposits). The Dolores River in southwest Colorado is a major tributary of the Colorado River that historically accounts for about 6 percent of the salinity load to the Upper Colorado River Basin with the Paradox Valley being the primary source of salinity to the Dolores River. The Paradox Valley, one of several salt-anticline valleys in the region, is a fault-bounded topographic basin aligned with and exposing an underlying salt-anticline core. Salt deposits in the Pennsylvanian Paradox Formation of the Hermosa Group form an elongated salt diapir oriented northwest to southeast that is up to 12,000 feet (ft) thick beneath the present valley floor. Surface erosion, groundwater circulation, and weathering during Tertiary and Quaternary valley formation contributed to development of a cap rock, collapse features, breccia, and brine at the top of the exposed salt diapir. Today (2023), brine occurring in the brecciated cap rock and underlying salt deposits is in hydraulic connection with an overlying freshwater alluvial aquifer, and depending on seasonal river stage and hydrologic conditions, the brine discharges to the Dolores River causing the observed increase in salinity as the river crosses the Paradox Valley. To reduce salinity concentrations in the Dolores River, the Bureau of Reclamation (Reclamation) operates the Paradox Valley Unit (PVU). The PVU project consists of nine shallow brine pumping wells near the Dolores River and one deep disposal well where the brine is injected for disposal. When operational, the PVU pumping wells extract brine from the base of the alluvial aquifer that is piped and injected into a deep disposal well about 3 miles southwest of the PVU. The PVU became fully operational July 1, 1996, and by 2015, operation of the PVU had reduced salinity concentrations in the Dolores River by as much as 70 percent compared to pre-PVU conditions. In response to a 4.5 magnitude earthquake, injection operations, and thus PVU pumping, were ceased from March 2019 to June 2022. A trial period of PVU operation began in June 2022 with a reduced injection rate, and thus PVU pumping rate, of about two-thirds capacity to gather additional information and guide future operational decisions. In cooperation with Reclamation, the U.S. Geological Survey (USGS) developed this report to present the current (2023) understanding of groundwater and brine occurrence and discharge to the Dolores River in the Paradox Valley. Results from the compilation of spatial datasets, groundwater sampling and age dating, and aquifer tests are presented to provide improved understanding of the Paradox Valley hydrogeology, to supply datasets for a numerical groundwater-flow and brine-transport model, and to support future operations of the PVU. The hydrogeologic data provided herein, along with the most recent loading analysis for the Dolores River in the Paradox Valley, and a previous conceptual model for brine discharge to the river are used to present a conceptual understanding of groundwater occurrence in the Paradox Valley.

Colorado

Geologic map of the Greater Antilles and Virgin Islands

Introduction The geologic map of the Greater Antilles and Virgin Islands is a compilation of information from the literature, integrated to provide a seamless geologic map of the region. This map was prepared to serve as a base map for a mineral resource assessment of the region. Several small-scale regional geologic maps of the region have been prepared in the past. This report supersedes an earlier version of the geologic map of the Greater Antilles and the Virgin Islands that was released as U.S. Geological Survey (USGS) Open-File Report 2019–1036. For this report, the regional geologic overview shown on sheet 1 covers Cuba; the island of Hispaniola, which includes Haiti and Dominican Republic; Jamaica; the Cayman Islands; Puerto Rico; and the U.S. and British Virgin Islands. Sheet 2 shows the geology of Cuba and Cayman Islands, sheet 3 the geology of Jamaica, sheet 4 the geology of Hispaniola, and sheet 5 the geology of Puerto Rico and Virgin Islands. Accompanying the maps are three tables: table 1 lists the map units (in alphabetical order by map-unit label), their ages, and the countries or territories in which they can be found. These are listed in alphabetical order for ease in locating unit names from the symbols on the maps. Table 2 is a list of sources for the geologic map compilation of Puerto Rico. Table 3 shows formations in the Devil’s Race Course Group of Jamaica.

Greater Antilles, Virgin Islands

New occurrences of the rare, REE minerals daqingshanite, törnebohmite, biraite, sahamalite, and ferriperbøeite from the Sheep Creek area, Montana, USA

Over 30 small, discontinuous, tabular carbonatite bodies are located in the Sheep Creek area, Ravalli County, southwest Montana. The age and origin of these REE-Nb-rich deposits are currently being investigated. The purpose of this paper is to document the occurrence of several rare minerals, including daqingshanite, törnebohmite, biraite, sahamalite, and ferriperbøeite, in two of the carbonatite bodies. These minerals are found in association with monazite, hydroxylbastnäsite, ferriallanite, calcite, dolomite, baryte, quartz, actinolite, apatite, celsian, and Sr-rich aragonite. Automated SEM-EDS was used to target the areas of interest in polished specimens for more detailed spot SEM-EDS and electron probe microanalysis. Raman spectra were also acquired for each of the rare minerals. The complex mineralogy of the Sheep Creek carbonatites is most likely due to several overlapping thermal events, including primary magmatic, overprinting hydrothermal, and supergene weathering stages. The rare minerals described in this study are believed to be hydrothermal and/or carbothermal in origin, although no estimates of temperature are available at this time.

Montana

Oligocene–Miocene development and evolution of the south Dome Rock Mountains basin, lower Colorado River corridor, Arizona, USA

Sedimentary basins in the Colorado River extensional corridor record large-magnitude Basin and Range extension and younger dextral shear deformation in the evolving Pacific−North America plate boundary. The south Dome Rock Mountains basin is located in west-central Arizona (USA), where the history of basin evolution, style of deformation, and timing of the transition between extension and dextral shear are not well constrained. We integrate new geologic mapping of the south Dome Rock Mountains basin with zircon U-Pb geochronology of six marker beds to characterize the timing of basin evolution and the slip history of the south Dome Rock Mountains normal fault. Structural analysis defines structures consistent with regional extension and younger dextral shear deformation. We use sedimentological and petrological analysis to interpret the depositional environments of three basin sequences. A lower basin sequence consists of fluvial strata deposited in an internally drained intermontane basin from ca. 35 Ma to 24.4 Ma. A sequence of volcanic rocks was emplaced in the basin between 24.4 Ma and 23.3 Ma, before the onset of local extensional faulting. An upper basin sequence of coarse conglomerate and sedimentary breccia was deposited on the hanging wall of a half graben structure after 23.3 Ma and through ca. 12.8 Ma, synchronous with 3−7 km of dip slip on the west-dipping, listric south Dome Rock Mountains fault. The basin was subsequently deformed by dextral shear after 12.8 Ma and before 4.8 Ma. These results document how distributed deformation related to the evolving Pacific−North America plate boundary occurred >100 km from the primary plate boundary.

Arizona, California