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.

Modelling the time-evolution of phytoplankton size spectra from satellite remote sensing

A dynamic size-structured model is developed for phytoplankton and nutrients in the oceanic mixed layer and applied to extract phytoplankton biomass at discrete size fractions from remotely sensed, ocean-colour data. General relationships between cell size and biophysical processes (such as sinking, grazing, and primary production) of phytoplankton were included in the model through a bottom–up approach. Time-dependent, mixed-layer depth was used as a forcing variable, and a sequential data-assimilation scheme was implemented to derive model trajectories. From a given time-series, the method produces estimates of size-structured biomass at every observation, so estimates seasonal succession of individual phytoplankton size, derived here from remote sensing for the first time. From these estimates, normalized phytoplankton biomass size spectra over a period of 9 years were calculated for one location in the North Atlantic. Further analysis demonstrated that strong relationships exist between the seasonal trends of the estimated size spectra and the mixed-layer depth, nutrient biomass, and total chlorophyll. The results contain useful information on the time-dependent biomass flux in the pelagic ecosystem.

CP-OQAM-OFDM based SC-FDMA: adjustable user bandwidth and space-time coding

The discrete Fourier transmission spread OFDM DFTS-OFDM) based single-carrier frequency division multiple access (SC-FDMA) has been widely adopted due to its lower peak-to-average power ratio (PAPR) of transmit signals compared with OFDM. However, the offset modulation, which has lower PAPR than general modulation, cannot be directly applied into the existing SC-FDMA. When pulse-shaping filters are employed to further reduce the envelope fluctuation of transmit signals of SC-FDMA, the spectral efficiency degrades as well. In order to overcome such limitations of conventional SC-FDMA, this paper for the first time investigated cyclic prefixed OQAMOFDM (CP-OQAM-OFDM) based SC-FDMA transmission with adjustable user bandwidth and space-time coding. Firstly, we propose CP-OQAM-OFDM transmission with unequally-spaced subbands. We then apply it to SC-FDMA transmission and propose a SC-FDMA scheme with the following features: a) the transmit signal of each user is offset modulated single-carrier with frequency-domain pulse-shaping; b) the bandwidth of each user is adjustable; c) the spectral efficiency does not decrease with increasing roll-off factors. To combat both inter-symbolinterference and multiple access interference in frequencyselective fading channels, a joint linear minimum mean square error frequency domain equalization using a prior information with low complexity is developed. Subsequently, we construct space-time codes for the proposed SC-FDMA. Simulation results confirm the powerfulness of the proposed CP-OQAM-OFDM scheme (i.e., effective yet with low complexity).

First contribution of mites (Acari) to the forensic analysis of hanged corpses: a case study from Spain

This case study from North Spain, highlights the importance of the collection of mites in addition to insects, from crime scenes or corpses subjected to environmental constraints that reduce or minimise insect activity, such as hanged corpses. In addition, this analysis highlights the relevance of arthropods’ collection in the field, even after the corpse has been moved away for autopsy. Four species of mites, phoretic on carrion (Silphidae) and rove (Staphylinidae) beetles, complemented and reinforced the autopsy analysis as well as the scarce information provided by insect activity. Poecilochirus carabi Canestrini & Canestrini, 1882 and Poecilochirus (Physoparasitus) davydovae Hyatt, 1980 (Mesostigmata: Parasitidae) were found in association with two Silphidae, Nicrophorus Fabricius, 1775 and Necrodes Leach, 1815, only when sampled in the autopsy room; this is suggestive of host-switching of mites and was likely due to the lack of availability of specific carriers in the field. The interpretation of the activity of Parasitidae mites both in the field and the autopsy room allows a better understanding of the timing and circumstances of decomposition. Phoretic deutonymphs of Pelzneria Scheucher 1957 (Astigmata: Histiostomatidae) were highly abundant, mostly P. crenulata Oudemans, 1909 and are reported for the first time on a Staphylinidae rove beetle, Creophilus maxillosus (L., 1758). Surprisingly, in this case study no Pelzneria were associated with the Silphidae found, which are their most common hosts, such as Necrodes littoralis (L., 1758) and Nicrophorus interruptus (Stephens, 1830). All histiostomatids were removed from the staphylinid (rove beetle) collected from the soil, at the scene of death, suggesting a recent arrival of the beetle. The occurrence of Staphylinidae beetles and their associated mites, such as Parasitidae and Pelzneria, and the information they provided would have been easily overlooked or lost if only the autopsy sampling would have been considered in the analysis of the case. The four mite species are reported for the first time for the Iberian Peninsula.

Identification of domestication-related loci associated with flowering time and seed size in soybean with the RAD-seq genotyping method

