Sunspot observation at a spatial resolution of 0.08 acrsec with GREGOR

The various mechanisms of magneto-convective energy transport determine the structure of sunspots. Using speckle-reconstructed images in the 486 nm band, we characterise the appearance of light bridges.
We find branched dark lanes extending along thin (≈1′′) light bridges in sunspots at various heliocentric angles. In thick (>2′′) light bridges the branches are disconnected from the central lane and have a Y shape with a bright grain toward the umbra. The images reveal that light bridges exist on varying intensity levels and that their small-scale features evolve on timescales of minutes. Faint-light bridges show dark lanes outlined by the surrounding bright features. Dark lanes are very common and are also found in the boundary of pores. They have a characteristic width of 0.1′′ or smaller. Intergranular dark lanes of that width are seen in active region granulation. We interpret our images in the context of magneto-convective simulations and findings: While central dark lanes in thin light bridges are elevated and associated with a density increase above upflows, the dark lane branches correspond to locations of downflows and are depressed relative to the adjacent bright plasma. Thick light bridges with central dark lanes show no projection effect. They have a flat elevated plateau that falls off steeply at the umbral boundary. There, Y-shaped filaments form as they do in the inner penumbra. This indicates the presence of inclined magnetic fields, meaning that the umbral magnetic field is wrapped around the convective light bridge.
Schlichenmaier, R., von der Lühe, O., Hoch, S., et al., 2016, Astronomy & Astrophysics, Volume 596, A7
DOI: 10.1051/0004-6361/201628561
Sunspot Decay

The stability of sunspots is one of the long-standing unsolved puzzles in the field of solar magnetism and the solar cycle. The thermal and magnetic structure of the sunspot beneath the solar surface is not accessible through observations, thus processes in these regions that contribute to the decay of sunspots can only be studied through theoretical and numerical studies.
New insights were gained taking advantage of a numerical simulation of a sunspot with a lateral extension of more than 98 000 km times 98 000 km extending almost 18 000 km beneath the solar surface. These simulations had been performed by Matthias Rempel at the High Altitude Observatory (HAO). A data set of 30 hours, showing a stable sunspot at the solar surface, was analysed with respect to sunspot decay as a PhD project at KIS, in collaboration with the HAO.
This simulation allows to study the non-local dynamics of fluting instability interacting with stabilising buoyance forces within in the magnetic flux tube that constitutes the spot. The simulation shows that the fluting instability causes a corrugation along the circumference of the flux tube (sunspot) which proceeds fastest at a depth around 7500 km beneath the solar surface. For a long time, the decay of the spot is suppressed by the stabilising effect close to the surface which holds the bundle of magnetic field lines together. Yet, ultimately, large intrusions of field-free plasma below the surface form convective bubbles which might form light-bridges in the spot at the solar surface. This can finally disrupt the spot.
The analysis of the simulation also allow to assess the results of a linear stability analysis developed in the 1970s (Meyer, Schmidt, & Weiss 1977). As predicted back then, the magnetic field configuration of a sunspot in the uppermost few 1000 km is stable, mostly due to buoyancy forces. Deeper layers in the spot are unstable to fluting (interchange) instability as it was predicted.
Strecker, H., Schmidt, W., Schlichenmaier, R., Rempel, M., Astronomy & Astrophysics, 2021, Volume 649, A123
DOI: 10.1051/0004-6361/202040199
Hanna Strecker, PhD Thesis, Universität Freiburg, https://freidok.uni-freiburg.de/data/165760
The magnetic field at the boundary between umbra and penumbra




