Earthworm distribution and abundance predicted by a process-based model

Earthworms are significant ecosystem engineers and are an important component of the diet of many vertebrates and invertebrates, so the ability to predict their distribution and abundance would have wide application in ecology, conservation and land management. Earthworm viability is known to be affected by the availability and quality of food resources, soil water conditions and temperature, but has not yet been modelled mechanistically to link effects on individuals to field population responses. Here we present a novel model capable of predicting the effects of land management and environmental conditions on the distribution and abundance of Aporrectodea caliginosa, the dominant earthworm species in agroecosystems. Our process-based approach uses individual based modelling (IBM), in which each individual has its own energy budget. Individual earthworm energy budgets follow established principles of physiological ecology and are parameterised for A. caliginosa from experimental measurements under optimal conditions. Under suboptimal conditions (e.g. food limitation, low soil temperatures and water contents) reproduction is prioritised over growth. Good model agreement to independent laboratory data on individual cocoon production and growth of body mass, under variable feeding and temperature conditions support our representation of A. caliginosa physiology through energy budgets. Our mechanistic model is able to accurately predict A. caliginosa distribution and abundance in spatially heterogeneous soil profiles representative of field study conditions. Essential here is the explicit modelling of earthworm behaviour in the soil profile. Local earthworm movement responds to a trade-off between food availability and soil water conditions, and this determines the spatiotemporal distribution of the population in the soil profile. Importantly, multiple environmental variables can be manipulated simultaneously in the model to explore earthworm population exposure and effects to combinations of stressors. Potential applications include prediction of the population-level effects of pesticides and changes in soil management e.g. conservation tillage and climate change.

The representation of the West-African Monsoon vertical cloud structure in the Met Office Unified Model: an evaluation with CloudSat

Weather and climate model simulations of the West African Monsoon (WAM) have generally poor representation of the rainfall distribution and monsoon circulation because key processes, such as clouds and convection, are poorly characterized. The vertical distribution of cloud and precipitation during the WAM are evaluated in Met Office Unified Model simulations against CloudSat observations. Simulations were run at 40-km and 12-km horizontal grid length using a convection parameterization scheme and at 12-km, 4-km, and 1.5-km grid length with the convection scheme effectively switched off, to study the impact of model resolution and convection parameterization scheme on the organisation of tropical convection. Radar reflectivity is forward-modelled from the model cloud fields using the CloudSat simulator to present a like-with-like comparison with the CloudSat radar observations. The representation of cloud and precipitation at 12-km horizontal grid length improves dramatically when the convection parameterization is switched off, primarily because of a reduction in daytime (moist) convection. Further improvement is obtained when reducing model grid length to 4 km or 1.5 km, especially in the representation of thin anvil and mid-level cloud, but three issues remain in all model configurations. Firstly, all simulations underestimate the fraction of anvils with cloud top height above 12 km, which can be attributed to too low ice water contents in the model compared to satellite retrievals. Secondly, the model consistently detrains mid-level cloud too close to the freezing level, compared to higher altitudes in CloudSat observations. Finally, there is too much low-level cloud cover in all simulations and this bias was not improved when adjusting the rainfall parameters in the microphysics scheme. To improve model simulations of the WAM, more detailed and in-situ observations of the dynamics and microphysics targeting these non-precipitating cloud types are required.

COCO-2: a model to assess the economic impact of an accident

COCO-2 is a model for assessing the potential economic costs likely to arise off-site following an accident at a nuclear reactor. COCO-2 builds on work presented in the model COCO-1 developed in 1991 by considering economic effects in more detail, and by including more sources of loss. Of particular note are: the consideration of the directly affected local economy, indirect losses that stem from the directly affected businesses, losses due to changes in tourism consumption, integration with the large body of work on recovery after an accident and a more systematic approach to health costs. The work, where possible, is based on official data sources for reasons of traceability, maintenance and ease of future development. This report describes the methodology and discusses the results of an example calculation. Guidance on how the base economic data can be updated in the future is also provided.

Sea-breeze dynamics and convection initiation: the influence of convective parameterization in weather and climate model biases

There are some long-established biases in atmospheric models that originate from the representation of tropical convection. Previously, it has been difficult to separate cause and effect because errors are often the result of a number of interacting biases. Recently, researchers have gained the ability to run multiyear global climate model simulations with grid spacings small enough to switch the convective parameterization off, which permits the convection to develop explicitly. There are clear improvements to the initiation of convective storms and the diurnal cycle of rainfall in the convection-permitting simulations, which enables a new process-study approach to model bias identification. In this study, multiyear global atmosphere-only climate simulations with and without convective parameterization are undertaken with the Met Office Unified Model and are analyzed over the Maritime Continent region, where convergence from sea-breeze circulations is key for convection initiation. The analysis shows that, although the simulation with parameterized convection is able to reproduce the key rain-forming sea-breeze circulation, the parameterization is not able to respond realistically to the circulation. A feedback of errors also occurs: the convective parameterization causes rain to fall in the early morning, which cools and wets the boundary layer, reducing the land–sea temperature contrast and weakening the sea breeze. This is, however, an effect of the convective bias, rather than a cause of it. Improvements to how and when convection schemes trigger convection will improve both the timing and location of tropical rainfall and representation of sea-breeze circulations.

