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.

Enzymatically-triggered, isothermally responsive polymers: re-programming poly(oligoethylene glycols) to respond to phosphatase

Polymers which can respond to externally applied stimuli have found much application in the biomedical field due to their (reversible) coil–globule transitions. Polymers displaying a lower critical solution temperature are the most commonly used, but for blood-borne (i.e., soluble) biomedical applications the application of heat is not always possible, nor practical. Here we report the design and synthesis of poly(oligoethylene glycol methacrylate)-based polymers whose cloud points are easily varied by alkaline phosphatase-mediated dephosphorylation. By fine-tuning the density of phosphate groups on the backbone, it was possible to induce an isothermal transition: A change in solubility triggered by removal of a small number of phosphate esters from the side chains activating the LCST-type response. As there was no temperature change involved, this serves as a model of a cell-instructed polymer response. Finally, it was found that both polymers were non cytotoxic against MCF-7 cells (at 1 mg·mL–1), which confirms promise for biomedical applications.

Shape-controlled continuous synthesis of metal nanostructures

A segmented flow-based microreactor is used for the continuous production of faceted nanocrystals. Flow segmentation is proposed as a versatile tool to manipulate the reduction kinetics and control the growth of faceted nanostructures; tuning the size and shape. Switching the gas from oxygen to carbon monoxide permits the adjustment in nanostructure growth from 1D (nanorods) to 2D (nanosheets). CO is a key factor in the formation of Pd nanosheets and Pt nanocubes; operating as a second phase, a reductant, and a capping agent. This combination confines the growth to specific structures. In addition, the segmented flow microfluidic reactor inherently has the ability to operate in a reproducible manner at elevated temperatures and pressures whilst confining potentially toxic reactants, such as CO, in nanoliter slugs. This continuous system successfully synthesised Pd nanorods with an aspect ratio of 6; thin palladium nanosheets with a thickness of 1.5 nm; and Pt nanocubes with a 5.6 nm edge length, all in a synthesis time as low as 150 s.

An accelerating cosmology without dark energy

The negative pressure accompanying gravitationally-induced particle creation can lead to a cold dark matter (CDM) dominated, accelerating Universe (Lima et al. 1996 [1]) without requiring the presence of dark energy or a cosmological constant. In a recent study, Lima et al. 2008 [2] (LSS) demonstrated that particle creation driven cosmological models are capable of accounting for the SNIa observations [3] of the recent transition from a decelerating to an accelerating Universe, without the need for Dark Energy. Here we consider a class of such models where the particle creation rate is assumed to be of the form Gamma = beta H + gamma H(0), where H is the Hubble parameter and H(0) is its present value. The evolution of such models is tested at low redshift by the latest SNe Ia data provided by the Union compilation [4] and at high redshift using the value of z(eq), the redshift of the epoch of matter - radiation equality, inferred from the WMAP constraints on the early Integrated Sachs-Wolfe (ISW) effect [5]. Since the contributions of baryons and radiation were ignored in the work of LSS, we include them in our study of this class of models. The parameters of these more realistic models with continuous creation of CDM are constrained at widely-separated epochs (z(eq) approximate to 3000 and z approximate to 0) in the evolution of the Universe. The comparison of the parameter values, {beta, gamma}, determined at these different epochs reveals a tension between the values favored by the high redshift CMB constraint on z(eq) from the ISW and those which follow from the low redshift SNIa data, posing a potential challenge to this class of models. While for beta = 0 this conflict is only at less than or similar to 2 sigma, it worsens as beta increases from zero.

Note on ""Continuous matter creation and the acceleration of the universe: the growth of density fluctuations""

Recently, de Roany and Pacheco (Gen Relativ Gravit, doi:10.1007/s10714-010-1069-2) performed a Newtonian analysis on the evolution of perturbations for a class of relativistic cosmological models with Creation of Cold Dark Matter (CCDM) proposed by the present authors (Lima et al. in JCAP 1011:027, 2010). In this note we demonstrate that the basic equations adopted in their work do not recover the specific (unperturbed) CCDM model. Unlike to what happens in the original CCDM cosmology, their basic conclusions refer to a decelerating cosmological model in which there is no transition from a decelerating to an accelerating regime as required by SNe type Ia and complementary observations.

Skull Modularity in Neotropical Marsupials and Monkeys: Size Variation and Evolutionary Constraint and Flexibility

An organism is built through a series of contingent factors, yet it is determined by historical, physical, and developmental constraints. A constraint should not be understood as an absolute obstacle to evolution, as it may also generate new possibilities for evolutionary change. Modularity is, in this context, an important way of organizing biological information and has been recognized as a central concept in evolutionary biology bridging on developmental, genetics, morphological, biochemical, and physiological studies. In this article, we explore how modularity affects the evolution of a complex system in two mammalian lineages by analyzing correlation, variance/covariance, and residual matrices (without size variation). We use the multivariate response to selection equation to simulate the behavior of Eutheria and Metharia skulls in terms of their evolutionary flexibility and constraints. We relate these results to classical approaches based on morphological integration tests based on functional/developmental hypotheses. Eutherians (Neotropical primates) showed smaller magnitudes of integration compared with Metatheria (didelphids) and also skull modules more clearly delimited. Didelphids showed higher magnitudes of integration and their modularity is strongly influenced by within-groups size variation to a degree that evolutionary responses are basically aligned with size variation. Primates still have a good portion of the total variation based on size; however, their enhanced modularization allows a broader spectrum of responses, more similar to the selection gradients applied (enhanced flexibility). Without size variation, both groups become much more similar in terms of modularity patterns and magnitudes and, consequently, in their evolutionary flexibility. J. Exp. Zool. (Mol. Dev. Evol.) 314B:663-683, 2010. (C) 2010 Wiley-Liss, Inc.

Cladistic analysis of continuous modularized traits provides phylogenetic signals in Homo evolution

Evolutionary novelties in the skeleton are usually expressed as changes in the timing of growth of features intrinsically integrated at different hierarchical levels of development(1). As a consequence, most of the shape- traits observed across species do vary quantitatively rather than qualitatively(2), in a multivariate space(3) and in a modularized way(4,5). Because most phylogenetic analyses normally use discrete, hypothetically independent characters(6), previous attempts have disregarded the phylogenetic signals potentially enclosed in the shape of morphological structures. When analysing low taxonomic levels, where most variation is quantitative in nature, solving basic requirements like the choice of characters and the capacity of using continuous, integrated traits is of crucial importance in recovering wider phylogenetic information. This is particularly relevant when analysing extinct lineages, where available data are limited to fossilized structures. Here we show that when continuous, multivariant and modularized characters are treated as such, cladistic analysis successfully solves relationships among main Homo taxa. Our attempt is based on a combination of cladistics, evolutionary- development- derived selection of characters, and geometric morphometrics methods. In contrast with previous cladistic analyses of hominid phylogeny, our method accounts for the quantitative nature of the traits, and respects their morphological integration patterns. Because complex phenotypes are observable across different taxonomic groups and are potentially informative about phylogenetic relationships, future analyses should point strongly to the incorporation of these types of trait.

Continuous and High-Level In Vivo Delivery of Endostatin From Recombinant Cells Encapsulated in TheraCyte (R) Immunoisolation Devices

Endostatin (ES) is a potent inhibitor of angiogenesis and tumor growth. Continuous ES delivery of ES improves the efficacy and potency of the antitumoral therapy. The TheraCyte (R) system is a polytetrafluoroethylene (PTFE) semipermeable membrane macroencapsulation system for implantation of genetically engineered cells specially designed for the in vivo delivery of therapeutic proteins, such as ES, which circumvents the problem of limited half-life and variation in circulating levels. In order to enable neovascularization at the tissues adjacent to the devices prior to ES secretion by the cells inside them, we designed a scheme in which empty TheraCyte (R) devices were preimplanted SC into immunodeficient mice. Only after healing (17 days later) were Chinese hamster ovary cells expressing ES injected into the preimplanted devices. In another model for device implantation, the cells expressing ES where loaded into the immunoisolation devices prior to implantation into the animals, and the TheraCyte (R) were then immediately implanted SC into the mice. Throughout the 2-month study, constant high ES levels of up to 3.7 mu g/ml were detected in the plasma of the mice preimplanted with the devices, while lower but also constant levels of ES (up to 2.1 mu g/ml plasma) were detected in the mice that had received devices preloaded with the ES-expressing cells. Immunohistochemistry using anti-ES antibody showed reaction within the device and outside it, demonstrating that ES, secreted by the confined recombinant cells, permeated through the membrane and reached the surrounding tissues.

Perinatal salt restriction: A new pathway to programming adiposity indices in adult female Wistar rats

Low birth weight has been associated with increased obesity in adulthood. It has been shown that dietary salt restriction during intrauterine life induces low birth weight and insulin resistance in adult Wistar rats. The present study had a two-fold objective: to evaluate the effects that low salt intake during pregnancy and lactation has on the amount and distribution of adipose tissue; and to determine whether the phenotypic changes in fat mass in this model are associated with alterations in the activity of the renin-angiotensin system. Maternal salt restriction was found to reduce birth weight in male and female offspring. In adulthood, the female offspring of dams fed the low-salt diet presented higher adiposity indices than those seen in the offspring of dams fed a normal-salt diet. This was attributed to the fact that adipose tissue mass (retroperitoneal but not gonadal, mesenteric or inguinal) was greater in those rats than in the offspring of dams fed a normal diet. The adult offspring of dams fed the low-salt diet, compared to those dams fed a normal-salt diet, presented the following: plasma leptin levels higher in males and lower in females; plasma renin activity higher in males but not in females; and no differences in body weight, mean arterial blood pressure or serum angiotensin-converting enzyme activity. Therefore, low salt intake during pregnancy might lead to the programming of obesity in adult female offspring. (c) 2009 Elsevier Inc. All rights reserved.

