Time-variable gravity signal in Greenland revealed by high-low satellite-to-satellite tracking

[1] In the event of a termination of the Gravity Recovery and Climate Experiment (GRACE) mission before the launch of GRACE Follow-On (due for launch in 2017), high-low satellite-to-satellite tracking (hl-SST) will be the only dedicated observing system with global coverage available to measure the time-variable gravity field (TVG) on a monthly or even shorter time scale. Until recently, hl-SST TVG observations were of poor quality and hardly improved the performance of Satellite Laser Ranging observations. To date, they have been of only very limited usefulness to geophysical or environmental investigations. In this paper, we apply a thorough reprocessing strategy and a dedicated Kalman filter to Challenging Minisatellite Payload (CHAMP) data to demonstrate that it is possible to derive the very long-wavelength TVG features down to spatial scales of approximately 2000 km at the annual frequency and for multi-year trends. The results are validated against GRACE data and surface height changes from long-term GPS ground stations in Greenland. We find that the quality of the CHAMP solutions is sufficient to derive long-term trends and annual amplitudes of mass change over Greenland. We conclude that hl-SST is a viable source of information for TVG and can serve to some extent to bridge a possible gap between the end-of-life of GRACE and the availability of GRACE Follow-On.

SURVIVAL PREDICTION FOR BRAIN TUMOR PATIENTS USING GENE EXPRESSION DATA

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. In order to evaluate this hypothesis, the general goal of this research is to build models for survival prediction of glioma patients using DNA molecular profiles (U133 Affymetrix gene expression microarrays) along with clinical information. First, a predictive Random Forest model is built for binary outcomes (i.e. short vs. long-term survival) and a small subset of genes whose expression values can be used to predict survival time is selected. Following, a new statistical methodology is developed for predicting time-to-death outcomes using Bayesian ensemble trees. Due to a large heterogeneity observed within prognostic classes obtained by the Random Forest model, prediction can be improved by relating time-to-death with gene expression profile directly. We propose a Bayesian ensemble model for survival prediction which is appropriate for high-dimensional data such as gene expression data. Our approach is based on the ensemble "sum-of-trees" model which is flexible to incorporate additive and interaction effects between genes. We specify a fully Bayesian hierarchical approach and illustrate our methodology for the CPH, Weibull, and AFT survival models. We overcome the lack of conjugacy using a latent variable formulation to model the covariate effects which decreases computation time for model fitting. Also, our proposed models provides a model-free way to select important predictive prognostic markers based on controlling false discovery rates. We compare the performance of our methods with baseline reference survival methods and apply our methodology to an unpublished data set of brain tumor survival times and gene expression data, selecting genes potentially related to the development of the disease under study. A closing discussion compares results obtained by Random Forest and Bayesian ensemble methods under the biological/clinical perspectives and highlights the statistical advantages and disadvantages of the new methodology in the context of DNA microarray data analysis.

Transcriptional regulation of the Borrelia burgdorferi antigenically variable VlsE surface protein.

The Lyme disease agent Borrelia burgdorferi can persistently infect humans and other animals despite host active immune responses. This is facilitated, in part, by the vls locus, a complex system consisting of the vlsE expression site and an adjacent set of 11 to 15 silent vls cassettes. Segments of nonexpressed cassettes recombine with the vlsE region during infection of mammalian hosts, resulting in combinatorial antigenic variation of the VlsE outer surface protein. We now demonstrate that synthesis of VlsE is regulated during the natural mammal-tick infectious cycle, being activated in mammals but repressed during tick colonization. Examination of cultured B. burgdorferi cells indicated that the spirochete controls vlsE transcription levels in response to environmental cues. Analysis of PvlsE::gfp fusions in B. burgdorferi indicated that VlsE production is controlled at the level of transcriptional initiation, and regions of 5' DNA involved in the regulation were identified. Electrophoretic mobility shift assays detected qualitative and quantitative changes in patterns of protein-DNA complexes formed between the vlsE promoter and cytoplasmic proteins, suggesting the involvement of DNA-binding proteins in the regulation of vlsE, with at least one protein acting as a transcriptional activator.

