Shadow Analysis of Soil Surface Roughness Compared to the Chain Set Method and Direct Measurement of Micro-relief.

Erosion potential and the effects of tillage can be evaluated from quantitative descriptions of soil surface roughness. The present study therefore aimed to fill the need for a reliable, low-cost and convenient method to measure that parameter. Based on the interpretation of micro-topographic shadows, this new procedure is primarily designed for use in the field after tillage. The principle underlying shadow analysis is the direct relationship between soil surface roughness and the shadows cast by soil structures under fixed sunlight conditions. The results obtained with this method were compared to the statistical indexes used to interpret field readings recorded by a pin meter. The tests were conducted on 4-m2 sandy loam and sandy clay loam plots divided into 1-m2 subplots tilled with three different tools: chisel, tiller and roller. The highly significant correlation between the statistical indexes and shadow analysis results obtained in the laboratory as well as in the field for all the soil?tool combinations proved that both variability (CV) and dispersion (SD) are accommodated by the new method. This procedure simplifies the interpretation of soil surface roughness and shortens the time involved in field operations by a factor ranging from 12 to 20.

Influence of tillage on soil surface roughness structure andsurface porosity.

Soil erosion is a complex phenomenon involving the detachment and transport of soil particles, storage and runoff of rainwater, and infiltration. The relative magnitude and importance of these processes depends on several factors being one of them surface microtopography, usually quantified trough soil surface roughness (SSR). Surface soil porosity and SSR can be altered by tillage operation. Even though the surface porosity is an important parameter of a tilled field, however, no practical technique for rapid and non-contact measurement of surface porosity has been developed yet.

Fractal analysis of laplacian pyramidal filters applied to segmentation of soil images

The laplacian pyramid is a well-known technique for image processing in which local operators of many scales, but identical shape, serve as the basis functions. The required properties to the pyramidal filter produce a family of filters, which is unipara metrical in the case of the classical problem, when the length of the filter is 5. We pay attention to gaussian and fractal behaviour of these basis functions (or filters), and we determine the gaussian and fractal ranges in the case of single parameter ?. These fractal filters loose less energy in every step of the laplacian pyramid, and we apply this property to get threshold values for segmenting soil images, and then evaluate their porosity. Also, we evaluate our results by comparing them with the Otsu algorithm threshold values, and conclude that our algorithm produce reliable test results.

Computer simulation of random parckings for self-similar particle size distributions in soil and granular materials: porosity and pore size distribution

A 2D computer simulation method of random packings is applied to sets of particles generated by a self-similar uniparametric model for particle size distributions (PSDs) in granular media. The parameter p which controls the model is the proportion of mass of particles corresponding to the left half of the normalized size interval [0,1]. First the influence on the total porosity of the parameter p is analyzed and interpreted. It is shown that such parameter, and the fractal exponent of the associated power scaling, are efficient packing parameters, but this last one is not in the way predicted in a former published work addressing an analogous research in artificial granular materials. The total porosity reaches the minimum value for p = 0.6. Limited information on the pore size distribution is obtained from the packing simulations and by means of morphological analysis methods. Results show that the range of pore sizes increases for decreasing values of p showing also different shape in the volume pore size distribution. Further research including simulations with a greater number of particles and image resolution are required to obtain finer results on the hierarchical structure of pore space.

Real-Time Dynamic PGD Calculation of Non-Linear Soil Behavior

