Removal of paracetamol on biomass-derived activated carbon: Modeling the fixed bed breakthrough curves using batch adsorption experiments

The remediation of paracetamol (PA), an emerging contaminant frequently found in wastewater treatment plants, has been studied in the low concentration range (0.3–10 mg L−1) using as adsorbent a biomass-derived activated carbon. PA uptake of up to 100 mg g−1 over the activated carbon has been obtained, with the adsorption isotherms being fairly explained by the Langmuir model. The application of Reichemberg and the Vermeulen equations to the batch kinetics experiments allowed estimating homogeneous and heterogeneous diffusion coefficients, reflecting the dependence of diffusion with the surface coverage of PA. A series of rapid small-scale column tests were carried out to determine the breakthrough curves under different operational conditions (temperature, PA concentration, flow rate, bed length). The suitability of the proposed adsorbent for the remediation of PA in fixed-bed adsorption was proven by the high PA adsorption capacity along with the fast adsorption and the reduced height of the mass transfer zone of the columns. We have demonstrated that, thanks to the use of the heterogeneous diffusion coefficient, the proposed mathematical approach for the numerical solution to the mass balance of the column provides a reliable description of the breakthrough profiles and the design parameters, being much more accurate than models based in the classical linear driving force.

Imaging and modeling of pore scale processes in porous media using X-ray computed tomography and lattice Boltzmann solver

In der Erdöl– und Gasindustrie sind bildgebende Verfahren und Simulationen auf der Porenskala im Begriff Routineanwendungen zu werden. Ihr weiteres Potential lässt sich im Umweltbereich anwenden, wie z.B. für den Transport und Verbleib von Schadstoffen im Untergrund, die Speicherung von Kohlendioxid und dem natürlichen Abbau von Schadstoffen in Böden. Mit der Röntgen-Computertomografie (XCT) steht ein zerstörungsfreies 3D bildgebendes Verfahren zur Verfügung, das auch häufig für die Untersuchung der internen Struktur geologischer Proben herangezogen wird. Das erste Ziel dieser Dissertation war die Implementierung einer Bildverarbeitungstechnik, die die Strahlenaufhärtung der Röntgen-Computertomografie beseitigt und den Segmentierungsprozess dessen Daten vereinfacht. Das zweite Ziel dieser Arbeit untersuchte die kombinierten Effekte von Porenraumcharakteristika, Porentortuosität, sowie die Strömungssimulation und Transportmodellierung in Porenräumen mit der Gitter-Boltzmann-Methode. In einer zylindrischen geologischen Probe war die Position jeder Phase auf Grundlage der Beobachtung durch das Vorhandensein der Strahlenaufhärtung in den rekonstruierten Bildern, das eine radiale Funktion vom Probenrand zum Zentrum darstellt, extrahierbar und die unterschiedlichen Phasen ließen sich automatisch segmentieren. Weiterhin wurden Strahlungsaufhärtungeffekte von beliebig geformten Objekten durch einen Oberflächenanpassungsalgorithmus korrigiert. Die Methode der „least square support vector machine” (LSSVM) ist durch einen modularen Aufbau charakterisiert und ist sehr gut für die Erkennung und Klassifizierung von Mustern geeignet. Aus diesem Grund wurde die Methode der LSSVM als pixelbasierte Klassifikationsmethode implementiert. Dieser Algorithmus ist in der Lage komplexe geologische Proben korrekt zu klassifizieren, benötigt für den Fall aber längere Rechenzeiten, so dass mehrdimensionale Trainingsdatensätze verwendet werden müssen. Die Dynamik von den unmischbaren Phasen Luft und Wasser wird durch eine Kombination von Porenmorphologie und Gitter Boltzmann Methode für Drainage und Imbibition Prozessen in 3D Datensätzen von Böden, die durch synchrotron-basierte XCT gewonnen wurden, untersucht. Obwohl die Porenmorphologie eine einfache Methode ist Kugeln in den verfügbaren Porenraum einzupassen, kann sie dennoch die komplexe kapillare Hysterese als eine Funktion der Wassersättigung erklären. Eine Hysterese ist für den Kapillardruck und die hydraulische Leitfähigkeit beobachtet worden, welche durch die hauptsächlich verbundenen Porennetzwerke und der verfügbaren Porenraumgrößenverteilung verursacht sind. Die hydraulische Konduktivität ist eine Funktion des Wassersättigungslevels und wird mit einer makroskopischen Berechnung empirischer Modelle verglichen. Die Daten stimmen vor allem für hohe Wassersättigungen gut überein. Um die Gegenwart von Krankheitserregern im Grundwasser und Abwässern vorhersagen zu können, wurde in einem Bodenaggregat der Einfluss von Korngröße, Porengeometrie und Fluidflussgeschwindigkeit z.B. mit dem Mikroorganismus Escherichia coli studiert. Die asymmetrischen und langschweifigen Durchbruchskurven, besonders bei höheren Wassersättigungen, wurden durch dispersiven Transport aufgrund des verbundenen Porennetzwerks und durch die Heterogenität des Strömungsfeldes verursacht. Es wurde beobachtet, dass die biokolloidale Verweilzeit eine Funktion des Druckgradienten als auch der Kolloidgröße ist. Unsere Modellierungsergebnisse stimmen sehr gut mit den bereits veröffentlichten Daten überein.

