Implications of converting native forest areas to agricultural systems on the dynamics of CO2 emission and carbon stock in a Cerrado soil, Brazil.

The conversion of native vegetation to agricultural areas leads to a natural process of carbon loss but these systems can stabilize in terms of carbon dynamics depending on the management and conversion time, presenting potential to both store and stabilize this carbon in the soil, resulting in lower soil respiration rates. In this context, this study aimed to investigate the effect of converting native Cerrado forest areas to agricultural systems with a forest planted with Eucalyptus camaldulensis and silvopastoral systems on the dynamics of CO2 emission and carbon stock at different soil depths. The experimental sites are located in the Midwest of Brazil, in the coordinates 20°22'31″ S and 51°24'12″ W. Were evaluated soil CO2 emission (FCO2), soil organic carbon, the degree of humification of soil organic matter (HLIFS), soil temperature, soil moisture, and soil chemical and physical attributes. The soil of the area is classified as an Oxisol (Haplic Acrustox). Soil samples were collected at depths of 0.00-0.10, 0.10-0.20, 0.20-0.30, and 0.30-0.40 m. The lowest FCO2 values were found in the silvopastoral system (1.05 µmol m-2 s-1), followed by the native forest (1.65 µmol m-2 s-1) and the eucalyptus system (1.96 µmol m-2 s-1), indicating a 36% reduction in FCO2 compared to the conversion of the native forest to the silvopastoral system and an increase of 19% when converting the native forest to the eucalyptus system. The soil chemical attributes (N, K+, Ca2+, H++Al3+, CEC, and organic carbon) showed a decrease along the profile. The shallowest depths (0.00-0.10 and 0.10-0.20 m) presented no differences between systems but the subsequent depths (0.20-0.30 and 0.30-0.40 m) had a difference (95% confidence interval), relative to N, Ca2+, H++Al3, CEC, and organic carbon stock (OCS), and the soil under silvopastoral system showed a higher concentration of these attributes than the native forest. The multivariate analysis showed that the eucalyptus and silvopastoral systems did not differ from the forest in the shallowest soil layer but differed from each other. This behavior changed from the second assessed depth (0.10-0.20 m), in which the silvopastoral system stands out, differing both from the eucalyptus system and from the native forest, and this behavior is maintained at the following depths (0.20-0.30 and 0.30-0.40 m). OCS, H++Al3, CEC, and nitrogen are strongly related to land use change for silvopastoral system. Regarding the behavior/relationship of attributes as a function of depth, the silvopastoral system contributed to soil carbon accumulation and stability over consecutive years.

Artificial neural networks and adaptive neuro-fuzzy inference systems for prediction of soil respiration in forested areas southern Brazil.

The purpose of this study was to estimate the temporal variability of CO2 emission (FCO2) from O2 influx into the soil (FO2) in a reforested area with native vegetation in the Brazilian Cerrado, as well as to understand the dynamics of soil respiration in this ecosystem. The database is composed of soil respiration data, agroclimatic variables, improved vegetation index (EVI), and soil attributes used to train machine learning algorithms: artificial neural network (ANN) and an adaptive neuro-fuzzy inference system (ANFIS). The predictive performance was evaluated based on the mean absolute error (MEA), root mean square error (RMSE), mean absolute percentage error (MAPE), agreement index (d), confidence coefficient (c), and coefficient of determination (R2). The best estimation results for validation were FCO2 with multilayer perceptron neural network (MLP) (R2 = 0.53, RMSE = 0.967 µmol m-2 s-1) and radial basis function neural network (RBF) (R2 = 0.54, RMSE = 0.884 µmol m-2 s-1) and FO2 with MLP (R2 = 0.45, RMSE = 0.093 mg m-2 s-1) and RBF (R2 = 0.74, 0.079 mg m-2 s-1). Soil temperature and macroporosity are important predictors of FCO2 and FO2. The best combination of variables for training the ANFIS was selected based on trial and error. The results were as follows: FCO2 (R2 = 16) and FO2 (R2 = 29). In all models, FCO2 outperformed FO2. A primary factor analysis was performed, and FCO2 and FO2 correlated best with the weather and soil attributes, respectively.

