The Origins of Lactase Persistence in Europe

Lactase persistence (LP) is common among people of European ancestry, but with the exception of some African, Middle Eastern and southern Asian groups, is rare or absent elsewhere in the world. Lactase gene haplotype conservation around a polymorphism strongly associated with LP in Europeans (-13,910 C/T) indicates that the derived allele is recent in origin and has been subject to strong positive selection. Furthermore, ancient DNA work has shown that the -13,910*T (derived) allele was very rare or absent in early Neolithic central Europeans. It is unlikely that LP would provide a selective advantage without a supply of fresh milk, and this has lead to a gene-culture coevolutionary model where lactase persistence is only favoured in cultures practicing dairying, and dairying is more favoured in lactase persistent populations. We have developed a flexible demic computer simulation model to explore the spread of lactase persistence, dairying, other subsistence practices and unlinked genetic markers in Europe and western Asia's geographic space. Using data on -13,910*T allele frequency and farming arrival dates across Europe, and approximate Bayesian computation to estimate parameters of interest, we infer that the -13,910*T allele first underwent selection among dairying farmers around 7,500 years ago in a region between the central Balkans and central Europe, possibly in association with the dissemination of the Neolithic Linearbandkeramik culture over Central Europe. Furthermore, our results suggest that natural selection favouring a lactase persistence allele was not higher in northern latitudes through an increased requirement for dietary vitamin D. Our results provide a coherent and spatially explicit picture of the coevolution of lactase persistence and dairying in Europe.

Linking the coevolutionary and population dynamics of host-parasitoid interactions

The interplay between coevolutionary and population or community dynamics is currently the focus of much empirical and theoretical consideration. Here, we develop a simulation model to study the coevolutionary and population dynamics of a hypothetical host-parasitoid interaction. In the model, host resistance and parasitoid virulence are allowed to coevolve. We investigate how trade-offs associated with these traits modify the system's coevolutionary and population dynamics. The most important influence on these dynamics comes from the incorporation of density-dependent costs of resistance ability. We find three main outcomes. First, if the costs of resistance are high, then one or both of the players go extinct. Second, when the costs of resistance are intermediate to low, cycling population and coevolutionary dynamics are found, with slower evolutionary changes observed when the costs of virulence are also low. Third, when the costs associated with resistance and virulence are both high, the hosts trade-off resistance against fecundity and invest little in resistance. However, the parasitoids continue to invest in virulence, leading to stable host and parasitoid population sizes. These results support the hypothesis that costs associated with resistance and virulence will maintain the heritable variation in these traits found in natural populations and that the nature of these trade-offs will greatly influence the population dynamics of the interacting species.

Constraint release in entangled binary blends of linear polymers: a molecular dynamics study

We present extensive molecular dynamics simulations of the dynamics of diluted long probe chains entangled with a matrix of shorter chains. The chain lengths of both components are above the entanglement strand length, and the ratio of their lengths is varied over a wide range to cover the crossover from the chain reptation regime to tube Rouse motion regime of the long probe chains. Reducing the matrix chain length results in a faster decay of the dynamic structure factor of the probe chains, in good agreement with recent neutron spin echo experiments. The diffusion of the long chains, measured by the mean square displacements of the monomers and the centers of mass of the chains, demonstrates a systematic speed-up relative to the pure reptation behavior expected for monodisperse melts of sufficiently long polymers. On the other hand, the diffusion of the matrix chains is only weakly perturbed by the diluted long probe chains. The simulation results are qualitatively consistent with the theoretical predictions based on constraint release Rouse model, but a detailed comparison reveals the existence of a broad distribution of the disentanglement rates, which is partly confirmed by an analysis of the packing and diffusion of the matrix chains in the tube region of the probe chains. A coarse-grained simulation model based on the tube Rouse motion model with incorporation of the probability distribution of the tube segment jump rates is developed and shows results qualitatively consistent with the fine scale molecular dynamics simulations. However, we observe a breakdown in the tube Rouse model when the short chain length is decreased to around N-S = 80, which is roughly 3.5 times the entanglement spacing N-e(P) = 23. The location of this transition may be sensitive to the chain bending potential used in our simulations.