El estudio sísmico en los últimos 50 años y el análisis del comportamiento dinámico del suelo revelan que el comportamiento del suelo es altamente no lineal e histéretico incluso para pequeñas deformaciones. El comportamiento no lineal del suelo durante un evento sísmico tiene un papel predominante en el análisis de la respuesta de sitio. Los análisis unidimensionales de la respuesta sísmica del suelo son a menudo realizados utilizando procedimientos lineales equivalentes, que requieren generalmente pocos parámetros conocidos. Los análisis de respuesta de sitio no lineal tienen el potencial para simular con mayor precisión el comportamiento del suelo, pero su aplicación en la práctica se ha visto limitada debido a la selección de parámetros poco documentadas y poco claras, así como una inadecuada documentación de los beneficios del modelado no lineal en relación al modelado lineal equivalente. En el análisis del suelo, el comportamiento del suelo es aproximado como un sólido Kelvin-Voigt con un módulo de corte elástico y amortiguamiento viscoso. En el análisis lineal y no lineal del suelo se están considerando geometrías y modelos reológicos más complejos. El primero está siendo dirigido por considerar parametrizaciones más ricas del comportamiento linealizado y el segundo mediante el uso de multi-modo de los elementos de resorte-amortiguador con un eventual amortiguador fraccional. El uso del cálculo fraccional está motivado en gran parte por el hecho de que se requieren menos parámetros para lograr la aproximación exacta a los datos experimentales. Basándose en el modelo de Kelvin-Voigt, la viscoelasticidad es revisada desde su formulación más estándar a algunas descripciones más avanzada que implica la amortiguación dependiente de la frecuencia (o viscosidad), analizando los efectos de considerar derivados fraccionarios para representar esas contribuciones viscosas. Vamos a demostrar que tal elección se traduce en modelos más ricos que pueden adaptarse a diferentes limitaciones relacionadas con la potencia disipada, amplitud de la respuesta y el ángulo de fase. Por otra parte, el uso de derivados fraccionarios permite acomodar en paralelo, dentro de un análogo de Kelvin-Voigt generalizado, muchos amortiguadores que contribuyen a aumentar la flexibilidad del modelado para la descripción de los resultados experimentales. Obviamente estos modelos ricos implican muchos parámetros, los asociados con el comportamiento y los relacionados con los derivados fraccionarios. El análisis paramétrico de estos modelos requiere técnicas numéricas eficientemente capaces de simular comportamientos complejos. El método de la Descomposición Propia Generalizada (PGD) es el candidato perfecto para la construcción de este tipo de soluciones paramétricas. Podemos calcular off-line la solución paramétrica para el depósito de suelo, para todos los parámetros del modelo, tan pronto como tales soluciones paramétricas están disponibles, el problema puede ser resuelto en tiempo real, porque no se necesita ningún nuevo cálculo, el solucionador sólo necesita particularizar on-line la solución paramétrica calculada off-line, que aliviará significativamente el procedimiento de solución. En el marco de la PGD, parámetros de los materiales y los diferentes poderes de derivación podrían introducirse como extra-coordenadas en el procedimiento de solución. El cálculo fraccional y el nuevo método de reducción modelo llamado Descomposición Propia Generalizada han sido aplicado en esta tesis tanto al análisis lineal como al análisis no lineal de la respuesta del suelo utilizando un método lineal equivalente. ABSTRACT Studies of earthquakes over the last 50 years and the examination of dynamic soil behavior reveal that soil behavior is highly nonlinear and hysteretic even at small strains. Nonlinear behavior of soils during a seismic event has a predominant role in current site response analysis. One-dimensional seismic ground response analysis are often performed using equivalent-linear procedures, which require few, generally well-known parameters. Nonlinear analyses have the potential to more accurately simulate soil behavior, but their implementation in practice has been limited because of poorly documented and unclear parameter selection, as well as inadequate documentation of the benefits of nonlinear modeling relative to equivalent linear modeling. In soil analysis, soil behaviour is approximated as a Kelvin-Voigt solid with a elastic shear modulus and viscous damping. In linear and nonlinear analysis more complex geometries and more complex rheological models are being considered. The first is being addressed by considering richer parametrizations of the linearized behavior and the second by using multi-mode spring-dashpot elements with eventual fractional damping. The use of fractional calculus is motivated in large part by the fact that fewer parameters are required to achieve accurate approximation of experimental data. Based in Kelvin-Voigt model the viscoelastodynamics is revisited from its most standard formulation to some more advanced description involving frequency-dependent damping (or viscosity), analyzing the effects of considering fractional derivatives for representing such viscous contributions. We will prove that such a choice results in richer models that can accommodate different constraints related to the dissipated power, response amplitude and phase angle. Moreover, the use of fractional derivatives allows to accommodate in parallel, within a generalized Kelvin-Voigt analog, many dashpots that contribute to increase the modeling flexibility for describing experimental findings. Obviously these rich models involve many parameters, the ones associated with the behavior and the ones related to the fractional derivatives. The parametric analysis of all these models require efficient numerical techniques able to simulate complex behaviors. The Proper Generalized Decomposition (PGD) is the perfect candidate for producing such kind of parametric solutions. We can compute off-line the parametric solution for the soil deposit, for all parameter of the model, as soon as such parametric solutions are available, the problem can be solved in real time because no new calculation is needed, the solver only needs particularize on-line the parametric solution calculated off-line, which will alleviate significantly the solution procedure. Within the PGD framework material parameters and the different derivation powers could be introduced as extra-coordinates in the solution procedure. Fractional calculus and the new model reduction method called Proper Generalized Decomposition has been applied in this thesis to the linear analysis and nonlinear soil response analysis using a equivalent linear method.

