The aim of the thesis is to find a method, in conjunction with the ordinary differential equation (ODE) based method of lines (MOL) solution of Richards’ equation, to model the steep wetting front infiltration in very dry soils, accurately and efficiently. Due to the steep pressure head or steep water volumetric content gradients, highly nonlinear soil hydraulic properties and the rapid movement of the wetting front, accurate solutions for infiltration into a dry soil are usually difficult to obtain. Additionally, such problems often require very small time steps and large computation times. As an enhancement to the used ODE/MOL approach, Higher Order Finite Differencing, Varying Order Finite Differencing, Vertical Scaling, Adaptive Schemes and Non-uniform Stretching Techniques have been implemented and tested in this thesis. Success has been found in the ability of Vertical Scaling to simulate very steep moving front solution for the Burgers’ equation. Unfortunately, the results also show that Vertical Scaling needs significant research and improvement before their full potential in routine applications for difficult nonlinear problems, such as Richard’s equation with very steep moving front solution, can be realized. However, we have also shown that the use of the composed form of RE and a 2nd order finite differencing for the first order derivative approximation is conducive for modelling steep moving front problem in a very dry soil. Additionally, with the combination of an optimal influx value at the edges of the inlet, the ODE/MOL approach is able to model a 2-D infiltration in very dry soils, effectively and accurately. Furthermore, one of the strengths of this thesis is the use of a MATLAB PDE template. Implementing the ODE/MOL approach via a MATLAB PDE template has shown to be most suitable for modelling of partial differential equations. The plug and play mode of modifying the PDE template for solving time-dependent partial differential equations is user-friendly and easy, as compared to more conventional approaches using Pascal, Fortran, C or C++. The template offers greater modularity, flexibility, versatility, and efficiency for solving PDE problems in both 1-D and 2-D spatial dimensions. Moreover, the 2-D PDE template has been extended for irregular shaped domains.

<p>For many surface waters, phosphorus (P) leaching is a serious problem that should be minimized to prevent eutrophication. In Sweden there is a demand for physical and technical development of high-performance P removal techniques to reduce phosphorus leaching from on-site wastewater treatment systems to the Baltic Sea. However, although these systems are designed to reduce eutrophication there are also other environmental impacts to be considered when implementing them in on-site systems; energy use and global warming potential are two examples. This study has investigated several bed filter materials (reactive media and natural soils) for their total environmental impact (in commercial applications) as well as for the predominating chemical phosphorus removal mechanisms. The use of life cycle assessment revealed that several reactive bed filters are relatively energy-consuming due to the material manufacturing process. Characterization of phosphorus compounds in used reactive media provided evidence for calcium phosphate precipitation as the predominating P removal mechanism in alkaline filter materials. However, in soil treatment systems with noncalcareous soils, batch experiments and extractions suggested that aluminium compounds were important for P removal. According to mass balance calculations that compared accumulated P with the estimated P load in a soil treatment system, the long term P removal capacity was very low; only 6.4 % of the applied phosphorus had been removed during 16 years of operation.</p>

Fosforutsläpp kan leda till övergödning i vattenrecipienten om kritiska nivåer överskrids. De enskilda avloppen står idag för en relativt stor andel av Sveriges fosforutsläpp jämfört med de kommunala avloppen. Detta kan sättas i bakgrund av att staten, alltsedan 1970-talet, har lagt ner stora resurser på att reducera utsläppen hos kommunala avlopp, medan en motsvarande satsning på de enskilda avloppen uteblivit. Merparten av dagens enskilda avlopp utgörs av markbaserade reningssystem, det vill säga avloppslösningar där avloppsvattnet får passera genom markmaterial för att renas. Fosforreningen i markbaserade reningssystem har forskningshistoriskt sett inte varit en prioriterad fråga, vilket har lett till kunskapsbrist inom området. Naturvårdsverket har i och med miljöbalkens införande skärpt sina krav på fosforreningen hos enskilda avlopp. Trots att Naturvårdsverket nu har satt upp exakta kravgränser på fosforutsläpp, exempelvis att en normal skyddsnivå motsvarar en rening på 70 %, saknas dock råd för hur gemene man, inför byggandet av ett markbaserat reningssystem, ska kunna uppskatta fosforreningsförmågan i markmaterialet. Emellertid rekommenderar Naturvårdsverket att man inför byggandet ska uppskatta markmaterialets hydrauliska konduktivitet, en uppskattning som ofta utförs genom att mäta markmaterialets kornstorleksfördelning. Det finns även forskning som tyder på att kornstorleksfördelningen skulle kunna ha samband med fosforreningen. Syftet med detta examensarbete var framförallt att utvärdera om kornstorleks- fördelningen kan användas som ett mått på fosforreningen. Som ett sidospår till detta utvärderades även om den hydrauliska konduktiviteten kan bestämmas ur kornstorleks- fördelningen. Empiriska formler har vid tidigare studier etablerats för detta, men frågan är om dylika empiriska formler gäller för de undersökta markbaserade reningssystemen. Undersökningens resultat tydde på ett samband mellan kornstorlek och hydraulisk konduktivitet medan det inte kunde påvisas något samband mellan kornstorleks- fördelning och fosforrening. Det är därmed önskvärt att finna ett annat mått på fosforreningen i markbaserade reningssystem. Kan vi inte uppskatta hur mycket de markbaserade reningssystemen bidrar till övergödningen kan vi inte heller väga denna negativa miljöpåverkan mot de markbaserade systemens fördelar, såsom relativt låg energiförbrukning och relativt små utsläpp av växthusgaser. / Release of phosphorus is one of the contributing factors of eutrophication in aquatic recipients. In Sweden, on-site wastewater treatment represents large emissions of phosphorous per person in comparison to municipal wastewater treatment. This can be put in the context of the Government’s large investments for reducing phosphorus in municipal wastewater treatment plants, while no corresponding investment has yet been made in the field of on-site wastewater treatment. In Sweden, the on-site wastewater treatment systems of today mainly consist of WSAS (Wastewater soil absorption systems), i.e. systems where wastewater is cleaned by passing though soil material. During the history of research about WSAS issues of phosphorous removal have not been first priority, which has led to a considerable lack of knowledge regarding these issues. Since 1999, when the Swedish regulations of environmental law were established, the SEPA (Swedish Environmental Protection Agency) has tightened up their demands of phosphorus removal in on-site wastewater treatment systems, for instance by stating that a phosphorous removal of 70% shall be reached to accomplish a “normal level of protection”. Despite this fact, there are no recommendations of how to estimate phosphorus removal in the soil material. However, the SEPA recommends estimation of hydraulic conductivity before building WSAS. This estimation is usually made by measuring grain size distribution in the soil material. One interesting fact is that there, according to some scientists, might be a connection between grain size distribution and phosphorus removal. The aim of this thesis was mainly to investigate if measure of grain size distribution can be used in order to estimate phosphorus removal. The aim was also to evaluate if grain size distribution is a good indicator of hydraulic conductivity in the investigated WSAS. It was concluded that grain size distribution is a suitable indicator of hydraulic conductivity, but not a suitable indicator of phosphorus removal. Thus, we still need to find an easy way to estimate the phosphorus removal in WSAS. If we do not know the extent of impact on eutrophication by WSAS, we cannot weigh this negative impact against the positive aspects of relatively low energy usage and low release of green house gases compared to other on-site wastewater treatment systems.

