Process-based modelling of mass transport in the subsurface

Modern societies rely on natural water pathways that include subsurface flow of water and dissolved chemicals. The thesis presents a range of numerical and analytical models for simulating physical, chemical and biological processes in the subsurface, including coastal aquifers, the near-surface vadose zone, and solute transport in biogeochemically active aquifers.

A three dimensional numerical investigation of hillslope flow processes

Physically based distributed models of catchment hydrology are likely to be made available as engineering tools in the near future. Although these models are based on theoretically acceptable equations of continuity, there are still limitations in the present modelling strategy. Of interest to this thesis are the current modelling assumptions made concerning the effects of soil spatial variability, including formations producing distinct zones of preferential flow. The thesis contains a review of current physically based modelling strategies and a field based assessment of soil spatial variability. In order to investigate the effects of soil nonuniformity a fully three dimensional model of variability saturated flow in porous media is developed. The model is based on a Galerkin finite element approximation to Richards equation. Accessibility to a vector processor permits numerical solutions on grids containing several thousand node points. The model is applied to a single hillslope segment under various degrees of soil spatial variability. Such variability is introduced by generating random fields of saturated hydraulic conductivity using the turning bands method. Similar experiments are performed under conditions of preferred soil moisture movement. The results show that the influence of soil variability on subsurface flow may be less significant than suggested in the literature, due to the integrating effects of three dimensional flow. Under conditions of widespread infiltration excess runoff, the results indicate a greater significance of soil nonuniformity. The recognition of zones of preferential flow is also shown to be an important factor in accurate rainfall-runoff modelling. Using the results of various fields of soil variability, experiments are carried out to assess the validity of the commonly used concept of `effective parameters'. The results of these experiments suggest that such a concept may be valid in modelling subsurface flow. However, the effective parameter is observed to be event dependent when the dominating mechanism is infiltration excess runoff.

Complementary methods to investigate the development of clogging within a horizontal sub-surface flow tertiary treatment wetland

A combination of experimental methods was applied at a clogged, horizontal subsurface flow (HSSF) municipal wastewater tertiary treatment wetland (TW) in the UK, to quantify the extent of surface and subsurface clogging which had resulted in undesirable surface flow. The three dimensional hydraulic conductivity profile was determined, using a purpose made device which recreates the constant head permeameter test in-situ. The hydrodynamic pathways were investigated by performing dye tracing tests with Rhodamine WT and a novel multi-channel, data-logging, flow through Fluorimeter which allows synchronous measurements to be taken from a matrix of sampling points. Hydraulic conductivity varied in all planes, with the lowest measurement of 0.1 md1 corresponding to the surface layer at the inlet, and the maximum measurement of 1550 md1 located at a 0.4m depth at the outlet. According to dye tracing results, the region where the overland flow ceased received five times the average flow, which then vertically short-circuited below the rhizosphere. The tracer break-through curve obtained from the outlet showed that this preferential flow-path accounted for approximately 80% of the flow overall and arrived 8 h before a distinctly separate secondary flow-path. The overall volumetric efficiencyof the clogged system was 71% and the hydrology was simulated using a dual-path, dead-zone storage model. It is concluded that uneven inlet distribution, continuous surface loading and high rhizosphere resistance is responsible for the clog formation observed in this system. The average inlet hydraulic conductivity was 2 md1, suggesting that current European design guidelines, which predict that the system will reach an equilibrium hydraulic conductivity of 86 md1, do not adequately describe the hydrology of mature systems.

Isotopic fingerprints in surficial waters : Stable isotope methods applied in hydrogeological studies

The driving force behind this study has been the need to develop and apply methods for investigating the hydrogeochemical processes of significance to water management and artificial groundwater recharge. Isotope partitioning of elements in the course of physicochemical processes produces isotopic variations to their natural reservoirs. Tracer property of the stable isotope abundances of oxygen, hydrogen and carbon has been applied to investigate hydrogeological processes in Finland. The work described here has initiated the use of stable isotope methods to achieve a better understanding of these processes in the shallow glacigenic formations of Finland. In addition, the regional precipitation and groundwater records will supplement the data of global precipitation, but as importantly, provide primary background data for hydrological studies. The isotopic composition of oxygen and hydrogen in Finnish groundwaters and atmospheric precipitation was determined in water samples collected during 1995 2005. Prior to this study, no detailed records existed on the spatial or annual variability of the isotopic composition of precipitation or groundwaters in Finland. Groundwaters and precipitation in Finland display a distinct spatial distribution of the isotopic ratios of oxygen and hydrogen. The depletion of the heavier isotopes as a function of increasing latitude is closely related to the local mean surface temperature. No significant differences were observed between the mean annual isotope ratios of oxygen and hydrogen in precipitation and those in local groundwaters. These results suggest that the link between the spatial variability in the isotopic composition of precipitation and local temperature is preserved in groundwaters. Artificial groundwater recharge to glaciogenic sedimentary formations offers many possibilities to apply the isotopic ratios of oxygen, hydrogen and carbon as natural isotopic tracers. In this study the systematics of dissolved carbon have been investigated in two geochemically different glacigenic groundwater formations: a typical esker aquifer at Tuusula, in southern Finland and a carbonate-bearing aquifer with a complex internal structure at Virttaankangas, in southwest Finland. Reducing the concentration of dissolved organic carbon (DOC) in water is a primary challenge in the process of artificial groundwater recharge. The carbon isotope method was used to as a tool to trace the role of redox processes in the decomposition of DOC. At the Tuusula site, artificial recharge leads to a significant decrease in the organic matter content of the infiltrated water. In total, 81% of the initial DOC present in the infiltrated water was removed in three successive stages of subsurface processes. Three distinct processes in the reduction of the DOC content were traced: The decomposition of dissolved organic carbon in the first stage of subsurface flow appeared to be the most significant part in DOC removal, whereas further decrease in DOC has been attributed to adsorption and finally to dilution with local groundwater. Here, isotope methods were used for the first time to quantify the processes of DOC removal in an artificial groundwater recharge. Groundwaters in the Virttaankangas aquifer are characterized by high pH values exceeding 9, which are exceptional for shallow aquifers on glaciated crystalline bedrock. The Virttaankangas sediments were discovered to contain trace amounts of fine grained, dispersed calcite, which has a high tendency to increase the pH of local groundwaters. Understanding the origin of the unusual geochemistry of the Virttaankangas groundwaters is an important issue for constraining the operation of the future artificial groundwater plant. The isotope ratios of oxygen and carbon in sedimentary carbonate minerals have been successfully applied to constrain the origin of the dispersed calcite in the Virttaankangas sediments. The isotopic and chemical characteristics of the groundwater in the distinct units of aquifer were observed to vary depending on the aquifer mineralogy, groundwater residence time and the openness of the system to soil CO2. The high pH values of > 9 have been related to dissolution of calcite into groundwater under closed or nearly closed system conditions relative to soil CO2, at a low partial pressure of CO2.

Construção e operação de wetland como polimento de lixiviado de aterro de resíduos sólidos urbanos.

O presente trabalho tem o objetivo de avaliar a metodologia de construção, de plantio, de manutenção, de operação e a eficiência de um sistema de tratamento por wetland construído como etapa de polimento da Estação de Tratamento de Lixiviado (ETC) de Aterro já encerrado de Resíduos Sólidos Urbanos na Região Metropolitana do Rio de Janeiro. Parte do efluente tratado por lodos ativados na ETC foi direcionado e tratado no wetland construído. Foi escolhido o projeto de Fluxo Horizontal Subsuperficial e a vegetação selecionada foi a taboa (Typha latifolia) que é nativa da área do aterro. Em média, foram feitas três amostragens mensais do afluente e do efluente do wetland, de maio a outubro de 2013. A eficiência do sistema foi avaliada por meio de parâmetros físico-químicos e de parâmetros coletivos específicos. A eficiência de remoção de Demanda Química de Oxigênio (DQO) foi de 60%, nitrogênio amoniacal de 67%, nitrito de 72% e nitrato de 57%. Outro parâmetro avaliado foi a toxicidade aguda, foram utilizados os organismos teste Danio rerio (peixe), a Daphnia similis (microcrustáceo) e a Aliivibrio fischeri (bactéria luminescente). Durante o período foram coletados diariamente as vazões de entrada e saída, a condutividade elétrica e o índice pluviométrico. Os resultados mostraram que o uso de wetland como etapa de polimento pode ser uma alternativa para o tratamento de lixiviado.

水文条件对紫色土坡面土壤侵蚀及养分流失的影响

采用人工模拟降雨试验,研究水文条件对紫色土坡面土壤侵蚀及氮和磷养分流失的影响。试验处理包括2个施肥水平(低肥和高肥水平),4个水文条件(自由下渗、土壤水分饱和、壤中流、壤中流+降雨)和一个降雨强度(60 mm/h,历时60 min)。结果表明:壤中流+降雨和土壤水分饱和条件下的土壤侵蚀量分别是自由下渗条件下的3.1和1.7倍,同自由下渗相比,壤中流、壤中流+降雨和土壤水分饱和条件下,地表径流中NO3-N、HPO4-P的浓度和流失量有显著增加;低肥水平条件下,自由下渗、土壤水分饱和、壤中流和壤中流+降雨地表径流中,NO3-N的浓度分别是0.88、58.90、698.41和87.80 mg/L,对应水文条件下地表径流中,HPO4-P的浓度分别是0.252、0.322、0.811和0.383 mg/L,高肥水平条件下,径流中的NO3-N和HPO4-P的浓度也有相同的趋势;土壤水分饱和条件下,地表径流中NO3-N和HPO4-P的流失量分别是自由下渗条件下的27~39和1.3倍,壤中流+降雨条件下,地表径流中NO3-N和HPO4-P的流失量分别是自由下渗条件下的100~114和1.5~1.7倍,同时,壤中流+降雨和土壤...

森林流域单元水文模型的构建与应用

目前全球缺水、水污染、洪涝灾害以及水土流失仍然非常严重，尤其在我国北方地区。流域水文模型可用来进行不同需水管理的情景分析，为解决我国水问题提供科学依据。分布式水文模型是流域水文模型的发展方向，具有显著特点：1）应用前景广泛，不仅可以模拟流域水文过程，还可以协助模拟泥沙或污染物的运移过程，为水利工程设计、水土保持、环境保护等领域提供技术支持；2）能够预测流域土地利用或气候变化下的流域水文响应过程变化，为管理部门提供决策支持；3）模型所需要的参数全部具有物理意义，可通过实际测量确定，适合模拟实测系列较短或是无观测流域的水文过程；4）对于目前国际水文界的前沿问题—水文尺度转换提供了一种有效的解决途径。 然而分布式水文模型还不完善，如1)真实性问题。对一些水文过程和边界条件还不确定。2)尺度转换问题。目前很少考虑尺度对参数有效性的影响。3)检验问题。还无法判断对有些难以测量的水文状态变量的模拟正确与否。4)计算时间和数据存储的问题。有些分布式水文模型虽然具有很强的水文物理基础和完善的模型结构，但是计算时间过长和(或)数据存储过大，难以应用。上述问题的核心就是对分布式水文模型的核心—单元水文模型的研究不够，需要为进一步完善单元水文模型进行研究。 本文采用饱和入渗理论、Saint-Venant方程、Richards方程、Penman-Monteith方程等等构建了以有限差分法求解的适用于森林流域的单元水文模型，并通过实验室模拟试验和坡地径流场资料进行了验证，主要结论为： 通过不同坡度和不同雨强下的室内坡面产汇流实验模拟，表明：该模型模拟的坡面流和壤中流过程与实测过程基本一致，峰现时间、径流历时、峰值流量、出流总量模拟值与实测值的相对误差均较小，基本小于10%。模型的模拟精度较高，实用性较强，为深入研究壤中流机制和改进流域降雨-径流模型提供了理论依据。 通过坡地径流观测场实测资料的验证，表明：该模型模拟的坡面流过程精度较高，累计流量的精度更高于小时过程的精度，离差系数、效率系数、确定系数均较理想，具有应用价值，有助于改善分布式水文模型在森林流域的模拟效果。

Connectivité hydrologique et signature géochimique à l'échelle événementielle dans un bassin versant forestier

Dans un bassin versant, la connectivité hydrologique constitue un état hydrologique qui lie le versant à la zone riveraine. Ses impacts sur la production du débit et le transfert des éléments dissous vers le cours d’eau sont présumés substantiels. L’étude vise à 1) détecter les hydrotopes et les connexions hydrologiques à l’aide d’un réseau de puits qui permet la mesure des fluctuations de la nappe phréatique (NP); 2) identifier la variabilité spatio-temporelle et la signature géochimique des sources potentielles en eau à l’aide des éléments majeurs et traces et 3) examiner la contribution spatio-temporelle respective des sources en eau du bassin lors d’un événement de précipitation. L’étude s’effectue dans un bassin versant forestier du Bouclier canadien (l’Hermine). Nous démontrons l’existence de quatre hydrotopes représentant un gradient de convergence de l’eau, soulignant la diversité de comportement de NP. Les connexions hydrologiques se caractérisent par des coefficients de Spearman élevés des relations entre la profondeur de la NP et le débit, dans leur partie en aval, et s’enclenchent par le fill and spill. Le comportement de NP est influencé par la distance aux limites du bassin, l’horizonation du sol et la topographie souterraine. En somme, trois sources en eau se connectent à partir du versant vers la zone riveraine durant l’événement pluvial de manière chronologique: 1) les horizons B et la NP de l’ensemble du bassin (Sr); 2) les horizons LFH des zones de convergence (Ba et Zn) et 3) une dépression de sol humide sur le versant nord (Co et Mn).

Phosphorus sorption and availability in soils amended with animal manures and sewage sludge

Soils that receive large applications of animal wastes and sewage sludge are vulnerable to releasing environmentally significant concentrations of dissolved P available to subsurface flow owing to the gradual saturation of the soil's P sorption capacity. This study evaluated P sorption (calculated from Langmuir isotherms) and availability of P (as CaCl2-P and resin P) in soils incubated for 20 d with poultry litter, poultry manure, cattle slurry, municipal sewage sludge, or KH2PO4, added on a P-equivalent basis (100 mg P kg(-1)). All the P sources had a marked negative effect on P sorption and a positive effect on P availability in all soils. In the cattle slurry- and KH2PO4- treated soils, the decreases in P sorption maximum (19-66%) and binding energy (25-89%) were consistently larger than the corresponding decreases (7-41% and 11-30%) in poultry litter-, poultry manure-, and sewage sludge-treated soils. The effects of cattle slurry and KH2PO4 on P availability were, in most cases, larger than those of the other P sources. In the poultry litter, poultry manure, and sewage sludge treatments, the increase in soil solution P was inversely related (R-2 = 0.75) to the input of Ca from these relatively high Ca (13.5-42 g kg(-1)) sources. Correlation analyses implied that the magnitude of the changes in P sorption and availability was not related to the water-extractable P content of the P sources. Future research on the sustainable application of organic wastes to agricultural soils needs to consider the non-P- as well as P-containing components of the waste.

Contribution of N Species and P Fractions to Stream Water Quality in Agricultural Catchments

The contribution from agricultural catchments to stream nitrogen and phosphorus concentrations was assessed by evaluation of the chemical composition of these nutrients in agricultural runoff for both surface and subsurface flow pathways. A range of land uses (grazed and ungrazed grassland, cereals, roots) in intensive agricultural systems was studied at scales from hillslope plots (0.5m2) to large catchment (>300km2). By fractionating the total nutrient load it was possible to establish that most of the phosphorus was transported in the unreactive (particulate and organic) fraction via surface runoff. This was true regardless of the scale of measurement. The form of the nitrogen load varied with land use and grazing intensity. High loads of dissolved inorganic nitrogen (with >90% transported as NH4-N) were recorded in surface runoff from heavily grazed land. In subsurface flow from small (2km2) subcatchments and in larger (>300 km2) catchments, organic nitrogen was found to be an important secondary constituent of the total nitrogen load, comprising 40% of the total annual load.

Modelling the impact of land use change on water quality in agricultural catchments

Export coefficient modelling was used to model the impact of agriculture on nitrogen and phosphorus loading on the surface waters of two contrasting agricultural catchments. The model was originally developed for the Windrush catchment where the highly reactive Jurassic limestone aquifer underlying the catchment is well connected to the surface drainage network, allowing the system to be modelled using uniform export coefficients for each nutrient source in the catchment, regardless of proximity to the surface drainage network. In the Slapton catchment, the hydrological path-ways are dominated by surface and lateral shallow subsurface flow, requiring modification of the export coefficient model to incorporate a distance-decay component in the export coefficients. The modified model was calibrated against observed total nitrogen and total phosphorus loads delivered to Slapton Ley from inflowing streams in its catchment. Sensitivity analysis was conducted to isolate the key controls on nutrient export in the modified model. The model was validated against long-term records of water quality, and was found to be accurate in its predictions and sensitive to both temporal and spatial changes in agricultural practice in the catchment. The model was then used to forecast the potential reduction in nutrient loading on Slapton Ley associated with a range of catchment management strategies. The best practicable environmental option (BPEO) was found to be spatial redistribution of high nutrient export risk sources to areas of the catchment with the greatest intrinsic nutrient retention capacity.

Influence of UV radiation on chlorophyll, and antioxidant enzymes of wetland plants in different types of constructed wetland

A surface- and vertical subsurface-flow-constructed wetland were designed to study the response of chlorophyll and antioxidant enzymes to elevated UV radiation in three types of wetland plants (Canna indica, Phragmites austrail, and Typha augustifolia). Results showed that (1) chlorophyll content of C. indica, P. austrail, and T. augustifolia in the constructed wetland was significantly lower where UV radiation was increased by 10 and 20 % above ambient solar level than in treatment with ambient solar UV radiation (p < 0.05). (2) The malondialdehyde (MDA) content, guaiacol peroxidase (POD), superoxide dismutase (SOD), and catalase (CAT) activities of wetland plants increased with elevated UV radiation intensity. (3) The increased rate of MDA, SOD, POD, and CAT activities of C. indica, P. australis, and T. angustifolia by elevated UV radiation of 10 % was higher in vertical subsurface-flow-constructed wetland than in surface-flow-constructed wetland. The sensitivity of MDA, SOD, POD, and CAT activities of C. indica, P. austrail, and T. augustifolia to the elevated UV radiation was lower in surface-flow-constructed wetland than in the vertical subsurface-flow-constructed wetland, which was related to a reduction in UV radiation intensity through the dissolved organic carbon and suspended matter in the water. C. indica had the highest SOD and POD activities, which implied it is more sensitive to enhanced UV radiation. Therefore, different wetland plants had different antioxidant enzymes by elevated UV radiation, which were more sensitive in vertical subsurface-flow-constructed wetland than in surface-flow-constructed wetland.

Multi‐method global sensitivity analysis (MMGSA) for modelling floodplain hydrological processes

When studying hydrological processes with a numerical model, global sensitivity analysis (GSA) is essential if one is to understand the impact of model parameters and model formulation on results. However, different definitions of sensitivity can lead to a difference in the ranking of importance of the different model factors. Here we combine a fuzzy performance function with different methods of calculating global sensitivity to perform a multi-method global sensitivity analysis (MMGSA). We use an application of a finite element subsurface flow model (ESTEL-2D) on a flood inundation event on a floodplain of the River Severn to illustrate this new methodology. We demonstrate the utility of the method for model understanding and show how the prediction of state variables, such as Darcian velocity vectors, can be affected by such a MMGSA. This paper is a first attempt to use GSA with a numerically intensive hydrological model.

Multi-method global sensitivity analysis (MMGSA) for modelling floodplain hydrological processes

When studying hydrological processes with a numerical model, global sensitivity analysis (GSA) is essential if one is to understand the impact of model parameters and model formulation on results. However, different definitions of sensitivity can lead to a difference in the ranking of importance of the different model factors. Here we combine a fuzzy performance function with different methods of calculating global sensitivity to perform a multi-method global sensitivity analysis (MMGSA). We use an application of a finite element subsurface flow model (ESTEL-2D) on a flood inundation event on a floodplain of the River Severn to illustrate this new methodology. We demonstrate the utility of the method for model understanding and show how the prediction of state variables, such as Darcian velocity vectors, can be affected by such a MMGSA. This paper is a first attempt to use GSA with a numerically intensive hydrological model