Dispersivity Determination Through a Modeling Approach From a Tracer Test Based on Total Br Concentration in Soil Samples

Determination of reliable solute transport parameters is an essential aspect for the characterization of the mechanisms and processes involved in solute transport (e.g., pesticides, fertilizers, contaminants) through the unsaturated zone. A rapid inexpensive method to estimate the dispersivity parameter at the field scale is presented herein. It is based on the quantification by the X-ray fluorescence solid-state technique of total bromine in soil, along with an inverse numerical modeling approach. The results show that this methodology is a good alternative to the classic Br− determination in soil water by ion chromatography. A good agreement between the observed and simulated total soil Br is reported. The results highlight the potential applicability of both combined techniques to infer readily solute transport parameters under field conditions.

(Table 1) Soil physical parameters of sediments obtained in the Lüchow-Dannenberg district, Lower Saxony

The mobility of the radionuclides of the elements Sr, I, Cs and Ce were investigated for three typical sands of Northern Germany under simulated natural, undersaturated flow conditions. The laboratory experiments include the determination of the flow parameters of the seepwater movement as well as the transport velocities of the radionuclides in the sands. For the three selected sands, the following mobility sequence/order has been established for the radionuclides: I < Sr < Cs < Ce

Chloride as a signature indicator of soil textural and hydrologic stratigraphies in variable charge deep profiles

Soil properties that influence water movement through profiles are important for determining flow paths, reactions between soil and solute, and the ultimate destination of solutes. This is particularly important in high rainfall environments. For highly weathered deep profiles, we hypothesize that abrupt changes in the distribution of the quotient [QT = (silt + sand)/clay] reflect the boundaries between textural units or textural (TS) and hydrologic (HS) stratigraphies. As a result, QT can be used as a parameter to characterize TS and as a surrogate for HS. Secondly, we propose that if chloride distributions were correlated with QT, under non-limiting anion exchange, then chloride distributions can be used as a signature indicator of TS and HS. Soil cores to a depth of 12.5 in were taken from 16 locations in the wet tropical Johnstone River catchment of northeast Queensland, Australia. The cores belong to nine variable charge soil types and were under sugarcane (Saccharun officinarum-S) production, which included the use of potassium chloride, for several decades. The cores were segmented at I m depth increments and subsamples were analysed for chloride, pH, soil water content (theta), clay, silt and sand contents. Selected bores were capped to serve as piezometers to monitor groundwater dynamics. Depth incremented QT, theta and chloride correlated, each individually, significantly with the corresponding profile depth increments, indicating the presence of textural, hydrologic and chloride gradients in profiles. However, rapid increases in QT down the profile indicated abrupt changes in TS, suggesting that QT can be used as a parameter to characterize TS and as a surrogate for HS. Abrupt changes in chloride distributions were similar to QT, suggesting that chloride distributions can be used as a signature indicator of QT (TS) and HS. Groundwater data indicated that chloride distributions depended, at least partially, on groundwater dynamics, providing further support to our hypothesis that chloride distribution can be used as a signature indicator of HS. Copyright (c) 2005 John Wiley & Sons, Ltd.

Circumpolar dataset of Soil Organic Carbon north of treeline derived from ENVISAT ASAR GM, link to GeoTIFF

