Pokroky matematiky, fyziky a astronomie
Marian Karlický Mohou na Slunci vzniknout supererupce? Pokroky matematiky, fyziky a astronomie, Vol. 59 (2014), No. 3, 187--193
Persistent URL: http://dml.cz/dmlcz/144024
Terms of use: © Jednota českých matematiků a fyziků, 2014 Institute of Mathematics of the Academy of Sciences of the Czech Republic provides access to digitized documents strictly for personal use. Each copy of any part of this document must contain these Terms of use. This paper has been digitized, optimized for electronic delivery and stamped with digital signature within the project DML-CZ: The Czech Digital Mathematics Library http://project.dml.cz
Mohou na Slunci vzniknout supererupce? Marian Karlick´y, Ondˇrejov ´ 1. Uvod Sluneˇcn´ı erupce je explozivn´ı proces ve sluneˇcn´ı atmosf´eˇre, kter´ y vznik´a rychlou pˇremˇenou magnetick´e energie v okol´ı sluneˇcn´ıch skvrn na ohˇrev a pohyby plazmatu, urychlov´an´ı ˇc´astic a z´aˇren´ı v ˇsirok´em oboru spektra od radiov´ ych vln pˇres optick´e, rentgenov´e aˇz po gama z´aˇren´ı [9], [3], [4]. Typick´e mnoˇzstv´ı energie uvolnˇen´e v tomto procesu, trvaj´ıc´ım nˇekolik minut aˇz nˇekolik hodin, je 1022 –1025 J [7].
Obr. 1. Sch´ema modelu mohutn´e sluneˇcn´ı erupce podle P. Gallaghera (http://solarmuri. ssl.berkeley.edu/∼hhudson/cartoons). Obr´ azek zn´ azorˇ nuje vyvrˇzen´ı filamentu, tj. magnetick´e smyˇcky, kterou teˇce elektrick´ y proud aˇz 1012 A. Pod n´ı se vytv´ aˇr´ı proudov´ a vrstva, kde doch´ az´ı k propojov´ an´ı magnetick´ ych siloˇc´ ar (rekonexi) se souˇcasn´ ym uvolˇ nov´ an´ım magnetick´e energie, tj. k prim´ arn´ımu erupˇcn´ımu procesu. Ve spodn´ı ˇc´ asti erupce se vytv´ aˇr´ı ark´ ada erupˇcn´ıch smyˇcek, jejichˇz zakotven´ı jsou pozorovateln´ a jako erupˇcn´ı vl´ akna. U nejmohutnˇejˇs´ıch erupc´ı jsou tato vl´ akna pozorovateln´ a i ve viditeln´em oboru spektra, viz tak´e b´ıl´ a vl´ akna na obr´ azku 5. Detailnˇejˇs´ı popis sluneˇcn´ı erupce lze naj´ıt v prac´ıch [9], [3], [4].
ˇ ´ , DrSc., Astronomick´ Doc. RNDr. Marian Karlicky yu ´stav, Akademie vˇed Cesk´ e republiky, 251 65 Ondˇrejov, e-mail:
[email protected]
Pokroky matematiky, fyziky a astronomie, roˇcn´ık 59 (2014), ˇc. 3
187
Obr. 2. Frekvence v´ yskytu sluneˇcn´ıch erupc´ı spolu s frekvenc´ı v´ yskytu stel´ arn´ıch supererupc´ı v z´ avislosti na energii erupce (pˇrevzato z pr´ ace [10]). 1 erg = 10−7 J
Sch´ema modelu sluneˇcn´ı erupce je zn´azornˇeno na obr´azku 1. K erupc´ım podobn´eho typu doch´az´ı i na hvˇezd´ach. V nˇekter´ ych pˇr´ıpadech energie, uvolnˇen´a v tˇechto stel´arn´ıch erupc´ıch, podstatnˇe pˇrevyˇsuje energii maxim´aln´ıch pozorovan´ ych sluneˇcn´ıch erupc´ı; dosahuje hodnot aˇz 1026 –1031 J [8]. K nejmohutnˇejˇs´ım erupc´ım doch´az´ı zejm´ena u mlad´ ych a dvojn´ ych hvˇezd. Tyto erupce, kter´e jsou tis´ıckr´at a v´ıcekr´at mohutnˇejˇs´ı neˇz dosud pozorovan´e nejmohutnˇejˇs´ı erupce na Slunci, se naz´ yvaj´ı supererupce. Z ned´avn´e doby je zn´am pˇr´ıpad mohutn´e sluneˇcn´ı erupce (intenzity X4.6 podle GOES klasifikace), kter´a zp˚ usobila silnou geomagnetickou bouˇri (∼540 nT) a vyˇradila z provozu celou elektrickou s´ıt’ v oblasti Quebec v Kanadˇe. Energie t´eto erupce byla odhadnuta na 1025 J. Na z´akladˇe t´eto ud´alosti se lze pt´at, zdali nem˚ uˇze doj´ıt k jeˇstˇe mohutnˇejˇs´ım sluneˇcn´ım erupc´ım s jeˇstˇe vˇetˇs´ımi ˇskodami na pozemsk´em prostˇred´ı a zaˇr´ızen´ıch modern´ı civilizace. Jak ukazuje obr´ azek 2, ˇcetnost v´ yskytu sluneˇcn´ıch erupc´ı N kles´a s r˚ ustem jejich energie E podle mocninn´e z´avislosti dN/dE ∼ E −1.8 . Napˇr´ıklad nejmohutnˇejˇs´ı sluneˇcn´ı erupce s energi´ı 1025 J (tj. 1032 erg) se vyskytuj´ı pˇribliˇznˇe jednou za rok. Na z´akladˇe tohoto grafu se d´a i extrapolovat ˇcetnost v´ yskytu moˇzn´ ych sluneˇcn´ıch supererupc´ı. Nen´ı vˇsak jasn´e, zda-li je takov´a extrapolace opr´avnˇen´a. Na druh´e stranˇe, doplnit statistiku tohoto grafu o skuteˇcn´a pozorov´an´ı sluneˇcn´ıch supererupc´ı je nere´aln´e. Ve vesm´ıru je vˇsak mnoho hvˇezd podobn´eho typu jako Slunce, a proto pozorov´an´ı erupc´ı na tˇechto hvˇezd´ach je cestou k ˇreˇsen´ı tohoto probl´emu.
188
Pokroky matematiky, fyziky a astronomie, roˇcn´ık 59 (2014), ˇc. 3
2. Supererupce na hvˇ ezd´ ach Maehara a kol. [6] hledali erupce na hvˇezd´ach sluneˇcn´ıho typu G hlavn´ı posloupnosti v datech mise Kepler NASA [5], pozorovan´ ych bˇehem 120 dn˚ u v roce 2009. Pro v´ ybˇer hvˇezd uvaˇzovali efektivn´ı teplotu Teff v rozmez´ı 5100–6000 K a povrchovou gravitaci log(g) ≥ 2, kde g je v jednotk´ach m s−2 . V anal´ yze v´ıce neˇz 83 000 hvˇezd nalezli 365 supererupc´ı na 148 hvˇezd´ach sluneˇcn´ıho typu G s energi´ı vˇetˇs´ı neˇz 1026 J. Typick´a d´elka supererupc´ı byla nˇekolik hodin a jejich amplituda byla 0.1–1 % hvˇezdn´e z´aˇrivosti. Chyba energi´ı supererupc´ı byla odhadnuta na 60 %. Pˇr´ıklad pozorov´ an´ı supererupc´ı na dvou hvˇezd´ach (KIC 9459362 a KIC 6034120) je na obr´azku 3. Perioda variac´ı svˇeteln´e kˇrivky hvˇezdy KIC 9459362 byla 12.5 dne a hvˇezdy KIC 6034120 5.7 dne. Obecnˇe tyto variace mohou b´ yt zp˚ usobeny rotac´ı hvˇezdy spolu s hvˇezdn´ ymi skvrnami, pohyby hvˇezd po eliptick´ ych drah´ach ve dvojn´em syst´emu, zakr´ yv´an´ım jedn´e hvˇezdy druhou a nebo pulzac´ı hvˇezdy. Vˇsechny tyto moˇznosti musely b´ yt analyzov´any. Bylo vˇsak zjiˇstˇeno, ˇze variace svˇeteln´ ych kˇrivek hvˇezd KIC 9459362 a KIC 6034120 jsou zp˚ usobeny rotac´ı hvˇezdy se skvrnami. Objev´ı-li se totiˇz na disku hvˇezdy tmav´a a chladn´a skvrna, z´ aˇrivost hvˇezdy klesne. Tato z´aˇrivost opˇet vzroste, kdyˇz skvrna pˇri rotaci hvˇezdy zapadne za okraj hvˇezdy. Perioda t´eto variace z´aˇrivosti proto odpov´ıd´a rotaˇcn´ı periodˇe hvˇezdy. Hloubka modulace t´eto variace je pak d´ana velikosti skvrny, coˇz umoˇzn ˇuje odhadnout jej´ı velikost. Pokud dojde na hvˇezdˇe k erupci, vyz´aˇren´a energie t´eto erupce je pak superponov´ana na z´akladn´ı rotaˇcn´ı variaci svˇeteln´e kˇrivky tak, jak se ukazuje na obr´azku 3. Z intenzity z´aˇren´ı supererupce a vzd´alenosti a parametr˚ u hvˇezdy se pak vypoˇc´ıt´a celkov´a energie supererupce. Variace svˇeteln´ ych kˇrivek na obr´azku 3 ukazuj´ı, ˇze na hvˇezd´ach KIC 9459362 a KIC 6034120 doˇslo nejenom k supererupc´ım, ale ˇze na tˇechto hvˇezd´ach existuj´ı velik´e skvrny. V analogii s erupˇcn´ımi procesy na Slunci lze supererupce na tˇechto hvˇezd´ach vysvˇetlit uvolnˇen´ım magnetick´e energie, naakumulovan´e v bl´ızkosti hvˇezdn´ ych skvrn. Data vˇsech zjiˇstˇen´ ych supererupc´ı byla zpracov´ana a v´ ysledky jsou shrnuty v obr´azc´ıch 2 a 4. Obr´azek 2 srovn´av´a ˇcetnost v´ yskytu hvˇezdn´ ych supererupc´ı (Maehara et al. [6]) a erupc´ı na Slunci. Histogram v´ yskytu supererupc´ı v podstatˇe odpov´ıd´a extrapolovan´e mocninn´e z´avislosti. Ukazuje se, ˇze analyzovan´e supererupce se vyskytuj´ı s ˇcetnost´ı jednou za 800 aˇz 5000 let. Obr´azek 4 pak ukazuje, ˇze energie supererupc´ı roste s rostouc´ı velikost´ı hvˇezdn´ ych skvrn. Zd´a se proto, ˇze nutnou podm´ınkou supererupc´ı jsou neobvykle velk´e skvrny, kter´e zauj´ımaj´ı plochu aˇz 10−2 –10−1 celkov´eho povrchu hvˇezdy. Proto, pokud by mˇelo doj´ıt k supererupci na Slunci, podm´ınkou jsou ohromn´e sluneˇcn´ı skvrny, viz obr´azek 5. Za norm´aln´ıch okolnost´ı je i u tˇech zat´ım nejmohutnˇejˇs´ıch sluneˇcn´ıch erupc´ı plocha skvrn podstatnˇe menˇs´ı, viz obr´azek 4. D´ale Maehara a kol. [6] zjistili, ˇze poˇcet supererupc´ı je vˇetˇs´ı u hvˇezd, kter´e maj´ı kratˇs´ı rotaˇcn´ı periodu, tj. u hvˇezd podstatnˇe mladˇs´ıch neˇz je Slunce. To znamen´a, ˇze v dalˇs´ım studiu t´eto problematiky bude nutn´e hledat supererupce na hvˇezd´ach typu G s pokud moˇzno stejnou periodou rotace jako Slunce (stˇredn´ı hvˇezdn´a perioda rotace Slunce je 25.38 dn´ı). Vˇsechna tato zjiˇstˇen´ı vedou k cel´e ˇradˇe teoretick´ ych ot´azek, zdali procesy na Slunci mohou skuteˇcnˇe v´est ke vzniku sluneˇcn´ı supererupce. Napˇr´ıklad ot´azka: Je dynamo proces, kter´ y na Slunci generuje magnetick´e pole ve sluneˇcn´ıch skvrn´ach dostateˇcn´ y Pokroky matematiky, fyziky a astronomie, roˇcn´ık 59 (2014), ˇc. 3
189
Obr. 3. Typick´e svˇeteln´e kˇrivky hvˇezd se supererupcemi. (a) Svˇeteln´ a kˇrivka se supererupc´ı na hvˇezdˇe G-typu KIC 9459362. BJD – barycentrick´e juli´ ansk´e datum. (b) Zvˇetˇsen´ a svˇeteln´ a kˇrivka supererupce, d znamen´ a dny od zaˇc´ atku supererupce. Odhadnut´ a energie supererupce je 5.63 · 1027 J. (c) Stejn´e jako (a), ale pro supererupci na hvˇezdˇe G-typu KIC 6034120. (d) Stejn´e jako (b). Odhadnut´ a energie t´eto supererupce je 3.03 · 1028 J (pˇrevzato z pr´ ace [6]).
190
Pokroky matematiky, fyziky a astronomie, roˇcn´ık 59 (2014), ˇc. 3
Obr. 4. Z´ avislost energie erupce na ploˇse skrvn: pln´e ˇctvereˇcky v obr´ azku nahoˇre jsou pro supererupce na hvˇezd´ ach sluneˇcn´ıho typu a pln´ a koleˇcka dole pro sluneˇcn´ı erupce. Obr´ azek je doplnˇen rentgenov´ ym (X) tokem pro pˇr´ıpad sluneˇcn´ıch erupc´ı a vypoˇcten´ ymi hodnotami magnetick´eho pole a magnetick´eho toku ve skvrn´ ach na hvˇezd´ ach (pˇrevzato z pr´ ace [10]). 1 maxwell (Mx) = 10−8 Wb
pro generaci vysok´ ych magnetick´ ych pol´ı a tok˚ u, nutn´ ych pro supererupci? Nebo: Jak dlouho by takov´ y proces trval? Je pro supererupce nutn´a pˇr´ıtomnost hork´e planety velikosti Jupitera, jak je nˇekdy diskutov´ano? Na tyto ot´azky se pokusili odpovˇedˇet ve sv´e pr´aci Shibata a kol. [10]. Autoˇri vyˇsli z Faradayovy indukˇcn´ı rovnice ∂B = ∇ × (v × B) + η∇2 B, ∂t
(1)
kde B je vektor magnetick´e indukce, ∇ · B = 0, v je rychlost proudˇen´ı sluneˇcn´ıho plazmatu a η je magnetick´a difuzivita. Pˇri zanedb´an´ı efektu magnetick´e difuzivity uk´azali, ˇze potˇrebn´ y magnetick´ y tok pro supererupci (1016 Wb) lze vygenerovat asi Pokroky matematiky, fyziky a astronomie, roˇcn´ık 59 (2014), ˇc. 3
191
Obr. 5. Pˇredstava moˇzn´e supererupce na Slunci (b´ıl´ a vl´ akna), doprov´ azen´ a skupinou neobvykle velk´ ych sluneˇcn´ıch skvrn, pozorovan´ a ve viditeln´em oboru spektra
za 40 let, coˇz je doba delˇs´ı neˇz je sluneˇcn´ı cyklus (11 let), ale podstatnˇe kratˇs´ı neˇz odhadovan´ y interval mezi supererupcemi. Z˚ ust´av´a ovˇsem ot´azka, nakolik m˚ uˇze ovlivnit tento proces zesilov´ an´ı magnetick´eho toku magnetick´a difuzivita. Autoˇri se zab´ yvali tak´e probl´emem role hork´ ych planet velikosti Jupitera, ale jejich pˇr´ıtomnost v bl´ızkosti studovan´ ych hvˇezd se neprok´azala. 3. Z´ avˇ er Na z´akladˇe vˇsech tˇechto v´ ysledk˚ u nelze vylouˇcit vznik supererupc´ı s energiemi 1027 28 aˇz 10 J a jejich ˇcetnosti jednou za 800–5000 let i na naˇsem Slunci. Pokud by k takov´e supererupci na Slunci doˇslo, lze zat´ım jen spekulovat, jak´e katastrof´aln´ı n´asledky by takov´a supererupce mohla m´ıt na podm´ınky ˇzivota na Zemi. V prvn´ı ˇradˇe by nejv´ıce byly ohroˇzeny vˇsechna elektrick´a a poˇc´ıtaˇcov´a zaˇr´ızen´ı, elektr´arny, elektrick´e a komunikaˇcn´ı s´ıtˇe, astronauti na obˇeˇzn´e dr´aze Zemˇe atd. Pokud by nav´ıc doˇslo i k devastaci ozonov´e vrsty Zemˇe, pak by vˇse bylo ohroˇzeno nebezpeˇcn´ ym z´aˇren´ım. V tˇechto souvislostech vznikaj´ı i dalˇs´ı ot´azky, napˇr´ıklad, zdali nˇekter´e jeˇstˇe vˇetˇs´ı supererupce, ke kter´ ym mohlo doj´ıt v minulosti na Slunci, nezp˚ usobily nˇekter´e diskutovan´e katastrofy v historii naˇs´ı Zemˇe. 192
Pokroky matematiky, fyziky a astronomie, roˇcn´ık 59 (2014), ˇc. 3
Tento ˇcl´anek je shrnut´ım v´ ysledk˚ u prac´ı prof. Kazunari Shibaty z Kyoto University, Japonsko, a jeho spolupracovn´ık˚ u. Tyto v´ ysledky prof. Kazunari Shibata prezentoval na ned´avn´e mezin´arodn´ı konferenci Solar and stellar flares: Observations, simulations and synergies konan´e v Praze 23–27. ˇcervna 2014 i ve sv´e pˇredn´aˇsce pro veˇrejnost dne ˇ v Praze. 26. ˇcervna 2014 v budovˇe Akademie vˇed CR Podˇ ekov´ an´ı. Autor dˇekuje prof. Kazunari Shibatovi za poskytnut´ı jeho ˇcl´ank˚ u ˇ a obr´azk˚ u. Tato pr´ace byla podpoˇrena grantem P209/12/0103 GA CR. Literatura [1] Aschwanden, M. J., Tarbell, T. D., Nightingale, R. W., Schrijver, C. J., Title, A., Kankelborg, C. C., Martens, P., Warren, H. P.: Time variability of the “quiet” Sun observed with TRACE, II. Physical parameters, temperature evolution, and energetics of extreme-ultraviolet nanoflares. Astrophys. J. 535 (2000), 1047–1065. [2] Crosby, N. B., Aschwanden, M. J., Dennis, B. R.: Frequency distributions and correlations of solar X-ray flare parameters. Solar Phys. 143 (1993), 275–299. ´ , M.: Solar flares: Radio and X-ray signatures of magnetic reconnection pro[3] Karlicky cesses. Res. Astronom. Astrophys. 14 (7) (2014), 753–772. ´ , M., Ba ´ rta, M.: Energetick´e kask´ [4] Karlicky ady v rekonexi magnetick´eho pole. Vesm´ır 92 (10) (2013), 558–561. [5] Koch, D. G., et al.: Kepler mission design, realized photometric performance, and early science. Astrophys. J. 713 (2010), L79–L86. [6] Maehara, H., Shibayama, T., Notsu, S., Notsu, Y., Nagao, T., Kusaba, S., Honda, S., Nogami, D., Shibata, K.: Superflares on solar-type stars. Nature 485 (2012), 478–481. [7] Priest, E. R.: Solar flare magnetohydrodynamics. Gordon and Breach Science Publishers, New York, 1981. [8] Schaefer, B. E.: Flashes from normal stars. Astrophys. J. 337 (1989), 927–933. [9] Shibata, K., Magara, T.: Solar flares: Magnetohydrodynamic processes. Living Rev. Sol. Phys. 8 (6) (2011), 1–99. [10] Shibata, K., Isobe, H., Hillier, A., Choudhuri, A. R., Maehara, H., Ishii, T., Shibayama, T., Notsu, S., Notsu, Y., Nagao, T., Honda, S., Nogami, D.: Can superflares occur on our Sun? Publ. Astron. Soc. Japan 65 (49), (2013), 1–8. [11] Shimizu, T.: Energetics and occurrence rate of active-region transient brightenings and implications for the heating of the active-region corona. Publ. Astron. Soc. Japan 47 (1995), 251–263.
Pokroky matematiky, fyziky a astronomie, roˇcn´ık 59 (2014), ˇc. 3
193