Impacto de la circulación ecuatorial en la zona mínima de oxógeno presente en el norte del Ecosistema de la Corriente del Humboldt

El Norte del Ecosistema de la Corriente Humboldt (NECH) constituye una de las mayores zonas de afloramiento, localizada en el borde oriental del Pacifico Sur, la cual presenta características particulares, entre las que destacan una alta producción primaria y la presencia de una de las zonas mínimas de oxigeno (ZMO) más intensas en el océano abierto. La ZMO presente en esta zona es producto de la alta demanda de oxígeno durante la remineralización de la materia orgánica, el largo tiempo de residencia de sus aguas y su poca ventilación. En el presente estudio nos enfocamos en estudiar la influencia de cambios en la ventilación en la ZMO del NECH, tomando como las principales fuentes de aporte de oxígeno en esta zona a la Corriente Sub-superficial Ecuatorial (EUC) y las Contracorrientes Sub-superficiales del Sur (SSCCs) o también conocidas como los Jets de Tsuchiya. Utilizamos el modelo acoplado físico-biogeoquímico ROMS-PISCES, para observar la sensibilidad de la ZMO a diferentes condiciones de la circulación ecuatorial provenientes de dos modelos oceánicos de circulación general (SODA y MERCATOR). Los resultados muestran que el flujo de oxígeno a los 88ºW disminuye latitudinalmente de la EUC a los SSCCs; además, se observa que la ZMO desaparece de los 4ºN - 4ºS en la simulación que presenta una circulación más intensa (RPSoda) por lo que se puede concluir que una intensificación de la circulación ecuatorial afectaría principalmente a la zona ecuatorial y no frente a Perú, debido a que una mayor ventilación sería compensada con un mayor consumo de oxigeno durante la remineralización, producto de una alta productividad generada por un mayor flujo de nutrientes.

Estudi per a la construcció de l'EDAR mitjançant filtres verds de la Nou de Gaià (Tarragonès)

El tractament de les aigües en nuclis menors de 2000 habitants es troba pendent de completar per part de l’Agència Catalana de l’Aigua més concretament al corresponent Pla de Sanejament d’Aigües Residuals Urbanes (PSARU). El nucli de La Nou de Gaià (al Tarragonès) es troba pendent de la construcció de la corresponent instal·lació de sanejament, projectada al 2007. Alternativament a les depuradores tradicionals basades en l’ús de formigó (o materials alternatius) i en la despesa elèctrica per assegurar una aeració i una evacuació dels fangs generats, existeixen tecnologies “toves”. Aquestes tecnologies, també conegudes com a “verdes”, es basen en imitar els sistemes naturals maximitzant el seu potencial d’autodepuració. A grans trets existeixen dos formes de depurar les aigües de forma ecològica”: llacunatges (existeix una capa d’aigua lliure) i filtres verds. El present estudi es basa en l’aplicació de filtres verds de morfologia vertical i flux subsuperficial, plantat amb canyes dels generes Scirpus o Phragmites. El resultat han estat 4 bases de 35*35 per a tractar un cabal de 150 m3/d i una població equivalent de 1272.

Numerical simulation of ambient seismic wavefield modification caused by pore-fluid effects in an oil reservoir

