Search NASA⌕ Search

SEARCH · Search NASA

Results for “Icing conditions”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Increasing large wildfires over the western United States linked to diminishing sea ice in the Arctic

Abstract The compound nature of large wildfires in combination with complex physical and biophysical processes affecting variations in hydroclimate and fuel conditions makes it difficult to directly connect wildfire changes over fire-prone regions like the western United States (U.S.) with anthropogenic climate change. Here we show that increasing large wildfires during autumn over the western U.S. are fueled by more fire-favorable weather associated with declines in Arctic sea ice during preceding months on both interannual and interdecadal time scales. Our analysis (based on observations, climate model sensitivity experiments, and a multi-model ensemble of climate simulations) demonstrates and explains the Arctic-driven teleconnection through regional circulation changes with the poleward-shifted polar jet stream and enhanced fire-favorable surface weather conditions. The fire weather changes driven by declining Arctic sea ice during the past four decades are of similar magnitude to other leading modes of climate variability such as the El Niño-Southern Oscillation that also influence fire weather in the western U.S.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Physicochemical characterization of free troposphere and marine boundary layer ice-nucleating particles collected by aircraft in the eastern North Atlantic

Abstract. Atmospheric ice nucleation impacts the hydrological cycle and climate by modifying the radiative properties of clouds. To improve our predictive understanding of ice formation, ambient ice-nucleating particles (INPs) need to be collected and characterized. Measurements of INPs at lower latitudes in a remote marine region are scarce. The Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign, in the region of the Azores islands, provided the opportunity to collect particles in the marine boundary layer (MBL) and free troposphere (FT) by aircraft during the campaign's summer and winter intensive operation period. The particle population in samples collected was examined by scanning transmission X-ray microscopy with near-edge X-ray absorption fine structure spectroscopy. The identified INPs were analyzed by scanning electron microscopy with energy-dispersive X-ray analysis. We observed differences in the particle population characteristics in terms of particle diversity, mixing state, and organic volume fraction between seasons, mostly due to dry intrusion events during winter, as well as between the sampling locations of the MBL and FT. These differences are also reflected in the temperature and humidity conditions under which water uptake, immersion freezing (IMF), and deposition ice nucleation (DIN) proceed. Identified INPs reflect typical particle types within the particle population on the samples and include sea salt, sea salt with sulfates, and mineral dust, all associated with organic matter, as well as carbonaceous particles. IMF and DIN kinetics are analyzed with respect to heterogeneous ice nucleation rate coefficients, Jhet, and ice nucleation active site density, ns, as a function of the water criterion Δaw. DIN is also analyzed in terms of contact angles following classical nucleation theory. Derived MBL IMF kinetics agree with previous ACE-ENA ground-site INP measurements. FT particle samples show greater ice nucleation propensity compared to MBL particle samples. This study emphasizes that the types of INPs can vary seasonally and with altitude depending on sampling location, thereby showing different ice nucleation propensities, which is crucial information when representing mixed-phase cloud and cirrus cloud microphysics in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Partitioning uncertainty in projections of Arctic sea ice

Abstract Improved knowledge of the contributing sources of uncertainty in projections of Arctic sea ice over the 21st century is essential for evaluating impacts of a changing Arctic environment. Here, we consider the role of internal variability, model structure and emissions scenario in projections of Arctic sea-ice area (SIA) by using six single model initial-condition large ensembles and a suite of models participating in Phase 5 of the Coupled Model Intercomparison Project. For projections of September Arctic SIA change, internal variability accounts for as much as 40%–60% of the total uncertainty in the next decade, while emissions scenario dominates uncertainty toward the end of the century. Model structure accounts for 60%–70% of the total uncertainty by mid-century and declines to 30% at the end of the 21st century in the summer months. For projections of wintertime Arctic SIA change, internal variability contributes as much as 50%–60% of the total uncertainty in the next decade and impacts total uncertainty at longer lead times when compared to the summertime. In winter, there exists a considerable scenario dependence of model uncertainty with relatively larger model uncertainty under strong forcing compared to weak forcing. At regional scales, the contribution of internal variability can vary widely and strongly depends on the calendar month and region. For wintertime SIA change in the Greenland-Iceland-Norwegian and Barents Seas, internal variability contributes 60%–70% to the total uncertainty over the coming decades and remains important much longer than in other regions. We further find that the relative contribution of internal variability to total uncertainty is state-dependent and increases as sea ice volume declines. These results demonstrate that internal variability is a significant source of uncertainty in projections of Arctic sea ice.

