Search

The Greenland and Antarctic ice sheets are important indicators of climate change and major contributors to sea level rise. Hence, precise, long-term observations of ice mass change are required to assess their contribution to sea level rise. Such observations can be achieved through three different methods. They can be achieved directly by measuring regional changes in the Earth's gravity field using the Gravity Recovery and Climate Experiment Follow-On (GRACE-FO) satellite system. Alternatively, they can be achieved indirectly by measuring changes in ice thickness using satellite altimetry or by estimating changes in the mass budget using a combination of regional climate model data and ice discharge across the grounding line, based on multi-sensor satellite radar observations of ice velocity (Hanna et al., 2013). Satellite radar altimetry has been used to measure elevation change since 1992 through a combination of various missions. In addition to the surface slope and complex topography, it has been shown that one of the most challenging issues concerns spatial and temporal variability in radar pulse penetration into the snowpack. This results in an inaccurate measurement of the true surface elevation and consequently affects surface elevation change (SEC) estimates. To increase the accuracy of surface elevation measurements retrieved by retracking the radar return waveform and thus reduce the uncertainty in the SEC, we developed a deep convolutional-neural-network architecture (AWI-ICENet1). AWI-ICENet1 is trained using a simulated reference data set with 3.8 million waveforms, taking into account different surface slopes, topography, and attenuation. The successfully trained network is finally applied as an AWI-ICENet1 retracker to the full time series of CryoSat-2 Low Resolution Mode (LRM) waveforms over both ice sheets. We compare the AWI-ICENet1-retrieved SEC with estimates from conventional retrackers, including the threshold first-maximum retracker algorithm (TFMRA) and the European Space Agency's (ESA) ICE1 and ICE2 products. Our results show less uncertainty and a great decrease in the effect of time-variable radar penetration, reducing the need for corrections based on its close relationship with backscatter and/or leading-edge width, which are typically used in SEC processing. This technique provides new opportunities to utilize convolutional neural networks in the processing of satellite altimetry data and is thus applicable to historical, recent, and future missions.

The 79° North Glacier (Nioghalvfjerdsbrae, 79NG) is one of three remaining glaciers with a floating tongue in Greenland. Although the glacier has been considered exceptionally stable in the past, earlier studies have shown that the ice tongue has thinned in recent decades. By conducting high-resolution ground-based and airborne radar measurements in conjunction with satellite remote-sensing observations, we find significant changes in the geometry of 79NG. In the vicinity of the grounding line, a 500 m high subglacial channel has grown since ∼ 2010 and has caused surface lowering of up to 7.6 m a−1. Our results show extreme basal melt rates exceeding 150 m a−1 over a period of 17 d within a distance of 5 km from the grounding line, where the ice has thinned by 32 % since 1998. We find a heterogeneous distribution of melt rates, likely due to variability in water column thickness and channelization of the ice base. Time series of melt rates show a decrease in basal melting since 2018, indicating an inflow of colder water into the cavity below 79NG. We discuss the processes that have led to the changes in geometry and conclude that the inflow of warm ocean currents has led to the extensive thinning of 79NG's floating ice tongue near the grounding line over the last 2 decades. In contrast, we hypothesize that the growth of the channel results from increased subglacial discharge due to a considerably enlarged area of summer surface melt due to the warming of the atmosphere.

A detailed understanding of how the Antarctic ice sheet (AIS) responds to a warming climate is needed because it will most likely increase the rate of global mean sea level rise. Time-variable satellite gravimetry, realized by the Gravity Recovery and Climate Experiment (GRACE) and Gravity Recovery and Climate Experiment Follow-On (GRACE-FO) missions, is directly sensitive to AIS mass changes. However, gravimetric mass balances are subject to two major limitations. First, the usual correction of the glacial isostatic adjustment (GIA) effect by modelling results is a dominant source of uncertainty. Second, satellite gravimetry allows for a resolution of a few hundred kilometres only, which is insufficient to thoroughly explore causes of AIS imbalance. We have overcome both limitations by the first global inversion of data from GRACE and GRACE-FO, satellite altimetry (CryoSat-2), regional climate modelling (RACMO2), and firn densification modelling (IMAU-FDM). The inversion spatially resolves GIA in Antarctica independently from GIA modelling jointly with changes of ice mass and firn air content at 50 km resolution. We find an AIS mass balance of −144 ± 27 Gt a−1 from January 2011 to December 2020. This estimate is the same, within uncertainties, as the statistical analysis of 23 different mass balances evaluated in the Ice sheet Mass Balance Inter-comparison Exercise (IMBIE; Otosaka et al., 2023b). The co-estimated GIA corresponds to an integrated mass effect of 86 ± 21 Gt a−1 over Antarctica, and it fits better with global navigation satellite system (GNSS) results than other GIA predictions. From propagating covariances to integrals, we find a correlation coefficient of −0.97 between the AIS mass balance and the GIA estimate. Sensitivity tests with alternative input data sets lead to results within assessed uncertainties.

