Original: http://math.ucr.edu/home/baez/this.week.html
11 квітня 2012 року
Джон Баез
На цьому тижні я намагаюся щось нове: в тому числі модель клімату, яка працює на вашому веб-браузері!
Це не реалістична модель; ми тільки почали. Але деякі програмісти в команді Азимут проекту зацікавлені в створенні більш таких моделей, особливо Аллан Ерскін (хто зробив цей), Джим Статтард (який допоміг мені змусити його працювати), Глин Еджі і Стаффан Лілейгерен. Це може бути цікавий спосіб для нас, щоб дізнатися і пояснити фізику клімату. З достатньою кількістю цих моделей, ми б весь курс онлайн! Якщо ви хочете допомогти нам, будь ласка сказати привіт.
Аллан буде говорити більше про проблеми програмування пізніше. Але по-перше, велика загадка: як невеликі зміни в орбіті Землі призводять до великих змін в кліматі Землі? Як я вже говорив минулого разу, здається, важко зрозуміти льодовикові цикли останніх декількох мільйонів років, відповідаючи на це.
Чи існують механізми зворотного зв’язку, які можуть ампліфікувати невеликі зміни температури? Так. Ось кілька очевидних з них:
- Зворотний зв’язок водяної пари. Коли стає тепліше, вода випаровується, і повітря стає вологим. Але водяна пара є парниковим газом, який викликає додаткове потепління. І навпаки, коли земля остигає, повітря стає сухішою, так що парниковий ефект стає слабкішим, що має тенденцію охолоджувати речі вниз.
- Зворотний зв’язок альбедо льоду. Сніг і лід відбивають більше світла, ніж рідкі океани або грунт. Коли Земля прогрівається, танення снігу і льоду, так що Земля стає темніше, поглинає більше світла, і, як правило, щоб отримати отримати ще тепліше. І навпаки, коли Земля охолоджується, більше снігу і льоду форма, так що Земля стає світліше, поглинає менше світла, і прагне отримати ще крутіше.
- Зворотний зв’язок розчинності діоксиду вуглецю. Холодна вода може містити більше вуглекислого газу, ніж теплою водою: саме тому відкриття теплий може соди може бути настільки вибухонебезпечна. Так що, коли океани Землі розігріти, вони виділяють вуглекислий газ. Але вуглекислий газ є парниковим газом, який викликає додаткове потепління. І навпаки, коли oceaans остигати, вони поглинають більше двоокису вуглецю, так що парниковий ефект стає слабкішим, що має тенденцію охолоджувати речі вниз.
Звичайно, є і негативні зворотні зв’язки: в іншому випадку клімат буде вкрай нестабільні! Існують також складні зворотні зв’язки, чий загальний ефект важче оцінити:
- Зворотний зв’язок Планка. Більш гарячий світ випромінює більше тепла, яка охолоджує його вниз. Це велика негативний зворотний зв’язок, яка зберігає всі позитивні відгуки від створення Землі шалено гарячої або шалено холодної.
- Зворотний зв’язок хмар. Тепліша Земля має більше хмар, які відбивають більше світла, але і збільшують парниковий ефект.
- Зворотний зв’язок градієнта зміни метеорологічного елемента з висотою. Підвищений парниковий ефект змінює вертикальний профіль температури атмосфери, яка має наслідки своєї власної, але це працює по-іншому у полюсів і поблизу екватора.
Див. тиждень302, щоб отримати докладнішу інформацію про зворотні зв’язки і наскільки великі вони, ймовірно, будуть.
На вершині всіх цих тонкощів, будь запропоноване рішення загадки льодовикових циклів потрібно тримати кілька інших речей на увазі, теж:
- Дійсно хороша теорія пояснює, чому не тільки у нас є льодовикові цикли в даний час, але чому у нас не було їх раніше. Як я вже пояснював у тиждень317, вони почали близько 5 мільйонів років тому, стало набагато холодніше, близько 2 мільйонів років тому, і перейшли з циклу приблизно 41000 років до циклу приблизно 100000 років складе близько 1 млн років тому.
- Скажімо, ми пофантазувати колосальні великий механізм позитивний результат впливу який робить велику роботу з підтримки тепла землі, коли тепло і холод, коли холодно. Якщо цей ефект досить сильний, Земля може бути бістабільна: вона матиме два стійких стани, теплий і холодний. На жаль, якщо ефект буде занадто сильним, він не буде легким для Землі, щоб пхати туди і назад між цими двома станами!Класичним прикладом бістабільної системи є автоматичне включення скажімо для електричного світла. Коли світло на ньому залишається включеним; коли світло вимкнений він залишається вимкненим. При натисканні перемикача дуже м’яко, нічого не відбудеться. Але якщо натиснути на нього досить важко, він раптово вискочити з включеного на вимкнений, або навпаки.Якщо ми намагаємося змоделювати льодовикові цикли, використовуючи цю ідею, нам потрібно перемикач, щоб мати досить драматичний ефект, але як і раніше чуйно реагувати на досить ніжний дотик. Для цієї роботи нам потрібно досить позитивного зворотного зв’язку… але не дуже багато.(Ми могли б також спробувати іншу ідею: може бути, Земля тримає себе в своєму крижаному льодовикової стані, або його теплого міжльодовикового стану, використовуючи якийсь механізм, який поступово використовує щось. Потім, коли Земля проходить з цього матеріалу, то, що це, клімат може легко перевернути в інший стан.)
- Ми завжди повинні пам’ятати, що в хорошому наближенні, загальна кількість сонячного світла, що вражає Землю щороку не змінюється в міру зміни орбіти Землі в так званих “циклах Миланковича”, які, здається, викликають льодовикові періоди. Я пояснив, чому минулого разу. Які зміни не загальна кількість сонячного світла, але щось набагато тонше: кількість сонячного світла на певних широтах в певні пори року! Зокрема, Миланкович стверджував, і більшість вчених вважають, що Земля має тенденцію ставати холодно, коли є трохи сонячного світла удару далеко північні широти влітку.
З цих та інших причин, будь-яке рішення льодовикового періоду головоломки зобов’язаний бути тонким. Замість того, щоб занурюватися прямо в цю складну трясовину, давайте спробуємо щось набагато простіше. Давайте просто думати про те, як ефект альбедо льоду може, теоретично, зробити Землю бістабільною.
Щоб зробити це, давайте подивимося на дуже простій моделі в цій великій поки ще не опублікованої книги:
- Джеральд Р. Норт, Прості моделі глобального клімату.
Це нульвимірна модель балансу енергії, а це означає, що вона включає в себе тільки середня температура Землі, середня сонячна радіація приходить в, а середнє інфрачервоне випромінювання виходить.
Середньою температурою буде T, вимірюване за Цельсієм. Ми будемо вважати, що Земля випромінює потужність квадратний метр, що дорівнює
де A = 218 Вт/м2, а B = 1,90 Вт/м2 на кожен градус за Цельсієм. Це лінійне наближення береться з супутникових даних на нашій Землі. Насправді, потужність випромінювання зростає швидше, ніж лінійно зі збільшенням температури.
Ми будемо вважати, що Земля поглинає сонячну енергію енергії на квадратний метр, що дорівнює
Ось:
- Q – це середня інсоляція: тобто, кількість сонячної енергії на квадратний метр удару у верхній частині атмосфери Землі, усереднене за місцем і часу року. У реальності Q дорівнює приблизно 341,5 Вт/м2. Це одна четверта сонячної постійної величини, тобто сонячної енергії на квадратний метр, що б вдарити панель, що завис у просторі над атмосферою Землі та були звернені прямо на Сонце.Чому чверть? Це хороша геометрія головоломка: ми одного разу працювали над нею на блозі Азимут.
- c(T) – коальбедо: частка сонячної енергії, яка поглинається. Коальбедо залежить від температури; ми повинні сказати, як.
Зважаючи на все вищенаведене, маємо таке:
де C – теплова потужність Землі в джоулях на градус за квадратний метр. Звичайно, це кумедна річ, тому що теплова енергія зберігається не тільки на поверхні, але і в повітрі і/або води, а також деталі сильно варіюються в залежності від того, де ми знаходимося. Але якщо розглядати єдину планету з сухим повітрям і не океан, Північний каже, ми можемо приблизно взяти С, рівну приблизно половині теплоємності при постійному тиску колони сухого повітря над квадратний метр, а саме 5 мільйонів джоулів на градус Цельсія.
Найпростіше зробити, це знайти рівноважні рішення, тобто рішення, в яких dTdt=0, таким чином
Тепер C не має значення! Ми хотіли б вирішити для T в залежності від сонячного випромінювання Q, але це легше вирішити для Q як функції Т:
Щоб йти далі, нам потрібно вгадати деяку формулу для коальбедо c(T). Коальбедо, пам’ятайте, що це частина сонячного світла, яка поглинається, коли воно потрапляє на Землю. Це 1 мінус альбедо, яка представляє собою фракцію, яка отримує відображення. Ось невелика діаграма альбедо:
Якщо переплутати між альбедо і коальбедо, просто пам’ятайте: вугілля має високе коальбедо.
Так як ми намагаємося тримати речі дуже просто право не, а не моделювати природу у всій своїй славної складності, давайте просто скажемо, середнє альбедо Землі становить 0,65, коли дуже холодно, і є багато снігу. Отже, нехай
буде “icy” (льодове) коальбедо, добре для дуже низьких температур. Точно так же, скажімо, середнє альбедо падає до 0,3, коли його дуже жарко, і Земля темніше. Отже, нехай
буде “ice-free” (без льоду) коальбедо, добре для високих температур, коли Земля темніше.
Потім нам потрібна функція температури, який інтерполює між ci та cf. Спробуймо цю:
Якщо ви не шанувальник функції гіперболічного тангенса Тань, це може здатися страшно. Але не лякайтеся!
Функція 12(1+tanh(γT)) це просто функція, яка плавно переходить від 0 при низьких температурах до 1 при високих температурах. Це гарантує, що коальбедо знаходиться поблизу крижаного значення ci при низьких температурах, а також поблизу його вільної від льоду значення cf при високих температурах. Але найцікавіше тут γ, параметр, який говорить, як швидко піднімається coalbedo як Земля стає тепліше. Залежно від цього, ми отримаємо різні ефекти!
Функція c(T) зростає найшвидше при T=0, оскільки ось де tanh(γT) має найбільший нахил. Ми просто пощастило, що в градусах Цельсія T = 0 є точкою плавлення льоду, так що це робить трохи сенсу.
Тепер програмування магії Алана Ерскін в гру вступає! Ви можете рухати повзунок для регулювання параметрів γ до різних значень між 0 і 1:
Можна побачити, як коальбедо c(T) змінюється як функція температури T. На цьому графіку температура коливається від -50 °C до 50 °C; графік залежить від того, яке значення γ ви вибираєте за допомогою повзунка:
Ось інсоляція Q, потрібна, щоб отримати задану температуру Т від -50 ° C до 50 ° C:
Це найпростіше вирішити для Q з точки зору Т. Але це більш інтуїтивним, щоб перевернути цей графік і побачити, що рівноважні температури T дозволено для даного інсоляції Q від 200 до 500 Вт на квадратний метр:
Вражає те, що коли γ значно зростає, три різні температури сумісні з тим же самим обсягом сонячної радіації! Це означає, що Земля може бути гарячим, холодним або щось проміжне з’єднання, навіть якщо кількість сонячного світла, наїзд фіксується. Проміжний стан нестійкий, виявляється. Тільки гарячі і холодні стану стабільні. Таким чином, ми говоримо, що Земля бістабільна в цій спрощеній моделі.
Бачите, наскільки велике γ має бути для цього бістабільність, щоб загнутися? Воно точно правильне, коли γ=0,05, оскільки тоді ми отримуємо такий графік:
Коли інсоляція становить менш ніж приблизно 385 Вт/м2 є тільки в холодному стані. Коли вона потрапляє 385 Вт/м2, як показано зеленою лінією, раптом є два можливих температур: холодний один і набагато більш гарячого. Коли інсоляція вище, як показано чорною лінією, є три можливих температур: холодний один, і один нестабільний проміжний, і гарячу. І коли інсоляція отримує вище 465 Вт/м2, як показано червоною лінією, є тільки гарячий стан!
Чому проміжний стан нестійкий, коли воно існує? Чому дві інших рівноваги стабільні? Щоб відповісти на ці питання, ми повинні були б повернутися і вивчити вихідне рівняння:
і подивитися, що відбувається, коли ми натискаємо T трохи в сторону від одного зі своїх рівноважних значень. Це дійсно цікаво, але ми не будемо робити це сьогодні. Замість цього давайте зробити деякі висновки з того, що ми тільки що бачили. Є принаймні три моралі: математична моральна, наука модель клімату, і програмне забезпечення морально.
Математично ця модель ілюструє теорію катастроф. Як ми повільно повернути вгору у, ми отримуємо різні криві, що показують, як температура є функцією інсоляції … поки раптом крива не є графіком функції більше: вона стає нескінченно крутий в одній точці! Після цього ми отримуємо бістабільність:
γ=0,00
γ=0,01
γ=0,02
γ=0,03
γ=0,04
γ=0,05
Це називається параболічною катастрофою, і ви можете візуалізувати ці криві, як скибочки поверхні в 3D, яка виглядає приблизно як зображення:
звідси:
- Вольфрам Метуолд, Параболічна катастрофа. (Містить пакет математичний.)
Параболічна катастрофа “структурно стабільна“, тобто малі обурення не змінюють своє якісне поведінку. Ця концепція точної в теорії катастроф. Це корисна концепція, оскільки вона акцентує нашу увагу на стійких особливостях моделей: особливості, які не йдуть геть, якщо модель трохи не так, як це завжди буває.
Наскільки наука про клімат йде, одна мораль така, що він платить, щоб провести деякий час, переконавшись, що ми розуміємо, прості моделі, перш ніж ми заглибимося в більш складні. Прямо зараз ми дивимося на дуже простій моделі, але ми вже бачимо деякі цікаві явища. Вид моделі ми розглядаємо тепер називається модель Будико-продавців. Вони були вивчені з кінця 1960-х років:
- М. І. Будико, Про походження льодовикової епохи (російською мовою), Метеор. Гідрол. 2 (1968), 3-8.
- М. І. Будико, Вплив варіацій сонячної радіації на клімат Землі, Скажи нам 21 (1969), 611-619.
- Уільям Д. Селлерз, Глобальна кліматична модель, заснована на енергетичному балансі системи Земля-атмосфера, Журнал. Прикл. метеор. 8 (1969), 392-400.
- Карл Кроуфорд і Ерланд Келлен, Зауваження про умови існування більш ніж одного стаціонарного рішення в моделях типу Будико-Селлерза, Ж. Атмос. Наук. 35 (1978), 1123-1125.
- Джеральд Р. Норт, Девід Поллард і Брюс Вілічкі, Варіаційне формулювання кліматичних моделей Будико-Селлерза, Ж. Атмос. Наук. 36 (1979), 255-259.
Він також платить, щоб порівняти наші моделі до реальності! Наприклад, графіки, які ми бачили показують деякі дивно гарячі і холодні температури для Землі. Це трохи нервує. Давайте досліджуємо. Припустимо, ми встановлюємо γ=0 на повзунку. Тоді коальбедо Землі стає незалежною від температури: це 0,525, на півдорозі між її зледенілих і вільних від льоду цінностей. Потім, коли інсоляція бере своє фактичне значення 342,5 Вт на квадратний метр, модель говорить температура Землі дуже холодно: близько -20 °C!
Чи означає це, що модель в корені помилковим? Можливо ні! Зрештою, вона заснована на дуже світлою Землі. Припустимо, що ми використовуємо фактичну альбедо Землі. Звичайно, це важко визначити, набагато менше визначити. Але давайте подивимося на деякий середнє значення альбедо Землі: мовляв, це близько 0,3. Це дає coalbedo c = 0,7. Якщо ми підключи, що в нашій формулі:
отримаємо 11 °C. Це не дуже далеко від фактичної середньої температури Землі, а саме близько 15 ° С Так зимно температура -20 ° C, здається, виходить від Землі, що дуже багато світліше, ніж у нас.
Наша модель включає в себе парниковий ефект, так як коефіціенти А і В були визначені супутникових вимірювань, скільки радіації насправді йде в атмосферу Землі і вистрілює в космос.Як подальшої перевірки нашої моделі, ми можемо поглянути на ще більш простої моделі балансу нульмерние енергії: повністю чорний Землю, без парникового ефекту. Інший учасник проекту Азимут написав про це:
- Тім ван Бік, Розміщення Землі в коробці, Азимут, 19 червня 2011 р.
- Тім ван Бік, Квант тепла, Азимут, 2 липня 2011 р.
Як він пояснює, ця модель дає Землю температуру 6 ° С. Він також показує, що в цій моделі, зниження альбедо до реалістичного значенням 0,3 знижує температуру допрохолоди -18° C. Для того, щоб отримати від того, до чого-то на кшталт нашої Землі, ми повинні прийняти парниковий ефект до уваги.
Такого роду махінацій навколо свого роду річ, яку ми повинні зробити, щоб вивчити недоліки і переваги моделі клімату. Звичайно, будь-яка реалістична модель клімату є набагато складнішим, ніж маленької іграшкою ми дивилися на, так що ‘дріб’язковий навколо’ також повинні бути більш складними. За допомогою більш складної моделі, ми також можемо бути більш вимогливими. Наприклад, коли я сказав, 11° С “не дуже далеко від фактичної середньої температури Землі, а саме близько 15° С”, я був дуже пересичений про те, що насправді велика розбіжність. Я тільки прийняв цю позицію, тому що розрахунки ми робимо зараз є попередніми.
Нарешті, ось що Аллан повинен сказати про програмне забезпечення, ви тільки що бачили, і деякі з вправнішим програмного забезпечення ви побачите в найближчі тижні:
Ваш початковий питання в форумі Азимут був “що це найпростіший спосіб писати прості програми такого роду, які можуть бути доступні і управляється неосвічених людей в Інтернеті?” А “проста програма” для моделі клімату ви запропонували необхідні два елементи: засіб для вирішення ЗДР (звичайного диференціального рівняння), яка описує модель, і засіб для взаємодії з і візуалізувати результати для (ясно) “неосвічених людей в Інтернеті”.
Деякі хороші пропозиції були зроблені членами форуму:
- використовувати повноцінний пакет чисельних обчислень, таких як шавлія або Matlab, які приходять до зубів навантажена із ЗДР решателей і інтерактивної побудови графіків;
- використовувати повнофункціональний мова програмування, як Java, який має бібліотеки, доступні для вирішення ЗДР і побудови графіків, і які можуть бути упаковані як аплет для Інтернету;
- робити всі обчислення і візуалізацію себе у Javascript.
У той час як перші дві пропозиції були вище для обчислення оді рішення, я знав, що з гіркого досвіду (як розробник програмного забезпечення), що по-справжньому неосвічені люди були нам сміливі учасники форуму, які беруть участь в цій новій онлайн підприємства: ніхто з нас не були експертами в цій інтерактивній/онлайн математики річ, і програмування нового програмного забезпечення майже завжди складніше, ніж ви очікуєте, що це буде.
Тоді насправді випускати нове програмне забезпечення ще складніше! Особливо перед аудиторією, як великий, як ваших читачів. Для того, щоб придумати інтерактивне рішення, яке буде працювати на багатьох різних комп’ютерах/браузерах, сама буденна і пішохід навіювання “робити все це в себе Javascript і мати їх запустити його в браузері” був також, швидше за все, буде успішним.
Проблема з Javascript, що не так багато людей використовують його для чисельних обчислень, і я впав на наші шанси на успіх, поки Стаффан не вказав на відмінне програмне забезпечення JSXGraph. JSXGraph має багато прикладів, доступних, щоб отримати і працює, має вирішувач ЗДР, і після копіювання/вставки або два і деякі настройки з мого боку ми були все готово.
Істинне виправдання для олл-Javascript, хоча в тому, що ви були згодом змогли зробити деякі копіювати/вставити свій власний безпосередньо в TWF без будь-яких серверів, які потребують налаштоване і т.д., або навіть будь-якої допомоги від мене! Графіки повинні бути доступні для перегляду ваших читачів до тих пір, поки браузери підтримують JavaScript (ознака хорошої версії програмного забезпечення є те, що ви не повинні думати про це потім).
Є деякі поліпшення я хотів би зробити, щоб, як ми обробляємо майбутні проекти, які ми обговорювали в форумі. Перш за все, за допомогою Javascript, щоб зробити всі наші чисельна робота не буде залучати кращі і яскраві уми з форуму (або в іншому місці), щоб допомогти з подальшими моделями. Моя особиста надія полягає в тому, що ми дозволяємо всім числова роботи належить зробити на будь-якій мові, люди відчувають себе продуктивним с, і що ми придумали пляма спосіб для вас, щоб вбудувати і взаємодіяти тільки з даними з цих моделей в ваших веб-сторінок. Глін Еджі і Джим Статтард, здається, є певний імпульс у цьому напрямку.
Або, можливо, створення і редагування інтерактивних математику онлайн в кінцевому підсумку стати настільки ж легко, як вікі-сторінок сьогодні, я знаю, що сказав Стаффан розробники Sage шукали, щоб зробити їх онлайн навчальні посібники більш інтерактивним. Крім того, яскраві люди за новою мовою Джулії обговорюють способи запуску (і, імовірно, взаємодіють з) Julia в хмарі. Так що, можливо, ми повинні просто затягли наші ноги над цим проектом протягом декількох років для всього цього цікавого матеріалу, щоб допомогти нам! (І давайте чекати Singularity поки ми на нього.)
Ні, давайте не будемо! Я сподіваюся, що ви програмісти там можуть допомогти нам знайти гарні рішення проблем, з якими стикаються Аллана. І я сподіваюся, що деякі з вас насправді приєднатися до команди.
До речі, Алан має дещо багатшу версію тієї ж моделі Будико-Селлерза тут.
Обговорення цього питання на “знахідки цього тижня” див. на Азимуті. І якщо ви хочете взяти участь в створенні інтернет-кліматичних моделей, зв’яжіться зі мною і/або приєднайтеся до форуму Азимут.
Таким чином, даний тепловий режим і оледенений Землі доводять, характеризується високою нестабільністю. Порівняно невеликі зміни радіації тільки на 1,0-1,5% – достатні для розвитку крижаного покриву на суші і в океані, який досягає помірних широт. – М. І. Будико
© 2012 John Baez
[email protected]