Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “NOAA Story Map”

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.

458 records · Page 8Linked to original sources
USGS2025

Alteration mapping in granitic gneiss using handheld geophysical and geochemical instruments: Implications for iron oxide-apatite and rare earth elements exploration

The Adirondack Mountains of New York, U.S.A. contain iron oxide-apatite (IOA) mineral deposits with variable concentrations of rare earth elements (REE). The IOA mineral deposits are typically hosted in the Lyon Mountain Granite Gneiss and are spatially correlated with extensive Na metasomatism (albitization) of the surrounding country rocks, although some mineral deposits also occur in metagabbro, paragneiss, and anorthosite hosts. The location of albitization zones is key to finding new IOA mineral deposits and better understanding the mineralization processes associated with their genesis. However, the Na alteration zones are generally not visibly identifiable in outcrop or hand specimens because the color and textural changes are subtle and are thus difficult to map. Here, we discuss the results from testing two spectroscopic techniques (handheld, field gamma-ray spectroscopy, and portable X-ray fluorescence (pXRF)) to characterize albitization zones surrounding IOA mineral deposits. The gamma-ray and pXRF spectrometers successfully distinguished the relative intensity of alteration in the Lyon Mountain Granite Gneiss based on the K abundance. The measured K content decreases towards the IOA mineral deposits, and the magnitude and width of the gradient are similar along the strike of each mineral deposit. Elevated Th and U values are present in host rocks adjacent to REE-bearing IOA mineral deposits. The pXRF and gamma-ray spectrometer K values are remarkably consistent with laboratory-based, whole-rock XRF compositional data and, therefore, useful for semi-quantitative analysis. Notably, albitization aureoles are consistent around REE-rich and REE-poor mineral deposits, suggesting that saline fluids are essential for the petrogenesis of IOA mineral deposits but may not be critical to REE mineralization. Ultimately, this study demonstrates the utility of handheld gamma-ray and pXRF spectrometry for identifying otherwise cryptic albitization gradients associated with IOA mineral deposits in granitic gneiss.

New York↗
USGS2025

Mapping potential sensitivity to hydrogeomorphic change in the UMRS riverscape

In 2020 the U.S. Geological Survey (USGS), as part of the Upper Mississippi River Restoration (UMRR) Program, began a new project to characterize potential hydrogeomorphic change associated with hydrogeomorphic units (HGUs) and their catenae (units linked by their association with sediment sources and flow origins). The goal of the project was to develop a geographic information system (GIS) database of HGUs for the Upper Mississippi River System (UMRS) available to both scientists and river managers working on UMRR studies and HREP planning and design studies. The characterization was based on a hydrogeomorphic change hierarchical classification developed previously for the UMRS. The products were generated with automated techniques in a GIS using systemic datasets. Landforms were mapped from the 2015 UMRS topobathymetric dataset with geomorphon (shape-based) tools tailored for the large riverscape, valley bottom environments in the UMRS. A clustering analysis was applied to the resulting landforms to identify HGUs associated with zones of perennial low flows, bankfull flows, and overbank floodplains. Catenae were assembled based on the proximity of the units to the main channel, tributary mouths, and side channels from previously published aquatic areas (USACE, 2018) coupled with least-cost flowpath linkages between potential sediment origins and planform change units developed by Rogala, Fitzpatrick, and Henderson (2020). These GIS-based analyses were successful at identifying a range of HGUs using an automated technique with available data across the entire riverscape, with emphasis on those that have the potential for hydrogeomorphic change. Most of the resulting features are depositional, as expected in a large river system. However, this is the first attempt of linking tributary inputs, side channel erosion and levee breaches with their depositional counterparts. The approach was successfully piloted in Pools 8 and 10 in the Upper Impounded Reach and Pool 14 in the Lower Impounded Reach, with next steps for application in reaches of the unimpounded section and Illinois River. This report emphasizes results from Pool 10, which was the focus of most of our attention during the pilot phase.

Illinois, Indiana, Iowa, Minnesota, Missouri, Sout↗
USGS2026

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2024