The European Beyond EPICA project aims to extract a continuous ice core of up to 1.5 Ma, with a maximum age density of 20 kyr m−1 at Little Dome C (LDC). We present a 1D numerical model which calculates the age of the ice around Dome C. The model inverts for basal conditions and accounts either for melting or for a layer of stagnant ice above the bedrock. It is constrained by internal reflecting horizons traced in radargrams and dated using the EPICA Dome C (EDC) ice core age profile. We used three different radar datasets ranging from a 10 000 km2 airborne survey down to 5 km long ground-based radar transects over LDC. We find that stagnant ice exists in many places, including above the LDC relief where the new Beyond EPICA drill site (BELDC) is located. The modelled thickness of this layer of stagnant ice roughly corresponds to the thickness of the basal unit observed in one of the radar surveys and in the autonomous phase-sensitive radio-echo sounder (ApRES) dataset. At BELDC, the modelled stagnant ice thickness is 198±44 m and the modelled oldest age of ice is 1.45±0.16 Ma at a depth of 2494±30 m. This is very similar to all sites situated on the LDC relief, including that of the Million Year Ice Core project being conducted by the Australian Antarctic Division. The model was also applied to radar data in the area 10–15 km north of EDC (North Patch), where we find either a thin layer of stagnant ice (generally m) or a negligible melt rate ( mm yr−1). The modelled maximum age at North Patch is over 2 Ma in most places, with ice at 1.5 Ma having a resolution of 9–12 kyr m−1, making it an exciting prospect for a future Oldest Ice drill site.

The largest floating tongue of Greenland’s ice sheet, Nioghalvfjerdsbræ, has been relatively stable with respect to areal retreat until 2022. Draining more than 6 % of the ice sheet, a disintegration of Nioghalvfjerdsbræ's floating tongue and subsequent acceleration due to loss in buttressing are likely to lead to sea level rise. Therefore, the stability of the floating tongue is a focus of this study. We employed a suite of observational methods to detect recent changes at the calving front. We found that the calving style has changed since 2016 at the southern part of the eastern calving front, from tongue-type calving to a crack evolution initiated at frontal ice rises reaching 5–7 km and progressing further upstream compared to 2010. The calving front area is further weakened by an area upstream of the main calving front that consists of open water and an ice mélange that has substantially expanded, leading to the formation of a narrow ice bridge. These geometric and mechanical changes may be a precursor of instability of the floating tongue. We complement our study by numerical ice flow simulations to estimate the impact of future ice-front retreat and complete ice shelf disintegration on the discharge of grounded ice. These idealized scenarios reveal that a loss of the south-eastern area of the ice shelf would lead to a 0.2 % increase in ice discharge at the grounding line, while a sudden collapse of the frontal area (46 % of the floating tongue area) will enhance the ice discharge by 5.1 % due to loss in buttressing. Eventually, a full collapse of the floating tongue increases the grounding line flux by 166 %.

The Northeast Greenland Ice Stream has recently seen significant change to its floating margins and has been identified as vulnerable to future climate warming. Inflow of warm Atlantic Intermediate Water (AIW) from the continental shelf has been observed in the vicinity of the Nioghalvfjerdsfjorden (79∘ N) Glacier calving front, but AIW penetration deep into the ice shelf cavity has not been observed directly. Here, we report temperature and salinity measurements from profiles in an epishelf lake, which provide the first direct evidence of AIW proximal to the grounding line of 79∘ N Glacier, over 50 km from the calving front. We also report evidence for partial un-grounding of the margin of 79∘ N Glacier taking place at the western end of the epishelf lake. Comparison of our measurements to those close to the calving front shows that AIW transits the cavity to reach the grounding line within a few months. The observations provide support for modelling studies that infer AIW-driven basal melt proximal to the grounding line and demonstrate that offshore oceanographic changes can be rapidly transmitted throughout the sub-ice-shelf cavity, with implications for near-future stability of the ice stream.

Ice losses from the Greenland and Antarctic ice sheets have accelerated since the 1990s, accounting for a significant increase in the global mean sea level. Here, we present a new 29-year record of ice sheet mass balance from 1992 to 2020 from the Ice Sheet Mass Balance Inter-comparison Exercise (IMBIE). We compare and combine 50 independent estimates of ice sheet mass balance derived from satellite observations of temporal changes in ice sheet flow, in ice sheet volume, and in Earth's gravity field. Between 1992 and 2020, the ice sheets contributed 21.0±1.9 mm to global mean sea level, with the rate of mass loss rising from 105 Gt yr−1 between 1992 and 1996 to 372 Gt yr−1 between 2016 and 2020. In Greenland, the rate of mass loss is 169±9 Gt yr−1 between 1992 and 2020, but there are large inter-annual variations in mass balance, with mass loss ranging from 86 Gt yr−1 in 2017 to 444 Gt yr−1 in 2019 due to large variability in surface mass balance. In Antarctica, ice losses continue to be dominated by mass loss from West Antarctica (82±9 Gt yr−1) and, to a lesser extent, from the Antarctic Peninsula (13±5 Gt yr−1). East Antarctica remains close to a state of balance, with a small gain of 3±15 Gt yr−1, but is the most uncertain component of Antarctica's mass balance. The dataset is publicly available at https://doi.org/10.5285/77B64C55-7166-4A06-9DEF-2E400398E452 (IMBIE Team, 2021).

Ice shelves play a key role in the stability of the Antarctic Ice Sheet due to their buttressing effect. A loss of buttressing as a result of increased basal melting or ice shelf disintegration will lead to increased ice discharge. Some ice shelves exhibit channels at the base that are not yet fully understood. In this study, we present in situ melt rates of a channel which is up to 330 m high and located in the southern Filchner Ice Shelf. Maximum observed melt rates are 2 m yr−1. Melt rates inside the channel decrease in the direction of ice flow and turn to freezing ∼55 km downstream of the grounding line. While closer to the grounding line melt rates are higher within the channel than outside, this relationship reverses further downstream. Comparing the modeled evolution of this channel under present-day climate conditions over 250 years with its present geometry reveals a mismatch. Melt rates twice as large as the present-day values are required to fit the observed geometry. In contrast, forcing the model with present-day melt rates results in a closure of the channel, which contradicts observations. The ice shelf experiences strong tidal variability in vertical strain rates at the measured site, and discrete pulses of increased melting occurred throughout the measurement period. The type of melt channel in this study diminishes in height with distance from the grounding line and is hence not a destabilizing factor for ice shelves.

