Mechanisms of the negative shortwave cloud feedback in mid to high latitudes

Increases in cloud optical depth and liquid water path (LWP) are robust features of global warming model simulations in high latitudes, yielding a negative shortwave cloud feedback, but the mechanisms are still uncertain. We assess the importance of microphysical processes for the negative optical depth feedback by perturbing temperature in the microphysics schemes of two aquaplanet models, both of which have separate prognostic equations for liquid water and ice. We find that most of the LWP increase with warming is caused by a suppression of ice microphysical processes in mixed-phase clouds, resulting in reduced conversion efficiencies of liquid water to ice and precipitation. Perturbing the temperature-dependent phase partitioning of convective condensate also yields a small LWP increase. Together, the perturbations in large-scale microphysics and convective condensate partitioning explain more than two-thirds of the LWP response relative to a reference case with increased SSTs, and capture all of the vertical structure of the liquid water response. In support of these findings, we show the existence of a very robust positive relationship between monthly-mean LWP and temperature in CMIP5 models and observations in mixed-phase cloud regions only. In models, the historical LWP sensitivity to temperature is a good predictor of the forced global warming response poleward of about 45°, although models appear to overestimate the LWP response to warming compared to observations. We conclude that in climate models, the suppression of ice-phase microphysical processes that deplete cloud liquid water is a key driver of the LWP increase with warming and of the associated negative shortwave cloud feedback.

On univoque points for self-similar sets

Let K⊆R be the unique attractor of an iterated function system. We consider the case where K is an interval and study those elements of K with a unique coding. We prove under mild conditions that the set of points with a unique coding can be identified with a subshift of finite type. As a consequence, we can show that the set of points with a unique coding is a graph-directed self-similar set in the sense of Mauldin and Williams (1988). The theory of Mauldin and Williams then provides a method by which we can explicitly calculate the Hausdorff dimension of this set. Our algorithm can be applied generically, and our result generalises the work of Daróczy, Kátai, Kallós, Komornik and de Vries.

Technical variation in a sample of high level judo players

Technical actions performed by two groups of judokas who won medals at World Championships and Olympic Games during the period 1995-2001 were analyzed. In the Super Elite group (n = 17) were the best athletes in each weight category. The Elite group (n = 16) were medal winners who were not champions and did not win more than three medals. Super Elite judokas used a greater number of throwing techniques which resulted in scores, even when expressed relative to the total number of matches performed, and these techniques were applied in more directions than those of Elite judokas. Further, the number of different throwing techniques and the variability of directions in which techniques were applied were significantly correlated with number of wins and the number of points and ippon scored. Thus, a greater number of throwing techniques and use of directions for attack seem to be important in increasing unpredictability during judo matches.

Invariant theory and reversible-equivariant vector fields

In this paper we present results for the systematic study of reversible-equivariant vector fields - namely, in the simultaneous presence of symmetries and reversing symmetries - by employing algebraic techniques from invariant theory for compact Lie groups. The Hilbert-Poincare series and their associated Molien formulae are introduced,and we prove the character formulae for the computation of dimensions of spaces of homogeneous anti-invariant polynomial functions and reversible-equivariant polynomial mappings. A symbolic algorithm is obtained for the computation of generators for the module of reversible-equivariant polynomial mappings over the ring of invariant polynomials. We show that this computation can be obtained directly from a well-known situation, namely from the generators of the ring of invariants and the module of the equivariants. (C) 2008 Elsevier B.V, All rights reserved.

Iteration of closed geodesics in stationary Lorentzian manifolds

Following the lines of Bott in (Commun Pure Appl Math 9:171-206, 1956), we study the Morse index of the iterates of a closed geodesic in stationary Lorentzian manifolds, or, more generally, of a closed Lorentzian geodesic that admits a timelike periodic Jacobi field. Given one such closed geodesic gamma, we prove the existence of a locally constant integer valued map Lambda(gamma) on the unit circle with the property that the Morse index of the iterated gamma(N) is equal, up to a correction term epsilon(gamma) is an element of {0,1}, to the sum of the values of Lambda(gamma) at the N-th roots of unity. The discontinuities of Lambda(gamma) occur at a finite number of points of the unit circle, that are special eigenvalues of the linearized Poincare map of gamma. We discuss some applications of the theory.

