Investigate different nutritional regimens on some forage quantitative and qualitative traits of Dracocephalum kotschyi Boiss under water deficit stress

مقدمه

مدیریت کود، یکی از مهم‌ترین عوامل در کشت موفق محصولات زراعی است که بر کیفیت و کمیت تولید اثر می‌گذارد (Tahmasebi et al., 2011). استفاده بیش از حد از کودهای شیمیایی، باعث ایجاد مشکلات متعددی در بخش کشاورزی از جمله تغییر در ساختار خاک، آلودگی آب‌های زیرزمینی و سمیت عناصر سنگین شده است (Behrooz et al., 2017) و استفاده از کودهای بیولوژیک مانند ورمی کمپوست و نیتروکسین که علاوه بر تأمین نیاز غذایی گیاه و تضمین عملکرد بالا، سبب افزایش کیفیت محصولات شده‌اند و اثرات مخرب زیست محیطی هم ندارند، می‌توانند جایگزین مناسبی برای کودهای شیمیایی باشند (Barker & Bryson, 2006). به دلیل اثرات مثبت مواد آلی بر خاک، به‌عنوان یکی از ارکان مهم بهره‌وری خاک شناخته شده‌اند (Behrooz et al., 2017). امروزه استفاده از کود کمپوست در اراضی کشاورزی، به‌طور عمومی مورد توجه است که از آن به‌عنوان بهترین تدبیر زیست محیطی عملی یاد شده است (Kabirinegad et al., 2009). ورمی‌کمپوست، کود آلی مناسبی برای تحقق این شرایط (Behrooz et al., 2017) و منبع مناسبی از عناصر غذایی و مواد تنظیم‌کننده رشد گیاهی (هیومیک اسیدها و تنظیم کننده رشد گیاهی مانند اکسین‌ها، جیبرلین‌ها و سیتوکنین‌ها) است که مسئول افزایش رشد گیاه و عملکرد بسیاری از محصولات زراعی هستند (Atiyeh et al., 2002). ورمی‌کمپوست به‌طور مؤثری سبب افزایش شکل‌گیری ریشه، افزایش طول ساقه و تولید زیست‌توده در علف لیمو می‌شود (Srinivas et al., 2017) و تأثیر مثبتی روی بسیاری از گیاهان معطر و دارویی دارد (Anwar et al., 2005). آزولا یک سرخس آبزی و شناور در آب است که همراه با جلبک سبز – آبی (Anabaena azollae) به صورت همزیست زندگی می‌کند و توانایی جذب نیتروژن اتمسفر و تثبیت آن را دارد و می‌تواند جایگزین بسیار خوبی برای کودهای شیمیایی باشد (Yousefzadeh et al., 2011). آزولا منبع بسیار مهمی از عناصر مورد نیاز گیاهان مانند نیتروژن، فسفر، پتاسیم و ... است که پس از اضافه شدن به خاک می‌تواند آن‌ها را به تدریج در اختیار گیاه قرار دهد (Pabby et al., 2003). یکی از روش‌های استفاده از آزولا در خاک‌های زراعی، تبدیل آن به کمپوست (آزوکمپوست) و سپس اضافه کردن آن به خاک است (Yousefzadeh et al., 2011). کودهای زیستی شامل انواع مختلف ریزجانداران آزادزی هستند که توانایی تبدیل عناصر غذای اصلی، از شکل غیرقابل دسترس به شکل قابل دسترس در فرایندهای زیستی دارند (Vessey, 2003) و سبب گسترش نظام‌های ریشه‌ای می‌شوند (Chen, 2006).

زرین‌گیاه با نام علمی Dracocephalum kotschyi Boiss از تیره نعناعیان^[1] می‌باشد که معمولاً در ارتفاعات کوهستانی استان‌های گرگان، مازندران، همدان، کرمانشاه، فارس (کوه دنا) و ارتفاعات تهران (کوه‌های البرز)، ارتفاعات شمال سمنان و اصفهان می‌روید (Mozaffarian, 2015). گونه D.kotschyi (2n=20) از گونه‌های انحصاری و بومی ایران است و با نام زرین‌گیاه یا بادرنجبویه دنایی شناخته می‌شود (Mozaffarian, 2003). به دلیل برداشت بی‌رویه این گیاه در مرحله گلدهی توسط افراد بومی و محدود بودن پراکنش جغرافیایی آن، یکی از گیاهان در معرض انقراض (^[2]EN) ایران می‌باشد (Abdi, 2008). بیشترین وزن خشک، درصد اسانس و عملکرد اسانس در گیاه بادرشبی (Dracocephalum moldavica  L.) از خانواده نعناعیان با کاربرد 100 کیلوگرم نیتروژن در هکتار به‌دست آمده است (Soroori et al., 2014).

ارزش غذایی علوفه در مورد گیاهان علوفه‌ای اهمیت بسیاری دارد. تعیین ماده خشک در گیاه، به‌تنهایی برای ارزیابی کیفی گیاهان علوفه‌ای کفایت نمی‌کند، بلکه میزان قابلیت هضم ماده خشک، کربوهیدارت‌های محلول در آب، الیاف نامحلول در شوینده خنثی، پروتئین و فیبر خام، الیاف نامحلول در شوینده اسیدی و درصد خاکستر نیز از لحاظ قابلیت هضم و میزان انرژی مورد نیاز می‌باشد. یکی دیگر از مهم‌ترین پارامترهای کیفیت علوفه و نماینده مهم‌ترین منبع انرژی در جیره، کربوهیدرات‌های محلول در آب می‌باشد (Coleman & Moore, 2003). الیاف نامحلول در شوینده خنثی[3] (NDF) و الیاف نامحلول در شوینده اسیدی[4] (ADF) به‌ترتیب نشان‌دهنده قابیلت مصرف علوفه توسط دام و قابلیت هضم می‌باشند (Contreras-Govea et al., 2009). محتوی الیاف نامحلول در شوینده اسیدی، نشان‌دهنده سهم دیواره سلولی در علوفه است که شامل سلولز و لیگنین می‌باشد و به‌طور معمول، با افزایش این شاخص، قابلیت هضم علوفه کاهش پیدا می‌کند (Albayrak et al., 2011). کاربرد مواد آلی و کودهای زیستی و شیمیایی روی گیاه سورگوم علوفه‌ای نشان داد که ارزش هضمی، قابلیت هضم ماده آلی و درصد خاکستر، تأثیر معنی‌داری نداشت، ولی پروتئین سورگوم تحت تأثیر قرار گرفت (Saeid-Nejad et al., 2012).

