Geostatistics allows the evaluation of the distribution pattern of data with high spatial variability in agricultural systems. This study aimed to evaluate the spatial variability of biological diversity indices of soil fauna under different land (agriculture and forest). Samples were collected in seven areas (millet, soybean, corn, eucalyptus, pasture crops, and preserved and disturbed Cerrado), in Maranhão state, Brazil. The soil fauna was caught trapped in pitfall traps, installed 3 m away from each other. In each area, 130 traps were maintained for seven days. After this period, they were removed and their content transferred to bottles and taken to the laboratory, where the insects were screened and identified at the level of orders and families. Eight indices were calculated, namely: individuals trap-1 day-1 , Jackknife richness estimator, the Simpson, McIntosh, Shannon, and total diversity, and Simpson dominance, and Pielou equitability indices. The spatial variability was derived from the semivariograms fitted to Gaussian, spherical, and exponential geostatistical models. Statistical analysis showed medium values of the coefficient of variation for millet, except for the indices individuals trap-1 day-1 and McIntosh diversity, which were considered high. The values of the correlation matrix were negative for some indices, suggesting an inverse relationship. For millet, corn, eucalyptus, disturbed Cerrado, and pasture areas, the Shannon diversity index exhibited a pure nugget effect. For the areas of millet, corn, disturbed Cerrado and pasture, the total diversity index was adjusted to the Gaussian model. The degree of spatial dependence was considered high for the individuals trap-1 day-1 and Pielou equitability indices for millet. Only for soybean and pasture similarity in the scaled semivariograms was observed for the spatial variability of the indices, indicating similarity of performance. Soil management and land use affect the patterns of soil fauna abundance, richness, and diversity. The presence of groups such as Araneae, Diplura, and Poduromorpha are related to ecological equilibrium, quality, and sustainability of the agricultural systems studied.

Soil resistance to penetration (PR) is an indirect measure of the state of soil compaction. Thus, the objective of this study was to characterize PR in vertical profiles in an area cultivated with sugarcane using multifractal models for different relief units. The experiment was carried out in an Oxisol with a clay texture, with 6.85 ha in the municipality of Coelho Neto (Maranhão state, Brazil), where 60 sampling points were demarcated. The area was divided into four relief units (Type A > 74 m, Type B from 71 to 74 m, Type C from 68 to 71 m and Type D from 65 to 68 m). The PR was measured at the 60 sampling points using an impact penetrometer, and the PR determined in the 0-0.60 m depth layer every 0.01 m. The multifractal analysis was performed considering the scale property of each profile and typified the singularity and Rènyi spectra estimated using the current method. Multifractal analysis allowed the identification of patterns at different scales and with high heterogeneity. The multifractal behavior was represented by the singularity spectrum (α), versus f(α), and the generalized dimension (Dq). The multifractal analysis allowed the differentiation between the profiles of the relief units (Types A, B, C and D), resulting in an important tool for studies of soil resistance to penetration.

Soil beetles' communities are responsible for many ecosystem services, and are very sensitive to environmental changes. Thus, this study aimed to evaluate the abundance and diversity of the soil coleoptera fauna under uses and management and also to identify relationships of the beetle community with soil's physical and chemical properties. The experiment had six experimental plots set up an Oxisol (Latossolo): corn (CO), soybean (SO), 7-year-old eucalyptus (EI), 4-year-old eucalyptus (EII), preserved Cerrado (PC), and disturbed Cerrado (DC). Soil beetles were sampled at 128 points for each experimental plot, where the soil physical and chemical properties were analyzed. The Coleoptera fauna organisms were identified at the family, subfamily, and gender level, and then, the number of individuals per day, richness, Shannon diversity indexes, and Pielou evenness were determined. The data were analyzed using multivariate techniques (hierarchical grouping and factor analysis). On total, 750 specimens of beetles were collected, distributed into 9 families, 14 subfamilies, and 27 genera. The most abundant family was Scarabaeidae (11 genera) with the highest occurrence in the PC (143 specimens) and DC (81 specimens). Cultivation with SO presented the greatest number of trap day individuals (ind trap -1 day -1 = 0.548); however, the highest diversity was found in the PC. (20 taxonomic groups) and CO (16 taxonomic groups). Shannon diversity was higher for the CO (H' = 3.107), followed by the PC (H' = 2.699), and the lowest value was found for the SO (H' = 1.530). The similarity dendrogram grouped the plots into two extracts, demonstrating how the intensity of land use influences the abundance and diversity of beetle fauna. The factor analysis grouped the Coleoptera and the physical and chemical soil properties in two factors: elements related to the state of aggregation and porous system's elements. The Coleoptera community was influenced by the intensity of land use and the portion with anthropized natural vegetation showed the highest richness, demonstrating that the Coleoptera fauna responds to environmental changes. Edaphic beetles in the different use and management systems were primarily related to soil physical properties, which explain the state of aggregation (pH, altitude, Ca 2+ , BD, clay, macroporosity, silt, K + , and microporosity) and the porous soil system (sand and total porosity).

In the global agribusiness, the herbicide use is a major problem for sustainable production, in this sense, it is necessary to better understand the interaction of weed species and floristic composition such as biodiversity indicators. The objective of this study was to analyze the spatial variability of weeds in an Oxisol under no-tillage system. Samples were taken in an area of 0.5 ha, in 50 sampling points with spacing of 5 m x 10 m. Data were analyzed by means of classical statistics, geostatistics, and spatial variability of the constructed maps by the interpolation by kriging technique. All the species of weeds presented in the study area showed spatial variability with the exception of Ipomoea triloba (L.) and Heliotropium indicum (L.), which showed pure nugget effect. The range values (a) shows that the spacing between samples can be extended to all species of weeds. The study was unable to determine specifics areas of management in the local since the different species of weed infested different plots of the area.

The objective of this study was to determine the multifractality of diversity indexes of edaphic fauna in areas with natural vegetation and in agricultural systems. Biological sampling was carried out in seven treatments (millet, maize, soybean, eucalyptus, preserved cerrado, disturbed cerrado and pasture), containing 130 pitfall traps, distributed in transects with 3 m of spacing between sampling points, totaling 390 m. The multifractal analysis was determined based on the moment method, where estimates such as the capacity dimensions, entropy and the correlation of the diversity indexes were calculated. The soybean area had greater Shannon diversity (2.69), however it had smaller abundance of individuals. The partition functions were adjusted with coefficient of determination > 0.90. The difference between D-10-D 10 , ranged from 0.080 to 1,707 for Pielou equitability in soybean cultivation, for richness in the area under the cultivation of eucalyptus. The singularity spectra expressed graphs with different degrees of heterogeneity for the soil fauna indexes, and the richness expressed the best structure. The area cultivated with soybean had a monofractal tendency, due to the homogeneous distribution of individuals of the edaphic fauna along the transect. The fractal analysis provided the description of patterns of variability that are not detected by classical methods.