Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ice physics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 559 records · Page 31
All sourcesNASA NTRSDOE 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.↗
DOE OSTI2024

The MOSAiC Distributed Network: Observing the coupled Arctic system with multidisciplinary, coordinated platforms

Central Arctic properties and processes are important to the regional and global coupled climate system. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) Distributed Network (DN) of autonomous ice-tethered systems aimed to bridge gaps in our understanding of temporal and spatial scales, in particular with respect to the resolution of Earth system models. By characterizing variability around local measurements made at a Central Observatory, the DN covers both the coupled system interactions involving the ocean-ice-atmosphere interfaces as well as three-dimensional processes in the ocean, sea ice, and atmosphere. The more than 200 autonomous instruments (“buoys”) were of varying complexity and set up at different sites mostly within 50 km of the Central Observatory. During an exemplary midwinter month, the DN observations captured the spatial variability of atmospheric processes on sub-monthly time scales, but less so for monthly means. They show significant variability in snow depth and ice thickness, and provide a temporally and spatially resolved characterization of ice motion and deformation, showing coherency at the DN scale but less at smaller spatial scales. Ocean data show the background gradient across the DN as well as spatially dependent time variability due to local mixed layer sub-mesoscale and mesoscale processes, influenced by a variable ice cover. The second case (May–June 2020) illustrates the utility of the DN during the absence of manually obtained data by providing continuity of physical and biological observations during this key transitional period. We show examples of synergies between the extensive MOSAiC remote sensing observations and numerical modeling, such as estimating the skill of ice drift forecasts and evaluating coupled system modeling. The MOSAiC DN has been proven to enable analysis of local to mesoscale processes in the coupled atmosphere-ice-ocean system and has the potential to improve model parameterizations of important, unresolved processes in the future.

54 ENVIRONMENTAL SCIENCES↗
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↗
DOE OSTI2023

Beyond the Stokes approximation: shallow visco-elastic ice-sheet models

Abstract The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.

58 GEOSCIENCES↗
DOE OSTI2022

Evidence for superionic H 2 O and diffusive He–H 2 O at high temperature and high pressure

In this work, we present the evidence of superionic phase formed in H 2 O and, for the first time, diffusive H 2 O–He phase, based on time-resolved x-ray diffraction experiments performed on ramp-laser-heated samples in diamond anvil cells. The diffraction results signify a similar bcc-like structure of superionic H 2 O and diffusive He–H 2 O, while following different transition dynamics. Based on time and temperature evolution of the lattice parameter, the superionic H 2 O phase forms gradually in pure H 2 O over the temperature range of 1350–1400 K at 23 GPa, but the diffusive He–H 2 O phase forms abruptly at 1300 K at 26 GPa. We suggest that the faster dynamics and lower transition temperature in He–H 2 O are due to a larger diffusion coefficient of interstitial-filled He than that of more strongly bound H atoms. This conjecture is then consistent with He disordered diffusive phase predicted at lower temperatures, rather than H-disordered superionic phase in He–H 2 O.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Quantum Spin Ice in Three-Dimensional Rydberg Atom Arrays

Quantum spin liquids are exotic phases of matter whose low-energy physics is described as the deconfined phase of an emergent gauge theory. With recent theory proposals and an experiment showing preliminary signs of Z 2 topological order [G. Semeghini , ], Rydberg atom arrays have emerged as a promising platform to realize a quantum spin liquid. In this work, we propose a way to realize a U(1) quantum spin liquid in three spatial dimensions, described by the deconfined phase of U(1) gauge theory in a pyrochlore lattice Rydberg atom array. We study the ground state phase diagram of the proposed Rydberg system as a function of experimentally relevant parameters. Within our calculation, we find that by tuning the Rabi frequency, one can access both the confinement-deconfinement transition driven by a proliferation of “magnetic” monopoles and the Higgs transition driven by a proliferation of “electric” charges of the emergent gauge theory. We suggest experimental probes for distinguishing the deconfined phase from ordered phases. This work serves as a proposal to access a confinement-deconfinement transition in three spatial dimensions on a Rydberg-based quantum simulator. Published by the American Physical Society 2025

Shah, Jeet (ORCID:0000000158738129)↗
DOE OSTI2021

Early detection and uncertainty quantification of rapid sea-level rise from Antarctica

Focal Areas: (i) classification and/or anomaly detection methods applied towards identifying critical climate system thresholds; improved emulator design for (ii) reducing computational costs relative to full-physics models and (iii) improved uncertainty quantification workflows. Rationale: The evolution of the Antarctic ice sheet in response to climate change remains the single largest uncertainty in projecting future sea-level rise (SLR), with risk-averse projections for 2100 spanning between zero and a half meter. During the past decade, DOE has made substantial investments in new ice sheet and Earth system models needed to improve both understanding and predictive capability in this area. Yet major challenges remain, including better understanding when and under what circumstances significant increases in SLR from Antarctica may be initiated and better quantifying uncertainties in model-based SLR projections. Here, we discuss the potential for transformative advances in these areas through the application of machine learning and artificial intelligence (ML and AI, respectively).

54 ENVIRONMENTAL SCIENCES↗
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.↗