با توجه به این‌که به زرین گیاه در منطقه آذربایجان، قوچ اوتی (علف گوسفند نر) می‌گویند و دلیل این نوع نام‌گذاری، رغبت گوسفند برای خوردن این گیاه می‌باشد، بررسی صفات کیفی علوفه‌ای این گیاه، به نظر لازم و ضروری می‌آید و از طرفی با توجه به این‌که امروزه استفاده از گیاهان دارویی در جیره غذایی دام تحت عنوان علوفه-دارو، موضوع جدیدی در تکمیل جیره غذایی دام به‌حساب می‌آید، اهمیت موضوع مورد بررسی بیشتر می‌شود؛ همچنین شرایط کمبود آب که در کشور حاکم است و پیداکردن گونه‌های مقاوم به شرایط تنش کم‌آبی، اهمیت موضوع را بالاتر می‌برد که در این باره، تحقیق درباره گیاهان بومی ایران را می‌توان در اولویت قرار داد. در این پژوهش، اثر سیستم‌های مختلف تغذیه‌ای بر برخی صفات کمی، کیفی و علوفه‌ای گیاه دارویی زرین گیاه تحت شرایط تنش کم‌آبی مورد بررسی قرار گرفت.

مواد و روش‌ها

این پژوهش به صورت اسپلیت پلات و در قالب بلوک‌های کامل تصادفی در سه تکرار در سال زراعی 97-1396 در مزرعه پژوهشی دانشکده کشاورزی دانشگاه تربیت مدرس، واقع در 17 کیلومتری غرب تهران، با مختصات جغرافیایی 35 درجه و 74 دقیقه عرض شمالی و 51 درجه و 16 دقیقه طول شرقی با ارتفاع 1275 متر از سطح دریا اجرا شد. عامل تنش کم‌آبی با چهار سطح (آبیاری مطلوب، آبیاری بعد از تخلیه 20 درصد آب قابل استفاده؛ آبیاری ملایم، آبیاری بعد از تخلیه 40 درصد آب قابل استفاده؛ آبیاری متوسط، آبیاری بعد از تخلیه 60 درصد آب قابل استفاده و آبیاری شدید، آبیاری بعد از تخلیه 80 درصد آب قابل استفاده) در کرت‌های اصلی و عامل سیستم‌های کودی با پنج سطح (کود شیمیایی اوره، کود زیستی نیتروکسین، ورمی‌کمپوست، آزوکمپوست و عدم کاربرد کود) در کرت‌های فرعی قرار گرفتند. قبل از انجام آزمایش، ویژگی‌های فیزیکی و شیمیایی خاک محل آزمایش در دو عمق صفر تا 30  سانتی­متری و 30 تا 60 سانتی‌متری تعیین شد (جدول 1).

مقدار کاربرد از هر نوع کود با توجه به نیاز کودی زرین گیاه محاسبه شد. کودهای آلی ورمی‌کمپوست و آزوکمپوست قبل از کاشت با خاک مخلوط شدند و میزان مصرف کودهای ورمی‌کمپوست و آزوکمپوست بر اساس میزان آزادسازی نیتروژن (54/30 درصد برای ورمی‌کمپوست و 40/31 درصد برای آزوکمپوست) در طول فصل زراعی و با توجه به خصوصیات نیتروژن موجود در کودهای آلی به‌طور جداگانه برآورد شد و برای مصرف در این پژوهش، 31476 کیلوگرم ورمی‌کمپوست در هکتار و 20276 کیلوگرم آزوکمپوست در هکتار تعیین شد. خصوصیات کودهای ورمی‌کمپوست و آزوکمپوست در جدول 2 قابل مشاهده است. کود زیستی نیتروکسین با توجه به توصیه شرکت سازنده آن (فناوری زیستی مهر آسیا) استفاده شد، به این ترتیب که هنگام کاشت نشاها، ریشه نشاها به محلول نیتروکسین (دو لیتر در 10 لیتر آب) در سایه آغشته و کشت شدند. همچنین کاربرد نیتروژن به روش شیمیایی به فرم اوره به صورت سرک و در سه مرحله ساقه‌دهی، شروع گلدهی و گلدهی، به‌ترتیب 25، 50 و 25 درصد کل نیتروژن صورت گرفت.

جدول 1- خصوصیات فیزیکی و شیمیایی خاک مزرعه

Table 1. Farm soil physicochemical properties

جدول 2- خصوصیات ورمی‌کمپوست و آزوکمپوست

Table 2. Vermicompost and Azocompost properties