Bonan, David B. (ORCID:0000000338676009)↗
DOE OSTI2022

Sea Ice Rheology Experiment (SIREx): 2. Evaluating Linear Kinematic Features in High-Resolution Sea Ice Simulations

Simulating sea ice drift and deformation in the Arctic Ocean is still a challenge because of the multiscale interaction of sea ice floes that compose the Arctic Sea ice cover. The Sea Ice Rheology Experiment (SIREx) is a model intercomparison project of the Forum of Arctic Modeling and Observational Synthesis (FAMOS). In SIREx, skill metrics are designed to evaluate different recently suggested approaches for modeling linear kinematic features (LKFs) to provide guidance for modeling small-scale deformation. These LKFs are narrow bands of localized deformation that can be observed in satellite images and also form in high resolution sea ice simulations. In this contribution, spatial and temporal properties of LKFs are assessed in 36 simulations of state-of-the-art sea ice models and compared to deformation features derived from the RADARSAT Geophysical Processor System. All simulations produce LKFs, but only very few models realistically simulate at least some statistics of LKF properties such as densities, lengths, or growth rates. All SIREx models overestimate the angle of fracture between conjugate pairs of LKFs and LKF lifetimes pointing to inaccurate model physics. The temporal and spatial resolution of a simulation and the spatial resolution of atmospheric boundary condition affect simulated LKFs as much as the model's sea ice rheology and numerics. Only in very high resolution simulations (≤2 km) the concentration and thickness anomalies along LKFs are large enough to affect air-ice-ocean interaction processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Riverbank temperatures on the Selawik River, Alaska 2010-2012, and Koyukuk River, Alaska June to July 2018.

The data package includes temperature measurements from riverbanks along the Selawik and Koyukuk Rivers in northwest Alaska. The Selawik River data also includes meteorological data collected locally and the Koyukuk River includes bulk density and ice content data from the riverbanks. The Selawik River data was collected between 2010 and 2013 with meteorological data collected for the entire time period and riverbank temperatures collected in shorter time intervals between 2010 to 2012. All of the Selawik River data were collected a single river bend. The Koyukuk River temperature data was collected from late June to early July 2018 at a total of five locations on three riverbanks. Ice content measurements on the Koyukuk were made at five riverbank locations. The data collection on both rivers was done in support of studies on the influence of permafrost on riverbank erosion. The temperature data was collected to monitor the thermal conditions of the river banks both seasonally and during periods of erosion. The ice content of the banks was measured to better understand how bank erosion rates are influenced by ice in the frozen riverbanks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modeling impacts of ice-nucleating particles from marine aerosols on mixed-phase orographic clouds during 2015 ACAPEX field campaign

Abstract. A large fraction of annual precipitation over the western United States comes from wintertime orographic clouds associated with atmospheric rivers (ARs). Transported African and Asian dust and marine aerosols from the Pacific Ocean may act as ice-nucleating particles (INPs) to affect cloud and precipitation properties over the region. Here we explored the effects of INPs from marine aerosols on orographic mixed-phase clouds and precipitation at different AR stages for an AR event observed during the 2015 ACAPEX field campaign under low dust (<0.02 cm−3) conditions. Simulations were conducted using the chemistry version of the Weather Research and Forecasting Model coupled with the spectral-bin microphysics at 1 km grid spacing, with ice nucleation connected with dust and marine aerosols. By comparing against airborne and ground-based observations, accounting for marine INP effects improves the simulation of AR-precipitation. The marine INPs enhance the formation of ice and snow, leading to less shallow warm clouds but more mixed-phase and deep clouds, as well as to a large spillover effect of precipitation after AR landfall. The responses of cloud and precipitation to marine INPs vary with the AR stages, with more significant effects before AR landfall and post-AR than after AR landfall, mainly because the moisture and temperature conditions change with the AR evolution. This work suggests weather and climate models need to consider the impacts of marine INPs since their contribution is notable under low dust conditions despite the much lower relative ice nucleation efficiency of marine INPs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Liquid-liquid phase separation in supercooled water from ultrafast heating of low-density amorphous ice

Recent experiments continue to find evidence for a liquid-liquid phase transition (LLPT) in supercooled water, which would unify our understanding of the anomalous properties of liquid water and amorphous ice. These experiments are challenging because the proposed LLPT occurs under extreme metastable conditions where the liquid freezes to a crystal on a very short time scale. Here, we analyze models for the LLPT to show that coexistence of distinct high-density and low-density liquid phases may be observed by subjecting low-density amorphous (LDA) ice to ultrafast heating. We then describe experiments in which we heat LDA ice to near the predicted critical point of the LLPT by an ultrafast infrared laser pulse, following which we measure the structure factor using femtosecond x-ray laser pulses. Consistent with our predictions, we observe a LLPT occurring on a time scale <100 ns and widely separated from ice formation, which begins at times >1 μs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

2D CFD of lean premixed hydrogen–air flame quenching under locomotive engine conditions

Hydrogen (H 2 ) is a promising fuel for reducing emissions in heavy-duty internal combustion engines (ICEs), but its low quenching distance increases the risk of flame propagation into narrow crevice regions, such as piston-liner gaps. This work uses detailed CFD simulations with finite-rate chemistry to investigate premixed H 2 –air flame quenching in a two-dimensional (2D) region consistent with the piston-liner gap of a diesel ICE. Model accuracy was assessed by comparison with experimentally measured quenching distances in an annular stepwise diverging tube (ASDT). A parametric study was conducted to assess the influence of crevice width (0.05–1.18 mm), pressure (50–150 bar), unburned gas temperature (431–573 K), and equivalence ratio ( Φ= 0.3–0.6). Results show that the critical Péclet number for flame survival is no greater than 3.25, consistent with prior literature, even under ultra-lean and high-pressure conditions (Φ ≤ 0.3, > 50 bar). Additionally, flames with Péclet numbers exceeding 6.15 exhibited front wrinkling, suggesting the onset of velocity-driven instabilities and enhanced flame robustness. These findings help define thresholds for flame quenching in confined geometries and support the safe design of H 2 fueled ICEs.

08 HYDROGEN↗
DOE OSTI2021

Facilitating Large‐Scale Snow Shedding from In‐Field Solar Arrays using Icephobic Surfaces with Low‐Interfacial Toughness

Abstract Large‐scale accrual of snow and ice on solar arrays in northern latitudes can cause significant power generation losses during winter. Depending on environmental conditions, snow can encompass a wide range in physical characteristics from dry snow (modulus ≈100 kPa and density ≈0.1 g cm −3 ) to bulk ice (modulus ≈8 GPa and density ≈0.9 g cm −3 ). This variation in snow morphology has made the development of a passive, broad‐spectrum, snow and ice‐shedding surface challenging. Here, the authors develop one of the first surfaces that simultaneously possesses both low‐interfacial strength ( τ˄ ice < 50 kPa) and toughness (Γ ice < 0.5 J m −2 ) with ice. These surfaces, fabricated via the addition of mobile polymer chains/oils to a thin polymeric coating, require extremely low detachment forces for ice, enabling its passive shedding at virtually any accretion length scale. Preliminary evidence that the new surfaces can shed different forms of snow and ice from field‐deployed solar arrays, over a range of subzero temperatures for several weeks, leading to significant increases in power generation is provided. The optically transparent surfaces are easily scalable and can be widely deployed by the solar industry in areas that see persistent snow. Other applications include automotive windshields, LIDAR covers for autonomous vehicles, and cold climate optical sensors.

Dhyani, Abhishek↗
DOE OSTI2022

Salinity tolerance mechanisms of an Arctic Pelagophyte using comparative transcriptomic and gene expression analysis

Abstract Little is known at the transcriptional level about microbial eukaryotic adaptations to short-term salinity change. Arctic microalgae are exposed to low salinity due to sea-ice melt and higher salinity with brine channel formation during freeze-up. Here, we investigate the transcriptional response of an ice-associated microalgae over salinities from 45 to 8. Our results show a bracketed response of differential gene expression when the cultures were exposed to progressively decreasing salinity. Key genes associated with salinity changes were involved in specific metabolic pathways, transcription factors and regulators, protein kinases, carbohydrate active enzymes, and inorganic ion transporters. The pelagophyte seemed to use a strategy involving overexpression of Na + -H + antiporters and Na + -Pi symporters as salinity decreases, but the K + channel complex at higher salinities. Specific adaptation to cold saline arctic conditions was seen with differential expression of several antifreeze proteins, an ice-binding protein and an acyl-esterase involved in cold adaptation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Supporting data for "Frozen flume experiments indicate rapid permafrost riverbank erosion depends on bank roughness"