Three-dimensional printing creates models for surgical planning of aortic valve replacement after previous coronary bypass grafting.

PURPOSE Resternotomy for aortic valve replacement in patients with previous coronary artery bypass grafting and an internal mammary artery graft may be a surgical problem. Thus, we are exploring the effect of using rapid prototyping techniques for surgical planning and intraoperative orientation during aortic valve replacement after previous coronary artery bypass grafting (CABG). DESCRIPTION As a proof of concept, we studied a patient who had undergone CABG 5 years earlier. At that time the patient received a left internal mammary artery graft to the left anterior descending artery and a venous graft to the right coronary artery. Now the patient required aortic valve replacement due to symptomatic aortic valve stenosis. The left internal mammary artery bypass and the right coronary artery bypass were patent and showed good flow in the angiography. The patient was examined by 128-slice computed tomography. The image data were visualized and reconstructed. Afterwards, a replica showing the anatomic structures was fabricated using a rapid prototyping machine. EVALUATION Using data derived from 128-slice computed tomography angiography linked to proprietary software, we were able to create three-dimensional reconstructions of the vascular anatomy after the previous CABG. The models were sterilized and taken to the operating theatre for orientation during the surgical procedure. CONCLUSIONS Stereolithographic replicas are helpful for choosing treatment strategies in surgical planning and for intraoperative orientation during reoperations of patients with previous CABG.

Statistical characterization and development of summary measures of non-normal distributions of nuclear morphometry data from prostate cancer cells for use as prognostic indicators

Nuclear morphometry (NM) uses image analysis to measure features of the cell nucleus which are classified as: bulk properties, shape or form, and DNA distribution. Studies have used these measurements as diagnostic and prognostic indicators of disease with inconclusive results. The distributional properties of these variables have not been systematically investigated although much of the medical data exhibit nonnormal distributions. Measurements are done on several hundred cells per patient so summary measurements reflecting the underlying distribution are needed.^ Distributional characteristics of 34 NM variables from prostate cancer cells were investigated using graphical and analytical techniques. Cells per sample ranged from 52 to 458. A small sample of patients with benign prostatic hyperplasia (BPH), representing non-cancer cells, was used for general comparison with the cancer cells.^ Data transformations such as log, square root and 1/x did not yield normality as measured by the Shapiro-Wilks test for normality. A modulus transformation, used for distributions having abnormal kurtosis values, also did not produce normality.^ Kernel density histograms of the 34 variables exhibited non-normality and 18 variables also exhibited bimodality. A bimodality coefficient was calculated and 3 variables: DNA concentration, shape and elongation, showed the strongest evidence of bimodality and were studied further.^ Two analytical approaches were used to obtain a summary measure for each variable for each patient: cluster analysis to determine significant clusters and a mixture model analysis using a two component model having a Gaussian distribution with equal variances. The mixture component parameters were used to bootstrap the log likelihood ratio to determine the significant number of components, 1 or 2. These summary measures were used as predictors of disease severity in several proportional odds logistic regression models. The disease severity scale had 5 levels and was constructed of 3 components: extracapsulary penetration (ECP), lymph node involvement (LN+) and seminal vesicle involvement (SV+) which represent surrogate measures of prognosis. The summary measures were not strong predictors of disease severity. There was some indication from the mixture model results that there were changes in mean levels and proportions of the components in the lower severity levels. ^

Theoretical Calving Rates from Glaciers Along Ice Walls Grounded in Water of Variable Depths

