Simulating the impact of four control strategies on the population dynamics of Neospora caninum infection in Swiss dairy cattle

A dynamic deterministic simulation model was developed to assess the impact of different putative control strategies on the seroprevalence of Neospora caninum in female Swiss dairy cattle. The model structure comprised compartments of "susceptible" and "infected" animals (SI-model) and the cattle population was divided into 12 age classes. A reference model (Model 1) was developed to simulate the current (status quo) situation (present seroprevalence in Switzerland 12%), taking into account available demographic and seroprevalence data of Switzerland. Model 1 was modified to represent four putative control strategies: testing and culling of seropositive animals (Model 2), discontinued breeding with offspring from seropositive cows (Model 3), chemotherapeutic treatment of calves from seropositive cows (Model 4), and vaccination of susceptible and infected animals (Model 5). Models 2-4 considered different sub-scenarios with regard to the frequency of diagnostic testing. Multivariable Monte Carlo sensitivity analysis was used to assess the impact of uncertainty in input parameters. A policy of annual testing and culling of all seropositive cattle in the population reduced the seroprevalence effectively and rapidly from 12% to <1% in the first year of simulation. The control strategies with discontinued breeding with offspring from all seropositive cows, chemotherapy of calves and vaccination of all cattle reduced the prevalence more slowly than culling but were still very effective (reduction of prevalence below 2% within 11, 23 and 3 years of simulation, respectively). However, sensitivity analyses revealed that the effectiveness of these strategies depended strongly on the quality of the input parameters used, such as the horizontal and vertical transmission factors, the sensitivity of the diagnostic test and the efficacy of medication and vaccination. Finally, all models confirmed that it was not possible to completely eradicate N. caninum as long as the horizontal transmission process was not interrupted.

Structural response and dynamics of fluid-structure-control interaction in wind turbine blades

Reducing the uncertainties related to blade dynamics by the improvement of the quality of numerical simulations of the fluid structure interaction process is a key for a breakthrough in wind-turbine technology. A fundamental step in that direction is the implementation of aeroelastic models capable of capturing the complex features of innovative prototype blades, so they can be tested at realistic full-scale conditions with a reasonable computational cost. We make use of a code based on a combination of two advanced numerical models implemented in a parallel HPC supercomputer platform: First, a model of the structural response of heterogeneous composite blades, based on a variation of the dimensional reduction technique proposed by Hodges and Yu. This technique has the capacity of reducing the geometrical complexity of the blade section into a stiffness matrix for an equivalent beam. The reduced 1-D strain energy is equivalent to the actual 3-D strain energy in an asymptotic sense, allowing accurate modeling of the blade structure as a 1-D finite-element problem. This substantially reduces the computational effort required to model the structural dynamics at each time step. Second, a novel aerodynamic model based on an advanced implementation of the BEM(Blade ElementMomentum) Theory; where all velocities and forces are re-projected through orthogonal matrices into the instantaneous deformed configuration to fully include the effects of large displacements and rotation of the airfoil sections into the computation of aerodynamic forces. This allows the aerodynamic model to take into account the effects of the complex flexo-torsional deformation that can be captured by the more sophisticated structural model mentioned above. In this thesis we have successfully developed a powerful computational tool for the aeroelastic analysis of wind-turbine blades. Due to the particular features mentioned above in terms of a full representation of the combined modes of deformation of the blade as a complex structural part and their effects on the aerodynamic loads, it constitutes a substantial advancement ahead the state-of-the-art aeroelastic models currently available, like the FAST-Aerodyn suite. In this thesis, we also include the results of several experiments on the NREL-5MW blade, which is widely accepted today as a benchmark blade, together with some modifications intended to explore the capacities of the new code in terms of capturing features on blade-dynamic behavior, which are normally overlooked by the existing aeroelastic models.

A Nonstationary Space-Time Gaussian Process Model for Partially Converged Simulations

