Assimilation of sea-ice concentration in a global climate model — physical and statistical aspects

We investigate the initialization of Northern-hemisphere sea ice in the global climate model ECHAM5/MPI-OM by assimilating sea-ice concentration data. The analysis updates for concentration are given by Newtonian relaxation, and we discuss different ways of specifying the analysis updates for mean thickness. Because the conservation of mean ice thickness or actual ice thickness in the analysis updates leads to poor assimilation performance, we introduce a proportional dependence between concentration and mean thickness analysis updates. Assimilation with these proportional mean-thickness analysis updates significantly reduces assimilation error both in identical-twin experiments and when assimilating sea-ice observations, reducing the concentration error by a factor of four to six, and the thickness error by a factor of two. To understand the physical aspects of assimilation errors, we construct a simple prognostic model of the sea-ice thermodynamics, and analyse its response to the assimilation. We find that the strong dependence of thermodynamic ice growth on ice concentration necessitates an adjustment of mean ice thickness in the analysis update. To understand the statistical aspects of assimilation errors, we study the model background error covariance between ice concentration and ice thickness. We find that the spatial structure of covariances is best represented by the proportional mean-thickness analysis updates. Both physical and statistical evidence supports the experimental finding that proportional mean-thickness updates are superior to the other two methods considered and enable us to assimilate sea ice in a global climate model using simple Newtonian relaxation.

Assimilation of sea-ice concentration in a global climate model — physical and statistical aspects

We investigate the initialisation of Northern Hemisphere sea ice in the global climate model ECHAM5/MPI-OM by assimilating sea-ice concentration data. The analysis updates for concentration are given by Newtonian relaxation, and we discuss different ways of specifying the analysis updates for mean thickness. Because the conservation of mean ice thickness or actual ice thickness in the analysis updates leads to poor assimilation performance, we introduce a proportional dependence between concentration and mean thickness analysis updates. Assimilation with these proportional mean-thickness analysis updates leads to good assimilation performance for sea-ice concentration and thickness, both in identical-twin experiments and when assimilating sea-ice observations. The simulation of other Arctic surface fields in the coupled model is, however, not significantly improved by the assimilation. To understand the physical aspects of assimilation errors, we construct a simple prognostic model of the sea-ice thermodynamics, and analyse its response to the assimilation. We find that an adjustment of mean ice thickness in the analysis update is essential to arrive at plausible state estimates. To understand the statistical aspects of assimilation errors, we study the model background error covariance between ice concentration and ice thickness. We find that the spatial structure of covariances is best represented by the proportional mean-thickness analysis updates. Both physical and statistical evidence supports the experimental finding that assimilation with proportional mean-thickness updates outperforms the other two methods considered. The method described here is very simple to implement, and gives results that are sufficiently good to be used for initialising sea ice in a global climate model for seasonal to decadal predictions.

Thermodynamic and statistical analysis of soybean oil extraction process using renewable solvent

The main goal of this work was to evaluate thermodynamic parameters of the soybean oil extraction process using ethanol as solvent. The experimental treatments were as follows: aqueous solvents with water contents varying from 0 to 13% (mass basis) and extraction temperature varying from 50 to 100 degrees C. The distribution coefficients of oil at equilibrium have been used to calculate enthalpy, entropy and free energy changes. The results indicate that oil extraction process with ethanol is feasible and spontaneous, mainly under higher temperature. Also, the influence of water level in the solvent and temperature were analysed using the response surface methodology (RSM). It can be noted that the extraction yield was highly affected by both independent variables. A joint analysis of thermodynamic and RSM indicates the optimal level of solvent hydration and temperature to perform the extraction process.

Crystal structure and statistical coupling analysis of highly glycosylated peroxidase from royal palm tree (Roystonea regia)