Sunspots are the longest-known manifestation of solar activity and their magnetic nature has been known for more than a century. Despite this, the boundary between umbrae and penumbrae, the two fundamental sunspot regions, has hitherto been solely defined by an intensity threshold.
Using Hinode satellite data, we can give statistical proof of an empirical law governing the boundary between umbra and penumbra in stable sunspots: The vertical component of the magnetic field, Bver, is invariant.
Contours of intensity at 50% of the granular quiet-sun value and contours of Bver=1867 G match. This is illustrated in the Fig. 1, in which contours for intensity (white) and Bver (red) are displayed for some of the spots that where included in the sample. Deviations were studied and shown to be non-significant (Jurcak et al. 2018).
This empirical law applies for all sunspots of the sample with different sizes, morphologies, evolutionary stages, and phases of the solar cycle. Thereby we have unveiled the magnetic property that discriminates between the umbral and penumbral modes of magneto-convection. This discovery carries strong consequences for the understanding of the fundamental processes of energy transport occurring in the magnetic Sun.
The validity of this empirical law and the value of this new constant is now being examined for different wavelengths and data sources. Figure 2 shows intensity images of sunspots observed with the GRIS@GREGOR spectograph. The value of Bver=1843 Gauss is consistent with the results from Jurcak et al 2018 with Hinode data.
We also studied the temporal evolution of Bver at the umbral boundary, in this case using HMI/SDO data. During the first disc passage, NOAA AR 11591, Bver remains constant at 1693 G with a root-mean-square deviation of 15 G, whereas the magnetic field strength varies substantially (mean 2171 G, rms of 48 G) and shows a long-term variation as shown in the residual in Fig. 3 after daily variations due to the satellite orbit are corrected. Hence, during the disc passage of a stable sunspot its umbral boundary can equivalently be defined by using the continuum intensity Ic or the vertical magnetic field component Bver. Contours of fixed magnetic field strength fail to outline the umbral boundary (Schmassmann et al. 2018).
Now, since observations have shown that in stable sunspots the umbral boundary is outlined by a critical value of the vertical magnetic field component, the key question is: What is the nature of the distinct magnetoconvection regimes in the umbra and penumbra? This is still unclear. To follow up on this, we analysed a sunspot simulation in an effort to understand the origin of the convective instabilities giving rise to the penumbral and umbral distinct regimes. We applied the criterion from Gough & Tayler (1966, MNRAS, 133, 85), accounting for the stabilising effect of the vertical magnetic field to investigate the convective instabilities in a MURaM sunspot simulation. The result is shown in Fig. 4: The deep subphotospheric sunspot trunk is magneto-convectivly stable. Upwards from some -5 Mm, columns of instabilty characterise the umbra of the sunspot. The subsphotospheric penumbra exhibits filamentary instabilities that are almost parallel to the solar surface. Hence, the Gough and Tayler criterion as a diagnostic tool reveals the tripartite nature of sunspot structure with distinct regimes of magneto-convection in the umbra, penumbra, and granulation operating in realistic MHD simulations. (Schmassmann et al. 2021)
Jurčák, J., Rezaei, R., Bello González, N. Schlichenmaier, R., Vomlel, J., 2018, Astronomy & Astrophysics, Volume 611, L4
DOI: 10.1051/0004-6361/201732528
Schmassmann, M., Schlichenmaier, R. 2018, Astronomy & Astrophysics, Volume 620, A104
DOI: 10.1051/0004-6361/201833441
Lindner, P., Schlichenmaier, R., Bello González N., 2020, Astronomy & Astrophysics, Volume 638, A25
DOI: 10.1051/0004-6361/202037716
Schmassmann, M., Rempel, M., Bello González, N., Schlichenmaier, R., Jurčák, J., 2021, Astronomy & Astrophysics, Volume 656, A92
DOI: 10.1051/0004-6361/202141607
Lösung der Strahlungstransportgleichung im geometrischen Maßstab

Das Vorhandensein von Magnetfeldern in der Sonne führt zu einer Fülle von magnetodynamischen Phänomenen, die das Interesse vieler Sonnenphysiker*innen geweckt haben (und noch wecken). Die Untersuchung dieser physikalischen Phänomene hängt von der genauen Ableitung der physikalischen Eigenschaften aus den aufgezeichneten Spektren ab.
In diesem Projekt wollen wir einen numerischen Code entwickeln, der die Strahlungstransportgleichung zusammen mit der magneto-hydrostatischen Gleichung für einen Bereich der Sonnenoberfläche selbstkonsistent löst. Zu diesem Zweck ist es zweckmäßig, den Löser der Strahlungstransportgleichung in einem geometrischen Maßstab zu implementieren, da dieser Punkt die Berechnung der räumlichen Ableitungen vereinfacht, die zur Lösung der oben genannten Gleichung erforderlich sind. Ein erster Schritt besteht darin, die Strahlungstransportgleichung zu lösen, sofern eine Höhenschichtung eines Satzes von physikalischen Parametern vorliegt. In der Abbildung demonstrieren wir die (parallelen) Möglichkeiten dazu, indem wir die synthetischen Polarisationskarten für die Spektrallinie Fe I 6173 Å für eine Momentaufnahme einer Sonnenfleckensimulation zeigen (Rempel, M. 2012, ApJ, 750, 62). Die Synthese benötigte 125,13 CPU-Stunden, die, aufgeteilt auf 40 Prozessoren, in 3,12 Stunden erledigt waren.
Der numerische Code FIRTEZ-DZ, der die Vorwärts- und Rückwärts-Strahlungstransportgleichung löst, wurde von Pastor Yabar et al. (2019), A&A, 629, 24. Der Code ist hier öffentlich zugänglich.
Thermische Struktur der Penumbra