In the context of expensive numerical experiments, a promising solution for alleviating the computational costs consists of using partially converged simulations instead of exact solutions. The gain in computational time is at the price of precision in the response. This work addresses the issue of fitting a Gaussian process model to partially converged simulation data for further use in prediction. The main challenge consists of the adequate approximation of the error due to partial convergence, which is correlated in both design variables and time directions. Here, we propose fitting a Gaussian process in the joint space of design parameters and computational time. The model is constructed by building a nonstationary covariance kernel that reflects accurately the actual structure of the error. Practical solutions are proposed for solving parameter estimation issues associated with the proposed model. The method is applied to a computational fluid dynamics test case and shows significant improvement in prediction compared to a classical kriging model.

Dynamics of vessel wall changes following the implantation of the Absorb everolimus-eluting bioresorbable vascular scaffold: a multi-imaging modality study at 6, 12, 24 and 36 months

Aims: To assess observations with multimodality imaging of the Absorb bioresorbable everolimus-eluting vascular scaffold performed in two consecutive cohorts of patients who were serially investigated either at 6 and 24 months or at 12 and 36 months. Methods and results: In the ABSORB multicentre single-arm trial, 45 patients (cohort B1) and 56 patients (cohort B2) underwent serial invasive imaging, specifically quantitative coronary angiography (QCA), intravascular ultrasound (IVUS), radiofrequency backscattering (IVUS-VH) and optical coherence tomography (OCT). Between one and three years, late luminal loss remained unchanged (6 months: 0.19 mm, 1 year: 0.27 mm, 2 years: 0.27 mm, 3 years: 0.29 mm) and the in-segment angiographic restenosis rate for the entire cohort B (n=101) at three years was 6%. On IVUS, mean lumen, scaffold, plaque and vessel area showed enlargement up to two years. Mean lumen and scaffold area remained stable between two and three years whereas significant reduction in plaque behind the struts occurred with a trend toward adaptive restrictive remodelling of EEM. Hyperechogenicity of the vessel wall, a surrogate of the bioresorption process, decreased from 23.1% to 10.4% with a reduction of radiofrequency backscattering for dense calcium and necrotic core. At three years, the count of strut cores detected on OCT increased significantly, probably reflecting the dismantling of the scaffold; 98% of struts were covered. In the entire cohort B (n=101), the three-year major adverse cardiac event rate was 10.0% without any scaffold thrombosis. Conclusions: The current investigation demonstrated the dynamics of vessel wall changes after implantation of a bioresorbable scaffold, resulting at three years in stable luminal dimensions, a low restenosis rate and a low clinical major adverse cardiac events rate.

The role of zinc dynamics in growth hormone secretion

Human growth hormone (GH) causes a variety of physiological and metabolic effects in humans and plays a pivotal role in postnatal growth. In somatotroph cells of the anterior pituitary, GH is stored in concentrated forms in secretory granules to be rapidly released upon GH-releasing hormone stimulation. During the process of secretory granule biogenesis, self-association of GH occurs in the compartments of the early secretory pathway (endoplasmic reticulum and Golgi complex). Since this process is greatly facilitated by the presence of zinc ions, it is of importance to understand the potential role of zinc transporters that participate in the fine-tuning of zinc homeostasis and dynamics, particularly in the early secretory pathway. Thus, the role of zinc transporters in supplying the secretory pathway with the sufficient amount of zinc required for the biogenesis of GH-containing secretory granules is essential for normal secretion. This report, illustrated by a clinical case report on transient neonatal zinc deficiency, focuses on the role of zinc in GH storage in the secretory granules and highlights the role of specific zinc transporters in the early secretory pathway.

Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

