نموذج «نيكوليت» لنمو الخس: أين تكمن أوجه القصور في المزارع العمودية

يظهر نموذج «نيكوليت» في معظم المراجع الأكاديمية المتعلقة بالتنبؤ بمحصول الخس، كما يتزايد استخدامه بشكل متزايد في وحدات التنبؤ التي تُباع للمشغلين. إنه نموذج جيد. لكنه يُستخدم بشكل روتيني بطرق تضمن خطأه، ونادرًا ما يكمن الفشل في الحسابات — بل في المعلمات والمدخلات وخطوة التحويل في النهاية.

NICOLET B3 هو نموذج ديناميكي غير خطي لزراعة الخس، يتنبأ بالمادة الجافة والمادة الطازجة استنادًا إلى الضوء ودرجة الحرارة وثاني أكسيد الكربون. وفي الممارسة العملية، يحدث انحراف في النتائج لثلاثة أسباب: أولًا، لا يتم أبدًا إعادة ضبط المعلمات التي تمت معايرتها في دفيئة شخص آخر لتناسب الظروف المحلية؛ وثانيًا، لا يتم تمثيل طيف مصابيح LED في النموذج؛ وثالثًا، يتم افتراض تحويل الوزن من الجاف إلى الطازج بدلًا من قياسه مقارنةً ببيانات الحصاد.

ما الذي تحاكيه «نيكوليت B3» فعليًّا

ينتمي نموذج «نيكوليت» إلى عائلة من نماذج توازن الكربون الآلية. تقسم صيغة سلالة Seginet النبات إلى حجرات افتراضية — «الهيكل»، و«الفجوة»، و«الكربون الزائد» — حيث تتكامل المركبات العضوية والمعدنية القابلة للذوبان داخل الفجوة مع بعضها البعض للحفاظ على ضغط تناضحي إجمالي ثابت. وهذه النقطة الأخيرة أكثر أهمية مما تبدو عليه: فهي الآلية التي تسمح للنموذج بتقديم معلومات عن محتوى النترات، وليس فقط عن الكتلة.

تتتبع صيغة «فان هنتن» المرتبطة ارتباطًا وثيقًا، والتي تُطبقها معظم البرامج ، متغيرين للحالة، هما الوزن الجاف غير الهيكلي والوزن الجاف الهيكلي، اللذين يتأثران بالإشعاع الساقط النشط ضوئيًا، ودرجة الحرارة، وثاني أكسيد الكربون. ثلاثة مدخلات، حالتان، وعدد قليل من المعادلات التفاضلية العادية. وهذا البساطة هي السبب في أنها تعمل على وحدة تحكم بدلاً من مجموعة وحدات معالجة الرسومات (GPU).

يعمل النموذج بشكل صحيح بعد معايرته. أشار تقييم نُشر عام 2022 في «المجلة الإسبانية للبحوث الزراعية» إلى أن معامل R² بلغ 0.9939 للمادة الطازجة و0.9858 للمادة الجافة بعد المعايرة باستخدام بيانات الدفيئة، ووجد أن درجة الحرارة كان لها التأثير الأكثر وضوحًا على أداء النموذج مقارنةً بالضوء النشط بيولوجيًا (PAR) وثاني أكسيد الكربون (CO₂)، كما أشار إلى أن النموذج يمكن أن يحقق أداءً جيدًا باستخدام مدخلين فقط.

اقرأ ذلك بعناية، لأنه يحتوي على تحذير. فقد جاءت تلك الملاءمات بعد معايرة المعلمات استنادًا إلى البيانات التجريبية المحلية. وقد حددت دراسة الحساسية التي أُجريت على نموذج NICOLET B3 ونُشرت في مجلة «Agrociencia» ثلاث معلمات لها تأثير كبير على سلوك النموذج: معامل التنفس الصياني، ومعامل انقراض الإشعاع عبر الغطاء النباتي، ومعامل معدل النمو غير المقيد للغطاء النباتي المغلق. وهذه المعلمات الثلاثة هي بالضبط تلك التي تتغير عند تغيير الصنف أو كثافة النباتات أو هندسة الرفوف.

لماذا ينحرف نموذج «نيكوليت» المعاير في مزرعة عمودية؟

لأن مجموعات المعلمات المنشورة تم تحديدها في ظل أشعة الشمس، وفي الصوبات الزراعية، مع تباعد بين النباتات مشابه لتلك الموجودة في الحقول — في حين أن المزرعة العمودية تخالف هذه الشروط الثلاثة جميعها في آن واحد.

