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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS2011

Vostok Subglacial Lake: A Review of Geophysical Data Regarding Its Discovery and Topographic Setting

Vostok Subglacial Lake is the largest and best known sub-ice lake in Antarctica. The establishment of its water depth (>500 m) led to an appreciation that such environments may be habitats for life and could contain ancient records of ice sheet change, which catalyzed plans for exploration and research. Here we discuss geophysical data used to identify the lake and the likely physical, chemical, and biological processes that occur in it. The lake is more than 250 km long and around 80 km wide in one place. It lies beneath 4.2 to 3.7 km of ice and exists because background levels of geothermal heating are sufficient to warm the ice base to the pressure melting value. Seismic and gravity measurements show the lake has two distinct basins. The Vostok ice core extracted >200 m of ice accreted from the lake to the ice sheet base. Analysis of this ice has given valuable insights into the lake s biological and chemical setting. The inclination of the ice-water interface leads to differential basal melting in the north versus freezing in the south, which excites circulation and potential mixing of the water. The exact nature of circulation depends on hydrochemical properties, which are not known at this stage. The age of the subglacial lake is likely to be as old as the ice sheet (approx.14 Ma). The age of the water within the lake will be related to the age of the ice melting into it and the level of mixing. Rough estimates put that combined age as approx.1 Ma.

Siegert, Martin J.↗

NASA NTRS2026

Satellite Optical Remote Sensing of Clouds and Aerosols: From Particle Single-Scattering and Gaseous Absorption Through Radiative Transfer to Retrieval Products

Clouds and aerosols are fundamental regulators of Earth’s radiation budget and climate system, influencing both solar and terrestrial radiation through scattering, absorption, and emission processes. Accurate characterization of their physical and radiative properties from space requires a rigorous understanding of particle single-scattering, gaseous absorption, and radiative transfer in the atmosphere, as well as reliable inversion methods. This review synthesizes the physical foundations and algorithmic implementations of satellite-based passive optical remote sensing of clouds and aerosols, spanning the ultraviolet to thermal infrared spectral range. Beginning with electromagnetic scattering theory and state-of-the-art methods for computing single-scattering by nonspherical particles and computationally efficient methods for accounting for atmospheric absorption, we discuss the radiative transfer framework underpinning cloud and aerosol retrievals. The connection between single-scattering and multiple-scattering is rigorously formulated. We then summarize operational and research-grade retrieval techniques, including cloud masking and thermodynamic phase determination, CO₂ slicing for cloud-top pressure, the Nakajima-King shortwave bi-spectral, and infrared split-window approaches for cloud optical thickness and effective particle size, inversion algorithms for determining aerosol properties from multi-spectral and/or multi-angle radiometric and polarimetric measurements, and active-passive sensing synergy. Examples of the global cloud and aerosol climatologies are illustrated using observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Multi-angle Imaging SpectroRadiometer (MISR). Furthermore, the unique strengths of active remote sensing techniques based on spaceborne lidar observations are briefly elaborated in the context of studying ice clouds composed of randomly and horizontally oriented ice crystals, which is a significant challenge for conventional passive remote sensing techniques. By connecting physical theory to practical retrievals, this review highlights both the maturity of current methodologies and the remaining challenges in reducing uncertainties in particle morphology, vertical structure, absorption, and aerosol-cloud interactions. Furthermore, the impact of artificial intelligence (AI) on atmospheric remote sensing is briefly addressed.

Aerosols↗

NASA NTRS2014

An Improved Cryosat-2 Sea Ice Freeboard Retrieval Algorithm Through the Use of Waveform Fitting