Correct estimation of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice core studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: outputs of a firn densification model, and measurements of δ15N of N2 in air trapped in ice core, assuming that δ15N is only affected by gravitational fractionation in the firn column. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi-coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available ice core air-δ15N measurements from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rates and temperature conditions. Our δ15N profiles reveal a heterogeneous response of the firn structure to glacial–interglacial climatic changes. While firn densification simulations correctly predict TALDICE δ15N variations, they systematically fail to capture the large millennial-scale δ15N variations measured at BI and the δ15N glacial levels measured at JRI and EDML – a mismatch previously reported for central East Antarctic ice cores. New constraints of the EDML gas–ice depth offset during the Laschamp event (~41 ka) and the last deglaciation do not favour the hypothesis of a large convective zone within the firn as the explanation of the glacial firn model–δ15N data mismatch for this site. While we could not conduct an in-depth study of the influence of impurities in snow for firnification from the existing datasets, our detailed comparison between the δ15N profiles and firn model simulations under different temperature and accumulation rate scenarios suggests that the role of accumulation rate may have been underestimated in the current description of firnification models.

Transient simulations of the carbon and nitrogen dynamics in northern peatlands: from the Last Glacial Maximum to the 21st century

The development of northern high-latitude peatlands played an important role in the carbon (C) balance of the land biosphere since the Last Glacial Maximum (LGM). At present, carbon storage in northern peatlands is substantial and estimated to be 500 ± 100 Pg C (1 Pg C = 1015 g C). Here, we develop and apply a peatland module embedded in a dynamic global vegetation and land surface process model (LPX-Bern 1.0). The peatland module features a dynamic nitrogen cycle, a dynamic C transfer between peatland acrotelm (upper oxic layer) and catotelm (deep anoxic layer), hydrology- and temperature-dependent respiration rates, and peatland specific plant functional types. Nitrogen limitation down-regulates average modern net primary productivity over peatlands by about half. Decadal acrotelm-to-catotelm C fluxes vary between −20 and +50 g C m−2 yr−1 over the Holocene. Key model parameters are calibrated with reconstructed peat accumulation rates from peat-core data. The model reproduces the major features of the peat core data and of the observation-based modern circumpolar soil carbon distribution. Results from a set of simulations for possible evolutions of northern peat development and areal extent show that soil C stocks in modern peatlands increased by 365–550 Pg C since the LGM, of which 175–272 Pg C accumulated between 11 and 5 kyr BP. Furthermore, our simulations suggest a persistent C sequestration rate of 35–50 Pg C per 1000 yr in present-day peatlands under current climate conditions, and that this C sink could either sustain or turn towards a source by 2100 AD depending on climate trajectories as projected for different representative greenhouse gas concentration pathways.

Methane emissions from floodplains in the Amazon Basin: challenges in developing a process-based model for global applications

