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.

Assessment of pedotransfer functions for estimating soil water retention curves for the amazon region

Knowledge of the soil water retention curve (SWRC) is essential for understanding and modeling hydraulic processes in the soil. However, direct determination of the SWRC is time consuming and costly. In addition, it requires a large number of samples, due to the high spatial and temporal variability of soil hydraulic properties. An alternative is the use of models, called pedotransfer functions (PTFs), which estimate the SWRC from easy-to-measure properties. The aim of this paper was to test the accuracy of 16 point or parametric PTFs reported in the literature on different soils from the south and southeast of the State of Pará, Brazil. The PTFs tested were proposed by Pidgeon (1972), Lal (1979), Aina & Periaswamy (1985), Arruda et al. (1987), Dijkerman (1988), Vereecken et al. (1989), Batjes (1996), van den Berg et al. (1997), Tomasella et al. (2000), Hodnett & Tomasella (2002), Oliveira et al. (2002), and Barros (2010). We used a database that includes soil texture (sand, silt, and clay), bulk density, soil organic carbon, soil pH, cation exchange capacity, and the SWRC. Most of the PTFs tested did not show good performance in estimating the SWRC. The parametric PTFs, however, performed better than the point PTFs in assessing the SWRC in the tested region. Among the parametric PTFs, those proposed by Tomasella et al. (2000) achieved the best accuracy in estimating the empirical parameters of the van Genuchten (1980) model, especially when tested in the top soil layer.

Three-dimensional interpolation of soil data: fertility and pedomorphological features in southern Brazil

The graphical representation of spatial soil properties in a digital environment is complex because it requires a conversion of data collected in a discrete form onto a continuous surface. The objective of this study was to apply three-dimension techniques of interpolation and visualization on soil texture and fertility properties and establish relationships with pedogenetic factors and processes in a slope area. The GRASS Geographic Information System was used to generate three-dimensional models and ParaView software to visualize soil volumes. Samples of the A, AB, BA, and B horizons were collected in a regular 122-point grid in an area of 13 ha, in Pinhais, PR, in southern Brazil. Geoprocessing and graphic computing techniques were effective in identifying and delimiting soil volumes of distinct ranges of fertility properties confined within the soil matrix. Both three-dimensional interpolation and the visualization tool facilitated interpretation in a continuous space (volumes) of the cause-effect relationships between soil texture and fertility properties and pedological factors and processes, such as higher clay contents following the drainage lines of the area. The flattest part with more weathered soils (Oxisols) had the highest pH values and lower Al3+ concentrations. These techniques of data interpolation and visualization have great potential for use in diverse areas of soil science, such as identification of soil volumes occurring side-by-side but that exhibit different physical, chemical, and mineralogical conditions for plant root growth, and monitoring of plumes of organic and inorganic pollutants in soils and sediments, among other applications. The methodological details for interpolation and a three-dimensional view of soil data are presented here.

Contribution of Topography and Incident Solar Radiation to Variation of Soil and Plant Litter at an Area with Heterogeneous Terrain

Natural processes that determine soil and plant litter properties are controlled by multiple factors. However, little attention has been given to distinguishing the effects of environmental factors from the effects of spatial structure of the area on the distribution of soil and litter properties in tropical ecosystems covering heterogeneous topographies. The aim of this study was to assess patterns of soil and litter variation in a tropical area that intercepts different levels of solar radiation throughout the year since its topography has slopes predominantly facing opposing geographic directions. Soil data (pH, C, N, P, H+Al, Ca, Mg, K, Al, Na, sand, and silt) and plant litter data (N, K, Ca, P, and Mg) were gathered together with the geographic coordinates (to model the spatial structure) of 40 sampling units established at two sites composed of slopes predominantly facing northwest and southeast (20 units each). Soil and litter chemical properties varied more among slopes within similar geographic orientations than between the slopes facing opposing directions. Both the incident solar radiation and the spatial structure of the area were relevant in explaining the patterns detected in variation of soil and plant litter. Individual contributions of incident solar radiation to explain the variation in the properties evaluated suggested that this and other environmental factors may play a particularly relevant role in determining soil and plant litter distribution in tropical areas with heterogeneous topography. Furthermore, this study corroborates that the spatial structure of the area also plays an important role in the distribution of soil and litter within this type of landscape, which appears to be consistent with the action of water movement mechanisms in such areas.