مراحل آماده­سازی زمین با اجرای عملیات شخم، دیسک­زنی برای خردکردن کلوخه­ها، تسطیح و ایجاد جوی و پشته توسط فاروئر قبل از کشت در اسفند 1396 و فروردین ماه 1397 انجام شد و کشت به صورت نشاکاری بود. برای این منظور، بذرهای زرین‌گیاه اکوتیپ فریدونشهر، سه ماه قبل از انتقال نشاء به مزرعه خریداری شدند و در گلخانه دانشکده کشاورزی دانشگاه تربیت مدرس، به تعداد کافی در محیط 80 درصد کوکوپیت و 20 درصد پیت ماس کاشته شدند. طول هر واحد آزمایشی، سه متر و عرض آن، دو متر بود و فاصله بین ردیف‌ 50 سانتی‌متر و فاصله روی ردیف‌ها 25 سانتی‌متر در نظر گرفته شد. کشت نشاها در اردیبهشت ماه سال 1397 در مزرعه صورت گرفت. به‌منظور جلوگیری از تداخل تیمارهای آزمایشی و سایر مراحل اجرای آزمایش، فاصله­ای به اندازه 5/1 متر بین بلوک­ها و کرت‌های اصلی و یک متر بین کرت­ها لحاظ شد. آبیاری به صورت قطره­ای و با استفاده از نوارهای آبیاری انجام گرفت. مبارزه با علف‌های‌هرز به‌وسیله وجین کن و با دست در طول فصل زراعی انجام شد و پس از استقرار کامل نشاها (هشت برگی)، تیمار­های آبیاری اعمال شد. برای تعیین سطوح مختلف آبیاری، از روابط ارائه شده توسط  Behera and Panda(2009) استفاده شد در این روش، برنامه زمان­بندی آبیاری بر اساس درصد تخلیه آب خاک در منطقه ریشه است. مقدار آب خاک با استفاده از دستگاه انعکاس سنجی زمانی (TDR^[5]) در عمق‌ گفته شده تعیین و از منحنی‌های کالیبراسیون رطوبتی خاک برای تعیین رابطه بین مقدار عددی ارائه شده توسط دستگاه و مقدار حجمی رطوبت خاک استفاده شد.

به­منظور بررسی علمی تأثیر تیمارهای اعمال شده، برداشت نهایی در گلدهی کامل انجام گرفت. بدین منظور، بوته­ها از فاصله پنج سانتی­متری سطح خاک و با رعایت حذف اثر حاشیه‌ای یه اندازه یک متر مربع کف­بر شدند و جهت اندازه­گیری صفات مورد بررسی به آزمایشگاه انتقال یافتند. بعد از برداشت، عملکرد تر اندام‌های هوایی اندازه­گیری شد و سپس نمونه‌ها سایه خشک شدند. در محله بعدی، نمونه‌های خشک شده به‌طور کامل آسیاب شدند و برای اندازه‌گیری صفات کیفی علوفه‌ای به آزمایشگاه مؤسسه تحقیقات جنگل‌ها و مراتع منتقل شدند. اندازه‌گیری صفات با استفاده از دستگاه طیف سنجی مادون قرمز ([6]NIR) و با روش ارائه شده توسط Jafari et al. (2003) انجام گرفت. صفات علوفه‌ای اندازه‌گیری شده با این دستگاه شامل درصد ماده خشک قابل هضم (DMD[7])، درصد قندهای محلول در آب (WSC[8])، درصد پروتئین خام ([9]CP)، درصد دیواره سلولی منهای همی سلولز یا درصد الیاف نامحلول در شونده اسیدی (ADF[10])، درصد خاکستر کل[11] (ASH)، درصد فیبر خام (CF[12]) و درصد الیاف نامحلول در شوینده خنثی (NDF[13]) بود.

برای تجزیه‌ و تحلیل داده‌ها، بر اساس روش تجزیه واریانس، از برنامه آماری (version 9.3) SAS ^[14]  و برای ترسیم نمودارها از نرم افزار13.0  Excel استفاده شد و همچنین مقایسه میانگین‌های صفات مورد بررسی به روش آزمون LSD^[15] در سطح احتمال پنج درصد انجام گرفت.

نتایج و بحث

وزن تر برگ، ساقه و کل

نتایج جدول تجزیه واریانس نشان داد که وزن تر برگ، ساقه و کل زرین گیاه، تحت تأثیر کاربرد رژیم‌های کودی در شرایط تنش کم‌آبی در سطح احتمال یک درصد قرار گرفت (جدول 3). علت عدم وجود تفاوت در نسبت برگ به ساقه را می‌توان به تأثیر تقریباً یکنواخت تنش اعمال شده بر تجمع ماده خشک و کاهش وزن در کلیه اندا‌م‌های هوایی گیاه نسبت داد. جدول مقایسه میانگین برهمکنش رژیم کودی و رژیم‌های آبیاری نشان داد که بیشترین وزن تر برگ با 1797 کیلوگرم در هکتار، در شرایط تنش ملایم و با کاربرد کود شیمیایی اوره به‌دست آمد (جدول4) وبیشترین وزن تر ساقه برگ از تیمارهای کود شیمیایی در شرایط تنش ملایم و کاربرد کود زیستی نیتروکسین در شرایط تنش شدید به‌ترتیب با 1161 و 1050 کیلوگرم در هکتار تولید شد (جدول 4). در بررسی جدول مقایسه میانگین برهمکنش عوامل مورد بررسی، مشاهده شد که کاربرد کود شیمیایی اوره در شرایط تنش کم‌آبی ملایم با میزان زیست‌توده 2957 کیلوگرم در هکتار، بیشترین وزن تر کل را تولید کرد (جدول 4).

به دلیل این‌که وزن تر کل در گیاهان دارویی، همبستگی مثبتی با عملکرد اسانس دارند، بنابراین باید در مرحله اول ،تولید زیست‌توده بیش‌تر با هدف افزایش عملکرد اقتصادی مدنظر قرار بگیرد (Ramezan & Abbaszadeh, 2016). زرین گیاه برای تولید حداکثر زیست‌توده، پاسخ مثبتی به شرایط تنش ملایم می‌دهد و همان‌گونه‌که مشاهده می‌شود، در شرایط تنش کم‌آبی ملایم، عملکرد بالاتری نسبت به شرایط آبیاری مطلوب داشت و نسبت به شرایط مطلوب، 12 درصد (210 کیلوگرم در هکتار) وزن تر کل بیش‌تری تولید کرد (جدول 4). افزایش تولید زیست‌توده در شرایط تنش ملایم و کاربرد کود شیمیایی اوره و نیتروکسین را می‌توان به در دسترس بودن اوره در شرایط رشدی گیاه و ایجاد رابطه همزیستی نیتروکسین و در ادامه، تثبیت نیتروژن به وسیله آن و در اختیار قرار دادن نیتروژن در دسترس گیاه نسبت داد. جدول همبستگی صفات مورد آزمایش نشان داد که وزن تر ساقه با وزن تر برگ (r=0.63**) و وزن تر کل (r=0.82**) و همچنین وزن تر برگ با وزن تر کل (r=0.80**) همبستگی مثبت و معنی‌داری را داشت (جدول 5).

 ماده خشک قابل هضم

بر اساس جدول تجزیه واریانس، برهمکنش رژیم‌های کودی و تنش کم‌آبی، تأثیر معنی‌داری در سطح احتمال پنج درصد روی ماده خشک قابل هضم داشتند (جدول 3). جدول مقایسه میانگین اثر متقایل رژیم کودی با تنش کم‌آبی نشان داد که بیشترین درصد ماده خشک قابل هضم، در تیمار آزوکمپوست در شرایط تنش کم‌آبی متوسط با 8/57 درصد تولید شد (جدول 4). افزایش درصد ماده خشک قابل هضم به نسبت محتویات داخل سلول مانند کربوهیدرات‌ها و پروتئین‌های محلول که قابلیت هضم بالایی دارند، به دیواره سلولی بستگی دارد و عوامل محیطی مانند تنش، دما، خاک و ... بر قابلیت هضم تأثیر دارند (Mehdinejad et al., 2017). بنابراین افزایش درصد ماده خشک قابل هضم را می‌توان به جذب عناصر غذایی در صورت تغذیه یا آزوکمپوست و در نتیجه آن افزایش کربوهیدارت‌های محلول و افزایش پروتئین در شرایط تنش کم‌آبی متوسط نسبت داد که با نتایج Fateh et al. (2008) مطابفت دارد. درصد ماده خشک قابل هضم، با الیاف نامحلول در شوینده اسیدی (r=-0.95**) و فیبر خام (r=-0.59**)، همبستگی معنی‌دار و منفی داشت، درحالی‌که با الیاف نامحلول در شوینده خنثی (r=0.35**)، همبستگی معنی‌دار و مثبت نشان داد (جدول 5) که با نتایج Ward et al. (2001) مطابقت دارد.

جدول 3- تجزیه واریانس برخی از صفات زرین گیاه، تحت تأثیر رژیم‌های کودی و تنش کم‌آبی

Table 3. Variance analysis of some Dracocephalum kotschyi Boiss traits affected by fertilizer and water deficit stress.

ns ، ٭ و ٭٭: به ترتیب نشانگر عدم معنی داری و معنی­دار بودن در سطح احتمال پنج و یک درصد می­باشد.

ns, * and **: non significant and significant  at 5%  and 1% of probability levels, respectively.

خاکستر کل

برهمکنش رژیم‌های تغذیه‌ای و تنش کم‌آبی، تأثیر معنی‌داری روی خاکستر کل در سطح احتمال پنج درصد داشت (جدول 3). جدول مقایسه میانگین برهمکنش رژیم کودی و تنش کم‌آبی نشان داد که درصد خاکستر کل در شرایط تنش کم‌آبی شدید با کاربرد اوره و عدم استفاده از کود با 10 درصد خاکستر، بیبیشتر از سایر تیمارها بود و با تیمارهای کاربرد آزوکمپوست و بدون استفاده از کود در شرایط آبیاری مطلوب اختلاف معنی‌داری نداشت (جدول 4). افزایش خاکستر در شرایط آبیاری مطلوب و استفاده از سیستم کودی آزوکمپوست می‌تواند به‌دلیل استفاده بهتر از منابع محیطی و تغذیه‌‌ای باشد که با افزایش جذب عناصر غذایی، درصد خاکستر  افزایش می‌یابد. جدول همبستگی صفات مورد بررسی نشان داد که خاکستر کل با صفت الیاف نامحلول در شوینده خنثی، همبستگی مثبت و معنی‌داری (r=0.39**) داشت (جدول 5). با توجه به این‌که درصد خاکستر، در واقع بیانگر مقدار مواد معدنی موجود در بافت‌های گیاهی است، بنابراین در شرایط تنش کم‌آبی، جذب این مواد توسط ریشه کاهش می‌یابد.

فیبر خام

جدول تجزیه واریانس نشان داد که رژیم‌‌های کودی و تنش کم‌آبی توانستند فیبر خام را در سطح احتمال پنج درصد تحت تأثیر قرار دهند (جدول 3). بیشترین درصد فیبر خام، از تیمار کاربرد ورمی‌کمپوست در شرایط تنش کم‌آبی ملایم با 5/31 درصد به‌دست آمد (جدول 4). جدول همبستگی صفات مورد بررسی نشان دادند که فیبر خام با پروتئین خام (r=-0.65**)، ماده خشک قابل هضم (r=-0.59**) و الیاف نامحلول در شوینده اسیدی (r=-0.46**)، همبستگی منفی و معنی‌داری داشت، اما با الیاف نامحلول در شوینده خنثی (r=0.67**)، همبستگی مثبت معنی‌داری نشان داد (جدول 5). از آن‌جا که فیبر خام شامل مواد غیرقابل هضم علوفه مانند لیگنین، سلولز و همی‌سلولز است، بالابودن آن سبب کاهش محتوی انرژی علوفه می‌شود (Saha et al., 2010). بالابودن درصد فیبر خام در گیاهان در شرایط تنش، به دلیل کاهش ساخته شدن دیواره سلولی است (Kramer, 1983).

