Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “snow contribution”

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 181 records · Page 10

Energetic Charged Particle Component or the NO(y) Budget of the Polar Middle Atmosphere

Analysis of nitrates measured in polar ice cap snow at a high resolution shows large variations in the nitrates. It has been shown that the nitrate signal may contain a signature of solar activity [Zeller and Dreschhoff, 19951. Reactive odd nitrogen production associated with solar particle events (SPEs) and auroral activity may be a source of some of the nitrate anomalies observed in the polar ice caps. Periods of large SPEs can lead to a production of polar atmospheric odd nitrogen in excess of the ambient sources in the polar stratosphere and mesosphere, and may leave a large nitrate signal stratified in the polar ice cap. Auroral electrons and photoelectrons produce odd nitrogen in the thermosphere, some of which may be transported to the polar (>50 degrees) mesosphere and stratosphere. Sources of odd nitrogen in the polar middle atmosphere associated with SPEs, galactic cosmic rays, and auroral electron precipitation have been quantified. The relative contributions by the energetic particles sources to the Noy budget of the polar middle atmosphere (from tropopause to 50 km, from 50 degrees to 90 degrees latitude) are compared with the nitrates observed in the polar ice sheets.

Vitt, F. M.↗

A Study of the Kinematic and Dynamic Processes associated with Mesoscale Snowbands in a Midwest United States Snowstorm on 5 February 2020

Mesoscale snowbands are frequently responsible for dumping large amounts of snow over a relatively small region and pose a considerable forecasting challenge. The present study focuses on the mesoscale snowbands observed on 5 February 2020 over the Midwest United States. In order to improve the forecasting of mesoscale snowbands, it is important to understand the processes that organize these features. To this end, we are investigating the kinematic and dynamic processes that shape the mesoscale snowbands using a combination of high-altitude airborne radar, in-situ aircraft measurements, surface observations, and model simulations. This presentation will primarily focus on the three-dimensional structure of the mesoscale snowbands in terms of EXRAD reflectivity and horizontal winds. The horizontal winds are retrieved from the EXRAD scanning-beam Doppler velocities using a VAD technique. These winds, in combination with virtual potential temperature taken from the HRRR operational analysis, are used to compute the bulk Richardson number in order to investigate any contributions to banding features that might be tied to Kelvin-Helmholtz instability. Our analysis does not support the presence of Kelvin-Helmholtz instability. That said, low, but not subcritical, Richardson numbers were found near the top of some banding features in this event and it is possible that the model analysis is smoothing out the virtual potential temperature fields. Finally, we present arguments for using a flight module designed to sample the along-band structures with particular emphasis on observing generating cells and low-level microphysics.

Charles N. Helms↗

The anomalous extinction curve in the direction of Rho Ophiuchi from 950 to 1180 A

We report the interstellar extinction from 950 to 1180 A in the direction of Rho Oph. These data were obtained with a moderate-resolution sounding rocket observation of Ro Oph, and two lightly reddened B stars of identical spectral type. Unlike Voyager observations in this wavelength region, our spectral resolution is sufficient to directly measure the contribution of molecular hydrogen along the line of sight and unambiguously identify the dust contribution. The derived curves compare favorably with the extrapolation of functional fits to the Rho Oph derived from longer wavelengths. In addition, we have fitted a model of dust extinction to the derived curve and find that the best fits indicate that the line of sight exhibits a larger than average grain size and a preferential concentration of graphite relative to silicates. The sharp rise seen at the shortest wavelengths (below 1000 A) cannot be explained by standard silicate-graphite dust models and may be indicative of a nonstandard grain shape or a missing component in the extinction models.

Green, James C.↗

Explicitly Resolving Lightning and Electrification Processes from the 10-12 April 2019 Thundersnow Outbreak

The 10-12 April 2019 thundersnow (i.e., lightning within snowfall) outbreak was examined via ground- and space-based lightning observations and was simulated using a numerical weather prediction model with an explicit electrification parameterization. When compared to observations, the simulation propagated the synoptic snowband two to six hours faster while also exaggerating the 3-D reflectivity structure. Throughout the event, the simulation produced 1,733 thundersnow flashes which was less than what was observed by ground- and space-based lightning sensors. In general, simulated thundersnow flashes were spatially offset from the largest reflectivities within the synoptic snowband and tended to occur within elevated convection that traversed isentropically along the top of mid-level frontogenesis. These simulated thundersnow flashes were associated with a tripole charge structure with ice/snow hydrometeors contributing most to the main negative charge region. Both simulated and observed thundersnow flashes initiated in conditionally unstable environments. Lastly, a conceptual model was developed to explain the spatial separation between the largest reflectivities in the snowband and the occurrence of thundersnow. It is hypothesized that the spatial offset of thundersnow initiation from the reflectivity cores within the synoptic snowband arose from a thermal circulation – induced by mid-level frontogenesis – that advects positively charged ice/snow hydrometeors towards the surface and creates a nearly homogeneous vertical charge structure.

lightning↗

Musculoskeletal adaptations to weightlessness and development of effective countermeasures

A Research Roundtable, organized by the American College of Sports Medicine with sponsorship from the National Aeronautics and Space Administration, met in November 1995 to define research strategies for effective exercise countermeasures to weightlessness. Exercise was considered both independently of, and in conjunction with, other therapeutic modalities (e.g., pharmacological nutritional, hormonal, and growth-related factors) that could prevent or minimize the structural and functional deficits involving skeletal muscle and bone in response to chronic exposure to weightlessness, as well as return to Earth baseline function if a degree of loss is inevitable. Musculoskeletal deficits and countermeasures are described with respect to: 1) muscle and connective tissue atrophy and localized bone loss, 2) reductions in motor performance, 3) potential proneness to injury of hard and soft tissues, and 4) probable interaction between muscle atrophy and cardiovascular alterations that contribute to the postural hypotension observed immediately upon return from space flight. In spite of a variety of countermeasure protocols utilized previously involving largely endurance types of exercise, there is presently no activity-specific countermeasure(s) that adequately prevent or reduce musculoskeletal deficiencies. It seems apparent that countermeasure exercises that have a greater resistance element, as compared to endurance activities, may prove beneficial to the musculoskeletal system. Many questions remain for scientific investigation to identify efficacious countermeasure protocols, which will be imperative with the emerging era of long-term space flight.

Non-NASA Center↗

Analysis of Multiple Precipitation Products and Preliminary Assessment of Their Impact on Global Land Data Assimilation System (GLDAS) Land Surface States

Land surface models (LSMs) are computer programs, similar to weather and climate prediction models, which simulate the stocks and fluxes of water (including soil moisture, snow, evaporation, and runoff) and energy (including the temperature of and sensible heat released from the soil) after they arrive on the land surface as precipitation and sunlight. It is not currently possible to measure all of the variables of interest everywhere on Earth with sufficient accuracy and space-time resolution. Hence LSMs have been developed to integrate the available observations with our understanding of the physical processes involved, using powerful computers, in order to map these stocks and fluxes as they change in time. The maps are used to improve weather forecasts, support water resources and agricultural applications, and study the Earth's water cycle and climate variability. NASA's Global Land Data Assimilation System (GLDAS) project facilitates testing of several different LSMs with a variety of input datasets (e.g., precipitation, plant type). Precipitation is arguably the most important input to LSMs. Many precipitation datasets have been produced using satellite and rain gauge observations and weather forecast models. In this study, seven different global precipitation datasets were evaluated over the United States, where dense rain gauge networks contribute to reliable precipitation maps. We then used the seven datasets as inputs to GLDAS simulations, so that we could diagnose their impacts on output stocks and fluxes of water. In terms of totals, the Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP) had the closest agreement with the US rain gauge dataset for all seasons except winter. The CMAP precipitation was also the most closely correlated in time with the rain gauge data during spring, fall, and winter, while the satellitebased estimates performed best in summer. The GLDAS simulations revealed that modeled soil moisture is highly sensitive to precipitation, with differences in spring and summer as large as 45% depending on the choice of precipitation input.

Gottschalck, Jon↗

Monitoring Atmospheric CO2 From Space: Challenge & Approach