Lattice Boltzmann modeling of fluid flow and solute transport in karst aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Lattice Boltzmann Modeling of Fluid Flow and Solute Transport in Karst Aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Inferring transport characteristics in a fractured rock aquifer by combining single-hole ground-penetrating radar reflection monitoring and tracer test data

Investigations of solute transport in fractured rock aquifers often rely on tracer test data acquired at a limited number of observation points. Such data do not, by themselves, allow detailed assessments of the spreading of the injected tracer plume. To better understand the transport behavior in a granitic aquifer, we combine tracer test data with single-hole ground-penetrating radar (GPR) reflection monitoring data. Five successful tracer tests were performed under various experimental conditions between two boreholes 6 m apart. For each experiment, saline tracer was injected into a previously identified packed-off transmissive fracture while repeatedly acquiring single-hole GPR reflection profiles together with electrical conductivity logs in the pumping borehole. By analyzing depth-migrated GPR difference images together with tracer breakthrough curves and associated simplified flow and transport modeling, we estimate (1) the number, the connectivity, and the geometry of fractures that contribute to tracer transport, (2) the velocity and the mass of tracer that was carried along each flow path, and (3) the effective transport parameters of the identified flow paths. We find a qualitative agreement when comparing the time evolution of GPR reflectivity strengths at strategic locations in the formation with those arising from simulated transport. The discrepancies are on the same order as those between observed and simulated breakthrough curves at the outflow locations. The rather subtle and repeatable GPR signals provide useful and complementary information to tracer test data acquired at the outflow locations and may help us to characterize transport phenomena in fractured rock aquifers.

Temperature effect in potassium and nitrate ions in soil transport

When doing researches on solute dynamics in porous medium, the knowledge of medium characteristics and percolating liquids, as well as of external factors is very important. An important external factor is temperature and, in this sense, our purpose was determining potassium and nitrate transport parameters for different values of temperature, in miscible displacement experiments. Evaluated parameters were retardation factor (R), diffusion/dispersion coefficient (D) and dispersivity, at ambient temperature (25 up to 28 ºC), 40 ºC and 50 ºC. Salts used were potassium nitrate and potassium chlorate, prepared in a solution made up of 5 ppm nitrate and 2.000 ppm potassium, with Red-Yellow Latosol porous medium. Temperature exhibited a positive influence upon porous medium solution and upon dispersion coefficient.

HYBRID ARTIFICIAL NEURAL NETWORK APPLIEDTO MODELING SCFE OF BASIL AND ROSEMARY OILS

This work presents the results of a Hybrid Neural Network (HNN) technique as applied to modeling SCFE curves obtained from two Brazilian vegetable matrices. A series Hybrid Neural Network was employed to estimate the parameters of the phenomenological model. A small set of SCFE data of each vegetable was used to generate an extended data set, sufficient to train the network. Afterwards, other sets of experimental data, not used in the network training, were used to validate the present approach. The series HNN correlates well the experimental data and it is shown that the predictions accomplished with this technique may be promising for SCFE purposes.

Avaliação da capacidade de adsorção do CO2 em zeólita 12X com gases sintéticos e originados da pirólise do lodo de esgoto