جدول 4- مقایسه میانگین برخی از صفات زرین گیاه تحت تأثیر برهمکنش رژیم‌‌های کودی و تنش کم‌آبی

Table 4. Mean comparison of some Dracocephalum kotschyi Boiss traits affected by interaction of fertilizer regimes and water deficit stress

میانگین­های دارای حروف مشترک در هر ستون، تفاوت معنی­داری در سطح احتمال یک درصد نشان ندارند.

Means with the same letters in the same column are not significantly different at the 1 % of probability level.

پروتئین خام

اثر متقایل رژیم‌های کودی و تنش کم‌آبی و اثر اصلی رژیم کودی نتوانستند پروتئین خام زرین گیاه را تحت تأثیر قرار دهند (جدول 3)، درحالی‌که اثر اصلی تنش کم‌آبی در سطح احتمال یک درصد بر پروتئین خام معنی‌دار بود (جدول 3). بیشترین درصد پروتئین خام در شرایط تنش کم‌آبی شدید با 72/19 درصد به‌دست آمد (شکل 1) و تنش کم‌آبی متوسط با 04/18 درصد در رتبه دوم قرار گرفت (شکل 1). رگرسیون خطی سطوح مختلف تنش کم‌آبی نشان داد که با کاهش میزان آب مصرفی (کاهش20 درصدی آب قابل استفاده گیاه در هر سطح تنش)، 27/1 درصد به درصد پروتئین خام گیاه زرین گیاه اضافه شد (شکل 1). پروتئین خام با فیبر خام همبستگی منفی و معنی‌دار (r=-0.65**) و با الیاف نامحلول در شوینده خنثی همبستگی مثبت و معنی‌داری (r=0.62**) داشت (جدول 5). محققان نشان دادند که درصد پروتئین در شرایط دیم و تنش کم‌آبی نسبت به شرایط آبیاری مطلوب افزایش پیدا می‌کند (Weichenthal et al., 2004; Ortega-Ochoa, 2005). زمانی که تنش کم‌آبی اتفاق می‌افتد، شدت تنفس افزایش می‌یابد و جذب مواد کاهش پیدا می‌کند، در نتیجه هیدرات‌های کرین ذخیره شده کم می‌شوند و پروتئین افزایش می‌یابد (Normohamadi et al., 2000). این نتایج با پژوهش‌های Nielsen et al. (2006) و Berenger & Faci (2000) مطابقت دارد.

شکل 1- تأثیر تنش کم‌آبی روی پروتئین خام (حروف مشترک نشان از عدم اختلاف معنی‌داری در سطح احتمال یک درصد می‌باشد)

Figure 1. Effect of water deficit stress on crude protein (similar letters indicate non significant differences at 1% of probability level).

الیاف نامحلول در شوینده اسیدی

الیاف نامحلول در شوینده اسیدی (ADF)، تحت تأثیر برهمکنش رژیم کودی و تنش کم‌آبی در سطح احتمال پنج درصد قرار گرفت (جدول 3). مقایسه میانگین برهمکنش رژیم‌های کودی و تنش کم‌آبی نشان داد که بیشترین درصد الیاف نامحلول در شوینده اسیدی، از کاربرد نیتروکسین در شرایط تنش کم‌آبی متوسط (37 درصد) به‌دست آمد، اما با کاربرد نیتروکسین، ورمی‌کمپوست و آزوکمپوست در شرایط تنش کم‌آبی ملایم و کاربرد اوره، ورمی‌کمپوست و آزوکمپوست در شرایط تنش کم‌آبی شدید، از لحاظ آماری تفاوت معنی‌داری نداشت (شکل 2). جدول ضرایب همبستگی صفات نشان داد که الیاف نامحلول در شوینده اسیدی، همبستگی معنی‌دار و منفی با صفت ماده خشک قابل هضم (r=-0.95**) و همبستگی مثبت و معنی‌داری با صفات قندهای محلول در آب (r=0.28**) و فیبر خام (r=0.67**) داشت (جدول 5).

قندهای محلول در آب

اثر برهمکنش رژیم کودی و تنش کم‌آبی بر قندهای محلول در آب در سطح احتمال یک درصد تحت معنی‌دار بود (جدول 3). مقایسه میانگین برهمکنش رژیم کودی و تنش کم‌آبی نشان داد که بیشترین درصد قندهای محلول در آب، از عدم استفاده از کود در شرایط تنش کم‌آبی متوسط (06/22 درصد) به‌دست آمد (شکل 3)، درحالی‌که با تیمارهای استفاده از کود شیمیایی اوره و کود زیستی نیتروکسین در شرایط تنش کم‌آبی متوسط و کاربرد کود زیستی نیتروکسین در شرایط تنش کم‌آبی ملایم، از لحاظ آماری تفاوت معنی‌داری نداشت (شکل 3).

شکل 2- تأثیر برهمکنش رژیم‌های تغذیه‌ای با تنش کم‌آبی روی الیاف نامحلول در شوینده اسیدی (حروف مشترک نشان‌دهنده عدم اختلاف معنی‌دار در سطح احتمال پنج درصد می‌باشد)

Figure 2. Interaction effect of fertilizers regimes and water deficit stress on acid detergent fiber (similar letters indicate non significant differences at 5% of probability level).

