L’hydrodynamique et la probabilité de minéralisation (R²AI) autour des intrusions magmatiques ont été étudiées par modélisation numérique couplant transfert de chaleur et circulation de fluide. L’objectif principal de ce travail est de tester la nature du lien génétique entre l’intrusion et le processus de minéralisation. La première série de résultats s’appuie sur une comparaison avec des exemples naturels de gisements bien connus : i) L’hydrodynamique et la localisation des zones probables de minéralisation sont fortement dépendantes de la profondeur de mise en place du pluton. Au-dessus de 4.5km de profondeur de mise en place, le seuil de perméabilité de 10-16 m² est atteint et les cellules convectives peuvent créer des zones de décharge additionnelles où des minéralisations peuvent avoir lieu ; ii) Pour toutes les profondeurs d’emplacement, la zone en dessous du pluton n'est pas favorable à la précipitation minérale ; iii) Les apophyses focalisent les fluides convectifs et les zones de minéralisation autour d’elles ; iv) La phase de refroidissement n'est pas la phase majeure de convection. La zone advective principale et celle de haute favorabilité peuvent se produire avant et pendant la phase la plus chaude d’emplacement, avant que le magma ne cristallise complètement; v) Les détachements sont capables de fortement modifier et de re-localiser les flux convectifs déclenchés par une intrusion syn-tectonique; vi) Les conditions physiques favorables à la minéraliser sont produites pendant une durée courte autour de la phase la plus chaude de l'intrusion. Même si les arguments chimiques sont absents, la circulation de fluide (induite par la mise en place de magma) joue un rôle principal dans la genèse des gisements d'or associés aux intrusions. De plus, la formation de ce type de gisement est favorisée par l'occurrence d'une auréole thermique fracturée autour de l'intrusion. La seconde série de résultats concerne l’étude du cas naturel de la minéralisation W-Au de Tighza (Jebel Aouam) au Maroc Central. Une campagne d’acquisition de données gravimétriques, l’inversion données et l’utilisation de logiciel 3D, ont permis d’obtenir la géométrie 3D complexe du pluton de Tighza. Les résultats sont les suivants : i) la zone probable de la minéralisation apparaît au début de la mise en place du magma dans la zone perméable (veine W1) et s'étend pour remplir W1 et couvrir la région autour du pluton pendant la phase la plus chaude de mise en place; ii) lors du refroidissement, la zone probable est réduite et limitée à la zone perméable (W1) pendant 0,6 Ma; iii) L’application de la température de fermeture isotopique de la muscovite et de la biotite avec la distribution du R²AI montre que les âges de refroidissement entre la minéralisation au niveau de la veine W1 et l'intrusion ne sont pas séparés de plus de 0,10 Ma. Ceci est confirmé par la datation absolue de la minéralisation de Tighza et permet de discuter la fiabilité des âges obtenus pour la minéralisation dans la veine W1. / Coupled hydro-thermal numerical modeling has been used to simulate the hydrothermal fluid flow regime and the mineralization probability (R²AI) around plutons. The main objective behind this work is to test the nature of the genetic link between mineralization and intrusions. The first series of results comes from comparison with well-constrained mineral deposits: i) Fluid circulation and mineralization patterns are strongly dependent of the emplacement depth of the pluton. Deep seated plutons emplaced below 10 km do not induce an advective heat dissipation. For emplacement depth less than 4.5 km, the permeability threshold of 10-16 m2 is reached and second order convection cells may create additional discharge zones where mineralization are expected; ii) For all emplacement depths, the pluton floor zone is not favorable for mineral deposition; iii) The apexes strongly modified the fluid flow patterns by focusing convective fluids and mineralization zones around them; iv) The cooling phase is not the main phase of convection for large pluton often associated with long-lived magma emplacement. Major advective heat dissipation and mineral deposition zone may also occur sometime before and during the hottest phase of emplacement; v) Extensional detachments faults are able to delocalize and strongly modify classical fluid flow patterns induced by coeval intrusion; vi) Favorable physical conditions for mineral deposition are encountered around middle crust pluton, during a short time span bracketing the hottest phase of intrusion. We conclude that, even if chemical arguments are absent, fluid circulation induced by granite emplacement plays a key role in the genesis of granite-related Au deposits. Moreover, formation of this type of deposit is promoted and controlled by the occurrence of a fractured thermal aureole around the intrusion. The second series of results deals with the W-Au granite related Tighza deposits (Jebel Aouam, Morocco). Based on gravimetric data, inversion, and 3D modeling software, we were able to construct the most probable complex geometry of the Tighza pluton. The 3D geometries of the pluton and major fractures (W1 vein) were injected in the hydro-thermal modeling procedure. The results are: i) the probable zone of mineralization appears at the beginning of magma emplacement within the permeable zone (W1 vein) and extends to fill up W1 and covers the area around pluton at the hottest phase; ii) During the cooling phase, the story was reversed; the probable zone was reduced and restricted in the permeable zone (W1) during 0.6 Myr of cooling; iii) Application of isotopic closure temperature of muscovite and biotite coupled with R²AI distribution shows that the cooling ages between mineralization in W1 veins and the intrusion are not separated by more than 0.10 Myr. This is confirmed by the absolute dating obtained for Tighza Au mineralization and allows discussing the significance of older ages obtained for the mineralization in W1 veins.

