SEARCH · Engineering Papers

Results for “ice”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

A LaGrangian Perspective of the Surface Energy Budget Contributions to Sea Ice Evolution in the ‘New Arctic’

Recent Arctic warming outpaces global warming and is characterized by dramatic changes in sea ice (record sea ice loss, thinner and more seasonal ice cover), warming SSTs, and a lengthening melt season. These changes impact the surface energy budget and the sensitivity of sea ice to radiative forcing, namely an increased sensitivity of thinner ice to local radiative forcings and a weaker buffering of surface-atmosphere exchanges of heat and moisture. However, these processes are a large source of uncertainty in global climate models with prior studies finding that poor representations of downwelling longwave radiation and sea ice albedo lead to model spread in sea ice conditions. A limitation of prior work is examining sea ice-atmosphere interactions through an Eulerian framework. This study utilizes a sea ice parcel dataset in a LaGrangian framework, allowing for tracking of surface energy budget-sea ice interactions over the lifespan of individual sea ice parcels; this framework also allows for the quantification of the sea ice response to episodic weather events along with seasonal drivers, including the drivers of rapid growth and melt events. First, lag-lead relationships between relevant surface energy budget terms and sea ice thickness changes are explored to understand the atmospheric conditions preceding rapid growth and melt. The impact of episodic events- moisture intrusion, warm air outbreaks, cyclones- on sea ice growth are also investigated by stratifying first-year parcels by sea ice thickness at the end of the growth season and analyzing the effect of atmospheric variability (number and intensity of episodic events) on parcel thickness change. Lastly, parcel growth and melt and the associated radiative anomalies are compared to snowfall events/snow thickness to determine how snowfall modifies the interaction between the local surface energy budget and sea ice mass balance.

Patrick Taylor↗

DOE OSTI2022

Simulating ice-shelf extent using damage mechanics

Abstract Inaccurate representations of iceberg calving from ice shelves are a large source of uncertainty in mass-loss projections from the Antarctic ice sheet. Here, we address this limitation by implementing and testing a continuum damage-mechanics model in a continental scale ice-sheet model. The damage-mechanics formulation, based on a linear stability analysis and subsequent long-wavelength approximation of crevasses that evolve in a viscous medium, links damage evolution to climate forcing and the large-scale stresses within an ice shelf. We incorporate this model into the BISICLES ice-sheet model and test it by applying it to idealized (1) ice tongues, for which we present analytical solutions and (2) buttressed ice-shelf geometries. Our simulations show that the model reproduces the large disparity in lengths of ice shelves with geometries and melt rates broadly similar to those of four Antarctic ice shelves: Erebus Glacier Tongue (length ~ 13 km), the unembayed portion of Drygalski Ice Tongue (~ 65 km), the Amery Ice Shelf (~ 350 km) and the Ross Ice Shelf (~ 500 km). These results demonstrate that our simple continuum model holds promise for constraining realistic ice-shelf extents in large-scale ice-sheet models in a computationally tractable manner.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Glacier ice archives nearly 15,000-year-old microbes and phages

Background Glacier ice archives information, including microbiology, that helps reveal paleoclimate histories and predict future climate change. Though glacier-ice microbes are studied using culture or amplicon approaches, more challenging metagenomic approaches, which provide access to functional, genome-resolved information and viruses, are under-utilized, partly due to low biomass and potential contamination. Results We expand existing clean sampling procedures using controlled artificial ice-core experiments and adapted previously established low-biomass metagenomic approaches to study glacier-ice viruses. Controlled sampling experiments drastically reduced mock contaminants including bacteria, viruses, and free DNA to background levels. Amplicon sequencing from eight depths of two Tibetan Plateau ice cores revealed common glacier-ice lineages including Janthinobacterium, Polaromonas, Herminiimonas, Flavobacterium, Sphingomonas, and Methylobacterium as the dominant genera, while microbial communities were significantly different between two ice cores, associating with different climate conditions during deposition. Separately, ~355- and ~14,400-year-old ice were subject to viral enrichment and low-input quantitative sequencing, yielding genomic sequences for 33 vOTUs. These were virtually all unique to this study, representing 28 novel genera and not a single species shared with 225 environmentally diverse viromes. Further, 42.4% of the vOTUs were identifiable temperate, which is significantly higher than that in gut, soil, and marine viromes, and indicates that temperate phages are possibly favored in glacier-ice environments before being frozen. In silico host predictions linked 18 vOTUs to co-occurring abundant bacteria (Methylobacterium, Sphingomonas, and Janthinobacterium), indicating that these phages infected ice-abundant bacterial groups before being archived. Functional genome annotation revealed four virus-encoded auxiliary metabolic genes, particularly two motility genes suggest viruses potentially facilitate nutrient acquisition for their hosts. Finally, given their possible importance to methane cycling in ice, we focused on Methylobacterium viruses by contextualizing our ice-observed viruses against 123 viromes and prophages extracted from 131 Methylobacterium genomes, revealing that the archived viruses might originate from soil or plants. Conclusions Together, these efforts further microbial and viral sampling procedures for glacier ice and provide a first window into viral communities and functions in ancient glacier environments. Such methods and datasets can potentially enable researchers to contextualize new discoveries and begin to incorporate glacier-ice microbes and their viruses relative to past and present climate change in geographically diverse regions globally.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2014