The prediction of moisture and contaminant transport through unsaturated soil to ground water is becoming increasingly important in the fields of hydrology, agriculture and environmental engineering. Computer aided simulation techniques enables one to conduct a series of systematic numerical experiments to analyze flow phenomenon in subsurface hydrology under various physical and chemical processes. The flow movement depends upon medium characteristics, initial and boundary conditions, which reflect, physical processes occurring below the ground. To understand the effects of physical process an efficient and accurate model is needed. Thus the model developed must be able to handle varied initial and boundary conditions. In this regard, infiltration into a very dry soil becomes a very important problem of study. Most of the earlier numerical models developed are concentrated on the development of an efficient algorithm or the modelling of a particular process which govern the flow in unsaturated or saturated-unsaturated homogeneous medium. Not much work has been done on the analysis of variably saturated flow in layered soil medium. Models to simulate unsaturated flow through dry soils, especially through layered soils with varied boundary conditions are very limited. Further, not much studies have been reported in the literature on the prediction of seepage face development and the phreatic surface movement in variably saturated media with layering. These aspects are very important in determining the flow field and the discharge from the domain. A detailed literature review covering above aspects has been made and is reported in this thesis. In the present work, two dimensional numerical models to predict the movement of wetting front in unsaturated domain and the movement of the phreatic surface in homogeneous and layered porous media under various initial and boundary conditions are developed based on finite difference and finite volume techniques. These models can handle flow in both rectangular flow domains and radial flow domains. The initial condition settings include the handling of very dry soil medium without any transformation of the governing equation, handling of infiltration and constant head conditions at the boundaries under steady state as well as transient scenarios. The models are also able to handle various soil moisture characteristics which depict the nonlinear behaviour between hydraulic conductivity, moisture content and pressure head in a soil media. A mixed form of the governing partial differential equation is used in the present study as it leads to better mass conservation. The finite difference model uses a central difference approximation for the space derivatives and an Eulerian backward difference approximation for the time derivative. The fully implicit formulation is solved iteratively using Strongly Implicit Procedure after making Picard approximation for the highly nonlinear coefficients. The process of infiltration into an initially dry soil leads to the development of a steep wetting front. As the finite volume technique is naturally an upwind method, it is expected to play a positive role in modelling such processes accurately. Hence, a finite volume model is also developed by approximating the convective part using a MUSCL approach and a fully implicit central difference method for the diffusive part of the governing equation. The models developed are validated using both experimental data and numerical solutions for problems reported in the literature. The validation problems cover a wide range of physical scenarios such as: infiltration into a very dry soil, infiltration into a dry soil column with gravity drainage, development of water table mound, steady state drainage in a sand filled wedge shaped tank with seepage face development and transient seepage face development in a rectangular domain. Five test problems are used for the validation of the models. The developed models perform very well for the test problems considered, indicating the models' capability in handling such situations. The results obtained by using the present models for simulating flow through highly unsaturated (very dry) soils show that the models perform very well when compared with models which use transformation techniques to handle such problems. The performance of the present models in comparison with the experimental data and numerical models available in the literature show the suitability of the present models in handling such situations. The present models are also used to analyse various types of unsaturated flow problems with varying initial and boundary conditions. The boundary conditions considered are no flow and /or free flow conditions along the left, right and bottom boundaries with infiltration condition along a part of the top boundary. For the various cases considered in the present study, infiltration rate varies from 5 cm/day to 50cm/day through an initially very dry soil of -15000 cm pressure head in homogeneous and layered soils. Different types of soil media considered vary from sandy loam, loam and clay with horizontal and vertical layering of these soils. A total number of 14 cases are analysed. The results are discussed in terms of pressure head variation in the flow domain along with moisture redistribution for all the cases under consideration. It is observed from these studies that the infiltration rate play an important role on the wetting front movement through layered soils depending on the type of layering and the boundary conditions considered. The soil properties of various layers affect the movement of wetting front by changing the direction of movement. Even though the wetting front movement is predominantly vertical, there is a tendency for the wetting front to move in the horizontal direction as it moves from a coarse soil to fine soil. It is also observed that the vertical layering of soils with different hydraulic conductivity helps in redirecting the flow towards the bottom boundary through the neighboring coarser layers. As finite volume method is more suitable for simulating sharp fronts, it is expected to perform better than finite difference method for simulating infiltration into very dry soils. So, a comparison is made between the performance of these two models by using the above test problems. It is observed from these studies that the performance of both the models are same except that the finite volume method takes much more CPU time than the finite difference model. Considering the type of problems tested, it is observed that for modelling unsaturated and saturated-unsaturated flows, finite difference method is better in comparison to finite volume method. It may be due to the predominant diffusive characteristics of the governing equation even while modelling flow through very dry soils. Proper estimation of the seepage height is an important aspect in finding the discharge through the porous medium. It is observed from the literature that the use of a saturated flow model in such situations can lead to an underestimation of the discharge through the porous medium. This effect is more important when dealing with small dimension problems. It is also observed that various parameters which govern the moisture movement through saturated-unsaturated regions affect proper estimation of the seepage face height and there by discharge. Various factors like effect of boundary conditions, type of soil layering, problem dimension and aspect ratio on seepage face development and the associated phreatic surface formation is studied in the present work. It is seen from the present study that the seepage face development is more in rectangular flow domain than in radial flow domain for both homogeneous and layered soils. It is also seen that the seepage face development in rectangular problems are more sensitive than radial flow problems for various factors considered. The seepage height is also influenced by the tail water level. It is seen from the present study that as the tail water level increases the seepage face reduces with no seepage face development for some of the cases studied. This influence is relatively less for radial flow problems. As the length of the domain increases the seepage height decreases. It is seen that for different cases with same horizontal dimension, as the height of the domain increases the seepage face height also increases. This phenomenon is observed for both homogeneous and layered soil medium. The influence of the aspect ratio, which is the ratio of the length to height of the domain indicate that as the aspect ratio increases the seepage height decreases. The type of the soil layering is observed to have a very strong influence on the seepage face development. The study for understanding the effect of soil layering on the development of seepage face and phreatic surface suggest that as the coarseness of the material increases, the phreatic surface become flatter and its steepness increases with the fineness of the soil. The present model is also used for studying the transient phreatic surface movement and the seepage face development. This is studied for homogeneous and layered rectangular soil medium. The present study is used to understand the effect of specific storage on the phreatic surface movement and the seepage face development. The studies indicate that the influence of specific storage on the seepage face development is insignificant in homogeneous soils with only very little effect in the early time for longer domains. It is also observed that the influence of the specific storage is significant in the case of layered soils. This effect depends on the type of layering and the problem dimension and is observed to have influence for relatively longer period. This observation suggests the importance of specific storage on transient seepage face development. When the specific storage effect is considered the drainage of the soil become faster resulting in a faster decline of the phreatic surface with time. The influence of specific storage is also studied considering the problem dimension effect. It is seen that as the aspect ratio increases, the effect of specific storage on the phreatic surface development decreases. The studies with change in the upstream boundary condition from a constant head to a no flow condition indicate that the effect of specific storage has no significant influence on the phreatic surface development for both homogeneous and layered soils.

