La cristallisation sous tension (SIC) du caoutchouc naturel (NR) a fait l’objet d’un nombre considérable d’études depuis sa découverte il y a près d’un siècle. Cependant, il existe peu d’informations dans la littérature concernant le comportement du caoutchouc à des vitesses de sollicitation proches des temps caractéristiques de cristallisation. L’objectif de cette thèse est alors de contribuer à la compréhension du phénomène de cristallisation sous tension grâce à des essais dynamiques à grandes vitesses. Pour répondre à cet objectif, nous avons développé une machine de traction permettant de déformer des échantillons d’élastomères à des vitesses de sollicitation pouvant aller jusqu’à 290s-1. Les essais ont été réalisés sur quatre NR avec des taux de soufre variables, deux NR chargés comportant des taux de noir de carbone différents. Nous avons également étudié un matériau synthétique à base de polyisoprène (IR) afin de comparer ses performances à celle du NR. Les essais dynamiques étant relativement difficiles à interpréter, un travail conséquent a donc été d’abord réalisé à basse vitesse. En outre, l’approche expérimentale proposée a été couplée à une approche thermodynamique de la SIC. Les mécanismes généraux associés à la cristallisation que nous identifions sont les suivants: lors d’une traction, la cristallisation consiste en l’apparition de populations cristallines conditionnée par l’hétérogénéité de réticulation des échantillons. Cette cristallisation semble nettement accélérée dès lors que ce cycle est réalisé au-dessus de la déformation de fusion. Nous attribuons ce phénomène à un effet mémoire dû à un alignement permanent des chaînes. Enfin, l’effet de la vitesse est décrit théoriquement en intégrant un terme de diffusion des chaînes dans la cinétique de SIC. Cette approche couplée à des essais mécaniques suggère que la SIC est essentiellement gouvernée par la cinétique de nucléation. Lors des tests dynamiques, la combinaison de l’effet mémoire et d’une accélération de la fusion pendant le cycle entraine une nette diminution voire une disparition de l’hystérèse cristalline. En outre, l’auto-échauffement, qui augmente progressivement avec la fréquence du cycle, tend à supprimer l’effet mémoire en provoquant le passage du cycle en dessous de la déformation de fusion. Lors de ces essais dynamiques, la SIC semble favorisée pour le matériau le moins réticulé. Nous attribuons cet effet au blocage d’enchevêtrements jouant le rôle de sites nucléants pour la SIC. Le matériau chargé semble avoir une moins bonne aptitude à cristalliser à hautes vitesses, par rapport à l’élastomère non chargé, en raison d’un auto-échauffement important à l’interface entre charges et matrice. Enfin, nous notons une convergence des cinétiques de cristallisation du caoutchouc naturel et synthétique à grande déformation et grande vitesse de sollicitation, que nous attribuons à la prédominance du terme énergétique d’origine entropique dans la cinétique de nucléation. / Strain induced crystallization (SIC) of Natural Rubber (NR) has been the subject of a large number of studies since its discovery in 1929. However, the literature is very poor concerning the study of SIC when samples are deformed with a stretching time in the range of the SIC characteristic time (around 10msec-100msec). Thus, the aim of this thesis is to contribute to the understanding of the SIC phenomenon thanks to dynamic tensile tests at high strain rates. To meet this goal, we have developed a dynamic tensile test machine allowing stretching samples of elastomers at strain rates up to 290 s-1. The tests are carried out on four NR with different sulphur amount, two NR with different carbon black filler amounts. We also studied a synthetic rubber made of polyisoprene chains (IR) able to crystallize under strain. Dynamic tests are relatively difficult to interpret; a significant work has thus been first performed at slow strain rate. Moreover, the experiments are coupled with a thermodynamic approach. First, the general mechanisms associated to the crystallization are identified as follows: during mechanical loading or during cooling in the deformed state, SIC is the result of successive appearance of crystallite populations whose nucleation and growth depend on the local network density. Crystallization is enhanced when the cycle is performed above the melting stretching ratio. This phenomenon is attributed to a memory effect due to a permanent alignment of the chains. Finally, the effect of the strain rate is theoretically described thanks to a diffusion term. This approach, coupled with experiments suggests that SIC is mainly governed by the nucleation kinetics. For the dynamic test, the combination of the memory effect and the acceleration of the melting during the cycle lead to a reduction or even disappearance of the crystalline hysteresis. In addition, self-heating, which progressively increases with the frequency of the cycle, causes the delay of the melting stretching ratio. This well explains why the crystallinity index decreases at the minimum stretching ratio of the dynamic cycles when the frequency increases. We finally compared the ability of our different rubbers to crystallize at high strain rates. SIC is enhanced for the weakly crosslinked rubber. This might be related to the dynamics of its free entanglements, these ones acting as supplementary crosslinks at high strain rates. Then, a filled rubber is compared to the unfilled one. We found that the filled sample has a lower ability to crystallize at high strain rates as compared to the unfilled one. This is likely due to the strong self-heating at the interface between the fillers and the rubbery matrix. Finally, we observe a convergence of crystallization kinetics in natural and synthetic rubbers at high strains and high strain rates. This is attributed to the predominance of the entropic energy in the nucleation kinetics in these experimental conditions.