We have modeled numerically the seismic response of a poroelastic inclusion with properties applicable to an oil reservoir that interacts with an ambient wavefield. The model includes wave-induced fluid flow caused by pressure differences between mesoscopic-scale (i.e., in the order of centimeters to meters) heterogeneities. We used a viscoelastic approximation on the macroscopic scale to implement the attenuation and dispersion resulting from this mesoscopic-scale theory in numerical simulations of wave propagation on the kilometer scale. This upscaling method includes finite-element modeling of wave-induced fluid flow to determine effective seismic properties of the poroelastic media, such as attenuation of P- and S-waves. The fitted, equivalent, viscoelastic behavior is implemented in finite-difference wave propagation simulations. With this two-stage process, we model numerically the quasi-poroelastic wave-propagation on the kilometer scale and study the impact of fluid properties and fluid saturation on the modeled seismic amplitudes. In particular, we addressed the question of whether poroelastic effects within an oil reservoir may be a plausible explanation for low-frequency ambient wavefield modifications observed at oil fields in recent years. Our results indicate that ambient wavefield modification is expected to occur for oil reservoirs exhibiting high attenuation. Whether or not such modifications can be detected in surface recordings, however, will depend on acquisition design and noise mitigation processing as well as site-specific conditions, such as the geologic complexity of the subsurface, the nature of the ambient wavefield, and the amount of surface noise.

Interactions between magnesium, calcium, and aluminum on soybean root elongation

Alleviation of Al rhizotoxicity by Ca and Mg can differ among species and genotypes. Root elongation of soybean [Glycine max (L.) Merr.] line N93-S-179 and cvs. Young and Ransom exposed to varying concentrations of Al, Ca and Mg were compared in two experiments using a vertically split root system. Roots extending from a surface compartment with limed soil grew for 12 days into a subsurface compartment with nutrient solution treatments maintained at pH 4.6 with either 0 or 15 µmol L-1 Al. Calcium and Mg concentrations in treatments ranging from 0 to 20 mmol L-1. Although an adequate supply of Mg was provided in the surface soil compartment for soybean top growth, an inclusion of Mg was necessary in the subsurface solutions to promote root elongation in both the presence and absence of Al. In the absence of Al in the subsurface solution, tap root length increased by 74 % and lateral root length tripled when Mg in the solutions was increased from 0 to either 2 or 10 mmol L-1. In the presence of 15 µmol L-1 Al, additions of 2 or 10 mmol L-1 Mg increased tap root length fourfold and lateral root length by a factor of 65. This high efficacy of Mg may have masked differences in Al tolerance between genotypes N93 and Young. Magnesium was more effective than Ca in alleviating Al rhizotoxicity, and its ameliorative properties could not be accounted for by estimated electrostatic changes in root membrane potential and Al3+ activity at the root surface. The physiological mechanisms of Mg alleviation of Al injury in roots, however, are not known.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

Mass conservative three-dimensional water tracer distribution from Markov chain Monte Carlo inversion of time-lapse ground-penetrating radar data

Time-lapse geophysical measurements are widely used to monitor the movement of water and solutes through the subsurface. Yet commonly used deterministic least squares inversions typically suffer from relatively poor mass recovery, spread overestimation, and limited ability to appropriately estimate nonlinear model uncertainty. We describe herein a novel inversion methodology designed to reconstruct the three-dimensional distribution of a tracer anomaly from geophysical data and provide consistent uncertainty estimates using Markov chain Monte Carlo simulation. Posterior sampling is made tractable by using a lower-dimensional model space related both to the Legendre moments of the plume and to predefined morphological constraints. Benchmark results using cross-hole ground-penetrating radar travel times measurements during two synthetic water tracer application experiments involving increasingly complex plume geometries show that the proposed method not only conserves mass but also provides better estimates of plume morphology and posterior model uncertainty than deterministic inversion results.

Genesis of the sombric horizon in ultisols (red argisols) in southern Santa Catarina, Brazil

Dark subsurface horizons, with properties similar to the sombric horizon characterized by the USA Soil Taxonomy, are frequent in Southern Brazil. The genesis of this horizon is controversial and poorly understood. This study aimed to describe the occurrence of sombric-like horizons in Ultisols in the South of Santa Catarina State, at low altitudes, and suggest possible processes of humus transference, accumulation and persistence in these horizons. Physical, chemical and mineralogical properties of four Ultisols were evaluated; three were sampled in a toposequence, and another representative one in an isolated profile (RSP). The dark subsurface horizons coincide with the AB and BA transitional genetic horizons; they are acid, low in base saturation, and have a similar clay mineralogy in all horizons. Very high amounts of Fe and Al extracted by ammonium oxalate and sodium pyrophosphate solution as well as maximum Al extracted by CuCl2 solution were observed in these dark subsurface horizons, indicating a possible migration of these elements in the form of organometallic complexes. The contents of Al plus ½ Fe extracted from the RSP soil horizons with ammonium oxalate indicated spodic materials in the sombric-like horizon, although the soil morphology was not compatible with Spodosols. Maximum contents of fine clay were also found in the sombric-like horizon, suggesting Fe and Al migration as clay-humic substances. However, the hypothesis that sombric-like horizons in these soils are a relict feature of a grass paleovegetation, different from the current dense seasonal forest, should not be discarded but investigated in further studies.

Integration of ground-penetrating radar and microgravimetric methods to map shallow caves

Ground-penetrating radar (GPR) and microgravimetric surveys have been conducted in the southern Jura mountains of western Switzerland in order to map subsurface karstic features. The study site, La Grande Rolaz cave, is an extensive system in which many portions have been mapped. By using small station spacing and careful processing for the geophysical data, and by modeling these data with topographic information from within the cave, accurate interpretations have been achieved. The constraints on the interpreted geologic models are better when combining the geophysical methods than when using only one of the methods, despite the general limitations of two-dimensional (2D) profiling. For example, microgravimetry can complement GPR methods for accurately delineating a shallow cave section approximately 10 X 10 mt in size. Conversely, GPR methods can be complementary in determining cavity depths and in verifying the presence of off-line features and numerous areas of small cavities and fractures, which may be difficult to resolve in microgravimetric data.

Availability of soil water under tillage systems, mulch management and citrus rootstocks

The increased availability of soil water is important for the management of non-irrigated orange orchards. The objective of this study was to evaluate the availability of soil water in a Haplorthox (Rhodic Ferralsol) under different tillage systems used for orchard plantation, mulch management and rootstocks in a "Pêra" orange orchard in northwest Paraná, Brazil. An experiment in a split-split-plot design was established in 2002, in an area cultivated with Brachiaria brizantha grass in which three tillage systems (no tillage, conventional tillage and strip-tillage) were used for orchard plantation. This grass was mowed twice a year between the rows, representing two mulch managements in the split plots (no mulching and mulching in the plant rows). The split-split-plots were represented by two rootstocks ("Rangpur" lime and "Cleopatra" mandarin). The soil water content in the plant rows was evaluated in the 0-20 cm layer in 2007 and at 0-20 and 20-40 cm in 2008-2009. The effect of soil tillage systems prior to implantation of orange orchards on soil water availability was less pronounced than mulching and the rootstocks. The soil water availability was lower when "Pêra" orange trees were grafted on "Cleopatra" mandarin than on "Rangpur" lime rootstocks. Mulching had a positive influence on soil water availability in the sandy surface layer (0-20 cm) and sandy clay loam subsurface (20-40 cm) of the soil in the spring. The production of B. brizantha between the rows and residue disposal in the plant rows as mulch increased water availability to the "Pêra" orange trees.

Origin of a restraining bend in an evolving strike-slip system: The Cenozoic As Pontes basin (NW Spain)

The As Pontes basin (12 km2), NW Iberian Peninsula, is bounded by a double restraining bend of a dextral strike-slip fault, which is related to the western onshore end of the Pyrenean belt. Surface and subsurface data obtained from intensive coal exploration and mining in the basin since the 1960s together with additional structural and stratigraphic sequence analysis allowed us to determine the geometric relationships between tectonic structures and stratigraphic markers. The small size of the basin and the large amount of quality data make the As Pontes basin a unique natural laboratory for improving our understanding of the origin and evolution of restraining bends. The double restraining bend is the end stage of the structural evolution of a compressive underlapping stepover, where the basin was formed. During the first stage (stepover stage), which began ca. 30 Ma ago (latest Rupelian) and lasted 3.4 My, two small isolated basins bounded by thrusts and normal faults were formed. For 1.3 My, the strike-slip faults, which defined the stepover, grew towards each other until joining and forming the double restraining bend, which bounds one large As Pontes basin (transition stage). The history of the basin was controlled by the activity of the double restraining bend for a further 3.4 My (restraining bend stage) and ended in mid-Aquitanian times (ca. 22 Ma).