الطيف غير مرئي بالنسبة للنموذج. يعتبر نموذج NICOLET معدل PAR كقيمة عددية. أظهرت مجموعة من جامعة ولاية ميشيغان وجامعة فرجينيا للتكنولوجيا تعمل مع مختبر إريك رانكل أنه مع وجود معامل كفاءة استخدام الضوء الثابت، لم يتنبأ نموذج النمو الديناميكي بدقة بالكتلة الجافة عبر معالجات الضوء خارج ظروف الدفيئة النموذجية، وأن ملاءمة معامل يعتمد على الطيف خفضت خطأ التنبؤ إلى ما يقارب 6–8% للأشكال البسيطة من الدرجة الأولى والثانية، وحوالي 2% كخطأ متوسط للشكل من الدرجة الرابعة. نفس النموذج، نفس المحصول؛ وكان الفرق يكمن في السماح لمصطلح الكفاءة بالتغير مع الطيف.

هندسة الغطاء النباتي تختلف عن هندسة الدفيئة. يعبر معامل انقراض الإشعاع عن كيفية تلاشي الضوء عبر الغطاء النباتي. وفي حالات الكثافة الخاصة بالزراعة الرأسية — مع ارتفاع ثابت للوحدات المعلقة وغياب مكون الإشعاع السماوي المنتشر — يتصرف الانقراض بشكل مختلف. إن الإبقاء على القيمة الواردة في الأدبيات دون تغيير يؤدي إلى تحيز منهجي في النصف الثاني من دورة المحصول، عندما يهيمن إغلاق الغطاء النباتي.

المناخ المحلي يختلف عن مناخ الغرفة. يعتمد النموذج على ظروف الغطاء النباتي. ويقوم معظم المشغلين بتزويده بقراءات مستشعر موجود في الغرفة على بعد مترين من النباتات. في رف مكون من ثمانية طوابق، لا تسود نفس درجة الحرارة أو نفس قيمة الضغط البخاري (VPD) في الطابق العلوي والطابق السفلي، وتعد درجة الحرارة هي المتغير الذي يتأثر به النموذج بشكل أكبر.

لا يتناول النموذج مسألة «احتراق الأطراف». فبرنامج «نيكوليت» (NICOLET) يتنبأ بالكتلة الحيوية، وليس بالقابلية للتسويق. ويقرن دليل كورنيل لتقييم تكلفة الإنتاج (CEA) الخاص بالخس هدفه المحدد لـ DLI البالغ 17 مول/م²/يوم بتدفق هواء هابط فوق النباتات، ويشير إلى أنه بدون هذا التدفق الهوائي، يتعين زراعة الخس في ظروف إضاءة مخفضة — على سبيل المثال 12 بدلاً من 17 مول/م²/يوم — مما يؤدي إلى إبطاء النمو. يمكن للنموذج أن يتنبأ برأس خس جميل يزن 160 غرامًا، لكن خط مراقبة الجودة الخاص بك يقوم بتخفيض تصنيفه.

ما هو الطراز الذي يتوافق مع أجهزة الاستشعار التي تمتلكها بالفعل؟

ينبغي أن يستند اختيار النموذج إلى الأدوات المستخدمة، وليس إلى المكانة. إن استخدام تركيبة تراعي مستويات النترات في مزرعة لا تُجرى فيها اختبارات الأنسجة يؤدي إلى دقة لا يمكن التحقق منها.

تقوم النماذج الميكانيكية بالاستقراء؛ بينما تقوم النماذج الإحصائية بالتقدير. إذا كنت بصدد إنشاء غرفة جديدة، فستحتاج إلى النماذج الميكانيكية. وهذه هي نفس المفاضلة التي تم تناولها في مقال «الزراعة في بيئة خاضعة للتحكم باستخدام التوأم الرقمي: ما يتطلبه الأمر فعليًّا».

تحويل الوزن الطازج هو النقطة التي عادةً ما تفشل عندها التوقعات

العملة الأساسية في نظام NICOLET هي الكربون. أما العملة التي يستخدمها المشتري فهي الوزن الطازج في عبوة قابلة للفتح. ويتم التحويل بينهما بناءً على افتراض محتوى المادة الجافة، وهو الرقم الأقل خضوعًا للتدقيق في معظم عمليات النشر.

استخدم معيار كورنيل كنقطة ارتكاز: فقد أدى الحفاظ على معدل 17 مول/م²/يوم إلى إنتاج رأس قابل للتسويق يبلغ وزنه الطازج 150 غرامًا مع ما يقرب من 7 غرامات من الوزن الجاف بعد 35 يومًا من البذر — أي ما يعادل 4.7٪ تقريبًا من المادة الجافة. والآن قم بإجراء الحسابات على نموذج يحدد الكتلة الجافة عند 7 غرامات ولكنه يستخدم معدل تحويل خاطئ. عند 5.5% من المادة الجافة، تبلغ التوقعات 127 غرامًا؛ وعند 4.0% تبلغ 175 غرامًا. تم توقع العوامل البيولوجية بشكل صحيح، لكن توقعات المحصول النهائي تختلف بمقدار الثلث عبر هذا النطاق. محتوى المادة الجافة ليس ثابتًا: فهو يتغير مع مؤشر DLI، ودرجة الحرارة، والصنف، ووقت الحصاد خلال اليوم.

