Evaluation of statistical and geostatistical models of digital soil properties mapping in tropical mountain regions

Soil properties have an enormous impact on economic and environmental aspects of agricultural production. Quantitative relationships between soil properties and the factors that influence their variability are the basis of digital soil mapping. The predictive models of soil properties evaluated in this work are statistical (multiple linear regression-MLR) and geostatistical (ordinary kriging and co-kriging). The study was conducted in the municipality of Bom Jardim, RJ, using a soil database with 208 sampling points. Predictive models were evaluated for sand, silt and clay fractions, pH in water and organic carbon at six depths according to the specifications of the consortium of digital soil mapping at the global level (GlobalSoilMap). Continuous covariates and categorical predictors were used and their contributions to the model assessed. Only the environmental covariates elevation, aspect, stream power index (SPI), soil wetness index (SWI), normalized difference vegetation index (NDVI), and b3/b2 band ratio were significantly correlated with soil properties. The predictive models had a mean coefficient of determination of 0.21. Best results were obtained with the geostatistical predictive models, where the highest coefficient of determination 0.43 was associated with sand properties between 60 to 100 cm deep. The use of a sparse data set of soil properties for digital mapping can explain only part of the spatial variation of these properties. The results may be related to the sampling density and the quantity and quality of the environmental covariates and predictive models used.

Revisiting Field Capacity (FC): variation of definition of FC and its estimation from pedotransfer functions

Taking into account the nature of the hydrological processes involved in in situ measurement of Field Capacity (FC), this study proposes a variation of the definition of FC aiming not only at minimizing the inadequacies of its determination, but also at maintaining its original, practical meaning. Analysis of FC data for 22 Brazilian soils and additional FC data from the literature, all measured according to the proposed definition, which is based on a 48-h drainage time after infiltration by shallow ponding, indicates a weak dependency on the amount of infiltrated water, antecedent moisture level, soil morphology, and the level of the groundwater table, but a strong dependency on basic soil properties. The dependence on basic soil properties allowed determination of FC of the 22 soil profiles by pedotransfer functions (PTFs) using the input variables usually adopted in prediction of soil water retention. Among the input variables, soil moisture content θ (6 kPa) had the greatest impact. Indeed, a linear PTF based only on it resulted in an FC with a root mean squared residue less than 0.04 m³ m-3 for most soils individually. Such a PTF proved to be a better FC predictor than the traditional method of using moisture content at an arbitrary suction. Our FC data were compatible with an equivalent and broader USA database found in the literature, mainly for medium-texture soil samples. One reason for differences between FCs of the two data sets of fine-textured soils is due to their different drainage times. Thus, a standardized procedure for in situ determination of FC is recommended.

ASPECTS OF THE SILVOPASTORAL SYSTEM CORRELATED WITH PROPERTIES OF A TYPIC QUARTZIPSAMMENT (ENTISOL) IN MATO GROSSO DO SUL, BRAZIL

In Brazil, grazing mismanagement may lead to soil and pasture degradation. To impede this process, integrated cropping systems such as silvopasture have been an effective alternative, allied with precision agriculture based on soil mapping for site-specific management. In this study, we aimed to define the soil property that best sheds light on the variability of eucalyptus and forage yield. The experiment was conducted in the 2011/12 crop year in Ribas do Rio Pardo, Mato Grosso do Sul State, Brazil. We analyzed linear and spatial correlations between eucalyptus traits and physical properties of a Typic Quartzipsamment at two depths (0.00-0.10 and 0.10-0.20 m). For that purpose, we set up a geostatistical grid for collection at 72 points. Gravimetric moisture in the 0.00-0.10 m layer is an important index of soil physical quality, showing correlation to eucalyptus circumference at breast height (CBH) in a Typic Quartzipsamment. With an increase in resistance to penetration in the soil surface layer, there is an increase in eucalyptus height and in neutral detergent fiber content in the forage crop. From a spatial point of view, the height of eucalyptus and the neutral detergent fiber of forage can be estimated by co-kriging analysis with soil resistance to penetration. Resistance to penetration values above 2.3 MPa indicated higher yielding sites.

Infrared thermometry to schedule irrigation of common bean

The objective of this work was to determine the critical irrigation time for common bean (Phaseolus vulgaris L. cv. Carioca) using infrared thermometry. Five treatments were analyzed. Canopy temperature differences between plants and a well-watered control about 1, 2, 3, 4, and 5±0.5ºC were tested. Physiological variables and plant growth were analyzed to establish the best time to irrigate. There was a significant linear correlation between the index and stomatal resistance, transpiration rate, and leaf water potential. Although significant linear correlation between the index and mean values of total dry matter, absolute growth rate, and leaf area index was found, no correlation was found with other growth index like relative growth rate, net assimilation rate, and leaf area ratio. Plants irrigated when their canopy temperature was 3±0.5ºC above the control had their relative growth rate mean value increased up to 59.7%, yielding 2,260.2 kg ha-1, with a reduction of 38.0% in the amount of water used. Plants irrigated when their canopy temperature was 4±0.5ºC yielded 1,907.6 kg ha-1, although their relative growth rate mean value was 4.0% below the control. These results show that the best moment to irrigate common bean is when their canopy temperature is between 3ºC and 4±0.5ºC above the control.