Fluctuating abundance of humpback whales (Megaptera novaeangliae) in a calving ground off coastal Brazil

The humpback whale (Megaptera novaeangliae) population that uses Abrolhos Bank, off the east coast of Brazil as a breeding ground is increasing. To describe temporal changes in the relative abundance of humpback whales around Abrolhos, seven years (1998-2004) of whale count data were collected during July through to November. During one-hour-scans, observers determined group size within 9.3 km (5 n.m.) of a land-based observing station. A total Of 930 scans, comprising 7996 sightings of adults and 2044 calves were analysed using generalized linear models that included variables for time of day, day of the season, years and two-way interactions as possible predictors. The pattern observed was the gradual build-up and decline in whale counts within seasons. Patterns and peaks of adult and calf counts varied among years. Although fluctuation was observed, there was generally an increasing trend in adult counts among years. Calf counts increased only in 2004. These fluctuations may have been caused by some environmental conditions in humpback whales` summering grounds and also by changes in spatial-temporal concentrations in Abrolhos Bank. The general pattern observed within the study area mirrored what was observed in the whole Abrolhos Bank. Knowledge of the consistency with which humpback whales use this important nursing area should prove beneficial for designing future monitoring programmes especially related to whale watching activities around Abrolhos Archipelago.

A linear optimization approach to the combined production planning model

Two fundamental processes usually arise in the production planning of many industries. The first one consists of deciding how many final products of each type have to be produced in each period of a planning horizon, the well-known lot sizing problem. The other process consists of cutting raw materials in stock in order to produce smaller parts used in the assembly of final products, the well-studied cutting stock problem. In this paper the decision variables of these two problems are dependent of each other in order to obtain a global optimum solution. Setups that are typically present in lot sizing problems are relaxed together with integer frequencies of cutting patterns in the cutting problem. Therefore, a large scale linear optimizations problem arises, which is exactly solved by a column generated technique. It is worth noting that this new combined problem still takes the trade-off between storage costs (for final products and the parts) and trim losses (in the cutting process). We present some sets of computational tests, analyzed over three different scenarios. These results show that, by combining the problems and using an exact method, it is possible to obtain significant gains when compared to the usual industrial practice, which solve them in sequence. (C) 2010 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.

Evolutionary model trees for handling continuous classes in machine learning

Model trees are a particular case of decision trees employed to solve regression problems. They have the advantage of presenting an interpretable output, helping the end-user to get more confidence in the prediction and providing the basis for the end-user to have new insight about the data, confirming or rejecting hypotheses previously formed. Moreover, model trees present an acceptable level of predictive performance in comparison to most techniques used for solving regression problems. Since generating the optimal model tree is an NP-Complete problem, traditional model tree induction algorithms make use of a greedy top-down divide-and-conquer strategy, which may not converge to the global optimal solution. In this paper, we propose a novel algorithm based on the use of the evolutionary algorithms paradigm as an alternate heuristic to generate model trees in order to improve the convergence to globally near-optimal solutions. We call our new approach evolutionary model tree induction (E-Motion). We test its predictive performance using public UCI data sets, and we compare the results to traditional greedy regression/model trees induction algorithms, as well as to other evolutionary approaches. Results show that our method presents a good trade-off between predictive performance and model comprehensibility, which may be crucial in many machine learning applications. (C) 2010 Elsevier Inc. All rights reserved.

Melting Regime of the Anionic Phospholipid DMPG: New Lamellar Phase and Porous Bilayer Model

Aqueous dispersions of the anionic phospholipid dimyristoyl phosphatidylglycerol (DMPG) at pH above the apparent pK of DMPG and concentrations in the interval 70-300 mM have been investigated by small (SAXS) and wide-angle X-ray scattering, differential scanning calorimetry, and polarized optical microscopy. The order. disorder transition of the hydrocarbon chains occurs along an interval of about 10 degrees C (between T(m)(on) similar to 20 degrees C and T(m)(off) similar to 30 degrees C). Such melting regime was previously characterized at lower concentrations, up to 70 mM DMPG, when sample transparency was correlated with the presence of pores across the bilayer. At higher concentrations considered here, the melting regime persists but is not transparent. Defined SAXS peaks appear and a new lamellar phase L(p) with pores is proposed to exist above 70 mM DMPG, starting at similar to 23 degrees C (similar to 3 degrees C above T(m)(on)) and losing correlation after T(m)(off). A new model for describing the X-ray scattering of bilayers with pores, presented here, is able to explain the broad band attributed to in-plane correlation between pores. The majority of cell membranes have a net negative charge, and the opening of pores across the membrane tuned by ionic strength, temperature, and lipid composition is likely to have biological relevance.

The double Caldeira-Leggett model: Derivation and solutions of the master equations, reservoir-induced interactions and decoherence

In this paper we analyze the double Caldeira-Leggett model: the path integral approach to two interacting dissipative harmonic oscillators. Assuming a general form of the interaction between the oscillators, we consider two different situations: (i) when each oscillator is coupled to its own reservoir, and (ii) when both oscillators are coupled to a common reservoir. After deriving and solving the master equation for each case, we analyze the decoherence process of particular entanglements in the positional space of both oscillators. To analyze the decoherence mechanism we have derived a general decay function, for the off-diagonal peaks of the density matrix, which applies both to common and separate reservoirs. We have also identified the expected interaction between the two dissipative oscillators induced by their common reservoir. Such a reservoir-induced interaction, which gives rise to interesting collective damping effects, such as the emergence of relaxation- and decoherence-free subspaces, is shown to be blurred by the high-temperature regime considered in this study. However, we find that different interactions between the dissipative oscillators, described by rotating or counter-rotating terms, result in different decay rates for the interference terms of the density matrix. (C) 2010 Elsevier B.V. All rights reserved.

Culture-area relation in Axelrod`s model for culture dissemination

