ДИНАМИКА УРОЖАЙНОСТИ УЛУЧШЕННЫХ СЕНОКОСОВ

Введение
Для функционального описания поведения территории нами вводятся новые понятия - активность и интенсивность растительного покрова. Причем территория понимается как простейшее геодезическое изображение ландшафта. А сам ландшафт, в свою очередь, является первым компонентом динамической геотриады «ландшафт + население + хозяйство».
Активность учитывается по доле площади растительного покрова (леса и древесно-кустарниковая растительность, луга и пастбища, особо охраняемые территории и болота) и этот экологический параметр позволяет хаpaктеризовать фактически образовавшиеся отклонения от территориального экологического равновесия [1] на конкретной территории.
Однако пока отсутствуют данные земельных кадастров по конкретным землепользователям. Поэтому приведены значения параметров растительного покрова по административным районам и городам субъекта федерации. В частности рассмотрены районы и города Республики Марий Эл (РМЭ) по состоянию распределения земель на 01.01.07 г.
В наиболее общем случае интенсивность проявляется как активность во времени. Физически интенсивность - это скорость изменений. А активность - это сами изменения в природной, природно‒техногенной или технической среде (по площади, урожайности растений, продуктивности почвы и пр.) в некотором срезе времени.
Ритмика жизнедеятельности растений четко предопределена циклами вращения Земли вокруг Солнца, поэтому единицей измерения любого показателя (параметр - это показатель, хаpaктеризующий исследуемую систему) растительного покрова и его компонент станет год. Активность и интенсивность поведения человека по отношению к процессу уничтожения ландшафта быстротечны. Например, планировалось атомными взрывами построить русло крупного канала для переброски воды северных рек на юг.
Растительный покров - это территориально распластанное живое существо, обладающее многими свойствами поведения живого вещества. Индикатором экологического равновесия на данной территории становится его активность. Она в простейшем случае исчисляется по доле занятой площади, причем в ходе конкуренции за плодородные почвы, с человеком. А интенсивность поведения (скорость реакции на внешние раздражения, в том числе и на антропогенные воздействия) растительного покрова за более чем 350‒450 млн. лет эволюции наземных растений пока не изменилась. За скоростью эволюции человека растения явно не успевают. Но они это возместят высокой своей активностью роста и развития сорняков, а также ликвидацией нужных человеку видов культур болезнями.
Травяной покров. Особое место в растительном покрове занимает травяной покров, и он более живуч, устойчив и стабилен в своем поведении по сравнению с древесными пологами. По данным американских ученых, трава появилась около 100 млн. лет назад. А деревья возникли 350‒450 млн. лет назад, то есть трава появилась как результат отклика древесных растений на поведение наземных животных. Поэтому трава биологически устойчивее деревьев, кустарников и кустарничков. При снижении плодородия почвы и ухудшении климата деревья исчезают, а их место обитания занимают степные травы.
Интенсивность воспроизводства травяного покрова. Годичная цикличность определяет урожайность травы в естественных условиях роста и развития, какими являются естественные луга, на которых заготовляется сено. Улучшенные луга уже меняют свое поведение, повышая урожайность на позитивное воздействие человека культурными и техническими мерами и внесением на такие луга удобрений.
Из статистического сборника [6, с. 181] были взяты сводные данные по сенокосным лугам Республики Марий Эл (РМЭ), приведенные в таблице 1.
Таблица 1
Урожайность сена в сельских хозяйствах на территории РМЭ, ц/га
|
Годы учета* |
в хозяйствах всех категорий |
в крупных и средних хозяйствах |
|||||
|
Время t, лет |
Естественные и улучшен. сенокосы |
Однолетние травы |
Много-летние травы |
Время t, лет |
Естественные сенокосы |
Улучшенные сенокосы |
|
|
1970 1975 1980 1985 1990 1995 1996 1997 1998 1999 2000 2001 2002 2003 |
0 5 10 15 20 25 26 27 28 29 30 31 32 33 |
14.5 8.0 12.0 17.1 19.0 3.0 12.2 16.0 11.3 9.9 13.6 14.6 9.1 12.5 |
12.4 7.9 20.5 29.9 26.9 16.1 30.4 18.8 21.4 11.1 33.7 20.6 17.7 21.1 |
19.9 10.1 18.6 35.1 33.6 16.8 18.8 21.3 16.1 15.5 18.2 19.0 12.1 18.6 |
- - 0 5 10 15 16 17 18 19 20 21 22 23 |
- - 11.0 14.7 16.7 12.1 10.0 13.0 10.2 9.5 12.0 14.1 8.7 12.8 |
- - 16.3 21.8 26.2 15.4 15.2 21.1 14.6 11.3 18.8 15.5 10.4 11.2 |
Примечание: * Недостатком динамического ряда является отсутствие данных по всем годам в промежутках между пятилетками, что не позволяет точно выявить циклы с полупериодом менее пяти лет.
Вначале рассмотрим статистические закономерности динамики урожайности травы после естественной сушки в сено на сенокосных лугах сельских хозяйств всех категорий. При этом будем выполнять все этапы статистического моделирования:
‒ эвристическая идентификация, когда выясняется логика каждого динамического ряда исходных количественных данных;
‒ структурная идентификация, когда по эвристической модели понимание процессов изменения значений изучаемого показателя выбирается основная (естественная составляющая) закономерность, то есть тренд изменения во времени исследуемого показателя;
‒ параметрическая идентификация, когда в решающей математической среде типа CurveExpert‒1.3 подбираются значения параметров статистической модели исходной конструкции, полученной на предыдущем этапе структурной идентификации.
После выполнения всех процедур и этапов идентификации биотехнического закона, предложенного проф. П.М. Мазуркиным [2‒5], сформируется готовая статистическая модель. Урожайность сена с естественных и улучшенных сенокосов со всей площади по РМЭ определяется рядом значений по точечному графику на рис. 1.
Рис. 1. Изменение урожайности сена естественных и улучшенных сенокосов
Республики Марий Эл, ц/га
Из роя точек резко выделяется одна (самая нижняя) точка.
Полнота динамического ряда. Пусть исходные данные [6, с.181] количественно достоверны и статистики адекватно описали табличными моделями типа табл. 1 картину динамики исследуемого явления или процесса. Тогда можно поставить задачу оценки полноты статистического ряда динамики значений изучаемого показателя.
Критерий полноты значений показателя. Представим динамический ряд значений урожайности сена с естественных и улучшенных сенокосов со всей площади по РМЭ в виде последовательности, показанной в таблице 2.
Весь этот ряд динамики можно разделить на две части:
а) динамика по пятилетним циклам, причем четко по хронологическому времени 1970, 1975, 1980 и т.д.;
б) динамика по годичным циклам, что соответствует действительности, то есть ритму вращения Земли вокруг Солнца.
Таблица 2
Полный динамический ряд урожайности сена, ц/га
|
Время t, лет |
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
|
|
Урож. q, ц/га |
14.5 |
- |
- |
- |
- |
8.0 |
- |
- |
- |
- |
|
|
10 |
11 |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
19 |
20 |
21 |
|
12.0 |
- |
- |
- |
- |
17.1 |
- |
- |
- |
- |
19.0 |
- |
|
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
33 |
|
- |
- |
- |
3.0 |
12.2 |
16.0 |
11.3 |
9.9 |
13.6 |
14.6 |
9.1 |
12.5 |
Введем новый показатель оценки добротности статистических данных динамического ряда - это коэффициент полноты ряда. Этот коэффициент k будет вычисляться как отношение количества m имеющихся значений показателя в ряду к общему количеству возможных члeнов n динамического ряда, то есть по формуле k=m/n .
Из данных табл. 2 получаем m=14 и n=34. Поэтому коэффициент полноты ряда будет k=14/34=0,412. Динамический ряд оказался сильно деформированным. Специальными уловками группировки исходных данных вовсю занимались советские экономисты, чтобы доказывать недоказуемое. Нужно статистическим службам запретить публикации неполных динамических рядов, так как это почти всегда приводит к ложной идентификации.
Картина распределения точек, показанная роем на рисунке 1, оказалась ложной из-за отсутствия 20 точек из требуемых 34 значений изучаемого показателя. При этом не ясно, какими по расположению относительно принятой системы координат были точки в промежутках между пятилетками. Поэтому может оказаться, что полный статистический ряд даст совершенно иную модель динамики (если урожайности 1970, 1975, 1980, 1985, 1990, 1995 годов окажутся не максимальными, а иными). Вот почему надо обращать серьезное внимание на полноту статистического ряда динамики значений изучаемого показателя.
Границы доверительного интервала статистического ряда. Применим известный прием выделения верхней и нижней границ доверительного интервала. Это обезопасит от ложной идентификации и неправильной тpaктовки найденной формулы, так как верхняя и нижняя границы показывают только возможное поле распределения значений показателя вне зависимости от полноты статистического ряда.
Перепишем первую часть показателей таблицы 1 с учетом границ доверительного интервала и результаты приведем в таблице 3.
Верхняя граница урожайности сена естественных и улучшенных сенокосов определится уравнением (рис. 2) вида
qmax=14,5092exp(-0,0091099t)+4,7361×10-8t9,23434exp(-0,44353t). (1)
Максимальная относительная погрешность уравнения (1) равна 4,59 %. Коэффициент корреляции по рис. 2 равен 0,9861.
Нижняя граница урожайности сена (рис. 3) определится статистическим уравнением
qmin = 7.64985exp(-0,0052882t)+0,000403119 t 5,45669 exp(-0,31879t). (2)
Максимальная относительная погрешность уравнения (2) равна 3,84 %, а коэффициент корреляции по рисунку 3 равен 1,0000, то есть формула (2) становится однозначно функциональной математической зависимостью.
Таблица 3
Урожайность сена в сельских хозяйствах всех категорий на территории РМЭ, ц/га
|
Годы учета |
Время t, лет |
Естественные и улучшенные сенокосы |
Однолетние травы |
Многолетние травы |
||||||
|
ряд |
max |
min |
ряд |
max |
min |
ряд |
max |
min |
||
|
1970 1975 1980 1985 1990 1995 1996 1997 1998 1999 2000 2001 2002 2003 |
0 5 10 15 20 25 26 27 28 29 30 31 32 33 |
14.5 8.0 12.0 17.1 19.0 3.0* 12.2 16.0 11.3 9.9 13.6 14.6 9.1 12.5 |
14.6
17.1 19.0
16.0
14.6
12.5 |
8.0 12.0
12.2
9.9
9.1 |
12.4 7.9 20.5 29.9 26.9 16.1 30.4 18.8 21.4 11.1 33.7 20.6 17.7 21.1 |
12.4
20.5 29.9
30.4
33.4
21.1 |
7.9
16.1
18.8
11.1
17.7 |
19.9 10.1 18.6 35.1 33.6 16.8 18.8 21.3 16.1 15.5 18.2 19.0 12.1 18.6 |
19.9
35.1 33.6
21.3
19.0
18.6 |
10.1
16.8 18.8
16.1 15.5
12.1 |
Примечание: * Такие резко отклоняющиеся точки исключаются из статистических рядов.
Рис. 2. Верхняя граница доверительного интервала урожайности сена естественных
и улучшенных сенокосов РМЭ, ц/га
Закономерности динамики урожайности сена. Конструкция уравнений для описания верхней и нижней границы доверительного интервал одинакова. Она содержит две составляющие:
а) закон гибели (спада) в упрощенной конструктивной форме, когда третий параметр статистической модели будет равен единице;
б) биотехнический закон проф. П.М. Мазуркина [2‒5] в упрощенной форме (когда четвертый параметр модели равен единице).
Когда известна общая конструкция статистической закономерности, то можем её применить ко всему динамическому ряду новых исходных статистических данных.
Рис. 3. Нижняя граница доверительного интервала урожайности сена естественных и улучшенных сенокосов РМЭ, ц/га
В результате идентификации был получен тренд (рис. 4) из двух составляющих в виде устойчивых законов, в котором знак экспоненциального закона изменился с отрицательного на положительный, и он стал показывать рост значений урожайности сена. При этом 1995 год не был исключен.
q=10,9017exp(0,00083530t)+1,1128×10-10 t 14,0597 exp(-0,88382t), (4)
После исключения резко отклоняющейся от других значений показателя точки, отмеченной в данных табл. 3 звёздочкой, была получена формула (рис. 5) статистической закономерности вида
q=11,4401exp(-0,0010562t)+4,9549×10-12 t 14,74586 exp(-0,80755t). (5)
Коэффициент корреляции увеличился до 0,7266, а закономерность стала одинаковой с границами доверительного интервала.
Рис. 4. Тренд урожайности сена естественных и улучшенных сенокосов РМЭ, ц/га
Рис. 5. Тренд урожайности сена естественных и улучшенных сенокосов РМЭ после исключения резко отклоняющейся точки, ц/га
При этом попытка наращивания первой составляющей до закона гибели в полной форме не дает результата. Остатки приведены на рис. 6.
Чередование знаков остатков ε показывает, что можно дополнить составляющей в виде закона колебательного движения.
Структурно-параметрическая идентификация выполняется относительно вейвлет‒функции проф. П.М. Мазуркина [2‒5] вида
q3 = Acos(πt/p0.5±a10), (6)
где A - амплитуда (половина) колебательного движения значений изучаемого показателя во времени по биотехническому закону
(6а)
p0.5 - половина периода колебательного возмущения показателя по биотехническому закону (с учетом предыстории в виде постоянного члeна)
, (6б)
a1...a10- параметры исходной конструкции статистической модели (6), являющейся обобщенным законом колебательного изменения значений изучаемого показателя.
Волновую закономерность нельзя искать по статистическим выборкам с пропусками данных. Поэтому принимается интервал времени с 1995 по 2003 гг. (причем результат 1995 года исключается).
После нескольких сеансов поиска значений параметров модели (6) в математической среде CurveExpert‒1.3 была получена конкретная закономерность в виде сложного уравнения
q=5,64537exp(0,022908t)+6,0178×10-13 t 14,8210 exp(-0,74130t)-A1cos(πt/p1+3,29024)+
+A2cos(πt/p2+1,31367), (7)
A1 = 2,09541 t 0,15870 exp(-0,0037692t), p1 = 2,79259 - 0,019856t,
A2 = 4,1925×10-9 t 8,82490 exp(-0,38095t), p2 = 1,12563 + 0,00014490t.
Рис. 7. Динамика урожайности сена естественных и улучшенных сенокосов РМЭ за период 1970‒2003 гг., ц/га
В статистической модели (7) из четырех составляющих две последние являются колебательными возмущениями. Из-за отрицательного знака первая волновая зависимость является кризисным возмущением сенокосных лугов на внешние раздражители (естественные и антропогенные), причем частота возмущения растет, а вторая волна с начальным периодом в 2×1,12563 = 2,25 лет показывает позитивное успокоение травяного покрова со снижением частоты адаптивного возмущения.
Сокращение ряда до 1996‒2003 гг. сужает прогнозные возможности готовой статистической модели. Исходя из высокой точности, когда максимальная относительная погрешность равна всего 0,47 %, уравнение (7) дает возможность принять горизонт прогноза, равный основанию прогноза, то есть в восемь лет до
2003 + 8 = 2011 года.
А при ряде со всеми 34 точками получили бы возможность прогноза до
2003 + 34 = 2037 года.
При этом остатки располагаются так, как это показано на рис. 8.
Расположение остатков показывает, что существует еще одна волновая составляющая с примерным периодом в пять лет. Погрешность измерений урожайности сена недостаточно высока, чтобы искать вейвлет‒функцию с амплитудой менее 0,08 ц/га (запись в справочнике [6, с.181] выполнена с округлением до одной десятой, поэтому точность измерений составляет ±0,05 ц/га).
Рис. 8. Остатки готовой статистической модели тренда (4.11) урожайности сена естественных и улучшенных сенокосов РМЭ, ц/га
Изменение шкалы времени. Если исключить пятилетние данные, а также резко отклоняющуюся точку 1995 г., то получим динамический ряд с новой шкалой абсцисс. Переход выполняется по формуле t1=t‒26.
Для новой шкалы получаем условие t1 = 0 для 1996 года. Такое сокращение времени приводит и к уменьшению количества составляющих модели.
После повторной параметрической идентификации модели (4.11) была получена формула (рис. 9) с тремя составляющими
q=12,3327exp(-0,00036540t) -A1cos(πt/p1+1,53433)+A2cos(πt/p2 -2,93294), (8)
A1 = 3,63948exp(0,0089534t), p1 = 1,91645-0,028705t,
A2 = 0,034478 t 8,51786 exp(-1,80140t), p2 = 1,16822 - 0,0086746t
Рис. 9. Динамика урожайности сена естественных и улучшенных сенокосов РМЭ за период 1996‒2003 гг., ц/га
Рис. 10. Остатки готовой статистической модели тренда (4.12) урожайности сена естественных и улучшенных сенокосов РМЭ за период 1996‒2003 гг., ц/га
По коэффициенту корреляции 1,000 и нулевым остаткам (рис. 10) модель (8) стала однозначной.
Даже по малым остаткам, находящихся в интервале пренебрежимо малых для пpaктики измерения массы сена чисел ǀɛǀ=7·10-12...7·10-10 на рис. 10 наблюдается новая волновая закономерность. Этот факт означает, что повышение точности измерения массы травы и сена с единицы площади позволит находить статистические закономерности динамики урожайности.
Список литературы
- Реймерс, Н.Ф. Природопользование: Словарь‒справочник / Н.Ф. Реймерс. - М.: Мысль, 1990. - 637 с.
- Мазуркин, П.М. Геоэкология: Закономерности современного естествознания: Научное изд. / П.М. Мазуркин. - Йошкар-Ола: МарГТУ, 2006. - 336 с.
- Мазуркин, П.М. Закономерности устойчивого развития / П.М. Мазуркин. ‒ Научное издание. - Йошкар-Ола: МарГТУ, 2002. - 302 с.
- Мазуркин, П.М. Рациональное природопользование: учебное пособие. В 3‒х ч. Ч. 1: Экологически ответственное землепользование / П.М. Мазуркин, С.Е. Анисимов, С.И. Михайлова; под ред. П.М. Мазуркина. - Йошкар-Ола: МарГТУ, 2006. - 176 с.
- Мазуркин, П.М. Математическое моделирование. Идентификация однофакторных статистических закономерностей: Учебное пособие / П.М. Мазуркин, А.С. Филонов. - Йошкар-Ола: МарГТУ, 2006. - 292 с.
- Республика Марий Эл: Статистический ежегодник / А.В. Целищев и др. - Часть 1. - Йошкар-Ола: ФСГС, 2004. - 342 с.
Статья подготовлена и опубликована при поддержке гранта 3.2.3/4603 МОН РФ
Статья в формате PDF
101 KB...
12 04 2026 15:53:53
Статья в формате PDF
133 KB...
11 04 2026 23:46:56
Статья в формате PDF
435 KB...
10 04 2026 1:17:32
Статья в формате PDF
130 KB...
09 04 2026 15:10:49
Одним из наиболее часто встречающихся осложнений после пластических операций остаются гипертрофические рубцы [1;6;10], этиология которых может быть обусловлена неадекватным образованием вазоактивных веществ. Репаративная регенерация операционной раны состоит из серии биохимических координированных реакций между различными типами клеток, регулируемых локальными медиаторами. В этом процессе участвуют не только клеточные элементы соединительной ткани, но и факторы, продуцируемые эндотелием [7]. При оперативных вмешательствах заполнение тканевого дефекта осуществляется грануляционной тканью, необходимым условием роста которой является развитие сети капилляров из эндотелиальных клеток (ангиогенез).
...
08 04 2026 10:45:36
Статья в формате PDF
117 KB...
06 04 2026 17:38:42
Статья в формате PDF
138 KB...
05 04 2026 8:29:46
Статья в формате PDF
464 KB...
04 04 2026 16:16:51
Статья в формате PDF
130 KB...
02 04 2026 11:58:48
Статья в формате PDF
125 KB...
01 04 2026 1:31:50
Статья в формате PDF
111 KB...
31 03 2026 6:27:30
Исследованы вопросы влияния давления, относительной влажности и температуры атмосферы на давление воздуха в шине 175/70R13 легкового автомобиля ВАЗ на основании данных ГУ «ВНИИГМИ-МЦД» по постам (станциям) о температуре воздуха, относительной влажности и атмосферном давлении на уровне станции по природно – климатическим поясам России. Вопросы влияния климатических хаpaктеристик на давление в автомобильных шинах рассмотрены для летнего периода, который является наиболее нагруженным в году периодом в плане эксплуатации автомобиля. Исследования выполнены методом случайной выборки с использованием данных срочных наблюдений по постам Федеральной службы по гидрометеорологии и мониторингу окружающей среды. Изменения давления в шине в течение рабочей смены значительно влияют на управляемость, надежность и экономическую эффективность эксплуатации автотрaнcпорта.
...
30 03 2026 4:22:43
Статья в формате PDF
269 KB...
29 03 2026 19:12:27
Статья в формате PDF
104 KB...
28 03 2026 21:42:28
Статья в формате PDF
285 KB...
27 03 2026 8:58:54
Статья в формате PDF
132 KB...
26 03 2026 8:15:28
Обсуждается проблема формирования структурных модулей, которые предназначены для конструирования невырожденных модулярных 3D структур кристаллов.
...
25 03 2026 11:52:22
Приведены методы ранжирования и рангового моделирования гидрологических параметров у множества крупных рек Земли по примеру статистических данных из учебника.
...
24 03 2026 22:34:45
Статья в формате PDF
100 KB...
23 03 2026 11:19:56
Статья в формате PDF
312 KB...
22 03 2026 16:53:48
Статья в формате PDF
127 KB...
21 03 2026 11:31:27
Статья в формате PDF
103 KB...
20 03 2026 13:14:34
Статья в формате PDF
113 KB...
17 03 2026 17:28:37
Статья в формате PDF
120 KB...
16 03 2026 11:48:53
Статья в формате PDF
125 KB...
15 03 2026 8:36:22
Статья в формате PDF
254 KB...
14 03 2026 1:51:24
Статья в формате PDF
112 KB...
13 03 2026 12:19:20
Статья в формате PDF
102 KB...
12 03 2026 23:41:38
Статья в формате PDF
124 KB...
11 03 2026 6:13:22
Статья в формате PDF
259 KB...
10 03 2026 11:14:33
Морфогенез лимфатической системы является результатом взаимодействия сосудов разного типа, растущих неравномерно. Его формы меняются так же, как строение и топография сосудов, их сочетания в связи с органогенезом. Поэтому морфогенез лимфатической системы протекает как процесс рекомбинации артерий и вен, а затем и лимфатических сосудов, служит проявлением самодифференциации сердечно-сосудистой системы, когда ее части вступают в повторное взаимодействие, в т.ч. и после их трaнcформации.
...
09 03 2026 13:22:32
Алтайский край разнообразен по рельефу, климату и почвам. Включает 5 природных зон – от сухой степи до увлажнённых предгорий. Гречиха посевная выращивается на всей территории края, однако её посевы наиболее продуктивны в условиях лесостепи, что связано с природными ресурсами и развитым пчеловодством. Применение зонального агротехнического комплекса в лесостепи позволяет получать высокий урожай зерна (1,5–2,0 т/га).
...
08 03 2026 18:28:17
Статья в формате PDF
141 KB...
07 03 2026 20:18:41
Статья в формате PDF
267 KB...
06 03 2026 11:17:26
Статья в формате PDF
122 KB...
04 03 2026 9:49:30
Еще:
Поддержать себя -1 :: Поддержать себя -2 :: Поддержать себя -3 :: Поддержать себя -4 :: Поддержать себя -5 :: Поддержать себя -6 :: Поддержать себя -7 :: Поддержать себя -8 :: Поддержать себя -9 :: Поддержать себя -10 :: Поддержать себя -11 :: Поддержать себя -12 :: Поддержать себя -13 :: Поддержать себя -14 :: Поддержать себя -15 :: Поддержать себя -16 :: Поддержать себя -17 :: Поддержать себя -18 :: Поддержать себя -19 :: Поддержать себя -20 :: Поддержать себя -21 :: Поддержать себя -22 :: Поддержать себя -23 :: Поддержать себя -24 :: Поддержать себя -25 :: Поддержать себя -26 :: Поддержать себя -27 :: Поддержать себя -28 :: Поддержать себя -29 :: Поддержать себя -30 :: Поддержать себя -31 :: Поддержать себя -32 :: Поддержать себя -33 :: Поддержать себя -34 :: Поддержать себя -35 :: Поддержать себя -36 :: Поддержать себя -37 :: Поддержать себя -38 ::