<p>This thesis uses results of systematic numerical models to argue that externally inactive salt structures, which are potential targets for radioactive waste disposal, might be internally active due to the presence of dense layers or blocks within a salt layer.</p><p>The three papers that support this thesis use the Gorleben salt diapir (NW Germany), which was targeted as a future final repository for high-grade radioactive waste, as a general guideline.</p><p>The first two papers present systematic studies of the parameters that control the development of a salt diapir and how it entrains a dense anhydrite layer. Results from these numerical models show that the entrainment of a dense anhydrite layer within a salt diapir depends on four parameters: sedimentation rate, viscosity of salt, perturbation width and the stratigraphic location of the dense layer. The combined effect of these four parameters, which has a direct impact on the rate of salt supply (volume/area of the salt that is supplied to the diapir with time), shape a diapir and the mode of entrainment. Salt diapirs down-built with sedimentary units of high viscosity can potentially grow with an embedded anhydrite layer and deplete their source layer (salt supply ceases). However, when salt supply decreases dramatically or ceases entirely, the entrained anhydrite layer/segments start to sink within the diapir. In inactive diapirs, sinking of the entrained anhydrite layer is inevitable and strongly depends on the rheology of the salt, which is in direct contact with the anhydrite layer. During the post-depositional stage, if the effective viscosity of salt falls below the threshold value of around 10¹⁸-10¹⁹ Pa s, the mobility of anhydrite blocks might influence any repository within the diapir. However, the internal deformation of the salt diapir by the descending blocks decreases with increase in effective viscosity of salt.</p><p>The results presented in this thesis suggest that it is highly likely that salt structures where dense and viscous layer/blocks are present undergo an internal deformation processes when these dense blocks start sinking within the diapir. Depending on size and orientation of these blocks, deformation pattern is significantly different within the diapir. Furthermore, model results applied to the Gorleben diapir show that the rate of descent of the entrained anhydrite blocks differs on different sides of the diapir. This suggests that if the anhydrite blocks descent within the Gorleben diapir, they initiate an asymmetric internal flow within it.</p>

Water managers must balance the need for clean and safe drinking water with ever-increasing amounts of waste-water. A technique for treating and storing surface water called “managed aquifer recharge” (MAR) is frequently used to help maintain this balance. When MAR is used to produce drinking water, water managers must ensure that disease-causing microbial contaminants are removed from the water prior to its distribution. This thesis examined the processes responsible for removing a specific class of microbial contaminants called “enteric viruses” during MAR. Viruses are naturally removed in groundwater through adsorption and inactivation mechanisms. This thesis investigated how these virus removal mechanisms were affected by ionic strength (IS), dissolved organic carbon (DOC), and the age of the sand used in a MAR infiltration basin. This was done using batch and flow-through column experiments designed to mimic conditions characteristic of a basin infiltration MAR scheme in Uppsala, Sweden. Bacteriophage MS2 was used as a proxy for enteric viruses. All of the experiments were conducted at 4°C. Experimental data were modeled to describe the fate and transport of viruses in the infiltrated groundwater. Conventional least-squares optimization and generalized likelihood uncertainty estimation (GLUE) were compared as model fitting-approaches in order to determine how data uncertainty affects parameter estimates and model predictions. Results showed that the sand used in the infiltration basins accumulates adsorbed organic matter as it is exposed to infiltrating surface waters. This reduced the amount of MS2 that was removed due to adsorption and inactivation. Results from GLUE indicated that MS2 is more likely to inactivate in a time-dependent manner when in the presence of sand with high concentrations of organic matter. Both model fitting techniques indicated that virus attachment rates were significantly lower for sand with high organic carbon content. Neither methodology was capable of adequately capturing the kinetics of virus adsorption. Uncertainties in the experimental data had a large effect on the conclusions that could be drawn from fitted models. This study showed that the presence of natural organic matter reduces the value of the infiltration basin as a microbial barrier.

