Diagnosing observation error correlations for Doppler radar radial winds in the Met Office UKV model using observation-minus-background and observation-minus-analysis statistics

With the development of convection-permitting numerical weather prediction the efficient use of high resolution observations in data assimilation is becoming increasingly important. The operational assimilation of these observations, such as Dopplerradar radial winds, is now common, though to avoid violating the assumption of un- correlated observation errors the observation density is severely reduced. To improve the quantity of observations used and the impact that they have on the forecast will require the introduction of the full, potentially correlated, error statistics. In this work, observation error statistics are calculated for the Doppler radar radial winds that are assimilated into the Met Office high resolution UK model using a diagnostic that makes use of statistical averages of observation-minus-background and observation-minus-analysis residuals. This is the first in-depth study using the diagnostic to estimate both horizontal and along-beam correlated observation errors. By considering the new results obtained it is found that the Doppler radar radial wind error standard deviations are similar to those used operationally and increase as the observation height increases. Surprisingly the estimated observation error correlation length scales are longer than the operational thinning distance. They are dependent on both the height of the observation and on the distance of the observation away from the radar. Further tests show that the long correlations cannot be attributed to the use of superobservations or the background error covariance matrix used in the assimilation. The large horizontal correlation length scales are, however, in part, a result of using a simplified observation operator.

Tests of sunspot number sequences: 1. Using ionosonde data

More than 70 years ago it was recognised that ionospheric F2-layer critical frequencies [foF2] had a strong relationship to sunspot number. Using historic datasets from the Slough and Washington ionosondes, we evaluate the best statistical fits of foF2 to sunspot numbers (at each Universal Time [UT] separately) in order to search for drifts and abrupt changes in the fit residuals over Solar Cycles 17-21. This test is carried out for the original composite of the Wolf/Zürich/International sunspot number [R], the new “backbone” group sunspot number [RBB] and the proposed “corrected sunspot number” [RC]. Polynomial fits are made both with and without allowance for the white-light facular area, which has been reported as being associated with cycle-to-cycle changes in the sunspot number - foF2 relationship. Over the interval studied here, R, RBB, and RC largely differ in their allowance for the “Waldmeier discontinuity” around 1945 (the correction factor for which for R, RBB and RC is, respectively, zero, effectively over 20 %, and explicitly 11.6 %). It is shown that for Solar Cycles 18-21, all three sunspot data sequences perform well, but that the fit residuals are lowest and most uniform for RBB. We here use foF2 for those UTs for which R, RBB, and RC all give correlations exceeding 0.99 for intervals both before and after the Waldmeier discontinuity. The error introduced by the Waldmeier discontinuity causes R to underestimate the fitted values based on the foF2 data for 1932-1945 but RBB overestimates them by almost the same factor, implying that the correction for the Waldmeier discontinuity inherent in RBB is too large by a factor of two. Fit residuals are smallest and most uniform for RC and the ionospheric data support the optimum discontinuity multiplicative correction factor derived from the independent Royal Greenwich Observatory (RGO) sunspot group data for the same interval.

An idealized study of coupled atmosphere-ocean 4D-Var in the presence of model error

Atmosphere only and ocean only variational data assimilation (DA) schemes are able to use window lengths that are optimal for the error growth rate, non-linearity and observation density of the respective systems. Typical window lengths are 6-12 hours for the atmosphere and 2-10 days for the ocean. However, in the implementation of coupled DA schemes it has been necessary to match the window length of the ocean to that of the atmosphere, which may potentially sacrifice the accuracy of the ocean analysis in order to provide a more balanced coupled state. This paper investigates how extending the window length in the presence of model error affects both the analysis of the coupled state and the initialized forecast when using coupled DA with differing degrees of coupling. Results are illustrated using an idealized single column model of the coupled atmosphere-ocean system. It is found that the analysis error from an uncoupled DA scheme can be smaller than that from a coupled analysis at the initial time, due to faster error growth in the coupled system. However, this does not necessarily lead to a more accurate forecast due to imbalances in the coupled state. Instead coupled DA is more able to update the initial state to reduce the impact of the model error on the accuracy of the forecast. The effect of model error is potentially most detrimental in the weakly coupled formulation due to the inconsistency between the coupled model used in the outer loop and uncoupled models used in the inner loop.

Dynamic portrait of the planetary 2/1 mean-motion resonance - I. Systems with a more massive outer planet

