УДК 633.11«324»:631.526.32:631.82:631.559 DOI 10.24411/0235-2516-2018-10018
ОПТИМИЗАЦИЯ ПИТАНИЯ СОРТОВ ОЗИМОЙ ПШЕНИЦЫ ПУТЕМ ВНЕСЕНИЯ РАСЧЕТНЫХ ДОЗ МИНЕРАЛЬНЫХ УДОБРЕНИЙ НА ПЛАНИРУЕМЫЙ УРОВЕНЬ УРОЖАЙНОСТИ
А.Н. Есаулко, д.с.-х.н., А.Ю. Ожередова, Н.В. Громова
Ставропольский государственный аграрный университет, e-mail: [email protected]
Представлены материалы по изучению влияния расчетных доз минеральных удобрений на планируемый уровень урожайности сортов озимой пшеницы. Результаты исследований (2010-2014 гг.) показали, что изучаемые методики расчета доз удобрений обеспечили планируемый уровень урожайности озимой пшеницы 4,0 т/га. Наибольшая достоверность (99%) была получена при внесении дозы N126P80K72 под планируемый урожай 6,0 т/га по методике расчета В.В. Агеева и А.Н. Есаулко. В 2015-2017 гг. на сортах Васса и Доля была достигнута планируемая урожайность 5,0 и 7,5 т/га. Урожайность 10 т/га была получена только в 2016 г. у сорта озимой пшеницы Доля.
Ключевые слова: озимая пшеница, минеральные удобрения, методика расчета доз удобрений, планируемая урожайность, программирование, чернозем выщелоченный, сорта озимой пшеницы.
OPTIMIZATION OF NUTRITION OF WINTER WHEAT VARIETIES BY INTRODUCING CALCULATING DOSES OF MINERAL FERTILIZERS ON THE PLANNED CROP LEVEL
Dr.Sci. A.N. Esaulko, A.Yu. Ozheredova, N.V. Gromova
Stavropol State Agrarian University, e-mail: [email protected]
Materials on the study of the effect of calculated doses of mineral fertilizers on the planned level of yields of winter wheat varieties are presented. The results of the studies (2010-2014) showed that the fertilizer dose calculation methods studied provided a planned level of winter wheat yield of 4.0 t/ha. The highest reliability (99%) was obtained by applying a dose of N126P80K72 for the planned yield of 6.0 t/ha according to the calculation method of V.V. Ageev and A.N. Esaulko. In 2015-2017 on varieties Vassa and Share, the planned yields were 5.0 and 7.5 t/ha. The yield of 10 tons per hectare was obtained only in 2016 in the variety of winter wheat.
Keywords: winter wheat, mineral fertilizers, method of calculation offertilizer doses, planned yield, programming, leached chernozem, varieties of winter wheat.
Озимая пшеница основная зерновая культура в Ставропольском крае. С 2010 по 2017 г. посевные площади под озимой пшеницей составили 1626,41717,5 тыс. га [1]. Однако валовые сборы зерна значительно колебались и находились в пределах 22,443,4 ц/га. Связано это с резкими колебаниями урожайности озимой пшеницы, которая составляла 34,4 ц/га в 2010 г., 2011 г. - 39,5 ц/га, 2012 г. - 22,4 ц/га, 2013 г. - 31,5 ц/га, 2014 г. - 40,2 ц/га, 2015 г. - 40,3 ц/га, 2016 г. - 42,8 ц/га, 2017 г. - 43,4 ц/га [2, 3].
Проблема увеличения производства зерна и других сельскохозяйственных продуктов решается за счет дальнейшего повышения продуктивности пашни. Этому способствуют направления - программирование, прогнозирование и планирование урожаев [4-7]. В основе его лежит требование удовлетворения потребности растений в жизненно важных ресурсах для формирования заданного урожая. В Ставропольском крае наблюдаются резкие колебания в урожайности озимой пшеницы, несмотря на стабильные посевные площади. Поэтому програм-
мированное выращивание озимой пшеницы предусматривает оптимизацию минерального питания растений в посевах обеспечением в соответствии с наличием тепла, влаги и света, для повышения урожайности и качества [9, 10]. В связи с этим с 2010 по 2017 г. были проведены исследования по изучению оптимизации питания сортов озимой пшеницы на основе внесения расчетных доз минеральных удобрений на планируемый уровень урожайности на сельскохозяйственной опытной станции Ставропольского ГАУ, которая расположена на северном склоне Ставропольской возвышенности в Грачев-ско-Калаусском ландшафте лесов и степей. Территория характеризуется неустойчивым увлажнением и неравномерностью выпадения осадков в течение года. Средняя многолетняя сумма осадков составляет 551 мм, за вегетационный период выпадает 350370 мм, среднегодовая температура воздуха 9,2°С. Гидротермический коэффициент - 1,1-1,3 [11].
Первый опыт был проведен в 2010-2014 гг. Цель исследований заключалась в оптимизации
применения минеральных удобрений на основе изучения балансовых методик расчета доз удобрений для достижения программируемого уровня урожайности озимой пшеницы. Объектом исследований был сорт озимой пшеницы Зустрич.
Второй опыт проведен в 2015-2017 гг. Цель исследований - изучение внесения расчетных доз минеральных удобрений для достижения максимальной продуктивности сортов озимой пшеницы на черноземе выщелоченном Ставропольской возвышенности. Объектом исследований были сорта озимой пшеницы Васса, Доля и Гром.
Методика исследований. Территория характеризуется неустойчивым увлажнением и неравномерностью выпадения осадков в течение года. Средняя многолетняя сумма осадков составляет 551 мм, за вегетационный период выпадает 350-370 мм, среднегодовая температура воздуха 9,2°С. Гидротермический коэффициент - 1,1-1,3 [11]. Почва -чернозем выщелоченный, мощный, малогумусный тяжелосуглинистый, который характеризуется в настоящее время средним содержанием гумуса (5,1-5,4%), подвижного фосфора (20-25 мг/кг по Мачигину) и обменного калия (220-270 мг/кг). Нитрификационная способность 16-30 мг/кг. Реакция почвенного раствора в верхних горизонтах почвы нейтральная, находится в пределах 6,1-6,5.
Опыт № 1. Размещение делянок по методу рен-домизированных повторений, повторность опыта трехкратная. Общая площадь опыта - 960 м2 (ширина - 12 м, длина - 80 м), площадь одной делянки
- 35 м2. Опыт двухфакторный: фактор А - планируемая урожайность озимой пшеницы - сорт Зустрич 4,0, 5,0 и 6,0 т/га; фактор В - способы расчета минеральных удобрений.
Схема опыта представлена следующим образом: 1. Контроль - без удобрений; 2. Рекомендованная -N60P60K30 [12]; 3. Планируемый урожай 4 т/га по методике В.В. Агеева и А.Н. Есаулко [13] -N60P34K34; 4. Планируемый урожай 4 т/га по методике Ставропольского НИИСХ и ГЦАС «Ставропольский» [14] - N68P44K24; 5. Планируемый урожай
5 т/га по методике В.В. Агеева и А.Н. Есаулко [13]
- N105P60K60; 6. Планируемый урожай 5 т/га по методике Ставропольского НИИСХ и ГЦАС «Ставропольский» [14] - N90P67K40; 7. Планируемый урожай 6 т/га по методике В.В. Агеева и А.Н. Есаулко [13] - N126P80K72; 8. Планируемый урожай
6 т/га по методике Ставропольского НИИСХ и ГЦАС «Ставропольский» [14] - N110P82K51.
Расчет доз минеральных удобрений на планируемую урожайность проводили по двум методикам. В соответствии с первым подходом, разработанным В.В. Агеевым и А.Н. Есаулко (2006), дозы фосфорных и калийных удобрений рассчитывали по формуле:
(В - В ■ Kn) Д Ку ■ 100 ,
где: Д - доза Р2О5 и К2О, кг/га; В - вынос Р2О5 и К2О с планируемым урожаем зерна, кг/га; Кп - коэффициент использования фосфора и калия из почвы от выноса с планируемым урожаем зерна с учетом содержания в почвах подвижных форм фосфора и калия и планируемой урожайности [13]; Ку -коэффициент использования фосфора и калия из удобрений (40 и 70% соответственно).
Дозы азотных удобрений рассчитывали по преобразованной формуле:
_ (В(К) - В(Р2О5) ■ Кп (Р2О5) ■ К) Д Ку ■ 100 '
где: К - отношение выноса N с планируемым урожаем зерна к выносу Р2О5 с планируемым урожаем зерна; Ку - коэффициент использования азота из удобрений (70%).
Согласно второй методике, разработанной специалистами Ставропольского НИИСХ и ГЦАС «Ставропольский» [14], дозы удобрений были рассчитаны по формуле:
Д = У ■ В ■ Кк, где: У - планируемая урожайность зерна, т/га; В -вынос N Р2О5 и К2О с 1 ц планируемого урожая зерна, кг; Кк - коэффициент компенсации выноса элементов питания за счет удобрений.
Опыт № 2. Делянки размещали по методу рен-домизированных повторений, повторность опыта трехкратная. Ширина - 3,6 м, длинна 5 м, общая Б делянки - 18 м2.
Опыт двухфакторный: фактор А - сорта озимой пшеницы Доля, Васса, Гром; фактор Б - дозы минеральных удобрений на продуктивность 5,0, 7,5 и 10,0 т/га. Схема опыта: 1. контроль - КбзР52 (фон); 2. планируемый урожай 5,0 т/га - КшР72Кзо; 3. планируемый урожай 7,5 т/га - КшР95К45; 4. планируемый урожай 10,0 т/га - К248РшКбо.
Расчет доз минеральных удобрений на планируемую урожайность озимой пшеницы 5,0, 7,5 и 10,0 т/га проводили по методике В.В. Агеева (2004). Предшественник - горох.
Результаты исследований. Опыт № 1. Наиболее благоприятные погодные условия для формирования урожая сложились в 2010-2011 г. Сумма осадков, выпавших за вегетацию культуры (580 мм), уступала норме 7%, однако их распределение способствовало оптимальной влагообеспеченности посевов и формированию наибольшей урожайности озимой пшеницы. Среднегодовая температура воздуха оказалась на 1,4°С выше многолетних значений, достигнув 10,6°С. Погодные условия в 20112012 г., сложились крайне неблагоприятно для формирования урожая, количество осадков оказалось меньше среднемноголетней годовой нормы на 27%. Погодные условия в 2012-2013 г. оказалось еще хуже, количество осадков было меньше многолетней годовой нормы на 15%. В 2013-2014 гг. выпало 705 мм осадков, что превысило среднемного-летние показатели на 13%. Однако их распределение
1. Расчет доз удобрений на планируемый уровень урожайности озимой пшеницы
Дозы удобрений Методика расчета Планируемая урожайность Урожайность, т/га
2011 г. 2012 г. 2013 г. 2014 г. средняя
0 Контроль - 3,12 2,63 3,74 3,25 3,19
Рекомендованная - 4,3 3,60 4,90 3,91 4,18
1 4,0 4,15 3,72 4,43 4,15 4,11
2 4,39 3,93 4,85 4,10 4,32
Nl05P60K60 1 5,0 4,63 4,34 5,57 4,90 4,86
N90P67K40 2 5,17 4,21 5,42 4,62 4,90
№26В80^2 1 6,0 6,02 4,91 6,91 5,85 5,92
N1^82^1 2 5,8 4,61 6,23 5,94 5,65
НСР05 - - 0,27 0,32 0,45 0,30 0,34
Sх, % - - 3,6 4,6 3,3 4,3 4,0
Примечание: 1. Методика расчета по В.В. Агееву и А.Н. Есаулко (2006); 2. Методика расчета ученых Ставропольского НИИСХ и ГЦАС «Ставропольский» (1987).
по вегетации растений было неравномерным, а наблюдавшийся дефицит осадков отмечался на фоне высоких температур воздуха и атмосферной засухи. При этом повышенный температурный режим был отмечен на протяжении всего вегетационного периода, среднегодовая температура превысила среднемноголетние показатели на 0,3°С.
Из данных таблицы 1 видно, что все изучаемые дозы минеральных удобрений достоверно увеличивали урожайность озимой пшеницы, и разница относительно контроля в 2011 г. составляла 1,03-2,9 т/га, в 2012 г. - 0,97-2,28 т/га, в 2013 г. - 0,69-3,17 т/га, в 2014 г. - 0,66-2,69 т/га.
При расчете доз минеральных удобрений на планируемый уровень урожайности озимой пшеницы 4,0 т/га установлено, что все изучаемые методы показали довольно высокую точность программирования урожайности, а отклонения от +3 до +7% был несущественными. Более высокий (4,85 т/га) уровень урожайности отмечен при внесении N68P44K24, рассчитанной по методике Ставропольского НИИСХ и ГЦАС «Ставропольский» (1987) урожайность оказалась выше на 18% по сравнению с расчетом доз минеральных удобрений по методике В.В. Агеева и А.Н. Есаулко (2006).
При планировании урожайности 5,0 т/га наблюдались аналогичные результаты. Обе расчетные методики обеспечили несущественное отклонение от планируемого уровня урожайности в сторону уменьшения -2,0 и -2,5%. При внесении удобрений на планируемый уровень урожайность 6,0 т/га точность методики В.В. Агеева и А.Н. Есаулко (2006) оказалась на 4,8% выше методики расчета, рекомендованной учеными Ставропольского НИИСХ и ГЦАС «Ставропольский». При внесении Nl26P8oK72 отмечалось существенное отклонение от планируемого уровня урожайности на 11%, чего нельзя сказать при внесении N1^82^1, рассчитанной по методике Ставропольского НИИСХ и ГЦАС «Ставропольский», где отклонение находилось в пределах ошибки опыта.
В среднем за годы исследований не удалось получить программированный уровень урожайности
озимой пшеницы 6,0 т/га независимо от методики расчета, но наибольший эффект был получен в 2012-2013 гг., когда все варианты обеспечивали достоверную прибавку урожайности (отклонение составило +4 и +13%). Сравнение изучаемых методик расчетных доз минеральных удобрений (N60P34K34 и N68P44K24) на программируемый уровень урожайности 4,0 и 6,0 т/га показало, что существенной разницы урожайности озимой пшеницы не выявлено. В среднем за 4 года оба метода расчета доз удобрений обеспечили программированный уровень 4,0 т/га. Программированный уровень 5,0 и 6,0 т/га достигнут не был, но наибольшая достоверность программирования 99% была получена при внесении дозы Nl26P8oK72 под планируемый урожай 6,0 т/га по методике В.В. Агеева и А.Н. Есаулко (2006).
Опыт № 2 был продолжением исследований. Расчет норм удобрений под планируемый урожай 5,0, 7,5, 10,0 т/га проводили по методике В.В. Агеева и А.Н. Есаулко (2006), так как данная методика оказалась более эффективной для исследуемой зоны.
Применение минеральных удобрений предусматривало два способа внесения: допосевное (под основную обработку почвы) и подкормки. Согласно схеме опыта, в среднем за два года дозы удобрений составляли: до посева - N^52 (рекомендованная), ^4Р72К30 (на планируемую урожайность 5,0 т/га), №бР95^5 (на планируемую урожайность 7,5 т/га), КшРшК60 (на планируемую урожайность 10,0 т/га). Подкормки проводили в фазу кущения (на контроле - Naaзo, при планировании урожайности 5,0 т/га - Naaзo, 7,5 т/га - Naa50 и 10,0 т/га -Naa5o), выхода в трубку (на контроле - Naa20, при планировании урожайности 5,0 т/га - Naa20, 7,5 т/га - Naaзo и 10,0 т/га - Naa5o) и колошения (при планировании урожайности 5,0 т/га - Nм20, 7,5 т/га -Ш20 и 10,0 т/га - Ш20). В качестве удобрений в опыте были использованы аммофос, калий хлористый, аммиачная селитра, мочевина.
Наиболее благоприятные погодные условия для формирования урожая озимой пшеницы сложились
в 2015-2016 гг. Сумма осадков, выпавших за вегетацию (654 мм), превышала норму на 19%, однако их распределение способствовало оптимальной влагообеспеченности посевов и формированию наибольшей урожайности. Среднегодовая температура воздуха оказалась на 0,9°С выше многолетних значений, достигнув 10,1°С.
Погодные условия в 2016-2017 гг. из-за неравномерного распределения осадков в весенне-летний период оказали неблагоприятное влияние на формирование урожая озимой пшеницы. В целом за 2016-2017 гг. выпало 661 мм осадков, что превысило среднемноголетние показатели на 20%. Повышенный температурный режим был отмечен на протяжении всей вегетации растений, среднегодовая температура превысила многолетние показатели на 0,2°С. Средняя урожайность по всем вариантам в 2015-2016 гг. составила 6,69 т/га, в 20162017 г. - 6,01 т/га (табл. 2).
Все изучаемые дозы минеральных удобрений достоверно увеличивали урожайность сортов озимой пшеницы относительно контроля. Максимальная урожайность была получена при внесении дозы К248РшКбо на планируемую урожайность 10,0 т/га и в среднем за два года составила 8,8 т/га. На кон-
троле максимальная урожайность зерна в среднем за два года сформировалось у сорта Доля 4,39 т/га. Планируемый уровень урожайности 5,0 т/га при внесении дозы КшР72Кзо был достигнут на всех сортах озимой пшеницы, а максимальный показатель так же обеспечил сорт Доля, но разница со сравниваемыми сортами несущественна. Уровень продуктивности на 7,5 т/га был достигнут при внесении К18бР95К45 у сортов Васса (7,58 т/га) и Доля (7,59 т/га). Планируемый уровень урожайности озимой пшеницы 10,0 т/га, при применении расчетной дозы №48РшКбо не был достигнут ни на одном сорте, а наибольшая достоверность 91% была получена на сорте Гром (9,06 т/га).
Изучаемые дозы минеральных удобрений оказали положительное влияние на качество зерна озимой пшеницы клейковины и белка (табл. 3).
По сравнению с контролем содержание клейковины у сорта Васса увеличилось на 1,6-4,0%, Гром -1,7-4,6% и Доля - на 2-4,9%, а увеличение содержание белка составило соответственно 1,0-1,8; 1,1-2,3 и 1,4-2,6%. Наибольшее (26%) содержание клейковины и белка в среднем за 2016-2017 гг. было получено при внесении ^8бР95К45 на варианте с планируемой продуктивностью 7,5 т/га. Применение всех изу-
2. Урожайность сортов озимой пшеницы на основе оптимизации минерального питания
Планируемая урожайность, т/га Дозы удобрений Сорта озимой пшеницы Урожайность, т/га
2016 г. 2017 г. средняя
Контроль N63P52 Васса 3,60 3,18 3,39
Гром 4,31 2,90 3,60
Доля 4,80 3,98 4,39
5,0 N124P72K30 Васса 4,71 5,45 5,08
Гром 5,32 4,66 4,99
Доля 6,13 5,23 5,68
7,5 N186P95K45 Васса 7,51 7,64 7,58
Гром 7,30 7,09 7,20
Доля 8,39 6,79 7,59
10,0 N248P133K60 Васса 8,23 9,23 8,73
Гром 9,46 8,65 9,06
Доля 10,47 7,30 8,89
НСР05 0,82 0,74 0,75
Планируемая Дозы Сорта озимой Содержание Показатель ИДК Белок
урожайность, т/га удобрений пшеницы клейковины, %
Васса 21,5 67 11,5
Контроль N63P52 Гром 21,5 67 11,5
Доля 21,4 65 11,4
Васса 23,2 74 12,6
5,0 N124P72K30 Гром 24,0 75 13,2
Доля 23,4 75 12,8
Васса 25,5 84 13,3
7,5 N186P95K45 Гром 26,1 83 13,8
Доля 26,3 85 14,0
Васса 23,1 12,5
10,0 N248P133K60 Гром 23,3 12,6
Доля 23,7 12,8
3. Качество зерна озимой пшеницы, среднее за 2016-2017 гг.
ченных доз минеральных удобрений способствовало и получению клейковины хорошего качества - показания прибора ИДК составили 74-85 ед. В итоге на контроле у всех изучаемых сортов было получено зерно 4 класса, а на остальных вариантах - 3 класса.
Таким образом, в ходе исследований 20102014 гг. было выявлено, что все изучаемые дозы минеральных удобрений существенно увеличивали урожайность озимой пшеницы по сравнению с контролем. Сравнение изучаемых методик расчетных доз минеральных удобрений на программируемый уровень урожайности 4,0 и 6,0 т/га показало, что существенной разницы в показателях урожайности озимой пшеницы не выявлено. В среднем за 4 года оба метода расчета доз удобрений обеспечили программированный уровень урожайности озимой пшеницы 4,0 т/га. Уровень урожайности 5,0 и 6,0 т/га достигнут не был, но наибольшая достоверность программирования 99% была получена при внесении дозы Nl26PsоK72 под планируемый урожай 6,0 т/га по методике расчета В.В. Агеева и А.Н. Есаулко (2006).
Все изучаемые дозы минеральных удобрений в
2016-2017 гг. существенно увеличивали урожайность озимой пшеницы по сравнению с контролем. Планируемый уровень урожайности 5,0 т/га при внесении дозы N12^7^30 был достигнут на всех сортах озимой пшеницы, а максимальный показатель обеспечил сорт Доля, но разница была несущественна. Уровень продуктивности на 7,5 т/га был достигнут при внесении NlS6P95K45 у сортов Васса (7,58 т/га) и Доля (7,59 т/га). Планируемый уровень 10,0 т/га при применении дозы N248Р1ззКбо не был достигнут ни на одном сорте, а наибольшая достоверность 91% была получена на сорте Гром (9,06 т/га).
Изучаемые дозы минеральных удобрений положительно влияли на качество зерна. Наибольшее (26,0%) содержание клейковины и белка в среднем за 2016-2017 гг. было получено при внесении NlS6P9sK4s на варианте с планируемой продуктивностью 7,5 т/га Существенного влияния изучаемые сорта на качественные показатели зерна не оказали, отмечается недостоверное увеличение содержания клейковины и белка у сорта Доля.
Литература
1. Трухачев В.И., Костюкова Е.И., Герасимов А.Н., Скрипниченко Ю.С. Аграрная экономика Ставропольского края: приоритеты и перспективы инновационного развития. - Ставрополь: СтавГАУ, 2016. - 272 с.
2. Сигида М.С. Влияние систем удобрений на продуктивность звена зернопропашного севооборота на выщелоченном черноземе: дисс. к.с.-х.н. - Ставрополь, 2008. - 237 с.
3. Esaulko A.N., Ozheredova A.Yu., Sigida M.S., Voskoboinikov A.V., Podkolzin O.A. Introduction of calculated doses of mineral fertilizers to achieve maximum productivity of winter wheat varieties on chernozem leached Stavropol upland // Research Journal of Pharmaceutical, Biological and Chemical Sciences, 2017, Т. 8, № 6. - Р. 778-781.
4. Минеев В.Г. Значение географической сети опытов с удобрениями в решении актуальных проблем агрохимии / Совершенствование организации и методологии агрохимических исследований в Географической сети опытов с удобрениями: материалы Всероссийской научно-методической конференции. - М.: ВНИИА, 2006. - С. 6-8.
5. Беличенко М.В., Рухович О.В., Романенков В.А. Использование результатов длительных полевых опытов с удобрениями для разработки стратегии обеспечения стабильных урожаев / 75 лет Географической сети опытов с удобрениями: сборник материалов Всероссийского совещания научных учреждений - участников Географической сети опытов с удобрениями. - М.: ВНИИА, 2016. - С. 23-27.
6. Бузов В.А., Гречишкина Ю.И. Продуктивность озимой пшеницы на черноземе при подкормке различными формами азотных удобрений // Плодородие, 2010, № 1. - С. 16-17.
7. Гречишкина Ю.И., Бузов В.А., Подорогин В.А. Эффективность форм азотных удобрений, при ранневесенней подкормке озимой пшеницы на черноземе выщелоченном // Плодородие, 2011, № 4. - С. 22-23.
8. Милащенко Н.З., Завалин А.А., Сычев В.Г., Самойлов Л.Н., Трушкин С.В. Факторы повышения эффективности удобрений в интенсивных технологиях возделывания пшеницы в России // Агрохимия, 2015, № 11. - С. 13-18.
9. Дорожко Г.Р., Власова О.И., Цховребов В.С. Развитие земледелия Ставрополья / Эволюция и деградация почвенного покрова: сборник научных статей по материалам V Международной научной конференции (г. Ставрополь, 19-22 сентября 2017 г.). - Ставрополь: ООО «СЕКВОЙЯ», 2017. - С. 249-251.
10. Трухачев В.И., Громов Е.И. Оценка уровня развития сельских территорий в разрезе регионов России // Экономика сельского хозяйства России, 2016, № 4. - С. 57-65.
11. Бадахова Г.Х., Кнутас А.В. Ставропольский край: современные климатические условия. - Ставрополь: Краевые сети связи, 2007. - 272 с.
12. Агеев В.В., Подколзин А.И. Системы удобрения в севооборотах Юга России: учебное пособие для студентов вузов агрономических специальностей (2 изд., перераб. и доп.). - Ставрополь: Ставропольский ГАУ, 2001. - 349 с.
13. Есаулко А.Н. Оптимизация систем удобрений в севооборотах центрального Предкавказья как фактор повышения плодородия почвы и продуктивности сельскохозяйственных культур: дисс. д.с.-х.н. - Ставрополь, 2006. - 515 с.
14. Петрова Л.Н., Чернов А.Я., Шустикова Е.П., Подколзин А.И., Карандашов Л.Г., Булавинов А.В. Методические указания для расчета потребности и распределения фондов минеральных удобрений в колхозах и совхозах Ставропольского края. - Ставрополь: Агрус, 1987. - 20 с.