Conditioning of the weak-constraint variational data assimilation problem for numerical weather prediction

4-Dimensional Variational Data Assimilation (4DVAR) assimilates observations through the minimisation of a least-squares objective function, which is constrained by the model flow. We refer to 4DVAR as strong-constraint 4DVAR (sc4DVAR) in this thesis as it assumes the model is perfect. Relaxing this assumption gives rise to weak-constraint 4DVAR (wc4DVAR), leading to a different minimisation problem with more degrees of freedom. We consider two wc4DVAR formulations in this thesis, the model error formulation and state estimation formulation. The 4DVAR objective function is traditionally solved using gradient-based iterative methods. The principle method used in Numerical Weather Prediction today is the Gauss-Newton approach. This method introduces a linearised `inner-loop' objective function, which upon convergence, updates the solution of the non-linear `outer-loop' objective function. This requires many evaluations of the objective function and its gradient, which emphasises the importance of the Hessian. The eigenvalues and eigenvectors of the Hessian provide insight into the degree of convexity of the objective function, while also indicating the difficulty one may encounter while iterative solving 4DVAR. The condition number of the Hessian is an appropriate measure for the sensitivity of the problem to input data. The condition number can also indicate the rate of convergence and solution accuracy of the minimisation algorithm. This thesis investigates the sensitivity of the solution process minimising both wc4DVAR objective functions to the internal assimilation parameters composing the problem. We gain insight into these sensitivities by bounding the condition number of the Hessians of both objective functions. We also precondition the model error objective function and show improved convergence. We show that both formulations' sensitivities are related to error variance balance, assimilation window length and correlation length-scales using the bounds. We further demonstrate this through numerical experiments on the condition number and data assimilation experiments using linear and non-linear chaotic toy models.

Mesh adaptation on the sphere using optimal transport and the numerical solution of a Monge-Ampère type equation

An equation of Monge-Ampère type has, for the first time, been solved numerically on the surface of the sphere in order to generate optimally transported (OT) meshes, equidistributed with respect to a monitor function. Optimal transport generates meshes that keep the same connectivity as the original mesh, making them suitable for r-adaptive simulations, in which the equations of motion can be solved in a moving frame of reference in order to avoid mapping the solution between old and new meshes and to avoid load balancing problems on parallel computers. The semi-implicit solution of the Monge-Ampère type equation involves a new linearisation of the Hessian term, and exponential maps are used to map from old to new meshes on the sphere. The determinant of the Hessian is evaluated as the change in volume between old and new mesh cells, rather than using numerical approximations to the gradients. OT meshes are generated to compare with centroidal Voronoi tesselations on the sphere and are found to have advantages and disadvantages; OT equidistribution is more accurate, the number of iterations to convergence is independent of the mesh size, face skewness is reduced and the connectivity does not change. However anisotropy is higher and the OT meshes are non-orthogonal. It is shown that optimal transport on the sphere leads to meshes that do not tangle. However, tangling can be introduced by numerical errors in calculating the gradient of the mesh potential. Methods for alleviating this problem are explored. Finally, OT meshes are generated using observed precipitation as a monitor function, in order to demonstrate the potential power of the technique.

Conditioning and preconditioning of the optimal state estimation problem

Optimal state estimation is a method that requires minimising a weighted, nonlinear, least-squares objective function in order to obtain the best estimate of the current state of a dynamical system. Often the minimisation is non-trivial due to the large scale of the problem, the relative sparsity of the observations and the nonlinearity of the objective function. To simplify the problem the solution is often found via a sequence of linearised objective functions. The condition number of the Hessian of the linearised problem is an important indicator of the convergence rate of the minimisation and the expected accuracy of the solution. In the standard formulation the convergence is slow, indicating an ill-conditioned objective function. A transformation to different variables is often used to ameliorate the conditioning of the Hessian by changing, or preconditioning, the Hessian. There is only sparse information in the literature for describing the causes of ill-conditioning of the optimal state estimation problem and explaining the effect of preconditioning on the condition number. This paper derives descriptive theoretical bounds on the condition number of both the unpreconditioned and preconditioned system in order to better understand the conditioning of the problem. We use these bounds to explain why the standard objective function is often ill-conditioned and why a standard preconditioning reduces the condition number. We also use the bounds on the preconditioned Hessian to understand the main factors that affect the conditioning of the system. We illustrate the results with simple numerical experiments.

