A Code for Automated Tracing of Coronal Loops Approaching Visual Perception

We develop a new numerical code with automated feature extraction, customized for tracing of coronal loops, a method we call Oriented Coronal CUrved Loop Tracing (OCCULT), which for the first time breaks even with the results of visual tracing. The method used is based on oriented-directivity tracing of curvi-linear features, but in contrast to other general feature-extraction algorithms, it is customized for solar EUV and SXR images by taking advantage of the specific property that coronal loops have large curvature radii compared with their widths. We evaluate the performance of this new code by comparing the cumulative distribution of loop lengths, the median and maximum loop lengths, the completeness of detection, and the congruency of the detected features with other numerical codes and visual tracings. We find that the new code closely approaches the results of visual perception and outperforms the other existing numerical codes. This algorithm is useful for the 3D reconstruction of the geometry, motion, and oscillations of coronal loops, with single or stereoscopic spacecraft, as well as for modeling of the loop hydrodynamics and the coronal magnetic field.     

Determining the 3D Structure of the Corona Using Vertical Height Constraints on Observed Active Region Loops

The corona associated with an active region is structured by high-temperature, magnetically dominated closed and open loops. The projected 2D geometry of these loops is captured in EUV filtergrams. In this study using SDO/AIA 171 Å filtergrams, we expand our previous method to derive the 3D structure of these loops, independent of heliostereoscopy. We employ an automated loop recognition scheme (Occult-2) and fit the extracted loops with 2D cubic Bézier splines. Utilizing SDO/HMI magnetograms, we extrapolate the magnetic field to obtain simple field models within a rectangular cuboid. Using these models, we minimize the misalignment angle with respect to Bézier control points to extend the splines to 3D (Gary, Hu, and Lee 2014). The derived Bézier control points give the 3D structure of the fitted loops. We demonstrate the process by deriving the position of 3D coronal loops in three active regions (AR 11117, AR 11158, and AR 11283). The numerical minimization process converges and produces 3D curves which are consistent with the height of the loop structures when the active region is seen on the limb. From this we conclude that the method can be important in both determining estimates of the 3D magnetic field structure and determining the best magnetic model among competing advanced magnetohydrodynamics or force-free magnetic-field computer simulations.     

Automated Temperature and Emission Measure Analysis of Coronal Loops and Active Regions Observed with the Atmospheric Imaging Assembly on the Solar Dynamics Observatory (SDO/AIA)

We developed numerical codes designed for automated analysis of SDO/AIA image datasets in the six coronal filters, including: i) coalignment test between different wavelengths with measurements of the altitude of the EUV-absorbing chromosphere, ii) self-calibration by empirical correction of instrumental response functions, iii) automated generation of differential emission measure [DEM] distributions with peak-temperature maps [T _p(x,y)] and emission measure maps [EM _p(x,y)] of the full Sun or active region areas, iv) composite DEM distributions [dEM(T)/dT] of active regions or subareas, v) automated detection of coronal loops, and vi) automated background subtraction and thermal analysis of coronal loops, which yields statistics of loop temperatures [T _e], temperature widths [σ _T], emission measures [EM], electron densities [n _e], and loop widths [w]. The combination of these numerical codes allows for automated and objective processing of numerous coronal loops. As an example, we present the results of an application to the active region NOAA 11158, observed on 15 February 2011, shortly before it produced the largest (X2.2) flare during the current solar cycle. We detect 570 loop segments at temperatures in the entire range of log(T _e)=5.7?–?7.0 K and corroborate previous TRACE and AIA results on their near-isothermality and the validity of the Rosner–Tucker–Vaiana (RTV) law at soft X-ray temperatures (T?2 MK) and its failure at lower EUV temperatures.     

Spatial and temporal variations of solar coronal loops

Skylab EUV observations of an active region near the solar limb were analyzed. Both cool (T < 10⁶ K) and hot (T > 10⁶ K) loops were observed in this region. For the hot loops the observed intensity variations were small, typically a few percent over a period of 30 min. The cool loops exhibited stronger variations, sometimes appearing and disappearing in 5 to 10 min. Most of the cool material observed in the loops appeared to be caused by the downward flow of coronal rain and by the upward ejection of chromospheric material in surges. The frequent EUV brightenings observed near the loop footpoints appear to have been produced by both in situ transient energy releases (e.g. subflares) and the infall/impact of coronal rain. The physical conditions in the loops (temperatures, densities, radiative and conducting cooling rates, cooling times) were determined. The mean energy required to balance the radiative and conductive cooling of the hot loops is approximately 3 × 10^–3 erg cm^–3 s^–1. One coronal heating mechanism that can account for the observed behavior of the EUV emission from McMath region 12634 is heating by the dissipation of fast mode MHD waves.     

