Numerical modeling of the effect of variation of boundary conditions on vadose zone hydraulic properties

An accurate estimation of hydraulic fluxes in the vadose zone is essential for the prediction of water, nutrient and contaminant transport in natural systems. The objective of this study was to simulate the effect of variation of boundary conditions on the estimation of hydraulic properties (i.e. water content, effective unsaturated hydraulic conductivity and hydraulic flux) in a one-dimensional unsaturated flow model domain. Unsaturated one-dimensional vertical water flow was simulated in a pure phase clay loam profile and in clay loam interlayered with silt loam distributed according to the third iteration of the Cantor Bar fractal object Simulations were performed using the numerical model Hydrus 1D. The upper and lower pressure heads were varied around average values of -55 cm for the near-saturation range. This resulted in combinations for the upper and lower constant head boundary conditions, respectively, of -50 and -60 cm, -40 and -70 cm, -30 and -80 cm, -20 and -90 cm, and -10 and -100 cm. For the drier range the average head between the upper and lower boundary conditions was set to -550 cm, resulting in the combinations -500 and -600 cm, -400 and -700 cm, -300 and -800 cm, -200 and -900 cm, and -100 and -1,000 cm, for upper and lower boundary conditions, respectively. There was an increase in water contents, fluxes and hydraulic conductivities with the increase in head difference between boundary conditions. Variation in boundary conditions in the pure phase and interlayered one-dimensional profiles caused significant deviations in fluxes, water contents and hydraulic conductivities compared to the simplest case (a head difference between the upper and lower constant head boundaries of 10 cm in the wetter range and 100 cm in the drier range).

Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities.

Quantifying the spatial configuration of hydraulic conductivity (K) in heterogeneous geological environments is essential for accurate predictions of contaminant transport, but is difficult because of the inherent limitations in resolution and coverage associated with traditional hydrological measurements. To address this issue, we consider crosshole and surface-based electrical resistivity geophysical measurements, collected in time during a saline tracer experiment. We use a Bayesian Markov-chain-Monte-Carlo (McMC) methodology to jointly invert the dynamic resistivity data, together with borehole tracer concentration data, to generate multiple posterior realizations of K that are consistent with all available information. We do this within a coupled inversion framework, whereby the geophysical and hydrological forward models are linked through an uncertain relationship between electrical resistivity and concentration. To minimize computational expense, a facies-based subsurface parameterization is developed. The Bayesian-McMC methodology allows us to explore the potential benefits of including the geophysical data into the inverse problem by examining their effect on our ability to identify fast flowpaths in the subsurface, and their impact on hydrological prediction uncertainty. Using a complex, geostatistically generated, two-dimensional numerical example representative of a fluvial environment, we demonstrate that flow model calibration is improved and prediction error is decreased when the electrical resistivity data are included. The worth of the geophysical data is found to be greatest for long spatial correlation lengths of subsurface heterogeneity with respect to wellbore separation, where flow and transport are largely controlled by highly connected flowpaths.

Risk Assessment of Modern Landfill Structures in Finland

The purpose of this thesis was to investigate environmental permits of landfills with respect to the appropriateness of risk assessments focusing on contaminant migration, structures capable to protect the environment, waste and leachate management and existing environmental impacts of landfills. According to the requirements, a risk assessment is always required to demonstrate compliance with environmental protection requirements if the environmental permit decision deviates from the set requirements. However, there is a reason to doubt that all relevant risk factors are identified in current risk assessment practices in order to protect people end environment. In this dissertation, risk factors were recognized in 12 randomly selected landfills. Based on this analysis, a structural risk assessment method was created. The method was verified with two case examples. Several development needs were found in the risk assessments of the environmental permit decisions. The risk analysis equations used in the decisions did not adequately take into account all the determining factors like waste prospects, total risk quantification or human delineated factors. Instead of focusing on crucial factors, the landfill environmental protection capability is simply expressed via technical factors like hydraulic conductivity. In this thesis, it could be shown, that using adequate risk assessment approaches the most essential environmental impacts can be taken into account by consideration of contaminant transport mechanisms, leachate effects, and artificial landfill structures. The developed structural risk analysing (SRA) method shows, that landfills structures could be designed in a more cost-efficient way taking advantage of recycled or by-products. Additionally, the research results demonstrate that the environmental protection requirements of landfills should be updated to correspond to the capability to protect the environment instead of the current simplified requirements related to advective transport only.