Le comportement thermique et thermo-mécanique des pieux énergétiques est étudié par plusieurs approches : mesures au laboratoire sur des éprouvettes de sol, modélisation physique en modèle réduit, expérimentations sur pieu en vraie grandeur, et calculs numériques/analytiques. D’abord, la conductivité thermique d’un loess à l’état non saturé est mesurée en fonction de la teneur en eau et de la succion. Les résultats montrent une relation univoque entre la conductivité thermique et la teneur en eau pendant un cycle d’humidification/séchage alors qu’une boucle d’hystérésis est observée pour la relation entre la conductivité thermique et la succion. Deuxièmement, des essais thermiques sont réalisés sur un pieu énergétique expérimental en vraie grandeur pour étudier le transfert thermique à l’échelle réelle. Troisièmement, une solution analytique est proposée pour simuler la conduction thermique d’un pieu énergétique vers le sol environnant pendant un chauffage. Les tâches mentionnées ci-dessus concernant le comportant thermique sont ensuite complétées par des études sur le comportement thermo-mécanique des pieux énergétiques. D’un côté, des expérimentations sont réalisées sur un modèle réduit de pieu installé dans un sable sec ou dans une argile saturée. Trente cycles thermiques, représentant trente cycles annuels, sont appliqués au pieu sous différentes charges axiales en tête. Les résultats montrent un tassement irréversible avec les cycles thermiques ; ce tassement est plus important sous une charge axiale plus grande. De plus, le tassement est plus marqué pendant les premiers cycles thermiques et devient négligeable pour les cycles suivants. De l’autre côté, les travaux expérimentaux sur le modèle réduit de pieu sont complétés par les calculs numériques utilisant la méthode des éléments finis. Cette approche est d’abord validée avec les résultats obtenus sur le pieu modèle avant d’être utilisée pour prédire les résultats des expérimentations en vraie grandeur / The thermal and thermo-mechanical behavior of energy piles is investigated by various approaches: laboratory measurement on small soil samples, physical modeling on small-scale pile, experiments on real-scale pile, and analytical/numerical calculations. First, the thermal conductivity of unsaturated loess is measured simultaneously with moisture content and suction. The results show a unique relationship between thermal conductivity and moisture content during a wetting/drying cycle while a clear hysteresis loop can be observed on the relationship between thermal conductivity and suction. Second, thermal tests are performed on a full-scale experimental energy pile to observe heat transfer at the real scale. Third, an analytical solution is proposed to simulate conductive heat transfer from an energy pile to the surrounding soil during heating. The above-mentioned tasks related to the thermal behavior are then completed by studies on the thermo-mechanical behavior of energy piles. On one hand, experiments are performed on a small-scale pile installed either in dry sand or in saturated clay. Thirty thermal cycles, representing thirty annual cycles, are applied to the pile under various constant pile head loads. The results show irreversible pile head settlement with thermal cycles; the settlement is higher at higher pile head load. In addition, the irreversible thermal settlement is the most significant during the first cycles; it becomes negligible at high number of cycles. On the other hand, the experimental work with small-scale pile is completed with numerical calculations by using the finite element method. This approach is first validated with the results on small-scale pile prior to be used to predict the results of full-scale experiments

