An Improved Cryosat-2 Sea Ice Freeboard Retrieval Algorithm Through the Use of Waveform Fitting

We develop an empirical model capable of simulating the mean echo power cross product of CryoSat-2 SAR and SAR In mode waveforms over sea ice covered regions. The model simulations are used to show the importance of variations in the radar backscatter coefficient with incidence angle and surface roughness for the retrieval of surfaceelevation of both sea ice floes and leads. The numerical model is used to fit CryoSat-2 waveforms to enable retrieval of surface elevation through the use of look-up tables and a bounded trust region Newton least squares fitting approach. The use of a model to fit returns from sea ice regions offers advantages over currently used threshold retrackingmethods which are here shown to be sensitive to the combined effect of bandwidth limited range resolution and surface roughness variations. Laxon et al. (2013) have compared ice thickness results from CryoSat-2 and IceBridge, and found good agreement, however consistent assumptions about the snow depth and density of sea ice werenot used in the comparisons. To address this issue, we directly compare ice freeboard and thickness retrievals from the waveform fitting and threshold tracker methods of CryoSat-2 to Operation IceBridge data using a consistent set of parameterizations. For three IceBridge campaign periods from March 20112013, mean differences (CryoSat-2 IceBridge) of 0.144m and 1.351m are respectively found between the freeboard and thickness retrievals using a 50 sea ice floe threshold retracker, while mean differences of 0.019m and 0.182m are found when using the waveform fitting method. This suggests the waveform fitting technique is capable of better reconciling the seaice thickness data record from laser and radar altimetry data sets through the usage of consistent physical assumptions

Model order reduction approaches for infinite horizon optimal control problems via the HJB equation

We investigate feedback control for infinite horizon optimal control problems for partial differential equations. The method is based on the coupling between Hamilton-Jacobi-Bellman (HJB) equations and model reduction techniques. It is well-known that HJB equations suffer the so called curse of dimensionality and, therefore, a reduction of the dimension of the system is mandatory. In this report we focus on the infinite horizon optimal control problem with quadratic cost functionals. We compare several model reduction methods such as Proper Orthogonal Decomposition, Balanced Truncation and a new algebraic Riccati equation based approach. Finally, we present numerical examples and discuss several features of the different methods analyzing advantages and disadvantages of the reduction methods

Rift in Antarctic Glacier: a Unique Chance to Study Ice Shelf Retreat

It happened again, but this time it was caught in the act. During the last week of September 2011 a large transverse rift developed across thefloating terminus of West Antarcticas PineIsland Glacier, less than 5 years after its lastlarge calving event, in 2007 (Figure 1). PineIsland Glaciers retreat has accelerated substantiallyin the past 2 decades, and it is nowlosing 50 gigatons of ice per year, or roughly 25 of Antarcticas total annual contributionto sea level rise [Rignot et al., 2008]. The glaciers recent accelerated retreat is likely triggered by ocean warming and increased submarine melting. As such, it is of significant interest to glaciologists and of heightened societal relevance

Influence of Persistent Wind Scour on the Surface Mass Balance of Antarctica

Accurate quantification of surface snow accumulation over Antarctica is a key constraint for estimates of the Antarctic mass balance, as well as climatic interpretations of ice-core records. Over Antarctica, near-surface winds accelerate down relatively steep surface slopes, eroding and sublimating the snow. This wind scour results in numerous localized regions (< or = 200 sq km) with reduced surface accumulation. Estimates of Antarctic surface mass balance rely on sparse point measurements or coarse atmospheric models that do not capture these local processes, and overestimate the net mass input in wind-scour zones. Here we combine airborne radar observations of unconformable stratigraphic layers with lidar-derived surface roughness measurements to identify extensive wind-scour zones over Dome A, in the interior of East Antarctica. The scour zones are persistent because they are controlled by bedrock topography. On the basis of our Dome A observations, we develop an empirical model to predict wind-scour zones across the Antarctic continent and find that these zones are predominantly located in East Antarctica. We estimate that approx. 2.7-6.6% of the surface area of Antarctica has persistent negative net accumulation due to wind scour, which suggests that, across the continent, the snow mass input is overestimated by 11-36.5 Gt /yr in present surface-mass-balance calculations

New Antarctic gravity anomaly grid for enhanced geodetic and geophysical studies in Antarctica

Gravity surveying is challenging in Antarctica because of its hostile environment and inaccessibility. Nevertheless, many ground-based, airborne, and shipborne gravity campaigns have been completed by the geophysical and geodetic communities since the 1980s. We present the first modern Antarctic-wide gravity data compilation derived from 13 million data points covering an area of 10 million km2, which corresponds to 73% coverage of the continent. The remove-compute-restore technique was applied for gridding, which facilitated leveling of the different gravity data sets with respect to an Earth gravity model derived from satellite data alone. The resulting free-air and Bouguer gravity anomaly grids of 10 km resolution are publicly available. These grids will enable new high-resolution combined Earth gravity models to be derived and represent a major step forward toward solving the geodetic polar data gap problem. They provide a new tool to investigate continental-scale lithospheric structure and geological evolution of Antarctica

Assessment of NASA Airborne Laser Altimetry Data Using Ground-Based GPS Data Near Summit Station, Greenland