Forms of inorganic phosphorus in soil under different long term soil tillage systems and winter crops

The cultivation of crops with different capacity of P uptake and use under long-term soil tillage systems can affect the distribution of P cycling and inorganic forms in the soil, as a result of higher or lower use efficiency of P applied in fertilizers. The purpose of this study was to evaluate the effect of long-term cultivation of different winter species under tillage systems on the distribution of inorganic P forms in the soil. In 1986, the experiment was initiated with six winter crops (blue lupin, hairy vetch, oat, oilseed radish, wheat and fallow) on a Rhodic Hapludox in southwestern Paraná, under no-tillage (NT) and conventional tillage (CT). The application of phosphate fertilizer in NT rows increased inorganic P in the labile and moderately labile forms, and soil disturbance in CT redistributed the applied P in the deeper layers, increasing the moderately labile P concentration in the subsurface layers. Black oat and blue lupin were the most efficient P-recyclers and under NT, they increased the labile P content in the soil surface layers.

Characterization and classification of soils in the Taquari river basin - Pantanal region, state of Mato Grosso do Sul, Brazil

Among the soils in the Mato Grosso do Sul, stand out in the Pantanal biome, the Spodosols. Despite being recorded in considerable extensions, few studies aiming to characterize and classify these soils were performed. The purpose of this study was to characterize and classify soils in three areas of two physiographic types in the Taquari river basin: bay and flooded fields. Two trenches were opened in the bay area (P1 and P2) and two in the flooded field (P3 and P4). The third area (saline) with high sodium levels was sampled for further studies. In the soils in both areas the sand fraction was predominant and the texture from sand to sandy loam, with the main constituent quartz. In the bay area, the soil organic carbon in the surface layer (P1) was (OC) > 80 g kg-1, being diagnosed as Histic epipedon. In the other profiles the surface horizons had low OC levels which, associated with other properties, classified them as Ochric epipedons. In the soils of the bay area (P1 and P2), the pH ranged from 5.0 to 7.5, associated with dominance of Ca2+ and Mg2+, with base saturation above 50 % in some horizons. In the flooded fields (P3 and P4) the soil pH ranged from 4.9 to 5.9, H+ contents were high in the surface horizons (0.8-10.5 cmol c kg-1 ), Ca2+ and Mg² contents ranged from 0.4 to 0.8 cmol c kg-1 and base saturation was < 50 %. In the soils of the bay area (P1 and P2) iron was accumulated (extracted by dithionite - Fed) and OC in the spodic horizon; in the P3 and P4 soils only Fed was accumulated (in the subsurface layers). According to the criteria adopted by the Brazilian System of Soil Classification (SiBCS) at the subgroup level, the soils were classified as: P1: Organic Hydromorphic Ferrohumiluvic Spodosol. P2: Typical Orthic Ferrohumiluvic Spodosol. P3: Typical Hydromorphic Ferroluvic Spodosol. P4: Arenic Orthic Ferroluvic Spodosol.

Phosphorus accumulation and pollution potential in a hapludult fertilized with pig manure

