Search

The geometry of englacial isochrones is a product of the past and present ice velocity field and is useful for our understanding of steady-state ice flow dynamics, flow regime re-organisation, and calibration of models. Ice rises contain various flow regimes (divide flow, flank flow, and grounding zones) on small spatial scales, meaning they are ideal locations to study ice-flow dynamics and stratigraphy to constrain model parameters. We run full Stokes, thermo-mechanically coupled simulations of Derwael Ice Rise in East Antarctica and simulate the three-dimensional stratigraphy of the ice rise and the surrounding ice shelf using the finite element model Elmer/Ice. Over the ice rise, we derive the accumulation rate from internal reflection horizons and use RACMO2.3 surface mass balance data over the surrounding ice shelf. Simulations are run for Glen's flow law exponents of n=3 and n=4 with appropriate values derived for the Arrhenius law. To calibrate the model, comparisons are made with the BedMachine surface elevation and density-adjusted internal reflection horizons observed in many transects recorded by AWI’s ultra-wide band radar covering the divide, the flanks, and the grounding zones. To understand ice flow dynamics where the velocity field of the ice rise and the ice shelf converge in the compressive and shear zones, we analyse the modelled englacial stress and strain rate fields. Our results allow us to investigate isochronal structures where observed internal reflection horizons are too steep or obscured to be adequately picked up by radar. A comparison between the model and observed fracturing can be used to infer threshold stress and strain rates for fracture initiation. These simulations are a blueprint for the full Stokes, three-dimensional modelling of ice rises and have further relevance in the study of three-dimensional influences on Raymond arch evolution, the constrained coupling of the anisotropy equations, comparisons with ice core data and the automated inference of ice flow parameters from internal reflection horizons.

The internal, isochronous layering of glaciers is shaped by accumulation and ice deformation. Information about these processes can be inferred from observing the layers using radar sounding. The reflectivity of the layers depends on density (permittivity) and acidity (conductivity) contrasts which tend to decrease with depth. At places like alpine glaciers where logistic limitations often only allow the deployment of lightweight and power-constrained ground-penetrating radar systems, it can therefore be challenging to illuminate the deeper radio-stratigraphy.The phase-sensitive Radio Echo Sounder (pRES) is a lightweight frequency modulated continuous wave radar which allows the use of coherent Synthetic Aperture Radar (SAR) processing techniques to improve the signal-to-noise ratio of internal reflection horizons. Using a mobile pRES we collected a radar profile on an alpine glacier (Colle Gnifetti, Italy/Switzerland). Here, we demonstrate how to apply layer-optimized SAR techniques to make deep internal layers visible, which could not be seen by a conventional pulsed radar. We evaluate the requirements on spatial resolution and positioning accuracy during data acquisition, necessary for applying layer-optimized SAR processing, as they constrain the feasibility of the method. We further discuss implications on how density and acidity contribute to decreasing dielectric contrasts.

Promontory ice rises are locally grounded features adjacent to ice shelves that are still connected to the ice sheet. Ice rises are an archive for the atmospheric and ice dynamic history of the respective outflow regions where the presence, absence, or migration of Raymond arches in radar stratigraphy represents a memory of the ice-rise evolution. However, ice rises and their inferred dynamic history are not yet used to constrain large-scale ice flow model spin-ups because matching the arch amplitudes includes many unknown parameters, e.g., those pertaining to ice rheology. In particular, anisotropic ice flow models predict gradients in ice fabric anisotropy on either side of an ice divide. However, this has thus far not been validated with observations. The ground-based phase-sensitive Radio Echo Sounder (pRES) has previously been used to infer ice fabric types for various flow regimes using the co-polarized polarimetric coherence phase as a metric to extract information from the birefringent radar backscatter. Here, we apply this technique using quad-polarimetric radar data along a 5 km transect across a ridge near the triple junction of Hammarryggen Ice Rise at the Princess Ragnhild Coast. A comparison with ice core data collected at the dome shows that the magnitude of ice fabric anisotropy can reliably be reconstructed from the quad-polarimetric data. We use the combined dataset also to infer the spatial variation of ice fabric orientations in the vicinity of the triple junction. The observations are integrated with airborne radar profiles and strain rates based on the shallow ice approximation. We then discuss whether estimated anisotropy from radar polarimetry on ice rises, in general, can be another observational constraint to better ice rises as an archive of ice dynamics.