La cristallisation induite par déformation des élastomères est un processus cinétique qui conduit à la formation de nano-cristallites thermodynamiquement stables. La présence de ces nano-cristallites au sein de la phase amorphe modifie considérablement les propriétés mécaniques des élastomères cristallisables. Ces élastomères ont en effet la propriété intéressante d'être auto-renforçants. L’objectif de ce travail est de développer un modèle physique capable de décrire localement l’évolution de la microstructure sous l’effet d’un champ de contrainte élastique durant la cristallisation. Dans ce but, un modèle de champ de phase est élaboré et mis en œuvre dans le cadre de la mécanique des milieux continus en couplant thermodynamique et mécanique avec une dynamique de transition de phase d’Allen-Cahn. La description thermodynamique de la cristallisation induite par déformation à petite échelle est basée sur la fonctionnelle d’énergie libre du système amorphe-cristal. Les conséquences du choix de cette formulation sont discutées, on étudie en particulier les effets de contraintes élastiques sur l’équilibre des phases en volumes ainsi que sur la cinétique de croissance des domaines cristallins au sein de l’amorphe. L’introduction de l’élasticité du réseau des contraintes topologiques induite par les enchevêtrements et/ou les nœuds de réticulation dans le modèle de champ de phase a permis de mettre en évidence l’existence d’un état stable de cristallites formées (modèle énergétique) mais aussi des instabilités de croissance (modèle cinétique). Sur la base de ces deux modèles, cinétique et énergétique, nous avons étudié systématiquement l’influence des contraintes topologiques sur la cinétique de croissance et nous montrons que cette cinétique est en effet contrôlée par l’accumulation de contraintes élastiques à l’interface. La prise en compte de l’élasticité du réseau des contraintes topologiques dans l’approche thermodynamique de cristallisation prédit une augmentation de la tension de surface et par conséquent un arrêt du mécanisme de croissance en donnant lieu à la formation de cristallites stables. Enfin, nous avons adopté le modèle énergétique pour modéliser le couplage entre nucléation, croissance et déformation cyclique. Pour valider le modèle local proposé une comparaison entre les résultats des simulations par champ de phase et les données expérimentales issues de la caractérisation d’un caoutchouc naturel réticulé est effectuée et nous montrons qualitativement l’accord entre l’expérience et le modèle. / Natural rubber NR and more generally elastomer presents unique physical properties that are very important for many engineering applications. Strain induced crystallization of elastomer presents a major interest because it improves considerably the mechanical properties. In fact, the presence of crystallites within the amorphous phase in a polymer network induces a strengthening of this material, giving NR a self-reinforcement character. In this thesis, we develop a mesoscopic model to describe the crystallization of elastomers under strain. In this context, we present a kinetic model using a new physical approach: a phase field model. This model combines the crystallization thermodynamics with the local stress field. The thermodynamic description of the phase transition is based on a Gibbs free energy functional F which contains all energy contributions of the system: the bulk contributions (enthalpy and entropy) and surface tension. To understand the experimental observation of nanometer size crystalites, an explicit account of the topological constraints induced by both entanglements and/or crosslinks is necessary. We investigated two limiting mechanisms, a kinetic limitation of the growth, and an energetic limitation. Based on both the kinetic and the energetic approaches, we have systematically studied the influence of topological constraints on the growth process. We have shown that the growth process is affected by the accumulation of elastic stress at the interface. The kinetic model predicts the existence of instabilities during the growth. These instabilities induce a heterogeneous dynamical growth which leads to the formation of dendrite like structures. On the contrary, the energetic approach predicts an exponential increase of the surface tension during the growth that limits the size of the crystallites very efficiently. In the last part we investigated elastomer crystallization under cyclic deformation. To this end, we coupled the previous energetic model with the nucleation process. Finally the simulation data are compared with experimental measurements.