We develop an empirical model capable of simulating the mean echo power cross product of CryoSat-2 SAR and SAR In mode waveforms over sea ice covered regions. The model simulations are used to show the importance of variations in the radar backscatter coefficient with incidence angle and surface roughness for the retrieval of surfaceelevation of both sea ice floes and leads. The numerical model is used to fit CryoSat-2 waveforms to enable retrieval of surface elevation through the use of look-up tables and a bounded trust region Newton least squares fitting approach. The use of a model to fit returns from sea ice regions offers advantages over currently used threshold retrackingmethods which are here shown to be sensitive to the combined effect of bandwidth limited range resolution and surface roughness variations. Laxon et al. (2013) have compared ice thickness results from CryoSat-2 and IceBridge, and found good agreement, however consistent assumptions about the snow depth and density of sea ice werenot used in the comparisons. To address this issue, we directly compare ice freeboard and thickness retrievals from the waveform fitting and threshold tracker methods of CryoSat-2 to Operation IceBridge data using a consistent set of parameterizations. For three IceBridge campaign periods from March 20112013, mean differences (CryoSat-2 IceBridge) of 0.144m and 1.351m are respectively found between the freeboard and thickness retrievals using a 50 sea ice floe threshold retracker, while mean differences of 0.019m and 0.182m are found when using the waveform fitting method. This suggests the waveform fitting technique is capable of better reconciling the seaice thickness data record from laser and radar altimetry data sets through the usage of consistent physical assumptions.

CryoSat-2↗

NASA NTRS2004

The Search for Young Planetary Systems And the Evolution of Young Stars

The Space Interferometer Mission (SIM) will provide a census of planetary systems by con- ducting a broad survey of 2,000 stars that will be sensitive to the presence of planets with masses as small as approx. 15 Earth masses (1 Uranus mass) and a deep survey of approx. 250 of the nearest, stars with a mass limit of approx.3 Earth masses. The broad survey will include stars spanning a wide range of ages, spectral types, metallicity, and other important parameters. Within this larger context, the Young Stars and Planets Key Project will study approx. 200 stars with ages from 1 Myr to 100 Myr to understand the formation and dynamical evolution of gas giant planets. The SIM Young Stars and Planets Project will investigate both the frequency of giant planet formation and the early dynamical history of planetary systems. We will gain insight into how common the basic architecture of our solar system is compared with recently discovered systems with close-in giant planets by examining 200 of the nearest (less than 150 pc) and youngest (1-100 Myr) solar-type stars for planets. The sensitivity of the survey for stars located 140 pc away is shown in the planet mass-separation plane. We expect to find anywhere from 10 (assuming that only the presently known fraction of stars. 5-7%, has planets) to 200 (all young stars have planets) planetary systems. W-e have set our sensitivity threshold to ensure the detection of Jupiter-mass planets in the critical orbital range of 1 to 5 AU. These observations, when combined with the results of planetary searches of mature stars, will allow us to test theories of planetary formation and early solar system evolution. By searching for planets around pre-main sequence stars carefully selected to span an age range from 1 to 100 Myr, we will learn a t what epoch and with what frequency giant planets are found at the water-ice snowline where they are expected to form. This will provide insight into the physical mechanisms by which planets form and migrate from their place of birth, and about their survival rate. With these data in hand, we will provide data, for the first time, on such important questions as: What processes affect the formation and dynamical evolution of planets? When and where do planets form? What is initial mass distribution of planetary systems around young stars? How might planets be destroyed? What is the origin of the eccentricity of planetary orbits? What is the origin of the apparent dearth of companion objects between planets and brown dwarfs seen in mature stars? The observational strategy is a compromise between the desire to extend the planetary mass function as low as possible and the essential need to build up sufficient statistics on planetary occurrence. About half of the sample will be used to address the "where" and "when" of planet formation. We will study classical T Tauri stars (cTTs) which have massive accretion disks and post- accretion, weak-lined T Tauri stars (wTTs). Preliminary estimates suggest the sample will consist of approx. 30% cTTs and approx. 70% wTTs, driven in part by the difficulty of making accurate astrometric measurements toward objects with strong variability or prominent disks.

Beichman, Charles A.↗

NASA NTRSDate not supplied

Age Estimates for Permanently Shadowed Craters in the VIPER Mission Area Based On Their Topography

