Patient-specific finite-element simulation of the human cornea: A clinical validation study on cataract surgery

The planning of refractive surgical interventions is a challenging task. Numerical modeling has been proposed as a solution to support surgical intervention and predict the visual acuity, but validation on patient specific intervention is missing. The purpose of this study was to validate the numerical predictions of the post-operative corneal topography induced by the incisions required for cataract surgery. The corneal topography of 13 patients was assessed preoperatively and postoperatively (1-day and 30-day follow-up) with a Pentacam tomography device. The preoperatively acquired geometric corneal topography – anterior, posterior and pachymetry data – was used to build patient-specific finite element models. For each patient, the effects of the cataract incisions were simulated numerically and the resulting corneal surfaces were compared to the clinical postoperative measurements at one day and at 30-days follow up. Results showed that the model was able to reproduce experimental measurements with an error on the surgically induced sphere of 0.38D one day postoperatively and 0.19D 30 days postoperatively. The standard deviation of the surgically induced cylinder was 0.54D at the first postoperative day and 0.38D 30 days postoperatively. The prediction errors in surface elevation and curvature were below the topography measurement device accuracy of ±5μm and ±0.25D after the 30-day follow-up. The results showed that finite element simulations of corneal biomechanics are able to predict post cataract surgery within topography measurement device accuracy. We can conclude that the numerical simulation can become a valuable tool to plan corneal incisions in cataract surgery and other ophthalmosurgical procedures in order to optimize patients' refractive outcome and visual function.

Climatic drivers of hourly to yearly tree radius variations along a 6°C natural warming gradient

Climate affects the timing, rate and dynamics of tree growth, over time scales ranging from seconds to centuries. Monitoring how a tree's stem radius varies over these time scales can provide insight into intra-annual stem dynamics and improve our understanding of climate impacts on tree physiology and growth processes. Here, we quantify the response of radial conifer stem size to environmental fluctuations via a novel assessment of tree circadian cycles. We analyze four years of sub-hourly data collected from 56 larch and spruce trees growing along a natural temperature gradient of ∼6 °C in the central Swiss Alps. During the growing season, tree stem diameters were greatest at mid-morning and smallest in the late evening, reflecting the daily cycle of water uptake and loss. Along the gradient, amplitudes calculated from the stem radius cycle were ∼50% smaller at the upper site (∼2200 m a.s.l.) relative to the lower site (∼800 m a.s.l.). We show changes in precipitation, temperature and cloud cover have a substantial effect on typical growing season diurnal cycles; amplitudes were nine times smaller on rainy days (>10 mm), and daily amplitudes are approximately 40% larger when the mean daily temperature is 15–20 °C than when it is 5–10 °C. We find that over the growing season in the sub-alpine forests, spruce show greater daily stem water movement than larch. However, under projected future warming, larch could experience up to 50% greater stem water use, which may severely affect future growth on already dry sites. Our data further indicate that because of the confounding influences of radial growth and short-term water dynamics on stem size, conventional methodology probably overstates the effect of water-linked meteorological variables (i.e. precipitation and relative humidity) on intra-annual tree growth. We suggest future studies use intra-seasonal measurements of cell development and consider whether climatic factors produce reversible changes in stem diameter. These study design elements may help researchers more accurately quantify and attribute changes in forest productivity in response to future warming.

Reliability and quality of water isotope data collected with a lowbudget rain collector

