Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data depth”

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.

At least 1,567 records · Page 87

Seismic Ambient Noise Analyses Reveal Changing Temperature and Water Signals to 10s of Meters Depth in the Critical Zone

The critical zone sustains terrestrial life, but we have few tools to explore it efficiently beyond the first few meters of the subsurface. Using analyses of high-frequency ambient seismic noise from densely spaced seismometers deployed in the forested Shale Hills subcatchment of the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO), we show that temporal changes in seismic velocities at depths from ~1 m to tens of m can be detected. Furthermore, these changes are driven by variations at the land surface. The Moving-Window Cross-Spectral (MWCS) method was employed to measure seismic-velocity changes in coda waves at hourly resolution in 10 different frequency bands. We observed a diurnal signal, a seasonal signal, and a meteorological-event-based signal. These signals were compared to time-series measurements of precipitation, well water levels, soil moisture, soil temperature, air temperature, latent heat flux, and air pressure in the heavily instrumented catchment. Most of the velocity changes can be explained by variations in temperature that result in thermoelastic strains that propagate to depth. But some double minima in seismic velocity time-series observed after large rain events were attributed in part to the effects of water infiltration. These results show that high-frequency ambient noise data may in some locations be used to detect changes in the critical zone from ~1 to ~100 m or greater depth with hourly resolution. But interpretation of such data requires multiple environmental data sets to deconvolve the complex interrelationships among thermoelastic and hydrological effects in the subsurface critical zone.

54 ENVIRONMENTAL SCIENCES↗

Interactive Visualization of High-Dimensional Petascale Ocean Data

We describe an application for interactive visualization of 5 petabytes of time-varying multivariate data from a high-resolution global ocean circulation model. The input data are 10311 hourly (ocean time) time steps of various 2D and 3D fields from a 22-billion point 1/48- degree “lat-lon cap” configuration of the MIT General Circulation Model (MITgcm). We map the global horizontal model domain onto our 128-screen (8x16) tiled display wall to produce a canonical tiling with approximately one MITgcm grid point per display pixel, and using this tiling we encode the entire time series for multiple native and computed scalar quantities at a collection of ocean depths. We reduce disk bandwidth requirements by converting the model’s floating point data to 16-bit fixed point values, and compressing those values with a lossless video encoder, which together allow synchronized playback at 24 time steps per second across all 128 displays. The application allows dynamic assignment of any two encoded tiles to any display, and has multiple interfaces for quickly specifying various orderly arrangements of tiles. All subsequent rendering is done on the fly, with run time control of colormaps, transfer functions, histogram equalization, and labeling. The two data streams on each screen can be rendered independently and combined in various ways, including blending, differencing, horizontal/ vertical wipes, and checkerboarding. The two data streams on any screen can optionally be displayed as a scatterplot in their joint attribute space. All scatterplots and map-view plots from the same x/y location and depth are linked so they all show the current brushable selection. Ocean scientists have used the system, and have found previously unidentified features in the data.

Interactive↗

Performance of the Lidar Design and Data Algorithms for the GLAS Global Cloud and Aerosol Measurements

The Geoscience Laser Altimeter System (GLAS) launched in early 2003 is the first polar orbiting satellite lidar. The instrument design includes high performance observations of the distribution and optical scattering cross sections of atmospheric clouds and aerosol. The backscatter lidar operates at two wavelengths, 532 and 1064 nm. For the atmospheric cloud and aerosol measurements, the 532 nm channel was designed for ultra high efficiency with solid state photon counting detectors and etalon filtering. Data processing algorithms were developed to calibrate and normalize the signals and produce global scale data products of the height distribution of cloud and aerosol layers and their optical depths and particulate scattering cross sections up to the limit of optical attenuation. The paper will concentrate on the effectiveness and limitations of the lidar channel design and data product algorithms. Both atmospheric receiver channels meet and exceed their design goals. Geiger Mode Avalanche Photodiode modules are used for the 532 nm signal. The operational experience is that some signal artifacts and non-linearity require correction in data processing. As with all photon counting detectors, a pulse-pile-up calibration is an important aspect of the measurement. Additional signal corrections were found to be necessary relating to correction of a saturation signal-run-on effect and also for daytime data, a small range dependent variation in the responsivity. It was possible to correct for these signal errors in data processing and achieve the requirement to accurately profile aerosol and cloud cross section down to 10-7 llm-sr. The analysis procedure employs a precise calibration against molecular scattering in the mid-stratosphere. The 1064 nm channel detection employs a high-speed analog APD for surface and atmospheric measurements where the detection sensitivity is limited by detector noise and is over an order of magnitude less than at 532 nm. A unique feature of the GLAS is a full acquisition of the surface return pulse, which has important application to the atmospheric transmission retrieval.