جدول 5- ضرایب همبستگی بین صفات موردبررسی

Table 5. Correlation coefficients between studied traits.

*  و **: معنی‌دار در سطوح پنج و یک درصد.

* and **: significant correlation in 5% and !5 of probability levels.

در بررسی اثر اصلی سطوح مختلف تنش کم‌آبی مشاهده شد که درصد قندهای محلول در آب، از آبیاری مطلوب تا تنش کم‌آبی متوسط رو به افزایش بود، به‌طوری‌که از 51/15 درصد به 48/20 درصد در شرایط تنش متوسط رسید (شکل 3)، اما با افزایش سطح تنش به تنش کم‌آبی شدید، میزان درصد قندهای محلول در آب کاهش پیدا کرد و به نظر می‌رسد که گیاه در سطوح بالای تنش، توانایی مقابله با شرایط تنش را نداشت. قندهای محلول در آب با الیاف نامحلول در شوینده اسیدی (r=0.28**)، همبستگی مثبت و معنی‌دار و با الیاف نامحلول در شوینده خنثی (r=-0.49**) همبستگی منفی و معنی‌داری داشت (جدول 5). گیاهان برای سازگاری به شرایط تنش کم‌آبی، پاسخ‌های فیزیولوژیکی، بیوشیمیایی و مورفولوژیکی مختلفی می‌دهند (Farooq et al., 2009; Mirzaee et al., 2013). در شرایط تنش، افزایش کربوهیدرات‌های محلول در آب، باعث حفظ تورژسانس سلول‌های برگ، حفاظت غشای سلولی و بازداری از تخریب پروتئین‌ها می‌شود؛ همچنین از طریق تأمین انرژی مورد نیاز گیاه، از مرگ حتمی آن جلوگیری می‌کند (Xue et al., 2008). گزارش شده است که کربوهیدرات‌های محلول، در تنظیم اسمزی و مکانیسم‌ها حفاظتی نقش دارند و قندهای محلول در شرایط تنش خشکی افزایش می‌یابد (Kamali et al., 2012).یکی از سازوکارهای سازگاری که گیاه در برابر تنش کم‌آبی از خود بروز می‌دهد، تنظیم اسمزی می‌باشد که با تجمع مواد محلول، سبب حفظ تورژسانس سلول‌ها در پتانسیل‌های پایین آب می‌باشد. کربوهیدرات‌ها نقش اصلی را در تنظیم اسمزی بازی می‌کنند؛ به همین دلیل در شرایط تنش کم‌آبی در گیاه، مقدار کربوهیدرات‌ها بیشتر شد. پاسخ گیاهان به تنش خشکی، به شدت، مدت و مرحله وقوع تنش بستگی دارد (Barnabas et al., 2008). در سطوح بالاتر تنش ، احتمال تخریب دیواره سلولی بالاتر می‌رود و در نتیجه می‌توان شاهد کاهش کربوهیدرات‌های محلول در آب شد. بنابراین کربوهیدرات‌های محلول در آب در شدت‌های بالاتر تنش، کاهش پیدا کرد که با نتایج Faroog et al. (2009) مطابقت داشت.

شکل 3 برهمکنش رژیم‌های تغذیه‌ای و تنش کم‌آبی روی کربوهیدرات‌های محلول در آب (حروف مشترک، نشان‌دهنده عدم اختلاف معنی‌دار در سطح احتمال یک درصد می‌باشند)

Figure 3. Interaction effect of fertilizer regimes and water deficit stress on water soluble carbohydrate (similar letters indicate non significant differences at 1% of probability level).

الیاف نامحلول در شوینده خنثی

الیاف نامحلول در شوینده خنثی (NDF) که در واقع همان دیواره سلولی می‌باشد، پتانسیل مصرف علوفه توسط دام را نشان می‌دهد و گیاهی که NDF پایین‌تری داشته باشد، برای دام خوشخوراک‌تر است (Contreras et al., 2009). بر اساس جدول تجزیه واریانس، الیاف نامحلول در شوینده خنثی، تحت تأثیر برهمکنش رژیم کودی و تنش کم‌آبی در سطح احتمال یک درصد قرار گرفت (جدول 3). مقایسه میانگین برهمکنش رژیم کودی با تنش کم‌آبی نشان داد که بیشترین درصد الیاف نامحلول در شوینده خنثی، از کاربرد کود اوره در شرایط تنش کم‌آبی شدید به‌دست آمد (شکل 4)، اما با تیمارهای نیتروکسین، ورمی‌کمپوست و عدم کاربرد کود در شرایط تنش مشابه، ورمی‌کمپوست در شرایط تنش کم‌آبی متوسط و کاربرد اوره در شرایط آبیاری مطلوب، تفاوت معنی‌داری نداشت (شکل 4). کم‌ترین درصد الیاف، به تیمارهای ورمی‌کمپوست، آزوکمپوست و نیتروکسین در شرایط تنش کم‌آبی شدید و کابرد آزوکمپوست و عدم استفاده از کود در شرایط آبیاری مطلوب تعلق داشت (شکل 4). به‌نظر می‌رسد که آزوکمپوست، شرایط تنش کم‌آبی شدید را تا حدودی کاهش داده است تا گیاه بتواند در شرایط مواجهه، NDF پایین‌تری تولید کند؛ پایین بودن NDF نشان‌دهنده خوشخوراکی گیاه می‌باشد. همچنین در بررسی سطوح محتلف تنش کم‌آبی مشاهد شد که هر چه سطح تنش کم‌آبی بیشتر بود و گیاه در معرض شرایط کم‌آبی بیشتری قرار گرفت، درصد الیاف نامحلول در شوینده خنثی نیز بالا رفت (شکل 4). در شرایط آبیاری مطلوب، به دلیل آبدار شدن سلول‌ها، تکامل و گسترش سلول‌ها کم‌تر شد و مقدار NDF و ADF کاهش یافت که نتیجه آن، افزایش کیفیت علوفه می‌باشد. الیاف نامحلول در شوینده خنثی، همبستگی مثبت و معنی‌داری با پروتئین خام (r=0.62**)، ماده خشک قابل هضم (r=0.35**) و خاکستر کل (r=0.39**) داشت، ولی با قندهای محلول در آب (r=-0.49**) و فیبر خام (r=-0.46**)، همبستگی منفی و معنی‌داری مشاهده شد (جدول 5).

