Selecting algorithms for Earth observation of climate within the European Space Agency Climate Change initiative: introduction to a special issue

This special issue is focused on the assessment of algorithms for the observation of Earth’s climate from environ- mental satellites. Climate data records derived by remote sensing are increasingly a key source of insight into the workings of and changes in Earth’s climate system. Producers of data sets must devote considerable effort and expertise to maximise the true climate signals in their products and minimise effects of data processing choices and changing sensors. A key choice is the selection of algorithm(s) for classification and/or retrieval of the climate variable. Within the European Space Agency Climate Change Initiative, science teams undertook systematic assessment of algorithms for a range of essential climate variables. The papers in the special issue report some of these exercises (for ocean colour, aerosol, ozone, greenhouse gases, clouds, soil moisture, sea surface temper- ature and glaciers). The contributions show that assessment exercises must be designed with care, considering issues such as the relative importance of different aspects of data quality (accuracy, precision, stability, sensitivity, coverage, etc.), the availability and degree of independence of validation data and the limitations of validation in characterising some important aspects of data (such as long-term stability or spatial coherence). As well as re- quiring a significant investment of expertise and effort, systematic comparisons are found to be highly valuable. They reveal the relative strengths and weaknesses of different algorithmic approaches under different observa- tional contexts, and help ensure that scientific conclusions drawn from climate data records are not influenced by observational artifacts, but are robust.

A Time-Series Analysis of the H alpha Emission Line in V3885 Sagitarii

Flickering is a phenomenon related to mass accretion observed among many classes of astrophysical objects. In this paper we present a study of flickering emission lines and the continuum of the cataclysmic variable V3885 Sgr. The flickering behavior was first analyzed through statistical analysis and the power spectra of lightcurves. Autocorrelation techniques were then employed to estimate the flickering timescale of flares. A cross-correlation study between the line and its underlying continuum variability is presented. The cross-correlation between the photometric and spectroscopic data is also discussed. Periodograms, calculated using emission-line data, show a behavior that is similar to those obtained from photometric datasets found in the literature, with a plateau at lower frequencies and a power-law at higher frequencies. The power-law index is consistent with stochastic events. The cross-correlation study indicates the presence of a correlation between the variability on Ha and its underlying continuum. Flickering timescales derived from the photometric data were estimated to be 25 min for two lightcurves and 10 min for one of them. The average timescales of the line flickering is 40 min, while for its underlying continuum it drops to 20 min.

Fluorescence spectroscopy to diagnose hepatic steatosis in a rat model of fatty liver

Steatosis is diagnosed on the basis of the macroscopic aspect of the liver evaluated by the surgeon at the time of organ extraction or by means of a frozen biopsy. In the present study, the applicability of laser-induced fluorescence (LIF) spectroscopy was investigated as a method for the diagnosis of different degrees of steatosis experimentally induced in rats. Rats received a high-lipid diet for different periods of time. The animals were divided into groups according to the degree of induced steatosis diagnosis by histology. The concentration of fat in the liver was correlated with LIF by means of the steatosis fluorescence factor (SFF). The histology classification, according to liver fat concentration was, Severe Steatosis, Moderate Steatosis, Mild Steatosis and Control (no liver steatosis). Fluorescence intensity could be directly correlated with fat content. It was possible to estimate an average of fluorescence intensity variable by means of different confidence intervals (P=95%) for each steatosis group. SFF was significantly higher in the Severe Steatosis group (P < 0.001) compared with the Moderate Steatosis, Mild Steatosis and Control groups. The various degrees of steatosis could be directly correlated with SFF. LIF spectroscopy proved to be a method capable of identifying the degree of hepatic steatosis in this animal model, and has the potential of clinical application for non-invasive evaluation of the degree of steatosis.

Skew-symmetric identities of octonions

Quadratic alternative superalgebras are introduced and their super-identities and central functions on one odd generator are described. As a corollary, all multilinear skew-symmetric identities and central polynomials of octonions are classified. (C) 2008 Elsevier B.V. All rights reserved.

