Mars est entourée d'une atmosphère ténue, composée à 95% de dioxyde de carbone (CO2). Au cours d'une année martienne, des transferts de masse (jusqu'à 30% du CO2 atmosphérique) entre l'atmosphère et les calottes polaires produisent des variations temporelles à très grande longueur d'onde du champ de gravité, notamment des harmoniques zonaux de son développement en harmoniques sphériques (de fait les coefficients ”composites” de degré 2 et 3). D'un autre côté, le potentiel gravitationnel du Soleil induit des déformations, dites de marée, du volume martien. Ces déformations produisent un potentiel perturbateur en tout point extérieur à la planète, proportionnel à son nombre de Love de degré 2 k2. k2 traduit la réponse élastique de la planète au potentiel solaire et permet de caractériser physiquement le noyau de Mars (sa nature, solide ou liquide, et son rayon). Une manière de quantifier les transferts de la masse atmosphérique et l'état du noyau est de déterminer les perturbations inhérentes sur le mouvement d'un satellite artificiel. Le cycle saisonnier du CO2 et l'état du noyau impliquent aussi des variations de la rotation de Mars. Une autre manière de quantifier les transferts de la masse atmosphérique et l'état du noyau est donc d'observer leurs effets sur la rotation. Des simulations d'observations de trajectographie de satellites (comme celles de Mars Global Surveyor, MGS, Odyssey, MODY) et/ou de la position d'un réseau de stations à la surface de Mars (comme dans l'expérience NEIGE) nous ont permis de voir s'il est possible de restituer précisément les variations des harmoniques zonaux de gravité de bas degré et/ou la rotation. Avec les observations réelles de trajectographie des missions américaines MGS et MODY, on a restitué les variations des harmoniques zonaux de gravité de bas degré et k2.

En 2009, le Grenelle de l’environnement a ouvert en France la possibilité d’intégrer des clauses environnementales au sein des cahiers des charges d’Indications Géographiques (IG). Cette opportunité nouvelle intéresse fortement les Agences de l’eau qui cherchent de plus en plus à agir en amont à la source des pollutions agricoles. Dans ce contexte, cette étude interroge la manière dont certaines démarches de qualité alimentaire intègrent dans leur dynamique des objectifs de préservation de l’eau. Pour y répondre, nous avons développé une approche globale permettant de resituer les enjeux de qualité de l’eau dans ceux plus larges d’environnement. Nous avons également dépassé la seule analyse des cahiers des charges des IG pour considérer l’ensemble des actions en faveur de l’environnement dans lesquelles elles s’impliquent. Le cadre d’analyse développé, celui de la résilience des Systèmes Alimentaires (SA) de qualité différenciée répond à ce besoin d’élargissement. Il nous permet d’étudier les dynamiques environnementales engagées pour faire face à des perturbations liées à la dégradation des ressources naturelles dont l’eau. Il s’agit autant de perturbations factuelles, mesurées, que de perturbations représentées lorsque les acteurs envisagent un risque potentiel. Ce cadre d’analyse a été appliqué à la viticulture et à l’arboriculture dans le Sud-Est de la France. Les résultats montrent la complexité des dynamiques environnementales engagées : elles reposent toutes sur plusieurs actions complémentaires et supposent la présence d’acteurs clés ayant la capacité à engager une dynamique collective et à envisager les SA à la fois dans leur dimension verticale les liant au marché et dans leur dimension horizontale les insérant dans le territoire. / In 2009, The Grenelle de l'Environnement (French environmental pact) introduced environmental clauses within Geographical Indications (GI) specifications. French Water agencies (agences de l’eau) perceive it as a good opportunity to address the issue of agricultural pollution, by takings actions at the upstream level. In this context, this study examines how certainfood quality processes are proactively setting targets for conserving water. This question led to develop a global approach in order to look at the issues of water quality from the wider angle of environmental preservation. It requires not only to focus on regulation tools of food quality products but also to take into account a larger set of environmental actions. Along the same lines, the theoretical framework of the resilience of food system of differentiated quality has also been developed. It allows to study proactive environmental dynamics set up by stakeholders to deal with perturbations linked to natural resources degradation, such as water. Those perturbations can be factual (measured) or represented, when stakeholders consider a potential risk. This analytical framework has been applied on case studies in viticulture and fruit production in the South-East of France. The results underscore the complexity of environmental dynamics: they are based on a set of complementary actions and require the mobilization of key players, able to engage in a collective dynamics but also to have a wider vision of food systems, both by looking at the vertical dimension (market oriented) and horizontal dimension (link with the local territory).

