Spatial and temporal variations of trace element distribution in soils and street dust of an industrial town in NW Spain: 15 years of study

Extensive spatial and temporal surveys, over 15 years, have been conducted in soil in urban parks and street dusts in one of the most polluted cities in western Europe, Avilés (NW Spain). The first survey was carried out in 1996, and since then monitoring has been undertaken every five years. Whilst the sampling site is a relatively small town, industrial activities (mainly the steel industry and Zn and Al metallurgy) and other less significant urban sources, such as traffic, strongly affect the load of heavy metals in the urban aerosol. Elemental tracers have been used to characterise the influence of these sources on the composition of soil and dust. Although PM10 has decreased over these years as a result of environmental measures undertaken in the city, some of the “industrial” elements still remain in concentrations of concern for example, up to 4.6% and 0.5% of Zn in dust and soil, respectively. Spatial trends in metals such as Zn and Cd clearly reflect sources from the processing industries. The concentrations of these elements across Europe have reduced over time, however the most recent results from Avilés revealed an upward trend in concentration for Zn, Cd, Hg and As. A risk assessment of the soil highlighted As as an element of concern since its cancer risk in adults was more than double the value above which regulatory agencies deem it to be unacceptable. If children were considered to be the receptors, then the risk nearly doubles from this element.

Enhanced radon emission in natural CO2 flows in Campo de Calatrava region (central Spain)

Activity of radon gas in natural soils is commonly low (in the order of few thousands of Bq·m-3) due to the fast decay (half-life= 3.8 days in the case of 222Rn) that prevents accumulation in soil pores. Exceptionally, high Rn soil activity (up to 430 KBq·m-3) is found around point sources of deep CO2 fluxes. These fluxes allow the transport of trace gases (including Rn) to long distances in the geosphere leading to a potential hazard as Rn accumulation in buildings. CO2 degassing is common in active or ancient volcanic fields and occurs as free gas fluxes or dissolved in groundwater. In this work, the occurrence of Rnbearing, CO2 fluxes from the Campo de Calatrava region in Central Spain has been studied in order to determine their (1) magnitude, (2) migration paths and (3) potential impact on the environment, and (4) methodologies to best detection and measurement.

Morphological Functions with Parallel Sets for the Pore Space of X-ray CT Images of Soil Columns

During the last few decades, new imaging techniques like X-ray computed tomography have made available rich and detailed information of the spatial arrangement of soil constituents, usually referred to as soil structure. Mathematical morphology provides a plethora of mathematical techniques to analyze and parameterize the geometry of soil structure. They provide a guide to design the process from image analysis to the generation of synthetic models of soil structure in order to investigate key features of flow and transport phenomena in soil. In this work, we explore the ability of morphological functions built over Minkowski functionals with parallel sets of the pore space to characterize and quantify pore space geometry of columns of intact soil. These morphological functions seem to discriminate the effects on soil pore space geometry of contrasting management practices in a Mediterranean vineyard, and they provide the first step toward identifying the statistical significance of the observed differences.

Availability and uptake of trace elements in a forage rotation under conservation and plough tillage

After 14 years under conventional plough tillage (CT) or conservation minimum tillage (MT), the soil available Al, Fe, Mn, Cu and Zn (0-5, 5-15 and 15-30 cm layers) and their plant uptake were evaluated during two years in a ryegrass-maize forage rotation in NW Spain (t emperate-humid region). The three-way ANOVA showed that trace element concentrations in soil were mainly influenced by sampling date, followed by soil depth and tillage system (35-73 %, 7-58 % and 3- 11 % of variance explained, respectively). Excepting for Fe (CT) and Al (CT and MT), the elemental concentrations decreased with depth, the stratification being stronger under MT. For soil available Al, Fe, Mn and Cu, the concentrations were higher in CT than in MT (5-15 and 15-30 cm layers) or were not affected by tillage system (0-5 cm). In contrast, the available Zn contents were higher in MT than CT at the soil surface and did not differ in deeper layers. The concentration of Al, Fe and Cu in crops were not influenced by tillage system, which explain 22 % of Mn variance in maize (CT > MT in the more humid year) and 18 % of Zn variance in ryegrass (MT > CT in both years). However, in the summer crop (maize) the concentrations of Fe, Mn and Zn tended to be higher in MT than in CT under drought conditions, while the opposite was true in the year without water limitation. Therefore, under the studied conditions of climate, soil, tillage and crop rotation, little influence of tillage system on crop nutritive value would be expected. To minimize the potential deficiency of Zn (maize) and Cu (maize and ryegrass) on crop yields the inclusion of these micro-nutrients in fertilization schedule is reco mmended, as well as liming to alleviate Al toxicity on maize crops.

Risk of Leaching in Soils Amended by Compost and digestate from Municipal Solid Waste

New European directives have proposed the direct application of compost and digestate produced from municipal solid wastes as organic matter sources in agricultural soils. Therefore information about phosphorus leaching from these residues when they are applied to the soil is increasingly mportant. Leaching experiments were conducted to determine the P mobility in compost and digestate mixtures, supplying equivalent amounts to 100 kg P ha?1 to three different types of soils. The tests were performed in accordance with CEN/TS 14405:2004 analyzing the maximum dissolved reactive P and the kinetic rate in the leachate. P biowaste fractionation indicated that digestate has a higher level of available P than compost has. In contrast, P losses in leaching experiments with soil-compost mixtureswere higher than in soil-digestate mixtures. For bothwastes, therewas no correlation between disolved reactive P lost and the water soluble P.The interaction between soil and waste, the long experimentation time, and the volume of leachate obtained caused the waste?s wettability to become an influential parameter in P leaching behavior. The overall conclusion is that kinetic data analysis provides valuable information concerning the sorption mechanism that can be used for predicting the large-scale behavior of soil systems.