Les processus hydrodynamiques et sédimentaires subtidaux sont étudiés en Manche au niveau du détroit du Pas-de-Calais et de son littoral méridional, et de la Mer d'Iroise et de la Baie de Douarnenez, à travers diverses modalités d'influence de l'hétérogénéité spatiale des sédiments de fond, de la houle via son interaction avec le courant en couche limite de fond et du vent. Une procédure nouvelle d'interpolation spatiale des distributions granulométriques permet de respecter l'hétérogénéité du substrat. L'approche méthodologique combine et confronte mesures de terrain et modélisation numérique. La modélisation, couplant des codes de circulation 3D (COHERENS) et de propagation de houle (SWAN), incorpore des modules de transport sédimentaire non-cohésif multiclasse, un module spécifique d'interaction houle-courant et diverses paramétrisations liées aux rides sableuses de fond (géométrie, coefficient de frottement), ainsi qu'un calcul optimisé du dépôt à l'interface eau-sédiment. La comparaison des prédictions numériques aux mesures recueillies le long du littoral méridional du détroit du Pas-de-Calais montre les interactions entre la bande littorale de sables fins envasés et la variabilité spatio-temporelle de l'hydrodynamique en conditions de courant de marée seul et superposé à une houle de Nord. La modélisation est ensuite appliquée à l'ensemble de la Manche orientale et centrale. En Baie de Douarnenez, l'analyse locale débouche sur une modélisation des influences respectives des courants de marée seuls, d'une houle de tempête superposée aux courants de marée et du vent. Enfin, une étude théorique des effets d'un cap sur la formation des bancs en drapeau est présentée.

Variabilité spatiale multiéchelle du zooplancton dans un lagoon récifal côtier - L'identification des changements dans les patrons écologiques selon l'échelle spatiale et la compréhension des processus qui génèrent ces changements sont d'une importance considérable en océanographie. Dans ce contexte, comprendre comment une communauté biologique répond à l'hétérogénéité environnementale requiert la connaissance des processus impliqués et l'échelle spatiale à laquelle ils opèrent. Les relations spatiales entre la variabilité du zooplancton et l'hétérogénéité environnementale sont encore imprécises dans les écosystèmes tropicaux côtiers. L'objectif de ce travail de thèse a donc été de déterminer les échelles de dépendance spatiale des patrons du zooplancton associé à un lagon récifal côtier et des processus environnementaux sous-jacents. Dans ce contexte, les intérêts de cette recherche ont été de quantifier les patrons de la variabilité du zooplancton dans l'espace, de comprendre comment ces patrons changent avec l'échelle spatiale et de déterminer les processus physiques et biologiques responsables de ces patrons spatiaux. L'échantillonnage, effectué le long d'un transect de la côte vers le large dans le lagon du Grand-Cul-de-Sac Marin (Guadeloupe), a concerné deux classes de taille du zooplancton (190–600 µm et > 600 µm) pour lesquelles la biomasse et l'abondance ont été estimées. L'analyse multiéchelle a été utilisée pour caractériser les patrons du zooplancton aux différentes échelles spatiales (de l'échelle de l'habitat à celle du lagon tout entier) et pour identifier les processus responsables de ces structures spatiales. Cette étude a montré que la variabilité du zooplancton est un phénomène multiéchelle dont l'amplitude et la dépendance spatiale dépendent de la taille des organismes, de leur motilité et de la variable-réponse considérée (biomasse ou abondance). La biomasse et l'abondance du zooplancton varient en réponse à la distribution spatiale du phytoplancton, au comportement du zooplancton, à l'hétérogénéité de l'habitat, à l'hydrodynamique et aux évènements météorologiques. La nature et les effets de ces processus sont dépendants de l'échelle spatiale. Ce travail a montré comment le changement d'échelle spatiale met en évidence différents niveaux d'organisation de la communauté zooplanctonique en réponse à l'hétérogénéité environnementale.