A new approach for the estimation of soil organic carbon (SOC) pools north of the tree line has been developed based on synthetic aperture radar (SAR; ENVISAT Advanced SAR Global Monitoring mode) data. SOC values are directly determined from backscatter values instead of upscaling using land cover or soil classes. The multi-mode capability of SAR allows application across scales. It can be shown that measurements in C band under frozen conditions represent vegetation and surface structure properties which relate to soil properties, specifically SOC. It is estimated that at least 29 Pg C is stored in the upper 30 cm of soils north of the tree line. This is approximately 25 % less than stocks derived from the soil-map-based Northern Circumpolar Soil Carbon Database (NCSCD). The total stored carbon is underestimated since the established empirical relationship is not valid for peatlands or strongly cryoturbated soils. The approach does, however, provide the first spatially consistent account of soil organic carbon across the Arctic. Furthermore, it could be shown that values obtained from 1 km resolution SAR correspond to accounts based on a high spatial resolution (2 m) land cover map over a study area of about 7 × 7 km in NE Siberia. The approach can be also potentially transferred to medium-resolution C-band SAR data such as ENVISAT ASAR Wide Swath with ~120 m resolution but it is in general limited to regions without woody vegetation. Global Monitoring-mode-derived SOC increases with unfrozen period length. This indicates the importance of this parameter for modelling of the spatial distribution of soil organic carbon storage.

(Tables 2-7) Soil characteristics and geochemistry of marshland in the North Frisian mainland, Schleswig-Holstein, Germany

Results from a large scale soil mapping on the North Frisian mainland indicate, that field characteristics, particularly the grain-size, bedding, and degree of compaction, with in general determine the soil units mapped, are closely correlated with each other and with other field and laboratory data. Exchangable ions and the Ca/Mg-ratio, however, indicate no explainable connections with the soil units and with most of the other field characteristics but are determined postsedimentarily by processes of the development of soil and landscape, such as desalting and decalcification, silicate weathering, fresh- and salt-water innundations, salty precipitations, salty groundwater and fertilization. Therefore the Ca/Mg-ratio is not suitable to differentiate between more clayey compacted Knick-marsh soils and less clayey permeable Klei-marsh soils. The results confirm that marsh-soils may only be classified and mapped by means of all available field-data which have to be supplemented by laboratory investigations.

Tab. 1+2+3: Selected principal soil properties in the oil exploitation region of KomiArcticOil (Usinsk) in the Russian tundra

Oil polluted and not oil polluted soils (crude oil hydrocarbons contents: 20-92500 mg/kg dry soil mass) under natural grass and forest vegetation and in a bog in the Russian tundra were compared in their principal soil ecological parameters, the oil content and the microbial indicators. CFE biomass-C, dehydrogenase and arylsulfatase activity were enhanced with the occurrence of crude oil. Using these parameters for purposes of controlling remediation and recultivation success it is not possible to distinguish bctween promotion of microbial activity by oil carbon or soil organic carbon (SOC). For this reason we think that these parameters are not appropriate to indicate a soil damage by an oil impact. In contrast the metabolie quotient (qC02), calculated as the ratio between soil basal respiration and the SIR biomass-C was adequate to indicate a high crude oil contamination in soil. Also, the ß-glucosidase activity (parameter ß-GL/SOC) was correlated negatively with oil in soil. The indication of a soil damage by using the stress parameter qCO, or the specific enzyme activities (activity/SOC) minimizes the promotion effect of the recent SOC content on microbial parameters. Both biomass methods (SIR, CFE) have technical problems in application for crude oil-contaminated and subarctic soils. CFE does not reflect the low C_mic level of the cold tundra soils. We recommend to test every method for its suitability before any data collection in series as well as application for cold soils and the application of ecophysiological ratios as R_mic/C_mic, C_mic/SOC or enzymatic activity/SOC instead of absolute data.

Soft computing based parameter identification in pavements and geomechanical systems

