Search

Precipitation is a major component of the hydrologic cycle and plays a significant role in the sea ice mass balance in the polar regions. Over the Southern Ocean, precipitation is particularly uncertain due to the lack of direct observations in this remote and harsh environment. Here we demonstrate that precipitation estimates from eight global reanalyses produce similar spatial patterns between 2000 and 2010, although their annual means vary by about 250 mm yr−1 (or 26% of the median values) and there is little similarity in their representation of interannual variability. ERA-Interim produces the smallest and CFSR produces the largest amount of precipitation overall. Rainfall and snowfall are partitioned in five reanalyses; snowfall suffers from the same issues as the total precipitation comparison, with ERA-Interim producing about 128 mm less snowfall and JRA-55 about 103 mm more rainfall compared to the other reanalyses. When compared to CloudSat-derived snowfall, these five reanalyses indicate similar spatial patterns, but differ in their magnitude. All reanalyses indicate precipitation on nearly every day of the year, with spurious values occurring on an average of about 60 days yr−1, resulting in an accumulation of about 4.5 mm yr−1. While similarities in spatial patterns among the reanalyses suggest a convergence, the large spread in magnitudes points to issues with the background models in adequately reproducing precipitation rates, and the differences in the model physics employed. Further improvements to model physics are required to achieve confidence in precipitation rate, as well as the phase and frequency of precipitation in these products.

Taking stock of our losses Earth's ice sheets are melting and sea levels are rising, so it behooves us to understand better which climate processes are responsible for how much of the mass loss. Smith et al. estimated grounded and floating ice mass change for the Greenland and Antarctic ice sheets from 2003 to 2019 using satellite laser altimetry data from NASA's ICESat and ICESat-2 satellites. They show how changing ice flow, melting, and precipitation affect different regions of ice and estimate that grounded-ice loss averaged close to 320 gigatons per year over that period and contributed 14 millimeters to sea level rise. Science , this issue p. 1239

NASA’s Multiple Altimeter Beam Experimental LiDAR (MABEL) is an aircraft-based photon-counting laser altimeter designed as a simulator to test measurement techniques and algorithms for Advanced Topographic Laser Altimeter System (ATLAS), the sole instrument on NASA’s Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) mission. By measuring the time of flight, pointing angle, and absolute position for individual photons, ICESat-2 provides detailed elevation measurements of earth’s surface. Calculating accurate and precise elevations requires an understanding of how photons interact with surfaces, and characterization of the photon distribution after returning from surfaces. Neither MABEL nor ATLAS records the transmitted laser pulse shape, relying instead on aggregating several pulses worth of photons, often using histograms, to characterize the pulse shape. In this paper, we assess the limitations of using histograms and propose a more robust method to describe MABEL’s system impulse-response function using an exponentially modified Gaussian distribution. We also provide standard error estimates for the arithmetic mean and standard deviation calculations, and for exponentially modified Gaussian parameters using a Monte Carlo sensitivity analysis. We apply this method to photon returns from a sea ice lead and from a dry salt lake bed as case studies for estimating the standard error associated with sample size for the arithmetic mean and standard deviation, and for the exponentially modified Gaussian parameters. We use these standard errors to calculate the minimum number of photons required to find both Gaussian and exponentially modified Gaussian distribution parameters within 3 cm of their parent population values.

The NASA Eulerian Snow On Sea Ice Model (NESOSIM) is a new, open-source snow budget model that is currently configured to produce daily estimates of the depth and density of snow on sea ice across the Arctic Ocean through the accumulation season. NESOSIM has been developed in a three-dimensional Eulerian framework and includes two (vertical) snow layers and several simple parameterizations (accumulation, wind packing, advection–divergence, blowing snow lost to leads) to represent key sources and sinks of snow on sea ice. The model is forced with daily inputs of snowfall and near-surface winds (from reanalyses), sea ice concentration (from satellite passive microwave data) and sea ice drift (from satellite feature tracking) during the accumulation season (August through April). In this study, we present the NESOSIM formulation, calibration efforts, sensitivity studies and validation efforts across an Arctic Ocean domain (100 km horizontal resolution). The simulated snow depth and density are calibrated with in situ data collected on drifting ice stations during the 1980s. NESOSIM shows strong agreement with the in situ seasonal cycles of snow depth and density, and shows good (moderate) agreement with the regional snow depth (density) distributions. NESOSIM is run for a contemporary period (2000 to 2015), with the results showing strong sensitivity to the reanalysis-derived snowfall forcing data, with the Modern-Era Retrospective analysis for Research and Applications (MERRA) and the Japanese Meteorological Agency 55-year reanalysis (JRA-55) forced snow depths generally higher than ERA-Interim, and the Arctic System Reanalysis (ASR) generally lower. We also generate and force NESOSIM with a consensus median daily snowfall dataset from these reanalyses. The results are compared against snow depth estimates derived from NASA's Operation IceBridge (OIB) snow radar data from 2009 to 2015, showing moderate–strong correlations and root mean squared errors of ∼ 10 cm depending on the OIB snow depth product analyzed, similar to the comparisons between OIB snow depths and the commonly used modified Warren snow depth climatology. Potential improvements to this initial NESOSIM formulation are discussed in the hopes of improving the accuracy and reliability of these simulated snow depths and densities.

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.

