Résumé : En raison de sa grande étendue, le Nord canadien présente plusieurs défis logistiques pour une exploitation rentable de ses ressources minérales. La TéléCartographie Prédictive (TCP) vise à faciliter la localisation de gisements en produisant des cartes du potentiel géologique. Des données altimétriques sont nécessaires pour générer ces cartes. Or, celles actuellement disponibles au nord du 60e parallèle ne sont pas optimales principalement parce qu’elles sont dérivés de courbes à équidistance variable et avec une valeur au mètre. Parallèlement, il est essentiel de connaître l'exactitude verticale des données altimétriques pour être en mesure de les utiliser adéquatement, en considérant les contraintes liées à son exactitude. Le projet présenté vise à aborder ces deux problématiques afin d'améliorer la qualité des données altimétriques et contribuer à raffiner la cartographie prédictive réalisée par TCP dans le Nord canadien, pour une zone d’étude située au Territoire du Nord-Ouest. Le premier objectif était de produire des points de contrôles permettant une évaluation précise de l'exactitude verticale des données altimétriques. Le second objectif était de produire un modèle altimétrique amélioré pour la zone d'étude. Le mémoire présente d'abord une méthode de filtrage pour des données Global Land and Surface Altimetry Data (GLA14) de la mission ICESat (Ice, Cloud and land Elevation SATellite). Le filtrage est basé sur l'application d'une série d'indicateurs calculés à partir d’informations disponibles dans les données GLA14 et des conditions du terrain. Ces indicateurs permettent d'éliminer les points d'élévation potentiellement contaminés. Les points sont donc filtrés en fonction de la qualité de l’attitude calculée, de la saturation du signal, du bruit d'équipement, des conditions atmosphériques, de la pente et du nombre d'échos. Ensuite, le document décrit une méthode de production de Modèles Numériques de Surfaces (MNS) améliorés, par stéréoradargrammétrie (SRG) avec Radarsat-2 (RS-2). La première partie de la méthodologie adoptée consiste à faire la stéréorestitution des MNS à partir de paires d'images RS-2, sans point de contrôle. L'exactitude des MNS préliminaires ainsi produits est calculée à partir des points de contrôles issus du filtrage des données GLA14 et analysée en fonction des combinaisons d’angles d'incidences utilisées pour la stéréorestitution. Ensuite, des sélections de MNS préliminaires sont assemblées afin de produire 5 MNS couvrant chacun la zone d'étude en totalité. Ces MNS sont analysés afin d'identifier la sélection optimale pour la zone d'intérêt. Les indicateurs sélectionnés pour la méthode de filtrage ont pu être validés comme performant et complémentaires, à l’exception de l’indicateur basé sur le ratio signal/bruit puisqu’il était redondant avec l’indicateur basé sur le gain. Autrement, chaque indicateur a permis de filtrer des points de manière exclusive. La méthode de filtrage a permis de réduire de 19% l'erreur quadratique moyenne sur l'élévation, lorsque que comparée aux Données d'Élévation Numérique du Canada (DNEC). Malgré un taux de rejet de 69% suite au filtrage, la densité initiale des données GLA14 a permis de conserver une distribution spatiale homogène. À partir des 136 MNS préliminaires analysés, aucune combinaison d’angles d’incidences des images RS-2 acquises n’a pu être identifiée comme étant idéale pour la SRG, en raison de la grande variabilité des exactitudes verticales. Par contre, l'analyse a indiqué que les images devraient idéalement être acquises à des températures en dessous de 0°C, pour minimiser les disparités radiométriques entre les scènes. Les résultats ont aussi confirmé que la pente est le principal facteur d’influence sur l’exactitude de MNS produits par SRG. La meilleure exactitude verticale, soit 4 m, a été atteinte par l’assemblage de configurations de même direction de visées. Par contre, les configurations de visées opposées, en plus de produire une exactitude du même ordre (5 m), ont permis de réduire le nombre d’images utilisées de 30%, par rapport au nombre d'images acquises initialement. Par conséquent, l'utilisation d'images de visées opposées pourrait permettre d’augmenter l’efficacité de réalisation de projets de SRG en diminuant la période d’acquisition. Les données altimétriques produites pourraient à leur tour contribuer à améliorer les résultats de la TCP, et augmenter la performance de l’industrie minière canadienne et finalement, améliorer la qualité de vie des citoyens du Nord du Canada. / Abstract : Due to its vast extent, Northern Canada faces several logistical challenges for a profitable exploitation of its mineral resources. Remote Predictive Mapping (RPM) aims to help in targeting mineral deposits through the production of geological potential maps. Elevation data is necessary for the generation of these maps. However, the currently available elevation data north of the 60th parallel are not optimal primarily because it has been derived from contours with values at a metric precision. Additionally, exact knowledge of the vertical accuracy of elevation data is essential to insure a suitable use, within its accuracy constraints. This project aimed to improve the quality of elevation data and to contribute to the refinement of RPM products for a study site located in the Northwest Territories. The first objective was to generate control points to evaluate vertical accuracy with precision. The second objective was to generate an improved elevation model for the study site. First, a filtering method for Global Land and Surface Altimetry Data (GLA14) from the ICESat (Ice, Cloud and land Elevation SATellite) mission is presented. This filtering is based on indicators, derived from information available in GLA14 data and terrain conditions, which are then applied successively to remove potentially contaminated elevation points. The points are filtered based on the attitude calculation, signal saturation, equipment noise, atmospheric conditions, slope and number of peaks. Next, a method to generate an improved Digital Surface Models (DSM) using StereoRadarGrammetry (SRG) with Radarsat-2 (RS-2) images is described. In the first part of the adopted methodology, DSM are stereorestituted from RS-2 image pairs, without control point. Then, the vertical accuracy of the DSM is calculated using the control points resulting from the filtering of GLA14 data, and analysed according to the incidence angles combination used for the stereorestitution. Next, selections from the preliminary DSM are assembled to generate 5 DSM, each covering entirely the study site. Finally, the DSM are analysed to identify the optimal selection for the area of interest. The selected indicators were found to be efficient and complementary, except for the indicator based on the noise/signal ratio. Otherwise, all indicators allowed to filter out points exclusively. A 19% reduction of the elevation mean square error was achieved with the filtering method, when compared to Canadian Digital Elevation Data (CDED). The initial density of the GLA14 allowed maintaining a spatially homogeneous distribution of the post-filtering elevation points despite a 69% rejection rate. From the analysis of the 136 preliminary DSM, no specific combination of the acquired RS-2 images incidence angles stood out as being ideal with SRG due to high variability in vertical accuracy. Nonetheless, the analysis showed that images should be ideally acquired at sub-zero temperatures to minimize radiometric discrepancies between scenes. Results also showed that the slope is the main factor influencing the accuracy of DSM generated with SRG. The best vertical accuracy (4 m) was achieved with same-side view configurations. Opposite-side view configurations, despite achieving a vertical accuracy of 5 m, allowed a 30% reduction in the amount of images initially acquired. Therefore, the use of opposite-side view configurations could help to improve the efficiency of SRG projects by reducing considerably the acquisition period. Elevation data generated using the proposed method could help to improve results from RPM and increase the efficiency of the mining industry in Northern Canada and finally contribute to the betterment of the lives of Northern Canada’s citizens.

