Search

Observations of glacier mass changes are key to understanding the response of glaciers to climate change and related impacts, such as regional runoff, ecosystem changes, and global sea-level rise. Spaceborne optical and radar sensors make it possible to quantify glacier elevation changes, and thus multi-annual mass changes, on a regional and global scale. However, estimates from a growing number of studies show a wide range of results with differences often beyond uncertainty bounds. Here, we present the outcome of a community-based inter-comparison experiment using spaceborne optical stereo (ASTER) and synthetic aperture radar interferometry (TanDEM-X) data to estimate elevation changes for defined glaciers and target periods that pose different assessment challenges. Using provided or self-processed digital elevation models (DEMs) for five test sites, 12 research groups provided a total of 97 spaceborne elevation-change datasets using various processing strategies. Validation with airborne data showed that using an ensemble estimate is promising to reduce random errors from different instruments and processing methods, but still requires a more comprehensive investigation and correction of systematic errors. We found that scene selection, DEM processing, and co-registration have the biggest impact on the results. Other processing steps, such as treating spatial data voids, differences in survey periods, or radar penetration, can still be important for individual cases. Future research should focus on testing different implementations of individual processing steps (e.g. co-registration) and addressing issues related to temporal corrections, radar penetration, glacier area changes, and density conversion. Finally, there is a clear need for our community to develop best practices, use open, reproducible software, and assess overall uncertainty in order to enhance inter-comparison and empower physical process insights across glacier elevation-change studies. [ABSTRACT FROM AUTHOR]

Unbiased estimation of causal effects with observational data requires adjustment for confounding variables that are related to both the outcome and treatment assignment. Standard variable selection techniques aim to maximize predictive ability of the outcome model, but they ignore covariate associations with treatment and may not adjust for important confounders weakly associated to outcome. We propose a novel method for estimating causal effects that simultaneously considers models for both outcome and treatment, which we call the bilevel spike and slab causal estimator (BSSCE). By using a Bayesian formulation, BSSCE estimates the posterior distribution of all model parameters and provides straightforward and reliable inference. Spike and slab priors are used on each covariate coefficient which aim to minimize the mean squared error of the treatment effect estimator. Theoretical properties of the treatment effect estimator are derived justifying the prior used in BSSCE. Simulations show that BSSCE can substantially reduce mean squared error over numerous methods and performs especially well with large numbers of covariates, including situations where the number of covariates is greater than the sample size. We illustrate BSSCE by estimating the causal effect of vasoactive therapy vs. fluid resuscitation on hypotensive episode length for patients in the Multiparameter Intelligent Monitoring in Intensive Care III critical care database. [ABSTRACT FROM AUTHOR]