This dataset contains supporting files detailing five frozen flume experiments conducted at the Caltech Earth Surface Dynamics Laboratory to investigate rates of ablation-limited permafrost riverbank erosion under controlled conditions. Water flowed past a bank of saturated, frozen sand and ice and gradually eroded the bank by thawing pore ice and immediately entraining sand and washing it downstream. Experiments were scaled for flow hydraulics and heat transfer allowing comparisons between our results and natural permafrost riverbanks. For each experiment, we measured the initial and final sand bank topography using a Keyence laser scanner, water surface slope at 3-min intervals throughout the experiment using a Massa sonar scanner, bank erosion using 10-sec overhead timelapse imagery taken by an overhead camera, water and bank temperature using thermistors frozen into the sand bank and sampling at 2 Hz, and water discharge using an in-line flow meter. We include calibration data for the carriage (engineered by the Saint Anthony Falls Laboratory) used to make sonar and laser topography measurements. We also include calibration data for temperature sensors, water discharge measurements, and images of a regular grid placed in the flume to align overhead camera images with the carriage datum. Grain size analysis for the channel bed (gravel) was produced using a pebble count and bank sand was measured using a Camsizer X2. In addition to the five frozen experiments, we include sonar scans of water surface slope and Keyence scans of bed and bank topography for calibration experiments ran with an immobile gravel bank and bed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Technical descriptions of the experimental dynamical downscaling simulations over North America by the CAM–MPAS variable-resolution model

Abstract. Comprehensive assessment of climate datasets is important for communicating model projections and associated uncertainties to stakeholders. Uncertainties can arise not only from assumptions and biases within the model but also from external factors such as computational constraint and data processing. To understand sources of uncertainties in global variable-resolution (VR) dynamical downscaling, we produced a regional climate dataset using the Model for Prediction Across Scales (MPAS; dynamical core version 4.0) coupled to the Community Atmosphere Model (CAM; version 5.4), which we refer to as CAM–MPAS hereafter. This document provides technical details of the model configuration, simulations, computational requirements, post-processing, and data archive of the experimental CAM–MPAS downscaling data. The CAM–MPAS model is configured with VR meshes featuring higher resolutions over North America as well as quasi-uniform-resolution meshes across the globe. The dataset includes multiple uniform- (240 and 120 km) and variable-resolution (50–200, 25–100, and 12–46 km) simulations for both the present-day (1990–2010) and future (2080–2100) periods, closely following the protocol of the North American Coordinated Regional Climate Downscaling Experiment. A deviation from the protocol is the pseudo-warming experiment for the future period, using the ocean boundary conditions produced by adding the sea surface temperature and sea-ice changes from the low-resolution version of the Max Planck Institute Earth System Model (MPI-ESM-LR) in the Coupled Model Intercomparison Project Phase 5 to the present-day ocean state from a reanalysis product. Some unique aspects of global VR models are evaluated to provide background knowledge to data users and to explore good practices for modelers who use VR models for regional downscaling. In the coarse-resolution domain, strong resolution sensitivity of the hydrological cycles exists over the tropics but does not appear to affect the midlatitude circulations in the Northern Hemisphere, including the downscaling target of North America. The pseudo-warming experiment leads to similar responses of large-scale circulations to the imposed radiative and boundary forcings in the CAM–MPAS and MPI-ESM-LR models, but their climatological states in the historical period differ over various regions, including North America. Such differences are carried to the future period, suggesting the importance of the base state climatology. Within the refined domain, precipitation statistics improve with higher resolutions, and such statistical inference is verified to be negligibly influenced by horizontal remapping during post-processing. Limited (≈50 % slower) throughput of the current code is found on a recent many-core/wide-vector high-performance computing system, which limits the lengths of the 12–46 km simulations and indirectly affects sampling uncertainty. Our experience shows that global and technical aspects of the VR downscaling framework require further investigations to reduce uncertainties for regional climate projection.

58 GEOSCIENCES↗
DOE OSTI2025

Arctic sea-ice loss drives a strong regional atmospheric response over the North Pacific and North Atlantic on decadal scales

Previous studies have suggested that Arctic sea-ice loss can have a profound influence on atmospheric circulation far away from the Arctic. However, there is little scientific consensus on the features of these remote responses, with the opposite impacts reported. Here we present a multi-model analysis of the decadal climate response to Arctic sea-ice loss using state-of-the-art energy conserving methodologies to isolate the impacts of sea-ice decline. We observe weakening of the Aleutian Low and development of a geopotential ridge in the North Pacific, associated with drier winter conditions over the southwest United States. Over the Atlantic, a negative NAO-like response drives wetter winter conditions across the western Mediterranean. These decadal-scale impacts substantially differ from reported centennial-scale responses to Arctic sea-ice loss simulated using non-energy conserving methodologies. Factors such as the timescale of the response and methodologies used to isolate the impacts of disappearing sea-ice cover should be carefully considered when consolidating scientific understanding on the future impacts of changing Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Micro-spectroscopic and freezing characterization of ice-nucleating particles collected in the marine boundary layer in the eastern North Atlantic