Experimentelle und numerische Untersuchungen zur Dispersion von Dichteströmungen in einem stochastischen Modellaquifer

Es ist bekannt, dass die Dichte eines gelösten Stoffes die Richtung und die Stärke seiner Bewegung im Untergrund entscheidend bestimmen kann. Eine Vielzahl von Untersuchungen hat gezeigt, dass die Verteilung der Durchlässigkeiten eines porösen Mediums diese Dichteffekte verstärken oder abmindern kann. Wie sich dieser gekoppelte Effekt auf die Vermischung zweier Fluide auswirkt, wurde in dieser Arbeit untersucht und dabei das experimentelle sowohl mit dem numerischen als auch mit dem analytischen Modell gekoppelt. Die auf der Störungstheorie basierende stochastische Theorie der macrodispersion wurde in dieser Arbeit für den Fall der transversalen Makodispersion. Für den Fall einer stabilen Schichtung wurde in einem Modelltank (10m x 1.2m x 0.1m) der Universität Kassel eine Serie sorgfältig kontrollierter zweidimensionaler Experimente an einem stochastisch heterogenen Modellaquifer durchgeführt. Es wurden Versuchsreihen mit variierenden Konzentrationsdifferenzen (250 ppm bis 100 000 ppm) und Strömungsgeschwindigkeiten (u = 1 m/ d bis 8 m/d) an drei verschieden anisotrop gepackten porösen Medien mit variierender Varianzen und Korrelationen der lognormal verteilten Permeabilitäten durchgeführt. Die stationäre räumliche Konzentrationsausbreitung der sich ausbreitenden Salzwasserfahne wurde anhand der Leitfähigkeit gemessen und aus der Höhendifferenz des 84- und 16-prozentigen relativen Konzentrationsdurchgang die Dispersion berechnet. Parallel dazu wurde ein numerisches Modell mit dem dichteabhängigen Finite-Elemente-Strömungs- und Transport-Programm SUTRA aufgestellt. Mit dem kalibrierten numerischen Modell wurden Prognosen für mögliche Transportszenarien, Sensitivitätsanalysen und stochastische Simulationen nach der Monte-Carlo-Methode durchgeführt. Die Einstellung der Strömungsgeschwindigkeit erfolgte - sowohl im experimentellen als auch im numerischen Modell - über konstante Druckränder an den Ein- und Auslauftanks. Dabei zeigte sich eine starke Sensitivität der räumlichen Konzentrationsausbreitung hinsichtlich lokaler Druckvariationen. Die Untersuchungen ergaben, dass sich die Konzentrationsfahne mit steigendem Abstand von der Einströmkante wellenförmig einem effektiven Wert annähert, aus dem die Makrodispersivität ermittelt werden kann. Dabei zeigten sich sichtbare nichtergodische Effekte, d.h. starke Abweichungen in den zweiten räumlichen Momenten der Konzentrationsverteilung der deterministischen Experimente von den Erwartungswerten aus der stochastischen Theorie. Die transversale Makrodispersivität stieg proportional zur Varianz und Korrelation der lognormalen Permeabilitätsverteilung und umgekehrt proportional zur Strömungsgeschwindigkeit und Dichtedifferenz zweier Fluide. Aus dem von Welty et al. [2003] mittels Störungstheorie entwickelten dichteabhängigen Makrodispersionstensor konnte in dieser Arbeit die stochastische Formel für die transversale Makrodispersion weiter entwickelt und - sowohl experimentell als auch numerisch - verifiziert werden.

The use of earth observation techniques to improve catchment-scale pollution predictions