Geometrical characterization of undisturbed soil samples using X-ray computed tomography image analysis. Effect of soil management on soil structure

El estudio de la estructura del suelo es de vital importancia en diferentes campos de la ciencia y la tecnología. La estructura del suelo controla procesos físicos y biológicos importantes en los sistemas suelo-planta-microorganismos. Estos procesos están dominados por la geometría de la estructura del suelo, y una caracterización cuantitativa de la heterogeneidad de la geometría del espacio poroso es beneficiosa para la predicción de propiedades físicas del suelo. La tecnología de la tomografía computerizada de rayos-X (CT) nos permite obtener imágenes digitales tridimensionales del interior de una muestra de suelo, proporcionando información de la geometría de los poros del suelo y permitiendo el estudio de los poros sin destruir las muestras. Las técnicas de la geometría fractal y de la morfología matemática se han propuesto como una poderosa herramienta para analizar y cuantificar características geométricas. Las dimensiones fractales del espacio poroso, de la interfaz poro-sólido y de la distribución de tamaños de poros son indicadores de la complejidad de la estructura del suelo. Los funcionales de Minkowski y las funciones morfológicas proporcionan medios para medir características geométricas fundamentales de los objetos geométricos tridimensionales. Esto es, volumen, superficie, curvatura media de la superficie y conectividad. Las características del suelo como la distribución de tamaños de poros, el volumen del espacio poroso o la superficie poro-solido pueden ser alteradas por diferentes practicas de manejo de suelo. En este trabajo analizamos imágenes tomográficas de muestras de suelo de dos zonas cercanas con practicas de manejo diferentes. Obtenemos un conjunto de medidas geométricas, para evaluar y cuantificar posibles diferencias que el laboreo pueda haber causado en el suelo. ABSTRACT The study of soil structure is of vital importance in different fields of science and technology. Soil structure controls important physical and biological processes in soil-plant-microbial systems. Those processes are dominated by the geometry of soil pore structure, and a quantitative characterization of the spatial heterogeneity of the pore space geometry is beneficial for prediction of soil physical properties. The technology of X-ray computed tomography (CT) allows us to obtain three-dimensional digital images of the inside of a soil sample providing information on soil pore geometry and enabling the study of the pores without disturbing the samples. Fractal geometry and mathematical morphological techniques have been proposed as powerful tools to analyze and quantify geometrical features. Fractal dimensions of pore space, pore-solid interface and pore size distribution are indicators of soil structure complexity. Minkowski functionals and morphological functions provide means to measure fundamental geometrical features of three-dimensional geometrical objects, that is, volume, boundary surface, mean boundary surface curvature, and connectivity. Soil features such as pore-size distribution, pore space volume or pore-solid surface can be altered by different soil management practices. In this work we analyze CT images of soil samples from two nearby areas with contrasting management practices. We performed a set of geometrical measures, some of them from mathematical morphology, to assess and quantify any possible difference that tillage may have caused on the soil.

Effects of biochar prepared from organic waste on soil properties

Biochar is a carbon-rich solid obtained by the thermal decomposition of organic matter under a limited supply of oxygen and at relatively low temperatures. Biochar can be prepared from the pyrolysis of different organic feed- stocks, such as wood and biomass crops, agricultural by-products, different types of waste or paper industry waste materials . The pyrolysis procedure of waste, i.e. sewage sludge, has mainly two advantages, firstly, it removes pathogens from waste and, secondly, biochar can reduce the leaching of heavy metals present in raw sewage sludge. This trend of the use of waste material as feedstocks to the preparation of biochar is increasing in the last years due to industrial development and economic growth imply an increase in waste generation. The application of biochar may have positive effects on soil physical properties as water holding capacity and structure or on soil biological activity and soil quality. Also, biochar can be used to remove water pollutants and can be used in multiple ways in soil remediation due to its adsorption of pesticides or metals. Also, biochar contribute to carbon sequestration due to carbon stability of biochar materials. The objective of this presentation is to review the positive effects of the biochar prepared from organic waste on soil properties.

Short-term response of the Ca cycle of a montane forest in Ecuador to low experimental CaCl2 additions

The tropical montane forests of the E Andean cordillera in Ecuador receive episodic Sahara- dust inputs particularly increasing Ca deposition. We added CaCl2 to isolate the effect of Ca deposition by Sahara dust to tropical montane forest from the simultaneously occurring pH effect. We examined components of the Ca cycle at four control plots and four plots with added Ca (2 × 5 kg ha?1 Ca annually as CaCl2) in a random arrangement. Between August 2007 and December 2009 (four applications of Ca), we determined Ca concentrations and fluxes in litter leachate, mineral soil solution (0.15 and 0.30 m depths), throughfall, and fine litterfall and Al con- centrations and speciation in soil solutions. After 1 y of Ca addition, we assessed fine-root bio- mass, leaf area, and tree growth. Only < 3% of the applied Ca leached below the acid organic layer (pH 3.5?4.8). The added CaCl2 did not change electrical conductivity in the root zone after 2 y. In the second year of fertilization, Ca retention in the canopy of the Ca treatment tended to decrease relative to the control. After 2 y, 21% of the applied Ca was recycled to soil with throughfall and litterfall. One year after the first Ca addition, fine-root biomass had decreased significantly. Decreasing fine-root biomass might be attributed to a direct or an indirect beneficial effect of Ca on the soil decomposer community. Because of almost complete association of Al with dissolved organic matter and high free Ca2+ : Al3+ activity ratios in solution of all plots, Al toxicity was unlikely. We conclude that the added Ca was retained in the system and had benefi- cial effects on some plants.

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.