This work depicts a study of the adsorption of carbon dioxide on zeolite 13X. The activities were divided into four stages: study batch adsorption capacity of the adsorbent with synthetic CO2 (4%), fixed bed dynamic evaluation with the commercial mixture of gases (4% CO2, 1.11% CO, 1 2% H2, 0.233% CH4, 0.1% C3, 0.0233% C4 argon as inert closing balance), fixed bed dynamic modeling and evaluation of the breakthrough curve of CO2 originated from the pyrolysis of sewage sludge. The sewage sludge and the adsorbent were characterized by analysis TG / DTA, SEM, XRF and BET. Adsorption studies were carried out under the following operating conditions: temperature 40 °C (for the pyrolysis of the sludge T = 600 °C), pressures of 0.55 to 5.05 bar (batch process), flow rate of the gaseous mixture between 50 - 72 ml/min and the adsorbent masses of 10, 15 and 20 g (fixed bed process). The time for the adsorption batch was 7 h and on the fixed bed was around 180 min. The results of this study showed that in batch adsorption process step with zeolite 13X is efficient and the mass of adsorbed CO2 increases with the increases pressure, decreases with temperature increases and rises due the increase of activation temperature adsorbent. In the batch process were evaluated the breakthrough curves, which were compared with adsorption isotherms represented by the models of Langmuir, Freündlich and Toth. All models well adjusted to the experimental points, but the Langmuir model was chosen in view of its use in the dynamic model does not have implications for adsorption (indeterminacy and larger number of parameters such as occurred with others) in solving the equation. In the fixed bed dynamic study with the synthetic gas mixture, 20 g of mass adsorbent showed the maximum adsorption percentage 46.7% at 40 °C temperature and 50 mL/min of flow rate. The model was satisfactorily fitted to the three breakthrough curves and the parameters were: axial dispersion coefficient (0.0165 dm2/min), effective diffusivity inside the particle (dm2/min 0.0884) and external transfer coefficient mass (0.45 dm/min). The breakthrough curve for CO2 in the process of pyrolysis of the sludge showed a fast saturation with traces of aerosols presents in the gas phase into the fixed bed under the reaction process

Recuperação e purificação de quitosanases usando adsorção em leito expandido com streamline DEAE com modelagem e simulação usando redes neurais

Expanded Bed Adsorption (EBA) is an integrative process that combines concepts of chromatography and fluidization of solids. The many parameters involved and their synergistic effects complicate the optimization of the process. Fortunately, some mathematical tools have been developed in order to guide the investigation of the EBA system. In this work the application of experimental design, phenomenological modeling and artificial neural networks (ANN) in understanding chitosanases adsorption on ion exchange resin Streamline® DEAE have been investigated. The strain Paenibacillus ehimensis NRRL B-23118 was used for chitosanase production. EBA experiments were carried out using a column of 2.6 cm inner diameter with 30.0 cm in height that was coupled to a peristaltic pump. At the bottom of the column there was a distributor of glass beads having a height of 3.0 cm. Assays for residence time distribution (RTD) revelead a high degree of mixing, however, the Richardson-Zaki coefficients showed that the column was on the threshold of stability. Isotherm models fitted the adsorption equilibrium data in the presence of lyotropic salts. The results of experiment design indicated that the ionic strength and superficial velocity are important to the recovery and purity of chitosanases. The molecular mass of the two chitosanases were approximately 23 kDa and 52 kDa as estimated by SDS-PAGE. The phenomenological modeling was aimed to describe the operations in batch and column chromatography. The simulations were performed in Microsoft Visual Studio. The kinetic rate constant model set to kinetic curves efficiently under conditions of initial enzyme activity 0.232, 0.142 e 0.079 UA/mL. The simulated breakthrough curves showed some differences with experimental data, especially regarding the slope. Sensitivity tests of the model on the surface velocity, axial dispersion and initial concentration showed agreement with the literature. The neural network was constructed in MATLAB and Neural Network Toolbox. The cross-validation was used to improve the ability of generalization. The parameters of ANN were improved to obtain the settings 6-6 (enzyme activity) and 9-6 (total protein), as well as tansig transfer function and Levenberg-Marquardt training algorithm. The neural Carlos Eduardo de Araújo Padilha dezembro/2013 9 networks simulations, including all the steps of cycle, showed good agreement with experimental data, with a correlation coefficient of approximately 0.974. The effects of input variables on profiles of the stages of loading, washing and elution were consistent with the literature

Migração de solutos em basalto fraturado: quantificação experimental em laboratório e validação matemática

