The trough-system algorithm and its application to spatial modeling of Greenland subglacial topography

This is the published version. Copyright 2014 Herzfeld et al.Dynamic ice-sheet models are used to assess the contribution of mass loss from the Greenland ice sheet to sea-level rise. Mass transfer from ice sheet to ocean is in a large part through outlet glaciers. Bed topography plays an important role in ice dynamics, since the acceleration from the slow-moving inland ice to an ice stream is in many cases caused by the existence of a subglacial trough or trough system. Problems are that most subglacial troughs are features of a scale not resolved in most ice-sheet models and that radar measurements of subglacial topography do not always reach the bottoms of narrow troughs. The trough-system algorithm introduced here employs mathematical morphology and algebraic topology to correctly represent subscale features in a topographic generalization, so the effects of troughs on ice flow are retained in ice-dynamic models. The algorithm is applied to derive a spatial elevation model of Greenland subglacial topography, integrating recently collected radar measurements (CReSIS data) of the Jakobshavn Isbræ, Helheim, Kangerdlussuaq and Petermann glacier regions. The resultant JakHelKanPet digital elevation model has been applied in dynamic ice-sheet modeling and sea-level-rise assessment

Practices in source code sharing in astrophysics

While software and algorithms have become increasingly important in astronomy, the majority of authors who publish computational astronomy research do not share the source code they develop, making it difficult to replicate and reuse the work. In this paper we discuss the importance of sharing scientific source code with the entire astrophysics community, and propose that journals require authors to make their code publicly available when a paper is published. That is, we suggest that a paper that involves a computer program not be accepted for publication unless the source code becomes publicly available. The adoption of such a policy by editors, editorial boards, and reviewers will improve the ability to replicate scientific results, and will also make the computational astronomy methods more available to other researchers who wish to apply them to their data.Comment: Accepted by Astronomy and Computing. 10 page

Algorithm for Detection of Ground and Canopy Cover in Micropulse Photon-Counting Lidar Altimeter Data in Preparation for the ICESat-2 Mission

NASA's Ice, Cloud and Land Elevation Satellite-II (ICESat-2) mission is a decadal survey mission (2016 launch). The mission objectives are to measure land ice elevation, sea ice freeboard, and changes in these variables, as well as to collect measurements over vegetation to facilitate canopy height determination. Two innovative components will characterize the ICESat-2 lidar: 1) collection of elevation data by a multibeam system and 2) application of micropulse lidar (photon-counting) technology. A photon-counting altimeter yields clouds of discrete points, resulting from returns of individual photons, and hence new data analysis techniques are required for elevation determination and association of the returned points to reflectors of interest. The objective of this paper is to derive an algorithm that allows detection of ground under dense canopy and identification of ground and canopy levels in simulated ICESat-2 data, based on airborne observations with a Sigma Space micropulse lidar. The mathematical algorithm uses spatial statistical and discrete mathematical concepts, including radial basis functions, density measures, geometrical anisotropy, eigenvectors, and geostatistical classification parameters and hyperparameters. Validation shows that ground and canopy elevation, and hence canopy height, can be expected to be observable with high accuracy by ICESat-2 for all expected beam energies considered for instrument design (93.01%-99.57% correctly selected points for a beam with expected return of 0.93 mean signals per shot (msp), and 72.85%-98.68% for 0.48 msp). The algorithm derived here is generally applicable for elevation determination from photoncounting lidar altimeter data collected over forested areas, land ice, sea ice, and land surfaces, as well as for cloud detection

Achieving superresolution with illumination-enhanced sparsity.

Recent advances in superresolution fluorescence microscopy have been limited by a belief that surpassing two-fold resolution enhancement of the Rayleigh resolution limit requires stimulated emission or the fluorophore to undergo state transitions. Here we demonstrate a new superresolution method that requires only image acquisitions with a focused illumination spot and computational post-processing. The proposed method utilizes the focused illumination spot to effectively reduce the object size and enhance the object sparsity and consequently increases the resolution and accuracy through nonlinear image post-processing. This method clearly resolves 70nm resolution test objects emitting ~530nm light with a 1.4 numerical aperture (NA) objective, and, when imaging through a 0.5NA objective, exhibits high spatial frequencies comparable to a 1.4NA widefield image, both demonstrating a resolution enhancement above two-fold of the Rayleigh resolution limit. More importantly, we examine how the resolution increases with photon numbers, and show that the more-than-two-fold enhancement is achievable with realistic photon budgets