Les phénomènes d’instabilité générés par les séismes dans les pentes représentent un risque naturel majeur à l’échelle mondiale. Ce risque a, de plus, tendance à croître du fait d’une urbanisation croissante dans des zones à fort aléa, comme en témoignent de nombreux évènements catastrophique à travers le passé. Pour garantir la sécurité des personnes et des biens face à ce type d’aléa, il faut améliorer les prédictions des déplacements générés par les séismes dans les pentes.Les premiers travaux scientifiques visant à corréler les caractéristiques des instabilités dans les pentes aux paramètres sismiques sont de nature empirique et remontent aux années 1980. Ces méthodes ne permettent pas d’expliquer tous les mouvements de terrain observés à travers le monde, notamment lorsque les effets de site modifient sensiblement la distribution des mouvements dans les pentes et génèrent des déplacements plus grands ou plus petits que ceux prédits par les lois empiriques qui négligent les effets de site.La première partie de cette thèse présente une nouvelle base de données mondiale de 277 glissements de terrain d’origine sismique ou non. L’analyse statistique conduite sur les données de cette base a montré que la forme géométrique moyenne des glissements de terrain reste stable lorsque le volume des instabilités augmente.La deuxième partie de la thèse est dédiée à l’analyse de la stabilité du glissement de terrain de Diezma (Espagne) au moyen de trois méthodes : la Méthode de NEWMARK (1965) basée sur le principe de l’équilibre limite en 2D qui suppose des conditions elasto-parfaitement-plastiques ; l’analyse par différences finies dans le domaine temporel avec le code FLAC en 2D et en 3D qui suppose des conditions visco-elasto-parfaitement-plastiques ; et l’analyse par recombinaison modale dans le domaine fréquentiel avec le code CESAR-LCPC en 2D et en 3D qui suppose des conditions visco-élastiques.Situé dans une région sismiquement active, le glissement de terrain de Diezma est susceptible d’être affecté par des séismes. Une analyse des déplacements induits dans ce versant par onze signaux sismiques différents a été menée dans l’objectif d’établir des corrélations entre les paramètres macro-sismiques des scenarios sismiques appliqués et les déformations calculées par les trois méthodes.Les résultats montrent des différences marquées à la fois qualitatives et quantitatives. Les conclusions principales sont : (a) Les déplacements obtenus par les trois méthodes ont des ordres de grandeur différents et peuvent donc conduire à une surestimation ou à une sous-estimation des déplacements ; (b) Les deux méthodes numériques montrent que les déplacements sont fortement conditionnés par la géométrie des modèles ; seule la méthode par différences finies fait apparaître une réponse en termes de déplacements qui dépend du scénario sismique considéré; (c) Les résultats des simulations2D et 3D ne sont pas comparables ; des analyses complémentaires doivent encore être menées pour guider l’utilisateur dans le choix de la méthode la plus appropriée; (d) La Méthode de NEWMARK (1965), dont l’utilisation reste très répandue de nos jours, est tout à fait appropriée à l’étude de cas simples mais elle peut se révéler inexacte lorsque la structure géologique / topographique du versant conduit à un fort effet de site car ce dernier n’est pas pris en compte par cette méthode / Seismically induced slope deformation is a worldwide common phenomenon that poses an increasing and considerable threat to fast expanding urbanization, and a great number of catastrophic events throughout the past attest thereof. For this reason, displacement predictions allowing for proper slope surveillance became a major concern.Early attempts to relate slope failures to seismic parameters are of empirical nature and date back to the 1980s. Although having proven stable, these relations are frequently disturbed by site effects causing outliers in terms of smaller or greater displacements than expected.The first part of this thesis presents a newly build chronological database of 277 globally distributed seismically and non-seismically induced landslides. A comprehensive statistical analysis was conducted on the data of this database with the main result that – statistically seen – the average geometrical shapes of landslides differ only proportionally.The second part of the thesis is dedicated to a comparative slope stability analysis of the Diezma Landslide (Spain) by means of three methods: the limit-equilibrium based NEWMARK-Method (1965) in 2D under elasto-perfectly-plastic conditions; finite-difference analysis in the time domain with the code FLAC in 2D and 3D under visco-elasto-perfectly-plastic conditions; and modal recombination analysis in the frequency domain with the finite-element code CESAR-LCPC in 2D and 3D under visco-elastic conditions.Located in a seismically active region, the Diezma Landslide is likely to be affected by earthquakes and was therefore considered as a suitable model case. A broad analysis of expected displacements was conducted using eleven strong-motion seismic scenarios. Moreover, the thesis searches for potential relations between macro-seismic parameters of the applied earthquake scenarios and the predicted deformations obtained from the three methods.It appeared that – due to the functionality of the methods – results differ quantitatively as well as qualitatively, and so does their suitability. Major findings are: (a) Results from the three methods are of different orders of magnitude and, thus, can easily lead to over- or under-estimations of displacements; (b) Both numerical methods reveal a strong influence of the model-geometry on the predicted displacements, whereas a scenario-dependent slope behavior manifested itself only within the finite difference analysis; (c) The switch from 2D to 3D does not necessarily result in a similar performance in each dimension and results must be critically judged before further use; (d) The NEWMARK-Method (1965) has proven itself once more to be appropriate for first slope assessments but not for sophisticated evaluations of ground response to seismic shaking