Development of a Database for Drilled Shaft Foundation Testing (DSHAFT), June 2012

Drilled shafts have been used in the US for more than 100 years in bridges and buildings as a deep foundation alternative. For many of these applications, the drilled shafts were designed using the Working Stress Design (WSD) approach. Even though WSD has been used successfully in the past, a move toward Load Resistance Factor Design (LRFD) for foundation applications began when the Federal Highway Administration (FHWA) issued a policy memorandum on June 28, 2000.The policy memorandum requires all new bridges initiated after October 1, 2007, to be designed according to the LRFD approach. This ensures compatibility between the superstructure and substructure designs, and provides a means of consistently incorporating sources of uncertainty into each load and resistance component. Regionally-calibrated LRFD resistance factors are permitted by the American Association of State Highway and Transportation Officials (AASHTO) to improve the economy and competitiveness of drilled shafts. To achieve this goal, a database for Drilled SHAft Foundation Testing (DSHAFT) has been developed. DSHAFT is aimed at assimilating high quality drilled shaft test data from Iowa and the surrounding regions, and identifying the need for further tests in suitable soil profiles. This report introduces DSHAFT and demonstrates its features and capabilities, such as an easy-to-use storage and sharing tool for providing access to key information (e.g., soil classification details and cross-hole sonic logging reports). DSHAFT embodies a model for effective, regional LRFD calibration procedures consistent with PIle LOad Test (PILOT) database, which contains driven pile load tests accumulated from the state of Iowa. PILOT is now available for broader use at the project website: http://srg.cce.iastate.edu/lrfd/. DSHAFT, available in electronic form at http://srg.cce.iastate.edu/dshaft/, is currently comprised of 32 separate load tests provided by Illinois, Iowa, Minnesota, Missouri and Nebraska state departments of transportation and/or department of roads. In addition to serving as a manual for DSHAFT and providing a summary of the available data, this report provides a preliminary analysis of the load test data from Iowa, and will open up opportunities for others to share their data through this quality–assured process, thereby providing a platform to improve LRFD approach to drilled shafts, especially in the Midwest region.

Diversity and structure of AMF communities as affected by tillage in a temperate soil.

Arbuscular mycorrhizal fungi (AMF) were studied in differently tilled soils from a long-term field experiment in Switzerland. Diversity and structure of AMF communities were surveyed either directly on spores isolated from the field soil or on spores isolated from trap cultures, planted with different host plants. Single-spore cultures were established from the AMF spores obtained from trap cultures. Identification of the AMF was made by observation of spore morphology and confirmed by sequencing of ITS rDNA. At least 17 recognised AMF species were identified in samples from field and/or trap cultures, belonging to five genera of AMF--Glomus, Gigaspora, Scutellospora, Acaulospora, and Entrophospora. Tillage had a significant influence on the sporulation of some species and non- Glomus AMF tended to be more abundant in the no-tilled soil. The community structure of AMF in the field soil was significantly affected by tillage treatment. However, no significant differences in AMF diversity were detected among different soil tillage treatments. AMF community composition in trap cultures was affected much more by the species of the trap plant than by the original tillage treatment of the field soil. The use of trap cultures for fungal diversity estimation in comparison with direct observation of field samples is discussed.

Fly Ash Soil Stabilization for Non-Uniform Subgrade Soils, Volume II: Influence of Subgrade Non-Uniformity on PCC Pavement Performance, TR-461, 2005

