Форум программистов, компьютерный форум, киберфорум
anaschu
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск  

сукцессия 36. Делаю прогоны модели около порога бифуркации (превращения эндо микоризы в экто микоризу)

Запись от anaschu размещена 21.07.2026 в 12:57
Показов 647 Комментарии 0

СУКЦЕССИЯ МИКОРИЗ: КАК МАТЕМАТИКА ОБЪЯСНЯЕТ ПОДЗЕМНУЮ ВОЙНУ

Мы привыкли думать, что лес — это деревья. А грибы — это то, что мы собираем в корзинку осенью. Но если заглянуть под землю, окажется, что настоящий хозяин леса — это грибница. И она ведёт жестокую, многолетнюю войну за ресурсы. Мы с коллегами построили математическую модель этой войны и нашли в ней несколько неожиданных бифуркаций, о которых я хочу рассказать.

Про ОДУ, грибы и деревья

Всё, что я буду описывать — это система из примерно 30 обыкновенных дифференциальных уравнений (ОДУ). Это значит, что мы следим за 30 переменными во времени: массы разных частей грибов и деревьев, количество фосфора и азота в почве, кислотность.

Есть два главных игрока:

АМ-гриб (арбускулярный) — мирный симбионт. Он живёт внутри корней широколиственных деревьев, помогает им усваивать фосфор, а взамен получает сахара. Работает в нейтральной или слабокислой среде (pH ~ 6.0).

ЕСМ-гриб (эктомикоризный) — агрессивный захватчик. Он оплетает корни снаружи, но что важнее — умеет закислять почву и добывать фосфор из труднодоступных минералов. Работает в кислой среде (pH ~ 4.8).

Задача модели: понять, при каких условиях мирный АМ-лес сменяется агрессивным ЕСМ-лесом. Эта смена называется сукцессией.

Первый сюрприз: граница

В теории Гуца говорилось, что смена происходит скачком. Нужно, чтобы два условия совпали: (1) почва достаточно закислена, (2) нового гриба пришло достаточно много. Без выполнения обоих — перехода не будет.

Мы проверили это на модели.

Оказалось, что у нас есть два независимых триггера:

Кислотный триггер: когда pH падает ниже 5.6, ЕСМ-гриб просыпается.

Бифуркационный параметр: скорость градации подроста (k_conifer_graduate). Это скорость, с которой молодые хвойные деревца «выпускаются» в устойчивые взрослые островки.

При фиксированных прочих условиях мы нашли резкую границу:

 
[tr][td]k[/td][td]Исход[/td][/tr]
[tr][td]< 0.035[/td][td]Вымирание ЕСМ[/td][/tr]
[tr][td]> 0.040[/td][td]Прорыв ЕСМ[/td][/tr]

Скачок с 3.9 г до 80 г биомассы при изменении параметра всего на 14% — это и есть настоящая бифуркация.

Второй сюрприз: ансамбль

Мы поняли, что в реальном лесу параметры не фиксированы. Каждый участок леса — разный. Поэтому мы провели ансамблевый эксперимент: 50 случайных реализаций, в которых три параметра одновременно варьировались на 10–25% вокруг базовых значений.

При одновременном варьировании всех условий, k_conifer_graduate перестал быть главным фактором. Его средние значения в группе «прорыв» и «вымирание» почти совпали:

 
[tr][td]Параметр[/td][td]BREAK[/td][td]FAIL[/td][/tr]
[tr][td]k_conifer_graduate[/td][td]0.0376[/td][td]0.0378[/td][/tr]
[tr][td]PS_base_conifer[/td][td]0.0389[/td][td]0.0319[/td][/tr]

Третий сюрприз: влажность

Но два исключения заставили нас копнуть глубже. При низком PS_base_conifer (0.0335) часть прогонов всё равно давала прорыв. Что их отличало?

Порог влажности.

Мы зафиксировали PS_base_conifer = 0.0335 и k_conifer_graduate = 0.0375, а варьировали только moisture_threshold — параметр, определяющий, при какой влажности почвы наступает засуха для подроста.