We analyse the global structure of the phase space of the planar planetary 2/1 mean-motion resonance in cases where the outer planet is more massive than its inner companion. Inside the resonant domain, we show the existence of two families of periodic orbits, one associated to the librational motion of resonant angle (sigma-family) and the other related to the circulatory motion of the difference in longitudes of pericentre (Delta pi-family). The well-known apsidal corotation resonances (ACR) appear as intersections between both families. A complex web of secondary resonances is also detected for low eccentricities, whose strengths and positions are dependent on the individual masses and spatial scale of the system. The construction of dynamical maps for various values of the total angular momentum shows the evolution of the families of stable motion with the eccentricities, identifying possible configurations suitable for exoplanetary systems. For low-moderate eccentricities, several different stable modes exist outside the ACR. For larger eccentricities, however, all stable solutions are associated to oscillations around the stationary solutions. Finally, we present a possible link between these stable families and the process of resonance capture, identifying the most probable routes from the secular region to the resonant domain, and discussing how the final resonant configuration may be affected by the extension of the chaotic layer around the resonance region.

DETECTABILITY AND ERROR ESTIMATION IN ORBITAL FITS OF RESONANT EXTRASOLAR PLANETS

We estimate the conditions for detectability of two planets in a 2/1 mean-motion resonance from radial velocity data, as a function of their masses, number of observations and the signal-to-noise ratio. Even for a data set of the order of 100 observations and standard deviations of the order of a few meters per second, we find that Jovian-size resonant planets are difficult to detect if the masses of the planets differ by a factor larger than similar to 4. This is consistent with the present population of real exosystems in the 2/1 commensurability, most of which have resonant pairs with similar minimum masses, and could indicate that many other resonant systems exist, but are currently beyond the detectability limit. Furthermore, we analyze the error distribution in masses and orbital elements of orbital fits from synthetic data sets for resonant planets in the 2/1 commensurability. For various mass ratios and number of data points we find that the eccentricity of the outer planet is systematically overestimated, although the inner planet`s eccentricity suffers a much smaller effect. If the initial conditions correspond to small-amplitude oscillations around stable apsidal corotation resonances, the amplitudes estimated from the orbital fits are biased toward larger amplitudes, in accordance to results found in real resonant extrasolar systems.

Alkaline phosphatase as a reporter of sigma(S) levels and rpoS polymorphisms in different E-coli strains

sigma(S) is responsible for the transcriptional regulation of genes related to protection against stresses and bacterial survival and it accumulates in the cell under conditions of stress, such as nutrient limitation. An increase in the levels of sigma(S) causes a reduction in the expression of genes that are transcribed by RNA polymerase associated with the principal sigma factor, sigma(70). phoA, that encodes alkaline phosphatase (AP) is expressed under phosphate shortage conditions, and is also repressed by sigma(S). Here we show that in a Pi-limited chemostat, accumulation of rpoS mutations is proportional to the intrinsic level of sigma(S) in the cells. Acquisition of mutations in rpoS relieves repression of the PHO genes. We also devised a non-destructive method based on the rpoS effect on AP that differentiates between rpo(S+) and rpoS mutants, as well as between high and low-sigma(S) producers. Using this method, we provide evidence that sigma(S) contributes to the repression of AP under conditions of Pi excess and that AP variation among different strains is at least partly due to intrinsic variation in sigma(S) levels. Consequently, a simple and non-destructive AP assay can be employed to differentiate between strains expressing different levels of sigma(S) on agar plates.

SIGMA THEORY AND TWISTED CONJUGACY CLASSES

Using Sigma theory we show that for large classes of groups G there is a subgroup H of finite index in Aut(G) such that for phi is an element of H the Reidemeister number R(phi) is infinite. This includes all finitely generated nonpolycyclic groups G that fall into one of the following classes: nilpotent-by-abelian groups of type FP(infinity); groups G/G `` of finite Prufer rank; groups G of type FP(2) without free nonabelian subgroups and with nonpolycyclic maximal metabelian quotient; some direct products of groups; or the pure symmetric automorphism group. Using a different argument we show that the result also holds for 1-ended nonabelian nonsurface limit groups. In some cases, such as with the generalized Thompson`s groups F(n,0) and their finite direct products, H = Aut(G).

(1)H NMR spectroscopy as a tool to determine accurate photoisomerization quantum yields of stilbene-like ligands coordinated to rhenium(I) polypyridyl complexes