Mutagenicity and DNA adduct formation of PAH, nitro-PAH, and oxy-PAH fractions of atmospheric particulate matter from Sao Paulo, Brazil

Urban particulate matter (UPM) contributes to lung cancer incidence. Here, we have studied the mutagenic activity and DNA adduct-forming ability of fractionated UPM extractable organic matter (EOM). UPM was collected with a high-volume sampler in June 2004 at two sites, one at street level adjacent to a roadway and the other inside a park within the urban area of the city of Sao Paulo, Brazil. UPM was extracted using dichloromethane, and the resulting EOM was separated by HPLC to obtain PAH, nitro-PAH, and oxy-PAH fractions which were tested for mutagenicity with the Salmonella strains TA98 and YG1041 with and without S9 metabolic activation. The PAH fraction from both sites showed negligible mutagenic activity in both strains. The highest mutagenic activity was found for the nitro-PAH fraction using YG1041 without metabolic activation; however, results were comparable for both sites. The nitro-PAH and oxy-PAH fractions were incubated with calf thymus DNA under reductive conditions appropriate for the activation of nitro aromatic compounds, then DNA adduct patterns and levels were determined with thin-layer chromatography (TLC) (32)p-postlabeling method using two enrichment procedures-nuclease PI digestion and butanol extraction. Reductively activated fractions from both sites produced diagonal radioactive zones (DRZ) of putative aromatic DNA adducts on thin layer plates with both enrichment procedures. No such DRZ were observed in control experiments using fractions from unexposed filters or from incubations without activating system. Total adduct levels produced by the nitro-PAH fractions were similar for both sites ranging from 30 to 45 adducts per 10(8) normal nucleotides. In contrast, the DNA binding of reductively activated oxy-PAH fractions was three times higher and the adduct pattern consisted of multiple discrete spots along the diagonal line on the thin layer plates. However, DNA adduct levels were not significantly different between the sampling sites. Both samples presented the same levels of mutagenic activity. The response in the Salmonella assay was typical of nitroaromatics. Although, more mutagenic activity was related to the nitro-PAH fraction in the Salmonella assay, the oxy-PAH fractions showed the highest DNA adduct levels. More studies are needed to elucidate the nature of the genotoxicants occurring in Sao Paulo atmospheric samples. (C) 2008 Elsevier B.V. All rights reserved.

Meteomap: Generation Of Meteorological Variables For Distributed Physically-Based Hydrological Modeling

Distributed energy and water balance models require time-series surfaces of the meteorological variables involved in hydrological processes. Most of the hydrological GIS-based models apply simple interpolation techniques to extrapolate the point scale values registered at weather stations at a watershed scale. In mountainous areas, where the monitoring network ineffectively covers the complex terrain heterogeneity, simple geostatistical methods for spatial interpolation are not always representative enough, and algorithms that explicitly or implicitly account for the features creating strong local gradients in the meteorological variables must be applied. Originally developed as a meteorological pre-processing tool for a complete hydrological model (WiMMed), MeteoMap has become an independent software. The individual interpolation algorithms used to approximate the spatial distribution of each meteorological variable were carefully selected taking into account both, the specific variable being mapped, and the common lack of input data from Mediterranean mountainous areas. They include corrections with height for both rainfall and temperature (Herrero et al., 2007), and topographic corrections for solar radiation (Aguilar et al., 2010). MeteoMap is a GIS-based freeware upon registration. Input data include weather station records and topographic data and the output consists of tables and maps of the meteorological variables at hourly, daily, predefined rainfall event duration or annual scales. It offers its own pre and post-processing tools, including video outlook, map printing and the possibility of exporting the maps to images or ASCII ArcGIS formats. This study presents the friendly user interface of the software and shows some case studies with applications to hydrological modeling.

Identification and estimation of interventions using changes in inequality measures

