Basal conditions of two Transantarctic Mountains outlet glaciers from observation-constrained diagnostic modelling

This is the published version. Copyright 2014 International Glaciological SocietyWe present a diagnostic glacier flowline model parameterized and constrained by new velocity data from ice-surface GPS installations and speckle tracking of TerraSAR-X satellite images, newly acquired airborne-radar data, and continental gridded datasets of topography and geothermal heat flux, in order to better understand two outlet glaciers of the East Antarctic ice sheet. Our observational data are employed as primary inputs to a modelling procedure that first calculates the basal thermal regime of each glacier, then iterates the basal sliding coefficient and deformation rate parameter until the fit of simulated to observed surface velocities is optimized. We find that the two glaciers have both frozen and thawed areas at their beds, facilitating partial sliding. Glacier flow arises from a balance between sliding and deformation that fluctuates along the length of each glacier, with the amount of sliding typically varying by up to two orders of magnitude but with deformation rates far more constant. Beardmore Glacier is warmer and faster-flowing than Skelton Glacier, but an up-glacier deepening bed at the grounding line, coupled with ice thicknesses close to flotation, lead us to infer a greater vulnerability of Skelton Glacier to grounding-line recession if affected by ocean-forced thinning and concomitant acceleration

Differential interferometric synthetic aperture radar for tide modelling in Antarctic ice-shelf grounding zones

Differential interferometric synthetic aperture radar (DInSAR) is an essential tool for detecting ice-sheet motion near Antarctica's oceanic margin. These space-borne measurements have been used extensively in the past to map the location and retreat of ice-shelf grounding lines as an indicator for the onset of marine ice-sheet instability and to calculate the mass balance of ice sheets and individual catchments. The main difficulty in interpreting DInSAR is that images originate from a combination of several SAR images and do not indicate instantaneous ice deflection at the times of satellite data acquisitions. Here, we combine the sub-centimetre accuracy and spatial benefits of DInSAR with the temporal benefits of tide models to infer the spatio-temporal dynamics of ice–ocean interaction during the times of satellite overpasses. We demonstrate the potential of this synergy with TerraSAR-X data from the almost-stagnant southern McMurdo Ice Shelf (SMIS). We then validate our algorithm with GPS data from the fast-flowing Darwin Glacier, draining the Antarctic Plateau through the Transantarctic Mountains into the Ross Sea. We are able to reconstruct DInSAR-derived vertical displacements to 7 mm mean absolute residual error and generally improve traditional tide-model output by up to 39 % from 10.8 to 6.7 cm RMSE against GPS data from areas where ice is in local hydrostatic equilibrium with the ocean and by up to 74 % from 21.4 to 5.6 cm RMSE against GPS data in feature-rich coastal areas where tide models have not been applicable before. Numerical modelling then reveals Young's modulus of E=1.0±0.56  GPa and an ice viscosity of ν=10±3.65  TPa s when finite-element simulations of tidal flexure are matched to 16 d of tiltmeter data, supporting the hypothesis that strain-dependent anisotropy may significantly decrease effective viscosity compared to isotropic polycrystalline ice on large spatial scales. Applications of our method include the following: refining coarsely gridded tide models to resolve small-scale features at the spatial resolution and vertical accuracy of SAR imagery, separating elastic and viscoelastic contributions in the satellite-derived flexure measurement, and gaining information about large-scale ice heterogeneity in Antarctic ice-shelf grounding zones, the missing key to improving current ice-sheet flow models. The reconstruction of the individual components forming DInSAR images has the potential to become a standard remote-sensing method in polar tide modelling. Unlocking the algorithm's full potential to answer multi-disciplinary research questions is desired and demands collaboration within the scientific community

High basal melting forming a channel at the grounding line of Ross Ice Shelf, Antarctica

Antarctica's ice shelves are thinning at an increasing rate, affecting their buttressing ability. Channels in the ice shelf base unevenly distribute melting, and their evolution provides insight into changing subglacial and oceanic conditions. Here we used phase-sensitive radar measurements to estimate basal melt rates in a channel beneath the currently stable Ross Ice Shelf. Melt rates of 22.2 ± 0.2 m a−1 (>2500% the overall background rate) were observed 1.7 km seaward of Mercer/Whillans Ice Stream grounding line, close to where subglacial water discharge is expected. Laser altimetry shows a corresponding, steadily deepening surface channel. Two relict channels to the north suggest recent subglacial drainage reorganization beneath Whillans Ice Stream approximately coincident with the shutdown of Kamb Ice Stream. This rapid channel formation implies that shifts in subglacial hydrology may impact ice shelf stability

Building a traceable climate model hierarchy with multi-level emulators

To study climate change on multi-millennial timescales or to explore a model’s parameter space, efficient models with simplified and parameterised processes are required. However, the reduction in explicitly modelled processes can lead to underestimation of some atmospheric responses that are essential to the understanding of the climate system. While more complex general circulations are available and capable of simulating a more realistic climate, they are too computationally intensive for these purposes. In this work, we propose a multi-level Gaussian emulation technique to efficiently estimate the outputs of steady-state simulations of an expensive atmospheric model in response to changes in boundary forcing. The link between a computationally expensive atmospheric model, PLASIM (Planet Simulator), and a cheaper model, EMBM (energy–moisture balance model), is established through the common boundary condition specified by an ocean model, allowing for information to be propagated from one to the other. This technique allows PLASIM emulators to be built at a low cost. The method is first demonstrated by emulating a scalar summary quantity, the global mean surface air temperature. It is then employed to emulate the dimensionally reduced 2-D surface air temperature field. Even though the two atmospheric models chosen are structurally unrelated, Gaussian process emulators of PLASIM atmospheric variables are successfully constructed using EMBM as a fast approximation. With the extra information gained from the cheap model, the multi-level emulator of PLASIM’s 2-D surface air temperature field is built using only one-third the amount of expensive data required by the normal single-level technique. The constructed emulator is shown to capture 93.2% of the variance across the validation ensemble, with the averaged RMSE of 1.33 °C. Using the method proposed, quantities from PLASIM can be constructed and used to study the effects introduced by PLASIM’s atmosphere

Ultra-light axions and the S8S_8 tension: joint constraints from the cosmic microwave background and galaxy clustering

We search for ultra-light axions as dark matter (DM) and dark energy particle candidates, for axion masses $10^{-32}\,\mathrm{eV} \leq m_\mathrm{a} \leq 10^{-24}\,\mathrm{eV}$, by a joint analysis of cosmic microwave background (CMB) and galaxy clustering data -- and consider if axions can resolve the tension in inferred values of the matter clustering parameter $S_8$. We give legacy constraints from Planck 2018 CMB data, improving 2015 limits on the axion density $\Omega_\mathrm{a} h^2$ by up to a factor of three; CMB data from the Atacama Cosmology Telescope and the South Pole Telescope marginally weaken Planck bounds at $m_\mathrm{a} = 10^{-25}\,\mathrm{eV}$, owing to lower (and theoretically-consistent) gravitational lensing signals. We jointly infer, from Planck CMB and full-shape galaxy power spectrum and bispectrum data from the Baryon Oscillation Spectroscopic Survey (BOSS), that axions are, today, $< 10\%$ of the DM for $m_\mathrm{a} \leq 10^{-26}\,\mathrm{eV}$ and $< 1\%$ for $10^{-30}\,\mathrm{eV} \leq m_\mathrm{a} \leq 10^{-28}\,\mathrm{eV}$. BOSS data strengthen limits, in particular at higher $m_\mathrm{a}$ by probing high-wavenumber modes ($k < 0.4 h\,\mathrm{Mpc}^{-1}$). BOSS alone finds a preference for axions at $2.7 \sigma$, for $m_\mathrm{a} = 10^{-26}\,\mathrm{eV}$, but Planck disfavours this result. Nonetheless, axions in a window $10^{-28}\,\mathrm{eV} \leq m_\mathrm{a} \leq 10^{-25}\,\mathrm{eV}$ can improve consistency between CMB and galaxy clustering data, e.g., reducing the $S_8$ discrepancy from $2.7 \sigma$ to $1.6 \sigma$, since these axions suppress structure growth at the $8 h^{-1}\,\mathrm{Mpc}$ scales to which $S_8$ is sensitive. We expect improved constraints with upcoming high-resolution CMB and galaxy lensing and future galaxy clustering data, where we will further assess if axions can restore cosmic concordance.Comment: 52 pages, 22 figure

Ratiometric electrochemical detection of Pd•••π interactions:application towards electrochemical molecular logic gates

<p>The widespread and large scale use of platinum group metals, especially palladium, in a wide variety of industrial applications has seen their levels in wastewater streams, roadside dust and even pharmaceuticals significantly rise over recent years. Due to the possible environmental damage and potential health risk this may cause, there is now substantial demand for inexpensive, efficient and robust methods for the detection of palladium. Based upon self-immolative linker technologies, we have designed and synthesised a number of allyl ether-functionalised electrochemical probes to determine the optimum probe structure required to deliver a ratiometric electrochemical detection method capable of achieving a limit of detection of <1 mg/mL within 20 min through the use of disposable screen-printed carbon electrodes. Combined with an enzymatic assay, this method was then used to achieve a proof-of-principle ratiometric electrochemical molecular logic gate.</p

Basal mass balance and prevalence of ice tongues in the Western ross sea

Ice tongues at the fringes of the Antarctic ice sheet lose mass primarily through both basal melting and calving. They are sensitive to ocean conditions which can weaken the ice both mechanically or through thinning. Ice tongues, which are laterally unconfined, are likely to be particularly sensitive to ocean-induced stresses. Here we examine ice tongues in the Western Ross Sea, by looking into the factors affecting their stability. We calculate the basal mass change of twelve Antarctic ice tongues using a flux gate approach, deriving thickness from ICESat-2 height measurements and ice surface velocities from Sentinel-1 feature-tracking over the same period (October 2018 to December 2021). The basal mass balance ranges between −0.14 ± 0.07 m yr−1 and −1.50 ± 1.2 m yr−1. The average basal mass change for all the ice tongues is −0.82 ± 0.68 m of ice yr−1. Low values of basal melt suggest a stable mass balance condition in this region, with low thermal ocean forcing, as other studies have shown. We found a heterogeneous basal melt pattern with no latitudinal gradient and no clear driver in basal melt indicating that local variables are important in the persistence of ice tongues in the absence of a strong oceanographic melting force. Moreover, thanks to the temporal resolution of the data we were able to resolve the seasonal variability of Drygalski and Aviator Ice Tongues, the two largest ice tongues studied