L’objectif de cette thèse est d’analyser le développement des instabilités secondaires bidimensionnelles et tridimensionnelles dans les couches de mélange à densité variable, incompressibles et à nombre de Froude infini. Dans ces conditions, la présence d’inhomogénéités de masse volumique modifie sensiblement la dynamique rotationnelle de l’écoulement et celle des instabilités secondaires sous l’action du couple barocline. Une analyse de stabilité linéaire non-modale est mise en oeuvre pour identifier les mécanismes physiques de croissance transitoire. Cette analyse permet également de prendre en compte le caractère instationnaire de la couche de mélange, absent dans l’analyse modale quasi-statique de Fontane (2005). Après établissement des équations de Navier–Stokes linéarisées directes et adjointes à densité variable, celles-ci sont utilisées dans une méthode d’optimisation itérative qui permet de déterminer les perturbations à croissance énergétique maximale. La première partie consiste en la description des perturbations optimales pour une couche de mélange homogène. Aux temps courts, lorsque la couche de mélange est quasi-parallèle, les perturbations optimales présentent de fortes amplifications transitoires, dont l’origine physique est due à la synergie des mécanismes classiques de Orr et de lift-up. Puis lorsque la couche s’enroule pour former un tourbillon de Kelvin–Helmholtz, les perturbations évoluent vers les instabilités tridimensionnelles elliptiques ou hyperboliques, selon le nombre d’onde latéral. Dans la deuxième partie, l’analyse est étendue aux couches de mélange à densité variable. Pendant la phase initiale de développement des perturbations optimales, les inhomogénéïtés de masse volumique ont une influence minime sur la croissance des perturbations. Ce n’est qu’une fois la couche de mélange enroulée que les effets de densité deviennent actifs, entraînant un supplément d’amplification significatif par rapport à la situation homogène. En particulier, le couple barocline favorise le développement des perturbations du côté du fluide léger du rouleau de Kelvin–Helmholtz. Enfin, lorsque le temps d’injection des perturbations est suffisamment retardé, la vorticité produite par le couple barocline favorise le développement d’une instabilité bidimensionnelle du type Kelvin-Helmholtz identifiée par Reinaud et al. (2000). / The purpose of this thesis is to analyse the development of two-dimensional and three-dimensional secondary instabilities in incompressible variable-density mixing layers, in the limit of infinite Froude number. Under these conditions, mass inhomogeneities alter significantly the rotational dynamics of the flow under the action of the baroclinic torque. A nonmodal stability analysis is implemented to identify the physical mechanisms of transient growth. This analysis allows to take into account the unsteady natureof the flow, which was absent in the quasi-static modal analysis (Fontane, 2005). After establishing of the direct and adjoint linearised Navier-Stokes equations for variable-density flows, they are used in an iterative optimization method to determine the perturbations that maximize their energy. The optimal perturbations are first obtained for a homogeneous time-evolving mixing layer. For times short enough, when the time-evolving mixing layer is almost parallel, optimal perturbations exhibit the largest transient growth. These amplifications arise from the synergy between the well-known Orr and liftup mechanisms. Once the mixing layer rolls up into a Kelvin–Helmholtz billow, the disturbances trigger the three-dimensional elliptical and hyperbolic instabilities. The analysis is then extended to variable-density mixing layers. During the initial development of optimal perturbations, mass inhomogeneities have no influence over the perturbations growth. Once the mixing layer has rolled up, the variable-density effects contribute significantly to the increase of the perturbation energy. In particular, the baroclinic torque enhances the development of perturbations in the light side of the Kelvin–Helmholtz billow. Finally, when the injection time of perturbations is delayed long enough, the baroclinic vorticity generation on the light side of the Kelvin–Helmholtz billow triggers a two-dimensional secondary Kelvin–Helmholtz instability, which has been identified by Reinaud et al. (2000).

Les nombres aléatoires ont été de tous temps utilisés pour des jeux de hasard, plus récemment pour créer des codes secrets et ils sont aujourd'hui nécessaire à l'exécution de programmes informatiques. Les générateurs de nombres aléatoires sont maintenant bien éloignés de simples dés à lancer et sont constitués de circuits électroniques ou d'algorithmes. Ceci pose des problèmes quant à la reconnaissance du caractère aléatoire des nombres générés. De plus, de la même manière ou autrefois les dés étaient pipés pour augmenter les chances de gagner, il est aujourd'hui possible d'influencer la sortie des générateurs de nombres aléatoires.Ce sujet est donc toujours d'actualité avec des exemples récents très médiatisés. Ceci concernait en effet la console de jeu PS3 qui génère un nombre aléatoire constant où la distribution de clefs secrètes redondantes sur internet.Ce mémoire présente l'étude de plusieurs générateurs ainsi que diverses manières de les perturber. Il montre ainsi des faiblesses inhérentes à leurs conceptions et des conséquences possibles de leur défaillance sur des composants de sécurité. Ces travaux ont de plus permis de mettre en évidence l'importance des problématiques concernant le test des nombres aléatoires ainsi que des retraitements corrigeant des biais dans ces nombres. / Random numbers have been used through the ages for games of chance, more recently for secret codes and today they are necessary to the execution of computer programs. Random number générators have now evolved from simple dices to electronic circuits and algorithms. Accordingly, the ability to distinguish between random and non-random numbers has become more difficult. Furthemore, whereas in the past dices were loaded in order to increase winning chances, it is now possible to influence the outcome of random number generators.In consequence, this subject is still very much an issue and has recently made the headlines. Indeed, there was talks about the PS3 game console which generates constant random numbers and redundant distribution of secret keys on the internet.This thesis presents a study of several generators as well as different means to perturb them. It shows the inherent defects of their conceptions and possible consequences of their failure when they are embedded inside security components. Moreover, this work highlights problems yet to be solved concerning the testing of random numbers and the post-processing eliminating bias in these numbers distribution.

L'un des enjeux actuel de la neutronique concerne la propagation rigoureuse des incertitudes d'entrée (e.g. données nucléaires, tolérances de fabrications, etc.) aux résultats finaux calculés par les codes (e.g. keff, taux de réaction, etc.). Pour propager les incertitudes, il est de coutume de faire l'hypothèse de petites variations autour d'une référence et de calculer, dans un premier temps, des profils de sensibilités. Or, les codes Monte-Carlo, qui se sont imposés comme des références de calcul, ne possèdent pas -- ou n'ont intégré que très récemment -- un moyen direct de calculer des sensibilités et donc de réaliser un calcul précis d'incertitudes. Les approches déterministes, elles, permettent le calcul de ces sensibilités mais introduisent parfois de très fortes hypothèses, notamment sur la géométrie.Le premier objectif de se travail de thèse est d'introduire dans le code Monte Carlo du CEA de transport des neutrons, TRIPOLI-4, des méthodes à même de calculer des profils de sensibilités du keff aux données nucléaires ou à toute autre perturbation. Pour cela, il a d'abord été nécessaire de mettre en place le calcul du flux adjoint d'un milieu critique. Pour la première fois, et grâce aux développements informatiques de ce travail, il a été possible de calculer dans un cas réel, concret, et applicatif, des spectres de flux adjoints en un point quelconque d'un réacteur. Ceci a été réalisé à l'aide de la probabilité itérée de fission (Iterated Fission Probability ou IFP) qui assimile le flux adjoint à l'importance d'un neutron dans un réacteur exactement critique. Ce calcul de flux adjoint a, par la suite, ouvert la porte au premier développement d'une méthode de calcul de « perturbations exacte » en Monte Carlo, théorie qui permet de s'affranchir des hypothèses de petites variations, et qui ouvre la porte à certaines applications jusqu'alors difficiles à analyser.Au delà de l'analyse poussée de la méthode IFP et de son application au calcul de flux adjoint, cette thèse propose également, d'obtenir dans le code TRIPOLI-4, les paramètres cinétiques d'un réacteur pondérés par le flux adjoint ou bien des aires de migration. A cette fin, l'implémentation reprend et améliore un algorithme déjà développé par la communauté scientifique pour estimer des perturbations au premier ordre. / A key issue in nowadays Reactor Physics is to propagate input data uncertainties (e.g. nuclear data, manufacturing tolerances, etc.) to nuclear codes final results (e.g. keff, reaction rate, etc.). In order to propagate uncertainties, one typically assumes small variations around a reference and evaluates at first sensitivity profiles. Problem is that nuclear Monte Carlo codes are not -- or were not until very recently -- able to straightforwardly process such sensitivity profiles, even thought they are considered as reference codes.First goal of this PhD thesis is to implement a method to calculate keff-sensitivity profiles to nuclear data or any perturbations in TRIPOLI-4, the CEA Monte Carlo neutrons transport code. To achieve such a goal, a method has first been developed to calculate the adjoint flux using the Iterated Fission Probability (IFP) principle that states that the adjoint flux at a given phase space point is proportional to the neutron importance in a just critical core after several power iterations. Thanks to our developments, it has been made possible, for the fist time, to calculate the continuous adjoint flux for an actual and complete reactor core configuration. From that new feature, we have elaborated a new method able to forwardly apply the exact perturbation theory in Monte Carlo codes. Exact perturbation theory does not rely on small variations which makes possible to calculate very complex experiments.Finally and after a deep analysis of the IFP method, this PhD thesis also reproduces and improves an already used method to calculate adjoint weighted kinetic parameters as well as reference migrations areas.

L’objectif de cette thèse est d’étudier et de donner des outils pour la compréhension des problèmes de perturbations singulières pour des modèles épidémiques et des problèmes de dynamiques de populations. Les modèles considérés sont des équations structurées en âge qui peuvent dans certains cas se réécrire comme des équations à retard. L’étude de ces classes d’exemples s’est faite avec succès et a permis de comprendre et de mettre en évidence toute la complexité et l’étendue de ces problèmes. Comme on peut le remarquer dans la littérature, l’une des clés fondamentales à la compréhension de ces problèmes est l’étude des variétés normalement hyperboliques en dimension infinie que nous avons largement étudiées dans cette thèse. L’approche utilisée est la méthode de Lyapunov-Perron. Ce qui nous a amené à étudier les problèmes de persistance et d’existence de trichotomie (dichotomie) exponentielle qui sont des éléments fondamentaux dans l’utilisation de cette méthode. / In this thesis we aim to give tools to understand singular perturbations in epidemic model sand population dynamic models. We study some singularly perturbed delay differential equation which does not enter into the class frame work of geometric singular perturbation for delay differential equations. An example of singularly perturbed age structured model is also studied. The study of these examples allowed us to understand and highlight some complexities of these problems. One of the main tools in understanding such questions is the normally hyperbolic manifolds theory which is our central focus in this thesis. The approach used here is the Lyapunov-Perron method. Therefore the problems of persistence and existence of exponential trichotomy (dichotomy) are also stressed since there are one of the mainingredients of this method.