Effects of Non-Aerobic Maximal Effort Exercise on Fatigue in Deconditioned Men and Women with Multiple Sclerosis

Multiple Sclerosis (MS) is a neurodegenerative disease of unknown etiology affecting women more frequently than men. Mental and physical fatigue complaints are often the most disabling symptoms for an MS patient. Both are multifactorial, potentially exacerbated by aerobic exercise, may prevent sustained physical functioning, and significantly interfere with activities of daily living1. A multi-center study was designed to investigate the effects of non-aerobic maximal effort exercise (MEE) for deconditioned persons with MS, with the expectation of minimizing fatigue. The IsoPUMP (Neuromuscular Engineering; Nashville, TN), is a specialized exercise and strength-sensing machine, designed to allow individuals to safely perform and record their non-aerobic MEE sessions. The Modified Fatigue Impact Scale (MFIS) and Multiple Sclerosis Functional Composite (MSFC) are common, accepted methods used to measure fatigue and function. The MFIS is a 21-item questionnaire which assesses the subjects’ perception of physical, cognitive, and psychosocial aspects of fatigue over a four-week period2. Each of the 21 items are scored on a scale from 0 (never) to 4 (almost always), and the total MFIS score is calculated by summing the circled number for each item. Total scores can range from 0 to 84; higher scores indicating a greater impact of fatigue on the person. The MFIS has three distinct subscales: (1) physical, (2) cognitive, and (3) psychosocial. These subscales can be scored independently by summing the questions that pertain to each subscale2. The MFIS physical subscale score can range from 0 – 36 and the MFIS cognitive subscale score can range from 0 – 40. The MSFC combines clinical measures used to assess lower limb function (Timed 25-Foot Walk [25-FW]), upper limb function (9-Hole Peg Test [9-HPT]), and cognition (Paced Auditory Serial Addition Test [PASAT-3”])3. The 25-FW is a quantitative measure of lower extremity function. The 9-HPT is a quantitative measure of arm and hand function where a subject inserts and then removes 9 pegs from a board, using one hand at a time. The time is recorded for each hand with the dominant hand trial first and the non-dominant hand trial second. The final score is recorded as the mean time for both hands. The PASAT-3” is a measure of cognitive function, specifically assessing auditory information processing speed, short-term memory, flexibility, and calculation ability. Cognitive dysfunction affects half of all MS patients; slowing ability to reason, concentrate, and recall5. In this test subjects listen to a series of 61 spoken numbers separated by 3 seconds and must add each number to the prior number. Their final PASAT-3” score is the number of correct additions in the series, with 60 reflecting a perfect score. The MSFC is then evaluated by creating Z-scores for each component, which compare each outcome with the average outcome of the study population. The three Z-scores are then averaged to create an overall composite score (the MSFC score) which represents change over time for that population of MS subjects3

An Analysis of Functional Status in Multiple Sclerosis Patients after Progressive Non-Aerobic High-Intensity Maximal Effort Exercise (MEE)

Background: Multiple Sclerosis (MS) is a disease with a wide-ranging impact on functional status. MS patient function has been assessed using Multiple Sclerosis Functional Composite Score (MSFCS). The MSFCS includes the standardized scores (Z-score) of three functional tests: the Paced Auditory Serial Addition Test (PASAT-3”) for cognitive function, 9-Hole Peg Test (9-HPT) for upper extremity function, and timed 25-foot walk (25-TW) for lower extremity function. One of the most common symptoms experienced by MS patients is severe fatigue, often brought on suddenly by aerobic exercise. Non-aerobic maximal effort exercise (MEE) is thought to increase strength without increasing fatigue. The IsoPUMP® (Neuromuscular Engineering; Nashville, TN) is a stationary exercise device designed for patient use to safely perform MEE leg presses and whole body lunges using isometric and eccentric exercises. The progressive functional changes of the MS patients were tracked using the MSFCs at specific intervals during the study

Global, regional, and national comparative risk assessment of 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks in 188 countries, 1990-2013: A systematic analysis for the Global Burden of Disease Study 2013

Background: The Global Burden of Disease, Injuries, and Risk Factor study 2013 (GBD 2013) is the first of a series of annual updates of the GBD. Risk factor quantification, particularly of modifiable risk factors, can help to identify emerging threats to population health and opportunities for prevention. The GBD 2013 provides a timely opportunity to update the comparative risk assessment with new data for exposure, relative risks, and evidence on the appropriate counterfactual risk distribution. Methods: Attributable deaths, years of life lost, years lived with disability, and disability-adjusted life-years (DALYs) have been estimated for 79 risks or clusters of risks using the GBD 2010 methods. Risk-outcome pairs meeting explicit evidence criteria were assessed for 188 countries for the period 1990-2013 by age and sex using three inputs: risk exposure, relative risks, and the theoretical minimum risk exposure level (TMREL). Risks are organised into a hierarchy with blocks of behavioural, environmental and occupational, and metabolic risks at the first level of the hierarchy. The next level in the hierarchy includes nine clusters of related risks and two individual risks, with more detail provided at levels 3 and 4 of the hierarchy. Compared with GBD 2010, six new risk factors have been added: handwashing practices, occupational exposure to trichloroethylene, childhood wasting, childhood stunting, unsafe sex, and low glomerular filtration rate. For most risks, data for exposure were synthesised with a Bayesian metaregression method, DisMod-MR 2.0, or spatial-temporal Gaussian process regression. Relative risks were based on meta-regressions of published cohort and intervention studies. Attributable burden for clusters of risks and all risks combined took into account evidence on the mediation of some risks such as high body-mass index (BMI) through other risks such as high systolic blood pressure and high cholesterol. Findings: All risks combined account for 57·2% (95% uncertainty interval [UI] 55·8-58·5) of deaths and 41·6% (40·1-43·0) of DALYs. Risks quantified account for 87·9% (86·5-89·3) of cardiovascular disease DALYs, ranging to a low of 0% for neonatal disorders and neglected tropical diseases and malaria. In terms of global DALYs in 2013, six risks or clusters of risks each caused more than 5% of DALYs: dietary risks accounting for 11·3 million deaths and 241·4 million DALYs, high systolic blood pressure for 10·4 million deaths and 208·1 million DALYs, child and maternal malnutrition for 1·7 million deaths and 176·9 million DALYs, tobacco smoke for 6·1 million deaths and 143·5 million DALYs, air pollution for 5·5 million deaths and 141·5 million DALYs, and high BMI for 4·4 million deaths and 134·0 million DALYs. Risk factor patterns vary across regions and countries and with time. In sub-Saharan Africa, the leading risk factors are child and maternal malnutrition, unsafe sex, and unsafe water, sanitation, and handwashing. In women, in nearly all countries in the Americas, north Africa, and the Middle East, and in many other high-income countries, high BMI is the leading risk factor, with high systolic blood pressure as the leading risk in most of Central and Eastern Europe and south and east Asia. For men, high systolic blood pressure or tobacco use are the leading risks in nearly all high-income countries, in north Africa and the Middle East, Europe, and Asia. For men and women, unsafe sex is the leading risk in a corridor from Kenya to South Africa. Interpretation: Behavioural, environmental and occupational, and metabolic risks can explain half of global mortality and more than one-third of global DALYs providing many opportunities for prevention. Of the larger risks, the attributable burden of high BMI has increased in the past 23 years. In view of the prominence of behavioural risk factors, behavioural and social science research on interventions for these risks should be strengthened. Many prevention and primary care policy options are available now to act on key risks

Combination of searches for heavy spin-1 resonances using 139 fb−1 of proton-proton collision data at s = 13 TeV with the ATLAS detector

A combination of searches for new heavy spin-1 resonances decaying into different pairings of W, Z, or Higgs bosons, as well as directly into leptons or quarks, is presented. The data sample used corresponds to 139 fb−1 of proton-proton collisions at = 13 TeV collected during 2015–2018 with the ATLAS detector at the CERN Large Hadron Collider. Analyses selecting quark pairs (qq, bb, , and tb) or third-generation leptons (τν and ττ) are included in this kind of combination for the first time. A simplified model predicting a spin-1 heavy vector-boson triplet is used. Cross-section limits are set at the 95% confidence level and are compared with predictions for the benchmark model. These limits are also expressed in terms of constraints on couplings of the heavy vector-boson triplet to quarks, leptons, and the Higgs boson. The complementarity of the various analyses increases the sensitivity to new physics, and the resulting constraints are stronger than those from any individual analysis considered. The data exclude a heavy vector-boson triplet with mass below 5.8 TeV in a weakly coupled scenario, below 4.4 TeV in a strongly coupled scenario, and up to 1.5 TeV in the case of production via vector-boson fusion

Global, regional, and national comparative risk assessment of 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks, 1990-2015: A systematic analysis for the Global Burden of Disease Study 2015

Background: The Global Burden of Diseases, Injuries, and Risk Factors Study 2015 provides an up-to-date synthesis of the evidence for risk factor exposure and the attributable burden of disease. By providing national and subnational assessments spanning the past 25 years, this study can inform debates on the importance of addressing risks in context. Methods: We used the comparative risk assessment framework developed for previous iterations of the Global Burden of Disease Study to estimate attributable deaths, disability-adjusted life-years (DALYs), and trends in exposure by age group, sex, year, and geography for 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks from 1990 to 2015. This study included 388 risk-outcome pairs that met World Cancer Research Fund-defined criteria for convincing or probable evidence. We extracted relative risk and exposure estimates from randomised controlled trials, cohorts, pooled cohorts, household surveys, census data, satellite data, and other sources. We used statistical models to pool data, adjust for bias, and incorporate covariates. We developed a metric that allows comparisons of exposure across risk factors—the summary exposure value. Using the counterfactual scenario of theoretical minimum risk level, we estimated the portion of deaths and DALYs that could be attributed to a given risk. We decomposed trends in attributable burden into contributions from population growth, population age structure, risk exposure, and risk-deleted cause-specific DALY rates. We characterised risk exposure in relation to a Socio-demographic Index (SDI). Findings: Between 1990 and 2015, global exposure to unsafe sanitation, household air pollution, childhood underweight, childhood stunting, and smoking each decreased by more than 25%. Global exposure for several occupational risks, high body-mass index (BMI), and drug use increased by more than 25% over the same period. All risks jointly evaluated in 2015 accounted for 57·8% (95% CI 56·6–58·8) of global deaths and 41·2% (39·8–42·8) of DALYs. In 2015, the ten largest contributors to global DALYs among Level 3 risks were high systolic blood pressure (211·8 million [192·7 million to 231·1 million] global DALYs), smoking (148·6 million [134·2 million to 163·1 million]), high fasting plasma glucose (143·1 million [125·1 million to 163·5 million]), high BMI (120·1 million [83·8 million to 158·4 million]), childhood undernutrition (113·3 million [103·9 million to 123·4 million]), ambient particulate matter (103·1 million [90·8 million to 115·1 million]), high total cholesterol (88·7 million [74·6 million to 105·7 million]), household air pollution (85·6 million [66·7 million to 106·1 million]), alcohol use (85·0 million [77·2 million to 93·0 million]), and diets high in sodium (83·0 million [49·3 million to 127·5 million]). From 1990 to 2015, attributable DALYs declined for micronutrient deficiencies, childhood undernutrition, unsafe sanitation and water, and household air pollution; reductions in risk-deleted DALY rates rather than reductions in exposure drove these declines. Rising exposure contributed to notable increases in attributable DALYs from high BMI, high fasting plasma glucose, occupational carcinogens, and drug use. Environmental risks and childhood undernutrition declined steadily with SDI; low physical activity, high BMI, and high fasting plasma glucose increased with SDI. In 119 countries, metabolic risks, such as high BMI and fasting plasma glucose, contributed the most attributable DALYs in 2015. Regionally, smoking still ranked among the leading five risk factors for attributable DALYs in 109 countries; childhood underweight and unsafe sex remained primary drivers of early death and disability in much of sub-Saharan Africa. Interpretation: Declines in some key environmental risks have contributed to declines in critical infectious diseases. Some risks appear to be invariant to SDI. Increasing risks, including high BMI, high fasting plasma glucose, drug use, and some occupational exposures, contribute to rising burden from some conditions, but also provide opportunities for intervention. Some highly preventable risks, such as smoking, remain major causes of attributable DALYs, even as exposure is declining. Public policy makers need to pay attention to the risks that are increasingly major contributors to global burden. Funding: Bill & Melinda Gates Foundation