Quality of surface water related to land use: a case study in a catchment with small farms and intensive vegetable crop production in southern Brazil

Water degradation is strongly related to agricultural activity. The aim of this study was to evaluate the influence of land use and some environmental components on surface water quality in the Campestre catchment, located in Colombo, state of Parana, Brazil. Physical and chemical attributes were analyzed (total nitrogen, ammonium, nitrate, total phosphorus, electrical conductivity, pH, temperature, turbidity, total solids, biological oxygen demand, chemical oxygen demand and dissolved oxygen). Monthly samples of the river water were taken over one year at eight monitoring sites, distributed over three sub-basins. Overall, water quality was worse in the sub-basin with a higher percentage of agriculture, and was also affected by a lower percentage of native forest and permanent preservation area, and a larger drainage area. Water quality was also negatively affected by the presence of agriculture in the riparian zone. In the summer season, probably due to higher rainfall and intensive soil use, a higher concentration of total nitrogen and particulate nitrogen was observed, as well as higher electrical conductivity, pH and turbidity. All attributes, except for total phosphorus, were in compliance with Brazilian Conama Resolution Nº 357/2005 for freshwater class 1. However, it should be noted that these results referred to the base flow and did not represent a discharge condition since most of the water samples were not collected at or near the rainfall event.

Soil fertility, nutrition and yield of maize and barley with gypsum application on soil surface in no-till

Annual crop yield and nutrition have shown differentiated responses to modifications in soil chemical properties brought about by gypsum application. The aim of this study was to evaluate the effect of gypsum application rates on the chemical properties of a Latossolo Bruno (Clayey Oxisol), as well as on the nutrition and yield of a maize-barley succession under no-till. The experiment was set up in November 2009 in Guarapuava, Parana, Brazil, applying gypsum rates of 0.0, 1.5, 3.0, 4.5, and 6.0 Mg ha-1 to the soil surface upon sowing maize, with crop succession of barley. Gypsum application decreased the levels of Al3+ and Mg2+ in the 0.0-0.1 m layer and increased soil pH in the layers from 0.2-0.6 m depth. Gypsum application has increased the levels of Ca2+ in all soil layers up to 0.6 m, and the levels of S-SO4(2-) up to 0.8 m. In both crops, the leaf concentrations of Ca and S were increased while Mg concentrations have decreased as a function of gypsum rates. There was also an effect of gypsum rates on grain yield, with a quadratic response of maize and a linear increase for barley. Yield increases were up to 11 and 12 % in relation to control for the maximum technical efficiency (MTE) rates of 3.8 and 6.0 Mg ha-1 of gypsum, respectively. Gypsum application improved soil fertility in the profile, especially in the subsurface, as well as plant nutrition, increasing the yields of maize and barley.

Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.

Effects of surface application of calcium-magnesium silicate and gypsum on soil fertility and sugarcane yield