Potentiometric-surface and depth-to-water maps for spring 2024 were created for the Mississippi River Valley alluvial aquifer using groundwater-altitude data from 1,151 wells completed in the Mississippi River Valley alluvial aquifer and from the altitude of the top of the water surface in area rivers from 160 U.S. Geological Survey and U.S. Army Corps of Engineers streamgages. The potentiometric-surface and depth-to-water maps for 2024 were created to support investigations to characterize the Mississippi River Valley alluvial aquifer as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the Mississippi River Valley alluvial aquifer for spring 2024 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 330 feet (ft) above the North American Vertical Datum of 1988. The regional direction of the groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient direction was into the depression, and near rivers, where the groundwater gradient direction was either from the aquifer to the river or from the river into the aquifer. Groundwater depressions in the potentiometric-surface map are in the lower one-half of the Cache region and in most of the Grand Prairie and Delta regions. Depth to water by well in the Mississippi River Valley alluvial aquifer in spring 2024 ranged from 0.27 ft above land surface to 145.65 ft below land surface.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗
USGS2025

Approximate inland extent of saltwater intrusion at the base of the Biscayne aquifer, Miami-Dade County, Florida, 2022

Miami-Dade County is part of a densely populated urban corridor in southeastern Florida. The Biscayne aquifer serves as Miami-Dade County’s primary drinking water source and is characterized by highly permeable karstic limestone and carbonate sand. The aquifer’s coastal location and permeable nature make it susceptible to saltwater intrusion. Monitoring the current inland extent and the rate of movement of the saltwater front in the aquifer can inform management strategies for conserving the long-term sustainability of the county’s water supply. In the 1950s, the U.S. Geological Survey published a map of the inland extent of saltwater intrusion in the Biscayne aquifer and has continued to update this map to monitor changes over time, with the most recent update published in 2018. An updated map has been created showing the approximate inland extent of saltwater intrusion in the Biscayne aquifer in eastern Miami-Dade County in 2022, with the 2018 extent shown for comparison. The inland extent of saltwater intrusion was mapped through the interpretation of borehole electromagnetic induction logs and measurements of chloride and specific conductance in groundwater samples. The location of the saltwater interface at the base of the Biscayne aquifer was represented by the 1,000-milligram-per-liter isochlor. This report describes changes in the location of the saltwater interface from 2018 to 2022. By 2022, the saltwater interface had moved farther inland in both the northern and southern parts of the county, advancing by as much as 0.3 kilometer in the north and up to 0.8 kilometer in the Model Land Area to the south. However, it remained relatively unchanged from its 2018 position in the east-central part of the county.

Florida↗
USGS2025

Methodology for compilation of previously published contour data showing the altitude of the base of Dakota Sandstone on the Colorado Plateau, Arizona, Colorado, New Mexico, and Utah

Structure contours and other geologic information from numerous published geologic maps were digitized and compiled into a digital dataset showing the configuration of a single stratigraphic datum, the base of the Dakota Sandstone and its equivalents across the Colorado Plateau. The principal maps compiled in digital form are a series of 1:250,000-scale 1 degree (°) × 2° quadrangle maps published by the U.S. Geological Survey, augmented by other geologic maps published at various map scales. The compiled digital dataset contains geologic map polygons of the Dakota Sandstone and regional stratigraphic equivalents, the location of faults and fold axes, structure contour lines that define the altitude of the base of the unit and bedding orientation data computed from the structure contour lines. This report provides the scientific rationale for compilation of these data and describes the compilation methodology for each of the data elements. This report provides an extended description of the data compilation in a companion U.S. Geological Survey digital data release of spatial data and attributes associated with the contoured surface and associated geologic data layers.

Arizona, Colorado, New Mexico, Utah↗
USGS2024

Field observations and logs from the Rose Hip trench exposure across a north-facing scarp within the Seattle Fault Zone, southern Bainbridge Island, Washington