Spatial statistical analysis and selection of genotypes in plant breeding

The objective of this study was to evaluate the efficiency of spatial statistical analysis in the selection of genotypes in a plant breeding program and, particularly, to demonstrate the benefits of the approach when experimental observations are not spatially independent. The basic material of this study was a yield trial of soybean lines, with five check varieties (of fixed effect) and 110 test lines (of random effects), in an augmented block design. The spatial analysis used a random field linear model (RFML), with a covariance function estimated from the residuals of the analysis considering independent errors. Results showed a residual autocorrelation of significant magnitude and extension (range), which allowed a better discrimination among genotypes (increase of the power of statistical tests, reduction in the standard errors of estimates and predictors, and a greater amplitude of predictor values) when the spatial analysis was applied. Furthermore, the spatial analysis led to a different ranking of the genetic materials, in comparison with the non-spatial analysis, and a selection less influenced by local variation effects was obtained.

Simulating maize phenology as a function of air temperature with a linear and a nonlinear model

The objective of this study was to adapt a nonlinear model (Wang and Engel - WE) for simulating the phenology of maize (Zea mays L.), and to evaluate this model and a linear one (thermal time), in order to predict developmental stages of a field-grown maize variety. A field experiment, during 2005/2006 and 2006/2007 was conducted in Santa Maria, RS, Brazil, in two growing seasons, with seven sowing dates each. Dates of emergence, silking, and physiological maturity of the maize variety BRS Missões were recorded in six replications in each sowing date. Data collected in 2005/2006 growing season were used to estimate the coefficients of the two models, and data collected in the 2006/2007 growing season were used as independent data set for model evaluations. The nonlinear WE model accurately predicted the date of silking and physiological maturity, and had a lower root mean square error (RMSE) than the linear (thermal time) model. The overall RMSE for silking and physiological maturity was 2.7 and 4.8 days with WE model, and 5.6 and 8.3 days with thermal time model, respectively.

Análise genética de escores de avaliação visual de bovinos com modelos bayesianos de limiar e linear

O objetivo deste trabalho foi comparar as estimativas de parâmetros genéticos obtidas em análises bayesianas uni-característica e bi-característica, em modelo animal linear e de limiar, considerando-se as características categóricas morfológicas de bovinos da raça Nelore. Os dados de musculosidade, estrutura física e conformação foram obtidos entre 2000 e 2005, em 3.864 animais de 13 fazendas participantes do Programa Nelore Brasil. Foram realizadas análises bayesianas uni e bi-características, em modelos de limiar e linear. De modo geral, os modelos de limiar e linear foram eficientes na estimação dos parâmetros genéticos para escores visuais em análises bayesianas uni-características. Nas análises bi-características, observou-se que: com utilização de dados contínuos e categóricos, o modelo de limiar proporcionou estimativas de correlação genética de maior magnitude do que aquelas do modelo linear; e com o uso de dados categóricos, as estimativas de herdabilidade foram semelhantes. A vantagem do modelo linear foi o menor tempo gasto no processamento das análises. Na avaliação genética de animais para escores visuais, o uso do modelo de limiar ou linear não influenciou a classificação dos animais, quanto aos valores genéticos preditos, o que indica que ambos os modelos podem ser utilizados em programas de melhoramento genético.

Growth performance and body composition of giant trahira fingerlings fed diets with different protein and energy levels

The objective of this work was to determine the proper levels of protein and energy in diets of Hoplias lacerdae fingerlings. The dietary crude protein (CP) and gross energy (GE) levels for fingerlings of giant trahira were evaluated in a completely randomized 4x3 factorial design with 35, 39, 43 and 47% CP and 4,100, 4,300 and 4,500 kcal kg-1 of GE, and four replicates. The survival rate was 99.22%, and a linear improvement on the performance parameters was detected after increasing diet crude protein levels. Feed conversion ratio decreased with increasing levels of dietary protein and energy in the diets. A significant interaction between crude protein and gross energy was observed over body protein and mineral matter. Body lipid has increased linearly as gross energy in the diet increased. The retention of crude protein and energy showed a linear increasing with rising of crude protein levels in the diet. Crude protein level at 47% provides the best performance and energy retention, independently of the gross energy levels in the diet.

Tamanho de amostra para estimação do coeficiente de correlação linear de Pearson entre caracteres de milho

O objetivo deste trabalho foi determinar o tamanho de amostra para a estimação do coeficiente de correlação linear de Pearson entre caracteres de três híbridos de milho. Para as análises, foram tomadas aleatoriamente 361, 373 e 416 plantas, respectivamente, de híbridos simples, triplo e duplo. Para cada planta, os seguintes caracteres foram mensurados: diâmetro maior e menor do colmo, altura da planta e altura, peso, comprimento e diâmetro da espiga, número de fileiras por espiga, peso e diâmetro de sabugo, massa de cem grãos, número de grãos por espiga, comprimento e produtividade de grãos. Para cada um dos 91 pares de caracteres e híbridos, foi determinado o tamanho de amostra a partir de "bootstrap", com reposição de 1.000 amostras, de cada tamanho de amostra simulado. Na estimação do coeficiente de correlação linear de Pearson com a mesma precisão, o tamanho de amostra (número de plantas) aumenta na direção de pares de caracteres com menor intensidade de relação linear, independentemente do tipo de híbrido. Para os 91 pares de caracteres estudados, 252 plantas são suficientes para a estimação do coeficiente de correlação linear de Pearson, no intervalo de confiança de "bootstrap" de 95%, igual a 0,30

Dry matter production of Tanzania grass as a function of agrometeorological variables

The objective of this work was to develop and validate linear regression models to estimate the production of dry matter by Tanzania grass (Megathyrsus maximus, cultivar Tanzania) as a function of agrometeorological variables. For this purpose, data on the growth of this forage grass from 2000 to 2005, under dry‑field conditions in São Carlos, SP, Brazil, were correlated to the following climatic parameters: minimum and mean temperatures, degree‑days, and potential and actual evapotranspiration. Simple linear regressions were performed between agrometeorological variables (independent) and the dry matter accumulation rate (dependent). The estimates were validated with independent data obtained in São Carlos and Piracicaba, SP, Brazil. The best statistical results in the development and validation of the models were obtained with the agrometeorological parameters that consider thermal and water availability effects together, such as actual evapotranspiration, accumulation of degree‑days corrected by water availability, and the climatic growth index, based on average temperature, solar radiation, and water availability. These variables can be used in simulations and models to predict the production of Tanzania grass.

Predicting the occurrence of alate aphids in Brassicaceae

The objective of this work was to predict the occurrence of alates of Brevicoryne brassicae, Lipaphis erysimi, and Myzus persicae (Hemiptera, Aphididae) in Brassicaceae. The alate aphids were collected in yellow water traps from July 1997 to August 2005. Aphid population peaks were predicted using a degree‑day model. The meteorological factors, temperature, air relative humidity, rainfall, and sunshine hours, were used to provide precision indexes to evaluate the best predictor for the date of the first capture of alate aphids by the traps. The degree‑day model indicated that the peak population of the evaluated aphid species can be predicted using one of the following biofix dates: January 1st, June 1st, and the date of the first capture of the alate aphid species by the yellow water traps. The best predictor of B. brassicae occurrence is the number of days with minimum temperature >15°C, and of L. erysimi and M. persicae, the number of days with rainfall occurrence.

Degree of multicollinearity and variables involved in linear dependence in additive-dominant models

The objective of this work was to assess the degree of multicollinearity and to identify the variables involved in linear dependence relations in additive-dominant models. Data of birth weight (n=141,567), yearling weight (n=58,124), and scrotal circumference (n=20,371) of Montana Tropical composite cattle were used. Diagnosis of multicollinearity was based on the variance inflation factor (VIF) and on the evaluation of the condition indexes and eigenvalues from the correlation matrix among explanatory variables. The first model studied (RM) included the fixed effect of dam age class at calving and the covariates associated to the direct and maternal additive and non-additive effects. The second model (R) included all the effects of the RM model except the maternal additive effects. Multicollinearity was detected in both models for all traits considered, with VIF values of 1.03 - 70.20 for RM and 1.03 - 60.70 for R. Collinearity increased with the increase of variables in the model and the decrease in the number of observations, and it was classified as weak, with condition index values between 10.00 and 26.77. In general, the variables associated with additive and non-additive effects were involved in multicollinearity, partially due to the natural connection between these covariables as fractions of the biological types in breed composition.

Random regression models to estimate genetic parameters for milk production of Guzerat cows using orthogonal Legendre polynomials

The objective of this work was to compare random regression models for the estimation of genetic parameters for Guzerat milk production, using orthogonal Legendre polynomials. Records (20,524) of test-day milk yield (TDMY) from 2,816 first-lactation Guzerat cows were used. TDMY grouped into 10-monthly classes were analyzed for additive genetic effect and for environmental and residual permanent effects (random effects), whereas the contemporary group, calving age (linear and quadratic effects) and mean lactation curve were analized as fixed effects. Trajectories for the additive genetic and permanent environmental effects were modeled by means of a covariance function employing orthogonal Legendre polynomials ranging from the second to the fifth order. Residual variances were considered in one, four, six, or ten variance classes. The best model had six residual variance classes. The heritability estimates for the TDMY records varied from 0.19 to 0.32. The random regression model that used a second-order Legendre polynomial for the additive genetic effect, and a fifth-order polynomial for the permanent environmental effect is adequate for comparison by the main employed criteria. The model with a second-order Legendre polynomial for the additive genetic effect, and that with a fourth-order for the permanent environmental effect could also be employed in these analyses.