Lime application recommendations for amendment of soil acidity in sugarcane were developed with a burnt cane harvesting system in mind. Sugarcane is now harvested in most areas without burning, and lime application for amendment of soil acidity in this system in which the sugarcane crop residue remains on the ground has been carried out without a scientific basis. The aim of this study was to evaluate the changes in soil acidity and stalk and sugar yield with different rates of surface application of calcium, magnesium silicate, and gypsum in ratoon cane. The experiment was performed after the 3rd harvest of the variety SP 81-3250 in a commercial green sugarcane plantation of the São Luiz Sugar Mill (47º 25' 33" W; 21º 59' 46" S), located in Pirassununga, São Paulo, in southeast Brazil. A factorial arrangement of four Ca-Mg silicate rates (0, 850, 1700, and 3400 kg ha-1) and two gypsum rates (0 and 1700 kg ha-1) was used in the experiment. After 12 months, the experiment was harvested and technological measurements of stalk and sugar yield were made. After harvest, soil samples were taken at the depths of 0.00-0.05, 0.05-0.10, 0.10-0.20, 0.20-0.40, and 0.40-0.60 m in all plots, and the following determinations were made: soil pH in CaCl2, organic matter, P, S, K, Ca, Mg, H+Al, Al, Si, and base saturation. The results show that the application of gypsum reduced the exchangeable Al3+ content and Al saturation below 0.05 m, and increased the Ca2+ concentration in the whole profile, the Mg2+ content below 0.10 m, K+ below 0.4 m, and base saturation below 0.20 m. This contributed to the effect of surface application of silicate on amendment of soil acidity reaching deeper layers. From the results of this study, it may be concluded that the silicate rate recommended may be too low, since the greater rates used in this experiment showed greater reduction in soil acidity, higher levels of nutrients at greater depths and an increase in stalk and sugar yield.

EFFECT OF SOIL TILLAGE AND PLANT RESIDUE ON SURFACE ROUGHNESS OF AN OXISOL UNDER SIMULATED RAIN

Surface roughness of the soil is formed by mechanical tillage and is also influenced by the kind and amount of plant residue, among other factors. Its persistence over time mainly depends on the fundamental characteristics of rain and soil type. However, few studies have been developed to evaluate these factors in Latossolos (Oxisols). In this study, we evaluated the effect of soil tillage and of amounts of plant residue on surface roughness of an Oxisol under simulated rain. Treatments consisted of the combination of the tillage systems of no-tillage (NT), conventional tillage (CT), and minimum tillage (MT) with rates of plant residue of 0, 1, and 2 Mg ha-1 of oats (Avena strigosa Schreb) and 0, 3, and 6 Mg ha-1 of maize (Zea mays L.). Seven simulated rains were applied on each experimental plot, with intensity of 60±2 mm h-1 and duration of 1 h at weekly intervals. The values of the random roughness index ranged from 2.94 to 17.71 mm in oats, and from 5.91 to 20.37 mm in maize, showing that CT and MT are effective in increasing soil surface roughness. It was seen that soil tillage operations carried out with the chisel plow and the leveling disk harrow are more effective in increasing soil roughness than those carried out with the heavy disk harrow and leveling disk harrow. The roughness index of the soil surface decreases exponentially with the increase in the rainfall volume applied under conditions of no tillage without soil cover, conventional tillage, and minimum tillage. The oat and maize crop residue present on the soil surface is effective in maintaining the roughness of the soil surface under no-tillage.

Organic Matter Fractions and Quality of the Surface Layer of a Constructed and Vegetated Soil After Coal Mining. I - Humic Substances and Chemical Characterization

After open coal mining, soils are “constructed”, which usually contain low levels and quality of organic matter (OM). Therefore, the use of plant species for revegetation and reclamation of degraded areas is essential. This study evaluated the distribution of carbon (C) in the chemical fractions as well as the chemical characteristics and humification degree of OM in a soil constructed after coal mining under cultivation of perennial grasses. The experiment was established in 2003 with the following treatments: Hemarthria altissima (T1), Paspalum notatum (T2), Cynodon dactilon (T3), Urochloa brizantha (T4), bare constructed soil (T5), and natural soil (T6). In 2009, soil samples were collected from the 0.00-0.03 m layer and the total organic carbon stock (TOC) and C stock in the chemical fractions: acid extract (CHCl), fulvic acid (CFA), humic acid (CHA), and humin (CHU) were determined. The humic acid (HA) fraction was characterized by infrared spectroscopy and the laser-induced fluorescence index (ILIF) of OM was also calculated. After six years, differences were only observed in the CHA stocks, which were highest in T1 (0.89 Mg ha-1) and T4 (1.06 Mg ha-1). The infrared spectra of HA in T1, T2 and T4 were similar to T6, with greater contribution of aliphatic organic compounds than in the other treatments. In this way, ILIF decreased in the sequence T5>T3>T4>T1>T2>T6, indicating higher OM humification in T3 and T5 and more labile OM in the other treatments. Consequently, the potential of OM quality recovery in the constructed soil was greatest in treatments T1 and T4.

Organic Matter Fractions and Quality of the Surface Layer of a Constructed and Vegetated Soil After Coal Mining. II - Physical Compartments and Carbon Management Index

Soils constructed after mining often have low carbon (C) stocks and low quality of organic matter (OM). Cover crops are decisive for the recovery process of these stocks, improving the quality of constructed soils. Therefore, the goal of this study was to evaluate the effect of cover crops on total organic C (TOC) stocks, C distribution in physical fractions of OM and the C management index (CMI) of a soil constructed after coal mining. The experiment was initiated in 2003 with six treatments: Hemarthria altissima (T1), Paspalum notatum (T2), Cynodon dactylon (T3), Urochloa brizantha (T4), bare constructed soil (T5), and natural soil (T6). Soil samples were collected in 2009 from the 0.00-0.03 m layer, and the TOC and C stocks in the physical particle size fractions (carbon in the coarse fraction - CCF, and mineral-associated carbon - MAC) and density fractions (free light fraction - FLF; occluded light fraction - OLF, and heavy fraction - HF) of OM were determined. The CMI components: carbon pool index (CPI), lability (L) and lability index (LI) were estimated by both fractionation methods. No differences were observed between TOC, CCF and MAC stocks. The lowest C stocks in FLF and OLF fractions were presented by T2, 0.86 and 0.61 Mg ha-1, respectively. The values of TOC stock, C stock in physical fractions and CMI were intermediate, greater than T5 and lower than T6 in all treatments, indicating the partial recovery of soil quality. As a result of the better adaptation of the species Hemarthria and Brizantha, resulting in greater accumulation of labile organic material, the CPI, L, LI and CMI values were higher in these treatments, suggesting a greater potential of these species for recovery of constructed soils.

SOIL PHOSPHORUS THRESHOLDS IN EVALUATING RISK OF ENVIRONMENTAL TRANSFER TO SURFACE WATERS IN SANTA CATARINA, BRAZIL

The State of Santa Catarina, Brazil, has agricultural and livestock activities, such as pig farming, that are responsible for adding large amounts of phosphorus (P) to soils. However, a method is required to evaluate the environmental risk of these high soil P levels. One possible method for evaluating the environmental risk of P fertilization, whether organic or mineral, is to establish threshold levels of soil available P, measured by Mehlich-1 extractions, below which there is not a high risk of P transfer from the soil to surface waters. However, the Mehlich-1 extractant is sensitive to soil clay content, and that factor should be considered when establishing such P-thresholds. The objective of this study was to determine P-thresholds using the Mehlich-1 extractant for soils with different clay contents in the State of Santa Catarina, Brazil. Soil from the B-horizon of an Oxisol with 800 g kg-1 clay was mixed with different amounts of sand to prepare artificial soils with 200, 400, 600, and 800 g kg-1 clay. The artificial soils were incubated for 30 days with moisture content at 80 % of field capacity to stabilize their physicochemical properties, followed by additional incubation for 30 days after liming to raise the pH(H2O) to 6.0. Soil P sorption curves were produced, and the maximum sorption (Pmax) was determined using the Langmuir model for each soil texture evaluated. Based on the Pmax values, seven rates of P were added to four replicates of each soil, and incubated for 20 days more. Following incubation, available P contents (P-Mehlich-1) and P dissolved in the soil solution (P-water) were determined. A change-point value (the P-Mehlich-1 value above which P-water starts increasing sharply) was calculated through the use of segmented equations. The maximum level of P that a soil might safely adsorb (P-threshold) was defined as 80 % of the change-point value to maintain a margin for environmental safety. The P-threshold value, in mg dm-3, was dependent on the soil clay content according to the model P-threshold = 40 + Clay, where the soil clay content is expressed as a percentage. The model was tested in 82 diverse soil samples from the State of Santa Catarina and was able to distinguish samples with high and low environmental risk.

Surface Complexation Modeling in Variable Charge Soils: Charge Characterization by Potentiometric Titration

ABSTRACT Intrinsic equilibrium constants of 17 representative Brazilian Oxisols were estimated from potentiometric titration measuring the adsorption of H+ and OH− on amphoteric surfaces in suspensions of varying ionic strength. Equilibrium constants were fitted to two surface complexation models: diffuse layer and constant capacitance. The former was fitted by calculating total site concentration from curve fitting estimates and pH-extrapolation of the intrinsic equilibrium constants to the PZNPC (hand calculation), considering one and two reactive sites, and by the FITEQL software. The latter was fitted only by FITEQL, with one reactive site. Soil chemical and physical properties were correlated to the intrinsic equilibrium constants. Both surface complexation models satisfactorily fit our experimental data, but for results at low ionic strength, optimization did not converge in FITEQL. Data were incorporated in Visual MINTEQ and they provide a modeling system that can predict protonation-dissociation reactions in the soil surface under changing environmental conditions.

Surface Complexation Modeling in Variable Charge Soils: Prediction of Cadmium Adsorption

ABSTRACT Intrinsic equilibrium constants for 22 representative Brazilian Oxisols were estimated from a cadmium adsorption experiment. Equilibrium constants were fitted to two surface complexation models: diffuse layer and constant capacitance. Intrinsic equilibrium constants were optimized by FITEQL and by hand calculation using Visual MINTEQ in sweep mode, and Excel spreadsheets. Data from both models were incorporated into Visual MINTEQ. Constants estimated by FITEQL and incorporated in Visual MINTEQ software failed to predict observed data accurately. However, FITEQL raw output data rendered good results when predicted values were directly compared with observed values, instead of incorporating the estimated constants into Visual MINTEQ. Intrinsic equilibrium constants optimized by hand calculation and incorporated in Visual MINTEQ reliably predicted Cd adsorption reactions on soil surfaces under changing environmental conditions.

Understanding and modulating the competitive surface-adsorption of proteins through coarse-grained molecular dynamics simulations

It is now well accepted that cellular responses to materials in a biological medium reflect greatly the adsorbed biomolecular layer, rather than the material itself. Here, we study by molecular dynamics simulations the competitive protein adsorption on a surface (Vroman effect), i.e. the non-monotonic behavior of the amount of protein adsorbed on a surface in contact with plasma as functions of contact time and plasma concentration. We find a complex behavior, with regimes during which small and large proteins are not necessarily competing between them, but are both competing with others in solution ("cooperative" adsorption). We show how the Vroman effect can be understood, controlled and inverted.