Academic literature on the topic 'Sous-sol (pédologie) – Teneur en eau'

A partir de résultats expérimentaux, des équations d'infiltration et de continuité, nous avons établi une équation permettant de prédire la profondeur maximale humidifiée par un goutteur à la surface du sol. Connaissant l'humidité initiale du sol Øi, la conductivité hydraulique Kf et la teneur en eau au niveau du front Øi ; l'équation établie permet de calculer la profondeur du front d'humectation Zf(t) à partir du suivi du rayon du front d'humectation Rf(t) à la surface du sol. Pour tester la validité de cette nouvelle approche, des essais ont été menés sur deux types de sol (sableux et limono-argileux) et avec trois débits différents. Les valeurs de Zf(t) mesurées (par deux techniques différentes) ont été comparées à celles calculées par cette équation, les résultats étaient, dans tous les cas, quasi-identiques (coefficient de corrélation R2, dans tous les cas, supérieur ou égal à 0.94).

Depuis plusieurs années, la Zone Non Saturée d'une parcelle expérimentale au nord de l'Italie (plaine du Pô) a fait l'objet de nombreuses études sur les caractéristiques physiques et géologiques du milieu, ainsi que sur les modalités du mouvement des eaux d'infiltration. Le site a été équipé de 6 tensiomètres à mercure, 6 bougies poreuses pour l'échantillonnage de l'eau et d'un tube d'accès pour la sonde à neutrons. L'application des différentes méthodes d'évaluation a fourni, pour la période d'étude (Septembre 1994-Août 1995), des valeurs de l'infiltration sensiblement différentes. La méthode physique, basée sur le calcul du bilan hydrique annuel, et considérant le sol comme un unique réservoir monocouche, a indiqué que 19% du total des précipitations arrive à s'infiltrer. D'autre part, un modèle de simulation intégrant les caractéristiques hydrodynamiques du sol a permis d'évaluer la recharge à 29% des précipitations annuelles. Le calcul de la répartition de la teneur en eau sur le profil porte sur la description des fluctuations temporelles du plan de flux nul, en relation avec les épisodes pluvieux et l'intensité de l'évaporation. La méthode géochimique du "profil en ions conservés" entre l'eau de pluie et du sol a permis, quant à elle, d'évaluer la lame d'eau infiltrée à environ 12% des précipitations annuelles. La composition isotopique des eaux extraites du sol manifeste un enrichissement global en isotopes lourds (2H et 18O) par rapport aux eaux météoriques, conséquence d'une évaporation intense. L'infiltration semble être de type "piston flow" dispersif. Elle suit, seulement pendant l'hiver, un flux où la composante convective dans la microporosité est prédominante sur la composante diffusive qui a lieu, de préférence, dans la macroporosité.

La matière organique (MO) est un élément important de la fertilité des sols cultivés. Les mesures de réflectance ont déjàété utilisées avec succès pour prédire la teneur en carbone organique (CO) des sols agricoles, soit au laboratoire surdes échantillons de sols séchés, soit directement au champ à l'aide de spectro-radiomètres de terrain ou grâce à desimages satellitaires. Cette étude propose une méthode de prédiction spatialisée du CO des sols, à partir d'imageriedrone. Une parcelle agricole de 13 ha « la Grande Borne », située dans la plaine de Versailles (78), a été suivie en avril2013, peu avant le semis alors qu'elle était encore en sol nu. Le drone AIRINOV® utilisé était équipé de la caméramultispectrale quadri-bande multiSPEC 4C® (550nm, 660nm, 735 nm et 790 nm) et évoluait à une altitude de 150 m.Vingt-trois placettes de contrôle au sol, réparties dans la parcelle, ont été délimitées par des cibles géolocaliséescentimétriquement par DGPS et ont fait l'objet de mesures synchrones au vol drone : mesures spectrales de terrain(ASD FieldSpec 3®), mesures de rugosité par une méthode photogrammétrique, prélèvements d'échantillons de solpour analyse de la composition chimique et de la teneur en eau. Les images obtenues d'une résolution de 15cm ont étégéoréférencées avec une précision infra-pixellaire. Les spectres de terrain ont été agrégés sur les bandes spectrales de2la multiSPEC® 4C afin de valider par régression linéaire la qualité des données de réflectance du drone. Les Robtenus, entre 0,8 et 0,9 selon les bandes, ainsi que l'erreur comprise entre 1 et 3% en réflectance, ont montrés lafiabilité du capteur multiSPEC 4C®. La teneur en CO a ensuite été estimée par régression des moindres carrés partiels(PLSR) à partir des mesures spectrales sur les 23 placettes. L'erreur quadratique moyenne de validation croisée(RMSECV) par la méthode LOO (Leave One Out) était de 1,97 g de CO par kg de sol. Un second modèle intégrant lacorrection des effets de l'humidité et de la rugosité sur la réflectance, a permis d'améliorer la qualité de la prédiction de18% avec une RMSECV de 1,61 g/kg. Le modèle a enfin été spatialisé sous SIG et comparé à de simples interpolationsspatiales des analyses de sol.

Cet article vise à démontrer l’intérêt scientifique et opérationnel de coupler de façon systématique carottages et méthodes géophysiques pour l’étude et la gestion de sédiments fluviaux. Deux techniques géophysiques sont présentées, car particulièrement adaptées au contexte fluvial : le radar géologique ou GPR (ground penetrating radar) est fondé sur l’émission d’ondes électromagnétiques et permet d’obtenir une image des structures sédimentaires superficielles à terre ou en eau peu profonde, tandis que le sondeur de sédiment (ou sub bottom profiler) est une méthode acoustique qui s’utilise en contexte aquatique uniquement, généralement dans des profondeurs d’eau plus importantes que le GPR. La mise en œuvre des méthodes géophysiques demande peu d’efforts logistiques et donne des informations sur la géométrie du sous-sol avec une très bonne résolution, mais ne permet pas une caractérisation des sédiments. Les carottes sédimentaires – prélevées à terre ou en eau – peuvent être échantillonnées et analysées, ce qui permet une caractérisation des sédiments (granulométrie, chimie, teneur en matière organique, datation, etc.) ponctuelle dans l’espace. Les deux approches sont donc complémentaires et leur couplage permet d’étendre la caractérisation des carottes à l’ensemble des structures identifiées par les sondages géophysiques. Pour illustrer ce propos, la démarche est appliquée à un ouvrage hydraulique, où une combinaison de cinq profils GPR et trois carottes sédimentaires suffit pour identifier et représenter en trois dimensions l’ensemble des structures sédimentaires, déterminer la chronologie de remplissage de l’aménagement et calculer et localiser les volumes de sédiments plus ou moins pollués en polychlorobiphényles indicateurs (PCBi).

Dans les périmètres irrigués du Maroc, la pollution nitrique diffuse des eaux souterraines augmente le risque de détérioration de la qualité des ressources hydriques et engendre un risque sanitaire pour la population rurale s'approvisionnant le plus souvent directement de l'aquifère. La présente étude, réalisée dans le périmètre du Tadla (3600 km2), située au centre du Maroc, vise à faire l'état des lieux de la contamination des eaux souterraines par les nitrates et à identifier les causes potentielles de cette pollution. Le suivi de la concentration des eaux souterraines en nitrates a été réalisé durant la période août 1996 - avril 1998 dans un réseau de cent puits. Les paramètres du milieu physique et des pratiques culturales ont été déterminés soit par des mesures au champs, des analyses au laboratoire ou par enquête auprès des agriculteurs dans un réseau de 40 sites. Les résultats sont présentés sous forme de cartes thématiques par le SIG. Le niveau actuel de la pollution nitrique des eaux souterraines dans le périmètre irrigué du Tadla commence à devenir critique. De part la pollution ponctuelle mise en évidence aux alentours de quelques agglomérations urbaines et à l'aval des sucreries, les zones dont les teneurs en nitrates des eaux souterraines dépassent 50 mg NO3- /l sont localisées en aval des Béni Amir et en aval et à l'Est des Béni Moussa Est. Béni Moussa Ouest reste en général moins polluée que Béni Moussa Est. Les apports excessifs en eau et en engrais azotés pour les cultures pratiquées (blé, betterave, cultures maraîchères, arboriculture et luzerne), sont les causes directes de la pollution diffuse. Par ailleurs, il est certain que la texture du sol (argileuse, argilo-limoneuse et limono-argileuse), la teneur en argile (29% à 55%) et tous les paramètres qui leurs sont corrélés notamment la réserve utile du sol, réduisent le risque de lixiviation des nitrates en profondeur et assurent en conséquence un contrôle naturel de la pollution nitrique des eaux souterraines. Néanmoins, tous les autres facteurs étudiés notamment le niveau des nappes faible par rapport au sol (<4m-12m), la vitesse d'infiltration du sol en milieu saturé généralement élevée (20-450 mm/h), les teneurs assez importantes du sol en matière organique (1.5%-3.5%) et en azote total (0.10%-0.15%), joueraient positivement sur la vulnérabilité de la nappe à cette pollution. De même, le pouvoir de minéralisation élevé des sols couplé aux conditions écologiques favorables à une minéralisation intense de la matière organique native du sol, assez caractéristique des régions irriguées à climat méditerranéen, seraient en partie responsable de cette pollution. Néanmoins, Aucune corrélation n'a été mise en évidence entre la teneur en nitrate des eaux souterraines et les combinaisons linéaires des variables du milieu physique. Ceci serait dû à l'homogénéité du milieu physique et à l'accumulation de la pollution suivant le sens d'écoulement des nappes de l'amont vers l'aval hydraulique. Vu le niveau d'intensification agricole que connaît cette région, la situation risquerait de s'aggraver rapidement induisant, en dehors de toute mesure de prévention, la détérioration de la qualité des eaux souterraines dans cette région.

Le teck (Tectona grandis L. f.) est une espèce ligneuse tropicale dont le bois est très demandé. Au Bénin, la provenance tanzanienne et la provenance locale sont les deux variétés de cette espèce les plus plantées. Pour contribuer à la caractérisa- tion du bois de teck du Bénin, nous avons déterminé et comparé, pour le bois de chaque provenance, les modules d’élasti- cité et de cisaillement, la longueur du fût, la circonférence et l’épaisseur d’écorce à hauteur d’abattage, le pourcentage de bois de cœur, la circonférence à une hau- teur de 1,30 m du sol et la masse volu- mique (ou la densité) à 12 % de teneur en eau. Les échantillons ont été préle- vés dans les plantations de la Lama sur des arbres de 20 à 25 ans (six billes de provenance locale et quatre billes de provenance tanzanienne). Les résultats obtenus sur la base de cet échantillon indiquent que, pour cette classe d’âge, les valeurs de ces propriétés pour le bois de teck de provenance locale sont supé- rieures à celles de la provenance tanza- nienne, sauf pour le pourcentage de bois de cœur. Sur la base de notre échantil- lonnage, le teck de provenance locale semble donc d’une certaine qualité, qu’il faut tendre à valoriser par la sélection bio- logique ou clonale des meilleurs sujets, aptes à s’adapter non seulement au chan- gement climatique mais aussi aux struc- tures des sols. Cela permettrait de favori- ser la production du bois de teck sur toute l’étendue du territoire national, voire de la sous-région. Toutefois, la mise en place d’une variété hybride conciliant les atouts des deux provenances pourrait aussi être une grande avancée dans l’amélioration des essences exploitées au Bénin.

