Search NASA⌕ Search

SEARCH · Search NASA

Results for “Icing”

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 523 records · Page 29

Expanded Understanding of the Western Antarctic Peninsula Sea‐Ice Environment Through Local and Regional Observations at Palmer Station

Abstract The Western Antarctic Peninsula (WAP) has been experiencing rapid regional warming since at least the 1950s, however, the impacts of this warming at the local scale are variable and nuanced. Previous studies that have linked sea‐ice variability to biogeochemical cycles and food web dynamics often combine local‐scale biogeochemical data with coarse‐resolution regional satellite sea‐ice data, which may not adequately capture local sea‐ice conditions. In this study, we analyzed local‐scale in situ sea‐ice observations collected as part of a 28‐year record (1992–2020) from the Palmer Long‐Term Ecological Research site at Anvers Island, mid‐WAP, in conjunction with isotopically‐derived sea‐ice meltwater (SIM) fractions and satellite‐derived sea‐ice motion and concentration, to quantify the variability and long‐term trends in local sea‐ice behavior. In situ sea ice observations at Palmer Station displayed higher variability than satellite observations and showed no significant declines over this time, despite region‐wide declines identified in prior studies. Higher spring SIM fractions were attributed to strong northward sea‐ice motion throughout the winter. Applying these local‐scale sea‐ice insights to similarly scaled stratification and chlorophyll‐ a measurements, we found that a longer‐lasting, more consistent sea‐ice pack led to greater water column stratification following the spring sea‐ice retreat. Greater sea‐ice persistence and stronger stratification led to larger peaks in chlorophyll‐ a , though sea‐ice metrics did not explain the positive temporal trends in either stratification strength or chlorophyll‐ a . Through this study, we identify how local sea‐ice observations and meltwater data can enhance satellite data to build an understanding of the intricate connections between ice, water column dynamics, and phytoplankton.

Goodell, E.↗

An integrated approach to the remote sensing of floating ice

Review article on remote sensing applications to glaciology. Ice parameters sensed include: ice cover vs open water, ice thickness, distribution and morphology of ice formations, vertical resolution of ice thickness, ice salinity (percolation and drainage of brine; flushing of ice body with fresh water), first-year ice and multiyear ice, ice growth rate and surface heat flux, divergence of ice packs, snow cover masking ice, behavior of ice shelves, icebergs, lake ice and river ice; time changes. Sensing techniques discussed include: satellite photographic surveys, thermal IR, passive and active microwave studies, microwave radiometry, microwave scatterometry, side-looking radar, and synthetic aperture radar. Remote sensing of large aquatic mammals and operational ice forecasting are also discussed.

Campbell, W. J.↗

Observation of pressure ridges in SAR images of sea ice: Scattering theory and comparison with observations

Ridges and keels (hummocks and bummocks) in sea ice flows are important in sea ice research for both scientific and practical reasons. Sea ice movement and deformation is driven by internal and external stresses on the ice. Ridges and keels play important roles in both cases because they determine the external wind and current stresses via drag coefficients. For example, the drag coefficient over sea ice can vary by a factor of several depending on the fluid mechanical roughness length of the surface. This roughness length is thought to be strongly dependent on the ridge structures present. Thus, variations in ridge and keel structure can cause gradients in external stresses which must be balanced by internal stresses and possibly fracture of the ice. Ridging in sea ice is also a sign of fracture. In a practical sense, large ridges form the biggest impediment to surface travel over the ice or penetration through sea ice by ice-strengthened ships. Ridges also play an important role in the damage caused by sea ice to off-shore structures. Hence, observation and measurement of sea ice ridges is an important component of sea ice remote sensing. The research reported here builds on previous work, estimating the characteristics of ridges and leads in sea ice from SAR images. Our objective is to develop methods for quantitative measurement of sea ice ridges from SAR images. To make further progress, in particular, to estimate ridge height, a scattering model for ridges is needed. Our research approach for a ridge scattering model begins with a survey of the geometrical properties of ridges and a comparison with the characteristics of the surrounding ice. For this purpose we have used airborne optical laser (AOL) data collected during the 1987 Greenland Sea Experiment. These data were used to generate a spatial wavenumber spectrum for height variance for a typical ridge - the typical ridge is the average over 10 large ridges. Our first-order model radar scattering includes both the quasi-specular and Bragg resonant scatter mechanisms. This model is extended to include contributions from volume scatter and scatter from discrete objects. Geometrical characteristics from the AOL survey and model calculations imply that for radar wavelengths and observation geometries that are dominated by the quasi-specular scattering mechanism radar backscatter from a ridge is a measure of peak ridge height. We present scattering model results and compare them with ridges observed during the LEADEX experiment of March-April 1992 when both X, C, and L-band aircraft SAR and the ERS-1 satellite SAR observed a region in the Beaufort Sea near 86 deg N, 10 deg W. Data were also collected documenting ridge characteristics on the surface. The surface data are used to generate a SAR signature via the scattering model described above. The predicted SAR signatures compare well with the SAR observations.