Abstract We use airborne accumulation radar data acquired over the Greenland Ice Sheet between 2002 and 2018 to identify changes in ice slab extent and thickness. We show that ice slabs several meters thick were already present at least as early as 2002. Between 2012 and 2018, they expanded by 13,400–17,600 km2 inland, or 37%–44%. Our results document that the extremely warm summer of 2012 produced near‐surface ice layers at higher elevations, enabling ice slabs to develop with only moderate melting in the following summers. With repeat flights along a transect in southwest Greenland, we show that moderate melting primarily causes slab thickening through uniform accretion on top of the ice slabs, while large melting events can also trigger localized accretion below existing ice slabs.

We present a high-resolution airborne radar data set (EGRIP-NOR-2018) for the onset region of the Northeast Greenland Ice Stream (NEGIS). The radar data were acquired in May 2018 with the Alfred Wegener Institute's multichannel ultra-wideband (UWB) radar mounted on the Polar 6 aircraft. Radar profiles cover an area of ∼24 000 km2 and extend over the well-defined shear margins of the NEGIS. The survey area is centered at the location of the drill site of the East Greenland Ice-Core Project (EastGRIP), and several radar lines intersect at this location. The survey layout was designed to (i) map the stratigraphic signature of the shear margins with radar profiles aligned perpendicular to ice flow, (ii) trace the radar stratigraphy along several flow lines, and (iii) provide spatial coverage of ice thickness and basal properties. While we are able to resolve radar reflections in the deep stratigraphy, we cannot fully resolve the steeply inclined reflections at the tightly folded shear margins in the lower part of the ice column. The NEGIS is causing the most significant discrepancies between numerically modeled and observed ice surface velocities. Given the high likelihood of future climate and ocean warming, this extensive data set of new high-resolution radar data in combination with the EastGRIP ice core will be a key contribution to understand the past and future dynamics of the NEGIS. The EGRIP-NOR-2018 radar data products can be obtained from the PANGAEA data publisher (https://doi.pangaea.de/10.1594/PANGAEA.928569; Franke et al., 2021a).

We derive recent surface mass balance (SMB) estimates from airborne radar observations along the iSTAR traverse (2013, 2014) at Pine Island Glacier (PIG), West Antarctica. Ground-based neutron probe measurements provide information of snow and firn density with depth at 22 locations and were used to date internal annual reflection layers. The 2005 layer was traced for a total distance of 2367 km to determine annual mean SMB for the period 2005–2014. Using complementary SMB estimates from two regional climate models, RACMO2.3p2 and MAR, and a geostatistical kriging scheme, we determine a regional-scale SMB distribution with similar main characteristics to that determined for the period 1985–2009 in previous studies. Local departures exist for the northern PIG slopes, where the orographic precipitation shadow effect appears to be more pronounced in our observations, and the southward interior, where the SMB gradient is more pronounced in previous studies. We derive total mass inputs of 79.9±19.2 and 82.1±19.2 Gt yr−1 to the PIG basin based on complementary ASIRAS–RACMO and ASIRAS–MAR SMB estimates, respectively. These are not significantly different to the value of 78.3±6.8 Gt yr−1 for the period 1985–2009. Thus, there is no evidence of a secular trend at decadal scales in total mass input to the PIG basin. We note, however, that our estimated uncertainty is more than twice the uncertainty for the 1985–2009 estimate on total mass input. Our error analysis indicates that uncertainty estimates on total mass input are highly sensitive to the selected krige methodology and assumptions made on the interpolation error, which we identify as the main cause for the increased uncertainty range compared to the 1985–2009 estimates.

Curvilinear channels on the surface of an ice shelf indicate the presence of large channels at the base. Modelling studies have shown that where these surface expressions intersect the grounding line, they coincide with the likely outflow of subglacial water. An understanding of the initiation and the ice–ocean evolution of the basal channels is required to understand the present behaviour and future dynamics of ice sheets and ice shelves. Here, we present focused active seismic and radar surveys of a basal channel, ∼950 m wide and ∼200 m high, and its upstream continuation beneath Support Force Glacier, which feeds into the Filchner Ice Shelf, West Antarctica. Immediately seaward from the grounding line, below the basal channel, the seismic profiles show an ∼6.75 km long, 3.2 km wide and 200 m thick sedimentary sequence with chaotic to weakly stratified reflections we interpret as a grounding line fan deposited by a subglacial drainage channel directly upstream of the basal channel. Further downstream the seabed has a different character; it consists of harder, stratified consolidated sediments, deposited under different glaciological circumstances, or possibly bedrock. In contrast to the standard perception of a rapid change in ice shelf thickness just downstream of the grounding line, we find a flat topography of the ice shelf base with an almost constant ice thickness gradient along-flow, indicating only little basal melting, but an initial widening of the basal channel, which we ascribe to melting along its flanks. Our findings provide a detailed view of a more complex interaction between the ocean and subglacial hydrology to form basal channels in ice shelves.