Axelrod`s model for culture dissemination offers a nontrivial answer to the question of why there is cultural diversity given that people`s beliefs have a tendency to become more similar to each other`s as they interact repeatedly. The answer depends on the two control parameters of the model, namely, the number F of cultural features that characterize each agent, and the number q of traits that each feature can take on, as well as on the size A of the territory or, equivalently, on the number of interacting agents. Here, we investigate the dependence of the number C of distinct coexisting cultures on the area A in Axelrod`s model, the culture-area relationship, through extensive Monte Carlo simulations. We find a non-monotonous culture-area relation, for which the number of cultures decreases when the area grows beyond a certain size, provided that q is smaller than a threshold value q (c) = q (c) (F) and F a parts per thousand yen 3. In the limit of infinite area, this threshold value signals the onset of a discontinuous transition between a globalized regime marked by a uniform culture (C = 1), and a completely polarized regime where all C = q (F) possible cultures coexist. Otherwise, the culture-area relation exhibits the typical behavior of the species-area relation, i.e., a monotonically increasing curve the slope of which is steep at first and steadily levels off at some maximum diversity value.

Mice and larvae tracking using a particle filter with an auto-adjustable observation model

This paper proposes a novel way to combine different observation models in a particle filter framework. This, so called, auto-adjustable observation model, enhance the particle filter accuracy when the tracked objects overlap without infringing a great runtime penalty to the whole tracking system. The approach has been tested under two important real world situations related to animal behavior: mice and larvae tracking. The proposal was compared to some state-of-art approaches and the results show, under the datasets tested, that a good trade-off between accuracy and runtime can be achieved using an auto-adjustable observation model. (C) 2009 Elsevier B.V. All rights reserved.

The trade-off between incentives and endogenous risk

Standard models of moral hazard predict a negative relationship between risk and incentives, but the empirical work has not confirmed this prediction. In this paper, we propose a model with adverse selection followed by moral hazard, where effort and the degree of risk aversion are private information of an agent who can control the mean and the variance of profits. For a given contract, more risk-averse agents suppIy more effort in risk reduction. If the marginal utility of incentives decreases with risk aversion, more risk-averse agents prefer lower-incentive contractsj thus, in the optimal contract, incentives are positively correlated with endogenous risk. In contrast, if risk aversion is high enough, the possibility of reduction in risk makes the marginal utility of incentives increasing in risk aversion and, in this case, risk and incentives are negatively related.

The risk-incentive trade-off in competitive search

Usando a abordagem de competitive search, modelo um mercado de trabalho com trabalhadores heterogêneos no qual há um problema de risco moral na relação entre firmas e trabalhadores. Nesse contexto, consigo prever como contratos reagem a mudanças nos parâmetros do mercado (em particular, o risco de produção), assim como a variação da probabilidade dos trabalhadores serem contratados. Minha contribuição principal é ver que, no nível individual, existe uma relação negativa entre risco e incentivos, mas efeitos de equilíbrio geral implicam que essa relação pode ser positiva no nível agregado. Esse resultado ajuda a esclarecer resultados empíricos contraditórios sobre a relação entre risco e incentivos.