Atmospheric CO2 is the key radiative forcing for the Earth's climate and may contribute a major part of the Earth's warming during the past 150 years. Advanced knowledge on the CO2 distributions and changes can lead considerable model improvements in predictions of the Earth's future climate. Large uncertainties in the predictions have been found for decades owing to limited CO2 observations. To obtain precise measurements of atmospheric CO2, certain challenges have to be overcome. For an example, global annual means of the CO2 are rather stable, but, have a very small increasing trend that is significant for multi-decadal long-term climate. At short time scales (a second to a few hours), regional and subcontinental gradients in the CO2 concentration are very small and only in an order of a few parts per million (ppm) compared to the mean atmospheric CO2 concentration of about 400 ppm, which requires atmospheric CO2 space monitoring systems with extremely high accuracy and precision (about 0.5 ppm or 0.125%) in spatiotemporal scales around 75 km and 10-s. It also requires a decadal-scale system stability. Furthermore, rapid changes in high latitude environments such as melting ice, snow and frozen soil, persistent thin cirrus clouds in Amazon and other tropical areas, and harsh weather conditions over Southern Ocean all increase difficulties in satellite atmospheric CO2 observations. Space lidar approaches using Integrated Path Differential Absorption (IPDA) technique are considered to be capable of obtaining precise CO2 measurements and, thus, have been proposed by various studies including the 2007 Decadal Survey (DS) of the U.S. National Research Council. This study considers to use the Intensity-Modulated Continuous-Wave (IM-CW) lidar to monitor global atmospheric CO2 distribution and variability from space. Development and demonstration of space lidar for atmospheric CO2 measurements have been made through joint adventure of NASA Langley Research Center and Exelis, Inc. As prototype space IPDA lidars, airborne laser absorption lidar systems operating in 1.57 CO2 absorption band have been developed and tested through lab, ground-based range, and flight campaigns. Very encouraging results have been obtained. The signal-to-noise ratio (SNR) for clear sky IPDA measurements of CO2 differential absorption optical depth (DAOD) for a 10-s integration over vegetated areas with about 10 km range was found to be as high as 1300, resulting in an error 0.077% or equivalent CO2 mixing ratio (XCO2) column precision of ~0.3 ppm. Precise range measurements using the IM-CW lidar approach were also achieved, and the uncertainties have been shown to be at sub meter level. Based on the airborne lidar development, space lidar and atmospheric CO2 observations are simulated. It shows that with the IM-CW approach, accurate atmospheric CO2 measurements can be achieved from space, and a space mission such as that proposed by the DS will meet science goals of atmospheric CO2 monitoring.

Lin, Bing↗

Polarimetric backscattering from sea ice

Polarimetric backscattering from sea ice is presented. Theoretical models for backscattering are first presented for various ice types. Then, theoretical results are compared with experimental data for new thin ice, first-year ice, and multi-year ice. Sea ice is modeled as a layer medium containing random scatterers and rough interfaces. For multi-year sea ice with snow cover, the sea ice layer is modeled as an ice background with embedded spheroidal air bubbles and the snow layer as air with ice grains. The hummocky topography on multi-year ice is characterized by a Gaussian distribution which has an averaging effect on backscattering coefficients. First-year sea ice is described by an ice medium hosting ellipsoidal brine inclusions. These inclusions are oriented preferentially in the vertical direction due to the columnar structure of first-year sea ice. Azimuthally, the orientation of the brine pockets are random corresponding to the random c-axes in sea ice, unless the axes are oriented by sea currents. For thin ice in newly opened leads, it has been observed that there exists a thin brine layer with very high salinity on the top surface of the new ice. This brine layer is depicted as a medium with high permittivity which can significantly affect electromagnetic scattering signatures from the lower thin ice layer, with a higher fractional volume of brine inclusions due to higher salinity as compared to thick first-year sea ice. The rough medium interfaces are described as Gaussian rough surfaces characterized by root-mean-square heights and surface correlation lengths. The contribution from rough surface, calculated under the Kirchhoff approximation or small perturbation method, is assumed to be independent from volume scattering. The total loss including absorption and scattering losses in the scattering media is represented by the imaginary part of effective permittivity obtained from the strong fluctuation theory. The polarimetric scattering coefficients for different ice types are then derived under the distorted Born approximation.