Calving has been studied for glaciers ranging from slow polar glaciers that calve on dry land, such as on Deception Island (63.0-degrees-S, 60.6-degrees-W) in Antarctica, through temperate Alaskan tide-water glaciers, to fast outlet glaciers that float in fiords and calve in deep water, such as Jakobshavns Isbrae (69.2-degrees-N, 49.9-degrees-W) in Greenland. Calving from grounded ice walls and floating ice shelves is the main ablation mechanism for the Antarctic and Greenland ice sheets, as it was along marine and lacustrine margins of former Pleistocene ice sheets, and is for tide-water and polar glaciers. Yet, the theory of ice calving is underdeveloped because of inherent dangers in obtaining field data to test and constrain calving models. An attempt is made to develop a calving theory for ice walls grounded in water of variable depth, and to relate slab calving from ice walls to tabular calving from ice shelves. A calving law is derived in which calving rates from ice walls are controled by bending creep behind the ice wall, and depend on wall height h, forward bending angle-theta, crevasse distance c behind the ice wall and depth d of water in front of the ice wall. Reasonable agreement with calving rates reported by Brown and others (1982) for Alaskan tide-water glaciers is obtained when c depends on wall height, wall height above water and water depth. More data are needed to determine which of these dependencies is correct. A calving ratio c/h is introduced to understand the transition from slab calving to tabular calving as water deepens and the calving glacier becomes afloat.

An application of combinatorial optimization in the printing industry

Offset printing is a common method to produce large amounts of printed matter. We consider a real-world offset printing process that is used to imprint customer-specific designs on napkin pouches. The print- ing technology used yields a number of specific constraints. The planning problem consists of allocating designs to printing-plate slots such that the given customer demand for each design is fulfilled, all technologi- cal and organizational constraints are met and the total overproduction and setup costs are minimized. We formulate this planning problem as a mixed-binary linear program, and we develop a multi-pass matching-based savings heuristic. We report computational results for a set of problem instances devised from real-world data.

Reliability of operational data from pig herds and performance ratings by veterinarians and pig farmers collected during telephone interviews for the evaluation of a PCV2 piglet vaccination

BackgroundThe aim of the present study was to evaluate the feasibility of using a telephone survey in gaining an understanding of the possible herd and management factors influencing the performance (i.e. safety and efficacy) of a vaccine against porcine circovirus type 2 (PCV2) in a large number of herds and to estimate customers¿ satisfaction.ResultsDatasets from 227 pig herds that currently applied or have applied a PCV2 vaccine were analysed. Since 1-, 2- and 3-site production systems were surveyed, the herds were allocated in one of two subsets, where only applicable variables out of 180 were analysed. Group 1 was comprised of herds with sows, suckling pigs and nursery pigs, whereas herds in Group 2 in all cases kept fattening pigs. Overall 14 variables evaluating the subjective satisfaction with one particular PCV2 vaccine were comingled to an abstract dependent variable for further models, which was characterized by a binary outcome from a cluster analysis: good/excellent satisfaction (green cluster) and moderate satisfaction (red cluster). The other 166 variables comprised information about diagnostics, vaccination, housing, management, were considered as independent variables. In Group 1, herds using the vaccine due to recognised PCV2 related health problems (wasting, mortality or porcine dermatitis and nephropathy syndrome) had a 2.4-fold increased chance (1/OR) of belonging to the green cluster. In the final model for Group 1, the diagnosis of diseases other than PCV2, the reason for vaccine administration being other than PCV2-associated diseases and using a single injection of iron had significant influence on allocating into the green cluster (P¿<¿0.05). In Group 2, only unchanged time or delay of time of vaccination influenced the satisfaction (P¿<¿0.05).ConclusionThe methodology and statistical approach used in this study were feasible to scientifically assess ¿satisfaction¿, and to determine factors influencing farmers¿ and vets¿ opinion about the safety and efficacy of a new vaccine.

Lake surface water temperatures of European Alpine lakes (1989–2013) based on the Advanced Very High Resolution Radiometer (AVHRR) 1 km data set