Ice shelves buttress ice flow from Antarctica’s interior and provide the interface for ice-ocean interactions. Here, we present a comprehensive dataset collected with autonomous phase-sensitive radio-echo sounders (ApRES) on the Ekström Ice Shelf in East Antarctica. The data include a one year time series of basal melt near the grounding zone and > 40 repeat and quad-polarimetric observations covering the entire ice shelf with the transition to grounded ice and multiple lateral shear zones. The inferred melt rates are put into context in terms of their seasonality and with respect to spatial patterns which were mapped in the internal ice stratigraphy. The polarimetric backscatter show signs of anisotropic ice fabric and spatial changes can be traced coherently from the grounding line to the ice-shelf front. We investigate those signatures in conjunction with the vertical strain rates inferred from ApRES and ice-flow modelling to learn more about ice-shelf rheology, particularly with respect to the stress exponent n.

Ice shelves surrounding the Antarctic perimeter moderate ice discharge towards the ocean through buttressing. Ice-shelf evolution and integrity depend on the local surface accumulation, basal melting and on the spatially variable ice-shelf viscosity. These components of ice-shelf mass balance are often poorly constrained by observations and introduce uncertainties in ice-sheet projections. Isochronal radar stratigraphy is an observational archive for the atmospheric, oceanographic and ice-flow history of ice shelves. Here, we predict the stratigraphy of locally accumulated ice on ice shelves with a kinematic forward model for a given atmospheric and oceanographic scenario. This delineates the boundary between local meteoric ice (LMI) and continental meteoric ice (CMI). A large LMI to CMI ratio hereby marks ice shelves whose buttressing strength is more sensitive to changes in atmospheric precipitation patterns. A mismatch between the steady-state predictions of the kinematic forward model and observations from radar can highlight inconsistencies in the atmospheric and oceanographic input data or be an indicator for a transient ice-shelf history not accounted for in the model. We discuss pitfalls in numerical diffusion when calculating the age field and validate the kinematic model with the full Stokes ice-flow model Elmer/Ice. The Roi Baudouin Ice Shelf (East Antarctica) serves as a test case for this approach. There, we find a significant east–west gradient in the LMI / CMI ratio. The steady-state predictions concur with observations on larger spatial scales (>10 km), but deviations on smaller scales are significant, e.g., because local surface accumulation patterns near the grounding zone are underestimated in Antarctic-wide estimates. Future studies can use these mismatches to optimize the input data or to pinpoint transient signatures in the ice-shelf history using the ever growing archive of radar observations of internal ice stratigraphy.

We collected surface snow samples from three different glaciers – Yala, Thana, and Sachin – in the central and western Himalayas to understand the spatial variability and radiative impacts of light-absorbing particles. The Yala and Thana glaciers in Nepal and Bhutan, respectively, were selected to represent the central Himalayas. The Sachin glacier in Pakistan was selected to represent the western Himalayas. The samples were collected during the pre- and post-monsoon seasons of the year 2016. The samples were analyzed for black carbon (BC) and water-insoluble organic carbon (OC) through the thermal optical method. The average mass concentrations (BC 2381 ng g−1; OC 3896 ng g−1; dust 101 µg g−1) in the western Himalayas (Sachin glacier) were quite high compared to the mass concentrations (BC 358 ng g−1, OC 904 ng g−1, dust 22 µg g−1) in the central Himalayas (Yala glacier). The difference in mass concentration may be due to the difference in elevation, snow age, local pollution sources, and meteorological conditions. BC in surface snow was also estimated through Weather Research and Forecasting (WRF) model coupled with Chemistry (WRF-Chem) simulations at the three glacier sites during the sampling periods. Simulations reasonably capture the spatial and seasonal patterns of the observed BC in snow but with a relatively smaller magnitude. Absolute snow albedo was estimated through the Snow, Ice, and Aerosol Radiative (SNICAR) model. The absolute snow albedo reduction ranged from 0.48 % (Thana glacier during September) to 24 % (Sachin glacier during May) due to BC and 0.13 % (Yala glacier during September) to 5 % (Sachin glacier during May) due to dust. The instantaneous radiative forcing due to BC and dust was estimated in the range of 0 to 96.48 and 0 to 25 W m−2, respectively. The lowest and highest albedo reduction and radiative forcing were observed in central and western Himalayan glaciers, respectively. The potential source regions of the deposited pollutants were inferred using WRF-Chem tagged-tracer simulations. Selected glaciers in the western Himalayas were mostly affected by long-range transport from the Middle East and central Asia; however, the central Himalayan glaciers were mainly affected by local and south Asia emissions (from Nepal, India, and China) especially during the pre-monsoon season. Overall, south Asia and west Asia were the main contributing source regions of pollutants.

