The effect of the equivalent-weights particle filter on dynamical balance in a primitive equation model

The disadvantage of the majority of data assimilation schemes is the assumption that the conditional probability density function of the state of the system given the observations [posterior probability density function (PDF)] is distributed either locally or globally as a Gaussian. The advantage, however, is that through various different mechanisms they ensure initial conditions that are predominantly in linear balance and therefore spurious gravity wave generation is suppressed. The equivalent-weights particle filter is a data assimilation scheme that allows for a representation of a potentially multimodal posterior PDF. It does this via proposal densities that lead to extra terms being added to the model equations and means the advantage of the traditional data assimilation schemes, in generating predominantly balanced initial conditions, is no longer guaranteed. This paper looks in detail at the impact the equivalent-weights particle filter has on dynamical balance and gravity wave generation in a primitive equation model. The primary conclusions are that (i) provided the model error covariance matrix imposes geostrophic balance, then each additional term required by the equivalent-weights particle filter is also geostrophically balanced; (ii) the relaxation term required to ensure the particles are in the locality of the observations has little effect on gravity waves and actually induces a reduction in gravity wave energy if sufficiently large; and (iii) the equivalent-weights term, which leads to the particles having equivalent significance in the posterior PDF, produces a change in gravity wave energy comparable to the stochastic model error. Thus, the scheme does not produce significant spurious gravity wave energy and so has potential for application in real high-dimensional geophysical applications.

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.

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.

Ecuador Paraiso Escondido virus, a new flavivirus isolated from New World sand flies in Ecuador, is the first representative of a novel clade in the genus flavivirus

A new flavivirus, Ecuador Paraiso Escondido virus (EPEV), named after the village where it was discovered, was isolated from sand flies (Psathyromyia abonnenci, formerly Lutzomyia abonnenci) that are unique to the New World. This represents the first sand fly-borne flavivirus identified in the New World. EPEV exhibited a typical flavivirus genome organization. Nevertheless, the maximum pairwise amino acid sequence identity with currently recognized flaviviruses was 52.8%. Phylogenetic analysis of the complete coding sequence showed that EPEV represents a distinct clade which diverged from a lineage that was ancestral to the nonvectored flaviviruses Entebbe bat virus, Yokose virus, and Sokoluk virus and also the Aedes-associated mosquito-borne flaviviruses, which include yellow fever virus, Sepik virus, Saboya virus, and others. EPEV replicated in C6/36 mosquito cells, yielding high infectious titers, but failed to reproduce either in vertebrate cell lines (Vero, BHK, SW13, and XTC cells) or in suckling mouse brains. This surprising result, which appears to eliminate an association with vertebrate hosts in the life cycle of EPEV, is discussed in the context of the evolutionary origins of EPEV in the New World.The flaviviruses are rarely (if ever) vectored by sand fly species, at least in the Old World. We have identified the first representative of a sand fly-associated flavivirus, Ecuador Paraiso Escondido virus (EPEV), in the New World. EPEV constitutes a novel clade according to current knowledge of the flaviviruses. Phylogenetic analysis of the virus genome showed that EPEV roots the Aedes-associated mosquito-borne flaviviruses, including yellow fever virus. In light of this new discovery, the New World origin of EPEV is discussed together with that of the other flaviviruses.

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.

Species-specific identification of Leishmania in naturally infected sand flies captured in Mato Grosso do Sul State, Brazil

An increase in cutaneous and visceral leishmaniasis cases has been reported in recent years in the state of Mato Grosso do Sul, Brazil, and little is known to date about their etiological agents. An investigation into natural Leishmania infection of sand flies captured in this state between December 2003 and August 2004 was carried out. Mini-exon sequences were used as targets to identify Leishmania, and an RFLP technique was employed for those identified as belonging to the Viannia subgenus. Calculation of the minimal infection rate (MR) revealed that 1.6% of sand flies captured in the forest, peridomicile and intradomicile were positive. Six species were found to be infected by Leishmania (V.) braziliensis. Interestingly, two of the six species. Lutzomyia longipalpis and Nyssomyia whitmani, were captured in anthropic environments. The findings of this study constitute a useful tool for planning control measures against this disease in the State of Mato Grosso do Sul. (C) 2010 Elsevier B.V. All rights reserved.

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.