Die thermische Struktur der Penumbra unterhalb ihrer sichtbaren Oberfläche hat wichtige Auswirkungen auf unser heutiges Verständnis von Sonnenflecken und ihrer Penumbrae: welche magneto-konvektive Mode transportiert Energie, und wie werden magneto-akustische Wellenmoden in der Sonnenfleckenseismologie umgewandelt.
Wir analysierten spektropolarimetrische Daten (d. h. Stokes-Profile) in drei Eisenlinien bei 1565 nm, die mit dem GRIS-Instrument am 1,5-Meter-Sonnenteleskop GREGOR beobachtet wurden. Die Daten werden um die Verschmierungseffekte des weitwinkligen Streulichts korrigiert und dann einem Inversionscode für die Strahlungstransportgleichung unterzogen, um unter anderem die Temperatur als Funktion der optischen Tiefe zu ermitteln. Wir stellen fest, dass der mittlere Temperaturunterschied zwischen den beiden durch etwa 75 km getrennten Schichten in den hellen Penumbra-Filamenten etwa 1057 K beträgt, während er in der ruhigen Sonne 1390 K beträgt. Dies deutet darauf hin, dass der Temperaturgradient unterhalb der sichtbaren Oberfläche der hellen Penumbra-Filamente kleiner ist als in der ruhigen Sonne, und es impliziert, dass die Temperatur in den subphotosphärischen Schichten von der der ruhigen Sonne abweicht. Wir interpretieren diese Ergebnisse als Beweis für eine dicke Penumbra, wobei die Magnetopause nicht in der Nähe ihrer sichtbaren Oberfläche liegt. Diese Ergebnisse wurden veröffentlicht von Borrero et al. (2017), A&A, 601, L8
Magnetischer Rückfluss in der Penumbra eines Sonnenflecks

Ein Höhepunkt der ersten wissenschaftlichen Ergebnisse von GREGOR ist die Erfassung der kleinräumigen Geometrie des Magnetfelds in einer Sonnenflecken-Penumbra. Spektralpolarimetrische Daten mit dem GRegor Infrared Spectropolarimeter (GRIS) zeigen, dass 35 % der Penumbra-Fläche mit Rückfluss bedeckt sind, d. h. mit Orten, an denen das Magnetfeld eine entgegengesetzte Polarität hat. Da die Feldstärke dieser nach unten gerichteten Komponente gering ist, macht dieser Rückfluss nur 10 % des gesamten magnetischen Flusses in der Penumbra aus. Eine sorgfältige Analyse der Tiefenabhängigkeit des Magnetfelds zeigt, dass der Rückfluss nur in den tiefsten Schichten der Photosphäre vorhanden ist, die mit den vier Eisenlinien in der spektralen Nähe von 1564,8 nm untersucht werden. Dies is der Grund dafür, dass das Spektropolarimeter an Bord von HINODE, das Eisenlinien höherer Ordnung um 630,2 nm verwendet, den Rückfluss nicht nachweisen kann, während frühere Beobachtungen bei 1564,8 nm mit dem VTT nicht die erforderliche räumliche Auflösung hatten.
Diese Erkenntnis hat wichtige Auswirkungen auf die Geschwindigkeits- und Magnetfeldstruktur der penumbralen Magnetokonvektion: Diese Flecken neigen dazu, mit den Orten zusammenzufallen, an denen der Evershed-Ausfluss am größten ist. Da die Evershed-Strömung magnetisiert ist, kann man davon ausgehen, dass sich die Strömung am Magnetfeld ausrichtet. Daher sollten diese Rückflussflecken mit Abflüssen verbunden sein. Diese Ergebnisse wurden von Borrero et al. (2016), A&A, 596, 2 veröffentlicht.
Ableitung des elektrischen Stromvektors in der Sonnenphotosphäre

