Статья: Математическая модель прогноза урожайности лекарственного сырья чистотела большого (Chelidonium majus L.)

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

2

Костромская государственная сельскохозяйственная академия

Математическая модель прогноза урожайности лекарственного сырья чистотела большого (Chelidonium majus L.)

УДК 633.8

Бородий С.А.

Аннотация

Предмет исследования -- дикорастущая ценопопуляция чистотела большого (Chelidonium majus L.). Цель исследований -- разработать математическую модель сезонного прогноза фитометрических параметров генеративного побега чистотела большого в дикорастущих ценопопуляциях для Костромской области. Исследования проводили в период 1995…2018 гг. (ареал распространения чистотела большого по территории Костромской области) и 2014…2015 гг. (сбор полевой информации, разработка модели прогноза). Установлено, что чистотел большой встречается на всей территории Костромской области как рудерально-сегетальный вид. Уборку надземной массы этого растения целесообразно проводить, когда первые плоды находятся в стадии начала созревания -- начала рассеивания семян, календарные сроки которых прогнозируются с точностью 99,77%. Прогнозы динамики фитометрических параметров стебля совпадали с эмпирическими данными на: 98,6…98,8 % (высота побега); 95,2…98,0 % (надземная общая фитомасса); 93,0…93,7 % (масса листьев); 94,5…97,8 % (масса стебля); 80,8…91,0 % (масса генеративных органов), что обеспечивает возможность применения модели в производственных условиях для прогноза урожайности дикорастущих ценопопуляций. Прогноз массы семян был менее точен (69,8…87,6 %), поэтому может применяться только для приблизительного расчета.

Ключевые слова: ЧИСТОТЕЛ БОЛЬШОЙ, CHELIDONIUM MAJUS L., ДИНАМИКА ПРОДУКТИВНОСТИ ПОБЕГА, ПРОГНОЗ УРОЖАЙНОСТИ, МАТЕМАТИЧЕСКАЯ МОДЕЛЬ

генеративный побег надземная масса чистотел

Введение

Надземную массу чистотела большого применяют в официальной и народной медицине как спазмолитическое, желчегонное, противовоспалительное и ранозаживляющее средство. Препараты чистотела эффективны для лечения папилломатоза, начальных стадий красной волчанки, задерживают рост злокачественных опухолей. В ветеринарной медицине свежей, смешанной с солью травой лечат тимпанию овец [1]. Настой свежей травы можно применять для защиты растений от вредителей и болезней [2], поскольку это растение обладает фунгистатическим и инсектицидным действием.

Запасы сырья чистотела значительно превышают его потребность. В СССР заготовительными конторами объединения «Союзлекраспром» ежегодно заготовлялось более 25 тыс. т дикорастущего сырья, часть которого экспортировалась. В настоящее время объемы сократились в десятки раз [2], несмотря на то, что потребность остается высокой. В связи с увеличением количества сельскохозяйственных предприятий, входящих в ассоциацию биологического земледелия, когда применение химических средств защиты растений запрещается, объемы заготовки травы чистотела будут возрастать.

Технология выращивания чистотела в культуре разработана [2], но основной объем сырья пока получают из дикорастущих ценопопуляций. Для планирования периода заготовительных работ и ожидаемой урожайности эксплуатационных участков требуется расчет сезонного прогноза динамики фенологических стадий развития растений и урожайности сырья стандартной влажности минимум за месяц до начала заготовки. В настоящее время такой прогноз для чистотела большого отсутствует. Таким образом, целью исследований ставилось разработать математическую модель сезонного прогноза фитометрических параметров генеративного побега чистотела большого в дикорастущих ценопопуляциях для Костромской области.

Исследования проводили в период 1995…2018 гг. (ареал распространения чистотела большого по территории Костромской области) и 2014…2015 гг. (сбор полевой информации, разработка модели прогноза). Ареал распространения изучали путем маршрутного обследования. Места произрастания отмечали на карте точками. Изучение основных фитометрических параметров растения (высота побега, масса листьев стандартной (13%) влажности) проводили на участке антропогенно преобразованной территории в окрестностях д. Апраксино Костромского района Костромской области. Наблюдения и учеты выполняли с интервалом 7 суток от весеннего отрастания побегов до начала рассеивания семян из сформировавшихся первыми плодов, поскольку в более поздние стадии развития нижние листья желтеют и отмирают, что снижает качество сырья. Выборка для учета включала 10 побегов, срезанных в случайном порядке на уровне почвы в четырехкратной повторности с последующим измерением высоты. Для изучения динамики массы фитоорганы отделяли от побега, высушивали до влажности 13% и взвешивали. Параметром времени в прогностической модели служила суммарная энтальпия воздуха, которую рассчитывали и прогнозировали на весенне-летний период текущего года по оригинальной методике [3]. Статистическая обработка результатов исследований и разработка модели прогноза выполнены с использованием табличного процессора Excel.

Результаты исследований