Lake water temperature (LWT) is an important driver of lake ecosystems and it has been identified as an indicator of climate change. Consequently, the Global Climate Observing System (GCOS) lists LWT as an essential climate variable. Although for some European lakes long in situ time series of LWT do exist, many lakes are not observed or only on a non-regular basis making these observations insufficient for climate monitoring. Satellite data can provide the information needed. However, only few satellite sensors offer the possibility to analyse time series which cover 25 years or more. The Advanced Very High Resolution Radiometer (AVHRR) is among these and has been flown as a heritage instrument for almost 35 years. It will be carried on for at least ten more years, offering a unique opportunity for satellite-based climate studies. Herein we present a satellite-based lake surface water temperature (LSWT) data set for European water bodies in or near the Alps based on the extensive AVHRR 1 km data record (1989–2013) of the Remote Sensing Research Group at the University of Bern. It has been compiled out of AVHRR/2 (NOAA-07, -09, -11, -14) and AVHRR/3 (NOAA-16, -17, -18, -19 and MetOp-A) data. The high accuracy needed for climate related studies requires careful pre-processing and consideration of the atmospheric state. The LSWT retrieval is based on a simulation-based scheme making use of the Radiative Transfer for TOVS (RTTOV) Version 10 together with ERA-interim reanalysis data from the European Centre for Medium-range Weather Forecasts. The resulting LSWTs were extensively compared with in situ measurements from lakes with various sizes between 14 and 580 km2 and the resulting biases and RMSEs were found to be within the range of −0.5 to 0.6 K and 1.0 to 1.6 K, respectively. The upper limits of the reported errors could be rather attributed to uncertainties in the data comparison between in situ and satellite observations than inaccuracies of the satellite retrieval. An inter-comparison with the standard Moderate-resolution Imaging Spectroradiometer (MODIS) Land Surface Temperature product exhibits RMSEs and biases in the range of 0.6 to 0.9 and −0.5 to 0.2 K, respectively. The cross-platform consistency of the retrieval was found to be within ~ 0.3 K. For one lake, the satellite-derived trend was compared with the trend of in situ measurements and both were found to be similar. Thus, orbital drift is not causing artificial temperature trends in the data set. A comparison with LSWT derived through global sea surface temperature (SST) algorithms shows lower RMSEs and biases for the simulation-based approach. A running project will apply the developed method to retrieve LSWT for all of Europe to derive the climate signal of the last 30 years. The data are available at doi:10.1594/PANGAEA.831007.

Planning of a make-to-order production process in the printing industry

Offset printing is a common method to produce large amounts of printed matter. We consider a real-world offset printing process that is used to imprint customer-specific designs on napkin pouches. The production equipment used gives rise to various technological constraints. The planning problem consists of allocating designs to printing-plate slots such that the given customer demand for each design is fulfilled, all technological and organizational constraints are met and the total overproduction and setup costs are minimized. We formulate this planning problem as a mixed-binary linear program, and we develop a multi-pass matching-based savings heuristic. We report computational results for a set of problem instances devised from real-world data.

Three-dimensional printing for perioperative planning of complex aortic arch surgery

Purpose In this study, we show the use of three-dimensional printing models for preoperative planning of surgery for patients with complex aortic arch anomalies. Description A 70-year-old man with an extensively arteriosclerotic aneurysm reaching from the ascending aorta to the descending aorta was referred to our center for complete aortic arch replacement. We visualized and reconstructed computed tomography data of the patient and fabricated a flexible three-dimensional model of the aortic arch including the aneurysm. Evaluation This model was very helpful for the preoperative decision making and planning of the frozen elephant trunk procedure owing to the exact and lifelike illustration of the native aortic arch. Conclusions Three-dimensional models are helpful in preoperative planning and postoperative evaluation of frozen elephant trunk procedures in patients with complex aortic anatomy.

Automated 3D Lumbar Intervertebral Disc Segmentation from MRI Data Sets

