Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ice physics”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Small field campaign: aerosol – ice formation closure pilot study. Final Report

Prediction of atmospheric ice formation from aerosol particles by heterogeneous nucleation represents one of the grand challenges in atmospheric science. Our insufficient predictive understanding of primary ice formation is the reason that climate models typically do not include heterogeneous ice nucleation with subsequent effects on climate uncertainty. Mixed-phase clouds, where supercooled water droplets and ice crystals coexist play globally an important role regulating climate. This is especially the case for the Arctic region that experiences the greatest warming due to climate change compared to other regions in the world. Immersion freezing initiated by ice-nucleating particles (INPs) in supercooled water droplets is recognized as the dominant primary ice formation pathway in mixed-phase cloud regimes. For this reason, it is crucial to evaluate our capability to predict immersion freezing for a given ambient aerosol population. The goal of this project is conducting a field-based pilot study at the U.S. DOE Atmospheric Radiation Measurement (ARM) user facility at Southern Great Plains (SGP) to evaluate our capability to predict the number concentration of aerosol particles that serve as INPs in the immersion freezing mode. Successful prediction of INP number concentrations is also termed “closure”. This field-observational approach represents a first-of-its kind attempt of an aerosol–ice formation closure study (AEROICESTUDY). Very few closure studies related to INPs have been conducted, and to our knowledge, none using robust size-resolved ambient aerosol composition measurements as a starting point. Achievement of aerosol–ice formation closure relies on our ability to characterize the ambient aerosol population with respect to particles size and composition and to determine INP number concentrations for specified freezing temperatures. This requires numerous online and offline instrumentation resulting in this pilot field campaign being a multi-institutional and community-collaborative effort. We chose the ARM SGP megasite for this first aerosol-ice formation closure pilot study due to its significant measurement capabilities available to obtain detailed physical characterization of the local aerosol population including size distribution, mass loading, and chemical composition of non-refractory aerosol particles. The overall objective of this project is to identify ice nucleation parameterizations that produce the most robust predictions of INP numbers and thus are best suited to be included in cloud and climate models. This objective includes the following goals for this field and laboratory-based project: i) What are the crucial aerosol physicochemical property measurements needed to accurately guide ice nucleation representations in models and long-term INP measurements? ii) What level of parameter details needs to be known to achieve aerosol–ice formation closure? iii) What are the leading causes for climate model bias in INP predictions? We found that the advances in our understanding of immersion freezing garnered over the last 20 years allowed us to yield partial and full closures of atmospheric immersion freezing from ambient aerosol particles. When the aerosol population is physicochemically complex and parameterizations for representative INP types are not yet available, we still struggle to accurately predict INP number concentrations. This project clearly demonstrates that with more laboratory and field measurements that are accompanied by particle composition analysis, the necessary datasets to achieve aerosol–ice formation closure for various locations will emerge, thus providing a robust foundation for guiding the representation of INPs in cloud and climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Epoch of Giant Planet Migration Planet Search Program. I. Near-infrared Radial Velocity Jitter of Young Sun-like Stars

We present early results from the Epoch of Giant Planet Migration program, a precise radial velocity (RV) survey of more than 100 intermediate-age (∼20–200 Myr) G and K dwarfs with the Habitable Zone Planet Finder spectrograph (HPF) at McDonald Observatory’s Hobby–Eberly Telescope. The goals of this program are to determine the timescale and dominant physical mechanism of giant planet migration interior to the water ice line of Sun-like stars. Here, we summarize results from the first 14 months of this program, with a focus on our custom RV pipeline for HPF, a measurement of the intrinsic near-infrared RV activity of young Solar analogs, and modeling the underlying population-level distribution of stellar jitter. We demonstrate on-sky stability at the sub-2 m s{sup −1} level for the K2 standard HD 3765 using a least-squares matching method to extract precise RVs. Based on a subsample of 29 stars with at least three RV measurements from our program, we find a median rms level of 34 m s{sup −1}. This is nearly a factor of 2 lower than the median rms level in the optical of 60 m s{sup −1} for a comparison sample with similar ages and spectral types as our targets. The observed near-infrared jitter measurements for this subsample are well reproduced with a log-normal parent distribution with μ = 4.15 and σ = 1.02. Finally, by compiling rms values from previous planet search programs, we show that near-infrared jitter for G and K dwarfs generally decays with age in a similar fashion to optical wavelengths, albeit with a shallower slope and lower overall values for ages ≲1 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A Formulation of the Many-Body Expansion (MBE) for Periodic Systems: Application to Several Ice Phases