Genetic heterogeneity of residual variance : estimation of variance components using double hierarchical generalized linear models

Background: The sensitivity to microenvironmental changes varies among animals and may be under genetic control. It is essential to take this element into account when aiming at breeding robust farm animals. Here, linear mixed models with genetic effects in the residual variance part of the model can be used. Such models have previously been fitted using EM and MCMC algorithms. Results: We propose the use of double hierarchical generalized linear models (DHGLM), where the squared residuals are assumed to be gamma distributed and the residual variance is fitted using a generalized linear model. The algorithm iterates between two sets of mixed model equations, one on the level of observations and one on the level of variances. The method was validated using simulations and also by re-analyzing a data set on pig litter size that was previously analyzed using a Bayesian approach. The pig litter size data contained 10,060 records from 4,149 sows. The DHGLM was implemented using the ASReml software and the algorithm converged within three minutes on a Linux server. The estimates were similar to those previously obtained using Bayesian methodology, especially the variance components in the residual variance part of the model. Conclusions: We have shown that variance components in the residual variance part of a linear mixed model can be estimated using a DHGLM approach. The method enables analyses of animal models with large numbers of observations. An important future development of the DHGLM methodology is to include the genetic correlation between the random effects in the mean and residual variance parts of the model as a parameter of the DHGLM.

Optimization Of Groundwater Remediation With New Efficient Parallel Algorithm For Global Optimization

Application of optimization algorithm to PDE modeling groundwater remediation can greatly reduce remediation cost. However, groundwater remediation analysis requires a computational expensive simulation, therefore, effective parallel optimization could potentially greatly reduce computational expense. The optimization algorithm used in this research is Parallel Stochastic radial basis function. This is designed for global optimization of computationally expensive functions with multiple local optima and it does not require derivatives. In each iteration of the algorithm, an RBF is updated based on all the evaluated points in order to approximate expensive function. Then the new RBF surface is used to generate the next set of points, which will be distributed to multiple processors for evaluation. The criteria of selection of next function evaluation points are estimated function value and distance from all the points known. Algorithms created for serial computing are not necessarily efficient in parallel so Parallel Stochastic RBF is different algorithm from its serial ancestor. The application for two Groundwater Superfund Remediation sites, Umatilla Chemical Depot, and Former Blaine Naval Ammunition Depot. In the study, the formulation adopted treats pumping rates as decision variables in order to remove plume of contaminated groundwater. Groundwater flow and contamination transport is simulated with MODFLOW-MT3DMS. For both problems, computation takes a large amount of CPU time, especially for Blaine problem, which requires nearly fifty minutes for a simulation for a single set of decision variables. Thus, efficient algorithm and powerful computing resource are essential in both cases. The results are discussed in terms of parallel computing metrics i.e. speedup and efficiency. We find that with use of up to 24 parallel processors, the results of the parallel Stochastic RBF algorithm are excellent with speed up efficiencies close to or exceeding 100%.

The role of no-arbitrage on forecasting: lessons from a parametric term structure model

Parametric term structure models have been successfully applied to innumerous problems in fixed income markets, including pricing, hedging, managing risk, as well as studying monetary policy implications. On their turn, dynamic term structure models, equipped with stronger economic structure, have been mainly adopted to price derivatives and explain empirical stylized facts. In this paper, we combine flavors of those two classes of models to test if no-arbitrage affects forecasting. We construct cross section (allowing arbitrages) and arbitrage-free versions of a parametric polynomial model to analyze how well they predict out-of-sample interest rates. Based on U.S. Treasury yield data, we find that no-arbitrage restrictions significantly improve forecasts. Arbitrage-free versions achieve overall smaller biases and Root Mean Square Errors for most maturities and forecasting horizons. Furthermore, a decomposition of forecasts into forward-rates and holding return premia indicates that the superior performance of no-arbitrage versions is due to a better identification of bond risk premium.