Is housing good for the economy?

Not surprisingly housing researchers and practitioners frequently call for more resources to be devoted to housing. But governments in recent years have devoted fewer resources to housing rather than more. One of the reasons is that housing expenditures have to be seen in terms of the overall resource constraints on the economy and in many instances this requires a macro‐economic perspective. This paper reviews the macro‐economic arguments for and against housing expenditures, particularly through the use of a quantitative policy simulation model.

Uncertainty and feasibility studies: an Italian case study

Uncertainty affects all aspects of the property market but one area where the impact of uncertainty is particularly significant is within feasibility analyses. Any development is impacted by differences between market conditions at the conception of the project and the market realities at the time of completion. The feasibility study needs to address the possible outcomes based on an understanding of the current market. This requires the appraiser to forecast the most likely outcome relating to the sale price of the completed development, the construction costs and the timing of both. It also requires the appraiser to understand the impact of finance on the project. All these issues are time sensitive and analysis needs to be undertaken to show the impact of time to the viability of the project. The future is uncertain and a full feasibility analysis should be able to model the upside and downside risk pertaining to a range of possible outcomes. Feasibility studies are extensively used in Italy to determine land value but they tend to be single point analysis based upon a single set of “likely” inputs. In this paper we look at the practical impact of uncertainty in variables using a simulation model (Crystal Ball ©) with an actual case study of an urban redevelopment plan for an Italian Municipality. This allows the appraiser to address the issues of uncertainty involved and thus provide the decision maker with a better understanding of the risk of development. This technique is then refined using a “two-dimensional technique” to distinguish between “uncertainty” and “variability” and thus create a more robust model.

How predation and landscape fragmentation affect vole population dynamics

Background: Microtine species in Fennoscandia display a distinct north-south gradient from regular cycles to stable populations. The gradient has often been attributed to changes in the interactions between microtines and their predators. Although the spatial structure of the environment is known to influence predator-prey dynamics of a wide range of species, it has scarcely been considered in relation to the Fennoscandian gradient. Furthermore, the length of microtine breeding season also displays a north-south gradient. However, little consideration has been given to its role in shaping or generating population cycles. Because these factors covary along the gradient it is difficult to distinguish their effects experimentally in the field. The distinction is here attempted using realistic agent-based modelling. Methodology/Principal Findings: By using a spatially explicit computer simulation model based on behavioural and ecological data from the field vole (Microtus agrestis), we generated a number of repeated time series of vole densities whose mean population size and amplitude were measured. Subsequently, these time series were subjected to statistical autoregressive modelling, to investigate the effects on vole population dynamics of making predators more specialised, of altering the breeding season, and increasing the level of habitat fragmentation. We found that fragmentation as well as the presence of specialist predators are necessary for the occurrence of population cycles. Habitat fragmentation and predator assembly jointly determined cycle length and amplitude. Length of vole breeding season had little impact on the oscillations. Significance: There is good agreement between our results and the experimental work from Fennoscandia, but our results allow distinction of causation that is hard to unravel in field experiments. We hope our results will help understand the reasons for cycle gradients observed in other areas. Our results clearly demonstrate the importance of landscape fragmentation for population cycling and we recommend that the degree of fragmentation be more fully considered in future analyses of vole dynamics.

Interactions between the physical soil environment and a horizontal ground coupled heat pump, for a domestic site in the UK