A new formulation of the many-body expansion (MBE) for periodic systems, based on their inherent translational symmetry, is introduced and applied to 7 ice polymorphs. This new formulation is built via a hierarchical procedure connecting gas phase unit cells over finite supercells to infinite solids. For finite supercells, it is demonstrated that this method successfully recovers the energetics while reducing the scaling of the calculation of the many-body terms by a factor of ~N, where N is the system size. Furthermore, the proposed framework delineates a straightforward way to obtain properties of macroscopic systems in the limit of an infinite cell. For periodic systems, the success of this approach is demonstrated by showing that the lattice energies computed (up to the 4-th order in the MBE) reproduce the lattice energies obtained using periodic boundary conditions with an Ewald summation for 7 polymorphs of ice (Ih, II, VIII, IX, XIII, XIV, XV). This development makes it possible to quantify, for the first time, the many-body contributions to the lattice energy of various ice polymorphs. The many-body (three-body and higher) interactions were found to vary significantly among the 7 ice phases, amounting between 7-24% of the total lattice energies. So, this development opens the door for obtaining insights into solid-state properties, while leveraging the computational benefits of the MBE.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A spatial evaluation of Arctic sea ice and regional limitations in CMIP6 historical simulations

The Arctic sea ice response to a warming climate is assessed in a subset of models participating in Phase 6 of the Coupled Model Intercomparison Project (CMIP6), using several metrics in comparison with satellite observations and results from the Pan-Arctic Ice Ocean Modeling and Assimilation System and the Regional Arctic System Model. Our study examines the historical representation of sea ice extent, volume, and thickness using spatial analysis metrics, such as the integrated ice-edge error, Brier score, and Spatial Probability Score. We find that the CMIP6 multi-model mean captures the mean annual cycle and 1979-2014 sea ice trends remarkably well. However, individual models experience a wide range of uncertainty in the spatial distribution of sea ice when compared against satellite measurements and reanalysis data. Our metrics expose common and individual regional model biases, which sea ice temporal analyses alone do not capture. We identify large ice edge and ice thickness errors in Arctic sub-regions, implying possible model specific limitations in or lack of representation of some key physical processes. We postulate that many of them could be related to the oceanic forcing, especially in the marginal and shelf seas, where seasonal sea ice changes are not adequately simulated. We therefore conclude that an individual model’s ability to represent the observed/reanalysis spatial distribution still remains a challenge. We propose the spatial analysis metrics as useful tools to diagnose model limitations, narrow down possible processes affecting them, and guide future model improvements critical to the representation and projections of Arctic climate change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Dipolar Spin Ice Regime Proximate to an All-In-All-Out Néel Ground State in the Dipolar-Octupolar Pyrochlore Ce 2 Sn 2 O 7

The dipolar-octupolar (DO) pyrochlores, R 2 M 2 O 7 ( R = Ce , Sm , Nd ), are key players in the search for realizable novel quantum spin liquid (QSL) states as a large parameter space within the DO pyrochlore phase diagram is theorized to host QSL states of both dipolar and octupolar nature. New single crystals and powders of Ce 2 Sn 2 O 7 , synthesized by hydrothermal techniques, present an opportunity for a new characterization of the exchange parameters in Ce 2 Sn 2 O 7 using the near-neighbor X Y Z model Hamiltonian associated with DO pyrochlores. Utilizing quantum numerical linked cluster expansion fits to heat capacity and magnetic susceptibility measurements, and classical Monte Carlo calculations to the diffuse neutron diffraction of the new hydrothermally grown Ce 2 Sn 2 O 7 samples, we place Ce 2 Sn 2 O 7 ’s ground state within the ordered dipolar all-in-all-out (AIAO) Néel phase, with quantum Monte Carlo calculations showing a transition to long-range order at temperatures below those accessed experimentally. Indeed, our new neutron diffraction measurements on the hydrothermally grown Ce 2 Sn 2 O 7 powders show a broad signal at low scattering wave vectors, reminiscent of a spin ice, in striking contrast from previous powder neutron diffraction on samples grown from solid-state synthesis, which found diffuse scattering at high scattering wave vectors associated with magnetic and suggested an octupolar quantum spin ice state. We conclude that new hydrothermally grown Ce 2 Sn 2 O 7 samples host a finite-temperature proximate dipolar spin ice phase, above the expected transition to AIAO Néel order. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Quantifying the Resolution Sensitivity of the Kain–Fritsch Scheme Across the Gray Zone by Isolating Interactions: A TWP‐ICE Case Study