La caractérisation d'un aquifère est classiquement réalisée à partir d'études géologiques, piézométriques, géochimiques et de pompages d'essai. Ces techniques d'étude présentent certaines limites, qui peuvent néanmoins être réduites par l'utilisation de méthodes complémentaires comme certaines méthodes géophysiques. L'objectif de cette thèse est de quantifier l'amélioration apportée par les informations géophysiques dans la caractérisation hydrogéologique des aquifères. Deux contextes géologiquessont sélectionnés pour les difficultés qu'ils posent à la caractérisation hydrogéologique, l'importance des questions sociétales qui s'y posent et leur pertinence dans le cadre des interventions humanitaires : les aquifères de socle cristallin et les lentilles d'eau douce des cordons sableux côtiers. Les résultats d'une étude numérique et leur validation sur deux sites expérimentaux en Inde montrent que les informations complémentaires apportées par la technique de panneau électrique et la méthode de sondage par résonance magnétique des protons (RMP) permettent de préciser l'interprétation de pompages d'essai en contexte de socle. L'incertitude sur la transmissivité est réduite en moyenne de 39 % et celle sur le coefficient d'emmagasinement de 34 %. Cette approche hydrogéophysique permet ainsi de mieux décrire le milieu souterrain, et de proposer des modèles conceptuels plus complets. En milieu côtier sableux dans le Sud de l'Inde, l'utilisation combinée d'observations piézométriques, de la loi d'Archie et des méthodes de sondage électromagnétiques en domaine temporel (TDEM) et RMP permet de quantifier un volume d'eau douce disponible de 510 litres en moyenne par mètre carré de surface au sol avant la mousson dans la partie inhabitée (contre 420 L/m2 dans la partie habitée), et 670 L/m2 après la mousson (contre 450 L/m2 dans la partie habitée). Ces volumes auraient été surestimés de près de 40 % si une approche basée uniquement sur le TDEM et la loi d'Archie avait été utilisée. Le suivi temporel des mesures géophysiques permet également d'estimer la recharge nette de l'année hydrologique 2008 à 140 mm sur ce cordon (pour 2030 mm de précipitations). De plus, ce cordon sableux présente un champ géomagnétique hétérogène, qui rend la mesure de signaux de Free induction decay, (FID, classiquement enregistrés en sondage RMP) difficile. Une nouvelle procédure de mesure et d'interprétation basée sur des signaux d'écho de spin a donc été testée et validée expérimentalement. La caractérisation du milieu avec ce nouveau protocole permet de préciser la teneur en eau du milieu, sous-estimée de 26 % par le protocole FID. Ce résultat ouvre de nouvelles perspectives dans la recherche d'une relation quantifiée entre la teneur en eau RMP et les différentes porosités hydrogéologiques.

Cette thèse porte sur les membranes biologiques, étudiées du point de vue de la physique théorique. Nous nous intéressons aux effets génériques de la présence d'une ou deux inclusions membranaires, par exemple des protéines, et à ceux d'une modification chimique locale de l'environnement de la membrane.<br /> Tout d'abord, nous étudions une interaction entre deux inclusions membranaires, qui est analogue à la force de Casimir : elle provient des contraintes que les inclusions imposent aux fluctuations thermiques de la forme de la membrane. Nous calculons les fluctuations de cette force entre deux inclusions ponctuelles. Nous clarifions la définition de la force exercée par un fluide corrélé sur une inclusion, dans un micro-état du fluide. Cette définition joue un rôle clé dans l'étude de la force de Casimir au-delà de sa valeur moyenne à l'équilibre. Nous étudions également les interactions de Casimir entre des inclusions membranaires de forme allongée.<br /> Ensuite, nous nous intéressons à l'élasticité membranaire à l'échelle nanométrique, qui est mise en jeu dans les déformations locales de l'épaisseur de la membrane à proximité de certaines protéines. Nous soulignons l'importance d'un terme énergétique qui a été négligé jusqu'à présent.<br /> Enfin, nous présentons une description théorique, développée à partir de principes fondamentaux, de la dynamique d'une membrane soumise à une perturbation chimique locale de son environnement. Nous comparons nos prévisions théoriques à de nouveaux résultats expérimentaux portant sur la déformation dynamique d'une membrane biomimétique soumise à une augmentation locale de pH.

