Forecasting fuel ethanol consumption in Brazil by time series models: 2006-2012

This article analysed scenarios for Brazilian consumption of ethanol for the period 2006 to 2012. The results show that if the country`s GDP sustains a 4.6% a year growth, domestic consumption of fuel ethanol could increase to 25.16 billion liters in this period, which is a volume relatively close to the forecasted gasoline consumption of 31 billion liters. At a lower GDP growth of 1.22% a year, gasoline consumption would be reduced and domestic ethanol consumption in Brazil would be no higher than 18.32 billion liters. Contrary to the current situation, forecasts indicated that hydrated ethanol consumption could become much higher than anhydrous consumption in Brazil. The former is being consumed in cars moved exclusively by ethanol and flex-fuel cars, successfully introduced in the country at 2003. Flex cars allow Brazilian consumers to choose between gasoline and hydrated ethanol and immediately switch to whichever fuel presents the most favourable relative price.

Pricing farm-level agricultural insurance: a Bayesian approach

This paper applies Hierarchical Bayesian Models to price farm-level yield insurance contracts. This methodology considers the temporal effect, the spatial dependence and spatio-temporal models. One of the major advantages of this framework is that an estimate of the premium rate is obtained directly from the posterior distribution. These methods were applied to a farm-level data set of soybean in the State of the Parana (Brazil), for the period between 1994 and 2003. The model selection was based on a posterior predictive criterion. This study improves considerably the estimation of the fair premium rates considering the small number of observations.

Spatio-temporal modeling of agricultural yield data with an application to pricing crop insurance contracts

This article presents a statistical model of agricultural yield data based on a set of hierarchical Bayesian models that allows joint modeling of temporal and spatial autocorrelation. This method captures a comprehensive range of the various uncertainties involved in predicting crop insurance premium rates as opposed to the more traditional ad hoc, two-stage methods that are typically based on independent estimation and prediction. A panel data set of county-average yield data was analyzed for 290 counties in the State of Parana (Brazil) for the period of 1990 through 2002. Posterior predictive criteria are used to evaluate different model specifications. This article provides substantial improvements in the statistical and actuarial methods often applied to the calculation of insurance premium rates. These improvements are especially relevant to situations where data are limited.

Spatial pattern of black spot incidence within citrus trees related to disease severity and pathogen dispersal

Guignardia citricarpa, the causal agent of citrus black spot, forms airborne ascospores on decomposing citrus leaves and water-spread conidia on fruits, leaves and twigs. The spatial pattern of diseased fruit in citrus tree canopies was used to assess the importance of ascospores and conidia in citrus black spot epidemics in Sao Paulo State, Brazil. The aggregation of diseased fruit in the citrus tree canopy was quantified by the binomial dispersion index (D) and the binary form of Taylor`s Power Law for 303 trees in six groves. D was significantly greater than 1 in 251 trees. The intercept of the regression line of Taylor`s Power Law was significantly greater than 0 and the slope was not different from 1, implying that diseased fruit was aggregated in the canopy independent of disease incidence. Disease incidence (p) and severity (S) were assessed in 2875 citrus trees. The incidence-severity relationship was described (R-2 = 88.7%) by the model ln(S) = ln(a) + bCLL(p) where CLL = complementary log-log transformation. The high severity at low incidence observed in many cases is also indicative of low distance spread of G. citricarpa spores. For the same level of disease incidence, some trees had most of the diseased fruit with many lesions and high disease severity, whereas other trees had most of the fruit with few lesions and low disease severity. Aggregation of diseased fruit in the trees suggests that splash-dispersed conidia have an important role in increasing the disease in citrus trees in Brazil.

Does using stepwise variable selection to build sequential path analysis models make sense?

Causal inference methods - mainly path analysis and structural equation modeling - offer plant physiologists information about cause-and-effect relationships among plant traits. Recently, an unusual approach to causal inference through stepwise variable selection has been proposed and used in various works on plant physiology. The approach should not be considered correct from a biological point of view. Here, it is explained why stepwise variable selection should not be used for causal inference, and shown what strange conclusions can be drawn based upon the former analysis when one aims to interpret cause-and-effect relationships among plant traits.

Spatial Modeling of a Soil Fertility Index using Visible-Near-Infrared Spectra and Terrain Attributes

Our objective was to develop a methodology to predict soil fertility using visible near-infrared (vis-NIR) diffuse reflectance spectra and terrain attributes derived from a digital elevation model (DEM). Specifically, our aims were to: (i) assemble a minimum data set to develop a soil fertility index for sugarcane (Sarcharum officinarum L.) (SFI-SC) for biofuel production in tropical soils; (ii) construct a model to predict the SFI-SC using soil vis-NIR spectra and terrain attributes; and (iii) produce a soil fertility map for our study area and assess it by comparing it with a green vegetation index (GVI). The study area was 185 ha located in sao Paulo State, Brazil. In total, 184 soil samples were collected and analyzed for a range of soil chemical and physical properties. Their vis-NIR spectra were collected from 400 to 2500 nm. The Shuttle Radar Topographic Mission 3-arcsec (90-m resolution) DEM of the area was used to derive 17 terrain attributes. A minimum data set of soil properties was selected to develop the SFI-SC. The SFI-SC consisted of three classes: Class 1, the highly fertile soils; Class 2, the fertile soils; and Class 3, the least fertile soils. It was derived heuristically with conditionals and using expert knowledge. The index was modeled with the spectra and terrain data using cross-validated decision trees. The cross-validation of the model correctly predicted Class 1 in 75% of cases, Class 2 in 61%, and Class 3 in 65%. A fertility map was derived for the study area and compared with a map of the GVI. Our approach offers a methodology that incorporates expert knowledge to derive the SFI-SC and uses a versatile spectro-spatial methodology that may be implemented for rapid and accurate determination of soil fertility and better exploration of areas suitable for production.

Using Imaging Spectroscopy to study soil properties

Imaging Spectroscopy (IS) is a promising tool for studying soil properties in large spatial domains. Going from point to image spectrometry is not only a journey from micro to macro scales, but also a long stage where problems such as dealing with data having a low signal-to-noise level, contamination of the atmosphere, large data sets, the BRDF effect and more are often encountered. In this paper we provide an up-to-date overview of some of the case studies that have used IS technology for soil science applications. Besides a brief discussion on the advantages and disadvantages of IS for studying soils, the following cases are comprehensively discussed: soil degradation (salinity, erosion, and deposition), soil mapping and classification, soil genesis and formation, soil contamination, soil water content, and soil swelling. We review these case studies and suggest that the 15 data be provided to the end-users as real reflectance and not as raw data and with better signal-to-noise ratios than presently exist. This is because converting the raw data into reflectance is a complicated stage that requires experience, knowledge, and specific infrastructures not available to many users, whereas quantitative spectral models require good quality data. These limitations serve as a barrier that impedes potential end-users, inhibiting researchers from trying this technique for their needs. The paper ends with a general call to the soil science audience to extend the utilization of the IS technique, and it provides some ideas on how to propel this technology forward to enable its widespread adoption in order to achieve a breakthrough in the field of soil science and remote sensing. (C) 2009 Elsevier Inc. All rights reserved.

Soil CO(2) emission estimated by different interpolation techniques

Soil CO(2) emissions are highly variable, both spatially and across time, with significant changes even during a one-day period. The objective of this study was to compare predictions of the diurnal soil CO(2) emissions in an agricultural field when estimated by ordinary kriging and sequential Gaussian simulation. The dataset consisted of 64 measurements taken in the morning and in the afternoon on bare soil in southern Brazil. The mean soil CO(2) emissions were significantly different between the morning (4.54 mu mol m(-2) s(-1)) and afternoon (6.24 mu mol m(-2) s(-1)) measurements. However, the spatial variability structures were similar, as the models were spherical and had close range values of 40.1 and 40.0 m for the morning and afternoon semivariograms. In both periods, the sequential Gaussian simulation maps were more efficient for the estimations of emission than ordinary kriging. We believe that sequential Gaussian simulation can improve estimations of soil CO(2) emissions in the field, as this property is usually highly non-Gaussian distributed.

Effects of dietary supplementation of rumen-protected conjugated linoleic acid to grazing cows in early lactation

Conjugated linoleic acids (CLA) are potent anticarcinogens in animal and in vitro models as well as inhibitors of fatty acid synthesis in mammary gland, liver, and adipose tissue. Our objective was to evaluate long-term CLA supplementation of lactating dairy cows in tropical pasture on milk production and composition and residual effects posttreatment. Thirty crossbred cows grazing stargrass (Cynodon nlemfuensis Vanderyst var. nlemfuensis) were blocked by parity and received 150 g/d of a dietary fat supplement of either Ca-salts of palm oil fatty acids (control) or a mixture of Ca-salts of CLA (CLA treatment). Supplements of fatty acids were mixed with 4 kg/d of concentrate. Grazing plus supplements were estimated to provide 115% of the estimated metabolizable protein requirements from 28 to 84 d in milk (treatment period). The CLA supplement provided 15 g/d of cis-9, trans-11 and 22 g of cis-10, trans-12. Residual effects were evaluated from 85 to 112 d in milk (residual period) when cows were fed an 18% crude protein concentrate without added fat. The CLA treatment increased milk production but reduced milk fat concentration from 2.90 to 2.14% and fat production from 437 to 348 g/d. Milk protein concentration increased by 11.5% (2.79 to 3.11%) and production by 19% (422 to 504 g/d) in the cows fed CLA. The CLA treatment decreased milk energy concentration and increased milk volume, resulting in unchanged energy output. Milk production and protein concentration and production were also greater during the residual period for the CLA-treated cows. The CLA treatment reduced production of fatty acids (FA) of all chain lengths, but the larger effect was on short-chain FA, causing a shift toward a greater content of longer chain FA. The CLA treatment increased total milk CLA content by 30% and content of the trans-10, cis-12 CLA isomer by 88%. The CLA treatment tended to decrease the number of days open, suggesting a possible effect on reproduction. Under tropical grazing conditions, in a nutritionally challenging environment, CLA-treated cows decreased milk fat content and secreted the same amount of milk energy by increasing milk volume and milk protein production.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Random regression models to estimate test-day milk yield genetic parameters Holstein cows in Southeastern Brazil

A total of 152,145 weekly test-day milk yield records from 7317 first lactations of Holstein cows distributed in 93 herds in southeastern Brazil were analyzed. Test-day milk yields were classified into 44 weekly classes of DIM. The contemporary groups were defined as herd-year-week of test-day. The model included direct additive genetic, permanent environmental and residual effects as random and fixed effects of contemporary group and age of cow at calving as covariable, linear and quadratic effects. Mean trends were modeled by a cubic regression on orthogonal polynomials of DIM. Additive genetic and permanent environmental random effects were estimated by random regression on orthogonal Legendre polynomials. Residual variances were modeled using third to seventh-order variance functions or a step function with 1, 6,13,17 and 44 variance classes. Results from Akaike`s and Schwarz`s Bayesian information criterion suggested that a model considering a 7th-order Legendre polynomial for additive effect, a 12th-order polynomial for permanent environment effect and a step function with 6 classes for residual variances, fitted best. However, a parsimonious model, with a 6th-order Legendre polynomial for additive effects and a 7th-order polynomial for permanent environmental effects, yielded very similar genetic parameter estimates. (C) 2008 Elsevier B.V. All rights reserved.

Spatial patterns of soil attributes and components in a mangrove system in Southeast Brazil (Sao Paulo)

Purpose Among environmental factors governing innumerous processes that are active in estuarine environments, those of edaphic character have received special attention in recent studies. With the objectives of determining the spatial patterns of soil attributes and components across different mangrove forest landscapes and obtaining additional information on the cause-effect relationships between these variables and position within the estuary, we analyzed several soil attributes in 31 mangrove soil profiles from the state of So Paulo (Guaruja, Brazil). Materials and methods Soil samples were collected at low tide along two transects within the CrumahA(0) mangrove forest. Samples were analyzed to determine pH, Eh, salinity, and the percentages of sand, silt, clay, total organic carbon (TOC), and total S. Mineralogy of the clay fraction (< 2 mm) was also studied by X-ray diffraction analysis, and partitioning of solid-phase Fe was performed by sequential extraction. Results and discussion The results obtained indicate important differences in soil composition at different depths and landscape positions, causing variations in physicochemical parameters, clay mineralogy, TOC contents, and iron geochemistry. The results also indicate that physicochemical conditions may vary in terms of different local microtopographies. Soil salinity was determined by relative position in relation to flood tide and transition areas with highlands. The proportions of TOC and total S are conditioned by the sedimentation of organic matter derived from vegetation and by the prevailing redox conditions, which clearly favored intense sulfate reduction in the soils (similar to 80% of the total Fe is Fe-pyrite). Particle-size distribution is conditioned by erosive/deposition processes (present and past) and probably by the positioning of ancient and reworked sandy ridges. The existing physicochemical conditions appear to contribute to the synthesis (smectite) and transformation (kaolinite) of clay minerals. Conclusions The results demonstrate that the position of soils in the estuary greatly affects soil attributes. Differences occur even at small scales (meters), indicating that both edaphic (soil classification, soil mineralogy, and soil genesis) and environmental (contamination and carbon stock) studies should take such variability into account.

Large-scale ecological restoration of high-diversity tropical forests in SE Brazil

The complex interactions among endangered ecosystems, landowners` interests, and different models of land tenure and use, constitute an important series of challenges for those seeking to maintain and restore biodiversity and augment the flow of ecosystem services. Over the past 10 years, we have developed a data-based approach to address these challenges and to achieve medium and large-scale ecological restoration of riparian areas on private lands in the state of Sao Paulo, southeastern Brazil. Given varying motivations for ecological restoration, the location of riparian areas within landholdings, environmental zoning of different riparian areas, and best-practice restoration methods were developed for each situation. A total of 32 ongoing projects, covering 527,982 ha, were evaluated in large sugarcane farms and small mixed farms, and six different restoration techniques have been developed to help upscale the effort. Small mixed farms had higher portions of land requiring protection as riparian areas (13.3%), and lower forest cover of riparian areas (18.3%), than large sugarcane farms (10.0% and 36.9%, respectively for riparian areas and forest cover values). In both types of farms, forest fragments required some degree of restoration. Historical anthropogenic degradation has compromised forest ecosystem structure and functioning, despite their high-diversity of native tree and shrub species. Notably, land use patterns in riparian areas differed markedly. Large sugarcane farms had higher portions of riparian areas occupied by highly mechanized agriculture, abandoned fields, and anthropogenic wet fields created by siltation in water courses. In contrast, in small mixed crop farms, low or non-mechanized agriculture and pasturelands were predominant. Despite these differences, plantations of native tree species covering the entire area was by far the main restoration method needed both by large sugarcane farms (76.0%) and small mixed farms (92.4%), in view of the low resilience of target sites, reduced forest cover, and high fragmentation, all of which limit the potential for autogenic restoration. We propose that plantations should be carried out with a high-diversity of native species in order to create biologically viable restored forests, and to assist long-term biodiversity persistence at the landscape scale. Finally, we propose strategies to integrate the political, socio-economic and methodological aspects needed to upscale restoration efforts in tropical forest regions throughout Latin America and elsewhere. (C) 2010 Elsevier BA/. All rights reserved.