Nghiem, S. V.↗

Elevation change (2000-2004) on the Malaspina Glacier, Alaska

The glaciers of the southeastern Alaska coastal region are the largest temperate glacier meltwater source on Earth and may contribute one third of the total glacier meltwater entering the global ocean. Since melt onset and refreeeze timing in this region show a tendency toward earlier onset and longer ablation seasons, accelerated glacier wastage may be occurring. In this study we focus on one of the largest temperate glacier systems on Earth, the Malaspina Glacier. This glacier, with a length of approximately 110 km and an area of approximately square 5,000 km, has the largest piedmont lobe of any temperate glacier. The entire lobe, which lies at elevations below 600 m, is within the ablation zone. We report and interpret ice elevation change between a digital elevation model (DEM) derived from the Shuttle Radar Topography Mission (SRTM C band) observations in Feb. 2000 and ICESat Laser 1-3 observations between Feb. 2003 and Nov. 2004. We use these elevation change results, along with earlier studies, to address the spatial and temporal variability in wastage of the piedmont lobe. Between 2000 and 2004 ice elevation changes of 10-30 meters occurred across the central Malaspina piedmont lobe. From 1972/73 (USGS DEM) to 1999 (SRTM corrected for estimated winter snow accumulation) Malaspina's (Agassiz, Seward Lobe, and Marvine) mean ice thinning was estimated at -47 m with maximum thinning on parts of the lobes to -160 m. The Malaspina's accumulation area is only slightly larger than its ablation area (2,575 km2 vs. 2,433 km2); unfortunately few glaciological observations are available from this source region. Snow accumulation rates have been largely inferred from low-altitude precipitation and temperature data. Comparing sequential ICESat observations in the Malaspina source region, we estimated short-term elevation increases of up to 5 meters during the winter of 2003/04.

Sauber, J.↗

Evaluation of the Surface Representation of the Greenland Ice Sheet in a General Circulation Model

Simulated surface conditions of the Goddard Earth Observing System model, version 5 (GEOS 5) atmospheric general circulation model (AGCM) are examined for the contemporary Greenland Ice Sheet (GrIS). A surface parameterization that explicitly models surface processes including snow compaction, meltwater percolation and refreezing, and surface albedo is found to remedy an erroneous deficit in the annual net surface energy flux and provide an adequate representation of surface mass balance (SMB) in an evaluation using simulations at two spatial resolutions. The simulated 1980-2008 GrIS SMB average is 24.7+/-4.5 cm yr(- 1) water-equivalent (w.e.) at.5 degree model grid spacing, and 18.2+/-3.3 cm yr(- 1) w.e. for 2 degree grid spacing. The spatial variability and seasonal cycle of the simulation compare favorably to recent studies using regional climate models, while results from 2 degree integrations reproduce the primary features of the SMB field. In comparison to historical glaciological observations, the coarser resolution model overestimates accumulation in the southern areas of the GrIS, while the overall SMB is underestimated. These changes relate to the sensitivity of accumulation and melt to the resolution of topography. The GEOS-5 SMB fields contrast with available corresponding atmospheric models simulations from the Coupled Model Intercomparison Project (CMIP5). It is found that only a few of the CMIP5 AGCMs examined provide significant summertime runoff, a dominant feature of the GrIS seasonal cycle. This is a condition that will need to be remedied if potential contributions to future eustatic change from polar ice sheets are to be examined with GCMs.

AGCM↗

Greenland Ice Sheet Surface Temperature, Melt, and Mass Loss: 2000-2006