2D Feature Recognition And 3d Reconstruction In Solar Euv Images

EUV images show the solar corona in a typical temperature range of T >rsim 1 MK, which encompasses the most common coronal structures: loops, filaments, and other magnetic structures in active regions, the quiet Sun, and coronal holes. Quantitative analysis increasingly demands automated 2D feature recognition and 3D reconstruction, in order to localize, track, and monitor the evolution of such coronal structures. We discuss numerical tools that “fingerprint” curvi-linear 1D features (e.g., loops and filaments). We discuss existing finger-printing algorithms, such as the brightness-gradient method, the oriented-connectivity method, stereoscopic methods, time-differencing, and space–time feature recognition. We discuss improved 2D feature recognition and 3D reconstruction techniques that make use of additional a priori constraints, using guidance from magnetic field extrapolations, curvature radii constraints, and acceleration and velocity constraints in time-dependent image sequences. Applications of these algorithms aid the analysis of SOHO/EIT, TRACE, and STEREO/SECCHI data, such as disentangling, 3D reconstruction, and hydrodynamic modeling of coronal loops, postflare loops, filaments, prominences, and 3D reconstruction of the coronal magnetic field in general.     

The coronal mass ejection associated with the loop eruption and coronal dimmings on 2009 December 13

By using the multi-wavelength observations from the Solar Terrestrial Relations Observatory (STEREO), the Solar and Heliospheric Observatory (SOHO), the Transition Region and Coronal Explorer (TRACE) and the HINODE, we study a coronal mass ejection (CME) and associated coronal dimming occurred on 2009 December 13, as a consequence of the expansion and eruption of EUV loops. The activities were probably triggered by the new flux emergence and the convergence motions, which were evident in the magnetograms from Michelson Doppler Imager (MDI) on the SOHO. The emergence led to the destabilization and eruption of the preexistent flux rope, which was highly-sheared over the polarity inversion line. Driven by the ejection of the flux rope, the overlying loops initially slowly rose and then erupted. As a result, two coronal dimming regions formed at the footprints of the loops and a B1.9 flare, about half an hour later, occurred in the eruptive region near one foot of the eruptive loops. The CME had a very close temporal and spatial relationship to the expanding loops, and it is very likely that the CME directly developed from the eruptive loops. The CME mass calculated from the EUV dimming was of the same order of the typical CME mass, suggesting that the dimming region supplied major of the mass for the CME. The kinematic evolution of the CME was basically consistent with the three-phases scenario: the initiation phase, impulsive acceleration phase, and propagation phase. The onset of the CME could be well explained by the emerging flux trigger mechanism.     

The Observational Evidence on the Loop–Loop Interaction in a Flare–CME Event on April 15, 1998

The evolution of the soft X-ray and EUV coronal loops related to the April 15, 1998 solar flare–CME event is studied with multiwavelength observations including hard X-rays (BATSE), microwaves (NoRP, CNAO) and magnetograms (SOHO/MDI), as well as images from Yohkoh/SXT and SOHO/EIT at 195 Å. It is shown that: (1) two soft X-ray and EUV loops rose, crossed and turned bright, (2) near one footpoint of these loops, the background magnetic field decreased, (3) there were similar quasi periodic oscillations in the time profiles of hard X-ray and microwave emissions, which characterized the loop–loop coalescence instability, (4) after the loop–loop reconnection, two new loops formed, the small one stayed at the original place, and the large one ejected out as part of the constructed prominence cloud. Based upon these observations, we argue that the decrease of the background magnetic field near these loops caused them to rise and approach each other, and in turn, the fast loop–loop coalescence instability took place and triggered the flare and the CME.     

Extreme Ultraviolet Late-Phase Flares: Before and During the Solar Dynamics Observatory Mission

The solar extreme-ultraviolet (EUV) observations from the Solar Dynamics Observatory (SDO) have revealed interesting characteristics of warm coronal emissions, such as Fe xvi 335 Å emission, which peak soon after the hot coronal X-ray emissions peak during a flare and then sometimes peak for a second time hours after the X-ray flare peak. This flare type, with two warm coronal emission peaks but only one X-ray peak, has been named the EUV late phase (Woods et al., Astrophys. J. 739, 59, 2011). These flares have the distinct properties of i) having a complex magnetic-field structure with two initial sets of coronal loops, with one upper set overlaying a lower set, ii) having an eruptive flare initiated in the lower set and disturbing both loop sets, iii) having the hot coronal emissions emitted only from the lower set in conjunction with the X-ray peak, and iv) having the first peak of the warm coronal emissions associated with the lower set and its second peak emitted from the upper set many minutes to hours after the first peak and without a second X-ray enhancement. The disturbance of the coronal loops by the eruption is at about the same time, but the relaxation and cooling down of the heated coronal loops during the post-flare reconnections have different time scales with the longer, upper loops being significantly delayed from the lower loops. The difference in these cooling time scales is related to the difference between the two peak times of the warm coronal emission and is also apparent in the decay profile of the X-ray emissions having two distinct decays, with the first decay slope being steeper (faster) and the delayed decay slope being smaller (slower) during the time of the warm-coronal-emission second peak. The frequency and relationship of the EUV late-phase decay times between the Fe xvi 335 Å two flare peaks and X-ray decay slopes are examined using three years of SDO/EUV Variability Experiment (EVE) data, and the X-ray dual-decay character is then exploited to estimate the frequency of EUV late-phase flares during the past four solar cycles. This study indicates that the frequency of EUV late-phase flares peaks before and after each solar-cycle minimum.     

Damped Oscillations of Coronal Loops

A mechanism of damped oscillations of a coronal loop is investigated. The loop is treated as a thin toroidal flux rope with two stationary photospheric footpoints, carrying both toroidal and poloidal currents. The forces and the flux-rope dynamics are described within the framework of ideal magnetohydrodynamics (MHD). The main features of the theory are the following: i) Oscillatory motions are determined by the Lorentz force that acts on curved current-carrying plasma structures and ii) damping is caused by drag that provides the momentum coupling between the flux rope and the ambient coronal plasma. The oscillation is restricted to the vertical plane of the flux rope. The initial equilibrium flux rope is set into oscillation by a pulse of upflow of the ambient plasma. The theory is applied to two events of oscillating loops observed by the Transition Region and Coronal Explorer (TRACE). It is shown that the Lorentz force and drag with a reasonable value of the coupling coefficient (c _d ) and without anomalous dissipation are able to accurately account for the observed damped oscillations. The analysis shows that the variations in the observed intensity can be explained by the minor radial expansion and contraction. For the two events, the values of the drag coefficient consistent with the observed damping times are in the range c _d ≈2 – 5, with specific values being dependent on parameters such as the loop density, ambient magnetic field, and the loop geometry. This range is consistent with a previous MHD simulation study and with values used to reproduce the observed trajectories of coronal mass ejections (CMEs).     

Automated Detection of Coronal Loops Using a Wavelet Transform Modulus Maxima Method

We propose and test a wavelet transform modulus maxima method for the automated detection and extraction of coronal loops in extreme ultraviolet images of the solar corona. This method decomposes an image into a number of size scales and tracks enhanced power along each ridge corresponding to a coronal loop at each scale. We compare the results across scales and suggest the optimum set of parameters to maximize completeness, while minimizing detection of noise. For a test coronal image, we compare the global statistics (e.g. number of loops at each length) to previous automated coronal-loop detection algorithms.     

The temperature structure and pressure balance of magnetic loops in active regions

EUV observations show many active region loops in lines formed at temperatures between 10⁴K and 2×l0⁶K. The brightest loops are associated with flux tubes leading to the umbrae of sunspots. It is shown that the high visibility of certain loops in transition region lines is due principallly to a sharp radial decrease of temperature to chromospheric values toward the loop axis. The plasma density of these cool loops is not significantly greater than in the hot gas immediately surrounding it. Consequently, the internal gas pressure of the cool material is clearly lower. The hot material immediately surrounding the cool loops is generally denser than the external corona by a factor 3–4. When the active region is examined in coronal lines, this hot high pressure plasma shows up as loops that are generally parallel to the cool loops but significantly displaced laterally. In general the loop phenomenon in an active region is the result of temperature variations by two orders of magnitude and density variations of around a factor five between adjacent flux tubes in the corona.     

A Study of Surges: II. On the Relationship between Chromospheric Surges and Coronal Mass Ejections