Flowering time and seed size are traits related to domestication. However, identification of domestication-related loci/genes of controlling the traits in soybean is rarely reported. In this study, we identified a total of 48 domestication-related loci based on RAD-seq genotyping of a natural population comprising 286 accessions. Among these, four on chromosome 12 and additional two on chromosomes 11 and 15 were associated with flowering time, and four on chromosomes 11 and 16 were associated with seed size. Of the five genes associated with flowering time and the three genes associated with seed size, three genes Glyma11g18720, Glyma11g15480 and Glyma15g35080 were homologous to Arabidopsis genes, additional five genes were found for the first time to be associated with these two traits. Glyma11g18720 and Glyma05g28130 were co-expressed with five genes homologous to flowering time genes in Arabidopsis, and Glyma11g15480 was co-expressed with 24 genes homologous to seed development genes in Arabidopsis. This study indicates that integration of population divergence analysis, genome-wide association study and expression analysis is an efficient approach to identify candidate domestication-related genes.

Second language influence on first language motion event encoding and categorization in Spanish-speaking children learning L2 English

Studies show cross-linguistic differences in motion event encoding, such that English speakers preferentially encode manner of motion more than Spanish speakers, who preferentially encode path of motion. Focusing on native Spanish speaking children (aged 5;00-9;00) learning L2 English, we studied path and manner verb preferences during descriptions of motion stimuli, and tested the linguistic relativity hypothesis by investigating categorization preferences in a non-verbal similarity judgement task of motion clip triads. Results revealed L2 influence on L1 motion event encoding, such that bilinguals used more manner verbs and fewer path verbs in their L1, under the influence of English. We found no effects of linguistic structure on non-verbal similarity judgements, and demonstrate for the first time effects of L2 on L1 lexicalization in child L2 learners in the domain of motion events. This pattern of verbal behaviour supports theories of bilingual semantic representation that postulate a merged lexico-semantic system in early bilinguals.

Solvatochromism in Binary Mixtures: First Report on a Solvation Free Energy Relationship between Solvent Exchange Equilibrium Constants and the Properties of the Medium

We have employed UV-vis spectroscopy in order to investigate details of the solvation of six solvatochromic indicators, hereafter designated as ""probes"", namely, 2,6-diphenyl-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (RB); 4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePM; 1-methylquinolinium-8-olate, QB; 2-bromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr, 2,6-dichloro-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (WB); and 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr,, respectively. These can be divided into three pairs, each includes two probes of similar pK(a) in water and different lipophilicity. Solvation has been studied in binary mixtures, BMs, of water, W, with 12 protic organic solvents, S, including mono- and bifunctional alcohols (2-alkoxyethanoles, unsaturated and chlorinated alcohols). Each medium was treated as a mixture of S, W, and a complex solvent, S-W, formed by hydrogen bonding. Values of lambda(max) (of the probe intramolecular charge transfer) were converted into empirical polarity scales, E(T)(probe) in kcal/mol, whose values were correlated with the effective mole fraction of water in the medium, chi w(effective). This correlation furnished three equilibrium constants for the exchange of solvents in the probe solvation shell; phi(W/S) (W substitutes S): phi(S-W/W) (S-W substitutes W), and phi(S-W/S) (S-W substitutes S), respectively. The values of these constants depend on the physicochemical properties of the probe and the medium. We tested, for the first time, the applicability of a new solvation free energy relationship: phi = constant + a alpha(BM) + b beta(BM) + s(pi*(BM) + d delta) + p log P(BM), where a, b, s, and p are regression coefficients alpha(BM), beta(BM), and pi*(BM) are solvatochromic parameters of the BM, delta is a correction term for pi*, and log P is an empirical scale of lipophilicity. Correlations were carried out with two-, three-, and four-medium descriptors. In all cases, three descriptors gave satisfactory correlations; use of four parameters gave only a marginal increase of the goodness of fit. For phi(W/S), the most important descriptor was found to be the lipophilicity of the medium; for phi(S-W/W) and phi(S-W/S), solvent basicity is either statistically relevant or is the most important descriptor. These responses are different from those of E(T)(probe) of many solvatochromic indicators in pure solvents, where the importance of solvent basicity is usually marginal, and can be neglected.

Risk of recurrent stroke and death after first stroke in long-distance ski race participants

Background: Physical activity is of benefit for primary prevention of cardiovascular diseases, but it appears to increase the risk for atrial fibrillation. We aimed to study a cohort of patients following a first stroke in individuals with previous high physical activity, compare them to the general population with respect to recurrent stroke and death, and relate these to atrial fibrillation. Methods and results: From the participants of the Vasaloppet, the world's largest ski-race, and matched individuals from the general population (n=708 604), we identified 5964 patients hospitalized with a first-time stroke between 1994 and 2010. Individuals with severe diseases were excluded. One half percent of skiers and 1% of nonskiers were hospitalized due to stroke. The incidence rate was 8.3 per 100 person-years among skiers and 11.1 among nonskiers. The hazard ratio (HR) for recurrent stroke or death between the 2 groups was 0.76 (95% CI 0.67 to 0.86). The result was consistent in subgroups. The HR for death was 0.66 (95% CI 0.56 to 0.78) and for recurrent stroke 0.82 (95% CI 0.70 to 0.96). After adjustment for smoking and socioeconomic factors, the HR for death was consistent at 0.70 (95% CI 0.56 to 0.87) while the HR for recurrent stroke was not statistically significant. Outcomes for skiers with atrial fibrillation tended to show a lower risk than for nonskiers. Conclusions: This large cohort study supports the hypothesis that patients with a stroke and with prior regular physical activity have a lower risk of death, while their risk for recurrent stroke is similar to that of nonskiers. The skiers had a higher incidence of atrial fibrillation, but still no increased risk of recurring stroke.