A primary objective of the VIPER [1] mission is to characterize the distribution and physical state of volatiles at the lunar poles, including within permanently shadowed regions (PSRs) where water ice has been inferred to be stable [e.g.. 2,3]. A mission area for VIPER has been defined that enables this scientific objective near Nobile crater (Fig. 1)[4]. This location enables a traverse that can both meet VIPER’s engineering constraints (Earth-direct communication, adequate power, etc.)as well as accomplish the planned scientific exploration. In this abstract, we describe observations of crater topography that provide insight into the age of several craters that host PSRs within the planned VIPER mission area. The role that the age of PSRs plays in controlling the presence or absence of polar volatiles is of substantial interest for discerning volatile history [e.g., 5-7]. The physical state, depth distribution, and spatial distribution of volatile deposits may also vary as a function of PSR age due to gardening and/or differing emplacement mechanisms [8]. Understanding the age of PSRs that VIPER may explore is thus a useful goal.

C I Fassett↗

NASA NTRS2009

Radar Differential Phase Signatures of Ice Orientation for the Prediction of Lightning Initiation and Cessation

The majority of lightning-related casualties typically occur during thunderstorm initiation (e.g., first flash) or dissipation (e.g., last flash). The physics of electrification and lightning production during thunderstorm initiation is fairly well understood. As such, the literature includes a number of studies presenting various radar techniques (using reflectivity and, if available, other dual-polarimetric parameters) for the anticipation of initial electrification and first lightning flash. These radar techniques have shown considerable skill at forecasting first flash. On the other hand, electrical processes and lightning production during thunderstorm dissipation are not nearly as well understood and few, if any, successful techniques have been developed to anticipate the last flash and subsequent cessation of lightning. One promising approach involves the use of dual-polarimetric radar variables to infer the presence of oriented ice crystals in lightning producing storms. In the absence of strong vertical electric fields, ice crystals fall with their largest (semi-major) axis in the horizontal associated with gravitational and aerodynamic forces. In thunderstorms, strong vertical electric fields (100-200 kV m(sup -1)) have been shown to orient small (less than 2 mm) ice crystals such that their semi-major axis is vertical (or nearly vertical). After a lightning flash, the electric field is typically relaxed and prior radar research suggests that ice crystals rapidly resume their preferred horizontal orientation. In active thunderstorms, the vertical electric field quickly recovers and the ice crystals repeat this cycle of orientation for each nearby flash. This change in ice crystal orientation from primarily horizontal to vertical during the development of strong vertical electric fields prior to a lightning flash forms the physical basis for anticipating lightning initiation and, potentially, cessation. Research has shown that radar reflectivity (Z) and other co-polar back-scattering radar measurements like differential reflectivity (Z(sub dr)) typically measured by operational dual-polarimetric radars are not sensitive to these changes in ice crystal orientation. However, prior research has demonstrated that oriented ice crystals cause significant propagation effects that can be routinely measured by most dual-polarimetric radars from X-band (3 cm) to S-band (10 cm) wavelengths using the differential propagation phase shift (often just called differential phase, phi(sub dp)) or its range derivative, the specific differential phase (K(sub dp)). Advantages of the differential phase include independence from absolute or relative power calibration, attenuation, differential attenuation and relative insensitivity to ground clutter and partial beam occultation effects (as long as the signal remains above noise). In research mode, these sorts of techniques have been used to anticipate initial cloud electrification, lightning initiation, and cessation. In this study, we develop a simplified model of ice crystal size, shape, orientation, dielectric, and associated radar scattering and propagation effects in order to simulate various idealized scenarios of ice crystals responding to a hypothetical electric field and their dual-polarimetric radar signatures leading up to lightning initiation and particularly cessation. The sensitivity of the K(sub dp) ice orientation signature to various ice properties and radar wavelength will be explored. Since K(sub dp) is proportional to frequency in the Rayleigh- Gans scattering regime, the ice orientation signatures should be more obvious at higher (lower) frequencies (wavelengths). As a result, simulations at radar wavelengths from 10 cm down to 1 cm (Ka-band) will be conducted. Resonance effects will be considered using the T-matrix method. Since most K(sub dp) Vbased observations have been shown at S-band, we will present ice orientation signatures from C-band (UAH/NASA ARMOR) and X-bd (UAH MAX) dual-polarimetric radars located in Northern Alabama. Issues related to optimal radar scanning for the detection of oriented ice will be discussed. Preliminary suggestions on how these differential phase signatures of oriented ice could contribute to lightning initiation and cessation algorithms will be presented.

