8. Bolotov A.G. Izmerenie temperatury pochv v polevykh usloviyakh // Sb. materialov II mezhd. nauch.-prak. konf. Antropogennoe vozdeystvie na lesnye ekosistemy. - Barnaul, 2002. - S. 148-150.
9. Bolotov A.G., Shein E.V., Makarychev S.V. Water Retention Capacity of Soils in the Altai Region. Eurasian Soil Science, 52 (2): 187-192 (2019). doi: 10.1134/S1064229319020030.
10. Bolotov G.A., Bekhovykh Yu.V., Seme-nov G.A. Opredelenie teplofizicheskikh svoystv kapil-lyarno-poristykh tel impulsnym metodom s ispol-zovaniem tekhnologii vizualnogo programmirovaniya // Vestnik Altayskogo gosudarstvennogo agrarnogo universiteta. - 2012. - No. 6 (68). - S. 37-40.
+ + +
I
УДК 631.445.4:504.53.06
О.В. Коваленко, Ю.Н. Сыромятников
O.V. Kovalenko, Yu.N. Syromyatnikov
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ТРАНСЛОКАЦИИ МИКРОЭЛЕМЕНТОВ-МЕТАБОЛИТОВ В СИСТЕМЕ ПОЧВА-РАСТЕНИЕ В УСЛОВИЯХ ЕЕ ПОЛИЭЛЕМЕНТНОГО ЗАГРЯЗНЕНИЯ ТЯЖЕЛЫМИ МЕТАЛЛАМИ
THE MATHEMATICAL MODELING OF TRACE ELEMENT METABOLITE TRANSLOCATION IN THE SOIL - PLANT SYSTEM UNDER THE CONDITIONS OF ITS MULTIPLE-ELEMENT POLLUTION WITH HEAVY METALS
Ключевые слова: микроэлементы, почва, растения, загрязнение почвы.
Доказана возможность использования метода математического моделирования как одного из подходов для изучения закономерностей транслокации микроэлементов-метаболитов ^п, Fе, Со, Мп, Си) в системе почва-растение при различных уровнях ее техногенной нагрузки Сd, РЬ, N Сг в условиях микрополевого опыта на черноземе подзолистом тяжелосуглинистом. Сформулированы основные принципы построения многофакторных регрессионных моделей и дана оценка транслокации микроэлементов под действием тяжелых металлов в системе почва-растение.
Keywords: trace elements, soil, plants, soil pollution.
This paper proves the possibility of using the method of mathematical modeling as one of the approaches to study the patterns of translocation of trace element metabolites (Zn, Fe, Co, Mn, Cu) in the soil - plant system at various levels of its technogenic load with Cd, Pb, Ni, Cr in a micro-field experiment in podzolic and heavily loamy chernozem soil. The basic principles of the construction of multivariate regression models are formulated and the translocation of trace elements under the action of heavy metals in the soil - plant system is evaluated.
Коваленко Ольга Владимировна, инженер, Харьковский национальный университет им. В.Н. Каразина, г. Харьков, Украина. E-mail: 7109905@gmail.com. Сыромятников Юрий Николаевич, аспирант, Харьковский национальный технический университет сельского хозяйства им. Петра Василенко, г. Харьков, Украина. E-mail: gara176@meta.ua.
Введение
Проблема полиэлементного загрязнения почв тяжелыми металлами и получения качественной продукции в зонах техногенного давления является актуальной [1-6]. Для эффективного анализа и выявления закономерностей накопления металлов-токсикантов, микроэлементов-метаболитов в трофических цепях необходимо, на наш взгляд, выявить взаимосвязь между факторами, определяющими направление процессов, и представить их в количественной форме - в виде математической модели.
Одной из задач исследований вопросов техногенного загрязнения системы почва-растение является разработка модели суммарного воздействия техногенного полиэлементного загрязнения почв тяжелыми металлами на процесс транслокации микроэлементов-метаболитов. Цель моделирования в исследованиях - нахождение математической формулы для имеющихся составляющих и их взаимной связи в системе почва-растение и дальнейшее их модельное и теоретическое объяснение.
Статистический уровень наших исследований имеет в основном эмпирический характер, так как разработки детальной гипотетической модели, которая объяснила бы причины вариации показателей перед началом исследования, не было. Кроме того, необходимость включения в совокупности условий большого количества показателей привела к значительным трудностям по сравнению с результатами статистического анализа, который показал, что данные содержания металлов в системе почва-растение можно описать уравнениями множественной регрессии. Это позволило выявить элементы, содержание которых выходит за пределы дефицита и избытка при полиэлементном загрязнении почвы тяжелыми металлами.
Kovalenko Olga Vladimirovna, Engineer, Kharkov National University named after V.N. Karazin, Kharkov, Ukraine. E-mail: 7109905@gmail.com. Syromyatnikov Yuriy Nikolayevich, post-graduate student, Kharkov National Technical University of Agriculture named after P.M. Vasilenko, Kharkov, Ukraine. E-mail: gara176@meta.ua.
v-
Методика и объекты
Зависимости транслокации микроэлементов-метаболитов в системе почва-растение изучались на групповых объектах методами математической статистики (дисперсионный, корреляционный, регрессивный и факторный анализы) с использованием пакета программ А.И. Новикова [7], усовершенствованных Ю.Н. Сыромятниковым [8]. Объектом исследований избран чернозем подзолистый тяжелосуглинистый Восточно-Лесосте-повой возвышенной провинции Украины с смоделированными уровнями техногенного загрязнения, превышающими естественный фон по содержанию Cd, Pb, Ni, Сг в 3, 5, 10, 15 раз. Как тестовые культуры использованы суданская трава (Soeghum vulgare sudanence), гречка (Fagoppyrum sagittayum), горох (Pisum sativum var. commune). Определение тяжелых металлов и микроэлементов-метаболитов в анализируемой системе проводилось с использованием атомно-абсорб-ционного метода.
Результаты и обсуждение
Использование статистических многомерных методов анализа позволило установить, что содержание цинка в растениях коррелирует с содержанием в них Cu (r = 0,79), Mn (r = 0,64), Cd (r = 0,58) и Cr (r = 0,5), а также с содержанием Cu (r = 0,62), Cr (r = 0,58), Mn (r = 0,5), Fe (r = 0,5) в почве. Лабильная форма цинка в почве связана с содержанием кислоторастворимых форм Ni (r = 0,54), Fe (r = 0,5).
Транслокация железа в растения обусловлена преимущественно фазой вегетации растений (r = 0,9), легкодоступными формами Cu (r = 0,69), Cr (r = 0,59), Pb (r = 0,5) и кислоторастворимыми формами Mn (r = 0,6), Zn (r = 0,53), Cr (r = 0,5), а также содержанием Cr (r = 0,92) в растениях.
Поступление кобальта в растения, изученное ранее, тесно связано с содержанием в них Fe (г = 0,97) и Сг (г = 0,9), фазой их вегетации (г = 0,98) и содержанием легкодоступных форм Со (г = 0,85), Мп (г = 0,73) в почве.
Накопление марганца тестовыми культурами связано с содержанием в них Zn (г = 0,76), Сг (г = 0,62), Си (г = 0,57), Cd (г = 0,52), фазой вегетации (г = 0,74), с уровнями поступления в почву Сг (г = 0,5), а также с содержанием в почве легкодоступных форм Мп (г = 0,62), Сг (г = 0,62) и кис-лоторастворимой формы Сг (г = 0,52). Количество марганца в почве коррелирует с содержанием в нем легкодоступных форм Сг (г = 0,99), Си (г = 0,72), Со (г = 0,69), N (г = 0,67), Zn (г = 0,65) и кислоторастворимых форм Fe (г = 0,65), Со (г = 0,5).
Поступление меди в растения преимущественно обусловлено фазой их вегетации (г = 0,83), содержанием в них Со (г = 0,82), Zn (г = 0,82), Мп (г = 0,5), а также содержанием в почве кислоторастворимых форм Си (г = 0,51), Сг (г = 0,43) и легкодоступных форм Си (г = 0,73), Сг (г = 0,48). Легкодоступные формы меди в почве связаны с кислоторастворимыми формами Fe (г = 0,66), Мп (г = 0,64), Со (г = 0,51), Си (г = 0,5) и легкодоступными формами Сг (г = 0,66), Со (г = 51), Си (г = 0,5), Мп и РЬ (г = 0,5) в почве.
По данным регрессионного анализа установлено, что урожайность растений зависит в большинстве случаев от содержания в них микроэлементов-метаболитов, то есть элементов, жизненно необходимых для нормального их развития: Мп (г = 0,71), Fe (г = 0,7), Со (г = 0,58), Си (г = 0,5), от уровня поступления в почву фоновых количеств загрязнителей: Сг (г = 0,59), Cd (г = 0,5) и их содержания в растениях (г = 0,55; г = 0,66). Особенно тесная связь выявлена с наличием в почве таких легкодоступных растениям форм металлов, как Сг (г = 0,62), Мп (г = 0,59), N (г = 0,58) и несколько меньше Cd (г = 0,5), а также кислоторастворимых форм N (г = 0,62), Сг (г = 0,54).
Поскольку в публикациях по изучению транслокации микроэлементов нет ссылок на выбор адекватной естественному процессу модели, эм-
пирическим путем установлено, что целесообразнее использовать многофакторные регрессионные модели [9, 10]. Результатом математической обработки полученных данных стала десятифак-торная полная квадратичная регрессионная модель [11], которая устанавливает формы зависимости, функции регрессии и дает оценку неизвестным значениям зависимой переменной в системе почва-растение при моделировании уровней загрязнения почвы тяжелыми металлами [12-15].
Полная квадратичная модель содержит в себе 65 параметров и 10 факторов РЬ, N Сг, Мд, д х2, Рр, Ъ Ки) и имеет такой общий вид:
У = a0 +Х aixi +Х ai 2 Х12 +Х
Q.-3Х-X ■.
I 3 l J
Для рассматриваемых в системе почва-растение микроэлементов-метаболитов и тяжелых металлов получены следующие уравнения (1)-(5):
ZnR = 101,41 + 2,68Znbuf + 9,55Pp - 83,3t + 12,79t2 + + 2,32Znbuft - l,71ZnbufKu - 0,203CdZnbuf -
- 0,785CrPp - 4,66ZnbufPp + 2,08Znbuft + 0,944Crt +
+ 3,89PpKu - 0,014СгМд; R = 0,804; F (13,188) = 26,57; Fst = 1,76; (1) Cor = 7,68 - 0,82Pp2 + 0,24tKu - 0,24Ku2 + 3,68Pp -- 0,72СонаРр - 1,83COhcl - 6,15t + + 0,77t2 + 1,3CoHCLt; R = 0,794; F (9,192) = 36,45; Fst = l,93; (2) FeR = 3066,39 - 2489,99Pp + + 485,32Pp2 + 6,53FebufKu; R = 0,91; F (3,198) = 318,45; Fst = 2,65; (3) MnR = -314,6 - 15,08Ppt - 18,22PpKu -
- 0,59MnHCLKu + 355,941Ku - 40,52Ku2 - 21,18tKu +
+ 0,27MnHCLt + 151,84Pp - 27,15Pp2; R = 0,708; F(9,192) = 21,55; Fst = l,93; (4) Cur = 9,62 - 5,34t + 1,12t2; R = 0,474; F (2,199) = 28,86; Fst = 3,04; (5) где Cd, Pb, Ni, Cr - превышение кларков исследуемых металлов по схеме опыта; t - время отбора образцов; Рр - органы растений; Мд - минеральные удобрения; Х1 - содержание исследуемого элемента в буферной вытяжке из почвы;
Х2 - содержание элемента в кислотораствори-мой вытяжке из почвы;
Ки - исследуемая культура.
Построенная нами модель позволяет выявить многочисленные явления взаимодействия металлов в почве и растениях, которые характеризуются изменчивостью и разнородностью проявления в системе, и описать их ограниченным количеством факторов влияния. Так, нами установлено, что содержание цинка в растениях определяется превышением фонового уровня в почве хрома и кадмия, содержанием подвижной формы цинка в почве, видом выращиваемой культуры, а также периодом ее вегетации. Содержание других исследованных микроэлементов (Со, Ре, Мп) в растениях связано с их содержанием в почве, видом выращиваемой культуры, распределением металлов в органах растений по времени (Си, Мп). Концентрация других металлов-токсикантов (РЬ, N1, С^ Сг) в растениях тесно связана с содержанием их в почве, распределением в органах растений, видом выращиваемой культуры. Прослеживается также связь содержания металлов в почве с внесенным их количеством.
С целью дальнейшей детализации полученной информации построили пятифакторную полную квадратичную регрессионную модель, которая содержит превышение местных кларков С^ РЬ, N1, Сг и фактор времени. Результатом стали такие уравнения (6)-(20), описывающие прочность, тесноту и форму связи между исследуемыми параметрами для суданской травы (г1), гречки (г2) и гороха (г3):
Znrl = 61,92 - 46,61 + 10,6912 + 5,77Сг - 0,73Сг2 +
+ 0,31Си№ - 0,38Си1 - 0,29№1;
R = 0,96; F (7,31) = 46,15; Fst = 2,31; (6) Znг2 = 87,15 + !,08Сг1 - 59,51 + 1412 - 0,11РЬМ -- 0,05Сг2 - 0,16Си1;
R = 0,93; F (6,38) = 41,4; Fst = 2,35; (7) Znrз = 106,6 - 56,81 + 10,312 - 3,19Си + 0,17Си2;
R = 0,9; F (1,40) = 45,0; Fst =4,08; (8) СОг1 = 2,98 + 0,92^ - 3,171;
R = 0,98; F (2,36) = 459,9; Fst = 3,26; (9) СОг2 = 6,76 - 4,271 + Ода + 0,008СиСг;
R = 0,94; F (3,41) = 113,13; Fst =2,83; (10)
Согз = 2,52 + 0,59^ - !,831 - 0,0Ш + 0,27Сг -- 0,07101 - 0,007Сг2; R = 0,71; F (6,38) = 6,35; Fst = 2,35; (11) Реп = 414,11 + 136,712 - 459,491; R = 0,98; F(2,36) = 491,84; Fst = 3,26; (12) Реп = 574,83 - 352,161 + 59,6712; R = 0,94; F (2,42) = 155,52; Fst = 3,22; (13) Рег3 = 259,36 - 165,271 + 31,4312 +
+ 1,09Сг1 - 0,46NiCг; R = 0,87; F(4,40) = 30,57; Fst = 2,61; (14) Мпг1 = 104,23 + 12,29Сг - !,04Сг2 + 22,3212 -- 84,251 - !,6Сг1; R = 0,88; F (5,33) = 22,29; Fst = 2,49; (15) Мпг2 = 437,74 - 286,091 + 50,2812 + 31,28РЬ -- !,28РЬ2 - 4,12РЬ1 - 0,76NiCг; R = 0,95; F (6,38) = 57,28; Fst = 2,35; (16) Мпг3 = 69,64 + 94,31 - 20,3112 + 0,56Сг1 -- 0,19Сг2 + 2,61Сг; R = 0,97; F (5,39) = 134,22; Fst = 2,45; (17) Сиг1 = 10,54 - 6,231 + 1,1312 + 0,03№Сг + 0,08Сг1; R = 0,88^ (4,34) = 30,66; Fst = 2,65; (18) Сиг2 = 6,34 - !,531 + 0,36Сг - 0,02Сг2;
R = 0,8; F (3,41) = 23,68; Fst = 2,83; (19) Сиг3 = 8,8! - 4,171 - 0,41Сг + 0,7712 + + 0,01Сг2 + 0,04СП; R = 0,89; F (5,39) = 28,39; Fst = 2,45. (20) При этом значение Р-критерия Фишера в деся-тифакторной и пятифакторной моделях взято при уровнях значимости а = 5%.
Анализируя полученные пятифакторные и де-сятифакторные модели, можно констатировать, что при повышении концентрации С^ РЬ, N и Сг в почве их поступление, а также транслокация микроэлементов-метаболитов ^п, Ре, Си, Со, Мп) в растении нелинейно зависят от превышения кларков металлов в почве, что подтверждается коэффициентами детерминации. Последние дают право считать, что показатели содержания тяжелых металлов и металлов-метаболитов в растениях достаточно тесно связаны (г > 0,5), а в случае поступления С^ N РЬ, Сг в растения - сильно связаны (г > 0,7) с варьированием уровней поступления металлов в почву и фактором времени, что наглядно подтверждает прочность влияния факторов на результативный признак. Кроме того,
четко прослеживается специфика выращиваемой культуры, которая проявляется в разном уровне накопления металлов-токсикантов и разных причинах, обуславливающих данное явление. Так, зависимость поступления того или иного металла-токсиканта, как и микроэлементов-метаболитов в тестовые культуры в каждом конкретном случае своя, что подтверждается полученными моделями.
Поступление микроэлементов-метаболитов в тестовые культуры связано преимущественно с разным уровнем поступления металлов-токсикантов РЬ, N Сг) в течение вегетационного периода, с взаимодействием металлов-токсикантов с микроэлементами в почвенном растворе.
Выводы
1. Благодаря использованию многомерных методов статистического анализа установлено, что содержание микроэлементов-метаболитов в системе почва-растение - это переменный признак, который представляет собой функцию многих переменных и множества факторов, влияющих на этот признак. Использованные виды анализа информации позволяют утверждать, что зависимости в анализируемой системе имеют не функциональный, а статистический характер.
2. С помощью метода математического моделирования разработаны десяти- и пятифакторные регрессионные модели суммарного воздействия техногенного полиэлементного загрязнения почвы тяжелыми металлами на процессы транслокации микроэлементов, что позволило найти математическую формулу для имеющихся составляющих и их взаимной связи в системе почва-растение, а также перейти от аналитического изучения отдельных функций и процессов к синтетическому познанию целостности системы почва-растение, предвидеть возможные изменения одного признака на основе известных изменений другого.
3. Результаты математического моделирования показывают, что зависимости транслокации микроэлементов-метаболитов, тяжелых металлов в тестовые культуры в условиях техногенеза не идентичные, а специфические. Тяжелые металлы
влияют не столько количеством элементов, которые поступают в почву, сколько изменением их соотношения в почве, которое стимулируется данным влиянием.
4. Определение связей в системе почва-растение при транслокации микроэлементов и тяжелых металлов имеет большое значение для теоретического обоснования возможных уровней техногенных нагрузок на агроэкосистемы в условиях техногенного прессинга с целью оценки уровней поступления тяжелых металлов, нормирования загрязнителей и получения продукции агроценозов надлежащего качества.
Библиографический список
1. Борисочкина Т.И., Водяницкий Ю.Н. Загрязнение агроландшафтов России тяжелыми металлами: источники, масштабы, прогнозы // Бюллетень Почвенного института им. В.В. Докучаева. -2007. - № 60.
2. Капралова О.А. и др. Влияние загрязнения тяжелыми металлами на эколого-биологические свойства почв г. Ростова-на-Дону // Научная мысль Кавказа. - 2012. - № 1. - С. 13-13.
3. Водяницкий Ю.Н. Формулы оценки суммарного загрязнения почв тяжелыми металлами и металлоидами // Почвоведение. - 2010. - № 10. -С. 1276-1280.
4. Фадеев А.И., Мирошниченко Н.Н., Само-хвалова В.Л. Миграция, транслокация и фитоток-сичность тяжелых металлов при полиэлементном загрязнении почв // Агрохимия. - 2001. - № 3. -С. 57-61.
5. Жигарева Т.Л. и др. Влияние природных мелиорантов и тяжелых металлов на урожайность зерновых культур и микрофлору дерново-подзолистой почвы // Агрохимия. - 2005. - № 11. - С. 60-65.
6. Карпухин М.М., Ладонин Д.В. Влияние компонентов почвы на поглощение тяжелых металлов в условиях техногенного загрязнения // Почвоведение. - 2008. - № 11. - С. 1388-1398.
7. Новиков А.И. Эконометрика: учебное пособие. - М.: Дашков и К°, 2017. - 223 с.
8. Сыромятников Ю.Н. - Режим доступа: URL: https://orcid.org/0000-0001-9502-626X.
9. Лакин Г.Ф. Биометрия. - М.: Высшая школа, 1990. - 352 с.
10. Математическое моделирование в биологии и химии. Новые подходы: сб. науч. тр. / отв. ред. Р.Г. Хлеборос; Рос. акад. наук. Сиб. отд. -Новосибирск: Наука; Сиб. отд., 1992. - 219 с.
11. Егоршин А.А., Малярец Л.М. Корреляционно-регрессионный анализ. - Харьков: Основа, 1998. - 201 с.
12. Задорин А.Д., Рубинштейн М.И., Юмагу-лова А.Н., Андрейчук А.Л. О моделировании связей в системе почва-растение // Почвоведение. -1982. - № 8. - С. 142-145.
13. Марчук Г.И. Математическое моделирование и охрана природы // Природа. - 1985. - № 7. -С. 6-12.
14. Математические модели рационального природопользования: Сб. науч. тр. / отв. ред. В.В. Пененко, И.Б. Токин; АН СССР, Сиб. отд. -Новосибирск: Наука, Сиб. отд., 1989. - 137 с.
15. Пененко В.В., Алоян А.Е. Модели и методы для задач охраны окружающей среды / отв. ред. М.М. Лаврентьев. - Новосибирск: Наука; Сиб. отд., 1985. - 256 с.
References
1. Borisochkina T.I., Vodyanitskiy Yu.N. Zagrya-znenie agrolandshaftov Rossii tyazhelymi metallami: istochniki, masshtaby, prognozy // Byulleten Poch-vennogo instituta im. V.V. Dokuchaeva. - 2007. -No. 60.
2. Kapralova O.A. i dr. Vliyanie zagryazneniya tyazhelymi metallami na ekologo-biologicheskie svoystva pochv g. Rostova-na-Donu // Nauchnaya mysl Kavkaza. - 2012. - No. 1. - S. 13-13.
3. Vodyanitskiy Yu.N. Formuly otsenki summar-nogo zagryazneniya pochv tyazhelymi metallami i metalloidami // Pochvovedenie. - 2010. - No. 10. -S. 1276-1280.
4. Fadeev A.I., Miroshnichenko N.N., Samokh-valova V.L. Migratsiya, translokatsiya i fitotoksichnost tyazhelykh metallov pri polielementnom zagryaznenii pochv // Agrokhimiya. - 2001. - No. 3. - S. 57-61.
5. Zhigareva T. L. i dr. Vliyanie prirodnykh melio-rantov i tyazhelykh metallov na urozhaynost zerno-vykh kultur i mikrofloru dernovo-podzolistoy pochvy // Agrokhimiya. - 2005. - No. 11. - S. 60-65.
6. Karpukhin M.M., Ladonin D.V. Vliyanie kom-ponentov pochvy na pogloshchenie tyazhelykh metallov v usloviyakh tekhnogennogo zagryazneniya // Pochvovedenie. - 2008. - No. 11. - S. 1388-1398.
7. Novikov A.I. Ekonometrika: uchebnoe posobie. - M.: Dashkov i K°, 2017. - 223 s.
8. Syromyatnikov Yu.N. [Elektronnyy resurs]. -URL: https://orcid.org/0000-0001-9502-626X.
9. Lakin G.F. Biometriya. - M.: Vysshaya shkola, 1990. - 352 s.
10. Matematicheskoe modelirovanie v biologii i khimii. Novye podkhody: Sb. nauch. tr. / Ros. akad. nauk. Sib. otd.; Otv. red. R.G. Khleboros. - Novosibirsk: Nauka, Sib. otd., 1992. - 219 s.
11. Yegorshin A.A., Malyarets L.M. Korrelya-tsionno-regressionnyy analiz. - Kharkov: Osnova, 1998. - 201 s.
12. Zadorin A.D., Rubinshteyn M.I., Yumagulo-va A.N., Andreychuk A.L. O modelirovanii svyazey v sisteme pochva-rastenie // Pochvovedenie. - 1982. -No. 8. - S. 142-145.
13. Marchuk G.I. Matematicheskoe modelirovanie i okhrana prirody // Priroda. - 1985. - No. 7. -S. 6-12.
14. Matematicheskie modeli ratsionalnogo pri-rodopolzovaniya: Sb. nauch. tr. / AN SSSR, Sib. otd.; otv. red. V.V. Penenko, I.B. Tokin. - Novosibirsk: Nauka, Sib. otd., 1989. - 137 s.
15. Penenko V.V., Aloyan A.Ye. Modeli i metody dlya zadach okhrany okruzhayushchey sredy / otv. red. M.M. Lavrentev. - Novosibirsk: Nauka, Sib. otd., 1985. - 256 s.
+ + +