The resolution sensitivity of the Kain–Fritsch (KF) convection scheme and the role of interactions between the physics and dynamics within the gray zone (<10 km) were investigated using the Separate Physics and Dynamics Experiment (SPADE) framework. Two groups of experiments were conducted using the Weather Research and Forecasting (WRF) model via traditional (Tradition) runs and SPADE runs with resolutions of 1, 2, 4, and 8 km during the wet period of the Tropical Warm Pool–International Cloud Experiment (TWP‐ICE). Results show that the KF scheme simulates the weakened convective processes well as the resolution increases in both groups, and the changes in the convective variables with resolution in SPADE are smaller than in the Tradition group. This indicates the important effects of interactions between model components on convection parameterizations as the resolution changes. Additionally, the microphysics variables remain nearly unchanged with resolution in SPADE and weaken slightly in Tradition as the resolution decreases, suggesting the relatively weaker influences of model interactions for the resolved‐cloud parameterization. Therefore, the scale‐aware behavior of KF scheme is further strengthened in Tradition runs, primarily through inhibiting the strength of stratiform processes through physics–dynamics interactions and physical components.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Heterogeneous Ice Growth in Micron-Sized Water Droplets Due to Spontaneous Freezing

Understanding how ice nucleates and grows into larger crystals is of crucial importance for many research fields. The purpose of this study was to shed light on the phase and structure of ice once a nucleus is formed inside a metastable water droplet. Wide-angle X-ray scattering (WAXS) was performed on micron-sized droplets evaporatively cooled to temperatures where homogeneous nucleation occurs. We found that for our weak hits ice grows more cubic compared to the strong hits that are completely hexagonal. Due to efficient heat removal caused by evaporation, we propose that the cubicity of ice at the vicinity of the droplet’s surface is higher than for ice formed within the bulk of the droplet. Moreover, the Bragg peaks were classified based on their geometrical shapes and positions in reciprocal space, which showed that ice grows heterogeneously with a significant population of peaks indicative of truncation rods and crystal defects. Frequent occurrences of the (100) reflection with extended in-planar structure suggested that large planar ice crystals form at the droplet surface, then fracture into smaller domains to accommodate to the curvature of the droplets. Planar faulting due to misaligned domains would explain the increased cubicity close to the droplet surface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Statistical emulation of a perturbed basal melt ensemble of an ice sheet model to better quantify Antarctic sea level rise uncertainties

Abstract. Antarctic ice shelves are vulnerable to warming ocean temperatures, and some have already begun thinning in response to increased basal melt rates. Sea level is therefore expected to rise due to Antarctic contributions, but uncertainties in its amount and timing remain largely unquantified. In particular, there is substantial uncertainty in future basal melt rates arising from multi-model differences in thermal forcing and how melt rates depend on that thermal forcing. To facilitate uncertainty quantification in sea level rise projections, we build, validate, and demonstrate projections from a computationally efficient statistical emulator of a high-resolution (4 km) Antarctic ice sheet model, the Community Ice Sheet Model version 2.1. The emulator is trained to a large (500-member) ensemble of 200-year-long 4 km resolution transient ice sheet simulations, whereby regional basal melt rates are perturbed by idealized (yet physically informed) trajectories. The main advantage of our emulation approach is that by sampling a wide range of possible basal melt trajectories, the emulator can be used to (1) produce probabilistic sea level rise projections over much larger Monte Carlo ensembles than are possible by direct numerical simulation alone, thereby providing better statistical characterization of uncertainties, and (2) predict the simulated ice sheet response under differing assumptions about basal melt characteristics as new oceanographic studies are published, without having to run additional numerical ice sheet simulations. As a proof of concept, we propagate uncertainties about future basal melt rate trajectories, derived from regional ocean models, to generate probabilistic sea level rise estimates for 100 and 200 years into the future.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Optical Absorption Spectroscopy Probes Water Wire and Its Ordering in a Hydrogen-Bond Network

Water wires, quasi-one-dimensional chains composed of hydrogen-bonded (H-bonded) water molecules, play a fundamental role in numerous chemical, physical, and physiological processes. Yet direct experimental detection of water wires has been elusive so far. Based on advanced many-body theory that includes electron-hole interactions, we report that optical absorption spectroscopy can serve as a sensitive probe of water wires and their ordering. In both liquid and solid water, the main peak of the spectrum is discovered to be a charge-transfer exciton. In water, the charge-transfer exciton is strongly coupled to the H-bonding environment where the exciton is excited between H-bonded water molecules with a large spectral intensity. In regular ice, the spectral weight of the charge-transfer exciton is enhanced by a collective excitation occurring on proton-ordered water wires, whose spectral intensity scales with the ordering length of water wire. The spectral intensity and excitonic interaction strength reaches its maximum in ice XI, where the long-range ordering length yields the most pronounced spectral signal. Our findings suggest that water wires, which widely exist in important physiological and biological systems and other phases of ice, can be directly probed by this approach. Published by the American Physical Society 2025

Tang, Fujie (ORCID:0000000197615359)↗
DOE OSTI2022

The understudied winter: Evidence of how precipitation differences affect stream metabolism in a headwater

Climate change is causing pronounced shifts during winter in the US, including shortening the snow season, reducing snowpack, and altering the timing and volume of snowmelt-related runoff. These changes in winter precipitation patterns affect in-stream freeze-thaw cycles, including ice and snow cover, and can trigger direct and indirect effects on in-stream physical, chemical, and biological processes in ~60% of river basins in the Northern Hemisphere. We used high-resolution, multi-parameter data collected in a headwater stream and its local environment (climate and soil) to determine interannual variability in physical, chemical, and biological signals in a montane stream during the winter of an El Niño and a La Niña year. We observed ~77% greater snow accumulation during the El Niño year, which caused the formation of an ice dam that shifted the system from a primarily lotic to a lentic environment. Water chemistry and stream metabolism parameters varied widely between years. They featured anoxic conditions lasting over a month, with no observable gross primary production (GPP) occurring under the ice and snow cover in the El Niño year. In contrast, dissolved oxygen and GPP remained relatively high during the winter months of the La Niña year. These redox and metabolic changes driven by changes in winter precipitation have significant implications for water chemistry and biological functioning beyond the winter. Our study suggests that as snow accumulation and hydrologic conditions shift during the winter due to climate change, hot-spots and hot-moments for biogeochemical processing may be reduced, with implications for the downstream movement of nutrients and transported materials.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Comparison of Shadowgraphy and X-Ray Phase Contrast Methods for Characterizing a DT Ice Layer in an Inertial Confinement Fusion Target