Географическое распространение. По опубликованным данным [1, 2], чистотел большой как евроазиатский вид широко распространен в европейской части России, на Кавказе, в Сибири (кроме Арктической зоны), на Дальнем Востоке. Экологические условия благоприятны в населенных пунктах, по окраинам полей, кустарникам, лесосекам, вдоль канав и др. на легко- и среднесуглинистых почвах с высоким содержанием гумуса.

Результаты обследования части территории Костромской области совпадают с данными публикаций, то есть и здесь вид относится к группе рудерально-сегетальных (рис. 1). Следовательно, присутствие чистотела большого в фитоценозах может служить индикатором присутствия антропогенного фактора на определенной территории.

Рис. 1. Географическое распространение дикорастущих ценопопуляций чистотела большого (Chelidonium majus L.) в Костромской области, 1995…2018 гг. (ориг.)

Фенологические стадии развития. Биологической особенностью чистотела большого является продолжительное по времени развитие. Стадия бутонизации как начало генеративного периода идентифицируется хорошо, но после начала формирования первых плодов последующие стадии наблюдаются практически одновременно, поскольку стебель с листьями и цветоносами продолжает рост из пазухи верхнего листа (рис. 2).

Рис. 2. Генеративный побег чистотела большого (Chelidonium majus L.) с генеративными органами различных стадий развития (11.07.2015) (ориг.)

По нашим наблюдениям цветение и плодоношение продолжаются до конца августа -- начала сентября, но после массового цветения (конец мая -- начало июня) в срезанной зеленой массе будут встречаться отмершие листья, что потребует дополнительных затрат на последующую ручную сортировку лекарственного сырья.

Следовательно, уборку надземной массы чистотела целесообразно проводить, когда первые плоды находятся в стадии начала созревания -- начала рассеивания семян. Для планирования календарного срока уборки мы рекомендуем применять прогноз суммарной энтальпии воздуха [3], для которого требуется знать теплоемкость стадий развития. Наблюдениями за динамикой развития генеративного побега чистотела в дикорастущей ценопопуляции было установлено значение теплоемкости фенологических стадий с высоким совпадением по годам (99,77%) (табл. 1).

Таблица 1. Теплоемкость фенологических стадий развития чистотела большого (Chelidonium majus L.) в дикорастущей ценопопуляции

Фенологическая стадия развития

Календарная дата

Суммарная энтальпия воздуха, КДж/кг

2014 г.

2015 г.

Среднее

2014 г.

2015 г.

Среднее

Бутонизация

10.05

09.05

09.05…10.05

423,37

394,65

409,01

Начало цветения

17.05

16.05

16.05…17.05

632,84

598,35

615,60

Начало образования семян

24.05

23.05

23.05…24.05

908,43

808,29

858,36

Начало созревания семян

31.05

30.05

30.05…31.05

1162,49

1085,09

1123,79

Начало рассеивания семян

07.06

06.06

06.06…07.06

1450,04

1344,59

1397,32

R2

0,9977

Рассчитанными значениями теплоемкости можно руководствоваться для прогноза календарных сроков стадий развития растений и, следовательно, планирования сроков заготовки лекарственного сырья чистотела.

Динамика высоты побега. Отрастание побегов дикорастущей ценопопуляции чистотела как многолетнего растения начиналось вскоре после схода снега, что наблюдалось в начале мая. В дальнейшем рост продолжался даже после стадии рассеивания семян из сформировавшихся первыми плодов, что объясняется формированием стебля с цветоносами из пазух верхних листьев. Таким образом, интенсивный рост продолжался до третьей декады июня (рис. 3).

Динамика высоты побега (hр) подчинялась линейной зависимости, поэтому для прогноза была разработана следующая модель:

(180,0?Ср ? 1900,0 КДж/кг)(1)

,

где: h'р -- корректирующий коэффициент; hfр -- эмпирическая высота побега на дату учёта, см; Cpf -- эмпирическая суммарная энтальпия воздуха от даты весеннего возобновления вегетации до даты учета высоты побега, КДж/кг; Cp -- прогнозируемая суммарная энтальпия воздуха, Кдж/кг.

Рис. 3. Динамика высоты побега чистотела большого (Chelidonium majus L.) в дикорастущей ценопопуляции (ориг.)

Ретроспективный расчет высоты побега показал высокое совпадение с эмпирическими данными: 98,8% в 2014 г. и 98,6% в 2015 г., что позволяет рекомендовать модель для прогноза высоты травостоя чистотела (табл. 2).

Таблица 2. Верификация модели прогноза высоты побега чистотела большого (Chelidonium majus L.) в дикорастущей ценопопуляции

Календарная дата учета

Суммарная энтальпия воздуха, кДж/кг

Высота побега, см

эмпирическая

расчетная

03.05.2014

321,98

14,25

14,25

10.05.2014

423,37

20,72

18,98

17.05.2014

632,84

32,10

28,77

24.05.2014

908,43

