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.

Parametric and nonparametric statistical modelling of crop yield: implications for pricing crop insurance contracts

This article considers alternative methods to calculate the fair premium rate of crop insurance contracts based on county yields. The premium rate was calculated using parametric and nonparametric approaches to estimate the conditional agricultural yield density. These methods were applied to a data set of county yield provided by the Statistical and Geography Brazilian Institute (IBGE), for the period of 1990 through 2002, for soybean, corn and wheat, in the State of Paran. In this article, we propose methodological alternatives to pricing crop insurance contracts resulting in more accurate premium rates in a situation of limited data.

Yield loss caused by huanglongbing in different sweet orange cultivars in So Paulo, Brazil

Citrus huanglongbing (HLB) reduces an affected orchard`s economic life. This work aimed to characterize yield loss due to HLB for different sweet orange cultivars and determine the relationship between disease severity and yield. Disease severity and yield were assessed on 949 individual trees distributed in 11 different blocks from sweet orange cultivars Hamlin, Westin, Pera and Valencia. In each block, plants showing a range of HLB severity levels and asymptomatic plants were selected. Total yield (weight of harvested fruit), mean weight of asymptomatic and symptomatic fruit, relative yield (symptomatic tree yield/mean yield of asymptomatic trees from the same block) and relative number of fruits (fruit number from symptomatic tree/mean number of fruits from asymptomatic trees from the same block) were determined. The weight of symptomatic fruit was lower than the weight of asymptomatic fruit, but the weights of asymptomatic and symptomatic fruit were not correlated with disease severity, indicating that the effects of HLB were restricted to symptomatic branches. The relationship of the relative yield with HLB severity can be satisfactorily described by a negative exponential model. The rates of yield decrease as a function of disease severity were similar for all assessed cultivars. A relative yield (up to 19%) was observed even for trees where disease severity was 100%. The strong linear relationship between relative number of fruits per tree and the relative yield per tree suggested that the yield reduction was due primarily to early fruit drop or lack of fruit set on affected branches.

Role of healthy nursery plants in orange yield during eight years of Citrus Variegated Chlorosis epidemics

The production and commercialization of citrus seedlings inspected and produced in protected screen-houses has become mandatory in Sao Paulo State, Brazil since January 2003. This law was intended to avoid the dispersion of Citrus Variegated Chlorosis (CVC), disease caused by Xylella fastidiosa. Our objective was to compare the yield over 8 years of `Natal` sweet orange trees grafted onto Rangpur lime obtained from healthy nursery plants and from plants artificially inoculated with X. fastidiosa. Yield was evaluated in an orchard planted in February 1999 with two treatments: (i) trees from healthy nursery plant, and (ii) trees from plants artificially inoculated with X. fastidiosa. The mean yield was 21% higher in trees from healthy nursery plants, as compared to trees from inoculated nursery plants. This difference represents a gain of approximately 203 boxes of 40.8 kg each, considering a planting density of 550 plants per hectare. (C) 2011 Elsevier B.V. All rights reserved.

Identification of microsatellites markers associated with yield components and quality parameters in sugarcane

Marker assisted selection depends on the identification of tightly linked association between marker and the trait of interest. In the present work, functional (EST-SSRs) and genomic (gSSRs) microsatellite markers were used to detect putative QTLs for sugarcane yield components (stalk number, diameter and height) and as well as for quality parameters (Brix, Pol and fibre) in plant cane. The mapping population (200 individuals) was derived from a bi-parental cross (IACSP95-3018 x IACSP93-3046) from the IAC Sugarcane Breeding Program. As the map is under construction, single marker trait association analysis based on the likelihood ratio test was undertaken to detect the QTLs. Of the 215 single dose markers evaluated (1:1 and 3:1), 90 (42%) were associated with putative QTLs involving 43 microsatellite primers (18 gSSRs and 25 EST-SSRs). For the yield components, 41 marker/trait associations were found: 20 for height, 6 for diameter and 15 for stalk number. An EST-SSRs marker with homology to non-phototropic hypocotyls 4 (NPH4) protein was associated with a putative QTL with positive effect for diameter as also with a negative effect for stalk number. In relation to the quality parameters, 18 marker trait associations were found for Brix, 19 for Pol, and 12 for fibre. For fibre, 58% of the QTLs detected showed a negative effect on this trait. Some makers associated with QTLs with a negative effect for fibre showed a positive effect for Pol, reflecting the negative correlation generally observed between these traits.

Analysis of genomic and functional RFLP derived markers associated with sucrose content, fiber and yield QTLs in a sugarcane (Saccharum spp.) commercial cross

Expressed sequence tags derived markers have a great potential to be used in functional map construction and QTL tagging. In the present work, sugarcane genomic probes and expressed sequence tags having homology to genes, mostly involved in carbohydrate metabolism were used in RFLP assays to identify putative QTLs as well as their epistatic interactions for fiber content, cane yield, pol and tones of sugar per hectare, at two crop cycles in a progeny derived from a bi-parental cross of sugarcane elite materials. A hundred and twenty marker trait associations were found, of which 26 at both crop cycle and 32 only at first ratoon cane. A sucrose synthase derived marker was associated with a putative QTL having a high negative effect on cane yield and also with a QTL having a positive effect on Pol at both crop cycles. Fifty digenic epistatic marker interactions were identified for the four traits evaluated. Of these, only two were observed at both crop cycles.

Tree performance and fruit yield and quality of `Okitsu` Satsuma mandarin grafted on 12 rootstocks

The citriculture in Brazil, as well as in other important regions in the world, is based on very few mandarin cultivars. This fact leads to a short harvest period and higher prices for off-season fruit. The `Okitsu` Satsuma (Citrus unshiu Marc.) is among the earliest ripening mandarin cultivars and it is considered to be tolerant to, citrus canker (Xanthomonas citri subsp. citri Schaad et al.) and to citrus variegated chlorosis (Xylella fastidiosa Wells et al.). Despite having regular fruit quality under hot climate conditions, the early fruit maturation and absence of seeds of `Okitsu` fruits are well suited for the local market in the summer(December through March), when the availability of citrus fruits for fresh consumption is limited. Yet, only a few studies have been conducted in Brazil on rootstocks for `Okitsu`. Consequently, a field trial was carried out in Bebeclouro, Sao Paulo State, to evaluate the horticultural performance of `Okitsu` Satsuma mandarin budded onto 12 rootstocks: the citrandarin `Changsha` mandarin (Citrus reticulata Blanco) x Poncirus trifoliata `English Small`: the hybrid Rangpur lime (Citrus limonia Osbeck) x `Swingle` citrumelo (P. trifoliata (L.) Raf. x Citrus paradisi Macfad.); the trifoliates (P. trifoliata (L) Raf)`Rubidoux`,`FCAV` and `Flying Dragon`(P. trifoliata var. monstrosa); the mandarins `Sun Chu Sha Kat`(C. reticulata Blanco) and `Sunki`(Citrus sunki (Hayata) Hort. ex. Tanaka); the Rangpur limes (C. limonia Osbeck) `Cravo Limeira` and `Cravo FCAV`;`Carrizo` citrange (Citrus sinensis x P. trifoliata), `Swingle` citrumelo (P. trifoliata x C. paradisi), and `Orlando` tangelo (C. paradisi x Citrus tangerina cv. `Dancy`). The experimental grove was planted in 2001, using a 6 m x 3 m spacing, in a randomized block design. No supplementary irrigation was applied. Fruit yield, canopy volume, and fruit quality were assessed for each rootstock. A cluster multivariate analysis identified three different rootstock pairs with similar effects on plant growth, yield and fruit quality of `Okitsu` mandarin. The `Flying Dragon `trifoliate had a unique effect over the `Okitsu` trees performance, inducing lower canopy volume and higher yield efficiency and fruit quality, and might be suitable for high-density plantings. The `Cravo Limeira` and `Cravo FCAV` Rangpur limes induced early-ripening of fruits, with low fruit quality. `Sun Chu Sha Kat` and `Sunki` mandarins and the `Orlando` tangelo conferred lower yield efficiency and less content of soluble solids for the latter rootstock. (C) 2009 Elsevier B.V. All rights reserved.

Water use and crop coefficient for hybrid DKB 390

Water use and crop coefficient for hybrid DKB 390. This work aims to characterize the water use of maize hybrid DKB 390 under suitable conditions of irrigation for both sufficient and below-optimal situations of nitrogen supply. Crop coefficient values for different stages are also presented as a result, in order to provide the basis for crop water budget and management throughout the cycle. A field experiment was carried Out during the main season, in which biomass, soil moisture, leaf area, climate data and light transmittance were evaluated. These have allowed deriving water balance, use and efficiency. The mentioned genotype requires around 600 nun for high yield targets, being less efficient when led under below-optimal nitrogen fertilization.

Modeling of transpiration reduction in van Genuchten-Mualem type soils

We derive an analytic expression for the matric flux potential (M) for van Genuchten-Mualem (VGM) type soils which can also be written in terms of a converging infinite series. Considering the first four terms of this series, the accuracy of the approximation was verified by comparing it to values of M estimated by numerical finite difference integration. Using values of the parameters for three soils from different texture classes, the proposed four-term approximation showed an almost perfect match with the numerical solution, except for effective saturations higher than 0.9. Including more terms reduced the discrepancy but also increased the complexity of the equation. The four-term equation can be used for most applications. Cases with special interest in nearly saturated soils should include more terms from the infinite series. A transpiration reduction function for use with the VGM equations is derived by combining the derived expression for M with a root water extraction model. The shape of the resulting reduction function and its dependency on the derivative of the soil hydraulic diffusivity D with respect to the soil water content theta is discussed. Positive and negative values of dD/d theta yield concave and convex or S-shaped reduction functions, respectively. On the basis of three data sets, the hydraulic properties of virtually all soils yield concave reduction curves. Such curves based solely on soil hydraulic properties do not account for the complex interactions between shoot growth, root growth, and water availability.

Seed development, yield and quality of two palm species growing in different tropical forest types in SE Brazil: implications for ecological restoration

Natural forest remnants have been set as seed production fields to supply seeds of native tree species for tropical forest restoration, but the effect of different forest types on seed production has not been accessed to date for palm species. In this work, we studied seed development, yield, and quality of two palm species in different tropical forest types in SE Brazil. Seed production of palmiteiro (Euterpe edulis) and queen-palm (Syagrus romanzoffiana), which are largely used in restoration efforts due to their importance for vertebrate frugivores, were studied in natural remnants of Atlantic Rainforest, Restinga Forest, Seasonally Dry Forest, and Cerrado Forest. We studied seed development, yield, size, and germination of seed lots produced in some of these forest types, including seeds harvested in 2008, 2009, and both years. Seed yield and quality, as well as seed dry mass in 2009, were higher for palmiteiro seeds produced in the Atlantic Rainforest, while queen-palm seeds produced at the Restinga Forest showed the higher mass and yield, but the lowest physiological potential. Consequently, these natural differences of seed yield and quality have to be taken into account for establishing standards for seed commercialization and analysis, seed pricing, and seedling production in forest nurseries.

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.

Milk yield, milk composition, eating behavior, and lamb performance of ewes fed diets containing soybean hulls replacing coastcross (Cynodon species) hay

The objective of this experiment was to evaluate the effects of replacing coastcross hay NDF by soybean hull (SH) NDF on the lactation performance and eating behavior of ewes and also on the performance of their lambs. Fifty-six Santa Ines lactating ewes (56.1 +/- 6.8 kg of initial BW; mean +/- SD) were penned individually and used in a randomized complete block design with 14 blocks and 4 treatments. Diets were formulated to provide similar concentrations of NDF (56%) and CP (16%). The SH NDF replaced 33 (SH33), 67 (SH67), or 100% (SH100) of the NDF contributed by coastcross hay in a 70% forage-based diet (SH0), resulting in SH inclusion rates of 0, 25, 54, and 85% of the dietary DM. Once a week, from the second to the eighth week of lactation (weaning time), ewes were separated from their lambs, stimulated by a 6-IU i.v. oxytocin injection, and hand milked to empty the udder. After 3 h, milk production was obtained after the same procedure. Quadratic effect for milk production (142.4, 179.8, 212.6, and 202.9 g/3 h) and cubic effect for DMI (2.27, 2.69, 3.25, and 3.00 kg/d) were observed as SH inclusion increased from 0 to 85% of the dietary DM. Milk fat (7.59, 7.86, 7.59, and 7.74%), protein (4.53, 4.43, 4.40, and 4.55%), and total solids (18.24, 18.54, 18.39, and 18.64%) did not differ among the 70% forage-based diet and diets with SH NDF replacing 33, 67, or 100% of the NDF. A linear increase in lactose concentration was observed with SH inclusion. Ewe BW gain during the trial showed a cubic response (0.37, 0.03, 4.80, and 2.80 kg) with SH inclusion. The preweaning ADG of lambs increased linearly, and ADG of lambs after weaning decreased linearly with SH inclusion. Final BW of lambs (2 wk after weaning) did not differ among treatments. Eating behavior observations were conducted with 44 ewes. The same facilities, experimental design, dietary treatments, and feeding management were used. Observations were visually recorded every 5 min for a 24-h period when ewes were 46 +/- 6.8 d in milk. Eating time (min/d, min/g of DMI, and min/g of NDF intake) and time expended in rumination and chewing activities (min/g of DMI and min/g of NDF intake) decreased linearly with the addition of SH in the diets. The inclusion of SH improved DMI and milk production, also reflecting on the BW of lambs at weaning. Milk performance was not affected when SH NDF replaced 100% of hay NDF.

Adapting the CROPGRO perennial forage model to predict growth of Brachiaria brizantha

Warm-season grasses are economically important for cattle production in tropical regions, and tools to aid in management and research of these forages would be highly beneficial. Crop simulation models synthesize numerous physiological processes and are important research tools for evaluating production of warm-season grasses. This research was conducted to adapt the perennial CROPGRO Forage model to simulate growth of the tropical species palisadegrass [Brachiaria brizantha (A. Rich.) Stapf. cv. Xaraes] and to describe model adaptation for this species. In order to develop the CROPGRO parameters for this species, we began with values and relationships reported in the literature. Some parameters and relationships were calibrated by comparison with observed growth, development, dry matter accumulation and partitioning during a 2-year experiment with Xaraes palisadegrass in Piracicaba, SP, Brazil. Starting with parameters for the bahiagrass (Paspalum notatum Flugge) perennial forage model, dormancy effects had to be minimized, and partitioning to storage tissue/root decreased, and partitioning to leaf and stem increased to provide for more leaf and stem growth and less root. Parameters affecting specific leaf area (SLA) and senescence of plant tissues were improved. After these changes were made to the model, biomass accumulation was better simulated, mean predicted herbage yield per cycle was 3573 kg ha(-1), with a RMSE of 538 kg DM ha(-1) (D-Stat = 0.838, simulated/observed ratio = 1.028). The results of the adaptation suggest that the CROPGRO model is an efficient tool to integrate physiological aspects of palisadegrass and can be used to simulate growth. (C) 2010 Elsevier B.V. All rights reserved.

SOYBEAN YIELD COMPONENTS IN DIFFERENT SOWING PERIODS IN THE WESTERN OF STATE OF BAHIA

Yield is closely linked to the plant yield components and depend directly of the genotype interaction with the environment. Essays were installed in the experimental field of the Bahia Foundation in the Maria Gabriela farm in the county of Sao Desiderio - BA in the year 2006-2007. The aim of this work was to evaluate yield components of five soybean cultivars with different maturation cycles indicated for the Western Region of Bahia in different sowing periods. The experimental design was in random blocks in 4 x 5 factorial scheme (four sowing periods: Ep1 first - 11/29/2006, Ep2 second - 12/14/2006, Ep3 third - 12/28/2007, Ep4 fourth - 01/12/2007 and five cultivars: M-SOY 8411, BRS Corisco, BRS 263 [Diferente], BRS Barreiras e M-SOY 9350) with four repetitions. The following characteristics were evaluated: total number of pods per plant, total number of beans per plant, mass of 1000 beans and yield. Plant yield components, total number of pods per plant, total number of beans per plant, mass of 1000 beans, reduced with sowing delay and showed compensation effect between cultivars and sowing periods. Late sowing Ep3 (28/12/2006) and Ep4 (12/01/2007) were not favourable to raise yield of soybean in the Western Region of Bahia.