Tree phenology and carbon dioxide fluxes - use of digital photography for process-based interpretation at the ecosystem scale

Vegetation phenology is an important indicator of climate change and climate variability and it is strongly connected to biospheric–atmospheric gas exchange. We aimed to evaluate the applicability of phenological information derived from digital imagery for the interpretation of CO2 exchange measurements. For the years 2005–2007 we analyzed seasonal phenological development of 2 temperate mixed forests using tower-based imagery from standard RGB cameras. Phenological information was jointly analyzed with gross primary productivity (GPP) derived from net ecosystem exchange data. Automated image analysis provided reliable information on vegetation developmental stages of beech and ash trees covering all seasons. A phenological index derived from image color values was strongly correlated with GPP, with a significant mean time lag of several days for ash trees and several weeks for beech trees in early summer (May to mid-July). Leaf emergence dates for the dominant tree species partly explained temporal behaviour of spring GPP but were also masked by local meteorological conditions. We conclude that digital cameras at flux measurement sites not only provide an objective measure of the physiological state of a forest canopy at high temporal and spatial resolutions, but also complement CO2 and water exchange measurements, improving our knowledge of ecosystem processes.

Topologically Accurate Dual Isosurfacing Using Ray Intersection

“Dual contouring” approaches provide an alternative to standard Marching Cubes (MC) method to extract and approximate an isosurface from trivariate data given on a volumetric mesh. These dual approaches solve some of the problems encountered by the MC methods. We present a simple method based on the MC method and the ray intersection technique to compute isosurface points in the cell interior. One of the advantages of our method is that it does not require us to use Hermite interpolation scheme, unlike other dual contouring methods. We perform a complete analysis of all possible configurations to generate a look-up table for all configurations. We use the look-up table to optimize the ray-intersection method to obtain minimum number of points necessarily sufficient for defining topologically correct isosurfaces in all possible configurations. Isosurface points are connected using a simple strategy.

Estimation of tool pose based on force-density correlation during robotic drilling

The application of image-guided systems with or without support by surgical robots relies on the accuracy of the navigation process, including patient-to-image registration. The surgeon must carry out the procedure based on the information provided by the navigation system, usually without being able to verify its correctness beyond visual inspection. Misleading surrogate parameters such as the fiducial registration error are often used to describe the success of the registration process, while a lack of methods describing the effects of navigation errors, such as those caused by tracking or calibration, may prevent the application of image guidance in certain accuracy-critical interventions. During minimally invasive mastoidectomy for cochlear implantation, a direct tunnel is drilled from the outside of the mastoid to a target on the cochlea based on registration using landmarks solely on the surface of the skull. Using this methodology, it is impossible to detect if the drill is advancing in the correct direction and that injury of the facial nerve will be avoided. To overcome this problem, a tool localization method based on drilling process information is proposed. The algorithm estimates the pose of a robot-guided surgical tool during a drilling task based on the correlation of the observed axial drilling force and the heterogeneous bone density in the mastoid extracted from 3-D image data. We present here one possible implementation of this method tested on ten tunnels drilled into three human cadaver specimens where an average tool localization accuracy of 0.29 mm was observed.

Implementation of a virtual correlative light and transmission electron microscope

In the long run, the widespread use of slide scanners by pathologists requires an adaptation of teaching methods in histology and cytology in order to target these new possibilities of image processing and presentation via the internet. Accordingly, we were looking for a tool with the possibility to teach microscopic anatomy, histology, and cytology of tissue samples which would be able to combine image data from light and electron microscopes independently of microscope suppliers. With the example of a section through the villus of jejunum, we describe here how to process image data from light and electron microscopes in order to get one image-stack which allows a correlation of structures from the microscopic anatomic to the cytological level. With commercially available image-presentation software that we adapted to our needs, we present here a platform which allows for the presentation of this new but also of older material independently of microscope suppliers.

