Document Type : Research Paper
Authors
1 Department of Agronomy and Plant Breeding, Faculty of Agriculture, University of Tehran, Karaj, Iran
2 Department of Agroecology, Dryland Agricultural Research Institute, Maragheh, Iran.
3 Department of Agricultural Economics, Faculty of Agricultural Economics and Development, University of Tehran, Karaj, Iran
Abstract
Keywords
. مقدمه
بشر در طول تاریخ، جهت تامین نیازهای اولیه غذایی خود، همواره مجبور به صرف انرژیهای مختلف بوده است
(Erdal et al., 2007). در سالهای اخیر، یکی از دغدغههای اساسی پیش روی فعالان بخش کشاورزی، افزایش مصرف انرژی به واسطه رشد روزافزون جمعیت و متعاقبا افزایش تقاضا برای محصولات کشاورزی میباشد (Javadi & Esfahani, 2023). افزایش جمعیت، تغییرات اقلیمی و کمبود منابع آب و زمین سبب ایجاد مخاطره در امنیت غذایی میشود (Tian et al., 2021). بر اساس گزارش مشترک سازمان خواربار و کشاورزی ملل متحد (FAO)، صندوق بینالمللی توسعه کشاورزی (IFAD)، برنامه جهانی غذای سازمان ملل (WFP) و سایر نهادهای ذیربط، حدود 828 میلیون نفر در سال 2021 تحت تاثیر گرسنگی قرار گرفتند.
این گزارش پیشبینی کرده حدود 670 میلیون نفر (هشت درصد از جمعیت جهان) همچنان در سال 2030 با گرسنگی روبرو خواهند بود. با توجه به بررسی سناریوهای مختلف از جمله ادامه بیثباتی و جنگ در اوکراین، از سوی دیگر اعمال سیاستهای محدودکننده به برخی کشورها همچون تحریم اقتصادی ایران، چگونگی تامین غذا در کشورهای مختلف نگرانکننده است (Unicef, 2022). گندم (Triticum aestivum L.) و جو (Hordeum vulgare L.) بهواسطه قرار گرفتن در جیره غذایی انسان و دام مهمترین محصولات مهم استراتژیک موثر بر امنیت غذایی در سراسر جهان و بهویژه در کشور ایران به شمار میآیند. بر اساس آمار فائو، تولید محصول استراتژیک گندم در سال 2019 در ایران 8/16 میلیون تن بود. نسبت وابستگی به واردات غلات در ایران در سال 2015، 7/35 درصد بوده که این میزان در سالهای 2016 تا 2018، (1/36 درصد) افزایش یافته است (FAO, 2021). بر اساس آمارهای وزارت جهاد کشاورزی، گندم با 46/49 درصد مهمترین محصول تولیدی کشور میباشد (Bayat et al., 2023; Hadi et al., 2017; Nasrabadi et al., 2015). جو پس از گندم، با سطح کشت بیش از 1660 هکتار و تولید 3800 هزار تُن عمدهترین محصول زراعی در کشور میباشد. این محصول بهواسطه سازگاری اکولوژیک در عمده نقاط کشور ایران مورد کاشت قرار میگیرد (Ahmadi et al., 2021). طبق اعلام مرکز آمار ایران در سال 1402 سهم استان تهران از مساحت اراضی زراعی آبی کشور 78/1 درصد (رتبه 20) میباشد؛ اما از جنبه میزان تولید محصولات زراعی سهم چهار درصدی داشته و در رتبه هفتم کشور قرار دارد. این آمار نشان میدهد بهرهبرداران بخش کشاورزی استان نسبت به سایر استانهای کشور در افزایش بهرهوری تولید موفق بودهاند. با این حال اراضی کشاورزی استان تهران با بهرهگیری از شرایط آبوهوایی متنوع، از پتانسیل عملکرد بالایی برخوردار بوده که این موضوع اهمیت نحوه مدیریت صحیح نهادههای تولید را بیشازپیش نشان میدهد. یکی از راهکارهای اصلی تامین و حفظ امنیت غذایی، مدیریت مصرف انرژی در فعالیتهای کشاورزی میباشد (Platis et al., 2019). بر اساس پیشبینیهای انجامشده، جمعیت جهان تا سال 2050 میلادی به 3/9 میلیارد نفر خواهد رسید (Heng et al., 2009). لذا، تامین نیازهای غذایی ناشی از افزایش جمعیت منجر به افزایش مصرف نهادههای شیمیایی و منابع انرژی از جمله سوختهای فسیلی، الکتریسیته، کودهای شیمیایی و آفتکشها میشود (Mohammadzadeh et al., 2017). افزایش مصرف این منابع سبب افزایش میزان انتشار گازهای گلخانهای و متعاقبا مخاطرات زیستمحیطی مختلف میشود (Yousefi et al., 2016). یکی از عوامل موثر بر افزایش بهرهوری انرژی و کاهش اثرات زیستمحیطی انرژیهای ورودی، شناسایی عوامل موثر بر کارایی و بهرهوری انرژی و در عین حال یافتن جایگزین با در نظر داشتن جوانب اقتصادی و فنی میباشد (Vafabakhsh et al., 2019). بر اساس گزارش فائو، در سال 2019 میزان انتشار گازهای گلخانهای در زمینهای کشاورزی ایران، 2/47 میلیون تن CO2 انتشاریافته میباشد (FAO, 2021).
مطالعات متعددی در ایران با رویکرد جریان انرژی، انتشار گازهای گلخانهای و تولید محصولات گندم و جو تاکنون انجام شده است (Heidarisoltanabadi, 2023, Kazemi et al., 2015; Mohammadzadeh et al., 2017, Alipour et al., 2012,
Omid-Mehr, 2015; Javadi & Esfahani, 2023; Vafabakhsh & Mohammadzadeh, 2019; Khodarezaie et al., 2022; Bayat et al., 2023, Khodaei Joghan et al., 2021). بر اساس سامانه اطلاعات آماری استان تهران، استان تهران با وسعتی حدود 13842 کیلومتر مربع و جمعیتی بالغ بر 13 میلیون نفر بین ۳۴ تا ۳۶٫۵ درجه عرض شمالی و ۵۰ تا ۵۳ درجه طول شرقی واقع شده است. این استان در نیمه شمالی ایران قرار دارد. از نظر پستی و بلندی استان تهران در دامنه جنوبی البرز مرکزی قرار دارد. در نواحی مختلف استان تهران به علت موقعیت ویژه جغرافیایی، آب و هوای متفاوتی شکل گرفته است. وجود رشتهکوههای البرز در شمال استان موجب تعدیل آب و هوای استان تهران شده؛ اما از سوی دیگر به دلیل مجاورت با مناطق خشک و کویری استانهای قزوین، قم و سمنان، تاثیرات حرارتی نامطلوبی بر اقلیم استان داشته است. بادهای مرطوب و بارانزای غربی نیز از دیگر عوامل تاثیرگذار بر اقلیم استان میباشد. بر اساس آمار سازمان جهاد کشاورزی استان تهران در سال 1400، مجموع سطح زیر کشت گندم و جو بهترتیب 2/44195 و 5/28709 هکتار میباشد. بهطور کلی، محصولات زراعی مورد کشت در استان تهران شامل گندم، جو، کلزا، ذرت، نباتات علوفهای، سبزی و صیفیجات میباشد که دو گیاه زراعی گندم و جو با توجه به اهمیت استراتژیکی در این استان دارای سهم اراضی قابل توجهی نسبت به سایر محصولات زراعی میباشند.
در این مطالعه، به محاسبه و مقایسه شاخصهای انرژی و میزان انتشار گازهای گلخانهای دو محصول زراعی گندم و جو در شهرستانهای استان تهران پرداخته شده است.
جمعآوری اطلاعات مورد نیاز این مطالعه برای محصولات گندم و جو شامل نوع و میزان نهادههای مصرفی (آب، الکتریسیته، کودهای شیمیایی و آلی زیستی، سموم شیمیایی، سوخت مصرفی) و میزان عملکرد محصول (عملکرد دانه و کاه) به روش پرسشنامه و انجام مصاحبه حضوری با 300 کشاورز در شهرستانهای استان تهران شامل شهرستانهای تهران، اسلامشهر، پاکدشت، دماوند، رباط کریم، ری، شهریار، فیروزکوه، ملارد و ورامین انجام شد. همچنین اطلاعات شهرستانهای پردیس، بهارستان، قدس، قرچک و پیشوا بهترتیب با توجه به اینکه زمانی بخشی از شهرستانهای تهران، رباط کریم، شهریار و ورامین بوده و از نظر کلیه جوانب دارای موقعیت نسبتا یکسانی با شهرستانهای مزبور میباشد نیز در میان اطلاعات شهرستانهای مرجع منظور شد. شهرستان شمیرانات نیز فاقد سطح زیر کشت زراعی میباشد. پس از جمعآوری دادهها، اطلاعات مربوط به هر شهرستان دستهبندی شد. سپس، شاخصهای انرژی شامل کارایی مصرف انرژی، بهرهوری انرژی، انرژی مخصوص و انرژی خالص (بهترتیب رابطه 1 تا 4) بر اساس انرژی نهادههای مصرفی طی عملیات زراعی شامل بذر، کود، آفتکشها، ماشینهای کشاورزی، آب و نیروی انسانی به همراه عملکرد محصول مطابق معادل انرژی آنها (جدول 1)، از طریق نرمافزار اکسل 2016 محاسبه شد
(Vafabakhsh & Mohammadzadeh, 2019).
(رابطه 1)
|
Energy use efficiency= |
(رابطه 2)
|
Energy use productivity= |
(رابطه 3)
|
Specific energy = |
(رابطه 4) |
Net energy = |
اشکال مختلف انرژی شامل انرژی مستقیم (نیروی انسانی، سوخت دیزلی، آب آبیاری و الکتریسیته)، انرژی غیر مستقیم (بذر، کودهای شیمیایی، آفتکشها و ماشینهای کشاورزی)، انرژی تجدیدپذیر (نیروی انسانی، بذر، آب آبیاری) و انرژی تجدیدناپذیر (الکتریسیته، کودهای شیمیایی، سوخت دیزلی، آفتکشها و ماشینهای کشاورزی) میباشد که بر اساس اطلاعات جمعآوریشده از طریق جمع انرژیهای مذکور بهدست آمد (Yilmaz et al., 2005). میزان انتشار گازهای گلخانهای بهواسطه مصرف نهادههای شیمیایی بر اساس ضرایب دیاکسید کربن (CO2)، متان (CH4) و نتیروز اکسید (N2O) در جدول 2 محاسبه شد
(Vafabakhsh & Mohammadzadeh, 2019).
جدول 1. معادل انرژی ورودی و خروجی در بوم نظامهای گندم و جو.
Inputs |
Unit |
Energy equivalents (MJ.unit-1) |
Reference |
Inputs |
Unit |
Energy equivalents (MJ.unit-1) |
Reference |
|
Human labour |
hr |
1.96 |
De et al. (2001) |
Herbicides |
l |
85 |
Kitani et al. (1999) |
|
Machinery |
hr |
62.7 |
Mandal et al. (2002) |
Insecticides |
l |
229 |
Kitani et al. (1999) |
|
Diesel |
l |
47.8 |
Kitani et al. (1999) |
Fungicides |
l |
115 |
Kitani et al. (1999) |
|
Gasoline |
l |
46.3 |
Kitani et al. (1999) |
Electricity |
KWh |
12 |
Kitani et al. (1999) |
|
Nitrogen |
kg |
66.14 |
Hatirli et al. (2006) |
Irrigation water |
m3 |
1.02 |
Acaroglu et al. (1998) |
|
Phosphorus (P2O5) |
kg |
12.44 |
Hatirli et al. (2006) |
Wheat seed |
kg |
15.7 |
Mohammadzadeh et al. (2017) |
|
Potassium (K2O) |
kg |
11.15 |
Hatirli et al. (2006) |
Barley seed |
kg |
14.7 |
Mohammadzadeh et al. (2017) |
|
Straw of wheat and barley |
kg |
12.5 |
Mohammadzadeh et al. (2017) |
Wheat and Barley grain |
kg |
14.7 |
Mohammadzadeh et al. (2017) |
|
جدول 2. انتشار گازهای گلخانهای به ازای مصرف هر واحد نهاده ورودی
Input |
CH4 |
N2O |
CO2 |
Reference |
Input |
CH4 |
N2O |
CO2 |
Reference |
Diesel (L) |
5.20 |
0.7 |
3560 |
Kramer et al. (1999) |
Electricity (kwh) |
0.02 |
8.82 |
61.20 |
Tzilivakis et al. (2005) |
Gasoline (L) |
* |
* |
2320 |
Koga & Tajima (2011) |
Herbicide (kg) |
* |
* |
6300 |
Lal (2004) |
Nitrogen (kg) |
3.70 |
0.03 |
3100 |
Snyder et al. (2009) |
Insecticide (kg) |
* |
* |
5100 |
Lal (2004) |
Phosphate (kg) |
1.80 |
0.02 |
1000 |
Snyder et al. (2009) |
Fungicide (kg) |
* |
* |
3900 |
Lal (2004) |
Potassium (kg) |
1.00 |
0.01 |
700 |
Snyder et al. (2009) |
equivalence factor CO2 |
21 |
30 |
1 |
IPCC (1995) |
مقادیر ورودی و خروجی به همراه انرژی معادل آن در بومنظامهای گندم و جو شهرستانهای استان تهران در جدول 3 نشان داده شده است. بر اساس نتایج، نیروی انسانی مورد نیاز در بومنظامهای گندم و جو در مناطق مختلف استان تهران متغیر بود. کمترین و بیشترین نیروی انسانی به کار گرفتهشده در کشتبومهای گندم بهترتیب در اراضی دیم شهرستان فیروزکوه با 48 ساعت در هکتار (معادل 1/94 مگاژول در هکتار انرژی) و اراضی شهرستان اسلامشهر با 300 ساعت در هکتار (معادل 588 مگاژول در هکتار انرژی) مشاهده شد. از میان مزارع گندم آبی استان تهران اراضی شهرستان ری از کمترین نیروی انسانی مورد نیاز برخوردار بود. همچنین، از لحاظ نیروی انسانی مورد نیاز اراضی دیم شهرستان فیروزکوه و دماوند بهترتیب با 40 و 70 ساعت معادل 1/94 و 149 مگاژول انرژی نیروی کار کمترین و شهرستانهای تهران و پاکدشت با 250 ساعت معادل 490 مگاژول در هکتار انرژی نیروی کار، بیشترین مقدار را به خود اختصاص دادند. نتایج سایر تحقیقات نیز تایید میکند که یکی از دلایل اصلی کاهش نیروی انسانی مورد نیاز برای تولید گندم در اراضی دیم کاهش عملیات خاکورزی، تعدیل کنترل آفات و بیماریها و حذف عملیات آبیاری مزارع میباشد (Kazemi et al., 2016). مقدار نیروی انسانی بهکاررفته برای تولید محصولات کشاورزی مختلف تحت تاثیر عوامل مختلفی از جمله نوع محصول، اندازه مزارع، نوع عملیات زراعی، ضریب مکانیزاسیون و سایر موارد میباشد
(Mohamadzadeh et al., 2017).
در بومنظام گندم، میزان ساعت کار ماشینهای کشاورزی در اراضی دیم شهرستان فیروزکوه، دماوند و اراضی آبی ری 15 ساعت و اراضی آبی شهرستان شهریار 16 ساعت بود و بیشترین مقدار آن در اراضی آبی به شهرستان تهران با 48 ساعت در هکتار معادل 6/3009 مگاژول در هکتار تعلق داشت. یکی از دلایل اصلی افزایش ساعت کار ماشینهای کشاورزی در شهرستان تهران خرد بودن و پراکندگی اراضی میباشد. در بومنظام جو، بیشترین ساعت کار ماشینهای کشاورزی در اراضی شهرستان تهران با 48 ساعت و کمترین میزان آن در اراضی دیم شهرستانهای فیروزکوه و دماوند با 12 و 15 ساعت معادل 4/752 و 2/1003 مگاژول در هکتار مشاهده شد. در میان اراضی آبی، شهرستانهای ری و ورامین با 12 و 14 ساعت معادل 4/752 و 8/877 مگاژول در هکتار عملیات کشاورزی کمتری نسبت به سایر شهرستانها داشتند. کمترین میزان مصرف سوخت به اراضی جو دیم و آبی شهرستانهای دماوند و فیروز کوه مربوط بود. بهطور کلی، خرد و پراکندهبودن اراضی شهرستان تهران و یکپارچگی اراضی شهرستانهای ری و ورامین، نسبت به سایر عوامل دخیل در عملیات کاشت، داشت و برداشت میتواند علت افزایش و کاهش ساعت کار ماشینهای کشاورزی در منطقه باشد. بر اساس نتایج، مصرف سوخت در اراضی گندم شهرستانهای تهران، اسلامشهر و ورامین به یک میزان بوده (300 لیتر در هکتار معادل 14340 مگاژول در هکتار)، از سوی دیگر اراضی شهرستانهای ملارد، شهریار و رباط کریم از شهرستانهای غرب استان تهران نیز با 320 لیتر در هکتار (معادل 15296 مگاژول در هکتار) بیشترین میزان مصرف سوخت را داشتند. کمترین میزان مصرف سوخت نیز بهترتیب متعلق به اراضی دیم شهرستانهای فیروزکوه (75 لیتر معادل 3585 مگاژول در هکتار) و دماوند (80 لیتر معادل 3824 مگاژول در هکتار) میباشد. مصرف سوخت در اراضی آبی جو شهرستانهای پاکدشت (100 لیتر معادل 4780 مگاژول در هکتار) کمترین مقدار و در شهرستان ملارد با 320 لیتر (معادل 16730 مگاژول در هکتار) بیشترین میزان را به خود اختصاص داد.
بیشترین میزان مصرف نیتروژن در بومنظام گندم در شهرستانهای اسلامشهر، تهران، رباطکریم و ری با 300 کیلوگرم معادل 13228 مگاژول در هکتار انرژی مشاهده شد و کمترین میزان آن به مزارع دیم شهرستانهای فیروزکوه و دماوند بهترتیب با 75 و 3/146 کیلوگرم معادل 1/4299 و 9921 مگاژول در هکتار مربوط بود. همچنین، از مزارع آبی استان تهران، شهرستان ملارد کمترین میزان مصرف نیتروژن (175 کیلوگرم) را در بومنظام گندم داشت؛ اما بیشترین میزان ورودی انرژی برای مصرف نیتروژن در مزارع جو به شهرستان ورامین با 6/208 کیلوگرم در هکتار معادل 8/13796 مگاژول انرژی در هکتار متعلق بود. در اراضی تولید جو، مقایسه میزان انرژی ورودی انواع کودهای شیمیایی در استان تهران نشان داد شهرستان شهریار در مجموع با 15587 مگاژول در هکتار انرژی، بیشترین میزان مصرف نیتروژن، فسفر و پتاسیم را به خود اختصاص داد. اراضی دیم شهرستان فیروزکوه کمترین میزان مصرف کودهای شیمیایی یا بهعبارتی کمترین انرژی ورودی (مجموعا 7/5019 مگاژول در هکتار) برای تولید جو را داشت.
میزان آب مصرفی درکشتبومهای گندم و جو استان مورد ارزیابی قرار گرفت. شهرستانهای اسلامشهر و پاکدشت بیشترین میزان آب مصرفی (8/7690 مگاژول در هکتار) را به خود اختصاص داد. کمترین میزان در اراضی فاریاب شهرستان فیروزکوه مشاهده شد.
در مزارع گندم، بیشترین میزان برق مصرفی چاهها در شهرستان اسلامشهر به میزان 2905 کیلووات ساعت معادل 7/34859 مگاژول انرژی در هکتار محاسبه شد. پس از آن شهرستانهای شهریار و رباط کریم بهترتیب با 6/29771و 6/22259 مگاژول انرژی در هکتار میزان انرژی الکتریسته بیشتری نسبت به سایر شهرستانها صرف کرد. کشتبومهای جو نیز همانند اراضی گندم استان از همین الگو تبعیت کرد. یکی از دلایل افزایش انرژی الکتریسیته در مناطق غربی استان تهران برقیبودن اکثر چاههای کشاورزی، همچنین نوع سیستم آبیاری میباشد. در این میان راهاندازی ایستگاههای پمپاژ در کانالهای فرعی آبرسانی نقاط غربی استان از جمله رباطکریم و اسلامشهر در میزان برق مصرفی برای تولید محصولات زراعی بیتاثیر نیست.
بیشترین میزان بذر مصرفی برای تولید گندم به میزان 250 کیلوگرم در هکتار معادل 3925 مگاژول در هکتار انرژی از میان شهرستانهای نیمه غربی استان تهران به شهرستان رباطکریم و از شهرستانهای واقع در شرق استان، به شهرستان فیروزکوه مربوط بود. این میزان میتواند به دلیل در دسترسبودن آب در این شهرستانها باشد. شهرستان فیروزکوه از نظر شرایط اقلیمی جزء مناطق پر آب استان محسوب میشود. شهرستان رباطکریم نیز به دلیل بهرهمندی از کانال آبرسانی محمدیه، ازآب قابل اطمینان برای تولید گیاهان زراعی برخوردار است. این مطالعه نشان داد در اراضی آبی و دیم شهرستان فیروزکوه و دماوند نسبت به سایر شهرستانهای استان کمترین میزان بذر مصرفی جو را داشتند. بیشترین میزان بذر مصرفی نیز متعلق به شهرستان رباطکریم با 220 کیلوگرم در هکتار معادل 3234 مگاژول در هکتار انرژی بود. از سویی کاهش در میزان بذر مصرفی جو در اراضی شرقی استان میتواند بهدلیل تفاوت در نوع رقم و شرایط جوی منطقه باشد. بیشترین مقادیر خروجی انرژی در خصوص عملکرد دانه به مزارع گندم شهرستان ورامین و ری بهترتیب با 6014 کیلوگرم معادل 88406 مگاژول انرژی و 5755 کیلوگرم معادل 84599 مگاژول انرژی تعلق داشت. بههمینترتیب بیشترین میزان تولید کاه گندم نیز متعلق به اراضی این شهرستانها بوده است. برآورد مقدار انرژی خروجی در بومنظام جو نشان داد بیشترین عملکرد دانه جو با 4908 کیلوگرم در هکتار معادل 6/72147 مگاژول در هکتار انرژی در اراضی شهرستان ورامین به دست آمد. غیر از اراضی دیم که به طور طبیعی میزان عملکرد محصول پایینی دارند، در اراضی آبی نیز شهرستان فیروزکوه با 2000 کیلوگرم در هکتار معادل 29400 مگاژول انرژی در هکتار، کمترین میزان عملکرد محصول را به خود اختصاص داد.
شاخصهای انرژی و اشکال انرژی در بومنظام تولید گندم و جو شهرستانهای استان تهران بهترتیب در جدول 4 و 5 نشـان داده شده است. محاسبه انرژیهای ورودی نشان داد بیشترین و کمترین انرژی ورودی در بومنظام گندم و جو بهترتیب برای شهرستان اسلامشهر و فیروزکوه بوده است. از سوی دیگر، بیشترین انرژی خروجی در بومنظامهای گندم متعلق به اراضی شهرستان ورامین و ری بهترتیب 148543 و 5/142148 مگاژول در هکتار و در کشتبومهای جو مربوط به اراضی شهرستانهای ورامین و شهریار بهترتیب 6/72147 و 6/62592 مگاژول در هکتار بود. کمترین میزان انرژی خروجی به اراضی دیم گندم و جو شهرستان فیروزکوه و دماوند تعلق داشت. محاسبه نسبت انرژی خروجی به انرژی ورودی در بومنظامهای گندم و جو شهرستانهای استان تهران نشان داد بیشترین کارآیی مصرف انرژی در مزارع آبی گندم شهرستان ملارد و فیروزکوه بود. همچنین مزارع شهرستانهای دماوند و ورامین بیشترین کارآیی مصرف را در بومنظام جو داشتند (جدول 4). مطالعه وضعیت انرژی خالص در کشتبومهای گندم و جو نشان داد بیلان مصرف انرژی تنها در مزارع جو آبی شهرستانهای اسلامشهر و شهریار و و دیمزارهای جو شهرستان دماوند منفی است. بهعبارتی، انرژی ورودی محصول جو در این شهرستانها بیش از انرژی خروجی بوده است. از سوی دیگر، بیشترین انرژی خالص در اراضی گندم و جو شهرستان ورامین بهترتیب 2/91711 و 9/25191 مگاژول در هکتار بوده است. در مزارع گندم آبی شهرستان اسلامشهر و جو دیم شهرستان دماوند به ازای تولید هر کیلوگرم محصول، انرژی بیشتری نسبت به سایر شهرستانها استفاده شده است.
بررسی اشکال مختلف انرژی در مزارع گندم و جو نشان داد در هر دو محصول مزارع شهرستان اسلامشهر بیشترین انرژی مستقیم و مزارع شهرستان تهران بیشترین انرژی غیر مستقیم را به خود اختصاص داد. بیشترین انرژی تجدیدپذیر و تجدیدناپذیر بهترتیب در مزارع گندم شهرستان رباطکریم (8/12006 مگاژول در هکتار) و شهرستان اسلامشهر (74426 مگاژول در هکتار) بود. در مجموع کمترین انرژی مستقیم، غیر مستقیم و تجدیدناپذیر در مزارع دیم گندم و جو شهرستان فیروزکوه مشاهده شد. همچنین کمترین انرژی تجدیدپذیر استان متعلق به مزارع دیم شهرستان دماوند بود. بهطور کلی، درکشتبومهای گندم و جو استان تهران سهم انرژیهای مستقیم (از جمله انرژی مصرفی نیروی انسانی، سوخت، آب مصرفی و الکتریسته)، بیش از انرژی غیر مستقیم (شامل میزان بذر، کود دامی، مکانیزاسیون، کنترلکنندههای زیستی و کودهای شیمیایی) میباشد (شکل 1 و 2).
جدول 3. مقادیر ورودی و خروجی در بومنظامهای گندم و جو شهرستانهای استان تهران.
Mean±SE
|
Varamin
|
Malard |
Firoozkuh- Rainfed |
Firoozkuh- Irrigated |
Shahriyar |
Ray |
Robat Karim |
Damavand- Rainfed |
Damavand- Irrigated |
Tehran |
Pakdasht |
Eslamshahr |
Product |
Unit |
Inputs |
|
174±72.2 |
164 |
150 |
48 |
170 |
200 |
136 |
200 |
76 |
144 |
250 |
250 |
300 |
Wheat |
hr |
Human labour |
|
144.8±66.8 |
120 |
145 |
40 |
105 |
200 |
120 |
200 |
73 |
107 |
250 |
250 |
128 |
Barley |
|||
20.66±9.2 |
25 |
22 |
15 |
17 |
16 |
15 |
20 |
15 |
15 |
48 |
20 |
20 |
Wheat |
hr |
Machinery |
|
18.8±9.6 |
14 |
18 |
12 |
16 |
16 |
12 |
20 |
15 |
15 |
48 |
20 |
20 |
Barley |
|||
220±108.1 |
300 |
320 |
75 |
120 |
320 |
285 |
320 |
80 |
120 |
300 |
100 |
300 |
Wheat |
L |
Diesel |
|
207.6±105 |
196 |
350 |
70 |
120 |
320 |
270 |
300 |
75 |
120 |
270 |
100 |
300 |
Barley |
|||
232.5±71.3 |
254.1 |
175 |
75 |
240 |
200 |
300 |
300 |
146.3 |
250 |
300 |
250 |
300 |
Wheat |
kg |
Nitrogen |
|
156.8±46.1 |
208.6 |
125 |
65 |
100 |
200 |
200 |
143 |
150 |
150 |
200 |
140 |
200 |
Barley |
|||
107±44.7 |
109.6 |
50 |
50 |
120 |
150 |
150 |
150 |
50 |
100 |
150 |
54 |
150 |
Wheat |
kg |
Phosphate (P2O5) |
|
68.6±24.9 |
45.4 |
50 |
40 |
50 |
100 |
100 |
50 |
50 |
77.8 |
100 |
60 |
100 |
Barley |
|||
55.9±37.8 |
16 |
50 |
20 |
50 |
150 |
50 |
20 |
50 |
100 |
50 |
40 |
75 |
Wheat |
kg |
Potassium (K2O) |
|
50.2±23.8 |
16 |
50 |
20 |
70 |
100 |
50 |
20 |
50 |
66.7 |
50 |
60 |
50 |
Barley |
|||
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
Wheat |
kg |
Calcium |
|
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
Barley |
|||
10±23.7 |
0 |
0 |
0 |
70 |
0 |
0 |
0 |
0 |
50 |
0 |
0 |
0 |
Wheat |
kg |
Sulphur |
|
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
Barley |
|||
2.3±1.3 |
1.5 |
5 |
1.5 |
4 |
2 |
2 |
1.5 |
0 |
4 |
2 |
2 |
2 |
Wheat |
L |
Herbicides |
|
1.2±1 |
1.5 |
3 |
0 |
0 |
2 |
2 |
0 |
0 |
1.6 |
2 |
2 |
1 |
Barley |
|||
4.4±4 |
0.5 |
4 |
0 |
12 |
2 |
8 |
0 |
0 |
5 |
8 |
6 |
8 |
Wheat |
kg |
Micronutrients |
|
0.25±0.9 |
0 |
3 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
Barley |
|||
1.2±1.1 |
1 |
1 |
0 |
3 |
2 |
0 |
1 |
0 |
3 |
0 |
2 |
1 |
Wheat |
L |
Insecticides |
|
0.5±0.7 |
1 |
1 |
0 |
0 |
1 |
0 |
0 |
0 |
0 |
0 |
2 |
1 |
Barley |
|||
0.4±0.5 |
0 |
1 |
0 |
1.5 |
0 |
0.5 |
0 |
0 |
0.5 |
0.5 |
0 |
0.6 |
Wheat |
L |
Fungicides |
|
0.25±0.3 |
0.5 |
0.5 |
0 |
0 |
0 |
0.5 |
0 |
0 |
0 |
0.5 |
0 |
1 |
Barley |
|||
ادامه جدول 3. مقادیر ورودی و خروجی در بومنظامهای گندم و جو شهرستانهای استان تهران.
Mean±SE
|
Varamin
|
Malard |
Firoozkuh- Rainfed |
Firoozkuh- Irrigated |
Shahriyar |
Ray |
Robat Karim |
Damavand- Rainfed |
Damavand- Irrigated |
Tehran |
Pakdasht |
Eslamshahr |
Product |
Unit |
Inputs |
1140±913.1 |
864.5 |
624.1 |
0 |
359.9 |
2481 |
1306.4 |
1855 |
0 |
1000.7 |
867.7 |
1416.2 |
2905 |
Wheat |
KWh |
Electricity |
1055.1±886.1 |
864 |
528.5 |
0 |
298.1 |
2066.4 |
1306.4 |
1855 |
0 |
775.2 |
646 |
1416.2 |
2905 |
Barley |
||
5017.2±2850.7 |
7539 |
5492 |
0 |
3018 |
5889 |
7539 |
7539 |
0 |
4508 |
3603 |
7540 |
7540 |
Wheat |
m3 |
Irrigation water |
4660.5±2937 |
7539 |
4650 |
0 |
2500 |
4905 |
7539 |
7539 |
0 |
3492 |
2682 |
7540 |
7540 |
Barley |
||
195.3±54.6 |
223.5 |
200 |
100 |
250 |
220 |
220 |
250 |
80.4 |
160 |
220 |
200 |
220 |
Wheat |
kg |
Seed |
172.8±45.7 |
217.1 |
200 |
85 |
165 |
180 |
180 |
220 |
80 |
166.7 |
180 |
200 |
200 |
Barley |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Outputs |
4457.4±1724 |
6014 |
5539 |
1000 |
4183 |
5209 |
5755 |
5521 |
896.3 |
4333.2 |
5205 |
4834 |
5000 |
Wheat |
kg |
Grain |
3563.1±1310.8 |
4811 |
4431 |
1050 |
4000 |
4167 |
4604 |
4416 |
1078 |
4000 |
4200 |
3500 |
2500 |
Wheat |
kg |
Straw |
3351.5±1367.6 |
4908 |
3710 |
796 |
2000 |
4258 |
4000 |
4105 |
800 |
4000 |
4000 |
3641 |
4000 |
Barley |
kg |
Grain |
2658.3±1054.2 |
3925 |
3600 |
875 |
2200 |
3406 |
3800 |
3284 |
960 |
3200 |
2500 |
2400 |
1750 |
Barley |
kg |
Straw |
شکل 1. سهم ورودیهای مختلف در کل انرژی ورودی در بومنظام گندم در شهرستانهای استان تهران.
شکل 2. سهم ورودیهای مختلف در کل انرژی ورودی در بومنظام جو در شهرستانهای استان تهران.
شکل 3. سهم ورودیهای مختلف در پتانسیل گرمایش جهانی در بومنظامهای گندم در شهرستانهای استان تهران.
جدول 4. شاخصهای انرژی در بومنظامهای گندم و جو شهرستانهای استان تهران.
Energy productivity (kg.MJ-1) |
Specific Energy (MJ.kg-1) |
Energy use Efficiency |
Net Energy (MJ.ha-1) |
Output Energy (MJ.ha-1) |
Input Energy (MJ.ha-1) |
Product |
|
0.06 |
17.23 |
1.21 |
18591.2 |
104750 |
86159 |
Wheat |
Eslamshahr |
0.1 |
19.2 |
0.8 |
-17993.9 |
58800 |
76793.9 |
Barley |
|
0.09 |
11.03 |
2.15 |
61459.8 |
114809.8 |
53350 |
Wheat |
Pakdasht |
0.1 |
12.5 |
1.2 |
8070.5 |
53522.7 |
45452.2 |
Barley |
|
0.09 |
11.3 |
2.19 |
70178.9 |
129013.5 |
58835 |
Wheat |
Tehran |
0.1 |
11.2 |
1.3 |
14003.4 |
58800 |
44796.6 |
Barley |
|
0.09 |
10.8 |
2.42 |
68871 |
113698 |
46827 |
Wheat |
Damavand Irrigated - |
0.1 |
8.5 |
1.7 |
24833.8 |
58800 |
33966.2 |
Barley |
|
0.05 |
19 |
1.56 |
9619.1 |
26650.6 |
17032 |
Wheat |
Damavand- Rainfed |
0 |
21.2 |
0.7 |
-5185.1 |
11760 |
16945.1 |
Barley |
|
0.07 |
13.24 |
1.86 |
63254.8 |
136358.7 |
73104 |
Wheat |
Robat Karim |
0.1 |
14.5 |
1 |
871.1 |
60343.5 |
59472.4 |
Barley |
|
0.09 |
11.31 |
2.18 |
77044.7 |
142148.5 |
65104 |
Wheat |
Ray |
0.1 |
13.8 |
1.1 |
3636.7 |
58800 |
55163.3 |
Barley |
|
0.07 |
14.12 |
1.75 |
55101.7 |
128659.8 |
73558 |
Wheat |
Shahriyar |
0.1 |
15.3 |
1 |
-2530 |
62592.6 |
65122.6 |
Barley |
|
0.1 |
9.29 |
2.86 |
72631.1 |
111490.1 |
38859 |
Wheat |
Firoozkuh- Irrigated |
0.1 |
11.8 |
1.3 |
5886.2 |
29400 |
23513.8 |
Barley |
|
0.08 |
12.12 |
2.29 |
15702.4 |
27825 |
12123 |
Wheat |
Firoozkuh -Rainfed |
0.1 |
13.1 |
1.1 |
1255.2 |
11701.2 |
10446 |
Barley |
|
0.11 |
8.5 |
2.9 |
89607.4 |
136810.8 |
47203 |
Wheat |
Malard |
0.1 |
11.5 |
1.3 |
12020.7 |
54537 |
42516.3 |
Barley |
|
0.1 |
9.45 |
2.61 |
91711.2 |
148543.3 |
56832 |
Wheat |
Varamin |
0.1 |
9.6 |
1.5 |
25191.9 |
72147.6 |
46955.7 |
Barley |
|
0.08±0.01 |
12.3±3.2 |
2.16±0.5 |
57814.4±28137.9 |
110063.2±40908.7 |
52415.5±22061.7 |
Wheat |
Mean±SE |
0.09±0.03 |
13.51±3.68 |
1.16±0.28 |
5838.4±12267.1 |
49267.05±20103.9 |
43428.7±19710.6 |
Barley |
شکل 3 و 4، سهم ورودیهای مختلف در پتانسیل گرمایش جهانی در بومنظامهای گندم و جو استان تهران را نشان میدهد. بر اساس نتایج، الکتریسیته و نیتروژن مصرفی بیشترین سهم را در انتشار گازهای گلخانهای استان دارند. نتایج حاصل از این تحقیق با مطالعات سایر محققان مطابقت دارد. در این مطالعات بیشترین میزان سهم تولید گازهای گلخانهای به نهادههای سوخت دیزل، کود نیتروژن و الکتریسیته تعلق داشت (Mondani et al., 2017; Taheri Rad et al., 2015; Khoshnevisan et al., 2013;
Javadi & Esfahani, 2023; Mohamadzadeh et al., 2017; Rajabi et al., 2012).
جدول 5. اشکال انرژی در بومنظامهای گندم و جو شهرستانهای استان تهران.
Non-renewable Energy (MJ.ha-1) |
Renewable Energy (MJ.ha-1) |
Indirect Energy (MJ.ha-1) |
Direct Energy (MJ.ha-1) |
Product |
|
74426 |
11732.8 |
28680.3 |
57478.5 |
Wheat |
Eslamshahr |
65912.2 |
10881.7 |
19652.5 |
57141.4 |
Barley |
|
42029.2 |
11320.8 |
23394.8 |
29955.2 |
Wheat |
Pakdasht |
34331.4 |
11120.8 |
15497 |
29955.2 |
Barley |
|
51215.6 |
7619.1 |
29916.6 |
28918 |
Wheat |
Tehran |
38925 |
5871.6 |
20912.6 |
23884 |
Barley |
|
39434.7 |
7392.4 |
24202 |
22625.1 |
Wheat |
Damavand Irrigated - |
27744.1 |
6222.1 |
15156.1 |
18810 |
Barley |
|
15620.3 |
1411.2 |
13058.6 |
3973 |
Wheat |
Damavand- Rainfed |
15626 |
1319.1 |
13217 |
3728.1 |
Barley |
|
61097.1 |
12006.8 |
27466.5 |
45637.4 |
Wheat |
Robat Karim |
48156.6 |
11315.8 |
14791 |
44681.4 |
Barley |
|
53693.5 |
11410.3 |
27847.5 |
37256.3 |
Wheat |
Ray |
44592.4 |
10571 |
18655.4 |
36507.9 |
Barley |
|
63705.3 |
9852.8 |
22091.7 |
51466.4 |
Wheat |
Shahriyar |
57081.5 |
8041.1 |
19635.2 |
45487.4 |
Barley |
|
31522.4 |
7336.6 |
25392.7 |
13466.3 |
Wheat |
Firoozkuh- Irrigated |
18332.5 |
5181.3 |
11445.2 |
12068.6 |
Barley |
|
10458.5 |
1664.1 |
8443.5 |
3679.1 |
Wheat |
Firoozkuh -Rainfed |
9118.1 |
1327.9 |
7021.6 |
3424.4 |
Barley |
|
38167.5 |
9035.8 |
18522.4 |
28681 |
Wheat |
Malard |
34549.1 |
7967.2 |
14417.1 |
28099.2 |
Barley |
|
45312 |
11520.2 |
24107.4 |
32724.7 |
Wheat |
Varamin |
35839.3 |
11116.4 |
19293.9 |
27661.8 |
Barley |
|
43890.2±18826.5 |
8525.2±3701.9 |
22760.3±6487.8 |
29655.1±17109 |
Wheat |
Mean±SE |
35850.7±16765.8 |
7578±3668.2 |
15807.9±4069.6 |
27620.8±16668.4 |
Barley |
بررسی نتایج این تحقیق نشان میدهد مناطق غربی استان تهران از جمله شهرستانهای اسلامشهر، شهریار و رباطکریم بیشترین سهم از انتشار گازهای گلخانهای را داشته است. از سوی دیگر، مناطق شرقی استان از جمله مزارع فاریاب و دیمزارهای واقع در شهرستان فیروزکوه و دماوند کمترین سهم را از انتشار گازهای گلخانهای داشتند.
شکل 4. سهم ورودیهای مختلف در پتانسیل گرمایش جهانی در بومنظام جو در شهرستانهای استان تهران.
بهطور کلی، نتایج نشان داد که کشتبومهای گندم و جو در شهرستانهای اسلامشهر و رباطکریم واقع در غرب استان، بیشترین مصرف نهادههای ورودی را به خود اختصاص دادند. پس از آن، در شهرستانهای واقع در جنوب استان از جمله شهرستانهای ری، پاکدشت و ورامین نسبت به مناطق شمالی و شرقی استان نهادههای ورودی بیشتری برای این محصولات صرف شده است. از سوی دیگر، بیشترین خروجی عملکرد متعلق به مزارع گندم و جو شهرستان ورامین و ری میباشد. اما مقایسه میان بومنظام گندم و جو استان تهران نشان داد مجموع انرژیهای ورودی در مزارع گندم بیش از جو است؛ بهعبارت دیگر، برای تولید گندم به میزان نهادههای بیشتری نسبت به جو نیاز است. اما نسبت انرژی خروجی به ورودی یا همان کارآیی مصرف انرژی در مزارع گندم استان بیش از مزارع جو بود. نتایج این مطالعه نشان داد الکتریسیته و نیتروژن مصرفی بیشترین سهم از انرژیهای ورودی و انتشار گازهای گلخانهای استان را دارند. افزایش بهرهوری انرژی و کاهش انتشار گازهای گلخانهای را میتوان از طریق مدیریت مصرف منابع و همچنین بهبود حاصلخیزی خاک محقق کرد. در مجموع میتوان نتیجه گرفت میزان انرژیهای ورودی در کشتبومهای گندم و جو شهرستانهای استان تهران به دلیل تنوع آب و هوایی، متفاوت بوده؛ به همین دلیل، مدیریت منابع تولید نیازمند شناخت کامل بومنظامهای زراعی میباشد.
Acaroglu, M. (1998). Energy from biomass, and applications. University of Selcuk, Graduate School of Natural and Applied Sciences, Turkey. 43 pp.
Ahmadi, K., Ebadzadeh, H.R., Hatami, F., Mohammadnia Afrozi, S., Esfandiaripour, A., & Abbas Taghani, R. (2021). Agricultural statistics for 2019-2020. Ministry of Jahad-e- Agriculture, Planning and Economic Deputy, Information and Communication Technology Center, 1, 20.
Alipour, A., Afshar, R.K., Nejad, M.K., Behbahani, A.G., & Mohammadi, V. (2014). Evaluation of energy flow in irrigated wheat agro ecosystems. A case study: Shahr-e-Rey City. Journal of Agriculture Science and Sustainable Production, 23, 59-69. (In Persian).
Bayat, F., Sarmadian, F., Jahansuz, M.R., & Sahbani, M. (2023). Estimation of irrigated winter wheat and barley production potential by FAO method. Journal of Climate Research, 1402(53), 21-34. (In Persian).
De, D., Singh, R., & Chandra, H. (2001). Technological impact on energy consumption in rainfed soybean cultivation in Madhya Pradesh. Applied Energy, 70, 193-213.
Erdal, G., Esengün, K., Erdal, H., & Gündüz, O. (2007). Energy use and economical analysis of sugar beet production in Tokat province of Turkey. Energy, 32(1), 35-41.
Hatirli, S.A., Ozkan, B., & Fert, C. (2006). Energy inputs and crop yield relationship in greenhouse tomato production. Renewable Energy, 31, 427-438.
Hadi, M., Jalili, M., & Heris, A. (2017). Assessing the wheat yield under irrigated and rainfed farming and evaluating the possibility of supplemental irrigation of rainfed by water stored in deficit irrigated farming. 3(11), 403-411. (In Persian).
FAO. (2021). World food and agriculture - statistical yearbook 2021. Rome.https://doi.org/10.4060/cb4477en.
Heidari Soltanabadi, M. (2023). Determination of energy production function in irrigated wheat of Isfahan province. Engineering and Energy Management, 11(1), 116-127. (In Persian).
Heng, L.K., Hsiao, T., Evett, S., Howell, T., & Steduto, P. (2009). Validating the FAO AquaCrop model for irrigated and water deficient field maize. Agronomy Journal, 101(3), 488-498.
Houshyar, E., Zareifard, H.R., Grundmann, P., & Smith, P. (2015). Determining efficiency of energy input for silage corn production: An econometric approach. Energy, 93, 2166-2174.
IPCC. (1995). Climate change, the science of climate change. In: Houghton, J.T., Meira Filho, L.G., Callander, B.A., Harris, N., Kattenberg, A., & Maskell, K. (Eds). Intergovernmental panel on climate change. Cambridge: Cambridge University Press.
Javadi, H., & Esfahani, S.M.J. (2023). Energy flow and greenhouse gases emission of crop production systems in South Khorasan province. (In Persian).
Kazemi, H., Alizadeh, P., & Nehbandani, A. (2016). Assessing energy flow in rainfed and irrigated wheat fields of Shahrekourd township under two tillage systems. Journal of Agroecology, 8(2), 281-295. (In Persian).
Khodarezaie, E., Khoshbakht, K., Veisi, H., & Nazari, M.R. (2022). Study of the energy flow and global warming potential (GWP) of alfalfa and maize silage production with different water supply sources (A case study: Qazvin plain). Environmental Sciences. (In Persian).
Kitani, O. (1999). CIGR handbook of agricultural engineering. American Society of Agricultural Engineers, United States of America.
Kramer, K.J., Moll, H.C., & Nonhebel, S. (1999). Total greenhouse gas emissions related to the dutch crop production system. Agriculture, Ecosystems and Environment, 72, 9-16.
Koga, N., & Tajima, R. (2011). Assessing energy efficiencies and greenhouse gas emissions under bioethanol-oriented paddy rice production in northern Japan. Journal of Environmental Management, 92, 967-973.
Khodaei, J.A., Taki, M., & Matoorian, H. (2022). Evaluating energy productivity, greenhouse gas emission, global warming potential and sustainability index of wheat and rapeseed agroecosystems in Khorramshahr. 309-324. (In Persian).
Khoshnevisan, B., Rafiee, S., Omid, M., Yousefi, M., & Movahedi, M. (2013). Modeling of energy consumption and GHG (greenhouse gas) emissions in wheat production in Esfahan province of Iran using artificial neural networks. Energy, 52, 333-338.
Lal, R. (2004). Carbon emission from farm operations. Environment International, 30, 981-990.
Mandal, K., Saha, K., Ghosh, P., Hati, K., & Bandyopadhyay, K. (2002). Bioenergy and economic analysis of soybeanbased crop production systems in central India. Biomass and Bioenergy, 23, 337-345.
Mohammadzadeh, A., Mahdavi Damghani, A., Vafabakhsh, J., & Deihimfard, R. (2018). Ecological– economic efficiency for alfalfa (Medicago sativa L.) and corn silage (Zea mays L.) production systems: Maragheh– Bonab plain, east Azerbaijan province. Journal of Agroecology, 10, 875-895. (In Persian).
Mohammadzadeh, A., Mahdavi, D.A., Vafabakhsh, J., & Deihimfard, R. (2018). Ecological–economic efficiency for alfalfa (Medicago sativa L.) and corn silage (Zea mays L.) production systems: Maragheh–Bonab plain, East Azerbaijan province. (In Persian).
Mondani, F., Aleagha, S., Khoramivafa, M., & Ghobadi, R. (2017). Evaluation of greenhouse gases emission based on energy consumption in wheat Agroecosystems. Energy Reports, 3, 37-45.
Mohammadi, A., Rafiee, S., Jafari, A., Keyhani, A., Mousavi-Avval, S.H., & Nonhebel, S. (2014). Energy use efficiency and greenhouse gas emissions of farming systems in north Iran. Renewable and Sustainable Energy Reviews, 30, 724-733.
Nasrabadi, T., Arab, E., & Pourasghar, F. (2015). Investigating the proportion of wheat planted area in Iran with wheat yield and water demand by focusing on virtual water approach. Journal of Environmental Studies, 41(3), 529-543. (In Persian).
Omid-Mehr, Z. (2016). Evaluating energy flow and greenhouse gas emissions in rainfed wheat production. Cereal Research, 6(3), 353-366. (In Persian).
Platis, D.P., Anagnostopoulos, C.D., Tsaboula, A.D., Menexes, G.C., Kalburtji, K.L., & Mamolos, A.P. (2019). Energy analysis, and carbon and water footprint for environmentally friendly farming practices in agroecosystems and agroforestry. Sustainability, 11(6), 1664.
Rajabi, M.H., Zeinali, E., & Soltani, E. (2012). Evaluation of energy use in wheat production in Gorgan. Journal of Plant Production Research, 19(3), 143-171. (In Persian).
TaheriRad, A.R., Nikkhah, A., Khojastehpour, M., & Nourozieh, S. (2015). Assessing GHG emissions, and energy and economic analysis of cotton production in the Golestan province. 428-445. (In Persian).
Snyder, C., Bruulsema, T., Jensen, T., & Fixen, P. (2009). Review of greenhouse gas emissions from crop production systems and fertilizer management effects. Agriculture, Ecosystems and Environment, 133, 247-266.
Tian, X., Engel, B.A., Qian, H., Hua, E., Sun, S., & Wang, Y. (2021). Will reaching the maximum achievable yield potential meet future global food demand? Journal of Cleaner Production, 294, 126285.
Vafabakhsh, J., & Mohammadzadeh, A. (2019). Energy flow and GHG emissions in major field and horticultural crop production systems (Case study: Sharif Abad plain). Journal of Agroecology, 11(2), 365-382. (In Persian).
Tzilivakis, J., Warner, D., May, M., Lewis, K., & Jaggard, K. (2005). An assessment of the energy inputs and greenhouse gas emissions in sugar beet (Beta vulgaris) production in the UK. Agricultural Systems, 85, 101-119.
Unicef. (2022). The State of Food Security and Nutrition in the World (SOFI) Report-2022.
Yousefi, M., Mahdavi Damghani, A., & Khoramivafa, M. (2016). Comparison greenhouse gas (GHG) emissions and global warming potential (GWP) effect of energy use in different wheat agroecosystems in Iran. Environmental Science and Pollution Research, 23, 7390-7397. (In Persian).
Yilmaz, I., Akcaoz, H., & Ozkan, B. (2005). An analysis of energy use and input costs for cotton production in Turkey. Renewable Energy, 30(2), 145-155.