SEARCH · USGS Science
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.
This CD contains imagery collected by the Advanced Very High Resolution Radiometer (AVHRR) on the NOAA polar-orbiting weather satellites. The AVHRR provides almost daily coverage of a site at a resolution of approximately 1 km. The data types included are sea surface temperature (SST), water reflectance (REF), and a false-color infrared overview (FCI).
Microbial processing of atmospheric nitrogen (N) deposition regulates the retention and mobilization of N in soils, with important implications for water quality. Understanding the links between N deposition, microbial communities, N transformations, and water quality is critical as N deposition shifts toward reduced N and remains persistently high in many regions. Here, we investigated these connections along an elevation transect in the Colorado Front Range. Although rates of N deposition and pools of extractable N increased down the elevation transect, soil microbial communities and N transformation rates did not follow clear elevational patterns. The subalpine microbial community was distinct, corresponding to a high C:N ratio and low pH, while the microbial communities at the lower elevation sites were all very similar. Net nitrification, mineralization, and nitrification potential rates were highest at the Plains (1,700 m) and Montane (2,527 m) sites, suggesting that these ecosystems mobilize N. In contrast, the net immobilization of N observed at the Foothills (1,978 m) and Subalpine (3,015 m) sites suggests that these ecosystems retain N deposition. The contrast in N transformation rates between the plains and foothills, both of which receive elevated N deposition, may be due to spatial heterogeneity not captured in this study and warrants further investigation. Stream N concentrations from the subalpine to the foothills were consistently low, indicating that these soils are currently able to process and retain N deposition, but this may be disrupted if drought, wildfire, or land-use change alter the ability of the soils to retain N.
Excessive nitrogen export from agricultural watersheds remains a critical water quality challenge, with the Upper Mississippi River Basin (UMRB) significantly contributing to downstream eutrophication and hypoxia in the Gulf. This study investigates the spatiotemporal dynamics of riverine nitrate plus nitrite (NO 3 – + NO 2 – -N) export across the UMRB at high spatial resolution (12-digit Hydrologic Unit Codes or HUC12 subwatershed scale) during 2001–2020 and quantifies the effects of anthropogenic activities and hydrological variability on riverine NO 3 – + NO 2 – -N export changes in the region between 2001–2005 and 2016–2020. Our results revealed hotspots of substantial increases in NO 3 – + NO 2 – -N yields across the UMRB, with distinct regional patterns in driving factors. Over the entire UMRB, NO 3 – + NO 2 – -N yields increased by 9.7 kg/ha/yr on average from 2001–2005 to 2016–2020, with anthropogenic activities contributing 4.8 kg/ha/yr and hydrological variability contributing 4.9 kg/ha/yr. The northern and western UMRB had combined influences from both anthropogenic activities and hydrological variability, while the east-central regions had predominantly hydrologically driven changes. Agricultural sources, including fertilizer, manure, and biological nitrogen fixation, collectively contributed over 80% of NO 3 – + NO 2 – -N loading throughout the basin. This framework for disentangling human and hydrological impacts provides critical insights for developing effective and targeted watershed management strategies to reduce nutrient losses and improve water quality.
Multiple generations of spring-fed streams traversed ∼800 km 2 of the Las Vegas Valley in southern Nevada between ca. 10.9 ka and 8.5 ka, depositing an extensive tufa network. The scale of this network and diversity of tufa morphologies is novel in North America and offers an opportunity to obtain quantitative paleoclimate data for the region during the early Holocene. We determined isotopic compositions and estimated past temperatures using clumped isotope data from early Holocene tufa on the valley floor (698 m) as well as tufa forming today at higher elevation in the nearby Spring Mountains at Cold Creek Spring (1856 m). Modern and fossil tufa yielded comparably low δ 18 O values, implying that source waters for both were derived from high-elevation winter precipitation. Clumped isotope temperatures of modern tufa average 15.8 ± 2.5 °C, aligning with mean summer temperatures of the emergent spring water, and indicate equilibrium conditions of tufa formation. The early Holocene tufa yielded similar clumped isotope temperatures, averaging 15.2 ± 3.9 °C, meaning it precipitated at temperatures that occur at much higher elevations today. The Las Vegas tufa record, combined with nearby and temporally correlative paleospring and lacustrine records, suggest that cool/wet conditions prevailed throughout the Mojave Desert during the early Holocene. These records also demonstrate that spring ecosystems responded to millennial-scale hydroclimate variations that supersede climate change driven solely by insolation. The previously unrecognized pattern of ecosystem response to hydroclimate documented here may assist in understanding climate drivers for the early Holocene and provide critical information for the fate of groundwater-dependent ecosystems in the southwestern United States.
As of October 2024, the U.S. Geological Survey (USGS) operated 588 sites across the United States and its territories as part of the Groundwater Climate Response Network (CRN). The CRN is comprised of wells selected to monitor the effects of climate variability, such as droughts, on groundwater levels nationwide. The CRN includes nearly 500 locations with real-time data and more than 100 sites with non-real-time data available to the public on the CRN web mapper and the USGS National Water Dashboard.
A marine controlled source electromagnetic (CSEM) campaign was carried out in the Gulf of Mexico to further develop marine electromagnetic techniques in order to aid the detection and mapping of gas hydrate deposits. Marine CSEM methods are used to obtain an electrical resistivity structure of the subsurface which can indicate the type of substance filling the pore space, such as gas hydrates which are more resistive. Results from the Walker Ridge 313 study (WR 313) are presented in this paper and compared with the Gulf of Mexico Gas Hydrate Joint Industry Project II (JIP2) logging while drilling (LWD) results and available seismic data. The hydrate, known to exist within sheeted sand deposits, is mapped as a resistive region in the two dimensional (2D) CSEM inversion models. This is consistent with the JIP2 LWD resistivity results. CSEM inversions that use seismic horizons provide more realistic results compared to the unconstrained inversions by providing sharp boundaries and architectural control on the location of the resistive and conductive regions in the CSEM model. The seismic horizons include: 1) the base of the gas hydrate stability zone (BGHSZ), 2) the top of salt, and 3) the top and bottom of a fine grained marine mud interval with near vertical hydrate filled fractures, to constrain the CSEM inversion model. The top of salt provides improved location for brines, water saturated salt, and resistive salt. Inversions of the CSEM data map the occurrence of a ‘halo’ of conductive brines above salt. The use of the BGHSZ as a constraint on the inversion helps distinguish between free gas and gas hydrate as well as gas hydrate and water saturated sediments.
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.
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.
Estimates of dynamic groundwater volumes supplying baseflow to streams are important for water availability projections during extended periods of drought. The primary goals of this study were to provide dynamic storage volume estimates, inferred from streamflow recession analysis, for baseflow regimes within seven gaged catchments within the Neversink Reservoir Watershed (NRW), a critical municipal water source for New York City. Additionally, geomorphological properties, surficial geology and hydro-meteorological processes were quantified and described in relation to time and spatially variable recession behaviour and storage estimates across the NRW. To explore these relationships, we (1) evaluated seasonal trends in streamflow recession behaviour in relation to modelled potential evapotranspiration (PET) and catchment runoff rates, (2) derived empirical streamflow models for cool-season runoff using both linear and nonlinear reservoir assumptions for baseflow and (3) calculated metrics related to the geology and geomorphology of each catchment and compared these metrics to area normalised baseflow dynamic storage estimates. Results show that baseflow recession behaves as a nonlinear reservoir, and applying linear groundwater reservoir assumptions may underestimate the total dynamic storage volumes compared to what would be predicted for a nonlinear reservoir. Increases in PET caused decreases in storage conditions that resulted in increased recession rates and nonlinearity in streamflow recession during the growing season. Additionally, we found that while no single physical catchment characteristic solely predicted catchment storage dynamics, sediment volume and stream gradients were stronger predictors of normalised storage volumes than catchment surface area or surface topography alone. Within the NRW, catchments with the highest sediment volume exhibited the lowest recession rates and higher dynamic storage volumes, while the smallest catchment, mostly devoid of sediment, had the fastest recession rate and lowest dynamic storage volume.
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.
As of October 2023, the U.S. Geological Survey (USGS) operated more than 660 sites across the United States and its territories as part of the Groundwater Climate Response Network (CRN). The CRN is comprised of wells and springs selected to monitor the effects of climate variability, such as droughts, on groundwater levels and spring discharge nationwide. The CRN includes more than 550 locations with realtime data and more than 100 sites with non-real-time data available to the public on the CRN web mapper and the USGS National Water Dashboard.
Pervasive chemical weathering on stable cratons may form thick regoliths and elemental enrichment, but constraining the age of regolith formation is challenging. In this study we utilize multiple geochronological techniques on different minerals from the world-class Mount Weld rare earth element (REE) deposit, formed by lateritic weathering of a carbonatite, to constrain the age of formation and provide insight into landscape evolution. The oldest dates, ca. 100 to 38 Ma, are from Lu-Hf dating of churchite [HREE(PO 4 )·2(H 2 O)], a heavy REE phosphate mineral. Growth bands on individual minerals show a younging outwards. 40 Ar/ 39 Ar geochronology of cryptomelane [K(Mn 4+ ,Mn 2+ )₈O₁₆] yielded dates from ca. 40 to 27 Ma. Similarly, (U-Th)/He geochronology of goethite [FeO(OH)] yielded dates ranging from ca. 45 to 19 Ma. Integrating results into regional constraints, suggests 1) churchite formed by mineral saturation in a karst-like setting below the water table from ca. 100 to 40 Ma, 2) with minor uplift and erosion, cryptomelane and goethite formed at or near the water table between ca. 45 and 19 Ma, 3) after ca. 15 to 10 Ma chemical weathering within the profile had ended. Other studies document that the region experienced minimal uplift and a wet, warm climate from ca. 100 Ma to 15 Ma. These conditions and the high carbonate content of the carbonatite promote extensive chemical weathering, a deep weathering profile, and the preservation of the weathered section. This study highlights the use of multiple geochronological techniques utilizing different minerals to provide insight into how laterites form and to constrain the timing and history of the formation of this important mineral deposit.
Alluvial wetland ecosystems are vital as biodiversity hotspots but are increasingly threatened by anthropogenic stressors and drought. These pressures are especially acute in arid and semi-arid regions, where eco-hydrologic connectivity is fragile and recovery is slow. This study quantifies the efficacy of nature-based solutions, particularly the ‘Zeedyk approach,’ which employs low-tech Natural Infrastructure in Dryland Streams (NIDS)—including rock detention structures—to slow surface water, raise groundwater tables, and restore wetland function at a spring-fed wetland in Cebolla Canyon, New Mexico, U.S.A. Our results depict a Restoration Feedback Loop that captures stages of change from a healthy wetland in 1935, altered by 20th-century agriculture and grazing, to the re-establishment of the historical flow regime by 2024 documented through an 89-year archive of aerial imagery (1935–2024). By the end of our study period, the Spring-Fed Wetland had expanded by roughly 229% of the original 1935 area, to 4.13 ha. Using 40 years of satellite data, we assess changes in vegetation and hydrology with remote sensing indices. Spatial and temporal analyses reveal significant increases in vegetation greenness and wetness, particularly in an Expanded Wetland subregion, which exhibited ∼3.5x higher wetness and ∼1.5x higher greenness trends compared to adjacent areas. Monthly metrics highlight seasonal variability, with increases in greenness linked to monsoonal rainfall and lateral water redistribution, indicating that restoration impacts extend beyond the primary wetland. This study demonstrates the utility of cloud-based platforms like Google Earth Engine and USGS EarthExplorer for long-term monitoring of wetland restoration, while quantifying the efficacy of the ‘Zeedyk approach’ and demonstrating its potential as a scalable method to restore and conserve wetland meadows in other arid and semi-arid landscapes.
Between the 1950s and 1980s, wastewater generated at the Idaho National Laboratory contained Iodine-129 ( 129 I); this wastewater was discharged directly into the eastern Snake River Plain (ESRP) aquifer through a deep disposal well, unlined infiltration ponds, or leaked from distribution systems below industrial facilities. During 2021–22, the U.S. Geological Survey, in cooperation with the U.S. Department of Energy and the Idaho Department of Environmental Quality Idaho National Laboratory Oversight Program, collected groundwater samples from 64 monitoring wells in the ESRP aquifer, 6 of which are part of a multilevel monitoring system, to determine the concentration of 129 I in the groundwater. These samples were analyzed by accelerator mass spectrometry as part of a long-term ongoing study to track trends and occurrences of this carcinogenic, long-lived radionuclide in the environment. Concentrations ranged from slightly above the locally determined background concentration of 5.4×10 −6 picocuries per liter, to just below the U.S. Environmental Protection Agency’s maximum contaminant level of 1 picocurie per liter. Discharge of wastewater containing 129 I has been discontinued to the aquifer, and long-term trends from a subset (n=15) of sampled wells show decreasing 129 I concentrations over the last three decades. Concentrations of 129 I in groundwater from monitoring wells near facilities at the Idaho National Laboratory are affected by episodic recharge from an ephemeral surface-water source and by the fracture-flow dominated hydrologic regime in the ESRP aquifer. The spatially focused sampling effort has also identified a low-level 129 I plume that affects long-term water quality near and downgradient from the Advanced Test Reactor Complex in the southwestern part of the facility that had not been clearly defined in previous sampling efforts, although the definition of the plume is somewhat limited by available data.
On July 9, 1926, monitoring well Oe-151 at Woodgate, New York (USGS ID 433112075091501) recorded its first groundwater data. Since then, the well has provided water data nearly continuously and has now reached a 100-year milestone for data collection. The well is part of the U.S. Geological Survey (USGS) Climate Response Network (CRN), which is a national network of wells selected to monitor natural groundwater conditions. Well Oe-151 is the first well in the network to reach a 100-year Centennial milestone.
What is the U.S. Geological Survey National Groundwater Climate Response Network? The U.S. Geological Survey’s (USGS) National Water Monitoring Network is a network of networks that includes the Climate Response Network (CRN). The CRN is a network of wells selected to monitor the effects of climate variability, including droughts, on groundwater systems within the United States, Puerto Rico, and the U.S. Virgin Islands. The primary purpose of the CRN is to monitor these effects on groundwater levels in unconfined or near-surface confined aquifers that are minimally affected by pumping or other human-induced stresses.
The Little Sequatchie River and Pryor Cove Branch, in southern Tennessee, drain the eastern escarpment of the Cumberland Plateau to the Sequatchie River near the southern end of the Sequatchie Valley. The Little Sequatchie River is the largest tributary to the Sequatchie River by drainage area, covering over 120 square miles. The hydrology of the two drainage areas has been largely altered by karst processes, which has caused the majority of the streams to sink at the contact between the Mississippian Pennington Formation and the underlying Mississippian Bangor Limestone. The U.S. Geological Survey, in cooperation with the U.S. Fish and Wildlife Service and Tennessee Department of Environment and Conservation, initiated a study in 2021 to map the karst groundwater pathways in both watersheds in order to delineate recharge areas for several springs. One of these springs, Sequatchie Cave, represents a significant habitat for two Species of Greatest Conservation Need, the Glyphopsyche sequatchie (Sequatchie caddisfly) and the federally endangered Marstonia ogmorhaphe (royal marstonia). Springs and springflow-dominated streams in the Little Sequatchie River valley and Pryor Cove also provide water for agricultural practices and serve as a drinking water source for nearby communities. During the study, a total of 25 dye injections were conducted over eight rounds from January 2022 through March 2023. Dye traces from these injections helped to delineate recharge areas for six major springs, ranging from 7.3 to 65.2 square miles in area. The majority of the dye traces remained subsurface (from sinkpoint to recovery site) for long distances, with karst groundwater travelling nearly 8 miles before resurfacing. The dye traces also had rapid traveltimes, often travelling hundreds to thousands of feet per hour. The goal of this project was to provide scientific data related to karst groundwater pathways and spring recharge areas to aid State and Federal agencies in making informed decisions to protect and preserve this unique and vulnerable karst system.
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.