Wolbachia in Anastrepha Fruit Flies (Diptera: Tephritidae)

Endosymbiotic bacteria of the genus Wolbachia are widespread among arthropods and cause a variety of reproductive abnormalities, such as cytoplasmic incompatibility, thelytokous parthenogenesis, male-killing, and host feminization. In this study, we used three sets of Wolbachia-specific primers (16S rDNA, ftsZ, and wsp) in conjunction with the polymerase chain reaction (PCR), cloning and sequencing to study the infection of fruit flies (Anastrepha spp. and Ceratitis capitata) by Wolbachia. The flies were collected at several localities in Brazil and at Guayaquil, Ecuador. All of the fruit flies studied were infected with Wolbachia supergroup A, in agreement with the high prevalence of this group in South America. Phylogenetic analysis showed that the wsp gene was the most sensitive gene for studying the relationships among Wolbachia strains. The Wolbachia sequences detected in these fruit flies were similar to those such as wMel reported for other fruit flies. These results show that the infection of Anastrepha fruit flies by Wolbachia is much more widespread than previously thought.

New genera of philopotine spider flies (Diptera, Acroceridae) with a key to living and fossil genera

In this paper we describe two new genera of philopotine Acroceridae: Schlingeriella irwini gen. et sp. n. (New Caledonia) and Quasi fisheri gen. et sp. n. (Mexico). The Baltic amber species Eulonchiella eocenica Meunier, 1912 is rediagnosed and a neotype designated based on a newly discovered specimen. We also provide a dichotomous key to the world genera of Philopotinae, both living and fossil.

The Muscini flies of the world (Diptera, Muscidae): identification key and generic diagnoses

Flies of the tribe Muscini (Diptera, Muscidae) are worldwide in distribution and are represented by some 350 species in 18 genera. The present study provides an identification key and diagnoses for all the genera of world Muscini: Biopyrellia Townsend, Curranosia Paterson, Dasyphora Robineau-Desvoidy, Deltotus Seguy, Hennigmyia Peris, Mesembrina Meigen, Mitroplatia Enderlein, Morellia Robineau-Desvoidy, Musca Linnaeus, Myiophaea Enderlein, Neomyia Walker, Neorypellia Pont, Polietes Rondani, Polietina Schnabl & Dziedzicki, Pyrellia Robineau-Desvoidy, Pyrellina Malloch, Sarcopromusca Townsend, Ziminellia Nihei & de Carvalho. Most infrageneric taxa are also represented, namely, the sub-genera of Dasyphora and Morellia. Comments on phylogeny support (whenever pertinent) and the major references containing revisions and regional identification keys to species are provided for each genus and subgenus.

New Trypanosoma (Duttonella) vivax genotypes from tsetse flies in East Africa

Salivarian trypanosomes pose a substantial threat to livestock, but their full diversity is not known. To survey trypanosomes carried by tsetse in Tanzania, DNA samples from infected proboscides of Glossina pallidipes and G. swynnertoni were identified using fluorescent fragment length barcoding (FFLB), which discriminates species by size polymorphisms in multiple regions of the ribosomal RNA locus. FELLB identified the trypanosomes in 65 of 105 (61.9%) infected proboscides, revealing 9 mixed infections. Of 7 different FFLB profiles, 2 were similar but not identical to reference West African Trypanosoma vivax; 5 other profiles belonged to known species also identified in fly midguts. Phylogenetic analysis of the glycosomal glyceraldehyde phosphate dehydrogenase gene revealed that the Tanzanian T. vivax samples fell into 2 distinct groups, both outside the main chide of African and South American T. vivax. These new T. vivax genotypes were common and widespread in tsetse in Tanzania. The T. brucei-like trypanosome previously described from tsetse midguts was also found in 2 proboscides, demonstrating a salivarian transmission route. Investigation of mammalian host range and pathogenicity will reveal the importance of these new trypanosomes for the epidemiology and control of animal trypanosomiasis in East Africa.