Remote sensing can potentially provide information useful in improving pollution transport modelling in agricultural catchments. Realisation of this potential will depend on the availability of the raw data, development of information extraction techniques, and the impact of the assimilation of the derived information into models. High spatial resolution hyperspectral imagery of a farm near Hereford, UK is analysed. A technique is described to automatically identify the soil and vegetation endmembers within a field, enabling vegetation fractional cover estimation. The aerially-acquired laser altimetry is used to produce digital elevation models of the site. At the subfield scale the hypothesis that higher resolution topography will make a substantial difference to contaminant transport is tested using the AGricultural Non-Point Source (AGNPS) model. Slope aspect and direction information are extracted from the topography at different resolutions to study the effects on soil erosion, deposition, runoff and nutrient losses. Field-scale models are often used to model drainage water, nitrate and runoff/sediment loss, but the demanding input data requirements make scaling up to catchment level difficult. By determining the input range of spatial variables gathered from EO data, and comparing the response of models to the range of variation measured, the critical model inputs can be identified. Response surfaces to variation in these inputs constrain uncertainty in model predictions and are presented. Although optical earth observation analysis can provide fractional vegetation cover, cloud cover and semi-random weather patterns can hinder data acquisition in Northern Europe. A Spring and Autumn cloud cover analysis is carried out over seven UK sites close to agricultural districts, using historic satellite image metadata, climate modelling and historic ground weather observations. Results are assessed in terms of probability of acquisition probability and implications for future earth observation missions. (C) 2003 Elsevier Ltd. All rights reserved.

Estudo da circulação hidrodinâmica e da dispersão de óleo na zona costeira entre Macau e Galinhos RN

This paper presents a study on coastal hydrodynamics and the spread of an oil spill in waters off Macau and Galinhos, on the east coast of the state of Rio Grande do Norte in Northeast Brazil. This area has a very marked coastal dynamic owing to the complexity of its geomorphological features, developed in a regime of semidiurnal mesotides involving reefs, spits, estuaries, mangroves, lakes and dunes. The region also plays an important role in the socioeconomic development of the state, given that the production of oil, natural gas, salt and shrimp is concentrated there. The series of oil platforms is interconnected by a pipeline system that carries oil to the local terminal. This pipeline could leak at any moment, causing immense ecological damage. To gauge the risks of an oil leak and resulting contamination of the coastal region, two hydrodynamic scenarios were simulated. The results obtained were used to implement a contaminant transport model with the creation of various oil leak scenarios modeled at different volumes (from small to large) and intensities (sporadic and continuous), at points considered critical for the model (on two platforms and at two pipeline intersections), under different wind (summer and winter) and tidal (high and low at new, full and quarter moon phases) conditions. The use of hydrodynamic circulation computer models as a tool for representing a real project design has been increasingly frequent in recent years, given that they enable the realistic simulation of the hydrodynamic circulation pattern in bodies of water and an analysis of the impacts caused by contaminants released into the water. This study used the computer models contained in SisBAHIA®, in continuous development in the area of Coastal Engineering and Oceanography at COPPE/UFRJ

Estimativa de alguns parâmetros de transporte de contaminantes para solos Brasileiros

This paper introduces an indirect estimate for the coefficients of distribution, hydrodynamic dispersion and retardation for contaminants commonly encountered in sanitary landfills and their liners, such as Cu2+ and K+; this estimate is based on the relationship between concentration and certain physical characteristics of typical Brazilian soils. The results of previous studies investigating the migration of contaminants were used to develop mathematical expressions from multiple non-linear regressions. Using minimal squares regression, this transport was linked to various combinations of contaminant concentration and both structural and textural characteristics of the porous medium. Various combinations of characteristics and concentrations were investigated, with a mathematical expression obtained for each. The relationship between percentage of clay and the contaminant content proved to be the most closely correlated with actual transport parameters, with coefficients close to one.

Tecniche avanzate per l'analisi di sensitività globale: applicazione alla contaminazione delle acque sotterranee.

The thesis is framed within the field of the stochastic approach to flow and transport themes of solutes in natural porous materials. The methodology used to characterise the uncertainty associated with the modular predictions is completely general and can be reproduced in various contexts. The theme of the research includes the following among its main objectives: (a) the development of a Global Sensitivity Analysis on contaminant transport models in the subsoil to research the effects of the uncertainty of the most important parameters; (b) the application of advanced techniques, such as Polynomial Chaos Expansion (PCE), for obtaining surrogate models starting from those which conduct traditionally developed analyses in the context of Monte Carlo simulations, characterised by an often not negligible computational burden; (c) the analyses and the understanding of the key processes at the basis of the transport of solutes in natural porous materials using the aforementioned technical and analysis resources. In the complete picture, the thesis looks at the application of a Continuous Injection transport model of contaminants, of the PCE technique which has already been developed and applied by the thesis supervisors, by way of numerical code, to a Slug Injection model. The methodology was applied to the aforementioned model with original contribution deriving from surrogate models with various degrees of approximation and developing a Global Sensitivity Analysis aimed at the determination of Sobol’ indices.

