Brevet FR3039035A1 — publié le 20170127

DESCRIPTION DE L’INVENTION TITRE PROCÉDÉ POUR LA PRÉCONISATION DES DOSES P (dP) D'ENGRAIS PHOSPHORIQUE POUR GRANDUES CULTURES AGRONOMIQUES EN PRÉSENCE DE RÉSIDUS DE CULTURE AU SOL DOMAINE TECHNIQUE DE L’INVENTION Agronomie conseil, et plus particulièrement l'analyse de terres et de végétaux destinés au pilotage des apports aux grandes cultures agronomiques d'engrais N, P et K. Il est aussi question de microbiologie des sols, et plus particulièrement de l'activité diazotrophe des populations azotobactériennes capables de fixation biologique non symbiotique de l'azote diatomique. ÉTAT DE LA TECHNIQUE Préconisation des apports d'engrais P La préconisation des doses d'engrais P et K pour grandes cultures agronomiques est imprécise. Elle s'appuie sur l'élaboration de courbes de réponses mitscherlichienne des rendements agronomiques de ces culture aux teneurs du sol en P dit échangeable ou encore (phyto)disponible renforcées de doses d'engrais P calibrées. Ces relations mitscherlichiennes (MITRDT) sont de la forme (eg. Ussiri et al. 1998, Morel et al. 1992) ; RDT = A ■ [1 - exp(-c ■ (b ■ P + dP))] (1 ) lorsque A, b et c sont des coefficients empiriques, RDT le rendement agronomique de la culture (t/ha, qx/ha ou encore g/m2), P la teneur en P "échangeable" du sol et dP la dose d'engrais P apportée, P et dP ayant comme de raison les même unités, généralement en mg-P/kg-sol ou kg-P / ha. Par transformation, il est ainsi possible de rétro calculer dP nécessaire à un rendement équivalent à 100% du potentiel de rendement, à savoir (cf. Ussiri et al. 1998) ; dP = (1 / c) ■ ln[c · A / p ■ (1 + R)] - b ■ P (2) lorsque p est le ratio entre les prix des engrais et de la récolte, et R le taux marginal de rendement économique (a lias taux d'intérêt ; marginal rate of return). Il est ainsi, en principe, possible de préconiser des dP marginalement de plus en plus faibles selon l'augmentation régulière des teneurs P des sols. Cette démarche à été appliquée en France par Morel et al. 1992. Cette approche MITRDT est néanmoins problématique. Premièrement, la réponse des RDT à P et dP n'est pas toujours très marquée (Figure 1 ). En effet, le P apporté au sols est rapidement fixé physico - chimiquement sous forme de composés calco - magnêsiques, ou encore alumino ferriques si les sols sont particulièrement acides. Il en résulte que les engrais P, à l'exception notable des engrais P localisés près des semences et/ou racines (engrais "starter"), ne servent pas à alimenter directement les racines, mais à enrichir le sol ; au delà d'une certaine richesse en P, il sera ainsi recommandé de faire l'impasse sur de tels apports d’engrais P (Comifer 2007). L'établissement d'un indice de rendement doit donc impérativement comprendre la compilation de résultats de nombreux sites et années différentes ; étant donné les différents potentiels de rendement en oeuvre, les rendements observés doivent être normalisés, généralement par rapport à un rendement maximal de 100% (eg. Morel et al. 1992). Cette normalisation sert donc à masquer, en sorte, la très grande variabilité interannuelle et inter - sites des réponses des rendements à P et dP. La fertilisation P, voire K, spécifique à une culture agronomique est fondée sur un certain CME-PK (coefficient multiplicatif d'exportation - phosphore/potassium ; Comifer 2007) appliqué au rendement agronomique moyen de la parcelle et ses teneurs en P et K pemettant de calculer une dose hectare de P et/ou K. Ces recommandations sont déclinées selon l 'exigence de la culture en P et K. Or, cette notion d'exigence est redondante par rapport à celle de rendement agronomiques, et donc d'exportation de P et K. De plus, la valeur du CME-PK dépend d'une relation RDT : dPKMg empirique, dite ici mitscherlichienne , en principe spécifique à une culture, mais aussi à un type de sol, à un régime hydroclimatique, etc. Or, ces relations n'existent pas pour l'ensemble de ces situations agronomiques particulières. Leur établissement à partir de données pluriannuelles et multisites étant très onéreux (Colomb ét al. 2007). Il est donc difficile d'établir le lien entre la disponibilité élémentaire du P sans avoir à élaborer une certaine relation mitscherlichienne entre RDT et la phytodisponibilité du P (MITRDT). Ces relations RDT : phystodisponibilité du P fonction des teneurs en apports de P firent l'objet de maintes études de calibrations (Morel et la. 1992), sans pour autant permettre à ce jour la création de grilles de fertilisation véritablement universelles, même au sein d'une même SAU nationale. Localement sur un site ou un sol en particulier il est très difficile de résoudre simultanément ce type d'équations mitscherlichiennes non linaires à deux variables, teneurs et apports en P. En effet, l'effet agronomique de tels apports de P sont a priori peu ou pas détectables ; il est surtout question d'enrichissement graduel et à moyen et long terme du sol dans l'espoir que par après de tels sols enrichis nécessitent moins, ou pas, d'engrais P. Il ne s'agit pas ici d'un problème expressément liée à l'agriculture de précision cherchant à adapter superficiellement par modulation en temps réel des dose P, K et/ou Mg selon les variations immédiates du terrain, mais plutôt du recalibrage des grilles de fertilisation P à ce jour trop dépendantes de fonctions MITRDT peu ou pas renseignées. Diazotrophie non symbiotique, P du sol et PUE Bien que la diazotrophie soit très énergivore, et donc grande utilisatrice de P via la phosphorylation de l'ADT en ATP, une telle fixation non symbiotique de l'azote en présence de résidus de culture au sol (RCS) n'a jamais, contrairement au potentiel RDT, servie à calibrer la fertilisation P. Augusto et al. 2013 rapportent pourtant une simili courbe MITAZB - i.e. elle ne résulte pas d'apports d'engrais P pré - calibrés, sans pourtant pouvoir donc traduire ces réponses mitscherlichiennes en préconisations P ; de plus il s'agit de diazotrophie symbiotique hors du domaine d'application de la présente invention. Autre constat. L'efficience agronomique du P ( alias ; phosphorus use efficiency, ou PUE) apporté sous forme d'engrais aux culture n'est apparemment pas améliorée en présence de RCS, du moins la base de connaissance publiée n'en fait pas mention. Malgré une contribution aux réserves de P du sol (Damon et al. 2014, Thibaud et al. 1988), ces RCS auraient surtout, semble-til, un effet dépressif sur PUE, vraisemblablement par entraînement suite à l'effet dépressif de ces RCS sur NUE (nitrogen use efficiency) en raison de l'immobilisation de l'azote minéral. Par contre, la contribution de ce P à l'augmentation de TARA (activité de réduction de l’acétylène) non symbiotique en proximité des RCS n'est pas documentée. Bien qu'étonnante, cette omission peut néanmoins être attribuée au fait que TARA in situ non symbiotique, y compris en proximité des RCS, fut jusqu'à ce jour perçue comme peu ou pas importante en termes agronomiques, contribuant aux sols arables tempérés que quelques kg N par hectare et par an. Depuis remise en question par quelques chercheurs (Gupta et al. 2006 ; Halsall et Gibson 1983), cette affirmation à néanmoins marqué les esprits ; TARA in situ non symbiotique est trop faible pour être intégrée aux méthodes de préconisations des doses d'engrais N, P et K pour grandes cultures. Enfin, notons que le PUE d'engrais P, hors engrais starter, est assez faible de l'ordre de 10 à 20% du fait d'un haut taux de fixation physico - chimique des ions phosphatés ainsi apportés directement au sol (eg. FAO 2008, Sattari et al. 2012). Toute méthode pour la préconisation de P sera donc la bienvenue. Cela est d'autant plus vrai que la mobilisation du P pour la culture et/ou la richesse en P des grains n'est pas en soi un objectif agronomique. Au contraire, dans un contexte d'épuisement des gisements de P - ressource non renouvelable, la réduction des taux de mobilisation du P par les graines (Rose et al. 2013) est proné. Le bouclage du cycle du P de manière à retourner au sol le P peu ou pas assimilables (acide phytique) pourrait globalement réduire de beaucoup les besoins en P des sols arables. Pour preuve, Sattari et al. 2012 (fig. 5 / p 6352) démontrent malgré une réduction global, notamment en Europe de l'ouest (op.cit. fig. 2 3 p 6350), les niveaux de mobilisation du P par les culture est stable, voire en légère augmentation ; selon les auteurs, cette tendance devrait ce maintenir jusqu'en voire au delà de 2050. Minorer les doses d'engrais P de manière à assurer ponctuellement, et non seulement globalement comme le rapportent Sattari et al. 2012, et donc d'actualité agronomique. Sigles et définitions AZB : azotobactérisation des résidus de culture cellulosique au sol (RCS) - eg. Claude et Fillon 2004. pAZB : potentiel l’activité diazotrophe en surface du sol arable, y compris celle attribuable à la susdite azotobactérisation des résidus de culture. L'pAZB est ici calculée comme la différence en l'pAZBn et ΓρΑΖΒρ, i.e. les rapports des ratios C sur N et C sur P (rCN, rCP) résidusphériques calculés ou mesurés et des ratios C sur N et C sur P critiques (rCNc, rCPc) en terme de libération de métabolites C, N et P assimilables par les microorganismes du sol. L'immobilisation de l'azote du fait d'un rCN supérieur à rCNc étant favorable à la diazotrophie, pAZBn est désigné comme positif ; in contrario, l'immobilisation du phosphore du fait d'un rCP supérieur au rCPc, pAZBp est désigné comme négatif, d'où la susdite différence pAZBn - pAZBp. MITRDT : relation de type mitscherlichienne (MIT), i.e. non linaire et multivariable, comportant comme variable dépendante le rendement (RDT) de la culture agronomique et comme variables indépendantes la teneur en P, ou tout autre élément fertilisant en principe, du sol, et divers niveaux d'apports d'engrais P au sol (dP). La relation atteint un plateau à titre d'optimum (cf. éq. 1) ; il est ainsi possible de rétro calculé une certaine dP selon la teneur en P du sol (cf. éq. 2). MITAZB : relation de type mitscherlichienne (MIT), i.e. non linaire et multivariable, comportant comme variable dépendante le niveau d'activité, potentiel ou réel, de la flore azotobactérienne (AZB) du sol, et en particulier en proximité des RCS, et comme variables indépendantes la quantité de P, ou tout autre élément fertilisant en principe, retourné au sol via les RCS et divers niveaux d'engrais P (dP) apportées directement aux RCS. La relation atteint elle aussi un plateau à titre d'pAZB optimum (cf. éq. 3) ; il est ainsi possible de rétro - calculé une certaine dP selon la susdite quantité de P retournée au sol via les RCS de la culture précédente (cf. éq. 4). RCS : résidus de culture cellulosiques, avantageusement pailleux, au sol, y compris de maïs - grain, bien que les RCS de cultures de tournesol et de colza sont s'avèrent aussi de bon substrats pour l'AZB. Il y a ici trois usages possible de ce sigle ; qRCS (kg de RCS par hectare), c CS (% de couverture du sol par des RCS), et iRCS (indice de dispersion, répartition ou amoncellement des RCS) iPAZ : méthode pour la préconisation de doses d'engrais phosphoriques (dP) selon la présente invention caractérisée en ce que (i) elle s'appuie sur MITAZB plutôt que MITRDT, et (ii) les dP sont appréciablement plus faibles que celles autrement préconisées à titre de fumure de fond PK. dose X (dX ; kg-N, ou "unités d'azote" / ha) : soit = [coefficient d'utilisation de ΓΝ spécifique à une variété de la culture x l'objectif de rendement], mais aussi très complexe comprenant une foule d 'écritures opérationnelles et un certain bilan prévisionnel d'azote (Comifer 2011), la plupart de ces écritures devant être grossièrement estimées à l'aide d'abaques plus ou moins bien renseignés. dose P (dP) : quantité d'engrais P, en équivalents P, apporté au sol, ou plus ingénieusement ici directement aux RCS avant leur enfouissement. Dans le cadre de la présente invention, il s'agit surtout de fumure de fonds, y compris si elles sont appliquée directement sur les RCS, et généralement pas d'engrais dits "starter" normalement appliqués en bandes et/ou près des semences. Ceo : quantité de substrat carboné, ici selon Poil et al. 2006 en kg/m3, apporté initialement au sol via le retour des RCS - pour fin de modélisation, j'ai supposé qu'il s'agit à court terme surtout de l'équivalant physicochimique (cf. coefficient de diffusion) du cellobiose. Cno : quantité de substrat azoté, ici selon Poil et al. 2006 en kg/m3, apporté initialement au sol via le retour des RCS - pour fin de modélisation, j'ai supposé qu'il s'agit à court terme surtout de l'équivalant physicochimique (cf. coefficient de diffusion) de l'ammonium. Cpo : quantité de substrat phosphorique, ici selon Poil et al. 2006 en kg/m3, apporté initialement au sol via le retour des RCS - pour fin de modélisation, j'ai supposé qu'il s'agit à court terme surtout de l'équivalant physicochimique (cf. coefficient de diffusion) des phosphates. Cpoi : quantité de P retournée au sol via la les RCS ; a priori, Cpoi = Cpo - la distinction est faite ici pour fin de modélisation mathématique via R dPo : fraction de Cpoi définissant un apport en équivalents de P, dCpoi, comme suite dCpoi = dPo x Cpoi dCpoi : quantité de P à apporter aux RCS, avantageusement sous formes solubles en solutions ou en granulés (voir texte et Tableau 7), de manière à optimiser pAZB au sens de la relation MITAZB (éq. 3) - dCpoi peut maintenant être traduite en dP selon l'équation (4) dP = -b · Cpoi DIVULGATION DE L’INVENTION Problème technique Le problème technique provient d'un calcul imprécis des apports d'engrais P préconisés sur la base de relations mitscherlichiennes de type MITRDT. Ces relations MITRDT supposent d'éventuelles augmentations de RDT suite à l'enrichissement en P du sol via l'apport d'engrais P soient appréciables et documentées pour un même site. Or, sur un même site, il y a trop peu d'observatons pour tracer intelligemment une courbe de type MITRDT. Il faut donc compiler des observations provenant de nombreux sites / années aux potentiels de rendement très différents, quitte à les normaliser, par exemple sur une base 100. C'est cette normalisation qui cause ici problème puisque les rendements escomptés sont eux absolus (eg. quintaux / hectare), et non relatifs, tout comme doivent l'être les préconisations d'engrais P. De plus, l'homme du métier parle de dP "renforcement", et de dP "enrichissement", l'impasse de l'une ou l'autre étant aussi possible si le sol est jugé suffisamment riche en P (Comifer 2007). Or, le passage du renforcement à l'impasse en P est dans une large mesure imprécis. J'attribue ce flottement à une application trop stricte du paradigme mitscherlichien, initialement prévu pour la fertilisation azotée et non phosphatée, à la préconisation de doses d'engrais P. La base de connaissance décrit donc des relations mitscherlichiennes par rapport aux rendements agronomiques (MITRDT). Or, et bien que l'activité azotobactérienne du sol (AZB), ou plus exactement le niveau d'ARA qu'elle suppose, puisse nécessairement bénéficier d'une plus grande disponibilité du P, des relations mitscherlichienne par rapport à AZB (MITAZB) de la même forme que MITRDT ne font pas partie de l'état de la technique, et ne peuvent donc pas se substituer aux usuelles relations de type MITRDT pourtant bien imparfaites. Dommage. In fine , toute cette problématique contribue à des niveau de PUE trop faibles pour être durables. Solution technique L'invention (iPAZ) affirme que ce sont les teneurs en P des RCS qui surtout alimentent en P, non les culture mais plutôt les flores azotobactériennes (AZB) résidusphériques. Il s'avère que TARA des ces AZB est toute aussi, voire plus, révélatrice de l'état phosphorique du sol que le potentiel de rendement des culture. La calibration des dP selon le potentiel d'activité des AZB résidusphériques (pAZB) s'avérera tout aussi utile en agronomie que la construction fastidieuse et tout relative des habituelles courbes MITRDT. En effet, le niveaux d'activité de la diazatrophie non symbiotique du sol en proximité des RCS enfouis et en cours de décomposition (résidusphère) est limité par la quantité de P disponible dans le sol, beaucoup plus que le sont les rendements agronomiques de grandes cultures. Les courbes de réponses MITRDT peuvent donc être remplacées par des courbes de réponses MITAZB décrivant non les augmentations des rendements agronomiques en réponse au renforcement de P par dP, mais plutôt les augmentations du pAZB en réponse au même renforcement de P par dP. Cette solution technique comprendra l'application, le cas échéant, de dP directement aux RCS et non plus au sol, réduisant ainsi notablement les risques de fixations physico - chimiques, calco-magnésiques et/ou aluminoferriques des engrais P apportés au sol. Comme nous le verrons, ces apports résidusphériques de P seront aussi appréciablement moins importants que les habituels apports de P sous forme de fumure de fonds destinés avant tout à enrichir en P les sols dans l'espoir de justifier ainsi une éventuelle impasse sur rapports d'ultimes dP. Il s'agit donc d'un procédé pour la calibration des apports aux cultures agronomiques d'engrais phosphoriques comprenant une courbe de réponse mitscherlichienne décrivant l'activité diazotrophe de la flore azotobactérienne dans un sol arable en proximité de résidus de la culture précédente enfouis. Cette courbe de réponse est dite ici de type MITAZB. Cette calibration ce fait en fonction, (i) des teneurs en C, P, voire N, desdit résidus de culture cellulosiques, avantageusement pailleux, au sol (RCS) avant leur enfouissement en pré - semis d'une culture agronomique suivante et, (ii) d'un équivalent de phosphore exogène, dit ici dose P ou dP, proportionnel aux susdites teneurs en C et P des résidus de culture au sol. Ce procédé est aussi caractérisé en ce qu'il nécessite ; > la détermination de la quantité desdits résidus de culture au sol (qRCS), soit indirectement via une relation empirique fonction du rendement de la culture précédente (RDT), soit plus directement au sol via la détermination du degré de couverture du sol par ces résidus (cRCS) ; > la détermination des teneurs en N (azote) et en P (phosphore) desdits RCS de manière à obtenir Cpoi, i.e. la quantité de P initialement retournée au sol par ces RCS > le calcul de dP en fonction de Cpoi capable d'assurer que le niveau d'activité diazoptrophe en proximité desdits résidus de culture une fois enfouis (pAZB) maximale. La détermination de la teneur en P peut ce faire lorsque la culture précédente est encore sur pieds de manière à pouvoir calculer dP avant la récolte et le retour des RCS. Cette détermination anticipe de dP est comme de raison logistiquement avantageuse ; elle est aussi fondée du fait qu'il exist une certaine corrélation entre les teneurs en P des grains et paille au stade BBCH 65 et celles à BBCH 89-91. Plus simplement, et sans avons à dépendre de cette corrélation, il et aussi possible d'échantillonner les RCS toujours sur pieds une fois la matruité physiologique, soit aux alentours de BBCH 75-77 ; le cas échéant, l'utilisateur dispose d'environs trois semaines, dans la cas d'une culture d'hiver telle que blé tendre en France, avant la récolte, voire de 4 à 5 semaines avant l'enfouissement des RCS, pour obtenir une valeurs %P rC s et de Cpoi. Le niveau d'activité diazoptrophe en proximité desdits résidus de culture une fois enfouis (pAZB) est défini comme l'équilibre entre l'immobilisation nette de l'azote du sol, phénomène favorable à ladite activité diazotrophe, et l'immobilisation nette du phosphore du sol, phénomène défavorable à ladite activité diazotrophe. Cet équilibre entre l'immobilisation nette de l'azote du sol, favorable à la activité diazotrophe, et l'immobilisation nette du phosphore du sol, défavorable à la activité diazotrophe est défini formellement comme le rapport entre pAZBn et pAZBp (pAZBn / pAZBp), et plus avantageusement comme la différence entre pAZBn et pAZBp (pAZB - pAZB). pAZB -n et -p sont eux, respectivement définis comme les ratios rCN/rCNc et rCP/rCPc, i.e. les rapports entre ratios C/N (ou C/P) de la résidusphère à l'équilibre et C/N (ou C/P) critique en dessus du quel l'azote (ou le phosphore) est libéré sous forme minéral (minéralisation). Le calcul de dP en fonction de Cpoi et d'un coefficient b issue d'une relation MITAZB de type pAZB = A ■ [1 - exp(-c ■ (b ■ Cpoi + dCpoi))] lorsque dCpoi et un équivalent du P apporté aux RCS et/ou au sol affectant pAZB, et dont la valeur est comprise entre 0,10 et 0,50, plus particulièrement entre 0,15 et 0,30, et avantageusement ici d'environs 0,20, peut être effectué comme suite ; dP = -b x Cpoi dP pouvant ainsi être exprimé en kg-P par m3 de sol, mg-P / kg-sol, ou encore en kg-P par hectare une fois les conversions unitaires nécessaire effectuées (eg. Figure 5). La détermination de la quantité des résidus de culture cellulosiques au sol (qRCS ; g/m 2 ) est elle déterminée indirectement à partir du rendement de la culture précédente (RDT ; g/m 2 ) selon une des relations empiriques suivantes ; qRCSsjH = [0,911 x RDT + 164] qRCSmais = [1 ,51 x RDT + 3,7 - RDT] qRCS s0 rgho = [2,56 x RDT + 0.4 - RDT] qRCScéréaies = [2, 13 x RDT + 1 , 1 - RDT] , Il est aussi possible de déterminer qRCS plus directement via la détermination préalable du degré de couverture du sol (cRCS ; %), la quantité de RCS (qRCS ; t/ha) pouvant être déterminé par l'une des relation empiriques suivantes ; qRCS BT H = 0,0633 x exp(0,0493 x cRCS) | R2 = 0,9275 qRCS M aïs = 0,0519 x exp (0,0554 x cRCS) | R2 = 0,8647 qRCSroumesoi = 0,0015 x cRCS 2 - 0,003 x cRCS + 0,378 | R2 = 0,9928 Les cultures agronomiques concernées sont des grandes cultures agronomiques d'hivers choisies parmi un groupe de cultures comprenant les céréales, y compris le blé, l'avoine, le seigle, le triticale, l'orge, le colza, le lin, ou encore toute autre grandes culture agronomique pouvant être semis en présence de résidus de culture cellulosiques, avantageusement pailleux au sol. dP est sous forme soluble est appliquée par pulvérisation liquide directement aux RCS avant leur enfouissement, avantageusement en pré - semis d'une grande culture agronomique. dP peut aussi être sous forme pulvérulente applicable par pulvérisation solide directement aux résidus de culture cellulosiques au sol avant leur enfouissement, avantageusement en pré - semis d'une grande culture agronomique. Dans le cadre de cette méthode de fertilisation phosphorique, il est aussi possible que les résidus de culture cellulosique au sol soient bactérisés à l'aide d'une préparation bactérienne réputée capable d'augmenter l'activité diazotrophe, voire l'activité de réduction de l'acétylène (ARA), du sol en proximité desdits résidus de culture un fois enfouis, par exemple selon Claude et Fillon 2004. Mieux, la préparation bactérienne contient des biomasses bactériennes capables de minéraliser partiellement en P inorganique le P organique, notamment sous forme de sels d'acide phytique ; plusieurs biomasses bactériennes de ce type existantes en ce sens sont rapportées dans la base de connaissance, à savoir par Idriss et al. 2002. Enfin, la préparation bactérienne peut être combinée avec dP au sein d'une même bouillie pulvérisable directement sur lesdits résidus de culture au sol. L'osmotolérance relative de ces biomasses, très recherchée en de telles circonstances, peut elle être assurée via la technologie décrite dans la demande de brevet EP2647695 (A1), à avoir un mode de préparation de bouillies inoculantes pulvérisables comprenant des populations bactériennes moyennement halophiles Avantages apportés et activité inventive Il n'est plus nécessaire d'obtenir une relation de type MITRDT reliant les apports d'engrais et teneur en P des sols aux rendements agronomiques. En effet, ce type de relation est trop souvent peu marquée et nécessairement relative du fait d'avoir à compiler des données de divers sites x années. Au contraire, l'activité diazotrophe de la flore azotobactérienne du sol est elle limitée par la disponibilité du P ; cette activité est donc un meilleur indicateur des besoins en P du sol, notamment s'il y azotobactérisation des RCS. Le potentiel agronomique de telles biofertilisations azotobactériennes est donc augmenté, et l'impacte agronomique de tels apports de P directement aux RCS sera plus apréciable que celle d'apports de P au sol sous forme de fumure de fond. De plus, ces doses P (dP) pulvérisables directement sur les RCS, seront appréciablement plus petites que les dP normalement préconisées sous forme de fumures de fonds PK. Cette réduction des dP, hors engrais "starter", est vraisemblablement le fait d'une meilleure utilisation du P par les flores AZB in situ en proximité immédiate de RCS. BRÈVE DESCRIPTION DES DESSINS ET FIGURES Figure 1 : Illustration d'un aspect du problème technique. Les corrélations, y compris mitscherlichiennes, entre le rendement d'une culture agronomique (RDT) et la disponibilité du P dans le sol, ici selon Dyer ou Olsen, sont a priori vagues. {Haut) 107 points variables RDT (blé tendre d'hiver, ou BTH) : P 2 05 Dye r· {Bas) idem, mais par rapport à P 2 0 5 oisen. Ces deux jeux de données RDT : P proviennent d'une compilation des rendements moyens pour l'ensemble des départements et régions françaises (données Agrestes, Ministère de l'Agriculture, France) adjoints aux teneurs P 0 isen et P Dye r correspondants. Ces valeurs proviennent de la Base de Données d'Analyses de Terre (BDAT : Gissol - Groupement d'intérêt scientifique pour les sols). Figure 2 : Coefficients de régression normalisés de pAZB (pAZBn - pAZBp) sur quantité et degré de couverture des RCS (kgRCS et cRCS, respectivement), l'indice iRCS, le ratio C/N et la teneur en N des résidus de culture (rCNo et %Nrcs, respectivement), ainsi que le reliquats d'azote minéral (Nm) à l'automne en ppm (mg/kg) d'Nm et, surtout, l'ampleur de l'apport de P au RCS en proportion de la quantité de P retournée au sol par ces RCS (dPo). Cette dernière proportion sera par la suite être reconvertie en dose d'engrais P (dP ; équation 4) soluble pulvérisable directement sur les RCS. Notons que ces coefficients normalisés comparent le poids relatif des variables ; pus la valeur d'un coefficient est élevée, plus le poids de la variable correspondante est important. Figure 3 : Module de rétroaction négative de la diazotrophie non symbiotique de la flore azotobactérienne du sol en présence de résidus de culture cellulosiques. Deux type de rétroactions sont en jeux ; la première tenant compte d'une accumulation graduelle de l'azote minéral au delà de 14 mg-N / L de solution du sol (1mM-N) capable de réduire appréciablement le fonctionnement de Nif (en abîme basse), et la deuxième fonction en puissance négative du ration NP des résidus de culture au moment de leur enfouissement (rNPo ; en abîme haute). Cette double rétroaction permet de maintenir la contribution de l'activité azotobactérienne (Ndiazo), ici sans l'action directe des éventuels inocula azotobactériens, à un régime a priori modeste en accord avec l'état de l'art. Figure 4 : Efficacité relative (er ; alias "potentiel") de la diazotrophie azotobactérienne (AZB) du sol selon l'augmentation de l'apport de P proportionnel à Cpoi (kg-P / m3 de sol), i.e. l'apport de P au sol provenant des résidus de culture cellulosiques provenant de la culture précédente. Cette apport supplémentaire de P, dCpoi, proportionnel à Cpoi, permet d'atteindre progressivement des pAZB égale à 0,00, sachant qu'pAZB = pAZBn - pAZBp, le pAZB négatif étant ici progressivement réduit en raison de la diminution du ratio rCP / rCPc dictant l'immobilisation du P. En effet, ici dCpoi réduit rCP sans affecter rCPc qui est lui strictement fonction de la qualité des résidus de culture avant l'application de l'équivalent dPo. Le texte de la description explicite plus en détails cette action de dCpoi sur pAZB. Figure 5 : Courbe de calibration des équivalents P, dCpoi (kg-P / ha), à apporter à Cpoi au moment de l'enfouissement des résidus de culture cellulosiques au sol de manière à ce que pAZB atteigne approximativement 0,00, soit ici le plein potentiel de l'activité azotobactérienne non symbiotique du sol du fait d'un certain équilibre entre l'immobilisation de l’azote et la non immobilisation du P ; cette notion est explicitée dans le texte. Notons que dCpoi représente ici constamment 20% de Cpoi, ce qui est facilement réalisable en agronomique, Cpoi ne dépassant pas 7 - 10 kg-P / ha / année. dCpoi ayant permit l'esetimation du coefficient b (équation 3), ces apports de P pourront être exprimée en dP (mg-P / kg-sol, ou kg-P / ha). MODE DE RÉALISATION PRÉFÉRÉ DE L’INVENTION Le potentiel e l'activité azotobactérienne non symbiotique (pAZB) en sols arables est limitée par la disponibilité du P (Reed et al. 2007, 2011 et 2013, Binkley et al. 2003, Augusto et al. 2003, Vitousek et al. 1991 , Eiselel et al. 1989, Howarth et Marino 1988). Le P apporté au niveau de la résidusphère représentera, ici à titre d'exemple non limitatif, de 10 à 25% du P normalement apporté au sol via les RCS. Ce P peut être maintenant apporté par pulvérisation liquide des RCS avant leur enfouissement ; notons que cette pulvérisation liquide d'une solution de P fut habituellement destinée au feuillage. J'ai prévue une gamme de scénarii comportant de 4 à 9 t de RCS, des rCNo de 33 à 125, ainsi que des reliquats d'azote minéral automnaux (Nm) de 10 à 25 mg-N/kg-sol. Ces divers scénarii au moment de l'enfouissement des résidus de culture, avantageusement à l'automne en pré - semis d'une culture hiver telle que le blé, sont décrit au Tableau 1. Ceux-ci, afin de rendre plus défendable la faisabilité de réalisation de l'invention, comportent différents taux de couverture au sol (cRCS ; %) de répartition plus ou moins aléatoire selon l'indice iRCS = [dcm · iEE-4], Notons aussi à ce stade qu'il n'est pas impératif de préciser l'espèce des RCS puisque les scénarii au Tableau 1 englobent a priori la plupart des cas de figures comprenant des RCS plus ou moins abondants, plus ou moins riches en N, et/ou plus ou moins amoncelés à la surface du sol. Il sera néanmoins possible à terme d'envisager des courbes de calibration de dP selon Cpoi (Figure 5) particulière aux différentes espèces produisant les RCS. Description de la modélisation iPAZ-R Entrée des data pour les différents scénarii (Tableau 1) Les statistiques descriptives pour les différents scénarii, 1 500 en tout, sont au Tableau 1. Ces données sont placées dans un fichier - data iPAZ.txt. En principe, tout les cas de figures concernant le retours de RCS par des grandes cultures agronomiques telles que le blé, le colza et le maïs sont ainsi couverts, bien que certains, c'est vrai, sont peut probables, eg. 9 t/ha de résidus à des rCN très faible ; il est ici surtout question de tester la robustesse de telles modélisations. Tableau 1 : Données initiales pour le modèle iPAZ-R ( alias input data) - data iPAZ.xls Echantillon Unités Minimum Maximum Médiane Moyenne écartype (n-1) Moynne géom. RCS kg/ha 3947 cRCS % 8,00 dmm mm 13,87 iEE-4 pixels 1 ,41 iRCS sans unités 26,00 rCNo sans unités 33,00 RDT g/m 2 42,39 Nrcs % 0,336 ppmNm mg/kg-sol 10,72 dPo kg/ha 0,010 9000 6421 6360 60.00 35,00 32,80 31,88 21,79 21,26 4,38 2,50 2,50 140.00 46,00 55,40 125.00 75,00 76,60 306,81 171,83 168,66 1,273 0,560 0,686 25,44 19,71 19,24 0,250 0,125 0,127 1791 18,51 5,64 0,88 32,61 33,18 93,70 0,340 4,16 0,083 6094 26,50 20,52 2,35 48,26 68,85 135,63 0,610 18,75 0,085 Concrètement, dix (10) niveaux de RCS allant de 4000 à 9000 kg/ha sont croisés avec six (6) différentes doses de P (dP) allant de l'équivalent de 1 à 25% la quantité de P retournée au sol par les RCS, et cinq (5) niveaux de richesse relative en N exprimées sous forme de ratios C/N (rNCo). Il y a donc a priori 10 x 6 x 5 = 300 observations simulées. Pour démontrer que l'indice iPAZ ne dépend pas des teneurs en N reliquats automnaux, j'ai aussi inclus cinq (5) valeurs aléatoires de tels reliquats (en mg-Nm / kg-sol) variant de 10 à 25. Or, comme nous le voyons à la Figure 2, une régression linéaire de l'efficacité relative de la flore azotobactérienne du sol (pAZB) sur ces variables démontrer qu'iPAZ n'est pas fonction, même partiellement, de tels reliquats. Il n'y a donc pas redondance avec d'autre type d'indices intégrant ces reliquats d'Nm. Les résultats de ces modélisations iPAZ-R (Tableaus 2-a, 2-b et 2-c) sont exprimés en terme de concentrations (kg élémentaire / m 3 de sol ; Poil et al. 2006 ; infra Tableau 4) de carbone, d'azote et de phosphore au moment de l'enfouissement soit, Ceo, Cno et Cpo, respectivement ; les ratios C/N, C/P et N/P à ce moment, soit respectivement rCNo, rCPo et rNPo, peuvent ainsi être calculés et comparés aux ratios critiques correspondant, rCNc, rCPc et rNPc, afin de calculer Ι'ρΑΖΒ ; pAZBn Les résultats de ces modélisation iPAZ-R comprennent aussi les concentrations (kg élémentaire / m 3 de sol ; Poil et al. 2006) de carbone, d'azote et de phosphore 90 jours post enfouissement (jpe) soit Cn, Cn et Cp, respectivement. Le modèle iPAZ-R tente de décrire l'équilibre entre r(im)mobilisation de l'azote, favorable, et du phosphore, défavorable, à l'activité azotobactérienne non symbiotique. Formellement, rCNc et rCPc sont, respectivement, les ratios C/N et C/P critiques pour la biomasse résidusphérique au delà desquels la minéralisation du N et P, respectivement, ne pourra libérer de N et/ou de P minéral ; ces valeurs critiques sont calculés à partir des équations empiriques proposées par Manzoni et al. 2008, à savoir ; rCNc = 2,3475 x cCNo 0,7452 rCPc = 2.3223 x rCPo 0,8430 rNPc Établissement des coefficients de diffusion à travers le sol pour C, N et P (Tableau 3) Les concentrations effectives de C, N et P en proximité des résidus de culture enfouis en cours de décomposition (résidusphère) dépendent des coefficients effectifs diffusion de ces divers substrats. Ces substrats vont migrer à des distances différentes, affectant ainsi leurs concentrations disponibles aux azotobacters capables d'AZB non symbiotique. Ce coefficient De (ou D * a selon Habib et Guennelon 1983) est lui calculé comme suite ; De = Do / (1 + ((p d A bdc) / Θ)) sachant que Do est le coefficient moléculaire du substrat dans l'eau, et b est un coefficient linéaire décrivant le taux d'adsorption du substrat in situ ; p d et Θ sont les densité apparentes du sol et la teneur volumétrique en eau (v/v) du sol. À ce stade p d et Θ sont fixés à 1 ,25 (g/cm 3 ) et 0,25 (v/v), respectivement. Plus précisément, les métabolites - substrats C, N et P sont ici, respectivement, la cellobiose, l'ammonium et les phosphates, soit les métabolites communément issus de la dégradation microbiologique in situ des résidus de culture cellulosiques, voire avantageusement ici pailleux. Les valeurs Do et De calculées selon Habib et Guennelon 1983 sont elles rapportées au Tableau 3. Tableau 2-a : Modèle iPAZ-R - Code R pour la modélisation du potentiel AZB (pAZB) ; Intégration des scénarii (Tableau 1) et paramétrage des équations de diffusion et de consommation (dégradation) des substrats C. N et P selon leurs distance de migration à l'équilibre. rm(list=ls()) data Doc # Habib et Guennelon 1983 éq. 12 (v/v) # Habib et Guennelon 1983 éq. 12 (g/cm3) # m2/jour (Wang et al. 2011 , tableau 1) # Neff et Asner 2001 , tableau 2, annexe # m2/jour (Ccokeborne et al. 1988, tableau 5) # NGWCLC 2003, tableau 3.2 # m2/jour (Habib et Guennelon 1983, page 118) # Habib et Guennelon 1983 éq. 7, figure 3 # Poil et al. 2006 (tableau 2, figure 2) : kD # Poil et al. 2006 (tableau 2, figure 2) : kD # Poil et al. 2006 (tableau 2, figure 2) : kD # Poil et al. 2006 (tableau 2, figure 2) : vD # Poil et al. 2006 (tableau 2, figure 2) : vD # Poil et al. 2006 (tableau 2, figure 2) : vD # kg de RCS # % de couverture du sol par les RCS # distance entre les foyers d'activités (cm) # i.e. mi-distance entre deux foyers d'activité # "intensité des foyers d'activité selon kde2d" - EEfor (n in 1 :length(data[,1])) { RCS 04 iRCS RDT Nrcs A 1.16 rNPo A 0.7452 rCPc A 0.8430 rNPc # Manzoni et al. 2010 (figure 3) # Manzoni et al. 2010 (figure 7c) # Manzoni et al. 2010 (figure 7d) Cco Cpoi # kg-C / m3 de sol sur 10 cm de profondeur # kg-N / m3 de sol sur 10 cm de profondeur # kg-P / m3 de sol sur 10 cm de profondeur # ... besoins en P pour atteindre pAZB ~ 0,00 Tableau 3 : Paramétrage du modèle iPAZ-R (cf. Tableaux 2-a, b et c) C N P unités / référence / cf. Métabolite Do De Bd De/Do cellobiose 1,5 0,207 2 0,138 ammonium 1,5 0,114 5 0,076 phosphates 1,5 0,039 10 0,026 m2 / jours x 10e6 m2 / jours x 10e6 Habib et Guennelon 1983, eq. 7 na coeff. ABC 3,33 (A) coeff. abc 3,65 (a) 5.00(B) 10.0(C) 3,65 (b) 3,65 (c) Au sens de Poil et al. 2006 Au sens de Poil et al. 2006 Références Neff et al. 2001 NGWCLC 2003 Habib et Guennelon 1983 Les coefficients de diffusion moléculaires (Do) sont essentiellement les mêmes pour les différents métabolites - substrats, tandis que les coefficients de diffusion effective (De) sont eux assez différents du fait des valeurs bd variables. Les valeurs de bd applicables à l'équation 7 de Habib et Guennelon 1983 ont été estimées, respectivement, pour les métabolites C, N et P à partir de données publiées dans Neff et al. 2001 (tableau 2 et annexe, op. cit.), NGWCLC 2003 (tableau 3.2, op. cit.) et Habib et Guennelon 1983 (eq. 15, page 1 17, op. cit.). Calcul des concentrations initiales de C, N et P (Tableau 2-a) Les concentrations initiales des sols en C, N et P provenant de l'enfouissement des RCS sont exprimés en kg par m 3 selon Poil et al. 2006. Pour obtenir cette concentration initiale du carbone, Ceo, il suffit de diviser les kg de résidus de culture au sol, RCS, par 0,42, soit la teneur approximative en C des RCS ; la teneur initiale en azote des RCS, Cno, est elle égale à Ceo / rCNo. Idem pour ce qui est de la teneur initiale des RCS en P des RCS égale à Ceo / rCPo. Il est important de distinguer entre Cpo et Cpoi, cette dernière valeur étant celle pour les RCS sans , ou avant, l'apport de dPo proportionnel à Cpo avant cet apport (i.e. Cpoi). Formellement, la modélisation iPAZ va donc comporter les calculs suivants ; Ceo Cno Cpoi Cpo L'apport dCpoi permettra d'augmenter pAZB de manière à atteindre un seuil où le P n'est plus limitatif ; les préconisations d'engrais P seront donc fortement proportionnelles à dCpoi. Calcul des concentrations finales de C, N et P (Tableau 2-b) Les équations différentielles de Poil et al. 2006 (éqs. 4 et 5) permettent de s’affranchir de la variable temps. En effet le ratios rCN et rCP (C/N et C/P) calculés pour les résidusphères sont elles fonction des susdites concentrations initiales des substrats C, N et P, Ce, Cn et Cp, une fois atteinte une certaine distance de migration, zc, zn et zc (mm), à partir de l'interface sol : RCS ; C(z) = C ( , exp(rz), sachant que la valeurs r est elle définie comme suit, et cela pour chaque type de substrats, rc, rn et rp (Tableau 2-b) ; et que D est le coefficient de diffusion effective (De), et que k et v sont des coefficients empiriques définis par Poil et al. 2006 (eq. 5). Dans le cas présent, k et v sont spécifiés pour chacun des substrats - métabolites C, N et P (i.e. k cnp et v cnp , selon le cas) comme suit ; kenp ■ De cnp # k/jour (Poil et al. 2008, tableau 2) Vcnp ■ De cnp # mm/j ; Convection (Poil et al. 2008, tableau 2) données publiées dans Neff et al. 2001 (tableau 2 et annexe, op. cit.), NGWCLC 2003 (tableau 3.2, op. cit.) et Habib et Guennelon 1983 (eq. 15, page 1 17, op. cit.). Calcul des concentrations initiales de C, N et P (Tableau 2-a) Les concentrations initiales des sols en C, N et P provenant de l'enfouissement des RCS sont exprimés en kg par m 3 selon Poil et al. 2006. Pour obtenir cette concentration initiale du carbone, Ceo, il suffit de diviser les kg de résidus de culture au sol, RCS, par 0,42, soit la teneur approximative en C des RCS ; la teneur initiale en azote des RCS, Cno, est elle égale à Ceo / rCNo. Idem pour ce qui est de la teneur initiale des RCS en P des RCS égale à Ceo / rCPo. Il est important de distinguer entre Cpo et Cpoi, cette dernière valeur étant celle pour les RCS sans , ou avant, l'apport de dPo proportionnel à Cpo avant cet apport (i.e. Cpoi). Formellement, la modélisation iPAZ va donc comporter les calculs suivants ; Ceo Cno Cpoi Cpo L'apport dCpoi permettra d'augmenter pAZB de manière à atteindre un seuil où le P n'est plus limitatif ; les préconisations d'engrais P seront donc fortement proportionnelles à dCpoi. Calcul des concentrations finales de C, N et P (Tableau 2-b) Les équations différentielles de Poil et al. 2006 (éqs. 4 et 5) permettent de s’affranchir de la variable temps. En effet le ratios rCN et rCP (C/N et C/P) calculés pour les résidusphères sont elles fonction des susdites concentrations initiales des substrats C, N et P, Ce, Cn et Cp, une fois atteinte une certaine distance de migration, zc, zn et zc (mm), à partir de l'interface sol : RCS ; C(z) = C ( , exp(rz), sachant que la valeurs r est elle définie comme suit, et cela pour chaque type de substrats, rc, rn et rp (Tableau 2-b) ; et que D est le coefficient de diffusion effective (De), et que k et v sont des coefficients empiriques définis par Poil et al. 2006 (eq. 5). Dans le cas présent, k et v sont spécifiés pour chacun des substrats - métabolites C, N et P (i.e. k cnp et v cnp , selon le cas) comme suit ; kenp ■ De cnp Vcnp ■ De cnp # k/jour (Poil et al. 2008, tableau 2) # mm/j ; Convection (Poil et al. 2008, tableau 2) sachant que les valeurs A, B, et C, ou a, b et c applicable au substrats C, N et P, respectivement ayant déjà été rapportées ci-dessus au Tableau 3. Tableau 2-b : Modèle iPAZ-R - Code R pour la modélisation du potentiel AZB (pAZB) ; Calculs des concentrations des substrats C. N et P. essentiellement de la cellobiose. de l'ammonium et des phosphates, une fois leurs distances de migrations maximale, z (mmV respectives atteintes. # Habib et Guennelon 1983 éq. 12 # k/jour (Poil et al. 2008, tableau 2) # mm/j ; Cnovection (Poil et al. 2008, tableau 2) # Poil et al. 2006 eq. 5 # kg/m3 ; Poil et al. 2006 eq. 4 # Concentration des substrats C-solubles ... Dec A bdc)/teta)) kc A 2 + kc/Dec) A 0.5 zc Ce # Concentration des substrats N-Nm ... Den A bdn)/teta)) kn A 2 + kn/Den) A 0.5 zn Cn Module diazotrophe ; Tableau 2-c # Concentration des substrats P-P04 ... Dep A bdp)/teta)) kp A 2 + kp/Dep) A 0.5 zp Cp # Habib et Guennelon 1983 éq. 12 # k/jour (Poil et al. 2008, tableau 2) # mm/j ; Cnovection (Poil et al. 2008, tableau 2) # Poil et al. 2006 eq. 5 # kg/m3 ; Poil et al. 2006 eq. 4 # Habib et Guennelon 1983 éq. 12 # k/jour (Poil et al. 2008, tableau 2) # mm/j ; Cnovection (Poil et al. 2008, tableau 2) # Poil et al. 2006 eq. 5 # kg/m3 ; Poil et al. 2006 eq. 4 Il est ainsi possible, toujours dans l'esprit de Poil et al. 2006, de relier le taux de consommation microbienne de ces substrats à leurs taux de diffusion ; plus le taux de diffusion est grand, et inversement plus le taux de consommation est faible, plus chacun des substrats pourra migrer loin de l'interface sol : résidus, soit ici d'approximativement 4 mm pour le cellobiose (C), 2 mm pour l'ammonium (N), et 1 mm pour les phosphates (P). Calcul du la contribution de la diazotrophie non symbiotique (Tableau 2-c) Pour ce faire, je propose un module de la diffusion de l'Nm entre les RCS et RAC (i.e. résidus de racines). À ce titre, j'ai pris en compte le rendement de la culture précédente, et donc la quantité de résidus de culture cellulosique aériens retournée au sol, mais aussi la teneur en N de ces résidus rapporté au rendement agronomique (RDT) de ladite culture précédente. Il s'avère que le RDT, avantageusement d'une céréale ou encore d'un maïs - grain, peut ici servir d'indicateur suffisamment précis du ratio entre les quantités post-récolte de résidus cellulosiques laissés au sol et racinaire, ratio dénommé ici rMS. De plus, le rapport entre RDT et la teneur en N des résidus de culture au sol et non celui des grains, ratio dénommé ici RVN s'avère être lui assez bien corrélé au ratio entre les teneurs en N des résidus aériens et souterrain, soit ici rNS ; rMS = RCS / RAC rNS = %Nrcs / %Nrac RVN = RDT / Nrcs Tableau 2-c : Modèle iPAZ-R - Code R pour la modélisation du potentiel AZB (pAZB) ; Module diazotrophe partie intégrante partie intégrante des calculs des concentrations des substrats C. N et P à z (mm). #CALCUL PE RMS, RNS, ET RCNRDS ET RCNRZS ; LES RDT PERMETTENT EUX P'ÉTABLIR MSRAC ET MSRCS rCNrcs A 0.7056 rNS A -0.2449 MSrcs N1 D A bdn)/teta)) # Habib et Guennelon 1983 éq. 12 for (tbbch in 1 :interculture) { INTERNE RDS. Bimmb A -0.2667) Bminr A -0.2667) IMMB MINR A -0.53)) DIAZO A -0.2493/100)) bimmb Nimmb Ndiazo Nminr N1 Ndiazo N1 Nminr N1 Bimmb A -0.2667) Bminr A -0.2667) IMMB MINR A -0.53)) DIAZO A -0.2493/100)) bimmb Nimmb Ndiazo Nminr N2 Ndiazo N2 Nminr N2 #CALCUL DE LA DYNAMIQUE DN1 /DT TELLE QU'AFFECTÉE PAR LA DIFFUSION Ni 2 ET N21. DIFF12 A 0.5)))) DIFF21 A 0.5)))) n12 Nrds Cn # ka/m3 : Poil et al. 2006 ea. 4 Notons qu'on peut aussi poser RVN = RCS / Nrcs. D'une manière ou l'autre, RVN a pour unités une masse par unité de surface. Or, ce "rendement" aérien fonction des rendements agronomiques et, étonnamment, des richesses en N des résidus cellulosique s'avère être dégressif par rapport à rNS, la richesse relative des RCS par rapport à celles des RAC. Or, plus ce rendement protéique, inversé en en sort par rapport à RDN, est important, plus l'azote résiduelle, i.e. non comprise dans RDN, des RAC sera important par rapport à celles des RCS. Il s'agit donc relation dégressive analogue à celle de l'azote critique (Ne ; %) par rapport à la biomasse totale aérienne (DM ; t/ha) décrite par Justes et al. 1994. Je propose d'établir des relations pour le calcul de rMS, rNS, ainsi que la quantité de résidus cellulosique produite en proportion des celle des racines au moment de la récolte de RDT. Plus précisément, et ce stade sur la base de données d'Allard et al. 2013 ; rMS = 0.0601 *RDT A 0.7056 rNS = 1.7481 * RVN A -0.2449 MSrcs = 1.911*RDT+164.15 iV = 0.0131*RDT+2.3159 | ratio sans unités | ratio sans unités I g / m2 I ratio sans unités Pour fin de validation, il faut aussi estimer les quantités d'azote (minéral ; mg/kg) immobilisé (Nimmb), minéralisé (Nminr) et d'origine diazotrophe (Ndiazo) au sein de la résidusphère (RDS) provenant de l'enfouissement des résidus de culture cellulosique au sol. Plus précisément et tel qu'incorporé au code R (Tableau 2-c) nous calculons ces quantités en fonction du temps (jours), voire avantageusement du temps phénologique (tbbch) ; Nimmb = N1/(1 +exp(-bimmb*tbbch)) Nminr = (Nrcs*MSrcs)+Ndiazo)/(1+exp(-bminr*tbbch))/(8.51*tbbch A -0.53) Ndiazo = 14/(1 +exp(-bdiazo*tbbch)))*(36*N1 A -0.2493/100) I mgN/kg I mgN/kg I mgN/kg Notons que la quantité Ndiazo et intégrée a priori dans l'azote à minéraliser (minéralisable). En effet, la diazotrophie produit avant tout de la biomasse cytosolique ammonifiable ; il est judicieux d'ajouter cette quantité à l'ensemble de l'azote minéralisable provenant des résidus de culture et/ou des biomasses racinaires résiduelles. Les valeurs des coefficients bimmb, bminr et bdizao étant eux directement fonction des rCNrcs, i.e. les ratios C/N des résidus de culture au sol, à savoir ; Bimmb = function(rCNrcs) return(-0,3096*rCNrcs A -0.2667) | coefficient Bminr = function(rCNrcs) return(+0.3096*rCNrcs A -0.2667) | coefficient Afin de tenir compte de la rétroaction négative de la diazotrophie (alias ; activité de réduction de l'acétylène - ARA) lorsque les rNP sont élevés, le coefficient bdiazo est calculé ainsi ; bdiazo Cette fonction est basée sur des données empiriques obtenues de Smith 1992 et Howarth et Marino 1988 et permet de plafonner la contribution de la diazotrophie (Ndiazo ; Figure 3) lors de l'accumulation d'Nm, soit généralement au delà d'une quinzaine de mg-N par kg-sol et/ou l'équivalent d'un mM en solution de sol. Les dynamiques de telles diazotrophies sont exemplifiées à la Figure 5. Plus formellement, cette dynamique de la diazotrophie comporte trois moments ; le calcul de bdiazo selon rNPo sur la base des données de Howarth et Marion 1988 ; > bdiazo = 21 ,569 x exp(-0.4338 x rNPo) | Figure 3, abîme A le calcul de DIAZO selon rNPo sur la base de bdiazo ; > DIAZO = 14 / (1 + exp(-bdiazo x tbbch) | Figure 3, abîme B le calcul de Ndiazo selon rNPo sur la base de DIAZO ; > Ndiazo = DIAZO x 36 x N m (" ° 2493/100) | Figure 3, courbe principale Il en résulte une courbe dégressive d'Ndiazo selon rNPo (Figure 3) fortement atténuée dès rNPo = 20 Pour ce qui est de l'estimation des quantités d'azote (minéral ; mg/kg) immobilisé (Nimmb), minéralisé (Nminr) et d'origine diazotrophe (Ndiazo) au sein de la rhizosphère (RZS), les calculs sont foncièrement les mêmes à l'exception des rCNrcs remplacé par rCNrzs, i.e. les ratios C/N des biomasses souterraines provenant des biomasses racinaires résiduelles, ainsi que les Nrcs et MSrcs remplacées par Nrac et MSrzs (cf. Tableau 2). Nous définissons par la suite les deux compartiments du sol, N1 et N2, contenant les quantités d'azote minéral provenant des susdits processus d'immobilisation, minéralisation et diazotrophie. Nous établissons donc que ; N1 = Nimmb + Nminr | RDS, en mgN/kg (de sol) N2 = Nimmb + Nminr |.RZS, en mgN/kg (de sol) Enfin, il faut tenir compte de la diffusion de l'azote minéral des RAC plutôt riches en N et donc en N minéralisé vers les RCS qui sont elles ici en principe plus pauvres en N. Plus précisément, et encore une fois tel qu'incorporé dans le code R (Tableau 1) ; n12 = (N2/2)*(1-erf(x/(2*(D*tbbch) A 0.5))) | mgN/kg n21 = (N1/2)*(1-erf(x/(2*(D*tbbch) A 0.5))) | mgN/kg Ces quantités d'azote minéral diffusables de N1 (RDS) ver N2 (RZS) sont fonction du coefficient de diffusion de l'ammonium, et non du nitrate (beaucoup plus important) puisque que c'est l'ammonium qui apparaît en premier lors de la minéralisation pendant l'inter-culture hivernale. Voir le Tableau 2-c pour le code R permettant de modéliser ces processus ; le cas échéant, les valeurs proposées pour les différents coefficients sont indiquées. Calcul du potentiel azotobactériens in situ (pAZB) (Tableau 2-d) Le pAZB maxium sera atteint que si l'azote minéral est suffisamment immobilisé, et si le phosphore du sol lui ne l'est pas. En effet, non seulement pAZB est très sujet à rétroaction négative par Nm, mais elle est aussi très énergivore requérant donc d'abondante réserves de ATP et donc aussi de P biodisponible. Rappelon que l'immobilisation de l'azote est donc ici un plus, et celle du P un moins ; c'est pour cela qu'ici l'efficacité relative (er) potentielle de l'activité diazotrophe du sol (AZB) est définie comme suit ; pAZB = pAZBn - pAZBp Pour faire simple, disons que l'pAZB potentielle maximum sera atteint lorsque les taux de formation des métabolites carbonés et azotés seront équilibrés. Par équilibre on veut dire ici que le rapport entre les rCN et P résidusphériques et les rCN et P critiques est approximativement unitaire. Tableau 2-d : Modèle iPAZ-R - Code R pour la modélisation du potentiel AZB (pAZB) ; Calcul du l'efficacité relative de l'azotobactérisation. ou de l'activité de la flore azotobactérienne du sol en cas de non inoculation des RCS, selon les concentration des substrats C, N et P à l'équilibre une fois z (mm) atteints au sens de Poil et al. 2006. # Calcul de l'pAZB selon rCNc et rCPc ... rCN * 1000) rZn Courbe de réponse MIT AZB, coefficient b et calcul de dP (Tableau 6) Le code R proposé (Tableau 2-d) crée un fichier contenant l'ensemble des données ainsi modélisées (Tableau 4), y compris donc pAZB et dCpoi qui peuvent ainsi être mise en relation (Figure 4). Au Tableau 5 sont rapporté les statistique descriptive pour ces résultats de modélisaton. Cette relation est clairement mitschlerichienne et dites ici MITAZB. Donc, à l'aide de la fonction (procédure) nls2 de R (Tableau 6), il est en principe possible d'établir les coefficients A, b et c d'une relation en touts points analogue à l'équation (1 ) ci-dessus, à savoir maintenant ; pAZB = A · [1 - exp(-c · (b ■ Cpoi + dCpoi))] (3) Notons que les teneurs P du sol sont ici remplacés par les concentration en P du sol provenant directement des résidus de culture enfouis dans le sol (Cpoi) et les apports d'engrais P, dP, par l'équivalent dCpoi proportionnel à une certaine fraction de Cpoi, soit ici de 1 à 25%. La réalisation de l'invention va maintenant dépendre de la combinaison de Cpoi et dCpoi. Or, ici cette équation est en ce sens sans solution, seule le coefficient b s'étant stabilisé (après 10 000 itérations) de l'algorithme nls2 via le code R proposé au Tableau 6. Cependant, il s'avère que cela suffira ici à calculer les dCpoi en équivalents kg-P par hectare (ha ; Figure 5). Tableau 4 : Résultats de la modélisation iPAZ-R comprenant une préconisation de dP dCpoi/ha, i.e. en fonction d'une fraction de l'apport de P initialement retournée au sol via les RCS et de la réponse MITAZB (éq. 3), y compris et plus particulièrement le coefficient b (éq. 4) permettant ainsi de calibrer cette dP de manière à obtenir une pAZB optimale. Ici, et selon le code R proposé au Tableau 5, le coefficient b est estimé à -0,1990 ; le graphe de cette fonction dCpoi = -b ■ Cpoi est rapporté à la Figure 5. Etc. N = 1 500 Préconisation des dP selon l'apport initiale de P via RCS (Figure 5) Pour obtenir les dCpoi nécessaire à la pleine réalisation d'pAZB il faut résoudre l'équation (3) lorsque pAZB = 0,00. Or, puisque pAZB doit ici égaler 0,00, dCpoi n'est plus que fonction du coefficient b, à savoir ; dP = -b ■ Cpoi (4 ; mg-P / kg-sol, ou kg-P / ha) Ce résultat surprenant provient du fait que le logarithme naturel (LN) de [0,00/A] que l'on rencontre lors du réarrangement le l’équation (3) est égale à 1 , et que le terme [0/-c] est lui défini comme étant égale à 0,00. In fine, y = -bx, ou ici dCpoi = -bCpoi. Il est ainsi possible de rétro - calculer dCpoi en fonction de Cpoi (Figure 5), ce qui s'avère être très pratique d'un point de vue agronomique. Tableau 5 : Données modélisées par iPAZ-R ( alias output data) - pAZB iPAZ.xls Echantillon Unités Minimum Maximum Médiane Moyenne 1-1) RCS, Ceo, Cno, Cpo, rCNo, rCPo, rNPo, Ce, Cn, Cp, pAZBn, pAZBp, pAZB, Cpoi dCpoi Cpoi/ha dCpoi /ha kg/ha kg/m 3 kg/m 3 kg/m 3 sans unités sans unités sans unités kg/m 3 kg/m 3 kg/m 3 Sans unités Sans unités Sans unités kg/m 3 sans unités kg/ha kg/ha 3947 1 ,658 0,013 0,001 33,000 586 17,759 2,686 0,030 0,001 0,899 0,944 -0,454 0,001 0,000 0,754 0,150 9000 3,780 0,115 0,008 125,000 2747 21,977 11,252 0,360 0,024 1,315 1,503 0,106 0,006 0,002 8,062 1,604 6421 2,697 0,035 0,002 75,000 1519 20,252 5,890 0,097 0,005 1 ,090 1 ,208 -0,117 0,002 0,000 2,198 0,437 6360 2,671 0,044 0,003 76,600 1584 20,035 6,696 0,121 0,006 1 ,099 1 ,213 -0,115 0,002 0,000 2,841 0,565 1791 0,752 0,026 0,002 33,181 781 1 ,516 3,043 0,080 0,005 0,114 0,138 0,111 0,001 0,000 1 ,873 0,373 Moyenne géom. 6094 2,559 0,037 0,002 68,848 1375 19.977 5.977 0,099 0,005 1 ,093 1 ,205 0,002 0,000 2,326 0,463 Tableau 6 : Proposition de code R pour la détermination des coefficients A, b et c rm(list=ls()) data Cpoi 1 5] kg/m3 nie * (1 - exp(-c*(b * Cpoi+dCpoi))), start = list(A=0.05, b=-.25, c=2500), control = nie, trace = TRUE) summary(datal) APPLICATIONS BIOINDUSTRIELLES ET AGRONOMIQUES L'agronome, l'agriculteur ou un quelconque utilisateur pour réaliser l'invention en déterminant la quantité de RCS. Pour ce faire, il peut le faire directement via la méthode dite LPTM (eg. Shelton et Jasa 2009) ou encore indirectement à partir du RDT de la culture précédente (eg. Steiner et la. 1994a, 1994b, 1995, et Steiner 1994) ; in fine, il peut vérifier si ces deux types de mesures concordent (cf. Steiner et al. 2000). Le dosage du P dans un échantillon représentatif de ces RCS servira à calculer Cpoi, et donc dP (dCpoi) selon l'équation (4) ci dessus, ou encore à partir de la Figure 5. Notons qu'à la Figure 5, dP est rapporté en unités de P, et non de P205 comme c'est souvent le cas pour des engrais phosphoriques commerciaux. Notons aussi que les dP ainsi préconisés sont non seulement bien inférieurs, "rationalisés" en sorte, mais aussi proportionnels à Cpoi. Cela détonnera auprès de l'homme du métier puisque normalement les retours de P sont comptabilisés entant que tels et contribuent généralement à réduire dP, et non à l'augmenter. Il existe dans le commerce de nombreuses préparation d'engrais P soluble dans l'eau pulvérisables, le plus souvent destinées à l'application foliaire ; leur application directement RCS est elle aussi ici assez inhabituelle - au Tableau 7 une liste de quelques une de ces produits ; Tableau 7 : Quelques exemples d'engais phosphoriques commercialisés pouvant être utilisés dans la réalisation de la présente invention. Produit Fabricant g-P205/L g-P/L Formulation Utilisation habituelle Capfol™ Florilège™ Ultra Solustar™ Mg Solustar™ PZ Starmax™ Mg Agronutrition, France Agronutrition, France Agronutrition, France Agronutrition, France Agronutrition, France 60 250 440 440 440 26,2 109 192 192 192 liquide foliaire liquide foliaire liquide sol (démarrage) liquide sol (démarrage) liquide foliaire YaraVita™ Magphos K YaraVita™ Pentaflo YaraVita™ Seniphos YaraVita™ Solatrel YaraVita™ Zeatrel Yara, France Yara, France Yara, France Yara, France Yara, France 440 240 310 440 440 192 liquide 105 liquide 135 liquide 192 liquide 192 liquide foliaire foliaire foliaire (fruits) foliaire foliaire Granoligo™ PB 29 Granoligo™ PZ 50 Microstar™ Céréales Agronutrition, France Agronutrition, France Agronutrition, France (14% p/p) 6,12% (16% p/p) 6,99% (47% p/p) 20,5% granulés granulés granulés sol sol sol Ces préparations n'ont pas été expressément conçues pour application sur RCS ; dans bien des cas elles ont probablement trop raffinées et coûteuses (de la "confiture aux cochons" en sorte). Elles démontrent néanmoins que ce type de préparations peuvent être fabriquées, quitte à en retirer les oligoéléments (sauf le Mo s'il y a), les adjuvants foliares, etc. Notons aussi les quelques exemples de formulations granulées (solides) ; suffisamment fines, pulvérulentes en sorte, elles peuvent aussi être répandues (saupoudrées) sur les RCS, bien que cette techniques nécessitera quelques mises au point afin d'éviter les problèmes sanitaires associés à l'utilisation pleins champs de telles poudre. A suivre. Les bouillies ainsi constituées peuvent être appliquées avec ou après des inocula azotobactériens capable, en soi, d'augmenter le potentiel AZB des sol ; ce type de valorisation bactériologique des RCS est maintenant en plein développement (eg. Claude et Fillon 2004, EP2845906, EP2730926, EP2728353, EP2345319, etc.). Si l'on devait combiner en une seule bouillie, inoculum azotobactérien et P soluble, la pression osmotique agissant sur les bactéries ainsi réintroduites aux RCS pourra être avantageusement gérée selon la technique proposée dans EP2647695 à savoir un mode de préparation de bouillies inoculantes pulvérisables comprenant des populations bactériennes moyennement halophiles. Réferences Allison, Steven D., Peter M. Vitousek. 2005. Responses of extracellular enzymes to simple and complex nutrient inputs. Soil Biology & Biochemistry 37 : 937-944 Augusto, Laurent, Florian Delerue, Anne Gallet-Budynek et David L. Acha. 2013. Global assessment of limitation to symbiotic nitrogen fixation by phosphorus availability in terrestrial ecosystems using a metaanalysis approach. Global Biogeochemical Cycles 27 : 804-815, doi:10.1002/gbc.20069 Binkley, Dan, Randy Senock , Kermit Cromack Jr. 2003. Phosphorus limitation on nitrogen fixation by Facaltaria seedlings. Forest Ecology and Management 186 : 171-176 Claude, PP et L Fillion. 2004. Effet de l’apport d’un inoculum bactérien aux résidus de culture de maïs-grain au sol sur le rendement et la qualité de blés d’hiver panifiables en France. Agrosol. 15 (1) : 23-29. Colomb, B, P. Debaeke, C. Jouany et J.M. Nolot. 2007. Phosphorus Management In Low Input Stockless Cropping Systems : Crop and Soil Responses to Contrasting P Régimes In A 36-Year Experiment In Southern France. Europ. J. Agronomy 26 (2007) 154-165 COMIFER 2011. Calcul de la fertilisation azotée Guide méthodologique pour l’établissement des prescriptions locales Cultures annuelles et prairies Édition 2011. Brochure éditée par le Comifer, Le Diamant A, 92909 Paris La Défense Cedex, ISBN 978-2-910393-09-0 COMIFER 2007. Fertilisation PK : grille de calcul de dose. Le groupe PKMg du COMIFER 2007 Damon, PM, B Bowden T Rose et Z Rengel. 2014. Crop residue contribution to phosphorus pools in agricultural soils : a review. Soil Biol. Biochem. 74 : 127-137 Eisele, KA, DS Schimel, LA Kapustka et WJ Parton. 1989. Effects of available P an N:P ratios on nonsymbiotic dinitrogen fixation in tallgrass prairie soils. Oecologia 79 : 471-474 FAO 2008. Efficiency of soil and fertilizer phosphorus use. FAO Fertilizer and Plant Nutrition Bulletin no. 18. Food and agriculture organization of the United Nations (FAO, Rome), 2008 Gupta, WSR, MM Roper et DK Roget. 2006. Potential for non symbiotic N2 fixation in different agroecological zones of Southern Australia. Aust. J. Soil Res. 44 : 343-354 Habib, Rober et Roger Guennelon. 1983. Evaluation du coefficient de diffusion du phosphore en sol de limon argileux calcaire. Agronomie 3 (2), pp.113-121. Halsall DM et AH Gibson. 1986. comparison of two Cellulomonas strains and their interaction with Azospirillum brasilense in dégradation of wheat straw and associated nitrogen fixation. Appl. Environ. Microbiol. 51 : 855-861 Howarth, Robert W. et Roxanne Marino. 1988. Nitrogen fixation in freshwater, estuarine, and marine ecosystems. 2. Riogeochemical control9. Limnol. Oceanogr., 33(4, part 2) : 688-701 Idriss, EE, O Makarewicz, A Farouk, K Rosner, R Greiner, H Bochow, T Richter et R Borriss. 2002. Extracellular phytase activity of Bacillus amyloliquefaciens FZB45 contributes to its plant-growth-promoting effect. Microbiology 148, 2097-2109 Manzoni, Stefano, John A. Trofymow, Robert B. Jackson et Amilcare Porporato. 2010. Stoichiometric contrais on carbon, nitrogen, and phosphorus dynamics in decomposing litter. Ecological Monographs, 80(1), pp. 89-106 McGill, W. B. et E. K. Christi. 1983. Biogeochemical Aspects of Nutrient Cycle Interactions in Soils and Organisms. In SCOPE 21 - The Major Biogeochemical Cycles and Their Interactions , chapitre 9 / http://www.scopenvironment.org/downloadpubs/scope21/chapter09.html Morel, C, C Plenchette et JC Fardeau. 1992. La fertilisation phosphatée raisonnée de la culture du blé. Agronomie 12 : 565-579 Neff, JC et GP Asner. 2001. Dissolved organic carbon in terrestrial ecosystems : synthesis and a model. Ecosystems 4 : 29-48. NGWCLC 2003. Review of ammonium atténuation in soil and groundwater. National Ground and Contaminated Land Centre (UK), report NC/02/49 Poil, C. J. Ingwersen ,M. Stemmer ,M.H. Gerzabek et E. Kandeler. 2006. Mechanisms of solut é transport affect small-scale abundance and functiono f soil microorganisms in the detritusphere. European Journal of Soil Science, August 2006, 57, 583-595 / doi: 10.1111/j.1365-2389.2006.00835.x Reed, Sasha C., Cory C. Cleveland et Alan R. Townsend. 2013. Relationships among phosphorus, molybdenum and free-living nitrogen fixation in tropical rain forests: results from observational and experimental analyses. Biogeochemistry DOI 10.1007/S10533-013-9835-3 Reed, Sasha C., Cory C. Cleveland, et Alan R. Townsend. 2011. Functional Ecology of Free-Living Nitrogen Fixation: A Contemporary Perspective. Annu. Rev. Ecol. Evol. Syst. 2011. 42:489-512 Reed, Sasha C., Timothy R. Seastedt , Colin M. Mann , Katharine N. Suding , Alan R. Townsend , Karie L. Cherwin. 2007. Phosphorus fertilization stimul â tes nitrogen fixation and increases inorganic nitrogen concentrations in a restored prairie, applied soil ecology 36 : 238-242 Rose TJ, L Liu et M Wissuwa. 2013. Improving phosphorus efficiency in cereal crops : is breeeding for reduced grains phosphorus concentration part of the solution ? Frontiers in Plant Science / November 2013 / Vol. 4 / Article 444. Sattari SZ, AF Bouwman, KE Giller et MK van Ittersum. 2012. Residual soil phosphorus as the missing piece in the global phosphorus crisis puzzle. PNAS 109 (16) : 6348-6353 Shelton DP et J Jasa 2009. Estimating percent residue cover using the line-transect method. NebGuide, University of Nebraska-Lincoln Extension, Institute of Agriculture and Natural Resources. G1931/ February 2009 / www.extension.unl.edu/publications Smith, Val H. 1992. Effects of Nitrogen: Phosphorus Supply Ratios on Nitrogen Fixation in Agricultural and Pastoral Ecosystems. Biogeochemistry, Vol. 18, No. 1. (1992), pp. 19-35. Steiner JL, HH Schomberg, PW Unger et J Cresap. 2000. Biomass and Residue Cover Relationships of Fresh and Decomposing Small Grain Residue. Soil Sci. Soc. Am. J. 64:2109-2114 Steiner, JL, HH Schomberg et JE Morrison. 1994a. Measuring surface residue and calculating losses from décomposition and redistribution, pp. 21-29 In W. C. Moldenhauer (Managing Ed.) and B. A. Stewart (Régional Ed.) Crop Residue Management to Reduce Erosion and Improve Soil Quality: Southern Great Plains. USDA-ARS Conservation Research Report No. 37. Steiner, JL, HH Schomberg et PW Unger. 1995. Residue Décomposition Submodel. pp. D1-D12, In L. J. Hagen, L. E. Wagner, and J. Tatarko. 1995. Wind Erosion Prédiction System (WEPS). Technical Documentation. Beta Release 95-08. USDA, Agricultural Research Service, Wind Erosion Research Unit, Throckmorton Hall, Manhattan, Kansas, 66506. Steiner, JL. 1994. Crop residue effects on water conservation, pp. 41-76, In P. W. Unger (ed.) Managing Agricultural Residues. Lewis Publ., Boca Raton, FL. Steiner, JL., HH Schomberg, CL Douglas Jr. et AL Black. 1994b. Standing stem persistence in no-tillage small grain fields. Agron. J. 86:76-81. Thibaud MC, C Morel et JC Fardeau. 1988. Contribution of phosphorus issued form crop residues to plant nutrition. Soil Sci. Plant Nutr. 34 : 481-491 Ussiri, D. A., P. N. S. Mnkeni , A. F. MacKenzie et J. M. R. Seraoka 1998. Soil test calibration studies for formulation of phosphorus fertilizer recommendations for maize in Morogoro district, Tanzania. II. estimation of optimum fertilizer rates, Communications in Soil Science and Plant Analysis, 29 (17-18) : 2815-2828 Vitousek, Peter M. et Robert W. Howarth. 1991. Nitrogen Limitation on Land and in the Sea: How Can It Occur? Biogeochemistry, Vol. 13, No. 2 (1991), pp. 87-115 Wang, Y.-P., B. Z. Houlton, and C. B. Field. 2007. A model of biogeochemical cycles of carbon, nitrogen, and phosphorus including symbiotic nitrogen fixation and phosphatase production. Global Biogeochemical Cycles, VOL. 21 , GB1018, doi:10.1029/2006GB002797, 2007 0 5 0 REVENDICATIONS 1. Procédé pour la calibration des apports aux cultures agronomiques d'engrais phosphoriques comprenant l'établissement d'une courbe de réponse mitscherlichienne décrivant l'activité diazotrophe de la flore azotobactérienne dans un sol arable en proximité de résidus de la culture précédente enfouis, courbe de réponse dite ici de type MITAZB, et cela en fonction, (i) des teneurs en C, P, voire N, desdit résidus de culture cellulosiques, avantageusement pailleux, au sol (RCS) avant leur enfouissement en pré - semis d'une culture agronomique suivante et, (ii) d'un équivalent de phosphore exogène, dit ici dose P ou dP, proportionnel aux susdites teneurs en C et P des résidus de culture au sol, caractérisé en ce qu'il nécessite ; > la détermination de la quantité desdits résidus de culture au sol (qRCS), soit indirectement via une relation empirique fonction du rendement de la culture précédente (RDT), soit plus directement au sol via la détermination du degré de couverture du sol par ces résidus (cRCS) ; la détermination des teneurs en N (azote) et en P (phosphore) desdits RCS de manière à obtenir Cpoi, i.e. la quantité de P initialement retournée au sol par ces RCS > le calcul de dP en fonction de Cpoi capable d'assurer que le niveau d'activité diazoptrophe en proximité desdits résidus de culture une fois enfouis (pAZB) maximale. 2. Procédé pour la calibration des apports aux cultures agronomiques d'engrais phosphoriques selon la revendication précédente caractérisé en ce que le niveau d'activité diazoptrophe en proximité desdits résidus de culture une fois enfouis (pAZB) est défini comme l'équilibre entre l'immobilisation nette de l'azote du sol, phénomène favorable à ladite activité diazotrophe, et l'immobilisation nette du phosphore du sol, phénomène défavorable à ladite activité diazotrophe. 3. Procédé pour la calibration des apports aux cultures agronomiques d'engrais phosphoriques selon la revendication précédente caractérisé en ce que l'équilibre entre l'immobilisation nette de l'azote du sol, phénomène favorable à ladite activité diazotrophe, et l'immobilisation nette du phosphore du sol, phénomène défavorable à ladite activité diazotrophe, est défini formellement comme le rapport entre pAZBn et pAZBp (pAZBn / pAZBp), et plus avantageusement comme la différence entre pAZBn et pAZBp (pAZB - pAZB), pAZB -n et -p étant, respectivement les ratios rCN/rCNc et rCP/rCPc, i.e. les rapports entre ratios C/N (ou C/P) de la résidusphère à l'équilibre et C/N (ou C/P) critique en dessus du quel l'azote (ou le phosphore) est libéré sous forme minéral (minéralisation). 4. Procédé pour la calibration des apports aux cultures agronomiques d'engrais phosphoriques selon une quelconque des revendications précédentes caractérisé en ce que le calcul de dP en fonction de Cpoi et d'un coefficient b issue d'une relation MITAZB de type pAZB = A ■ [1 - exp(-c ■ (b · Cpoi + dCpoi))] 0 lorsque dCpoi est un équivalent du P apporté aux RCS et/ou au sol affectant pAZB, dont la valeur est comprise entre 0,10 et 0,50, plus particulièrement entre 0,15 et 0,30, et avantageusement ici d'environs 0,20, peut être effectué comme suite ; dP = -b x Cpoi dP pouvant ainsi être exprimé en kg-P par m3 de sol, mg-P / kg-sol, ou encore en kg-P par hectare une fois les conversions unitaires nécessaire effectuées. 5. Procédé pour la calibration des apports aux cultures agronomiques d’engrais phosphoriques selon une quelconque des revendications précédentes caractérisé en ce que la détermination de la quantité des résidus de culture cellulosiques au sol (qRCS ; g/m 2 ) est déterminée indirectement à partir du rendement de la culture précédente (RDT ; g/m 2 ) selon une des relations empiriques suivantes ; qRCS B TH = [0,911 x RDT + 164] qRCSmaïs = [1 ,51 x RDT + 3,7 - RDT] qRCS S orgho = [2,56 x RDT + 0.4 - RDT] qRCSoéréaies = [2,13 x RDT + 1 ,1 - RDT] ou encore plus directement via la détermination préalable du degré de couverture du sol (cRCS ; %), la quantité de RCS (qRCS ; t/ha) pouvant être déterminé par l’une des relation empiriques suivantes ; qRCS B TH = 0,0633 x exp(0,0493 x cRCS) | R2 = 0,9275 qRCS M a(s = 0,0519 x exp (0,0554 x cRCS) | R2 = 0,8647 qRCSrournesoi = 0,0015 x cRCS 2 - 0,003 x cRCS + 0,378 | R2 = 0,9928 6. Procédé pour la calibration des apports aux cultures agronomiques d’engrais phosphoriques selon une quelconque des revendications précédentes caractérisé en ce que les cultures agronomiques sont des grandes cultures agronomiques d’hivers choisies parmi un groupe de cultures comprenant les céréales, y compris le blé, l’avoine, le seigle, le triticale, l’orge, le colza, le lin, ou encore toute autre grandes culture agronomique pouvant être semis en présence de résidus de culture cellulosiques, avantageusement pailleux au sol. 7. Méthode de fertilisation phosphorique pour cultures agronomiques comprenant l’application d’une dose dP de phosphore calibrée selon une quelconque des revendications précédentes caractérisée en ce que dP est sous forme soluble est appliquée par pulvérisation liquide directement aux RCS avant leur enfouissement, avantageusement en pré - semis d’une grande culture agronomique. 8. Méthode de fertilisation phosphorique pour cultures agronomiques comprenant l’application d’une dose dP de phosphore selon la revendication 7 caractérisée en ce que la dP sous forme pulvérulente est appliquée par pulvérisation solide directement aux résidus de culture cellulosiques au sol avant leur enfouissement, avantageusement en pré - semis d’une grande culture agronomique telle que le blé tendre, le colza, ou encore toute autre grandes cultures agronomiques pouvant être semis en présence de résidus de culture cellulosiques, avantageusement pailleux, au sol. 9. Méthode de fertilisation phosphorique pour cultures agronomiques comprenant l'application d'une dose dP de phosphore selon les deux revendications précédentes caractérisée en ce que les résidus de culture cellulosique au sol sont bactérisés à l'aide d'une préparation bactérienne réputée capable d'augmenter l'activité diazotrophe, voire l'activité de réduction de l'acétylène (ARA), du sol en proximité desdits résidus de culture un fois enfouis. 10. Méthode de fertilisation phosphorique pour cultures agronomiques comprenant l'application d'une dose dP de phosphore selon la revendication précédente caractérisée en ce que la préparation bactérienne contient des biomasses bactériennes capables de minéraliser partiellement en P inorganique le P organique, notamment sous forme de sels d'acide phytique. 11. Méthode de fertilisation phosphorique pour cultures agronomiques comprenant l'application d’une dose dP de phosphore selon les deux revendications précédentes caractérisée en ce que la préparation bactérienne est combinée avec dP au sein d'une même bouillie pulvérisable directement sur lesdits résidus de culture au sol.