RATIONALELow-budget rain collectors for water isotope analysis, such as the `ball-in-funnel type collector' (BiFC), are widely used in studies on stable water isotopes of rain. To date, however, an experimental quality assessment of such devices in relation to climatic factors does not exist. METHODSWe used Cavity Ring-Down Spectrometry (CRDS) to quantify the effects of evaporation on the O-18 values of reference water under controlled conditions as a function of the elapsed time between rainfall and collection for isotope analysis, the sample volume and the relative humidity (RH: 31% and 67%; 25 degrees C). The climate chamber conditions were chosen to reflect the warm and dry end of field conditions that favor evaporative enrichment (EE). We also tested the performance of the BiFC in the field, and compared our H-2/O-18 data obtained by isotope ratio mass spectrometry (IRMS) with those from the Swiss National Network for the Observation of Isotopes in the Water Cycle (ISOT). RESULTSThe EE increased with time, with a 1 increase in the O-18 values after 10days (RH: 25%; 25 degrees C; 35mL (corresponding to a 5mm rain event); p <0.001). The sample volume strongly affected the EE (max. value +1.5 parts per thousand for 7mL samples (i.e., 1mm rain events) after 72h at 31% and 67% RH; p <0.001), whereas the relative humidity had no significant effect. Using the BiFC in the field, we obtained very tight relationships of the H-2/O-18 values (r(2) 0.95) for three sites along an elevational gradient, not significantly different from that of the next ISOT station. CONCLUSIONSSince the chosen experimental conditions were extreme compared with the field conditions, it was concluded that the BiFC is a highly reliable and inexpensive collector of rainwater for isotope analysis. Copyright (c) 2014 John Wiley & Sons, Ltd.

The effect of abrupt climatic warming on biogeochemical cycling and N2O emissions in a terrestrial ecosystem

The large, rapid increase in atmospheric N2O concentrations that occurred concurrent with the abrupt warming at the end of the Last Glacial period might have been the result of a reorganization in global biogeochemical cycles. To explore the sensitivity of nitrogen cycling in terrestrial ecosystems to abrupt warming, we combined a scenario of climate and vegetation composition change based on multiproxy data for the Oldest Dryas–Bølling abrupt warming event at Gerzensee, Switzerland, with a biogeochemical model that simulates terrestrial N uptake and release, including N2O emissions. As for many central European sites, the pollen record at the Gerzensee is remarkable for the abundant presence of the symbiotic nitrogen fixer Hippophaë rhamnoides (L.) during the abrupt warming that also marks the beginning of primary succession on immature glacial soils. Here we show that without additional nitrogen fixation, climate change results in a significant increase of N2O emissions of approximately factor 3.4 (from 6.4 ± 1.9 to 21.6 ± 5.9 mg N2O–N m− 2 yr− 1). Each additional 1000 mg m− 2 yr− 1 of nitrogen added to the ecosystem through N-fixation results in additional N2O emissions of 1.6 mg N2O–N m− 2 yr− 1 for the time with maximum H. rhamnoides coverage. Our results suggest that local reactions of emissions to abrupt climate change could have been considerably faster than the overall atmospheric concentration changes observed in polar ice. Nitrogen enrichment of soils due to the presence of symbiotic N-fixers during early primary succession not only facilitates the establishment of vegetation on soils in their initial stage of development, but can also have considerable influence on biogeochemical cycles and the release of reactive nitrogen trace gases to the atmosphere.

Forest dynamics during the transition from the Oldest Dryas to the Bølling–Allerød at Gerzensee—a simulation study

The aim of this study was to explore potential causes and mechanisms for the sequence and temporal pattern of tree taxa, specifically for the shift from shrub-tundra to birch–juniper woodland during and after the transition from the Oldest Dryas to the Bølling–Allerød in the region surrounding the lake Gerzensee in southern Central Europe. We tested the influence of climate, forest dynamics, community dynamics compared to other causes for delays. For this aim temperature reconstructed from a δ18O-record was used as input driving the multi-species forest-landscape model TreeMig. In a stepwise scenario analysis, population dynamics along with pollen production and transport were simulated and compared with pollen-influx data, according to scenarios of different δ18O/temperature sensitivities, different precipitation levels, with/without inter-specific competition, and with/without prescribed arrival of species. In the best-fitting scenarios, the effects on competitive relationships, pollen production, spatial forest structure, albedo, and surface roughness were examined in more detail. The appearance of most taxa in the data could only be explained by the coldest temperature scenario with a sensitivity of 0.3‰/°C, corresponding to an anomaly of − 15 °C. Once the taxa were present, their temporal pattern was shaped by competition. The later arrival of Pinus could not be explained even by the coldest temperatures, and its timing had to be prescribed by first observations in the pollen record. After the arrival into the simulation area, the expansion of Pinus was further influenced by competitors and minor climate oscillations. The rapid change in the simulated species composition went along with a drastic change in forest structure, leaf area, albedo, and surface roughness. Pollen increased only shortly after biomass. Based on our simulations, two alternative potential scenarios for the pollen pattern can be given: either very cold climate suppressed most species in the Oldest Dryas, or they were delayed by soil formation or migration. One taxon, Pinus, was delayed by migration and then additionally hindered by competition. Community dynamics affected the pattern in two ways: potentially by facilitation, i.e. by nitrogen-fixing pioneer species at the onset, whereas the later pattern was clearly shaped by competition. The simulated structural changes illustrate how vegetation on a larger scale could feed back to the climate system. For a better understanding, a more integrated simulation approach covering also the immigration from refugia would be necessary, for this combines climate-driven population dynamics, migration, individual pollen production and transport, soil dynamics, and physiology of individual pollen production.

A model-data comparison of Holocene timberline changes in the Swiss Alps reveals past and future drivers of mountain forest dynamics

Mountain vegetation is strongly affected by temperature and is expected to shift upwards with climate change. Dynamic vegetation models are often used to assess the impact of climate on vegetation and model output can be compared with paleobotanical data as a reality check. Recent paleoecological studies have revealed regional variation in the upward shift of timberlines in the Northern and Central European Alps in response to rapid warming at the Younger Dryas/Preboreal transition ca. 11700years ago, probably caused by a climatic gradient across the Alps. This contrasts with previous studies that successfully simulated the early Holocene afforestation in the (warmer) Central Alps with a chironomid-inferred temperature reconstruction from the (colder) Northern Alps. We use LandClim, a dynamic landscape vegetation model to simulate mountain forests under different temperature, soil and precipitation scenarios around Iffigsee (2065m a.s.l.) a lake in the Northwestern Swiss Alps, and compare the model output with the paleobotanical records. The model clearly overestimates the upward shift of timberline in a climate scenario that applies chironomid-inferred July-temperature anomalies to all months. However, forest establishment at 9800 cal. BP at Iffigsee is successfully simulated with lower moisture availability and monthly temperatures corrected for stronger seasonality during the early Holocene. The model-data comparison reveals a contraction in the realized niche of Abies alba due to the prominent role of anthropogenic disturbance after ca. 5000 cal. BP, which has important implications for species distribution models (SDMs) that rely on equilibrium with climate and niche stability. Under future climate projections, LandClim indicates a rapid upward shift of mountain vegetation belts by ca. 500m and treeline positions of ca. 2500m a.s.l. by the end of this century. Resulting biodiversity losses in the alpine vegetation belt might be mitigated with low-impact pastoralism to preserve species-rich alpine meadows.

High-resolution dynamic computer simulation analysis of the behavior of sample components with pI values outside the pH gradient established by carrier ampholyte CIEF

The behavior of sample components whose pI values are outside the pH gradient established by 101 hypothetical biprotic carrier ampholytes covering a pH 6-8 range was investigated by computer simulation under constant current conditions with concomitant constant electroosmosis toward the cathode. Data obtained with the sample being applied between zones of carrier ampholytes and on the anodic side of the carrier ampholytes were studied and found to evolve into zone structures comprising three regions between anolyte and catholyte. The focusing region with the pH gradient is bracketed by two isotachopheretic zone structures comprising selected sample and carrier components as isotachophoretic zones. The isotachophoretic structures electrophoretically migrate in opposite direction and their lengths increase with time due to the gradual isotachophoretic decay at the pH gradient edges. Due to electroosmosis, however, the overall pattern is being transported toward the cathode. Sample components whose pI values are outside the established pH gradient are demonstrated to form isotachophoretic zones behind the leading cation of the catholyte (components with pI values larger than 8) and the leading anion of the anolyte (components with pI values smaller than 6). Amphoteric compounds with appropriate pI values or nonamphoteric components can act as isotachophoretic spacer compounds between sample compounds or between the leader and the sample with the highest mobility. The simulation data obtained provide for the first time insight into the dynamics of amphoteric sample components that do not focus within the established pH gradient.

Sampling strategies for capillary isoelectric focusing with electroosmotic zone mobilization assessed by high-resolution dynamic computer simulation

The impact of initial sample distribution on separation and focusing of analytes in a pH 3–11 gradient formed by 101 biprotic carrier ampholytes under concomitant electroosmotic displacement was studied by dynamic high-resolution computer simulation. Data obtained with application of the analytes mixed with the carrier ampholytes (as is customarily done), as a short zone within the initial carrier ampholyte zone, sandwiched between zones of carrier ampholytes, or introduced before or after the initial carrier ampholyte zone were compared. With sampling as a short zone within or adjacent to the carrier ampholytes, separation and focusing of analytes is shown to proceed as a cationic, anionic, or mixed process and separation of the analytes is predicted to be much faster than the separation of the carrier components. Thus, after the initial separation, analytes continue to separate and eventually reach their focusing locations. This is different to the double-peak approach to equilibrium that takes place when analytes and carrier ampholytes are applied as a homogenous mixture. Simulation data reveal that sample application between two zones of carrier ampholytes results in the formation of a pH gradient disturbance as the concentration of the carrier ampholytes within the fluid element initially occupied by the sample will be lower compared to the other parts of the gradient. As a consequence thereof, the properties of this region are sample matrix dependent, the pH gradient is flatter, and the region is likely to represent a conductance gap (hot spot). Simulation data suggest that sample placed at the anodic side or at the anodic end of the initial carrier ampholyte zone are the favorable configurations for capillary isoelectric focusing with electroosmotic zone mobilization.

Dynamic high-resolution computer simulation of electrophoretic enantiomer separations with neutral cyclodextrins as chiral selectors.

GENTRANS, a comprehensive one-dimensional dynamic simulator for electrophoretic separations and transport, was extended for handling electrokinetic chiral separations with a neutral ligand. The code can be employed to study the 1:1 interaction of monovalent weak and strong acids and bases with a single monovalent weak or strong acid or base additive, including a neutral cyclodextrin, under real experimental conditions. It is a tool to investigate the dynamics of chiral separations and to provide insight into the buffer systems used in chiral capillary zone electrophoresis (CZE) and chiral isotachophoresis. Analyte stacking across conductivity and buffer additive gradients, changes of additive concentration, buffer component concentration, pH, and conductivity across migrating sample zones and peaks, and the formation and migration of system peaks can thereby be investigated in a hitherto inaccessible way. For model systems with charged weak bases and neutral modified β-cyclodextrins at acidic pH, for which complexation constants, ionic mobilities, and mobilities of selector-analyte complexes have been determined by CZE, simulated and experimentally determined electropherograms and isotachopherograms are shown to be in good agreement. Simulation data reveal that CZE separations of cationic enantiomers performed in phosphate buffers at low pH occur behind a fast cationic migrating system peak that has a small impact on the buffer composition under which enantiomeric separation takes place.

SLA-Driven Simulation of Multi-Tenant Scalable Cloud-Distributed Enterprise Information Systems

Cloud Computing is an enabler for delivering large-scale, distributed enterprise applications with strict requirements in terms of performance. It is often the case that such applications have complex scaling and Service Level Agreement (SLA) management requirements. In this paper we present a simulation approach for validating and comparing SLA-aware scaling policies using the CloudSim simulator, using data from an actual Distributed Enterprise Information System (dEIS). We extend CloudSim with concurrent and multi-tenant task simulation capabilities. We then show how different scaling policies can be used for simulating multiple dEIS applications. We present multiple experiments depicting the impact of VM scaling on both datacenter energy consumption and dEIS performance indicators.

A compilation of Western European terrestrial records 60–8 ka BP: towards an understanding of latitudinal climatic gradients

Terrestrial records of past climatic conditions, such as lake sediments and speleothems, provide data of great importance for understanding environmental changes. However, unlike marine and ice core records, terrestrial palaeodata are often not available in databases or in a format that is easily accessible to the non-specialist. As a consequence, many excellent terrestrial records are unknown to the broader palaeoclimate community and are not included in compilations, comparisons, or modelling exercises. Here we present a compilation of Western European terrestrial palaeo-records covering, entirely or partially, the 60–8-ka INTIMATE time period. The compilation contains 56 natural archives, including lake records, speleothems, ice cores, and terrestrial proxies in marine records. The compilation is limited to include records of high temporal resolution and/or records that provide climate proxies or quantitative reconstructions of environmental parameters, such as temperature or precipitation, and that are of relevance and interest to a broader community. We briefly review the different types of terrestrial archives, their respective proxies, their interpretation and their application for palaeoclimatic reconstructions. We also discuss the importance of independent chronologies and the issue of record synchronization. The aim of this exercise is to provide the wider palaeo-community with a consistent compilation of high-quality terrestrial records, to facilitate model-data comparisons, and to identify key areas of interest for future investigations. We use the compilation to investigate Western European latitudinal climate gradients during the deglacial period and, despite of poorly constrained chronologies for the older records, we summarize the main results obtained from NW and SW European terrestrial records before the LGM.

A stratigraphic framework for abrupt climatic changes during the Last Glacial period based on three synchronized Greenland ice-core records: refining and extending the INTIMATE event stratigraphy

Due to their outstanding resolution and well-constrained chronologies, Greenland ice-core records provide a master record of past climatic changes throughout the Last Interglacial–Glacial cycle in the North Atlantic region. As part of the INTIMATE (INTegration of Ice-core, MArine and TErrestrial records) project, protocols have been proposed to ensure consistent and robust correlation between different records of past climate. A key element of these protocols has been the formal definition and ordinal numbering of the sequence of Greenland Stadials (GS) and Greenland Interstadials (GI) within the most recent glacial period. The GS and GI periods are the Greenland expressions of the characteristic Dansgaard–Oeschger events that represent cold and warm phases of the North Atlantic region, respectively. We present here a more detailed and extended GS/GI template for the whole of the Last Glacial period. It is based on a synchronization of the NGRIP, GRIP, and GISP2 ice-core records that allows the parallel analysis of all three records on a common time scale. The boundaries of the GS and GI periods are defined based on a combination of stable-oxygen isotope ratios of the ice (δ18O, reflecting mainly local temperature) and calcium ion concentrations (reflecting mainly atmospheric dust loading) measured in the ice. The data not only resolve the well-known sequence of Dansgaard–Oeschger events that were first defined and numbered in the ice-core records more than two decades ago, but also better resolve a number of short-lived climatic oscillations, some defined here for the first time. Using this revised scheme, we propose a consistent approach for discriminating and naming all the significant abrupt climatic events of the Last Glacial period that are represented in the Greenland ice records. The final product constitutes an extended and better resolved Greenland stratotype sequence, against which other proxy records can be compared and correlated. It also provides a more secure basis for investigating the dynamics and fundamental causes of these climatic perturbations.

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.

Electron and photon energy calibration with the ATLAS detector using LHC Run 1 data

This paper presents the electron and photon energy calibration achieved with the ATLAS detector using about 25 fb−1 of LHC proton–proton collision data taken at centre-of-mass energies of √s = 7 and 8 TeV. The reconstruction of electron and photon energies is optimised using multivariate algorithms. The response of the calorimeter layers is equalised in data and simulation, and the longitudinal profile of the electromagnetic showers is exploited to estimate the passive material in front of the calorimeter and reoptimise the detector simulation. After all corrections, the Z resonance is used to set the absolute energy scale. For electrons from Z decays, the achieved calibration is typically accurate to 0.05% in most of the detector acceptance, rising to 0.2% in regions with large amounts of passive material. The remaining inaccuracy is less than 0.2–1% for electrons with a transverse energy of 10 GeV, and is on average 0.3% for photons. The detector resolution is determined with a relative inaccuracy of less than 10% for electrons and photons up to 60 GeV transverse energy, rising to 40% for transverse energies above 500 GeV.

Measurement of the muon reconstruction performance of the ATLAS detector using 2011 and 2012 LHC proton-proton collision data

This paper presents the performance of the ATLAS muon reconstruction during the LHC run with pp collisions at √s = 7–8 TeV in 2011–2012, focusing mainly on data collected in 2012. Measurements of the reconstruction efficiency and of the momentum scale and resolution, based on large reference samples of J/ψ → μμ, Z → μμ and ϒ → μμ decays, are presented and compared to Monte Carlo simulations. Corrections to the simulation, to be used in physics analysis, are provided. Over most of the covered phase space (muon |η| < 2.7 and 5 ≲ pT ≲ 100 GeV) the efficiency is above 99% and is measured with per-mille precision. The momentum resolution ranges from 1.7% at central rapidity and for transverse momentum pT ≅ 10 GeV, to 4% at large rapidity and pT ≅ 100 GeV. The momentum scale is known with an uncertainty of 0.05% to 0.2% depending on rapidity. A method for the recovery of final state radiation from the muons is also presented.