The Seattle Fault Zone is an approximately 70-km-long, east-west-trending zone of south-dipping blind reverse faults within the Puget lowland region in Washington. Because of the proximity, the Seattle Fault Zone poses a significant earthquake hazard to the Puget sound and Seattle metropolitan regions. We present preliminary mapping and trench-site information from a paleoseismic investigation across a newly identified active fault scarp located within the hanging wall of the Seattle Fault Zone on southern Bainbridge Island, Washington. The trench exposed monoclinally folded Miocene bedrock, fractured and faulted glacial-related deposits, and laminated lacustrine deposits capped by slope-derived colluvium. The observations from this investigation record late Pleistocene to Holocene north-vergent folding and faulting along this new fault scarp.

Washington↗
USGS2026

Final project peport for “Mapping riparian vegetation response to climate change on the San Carlos Apache Reservation and Upper Gila River watershed to inform restoration priorities: 1935 to present – Phase 2: Focus on tamarisk vegetation”

Riparian ecosystems play a critical role in supporting wildlife habitats, maintaining water quality, and sustaining ecological resilience in arid landscapes. In the Upper Gila River level-4 Hydrologic Unit Code (HUC-4; Identification Number – 1504) watershed of Arizona and New Mexico, and across areas of the San Carlos Apache Tribe of the San Carlos Apache Reservation (hereafter, Tribe/Tribal – entity; or Reservation - place), riparian ecosystems have been substantially altered by the widespread expansion of tamarisk ( Tamarix spp.), an invasive and non-native riparian species. Tamarisk has high water use and increased flammability, making it a growing concern as droughts intensify across the southwestern United States (U.S.). Recent research indicates that tamarisk is increasingly stressed under prolonged drought conditions, which can elevate wildfire risk and further degrade riparian habitat. In response, the U.S. Geological Survey (USGS) Western Geographic Science Center, in collaboration with the Tribe, developed remote sensing–based tools to map riparian vegetation composition and monitor vegetation condition over time. These tools enable accurate identification of tamarisk extent, detection of vegetation stress, and comparison with native species such as cottonwood and willow. This information directly supports restoration and management actions, including targeted tamarisk removal, protection of endangered species habitat, and prioritization of areas for native vegetation recovery. The tools also provide insight into how riparian vegetation responds to changing climate and hydrologic conditions, strengthening long-term planning for riparian forest management. Findings from this work demonstrate that riparian vegetation responses vary across river systems, indicating differences in plant composition and hydro-climatic conditions. Tamarisk on both the Gila and San Carlos Rivers generally exhibit greater declines in greenness in response to higher temperature as well as lower precipitation and river flow, while cottonwood ( Populus fremontii ) and willow (Salix spp.) are particularly sensitive to changes in discharge along the San Carlos River and respond more directly to both high and low flows, with more moderate responses to precipitation and temperature. By providing actionable, science-based information, this research supports climate adaptation planning, wildfire risk reduction, and improved riparian ecosystem health, with benefits that extend beyond the Reservation to regional water and wildlife conservation efforts.

Arizona, New Mexico↗
USGS2023

Sediment thickness map of United States Atlantic and Gulf Coastal Plain Strata, and their influence on earthquake ground motions

With the recent successful accounting of basin depth ground-motion adjustments in seismic hazard analyses for select areas of the western United States, we move toward implementing similar adjustments in the Atlantic and Gulf Coastal Plains by constructing a sediment thickness model and evaluating multiple relevant site amplification models for central and eastern United States seismic hazard analyses. We digitize and combine existing sediment thickness data sets into a composite surface that delineates the base of Cretaceous sediments under the Atlantic Coastal Plain and the base of Mesozoic sediments under the Gulf Coastal Plain. Amplification models dependent on sediment thickness, site natural period, and source-to-site path length are compared with data sets of observed ground motions to evaluate the ability of the new models to improve ground motion estimates. We find that the amplification models can account for observed trends in sediment-thickness and period-dependent residuals, but some tuning is required. For example, the model of Chapman and Guo requires a reference V S 30 , the time-averaged shear-wave velocity within 30 m of the Earth’s surface, for non-Coastal Plain sites, which we estimate to be between about 1 and 2 km/s. Along with our sediment thickness model, we estimate a velocity profile for application to the Harmon et al. site-natural-period-based model in order to best match the Chapman and Guo period dependence for a broad range of sediment thicknesses. The Next Generation of Attenuation models for the eastern United States Gulf Coast path-based adjustment models can also account for seismic attenuation in the Coastal Plain sediments and reduce the standard deviation of total residuals. If enacted in the U.S. Geological Survey National Seismic Hazard Model, these amplification models will reduce predicted short-period (<1 s) and increase predicted long-period (>1 s) ground motions in the Coastal Plains appreciably.