In-depth search focused on furans, lactones, volatile phenols, and acetals as potential age markers of Madeira wines by comprehensive two-dimensional gas chromatography with time-of-flight mass spectrometry combined with solid phase microextraction

The establishment of potential age markers of Madeira wine is of paramount significance as it may contribute to detect frauds and to ensure the authenticity of wine. Considering the chemical groups of furans, lactones, volatile phenols, and acetals, 103 volatile compounds were tentatively identified; among these, 71 have been reported for the first time in Madeira wines. The chemical groups that could be used as potential age markers were predominantly acetals, namely, diethoxymethane, 1,1-diethoxyethane, 1,1-diethoxy-2-methyl-propane, 1-(1-ethoxyethoxy)-pentane, trans-dioxane and 2-propyl-1,3-dioxolane, and from the other chemical groups, 5-methylfurfural and cis-oak-lactone, independently of the variety and the type of wine. GC × GC-ToFMS system offers a more useful approach to identify these compounds compared to previous studies using GC−qMS, due to the orthogonal systems, that reduce coelution, increase peak capacity and mass selectivity, contributing to the establishment of new potential Madeira wine age markers. Remarkable results were also obtained in terms of compound identification based on the organized structure of the peaks of structurally related compounds in the GC × GC peak apex plots. This information represents a valuable approach for future studies, as the ordered-structure principle can considerably help the establishment of the composition of samples. This new approach provides data that can be extended to determine age markers of other types of wines.

First report of a B chromosome in a natural population of Astyanax altiparanae (Characiformes, Characidae)

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

A first case of protease codon 35 amino acid insertion in a HIV-1 subtype B sequence detected in the Bauru region, state of São Paulo, Brazil: case report

Inserções de aminoácidos na protease têm sido raramente descritas em pacientes infectados pelo HIV. Uma destas inserções foi, recentemente, descrita no codon 35, embora seu impacto na resistência mantém-se pouco conhecido. Este trabalho apresenta um caso de uma variante viral com inserção no codon 35 da protease, descrita pela primeira vez em Bauru, São Paulo, Brasil, circulante em um homem, caucasiano, com 38 anos, o qual apresenta infecção assintomática pelo HIV desde 1997. A variante isolada mostrou uma inserção no codon 35 da protease de dois aminoácidos: uma treonina e um ácido aspártico, resultando na sequência de aminoácidos E35E_TD.

The first reported case of human tick paralysis in Brazil: a new induction pattern by immature stages

Tick paralysis (TP) is a rare disease with rapid progression and potential fatal evolution. Immediately after the diagnosis, removal of all ticks from the body of the patient is mandatory. The present study reports for the first time a human case of the disease in Brazil. The patient had loss of muscle strength, decreased reflexes and marked palpebral ptosis. Six hours after removal of the last tick, the ptosis improved and on the following day, the patient had near total regression of the symptoms. This report emphasizes the possible presence of similar cases that should be promptly diagnosed and quickly treated. A new induction pattern for TP in humans associated with immature stages of ticks is also presented.

Intraspecific crosses resulting in the first occurrence of eight and nine B chromosomes in Prochilodus lineatus (Characiformes, Prochilodontidae)

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

Occurrence of sibling species of Lutzomyia longipalpis (Diptera : Psychodidae) in Venezuela: First evidence from reproductively isolated sympatric populations

The delimitation of cryptic species within the main vector of the American visceral leishmaniasis, Lutzomyia longipalpis, remains a topic of controversy. An analysis of generic variability based on 8 enzymatic loci revealed fixed differences in 2 diagnostic loci, adenylate kinase (Ak) and hexokinase (Hk), between sympatric and allopatric populations at 4 localities in Venezuela. The absence of heterozygotes for these 2 loci within 1 locality indicates, for the first time, the presence of 2 sympatric reproductively isolated populations or cryptic species within L. longipalpis. Significant differences were also detected between these cryptic species in the allele frequencies of glucose-6-phosphate isomerase (Gpi) and malate dehydrogenase, decarboxylating (Me). One species showed mean heterozygosities that ranged between 6.6% and 6.7%, with 1.6-1.9 alleles detected per locus, while the other had mean heterozygosities that ranged from 4.3% to 6.3%, with 1.3-1.6 alleles per locus. Comparisons of isozyme profiles with published data suggests that 1 species is similar to the L. longipalpis described in Colombian and Brazilian populations, whereas the other has not been previously reported.

Chromosomal mapping of repetitive DNAs in the beetle Dichotomius geminatus provides the first evidence for an association of 5S rRNA and histone H3 genes in insects, and repetitive DNA similarity between the B chromosome and A complement

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