Quantitative 3D strain analysis in analogue experiments: Integration of X-ray computed tomography and digital volume correlation techniques

The combination of scaled analogue experiments, material mechanics, X-ray computed tomography (XRCT) and Digital Volume Correlation techniques (DVC) is a powerful new tool not only to examine the 3 dimensional structure and kinematic evolution of complex deformation structures in scaled analogue experiments, but also to fully quantify their spatial strain distribution and complete strain history. Digital image correlation (DIC) is an important advance in quantitative physical modelling and helps to understand non-linear deformation processes. Optical non-intrusive (DIC) techniques enable the quantification of localised and distributed deformation in analogue experiments based either on images taken through transparent sidewalls (2D DIC) or on surface views (3D DIC). X-ray computed tomography (XRCT) analysis permits the non-destructive visualisation of the internal structure and kinematic evolution of scaled analogue experiments simulating tectonic evolution of complex geological structures. The combination of XRCT sectional image data of analogue experiments with 2D DIC only allows quantification of 2D displacement and strain components in section direction. This completely omits the potential of CT experiments for full 3D strain analysis of complex, non-cylindrical deformation structures. In this study, we apply digital volume correlation (DVC) techniques on XRCT scan data of “solid” analogue experiments to fully quantify the internal displacement and strain in 3 dimensions over time. Our first results indicate that the application of DVC techniques on XRCT volume data can successfully be used to quantify the 3D spatial and temporal strain patterns inside analogue experiments. We demonstrate the potential of combining DVC techniques and XRCT volume imaging for 3D strain analysis of a contractional experiment simulating the development of a non-cylindrical pop-up structure. Furthermore, we discuss various options for optimisation of granular materials, pattern generation, and data acquisition for increased resolution and accuracy of the strain results. Three-dimensional strain analysis of analogue models is of particular interest for geological and seismic interpretations of complex, non-cylindrical geological structures. The volume strain data enable the analysis of the large-scale and small-scale strain history of geological structures.

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.

Divergent selection during speciation of Lake Malawi cichlid fishes inferred from parallel radiations in nuptial coloration

Repeated evolution of the same phenotypic difference during independent episodes of speciation is strong evidence for selection during speciation. More than 1,000 species of cichlids, >10% of the world's freshwater fish species, have arisen within the past million years in Lakes Malawi and Victoria in eastern Africa. Many pairs of closely related sympatric species differ in their nuptial coloration in very similar ways. Nuptial coloration is important in their mate choice, and speciation by sexual selection on genetically or ecologically constrained variation in nuptial coloration had been proposed, which would repeatedly produce similar nuptial types in different populations, a prediction that was difficult to test in the absence of population-level phylogenies. We measured genetic similarity between individuals within and between populations, species, and lake regions by typing 59 individuals at >2,000 polymorphic genetic loci. From these data, we reconstructed, to our knowledge, the first larger species level phylogeny for the most diverse group of Lake Malawi cichlids. We used the genetic and phylogenetic data to test the divergent selection scenario against colonization, character displacement, and hybridization scenarios that could also explain diverse communities. Diversity has arisen by replicated radiations into the same color types, resulting in phenotypically very different, yet closely related, species within and phenotypically highly similar yet unrelated sets of species between regions, which is consistent with divergent selection during speciation and is inconsistent with colonization and character displacement models.

Scanning electron microscope study of the developing microvasculature in the postnatal rat lung.