Accurate estimation of road pavement geometry and layer material properties through the use of proper nondestructive testing and sensor technologies is essential for evaluating pavement’s structural condition and determining options for maintenance and rehabilitation. For these purposes, pavement deflection basins produced by the nondestructive Falling Weight Deflectometer (FWD) test data are commonly used. The nondestructive FWD test drops weights on the pavement to simulate traffic loads and measures the created pavement deflection basins. Backcalculation of pavement geometry and layer properties using FWD deflections is a difficult inverse problem, and the solution with conventional mathematical methods is often challenging due to the ill-posed nature of the problem. In this dissertation, a hybrid algorithm was developed to seek robust and fast solutions to this inverse problem. The algorithm is based on soft computing techniques, mainly Artificial Neural Networks (ANNs) and Genetic Algorithms (GAs) as well as the use of numerical analysis techniques to properly simulate the geomechanical system. A widely used pavement layered analysis program ILLI-PAVE was employed in the analyses of flexible pavements of various pavement types; including full-depth asphalt and conventional flexible pavements, were built on either lime stabilized soils or untreated subgrade. Nonlinear properties of the subgrade soil and the base course aggregate as transportation geomaterials were also considered. A computer program, Soft Computing Based System Identifier or SOFTSYS, was developed. In SOFTSYS, ANNs were used as surrogate models to provide faster solutions of the nonlinear finite element program ILLI-PAVE. The deflections obtained from FWD tests in the field were matched with the predictions obtained from the numerical simulations to develop SOFTSYS models. The solution to the inverse problem for multi-layered pavements is computationally hard to achieve and is often not feasible due to field variability and quality of the collected data. The primary difficulty in the analysis arises from the substantial increase in the degree of non-uniqueness of the mapping from the pavement layer parameters to the FWD deflections. The insensitivity of some layer properties lowered SOFTSYS model performances. Still, SOFTSYS models were shown to work effectively with the synthetic data obtained from ILLI-PAVE finite element solutions. In general, SOFTSYS solutions very closely matched the ILLI-PAVE mechanistic pavement analysis results. For SOFTSYS validation, field collected FWD data were successfully used to predict pavement layer thicknesses and layer moduli of in-service flexible pavements. Some of the very promising SOFTSYS results indicated average absolute errors on the order of 2%, 7%, and 4% for the Hot Mix Asphalt (HMA) thickness estimation of full-depth asphalt pavements, full-depth pavements on lime stabilized soils and conventional flexible pavements, respectively. The field validations of SOFTSYS data also produced meaningful results. The thickness data obtained from Ground Penetrating Radar testing matched reasonably well with predictions from SOFTSYS models. The differences observed in the HMA and lime stabilized soil layer thicknesses observed were attributed to deflection data variability from FWD tests. The backcalculated asphalt concrete layer thickness results matched better in the case of full-depth asphalt flexible pavements built on lime stabilized soils compared to conventional flexible pavements. Overall, SOFTSYS was capable of producing reliable thickness estimates despite the variability of field constructed asphalt layer thicknesses.

Effect of nitrogen fertilization on soil microbial biomass in an Oxisol cultivated with irrigated barley in the Brazilian Cerrado.