Graphical models and point set matching

Point pattern matching in Euclidean Spaces is one of the fundamental problems in Pattern Recognition, having applications ranging from Computer Vision to Computational Chemistry. Whenever two complex patterns are encoded by two sets of points identifying their key features, their comparison can be seen as a point pattern matching problem. This work proposes a single approach to both exact and inexact point set matching in Euclidean Spaces of arbitrary dimension. In the case of exact matching, it is assured to find an optimal solution. For inexact matching (when noise is involved), experimental results confirm the validity of the approach. We start by regarding point pattern matching as a weighted graph matching problem. We then formulate the weighted graph matching problem as one of Bayesian inference in a probabilistic graphical model. By exploiting the existence of fundamental constraints in patterns embedded in Euclidean Spaces, we prove that for exact point set matching a simple graphical model is equivalent to the full model. It is possible to show that exact probabilistic inference in this simple model has polynomial time complexity with respect to the number of elements in the patterns to be matched. This gives rise to a technique that for exact matching provably finds a global optimum in polynomial time for any dimensionality of the underlying Euclidean Space. Computational experiments comparing this technique with well-known probabilistic relaxation labeling show significant performance improvement for inexact matching. The proposed approach is significantly more robust under augmentation of the sizes of the involved patterns. In the absence of noise, the results are always perfect.

Method of numerical correction of errors occasioned by delay of records during the monitoring of environmental variables of interest for animal production

As condições de ambiente térmico e aéreo, no interior de instalações para animais, alteram-se durante o dia, devido à influência do ambiente externo. Para que análises estatísticas e geoestatísticas sejam representativas, uma grande quantidade de pontos distribuídos espacialmente na área da instalação deve ser monitorada. Este trabalho propõe que a variação no tempo das variáveis ambientais de interesse para a produção animal, monitoradas no interior de instalações para animais, pode ser modelada com precisão a partir de registros discretos no tempo. O objetivo deste trabalho foi desenvolver um método numérico para corrigir as variações temporais dessas variáveis ambientais, transformando os dados para que tais observações independam do tempo gasto durante a aferição. O método proposto aproximou os valores registrados com retardos de tempo aos esperados no exato momento de interesse, caso os dados fossem medidos simultaneamente neste momento em todos os pontos distribuídos espacialmente. O modelo de correção numérica para variáveis ambientais foi validado para o parâmetro ambiental temperatura do ar, sendo que os valores corrigidos pelo método não diferiram pelo teste Tukey, a 5% de probabilidade dos valores reais registrados por meio de dataloggers.

Tornar-se, ser e viver do professorado: entre regularidades e variações identitárias

The main aim of this research has been to analyze the identity patterns of the teacher s staff of fundamental education public schools in the Metropolitan Area of Natal-RN. It sets out from the hypothesis that being a teacher within this context grows out of the regularities of a specific habitus, which, according to Bourdieu, develops into mental schemes of thought and action within a specific social group. This habitus forms the basis on which is built the social representation of being a teacher prevailing in the group, as well as the symbolic differences that typify its identity variations. Three data sources have been fundamental in building up this thesis: (a) formative essays of students graduating from a Higher Teacher s Formation Course, as well as observing some of the public defense of these essays during field work; (b) a questionnaire aimed at classifying economically, socially, and culturally a sample of public teachers of the Natal-RN county; and (c) submitting a sub-sample of this group to the process of Multiple Classification Procedures (MCP). The analysis of data was done according to the multidimensional, non-parametric statistical procedures of both the Category Content Analysis and Enunciation Analysis methods. The results of the analysis took into account an ample set of variables, its associations and implications, the cultural and social profile of the population under scrutiny, their life styles, as well as the strategies they developed in the process of becoming a teacher, and the social representation of being a teacher . We came to the conclusion that the social identity of the teachers corps, or as we prefer to say it being a teacher , is a result of a set of regularities produced by the habitus that gives social shape and meaning to the existence of the group proper. We note the existence of identity variations caused by the variables (a) educational level; and (b) mode of action in fundamental education (if these are the first or last grades where the subjects operate). However, these variations will not break the power of the regularities that give shape, meaning, and social visibility to the group. The social representation of being a teacher points to the tensions, ambiguities, and trends inherent to common sense, as well as to a strong tendency to reassign a new meaning to being a teacher. Our thesis, therefore, is that the identity configuration of the teachers corps under scrutiny is characterized by an integrative synthesis, by-product of a habitus that is superimposed, and at the same time co-exists with different identity variations

Caracterização e avaliação do macro sistema urbano de distribuição de água de Natal RN

The aim of this study is to characterize and evaluate the Macro System of Regional Water Distribution Natal North (RNN) and Southern Regional Natal (RNS), covering 35% and 65% respectively of the Natal-RN City. The terms of the quality and quantity of water (surface and groundwater) were also evaluated in order to adjust the parameters that contribute to proper distribution and control in water reserves. The methodology of the work took place from collecting volumetric data of production capacity and distribution of the two treatment plants for Regional as well as the flow rates of wells. Yet the quantitative capacity of reservation, distribution and consumption of the main reservoirs, population numbers and consumption of members neighborhoods were collected. Data were tabulated and used in computational simulator EPANET to diagnose possible through the water balance, the offers and demands on the water supply system in the neighborhoods of the capital, linking them to specific distribution points. We also evaluated the wells in the levels of nitrate in water consumed. As a result it was found that some neighborhoods in the South Regional Natal, was ranked as critical supply situation: City of Hope, Lagoa Nova and Nova Descoberta, where demand exceeds supply. While in most Northern Regional Natal present deficiency in the supply system as: Lagoa Azul, the Parque dos Coqueiros, igapó, Amarante and Salinas. The rates of nitrate in the city were significant, but manageable with corrective and preventive measures. The averages were 12 mg /l-N in Candelária, 10 mg/l-N in Lagoa Nova, 9 mg/l-N in Satelite, 20 mg/l-N in Gramore and 15 mg/l-N in N. Sra. Apresentação. Therefore proper distribution of water abstracted and implementation of quality control ensures the supply required by the system, associated with preservation of Water Resources of the Metropolitan Region of Natal

Prodrugs for the treatment of neglected diseases

Recently, World Health Organization ( WHO) and Medicins San Frontieres (MSF) proposed a classification of diseases as global, neglected and extremely neglected. Global diseases, such as cancer, cardiovascular and mental (CNS) diseases represent the targets of the majority of the R&D efforts of pharmaceutical companies. Neglected diseases affect millions of people in the world yet existing drug therapy is limited and often inappropriate. Furthermore, extremely neglected diseases affect people living under miserable conditions who barely have access to the bare necessities for survival. Most of these diseases are excluded from the goals of the R&D programs in the pharmaceutical industry and therefore fall outside the pharmaceutical market. About 14 million people, mainly in developing countries, die each year from infectious diseases. From 1975 to 1999, 1393 new drugs were approved yet only 1% were for the treatment of neglected diseases [ 3]. These numbers have not changed until now, so in those countries there is an urgent need for the design and synthesis of new drugs and in this area the prodrug approach is a very interesting field. It provides, among other effects, activity improvements and toxicity decreases for current and new drugs, improving market availability. It is worth noting that it is essential in drug design to save time and money, and prodrug approaches can be considered of high interest in this respect. The present review covers 20 years of research on the design of prodrugs for the treatment of neglected and extremely neglected diseases such as Chagas' disease ( American trypanosomiasis), sleeping sickness ( African trypanosomiasis), malaria, sickle cell disease, tuberculosis, leishmaniasis and schistosomiasis.