During postnatal growth the parenchymal septa of rat lung undergo an impressive restructuring. While immature septa are thick and contain two capillary layers, mature septa are slender and contain a single microvascular network. Using the Mercox casting technique and scanning electron microscopy, we investigated the mode and the timing of the transformation of the pulmonary capillary bed. During the third postnatal week the parenchymal septa rapidly mature to match adult morphology. Even in adult lungs, however, remnants of the immature status are present: A capillary bilayer is regularly found at the base and the tip of the septa. Our observations support the concept that reduction of intervening tissue, partial fusion of the two capillary networks, and preferential growth lead to the mature vascular arrangement. The fact that true mature interalveolar septa show a denser capillary network than alveolar walls abutting onto pleura, bronchi, or larger vessels is consonant with the fusion theory. Towards the nonparenchyma, the capillary network surrounding every airspace had no counterpart to fuse with. From quantitative data it can be calculated that owing to lung growth, mesh size should increase more than four times between birth and adult age. The adult lung network, however, is denser than the one in young animals. This means that new meshes must be added during growth. We propose that small holes observed in sheet-like regions of the microvasculature enlarge to form new capillary meshes. With this mechanism of in-itself or intussusceptional growth, sprouting of individual capillary segments to increase network size is no longer needed.

Sulfur Isotopic Measurements from a West Antarctic Ice Core: Implications for Sulfate Source and Transport

Measurements of delta(34)S covering the years 1935-76 and including the 1963 Agung (Indonesia) eruption were made on a West Antarctic firn core, RIDSA (78.73 degrees S, 116.33 degrees W; 1740m a.s.l.), and results are used to unravel potential source functions in the sulfur cycle over West Antarctica. The delta(34)S values Of SO42- range from 3.1 parts per thousand to 9.9 parts per thousand. These values are lower than those reported for central Antarctica, from near South Pole station, of 9.3-18.1 parts per thousand (Patris and others, 2000). While the Agung period is isotopically distinct at South Pole, it is not in the RIDSA dataset, suggesting differences in the source associations for the sulfur cycle between these two regions. Given the relatively large input of marine aerosols at RIDSA (determined from Na+ data and the seasonal SO42- cycle), there is likely a large marine biogenic SO42- influence. The delta(34)S values indicate, however, that this marine biogenic SO42-, with a well-established delta(34)S of 18 parts per thousand, is mixing with SO42- that has extremely negative delta(34)S values to produce the measured isotope values in the RIDSA core. We suggest that the transport and deposition of stratospheric SO42- in West Antarctica, combined with local volcanic input, accounts for the observed variance in delta(34)S values.

Atmospheric Dynamics of Hot Exoplanets

The characterization of exoplanetary atmospheres has come of age in the last decade, as astronomical techniques now allow for albedos, chemical abundances, temperature profiles and maps, rotation periods and even wind speeds to be measured. Atmospheric dynamics sets the background state of density, temperature and velocity that determines or influences the spectral and temporal appearance of an exoplanetary atmosphere. Hot exoplanets are most amenable to these characterization techniques; in the present review, we focus on highly-irradiated, large exoplanets (the "hot Jupiters"), as astronomical data begin to confront theoretical questions. We summarize the basic atmospheric quantities inferred from the astronomical observations. We review the state of the art by addressing a series of current questions and look towards the future by considering a separate set of exploratory questions. Attaining the next level of understanding will require a concerted effort of constructing multi-faceted, multi-wavelength datasets for benchmark objects. Understanding clouds presents a formidable obstacle, as they introduce degeneracies into the interpretation of spectra, yet their properties and existence are directly influenced by atmospheric dynamics. Confronting general circulation models with these multi-faceted, multi-wavelength datasets will help us understand these and other degeneracies. The coming decade will witness a decisive confrontation of theory and simulation by the next generation of astronomical data.

Patient-Specific Semi-Supervised Learning for Postoperative Brain Tumor Segmentation

In contrast to preoperative brain tumor segmentation, the problem of postoperative brain tumor segmentation has been rarely approached so far. We present a fully-automatic segmentation method using multimodal magnetic resonance image data and patient-specific semi-supervised learning. The idea behind our semi-supervised approach is to effectively fuse information from both pre- and postoperative image data of the same patient to improve segmentation of the postoperative image. We pose image segmentation as a classification problem and solve it by adopting a semi-supervised decision forest. The method is evaluated on a cohort of 10 high-grade glioma patients, with segmentation performance and computation time comparable or superior to a state-of-the-art brain tumor segmentation method. Moreover, our results confirm that the inclusion of preoperative MR images lead to a better performance regarding postoperative brain tumor segmentation.