A modelagem matemática do escoamento em duas dimensões é uma ferramenta importante para a pesquisa hidráulica e o desenvolvimento de projetos. Nesta dissertação é apresentado um estudo da hidrodinâmica bidimensional através de modelos matemáticos, sendo inicialmente introduzidas as equações básicas da mecânica dos fluídos que regem o movimento em duas direções. Em seguida são discutidas alguns modelos apresentados pela bibliografia, de onde podem ser obtidas valiosas informações sobre a técnica de modelação. Após a apresentação do método das diferenças finitas, um modelo matemático implícito é proposto para a aplicação em estudos bidimensionais onde é necessário o cálculo detalhado do campo de velocidades, evitando-se as restrições devido a estabilidade e acuracidade do método de diferenças. Como conclusão, são apresentados alguns exemplos de aplicação e comparações dos resultados com outros modelos. / Numerical modelling of two-dimensional flow is very important as a tool for hydraulic research and design. This work starts with a review of basic fluid mechanics equations of the wave motion in two-directions. Then some aspects of models proposed by several authors are discussed. The Finite Diference method is presented and na implicit model is proposed for aplication in detailed studies. The stability and acuraccy are investigated to recognize the range of model aplication. Some practical examples and comparisons are shown to test model capabilities.

O termo macro-rugosidade é empregado para caracterizar os elementos contínuos ou esparsos de formas diversas, que são aplicados artificialmente no fundo e paredes dos condutos hidráulicos, em oposição à rugosidade de grão, relacionada às asperezas das paredes destes condutos. Neste trabalho apresenta-se o estudo teórico-experimental do comportamento do escoamento livre com a presença de macro-rugosidade constituída de elementos prismáticos de seção transversal retangular, denominados travessões, que encontram aplicações práticas na estabilização de leitos e regularização de canais para navegação. O emprego da macro-rugosidade resulta na ocorrência de maiores perdas de energia quando as profundidades são pequenas e na redução progressiva deste efeito à medida que a lâmina d\'água se eleva. O estudo teórico compreende o tratamento dado pela literatura ao problema da resistência ao escoamento em canais, desenvolvido a partir da integração da tradicional lei logarítmica para o para o perfil bidimensional de velocidades. Também foi abordado o modelo de distribuição probabilística de velocidades de CHIU, baseado na maximização da entropia. Estudos específicos sobre macro-rugosidade de fundo, desenvolvidos a partir da década de 1960, foram pesquisados para orientação dos estudos experimentais desenvolvidos, notadamente os trabalhos de MORRIS & WIGGERT, SAYRE & ALBERTSON e KNIGHT & MACDONALD. Ensaios experimentais foram desenvolvidos em canal de laboratório para acaracterização das perdas de carga dos travessões de fundo, em diferentes condições de profundidade e espaçamento longitudinal. Nestes ensaios desenvolveu-se uma técnica para medição das perdas de energia e cálculo das tensões de atrito em escoamentos permanentes e gradualmente variados no espaço, situação esta usual nos canais empregados em engenharia. ) Foi também considerado o efeito de parede, resultante do emprego de um canal de laboratório de pequenas dimensões, de modo a poder-se generalizar os resultados para o caso bidimensional. Os dados experimentais permitiram o teste dos modelos para cálculo de perdas de carga de Manning e do Perfil Logarítmico, estabelecendo-se correlações entre as perdas e as características geométricas da macro-rugosidade. Como contribuição ao tema foi proposta uma adaptação do modelo logarítmico considerando a hipótese de MORRIS & WIGGERT da separação do escoamento em uma zona de interferência, denominada de zona de hiper-turbulência e outra, de turbulência normal. Os resultados obtidos indicaram a validade do modelo proposto. Ainda como contribuição ao tema, a lei de distribuição probabilística das velocidades foi aplicada, obtendo-se uma estimativa dos parâmetros para as diferentes configurações geométricas. / The designation large-scale roughness is generally employed to describe spot or continuous elements artificially disposed in channel bottom, different from wall roughness or grain roughness, associated to the wall material. This work presents a theoretical and experimental study of the free flow in channels with prismatic rectangular elements regularly disposed on the bottom, named traverse blocks, which have large application in river stabilization and navigation works. The use of the large-scale roughness results in a great head loss for small depths and a progressive reduction of these losses as the depth increases. The theoretical analysis initially focused the traditional approach considering the logarithmic velocity distribution law. It was also analyzed the probabilistic velocity distribution law by CHIU, based on the maximum entropy principle. Specific studies about large-scale roughness, developed since 1960, where considered to guide the model testes, especially those by MORRIS & WIGGERT, SAYRE & ALBERTSON and KNIGHT & MACDONALD. Experimental data where collected from a laboratory channel and a technique for computing bottom shear stress were developed allowing to consider permanent gradually varied flows, a regular situation in engineering applications. The wall interference in total channel shear stress where considered in order to make data comparable with those from real two-dimensional flow. Practical correlations using geometriccharacteristics of the traverse block roughness and flow parameters where established considering Manning\'s Formula and the Logarithmic Law model. In order to achieve a better modeling of the phenomenon, a complementary term to the Logarithmic law was proposed, considering the hyper-turbulent layer over the roughness, and fit to the data. Finally, parameter estimation for the use of the CHIU\'s law with large-scale roughness where developed as another contribution to the matter.