A avaliação do risco a contaminação e a escolha de técnicas de remediação de poluentes em aquíferos fraturados depende da quantificação dos fenômenos envolvidos no transporte de solutos. A geometria da fratura, usualmente caracterizada pela abertura, é o principal parâmetro que indiretamente controla o transporte nos aquíferos fraturados. A simplificação mais comum desse problema é assumir que as fraturas são um par de placas planas e paralelas, isto é, com uma abertura constante. No entanto, por causa do limitado número de trabalhos experimentais, não está esclarecida a adequabilidade do uso de uma abertura constante para simular o transporte conservativo em fraturas do Aquífero Serra Geral (ASG), Brasil. O objetivo deste trabalho é avaliar a influência da abertura de uma fratura natural do Aquífero Serra Geral sob o transporte conservativo de solutos. Uma amostra natural de basalto fraturado foi usada em um experimento hidráulico e de transporte de um traçador conservativo (escala de laboratório). O campo de abertura foi medido usando a técnica avançada, de alta resolução e tridimensional, chamada microtomografia computadorizada de raios-X. A concentração de traçador medida foi utilizada para validar uma solução analítica unidimensional da Equação de Advecção-dispersão (ADE). O desemprenho do ajuste da ADE às curvas de passagem experimentais foi avaliado para quatro diferentes tipos de aberturas constantes. Os resultados mostraram que o escoamento de água e o transporte de contaminantes pode ocorrer através de fraturas micrométricas, ocasionando, eventualmente, a contaminação do ASG. A abertura de balanço de massa é a única que pode ser chamada propriamente de \"abertura equivalente\". O uso de aberturas constantes na ADE não permitiu representar completamente o formato das curvas de passagem porque o campo de velocidade não é uniforme e intrinsicamente bidimensional. Portanto, na simulação do transporte deve-se incorporar a heterogeneidade da abertura da fratura.

Purificação, modelagem e simulação da adsorção de bioprodutos por troca iônica em coluna de leito expandido

A produção de proteínas através de microrganismos tornou-se uma técnica muito importante na obtenção de compostos de interesse da indústria farmacêutica e alimentícia. Extratos brutos nos quais as proteínas são obtidas são geralmente complexos, contendo sólidos e células em suspensão. Usualmente, para uso industrial destes compostos, é necessário obtê-los puros, para garantir a sua atuação sem interferência. Um método que vem recebendo destaque especialmente nos últimos 10 anos é o uso da cromatografia de troca iônica em leito expandido, que combina em uma única etapa os passos de clarificação, concentração e purificação da molécula alvo, reduzindo assim o tempo de operação e também os custos com equipamentos para realização de cada etapa em separado. Combinado a este fato, a última década também é marcada por trabalhos que tratam da modelagem matemática do processo de adsorção de proteínas em resinas. Está técnica, além de fornecer informações importantes sobre o processo de adsorção, também é de grande valia na otimização da etapa de adsorção, uma vez que permite que simulações sejam feitas, sem a necessidade de gasto de tempo e material com experimentos em bancada, especialmente se é desejado uma ampliação de escala. Dessa forma, o objetivo desta tese foi realizar a modelagem e simulação do processo de adsorção de bioprodutos em um caldo bruto na presença de células, usando inulinase e C-ficocianina como objeto de estudo e purificar C-ficocianina utilizando resina de troca iônica em leito expandido. A presente tese foi então dividida em quatro artigos. O primeiro artigo teve como objeto de estudo a enzima inulinase, e a otimização da etapa de adsorção desta enzima em resina de troca iônica Streamline SP, em leito expandido, foi feita através da modelagem matemática e simulação das curvas de ruptura em três diferentes graus de expansão (GE). As máximas eficiências foram observadas quando utilizadas maiores concentrações de inulinase (120 a 170 U/mL), e altura de leito entre 20 e 30 cm. O grau de expansão de 3,0 vezes foi considerado o melhor, uma vez que a produtividade foi consideravelmente superior. O segundo artigo apresenta o estudo das condições de adsorção de C-ficocianina em resina de troca iônica, onde foi verificado o efeito do pH e temperatura na adsorção e após construída a isoterma de adsorção. A isoterma de adsorção da C-ficocianina em resina Streamline Q XL feita em pH 7,5 e a 25°C (ambiente), apresentou um bom ajuste ao modelo de Langmuir (R=0,98) e os valores qm (capacidade máxima de adsorção) e Kd (constante de equilíbrio) estimados pela equação linearizada da isoterma, foram de 26,7 mg/mL e 0,067mg/mL. O terceiro artigo aborda a modelagem do processo de adsorção de extrato não clarificado de C-ficocianina em resina de troca iônica Streamline Q XL em coluna de leito expandido. Três curvas de ruptura foram feitas em diferentes graus de expansão (2,0, 2,5 e 3,0). A condição de adsorção de extrato bruto não clarificado de C-ficocianina que se mostrou mais vantajosa, por apresentar maior capacidade de adsorção, é quando se alimenta o extrato até atingir 10% de saturação da resina, em grau de expansão 2,0, com uma altura inicial de leito de 30 cm. O último artigo originado nesta tese foi sobre a purificação de C-ficocianina através da cromatografia de troca iônica em leito expandido. Uma vez que a adsorção já havia sido estudada no artigo 2, o artigo 4 enfoca na otimização das condições de eluição, visando obter um produto com máxima pureza e recuperação. A pureza é dada pela razão entre a absorbância a 620 nm pela absorbância a 280 nm, e dizse que quando C-ficocianina apresenta pureza superior a 0,7 ela pode ser usada em como corante em alimentos. A avaliação das curvas de contorno indicou que a faixa de trabalho deve ser em pH ao redor de 6,5 e volumes de eluição próximos a 150 mL. Tais condições combinadas a uma etapa de pré-eluição com 0,1M de NaCl, permitiu obter C-ficocianina com pureza de 2,9, concentração 3 mg/mL, e recuperação ao redor de 70%.