Shadowgraphy and x-ray phase contrast (XPC) imaging are two techniques that are used for characterizing the deuterium-tritium ice layer in Inertial Confinement Fusion targets. Each technique has limitations that affect how accurately they can characterize small crystalline defects and measure the ice-thickness non-uniformities that may be only a few micrometers in height. The concern is that shadowgraphy may be overly sensitive to the shape and depth of defects in the surface and insufficiently sensitive to the shape of longer-wavelength roughness, while x-ray phase contrast may be too insensitive to defects in the ice surface. Multiple ice layers with different thicknesses (40 to 63 μm), thickness uniformities (peak-valley variations that range from <2- to 12 μm), and crystal defects were analyzed using shadowgraphy and x-ray phase contrast techniques. The results from each method agree when the ice layer is uniformly thick and the crystal lacks defects. That agreement worsens as the number of defects in the surface of the ice layer increases and the roughness that is determined from a shadowgram image of the target’s limb is greater than can be justified by the number of defects that are seen in the target’s front and rear surfaces. The XPC technique is considerably less sensitive to surface defects, in part due to the poorer dynamic range and image resolution compared to shadowgraphy. Localized regions of the ice layer that are thicker or thinner than the average thickness of the layer are reported by shadowgraphy to be smaller in height and footprint (by up to 30%) than by XPC. As a result, the two techniques report different ice layer thicknesses that can vary by up to 10%. Shadowgraphy results from two caustics that trace different paths through the target, and, in theory, image the same ice-vapor surface (but reflect from either the vapor or ice side of the interface) did not consistently characterize the size or shape of ice features to be the same magnitude. Here, the XPC technique provides the best assessment of low-mode (l < 7) roughness in the ice layer. Shadowgraphy results using the strongest caustic is best for detecting the presence of grooves in the ice, though not for quantifying the size of them. If multiple grooves are present it is best to discard and reform the ice layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dissecting Anvil Cloud Response to Sea Surface Warming

Abstract We derive an anvil cloud diagnostic from the continuity equation of cloud ice and apply it to the output of convection‐permitting Energy Exascale Earth System Model (E3SM) simulations run in radiative‐convective equilibrium mode. This diagnostic shows that anvil cloud fraction can be reliably diagnosed as a product of cloud detrainment and lifetime. Detrainment is found to be approximated well by a product of clear sky convergence and cloud ice mixing ratio, while cloud lifetime is dominated by sedimentation. Taken together, this diagnostic expresses anvil cloud fraction as a function of five physically measurable quantities. Of these, clear‐sky convergence changes drive the anvil cloud reduction with warming while an increase in cloud ice mixing ratio buffers the decrease. Accordingly, this study provides a theoretical foundation upon which the Stability‐Iris hypothesis can be tested.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Imaging the magnetic nanowire cross section and magnetic ordering within a suspended 3D artificial spin-ice

Artificial spin-ice systems are patterned arrays of magnetic nanoislands arranged into frustrated geometries and provide insight into the physics of ordering and emergence. The majority of these systems have been realized in two-dimensions, mainly due to the ease of fabrication, but with recent developments in advanced nanolithography, three-dimensional artificial spin ice (ASI) structures have become possible, providing a new paradigm in their study. Such artificially engineered 3D systems provide new opportunities in realizing tunable ground states, new domain wall topologies, monopole propagation, and advanced device concepts, such as magnetic racetrack memory. Direct imaging of 3DASI structures with magnetic force microscopy has thus far been key to probing the physics of these systems but is limited in both the depth of measurement and resolution, ultimately restricting measurement to the uppermost layers of the system. In this work, a method is developed to fabricate 3DASI lattices over an aperture using two-photon lithography, thermal evaporation, and oxygen plasma exposure, allowing the probe of element-specific structural and magnetic information using soft x-ray microscopy with x-ray magnetic circular dichroism (XMCD) as magnetic contrast. The suspended polymer–permalloy lattices are found to be stable under repeated soft x-ray exposure. Analysis of the x-ray absorption signal allows the complex cross section of the magnetic nanowires to be reconstructed and demonstrates a crescent-shaped geometry. Measurement of the XMCD images after the application of an in-plane field suggests a decrease in magnetic moment on the lattice surface due to oxidation, while a measurable signal is retained on sub-lattices below the surface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Experimental realization of classical Z 2 spin liquids in a programmable quantum device

We build and probe a Z 2 spin liquid in a programmable quantum device, the D-Wave DW-2000Q. Specifically, we observe the classical eight-vertex and six-vertex (spin ice) states and transitions between them. To realize this state of matter, we design a Hamiltonian with combinatorial gauge symmetry using only pairwise-qubit interactions and a transverse field, i.e., interactions which are accessible in this quantum device. The combinatorial gauge symmetry remains exact along the full quantum annealing path, landing the system onto the classical eight-vertex model at the endpoint of the path. The output configurations from the device allow us to directly observe the loop structure of the classical model. Moreover, we deform the Hamiltonian so as to vary the weights of the eight vertices and show that we can selectively attain the classical six-vertex (ice) model, or drive the system into a ferromagnetic state. Additionally, we present studies of the classical phase diagram of the system as a function of the eight-vertex deformations and effective temperature, which we control by varying the relative strengths of the programmable couplings, and we show that the experimental results are consistent with theoretical analysis. Finally, we identify additional capabilities that, if added to these devices, would allow us to realize Z 2 quantum spin liquids on which to build topological qubits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Reconfigurable Pinwheel Artificial-Spin-Ice and Superconductor Hybrid Device