Investigating ice nucleation in cirrus clouds with an aerosol-enabled Multiscale Modeling Framework

In this study, an aerosol-dependent ice nucleation scheme [Liu and Penner, 2005] has been implemented in an aerosol-enabled multi-scale modeling framework (PNNL MMF) to study ice formation in upper troposphere cirrus clouds through both homogeneous and heterogeneous nucleation. The MMF model represents cloud scale processes by embedding a cloud-resolving model (CRM) within each vertical column of a GCM grid. By explicitly linking ice nucleation to aerosol number concentration, CRM-scale temperature, relative humidity and vertical velocity, the new MMF model simulates the persistent high ice supersaturation and low ice number concentration (10 to 100/L) at cirrus temperatures. The low ice number is attributed to the dominance of heterogeneous nucleation in ice formation. The new model simulates the observed shift of the ice supersaturation PDF towards higher values at low temperatures following homogeneous nucleation threshold. The MMF models predict a higher frequency of midlatitude supersaturation in the Southern hemisphere and winter hemisphere, which is consistent with previous satellite and in-situ observations. It is shown that compared to a conventional GCM, the MMF is a more powerful model to emulate parameters that evolve over short time scales such as supersaturation. Sensitivity tests suggest that the simulated global distribution of ice clouds is sensitive to the ice nucleation schemes and the distribution of sulfate and dust aerosols. Simulations are also performed to test empirical parameters related to auto-conversion of ice crystals to snow. Results show that with a value of 250 μm for the critical diameter, Dcs, that distinguishes ice crystals from snow, the model can produce good agreement to the satellite retrieved products in terms of cloud ice water path and ice water content, while the total ice water is not sensitive to the specification of Dcs value.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling

The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer ice-philicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface–ice contact at the expense of surface–water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Furthermore, our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Sea ice heat and mass balance measurements from four autonomous buoys during the MOSAiC drift campaign

As part of the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC), four autonomous seasonal ice mass balance buoys were deployed in first- and second-year ice. These buoys measured position, barometric pressure, snow depth, ice thickness, ice growth, surface melt, bottom melt, and vertical profiles of temperature from the air, through the snow and ice, and into the upper ocean. Observed air temperatures were similar at all four sites; however, snow–ice interface temperatures varied by as much as 10°C, primarily due to differences in snow depth. Observed winter ice growth rates (November to May) were <1 cm day -1 , with summer melt rates (June to July) as large as 5 cm day -1 . Air temperatures changed as much as 2°C hour -1 but were dampened to <0.3°C hour -1 at the snow–ice interface. Initial October ice thicknesses ranged from 0.3 m in first-year ice to 1.2 m in second-year ice. By February, this range was only 1.20–1.46 m, due in part to differences in the onset of basal freezing. In second-year ice, this delay was due to large brine-filled voids in the ice; propagating the cold front through this ice required freezing the brine. Mass balance results were similar to those measured by autonomous buoys deployed at the North Pole from 2000 to 2013. Winter average estimates of the ocean heat flux ranged from 0 to 3 W m -2 , with a large increase in June 2020 as the floe moved into warmer water. Estimates of average snow thermal conductivity measured at two buoys during periods of linear temperature profiles were 0.41 and 0.42 W m -1 °C -1 , higher than previously published estimates. Results from these ice mass balance buoys can contribute to efforts to close the MOSAiC heat budget.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2012