Une grande quantité d'eau est perdue dans les zones de pâturage et prairies naturelles du fait de la présence dans ces régions de plantes à forte transpiration. La gestion du couvert végétal et des bassins versants a été proposée comme moyen pour augmenter la disponibilité des ressources en eau. Des efforts accrûs ont été consacrés au développement de pratiques de gestion et d'outils pour évaluer le potentiel d'augmentation de la ressource en eau. La modélisation hydrologique joue un rôle clé dans ces efforts. Un des outils les plus complets pour la modélisation dans les zones de pâturage et prairies naturelles est le modèle SPUR. Il s'agit d'un modèle de bassin versant, spatialement semi-distribué. Le modèle est constitué de cinq modules principaux qui incluent les aspects suivants: climat, hydrologie, plantes, animaux et économie. La composante hydrologique du modèle prend en compte à un pas de temps journalier les phénomènes de ruissellement, évapotranspiration, percolation et écoulement latéral. Le ruissellement est calculé à partir du numéro de courbe qui dépend du couvert végétal, des pratiques culturales, ainsi que des conditions hydrologiques. Cependant l'utilisation par le modèle de la méthode du numéro de courbe pour déterminer le ruissellement pose de sérieux problèmes quant à l'efficacité du modèle. Dans notre recherche, nous avons substitué la méthode des numéros de courbe par l'équation de Green et Ampt. Un avantage majeur de cette approche est l'utilisation de l'intensité de la pluie comme variable de forçage au lieu de la pluie journalière. De plus, cette équation d'infiltration utilise des paramètres physiques comme la conductivité hydraulique à saturation. L'objectif de cet article est de présenter les performances du modèle SPUR original et modifié sur trois types de couvert végétal : sol nu, sol enherbé et buissons. Trois années de mesures collectées sur le bassin versant de Throckmorton (Texas, Etats Unis d'Amérique) ont été utilisées pour la calibration et la validation des modèles. Les performances des modèles ont été évaluées en utilisant le coefficient d'efficience de Nash et Sutcliff. Le calage a porté sur la première année de mesure. Pour le modèle original, le calage a consisté à ajuster les numéros de courbe de manière à optimiser l'efficience. Pour le modèle modifié, il n'a été procédé à aucun calage. Les valeurs de conductivité hydraulique à saturation ont été estimées en utilisant des équations de pédotransfert en se basant sur les propriétés texturales et structurales des sols. L'introduction dans le modèle SPUR de l'équation de Green et Ampt a considérablement amélioré la performance du modèle pour la prévision du ruissellement sur tous les sites. L'efficience moyenne pour les prévisions du ruissellement mensuel sur sol nu est de 0.16, alors que celle ci est négative pour le modèle original (-0.11). Pour les sites enherbés l'efficience du modèle modifié est de 0.48 alors qu'elle est négative pour le modèle original. L'utilisation du numéro de courbe a résulté en une surestimation systématique du ruissellement sur tous les sites. De manière générale, le modèle original et le modèle modifié présentent de meilleures performances sur les sites non nus. Ceci est dû au fait que les deux modèles sous-estiment de manière significative l'évaporation sur les sols nus. Un des désavantages des deux modèles est en effet de limiter l'évaporation aux premiers 15 cm du sol. L'introduction de l'équation de Green et Ampt a amélioré les performances du modèle pour la prévision du ruissellement aussi bien à l'échelle mensuelle qu'annuelle. De plus, le modèle modifié est sensible au type d'occupation du sol et est donc adapté comme outil pour l'analyse de scénario en vue de préserver les ressources en eau. Une analyse de sensibilité a été conduite afin d'évaluer l'impact des paramètres d'entrée sur les sorties des deux modèles. L'analyse de sensibilité a consisté à modifier systématiquement les paramètres d'entrée de plus ou moins 10%. Pour le modèle original, l'analyse a porté sur l'influence du numéro de courbe, et pour le modèle modifié celle ci a porté sur l'étude de l'impact lié aux paramètres utilisés pour calculer la conductivité hydraulique à saturation. Concernant le modèle original, une augmentation du numéro de courbe de 10% entraîne une augmentation du ruissellement de 120% pour le sol nu, et aux alentours de 100% pour les autres sites. L'impact de ces variations sur l'évaporanspiration est minimal, avec une variation maximale de 16% pour le sol nu. Concernant le modèle modifié, la teneur du sol en sable est le paramètre ayant la plus grande influence sur la quantité d'eau ruisselée pour le sol nu. Par contre, pour les lysimètres ayant un couvert végétal, le pourcentage de sol couvert par la canopée est le factor majeur contrôlant la quantité d'eau ruisselée. Les paramètres liés au couvert végétal ont un plus grand impact sur le ruissellement que les paramètres liés aux propriétés intrinsèques du sol. Globalement l'introduction de l'équation de Green et Ampt a amélioré les capacités prédictives du modèle. Outre le fait que le modèle modifié ne nécessite pas un calage particulier pour la détermination des paramètres de transfert de l'eau dans le sol, il se base sur l'intensité de la pluie pour la détermination du ruissellement. Il a été montré que le modèle modifié est sensible aux changements de type d'occupation du sol. Il peut donc être donc utilisé comme outil pour évaluer l'impact de différents scénarios d'occupation du sol sur les ressources en eau dans les zones de pâturage et prairies naturelles. Toutefois, des améliorations, telles que l'introduction de l'impact du développement de fissures sur l'infiltration (écoulement préférentiel) ainsi que sur le phénomène d'évaporation devraient être prises en compte afin d'améliorer les prévisions du bilan hydrologiques, notamment sur sol nu.

Le radar géologique est utilisé pour déterminer différentes caractéristiques du sous-sol. Cependant, l'analyse des signaux est parfois compliquée par la structure hétérogène du milieu. De plus, la méthode est très sensible à la présence d'eau dans le sous-sol. Afin de pallier ces problèmes, nous utilisons un programme de modélisation en domaine temporel par différences finies (FDTD). Nous ajustons l'impulsion source pour nous rapprocher le plus possible des résultats réels. Nous commençons par modéliser la structure d'un bassin à flanc de falaise dans les Alpes et y découvrons une structure géologique qui peut indiquer la présence de mouvements d'enfoncement du bassin. Une étude sur des cavités dans des calcaires nous permet ensuite de discriminer les signaux dus à des cavités remplies d'argile et des cavités vides. Enfin, une dernière étude nous permet de suivre l'évolution d'une couche de neige au cours du temps. Dans un bac, un sable dans lequel sont enterrés des objets est soumis à des variations de hauteur de nappe phréatique. Nous notons qu'aucune réflexion sur le toit de la nappe phréatique n'est enregistré sur les données réelles. Nous montrons que les 3 formules utilisées fréquemment en radar pour relier la permittivité diélectrique et la teneur en eau du milieu donnent des résultats comparables, tant sur les valeurs de teneur en eau que sur les variations de volume d'eau. Une modélisation par FDTD nous permet de donner une estimation sur les paramètres diélectriques du sable. Nous simulons dans le bac une pollution au gasoil à partir de la surface. Nous montrons que le gasoil ne se répartit pas de façon homogène dans le bac. Peu de signes permettent de caractériser la pollution en-dehors de la variation de la vitesse de propagation des ondes dans le milieu. Une modélisation permet de donner une forme au panache de pollution. Nous montrons également que des déplacements d'eau importants ont lieu lors de la mise en place du gasoil
Ground Penetrating Radar is used in near surface acquisitions to determine several characteristics from the soil. The analysis of the data is generally complicated by the heterogeneity of the near surface layer. The influence of water is also very strong and a change is the water content of the soil can change the results. In order to improve the data analysis, we use a modelling program with a finite differences in time domain (FDTD) algorithm. We determine the shape of the source impulse to be as near as possible from the real results. At first, we use this program to model the structure in a peat basin in the French Alps. We discover a geological structure which may result from a progressive sinking of the basin nearby the cliff. Another study helps us to discriminate between karsts filled with clay and hollow karsts in limestone. Finally, we can follow the changes affecting a snow layer at several times during the winter. We mounted a laboratory experiment in order to reproduce a groundwater level in a sand layer at varying depths. We had no signal reflection on the top of the saturated sand layer. We show that the 3 commonly used relations linking the water content to the relative dielectric permittivity give similar results, both on the water contents and on the variations of the water volume in the sand. With the modelling of the profiles from the sans box, we can give some estimations for the dielectric parameters of the sand. At last, we simulated an accidental fuel pollution from the surface of the sand box. We show that the fuel is not homogeneously distributed in the sand. We have no reflection on the fuel saturated layer, and only the wave velocity is clearly affected by the presence of fuel. The plume is more concentrated in the vicinity of the injection point and is the origin of water movements in the sand. New signals appear after a longer time due to these changes. A modelling gives a good idea from the shape of the plume

Cette recherche doctorale porte sur la modélisation 3D de la dynamique des teneurs en eau dans le sol et le soussol lorsque les écoulements sont à saturation variable. La modélisation est basée sur une version généralisée de la loi de Darcy-Buckingham et de l’équation de Richards multidimensionnelle. Les recherches présentées dans cette thèse concernent différents volets, présentés ci-dessous, dont le fil conducteur est l’analyse de phénomènes d’écoulements en milieux poreux, contenant de fortes hétérogénéités et/ou perturbés par des fluctuations temporelles. Dans le cadre d’un partenariat de recherche entre l’IMFT et l'IRSN sur la problématique du stockage souterrain de déchets radioactifs, nous avons modélisé en 3D la dynamique du front de désaturation d’une couche argileuse autour d’une galerie souterraine ventilée, à l’aide du code volumes finis BIGFLOW 3D. Ce travail a permis de mettre au point une approche d’immersion pour la modélisation des écoulements en milieux composites. Un autre volet de cette thèse concerne une étude analytique et numérique des profils verticaux de succion et de teneur en eau lors d’une infiltration verticale, non-saturante, dans un sol hétérogène finement stratifié. Des solutions analytiques exactes et approchées sont obtenues en régime permanent, à l’aide de transformations de variables, et sont comparées avec des solutions numériques pour différents degrés d’hétérogénéité. De même, les phénomènes de « barrière capillaire » sont étudiés d’une part analytiquement et d’autre part, par expérimentations numériques transitoires d’infiltration sur des systèmes bicouches, en présence d’une nappe plus ou moins profonde. D’autre part, les écoulements nonsaturés sont étudiés, cette fois, en présence d’un forçage transitoire fortement oscillatoire, dans une colonne de sable fin homogène, sous l’effet de fluctuations périodiques du niveau de la surface libre. Une méthode multi-front mise au point lors d’une thèse précédente à l’IMFT est validée numériquement en montrant que cette méthode simule efficacement la dynamique oscillatoire des flux et des profils de pression avec un nombre limité de « fronts ». De plus, nous avons analysé la phénoménologie capillaire/gravitaire des écoulements oscillatoires dans la colonne grâce à un suivi dynamique du plan de flux nul. Enfin, nous présentons, comme extension aux travaux précédents, une étude préliminaire des phénomènes 3D d’infiltration et de redistribution d’eau, notamment lorsque l’hétérogénéité du sol est de type aléatoire. Les premiers essais d’infiltration permettent de tester, en réplique unique, des méthodes de prises de moyennes spatiales des champs de succion et teneur en eau et évaluer la taille minimum du domaine de calcul 3D permettant d’obtenir des profils verticaux moyens représentatifs de l’infiltration dans une réplique unique du sol aléatoire
This doctoral research, defended at the Institut de Mécanique des Fluides de Toulouse, is devoted to modeling water pressure and water content in soils and in subsurface geologic formations, in the case of variably saturated flow. One of the main scientific objectives of this work is to analyze the response of unsaturated flow systems, and particularly capillary and gravitational effects, in the presence of material heterogeneities, discontinuities, and/or space-time fluctuations. Modeling is based on a generalized version of Darcy- Buckingham’s law, and of Richards’ flow equation. Overall, the topics developed in this PhD thesis focus on several related aspects of variably saturated water flow in the subsurface. These aspects all occur at once in most applications (drying/wetting,heterogeneity, temporal forcing), but they are "decoupled" here for convenience. A preliminary research (collaborative project IMFT / IRSN) was developed to study the 2D/3D drying process at the porous wall of a deep cylindrical excavation in response to hydrometeorological signals. This project has motivated the design and testing of a novel approach to include cavities in the modeling domain. A detailed study of steady state infiltration was developed for the case of finely stratified soils, with parameters that vary continuously and cyclically with depth. Exact and approximate analytical solutions are calculated based on variable transformation methods and on perturbation type approximations, and they are tested numerically using a finite volume code (BIGFLOW 3D). The sensitivity of suction fluctuations vs. stratification wavelength is investigated, as well as the effect of the degree of heterogeneity, and of water table depth. Capillary barrier effects are studied for the case of unsaturated infiltration in multilayer soil systems characterized by a discontinuity of soil properties at interfaces. Numerical experiments are developed for transient infiltration towards a water table through a two-layer system, the goal being to analyze possible capillary barrier effects under various scenarios. On the other hand, we study numerically the case of a partially saturated / unsaturated soil column submitted to highly variable oscillatory pressure at the bottom of the column: this leads to vertical flow oscillations in the unsaturated zone above the water table influenced by tides (coastal beach sand). We analyze the dynamics of this oscillatory flow, where capillary and gravitational effects compete; for this purpose we use a novel method that tracks the positions of the zero flux plane in the unsaturated column. Finally, we also present, as an extension to the previous studies, a preliminary investigation of multidimensional infiltration/redistribution phenomena, particularly for the case of fully 3D random-type soil heterogeneity. The first numerical experiments of 3D infiltration are undertaken based on the single realization approach to soil heterogeneity, and assuming a uniform distribution of wetting at soil surface

Stocker du carbone (C) dans le sol représente à la fois des enjeux agronomiques (amélioration des propriétés physico-chimiques du sol) et environnementaux (réduction de la quantité de C02 atmosphérique). A Madagascar, le stockage et la protection de C dans le sol, sous systèmes de culture en semis direct avec couverture végétale (SCY) sont étudiés sur un dispositif de longue durée (11 ans) de la région d'Antsirabe, Hautes Terres malgaches). Quatre systèmes sont étudiés: un système en labour conventionnel (Conventional-tillage) [CT mis, rotation maïs {Zea mays L.)-soja (Glycine max. L.)]. et trois systèmes en SCY sans travail du sol (No-tillage, NT) [NT mis, rotation maïs-soja; NT rn/m-d. rotation maïs-maïs avec une couverture végétale de Desmodium uncinatum ; et NT h/s-k. rotation haricot (Phaseolus vulgaris)-soja avec une couverture de Pennissetum c/andestinum]. Les résidus de récoltes sont exportés de la parcelle pour CT mis alors qu'ils sont restitués pour les NT. Les résultats montrent des teneurs en C significativement plus élevées sous systèmes NT par rapport à CT à 0-5 et 5-10 cm de profondeur. Au-delà de 10 cm de profondeur. les deux modes de gestion ne sont plus différents. Les valeurs de stocks de C à masse équivalente indiquent des stockages annuels élevés pour les systèmes NT étudiés. de 0.7 à 1,0 Mg C.haI.anI pour la couche de 0-20 cm. Ce stockage annuel élevé sous systèmes NT est en grande partie attribué à l'importante quantité de biomasse restituée par ces systèmes. par rapport au système labouré. Les teneurs en C élevées sous NT s'accompagnent d'une teneur élevée en macroagrégats (MA) stables à l'eau, de 16 à 33 % plus élevée par rapport à CT, à 0-5 et 5-10 cm. Ces macroagrégats stables contribuent à plus de 80 % à la différence de teneur en Centre systèmes NT et CT. Ceci pourrait induire une protection physique du C se trouvant à l'intérieur des agrégats contre la minéralisation microbienne. Les fractionnements granulo-densimétriques effectués sur ces échantillons permettent de localiser le C stocké sous NT. A 0-5 cm, 30 à 40 % du C stocké sont associés aux matières organiques particulaires internes aux agrégats stables (MOPi) de taille supérieure à 50 urn et 40 à 60 % aux fractions fines de taille inférieure à 50 urn, A 5-10 cm, 60 à 90 % du C stocké sont localisés dans les fractions fines " 50 urn) du sol. Si la matière organique associée aux fractions fines du sol est protégée physico-chimiquement contre la minéralisation microbienne par son adsorption sur les particules minérales du sol. les MOPi sont protégées physiquement par leur emprisonnement à l'intérieur des agrégats. Par ailleurs, la déprotection de ces MOPi, par broyage grossier " 200 urn) du sol total ou par broyage fin " 50 urn) des macro- et mésoagrégats, ne révèle. qu'une faible quantité de C protégé physiquement (de 26,2 à 54,0 I-lgC.g-1 sol). Ceci pourrait s'expliquer par le fait que ces MOPi, comme le montrent les observations au microscopie électronique à balayage, sont recouvertes et imprégnées de fractions fines qui les protègent physico-chimiquement contre la minéralisation microbienne malgré la destruction des agrégats. En conclusion. pour le sol argileux considéré ici, les systèmes SCY testés s'avèrent très efficaces pour stocker du C dans le sol. Ce C stocké est relativement stabilisé puisqu'H est protégé, au moins physico-chimiquement. contre la minéralisation microbienne.