To provide insight into subgrade non-uniformity and its effects on pavement performance, this study investigated the influence of non-uniform subgrade support on pavement responses (stress and deflection) that affect pavement performance. Several reconstructed PCC pavement projects in Iowa were studied to document and evaluate the influence of subgrade/subbase non-uniformity on pavement performance. In situ field tests were performed at 12 sites to determine the subgrade/subbase engineering properties and develop a database of engineering parameter values for statistical and numerical analysis. Results of stiffness, moisture and density, strength, and soil classification were used to determine the spatial variability of a given property. Natural subgrade soils, fly ash-stabilized subgrade, reclaimed hydrated fly ash subbase, and granular subbase were studied. The influence of the spatial variability of subgrade/subbase on pavement performance was then evaluated by modeling the elastic properties of the pavement and subgrade using the ISLAB2000 finite element analysis program. A major conclusion from this study is that non-uniform subgrade/subbase stiffness increases localized deflections and causes principal stress concentrations in the pavement, which can lead to fatigue cracking and other types of pavement distresses. Field data show that hydrated fly ash, self-cementing fly ash-stabilized subgrade, and granular subbases exhibit lower variability than natural subgrade soils. Pavement life should be increased through the use of more uniform subgrade support. Subgrade/subbase construction in the future should consider uniformity as a key to long-term pavement performance.

DNA barcodes for soil animal taxonomy

The biodiversity of soil communities remains very poorly known and understood. Soil biological sciences are strongly affected by the taxonomic crisis, and most groups of animals in that biota suffer from a strong taxonomic impediment. The objective of this work was to investigate how DNA barcoding - a novel method using a microgenomic tag for species identification and discrimination - permits better evaluation of the taxonomy of soil biota. A total of 1,152 barcode sequences were analyzed for two major groups of animals, collembolans and earthworms, which presented broad taxonomic and geographic sampling. Besides strongly reflecting the taxonomic impediment for both groups, with a large number of species-level divergent lineages remaining unnamed so far, the results also highlight a high level (15%) of cryptic diversity within known species of both earthworms and collembolans. These results are supportive of recent local studies using a similar approach. Within an impeded taxonomic system for soil animals, DNA-assisted identification tools can facilitate and improve biodiversity exploration and description. DNA-barcoding campaigns are rapidly developing in soil animals and the community of soil biologists is urged to embrace these methods.

Prediction of soil orders with high spatial resolution: response of different classifiers to sampling density

The objective of this work was to evaluate sampling density on the prediction accuracy of soil orders, with high spatial resolution, in a viticultural zone of Serra Gaúcha, Southern Brazil. A digital elevation model (DEM), a cartographic base, a conventional soil map, and the Idrisi software were used. Seven predictor variables were calculated and read along with soil classes in randomly distributed points, with sampling densities of 0.5, 1, 1.5, 2, and 4 points per hectare. Data were used to train a decision tree (Gini) and three artificial neural networks: adaptive resonance theory, fuzzy ARTMap; self‑organizing map, SOM; and multi‑layer perceptron, MLP. Estimated maps were compared with the conventional soil map to calculate omission and commission errors, overall accuracy, and quantity and allocation disagreement. The decision tree was less sensitive to sampling density and had the highest accuracy and consistence. The SOM was the less sensitive and most consistent network. The MLP had a critical minimum and showed high inconsistency, whereas fuzzy ARTMap was more sensitive and less accurate. Results indicate that sampling densities used in conventional soil surveys can serve as a reference to predict soil orders in Serra Gaúcha.

Yield mapping, soil fertility and tree gaps in an orange orchard