Over the past three decades, technological innovations with respect to cemented paste backfill (CPB) as a means of ground support has allowed for increased production within the mining industry, management mine waste costs, as well as the improvement of the overall health and safety of underground mining operations. Despite the extensive use of this relatively new ground support material, many fundamental factors affecting the design of safe and economical CPB structures are still not well understood.Recently, a significant amount of academic and industry research has been conducted to better understanding the distribution of stress with respect to primary-secondary extraction sequencing for stope-and-fill mining operations. While current, as well as past research, as provided a wealth of knowledge on the distribution of stress through the fill material itself, it lacks in providing an examination into the mechanism by which stress is able to redistribute itself through the backfill material as well as within the surrounding rockmass. The scope of this work is to optimize stope-and-fill extraction sequencing through the analysis of stress distributions as well as local and global stability of multiple narrow verticalfully-drained backfilled stopes. Scientific investigations into the behavior of the CPB material and surrounding rockmass will result in animproved understanding of how to better implement engineered paste-fill materials as a means of ground support for underground mining operations. Numerical simulations (FLAC3D and RocScience) were utilized in analyzing hypothetical (literature) as well as site-specific (field) case studies. While these simulations confirm generalized stress behaviors within the backfill material for single and adjacent stopes, stress redistributions within the surrounding rockmass as well as the rock-pillarindicate the development of tensile and compressive zones. From these results, one is able to better approximate ground and CPB instability with respect to site-specific conditions, geometries, and material properties. These simulations have been validated with respect to published analytical solutions, numerical simulations, and site-measurements for single (isolated) and adjacent narrow vertical fully-drained backfilled stopes.

The salinization of groundwater resources is a widespread problem in arid agricultural environments. In Mewat, the amount of solutes dissolved in the water has become too high to use for drinking or agriculture. The only fresh water recharge to this bowl-shaped region is through precipitation, which is focused at the foothills of the mountain. The freshest water is found closest to the mountains and the salinity of the groundwater increases as the distance from the mountains increase. The pumps that supply the region with fresh water are located in the shrinking freshwater zone. Locally-monitored wells show the movement of salinity in the region, as the saline water encroaches upon the freshwater. This study aims to answer the following questions: How long until the region runs out of fresh water? What would it take to have sustainable fresh water supplies? Is it even possible to have sustainable fresh water supplies in this environment? In order to answer these questions, we will quantify potential futures for an arid, groundwater-dependent location in rural India, using numerical groundwater modeling to quantify interactions between human water use, infrastructure, and climate. Outcomes of this modeling study will inform sustainable management of groundwater resources

The mountainous areas in the Intermountain West (IMW) of the North America are considered as the major water reservoir for the Western US. Summer evapotranspiration (ET) and soil moisture are key factors affecting the annual water yield in the montane region of the IMW. This research estimated ET of four common vegetation types (aspen, conifer, grass, and sage) and areal soil moisture in an advanced instrumentation site located at the T.W. Daniel Experimental Forest (TWDEF). Among instrumented forest research sites worldwide, TWDEF is one of a few with triplicate measures of meteorological parameters, radiation, and soil moisture within four common vegetation types in the IMW. This unique dataset enables study and understanding of the ecological and hydrological responses to climate change in Utah and the IMW region. In a second phase of this study, summer water uses from the four common vegetation types were simulated using a numerical simulation model, Hydrus-1D. The simulation was informed by soil moisture measurements at three depths (0.1 m, 0.25 m, and 0.5 m) and by ET measured from an eddy covariance tower. The results confirmed the value of numerical simulations as a viable alternate method to estimated ET where no direct ET measurements are available. It also provided comparison of water use by these vegetation species including both high and low water years. In the third phase of this study, a comparison was made between the intermediate-scale areal soil moisture measured by a Cosmic-ray neutron probe (CRNP) and the in situ TDT soil moisture network at the TWDEF site. Improved correlations were obtained, especially after shallow rainfall events, by including numerically simulated soil moisture above 0.1 m where no measurements were available. The original CRNP calibration exhibited a dry bias during spring/early summer, leading to the need for a site-specific enhanced calibration, which improved the accuracy of the CRNP soil moisture estimate at the TWDEF site.