This paper proposed an automated 3D lumbar intervertebral disc (IVD) segmentation strategy from MRI data. Starting from two user supplied landmarks, the geometrical parameters of all lumbar vertebral bodies and intervertebral discs are automatically extracted from a mid-sagittal slice using a graphical model based approach. After that, a three-dimensional (3D) variable-radius soft tube model of the lumbar spine column is built to guide the 3D disc segmentation. The disc segmentation is achieved as a multi-kernel diffeomorphic registration between a 3D template of the disc and the observed MRI data. Experiments on 15 patient data sets showed the robustness and the accuracy of the proposed algorithm.

Quick divergence but slow convergence during ecotype formation in lake and stream stickleback pairs of variable age

When genetic constraints restrict phenotypic evolution, diversification can be predicted to evolve along so-called lines of least resistance. To address the importance of such constraints and their resolution, studies of parallel phenotypic divergence that differ in their age are valuable. Here, we investigate the parapatric evolution of six lake and stream threespine stickleback systems from Iceland and Switzerland, ranging in age from a few decades to several millennia. Using phenotypic data, we test for parallelism in ecotypic divergence between parapatric lake and stream populations and compare the observed patterns to an ancestral-like marine population. We find strong and consistent phenotypic divergence, both among lake and stream populations and between our freshwater populations and the marine population. Interestingly, ecotypic divergence in low-dimensional phenotype space (i.e. single traits) is rapid and seems to be often completed within 100 years. Yet, the dimensionality of ecotypic divergence was highest in our oldest systems and only there parallel evolution of unrelated ecotypes was strong enough to overwrite phylogenetic contingency. Moreover, the dimensionality of divergence in different systems varies between trait complexes, suggesting different constraints and evolutionary pathways to their resolution among freshwater systems.

Eclogitic diamonds from variable crustal protoliths in the northeastern Siberian craton: trace elements and coupled δ13C – δ18O signatures in diamonds and garnet inclusions

Diamonds of eclogitic assemblages are dominant in the placer diamond deposits of the northeastern Siberian platform. In this study we present new trace elements and stable isotopes (δ13C and δ18O) data for alluvial diamonds and their garnet inclusions from this locality. Cr-rich garnets of peridotitic affinity in the studied diamonds have a narrow range of δ18O values from 5.7‰ to 6.2‰, which is largely overlapping with the accepted mantle range. This narrow range suggests that the garnet inclusions showing different REE patterns and little variations in oxygen isotopes may have formed by different processes involving fluid/melts that, however, were in oxygen isotopic equilibrium with the mantle. The trace element composition of the eclogitic garnet inclusions supports a crustal origin for at least the high-Ca garnets, which show flat HREE patterns and in some cases a positive Eu-anomaly. High-Ca eclogitic garnets generally show heavier oxygen isotope compositions (δ18O 6.5–9.6‰) than what is observed in low-Ca garnets (δ18O 5.7–7.4‰). The variability in oxygen isotopes and trace elements is suggested to be inherited from contrasting crustal protoliths. The relationship between the high δ18O values of inclusions and the low δ13C values of the host diamonds implies that the high-Ca garnet inclusions were derived from intensely hydrated (e.g., δ18O > 7‰) and typically oxidised basaltic rock close to the seawater interface, and that the carbon for diamonds was closely associated with this protolith.

Automated 3D Lumbar Intervertebral Disc Segmentation from MRI Data Sets

This paper proposed an automated three-dimensional (3D) lumbar intervertebral disc (IVD) segmentation strategy from Magnetic Resonance Imaging (MRI) data. Starting from two user supplied landmarks, the geometrical parameters of all lumbar vertebral bodies and intervertebral discs are automatically extracted from a mid-sagittal slice using a graphical model based template matching approach. Based on the estimated two-dimensional (2D) geometrical parameters, a 3D variable-radius soft tube model of the lumbar spine column is built by model fitting to the 3D data volume. Taking the geometrical information from the 3D lumbar spine column as constraints and segmentation initialization, the disc segmentation is achieved by a multi-kernel diffeomorphic registration between a 3D template of the disc and the observed MRI data. Experiments on 15 patient data sets showed the robustness and the accuracy of the proposed algorithm.