Extensive melt on the Greenland Ice Sheet has been documented by a variety of ground and satellite measurements in recent years. If the well-documented warming continues in the Arctic, melting of the Greenland Ice Sheet will likely accelerate, contributing to sea-level rise. Modeling studies indicate that an annual or summer temperature rise of 1 C on the ice sheet will increase melt by 20-50% therefore, surface temperature is one of the most important ice-sheet parameters to study for analysis of changes in the mass balance of the ice-sheet. The Greenland Ice Sheet contains enough water to produce a rise in eustatic sea level of up to 7.0 m if the ice were to melt completely. However, even small changes (centimeters) in sea level would cause important economic and societal consequences in the world's major coastal cities thus it is extremely important to monitor changes in the ice-sheet surface temperature and to ultimately quantify these changes in terms of amount of sea-level rise. We have compiled a high-resolution, daily time series of surface temperature of the Greenland Ice Sheet, using the I-km resolution, clear-sky land-surface temperature (LST) standard product from the Moderate-Resolution Imaging Spectroradiometer (MODIS), from 2000 - 2006. We also use Gravity Recovery and Climate Experiment (GRACE) data, averaged over 10-day periods, to measure change in mass of the ice sheet as it melt and snow accumulates. Surface temperature can be used to determine frequency of surface melt, timing of the start and the end of the melt season, and duration of melt. In conjunction with GRACE data, it can also be used to analyze timing of ice-sheet mass loss and gain.

Hall, Dorothy K.↗

Design, Evaluation and GCM-Performance of a New Parameterization for Microphysics of Clouds with Relaxed Arakawa-Schubert Scheme (McRas)

A prognostic cloud scheme named McRAS (Microphysics of clouds with Relaxed Arakawa-Schubert Scheme) was developed with the aim of improving cloud-microphysics, and cloud-radiation interactions in GCMs. McRAS distinguishes convective, stratiform, and boundary-layer clouds. The convective clouds merge into stratiform clouds on an hourly time-scale, while the boundary-layer clouds do so instantly. The cloud condensate transforms into precipitation following the auto-conversion relations of Sundqvist that contain a parametric adaptation for the Bergeron-Findeisen process of ice crystal growth and collection of cloud condensate by precipitation. All clouds convect, advect, and diffuse both horizontally and vertically with a fully active cloud-microphysics throughout its life-cycle, while the optical properties of clouds are derived from the statistical distribution of hydrometeors and idealized cloud geometry. An evaluation of McRAS in a single column model (SCM) with the GATE Phase III data has shown that McRAS can simulate the observed temperature, humidity, and precipitation without discernible systematic errors. An evaluation with the ARM-CART SCM data in a cloud model intercomparison exercise shows reasonable but not an outstanding accurate simulation. Such a discrepancy is common to almost all models and is related, in part, to the input data quality. McRAS was implemented in the GEOS II GCM. A 50 month integration that was initialized with the ECMWF analysis of observations for January 1, 1987 and forced with the observed sea-surface temperatures and sea-ice distribution and vegetation properties (biomes, and soils), with prognostic soil moisture, snow-cover, and hydrology showed a very realistic simulation of cloud process, incloud water and ice, and cloud-radiative forcing (CRF). The simulated ITCZ showed a realistic time-mean structure and seasonal cycle, while the simulated CRF showed sensitivity to vertical distribution of cloud water which can be easily altered by the choice of time constant and incloud critical cloud water amount regulators for auto-conversion. The CRF and its feedbacks also have a profound effect on the ITCZ. Even though somewhat weaker than observed, the McRAS-GCM simulation produces robust 30-60 day oscillations in the 200 hPa velocity potential. Two ensembles of 4-summer (July, August, September) simulations, one each for 1987 and 1988 show that the McRAS-GCM simulates realistic and statistically significant precipitation differences over India, Central America, and tropical Africa. Several seasonal simulations were performed with McRAS-GEOS II GCM for the summer (June-July- August) and winter (December-January-February) periods to determine how the simulated clouds and CRFs would be affected by: i) advection of clouds; ii) cloud top entrainment instability, iii) cloud water inhomogeneity correction, and (iv) cloud production and dissipation in different cloud-processes. The results show that each of these processes contributes to the simulated cloud-fraction and CRF.

Sud, Y. C.↗

ICESat-2 Mission: Status, Outlook, and Contributions to Polar Science

The Ice, Cloud and land Elevation Satellite-2 (ICESat-2) mission carries the Advanced Topographic Laser Altimeter System (ATLAS) lidar to measure the changing height of Earth surface. After more than 3 years of science data collection, ATLAS has emitted well over a trillion laser pulses and continues to operate nominally. The ICESat-2 data products were initially released in May 2019, and have been used in over 150 peer-reviewed publications to date. Recent community white papers have noted the importance of ice elevation measurements for the coming decades, and as CryoSat-2 has entered it’s second decade and future missions are uncertain, ICESat-2 will be a critical part of the Earth System Observatory for the 2020s. In the polar regions, ICESat-2 has enabled year-round sea ice freeboard measurements as well as ice sheet elevation changes on seasonal, annual, and decadal time scales when combined with other missions. These data enable measurement of change at fine spatial and temporal scales to help understand or characterize the processes driving these observations. As of this presentation, Release 005 is the current version of the along-track data products, including surface-specific products for sea ice, land ice, ground and canopy height, ocean, and inland water heights. In addition, first versions of gridded data products for land ice, sea ice, and ocean are all available. Data spans the start of the mission (14 October 2018) through early 2022. The ICESat-2 project has also begun to produce and distribute QuickLook products with a nominal 3-day latency (as opposed to the ~45-day latency of the final data products). Initial versions of these QuickLook products have elevation accuracy of ~3 meters, and geolocation accuracy of ~50 meters. All products are available through the National Snow and Ice Data Center (nsidc.org). This presentation will summarize the current state and health of the observatory, current state of the ICESat-2 data products, data product quality assessment, and outlook for the future.

ICESat-2↗

Contribution of National near Real Time MODIS Forest Maximum Percentage NDVI Change Products to the U.S. ForWarn System

This presentation reviews the development, integration, and testing of Near Real Time (NRT) MODIS forest % maximum NDVI change products resident to the USDA Forest Service (USFS) ForWarn System. ForWarn is an Early Warning System (EWS) tool for detection and tracking of regionally evident forest change, which includes the U.S. Forest Change Assessment Viewer (FCAV) (a publically available on-line geospatial data viewer for visualizing and assessing the context of this apparent forest change). NASA Stennis Space Center (SSC) is working collaboratively with the USFS, ORNL, and USGS to contribute MODIS forest change products to ForWarn. These change products compare current NDVI derived from expedited eMODIS data, to historical NDVI products derived from MODIS MOD13 data. A new suite of forest change products are computed every 8 days and posted to the ForWarn system; this includes three different forest change products computed using three different historical baselines: 1) previous year; 2) previous three years; and 3) all previous years in the MODIS record going back to 2000. The change product inputs are maximum value NDVI that are composited across a 24 day interval and refreshed every 8 days so that resulting images for the conterminous U.S. are predominantly cloud-free yet still retain temporally relevant fresh information on changes in forest canopy greenness. These forest change products are computed at the native nominal resolution of the input reflectance bands at 231.66 meters, which equates to approx 5.4 hectares or 13.3 acres per pixel. The Time Series Product Tool, a MATLAB-based software package developed at NASA SSC, is used to temporally process, fuse, reduce noise, interpolate data voids, and re-aggregate the historical NDVI into 24 day composites, and then custom MATLAB scripts are used to temporally process the eMODIS NDVIs so that they are in synch with the historical NDVI products. Prior to posting, an in-house snow mask classification product is computed for the current compositing period and integrated into the change images to account for snow related NDVI drops. The supplemental snow classification product was needed because other available QA cloud/snow mask typically underestimates snow cover. MODIS true and false color composites were also computed from eMODIS reflectance data and the true color RGBs are also posted on ForWarn?s FCAV; this data is used for assessing apparent occasional quality issues on the change products due to residual unmasked cloud cover. New forest change products are posted with typical latencies of 1-2 days after the last input eMODIS data collection date for a given 24 day compositing period.

Spruce, Joseph P.↗

Photometric Properties of Quiet Sun Observed With the DKIST and HINODE

The variability of irradiance is primarily influenced by surface magnetism, such as sunspots, faculae, and network structures. Each of these features contributes differently to the overall variability. However, the impact of small magnetic elements, which are commonly observed in high resolution images of the solar photosphere but remain unresolved in full-disk images used for irradiance modeling, is still a topic of debate. Understanding the influence of these small magnetic elements to the overall brightness, particularly in quiet regions, is fundamental to understanding irradiance fluctuations over decadal and longer temporal scales. We investigate the brightness of small-size magnetic elements using very high spatial resolution spectro-polarimetric observations acquired at different spectral ranges with the National Science Foundations Daniel K. Inouye Solar Telescope (DKIST) and with HINODE/SOT. Our results are compared to 3D MHD Small Solar Dynamo simulations obtained with the MURaM code.

SMD↗

Intercomparison of Precipitation Estimates over the Arctic Ocean and Its Peripheral Seas from Reanalyses

Precipitation over the Arctic Ocean has a significant impact on the basin-scale freshwater and energy budgets but is one of the most poorly constrained variables in atmospheric reanalyses. Precipitation controls the snow cover on sea ice, which impedes the exchange of energy between the ocean and atmosphere, inhibiting sea ice growth. Thus, accurate precipitation amounts are needed to inform sea ice modeling, especially for the production of thickness estimates from satellite altimetry freeboard data. However, obtaining a quantitative estimate of the precipitation distribution in the Arctic is notoriously difficult because of a number of factors, including a lack of reliable, long-term in situ observations; difficulties in remote sensing over sea ice; and model biases in temperature and moisture fields and associated uncertainty of modeled cloud microphysical processes in the polar regions. Here, we compare precipitation estimates over the Arctic Ocean from eight widely used atmospheric reanalyses over the period 2000–16 (nominally the ‘‘new Arctic’’). We find that the magnitude, frequency, and phase of precipitation vary drastically, although interannual variability is similar. Reanalysis-derived precipitation does not increase with time as expected; however, an increasing trend of higher fractions of liquid precipitation (rainfall) is found. When compared with drifting ice mass balance buoys, three reanalyses (ERA-Interim, MERRA, and NCEP R2) produce realistic magnitudes and temporal agreement with observed precipitation events, while two products [MERRA, version 2 (MERRA-2), and CFSR] show large, implausible magnitudes in precipitation events. All the reanalyses tend to produce overly frequent Arctic precipitation. Future work needs to be undertaken to determine the specific factors in reanalyses that contribute to these discrepancies in the new Arctic.

Mcpartland, Linette N.↗

Radiometric Calibration of A Non-Imaging Airborne Spectrometer To Measure the Greenland Ice Sheet Surface

Methods to radiometrically calibrate a non-imaging airborne visible-to-shortwave infrared (VSWIR) spectrometer to measure the Greenland ice sheet surface are presented. Airborne VSWIR measurement performance for bright Greenland ice and dark bare rock/soil targets is compared against the MODerate resolution atmospheric TRANsmission (MODTRAN ® ) radiative transfer code (version 6.0), and a coincident Landsat 8 Operational Land Imager (OLI) acquisition on 29 July 2015 during an in-flight radiometric calibration experiment. Airborne remote sensing flights were carried out in northwestern Greenland in preparation for the Ice, Cloud, and land Elevation Satellite 2 (ICESat-2) laser altimeter mission. A total of nine science flights were conducted over the Greenland ice sheet, sea ice, and open-ocean water. The campaign's primary purpose was to correlate green laser pulse penetration into snow and ice with spectroscopic-derived surface properties. An experimental airborne instrument configuration that included a nadir-viewing (looking downward at the surface) non-imaging Analytical Spectral Devices (ASD) Inc. spectrometer that measured upwelling VSWIR (0.35 to 2.5 µm) spectral radiance (Wm −2 sr −1 µm −1 ) in the two-color Slope Imaging Multi-polarization Photon-Counting Lidar's (SIMPL) ground instantaneous field of view, and a zenith-viewing (looking upward at the sky) ASD spectrometer that measured VSWIR spectral irradiance (W m −2 nm −1 ) was flown. National Institute of Standards and Technology (NIST) traceable radiometric calibration procedures for laboratory, in-flight, and field environments are described in detail to achieve a targeted VSWIR measurement requirement of within 5 % to support calibration/validation efforts and remote sensing algorithm development. Our MODTRAN predictions for the 29 July flight line over dark and bright targets indicate that the airborne nadir-viewing spectrometer spectral radiance measurement uncertainty was between 0.6 % and 4.7 % for VSWIR wavelengths (0.4 to 2.0 µm) with atmospheric transmittance greater than 80 %. MODTRAN predictions for Landsat 8 OLI relative spectral response functions suggest that OLI is measuring 6 % to 16 % more top-of-atmosphere (TOA) spectral radiance from the Greenland ice sheet surface than was predicted using apparent reflectance spectra from the nadir-viewing spectrometer. While more investigation is required to convert airborne VSWIR spectral radiance into atmospherically corrected airborne surface reflectance, it is expected that airborne science flight data products will contribute to spectroscopic determination of Greenland ice sheet surface optical properties to improve understanding of their potential influence on ICESat-2 measurements.