Royal palm tree peroxidase (RPTP) is a very stable enzyme in regards to acidity, temperature, H(2)O(2), and organic solvents. Thus, RPTP is a promising candidate for developing H(2)O(2)-sensitive biosensors for diverse applications in industry and analytical chemistry. RPTP belongs to the family of class III secretory plant peroxidases, which include horseradish peroxidase isozyme C, soybean and peanut peroxidases. Here we report the X-ray structure of native RPTP isolated from royal palm tree (Roystonea regia) refined to a resolution of 1.85 angstrom. RPTP has the same overall folding pattern of the plant peroxidase superfamily, and it contains one heme group and two calcium-binding sites in similar locations. The three-dimensional structure of RPTP was solved for a hydroperoxide complex state, and it revealed a bound 2-(N-morpholino) ethanesulfonic acid molecule (MES) positioned at a putative substrate-binding secondary site. Nine N-glycosylation sites are clearly defined in the RPTP electron-density maps, revealing for the first time conformations of the glycan chains of this highly glycosylated enzyme. Furthermore, statistical coupling analysis (SCA) of the plant peroxidase superfamily was performed. This sequence-based method identified a set of evolutionarily conserved sites that mapped to regions surrounding the heme prosthetic group. The SCA matrix also predicted a set of energetically coupled residues that are involved in the maintenance of the structural folding of plant peroxidases. The combination of crystallographic data and SCA analysis provides information about the key structural elements that could contribute to explaining the unique stability of RPTP. (C) 2009 Elsevier Inc. All rights reserved.

Reduction capability of soil humic substances from the Rio Negro basin, Brazil, towards Hg(II) studied by a multimethod approach and principal component analysis (PCA)

This paper characterizes humic substances (HS) extracted from soil samples collected in the Rio Negro basin in the state of Amazonas, Brazil, particularly investigating their reduction capabilities towards Hg(II) in order to elucidate potential mercury cycling/volatilization in this environment. For this reason, a multimethod approach was used, consisting of both instrumental methods (elemental analysis, EPR, solid-state NMR, FIA combined with cold-vapor AAS of Hg(0)) and statistical methods such as principal component analysis (PCA) and a central composite factorial planning method. The HS under study were divided into groups, complexing and reducing ones, owing to different distribution of their functionalities. The main functionalities (cor)related with reduction of Hg(II) were phenolic, carboxylic and amide groups, while the groups related with complexation of Hg(II) were ethers, hydroxyls, aldehydes and ketones. The HS extracted from floodable regions of the Rio Negro basin presented a greater capacity to retain (to complex, to adsorb physically and/or chemically) Hg(II), while nonfloodable regions showed a greater capacity to reduce Hg(II), indicating that HS extracted from different types of regions contribute in different ways to the biogeochemical mercury cycle in the basin of the mid-Rio Negro, AM, Brazil. (c) 2007 Published by Elsevier B.V.

Statistical methods for assessing agreement between double readings of clinical measurements

Statistical analysis of data is crucial in cephalometric investigations. There are certainly excellent examples of good statistical practice in the field, but some articles published worldwide have carried out inappropriate analyses. Objective: The purpose of this study was to show that when the double records of each patient are traced on the same occasion, a control chart for differences between readings needs to be drawn, and limits of agreement and coefficients of repeatability must be calculated. Material and methods: Data from a well-known paper in Orthodontics were used for showing common statistical practices in cephalometric investigations and for proposing a new technique of analysis. Results: A scatter plot of the two radiograph readings and the two model readings with the respective regression lines are shown. Also, a control chart for the mean of the differences between radiograph readings was obtained and a coefficient of repeatability was calculated. Conclusions: A standard error assuming that mean differences are zero, which is referred to in Orthodontics and Facial Orthopedics as the Dahlberg error, can be calculated only for estimating precision if accuracy is already proven. When double readings are collected, limits of agreement and coefficients of repeatability must be calculated. A graph with differences of readings should be presented and outliers discussed.

Cytogenetic analysis of A- and B-chromosomes of Prochilodus lineatus (Teleostei, Prochilodontidae) using different restriction enzyme banding and staining methods

Different cytogenetic techniques were used to analyse the chromosomes of Prochilodus lineatus with the main objective of comparing the base composition of A- and B-chromosomes. The results of digestion of chromosomes with 10 different restriction endonucleases (REs), silver staining, CMA(3) staining and C-banding indicated the existence of different classes of highly repetitive DNA in the A-set and also suggested the existence of compositional differences between the chromatin of A- and B-chromosomes. The 5-BrdU incorporation technique showed a late replicating pattern in all B-chromosomes and in some heterochromatic pericentromeric regions of A-chromosomes. The cleavage with RE BamHI produced a band pattern in all chromosomes of P. lineatus which permitted the tentative pairing of homologues in the karyotype of this species. We concluded that the combined use of the above techniques can contribute to the correct identification of chromosomes and the karyotypic analysis in fishes. on the basis of the results, some aspects of chromosome structure and the origin of the B-chromosomes in P. lineatus are discussed.