Dans cette thèse, nous nous intéressons aux propriétés spectrales des opérateurs non-auto-adjoints aléatoires. Nous allons considérer principalement les cas des petites perturbations aléatoires de deux types des opérateurs non-auto-adjoints suivants :1. une classe d’opérateurs non-auto-adjoints h-différentiels Ph, introduite par M. Hager [32],dans la limite semiclassique (h→0); 2. des grandes matrices de Jordan quand la dimension devient grande (N→∞). Dans le premier cas nous considérons l’opérateur Ph soumis à de petites perturbations aléatoires. De plus, nous imposons que la constante de couplage δ vérifie e (-1/Ch) ≤ δ ⩽ h(k), pour certaines constantes C, k > 0 choisies assez grandes. Soit ∑ l’adhérence de l’image du symbole principal de Ph. De précédents résultats par M. Hager [32], W. Bordeaux-Montrieux [4] et J. Sjöstrand [67] montrent que, pour le même opérateur, si l’on choisit δ ⪢ e(-1/Ch), alors la distribution des valeurs propres est donnée par une loi de Weyl jusqu’à une distance ⪢ (-h ln δ h) 2/3 du bord de ∑. Nous étudions la mesure d’intensité à un et à deux points de la mesure de comptage aléatoire des valeurs propres de l’opérateur perturbé. En outre, nous démontrons des formules h-asymptotiques pour les densités par rapport à la mesure de Lebesgue de ces mesures qui décrivent le comportement d’un seul et de deux points du spectre dans ∑. En étudiant la densité de la mesure d’intensité à un point, nous prouvons qu’il y a une loi de Weyl à l’intérieur du pseudospectre,une zone d’accumulation des valeurs propres dûe à un effet tunnel près du bord du pseudospectre suivi par une zone où la densité décroît rapidement. En étudiant la densité de la mesure d’intensité à deux points, nous prouvons que deux valeurs propres sont répulsives à distance courte et indépendantes à grande distance à l’intérieur de ∑. Dans le deuxième cas, nous considérons des grands blocs de Jordan soumis à des petites perturbations aléatoires gaussiennes. Un résultat de E.B. Davies et M. Hager [16] montre que lorsque la dimension de la matrice devient grande, alors avec probabilité proche de 1, la plupart des valeurs propres sont proches d’un cercle. De plus, ils donnent une majoration logarithmique du nombre de valeurs propres à l’intérieur de ce cercle. Nous étudions la répartition moyenne des valeurs propres à l’intérieur de ce cercle et nous en donnons une description asymptotique précise. En outre, nous démontrons que le terme principal de la densité est donné par la densité par rapport à la mesure de Lebesgue de la forme volume induite par la métrique de Poincaré sur la disque D(0, 1). / In this thesis we are interested in the spectral properties of random non-self-adjoint operators. Weare going to consider primarily the case of small random perturbations of the following two types of operators: 1. a class of non-self-adjoint h-differential operators Ph, introduced by M. Hager [32], in the semiclassical limit (h→0); 2. large Jordan block matrices as the dimension of the matrix gets large (N→∞). In case 1 we are going to consider the operator Ph subject to small Gaussian random perturbations. We let the perturbation coupling constant δ be e (-1/Ch) ≤ δ ⩽ h(k), for constants C, k > 0 suitably large. Let ∑ be the closure of the range of the principal symbol. Previous results on the same model by M. Hager [32], W. Bordeaux-Montrieux [4] and J. Sjöstrand [67] show that if δ ⪢ e(-1/Ch) there is, with a probability close to 1, a Weyl law for the eigenvalues in the interior of the pseudospectrumup to a distance ⪢ (-h ln δ h) 2/3 to the boundary of ∑. We will study the one- and two-point intensity measure of the random point process of eigenvalues of the randomly perturbed operator and prove h-asymptotic formulae for the respective Lebesgue densities describing the one- and two-point behavior of the eigenvalues in ∑. Using the density of the one-point intensity measure, we will give a complete description of the average eigenvalue density in ∑ describing as well the behavior of the eigenvalues at the pseudospectral boundary. We will show that there are three distinct regions of different spectral behavior in ∑. The interior of the of the pseudospectrum is solely governed by a Weyl law, close to its boundary there is a strong spectral accumulation given by a tunneling effect followed by a region where the density decays rapidly. Using the h-asymptotic formula for density of the two-point intensity measure we will show that two eigenvalues of randomly perturbed operator in the interior of ∑ exhibit close range repulsion and long range decoupling. In case 2 we will consider large Jordan block matrices subject to small Gaussian random perturbations. A result by E.B. Davies and M. Hager [16] shows that as the dimension of the matrix gets large, with probability close to 1, most of the eigenvalues are close to a circle. They, however, only state a logarithmic upper bound on the number of eigenvalues in the interior of that circle. We study the expected eigenvalue density of the perturbed Jordan block in the interior of thatcircle and give a precise asymptotic description. Furthermore, we show that the leading contribution of the density is given by the Lebesgue density of the volume form induced by the Poincarémetric on the disc D(0, 1).

