A Hierarchical Spatio-Temporal Statistical Model Motivated by Glaciology

In this paper, we extend and analyze a Bayesian hierarchical spatio-temporal model for physical systems. A novelty is to model the discrepancy between the output of a computer simulator for a physical process and the actual process values with a multivariate random walk. For computational efficiency, linear algebra for bandwidth limited matrices is utilized, and first-order emulator inference allows for the fast emulation of a numerical partial differential equation (PDE) solver. A test scenario from a physical system motivated by glaciology is used to examine the speed and accuracy of the computational methods used, in addition to the viability of modeling assumptions. We conclude by discussing how the model and associated methodology can be applied in other physical contexts besides glaciology.Comment: Revision accepted for publication by the Journal of Agricultural, Biological, and Environmental Statistic

Spring warming in Yukon mountains is not amplified by the snow albedo feedback

Decreasing spring snow cover may amplify Arctic warming through the snow albedo feedback. To examine the impact of snowmelt on increasing temperature we used a 5,000 m elevation gradient in Yukon, Canada, extending from valley-bottom conifer forests, through middle elevation tundra, to high elevation icefields, to compare validated downscaled reanalysis air temperature patterns across elevational bands characterized by different patterns of spring snowmelt. From 2000 to 2014 we observed surface warming of 0.01 °C/a·1,000 m in May (0.14 °C/a at 1,000 m to 0.19 °C/a at 5,000 m), and uniform cooling of 0.09 °C/a in June at all elevations. May temperature trends across elevationally dependent land cover types were highly correlated with each other despite large variations in albedo and snow cover trends. Furthermore, a clear dependency of infrared skin temperature on snow cover mediated albedo decline was observed in tundra, but this was insufficient to influence average diurnal air temperature. We observed negative June temperature trends which we attribute to increasing daytime cloud cover because albedo and snow cover trends were unchanging. We conclude that 8-day and monthly averaged Spring air temperature trends are responding to a synoptic external forcing that is much stronger than the snow albedo feedback in sub-Arctic mountains

A Bayesian hierarchical model for glacial dynamics based on the shallow ice approximation and its evaluation using analytical solutions

Bayesian hierarchical modeling can assist the study of glacial dynamics and ice flow properties. This approach will allow glaciologists to make fully probabilistic predictions for the thickness of a glacier at unobserved spatiotemporal coordinates, and it will also allow for the derivation of posterior probability distributions for key physical parameters such as ice viscosity and basal sliding. The goal of this paper is to develop a proof of concept for a Bayesian hierarchical model constructed, which uses exact analytical solutions for the shallow ice approximation (SIA) introduced by Bueler et al. (2005). A suite of test simulations utilizing these exact solutions suggests that this approach is able to adequately model numerical errors and produce useful physical parameter posterior distributions and predictions. A byproduct of the development of the Bayesian hierarchical model is the derivation of a novel finite difference method for solving the SIA partial differential equation (PDE). An additional novelty of this work is the correction of numerical errors induced through a numerical solution using a statistical model. This error-correcting process models numerical errors that accumulate forward in time and spatial variation of numerical errors between the dome, interior, and margin of a glacier.The Icelandic Research Fund (RANNIS) is thanked for funding this research.Peer Reviewe

Delivery of stable ultra-thin liquid sheets in vacuum for biochemical spectroscopy

The development of ultra-thin flat liquid sheets capable of running in vacuum has provided an exciting new target for X-ray absorption spectroscopy in the liquid and solution phases. Several methods have become available for delivering in-vacuum sheet jets using different nozzle designs. We compare the sheets produced by two different types of nozzle; a commercially available borosillicate glass chip using microfluidic channels to deliver colliding jets, and an in-house fabricated fan spray nozzle which compresses the liquid on an axis out of a slit to achieve collision conditions. We find in our tests that both nozzles are suitable for use in X-ray absorption spectroscopy with the fan spray nozzle producing thicker but more stable jets than the commercial nozzle. We also provide practical details of how to run these nozzles in vacuum

Direct observation of ultrafast exciton localization in an organic semiconductor with soft X-ray transient absorption spectroscopy

The localization dynamics of excitons in organic semiconductors influence the efficiency of charge transfer and separation in these materials. Here we apply time-resolved X-ray absorption spectroscopy to track photoinduced dynamics of a paradigmatic crystalline conjugated polymer: poly(3-hexylthiophene) (P3HT) commonly used in solar cell devices. The π→π* transition, the first step of solar energy conversion, is pumped with a 15 fs optical pulse and the dynamics are probed by an attosecond soft X-ray pulse at the carbon K-edge. We observe X-ray spectroscopic signatures of the initially hot excitonic state, indicating that it is delocalized over multiple polymer chains. This undergoes a rapid evolution on a sub 50 fs timescale which can be directly associated with cooling and localization to form either a localized exciton or polaron pair

A Hierarchical Spatiotemporal Statistical Model Motivated by Glaciology

This is a post-peer-review, pre-copyedit version of an article published in Journal of Agricultural, Biological and Environmental Statistics. The final authenticated version is available online at: http://dx.doi.org/10.1007/s13253-019-00367-1In this paper, we extend and analyze a Bayesian hierarchical spatiotemporal model for physical systems. A novelty is to model the discrepancy between the output of a computer simulator for a physical process and the actual process values with a multivariate random walk. For computational efficiency, linear algebra for bandwidth limited matrices is utilized, and first-order emulator inference allows for the fast emulation of a numerical partial differential equation (PDE) solver. A test scenario from a physical system motivated by glaciology is used to examine the speed and accuracy of the computational methods used, in addition to the viability of modeling assumptions. We conclude by discussing how the model and associated methodology can be applied in other physical contexts besides glaciology.Icelandic Centre for Research (152457).Peer reviewe

Non-surface mass balance of glaciers in Iceland

Publisher's version (útgefin grein)Non-surface mass balance is non-negligible for glaciers in Iceland. Several Icelandic glaciers are in the neo-volcanic zone where a combination of geothermal activity, volcanic eruptions and geothermal heat flux much higher than the global average lead to basal melting close to 150 mm w.e. a−1 for the Mýrdalsjökull ice cap and 75 mm w.e. a−1 for the largest ice cap, Vatnajökull. Energy dissipation in the flow of water and ice is also rather large for the high-precipitation, temperate glaciers of Iceland resulting in internal and basal melting of 20–150 mm w.e. a−1. The total non-surface melting of glaciers in Iceland in 1995–2019 was 45–375 mm w.e. a−1 on average for the main ice caps, and was largest for Mýrdalsjökull, the south side of Vatnajökull and Eyjafjallajökull. Geothermal melting, volcanic eruptions and the energy dissipation in the flow of water and ice, as well as calving, all contribute, and thus these components should be considered in mass-balance studies. For comparison, the average mass balance of glaciers in Iceland since 1995 is −500 to −1500 mm w.e. a−1. The non-surface mass balance corresponds to a total runoff contribution of 2.1 km3 a−1 of water from Iceland.Financial support for lidar mapping of glaciers in Iceland in 2008–2012 was provided by the Icelandic Research Fund (163391-052), the Landsvirkjun (National Power Company of Iceland) Research Fund, the Icelandic Road Administration, the Reykjavík Energy Environmental and Energy Research Fund, the National Land Survey of Iceland, the Klima- og Luftgruppen (KoL) research fund of the Nordic Council of Ministers, and the Vatnajökull National Park. The acquisition of the Hofsjökull 2013 DEM was funded by AlpS GmbH and the University of Innsbruck. The acquisition of the Langjökull 2013 DEM was funded by NERC grant IG 13/12 and the DEM was provided by Ian Willis at the Scott Polar Research Institute. The work on estimating geothermal and volcanic power is based on funding from many sources, including the Research Fund of the University of Iceland, ISAVIA (the Icelandic Aviation Service), the Icelandic Road Administration and Landsvirkjun; logistical support has been provided by the Iceland Glaciological Society.Peer Reviewe