Determination of Ice Cloud Models Using MODIS and MISR Data

Representation of ice clouds in radiative transfer simulations is subject to uncertainties associated with the shapes and sizes of ice crystals within cirrus clouds. In this study, we examined several ice cloud models consisting of smooth, roughened, homogeneous and inhomogeneous hexagonal ice crystals with various aspect ratios. The sensitivity of the bulk scattering properties and solar reflectances of cirrus clouds to specific ice cloud models is investigated using the improved geometric optics method (IGOM) and the discrete ordinates radiative transfer (DISORT) model. The ice crystal habit fractions in the ice cloud model may significantly affect the simulations of cloud reflectances. A new algorithm was developed to help determine an appropriate ice cloud model for application to the satellite-based retrieval of ice cloud properties. The ice cloud particle size retrieved from Moderate Resolution Imaging Spectroradiometer (MODIS) data, collocated with Multi-angle Imaging Spectroradiometer (MISR) observations, is used to infer the optical thicknesses of ice clouds for nine MISR viewing angles. The relative differences between view-dependent cloud optical thickness and the averaged value over the nine MISR viewing angles can vary from −0.5 to 0.5 and are used to evaluate the ice cloud models. In the case for 2 July 2009, the ice cloud model with mixed ice crystal habits is the best fit to the observations (the root mean square (RMS) error of cloud optical thickness reaches 0.365). This ice cloud model also produces consistent cloud property retrievals for the nine MISR viewing configurations within the measurement uncertainties.

MODIS↗

NASA NTRSDate not supplied

Effect of Geometric Fidelity on the Aerodynamics of a Swept Wing with Glaze Ice Accretion

Aerodynamic assessment of icing effects on swept wings is an important component of a larger effort to improve three-dimensional icing simulation capabilities. An understanding of ice-shape geometric fidelity on iced-wing aerodynamics and the associated flowfield features are needed to guide the development and validation of ice-accretion simulation tools. To this end, wind-tunnel testing was carried out for 8.9% and 13.3% scale semispan wing models based upon the Common Research Model airplane configuration. Various levels of geometric fidelity of an artificial ice shape representing a glaze-ice accretion on a swept wing were investigated. The highest fidelity artificial ice shape reproduced all of the three-dimensional features associated with the glaze ice accretion. The lowest fidelity artificial ice shapes were simple, spanwise-varying horn ice geometries intended to represent the maximum ice thickness on the wing upper surface. The results presented in this paper show that the addition of grit roughness to some lower-fidelity artificial ice shapes resulted in favorable lift and pitching moment comparisons to the wing with the highest fidelity artificial ice shape. In the range of 4.3 to 7.4 deg. angle of attack, surface oil flow visualization and pressure data show that the wing with the two lower fidelity simulations clearly demonstrated a leading edge vortex dominated flowfield, referred to as type I. For the wing with the high fidelity ice shape, the flowfield at lower angles of attack was characterized by streamwise-running, counter-rotating vortical flow referred to as type II. Between 6.4 and 7.4 deg. angle of attack, the effect of the type II flow structures was significantly altered and gave way to the type I leading edge vortex. This means that for angles of attack 7.4 deg. and higher, the wing with all three configurations exhibited the same type of flowfield. This helps to explain why there is reasonably good agreement in the lift and pitching moment coefficients among these configurations.

Aircraft icing↗

NASA NTRS1947

Laboratory Investigation of Ice Formation and Elimination in the Induction System of a Large Twin-engine Cargo Aircraft