Vesecky, J. F.↗

Light transmission and reflection in perennially ice-covered Lake Hoare, Antarctica

We have investigated the transmission and albedo of the perennial ice cover on Lake Hoare, Antarctica. Our database consists of year-round measurements of the photosynthetically active radiation (400-700 nm) under the ice, measurements of the spatial variation of the under-ice light in midsummer, and spectrally resolved measurements from 400 to 700 nm of the albedo and transmission of the ice cover in early (November) and in midsummer (January). Our results show that the transmission decreases in the first part of summer, dropping by a factor of approximately 4 from November to January. We suggest that this is due to heating in the upper layers of the ice cover and the formation of Tyndall figures. Later in the summer when a significant liquid water fraction occurs within the ice cover, the transmission increases. In the fall when the ice cover freezes solid the transmission drops markedly. The spectrally resolved measurements from 400 to 700 nm show that approximately 2-5% of the incident light in this spectral region penetrates the 3.5-m thick ice cover. We have analyzed the spectral data using a two-stream scattering solution to the radiative transfer equation with three vertical layers in the ice cover. A surficial glaze of scattering ice 1 cm thick overlies a layer of sandy, bubbly ice about a meter thick, and below this is a thick layer of sand-free ice with bubbles. We find that the ice cover is virtually opaque at wavelengths longer than 800 nm. The net transmission of solar energy is approximately 2%. Significant changes in the thickness of the ice cover have been reported at Lake Hoare. These are due primarily to changes in the thickness of the bottom layer only. Because this layer is relatively clear, the effect on the transmission through the ice cover from these changes is less than would be predicted assuming a homogeneous ice cover.

Mckay, C. P.↗

Spatial patterns in the length of the sea ice season in the Southern Ocean, 1979-1986

The length of the sea ice season summarizes in one number the ice coverage conditions for an individual location for an entire year. It becomes a particularly valuable variable when mapped spatially over a large area and examined for regional and interannual differences, as is done here for the Southern Ocean over the years 1979-1986, using the satellite passive microwave data of the Nimbus 7 scanning multichannel microwave radiometer. Three prominent geographic anomalies in ice season lengths occur consistently in each year of the data set, countering the general tendency toward shorter ice seasons from south to north: (1) in the Weddell Sea the tendency is toward shorter ice seasons from southwest to northeast, reflective of the cyclonic ice/atmosphere/ocean circulations in the Weddell Sea region. (2) Directly north of the Ross Ice Shelf anomalously short ice seasons occur, lasting only 245-270 days, in contrast to the perennial ice coverage at comparable latitudes in the southern Bellingshausen and Amundsen Seas and in the western Weddell Sea. The short ice season off the Ross Ice Shelf reflects the consistently early opening of the ice cover each spring, under the influence of upwelling along the continental slope and shelf and atmospheric forcing from winds blowing off the Antarctic continent. (3) In the southern Amundsen Sea, anomalously short ice seasons occur adjacent to the coast, owing to the frequent existence of coastal polynyas off the many small ice shelves bordering the sea. Least squares trends in the ice season lengths over the 1979-1986 period are highly coherent spatially, with overall trends toward shorter ice seasons in the northern Weddell and Bellingshausen seas and toward longer ice seasons in the Ross Sea, around much of East Antarctica, and in a portion of the south central Weddell Sea.