Characterization of granite matrix porosity and pore-space geometry by in situ and laboratory methods

Most available studies of interconnected matrix porosity of crystalline rocks are based on laboratory investigations; that is, work on samples that have undergone stress relaxation and were affected by drilling and sample preparation. The extrapolation of the results to in situ conditions is therefore associated with considerable uncertainty, and this was the motivation to conduct the ‘in situ Connected Porosity’ experiment at the Grimsel Test Site (Central Swiss Alps). An acrylic resin doped with fluorescent agents was used to impregnate the microporous granitic matrix in situ around an injection borehole, and samples were obtained by overcoring. The 3-D structure of the porespace, represented by microcracks, was studied by U-stage fluorescence microscopy. Petrophysical methods, including the determination of porosity, permeability and P -wave velocity, were also applied. Investigations were conducted both on samples that were impregnated in situ and on non-impregnated samples, so that natural features could be distinguished from artefacts. The investigated deformed granites display complex microcrack populations representing a polyphase deformation at varying conditions. The crack population is dominated by open cleavage cracks in mica and grain boundary cracks. The porosity of non-impregnated samples lies slightly above 1 per cent, which is 2–2.5 times higher than the in situ porosity obtained for impregnated samples. Measurements of seismic velocities (Vp ) on spherical rock samples as a function of confining pressure, spatial direction and water saturation for both non-impregnated and impregnated samples provide further constraints on the distinction between natural and induced crack types. The main conclusions are that (1) an interconnected network of microcracks exists in the whole granitic matrix, irrespective of the distance to ductile and brittle shear zones, and (2) conventional laboratory methods overestimate the matrix porosity. Calculations of contaminant transport through fractured media often rely on matrix diffusion as a retardation mechanism.

Proceedings of the 22nd Annual Biochemical Engineering Symposium

This is the twenty-second of a series of symposia devoted to talks and posters by students about their biochemical engineering research. The first, third, fifth, ninth, twelfth, sixteenth, and twenti~th were hosted by Kansas State University, the second and fourth by the University of Nebraska- Lincoln, the sixth, seventh, tenth, thirteenth, seventeenth, and twenty-second by Iowa State University, the eighth, fourteenth, and nineteenth by the University of Missouri-Columbia, the eleventh, fifteenth, and twenty-first by Colorado State University, and the eighteenth by the University of Colorado. Next year's symposium will be at the University of Oklahoma. Symposium proceedings are edited and issued by faculty of the host institution. Because final publication usually takes place in refereed journals, articles included here are brief and often cover work in progress. ContentsC. A. Baldwin, J.P. McDonald, and L. E. Erickson, Kansas State University. Effect of Hydrocarbon Phase on Kinetic and Transport Limitations for Bioremediation of Microporous Soil J. C. Wang, S. K. Banerji, and Rakesh Bajpai, University of Missouri-Columbia. Migration of PCP in Soil-Columns in Presence of a Second Organic Phase Cheng-Hsien Hsu and Roger G. Harrison, University of Oklahoma. Bacterial Leaching of Zinc and Copper from Mining Wastes James A. Searles, Paul Todd, and Dhinakar S. Kompala, University of Colorado. Suspension Culture of Chinese Hamster Ovary Cells Utilizing Inclined Sedimentation Ron Beyerinck and Eric H. Dunlop, Colorado State University. The Effect of Feed Zone Turbulence as Measured by Laser Doppler Velocimetry on Baker's Yeast Metabolism in a Chemostat Paul Li-Hong Yeh, GraceY. Sun, Gary A. Weisman, and Rakesh Bajpai, University of Missouri-Columbia. Effect of Medium Constituents upon Membrane Composition of Insect Cells R. Shane Gold, M. M. Meagher, R. Hutkins, and T. Conway, University of Nebraska-Lincoin. Ethanol Tolerance and Carbohydrate Metabolism in Lactobacilli John Sargantanis and M. N. Karim, Colorado State University. Application of Kalman Filter and Adaptive Control in Solid Substrate Fermentation D. Vrana, M. Meagher, and R. Hutkins, University of Nebraska-Lincoln. Product Recovery Optimization in the ABE Fermentation Kalyan R. Tadikonda and Robert H. Davis, University of Colorado. Cell Separations Using Targeted Monoclonal Antibodies Against Surface Proteins Meng H. Heng and Charles E. Glatz, Iowa State University. Charged Fusion for Selective Recovery of B-Galactosidase from Cell Extract Using Hollow Fiber Ion-Exchange Membrane Adsorption Hsiu-Mei Chen, Peter J. Reilly, and Clark Ford, Iowa State University. Site-Directed Mutagenesis to Enhance Thermostability of Glucoamylase from Aspergillus: A Rational Approach P. Tuitemwong, L. E. Erickson, and D. Y. C. Fung, Kansas State University. Applications of Enzymatic Hydrolysis and Fermentation on the Reduction of Flatulent Sugars in the Rapid Hydration Hydrothermal Cooked Soy Milk Sanjeev Redkar and Robert H. Davis, University of Colorado. Crossflow Microfiltration of Yeast Suspensions Linda Henk and James C. Linden, Colorado State University, and Irving C. Anderson, Iowa State University. Evaluation of Sorghum Ensilage as an Ethanol Feedstock Marc Lipovitch and James C. Linden, Colorado State University. Stability and Biomass Feedstock Pretreatability for Simultaneous Saccharification and Fermentation Ali Demirci, Anthony L. Pometto Ill, and Kenneth E. Johnson, Iowa State University. Application of Biofilm Reactors in Lactic Acid Fermentation Michael K. Dowd, Peter I. Reilly, and WalterS. Trahanovsky, Iowa State University. Low Molecular-Weight Organic Composition of Ethanol Stillage from Corn Craig E. Forney, Meng H. Heng, John R. Luther, Mark Q. Niederauer, and Charles E. Glatz, Iowa State University. Enhancement of Protein Separation Using Genetic Engineering J. F. Shimp, J. C. Tracy, E. Lee, L. C. Davis, and L. E. Erickson, Kansas State University. Modeling Contaminant Transport, Biodegradation and Uptake by Plants in the Rhizosphere Xiaoqing Yang, L. E. Erickson, and L. T. Fan, Kansas State University. Modeling of Dispersive-Convective Characteristics in Bioremediation of Contaminated Soil Jan Johansson and Rakesh Bajpai, University of Missouri-Columbia. Fouling of Membranes J. M. Wang, S. K. Banerji, and R. K. Bajpai, University of Missouri-Columbia. Migration of Sodium-Pentachorophenol (Na-PCP) in Unsaturated and Saturated Soil-Columns J. Sweeney and M. Meagher, University of Nebraska-Lincoln. The Purification of Alpha-D-Glucuronidase from Trichoderma reesei