Tropical wetlands are estimated to represent about 50% of the natural wetland methane (CH4) emissions and explain a large fraction of the observed CH4 variability on timescales ranging from glacial–interglacial cycles to the currently observed year-to-year variability. Despite their importance, however, tropical wetlands are poorly represented in global models aiming to predict global CH4 emissions. This publication documents a first step in the development of a process-based model of CH4 emissions from tropical floodplains for global applications. For this purpose, the LPX-Bern Dynamic Global Vegetation Model (LPX hereafter) was slightly modified to represent floodplain hydrology, vegetation and associated CH4 emissions. The extent of tropical floodplains was prescribed using output from the spatially explicit hydrology model PCR-GLOBWB. We introduced new plant functional types (PFTs) that explicitly represent floodplain vegetation. The PFT parameterizations were evaluated against available remote-sensing data sets (GLC2000 land cover and MODIS Net Primary Productivity). Simulated CH4 flux densities were evaluated against field observations and regional flux inventories. Simulated CH4 emissions at Amazon Basin scale were compared to model simulations performed in the WETCHIMP intercomparison project. We found that LPX reproduces the average magnitude of observed net CH4 flux densities for the Amazon Basin. However, the model does not reproduce the variability between sites or between years within a site. Unfortunately, site information is too limited to attest or disprove some model features. At the Amazon Basin scale, our results underline the large uncertainty in the magnitude of wetland CH4 emissions. Sensitivity analyses gave insights into the main drivers of floodplain CH4 emission and their associated uncertainties. In particular, uncertainties in floodplain extent (i.e., difference between GLC2000 and PCR-GLOBWB output) modulate the simulated emissions by a factor of about 2. Our best estimates, using PCR-GLOBWB in combination with GLC2000, lead to simulated Amazon-integrated emissions of 44.4 ± 4.8 Tg yr−1. Additionally, the LPX emissions are highly sensitive to vegetation distribution. Two simulations with the same mean PFT cover, but different spatial distributions of grasslands within the basin, modulated emissions by about 20%. Correcting the LPX-simulated NPP using MODIS reduces the Amazon emissions by 11.3%. Finally, due to an intrinsic limitation of LPX to account for seasonality in floodplain extent, the model failed to reproduce the full dynamics in CH4 emissions but we proposed solutions to this issue. The interannual variability (IAV) of the emissions increases by 90% if the IAV in floodplain extent is accounted for, but still remains lower than in most of the WETCHIMP models. While our model includes more mechanisms specific to tropical floodplains, we were unable to reduce the uncertainty in the magnitude of wetland CH4 emissions of the Amazon Basin. Our results helped identify and prioritize directions towards more accurate estimates of tropical CH4 emissions, and they stress the need for more research to constrain floodplain CH4 emissions and their temporal variability, even before including other fundamental mechanisms such as floating macrophytes or lateral water fluxes.

Intra-annual dynamics of non-structural carbohydrates in the cambium of mature conifer trees reflects radial growth demands

The presence of soluble carbohydrates in the cambial zone, either from sugars recently produced during photosynthesis or from starch remobilized from storage organs, is necessary for radial tree growth. However, considerable uncertainties on carbohydrate dynamics and the consequences on tree productivity exist. This study aims to better understand the variation in different carbon pools at intra-annual resolution by quantifying how cambial zone sugar and starch concentrations fluctuate over the season and in relation to cambial phenology. A comparison between two physiologically different species growing at the same site, i.e., the evergreen Picea abies Karst. and the deciduous Larix decidua Mill., and between L. decidua from two contrasting elevations, is presented to identify mechanisms of growth limitation. Results indicate that the annual cycle of sugar concentration within the cambial zone is coupled to the process of wood formation. The highest sugar concentration is observed when the number of cells in secondary wall formation and lignification stages is at a maximum, subsequent to most radial growth. Starch disappears in winter, while other freeze-resistant non-structural carbohydrates (NSCs) increase. Slight differences in NSC concentration between species are consistent with the differing climate sensitivity of the evergreen and deciduous species investigated. The general absence of differences between elevations suggests that the cambial activity of trees growing at the treeline was not limited by the availability of carbohydrates at the cambial zone but instead by environmental controls on the growing season duration.

Terahertz macrospin dynamics in insulating ferrimagnets

We investigate numerically the excitation of nonlinear magnetic interactions in a ferrite material by an energetic pump pulse of terahertz (THz) radiation. The calculations are performed by solving the coupled Maxwell and Landau-Lifshitz-Gilbert differential equations. In a time-resolved THz pump/THz probe scheme, it is demonstrated that Faraday rotation of a delayed THz probe pulse can be used to map these interactions. Our study is motivated by the ability of soft x-ray free electron lasers to perform time-resolved imaging of the magnetization process at the submicrometer and subpicosecond length and time scales.

Coupled human and natural system dynamics as key to the sustainability of Lake Victoria's ecosystem services