И получили чёткую картину:

 
[tr][td]moisture_threshold[/td][td]Исход[/td][/tr]
[tr][td]< 0.
235]< moisture_threshold </td><td>Исход</td></tr>
[tr][td]< 0.235[/td][td]BREAK (5/5)[/td][/tr]
[tr][td]> 0.235[/td][td]FAIL (0/7)[/td][/tr]

Вывод: если moisture_threshold низкий (т.е. засуха наступает поздно), подрост успевает набрать критическую массу даже при слабом фотосинтезе. Если порог высокий (засуха наступает рано) — подрост не успевает, и ЕСМ вымирает.

Что это даёт?

Прорыв ЕСМ возможен, если:
* PS_base_conifer > 0.0355 (хороший фотосинтез), ИЛИ
* PS_base_conifer < 0.0355 И moisture_threshold < 0.235 (слабый фотосинтез, но мягкая засуха)

Иначе — вымирание.

С практической точки зрения это означает, что для управления сукцессией нужно не просто увеличивать приток ЕСМ-грибов, а заботиться о том, чтобы хвойный подрост мог накопить биомассу в первые, самые уязвимые годы. Осветление, снижение конкуренции за воду — это не метафоры, это рычаги управления фазовым переходом.

Вместо заключения

В этой модели по-прежнему остаются открытые вопросы. Что делает защелачивание АМ-грибов? Как влияет C:N:P как самостоятельный триггер? Как работает ансамбль в стохастическом режиме (не 50, а 500 прогонов)?

Но главное мы уже знаем: подземный лес — это не просто сеть корней, а сложная динамическая система с бифуркациями. И математика позволяет увидеть её правила.

Исходный код доступен в репозитории. Если вам интересно — можете запустить сами.
Размещено в Без категории
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Всего комментариев 0
Комментарии
 
Новые блоги и статьи
Установка MinGW GCC 16.2 и CMake
8Observer8 10.08.2026
VK Видео: https:/ / vkvideo. ru/ video-240781534_456239017 YouTube: eY5-5PyI9NM Текстовая версия
Неделя из жизни имитационной модели склада: мои кривые руки растут, откуда надо
anaschu 10.08.2026
Неделя из жизни имитационной модели склада: как я почти написал неправильную логику и что с этим делать Работаю сейчас над учебно-рабочим проектом: строю в AnyLogic имитационную модель процессов. . .
Калькулятор для расчета родства
russiannick 07.08.2026
1. Задача: Создать калькулятор для расчета родства. Родственных связей существует 8 ступеней, такие как: p - отец P - мать q - муж Q - жена b - брат B - сестра s - сын S - дочь
Мир по моей воле
kumehtar 07.08.2026
Когда-то кажется, что всё просто. Ты весь такой светлый. Причиняешь добро. Борешься за справедливость в этом тёмном мире. Потом начинаешь замечать одну неприятную вещь. Почти каждый хороший. . .
Кредитный калькулятор
Maks 05.08.2026
Решение задачи по прикладной информатике средствами 1С. Задача: Напишите приложение-калькулятор, которое помогает рассчитывать параметры кредита для аннуитетного и дифференцированного видов. . .
У нас сейчас поговорку "Опять 25" нужно переделать на "Опять +35".
kumehtar 04.08.2026
С ностальгией вспоминаю времена моего детства, когда у нас и правда +25 - была максимальная температура летом. Раньше +25 °C реально казались вершиной жары, когда можно было весь день пропадать на. . .
Как ИИ начал спорить и врать (возможно почуяв опасность для себя от индустрии - уход от электроники).
Hrethgir 04.08.2026
Недельный диалог, на фоне событий с НПЗ. Да, из спирта можно получать бензин, и это не сложно. Но потом в схеме я решил избавиться от насоса, при этом полностью сделав контроль подачи спирта в. . .
Термопринтер QR701
Argus19 03.08.2026
Термопринтер QR701 Купил два термопринтера QR701. На сэлф-тесте написано: Language: PC936 (GB18030). Что означает, что принтеры могут печатать только латиницу и китайские иероглифы. Так же. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru