Indices representing landscape as a support for pedological survey in GIS environment

Mapping of soil has been highlighted in the scientific community, because as alertness about the environment increases, it is necessary to understand more and more about the distribution of the soil in the landscape, as well as its potential and its limitations for the use. In that way the main aim of this study was to apply indices representing landscape with the use of geoprocessing to give support in the delimitation of different compartments of landscape. Primary indices used were altitude above channel network (AACN) and secondary channel network base level (CNBL), multiresolution index of valley bottom flatness (MRVBF) and Wetness index (ITW), having as object of study the Canguiri Experimental Farm, located in Pinhais, Curitiba's Metropolitan region. To correlate the chemical attributes and granulometric ones in sampling groups, totalizing 17 points (Sugamosto, 2002), a matrix of a simple linear correlation (Pearson) with the indices of the landscape were generated in the Software Statistica. The conclusion is that the indices representing the landscape used in the analysis of groupings were efficient as support to map soil at the level of suborder of Brazilian Soil Classification System.

Índices para a representação da paisagem como apoio para levantamento pedológico em ambiente de geoprocessamento

O mapeamento de solos tem ganhado destaque dentro da comunidade científica, pois, à medida em que a preocupação com o meio ambiente aumenta, ocorre a necessidade de se entender cada vez mais a distribuição dos solos na paisagem, tal como seu potencial e limitações de uso. Desta forma, o principal objetivo do trabalho foi aplicar índices de representação da paisagem com o apoio de geoprocessamento, suporte na delimitação dos diferentes compartimentos da paisagem. Foram utilizados índices primários Altitude above channel network (AACN) e secundários Channel network base level (CNBL, Multiresolution índex of valley bottom flatness (MRVBF) e Wetness índex (ITW), cujo estudo foi a Fazenda Experimental Canguiri, no município de Pinhais, região Metropolitana de Curitiba. Para correlacionar os atributos químicos e granulométricos por grupo de amostragem, totalizando 17 pontos (Sugamosto, 2002) foi gerada, no Software Statistica, uma matriz de correlação linear simples (Pearson) com os índices de representação da paisagem, os quais foram não apenas utilizados na análise de agrupamentos mas também eficientes no apoio do mapeamento dos solos, a nível de subordem do Sistema Brasileiro de Classificação de Solo.

Numerical analysis of the impact of charcoal production on soil hydrological behavior, runoff response and erosion susceptibility

The impact of charcoal production on soil hydraulic properties, runoff response and erosion susceptibility were studied in both field and simulation experiments. Core and composite samples, from 12 randomly selected sites within the catchment of Kotokosu were taken from the 0-10 cm layer of a charcoal site soil (CSS) and adjacent field soils (AFS). These samples were used to determine saturated hydraulic conductivity (Ksat), bulk density, total porosity, soil texture and color. Infiltration, surface albedo and soil surface temperature were also measured in both CSS and AFS. Measured properties were used as entries in a rainfall runoff simulation experiment on a smooth (5 % slope) plot of 25 x 25 m grids with 10 cm resolutions. Typical rainfall intensities of the study watershed (high, moderate and low) were applied to five different combinations of Ks distributions that could be expected in this landscape. The results showed significantly (p < 0.01) higher flow characteristics of the soil under charcoal kilns (increase of 88 %). Infiltration was enhanced and runoff volume reduced significantly. The results showed runoff reduction of about 37 and 18 %, and runoff coefficient ranging from 0.47-0.75 and 0.04-0.39 or simulation based on high (200 mm h-1) and moderate (100 mm h-1) rainfall events over the CSS and AFS areas, respectively. Other potential impacts of charcoal production on watershed hydrology were described. The results presented, together with watershed measurements, when available, are expected to enhance understanding of the hydrological responses of ecosystems to indiscriminate charcoal production and related activities in this region.

Numerical Dominance of the Argentine Ant vs Native Ants and Consequences on Soil Resource Searching in Mediterranean Cork-Oak Forests (Hymenoptera: Formicidae)

This study is focused on the dominance exerted by the invasive Argentine ant over native ants in a coastal Mediterranean area. Theimpact of this invasive ant on native ant assemblages and its consequences on total ant biomass and on the intensity of habitat explorationwere evaluated. Foraging ants were observed and their trajectories recorded during 5-minute periods in two study zones, one invaded andthe other non-invaded. Ant species detected, ant worker abundance, ant biomass and the intensity of soil surface searching done by antswere compared between the two zones. The Argentine ant invasion provoked a drastic reduction of the ant species richness. Apparentlyonly one native ant species is able to coexist with the Argentine ant, the cryptic Plagiolepis pygmaea. Ant worker abundance was also modified after the invasion: the number of Argentine ant workers detected, which represented 92% of the invaded zone, was two times higher than the number of native ant workers detected in the non-invaded zone. The total ant biomass was inversely affected, becoming four times lower in the invaded zone highly dominated by Linepithema humile. The higher number of Argentine ant workers and their fast tempo of activity implied an alteration of the intensity of soil surface searching: scanning by the Argentine ants in the invaded zone was higher than that done by the native ants in the non-invaded zone, and the estimated time for a complete soil surface scan was 64 minutes in the invaded zone and 108 minutes in the non-invaded zone. Consequently, resources will be discovered faster by ants in the invaded zone than in the non-invaded zone. The increase of the mean temperature and the decrease of the relative humidity from May to August reduced the ant activity in the two study zones but this reduction was greater in the invaded zone

Numerical Investigations on Soil Structure Interaction of Multistorey Frames

Frames are the most widely used structural system for multistorey buildings. A building frame is a three dimensional discrete structure consisting of a number of high rise bays in two directions at right angles to each other in the vertical plane. Multistorey frames are a three dimensional lattice structure which are statically indeterminate. Frames sustain gravity loads and resist lateral forces acting on it. India lies at the north westem end of the Indo-Australian tectonic plate and is identified as an active tectonic area. Under horizontal shaking of the ground, horizontal inertial forces are generated at the floor levels of a multistorey frame. These lateral inertia forces are transferred by the floor slab to the beams, subsequently to the columns and finally to the soil through the foundation system. There are many parameters that affect the response of a structure to ground excitations such as, shape, size and geometry of the structure, type of foundation, soil characteristics etc. The Soil Structure Interaction (SS1) effects refer to the influence of the supporting soil medium on the behavior of the structure when it is subjected to different types of loads. Interaction between the structure and its supporting foundation and soil, which is a complete system, has been modeled with finite elements. Numerical investigations have been carried out on a four bay, twelve storeyed regular multistorey frame considering depth of fixity at ground level, at characteristic depth of pile and at full depth. Soil structure interaction effects have been studied by considering two models for soil viz., discrete and continuum. Linear static analysis has been conducted to study the interaction effects under static load. Free vibration analysis and further shock spectrum analysis has been conducted to study the interaction effects under time dependent loads. The study has been extended to four types of soil viz., laterite, sand, alluvium and layered.The structural responses evaluated in the finite element analysis are bending moment, shear force and axial force for columns, and bending moment and shear force for beams. These responses increase with increase in the founding depth; however these responses show minimal increase beyond the characteristic length of pile. When the soil structure interaction effects are incorporated in the analysis, the aforesaid responses of the frame increases upto the characteristic depth and decreases when the frame has been analysed for the full depth. It has been observed that shock spectrum analysis gives wide variation of responses in the frame compared to linear elastic analysis. Both increase and decrease in responses have been observed in the interior storeys. The good congruence shown by the two finite element models viz., discrete and continuum in linear static analysis has been absent in shock spectrum analysis.

Delineation of specific management areas for coffee cultivation based on the soil-relief relationship and numerical classification

Predicting and mapping productivity areas allows crop producers to improve their planning of agricultural activities. The primary aims of this work were the identification and mapping of specific management areas allowing coffee bean quality to be predicted from soil attributes and their relationships to relief. The study area was located in the Southeast of the Minas Gerais state, Brazil. A grid containing a total of 145 uniformly spaced nodes 50 m apart was established over an area of 31. 7 ha from which samples were collected at depths of 0. 00-0. 20 m in order to determine physical and chemical attributes of the soil. These data were analysed in conjunction with plant attributes including production, proportion of beans retained by different sieves and drink quality. The results of principal component analysis (PCA) in combination with geostatistical data showed the attributes clay content and available iron to be the best choices for identifying four crop production environments. Environment A, which exhibited high clay and available iron contents, and low pH and base saturation, was that providing the highest yield (30. 4l ha-1) and best coffee beverage quality (61 sacks ha-1). Based on the results, we believe that multivariate analysis, geostatistics and the soil-relief relationships contained in the digital elevation model (DEM) can be effectively used in combination for the hybrid mapping of areas of varying suitability for coffee production. © 2012 Springer Science+Business Media New York.

Soil-structure interaction during tunnelling in urban area: observations and 3D numerical modelling

This work illustrates a soil-tunnel-structure interaction study performed by an integrated,geotechnical and structural,approach based on 3D finite element analyses and validated against experimental observations.The study aims at analysing the response of reinforced concrete framed buildings on discrete foundations in interaction with metro lines.It refers to the case of the twin tunnels of the Milan (Italy) metro line 5,recently built in coarse grained materials using EPB machines,for which subsidence measurements collected along ground and building sections during tunnelling were available.Settlements measured under freefield conditions are firstly back interpreted using Gaussian empirical predictions. Then,the in situ measurements’ analysis is extended to include the evolving response of a 9 storey reinforced concrete building while being undercrossed by the metro line.In the finite element study,the soil mechanical behaviour is described using an advanced constitutive model. This latter,when combined with a proper simulation of the excavation process, proves to realistically reproduce the subsidence profiles under free field conditions and to capture the interaction phenomena occurring between the twin tunnels during the excavation. Furthermore, when the numerical model is extended to include the building, schematised in a detailed manner, the results are in good agreement with the monitoring data for different stages of the twin tunnelling. Thus, they indirectly confirm the satisfactory performance of the adopted numerical approach which also allows a direct evaluation of the structural response as an outcome of the analysis. Further analyses are also carried out modelling the building with different levels of detail. The results highlight that, in this case, the simplified approach based on the equivalent plate schematisation is inadequate to capture the real tunnelling induced displacement field. The overall behaviour of the system proves to be mainly influenced by the buried portion of the building which plays an essential role in the interaction mechanism, due to its high stiffness.

Long-term behaviour of railway transitions under dynamic loading application to soft soil sites

Dissertação para obtenção do Grau de Doutor em Engenharia Civil

Desarrollo sustentable de geoestructuras compuestas por suelo loéssico compactado mejorado con agregado de materiales estabilizantes. Sustainable development of compacted loess soil geo-structures with addition of stabilizing materials.

Los suelos estabilizados mediante compactación, permiten obtener materiales con ventajas ténicas y economicas en diferentes tipos de obras de ingeniería. Ejemplos de su uso se tiene en bases viales de autopistas, rutas o calles urbanas, pistas de aterrizaje, barreras de contención para enterramientos sanitarios o lagunas de estabilización, apoyos de plateas para fundación de edificios, losas industriales, entre otras aplicaciones. Las fallas en este tipo de construcciones pueden resultar en catástrofes ambientales, sociales y elevadas pérdidas económicas, por lo que resulta de gran importancia optimizar el diseño e incrementar la seguridad de este tipo de construcciones. Las obras con estas características involucran grandes volúmenes y/o superficies que requieren controles sistemáticos durante su desarrollo, a los fines de garantizar el cumplimiento de las propiedades de los materiales establecidos en la etapa de diseño. De esta forma, es necesario contar con ensayos de campo sencillos, confiables y eficientes que permitan identificar propiedades físicas, mecánicas e hidráulicas. Las geoestructuras generadas mediante la compactación del suelo próximo al sector de construcción pueden funcionar adecuadamente, con reducidos costos de material y transporte. Su estabilización puede ejecutarse en forma natural, o con la incorporación de agregados minerales como bentonita, cal o cemento. Estas incorporaciones mejoran las propiedades hidráulicas y mecánicas del material, optimizando el comportamiento requerido para la obra. Para establecer la forma en la que estos minerales modifican el comportamiento del suelo local compactado deben realizarse investigaciones especiales con los materiales involucrados. En el ámbito internacional existen numerosas investigaciones sobre comportamiento de suelos compactados, no obstante, si bien aportan antecedentes para la planificación de estudios locales, sus resultados no pueden trasladarse de manera directa. Las características propias del suelo local constituye la principal variable debido a la diversidad en las propiedades geotécnicas de cada Región. Esta investigación, se focaliza en el empleo de suelos limosos de la formación loéssica de la zona central de Argentina. Los suelos de la llanura cordobesa poseen comportamientos particulares, los cuales son contemplados en los diseños presentados como resutado de las investigaciones internacionales. Esta particularidad se relaciona con su inestabilidad, lo que los clasifica como suelos colapsables. Los resultados obtenidos en este trabajo podrán ser extendidos a una gran superficie de la Provincia de Córdoba y a la Región Pampeana en general, a los fines de establecer recomendaciones de diseño y construcción para la confección de Pliegos de Especificaciones Técnicas de diferentes tipos de obras públicas y privadas. El estudio contempla la ejecución de un plan experimental a escala de laboratorio y campo. Los materiales corresponden a suelo limosos puros, y diferentes agregados tales como bentonita, cal y cemento. Se planifican ensayos para evaluar el desempeño del material, a partir de la confección de muestras preparadas con diferentes condiciones de compactación (energía, humedad y método), y en forma de mezcla con los distintos tipos de agregados. Se realizarán ensayos de permeabilidad en celdas de pared rígida y flexible, junto a ensayos mecánicos de compresión confinada, simple y triaxial. Para el trabajo experimental de campo se prevé la ejecución de terraplenes de prueba instrumentados con tensiómetros e infiltrómetros para evaluar el comportamiento hidraúlico en el tiempo, junto con ensayos de penetración y plato de carga para la caracterización mecánica. En forma conjunta se propone el desarrollo de modelos numéricos de caracterización hidromecánica. Stabilized soils by compaction, produce materials technical and economic advantages in different types of engineering works. For example, road bases in highways, roads or city streets, containment barriers for sanitary landfill or stabilization ponds, foundation support of building, industrial flat, and other applications. Failures can result in environmental catastrophes, social, and economic loss, so it is important to optimize the design and increase the safety of such buildings. These works involve large surfaces that require systematic tests during construction, so it is necessary to have simple field tests, reliable and efficient to identify physical, mechanical and hydraulic properties. The geo-structures generated by local soil compaction have reduced material and transportation costs. Stabilization can be naturally, or with the addition of mineral aggregates as bentonite, lime and cement. These additions improve the hydraulic and mechanical properties of the material. So, special investigations should be conducted with the materials involved. There are many international studies on compacted soils behavior but their results can not be transferred directly due to the particularities of regional soils. For this research silty soils of central Argentina are the main focus. The soils of Córdoba plains are instability, so are classified as collapsible soils. The results obtained in this work may be extended to a large area of the Province of Cordoba and the Pampas region in general, in order to establish design and construction recommendations. The study includes laboratory and field tests. The materials are pure silty soil, and different aggregates such as bentonite, lime and cement. Tests are planned to evaluate the performance. Laboratory includes rigid and flexible wall cells, confined, triaxial and simple compression tests. For field experimental instrumented embankments will be constructed. A numerical hydromechanical model will be developed.

3D geoelectrical monitoring to follow hydric soil behaviour during rainfall

RESUME Durant les dernières années, les méthodes électriques ont souvent été utilisées pour l'investigation des structures de subsurface. L'imagerie électrique (Electrical Resistivity Tomography, ERT) est une technique de prospection non-invasive et spatialement intégrée. La méthode ERT a subi des améliorations significatives avec le développement de nouveaux algorithmes d'inversion et le perfectionnement des techniques d'acquisition. La technologie multicanale et les ordinateurs de dernière génération permettent la collecte et le traitement de données en quelques heures. Les domaines d'application sont nombreux et divers: géologie et hydrogéologie, génie civil et géotechnique, archéologie et études environnementales. En particulier, les méthodes électriques sont souvent employées dans l'étude hydrologique de la zone vadose. Le but de ce travail est le développement d'un système de monitorage 3D automatique, non- invasif, fiable, peu coûteux, basé sur une technique multicanale et approprié pour suivre les variations de résistivité électrique dans le sous-sol lors d'événements pluvieux. En raison des limitations techniques et afin d'éviter toute perturbation physique dans la subsurface, ce dispositif de mesure emploie une installation non-conventionnelle, où toutes les électrodes de courant sont placées au bord de la zone d'étude. Le dispositif le plus approprié pour suivre les variations verticales et latérales de la résistivité électrique à partir d'une installation permanente a été choisi à l'aide de modélisations numériques. Les résultats démontrent que le dispositif pôle-dipôle offre une meilleure résolution que le dispositif pôle-pôle et plus apte à détecter les variations latérales et verticales de la résistivité électrique, et cela malgré la configuration non-conventionnelle des électrodes. Pour tester l'efficacité du système proposé, des données de terrain ont été collectées sur un site d'étude expérimental. La technique de monitorage utilisée permet de suivre le processus d'infiltration 3D pendant des événements pluvieux. Une bonne corrélation est observée entre les résultats de modélisation numérique et les données de terrain, confirmant par ailleurs que le dispositif pôle-dipôle offre une meilleure résolution que le dispositif pôle-pôle. La nouvelle technique de monitorage 3D de résistivité électrique permet de caractériser les zones d'écoulement préférentiel et de caractériser le rôle de la lithologie et de la pédologie de manière quantitative dans les processus hydrologiques responsables d'écoulement de crue. ABSTRACT During the last years, electrical methods were often used for the investigation of subsurface structures. Electrical resistivity tomography (ERT) has been reported to be a useful non-invasive and spatially integrative prospecting technique. The ERT method provides significant improvements, with the developments of new inversion algorithms, and the increasing efficiency of data collection techniques. Multichannel technology and powerful computers allow collecting and processing resistivity data within few hours. Application domains are numerous and varied: geology and hydrogeology, civil engineering and geotechnics, archaeology and environmental studies. In particular, electrical methods are commonly used in hydrological studies of the vadose zone. The aim of this study was to develop a multichannel, automatic, non-invasive, reliable and inexpensive 3D monitoring system designed to follow electrical resistivity variations in soil during rainfall. Because of technical limitations and in order to not disturb the subsurface, the proposed measurement device uses a non-conventional electrode set-up, where all the current electrodes are located near the edges of the survey grid. Using numerical modelling, the most appropriate arrays were selected to detect vertical and lateral variations of the electrical resistivity in the framework of a permanent surveying installation system. The results show that a pole-dipole array has a better resolution than a pole-pole array and can successfully follow vertical and lateral resistivity variations despite the non-conventional electrode configuration used. Field data are then collected at a test site to assess the efficiency of the proposed monitoring technique. The system allows following the 3D infiltration processes during a rainfall event. A good correlation between the results of numerical modelling and field data results can be observed since the field pole-dipole data give a better resolution image than the pole-pole data. The new device and technique makes it possible to better characterize the zones of preferential flow and to quantify the role of lithology and pedology in flood- generating hydrological processes.

Numerical modeling of the effect of variation of boundary conditions on vadose zone hydraulic properties

An accurate estimation of hydraulic fluxes in the vadose zone is essential for the prediction of water, nutrient and contaminant transport in natural systems. The objective of this study was to simulate the effect of variation of boundary conditions on the estimation of hydraulic properties (i.e. water content, effective unsaturated hydraulic conductivity and hydraulic flux) in a one-dimensional unsaturated flow model domain. Unsaturated one-dimensional vertical water flow was simulated in a pure phase clay loam profile and in clay loam interlayered with silt loam distributed according to the third iteration of the Cantor Bar fractal object Simulations were performed using the numerical model Hydrus 1D. The upper and lower pressure heads were varied around average values of -55 cm for the near-saturation range. This resulted in combinations for the upper and lower constant head boundary conditions, respectively, of -50 and -60 cm, -40 and -70 cm, -30 and -80 cm, -20 and -90 cm, and -10 and -100 cm. For the drier range the average head between the upper and lower boundary conditions was set to -550 cm, resulting in the combinations -500 and -600 cm, -400 and -700 cm, -300 and -800 cm, -200 and -900 cm, and -100 and -1,000 cm, for upper and lower boundary conditions, respectively. There was an increase in water contents, fluxes and hydraulic conductivities with the increase in head difference between boundary conditions. Variation in boundary conditions in the pure phase and interlayered one-dimensional profiles caused significant deviations in fluxes, water contents and hydraulic conductivities compared to the simplest case (a head difference between the upper and lower constant head boundaries of 10 cm in the wetter range and 100 cm in the drier range).

Using numerical classification of profiles based on Vis-NIR spectra to distinguish soils from the Piracicaba Region, Brazil

Considering that information from soil reflectance spectra is underutilized in soil classification, this paper aimed to evaluate the relationship of soil physical, chemical properties and their spectra, to identify spectral patterns for soil classes, evaluate the use of numerical classification of profiles combined with spectral data for soil classification. We studied 20 soil profiles from the municipality of Piracicaba, State of São Paulo, Brazil, which were morphologically described and classified up to the 3rd category level of the Brazilian Soil Classification System (SiBCS). Subsequently, soil samples were collected from pedogenetic horizons and subjected to soil particle size and chemical analyses. Their Vis-NIR spectra were measured, followed by principal component analysis. Pearson's linear correlation coefficients were determined among the four principal components and the following soil properties: pH, organic matter, P, K, Ca, Mg, Al, CEC, base saturation, and Al saturation. We also carried out interpretation of the first three principal components and their relationships with soil classes defined by SiBCS. In addition, numerical classification of the profiles based on the OSACA algorithm was performed using spectral data as a basis. We determined the Normalized Mutual Information (NMI) and Uncertainty Coefficient (U). These coefficients represent the similarity between the numerical classification and the soil classes from SiBCS. Pearson's correlation coefficients were significant for the principal components when compared to sand, clay, Al content and soil color. Visual analysis of the principal component scores showed differences in the spectral behavior of the soil classes, mainly among Argissolos and the others soils. The NMI and U similarity coefficients showed values of 0.74 and 0.64, respectively, suggesting good similarity between the numerical and SiBCS classes. For example, numerical classification correctly distinguished Argissolos from Latossolos and Nitossolos. However, this mathematical technique was not able to distinguish Latossolos from Nitossolos Vermelho férricos, but the Cambissolos were well differentiated from other soil classes. The numerical technique proved to be effective and applicable to the soil classification process.

Numerical procedure for estimating temperature in solarized soils

The objective of this work was to develop a simplified numerical procedure for the estimation of accumulated monthly hours of solarized soil temperatures. The proposed model requires monthly means of daily solar radiation and maximum air temperature as input data, and a daily pattern of temperature variation assumed to be sine-shaped. The procedure was verified using observations made during the years 1992 and 1993 in Jaguariúna, SP. The proposed procedure can predict monthly temperature hours at 10 cm depth in the solarized soil, with acceptable accuracy, in the region for which it was developed.

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.

Application of x-ray computed tomography in the evaluation of soil porosity in soil management systems

The study aimed to evaluate a methodology to quantify the porosity of the soil using computed tomography in areas under no-tillage, conventional tillage and native forest. Three soil management systems were selected for the study: forest, conventional tillage and no-tillage. In each soil management system, undisturbed soil samples were collected in the surface layer (0.0 to 0.10 m). The tomographic images were obtained using a X-ray microtomography. After obtaining the images, they were processed, and a methodology was evaluated for image conversion into numerical values. The statistical method which provided the greatest accuracy was the percentile method. The methodology used to analyze the tomographic image allowed quantifying the porosity of the soil under different soil management. The method enabled the characterization of soil porosity in a non-evasive and non-destructive way.