La prédiction de l'émission infrarouge des jets de propulseurs composites aluminisés nécessite principalement trois étapes : le calcul des grandeurs aérothermochimiques du jet, la conversion de ces grandeurs en propriétés optiques (coefficient d'absorption, coefficient de diffusion, fonction de phase) puis la résolution de l'équation de transfert radiatif. Cette thèse,essentiellement consacrée à cette troisième étape, propose de nouvelles voies pour l'application des modèles de bande aux cas de milieux biphasiques et diffusants.D'une part, nous avons étendu ces modèles aux cas de milieux caractérisés par un déséquilibre thermique entre gaz et particules. D'autre part, nous avons proposé une méthode de séparation de la luminance en deux contributions, appelées respectivement non diffusée et diffusée, à partir d'une idée originale de Liu et al. La contribution non diffusée est solution de l'équation de transfert radiatif obtenue en ignorant les effets de la diffusion. Par conséquent, elle peut être résolue par une formulation en modèles de bande. Cette approche permet de réduire les erreurs de corrélations spectrales au seul terme de luminance diffusée.Nous avons montré l'intérêt de ces approches par comparaison avec une résolution de l'équation de transfert radiatif en raie par raie, sur des milieux représentatifs de situations de télédétection de jets. / Prediction of infrared emission of exhaust plumes from aluminized composite rocket, follows mainly three steps : calculating aero-thermo-chemical values in the plume, converting those valuesto optical properties (absorption coefficient, scattering coefficient and phase function) and resolving the radiative transfer equation. This thesis is mostly devoted to this third step, and proposes new ways to use band models on two-phases and scattering media.Firstly, we extended band models to cases with thermic non equilibrium between gas and particles. Secondly, we proposed a method consisting in splitting radiance in two parts, one called un-scattered and the other scattered, from an original idea of Liu et al. The un-scattered part is solution of the radiative transfer equation obtained by ignoring scattering. As a result, the unscattered radiance can be found by using band models. By this approach, errors on spectral correlations are only present on the scattered radiance.We show the interest of thoses approches by comparing them with a line by line resolutionof the radiative transfer equation, on media representative of remote sensing cases of rocket exhaust plumes.

Le processus de formation d’étoiles se déroule selon plusieurs phases. Tout d’abord une phase à grande échelle, durant laquelle le nuage moléculaire se fragmente sous l’action de sa propre gravité et de la turbulence en coeurs denses gravitationnellement instables. Dans ces fragments le milieu est optiquement mince, l’énergie libérée par le travail de compression s’échappe sous forme de rayonnement, d’où un processus quasi isotherme. Lorsque le nuage devient optiquement épais à son propre rayonnement, la matière en effondrement forme un objet en équilibre hydrostatique appelé premier cœur dit de Larson. S’ensuit une phase d’accrétion, qui conduit ultimement à la dissociation du dihydrogène. Une partie du travail de compression est alors absorbée par l'énergie de dissociation de la molécule, et non plus convertie en énergie thermique, permettant à l'effondrement de recommencer. Lorsque que toutes les molécules de dihydrogene ont été dissociées, la phase adiabatique recommence et le second cœur de Larson (proto-étoile) est formé.L'ajout des éléments nécessaires au traitement de la magnétohydrodynamique (MHD) non-idéale dans le code à grille adaptative RAMSES constitue la première partie de la thèse. L'étude détaillée des stades ultimes (premier et second cœur de Larson) de la formation des étoiles constitue la seconde partie de la thèse. Cette étude a pu mettre en évidence des effets importants de la MHD non-idéale sur la répartition du champ magnétique et l'efficacité du transport de moment angulaire. / Stars formation occurs in several steps. First a large scale phase during which the molecular cloud undergo fragmentation due to its self-gravity and turbulence. In the gravitationally unstable fragments the medium is optically thin causing all the energy generated by the collapse to escape freely. This is called the isothermal compression phase. When the cloud becomes optically thick to its own radiation, an hydrostatic core forms: the first Larson core. Follow an adiabatic accretion phase ending up ultimately in the dissociation of dihydrogen molecules. Part of the energy from the gravitational collapse is absorbed by the chemical process allowing for another quasi isothermal collapse to start until depletion of dihydrogen molecules. When the adiabatic phase is restored, the second Larson core (proto-star) is formed.Coding the non-ideal magnetohydrodynamics (MHD) solver in the adaptive mesh refinement code RAMSES has been the focus for the first part of the thesis. The precise study of the last steps (first and second Larson core) of star formation is the second part of the thesis. This study highlighted the impact of non-ideal MHD on the magnetic field repartition and the efficiency of the angular momentum transport.

