Про хвилю популяції

Ця симуляція моделює просторове поширення інвазивного виду у 2D ландшафті за допомогою реакційно-дифузійного рівняння Фішера-КПП. Починаючи з невеликої центральної колонії, спостерігайте, як хвиля популяції розширюється назовні з передбачуваною швидкістю, що визначається лише двома параметрами: дифузією D та швидкістю росту r.

Рівняння Фішера-КПП

Рівняння, що керує процесом: ∂u/∂t = D∇²u + ru(1−u), де u(x,y,t) — локальна щільність популяції, нормована до [0,1]. Перший доданок описує випадкове розселення (дифузія), другий — логістичний ріст популяції. Фішер довів у 1937 році, що будь-яка локалізована початкова популяція поширюється як хвиля, що рухається, з мінімальною швидкістю v = 2√(Dr).

Ефект Аллі

Увімкнення ефекту Аллі змінює доданок росту на ru(u−a)(1−u). Коли локальна щільність падає нижче порогу a, швидкість росту стає від'ємною — популяція скорочується замість зростання. Це створює мінімальну життєздатну щільність популяції для інвазії, уповільнюючи або зупиняючи хвильовий фронт.

Чисельний метод

Сітка 150×150 еволюціонує за допомогою явних кінцевих різниць з граничними умовами нульового потоку (Неймана). Крок за часом автоматично обмежується умовою КФЛ dt ≤ dx²/(4D) для чисельної стійкості. Колір кодує щільність: синій (u=0, порожньо) → зелений (u=0.5, половина ємності) → червоний (u=1, насичено). Нижня смуга показує загальну популяцію в часі.

Часті запитання