كما أن الجانب التناضحي لنظام «نيكوليت» يتمتع بميزة تجارية أيضًا. ونظرًا لأن النترات الفجوية والكربوهيدرات تحل إحداهما محل الأخرى، فإن الإنتاج في ظروف الإضاءة المنخفضة يميل إلى ارتفاع مستويات النترات. إذا كنت تقوم بالشحن إلى الاتحاد الأوروبي، فإن لائحة المفوضية (EU) رقم 1258/2011 تحدد الحد الأقصى لمستوى النترات في الخس المزروع تحت غطاء بـ 5,000 ملغ NO3/كجم للمنتجات التي يتم حصادها في الفترة من 1 أكتوبر إلى 31 مارس، و4,000 ملغ/كجم في الفترة من 1 أبريل إلى 30 سبتمبر. إن النموذج الذي يتتبع الفجوة الخلوية يوفر لك معلومات حول مخاطر عدم الامتثال، وليس فقط حول المحصول.

الحل ليس مثيرًا للإعجاب: قم بوزن وتجفيف عينة من كل محصول، وقم بتسجيلها مقارنةً بالسجل المناخي الفعلي للدفعة، ثم أعد المعايرة. فالنموذج الذي يفتقر إلى حلقة المعايرة ما هو إلا جدول بيانات يحتوي على معادلات تفاضلية عادية.

كيف تتعامل شركة AGEYE مع هذا الأمر

يُعد CultivAid Predict (predict.ageyetech.com) التوأم الرقمي الزراعي ومحرك التنبؤ الخاص بشركة AGEYE، ويضم ستة نماذج معادلات تفاضلية عادية (ODE) للمحاصيل بدلاً من نموذج واحد: نموذج NICOLET B3 للخس، ونموذج Van Henten للخس، ونموذج Xu-Du المراعي لـ VPD للخس، ونموذج TOMGRO للطماطم، ونموذج SIMPLE الشامل، ونموذج الأعشاب اللوجستي. يتم اختيار النموذج من خلال توجيه يعتمد على أجهزة الاستشعار، حيث يصنف التدفقات الواردة إلى أجهزة استشعار كاملة، أو أجهزة استشعار أساسية، أو تحسين المناخ، ثم يختار أفضل نموذج لكل محصول — وبذلك لا تُجبر المزرعة التي تمتلك بيانات VPD على مستوى الغطاء النباتي على استخدام نفس الصيغة التي تستخدمها المزرعة التي تسجل درجة الحرارة فقط.

تسجل نقطة نهاية المعايرة البيانات الفعلية للحصاد وتتتبع معدل الخطأ المطلق (MAPE) مقارنةً بالتنبؤات السابقة، وهو ما يحول النموذج النظري الثابت إلى نموذج يعكس ظروف غرفتك. تتضمن التوقعات حدود ثقة بارامترية وكشفًا عن الحالات الشاذة، وتتيح لك نقطة نهاية السيناريو اختبار التغييرات الافتراضية في الظروف قبل تطبيق الوصفة. تدعم المنصة 16 صنفًا من المحاصيل تشمل الخس والأعشاب والكرنب والطماطم والفراولة والخضروات الصغيرة، بالإضافة إلى محاكاة مزرعة متعددة العوامل تضم 6 عوامل و18 نوعًا من الإجراءات للتخطيط التشغيلي.

ماذا يعني ذلك

تتجه نماذج الخس نحو التصاميم الهجينة: نوى آلية للاستقراء، مع طبقات تصحيح تجريبية للطيف الضوئي والرطوبة ومساحة الورقة النوعية — وهي عوامل لم تكن الصيغ الأصلية التي وُضعت في التسعينيات مضطرة للتعامل معها قط. والأبحاث المنشورة التي تشير إلى هذا الاتجاه موجودة بالفعل. وما تفتقر إليه معظم المزارع ليس معادلة أفضل — بل الانضباط في تسجيل بيانات الحصاد الفعلية الذي يتيح تصحيح أي معادلة. سيحصل المشغلون الذين يقيسون ظروف الغطاء النباتي ويسجلون المادة الجافة في كل حصاد على تنبؤات مفيدة من نموذج عمره ثلاثون عامًا. أما المشغلون الذين لا يفعلون ذلك، فسيحصلون على أرقام موثوقة من أي نموذج يشترونه لاحقًا، وستظل هذه الأرقام خاطئة. تعد دقة التنبؤ مدخلاً مباشرًا لتخطيط العمالة والتعبئة الذي تمت مناقشته في مقال «كيفية خفض التكلفة لكل رطل في الزراعة الداخلية: العوامل الأربعة المهمة».

إذا كنت ترغب في التحقق مما إذا كانت بيئة الإضاءة الفعلية لديك تتطابق مع الافتراضات التي يستند إليها نموذج النمو الخاص بك، فابدأ باستخدام «حاسبة DLI» المجانية من AGEYE.