شکل 4- برهمکنش رژیم‌های تغذیه‌ای و تنش کم‌آبی روی الیاف نامحلول در شوینده خنثی (حروف مشترک، نشان‌دهنده عدم اختلاف معنی‌دار در سطح احتمال یک درصد می‌باشد)

Figure 4. Interaction effect of fertilizer regimes and water deficit stress on neutral detergent fiber (similar letters indicate non significant differences at 1% of probability level).

نتیجه‌گیری کلی

نتایج این پژوهش نشان داد که زرین گیاه، به‌دلیل داشتن نزدیک به 20 درصد پروتئین خام می‌تواند گزینه مناسبی به‌عنوان مکمل در جیره غذایی دام قرار باشد. بیشترین درصد ماده خشک قابل هضم از کاربرد آزوکمپوست در شرایط تنش‌کم‌آبی متوسط (8/57 درصد) به‌دست آمد و درصد قندهای محلول در آب در تیمارهایی با سیستم‌های تغذیه‌ای اوره، نیتروکسین و عدم استفاده از کود در شرایط تنش کم‌آبی متوسط و تیمار نیتروکسین در شرایط تنش کم‌آبی ملایم به بیشترین میزان خود رسید. با توجه به نتایج این پژوهش، زرین گیاه توانست شرایط تنش‌کم‌آبی را تحمل کند. با توجه به این‌که علاوه بر انتخاب علوفه‌ای با کیفیت بالا، باید به تولید اندام هوایی مناسب نیز توجه شود، کاربرد کود شیمیایی اوره در شرایط تنش ملایم  با 32 درصد تولید بیشتر نسبت به شاهد، به‌عنوان تیمار برتر انتخاب شد.

REFERENCES

