Title: Weekly Mapping of Sea Ice Freeboard in the Ross Sea from ICESat-2
NASA’s ICESat-2 has been providing sea ice freeboard measurements across the polar regions since October 2018. In spite of the outstanding spatial resolution and precision of ICESat-2, the spatial sparsity of the data can be a critical issue for sea ice monitoring. This study employs a geostatistical approach (i.e., ordinary kriging) to characterize the spatial autocorrelation of the ICESat-2 freeboard measurements (ATL10) to estimate weekly freeboard variations in 2019 for the entire Ross Sea area, including where ICESat-2 tracks are not directly available. Three variogram models (exponential, Gaussian, and spherical) are compared in this study. According to the cross-validation results, the kriging-estimated freeboards show correlation coefficients of 0.56–0.57, root mean square error (RMSE) of ~0.12 m, and mean absolute error (MAE) of ~0.07 m with the actual ATL10 freeboard measurements. In addition, the estimated errors of the kriging interpolation are low in autumn and high in winter to spring, and low in southern regions and high in northern regions of the Ross Sea. The effective ranges of the variograms are 5–10 km and the results from the three variogram models do not show significant differences with each other. The southwest (SW) sector of the Ross Sea shows low and consistent freeboard over the entire year because of the frequent opening of wide polynya areas generating new ice in this sector. However, the southeast (SE) sector shows large variations in freeboard, which demonstrates the advection of thick multiyear ice from the Amundsen Sea into the Ross Sea. Thus, this kriging-based interpolation of ICESat-2 freeboard can be used in the future to estimate accurate sea ice production over the Ross Sea by incorporating other remote sensing data.  more » « less
Award ID(s):
1835784
PAR ID:
10308731
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Remote Sensing
Volume:
13
Issue:
16
ISSN:
2072-4292
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract. Ocean–sea-ice coupled models constrained by various observations provide different ice thickness estimates in the Antarctic. We evaluatecontemporary monthly ice thickness from four reanalyses in the Weddell Sea: the German contribution of the project Estimating the Circulation and Climate ofthe Ocean Version 2 (GECCO2), the Southern Ocean State Estimate (SOSE), the Ensemble Kalman Filter system based on the Nucleus for European Modelling of the Ocean (NEMO-EnKF) and the Global Ice–Ocean Modeling and Assimilation System (GIOMAS). The evaluation is performed againstreference satellite and in situ observations from ICESat-1, Envisat, upward-looking sonars and visual ship-based sea-ice observations. Compared withICESat-1, NEMO-EnKF has the highest correlation coefficient (CC) of 0.54 and lowest root mean square error (RMSE) of 0.44 m. Compared within situ observations, SOSE has the highest CC of 0.77 and lowest RMSE of 0.72 m. All reanalyses underestimate ice thickness near the coast ofthe western Weddell Sea with respect to ICESat-1 and in situ observations even though these observational estimates may be biased low. GECCO2 andNEMO-EnKF reproduce the seasonal variation in first-year ice thickness reasonably well in the eastern Weddell Sea. In contrast, GIOMAS ice thicknessperforms best in the central Weddell Sea, while SOSE ice thickness agrees most with the observations from the southern coast of the Weddell Sea. Inaddition, only NEMO-EnKF can reproduce the seasonal evolution of the large-scale spatial distribution of ice thickness, characterized by the thickice shifting from the southwestern and western Weddell Sea in summer to the western and northwestern Weddell Sea in spring. We infer that the thickice distribution is correlated with its better simulation of northward ice motion in the western Weddell Sea. These results demonstrate thepossibilities and limitations of using current sea-ice reanalysis for understanding the recent variability of sea-ice volume in the Antarctic. 
    more » « less
  2. Assimilation of remote-sensing products of sea ice thickness (SIT) into sea ice–ocean models has been shown to improve the quality of sea ice forecasts. Key open questions are whether assimilation of lower-level data products such as radar freeboard (RFB) can further improve model performance and what performance gains can be achieved through joint assimilation of these data products in combination with a snow depth product. The Arctic Mission Benefit Analysis system was developed to address this type of question. Using the quantitative network design (QND) approach, the system can evaluate, in a mathematically rigorous fashion, the observational constraints imposed by individual and groups of data products. We demonstrate the approach by presenting assessments of the observation impact (added value) of different Earth observation (EO) products in terms of the uncertainty reduction in a 4-week forecast of sea ice volume (SIV) and snow volume (SNV) for three regions along the Northern Sea Route in May 2015 using a coupled model of the sea ice–ocean system, specifically the Max Planck Institute Ocean Model. We assess seven satellite products: three real products and four hypothetical products. The real products are monthly SIT, sea ice freeboard (SIFB), and RFB, all derived from CryoSat-2 by the AlfredWegener Institute. These are complemented by two hypothetical monthly laser freeboard (LFB) products with low and high accuracy, as well as two hypothetical monthly snow depth products with low and high accuracy. On the basis of the per-pixel uncertainty ranges provided with the CryoSat-2 SIT, SIFB, and RFB products, the SIT and RFB achieve a much better performance for SIV than the SIFB product. For SNV, the performance of SIT is only low, the performance of SIFB is higher and the performance of RFB is yet higher. A hypothetical LFB product with low accuracy (20 cm uncertainty) falls between SIFB and RFB in performance for both SIV and SNV. A reduction in the uncertainty of the LFB product to 2 cm yields a significant increase in performance. Combining either of the SIT or freeboard products with a hypothetical snow depth product achieves a significant performance increase. The uncertainty in the snow product matters: a higher-accuracy product achieves an extra performance gain. Providing spatial and temporal uncertainty correlations with the EO products would be beneficial not only for QND assessments, but also for assimilation of the products. 
    more » « less
  3. null (Ed.)
    The measurement of sea ice elevation above sea level or the “freeboard” depends upon an accurate retrieval of the local sea level. The local sea level has been previously retrieved from altimetry data alone by the lowest elevation method, where the percentage of the lowest elevations over a particular segment length scale was used. Here, we provide an evaluation of the scale dependence on these local sea level retrievals using data from NASA Operation IceBridge (OIB) which took place in the Ross Sea in 2013. This is a unique dataset of laser altimeter measurements over five tracks from the Airborne Topographic Mapper (ATM), with coincidently high-spatial resolution images from the Digital Mapping System (DMS), that allows for an independent sea level validation. The local sea level is first calculated by using the mean elevation of ATM L1B data over leads identified by using the corresponding DMS imagery. The resulting local sea level reference is then used as ground truth to validate the local sea levels retrieved from ATM L2 by using nine different percentages of the lowest elevation (0.1%, 0.5%, 1%, 1.5%, 2%, 2.5%, 3%, 3.5%, and 4%) at seven different segment length scales (1, 5, 10, 15, 20, 25, and 50 km) for each of the five ATM tracks. The closeness to the 1:1 line, R2, and root mean square error (RMSE) is used to quantify the accuracy of the retrievals. It is found that all linear least square fits are statistically significant (p < 0.05) using an F test at every scale for all tested data. In general, the sea level retrievals are farther away from the 1:1 line when the segment length scale increases from 1 or 5 to 50 km. We find that the retrieval accuracy is affected more by the segment length scale than the percentage scale. Based on our results, most retrievals underestimate the local sea level; the longer the segment length (from 1 to 50 km) used, especially at small percentage scales, the larger the error tends to be. The best local sea level based on a higher R2 and smaller RMSE for all the tracks combined is retrieved by using 0.1–2% of the lowest elevations at the 1–5 km segment lengths. 
    more » « less
  4. Abstract The West Antarctic Ice Sheet (WAIS) contains an ice volumetric equivalent of 5.3 m of global mean sea‐level rise. Understanding the climate mechanisms that influence WAIS surface mass balance (SMB) can help reduce uncertainties in sea‐level rise projections. Previous work has focused on trends related to temperature and the roles of well‐known climate phenomena in the Pacific sector, including the El Niño Southern Oscillation (ENSO) and the Southern Annular Mode. Here, the focus is on interannual variability of precipitation associated with Rossby wave trains emanating from tropical convection regions in the Indian Ocean and the central tropical Pacific. Objective clustering of SMB measurements from 25 WAIS ice and firn cores from 1900 through 1999 indicated oppositely signed mean SMB trends in the East and West sectors of WAIS, and these trends were removed to focus on year‐to‐year variations. Climatically important months for each sector were identified by correlating reanalysis precipitation with ice core SMB. Tropical Rossby wave source regions that produced anomalous WAIS SMB were identified by analyzing composite convection, geopotential height, and wave activity flux. East WAIS SMB variability was controlled by anomalously low geopotential height over the Ross and Amundsen Seas. West WAIS SMB variability was driven by ENSO and non‐uniformity in the phase distribution of the Madden‐Julian Oscillation. These results provide a deeper understanding of the drivers of WAIS SMB variability and may improve projections of future sea‐level rise. 
    more » « less
  5. As climate warms and the transition from a perennial to a seasonal Arctic sea-ice cover is imminent, understanding melt ponding is central to understanding changes in the new Arctic. National Aeronautics and Space Administration (NASA)’s Ice, Cloud and land Elevation Satellite (ICESat-2) has the capacity to provide measurements and monitoring of the onset of melt in the Arctic and on melt progression. Yet ponds are currently not identified on the ICESat-2 standard sea-ice products, in which only a single surface is determined. The objective of this article is to introduce a mathematical algorithm that facilitates automated detection of melt ponds in the ICESat-2 Advanced Topographic Laser Altimeter System (ATLAS) data, retrieval of two surface heights, pond surface and bottom, and measurements of depth and width of melt ponds. With ATLAS, ICESat-2 carries the first spaceborne multibeam micropulse photon-counting laser altimeter system, operating at 532-nm frequency. ATLAS data are recorded as clouds of discrete photon points. The Density-Dimension Algorithm for bifurcating sea-ice reflectors (DDA-bifurcate-seaice) is an autoadaptive algorithm that solves the problem of pond detection near the 0.7-m nominal along-track spacing of ATLAS data, utilizing the radial basis function for calculation of a density field and a threshold function that automatically adapts to changes in the background, apparent surface reflectance, and some instrument effects. The DDA-bifurcate-seaice is applied to large ICESat-2 datasets from the 2019 and 2020 melt seasons in the multiyear Arctic sea-ice region. Results are evaluated by comparison with those from a manually forced algorithm. 
    more » « less