East Africa’s Lake Victoria provides resources and services to millions of people on the lake’s shores and abroad. In particular, the lake’s fisheries are an important source of protein, employment, and international economic connections for the whole region. Nonetheless, stock dynamics are poorly understood and currently unpredictable. Furthermore, fishery dynamics are intricately connected to other supporting services of the lake as well as to lakeshore societies and economies. Much research has been carried out piecemeal on different aspects of Lake Victoria’s system; e.g., societies, biodiversity, fisheries, and eutrophication. However, to disentangle drivers and dynamics of change in this complex system, we need to put these pieces together and analyze the system as a whole. We did so by first building a qualitative model of the lake’s social-ecological system. We then investigated the model system through a qualitative loop analysis, and finally examined effects of changes on the system state and structure. The model and its contextual analysis allowed us to investigate system-wide chain reactions resulting from disturbances. Importantly, we built a tool that can be used to analyze the cascading effects of management options and establish the requirements for their success. We found that high connectedness of the system at the exploitation level, through fisheries having multiple target stocks, can increase the stocks’ vulnerability to exploitation but reduce society’s vulnerability to variability in individual stocks. We describe how there are multiple pathways to any change in the system, which makes it difficult to identify the root cause of changes but also broadens the management toolkit. Also, we illustrate how nutrient enrichment is not a self-regulating process, and that explicit management is necessary to halt or reverse eutrophication. This model is simple and usable to assess system-wide effects of management policies, and can serve as a paving stone for future quantitative analyses of system dynamics at local scales.

The Emergence of Formal Control Specificity in Information Systems Outsourcing: A Process-View

Information systems (IS) outsourcing projects often fail to achieve initial goals. To avoid project failure, managers need to design formal controls that meet the specific contextual demands of the project. However, the dynamic and uncertain nature of IS outsourcing projects makes it difficult to design such specific formal controls at the outset of a project. It is hence crucial to translate high-level project goals into specific formal controls during the course of a project. This study seeks to understand the underlying patterns of such translation processes. Based on a comparative case study of four outsourced software development projects, we inductively develop a process model that consists of three unique patterns. The process model shows that the performance implications of emergent controls with higher specificity depend on differences in the translation process. Specific formal controls have positive implications for goal achievement if only the stakeholder context is adapted, while they are negative for goal achievement if in the translation process tasks are unintendedly adapted. In the latter case projects incrementally drift away from their initial direction. Our findings help to better understand control dynamics in IS outsourcing projects. We contribute to a process theoretic understanding of IS outsourcing governance and we derive implications for control theory and the IS project escalation literature.

Mechanisms of recognition and binding of α-TTP to the plasma membrane by multi-scale molecular dynamics simulations

We used multiple sets of simulations both at the atomistic and coarse-grained level of resolution to investigate interaction and binding of α-tochoperol transfer protein (α-TTP) to phosphatidylinositol phosphate lipids (PIPs). Our calculations indicate that enrichment of membranes with such lipids facilitate membrane anchoring. Atomistic models suggest that PIP can be incorporated into the binding cavity of α-TTP and therefore confirm that such protein can work as lipid exchanger between the endosome and the plasma membrane. Comparison of the atomistic models of the α-TTP-PIPs complex with membrane-bound α-TTP revealed different roles for the various basic residues composing the basic patch that is key for the protein/ligand interaction. Such residues are of critical importance as several point mutations at their position lead to severe forms of ataxia with vitamin E deficiency (AVED) phenotypes. Specifically, R221 is main residue responsible for the stabilization of the complex. R68 and R192 exchange strong interactions in the protein or in the membrane complex only, suggesting that the two residues alternate contact formation, thus facilitating lipid flipping from the membrane into the protein cavity during the lipid exchange process. Finally, R59 shows weaker interactions with PIPs anyway with a clear preference for specific phosphorylation positions, hinting a role in early membrane selectivity for the protein. Altogether, our simulations reveal significant aspects at the atomistic scale of interactions of α-TTP with the plasma membrane and with PIP, providing clarifications on the mechanism of intracellular vitamin E trafficking and helping establishing the role of key residue for the functionality of α-TTP.