(Table 1) Concentration of PCB and other contaminants in blood plasma of polar bears (Ursus maritimus) from Svalbard

Persistent chemicals accumulate in the arctic environment due to their chemical reactivity and physicochemical properties and polychlorinated biphenyls (PCBs) are the most concentrated pollutant class in polar bears (Ursus maritimus). Metabolism of PCB and polybrominated biphenyl ether (PBDE) flame-retardants alter their toxicological properties and these metabolites are known to interfere with the binding of thyroid hormone (TH) to transthyretin (TTR) in rodents and humans. In polar bear plasma samples no binding of [125I]-T4 to TTR was observed after incubation and PAGE separation. Incubation of the plasma samples with [14C]-4-OH-CB107, a compound with a higher binding affinity to TTR than the endogenous ligand T4 resulted in competitive binding as proven by the appearance of a radio labeled TTR peak in the gel. Plasma incubation with T4 up to 1 mM, a concentration that is not physiologically relevant anymore did not result in any visible competition. These results give evidence that the binding sites on TTR for T4 in wild living polar bears are completely saturated. Such saturation of binding sites can explain observed lowered levels of THs and could lead to contaminant transport into the developing fetus.

Método multigrid algébrico: reutilização das estruturas multigrid no transporte de contaminantes

A necessidade de obter solução de grandes sistemas lineares resultantes de processos de discretização de equações diferenciais parciais provenientes da modelagem de diferentes fenômenos físicos conduz à busca de técnicas numéricas escaláveis. Métodos multigrid são classificados como algoritmos escaláveis.Um estimador de erros deve estar associado à solução numérica do problema discreto de modo a propiciar a adequada avaliação da solução obtida pelo processo de aproximação. Nesse contexto, a presente tese caracteriza-se pela proposta de reutilização das estruturas matriciais hierárquicas de operadores de transferência e restrição dos métodos multigrid algébricos para acelerar o tempo de solução dos sistemas lineares associados à equação do transporte de contaminantes em meio poroso saturado. Adicionalmente, caracteriza-se pela implementação das estimativas residuais para os problemas que envolvem dados constantes ou não constantes, os regimes de pequena ou grande advecção e pela proposta de utilização das estimativas residuais associadas ao termo de fonte e à condição inicial para construir procedimentos adaptativos para os dados do problema. O desenvolvimento dos códigos do método de elementos finitos, do estimador residual e dos procedimentos adaptativos foram baseados no projeto FEniCS, utilizando a linguagem de programação PYTHONR e desenvolvidos na plataforma Eclipse. A implementação dos métodos multigrid algébricos com reutilização considera a biblioteca PyAMG. Baseado na reutilização das estruturas hierárquicas, os métodos multigrid com reutilização com parâmetro fixo e automática são propostos, e esses conceitos são estendidos para os métodos iterativos não-estacionários tais como GMRES e BICGSTAB. Os resultados numéricos mostraram que o estimador residual captura o comportamento do erro real da solução numérica, e fornece algoritmos adaptativos para os dados cuja malha retornada produz uma solução numérica similar à uma malha uniforme com mais elementos. Adicionalmente, os métodos com reutilização são mais rápidos que os métodos que não empregam o processo de reutilização de estruturas. Além disso, a eficiência dos métodos com reutilização também pode ser observada na solução do problema auxiliar, o qual é necessário para obtenção das estimativas residuais para o regime de grande advecção. Esses resultados englobam tanto os métodos multigrid algébricos do tipo SA quanto os métodos pré-condicionados por métodos multigrid algébrico SA, e envolvem o transporte de contaminantes em regime de pequena e grande advecção, malhas estruturadas e não estruturadas, problemas bidimensionais, problemas tridimensionais e domínios com diferentes escalas.