Parkinson, Claire L.↗

Modern Airfoil Ice Accretions

This report presents results from the first icing tests performed in the Modem Airfoils program. Two airfoils have been subjected to icing tests in the NASA Lewis Icing Research Tunnel (IRT). Both airfoils were two dimensional airfoils; one was representative of a commercial transport airfoil while the other was representative of a business jet airfoil. The icing test conditions were selected from the FAR Appendix C envelopes. Effects on aerodynamic performance are presented including the effects of varying amounts of glaze ice as well as the effects of approximately the same amounts of glaze, mixed, and rime ice. Actual ice shapes obtained in these tests are also presented for these cases. In addition, comparisons are shown between ice shapes from the tests and ice shapes predicted by the computer code, LEWICE for similar conditions. Significant results from the tests are that relatively small amounts of ice can have nearly as much effect on airfoil lift coefficient as much greater amounts of ice and that glaze ice usually has a more detrimental effect than either rime or mixed ice. LEWICE predictions of ice shapes, in general, compared reasonably well with ice shapes obtained in the IRT, although differences in details of the ice shapes were observed.

Addy, Harold E., Jr.↗

A Numerical Evaluation of Icing Effects on a Natural Laminar Flow Airfoil

As a part of CFD code validation efforts within the Icing Branch of NASA Glenn Research Center, computations were performed for natural laminar flow (NLF) airfoil, NLF-0414. with 6 and 22.5 minute ice accretions. Both 3-D ice castings and 2-D machine-generated ice shapes were used in wind tunnel tests to study the effects of natural ice is well as simulated ice. They were mounted in the test section of the Low Turbulence Pressure Tunnel (LTPT) at NASA Langley that the 2-dimensionality of the flow can be maintained. Aerodynamic properties predicted by computations were compared to data obtained through the experiment by the authors at the LTPT. Computations were performed only in 2-D and in the case of 3-D ice, the digitized ice shape obtained at one spanwise location was used. The comparisons were mainly concentrated on the lift characteristics over Reynolds numbers ranging from 3 to 10 million and Mach numbers ranging from 0.12 to 0.29. WIND code computations indicated that the predicted stall angles were in agreement with experiment within one or two degrees. The maximum lift values obtained by computations were in good agreement with those of the experiment for the 6 minute ice shapes and the minute 3-D ice, but were somewhat lower in the case of the 22.5 minute 2-D ice. In general, the Reynolds number variation did not cause much change in the lift values while the variation of Mach number showed more change in the lift. The Spalart-Allmaras (S-A) turbulence model was the best performing model for the airfoil with the 22.5 minute ice and the Shear Stress Turbulence (SST) turbulence model was the best for the airfoil with the 6 minute ice and also for the clean airfoil. The pressure distribution on the surface of the iced airfoil showed good agreement for the 6 minute ice. However, relatively poor agreement of the pressure distribution on the upper surface aft of the leading edge horn for the 22.5 minute ice suggests that improvements are needed in the grid or turbulence models.

Chung, James J.↗

Satellite-Derived Dynamics of Southern Ocean Sea Ice