Determination of Hg and diet identification in otter (Lontra longicaudis) feces

An analytical procedure for the determination of Hg in otter (Lontra longicaudis) feces was developed, to separate fish scales for the identification of the animal diet. Samples were washed with ultra-pure water and the suspension was sampled and transferred for digestion. The solubilization was performed with nitric-perchloric acid mixture, and detection carried out by the atomic fluorescence spectrometry (AFS). The quality of the analytical procedure was assessed by analyzing in-house standard solutions and certified reference materials. Total Hg concentrations were in the range of 7.6-156 ng g(-1) (July 2004), 25.6-277 ng g(-1) (January 2005) and 14.6-744 ng g(-1) (May 2005) that is approximately the same order of magnitude for all samples collected in two reservoirs at the Tiete River, Brazil. Although Hg concentrations varied with sampling periods and diet, high levels were correlated to the percentage of carnivorous fish scales present in the otter feces. (c) 2007 Elsevier Ltd. All rights reserved.

Spatial heterogeneity of leaf wetness duration in apple trees and its influence on performance of a warning system for sooty blotch and flyspeck

To determine the effect of sensor placement on the performance of a disease-warning system for sooty blotch and flyspeck (SBFS), we measured leaf wetness duration (LWD) at 12 canopy positions in apple trees, then simulated operation of the disease-warning system using LWD measurements from different parts of the canopy. LWD sensors were placed in four trees within one Iowa orchard during two growing seasons, and in one tree in each of four orchards during a single growing season. The LWD measurements revealed substantial heterogeneity among sensor locations. In all data sets, the upper, eastern portion of the canopy had the longest mean daily LWD, and was the first site to form dew and the last to dry. The lower, western portion of the canopy averaged about 3 It less LWD per day than the top of the canopy, and was the last zone where dew formed and the first to dry off. On about 25% of nights when dew occurred in the top of the canopy, no dew formed in the lower, western canopy. Intracanopy variability of LWD was more pronounced when dew was the sole source of wetness than on days when rainfall occurred. Daily LWD in the upper, eastern portion of the canopy was slightly less than reference measurements made at a 0.7-m height over turfgrass located near the orchard. When LWD measurements from several canopy positions were input to the SBFS warning system, timing of occurrence of a fungicide-spray threshold varied by as much as 30 days among canopy positions. Under Iowa conditions, placement of an LWD sensor at an unobstructed site over turfgrass was a fairly accurate surrogate for the wettest part of the canopy. Therefore, such an extra-canopy LWD sensor might be substituted for a within-canopy sensor to enhance operational reliability of the SBFS warning system.