The icing characteristics, the de-icing rate with hot air, and the effect of impact ice on fuel metering and mixture distribution have been determined in a laboratory investigation of that part of the engine induction system consisting of a three-barrel injection-type carburetor and a supercharger housing with spinner-type fuel injection from an 18-cylinder radial engine used on a large twin-engine cargo airplane. The induction system remained ice-free at carburetor-air temperatures above 36 F regardless of the moisture content of the air. Between carburetor-air temperatures of 32 F and 36 F with humidity ratios in excess of saturation, serious throttling ice formed in the carburetor because of expansion cooling of the air; at carburetor-air temperatures below 32 F with humidity ratios in excess of saturation, serious impact-ice formations occurred, Spinner-type fuel injection at the entrance to the supercharger and heating of the supercharger-inlet elbow and the guide vanes by the warn oil in the rear engine housing are design features that proved effective in eliminating fuel-evaporation icing and minimized the formation of throttling ice below the carburetor. Air-flow recovery time with fixed throttle was rapidly reduced as the inlet -air wet -bulb temperature was increased to 55 F; further temperature increase produced negligible improvement in recovery time. Larger ice formations and lower icing temperatures increased the time required to restore proper air flow at a given wet-bulb temperature. Impact-ice formations on the entrance screen and the top of the carburetor reduced the over-all fuel-air ratio and increased the spread between the over-all ratio and the fuel-air ratio of the individual cylinders. The normal spread of fuel-air ratio was increased from 0.020 to 0.028 when the left quarter of the entrance screen was blocked in a manner simulating the blocking resulting from ice formations released from upstream duct walls during hot-air de-icing.

ICE PREVENTION - AIR HEATING↗

NASA NTRS2021

Sea Ice Extents Continue to Set New Records: Arctic, Antarctic, and Global Results

The multi-channel satellite passive-microwave record of Earth’s sea ice coverage, extending back to the late 1970s, has long revealed declining sea ice coverage in the Arctic but through 2015 revealed an overall increase rather than decrease in Antarctic sea ice coverage. Following major decreases in Antarctic sea ice since 2015, the 42-year 1979-2020 satellite dataset now shows losses in sea ice coverage in both hemispheres, and this is convincingly demonstrated by the enumeration of monthly and yearly record high and record low sea ice extents experienced over the course of the 42 years. In fact, one of the most convincing statistics on the declining Arctic sea ice cover is the fact that since 1986 the Arctic has not experienced a single monthly record high sea ice extent in any month but has experienced 93 monthly record lows. In contrast, all 12 calendar months have their 42-year Antarctic monthly record high sea ice extents in the period 2007-2015, while 8 of the 12 calendar months have had Antarctic record lows since 2015. Globally, every calendar month has registered a new monthly record low within the past 5 years. These results are complemented (and somewhat tempered) by quantification of the range of monthly and yearly sea ice extent values over the 42 years. For instance, although the Arctic’s lowest September monthly average sea ice extent (in 2012) is 53% lower than its highest September monthly average sea ice extent (in 1980), the other months have far smaller percent differences between their lowest and highest Arctic values. For yearly average sea ice extents, the Arctic’s lowest value (in 2020) is 18% lower than its highest value (in 1982), the Antarctic’s lowest value (in 2017) is 16% lower than its highest value (in 2014), and the global lowest value (in 2019) is only 12% lower than its highest value (in 1982).

sea ice↗

DOE OSTI2019

Deglaciation of the Greenland and Laurentide ice sheets interrupted by glacier advance during abrupt coolings

Understanding patterns of ice-sheet deglaciation is key for predicting the rate of future ice-sheet melt, yet the processes underlying deglaciation remain elusive. The early Holocene (11.7 ka to 8.2 ka; Greenlandian Stage) represents the most recent period when the Laurentide and Greenland ice sheets underwent large-scale recession. Moreover, this ice-sheet recession occurred under the backdrop of regional temperatures that were similar to or warmer than today, and comparable to those projected for the upcoming centuries. Reconstructing Laurentide and Greenland ice sheet behavior during the early Holocene, and elucidating the mechanisms dictating this behavior may serve as a partial analog for future Greenland ice-sheet change in a warming world. Here, we present 123 new 10 Be surface exposure ages from two sites on Baffin Island and southwestern Greenland that constrain the behavior of the Laurentide and Greenland ice sheets, and an independent alpine glacier during the early Holocene. On Baffin Island, we focus on a unique area where moraines deposited by the Laurentide Ice Sheet rest directly adjacent to moraines deposited by an independent alpine glacier. Sixty-one 10 Be ages reveal that advances and/or stillstands of the Laurentide Ice Sheet and an alpine glacier occurred in unison around 11.8 ka, 10.3 ka, and 9.2 ka. Sixty-two 10 Be ages from southwestern Greenland indicate that the Greenland Ice Sheet margin experienced re-advances or stillstands around 11.6 ka, 10.4 ka, 9.1 ka, 8.1 ka, and 7.3 ka. Our results reveal that ice sheets respond to climate perturbations on the same centennial timescale as small alpine glaciers. Further, we hypothesize that during the warming climate of the early Holocene, freshening of the North Atlantic Ocean induced by ice-sheet melt resulted in regional cooling and brief periods of ice-sheet stabilization superimposed on net glacier recession. These observations point to a negative feedback mechanism inherent to melting ice sheets in the Baffin Bay region that slows ice-sheet recession during intervals of otherwise rapid deglaciation.