This paper presents semiparametric estimators of changes in inequality measures of a dependent variable distribution taking into account the possible changes on the distributions of covariates. When we do not impose parametric assumptions on the conditional distribution of the dependent variable given covariates, this problem becomes equivalent to estimation of distributional impacts of interventions (treatment) when selection to the program is based on observable characteristics. The distributional impacts of a treatment will be calculated as differences in inequality measures of the potential outcomes of receiving and not receiving the treatment. These differences are called here Inequality Treatment Effects (ITE). The estimation procedure involves a first non-parametric step in which the probability of receiving treatment given covariates, the propensity-score, is estimated. Using the inverse probability weighting method to estimate parameters of the marginal distribution of potential outcomes, in the second step weighted sample versions of inequality measures are computed. Root-N consistency, asymptotic normality and semiparametric efficiency are shown for the semiparametric estimators proposed. A Monte Carlo exercise is performed to investigate the behavior in finite samples of the estimator derived in the paper. We also apply our method to the evaluation of a job training program.

Optimal two-sided tests for instrumental variables regression with heteroskedastic and autocorrelated errors

This paper considers two-sided tests for the parameter of an endogenous variable in an instrumental variable (IV) model with heteroskedastic and autocorrelated errors. We develop the nite-sample theory of weighted-average power (WAP) tests with normal errors and a known long-run variance. We introduce two weights which are invariant to orthogonal transformations of the instruments; e.g., changing the order in which the instruments appear. While tests using the MM1 weight can be severely biased, optimal tests based on the MM2 weight are naturally two-sided when errors are homoskedastic. We propose two boundary conditions that yield two-sided tests whether errors are homoskedastic or not. The locally unbiased (LU) condition is related to the power around the null hypothesis and is a weaker requirement than unbiasedness. The strongly unbiased (SU) condition is more restrictive than LU, but the associated WAP tests are easier to implement. Several tests are SU in nite samples or asymptotically, including tests robust to weak IV (such as the Anderson-Rubin, score, conditional quasi-likelihood ratio, and I. Andrews' (2015) PI-CLC tests) and two-sided tests which are optimal when the sample size is large and instruments are strong. We refer to the WAP-SU tests based on our weights as MM1-SU and MM2-SU tests. Dropping the restrictive assumptions of normality and known variance, the theory is shown to remain valid at the cost of asymptotic approximations. The MM2-SU test is optimal under the strong IV asymptotics, and outperforms other existing tests under the weak IV asymptotics.

Generalized linear mixed models for the genetic evaluation of binary reproductive traits: a simulation study

The objective of this study was to evaluate the use of probit and logit link functions for the genetic evaluation of early pregnancy using simulated data. The following simulation/analysis structures were constructed: logit/logit, logit/probit, probit/logit, and probit/probit. The percentages of precocious females were 5, 10, 15, 20, 25 and 30% and were adjusted based on a change in the mean of the latent variable. The parametric heritability (h²) was 0.40. Simulation and genetic evaluation were implemented in the R software. Heritability estimates (ĥ²) were compared with h² using the mean squared error. Pearson correlations between predicted and true breeding values and the percentage of coincidence between true and predicted ranking, considering the 10% of bulls with the highest breeding values (TOP10) were calculated. The mean ĥ² values were under- and overestimated for all percentages of precocious females when logit/probit and probit/logit models used. In addition, the mean squared errors of these models were high when compared with those obtained with the probit/probit and logit/logit models. Considering ĥ², probit/probit and logit/logit were also superior to logit/probit and probit/logit, providing values close to the parametric heritability. Logit/probit and probit/logit presented low Pearson correlations, whereas the correlations obtained with probit/probit and logit/logit ranged from moderate to high. With respect to the TOP10 bulls, logit/probit and probit/logit presented much lower percentages than probit/probit and logit/logit. The genetic parameter estimates and predictions of breeding values of the animals obtained with the logit/logit and probit/probit models were similar. In contrast, the results obtained with probit/logit and logit/probit were not satisfactory. There is need to compare the estimation and prediction ability of logit and probit link functions.

Covariance functions for body weight from birth to maturity in Nellore cows

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Estimação de componentes de co-variância para pesos corporais do nascimento aos 365 dias de idade de bovinos Guzerá empregando-se modelos de regressão aleatória

Um total de 19.770 pesos corporais de bovinos Guzerá, do nascimento aos 365 dias de idade, pertencentes ao banco de dados da Associação Brasileira dos Criadores de Zebu (ABCZ) foi analisado com os objetivos de comparar diferentes estruturas de variâncias residuais, considerando 1, 18, 28 e 53 classes residuais e funções de variância de ordens quadrática a quíntica; e estimar funções de co-variância de diferentes ordens para os efeitos genético aditivo direto, genético materno, de ambiente permanente de animal e de mãe e parâmetros genéticos para os pesos corporais usando modelos de regressão aleatória. Os efeitos aleatórios foram modelados por regressões polinomiais em escala de Legendre com ordens variando de linear a quártica. Os modelos foram comparados pelo teste de razão de verossimilhança e pelos critérios de Informação de Akaike e de Informação Bayesiano de Schwarz. O modelo com 18 classes heterogêneas foi o que melhor se ajustou às variâncias residuais, de acordo com os testes estatísticos, porém, o modelo com função de variância de quinta ordem também mostrou-se apropriado. Os valores de herdabilidade direta estimados foram maiores que os encontrados na literatura, variando de 0,04 a 0,53, mas seguiram a mesma tendência dos estimados pelas análises unicaracterísticas. A seleção para peso em qualquer idade melhoraria o peso em todas as idades no intervalo estudado.

Sparse one-dimensional discrete Dirac operators II: Spectral properties

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

Molecular typing and virulence markers of Yersinia enterocolitica strains from human, animal and food origins isolated between 1968 and 2000 in Brazil

Molecular typing and virulence markers were used to evaluate the genetic profiles and virulence potential of 106 Yersinia enterocolitica strains. of these strains, 71 were bio-serotype 4/O: 3, isolated from human and animal clinical material, and 35 were of biotype 1 A or 2 and of diverse serotypes, isolated from food in Brazil between 1968 and 2000. Drug resistance was also investigated. All the strains were resistant to three or more drugs. The isolates showed a virulence-related phenotype in the aesculin, pyrazinamidase and salicin tests, except for the food isolates, only two of which were positive for these tests. For the other phenotypic virulence determinants (autoagglutination, Ca++ dependence and Congo red absorption), the strains showed a diverse behaviour. The inv, ail and ystA genes were detected in all human and animal strains, while all the food isolates were positive for inv, and 3% of them positive for ail and ystA. The presence of virF was variable in the three groups of strains. The strains were better discriminated by PFGE than by enterobacterial repetitive intergenic consensus PCR (ERIC-PCR). A higher genomic similarity was observed among the 4/O: 3 strains, isolated from human and animal isolates, than among the food strains, with the exception of two food strains possessing the virulence genes and grouped close to the 4/O: 3 strains by ERIC-PCR. Unusually, the results revealed the virulence potential of a bio-serotype 1 A/O: 10 strain, suggesting that food contaminated with Y. enterocolitica biotype 1 A may cause infection. This also suggests that ERIC-PCR may be used as a tool to reveal clues about the virulence potential of Y. enterocolitica strains. Furthermore, the results also support the hypothesis that animals may act as reservoirs of Y. enterocolitica for human infections in Brazil, an epidemiological aspect that has not been investigated in this country, confirming data from other parts of the world.

Stress induces behavioral sensitization, increases nicotine-seeking behavior and leads to a decrease of CREB in the nucleus accumbens

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

Stability of neutral functional differential equations in terms of measures

In this paper we investigate the relationships between different concepts of stability in measure for the solutions of an autonomous or periodic neutral functional differential equation.