Antarctic ERS-2, RADARSAT Synthetic Aperture Radar and ERS-1/2 Scatterometer images were analyzed with SMMI radiometer image time-series data to investigate seasonal variability in satellite-tracked sea-ice dynamics in the Southern Ocean during 1992. Supporting field data were acquired during 'in-situ' experiments including the winter 1992 Ice Station Weddell and Winter Weddell Gyre studies. A variety of surface measurements were made during these experiments including Argos-buoy deployment and GPS drift measurements. These are used in conjunction with International Program for Antarctic Buoys drift trajectories for ice-motion tracking validation. Comparisons between gridded Special Sensor Microwave Imager (SSMI) ice-motion vectors and European Center for Medium Range Weather Forecasts/National Centers for Environmental Prediction (ECMWF/NCEP) analyses indicate that large-scale drift is forced predominantly by the long-term mean, large-scale synoptic pressure field. Only sub-daily SAR sea-ice tracking can capture high-frequency fluctuations, driven by polar lows or tidal forcing. In these cases, sea-ice drift can respond rapidly to changes in forcing on semi-diurnal time scales depending on the location with respect to the coastline. Seasonality of ice drift, particularly in the Weddell and Ross Seas, is linked to ice extent and compactness, and internal ice stresses transmitted through the pack ice from the coast. Three-monthly seasonal climatologies are presented of austral winter of ice drift in the Southern Ocean. The large Weddell and Ross Sea gyres are clearly resolved along with key seasonal and spatial attributes of their cyclonic circulation. Regional time series of ice dynamics parameters are used to illustrate correlations with meteorological forcing. Persistent divergence such as that occurring in the Ronne-Filchner polynya system results in large fractions of new ice. Similarly, convergence zones produce large fractions of deformed ice and characterize the dynamics of regions where perennial ice is observed. High shear strains also help delineate the axis of the Antarctic divergence in many places in the ice cover. In these regions, the separation between the coastal 'east wind drift'and ACC-dominated drift regimes is characterized by zonally extended regions of intense shear.

Drinkwater, Mark R.↗

Sea Ice Motion from Wavelet Analysis of Satellite Data

Wavelet analysis of NASA scatterometer (NSCAT) backscatter and Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave/Imager (SSM/I) radiance data can be used to obtain daily sea ice drift information for the Arctic region. This technique provides improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing data from satellite synthetic aperture radars. Comparisons with ice motion derived from ocean buoys give good quantitative agreement. Both comparison results from NSCAT and SSM/I are compatible, and the results from NSCAT can definitely complement that from SSM/I when there are cloud or surface effects. Then three sea-ice drift daily results from NSCAT, SSM/I, and buoy data can be merged as a composite map by some data fusion techniques. The ice flow streamlines are highly correlated with surface air pressure contours. Examples of derived ice-drift maps in December 1996 illustrate large-scale circulation reversals over a period of four days. A method for deriving divergence and shear at the large-scale has been developed and comparison between buoys and satellite results shows a good agreement. These calibrated/validated results indicate that NSCAT, SSM/I merged daily ice motion are suitably accurate to identify and closely locate sea ice processes, and to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model. For demonstration purpose, the ice velocities derived from satellite data are compared with the ice velocities derived from a coupled ice-ocean interaction model. The comparison reveals that the general circulation patterns of the two are quite similar but the ice velocity differences between the two are quite significant. In order to quantify the wind effects on ice motion, empirical orthogonal functions (EOF) are used in the principal component analysis for both ice motion and pressure field. Some preliminary results of sea-ice motion from QuikScat will also be presented.

Liu, Antony K.↗

Space/Time Statistics of Polar Ice Motion

Ice motions have been computed from passive microwave imagery (SMMR and SSM/I) on a daily basis for both Polar Regions. In the Arctic these daily motions have been merged with daily motions from AVHRR imagery and the Arctic buoy program. In the Antarctic motion only from the AVHRR were available for merging with the passive microwave vectors. Long-term means, monthly means and weekly means have all been computed from the resulting 22-year time series of polar ice motion. Papers are in preparation that present the long term (22 year) means, their variability and show animations of the monthly means over this time period for both Polar Regions. These papers will have links to "enhanced objects" that allow the reader to view the animations as part of the paper. The first paper presents the ice motion results from each of the Polar Regions. The second paper looks only at ice motion in the Arctic in order to develop a time series of ice age in the Arctic. Starting with the first full SMMR year in 1979 we keep track of each individual "ice element" (resolution of the sensor) and track it in the subsequent monthly time series. After a year we "age" each "particle" and we thus can keep track of the age of the ice starting in 1979. We keep track of ice age classes between one and five years and thus we can see the evolution of the ice as it ages after the initial 5-year period. This calculation shows how we are losing the older ice through Fram Strait at a rather alarming rate particularly in the past 15 years. This loss of older ice has resulted in an overall decrease in the thickest, oldest ice, which is now limited to a region just north of the Canadian Archipelago with tongues extending out across the pole towards the Siberian Shelf. This loss of old ice is consistent with the effects of global warming which provides the heat needed to melt, move and disperse this oldest ice through Fram Strait. This is the first step in a progression that may eventually open the Arctic ice pack and lead to an ice-free Arctic Ocean.

Emery, William J.↗

Airframe Icing Research Gaps: NASA Perspective

qCurrent Airframe Icing Technology Gaps: Development of a full 3D ice accretion simulation model. Development of an improved simulation model for SLD conditions. CFD modeling of stall behavior for ice-contaminated wings/tails. Computational methods for simulation of stability and control parameters. Analysis of thermal ice protection system performance. Quantification of 3D ice shape geometric characteristics Development of accurate ground-based simulation of SLD conditions. Development of scaling methods for SLD conditions. Development of advanced diagnostic techniques for assessment of tunnel cloud conditions. Identification of critical ice shapes for aerodynamic performance degradation. Aerodynamic scaling issues associated with testing scale model ice shape geometries. Development of altitude scaling methods for thermal ice protections systems. Development of accurate parameter identification methods. Measurement of stability and control parameters for an ice-contaminated swept wing aircraft. Creation of control law modifications to prevent loss of control during icing encounters. 3D ice shape geometries. Collection efficiency data for ice shape geometries. SLD ice shape data, in-flight and ground-based, for simulation verification. Aerodynamic performance data for 3D geometries and various icing conditions. Stability and control parameter data for iced aircraft configurations. Thermal ice protection system data for simulation validation.

Potapczuk, Mark↗

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.↗

Examination of Icing Induced Loss of Control and Its Mitigations

Factors external to the aircraft are often a significant causal factor in loss of control (LOC) accidents. In today s aviation world, very few accidents stem from a single cause and typically have a number of causal factors that culminate in a LOC accident. Very often the "trigger" that initiates an accident sequence is an external environment factor. In a recent NASA statistical analysis of LOC accidents, aircraft icing was shown to be the most common external environmental LOC causal factor for scheduled operations. When investigating LOC accident or incidents aircraft icing causal factors can be categorized into groups of 1) in-flight encounter with super-cooled liquid water clouds, 2) take-off with ice contamination, or 3) in-flight encounter with high concentrations of ice crystals. As with other flight hazards, icing induced LOC accidents can be prevented through avoidance, detection, and recovery mitigations. For icing hazards, avoidance can take the form of avoiding flight into icing conditions or avoiding the hazard of icing by making the aircraft tolerant to icing conditions. Icing detection mitigations can take the form of detecting icing conditions or detecting early performance degradation caused by icing. Recovery from icing induced LOC requires flight crew or automated systems capable of accounting for reduced aircraft performance and degraded control authority during the recovery maneuvers. In this report we review the icing induced LOC accident mitigations defined in a recent LOC study and for each mitigation describe a research topic required to enable or strengthen the mitigation. Many of these research topics are already included in ongoing or planned NASA icing research activities or are being addressed by members of the icing research community. These research activities are described and the status of the ongoing or planned research to address the technology needs is discussed

Reehorst, Andrew L.↗

Independent Effects of Reynolds and Mach Numbers on the Aerodynamics of an Iced Swept Wing

Aerodynamic assessment of icing effects on swept wings is an important component of a larger effort to improve three-dimensional icing simulation capabilities. An understanding of ice-shape geometric fidelity and Reynolds and Mach number effects on the iced-wing aerodynamics is needed to guide the development and validation of ice-accretion simulation tools. To this end, wind-tunnel testing was carried out for a 13.3%-scale semispan wing based upon the Common Research Model airplane configuration. The wind-tunnel testing was conducted at the ONERA F1 pressurized wind tunnel with Reynolds numbers of 1.6×10(exp 6) to 11.9×10(exp 6) and Mach numbers of 0.09 to 0.34. Five different configurations were investigated using fully 3D, high-fidelity artificial ice shapes that maintain nearly all of the 3D ice accretion features documented in prior icing-wind tunnel tests. These large, leadingedge ice shapes were nominally based upon airplane holding in icing conditions scenarios. For three of these configurations, lower-fidelity simulations were also built and tested. The results presented in this paper show that while Reynolds and Mach number effects are important for quantifying the clean-wing performance, there is very little to no effect for an iced-wing with 3D, high-fidelity artificial ice shapes or 3D smooth ice shapes with grit roughness. These conclusions are consistent with the large volume of past research on icedairfoils. However, some differences were also noted for the associated stalling angle of the iced swept wing and for various lower-fidelity versions of the leading-edge ice accretion. More research is planned to further investigate the key features of ice accretion geometry that must be simulated in lower-fidelity versions in order to capture the essential aerodynamics.