Post-combustion CO2 adsorption on activated carbons with different textural properties

Fixed bed CO2 adsorption tests were carried out in model flue-gas streams onto two commercial activated carbons, namely Filtrasorb 400 and Nuchar RGC30, at 303 K, 323 K and 353 K. Thermodynamic adsorption results highlighted that the presence of a narrower micropore size distribution with a prevailing contribution of very small pore diameters, observed for Filtrasorb 400, is a key factor in determining a higher CO2 capture capacity, mostly at low temperature. These experimental evidences were also corroborated by the higher value of the isosteric heat derived for Filtrasorb 400, testifying stronger interactions with CO2 molecules with respect to Nuchar RGC30. Dynamic adsorption results on the investigated sorbents highlighted the important role played by both a greater contribution of mesopores and the presence of wider micropores for Nuchar RGC30 in establishing faster capture kinetics with respect to Filtrasorb 400, in particular at 303 K. Furthermore, the modeling analysis of 15% CO2 breakthrough curves allowed identifying intraparticle diffusion as the rate-determining step of the process.

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

Effects of Soil Compaction and Organic Carbon Content on Preferential Flow in Loamy Field Soils

Abstract: Preferential flow and transport through macropores affect plant water use efficiency and enhance leaching of agrochemicals and the transport of colloids, thereby increasing the risk for contamination of groundwater resources. The effects of soil compaction, expressed in terms of bulk density (BD), and organic carbon (OC) content on preferential flow and transport were investigated using 150 undisturbed soil cores sampled from 15 × 15–m grids on two field sites. Both fields had loamy textures, but one site had significantly higher OC content. Leaching experiments were conducted in each core by applying a constant irrigation rate of 10 mm h−1 with a pulse application of tritium tracer. Five percent tritium mass arrival times and apparent dispersivities were derived from each of the tracer breakthrough curves and correlated with texture, OC content, and BD to assess the spatial distribution of preferential flow and transport across the investigated fields. Soils from both fields showed strong positive correlations between BD and preferential flow. Interestingly, the relationships between BD and tracer transport characteristics were markedly different for the two fields, although the relationship between BD and macroporosity was nearly identical. The difference was likely caused by the higher contents of fines and OC at one of the fields leading to stronger aggregation, smaller matrix permeability, and a more pronounced pipe-like pore system with well-aligned macropores.

Geological structure, recharge processes and underground drainage of a glacierised karst aquifer system, Tsanfleuron-Sanetsch, Swiss Alps

The relationships between stratigraphic and tectonic setting, recharge processes and underground drainage of the glacierised karst aquifer system `Tsanfleuron-Sanetsch' in the Swiss Alps have been studied by means of various methods, particularly tracer tests (19 injections). The area belongs to the Helvetic nappes and consists of Jurassic to Palaeogene sedimentary rocks. Strata are folded and form a regional anticlinorium. Cretaceous Urgonian limestone constitutes the main karst aquifer, overlain by a retreating glacier in its upper part. Polished limestone surfaces are exposed between the glacier front and the end moraine of 1855/1860 (Little Ice Age); typical alpine karrenfields can be observed further below. Results show that (1) large parts of the area are drained by the Glarey spring, which is used as a drinking water source, while marginal parts belong to the catchments of other springs; (2) groundwater flow towards the Glarey spring occurs in the main aquifer, parallel to stratification, while flow towards another spring crosses the entire stratigraphic sequence, consisting of about 800 m of marl and limestone, along deep faults that were probably enlarged by mass movements; (3) the variability of glacial meltwater production influences the shape of the tracer breakthrough curves and, consequently, flow and transport in the aquifer.