Elektrische Ströme spielen eine sehr wichtige Rolle für die Energiebilanz und den Transport im Plasma der Sonnenatmosphäre und insbesondere in der Sonnenkorona. Sie sind auch ein Indikator für die magnetische Rekonnexion, die als Hauptursache für explosive und transiente Ereignisse in der Chromosphäre und Korona gilt. Folglich ist die Messung dieser elektrischen Ströme eines der Hauptziele der Sonnenphysikgemeinschaft.
Trotz ihrer Bedeutung waren routinemäßige und direkte Messungen der elektrischen Ströme j in der Sonnenatmosphäre bisher nicht allgemein möglich. Mit unserem neuen Inversionscode FIRTEZ-DZ, der an unserem Institut entwickelt wurde, sind wir in der Lage, genaue Rückschlüsse auf den vollständigen elektrischen Stromvektor j zu ziehen. Dies ist möglich, weil FIRTEZ das magneto-hydrostatische Gleichgewicht (MHS) anstelle des weit verbreiteten hydrostatischen Gleichgewichts (HE)
berücksichtigt, das von anderen Inversionscodes verwendet wird. Diese Studie ist zur Veröffentlichung in Astronomy & Astrophysics Letters angenommen worden und wird demnächst veröffentlicht werden. Einen Vorabdruck finden Sie hier: arxiv.org/abs/2112.04356
Stokes-Profil-Inversionen im magnetisch-hydrostatischen Gleichgewicht

Inversionscodes für die Strahlungstransportgleichung sind Werkzeuge, mit denen wir die Temperatur, das Magnetfeld und die Plasmageschwindigkeit in der Sonnenatmosphäre aus spektralpolarimetrischen Beobachtungen (d. h. dem Stokes-Vektor) in Spektrallinien als Funktion der optischen Tiefe bestimmen können.
Zur Umrechnung von der optischen Tiefe in die geometrische Höhe werden zusätzliche Näherungen verwendet, typischerweise die Annahme des hydrostatischen Gleichgewichts (HE). Das hydrostatische Gleichgewicht ist jedoch keine gute Näherung, wenn starke und inhomogene Magnetfelder vorhanden sind. Dies ist zum Beispiel bei Sonnenflecken der Fall. Wir haben eine neue Methode entwickelt, die stattdessen die Annäherung des magneto-hydrostatischen Gleichgewichts (MHS) verwendet, das die Auswirkungen der magnetischen Druck- und Spannungskräfte berücksichtigt.
Um unsere neue Methode zu vergleichen, haben wir die so genannte Wilson-Depression (geometrische Höhe des Ortes für die Einheit der optischen Tiefe durch einen Sonnenfleck) mit dem alten hydrostatischen Ansatz (untere Tafel) und unserer neuen magneto-hydrostatischen Methode (mittlere Tafel) bestimmt und mit den Ergebnissen dreidimensionaler numerischer magnetohydrodynamischer Simulationen (obere Tafel) verglichen. Wie man sieht, liefert der magneto-hydrostatische Ansatz Werte für die Wilson-Depression, die in viel besserer Übereinstimmung mit den numerischen Simulationen sind als der hydrostatische Ansatz.
Diese Arbeit wurde veröffentlicht von Borrero et al. (2019), A&A, 632, 111.
Seismische Untersuchung der Sonnenatmosphäre

Da die Sonnenatmosphäre geschichtet ist, d. h. aufgrund ihrer unterschiedlichen Dichte und Temperatur eine geschichtete Struktur aufweist, werden Wellen, die in die Atmosphäre wandern, gebrochen. Abhängig von der Frequenz können Wellen unterschiedliche Höhen erreichen, bevor ihr Laufweg "abgeschnitten" wird.
Mit Daten des Interferometers HELLRIDE (HELioseismological Large Regions Interferometric DEvice), das am Vakuumturm-Teleskop auf Teneriffa betrieben wird, messen wir diese Grenzfrequenz in Abhängigkeit von der Höhe.
Diese Messungen können mit theoretischen Modellen verglichen werden und geben Einblicke in die Physik dieses inhomogenen Mediums.
Bezug:
Wisniewska, A., Musielak, Z. E., Staiger, J., Roth, M., Astrophys. J. Lett. 819, L23
Die Jagd nach Strömungen im tiefen Sonneninnern