Un bon dimensionnement d’une fondation superficielle nécessite une bonne caractérisation du sol sur lequel elle reposera. Pour le cas des sols granulaires, les essais de reconnaissances destructifs (SPT, CPT, etc.) impliquent la conversion des profils obtenus du sol exploré en un profil de rigidité afin de pouvoir estimer sa portance. Une telle tâche peut induire des erreurs au niveau de l’estimation étant donné qu’elle nécessite le passage d’un nombre de paramètres (N, qc) n’illustrant pas l’état réel du sol (rigidité) à des paramètres caractérisant sa performance (module élastique, coefficient de Poisson). La caractérisation des sols granulaires en utilisant la vitesse des ondes de cisaillement (non destructive) nous permet de contourner ce problème étant donné qu’elle nous donne directement la variation d’un paramètre intrinsèque au sol (Vs) à partir duquel il est possible d’effectuer le dimensionnement. La disponibilité d’une technique de mesure de vitesses de cisaillement (P-RAT) sur des échantillons de sols au laboratoire géotechnique de l’UdeS, nous a permis de mettre en œuvre l’application de Vs dans l’estimation de la portance des sols granulaires. Au cours de ce mémoire, une procédure servant à estimer le tassement des fondations superficielles reposant sur des dépôts granulaires est présenté. Les paramètres de cisaillement (Vs, G) d’un sable déterminés à l’aide des essais de caractérisation (P-RAT) ont été utilisés pour construire un modèle numérique de fondation. La consistance des résultats numériques du modèle a été confirmée à travers des essais de chargement de plaque. Ce travail montre donc l’importance des paramètres de cisaillement dans la modélisation d’une application fréquente en géotechnique. Il constitue un départ considérable pour la construction de relations constitutives basées sur le paramètre vitesse de cisaillement.