Non-bee insects are important contributors to global crop pollination

Wild and managed bees are well documented as effective pollinators of global crops of economic importance. However, the contributions by pollinators other than bees have been little explored despite their potential to contribute to crop production and stability in the face of environmental change. Non-bee pollinators include flies, beetles, moths, butterflies, wasps, ants, birds, and bats, among others. Here we focus on non-bee insects and synthesize 39 field studies from five continents that directly measured the crop pollination services provided by non-bees, honey bees, and other bees to compare the relative contributions of these taxa. Non-bees performed 25–50% of the total number of flower visits. Although non-bees were less effective pollinators than bees per flower visit, they made more visits; thus these two factors compensated for each other, resulting in pollination services rendered by non-bees that were similar to those provided by bees. In the subset of studies that measured fruit set, fruit set increased with non-bee insect visits independently of bee visitation rates, indicating that non-bee insects provide a unique benefit that is not provided by bees. We also show that non-bee insects are not as reliant as bees on the presence of remnant natural or seminatural habitat in the surrounding landscape. These results strongly suggest that non-bee insect pollinators play a significant role in global crop production and respond differently than bees to landscape structure, probably making their crop pollination services more robust to changes in land use. Non-bee insects provide a valuable service and provide potential insurance against bee population declines.

Plastic and evolutionary responses to heat stress in a temperate dung fly: negative correlation between basal and induced heat tolerance?

Extreme weather events such as heat waves are becoming more frequent and intense. Populations can cope with elevated heat stress by evolving higher basal heat tolerance (evolutionary response) and/or stronger induced heat tolerance (plastic response). However, there is ongoing debate about whether basal and induced heat tolerance are negatively correlated and whether adaptive potential in heat tolerance is sufficient under ongoing climate warming. To evaluate the evolutionary potential of basal and induced heat tolerance, we performed experimental evolution on a temperate source 4 population of the dung fly Sepsis punctum. Offspring of flies adapted to three thermal selection regimes (Hot, Cold and Reference) were subjected to acute heat stress after having been exposed to either a hot-acclimation or non-acclimation pretreatment. As different traits may respond differently to temperature stress, several physiological and life history traits were assessed. Condition dependence of the response was evaluated by exposing juveniles to different levels of developmental (food restriction/rearing density) stress. Heat knockdown times were highest, whereas acclimation effects were lowest in the Hot selection regime, indicating a negative association between basal and induced heat tolerance. However, survival, adult longevity, fecundity and fertility did not show such a pattern. Acclimation had positive effects in heat-shocked flies, but in the absence of heat stress hot-acclimated flies had reduced life spans relative to nonacclimated ones, thereby revealing a potential cost of acclimation. Moreover, body size positively affected heat tolerance and unstressed individuals were less prone to heat stress than stressed flies, offering support for energetic costs associated with heat tolerance. Overall, our results indicate that heat tolerance of temperate insects can evolve under rising temperatures, but this response could be limited by a negative relationship between basal and induced thermotolerance, and may involve some but not other fitness-related traits.

Detection of Leishmania (Leishmania) infantum RNA in fleas and ticks collected from naturally infected dogs

The occurrence of the insect vector (sand flies) with low rates of Leishmania infection, as well as autochthonous transmission in the absence of the natural vector in dogs, have been reported. These unexpected data suggest a hypothesis of other arthropods as a possible way of Leishmania transmission. The prevalence of Leishmania (Leishmania) infantum in fleas and ticks collected from dogs with canine visceral leishmaniasis (CVL), as well as parasite viability, were evaluated herein. The presence of L. (L.) infantum was assayed by PCR and ELISA in ectoparasites and biological samples from 73 dogs living in a Brazilian endemic area. As the occurrence of Leishmania DNA in ticks and fleas is expected given their blood-feeding habits, we next investigated whether parasites can remain viable inside ticks. PCR and ELISA confirmed that 83% of the dogs had CVL. Fleas and ticks (nymphs, male and female adults) were collected in 55% and 63% of the 73 dogs, respectively. Out of the 60 dogs with CVL, 80% harbored ectoparasites infected with L. (L.) infantum. The infection rates of the ectoparasites were 23% and 50% for fleas and ticks, respectively. The RNA analysis of the extract from ticks left in laboratory conditions during 7 to 10 days after removal from CVL dogs showed that parasites were alive. In addition, live parasites were also detected inside adult ticks recently molted in laboratory conditions. These findings indicate a higher infection rate of L. (L.) infantum in ticks and fleas, but they do not conclusively demonstrate whether these ticks can act as vectors of CVL, despite the fact that their rates were higher than those previously described in Lutzomyia longipalpis. The presence of viable L. (L.) infantum in ticks suggests the possible importance of dog ectoparasites in CVL dissemination.