Liu et al. (Astrophys. J. 628, 1056, 2005a) described one surge – coronal mass ejection (CME) event showing a close relationship between solar chromospheric surge ejection and CME that had not been noted before. In this work, large Hα surges (>72 Mm, or 100 arcsec) are studied. Eight of these were associated with CMEs. According to their distinct morphological features, Hα surges can be classified into three types: jetlike, diffuse, and closed loop. It was found that all of the jetlike surges were associated with jetlike CMEs (with angular widths ≤30 degrees); the diffuse surges were all associated with wide-angle CMEs (e.g., halo); the closed-loop surges were not associated with CMEs. The exclusive relation between Hα surges and CMEs indicates difference in magnetic field configurations. The jetlike surges and related narrow CMEs propagate along coronal fields that are originally open. The unusual transverse mass motions in the diffuse surges are suggested to be due to magnetic reconnections in the corona that produce wide-angle CMEs. For the closed-loop surges, their paths are just outlining stable closed loops close to the solar surface. Thus no CMEs are associated with them.     

Global Coronal Seismology

Following the observation and analysis of large-scale coronal-wave-like disturbances, we discuss the theoretical progress made in the field of global coronal seismology. Using simple mathematical techniques we determine average values for the magnetic field together with a magnetic map of the quiet Sun. The interaction between global coronal waves and coronal loops allows us to study loop oscillations in a much wider context, i.e. we connect global and local coronal oscillations.     

Forward-Modeling of Doppler Shifts in EUV Spectral Lines

The interpretation of red- and blueshifts in EUV spectral observations remains a challenge that could provide important clues to the heating processes in the solar atmosphere. Hinode/EUV Imaging Spectrometer (EIS) observations near the footpoints of coronal loops show blueshifts for emission lines with temperatures above 1 MK and redshifts for lines below 1 MK. The implications are addressed through numerical modeling of loop dynamics. The simulation results are converted into synthetic EIS observations. A single one-dimensional loop cannot reproduce the observed behavior. However, persistent red- and blueshifts can be understood as a collective spectral signature of a bundle of 10 or more loops that have an average temperature of around 1 MK and evolve in a similar way: small-scale heating events occur randomly along each loop on a timescale of several minutes. Strong blueshifts are accompanied by low intensities. The power-law index of the energy distribution has a minor role in determining the average Doppler shifts.     

Solar Flare and CME Observations with STEREO/EUVI

STEREO/EUVI observed 185 flare events (detected above the GOES class C1 level or at >?25 keV with RHESSI) during the first two years of the mission (December 2006?–?November 2008), while coronal mass ejections (CMEs) were reported in about a third of these events. We compile a comprehensive catalog of these EUVI-observed events, containing the peak fluxes in soft X rays, hard X rays, and EUV, as well as a classification and statistics of prominent EUV features: 79% show impulsive EUV emission (coincident with hard X rays), 73% show delayed EUV emission from postflare loops and arcades, 24% represent occulted flares, 17% exhibit EUV dimming, 5% show loop oscillations or propagating waves, and at least 3% show erupting filaments. We analyze an example of each EUV feature by stereoscopic modeling of its 3D geometry. We find that EUV emission can be dominated by impulsive emission from a heated, highly sheared, noneruptive filament, in addition to the more common impulsive EUV emission from flare ribbons or the delayed postflare EUV emission that results from cooling of the soft-X-ray-emitting flare loops. Occulted flares allow us to determine CME-related coronal dimming uncontaminated from flare-related EUV emission. From modeling the time evolution of EUV dimming we can accurately quantify the initial expansion of CMEs and determine their masses. Further, we find evidence that coronal loop oscillations are excited by the rapid initial expansion of CMEs. These examples demonstrate that stereoscopic EUV data provide powerful new methods to model the 3D aspects in the hydrodynamics of flares and kinematics of CMEs.     

The pre-onset morphology of the 5 September 1973 flare

Simultaneous visible, EUV, and X-ray observations of magnetic structures before and during the onset of the flare of 5 September 1973 are co-registered and interpreted. Ninety minutes before the flare, intense EUV knots fluctuate near the loops which subsequently flare. The pre-flare loop is observed in O IV 554, but not in X-rays, which show instead a parallel structure which is related either to a darkening filament or the subsequent flare kernels. As the full disk X-ray emission increases, first the EUV flare loop appears, then X-ray kernels form at the feet of two EUV loops, one of which overlies the activated filament. The flaring, at any given time, is confined to a single loop (or bundle of loops) whose long axis (barely) crosses the neutral line. As time progresses, the flaring moves to other (probably higher) loops sharing the off-band H footpoints but whose axes are rotated relative to the earlier loops by angles of about 30°. Previous interpretations of single-telescope observations are revised in this joint investigation.     