Abstract. Formation of atmospheric ice plays a crucial role in the microphysical evolution of mixed-phase and cirrus clouds and thus climate. How aerosol particles impact ice crystal formation by acting as ice-nucleating particles (INPs) is a subject of intense research activities. To improve understanding of atmospheric INPs, we examined daytime and nighttime particles collected during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) field campaign conducted in summer 2017. Collected particles, representative of a remote marine environment, were investigated for their propensity to serve as INPs in the immersion freezing (IMF) and deposition ice nucleation (DIN) modes. The particle population was characterized by chemical imaging techniques such as computer-controlled scanning electron microscopy with energy-dispersive X-ray analysis (CCSEM/EDX) and scanning transmission X-ray microscopy with near-edge X-ray absorption fine-structure spectroscopy (STXM/NEXAFS). Four major particle-type classes were identified where internally mixed inorganic–organic particles make up the majority of the analyzed particles. Following ice nucleation experiments, individual INPs were identified and characterized by SEM/EDX. The identified INP types belong to the major particle-type classes consisting of fresh sea salt with organics or processed sea salt containing dust and sulfur with organics. Ice nucleation experiments show IMF events at temperatures as low as 231 K, including the subsaturated regime. DIN events were observed at lower temperatures of 210 to 231 K. IMF and DIN observations were analyzed with regard to activated INP fraction, ice-nucleation active site (INAS) densities, and a water activity-based immersion freezing model (ABIFM) yielding heterogeneous ice nucleation rate coefficients. Observed IMF and DIN events of ice formation and corresponding derived freezing rates demonstrate that the marine boundary layer aerosol particles can serve as INPs under typical mixed-phase and cirrus cloud conditions. The derived IMF and DIN parameterizations allow for implementation in cloud and climate models to evaluate predictive effects of atmospheric ice crystal formation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

NED_Asteroid_Energy_Deposition

In the event of a potentially catastrophic asteroid impact, with sufficient warning time, deploying a nuclear device remains a powerful option for planetary defense if a kinetic impactor proves insufficient. Predicting the effectiveness of a potential nuclear deflection or disruption mission depends on accurate multiphysics simulations of the device's x-ray energy deposition into the asteroid and the resulting material ablation. These simulations span many orders of magnitude, require a variety of different complex physics packages, and are computationally expensive. Having an efficient and accurate way of modeling this system is necessary for exploring a mission's sensitivity to the asteroid's range of physical properties. To expedite future simulations, we present a completed x-ray energy deposition model developed using the radiation-hydrodynamics code Kull which can be used to initiate a nuclear mitigation mission hydrocode. The model spans a wide variety of possible mission initial conditions: four different asteroid-like materials (Silicon Dioxide, Forsterite, Iron, and Ice), two different source spectra (1 and 2 keV blackbodies), and then a broad range of radiation fluences (0.0001 to 1 kt per square meter), source durations (10 to 100 ns), and asteroid porosities (0 to 80 percent). Using blowoff momentum as the primary metric, the modelinitiated simulation results match the full radiation-hydrodynamics results to within 10 percent. Please reference the journal article: Burkey et al., X-Ray Energy Deposition Model for Simulating Asteroid Response to a Nuclear Planetary Defense Mitigation Mission, Planetary Science Journal, (2023) for more information.

Managan, RobertA↗
DOE OSTI2025

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: II. Time evolution

Dedicated DIII-D experiments studied instabilities below the ion cyclotron frequency f ci and at cyclotron harmonics by injecting thirteen different neutral beam populations into plasmas with different thermal compositions of H and D at five different values of toroidal field. The resulting database of 2529 conditions is described in Part 1: Stability. This paper focuses on the temporal evolution. For some beam geometries, the ion cyclotron emission (ICE) is much stronger at the beginning of injection, suggesting the importance of ‘bump-on-tail’ energy drive. Sawtooth crashes cause the emission to cease for some beam types and emission bands, but the opposite effect occurs for others. When unstable, for virtually all conditions and frequency bands, the emission bursts in a cycle reminiscent of predator-prey oscillations. Evidence for nonlinear interaction between ICE harmonics is observed. Strongly unstable modes typically appear within 1 ms from the start of beam injection and often decay within 0.1 ms when the beam turns off, but both the initial rise and final decay of the mode power exhibit a wide variety of behavior.

electromagnetic ion cyclotron instabilities↗
DOE OSTI2020

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