Machine learning for prediction of soil CO2 emission in tropical forests in the Brazilian Cerrado.

Soil CO2 emission (FCO2) is a critical component of the global carbon cycle, but it is a source of great uncertainty due to the great spatial and temporal variability. Modeling of soil respiration can strongly contribute to reducing the uncertainties associated with the sources and sinks of carbon in the soil. In this study, we compared five machine learning (ML) models to predict the spatiotemporal variability of FCO2 in three reforested areas: eucalyptus (RE), pine (RP) and native species (RNS). The study also included a generalized scenario (GS) where all the data from RE, RP and RNS were included in one dataset. The ML models include generalized regression neural network (GRNN), radial basis function neural network (RBFNN), multilayer perceptron neural network (MLPNN), adaptive neuro-fuzzy inference system (ANFIS) and random forest (RF). Initially, we had 32 attributes and after pre-processing, including Pearson's correlation, canonical correlation analysis (CCA), and biophysical justification, only 21 variables remained. We used as input variables 19 soil properties and climate variables in reforested areas of eucalyptus, pine and native species. RF was the best model to predict soil respiration to RE [adjusted coefficient of determination (R2 adj): 0.70 and root mean square error (RMSE): 1.02 µmol m-2 s-1], RP (R2 adj: 0.48 and RMSE: 1.07 µmol m-2 s-1) and GS (R2 adj: 0.70 and RMSE: 1.05 µmol m-2 s-1). Our findings support that RF and GRNN are promising for predicting soil respiration of reforested areas which could help to identify and monitor potential sources and sinks of the main additional greenhouse gas over ecosystems.

The fewer, the better fare: Can the loss of vegetation in the Cerrado drive the increase in dengue fever cases infection?

Several studies have reported the relationship of deforestation with increased incidence of infectious diseases, mainly due to the deregulation caused in these environments. The purpose of this study was to answer the following questions: a) is increased loss of vegetation related to dengue cases in the Brazilian Cerrado? b) how do different regions of the tropical savanna biome present distinct patterns for total dengue cases and vegetation loss? c) what is the projection of a future scenario of deforestation and an increased number of dengue cases in 2030? Thus, this study aimed to assess the relationship between loss of native vegetation in the Cerrado and dengue infection. In this paper, we quantify the entire deforested area and dengue infection cases from 2001 to 2019. For data analyses, we used Poisson generalized linear model, descriptive statistics, cluster analysis, non-parametric statistics, and autoregressive integrated moving average (ARIMA) models to predict loss of vegetation and fever dengue cases for the next decade. Cluster analysis revealed the formation of four clusters among the states. Our results showed significant increases in loss of native vegetation in all states, with the exception of Piauí. As for dengue cases, there were increases in the states of Minas Gerais, São Paulo, and Mato Grosso. Based on projections for 2030, Minas Gerais will register about 4,000 dengue cases per 100,000 inhabitants, São Paulo 750 dengue cases per 100,000 inhabitants, and Mato Grosso 500 dengue cases per 100,000 inhabitants. To reduce these projections, Brazil will need to control deforestation and implement public health, environmental and social policies, requiring a joint effort from all spheres of society.

Contribution of morphological traits for grain yield in common bean / Contribuição de caracteres morfoagronômicos para produtividade de grãos em feijoeiro comum

In Brazil, the common bean crop has representative agricultural exploitation, not only because of its production economic value, but also because there is a large cultivation area. This research aimed to investigate the direct and indirect relationship of morphological components on grain yield in common bean plants. This study was carried out in a Quartzarenic Neosol in the municipality of Cassilândia, Mato Grosso do Sul State, Brazil, in the agricultural year of 2016/2017. The evaluated traits were: grains yield (GY) with its primary components; mass of one hundred grains (HG); number of grains per plant (GP); number of grains per pod (GPP); dry weight of aerial parts (DWA); number of pod per plant (PP); plant dry mass (DM); plant high (PH); and stem diameter (SD). Initially, the Pearson's correlation among these traits was estimated and the correlation network was used to graphically express the obtained results. Analysis of these data through the statistical techniques of multicollinearity diagnosis followed by path analysis enabled to verify that the number of pod per plant, the mass of one hundred grains, and the number of grains per plant, among the primary components of grain yield, are the traits of greater potential to select and identify superior genotypes for grain productivity yield, and that dry matter and stem diameter traits showed a negative correlation with grain yield in common bean grains.