41,15

41,64

31.05.2014

1162,49

62,40

53,50

07.06.2014

1450,04

71,62

66,93

14.06.2014

1711,11

80,32

79,12

21.06.2014

1906,81

89,15

88,26

R2

0,9883

02.05.2015

186,82

13,29

13,29

09.05.2015

394,65

19,06

23,00

16.05.2015

598,35

25,17

32,51

23.05.2015

808,29

39,33

42,31

30.05.2015

1085,09

51,18

55,24

06.06.2015

1344,59

69,46

67,36

13.06.2015

1546,36

70,25

76,78

20.06.2015

1801,54

83,69

88,70

R2

0,9859

Динамика надземной общей фитомассы побега. Стабильное увеличение фиотомассы побега наблюдалось до начала рассеивания семян, далее процесс имел тенденцию к замедлению и снижению, что особенно заметно в 2014 году (рис. 4).

Рис. 4. Динамика надземной общей фитомассы побега чистотела большого (Chelidonium majus L.) в дикорастущей ценопопуляции (ориг.)

Нестабильность динамики после стадии начала рассеивания семян объясняется отмиранием и опадением нижних листьев и созревших семян из коробочек. Известно, что лекарственным сырьем является вся зеленая надземная фитомасса, причем в стадию массового цветения созревшие семена нижних коробочек используются как посевной материал для закладки культурных плантаций [2]. Следовательно, для прогноза количества заготовляемой фитомассы диапазон работы модели можно ограничить стадией созревания -- начала рассеивания семян, которая в годы исследования наблюдалась в июне.

Для дикорастущих ценопопуляций разработана модель на базе системы уравнений:

(180,0 ? Cp ? 630,0 КДж/кг),

(180,1 ? Cp ? 1900,0 КДж/кг),(2)

где Mo -- общая надземная фитомасса побега, г.

Прогноз Mо в зависимости от Cp рассчитывается по модели:

(3)

,

где: Mо' -- корректирующий коэффициент; Mоf -- эмпирическая общая надземная фитомасса на календарную дату учёта, г/побег; f(Cp), f(Cpf) -- функция динамики воздушно-сухой (13%) общей надземной фитомассы побега (см. (2)).

Корреляционный анализ показал хорошее совпадение эмпирических и прогнозируемых значений (табл. 3), что позволяет рекомендовать модель (3) для прогноза динамики общей надземной фитомассы побега.

Таблица 3. Верификация модели прогноза массы (13% влажность) побега чистотела большого (Chelidonium majus L.) в дикорастущей ценопопуляции

Календарная дата учета

Суммарная энтальпия воздуха, кДж/кг

Общая надземная фитомасса (Mо), г/побег

Масса листьев (Ml), г/побег

Масса стебля (Ms), г/побег

Масса соцветий и плодов (Mg), г/побег

Масса семян (Mgs), г/побег

Э*

П*

Э*

П*

Э*

П*

Э*

П*

Э*

П*

03.05.2014

321,98

0,29

0,29

0,16

0,15

0,13

0,13

10.05.2014

423,37

0,79

0,62

0,36

0,29

0,26

0,20

0,17

0,13

17.05.2014

632,84

1,55

1,70

0,45

0,49

0,83

0,91

0,28

0,28

24.05.2014

908,43

3,12

2,64

1,26

1,07

1,42

1,18

0,44

0,36

31.05.2014

1162,49

3,53

3,47

1,14

1,17

1,71

1,62

0,36

0,41

0,32

0,33

07.06.2014

1450,04

4,42

4,10

1,21

1,19

2,07

1,99

**

0,43

0,54

0,50

14.06.2014

1711,11

3,74

4,38

1,03

1,23

1,77

2,18

0,40

0,42

0,54

0,49

21.06.2014

1906,81

3,99

4,41

1,11

1,28

2,00

2,21

0,44

0,40

0,43

0,47

R2

0,9515

0,9303

0,9445

0,8079

0,8758

02.05.2015

186,82

0,25

0,25

0,12

0,12

0,13

0,13

09.05.2015

394,65

0,63

0,72

0,26

0,34

0,22

0,27

0,15

0,16

16.05.2015

598,35

1,80

1,69

0,64

0,55

0,96

0,88

0,20

0,29

23.05.2015

808,29

2,29

2,45

0,80

0,86

1,22

1,27

0,27

0,36

30.05.2015

1085,09

2,87

3,45

1,10

1,22

1,26

1,58

*

0,43

0,28

0,31

06.06.2015

1344,59

3,67

4,11

1,14

1,25

1,83

1,97

0,39

0,45

0,31

0,46

13.06.2015

1546,36

4,28

4,44

1,32

1,26

2,10

2,18

0,40

0,45

0,46

0,52

20.06.2015

1801,54

4,73

4,62

1,57

1,31

2,35

2,30

0,41

0,43

0,40

0,50

R2

0,9797

0,9367

0,9776

0,9103

0,6984