The current high competition on Citrus industry demands from growers new management technologies for superior efficiency and sustainability. In this context, precision agriculture (PA) has developed techniques based on yield mapping and management systems that recognize field spatial variability, which contribute to increase profitability of commercial crops. Because spatial variability is often not perceived the orange orchards are still managed as uniform and adoption of PA technology on citrus farms is low. Thus, the objective of the present study was to characterize the spatial variability of three factors: fruit yield, soil fertility and occurrence of plant gaps caused by either citrus blight or huanglongbing (HLB) in a commercial Valencia orchard in Brotas, São Paulo State, Brazil. Data from volume, geographic coordinates and representative area of the bags used on harvest were recorded to generate yield points that were then interpolated to produce the yield map. Soil chemical characteristics were studied by analyzing samples collected along planting rows and inter-rows in 24 points distributed in the field. A map of density of tree gaps was produced by georeferencing individual gaps and later by counting the number of gaps within 500 m² cells. Data were submitted to statistical and geostatistical analyses. A t test was used to compare means of soil chemical characteristics between sampling regions. High variation on yield and density of tree gaps was observed from the maps. It was also demonstrated overlapping regions of high density of plant absence and low fruit yield. Soil fertility varied depending on the sampling region in the orchard. The spatial variability found on yield, soil fertility and on disease occurrence demonstrated the importance to adopt site specific nutrient management and disease control as tools to guarantee efficiency of fruit production.

Georeferenced database generation with the purpose of hydrologic molding in reservoirs of the hydrographic basin of Jaguaribe river in the state of Ceará, Brazil

The edafoclimatic conditions of the Brazilian semiarid region favor the water loss by surface runoff. The state of Ceará, almost completely covered by semiarid, has developed public policies for the construction of dams in order to attend the varied water demand. Several hydrological models were developed to support decisive processes in the complex management of reservoirs. This study aimed to establish a methodology for obtaining a georeferenced database suitable for use as input data in hydrological modeling in the semiarid of Ceará. It was used images of Landsat satellite and SRTM Mission, and soil maps of the state of Ceará. The Landsat images allowed the determination of the land cover and the SRTM Mission images, the automatic delineation of hydrographic basins. The soil type was obtained through the soil map. The database was obtained for Jaguaribe River hydrographic basin, in the state of Ceará, and is applicable to hydrological modeling based on the Curve Number method for estimating the surface runoff.

Mapping of use and occupation of the soil and irrigation water quality in the city of Salto do Lontra-Paraná, Brazil

The objective of this study consisted on mapping the use and soil occupation and evaluation of the quality of irrigation water used in Salto do Lontra, in the state of Paraná, Brazil. Images of the satellite SPOT-5 were used to perform the supervised classification of the Maximum Likelihood algorithm - MAXVER, and the water quality parameters analyzed were pH, EC, HCO3-, Cl-, PO4(3-), NO3-, turbidity, temperature and thermotolerant coliforms in two distinct rainfall periods. The water quality data were subjected to statistical analysis by the techniques of PCA and FA, to identify the most relevant variables in assessing the quality of irrigation water. The characterization of soil use and occupation by the classifier MAXVER allowed the identification of the following classes: crops, bare soil/stubble, forests and urban area. The PCA technique applied to irrigation water quality data explained 53.27% of the variation in water quality among the sampled points. Nitrate, thermotolerant coliforms, temperature, electrical conductivity and bicarbonate were the parameters that best explained the spatial variation of water quality.

Design of a Knowledge-Based Geographic Information System

A GIS has been designed with limited Functionalities; but with a novel approach in Aits design. The spatial data model adopted in the design of KBGIS is the unlinked vector model. Each map entity is encoded separately in vector fonn, without referencing any of its neighbouring entities. Spatial relations, in other words, are not encoded. This approach is adequate for routine analysis of geographic data represented on a planar map, and their display (Pages 105-106). Even though spatial relations are not encoded explicitly, they can be extracted through the specially designed queries. This work was undertaken as an experiment to study the feasibility of developing a GIS using a knowledge base in place of a relational database. The source of input spatial data was accurate sheet maps that were manually digitised. Each identifiable geographic primitive was represented as a distinct object, with its spatial properties and attributes defined. Composite spatial objects, made up of primitive objects, were formulated, based on production rules defining such compositions. The facts and rules were then organised into a production system, using OPS5

