Gaussian anamorphosis in the analysis step of the EnKF: a joint state-variable/observation approach

The analysis step of the (ensemble) Kalman filter is optimal when (1) the distribution of the background is Gaussian, (2) state variables and observations are related via a linear operator, and (3) the observational error is of additive nature and has Gaussian distribution. When these conditions are largely violated, a pre-processing step known as Gaussian anamorphosis (GA) can be applied. The objective of this procedure is to obtain state variables and observations that better fulfil the Gaussianity conditions in some sense. In this work we analyse GA from a joint perspective, paying attention to the effects of transformations in the joint state variable/observation space. First, we study transformations for state variables and observations that are independent from each other. Then, we introduce a targeted joint transformation with the objective to obtain joint Gaussianity in the transformed space. We focus primarily in the univariate case, and briefly comment on the multivariate one. A key point of this paper is that, when (1)-(3) are violated, using the analysis step of the EnKF will not recover the exact posterior density in spite of any transformations one may perform. These transformations, however, provide approximations of different quality to the Bayesian solution of the problem. Using an example in which the Bayesian posterior can be analytically computed, we assess the quality of the analysis distributions generated after applying the EnKF analysis step in conjunction with different GA options. The value of the targeted joint transformation is particularly clear for the case when the prior is Gaussian, the marginal density for the observations is close to Gaussian, and the likelihood is a Gaussian mixture.

An evaluation of decadal probability forecasts from state-of-the-art climate models

While state-of-the-art models of Earth's climate system have improved tremendously over the last 20 years, nontrivial structural flaws still hinder their ability to forecast the decadal dynamics of the Earth system realistically. Contrasting the skill of these models not only with each other but also with empirical models can reveal the space and time scales on which simulation models exploit their physical basis effectively and quantify their ability to add information to operational forecasts. The skill of decadal probabilistic hindcasts for annual global-mean and regional-mean temperatures from the EU Ensemble-Based Predictions of Climate Changes and Their Impacts (ENSEMBLES) project is contrasted with several empirical models. Both the ENSEMBLES models and a “dynamic climatology” empirical model show probabilistic skill above that of a static climatology for global-mean temperature. The dynamic climatology model, however, often outperforms the ENSEMBLES models. The fact that empirical models display skill similar to that of today's state-of-the-art simulation models suggests that empirical forecasts can improve decadal forecasts for climate services, just as in weather, medium-range, and seasonal forecasting. It is suggested that the direct comparison of simulation models with empirical models becomes a regular component of large model forecast evaluations. Doing so would clarify the extent to which state-of-the-art simulation models provide information beyond that available from simpler empirical models and clarify current limitations in using simulation forecasting for decision support. Ultimately, the skill of simulation models based on physical principles is expected to surpass that of empirical models in a changing climate; their direct comparison provides information on progress toward that goal, which is not available in model–model intercomparisons.

Estimating Lorenz's reference state in an ocean with a nonlinear equation of state for seawater

The study of the mechanical energy budget of the oceans using Lorenz available potential energy (APE) theory is based on knowledge of the adiabatically re-arranged Lorenz reference state of minimum potential energy. The compressible and nonlinear character of the equation of state for seawater has been thought to cause the reference state to be ill-defined, casting doubt on the usefulness of APE theory for investigating ocean energetics under realistic conditions. Using a method based on the volume frequency distribution of parcels as a function of temperature and salinity in the context of the seawater Boussinesq approximation, which we illustrate using climatological data, we show that compressibility effects are in fact minor. The reference state can be regarded as a well defined one-dimensional function of depth, which forms a surface in temperature, salinity and density space between the surface and the bottom of the ocean. For a very small proportion of water masses, this surface can be multivalued and water parcels can have up to two statically stable levels in the reference density profile, of which the shallowest is energetically more accessible. Classifying parcels from the surface to the bottom gives a different reference density profile than classifying in the opposite direction. However, this difference is negligible. We show that the reference state obtained by standard sorting methods is equivalent, though computationally more expensive, to the volume frequency distribution approach. The approach we present can be applied systematically and in a computationally efficient manner to investigate the APE budget of the ocean circulation using models or climatological data.

Histological Outcomes on the Development of New Space-making Devices for Maxillary Sinus Floor Augmentation

Background: Previous studies have pointed out that the mere elevation of the maxillary sinus membrane promotes bone formation without the use of augmentation materials. Purpose: This experimental study aimed at evaluating if the two-stage procedure for sinus floor augmentation could benefit from the use of a space-making device in order to increase the bone volume to enable later implant installation with good primary stability. Materials and Methods: Six male tufted capuchin primates (Cebus apella) were subjected to extraction of the three premolars and the first molar on both sides of the maxilla to create an edentulous area. The sinuses were opened using the lateral bone-wall window technique, and the membrane was elevated. One resorbable space-making device was inserted in each maxillary sinus, and the bone window was returned in place. The animals were euthanatized after 6 months, and biopsy blocks containing the whole maxillary sinus and surrounding soft tissues were prepared for ground sections. Results: The histological examination of the specimens showed bone formation in contact with both the schneiderian membrane and the device in most cases even when the device was displaced. The process of bone formation indicates that this technique is potentially useful for two-stage sinus floor augmentation. The lack of stabilization of the device within the sinus demands further improvement of space-makers for predictable bone augmentation. Conclusions: It is concluded that (1) the device used in this study did not trigger any important inflammatory reaction; (2) when the sinus membrane was elevated, bone formation was a constant finding; and (3) an ideal space-making device should be stable and elevate the membrane to ensure a maintained connection between the membrane and the secluded space.

Bjorling problem for timelike surfaces in the Lorentz-Minkowski space

We solve the Bjorling problem for timelike surfaces in the Lorentz-Minkowski space through a split-complex representation formula obtained for this kind of surface. Our approach uses the split-complex numbers and natural split-holomorphic extensions. As applications, we show that the minimal timelike surfaces of revolution as well as minimal ruled timelike surfaces can be characterized as solutions of certain adequate Bjorling problems in the Lorentz-Minkowski space. (C) 2010 Elsevier Inc. All rights reserved.

Dynamic whole system testing of combined renewable heating systems : the current state of the art

Objective: For the evaluation of the energetic performance of combined renewable heating systems that supply space heat and domestic hot water for single family houses, dynamic behaviour, component interactions, and control of the system play a crucial role and should be included in test methods. Methods: New dynamic whole system test methods were developed based on “hardware in the loop” concepts. Three similar approaches are described and their differences are discussed. The methods were applied for testing solar thermal systems in combination with fossil fuel boilers (heating oil and natural gas), biomass boilers, and/or heat pumps. Results: All three methods were able to show the performance of combined heating systems under transient operating conditions. The methods often detected unexpected behaviour of the tested system that cannot be detected based on steady state performance tests that are usually applied to single components. Conclusion: Further work will be needed to harmonize the different test methods in order to reach comparable results between the different laboratories. Practice implications: A harmonized approach for whole system tests may lead to new test standards and improve the accuracy of performance prediction as well as reduce the need for field tests.

Space Window*: an exhibition devoted to Man in Space and Space in Art

Call for submissions to participate in a show in the Bell Gallery, List Art Building at Brown University. Co-sponsored by the RI State Council on the Arts and the Providence Inner City Arts Association.

Alcoolemia em vítimas fatais de acidentes de trânsito no Rio Grande do Norte empregando cromatografia em fase gasosa - head space

Alcohol is one of the few psychotropic drugs that their consumption has admitted legally and sometimes encouraged by the society. Studies show alcohol as the highest consumption of drugs among young people and society in general, probably because of its availability and easy access. The abuse causes public health problems, which was closely related to the violence, socioeconomic problems and the high number of automobile accidents. Transit is one of the main sectors affected by the effects of alcohol, observing a high incidence in the studies. About half of automobile accidents occurs after the consumption of alcoholic beverage, and the vast majority of cases related to high concentrations of alcohol in the bloodstream. The relationship of drunk with traffic accidents is in fact evident everywhere in the world, including Brazil, where studies have shown a high relationship between alcohol consumption and traffic accidents. This study determined the alcohol in fatal victims of traffic accidents in the state of Rio Grande do Norte and established the profile of this population compared with those found in Brazil and other countries. Samples of blood of ethanol added to fulfillment of the standardization of chromatographic conditions and procedures for the analysis, being employed in the determination of alcohol in blood samples of 277 victims of traffic accidents, collected at the Institute of Scientific Technical Police of Rio Grande do North (ITEP) in the year 2007. The blood alcohol level was determined in these samples correlated with the sex, age and marital status of the victim and the location, day of week and month when the accident occurred, is doing a statistical analysis and outlining a profile of the victims of an accident at transit in the state of Rio Grande do Norte. The parameters of standardization studied ensured the quality of the analytical method and, consequently, to obtain reliable laboratory results. Being given the best temperature for injector (150 ºC), detector (250 ºC) and column (50 ºC) with a flow of gas in the column of 2mL/minutos and analysis of time of 12 minutes. The method was linear in the range of 0.01 to 3.2 g / L (r2 = 0.9989) with average recovery of 100.2% and precision with coefficient of variation less than 15%. The analysis carried out on victims of fatal road traffic accidents, ethanol detected in the blood in 66.43% of the victims and these, 96% showed concentration ≥ 0.2 g / L, 87.73% of victims were male, while 12.27% female. The younger age group (1535 years) was the most involved (52,35%) and most single (55.60%). The accidents occurred with greater prevalence in the day on Monday (27%) followed by Sunday (24,19%) and Saturday (15,52%) and it was found that the prevalence of injuries varied between the different months of the year, and in February (14.4%) and April (10.47%) the months that had a higher number of accidents, however this oscillation showed no statistically significant difference. Also no significant difference was observed between the tracks of concentration found in men and women. The standardized method showed to be efficient, given satisfactorily to the goals of this work, and the high levels of alcohol found in victims of fatal road traffic accidents are consistent with several studies of literature, and the profile of the victim also supported by presenting in its most young adults, male and single

Genetic diversity of natural populations of Anemopaegma arvense (Bignoniaceae) in the Cerrado of São Paulo State, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A new method for real time computation of power quality indices based on instantaneous space phasors

One of the important issues about using renewable energy is the integration of dispersed generation in the distribution networks. Previous experience has shown that the integration of dispersed generation can improve voltage profile in the network, decrease loss, etc. but can create safety and technical problems as well. This work report the application of the instantaneous space phasors and the instantaneous complex power in observing performances of the distribution networks with dispersed generators in steady state. New IEEE apparent power definition, the so-called Buchholz-Goodhue effective apparent power, as well as new proposed power quality (oscillation) index in the three-phase distribution systems with unbalanced loads and dispersed generators, are applied. Results obtained from several case studies using IEEE 34 nodes test network are presented and discussed. (C) 2006 Elsevier B.V. All rights reserved.

Configuration-dependent diffusion can shift the kinetic transition state and barrier height of protein folding

We show that diffusion can play an important role in protein-folding kinetics. We explicitly calculate the diffusion coefficient of protein folding in a lattice model. We found that diffusion typically is configuration- or reaction coordinate-dependent. The diffusion coefficient is found to be decreasing with respect to the progression of folding toward the native state, which is caused by the collapse to a compact state constraining the configurational space for exploration. The configuration- or position-dependent diffusion coefficient has a significant contribution to the kinetics in addition to the thermodynamic free-energy barrier. It effectively changes (increases in this case) the kinetic barrier height as well as the position of the corresponding transition state and therefore modifies the folding kinetic rates as well as the kinetic routes. The resulting folding time, by considering both kinetic diffusion and the thermodynamic folding free-energy profile, thus is slower than the estimation from the thermodynamic free-energy barrier with constant diffusion but is consistent with the results from kinetic simulations. The configuration- or coordinate-dependent diffusion is especially important with respect to fast folding, when there is a small or no free-energy barrier and kinetics is controlled by diffusion. Including the configurational dependence will challenge the transition state theory of protein folding. The classical transition state theory will have to be modified to be consistent. The more detailed folding mechanistic studies involving phi value analysis based on the classical transition state theory also will have to be modified quantitatively.

Low-energy three-body charge transfer reactions with Coulomb interaction in the final state

Three-body charge transfer reactions with Coulomb interaction in the final state are considered within the framework of coordinate-space integro-differential Faddeev-Hahn-type equations within two- and six-state close-coupling approximations. The method is employed to study direct muon transfer in low-energy collisions of the muonic hydrogen H-mu by helium (He2+) and lithium (Li3+) nuclei. The experimentally observed isotopic dependence is reproduced.

Low-energy direct muon transfer from H to Ne10+, S16+ and Ar18+ using the two-state close-coupling approximation to the Faddeev-Hahn-type equation

We perform a three-body calculation of direct muon-transfer rates from thermalized muonic hydrogen isotopes to bare nuclei Ne10+, S16+ and Ar18+ employing integro-differential Faddeev-Hahn-type equations in configuration space with a two-state close-coupling approximation scheme. All Coulomb potentials including the strong final-state Coulomb repulsion are treated exactly. A long-range polarization potential is included in the elastic channel to take into account the high polarizability of the muonic hydrogen. The transfer rates so-calculated are in good agreement with recent experiments. We find that the muon is captured predominantly in the n = 6, 9 and 10 states of muonic Ne10+, S16+ and Ar18+, respectively.

Spacelike and timelike pion electromagnetic form factor and Fock state components within the light-front dynamics

The simultaneous investigation of the pion electromagnetic form factor in the space- and timelike regions within a light-front model allows one to address the issue of nonvalence components of the pion and photon wave functions. Our relativistic approach is based on a microscopic vector-meson-dominance model for the dressed vertex where a photon decays in a quark-antiquark pair, and on a simple parametrization for the emission or absorption of a pion by a quark. The results show an excellent agreement in the space like region up to -10 (GeV/c)(2), while in timelike region the model produces reasonable results up to 10 (GeV/c)(2).

Solutions of the bound-state Faddeev-Yakubovsky equations in three dimensions by using NN and 3N potential models

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)