Issue |
A&A
Volume 679, November 2023
|
|
---|---|---|
Article Number | A100 | |
Number of page(s) | 10 | |
Section | The Sun and the Heliosphere | |
DOI | https://doi.org/10.1051/0004-6361/202347180 | |
Published online | 15 November 2023 |
The short-term stability and tilting motion of a well-observed low-latitude solar coronal hole
1
Department of Physics, University of Helsinki, PO Box 64 00014 Helsinki, Finland
e-mail: stephan.heinemann@hmail.at
2
Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany
3
Predictive Science Inc., 9990 Mesa Rim Road, Suite 170, San Diego, CA 92121, USA
4
NASA/Marshall Space Flight Center, Huntsville, AL 35812, USA
5
Leibniz-Institut für Sonnenphysik (KIS), Schöneckstr. 6, 79104 Freiburg, Germany
Received:
14
June
2023
Accepted:
15
September
2023
Context. Our understanding of the solar magnetic coronal structure is tightly linked to the shape of open field regions, specifically coronal holes. A dynamically evolving coronal hole coincides with the local restructuring of open to closed magnetic field, which leads to changes in the interplanetary solar wind structure.
Aims. By investigating the dynamic evolution of a fast-tilting coronal hole, we strive to uncover clues about what processes may drive its morphological changes, which are clearly visible in extreme ultraviolet (EUV) filtergrams.
Methods. Using combined 193 Å and 195 Å EUV observations by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory and the Extreme UltraViolet Imager on board the Solar Terrestrial Relations Observatory-Ahead, in conjunction with line-of-sight magnetograms taken by the Helioseismic and Magnetic Imager, also on board the Solar Dynamics Observatory, we tracked and analyzed a coronal hole over 12 days to derive changes in morphology, area, and magnetic field. We complemented this analysis by potential field source surface modeling to compute the open field structure of the coronal hole.
Results. We find that the coronal hole exhibits an apparent tilting motion over time that cannot solely be explained by solar differential rotation. It tilts at a mean rate of ∼3.2° day−1 that accelerates up to ∼5.4° day−1. At the beginning of May the area of the coronal hole decreased by more than a factor of three over four days (from ∼13 × 109 km2 to ∼4 × 109 km2), but its open flux remained constant (∼2 × 1020 Mx). Furthermore, the observed evolution is not reproduced by modeling that assumes the coronal magnetic field to be potential.
Conclusions. In this study we present a solar coronal hole that tilts at a rate that has yet to be reported in literature. The rate exceeds the effect of the coronal hole being advected by either photospheric or coronal differential rotation. Based on the analysis we find it likely that this is due to morphological changes in the coronal hole boundary caused by ongoing interchange reconnection and the interaction with a newly emerging ephemeral region in its vicinity.
Key words: Sun: rotation / Sun: corona / Sun: magnetic fields / Sun: photosphere / methods: observational / methods: analytical
© The Authors 2023
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
The topology and dynamics of the solar atmosphere are governed by the complex interplay of open and closed magnetic fields. Open fields link the solar surface to interplanetary space, allowing outflowing plasma to be accelerated to supersonic speeds, and are the source of the heliospheric open flux (see review by Cranmer & Winebarger 2019). Open field structures lead to the formation of large-scale, low-density, low-temperature regions in the solar corona that can be observed as structures of reduced emission in the Extreme UltraViolet (EUV) and X-rays (on-disk) and also in white-light (off-limb), known as the coronal holes (see review by Cranmer 2009, and references therein).
Due to plasma outflow along open magnetic field lines, a coronal hole features an electron density that can be more than 30% lower than the surrounding quiet (closed) corona ((1 − 2.5)×108 cm−3; Warren & Hassler 1999; Hahn et al. 2011; Heinemann et al. 2021) and a temperature around 0.9 MK (e.g., Fludra et al. 1999; Saqri et al. 2020). The magnetic field structure of a coronal hole consists of a mix of open and closed magnetic fields. The open fields form magnetic funnels (Tu et al. 2005) that are rooted in photospheric footpoints. These footpoints are unipolar magnetic elements of the coronal hole’s dominant polarity. The magnetic elements contain the majority of the coronal hole’s signed flux (38 − 72%) in only a fraction of the coronal hole’s area (1 − 8%; Hofmeister et al. 2017, 2019; Heinemann et al. 2019). Most of the area of the coronal hole is covered by small closed loops whose average height is significantly lower than in quiet-Sun regions (Wiegelmann & Solanki 2004). The interplay of these open and closed fields within the coronal hole and the surrounding and global magnetic field configuration determine the properties of the resulting high-speed solar wind stream in terms of peak velocity, longitudinal and latitudinal extent, shape of the solar wind profile, and whether the stream is geoeffective (see, e.g., Garton et al. 2018; Geyer et al. 2021; Hofmeister et al. 2020, 2022; Samara et al. 2022). During solar maximum, well-defined low-latitude coronal holes, which often but not exclusively emerge from the “ashes” of decaying active regions (Karachik et al. 2010; Petrie & Haislmaier 2013), dominate the solar corona, whereas during solar minimum large-scale polar coronal holes appear as a result of flux transport to the poles and the ensuing strong dipole field (Cranmer 2009, and references therein). Large fragmented coronal holes with diffuse boundaries (i.e., “patchy” coronal holes) also often appear at low latitudes during solar minimum (Heinemann et al. 2020; Samara et al. 2022).
To understand and model the coronal structure in a time-dependent fashion, the evolution of coronal holes must be analyzed. The evolution of the area of long-lived coronal holes typically exhibits a pattern of growing to a maximum before decaying over multiple solar rotations with lifetimes from a few solar rotations up to more than two years (Heinemann et al. 2018b, 2020; Hewins et al. 2020). Wang & Sheeley (1990) and Gosling (1996) argue that changes in the underlying photospheric magnetic field affect the appearance of coronal holes. Heinemann et al. (2020), however, showed that there is no systematic correlation between the evolution of the underlying magnetic field and the area of the coronal hole, suggesting that the primary driver of the coronal hole’s evolution is the interaction with the local magnetic field configuration surrounding the coronal hole and the global magnetic field.
On small scales the coronal hole boundary changes dynamically through interchange reconnection, also called “footpoint switching”, where surrounding closed fields interact with the coronal hole’s open fields. This reconfiguration may cause a coronal hole to grow or decay. The rate of this process seems to be largely reliant on the surrounding fields, specifically the abundance of large-scale arcades and nearby active regions (see, e.g., Shelke & Pande 1984; Wang & Sheeley 2004; Madjarska et al. 2004; Madjarska & Wiegelmann 2009; Yang et al. 2011; Ma et al. 2014; Kong et al. 2018). Studies show that coronal holes evolve continuously (as long as external drivers, such as a rapid global reconfiguration of the coronal magnetic field, are not apparent), which suggests that interchange reconnection processes occur regularly at the coronal hole boundaries (Bohlin 1977; Kahler & Moses 1990). Interchange reconnection due to emerging flux within the coronal hole may, however, further aid in the decay of coronal holes (Zhang et al. 2007). Recently, interchange reconnection has also been suggested as a driver of the outflowing fast solar wind (Bale et al. 2023).
Coronal holes may also change their shape over the course of their lifetime. Although often considered to be rigid structures (Timothy et al. 1975), studies have shown that solar differential rotation can have a significant effect on the coronal manifestation of a coronal hole (Caplan et al. 2016; Heinemann et al. 2020). Wang & Sheeley (1993) proposed that to maintain a force-free state of open field structures that makes up the coronal holes, the rotation needs to be quasi-rigid to avoid field line twisting, but they neglected interchange reconnection. Although many coronal holes maintain their rough shape for some time, solar differential rotation often results in an apparent counterclockwise (in the northern hemisphere) or clockwise (in the southern hemisphere) tilting motion of coronal holes with a north–south extension. Prominent examples are the 1996 Elephant’s Trunk coronal hole (as studied by Del Zanna & Bromage 1999, and others) and the well-studied 2012 long-lived coronal hole (described by Heinemann et al. 2018a,b).
While coronal holes have been intensively studied, the evolution of the connection between the magnetic field in the photosphere and the resulting observed coronal hole morphology is not yet understood. Usually, coronal holes are treated as force-free and potential open fields rooted in the photosphere, and yet studies show significant differences in modeled open fields and coronal holes observed in the EUV (Asvestari et al. 2019; Linker et al. 2021). For this work we analyzed the evolution of one well-observed on-disk small coronal hole, using EUV images and magnetograms as well as magnetic field extrapolations. We observed an especially rapid rate of coronal hole tilting motion that has not yet been reported in the literature. This rapid evolution is inconsistent with changes caused by the motion of magnetic flux elements in the photosphere, and exceeds the effects of differential rotation. Our results provide insight into the connection between the photospheric and coronal field signatures in open-field regions, and highlight processes that may be involved in coronal hole evolution and the possible subsequent restructuring of the global field. In Sect. 2 we present the coronal hole and its evolution over 12 days. In Sect. 3 we explore the results, and propose two possible physical mechanisms that may be involved in the apparent accelerated tilting motion. The results are then summarized in Sect. 4.
2. Observational results
2.1. Data and preparation
In this study we use data from the Extreme UltraViolet Imager (EUVI; Wuelser et al. 2004) and the Atmospheric Imaging Assembly (AIA; Lemen et al. 2012), which are EUV imagers on board the Solar TErrestrial RElations Observatories-Ahead (STEREO-A; Kaiser et al. 2008) and the Solar Dynamics Observatory (SDO; Pesnell et al. 2012) spacecrafts. Since STEREO-A was situated ∼53 deg east of SDO, combining the two datasets increased the observational coverage of the corona during the time range of interest. For EUVI and AIA, we used the 195 Å and 193 Å channels because of the reduced emissivity of coronal holes in contrast to the surrounding corona in these narrow bands. This leads to a significant contrast. The primary contribution to these filters is the coronal emission line of 11 times ionized iron (Fe XII), and they are commonly used for coronal hole extraction (e.g., Boucheron et al. 2016; Hofmeister et al. 2017; Heinemann et al. 2019).
The EUV images used for the analysis and the coronal hole boundaries were prepared and extracted using the Coronal Hole Mapping and Analysis Pipeline (CHMAP) open-source Python software package1. CHMAP integrates all available instruments, in this case AIA/SDO and EUVI/STEREO-A into a single synchronic Carrington map with coronal hole detections. Prior to mapping, each image undergoes several preparation steps: (1) GPU-enabled PSF deconvolution, (2) data-driven limb brightening correction, and (3) inter-instrument transformation. The deconvolution step is instrument specific and tends to sharpen features in the disk image and removes a portion of the long-distance scattered light, increasing the dynamic contrast in the images. The limb brightening correction uses the equatorial band of images averaged over a six-month period to fit a correction factor as a function of center-to-limb angle and pixel intensity. This correction has the effect of “flattening” image intensity from center to limb, while leaving active regions and holes intact. The third correction is a linear transform between instruments. This is again calculated from an equatorial band that is time-averaged over six months and brings images from the two instruments to a similar overall brightness. Once the images from different instruments are brought to a consistent format, coronal hole detections are performed by an iterative, kernel-based scheme with one lower intensity threshold for seeding and a second higher threshold for boundary definitions
(Caplan et al. 2016). Thus, a coronal hole must have at least one pixel less than the lower threshold. The chosen values are shown in Table 1. The original detection parameters in Caplan et al. (2016) were trained on qualitative expert evaluations for the period of 2012–2014. For this study, the seeding parameter was increased from 0.95 to 1.105 due to increased overall brightness of the corrected images relative to the training period causing some apparent coronal holes to be missed. The algorithm then iteratively grows the hole from these seeds until the kernel and/or upper threshold are met at all coronal hole edges. The detection is performed at full disk-image resolution, and the images are then interpolated to Carrington maps (4096 × 2048 in longitude and latitude, respectively). For each time step we now have a map for each imaging instrument. These maps are combined into a single synchronic map by a minimum-intensity-merge procedure with consideration for viewing angles. A map pixel that is on the limb for AIA, but nearly at disk center for STEREO-A, will always have the STEREO-A value. Following the merge, the map resolution was reduced by integration to 1080 × 540, which corresponds to a spatial resolution of 0.333 degrees. We used EUV images prepared by the CHMAP pipeline from April 23, 2021 (when the coronal hole rotated into STEREO-A’s field-of-view), to May 5, 2021 (when the coronal hole rotated out of SDO’s field of view) at a five-minute cadence.
CHMAP extraction parameters.
During the time period when the coronal hole was visible in SDO observations, we used 720 s line-of-sight (LOS) observations from the Helioseismic and Magnetic Imager (HMI; Schou et al. 2012) on board SDO to derive the underlying magnetic flux density and its signed and unsigned magnetic flux. The magnetic field data were projected to heliographic coordinates and corrected for solar curvature assuming radial fields, following Hofmeister et al. (2017) and Heinemann et al. (2019).
2.2. A tilting coronal hole
Figure 1 shows snapshots of the evolution of the coronal hole over 12 days. Even without sophisticated analysis, it is clear that the coronal hole undergoes significant changes during this period of time: The coronal hole displays a strong clockwise tilting motion with time over the duration of the observations and a decrease in its area. Although the complete dissolution of the coronal hole cannot be fully observed because it rotates out of SDOs field of view, we note that (1) the coronal hole is barely visible in the EUV after May 5, while it is still in the field of view of SDO (even when taking into account projection effects) and (2) no trace of it is found when the corresponding Carrington longitude rotates into the STEREO-A field of view on May 19. In addition to the morphological evolution of the coronal hole itself, we observe changes in the surrounding corona: the decay of an ephemeral region at the northwestern edge (around 155° Carrington longitude and −6° latitude) that was present at the start of the observations, and the emergence of a second ephemeral region close-by a few days later (around 142° Carrington longitude and −8° latitude). The evolution of this region is described in Diercke et al. (in prep.).
![]() |
Fig. 1. Daily EUV images featuring the coronal hole. The consecutive panels show EUV images that highlight the evolution of the coronal hole between April 24 and May 5, 2021. The images were prepared by merging STEREO-A/EUVI 195 Å and SDO/AIA 193 Å filtergrams in the CHMAP pipeline. The coronal hole is shown in the heliographic Carrington frame on a 10° grid (white guidelines) from −50° to −10° latitude and 120°–160° longitude. |
2.3. Morphological evolution
Using the extracted boundaries of the coronal hole, we calculated a tilting rate of the coronal hole by defining the tilt as the angle (α) between the major axis of the coronal hole and the north–south axis of the Sun. We determined the major axis by calculating a longitudinal geometric center of mass of the coronal hole extraction at every latitude and then applying a linear fit. An example is shown in Fig. 2a. From this the tilting rate of the coronal hole is defined as
![]() |
Fig. 2. Tilt angle calculation. Panel a shows an example of the tilt angle calculation for the coronal hole on April 30, 2021, at 23:50. The red crosses indicate the latitudinal geometric centers of mass, which are then fitted (red line). The angle between the fit and the north–south direction (blue line) is defined as the tilt angle of the coronal hole. In panel b a schematic of the tilt angle calculation solely caused by differential rotation is shown. The light blue arrows highlight the solar differential rotation rates at different latitudes (not to scale). |
To compare the observed tilting rate of the coronal hole to the differential rotation, we considered an idealized vertical coronal hole (approximated by a line ≡ tilt axis) with a given latitudinal extent (h). Due to differential rotation, the axis of such a coronal hole should tilt approximately according to the difference in rotation rate between the lower and upper part. The difference in rotation rate is given by
with ω⊙ denoting the differential rotation rate at a given latitude (λ). With this relative rotation rate between the lower and upper part of the theoretical coronal hole, we can calculate the relative longitudinal shift dr as a function of time (t)
and subsequently the tilt angle
This is schematically shown in Fig. 2b.
To compare our coronal hole tilting rate to an established photospheric rotation rate we use the differential rotation profile given by Snodgrass & Ulrich (1990), derived from the magnetic feature correlation of the Mount Wilson magnetograms:
We approximate the coronal differential rotation by using the profile for coronal hole rotation rates in the southern hemisphere derived by Bagashvili et al. (2017):
By assuming a longitudinal extent of 20° stretching from latitudes λmin = −20° to λmax = −40°, we can calculate the tilt angle of a hypothetical coronal hole under the influence of either the photospheric or the coronal differential rotation rate. By substituting the differential rotation profiles (Eqs. (5) and (6), respectively) into Eq. (2), we obtain
Using Eqs. (3) and (4), we derive the respective tilt angles (α(ωphot) and α(ωch)) as a function of time. These are shown in Fig. 3 as the orange and green lines. The tilting rate can then be derived using Eq. (1), giving hypothetical coronal hole tilting rates solely caused by advection due to differential rotation of and
.
![]() |
Fig. 3. Coronal hole tilt angle (α) as a function of time from April 24 to May 5. The three-hour running mean is highlighted in black ( |
Figure 3 shows the coronal hole tilt angle as a function of time between April 24 and May 5. On average, we find that the coronal hole shows an apparent clockwise tilting motion, as expected for a north–south elongated coronal hole located in the southern hemisphere (Heinemann et al. 2018b). Between April 24 and May 1 the coronal hole tilts at a mean rate of ∼3.2° day−1, although periods of fast counterclockwise followed by immediate fast clockwise tilting are present (indicated in Fig. 3 by the red arrows). After visual inspection of the coronal hole boundaries (see Fig. 1), we conclude that these short-term changes of the tilt angle are caused by localized short-term changes in the coronal hole boundary and ambiguity in the determination of the tilt angle rather than the long-term evolution of the coronal hole. Interestingly, between May 2 and May 4 the tilt angle seems to change at a higher rate of ∼5.4° day−1 (yellow line), but in contrast to the earlier fast changes, it does not seem to be caused by short-term boundary displacement, but rather by an apparent tilting motion of the whole coronal hole structure.
We find that the mean tilting rate of the coronal hole (∼3.2° day−1) between April 24 and May 1 is slightly higher than the photospheric rate (2.6° day−1), and significantly higher than the coronal rate (1.2° day−1). Between May 2 and May 4, where the coronal hole tilts rapidly, the rate of ∼5.4° day−1 is more than double the photospheric rate and over four times the coronal rate.
Figure 4 shows the coronal hole area as a function of time. We observe that the area changes significantly between April 24 and April 26 by increasing in size by over 50% from 10 × 109 km2 to 15 × 109 km2 and then decaying back to 10 × 109 km2 within a day. Subsequently, until around April 30 the area seems to remain constant, even though the data in the time period between April 28 21:00 and April 29 15:00 is unreliable because SDO was pointing off-limb, due to a special observation campaign and the STEREO-A data showing strong projection effects. Starting around April 29 06:00, we observe a steady decline to ∼4 × 109 km2 until May 4 . Finally, we find a slight increase in the area up to ∼5 × 109 km2 in the hours before the coronal hole rotated out of SDO’s field of view.
![]() |
Fig. 4. Coronal hole area (A) as a function of time from April 24 to May 5. The three-hour running mean is highlighted in black ( |
2.4. Magnetic evolution
Here we present the evolution of the coronal hole’s magnetic properties during the time period when the coronal hole was observed by SDO. In Fig. 5 the area, magnetic flux density, and magnetic flux are shown for the time period April 30 to May 5. The magnetic field properties were calculated by integrating each magnetogram over the projected coronal hole boundary and correcting for radial fields (for details, see Heinemann et al. 2019). The area decreases continuously from ∼12.8 × 109 km2 on April 29 22:00 to ∼3.9 × 109 km2 on May 3 22:00 (Fig. 5, top panel). The area-averaged (signed) magnetic flux density of the underlying photospheric field within the projected coronal hole area (Fig. 5, middle panel) increases slowly from 1.46 ± 0.32 G on April 29 to 2.26 ± 0.42 G on May 3 (the time-averaged value was found to be 2.03 ± 0.47 G). The respective signed (or open) flux remains approximately constant over most of the observed time period between April 29 and May 3 at (1.82 ± 0.43)×1020 Mx.
![]() |
Fig. 5. Coronal hole properties as functions of time from April 29 15:00 to May 4 00:00. The coronal hole area (A; red), signed magnetic flux density (Bs;green), unsigned magnetic flux (Φus; light blue), and signed magnetic flux (Φs; blue) are shown from top to bottom with a three-hour running mean highlighted in black ( |
Figure 3 shows a tilting rate of the coronal hole (up to May 1) that exceeds the rate expected from a coronal hole being advected by the photospheric differential rotation. This disparity between the rotation rate of the coronal hole and photosphere manifests as magnetic elements rotating out of the projected eastern (left) coronal hole boundary. In Fig. 6 we show the magnetic field evolution of a subregion of the coronal hole over a few days. The overlaid coronal hole boundary highlights the respective area evolution in that subfield on the coronal hole’s east side. In the first time step shown in panel a a small bipolar region, which corresponds to a coronal bright point in AIA 193 Å, is excluded from the extracted coronal hole boundary. This bright point cannot be observed in the next image shown, likely due to flux cancellation with the opposite polarity region next to it within less than 14 h, and only a large unipolar magnetic element remains as a footpoint of the open field and part of the coronal hole. Over the next three days this specific magnetic element moves eastward relative to the boundary and can be seen located outside the coronal hole boundary at later times. This is consistent with a motion of the coronal hole, faster than the photospheric differential rotation. Another possible explanation is that the apparent westward motion of the coronal hole in this subfield might be caused by the evolution (in this case closing or shrinking) of the boundary. We find that other long-lived unipolar magnetic elements also show an eastward motion relative to the overlying coronal hole.
![]() |
Fig. 6. Evolution of the LOS magnetic field underlying the coronal hole at a cadence of 14 h. A subregion is highlighted to show the motion of a long-lived unipolar magnetic element (blue contour). The coronal hole boundary derived from EUV observations is shown in red. The green contour in panel a gives the location of a coronal bright point within the coronal hole (which can be seen in 193 Å), which no longer exists 14 h later. |
2.5. Topological evolution in magnetic field extrapolations
Coronal holes are commonly defined by their reduced emission in the corona and their open magnetic field configuration. They are often reconstructed as open field areas derived from magnetic field modeling. The most widespread approach is to consider the coronal hole approximately force free and potential. This approach uses either potential field source surface (PFSS; Altschuler & Newkirk 1969) models or adaptations, such as the Wang–Sheeley–Arge model (Arge & Pizzo 2000) that includes a Schatten current sheet (Schatten et al. 1968). To evaluate the conformity of the open fields derived from a potential magnetic field model with the observed coronal hole and to determine if the modeled open field region shows an evolutionary behavior similar to the observed structure, we model multiple instances over time using a global PFSS model. We employ a finite difference scheme (Stansby et al. 2020) with synchronic HMI magnetograms (hmi.mrdailysynframe_polfil_720s).
In Fig. 7 we present the modeled open fields calculated for different source surface heights Rss (i.e., the height at which the field lines are forced to be radial). We found, that for the standard height for the source surface (Rss = 2.5 R⊙), barely any open field areas are reconstructed. For that source surface height and from April 27 to April 30, only a tiny patch of open field can be seen near the coronal hole’s northern edge. Following this time period, no open fields are produced for this height. When reducing the value for the source surface height to around Rss = 2.1 R⊙, we find some open fields roughly at the location of the coronal hole for all times considered. The shape and precise location of these modeled open fields, however, do not match the observed coronal hole shape. The surrounding modeled closed field structure also does not exhibit any sudden or large-scale reconfiguration during this time period (not shown here), and although it is common to observe some disparity between observed coronal hole areas and modeled open fields as a function of source surface height and in terms of coronal hole shape (Asvestari et al. 2019, 2020; Linker et al. 2021), our model results lack any evidence of the evolution observed in EUV. There are no changes in the calculated open field area that resemble the changes in the observed morphology nor the observed apparent shear or tilting motion of the coronal hole presented in Figs. 3 and 5. We note that the choice of magnetogram may influence the resulting open field regions (e.g., Caplan et al. 2021; Asvestari et al. 2023). To account for this, we calculated the PFSS open field regions for multiple magnetograms (using the GONG and ADAPT maps as well as different resolutions), but found none that matched the observations to a significantly better degree. However, this does not exclude the possibility that this discrepancy may stem from inadequate boundary conditions in the model.
![]() |
Fig. 7. Coronal hole area overlaid with open field footpoints. The open fields were calculated at 1.05 R⊙ using PFSS modeling and different source surface heights (RSS). The colored contours represent the open field area for respective source surface heights as noted on the colorbar. |
2.6. Dynamic evolution
To investigate the dynamics of coronal hole evolution, we looked for small-scale short-term EUV brightenings that are commonly regarded as tracers for reconnection events (e.g., Shokri et al. 2022). In our dataset we defined EUV brightenings as pixels where the intensity increases by more than 100% in one time step and whose enhancement does not appear to last longer than 10 min. In addition, the intensity enhancement must be 3σ above the mean coronal hole intensity to avoid the inclusion of noise-dominated changes.
In Fig. 8 we indicate the location of all brightenings detected in and around the coronal hole, color-coded as a function of time. In each panel all tracer locations detected within ±12 h of the times when the shown EUV image are displayed. We find that these brightenings occur continuously throughout the observed time period. Preferred locations are where changes in the coronal hole are seen in subsequent images. From around May 3 to May 5 we find enhancements in the northern part of the coronal hole in close proximity to the emerging ephemeral region. We detected 205 such brightenings in the observed time period. Between April 24 and 26 we find an occurrence rate of just above one per hour, while afterwards the rate drops to one every two to three hours. We suspect that the higher initial rate of brightenings (during times that are dominated by STEREO-A observerations) might be due to some instrumental differences between STEREO and SDO that remains even after the inter-instrument calibration.
![]() |
Fig. 8. Daily images of the evolution of the coronal hole from April 24, 2021, to May 5, 2021. The images show merged STEREO-A/EUVI 195 Å and SDO/AIA 193 Å filtergrams as prepared by the CHMAP pipeline. The coronal hole is shown in the heliographic Carrington frame on a 10° grid (white guidelines) from −50° to −10° latitude and 120°–160° longitude. In each image the locations of short-term EUV brightenings that occurred in the time period of ±12 h around the shown frame are labeled “x”. The labeled locations are color-coded according to the time of detection. |
3. Discussion
In Sect. 2 we presented the evolution of a coronal hole over 12 days. Until May 1 the coronal hole shows a clockwise tilting motion that appears to develop at a rate slightly higher than the expected value if the sole source of the tilting was advection of its photospheric footpoints according to a differential rotation profile. Around May 2 to May 4 the coronal hole starts tilting rapidly at nearly double the previous rate. In addition, during the same time period from April 30 to May 3 the area decreases by more than a factor of three within a few days, while at the same time the majority of the coronal hole’s open flux is conserved.
It is believed that the open field that is seen as a dark region in the solar corona is rooted in unipolar magnetic elements in the photosphere (Hofmeister et al. 2019). According to the general understanding of the expansion of these open fields (Wedemeyer-Böhm et al. 2009; Cranmer & Winebarger 2019), it can be assumed that their coronal representation (i.e., the coronal hole) is formed by the sum of the open field lines extending from their footpoints and the interaction with the surrounding and global fields. It is well known that the photosphere, and thus the magnetic elements on large scales, rotate differentially, which would suggest that a coronal hole should rotate in the same manner as those magnetic elements if an inflexible connection between footpoints and the corona is assumed. However, this contradicts the observed rigid behavior (e.g., Timothy et al. 1975) of coronal holes not rotating entirely differentially. Heinemann et al. (2018b) found a tilting rate for one specific coronal hole in 2012 of 0.42° day−1, which is at least a factor of four smaller than the photospheric rotation rate for the same latitudinal range due to differential rotation (as given by Snodgrass & Ulrich 1990), and around half of the rotation rate of coronal holes derived by Bagashvili et al. (2017). A rotation rate slower than differential rotation could be caused by the faster rotating underlying photosphere applying a twist on the more uniform magnetic coronal structure causing footpoint switching near the boundary to retain the coronal hole’s integrity (e.g., Kong et al. 2018).
In our study we show first that the tilting rate appears to exceed the rate predicted by a coronal hole solely advected by a photospheric differential rotation profile. The excess tilting is even greater when comparing to a coronal rotation profile (see Fig. 3). Second, we show that the suspected coronal hole footpoints (unipolar magnetic elements) move eastward in comparison to the coronal hole boundary, which indicates that parts of the coronal hole (as seen in the corona) rotate faster than the photospheric differential rotation (see Fig. 6). Finally, we find that the PFSS modeled open-field structure does not show evolutionary behavior that is similar to the structure observed in EUV. The noticeable tilting motion cannot be seen in the modeled open fields (see Fig. 7). From these three observations, it follows that the rate at which the coronal hole tilts cannot be explained by photospheric and/or coronal differential rotation acting on the coronal hole. Thus, it becomes necessary to identify alternative physical mechanisms responsible for this phenomenon. We propose two possibilities:
First, we cannot discern whether the strong apparent tilting motion, starting around May 2, was caused by coronal hole substructures evolving individually. The evolution of the northernmost parts of the coronal hole seems to be strongly influenced by the emergence of a nearby ephemeral active region. The sudden flux emergence near the boundary may have led to a westward shifting of a part of the coronal hole. This can occur due to reconnection to an energetically lower state or due to being pushed by magnetic pressure. Terradas (2023) showed such a case using 2D magnetostatic simulations.
Second, we suspect that the continuous interchange reconnection in and around the coronal hole (Madjarska & Wiegelmann 2009) caused dynamic changes in the open field configuration which could account for the observed apparent tilt motion. Indicators for this are the magnetic elements that are seen moving in and out of the coronal hole in an eastward direction and that the open flux in the coronal hole remains roughly constant until after May 3, although the area has decreased by approximately a factor of three. This suggests that the majority of the morphological evolution and area decay is caused by a flux conserving process, likely interchange reconnection in the corona, rather than changes in the underlying photospheric magnetic field structure due to flux emergence and cancellation. We support this by showing that localized short-term EUV brightenings can be seen during the coronal hole evolution (see Fig. 8). Those brightenings are primarily located around areas that show subsequent morphological changes. They may indicate coronal transient events observed at these spatial scales, such as coronal jets (Sterling et al. 2015; Panesar et al. 2018), coronal bright points (Habbal & Withbroe 1981; Matković et al. 2023), or blinkers (Subramanian et al. 2012; Shokri et al. 2022); furthermore, it is generally acknowledged that reconnection plays a major role in their formation (e.g., Shokri et al. 2022). Additionally, a coronal hole features a primarily open-close magnetic configuration with lower-lying closed loops and open fields stretching out into the heliosphere (e.g., Wiegelmann & Solanki 2004; Bale et al. 2023). Thus, it is reasonable to speculate that these detected brightenings are connected to interchange reconnection with associated small-scale restructuring of the open-closed field.
4. Conclusions
In this study we presented the evolution of a well-observed coronal hole over 12 consecutive days. Using EUV and LOS magnetic field observations in combination with PFSS modeling, we tracked and analyzed the morphological and magnetic evolution of the coronal hole from April 23 to May 5, 2021. Our findings can be summarized as follows:
-
For the first time, we describe a coronal hole that tilts at a rate exceeding the rate expected from differential rotation. Our coronal hole tilts at a mean rate (3.2° day−1) that exceeds the expected effect of solar plasma being advected by either photospheric (2.6° day−1) or coronal hole differential rotation (1.2° day−1).
-
In the time period from May 1 to May 3 we observed an acceleration of the tilting rate to 5.4° day−1.
-
While the area of the coronal hole decays by more than a factor of three in approximately four days (from ∼12.8 × 109 km2 to ∼3.9 × 109 km2), its signed flux remains conserved at a value of (1.82 ± 0.43)×1020 Mx.
-
We find that photospheric unipolar magnetic elements move eastward relative to the overlying coronal hole boundary. This hints that parts of the coronal hole rotate (or evolve) faster than the underlying photosphere.
-
The prominent tilting motion of the coronal hole cannot be reproduced by potential field extrapolations of the coronal hole’s open field structure.
-
Throughout the observed time period we find small-scale brightenings occurring in and around the coronal hole at a rate of 0.3–1 per hour. These are suspected to be tracers for interchange reconnection.
-
We suggest that the rapid morphological evolution of the coronal hole may be primarily driven by interchange reconnection at the boundary (which is a flux conserving process) and the interaction with a newly emerging ephemeral region in the northern part of the coronal hole.
This coronal hole highlights the dynamic nature of the solar corona even in regions traditionally believed to be rigid and stable. The dynamic evolution may have a significant effect on the outflowing solar wind and could also be a reason for the poor agreement between EUV observations and PFSS modeled open fields.
Acknowledgments
This paper resulted from discussions with the team of the extensive multi-instrument observation campaign lead by S. J. Hofmeister from April 28th 2021. Thus, the authors would like to thank P. Bryans, M. Cantoresi, D. J. Christian, M. Globa, M. Hahn, N. Huang, C. Kuckein, D. Lacatus, G. de Toma, K. Tziotziou, E. Samara, C. Vocks and D. W. Savin for fruitful discussions and helpful comments. The authors would like to thank J. Linker and R. Caplan from Predictive Science Inc. for providing the coronal hole extractions. S.G.H. expresses gratitude to K. Horaites for engaging in inspirational listening. This research was funded in whole, or in part, by the Austrian Science Fund (FWF) Erwin-Schrödinger fellowship J-4560. S.J.H. acknowledges support by the German Science Fund (DFG), grant No. 448336908. A.C.S. was supported with funding from the Heliophysics Division of NASA’s Science Mission Directorate through the Heliophysics Supporting Research (HSR, grant No. 20-HSR20 2-0124) Program, and through the Heliophysics System Observatory Connect (HSOC, grant No. 80NSSC20K1285) Program. J.A.T. and C.D. acknowledge support in the creation of CHMAP from the NASA Guest Investigators Program (NNX17AB78G). J.P. acknowledges the Academy of Finland Project 343581. E.A. acknowledges support from the Academy of Finland (Postdoctoral Researcher Grant 322455) and the ERC under the European Union’s Horizon 2020 Research and Innovation Programme Project 724391 (SolMAG). The authors acknowlege support from NASA Heliophysics Living with a Star grant 80NSSC20K0183. The SDO and STEREO image data are available by courtesy of NASA and the respective science teams.
References
- Altschuler, M. D., & Newkirk, G. 1969, Sol. Phys., 9, 131 [NASA ADS] [CrossRef] [Google Scholar]
- Arge, C. N., & Pizzo, V. J. 2000, J. Geophys. Res., 105, 10465 [Google Scholar]
- Asvestari, E., Heinemann, S. G., Temmer, M., et al. 2019, J. Geophys. Res. (Space Phys.), 124, 8280 [NASA ADS] [CrossRef] [Google Scholar]
- Asvestari, E., Heinemann, S. G., Temmer, M., et al. 2020, J. Phys. Conf. Ser., 1548, 012004 [NASA ADS] [CrossRef] [Google Scholar]
- Asvestari, E., Arge, C. N., Caplan, R. M., et al. 2023, ApJ, submitted [Google Scholar]
- Bagashvili, S. R., Shergelashvili, B. M., Japaridze, D. R., et al. 2017, A&A, 603, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bale, S. D., Drake, J. F., McManus, M. D., et al. 2023, Nature, 618, 252 [NASA ADS] [CrossRef] [Google Scholar]
- Bohlin, J. D. 1977, Sol. Phys., 51, 377 [NASA ADS] [CrossRef] [Google Scholar]
- Boucheron, L. E., Valluri, M., & McAteer, R. T. J. 2016, Sol. Phys., 291, 2353 [NASA ADS] [CrossRef] [Google Scholar]
- Caplan, R. M., Downs, C., & Linker, J. A. 2016, ApJ, 823, 53 [CrossRef] [Google Scholar]
- Caplan, R. M., Downs, C., Linker, J. A., & Mikic, Z. 2021, ApJ, 915, 44 [NASA ADS] [CrossRef] [Google Scholar]
- Cranmer, S. R. 2009, Liv. Rev. Sol. Phys., 6, 3 [Google Scholar]
- Cranmer, S. R., & Winebarger, A. R. 2019, ARA&A, 57, 157 [Google Scholar]
- Del Zanna, G., & Bromage, B. J. I. 1999, J. Geophys. Res., 104, 9753 [NASA ADS] [CrossRef] [Google Scholar]
- Fludra, A., Del Zanna, G., & Bromage, B. J. I. 1999, Space. Sci. Rev., 87, 185 [NASA ADS] [CrossRef] [Google Scholar]
- Garton, T. M., Murray, S. A., & Gallagher, P. T. 2018, ApJ, 869, L12 [NASA ADS] [CrossRef] [Google Scholar]
- Geyer, P., Temmer, M., Guo, J., & Heinemann, S. G. 2021, A&A, 649, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gosling, J. T. 1996, ARA&A, 34, 35 [NASA ADS] [CrossRef] [Google Scholar]
- Habbal, S. R., & Withbroe, G. L. 1981, Sol. Phys., 69, 77 [Google Scholar]
- Hahn, M., Landi, E., & Savin, D. W. 2011, ApJ, 736, 101 [NASA ADS] [CrossRef] [Google Scholar]
- Heinemann, S. G., Hofmeister, S. J., Veronig, A. M., & Temmer, M. 2018a, ApJ, 863, 29 [NASA ADS] [CrossRef] [Google Scholar]
- Heinemann, S. G., Temmer, M., Hofmeister, S. J., Veronig, A. M., & Vennerstrøm, S. 2018b, ApJ, 861, 151 [NASA ADS] [CrossRef] [Google Scholar]
- Heinemann, S. G., Temmer, M., Heinemann, N., et al. 2019, Sol. Phys., 294, 144 [NASA ADS] [CrossRef] [Google Scholar]
- Heinemann, S. G., Jerčić, V., Temmer, M., et al. 2020, A&A, 638, A68 [EDP Sciences] [Google Scholar]
- Heinemann, S. G., Saqri, J., Veronig, A. M., Hofmeister, S. J., & Temmer, M. 2021, Sol. Phys., 296, 18 [NASA ADS] [CrossRef] [Google Scholar]
- Hewins, I. M., Gibson, S. E., Webb, D. F., et al. 2020, Sol. Phys., 295, 161 [Google Scholar]
- Hofmeister, S. J., Veronig, A., Reiss, M. A., et al. 2017, ApJ, 835, 268 [Google Scholar]
- Hofmeister, S. J., Utz, D., Heinemann, S. G., Veronig, A., & Temmer, M. 2019, A&A, 629, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hofmeister, S. J., Veronig, A. M., Poedts, S., Samara, E., & Magdalenic, J. 2020, ApJ, 897, L17 [NASA ADS] [CrossRef] [Google Scholar]
- Hofmeister, S. J., Asvestari, E., Guo, J., et al. 2022, A&A, 659, A190 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kahler, S. W., & Moses, D. 1990, ApJ, 362, 728 [NASA ADS] [CrossRef] [Google Scholar]
- Kaiser, M. L., Kucera, T. A., Davila, J. M., et al. 2008, Space. Sci. Rev., 136, 5 [Google Scholar]
- Karachik, N. V., Pevtsov, A. A., & Abramenko, V. I. 2010, ApJ, 714, 1672 [NASA ADS] [CrossRef] [Google Scholar]
- Kong, D. F., Pan, G. M., Yan, X. L., Wang, J. C., & Li, Q. L. 2018, ApJ, 863, L22 [Google Scholar]
- Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17 [Google Scholar]
- Linker, J. A., Heinemann, S. G., Temmer, M., et al. 2021, ApJ, 918, 21 [CrossRef] [Google Scholar]
- Ma, L., Qu, Z.-Q., Yan, X.-L., & Xue, Z.-K. 2014, RAA, 14, 221 [NASA ADS] [Google Scholar]
- Madjarska, M. S., & Wiegelmann, T. 2009, A&A, 503, 991 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Madjarska, M. S., Doyle, J. G., & van Driel-Gesztelyi, L. 2004, ApJ, 603, L57 [NASA ADS] [CrossRef] [Google Scholar]
- Matković, F., Brajša, R., Temmer, M., et al. 2023, A&A, 670, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Panesar, N. K., Sterling, A. C., & Moore, R. L. 2018, ApJ, 853, 189 [NASA ADS] [CrossRef] [Google Scholar]
- Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, Sol. Phys., 275, 3 [Google Scholar]
- Petrie, G. J. D., & Haislmaier, K. J. 2013, ApJ, 775, 100 [NASA ADS] [CrossRef] [Google Scholar]
- Samara, E., Magdalenić, J., Rodriguez, L., et al. 2022, A&A, 662, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Saqri, J., Veronig, A. M., Heinemann, S. G., et al. 2020, Sol. Phys., 295, 6 [Google Scholar]
- Schatten, K. H., Ness, N. F., & Wilcox, J. M. 1968, Sol. Phys., 5, 240 [NASA ADS] [CrossRef] [Google Scholar]
- Schou, J., Scherrer, P. H., Bush, R. I., et al. 2012, Sol. Phys., 275, 229 [Google Scholar]
- Shelke, R. N., & Pande, M. C. 1984, Bull. Astron. Soc. India, 12, 404 [NASA ADS] [Google Scholar]
- Shokri, Z., Alipour, N., Safari, H., et al. 2022, ApJ, 926, 42 [NASA ADS] [CrossRef] [Google Scholar]
- Snodgrass, H. B., & Ulrich, R. K. 1990, ApJ, 351, 309 [NASA ADS] [CrossRef] [Google Scholar]
- Stansby, D., Yeates, A., & Badman, S. T. 2020, J. Open Source Softw., 5, 2732 [Google Scholar]
- Sterling, A. C., Moore, R. L., Falconer, D. A., & Adams, M. 2015, Nature, 523, 437 [NASA ADS] [CrossRef] [Google Scholar]
- Subramanian, S., Madjarska, M. S., Doyle, J. G., & Bewsher, D. 2012, A&A, 538, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Terradas, J. 2023, Physics, 5, 276 [NASA ADS] [CrossRef] [Google Scholar]
- Timothy, A. F., Krieger, A. S., & Vaiana, G. S. 1975, Sol. Phys., 42, 135 [NASA ADS] [CrossRef] [Google Scholar]
- Tu, C.-Y., Zhou, C., Marsch, E., et al. 2005, Science, 308, 519 [Google Scholar]
- Wang, Y.-M., & Sheeley, N. R., Jr. 1990, ApJ, 355, 726 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, Y.-M., & Sheeley, N. R., Jr. 1993, ApJ, 414, 916 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, Y.-M., & Sheeley, N. R., Jr. 2004, ApJ, 612, 1196 [NASA ADS] [CrossRef] [Google Scholar]
- Warren, H. P., & Hassler, D. M. 1999, J. Geophys. Res., 104, 9781 [NASA ADS] [CrossRef] [Google Scholar]
- Wedemeyer-Böhm, S., Lagg, A., & Nordlund, Å. 2009, Space. Sci. Rev., 144, 317 [Google Scholar]
- Wiegelmann, T., & Solanki, S. K. 2004, Sol. Phys., 225, 227 [NASA ADS] [CrossRef] [Google Scholar]
- Wuelser, J. P., Lemen, J. R., Tarbell, T. D., et al. 2004, in Telescopes and Instrumentation for Solar Astrophysics, eds. S. Fineschi, & M. A. Gummin, SPIE Conf. Ser, 5171, 111 [NASA ADS] [CrossRef] [Google Scholar]
- Yang, S., Zhang, J., Li, T., & Liu, Y. 2011, ApJ, 732, L7 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, J., Zhou, G., Wang, J., & Wang, H. 2007, ApJ, 655, L113 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
![]() |
Fig. 1. Daily EUV images featuring the coronal hole. The consecutive panels show EUV images that highlight the evolution of the coronal hole between April 24 and May 5, 2021. The images were prepared by merging STEREO-A/EUVI 195 Å and SDO/AIA 193 Å filtergrams in the CHMAP pipeline. The coronal hole is shown in the heliographic Carrington frame on a 10° grid (white guidelines) from −50° to −10° latitude and 120°–160° longitude. |
In the text |
![]() |
Fig. 2. Tilt angle calculation. Panel a shows an example of the tilt angle calculation for the coronal hole on April 30, 2021, at 23:50. The red crosses indicate the latitudinal geometric centers of mass, which are then fitted (red line). The angle between the fit and the north–south direction (blue line) is defined as the tilt angle of the coronal hole. In panel b a schematic of the tilt angle calculation solely caused by differential rotation is shown. The light blue arrows highlight the solar differential rotation rates at different latitudes (not to scale). |
In the text |
![]() |
Fig. 3. Coronal hole tilt angle (α) as a function of time from April 24 to May 5. The three-hour running mean is highlighted in black ( |
In the text |
![]() |
Fig. 4. Coronal hole area (A) as a function of time from April 24 to May 5. The three-hour running mean is highlighted in black ( |
In the text |
![]() |
Fig. 5. Coronal hole properties as functions of time from April 29 15:00 to May 4 00:00. The coronal hole area (A; red), signed magnetic flux density (Bs;green), unsigned magnetic flux (Φus; light blue), and signed magnetic flux (Φs; blue) are shown from top to bottom with a three-hour running mean highlighted in black ( |
In the text |
![]() |
Fig. 6. Evolution of the LOS magnetic field underlying the coronal hole at a cadence of 14 h. A subregion is highlighted to show the motion of a long-lived unipolar magnetic element (blue contour). The coronal hole boundary derived from EUV observations is shown in red. The green contour in panel a gives the location of a coronal bright point within the coronal hole (which can be seen in 193 Å), which no longer exists 14 h later. |
In the text |
![]() |
Fig. 7. Coronal hole area overlaid with open field footpoints. The open fields were calculated at 1.05 R⊙ using PFSS modeling and different source surface heights (RSS). The colored contours represent the open field area for respective source surface heights as noted on the colorbar. |
In the text |
![]() |
Fig. 8. Daily images of the evolution of the coronal hole from April 24, 2021, to May 5, 2021. The images show merged STEREO-A/EUVI 195 Å and SDO/AIA 193 Å filtergrams as prepared by the CHMAP pipeline. The coronal hole is shown in the heliographic Carrington frame on a 10° grid (white guidelines) from −50° to −10° latitude and 120°–160° longitude. In each image the locations of short-term EUV brightenings that occurred in the time period of ±12 h around the shown frame are labeled “x”. The labeled locations are color-coded according to the time of detection. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.