Abstract:The aim of this study was to evaluate the effect of different nitrogen doses and five period of sample collection, on soil microbial biomass - nitrogen (SMB-N), total nitrogen (total N) and percentual ratio of the microbial biomass and total N (SMB-N/total N) in a Oxisol cultivated with barley (Hordeum vulgare L.). The experiment was installed in June, 2005, in an area located at Embrapa Cerrados, Federal District. The experimental design was a randomized block, with three replicates. The plots received doses of nitrogen: 20 - 40 - 80 kg ha-1 N and a control without it and the subplots were period of soil sample. Three applications of N were realized: 10 kg ha-1 on the 5th day (06/14) after sowing; the rest of N was parceled in two applications with fertigation, on tillage, on the 27th (07/08) DAP, e no 43rd (07/22) DAP. Soil samples layer (0 - 10 cm deep) were collected for (SMB-N) determination and total N in six periods: 02 days before of the first fertigation; 02 days after of the first fertigation; 04 days before of the last fertigation and 04 days after of the last fertigation; on flowering stage and after harvesting. There was effect of the doses of N and the period of soil collection on the SMB-N, total N and in the ratio SMB-N/total N. The average values of total N revealed steadier in short-term (cycle of the culture) and this was not a good parameter to evaluate the behavior and N transformations in the soil-plant system. Resumen: El objetivo de este estudio fue evaluar el efecto de diferentes dosis de nitrógeno y cinco período de muestreo en la biomasa microbiana del complejo suelo - nitrogeno (BMS-N), nitrógeno total (N total) y la relación porcentual de la biomasa microbiana y N total (BMS-N/N total) en un Oxisol cultivado con cebada (Hordeum vulgare L.). El estudio se inició en junio de 2005 en la estación experimental de la Empresa de Pesquisa Agropecuaria (Embrapa-Cerrados), Distrito Federal, Brazil. El experimento se dispuso en bloques al azar con tres repeticiones. Las parcelas recibieron dosis de nitrógeno: 20 - 40 - 80 kg/ha de N más un control sin N, y las subparcelas fueron el periodo de muestro. Las aplicaciones de N se realizaron de la forma siguiente: cinco días después de la siembra (dds) se aplicaron 10 kg/ha y el resto de la dosis se aplicó con fertirrigación en dos dosis 27 y 43 dds. Las muestras de suelo (0-10 cm de profundidad) para determinar BMS-N y N total fueron tomadas, 2 días antes e igual número de días después de la primera fertirrigación y 4 días antes y después de la última, en floración y después de la cosecha. No se encontró efecto de las dosis de N y el período de muestreo en el BMS-N, N total y en la relación BMS-N/N total. Los valores medios de N total fueron más estable en el corto plazo (ciclo de la cultivo) lo que indica que éste no es un buen parámetro para evaluar la dinámica del N y sus transformaciones en el sistema suelo-planta.

Changes In The Bacterial Community Of Soil From A Neutral Mine Drainage Channel.

Mine drainage is an important environmental disturbance that affects the chemical and biological components in natural resources. However, little is known about the effects of neutral mine drainage on the soil bacteria community. Here, a high-throughput 16S rDNA pyrosequencing approach was used to evaluate differences in composition, structure, and diversity of bacteria communities in samples from a neutral drainage channel, and soil next to the channel, at the Sossego copper mine in Brazil. Advanced statistical analyses were used to explore the relationships between the biological and chemical data. The results showed that the neutral mine drainage caused changes in the composition and structure of the microbial community, but not in its diversity. The Deinococcus/Thermus phylum, especially the Meiothermus genus, was in large part responsible for the differences between the communities, and was positively associated with the presence of copper and other heavy metals in the environmental samples. Other important parameters that influenced the bacterial diversity and composition were the elements potassium, sodium, nickel, and zinc, as well as pH. The findings contribute to the understanding of bacterial diversity in soils impacted by neutral mine drainage, and demonstrate that heavy metals play an important role in shaping the microbial population in mine environments.

Bacterial Diversity Assessment In Soil Of An Active Brazilian Copper Mine Using High-throughput Sequencing Of 16s Rdna Amplicons.

Mining activities pose severe environmental risks worldwide, generating extreme pH conditions and high concentrations of heavy metals, which can have major impacts on the survival of organisms. In this work, pyrosequencing of the V3 region of the 16S rDNA was used to analyze the bacterial communities in soil samples from a Brazilian copper mine. For the analysis, soil samples were collected from the slopes (geotechnical structures) and the surrounding drainage of the Sossego mine (comprising the Sossego and Sequeirinho deposits). The results revealed complex bacterial diversity, and there was no influence of deposit geographic location on the composition of the communities. However, the environment type played an important role in bacterial community divergence; the composition and frequency of OTUs in the slope samples were different from those of the surrounding drainage samples, and Acidobacteria, Chloroflexi, Firmicutes, and Gammaproteobacteria were responsible for the observed difference. Chemical analysis indicated that both types of sample presented a high metal content, while the amounts of organic matter and water were higher in the surrounding drainage samples. Non-metric multidimensional scaling (N-MDS) analysis identified organic matter and water as important distinguishing factors between the bacterial communities from the two types of mine environment. Although habitat-specific OTUs were found in both environments, they were more abundant in the surrounding drainage samples (around 50 %), and contributed to the higher bacterial diversity found in this habitat. The slope samples were dominated by a smaller number of phyla, especially Firmicutes. The bacterial communities from the slope and surrounding drainage samples were different in structure and composition, and the organic matter and water present in these environments contributed to the observed differences.