L’étude présentée s’inscrit dans un programme de recherche mené par l’Université de Sherbrooke (UdeS) pour le ministère des Transports du Québec (MTQ) sur la détection de l’endommagement par des mesures de vibrations ambiantes sur un pont routier en service. La méthode s’appuie sur la mise à jour automatique des propriétés physiques d’un modèle numérique d’un pont de façon à reproduire les fréquences et modes propres de la structure mesurés expérimentalement. Pour cette étude, le logiciel commercial FEMtools a été utilisé. L’ouvrage ciblé est le pont Est de la Rivière-aux-Mulets situé à Sainte-Adèle (Qc) sur l’autoroute des Laurentides (autoroute 15) sur lequel des essais d’analyse modale ont été effectués sous vibrations ambiantes en novembre 2013. L’étude présentée porte essentiellement sur la phase d’ingénierie la plus délicate du processus de détection d’endommagement par mise à jour d’un modèle numérique, soit la définition et le recalage du modèle initial. Pour cela, un modèle du pont par éléments finis a été développé sous le logiciel Abaqus, puis a été exporté vers FEMtools selon une procédure qui a été établie. Avant d’utiliser le logiciel FEMtools pour le recalage automatique du modèle, une brève étude de vérification des capacités du logiciel à recaler un modèle et à détecter de l’endommagement est effectuée à partir d’un cas d’étude connu étudié antérieurement à l’UdeS (cadre de deux étages en béton armé endommagé sismiquement). Puisque les algorithmes de recalage automatique de modèles par éléments finis sont très sensibles à plusieurs variables et choix définis par l’utilisateur et que le problème de recalage est fortement sous-déterminé, le risque d’obtenir des résultats incohérents est grand si on ne suit pas une démarche rigoureuse. Une méthodologie en cinq étapes adaptée aux problématiques de pont est donc proposée et appliquée au cas du pont de la Rivière-aux-Mulets. Les résultats du recalage en termes de réductions de rigidité des sections de béton par rapport aux rigidités initiales supposées sont cohérents avec les comportements statique et dynamique attendus d’un pont routier en service. Pour en arriver à ce résultat, plusieurs tentatives de recalage ont été effectuées. Un échantillon de ces cas d’études est présenté, montrant notamment l’influence dans le processus de la définition et du nombre des paramètres de recalage, de leurs valeurs initiales, de la complexité des modes expérimentaux et du critère de convergence utilisé. Il est aussi montré que des résultats de recalage manifestement erronés permettent parfois d’identifier des erreurs majeures de modélisation, ce qui est un atout important de la méthode. Finalement, plusieurs scénarios d’endommagement fictifs imposés numériquement sont étudiés pour évaluer brièvement la capacité de FEMtools à détecter de l’endommagement sur le pont. Les résultats montrent notamment que la précision de la détection est grandement améliorée par l’utilisation d’un nombre croissant de modes expérimentaux.

Ce travail de recherche s'inscrit dans le cadre de l'analyse vibroacoustique des structures complexes. Il propose une nouvelle approche énergétique utilisant le champ de densité d'énergie afin de simplifier une telle analyse. Cette approche est basée sur un modèle numérique stochastique construit à l'aide de l'approche probabiliste non paramétrique des incertitudes de modélisation et de paramètres. L'approche énergétique stochastique développée correspond à une nouvelle représentation du système vibroacoustique en terme des grandeurs énergétiques aléatoires. Un modèle vibroacoustique énergétique moyen est alors construit en prenant la moyenne statistique des grandeurs énergétiques. On dispose alors d'un modèle énergétique moyen pour analyser la vibroacoustique des systèmes complexes dans la bande des basses et des moyennes fréquences alors que la méthode SEA ne permet pas d'analyser cette bande de fréquence. L'analyse des propriétés des grandeurs énergétiques moyennes utilisées pour la construction du modèle vibroacoustique énergétique permet de construire une version simplifiée conduisant à un modèle énergétique simplifié pour lequel une méthodologie de partition structurale par zone est établie. Une application de cette approche énergétique et de la méthodologie de partition structurale par zone est présentée pour une voiture constituée d'une structure automobile couplée avec sa cavité acoustique interne