Christopher J. Crawford↗

Sources of Sulfate Found in Mounds and Lakes at the Lewis Cliffs Ice Tongue, Transantarctic

Massive but highly localized Na-sulfate mounds (mirabilite, Na2SO4.10H2O) have been found at the terminal moraine of the Lewis Cliffs Ice Tongue (LCIT), Antarctica. (Sigma)34S and (Sigma)18O values of LCIT mirabilite range from +48.8 to +49.3% (CDT), and -16.6 to -17.1% (V-SMOW), respectively, while (Delta)17O average -0.37% (V-SMOW). LCIT mirabilite mounds are isotopically different from other mirabilite mounds found in coastal regions of Antarctica, which have isotope values close to seawater compositions. (Sigma)18O and (Delta)17O values suggest the incorporation of isotopically light glacial water. Data point to initial sulfate formation in an anoxic water body, either as a stratified anoxic deep lake on the surface, a sub-glacial water reservoir, or a sub-glacial lake. Several surface lakes of varying size are also present within this region of the LCIT, and in some cases are adjacent to the mirabilite mounds. O and D isotope compositions of surface lakes confirm they are derived from a mixture of glacial ice and snow that underwent moderate evaporation. (Sigma)18O and (Sigma)D (V-SMOW) values of snow, ice, and lake water range from -64.2 to -29.7%, and -456.0 to -231.7%, respectively. However, the isotope chemistry of these surface lakes is extremely different from the mounds. Dissolved SO4-2 (Sigma)34S and (Sigma)18O values range from +12.0 to +20.0% and -12.8 to -22.2% (the most negative (Sigma)18O of terrestrial sulfate ever reported), respectively, with sulfate (Delta)17O ranging from +0.93 to 2.24%. Ion chromatography data show that lake water is fresh to brackish in origin, with TDS less than 1500 ppm, and sulfate concentration less than 431 ppm. Isotope and chemical data suggest that these lakes are unlikely the source of the mirabilite mounds. We suggest that lake water sulfate is potentially composed of a mixture of atmospheric sulfate and minor components of sulfate of weathering origin, much like the sulfate in the polar plateau soils of the McMurdo Dry Valleys. A simple model explains mirabilite mound formation at the LCIT. Sulfur redox processes could occur sub-glacially as a result of liquid-water-based glacial conditions (Alpine style glacier), most likely formed by pressure melting of overlying ice (Aharon, GCA, 52, 2321-2331). We suggest that the aqueous base of the LCIT contains dissolved SO42- and is anoxic where sulfate reduction to H2S, HS-, or native sulfur takes place. Sulfide is removed by either precipitation as sulfide minerals or by escape of H2S (neither of which have been observed). Mirabilite precipitation is likely the result of evaporation or freezing of sulfate-rich brines as they reach the surface where they manifest themselves as mounds. Pressure from the overlying ice contributing to a pressure-melting scenario that creates the sub-glacial aqueous environment also contributes to the mechanism of upward transport of the sulfate-rich fluids. Further evidence to support this upward transport model comes from the nature of ice motion at the LCIT. Cassidy et al (Meteoritics, 27, 490-525, 1992) pointed out that it is the vertical ice motion in this area that creates the meteorite-stranding surface that could also account for upward transport of sulfate-rich fluids. Alternatively, mirabilite was deposited in a similar condition as present-day coastal Antarctica when the LCIT was wetter and warmer

Socki, Richard↗