The potential of deriving snow depth estimates using differences in freeboard heights from CryoSat-2 (CS-2) and ICESat-2 (IS-2) is examined. In our analysis, we use lidar freeboard from the Airborne Topographic Mapper (ATM) on Operation IceBridge (OIB) as proxy of IS-2 total (snow + ice) freeboard. Snow depths are estimates from the OIB snow radar. Differences in height between the total (ATM) and ice (CS-2) freeboards are related to snow depth by the refractive index of the snow layer ( η s ) , which is dependent on snow density. For two years (2014 and 2015), regression of the ATM and CS-2 freeboard differences against OIB snow depth gives correlations of ∼0.80, estimated η s of ∼1.21, and standard errors of ∼8 cm. The resulting refractive index, η s , can be compared to that expected of the Arctic snow cover in early spring (1.25 ± 0.05). The expected biases and variability in the regression analysis are discussed. Results suggest that snow depth can be estimated from the freeboard differences. The benefits of adjusting the orbit of CS-2 for providing more optimized overlaps between IS-2 and CS-2 are considered.

NASA's Ice, Cloud, and land Elevation Satellite (ICESat-2), to be launched in 2018, will measure the height of Earth from space using lasers, collecting the most precise and detailed account yet of our planet's elevation. The mission will allow scientists to investigate how global warming is changing the planet's icy polar regions and to take stock of Earth's vegetation. ICESat-2's emphasis on polar ice, as well as its unique measurement approach, will provide an intriguing and accessible focus for the mission's education and outreach programs. Sea ice and land ice are areas that have experienced significant change in recent years. It is key to communicate why we are measuring these areas and their importance. ICESat-2 science data will provide much-needed answers to climate change questions such as, “Is the ice really melting in the polar regions?” and “What does studying Earth's frozen regions tell us about our changing climate?” In this paper, lessons-learned and novel techniques for engaging and educating all audiences in the mission will be discussed, such as including results of a unique collaboration with art design school the Savannah College of Art Design (SCAD) to create fun and exciting products such as animated characters and interactive stories. Future collaborations with wildlife researchers, a new citizen science program in collaboration with GLOBE, and evidence from other STEAM (Science, Technology, Engineering, Arts, Math) education approaches will also be detailed in this paper.

Sea-ice thickness in the Sea of Okhotsk is estimated for 2004–2008 from ICESat derived freeboard under the assumption of hydrostatic balance. Total ice thickness including snow depth (htot) averaged over 2004–2008 is 95 cm. The interannual variability of htot is large; from 77.5 cm (2008) to 110.4 cm (2005). The mode of htot varies from 50–60 cm (2007 and 2008) to 70–80 cm (2005). Ice thickness derived from ICESat data is validated from a comparison with that observed by Electromagnetic Induction Instrument (EM) aboard the icebreaker Soya near Hokkaido, Japan. Annual maps of htot reveal that the spatial distribution of htot is similar every year. Ice volume of 6.3 × 1011 m3 is estimated from the ICESat derived htot and AMSR-E derived ice concentration. A comparison with ice area demonstrates that the ice volume cannot always be represented by the area solely, despite the fact that the area has been used as a proxy of the volume in the Sea of Okhotsk. The ice volume roughly corresponds to that of annual ice production in the major coastal polynyas estimated based on heat budget calculations. This also supports the validity of the estimation of sea-ice thickness and volume using ICESat data.