No Brasil, a cultura do feijoeiro tem uma exploração agrícola representativa, não só por seu valor econômico de produção, mas também porque existe grande área de cultivo. Este trabalho teve como objetivo investigar a relação direta e indireta de componentes morfológicos no rendimento de grãos nas plantas de feijão comum. Este estudo foi realizado em Neossolo Quartzarênico no município de Cassilândia, Estado de Mato Grosso do Sul, no ano agrícola de 2016/2017. As características avaliadas foram: produtividade de grãos (GY) com seus componentes primários; massa de cem grãos (HG); número de grãos por planta (GP); número de grãos por vagem (GPP); massa seca da parte aérea (DWA); número de vagens por planta (PP); massa seca da planta (DM); altura de planta (PH); e diâmetro do caule (SD). Inicialmente, a correlação de Pearson entre esses caracteres foi estimada e a rede de correlação foi usada para expressar graficamente os resultados obtidos. A análise desses dados através do diagnóstico de multicolinearidade seguida de análise de trilha permitiu verificar que o número de vagens por planta, a massa de cem grãos e o número de grãos por planta, entre os principais componentes do rendimento de grãos, são os caracteres de maior potencial para selecionar e identificar genótipos superiores para a produtividade de grãos no feijoeiro e que a matéria seca e o diâmetro do caule mostraram uma correlação negativa com a produtividade de grãos em grãos.

Produtividade da cana-de-açúcar: variabilidade linear e espacial entre componentes tecnológicos e da produção / Productivity of cane sugar: linear and spatial variability between components and technological production

Na atualidade, a cultura da cana-de-açúcar desempenha frente à realidade brasileira, sobretudo no aspecto relacionado às fontes energéticas alternativas, importantíssimo papel. No ano de 2009, no município de Suzanápolis (SP), no Cerrado Brasileiro, foi instalado um experimento com a cultura da cana-de-açúcar em um Argissolo Vermelho eutrófico, com o objetivo de selecionar, por meio de coeficientes de correlação de Pearson, da modelagem de regressões lineares simples e de múltiplas e da correlação espacial, os melhores componentes tecnológicos e os de produção, para explicar a variabilidade da produtividade da cana-de-açúcar. Foi instalada a malha geoestatística, para a coleta de dados, com 120 pontos amostrais, numa área de 14,53 ha. Para as regressões lineares simples, a população de plantas é o componente da produção que apresenta a melhor correlação quadrática com a produtividade da cana-de-açúcar, dada por: PRO = -0,553**xPOP2+16,14*xPOP-15,77. Entretanto, para as regressões lineares múltiplas, a equação PRO = - 21,11+4,92xPOP**+0,76xPUR** é a que melhor se apresenta para estimar a referida produtividade. Espacialmente, a melhor correlação com a produtividade da cana-de-açúcar é determinada, também, pelo componente da produção da população de plantas.

Nowadays, the culture of the sugarcane plays an important role regarding the Brazilian reality, especially in the aspect related to the alternative energy sources. In 2009, the municipality of Suzanápolis (SP), in the Brazilian Cerrado, an experiment was conducted with the culture of the sugarcane in a Red eutrophic, with the aim of selecting, using Pearson correlation coefficients, modeling, simple, linear and multiple regressions and spatial correlation, and also the best technological and productive components, to explain the variability of the productivity of the sugarcane. The geostatistical grid was installed in order to collect the data, with 120 sampling points, in an area of 14.53 ha. For the simple linear regressions, the plants population is the component of production that presents the best quadratic correlation with the productivity of the sugarcane, given by: PRO = -0.553**xPOP2+16.14*xPOP-15.77. However, for multiple linear regressions, the equation PRO = -21.11+4.92xPOP**+0.76xPUR** is the one that best presents in order to estimate that productivity. Spatially, the best correlation with yield of the sugarcane is also determined by the component of the production population of plants.