Quantitative characterization of physico-chemical properties of the active layers of reverse osmosis and nanofiltration membranes, and their relation to membrane performance

The main objectives of this dissertation were: (i) to develop experimental and analytical procedures to quantify different physico-chemical properties of the ultra-thin (~ 100 nm) active layers of reverse osmosis (RO) and nanofiltration (NF) membranes and their interactions with contaminants; (ii) to use such procedures to evaluate the similarities and differences between the active layers of different RO/NF membranes; and (iii) to relate characterization results to membrane performance. Such objectives were motivated by the current limited understanding of the physico-chemical properties of active layers as a result of traditional characterization techniques having limitations associated with the nanometer-scale spatial resolution required to study these ultra-thin films. Functional groups were chosen as the main active layer property of interest. Specific accomplishments of this study include the development of procedures to quantify in active layers as a function of pH: (1) the concentration of both negatively and positively ionized functional groups; (2) the stoichiometry of association between ions (i.e., barium) and ionized functional groups (i.e., carboxylate and sulfonate); and (3) the steric effects experienced by ions (i.e., barium). Conceptual and mathematical models were developed to describe experimental results. The depth heterogeneity of the active layer physico-chemical properties and interactions with contaminants studied in this dissertation was also characterized. Additionally, measured concentrations of ionized functional groups in the polyamide active layers of several commercial RO/NF membranes were used as input in a simplified RO/NF transport model to predict the rejection of a strong electrolyte (i.e., potassium iodide) and a weak acid (i.e., arsenious acid) at different pH values based on rejection results at one pH condition. The good agreement between predicted and experimental results showed that the characterization procedures developed in this study serve as useful tools in the advancement of the understanding of the properties and structure of the active layers of RO/NF membranes, and the mechanisms of contaminant transport through them.

Solução semi-analítica para modelagem das condições inicial e de contorno aplicáveis ao ensaio de difusão pura

Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Civil e Ambiental, 2015.

BeachWin: modelling groundwater effects on swash sediment transport and beach profile changes

Field and laboratory observations have shown that a relatively low beach groundwater table enhances beach accretion. These observations have led to the beach dewatering technique (artificially lowering the beach water table) for combating beach erosion. Here we present a process-based numerical model that simulates the interacting wave motion on the beach. coastal groundwater flow, swash sediment transport and beach profile changes. Results of model simulations demonstrate that the model replicates accretionary effects of a low beach water table on beach profile changes and has the potential to become a tool for assessing the effectiveness of beach dewatering systems. (C) 2002 Elsevier Science Ltd. All rights reserved.