Carey, L.D.↗

NASA NTRSDate not supplied

VIPER: Systems Integration Status

NASA’s Artemis Program plans to return humans to the Moon for an extended stay. To do so will require substantial resources to sustain that continued human presence, including continuous supplies delivered from the Earth. Given the expense and complexity of resource deliveries from Earth, if some resources were indigenously available, substantial logistical savings could be available by “living off the land”, wherever possible. The LCROSS[1] , LRO and other missions have confirmed the presence of resources such as water-ice and other volatiles in lunar polar regions, so the next step is to understand the scientific nature and physical distribution of those candidate resources. Those local volatiles could be processed into propellants and human life-supporting needs, reducing risk of maintaining a permanent human presence on the Moon. The Volatiles Investigating Polar Exploration Resource (VIPER) is a surface mobility scientific platform, designed to spend ~100 days mapping and surveying four different Ice Stability Regions to understand the scientific nature and distribution of water and other volatiles. VIPER will also provide scientific mineralogical context of the lunar regolith, such as the presence of silicon and light metals in lunar regolith, providing a composite picture of resource availability and sustainment. This paper will discuss the latest systems-level integration activities by the VIPER team, following our initial introduction to this mission at IAC2021[2] . The VIPER team successfully passed its Systems Integration Review (SIR) in late-2022, and in early 2023, began system-level surface segment (rover) flight hardware assembly. VIPER is managed within NASA’s Science Mission Directorate (SMD), utilizing the Commercial Lunar Payload Services (CLPS) lunar delivery model with partner, Astrobotic, Inc.

VIPER↗

NASA NTRS2015

Model Forecast Skill and Sensitivity to Initial Conditions in the Seasonal Sea Ice Outlook

We explore the skill of predictions of September Arctic sea ice extent from dynamical models participating in the Sea Ice Outlook (SIO). Forecasts submitted in August, at roughly 2 month lead times, are skillful. However, skill is lower in forecasts submitted to SIO, which began in 2008, than in hindcasts (retrospective forecasts) of the last few decades. The multimodel mean SIO predictions offer slightly higher skill than the single-model SIO predictions, but neither beats a damped persistence forecast at longer than 2 month lead times. The models are largely unsuccessful at predicting each other, indicating a large difference in model physics and/or initial conditions. Motivated by this, we perform an initial condition sensitivity experiment with four SIO models, applying a fixed −1 m perturbation to the initial sea ice thickness. The significant range of the response among the models suggests that different model physics make a significant contribution to forecast uncertainty.

Arctic Sea↗

NASA NTRS2011

Understanding Effective Diameter and Its Application to Terrestrial Radiation in Ice Clouds