La réfrigération magnétique est une technologie émergente grâce à des avantages considérables par rapport aux technologies de réfrigération classiques. Cette technologie basée sur l’effet magnétocalorique offre d’importants avantages environnementaux car d’une part l’efficacité théorique des cycles utilisés est supérieure à celle des technologies classiques et d’autre part son fonctionnement ne nécessite pas une utilisation de gaz/vapeur `a fort effet de serre. En revanche des verrous scientifiques restent à lever, Le modèle multi-physique proposé dans cette thèse à pour but d’améliorer la précision de calcul. Il consiste à coupler un modèle 3D magnétostatique résolu par la méthode des éléments finis, un modèle magnétocalorique analytique et un modèle thermo-fluidique résolu par méthode des différences finies. Parallèlement, un banc d’essais a été conçu, optimisé et réalisé, ce banc permettra de faire des mesures fines des différents phénomènes qui interagissent dans la réfrigération magnétique. / The magnetic refrigeration technology is a promising alternative technology to the production of cold. The work of this thesis deals with studying and designing a magnetic refrigeration prototype. A multiphysic model is developed, this model taking into account several magnetic and magnetocaloric aspect that never dealt in the literature. It is used to investigate the influence of a range of parameters on the performance of the AMR. A new test bench of magnetic refrigeration is also designed, it is based on a particular electromagnet that was optimally realized. The magnetic performances are showed and concord with design prevision. Finally, suggestions for future works are provided based on the knowledge presented here.

Foi realizado um estudo com dados físicos em superfície provenientes de implementação do modelo numérico POM e dados de clorofila a do sensor MODIS, no Atlântico Sul e Tropical, tendo como finalidade analisar as variações espaciais de suas correlações. O modelo hidrodinâmico foi submetido a uma validação com dados do Projeto PIRATA e do sensor MODIS, obtendo-se resultados satisfatórios, principalmente em superfície. A clorofila a apresenta altos valores de correlação com variáveis físicas em superfície, porém esses valores são influenciados pela presença dos sinais anual e semianual nas séries temporais. Consequentemente, as análises foram desenvolvidas após a remoção desses sinais das séries. As correlações cruzadas entre as variáveis físicas e clorofila a tiveram valores altos em módulo apresentando, de maneira geral, padrão espacial característico de correlações negativas no interior do Giro Subtropical e positivas nas demais regiões para temperatura e salinidade em superfície e elevação do nível do mar. Esse padrão espacial é mais evidente para elevação e temperatura em superfície, apresentando defasagem no tempo de resposta biológica de quase zero no interior do Giro Subtropical e Equador, com predomínio de atrasos nessa resposta nas demais regiões. A covariância cruzada foi calculada para alguns pontos ao longo do meridiano de 20ºW, para elevação de superfície e temperatura, verificando-se concordância com os resultados obtidos pela correlação cruzada. A análise espectral dessa covariância salientou a ocorrência de alguns sinais, dentre eles os sinais com períodos de 3,5 e 0,7 anos (que provavelmente estão associados aos fenômenos ENSO) e o período de 0,33 anos (que pode estar relacionado a ressonâncias do sinal anual). A assinatura desses sinais apresentou variabilidade latitudinal. / A study was performed with physical data at the surface produced by the implementation of the numerical model POM and chlorophyll a data from the MODIS sensor, in the Tropical and South Atlantic, with the aim of to analyze the spatial variations of their correlations. The hydrodynamic model was subjected to a validation with data of Project PIRATE and MODIS sensor, obtaining satisfactory results, especially at the surface. The chlorophyll a values have high correlation with physical variables at the surface, but these values are influenced by the presence of annual and semiannual signals in the time series. Then, the analyses were carried out after removal of these signals from the series. The cross-correlations between physical variables and chlorophyll a had high magnitudes showing, in general, characteristic spatial pattern of negative correlations within the Subtropical Gyre and positive in other regions, for surface temperature and salinity and sea surface level. This spatial pattern is most evident in surface elevation and temperature, presenting lag in the biological response time nearly zero within the Subtropical Gyre and Equator, with a predominance of delays of the responses in other regions. The cross-covariance is calculated for several points along the meridian of 20º W for surface elevation and surface temperature, verifying agreement with the results obtained by the cross correlation. Spectral analysis of this covariance shows the occurrence of some signals, including signals with periods of 3.5 and 0.7 years (which are probably associated with ENSO phenomena) and the period of 0.33 years (that can be related to resonances of the annual signal). The signature of these signals showed latitudinal variability.

Submitted by Celso Magalhaes (celsomagalhaes@ufrrj.br) on 2017-08-29T12:36:36Z No. of bitstreams: 1 2017 - Iris Cristiane Magistrali.pdf: 5586295 bytes, checksum: cf90d85c09628780240eede5db4ec99e (MD5) / Made available in DSpace on 2017-08-29T12:36:36Z (GMT). No. of bitstreams: 1 2017 - Iris Cristiane Magistrali.pdf: 5586295 bytes, checksum: cf90d85c09628780240eede5db4ec99e (MD5) Previous issue date: 2017-02-20 / Coordena??o de Aperfei?oamento de Pessoal de N?vel Superior - CAPES / The Intergovernmental Panel on Climate Change (IPCC) has published reports on issues related to global climate changes. As a result, the climatic suitability of various agricultural crops in different regions of the world may change. Sugarcane is an important agricultural species and it is grown in several regions of the world, mainly for the generation of sugar and ethanol. Given this, there is a need for studies that assess the consequences of climate change in this crop in future scenarios. Therefore, the present work has as objectives: 1) to evaluate the performance of two climatic models CCCma (Canadian Center for Climate Modeling and Analysis) and GFDL (Geophysical Fluid Dynamic Laboratory) in the baseline period (1961-2000), and in the Future scenario A1B (2046-2065); 2) to elaborate the zoning of sugarcane for the State of Rio de Janeiro due to the A1B scenario considered intermediate. For the evaluation of the models performance, descriptive statistical analysis and spatial interpolation (Kriging and Co- Kriging ordinary) was used based on temperature and rainfall variables in water balance of the crop for the elaboration of agro climatic zoning, calculated by Thornthwaite & Mather?s (1955) method. The zoning was elaborated based on the reclassified thematic maps of water deficit and temperature. We used images from SRTM (Shuttle Radar Topography Mission) radar for altimetry, and map of use and ground cover. The results referring to the descriptive measures in the different points of the model showed the average rainfall ranged from 1118.02 mm to 1447.47 mm for the CCCma model. However, for the GFDL model average rainfall was lower than the CCCma, ranging from 804.82 mm to 1013.72 mm. The temperature averages presented similar behaviour and presented little oscillation during the year, when compared with the rainfall. For the CCCma model, the minimum average temperature was 21.02 ?C and the maximum was 22.99 ?C. For the GFDL model, the mean minimum temperature was 20.62 ?C and the maximum was 23.47 ?C. The CCCma model presented the data closest to the National Oceanic & Atmospheric Administration (NOOA) reanalysis data, and the GFDL model underestimated most rainfall and temperature data. The average rainfall of the NOOA data ranged from a minimum to a maximum of 947.74 mm at 1689.77 mm, respectively, and the minimum and maximum average temperature of 19.62 ? C and 21.62 ? C, respectively. Through geostatistics, it was possible to verify that, for the past and future scenarios, the exponential transitive model presented in the majority of cases the least degree of spatial dependence (GDE), being therefore considered the best. When comparing the GDE of the two climate models, noticed that the CCCma presented the best geostatistical performance. Due to the high GDE values, it was not possible to interpolate the GFDL model data regarding temperature for the future scenario. Through the product MCD12Q1, it was possible to register relevant differences in rainfall and temperature. Urban areas in general have lower amounts of rainfall and higher temperatures compared to areas with forest occurrence. In relation to the zoning, the results showed that for the past conditions, as well as for future projections, there is no water restriction for the development of the sugarcane crop. The temperature showed an increase in the future scenario adopted and, consequently, expansion of areas suitable for cultivation / O Painel Intergovernamental sobre Mudan?as Clim?ticas (IPCC) tem publicado relat?rios sobre quest?es ligadas as altera??es no clima global. Em fun??o disso, a aptid?o clim?tica de diversas culturas agr?colas em diferentes regi?es do mundo pode sofrer altera??es. A cana-de-a??car ? uma esp?cie agr?cola importante e, ? cultivada em v?rias regi?es do mundo, principalmente para a gera??o de a??car e etanol. Diante disso, existe a necessidade de estudos que avaliem as consequ?ncias das mudan?as clim?ticas nessa cultura em cen?rios futuros. Portanto, o presente trabalho tem como objetivos: 1) avaliar o desempenho de dois modelos clim?ticos CCCma (Canadian Centre for Climate Modelling and Analysis) e GFDL (Geophysical Fluid Dynamic Laboratory) no per?odo ?baseline? (1961-2000), e no cen?rio futuro A1B (2046-2065); 2) elaborar o zoneamento da cana-de-a??car para o Estado do Rio de Janeiro em fun??o do cen?rio de emiss?o A1B considerado intermedi?rio. Para a avalia??o do desempenho dos modelos foi utilizado an?lise estat?stica descritiva e interpola??o espacial (Krigagem e Co-Krigagem ordin?ria) para as vari?veis temperatura e chuva. Para a etapa da elabora??o do zoneamento agroclim?tico procedeu-se o c?lculo dos balan?os h?dricos da cultura, pelo m?todo de Thornthwaite e Mather (1955). Baseados nos mapas tem?ticos reclassificados de defici?ncia h?drica e temperatura elaborou-se o zoneamento. Utilizou-se imagens provenientes do radar SRTM (Shuttle Radar Topography Mission) para a altimetria, e, mapa de uso e cobertura do solo. Os resultados referentes as medidas descritivas nos diferentes pontos do modelo, mostraram que a m?dia da chuva variou de 1118,02 mm ? 1447,47 mm, para o modelo CCCma. Entretanto, para o modelo GFDL a m?dia da chuva foi inferior ao CCCma, variando de 804,82 mm ? 1013,72 mm. As m?dias de temperatura apresentaram comportamento semelhante e apresentaram pouca oscila??o durante o ano, se comparada a chuva. Para o modelo CCCma a temperatura m?dia m?nima foi de 21,02 ?C e a m?xima de 22,99 ?C. Para o modelo GFDL temperatura m?dia m?nima foi de 20,62 ?C e a m?xima de 23,47 ?C. O modelo CCCma apresentou os dados mais pr?ximos aos dados de rean?lise do National Oceanic & Atmospheric Administration (NOOA), sendo que o modelo GFDL subestimou a maioria dos dados de chuva e temperatura. A chuva m?dia dos dados do NOOA variou entre a m?nima e m?xima de 947,74 mm ? 1689,77 mm, respectivamente, e a temperatura m?dia m?nima e m?xima de 19,62 ?C e 21,62 ?C, respectivamente. Atrav?s da geoestat?stica, foi poss?vel constatar que tanto para o cen?rio passado quanto o futuro, o modelo transitivo exponencial apresentou na maioria dos casos o menor grau de depend?ncia espacial (GDE), sendo portanto, considerado o melhor. Ao comparar o GDE dos dois modelos clim?ticos, verifica-se que o CCCma apresentou o melhor desempenho geoestat?stico. Devido ao altos valores de GDE n?o foi poss?vel interpolar os dados do modelo GFDL referente a temperatura para o cen?rio futuro. Atrav?s do produto MCD12Q1 foi poss?vel registrar diferen?as relevantes em rela??o a chuva e a temperatura. ?reas urbanas em geral apresentam menores quantidades de chuvas e maiores temperaturas se comparadas ? ?reas com ocorr?ncia de floresta. Em rela??o ao zoneamento os resultados mostraram que tanto para as condi??es do passado, quanto para proje??es futuras, n?o h? restri??o h?drica para o desenvolvimento da cultura da cana-de-a??car. A temperatura apresentou um aumento no cen?rio futuro adotado e consequentemente expans?o de ?reas aptas ao cultivo. Praticamente todo o estado apresenta ?reas aptas em rela??o a topografia, exceto as regi?es serranas. Atrav?s do zoneamento final, ix foi poss?vel constatar que no cen?rio futuro adotado haver? uma expans?o da cana-de-a??car em torno de 58%, podendo esse resultado estar aliado ao aumento da temperatura.

Numerical models are becoming fundamental tools to predict a range of complex problems faced by geotechnical and geo-environmental engineers. However, to render the model reliable for future predictions, the model input parameters must be determined with consideration of the scale effects. If there is a difference of scales between the observation and the model scales there are two possible ways to consider it: or models are constructed with elements of a size similar to that at which the data were measured, or some upscaling rules must be defined. In this context, this thesis focuses on upscaling of water flow and mass transport in a tropical soil by means of numerical, laboratory and field studies. This thesis is organized in four parts. First, the heterogeneity, correlation and cross-correlation between solute transport parameters (dispersivity, α, and partition coefficient, Kd) and soil properties are studied in detail. In this part, it is verified that the hydraulic conductivity (K) and solute transport parameters are highly heterogeneous, while soil properties are not. Spatial correlation of α, K, and statistically significant variables are studied, and it would probably improve the estimation only in a small-scale study, since the spatial correlation are only observed up to 2.5 m. This study is a first attempt to evaluate the spatial variation in the correlation coefficient of transport parameters of a reactive and a nonreactive solute, indicating the more relevant variables and the one that should be included in future studies. In the second part, scale effect on K, dispersivity and partition coefficient of potassium and chloride are studied experimentally by means of laboratory and field experiments. The purpose is to contribute to the discussion about scale effects on K, α and Kd and understanding how these parameters behave with the change in the scale of measurement. Results show that K values increases with scale, regardless of the method of measurement, except for the results obtained from double-ring infiltrometer tests. Dispersivity trends to increase exponentially with the sample height. Partition coefficient tends to increase with sample length, diameter and volume. These differences in the parameters according to the scale of measurement must be considered when these observations are later used as input to numerical models, otherwise the responses can be misrepresented. Third, stochastic analysis of three-dimensional hydraulic conductivity upscaling is performed using a simple average and the Laplacian-with-skin methods for a variety of block sizes based on real K measurements. In this part it is demonstrated the errors that can be introduced by using a deterministic upscaling using simple averages of the measured K without accounting for the spatial correlation. Results show that K heterogeneity can be incorporated in the daily practice of the geotechnical modeler. The aspects to consider when performing the upscaling are also discussed. Finally, the dependence of the exponent of the p-norm as a function of the block size is analyzed. In the last part, stochastic upscaling of hydrodynamic dispersion coefficient (D) and retardation factor (R) is performed using real data aiming to reduce the lack in experimental upscaling of reactive solute transport research. The enhanced macrodispersion coefficient approach is used to upscale the local scale hydrodynamic dispersion (D) and, as a novelty, the impact of heterogeneity of local dispersivity is also taken into account. To upscale retardation factor (R), a p-norm is used to compute an equivalent R. Uncertainty analyses are also performed and a good propagation of the uncertainties is achieved after upscaling. Simple upscaling methods can be incorporated to the modeling practice using commercial transport codes and properly reproduce de transport at coarse scale but may require corrections to reduce smoothing of the heterogeneity caused by the upscaling procedure. / Modelos numéricos estão se tornando ferramentas fundamentais para prever uma série de problemas complexos enfrentados por engenheiros geotécnicos e geoambientais. No entanto, para que o modelo seja confiável para previsões futuras, seus parâmetros de entrada devem ser determinados com a consideração do efeito da escala. Se há uma diferença de escalas entre a escala da observação e a escala do modelo, existem duas maneiras possíveis de considerá-la: ou constrói-se modelos com elementos de tamanhos semelhantes àqueles em que os dados foram medidos, ou definem-se algumas regras de mudança de escala. Neste contexto, esta tese enfoca a mudança de escala do fluxo de água e do transporte de massa em um solo tropical, por meio de estudos numéricos, laboratoriais e de campo. Esta tese é organizada em quatro partes. Em primeiro lugar, estudou-se em detalhe a heterogeneidade, a correlação e a correlação cruzada entre os parâmetros de transporte de soluto (dispersividade, α, e coeficiente de partição, Kd) e as propriedades do solo. Nesta parte, verificou-se que a condutividade hidráulica (K) e os parâmetros de transporte de soluto são altamente heterogêneos, enquanto as propriedades do solo não o são. A correlação espacial de α, K e das variáveis estatisticamente significativas foi estudada, e, provavelmente, melhoraria a estimativa apenas em um estudo em pequena escala, uma vez que a correlação espacial só foi observada até 2,5 m. Este estudo foi uma primeira tentativa de avaliar a variação espacial no coeficiente de correlação dos parâmetros de transporte de um soluto reativo e não-reativo, indicando as variáveis mais relevantes e as que devem ser incluídas em estudos futuros. Na segunda parte, o efeito de escala em K, na dispersividade e no coeficiente de partição de potássio e cloreto é estudado experimentalmente por meio de ensaios laboratoriais e de campo. O objetivo foi contribuir com a discussão sobre os efeitos de escala em K, α e Kd e entender como esses parâmetros se comportam com a mudança na escala da medição. Os resultados mostram que K aumenta com a escala, independentemente do método de medição. A dispersão tende a aumentar de maneira exponencial com a altura da amostra. O coeficiente de partição tende a aumentar tanto com o comprimento, quanto com o diâmetro e o volume da amostra. Essas diferenças nos parâmetros de acordo com a escala de medida devem ser consideradas quando essas observações são posteriormente usadas como entrada para modelos numéricos, caso contrário, as respostas podem ser mal representadas. Em terceiro lugar, uma análise estocástica tridimensional da mudança de escala da condutividade hidráulica foi realizada usando tanto média simples quanto o método Laplaciano-com-pele para vários tamanhos de blocos usando medidas K reais. Nesta parte, foram demonstrados os erros que podem ser introduzidos ao se usar métodos determinísticos de mudança de escala, usando médias simples das medições de K sem se considerar a correlação espacial. A aplicação das técnicas de mudança de escala mostra que a heterogeneidade de K pode ser incorporada na prática diária do modelador geotécnico. Os aspectos a serem considerados ao realizar a mudança de escala também foram discutidos. Finalmente, analisou-se a dependência do expoente da norma p em função do tamanho do bloco. Na última parte, uma aplicação de mudança de escala estocástica do coeficiente de dispersão hidrodinâmica (D) e do fator de retardo (R) foi realizada usando dados reais visando reduzir a falta pesquisas no tema de mudança de escala do transporte de soluto reativo. A mudança de escala do D foi feito usando o método de macrodispersão. O método da média simples baseado na norma p foi usado para executar a mudança de escala de R. A incerteza foi propagada satisfatoriamente. Métodos simples de mudança de escala podem ser incorporados à prática de modelagem usando programas comerciais, e reproduzir corretamente o transporte em escala grossa, mas podem exigir correções para reduzir o efeito suavizado da heterogeneidade causada pelo procedimento de mudança de escala.

La déformation actuelle et la sismicité en domaine intraplaque sont des phénomènes mal compris. En particulier le rôle de l'héritage structural sur les variations spatiales de sismicité, et donc d'aléa sismique, reste une problématique majeure à la fois académique et sociétale. Plusieurs études semblent montrer que l'adoucissement rhéologique créé lors de la formation de grandes structures tectoniques aujourd'hui héritées, est un facteur indispensable pour expliquer la localisation et la magnitude de la déformation actuelle. Cependant, la quantification de l'impact de cet l’héritage structural sur la déformation de la croûte supérieure est toujours mal contrainte voir inconnue.Ces travaux de thèse ont pour but de quantifier l'impact de l'héritage structural sur la déformation actuelle des premiers kilomètres de la croûte. Il s’articule autour de deux axes de recherche. Le premier est la mesure précise des taux de déformations dans la vallée du Saint Laurent, à l’Est du Canada, à l’aide du Global Positioning System (GPS). Cette zone présente un fort héritage structural associé au paléo-rift Iapetus (~700 Ma) et une forte densité de station GPS permettant la mesure des taux de déformation tant à l’extérieur qu’à l’intérieur de la zone d’héritage structural. Une des nouveautés de cette étude est la combinaison de stations GPS continues et de campagnes (réseau de campagne dédié, pour la première fois utilisé pour contraindre les taux de déformation dans cette zone). L’analyse des taux de déformation obtenus à l’aide du GPS a montré qu’avec des taux de déformation de 2 à 11 fois supérieurs aux régions limitrophes, la zone d’héritage structural concentre significativement la déformation de surface. La comparaison des taux de déformation GPS avec ceux prédits par les modèles de réajustement isostatique glaciaire (GIA, processus majeur actif dans la région) a mis en évidence que l’héritage structural amplifie les taux de déformation générés par le GIA d’un facteur 6 à 28.Le deuxième axe de recherche consiste à développer des modèles numériques thermo-mécaniques 2D prenant en compte l’adoucissement rhéologique associé à l’héritage structural. La quantification de l’amplification des taux de déformation de la croûte supérieure est effectuée pour différentes localisations de l'adoucissement au sein de la lithosphère. Les modèles sont d’abord réalisés pour un cas intraplaque générique, puis appliqués à la vallée du St Laurent. Une des nouveautés de cette étude est d’analyser les modèles pour une force fixée, représentative du bilan des forces agissant sur la lithosphère, et supposée être à l’équilibre avec la résistance intégrée de la lithosphère. Cela permet d’appréhender le comportement mécanique de la lithosphère pour un état non stationnaire. Pour les modèles intraplaques génériques, l’amplification des taux de déformation varie de 1 à 27 en fonction des paramètres testés (lithologie crustale, quantité de déformation héritée, géotherme et force). Ces amplifications modélisées sont en accord avec celles mesurées par le GPS ou la sismicité à grande échelle spatiale (100s km) dans l’est de l’Amérique du Nord. Les modèles numériques associés à la vallée du St Laurent intègrent l’impact de la géométrie de la zone adoucie et des variations latérales de la quantité d’adoucissement et du géotherme. Les amplifications associées vont de 1 à 11, en accord avec les amplifications GPS calculées dans cette zone. Un résultat majeur mis en évidence par cette étude est la présence d’un noyau élastique conservé dans le manteau lithosphérique supérieur malgré le fort adoucissement imposé, limitant les effets d’amplification.Nos résultats reflètent un comportement mécanique de la lithosphère et de la déformation associée, supposés long-terme (10^3 à 10^5 ans). Des travaux restent à entreprendre pour quantifier l’impact des processus transitoires pouvant être à l’origine de la réactivation des zones d’héritage structural. / Current deformation and seismicity are poorly understood phenomena in intraplate domains. In particular, the role of structural inheritance on spatial variations of seismicity, and hence seismic hazard, remains a major social and academic issue. Previous studies tend to show that rheology weakening created during the formation of large tectonic structures, currently inherited, is essential to explain the current deformation localization and magnitude. However, quantifying the structural inheritance impact on upper crustal deformation is still poorly constrained, if not unknown.The goal of this thesis is to quantify the impact of structural inheritance on present-day upper crustal deformation. It is performed through two research axes. The first one is the measure of precise strain rates in the Saint Lawrence Valley, eastern Canada, using the Global Positioning System (GPS). This region presents a strong structural inheritance related to the Iapetus paleo-rift (~700 Ma) and a dense GPS network, allowing a high resolution in the strain rate measurements both outside and inside the structural inheritance area. One of the novelties in this study is the combination of continues and campaign GPS stations (dedicated network, used for the first time to constrain strain rates in this region). GPS strain rate analysis show that structural inheritance significantly localizes the deformation, increasing 2-11 times the average surface strain rates in the structural inheritance area compared to surrounding regions. The comparison of GPS strain rates and glacial isostatic adjustment models (GIA, the main active process in the St Lawrence Valley) shows that GIA-induced strain rates are amplified inside the inheritance area about a factor of 6 to 28, whereas they are not amplified outside the inherited domain.The second research axe involves developing 2D thermo-mechanical numerical models including rheology weakening related to structural inheritance. Quantification of the upper crustal strain rate amplification is carried out following different weak zone localization in the lithosphere. Models are first produced for a generic intraplate domain and then, dedicated to the St Lawrence Valley. One of the novelties of this study is to analyze the models for a fixed net driving force (including the impact of regional and local forces) assumed to be at equilibrium with the lithospheric integrated strength. This approach allows studying the mechanical behavior of the lithosphere for a non-steady state. For a generic intraplate model, the strain rate amplification varies between 1 to 27 depending on the tested parameters (crustal strength, amount of inherited strain, geotherm and net driving force). These model amplification factors are in good agreement with observed GPS and seismic strain rates calculated on a large spatial scale (100s km) in Eastern North America. Representatives models of the St Lawrence Valley investigate the impact of the weak zone geometry and lateral variations of both the amount of weakening and the geotherm. The associated amplification range between 1-11 and is in good agreement with GPS amplifications calculated in this region. A major feature highlighted through this study is the presence of a preserved elastic core in the upper lithospheric mantle despite the high imposed weakening, limiting amplification effects.Our results reflect long-term behavior (10^3 to 10^5 yrs) of the lithosphere deformation and mechanical behavior. The next step is to include and quantify the impact of transient processes, which could be at the origin of structural inheritance reactivation.

A erosão costeira é um processo de grande poder destrutivo e que afeta pelo menos 70% da zona costeira no mundo. Entender as causas e os processos que levam à erosão, bem como as regiões mais ou menos sensíveis a ela, é essencial para o correto manejo costeiro. O presente trabalho visa identificar as áreas vulneráveis à erosão através da aplicação de um Índice de Vulnerabilidade Costeira em praias do norte do Espírito Santo (Pontal do Ipiranga, Conceição da Barra) e sul da Bahia (Mucuri, Nova Viçosa, Caravelas, Prado, Cumuruxatiba, Corumbau, Arraial d\'Ajuda, Porto Seguro). O cálculo do índice se deu pela análise de treze indicadores ambientais: número de frentes frias, regime de tempestade, força de ondas, ângulo de incidência das ondas, estimativa de deriva potencial, morfologia da praia, exposição às ondas, presença de rios e/ou desembocaduras, elevação do terreno, vegetação, taxa de ocupação, obras de engenharia costeira e evidências de erosão. Os resultados mostraram que os indicadores \"morfologia da praia\" e \"exposição às ondas\" parecem ser os principais contribuintes para a erosão costeira nos municípios em que foram observados grandes prejuízos na infraestrutura urbana e perigo iminente à população local. Já nos demais municípios que também apresentaram focos de erosão, a combinação entre os indicadores \"estimativa de deriva potencial\" e \"regimes de tempestade\" parece ser a principal influência ao processo erosivo. Com exceção de Pontal do Ipiranga (que obteve um resultado de baixa vulnerabilidade), todas as praias apresentaram um IVC correspondente à vulnerabilidade moderada à erosão. O estudo se mostra eficiente para a determinação do nível de vulnerabilidade relativo entre as praias estudadas. / Coastal erosion is a process with great destructive power and affects more than 70% of worldwide coastlines. In order to achieve an effective coastal management, it is important to find and understand not only the causes and process that leads to erosion, but also areas that are more or less vulnerable to it. This study aims to identify erosion vulnerable areas through application of a Coastal Vulnerability Index in Brazilian beaches in north Espírito Santo (Pontal do Ipiranga, Conceição da Barra) and south Bahia (Mucuri, Nova Viçosa, Caravelas, Prado, Cumuruxatiba, Corumbau, Arraial d\'Ajuda, Porto Seguro). The index was calculated using 13 indicators: number of cold fronts, storm impact scale, wave power, angle of wave approach, longshore drift estimation, beach morphology, wave exposure, presence of inlets or river mouths, terrain elevation, vegetation, occupation rate, coastal engineering structures and erosion evidences. In places where great urban infrastructure destruction and imminent danger to local community were observed, both \"beach morphology\" and \"wave exposure\" seem to be the main contributors to coastal erosion. In other locations where erosion evidences were also observed, the combination between the indicators \"longshore drift estimation\" and \"storm impact scale\" seems to be the main influence to erosive process. IVC results showed that all the beaches are in a moderate vulnerability context, except in Pontal do Ipiranga, where the result was for low vulnerability. The study has proven to be effective in defining the relative vulnerability of the studied beaches.

Les zones de convergence sont des objets clés à la compréhension de la dynamique de la lithosphère. Elles sont le siège de déformations majeures comme en témoignent la concentration et l’intensité des séismes qui leur sont associées. A plus grande échelle de temps et d'espace, ces déformations résultent généralement en un empilement d’unités dont l’étude offre l'accès aux différents niveaux structuraux mis en jeu dans l'organisation de la chaîne de montagnes. Caractériser précisément la dynamique de structuration de ces unités ainsi que les paramètres contrôlant cette dynamique constitue une étape cruciale permettant d’aller plus loin, notamment quantitativement, dans notre connaissance de la dynamique lithosphérique. Cette thèse a comme double objectif (1) de préciser le comportement rhéologique des roches aux conditions Pressions-Température P-T du faciès éclogitique au sein de l' interface de subduction et (2) de caractériser et quantifier l' influence de la rhéologie de la plaque supérieure, et plus spécifiquement, le rôle de sa partie crustale, sur l'évolution spatio-temporelle des zones de convergence. Pour cela, une approche multi-disciplinaire a été utilisée. Dans un premier temps, je présente une analyse multi-échelle couplant travail de terrain et de pétrologie métamorphique qui a permis d'étudier la déformation au sein de roches HP-BT à l' interface de subduction au niveau de la klippe du Mont-Emilius (Alpes occidentales, Italie). Je détaille ensuite une étude quantitative alliant modélisation numérique thermo-mécanique 3D et 2D de zones de convergence. L'ensemble des modèles a permis d'analyser de nombreux paramètres influençant la structure rhéologique de la plaque supérieure tels que le géotherme initial, l'épaisseur de la lithosphère et de la croûte et la nature des matériaux impliqués. L'ensemble des modèles réalisés sont contraints/confrontés par/aux des données issues d'exemples naturels. Les résultats de l'étude sur des roches déformées au sein de l’interface de subduction mettent en évidence le possible comportement cassant des roches à des conditions de pression et température de l'ordre de 2.15-2.40 GPa, 500-550 °C, i.e., dans le faciès éclogitique. L’enregistrement d’un tel mode de déformation est d'une importance capitale car il remet en question le paradigme d’un comportement ductile sans résistance au niveau de l’ interface de subduction. Les résultats obtenus via les modèles numériques montrent par ailleurs que la rhéologie de la plaque supérieure, ainsi que celle de sa seule partie crustale, a une influence de premier ordre sur la dynamique globale des zones de convergence en modifiant notamment le mode de subduction, la cinématique de la fosse, le mode d'exhumation lors d'une collision, le timing de la déchirure du slab et de la formation de bassins d'arrière-arc, la répartition et l'intensité de la déformation au sein de la plaque supérieure. La combinaison des méthodes de pétrologie et de modélisation numérique a permis d'obtenir une analyse quantifiée de l' influence de la rhéologie des lithosphères impliquées dans les zones de convergence sur la dynamique de ces zones. Cette thèse présente ainsi de nouvelles contraintes à notre compréhension de la réponse mécanique de la lithosphère, en fonction de sa structuration rhéologique et de sa place au sein des zones de convergence à petite et grande échelle. Les nouvelles données présentées révèlent l' impact majeur de la rhéologie de la lithosphère dans les zones de convergence. Ce paramètre amène ainsi à reconsidérer notre vision actuelle des zones de convergence. / Convergence zones are key objects to the understanding of the lithosphere dynamics. They are the location correspond to places of intense deformation as evidenced by the concentration and magnitude of recorded earthquakes. On a larger scale of time and space, these deformations generally result in nappes stacking whose study offers access to the different structural levels involved in the mountain belt structuration. Precisely characterizing the structuring dynamics of these units as well as the parameters controlling this dynamic is a crucial step that would allow in particular a better quantification of lithospheric dynamics. The aim of this thesis is twofold: (1) clarifying the rheological behavior of rocks in the Pressure-Temperature (P-T) conditions of the eclogitic facies at the subduction interface and (2) characterizing and quantifying the influence of the overriding plate rheology, and more specifically, the role of its crustal part, on the spatio-temporal evolution of convergence zones. For this, I used a multi-disciplinary approach. First, I present a multi-scale analysis combining fieldwork and metamorphic petrology, which allowed me to study the deformation within High Pressure-Low Temperature (HP-LT) rocks at the subduction interface in the Mont-Emilius klippe (Western Alps, Italy). Then, I show the results of a quantitative study combining 3D and 2D thermo-mechanical modeling of convergence zones. The entire set of models allowed me to analyze different parameters influencing the rheological structure of the overriding plate, such as the initial geotherm, the thicknesses of the lithosphere and the crust, and the nature of the involved materials. All the performed models are constrained/compared by/with data from natural examples. The results of the study on deformed rocks within the subduction interface highlight the possible brittle behavior of rocks at pressure and temperature conditions on the order of 2.15-2.40 GPa and 500-550 °C, i.e., in the eclogitic facies. The recording of such a deformation mode is of paramount importance because it challenges the paradigm of subduction interface caracterized by ductile behavior without resistance. The results obtained with the numerical models show that the rheology of the overriding plate, as well as that of only its crustal part, has a first-order influence on the overall dynamics of the convergence zones by modifying the mode of subduction, trench kinematics, the mode of exhumation during collision, the timing for slab break-off and back-arc basin formation, the location and intensity of deformation within the overriding plate. The combination of petrology and numerical modeling methods allowed me to obtain a quantified analysis of the influence of the rheology of the lithospheres involved in convergence zones on the dynamics of these zones. This thesis presents new constraints for our understanding of the mechanical response of the lithosphere at different spatial scales as a function of its rheological structure. The new data presented here reveal the major impact of the lithosphere rheology in convergence zones. This parameter leads us to reconsider our current view of the convergence zones.

This dissertation investigates unexplained vapor intrusion field data sets that have been observed at hazardous waste sites, including: 1) non-linear soil gas concentration trends between the VOC source (i.e. contaminated groundwater plume) and the ground surface; and, 2) alternative pathways that serve as entry points for vapors to infiltrate into buildings and serve to increase VOC exposure risks as compared to the classic vapor intrusion model, which primarily considered foundation cracks as the route for vapor entry. The overall hypothesis of this research is that theoretical knowledge of fate and transport processes can be systematically applied to vapor intrusion field data using a multiple lines of evidence approach to improve the science-based understanding of how and when vapor intrusion exposure risks will pose increased exposure risk; and, ultimately this knowledge can be used to develop policies that reduce exposure risks. The first objective of this research involved numerical modeling, field sampling and laboratory tests to investigate which factors influence soil gas transport within the subsurface. Combining results of all of these studies provide improved understanding of which factors influence VOC fate and transport within the subsurface. Importantly, the results demonstrate a non-linear trend between the VOC source concentration in the subsurface and the ground surface concentration at the study site, which disagrees with many vapor intrusion conceptual models. Ultimately, the source concentration may not be a good predictor of shallow soil gas concentrations. Laboratory tests described the effect of soil characteristics such as the soil water content on VOC vapor diffusion. The numerical model was able to explain specific conditions that could not be described by the field and laboratory data alone. A paper was published that summarizes the major outcomes from this objective (Pennell et al, 2016). The second objective of this research investigated preferential pathways for VOC vapor migration into buildings. Sewer systems can act as important pathways for vapor intrusion. The research objective is to evaluate conditions that increase the potential for inhalation exposure risks via vapor intrusion thorough sewer systems into indoor spaces. A field study was conducted in California over a 4-year period to investigate the spatial and temporal variability of alternative pathways (e.g. aging infrastructure piping systems) within the context of vapor intrusion exposure risks. A paper was published that summarizes the major outcomes from the field study (Roghani et al. 2018). The final research objective involved the development of a numerical model to describe VOC fate and transport within a sewer system. The numerical model predicts VOC mass transport. The model results were compared to the field data and provides insight about the role preferential pathways play in increasing VOC exposure risks.