THE POINCARE-BENDIXSON THEOREM ON THE KLEIN BOTTLE FOR CONTINUOUS VECTOR FIELDS

We present a version of the Poincare-Bendixson Theorem on the Klein bottle K(2) for continuous vector fields. As a consequence, we obtain the fact that K(2) does not admit continuous vector fields having a omega-recurrent injective trajectory.

THE IMPLICIT FUNCTION THEOREM FOR CONTINUOUS FUNCTIONS

In the present paper we obtain a new homological version of the implicit function theorem and some versions of the Darboux theorem. Such results are proved for continuous maps on topological manifolds. As a consequence. some versions of these classic theorems are proved when we consider differenciable (not necessarily C-1) maps.

An integer linear programming approach for approximate string comparison

We introduce a problem called maximum common characters in blocks (MCCB), which arises in applications of approximate string comparison, particularly in the unification of possibly erroneous textual data coming from different sources. We show that this problem is NP-complete, but can nevertheless be solved satisfactorily using integer linear programming for instances of practical interest. Two integer linear formulations are proposed and compared in terms of their linear relaxations. We also compare the results of the approximate matching with other known measures such as the Levenshtein (edit) distance. (C) 2008 Elsevier B.V. All rights reserved.

Aerosol inorganic composition at a tropical site: Discrepancies between filter-based sampling and a semi-continuous method

The concentrations of the water-soluble inorganic aerosol species, ammonium (NH4+), nitrate (NO3-), chloride (Cl-), and sulfate (SO42-), were measured from September to November 2002 at a pasture site in the Amazon Basin (Rondnia, Brazil) (LBA-SMOCC). Measurements were conducted using a semi-continuous technique (Wet-annular denuder/Steam-Jet Aerosol Collector: WAD/SJAC) and three integrating filter-based methods, namely (1) a denuder-filter pack (DFP: Teflon and impregnated Whatman filters), (2) a stacked-filter unit (SFU: polycarbonate filters), and (3) a High Volume dichotomous sampler (HiVol: quartz fiber filters). Measurements covered the late dry season (biomass burning), a transition period, and the onset of the wet season (clean conditions). Analyses of the particles collected on filters were performed using ion chromatography (IC) and Particle-Induced X-ray Emission spectrometry (PIXE). Season-dependent discrepancies were observed between the WAD/SJAC system and the filter-based samplers. During the dry season, when PM2.5 (D-p <= 2.5 mu m) concentrations were similar to 100 mu g m(-3), aerosol NH4+ and SO42- measured by the filter-based samplers were on average two times higher than those determined by the WAD/SJAC. Concentrations of aerosol NO3- and Cl- measured with the HiVol during daytime, and with the DFP during day- and nighttime also exceeded those of the WAD/SJAC by a factor of two. In contrast, aerosol NO3- and Cl- measured with the SFU during the dry season were nearly two times lower than those measured by the WAD/SJAC. These differences declined markedly during the transition period and towards the cleaner conditions during the onset of the wet season (PM2.5 similar to 5 mu g m(-3)); when filter-based samplers measured on average 40-90% less than the WAD/SJAC. The differences were not due to consistent systematic biases of the analytical techniques, but were apparently a result of prevailing environmental conditions and different sampling procedures. For the transition period and wet season, the significance of our results is reduced by a low number of data points. We argue that the observed differences are mainly attributable to (a) positive and negative filter sampling artifacts, (b) presence of organic compounds and organosulfates on filter substrates, and (c) a SJAC sampling efficiency of less than 100%.

New mechanical samples positioning system for irradiations on a radial channel at nuclear research reactor in a full-power continuous operation

This paper describes a new mechanical samples positioning system that allows the safe placement and removal of biological samples for prolonged irradiation, in a nuclear reactor during full-power continuous operation. Also presented herein the materials of construction and operating principles. Additionally, this sample positioning system is compared with an existing pneumatic and automated transfer system, already available at the research reactors. The system consists of a mechanical arm with a claw, which can deliver the samples for irradiations without reactor shutdown. It was installed in the lEA-R1 research reactor at Instituto de Pesquisas Energeticas e Nucleares (IPEN), Sao Paulo, Brazil, and for the past 5 years, the system has successfully operated and allowed the conducting of important experiments. As a result of its introduction, the facility has been in a position to positively respond to the increased demand in studies of biology, medicine, physics, engineering, detector/dosimeter calibrations, etc. It is one example of the appropriated technologies that save energy and resources. (C) 2010 Elsevier Ltd. All rights reserved.