Contexte: Les champignons mycorhiziens à arbuscules (AMF) établissent des relations symbiotiques avec la plupart des plantes grâce à leurs réseaux d’hyphes qui s’associent avec les racines de leurs hôtes. De précédentes études ont révélé des niveaux de variation génétique extrêmes pour des loci spécifiques permettant de supposer que les AMF peuvent contenir des milliers de noyaux génétiquement divergents dans un même cytoplasme. Si aucun processus de reproduction sexuée n’a jusqu’ici été observé chez ces mycorhizes, on constate cependant que des niveaux élevés de variation génétique peuvent être maintenus à la fois par l’échange de noyaux entre hyphes et par des processus fréquents de recombinaison entre noyaux. Les AMF se propagent par l’intermédiaire de spores qui contiennent chacune un échantillon d’une population initiale de noyaux hétérogènes, directement hérités du mycélium parent. À notre connaissance les AMF sont les seuls organismes qui ne passent jamais par un stade mononucléaire, ce qui permet aux noyaux de diverger génétiquement dans un même cytoplasme. Ces aspects singuliers de la biologie des AMF rendent l’estimation de leur diversité génétique problématique. Ceci constitue un défi majeur pour les écologistes sur le terrain mais également pour les biologistes moléculaires dans leur laboratoire. Au-delà même des problématiques de diversité spécifique, l’amplitude du polymorphisme entre noyaux mycorhiziens est mal connue. Le travail proposé dans ce manuscrit de thèse explore donc les différents aspects de l’architecture génomique singulière des AMF. Résultats L’ampleur du polymorphisme intra-isolat a été déjà observée pour la grande sous-unité d’ARN ribosomal de l’isolat Glomus irregulare DAOM-197198 (précédemment identifié comme G. intraradices) et pour le gène de la polymerase1-like (PLS) de Glomus etunicatum isolat NPI. Dans un premier temps, nous avons pu confirmer ces résultats et nous avons également pu constater que ces variations étaient transcrites. Nous avons ensuite pu mettre en évidence la présence d’un goulot d’étranglement génétique au moment de la sporulation pour le locus PLS chez l’espèce G. etunicatum illustrant les importants effets d’échantillonnage qui se produisaient entre chaque génération de spore. Enfin, nous avons estimé la différentiation génétique des AMF en utilisant à la fois les réseaux de gènes appliqués aux données de séquençage haut-débit ainsi que cinq nouveaux marqueurs génomiques en copie unique. Ces analyses révèlent que la différenciation génomique est présente de manière systématique dans deux espèces (G. irregulare et G. diaphanum). Conclusions Les résultats de cette thèse fournissent des preuves supplémentaires en faveur du scénario d’une différenciation génomique entre noyaux au sein du même isolat mycorhizien. Ainsi, au moins trois membres du genre Glomus, G. irregulare, G. diaphanum and G. etunicatum, apparaissent comme des organismes dont l’organisation des génomes ne peut pas être décrit d’après un modèle Mendélien strict, ce qui corrobore l’hypothèse que les noyaux mycorhiziens génétiquement différenciés forment un pangenome. / Background: Arbuscular mycorrhizal fungi (AMF) are root-inhabiting fungi whose hyphal networks form symbioses with plants. Previous studies have revealed extremely high levels of genetic variation for some loci, which has lead to the proposition that AMF contain thousands of genetically divergent nuclei that share the same cytoplasm, i.e. they are heterokaryotic coenocytes. No reproductive stage has as yet been observed in AMF, yet evidence is accumulating that the observed high levels of diversity could be maintained by the exchange of nuclei between hyphal systems and (meiotic) recombination. AMF spores contain varying fractions of this heterogeneous population of nuclei, which migrate directly from the parent mycelium. To our knowledge, AMF are the only organisms that never pass through a single nucleus stage in their life cycle, which allows nuclei to diverge into genetically distinct nuclei within the same cytoplasm. Thus, estimating genetic diversity in arbuscular mycorrhizal fungi (AMF) is a major challenge, not only for ecologists in the field but also for molecular biologists in the lab. It is unclear what the extent of polymorphism is in AMF genomes. The present thesis investigates different aspects of this peculiar genome organization. Results The second chapter in this thesis confirms the extensive intra-isolate polymorphism that was previously observed for large subunit rDNA (in G. irregulare DAOM-197198) and the polymerase1-like gene, PLS (in G. etunicatum), and shows that this polymorphism is transcribed. In the third chapter I report the presence of a bottleneck of genetic variation at sporulation for the PLS locus, in G. etunicatum. Analyses in the fourth chapter, based on a conservative network-based clustering approach and five novel single copy genomic markers, reveal extensive genome-wide patterns of diversity in two different AMF species (G. irregulare and G. diaphanum). Conclusions The results from this thesis provide additional evidence in favor of genome differentiation between nuclei in the same isolate for AMF. Thus, at least three members of the Glomus genus, G. irregulare, G. diaphanum and G. etunicatum appear to be organisms whose genome organization cannot be described by a single genome sequence: genetically differentiated nuclei in AMF form a pangenome.

No description available.

No description available.

Les bassins d’infiltration font partie intégrante des techniques alternatives de gestion des eaux pluviales en milieu urbain. Néanmoins, la potentialité de transfert de polluants vers la nappe est accrue en cas d’écoulements préférentiels dans les sols sous-jacents. Une bonne compréhension du couplage entre processus d’écoulements préférentiels en zone vadose et mécanismes géochimiques est requise. Cette thèse s’inscrit dans le cadre du suivi d’un bassin d’infiltration depuis plusieurs dizaines d’années de fonctionnement. Le site d’étude est situé sur le dépôt fluvioglaciaire hétérogène couvrant une grande partie de la région lyonnaise. Des auscultations sur une fosse sous le bassin (section 13.5m*2.5m) ont mis en évidence une régionalisation particulière de la pollution dans le sol. Cette étude s’appuie sur une étude numérique visant à identifier l’origine de la régionalisation des polluants et à la relier aux écoulements préférentiels résultant des hétérogénéités lithologiques. En amont de l’étude numérique, les lithofaciès sont complètement caractérisés aux regards de leurs propriétés hydrodynamiques, hydrodispersives et géochimiques. La modélisation numérique permet de souligner l’établissement de cheminements préférentiels en lien avec le contraste de propriétés hydrodynamiques, notamment lorsque de faibles débits sont appliqués en surface. Le rôle de chaque lithofaciès et de l’architecture du dépôt (stratification et inclusions) est clairement identifié. Les répercussions de tels écoulements sur les transferts non réactifs sont ensuite investiguées en combinant l’influence des écoulements préférentiels et le fractionnement de l’eau en fractions mobile et immobile résultant de l’hétérogénéité intrinsèque au sein de chaque lithofaciès. Enfin, ces processus physiques sont couplés à la réactivité géochimique pour le cas d’un polluant modèle (le cuivre) en prenant en compte la réactivité différentielle des lithofaciès. Ces résultats permettent de générer un modèle conceptuel d’écoulements préférentiels et de transfert de masse en milieu fortement hétérogène. / An infiltration basin is a stormwater best management practice (BMP) designed to infiltrate runoff volumes in urban areas. Nevertheless, preferential flow paths in the underlying soil may cause rapid migration of pollutants, thus contributing to groundwater contamination. Understanding the coupling between preferential flow processes in the vadose zone and geochemical mechanisms is then required. This thesis is a part of the follow-up of an infiltration basin for several decades of exploitation. The study site was settled over a highly heterogeneous glaciofluvial deposit covering much of the Lyon region. The investigation of an excavated section of the basin (13.5m long and 2.5m deep) pointed out a specific regionalization of pollution in the soil. This research is based on a numerical study to identify the origin of such a pollutant pattern and link this with preferential flow resulting from lithological heterogeneities. Different lithofacies were fully characterized regarding their hydraulic, hydrodispersive and geochemical properties. The numerical study proves that the high contrast in hydraulic properties triggers the establishment of preferential flow (capillary barriers and funneled flow). Preferential flow develops mainly for low initial water contents and low fluxes imposed at surface. The role of each lithofacies and architecture of deposit (stratification and inclusions) is clearly identified. The impact of such flows on non-reactive transfers is then investigated by combining the influence of preferential flow and pore water fractionation info into mobile and immobile fractions, resulting from the intrinsic heterogeneity within each lithofacies. Finally, these physical processes are coupled to the geochemical reactivity for a pollutant model (copper), taking into account the differential reactivity of lithofacies. These results generate a conceptual model of preferential flow and mass transfer in strongly heterogeneous media.