58 GEOSCIENCES↗

DOE OSTI2021

Chirping ion cyclotron emission (ICE) on NSTX-U

In this work, we report on the discovery on NSTX-U of a qualitatively new type of non-thermal emission from plasmas in the ion-cyclotron frequency range (ICE). This new type of ICE is weakly damped, γ damp /ω ICE ≈ 0.14%, low order eigenmodes, excited by a fast ion population radially located near an internal transport barrier. The ICE appears as repetitive bursts with durations of ≈1 ms, and the bursts exhibit frequency chirps upwards and/or downwards with maximum δf/f ≈ 1%. We refer to this chirping ICE as ch-ICE. The bursts are longer than the typical ICE bursts of ≈100 μs reported previously for NSTX and NSTX-U, for simplicity, NSTX(-U). The chirping resembles that seen in modeling of weakly unstable fast particle driven instabilities (cf, Berk et al 1997 Phys. Lett. A 234 213). These data also provide some preliminary indications of non-linear coupling between the ch-ICE harmonics. Despite extensive experimental and theoretical studies of ICE, theoretical understanding of ICE remains qualitative at best. The detailed experimental observations of ICE from NSTX(-U) find a wide variation in the ICE characteristics. The new observations of ch-ICE presented here, combined with the characteristics of ICE on NSTX(-U) previously reported (Fredrickson et al 2019 Phys. Plasmas 26 032111), could guide the development of more complete theoretical models of ICE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

A Comprehensive Parameterization of Heterogeneous Ice Nucleation of Dust Surrogate: Laboratory Study with Hematite Particles and Its Application to Atmospheric Models

A new heterogeneous ice nucleation parameterization that covers a wide temperature range (-36 °C to -78 °C) is presented. Developing and testing such an ice nucleation parameterization, which is constrained through identical experimental conditions, is critical in order to accurately simulate the ice nucleation processes in cirrus clouds. The surface-scaled ice nucleation efficiencies of hematite particles, inferred by n s , were derived from AIDA (Aerosol Interaction and Dynamics in the Atmosphere) cloud chamber measurements under water subsaturated conditions that were realized by continuously changing temperature (T) and relative humidity with respect to ice (RH ice ) in the chamber. Our measurements showed several different pathways to nucleate ice depending on T and RH ice conditions. For instance, almost independent freezing was observed at -60 °C < T < -50 °C, where RH ice explicitly controlled ice nucleation efficiency, while both T and RH ice played roles in other two T regimes: -78 °C < T < -60 °C and -50 °C < T < -36 °C. More specifically, observations at T colder than -60 °C revealed that higher RHice was necessary to maintain constant n s , whereas T may have played a significant role in ice nucleation at T warmer than -50 °C. We implemented new n s parameterizations into two cloud models to investigate its sensitivity and compare with the existing ice nucleation schemes towards simulating cirrus cloud properties. Our results show that the new AIDA-based parameterizations lead to an order of magnitude higher ice crystal concentrations and inhibition of homogeneous nucleation in colder temperature regions. Our cloud simulation results suggest that atmospheric dust particles that form ice nuclei at lower temperatures, below -36 °C, can potentially have stronger influence on cloud properties such as cloud longevity and initiation when compared to previous parameterizations.

Hiranuma, Naruki↗

DOE OSTI2019

Laboratory study of the heterogeneous ice nucleation on black-carbon-containing aerosol