The diurnal variation in stratospheric ozone from the MACC reanalysis, the ERA-Interim reanalysis, WACCM and Earth observations: characteristics and intercomparison

In this study we compare the diurnal variation in stratospheric ozone derived from free-running simulations of the Whole Atmosphere Community Climate Model (WACCM) and from reanalysis data of the atmospheric service MACC (Monitoring Atmospheric Composition and Climate) which both use a similar stratospheric chemistry module. We find good agreement between WACCM and the MACC reanalysis for the diurnal ozone variation in the high-latitude summer stratosphere based on photochemistry. In addition, we consult the ozone data product of the ERA-Interim reanalysis. The ERA-Interim reanalysis ozone system with its long-term ozone parametrization can not capture these diurnal variations in the upper stratosphere that are due to photochemistry. The good dynamics representations, however, reflects well dynamically induced ozone variations in the lower stratosphere. For the high-latitude winter stratosphere we describe a novel feature of diurnal variation in ozone where changes of up to 46.6% (3.3 ppmv) occur in monthly mean data. For this effect good agreement between the ERA-Interim reanalysis and the MACC reanalysis suggest quite similar diurnal advection processes of ozone. The free-running WACCM model seriously underestimates the role of diurnal advection processes at the polar vortex at the two tested resolutions. The intercomparison of the MACC reanalysis and the ERA-Interim reanalysis demonstrates how global reanalyses can benefit from a chemical representation held by a chemical transport model. The MACC reanalysis provides an unprecedented description of the dynamics and photochemistry of the diurnal variation of stratospheric ozone which is of high interest for ozone trend analysis and research on atmospheric tides. We confirm the diurnal variation in ozone at 5 hPa by observations of the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) experiment and selected sites of the Network for Detection of Atmospheric Composition Change (NDACC). The latter give valuable insight even to diurnal variation of ozone in the polar winter stratosphere.

Fully Automatic CT Segmentation for Computer-Assisted Pre-operative Planning of Hip Arthroscopy

Extraction of both pelvic and femoral surface models of a hip joint from CT data for computer-assisted pre-operative planning of hip arthroscopy is addressed. We present a method for a fully automatic image segmentation of a hip joint. Our method works by combining fast random forest (RF) regression based landmark detection, atlas-based segmentation, with articulated statistical shape model (aSSM) based hip joint reconstruction. The two fundamental contributions of our method are: (1) An improved fast Gaussian transform (IFGT) is used within the RF regression framework for a fast and accurate landmark detection, which then allows for a fully automatic initialization of the atlas-based segmentation; and (2) aSSM based fitting is used to preserve hip joint structure and to avoid penetration between the pelvic and femoral models. Validation on 30 hip CT images show that our method achieves high performance in segmenting pelvis, left proximal femur, and right proximal femur surfaces with an average accuracy of 0.59 mm, 0.62 mm, and 0.58 mm, respectively.

The mass-hierarchy and CP-violation discovery reach of the LBNO long-baseline neutrino experiment

The next generation neutrino observatory proposed by the LBNO collaboration will address fundamental questions in particle and astroparticle physics. The experiment consists of a far detector, in its first stage a 20 kt LAr double phase TPC and a magnetised iron calorimeter, situated at 2300 km from CERN and a near detector based on a highpressure argon gas TPC. The long baseline provides a unique opportunity to study neutrino flavour oscillations over their 1st and 2nd oscillation maxima exploring the L/E behaviour, and distinguishing effects arising from δCP and matter. In this paper we have reevaluated the physics potential of this setup for determining the mass hierarchy (MH) and discovering CP-violation (CPV), using a conventional neutrino beam from the CERN SPS with a power of 750 kW. We use conservative assumptions on the knowledge of oscillation parameter priors and systematic uncertainties. The impact of each systematic error and the precision of oscillation prior is shown. We demonstrate that the first stage of LBNO can determine unambiguously the MH to > 5δ C.L. over the whole phase space. We show that the statistical treatment of the experiment is of very high importance, resulting in the conclusion that LBNO has ~ 100% probability to determine the MH in at most 4-5 years of running. Since the knowledge of MH is indispensable to extract δCP from the data, the first LBNO phase can convincingly give evidence for CPV on the 3δ C.L. using today’s knowledge on oscillation parameters and realistic assumptions on the systematic uncertainties.