Espacialização do banco de sementes de plantas daninhas sob diferentes manejos de cana-de-açúcar em rio brilhante, MS / Geographical distribution of the weeds seed bank of under different managements of sugar cane in rio brilhante, MS

O setor sucroalcooleiro está em pleno crescimento no Estado de Mato Grosso do Sul e, com isso a preocupação com os sistemas de colheita passam a ser estudadas, pois estes sistemas podem influenciar na interferência da comunidade infestante de plantas daninhas na cultura da cana-de-açúcar. O manejo integrado unido à ferramenta geoestatística visa evitar perdas de produtividade devido à interferência de plantas daninhas. O objetivo deste trabalho foi estudar a variabilidade espacial do banco de sementes de plantas daninhas em função do sistema de colheita de cana-de-açúcar (sem queima e com queima). O experimento foi realizado na área de cultivo comercial da Usina ETH Bioenergia S/A Unidade Eldorado. As amostras de solo foram retiradas com trado da camada de 0,00 a 0,40 m de profundidade, nos dois sistemas de colheita. A parcela experimental foi composta por uma malha constituída por 50 pontos georreferenciados com distâncias irregulares. As amostras de solo foram levadas para a casa de vegetação para emergência das plântulas. O número de espécies daninhas foi analisado por meio de análises estatísticas descritivas e de técnicas geoestatísticas. As sementes de picão-preto, dicotiledôneas, capim-amargoso, tiririca, trapoeraba e monocotiledôneas apresentaram dependência espacial do banco de sementes no sistema de colheita com queima da cana-de-açúcar. Para o sistema de colheita da cana queimada apenas a espécie tiririca apresentou dependência espacial de distribuição no banco de sementes. Na colheita de cana sem queima o mapeamento dessas espécies possibilitará por meio dos mapas de krigagem produzidos, aplicações localizadas de herbicidas no manejo integrado de Cyperus rotundus.

The State of Mato Grosso do Sul is in full growth of this sector, thus the concern about harvesting systems are being studied, and these systems may influence the weed community interference of weeds in the cane sugar. The integrated management tool attached to geostatistics is to avoid productivity losses due to weed interference. The objective of this work was to study the spatial variability of the seed bank of weeds depending on the system for collecting cane sugar (raw and burning). The experiment was conducted in the area of commercial cultivation of the plant ETH Bioenergy S/A Eldorado Unity. Soil samples were taken with auger layer from 0.00 to 0.40 m depth in both cropping systems. The experimental plot was composed by a mesh consisting of 50 points georeferenced with irregular distances. Soil samples were taken to the greenhouse for germination. The number of weed species was analyzed using descriptive statistics and geostatistical techniques. The seeds of B. pilosa, dicots, bitter grass, nutsedge, dayflower monocots and spatial dependence of the seed bank in the collection system with burning of cane sugar. For the system of harvest only the raw sedge species present spatial dependence of distribution in the seed bank. In the harvest green cane enable the mapping of these species through the kriging maps produced, spot applications of herbicides in integrated management of Cyperus rotundus.

Produtividade da soja correlacionada com a porosidade e a densidade de um Latossolo Vermelho do cerrado brasileiro / Soybean yield correlated with the porosity and the density of a red Latosol of the Brazilian Savanna