The cloud property known as "effective diameter" or "effective radius", which in essence is the cloud particle size distribution (PSD) volume at bulk density divided by its projected area, is used extensively in atmospheric radiation transfer, climate modeling and remote sensing. This derives from the assumption that PSD optical properties can be uniquely described in terms of their effective diameter, D(sub e), and their cloud water content (CWC), henceforth referred to as the D(sub e)-CWC assumption. This study challenges this assumption, showing that while the D(sub e)-CWC assumption appears generally valid for liquid water clouds, it appears less valid for ice clouds in regions where (1) absorption is not primarily a function of either the PSD ice water content (IWC) or the PSD projected area, and (2) where wave resonance (i.e. photon tunneling) contributes significantly to absorption. These two regions often strongly coincide at terrestrial wavelengths when De less than ~60 m, which is where this D(sub e)-CWC assumption appears poorest. Treating optical properties solely in terms of D(sub e) and IWC may lead to errors up to 24%, 26% and 20% for terrestrial radiation in the window region regarding the absorption and extinction coefficients and the single scattering albedo, respectively. Outside the window region, errors may reach 33% and 42% regarding absorption and extinction. The magnitude and sign of these errors can change rapidly with wavelength, which may produce significant errors in climate modeling, remote sensing and other applications concerned with the wavelength dependence of radiation. Where the D(sub e)-CWC assumption breaks down, ice cloud optical properties appear to depend on D(sub e), IWC and the PSD shape. Optical property parameterizations in climate models and remote sensing algorithms based on historical PSD measurements may exhibit errors due to previously unknown PSD errors (i.e. the presence of ice artifacts due to the shattering of larger ice particles on the probe inlet tube during sampling). More recently developed cloud probes are designed to mitigate this shattering problem. Using realistic PSD shapes for a given temperature (and/or IWC) and cloud type may minimize errors associated with PSD shape in ice optics parameterizations and remote sensing algorithms. While this topic was investigated using two ice optics schemes (the Yang et al., 2005 database and the modified anomalous diffraction approximation, or MADA), a physical understanding of the limitations of the D(sub e)-IWC assumption was made possible by using MADA. MADA allows one to approximate the contribution of photon tunneling to absorption relative to other optical processes, which reveals that part of the error regarding the D(sub e)-IWC assumption can be associated with tunneling. By relating the remaining error to the radiation penetration depth in bulk ice (DELTA L) due to absorption, the domain where the D(sub e)-IWC assumption is weakest was described in terms of D(sub e) and DELTA L.

Mitchell, D. L.↗

NASA NTRS1994

H2-rich interstellar grain mantles: An equilibrium description

Experiments simulating the codeposition of molecular hydrogen and water ice on interstellar grains demonstrate that amorphous water ice at 12 K can incorporate a substantial amount of H2, up to a mole ratio of H2/H2O = 0.53. We find that the physical behavior of approximately 80% of the hydrogen can be explained satisfactorily in terms of an equilibrium population, thermodynamically governed by a wide distribution of binding site energies. Such a description predicts that gas phase accretion could lead to mole fractions of H2 in interstellar grain mantles of nearly 0.3; for the probable conditions of WL5 in the rho Ophiuchi cloud, an H2 mole fraction of between 0.05 and 0.3 is predicted, in possible agreement with the observed abundance reported by Sandford, Allamandola, & Geballe. Accretion of gas phase H2 onto grain mantles, rather than photochemical production of H2 within the ice, could be a general explanation for frozen H2 in interstellar ices. We speculate on the implications of such a composition for grain mantle chemistry and physics.

Dissly, Richard W.↗

NASA NTRS1998

35 GHz Measurements of CO2 Crystals for Simulating Observations of the Martian Polar Caps

In order to learn more about the Martian polar caps, it is important to compare and contrast the behavior of both frozen H2O and CO2 in different parts of the electromagnetic spectrum. Relatively little attention has been given, thus far, to observing the thermal microwave part of the spectrum. In this experiment, passive microwave radiation emanating from within a 33 cm snowpack was measured with a 35 GHz hand-held radiometer, and in addition to the natural snow measurements, the radiometer was used to measure the microwave emission and scattering from layers of manufactured CO2 (dry ice). A 1 m x 2 m plate of aluminum sheet metal was positioned beneath the natural snow so that microwave emissions from the underlying soil layers would be minimized. Compared to the natural snow crystals, results for the dry ice layers exhibit lower' microwave brightness temperatures for similar thicknesses, regardless of the incidence angle of the radiometer. For example, at 50 degree H (horizontal polarization) and with a covering of 21 cm of snow and 18 cm of dry ice, the brightness temperatures were 150 K and 76 K, respectively. When the snow depth was 33 cm, the brightness temperature was 144 K, and when the total thickness of the dry ice was 27 cm, the brightness temperature was 86 K. The lower brightness temperatures are due to a combination of the lower physical temperature and the larger crystal sizes of the commercial CO2 Crystals compared to the snow crystals. As the crystal size approaches the size of the microwave wavelength, it scatters microwave radiation more effectively, thus lowering the brightness temperature. The dry ice crystals in this experiment were about an order of magnitude larger than the snow crystals and three orders of magnitude larger than the CO2 Crystals produced in the cold stage of a scanning electron microscope. Spreading soil, approximately 2 mm in thickness, on the dry ice appeared to have no effect on the brightness temperatures.