In this work, the use of proton nuclear magnetic resonance, (1)H NMR, was fully described as a powerful tool to follow a photoreaction and to determine accurate quantum yields, so called true quantum yields (Phi(true)), when a reactant and photoproduct absorption overlap. For this, Phi(true) for the trans-cis photoisomerization process were determined for rhenium(I) polypyridyl complexes, fac-[Re(CO)(3)(NN)(trans-L)](+) (NN = 1,10-phenanthroline, phen, or 4,7-diphenyl-1,10-phenanthroline, ph(2)phen, and L = 1,2-bis(4-pyridyl) ethylene, bpe, or 4-styrylpyridine, stpy). The true values determined at 365 nm irradiation (e. g. Phi(NMR) = 0.80 for fac-[Re(CO)(3)(phen)(trans-bpe)](+)) were much higher than those determined by absorption spectral changes (Phi(UV-Vis) = 0.39 for fac-[Re(CO)(3)(phen)(trans-bpe)](+)). Phi(NMR) are more accurate in these cases due to the distinct proton signals of trans and cis-isomers, which allow the actual determination of each component concentration under given irradiation time. Nevertheless when the photoproduct or reactant contribution at the probe wavelength is negligible, one can determine Phi(true) by regular absorption spectral changes. For instance, Phi(313) nm for free ligand photoisomerization determined both by absorption and (1)H NMR variation are equal within the experimental error (bpe: Phi(UV-Vis) = 0.27, Phi(NMR) = 0.26; stpy: Phi(UV-Vis) = 0.49, Phi(NMR) = 0.49). Moreover, (1)H NMR data combined with electronic spectra allowed molar absorptivity determination of difficult to isolate cis-complexes. (C) 2009 Elsevier B. V. All rights reserved.

What happens after the central bank of Brazil increases the target interbank rate by 1%?

I start presenting an explicit solution to Taylorís (2001) model, in order to illustrate the link between the target interest rate and the overnight interest rate prevailing in the economy. Next, I use Vector Auto Regressions to shed some light on the evolution of key macroeconomic variables after the Central Bank of Brazil increases the target interest rate by 1%. Point estimates show a four-year accumulated output loss ranging from 0:04% (whole sample, 1980 : 1-2004 : 2; quarterly data) to 0:25% (Post-Real data only) with a Örst-year peak output response between 0:04% and 1:0%; respectively. Prices decline between 2% and 4% in a 4-year horizon. The accumulated output response is found to be between 3:5 and 6 times higher after the Real Plan than when the whole sample is considered. The 95% confidence bands obtained using bias-corrected bootstrap always include the null output response when the whole sample is used, but not when the data is restricted to the Post-Real period. Innovations to interest rates explain between 4:9% (whole sample) and 9:2% (post-Real sample) of the forecast error of GDP.

Práticas de avaliação nas escolas oficiais de 1. grau, em Manaus

Este trabalho visa a examinar as práticas de avaliação das escolas oficiais da cidade de Manaus que oferecem as quatro primeiras séries do Ensino de 1o Grau, com vistas a averiguar se a execução da avaliação está de acordo com as prescrições legais em vigor. Trata-se de estudo piloto exploratório feito no período de agosto a setembro de 1977, em escolas selecionadas por amostragem estratificada proporcional, tendo sido estudadas 50% das escolas da população alvo (em número de 25) nas quais foram envolvidos os diretores e os supervisores, bem como, 300 professores selecionados pelos diretores. Na coleta de dados foram utilizados dois questionários diferentes, um aplicado aos professores e outro aos supervisores, fazendo-se ainda uma entrevista individual semi-estruturada com os diretores das escolas envolvidas. Para reduzir as fontes de erro, a própria autora entregou e recolheu pessoalmente os questionários (obtendo uma taxa de retorno de 77%) e conduziu todas as entrevistas. A análise de dados focalizou as seguintes variáveis: o preparo do professor para realizar a avaliação prescrita pela legislação; as práticas de avaliação empregadas pelos professores; o tipo de orientação dada pelo supervisor ao professor no que se refere à avaliação da aprendizagem; e as condições oferecidas ao professor, ao supervisor e às escolas para a prática da avaliação. Os resultados indicaram que, de modo geral, as práticas de avaliação escolar não correspondem aos padrões exigidos pelos textos legais. A deficiência de informações adequadas sobre a teoria atual de avaliação e sobre instrumentos e procedimentos, bem como a insuficiência de condições de prática necessárias à execução da avaliação nos moldes recomendados destacam-se como questões chave do problema.

Modelagem não-paramétrica da dinâmica da taxa de juros instantânea utilizando contratos futuros da taxa média dos depósitos interfinanceiros de 1 dia (DI1)

Modelos de predição baseados em estimações não-paramétricas continuam em desenvolvimento e têm permeado a comunidade quantitativa. Sua principal característica é que não consideram a priori distribuições de probabilidade conhecidas, mas permitem que os dados passados sirvam de base para a construção das próprias distribuições. Implementamos para o mercado brasileiro os estimadores agrupados não-paramétricos de Sam e Jiang (2009) para as funções de drift e de difusão do processo estocástico da taxa de juros instantânea, por meio do uso de séries de taxas de juros de diferentes maturidades fornecidas pelos contratos futuros de depósitos interfinanceiros de um dia (DI1). Os estimadores foram construídos sob a perspectiva da estimação por núcleos (kernels), que requer para a sua otimização um formato específico da função-núcleo. Neste trabalho, foi usado o núcleo de Epanechnikov, e um parâmetro de suavizamento (largura de banda), o qual é fundamental para encontrar a função de densidade de probabilidade ótima que forneça a estimação mais eficiente em termos do MISE (Mean Integrated Squared Error - Erro Quadrado Integrado Médio) no momento de testar o modelo com o tradicional método de validação cruzada de k-dobras. Ressalvas são feitas quando as séries não possuem os tamanhos adequados, mas a quebra estrutural do processo de difusão da taxa de juros brasileira, a partir do ano 2006, obriga à redução do tamanho das séries ao custo de reduzir o poder preditivo do modelo. A quebra estrutural representa um processo de amadurecimento do mercado brasileiro que provoca em grande medida o desempenho insatisfatório do estimador proposto.