Propagation of the Slow Magnetoacoustic Waves in Coronal Loops Above Sunspots

The observations from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) have revealed the weak dis- turbances (WDs) propagating in the fan-like coronal loops of the active region (AR 11092) at 171 ?A, 193 ?A, and 211 ?A. These WDs seem to be a common phenomenon in this part of the active region. The disturbances originate from the bright loop foot, and propagate along the loops. The observed propagation speed decreases with the increasing temperature, and varies between 40 km/s and 121 km/s, close to and less than the sound speed in coronal loops. Consid- ering the projection effect and the different angles of the loops with respect to the line of sight, this is exactly what the slow-wave model expects. The wavelet analysis shows that the periods of the WDs observed in different wavebands have no signi?cant difference, the two distinct periods, 3 min and more than 10 min, are all detected in the three EUV wavebands. Not only the coronal loops but also the sunspot region in the chromosphere exhibit intensity oscillations with a period of the order of 3 min. This result suggests that the sunspot oscillations can propagate into the corona through the chromosphere and transition region.     

Dynamical Features and Evolutional Characteristics of Brightening Coronal Loops

We present a detailed study of coronal loop brightenings observed in an active region on the solar limb. These brightening loops show expanding and shrinking motions in EUV coronal line images and also show downflow along the loops in Lα and Hα images. By means of time-slice analysis of the images, we found that both the expanding and shrinking motions of the loops are not real motions of plasma but apparent motions like post-flare loops, where the loops at the different height are successively heated and cooled. From a temperature analysis, the time delay between the brightenings of hot 195 Å and cool Lα loops is found to be nearly equal to the time-scale of the conduction cooling. We conclude that these loop brightenings are sources of so called Hα coronal rains.     

Culgoora radio and skylab euv observations of emerging magnetic flux in the lower corona

Detailed comparisons of Culgoora 160 MHz radioheliograms of solar noise storms and Skylab EUV spectroheliograms of coronal loop structures are presented. It is concluded that: (1) there is a close association between changes in large-scale magnetic fields in the corona and the onset or cessation of noise storms; (2) these coronal changes result from the emergence of new magnetic flux at the photospheric level; (3) although new magnetic flux at the photospheric level is often accompanied by an increase in flare activity the latter is not directly responsible for noise storm activity; rather the new magnetic flux diffuses slowly outwards through the corona at rates 1–2 km s^–1 and produces noise storms at 160 MHz 1–2 days later; (4) the coronal density above or in large-scale EUV loop systems is sufficiently dense to account for noise storm emission at the fundamental plasma frequency; (5) the scatter in noise storm positions can be accounted for by the appearance and disappearance of individual loops in a system.     

On the Coronal Mass Ejection onset and Coronal Dimming

A comprehensive case and statistical study of CME onsets has been conducted on the solar limb using the CDS, LASCO and EIT instruments aboard the SOHO spacecraft. This is the first dedicated campaign to establish firmly the EUV signatures of CME onsets and is based on a series of low-corona observing campaigns made in 2002. The event database consisted of 36 multiple emission line sequences observed with CDS and the study builds, in particular, on studies of EUV coronal dimming which have been associated with CME onsets. We witness a range of dimming events in EUV coronal emission line data. Shorter events, commonly of duration < 4 hours, we find are indirectly associated with CME onsets whereas longer-duration dimmings (> 4 hours) appear to be either due to coronal evolution or rotational effects. However, for some CME onsets, where the CDS pointing was appropriate, no dimming was observed. Dimming observed in EIT typically occurred immediately after the launch of a loop or prominence, and in 5 out of 9 events there is evidence of a matter buildup within the loop before launch. A total of 10 events occurred where CDS was used to directly observe the CME footprint, but no relationship between these events was found. The results suggest that the response of the corona to a CME launch differs between the low (1.0 R _⊙≤R≤1.2 R _⊙) and middle (1.2 R _⊙<R≤2.0 R _⊙) corona regions, hence implying a difference between dimming observations conducted with different instruments.