Both primary and secondary liver cancers are common and the majority of patients are not eligible for surgical resection or a liver transplant, which are considered the only hope of cure. Mortality rates are high and there is a need for alternative treatment options. New forms of local treatment work best on small tumours; large ones, however, remain difficult to treat. Hyperthermia involves heating tumours to 40??-44?? C. The aim is to heat the entire tumour without damaging the surrounding normal tissue. Treating deep seated tumours is technically challenging. Ferromagnetic embolisation hyperthermia (FEH) is a novel method of treating liver tumours. Magnetic microspheres are infused into the hepatic artery and lodge primarily in the tumour periphery. An applied alternating-current magnetic field causes the microspheres to heat. Animal experiments have shown that this is a promising technique. There is a need for modelling of FEH prior to commencement of clinical trials. Analytical and numerical models of tumour heating during FEH treatment are presented here. The models help predict the temperature distributions that are likely to arise during treatment and give insight into the factors affecting tumour and liver heating. The models incorporate temperature-dependent thermal properties and blood perfusion rates of the tissues and a heterogeneous clustering of microspheres in the tumour periphery. Simulations show that the poorly perfused tumours heat preferentially while the liver is effectively cooled by blood flow from the portal vein. A peripheral distribution of heat sources produces a more even temperature field throughout the tumour, compared to a heat source that is centred within the tumour core. Large tumours reach higher temperatures and have higher heating rates, supporting experimental findings. Using temperature-dependent, rather than constant, values for thermal conductivities and blood perfusion rates results in higher temperatures within the tumour. The uneven clustering of microspheres in the tumour periphery leads to a more heterogeneous temperature distribution in the core, but it has less of an effect on the wellperfused liver. The results show that FEH has the potential to effectively treat liver tumours and the technique merits further investigation.

This thesis uses results of systematic numerical models to argue that externally inactive salt structures, which are potential targets for radioactive waste disposal, might be internally active due to the presence of dense layers or blocks within a salt layer. The three papers that support this thesis use the Gorleben salt diapir (NW Germany), which was targeted as a future final repository for high-grade radioactive waste, as a general guideline. The first two papers present systematic studies of the parameters that control the development of a salt diapir and how it entrains a dense anhydrite layer. Results from these numerical models show that the entrainment of a dense anhydrite layer within a salt diapir depends on four parameters: sedimentation rate, viscosity of salt, perturbation width and the stratigraphic location of the dense layer. The combined effect of these four parameters, which has a direct impact on the rate of salt supply (volume/area of the salt that is supplied to the diapir with time), shape a diapir and the mode of entrainment. Salt diapirs down-built with sedimentary units of high viscosity can potentially grow with an embedded anhydrite layer and deplete their source layer (salt supply ceases). However, when salt supply decreases dramatically or ceases entirely, the entrained anhydrite layer/segments start to sink within the diapir. In inactive diapirs, sinking of the entrained anhydrite layer is inevitable and strongly depends on the rheology of the salt, which is in direct contact with the anhydrite layer. During the post-depositional stage, if the effective viscosity of salt falls below the threshold value of around 1018-1019 Pa s, the mobility of anhydrite blocks might influence any repository within the diapir. However, the internal deformation of the salt diapir by the descending blocks decreases with increase in effective viscosity of salt. The results presented in this thesis suggest that it is highly likely that salt structures where dense and viscous layer/blocks are present undergo an internal deformation processes when these dense blocks start sinking within the diapir. Depending on size and orientation of these blocks, deformation pattern is significantly different within the diapir. Furthermore, model results applied to the Gorleben diapir show that the rate of descent of the entrained anhydrite blocks differs on different sides of the diapir. This suggests that if the anhydrite blocks descent within the Gorleben diapir, they initiate an asymmetric internal flow within it.