Spinhirne, James D.↗

Global Variations in Regolith Properties on Asteroid Vesta from Dawn's Low-Altitude Mapping Orbit

We investigate the depth, variability, and history of regolith on asteroid Vesta using data from the Dawn spacecraft. High-resolution (15-20 m pixel(sup -1)) Framing Cameraimages are used to assess the presence of morphologic indicators of a shallow regolith,including the presence of blocks in crater ejecta, spur-and-gully-type features in crater walls,and the retention of small (less than 300 m) impact craters. Such features reveal that the broad,regional heterogeneities observed on Vesta in terms of albedo and surface composition extend to the physical properties of the upper approx. 1 km of the surface. Regions of thin regolithare found within the Rheasilvia basin and at equatorial latitudes from approx. 0-90 deg. E and approx.260-360 deg. E. Craters in these areas that appear to excavate material from beneath the regolithhave more diogenitic (Rheasilvia, 090 deg. E) and cumulate eucrite (260-360 deg. E) compositions.A region of especially thick regolith, where depths generally exceed 1 km, is found from approx.100-240 deg. E and corresponds to heavily cratered, low-albedo surface with a basaltic eucritecomposition enriched in carbonaceous chondrite material. The presence of a thick regolithin this area supports the idea that this is an ancient terrain that has accumulated a larger component of exogenic debris. We find evidence for the gardening of crater ejecta towardmore howarditic compositions, consistent with regolith mixing being the dominant form of "weathering" on Vesta.

asteroid↗

Use of ERTS-1 data: Summary report of work on ten tasks

The author has identified the following significant results. Depth mapping's for a portion of Lake Michigan and at the Little Bahama Bank test site have been verified by use of navigation charts and on-site visits. A thirteen category recognition map of Yellowstone Park has been prepared. Model calculation of atmospheric effects for various altitudes have been prepared. Radar, SLAR, and ERTS-1 data for flooded areas of Monroe County, Michigan are being studied. Water bodies can be reliably recognized and mapped using maximum likelihood processing of ERTS-1 digital data. Wetland mapping has been accomplished by slicing of single band and/or ratio processing of two bands for a single observation date. Both analog and digital processing have been used to map the Lake Ontario basin using ERTS-1 data. Operating characteristic curves were developed for the proportion estimation algorithm to determine its performance in the measurement of surface water area. The signal in band MSS-5 was related to sediment content of waters by modelling approach and by relating surface measurements of water to processed ERTS data. Radiance anomalies in ERTS-1 data could be associated with the presence of oil on water in San Francisco Bay, but the anomalies were of the same order as those caused by variations in sediment concentration and tidal flushing.

Thomson, F. J.↗

Mars 2020 Entry, Descent and Landing Instrumentation (MEDLI2)

This paper will introduce Mars Entry Descent and Landing Instrumentation (MEDLI2) on NASA's Mars2020 mission. Mars2020 is a flagship NASA mission with science and technology objectives to help answer questions about possibility of life on Mars as well as to demonstrate technologies for future human expedition. Mars2020 is scheduled for launch in 2020. MEDLI2 is a suite of instruments embedded in the heatshield and backshell thermal protection systems of Mars2020 entry vehicle. The objectives of MEDLI2 are to gather critical aerodynamics, aerothermodynamics and TPS performance data during EDL phase of the mission. MEDLI2 builds up the success of MEDLI flight instrumentation on Mars Science Laboratory mission in 2012. MEDLI instrumentation suite measured surface pressure and TPS temperature on the heatshield during MSL entry into Mars. MEDLI data has since been used for unprecedented reconstruction of aerodynamic drag, vehicle attitude, in-situ atmospheric density, aerothermal heating, transition to turbulence, in-depth TPS performance and TPS ablation. [1,2] In addition to validating predictive models, MEDLI data has highlighted extra margin available in the MSL forebody TPS, which can potentially be used to reduce vehicle parasitic mass. MEDLI2 expands the scope of instrumentation by focusing on quantities of interest not addressed in MEDLI suite. The type the sensors are expanded and their layout on the TPS modified to meet these new objectives. The paper will provide key motivation and governing requirements that drive the choice and the implementation of the new sensor suite. The implementation considerations of sensor selection, qualification, and demonstration of minimal risk to the host mission will be described. The additional challenges associated with mechanical accommodation, electrical impact, data storage and retrieval for MEDLI2 system, which extends sensors to backshell will also be described.

Thermal Protection System↗

Initial Validation of a Gas-Granular Flow Solver Using a Subscale, Reduced Pressure Plume Surface Interaction Ground Test

With NASA’s goal to land the next humans on the lunar surface in the next few years, it is vitally important to have a better understanding of the plume surface interaction (PSI) between the landing vehicles and the lunar regolith. The Fluid Dynamics Branch at NASA/MSFC has previously used the gas-granular flow solver Loci/GGFS to qualitatively predict crater formation due to PSI effects in a lunar (near vacuum) ambient environment. In this paper, validation of Loci/GGFS crater width and depth predictions in ambient near-lunar conditions are provided using experimental data collected at MSFC during the Physics-Focused Ground Test 1 (PFGT-1) campaign in 2022. To observe sensitivity to soil models, simulations were conducted with both monodisperse glass bead (MGB) and BP-1 lunar regolith simulant soil models in Loci/GGFS. Crater depth and width comparisons are made with PFGT-1 Run 56, which used BP-1 soil. The Loci/GGFS BP-1 soil model performed slightly better with a mean predicted crater depth within 10% of the experiment. Both soil models predicted crater width within 10%. Cratering occurred more quickly with the MGB soil model. Mesh and spatial order sensitivity are also examined for the BP-1 soil model.

Validation↗

Tropical Ocean Global Atmosphere (TOGA) Meteorological and Oceanographic Data Sets for 1985 and 1986

The Tropical Ocean Global Atmosphere (TOGA) Program is a component of the World Meteorological Organization (WMO)/International Council of Scientific Unions (ICSU) World Climate Research Program (WCRP). One of the objectives of TOGA, which began in 1985, is to determine the limits of predictability of monthly mean sea surface temperature variations in tropical regions. The TOGA program created a raison d'etre for an explosive growth of the tropical ocean observing system and a substantial improvement in numerical simulations from atmospheric and oceanic general circulation models. Institutions located throughout the world are involved in the TOGA-distributed active data archive system. The diverse TOGA data sets for 1985 and 1986, including results from general circulation models, are included on a CD-ROM. Variables on the CD-ROM are barometric pressure, surface air temperature, dewpoint temperature Cartesian components of surface wind, surface sensible and latent heat fluxes,Cartesian components of surface wind stress and of an index of surface wind stress, sea level, sea surface temperature, and depth profiles of temperature and current in the upper ocean. Some data sets are global in extent, some are regional and cover portions of an ocean basin. Data on the CD-ROM can be extracted with an Apple Macintosh or an IBM PC.

Halpern, D.↗

Implications of high-spatial-resolution thermal infrared (Termoskan) data for Mars landing site selection

Thermal infrared observations of Mars from spacecraft provide physical information about the upper thermal skin depth of the surface, which is on the order of a few centimeters in depth and thus very significant for lander site selection. The Termoskan instrument onboard the Soviet Phobos '88 spacecraft acquired the highest spatial-resolution thermal infrared data obtained for Mars, ranging in resolution from 300 m to 3 km per pixel. It simultaneously obtained broadband reflected solar flux data. Although the 6 deg N - 30 deg S Termoskan coverage only slightly overlaps the nominal Mars Pathfinder target range, the implications of Termoskan data for that overlap region and the extrapolations that can be made to other regions give important clues for optimal landing site selection.

Betts, Bruce H.↗

Subsurface Temperature, Moisture, Thermal Conductivity and Heat Flux, Area A, B, C, and D, Utqiagvik (Barrow), Alaska, Ongoing from 2012

Subsurface soil monitoring stations are located across each of the four NGEE Arctic, Utqiagvik (Barrow), Alaska, intensive study Areas: (A) low-centered polygons with well-defined troughs, (B) high-centered polygons, (C) flat-centered or transitional polygons, and (D) inundated low-centered polygons with no troughs. Two data products are being generated - both with high temporal (hourly) and spatial resolution (within polygon and with depth): (1) soil profile temperature measurements and (2) active layer and permafrost temperature and moisture measurements from instrumented soil pits. This data package currently contains seventy-nine .csv files, one .xlsx (with complimentary .csv file), and one .pdf. Data from September 2012 through August 2022 are available as annual files. Additional data will be added annually after quality checks have been completed. These data are being collected to better understand the thermal dynamics of the active layer and permafrost. Together this information serves as calibration and evaluation data for fine-scale simulations of polygons. 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).

54 ENVIRONMENTAL SCIENCES↗

Voyager disk-integrated photometry of Triton

Hapke's (1981) photometric model has been combined with a plane-parallel thin atmospheric haze model to describe Voyager whole-disk observations of Triton, in the violet, blue, and green wavelength bands, in order to obtain estimates of Triton's geometric albedo, phase integral, and Bond albedo. Phase angle coverage in these filters ranging from about 12 to 159 deg was obtained by combining narrow- and wide-angle camera images. An upturn in the data at the highest phase angles observed can be explained by including scattering in a thin atmospheric haze layer with optical depths systematically decreasing with wavelength from about 0.06 in the violet to 0.03 for the green filter data.

Hillier, J.↗

Lidar measurements of boundary layers, aerosol scattering and clouds during project FIFE

A detailed account of progress achieved under this grant funding is contained in five journal papers. The titles of these papers are: The calculation of area-averaged vertical profiles of the horizontal wind velocity using volume imaging lidar data; Volume imaging lidar observation of the convective structure surrounding the flight path of an instrumented aircraft; Convective boundary layer mean depths, cloud base altitudes, cloud top altitudes, cloud coverages, and cloud shadows obtained from Volume Imaging Lidar data; An accuracy analysis of the wind profiles calculated from Volume Imaging Lidar data; and Calculation of divergence and vertical motion from volume-imaging lidar data. Copies of these papers form the body of this report.

Eloranta, Edwin W.↗

Snow Distribution Patterns Revisited: A Physics-Based and Machine Learning Hybrid Approach to Snow Distribution Mapping in the Sub-Arctic

Snowpack distribution in Arctic and alpine landscapes often occurs in repeating, year-to-year patterns due to local topographic, weather, and vegetation characteristics. Previous studies have suggested that with years of observational data, these snow distribution patterns can be statistically integrated into a snow process modeling workflow. Recent advances in snow hydrology and machine learning (ML) have increased our ability to predict snowpack distribution using in-situ observations, remote sensing data sets, and simple landscape characteristics that can be easily obtained for most environments. Here, we propose a hybrid approach to couple a ML snow distribution pattern (MLSDP) map with a physics-based, snow process model. We trained a random forest ML algorithm on tens of thousands of snow survey observations from a subarctic study area on the Seward Peninsula, Alaska, collected during peak snow water equivalent (SWE). We validated hybrid model outputs using in-situ snow depth and SWE observations, as well as a light detection and ranging data set and a distributed temperature profiling sensor data set. When the hybrid results were compared with the physics-based method, the hybrid method more accurately depicted the spatial patterns of the snowpack, areas of drifting snow, and years when no in-situ observations were used in the random forest ML training data set. The hybrid method also showed improvements in root mean squared error at 61% of locations where time-series estimations of snow depth were observed. These results can be applied to any physics-based model to improve the snow distribution patterning to reflect observed conditions in high latitude and high elevation cold region environments.

54 ENVIRONMENTAL SCIENCES↗

Soil and groundwater environmental sensor data, Wax Lake Delta, Louisiana, March 2023 - March 2024

This study evaluates how environmental parameters that integrate biogeochemical processes vary with water table fluctuations in the freshwater Wax Lake Delta (WLD) in Louisiana, U.S.A. This data package contains seven *.csv files and one Excel file that compiles all the data from the individual .csv files. This dataset reports high frequency (15-min) observations of water level, soil redox potential, specific conductance, and pH made for one year along elevation transects located on the older, proximal (OT) and younger, distal (YT) ends of a deltaic island. Water depth relative to the ground surface (cm; HOBO U20L-04; error ± 0.4 cm), water pH and temperature (HOBO MX2501), and specific conductance and temperature (HOBO U24-001) sensors were installed in March 2023. Water depth was corrected for barometric pressure recorded by a separate logger secured to a platform above the highest water level. Soil redox probes (SWAP ORP-40-4-B) were also installed in March 2023. Each probe had four Pt sensors (2 mm width) placed at 10 cm, 20 cm, 30 cm, and 40 cm below the ground surface. Redox data were referenced to an external Ag0/AgCl (3M KCl) reference probe placed in saturated ground and recorded on CR1000X dataloggers (Campbell Scientific) powered by solar panels. A second reference probe was positioned near the primary reference probe for backup and data correction. The tops of the soil redox probes and soil moisture probes were flush with the soil surface so that sensors are reported at their indicated depths below ground surface. Here, we report data collected between 15 March 2023 to 15 March 2024 for all sensors, with some differences due to exact dates of sensor placement or data gaps associated with sensor malfunction. For example, water depth at OT4 was not recorded between March to November 2023. Data flags indicate whether a value is valid (1) or was excluded from data analysis in the associated manuscript (-1).

EARTH SCIENCE > LAND SURFACE > SOILS↗

Microbial Community Metabolic Modeling: A Community Data-Driven Network Reconstruction: COMMUNITY DATA-DRIVEN METABOLIC NETWORK MODELING

Metabolic network modeling of microbial communities provides an in-depth understanding of community-wide metabolic and regulatory processes. Compared to single organism analyses, community metabolic network modeling is more complex because it needs to account for interspecies interactions. To date, most approaches focus on reconstruction of high-quality individual networks so that, when combined, they can predict community behaviors as a result of interspecies interactions. However, this conventional method becomes ineffective for communities whose members are not well characterized and cannot be experimentally interrogated in isolation. Here, we tested a new approach that uses community-level data as a critical input for the network reconstruction process. This method focuses on directly predicting interspecies metabolic interactions in a community, when axenic information is insufficient. We validated our method through the case study of a bacterial photoautotroph-heterotroph consortium that was used to provide data needed for a community-level metabolic network reconstruction. Resulting simulations provided experimentally validated predictions of how a photoautotrophic cyanobacterium supports the growth of an obligate heterotrophic species by providing organic carbon and nitrogen sources.

compartmentalized metabolic network↗

Dare Mighty Things

Winston Churchill once said, “To improve is to change; to be perfect is to change often.” JPL’s Property Accountability objective is to provide superior services related to property accountability, reutilization, and disposition. JPL must report property reutilized and disposed through either sales, donation, or scrap. Implementing the JPL designed and built Property Information Reporting System (PIRS), opened a new perspective with our data and what are we reporting to NASA. With this new depth of visibility, we self-implemented as assessment to validate the fiduciary and stewardship responsibilities of what is being reporting to NASA. JPL strives to perform at a level beyond the basic primary expectations of our NASA requirements. JPL developed the PIRS - which rolled out in 2018 - to simplify the delivery of JPL’s Personal Property and Equipment (PP&E) reporting: PIRS takes the raw data from Oracle to generate the NASA Form (NF) 1018 and Contractor Held Asset Reporting System (CHATS) Reports, and to meet the requirements established for AS9100 compliance. The data presented from PIRS is exportable and used to analyze our property records. JPL is also embarking on an effort to define the JPL of the future with “Enterprise 2.0”. Enterprise 2.0 is the driving force of JPL’s strategy to digitally transform the Laboratory from the current legacy processes and systems to a modern, integrated, information-driven highway of business transactions. This should include automation of routine processes and data-related tasks, integration of data and systems, advanced search and analytics, and improved information sharing and collaboration. As JPL’s evolving landscape of digital technologies advance to meet the future, JPL Property Accountability is making strides to surpass these expectations. For every cause there is an effect: a revealing of something that aids in the further development and transformation of JPL best business practices.

Sucher, Jay M↗

Non-Tidal Non-Seasonal Oceanic Mass Redistribution Estimated from the TOPEX/Poseidon Observation

Topex/Poseidon altimetry data are used to estimate the non-tidal mass redistribution as a function of space-time. The goal is to study the contribution of ocean circulations in the geodynamic effects including Earth's rotational and gravitational variations. We examine the non-seasonal anomalies at monthly sampling rate over the T/P span of eight years, concentrating especially on interannual variabilities. Since the sea-surface height data obtained from altimetry is the combined effect of steric change (primarily thermal effect) and the mass flux, and because the former has zero contribution to the geodynamic effects, one needs to do a so-called steric correction by removing from the altimetry data the steric contributions. We achieve it using multiyear monthly sea-surface temperature maps, together with monthly "climatology" for mixed-layer depth maps for mean-months of the year. We analyze both sets of the altimetry data and the steric correction data using the empirical orthogonal function/principal component analysis (in which we take care of issues associated with the area-weighting and non-zero mean), and examine the most important modes, either globally or regionally. In particular, the ENSO in the tropical Pacific/Indian Oceans exhibits the most prominent pattern. The net mass transport after the steric correction can then be compared with: (1) ocean general circulation model outputs for the same period of time (such as from POCM-4B); (2) non-atmospheric Earth rotation variations obtained from space geodesy data and atmospheric angular momentum data; (3) non-atmospheric low-degree gravitational variations from satellite-laser-ranging observations and global atmospheric data.