Durch die turbulenten Bewegungen des Plasmas werden nahe der Sonnenoberfläche seismische Wellen erzeugt. Dort durchqueren sie das Sonneninnere und untersuchen die Bedingungen im Inneren der Sonne.
Diese seismischen Wellen interagieren mit physikalischen Größen im Inneren der Sonne. Dies kann ausgenutzt werden, um beispielsweise die Struktur und Amplitude von Strömungen im Inneren der Sonne abzuschätzen.
KIS-Forscher*innen entwickeln mit ihren Partnern im In- und Ausland eine neue Methode zur Modellierung des Einflusses von Strömungen auf die Laufzeiten von Schallwellen im Sonneninneren. Diese Entwicklung ist notwendig, um unser Verständnis großskaliger Strömungen im tiefen Sonneninneren voranzutreiben und dadurch Erkenntnisse über den Ursprung der Sonnenaktivität zu gewinnen.
Bezug:
Böning, V.G.A., Roth, M., Zima, W., Birch, A.C., and Gizon, L. , 2016, ApJ, 824, 49
Wie Magnetfelder den Klang der Sonne verändern

Der solare 11-Jahres-Zyklus beeinflusst viele der beobachtbaren Werte der Sonne. Im Laufe des Sonnenzyklus ändern sich die Schwingungsfrequenzen der Sonne und die Amplituden dieser Schwingungen geringfügig aber messbar. Sonnenschwingungen sind stehende Schallwellen im Inneren der Sonne. Ähnlich wie bei Schallwellen hier auf der Erde ist der Druck innerhalb des Mediums der Schlüssel für die Art und Weise, wie sich diese Wellen ausbreiten. Aus diesem Grund werden diese Schwingungen auch als p-Moden (p für Druck) bezeichnet. Änderungen der p-Mode-Eigenschaften während des Sonnenzyklus sind auf die unterschiedliche Stärke des Magnetfelds im Laufe der Zeit in verschiedenen Tiefen innerhalb der Sonne zurückzuführen.
Zwischen 2009 und 2013 hat das NASA-Weltraumteleskop Kepler kontinuierlich über 150 000 Sterne beobachtet. Der Hauptzweck dieser Mission war die Entdeckung von Planeten außerhalb unseres Sonnensystems. Dieselben Daten können jedoch auch verwendet werden, um sich mit der Asteroseismologie zu beschäftigen, d. h. um stellare Eigenschaften zu erforschen, indem man ihre inhärenten Schwingungen untersucht.
Diese Daten wurden verwendet, um die Eigenschaften der Schwingungen von 24 sonnenähnlichen Sternen zu untersuchen, die von Kepler beobachtet wurden. Von den 24 untersuchten Sternen zeigten 23 zeitliche Variationen ihrer Schwingungsfrequenz. Bei 6 dieser 23 Sterne war jede Zunahme der Oszillationsfrequenz mit einer Abnahme der Amplitude verbunden, während jede Abnahme der Frequenz mit einer Zunahme der Amplitude verbunden war. Dieser Zusammenhang zwischen Frequenz und Amplitude ist von der Sonne her bekannt und gilt als Hinweis darauf, dass die beobachteten Veränderungen der Schwingungseigenschaften auf zeitliche Variationen der Magnetfelder in diesen Sternen zurückzuführen sind.
Das deutlichste Beispiel dafür zeigt das Diagramm für star KIC 8006161, auch Doris genannt. Doris ist ein sonnenähnlicher Stern. Seine Masse beträgt das 1,04-fache der Sonnenmasse, sein Radius das 0,947-fache des Sonnenradius und sein Alter 5,04 Milliarden Jahre. Während des Beobachtungszeitraums erhöhte sich die Schwingungsfrequenz dieses Sterns um fast 1µHz. Im gleichen Zeitraum nahmen die Schwingungsamplituden deutlich ab. Zum Vergleich: Während ihres 11-jährigen Zyklus der magnetischen Aktivität ändern sich die Schwingungsfrequenzen der Sonne um etwa 0,5 µHz.
Die Beobachtung der Aktivitätszyklen in einer größeren Anzahl von Sternen hilft, die Prozesse, die diese Zyklen steuern, besser zu verstehen. Selbst im Fall der Sonne kennen wir die Funktionsprinzipien des Dynamos, der den Sonnenzyklus aufrechterhält und antreibt, noch nicht vollständig. In dieser Hinsicht können asteroseismische Analysen anderer Sterne – wie die hier vorgestellte – wertvolle Zusatzinformationen liefern.
Verweise:
Stellare magnetische Aktivität und Variabilität der Schwingungsparameter: Eine Untersuchung von 24 sonnenähnlichen Sternen, beobachtet von Kepler, R. Kiefer, A. Schad, G. Davies, M. Roth, Astronomy & Astrophysics DOI: 10.1051/0004-6361/201628469
Link: NASA ADS
Helioseismische Messung der Meridionalen Strömung


