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.

A unification of models for meta-analysis of diagnostic accuracy studies

Studies of diagnostic accuracy require more sophisticated methods for their meta-analysis than studies of therapeutic interventions. A number of different, and apparently divergent, methods for meta-analysis of diagnostic studies have been proposed, including two alternative approaches that are statistically rigorous and allow for between-study variability: the hierarchical summary receiver operating characteristic (ROC) model (Rutter and Gatsonis, 2001) and bivariate random-effects meta-analysis (van Houwelingen and others, 1993), (van Houwelingen and others, 2002), (Reitsma and others, 2005). We show that these two models are very closely related, and define the circumstances in which they are identical. We discuss the different forms of summary model output suggested by the two approaches, including summary ROC curves, summary points, confidence regions, and prediction regions.

Impacts of changing climate and land use on vegetation dynamics in a Mediterranean ecosystem: insights from paleoecology and dynamic modeling

Forests near the Mediterranean coast have been shaped by millennia of human disturbance. Consequently, ecological studies relying on modern observations or historical records may have difficulty assessing natural vegetation dynamics under current and future climate. We combined a sedimentary pollen record from Lago di Massacciucoli, Tuscany, Italy with simulations from the LandClim dynamic vegetation model to determine what vegetation preceded intense human disturbance, how past changes in vegetation relate to fire and browsing, and the potential of an extinct vegetation type under present climate. We simulated vegetation dynamics near Lago di Massaciucoli for the last 7,000 years using a local chironomid-inferred temperature reconstruction with combinations of three fire regimes (small infrequent, large infrequent, small frequent) and three browsing intensities (no browsing, light browsing, and moderate browsing), and compared model output to pollen data. Simulations with low disturbance support pollen-inferred evidence for a mixed forest dominated by Quercus ilex (a Mediterranean species) and Abies alba (a montane species). Whereas pollen data record the collapse of A. alba after 6000 cal yr bp, simulated populations expanded with declining summer temperatures during the late Holocene. Simulations with increased fire and browsing are consistent with evidence for expansion by deciduous species after A. alba collapsed. According to our combined paleo-environmental and modeling evidence, mixed Q. ilex and A. alba forests remain possible with current climate and limited disturbance, and provide a viable management objective for ecosystems near the Mediterranean coast and in regions that are expected to experience a mediterranean-type climate in the future.

Grazing response patterns indicate isolation of semi-natural European grasslands

Identifying drivers of species diversity is a major challenge in understanding and predicting the dynamics of species-rich semi-natural grasslands. In particular in temperate grasslands changes in land use and its consequences, i.e. increasing fragmentation, the on-going loss of habitat and the declining importance of regional processes such as seed dispersal by livestock, are considered key drivers of the diversity loss witnessed within the last decades. It is a largely unresolved question to what degree current temperate grassland communities already reflect a decline of regional processes such as longer distance seed dispersal. Answering this question is challenging since it requires both a mechanistic approach to community dynamics and a sufficient data basis that allows identifying general patterns. Here, we present results of a local individual- and trait-based community model that was initialized with plant functional types (PFTs) derived from an extensive empirical data set of species-rich grasslands within the `Biodiversity Exploratories' in Germany. Driving model processes included above- and belowground competition, dynamic resource allocation to shoots and roots, clonal growth, grazing, and local seed dispersal. To test for the impact of regional processes we also simulated seed input from a regional species pool. Model output, with and without regional seed input, was compared with empirical community response patterns along a grazing gradient. Simulated response patterns of changes in PFT richness, Shannon diversity, and biomass production matched observed grazing response patterns surprisingly well if only local processes were considered. Already low levels of additional regional seed input led to stronger deviations from empirical community pattern. While these findings cannot rule out that regional processes other than those considered in the modeling study potentially play a role in shaping the local grassland communities, our comparison indicates that European grasslands are largely isolated, i.e. local mechanisms explain observed community patterns to a large extent.

Effects of climate, fire, and humans on forest dynamics: forest simulations compared to the palaeological record

In order to find out which factors influenced the forest dynamics in northern Italy during the Holocene, a palaeoecological approach involving pollen analysis was combined with ecosystem modelling. The dynamic and distribution based forest model DisCForm was run with different input scenarios for climate, species immigration, fire, and human impact and the similarity of the simulations with the original pollen record was assessed. From the comparisons of the model output and the pollen core, it appears that immigration was most important in the first part of the Holocene, and that fire and human activity had a major influence in the second half. Species not well represented in the simulation outputs are species with a higher abundance in the past than today (Corylus), with their habitat in riparian forests (Alnus) or with a strong response to human impact (Castanea).

A model for radiofrequency electromagnetic field predictions at outdoor and indoor locations in the context of epidemiological research

We present a geospatial model to predict the radiofrequency electromagnetic field from fixed site transmitters for use in epidemiological exposure assessment. The proposed model extends an existing model toward the prediction of indoor exposure, that is, at the homes of potential study participants. The model is based on accurate operation parameters of all stationary transmitters of mobile communication base stations, and radio broadcast and television transmitters for an extended urban and suburban region in the Basel area (Switzerland). The model was evaluated by calculating Spearman rank correlations and weighted Cohen's kappa (kappa) statistics between the model predictions and measurements obtained at street level, in the homes of volunteers, and in front of the windows of these homes. The correlation coefficients of the numerical predictions with street level measurements were 0.64, with indoor measurements 0.66, and with window measurements 0.67. The kappa coefficients were 0.48 (95%-confidence interval: 0.35-0.61) for street level measurements, 0.44 (95%-CI: 0.32-0.57) for indoor measurements, and 0.53 (95%-CI: 0.42-0.65) for window measurements. Although the modeling of shielding effects by walls and roofs requires considerable simplifications of a complex environment, we found a comparable accuracy of the model for indoor and outdoor points.

Simultaneous Multiscale Polyaffine Registration by Incorporating Deformation Statistics

Locally affine (polyaffine) image registration methods capture intersubject non-linear deformations with a low number of parameters, while providing an intuitive interpretation for clinicians. Considering the mandible bone, anatomical shape differences can be found at different scales, e.g. left or right side, teeth, etc. Classically, sequential coarse to fine registration are used to handle multiscale deformations, instead we propose a simultaneous optimization of all scales. To avoid local minima we incorporate a prior on the polyaffine transformations. This kind of groupwise registration approach is natural in a polyaffine context, if we assume one configuration of regions that describes an entire group of images, with varying transformations for each region. In this paper, we reformulate polyaffine deformations in a generative statistical model, which enables us to incorporate deformation statistics as a prior in a Bayesian setting. We find optimal transformations by optimizing the maximum a posteriori probability. We assume that the polyaffine transformations follow a normal distribution with mean and concentration matrix. Parameters of the prior are estimated from an initial coarse to fine registration. Knowing the region structure, we develop a blockwise pseudoinverse to obtain the concentration matrix. To our knowledge, we are the first to introduce simultaneous multiscale optimization through groupwise polyaffine registration. We show results on 42 mandible CT images.

Closed-loop control of mean arterial blood pressure during surgery with alfentanil: clinical evaluation of a novel model-based predictive controller

BACKGROUND: In contrast to hypnosis, there is no surrogate parameter for analgesia in anesthetized patients. Opioids are titrated to suppress blood pressure response to noxious stimulation. The authors evaluated a novel model predictive controller for closed-loop administration of alfentanil using mean arterial blood pressure and predicted plasma alfentanil concentration (Cp Alf) as input parameters. METHODS: The authors studied 13 healthy patients scheduled to undergo minor lumbar and cervical spine surgery. After induction with propofol, alfentanil, and mivacurium and tracheal intubation, isoflurane was titrated to maintain the Bispectral Index at 55 (+/- 5), and the alfentanil administration was switched from manual to closed-loop control. The controller adjusted the alfentanil infusion rate to maintain the mean arterial blood pressure near the set-point (70 mmHg) while minimizing the Cp Alf toward the set-point plasma alfentanil concentration (Cp Alfref) (100 ng/ml). RESULTS: Two patients were excluded because of loss of arterial pressure signal and protocol violation. The alfentanil infusion was closed-loop controlled for a mean (SD) of 98.9 (1.5)% of presurgery time and 95.5 (4.3)% of surgery time. The mean (SD) end-tidal isoflurane concentrations were 0.78 (0.1) and 0.86 (0.1) vol%, the Cp Alf values were 122 (35) and 181 (58) ng/ml, and the Bispectral Index values were 51 (9) and 52 (4) before surgery and during surgery, respectively. The mean (SD) absolute deviations of mean arterial blood pressure were 7.6 (2.6) and 10.0 (4.2) mmHg (P = 0.262), and the median performance error, median absolute performance error, and wobble were 4.2 (6.2) and 8.8 (9.4)% (P = 0.002), 7.9 (3.8) and 11.8 (6.3)% (P = 0.129), and 14.5 (8.4) and 5.7 (1.2)% (P = 0.002) before surgery and during surgery, respectively. A post hoc simulation showed that the Cp Alfref decreased the predicted Cp Alf compared with mean arterial blood pressure alone. CONCLUSION: The authors' controller has a similar set-point precision as previous hypnotic controllers and provides adequate alfentanil dosing during surgery. It may help to standardize opioid dosing in research and may be a further step toward a multiple input-multiple output controller.

ECG-triggered muscular counterpulsation for treatment of low cardiac output

BACKGROUND: Skeletal muscular counterpulsation (MCP) has been used as a new noninvasive technique for treatment of low cardiac output. The MCP method is based on ECG-triggered skeletal muscle stimulation. The purpose of the present study was to evaluate acute hemodynamic changes induced by MCP in the experimental animal. METHODS: Eight anaesthetized pigs (43+/-4 kg) were studied at rest and after IV â-blockade (10 mg propranolol) before and after MCP. Muscular counterpulsation was performed on both thighs using trains (75 ms duration) of multiple biphasic electrical impulses with a width of 1 ms and a frequency of 200 Hz at low (10 V) and high (30 V) amplitude. ECG-triggering was used to synchronize stimulation to a given time point. LV pressure-volume relations were determined using the conductance catheter. After baseline measurements, MCP was carried out for 10 minutes at low and high stimulation amplitude. The optimal time point for MCP was determined from LV pressure-volume loops using different stimulation time points during systole and diastole. Best results were observed during end-systole and, therefore, this time point was used for stimulation. RESULTS: Under control conditions, MCP was associated with a significant decrease in pulmonary vascular resistance (-18%), a decrease in systemic vascular resistance (-11%) and stroke work index (-4%), whereas cardiac index (+2%) and ejection fraction (+6%) increased slightly. Pressure-volume loops showed a leftward shift with a decrease in end-systolic volume. After â-blockade, cardiac function decreased (HR, MAP, EF, dP/dt max), but it improved with skeletal muscle stimulation (HR +10% and CI +17%, EF +5%). There was a significant decrease in pulmonary (-19%) and systemic vascular resistance (-29%). CONCLUSIONS: In the animal model, ECG-triggered skeletal muscular counterpulsation is associated with a significant improvement in cardiac function at baseline and after IV â-blockade. Thus, MCP represents a new, non-invasive technique which improves cardiac function by diastolic compression of the peripheral arteries and veins, with a decrease in systemic vascular resistance and increase in cardiac output.

Human temporal bones versus mechanical model to evaluate three middle ear transducers

A life-size mechanical middle ear model and human temporal bones were used to evaluate three different middle ear transducers for implantable hearing aids: the driving rod transducer (DRT), the floating mass transducer (FMT) or vibrant sound bridge, and the contactless transducer (CLT). Results of the experiments with the mechanical model were within the range of the results for human temporal bones. However, results with the mechanical model showed better reproducibility. The handling of the mechanical model was considerably simpler and less time-consuming. Systematic variations of mounting parameters showed that the angle of the rod has virtually no effect on the output of the DRT, the mass loading on the cable of the FMT has a larger impact on the output than does the tightness of crimping, and the output level of the CLT can be increased by 10 dB by optimizing the mounting parameters.

Effects of an aging vascular model on healthy and diseased hearts

The vitamin D(3) and nicotine (VDN) model is a model of isolated systolic hypertension (ISH) due to arterial calcification raising arterial stiffness and vascular impedance similar to an aged and stiffened arterial tree. We therefore analyzed the impact of this aging model on normal and diseased hearts with myocardial infarction (MI). Wistar rats were treated with VDN (n = 9), subjected to MI by coronary ligation (n = 10), or subjected to a combination of both MI and VDN treatment (VDN/MI, n = 14). A sham-treated group served as control (Ctrl, n = 10). Transthoracic echocardiography was performed every 2 wk, whereas invasive indexes were obtained at week 8 before death. Calcium, collagen, and protein contents were measured in the heart and the aorta. Systolic blood pressure, pulse pressure, thoracic aortic calcium, and end-systolic elastance as an index of myocardial contractility were highest in the aging model group compared with MI and Ctrl groups (P(VDN) < 0.05, 2-way ANOVA). Left ventricular wall stress and brain natriuretic peptide (P(VDNxMI) = not significant) were highest, while ejection fraction, stroke volume, and cardiac output were lowest in the combined group versus all other groups (P(VDNxMI) < 0.05). The combination of ISH due to this aging model and MI demonstrates significant alterations in cardiac function. This model mimics several clinical phenomena of cardiovascular aging and may thus serve to further study novel therapies.

A multiple source model for 6 MV photon beam dose calculations using Monte Carlo

A multiple source model (MSM) for the 6 MV beam of a Varian Clinac 2300 C/D was developed by simulating radiation transport through the accelerator head for a set of square fields using the GEANT Monte Carlo (MC) code. The corresponding phase space (PS) data enabled the characterization of 12 sources representing the main components of the beam defining system. By parametrizing the source characteristics and by evaluating the dependence of the parameters on field size, it was possible to extend the validity of the model to arbitrary rectangular fields which include the central 3 x 3 cm2 field without additional precalculated PS data. Finally, a sampling procedure was developed in order to reproduce the PS data. To validate the MSM, the fluence, energy fluence and mean energy distributions determined from the original and the reproduced PS data were compared and showed very good agreement. In addition, the MC calculated primary energy spectrum was verified by an energy spectrum derived from transmission measurements. Comparisons of MC calculated depth dose curves and profiles, using original and PS data reproduced by the MSM, agree within 1% and 1 mm. Deviations from measured dose distributions are within 1.5% and 1 mm. However, the real beam leads to some larger deviations outside the geometrical beam area for large fields. Calculated output factors in 10 cm water depth agree within 1.5% with experimentally determined data. In conclusion, the MSM produces accurate PS data for MC photon dose calculations for the rectangular fields specified.