L'un des plus grands problèmes ouverts de la cosmologie moderne est l'origine de l'expansion accélérée de l'Univers, découverte en 1998. L'explication théorique la plus simple repose sur l'introduction d'une constante cosmologique $Lambda$. Ce modèle, connu sous le nom de $Lambda$CDM, est en accord avec les différentes observations liées à l'expansion accélérée. Cependant, il présente des problèmes d'ordre théorique. Par conséquent, plusieurs alternatives, connues collectivement sous le nom de {it modèles d'énergie noire}, ont été proposées pour expliquer cette accélération. Plusieurs d'entre eux restent viables, car leurs {it backgrounds} cosmologiques ne présentent pas de signatures identifiables. Par contre, les effets sur les phénomènes perturbatifs sont plus spécifiques à chacun de ces modèles. Dans cette thèse, nous explorons les caractéristiques particulières de la croissance des perturbations de matière à l'ordre linéaire dans les théories $f(R)$ avec un regard complémentaire sur les modèles chameleon. La paramétrisation du taux de croissance de la matière en termes d'une fonction $gamma$ permet d'identifier une signature très spécifique de ces modèles en comparaison avec le modèle $Lambda$CDM. Une étude supplémentaire a permis de trouver une dépendance en échelle explicite, nommée {it dispersion}, dans la croissance des perturbations. Des observations plus précises pourraient permettre de faire la différence entre ces différents modèles selon la présence de ces caractéristiques. / One of the most important open issues in modern cosmology is the origin of the accelerated expansion of the Universe, observed in 1998. The simplest theoretical explanation relies on the introduction of a cosmological constant $Lambda$. This model, known as $Lambda$CDM, agrees with all the different observations connected to the accelerated expansion. However, it presents some theoretical issues. As a result, several alternatives, known collectively under the name of {it dark energy models}, have been proposed to explain this acceleration. Several among them remain viable, since their cosmological backgrounds do not show any identifiable signature. On the other hand, effects on the perturbative level are more specific to each model. In this thesis, we explore the particular characteristics of the growth of linear matter perturbations in $f(R)$ theories, with a complementary look on chameleon models. The parameterization of the growth rate in terms of a $gamma$ function allows us to identify a very specific signature of these models in comparison with the $Lambda$CDM model. A subsequent study allowed us to find an explicit scale dependance, known as {it dispersion}, in the growth of perturbations. More precise observations could enable us to distinguish between dark energy models according to the presence of this type of feature.