No ano agrícola de 2004/05, na Fazenda Experimental da FE/Unesp-Campus de Ilha Solteira, localizada no Município de Selvíria, Mato Grosso do Sul - (MS), foram estudadas a variabilidade e a dependência espacial entre a macroporosidade (MA), a microporosidade (MI), a porosidade total (PT) e a densidade do solo (DS), em três profundidades, com a produtividade da soja (PG), em sistema plantio direto, sobre um Latossolo Vermelho distroférrico. O objetivo foi pesquisar a variabilidade e as correlações lineares e espaciais entre os atributos da planta e do solo, visando a selecionar um indicador da qualidade física do solo para a produtividade da soja. Foi instalada a malha geoestatística, para a coleta de dados do solo e da planta, contendo 124 pontos amostrais, numa área de 4.000m². A correlação linear entre a produtividade de grãos de soja com os atributos do solo estudados foi baixa. Houve correlação espacial direta entre a PG e a MA nas camadas de zero-0,10m e 0,10-0,20m. A macroporosidade do solo (MA1), avaliada na camada de zero-0,10m, apresentou-se como o melhor indicador da qualidade física do solo quando destinada à avaliação da produtividade de soja, nas condições de plantio direto, em solo de cerrado.

The variability and the spatial dependence attributes of the following plant and soil: a) the soybean yield (SY) in no-tillage and irrigated, b) the macroporosity (MA), c) microporosity (MI), d) total porosity (TP), e) and bulk density (BD), in depths of 1 (0-0.10m), 2 (0.10-0.20m) and 3 (0.20-0.30 m), in a Red Latosol of the Engineering college of the São Paulo State University (FE/UNESP) in the Ilha Solteira campus, located in Selvíria County, Mato Grosso do Sul state, Brazil, were analyzed in the agricultural year of 2004/05. The main goal was to study the variability and the linear and spatial correlations among the attributes (plant and soil), aiming to select an index of soil physical quality with good capacity to represent the soybean yield. A geostatistical grid to collect soil and plant data was installed, with one hundred and twenty four sample points, in an area of 4000m². The linear relationship between soybean yield and physical attributes of soil studied was low. So, there was a direct spatial relationship between SY and MA1 (0-0.10m), as well as between SY and MA2 (0.10-0.20m). The MA1 showed its best index of the soil physical quality, when it was destined to estimate the soybean yield on no-tillage.

Variabilidade espacial de atributos químicos em latossolo e argissolos / Spatial variability of chemical attributes in latosol and argisol

O trabalho teve como objetivo caracterizar a variabilidade espacial de atributos químicos de Latossolos e Argissolos, sob cultivo de cana-de-açúcar em áreas com variações na forma do relevo. No presente estudo utilizou-se duas áreas, sendo uma em Latossolo em pedoforma convexa (158ha) e a outra em Argissolo na pedoforma linear (172ha). Foi coletada amostra de solo em malha na profundidade de 0,00-0,50m, realizando-se análise química de cada ponto amostrado. Os maiores coeficientes de variação e alcances foram observados na pedoforma convexa (Latossolo). Portanto, o Latossolo inserido na pedoforma convexa apresentou maior variabilidade espacial para os atributos químicos em relação ao Argissolo na pedoforma linear. O latossolo inserido pedoforma convexa necessita de maior número de pontos de coleta por apresentar maior variabilidade espacial. Recomenda-se que o intervalo de amostragem seja igual ao alcance da dependência espacial, para associar menor esforço de amostragem com maior representatividade.

The objective of this research was to characterize the spatial variability of chemical attributes of soils cultivated with sugarcane in areas with variations in relief forms. The two soils studied were a Latosol, placed in a convex landform (158ha), and an Argisol, placed in a linear landform (172ha). Soil samples for chemical analysis were collected in a regular grid, at the depth of 0.0-0.5m. The convex landform presented the highest coefficient variation and ranges. Therefore, the Latosol inserted in the landform convex presented greater spatial variability for chemical attributes than the Argisol in the linear landform. A larger number of points must be sampled in order to represent the spatial variability of the Latosol from the convex landform. We recommend the sampling interval to be equal to the range of the spatial dependency in order to associate a less intensive sampling effort with a higher representativity.

Influência de formas do relevo em atributos físicos de um latossolo sob cultivo de cana-de-açúcar / Influence of the relief forms on physical attributes of an oxisol cropped with sugarcane