Earthquake Spectra↗
USGS2024

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

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

Fire↗
USGS2024

Geologic framework and hydrostratigraphy of the Edwards and Trinity aquifers within parts of Bandera and Kendall Counties, Texas

The karstic Edwards and Trinity aquifers are classified as major sources of water in south-central Texas by the Texas Water Development Board. During 2019–23 the U.S. Geological Survey, in cooperation with the Edwards Aquifer Authority, mapped and described the geology and hydrostratigraphy of the rocks composing the Edwards and Trinity aquifers within parts of Bandera and Kendall Counties from field observations of the surficial expressions of the rocks. The thicknesses of the mapped lithostratigraphic and hydrostratigraphic units were also estimated from field observations in the study area. The Cretaceous rocks in the study area are part of the Trinity Group and Edwards Group. The groups, formations, and members are composed primarily of layers of marls, shales, and limestones. The limestones are composed of mudstone through grainstone, framestone and boundstone, dolomite, and argillaceous and evaporitic rocks. The principal structural feature in the study area is the Balcones fault zone. The Balcones fault zone is the result of late Oligocene and early Miocene extensional faulting and fracturing that was a result of the eastern Edwards Plateau uplift. In the Balcones fault zone, most of the faults in the study area are high-angle to vertical, en echelon, normal faults that are predominantly downthrown to the southeast. Hydrostratigraphically, the rocks exposed in the study area are those that contain the Edwards aquifer, the upper zone of the Trinity aquifer, and the middle zone of the Trinity aquifer. Descriptions of the hydrostratigraphic units, thicknesses, hydrologic function, porosity types, and field identification and observations are provided, including those for the informal Bandera and Love Creek hydrostratigraphic units of the Edwards aquifer, which were identified through the mapping for this study.

Texas↗
USGS2024

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↗
USGS2025

Cenozoic stratigraphy of Colorado

As a successor to previous Colorado stratigraphy charts (MS-53 Colorado Stratigraphic Chart and MS-54 Cretaceous Stratigraphy of Colorado), this Colorado Geological Survey (CGS) publication resulted from a collaboration between the CGS, USGS, and the Denver Museum of Nature and Science (DMNS). The chart was designed to illustrate Cenozoic stratigraphy spanning the state’s many sedimentary basins. It builds upon the work of dozens of colleagues and updates Richard Pearl’s seminal 1974 stratigraphy chart. The chart leverages the community’s stratigraphic work in both the subsurface and outcrop, and depicts new geochronologic constraints for many units. To facilitate comparison of strata to external forcing factors, the chart employs a linear timescale. Each unit’s dominant depositional environment is depicted, as are major mountain building events, erosional events, and regional unconformities.

Colorado↗
USGS2025

Footprints of past mining in Alaska (USA) derived from high-resolution satellite imagery

Mapping the land area used for mining in the past is essential for guiding the remediation of affected landscapes and assessing the resource potential of related waste products. Despite significant recent progress delineating footprints of active and inactive mining globally, the known inventory of such mine lands remains incomplete. Here, I describe a new map dataset of footprints of land surface disturbance and waste at sites of past mining in Alaska (USA) based on visual interpretation of satellite imagery. This dataset maps 6–14 times the area of previous regional and global mine footprint maps in Alaska and is the first in the region to explicitly delineate mine waste landforms (e.g., tailings piles). The data are publicly available from the U.S. Geological Survey under a “no rights reserved” Creative Commons (CC0) license agreement.

Alaska↗
USGS2025

Potentiometric surface maps and groundwater-level hydrographs for confined aquifers of the New Jersey Coastal Plain, 2018