A series of NASA airborne lidars have been used in support of satellite laser altimetry missions. These airborne laser altimeters have been deployed for satellite instrument development, for spaceborne data validation, and to bridge the data gap between satellite missions. We used data from ground-based Global Positioning System (GPS) surveys of an 11 km long track near Summit Station, Greenland, to assess the surface–elevation bias and measurement precision of three airborne laser altimeters including the Airborne Topographic Mapper (ATM), the Land, Vegetation, and Ice Sensor (LVIS), and the Multiple Altimeter Beam Experimental Lidar (MABEL). Ground-based GPS data from the monthly ground-based traverses, which commenced in 2006, allowed for the assessment of nine airborne lidar surveys associated with ATM and LVIS between 2007 and 2016. Surface–elevation biases for these altimeters – over the flat, ice-sheet interior – are less than 0.12 m, while assessments of measurement precision are 0.09 m or better. Ground-based GPS positions determined both with and without differential post-processing techniques provided internally consistent solutions. Results from the analyses of ground-based and airborne data provide validation strategy guidance for the Ice, Cloud, and land Elevation Satellite 2 (ICESat-2) elevation and elevation-change data products

Evidence for Extending Anomalous Miocene Volcanism at the Edge of the East Antarctic Craton

Using field observations followed by petrological, geochemical, geochronological, and geophysical data, we infer the presence of a previously unknown Miocene subglacial volcanic center ~230 km from the South Pole. Evidence of volcanism is from boulders of olivine‐bearing amygdaloidal/vesicular basalt and hyaloclastite deposited in a moraine in the southern Transantarctic Mountains. 40Ar/39Ar ages from five specimens plus U‐Pb ages of detrital zircon from glacial till indicate igneous activity 25–17 Ma. The likely source of the volcanism is a circular −735 nT magnetic anomaly 60 km upflow from the sampling site. Subaqueous textures of the volcanics indicate eruption beneath ice or into water at the margin of an ice mass during the early Miocene. These rocks record the southernmost Cenozoic volcanism in Antarctica and expand the known extent of the oldest lavas associated with West Antarctic Rift system. They may be an expression of lithospheric foundering beneath the southern Transantarctic Mountains

Understanding biases in ICESat-2 data due to subsurface scattering using Airborne Topographic Mapper waveform data

The process of laser light reflecting from surfaces made of scattering materials that do not strongly absorb at the wavelength of the laser can involve reflections from hundreds or thousands of individual grains, which can introduce delays in the time between light entering and leaving the surface. These time-of-flight biases depend on the grain size and density of the medium, and thus they can result in spatially and temporally varying surface height biases estimated from laser altimeters, such as NASA's ICESat-2 (Ice Cloud, and land Elevation Satellite-2) mission. Modeling suggests that ICESat-2 might experience a bias difference as large as 0.1–0.2 m between coarse-grained melting snow and fine-grained wintertime snow (Smith et al., 2018), which exceeds the mission's requirement to measure seasonal height differences to an accuracy better than 0.1 m (Markus et al., 2017). In this study, we investigate these biases using a model of subsurface scattering, laser altimetry measurements from NASA's ATM (Airborne Topographic Mapper) system, and grain size estimates based on optical imagery of the ice sheet. We demonstrate that distortions in the shapes of waveforms measured using ATM are related to the optical grain size of the surface estimated using optical reflectance measurements and show that they can be used to estimate an effective grain radius for the surface. Using this effective grain radius as a proxy for the severity of subsurface scattering, we use our model with grain size estimates from optical imagery to simulate corrections for biases in ICESat-2 data due to subsurface scattering and demonstrate that, on the basis of large-scale averages, the corrections calculated based on the satellite optical imagery match the biases in the data. This work demonstrates that waveform-based altimetry data can measure the optical properties of granular surfaces and that corrections based on optical grain size estimates can correct for subsurface-scattering biases in ICESat-2 data.</p

Assessment of NASA airborne laser altimetry data using ground-based GPS data near Summit Station, Greenland

A series of NASA airborne lidars have been used in support of satellite laser altimetry missions. These airborne laser altimeters have been deployed for satellite instrument development, for spaceborne data validation, and to bridge the data gap between satellite missions. We used data from ground-based Global Positioning System (GPS) surveys of an 11 km long track near Summit Station, Greenland, to assess the surface–elevation bias and measurement precision of three airborne laser altimeters including the Airborne Topographic Mapper (ATM), the Land, Vegetation, and Ice Sensor (LVIS), and the Multiple Altimeter Beam Experimental Lidar (MABEL). Ground-based GPS data from the monthly ground-based traverses, which commenced in 2006, allowed for the assessment of nine airborne lidar surveys associated with ATM and LVIS between 2007 and 2016. Surface–elevation biases for these altimeters – over the flat, ice-sheet interior – are less than 0.12 m, while assessments of measurement precision are 0.09 m or better. Ground-based GPS positions determined both with and without differential post-processing techniques provided internally consistent solutions. Results from the analyses of ground-based and airborne data provide validation strategy guidance for the Ice, Cloud, and land Elevation Satellite 2 (ICESat-2) elevation and elevation-change data products