The ability to control the potential landscape in a medium of interacting particles could lead to intriguing collective behavior and innovative functionalities. In this work, we utilize spatially reconfigurable magnetic potentials of a pinwheel artificial-spin-ice (ASI) structure to tailor the motion of superconducting vortices. The reconstituted chain structures of the magnetic charges in the pinwheel ASI and the strong interaction between magnetic charges and superconducting vortices allow significant modification of the transport properties of the underlying superconducting thin film, resulting in a reprogrammable resistance state that enables a reversible and switchable vortex Hall effect. Our results highlight an effective and simple method of using ASI as an in situ reconfigurable nanoscale energy landscape to design reprogrammable superconducting electronics, which could also be applied to the in situ control of properties and functionalities in other magnetic particle systems, such as magnetic skyrmions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Diagnosing the sensitivity of grounding-line flux to changes in sub-ice-shelf melting

Abstract. Using a numerical ice flow model, we study changes in ice shelf buttressing and grounding-line flux due to localized ice thickness perturbations, a proxy for localized changes in sub-ice-shelf melting. From our experiments, applied to idealized (MISMIP+) and realistic (Larsen C) ice shelf domains, we identify a correlation between a locally derived buttressing number on the ice shelf, based on the first principal stress, and changes in the integrated grounding-line flux. The origin of this correlation, however, remains elusive from the perspective of a theoretical or physically based understanding. This and the fact that the correlation is generally much poorer when applied to realistic ice shelf domains motivate us to seek an alternative approach for predicting changes in grounding-line flux. We therefore propose an adjoint-based method for calculating the sensitivity of the integrated grounding-line flux to local changes in ice shelf geometry. We show that the adjoint-based sensitivity is identical to that deduced from pointwise, diagnostic model perturbation experiments. Based on its much wider applicability and the significant computational savings, we propose that the adjoint-based method is ideally suited for assessing grounding-line flux sensitivity to changes in sub-ice-shelf melting.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Mode structure measurements of ion cyclotron emission and sub-cyclotron modes on DIII-D

We report mode structure measurements of coherent ion cyclotron emission (ICE) and sub-cyclotron modes on DIII-D. Through a dedicated experiment, we aimed to characterize a variety of modes via the upgraded ICE diagnostic in both L- and H-mode plasmas. In the L-mode plasmas, autopower spectrum peaks at harmonics of the ion cyclotron frequency fci were observed and are localized in the core of the plasma. Sub-cyclotron modes (f ~ 0.5f ci ) were also observed in L-mode plasmas when the toroidal magnetic field strength was lowered from B T = 2.17 T. In H-mode plasmas, many ICE harmonics localized to the plasma edge were observed, with some exceeding the nominal ICE diagnostic bandwidth of f ϵ [0, 100] MHz. Polarization estimates made using an orthogonal pair of vertical and horizontal magnetic pickup loops on the outer wall of the machine indicate that ICE and sub-cyclotron modes have compressional polarization at the plasma edge, the latter being consistent with simulation efforts for comparable DIII-D plasmas. For all modes, the same harmonics are observed on both the centerpost and outer wall loops, indicating that ICE is poloidally extended. Finally, toroidal mode numbers were calculated using three outer wall loops for both core ICE and sub-cyclotron modes in L-mode plasmas. Furthermore, the sub-cyclotron case served as a benchmark for our calculation method, with measured numbers of n ϵ [–24, –18] roughly agreeing with heuristic estimates of n ϵ [–20, –13]. Core ICE mode numbers were measured to be n ϵ [–10, 5], which is comparable to measurements made on other machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