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.
There are thousands of thaw (thermokarst) lakes on the North Slope of Alaska, where they cover as much as 40% of the land area. Their very name recognizes the fact that they owe their origin to the impact they have on the ground thermal regime, but there have been few quantitative studies of the impact of the lakes on atmosphere-land interactions in this tundra region.
Gas hydrate-bearing sands on the Alaska North Slope (ANS) host minor volumes of crude oil whose impact on formation permeability has never been quantified relative to a gas hydrate reservoir system. Here, we combine in situ pressure-core microcomputed-tomography (μ-CT), thermogravimetric analysis (TGA), gas-chromatography–mass-spectrometry (GC–MS), scanning-electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), and Stokes-flow simulation to (i) quantify oil saturation, (ii) infer its source, and (iii) evaluate its influence on permeability after gas-hydrate dissociation. μ-CT and SEM-EDS imaging identify the host sediment as silty, with some grain-coating clay present. Eight representative 3003-voxel subvolumes extracted from a preserved pressure core (Core 15P-3c, from the B1 sand (unit B), HYDRATE 02 Geo Data Well (GDW), 886.96–887.29 m measured depth, MD) exhibit porosities of 0.34–0.37 and an average oil saturation of 0.09 ± 0.03, in comparison to an independent TGA-based estimate of 0.16. GC-MS fingerprints obtained from an adjacent depressurized core (Core 17P-1, 891.37–891.44 m MD), together with oil pore habit revealed by μ-CT, suggest that the oil is partly native but also includes contributions from mineral oil-based drilling fluid contamination. Further analysis of oil–porewater interactions and flushing tests confirms that native oil saturation in the analyzed samples ranges from 0.04 to 0.08 with an average of approximately 0.06, and the oil is highly immobile. Flow simulations demonstrate that the native oil saturation of ≈0.06 reduces permeability by approximately 50% due to both pore blockage and increased flow-path tortuosity. These findings reveal the potential presence of native oil within the ANS gas hydrate reservoir and suggest that future simulation models may need to account for its impact on permeability to improve long-term performance predictions of gas and water production.
Gas hydrate, a solid combination of gas (mostly methane in nature) and water molecules stable at low temperatures and elevated pressures, occurs naturally in marine and permafrost-associated environments. Gas hydrate reservoirs, such as those in the Alaska North Slope, have been considered potential energy resources for gas production. To understand the petrophysical and geo-mechanical characteristics of the reservoir, core samples retrieved from the site of the JOGMEC-DOE-USGS collaborative gas hydrate R&D project have been analyzed in the laboratory for their hydraulic and mechanical properties. This paper focuses on both seal and reservoir samples associated with the B1 and D1 sands, which are evaluated for index properties (including porosity, grain size distribution, liquid and plastic limits, specific surface area, and specific gravity), consolidation, permeability, and water retention. Furthermore, the reservoir core samples were tested with pore-filling, laboratory-grown tetrahydrofuran hydrate, in order to assess reservoir behavior during gas production from hydrates. Under simulated in situ stress conditions, the seal and hydrate-free reservoir cores had a permeability anisotropy ratio of k h /k v = 3.0−5.0, and k h /k v = 2.4−3.0 for the reservoir tetrahydrofuran hydrate-bearing cores. The data suggest that depressurizing the reservoir to induce hydrate dissociation alters the reservoir effective permeability in three ways: permeabilities decrease due to porosity lost (e.g., the initial reservoir thickness can decrease by up to 5% upon 7 MPa depressurization), permeability increases due to the loss of solid hydrate in the pore space, and permeability anisotropy k h /k v decreases in response to the evolving pore-space geometry. We show that given the simulated in situ gas hydrate saturations (i.e., S h = 32% in core 7P-2E and S h = 21% in core 20P-4), gas production from the dissociation of tetrahydrofuran hydrate in the two tested cores results in a net increase in effective permeability and a decrease in k h /k v . This study highlights the importance of investigating seal and reservoir sediments and the impacts of depressurization on the porosity and permeability responses during production.
Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.
Here, we investigated numerical simulation strategies for a long-term test of depressurization-induced gas production from the B1 Sand of Unit B at the Hydrate-01 Stratigraphic Test Well. The main objective of this study was to estimate fluid production rates (with emphasis on water production) under a variety of conditions and production scenarios and contribute new insights to the design and management of the field test. In the first part of the study, we investigated the system response to a three-step depressurization process using two limiting sets of flow properties─the expected maximum and minimum intrinsic and effective permeabilities─for the very heterogeneous reservoir. In the second part, we investigated the effect of the production interval length and placement within the formation relative to the boundaries of the hydrate-bearing unit. The best performing well configuration was used in the third part of the study, which used the most representative subsurface flow properties to investigate the effect of the depressurization strategy on the production performance. The best overall performance (largest gas production with modest water production and a strong response at the observation wells) was obtained with a 10 m-long well situated 3 m below the top of the formation and a three-step depressurization scheme at 15-day intervals to a terminal bottomhole pressure of 2.8 MPa. The overall production performance was enhanced by a faster rate of depressurization. Estimated water production rates in all cases were limited and easily manageable. None of the tested well configurations or depressurization strategies significantly reduced water production without also severely reducing gas production. In all the investigated cases, 95% of the long-term fraction of produced water was replenished by inflows from the boundaries and could not be reduced. These substantial water inflows are an unavoidable feature of HU-B and cannot be easily mitigated by a hydraulic control.
Explore the source record for details and available documents.
Three physiographic provinces comprise the North Slope of Alaska: the Arctic Mountains, the Arctic Foothills and the Arctic Coastal Plain Provinces. The features and processes in the Arctic Coastal Plain, a zone of continuous permafrost, are stressed in this paper. The evidence for and mechanisms of the geomorphic cycle are discussed starting with frost cracks. Frost cracks may form polygonal ground which leads to low-centered ice wedge polygons in areas having ice-rich permafrost. As the low-centered ice wedge polygons enlarge due to thermal erosion they may evolve into thaw lakes which are largely oriented in a northwest-southeast direction on the Arctic Coastal Plain. Eventual drainage of a deep lake may result in a closed-system pingo. Evidence of the various stages of the geomorphic cycle is ubiquitous on the Alaskan Arctic Coastal Plain and indicates the ice content of the permafrost in some areas.
Measurements of the C-band (wavelength = 5 cm) radar cross section of an area in the Brooks Range foothills on the North Slope of Alaska using images from the ERS-1 satellite show significant temporal changes. These changes are strongly correlated with elevation and hillslope orientation and are greatest on some of the elevated areas and weaker in river drainages. By constructing 'difference images' using various image pairs, and by analyzing climatological and hydrological data from the site, we conclude that the radar backscatter changes are largely due to changes in soil and vegetation liquid water content induced by freeze/thaw events. The correlation with topography in the difference images arises from the dependence of vegetation, organic layer thickness, and volumetric water content on hillslope position and orientation. These results demonstrate the viability of radar backscatter intensity comparisons using repeat-pass images as a means of change detection.
The North Slope of Alaska (NSA) atmospheric observatory, operated by the U.S. Department of Energy (DOE)’s Atmospheric Radiation Measurement (ARM) User Facility, is a measurement site in the Arctic that has been collecting crucial atmospheric data for more than 25 years. The central facility located in Utqiaġvik, Alaska (formerly known as Barrow) hosts a suite of instruments that are used to better understand arctic processes, which are often not well represented in earth system models. The NSA site sits only a few kilometers from the Arctic Ocean, which also makes it a prime location to study complex ocean-atmosphere-ice interactions. Arctic cloud and precipitation processes are also of scientific interest, and remote-sensing instruments including radars are a key component of the NSA instrument suite. One of the radars at NSA is the X-band Scanning ARM Precipitation Radar (XSAPR). This report evaluates one year of recent XSAPR data from April 2024 through April 2025 and details the process of generating b1-level data. This analysis marks the first effort by ARM staff to quality-control NSA XSAPR data with the goal of routinely producing b1-level data in the future depending on radar operations.
The Examining the Ice-Nucleating Particles from the North Slope of Alaska (ExINP-NSA) campaign was conducted at the National Oceanic and Atmospheric Administration’s (NOAA) Barrow Atmospheric Baseline Observatory (71.3230° N, 156.6114° W, “BRW” hereafter), next to the Atmospheric Radiation Measurement (ARM) user facility NSA site and ~ 6 km northeast of the town of Utqiaġvik. The location of the NSA site is shown in Figure 1. Our observing period began in October 2021 and continued until May 2024. This campaign was funded by the U.S. Department of Energy (DOE) Office of Science Early Career Research Program through grant number DE-SC001879. This grant contains funding for three ARM field campaigns. The previous field campaigns were performed at ARM’s Southern Great Plains atmospheric observatory in Oklahoma (SGP; 36.6073° N, 97.4876° W) and the ARM’s Eastern North Atlantic atmospheric observatory on Graciosa Island, Azores (ENA; 39.0916° N, 28.0257° W). The ExINP-NSA campaign aims included: • Determining the number concentration of ice-nucleating particles (n INPs ) active at temperatures spanning the range of heterogeneous freezing processes (from ≈−30 °C to 0 °C) using a combination of online and offline measurements, • Determining whether local meteorological conditions and/or synoptic scale air mass transport impact INP abundance and/or ice nucleation efficiency, • Examining if the physicochemical properties of INPs relate to aerosol chemistry, and • Assessing if there is a similarity in INP properties across three ARM sites Multi-seasonal datasets of INP abundance in the NSA region were delivered from this campaign. This campaign also allowed researchers to perform a comprehensive analysis of atmospheric INPs based on long- term ground-based measurements in the Alaskan Arctic. Our data and results from the ExINP-NSA campaign will help refine current earth system models. One of the stated goals of ARM is to advance aerosol-cloud ice interaction, which will be a direct result of this campaign. Current earth system models poorly represent INPs, and the 15-minute time resolution data over several seasons generated during this campaign will provide an invaluable resource, especially combined with the datasets generated during two previous ARM ExINP campaigns. These datasets will allow for a greater understanding of ice nucleation processes as they may (or may not) relate to local meteorological processes and aerosol chemistry and will eventually help further the understanding of the Earth’s atmospheric processes and energy balance.
This dataset provides shapefiles and trained machine learning models used for aufeis detection at four sites on the North Slope of Alaska. It includes reference data for evaluating Landsat-based detection methods, supporting research on remote sensing approaches for identifying aufeis. The ReferenceData folder contains ArcGIS shapefiles of semi-automated land cover classifications for 217 Landsat Collection 2 images, categorizing pixels into six classes: aufeis, snow, ground, none, water, and cloud. The SiteBuffers.zip file includes 10-kilometer buffer shapefiles defining regions of interest around four aufeis fields (Canning21, FH1, Firth, and Kuparuk), used to test three detection techniques. Additionally, the TrainedRFModels folder contains six pre-trained Scikit-Learn Random Forest classifiers (100 trees, max depth = 30) designed to predict aufeis presence in Landsat Collection 2 Surface Reflectance images using Red, Blue, SWIR2, NDVI, and NDWI bands. This dataset supports the development and validation of remote sensing methods for mapping aufeis in Arctic environments.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).
Particles are removed from the atmosphere through both wet and dry deposition. These processes are poorly understood, though they constitute important uncertainties in climate and air quality models. This project aims to use observational constraints on particle fluxes to improve model representations of dry deposition. Within this data set, we measured fluxes of size resolved particles with an ultrahigh sensitivity aerosol spectrometer (UHSAS) at the ARM facility’s North Slope Alaska site. Measurements were made at the meteorological tower between 9 Sept and 25 Oct 2021. This site is a coastal tundra location and received snow during the project.
A single-channel (56.7 GHz) scanning radiometer was deployed in August 2022 at the Atmospheric Radiation Measurement (ARM) North Slope of Alaska site near Utqiaġvik. The radiometer is designed to provide temperature profiles between 0 and 1 km every 5 min. Averaging kernels show that this single-channel radiometer, taking observations at 10 discrete elevation angles, yields approximately the same information as a seven-channel V-band radiometer scanning three elevation angles. The instrument is able to reproduce the occurrence of temperature inversions between the surface and 1 km and their strength showing a correlation of 0.85, bias of −0.6 K, and slope of 1.04 with respect to radiosondes. Uncertainty in the inversion base height varies from 50 m near the surface to ∼300 m above 0.4 km when compared with radiosondes. Conversely, the inversion top height is overestimated and has higher uncertainty due to the degrading effects of the averaging kernels on the vertical resolution of the retrievals. Thanks to the high temporal resolution of the retrievals, the diurnal cycle of boundary layer temperature was evaluated showing that the radiometer can capture some aspects of the boundary layer thermal structure. The present analysis provides an overview of the capabilities of this simple observing configuration for selected applications.
The U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) User Facility operates an X-band Scanning ARM Precipitation Radar (XSAPR) at the North Slope of Alaska (NSA) atmospheric observatory in Utqiaġvik. This radar provides year-round dual-polarization measurements of snow and ice in an arctic environment. Over the past few years, there has been greater focus on the NSA XSAPR, and recent efforts by ARM staff have resulted in quality-controlled, b-level data (Rocque et al. 2026b). This report outlines the processing of the second set of b-level data released for this radar from May 2025 through April 2026 following similar methods from Rocque et al. (2026b).
Arrays of instruments for monitoring winter precipitation (snowfall) and snow cover on the ground installed in 2017 at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s North Slope of Alaska (NSA) C1 and at Oliktok Point sites became operational in 2018. In 2022, the instruments from Oliktok Point were moved to NSA E12, about 5.4 km south of NSA C1, where two arrays now operate. The instrument arrays monitor wind speed and direction, snow depth at multiple locations, the horizontal flux of blowing snow, and the number and fall speed of hydrometeors. The arrays are monitored using digital cameras. Collectively, the instruments produce a wealth of data on falling and accumulated snow, but, as with any instrument array, some data are more reliable and accurate than others. In this document, we present our findings on the efficacy, accuracy, and reliability of each type of instrument. Overall, a key finding is that no single instrument provides sufficient information to determine the source of falling snow particles nor the cause of changes in snow depth. However, if used in concert, the instruments produce a reliable understanding of the processes affecting the snow cover depth distribution and the true winter precipitation.
The Global Electric Circuit (GEC) of the atmosphere provides a unique perspective on Earth’s changing climate. Monitoring this global electrical signature provides details of the global nature of electrified clouds and thunderstorms. The North Slope of Alaska (NSA) is a unique location for collecting these electric field measurements. Besides being at the opposite pole from many previous electric field measurements taken at Russia’s Vostok station in Antarctica, this site provides a rare opportunity to use other instruments at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s NSA observatory at Utqiaġvik (formerly known as Barrow), such as the Ka-band ARM Zenith Radar (KAZR), upward-facing micropulse lidar (MPL), vertical profile of meteorological measures, and other aerosol measurements. Since 2017, Texas A&M University at Corpus Christi (TAMUCC) has established an observational instrument package, including CS110 electric fields and anemometer, at ARM’s NSA observatory at Utqiaġvik. This enabled the unique opportunity to not only provide information about the global signature of the GEC, but also the physical inputs to the local electric field, by analyzing the physical properties of the simultaneous cloud, wind, and aerosol properties occurring with the vertical electric field.
The Global Electric Circuit (GEC) of the atmosphere provides a unique perspective on Earth’s changing climate. Monitoring this global electrical signature provides details of the global nature of electrified clouds and thunderstorms. The North Slope of Alaska (NSA) is a unique location for collecting these electric field measurements. Besides being at the opposite pole from many previous electric field measurements taken at Russia’s Vostok station in Antarctica, this site provides a rare opportunity to use other instruments at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s NSA observatory at Utqiaġvik (formerly known as Barrow), such as the Ka-band ARM Zenith Radar (KAZR), upward-facing micropulse lidar (MPL), vertical profile of meteorological measures, and other aerosol measurements. Since 2017, Texas A&M University at Corpus Christi (TAMUCC) has established an observational instrument package, including CS110 electric fields and anemometer, at ARM’s NSA observatory at Utqiaġvik. This enabled the unique opportunity to not only provide information about the global signature of the GEC, but also the physical inputs to the local electric field, by analyzing the physical properties of the simultaneous cloud, wind, and aerosol properties occurring with the vertical electric field.
In this paper, we explore scaling and data assimilation-related issues associated with utilizing passive microwave satellite observations of Cold Lands-in this case, the climatologically and ecologically sensitive arctic tundra. Our approach expands on our earlier work using a one-year dataset from the Radiobrightness Energy Balance Experiment-3 (REBEX-3). REBEX-3 featured a tower-based SSM/I (Special Sensor Microwave/Imager) simulator deployed on the North Slope of Alaska in 1994-95. Two findings are significant here. First, a comparison of tower and satellite signatures at 19 and 37 GHz strongly suggested that the North Slope is radiometrically homogeneous for spatial scales up to SSM/I footprints (approximately 25 km), an unusual and valuable characteristic for monitoring and retrieving land surface conditions. And second, at the plot scale, signatures of snow/no-snow and freeze/thaw transitions were identifiable for tussock tundra land cover, so that even snow-free frozen tundra could be unambiguously distinguished from tundra covered with dry snow, another unusual and valuable characteristic. We present results from analyzing satellite brightness signatures of selected North Slope pixels corresponding to instrumented sites along a transect from the Brooks Range to the Arctic Ocean. A custom EASE (Equal Area Scalable Earth)-Grid processor was used to extract SSMJI data for every orbit with observations of this region during the 1994-95 year. The resulting high temporal-resolution (4-8 points/day), gridded data were then analyzed for evidence of the same diurnal and seasonal signatures seen at the plot scale (through micrometeorological and/or brightness data). Differences between satellite and tower brightness observations are quantified for various conditions at the REBEX-3 site. Such differences from the less-frequent and/or larger-scale satellite observations represent a form of input 'noise' in data assimilation applications. For the other sites, the performance of snow/no-snow and freeze/thaw discriminators vs. ground truth represents an opportunity to gauge the homogeneity of other pixels.