Chao, Benjamin F.↗

Final Scientific/Technical Report of Gas Hydrate Dynamics on the Alaskan Beaufort Continental Slope: Modeling and Field Characterization

Methane Hydrates, a solid form of methane and water, exist at high pressures and low temperatures, occurs on every continental margin on Earth, represents one of the largest reservoirs of carbon on the planet, and, if destabilized, may play an important role in both slope stability and climate change. For decades, researchers have studied methane hydrates with the hope of determining if methane hydrates are destabilizing, and if so, how this destabilization might impact slope stability and ocean/atmosphere carbon budgets. In the past ~5 years, it has become well established that the upper “feather-edge” of methane hydrate stability (intermediate water depths of ~200-500 meters below sea level) represents an important frontier for methane hydrates stability research, as this zone is most susceptible to destabilization due to minor fluctuations in ocean temperature in space and time. The Arctic Ocean—one of the fastest warming regions on Earth—is perhaps the best place to study possible changes to methane hydrate stability due to ocean warming. To address the stability of methane hydrates at intermediate ocean depths, Southern Methodist University in partnership with Oregon State University and The United State Geological Survey at Woods Hole began investigating methane hydrate stability in intermediate water depths below both the US Beaufort Sea and the Atlantic Margin, from 2012-2017. The work was funded by the Department of Energy’s (DOE) National Energy Technology Laboratory (NETL). The key goal of the SMU component of this study was to collect the first ever heat flow data in the Beaufort Sea and compare measured shallow (probe-based1) heat flow values with deeper (BSR-derived2) heat flow values, and from this, determine whether hydrates were in thermal equilibrium. In September 2016, SMU/OSU collected the first ever heat flow measurements in the US Beaufort Sea. Despite poor weather and rough seas, the cruise was a success, with 116 heat flow measurements acquired across the margin, spanning 4 transects separated by more than 400 km. Useable heat flow data exists for 97% (113) of probe heat flow measurements, revealing a clear picture of regional heat flow across the basin. During the past 8 months since the cruise, SMU researchers have processed the heat flow and thermal conductivity measurements and compared results to deeper heat flow estimates obtained from seismic data. The analysis reveals clear, consistent trends: All probe heat flow measurements in depths greater than 800 mbsl are consistent with BSR-derived values; heat flow measurements obtained in water depths between ~250-750 mbsl are systematically lower than those estimated from BSRs; and heat flow estimates in water depths shallower than ~250 mbsl are systematically warmer than deeper estimates. The consistency between shallow (probe) and deep (BSR) heat flow measurements at depths greater than ~750 m where ocean temperature changes are minimal supports the premise that the hydrates consist primarily of methane and represent a valuable tool for estimating heat flow. The anomalous cooling trend observed in the upper 250 m is consistent with expected seasonal effects observed in shallow ocean buoy measurements in the arctic, when cold, less dense melting sea ice cools the upper 200 m of the ocean during the summer as ice melting occurs. The discrepancy in heat flow at intermediate water depths is best explained via recent intermediate ocean temperature warming, where long-term (annual or longer) warming intermediate ocean bottom waters result in an anomalously low heat flow in shallow heat flow measurements. Using the characteristic 1D time-length scale for diffusion, we estimate that ocean temperature warming began no later than ~1200 years ago but arguably much more recently as results are limited by seismic resolution. More importantly, our analysis indicates methane hydrate is destabilizing not only in the upper feather edge (200-500 mbsl) but at depths as great as 750 mbsl. The intermediate ocean warming rate supports previous studies suggesting geologically rapid warming (>0.1 deg C/decade) at intermediate ocean depths in the Beaufort Sea. Assuming no further changes or additional warming, our analysis indicates methane hydrates will destabilize at seafloor depths shallower than 750 mbsl in the Beaufort Sea within the next ~3000 years. 1 Probe outfitted with sensors inserted into the seafloor sediment 2 Bottom-simulating reflector (BSR) seismic data indicates presence of hydrate deposits

58 GEOSCIENCES↗