Engineering Papers⌕ Search

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.

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

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Characterization of Ice Roughness Variations in Scaled Glaze Icing Conditions

Because of the significant influence of surface tension in governing the stability and breakdown of the liquid film in flooded stagnation regions of airfoils exposed to glaze icing conditions, the Weber number is expected to be a significant parameter governing the formation and evolution of ice roughness. To investigate the influence of the Weber number on roughness formation, 53.3-cm (21-in.) and 182.9-cm (72-in.) NACA 0012 airfoils were exposed to flow conditions with essentially the same Weber number and varying stagnation collection efficiency to illuminate similarities of the ice roughness created on the different airfoils. The airfoils were exposed to icing conditions in the Icing Research Tunnel (IRT) at the NASA Glenn Research Center. Following exposure to the icing event, the airfoils were then scanned using a ROMER Absolute Arm scanning system. The resulting point clouds were then analyzed using the self-organizing map approach of McClain and Kreeger (2013) to determine the spatial roughness variations along the surfaces of the iced airfoils. The roughness characteristics on each airfoil were then compared using the relative geometries of the airfoil. The results indicate that features of the ice shape and roughness such as glaze-ice plateau limits and maximum airfoil roughness were captured well by Weber number and collection efficiency scaling of glaze icing conditions. However, secondary ice roughness features relating the instability and waviness of the liquid film on the glaze-ice plateau surface are scaled based on physics that were not captured by the local collection efficiency variations.

Roughness Formation↗
NASA NTRS2017

Ice Crystal Icing Research at NASA

Ice crystals found at high altitude near convective clouds are known to cause jet engine power-loss events. These events occur due to ice crystals entering a propulsion system's core flowpath and accreting ice resulting in events such as uncommanded loss of thrust (rollback), engine stall, surge, and damage due to ice shedding. As part of a community with a growing need to understand the underlying physics of ice crystal icing, NASA has been performing experimental efforts aimed at providing datasets that can be used to generate models to predict the ice accretion inside current and future engine designs. Fundamental icing physics studies on particle impacts, accretion on a single airfoil, and ice accretions observed during a rollback event inside a full-scale engine in the Propulsion Systems Laboratory are summarized. Low fidelity code development using the results from the engine tests which identify key parameters for ice accretion risk and the development of high fidelity codes are described. These activities have been conducted internal to NASA and through collaboration efforts with industry, academia, and other government agencies. The details of the research activities and progress made to date in addressing ice crystal icing research challenges are discussed.

crysta↗
DOE OSTI2015

Impacts of a mushy‐layer thermodynamic approach in global sea‐ice simulations using the CICE sea‐ice model

Abstract We perform global simulations of the Los Alamos sea‐ice model, CICE, with a new thermodynamics component that has a fully prognostic, variable bulk salinity vertical profile based on mushy layer physics. The processes of gravity drainage, melt‐water flushing and snow‐ice formation are parameterized to allow the bulk salinity to evolve with time. We analyze the seasonal and spatial variation of sea‐ice bulk salinity, area, volume and thickness and compare these quantities to simulations using the previous thermodynamic component. Adjusting one of the gravity drainage parameters, we find good agreement between simulation results and fieldwork ice‐core observations of sea‐ice bulk salinity. As expected, bulk salinity is highest during periods of ice growth and lowest after periods of ice melt. In the northern hemisphere the new thermodynamics component produces thicker ice than the previous thermodynamics component. Of the nine major differences between the two models, differences in how salinities are calculated and how melt‐pond flushing is parameterized are the principal causes of this thickness difference. Thickness differences are smaller in the southern hemisphere than in the northern hemisphere since a greater fraction of ice melts, and differences cannot accumulate year‐on‐year. Model differences in how ice thickness changes and snow‐ice formation are calculated are the most important causes of the different thickness between the two thermodynamic components in the southern hemisphere. The melt‐pond area and volume are found to be highly sensitive to a parameter choice controlling drainage through macroscopic holes in the ice, in both hemispheres.