Bilateral anomaly in the style of Sciopemyia sordellii (Shannon & Del Ponte) (Diptera, Psychodidae)

This paper presents a male specimen of Sciopemyia sordellii with a rare bilateral anomaly, consisting in eight spines in a style and five in the other. This species has four spines in each style as its normal number. The specimen was captured using a CDC light trap, in a forested area in the State Park ""Floresta Estadual Edmundo Navarro de Andrade"", in May 2004, located in the city of Rio Claro, Sao Paulo State, Brazil. Similar anomaly was once described but this is the first specimen found with a bilateral alteration. It may cause confusion in taxonomic identification and even lead to description of new species, increasing the number of synonymies.

Survival, Population Size, and Gonotrophic Cycle Duration of Nyssomyia neivai (Diptera: Psychodidae) at an Endemic Area of American Cutaneous Leishmaniasis in Southeastern Brazil

The survival, absolute population size, gonotrophic cycle duration, and temporal and spatial abundance of Nyssomyia neivai (Pinto) were studied in a rural area endemic for American cutaneous leishmaniasis (ACL) in Conchal, Sao Paulo State, southeastern Brazil, using mark-release-recapture techniques and by monitoring population fluctuation. The monthly abundance exhibited a unimodal pattern, with forest and domicile habitats having the highest relative abundances. A total of 1,873 males and 3,557 females were marked and released during the six experiments, of which 4.1-13.0% of males and 4.1-11.8% of females were recaptured. Daily survivorship estimated from the decline in recaptures per day was 0.681 for males and 0.667 for females. Gonotrophic cycle duration was estimated to be 4.0 d. Absolute population size was calculated using the Lincoln Index and ranged from 861 to 4,612 males and from 2,187 to 19,739 females. The low proportion of females that reach the age when they are potentially infective suggests that N. neivai has a low biological capacity to serve as a vector and that factors such as high biting rates and opportunistic feeding behavior would be needed to enable Leishmania (Viannia) braziliensis Vianna transmission. This agreed with the epidemiological pattern of ACL in southeastern Brazil that is characterized by low incidence, with isolated cases acquired principally within domiciliary habitats.

Observations on the feeding habits of Lutzomyia longipalpis (Lutz & Neiva, 1912) (Diptera: Psychodidae: Phlebotominae) in Campo Grande, an endemic area of visceral leishmaniasis in Mato Grosso do Sul, Brazil

Sand flies were captured weekly with CDC light traps from December 2003 to November 2005 in three areas of Campo Grande, in the Brazilian state of Mato Grosso do Sul. These areas incorporated two patches of remnant forest and five houses. The blood meals of engorged female sand flies were identified using the avidin-biotin system of immunoenzymatic ELISA capture. Most (327/355) of the females analysed were Lutzomyia longipalpis, of which 66.4% reacted with human blood, 64.8% with that of birds and 8.9% with that of dogs. Females that had taken human blood predominated in the residential areas and two forest patches. The following combinations of blood were also detected for L. longipalpis in some of the samples analysed: bird + human (43.4%), bird + human + dog (6.1%). The combination bird + human + dog + pig was also found for Nyssomyia whitmani. Dogs and pigs appear to have little attractiveness for L. longipalpis. The results obtained demonstrate the eclecticism and high anthropophily of L. longipalpis and raise new questions with regard to the importance of dogs in VL epidemiology and the possible role of man as a source of infection for sand flies. (C) 2008 Elsevier B.V. All rights reserved.