Existe trade-off entre eficiência e qualidade nas organizações hospitalares?

A mudança do perfil demográfico e epidemiológico das populações, com progressivo envelhecimento populacional e aumento de portadores de doenças crônicas não transmissíveis, somado a necessidade da ampliação da oferta de serviços de saúde e crescentes custos em saúde, impõe enormes desafios aos sistemas e serviços de saúde. A eficiência organizacional dos serviços de saúde tem papel importante tanto na racionalização dos custos quanto na melhoria da qualidade e segurança assistencial. Tendo papel central nos sistemas de saúde como centros difusores de conhecimento, capacitação profissional, incorporação de tecnologias, prestação de serviços de maior complexidade aos pacientes e, consequentemente, elevados custos destes serviços, aos hospitais é fundamental a busca por essa eficiência. Este estudo buscou analisar se existe trade-off entre eficiência e qualidade em organizações hospitalares e identificar quais determinantes poderiam estar associados com maiores ou menores escores de eficiência. Utilizou-se dois modelos de análise de envelopamento de dados (data envelopment analysis, DEA), sem e com variáveis de qualidade, com retornos variáveis de escala e orientados para resultado. Foram estudados 47 hospitais gerais públicos do estado de São Paulo. No modelo sem variáveis de qualidade 14 deles foram considerados eficientes, enquanto que 33 no modelo com estas variáveis. O coeficiente de correlação de Spearman entre os dois modelos foi de 0,470 (correlação moderada). Não há evidências de que haja trade-off entre eficiência e qualidade nestas organizações hospitalares. Hospitais eficientes no modelo sem variáveis de qualidade, também o foram com variáveis de qualidade, assim como houve hospitais ineficientes no modelo sem variáveis de qualidade que foram eficientes com estas variáveis. Não foram encontradas associações estatisticamente significantes (p<0,05) entre eficiência e as características dos hospitais estudados, como acreditação, modelos de gestão, porte hospitalar e atividades de ensino, apesar de alguns achados de maior ou menor escore de eficiência para alguns determinantes. Desta maneira, concluiu-se que a utilização de variáveis de qualidade é um fator fundamental na determinação da eficiência de organizações de saúde, e não podem estar dissociadas. Gestões eficientes também estão relacionadas à obtenção de melhores resultados assistenciais sem a necessidade que se tenha de optar em alcançar melhores resultados econômico-financeiros ou melhores resultados assistenciais.

Is the meander growth in the Brazil Current system off Southeast Brazil due to baroclinic instability?

Temporally-growing frontal meandering and occasional eddy-shedding is observed in the Brazil Current (BC) as it flows adjacent to the Brazilian Coast. No study of the dynamics of this phenomenon has been conducted to date in the region between 22 degrees S and 25 degrees S. Within this latitude range, the flow over the intermediate continental slope is marked by a current inversion at a depth that is associated with the Intermediate Western Boundary Current (IWBC). A time series analysis of 10-current-meter mooring data was used to describe a mean vertical profile for the BC-IWBC jet and a typical meander vertical structure. The latter was obtained by an empirical orthogonal function (EOF) analysis that showed a single mode explaining 82% of the total variance. This mode structure decayed sharply with depth, revealing that the meandering is much more vigorous within the BC domain than it is in the IWBC region. As the spectral analysis of the mode amplitude time series revealed no significant periods, we searched for dominant wavelengths. This search was done via a spatial EOF analysis on 51 thermal front patterns derived from digitized AVHRR images. Four modes were statistically significant at the 95% confidence level. Modes 3 and 4, which together explained 18% of the total variance, are associated with 266 and 338-km vorticity waves, respectively. With this new information derived from the data, the [Johns, W.E., 1988. One-dimensional baroclinically unstable waves on the Gulf Stream potential vorticity gradient near Cape Hatteras. Dyn. Atmos. Oceans 11, 323-350] one-dimensional quasi-geostrophic model was applied to the interpolated mean BC-IWBC jet. The results indicated that the BC system is indeed baroclinically unstable and that the wavelengths depicted in the thermal front analysis are associated with the most unstable waves produced by the model. Growth rates were about 0.06 (0.05) days(-1) for the 266-km (338-km) wave. Moreover, phase speeds for these waves were low compared to the surface BC velocity and may account for remarks in the literature about growing standing or stationary meanders off southeast Brazil. The theoretical vertical structure modes associated with these waves resembled very closely to the one obtained for the current-meter mooring EOF analysis. We interpret this agreement as a confirmation that baroclinic instability is an important mechanism in meander growth in the BC system. (C) 2008 Elsevier B.V. All rights reserved.