Foster, J. L.↗

NASA NTRS2002

Simulating Snow Over Sea Ice In Climate Models

We have evaluated two methods of simulating the seasonal cycle of snow over sea ice in and around the Arctic: The NCAR global climate model CCM3, with its standard snow hydrology, and the snow pack model SNTHERM, forced with hourly atmospheric output from CCM3. A new dataset providing dates for the onset of snow melt over Arctic sea ice provides a means for assessing basin-wide how well the models simulate melt onset, but contains no information on how long it then takes for all the snow to melt. Use of data from the SHEBA site provides very detailed information on the behavior of the snow before and during the melt season, but only for a very limited area. Russian drift data provide climatological data on the seasonal cycle of snow water equivalent and snow density, over multi-year sea ice in the central Arctic basin. These datasets are used to compare the two modeling methods, and to see if use of the more physically-realistic SNTHERM provides any significant improvements. Conclusions obtained so far include: 1. Both CCM3 and CCM3/SNTHERM do a good job overall of matching the onset of snow melt dataset; although CCM3/SNTHERM consistently trends to underestimate the date and CCM3 to overestimate it. 2. SHEBA and ice drift data for the Arctic show that CCM3/ SNTHERM does a better job than CCM3 at simulating the total melt period. 3. Ice drift snow density and accumulation data suggest that while providing superior results, CCM3/SNTHERM may still suffer from overly vigorous melting. 4. Both the large-scale atmospheric forcing and snow pack physical processes are important in proper simulation of the snow seasonal cycle. Ongoing work includes further diagnosis of CCM3/SNTHERM, use of more observational datasets, especially from marginal seas in the pan-Arctic, and full coupling of SNTHERM into CCM3 (work to date has all been off-line simulations).

Arnold, James E.↗

NASA NTRS1985

Progress in development of a Navier-Stokes solver for evaluation of iced airfoil performance

A method is being developed for evaluation of the flow field behavior about an airfoil with significant ice accretion on the leading edge. The computer code, being evaluated for this purpose, solves the Navier-Stokes equations in a body-fitted curvilinear coordinate system. This requires the use of a grid generation code to transform the x-y coordinates of the physical space into xi-eta coordinates of the computational space. Evaluation of the suitability of these two codes for predicting iced airfoil performance is presently being carried out in anticipation of use in an overall icing analysis effort. Results of this evaluation to date indicate good correlation with known information on clean airfoils. Preliminary results for rime and glaze, iced airfoil shapes are also presented.

Potapczuk, M. G.↗

NASA NTRS2002

[Equipment Inventory Report of Inventions and Subcontracts]

This was a joint award to the PI, then of Texas A&M University, and Gerald Heymsfield of NASA GSFC, with co-I's Robbie Hood and Richard Blakeslee of NASA Marshall Space Flight Center (MSFC). Although the PI moved to the University of Utah during year-2 of this grant, it continued to be administered through Texas A&M for the convenience of the Ph.D. student Daniel Cecil until after he defended his dissertation in late 2000. Cecil has assembled a subset of this database for 261 Tropical Rainfall Measuring Mission (TRMM) passes over 45 hurricanes or tropical cyclones. It consists of TRMM Microwave Imager (TMI) ice scattering signatures and TRMM Precipitation Radar (PR) reflectivity profiles and Lightning Imaging Sensor (LIS) flash locations for the rain features associated with all tropical cyclones observed by TRMM from Dec. 1997 through Dec. 1998. Each rain feature has been subjectively cataloged as belonging to a hurricane eyewall, inner rainband, outer rainband, or the inner or outer regions of a developing or non-developing tropical cyclone. This database quantifies the relative abundance of lightning in outer rainbands compared to eyewalls and inner rainbands. The TRMM data set permits studies that can seek more specific physical relationships than Cecil and Zipser were able to find when they related some SSM/I measures of ice scattering to hurricane intensity change.

