Evaluation of FAO Penman-Monteith and alternative methods for estimating reference evapotranspiration with missing data in Southern Ontario, Canada

Grass reference evapotranspiration (ETo) is an important agrometeorological parameter for climatological and hydrological studies, as well as for irrigation planning and management. There are several methods to estimate ETo, but their performance in different environments is diverse, since all of them have some empirical background. The FAO Penman-Monteith (FAD PM) method has been considered as a universal standard to estimate ETo for more than a decade. This method considers many parameters related to the evapotranspiration process: net radiation (Rn), air temperature (7), vapor pressure deficit (Delta e), and wind speed (U); and has presented very good results when compared to data from lysimeters Populated with short grass or alfalfa. In some conditions, the use of the FAO PM method is restricted by the lack of input variables. In these cases, when data are missing, the option is to calculate ETo by the FAD PM method using estimated input variables, as recommended by FAD Irrigation and Drainage Paper 56. Based on that, the objective of this study was to evaluate the performance of the FAO PM method to estimate ETo when Rn, Delta e, and U data are missing, in Southern Ontario, Canada. Other alternative methods were also tested for the region: Priestley-Taylor, Hargreaves, and Thornthwaite. Data from 12 locations across Southern Ontario, Canada, were used to compare ETo estimated by the FAD PM method with a complete data set and with missing data. The alternative ETo equations were also tested and calibrated for each location. When relative humidity (RH) and U data were missing, the FAD PM method was still a very good option for estimating ETo for Southern Ontario, with RMSE smaller than 0.53 mm day(-1). For these cases, U data were replaced by the normal values for the region and Delta e was estimated from temperature data. The Priestley-Taylor method was also a good option for estimating ETo when U and Delta e data were missing, mainly when calibrated locally (RMSE = 0.40 mm day(-1)). When Rn was missing, the FAD PM method was not good enough for estimating ETo, with RMSE increasing to 0.79 mm day(-1). When only T data were available, adjusted Hargreaves and modified Thornthwaite methods were better options to estimate ETo than the FAO) PM method, since RMSEs from these methods, respectively 0.79 and 0.83 mm day(-1), were significantly smaller than that obtained by FAO PM (RMSE = 1.12 mm day(-1). (C) 2009 Elsevier B.V. All rights reserved.

Definition of priority areas for forest conservation through the ordered weighted averaging method

The general objective of this study was to evaluate the ordered weighted averaging (OWA) method, integrated to a geographic information systems (GIS), in the definition of priority areas for forest conservation in a Brazilian river basin, aiming at to increase the regional biodiversity. We demonstrated how one could obtain a range of alternatives by applying OWA, including the one obtained by the weighted linear combination method and, also the use of the analytic hierarchy process (AHP) to structure the decision problem and to assign the importance to each criterion. The criteria considered important to this study were: proximity to forest patches; proximity among forest patches with larger core area; proximity to surface water; distance from roads: distance from urban areas; and vulnerability to erosion. OWA requires two sets of criteria weights: the weights of relative criterion importance and the order weights. Thus, Participatory Technique was used to define the criteria set and the criterion importance (based in AHP). In order to obtain the second set of weights we considered the influence of each criterion, as well as the importance of each one, on this decision-making process. The sensitivity analysis indicated coherence among the criterion importance weights, the order weights, and the solution. According to this analysis, only the proximity to surface water criterion is not important to identify priority areas for forest conservation. Finally, we can highlight that the OWA method is flexible, easy to be implemented and, mainly, it facilitates a better understanding of the alternative land-use suitability patterns. (C) 2008 Elsevier B.V. All rights reserved.

Suitability of relative humidity as an estimator of leaf wetness duration

Leaf wetness duration (LWD) is a key parameter in agricultural meteorology since it is related to epidemiology of many important crops, controlling pathogen infection and development rates. Because LWD is not widely measured, several methods have been developed to estimate it from weather data. Among the models used to estimate LWD, those that use physical principles of dew formation and dew and/or rain evaporation have shown good portability and sufficiently accurate results, but their complexity is a disadvantage for operational use. Alternatively, empirical models have been used despite their limitations. The simplest empirical models use only relative humidity data. The objective of this study was to evaluate the performance of three RH-based empirical models to estimate LWD in four regions around the world that have different climate conditions. Hourly LWD, air temperature, and relative humidity data were obtained from Ames, IA (USA), Elora, Ontario (Canada), Florence, Toscany (Italy), and Piracicaba, Sao Paulo State (Brazil). These data were used to evaluate the performance of the following empirical LWD estimation models: constant RH threshold (RH >= 90%); dew point depression (DPD); and extended RH threshold (EXT_RH). Different performance of the models was observed in the four locations. In Ames, Elora and Piracicaba, the RH >= 90% and DPD models underestimated LWD, whereas in Florence these methods overestimated LWD, especially for shorter wet periods. When the EXT_RH model was used, LWD was overestimated for all locations, with a significant increase in the errors. In general, the RH >= 90% model performed best, presenting the highest general fraction of correct estimates (F(C)), between 0.87 and 0.92, and the lowest false alarm ratio (F(AR)), between 0.02 and 0.31. The use of specific thresholds for each location improved accuracy of the RH model substantially, even when independent data were used; MAE ranged from 1.23 to 1.89 h, which is very similar to errors obtained with published physical models for LWD estimation. Based on these results, we concluded that, if calibrated locally, LWD can be estimated with acceptable accuracy by RH above a specific threshold, and that the EXT_RH method was unsuitable for estimating LWD at the locations used in this study. (C) 2007 Elsevier B.V. All rights reserved.

Characterization of fungal soil communities by F-RISA and arbuscular mycorrhizal fungi from Araucaria angustifolia forest soils after replanting and wildfire disturbances

The Fungal Ribosomal Intergenic Spacer Analysis (F-RISA) was used to characterize soil fungal communities from three ecosystems of Araucaria angustifolia from Brazil: a native forest and two replanted forest ecosystems, one of them with a past history of wildfire. The arbuscular mycorrhizal fungi (AMF) infection was evaluated in Araucaria roots of 18-month-old axenic plants previously inoculated with soils collected from those areas in a greenhouse experiment. The principal component analysis of F-RISA profiles showed different soil fungal community between the three studied areas. Sixty three percent of F-RISA fragments amplified in the soil and the substrate samples presented lengths between 500 and 700 bp. The number of Operational Taxonomic Units (OTUs) was 34 for soil and 38 for substrate, however, more fragments were detected in soil (214) than in substrate (163). An in silico F-RISA analysis to compare our data with ITS1-5.8S-ITS2 sequences from NCBI database showed the presence of Ascomycota, Basidiomycota and Glomeromycota among the soil and substrate fungal communities. AMF infection was higher in plants inoculated with soil from the native forest and the replanted forest with wildfire, both presenting similar chemical characteristics but with different disturbance levels. These results indicate that soil chemical composition may influence the soil fungal community structures rather than the anthropogenic or fire disturbances.

Micronutrient dynamics in the forest biomass in succession stages of a seasonal deciduous forest, RS, Brazil

There is a great need of research to assess the behavior of micronutrients in natural forests of southern Brazil. Do to this need, the objective of this work was to study the levels and amounts of micronutrients in forest above ground biomass of the forest, in a comparative way, in two secondary succession stages (SSS) in a Seasonal Deciduous Forest in Rio Grande do Sul, Brazil. The SSS had enjoyed 35 and 55 years of regeneration since the end of agricultural use, respectively for initial secondary forest (ISF) and late secondary forest (LSF). The above-ground biomass was collected and separated into vegetative strata and these in fractions, thereafter chemically analyzed for the levels of B, Fe, Zn, Mn and Cu. Leaf fractions of arboreal, shrubs and herbaceous strata showed the highest levels for most nutrients. Only the levels of iron and manganese were higher in the bark fraction, for both sucession stages. In the LSF, the herbaceous stratum also showed high levels of Fe. The average levels of micronutrients showed differences between the two sucession stages only in relation to Fe and Mn, with higher levels in LSF biomass. The amount of nutrients stored was always higher in LSF, because of the largest biomass and the higher levels of Fe and Mn in the biomass of this SSS. The quantitative order of nutrient storage in biomass was Fe> Mn> Zn> B> Cu.

Within-stand and seasonal variations of specific leaf area in a clonal Eucalyptus plantation in the Republic of Congo

Specific leaf area (SLA; m(leaf)(2) kg(leaf)(-1)) is a key ecophysiological parameter influencing leaf physiology, photosynthesis, and whole plant carbon gain. Both individual tree-based models and other forest process-based models are generally highly sensitive to this parameter, but information on its temporal or within-stand variability is still scarce. In a 2-4-year-old Eucalyptus plantation in Congo, prone to seasonal drought, the within-stand and seasonal variability in SLA were investigated by means of destructive sampling carried out at 2-month intervals, over a 2-year period. Within-crown vertical gradients of SLA were small. Highly significant relationships were found between tree-average SLA (SLA(t)) and tree size (tree height, H(t), or diameter at breast height, DBH): SLA(t) ranged from about 9 m(2) kg(-1) for dominant trees to about 14-15 m(2) kg(-1) for the smallest trees. The decrease in SLA(t) with increasing tree size was accurately predicted from DBH using power functions. Stand-average SLA varied by about 20% during the year, with lowest values at the end of the 5-month dry season, and highest values about 2-3 months after the onset of the wet season. Variability in leaf water status according to tree size and season is discussed as a possible determinant of both the within-stand and seasonal variations in SM. (C) 2009 Elsevier B.V. All rights reserved.

Natural regeneration of Atlantic Forest species in the understory of Eucalyptus saligna Smith. in a former forest production unit at the Parque das Neblinas, Bertioga, SP

This study aimed at characterizing the potential for natural regeneration of native vegetation in the under-story of an earlier Eucalyptus saligna Smith production stand. The study was carried out at the Parque das Neblinas, Bertioga municipality, SP, in a 45 ha third rotation stand; which had been abandoned 15 years ago for natural regeneration to occur. The sampling was done in 24 plots of 20 x 40 m. The sampled area was of 19,200 m(2), with inventory made of 100% of the eucalyptus trees. All regeneration trees with a height >= 1.30 m and DBH >= 5.0 cm were measured, as well as adult individuals with DBH >= 5.0 cm; surveyed in two size classes. 1,417 individuals of E. saligna were measured, with a density of 738,02 individuals/ha and a basal area of 22.69 m(2)/ha. Among 2,763 natural regeneration individuals, 111 species belonged to 66 genera and 34 botanical families. The species represented 43.7% of the tree richness of neighboring native forest fragments. The total estimated density and the basal area were respectively 1,052.6 individuals/ha and 6.4 m(2)/ha of autochthonous trees with DBH >= 5.0 cm (Class 1); while for regeneration there were 3,864.58 individuals/ha, and 2.76 m(2)/ha of individuals with a height >= 1.30 m and DBH <5.0 cm (Class 2). Shannon diversity (H`) was 2.83 and 3.68, respectively, for Classes 1 and 2, and the corrected species richness for a 1000-individual sample (R(1000)) were 75.6 and 87.29 (Fisher`s a index) for the same classes. The majority of the species (34.84%) was typical from the understory of wet tropical forest and had zoochoric fruit dispersal (67.57%). The results indicate that, under these conditions, a eucalyptus forest is able to provide adequate regeneration niches for native vegetation, and may represent a sink habitat for local populations.

Methods for detection of anatoxin-a(s) by liquid chromatography coupled to electrospray ionization-tandem mass spectrometry

Anatoxin-a(s) is a potent irreversible inhibitor of the enzyme acetylcholinesterase with a unique N-hydroxyguanidine methylphosphate ester chemical structure. Determination of this toxin in environmental samples is hampered by the lack of specific methods for its detection. Using the toxic strain of Anabaena lemmermani PH-160 B as positive control, the fragmentation characteristics of anatoxin-a(s) under collision-induced dissociation conditions have been investigated and new LC-MS/MS methods proposed. Recommended ion transitions for correct detection of this toxin are 253 > 58, 253 > 159, 235 > 98 and 235 > 96. Chromatographic separation is better achieved under HILIC conditions employing a ZIC-HILIC column. This method was used to confirm for the first time the production of anatoxin-a(s) by strains of Anabaena oumiana ITEP-025 and ITEP-026. Considering no standard solutions are commercially available, our results will be of significant use for the correct identification of this toxin by LC-MS/MS. (C) 2009 Elsevier Ltd. All rights reserved.

In vitro antioxidant activity and in vivo efficacy of topical formulations containing vitamin C and its derivatives studied by non-invasive methods

Background/purpose: Vitamins C and its derivatives, mainly due to their antioxidant properties, are being used in cosmetic products to protect and to reduce the signs of ageing. However, there are no studies comparing the effects of vitamin C [ascorbic acid (AA)] and its derivatives, magnesium ascorbyl phosphate (MAP) and ascorbyl tetra-isopalmitate (ATIP), when vehiculated in topical formulations, mainly using objective measurements, which are an important tool in clinical efficacy studies. Thus, the objective of this study was to determine the in vitro antioxidant activity of AA and its derivatives, MAP and ATIP, as well as their in vivo efficacy on human skin, when vehiculated in topical formulations. Methods: The study of antioxidant activity in vitro was performed with an aqueous and a lipid system. The in vivo methodology consisted of the application of these formulations on human volunteers` forearm skin and the analysis of the skin conditions after 4-week period daily applications in terms of transepidermal water loss (TEWL), stratum corneum moisture content and viscoelasticity using a Tewameter (R), Corneometer (R) and Cutometer (R), respectively. Results: In vitro experiments demonstrated that in an aqueous system, AA had the best antioxidant potential, and MAP was more effective than ATIP, whereas in the lipid system ATIP was more effective than MAP. In in vivo studies, all formulations enhanced stratum corneum moisture content after a 4-week period daily applications when compared with baseline values; however, only the formulation containing AA caused alterations in TEWL values. The formulations containing MAP caused alterations in the viscoelastic-to-elastic ratio, which suggested its action in the deeper layers of the skin. Conclusion: AA and its derivates presented an in vitro antioxidant activity but AA had the best antioxidant effect. In in vivo efficacy studies, only the formulation containing AA caused alterations in TEWL values and the formulation containing MAP caused alterations in the viscoelastic-to-elastic ratio. This way, vitamin C derivatives did not present the same effects of AA on human skin; however, MAP showed other significant effect-improving skin hydration, which is very important for the normal cutaneous metabolism and also to prevent skin alterations and early ageing.

Quercetin inhibits UV irradiation-induced inflammatory cytokine production in primary human keratinocytes by suppressing NF-kappa B pathway

Background: Topical flavonoids, such as quercetin, have been shown to reduce ultraviolet (UV) irradiation-mediated skin damage. However, the mechanisms and signaling pathways involved in this protective effect are not clear. UV irradiation leads to activation of two major signaling pathways, namely nuclear factor kappa B (NF-kappa B) and activator protein-1 (AP-1) pathways. Activation of NF-kappa B pathway by UV irradiation stimulates inflammatory cytokine expression, whereas activation of AP-1 pathway by UV irradiation promotes matrix metalloproteinase (MMP) production. Both pathways contribute to UV irradiation-induced skin damage, such as photoaging and skin tumor formation. Objective: To elucidate the underlying mechanism, we examined the effect of quercetin on UV irradiation induced activation of NF-kappa B and AP-1 pathways. Methods: Primary human keratinocytes, the major skin cell type subjected to physiological solar UV irradiation, were used to study the effects of quercetin on UV irradiation-induced signal transduction pathways. Results: Quercetin decreased UV irradiation-induced NF-kappa B DNA-binding by 80%. Consequently, quercetin suppressed UV irradiation-induced expression of inflammatory cytokines IL-1 beta (similar to 60%), IL-6 (similar to 80%), IL-8 (similar to 76%) and TNF-alpha (similar to 69%). In contrast, quercetin had no effect on UV irradiation activation of three MAP kinases, ERK, JNK, or p38. Accordingly, induction of AP-1 target genes such as MMP-1 and MMP-3 by UV irradiation was not suppressed by quercetin. Conclusion: Our data indicate that the ability of quercetin to block UV irradiation-induced skin inflammation is mediated, at least in part, by its inhibitory effect on NF-kappa B activation and inflammatory cytokine production. (C) 2011 Japanese Society for Investigative Dermatology. Published by Elsevier Ireland Ltd. All rights reserved.

Determination of mitragynine in rat plasma by LC-MS/MS: Application to pharmacokinetics

This study used for the first time LC-MS/MS for the analysis of mitragynine (MIT), a mu-opioid agonist with antinociceptive and antitussive properties, in rat plasma. Mitragynine and the internal standard (amitriptyline) were extracted from plasma with hexane-isoamyl alcohol and resolved on a Lichrospher (R) RP-SelectB column (9.80 and 12.90 min, respectively). The quantification limit was 0.2 ng/mL within a linear range of 0.2-1000 ng/mL The method was applied to quantify mitragynine in plasma samples of rats (n = 8 per sampling time) treated with a single oral dose of 20 mg/kg. The following pharmacokinetic parameters were obtained (mean): maximum plasma concentration: 424 ng/mL; time to reach maximum plasma concentration: 1.26 h; elimination half-life: 3.85 h, apparent total clearance: 6.35 L/h/kg, and apparent volume of distribution: 37.90 L/kg. (C) 2009 Elsevier B.V. All rights reserved.

A chemometric study on the analgesic activity of cannabinoid compounds using SDA, KNN and SIMCA methods

The supervised pattern recognition methods K-Nearest Neighbors (KNN), stepwise discriminant analysis (SDA), and soft independent modelling of class analogy (SIMCA) were employed in this work with the aim to investigate the relationship between the molecular structure of 27 cannabinoid compounds and their analgesic activity. Previous analyses using two unsupervised pattern recognition methods (PCA-principal component analysis and HCA-hierarchical cluster analysis) were performed and five descriptors were selected as the most relevants for the analgesic activity of the compounds studied: R (3) (charge density on substituent at position C(3)), Q (1) (charge on atom C(1)), A (surface area), log P (logarithm of the partition coefficient) and MR (molecular refractivity). The supervised pattern recognition methods (SDA, KNN, and SIMCA) were employed in order to construct a reliable model that can be able to predict the analgesic activity of new cannabinoid compounds and to validate our previous study. The results obtained using the SDA, KNN, and SIMCA methods agree perfectly with our previous model. Comparing the SDA, KNN, and SIMCA results with the PCA and HCA ones we could notice that all multivariate statistical methods classified the cannabinoid compounds studied in three groups exactly in the same way: active, moderately active, and inactive.

Sampling phylogenetic tree space with the generalized Gibbs sampler

The generalized Gibbs sampler (GGS) is a recently developed Markov chain Monte Carlo (MCMC) technique that enables Gibbs-like sampling of state spaces that lack a convenient representation in terms of a fixed coordinate system. This paper describes a new sampler, called the tree sampler, which uses the GGS to sample from a state space consisting of phylogenetic trees. The tree sampler is useful for a wide range of phylogenetic applications, including Bayesian, maximum likelihood, and maximum parsimony methods. A fast new algorithm to search for a maximum parsimony phylogeny is presented, using the tree sampler in the context of simulated annealing. The mathematics underlying the algorithm is explained and its time complexity is analyzed. The method is tested on two large data sets consisting of 123 sequences and 500 sequences, respectively. The new algorithm is shown to compare very favorably in terms of speed and accuracy to the program DNAPARS from the PHYLIP package.