The U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection (NJDEP), prepared potentiometric surface maps for 10 confined aquifers of the New Jersey Coastal Plain physiographic province based on water-level measurements collected during late 2018 and early 2019 from 951 wells in New Jersey and parts of Pennsylvania and Delaware. Maps were prepared for the confined Cohansey aquifer, Rio Grande water-bearing zone, Atlantic City 800-foot sand, Piney Point aquifer, Vincentown aquifer, Wenonah-Mount Laurel aquifer, Englishtown aquifer system, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system. Potentiometric surface maps indicate regional cones of depression in the following aquifers and the counties in which they are centered: Atlantic City 800-foot sand in Atlantic County, the Piney Point aquifer in Cumberland County, the Wenonah-Mount Laurel aquifer and Englishtown aquifer system in Monmouth and Ocean Counties, the Wenonah-Mount Laurel aquifer in Camden and Gloucester County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Ocean County, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system in Camden County. Cones of depression with smaller areal extents were in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Atlantic City 800-foot sand centered in Cape May County, the Piney Point aquifer centered in Ocean County, the Wenonah-Mount Laurel aquifer in Salem and Burlington Counties, the Englishtown aquifer system in Camden County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth County, and the middle aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth, Ocean, and Salem Counties. No cone of depression was interpreted in the Vincentown aquifer. Long-term hydrographs are presented for 75 wells spanning each of the 10 confined aquifers, and contain a mix of discrete water-level measurements and daily mean water levels based on continuously recorded 15-minute data. Changes of water levels during 2014–19, as indicated by the hydrographs, were compared with those of previous periods to assess any departures from historical data. During 2014–19, water levels were stable and fluctuated within similar ranges as previous periods in the following aquifers and locations: all wells in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Vincentown aquifer, the Englishtown aquifer system, the Piney Point aquifer wells in Burlington and Ocean Counties, six of eight wells in the Wenonah-Mount Laurel aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system except those within NJDEP Critical Area II. Increasing water levels in 2014–19, ongoing since historical periods, were indicated in the following aquifers and locations: Atlantic County wells in the Piney Point aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system within NJDEP Critical Area II. Water levels in the Atlantic City 800-foot sand also increased during 2014–19 in wells in Atlantic County and northern Cape May County closer to the center of the cone of depression in that aquifer, which is a response unique to this period and absent from previous periods. During 2013–19, continued decreasing water levels, ongoing since previous periods, were indicated by hydrographs of Atlantic City 800-foot sand wells in southern Cape May County, Piney Point aquifer wells in Cumberland County where the regional cone of depression is located, and two wells in the Wenonah-Mount Laurel aquifer—070478, which in 2014–19 departed from previous periods, and 330020, which continued a gradual decrease throughout its period of record.

New Jersey↗
USGS2025

Lunar grid systems, coordinate systems, and map projections for the Artemis missions and lunar surface navigation