Broeren, Andy P.↗

Antarctic Ice Shelf Thickness Change from Multimission Lidar Mapping

We calculate rates of ice thickness change and bottom melt for ice shelves in West Antarctica and the Antarctic Peninsula from a combination of elevation measurements from NASA–CECS Antarctic ice mapping campaigns and NASA Operation IceBridge corrected for oceanic processes from measurements and models, surface velocity measurements from synthetic aperture radar, and high-resolution outputs from regional climate models. The ice thickness change rates are calculated in a Lagrangian reference frame to reduce the effects from advection of sharp vertical features, such as cracks and crevasses, that can saturate Eulerian-derived estimates. We use our method over different ice shelves in Antarctica, which vary in terms of size, repeat coverage from airborne altimetry, and dominant processes governing their recent changes. We find that the Larsen-C Ice Shelf is close to steady state over our observation period with spatial variations in ice thickness largely due to the flux divergence of the shelf. Firn and surface processes are responsible for some short-term variability in ice thickness of the Larsen-C Ice Shelf over the time period. The Wilkins Ice Shelf is sensitive to short-timescale coastal and upper-ocean processes, and basal melt is the dominant contributor to the ice thickness change over the period. At the Pine Island Ice Shelf in the critical region near the grounding zone, we find that ice shelf thickness change rates exceed 40m/yr, with the change dominated by strong submarine melting. Regions near the grounding zones of the Dotson and Crosson ice shelves are decreasing in thickness at rates greater than 40m/yr, also due to intense basal melt. NASA–CECS Antarctic ice mapping and NASA Operation IceBridge campaigns provide validation datasets for floating ice shelves at moderately high resolution when coregistered using Lagrangian methods.

Sutterley, Tyler C.↗

Analyzing the Ages of South Polar Craters on the Moon: Implications for the Sources and Evolution of Surface Water Ice

Surface water ice in the permanently shadowed polar regions of the Moon has a patchy surficial distribution and is not found within all available cold-trapping areas. To date it is not well understood when the ice was delivered, which has important implications for the surficial characteristics of the ice as well as for possible delivery mechanisms. Here we present absolute model ages for 20 south polar craters that host surface water ice, providing maximum estimates of the ages of surface ice contained within these craters. We quantify the amount of available cold-trapping surface area that is occupied by water ice in order to examine the relationship between the patchiness of ice within each crater and the age of each host crater. The majority of surface ice is contained in old craters ≥∼3.1 Gyr, where the majority of cold-trapping area on the pole exists. The ice is these ancient craters is very patchy in surficial distribution, occupying <11.5% of cold-trapping surface area available in individual craters. This patchy distribution of ice in old craters is likely to be due to impact bombardment and regolith overturn within the polar regions. Interestingly, surface ice is also located within smaller craters (<15 km in diameter), whose sharp crater rim crest morphologies suggest that they may be relatively young. Ice in fresh-looking craters suggests that ice has been delivered to the lunar surface more recently, perhaps from micrometeorites or through solar wind interactions with the lunar regolith. Finally, we also analyze a group of ancient craters that does not host surface water ice, even though these craters are present-day cold traps. These specific ancient craters would not have been thermally stable for the cold-trapping of water ice before the onset of true polar wander suggested by Siegler et al. (2016). If true polar wander did occur on the Moon, then the ages of ice-bearing craters presented here set an upper limit for the age of post-true polar wander hydrogen emplacement of 4.1 ± 0.1 Gyr.

Deutsch, Ariel N.↗

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

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

Patrick Taylor↗

Determination of Ice Cloud Models Using MODIS and MISR Data

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

MODIS↗