The internal ice stratigraphy as imaged by radar is an integrated archive of the atmospheric- oceanographic, and ice-dynamic history that the ice sheet has experienced. It provides an observational constraint for ice flow modeling that has been used for instance to predict age-depth relationships at prospective ice-coring sites in Antarctica’s interior. The stratigraphy is typically more disturbed and more difficult to image in coastal regions due to faster ice flow. Yet, knowledge of ice stratigraphy across ice shelf grounding lines and further seawards is important to help constrain ocean-induced melting and associated stability.Here, we present preliminary results of synthesizing information from radar stratigraphic characteristics from airborne and ground-based radar surveys that have been collected for specific projects starting from the 1990s onwards focusing on ice marginal zones of Antarctica. The key data is based on airborne surveys from the German Alfred Wegener Institute’s polar aircrafts equipped with a 150 MHz radar. In the meantime this system has been replaced by an ultra-wide band 150-520 MHz radar. The older data will provide a baseline with extensive coverage that can be used for model calibration and change detection over time. We aim to provide metrics of the radio stratigraphy (e.g. shape and slope of internal reflection horizons) as well as classified prevalent stratigraphy types that can be used to calibrate machine learning approaches such as simulation based inference. The data obtained will be integrated in coordination efforts within the SCAR AntArchitecture Action Group.

The internal stratigraphy of alpine glaciers entails information about its past dynamics and accumulation rates. It further can be used for intercalibrating the age-depth scales of ice cores. The internal ice stratigraphy is often imaged using radar, but similar to polar ice sheets the deeper stratigraphy is often difficult to resolve with classical pulsed radar systems. For polar ice sheets, the introduction of phase coherent radars has illuminated this former echo-free zone (EFZ) and now patterns of folded, buckled and disrupted ice stratigraphy are clearly visible. Unfortunately, the new airborne and ground-based radar systems applied in polar regions are typically too heavy to be deployed in an alpine environment.Here, we transfer the lightweight autonomous phase-sensitive radio-echo sounder (ApRES) to an alpine glacier targeting its echo-free zone (Colle Gnifetti, Italy/Switzerland). The ApRES is a coherent frequency modulated continuous wave radar with an integration time of 1 s per trace which we deployed in combination with a GNSS used in real time kinematic (RTK) mode. The latter allows repositioning of the antennas with sub-wavelength accuracy (approximately 5 cm) required to exploit the coherent signal. Like this, the radio-stratigraphy of the former EFZ at this site could be imaged using a matched filtering SAR method. The resulting radargrams cover former ice core sites (e.g., Ice Memory and KCC) and can be used to harmonize conflicting age-depth scales. This dataset will be analysed further in conjunction with ice-fabric measurements from ice cores to reveal how the anisotropic ice rheology imprints on the flow field of glaciers.