Ce(III) and Ce(IV) (re)distribution and fractionation in a laterite profile from Madagascar: Insights from in situ XANES spectroscopy at the Ce LIII-edge

The distribution of trivalent and tetravalent cerium, Ce(III) and Ce(IV) respectively, in a lateritic profile from Madagascar, has been characterized by X-ray-absorption near-edge structure (XANES) spectroscopy at the Ce LIII-edge on the LUCIA beamline (SOLEIL synchrotron, France). XANES spectra were acquired on bulk-rock samples as well as on specific lateritic minerals or polymineral zones (in-situ measurements) of the tonalite bedrock and the three overlying weathered horizons (C-, B- and A-horizons). Geochemically, the bedrock, and the A- and C-horizons show similar rare earth element content (REE = 363–405 mg/kg). They also display the same positive Ce-anomaly (CeCN/Ce∗ = 1.12–1.45), which is therefore likely to be inherited from the bedrock. In the B-horizon, the higher REE content (REE = 2194 mg/kg) and the larger Ce-anomaly (CeCN/Ce∗ = 4.26) are consistent with an accumulation zone caused by the evaporation of groundwater during the dry season. There is a good agreement between the Ce(III)/Cetotal ratio (XCe(III)) deduced from the positive Ce-anomaly (bulk-rock geochemical data) and that derived from XANES spectroscopy on the same bulk-rock samples (BR-XCe(III)-XANES) in the bedrock, and the C- and B-horizons. In the A-horizon, XANES measurements on bulk rock and minerals revealed a higher BR-XCe(III)-XANES (up to 100%) compared to the XCe(III) deduced from geochemical data (XCe(III) = 79%). The preservation of a positive Ce-anomaly in the A-horizon suggests that the Ce mobilization and redistribution during weathering occurred with no significant Ce fractionation from other trivalent REE. Remarkably, the only investigated sample where cerianite is observed belongs to the B-horizon. Within this horizon, Ce oxidation state varies depending on the microstructural position (porosity, cracks, clay-rich groundmass). The highest Ce(IV) concentrations are measured in cerianite (and aluminophosphates) localized in pores at the vicinity of Mn-rich domains (XCe(III)-XANES = 30–51%). Therefore, Ce fractionation from other REE is attributed to a Ce oxidation and precipitation potentially assisted by oxyhydroxide scavenging. In the C-horizon, Ce(III) and Ce(IV) are mainly distributed in REE-minerals of the rhabdophane group found in pores and cracks. The similarity between the Ce(III) proportion of rhabdophane grains (XCe(III)-XANES = 74–89%) with that of the bedrock (BR-XCe(III)-XANES = 79%) suggests no significant fractionation of Ce(III) and Ce(IV) between solution and mineral during the successive stages of primary REE-mineral alteration, transport in solution and secondary precipitation in the incipient stages of weathering. Overall, our novel spectroscopic approach shows that Ce is not necessarily oxidized nor fractionated from other REE during weathering in lateritic conditions. This implies that like Ce(III), Ce(IV) can be mobilized in aqueous fluids during weathering, possibly thanks to complexation with organic molecules, and can precipitate together with Ce(III) in secondary REE-bearing minerals. The corollary is that (paleo)redox reconstructions in soils and/or sediments based on Ce-anomaly in weathered rocks or minerals must be interpreted with caution.