Що таке рівняння Фішера-КПП?
Рівняння Фішера-КПП — це диференціальне рівняння з частинними похідними, що описує просторове поширення популяції у середовищі. Запропоноване незалежно Рональдом Фішером та Колмогоровим, Петровським і Піскуновим у 1937 році, воно поєднує дифузію (поширення особин через випадковий рух) з логістичним зростанням (збільшення популяції, обмежене ємністю середовища). Рівняння: ∂u/∂t = D∇²u + ru(1−u), де u — локальна щільність популяції (від 0 до 1), D — дифузія, r — внутрішня швидкість росту.
Що означає формула швидкості хвилі v=2√(Dr)?
Фішер довів, що для рівняння КПП з будь-якою локалізованою початковою умовою фронт популяції рухається з мінімальною швидкістю v = 2√(Dr). Цей елегантний результат пов'язує два біологічні параметри: D (як швидко розселяються особини) та r (як швидко зростає популяція). Більша швидкість розселення або більша швидкість росту збільшують швидкість інвазії пропорційно — подвоєння D або r збільшує швидкість хвилі в √2 разів. Хвильовий фронт має характерну сигмоїдну форму, переходячи від u=1 (заселений ареал) ззаду до u=0 (порожній ареал) спереду.
Що таке ефект Аллі і як він уповільнює інвазію?
Ефект Аллі виникає, коли особини при малій щільності популяції мають нижчу пристосованість, ніж при помірній — наприклад, через труднощі з пошуком партнерів, зниження ефекту розбавлення хижаками або провал кооперативного пошуку їжі при малій щільності. У модифікованому рівнянні ∂u/∂t = D∇²u + ru(u−a)(1−u) поріг a > 0 створює штраф за швидкість росту нижче щільності a. Популяція з початковою щільністю нижче a скорочуватиметься, а не зростатиме, унеможливлюючи встановлення колонії з однієї-двох особин. Це різко уповільнює або зупиняє хвильові фронти інвазії.
Як чисельно симулюється 2D дифузія?
Симуляція використовує метод кінцевих різниць із явною схемою Ейлера. Лапласіан ∇²u апроксимується за допомогою п'ятиточкового шаблону: ∇²u ≈ (u[i+1,j] + u[i−1,j] + u[i,j+1] + u[i,j−1] − 4u[i,j]) / Δx². При кроці сітки Δx=1 та кроці за часом Δt=0.1 схема стійка за умови Δt ≤ Δx²/(4D) (умова КФЛ для 2D дифузії). Граничні умови — нульовий потік (Неймана): градієнт щільності нульовий на краях сітки, тому популяція не може покинути область симуляції.
Чому хвильовий фронт має певну ширину?
Розв'язок хвилі, що рухається, для рівняння Фішера-КПП має ширину фронту, що масштабується як √(D/r): висока дифузія розширює фронт, тоді як висока швидкість росту утримує його вузьким. Фізично ширина фронту представляє перехідну зону між усталеною популяцією (u≈1) та незаселеним ареалом (u≈0). У цій зоні особини активно ростуть і розселяються. Профіль фронту приблизно є логістичною сигмоїдою, а його точна ширина залежить від швидкості хвилі відносно мінімальної швидкості 2√(Dr).
Що таке неоднорідний ландшафт і як він впливає на поширення?
Неоднорідний ландшафт має просторово змінну дифузію D(x,y) — одні ділянки дозволяють швидкий рух (відкриті поля), інші — повільний (густий ліс, річки). У симуляції режим неоднорідності множить D на випадковий коефіцієнт від 0.5 до 1.5 у кожній комірці сітки. Це спричиняє швидке просування фронту інвазії через ділянки з високим D і затримку на ділянках з низьким D, формуючи нерівний, «пальцеподібний» фронт замість гладкої кругової хвилі. Це відповідає польовим спостереженням за інвазивними видами, що поширюються швидше вздовж доріг, річок або коридорів середовища існування.
Які реальні інвазії описала модель Фішера-КПП?
Модель Фішера-КПП була застосована до поширення мускусних щурів у Європі (завезені 1905 р., ~1 км/рік, підібрано D≈51 км²/рік), поширення каштанової гнилі в Північній Америці, поширення ВІЛ у лімфовузлах та історичного поширення землеробства в неолітичній Європі (~1 км/рік протягом тисячоліть). В усіх випадках ключове передбачення — постійна швидкість хвилі пропорційна √(Dr) — добре узгоджується з польовими даними D та r.
Як симуляція вимірює швидкість хвилі?
Симуляція відстежує радіус, при якому щільність популяції дорівнює 0.5 (контур напівмаксимуму), скануючи від центру назовні. При кожному кадрі зміна цього радіуса, поділена на минулий час, дає миттєву швидкість хвилі в комірках за одиницю часу. У режимі Фішера-КПП виміряна швидкість збігається з теоретичною 2√(Dr) після короткого перехідного періоду, поки фронт загострюється від початкового гладкого профілю колонії до справжньої форми хвилі, що рухається.
Чи можна зупинити або розвернути інвазію?
З увімкненим ефектом Аллі та достатньо високим порогом a хвильовий фронт інвазії може зупинитися — щільність на передньому краї падає нижче a і локальна швидкість росту стає від'ємною, спричиняючи відступ. Ця бістабільність (або повна інвазія, або вимирання) є ключовою рисою моделі сильного ефекту Аллі. Без ефекту Аллі (стандартний Фішер-КПП) будь-яка позитивна початкова щільність врешті поширюється по всьому ареалу — тривіальний стан u=0 нестійкий. Стратегії контролю, такі як біологічний контроль (зменшення r) або створення бар'єрів (локальне зменшення D), можуть уповільнити, але зазвичай не зупинити хвилю Фішера-КПП.
Які інші явища описує реакційно-дифузійне моделювання?
Реакційно-дифузійні рівняння, як-от модель Фішера-КПП, описують широке коло явищ, окрім біологічної інвазії. Поширення хімічних хвиль у реакції Бєлоусова-Жаботинського, інвазія пухлин у тканинах, кальцієві хвилі у клітинах, формування візерунків на шкурі тварин (нестійкість Тюрінга) та поширення епідемій у мережах — усі ці явища підпорядковуються реакційно-дифузійній динаміці. Ключовий висновок полягає в тому, що поєднання локальної реакції (ріст/загасання) та просторової дифузії породжує складні просторові патерни та хвилі, що поширюються, які принципово відрізняються від кожного процесу окремо.

Про цю симуляцію

Ця симуляція моделює просторове поширення інвазивного виду за допомогою реакційно-дифузійного рівняння Фішера-КПП, ∂u/∂t = D∇²u + ru(1−u), розв'язаного чисельно на сітці 150×150. Невелика центральна колонія розширюється назовні у вигляді хвилі, що біжить, чия швидкість збігається до теоретичного мінімуму v = 2√(Dr), а необов'язковий доданок ефекту Аллі може сповільнити або зупинити інвазію нижче порогу щільності.

🔬 Що демонструє

Щільність населення u(x,y), що еволюціонує на сітці 150×150 за допомогою явних скінченних різниць із п'ятиточковим шаблоном Лапласіана та нульовопотоковими (Неймана) межами. Синій позначає порожній ареал (u=0), зелений — половину місткості (u=0.5), а червоний — насичену колонію (u=1). Смуга під сіткою показує загальне населення в часі, а живі показники порівнюють виміряну швидкість хвилі з теоретичним значенням 2√(Dr).