Zipser, Edward J.↗

NASA NTRS2021

Physics-Based SNOWPACK Model Improves Representation of Near-Surface Antarctic Snow and Firn Density

Estimates of snow and firn density are required for satellite-altimetry-based retrievals of ice sheet mass balance that rely on volume-to-mass conversions. Therefore, biases and errors in presently used density models confound assessments of ice sheet mass balance and by extension ice sheet contribution to sea level rise. Despite this importance, most contemporary firn densification models rely on simplified semi-empirical methods, which are partially reflected by significant modeled density errors when compared to observations. In this study, we present a new drifting-snow compaction scheme that we have implemented into SNOWPACK, a physics-based land surface snow model. We show that our new scheme improves existing versions of SNOWPACK by increasing simulated near-surface (defined as the top 10 m) density to be more in line with observations (near-surface bias reduction from −44.9 to −5.4 kg m−3). Furthermore, we demonstrate high-quality simulation of near-surface Antarctic snow and firn density at 122 observed density profiles across the Antarctic ice sheet, as indicated by reduced model biases throughout most of the near-surface firn column when compared to two semi-empirical firn densification models (SNOWPACK mean bias=−9.7 kg m−3, IMAU-FDM mean bias=−32.5 kg m−3, GSFC-FDM mean bias=15.5 kg m−3). Notably, our analysis is restricted to the near surface where firn density is most variable due to accumulation and compaction variability driven by synoptic weather and seasonal climate variability. Additionally, the GSFC-FDM exhibits lower mean density bias from 7–10 m (SNOWPACK bias=−22.5 kg m−3, GSFC-FDM bias=10.6 kg m−3) and throughout the entire near surface at high-accumulation sites (SNOWPACK bias=−31.4 kg m−3, GSFC-FDM bias=−4.7 kg m−3). However, we found that the performance of SNOWPACK did not degrade when applied to sites that were not included in the calibration of semi-empirical models. This suggests that SNOWPACK may possibly better represent firn properties in locations without extensive observations and under future climate scenarios, when firn properties are expected to diverge from their present state.

Eric Keenan↗

NASA NTRS2023

A 1D Model for Nucleation of Ice From Aerosol Particles: An Application to a Mixed-Phase Arctic Stratus Cloud Layer

Mixed-phase clouds (MPCs) have been identified as significant contributors to uncertainties in climate projections, attributable to model representation of processes controlling the formation and loss of supercooled water droplets and ice particles from the atmosphere. Arctic MPCs are commonly widespread and long-lived, with sustained ice crystal formation processes that challenge current understanding. This study examines the ice-nucleating particle (INP) reservoir dynamics governing immersion-mode heterogeneous freezing in an observed case of Arctic MPCs using a simplified 1D aerosol-cloud model. The model setup includes prescribed dynamical forcings and thermodynamic profiles, and represents INPs as multicomponent and polydisperse particle size distributions. Diagnostic and prognostic approaches to immersion freezing parameterization are compared, including time-independent (singular) number- and surface area-based descriptions and a time-dependent description following classical nucleation theory (CNT). The choice of freezing parameterization defines the size of the INP reservoir. The CNT-based description yields an orders of magnitude larger INP reservoir than the singular parameterizations, which is the dominant factor for sustained ice crystal formation. The efficiency of the freezing process and cloud cooling are of secondary importance. A diagnostic treatment neglecting INP loss is only accurate when the INP reservoir size is large and INP depletion weak. Since a larger INP reservoir sustains ice crystal formation substantially longer, and ice water path scales with ice crystal concentrations for the conditions considered, resolving the source of differences in INP reservoir dynamics due to model implementation is a high priority for advancing climate model physics.

aerosol↗

NASA NTRS2005

Phytoplankton in the Beaufort and Chukchi Seas: Distributions, Dynamics and Environmental Forcing

Time-series of remotely sensed distributions of phytoplankton, sea ice, surface temperature, albedo, and clouds were examined to evaluate the impact of the variability of environmental conditions and physical forcing on the phytoplankton distribution in the Beaufort and Chukchi Seas. Large-scale distributions of these parameters were studied for the first time using weekly and monthly composites from April 1998 through September 2002. The basic data set used in this study are phytoplankton pigment concentration derived from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS), ice concentration obtained from the Special Sensor Microwave Imager (SSM/I) and surface temperature, cloud cover, and albedo derived from the Advanced Very High Resolution Radiometer (AVHRR). Seasonal variations of the sea ice cover was observed to be the dominant environmental factor as the ice edge blooms followed the retreating marginal ice zones northward. Blooms were most prominent in the southwestern Chukchi Sea, and were especially persistent immediately north of the Bering Strait in nutrient- rich Anadyr water and in some fronts. Chlorophyll concentrations are shown to increase from a nominal value during onset of melt in April to a maximum value in mid-spring or summer depending on location. Large interannual variability of ice cover and phytoplankton distributions was observed with the year 1998 being uniquely associated with an early season occurrence of a massive bloom. This is postulated to be caused in part by a rapid response of phytoplankton to an early retreat of the sea ice cover in the Beaufort Sea region. Correlation analyses showed relatively high negative correlation between chlorophyll and ice concentration with the correlation being highest in May, the correlation coefficient being -0.45. 1998 was also the warmest among the five years globally and the sea ice cover was least extensive in the Beaufort-Khukchi Sea region, partly because of the 1997-98 El Nino. Strong correlations were noted between ice extent and surface temperature, the correlation coefficient being highest at - 0.79 in April, during the onset of the bloom period

Wang, Jian↗

NASA NTRS2005

Neptune Polar Orbiter with Probes

The giant planets of the outer solar system divide into two distinct classes: the gas giants Jupiter and Saturn, which consist mainly of hydrogen and helium; and the ice giants Uranus and Neptune, which are believed to contain significant amounts of the heavier elements oxygen, nitrogen, and carbon and sulfur. Detailed comparisons of the internal structures and compositions of the gas giants with those of the ice giants will yield valuable insights into the processes that formed the solar system and, perhaps, other planetary systems. By 2012, Galileo, Cassini and possibly a Jupiter Orbiter mission with microwave radiometers, Juno, in the New Frontiers program, will have yielded significant information on the chemical and physical properties of Jupiter and Saturn. A Neptune Orbiter with Probes (NOP) mission would deliver the corresponding key data for an ice giant planet. Such a mission would ideally study the deep Neptune atmosphere to pressures approaching and possibly exceeding 1000 bars, as well as the rings, Triton, Nereid, and Neptune s other icy satellites. A potential source of power would be nuclear electric propulsion (NEP). Such an ambitious mission requires that a number of technical issues be investigated, however, including: (1) atmospheric entry probe thermal protection system (TPS) design, (2) probe structural design including seals, windows, penetrations and pressure vessel, (3) digital, RF subsystem, and overall communication link design for long term operation in the very extreme environment of Neptune's deep atmosphere, (4) trajectory design allowing probe release on a trajectory to impact Neptune while allowing the spacecraft to achieve a polar orbit of Neptune, (5) and finally the suite of science instruments enabled by the probe technology to explore the depths of the Neptune atmosphere. Another driving factor in the design of the Orbiter and Probes is the necessity to maintain a fully operational flight system during the lengthy transit time from launch through Neptune encounter, and throughout the mission. Following our response to the recent NASA Research Announcement (NRA) for Space Science Vision Missions for mission studies by NASA for implementation in the 2013 or later time frame, our team has been selected to explore the feasibility of such a Neptune mission.

Bienstock, Bernard↗