Turner, Adrian K.↗
DOE OSTI2014

Micro-Spectroscopic Imaging and Characterization of Individually Identified Ice Nucleating Particles from a Case Field Study

The effect of anthropogenic and biogenic organic particles on atmospheric glaciation processes is poorly understood. We use an optical microscopy (OM) setup to identify the location of ice nuclei (IN) active in immersion freezing and deposition ice nucleation for temperatures of 200-273 K within a large population of particles sampled from an ambient environment. Applying multi-modal micro-spectroscopy methods we characterize the physicochemical properties of individual IN in particle populations collected in central California. Chemical composition and mixing state analysis of particle populations are performed to identify characteristic particle-type classes. All particle-types contained organic material. Particles in these samples take up water at subsaturated conditions, induce immersion freezing at subsaturated and saturated conditions above 226 K, and act as deposition IN below 226 K. The identified IN belong to the most common particle-type classes observed in the field samples: organic coated sea salt, Na-rich, and secondary and refractory carbonaceous particles. Based on these observations, we suggest that the IN are not always particles with unique chemical composition and exceptional ice nucleation propensity; rather, they are common particles in the ambient particle population. Thus, particle composition and morphology alone are insufficient to assess their potential to act as IN. The results suggest that particle-type abundance is also a crucial factor in determining the ice nucleation efficiency of specific IN types. These findings emphasize that ubiquitous organic particles can induce ice nucleation under atmospherically relevant conditions and that they may play an important role in atmospheric glaciation processes.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2016

Ice nucleation activity of diesel soot particles at Cirrus relevant conditions: Effects of hydration, secondary organics coating, hydration, soot morphology, and coagulation

The role of atmospheric relevant soot particles that are processed in the atmosphere toward ice nucleation at cirrus cloud condition is poorly understood. In this study, the ice nucleating properties of diesel soot particles subjected to various physical and chemical aging treatments were investigated at temperatures ranging from -40 to -50 °C. We show that bare soot particles nucleate ice in deposition mode, but coating with secondary organics suppresses the heterogeneous ice nucleation potential of soot particles requiring homogeneous freezing threshold conditions. However, the ice nucleation efficiency of soot particles coated with an aqueous organic layer was similar to bare soot particles. Hydration of bare soot particles slightly enhanced the ice nucleation efficiency, and the IN abilities of compact soot particles (roundness = ~ 0.6) were similar to bare lacey soot particles (roundness = ~ 0.4). These results indicate that ice nucleation properties are sensitive to the various aging treatments.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2025

On the Role of Airborne Ice Nucleating Particles in Primary and Secondary Ice Formation Processes in Convective Midlatitude Clouds

Abstract Formation of ice in mixed-phase clouds may proceed via primary nucleation and secondary ice production (SIP). Primary nucleation involves the presence of ice nucleating particles (INPs), while SIP follows primary nucleation, resulting in ice crystal number concentrations N ice exceeding the number of INPs. The ice formation pathways in young congestus clouds are not well understood nor are the environmental conditions favorable for the onset of SIP. Coincident airborne measurements of INPs and N ice collected during the Secondary Production of Ice in Cumulus Experiment (SPICULE) campaign over the U.S. central Great Plains are reported for young congestus clouds. Number concentrations of INPs and N ice at temperatures between −10° and −20°C were used to categorize whether SIP was active in each analyzed cloud pass. Observational data suggested that fragmentation of freezing droplets (FFDs) was occurring and the most likely SIP mechanism for ice enhancement during the early stages of the sampled clouds, where higher cloud-base temperatures and stronger updrafts produced conditions favorable for the onset of FFD. Simple model simulations of a congestus cloud were used to compare predictions of N ice from newly implemented SIP mechanisms (including FFD) alongside those from the existing Hallett–Mossop (HM) mechanism. Although predicted N ice , when all SIP mechanisms were active, was in good agreement with observations, the dominant SIP mechanism in the model was predicted to be HM. The rate of heterogeneous freezing in raindrops was likely inhibiting realistic FFD rates in the simulated cloud. These discrepancies between the observations and simulations underscore the need for extended laboratory studies to formulate improved model representations of SIP in convective clouds.