Dois critérios de seleção na pré-desmama em bovinos da raça Gir. 1. Estimativas de parâmetros genéticos

Dados de 16.592 animais da raça Gir, provenientes do controle de desenvolvimento ponderal da Associação Brasileira dos Criadores de Zebu, nascidos no período de 1978 a 1994, criados em diversas condições de ambiente no Brasil, foram usados para estimar parâmetros genéticos e ambientais das características ganho médio diário do nascimento à desmama (GMD) e dias para atingir 160 kg do nascimento à desmama (D160). As estimativas dos componentes de variância e herdabilidade foram obtidas pelo método da máxima verossimilhança restrita por modelos animais uni-característica. Foram utilizados dois modelos o modelo 1, com grupo de contemporâneos (GC) como efeito fixo e com efeitos genético aditivo de animal e materno, efeito de ambiente permanente materno e erro como aleatórios; e modelo 2, com os mesmos efeitos do modelo 1, mas com covariância genética entre os efeitos direto e materno igual a zero (sigmaam=0). As estimativas de herdabilidade, utilizando-se o Modelo 1, foram iguais a 0,11; 0,04; 0,11 e 0,07 para GMD (efeito direto), GMD (efeito materno), D160 (efeito direto) e D160 (efeito materno), respectivamente. Para o Modelo 2, as estimativas de herdabilidade foram: 0,12; 0,05; 0,10 e 0,05 para GMD (efeito direto), GMD (efeito materno), D160 (efeito direto) e D160 (efeito materno). A diferença entre os modelos pela não inclusão da estrutura de covariância entre os efeitos direto e materno (sigmaam=0), neste trabalho, não foi significativa. A contribuição da variância residual (como porcentagem da variação fenotípica total) foi de 76% para GMD e de 80% para D160, para os dois modelos, indicando a necessidade da melhoria do ambiente de criação, assim como a padronização e formação de lotes de manejo, o que reflete na possibilidade de melhor definição de grupo de contemporâneos.

Improving the GNSS positioning stochastic model in the presence of ionospheric scintillation

Ionospheric scintillations are caused by time-varying electron density irregularities in the ionosphere, occurring more often at equatorial and high latitudes. This paper focuses exclusively on experiments undertaken in Europe, at geographic latitudes between similar to 50 degrees N and similar to 80 degrees N, where a network of GPS receivers capable of monitoring Total Electron Content and ionospheric scintillation parameters was deployed. The widely used ionospheric scintillation indices S4 and sigma(phi) represent a practical measure of the intensity of amplitude and phase scintillation affecting GNSS receivers. However, they do not provide sufficient information regarding the actual tracking errors that degrade GNSS receiver performance. Suitable receiver tracking models, sensitive to ionospheric scintillation, allow the computation of the variance of the output error of the receiver PLL (Phase Locked Loop) and DLL (Delay Locked Loop), which expresses the quality of the range measurements used by the receiver to calculate user position. The ability of such models of incorporating phase and amplitude scintillation effects into the variance of these tracking errors underpins our proposed method of applying relative weights to measurements from different satellites. That gives the least squares stochastic model used for position computation a more realistic representation, vis-a-vis the otherwise 'equal weights' model. For pseudorange processing, relative weights were computed, so that a 'scintillation-mitigated' solution could be performed and compared to the (non-mitigated) 'equal weights' solution. An improvement between 17 and 38% in height accuracy was achieved when an epoch by epoch differential solution was computed over baselines ranging from 1 to 750 km. The method was then compared with alternative approaches that can be used to improve the least squares stochastic model such as weighting according to satellite elevation angle and by the inverse of the square of the standard deviation of the code/carrier divergence (sigma CCDiv). The influence of multipath effects on the proposed mitigation approach is also discussed. With the use of high rate scintillation data in addition to the scintillation indices a carrier phase based mitigated solution was also implemented and compared with the conventional solution. During a period of occurrence of high phase scintillation it was observed that problems related to ambiguity resolution can be reduced by the use of the proposed mitigated solution.

A Discussion Related to Orbit Determination Using Nonlinear Sigma Point Kalman Filter

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