Advances in sequential data assimilation and numerical weather forecasting: an Ensemble Transform Kalman-Bucy Filter, a study on clustering in deterministic Ensemble Square Root Filters, and a test of a new time stepping scheme in an atmospheric model

This dissertation deals with aspects of sequential data assimilation (in particular ensemble Kalman filtering) and numerical weather forecasting. In the first part, the recently formulated Ensemble Kalman-Bucy (EnKBF) filter is revisited. It is shown that the previously used numerical integration scheme fails when the magnitude of the background error covariance grows beyond that of the observational error covariance in the forecast window. Therefore, we present a suitable integration scheme that handles the stiffening of the differential equations involved and doesn’t represent further computational expense. Moreover, a transform-based alternative to the EnKBF is developed: under this scheme, the operations are performed in the ensemble space instead of in the state space. Advantages of this formulation are explained. For the first time, the EnKBF is implemented in an atmospheric model. The second part of this work deals with ensemble clustering, a phenomenon that arises when performing data assimilation using of deterministic ensemble square root filters in highly nonlinear forecast models. Namely, an M-member ensemble detaches into an outlier and a cluster of M-1 members. Previous works may suggest that this issue represents a failure of EnSRFs; this work dispels that notion. It is shown that ensemble clustering can be reverted also due to nonlinear processes, in particular the alternation between nonlinear expansion and compression of the ensemble for different regions of the attractor. Some EnSRFs that use random rotations have been developed to overcome this issue; these formulations are analyzed and their advantages and disadvantages with respect to common EnSRFs are discussed. The third and last part contains the implementation of the Robert-Asselin-Williams (RAW) filter in an atmospheric model. The RAW filter is an improvement to the widely popular Robert-Asselin filter that successfully suppresses spurious computational waves while avoiding any distortion in the mean value of the function. Using statistical significance tests both at the local and field level, it is shown that the climatology of the SPEEDY model is not modified by the changed time stepping scheme; hence, no retuning of the parameterizations is required. It is found the accuracy of the medium-term forecasts is increased by using the RAW filter.

Effect of a methyl group on the spontaneous resolution of a square-pyramidal coordination compound: crystal packing and conglomerate formation

Reaction of [Cu(pic)2]·2H2O (where pic stands for 2-picolinato) with 2-({[2-(dimethylamino)ethyl]amino}methyl)phenol (HL1) produces the square-pyramidal complex [CuL1(pic)] (1), which crystallizes as a conglomerate (namely a mixture of optically pure crystals) in the Sohncke space group P212121. The use of the methylated ligand at the benzylic position, i.e. (±)-2-(1-{[2-(dimethylamino)ethyl]amino}ethyl)phenol (HL2), yields the analogous five-coordinate complex [CuL2(pic)] (2) that crystallizes as a true racemate (namely the crystals contain both enantiomers) in the centrosymmetric space group P21/c. Density functional theory (DFT) calculations indicate that the presence of the methyl group indeed leads to a distinct crystallization behaviour, not only by intramolecular steric effects, but also because its involvement in non-covalent C–H···π and hydrophobic intermolecular contacts appears to be an important factor contributing to the crystal-lattice (stabilizing) energy of 2

Sparse probability density function estimation using the minimum integrated square error

We develop a new sparse kernel density estimator using a forward constrained regression framework, within which the nonnegative and summing-to-unity constraints of the mixing weights can easily be satisfied. Our main contribution is to derive a recursive algorithm to select significant kernels one at time based on the minimum integrated square error (MISE) criterion for both the selection of kernels and the estimation of mixing weights. The proposed approach is simple to implement and the associated computational cost is very low. Specifically, the complexity of our algorithm is in the order of the number of training data N, which is much lower than the order of N2 offered by the best existing sparse kernel density estimators. Numerical examples are employed to demonstrate that the proposed approach is effective in constructing sparse kernel density estimators with comparable accuracy to those of the classical Parzen window estimate and other existing sparse kernel density estimators.

Chemistry and structure by design: ordered CuNi(CN)4 sheets with copper(II) in a square-planar environment

Layered copper–nickel cyanide, CuNi(CN)4, a 2-D negative thermal expansion material, is one of a series of copper(II)-containing cyanides derived from Ni(CN)2. In CuNi(CN)4, unlike in Ni(CN)2, the cyanide groups are ordered generating square-planar Ni(CN)4 and Cu(NC)4 units. The adoption of square-planar geometry by Cu(II) in an extended solid is very unusual.

Least square projection: A fast high-precision multidimensional projection technique and its application to document mapping

The problem of projecting multidimensional data into lower dimensions has been pursued by many researchers due to its potential application to data analyses of various kinds. This paper presents a novel multidimensional projection technique based on least square approximations. The approximations compute the coordinates of a set of projected points based on the coordinates of a reduced number of control points with defined geometry. We name the technique Least Square Projections ( LSP). From an initial projection of the control points, LSP defines the positioning of their neighboring points through a numerical solution that aims at preserving a similarity relationship between the points given by a metric in mD. In order to perform the projection, a small number of distance calculations are necessary, and no repositioning of the points is required to obtain a final solution with satisfactory precision. The results show the capability of the technique to form groups of points by degree of similarity in 2D. We illustrate that capability through its application to mapping collections of textual documents from varied sources, a strategic yet difficult application. LSP is faster and more accurate than other existing high-quality methods, particularly where it was mostly tested, that is, for mapping text sets.

A group topology on the free abelian group of cardinality c that makes its square countably compact

Under p = c, we prove that it is possible to endow the free abelian group of cardinality c with a group topology that makes its square countably compact. This answers a question posed by Madariaga-Garcia and Tomita and by Tkachenko. We also prove that there exists a Wallace semigroup (i.e., a countably compact both-sided cancellative topological semigroup which is not a topological group) whose square is countably compact. This answers a question posed by Grant.

Development of a sequential injection-square wave voltammetry method for determination of paraquat in water samples employing the hanging mercury drop electrode

This work describes the development and optimization of a sequential injection method to automate the determination of paraquat by square-wave voltammetry employing a hanging mercury drop electrode. Automation by sequential injection enhanced the sampling throughput, improving the sensitivity and precision of the measurements as a consequence of the highly reproducible and efficient conditions of mass transport of the analyte toward the electrode surface. For instance, 212 analyses can be made per hour if the sample/standard solution is prepared off-line and the sequential injection system is used just to inject the solution towards the flow cell. In-line sample conditioning reduces the sampling frequency to 44 h(-1). Experiments were performed in 0.10 M NaCl, which was the carrier solution, using a frequency of 200 Hz, a pulse height of 25 mV, a potential step of 2 mV, and a flow rate of 100 mu L s(-1). For a concentration range between 0.010 and 0.25 mg L(-1), the current (i(p), mu A) read at the potential corresponding to the peak maximum fitted the following linear equation with the paraquat concentration (mg L(-1)): ip = (-20.5 +/- 0.3) Cparaquat -(0.02 +/- 0.03). The limits of detection and quantification were 2.0 and 7.0 mu g L(-1), respectively. The accuracy of the method was evaluated by recovery studies using spiked water samples that were also analyzed by molecular absorption spectrophotometry after reduction of paraquat with sodium dithionite in an alkaline medium. No evidence of statistically significant differences between the two methods was observed at the 95% confidence level.

Square wave adsorptive cathodic stripping voltarnmetry automated by sequential injection analysis - Potentialities and limitations exemplified by the determination of methyl parathion in water samples

This paper describes the development and evaluation of a sequential injection method to automate the determination of methyl parathion by square wave adsorptive cathodic stripping voltammetry exploiting the concept of monosegmented flow analysis to perform in-line sample conditioning and standard addition. Accumulation and stripping steps are made in the sample medium conditioned with 40 mmol L-1 Britton-Robinson buffer (pH 10) in 0.25 mol L-1 NaNO3. The homogenized mixture is injected at a flow rate of 10 mu Ls(-1) toward the flow cell, which is adapted to the capillary of a hanging drop mercury electrode. After a suitable deposition time, the flow is stopped and the potential is scanned from -0.3 to -1.0 V versus Ag/AgCl at frequency of 250 Hz and pulse height of 25 mV The linear dynamic range is observed for methyl parathion concentrations between 0.010 and 0.50 mgL(-1), with detection and quantification limits of 2 and 7 mu gL(-1), respectively. The sampling throughput is 25 h(-1) if the in line standard addition and sample conditioning protocols are followed, but this frequency can be increased up to 61 h(-1) if the sample is conditioned off-line and quantified using an external calibration curve. The method was applied for determination of methyl parathion in spiked water samples and the accuracy was evaluated either by comparison to high performance liquid chromatography with UV detection, or by the recovery percentages. Although no evidences of statistically significant differences were observed between the expected and obtained concentrations, because of the susceptibility of the method to interference by other pesticides (e.g., parathion, dichlorvos) and natural organic matter (e.g., fulvic and humic acids), isolation of the analyte may be required when more complex sample matrices are encountered. (C) 2007 Elsevier B.V. All rights reserved.