Successive applications of pig litter to the soil surface can increase the phosphorus (P) content and alter its adsorption, promoting P transfer to surface or subsurface waters. The purpose of this study was to evaluate P accumulation and the pollution potential of a soil after application of pig litter. In March 2010, eight years after the installation of an experiment in Braço do Norte, Santa Catarina, SC, Brazil, on a Typic Hapludult, soil was sampled (layers 0-2.5, 2.5-5, 5-10, 10-15, 15-20 and 20-30 cm) after the following fertilization treatments: no pig litter fertilization, pig slurry application and pig manure application. In this period, 694 and 1,890 kg P2O5 ha-1 were applied in the treatments with pig slurry and pig manure, respectively. The P content was determined, based on Mehlich-1, anion exchange resin (AER), 0.01 mol L-1 CaCl2 and total P in the samples. The adsorption isotherm parameters were also determined by the Langmuir and Koski-Vähälä & Hartikainem models in the layers 0-2.5 and 20-30 cm. The application of 1,890 kg P2O5 ha-1 in the form of pig manure led to P accumulation, as evidenced by Mehlich-1, down to a depth of 15 cm, by AER and 0.01 mol L-1 CaCl2 down to 20 cm and by total P to 30 cm. After application of 1,890 kg P2O5 ha-1 in the form of pig manure, the values of maximum P adsorption capacity were lowest in the deepest layer (20-30 cm), indicating the occupation of part of the adsorption sites of the particles. The application of swine manure to the soil over eight years increased the P quantity in the soil solution of the surface layer, indicating environmental contamination risk for surface and subsurface waters.

Heavy metals in soils and plants in mango orchards in Petrolina, Pernambuco, Brazil

The monitoring of heavy metal concentrations in areas under intensive agriculture is essential for the agricultural sustainability and food safety. This paper evaluates the total contents of heavy metals in soils and mango trees in orchards of different ages (6, 7, 8, 9, 10, 11, 14, 16, 17, 19, and 26 years) in Petrolina, Pernambuco, Brazil. Soil samples were taken from the layers 0-20 cm and 20-40 cm, and mango leaves were collected in the growth stage. Areas of native vegetation (Caatinga) adjacent to the cultivated areas were used for comparison. The total concentrations of heavy metals (Cu, Cr, Fe, Zn, Mn, Ni, and Pb) were determined in soils and leaves. In general, mango cultivation led to Cu and Zn accumulation in the soil surface and to a reduction in the contents of Ni, Pb, Mn, and Fe in surface and subsurface. Since contamination by Cu, Zn, and Cr was detected, these areas must be monitored to prevent negative environmental impacts. For instance, the presence of Cr in mango tree leaves indicates the need to investigate the source of the element in these orchards. The management strategies of the different companies led to deficiency or excess of some metals in the evaluated areas. However, the Fe and Mn levels were adequate for the mineral nutrition of mango in all areas.

Estimating subsoil resistance to nitrate leaching from easily measurable pedological properties

Leaching of nitrate (NO3-) can increase the groundwater concentration of this anion and reduce the agronomical effectiveness of nitrogen fertilizers. The main soil property inversely related to NO3- leaching is the anion exchange capacity (AEC), whose determination is however too time-consuming for being carried out in soil testing laboratories. For this reason, this study evaluated if more easily measurable soil properties could be used to estimate the resistance of subsoils to NO3- leaching. Samples from the subsurface layer (20-40 cm) of 24 representative soils of São Paulo State were characterized for particle-size distribution and for chemical and electrochemical properties. The subsoil content of adsorbed NO3- was calculated from the difference between the NO3- contents extracted with 1 mol L-1 KCl and with water; furthermore, NO3- leaching was studied in miscible displacement experiments. The results of both adsorption and leaching experiments were consistent with the well-known role exerted by AEC on the nitrate behavior in weathered soils. Multiple regression analysis indicated that in subsoils with (i) low values of remaining phosphorus (Prem), (ii) low soil pH values measured in water (pH H2O), and (iii) high pH values measured in 1 moL L-1 KCl (pH KCl), the amounts of surface positive charges tend to be greater. For this reason, NO3- leaching tends to be slower in these subsoils, even under saturated flow condition.