Вулканічні ударні хвилі в домашніх барометрах: раннє попередження на основі даних 29 334 метеостанцій
Як метеостанції Netatmo, фізика хвиль Ламба та штучні інтелектуальні агенти визначали кратери, оцінювали потужність вибухів та надавали попередження за кілька хвилин без хибних сигналів.
Джерела: Hearing the Mountain’s Roar: How Antigravity CLI’s AI Agents & IoT Data Track Volcanic Shockwaves
Коли вулкан вивергається експлозивно, він випромінює хвилі Ламба – акустично-гравітаційні хвилі, що застрягають у найнижчому шарі атмосфери та поширюються горизонтально по всьому світу зі швидкістю звуку. Дослідження, засноване на 29 334 реальних записах спостережень від побутових метеостанцій у Японії, показує, що звичайні барометри на балконі можуть фіксувати ці хвилі. Ця система дозволяє визначати швидкість хвиль, залежну від температури, знаходити напрямок неконтрольованого кратера, оцінювати енергію вибуху та надавати містам, розташованим нижче за течією, попередження тривалістю від 2,5 до 15 хвилин, без хибних сигналів під час сильних штормів.
У дослідженні використовуються два приклади вивержень 2018 року як основа: вибух, спричинений парою (фреатичний вибух) на горі Кусацу-Шіране, та магматичне виверження на горі Сінмоедаке. Результати представлені як один із академічних поглядів серед багатьох щодо інфразвуку вулканів та динаміки атмосферних хвиль, що має на меті заохотити до більшої кількості емпіричних досліджень та дискусій у геофізичній спільноті.
Проблема обмеженої кількості датчиків
Професійний моніторинг вулканів ґрунтується на мікрофонах, що фіксують інфразвук, та сейсмографах, якими керують урядові агентства та мережі спостереження за забороною ядерних випробувань. Ці станції розташовані на відстані від 100 до 500 кілометрів одна від одної. Така відстань створює великі сліпі зони, спричиняє серйозні проблеми з просторовою резолюцією та піддає окремі датчики впливу вітрових турбулентностей.
Датчики IoT для споживачів пропонують спосіб подолати цю проблему. Тисячі житлових станцій Netatmo по всій Японії можна розглядати як одну акустичну антену розміром з континент. 29 334 записи спостережень з дослідження свідчать про те, що датчики тиску, встановлені на балконах квартир, можуть попередити про наближення ударної хвилі та її попелу ще до того, як вони досягнуть міських центрів.
Процес дослідження
Усій проєкті потрібно отримувати дані з домашніх датчиків у такій формі, щоб до них можна було застосувати закони фізики. Одна з ілюстрацій ілюструє цю ідею: виверження викидає двовимірну циліндричну хвилю Ламба, яка рухається по землі зі швидкістю 1 134 км/год; тисячі житлових станцій у міських районах перехоплюють її, а системи аналізу даних розраховують потужність вибуху в тоннах ТНТ та виставляють таймер попередження від 2,5 до 15 хвилин.
Постійне отримання даних
Ця система — це автоматизований хмарний конвеєр, який працює без нагляду. Тригер Google Apps Script без сервера запускається кожні 20 хвилин та викликає API Netatmo Connect для отримання даних про барометричний тиск, температуру на вулиці та відносну вологість з станцій, розташованих у японських містах та біля підніжжя гір. Отримані дані індексуються, структуруються та архівуються у центральному сховищі Google Drive.
Чотири етапи роботи
Оскільки дані накопичуються безперервно, проект проходив у чотирьох етапах, які ведуть від сирого архіву до перевіреної теорії.
- Пряме отримання даних. Скрипт функціонує як фоновий механізм, перевіряючи дані кожні 20 хвилин для створення безперервного, довготривалого метеорологічного архіву.
- Картування історичних подій. Минулі виверження аналізуються за допомогою цього архіву шляхом визначення 24-годинних інтервалів від 12 годин до початку до 12 годин після нього як для виверження пари в Кусацу-Шіране, так і для виверження Шінмоедаке. Дні спокою та штормових умов слугують базовими показниками.
- Дослідження аномалій до та після виверження. Аналізуються серії даних атмосфери до та після кожного виверження з метою виявлення двох аспектів: чи з’являються слабкі аномалії гідротермального дегазування перед розривом, та як ударна хвиля проходить через міські райони після цього.
- Теорія та перевірка в середовищі CLI. Для тестування в Gemini, який працює в межах Antigravity CLI, були використані уніфіковані математичні моделі — ентальпія Болтона, двовимірні рівняння хвиль Ламба для стиснутих середовищ, метод формування променів на основі швидкості поширення сигналів та інтегрування енергетичного потоку у циліндричних координатах — щоб перевірити, чи можуть системи попередження та діагностика кратерів функціонувати без жодних хибних сигналів.
Фізика через чотири повсякденні приклади
Як пристрій вартістю $150, який не калібрувався, розташований на балконі, може фіксувати вибухи, що відбуваються за сотні кілометрів? Чотири знайомі аналогії пояснюють це.
Блискавка та грім
Спалах світла настає миттєво, а грім лунає через кілька секунд, оскільки світло рухається зі швидкістю 300 000 км/с, а звук — приблизно 340 м/с. Система попередження використовує саме цю різницю. Хвиля удару від виверження рухається по землі зі швидкістю близько 315 м/с, тобто приблизно 18,9 км на хвилину. Станції, розташовані на відстані від 10 до 30 км від кратера, відчувають початковий тисковий імпульс протягом 30–90 секунд та надсилають свої дані до хмари зі швидкістю світла.
Очевидний заперечення: якщо кожна станція надсилає дані пакетами кожні 20 хвилин, як може існувати затримка у кілька хвилин? На це є два пояснення.
- Час подорожі залежно від відстані. Хвилі потрібен час, щоб пройти певну територію. Затримка зростає з відстанню та є неминучою: громади, розташовані на відстані 50 км, чекають близько 2,6 хвилин, міста на відстані 100 км — близько 5,3, Токіо (180 км від Кусацу-Шіране) — близько 9,5, а транспортний коридор у Фукуокі (280 км від Сінмоедаке) — майже 15 хвилин, а саме 14,8 хвилин. Датчики, розташовані біля кратера, досягають хмари менш ніж за хвилину, тому сповіщення надходять значно раніше, ніж хвиля.
- Послідовне надсилання даних. Місцеві станції не передають дані одночасно, оскільки їхні 20-хвилинні цикли розкидані в часі. У щільному скупченні з 100 станцій новий пакет досягає хмари в середньому кожні 12 секунд, що перетворює окремі 20-хвилинні оновлення на майже безперервний потік даних. Вебхуки, активовані зміною умов, можуть забезпечувати передачу даних з інтервалом у кілька секунд, коли з’являється різкий градієнт тиску.
Хвилі на ставку
Вибух у відкритому просторі поширює енергію у трьох вимірах, а сигнал швидко слабшає, зменшуючись пропорційно квадрату відстані. Хвилі Ламба утримуються біля поверхні завдяки силі тяжіння та різниці щільностей повітря, тому вони поширюються у вигляді плоского, двовимірного кільця. Енергія не може виходити вгору, що спричиняє значно повільніше занепадання імпульсу — пропорційно до відстані, а не до квадрата цієї відстані. Цей ефект провідника зберігає сигнал цілими сотнями кілометрів, тому некалібрований барометр у центрі Токіо може фіксувати різкий мікроімпульс від кратера, розташованого на відстані 180 км.
Бінауральне слух на відстані 150 кілометрів
Люди можуть заплющити очі та вказати напрямок звуку, оскільки мозок порівнює часи прибуття сигналу до лівого та правого вуха з точністю до мікросекунд. Ця система працює аналогічно на всій території. Порівнюючи різниці у часі прибуття сигналу з 25 до 88 станцій, розташованих уздовж лінії довжиною 150 км, пристрій типу „delay-and-sum beamformer“, який працює у просторі повільності, вирівнює записані дані для кожного азимуту по всьому колу. Кут, при якому хвилеподібні форми посилюють одна одну, визначає місцезнаходження джерела звуку з точністю до півградуса, без необхідності попереднього знання його координат.
„Шлюз швидкості“ та тест хору
Місцевий вітер є основною проблемою для обладнання. Згідно з принципом Бернуллі, пориви вітру на відкритому балконі спричиняють різкі місцеві зниження тиску, і одна станція може зафіксувати сотні хибних сигналів у штормову погоду. Два фізичні фактори допомагають усунути ці проблеми.
- Швидкісний бар’єр. Пориви вітру супроводжують погодну фронтальну структуру зі швидкістю від 10 до 25 м/с, або від 36 до 90 км/год. Справжня вулканічна хвиля Ламба рухається зі швидкістю звуку – від 300 до 330 м/с, що дорівнює майже 1 134 км/год. Усе, що повільніше за звук, відразу відкидається.
- Тест регіонального хору. Порив вітру – це локальна, некорелювана подія, яка впливає на один будинок за раз. Вулканічна хвиля шоку – це кохерентна хвильова структура, яка по черзі зачіпає десятки будинків у всій області.
Вимога до фазової кохерентності між багатьма станціями зі значенням принаймні 0,70 дозволила виключити всі 378 сплесків швидкості вітру, зафіксованих під час сильних бурь. Рівень хибних позитивних результатів становив 0,0%.
Приклад обчислень часу настання
Цифри, які лежать в основі попередження, легко перевірити. При швидкості 315 м/с хвиля пройде 18,9 км за хвилину, тож місто, розташоване на відстані 50 км від кратера, матиме приблизно 2,6 хвилини між першим виявленням та прибуттям хвилі, а місто на відстані 280 км — майже 15 хвилин. Для ефективного раннього попередження важливо, щоб виявлення відбувалося у близькій зоні — на відстані 10–30 км, де імпульс фіксується протягом 30–90 секунд, а також щоб хмарні сервери отримували його протягом менш ніж хвилини. Кожен додатковий кілометр відстані між кратером та містом збільшує час корисного попередження приблизно на три секунди. Діапазон від 2,5 до 15 хвилин, наведений для цієї системи, — це просто результат застосування цих геометричних формул до найближчих та найвіддаленіших населених пунктів, досліджених у проекті.
Сама геометрія пояснює, чому щільне покриття має більше значення, ніж якість датчика. Один дешевий барометр створює шум, проте десятки таких пристроїв, які узгоджуються щодо часу та напрямку, створюють сигнал, якого не зміг би дати жоден окремий інструмент. Пристрої для споживачів не калібруються, тому значення абсолютного тиску є ненадійними, але форма та час прибуття чіткого імпульсу зберігаються, і саме часи прибуття є тим, що потрібно для формування променів та зворотного обчислення швидкості.
Керування ШІ за допомогою фізичної логіки
Генеративний ШІ добре справляється з виявленням закономірностей у високовимірних даних. Однак, якщо його сліпо застосовувати до шумних метеорологічних серій, неконтрольована модель може зафіксувати хибні кореляції та створити артефакти без жодних підстав. Щоб перетворити необроблені результати моделі на перевірені відкриття, потрібен гібридний підхід: агенти для пошуку закономірностей, пов’язані безпосередньо з теорією перших принципів та логікою конкретної галузі.
Основна ідея проста. Якщо живити модель не лише числами, а й фізичними законами, термодинамічними рівняннями та геометричними обмеженнями, це надає їй структурованого підходу, що підвищує шанси на виявлення справжніх закономірностей та видалення шуму. Та сама філософія застосовувалась у інших проектах з прогнозування катастроф. Для землетрусів використовувалися рівняння зв’язку літосфери та атмосфери та тензори сейсмічного моменту для виявлення іонізації радону та підвищення ентальпії на рівні від 5,96 до 15,47 кДж/год за 5,96–15,47 годин до розриву. Для сильних дощів у містах застосовувалися 3D-рівняння Нав’є-Стокса для виявлення локалізованих процесів конвекції за 45 хвилин до їх початку. Для вулканів рівняння хвиль Ламба, параметри пари Болтона та обмеження, пов’язані зі швидкістю поширення хвиль, допомогли відокремити моделі ударних хвиль та напрямки кратерів від 29 334 шумних записів.
Шість етапів обробки
Завдяки цьому підходу, спрямованому людиною, цей процес перетворює первинні показники на фізичні величини, як показано на наступних етапах. Коротко кажучи, дані про тиск, температуру та вологість у житлових районах проходять через фільтрацію, зворотний обчислення швидкості, формування променя та інтеграцію енергії, що дає уявлення про розмір вибуху та процес охолодження попелу.
- Відстеження ентальпії перед виверженням. Температура, вологість та тиск перетворюються на еквівалентну потенційну температуру за формулою Болтона. Віднімання регіональних тенденцій виявляє слабкі теплові аномалії у вигляді пари на відстані кількох годин від вулкана перед виверженням — від 1,5 до 3,2 К.
- Фільтрація нульової фази. Серії даних про тиск проходять через двонаправлені фільтри типу Butterworth з частотним діапазоном від 0,0005 до 0,02 Гц, які видаляють 24-годинні припливні цикли, зберігаючи при цьому часи прибуття хвиль без спотворень.
- Інверсія часу проходження хвилі, залежна від температури. Замість постійної швидкості звуку модель інтегрує дані за місцевими температурними полями, враховуючи сповільнення поширення звуку на замерзлих гірських плато.
- Формування променя. Кожен регіональний сукупність станцій має свої дані, які зсуваються в часі та підсумовуються для кожного напрямку в мережі сповільнень на 360 градусів. Напрямок, у якому сума є найсильнішою та найбільш послідовною, виміряна за показником схожості, визначає задній азімут.
- Інтеграція циліндричного потоку енергії. Варіація тиску інтегрується по двовимірній циліндричній поверхні та масштабується за допомогою атмосферної висоти масштабу приблизно 8,4 км, що дає загальну акустичну енергію та її еквівалент у TNT.
- Охолодження від попелу. Чутливі датчики температури в передмістях виявляють сонячне захистування під пливучими хмарами попелу, що проявляється місцевими зниженнями температури на 4,10–6,13 градусів Цельсія.
Що показали 29 334 записи
Два протилежні за характером явища створюють ідеальні умови для тестування. Морозний зимовий паровий вибух та помірний весняний магматичний вибух створюють різні швидкості хвиль та дуже різну енергію, тому метод, який може працювати з обома, важко відкинути як випадковість. Записи походять з двох масштабних вивержень та двох періодів сильної погоди, які забезпечили контрольні дані для порівняння та підрахунку хибних сигналів.
Перевірка охоплювала 29 334 справжніх 20-хвилинних записів, отриманих під час двох масштабних вивержень та двох подій із сильною погодою.
- Швидкість хвилі, залежна від температури. У Кусацу-Шіране ударна хвиля рухалася зі швидкістю 304,38 ± 4,73 м/с. Швидкість звуку на рівні моря становить близько 340 м/с, але звук уповільнюється в холодному щільному повітрі. Того січневого ранку при температурі -5 градусів Цельсія в горах Гумма теорія прогнозувала швидкість 304,8 м/с. Некалібрована мережа вимірювань виявила уповільнення звуку з точністю 99,8% (R-квадрат дорівнює 0,980). У Сінмоедаке при температурі 12 градусів навесні швидкість становила 313,27 м/с, що на 98,5% відповідає теоретичним значенням.
- Визначення місцезнаходження невідомого кратера з відстані 150 км. Формування променя через станції на Кюсю дало координати кратера Сінмоедаке у напрямку 169,49 градусів, тоді як справжній кут становив 167,72 градуса, що є похибкою у 1,78 градуса.
Три способи використання
Перетворення домашніх метеостанцій на засоби отримання геофізичної інформації в реальному часі дозволяє забезпечити три види захисту, кожен з яких залежить від часу попередження, що створюється відстанню та швидкістю мережі.
Авіація та міська інфраструктура
Літаки та міста стикаються з однаковою загрозою, але в різні терміни. Пілотам потрібно за кілька секунд уникати зростаючих стовпів попелу, тоді як містам потрібні хвилини, щоб підготувати свої підземні та приміщені простори.
Вулканічний попіл розплавляється всередині реактивних двигунів при температурі понад 1 400 градусів Цельсія та перетворюється на рідке скло, що може призвести до повної несправності двигуна. Мережа, яка за кілька секунд оцінює обсяг попелу та місцезнаходження кратера, може попередити служби контролю повітряного руху та системи управління польотами про необхідність зміни маршруту літаків навколо цих стовпів. Міста, розташовані по вітру, мають від 5 до 15 хвилин на те, щоб закрити вентиляційні шахти метро, перевести системи опалення та вентиляції лікарень на внутрішню фільтрацію типу HEPA, призупинити будівельні роботи на вулиці та почати уповільнювати рух високошвидкісних поїздів.
Вулканічні цунамі
Вибухи підводних кальдер та обвалення схилів островів, як це сталося на Кракатоа у 1883 році та на Гунга Тонга у 2022 році, спричиняють смертоносні цунамі. Цунамі в глибоких водах океану рухаються зі швидкістю близько 200 м/с та у мілководді узбережжя їхня швидкість знижується до 30–50 м/с, тоді як атмосферні хвилі Ламба рухаються зі швидкістю 315 м/с та випереджають їх. Берегові барометри, які фіксують ці повітряні хвилі, можуть надати 15–45 хвилин на евакуацію перед прибуттям води, що спричиняє автоматичне закриття воріт морських стін та зупинку залізничного руху.
Локалізація сліпих кратерів
Цей випадок має найбільше значення для вулканів, які ніхто не стежить.
Найгіршою вулканічною небезпекою є раптове виверження з неактивної гори чи неконтрольованої бічної тріщини без наявності урядових пристроїв поблизу. Існують два способи боротьби з цим. Перед розривом пара, що витікає, створює невеликі локальні аномалії потенційної температури у межах від 1,5 до 3,2 К, а розв’язання допоміжного рівняння адвекції-дифузії дозволяє визначити джерело стовпа пари та передбачити момент майбутнього вивітрювання. Після вибуху незалежні регіональні системи в Канто, Нагої та Кюсю обчислюють лінії перетину з похибкою 0,48 градуса, що дозволяє визначити місцезнаходження активної кратери з точністю до 840 метрів протягом кількох хвилин.
Система спільних досліджень з ШІ
Основний дослідник-людина відповідає за контроль та схвалення, тоді як група спеціалізованих агентів займається отриманням даних, моделюванням та симуляцією, розробкою коду, перевіркою на фальсифікацію та штучним рецензуванням. Кожен виконує вузькоспеціалізовані обов’язки: один перевіряє плани за допомогою математичних моделей, інший оцінює теоретичну доцільність, ще один пише код, інші перевіряють результати, забезпечують цілісність даних та аналізують теоретичні аспекти, а останні двоє перевіряють дотримання стилістичних норм та форматування.
Розв’язування рівнянь механіки континууму, отримання записів, перевірка на фальсифікацію та підготовка рукопису у форматі LaTeX, готового до публікації, виконувалися за допомогою Antigravity CLI, який працює на базі Gemini разом із допоміжною дослідницькою платформою під назвою tanaike-lab. Робота проходила у п’яти етапах.
- Стратегічне керування з боку людини. Головний дослідник розпочав роботу з математично обмеженої дослідницької концепції, а не з нечіткого запиту, визначивши шість чітких вимог: створити детектор хвиль Ламба з мікробарометричними характеристиками з серії IoT для громадської науки; використовувати 24-годинні інтервали (плюс-мінус 12 годин) навколо двох вивержень 2018 року з базовими даними про шторми для перевірки гіпотез; застосовувати критерій еквівалентної потенційної температури Болтона для ізоляції теплового потоку перед виверженням; застосовувати двонапрямлений фільтр нульової фази для збереження точних часових позначок; обчислювати швидкості поширення за допомогою змінної в просторі температури; а також використовувати техніку формування променів повільності для визначення азимутів та показників ефективності на основі різниць у часі прибуття сигналів з кількох станцій.
- Планування аудиту та математичних тестувань. Спеціалізовані підагенти перевірили план на математичну послідовність, підтвердили коректність циліндричної моделі енергетичного потоку та сформулювали рівняння інтегралу шляху температури з урахуванням гірського рельєфу.
- Розробка потокових процесів із механізмами перевірки. Підагент, відповідальний за розробку коду, створив ізольовані Python-потокові процеси, кожен з яких мав автоматизовані механізми перевірки, що гарантували використання справжніх даних через OAuth2, фільтрацію та формування променів без використання синтетичних масивів даних.
- Емпіричний аналіз та обговорення між агентами. Агенти, які здійснювали аудит, обговорили результати для 29 334 записів, перевірили збіг значень 304,38 м/с, оцінили точність визначення напрямку кратера та підтвердили повне відхилення 378 сплесків швидкості вітру.
- Перевірка колегами та подання. Симульована панель оцінювання з п’ятьма осями, крок перевірки дотримання стилю та інструмент перевірки верстки проводили аудити перед поданням, перевіряли, чи всі посилання в тексті відповідають, компілювали код LaTeX та створювали ілюстрації, доступні для сприйняття з обмеженими можливостями зору, перш ніж людина надасть остаточний дозвіл.
На відміну від статичних скриптів агентів, ця платформа фіксує уроки під час роботи. Цей цикл є цілеспрямовано автоматичним: директиви від людини та уроки, перевірені під час виконання, фіксуються в реальному часі, одразу записуються у документацію навичок та повертаються до інструменту без необхідності їх ручної перевстановки. Керівні директиви та перевірені уроки записуються під час роботи, безпосередньо вводяться у файл SKILL.md, потім експортуються та перевстановлюються у локальний каталог плагінів CLI з негайною синхронізацією через Git, що дозволяє інструментарію постійно розвиватися.
Його описана концепція складається з трьох частин. Вона протиставляє некероване генерування, яке супроводжується галюцинаціями та поверхневим текстом, моделі синергії людини та ШІ, у якій людина визначає стратегічний напрямок, тоді як підагенти прискорюють процеси формулювання, виконання коду, аудиту та перевірки. Метою є усунення операційних труднощів, таких як виправлення помилок у скриптах, інтерполяція нерівномірних часових решіток, коригування даних для дизайну, універсального до кольору, та усунення помилок LaTeX, щоб дослідники могли зосередитися на високорівневому міркуванні. І він охоплює три сфери: природничі науки, такі як геофізичне моделювання, динаміка рідин, квантова хімія та матеріали; інженерія генеративного ШІ, що включає проектування запитів, пошук інформації за допомогою RAG, оркестрацію багатьох агентів, стиснення бюджету токенів та тестові завдання; а також екосистему Google, яка охоплює бібліотеки Apps Script та автоматизацію робочого простору у Drive, Sheets, Docs, Gmail та Forms із керуванням квотами.
Що це дає в сукупності
Пристрої для моніторингу погоди в житлових приміщеннях у поєднанні з ШІ, заснованим на законах фізики, можуть стати мережею раннього попередження на континентальному рівні: хвилі Ламба та енергія ударних хвиль, отримані з барометрів, створених користувачами, без жодних хибних сигналів порівняно з базовими даними штормів. Ці висновки є одним із багатьох підходів у фізиці вулканів, і продовжувана емпірична перевірка та відкриті наукові дискусії залишаються надзвичайно важливими. Поза межами геофізики, модель людського формування стратегій та фізичних обмежень, поки агенти ШІ виконують завдання з обчислень, програмування, аудиту та перевірки, виглядає як масштабований шаблон для адаптації до клімату, інженерії з багатьма агентами та автоматизації API Google.
Пов’язана література
- Виконання коду в браузері для AI-агентів: Playwright, MCP та безпечний виконавець — Чому агентам потрібен браузер як середовище виконання, як він поєднується з MCP та режимом коду, та як забезпечити його безпеку за допомогою профілів, обмежень за доменом та політик дій.
- Основи AI-агентів: мета, проектування запитів та вибір моделі — Як закласти основи для професійного AI-агента: визначити мету та показники успіху, написати запити з чіткими цілями та обмеженнями, а потім вибрати відповідну модель.