Foreward This document contains design specifications of a navigational standard for the Moon, including a Lunar Transverse Mercator system, a Lunar Polar Stereographic system, a Lunar Grid Reference System, and a unique coordinate structure, Artemis Condensed Coordinates, for Artemis mission navigation and lunar surface science. The National Aeronautics and Space Administration (NASA) Artemis campaign seeks to place humans on the Moon for the first time since the Apollo missions. Early Artemis missions are heavily focused on the lunar south pole, which promises to return valuable data on the Moon’s geologic record, amongst other mission objectives. Coordinate systems in use today for the lunar south pole provides crew members on the surface neither an efficient nor intuitive means to communicate their position and orientation. A novel grid coordinate system, the Lunar Grid Reference System, is proposed to address these concerns for use in real-time extravehicular activity operations on the lunar surface. The many stakeholders involved in the Artemis missions will need a common system to communicate position and orientation while astronauts are operating on the lunar surface. To that end, Artemis crew members will need that system to be efficient and intuitive to promote efficient extravehicular activity timelines and reduce confusion. In the context of this document, these characteristics are addressed on the design of lunar coordinate systems: Efficient.—The number of characters required to communicate a location within a desired precision level in both local and global contexts, and how many steps are required for a recipient or sender to interpret a location. Intuitive.—How well the system aligns with human perceptual abilities, and whether the system yields distances that have the same relationship to actual lunar surface distance in all directions from the point where a person is located. Technological systems are currently being investigated to supplement the crew members’ ability to locate and orient themselves and other assets on the lunar surface; however, it is unlikely that those systems will be fully operational for the first few landed missions. Even with future positional aids, crew members will still need an efficient and intuitive means to communicate position and orientation. In addition, if technological systems fail, the crews will require land navigation skills and have maps available, thus providing further motivation for a crew-centric coordinate system. The contents of this U.S. Geological Survey (USGS) document detail a comprehensive framework for standardizing lunar crewed surface navigation within NASA and outlines the protocols, methods, and designs necessary for achieving consistency and interoperability across relevant space mission teams and lunar surface navigators. Key components of this document include designs of map projections, projected coordinate reference systems (Lunar Transverse Mercator and Lunar Polar Stereographic systems), and a grid system (Lunar Grid Reference System and Artemis Condensed Coordinates) for the Moon. The work proposed in this document seeks to accomplish something similar to the National Geospatial-Intelligence Agency (NGA) document SIG 0012 (NGA, 2014a), but for using grid systems for the Moon. This report incorporates initial feedback and input from NASA’s Artemis Geospatial Data Team, NASA’s Flight Operations Directorate, National Geodetic Survey, USGS Astrogeology Science Center, and NGA and is intended to serve as a resource for all involved with the Artemis missions, as well as for engineers designing and operating lunar infrastructure.

Techniques and Methods↗
USGS2026

Geologic framework, hydrostratigraphy, and ichnology of the Blanco, Payton, and Rough Hollow 7.5-minute quadrangles, Blanco, Comal, Hays, and Kendall Counties, Texas

During 2023, the U.S. Geological Survey, in cooperation with the Edwards Aquifer Authority, revised an initial characterization completed during 2015–16 of the geologic framework, hydrostratigraphy, and ichnology of the Edwards and Trinity aquifers in the Blanco, Payton, and Rough Hollow 7.5-minute quadrangles in Blanco, Comal, Hays, and Kendall Counties, Texas. The purpose of this report is to present the updated geologic framework, hydrostratigraphy, and ichnology of the Trinity and Edwards Groups in those quadrangles. Rocks exposed in the study area are of the Lower Cretaceous Trinity Group and lower part of the Fort Terrett Formation of the Lower Cretaceous Edwards Group. The faulting and fracturing in the study area are part of the Balcones Fault Zone, an extensional system of faults active during the Paleocene to middle Eocene that generally trends southwest to northeast in south-central Texas. Hydrostratigraphically, the rocks exposed in the study area contain a section of the Edwards (Balcones Fault Zone) aquifer, upper zone of the Trinity aquifer, and middle zone of the Trinity aquifer. In the study area, the only hydrostratigraphic units of the Edwards (Balcones Fault Zone) aquifer remaining are (from top to bottom) VII and VIII, which cap several hills. The mapped hydrostratigraphic units of the upper zone of the Trinity aquifer are (from top to bottom) the cavernous, Camp Bullis, upper evaporite, fossiliferous, and lower evaporite. The mapped hydrostratigraphic units of the middle zone of the Trinity aquifer are (from top to bottom) the Bulverde, Little Blanco, Twin Sisters, Doeppenschmidt, Herff Falls (where present), Rust, Honey Creek, Hensell, and Cow Creek hydrostratigraphic units.

Texas↗
USGS2026

Lunar analog study using portable gamma-ray neutron detector: Radiochemical mapping of silicic-to-basaltic volcanic terrains in the San Francisco Volcanic Field, Arizona