EthoSeq: A tool for phylogenetic analysis and data mining in behavioral sequences

This article introduces the software program called EthoSeq, which is designed to extract probabilistic behavioral sequences (tree-generated sequences, or TGSs) from observational data and to prepare a TGS-species matrix for phylogenetic analysis. The program uses Graph Theory algorithms to automatically detect behavioral patterns within the observational sessions. It includes filtering tools to adjust the search procedure to user-specified statistical needs. Preliminary analyses of data sets, such as grooming sequences in birds and foraging tactics in spiders, uncover a large number of TGSs which together yield single phylogenetic trees. An example of the use of the program is our analysis of felid grooming sequences, in which we have obtained 1,386 felid grooming TGSs for seven species, resulting in a single phylogeny. These results show that behavior is definitely useful in phylogenetic analysis. EthoSeq simplifies and automates such analyses, uncovers much of the hidden patterns of long behavioral sequences, and prepares this data for further analysis with standard phylogenetic programs. We hope it will encourage many empirical studies on the evolution of behavior.

Bayesian and multivariate methods applied to favorability quantification in recôncavo Basin and Ribeira Belt, Brazil

A methodology to define favorable areas in petroleum and mineral exploration is applied, which consists in weighting the exploratory variables, in order to characterize their importance as exploration guides. The exploration data are spatially integrated in the selected area to establish the association between variables and deposits, and the relationships among distribution, topology, and indicator pattern of all variables. Two methods of statistical analysis were compared. The first one is the Weights of Evidence Modeling, a conditional probability approach (Agterberg, 1989a), and the second one is the Principal Components Analysis (Pan, 1993). In the conditional method, the favorability estimation is based on the probability of deposit and variable joint occurrence, with the weights being defined as natural logarithms of likelihood ratios. In the multivariate analysis, the cells which contain deposits are selected as control cells and the weights are determined by eigendecomposition, being represented by the coefficients of the eigenvector related to the system's largest eigenvalue. The two techniques of weighting and complementary procedures were tested on two case studies: 1. Recôncavo Basin, Northeast Brazil (for Petroleum) and 2. Itaiacoca Formation of Ribeira Belt, Southeast Brazil (for Pb-Zn Mississippi Valley Type deposits). The applied methodology proved to be easy to use and of great assistance to predict the favorability in large areas, particularly in the initial phase of exploration programs. © 1998 International Association for Mathematical Geology.

Growth and yield analysis in different methods of beet (Cv. 'tall top early wonder') plantation

This objective of this study was to evaluate the effects of different methods of red beet seedling production and direct sowing on the development of the plant in field conditions. 'Tall Top Early Wonder' was used as the cultivar. The experimental design was a randomized blocks with 4 replications of each treatment: T1 - seedling produced in trays of 288 cells; T2 - 200 cells; T3 - 128 cells; T4 - 128 cells; T5 - direct sowing and T6 - bare-root seedling. The seedlings produced in trays and bare root seedlings were transplanted in the field 28 DAS with spacing of 0.20 x 0.10m. The mean height of plants, leaf area, leaf dry matter, petiole dry matter, shoot dry matter, root dry matter, relationship between shoot dry matter and root dry matter, absolute growth rate, relative growth rate, net assimilation rate, leaf area ratio and specific leaf area were determined. Storage root fresh matter and mean storage root diameter were determined starting from 77 DAS. Initial growth of the plants were superior for the direct sowing, resulting in smaller RGR and NAR than the other treatments. Independent of the production method, an increase of the cycle of the crop was verified. T6 had larger delay in the initial development. There was no difference for productivity. Method T1, had less expenses with respect to substrate and space in the vegetation home, without reduction in production.

A comparison of mycolic acid analysis for nontuberculous mycobacteria identification by thin-layer chromatography and molecular methods