There is currently an increased interest of Government and Industry in the UK, as well as at the European Community level and International Agencies (i.e. Department of Energy, American International Energy Agency), to improve the performance and uptake of Ground Coupled Heat Pumps (GCHP), in order to meet the 2020 renewable energy target. A sound knowledge base is required to help inform the Government Agencies and advisory bodies; detailed site studies providing reliable data for model verification have an important role to play in this. In this study we summarise the effect of heat extraction by a horizontal ground heat exchanger (installed at 1 m depth) on the soil physical environment (between 0 and 1 m depth) for a site in the south of the UK. Our results show that the slinky influences the surrounding soil by significantly decreasing soil temperatures. Furthermore, soil moisture contents were lower for the GCHP soil profile, most likely due to temperature-gradient related soil moisture migration effects and a decreased hydraulic conductivity, the latter as a result of increased viscosity (caused by the lower temperatures for the GCHP soil profile). The effects also caused considerable differences in soil thermal properties. This is the first detailed mechanistic study conducted in the UK with the aim to understand the interactions between the soil, horizontal heat exchangers and the aboveground environment. An increased understanding of these interactions will help to achieve an optimum and sustainable use of the soil heat resources in the future. The results of this study will help to calibrate and verify a simulation model that will provide UK-wide recommendations to improve future GCHP uptake and performance, while safeguarding the soil physical resources.

Nonlinear set-membership estimation: a support vector machine approach

In this paper a support vector machine (SVM) approach for characterizing the feasible parameter set (FPS) in non-linear set-membership estimation problems is presented. It iteratively solves a regression problem from which an approximation of the boundary of the FPS can be determined. To guarantee convergence to the boundary the procedure includes a no-derivative line search and for an appropriate coverage of points on the FPS boundary it is suggested to start with a sequential box pavement procedure. The SVM approach is illustrated on a simple sine and exponential model with two parameters and an agro-forestry simulation model.

Changing climate and the Irish landscape

The impacts of current and future changes in climate have been investigated for Irish vegetation. Warming has been observed over the last two decades, with impacts that are also strongly influenced by natural oscillations of the surrounding ocean, seen as fluctuations in the North Atlantic Oscillation and the Atlantic Multidecadal Oscillation. Satellite observations show that vegetation greenness increases in warmer years, a feature mirrored by increases in net ecosystem production observed for a grassland and a plantation forest. An ensemble of general circulation model simulations of future climates indicate temperature rises over the twenty-first century ranging from 1°C to 7°C, depending on future scenarios of greenhouse gas emissions. Net primary production is simulated to increase under all scenarios, due to the positive impacts of rising temperature, a modest rise of precipitation and rising carbon dioxide concentrations. In an optimistic scenario of reducing future emissions, CO2 concentration is simulated to flatten from about 2070, although temperatures continue to increase. Under this scenario Ireland could become a source of carbon, whereas under all other emission scenarios Ireland is a sink for carbon that may increase by up to three-fold over the twenty-first century. A likely and unavoidable impact of changing climate is the arrival of alien plant species, which may disrupt ecosystems and exert negative impacts on native biodiversity. Alien species arrive continually, with about 250 dated arrivals in the twentieth century. A simulation model indicates that this rate of alien arrival may increase by anything between two and ten times, dependent on the future climatic scenario, by 2050. Which alien species may become severely disruptive is, however, not known.

Modelling predicts that heat stress, not drought, will increase vulnerability of wheat in Europe

New crop cultivars will be required for a changing climate characterised by increased summer drought and heat stress in Europe. However, the uncertainty in climate predictions poses a challenge to crop scientists and breeders who have limited time and resources and must select the most appropriate traits for improvement. Modelling is a powerful tool to quantify future threats to crops and hence identify targets for improvement. We have used a wheat simulation model combined with local-scale climate scenarios to predict impacts of heat stress and drought on winter wheat in Europe. Despite the lower summer precipitation projected for 2050s across Europe, relative yield losses from drought is predicted to be smaller in the future, because wheat will mature earlier avoiding severe drought. By contrast, the risk of heat stress around flowering will increase, potentially resulting in substantial yield losses for heat sensitive cultivars commonly grown in northern Europe.

Modelling of control options for an outbreak of Mycoplasma gallisepticum in egg production: a decision support tool