Patnaude, Ryan J. [Department of Atmospheric Scien↗
DOE OSTI2024

Probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under conditions of high ice-shelf basal melt

Abstract. Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter “the Amery sector”) to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1983

A coupled ice-ocean model of upwelling in the marginal ice zone

A dynamical coupled ice-ocean numerical model for the marginal ice zone (MIZ) is suggested and used to study upwelling dynamics in the MIZ. The nonlinear sea ice model has a variable ice concentration and includes internal ice stress. The model is forced by stresses on the air/ocean and air/ice surfaces. The main coupling between the ice and the ocean is in the form of an interfacial stress on the ice/ocean interface. The ocean model is a linear reduced gravity model. The wind stress exerted by the atmosphere on the ocean is proportional to the fraction of open water, while the interfacial stress ice/ocean is proportional to the concentration of ice. A new mechanism for ice edge upwelling is suggested based on a geostrophic equilibrium solution for the sea ice medium. The upwelling reported in previous models invoking a stationary ice cover is shown to be replaced by a weak downwelling due to the ice motion. Most of the upwelling dynamics can be understood by analysis of the divergence of the across ice edge upper ocean transport. On the basis of numerical model, an analytical model is suggested that reproduces most of the upwelling dynamics of the more complex numerical model.

Roed, L. P.↗
NASA NTRS1992

Effects of weather on the retrieval of sea ice concentration and ice type from passive microwave data

Effects of wind, water vapor, and cloud liquid water on ice concentration and ice type calculated from passive microwave data are assessed through radiative transfer calculations and observations. These weather effects can cause overestimates in ice concentration and more substantial underestimates in multi-year ice percentage by decreasing polarization and by decreasing the gradient between frequencies. The effect of surface temperature and air temperature on the magnitudes of weather-related errors is small for ice concentration and substantial for multiyear ice percentage. The existing weather filter in the NASA Team Algorithm addresses only weather effects over open ocean; the additional use of local open-ocean tie points and an alternative weather correction for the marginal ice zone can further reduce errors due to weather. Ice concentrations calculated using 37 versus 18 GHz data show little difference in total ice covered area, but greater differences in intermediate concentration classes. Given the magnitude of weather-related errors in ice classification from passive microwave data, corrections for weather effects may be necessary to detect small trends in ice covered area and ice type for climate studies.

Maslanik, J. A.↗
NASA NTRS1994

The effect of sea ice on the solar energy budget in the astmosphere-sea ice-ocean system: A model study

A coupled one-dimensional multilayer and multistream radiative transfer model has been developed and applied to the study of radiative interactions in the atmosphere, sea ice, and ocean system. The consistent solution of the radiative transfer equation in this coupled system automatically takes into account the refraction and reflection at the air-ice interface and allows flexibility in choice of stream numbers. The solar radiation spectrum (0.25 micron-4.0 micron) is divided into 24 spectral bands to account adequately for gaseous absorption in the atmosphere. The effects of ice property changes, including salinity and density variations, as well as of melt ponds and snow cover variations over the ice on the solar energy distribution in the entire system have been studied quantitatively. The results show that for bare ice it is the scattering, determined by air bubbles and brine pockets, in just a few centimeters of the top layer of ice that plays the most important role in the solar energy absorption and partitioning in the entire system. Ice thickness is important to the energy distribution only when the ice is thin, while the absorption in the atmosphere is not sensitive to ice thickness exceeds about 70 cm. The presence of clouds moderates all the sensitivities of the absorptive amounts in each layer to the variations in the ice properties and ice thickness. Comparisons with observational spectral albedo values for two simple ice types are also presented.

Jin, Z.↗
NASA NTRS2000

Ice Accretions and Icing Effects for Modern Airfoils

Icing tests were conducted to document ice shapes formed on three different two-dimensional airfoils and to study the effects of the accreted ice on aerodynamic performance. The models tested were representative of airfoil designs in current use for each of the commercial transport, business jet, and general aviation categories of aircraft. The models were subjected to a range of icing conditions in an icing wind tunnel. The conditions were selected primarily from the Federal Aviation Administration's Federal Aviation Regulations 25 Appendix C atmospheric icing conditions. A few large droplet icing conditions were included. To verify the aerodynamic performance measurements, molds were made of selected ice shapes formed in the icing tunnel. Castings of the ice were made from the molds and placed on a model in a dry, low-turbulence wind tunnel where precision aerodynamic performance measurements were made. Documentation of all the ice shapes and the aerodynamic performance measurements made during the icing tunnel tests is included in this report. Results from the dry, low-turbulence wind tunnel tests are also presented.

Addy, Harold E., Jr.↗
NASA NTRS2006

Recent Changes in High-Latitude Glaciers, Ice Caps, and Ice Sheets

The glaciers and ice sheets of the world contain enough ice to raise sea level by approximately 70 meters if they were to disappear entirely, and most of this ice is located in the climatically sensitive polar regions. Fortunately changes of this magnitude would probably take many thousands of years to occur, but recent discoveries indicate that these ice masses are responding to changes in today s climate more rapidly than previously thought. These responses are likely to be of great societal significance, primarily in terms of their implications for sea level, but also in terms of how their discharge of freshwater, through melting or calving, may impact ocean circulation. For millions of years, oceans have risen and fallen as the Earth has warmed and cooled, and ice on land has shrunk and grown. Today is no different in that respect, as sea levels have been rising at a rate of nearly 2 m per year during the last century (Miller and Douglas 2004), and 3 mm/yr in the last 12 years (Leuliette et al. 2004). What is different today, however, is that tens - perhaps hundreds - of millions of people live in coastal areas that are vulnerable to changes in sea level. Rising seas erode beaches, increase flood potential, and reduce the ability of barrier islands and coastal wetlands to mitigate the effects of major storms and hurricanes. The costs associated with a one-meter rise in sea level are estimated to be in the hundreds of billions of dollars in the United States alone. The worldwide costs in human terms would be far greater as some vulnerable low-lying coastal regions would become inundated, especially in poorer nations that do not have the resources to deal with such changes. Such considerations are particularly important in light of the fact that a one meter sea level rise is not significantly outside the 0.09 to 0.88 range of predictions for this century (IPCC 2001), and rises of this magnitude have occurred in the past in as little as 20 years (Fairbanks 1989). While the expansion of the warming oceans is estimated to be about a third of recent sea level rise, (Miller and Douglas 2004) the greatest potential for significantly increasing sea level lies in the Greenland and Antarctic ice sheets. For different reasons, each exhibits characteristics that suggest they are potentially unstable. In Antarctica, large portions of the ice cover rest on a soft bed that lies below sea level, making it vulnerable to runaway retreat. The Greenland ice sheet experiences considerable melt, which has the potential to rapidly accelerate the flow of ice toward the sea. While smaller ice masses, such as the Alaskan Glaciers and the Canadian ice caps, do not have anywhere near the same potential to impact sea level as the vast ice sheets do, many are melting rapidly, posing a significant near-term threat.

Abdalati, Waleed↗
NASA NTRS2011

Regional Changes in the Sea Ice Cover and Ice Production in the Antarctic

Coastal polynyas around the Antarctic continent have been regarded as sea ice factories because of high ice production rates in these regions. The observation of a positive trend in the extent of Antarctic sea ice during the satellite era has been intriguing in light of the observed rapid decline of the ice extent in the Arctic. The results of analysis of the time series of passive microwave data indicate large regional variability with the trends being strongly positive in the Ross Sea, strongly negative in the Bellingshausen/Amundsen Seas and close to zero in the other regions. The atmospheric circulation in the Antarctic is controlled mainly by the Southern Annular Mode (SAM) and the marginal ice zone around the continent shows an alternating pattern of advance and retreat suggesting the presence of a propagating wave (called Antarctic Circumpolar Wave) around the circumpolar region. The results of analysis of the passive microwave data suggest that the positive trend in the Antarctic sea ice cover could be caused primarily by enhanced ice production in the Ross Sea that may be associated with more persistent and larger coastal polynyas in the region. Over the Ross Sea shelf, analysis of sea ice drift data from 1992 to 2008 yields a positive rate-of-increase in the net ice export of about 30,000 km2 per year. For a characteristic ice thickness of 0.6 m, this yields a volume transport of about 20 km3/year, which is almost identical, within error bars, to our estimate of the trend in ice production. In addition to the possibility of changes in SAM, modeling studies have also indicated that the ozone hole may have a role in that it causes the deepening of the lows in the western Antarctic region thereby causing strong winds to occur offthe Ross-ice shelf.

Comiso, Josefino C.↗
NASA NTRS2017

Evaluation of Alternative Altitude Scaling Methods for Thermal Ice Protection System in NASA Icing Research Tunnel

A test was conducted at NASA Icing Research Tunnel to evaluate altitude scaling methods for thermal ice protection system. Two new scaling methods based on Weber number were compared against a method based on Reynolds number. The results generally agreed with the previous set of tests conducted in NRCC Altitude Icing Wind Tunnel where the three methods of scaling were also tested and compared along with reference (altitude) icing conditions. In those tests, the Weber number-based scaling methods yielded results much closer to those observed at the reference icing conditions than the Reynolds number-based icing conditions. The test in the NASA IRT used a much larger, asymmetric airfoil with an ice protection system that more closely resembled designs used in commercial aircraft. Following the trends observed during the AIWT tests, the Weber number based scaling methods resulted in smaller runback ice than the Reynolds number based scaling, and the ice formed farther upstream. The results show that the new Weber number based scaling methods, particularly the Weber number with water loading scaling, continue to show promise for ice protection system development and evaluation in atmospheric icing tunnels.

upwind schemes↗
NASA NTRS2019

Ice-Crystal Icing Accretion Studies at the NASA Propulsion Systems Laboratory

This paper describes an ice-crystal icing experiment conducted at the NASA Propulsion System Laboratory during June 2018. This test produced ice shape data on an airfoil for different test conditions similar to those inside the compressor region of a turbo-fan jet engine. Mixed-phase icing conditions were generated by partially freezing out a water spray using the relative humidity of flow as the primary parameter to control freeze-out. The paper presents the ice shape data and associated conditions which include pressure, velocity, temperature, humidity, total water content, melt ratio, and particle size distribution. The test featured a new instrument traversing system which allowed surveys of the flow and cloud. The purpose of this work was to provide experimental ice shape data and associated conditions to help develop and validate ice-crystal icing accretion models. The results support previous experimental observations of a minimum melt-ratio threshold for accretion to occur as well as the existence of a plateau region where the icing severity is high for a range of melt ratios. However, a maximum limit for melt ratio, which is suggested in the ice crystal icing literature, was not observed perhaps complicated by the potential for some supercooling of the water at these conditions.

Engine Icing↗
NASA NTRS2020

Analyzing the Impact of CryoSat-2 Ice Thickness Initialization on Seasonal Arctic Sea Ice Prediction

Twin 5-month seasonal forecast experiments are performed to predict the September 2018 minimum ice extent using the fully coupled Navy Earth System Prediction Capability (ESPC). In the control run, ensemble forecasts are initialized from the operational US Navy Global Ocean Forecasting System (GOFS) 3.1 for the ocean and sea ice but do not assimilate ice thickness data. Another set of forecasts are initialized from the same GOFS 3.1 fields but with sea ice thickness derived from CryoSat-2 (CS2). The Navy ESPC ensemble mean September 2018 minimum sea ice extent initialized with GOFS 3.1 ice thickness was over-predicted by 0.68 M sq.km (5.27 M sq.km) versus the ensemble set of forecasts initialized with CS2 ice thickness that had an error of 0.40 M sq.km (4.99 M sq.km), a 56% reduction in error. The September mean Integrated Ice Edge Error (IIEE) shows a 19% improvement for the entire Arctic with the CS2 data versus the control run. Comparison against Upward Looking Sonar (ULS) ice thickness in the Beaufort Sea reveals a lower bias and RMSE with the CS2 forecasts at all three moorings. Ice concentration at these locations is also improved, but neither set of experiments show ice free conditions as observed at moorings A and D.

time-lagged ensembles↗
NASA NTRSDate not supplied

Predictability of the Minimum Sea Ice Extent from Late Winter Fram Strait Ice Export: Model vs Observations

Late winter coastal divergence along the Eurasian coastline (referred to as the ice factory) – or the Fram Strait ice export, a proxy for coastal divergence in the ice factory – is a skillful predictor of the minimum sea ice extent (Williams et al., 2016; Brunette et al., 2018; Sesternikov and, 1979). Coastal divergence leads to the formation of coastal polynya where new ice grows but to a thickness that is not large enough to survive the following summer melt. This signal is then amplified by the ice albedo feedback and leads to more open water at the end of the summer melt season. In this thesis project, we will identify if this source of predictive skill in the seasonal forecast of the minimum sea ice extent is present in General Circulation Models, more precisely in the CESM2-LE, GISS-E2.1-G, GFDL FLOR-LE, CNRM-CM6-1 and the CanESM5. Since even small biases in the large-scale atmospheric circulation simulated by a model can short circuit this coupling between dominant modes of atmospheric variability (NAO and AO), coastal divergence along the Eurasian coastline, Fram Strait ice export and therefore seasonal forecasting skill of the model, failure to reproduce this coupling observed in the real Arctic will be used to identify biases in GCM. Preliminary results show that subtle changes in the large-scale atmospheric circulation leads to opposite statistical relationship between Fram Strait ice area export, coastal divergence along the Eurasian coastline, and seasonal predictability of the minimum sea ice extent both in the CESM2-LE and the GISS-E2.1-G models. Differences are linked with the partitioning between recirculation within the Beaufort Gyre, ice exported through Fram Strait and ridging north of the Canadian Arctic Archipelago.

Sea ice↗
NASA NTRSDate not supplied

Improved Treatment of Snow Over Ice Sheets in the NASA GISS Climate Model: Towards Ice Sheet–Climate Coupling

Representing the interactions between ice sheets and climate is essential for more accurate prediction of climate change and sea level rise. Ice sheets interact with the overlying atmosphere via the accumulation of snow and its compaction into firn, then ice, as well as the melting of surface snow and ice and the creation of runoff water. Getting an adequate representation of heat transfer, compaction, and melting processes is essential for an accurate representation of snow on land ice in global climate models. We are implementing an improved snow model on top of land ice as part of an effort to couple the NASA GISS climate model with the PISM ice sheet model. The new snow model includes additional layers and processes that are not currently incorporated (e.g., liquid water retention, percolation and refreezing, and snow densification), and mass and energy transfer methods that are consistent with both static ice sheets (with implicit iceberg fluxes) and interactive ice sheets (with explicit dynamics). We are tuning the densification scheme of this snow model with temperature and density data from common FirnCover and SumUp observations at locations in the accumulation zone of Greenland, and we compare the resulting density profiles to other SumUp density profiles in Greenland and Antarctica. We will assess the impact of this new snow model in climate model simulations with a static ice sheet compared with the previous (simpler) 2-layer snow model. Finally, we aim to use the non-coupled simulations as a baseline to assess the impact of dynamic coupling with an interactive ice sheet model.

snow↗