Abstract. Soot and black carbon (BC) particles are generated in the incomplete combustion of fossil fuels, biomass, and biofuels. These airborne particles affect air quality, human health, aerosol–cloud interactions, precipitation formation, and climate. At present, the climate effects of BC particles are not well understood. Their role in cloud formation is obscured by their chemical and physical variability and by the internal mixing states of these particles with other compounds. Ice nucleation in field studies is often difficult to interpret. Nonetheless, most field studies seem to suggest that BC particles are not efficient ice-nucleating particles (INPs). On the other hand, laboratory measurements show that in some cases, BC particles can be highly active INPs under certain conditions. By working with well-characterized BC particles, our aim is to systematically establish the factors that govern the ice nucleation activity of BC. The current study focuses on laboratory measurements of the effectiveness of BC-containing aerosol in the formation of ice crystals in temperature and ice supersaturation conditions relevant to cirrus clouds. We examine ice nucleation on BC particles under water-subsaturated cirrus cloud conditions, commonly understood as deposition-mode ice nucleation. We study a series of well-characterized commercial carbon black particles with varying morphologies and surface chemistries as well as ethylene flame-generated combustion soot. The carbon black particles used in this study are proxies for atmospherically relevant BC aerosols. These samples were characterized by electron microscopy, mass spectrometry, and optical scattering measurements. Ice nucleation activity was systematically examined in temperature and saturation conditions in the ranges of 217≤T≤235 K and 1.0≤Sice≤1.5 and 0.59≤Swater≤0.98, respectively, using a SPectrometer for Ice Nuclei (SPIN) instrument, which is a continuous-flow diffusion chamber coupled with instrumentation to measure light scattering and polarization. To study the effect of coatings on INPs, the BC-containing particles were coated with organic acids found in the atmosphere, namely stearic acid, cis-pinonic acid, and oxalic acid. The results show significant variations in ice nucleation activity as a function of size, morphology, and surface chemistry of the BC particles. The measured ice nucleation activity dependencies on temperature, supersaturation conditions, and the physicochemical properties of the BC particles are consistent with an ice nucleation mechanism of pore condensation followed by freezing. Coatings and surface oxidation modify the initial formation efficiency of pristine ice crystals on BC-containing aerosol. Depending on the BC material and the coating, both inhibition and enhancement in INP activity were observed. Our measurements at low temperatures complement published data and highlight the capability of some BC particles to nucleate ice under low ice supersaturation conditions. These results are expected to help refine theories relating to soot INP activation in the atmosphere.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Primary and secondary ice production: interactions and their relative importance

A discrepancy of up to 5 orders of magnitude between ice crystal and ice nucleating particle (INP) number concentrations was found in the measurements, indicating the potentially important role of secondary ice production (SIP) in the clouds. However, the interactions between primary and SIP processes and their relative importance remain unexplored. In this study, we implemented five different ice nucleation schemes as well as physical representations of SIP processes (i.e., droplet shattering during rain freezing, ice-ice collisional break-up, and rime splintering) in the Community Earth System Model version 2 (CESM2). We ran CESM2 in the single column mode for model comparisons with the DOE Atmospheric Radiation Measurement (ARM) Mixed-Phase Arctic Cloud Experiment (M-PACE) observations. We found that the model experiments with aerosol-aware ice nucleation schemes and SIP processes yield the best simulation results for the M-PACE single-layer mixed-phase clouds. We further investigated the relative importance of ice nucleation and SIP to ice number and cloud phase as well as interactions between ice nucleation and SIP in the M-PACE single-layer mixed-phase clouds. Our results show that SIP contributes 80 % to the total ice formation and transforms ~30 % of pure liquid-phase clouds simulated in the model experiments without considering SIP into mixed-phase clouds. The SIP is not only a result of ice crystals produced from ice nucleation, but also competes with the ice nucleation by reducing the number concentrations of cloud droplets and cloud-borne dust INPs. Conversely, strong ice nucleation also suppresses SIP by glaciating mixed-phase clouds and thereby reducing the amount of precipitation particles (rain and graupel).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD

We present a new capability of the ice sheet model SICOPOLIS that enables flexible adjoint code generation via source transformation using the open-source algorithmic differentiation (AD) tool OpenAD.The adjoint code enables efficient calculation of the sensitivities of a scalar-valued objective function or quantity of interest (QoI) to a range of important, often spatially varying and uncertain model input variables, including initial and boundary conditions, as well as model parameters.Compared to earlier work on the adjoint code generation of SICOPOLIS, our work makes several important advances:(i) it is embedded within the up-to-date trunk of the SICOPOLIS repository – accounting for 1.5 decades of code development and improvements – and is readily available to the wider community;(ii) the AD tool used, OpenAD, is an open-source tool;(iii) the adjoint code developed is applicable to both Greenland and Antarctica, including grounded ice as well as floating ice shelves, with an extended choice of thermodynamical representations.A number of code refactorization steps were required. They are discussed in detail in an Appendix as they hold lessons for the application of AD to legacy codes at large.As an example application, we examine the sensitivity of the total Antarctic Ice Sheet volume to changes in initial ice thickness, austral summer precipitation, and basal and surface temperatures across the ice sheet.Simulations of Antarctica with floating ice shelves show that over 100 years of simulation the sensitivity of total ice sheet volume to the initial ice thickness and precipitation is almost uniformly positive, while the sensitivities to surface and basal temperature are almost uniformly negative. Sensitivity to austral summer precipitation is largest on floating ice shelves from Queen Maud to Queen Mary Land. The largest sensitivity to initial ice thickness is at outlet glaciers around Antarctica. Comparison between total ice sheet volume sensitivities to surface and basal temperature shows that surface temperature sensitivities are higher broadly across the floating ice shelves, while basal temperature sensitivities are highest at the grounding lines of floating ice shelves and outlet glaciers. A uniformly perturbed region of East Antarctica reveals that, among the four control variables tested here, total ice sheet volume is the most sensitive to variations in austral summer precipitation as formulated in SICOPOLIS.Comparison between adjoint- and finite-difference-derived sensitivities shows good agreement, lending confidence that the AD tool is producing correct adjoint code.The new modeling infrastructure is freely available at http://www.sicopolis.net (last access: 2 April 2020) under the development trunk.

58 GEOSCIENCES↗

DOE OSTI2020

ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century

Abstract. Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between −7.8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between −6.1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

SNICAR-ADv4: a physically based radiative transfer model to represent the spectral albedo of glacier ice

Abstract. Accurate modeling of cryospheric surface albedo is essential for our understanding of climate change as snow and ice surfaces regulate the global radiative budget and sea-level through their albedo and mass balance. Although significant progress has been made using physical principles to represent the dynamic albedo of snow, models of glacier ice albedo tend to be heavily parameterized and not explicitly connected with physical properties that govern albedo, such as the number and size of air bubbles, specific surface area (SSA), presence of abiotic and biotic light absorbing constituents (LACs), and characteristics of any overlying snow. Here, we introduce SNICAR-ADv4, an extension of the multi-layer two-stream delta-Eddington radiative transfer model with the adding–doubling solver that has been previously applied to represent snow and sea-ice spectral albedo. SNICAR-ADv4 treats spectrally resolved Fresnel reflectance and transmittance between overlying snow and higher-density glacier ice, scattering by air bubbles of varying sizes, and numerous types of LACs. SNICAR-ADv4 simulates a wide range of clean snow and ice broadband albedo (BBA), ranging from 0.88 for (30 µm) fine-grain snow to 0.03 for bare and bubble-free ice under direct light. Our results indicate that representing ice with a density of 650 kg m−3 as snow with no refractive Fresnel layer, as done previously, generally overestimates the BBA by an average of 0.058. However, because most naturally occurring ice surfaces are roughened “white ice”, we recommend modeling a thin snow layer over bare ice simulations. We find optimal agreement with measurements by representing cryospheric media with densities less than 650 kg m−3 as snow and larger-density media as bubbly ice with a Fresnel layer. SNICAR-ADv4 also simulates the non-linear albedo impacts from LACs with changing ice SSA, with peak impact per unit mass of LACs near SSAs of 0.1–0.01 m2 kg−1. For bare, bubble-free ice, LACs actually increase the albedo. SNICAR-ADv4 represents smooth transitions between snow, firn, and ice surfaces and accurately reproduces measured spectral albedos of a variety of glacier surfaces. This work paves the way for adapting SNICAR-ADv4 to be used in land ice model components of Earth system models.

54 ENVIRONMENTAL SCIENCES↗