Mycoplasma gallisepticum (MG) is a bacterium that causes respiratory disease in chickens, leading to reduced egg production. A dynamic simulation model was developed that can be used to assess the costs and benefits of control using antimicrobials or vaccination in caged or free range systems. The intended users are veterinarians and egg producers. A user interface is provided for input of flock specific parameters. The economic consequence of an MG outbreak is expressed as a reduction in expected egg output. The model predicts that either vaccination or microbial treatment can approximately halve potential losses from MG in some circumstances. Sensitivity analysis is used to test assumptions about infection rate and timing of an outbreak. Feedback from veterinarians points to the value of the model as a discussion tool with producers.

Electron acceleration and parallel electric fields due to kinetic Alfvén waves in plasma with similar thermal and Alfvén speeds

We investigate electron acceleration due to shear Alfven waves in a collissionless plasma for plasma parameters typical of 4–5RE radial distance from the Earth along auroral field lines. Recent observational work has motivated this study, which explores the plasma regime where the thermal velocity of the electrons is similar to the Alfven speed of the plasma, encouraging Landau resonance for electrons in the wave fields. We use a self-consistent kinetic simulation model to follow the evolution of the electrons as they interact with a short-duration wave pulse, which allows us to determine the parallel electric field of the shear Alfven wave due to both electron inertia and electron pressure effects. The simulation demonstrates that electrons can be accelerated to keV energies in a modest amplitude sub-second period wave. We compare the parallel electric field obtained from the simulation with those provided by fluid approximations.

Modelling capability and affordance as properties of human/machine resource systems

Understanding how and why the capability of one set of business resources, its structural arrangements and mechanisms compared to another works can provide competitive advantage in terms of new business processes and product and service development. However, most business models of capability are descriptive and lack formal modelling language to qualitatively and quantifiably compare capabilities, Gibson’s theory of affordance, the potential for action, provides a formal basis for a more robust and quantitative model, but most formal affordance models are complex and abstract and lack support for real-world applications. We aim to understand the ‘how’ and ‘why’ of business capability, by developing a quantitative and qualitative model that underpins earlier work on Capability-Affordance Modelling – CAM. This paper integrates an affordance based capability model and the formalism of Coloured Petri Nets to develop a simulation model. Using the model, we show how capability depends on the space time path of interacting resources, the mechanism of transition and specific critical affordance factors relating to the values of the variables for resources, people and physical objects. We show how the model can identify the capabilities of resources to enable the capability to inject a drug and anaesthetise a patient.

Folding of graphene slit like pore walls—a simple method of improving CO2 separation from mixtures with CH4 or N2

We report for the first time a detailed procedure for creating a simulation model of energetically stable, folded graphene-like pores and simulation results of CO2/CH4 and CO2/N2 separation using these structures. We show that folding of graphene structures is a very promising method to improve the separation of CO2 from mixtures with CH4 and N2. The separation properties of the analyzed materials are compared with carbon nanotubes having similar diameters or S/V ratio. The presented results have potential importance in the field of CO2 capture and sequestration.

Simulation of the diurnal cycle in a climate model and its evaluation using data from Meteosat 7

The representation of the diurnal cycle in the Hadley Centre climate model is evaluated using simulations of the infrared radiances observed by Meteosat 7. In both the window and water vapour channels, the standard version of the model with 19 levels produces a good simulation of the geographical distributions of the mean radiances and of the amplitude of the diurnal cycle. Increasing the vertical resolution to 30 levels leads to further improvements in the mean fields. The timing of the maximum and minimum radiances reveals significant model errors, however, which are sensitive to the frequency with which the radiation scheme is called. In most regions, these errors are consistent with well documented errors in the timing of convective precipitation, which peaks before noon in the model, in contrast to the observed peak in the late afternoon or evening. When the radiation scheme is called every model time step (half an hour), as opposed to every three hours in the standard version, the timing of the minimum radiance is improved for convective regions over central Africa, due to the creation of upper-level layer-cloud by detrainment from the convection scheme, which persists well after the convection itself has dissipated. However, this produces a decoupling between the timing of the diurnal cycles of precipitation and window channel radiance. The possibility is raised that a similar decoupling may occur in reality and the implications of this for the retrieval of the diurnal cycle of precipitation from infrared radiances are discussed.