Mapping the ice bed interface with radar is challenging in many alpine glaciers where the ice is temperate, and in-ice absorption is high. It is also difficult in selected regions of polar ice sheets such as near grounding zones and in ice streams where clutter and rough beds increase incoherent volume scattering. The lack of information for the ice geometry impedes our process understanding, e.g., basal sliding (requires knowledge about the basal roughness) and the routing of subglacial water flow (requires knowledge on basal smoothness). The lack of observations to constrain variations in ice thickness on the sub-kilometre scale is thus still a bottleneck to confidently predict ice dynamics and expected rates of sea-level rise.A recent development in radioglaciology, namely the application of phase-coherent polarimetric radar, provides an excellent opportunity to overcome these limitations. Radar polarimetry has made significant strides in the last few years to constrain internal ice structure and their impact on the deformation of ice sheets, including the reconstruction of ice micro-structure parameters previously obtained from ice cores. Here, we suggest that the ice-bed interface can be identified in characteristic patterns of the polarimetric coherence phase. This new metric provides information in areas where the backscattered power amplitude does not show any signatures of the ice-bed interface. We provide examples for this across a wide range of glaciological settings, including cold (Colle Gnifetti, Switzerland) and temperate (Hintereisferner, Austria) alpine glaciers, thin grounding zones (Ekström Ice Shelf, East Antarctica) and thick ice domes (Dome C, East Antarctica). If this holds, then the ice thickness mapping in challenging glaciological settings should preferably be done using a quad-polarimetric acquisition geometry. For ground-based surveys, this can be done using an autonomous ice rover, for which we provide a proof-of-concept study on the Ekström Ice Shelf in Antarctica.

Ice shelves surrounding the Antarctic perimeter buttress ice flow from the continent towards the ocean, and their disintegration leads to an increase in ice discharge and sea level rise. The evolution and integrity of ice shelves is governed by surface accumulation, basal melting, and ice dynamics. We find history of these processes imprinted in the ice-shelf stratigraphy, which is mapped using isochrones imaged with radar. As an observational archive, the radar obtained stratigraphy combined with ice flow modeling has high potential to assist model calibration and reduce uncertainties in projections for the ice-sheet evolution. In this study we use a simplistic and observationally driven ice-dynamic forward model to predict the ice-shelf stratigraphy. We validate this approach with the full Stokes ice-flow model Elmer/Ice, and present a test-case for the Roi Baudouin Ice Shelf (East Antarctica) - where our model predictions agree well with radar obtained observations. The presented method enables us to investigate whether ice shelves are in steady-state, as well as to map spatial variations of how much of the ice-shelf volume is determined by its local surface mass balance. In the case of Roi Baudouin, we find the ice-shelf volume in the western part to be dominated by ice inflowing from the ice sheet, while the eastern part of the ice shelf is dominated by ice locally accumulated on the shelf. Such analysis serves as a metric for the susceptibility of ice shelves to climate change. We further apply our approach to other ice shelves in Antarctica.

Acquisition of quad-polarimetric radar data on ice sheets gives insights about the ice-fabric variability with depth and consequently can deliver essential constraints on the spatially variable ice rheology. Polarimetric measurements are collected manually in most ground-based surveys, discretely sampling a limited profile range. Measurements are time-intensive and often do not cover critical areas such as shear zones where field safety is a concern. Autonomous rovers can provide an alternative that optimizes for time, sampling resolution and safety. Here, we present an autonomous acquisition technique of quad-polarimetric radar data using a rover. This technique is based on a previous layout that has proven its capacity to navigate in various snow conditions but did not yet actively trigger the geophysical instruments attached. We upgraded the rover with a novel Robotic Operating System (ROS2) that interfaces simultaneously with a real-time positioning GPS and an automatic phase-sensitive radio-echo sounder (ApRES) with multiple transmitters multiple receivers. Like this, the rover can autonomously steer to pre-destined waypoints and then take static measurements at those locations also in areas where field safety might be compromised. We demonstrate this proof-of-concept on the Ekström Ice Shelf Antarctica, where we acquired densely spaced polarimetric radar data measurements. The rover’s operating system offers many opportunities for other measurement principles, e.g., densely spaced co-polarized data suitable for synthetic aperture radar (SAR) processing.

Summary Precipitation is critical to the water towers of the Hindu Kush-Himalaya-Karakoram region, exerting an important control on glacier mass balance and the water resources for 1.65 billion people. Given that hydroclimatic extremes and water stress have emerged as key hazards in the context of climate change, Nepal's Khumbu region overlaps key vulnerabilities. Here, we investigate the region's precipitation characteristics and moisture sources through analysis of data from a new high-altitude network of automatic weather stations, which allows for a more complete understanding of the climatological precipitation data that are critical information for local communities in the Khumbu region, visitors, and downstream populations. Our findings demonstrate that the northern Bay of Bengal is potentially an important moisture source during the monsoon period (June to August) and that westerly trajectories over land predominate for precipitation events during the post-monsoon, winter, and pre-monsoon seasons.