Seasonal variation of Lutzomyia longipalpis (Lutz & neiva, 1912) (Diptera : Psychodidae : Phlebotominae) in endemic area of visceral leishmaniasis, Campo Grande, state of Mato Grosso do Sul, Brazil

The seasonal distribution of Lutzomyia longipalpis was studied in two forested and five domiciliary areas of the urban area of Campo Grande; MS, from December 2003 to November 2005. Weekly captures were carried out with CDC light traps positioned on ground and in the canopy inside a residual forest and on the edge (ground) of a woodland and in at least one of the following ecotopes in peridomiciles-a cultivated area, a chicken coop, a pigsty, a kennel, a goat and sheep shelter and an intradomicile. A total of 9519 sand flies were collected, 2666 during the first year and 6853 during the second. L. longipalpis was found throughout the 2-year period, presenting smaller peaks at intervals of 2-3 months and two greater peaks, the first in February and the second in April 2005, soon after periods of heavy rain. Only In one of the woodlands was a significant negative correlation (p < 0.05) between the number of insects and temperature during the first year and the climatic factors (temperature, RHA and rain) was observed. In the domiciliary areas in four domiciles some positive correlations (p < 0.05) occurred in relation to one or more climatic factors; however, the species shows a clear tendency to greater frequency (72%) in the rainy season than in the dry (28%). Thus, we recommend an intensification of the VL control measures applied in Campo Grande, MS, during the rainy season with a view to reducing the risk of the transmission of the disease. (C) 2007 Elsevier B.V. All fights reserved.

Abundance and seasonality of Diptera (Insecta) in a poultry house in the northeast region of the state of Sao Paulo, Brazil

Modem production systems accommodate broody hens in high densities, leading to the accumulation of excrement under the cages. This substrate is excellent for the development of sinantropic flies. Thus, the accomplishment of surveys in these places becomes essential, in order to plan better strategies of control. The present work aimed at studying the entornofauna and the seasonality of the species of dipterous present in the Crisdan poultry house located in the Municipality of Sao Joao da Boa Vista, the State of Sao Paulo, Brazil. In the period of January of 2001 to December of 2002, 1,012,595 flies were captured using the ""jug-trap"". The species were identified: Drosophi-la repleta (Wollaston, 1858), Musca domestica (Linnaeus, 1758), Ophyra spp., Hennetria illucens (Linnaeus, 1758), Fannia canicularis (Linnaeus, 1761), Chrysomya megacephala (Fabricius, 1794), and Sepsidae. More frequently D. repleta and M. domestica had added 99.47% of the dipterous. Increased rainfall and the collection months influenced the sampling of dipterous (P < 0.05). Drosophila repleta was the most abundant species, representing 91% of all captured flies. However, this diptera did not develop at the surveyed site since immatures were not captured therein.

Spider-fed bromeliads: seasonal and interspecific variation in plant performance

Background and Aims Several animals that live on bromeliads can contribute to plant nutrition through nitrogen provisioning (digestive mutualism). The bromeliad-living spider Psecas chapoda (Salticidae) inhabits and breeds on Bromelia balansae in regions of South America, but in specific regions can also appear on Ananas comosus (pineapple) plantations and Aechmea distichantha. Methods Using isotopic and physiological methods in greenhouse experiments, the role of labelled ((15)N) spider faeces and Drosophila melanogaster flies in the nutrition and growth of each host plant was evaluated, as well as seasonal variation in the importance of this digestive mutualism. Key Results Spiders contributed 0.6 +/- 0.2% (mean +/- s.e.; dry season) to 2.7 +/- 1% (wet season) to the total nitrogen in B. balansae, 2.4 +/- 0.4% (dry) to 4.1 +/- 0.3% (wet) in An. comosus and 3.8 +/- 0.4% (dry) to 5 +/- 1% (wet) in Ae. distichantha. In contrast, flies did not contribute to the nutrition of these bromeliads. Chlorophylls and carotenoid concentrations did not differ among treatments. Plants that received faeces had higher soluble protein concentrations and leaf growth (RGR) only during the wet season. Conclusions These results indicate that the mutualism between spiders and bromeliads is seasonally restricted, generating a conditional outcome. There was interspecific variation in nutrient uptake, probably related to each species` performance and photosynthetic pathways. Whereas B. balansae seems to use nitrogen for growth, Ae. distichantha apparently stores nitrogen for stressful nutritional conditions. Bromeliads absorbed more nitrogen coming from spider faeces than from flies, reinforcing the beneficial role played by predators in these digestive mutualisms.

On the Relationship between Zaprionus spinipilus Chassagnard & McEvey and Z. vittiger Coquillett, the Type Species of the Genus Zaprionus (Diptera: Drosophilidae)

Zaprionus vittiger Coquillett is the type species of the genus Zaprionus Coquillett. However, the species is only known from five old museum specimens collected from South Africa and Malawi. It has often been confused with many other Zaprionus species, especially with Z. spinipilus Chassagnard & McEvey, a widespread species in Africa known from Madagascar, Malawi, Ethiopia and Cameroon. We have recently collected flies from the type localities of both species (South Africa and Madagascar, respectively). This has prompted us to test the taxonomic boundaries of these two nominal species using molecular (the mitochondrial COII and the nuclear Amyrel genes), chromosomal, morphological (internal and external genitalia), and reproductive isolation analyses. The results suggest Z. spinipilus to be a junior synonym to Z. vittiger.

Sequence and function of lysosomal and digestive cathepsin D-like proteinases of Musca domestica midgut

Musca domestica larvae display in anterior and middle midgut contents, a proteolytic activity with pH optimum of 3.0-3.5 and kinetic properties like cathepsin D. Three cDNAs coding for preprocathepsin D-like proteinases (ppCAD 1, ppCAD 2, ppCAD 3) were cloned from a M. domestica midgut cDNA library. The coded protein sequences included the signal peptide, propeptide and mature enzyme that has all conserved catalytic and substrate binding residues found in bovine lysosomal cathepsin D. Nevertheless, ppCAD 2 and ppCAD 3 lack the characteristic proline loop and glycosylation sites. A comparison among the sequences of cathepsin D-like enzymes from some vertebrates and those found in M. domestica and in the genomes of Aedes aegypti, Drosophila melanogaster, Tribolium castaneum, and Bombyx mori showed that only flies have enzymes lacking the proline loop (as defined by the motif: DxPxPx(G/A)P), thus resembling vertebrate pepsin. ppCAD 3 should correspond to the digestive cathepsin D-like proteinase (CAD) found in enzyme assays because: (1) it seems to be the most expressed CAD, based on the frequency of ESTs found. (2) The mRNA for CAD 3 is expressed only in the anterior and proximal middle midgut. (3) Recombinant procathepsin D-like proteinase (pCAD 3), after auto-activation has a pH optimum of 2.5-3.0 that is close to the luminal pH of M. domestica midgut. (4) Immunoblots of proteins from different tissues revealed with anti-pCAD 3 serum were positive only in samples of anterior and middle midgut tissue and contents. (5) CAD 3 is localized with immunogold inside secretory vesicles and around microvilli in anterior and middle midguit cells. The data support the view that on adapting to deal with a bacteria-rich food in an acid midgut region, M. domestica digestive CAD resulted from the same archetypical gene as the intracellular cathepsin D, paralleling what happened with vertebrates. The lack of the proline loop may be somehow associated with the extracellular role of both pepsin and digestive CAD 3. (C) 2009 Elsevier Ltd. All rights reserved.

The DNA puff 4C expresses a salivary secretion protein in Trichosia pubescens (Diptera; Sciaridae)

DNA puffs are genomic regions of polytene chromosomes that undergo developmentally controlled DNA amplification and transcription in salivary glands of sciarid flies. Here, we tested the hypothesis that DNA puff genes code for salivary proteins in Trichosia pubescens. To do that, we generated antibodies against saliva and immunoscreened a cDNA library made from salivary glands. We isolated clones corresponding to DNA puff regions, including clone D-50 that contained the entire coding sequence of the previously isolated C4B1 gene from puff 4C. Indeed, we showed that puff 4C is a DNA puff region detecting its local transcription and its extra rounds of DNA incorporation compared to neighboring regions. We further confirmed D-50 clone identity in Western blots reacted with the anti-saliva anitiserum. We detected a recombinant protein expressed by this clone that had the expected size for a full-length product of the gene. We end with a discussion of the relationship between DNA puff genes and their products.