SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Develop an efficient parallel implementation and assess the performance of local time-stepping (LTS) schemes for the shallow water equations in the Model for Prediction Across Scales (MPAS).
The US Department of Energy (DOE) conducts climate simulations on some of the world’s largest supercomputers. These exascale machines use heterogeneous architectures with both CPUs and GPUs, and scientific codes must adapt to make full use of this computing power. Los Alamos National Lab is developing Omega: The Ocean Model for E3SM Global Applications, which is specifically designed for modern exascale computers. It uses external libraries that have been optimized for a variety of architectures to run on different supercomputers. Omega is an unstructured-mesh ocean model based on TRiSK numerical methods. It will be the new ocean component of the DOE’s Energy Exascale Earth System Model (E3SM). The algorithms in Omega follow those of the current ocean component, MPAS-Ocean, but it will be written in C++ rather than Fortran to take advantage of the Kokkos performance portability library. Omega spatial operators are written as Kokkos kernels to run efficiently on both CPUs and GPUs. Work on Omega began in 2023 with a new C++ framework for unstructured mesh partitioning, halo exchanges, parallel IO, and Kokkos interfaces. The current version, Omega-0, is being developed to solve the shallow water equations and at present includes all of the tendency terms but not time stepping. Here we share the results of Omega-0 verification and performance testing. Verification includes unit tests implemented with CTest as well as convergence tests in Polaris, an in-house python package with a large suite of test problems. Performance tests compare simulations conducted on CPUs versus GPUs and across different architectures: tests are run on Frontier, which has AMD “Optimized 3rd Gen EPYC” CPUs and AMD MI250X GPUs, as well as Perlmutter, which is composed of AMD EPYC 7763 CPUs and NVIDIA A100 GPUs.
Benthic invertebrates play critical ecosystem roles including the breakdown of organic matter, sediment mixing, and nutrient cycling. In the lower Columbia River and estuary, benthic invertebrates provide a foraging resource for threatened and endangered juvenile salmon. Driven by the goal to create low velocity, shallow water, and riparian shrub habitats to benefit juvenile salmon, the USACE placed 237,000 CY of dredged material, resulting in the deposition of 13.5 acres of sand on the off-channel margins at Woodland Islands. The Pacific Northwest National Laboratory (PNNL) designed and implemented an action effectiveness research study to understand how dredged material placement at Woodland Islands affected sediment conditions and benthic invertebrate assemblages. A Before-After-Control-Impact, or BACI, study design was used to evaluate the response of environmental conditions and benthic invertebrates to dredged material placement. The design included one impact site and two control sites. Sampling occurred for two years prior to dredged material placement and two years after. Our study found that spatial variation was a significant factor for both environmental and biological response variables which suggests local conditions are important considerations for mechanisms affecting benthic assemblages. We found that off-channel habitats across all locations sampled were producing benthic invertebrates, many of which are common prey items for juvenile salmon and steelhead—e.g., insects, chironomids, crustaceans, and corophium. There was a significant BACI effect (i.e., an effect at the dredged material placement site, relative to conditions at the control sites) for concentration of carbon and ammonium in sediment, but not for phosphorous concentrations or for percent sand. The estimated abundance for the three invertebrate response variables—total abundance, total chironomid abundance, and total corophium abundance—was significantly lower at the impact site after dredged material placement, compared to the control sites. At Woodland Island, the estimated mean abundance for all invertebrates combined decreased 28% after dredged material placement. Estimated mean abundances of chironomid and corophium decreased by 8% and 88%, respectively. While invertebrate abundances were lower after dredged placement, the composition of benthic invertebrates was similar before and after placement suggesting that as the new habitat feature evolves recolonization will likely follow. These findings provide a foundation for understanding potential benefits and consequences of repurposing dredged material for habitat creation in the LCRE.
The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.
The research goal of this Department of Energy (DOE) Office of Fusion Energy Science (OFES) project was to investigate the electrified liquid surface dynamics using a shallow water model. The project was an extension of the Princeton Collaborative Research Facility (PCRF) project to work with Dr. Mikhail Shneider at Princeton University.
With changing demands in regulation, understanding methane emissions from offshore oil and gas production infrastructure has become increasingly important. Reported emissions from facilities in the Gulf of Mexico range from zero to thousands of tons of methane per hour, but these is currently no clear understanding of how this range compares to expected emissions from normally operating facilities. To generate realistic emission estimates, we create two bottom-up models that simulate emissions from facilities operating in the Gulf of Mexico. We estimate type 1 prototypical facilities (typically unmanned, older, lower-producing platforms in shallow water with little processing equipment, compressors, or storage tanks) to emit an average of 13 kg CH 4 h −1 , which corresponds to a loss of 2.7% of the average facility production. Type 2 prototypical facilities (continuously manned, higher production and operate in deeper water with processing equipment, oil storage tanks, compressors and power generation) emit an average of 88 kg CH 4 h −1 , which corresponds to a loss of 2.5% of production. The average measured emission from type 1 facilities was 18 kg CH 4 h −1 with a median production loss estimated at 8%. The average measured emission from type 2 facilities was 36 kg CH 4 h −1 with a median production loss estimated at 2.4%. Using emission factors that consider the long-tail emission distribution partly reconciles the difference between modelled and measured emission estimates, but we suggest the current the fugitive emission estimate may be an underestimate and more data on the number and size of fugitive emissions could explain differences between the modelled and measured emission estimate. We suggest the bottom-up approach described here that uses production data coupled with facility equipment could be used to identify facilities that have abnormally large measured emissions, caused by methodological failure or larger than expected fugitive emissions, which should be targeted for further evaluation resulting in remeasurement or identification of source type so that a more accurate estimates can be made on the absolute emission.
One of the biggest risks to safety on offshore platform safety is the ignition of high-pressure natural gas streams. Currently, the size and number of fugitive emissions on offshore platforms is unknown and methods used to detect fugitives have significant shortcomings. To investigate the frequency, size, and potential impact of fugitives, a data collection exercise was conducted using incidents reported, leak survey data, and independent measurements. The size and number of fugitives on offshore facilities were simulated to investigate likely areas of safety concern. Incident reports indicate in 2021 there were 113 reports of gas leaks on 1119 offshore facilities, suggesting 0.02 fugitives per Type 1 facility (older, shallow-water platforms) and 0.31 fugitives per Type 2 facility (larger deeper-water facilities). Leak survey data report 12 fugitives per Type 1 facility (average emission 0.6 kg CH 4 h −1 leak −1 ) and 15 fugitives per Type 2 facility (average emission 1.5 kg CH 4 h −1 leak −1 ). Reconciliation of direct measurements with a bottom-up model suggests that the number of fugitive emissions generated from the leak report data is an underestimate for Type 1 platforms (44 fugitives facility −1 ; average emission 0.6 kg CH 4 h −1 leak −1 ) and in general agreement for the Type 2 platforms (15 fugitives facility −1 ; average emission 1.5 kg CH 4 h −1 leak −1 ). Analysis of the fugitive emission rates on an offshore platform suggests that gas will not collect to explosive concentration if any air movement is present (>0.36 mph); however, large volumes of air (~600 m 3 ) near representative leaks on the working deck could become explosive in hour-long zero-wind conditions. We suggest that wearable technology could be employed to indicate gas build up, safety regulations amended to consider low-wind conditions and real-world experiments are conducted to test assumptions of air mixing on the working deck.
A new sweep-search algorithm (SSA) is developed and tested to identify the channel geometry transitions responsible for numerical convergence failure in a Saint-Venant equation (SVE) simulation of a large-scale open-channel network. Numerical instabilities are known to occur at “sharp” transitions in discrete geometry, but the identification of problem locations has been a matter of modeler’s art and a roadblock to implementing large-scale SVE simulations. The new method implements techniques from graph theory applied to a steady-state 1D shallow-water equation solver to recursively examine the numerical stability of each flowpath through the channel network. The SSA is validated with a short river reach and tested by the simulation of ten complete river systems of the Texas–Gulf Coast region by using the extreme hydrological conditions recorded during hurricane Harvey. The SSA successfully identified the problematic channel sections in all tested river systems. Subsequent modification of the problem sections allowed stable solution by an unsteady SVE numerical solver. The new SSA approach permits automated and consistent identification of problem channel geometry in large open-channel network data sets, which is necessary to effectively apply the fully dynamic Saint-Venant equations to large-scale river networks or for city-wide stormwater networks.
Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.
Abstract. Atmospheric rivers (ARs) are extreme weather events that can alleviate drought or cause billions of US dollars in flood damage. By transporting significant amounts of latent energy towards the poles, they are crucial to maintaining the climate system's energy balance. Since there is no first-principle definition of an AR grounded in geophysical fluid mechanics, AR identification is currently performed by a multitude of expert-defined, threshold-based algorithms. The variety of AR detection algorithms has introduced uncertainty into the study of ARs, and the thresholds of the algorithms may not generalize to new climate datasets and resolutions. We train convolutional neural networks (CNNs) to detect ARs while representing this uncertainty; we name these models ARCNNs. To detect ARs without requiring new labeled data and labor-intensive AR detection campaigns, we present a semi-supervised learning framework based on image style transfer. This framework generalizes ARCNNs across climate datasets and input fields. Using idealized and realistic numerical models, together with observations, we assess the performance of the ARCNNs. We test the ARCNNs in an idealized simulation of a shallow-water fluid in which nearly all the tracer transport can be attributed to AR-like filamentary structures. In reanalysis and a high-resolution climate model, we use ARCNNs to calculate the contribution of ARs to meridional latent heat transport, and we demonstrate that this quantity varies considerably due to AR detection uncertainty.
This paper introduces Omega, the Ocean Model for E3SM Global Applications. Omega is a new ocean model designed to run efficiently on high performance computing (HPC) platforms, including exascale heterogeneous architectures with accelerators, such as Graphics Processing Units (GPUs). Omega is written in C and uses the Kokkos performance portability library. These were chosen because they are well-supported and will help future-proof Omega for upcoming HPC architectures. Omega will eventually replace the Model for Prediction Across Scales-Ocean (MPAS-Ocean) in the US Department of Energy's (DOE's) Energy Exascale Earth System Model (E3SM). Omega runs on unstructured horizontal meshes with variable-resolution capability and implements the same horizontal discretization as MPAS-Ocean. This work documents the design and performance of Omega Version 0.1.0 (Omega-V0), which solves the shallow water equations with passive tracers and is the first step towards the full primitive equation ocean model. On Central Processing Units (CPUs), Omega-V0 is 1.4 times faster than MPAS-Ocean with the same configuration. Omega-V0 is more efficient on GPUs than CPUs on a per-watt basis – by a factor of 5.3 on Frontier and 3.6 on Aurora, two of the world's fastest exascale computers.
The Technology Commercialization Fund (TCF) award is built upon NREL's previous research exploring a bottom-fixed variable-geometry oscillating surge wave energy converter (VGOSWEC); however, unlike previous investigations, the VGOSWEC will be raised off the sea floor. Researchers believe that advances in offshore pilings and foundations, such as 3D printing, will provide low cost and easily deployable structures. There are several benefits to designing a raised VGOSWEC such as 1) reducing issues with sediment transport and environmental impacts near shore, 2) less expensive deployment and installation as work boats will not risk grounding in shallow water, and 3) deeper water deployments improve the wave resource.
Dryland ecosystems occupy a vast swath of the terrestrial land surface and exert a sizeable impact on the cycling of carbon and water globally. These biomes are characterized by tightly coupled carbon and water cycles that respond rapidly to transient pulses in water availability. However, there exist many mechanistic uncertainties regarding the environmental drivers of, and linkages between, plant and ecosystem processes. Thus, drylands are often poorly represented in many vegetation and land surface models. An enhanced understanding of dryland ecosystem function is limited by the lack of long-term, co-located, and frequent measurements of plant and ecosystem processes. At a piñon-juniper woodland in southeastern Utah, USA, we collected a continuous dataset of meteorological conditions, soil water potential from surface to bedrock, tree water potential, and ecosystem carbon and water fluxes from eddy covariance. We found that predawn and midday tree water potential and daily ecosystem fluxes were highly sensitive to fluctuations in soil water availability, particularly in shallower layers, and that daytime variability in atmospheric drivers only loosely controlled these processes. The strong connections between shallow soil water potential, tree water potential, and ecosystem fluxes occurred because of the dominant role of precipitation pulses in driving vegetation activity, as even small pulses of moisture stimulated shallow soil water potential, tree water potential, and evapotranspiration for between 1 and 2 weeks. Carbon fluxes (net ecosystem exchange and gross primary productivity) were sensitive to precipitation pulses for longer, up to 3 weeks. Our results highlight that improved monitoring and sensing of shallow soil moisture can greatly enhance our understanding of dryland ecosystem function. Here, a better mechanistic understanding of the impacts of precipitation pulses is also needed to improve vegetation modeling of dryland ecosystems.
The profile of the liquid water content (LWC) in clouds provides fundamental information for understanding the internal structure of clouds, their radiative effects, propensity to precipitate, and degree of entrainment and mixing with the surrounding environment. In principle, differential absorption techniques based on coincident dual-frequency radar reflectivity observations have the potential to provide the LWC profile. Previous differential frequency radar reflectivity (DFR) efforts were challenged by the fact that the measurable differential attenuation for small quantities of LWC is usually comparable to the system measurement error. This typically renders the retrieval impractical, as the uncertainty can become many times greater than the retrieved value itself. Theoretically, this drawback can be mitigated following two interconnected approaches: (1) increasing the frequency separation between the dual-frequency radar system to measure greater differential attenuation and (2) increasing the radar operating frequency to reduce the instrument measurement random error. Our recently developed 239 GHz radar was deployed during the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) along with a variety of collocated remote sensing and in situ instruments. We have combined Ka-band (35 GHz) and G-band (239 GHz) observations to retrieve the LWC from more than 100 vertical profiles of shallow clouds with typical amounts of LWC smaller than 1 g m -3 . We theoretically and experimentally demonstrate that the Ka-band and G-band pair of frequencies offers at least a 65 % relative improvement in the LWC retrieval sensitivity compared to previous works reported in the literature using lower-frequency radars. This new technique provides a missing capability to determine the LWC in the challenging low liquid water path (LWP) range (< 200 g m -2 ) and suggests a way forward to characterize microphysical and dynamical processes more precisely in shallow clouds.
Based on the projected increase in hurricane landfall frequency on the middle to lower US east coast, we examined the crucial role of the estuarine wind field in exacerbating coastal flooding. A regionally refined atmospheric and two high-resolution hydrology and ocean models are integrated to provide plausible and physically consistent ensembles of hurricane events and the associated flooding inside the Delaware Bay and River, a US mid-Atlantic estuary. Model results show that the hurricane propagation direction, estuarine geometry, remote surge from the open ocean, and direct nearshore upwind stress could magnify the flood magnitude. More specifically, inland-bound tracks that make landfall before reaching the mid-Atlantic coast produce a more significant surge within Delaware Bay than the shore-parallel tracks, where the estuarine wind direction plays the primary role in surge amplification. Ultimately, this study emphasized the need for integrated models to capture the nonlinear dynamics and interactions in flood hazard modeling.
In seasonally dry tropical forests, plant functional type can be classified as deciduous low wood density, deciduous high wood density, or evergreen high wood density species. While deciduousness is often associated with drought-avoidance and low wood density is often associated with tissue water storage, the degree to which these functional types may correspond to diverging and unique water use strategies has not been extensively tested. We examined (a) tolerance to water stress, measured by predawn and mid-day leaf water potential; (b) water use efficiency, measured via foliar δ 13 C; and (c) access to soil water, measured via stem water δ 18 O. We found that deciduous low wood density species maintain high leaf water potential and low water use efficiency. Deciduous high wood density species have lower leaf water potential and variable water use efficiency. Both groups rely on shallow soil water. Evergreen high wood density species have low leaf water potential, higher water use efficiency, and access alternative water sources. These findings indicate that deciduous low wood density species are drought avoiders, with a specialized strategy for storing root and stem water. Deciduous high wood density species are moderately drought tolerant, and evergreen high wood density species are the most drought tolerant group. Our results broadly support the plant functional type framework as a way to understand water use strategies, but also highlight species-level differences.
Abstract A new infaunal pinnotherid crab, Pinnixulala heardin. sp., is described from muddy substrates of the middle to inner continental shelf of the northern Gulf of Mexico. Molecular genetic analysis supports its separation from P. retinens (Rathbun, 1918), the only other western Atlantic member of the genus and a species most commonly found living as a symbiont of large polychaete worms. Both western Atlantic species are clearly separated from several described and undescribed populations of congeners occurring along shorelines of the eastern Pacific. We herein include a rediagnosis of P. retinens to facilitate morphological comparisons with P. heardin. sp. Unlike P. retinens, which is more common in shallow inshore waters and reaches larger sizes, the new species lacks development of a strong distally directed spine on the fourth pereopod ischium and has chelae that are at least as heavily developed as the second and third pereopods. In addition, the mature male first gonopods in P. heardin. sp. are tapered over their terminal one-third and directed distally, rather than being curved in their terminal one-third with the tips directed mesially, as they are in P. retinens. The new species is invariably collected from samples including infaunal worms and other invertebrates, but the specific hosts remain to be determined
In this article, federal offshore waters in the Gulf of Mexico are of interest for large-scale geologic carbon storage (GCS). However, more than 80,000 offshore oil and gas wells exist in the region, which could impact the integrity of sealing intervals. In this study, we propose a screening methodology for ranking offshore legacy wells based on the challenge they may present to GCS. The methodology relies on the review of well regulatory records to 1) identify leakage pathways and assess the potential hazards that wells pose to planned GCS operations, 2) evaluate well features that impact the accessibility of wells to determine the feasibility of potential corrective actions, and 3) rank wells based on the overall challenge they may pose for GCS. We demonstrate our framework by evaluating the construction and abandonment of 156 wells across eight areas of interest (AOIs) in shallow federal waters along the Texas Gulf Coast. The majority (99.3 %) of wells considered were constructed and plugged in a manner that did not isolate prospective GCS targets in the Upper and Lower Miocene formations and may potentially require a challenging or uncertain corrective action prior to GCS. Dataset trends suggest that the observed well construction and plugging designs may be common in shallow offshore federal waters along the Texas Gulf Coast. Consequently, operators pursuing offshore GCS projects in the region may consider selecting areas that avoid challenging wells or performing robust evaluations of legacy well leakage risks to plan corrective action prior to CO 2 injection.