Fast parallel kriging-based stepwise uncertainty reduction with application to the identification of an excursion set

Stepwise uncertainty reduction (SUR) strategies aim at constructing a sequence of points for evaluating a function  f in such a way that the residual uncertainty about a quantity of interest progressively decreases to zero. Using such strategies in the framework of Gaussian process modeling has been shown to be efficient for estimating the volume of excursion of f above a fixed threshold. However, SUR strategies remain cumbersome to use in practice because of their high computational complexity, and the fact that they deliver a single point at each iteration. In this article we introduce several multipoint sampling criteria, allowing the selection of batches of points at which f can be evaluated in parallel. Such criteria are of particular interest when f is costly to evaluate and several CPUs are simultaneously available. We also manage to drastically reduce the computational cost of these strategies through the use of closed form formulas. We illustrate their performances in various numerical experiments, including a nuclear safety test case. Basic notions about kriging, auxiliary problems, complexity calculations, R code, and data are available online as supplementary materials.

Radiological anatomy of the obturator nerve and its articular branches: basis to develop a method of radiofrequency denervation for hip joint pain

OBJECTIVE: A previous study of radiofrequency neurotomy of the articular branches of the obturator nerve for hip joint pain produced modest results. Based on an anatomical and radiological study, we sought to define a potentially more effective radiofrequency method. DESIGN: Ten cadavers were studied, four of them bilaterally. The obturator nerve and its articular branches were marked by wires. Their radiological relationship to the bone structures on fluoroscopy was imaged and analyzed. A magnetic resonance imaging (MRI) study was undertaken on 20 patients to determine the structures that would be encountered by the radiofrequency electrode during different possible percutaneous approaches. RESULTS: The articular branches of the obturator nerve vary in location over a wide area. The previously described method of denervating the hip joint did not take this variation into account. Moreover, it approached the nerves perpendicularly. Because optimal coagulation requires electrodes to lie parallel to the nerves, a perpendicular approach probably produced only a minimal lesion. In addition, MRI demonstrated that a perpendicular approach is likely to puncture femoral vessels. Vessel puncture can be avoided if an oblique pass is used. Such an approach minimizes the angle between the target nerves and the electrode, and increases the likelihood of the nerve being captured by the lesion made. Multiple lesions need to be made in order to accommodate the variability in location of the articular nerves. CONCLUSIONS: The method that we described has the potential to produce complete and reliable nerve coagulation. Moreover, it minimizes the risk of penetrating the great vessels. The efficacy of this approach should be tested in clinical trials.

Geochemical and isotopic characterization of the Bodélé Depression dust source and implications for transatlantic dust transport to the Amazon Basin

The Bodélé Depression (Chad) in the central Sahara/Sahel region of Northern Africa is the most important source of mineral dust to the atmosphere globally. The Bodélé Depression is purportedly the largest source of Saharan dust reaching the Amazon Basin by transatlantic transport. Here, we have undertaken a comprehensive study of surface sediments from the Bodélé Depression and dust deposits (Chad, Niger) in order to characterize geochemically and isotopically (Sr, Nd and Pb isotopes) this dust source, and evaluate its importance in present and past African dust records. We similarly analyzed sedimentary deposits from the Amazonian lowlands in order to assess postulated accumulation of African mineral dust in the Amazon Basin, as well as its possible impact in fertilizing the Amazon rainforest. Our results identify distinct sources of different ages and provenance in the Bodélé Depression versus the Amazon Basin, effectively ruling out an origin for the Amazonian deposits, such as the Belterra Clay Layer, by long-term deposition of Bodélé Depression material. Similarly, no evidence for contributions from other potential source areas is provided by existing isotope data (Sr, Nd) on Saharan dusts. Instead, the composition of these Amazonian deposits is entirely consistent with derivation from in-situ weathering and erosion of the Precambrian Amazonian craton, with little, if any, Andean contribution. In the Amazon Basin, the mass accumulation rate of eolian dust is only around one-third of the vertical erosion rate in shield areas, suggesting that Saharan dust is “consumed” by tropical weathering, contributing nutrients and stimulating plant growth, but never accumulates as such in the Amazon Basin. The chemical and isotope compositions found in the Bodélé Depression are varied at the local scale, and have contrasting signatures in the “silica-rich” dry lake-bed sediments and in the “calcium-rich” mixed diatomites and surrounding sand material. This unexpected finding implies that the Bodélé Depression material is not “pre-mixed” at the source to provide a homogeneous source of dust. Rather, different isotope signatures can be emitted depending on subtle vagaries of dust-producing events. Our characterization of the Bodélé Depression components indicate that the Bodélé “calcium-rich” component, identified here, is most likely released via eolian processes of sand grain saltation and abrasion and may be significant in the overall global budget of dusts carried out by the Harmattan low-level jet during the winter.

The importance of biotic interactions and local adaptation for plant response to environmental changes: field evidence along an elevational gradient

Predicting the response of species to environmental changes is a great and on-going challenge for ecologists, and this requires a more in-depth understanding of the importance of biotic interactions and the population structuration in the landscape. Using a reciprocal transplantation experiment, we tested the response of five species to an elevational gradient. This was combined to a neighbour removal treatment to test the importance of local adaptation and biotic interactions. The trait studied was performance measured as survival and biomass. Species response varied along the elevational gradient, but with no consistent pattern. Performance of species was influenced by environmental conditions occurring locally at each site, as well as by positive or negative effects of the surrounding vegetation. Indeed, we observed a shift from competition for biomass to facilitation for survival as a response to the increase in environmental stress occurring in the different sites. Unlike previous studies pointing out an increase of stress along the elevation gradient, our results supported a stress gradient related to water availability, which was not strictly parallel to the elevational gradient. For three of our species, we observed a greater biomass production for the population coming from the site where the species was dominant (central population) compared to population sampled at the limit of the distribution (marginal population). Nevertheless, we did not observe any pattern of local adaptation that could indicate adaptation of populations to a particular habitat. Altogether, our results highlighted the great ability of plant species to cope with environmental changes, with no local adaptation and great variability in response to local conditions. Our study confirms the importance of taking into account biotic interactions and population structure occurring at local scale in the prediction of communities’ responses to global environmental changes.

Identification of gunshots to the head by detection of RNA in backspatter primarily expressed in brain tissue

Traces of backspatter recovered from the inside of the barrel of a gun that was used to deliver suicidal or homicidal contact shots may be a source of valuable forensic evidence and first systematic investigations of the persistence of victim DNA from inside firearms have been presented. The aim of the present study was to include victim RNA in such analyses to determine the origin of tissues in addition and parallel to standard DNA profiling for forensic identification purposes. In a first step, suitable mRNA (C1orf61) and micro-RNAs (miR-124a and miR-124*) that are primarily expressed in brain tissue were selected from potential candidates and confirmed using quantitative PCR (qPCR). Secondly, a co-extraction procedure for RNA and DNA was established and brain differentiability of the selected RNAs was demonstrated via qPCR using samples from experimental shots at ballistic models. In a third step, this procedure was successfully applied to analyse samples from real casework comprising eight cases of suicidal contact shots. In this pilot study, we are first to report the possibility of co-extracting mRNA, miRNA and DNA from ballistic trace samples collected from the inside of firearms and we demonstrate that RNA and DNA based analyses can be performed in parallel to produce informative and highly complementary evidence.

Persistence of endocrine disruption in zebrafish (Danio rerio) after discontinued exposure to the androgen 17β-trenbolone.

The aim of the present study was to investigate the effects of the androgenic endocrine disruptor 17β-trenbolone on the sexual development of zebrafish (Danio rerio) with special emphasis on the question of whether adverse outcomes of developmental exposure are reversible or persistent. An exposure scenario including a recovery phase was chosen to assess the potential reversibility of androgenic effects. Zebrafish were exposed to environmentally relevant concentrations of 17β-trenbolone (1 ng/L-30 ng/L) from fertilization until completion of gonad sexual differentiation (60 d posthatch). Thereafter, exposure was either followed by 40 d of recovery in clean water or continued until 100 d posthatch, the age when zebrafish start being able to reproduce. Fish exposed for 100 d to 10 ng/L or 30 ng/L 17β-trenbolone were masculinized at different biological effect levels, as evidenced from a concentration-dependent shift of the sex ratio toward males as well as a significantly increased maturity of testes. Gonad morphological masculinization occurred in parallel with decreased vitellogenin concentrations in both sexes. Changes of brain aromatase (cyp19b) mRNA expression showed no consistent trend with respect to either exposure duration or concentration. Gonad morphological masculinization as well as the decrease of vitellogenin persisted after depuration over 40 d in clean water. This lack of recovery suggests that androgenic effects on sexual development of zebrafish are irreversible.

Dust measurements in the coma of comet 67P/Churyumov-Gerasimenko inbound to the Sun

Critical measurements for understanding accretion and the dust/gas ratio in the solar nebula, where planets were forming 4.5 billion years ago, are being obtained by the GIADA (Grain Impact Analyser and Dust Accumulator) experiment on the European Space Agency's Rosetta spacecraft orbiting comet 67P/Churyumov-Gerasimenko. Between 3.6 and 3.4 astronomical units inbound, GIADA and OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) detected 35 outflowing grains of mass 10(-10) to 10(-7) kilograms, and 48 grains of mass 10(-5) to 10(-2) kilograms, respectively. Combined with gas data from the MIRO (Microwave Instrument for the Rosetta Orbiter) and ROSINA (Rosetta Orbiter Spectrometer for Ion and Neutral Analysis) instruments, we find a dust/gas mass ratio of 4 +/- 2 averaged over the sunlit nucleus surface. A cloud of larger grains also encircles the nucleus in bound orbits from the previous perihelion. The largest orbiting clumps are meter-sized, confirming the dust/gas ratio of 3 inferred at perihelion from models of dust comae and trails.

Low current performance of the Bern medical cyclotron down to the pA range

A medical cyclotron accelerating H- ions to 18 MeV is in operation at the Bern University Hospital (Inselspital). It is the commercial IBA 18/18 cyclotron equipped with a specifically conceived 6 m long external beam line ending in a separate bunker. This feature is unique for a hospital-based facility and makes it possible to conduct routine radioisotope production for PET diagnostics in parallel with multidisciplinary research activities, among which are novel particle detectors, radiation biophysics, radioprotection, radiochemistry and radiopharmacy developments. Several of these activities, such as radiobiology experiments for example, require low current beams down to the pA range, while medical cyclotrons are designed for high current operation above 10 μA. In this paper, we present the first results on the low current performance of a PET medical cyclotron obtained by ion source, radio-frequency and main coil tuning. With this method, stable beam currents down to (1.5+/- 0.5 ) pA were obtained and measured with a high-sensitivity Faraday cup located at the end of the beam transport line.