As the highest mountain on Earth, Mount Everest is an iconic peak that offers an unrivalled natural platform for measuring ongoing climate change across the full elevation range of Asia’s water towers. However, Everest’s extreme environment challenges data collection, particularly on the mountain’s upper slopes, where glaciers accumulate mass and mountaineers are most exposed. Weather stations have operated on Everest before, including the world’s previous highest, but coverage has been sparse in space and time. Here we describe the installation of a network of five automatic weather stations (AWSs), including the two highest stations on Earth (8,430 and 7,945 m MSL) which greatly improves monitoring of this iconic mountain. We highlight sample applications of the new data, including an initial assessment of surface energy fluxes at Camp II (6,464 m MSL) and the South Col (7,945 m MSL), which suggest melt occurs at both sites, despite persistently below-freezing air temperatures. This analysis indicates that melt may even be possible at the 8,850 m MSL summit, and prompts a reevaluation of empirical temperature index models used to simulate glacier melt in the Himalayas that focus only on air temperature. We also provide the first evaluation of numerical weather forecasts at almost 8,000 m MSL and use of model output statistics to reduce forecast error, showcasing an important opportunity to improve climber safety on Everest. Looking forward, we emphasize the considerable potential of these freely available data for understanding weather and climate in the Himalayas and beyond, including tracking the behavior of upper-atmosphere winds, which the AWS network is uniquely positioned to monitor.

We collected surface snow samples from three different glaciers: Yala, Thana, and Sachin in the central and western Himalayas to understand the spatial variability and radiative impacts of light-absorbing particles. The Yala and Thana glaciers in Nepal and Bhutan, respectively, were selected to represent the central Himalayas. The Sachin glacier in Pakistan was selected to represent the western Himalayas. The samples were collected during the pre-and post-monsoon seasons of the year 2016. The samples were analysed for black carbon (BC) and water-insoluble organic carbon (OC) through the thermal optical method. The average mass concentrations (BC 2381.39 ng g−1; OC 3896.00 ng g−1; dust 101.05 µg g−1) in the western Himalaya (Sachin glacier) were quite higher compared to the mass concentrations (BC 357.93 ng g−1, OC 903.86 ng g−1, dust 21.95 µg g−1) at the central Himalaya (Yala glacier). The difference in mass concentration may be due to the difference in elevation, snow age, local pollution sources, and difference in meteorological conditions. BC in surface snow was also estimated through WRF-Chem simulations at the three glacier sites during the sampling periods. Simulations reasonably capture the spatial and seasonal patterns of the observed BC in snow but with a relatively smaller magnitude. Absolute snow albedo was estimated through the Snow, Ice, and Aerosol Radiation (SNICAR) model. The absolute snow albedo reduction was ranging between 0.48 % (Thana glacier during September) to 24 % (Sachin glacier during May) due to BC and 0.13 % (Yala glacier during September) to 5 % (Sachin glacier during May) due to dust. The instantaneous radiative forcing due to BC and dust was estimated in the range of 0 to 96.48 W m−2 and 0 to 25 W m−2 respectively. The lowest and highest albedo reduction and radiative forcing were observed in central and western Himalayan glaciers, respectively. The potential source regions of the deposited pollutants were inferred using WRF-Chem tagged-tracer simulations. Selected glaciers in the western Himalayas were mostly affected by long-range transport from the Middle East and Central Asia; however, the central Himalayan glaciers were mainly affected by local and South Asia emissions (from Nepal, India, and China) especially during the pre-monsoon season. Overall, South Asia and West Asia were the main contributing source regions of pollutants.