The development of fast, inexpensive, and reliable tests to identify nontuberculous mycobacteria (NTM) is needed. Studies have indicated that the conventional identification procedures, including biochemical assays, are imprecise. This study evaluated a proposed alternative identification method in which 83 NTM isolates, previously identified by conventional biochemical testing and in-house M. avium IS1245-PCR amplification, were submitted to the following tests: thin-layer chromatography (TLC) of mycolic acids and PCR-restriction enzyme analysis of hsp65 (PRA). High-performance liquid chromatography (HPLC) analysis of mycolic acids and Southern blot analysis for M. avium IS1245 were performed on the strains that evidenced discrepancies on either of the above tests. Sixty-eight out of 83 (82%) isolates were concordantly identified by the presence of IS1245 and PRA and by TLC mycolic acid analysis. Discrepant results were found between the phenotypic and molecular tests in 12/83 (14.4%) isolates. Most of these strains were isolated from non-sterile body sites and were most probably colonizing in the host tissue. While TLC patterns suggested the presence of polymycobacterial infection in 3/83 (3.6%) cultures, this was the case in only one HPLC-tested culture and in none of those tested by PRA. The results of this study indicated that, as a phenotypic identification procedure, TLC mycolic acid determination could be considered a relatively simple and cost-effective method for routine screening of NTM isolates in mycobacteriology laboratory practice with a potential for use in developing countries. Further positive evidence was that this method demonstrated general agreement on MAC and M. simiae identification, including in the mixed cultures that predominated in the isolates of the disseminated infections in the AIDS patients under study. In view of the fact that the same treatment regimen is recommended for infections caused by these two species, TLC mycolic acid analysis may be a useful identification tool wherever molecular methods are unaffordable.

Comparative analysis between three methods of bone estimating age in individuals with down syndrome by mode of the hand and wrist ray

The wrist and hand region has been the most commonly used for estimating age and osseous development due to the great number of ossification centers. The aim was to determine which method, Tanner & Whitehouse's (TW3), Greulich & Pyle's (GP) or Eklof & Ringertz's, more closely relates to the chronological age in subjects with Down syndrome with chronological ages between 61 and 180 months, using wrist and hand radiographs. The sample consisted of 85 radiographs, 52 of males and 33 of females. Eklof & Ringertz's method was computerized (Radiomemory). Greulich & Pyle's atlas was used and compared with the wrist and hand radiographs. For the TW3 method, 13 ossification centers were evaluated; for each one of them, there are seven or eight development stages to which scores are assigned; these scores are then added and the results are transformed into osseous age values. No statistically significant differences were observed between the male and female genders for methods TW3 and GP, contrasting with the observed differences for the Eklof & Ringertz method. Correlation (r2) between osseous and chronological ages was 0.8262 for TW3 and 0.7965 for GP, while for the method of Eklof & Ringertz, it was 0.7656 for females and 0.8353 for males. The author concluded that the osseous age assessment method that better related to the chronological age was the TW3, followed by Greulich & Pyle's and Eklof & Ringertz's.

A study of neolignan compounds with biological activity against Paracoccidioides brasiliensis by using quantum chemical and chemometric methods

Métodos quimiométricos (estatísticos) são empregados para classificar um conjunto de compostos derivados de neolignanas com atividade biológica contra a Paracoccidioides brasiliensis. O método AM1 (Austin Model 1) foi utilizado para calcular um conjunto de descritores moleculares (propriedades) para os compostos em estudo. A seguir, os descritores foram analisados utilizando os seguintes métodos de reconhecimento de padrões: Análise de Componentes Principais (PCA), Análise Hierárquica de Agrupamentos (HCA) e o método de K-vizinhos mais próximos (KNN). Os métodos PCA e HCA mostraram-se bastante eficientes para classificação dos compostos estudados em dois grupos (ativos e inativos). Três descritores moleculares foram responsáveis pela separação entre os compostos ativos e inativos: energia do orbital molecular mais alto ocupado (EHOMO), ordem de ligação entre os átomos C1'-R7 (L14) e ordem de ligação entre os átomos C5'-R6 (L22). Como as variáveis responsáveis pela separação entre compostos ativos e inativos são descritores eletrônicos, conclui-se que efeitos eletrônicos podem desempenhar um importante papel na interação entre receptor biológico e compostos derivados de neolignanas com atividade contra a Paracoccidioides brasiliensis.

The effect of power-ultrasound on the pretreatment of acidified aqueous solutions of banana flower-stalk: Structural, chemical and statistical analysis

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