Afin d'améliorer la qualité des produits et tendre vers des conceptions fiables et robustes, la simulation numérique joue de nos jours un rôle clé dans de nombreux secteurs de l'ingénierie. Malgré l'utilisation de modèles de plus en plus complexes et réalistes, les corrélations entre les mesures expérimentales et les simulations déterministes ne s'avèrent pas toujours évidentes, en particulier, si le phénomène observé est de nature fugace. Afin de prendre en compte les variations possibles de comportement, des techniques de tirages multiples comme les plans d'expériences, les analyses de sensibilité ou les approches non déterministes peuvent être exploitées. Cependant, ces simulations avancées conduisent inévitablement à des temps de calcul prohibitifs qui ne sont pas en adéquation avec des phases de conception de plus en plus courtes.L'objectif de cette thèse est d'explorer de nouvelles stratégies de résolution pour les problèmes mécaniques, où la non-linéarité de contact frottant et des variations sur les paramètres du modèle numérique sont considérés en même temps. Pour y parvenir, nous avons, dans un premier temps, étudié l'intégration de contrôleurs, basés sur la logique floue, pour résoudre un problème de contact frottant. L'idée proposée est de transformer le problème non linéaire en un ensemble de problèmeslinéaires de tailles réduites que l'on peut réanalyser grâce des développements homotopiques et des techniques de projection. Dans un second temps, nous avons étendu la démarche proposée au cas des problèmes de vibrations induites par le frottement comme le crissement. / To improve the quality of products and tend to reliable and robust designs, numerical simulations have nowadays taken a key role in many engineering domains. In spite of more complex and realistic numerical models, the correlation between a deterministic simulation and experimentations are not obvious, especially if the observed phenomenon have a fugitive nature. To take into account possible evolutions of behaviour, multiple samplings techniques such as designs of experiments, sensitivity analyses or non-deterministic approaches are currently performed. Nevertheless, these advanced simulations necessarily generate prohibitive computational times, which are not compatible with more and more shorter design steps.The aim of this work is to explore new numerical ways to solve mechanical problems including both the contact nonlinearity, the friction and several variability on model parameters. To achieve this objective, the integration of Fuzzy Logic Controllers has been first studied in the case of static frictional contact problems. The proposed idea is to decompose the non linear problem in a set of reduced linear problems. These last ones can be reanalyzed thanks to homotopy developments and projection techniques as a function of introduced perturbations. Second, the proposed strategy has been extended to the case of friction induced vibrations problems such as squeal.

De longs délais temporels peuvent exister entre les perturbations survenant à un site et les réponses des espèces de la communauté occupant ce site. Ces changements écosystémiques transitoires lents en comparaison d’une vie humaine peuvent être difficiles à détecter ou à comprendre sans la prise en compte du contexte temporel par l’utilisation de données historiques. C’est souvent le cas pour les processus d’invasion de plantes ou d’animaux non-indigènes dans les communautés forestières. Bien que les perturbations anthropiques soient reconnues pour faciliter les invasions dans ce type de communauté, des temps de latence importants peuvent tout de même survenir à chaque étape du processus. On peut donc penser que les patrons d’invasion d’un site représenteraient davantage les legs écologiques des perturbations plus anciennes que ceux des perturbations plus récentes. Ce mémoire vise à évaluer si l’on gagne à tenir compte de l’âge des perturbations pour expliquer les invasions biologiques établies dans les sites environnants. Pour vérifier cela, nous avons inventorié l’abondance et la composition des communautés de plantes et de vers de terre non-indigènes dans 85 parcelles dispersées dans l’ensemble de la Réserve Naturelle Gault. Suite à une reconstitution fouillée de l’historique des perturbations anthropiques, nous avons comparé le pouvoir explicatif des perturbations initiées avant 1910 avec celui des perturbations survenues après 1910. Nous avons aussi pris en compte les variables environnementales reconnues comme importantes. Les résultats suggèrent que dans les forêts de la Réserve Naturelle Gault, la présence et l’intensité des invasions de plantes et de vers de terre non-indigènes sont en premier lieu expliquées par la proximité aux perturbations anthropiques initiées il y a plus d’un siècle. Les conclusions de ce mémoire soulignent l’importance de l’utilisation d’une approche historique pour comprendre pleinement la durée des impacts écologiques des perturbations anthropiques quant à la facilitation des processus d’invasion biologique.