FR 3027187 A1 20160422 FR 1459945 A 20141016 - 1 - Procédé de détermination de la quantité d'azote à fournir à une culture La présente invention concerne un procédé de détermination de la quantité d'azote à fournir à une culture végétale. Le colza est une plante annuelle à fleurs jaunes de la famille des Brassicacées, famille anciennement nommée Crucifères, qui est largement cultivée pour la production d'huile alimentaire de tourteaux et d'agro/bio carburants. En France, le colza d'hiver se sème en fin d'été (du 15 août au 15 septembre), le peuplement recherché pour la culture étant d'environ 30 pieds/m2 pour les lignées et 20 pieds/m2 pour les hybrides. Le colza peut présenter une croissance très forte à l'automne, beaucoup plus que les céréales, et absorber jusqu'à plus de 250 kgN/ha à cette période suivant la date de semis et la disponibilité en azote du sol.. Le colza a un besoin élevé en azote (environ 7 kg/q).. La fertilisation totale du colza doit apporter entre 0 et 250 unités d'azote avec des apports d'engrais chimiques ou 15 organiques. Un colza semé tôt peut mobiliser des quantités importantes d'azote à l'automne et durant l'hiver, c'est-à-dire pendant la période où les risques de pertes de nitrate par lixiviation sont les plus élevés. Il va donc contribuer à appauvrir le sol en azote minéral, donc à réduire ces pertes 20 Aussi, au vu de ce fort besoin en azote et de l'impact environnemental que peut générer l'ajout d'engrais azoté aux cultures, il est nécessaire de disposer de modes de calcul précis des quantités d'azote à apporter aux cultures, notamment au printemps. En France, en zones vulnérables vis-à-vis de la pollution par les nitrates, le mode de calcul de la quantité d'engrais azoté autorisé est défini par arrêtés préfectoraux 25 régionaux. Ces arrêtés donnent des préconisations spécifiques en fonction des conditions de culture (sol, climat...). Les méthodes préconisées par lesdits arrêtés se basent généralement sur les préconisations du comité français d'études et de développement de la fertilisation raisonnée (COMIFER), reposant sur le bilan de masse prévisionnel et ses adaptations. 30 Les préconisations du COMIFER sont relativement complexes car très complètes et sont difficiles à mettre en oeuvre par l'exploitant. Aussi, il existe un besoin de fournir une méthode simple permettant à l'exploitant de calculer la dose la plus juste d'engrais azoté à fournir à une culture d'azote au printemps. 35 Afin de pallier ces inconvénients, la demanderesse a proposé en 1998 une règle de calcul très simplifiée de la quantité d'azote à fournir à une culture de colza, en tenant compte de la nature des sols, des apports azotés sous forme organique avant la - 2 - culture, et du rendement espéré pour la culture. Si cette méthode a constitué un progrès considérable à la fin des années 1990 et au cours de années 2000, elle nécessite toutefois d'être améliorée grâce au progrès de la connaissance depuis lors, pour prendre en compte de façon plus précise des situations culturales particulières et limiter les hétérogénéités de paramétrage entre régions puisque la méthode de 1998 distinguait 9 régions. Aussi, il existe un réel besoin d'améliorer l'aide à la détermination de la quantité d'azote à ajouter à une culture de colza, qui soit plus précise, plus uniforme sur le territoire français, et qui permette toujours à l'exploitant de déterminer la dose la plus 10 juste en respectant l'environnement et les contraintes réglementaires. Un des buts de l'invention est donc de fournir une méthode simplifiée permettant à un exploitant de fertiliser une culture de colza. Un autre but de l'invention est de fournir des moyens de mise en oeuvre de cette méthode. 15 Ainsi, l'invention concerne une méthode de fertilisation d'une culture de colza d'hiver en France métropolitaine à l'aide d'une quantité optimale d'engrais azoté, au printemps, dans le but de stimuler la croissance de ladite culture de colza, ladite méthode comprenant les étapes suivantes : a) une étape de mesure d'au moins la masse de matière fraiche aérienne à la sortie 20 de l'hiver MVSH d'au moins deux échantillons de ladite culture, éventuellement la masse de matière verte aérienne à la sortie de l'hiver MVSH et la masse de matière verte aérienne à l'entrée de l'hiver MVEH, b) une étape de collecte des paramètres suivants : - la localisation géographie L de ladite culture, 25 - la nature N du sol sur lequel se trouve ladite culture, - les pratiques culturales suivantes ayant un effet sur la fourniture d'azote à ladite culture: - la présence sur la parcelle de pois protéagineux pendant la campagne précédent ladite culture, 30 - la co-culture de ladite culture de colza avec des plants de légumineuses gélifs, - les habitudes d'apport de produits résiduaires organiques sur le sol où se trouve ladite culture pendant les années précédent ladite culture, 35 - l'apport de produits organiques pendant l'été précédent la culture et sur ladite culture pendant l'été et l'automne, et - 3 - l'objectif de rendement y pour ladite culture, c) une étape de calcul, notamment par un produit programme d'ordinateur mis en oeuvre sur un support approprié, de la quantité d'engrais azoté à ajouter X à ladite culture selon la méthode des bilans prévisionnels simplifiée et adaptée au colza, de sorte que : (i) si la culture de colza est localisée dans une région de France métropolitaine qui n'est pas la région Nord-Est ou la région Auvergne-Rhône-Alpes, la quantité de d'azote à ajouter X, au printemps, à ladite culture, exprimée en kilogramme d'azote par hectare (ou en unité d'azote), est calculée selon la formule suivante ; X = [(b.y) + Rf] - (Pi + Ri + M + Mha + Mpro1 + Fleg + Fass), où, b correspond au besoin unitaire en kg d'azote absorbé à la fermeture du bilan sur la quantité de graines produites par hectare, y est l'objectif de rendement exprimé en quintaux par hectare aux normes de commercialisation, Pi est la quantité d'azote absorbé par la culture à l'ouverture du bilan (kgN/ha), tel que si le produit MVSH x Coeff SH est supérieur ou égal à MVEH x Coeff EH, alors Pi= MVSH x Coeff SH, sinon, Pi = (MVSH x Coeff SH) + 1/2((MVEH x Coeff EH) - (MVSH x CoeffSH) / 1,35) où Coeff EH est le coefficient de conversion de la biomasse verte aérienne en kg/m2 en quantité d'azote absorbé en kg d'azote par hectare à l'entrée de l'hiver, et Coeff SH est le coefficient de conversion de la biomasse verte aérienne en kg/m2 en quantité d'azote absorbé kg d'azote par hectare à la sortie de l'hiver, et où Rf est la quantité d'azote minéral dans le sol à la fermeture du bilan exprimée en kg d'azote par hectare, Ri correspond à la quantité d'azote minéral dans le sol à l'ouverture du bilan exprimée en kg d'azote par hectare, M correspond à la minéralisation nette de l'azote de l'humus et des résidus de récolte du précédent exprimée en kg d'azote par hectare, Mpro1 est la minéralisation nette de l'azote organique des PRO apporté l'année de la culture avant l'ouverture du bilan exprimée en kg d'azote par hectare, Mpro1=qPRO x NPRO x Keq - 4 - où qPRO correspond à la quantité de produit résiduaire organique épandu, NPRO correspond à la teneur en azote desdits produits résiduaires, et Keq correspond au coefficient d'équivalence engrais Mha correspond à la fourniture d'azote liée à l'historique d'apport de produits résiduaires organiques Fleg correspond au supplément de fourniture d'azote lié à la présence sur la parcelle de pois protéagineux à la campagne précédent la culture de colza exprimée en kg d'azote par hectare, Fass correspond au supplément de fourniture d'azote lié à des cultures compagnes de légumineuses gélives exprimé en kg d'azote par hectare, où N est un sol superficiel (Nss), un sol profond (Nsp) ou un sol de craie de champagne (Nsc), et Rf, Ri et M sont dépendants de N et de L, (ii) si la culture de colza est localisée dans la région Nord-Est ou dans la région Auvergne-Rhône-Alpes, la quantité de d'azote à ajouter X, au printemps, à ladite culture, exprimée en unité d'azote, est calculée selon la formule suivante : X = (((b.y) - (Pi + Np)) / CAU) - (Mha + Mpro1+ Fleg + Fass), où, b correspond au besoin unitaire en kg d'azote absorbé à la fermeture du bilan sur la quantité de graines produites par hectare, y est l'objectif de rendement (q/ha aux normes), Pi est la quantité d'azote absorbé par la culture à l'ouverture du bilan (kgN/ha), tel que si le produit MVSH x Coeff SH est supérieur ou égal à MVEH x Coeff EH, alors Pi= MVSH x Coeff SH, sinon, Pi = (MVSH x Coeff SH) + 1/2((MVEH x Coeff EH) - (MVSH x CoeffSH) / 1.35) où Coeff EH est le coefficient de conversion de la biomasse verte aérienne (kg/m2) en quantité d'azote absorbé (kgN/ha) à l'entrée de l'hiver, et Coeff SH est le coefficient de conversion de la biomasse verte aérienne (kg/m2) en quantité d'azote absorbé (kgN/ha) à la sortie de l'hiver, et où Np (= P0: fourniture d'azote par le sol) correspond la fourniture - 5 - d'azote par le sol pendant la période d'ouverture du bilan en kg d'azote par hectare, Mpro1 est la minéralisation nette de l'azote organique des PRO apporté l'année de la culture avant l'ouverture du bilan exprimée en kg d'azote par hectare, Mpro1=qPRO x NPRO x Keq où qPRO correspond à la quantité de produit résiduaire organique épandu, NPRO correspond à la teneur en azote desdits produits résiduaires, et 10 Keq correspond au coefficient d'équivalence engrais Mha correspond à la fourniture d'azote liée à l'historique d'apport de produits résiduaires organiques, Fleg correspond au supplément d'absorption d'azote par la culture exprimée en kg d'azote par hectare, ledit supplément étant lié à la 15 présence sur la parcelle de pois protéagineux lors de la campagne précédent la culture de colza, Fass correspond au supplément de fourniture d'azote par la culture lié à des cultures compagnes de légumineuses gélives exprimé en kg d'azote par hectare, 20 CAU est le coefficient apparent d'utilisation variant de 0,6 à 1, où N est un sol superficiel (Nss) ou un sol profond (Nss) et Np sont dépendants de N et L, d) une étape de correction du calcul tel que : si le produit (b.y) calculé est supérieur à 330, alors le produit (b.y) est 25 fixé à 330, si X est supérieur à 250 unités d'azote, alors X est fixé à 250, b variant de 6 à 7, e) et éventuellement une étape de pesées des engrais azotés apportés correspondant au total à la valeur X obtenue si celle-ci est non nulle. 30 L'invention est basée sur la constatation faite par les inventeurs que les nouvelles règles de calcul permettent d'obtenir de manière simplifiée une estimation de la quantité d'engrais azoté à une culture de colza qui respecte les l'esprit de la méthode préconisée par le COMIFER, et pourrait à terme servir de base dans le cadre des arrêtés préfectoraux définissant la réglementation de 35 l'azote dans les différents département de France métropolitaine. Cette nouvelle méthode de calcul de la quantité d'azote a ajouter, sous forme d'engrais à une culture de colza, à partir de la reprise de la végétation au printemps est - 6 - désormais basée sur un calcul simplifié, sans trop de spécificités régionales, et tient compte des fournitures d'azote par le sol lié au système de culture et des apports azotés sous forme organique fournis à la culture de colza avant l'hiver. La méthode de l'invention est simple. L'exploitant qui souhaite mettre en oeuvre la méthode de l'invention doit donc dans un premier temps au moins déterminer la masse de matière fraiche aérienne de la culture de colza à la sortie de l'hiver. Cela signifie que l'exploitant, qui a semé sa culture de colza au mois d'août de l'année n, laisse passer l'hiver, et à la reprise de la végétation à la fin de l'hiver l'année n+1, va déterminer la masse de matière fraiche aérienne. Par « la masse de matière fraiche aérienne » on entend la masse végétale fraiche, c'est-à-dire non sèche (vivante), de colza correspondant à toute la partie verte de la plante qui n'est pas souterraine : la tige, les feuilles (hors feuilles sénescentes non vertes),. La mesure de cette masse de matière fraiche aérienne est réalisable selon plusieurs techniques décrites dans l'exemple 2 ci-après. Il est préconisé de réaliser une véritable pesée, en coupant les plants de colza au niveau du collet au ras du sol sur plusieurs échantillons de la parcelle. Avantageusement ces échantillons ou placettes ont une dimension d'environ 1 m2 précisément mesuré. Avantageusement, l'exploitant prélèvera la totalité des parties ariennes du colza sur au moins deux placettes, avantageusement sur 5 placettes, lesdites placettes étant éparpillées sur une zone représentative de l'état de croissance de la culture sur la parcelle. Il est également possible de mesurer cette masse de manière indirecte soit en s'aidant de la règle d'équivalence indiquée aux figures 32A à 32E, soit à l'aide de techniques d'imageries qui peuvent être obtenues par différents moyens de capture d'image (depuis un appareil photo, possiblement installé sur un appareil portatif, un appareil volant comme un drone, ou encore au moyen d'images satellites). Grâce à un traitement de l'image, il sera possible d'estimer la masse de matière fraiche aérienne. L'exploitant, s'il souhaite disposer d'une réponse plus précise quant à la quantité d'azote à ajouter à la culture devrait, outre la mesure de la MVSH, effectuer la même 30 mesure sur ladite culture de colza à l'entrée de l'hiver de l'année n où il a planté le colza. Les techniques de mesure décrites ci-dessus sont applicables. Une fois cette mesure effectué, l'exploitant a besoin, pour mettre en oeuvre la méthode de l'invention, de disposer d'informations relatives à sa culture : - sa localisation géographique ainsi que la nature du sol sur lequel ladite culture est 35 implantée, et - l'objectif de rendement pour sa culture. Cet objectif de rendement est un objectif de rendement aux normes de commercialisation (lot de graines renfermant 2% - 7 - d'impuretés et 9% d'eau). Les impuretés correspondent à tout ce qui n'est pas des graines : petites mottes de terres, fragments de tiges et de siliques qui n'ont pas pu être séparés des graines lors du battage et du nettoyage dans la moissonneuse batteuse. Le teneur en eau s'applique sur les graines pures (après séparation des impuretés). A la récolte, les lots de graines ne présentant jamais exactement ces caractéristiques, on mesure donc le taux d'impuretés et la teneur en eau réelle et par calcul, puis on détermine un équivalent de rendement aux normes Ces informations, outre le MVSH, sont indispensables pour mettre en oeuvre la méthode de l'invention. Il est également possible d'affiner le résultat en disposant d'informations concernant l'azote minéral présent dans le sol où est implantée ladite culture de colza. En effet, en disposant de telles informations, il est possible de tenir compte de l'azote du sol susceptible d'être utilisé par la culture, et donc de réduire ou d'augmenter la quantité 15 d'engrais azoté à ajouter, tout en maintenant une fertilisation optimale. L'azote présent dans le sol peut être apporté de différentes manières, notamment du fait des pratiques culturales utilisées sur la parcelle. En effet, si la culture de colza est cultivée sur une parcelle sur laquelle il y a eu lors des campagnes de cultures précédentes des apports d'engrais azoté organiques, il est 20 opportun d'en tenir compte, car une partie de l'azote est disponible dans le sol et pourra, outre l'apport supplémentaire, participer à la fertilisation du colza. De plus, les techniques de culture du colza avec des légumineuses gélives, qui capteront l'azote à la fin de l'été et à l'automne, mais gèleront en hiver, permettront d'améliorer la fourniture d'azote par le sol en une certaine quantité, cette certaine 25 quantité d'azote correspondant à autant d'engrais en moins à ajouter à la culture à la sortie de l'hiver. Dans l'invention on parle donc de plante gélive lorsque ladite plante est sensible au gel. Il est également tenu compte de la nature de la culture sur la parcelle pendant la 30 compagne précédente, c'est-à-dire la culture qui était présente sur la parcelle juste avant la culture du colza, dans le cas où il s'agit de pois protéagineux qui est également susceptible d'augmenter la fourniture d'azote par le sol à la culture de colza. A l'aide des informations préalablement citées, il est alors possible pour l'exploitant de calculer la quantité d'engrais à ajouter pour fertiliser la culture à l'aide des formules 35 susmentionnées. Le calcul peut être manuel ou être mis en oeuvre au moyen d'un programme d'ordinateur installé sur un ordinateur, ou en connexion via l'internet sur un serveur à - 8 - distance, en particulier le serveur de la demanderesse dédié à cet effet. Selon la région envisagée, la méthode de calcul sera différente. Si la culture a lieu dans la région Nord-Est, définie dans l'invention comme recouvrant les départements de la Haute-Marne, de la Meuse, de la Meurthe-et5 Moselle, de la Moselle, des Vosges, du Bas-Rhin et du Haut-Rhin, ou dans la région Auvergne-Rhône Alpes, définie dans l'invention comme recouvrant les départements de l'Ain, de l'Allier, de l'Ardèche, du Cantal, de la Drôme, de la Haute-Loire, de l'Isère, de la Loire, du Puy-de-Dôme, du Rhône, de la Savoie et de la Haute-Savoie, la formule de calcul utilisée sera la suivante : X = (((b.y) - (Pi + Np)) / CAU) - (Mha + Mpro1 + Fleg + 10 Fass). Cette équation correspond à la variante CAU de la méthode de l'invention. Si par contre la culture est réalisée sur le reste du territoire français métropolitain, c'est-à-dire en excluant les territoires et département d'outre-mer et la Corse, la formule de calcul utilisée sera la suivante X = [(b.y) + Rf] - (Pi + Ri + M + Mha + Mpro1 + Fleg + Fass), Cette équation correspond à la variante additive de la méthode de l'invention. 15 Dans le cadre de la variante additive de la méthode selon l'invention, il sera pris en considération deux régions spécifiques : la région Champagne définie dans l'invention comme regroupant les départements de l'Aisne, des Ardennes, de l'Aube et de la Marne, et de la région Grand-Ouest, définie dans l'invention comme regroupant les départements du Finistère, des Côtes d'Armor, du Morbihan, de la Loire Atlantique, 20 d'Ille et Vilaine, la Manche, la Mayenne, le Maine et Loire, la Sarthe et l'Orne. La région Champagne selon l'invention est importante car outre les types de sol retrouvés dans les autres départements (sols superficiels ou sols profonds), il existe dans cette région un autre type de sol : les craies de Champagne. La région Grand-Ouest est également importante du fait de sa spécificité en terme 25 de fourniture d'azote par le sol (M et Ri) Dans l'invention on entend par « ouverture du bilan » la période se situant aux alentours de l'acquisition du stade de développement de la culture de colza défini comme stade 30 BBCH qui correspond au début d'élongation de la tige principale lorsqu'il n'y a pas d'entrenceud (rosette). 30 On entend par ailleurs dans l'invention par « fermeture du bilan » la date ou la période correspondant à la fin d'absorption d'azote par la culture de colza. Dans l'invention, on définit par « quantité optimale d'engrais azoté », la quantité la plus appropriée d'engrais à ajouter à la culture de colza pour l'élaboration de son rendement. 35 Grace aux règles de calcul susmentionné, qui tiennent compte des paramètres susmentionnés, l'exploitant obtient une dose d'azote à apporter, à la fin de l'hiver, afin de stimuler la croissance de ladite culture et d'espérer obtenir l'objectif de rendement -9 indiqué. En outre, une fois cette valeur obtenue, l'exploitant peut déterminer, à l'aide de ses connaissances générales, le ou les types d'engrais à utiliser pour fournir la quantité d'azote obtenue par la méthode selon l'invention. Il a alors la possibilité de peser le ou 5 les engrais qui sera ou seront épandu(s) sur la culture. Dans l'hypothèse où la valeur X obtenue grâce au calcul susmentionné est nulle, cela signifie que le sol fournira suffisamment d'azote pour permettre le rendement espéré, et l'exploitant n'aura donc pas à ajouter d'azote sous forme d'engrais à la culture. 10 Avantageusement, l'invention concerne une méthode de fertilisation telle que définie ci-dessus, comprenant en outre une étape d'épandage, au printemps de la quantité X d'engrais à ladite culture de colza. Plus avantageusement, l'invention concerne une méthode de fertilisation susmentionnée, dans laquelle 15 Rf varie de 10 à 60, Ri varie de 10 à 50, M varie de 10 à 100, Mha varie de 0 à 100, Keq varie de 0 à 0,4 20 Fleg, et Fass varient de 0 à 50, Coeff SH vaut 65, et coeff EH vaut 50. Dans ce mode de réalisation avantageux, les différents paramètres susmentionnés sont définis dans les gammes indiquées. Rf varie de 10 à 60, ce qui signifie que Rf peut prendre les valeurs suivantes : 10, 25 30 35 Mha varie de 0 à 100, ce qui signifie que Mha peut prendre les valeur suivantes : 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59 ou 60. Ri varie de 10 à 50 ce qui signifie que Ri peut prendre les valeurs suivantes : 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49 ou 50. M varie de 10 à 100, ce qui signifie que M peut prendre les valeurs suivantes : 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99 ou 100. - 10 - 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99 ou 100. Keq varie de 0 à 0,4 signifie que Keq peut prendre les valeurs suivantes : 0, 0,01, 0,02, 0,03, 0,04, 0,05, 0,06, 0,07, 0,08, 0,09, 0,1, 0,11, 0,12, 0,13, 0,14, 0,15, 0,16, 0,17, 0,18, 0,19, 0,2, 0,21, 0,22, 0,23, 0,24, 0,25, 0,26, 0,27, 0,28, 0,29, 0,3, 0,31, 0,32, 0,33, 0,34, 0,35, 0,36, 0,37, 0,38, 0,39 ou 0,4. Fleg, Fass varient de 0 à 50, ce qui signifie qu'indépendamment les uns des autres ils peuvent prendre les valeurs suivantes : 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49 ou 50. Dans un mode de réalisation avantageux l'invention concerne une méthode telle que définie ci-dessus, d'une culture de colza localisée dans une région de France 15 métropolitaine qui n'est pas à la région Nord-Est ou la région Auvergne-Rhône-Alpes, dans laquelle : -si le sol est un sol Nss, Rf vaut 15, Ri vaut 30 pour le Grand Ouest [et 20 pour le reste de la France, M vaut 30 le Grand Ouest et 20 pour le reste de la France, -si le sol est un sol Nsp, Rf vaut 30, Ri vaut 40 pour le Grand Ouest et 30 pour le 20 reste de la France, M vaut 50 le Grand Ouest et 40 pour le reste de la France -si le sol est un sol Nsc, Rf vaut 30, Ri vaut 40, et M vaut 20. Dans un mode de réalisation avantageux l'invention concerne une méthode telle que définie ci-dessus, d'une culture de colza localisée dans la région Nord-Est ou dans la région Auvergne-Rhône-Alpes 25 si le sol est un sol Nss, Np vaut 30 pour la région Nord-Est ou vaut 40 la région Auvergne-Rhône-Alpes, si le sol est un sol Nsp, Np vaut 50 pour la région Nord-Est ou vaut 60 la région Auvergne-Rhône-Alpes. Dans un autre mode de réalisation, l'invention concerne une méthode telle que 30 définie ci-dessus, dans laquelle si la culture de colza est implantée sur un terrain ayant porté une culture de pois protéagineux la campagne précédente, Fleg vaut 25, sinon Fleg vaut 0, et si la culture de colza est implantée avec des légumineuses gélives associées, Fass faut 30, sinon Fass vaut 0. Avantageusement, l'invention concerne une méthode telle que définie ci-dessus, 35 dans laquelle Mha est déterminé comme indiqué dans le tableau suivant : Produits Résiduaires Organiques utilisé tous les 2 années tous les tous les ans sur 3 2 ans 3 ans et plus Fumier de bovins 30 20 15 5 Fumier de volailles 15 10 5 5 Fumier de porcs 20 15 5 0 Fumier de cheval 40 25 20 10 Fumier d'ovins 35 20 15 10 Lisier de porcs 20 15 10 5 Lisier de bovins 25 15 10 5 Lisier de bovins dilué 10 5 5 0 Lisier, fientes de volailles (60%MS) 10 5 5 0 Fientes de volailles séches (80%MS1 20 10 10 0 Boues urbaines liquides et pâteuses (É20%MS) 15 15 10 5 Boues urbaines sèchées (90%MS) 100 60 50 25 Compost de déchets verts 30 20 15 10 Compost de fumier de bovins, boues STEP + déchets verts 30 20 15 5 Compost de fumier de volailles 25 15 10 10 Vinasse de betterave concentrée 0 0 0 0 Ecu mes de sucreries surpressées 5 5 5 0 En d'autres termes, selon le type de produit résiduaire utilisé, et la récurrence d'épandage, l'homme du métier, à l'aide du tableau ci-dessus sera capable de déterminer la valeur Mha. Par exemple si du lisier de bovin est ajouté sur la parcelle tous les deux ans, Mha 5 aura une valeur de 10. De la même manière, si du compost de déchets verts est épandu tous les ans, Mha aura une valeur de 30. L'homme du métier comprendra ainsi comment utiliser le tableau ci-dessus. Dans un autre mode de réalisation avantageux, l'invention concerne une méthode 10 telle que définie ci-dessus, dans laquelle Mpro1 = teneur en azote des produits résiduaires (NPRO) x quantité de produits résiduaires épandue (qPRO) x le coefficient d'équivalence engrais (keq), où NPRO, qPRO et Keq sont définis dans le tableau suivant : - 12 - Produit Résiduaire Organique NPRO qPRO unité (t ou m3 / ha) Keq Fumiers Fumier de bovins 5.8 20 t / ha 0.1 Fumier de volailles 25 5 t / ha 0.2 Fumier de porcs 8 20 t / ha 0.1 Fumier de cheval 8 20 t / ha 0.1 Fumier d'ovins 7 20 t / ha 0.1 Lisiers Lisier de porcs à l'engrais 5.8 20 m3 / ha 0.05 Lisier de porcs mixtes 3.5 30 m3 / ha 0.05 Lisier de bovins 4.5 30 m3 / ha 0.15 Lisier de bovins dilué 1.6 30 m3 / ha 0.15 Lisier de canards 7.7 20 m3 / ha 0.05 Fientes Fientes de volailles séches (80%MS) 40 5 t /ha 0.05 Fientes de volailles (60%MS) 24 5 t /ha 0.05 Boues urbaines Boues urbaines liquides 3 30 m3 / ha 0.15 Boues urbaines pâteuses (20%MS) 10 10 t! ha 0.1 Boues urbaines sèchées (90%MS) 40 10 t! ha 0.1 Composts Compost de déchets verts 10 10 t / ha 0.05 Compost de fumier de bovins 6.3 15 t / ha 0.12 Compost de fumier de volailles 23 3 t / ha 0.12 Autres Vinasse de betterave concentrée 20 3 t / ha 0.15 Ecumes de sucreries surpressées 4 10 t / ha o Autre produit organique "champ libre" "champ libre" t / ha 0.1 et où champ libre correspond à des valeurs fournies par l'exploitant. En d'autres termes, si le produit résiduaire est du fumier de porc, épandu à une quantité de 20t/ha, MPRO sera alors égal à 8x20x0,1= 16. L'homme du métier sera capable à l'aide du tableau précédent de déterminer la 5 valeur de MPRO. Dans encore un autre mode de réalisation avantageux, l'invention concerne une méthode telle que définie ci-dessus, dans laquelle (i) si la culture de colza est localisée dans une région de France métropolitaine qui n'est pas la région Nord Est ou la région Auvergne- 10 Rhône-Alpes, la quantité de d'azote à apporter X, au printemps, à ladite culture, exprimée en kg d'azote par hectare, est calculée selon la formule suivante ; - 13 - O si la région est le Grand Ouest - si le sol est un sol Nss X = [(7y) + 15] - (Pi + 60 + Mha + Mpro1 + Fleg + Fass), - si le sol est un sol NSP X = [(7y) + 30] - (Pi + 90 + Mha + Mpro1 + Fleg + Fass), O si la région n'est pas le Grand Ouest - si le sol est un sol Nss X = [(7y) + 15] - (Pi + 40 + Mha + Mpro1 + Fleg + Fass), - si le sol est un sol NsP X = [(7y) + 30] - (Pi + 70 + Mha + Mpro1 + Fleg + Fass), - si le sol est un sol Nsc X = [(7y) + 30] - (Pi + 60 + Mha + Mpro1 + Fleg + Fass), (ii) si la culture de colza est localisée dans la région Nord-Est ou la région Auvergne-Rhône-Alpes, la quantité de d'azote à ajouter X, au printemps, à ladite culture, exprimée en unité d'azote, est calculée selon la formule suivante ; O si la région est la région Nord-Est - si le sol est un sol Nss X = (((7y) - (Pi + 30)) / 0.8) - (Mha + Mpro1+ Fleg + Fass), - si le sol est un sol NSP X = (((7y) - (Pi + 50)) / 0.8) - (Mha + Mprol+ Fleg + Fass), O si la région est la région Auvergne-Rhône-Alpes - si le sol est un sol Nss X = (((7y) - (Pi + 40))! 0.8) - (Mha + Mpro1+ Fleg + Fass), - si le sol est un sol Np X = (((7y) - (Pi + 60)) / 0.8) - (Mha + Mpro1+ Fleg + Fass). Ce mode de réalisation est le plus avantageux. Il permet à l'exploitant de déterminer au plus juste, en tenant compte des contraintes réglementaires, la quantité d'azote à ajouter à la culture de colza pour obtenir le rendement espéré. L'invention concerne en outre un produit programme d'ordinateur sur un support approprié pour la mise en oeuvre de la méthode définie précédemment. Un autre aspect de l'invention porte sur un logiciel ou un produit programme d'ordinateur conçu pour mettre en oeuvre le procédé décrit ci-dessus et/ou comprenant des portions/moyens/instructions de code de programme pour l'exécution dudit procédé lorsque ledit programme est exécuté sur un ordinateur. Avantageusement ledit programme est compris dans un support d'enregistrement de données lisible par ordinateur. Un tel support n'est pas limité à un support d'enregistrement portable tel - 14 - qu'un CD-ROM mais peut également faire part d'un dispositif comprenant une mémoire interne dans un ordinateur (par exemple des RAM et/ou ROM), ou de dispositif à mémoires externes tels des disques durs ou des clefs USB, ou un serveur à proximité ou à distance. L'invention sera mieux comprise à la lumière des 32 figures et des exemples suivants : Brève description des figures : Les figures 1A et 1B représentent des graphiques montrant le rendement aux 10 normes. La figure 1A représente le rendement aux normes en fonction du poids de matière fraîche aérienne à la sortie de l'hiver dans le référentiel d'essais doses d'azote 2008 à 2012. La figure 1B représente le rendement aux normes en fonction de la quantité d'azote 15 absorbé à la sortie de l'hiver dans le référentiel « observatoire région Centre ». Les figures 2A et 2B sont des graphiques montrant un exemple d'essai satisfaisant les critères de sélection pour le calcul du besoin unitaire. La figure 2A représente le rendement aux normes en fonction de la dose d'azote. Il 20 existe une différence significative entre les doses N dans la suite 2, une absence de différence significatives dans la suite 1 , l'absence de rendement classé « non a » au sein de la suite 1 et l'absence de rendement classé « a » dans la suite 2. La figure 2B représente l'azote absorbé (Nabs pe G4 = quantité d'azote absorbé dans les plantes entières au stade G4) en fonction de la dose d'azote. On observe une 25 linéarité de la réponse de Nabs pe G4 pour les doses 0 à 280 kgN/ha. Les figures 3A à 3F représentent les effectifs pour quelques types de segmentation de la base de données retenue pour l'estimation du besoin unitaire (22 essais) (MV SH : biomasse fraîche aérienne à la sortie de l'hiver ; Rdt max : rendement maximal 30 enregistré sur l'essai) La figure 3A représente la répartition selon l'année. La figure 3B représente la répartition selon la région, La figure 3C représente la répartition selon le type variétal. La figure 3D représente la répartition selon MV SH. 35 La figure 3E représente la répartition selon le rendement maximum. La figure 3F représente la répartition selon le type de sol. - 15 - Les figures 4A et 4B représentent des graphiques montrant des exemples de calcul du besoin unitaire sur un essai et pour deux hypothèses de niveau d'atteinte du rdtmaxQP. (Figure 4A: scénario 1 ; figure 4B scénario 2). Les figures 5A à 5C représentent le besoin unitaire en fonction du rendement maximal pour différents niveaux de rendement maximal (rdtmaxQP, 99% (figure 5A) de rdtmaxQP (figure 5B) et 95% de rdtmaxQP (figure 5C)). Référentiel de base = 23 essais ; 1 essai éliminé car valeur de « b » anormalement élevée (point barré) ; restent 22 essais sur lesquels ont été ajustées les droites de régression (Excel 2000). La figure 6 représente le besoin unitaire en fonction du rendement maximal (rdtmaxQP) et valeurs de statistique descriptives (déciles). Référentiel de 22 essais sur lesquels ont été calculés les déciles (Excel 2010). A noter, le changement d'échelle sur l'axe des ordonnées par rapport aux figures 5A à 5C). Les figures 7A et 7B représentent des graphiques illustrant deux catégories d'essai : rendement maximal modélisé (plateau) atteint pour une quantité d'azote absorbé à G4 plus faible (figure 7A) ou plus élevée (figure 7B) que la quantité maximale d'azote absorbé observée dans l'essai (courbe : relation modélisée ; points : données observées). Les figures 8A à 8F représentent les effectifs (nombre d'essais) pour quelques types de segmentation de la base de données retenue pour l'estimation du plafond de besoin (44 essais) (MV SH : biomasse fraîche aérienne à la sortie de l'hiver ; Rdt max : 25 rendement maximal enregistré sur l'essai). La figure 8A représente la répartition selon l'année. La figure 8B représente la répartition selon la région, La figure 8C représente la répartition selon le type variétal. La figure 8D représente la répartition selon MV SH. 30 La figure 8E représente la répartition selon le rendement maximum. La figure 8F représente la répartition selon le type de sol. Les figures 9A et 9B représentent des graphiques montrant les courbes de réponse du rendement relatif à la quantité d'azote absorbé dans les plantes entières au stade 35 G4 ; (figure 9A), 31 essais pour lesquels le rendement maximal est obtenu avec N absorbé à G4 inférieur ou égal à N absorbé à G4 maximal de l'essai ; (Figure 9B) : 10 essais pour lesquels le rendement maximal est obtenu avec N absorbé à G4 - 16 - légèrement supérieur (écart La figure 10 représente la quantité d'azote absorbé dans les plantes entières (kgN/ha) en fonction du poids de matière fraîche aérienne (g/m2) à l'entrée de l'hiver et pentes des régressions linéaires (ordonnée à origine = 0) par région. La figure 11 représente la quantité d'azote absorbé dans les plantes entières 10 (kgN/ha) en fonction du poids de matière fraîche aérienne (g/m2) à l'entrée de l'hiver et pentes des régressions linéaires (ordonnée à origine = 0) par année. La figure 12 représente la quantité d'azote absorbé dans les plantes entières (kgN/ha) en fonction du poids de matière fraîche aérienne (g/m2) à la sortie de l'hiver et 15 pentes des régressions linéaires (ordonnée à origine = 0) par région. La figure 13 représente la quantité d'azote absorbé dans les plantes entières (kgN/ha) en fonction du poids de matière fraîche aérienne (g/m2) à la sortie de l'hiver et pentes des régressions linéaires (ordonnée à origine = 0) par année. 20 La figure 14 représente la quantité d'azote absorbé dans les plantes entières (kgN/ha) en fonction du poids de matière fraîche aérienne (g/m2) à l'entrée (EH) et à la sortie (SH) de l'hiver et pentes des régressions linéaires (ordonnée à origine = 0) par région. 25 Les figures 15A à 151 représentent la quantité d'azote absorbé dans les plantes entières (kgN/ha) en fonction du poids de matière fraîche aérienne (g/m2) à l'entrée de l'hiver par région, toutes années confondues et coefficients de détermination (r2) des régressions linéaires et quadratiques (ordonnée à origine = 0). 30 La figure 15A représente les données pour la France. La figure 15B représente les données pour la Bourgogne. La figure 15C représente les données pour la région Centre. La figure 15D représente les données pour la région Champagne-Ardenne Nord- Picardie. 35 La figure 15E représente les données pour la région Nord-Est La figure 15F représente les données pour la région Poitou-Charentes Vendée Limousin. - 17 - La figure 15G représente les données pour la région Rhône-Alpes. La figure 15H représente les données pour la région Sud. La figure 151 est un tableau récapitulatif des données des figures 15A à 15H. Les figures 16A à 16X représentent des graphiques montrant la quantité d'azote absorbé dans les plantes entières (kgN/ha) en fonction du poids de matière fraîche aérienne (g/m2) à l'entrée de l'hiver par « région x année » et régressions linéaires et quadratiques toutes années confondues (lignes noires) et par année (ordonnée à origine = 0) (points noirs : années concernée ; points gris : autres années). La figure 17 représente la quantité d'azote absorbé dans les plantes entières (kgN/ha) en fonction du poids de matière fraîche aérienne (g/m2) à la sortie de l'hiver toutes régions, toutes années confondues et coefficients de détermination (r2) des régressions linéaires et quadratiques (ordonnée à origine = 0). Les figures 18A à 181 représentent la quantité d'azote absorbé dans les plantes entières (kgN/ha) en fonction du poids de matière fraîche aérienne (kg/m2) à la sortie de l'hiver par région, toutes années confondues et coefficients de détermination (r2) des régressions linéaires et quadratiques (ordonnée à origine = 0). La figure 18A représente les données pour la région Auvergne. La figure 18B représente les données pour la Bourgogne. La figure 18C représente les données pour la région Centre. La figure 18D représente les données pour la région Champagne-Ardenne Nord- Picardie. La figure 18E représente les données pour la région Nord-Est La figure 18F représente les données pour la région Poitou-Charentes Vendée Limousin. La figure 18G représente les données pour la région Rhône-Alpes. La figure 18H représente les données pour la région Sud. La figure 181 est un tableau récapitulatif des données des figures 18A à 18H. Les figures 19A à 19AQ représentent des graphiques représentant la quantité d'azote absorbé dans les plantes entières (kgN/ha) en fonction du poids de matière fraîche aérienne (kg/m2) à la sortie de l'hiver par « région x année » et régressions linéaires et quadratiques toutes années confondues (lignes noires) et par année (ordonnée à origine = 0) (points noirs : années concernée ; points gris : autres années). - 18 - La figure 20 est un graphique montrant la teneur en azote des parties aériennes (% de matière sèche) en fonction du poids de matière sèche aérienne (g/m2) à l'entrée de l'hiver toutes régions confondues et courbe critique de teneur en azote (d'après Colnenne, 1998). La figure 21 est un graphique montrant la teneur en azote des parties aériennes (% de matière sèche) en fonction du poids de matière sèche aérienne (g/m2) à la sortie de l'hiver toutes régions confondues et courbe critique de teneur en azote (d'après Colnenne, 1998). Les figures 22A et 22B représentent des graphiques montrant la relation entre la quantité d'azote absorbé dans les parties aériennes (Nabs aérien) et le poids de matière sèche aérienne (MS) (figure 22A) et entre le poids de matière sèche aérienne (MS) et le poids de matière fraîche aérienne (MV) (figure 22B) à l'entrée (EH) et à la sortie de l'hiver (SH) (base de données : BDD EH_SH). La figure 23 représente un graphique montrant la quantité d'azote dans les plantes entières (kgN/ha) en fonction du poids de matière fraîche aérienne (kg/m2) à l'entrée de l'hiver ; droites et pentes des droites correspondant aux déciles 1, 3, 5, 7 et 9. La figure 24 représente un graphique montrant la quantité d'azote dans les plantes entières (kgN/ha) en fonction du poids de matière fraîche aérienne (kg/m2) à la sortie de l'hiver ; droites et pentes des droites correspondant aux déciles 1, 3, 5, 7 et 9. Les figures 25A à 25C représentent des histogrammes montrant les principales caractéristiques statistiques des bases de données nettoyées (SAS figure 25A, CA28 figure 25B et CAMA figure 25C) concernant le reliquat d'azote minéral du sol à l'ouverture du bilan (kgN/ha) La figure 26 représente le nombre d'analyses par type de sol et par département dans la base de données « reliquat azoté à la sortie de l'hiver sous colza » du SAS. La figure 27 est un graphique montrant le rendement du colza en fonction du précédent cultural et de la dose d'engrais N : exemple d'une expérimentation conduite 35 par la demanderesse dans le département de la Somme en 2009/2010. Les figures 28A à 28F représentent les effectifs (nombre d'essais) pour quelques - 19 - types de segmentation de la base de données retenue pour l'estimation du CAU (25 essais) (MV SH : biomasse fraîche aérienne à la sortie de l'hiver ; Rdt max : rendement maximal enregistré sur l'essai). La figure 28A représente la répartition selon l'année. La figure 28B représente la répartition selon la région, La figure 28C représente la répartition selon le type variétal. La figure 28D représente la répartition selon MV SH. La figure 28E représente la répartition selon le rendement maximum. La figure 28F représente la répartition selon le type de sol. La figure 29 représente un histogramme montrant la distribution des coefficients apparents d'utilisation de l'azote de l'engrais minéral et données statistiques simples. Les figures 30A et 30B sont des graphiques représentant des courbes de réponse modélisées (ajustement quadratique + plateau) du rendement relatif (rendement modélisé à la dose N / rendement modélisé à la dose d'azote apporté la plus élevée de chaque essai) en fonction de la dose d'azote apporté pour les essais avec (Figure 30A; 43 essais) ou sans atteinte d'un plateau de rendement (Figure 30B; 25 essais). La figure 31 représente les types de méthodes de l'invention utilisées selon les régions. Les figures 32A à 32E représentent une corrélation entre le colza observé à la sortie de l'hier, et le poids frais aérien estimé en kg/m2. Exemples : Exemple 1 : Exemple de comparaisons des paramètres de la réglette 1998 et de la méthode selon l'invention. L'objectif de la demanderesse est de construire un outil d'aide à la fertilisation azotée 30 du colza basé sur les principes du bilan prévisionnel (approche statique et non dynamique), mais avec un certain nombre de simplifications pour le rendre plus opérationnel. De nos jours, compte-tenu de l'obligation d'établir des plans prévisionnels de fertilisation azotée multicultures dans la majeure partie des zones de grandes cultures 35 (zones vulnérables), les outils d'aide à la décision multi-espèces semblent les mieux adaptés à une utilisation par les agriculteurs et leurs techniciens (appropriation d'un seul outil pour toutes les cultures). La demanderesse a toutefois fait le choix de - 20 - développer un outil spécifique au colza pour maintenir un affichage clair sur les cultures dont elle s'occupe et pour donner à l'agriculteur la possibilité d'utiliser un outil simple. Par rapport à l'outil préexistant de calcul de la quantité d'azote à ajouter, (appelée « réglette 1998 »), l'objectif est : - d'actualiser un certain nombre de paramètres avec des bases données récentes et étoffées : o le besoin unitaire en azote (b), o le coefficient de conversion de la biomasse fraîche en quantité d'azote absorbé (coeff), o les reliquats d'azote minéral du sol à l'ouverture et à la fermeture du bilan (Ri et Rf), o le coefficient apparent d'utilisation (CAU) de l'azote des feuilles vertes gelées pendant l'hiver ; - de réduire les différences de paramètres entre régions (homogénéiser ou justifier les différences) : b, coeff, Ri, Rf, minéralisation de l'humus et des résidus du précédent (Mn), fourniture d'azote par le sol (Np) ; - de mieux prendre en compte les situations particulières : apports de produits résiduaires organiques (PRO) ponctuels sur le colza et réguliers dans la rotation (diversité des produits, quantités et fréquence), association de cultures (cultures compagnes de type légumineuses gélives), précédent légumineuse (restreint pour le moment au pois protéagineux) ; - de compléter le référentiel de sols (cas des craies de Champagne), - de prendre en compte les dernières recommandations du COMIFER pour limiter le risque de volatilisation et ajuster les doses d'engrais si besoin. 1-Types de sols Les regroupements de types de sols en deux catégories « sols superficiels » et « sols profonds » sont maintenus à l'identique de la réglette 1998 : - sols superficiels : caractérisés par une faible minéralisation au printemps et une réserve en eau faible, susceptible de limiter le rendement ; - sols profonds : caractérisés par une minéralisation « moyenne » à « élevée » au printemps et une bonne réserve en eau ne limitant pas (ou peu) le rendement. Un nouveau type de sol « sols de craie de Champagne » est créé pour prendre en 35 compte les caractéristiques spécifiques de ces sols : minéralisation printanière moyenne et bonne réserve en eau ne limitant pas (ou peu) le rendement, même en situation de Pi faible. -21 - 2-Nouveaux paramétrages 2.1 Objectif de rendement (y) Dans les zones vulnérables, l'objectif de rendement doit être fixé en respectant les arrêtés préfectoraux établissant les référentiels régionaux de mise en oeuvre de l'équilibre de la fertilisation azotée. Ailleurs, la règle proposée par la demanderesse est la suivante : moyenne des rendements réalisés en colza sur l'exploitation au cours des 5 dernières campagnes successives, si possible pour des conditions comparables de sol et de précédents de culture, en excluant la valeur minimale et la valeur maximale. La réglette azote 1998 introduit une limitation de l'objectif de rendement lorsque la biomasse aérienne fraîche à l'ouverture du bilan est trop faible. Les objectifs de rendement supérieurs à ces limites sont des situations dites « non réalistes » comme indiqué dans le tableau suivant : MV Objectif de rendement Sol superficiel Sol profond (kg/m2) Bourgogne 0,2 40 NR NR 0,4 40 NR NR 0,6 40 NR Champagne 0,2 40 NR NR 0,4 40 NR Centre 0,2 40 NR NR 0,4 40 NR NR 0,6 40 NR 0,8 40 NR (Sans PRO) Nord-Est 0,2 35 NR 0,2 40 NR NR 0,4 40 NR NR 0,6 40 NR 0,8 40 NR RA 0,2 35 NR (Sans PRO) 0,2 40 NR NR 0,4 40 NR NR 0,6 40 NR (Sans PRO) 15 - 22 - MV Objectif de Sol superficiel Sol profond (kg/m2) rendement PC 0,2 35 NR 0,2 40 NR NR 0,4 40 NR NR 0,6 40 NR NR 0,8 40 NR Auvergne 0,2 35 NR 0,2 40 NR NR 0,4 40 NR NR 0,6 40 NR (Sans PRO) Sud 0,2 35 NR NR 0,2 40 NR NR 0,4 40 NR NR 0,6 40 NR NR 0,8 40 NR Bretagne 0,2 40 NR NR 0,4 40 NR NR 0,6 40 NR NR 0,8 40 NR Le tableau 1 fait un inventaire des situations « non réalistes (plafond d'objectif de rendement) avec la réglette 1998(NR = non réaliste ; NR sans PRO: NR seulement dans les situations sans apport de produits résiduaires organiques ; case vide : calcul 5 possible) Le calcul n'est pas effectué pour ces situations. Ces plafonds sont très divers en fonction des régions et ils sont jugés plus ou moins restrictifs, voire trop restrictifs au regard des références aujourd'hui disponibles (essais doses d'azote et observatoire en 10 région Centre en (Figures 1A et 1B). Il est en effet possible d'atteindre des rendements élevés avec de faibles biomasses, mais on constate toutefois une limitation du potentiel en dessous de 0.4 kg/m2. Pour la méthode selon l'invention, il est donc proposé de ne jamais empêcher le calcul, mais d'alerter l'utilisateur sur la cohérence de sa saisie sur la base des plafonds 15 de rendement proposés en fonction du poids de matière fraîche à la sortie hiver et du type de sol (voir tableau 2). Il lui reviendrait ensuite de modifier sa valeur d'objectif de - 23 - rendement ou de la confirmer avant de lancer le calcul. Gamme de biomasse Sols Sols profonds Sols de craie fra îche aérienne superficiels de Champagne PV 200 g/m2 35q/ha 40 q/ha 40 q/ha 200 g/m2 Dans la réglette azote 1998, deux valeurs étaient retenues : 7.0 pour les régions Centre, Poitou-Charentes, Vendée, Limousin et Sud et 6.5 pour toutes les autres régions. Ce besoin unitaire est calculé à l'optimum de fertilisation azotée vis-à-vis du rendement. Pour l'atteindre, il convient donc de disposer d'essais « doses d'azote » permettant d'identifier cet optimum (courbe de réponse du rendement à la dose avec une partie croissante pour les doses les plus faibles et un plateau pour les doses les plus élevées). 2.2.1 Matériel et méthodes Parmi les 69 essais « doses d'azote » conduits par la demanderesse de 2008 à 20 2012 (thèmes MFN, MFO et FAZ), seuls ont été sélectionnés ceux présentant les caractéristiques suivantes (Figures 2A et 2B présentent un exemple d'essai retenu) : Dose optimale de fertilisation identifiable vis-à-vis du rendement, c'est à dire : o rendements significativement différents entre dose d'azote d'après le test Student - Newman - Keuls (snk) au seuil de 10 %; 25 o existence d'une suite (suite 1) d'au moins 2 doses consécutives dont les rendements sont classés « a » par le test snk au seuil de 10 °A parmi les doses les plus élevées ; o existence d'une suite (suite 2) d'au moins 2 doses consécutives (inférieures à la plus faible dose de la suite de doses élevées classée « a ») dont les 30 rendements ne sont pas classés « a » par le test snk au seuil de 10 %; o pas de dose dont le rendement est classé « a » au sein de la suite « 2 » o ajustement possible des points à une fonction polynomiale d'ordre 2 (quadratique) + plateau (logiciel SAS) avec dose de début de plateau de rendement -24 - modélisée inférieure à la dose maximale de l'essai mesures de quantité d'azote absorbé dans les parties aériennes à la fermeture du bilan (NabsG4 aérien) disponibles pour chaque modalité avec les contraintes suivantes o dans le cas du colza, c'est le stade G4 « azote » (de l'ordre de 2 à 3 semaines après fin floraison) qui est retenu pour cette estimation ; o il doit exister une relation linéaire correcte entre la quantité d'azote absorbé dans les parties aériennes au stade G4 et la dose d'azote lorsque l'on ne conserve que les doses d'azote apportées inférieures ou égales et celle immédiatement supérieure à la dose correspondant au début du plateau du rendement ; o la quantité d'azote absorbé dans les plantes entières au stade G4 est calculée de la façon suivante : NabsG4 = Nabs G4 aérien x 1.25. Après sélection sur les critères précédemment cités, 22 essais sont retenus pour le calcul du besoin unitaire. Quelques caractéristiques de la base de données sont présentées en Figures 3A à 3F. Il convient de noter la forte représentation des années 2008 et 2009, des régions Centre et Poitou-Charentes et des sols profonds. Le calcul du besoin unitaire est réalisé de la façon suivante pour chaque essai (Figures 4A et 4B) : - 1 : calcul du plateau de rendement - 2 : calcul de la dose d'azote (X) la plus faible permettant d'atteindre de rendement maximum (rdtmaxQP = rendement plateau) d'après l'ajustement quadratique + plateau 25 ; - 3: calcul de la quantité d'azote absorbé plante entière au stade G4 (NabsG4_X) pour cette dose X d'après l'ajustement linéaire de la relation NabsG4 = f(doseN). - calcul du besoin unitaire : b = NabsG4_X / X. 30 En phase exploratoire, ce calcul a également été réalisé pour différents niveaux de rendements inférieurs au rendement maximal : 99 % de rdtmaxQP et 95% de rdtmaxQP. 3.2.2 Résultats Les valeurs de besoin unitaires calculées pour les différents essais ont ensuite été 35 représentées en fonction du rendement maximal de référence (rdtmaxQP ou 99% de rdtmaxQP ou 95% de rdtmaxQP) afin de tester l'hypothèse d'une relation entre le besoin et le potentiel de rendement (Figures 5A à 5C). -25- 23 essais ont été retenus sur les critères précédemment mentionnés. Parmi eux, 1 essai a été écarté car il présente une valeur de besoin unitaire anormalement élevée (plus de 10). Restent donc 22 essais sur lesquels ont porté l'ajustement linéaire. Le besoin unitaire et le rendement maximal sont négativement corrélés, mais les coefficients de détermination des droites de régression sont faibles (maximum 0,2) et la pente négative est fortement influencée par 3 points (2 points à gauche dans la partie haute du nuage et le point plus à droite). En conséquence, Les inventeurs ont considéré que le besoin était indépendant du rendement. Il n'a pas non plus été identifié d'effet régional sur la valeur du besoin, ce qui amène à proposer une seule valeur nationale. Comme attendu, les valeurs de besoin unitaire diminuent lorsque l'on prend comme référence des niveaux de rendement qui s'éloigne du rendement maximal (moyennes respectives de 6,45 pour rdtmaxQP, 5,97 pour 99% rdtmaxQP et 5,45 pour 95% de 15 rdtmaxQP). Parmi ces trois scénarios de rendement maximal, celui avec rdtmaxQP a été retenu car l'optimum économique est proche de cette valeur pour une large gamme de prix de la graine et de l'engrais N (résultats non présentés). Les valeurs de besoin unitaire présentent une dispersion (4,82 à 9,13). Dans une 20 situation donnée, la valeur médiane de 6.45 fait donc prendre un risque de surestimation ou de sous-estimation de ce besoin. C'est pourquoi différents niveaux de risque ont été testés en utilisant les déciles 1, 3, 5, 7 et 9 des valeurs de besoin unitaire (calcul des déciles avec Excel 2010 ; Figure 6). Les 5 valeurs correspondantes ont été introduites dans autant de scénarios de paramétrage dans l'étape d'évaluation de la 25 règle de décision globale. Il convient de souligner que les valeurs des déciles 5 et 7 sont proches des valeurs de besoin diffusées jusqu'à maintenant (6,5 et 7,0). 2.3-Plafonnement du besoin de la culture (b.y) 30 Le produit du besoin unitaire et de l'objectif de rendement peut conduire à des valeurs élevées, surtout si la valeur de besoin unitaire retenue est élevée et/ou si l'objectif de rendement est surestimé. Or la réponse du rendement à la quantité d'azote absorbé par la culture à la fermeture du bilan présente la plupart du temps un plateau : au-delà d'une certaine quantité d'azote absorbé, le rendement ne progresse plus. Il 35 s'agit donc ici de définir cette valeur plafond de quantité d'azote absorbé. Dans la réglette azote 1998, le plafond ne concernait que l'écriture additive et était appliqué à la somme du besoin de la culture (b.y) et du reliquat d'azote minéral dans le - 26 - sol à la fermeture du bilan (Rf). Sa valeur était fixée à 330 kgN/ha. Dans le cadre de la mise au point de la méthode de l'invention, les inventeurs ont ré-estimé ce plafond à l'aide de références expérimentales récentes. Les inventeurs ont fait le choix de travailler sur b.y et non plus sur b.y + Rf. 2.3.1 Matériel et méthodes La base de données utilisée est constituée des essais « doses d'azote » conduits par la demanderesse de 2008 à 2012. La relation étudiée est celle exprimant le rendement aux normes en fonction de la quantité d'azote absorbé au stade G4 « azote » dans les plantes entières. N'ont été retenus que les essais pour lesquels cette relation s'ajuste correctement à une fonction quadratique + plateau ou linéaire. Le critère de sélection retenu est la valeur de somme des carrés des erreurs de prédiction (sse) qui doit être inférieure à 25 (cette valeur correspond à la valeur maximale pour laquelle la relation « plait à l'ceil »). 44 essais ont ainsi été retenus sur les 69 constituant la base de données. Sur les 25 essais non retenus, 16 n'ont pas fait l'objet de mesures d'azote absorbé à G4, 3 présentent une courbe de réponse plate et 5 une valeur de « sse » supérieure à 25. Sur les 44 essais retenus, 31 atteignent le plateau de rendement calculé par l'ajustement mathématique pour une quantité d'azote absorbé inférieure à la valeur maximale enregistrée sur l'essai. Pour les 13 autres, le plateau calculé de rendement est atteint pour des valeurs d'azote absorbé supérieures à la valeur maximale observée. Une illustration de ces deux catégories d'essai est présentée en Figures 7A et 7B. La quantité minimale d'azote absorbé permettant d'atteindre le rendement maximal est calculée essai par essai en utilisant les fonctions modélisées (quadratique + plateau ou linéaire). Le rendement maximal correspond au rendement plateau calculé avec le modèle quadratique + plateau pour les essais qui atteignent ce plateau pour une valeur inférieure à la quantité d'azote absorbé maximale de l'essai (31 essais) ou légèrement supérieure (écart de moins de 25 kgN/ha) à la quantité d'azote absorbé maximale de l'essai (10 essais). Quatre essais ne satisfont pas ces critères et ont été étudiées visuellement. Quelques caractéristiques de la base de données sont présentées en Figures 8A à 8F. Il convient de noter la forte représentation des années 2008 et 2009, de la région Centre et des sols profonds. Afin de pouvoir comparer plus facilement les courbes, les rendements de chaque essai ont été exprimés en rendement relatif (= rendement / rendement maximal). Les courbes de réponse de chacun des essais correspondant à ces types d'allure sont présentées en Figures 9A et 9B. La quantité plafond d'azote absorbé au stade G4 a ensuite été déterminée sur la - 27 - base de statistiques descriptives simples de type fréquence d'atteinte du rendement maximal pour différentes hypothèses de quantité d'azote absorbé et déciles d'absorption d'azote (calculés avec Excel 2010). 2.3.2 Résultats Sur la base des 41 essais répondant aux critères de sélection et pour les différentes hypothèses croissantes de quantité plafond d'azote absorbé testées (200 à 360 kgN/ha), la fréquence d'atteinte du rendement maximal augmente, d'abord rapidement jusqu'à 330 kgN/ha (tableau 3 suivant). Avec 330 kgN/ha, le rendement maximal est atteint dans 85 % des essais (35 essais sur 41 essais répondant aux critères de sélection). Au-delà de 330 kgN/ha, cette fréquence n'augmente plus que très lentement. Visuellement, il apparaît clairement que les courbes de réponse des 3 essais qui ne répondant pas aux critères de sélection atteignent le rendement maximal pour une valeur supérieure à 330 kgN/ha absorbé. L'intégration de ces trois essais ne modifie donc pas la valeur plafond d'absorption au-delà de laquelle la fréquence d'atteinte du rendement maximal commence à plafonner. Cette valeur de 330 kgN/ha se situe légèrement au-dessus du 8ème décile des valeurs d'absorption de notre référentiel. minimum 76 décile 1 197 décile 2 215 quartile 1 218 décile 3 221 décile 4 249 décile 5 263 décile 6 278 décile 7 291 quartile 3 314 décile 8 324 décile 9 354 504 maximum NabsG4 pour rendement max (kgN/ha) NabsG4 Fréquence d'atteinte du rendement max (kgN/ha) 200 0.15 250 0.41 280 0.63 290 0.68 300 0.73 310 0.73 320 0.76 330 0.85 340 0.85 350 0.88 360 0.90 Les tableaux 3 représentent la fréquence d'atteinte du rendement maximal pour différentes hypothèses de quantité d'azote absorbé au stade G4 dans les plantes entières (tableau de gauche) et statistiques descriptives simples de la distribution des quantités d'azote absorbé au stade G4 (excel 2010; tableau de droite) (31 + 10 essais répondant aux critères de sélection) Les inventeurs proposent donc de retenir 330 kgN/ha comme valeur plafond de - 28 - besoin (b.y). Le formalisme utilisé sera le suivant : si b.y 30% L 45% CaCO3 40% Sol superficiel (0 à 30 cm) 5 10 15* 15 Peu profond (0 à 60 cm) 10* 15 20 20*** Profond (0 à 90cm) 15* 20 30** 30*** Très profond (0 à 120 cm) 20* 30** 40 40*** Rf : quantité d'azote minéral dans le sol a la fermeture du bilan (kgN/ha) Source : AZOBIL © INRA, 2012 Le tableau 4 regroupe les valeurs de quantité d'azote minéral dans le sol à la fermeture du bilan (sols superficiels avec * ; sols profonds avec ** et sols de craie de Champagne avec ***). Ce tableau concerne des plants à enracinement profond explorant bien le sol (cas des principales grandes cultures : céréales, colza, maïs). Ce tableau peut être adapté en fonction des capacités d'interception de l'azote par les plantes considérées Faute de données nouvelles disponibles pour la demanderesse, les inventeurs se sont appuyés sur les références du COMIFER en les simplifiant à dire d'expert de la façon suivante : - sols superficiels = sols légers peu profonds à profond et sols argileux superficiels - sols profonds = sols argileux profonds et sols limoneux très profonds ; - sols de craie de Champagne = sols de craie peu profonds à très profonds. Il n'est pas introduit de déclinaison régionale. Pour chacun des trois types de sol, les valeurs médianes des cases - 29 - correspondantes sont retenues (tableau 5). Rf (kqN/ha) Réglette 2014 Réglette 1998 Type de sol France Bourgogne, Reste France Bretagne, Centre Sols superficiels 15 20 30 Sols profonds 30 30 40 Sols de craie de Champagne 30 . Le tableau 5 est une table simplifiée des valeurs forfaitaires de quantité d'azote minéral 5 dans le sol à la fermeture du bilan (Rf) proposées pour la méthode de l'invention (la réglette 2014) et rappel des valeurs retenues pour la réglette 1998. 2.5 Quantité d'azote dans la culture en place à l'ouverture du bilan (Pi) 2.5.1 Un indicateur simple : le poids de matière verte aérienne 10 Le principe de l'estimation de la quantité d'azote présent dans la culture à partir du poids de matière verte aérienne a été conservé, faute de disponibilité d'une méthode simple permettant d'accéder plus directement à cette grandeur. 2.5.2 Estimation de la quantité d'azote absorbé par la culture à l'entrée et à la sortie de l'hiver 15 Dans la réglette 1998, le coefficient de conversion de la biomasse fraîche aérienne (kg/m2) en quantité d'azote dans les plantes entières à l'ouverture du bilan (en kgN/ha) est le même à l'entrée et à la sortie de l'hiver. Ce coefficient peut prendre 3 valeurs selon la région : 65, 70 ou 75. Dans le cadre de la mise au point d'une nouvelle règle de décision, la pertinence de 20 ces choix (invariance du coefficient pendant l'hiver et segmentation régionale) a été remise en question. 2.5.2.1 Matériel et méthodes Deux bases de données ont été constituées : une pour l'entrée de l'hiver (BDD_EH) et l'autre pour la sortie de l'hiver (BDD_SH). Chaque base de données regroupe des 25 valeurs de poids de matière verte aérienne (MVaer), de poids de matière sèche aérienne et de teneur en azote total mesurées au champ sur les mêmes échantillons de quelques m2 issus de différentes régions pendant plusieurs années. Les bases de données ont été constituées à partir des dispositifs RCA (Farmstar) 2006 à 2009 et des essais MFN (doses d'azote) 2008 à 2012 ainsi que, seulement pour la sortie hiver, des 30 données de paramétrage de la réglette 1998 (1994 à 1998).- Ces bases de données sont déséquilibrées : il n'y pas le même nombre - 30 - d'échantillons chaque année dans chaque région. Le nombre de références est nettement plus élevé à la sortie qu'à l'entrée de l'hiver en raison du gros effort qui avait été réalisé lors de la construction de la réglette 1998 (travail réalisé seulement à la sortie de l'hiver à l'époque). Toutes les situations ayant fait l'objet de mesures à l'entrée de l'hiver n'en ont pas nécessairement fait l'objet à la sortie de l'hiver. La quantité d'azote absorbé dans les parties aériennes a été calculée en multipliant le poids de matière sèche aérienne par la teneur en azote total. La quantité d'azote absorbé dans les plantes entières (Nabspe) a été obtenue en multipliant la quantité d'azote absorbé dans les parties aériennes par un coefficient égal 10 à 1.35. Ce coefficient établit au cours des années 1990 n'a pas été remis en cause. Chaque base de données a été « nettoyée » en suivant la procédure suivante : - calcul de la pente de la droite de régression passant par (0,0) de la relation exprimant la quantité d'azote absorbé dans les plantes entières en fonction de la biomasse aérienne verte toutes années et régions confondues ; 15 - calcul des valeurs modélisées de quantité d'azote absorbé dans les plantes entières grâce à la pente calculée précédemment et aux valeurs de biomasse aérienne verte ; - calcul des écarts relatifs entre valeurs réelles et valeurs modélisées de quantité d'azote absorbé : (modélisé - réel)/ modélisé ; 20 - calcul de l'écart-type de ces écarts relatifs ; - élimination des points pour lesquels l'écart relatif entre valeurs observées et modélisées est supérieur à 2 fois l'écart-type des écarts relatifs. Le nombre de situations écartées lors de ce nettoyage est de : - 18 sur 319 situations dans la base de donnes initiales à l'entrée de l'hiver (restent 25 301 situations) ; - 16 sur 1253 à la sortie de l'hiver (restent 1237). La répartition des effectifs des bases de données nettoyées en fonction des années et des régions est présentée dans le tableau 6 pour l'entrée de l'hiver et dans le tableau 7 pour la sortie de l'hiver. Pour chaque analyse statistique, les situations présentant un 30 effectif inférieur à 10 ou trop particulières (par exemple, à la sortie de l'hiver, l'année 1998 dans la comparaison des années entre-elles toutes région confondues car une seule région représentée) n'ont pas été retenues (valeurs soulignées dans les tableaux 3.6 et 3.7). - 31 - Effectif 2006 2007 2008 2009 2010 2011 2012 Total Bourgogne 9 14 1 24 Centre 12 11 16 21 1 3 64 Champagne NP 10 10 18 10 2 50 Nord-Est 4 7 1 1 13 Poitou-Charentes 20 10 12 12 54 Rhône-Alpes 14 13 11 9 47 Sud 14 14 11 8 1 1 49 Total 70 58 81 81 2 3 6 301 bleau 6 représente la réoartition des effectifs de la base de données « entrée I » (BDD_EH) en fonction des années et des régions. Le ta iver - 32 - Effecti f a) LO CO h- CO a) CD h- CO CS) 0 n-- (N 73 O) a) a) a) a) a) a) a a 0 a a a a , a .-- a .-- a .,._. .- .- , a) .-- CN a CN CV CV CN CN o .-- C'I H Auvergne 19 28 31 22 19 119 Bourgogne 23 34 44 20 11 15 1 1 149 Centre 56 12 12 21 21 2 1 7 132 Champa gne NP 71 93 14 10 20 7 215 Nord-Est 36 37 29 5 7 1 115 Poitou- 70 56 23 10 12 13 1 185 Charente ci) 24 65 37 49 15 13 12 9 224 c cf) «D 0 _c ca- rz Maximum Mvaer 1994 1995 1996 1997 1998 2006 2007 2008 2009 2010 2011 2012 TOTAL AUVERGNE 1.296 3.869 2.783 2.206 2.352 3.869 BOURGOGNE 0.920 3.793 - 2.394 , 1.852 1.193 2.398 2.000 3.793 2.083 CENTRE 4.379 1.716 2.071 2.292 1.896 1.140 1.174 1.834 4.379 CHAMP_NP 3.058 2.260 , 2.005 1.791 2.075 1.593 3.058 NORD-EST 2.020 1.452 0.901 , 1.520 0.685 1.425 POITOU-CHARENTES 3.178 3.987 . 2.841 1.785 1.780 1.063 3.987 1.360 RHONE-ALPES 1.492 3.869 3.550 3.451 2.179 2.376 2.335 0.772 3.869 SUD 1.942 3.789 2.859 2.298 2.532 2.306 4.168 3.344 3.290 4.168 TOTAL 3.178 4.379 3.550 3.451 2.352 2.859 2.841 2.532 2.306 4.168 3.344 3.290 4.379 Le tableau 9 représente la biomasse verte aérienne maximale de chaque combinaison « année x région » à la sortie de l'hiver (BDD_SH). - 35 - Les analyses statistiques (logiciel SAS) réalisées sur chaque base de données ont porté sur les pentes des relations entre la quantité d'azote absorbé dans les plantes entières et la biomasse fraîche aérienne. Elles sont les suivantes : - à l'entrée de l'hiver o analyse de variance et comparaison des moyennes par région toutes années confondues et par année toutes régions confondues par rapport à la moyenne générale ; o analyse de variance et comparaison entre-elles des moyennes par région, toutes années confondues et des moyennes par année, toutes régions confondues 10 (méthode des contrastes sur les régions ou sur les années 2 à 2 au seuil de 5 % et consolidation de proche en proche sur l'ensemble des régions ou des années) ; o analyse de variance et comparaison des moyennes par région par rapport à la moyenne générale pour chaque année ; o analyse de variance et comparaison entre-elles des moyennes par région 15 pour chaque année (méthode des contrastes sur les régions 2 à 2 par année au seuil de 5 % et consolidation de proche en proche sur l'ensemble des régions par année) ; - à la sortie de l'hiver : mêmes analyses qu'à l'entrée de l'hiver ; - à l'entrée et à la sortie de l'hiver, sur les situations où les mesures ont été réalisées à ces 2 périodes : 20 o analyse de variance et comparaison des moyennes des deux périodes de mesure entre-elles (entrée et sortie hiver) toutes années et régions confondues ; o analyse de variance et comparaison des moyennes des deux périodes de mesure entre-elles, par région, toutes années confondues ; o analyse de variance et comparaison des moyennes des deux périodes de 25 mesure entre-elles, par année, toutes régions confondues. L'hypothèse de non linéarité de la relation entre la biomasse verte aérienne et la quantité d'azote absorbé dans les plantes entières a également été étudiée en comparant les coefficients de détermination de deux types d'ajustement : - ajustement linéaire avec ordonnée à origine = 0 (Nabs=a.MVaer) ; 30 - ajustement quadratique avec ordonnées à origine = 0 (Nabs = a.MVaer2 + b.MVaer). 2.5.2.2 Comparaison des pentes de la relation Nabspe=f(MVaer) à l'entrée de l'hiver 35 Toutes années confondues, la relation entre la quantité d'azote absorbé dans les plantes entières et la biomasse fraîche aérienne à l'entrée de l'hiver présente des pentes différentes selon les régions : les valeurs vont de 54.4 pour le Sud à 62.1 pour - 36 - le Nord-Est (Figure 10). Ces différences sont souvent significatives au seuil de 5%. Ainsi, les valeurs obtenues dans le Nord-Est, le Centre et la Bourgogne sont significativement supérieure à la valeur moyenne France et aux valeurs obtenues en Rhône-Alpes, Poitou-Charentes et Sud (opposition Nord - Sud ; tableau 11). Pente Test F (5%) Pr> RI (écart par n (tests rapport à la « contraste moyenne Tous) 2 à 2» consolidés) Nord-Est 62,1 a 0,0664 * 13 Bourgogne 61,7 a 0,0416 ** 24 Centre 60,2 a 0,0171 ** 64 Champagne Nord Picardie 57,1 ab 0,844 NS 50 Rhône Alpes 55,7 b 0,4448 NS 47 Poitou Charente Vendée Limousin 54,5 b 0,1121 NS 54 Sud 54,4 b 0,0628 * 49 Tous 56,8 301 Le tableau 11 représente les tests de comparaison des pentes régionales des régression linéaires des relations exprimant la quantité d'azote absorbé dans les plantes entières et la biomasse fraîche aérienne à l'entrée de l'hiver, toutes années confondues (NS : non significatif au seuil de 10 % ; * : significatif au seuil de 10% ; ** : significatif au seuil de 5% ; ***: significatif au seuil de 1 % ; n : effectif). Les pentes sont également variables dans le cas d'une segmentation par année toutes régions confondues, les valeurs allant de 54.4 en 2007 à 62.4 en 2008 (Figure 11). La valeur moyenne toutes années et régions confondues est de 56.8. L'année 2008 se distingue significativement de 2006, 2007 et 2009 par sa valeur élevée (tableau 12). Pente Test F (5%) (tests Pr> iti (écart par rapport à la moyenne Tous) "contraste 2 à 2" consolidés) 2006 54.8 b 0.0903 * 70 2007 54.4 b 0.0466 ** 58 2008 62.4 a L'analyse annuelle ne met en évidence aucune différence significative entre régions en 2006 et 2007. En 2008, seule la région Centre se distingue par la valeur élevée de la pente alors qu'en 2009, c'est la région Poitou-Charentes qui se distingue par une pente faible (tableau 13). Il n'y a pas de stabilité interannuelle des valeurs régionales de pentes ou du classement entre région. Il y a au moins autant de variabilité entre années pour une région donnée qu'entre régions pour une année donnée (tableaux 13 et 14). Une analyse de variance réalisée sur la variable « pente » met d'ailleurs en évidence une forte interaction entre les effets années et régions : Pr. >F = 0.0002 (après élimination des combinaisons « région x année » présentant un effectif inférieur à 10). 2006 Pente Test F (5%) (tests Pr> RI (écart par à la moyenne rapport Tous) "contraste 2 à 2" consolidés) Poitou-Charentes, Vendée, Limousin 57.1 a 0.4548 NS 20 Sud 55.7 a 0.6826 NS 14 Centre 55.6 a 0.7964 NS 12 Champagne Nord Picardie 53.6 a 0.6377 NS 10 Rhône-Alpes 53.3 a 0.5370 NS 14 Tous 54.8 70 2007 Centre 55.5 a 0.6327 NS 10 Champagne Nord Picardie 55.2 a 0.6664 NS 13 Rhône-Alpes 53.8 a 0.8199 NS 14 Poitou-Charentes, Vendée, Limousin 52.5 a 0.569 NS 11 Sud 51.3 a 0.4094 NS 10 Tous 54.4 58 2008 Centre 70.4 a 0.0014 *** 16 Rhône-Alpes 62.4 b 0.9963 NS 11 Champagne Nord Picardie 60.5 b 0.395 NS 18 Sud 59.5 b 0.5202 NS 11 Poitou-Charentes, Vendée, Limousin 59.2 b 0.1316 NS 12 Tous 62.4 68 2009 Champagne Nord Picardie 60.7 a 0.0438 10 Bourgogne 56.3 a 0.4003 NS 14 Centre 54.8 a 0.4857 NS 21 Poitou-Charentes, Vendée, Limousin 46.0 b 0.0073 **. 12 Tous 53.2 57 Le tableau 13 représente les tests de comparaison des pentes régionales des régressions linéaires des relations exprimant la quantité d'azote absorbé dans les plantes entières et la biomasse fraîche aérienne à l'entrée de l'hiver, par année (NS : non significatif au seuil de 10%; * : significatif au seuil de 10%; ** : significatif au seuil de 5% ; *** : significatif au seuil de 1 ; n : effectif). - 38 - Centre Pente Test F (5%) (tests Pr> Itl (écart par à la moyenne rapport Tous) "contraste 2 à 2" consolidés) 2008 70.4 a 0.0003 HS 16 2006 55.6 b 0.5491 NS 12 2009 54.8 b 0.2664 NS 21 2007 52.5 b 0.1567 NS 11 Tous 57.6 60 Champagne-Nord-Picardie 2009 60.7 a 0.0754 10 2008 60.5 a 0.0226 ** 18 2006 53.6 b 0.1059 NS 10 2007 51.3 b 0.0436 ** 10 Tous 56.2 48 Poitou-Charentes-Limousin 2008 59.2 a 0.0657 12 2006 57.1 a 0.3828 NS 20 2007 55.5 a 0.6704 NS 10 2009 46.0 b 0.0017 *** 12 Tous 54.5 54 Rhône-Alpes 2008 62.4 a 0.0962 11 2007 55.2 ab 0.8537 NS 13 2006 53.3 b 0.4099 NS 14 Tous 55.7 38 Sud 59.5 a 0.4686 NS 11 2008 2006 55.7 a 0.8556 NS 14 2007 53.8 a 0.6107 NS 14 Tous 55.3 39 Le tableau 14 représente les tests de comparaison des pentes annuelles des régression linéaires des relations exprimant la quantité d'azote absorbé dans les plantes entières et la biomasse fraîche aérienne à l'entrée de l'hiver, par région (NS: non significatif au seuil de 10 %; *: significatif au seuil de 10%; ** : significatif au seuil de 5%; ***: significatif au seuil de 1 %; n : effectif). Tous ces résultats ne permettent donc pas de mettre en évidence des différences systématiques de pentes entre régions et entre années à l'entrée de l'hiver. 2.5.2.3 Comparaison des pentes de la relation Nabspe=f(MVaer) à la sortie de l'hiver Toutes années confondues, la relation entre la quantité d'azote absorbé dans les plantes entières et la biomasse fraîche aérienne à la sortie de l'hiver présente également des pentes différentes selon les régions : les valeurs vont de 65,1 en Poitou- Charentes à 76.3 en Auvergne (Figure 12). Ces différences sont souvent significatives au seuil de 5%. Ainsi, les valeurs obtenues en Auvergne sont significativement supérieures à toutes les autres ainsi qu'à la moyenne générale, alors que les valeurs issues de Rhône-Alpes et Poitou-Charentes sont significativement plus faibles que toutes les autres ainsi qu'à la moyenne (tableau 15). La pente moyenne toutes régions et années confondues est de 68,9. - 39 - Pente Test F (5%) Pr> Itj (écart par n (tests rapport à la « contraste moyenne Tous) 2 à 2» consolidés) Auvergne 76,2 a Iti (écart par rapport à la moyenne Tous "contraste 2 à 2" consolidés) 2008 72.7 a 0.0005 *** 93 1997 72.1 a 0.0013 *** 120 1996 72.0 a L'analyse annuelle met en évidence des différences significatives entre régions tous les ans sauf en 2009, année où il apparaît néanmoins de fortes différences au niveau des pentes régionales, mais où la variabilité au sein de chaque région est trop élevée pour mettre en évidence des différences statistiques (tableau 17). Le classement des régions n'est pas le même chaque année. Par exemple, la région Centre est en tête de classement avec des pentes élevées en 2008 et 2009 alors qu'elle est en queue de classement avec des pentes faibles en 2006 et 2007. L'Auvergne est souvent en haut de classement (1994, 1995, 1997), sauf en 1996, mêmes si ses valeurs sont toujours élevées. Il apparaît également des différences significatives entre années pour chaque région. Le classement des années est variable d'une région à l'autre. Certaines années se situent toutefois plus régulièrement en haut de classement que d'autres. Ainsi, 1997 présente souvent des valeurs élevées et 2007 des valeurs basses. Mais l'année 1995 par exemple présente un comportement beaucoup plus irrégulier (tableau 18). Il y a presque autant de variabilité entre années pour une région donnée qu'entre régions pour une année donnée (tableaux 17 et 18). Une analyse de variance réalisée sur la variable « pente » met d'ailleurs en évidence une forte interaction entre les effets années et régions : Pr. >F : Itl (écart par rapport à la moyenne Tous) "contraste 2 à 2" consolidés) Bourgogne 78.9 a (écart par rapport Ainm:grue "contraste 2 à 2" à la moyenne Tous) consolidés) 1997 85.3 a 0.0032 *** 22 1994 78.4 ab 0.5764 NS 19 1998 78.4 ab 0.3452 NS 19 1996 74.9 0.4444 NS 31 1995 73.3 0.117 NS 28 Tous 76.2 119 Bourgogne 1994 78.9 a 0.0614 23 1997 75.6 a 0.0194 ** 20 1996 75.0 a 0.0163 ** 44 2008 71.0 ab 0.8034 NS 11 2009 67.5 ab 0.5967 NS 15 1995 65.2 0.0014 *** 34 Tous 70.2 147 Centre 2008 83.7 a Iti (écart par rapport à la moyenne Tous n "contraste") EH 56.0 b iti (écart par rapport à la moyenne Tous) Test F (5%) (tests "contraste 2 à 2" consolidés) EH 54.6 b 0.0367 ** 64 SH 62.1 a 0.0085 *** 64 Tous 57.7 128 2007 EH 54.4 b 0.005 *** 58 SH 62.4 a 0.002 *** 58 Tous 58.1 116 2008 EH 62.2 b 0.1058 NS 66 SH 66.3 a 0.0971 * 66 Tous 64.2 132 2009 EH 53.9 b 0.0047 *** 60 SH 65.0 a iti (écart par rapport à la moyenne Tous) "contraste 2 à 2" consolidés) EH 59.7 b 0.0481 ** 22 SH 68.9 a 0.0126 ** 22 Tous 63.5 44 Centre EH 54.6 a 0.3136 NS 42 SH 58.8 a 0.2109 NS 42 Tous 56.3 84 Champagne-Nord-Picardie EH 56.2 b 0.0001 *** 40 SH 71.0 a Iti (écart par rapport à la moyenne Tous) Bourgogne - 2009 "contraste 2 à 2'' consolidés) EH 56.1 0.0327 ** 13 SH 68.7 a 0.0015 *** 13 Tous 60.5 26 Centre -2006 EH 55.6 a 0.8247 NS 12 SH 57.2 a 0.7940 NS 12 Tous 56.3 24 Centre - 2007 EH 52.5 a 0.9616 NS 11 SH 52.8 a 0.9580 NS 11 Tous 52.7 22 Centre -2009 EH 51.9 0.1825 NS 10 SH 66.2 a 0.0319 ** 10 Tous 56.5 20 Cham - a - ne-Nord-Picardie - 2007 EH 51.3 0.0348 ** 10 SH 65.7 a 0.0387 ** 10 Tous 58.6 20 Champagne-Nord-Picardie - 2008 EH 61.2 0.0172 ** 17 SH 70.7 a 0.0292 ** 17 Tous 66.3 34 Poitou-Charentes - 2006 EH 57.0 0.2284 NS 19 SH 64.1 a 0.1580 NS 19 Tous 60.1 38 Poitou-Charentes - 2007 EH 55.5 0.0996 10 SH 65.9 a 0.0789 10 Tous 60.5 20 Poitou-Charentes - 2008 EH 61.2 a 0.6602 NS 10 SH 64.3 a 0.5916 NS 10 Tous 62.5 20 Poitou-Charentes - 2009 EH 45.7 0.1086 NS 10 SH 57.4 a 0.0086 *** 10 Tous 49.3 20 Rhône-Alpes- 2006 EH 53.3 a 0.3956 NS 14 SH 61.1 a 0.2248 NS 14 Tous 56.2 28 Rhône-Alpes- 2007 EH 55.2 0.1618 NS 13 SH 63.7 a 0.0869 13 Tous 58.8 26 Rhône-Alpes- 2008 EH 62.3 a 0.7848 NS 10 SH 63.8 a 0.7861 NS 10 Tous 63.0 20 Sud - 2006 EH 55.8 a 0.5157 NS 13 SH 60.1 a 0.4900 NS 13 Tous 57.8 26 Sud - 2007 EH 53.8 a 0.3031 NS 14 SH 59.8 a 0.2965 NS 14 Tous 56.8 28 -46 - Le tableau 22, ci-dessus, représente des tests de comparaison des pentes des régressions linéaires des relations exprimant la quantité d'azote absorbé dans les plantes entières et la biomasse fraîche aérienne à l'entrée et à la sortie de l'hiver par situation « région x année » (NS : non significatif au seuil de 10 ; * : significatif au seuil de 10% ; ** : significatif au seuil de 5% ; *** : significatif au seuil de 1 ; n : effectif). 2.5.2.5 Discussion de l'hypothèse de linéarité de la relation Nabspe=f(MVaer) La question posée est la suivante : la relation entre le poids de matière fraîche aérienne et la quantité d'azote absorbé dans les plantes entières est-elle linéaire à l'entrée et à la sortie de l'hiver (coefficient de conversion constant quel que soit le niveau de biomasse) ou existe-il un phénomène de dilution de l'azote qui conduirait à une réduction du coefficient de conversion pour les biomasses élevées ? A l'entrée de l'hiver, l'effet de dilution n'apparaît pas. Les coefficients de 15 détermination des régressions linéaires et quadratiques sont très voisins (Figures 15A à 151). Il convient de préciser que dans notre base de données, les valeurs de poids de matière fraîche aérienne supérieurs à 3 kg/m2 sont rares. L'étude des résultats par région et année (Figures 16A à 16X) fait apparaître visuellement quelques situations qui pourraient être mieux représentées par une 20 régression quadratique (Centre 2007 et 2009), mais l'effet de dilution s'il existe n'est généré que par 1 ou 2 points. A la sortie de l'hiver, l'effet de dilution ne semble apparaitre que pour les régions Champagne-Nord-Picardie et Nord-Est. Dans ces deux régions, les biomasses maximales observées sont au plus de 3 kg/m2 en Champagne et de 2 kg/m2 dans le 25 Nord-Est. Pour les autres régions et pour la France entière, les coefficients de détermination des régressions linéaires et quadratiques sont très voisins malgré la présence de situations à très forte croissance (jusqu'à 4 kg/m2) (Figures 17 et Figures 18A à 181). L'étude des résultats par région et année (Figures 19A à 19AQ) montre que pour 30 Champagne-Nord-Picardie et le Nord-Est, l'effet de dilution est principalement causé par un effet année, 1995 présentant une pente plus faible que celle des autres années. Pour cette année 1995, les régressions linéaire et quadratique présentent d'ailleurs des coefficients de détermination comparables. Pour les autres régions, il n'apparaît pas d'effet de dilution à l'échelle annuelle. 35 Cette absence d'effet de dilution pour les biomasses élevées est confirmée par le positionnement des données sur un graphique qui met en relation la teneur en azote et la biomasse sèche (Figures 20 et 21) : en moyenne, la teneur en azote est stable dans -47 - la gamme de biomasses étudiée. En revanche, la courbe critique de dilution N, qui représente le %N minimum pour que cet élément ne soit pas limitant pour la croissance, présente comme son nom l'indique une diminution de la teneur en azote quand la biomasse augmente (Colnenne, 1998). La stabilité de la teneur en azote, qui 5 est observée en moyenne dans nos données, indique que les petits colzas sont davantage carencés en azote que les gros (cf. la plus grande proportion de points en dessous de la courbe critique de dilution N). Cela suggère que les petits colzas sont en partie le résultat d'une faible disponibilité en azote qui conduit à des carences et à des faibles biomasses. C'est pourquoi on observe une corrélation linéaire entre la quantité 10 d'azote dans les plantes et la biomasse. Les tendances sont similaires en cas de segmentation des bases de données par région (non présenté). 3.5.2.6 Règle de calcul de Pi En l'absence de perte de feuilles vertes due au gel pendant l'hiver, le formalisme de 15 calcul de Pi de la réglette 1998 est conservé. Seul le poids de matière fraîche aérienne à la sortie de l'hiver (MV SH) est mobilisé et il est multiplié par le coefficient de conversion de la biomasse aérienne fraîche en quantité d'azote absorbé dans les plantes entières spécifique à la sortie de l'hiver (coeffSH). En cas de perte de feuilles vertes due au gel pendant l'hiver, les inventeurs 20 proposent une modification du calcul du Pi par rapport à la réglette 1998. La règle de décision de la réglette 1998 consiste à comparer les poids frais aériens à l'entrée et à la sortie de l'hiver. En cas de diminution du poids frais aérien pendant l'hiver, la quantité correspondante d'azote perdu par la culture (dNabs) est estimée en multipliant l'écart de poids frais par le coefficient de conversion de la biomasse fraîche 25 aérienne et quantité d'azote absorbé dans les plantes entières déterminé à la sortie de l'hiver (coeff). La quantité d'azote issue de cette perte de feuille réabsorbée par la culture après l'ouverture du bilan est estimée à 50 % de l'azote perdu soit Nabs x 0.5. Elle est donc ajoutée à la quantité d'azote présent dans la culture à l'ouverture du bilan, elle-même calculée en multipliant le poids frais aérien à la sortie de l'hiver et le 30 coefficient de conversion (estimé à la sortie de l'hiver). Le formalisme est donc le suivant : si MV EH > MV_SH, alors Pi = (MV_SH x coeff) + % [(MV_EH - MV_SH) x coeff] = (MV_SH + MV_EH) /2 x coeff. 35 Cette règle de décision repose sur l'hypothèse que le coefficient de conversion est le même à l'entrée et à la sortie de l'hiver. De plus, comme ce coefficient de conversion a été calculé en intégrant la partie racinaire, on considère qu'une partie de l'azote perdu - 48 - provient des racines, ce qui est faux car on n'observe jamais de décroissance du poids de matière fraîche, du poids de matière sèche et de la quantité d'azote des racines pendant l'hiver (résultats des essais APN non présentés).. Les inventeurs ont démontré que coeffEH est significativement plus faible que 5 coeffSH. Dans la méthode selon l'invention, si les inventeurs avaient continué à comparer les poids frais aériens, il y aurait eu des situations où le poids frais aérien aurait diminué pendant l'hiver sans qu'il y ait eu pour autant d'augmentation de Pi lié à la réabsorption après l'ouverture du bilan de l'azote perdu pendant l'hiver, ce qui est difficilement justifiable (exemple : en faisant l'hypothèse que MV EH = 3 kg/m2, MV SH 10 = 2.5 kg/m2, coeffEH = 57 et coeffSH = 69 (coeffEH et coeffSH = pentes des droites de régression) alors NabsEH = 171 kg/ha et NabsSH = 173 kg/ha. Faire la moyenne de NabsEH et NabsSH ne changerait pratiquement rien par rapport à NabsSH seul alors que la perte de biomasse fraîche serait de 0.5 kg/m2). Il convient donc de proposer une écriture qui permette de n'affecter la diminution de 15 biomasse qu'aux parties aériennes et de provoquer une augmentation de Pi en cas de perte de feuilles vertes pendant l'hiver, tout en ne mobilisant pas d'autre variable observée que le poids de matière fraîche aérienne. Il est tout d'abord présenté la formulation qui aurait été la plus rigoureuse pour estimer la perte d'azote liée au gel des feuilles vertes (dNabs en kg/ha), sachant que 20 les inventeurs ne l'ont pas retenue finalement car ils ont considéré qu'elle était trop complexe pour être facilement vulgarisée. Cette formulation serait la suivante : dNabs = ((MSa EH x %Na EH /100) - (MSa SH x %Na SH /100)) x 10000 / 1000 25 équation 1; avec MSa = poids de matière sèche aérienne (g/m2) et %Na = teneur en azote de la matière sèche des parties aériennes (en %) L'analyse de la base de données « BDD_EH_SH » décrite au chapitre 2.5.2.1 30 (parcelles avec mesures de biomasse fraîche, de biomasse sèche et de teneur en azote à l'entrée et à la sortie de l'hiver) montre que la teneur en azote moyenne de la matière sèche aérienne est la même (pas de différence significative) à l'entrée et à la sortie de l'hiver et que sa valeur est de 3.02 % (moyenne arithmétique de toutes les valeurs de teneur en azote en entrée et sortie d'hiver). Ceci est illustré en Figure 22A 35 et 22B. L'équation 1 devient donc : dNabs = (MSa EH - MSa SH) x 3.02 / 10 équation 2. -49 - Il reste à exprimer MSa EH et MSaSH en fonction les variables observées MV EH et MV SH en kg/m2. L'analyse de BDD EH_SH permet de déterminer une teneur en matière sèche moyenne de la matière fraîche à l'entrée (14.5 %) et à la sortie de l'hiver (16.2 %). Ces valeurs sont obtenues en réalisant la moyenne arithmétique des valeurs unitaires de teneur en matière sèche. Elles sont significativement différentes l'une de l'autre. L'équation 2 devient donc : dNabs = ((MV EH x 14.5 / 100) -(MV SH x 16.2 / 100)) x 3.02 x 100 équation 3. Du fait de l'absence de données nouvelles permettant de ré-estimer la part de l'azote issu des feuilles vertes tombées pendant l'hiver et réabsorbé par le colza en place au printemps, la valeur de la réglette 1998 a été conservée, soit 50 %. Cette valeur est issue des travaux de Dejoux (thèse 1999) qui indique que « le colza a bénéficié au cours de sa croissance de printemps d'environ 40 % de l'azote apporté par les feuilles gelées ». Ce CAU de 40 % conduit à une équivalence engrais N minéral de 50 % pour l'azote des feuilles. L'équivalence engrais N minéral d'un produit organique est en effet égale au rapport entre le CAU du produit organique et le CAU de l'engrais N minéral. Si l'on considère que le CAU de l'engrais N minéral est en moyenne de 80 c/o, on obtient en effet une équivalence de 50 % (40 / 80). Comme dans la réglette 1998, cette valeur de 50% est calculée en réalisant la moyenne des quantités d'azote présent dans les plantes à l'entrée et à la sortie de l'hiver. On fait alors l'hypothèse que la croissance des plantes a été négligeable pendant l'hiver, donc que la différence de biomasse est totalement imputable aux pertes de feuilles. Ceci est une approximation, car il est fort probable que dans de nombreux cas la croissance ait été significative. Le formalisme de calcul de Pi en cas de perte de feuilles vertes due au gel est donc le suivant : Si MSaEH > MSaSH, alors Pi = (MV SH x coeff SH) + % [((MV EH x 14.5 / 100) - (MV SH x 16.2 / 100)) x 3.02 x 100]. Le test d'égalité des biomasses ne porterait plus sur la biomasse fraîche comme dans la réglette 1998, mais sur la biomasse sèche aérienne (calculée à partir de la biomasse fraîche). Ce serait plus rigoureux et cela réglerait les problèmes de diminution de biomasse fraîche non accompagnée d'augmentation du Pi que l'on rencontre à cause de la différence entre coeffEH et coeffSH lorsque le test de - 50 - comparaison est réalisé sur les biomasses fraîches. Ce formalisme présente en revanche l'inconvénient de mobiliser des paramètres issus de bases de données différentes à la sortie de l'hiver (coeff_SH issu de la base de données BDD_SH complète et teneur en matière sèche et teneur en azote de la matière sèche issue de BDD_EH_SH). Pour une raison de simplification de la règle de décision pour la rendre plus facilement vulgarisable, les inventeurs ont finalement opté pour un formalisme consistant à comparer non pas les biomasses sèches, mais les quantités d'azote absorbé à l'entrée et la sortie de l'hiver calculées à partir du poids de matière fraîche aérienne à l'entrée et à la sortie de l'hiver et des coefficients de conversion déterminés à l'entrée et à la sortie de l'hiver : - Nabs_EH (plantes entières) = MVaer_EH x coeff EH - Nabs_SH (plante entières) = MVaer_SH x coeff_SH. Comparer Nabs EH et SH sur plantes entières revient à comparer les Nabs EH et SH sur parties aériennes puisqu'ils sont proportionnels : - Nabs_EH = Nabs_EH aérien x 1.35 - Nabs_SH = Nabs_SH aérien x 1.35. Cela revient également à comparer les poids de matière sèche aérienne si l'on fait l'hypothèse que les teneurs moyennes en azote de la matière sèche aérienne sont équivalentes à l'entrée et à la sortie de l'hiver comme cela a été montré précédemment sur la base de données BDD_EH_SH. Le formalisme devient donc : Si Nabs_EH > Nabs_SH alors dNabs = (Nabs_EH - Nabs_SH) /1.35 La diminution de quantité d'azote dans la culture pendant l'hiver (dNabs) n'est calculée que sur la base des parties aériennes grâce à la division par 1.35 de l'écart d'absorption sur plantes entières. Le calcul de Pi est donc le suivant : Pi = Nabs_SH + 1/2 (dNabs) = Nabs_SH + (1/2 30 (Nabs_EH - NabsSH) / 1.35). 2.5.2.7 Détermination des valeurs de coeff_EH et coeff_SH Au regard des éléments précédents, les inventeurs proposent de ne plus avoir qu'une valeur unique nationale pour chacun des coefficients de conversion. 35 Plutôt que de choisir d'emblée les valeurs de ces coefficients correspondant aux pentes des régressions linéaires obtenues sur nos jeux de données, les inventeurs ont fait le choix, dans un premier temps, de proposer différentes valeurs correspondant à - 51 - différents déciles (1, 3, 5, 7 et 9), c'est-à-dire à différents niveaux de risque de surévaluer ou sous-évaluer la quantité d'azote absorbé. Par exemple, le coefficient dit « décile 1 » revient à prendre le risque de sous-estimer de la quantité d'azote absorbé dans 9 cas sur 10. Le choix final des valeurs sera réalisé lors de l'évaluation globale de la règle de décision présentée plus loin dans cette synthèse technique. Les valeurs retenues pour l'évaluation sont présentées sur les Figures 23 et 24. Elles correspondent aux pentes pour lesquelles 10%, 30%, 50%, 70% ou 90% des points sont situées sous la droite (calcul avec Excel 2010). 10 2.6 Quantité d'azote minéral dans le sol à l'ouverture du bilan (Ri) 2.6.1 Matériel et Méthodes Trois bases de données de reliquats en sortie d'hiver sous colza mises à disposition par nos partenaires sont mobilisées : - base de données de SAS Laboratoire : références acquises de 2009 à 2013 sur 15 4011 parcelles sur tous types de sol (dont 119 en craies de Champagne) ; - base de données de la Chambre d'Agriculture d'Eure-et-Loir (CA28) : références acquises de 2007 à 2013 sur 1406 parcelles ; - base de données de la CAMA (51) : références acquises de 2005 à 2014 sur 1050 parcelles en sols de craie de Champagne. 20 Les mesures ont été réalisées sur 1 à 3 couches de 20 à 30 cm (le plus souvent 30 cm). Le calcul de la quantité d'azote minéral (Ri) est le suivant : somme N-NO3- sur les couches mesurées + N-NH4+ sur couche 1 + 1/3 N-NH4+ de la couche 2. Les bases de données sont nettoyées selon les critères suivants : 25 - suppression des parcelles avec une quantité de N- NH4 + supérieure à 20 kgN/ha sur la couche 1; - suppression des parcelles avec une quantité d'azote minéral total supérieure à 250 kg N/ha. La distribution des Ri ainsi calculés dans chacune des bases est présentée sur les 30 Figures 25A à 25C (logiciel SAS). La base de données du SAS a été utilisée pour estimer des valeurs forfaitaires de reliquats sur les sols profonds et superficiels. La base de la CA28 a été utilisée pour évaluer la pertinence de ces valeurs. La base de données du SAS ne couvre pas l'intégralité du territoire national mais les 35 principales zones de production de colza sont représentées (moitié Nord de la France). Le Sud de la France est probablement sous-représenté mais les références sous colza sont rares, y compris dans les organismes de développement (Figure 26). - 52 - Il convient de signaler que le nombre de couches prélevées est variable d'une parcelle à l'autre. Cela a un impact sur le reliquat d'azote minéral total qui est calculé à partir des couches ayant fait l'objet d'une analyse. Pour la base de données du SAS, on retrouve logiquement une faible proportion de parcelles avec 3 couches prélevées 5 parmi les sols argilo-calcaires superficiels. Dans les sols de limons, seulement deux tiers des parcelles ont fait l'objet d'un prélèvement sur 3 couches et 10 % n'ont fait l'objet d'un prélèvement que sur une couche (tableau 23). Dans la base de données de la CA28, aucun prélèvement n'a été réalisé que sur une couche. Dans la base de données de la CAMA, 31 % des parcelles ont fait l'objet d'une analyse sur 2 couches et 10 69 % sur 3 couches. SAS Nombre de couches prélevées Nombre de parcelles 1 2 I 3 Types de sois % de parcelles Argiles 18 62 20 254 Argilo-calcaires sup. 47 46 7 921 Craies 12 56 32 119 Limons 12 51 37 1997 Sables 18 54 28 384 Autres sols 9 63 28 65 CA 28 Nombre de couches prélevées Nombre de parcelles 2 3 Types de sols % de parcelles Argiles 0 29 71 7 Argilo-calcaires 0 84 16 89 Limons 0 30 70 1182 Sables 0 36 64 22 Sols caillouteux 0 86 14 105 Le tableau 23 résume la distribution des parcelles en fonction du nombre de 15 couches prélevées et du type de sol (base de données du SAS en haut et de la CA28 en bas). 2.6.2 Résultats 2.6.2.1 Tous types de sols hors craie de Champagne Le traitement des données fournies par le SAS (hors sols de craie de Champagne) 20 met en évidence un effet significatif régional sur les moyennes de reliquat d'azote minéral : les valeurs obtenues dans l'Ouest de la France sont significativement supérieures (certainement en raison des systèmes de production) à celles obtenues dans les autres régions qui ne se distinguent pas entre-elles. Les inventeurs proposent - 53 - un découpage du territoire national en 2 zones : - Grand Ouest (région Bretagne, région Pays de la Loire hors Vendée et région Basse-Normandie hors Calvados) - reste France. Les principales statistiques calculées pour ces deux zones sont présentées dans les tableaux 24 et 25. Type de sol Nombre Moyenne Médiane .10% '5% 75% 90% Profond 104 45 39 22 30 51 78 Superficiel 2! 28 28 26 26 30 30 Le tableau 24 représente les principales caractéristiques statistiques du reliquat à 10 l'ouverture du bilan (kgN/ha) pour le « grand Ouest » (base de données du SAS). Type de sol Nombre Moyenne Médiane 10% 25% 75% 90% Profond 1922 35 31 17 23 41 57 Superficiel 426 22 19 10 14 25 35 Le tableau 25 représente les principales caractéristiques statistiques du reliquat à l'ouverture du bilan (kgN/ha) pour le « grand Ouest » (base de données du SAS). 15 Les définitions retenues pour les catégories de sol sont les suivantes : - sols « profonds » : parcelles avec 2 ou 3 couches prélevées en excluant les sols argilo-calcaire superficiels, les craies et les sables ; - sols «superficiels» : sols dénommés argilo-calcaire superficiel avec une couche 20 prélevée. Les moyennes étant fortement influencées par les valeurs hautes des échantillons, il a semblé préférable aux inventeurs de retenir les médianes. Les valeurs forfaitaires retenues sont présentées dans le tableau 26. Ri forfaitaire (kgN/ha Grand Ouest Reste France Sols profonds 40 30 Sols superficiels 30* 20 25 Le tableau 26 représente les valeurs forfaitaires de Ri (kgN/ha) selon la situation géographique et le type de sol (base : médiane). * peu de références : valeur estimée 30 (médiane sols profonds -10 kgN/ha). Le choix de la médiane sous-entend une certaine prise de risque de surestimer 30 (50% des parcelles présentent des Ri inférieurs) ou sous-estimer le reliquat (50% des parcelles présentent des Ri supérieurs). Au regard de la forte variabilité interannuelle - 54 - associée aux valeurs de Ri, ce choix n'est pas très « sécuritaire ». On pourrait envisager de réviser à la hausse ou à la baisse ces valeurs pour limiter les risques de surfertiliser ou sous-fertiliser. Par exemple, pour limiter le risque de sous-fertilisation, les inventeurs pourraient retenir les valeurs du premier quartile (25% des valeurs supérieures). Cela conduirait aux valeurs forfaitaires présentées dans le tableau 27. Le tableau 27 représente les valeurs forfaitaires de Ri (kgN/ha) selon la situation géographique et le type de sol (base : premier quartile). Plusieurs hypothèses et incertitudes entourent le traitement de la base de données du SAS : les sols bien pourvus sont-ils sur-représentés ? Les critères de nettoyage de la base de données sont-ils suffisants pour éliminer les valeurs aberrantes ? Pour conforter ou nuancer les résultats présentés ci-dessus, les inventeurs se sont livrés au même traitement sur la base de données de la CA28. Il apparaît que les valeurs médianes obtenues en sols superficiels et profonds sont plus faibles que celles issues de la base de données SAS, zone « reste France » (tableau 28). Elles correspondent davantage aux valeurs du premier quartile de cette base de données. Types de sols Effectif Moyenne Médiane 10% 25% 75% 90% Argiles 7 26 24 3 7 50 60 Argilo-calcaires 89 24 17 9 13 33 44 Limons 1182 25 20 9 13 32 48 Sables 22 23 17 9 11 23 48 Sols caillouteux 105 20 17 7 12 24 35 Le tableau 28 représente les principales caractéristiques du reliquat à l'ouverture du bilan (kgN/ha) pour la base de données de la CA28. Les inventeurs considèrent cependant qu'il n'est pas nécessaire de chercher à prendre une marge de sécurité sur ce paramètre, car l'utilisateur pourra substituer une valeur mesurée (mesurée sur la parcelle ou issue d'un réseau régional) à la valeur forfaitaire. L'hypothèse d'un effet d'apports réguliers de produits résiduaires organiques (PRO) sur le Ri n'est pas vérifiée. Les parcelles avec apport de PRO ne présentent des valeurs supérieures que d'environ 3 kgN/ha à celles des autres parcelles (tableau 29) 30 sur le critère de la moyenne ajustée. Le paramétrage n'en tient donc pas compte. Ri forfaitaire (kgN/ha Grand Ouest Reste France Sols profonds Sols superficiels 25 20* 15 - 55 - Test Apport SNK régulier (seul en Effectif Moyenne Médiane 10% 25% 75% 90% 5%) sur PRO moyenne ajoutée NON 2563 32 28 14 19 39 54 a OUI 1176 35 31 16 22 43 58 a Le tableau 29 représente les principales statistiques liées à l'effet des apports réguliers de PRO sur le reliquat d'azote minéral à l'ouverture du bilan (kgN/ha) (base de données du SAS ; toutes situations). L'effet des apports de PRO juste avant la culture de colza est quant à lui significatif, mais le surplus de reliquat associé est faible (tableau 30). Ce facteur n'a donc pas non plus été pris en compte dans e paramétrage. Apport Effectif Moyenne Médiane 10% 25% 75% 90% Test de SNK PRO (seul 5%) sur moyenne ajoutée sur le colza NON 2456 31 27 14 19 37 50 a OUI 1284 38 32 17 22 45 65 b Le tableau 30 représente les principales statistiques liées à l'effet de I apport de PRO sur le colza PRO sur le reliquat d'azote minéral à l'ouverture du bilan (kgN/ha) 10 (base de données du SAS ; toutes situations). 2.6.2.2 Sols de craie de Champagne La base de données utilisée regroupe toutes les références de la CAMA et les 119 références en sols de craie de Champagne issues du SAS, quel que soit le nombre de 15 couches. L'effet des apports de PRO sur le colza avant l'ouverture du bilan sur le reliquat d'azote minéral du sol médian est très faible (+ 4 kgN/ha). Les inventeurs ont donc choisi de ne pas l'introduire dans notre règle de décision. 20 La base de données ne permet pas de tester l'effet de l'historique d'apport de PRO. - 56 - Cet éventuel effet ne sera donc pas pris en compte non plus. Les principales caractéristiques statistiques de la base de données sont présentées dans le tableau 31. Pour les raisons présentée précédemment, les inventeurs choisiront de retenir la médiane plutôt que la moyenne ou le premier quartile. La valeur 5 forfaitaire retenue pour le Ri en sols de craie de Champagne est donc de 40 kgN/ha. Nombre Moyenne Médiane 10% 25% 75% 90% 1169 47 41 23 30 56 79 Le tableau 31 représente les principales caractéristiques du reliquat à 'ouverture du bilan (kgN/ha) pour les sols de craie de Champagne (base de données CAMA + SAS). 10 3.6.2.3 Conclusion Le paramétrage choisi pour l'étape d'évaluation est celui présenté dans le tableau 32, sur la base de la médiane de la base de données du SAS pour deux zones (zone Ouest et reste de la France) sans effet spécifique de l'apport de PRO (régulier ou sur la 15 culture). Méthode de l'invention Région Sols superficiels Sols profonds Sols de craie de Champagne Grand Ouest 30 40 S.S. (Bretagne, Pays de la Loire hors Vendée et Basse Normandie hors Calvados) Champagne 20 30 40 Ardenne, Aines Reste de la France 20 30 s.s. Réglette 1998 Région Sols superficiels Sols profonds Bretagne 20 30 Centre 10 20 Reste de la France 10 30 Le tableau 32 représente les valeurs forfaitaires de quantité d'azote minéral dans le sol à l'ouverture du bilan (Ri) proposées pour la réglette azote 2014 (en haut) et rappel des valeurs retenues dans la réglette 1998 (en bas) ; valeurs exprimées en kgN/ha. 20 2.7 Minéralisation de l'humus et des résidus de récolte (M) Ce poste ne concerne que les régions en écriture additive du bilan. En l'absence de références nouvelles, les inventeurs ont fait le choix de reprendre le référentiel des - 57 - valeurs de minéralisation nette (Mn) de la réglette 1998, en simplifiant le découpage régional de façon à le rendre homogène avec celui proposé pour les reliquats d'azote minéral du sol à l'ouverture du bilan. Le précédent du colza est une céréale à paille sur plus de 92 % des surfaces françaises (enquête CETIOM 2012). Les résidus de culture du précédent sont ramassés sur 63 % des surfaces, enfouis sur 29 % et laissés en surface sur 7 %. Dans la base de données d'évaluation, les pailles sont enfouies sur 57% des parcelles et enlevées sur 43 % des parcelles. Les valeurs de minéralisation proposées correspondent aux situations avec comme précédent une céréale à paille quel que soit le mode de gestion des pailles. La valeur de 20 kgN/ha annoncée comme effet de la gestion des pailles devrait en effet être modulée en fonction des conditions pédoclimatiques, ce qui n'est pas compatible avec notre objectif d'outil simple, et des simulations avec AZOFERT de l'effet des conditions pédoclimatiques montrent que la valeur moyenne de 20 kgN/ha est surestimée. Les inventeurs ont donc choisi de ne pas prendre en compte cet effet. Pour respecter leur choix de simplicité du formalisme de la réglette et pour pouvoir prendre en compte les apports très fréquents de produits résiduaires organiques, les inventeurs ont choisi de ne pas reprendre le formalisme du COMIFER faisant appel à un coefficient correcteur de la minéralisation (Fsyst) fonction du système de culture (gestion des résidus de récolte, type et fréquence des apports de PRO). L'historique d'apport des PRO sera pris en compte au travers d'un surplus de minéralisation (paragraphe 2.10). Les valeurs forfaitaires de M sont présentées dans le tableau 33 qui rappelle également les valeurs utilisées dans la réglette azote 1998. Méthode de l'invention Région Sols superficiels Sols profonds Sols de craie de Champagne Grand Ouest 30 50 s.s. Champagne Ardenne, Aines 20 40 20 Reste de la France 20 40 s.s. Réglette 1998 Région Sols superficiels Sols profonds Bretagne 30 50 Reste de la France 20 40 Le tableau 33 représente les valeurs forfaitaires de minéralisation de l'humus et des résidus de culture (M) proposés pour la réglette azote 2014 (à gauche) et valeurs de Mn de la réglette 1998 (à droite) ; valeurs exprimées en kgN/ha. 3.8 Fourniture d'azote par le sol (PO) Ce poste ne concerne que les régions en écriture « CAU » du bilan. - 58 - En l'absence de références nouvelles, les inventeurs ont fait le choix de reprendre le référentiel des valeurs de fourniture d'azote par le sol (Np) de la réglette 1998. Les valeurs proposées (tableau 34) valent pour tous les modes de gestion des pailles du précédent confondus (hors précédent légumineuse). Région Sols Sols superficiels profonds Auvergne et Rhône-Alpes 40 60 Nord-Est 30 50 Le tableau 34 représente les valeurs forfaitaires de fourniture d'azote par le sol (PO) proposés pour la réglette azote 2014 ; valeurs exprimées en kgN/ha. 2.9 Fourniture d'azote issu des PRO épandus en fin d'été ou à l'automne juste 10 avant ou sur le colza en place (Mpro1) Une des avancées majeure attendue dans la Réglette 2014 consiste à mieux prendre en compte les apports de produits organiques sur colza, le plus souvent appliqués à la fin de l'été ou à l'automne, dans le calcul de la fertilisation azotée de 15 printemps. A cette fin, Les inventeurs ont décidé d'établir une table rassemblant les éléments nécessaires au calcul par défaut de la fourniture d'azote issu de ces apports. 2.9.1 Matériel et méthodes 20 Plusieurs sources d'information ont été consultées : les références du COMIFER, du CORPEN, les arrêtés préfectoraux parus à l'été 2012 qui déclinent les référentiels agronomiques faisant foi dans les zones vulnérables des principales régions françaises d'élevage (pour l'Ouest : la Haute-Normandie, la Basse-Normandie et les Pays de la Loire ; pour l'Est : la Lorraine, la Bourgogne, la Franche-Comté et aussi la Champagne 25 Ardenne et la Picardie pour des références spécifiques ; pour le Sud : Midi-Pyrénées), le Guide des Matières Organiques Tome 2 édité par l'ITAB en 2001 et la base de données Arvalis-Fertiweb en cours de construction (juillet 2013). Les références se comptent par dizaines voire centaine en fonction du type d'animaux et du type d'exploitation. Cela s'accompagne d'une forte divergence des 30 valeurs de référence proposées pour les teneurs en azote des produits (kg N / t ou m3) et surtout pour les coefficients d'équivalence engrais minéral (Keq). En conséquence, dans un souci de faciliter l'usage du futur outil d'aide à la décision et d'apporter une information pertinente pour le calcul, les inventeurs proposent de produire une table par défaut restreinte (20 produits décrits + 1 référence vierge) et de5 - 59 - laisser à l'utilisateur la possibilité de modifier la teneur en azote et la quantité de produit apporté. Cela se justifie d'autant plus qu'il existe une forte variabilité dans la composition des produits et dans les pratiques d'épandage. Les PRO référencés dans la table sont ceux qui sont le plus fréquemment appliqués sur colza. Ils ont été identifiés à l'aide d'une analyse de représentativité (fréquence d'utilisation) sur les bases de données du SAS, de la CA28 (bases de données des reliquats d'azote minéral du sol à l'ouverture du bilan sous colza) et de l'enquête SSP 2006 (annexe 4). Ils couvrent environ 85% des pratiques d'apport. Les quantités moyennes apportées ont été extrapolées de ces mêmes bases de données. Afin d'utiliser des références de quantité d'azote par unité de produit brut et de coefficient d'équivalence engrais reconnues, Les inventeurs ont toujours privilégié les valeurs publiées par le COMIFER (mai 2013). Lorsque les références n'étaient pas disponibles dans les publications du COMIFER, les inventeurs ont choisi des valeurs issues de la base de données d'Arvalis ou de référentiels régionaux qui leur semblaient pertinents et cohérents. 2.9.2 Résultats Le tableau 35 présente la table proposée par défaut. Deux mises en alerte devront par ailleurs être exercées vis-à-vis de l'utilisateur : - compte-tenu de la forte variabilité des valeurs observées autour des compositions moyennes des produits, il est toujours préférable de disposer de mesures réalisées sur le produit épandu ; - en cas d'apport de produits organiques, la mesure du reliquat d'azote minéral 25 dans le sol à la sortie de l'hiver est fortement conseillée. Les valeurs de Keq correspondent à la fourniture d'azote issue du PRO pendant la période d'ouverture du bilan et non pas sur tout le cycle. On considère que la partie minéralisée avant l'ouverture du bilan est le plus souvent absorbée par la culture et donc prise en compte dans le poste Pi (quantité d'azote absorbé par la culture à 30 l'ouverture du bilan). Cette hypothèse justifie toutefois le conseil de mesure du reliquat d'azote minéral du sol à l'ouverture du bilan dans ces situations d'apport de PRO. Le calcul du poste Mpro1 est effectué de la façon suivante : Mpro1 (kgN/ha) = quantité PRO apporté (t ou m3 / ha) x teneur N total PRO (kg/t ou kg/m3) x Keq. Produit résiduaire organique Teneur en azote total (Kg N/t ou Quantité Unité Coefficient d'équivalence engrais (Keq) sur la période du bilan colza Hypothèse produits, période et condition d'application du colza m3) de produit (t ou épendue m3/ha) par hectare Fumier - de bovins 5,8 20 t/ha 0,1 Fumier de bovin pailleux de litières accumulée, fin été - de volailles 25 5 t/ha 0,2 Fin été, incorporation immédiate - de porcs 8 20 t/ha 0,1 Fin été - de de cheval 8 20 t/ha 0,1 Fin été - d'ovins 7 20 t/ha 0,1 Fin été Lisiers - de porcs à l'engrais 5,8 * 20 m3/ha 0,05 Fin été, incorporation immédiate - de porcs mixtes 3,5 30 m3/ha 0,05 Fin été, incorporation immédiate - de bovins 4,5 30 m3/ha 0,15 Fin été, incorporation immédiate - de bovins dilué 1,6 30 m3/ha 0,15 Fin été, incorporation immédiate - de canard 7,7 * 20 m3/ha 0,05 ** Fientes - de volailles sèches (80%MS) 40 5 t/ha 0,05 ** - de volailles (60%MS) 24 5 t/ha 0,05 ** Fientes de volailles 60%MS Boues urbaines - liquides 3* 30 m3/ha 0,15 * Boues urbaines liquides biologiques - pâteuses (20%MS) 10 * 10 t/ha 0,1 ** Boues pâteuses physico-chimiques ou biologiques (20%MS) - séchées (90%MS) 40 * 10 t/ha 0,1* Composts - de déchets verts 10 10 t/ha 0,05 Fin été, déchet vert de plus de 6 mois - du fumier de bovins 6,3 15 t/ha 0,12 Fin été, compost de fumier de bovins jeune de moins de 6 mois - de fumier de volailles 23 3 t/ha 0,12 Fin été Autres Vinasse de betteraves concentrées 20 3 t/ha 0,15 Fin été Ecumes de sucrerie supressées 4 * 10 t/ha 0** Autre produit « champ libre » t ou 0,1 Fin été m3/ha Tableau 35 représente les valeurs forfaitaires de teneur en azote total sur produit brut, de quantité de produit brut épandu et de coefficient d'équivalence engrais minéral pour les 20 produits les plus épandus avant colza. Encadré : Expertise, synthése des informations déclarées lors des analyses de reliquat s sortie d'hiver sous colza (bases SAS et CA28) ; * : ARVALIS, base PRO juillet 2013 ; ** : Arrêté du 20 juillet 2012 établissant le référentiel régional de mise en oeuvre de l'équilibre de la fertilisation azotée en Champagne-Ardennes ; Reste : COMIFER - Groupe azote mai 2013. -61- 2.10 Supplément de fourniture d'azote lié à l'historique d'apport de PRO (Mha) La règle de décision de la Réglette azote 1998 tient compte du surplus de minéralisation engendré par des apports organiques réguliers sur la parcelle sans distinction de produit ou de fréquence d'apport. Dès lors que la parcelle reçoit au moins tous les 3 ou 4 ans, depuis une longue période, des effluents d'élevage, une réduction forfaitaire de 40 kgN/ha (50 kgN/ha en sol profond en régions Bretagne et Pays de la Loire) est appliquée à la dose d'azote conseillée. Or, d'après des travaux récents publiés par le COMIFER, il apparait que le surplus 10 de minéralisation pris en compte dans la Réglette 1998 est nettement surestimé pour des pratiques d'apport moyennes. Pour preuve, un surplus de minéralisation de l'ordre de 40 kg N/ha/an correspond, en tendance, à un apport de 20 à 25 t/ha de fumier de bovins, chaque année, pendant une dizaine d'années (Arvalis, lnra Rennes, CRA Bretagne dans Perspectives Agricoles nc392 - septem bre 2012). De plus, l'incidence 15 des apports réguliers de PRO est variable selon la quantité apportée, la nature du produit et la fréquence des apports sur la parcelle. L'arrière-effet des PRO ne peut donc pas se limiter à une valeur unique. La révision de l'effet à long terme des apports de produits résiduaires organiques constitue donc une voie de progrès importante pour la Réglette 2014. 20 2.10.1 Matériel et méthodes L'impact à long terme des apports organiques exogènes est pris en compte de manière très variable dans les méthodes et les outils de calcul de fertilisation. Deux méthodes font aujourd'hui référence en France : le coefficient « facteur système » 25 (FSyst) proposé par le COMIFER et le « référentiel breton » sur la contribution des apports historiques de produits organiques, largement utilisé dans l'Ouest de la France. Méthode dite du FSyst Dans la méthode de calcul proposée par le COMIFER (version 2011), les effets à long terme des PRO sont pris en compte et estimés à partir de la minéralisation du 30 stock d'azote organique du sol modulé par un coefficient facteur système (FSyst). Ce coefficient est utilisé pour moduler la vitesse potentielle de minéralisation de l'azote du sol en fonction de la fréquence de restitution des apports organiques et de la gestion des résidus de récolte (tableau 36). Les valeurs de FSyst sont proposées à dires d'experts, mais n'ont pas fait l'objet de validation à partir de données expérimentales. 35 (Communication personnelle Alain Bouthier, Arvalis). L'incidence des apports réguliers de produits organiques exogènes sur la dose d'azote est relativement faible avec cette méthode de calcul. - 62 - Fréquence des apports organiques exogènes et type de produit Jamais 5-10ans 3-4ans 1-2ans Facteurs multiplicateurs en + Résidus de récolte A BC A BC A BC Retour-prairie Cl Enlevés-brûlés 0,80 0,95 0,90 1,00 0,95 1,05 1,00 1,1 En cours d'étude Enfouis 1/2 0,90 1,00 0,95 1,05 1,00 1,10 1,02 1,1 Enfouis 1/1 1,00 1,05 1,00 1,10 1,02 1,20 1,05 1,1 Le Tableau 36 est une table des facteurs systèmes (FSyst) proposée par le COMIFER (extrait du guide COMIFER de calcul de la fertilisation azotée, 2011). Types de produits : A = fumiers (décomposition lente) ; B et C = autres, ainsi que les fumiers de volaille (décomposition rapide). Dans le cas où plusieurs types de produits sont apportés (des A et des BC), alors on privilégie les types A. Cl = couverts intermédiaires Cette méthode présente pour les inventeurs les principaux avantages et inconvénients suivants : - avantage : être reconnue, « estampillée COMIFER » ; elle prend en compte une 5 nouvelle conception de l'effet long terme des PRO (augmentation du pool d'azote organique humifié du sol) ; - inconvénients : (i) pas de prise en compte spécifique des apports très fréquents (tous les ans ou 2 années sur 3), (ii) complexité de l'adaptation de l'utilisation de FSyst dans le formalisme de la Réglette de l'invention, (iii) problème pour identifier la situation 10 correspondant aux valeurs forfaitaires de minéralisation : est-ce que les valeurs de M correspondent à l'absence d'apport organique et à des résidus enfouis une fois sur deux (dans ce cas, FSyst = 0.9 pour des valeurs de M = 20 kgN/ha en sols superficiels et 40kgN/ha en sols profonds) ou est-ce qu'elles correspondent à une autre situation ? « Référentiel Mha breton » 15 L'effet à long terme des PRO est pris en compte dans le poste du bilan historiquement appelé arrière-effets d'apports de PRO (Mha). Dans cette méthode, le terme de minéralisation spécifique lié à des arrières-effets d'apports de PRO n'a pas été remplacé par l'utilisation du facteur système FSyst. Les valeurs proposées dans la table sont issues de références régionales et de l'expertise des spécialistes. Cette table 20 fait référence dans le grand Ouest de la France (tableau 37). Nature des apport) ;«. Fee Ce«rw«t. f.re- tE:c e" Cernr ,>:«itel,s 3es, r..,ektis 25 10 4 15 15 o ans 2 anrées sur s les 3 ans et - 4 35 20 20 Le tableau 37 montre la contribution des apports organiques des années précédant le semis (« Référentiel Mha breton »). -64 - Cette méthode présente pour les inventeurs les principaux avantages et inconvénients suivants : avantages : o intégration aisée de ce poste dans le formalisme de la Réglette 2014; o prise en compte des apports très fréquents qui peuvent être réalisés dans les zones d'élevage du Grand Ouest ; o prise en compte des produits à C/N élevé de type boues et composts qui représentent près d'un quart des apports de PRO sur colza ; inconvénients : o références faisant appel à un formalisme ancien de la méthode du bilan ; o liste des PRO assez restrictive ; o quantités d'azote apportées ne représentant pas toujours la pratique d'apport moyenne (mais règle de trois possible) ; o utilisation de la table assez complexe en l'état dans un outil d'aide à la décision. 15 2.10.2 Méthode retenue Faute de disposer d'une méthode qui fasse consensus et après analyse des deux méthodes précédemment mentionnées, les inventeurs proposent de s'orienter vers une adaptation de la table « Contribution des apports organiques des années précédant le 20 semis (Mha) » utilisée dans le grand Ouest. Toutefois, la prise en compte de ce poste du bilan dans l'outil « Réglette azote colza » pourra être reconsidérée à moyen terme si les recommandations du COMIFER évoluent. Des expérimentations sur les effets à long terme des PRO et des discussions 25 entre spécialistes de la question sont en effet toujours en cours. On peut espérer que les 2 méthodes se complètent et fusionnent. Adaptation du référentiel Mha breton pour une utilisation nationale dans la Réglette azote colza Afin d'élaborer une règle de décision qui convienne pour l'ensemble du territoire 30 national, la table Mha issue du référentiel breton est adaptée de la façon suivante : - Elargissement de la gamme des PRO avec la prise en compte des 20 références couvrant 85% des pratiques d'apports sur colza (cf. chapitre sur Mpro1). L'établissement d'une concordance entre types de produits a, dans certains cas, été nécessaire. Ce travail a été réalisé en collaboration avec Christine Le Souder (Arvalis). 35 Ainsi, le rattachement des produits aux références bretonnes est-il le même dans Fertiweb. -65- - Pondération des effets en fonction des quantités moyennes apportées et des teneurs en azote définies par défaut dans notre table des PRO (application d'une règle de trois) (tableau 38). Compte-tenu de l'enjeu sur la dose d'azote et de l'ambition de conserver un outil facile d'utilisation, l'effet long terme des PRO est calculé avec les 5 valeurs moyennes des pratiques d'apport et des teneurs en azote. L'utilisateur ne sera pas interrogé sur les quantités épandues et les teneurs en azote des produits organiques apportés au cours des dernières années. Ces informations sont en effet difficiles d'accès, contrairement aux pratiques d'apport réalisées en fin d'été ou à l'automne sur le colza. Pour information, le détail de l'historique des pratiques d'apports 10 organiques n'est pas demandé dans la méthode du bilan du Comifer et dans Fertiweb. Dans le cas où plusieurs types de produits sont apportés, les produits à décomposition lente comme les fumiers et les composts sont privilégiés. - Arrondi des valeurs de Mha par pas de 5 et rapprochement des références 15 lorsque cela est possible pour aboutir à une table simplifiée (tableau 3.39). Liste PRO Réglette 2014 Correspondance Teneur N total Qté épandue N total Fourniture d'azote au colza (kgN/ha) en fonction de la fréquence d'apport référentiel "breton" apporté (kg/t ou m de produit brut) (t ou m3/ha de produit brut (k gN/ha) tous les ans 2 années sur 3 tous les 2 ans tous les 3 ans et+ Fumier de bovins Fumier de bovins 5.8 20 116 28 18 14 7 Lisier de bovins Lisier de bovins 4.5 30 135 23 17 11 6 Lisier de bovins dilué Lisier de bovins 1.6 30 48 8 6 4 2 Fumier devolailles Fumier de volailles 25 5 125 13 8 5 3 Lisier devolailles Lisier de volailles 7 15 105 10 5 5 0 Lisier de canards Lisier de volailles 7.7 20 154 15 7 7 0 Fientes de volailles sèches Lisier de volailles 40 5 200 19 10 10 0 Fientes de volailles (60%MS) Lisier de volailles 24 5 120 11 6 6 0 Lisier de porcs mixte Lisier de porcs 3.5 30 105 20 15 10 5 Lisier de porcs à l'engrais Lisier de porcs 5.8 20 116 22 17 11 6 Fumier de porcs Fumier deporcs 8 20 160 20 13 7 0 Fumier de chevaux Fumier de bovins 8 20 160 39 24 19 10 Fumier d'ovins Fumier de bovins 8 20 160 34 21 17 8 Compost fumier bov, et boues Compost fumier bovins 6.3 15 94.5 29 22 15 7 Compost fumier volailles Compost fumier volailles 23 3 69 24 15 11 9 Compost de déchets verts Compost mûr de déchets 10 10 100 29 19 15 8 Boues urbaines liquides Lisier de bovins 3 30 90 15 11 8 4 Boues urbaines pâteuses Lisier de bovins 10 10 100 17 13 8 4 Boues urbaines séchées (90%MS) Fumier de bovins 40 10 400 97 61 48 24 Vinasse de betteraves concentrée Effet nul 20 3 60 0 0 0 0 Ecumes de sucrerie surpressées Lisier de bovins 4 10 40 7 5 3 2 Le tableau 38 est une mise en correspondance des tables Mpro1 « réglette 2014 » et « référentiel Mha breton » avant simplification Fourniture d'azote au colza (kgN/ha)en fonction de la fréquence tous les 2 ans d'apport tous les 3 ans et + Type PRO tous les ans 2 années sur 3 Fumier de bovins 30 20 15 5 Fumier de volailles 15 10 5 5 Fumier de porcs 20 15 5 0 Fumier de cheval 40 25 20 10 Fumier d'ovins 35 20 15 10 Lisier de porcs 20 15 10 5 Lisier de bovins 25 15 10 5 Lisier de bovins dilué 10 5 5 0 Lisier, fientes de volailles (É60%MS) 10 5 5 0 Fientes de volailles séches (80%MS) 20 10 10 0 Boues urbaines liquides et pâteuses (É20%MS) 15 15 10 5 Boues urbaines sèchées (90%MS) 100 60 50 25 Compost de déchets verts 30 20 15 10 Compost de fumier de bovins, boues STEP + déchets verts 30 20 15 5 Compost de fumier de vol ai I les 25 15 10 10 Vinasse de betterave concentrée 0 0 0 0 Ecumes de sucreries surpréssées 5 5 5 0 Le tableau 39 est une table des suppléments de fourniture d'azote liés aux effets des apports réguliers de PRO » proposée pour la Réglette 2014. - 68 - 2.11 Supplément de fourniture d'azote lié à une culture précédente légumineuse (Fleg) 2.11.1 Matériel et méthodes L'effet du précédent pois protéagineux sur la fourniture d'azote au colza suivant a été évalué par le CETIOM et les Chambres d'Agriculture de la Mayenne, de la Moselle, de la Nièvre et de l'Yonne. Douze expérimentations ont été conduites de 2007/08 à 2009/10 pour comparer les courbes de réponse du rendement du colza à la dose d'engrais N minéral en fonction du précédent cultural (pois protéagineux vs. céréale à pailles). Ces essais ont fait l'objet d'une synthèse (rapport de synthèse de l'action 2 du projet CASDAR n°7175, F. Flénet, E. Bonnin, M. Gel oen, 2012) et d'une exploitation complémentaire formalisée dans une note (F. Flénet, octobre 2013). L'étude se concentre sur la comparaison d'un précédent pois protéagineux et d'un précédent céréale à pailles enfouies, en l'absence d'apport de matière organique sur colza, toutes choses égales par ailleurs. La dose d'engrais azoté permettant d'atteindre la marge brute maximale a été étudiée pour des gammes de prix de la graine de colza et d'engrais azoté minéral (20 à 50 E/q graines et 0,50 à 1,25 E/unité N). Ce cahier des charges (pailles enfouies, pas d'apport organique et détermination de la marge brute maximale pour toutes les hypothèses de prix et pour les deux précédents culturaux) a conduit à exclure 4 expérimentations. Les 8 qui ont été utilisées sont recensées dans le tableau 40 (l'expérimentation CA89_2008, pour laquelle cette condition d'identification de la marge brute maximale est respectée à l'exception d'une hypothèse de combinaison de prix (50 E/q graines et 0,5 E / unité), a été « repêchée »).. La part de l'effet N du précédent pois qui provient d'une différence d'absorption d'azote à la sortie de l'hiver a pu être distinguée de celle qui provient de différences de reliquats N à la sortie de l'hiver et/ou de minéralisation de printemps. La première est prise en compte dans le poste Pi (quantité d'azote absorbé par la culture à l'ouverture du bilan). La deuxième correspond à Fleg. Les quantités d'azote absorbé par le colza (Nabs) ont été calculées en appliquant aux observations de biomasse les coefficients de 57 et de 68 kg N / ha / (kg MF / m2), respectivement pour l'entrée (EH) et la sortie de l'hiver (SH) (valeurs moyennes définies dans le §35). La quantité d'azote à prendre en compte dans le bilan N (paramètre Nh) a été calculée comme suit : si Nabs_SH > Nabs_EH, ou si absence de mesure à EH, alors Nh = Nabs_SH, Sinon, Nh = (Nabs_EH + Nabs_SH) / 2. - 69 - Année Code de l'essai Organisme Lieu Type de sol 2007/08 CA58_2008 Chambre d'Agriculture Nièvre Limon argileux CA89 2008 Chambre d'Agriculture Yonne Argilo-calcaire 2008/09 CETIOMchât 2009 CETIOM Indre Moselle Somme Yonne Sablo-argilo-limoneux CA57 2009 Chambre d'Agriculture Limon - CETIOM Limon moyen CETIOMmons 2009 Chambre d'Agriculture Limon à silex _ CA89 2009 2009/10 CETIOMchât 2010 CETIOM Indre Limon argileux CETIOMmons 2010 CETIOM Somme Limon sableux Le tableau 40 représente les principales caractéristiques des essais retenus pour l'estimation de Fleg. 2.11.2 Résultats et discussions La synthèse des résultats montre qu'en précédent pois, la disponibilité en N est plus élevée qu'en précédent céréale à pailles : rendement du colza plus élevé sans apport d'engrais N et rendement maximal atteint pour une dose d'engrais N plus faible (Figure 27). En revanche, il n'y a pas de différence de rendement maximal du colza entre les précédents culturaux. En moyenne la marge brute maximale est obtenue pour 41 kg N / ha de moins en précédent pois et le rendement en graines obtenu à cette dose est en moyenne équivalent à celui du précédent céréale à pailles. Les hypothèses de prix de la graine et de l'engrais ont peu d'effet sur cet écart de dose N entre les précédents. Le tableau 40 présente les résultats pour 40 E / q et 0,75 E / unité. Pour ce scénario, la différence de dose d'engrais N n'est pas significative au seuil de 5%, en raison de la variabilité entre les expérimentations qui s'explique notamment par le dispositif expérimental. La différence, entre les précédents, de dose d'engrais N pour atteindre la marge brute maximale provient pour partie d'une disponibilité en azote plus élevée à l'« automne » en précédent pois, qui se traduit par une biomasse de colza plus élevée à l'entrée et à la sortie de l'hiver (tableau 41). En moyenne, cette différence de biomasse correspond à 14 kg N / ha. La différence entre les deux précédents n'est pas significative au seuil de 5%, en raison d'une variabilité importante, notamment parce que l'absorption d'azote a été très faible en précédent pois dans l'expérimentation CETIOM mons_2010. Une mauvaise levée, causée par un dessèchement du lit de semence par le pois, a en effet été observée dans cette expérimentation. - 70 - Code essai Précédent Dose engrais Biomasse Biomasse Nh (kg N / ha) entrée hiver sortie hiver (kg N / ha (kg MF/m2) (kg MF/m2 CA89_2008 Pois 141 nc* 1.6 108.8 CA89 2008 Céréale 120 nc 1.5 102 CA57_ 2008 Pois 94 nc 1.25 85 CA57 2008 Céréale 186 nc 1.025 69.7 CETIOMmons_2009 Pois 67 0.4 0.4 27.2 CETIOMmons 2009 Céréale 296 0.4 0.4 27.2 CETIOMchât_2009 Pois 124 2.56 2.19 148.9 CETIOMchât 2009 Céréale 107 2.155 0.98 94.7 CA57_2009 Pois 213 1.45 1.1 78.7 CA57 2009 Céréale 152 1.18 0.75 59.1 CA89_2009 Pois 158 2 2 136 CA89 2009 Céréale 191 1.35 1.025 73.3 CETIOMmons_2010 Pois 217 0.25 0.25 17 CETIOMmons 2010 Céréale 255 1.2 0.8 61.4 CETIOMchât_2010 Pois 154 1.69 1.25 90.7 CETIOMchât 2010 Céréale 192 1.54 1.36 92.5 Moyenne Pois 146 1.26 86.5 Céréale 187 0.98 72.5 Test de Student P = 0.115 P = 0.138 Le tableau 41 represente la dose d'engrais azoté pour atteindre la marge brute maximale (exemple du prix de la graine de 40 E / q et du prix de l'engrais N de 0,75 E / unité), biomasses fraîches aériennes à l'entrée et à la sortie de l'hiver et quantité d'azote déjà absorbé par le colza à l'ouverture du bilan (Nh) pour les 8 expérimentations retenues pour comparer l'effet du précédent cultural (* nc = non connu). La synthèse des résultats permet de conclure qu'en moyenne, le pois protéagineux permet de réduire la fertilisation azotée du colza d'environ 40 unités, dont 15 sont déjà prises en compte par l'azote déjà absorbé à l'ouverture du bilan. En conséquence, 25 unités correspondent au supplément de fourniture en azote lié au précédent pois par rapport au précédent céréales à pailles (Fleg). Elles doivent être retranchées du calcul de la dose d'azote en écriture additive La valeur de Fleg utilisée en écriture additive a été convertie par calcul en valeur de Fleg' pour l'écriture CAU grâce au formalisme suivant : Fleg' = (effet « légumineuse » moyen sur la dose N x CAU) - écart de Pi d'où Fleg' = (40 x 0.8) - 15 = 17 kgN/ha. Ces résultats proviennent de 8 expérimentations seulement. Il est nécessaire 20 d'acquérir de nouvelles références expérimentales pour les consolider et, si possible, pour les préciser en prenant en compte de l'état de croissance du pois protéagineux. - 71 - De la même manière, les expérimentations en cours et à venir seront nécessaires pour prendre en compte d'autres précédents protéagineux, notamment la féverole. 2.12 Supplément de fourniture d'azote lié à l'association (Fass) Extrait de la synthèse technique "Colza d'hiver associé à un couvert automnal" (Cadoux, 2014). Les colzas seuls contiennent en moyenne plus d'azote en entrée hiver que les colzas associés (+5.5 kgN/ha ; effectif : 45 en colza seul et 456 en colza associé). La 10 différence est non significative pour les couverts à base de légumineuses pures. Qu'il soit associé ou non, le colza présente des quantités d'azote similaires à biomasse équivalente à l'entrée et à la sortie de l'hiver, qu'il soit associé ou pas. Les paramètres coeffEH et coeffSH du colza seul peuvent donc être appliqués sans modification aux colzas associés. 15 Au stade F2, les biomasses et les quantités d'azote absorbé sont proches entre les colzas seuls fertilisés à la dose X selon la réglette 1998 et associés avec des légumineuses ou des mélanges de légumineuses et non légumineuses fertilisés à la dose X - 30 kgN/ha. Dans les essais non fertilisés les colzas associés à des légumineuses contiennent en moyenne 8 kg N de plus que les colzas seuls. 20 À G4, dans les modalités fertilisées, les colzas seuls fertilisés à la dose X contiennent en moyenne 15 kgN/ha de plus que les colzas associés avec des légumineuses fertilisés à la dose X - 30 kgN/ha. Dans les essais non fertilisés les colzas associés à des légumineuses contiennent en moyenne 8 kg N de plus que les colzas seuls. 25 A la récolte, malgré une réduction de fertilisation de 30kgN/ha d'azote pour les colzas associés par rapport aux colzas seuls, le rendement est équivalent, dans la plupart des cas. Il est toutefois légitime de s'interroger sur la pertinence de la dose X utilisée. En effet, si celle-ci était supérieure à la dose optimale a posteriori d'une trentaine d'unités ou plus, il ne serait pas étonnant que le rendement du colza associé 30 ne soit pas affecté à la dose X-30 kgN/ha par rapport au colza seul à la dose X. Pour écarter au moins en partie ce doute, les écarts de rendement ont été recalculés sur les situations ou l'indice de nutrition azotée (INN) au stade F2 étaient inférieurs à 1 sur les colzas seuls (situation probablement non sur-fertilisées). Il apparaît alors que les écarts de rendements entre les colzas seuls à la dose X et les colzas associés à la dose X-30 35 kgN/ha restent non significatifs. Ceci conforte l'hypothèse d'équivalence de rendement pour une fertilisation réduite de 30 kgN/ha sur les colzas associés. - 72 - En résumé, les couverts associés semblent fournir au colza, un supplément d'azote moyen mais très variable d'environ 10 kg/ha ainsi que d'autres bénéfices qui auraient permis dans nos essais de compenser l'écart d'environ 20kg/ha d'azote disponible (X30kg de fertilisation + 10kg de bénéfices des couverts associés). Il est alors possible de recommander, un abattement forfaitaire de fertilisation par rapport à la dose X de 30 kgN/ha en visant un rendement équivalent. Le paramètre Fass prend donc la valeur de 30 kgN/ha. Il peut être cumulé à Fleg en cas de précédent pois protéagineux. Dans tous les cas, l'abattement ne sera appliqué que si le couvert associé est à 10 base de légumineuses pures et s'il a levé. L'abattement ne sera dépendant ni de la croissance ni du gel du couvert associé. 2.13 Coefficient apparent d'utilisation de l'azote issu de l'engrais (CAU) 2.13.1 Matériel et méthodes 15 Le coefficient apparent d'utilisation de l'azote de l'engrais a été ré-estimé à partir de la base de données des essais « doses d'azote » 2008 à 2012 déjà mentionnée précédemment. La grande majorité des essais a été fertilisée avec de l'ammonitrate 33.5%. Le calcul du CAU a été réalisé de la façon suivante : calcul de la pente de la 20 régression linéaire de la relation exprimant la quantité d'azote absorbé dans les plantes entières à G4 en fonction de la dose d'azote apporté à G4. Ce calcul a été réalisé en ne conservant que les doses d'azote permettant d'atteindre un rendement inférieur au plateau de rendement, plus la dose la plus faible permettant d'atteindre ce plateau de rendement. Dans les essais où le plateau de rendement n'est jamais atteint, le calcul a 25 été réalisé sur toutes les doses. Dans cette étude, le plafond de rendement est défini comme le rendement observé le plus faible dans le groupe statistique homogène le plus élevé au sens du test de Student, Neuman et Keuls au seuil de 10% (logiciel SAS). Seuls les essais présentant une régression linéaire avec un coefficient de détermination supérieur à 0.8 ont été retenus. Cela représente 31 essais. Ont 30 également été écartés les essais pour lesquels seulement deux modalités de dose d'azote satisfaisaient le critère lié au test de Student, Neuman et Keuls (6 essais). 25 essais ont donc été exploités. Les valeurs des déciles ont été calculées avec Excel 2010. Quelques caractéristiques de la base de données sont présentées sur les Figures 35 28A à 28F. Il convient de noter une nouvelle fois la forte représentation des années 2008 et 2009 et des régions Centre et Poitou-Charentes. - 73 - 2.13.2 Résultats Les valeurs CAU vont de 0,34 à 1,48 (Figure 29). Sept essais présentent des valeurs de CAU supérieures à 1, ce qui traduit la difficulté d'estimer ce coefficient. La valeur moyenne est de 0,85 et la valeur médiane est de 0,87. Ces valeurs ne sont que légèrement supérieures à la valeur retenue pour la Réglette azote 1998 (0,8). En conséquence, les inventeurs proposent de conserver la valeur retenue pour la réglette azote 1998, à savoir 0,8. 2.14 Plafonnement de la dose totale d'azote Sur la base des valeurs de paramètres retenues dans les chapitres précédents, le calcul de la dose d'azote à appliquer peut conduire à des valeurs très élevées, surtout en méthode CAU (tableau 42). Il a donc semblé judicieux aux inventeurs d'introduire une notion de dose d'azote plafond. 2.14.1 Matériel et méthodes La base de données des essais « doses d'azote » 2008 à 2012, déjà mentionnée précédemment, a été mobilisée. Deux catégories d'essais ont été retenues : - 43 essais pour lesquels un plafond de rendement a été atteint selon l'ajustement quadratique + plateau. Le plafond de dose d'azote a été estimé à partir du calcul des centiles (Excel 2010) de la variable « dose N » permettant d'atteindre le rendement maximal d'après l'ajustement quadratique + plateau. - 25 essais pour lesquels un plafond de rendement n'a pas été atteint selon l'ajustement quadratique + plateau. 2.14.2 Résultats Parmi les 43 essais présentant un plateau de rendement, seul un essai présente une 25 dose minimale permettant d'atteindre le rendement maximal plus élevée que 230 kgN/ha (tableau 43 et figure 30A). Hypothèse de Nombre d'essais % d'essais ayant dose N plafond ayant atteint le atteint le rendement maxi rendement maxi 150 24 55.8 160 29 67.4 30 170 34 79.1 180 39 90.7 190 40 93.0 200 41 95.3 210 41 95.3 220 41 95.3 230 42 97.7 240 42 97.7 35 250 42 97.7 260 43 100.0 Le tableau 43 représente les nombre et pourcentage d'essais pour lesquels le rendement maximal est atteint pour différentes doses d'azote selon un ajustement quadratique + plateau (base de données de 43 essais) - 74 - Parmi les essais ne présentant pas de plateau de rendement, très peu comprennent une dose d'azote appliquée maximale nettement supérieure à 250 kgN/ha (figure 30B droite). Trois groupes d'essais peuvent être distingués. L'examen visuel des courbes noires (20 essais) combiné à l'analyse des doses d'azote théoriques permettant 5 d'atteindre les plateaux de rendement (doses supérieures à la dose d'azote maximale appliquée sur l'essai) permettent de penser que le plateau de rendement aurait été atteint pour une dose inférieure à 250 kgN/ha. Pour les 4 essais correspondant aux courbes roses, la dose d'azote qui permet d'atteindre le rendement maximal est clairement supérieure à 250 kgN/ha. Enfin l'examen visuel de la courbe bleue, laisse 10 également penser que le plateau de rendement aurait été atteint pour une dose d'azote supérieure à 250 kgN/ha. De plus, il est rappelé que la base de données est biaisée en raison d'une surreprésentation des situations avec des colzas moyens à gros et en sols profonds. 15 C'est pour toutes ces raisons que les inventeurs se sont orientés vers un plafond de dose relativement élevé, à savoir 250 kgN/ha. Essai objectif MV Dose prévisionnelle Réglette 1998 (kgN/ha) Dose prévisionnelle Réglette 2014 (kgN/ha) Dose optimale Type de sol rendement (q/ha) réelle (courbe de réponse) (kgN/ha) ouverture Tables (plafond Calcul équation Déciles 5-5 Déciles 7-3 Déciles 9-1 bilan objectif rendement) réglette avec objectif (kg/m2) de rdt réel C09MFN54028 39.86 0.620 190 (35 q/ha) 232 235 270 291 69 sup C09MFN54029 49.22 0.520 170 (40 q/ha) 267 295 314 318 211 prof CA57_2009 50.34 0.965 170 (40 q/ha et +) 228 268 282 290 126 prof Le tableau 42 montre un exemple de 3 essais en méthode CAU présentant des doses d'azote conseillées élevées surtout avec les paramètres de la réglette 2014, malgré le plafonnement du besoin à 330 kgN/ha. - 76 - Exemple 2: Exemple de mise en oeuvre de la méthode selon l'invention (Réglette 2014) 1- Département / Choix de la Réglette L'utilisateur déclare son département pour définir quelle « Réglette » sera utilisée. Il choisit son département dans un menu déroulant (n° dpt - nom département, par exemple) L'essentiel est d'avoir un système souple de passage du département à la réglette en cas de modification. Le choix du département, contrairement à une carte cliquable, 10 limite également le questionnement des utilisateurs sur les différences entre les réglettes régionales. Pour information, les limites géographiques des « régions réglette », dans l'état actuel de nos travaux, sont représentées sur la carte ci-après. La Figure 31 représente une carte de France montrant les types de réglettes 15 utilisées. 2- Type de sol L'utilisateur choisit son type de sol parmi les propositions suivantes : superficiel, profond (pour toutes les réglettes), craies de Champagne (uniquement pour la réglette 20 Champagne-Ardenne Aisne). Cela revient à proposer une réponse supplémentaire pour la réglette Champagne-Ardenne Aisne et à adapter les fenêtres pop-up d'aide. Sols superficiels : Ils se caractérisent par une faible minéralisation au printemps et une réserve en eau faible, susceptible de limiter le rendement. Exemples : (se 25 référer au tableau ci-dessous) Sols profonds : Ils se caractérisent par une minéralisation moyenne à élevée au printemps et une bonne réserve en eau ne limitant pas (ou peu) le rendement. Exemples : (se référer au tableau ci-dessous) Sols de craies de Champagne : Ils se caractérisent par une minéralisation moyenne 30 au printemps et une bonne réserve en eau ne limitant pas (ou peu) le rendement. Exemples : craies profondes, craies moyennement profondes, craies à poches, rendzines grises. Les exemples de sols superficiels et profonds cités dans la fenêtre d'aide varient selon le département choisi. Les correspondances géographiques et les dénominations 35 des sols sont présentées dans le tableau 43 ci-dessous. - 77 - Départements Exemples de type de sol Exemples de type de sol associés associés aux « sols associés aux « sols superficiels » profonds » 59, 62, 76, 60, 80 Sables, Biefs, Argiles à silex, Limons, Limons argileux, Argilo-calcaires superficiels et Argilo-calcaires profonds moyens, Cranettes 08, 10, 51, 02, 77 Sols sableux, Argiles à silex, Argiles, Limons profonds, Argilo-calcaires superficiels et Limons moyens, Limons moyens, Cranettes, argileux, Sols argilo- Graveluches, Craies limoneux, Argilo-calcaires superficielles, Gréves profonds, Rendzines colorées, Terres de vallées, Terres humifères 22, 29, 35, 56, 14, 50, Argilo-calcaires, Sols Limons 61, 44, 49, 53, 72, superficiels sur schistes dégradés 18, 28, 36, 37, 41, 45, Sables, Argilo-calcaires Limons, Limons argileux, 78, 91, 92, 93, 94, 95, superficiels et moyens, Limons sableux sains, Argilo- 27 Limons caillouteux calcaires profonds superficiels, Limons humides 67, 68, 52, 54, 55, 57, Argilo-calcaires superficiels et Argiles, Argilo-calcaires 88 moyens, Sables profonds, Limons, limons argileux 19, 23, 87, 85, 16, 17, Sables, Groles légères Groles profondes, 79, 86 Champagne, Terres rouges, Bornais, Varennes, Au bues, Doucins 03, 15, 43, 63 Argilo-calcaires plus ou Terres noires, Argilo-calcaires moins superficiels, Sols profonds, Sols limons granitiques, Alluvions, Sols sableux hydromorphes sableux, Argilo-sableux profonds et drainés 21, 58, 71, 89, 25, 39, Sables, Argilo-calcaires Limons blancs, Limons 70, 90, superficiels et moyens, argileux, Limons sableux, Argiles à silex Argilo-calcaires profonds, Alluvions 01, 07, 26, 38, 42, 69, Sols sableux, Argilo-calcaires Limons sablo limoneux sains, 73, 74 superficiels, Graviers plus ou Limons argileux profonds et moins profonds sains, Limons humides, Argiles ou Limons argileux profonds, Argiles humide, Argilo-calcaires profonds, Alluvions organiques sains, Marais humides 24, 33, 40, 47, 64, 11, Coteaux et plateaux argilo- Alluvions fertiles de vallées 30, 34, 48, 66, 09, 12, calcaires superficiels ou peu (Garonne, Hers, Girou, Tarn, 31, 32, 46, 65, 81, 82, profonds, Sols caillouteux vallées du Gers, Orb, Hérault, 04, 05, 06, 13, 83, 84, (Grausses, Costières...), Rhône, Durance...), Calcaires superficiels et Boulbènes profondes, Groies, Boulbènes Terreforts profonds des bas superficielles ou de coteaux du Sud-Ouest, caillouteuses, Chalosse, Limons du littoral Razès, Sables des Landes, méditerranéen Ségala - 78 - Chaque interlocuteur régional s'est exprimé sur les exemples de type de sol. Sur la région Poitou-Charentes, Vendée, Limousin, les appellations régionales sont préférées aux dénominations communes des sols. Pour la zone Sud, c'est une combinaison des grands types de sols, des noms locaux et des indications géographiques qui est retenue. Ailleurs, ce sont principalement les dénominations par grand type de sol qui sont privilégiées dans un souci d'harmonie avec les référentiels agronomiques des arrêtés Directive Nitrate. Dans tous les cas, ils demeurent des exemples pour aider l'utilisateur à se positionner mais ils ne constituent pas une liste exhaustive des sols qu'il est possible de rencontrer dans le département. 3-1- Historique des apports organiques sur la parcelle L'utilisateur indique s'il apporte régulièrement des produits résiduaires organiques (PRO) sur la parcelle. La réponse « pas d'apport organique » est affichée par défaut. Si oui, il doit choisir un PRO dans une liste (tableau 44 ci-dessous) et définir une 15 fréquence d'apport (tous les ans, 2 années sur 3, tous les 2 ans, tous les 3 ans et plus). L'utilisateur ne peut choisir qu'une seule combinaison PRO x fréquence d'apport. Table des arrières effets des PRO utilisée dans la nouvelle réglette azote colza. Produits Résiduaires Organiques tous les 2 années tous les tous les ans sur 3 2 ans 3 ans et plus Fumier de bovins 30 20 15 5 Fumier de volailles 15 10 5 5 Fumier de porcs 20 15 5 o Fumier de cheval 40 25 20 10 Fumier d'ovins 35 20 15 10 Lisier de porcs 20 15 10 5 Lisier de bovins 25 15 10 5 Lisier de bovins dilué 10 5 5 o Lisier, fientes de volailles (560%MS) 10 5 5 0 Fientes de volailles séches (80%MS) 20 10 10 o Boues urbaines liquides et pâteuses (520%MS) 15 15 10 5 Boues urbaines sèchées (90%MS) 100 60 50 25 Compost de déchets verts 30 20 15 10 Compost de fumier de bovins boues STEP + déchets verts 30 20 15 5 - 79 - Produits Résiduaires Organiques tous les 2 années tous les tous les ans sur 3 2 ans 3 ans et plus Compost de fumier de volailles 25 15 10 10 Vinasse de betterave concentrée 0 0 0 0 Ecumes de sucreries surpréssées 5 5 5 0 Tableau 44 Pour la méthode additive : Historique des apports organiques sur la parcelle : On considère qu'il y a un « arrière effet » des apports réguliers de produits organiques si la parcelle reçoit régulièrement, depuis une longue période, du fumier, du lisier, des fientes, des boues, des composts, etc. Dans le cas où plusieurs types de produits sont apportés, alors on privilégie les produits à décomposition lente comme les fumiers et les composts. La mesure du reliquat sortie hiver est conseillée dans les parcelles recevant régulièrement de la matière organique si le colza est très vert ou petit. Pour la méthode CAU : Historique des apports organiques sur la parcelle : On considère qu'il y a un « arrière effet » des apports réguliers de produits organiques si la parcelle reçoit régulièrement, depuis une longue période, du fumier, du lisier, des fientes, des boues, des composts, etc. Dans le cas où plusieurs types de produits sont apportés, alors on privilégie les produits à décomposition lente comme les fumiers et les composts. 3-2- Apport de produit organique à l'automne sur colza L'utilisateur indique s'il a apporté un produit résiduaire organique (PRO) à l'automne (avant le semis du colza ou en végétation). La réponse « pas d'apport organique » est affichée par défaut. Le cas échéant, il doit choisir un PRO dans une liste (table des PRO ci-dessous). On 25 lui indique qu'il peut modifier les valeurs par défaut proposées pour la teneur en azote et la quantité épandue. Un contrôle de saisie est prévu sur la teneur en azote et la quantité de produit épandue : l'utilisateur est alerté si valeur 100. L'utilisateur peut saisir jusqu'à 2 PRO appliqués à l'automne. Si deux saisies sont réalisées, les effets se cumulent. 30 Tableau 45 des PRO proposée par défaut dans l'outil pour calculer l'effet d'un apport - 80 - organique à l'automne. NB : Le coefficient d'équivalence « Keq sur la période du bilan colza » est caché pour l'utilisateur Produit Résiduaire Organique Teneur en azote total (Kg N i t ou m3) Quantité de . Coefficient épproduit andue par unite . . hectare (tou m3 / ha) d'equivalence engrais (Keq) sur la période du bilan colza Fumiers Fumier de bovins 5.8 20 t / ha 0.1 Fumier de volailles 25 5 t / ha 0.2 Fumier de porcs 8 20 t / ha 0.1 Fumier de cheval 8 20 t / ha 0.1 Fumier d'ovins 7 20 t / ha 0.1 Lisiers Lisier de porcs à l'engrais 5.8 20 m3 / ha 0.05 Lisier de porcs mixtes 3.5 30 m3 / ha 0.05 Lisier de bovins 4.5 30 m3 / ha 0.15 Lisier de bovins dilué 1.6 30 m3 / ha 0.15 Lisier de canards 7.7 20 m3 / ha 0.05 Fientes Fientes de volailles séches (80%MS) 40 5 t /ha 0.05 Fientes de volailles (60%MS) 24 5 t /ha 0.05 Boues urbaines m3 / ha Boues urbaines liquides 3 30 0.15 Boues urbaines pâteuses (20%MS) 10 10 t! ha 0.1 Boues urbaines sèchées (90%MS) 40 10 t! ha 0.1 Composts t! ha Compost de déchets verts 10 10 0.05 Compost de fumier de bovins 6.3 15 t / ha 0.12 Compost de fumier de volailles 23 3 t / ha 0.12 Autres 3 0.15 Vinasse de betterave concentrée 20 t! ha Ecumes de sucreries surpréssées 4 10 t / ha 0 Autre produit organique "champ libre" "champ libre" t / ha 0.1 Tableau 45 Pour la méthode additive : Apport de produit organique à l'automne sur colza : En cas d'apport de produit organique sur le colza en fin d'été ou à l'automne, que ce soit avant le semis ou en végétation, indiquez le type de produit épandu. Compte-tenu de la forte variabilité des compositions moyennes des produits et des pratiques d'épandage, il est préférable d'utiliser des valeurs mesurées. Si vous en disposez, modifiez les valeurs proposées par défaut. La mesure du reliquat sortie hiver est conseillée dans les parcelles recevant de la - 81 - matière organique si le colza est très vert ou petit. Pour la méthode CAU : Apport de produit organique à l'automne sur colza : En cas d'apport de produit organique sur le colza en fin d'été ou à l'automne, que ce soit avant le semis ou en végétation, indiquez le type de produit épandu. Compte-tenu de la forte variabilité des compositions moyennes des produits et des pratiques d'épandage, il est préférable d'utiliser des valeurs mesurées. Si vous en disposez, modifiez les valeurs proposées par défaut. 4- Précédent pois protéagineux L'utilisateur indique si le colza est implanté après un pois protéagineux en répondant par oui ou par non (ou case à cocher). La réponse « non » est affichée par défaut car le précédent céréale à paille est largement majoritaire. Précédent pois protéagineux : Indiquer si le colza est implanté après un pois 15 protéagineux. 5- Colza associé à l'automne avec un couvert de légumineuses gélif L'utilisateur indique si le colza est implanté à l'automne avec un couvert de légumineuses gélif en répondant par oui ou par non (ou case à cocher). La réponse « 20 non » est affichée par défaut car les associations sont pratiquées sur moins de 3% des surfaces à ce jour. Colza associé : Indiquer si le colza est implanté avec un couvert de légumineuses gélif (vesces, gesse, lentille, féverole, fenugrec...). Si le couvert n'a pas levé, l'effet de l'association est négligé. Dans ce cas, ne cocher pas la case « colza associé ». 25 6- Objectif de rendement L'utilisateur indique son objectif de rendement exprimé en q /ha en saisissant une valeur entière (attention au respect du format requis et aux erreurs de saisie). Il doit pourvoir saisir sa propre valeur sans plancher ni plafond de rendement. 30 Objectif de rendement : Ne choisissez pas le rendement record de la parcelle. Pour calculer l'objectif de rendement, prenez la moyenne des rendements colza réalisés sur l'exploitation pour des conditions comparables de sol au cours des 5 dernières années, en enlevant l'année la plus haute et l'année la plus basse. Si les références disponibles sont insuffisantes, reportez-vous au référentiel de mise en 35 oeuvre de l'équilibre de la fertilisation azotée publié dans chaque région. - 82 - 7- Poids frais du colza : entrée hiver et sortie hiver L'utilisateur indique le poids sortie hiver ou les poids entrée hiver et sortie hiver exprimés en kg/m2 en saisissant une ou des valeurs avec décimale. Il doit pouvoir saisir sa propre valeur sans plancher ni plafond de biomasse. Attention, le risque d'erreur de saisie est important (faute de frappe, séparateur décimal et surtout erreur d'unité g/m2, kg/m2 voire t/ha). Prévoir des contrôles de cohérence. Méthodes d'estimation du poids frais du colza 1- méthode par pesée (recommandée) + choisir 2 à 4 placettes de 1 m2 chacune, représentatives de la parcelle + délimiter chaque placette, puis prélever les plantes, lorsque la végétation est ressuyée (en absence de rosée ou de pluie) + couper les plantes au niveau du collet, au ras du sol + prélever en entrée et en sortie d'hiver. Dans les régions froides, faire la pesée 15 entrée d'hiver avant la destruction des feuilles par le gel. A la sortie d'hiver, prélever juste avant la date prévue du 1er apport d'azote au printemps : en décembre/janvier dans le Sud, en janvier/février dans le Nord. - peser les plantes fraîchement coupées sur chaque placette sans séchage. + calculer ensuite la moyenne des pesées réalisées sur les placettes et reporter ce 20 poids moyen dans le champ "Poids frais du colza en kg/m2" en arrondissant à la valeur la plus proche. Si la parcelle comprend plusieurs zones avec des densités ou des niveaux de croissance très différents, il est intéressant de réaliser la même opération sur chacune de ces zones (2 à 4 placettes par zone). Les outils d'imagerie satellitaire facilitent 25 l'identification de ces différentes zones. Vous calculerez donc des doses d'azote différentes sur chaque zone. Si possible, vous appliquerez ces doses différentes sur chacune d'elles. Sinon, vous devrez calculer la dose à apporter uniformément sur la parcelle. 30 2- Méthode visuelle L'observation des parcelles à la sortie d'hiver et les photos suivantes permettent d'estimer le poids frais du colza. Toutefois, cette méthode est moins précise que la méthode par pesée, notamment au delà d'1 kg/m2. 35 Les aides sont représentées aux figures 32A à 32E. - 83 - 3- Méthode Image IT ImagelT est une application pour smartphones proposée par YARA avec la collaboration de la demanderesse. Elle évalue la biomasse et la quantité d'azote absorbée par vos cultures de colza en entrée et en sortie d'hiver grâce au traitement 5 informatique de photos numériques. Compte-tenu de l'état actuel de son développement, la méthode requiert le respect de certaines conditions d'utilisation et n'est pas adaptée pour les colzas dépassant 1.2 kg/m2 (se référer aux consignes d'utilisation de l'outil). Pour calculer la dose d'azote à apporter reporter le poids frais estimé par image IT à l'entrée de l'hiver et/ou à la sortie de l'hiver dans les champs 10 concernés. Attention, si le poids frais est exprimé en t/ha dans le rapport ImagelT, il est nécessaire de le convertir en kg/m2 en divisant par 10 la valeur. Exemple : Le rapport image IT indique une biomasse sortie hiver d'environ 11 t/ha. Saisissez dans le champ « poids frais du colza sortie hiver », la valeur de 1.1 kg/m2. 15 8- Lancement du calcul de la dose d'azote et gestion des situations non réalistes Après avoir renseigné l'ensemble des champs précédemment cités, l'utilisateur commande le calcul de la dose d'azote en cliquant sur un bouton d'action « Calculer la dose d'azote ». Sur l'interface d'entrée, avant que le calcul soit lancé, un système expert s'assure de 20 la cohérence d'ensemble. Il est souhaité, pour cette nouvelle version de la réglette azote colza, que l'outil ne bloque pas le calcul mais alerte l'utilisateur sur la cohérence de sa saisie (warning). A lui ensuite de modifier ses paramètres ou de les confirmer et de lancer le calcul. Le système expert s'applique à plusieurs niveaux sur l'interface d'entrée :Biomasses 25 EH et SH : Si l'utilisateur saisit une biomasse 5 kg, l'outil affiche une alerte: l'utilisateur est invité à vérifier la valeur saisie, il la modifie ou la confirme. Objectif de rendement : Si l'utilisateur saisit un rendement 50 q/ha, l'outil affiche une alerte l'utilisateur est invité à vérifier la valeur saisie, il la modifie ou la confirme. 30 Biomasse colza SH x Objectif de rendement x Type de sol : Si l'utilisateur saisit un objectif de rendement non réaliste au regard de la biomasse et du type de sol, l'outil affiche une alerte l'utilisateur est invité à ré-évaluer l'objectif de rendement de sa parcelle. Les plafonds de rendement proposés en fonction du poids vert SH et du type de sol sont présentés dans le tableau 46 suivant : 35 -84 - Sol sup Sol prof Craie PV É 200g É 35 q/ha É 40 q/ha É 40 q/ha 200 g 330, alors b x y = 330 Plafond de X X = MV EH x coeff EH, alors Pi = MV SH x coeff SH, sinon Pi=Nabs SH + (1/2 (Nabs EH - Nabs SH) /1.35) avec Nabs = MV x coeff 25 Si MV EH pas disponible, alors Pi = MV SH x coeff SH Si MV SH pas disponible, le calcul de Pi est impossible Calcul de Mpro1 : Mpro1 = teneur N PRO (kg/t MB) x qté PRO épandue (t ou m3 MB/ha) x keq 30 Si 2 PRO sont appliqués à l'automne, les valeurs de Mpro1 s'additionnent. X: Dose d'azote à apporter (kgN/ha) b : besoin unitaire en kgN absorbé à la fermeture du bilan / q de graines produites (/ha) 35 y: objectif de rendement (q/ha aux normes) Pi : quantité d'azote absorbé par la culture à l'ouverture du bilan (kgN/ha) Np (= PO) fourniture d'azote par le sol pendant la période d'ouverture du bilan - 87 - (kgN/ha) Rf : quantité d'azote minéral dans le sol à la fermeture du bilan (kgN/ha) Ri : quantité d'azote minéral dans le sol à l'ouverture du bilan (kgN/ha) M : minéralisation nette de l'humus et des résidus de récolte du précédent (kgN/ha) Mpro1 : minéralisation nette de l'azote organique des PRO apporté l'année de la culture avant l'ouverture du bilan (kgN/ha) Mha : fourniture d'azote liée à l'historique d'apport de PRO Fleg : supplément de fourniture d'azote lié à un précédent pois protéagineux (kgN/ha) Fass : supplément de fourniture d'azote lié à des cultures compagnes (kgN/ha) CAU : coefficient apparent d'utilisation MV : poids de matière verte aérienne (kg/m2) EH : entrée hiver SH : sortie hiver (ouverture du bilan) Coeff : coefficient de conversion de la biomasse verte aérienne (kg/m2) en quantité d'azote absorbé (kgN/ha) Valeurs des paramètres : b = 7.0 kgN absorbé plantes entières / q de graines produites Rf : - sols superficiels : 15 - sols profonds : 30 sols de craie de Champagne : 30 Coeff EH: 50 Coeff SH :65 Ri : table 30 - Grand Ouest o Sols superficiels : 30 o Sols profonds : 40 - Reste France o Sol superficiels : 20 35 o Sols profonds : 30 o Sols de craie de Champagne : 40 25 - 88 - M : table - Grand Ouest o Sols superficiels : 30 o Sols profonds : 50 - Reste France o Sol superficiels : 20 o Sols profonds : 40 o Sols de craie de Champagne : 20 Np : table Nord-Est : o Sols superficiels : 30 o Sols profonds : 50 - Auvergne et Rhône-Alpes : o Sols superficiel : 40 o Sols profonds : 60 Mha : table du chapitre 3.1 Fleg: Si précédent pois protéagineux, alors Fleg= 25 Sinon Fleg =0 Fass: Si colza associé à un couvert de légumineuses gélif ; alors Fass= 30 Sinon Fass =0 CAU : 0,8. REVENDICATIONS1. Méthode de fertilisation d'une culture de colza d'hiver en France métropolitaine à l'aide d'une quantité optimale d'engrais azoté, au printemps, dans le but de stimuler la croissance de ladite culture de colza, ladite méthode comprenant les étapes suivantes : a) une étape de mesure d'au moins la masse de matière fraiche aérienne à la sortie de l'hiver MVSH d'au moins deux échantillons de ladite culture, éventuellement la masse de matière verte aérienne à la sortie de l'hiver MVSH et la masse de matière verte aérienne à l'entrée de l'hiver MVEH, 10 b) une étape de collecte des paramètres sietants : la localisation géographie L de ladite culture, - la nature N du sol sur lequel se trouve ladite culture, - les pratiques culturales suivantes ayant un effet sur la fourniture d'azote à ladite culture: 15 - la présence sur la parcelle de pois protéagineux pendant la campagne précédent ladite culture, - la co-culture de ladite culture de colza avec des plants de légumineuses gélifs, - les habitudes d'apport de produits résiduaires organiques sur le sol 20 où se trouve ladite culture pendant les années précédent ladite culture, - l'apport de produits organiques pendant l'été précédent la culture et sur ladite culture pendant l'été et l'automne, et 25 l'objectif de rendement y pour ladite culture, c) une étape de calcul de la quantité d'engrais azoté à ajouter X à ladite culture selon la méthode des bilans prévisionnels simplifiée et adaptée au colza, de sorte que : (i) si la culture de colza est localisée dans une région de France métropolitaine qui n'est pas la région Nord-Est ou la région 30 Auvergne-Rhône-Alpes, la quantité de d'azote à ajouter X, au printemps, à ladite culture, exprimée kilogramme d'azote par hectare, est calculée selon la formule suivante ; X = [(b.y) + Rf] - (Pi + Ri + M + Mha + Mpro1 + Fleg + Fass), où, 35 b correspond au besoin unitaire en kg d'azote absorbé à la fermeture- 90 - du bilan sur la quantité de graines produites par hectare, y est l'objectif de rendement exprimé en quintaux par hectare aux normes de commercialisation, Pi est la quantité d'azote absorbé par la culture à l'ouverture du bilan (kgN/ha), tel que si le produit MVSH x Coeff SH est supérieur ou égal à MVEH x Coeff EH, alors Pi= MVSH x Coeff SH, sinon, Pi = (MVSH x Coeff SH) +1/2(((MVEH x Coeff EH) - (MVSH x CoeffSH)) / 1,35) où Coeff EH est le coefficient de conversion de la biomasse verte aérienne en kg/m2 en quantité d'azote absorbé en kg d'azote 10 par hectare à l'entrée de l'hiver, et Coeff SH est le coefficient de conversion de la biomasse verte aérienne en kg/ne en quantité d'azote absorbé kg d'azote par hectare à la sortie de l'hiver, et où Rf est la quantité d'azote minéral dans le sol à la fermeture du bilan 15 exprimée en kg d'azote par hectare, Ri correspond à la quantité d'azote minéral dans le sol à l'ouverture du bilan exprimée en kg d'azote par hectare, M correspond à la minéralisation nette de l'azote de l'humus et des résidus de récolte du précédent exprimée en kg d'azote par hectare, 20 Mpro1 est la minéralisation nette de l'azote organique des PRO apporté l'année de la culture avant l'ouverture du bilan exprimée en kg d'azote par hectare, Mpro1=qPRO x NPRO x Keq où qPRO correspond à la quantité de produit résiduaire organique épandu, 25 NPRO correspond à la teneur en azote desdits produits résiduaires, et Keq correspond au coefficient d'équivalence engrais Mha correspond à la fourniture d'azote liée à l'historique d'apport de produits résiduaires organiques 30 Fleg correspond au supplément de fourniture d'azote lié à la présence sur la parcelle de pois protéagineux la campagne précédant la culture de colza exprimée en kg d'azote par hectare, Fass correspond au supplément de fourniture d'azote par la culture lié à des cultures compagnes de légumineuses gélives exprimé en kg 35 d'azote par hectare,- 91 - où N est un sol superficiel (Nss), un sol profond (Np) ou un sol de craie de champagne (Nec), et Rf, Ri et M sont dépendants de N et de L, ladite quantité de d'azote à apporter X, au printemps, à ladite culture, exprimée en kg d'azote par hectare, est calculée selon la formule suivante ; 0 si la région est le Grand Ouest - si le sol est un sol Nss X = [(7y) + 15] - (Pi + 60 + Mha + Mpro1 + Fleg + Fass), - si le sol est un sol Np 10 X = [(7y) 30] - (Pi + 90 + Mha + Mpro1 + Reg + Fass), 0 si la région n'est pas le Grand Ouest - si le sol est un sol Nss X = [(7y) + 15] - (Pi + 40 + Mha + Mpro1 + Fleg + Fass), - si le sol est un sol Nsp 15 X = [(7y) + 30] - (Pi + 70 + Mha + Mpro1 + Fleg + Fass), - si le sol est un sol Nsc X = [(7y) + 30] - (Pi + 60 + Mha + Mpro1 + Fleg + Fass), (ii) si la culture de colza est localisée dans la région Nord-Est ou dans la région Auvergne-Rhône-Alpes, la quantité de d'azote à ajouter X, 20 au printemps, à ladite culture, exprimée en unité d'azote, est calculée selon la formule suivante : X = (((b.y) - (Pi + Np)) / CAU) - (Mha + Mpro1 + Fleg + Fass), où, b correspond au besoin unitaire en kg d'azote absorbé à la fermeture 25 du bilan sur la quantité de graines produites par hectare, y est l'objectif de rendement (q/ha aux normes), Pi est la quantité d'azote absorbé par la culture à l'ouverture du bilan (kgN/ha), tel que si le produit MVSH x Coeff SH est supérieur ou égal à MVEH x Coeff EH, alors Pi= MVSH x Coeff SH, sinon, Pi = 30 (MVSH x Coeff SH) + 1/2(((MVEH x Coeff EH) - (MVSH x CoeffSH)) /1,35) où Coeff EH est le coefficient de conversion de la biomasse verte aérienne (kg/m2) en quantité d'azote absorbé (kgN/ha) à l'entrée de l'hiver, et Coeff SH est le coefficient de conversion de la biomasse 35 verte aérienne (kg/m2) en quantité d'azote absorbé (kgN/ha) à la sortie -92 - de l'hiver, et où Np correspond la fourniture d'azote par le sol pendant la période d'ouverture du bilan en kg d'azote par hectare, Mpro1 est la minéralisation nette de l'azote organique des PRO apporté l'année de la culture avant l'ouverture du bilan exprimée en kg d'azote par hectare, Mpro1=qPRO x NPRO x Keq où qPREO correspond à la quantité de produit résiduaire organique épendu, NPRO correspond à la teneur en azote desdits produits résiduaires, et Keq correspond au coefficient d'équivalence engrais Mha correspond à la fourniture d'azote liée à l'historique d'apport de produits résiduaires organiques, Fleg correspond au supplément d'absorption d'azote par la culture exprimée en kg d'azote par hectare, ledit supplément étant lié à la présence sur la parcelle de pois protéagineux lors de la campagne précédent la culture de colza, Fass correspond au supplément de fourniture d'azote par la culture lié à des cultures compagnes de légumineuses gélives exprimé en kg d'azote par hectare, CAU est le coefficient apparent d'utilisation variant de 0,6 à 1, où N est un sol superficiel (Nss) ou un sol profond (Nss) et Np sont dépendants de N et L, ladite quantité de d'azote à ajouter X, au printemps, à ladite culture, exprimée en unité d'azote, est calculée selon la formule suivante ; o si la région est la région Nord-Est - si le sol est un sol Nss X = (((7y) - (Pi + 30))! 0.8) - (Mha + Mpro1+ Fleg + Fass), - si le sol est un sol NsP X = (((7y) - (Pi + 50))! 0.8) - (Mha + Mpro1+ Fleg + Fass), o si la région est la région Auvergne-Rhône-Alpes - site sol est un sol Nss X = (((7y) - (Pi + 40))! 0.8) - (Mha + Mpro1+ Fleg + Fass), 3027187 - 93 - si le sol est un sol Nsp X = (((7y) - (Pi + 60))! 0.8) - (Mha + Moral + Fleg + Fess), d) une étape de correction du calcul tel que : si le produit (b y) calculé est supérieur à 330, alors le produit (b y) est fixé à 330, si X est supérieur à 250 unités d'azote, alors X est fixé à 250. e) et éventuellement une étape de pesées des engrais azotés apportés correspondant au total à la valeur X obtenue si celle-ci est non nulle. 10 2. Méthode de fertilisation selon la revendication 1, çomprenant en outre une étape d'épandage, au printemps de la quantité X d'engrais à ladite culture de colza. 3. Méthode de fertilisation selon la revendication 1 ou 2, dans laquelle Mha varie de 0 à 100, 15 Keq varie de 0 à 0,4 Fleg et Fass varient de 0à 50, Coeff SH vaut 65, et coeff EH vaut 50. 4. Méthode selon l'une quelconque des revendications 1 à 3, dans laquelle si la culture de colza est implantée sur un terrain ayant porté une culture de pois 20 protéagineux la campagne précédente, Fleg vaut 25, sinon Fleg vaut 0, et si la culture de colza est implantée avec des légumineuses gélives associées, Fass faut 30, sinon Fass vaut 0. 5. Méthode selon l'une quelconque des revendications 1 à 4, dans laquelle Mha est 25 déterminé _comme indiqué dans le tableau suivant : Produits Résiduaires Organiques utilisé tous les 2 années tous les tous les ans sur 3 2 ans 3 ans et . plus Fumier de bovins 30 20 15 5 Fumier de volailles 15 10 5 5 Fumier de porcs 20 15 5 0 Fumier de cheval 40 25 20 10 Fumier d'ovins 35 20 15 10 Lisier de porcs 20 15 10 5 Lisier de bovins 25 15 10 5- 94 - Lisier de bovins dilué 10 5 5 o Lisier, fientes de volailles (560%MS) 10 5 5 o Fientes de volailles séches (80%MS) 20 10 10 o Boues urbaines liquides et pâteuses (520%MS) 15 15 10 5 Boues urbaines sèchées (90%MS) 100 60 50 25 Compost de déchets verts 30 20 15 10 Compost de fumier de bovins, boues STEP + déchets verts 30 20 15 5 Compost de fumier de volailles 25 15 10 10 Vinasse de betterave concentrée o o o o Ecumes de sucreries surpressées 5 5 5 o 6. Méthode selon l'une quelconque des revendications 1 à 5, dans laquelle Mpro1 = teneur en azote des produits résiduaires (NPRO) x quantité de produits résiduaires épandue (qPRO) le coefficient d'équivalence engrais (keq), où NPRO, qPRO et Keq sont définis dans le tableau suivant : Produit Résiduaire Organique NPRO qPRO unité ha) (t ou Keq m3 / Fumiers 5,8 Fumier de bovins 20 t / ha 0,1 Fumier de volailles 25 5 t / ha 0,2 Fumier de porcs 8 20 t / ha 0,1 Fumier de cheval 8 20 t / ha 0,1 Fumier d'ovins 7 20 t / ha 0,1 Lisiers Lisier de porcs à l'engrais 5,8 20 m3 / ha 0,05 Lisier de porcs mixtes 3,5 30 m3 / ha 0,05 Lisier de bovins 4,5 30 m3 / ha 0,15 Lisier de bovins dilué 1,6 30 m3 / ha 0,15 Lisier de canards 7,7 20 m3 / ha 0,05 Fientes Fientes de volailles séches (80%MS) 40 5 t /ha 0,05 Fientes de volailles (60%MS) 24 5 t /ha 0,05 Boues urbaines Boues urbaines liquides 3 30 m3 / ha 0,15 Boues urbaines pâteuses (20%MS) 10 10 t! ha 0,1 Boues urbaines sèçhées (90%MS) 40 10 t! ha 0,1- 95 - Composts Compost de déchets verts 10 10 t / ha 0,05 Compost de fumier de bovins 6,3 15 t / ha 0,12 Compost de fumier de volailles 23 3 t/ha 0,12 Autres Vinasse de betterave concentrée 20 3 t! ha 0,15 Ecumes de sucreries surpressées 4 10 t / ha 0 Autre produit organique "champ libre" "champ libre" t / ha 0,1 et ou champ libre correspond a des valeurs fournies par l'exploitant. 7. Produit programme d'ordinateur sur un support approprié pour la mise en oeuvre de la méthode définie dans l'une quelconque des revendications 1 à 6.