In den oberflächennahem Schichten der Sonne kann man eine meridionale Strömung beobachten, die sich vom Äquator bis zu den Polen erstreckt. In unserer Gruppe haben wir eine neue helioseismische Analysemethode entwickelt und angewendet, um die Ausdehnung dieser meridionalen Strömung im tieferen Sonneninneren zu messen. Das Verfahren nutzt die Verformung tief eindringender solarer Schallwellen aufgrund der Strömung. Wir haben sechs aufeinanderfolgende Jahre von Doppler-Geschwindigkeitsdaten mit dieser Methode analysiert. Die Daten wurden zwischen 2004-2010 vom Michelson Doppler Imager (MDI) an Bord des Solar and Heliospheric Observatory (SoHO), einer Raumsonde der ESA/NASA, gewonnen. Abbildung 1 zeigt einen Querschnitt der meridionalen Strömung, gemessen unterhalb der Sonnenoberfläche zwischen 0,82-0,97 Sonnenradien in der Radius-Breiten-Ebene. Die gestrichelte Linie markiert die Sonnenoberfläche; strich-punktierte Linien markieren heliographische Breitengrade von +/-60 Grad. Positive (negative) Geschwindigkeiten entsprechen einer nordwärts (südwärts) gerichteten Strömung. Unsere Messung zeigt ein komplexes Strömungsprofil im Solarinnern mit Horizontalgeschwindigkeiten unter 50 m/s. Es besteht aus mehreren größeren und kleineren Flusszellen.
Die meridionale Strömung ist in größeren Tiefen vorhanden als in Abbildung 1 dargestellt, aber die aktuelle Analyse erlaubt noch keine Messungen der großräumigen Strömungskomponente in Tiefen über 126 Mm.
In Abbildung 2 (unten) zeigen wir die Strömungskomponente bestehend aus kleinen Strömungszellen als Funktion der Breite und Tiefe unter der Sonnenoberfläche; positive (negative) Werte entsprechen wiederum einer nordwärts (südwärts) gerichteten Strömung. Für diese Strömungskomponente erhalten wir signifikante Geschwindigkeiten bis in eine Tiefe von etwa 200 Mm. Somit liefert das Ergebnis den Beweis für eine tief eindringende meridionale Strömung, die die gesamte solare Konvektionszone durchdringt. Wir haben unsere Messung mit einer oberflächennahemn Messung der kleinräumigen Strömungskomponente (Bild 2, oben) von Komm et. al. (2004) verglichen. Diese Messung reicht von 0,5-14 Mm Tiefe. Nahe der Oberfläche stimmen beide Ergebnisse gut überein. Die bereits oberflächennah beobachteten kleinskaligen Fließzellen reichen bis tief ins Sonneninnere und ändern ihre Fließrichtung in etwa 100 Mm Tiefe.
Verweise:
Schad A., Timmer J., Roth M.: "Globaler helioseismischer Beweis für einen tief eindringenden Meridianfluss bestehend aus mehreren Flusszellen", ApJL, 778, L38-L44 (2013)
Komm, R., Howe, R., Hill, F., González Hernández, I. & Toner, C., ApJ, 631, 636 (2005)