Records from new high altitude weather stations reveal the meteorological conditions on Mt Everest during the 2019 monsoon. Using data from June-October, we show that the temperature, humidity, and winds announce the arrival of the monsoon, with changes that amplify with elevation. The largest change is therefore at the summit, where we estimate that monthly mean air temperature increased by 5.5 °C between June and July to average -19.1 °C over the monsoon. Such warming takes temperatures into the realm of winter conditions on much lower mountains of the mid-latitudes, illustrated with the well-known Mount Washington observatory (1,916 m; New Hampshire, USA). Although other dangers of climbing Everest may be enhanced during the monsoon, the cold induced hazard is much reduced.

Glacier-based mountain water systems are the water towers of the world, and billions are dependent on the vital services they provide. In an effort to better understand the effects of climate change on the glaciers of the Hindu Kush-Himalaya, National Geographic and Tribhuvan University conducted a multidisciplinary expedition to Mount Everest. These are the voices from that journey.

Accurate snow depth observations are critical to assess water resources. More than a billion people rely on water from snow, most of which originates in the Northern Hemisphere mountain ranges. Yet, remote sensing observations of mountain snow depth are still lacking at the large scale. Here, we show the ability of Sentinel-1 to map snow depth in the Northern Hemisphere mountains at 1 km² resolution using an empirical change detection approach. An evaluation with measurements from ~4000 sites and reanalysis data demonstrates that the Sentinel-1 retrievals capture the spatial variability between and within mountain ranges, as well as their inter-annual differences. This is showcased with the contrasting snow depths between 2017 and 2018 in the US Sierra Nevada and European Alps. With Sentinel-1 continuity ensured until 2030 and likely beyond, these findings lay a foundation for quantifying the long-term vulnerability of mountain snow-water resources to climate change., Remote sensing observations of mountain snow depth are still lacking for the Northern Hemisphere mountains. Here authors use Sentinel-1 satellite radar measurements to assess the snow depth in mountainous areas at 1 km² resolution and show that the Sentinel-1 retrievals capture the spatial variability between and within mountain ranges, as well as their inter-annual differences.

Melting snow and glacier ice in the Himalaya forms an important source of water for people downstream. Incoming longwave radiation (LWin) is an important energy source for melt, but there are only few measurements of LWin at high elevation. For the modelling of snow and glacier melt, the LWin is therefore often represented by parameterizations that were originally developed for lower elevation environments. With LWin measurements at eight stations in three catchments in the Himalaya, with elevations between 3,980 and 6,352 m.a.s.l., we test existing LWin parameterizations. We find that these parameterizations generally underestimate the LWin, especially in wet (monsoon) conditions, where clouds are abundant and locally formed. We present a new parameterization based only on near‐surface temperature and relative humidity, both of which are easy and inexpensive to measure accurately. The new parameterization performs better than the parameterizations available in literature, in some cases halving the root‐mean‐squared error. The new parameterization is especially improving existing parameterizations in cloudy conditions. We also show that the choice of longwave parameterization strongly affects melt calculations of snow and ice., We assess the performance of incoming longwave radiation models in the Himalaya, using data from eight automatic weather stations. In the Himalaya, clouds often form locally, and we show that a new model, which uses relative humidity as a parameter, performs better than the ones available in previous literature. Source: Usmar Helleman.

We derive the surface and basal radar reflectance and backscatter coefficients of the southern McMurdo Ice Shelf (SMIS) and part of the nearby Ross Ice Shelf (RIS), Antarctica, from radar statistical reconnaissance using a 60-MHZ airborne survey. The surface coefficients are further inverted in terms of snow density and roughness, providing a spatial distribution of the processes contributing to the surface boundary conditions. We disentangle the basal coefficients from surface transmission losses, and we provide the basal coherent content, an indicator of the boundary geometric disorder that is also self-corrected from englacial attenuation. The basal radar properties exhibit sharp gradients along specific iso-depths, suggesting an abrupt modification of the ice composition and geometric structure. We interpret this behavior as locations where the pressure-melting point is reached, outlining fields of freezing and melting ice. Basal steps are observed at both SMIS and RIS, suggesting a common geometric expression of widespread basal processes. This technique offers a simultaneous view of both the surface and basal boundary conditions to help investigate the ice-shelf stability, while its application to airborne data significantly improves coverage of the difficult-to-observe ice–ocean boundary. It also provides constraints on thermohaline circulation in ice shelves cavities, which are analogs for ice-covered ocean worlds.