CoordenaÃÃo de AperfeiÃoamento de Pessoal de NÃvel Superior / FundaÃÃo Cearense de Apoio ao Desenvolvimento Cientifico e TecnolÃgico / A erosÃo em parcelas irrigadas por sulcos à um fenÃmeno de especial importÃncia agrÃcola e ambiental que, tradicionalmente, no Brasil em especial, nÃo tem sido reconhecida nem tampouco estudada com detalhes. Com o objetivo de avaliar quantidade de solo erodido nesses sistemas, uma pesquisa foi conduzida em duas unidades experimentais da Fazenda Experimental Vale do Curu no municÃpio de Pentecoste (CearÃ- Brasil). As unidades experimentais foram escolhidas de acordo com texturas do solo representativas dos PerÃmetros Irrigados por superfÃcie no Nordeste brasileiro: franco arenosa e franco argilosiltosa , cada uma com 0,15 ha. Os experimentos foram conduzidos no perÃodo de julho a outubro de 2010 e de julho a outubro de 2012. Em cada unidade, nove sulcos foram analisados, correspondendo a trÃs tratamentos de vazÃo com trÃs repetiÃÃes para cada vazÃo, com um sulco-bordadura entre os sulcos-teste e dois sulcos sem irrigar delimitando os tratamentos de vazÃo. As vazÃes foram diferenciadas em funÃÃo dos diÃmetros dos sifÃes, sendo utilizados os de 1â, 1 Ââ e 2â de diÃmetro. As irrigaÃÃes foram realizadas em intervalos de no mÃximo sete dias. Durante o perÃodo do experimento foram avaliados: a concentraÃÃo de sedimentos na Ãgua de irrigaÃÃo; o efeito da perda de solo nas seguintes medidas de desempenho do sistema: eficiÃncia de aplicaÃÃo e coeficiente de uniformidade de Christiansen e os efeitos da erosÃo sobre a fertilidade do solo e a geometria da secÃÃo transversal ao fluxo de Ãgua. A pesquisa revelou que o excesso de Ãgua proveniente da irrigaÃÃo por sulcos carreou consigo altas concentraÃÃes de sedimentos, valores estes que variaram de 0,1 g L-1 a 7,0 g L-1, gerando uma perda de solo de atà 20 t ha-1ano-1. A perda de solo aumentou com a vazÃo aplicada na cabeceira do sulco segundo um modelo potencial a uma taxa crescente, entretanto, outros fatores igualmente importantes â comprimento do sulco, declividade e textura do solo â devem ser acrescidos à anÃlise, visando melhorar as prÃticas de manejo para minorar os problemas de erosÃo, otimizando assim a eficiÃncia hidrÃulica das irrigaÃÃes e a conservaÃÃo do solo. / Erosion in furrow irrigated plots is a phenomenon of agricultural and environmental importance that, particularly in Brazil, has not received the attention that deserves. Experiments were conducted in two experimental units of an experimental farm in Pentecoste (Cearà - Brazil) aiming to quantify the amount of eroded soil provoked by furrow irrigation. The experimental units were chosen to represent typical soil textures in the irrigated perimeters areas of Northeastern Brazil: sandy loam (B -T2) and silty clay loam (C - S2). Each unit had 0.15 ha. The experiments were conducted from July to October 2010 and from July to October 2012. In each unit, nine furrows were monitored, corresponding to three inflow rate treatments with three replications for each inflow rate, leaving a non-monitored furrow between the test furrows and two furrows without irrigation delimiting the flow rates treatments. Flow rates were differentiated using different diameters of siphons (1 ", 1 Â" and 2 "). Irrigations were performed at intervals of no more than seven days. During the experiments, the following variables were analyzed: sediment concentration in the runoff water, variation of the furrow cross section geometry, the relation of soil loss and system performance measures (application efficiency and ChristiansenÂs uniformity coefficient) and the effects of erosion on soil fertility. The research revealed that the excess water from the furrow irrigation carried high sediment concentrations, ranging from 0.1 g L-1 to 7.0 g L-1,resulting in a soil loss of up to 20 Mg ha-1 yr-1. Soil loss increased with the inflow rate following a potential function. Other important factors - furrow length, slope and soil texture - must be also considered to improve management that minimizes erosion and optimizes the irrigation efficiency.

For many surface waters, phosphorus (P) leaching is a serious problem that should be minimized to prevent eutrophication. In Sweden there is a demand for physical and technical development of high-performance P removal techniques to reduce phosphorus leaching from on-site wastewater treatment systems to the Baltic Sea. However, although these systems are designed to reduce eutrophication there are also other environmental impacts to be considered when implementing them in on-site systems; energy use and global warming potential are two examples. This study has investigated several bed filter materials (reactive media and natural soils) for their total environmental impact (in commercial applications) as well as for the predominating chemical phosphorus removal mechanisms. The use of life cycle assessment revealed that several reactive bed filters are relatively energy-consuming due to the material manufacturing process. Characterization of phosphorus compounds in used reactive media provided evidence for calcium phosphate precipitation as the predominating P removal mechanism in alkaline filter materials. However, in soil treatment systems with noncalcareous soils, batch experiments and extractions suggested that aluminium compounds were important for P removal. According to mass balance calculations that compared accumulated P with the estimated P load in a soil treatment system, the long term P removal capacity was very low; only 6.4 % of the applied phosphorus had been removed during 16 years of operation. / <p>QC 20110413</p>

Modern pig farming production may over burden the environment with organic substances, exposure of bacterial pathogens and introduction of resistance gene. This may be caused by the pig’s droppings, lack of seepage management or accidental spillage of seepage which may impact on the environment and its physicochemical parameters. The objective of this study is to determine and assess the level of bacteriological pollution emanating from the pig farm and their impact on the physicochemical parameters of soil and water as well as to identify the presence of antibiotic resistance gene of these prevailing bacteria. Soil and water samples were collected monthly for a period of six months (March- August 2013). Samples were collected at pig enclosures, soil 20 m and 100 m away from pig enclosures, constructed wetland used for treating pig farm wastewater, soil 20m and 100 m away from constructed wetland. Procedure followed for analysing soil and water samples includes physicochemical analyses, viable cell counts of 10-1 to 10-8 dilutions, identification of bacteria using API 20E test kit, antibiotic susceptibility analyses, and identification of resistance gene using molecular procedures. The media that were used for viable cell counts were, Nutrient agar, MacConkey Agar, Xylose Lysine Deoxycholate agar (XLD agar), and Eosin Methylene Blue (EMB). Physicochemical parameters of water showed unacceptable high levels of analysed parameters for BOD (163 mg/L to 3350 mg/L), TDS (0.77 g/L to 6.48 mg/L), COD (210 mg/L to 9400 mg/L), NO3 (55 mg/L to 1680 mg/L), NO2 (37.5 mg/L to 2730 mg/L), and PO43− (50 mg/L to 1427 mg/L) were higher than the maximum permissible limits set by Department of Water Affairs and Forestry (DWAF). For soil samples TDS (0.01g/L to 0.88 g/L), COD (40 mg/L to 304 mg/L), NO3 (32.5 mg/L to 475 mg/L), and NO2 (7.35 mg/L to 255 mg/L) and PO43- (32.5 mg/L to 475 mg/L ) were observed to be higher than recommended limits set by Federal Ministry for the Environmental (FME). The viable cells in soil samples 30cm depth ranged from 0 cfu/mL to 2.44 x 1010cfu/mL, in soil 5cm depth ranged from 1.00 x 101 cfu/mL to 1.91 x 1010 cfu/mL, and in water samples viable cells ranged from 5.00 x 101 to 5.05 x 109. Pseudomonas luteola (Ps. luteola), Escherichia vulneris (E. vulneris), Salmonella choleraesuis spp arizonae, Escherichia coli 1(E. coli 1), Enterobacter cloacae, Pseudomonas flourescens/putida (Ps. flourescens/putida), Enterobacter aerogenes, Serratia ordoriferal, Pasteurella pneumotropica, Ochrobactrum antropi, Proteus vulgaris group, Proteus vulgaris, Salmonella spp, Aeromonas Hydrophila/caviae/sobria1, Proteus Mirabillis, Vibrio fluvials, Rahnella aquatillis, Pseudomonas aeruginosa (Ps. aeruginosa), Burkholderia Cepacia, Stenotrophomonas maltophilia (St. maltophilia), Shwenella putrefaciens, Klebsiela pneumonia, Cedecea davisa, Serratia liquefaciens, Serratia plymuthica, Enterobacter sakaziki, Citrobacter braakii, Enterobacter amnigenus 2, Yersinia pestis, Serratia ficaria, Enterobacter gergoriae, Enterobacter amnigenus 1, Serratia marcescens, Raoutella terrigena, Hafnia alvei 1, Providencia rettgeri, and Pantoa were isolated from soil and water samples from the pig farm. Isolates were highly resistant to Penicillin G, Sulphamethaxazole, Vancomycin, Tilmocozin, Oxytetracycline, Spectinomycin, Lincomycin, and Trimethoprim. The most resistance genes detected in most isolates were aa (6’)-le-aph (2”)-la, aph (2”)-lb, aph (3”)-llla, Van A, Van B, Otr A and Otr B. Pig farm seepage is causing bacterial pollution which is impacting negatively on the natural environment in the vicinity of pig farm by introducing bacterial pathogens that have an antibiotic resistance gene and is increasing the physicochemical parameters for soil and water in the natural environment at the pig farm. It is therefore recommended that pig farms should consider the need to implement appropriate regulatory agencies that may include the regular monitoring of the qualities of final effluents from waste water treatment facilities. In addition there is a need to limit soil pollution in order to safe guard the natural environment in the vicinity of pig farm from bacteriological pollution and introduction of antibiotic resistance gene. It is also recommended that more advanced technologies should be introduced that will assist pig farms to manages the seepage properly. / Environmental Sciences / M. Sc. (Environmental Sciences)

Modern pig farming production may over burden the environment with organic substances, exposure of bacterial pathogens and introduction of resistance gene. This may be caused by the pig’s droppings, lack of seepage management or accidental spillage of seepage which may impact on the environment and its physicochemical parameters. The objective of this study is to determine and assess the level of bacteriological pollution emanating from the pig farm and their impact on the physicochemical parameters of soil and water as well as to identify the presence of antibiotic resistance gene of these prevailing bacteria. Soil and water samples were collected monthly for a period of six months (March- August 2013). Samples were collected at pig enclosures, soil 20 m and 100 m away from pig enclosures, constructed wetland used for treating pig farm wastewater, soil 20m and 100 m away from constructed wetland. Procedure followed for analysing soil and water samples includes physicochemical analyses, viable cell counts of 10-1 to 10-8 dilutions, identification of bacteria using API 20E test kit, antibiotic susceptibility analyses, and identification of resistance gene using molecular procedures. The media that were used for viable cell counts were, Nutrient agar, MacConkey Agar, Xylose Lysine Deoxycholate agar (XLD agar), and Eosin Methylene Blue (EMB). Physicochemical parameters of water showed unacceptable high levels of analysed parameters for BOD (163 mg/L to 3350 mg/L), TDS (0.77 g/L to 6.48 mg/L), COD (210 mg/L to 9400 mg/L), NO3 (55 mg/L to 1680 mg/L), NO2 (37.5 mg/L to 2730 mg/L), and PO43− (50 mg/L to 1427 mg/L) were higher than the maximum permissible limits set by Department of Water Affairs and Forestry (DWAF). For soil samples TDS (0.01g/L to 0.88 g/L), COD (40 mg/L to 304 mg/L), NO3 (32.5 mg/L to 475 mg/L), and NO2 (7.35 mg/L to 255 mg/L) and PO43- (32.5 mg/L to 475 mg/L ) were observed to be higher than recommended limits set by Federal Ministry for the Environmental (FME). The viable cells in soil samples 30cm depth ranged from 0 cfu/mL to 2.44 x 1010cfu/mL, in soil 5cm depth ranged from 1.00 x 101 cfu/mL to 1.91 x 1010 cfu/mL, and in water samples viable cells ranged from 5.00 x 101 to 5.05 x 109. Pseudomonas luteola (Ps. luteola), Escherichia vulneris (E. vulneris), Salmonella choleraesuis spp arizonae, Escherichia coli 1(E. coli 1), Enterobacter cloacae, Pseudomonas flourescens/putida (Ps. flourescens/putida), Enterobacter aerogenes, Serratia ordoriferal, Pasteurella pneumotropica, Ochrobactrum antropi, Proteus vulgaris group, Proteus vulgaris, Salmonella spp, Aeromonas Hydrophila/caviae/sobria1, Proteus Mirabillis, Vibrio fluvials, Rahnella aquatillis, Pseudomonas aeruginosa (Ps. aeruginosa), Burkholderia Cepacia, Stenotrophomonas maltophilia (St. maltophilia), Shwenella putrefaciens, Klebsiela pneumonia, Cedecea davisa, Serratia liquefaciens, Serratia plymuthica, Enterobacter sakaziki, Citrobacter braakii, Enterobacter amnigenus 2, Yersinia pestis, Serratia ficaria, Enterobacter gergoriae, Enterobacter amnigenus 1, Serratia marcescens, Raoutella terrigena, Hafnia alvei 1, Providencia rettgeri, and Pantoa were isolated from soil and water samples from the pig farm. Isolates were highly resistant to Penicillin G, Sulphamethaxazole, Vancomycin, Tilmocozin, Oxytetracycline, Spectinomycin, Lincomycin, and Trimethoprim. The most resistance genes detected in most isolates were aa (6’)-le-aph (2”)-la, aph (2”)-lb, aph (3”)-llla, Van A, Van B, Otr A and Otr B. Pig farm seepage is causing bacterial pollution which is impacting negatively on the natural environment in the vicinity of pig farm by introducing bacterial pathogens that have an antibiotic resistance gene and is increasing the physicochemical parameters for soil and water in the natural environment at the pig farm. It is therefore recommended that pig farms should consider the need to implement appropriate regulatory agencies that may include the regular monitoring of the qualities of final effluents from waste water treatment facilities. In addition there is a need to limit soil pollution in order to safe guard the natural environment in the vicinity of pig farm from bacteriological pollution and introduction of antibiotic resistance gene. It is also recommended that more advanced technologies should be introduced that will assist pig farms to manages the seepage properly. / Environmental Sciences / M. Sc. (Environmental Sciences)

Several decades of intensive dryland-farming in the Gadarif Region, located in the Eastern part of Sudan, has led to rapid landuse/landcover (LULC) changes mainly due to agricultural expansion, government policies and environmental calamities such as drought. The study area represents part of the African Sahel. The fundamental goal of the thesis was to assess land degradation and the impact of agriculture expansion on land cover, soil and crops production. To analyse and to monitor the LULC changes, multi-temporal Landsat data of the years 1979, 1989 and 1999 and ASTER data of the year 2009 covering an area of approximately 1200 km² were used. For this a post-classification comparison technique was applied to detect LULC changes from satellite images. Six LULC classes were identified during the classification scheme, namely cultivated land, fallow land, woodland, bare land, settlement and water. For the four dates of satellite images the overall classification accuracy ranged from 86 % to 92 %. During the three decades of the study period an extensive change of LULC patterns occurred. The cultivated areas increased significantly, covering 81 % of the previous woodland in the period 1979 – 2009. Fallow land only increased during the period 1989 – 1999. Over the three decades, urban expansion continuously increased covering an area of 23, 21 and 27 km² for the periods 1979 – 1989, 1989 – 1999 and 1999 – 2009 respectively. The detailed LULC map of the study area was obtained by using a dual polarisation (HH and HV) TerraSAR-X data of the year 2009. The different LULCs of the study area were analysed by employing an object-oriented classification approach. For that purpose, multi-resolution segmentation of the Definiens Software was used for creating the image objects. Using the feature-space optimisation tool the attributes of the TerraSAR-X images were optimised in order to obtain the best separability among classes for the LULC mapping. In addition to the classes that have been obtained by the optical data, the following LULCs resulted from SAR data: harvested land, rock, settlement 1 (local-roof buildings) and settlement 2 (concrete roof buildings). The backscattering coefficients for some classes were different along HH and HV polarisation. The best separation distance of the tested spectral, shape and textural features showed different variations among the discriminated LULC classes. An overall accuracy of 84 % with a kappa value of 0.82 was resulted from the classification scheme. Accuracy differences among the classes were kept minimal. For more than six decades in the Gadarif Region mechanised dryland farming is practised. As a result, due to continuous conventional tillage, extensive woodcutting and over-grazing, serious soil degradation occurred. To discuss the impact of LULC changes on the selected soil properties, three main LULC types were chosen to be investigated, namely: cultivated land, fallow land and woodland. In addition to the reference soil profiles, soil samples were also collected at two depths from ten sample plots for each of the LULC type. For these soil samples, various soil properties such as texture, bulk density (BD), organic matter (OM), soil pH, electrical conductivity (EC), sodium adsorption ratio (SoAR), phosphorous (P) and potassium (K) were analysed. Laboratory tests proved that soil properties were significantly affected by LULC changes. Within the different LULC types, clay content in the surface layers (0 – 5 and 5 – 15 cm) varied from 59 % to 65 %, whereas silt fractions ranged from 27 % to 37 %. Soil BD, OM and P were significantly different (p < 0·05) across the three LULC types. Soil pH was significantly different between cultivated land and woodland on one side and between fallow land and woodland on the other side. EC and SoAR values of fallow land were found to be significantly different (p < 0·05) from woodland. The dryland vertisol of the Gadarif Region in Sudan produced more than one-third of the national production of sorghum – the main food stuff in the country. Soil compaction has been recognised as one of the major problems in crop production worldwide. Soil strength and infiltration rate are important variables for understanding and predicting the soil processes. The effects of three different landuse systems (cultivated land, fallow land and woodland) on soil compaction and infiltration rate were investigated at two sites of the study area. Site 1 represents the older one of the two. The soil penetration resistance (SPR) was measured in three depths using a manually operated cone penetrometer. Infiltration rate was measured in the field using a double-ring infiltrometer. Following the cone-penetrometer sampling, soil samples were collected to determine the variables that affect SPR and infiltration rate vs. particle size, dry BD, volumetric moisture content (VMC) and organic carbon (OC) content. Field measurements and soil samples were collected for each landuse type. The measured infiltration rate data were inserted into the Kostiakov Model in order to predict the cumulative soil water infiltration. Soil compaction for the cultivated land was 65 % larger in comparison to woodland. Woodland areas showed an increase in the infiltration rate by 87 % and 74 % compared to cultivated and fallow land respectively. Both study sites showed an increase in the dry BD when SPR is increasing, while VMC decreases with increasing SPR. Also, low OC contents were observed to be associated with high SPR values. For Site 1 the average coefficient of determination (R²) for the infiltration data fit to the Kostiakov Model were 0.65, 0.73 and 0.84 for cultivated land, fallow land and woodland respectively. However, for Site 2 they were 0.63, 0.76 and 0.78. In the Gadarif Region agriculture is the main activity and practised in many forms with a variety of environmental effects and consequences. Continuous ploughing of the cultivated land coupled with inproper soil management has contributed to soil deterioration when the landuse changed from woodland to cultivated and fallow land. Therefore, the development of sustainable landuse practises in the dryland-farming of the study area need to be improved in order to reduce the amount of soil degradation in the future. / Mehrere Jahrzehnte intensiven Trockenfeldbaus in der Region von Gadarif, welche sich im östlichen Teil des Sudans befindet, führten hauptsächlich aufgrund von landwirtschaftlicher Expansion, politischen Beschlüssen der Regierung und Naturkatastrophen wie Trockenheit zu einer raschen Veränderung der Landnutzung und Landbedeckung. Das wesentliche Ziel dieser Dissertation war es, die Degradation des Landes, sowie die Auswirkungen von landwirtschaftlicher Expansion auf die Landbedeckung, den Boden und den Pflanzenbau im Untersuchungsgebiet, welches Teile der afrikanischen Sahelzone beinhaltet, abzuschätzen. Zur Analyse und Beobachtung der Veränderungen der Landnutzung und Landbedeckung wurden multi-temporale Landsat-Daten der Jahre 1979, 1989 und 1999 sowie ASTER-Daten aus dem Jahr 2009 genutzt, welche eine Fläche von schätzungsweise 1200 km² abdecken. Um Veränderungen von Landnutzung und Landbedeckung aus Satellitenbilddaten zu bestimmen, wurde ein auf Post-Klassifikation basierendes Vergleichsverfahren angewandt. Sechs Landnutzungs- und Landbedeckungsklassen, welche die Namen bewirtschaftetes Land, brach liegendes Land, Waldgebiet, Ödland, besiedeltes Land und Wasserfläche tragen, wurden während des Klassifikationsprozesses bestimmt. Für die vier Aufnahmezeitpunkte der Satellitendaten lag die allgemeine Klassifikationsgenauigkeit zwischen 86 % und 92 %. Während des dreißigjährigen Untersuchungszeitraums fand eine beträchtliche Veränderung der Landnutzungs- und Landbedeckungsstruktur statt. Bewirtschaftete Flächen nahmen in ihrem Anteil signifikant zu und bedeckten innerhalb des Zeitraums von 1979 bis 2009 81 % der früheren Waldgebiete. Der Anteil von brach liegendem Land nahm lediglich während des Zeitraums von 1989 bis 1999 zu. Besiedelte Gebiete breiteten sich über die drei Jahrzehnte kontinuierlich aus und wuchsen innerhalb des Zeitraums von 1979 bis 1989 um eine Fläche von 23 km², sowie um 21 km² zwischen 1989 und 1999 und um 27 km² in dem Zeitabschnitt 1999 – 2009. Eine detaillierte Karte zur Landnutzung und Landbedeckung des Untersuchungsgebiets wurde mittels der Nutzung dual polarisierter (HH und HV) TerraSAR-X Daten aus dem Jahr 2009 erzeugt. Die verschiedenen Landnutzungen und Landbedeckungen im Beobachtungsgelände wurden durch die Anwendung eines objektorientierten Klassifikationsansatzes analysiert. Um Bildobjekte zu erzeugen, wurde für diesen Zweck die auf einer mehrfachen Auflösung basierende Segmentierung der Software Definiens genutzt. Das Werkzeug Feature Space Optimisation wurde für die Optimierung der Attribute der TerraSAR-X Bilder angewandt, damit eine ideale Unterscheidungsfähigkeit entlang der Klassen für die Kartierung der Landnutzungen und Landbedeckungen erreicht werden kann. Zusätzlich zu jenen Klassen, welche mittels optischer Daten abgeleitet wurden, ergaben sich aus SAR-Daten noch die nachfolgenden Landnutzungen und Landbedeckungen: Abgeerntetes Land, Fels, Besiedlung 1 (Gebäude mit landestypischer Bedachung) und Besiedlung 2 (Gebäude mit Betondach). Die Koeffizienten der Rückstreuung entlang der Polarisationen HH und HV waren für einige Klassen unterschiedlich. Der günstigste Trennungsabstand der getesteten spektralen, formgebenden und texturalen Features ergab verschiedene Abweichungen zwischen den bestimmten Klassen der Landnutzung und Landbedeckung. Die Klassifikationsmaßnahmen ergaben eine Gesamtgenauigkeit von 84 % mit einem Kappa-Wert von 0.82. Genauigkeitsunterschiede entlang der Klassen wurden minimal gehalten. Seit über sechs Jahrzehnten wird in der Region Gadarif maschinenbetriebener Trockenfeldbau ausgeübt. In Folge dessen fand eine beträchtliche Abholzung und Überweidung sowie eine schwerwiegende Bodendegradation aufgrund des stetigen konventionellen Feldbaus statt. Um die Auswirkungen der Veränderung von Landnutzung und Landbedeckung auf die ausgewählten Bodenbeschaffenheiten auszuwerten, wurden drei Haupttypen der Landnutzung und Landbedeckung für die weitere Untersuchung ausgewählt: Bewirtschaftetes Land, brach liegendes Land, und Waldgebiet. Zusätzlich zu den Referenzbodenprofilen wurden außerdem für jeden Landnutzungs- und Landbedeckungstyp auf je zehn Probeflächen Bodenproben in zwei Tiefen entnommen. Bei diesen Bodenproben wurden zahlreiche Bodeneigenschaften analysiert, wie etwa Textur, Bodendichte (BD), organischer Materialgehalt (OM), pH-Wert des Bodens, elektrische Leitfähigkeit (EC), Adsorptionsgeschwindigkeit von Natrium (SoAR), Phosphorgehalt (P) sowie Kaliumgehalt (K). Labortests ergaben, dass die Bodeneigenschaften signifikant durch die Veränderungen der Landnutzung und Landbedeckung beeinflusst werden. Innerhalb der verschiedenen Landnutzungs- und Landbedeckungstypen variierte der Tongehalt in den Deckschichten (0 – 5 cm und 5 – 15 cm) zwischen 59 % und 65 %, wohin gegen sich die Lehmanteile von 27 % bis 37 % bewegten. Bodendichte, organischer Materialgehalt und Phosphorgehalt zeigten signifikant unterschiedliche Werte bei den drei Typen der Landnutzung und Landbedeckung (p < 0.05). Der pH-Wert des Bodens war signifikant verschieden zwischen bewirtschaftetem Land und Waldgebiet zum einen, und zwischen brach liegendem Land und Waldgebiet zum anderen. Die Werte der elektrischen Leitfähigkeit und der Adsorptionsgeschwindigkeit von Natrium bei brach liegendem Land erwiesen sich als maßgeblich verschieden zu jenen von Waldgebieten (p < 0.05). Auf dem Trockenland-Vertisolboden der Region Gadarif im Sudan wurde mehr als ein Drittel der nationalen Hirseproduktion erwirtschaftet – dem Haupternährungserzeugnis des Landes. Bodenverdichtung erwies sich als eines der weltweiten Hauptprobleme für den Pflanzenbau. Bodenfestigkeit und Versickerungsrate sind wichtige Variabeln, um Bodenprozesse verstehen und vorhersagen zu können. Die Auswirkungen der drei verschiedenen Landnutzungssysteme (bewirtschaftetes Land, brach liegendes Land und Waldgebiet) auf die Bodenverdichtung und Versickerungsrate wurden an zwei Standorten im Beobachtungsgebiet untersucht. Standort 1 ist der ältere der beiden. Der Widerstand der Bodenpenetration (SPR) wurde in drei Tiefen durch eine manuell angewandte Rammsonde gemessen. Mittels der Nutzung eines Doppelring-Infiltrometers ist die Versickerungsrate im Feld gemessen worden. Im Anschluss an die Probenentnahme mittels Rammsonden wurden Bodenproben gesammelt, um jene Variabeln bestimmen zu können, welche den Widerstand der Bodenpenetration sowie der Versickerungsrate im Vergleich zur Partikelgröße, zur trockenen Bodendichte, zum volumetrischen Feuchtigkeitsgehalt (VMC) und zum organischen Karbongehalt (OC) beeinflussen. Für jeden Landnutzungstypen wurden Feldmessungen durchgeführt und Bodenproben entnommen. Die gemessenen Daten der Versickerungsrate wurden in das Kostiakov-Modell eingespeist, um die gesamte Bodenwasserversickerung vorhersagen zu können. Die Bodenverdichtung bei bewirtschaftetem Land war 65 % stärker als bei Waldgebiet. Für Waldgebietsflächen wurde eine Zunahme der Versickerungsrate um 87 % verglichen mit bewirtschaftetem Land und um 74 % im Vergleich zu brach liegendem Land aufgezeigt. Beide Untersuchungsstandorte zeigten eine Zunahme in der trockenen Bodendichte für den Fall, dass der Widerstand der Bodenpenetration zunimmt, während der volumetrische Feuchtigkeitsgehalt mit zunehmendem Bodenpenetrationswiderstand abnimmt. Ebenso wurde beobachtet, dass ein geringer organischer Karbongehalt in Verbindung zu hohen Widerstandswerten der Bodenpenetration steht. Bei Standort 1 passte der durchschnittliche Bestimmungskoeffizient (R²) der Versickerungsrate zum Kostiakov-Modell mit den Werten 0.65 für bewirtschaftetes Land, 0.73 für brach liegendes Land und 0.84 für Waldgebiet. Für Standort 2 indessen ergaben die Werte 0.63, 0.76 und 0.78. Landwirtschaft, die in vielen Formen ausgeübt wird, ist die Haupttätigkeit in der Region Gadarif, und geht mit verschiedenartigsten Umweltauswirkungen und Konsequenzen einher. Kontinuierliche Feldbestellung des bewirtschafteten Landes, verbunden mit ungeeigneter Bodenbewirtschaftung, hat sich seit jenem Zeitpunkt, als sich die Landnutzung von Waldgebiet zu bewirtschaftetem und brach liegendem Land änderte, zu Bodenschädigung geführt. Daher muss die Entwicklung nachhaltiger Landnutzungspraktiken beim Trockenfeldbau im Untersuchungsgebiet verbessert werden, damit in Zukunft der Umfang der Bodendegradation verringert werden kann.

Les routes sont normalement pourvues de systèmes de drainage dimensionnés et exécutés conformément aux règles de l’art en vue d’évacuer, le plus rapidement possible de l’emprise de la route, les eaux des précipitations à considérer sur la base des périodes de retour prises en compte. Toutefois, il subsiste souvent des écoulements d’eau indésirables au niveau des talus et parfois dans les accotements et/ou les surfaces de roulement non imperméabilisées. Une succession de pluies entraîne des infiltrations d’eau qui varient notamment en fonction des conditions climatiques et suivant la texture et la structure du sol. De telles infiltrations ont pour conséquence la réduction des coefficients de sécurité des talus.<p><p>Bien qu’il existe plusieurs publications scientifiques traitant de pluies qui ont conduit à des glissements de terrain (Lim et al. 1996 ;Cho et al. 2001 ;Kim et al. 2004 ;Xue et al. 2007 ;Gavin et al. 2008), les incidences des infiltrations résultant de pluies successives sur le comportement des couches superficielles des sols non saturés ne sont généralement pas prises en compte. Les modèles permettant le calcul de la stabilité de talus des massifs de sols non saturés exigent beaucoup de paramètres parfois difficiles à évaluer et se rapportent ordinairement aux cas d’instabilité provoqués par une remontée du niveau piézométrique des nappes phréatiques.<p><p>Sur la base des essais réalisés en laboratoire, une méthodologie adaptée permettant d’évaluer la variation spatio-temporelle de la teneur en eau du sol en fonction d’une suite de pluies a été développée. Cette méthodologie facilite la prise en compte des effets cumulés des taux d’infiltration associés aux évènements pluvieux et permet d’en déduire le profil de succion ainsi que celui de la cohésion apparente à utiliser en vue de calculer, pour une inclinaison β du talus, l’intervalle de variation du coefficient de sécurité Fs. La méthodologie développée présente un intérêt particulier dans le cas de budgets et infrastructures limités.<p><p>ABSTRACT<p><p>Roads are normally equipped with drainage systems sized and implemented in accordance with the rules of art to evacuate as quickly as possible to the right of way, water precipitation to be considered on the basis of return periods taken into account. However, there are often water flows at the slope side and sometimes in the shoulders and / or running surfaces that are not waterproof. A succession of rain causes a certain amount of water infiltration, which varies according to climatic conditions and depending on the soil texture and structure. Such infiltrations have resulted in reduced safety factor of slopes.<p><p>Although there are several scientific publications on rainfall leading to landslides (Lim et al. 1996; Cho et al. 2001, Kim et al. 2004; Xue et al. 2007; Gavin and al. 2008), impacts resulting from infiltration of successive rains on the behavior of surface layers of unsaturated soils are usually not taken into account. Models for calculating the slope stability of unsaturated soils require many parameters that can be, in certain circumstances, difficult to assess and refer generally to cases of instability caused by a rise in groundwater level.<p><p>Based on laboratory tests, a suitable methodology for assessing the spatial and temporal variation of soil water content induced by a set of rains has been developed. This methodology facilitates the inclusion of the cumulative effects of the infiltration rates associated with rain events and infers from them the profile of suction and that of the apparent cohesion to be used to calculate, for a slope angle β, the range of variation of the safety factor Fs. This methodology is particularly relevant in the case of limited budgets and infrastructures. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished