Twin experiments with the equivalent weights particle filter and HadCM3

This paper investigates the use of a particle filter for data assimilation with a full scale coupled ocean–atmosphere general circulation model. Synthetic twin experiments are performed to assess the performance of the equivalent weights filter in such a high-dimensional system. Artificial 2-dimensional sea surface temperature fields are used as observational data every day. Results are presented for different values of the free parameters in the method. Measures of the performance of the filter are root mean square errors, trajectories of individual variables in the model and rank histograms. Filter degeneracy is not observed and the performance of the filter is shown to depend on the ability to keep maximum spread in the ensemble.

Sparse space-time boundary element methods for the heat equation

The goal of this work is the efficient solution of the heat equation with Dirichlet or Neumann boundary conditions using the Boundary Elements Method (BEM). Efficiently solving the heat equation is useful, as it is a simple model problem for other types of parabolic problems. In complicated spatial domains as often found in engineering, BEM can be beneficial since only the boundary of the domain has to be discretised. This makes BEM easier than domain methods such as finite elements and finite differences, conventionally combined with time-stepping schemes to solve this problem. The contribution of this work is to further decrease the complexity of solving the heat equation, leading both to speed gains (in CPU time) as well as requiring smaller amounts of memory to solve the same problem. To do this we will combine the complexity gains of boundary reduction by integral equation formulations with a discretisation using wavelet bases. This reduces the total work to O(h

A Variance-Minimizing Filter for Large-Scale Applications

A truly variance-minimizing filter is introduced and its per for mance is demonstrated with the Korteweg– DeV ries (KdV) equation and with a multilayer quasigeostrophic model of the ocean area around South Africa. It is recalled that Kalman-like filters are not variance minimizing for nonlinear model dynamics and that four - dimensional variational data assimilation (4DV AR)-like methods relying on per fect model dynamics have dif- ficulty with providing error estimates. The new method does not have these drawbacks. In fact, it combines advantages from both methods in that it does provide error estimates while automatically having balanced states after analysis, without extra computations. It is based on ensemble or Monte Carlo integrations to simulate the probability density of the model evolution. When obser vations are available, the so-called importance resampling algorithm is applied. From Bayes’ s theorem it follows that each ensemble member receives a new weight dependent on its ‘ ‘distance’ ’ t o the obser vations. Because the weights are strongly var ying, a resampling of the ensemble is necessar y. This resampling is done such that members with high weights are duplicated according to their weights, while low-weight members are largely ignored. In passing, it is noted that data assimilation is not an inverse problem by nature, although it can be for mulated that way . Also, it is shown that the posterior variance can be larger than the prior if the usual Gaussian framework is set aside. However , i n the examples presented here, the entropy of the probability densities is decreasing. The application to the ocean area around South Africa, gover ned by strongly nonlinear dynamics, shows that the method is working satisfactorily . The strong and weak points of the method are discussed and possible improvements are proposed.

Analysis Scheme in the Ensemble Kalman Filter

This paper discusses an important issue related to the implementation and interpretation of the analysis scheme in the ensemble Kalman filter . I t i s shown that the obser vations must be treated as random variables at the analysis steps. That is, one should add random perturbations with the correct statistics to the obser vations and generate an ensemble of obser vations that then is used in updating the ensemble of model states. T raditionally , this has not been done in previous applications of the ensemble Kalman filter and, as will be shown, this has resulted in an updated ensemble with a variance that is too low . This simple modification of the analysis scheme results in a completely consistent approach if the covariance of the ensemble of model states is interpreted as the prediction error covariance, and there are no further requirements on the ensemble Kalman filter method, except for the use of an ensemble of sufficient size. Thus, there is a unique correspondence between the error statistics from the ensemble Kalman filter and the standard Kalman filter approach

Assimilation of Geosat Altimeter Data for the Agulhas Current Using the Ensemble Kalman Filter with a Quasigeostrophic Model

The ring-shedding process in the Agulhas Current is studied using the ensemble Kalman filter to assimilate geosat altimeter data into a two-layer quasigeostrophic ocean model. The properties of the ensemble Kalman filter are further explored with focus on the analysis scheme and the use of gridded data. The Geosat data consist of 10 fields of gridded sea-surface height anomalies separated 10 days apart that are added to a climatic mean field. This corresponds to a huge number of data values, and a data reduction scheme must be applied to increase the efficiency of the analysis procedure. Further, it is illustrated how one can resolve the rank problem occurring when a too large dataset or a small ensemble is used.

Noether-type discrete conserved quantities arising from a finite element approximation of a variational problem

In this work, we prove a weak Noether-type Theorem for a class of variational problems that admit broken extremals. We use this result to prove discrete Noether-type conservation laws for a conforming finite element discretisation of a model elliptic problem. In addition, we study how well the finite element scheme satisfies the continuous conservation laws arising from the application of Noether’s first theorem (1918). We summarise extensive numerical tests, illustrating the conservation of the discrete Noether law using the p-Laplacian as an example and derive a geometric-based adaptive algorithm where an appropriate Noether quantity is the goal functional.

Thyroid hormone can increase estrogen-mediated transcription from a consensus estrogen response element in neuroblastoma cells

Thyroid hormones (T) and estrogens (E) are nuclear receptor ligands with at least two molecular mechanisms of action: (i) relatively slow genomic effects, such as the regulation of transcription by cognate T receptors (TR) and E receptors (ER); and (ii) relatively rapid nongenomic effects, such as kinase activation and calcium release initiated at the membrane by putative membrane receptors. Genomic and nongenomic effects were thought to be disparate and independent. However, in a previous study using a two-pulse paradigm in neuroblastoma cells, we showed that E acting at the membrane could potentiate transcription from an E-driven reporter gene in the nucleus. Because both T and E can have important effects on mood and cognition, it is possible that the two hormones can act synergistically. In this study, we demonstrate that early actions of T via TRalpha1 and TRbeta1 can potentiate E-mediated transcription (genomic effects) from a consensus E response element (ERE)-driven reporter gene in transiently transfected neuroblastoma cells. Such potentiation was reduced by inhibition of mitogen-activated protein kinase. Using phosphomutants of ERalpha, we also show that probable mitogen-activated protein kinase phosphorylation sites on the ERalpha, the serines at position 167 and 118, are important in TRbeta1-mediated potentiation of ERalpha-induced transactivation. We suggest that crosstalk between T and E includes potential interactions through both nuclear and membrane-initiated molecular mechanisms of hormone signaling.

Implicit equal-weights particle filter

Filter degeneracy is the main obstacle for the implementation of particle filter in non-linear high-dimensional models. A new scheme, the implicit equal-weights particle filter (IEWPF), is introduced. In this scheme samples are drawn implicitly from proposal densities with a different covariance for each particle, such that all particle weights are equal by construction. We test and explore the properties of the new scheme using a 1,000-dimensional simple linear model, and the 1,000-dimensional non-linear Lorenz96 model, and compare the performance of the scheme to a Local Ensemble Kalman Filter. The experiments show that the new scheme can easily be implemented in high-dimensional systems and is never degenerate, with good convergence properties in both systems.

NEONATAL HANDLING AND THE MATERNAL ODOR PREFERENCE IN RAT PUPS: INVOLVEMENT OF MONOAMINES AND CYCLIC AMP RESPONSE ELEMENT-BINDING PROTEIN PATHWAY IN THE OLFACTORY BULB

Early-life environmental events, such as the handling procedure, can induce long-lasting alterations upon several behavioral and neuroendocrine systems. However, the changes within the pups that could be causally related to the effects in adulthood are still poorly understood. In the present study, we analyzed the effects of neonatal handling on behavioral (maternal odor preference) and biochemical (cyclic AMP response element-binding protein (CREB) phosphorylation, noradrenaline (NA), and serotonin (5-HT) levels in the olfactory bulb (OB)) parameters in 7-day-old male and female rat pups. Repeated handling (RH) abolished preference for the maternal odor in female pups compared with nonhandled (NH) and the single-handled (SH) ones, while in RH males the preference was not different than NH and SH groups. In both male and female pups, RH decreased NA activity in the OB, but 5-HT activity increased only in males. Since preference for the maternal odor involves the synergic action of NA and 5-HT in the OB, the maintenance of the behavior in RH males could be related to the increased 5-HT activity, in spite of reduction in the NA activity in the OB. RH did not alter CREB phosphorylation in the OB of both male and females compared with NH pups. The repeated handling procedure can affect the behavior of rat pups in response to the maternal odor and biochemical parameters related to the olfactory learning mechanism. Sex differences were already detected in 7-day-old pups. Although the responsiveness of the hypothalamic-pituitary-adrenal axis to stressors is reduced in the neonatal period, environmental interventions may impact behavioral and biochemical mechanisms relevant to the animal at that early age. (C) 2009 IBRO. Published by Elsevier Ltd. All rights reserved.

Characterization of Metal and Trace Element Contents of Particulate Matter (PM10) Emitted by Vehicles Running on Brazilian FuelsHydrated Ethanol and Gasoline with 22% of Anhydrous Ethanol

Emission of fine particles by mobile sources has been a matter of great concern due to its potential risk both to human health and the environment. Although there is no evidence that one sole component may be responsible for the adverse health outcomes, it is postulated that the metal particle content is one of the most important factors, mainly in relation to oxidative stress. Data concerning the amount and type of metal particles emitted by automotive vehicles using Brazilian fuels are limited. The aim of this study was to identify inhalable particles (PM10) and their trace metal content in two light-duty vehicles where one was fueled with ethanol while the other was fueled with gasoline mixed with 22% of anhydrous ethanol (gasohol); these engines were tested on a chassis dynamometer. The elementary composition of the samples was evaluated by the particle-induced x-ray emission technique. The experiment showed that total emission factors ranged from 2.5 to 11.8 mg/km in the gasohol vehicle, and from 1.2 to 3 mg/km in the ethanol vehicle. The majority of particles emitted were in the fine fraction (PM2.5), in which Al, Si, Ca, and Fe corresponded to 80% of the total weight. PM10 emissions from the ethanol vehicle were about threefold lower than those of gasohol. The elevated amount of fine particulate matter is an aggravating factor, considering that these particles, and consequently associated metals, readily penetrate deeply into the respiratory tract, producing damage to lungs and other tissues.

Reappraisal of the effective elastic thickness for the sub-Andes using 3-D finite element flexural modelling, gravity and geological constraints

P>Estimates of effective elastic thickness (T(e)) for the western portion of the South American Plate using, independently, forward flexural modelling and coherence analysis, suggest different thermomechanical properties for the same continental lithosphere. We present a review of these T(e) estimates and carry out a critical reappraisal using a common methodology of 3-D finite element method to solve a differential equation for the bending of a thin elastic plate. The finite element flexural model incorporates lateral variations of T(e) and the Andes topography as the load. Three T(e) maps for the entire Andes were analysed: Stewart & Watts (1997), Tassara et al. (2007) and Perez-Gussinye et al. (2007). The predicted flexural deformation obtained for each T(e) map was compared with the depth to the base of the foreland basin sequence. Likewise, the gravity effect of flexurally induced crust-mantle deformation was compared with the observed Bouguer gravity. T(e) estimates using forward flexural modelling by Stewart & Watts (1997) better predict the geological and gravity data for most of the Andean system, particularly in the Central Andes, where T(e) ranges from greater than 70 km in the sub-Andes to less than 15 km under the Andes Cordillera. The misfit between the calculated and observed foreland basin subsidence and the gravity anomaly for the Maranon basin in Peru and the Bermejo basin in Argentina, regardless of the assumed T(e) map, may be due to a dynamic topography component associated with the shallow subduction of the Nazca Plate beneath the Andes at these latitudes.

The Diffusion Kernel Filter

A particle filter method is presented for the discrete-time filtering problem with nonlinear ItA ` stochastic ordinary differential equations (SODE) with additive noise supposed to be analytically integrable as a function of the underlying vector-Wiener process and time. The Diffusion Kernel Filter is arrived at by a parametrization of small noise-driven state fluctuations within branches of prediction and a local use of this parametrization in the Bootstrap Filter. The method applies for small noise and short prediction steps. With explicit numerical integrators, the operations count in the Diffusion Kernel Filter is shown to be smaller than in the Bootstrap Filter whenever the initial state for the prediction step has sufficiently few moments. The established parametrization is a dual-formula for the analysis of sensitivity to gaussian-initial perturbations and the analysis of sensitivity to noise-perturbations, in deterministic models, showing in particular how the stability of a deterministic dynamics is modeled by noise on short times and how the diffusion matrix of an SODE should be modeled (i.e. defined) for a gaussian-initial deterministic problem to be cast into an SODE problem. From it, a novel definition of prediction may be proposed that coincides with the deterministic path within the branch of prediction whose information entropy at the end of the prediction step is closest to the average information entropy over all branches. Tests are made with the Lorenz-63 equations, showing good results both for the filter and the definition of prediction.

Copia Retrotransposon in the Zaprionus Genus: Another Case of Transposable Element Sharing with the Drosophila melanogaster Subgroup

Copia is a retrotransposon that appears to be distributed widely among the Drosophilidae subfamily. Evolutionary analyses of regulatory regions have indicated that the Copia retrotransposon evolved through both positive and purifying selection, and that horizontal transfer (HT) could also explain its patchy distribution of the among the subfamilies of the melanogaster subgroup. Additionally, Copia elements could also have transferred between melanogaster subgroup and other species of Drosophilidae-D. willistoni and Z. tuberculatus. In this study, we surveyed seven species of the Zaprionus genus by sequencing the LTR-ULR and reverse transcriptase regions, and by using RT-PCR in order to understand the distribution and evolutionary history of Copia in the Zaprionus genus. The Copia element was detected, and was transcriptionally active, in all species investigated. Structural and selection analysis revealed Zaprionus elements to be closely related to the most ancient subfamily of the melanogaster subgroup, and they seem to be evolving mainly under relaxed purifying selection. Taken together, these results allowed us to classify the Zaprionus sequences as a new subfamily-ZapCopia, a member of the Copia retrotransposon family of the melanogaster subgroup. These findings indicate that the Copia retrotransposon is an ancient component of the genomes of the Zaprionus species and broaden our understanding of the diversity of retrotransposons in the Zaprionus genus.