Ces dernières années, les robots manipulateurs ont été utilisés pour des tâches, telles que le polissage, le fraisage, le placement des fibres, entre autres, auparavant réalisées par des machines spécialisées. Cependant, l’utilisation des robots pour ces tâches en contact, présente des problèmes de précision. Ce contact provoque des flexibilités dans le robot, l’outil et l’environnement. Cette thèse cible l’amélioration de la précision de la tâche robotisée de dépose de fibres qui nécessite un compactage homogène selon le matériau et la surface. La précision du drapage peut être améliorée par la correction en ligne de la trajectoire du robot et du point de centrage dynamique de l’outil (DTCP). Cette thèse propose l’asservissement en position/force-pression en modifiant le TCP, pour commander le compactage selon la surface de contact et prend en compte les flexibilités dues au contact. La pression et les flexibilités sont calculées à partir d’un modèle mécanique numérique réduit grâce à des réseaux de neurones. Cet asservissement a besoin d’une sélection appropriée des correcteurs pour garantir un bon fonctionnement et pour travailler avec divers types de moules. La sélection est réalisée dynamiquement par logique floue. Le DTCP est calculé dynamiquement grâce aux flexibilités de l’outil et il est redéfini lors de l’exécution de la trajectoire du robot pendant le drapage. Un dispositif expérimental est mis en oeuvre en utilisant des éléments similaires aux cellules robotisées de placement des fibres. Une application logicielle intègre l’asservissement proposé dans la baie de commande. Cet asservissement permet d’améliorer la précision de dépose de fibres. / In the last years, the robots have been used in the tasks that were previously performed by specialized machines such as polishing, milling, drilling, fibers lay up, among others. However, the use of industrial robots for these contact tasks presents problems of precision. This contact causes flexibility in the robot, in the tool and in the environment. This thesis is aimed at improving the accuracy of the robotic fiber placement task. This task requires a homogeneous and adapted compaction according to the material and the surface. The lay-up precision can be improved by online robot trajectory correction and the dynamic tool center point (DTCP) definition. This thesis presents a position/force-pressure control by modifying the TCP, to control the compaction according to the contact surface and which takes into account the flexibility due to the compaction force. Pressure and flexibility are calculated from a mechanical numerical model reduced using neural networks. The force-pressure control needs a selection of corrector to ensure a good performance and to work with several types of molds. This dynamic correction selection is performed by a fuzzy logic approach. The DTCP has been calculated dynamically thanks to the tool flexibility and it is redefined when the robot’s trajectory is running during the draping. An experimental setup has been implemented using industrial elements similar to robotic fiber placement cells. A software application has been developed which integrates the control loop into the robot controller. This loop improves the accuracy of fiber placement.

La thèse s'attache à l'étude des processus hydrodynamiques et hydro-sédimentaires ``cross-shore'' affectant la morphodynamique des plages sableuses dissipatives soumises à l'action des vagues, de la marée et des rechargements artificiels. La méthode choisie, qui s'appuie sur la modélisation numérique déterministe, couple le modèle de propagation de vagues FUNWAVE basé sur les équations complètement non linéaires de Wei et al. avec deux modèles de transport sédimentaire et un module d'évolution morphologique. Le premier modèle de transport, 1DH, intègre la vitesse près du fond corrigée par Lynett et la formule de transport sédimentaire total de Bailard. Le second modèle, pseudo-2DV, est une juxtaposition de modèles 1DV de turbulence exprimée en coordonnée sigma. Ce double outil est validé par des mesures de laboratoire (principalement dans le canal Delta Flume du Delft Hydraulics aux Pays-Bas) et in situ sur la plage de Pentrez en baie de Douarnenez (Finistère). Il reproduit bien la propagation des vagues, la distribution verticale des vitesses et des concentrations de sédiments en suspension, et de manière satisfaisante en dépit d'un lissage excessif l'évolution morphologique comme le mouvement des barres. La version 1DH, plus stable et moins coûteuse en temps calcul que la version pseudo-2DV, a été appliquée pour quantifier l'influence des différents paramètres sur la morphodynamique de plages planes soumises à l'action conjuguée des vagues et de la marée. Les résultats des simulations sont conformes au modèle descriptif de Masselink et Short. Finalement, ce modèle 1DH a été mis en oeuvre pour explorer l'évolution de rechargements artificiels et esquisser quelques recommandations préliminaires sur leur pratique.

Ce travail est consacré au modèle rapide du procédé de déposition par projection thermique en plasma d'arc soufflé. Il est plus particulièrement consacré au traitement d'une particule isolée ou d'un lot de particules échantillon représentatif d'une poudre.<br />Le procédé de projection par plasma d'arc soufflé et les principaux phénomènes qui régissent la formation du jet de plasma, les échanges thermiques et dynamiques entre la particule et l'écoulement du jet gazeux et la construction du dépôt sont examinés. Nous avons fait une étude bibliographique des modèles développées par ailleurs pour simuler les différentes fonctionnalités du procédé.<br />Les fondements et les caractéristiques du logiciel « Jets&Poudres » sont exposés, avant de présenter le modèle du transfert plasma-particule qui prend en compte la conduction interne à la particule et les déplacements des fronts de changement de phase. Dans ce modèle le calcul de l'évaporation de la particule est découplé du problème de la dynamique du gaz dans le jet et nous exploitons les résultats de J. C. Knight et le modèle de « Pression en retour » (Back Pressure) qu'il a développé. Une étude qualitative est consacrée à l'effet des paramètres de dispersion de la poudre en sortie de l'injecteur sur le traitement d'une particule isolée et sur la construction du dépôt. Ces effets mis en évidences, nous exposons un modèle de transport d'un lot de particules représentatif d'une poudre dans l'injecteur afin d'évaluer la dispersion en masse, taille et vitesse avant l'entrée dans le jet. Ce modèle complexe prend ent compte les collisions particule-parois et les colisions binaires particule- particule.<br />Les résultats de ces deux modèles sont discutés. Le traitement dynamique et thermique de la particule isolée est en bon accord avec ceux de la littérature. Ce qui autorise l'étude des conditions opératoires et des paramètres des différentes composantes fonctionnelles du procédé (torche, gaz plasmagène, injecteur, poudre...) sur le traitement des particules et leurs histoires thermique et dynamique.<br />Le modèle a été utilisé pour évaluer la quantité de matière nanostructurée conservée en fin de traitement et avant impact sur le substrat. Nous avons également exploité le modèle de dispersion de poudre pour la simuler la tache-dépôt formée par la projection ‘statique' de poudres de différents matériaux et de différentes granulométries.

Cette thèse se place dans un cadre de reconstruction urbaine et propose un corpus algorithmique pour extraire des formes simples sur les Modèles Numériques d'Elévation. Ce type de données décrit le relief d'une zone urbaine par une grille régulière de points à chacun desquels est associée une information de hauteur. Les modèles utilisés reposent sur l'utilisation de processus ponctuels marqués. Il s'agit de variables aléatoires dont les réalisations sont des configurations d'objets géométriques. Ces modèles permettent d'introduire des contraintes sur la forme des objets recherchés dans une image ainsi qu'un terme de régularisation modélisé par des interactions entre les objets. Une énergie peut être associée aux configurations d'objets et la configuration minimisant cette énergie trouvée au moyen d'un recuit-simulé couplé à un échantillonneur de type Monte Carlo par Chaîne de Markov à sauts réversibles (RJMCMC). Nous proposons quatre modèles pour extraire des caricatures de bâtiments à partir de descriptions altimétriques de zones urbaines denses. Chaque modèle est constitué par une forme d'objet, une énergie d'attache aux données et une énergie de régularisation. Les deux premiers modèles permettent d'extraire des formes simples (rectangles) en utilisant une contrainte d'homogénéité pour l'un et une détection des discontinuités pour l'autre. Le troisième modèle modélise les bâtiments par une forme polyhédrique. Le dernier modèle s'intéresse à l'apport d'une coopération entre des objets simples. Les algorithmes obtenus, automatiques, sont évalués sur des données réelles fournies par l'IGN (MNE Laser et optiques de différentes qualités).

Les villes regroupent la moitié de la population mondiale. Le suivi de leur croissance est devenu un enjeu majeur en aménagement du territoire, en écologie ou pour la défense. Ce suivi nécessite une connaissance fine de l'occupation du sol et de la géométrie tridimensionnelle du bâti. L'interférométrie par radar à synthèse d'ouverture (RSO) a fait ses preuves pour mesurer l'élévation du milieu naturel à basse résolution, mais il n'est pas certain que les progrès récents en imagerie RSO à haute résolution offrent des capacités suffisantes pour rendre compte de la complexité du milieu urbain. L'objectif de cette thèse était de répondre à cette question. En raison des réflexions multiples, des discontinuités et de la grande diversité des matériaux du milieu urbain, les images RSO et, a fortiori les interférogrammes, sont difficiles à interpréter. Une inversion classique de l'interférogramme ne fournit pas un Modèle Numérique de Surface exploitable. Nous avons donc proposé un ensemble de traitements de plus haut niveau qui prend en compte toute l'information disponible, aussi bien dans les images que dans leur interférogramme. Notre travail comporte une étape d'analyse du signal RSO et une chaîne de reconstruction 3D. L'analyse est faite à trois échelles: au niveau de la cellule de résolution par simulation électromagnétique, au niveau de l'objet par un découpage en sous-bandes et au niveau de la texture par modélisation statistique. Elle aboutit à l'identification de plusieurs caractéristiques de la scène et au calcul joint de l'élévation en chaque pixel et d'une classification.

Le but de ce travail est l'étude de processus à mésoéchelle de la mer Ligure. L'étude concerne les tourbillons et leur effet sur les variables biogéochimiques liées au plancton. Le modèle implémenté pour les simulations est ROMS, Regional Ocean Modeling System (Shchepetkin and McWilliams, 2005). La simulation à haute résolution du courant Ligure-Provençal-Catalan dans la Mer Ligure a révélé la présence d'instabilités qui induit la formation de structures tourbillonnaires. Ces tourbillons sont particulièrement intenses lorsque le forçage du vent provient d'un modèle météorologique à haute résolution spatiale et temporelle. Au contraire, l'utilisation de vents climatologiques à basse résolution comme forçage, induit des structures à mésoéchelle moins intenses et moins nombreuses. Les simulations forcées avec des vents réels permettent la formation d'intenses tourbillons anticycloniques entre le courant Ligure-Provençal et la côte. Ces tourbillons se déplacent lentement avec le courant. L'augmentation de la résolution permet de faire apparaître de petites structures tourbillonnaires autour du vortex principal. Les tourbillons observés sont associés à des structures de vitesse verticale complexes. La simulation montre d'importantes vitesses verticales positives (upwelling) près de la côte qui peuvent induire la remontée de masses d'eaux riches en nutriments, et des vitesses verticales négatives au centre du tourbillon. Les vitesses verticales augmentent avec la résolution de la simulation et peuvent atteindre des valeurs de l'ordre de la dizaine de mètres par jour. L'augmentation de la résolution permet de résoudre plus en détail la distribution du champ de vitesse et on observe ainsi des structures complexes de vitesses positives et négatives que s'alternent en correspondance d'un tourbillon. Pour ce qui concerne les vitesses horizontales, la résolution n'a pas un effet significatif sur leur intensité et distribution. Cela signifie que le processus est suffisamment bien résolu dans la simulation à plus basse résolution. Les vitesses verticales associées V au tourbillon peuvent avoir un rôle significatif pour la dynamique de l'écosystème marin de la Mer Ligure. Ces structures peuvent apporter en surface des masses d'eaux riches en nutriments qui peuvent mettre en route la chaîne trophique liée à la production primaire. Tous les vivants sont liés directement ou indirectement à cette production. Pour étudier l'impact des structures à mésoéchelle sur les variables biologiques nous avons appliqué un module biogéochimique couplé au modèle hydrodynamique. Les résultats mettent en évidence la complexité du comportement des variables biogéochimiques en présence des structures dynamiques, ainsi que leur variabilité en fonction du type de tourbillons (cyclonique ou anticyclonique). Les tourbillonnes cycloniques représentent essentiellement des zones de croissance de biomasse. Les tourbillonnes anticycloniques représentent aussi des zones d'accumulation de biomasse mais aussi des zones de croissance de biomasse due à la remontée d'eaux riches en nutriment au bord du tourbillon. Les résultats du travail de thèse et décrits dans ce document, montrent l'importance des processus à la mesoéchelle, l'impact de tourbillons sur la production biologique et la nécessité d'adopter une résolution spatiotemporelle suffisamment haute pour le modèle océanique et pour les forçages atmosphériques.