The surface reflection of solar radiation comprises an important boundary condition for solar radiative transfer simulations. In polar regions above snow surfaces, the surface reflection is particularly anisotropic due to low Sun elevations and the highly anisotropic scattering phase function of the snow crystals. The characterization of this surface reflection anisotropy is essential for satellite remote sensing over both the Arctic and Antarctica. To quantify the angular snow reflection properties, the hemispherical-directional reflectance factor (HDRF) of snow surfaces was derived from airborne measurements in Antarctica during austral summer in 2013/14. For this purpose, a digital 180∘ fish-eye camera (green channel, 490–585 nm wavelength band) was used. The HDRF was measured for different surface roughness conditions, optical-equivalent snow grain sizes, and solar zenith angles. The airborne observations covered an area of around 1000 km × 1000 km in the vicinity of Kohnen Station (75∘0′ S, 0∘4′ E) at the outer part of the East Antarctic Plateau. The observations include regions with higher (coastal areas) and lower (inner Antarctica) precipitation amounts and frequencies. The digital camera provided upward, angular-dependent radiance measurements from the lower hemisphere. The comparison of the measured HDRF derived for smooth and rough snow surfaces (sastrugi) showed significant differences, which are superimposed on the diurnal cycle. By inverting a semi-empirical kernel-driven bidirectional reflectance distribution function (BRDF) model, the measured HDRF of snow surfaces was parameterized as a function of solar zenith angle, surface roughness, and optical-equivalent snow grain size. This allows a direct comparison of the HDRF measurements with the BRDF derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite product MCD43. For the analyzed cases, MODIS observations (545–565 nm wavelength band) generally underestimated the anisotropy of the surface reflection. The largest deviations were found for the volumetric model weight fvol (average underestimation by a factor of 10). These deviations are likely linked to short-term changes in snow properties.

The ESA Earth Explorer CryoSat-2 was launched on 8 April 2010 to monitor the precise changes in the thickness of terrestrial ice sheets and marine floating ice. To do that, CryoSat orbits the planet at an altitude of around 720 km with a retrograde orbit inclination of 92∘ and a quasi repeat cycle of 369 d (30 d subcycle). To reach the mission goals, the CryoSat products have to meet the highest quality standards to date, achieved through continual improvements of the operational processing chains. The new CryoSat Ice Baseline-D, in operation since 27 May 2019, represents a major processor upgrade with respect to the previous Ice Baseline-C. Over land ice the new Baseline-D provides better results with respect to the previous baseline when comparing the data to a reference elevation model over the Austfonna ice cap region, improving the ascending and descending crossover statistics from 1.9 to 0.1 m. The improved processing of the star tracker measurements implemented in Baseline-D has led to a reduction in the standard deviation of the point-to-point comparison with the previous star tracker processing method implemented in Baseline-C from 3.8 to 3.7 m. Over sea ice, Baseline-D improves the quality of the retrieved heights inside and at the boundaries of the synthetic aperture radar interferometric (SARIn or SIN) acquisition mask, removing the negative freeboard pattern which is beneficial not only for freeboard retrieval but also for any application that exploits the phase information from SARIn Level 1B (L1B) products. In addition, scatter comparisons with the Beaufort Gyre Exploration Project (BGEP; https://www.whoi.edu/beaufortgyre, last access: October 2019) and Operation IceBridge (OIB; Kurtz et al., 2013) in situ measurements confirm the improvements in the Baseline-D freeboard product quality. Relative to OIB, the Baseline-D freeboard mean bias is reduced by about 8 cm, which roughly corresponds to a 60 % decrease with respect to Baseline-C. The BGEP data indicate a similar tendency with a mean draft bias lowered from 0.85 to −0.14 m. For the two in situ datasets, the root mean square deviation (RMSD) is also well reduced from 14 to 11 cm for OIB and by a factor of 2 for the BGEP. Observations over inland waters show a slight increase in the percentage of good observations in Baseline-D, generally around 5 %–10 % for most lakes. This paper provides an overview of the new Level 1 and Level 2 (L2) CryoSat Ice Baseline-D evolutions and related data quality assessment, based on results obtained from analyzing the 6-month Baseline-D test dataset released to CryoSat expert users prior to the final transfer to operations.

We developed a multi-mission satellite altimetry analysis over the Antarctic Ice Sheet which comprises Seasat, Geosat, ERS-1, ERS-2, Envisat, ICESat and CryoSat-2. After a consistent reprocessing and a stepwise calibration of the inter-mission offsets, we obtained monthly grids of multi-mission surface elevation change (SEC) with respect to the reference epoch 09/2010 (in the format of month/year) from 1978 to 2017. A validation with independent elevation changes from in situ and airborne observations as well as a comparison with a firn model proves that the different missions and observation modes have been successfully combined to a seamless multi-mission time series. For coastal East Antarctica, even Seasat and Geosat provide reliable information and, hence, allow for the analysis of four decades of elevation changes. The spatial and temporal resolution of our result allows for the identification of when and where significant changes in elevation occurred. These time series add detailed information to the evolution of surface elevation in such key regions as Pine Island Glacier, Totten Glacier, Dronning Maud Land or Lake Vostok. After applying a density mask, we calculated time series of mass changes and found that the Antarctic Ice Sheet north of 81.5∘ S was losing mass at an average rate of -85±16 Gt yr−1 between 1992 and 2017, which accelerated to -137±25 Gt yr−1 after 2010.

A key objective in palaeo-climatology is the retrieval of a continuous Antarctic ice-core record dating back 1.5 Ma. The identification of a suitable Antarctic site requires sufficient knowledge of the subglacial landscape beneath the Antarctic Ice Sheet. Here, we present new ice thickness information from the Dome Fuji region, East Antarctica, based on airborne radar surveys conducted during the 2014/15 and 2016/17 southern summers. Compared to previous maps of the region, the new dataset shows a more complex landscape with networks of valleys and mountain plateaus. We use the new dataset as input in a thermokinematic model that incorporates uncertainties in geothermal heat flux values in order to improve the predictions of potential ice-core sites. Our results show that especially the region immediately south of Dome Fuji station persists as a good candidate site for obtaining an old ice core. An initial assessment of basal conditions revealed the existence of what appears to be subglacial lakes. Further radar data analysis shows overall high continuity of layer stratigraphy in the region. This indicates that extending the age–depth information from the Dome Fuji ice core to a new ice-core drill site is a viable option.

Ice shelves control the dynamic mass loss of ice sheets through buttressing and their integrity depends on the spatial variability of their basal mass balance (BMB), i.e. the difference between refreezing and melting. Here, we present an improved technique – based on satellite observations – to capture the small-scale variability in the BMB of ice shelves. As a case study, we apply the methodology to the Roi Baudouin Ice Shelf, Dronning Maud Land, East Antarctica, and derive its yearly averaged BMB at 10 m horizontal gridding. We use mass conservation in a Lagrangian framework based on high-resolution surface velocities, atmospheric-model surface mass balance and hydrostatic ice-thickness fields (derived from TanDEM-X surface elevation). Spatial derivatives are implemented using the total-variation differentiation, which preserves abrupt changes in flow velocities and their spatial gradients. Such changes may reflect a dynamic response to localized basal melting and should be included in the mass budget. Our BMB field exhibits much spatial detail and ranges from −14.7 to 8.6 m a−1 ice equivalent. Highest melt rates are found close to the grounding line where the pressure melting point is high, and the ice shelf slope is steep. The BMB field agrees well with on-site measurements from phase-sensitive radar, although independent radar profiling indicates unresolved spatial variations in firn density. We show that an elliptical surface depression (10 m deep and with an extent of 0.7 km × 1.3 km) lowers by 0.5 to 1.4 m a−1, which we tentatively attribute to a transient adaptation to hydrostatic equilibrium. We find evidence for elevated melting beneath ice shelf channels (with melting being concentrated on the channel's flanks). However, farther downstream from the grounding line, the majority of ice shelf channels advect passively (i.e. no melting nor refreezing) toward the ice shelf front. Although the absolute, satellite-based BMB values remain uncertain, we have high confidence in the spatial variability on sub-kilometre scales. This study highlights expected challenges for a full coupling between ice and ocean models.

The formation mechanisms of ice-shelf channels remain poorly understood. Here, using ice-penetrating radar data, the authors propose that ice-shelf channel morphology in the Roi Baudouin Ice Shelf, East Antarctica, is seeded by esker ridges, indenting the ice from below.

We report annual snow accumulation rates from 1959 to 2004 along a 250 km segment of the Expéditions Glaciologiques Internationales au Groenland (EGIG) line across central Greenland using Airborne SAR/Interferometric Radar Altimeter System (ASIRAS) radar layers and high resolution neutron-probe (NP) density profiles. ASIRAS-NP-derived accumulation rates are not statistically different (95 % confidence interval) from in situ EGIG accumulation measurements from 1985 to 2004. ASIRAS-NP-derived accumulation increases by 20 % below 3000 m elevation, and increases by 13 % above 3000 m elevation for the period 1995 to 2004 compared to 1985 to 1994. Three Regional Climate Models (PolarMM5, RACMO2.3, MAR) underestimate snow accumulation below 3000 m by 16–20 % compared to ASIRAS-NP from 1985 to 2004. We test radar-derived accumulation rates sensitivity to density using modeled density profiles in place of NP densities. ASIRAS radar layers combined with Herron and Langway (1980) model density profiles (ASIRAS-HL) produce accumulation rates within 3.5 % of ASIRAS-NP estimates in the dry snow region. We suggest using Herron and Langway (1980) density profiles to calibrate radar layers detected in dry snow regions of ice sheets lacking detailed in situ density measurements, such as those observed by the Operation IceBridge campaign.

Knowledge about Antarctic sea-ice volume and its changes over the past decades has been sparse due to the lack of systematic sea-ice thickness measurements in this remote area. Recently, first attempts have been made to develop a sea-ice thickness product over the Southern Ocean from space-borne radar altimetry and results look promising. Today, more than 20 years of radar altimeter data are potentially available for such products. However, the characteristics of individual radar types differ for the available altimeter missions. Hence, it is important and our goal to study the consistency between single sensors in order to develop long and consistent time series. Here, the consistency between freeboard measurements of the Radar Altimeter 2 on board Envisat and freeboard measurements from the Synthetic-Aperture Interferometric Radar Altimeter on board CryoSat-2 is tested for their overlap period in 2011. Results indicate that mean and modal values are in reasonable agreement over the sea-ice growth season (May–October) and partly also beyond. In general, Envisat data show higher freeboards in the first-year ice zone while CryoSat-2 freeboards are higher in the multiyear ice zone and near the coasts. This has consequences for the agreement in individual sectors of the Southern Ocean, where one or the other ice class may dominate. Nevertheless, over the growth season, mean freeboard for the entire (regionally separated) Southern Ocean differs generally by not more than 3 cm (8 cm, with few exceptions) between Envisat and CryoSat-2, and the differences between modal freeboards lie generally within ±10 cm and often even below.

In the context of quantifying Arctic ice-volume decrease at global scale, the CryoSat-2 satellite was launched in 2010 and is equipped with the Ku band synthetic aperture radar altimeter SIRAL (Synthetic Aperture Interferometric Radar Altimeter), which we use to derive sea-ice freeboard defined as the height of the ice surface above the sea level. Accurate CryoSat-2 range measurements over open water and the ice surface of the order of centimetres are necessary to achieve the required accuracy of the freeboard-to-thickness conversion. Besides uncertainties of the actual sea-surface height and limited knowledge of ice and snow properties, the composition of radar backscatter and therefore the interpretation of radar echoes is crucial. This has consequences in the selection of retracker algorithms which are used to track the main scattering horizon and assign a range estimate to each CryoSat-2 measurement. In this study we apply a retracker algorithm with thresholds of 40, 50 and 80% of the first maximum of radar echo power, spanning the range of values used in the current literature. By using the selected retrackers and additionally results from airborne validation measurements, we evaluate the uncertainties of sea-ice freeboard and higher-level products that arise from the choice of the retracker threshold only, independent of the uncertainties related to snow and ice properties. Our study shows that the choice of retracker thresholds does have a significant impact on magnitudes of estimates of sea-ice freeboard and thickness, but that the spatial distributions of these parameters are less affected. Specifically we find mean radar freeboard values of 0.121 m (0.265 m) for the 40% threshold, 0.086 m (0.203 m) for the 50% threshold and 0.024 m (0.092 m) for the 80% threshold, considering first-year ice (multiyear ice) in March 2013. We show that the main source of freeboard and thickness uncertainty results from the choice of the retracker and the unknown penetration of the radar pulse into the snow layer in conjunction with surface roughness effects. These uncertainties can cause a freeboard bias of roughly 0.06–0.12 m. Furthermore we obtain a significant rise of 0.02–0.15 m of freeboard from March 2013 to November 2013 in the area for multiyear sea ice north of Greenland and Canada. Since this is unlikely, it gives rise to the assumption that applying different retracker thresholds depending on seasonal properties of the snow load is necessary in the future.

This study focuses on the present-day surface elevation of the Greenland and Antarctic ice sheets. Based on 3 years of CryoSat-2 data acquisition we derived new elevation models (DEMs) as well as elevation change maps and volume change estimates for both ice sheets. Here we present the new DEMs and their corresponding error maps. The accuracy of the derived DEMs for Greenland and Antarctica is similar to those of previous DEMs obtained by satellite-based laser and radar altimeters. Comparisons with ICESat data show that 80% of the CryoSat-2 DEMs have an uncertainty of less than 3 m ± 15 m. The surface elevation change rates between January 2011 and January 2014 are presented for both ice sheets. We compared our results to elevation change rates obtained from ICESat data covering the time period from 2003 to 2009. The comparison reveals that in West Antarctica the volume loss has increased by a factor of 3. It also shows an anomalous thickening in Dronning Maud Land, East Antarctica which represents a known large-scale accumulation event. This anomaly partly compensates for the observed increased volume loss of the Antarctic Peninsula and West Antarctica. For Greenland we find a volume loss increased by a factor of 2.5 compared to the ICESat period with large negative elevation changes concentrated at the west and southeast coasts. The combined volume change of Greenland and Antarctica for the observation period is estimated to be −503 ± 107 km3 yr−1. Greenland contributes nearly 75% to the total volume change with −375 ± 24 km3 yr−1.

We have derived digital elevation models (DEMs) over the western part of the Devon Ice Cap in Nunavut, Canada, using "swath processing" of interferometric data collected by Cryosat between February 2011 and January 2012. With the standard ESA (European Space Agency) SARIn (synthetic aperture radar interferometry) level 2 (L2) data product, the interferometric mode is used to map the cross-track position and elevation of the "point-of-closest-approach" (POCA) in sloping glacial terrain. However, in this work we explore the extent to which the phase of the returns in the intermediate L1b product can also be used to map the heights of time-delayed footprints beyond the POCA. We show that there is a range of average cross-track slopes (~ 0.5 to ~ 2°) for which the returns will be dominated by those beneath the satellite in the main beam of the antenna so that the resulting interferometric phase allows mapping of heights in the delayed range window beyond the POCA. In this way a swath of elevation data is mapped, allowing the creation of DEMs from a sequence of L1b SARIn Cryosat data takes. Comparison of the Devon results with airborne scanning laser data showed a mean difference of order 1 m with a standard deviation of about 1 m. The limitations of swath processing, which generates almost 2 orders of magnitude more data than traditional radar altimetry, are explored through simulation, and the strengths and weaknesses of the technique are discussed.

Sea ice thickness is one of the most sensitive variables in the Arctic climate system. In order to quantify changes in sea ice thickness, CryoSat-2 was launched in 2010 carrying a Ku-band radar altimeter (SIRAL) designed to measure sea ice freeboard with a few centimeters accuracy. The instrument uses the synthetic aperture radar technique providing signals with a resolution of about 300 m along track. In this study, airborne Ku-band radar altimeter data over different sea ice types have been analyzed. A set of parameters has been defined to characterize the differences in strength and width of the returned power waveforms. With a Bayesian-based method, it is possible to classify about 80% of the waveforms from three parameters: maximum of the returned power waveform, the trailing edge width and pulse peakiness. Furthermore, the maximum of the power waveform can be used to reduce the number of false detections of leads, compared to the widely used pulse peakiness parameter. For the pulse peakiness the false classification rate is 12.6% while for the power maximum it is reduced to 6.5%. The ability to distinguish between different ice types and leads allows us to improve the freeboard retrieval and the conversion from freeboard into sea ice thickness, where surface type dependent values for the sea ice density and snow load can be used.

We present estimates of accumulation rate along a 200 km transect ranging in elevation from 2750 to 3150 m in the dry snow zone on the western slope of the Greenland Ice Sheet. An airborne radar altimeter is used to estimate the thickness of annual internal layers and, in conjunction with ground based snow/firn density profiles, annual accumulation rates between 1998 and 2003 are derived. A clear gradient in the thickness of each layer observed by the radar altimeter and in the associated estimates of annual accumulation is seen along the transect, with a 33.6% ± 16% mean decrease in accumulation from west to east. The observed inter-annual variability is high, with the annual mean accumulation rate estimated at 0.359 m.w.e. yr−1 (s.d. ± 0.049 m.w.e. yr−1). Mean accumulation rates modelled using meteorological models overestimate our results by 16% on average, but by 32% and 42% in the years 2001 and 2002. The methodology presented here demonstrates the potential to obtain accurate and spatially extensive accumulation rates from radar altimeters in regions of ice sheets where field observations are sparse, and accumulation rates greater than several tens of cm.

Abstract In this work we have assessed the decolorization of textile effluents throughout their treatment in a solid-state fermentation (SSF) system. SSF assays were conducted with peach-palm (Bactris gasipaes) residue using the white rot fungus Ganoderma lucidum EF 31. The influence of the dye concentration and of the amounts of peach-palm residue and liquid phase on both the discoloration efficiency and enzyme production was studied. According to our results, independently of experimental conditions employed, laccase was the main ligninolytic enzyme produced by G. lucidum. The highest laccase activity was obtained at very low effluent concentrations, suggesting the existence of an inhibitory effect of higher concentrations on fungal metabolism. The highest percentage of color removal was reached when 10 grams of peach palm residue was moistened with 60 mL of the final effluent. In control tests carried out with the synthetic dye Remazol Brilliant Blue R (RBBR) decolorization efficiencies about 20% higher than that achieved with the industrial effluent were achieved. The adsorption of RBBR on peach-palm residue was also investigated. Equilibrium tests showed that the adsorption of this dye followed both Langmuir and Freundlich isotherms. Hence, our experimental results indicate that peach-palm residue is suitable substrate for both laccase production and color removal in industrial effluents.

The production of ethanol from lignocellulosic biomass is referred as a second generation biofuel, whose processing is one of the most promising technologies under development. There are few available studies on the use of enzymes produced by fungi as active for the biodegradation of lignocellulosic biomass. However, the manganese peroxidase (MnP) enzyme presents high potential to degrade lignin and the basidiomycetes are the major producers of this oxidase. Thus, this study aimed at evaluating the ability of fungi Lentinula edodes and Lentinula boryana to produce this enzyme when cultivated in submerged fermentation system (SS) and also in solid-state fermentation system (SSF) containing Eucalyptus benthamii sawdust with or without corn cob meal. In the SS the greatest MnP expression occurred on the 25th day, being of 70 UI.L–1 for L. boryana and of 20 UI.L–1 for L. edodes. In the SSF, the best results were obtained on the 10th day for L. edodes, while for L. boryana it happened between the 20th and the 25th days, despite both species presented values close to 110 UI.L–1. Therefore, the results indicated that the studied fungi express the enzyme of interest and that its production is enhanced when cultivated in solid system.

Abstract Lignocellulose is the most abundant environmental component and a renewable organic resource in soil. There are some filamentous fungi which developed the ability to break down and use cellulose, hemicellulose and lignin as an energy source. The objective of this research was to analyze the effect of three nitrogen resources (ammonium sulfate, saltpetre, soybean) in the holocellulolitic activity of Lentinula edodes EF 50 using as substrate sawdust E. benthamii. An experimental design mixture was applied with repetition in the central point consisting of seven treatments (T) of equal concentrations of nitrogen in ammonium sulfate, potassium nitrate and soybean. The enzymatic activity of avicelase, carboxymetilcellulase, β-glucosidase, xylanases and manganese peroxidase was determined. The humidity, pH, water activity (aw) and qualitative analysis of mycelial growth in 8 times of cultivation were evaluated. The results showed negative effect on enzyme production in treatments with maximum concentration of ammonium sulfate and potassium nitrate. The treatments with cooked soybean flour expressed higher enzymatic activities in times of 3, 6 and 9 days of culture, except in the activity of manganese peroxidase. The highest production was observed in the treatment with ammonium sulfate, and soybean (83.86 UI.L–1) at 20 days of cultivation.

We present a dataset of reconstructed three-dimensional (3D) englacial stratigraphic horizons in northern Greenland. The data cover four different regions representing key ice-dynamic settings in Greenland: (i) the onset of Petermann Glacier, (ii) a region upstream of the 79° North Glacier (Nioghalvfjerdsbræ), near the northern Greenland ice divide, (iii) the onset of the Northeast Greenland Ice Stream (NEGIS) and (iv) a 700 km wide region extending across the central ice divide over the entire northern part of central Greenland. In this paper, we promote the advantages of a 3D perspective of deformed englacial stratigraphy and explain how 3D horizons provide an improved basis for interpreting and reconstructing the ice-dynamic history. The 3D horizons are provided in various formats to allow a wide range of applications and reproducibility of results., (© 2023. Springer Nature Limited.)

The dynamic mass loss of ice sheets constitutes one of the biggest uncertainties in projections of ice-sheet evolution. One central, understudied aspect of ice flow is how the bulk orientation of the crystal orientation fabric translates to the mechanical anisotropy of ice. Here we show the spatial distribution of the depth-averaged horizontal anisotropy and corresponding directional flow-enhancement factors covering a large area of the Northeast Greenland Ice Stream onset. Our results are based on airborne and ground-based radar surveys, ice-core observations, and numerical ice-flow modelling. They show a strong spatial variability of the horizontal anisotropy and a rapid crystal reorganisation on the order of hundreds of years coinciding with the ice-stream geometry. Compared to isotropic ice, parts of the ice stream are found to be more than one order of magnitude harder for along-flow extension/compression while the shear margins are potentially softened by a factor of two for horizontal-shear deformation., (© 2023. The Author(s).)

Over the past two decades, ice loss from the Greenland ice sheet (GrIS) has increased owing to enhanced surface melting and ice discharge to the ocean 1-5 . Whether continuing increased ice loss will accelerate further, and by how much, remains contentious 6-9 . A main contributor to future ice loss is the Northeast Greenland Ice Stream (NEGIS), Greenland's largest basin and a prominent feature of fast-flowing ice that reaches the interior of the GrIS 10-12 . Owing to its topographic setting, this sector is vulnerable to rapid retreat, leading to unstable conditions similar to those in the marine-based setting of ice streams in Antarctica 13-20 . Here we show that extensive speed-up and thinning triggered by frontal changes in 2012 have already propagated more than 200 km inland. We use unique global navigation satellite system (GNSS) observations, combined with surface elevation changes and surface speeds obtained from satellite data, to select the correct basal conditions to be used in ice flow numerical models, which we then use for future simulations. Our model results indicate that this marine-based sector alone will contribute 13.5-15.5 mm sea-level rise by 2100 (equivalent to the contribution of the entire ice sheet over the past 50 years) and will cause precipitous changes in the coming century. This study shows that measurements of subtle changes in the ice speed and elevation inland help to constrain numerical models of the future mass balance and higher-end projections show better agreement with observations., (© 2022. The Author(s).)

We use satellite and airborne altimetry to estimate annual mass changes of the Greenland Ice Sheet. We estimate ice loss corresponding to a sea-level rise of 6.9 ± 0.4 mm from April 2011 to April 2020, with a highest annual ice loss rate of 1.4 mm/yr sea-level equivalent from April 2019 to April 2020. On a regional scale, our annual mass loss timeseries reveals 10-15 m/yr dynamic thickening at the terminus of Jakobshavn Isbræ from April 2016 to April 2018, followed by a return to dynamic thinning. We observe contrasting patterns of mass loss acceleration in different basins across the ice sheet and suggest that these spatiotemporal trends could be useful for calibrating and validating prognostic ice sheet models. In addition to resolving the spatial and temporal fingerprint of Greenland's recent ice loss, these mass loss grids are key for partitioning contemporary elastic vertical land motion from longer-term glacial isostatic adjustment (GIA) trends at GPS stations around the ice sheet. Our ice-loss product results in a significantly different GIA interpretation from a previous ice-loss product., (© 2022 The Authors.)

The Greenland Ice Sheet is the largest land ice contributor to sea level rise. This will continue in the future but at an uncertain rate and observational estimates are limited to the last few decades. Understanding the long-term glacier response to external forcing is key to improving projections. Here we use historical photographs to calculate ice loss from 1880-2012 for Jakobshavn, Helheim, and Kangerlussuaq glacier. We estimate ice loss corresponding to a sea level rise of 8.1 ± 1.1 millimetres from these three glaciers. Projections of mass loss for these glaciers, using the worst-case scenario, Representative Concentration Pathways 8.5, suggest a sea level contribution of 9.1-14.9 mm by 2100. RCP8.5 implies an additional global temperature increase of 3.7 °C by 2100, approximately four times larger than that which has taken place since 1880. We infer that projections forced by RCP8.5 underestimate glacier mass loss which could exceed this worst-case scenario.

We report the discovery of a large impact crater beneath Hiawatha Glacier in northwest Greenland. From airborne radar surveys, we identify a 31-kilometer-wide, circular bedrock depression beneath up to a kilometer of ice. This depression has an elevated rim that cross-cuts tributary subglacial channels and a subdued central uplift that appears to be actively eroding. From ground investigations of the deglaciated foreland, we identify overprinted structures within Precambrian bedrock along the ice margin that strike tangent to the subglacial rim. Glaciofluvial sediment from the largest river draining the crater contains shocked quartz and other impact-related grains. Geochemical analysis of this sediment indicates that the impactor was a fractionated iron asteroid, which must have been more than a kilometer wide to produce the identified crater. Radiostratigraphy of the ice in the crater shows that the Holocene ice is continuous and conformable, but all deeper and older ice appears to be debris rich or heavily disturbed. The age of this impact crater is presently unknown, but from our geological and geophysical evidence, we conclude that it is unlikely to predate the Pleistocene inception of the Greenland Ice Sheet.

Ice-shelf channels are long curvilinear tracts of thin ice found on Antarctic ice shelves. Many of them originate near the grounding line, but their formation mechanisms remain poorly understood. Here we use ice-penetrating radar data from Roi Baudouin Ice Shelf, East Antarctica, to infer that the morphology of several ice-shelf channels is seeded upstream of the grounding line by large basal obstacles indenting the ice from below. We interpret each obstacle as an esker ridge formed from sediments deposited by subglacial water conduits, and calculate that the eskers' size grows towards the grounding line where deposition rates are maximum. Relict features on the shelf indicate that these linked systems of subglacial conduits and ice-shelf channels have been changing over the past few centuries. Because ice-shelf channels are loci where intense melting occurs to thin an ice shelf, these findings expose a novel link between subglacial drainage, sedimentation and ice-shelf stability.