Os atributos físicos do solo variam em função das formas do relevo e sofrem influência da mineralogia da fração argila e do manejo da cultura de cana-de-açúcar, podendo interferir no processo de compactação do solo. Este trabalho foi realizado com o objetivo de avaliar o estado dos atributos físicos em diferentes formas do relevo em um Latossolo Vermelho eutroférrico argiloso sob cultivo de cana-de-açúcar. A área apresenta duas formas de relevo: uma côncava, que ocorre nas posições mais elevadas, e uma linear, constituída pelos segmentos ombro, escarpa, meia encosta e encosta inferior. As amostras de solo foram coletadas durante o ciclo da cultura, nas camadas de 0,00-0,15m, 0,15-0,30m e 0,30-0,45m, para determinação do teor de matéria orgânica e dos seguintes atributos físicos: densidade do solo, porosidade total, macroporosidade, microporosidade, resistência do solo à penetração e teor de água no solo. A mineralogia mais gibbsítica e o maior teor de matéria orgânica, encontrada na forma de relevo côncava e de segmento ombro, proporciona menores valores de densidade do solo, de resistência do solo à penetração e de microporosidade e maiores valores de macroporosidade e de porosidade total do que a mineralogia mais caulinítica, encontrada nos demais segmentos da forma linear.

The physical attributes of the soil vary according to the relief forms and suffer from the mineralogy of the clay fraction and the management of the sugarcane culture, being able to influence the soil compaction process. This research was developed in Jaboticabal-SP, aiming at evaluating the behavior of physical attributes in different relief forms in an oxisol cropped with sugarcane. The area presents two relief forms (concave and linear), the concave one, that occurs at the highest positions, and the linear one, consisting of shoulder, scarp, stocking lean and inferior lean segments. The samples soil were collected during the sugarcane cycle, in the 0.00-0.15, 0.15-0.30 and 0.30-0.45m layers, for the determination of the organic matter and the following physical attributes: bulk density, total porosity, macro and microporosity, soil penetration resistance and soil moisture. The mineralogy most gibbsitic and the greatest value of organic matter, found in the concave relief form and in the shoulder segment, provide lower bulk density, soil penetration resistance and microporosity values and higher macroporosity and total porosity values when compared with the other segments of the linear form (caulinitic mineralogy).

Otimização amostral de atributos de latossolos considerando aspectos solo-relevo / Sampling optimization of oxisol attributes considering soil-relief relationship

O custo de amostragem do solo para caracterizar a variabilidade espacial de atributos químicos e físicos é um problema em agricultura de precisão. O objetivo deste trabalho foi utilizar a forma da paisagem como critério auxiliar na otimização do esquema amostral, a fim de estimar os atributos químicos e físicos de latossolos em áreas sob cultivo de cana-de-açúcar. O solo foi amostrado na profundidade de 0,0-0,2m nos pontos de cruzamento de uma malha, com intervalos regulares de 10m, perfazendo um total de 100 pontos. Todos os atributos químicos e físicos tiveram mais de 84 por cento dos pontos plotados na pedoforma côncava e convexa em relação aos da pedoforma linear. A aplicação do programa Sanos 0.1 na malha amostral permitiu verificar que, na pedoforma côncava e convexa, há maior variabilidade espacial para os atributos químicos e físicos do solo, em comparação à pedoforma linear.

The cost of sampling for spatial variability characterization of soil chemical and physical attributes is a fundamental problem in precision agriculture. The objective of this work was to use the landscape form as an auxiliary criterion for sampling optimization in an Oxisol under sugarcane cultivation, aiming at the estimation of their chemical and physical attributes. Soil samples were collected in the depths of 0.0-0.2m in the intersections of a 10m regular interval grid, comprising a total of 100 points. All chemical and physical attributes had more than 84 percent of the points plotted in the concave and convex landforms compared to the linear landform. Larger spatial variability in the soil of chemical and physical attributes was observed in the concave landform than in the linear one according to analysis performed with the program Sanos 0.1.