A fully adaptive multiresolution scheme for shock computations

The scheme is based on Ami Harten's ideas (Harten, 1994), the main tools coming from wavelet theory, in the framework of multiresolution analysis for cell averages. But instead of evolving cell averages on the finest uniform level, we propose to evolve just the cell averages on the grid determined by the significant wavelet coefficients. Typically, there are few cells in each time step, big cells on smooth regions, and smaller ones close to irregularities of the solution. For the numerical flux, we use a simple uniform central finite difference scheme, adapted to the size of each cell. If any of the required neighboring cell averages is not present, it is interpolated from coarser scales. But we switch to ENO scheme in the finest part of the grids. To show the feasibility and efficiency of the method, it is applied to a system arising in polymer-flooding of an oil reservoir. In terms of CPU time and memory requirements, it outperforms Harten's multiresolution algorithm.The proposed method applies to systems of conservation laws in 1Dpartial derivative(t)u(x, t) + partial derivative(x)f(u(x, t)) = 0, u(x, t) is an element of R-m. (1)In the spirit of finite volume methods, we shall consider the explicit schemeupsilon(mu)(n+1) = upsilon(mu)(n) - Deltat/hmu ((f) over bar (mu) - (f) over bar (mu)-) = [Dupsilon(n)](mu), (2)where mu is a point of an irregular grid Gamma, mu(-) is the left neighbor of A in Gamma, upsilon(mu)(n) approximate to 1/mu-mu(-) integral(mu-)(mu) u(x, t(n))dx are approximated cell averages of the solution, (f) over bar (mu) = (f) over bar (mu)(upsilon(n)) are the numerical fluxes, and D is the numerical evolution operator of the scheme.According to the definition of (f) over bar (mu), several schemes of this type have been proposed and successfully applied (LeVeque, 1990). Godunov, Lax-Wendroff, and ENO are some of the popular names. Godunov scheme resolves well the shocks, but accuracy (of first order) is poor in smooth regions. Lax-Wendroff is of second order, but produces dangerous oscillations close to shocks. ENO schemes are good alternatives, with high order and without serious oscillations. But the price is high computational cost.Ami Harten proposed in (Harten, 1994) a simple strategy to save expensive ENO flux calculations. The basic tools come from multiresolution analysis for cell averages on uniform grids, and the principle is that wavelet coefficients can be used for the characterization of local smoothness.. Typically, only few wavelet coefficients are significant. At the finest level, they indicate discontinuity points, where ENO numerical fluxes are computed exactly. Elsewhere, cheaper fluxes can be safely used, or just interpolated from coarser scales. Different applications of this principle have been explored by several authors, see for example (G-Muller and Muller, 1998).Our scheme also uses Ami Harten's ideas. But instead of evolving the cell averages on the finest uniform level, we propose to evolve the cell averages on sparse grids associated with the significant wavelet coefficients. This means that the total number of cells is small, with big cells in smooth regions and smaller ones close to irregularities. This task requires improved new tools, which are described next.

Structural aspects of the fermion-boson mapping in two-dimensional gauge and anomalous gauge theories with massive fermions

Using a synthesis of the functional integral and operator approaches we discuss the fermion-buson mapping and the role played by the Bose field algebra in the Hilbert space of two-dimensional gauge and anomalous gauge field theories with massive fermions. In QED, with quartic self-interaction among massive fermions, the use of an auxiliary vector field introduces a redundant Bose field algebra that should not be considered as an element of the intrinsic algebraic structure defining the model. In anomalous chiral QED, with massive fermions the effect of the chiral anomaly leads to the appearance in the mass operator of a spurious Bose field combination. This phase factor carries no fermion selection rule and the expected absence of Theta-vacuum in the anomalous model is displayed from the operator solution. Even in the anomalous model with massive Fermi fields, the introduction of the Wess-Zumino field replicates the theory, changing neither its algebraic content nor its physical content. (C) 2002 Elsevier B.V. (USA).