1. Abdi, A., (2008). Evalution of red plant biodiversity in Markazi province. Iranian Journal of Rangelands and Forests Plant Breeding and Genetic Research, 16(1), 50-74. (In Persian)
2. Albayrak, S., M. Turk, O. Yuksel, & Yilmaz, M. (2011). Forage yield and the quality of perennial legume-grass mixtures under rainfed conditions. Notulae Botanicae Horti Agrobotanici Cluj-Napoca, 39(1), 114-118.
3. Anwar, M., Patra, D. D., Chand, S., Kumar, A., Naqvi, A. A. & Khanuja, S. P. S. (2005). Effect of organic manures and inorganic fertilizer on growth, herb and oil yield, nutrient accumulation, and oil quality of French basil. Communications in Soil Science and Plant Analysis, 36(14), 1737-1746.
4. Barker, A. V. & Bryson, G. M. (2006). Comparisons of compost with low or high nutrient status for the growth of plants in containers. Communications in Soil Science and Plant Analysis, 37, 1303-1319.
5. Barnabas, B., Jager, K. & Feher, A. (2008). The effect of drought and heat stress on reproductive processes in cereals. Plant Cell Environ, 31, 11–38.
6. Behera S. K., & Panda R. K. (2009). Effect of fertilization and irrigation schedule on water and fertilizer solute transport for wheat crop in a sub-humid sub-tropical region. Agriculture, Ecosystems and Environment, 130, 141-155.
7. Behrooz, E., Rahmanian, M., Heidarpour, O. & Shahriari, M. H. (2017). Effect of vermicompost and spent mushroom compost on the nutrient and essential oil composition of basil (Ocimum basilicum). Journal of Essential Oil Bearing Plants, 20(5), 1283-1292.
8. Chen, J. (2006). The combined use of chemical and organic fertilizers and /or biofertilizer for crop growth and soil fertility. International Workshop on Sustained Management of the Soil–Rhizosphere System for Efficient Crop Production and Fertilizer Use, 11, 16-20.
9. Contreras-Govea, F. E., Muck, R. E., Armstrong, K. L. & Albrecht, K. A. (2009). Nutritive value of corn silage in mixture with climbing beans. Animal Feed Science and Technology, 150, 1-8.
10. Farooq, M., Wahid, A., Kobayashi, N., Fujita, D. and Basra, S. M. A. (2009). Plant drought stress: effects, mechanisms and management. Agronomy for Sustainable Development, 29, 185– 212.
11. Fateh, E. (2008). Effects of organic and chemical fertilizers on forage yield and quality of Globe artichoke (Cynara scolymus),D. thesis, University of Tehran.
12. Jafari, A., Connolly, V., Frolich, A. & Walsh, E. K. (2003). A note on estimation of quality in perennial ryegrass by near infrared spectroscopy. Irish journal of agricultural and food research, 42, 293-299.
13. Kabirinegad, S. H., Abtahi, A. & Hodji, M. (2009). Short term effects of compost fertilizer on fertigation and concentration of cadmium and iron in two calcareous soils and corn. Proceedings of the Eleventh Iranian Soil Science Congress, Gorgan, 1245-1247. (In Persian)
14. Kamali, M., Goldani, M., & Farzaneh, A. (2012). The effect of that different irrigation levels on growth parameters and photosynthesis and hydrogen peroxide in Amaranthus tricolor. Journal of Soil and Water (Agricultural Science and Technology), 26, 318-309. (In Persian)
15. Kramer, P. J. (1983). Water relation of plants. Academic press, New York.
16. Krishnamoorthy, R. V. & Vajrabhiah, S. N. (1986). Biological activity of earthworm casts: an assessment of plant growth promoter levels in casts. Proceedings of the Indian Academy of Sciences, 95, 341–
17. Mehdinejad, N., Mohammakhani, M. & Fakheri, B. A. (2017). The effect of foliar application manganese on quantitative and qualitative traits of barley (Hordeum vulgare). Journal of Agronomy and Plant Breeding, 13(4), 13-29. (In Persian)
18. Minson, D. J. (1990). Forage in ruminant nutrition. Academic Press, 734 Pp.
19. Mirzaee, M., Moieni, A. & Ghanati, F. (2013). Effects of drought stress on the lipid peroxidation and antioxidant enzyme activities in two canola (Brassica napus ) cultivar. Journal of Agricultural Science and Technology, 1, 593-602.
20. Mozafarian, V. (2003). Dictionary of Iranian Plant Name. Farhang Moaser, 740Pp. (In Persian)
21. Mozafarian, V. (2015). Identification of medicinal and aromatic plants of Iran. Farhang Moaser, 1444Pp. (In Persian)
22. Nielsen, D. C., Vigil, M. F. & Benjamin, J.G. (2006). Forage yield response to water use for dry land corn, millet and triticale in the central Great Plains. Agronomy Journal, 98, 992-998.
23. Normohamadi, G h., Siyadat, S. A. & Kashani, A. (2000). The cultivation of crops. First volume. Shahid Chamran University Press. 446Pp.
24. Ortega-Ochoa, C. (2005). Effect of levels of irrigation on forage standing crop and quality of WW-B. Dahl (Bothriochloa bladhii) pasture under summer grazing. Ph.D. Thesis.
25. Pabby, A., Prasana, R. & Singh, P. K. (2003). Azolla-anabaena symbiosis- from traditional agriculture to biotechnology. Indian Journal of Biotechnology, 2, 26-37.
26. Ramezan, G. & Abbaszadeh, B. (2016). The effect of drought stress on yield, content and percentage of essential oil of Nepeta pogonosperma Jamzad e Assadi under different plant density. Iranian Journal of Medicinal and Aromatic Plants, 31(6), 1071-1085. (In Persian)
27. Saeid-Nejad, A. H., Rezvani-Moghaddam, P. Khazaeii, H. R. & Nasiri-Mahallati, M. (2012). Effect of organic fertilizers, bio-fertilizers and chemical fertilizers on the dry matter digestibility and crude protein of forage sorghum (Speedfeed cv). Iranian Journal of Field Crops Research, 9(4), 623-630. (In Persian)
28. Saha, U. K., Sonon, L. S., Hancock, D. W., Hill, N. S., Stewart, L., Heusner, G. L. & Kissel, D. E. (2010). Common terms used in animal feeding and nutrition. The University of Georgia, College of Agriculture and Environmental. Sciences, Virginia, USA.
29. Soroori, S., Moghaddam, M. & Hashemi Moghaddam, H. (2014). Effect of N fertilizer of different levels and cultivated density on yield and essential oil content on Moldavian dragonhead. Journal of Crops Improvment, 15(4), 179-194. (In Persian)
30. Srinivas, P., Vijaya-Padma, S. S., Pandu-Sastry, K. & Sunitha-Devi, K. B. (2017). Analysis the Effect of Fly Ash and Vermicompost Combination on Herb Yield, Oil Content and Oil Composition of Lemongrass (Cymbopogon flexuosus Nees). International Journal of Pure and Applied Bioscience, 5(4), 1710-1717.
31. Tahmasebi, D., Zarghami, R., Azghandi, A. V. & Chaichi, M. (2011). Effects of nanosilver and nitroxin biofertilizer on yield and yield components of potato minitubers. International Journal of Agriculture and Biology, 13, 986–990.
32. Vessey, J. K. (2003). Plant growth promoting rhizobacteria as biofertilizer. Plant Soil, 255(2), 571-586.
33. Ward, J. D., Redfearn, D. D., McCormick, M. E. & Cuomo, G. J. (2001). Chemical composition, ensiling characteristics, and apparent digestibility of summer annual forages in a subtropical double-cropping system with annual Journal of Dairy Science, 84, 177-182.
34. Weichenthal, B. A., Baltensperger, D. D. & Voge, K. P. (2004). Feed values for annual forages in Western Nebraska. http://digitalcommons. unl.edu/animalscinbcr/21.
35. Xue, G., Lynne McIntyre, C., Glassop, D. & Shorter, R. (2008). Use of expression analysis to dissect alterations in carbohydrate metabolism in wheat leaves during drought stress. Plant Molecular Biology, 67, 197–214.
36. Yousefzadeh, S., Modarres-Sanavy, S. A. M., Sefidkon, F., Ghalavand, A., Roshdi, M. & Safaralizadeh, A. (2011). Effect of biofertilizer, azocompost and nitrogen on morphologic traits and essential oil content of Dracocephalum moldavica in two regions of Iran. Iranian Journal of Medicinal and Aromatic Plants, 29(2), 438-459. (In Persian)

[1] Lamiaceae

[2] Endangered

[3] Neutral detergent fiber

[4] Acid Detergent Fiber

[5] Time-Domain Reflectometer

² Near Infrared Spectroscopy

[7] Dry Matter Digestibility

[8] Water Soluble Carbohydrates

[9] Crude Protein

[10] Acid Detergent Fiber

[11] Total ASH

[12] Crude fiber

[13] Neutral detergent fiber

[14] Statistical Analysis System

[15] Least Significant Difference