Les progrès effectués dans le domaine des instruments optiques permettent d’envisager désormais de placer des capteurs à moyenne résolution spatiale et haute résolution spectrale sur des orbites géostationnaires ce qui permettrait d’observer la Terre à la fois avec une moyenne résolution spatiale, une bonne résolution spectrale et une haute résolution temporelle, ce que ne permettent pas les satellites héliosynchrones. L'objectif de cette thèse est donc de simuler les images en luminance qui seraient acquises par un capteur géostationnaire hypothétique dédié à l'observation de la couleur de l'océan, afin de caractériser les images qui seraient acquises, de spécifier les futurs instruments et de valider des algorithmes de traitement sur des images réalistes. Pour parvenir à simuler les luminances spectrales qui seraient mesurées par le capteur, nous avons combiné les modèles de transfert radiatif directs de l'eau (Hydrolight) et de l'atmosphère (MODTRAN). Les données d’entrée d’Hydrolight ont été déduites des données provenant des cartes climatologiques de GlobColour (chlorophylle, matières minérales et organiques en suspension) pour simuler des images de couverture globale (Global Area Coverage) d'un capteur géostationnaire à basse résolution spatiale, moyenne résolution spectrale et haute résolution temporelle. L'analyse de ces images permet de déterminer la gamme de luminances pour chaque bande spectrale et pour différentes acquisitions au cours de la journée. Les résultats indiquent aussi les limites angulaires pour éviter le phénomène du sunglint et pour assurer un niveau minimum de luminance provenant de l'eau, contenu dans le signal mesuré. Les images en luminance ainsi simulées ont permis d’évaluer aussi la performance d’un algorithme d'estimation de concentration en chlorophylle de rapport de bandes dans les eaux au large. Les résultats montrent que l'erreur d'estimation est moyennement sensible aux angles d'observation et d’éclairement, importants dans les images géostationnaires. / Improvements of optical sensors now make it possible to launch medium spatial and high spectral resolution sensors on the Geostationary Earth Orbit (GEO). The GEO allows sensors to observe the Earth with high temporal resolution, unlike Low Earth Orbit (LEO) satellites. The objective of this thesis is to simulate radiance images that would be provided by a theoretical ocean color sensor, on a geostationary orbit, in order to characterize the images to be acquired, to specify the future instruments and to validate the image processing algorithms. For this purpose, we have combined a radiative transfer model for ocean waters (Hydrolight) and a radiative transfer model for the atmosphere (MODTRAN) to obtain simulated radiance images at the sea surface level and at the Top-Of-Atmosphere (TOA) level. The inputs of Hydrolight have been deduced from GlobColour Project data (chlorophyll, mineral suspended matter and organic matter) to simulate the Global Area Coverage (GAC) images of a geostationary sensor with low spatial resolution, medium spectral resolution and high temporal resolution. The analysis of these images enables to determine the radiance dynamic range for each spectral band over the day. The results indicate the angular limits to avoid the sunglint phenomenon and to ensure a minimum proportion of water leaving radiance in the measured TOA radiance. The simulated radiance images allow to evaluate the performance of an ocean color algorithm for chlorophyll concentration estimation in open waters. The results show that the error of chlorophyll estimation is moderately sensitive to observation and lighting angles, which are so important in geostationary images.

L'atténuation est un paramètre clé dans l'évaluation de l'aléa sismique car elle contrôle l'amplitude et la durée du mouvement du sol. Deux mécanismes contribuent à l'atténuation des ondes sismiques courte-période (f>1Hz, avec f la fréquence) : (1) l'absorption, quantifiée par son facteur de qualité Qi(f); et (2) le scattering, quantifié par Qsc(f). L'objectif principal de cette thèse est de cartographier l'atténuation sismique en France Métropolitaine en déterminant l'importance relative de ces deux processus. Pour mener à bien cette tâche, nous avons modélisé le transport de l'énergie sismique multi-diffusée - aussi appelée coda sismique - à l'aide de l'équation de transfert radiatif dans un milieu présentant des variations latérales des propriétés de diffusion et d'absorption. En utilisant une approche perturbative, nous avons calculé les noyaux de sensibilité de l'intensité de la coda sismique à des variations spatiales de Qi(f) et Qsc(f) à 2-D, dans le cas où la diffusion est isotrope. Dans un deuxième temps, nous avons mesuré la vitesse de décroissance énergétique de la coda sismique, quantifiée par le facteur de qualité de la coda Qc(f), sur plus de 120000 formes d'onde collectées en France, en Belgique et sur l'ensemble de l'arc alpin. L'approche théorique développée dans cette thèse nous a permis d'établir une relation approximative entre Qc et Qi. Cette relation prend la forme d'une intégrale de Qi sur le rai direct reliant la source à la station. Cette approximation a été utilisée afin de cartographier pour la première fois les variations régionales de l'absorption entre 1 et 32 Hz dans les Alpes et en France Métropolitaine. Nos résultats mettent en évidence de très fortes variations latérales (±30%) de l'absorption. À basse fréquence (f~1Hz), une corrélation claire apparaît entre la géologie de surface et les structures d'absorption : les zones de forte absorption se localisent sur les séries sédimentaires peu consolidées tandis que les régions de faible absorption correspondent au socle affleurant dans le Massif Armoricain, le Massif Central et dans les chaînes de montagne comme les Pyrénées et les Alpes. À haute fréquence (f~24Hz), la corrélation entre géologie de surface et atténuation disparaît. Nous avons formulé l'hypothèse que la dépendance fréquentielle de la structure en atténuation est due à un changement du contenu ondulatoire de la coda sismique avec la fréquence. Ainsi la sensibilité de la coda aux structures profondes de la croûte augmenterait avec la fréquence. Enfin, nous avons initié la mise en oeuvre des noyaux de sensibilité isotrope exacts à 2-D afin d'obtenir une tomographie d'absorption. L'inversion des mesures de Qc par décomposition en valeurs singulières nous a permis de construire la première tomographie d'absorption qui tienne compte de toute la sensibilité des trajets des ondes multi-diffusées. La carte préliminaire d'absorption obtenue pour les Pyrénées est prometteuse. Nos résultats fournissent des bases théoriques solides pour l'inversion linéarisée de Qsc et Qi à partir de la coda sismique. Ils offrent également des informations complémentaires aux tomographies de vitesse et des perspectives d'amélioration de la régionalisation des calculs de l'aléa sismique en France Métropolitaine. / Attenuation is a key parameter for seismic hazard assessment. It plays an important role in the observed variability of ground motion amplitude and duration. There are two main causes for attenuation of short period seismic waves (f>1Hz, with f the frequency) : (1) absorption, quantified by the quality factor Qi ; and (2) scattering, quantified by Qsc. The main objective of this thesis is to map the seismic attenuation in Metropolitan France and to determine the relative importance of these two processes. To model the transport of multiply-scattered seismic waves - also known as the seismic coda - we employ a scalar version of the radiative transfer equation with spatially dependent absorption and scattering properties. The sensitivity kernels of the coda intensity to spatial variations of Qi and Qsc are computed in 2-D isotropically and anisotropically scattering media. The coda quality factor Qc (f) - quantifying the decay rate of the seismic coda energy - have been estimated on a collection of 120000 waveforms recorded in France, Belgium and in the Alpine range. According to the theory developed in this thesis, we establish a linearized approximate relation between the coda quality factor Qc (f) and the absorption quality factor Qi. This relation is expressed as an integral along the direct ray path connecting the source to the station. This approximation is used to map regional variations of absorption in the Alps and in Metropolitan France between 1 and 32Hz. Our maps reveal strong lateral variations (±30%) of absorption in the crust. At low frequency (f~1Hz), the correlation between sedimentary deposits and high absorption regions is clear : strong absorption zones are localized on the poorly consolidated sedimentary series while low absorption regions correspond to the basement which outcrops in the Massif Armoricain, the Massif Central and in the mountain ranges as the Pyrenees or the Alps. At high frequency (f~24Hz), the correlation between surface geology and absorption structures tends to disappear. We hypothesize that the frequency dependence of the attenuation structure is caused by a change of the wavefield composition which accentuates the sensitivity of the coda to the deeper parts of the medium as the frequency increases. Finally, we initiate the implementation of the exact 2D isotropic sensitivity kernels to retrieve the crustal absorption structures. Using the linear relation between Qc and Qi, we obtain the first absorption map which takes into account the precise spatio-temporal sensitivity of coda waves. The inverse problem is solved with a singular value decomposition approach. The preliminary map of Qi for the Pyrenees is promising. Our results constitute a solid theoretical basis to develop linearized inversions of Qsc and Qi from the analysis of the seismic coda. They also significantly improve the knowledge of the regional variations of seismic attenuation in Metropolitan France. Our maps provide new insights on the crustal structure of the Alpine Range in complement to seismic velocity images. These attenuation maps have direct implication for the design of future seismic hazard maps.

La présente thèse, concernant l’ingénierie de la photosynthèse naturelle et notamment l’optimisation des procédés de production de microalgues, a été menée à l'Institut Pascal. L'approche suivie repose sur la construction de modèles de connaissance des photobioréacteurs, capables d'en prédire les performances quelle que soit leur géométrie, les conditions d'éclairement, ou le type de microalgues cultivées. Ces modèles de connaissance permettent de dégager des stratégies de conception et de conduite du procédé, qui sont utilisées pour développer et réaliser des démonstrateurs performants, à l'échelle pilote. Cette thèse a pour objet d'étude un photobioréacteur pilote de 24 litres utiles à agitation pneumatique dans lequel la lumière est apportée dans le volume réactionnel par 1000 fibres optiques à diffusion latérale. Sa conception repose sur le principe innovant de Dilution Contrôlée du Flux solaire en Volume (DiCoFluV) dont l'objectif est d'atteindre l'efficacité thermodynamique maximale de la photosynthèse. Au cours de cette thèse, le réacteur a été étudié en conditions d’éclairement parfaitement contrôlées, grâce à des lampes à décharges. Il a été rendu fonctionnel (en incluant sa régulation et son amélioration), puis caractérisé en termes d'hydrodynamique et de transfert radiatif. Notamment, le flux radiatif incident à la surface des fibres optiques a été déterminé par des expériences d'actinométrie (sel de Reinecke), grâce à un traitement original permettant l'analyse de situations avec absorption partielle du rayonnement et en géométrie complexe. Ensuite, une année de culture ininterrompue de la cyanobactérie Arthrospira platensis (dont six phases de fonctionnement continu) ont permis de mesurer la vitesse volumétrique moyenne de production de biomasse et l'efficacité thermodynamique du pilote. Pour les faibles densités de flux testées (dilution), les résultats expérimentaux ont montré l’apparition d’un phénomène de photorespiration qui a été intégré dans le modèle de couplage thermocinétique. De plus, ce manuscrit présente également l'utilisation de la méthode de Monte Carlo pour la résolution des modèles radiatifs dans la géométrie réelle du procédé, définie à partir d’une conception assistée par ordinateur, ce qui est une nouvelle avancée pour le traitement des géométries complexes. Après validation de l'adéquation entre les mesures expérimentales et l'estimation prédictive par le modèle, des premières simulations du pilote en fonctionnement solaire ont été menées, à partir de bases de données solaires (DNI). Les résultats obtenus donnent des premières indications quant aux paramètres d'ingénierie (en particulier le facteur de dilution) menant à une productivité surfacique moyenne annuelle maximale, en fonction de l'implantation géographique. / The present PhD dissertation deals with photobioreaction engineering for efficient microalgae production. The approach is based on the construction of knowledge models that permit predicting the performances of the process, whatever its design, the illumination conditions, or the microalgae species cultivated. These models are used to establish optimal design and control strategies that are implemented to construct and operate pilot-scale plants. Here, a 24 liters air-lift photobioreactor is studied, that is based on the principle of incident solar light-flux dilution for approaching the maximum thermodynamic efficiency of natural photosynthesis. For that purpose, the culture volume is internally illuminated by 1000 light-diffusing optical fibers. As a first step toward solar production, this PhD work focuses on perfectly controlled illumination conditions ensured by discharge lamps. First, the reactor hydrodynamics and radiative transfer are characterized. In particular, incident light-flux density at the fiber surfaces is measured by actinometry (with Reinecke salt), thanks to a novel treatment enabling analyses of situations with partial light absorption and complex geometries. Then, the mean volumetric rate of biomass production and the process efficiency are measured based on one year of continuous Arthrospira platensis culture. For the low radiative flux densities tested (dilution), photorespiration by the cyanobacterium is observed and included in the thermokinetic model. Moreover, this dissertation includes a presentation of the Monte Carlo method for solving the radiative transfer equation withinthe complex geometry of the computer aided design used for manufacturing the reactor. After validation against experimental results, the model is predictively used to simulate the pilot operating with natural solar light, based on solar DNI databases. These results indicate engineering parameters (in particular the dilution factor) for optimal yearly-averaged surface productivity, as a function of Earth location.

La caractérisation du couvert nival en forêt boréale est un élément important pour la compréhension des régimes climatiques et hydrologiques. Depuis plusieurs années, l’utilisation des micro-ondes passives est étudiée pour l’estimation de l’équivalent en eau de la neige (SWE : Snow Water Equivalent) à partir de capteurs satellitaires. Les algorithmes empiriques traditionnels étant limités en forêt boréale, le couplage d’un modèle de transfert radiatif (MTR) micro-onde passive (qui prend en compte les contributions du sol, de la neige, de la végétation et de l’atmosphère) avec un modèle de neige pour l’inversion du SWE semble une avenue prometteuse. La thèse vise donc à coupler un MTR avec le schéma de surface du modèle climatique canadien (CLASS) dans une perspective d’application opérationnelle pour les estimations de SWE à partir de données satellitaires micro-onde à 10.7, 19 et 37 GHz. Dans ce contexte, certains aspects centraux du MTR, dont l’effet de la taille des grains ainsi que la contribution de la végétation sont développés et quantifiés. Le premier aspect étudié dans la thèse concerne l’adaptation du modèle d’émission micro-onde passive DMRT-ML (Dense media radiative transfer theory – multi layer) pour l’intégration d’une nouvelle métrique représentant la taille des grains (surface spécifique des grains de neige: SSA). L’étude basée sur des mesures radiométriques et de neige in situ, montre la pertinence de l’utilisation de la SSA dans DMRT-ML et permet d’analyser le sens physique de l’adaptation nécessaire pour amener le modèle à simuler les températures de brillance (T[indice inférieur B) de la neige avec une erreur quadratique moyenne minimale de l’ordre de 13 K. Dans un contexte du couplage entre le modèle de neige de CLASS et DMRT-ML, un modèle d’évolution de la SSA est ensuite implémenté dans CLASS. Les SSA simulées par le module développé sont validées avec des données in situ basées sur la réflectance de la neige dans l’infrarouge à courte longueur d’onde pour différents types d’environnement. Au niveau de la contribution de la végétation, le modèle γ-ω a été étudié à partir de différentes bases de données (satellite, avion et au sol) en forêt boréale dense. L’étude montre l’importance de la considération de la diffusion (ω) pour l’estimation de l’émission de la végétation, paramètre auparavant généralement négligé aux hautes fréquences. Ensuite, des relations entre les transmissivités et certains paramètres structuraux de la forêt, dont l’indice de surface foliaire (LAI), ont été établies pour des forêts boréales en été. Des valeurs d’albédo de diffusion (ω) ainsi que les paramètres définissant la réflectivité du sol (QH) en forêt boréale ont aussi été inversées. Finalement, les simulations de T [indice inférieur] B issues du couplage du MTR (DMRT-ML, modèle γ-ω, et modèle atmosphérique) avec CLASS (dont les SSA simulées) ont été comparées avec les données AMSR-E sur une série temporelle continue de sept ans. Les premières comparaisons montrent une différence entre les paramètres de végétation (γ-ω) d’été et d’hiver, ainsi qu’une importante contribution des croûtes de glace dans la neige au signal. Les simulations du modèle ajusté montrent une bonne correspondance avec les observations d’AMSR-E (de l’ordre de 3 à 7 K selon la fréquence et la polarisation). Des tests de sensibilité montrent par contre une faible sensibilité du MTR/CLASS au SWE pour des forêts denses et des couverts nivaux épais. Le MTR-CLASS développé pourrait permettre l’assimilation de températures de brillance satellitaires en forêt boréale dans des systèmes opérationnels pour l’amélioration de paramètres de surface, dont la neige, dans les modèles météorologiques et climatiques.

The feasibility of a radiance assimilation using a multi-layered snow physical model to estimate snow physical parameters is studied.The work is divided in five parts.The first two chapters are dedicated to the literature review. In the third chapter, experimental work was conducted in the alpine snow to estimate snow correlation (for microwave emission modelling) using near-infrared digital photography. We made microwave radiometric and near-infrared reflectance measurements of snow slabs under different experimental conditions. We used an empirical relation to link near-infrared reflectance of snow to the specific surface area (SSA), and converted the SSA into the correlation length. From the measurements of snow radiances at 21 and 35 GHz, we derived the microwave scattering coefficient by inverting two coupled radiative transfer models (RTM) (the sandwich and six-flux model).The correlation lengths found are in the same range as those determined in the literature using cold laboratory work.The technique shows great potential in the determination of the snow correlation length under field conditions. In the fourth chapter, the performance of the ensemble Kalman filter (EnKF) for snow water equivalent (SWE) estimation is assessed by assimilating synthetic microwave observations at Ground Based Microwave Radiometer (GBMR-7) frequencies (18.7, 23.8, 36.5, 89 vertical and horizontal polarization) into a snow physics model, CROCUS. CROCUS has a realistic stratigraphic and ice layer modelling scheme. This work builds on previous methods that used snow physics model with limited number of layers. Data assimilation methods require accurate predictions of the brightness temperature (Tb) emitted by the snowpack. It has been shown that the accuracy of RTMs is sensitive to the stratigraphic representation of the snowpack. However, as the stratigraphic fidelity increases, the number of layers increases, as does the number of state variables estimated in the assimilation. One goal of the present study is to investigate whether passive microwave measurements can be used in a radiance assimilation (RA) scheme to characterize a more realistic stratigraphy.The EnKF run was performed with an ensemble size of 20 using artificially biased meteorological forcing data.The snow model was given biased precipitation to represent systematic errors introduced in modelling, yet the EnKF was still able to recover the"true" value of SWE with a seasonally-integrated RMSE of only 1.2 cm (8.1%).The RA was also able to extract the grain size profile at much higher dimensionality which shows that the many-to-one problem of SWE-Tb relationship can be overcome by assimilation, even when the grain size profile varies constantly with depth.The last chapter was on the validation of the data assimilation system using a point-scale radiance observations from the CLPX-1 GBMR-7. We first predicted snow radiance by coupling the snow model CROCUS to the snow emission model (MEMLS). Significant improvement of Tb simulation was achieved for the late February window for all three frequencies.The range of the underestimation of the polarization difference is between 25% and 75%. We then assimilated all six channels measurements of the GBMR-7.The filter was able to accurately retrieve the SWE for periods of time when the Tb measurements were available.The results show that RA using EnKF with a multi-layered snow model can be used to determine snow physical parameters even with a biased precipitation forcing.

This thesis takes its place in a context where information on the biophysical state of forest ecosystems at spatial scales only remote sensing can retrieve is in demand more than ever. In order to provide reliable information using validated approaches, the remote sensing research community recognises the need for new and innovative methods, especially in heterogeneous environments like savannas. The recent emergence of terrestrial LiDAR scanners (TLS) and the increase in the computational capability of computers which allow running ray tracing model simulations with a high level of realism hold great potential to improve our understanding of the processes influencing the radiance measured by satellite sensors. This thesis makes use of these two cutting edge technologies for estimating the spatial distribution of tree leaf area, a key element of modeling radiative transfer processes. The first part of the thesis concerns the development of methods for estimating tridimensional leaf area distribution in a savanna environment from TLS measurements. The methods presented address certain issues related to TLs measures affecting the application of classical theories (the probability of light transmission and the contact frequency) to the estimation of leaf area through indirect means. These issues pertain to the cross-section of laser pulses emitted by a TLS and the occlusion effects caused by the interception of laser pulses by material inside the crown. The developed methods also exploit additional information provided by the active nature of the TLS sensor that is not available to passive sensors like hemispherical photography, i.e. the intensity of a pulse return offers the possibility to distinguish between energy interception by wood and foliage. A simplified approach of this method is presented to promote its use by other research groups. This approach consists of a series of parameterisations and represents a significant gain in terms of the required resources to produce the leaf area, estimates. The second part of the thesis explores the combination of the tree representations generated in the first part with a ray tracing model to simulate the interactions of light with tree crowns. This approach is highly innovative and our study showed its potential to improve our understanding of the factors influencing the radiative environment in a savanna. The methods presented offer a solution to map leaf area at the individual tree scale over large areas from very high spatial resolution imagery.

La réalisation de nouveaux instruments de télédétection à très haute résolution spatiale offre la possibilité d'étudier plus précisément les villes. Pour ces études, la connaissance de l'atmosphère et plus particulièrement des aérosols peut s'avérer essentielle. Le but de cette thèse est donc de développer une méthode de caractérisation des aérosols adaptée aux images de télédétection des milieux urbains à l'échelle métrique dans les domaines visible et proche-infrarouge. Dans un premier temps, les propriétés optiques de ces particules ont été étudiées en utilisant les données fournies par 68 stations urbaines du réseau AERONET. Ensuite, afin de pouvoir évaluer l'impact des aérosols présents dans les villes sur le signal, un code de transfert radiatif 3D a été réalisé : AMARTIS v2. L'utilisation de cet outil pour une scène urbaine typique a permis de quantifier l'impact des particules sur le signal, à l'ombre et au soleil, en fonction de leurs propriétés optiques. Enfin, une méthode de télédétection des aérosols a été définie, basée sur l'observation de transitions ombre/soleil. Afin de mettre en oeuvre cette méthode, un code d'inversion a été développé : OSIS. Une étude de sensibilité d'OSIS a alors été menée à partir d'images synthétiques générées avec AMARTIS v2 et une utilisation expérimentale a été effectuée sur des acquisitions PELICAN obtenues lors de la campagne aéroportée MUSARDE sur la ville de Toulouse. Ces études ont notamment permis de quantifier la précision intrinsèque d'OSIS, comparable aux précisions obtenues avec les produits satellitaires pour l'inversion des épaisseurs optiques, et de montrer que cette procédure est applicable à tout instrument à très haute résolution spatiale, multispectral ou hyperspectral, aéroporté ou satellitaire.

L'objectif de ce travail est de développer un algorithme utilisant des réseaux neuronaux pour estimer des variables biophysiques des couverts végétaux à partir des données du capteur VEGETATION : la fraction de trou dans différentes directions de visée, l'indice foliaire, la fraction de rayonnement photosynthétiquement actif absorbé par le couvert, l'angle moyen d'inclinaison des feuilles, et la teneur en chlorophylle intégrée au niveau du couvert. Une base de données a été élaborée à partir de simulations de modèles de transfert radiatif dans la végétation. Ceci nous a permis de développer des réseaux de neurones pour estimer les différentes variables biophysiques considérées. L'estimation se fait à partir des réflectances bidirectionnelles mesurées par VEGETATION durant un cycle orbital, associées à l'angle zénithal solaire lors de l'acquisition de ces données. L'analyse de sensibilité à l'hétérogénéité intrapixellaire montre que les variables de type fraction de trou et fAPAR sont beaucoup moins sensibles que les variables de type LAI. La validation de l'algorithme sur des données expérimentales montre la robustesse de la méthode lorsque le nombre de données bidirectionnelles acquises durant le cycle orbital est faible. Par ailleurs, elle confirme que les performances d'estimation des variables biophysiques sont affectées non seulement par la représentativité du couvert et des conditions d'observation qui lui sont associées dans la base d'apprentissage, mais aussi par les hypothèses inhérentes aux modèles de transfert radiatif utilisés pour générer cette base. Nous mettons également en évidence la sensibilité de l'algorithme à des bruits de natures diverses (mesure, correction atmosphérique, ...). Les résultats obtenus sur les données expérimentales sont satisfaisants et aboutissent à un écart quadratique moyen inférieur à 0,1. L'intérêt des réseaux de neurones apparaît donc clairement, en particulier vis à vis des approches basées sur les indices de végétation.