Ice shelves are widely known to slow the transfer of Antarctic grounded ice to the ocean, especially if their flow is decelerated by local pinning points. Their longevity is influenced by variations in ice dynamics, surface accumulation and oceanic conditions in the ice shelf cavity. This is reflected in the ice shelf structure, which can be characterized by the shape of internal radar reflection horizons.We aim to map the internal ice shelf stratigraphy of ice shelves, starting with the narrow belt of ice-shelves in the Dronning Maud Land area. The final goal will be to evaluate these as a spatiotemporal archive of ice provenance and ice dynamics. The bulk of the data presented here were collected with AWI’s airborne multi frequency ultra-wideband radar and we combine these new observations with airborne and ground-based radar surveys from previous years. We present a consistent set of internal radar isochrones on a catchment scale for the Roi Baudoin area including the Ragnhild ice streams, the grounding-zone, the iceshelf and multiple ice rises. Using pattern matching technique we can link isochrones across different ice rises in the area, and hence provide first observational constraints on how ice rises jointly react to changes in atmospheric and oceanographic forcings. We also find a number of interesting features including dynamically induced dips in shear zones, truncating layers at the ice-shelf base, and the development of a meteoric ice layer distinguishing advected from newly accumulated ice in the iceshelf. The time series provided by radar observations over the last 10 years also offers the potential to map temporal changes. We use ice-flow modelling to provide age constraints for some internal layers and delineate portions within the shelf as a function of their advection history, hence marking areas of differing rheologies within the shelf. Taken together, this case study on a catchment scale is a primer to unravel the information stored in the isochronal stratigraphy of coastal Antarctica and contributes to international efforts (e.g., SCAR AntArchitecture) of mapping stratigraphy on ice sheet scales.

Ice shelves buttress ice flow from the continent towards the ocean, and their disintegration results in increased ice discharge. Ice-shelf evolution and integrity is influenced by surface accumulation, basal melting, and ice dynamics. We find signals of all of these processes imprinted in the ice-shelf stratigraphy that can be mapped using isochrones imaged with radar.Our aim is to develop an inverse approach to infer ice shelf basal melt rates using radar isochrones as observational constraints. Here, we investigate the influence of basalt melt rates on the shape of isochrones using combined insights from both forward and inverse modeling. We use the 3D full Stokes model Elmer/Ice in our forward simulations, aiming to reproduce isochrone patterns observed in our data. Moreover we develop an inverse approach based on the shallow shelf approximating, aiming to constrain basal melt rates using isochronal radar data and surface velocities. Insights obtained from our simulations can also guide the collection of new radar data (e.g., profile lines along vs. across-flow) in a way that ambiguities in interpreting the ice-shelf stratigraphy can be minimized. Eventually, combining these approaches will enable us to better constrain the magnitude and history of basal melting, which will give valuable input for ocean circulation and sea level rise projections.

This case study provides a framework for future monitoring and evidence for human source pollution in the Khumbu region, Nepal. We analyzed the chemical composition (major ions, major/trace elements, black carbon, and stable water isotopes) of pre-monsoon stream water (4300-5250 m) and snow (5200-6665 m) samples collected from Mt. Everest, Mt. Lobuche, and the Imja Valley during the 2019 pre-monsoon season, in addition to a shallow ice core recovered from the Khumbu Glacier (5300 m). In agreement with previous work, pre-monsoon aerosol deposition is dominated by dust originating from western sources and less frequently by transport from southerly air mass sources as demonstrated by evidence of one of the strongest recorded pre-monsoon events emanating from the Bay of Bengal, Cyclone Fani. Elevated concentrations of human-sourced metals (e.g., Pb, Bi, As) are found in surface snow and stream chemistry collected in the Khumbu region. As the most comprehensive case study of environmental chemistry in the Khumbu region, this research offers sufficient evidence for increased monitoring in this watershed and surrounding areas.