Electroanalytical determination of promethazine hydrochloride in pharmaceutical formulations on highly boron-doped diamond electrodes using square-wave adsorptive voltammetry

The electrochemical oxidation of promethazine hydrochloride was made on highly boron-doped diamond electrodes. Cyclic voltammetry experiments showed that the oxidation mechanisms involved the formation of an adsorbed product that is more readily oxidized, producing a new peak with lower potential values whose intensity can be increased by applying the accumulation potential for given times. The parameters were optimized and the highest current intensities were obtained by applying +0.78 V for 30 seconds. The square-wave adsorptive voltammetry results obtained in BR buffer showed two well-defined peaks, dependent on the pH and on the voltammetric parameters. The best responses were obtained at pH 4.0, frequency of 50 s(-1), step of 2 mV, and amplitude of 50 mV. Under these conditions, linear responses were obtained for concentrations from 5.96 x 10(-7) to 4.76 x 10(-6) mol L-1, and calculated detection limits of 2.66 x 10(-8) mol L-1 (8.51 mu g L-1) for peak 1 and of 4.61 x 10(-8) mol L-1 (14.77 mu g L-1) for peak 2. The precision and accuracy were evaluated by repeatability and reproducibility experiments, which yielded values of less than 5.00% for both voltammetric peaks. ne applicability of this procedure was tested on commercial formulations of promethazine hydrochloride by observing the stability, specificity, recovery and precision of the procedure in complex samples. All results obtained were compared to recommended procedure by British Pharmacopeia. The voltammetric results indicate that the proposed procedure is stable and sensitive, with good reproducibility even when the accumulation steps involve short times. It is therefore very suitable for the development of the electroanalytical procedure, providing adequate sensitivity and a reliable method.

A simple demonstration of the existence of the jordan canonical form for any square matrix

All the demonstrations known to this author of the existence of the Jordan Canonical Form are somewhat complex - usually invoking the use of new spaces, and what not. These demonstrations are usually too difficult for an average Mathematics student to understand how he or she can obtain the Jordan Canonical Form for any square matrix. The method here proposed not only demonstrates the existence of such forms but, additionally, shows how to find them in a step by step manner. I do not claim that the following demonstration is in any way “elegant” (by the standards of elegance in fashion nowadays among mathematicians) but merely simple (undergraduate students taking a fist course in Matrix Algebra would understand how it works).

Minimum size for the occurrence of vortex matter in a square mesoscopic superconductor

The study of superconducting samples in mesoscopic scale presented a remarkable improvement during the last years. Certainly, such interest is based on the fact that when the size of the samples is close to the order of the temperature dependent coherence length xi(T), and/or the size of the penetration depth lambda(T), there are some significant modifications on the physical properties of the superconducting state. This contribution tests the square cross-section size limit for the occurrence (or not) of vortices in mesoscopic samples of area L-2, where L varies discretely from 1 xi(0) to 8 xi(0).The time dependent Ginzburg-Landau (TDGL) equations approach is used upon taking the order parameter and the local magnetic field invariant along the z-direction. The vortex configurations at the equilibrium can be obtained from the TDGL equations for superconductivity as the system relaxes to the stationary state.The obtained results show that the limit of vortex penetration is for the square sample of size 3 xi(0) x 3 xi(0) in which only a single vortex are allowed into the sample. For smaller specimens, no vortex can be formed and the field entrance into the sample is continuous and the total flux penetration occurs at higher values of H/H-c2(0), where H-c2(T) is the upper critical field. Otherwise, for larger samples different vortices patterns can be observed depending on the sample size. (c) 2007 Elsevier B.V. All rights reserved.

Vortex dynamics in mesoscopic superconducting square of variable surface

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

Vortex-antivortex annihilation dynamics in a square mesoscopic superconducting cylinder

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