Studying lunar silicic volcanism provides key insights into the Moon’s crustal evolution, magmatic processes, and volcanic history. The Lunar-VISE (Lunar Vulkan Imaging and Spectroscopy Explorer) mission will investigate the Gruithuisen domes, a unique lunar region hypothesized to have formed through silicic volcanism. Using instruments on a Firefly Aerospace lander and a Honeybee Robotics rover, Lunar-VISE will analyze mineralogy, geochemistry, and surface properties to determine the origin and evolution of the domes, with a gamma-ray and neutron spectrometer (LV-GRNS) among its payload instruments. In preparation for this mission, we conducted preliminary fieldwork using a handheld gamma-ray neutron detector with NaI(Tl) and Cs LiYCl :Ce (CLYC) scintillators. This study focuses on various rhyolitic and basaltic volcanic centers in the San Francisco Volcanic Field (SFVF) near Flagstaff, Arizona—specifically Sugarloaf Peak, Bonito Lava Flow, and Robinson Mountain, a region containing several well-characterized lunar analog sites. The SFVF was selected for this study due to its broad compositional diversity spanning basaltic to rhyolitic compositions, its well-preserved volcanic morphologies, and its extensive use in previous NASA field campaigns and astronaut training exercises, making it an ideal terrestrial laboratory for testing planetary exploration techniques. We measured natural radioactivity, specifically from potassium (K), thorium (Th), and uranium (U) and other elements in their decay chains, which serve as diagnostic tracers of magmatic differentiation and crustal evolution processes, to assess geochemical variability across compositionally diverse terrains. The detector’s sensitivity was assessed across varying concentration levels. Regionally averaged concentrations of radioisotopes were determined by gamma-ray spectroscopy at selected sites. Our measurements reveal significant compositional variations between sites, with Sugarloaf Peak (rhyolitic, silicic dome) exhibiting the highest average radioisotope concentrations (K: 3.29 ± 0.24 wt%, U: 14.35 ± 1.81 ppm, Th: 27.14 ± 2.43 ppm), while Bonito Lava Flow (K: 1.13 ± 0.08 wt%, U: 3.86 ± 0.52 ppm, Th: 7.52 ± 1.02 ppm), and Robinson Mountain (K: 1.46 ± 0.16 wt%, U: 5.59 ± 1.06 ppm, Th: 11.75 ± 1.98 ppm) show lower concentrations aligned with basaltic compositions. Gamma-ray fluxes were elevated by a factor of approximately three to four at Sugarloaf Peak relative to nearby basaltic terrains, consistent with expected geochemical differentiation patterns. Furthermore, at meter scales, proximity to geological features significantly affects measurements. Th concentrations adjacent to a cliff face at Sugarloaf’s base were 28% higher than values measured 3-5 meters away from the same feature, demonstrating localized compositional heterogeneity. Analysis of station-by-station measurements reveals that K, U, and Th concentrations show an elevation trend at Robinson Mountain and generally higher concentrations near Sugarloaf Peak’s summit, suggesting progressive magmatic differentiation and/or the presence of more evolved lithologies at higher elevations. By mapping these elements on Earth using GRNS, we aim to optimize measurement techniques for Lunar-VISE and similar rover-borne missions by demonstrating the applicability of portable GRNS instruments for planetary surface exploration in an analog environment.

Arizona↗
USGS2026

Methodology for construction of a three-layer geologic model of the conterminous United States using land surface, top of bedrock, and top of basement

This report describes the methodology used for the construction of a digital three-layer geologic model of the conterminous United States by mapping the altitude of three surfaces: land surface, the top of bedrock, and the top of basement. These surfaces are mapped through the compilation and synthesis of published stratigraphic horizons from numerous topical studies. The mapped surfaces create a three-layer geologic model with three geomaterial-based subdivisions: unconsolidated to weakly consolidated sediment; layered consolidated rock strata that constitute bedrock; and crystalline rocks that are described as “basement,” consisting of either igneous, metamorphic, or highly deformed rocks. The data compilation and synthesis are highly dependent on the definition of the informal terms “bedrock” and “basement,” which may describe different ages or types of rock in different parts of the conterminous United States. This report presents the conceptualization of the three mapped layers, describes the datasets used, and summarizes the decisions made while compiling the three-layer model from the various sources. This digital dataset was created as part of efforts by the U.S. Geological Survey to develop subsurface geologic data in geospatial form as part of a broad directive to develop two-dimensional and three-dimensional geologic information at detailed, national, and continental scales. This digital dataset partly fulfills the goal of the U.S. Geological Survey’s National Cooperative Geologic Mapping Program to construct a national-scale three-dimensional geologic model.

conterminous United States↗