Snow dynamics play a crucial role in the hydrology of alpine catchments in the Himalaya. However, studies based on in-situ observations that elucidate the energy and mass balance of the snowpack at high altitude in this region are scarce. In this study, we use meteorological and snow observations at two high-altitude sites in the Nepalese Himalaya to quantify the mass and energy balance of the seasonal snowpack. Using a data driven experimental set-up we aim to understand the main meteorological drivers of snowmelt, illustrate the importance of accounting for the cold content dynamics of the snowpack, and gain insight into the role that snow meltwater refreezing plays in the energy and mass balance of the snowpack. Our results show an intricate relation between the sensitivity of melt and refreezing on the albedo, the importance of meltwater refreezing, and the amount of positive net energy used to overcome the cold content of the snowpack. The net energy available at both sites is primarily driven by the net shortwave radiation, and is therefore extremely sensitive to snow albedo measurements. We conclude that, based on observed snowpack temperatures, 21% of the net positive energy is used to overcome the cold content build up during the night. We also show that at least 32–34% of the snow meltwater refreezes again for both sites. Even when the cold content and refreezing are accounted for, excess energy is available beyond what is needed to melt the snowpack. We hypothesize that this excess energy may be explained by uncertainties in the measurement of shortwave radiation, an underestimation of refreezing due to a basal ice layer, a cold content increase due to fresh snowfall and the ground heat flux. Our study shows that in order to accurately simulate the mass balance of seasonal snowpacks in Himalayan catchments, simple temperature index models do not suffice and refreezing and the cold content needs to be accounted for.

A layer of debris cover often accumulates across the surface of glaciers in active mountain ranges with exceptionally steep terrain, such as the Andes, Himalaya, and New Zealand Alps. Such a supraglacial debris layer has a major influence on a glacier’s surface energy budget, enhancing radiation absorption, and melt when the layer is thin, but insulating the ice when thicker than a few cm. Information on spatially distributed debris surface temperature has the potential to provide insight into the properties of the debris, its effects on the ice below and its influence on the near-surface boundary layer. Here, we deploy an unmanned aerial vehicle (UAV) equipped with a thermal infrared sensor on three separate missions over one day to map changing surface temperatures across the debris-covered Lirung Glacier in the Central Himalaya. We present a methodology to georeference and process the acquired thermal imagery, and correct for emissivity and sensor bias. Derived UAV surface temperatures are compared with distributed simultaneous in situ temperature measurements as well as with Landsat 8 thermal satellite imagery. Results show that the UAV-derived surface temperatures vary greatly both spatially and temporally, with −1.4 ± 1.8, 11.0 ± 5.2, and 15.3 ± 4.7◦C for the three flights (mean ± sd), respectively. The range in surface temperatures over the glacier during the morning is very large with almost 50 ◦C. Ground-based measurements are generally in agreement with the UAV imagery, but considerable deviations are present that are likely due to differences in measurement technique and approach, and validation is difficult as a result. The difference in spatial and temporal variability captured by the UAV as compared with much coarser satellite imagery is striking and it shows that satellite derived temperature maps should be interpreted with care. We conclude that UAVs provide a suitable means to acquire surface temperature maps of debris-covered glacier surfaces at high spatial and temporal resolution, but that there are caveats with regard to absolute temperature measurement.

Marine ice accretes at the base of ice shelves, often infilling open structural weaknesses and increasing ice-shelf stability. However, the timing and location of marine ice formation remain poorly understood. This study determines marine ice source water composition and origin by examining marine ice crystal morphology, water isotope and solute chemistry in ice samples collected from the southern McMurdo Ice Shelf (SMIS), Antarctica. The measured co-isotopic record together with the output of a freezing model for frazil crystals indicate a spatio-temporally varying water source of sea water and relatively fresher water, such as melted meteoric or marine ice. This is in agreement with the occurrence of primarily banded and granular ice crystal facies typical for frazil ice crystals that nucleate in a supercooled mixture of water masses. We propose that marine ice exposed at the surface of SMIS, which experiences summer melt, is routed to the ice-shelf base via the tide crack. Here frazil crystals nucleate in a double diffusion mechanism of heat and salt between two water masses at their salinity-dependent freezing point. Recycling of previously formed marine ice facilitates ice-shelf self-sustenance in a warming climate.