The Representative Soil Sampling Scheme of England and Wales: the spatial variation of topsoil nutrient status and pH between 1971 and 2001

The Representative Soil Sampling Scheme (RSSS) has monitored the soil of agricultural land in England and Wales since 1969. Here we describe the first spatial analysis of the data from these surveys using geostatistics. Four years of data (1971, 1981, 1991 and 2001) were chosen to examine the nutrient (available K, Mg and P) and pH status of the soil. At each farm, four fields were sampled; however, for the earlier years, coordinates were available for the farm only and not for each field. The averaged data for each farm were used for spatial analysis and the variograms showed spatial structure even with the smaller sample size. These variograms provide a reasonable summary of the larger scale of variation identified from the data of the more intensively sampled National Soil Inventory. Maps of kriged predictions of K generally show larger values in the central and southeastern areas (above 200 mg L-1) and an increase in values in the west over time, whereas Mg is fairly stable over time. The kriged predictions of P show a decline over time, particularly in the east, and those of pH show an increase in the east over time. Disjunctive kriging was used to examine temporal changes in available P using probabilities less than given thresholds of this element. The RSSS was not designed for spatial analysis, but the results show that the data from these surveys are suitable for this purpose. The results of the spatial analysis, together with those of the statistical analyses, provide a comprehensive view of the RSSS database as a basis for monitoring the soil. These data should be taken into account when future national soil monitoring schemes are designed.

Predicted soil organic carbon stocks and changes in the Brazilian Amazon between 2000 and 2030

Currently we have little understanding of the impacts of land use change on soil C stocks in the Brazilian Amazon. Such information is needed to determine impacts'6n the global C cycle and the sustainability of agricultural systems that are replacing native forest. The aim of this study was to predict soil carbon stocks and changes in the Brazilian Amazon during the period between 2000 and 2030, using the GEFSOC soil carbon (C) modelling system. In order to do so, we devised current and future land use scenarios for the Brazilian Amazon, taking into account: (i) deforestation, rates from the past three decades, (ii) census data on land use from 1940 to 2000, including the expansion and intensification of agriculture in the region, (iii) available information on management practices, primarily related to well managed pasture versus degraded pasture and conventional systems versus no-tillage systems for soybean (Glycine max) and (iv) FAO predictions on agricultural land use and land use changes for the years 2015 and 2030. The land use scenarios were integrated with spatially explicit soils data (SOTER database), climate, potential natural vegetation and land management units using the recently developed GEFSOC soil C modelling system. Results are presented in map, table and graph form for the entire Brazilian Amazon for the current situation (1990 and 2000) and the future (2015 and 2030). Results include soil organic C (SOC) stocks and SOC stock change rates estimated by three methods: (i) the Century ecosystem model, (ii) the Rothamsted C model and (iii) the intergovernmental panel on climate change (IPCC) method for assessing soil C at regional scale. In addition, we show estimated values of above and belowground biomass for native vegetation, pasture and soybean. The results on regional SOC stocks compare reasonably well with those based on mapping approaches. The GEFSOC system provided a means of efficiently handling complex interactions among biotic-edapho-climatic conditions (> 363,000 combinations) in a very large area (similar to 500 Mha) such as the Brazilian Amazon. All of the methods used showed a decline in SOC stock for the period studied; Century and RothC simulated values for 2030 being about 7% lower than those in 1990. Values from Century and RothC (30,430 and 25,000 Tg for the 0-20 cm layer for the Brazilian Amazon region were higher than those obtained from the IPCC system (23,400 Tg in the 0-30 cm layer). Finally; our results can help understand the major biogeochemical cycles that influence soil fertility and help devise management strategies that enhance the sustainability of these areas and thus slow further deforestation. (C) 2007 Elsevier B.V. All rights reserved.