🎮 Як користуватися

Регулюйте коефіцієнт дифузії D і темп росту r, щоб змінити, наскільки швидко особини розселюються й розмножуються — спостерігайте, як теоретична швидкість хвилі оновлюється миттєво. Увімкніть ефект Аллі та підвищте поріг a, щоб побачити, як інвазія сповільнюється або повністю провалюється. Перемкніть неоднорідний ландшафт, щоб рандомізувати локальну дифузію на ±50%, отримуючи нерівний, розгалужений фронт замість гладкого кола. Скидання перезапускає з чистої центральної колонії; пауза/відтворення зупиняє анімацію.

💡 Чи знали ви?

Та сама формула v = 2√(Dr), що керує цією змодельованою колонією, була підігнана до реальних інвазій — ондатри, що поширювалися Європою після їх завезення у 1905 році, просувалися приблизно на 1 км/рік, що відповідає хвилі Фішера-КПП з D≈51 км²/рік. Та сама математика описує й поширення каштанового раку, ВІЛ у лімфатичних вузлах і навіть історичне поширення землеробства неолітичною Європою.

Поширені запитання

Які повзунки керують симуляцією і що вони змінюють?

Коефіцієнт дифузії D (0.1–3.0) задає, наскільки швидко особини розселюються за рахунок випадкового руху; темп росту r (0.1–2.0) задає логістичну швидкість розмноження; а поріг Аллі a (0–0.5) задає мінімальний штраф щільності, що використовується при увімкненому ефекті Аллі. І D, і r напряму впливають на теоретичну швидкість хвилі, показану як 2√(Dr), тож зміна будь-якого повзунка миттєво оновлює цей показник, а скидання відновлює свіжу колонію радіусом 8 клітин у центрі сітки 150×150.

Як симуляція фактично обчислює хвилю на кожному кадрі?

Кожен кадр виконує чотири явні кроки Ейлера рівняння Фішера-КПП. Лапласіан апроксимується п'ятиточковим шаблоном (сума чотирьох сусідів мінус чотири рази центральна клітина, поділена на квадрат кроку сітки Δx²=1), крок часу обмежений dt = min(0.1, 1/(4D)) для задоволення умови стійкості CFL двовимірної дифузії, а ріст використовує або ru(1−u), або, з увімкненим ефектом Аллі, ru(u−a)(1−u). Межі нульовопотокові, тож щільність не може вийти за межі сітки.

Чим виміряна швидкість хвилі відрізняється від показника «Теорія 2√(Dr)»?

Симуляція відстежує радіус, на якому щільність вперше падає до 0.5, скануючи назовні від центру, а потім ділить зміну цього радіусу на минулий змодельований час, щоб отримати виміряну швидкість, що оновлюється приблизно кожні півтори одиниці часу. Показник «Теорія 2√(Dr)» — це замкнена формула Фішера, обчислена безпосередньо з поточних значень повзунків D і r. На початку прогону виміряна швидкість нижча за теоретичну, поки фронт стає крутішим від початкового гладкого профілю; після короткого перехідного періоду обидва значення збігаються.

Що насправді змінює увімкнення ефекту Аллі в моделі?

З вимкненим перемикачем Аллі ріст є стандартним логістичним доданком ru(1−u), і будь-яка додатна щільність врешті-решт захоплює всю сітку. Увімкнення змінює ріст на ru(u−a)(1−u): щойно локальна щільність падає нижче порогу Аллі a, цей доданок стає від'ємним, тож ділянки нижче a зменшуються замість росту. Підвищення a повзунком робить ефективний поріг інвазії суворішим, і достатньо високе значення може повністю зупинити фронт або призвести до колапсу колонії замість поширення.

Що відбувається при увімкненні неоднорідного ландшафту?

Увімкнення будує карту дифузії для кожної клітини, де кожна з 150×150 клітин отримує D, помножений на випадковий множник від 0.5 до 1.5, згенерований один раз при скиданні. Під час обчислення симуляція використовує це локальне значення замість однорідного D усюди, де активний перемикач, тож хвильовий фронт просувається швидше через ділянки з високою дифузією і відстає в ділянках з низькою дифузією — створюючи той самий тип нерівного, розгалуженого поширення, що спостерігається, коли реальні інвазії йдуть уздовж доріг, річок чи лісових коридорів, а не по відкритій, однорідній місцевості.