Fractal and Multifractal Properties of Active Regions as Flare Precursors: A Case Study Based on SOHO/MDI and SDO/HMI Observations

Several studies indicate that fractal and multifractal parameters inferred from solar photospheric magnetic field measurements may help assessing the eruptive potential of Active Regions (ARs) and also predicting their flare activity. We further investigate this topic, by exploring the sensitivity of some parameters already used in the literature on data and methods employed for their estimation. In particular, we measured the generalized fractal dimensions D ₀ and D ₈, and the multifractal parameters C _div and D _div, on the time series of photospheric magnetograms of the flaring AR NOAA 11158 obtained with the SOHO/MDI and SDO/HMI. The observations by the latter instrument are characterized by a higher spatial and temporal resolution, as well as higher flux sensitivity, than the ones obtained from SOHO/MDI, which were widely employed in earlier studies. We found that the average and peak values of complexity parameters measured on the two data sets agree within measurement uncertainties. The temporal evolution of the parameters measured on the two data sets show rather similar trends, but the ones derived from the SOHO/MDI observations show larger and spurious variations over time than those deduced from analysis of the corresponding SDO/HMI data. We also found a larger sensitivity of these measurements to characteristics of the data analyzed than reported by earlier studies. In particular, analysis of the higher resolution and higher cadence SDO/HMI data allows us also to detect slight variations of the complexity indicators that cannot be derived from the analysis of the SOHO/MDI data. These variations occur right after the major events in the analyzed AR. They may be the signature of photospheric effects of coronal magnetic field re-arrangement.     

Long-term oscillations of sunspots and a special class of artifacts in SOHO/MDI and SDO/HMI data

A specific type of artifacts (named as “p2p”), that originate due to displacement of the image of a moving object along the digital (pixel) matrix of receiver are analyzed in detail. The criteria of appearance and the influence of these artifacts on the study of long-term oscillations of sunspots are deduced. The obtained criteria suggest us methods for reduction or even elimination of these artifacts. It is shown that the use of integral parameters can be very effective against the “p2p” artifact distortions. The simultaneous observations of sunspot magnetic field and ultraviolet intensity of the umbra have given the same periods for the long-term oscillations. In this way the real physical nature of the oscillatory process, which is independent of the artifacts have been confirmed again. A number of examples considered here confirm the dependence between the periods of main mode of the sunspot magnetic field long-term oscillations and its strength. The dependence was derived earlier from both the observations and the theoretical model of the shallow sunspot. The anti-phase behavior of time variations of sunspot umbra area and magnetic field of the sunspot demonstrates that the umbra of sunspot moves in long-term oscillations as a whole: all its points oscillate with the same phase.     

Multi-height Scattering Phase Shifts from Sunspots and Phase Differences Between Photospheric Heights with SDO/HMI and AIA Data

Following Couvidat (Solar Phys. 282, 15, 2013), we analyze data from the Helioseismic and Magnetic Imager (HMI) and the Atmospheric Imaging Assembly (AIA) instruments onboard the Solar Dynamics Observatory. Doppler velocity and continuum intensity at 6173 Å as well as intensity maps at 1600 Å and 1700 Å are studied on 14 active regions and four quiet-Sun regions at four heights in the solar photosphere. A Hankel–Fourier analysis is performed around these regions of interest. Outgoing–ingoing scattering phase shifts at a given atmospheric height are computed, as well as ingoing–ingoing and outgoing–outgoing phase differences between two atmospheric heights. The outgoing–ingoing phase shifts produced by sunspots show little dependence on measurement height, unlike the acoustic power reduction measured in Couvidat (2013). Phenomena happening above the continuum level, like acoustic glories, appear not to have a significant effect on the phases. However, there is a strong dependence on sunspot size, line-of-sight magnetic flux, and intensity contrast. As previously suggested by other groups, the region of wave scattering appears both horizontally smaller and vertically less extended than the region of acoustic power suppression, and occurs closer to the surface. Results presented here should help refine theoretical models of acoustic wave scattering by magnetic fields. Ingoing–ingoing phase differences between two measurement heights are also computed and show a complex pattern highly dependent on atmospheric height. In particular, a significant sensitivity of AIA signals to lower chromospheric layers is shown. Finally, ingoing–ingoing minus outgoing–outgoing phase differences between various measurement heights are discussed.     

Resolving Differences in Absolute Irradiance Measurements Between the SOHO/CELIAS/SEM and the SDO/EVE

The Solar EUV Monitor (SEM) onboard SOHO has measured absolute extreme ultraviolet (EUV) and soft X-ray solar irradiance nearly continuously since January 1996. The EUV Variability Experiment (EVE) on SDO, in operation since April of 2010, measures solar irradiance in a wide spectral range that encompasses the band passes (26?–?34 nm and 0.1?–?50 nm) measured by SOHO/SEM. However, throughout the mission overlap, irradiance values from these two instruments have differed by more than the combined stated uncertainties of the measurements. In an effort to identify the sources of these differences and eliminate them, we investigate in this work the effect of reprocessing the SEM data using a more accurate SEM response function (obtained from synchrotron measurements with a SEM sounding-rocket clone instrument taken after SOHO was already in orbit) and time-dependent, measured solar spectral distributions – i.e., solar reference spectra that were unavailable prior to the launch of the SDO. We find that recalculating the SEM data with these improved parameters reduces mean differences with the EVE measurements from about 20 % to less than 5 % in the 26?–?34 nm band, and from about 35 % to about 15 % for irradiances in the 0.1?–?7 nm band extracted from the SEM 0.1?–?50 nm channel.     

Spectro-Polarimetric Observations of Solar Magnetic Fields and the SOHO/MDI Calibration Issue

Comparisons of solar magnetic-field measurements made in different spectral lines are very important, especially in those lines in which observations have a long history or (and) specific diagnostic significance. The spectral lines Fe i 523.3 nm and Fe i 525.0 nm belong to this class. Therefore, this study is devoted to a comprehensive analysis using new high-precision Stokes-meter full-disk observations. The disk-averaged magnetic-field strength ratio R=B(523.3)/B(525.0) equals 1.97±0.02. The center-to-limb variation (CLV) is R=1.74−2.43μ+3.43μ ², where μ is the cosine of the center-to-limb angle. For the disk center, we find R=2.74, and for near-limb areas with μ=0.3, R equals 1.32. There is only a small dependence of R on the spatial resolution. Our results are rather close to those published three decades ago, but differ significantly from recent magnetographic observations. An application of our results to the important SOHO/MDI magnetic data calibration issue is discussed. We conclude that the revision of the SOHO/MDI data, based only on the comparison of magnetic-field measurements in the line pair Fe i 523.3 nm and Fe i 525.0 nm (increasing by a factor of 1.7 or 1.6 on average according to recent publications) is not obvious and new investigations are urgently needed.     

Dynamics of the solar magnetic field from SOHO/MDI

The investigation of the dynamics of magnetic fields from small scales to the large scales is very important for the understanding of the nature of solar activity. It is also the base for producing adequate models of the solar cycle with the purpose to predict the level of solar activity. Since December 1995 the Michelson Doppler Imager (MDI) on board of the Solar and Heliospheric Observatory (SOHO) provides full disk magnetograms and synoptic maps which cover the period of solar cycle 23 and the current minimum. In this paper, I review the following important topics with a focus on the dynamics of the solar magnetic field. The synoptic structure of the solar cycle; the birth of the solar cycle (overlapping cycles 23 and 24); the relationship of the photospheric magnetic activity and the EUV solar corona, polar magnetic fields and dynamo theory (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Comparison Between Nonlinear Force-Free Field and Potential Field Models Using Full-Disk SDO/HMI Magnetogram

Measurements of magnetic fields and electric currents in the pre-eruptive corona are crucial to the study of solar eruptive phenomena, like flares and coronal mass ejections (CMEs). However, spectro-polarimetric measurements of certain photospheric lines permit a determination of the vector magnetic field only at the photosphere. Therefore, there is considerable interest in accurate modeling of the solar coronal magnetic field using photospheric vector magnetograms as boundary data. In this work, we model the coronal magnetic field above multiple active regions with the help of a potential field and a nonlinear force-free field (NLFFF) extrapolation code over the full solar disk using Helioseismic and Magnetic Imager (SDO/HMI) data as boundary conditions. We compare projections of the resulting magnetic field lines with full-disk coronal images from the Atmospheric Imaging Assembly (SDO/AIA) for both models. This study has found that the NLFFF model reconstructs the magnetic configuration closer to observation than the potential field model for full-disk magnetic field extrapolation. We conclude that many of the trans-equatorial loops connecting the two solar hemispheres are current-free.     

The Energetics of White-light Flares Observed by SDO/HMI and RHESSI

White-light(WL) flares have been observed and studied for more than a century since their first discovery. However, some fundamental physics behind the brilliant emission remains highly controversial.One of the important facts in addressing the flare energetics is the spatio-temporal correlation between the WL emission and the hard X-ray(HXR) radiation, presumably suggesting that energetic electrons are the energy sources. In this study, we present a statistical analysis of 25 strong flares(≥M5) observed simultaneously by the Helioseismic and Magnetic Imager(HMI), on board the Solar Dynamics Observatory(SDO),and the Reuven Ramaty High Energy Solar Spectroscopic Imager(RHESSI). Among these events, WL emission was detected by SDO/HMI in 13 flares, associated with HXR emission. To quantitatively describe the strength of WL emission, equivalent area(EA) is defined as the integrated contrast enhancement over the entire flaring area. Our results show that the EA is inversely proportional to the HXR power-law index,indicating that stronger WL emission tends to be associated with a larger population of high energy electrons. However, no obvious correlation is found between WL emission and flux of non-thermal electrons at50 ke V. For the other group of 13 flares without detectable WL emission, the HXR spectra are softer(larger power-law index) than those flares with WL emission, especially for the X-class flares in this group.     

Analysis of Supergranule Sizes and Velocities Using Solar Dynamics Observatory (SDO)/Helioseismic Magnetic Imager (HMI) and Solar and Heliospheric Observatory (SOHO)/Michelson Doppler Imager (MDI) Dopplergrams

Co-temporal Doppler images from Solar and Heliospheric Observatory (SOHO)/Michelson Doppler Imager (MDI) and Solar Dynamics Observatory (SDO)/Helioseismic Magnetic Imager (HMI) have been analyzed to extract quantitative information about global properties of the spatial and temporal characteristics of solar supergranulation. Preliminary comparisons show that supergranules appear to be smaller and have stronger horizontal velocity flows within HMI data than was measured with MDI. There appears to be no difference in their evolutionary timescales. Supergranule sizes and velocities were analyzed over a ten-day time period at a 15-minute cadence. While the averages of the time-series retain the aforementioned differences, fluctuations of these parameters first observed in MDI data were seen in both MDI and HMI time-series, exhibiting a strong cross-correlation. This verifies that these fluctuations are not instrumental, but are solar in origin. The observed discrepancies between the averaged values from the two sets of data are a consequence of instrument resolution. The lower spatial resolution of MDI results in larger observed structures with lower velocities than is seen in HMI. While these results offer a further constraint on the physical nature of supergranules, they also provide a level of calibration between the two instruments.     

An Early Diagnostics of the Geoeffectiveness of Solar Eruptions from Photospheric Magnetic Flux Observations: The Transition from SOHO to SDO

In our previous articles (Chertok et al. in Solar Phys. 282, 175, 2013; Chertok et al. in Solar Phys. 290, 627, 2015), we presented a preliminary tool for the early diagnostics of the geoeffectiveness of solar eruptions based on the estimate of the total unsigned line-of-sight photospheric magnetic flux in accompanying extreme ultraviolet (EUV) arcades and dimmings. This tool was based on the analysis of eruptions observed during 1996?–?2005 with the Extreme-ultraviolet Imaging Telescope (EIT) and the Michelson Doppler Imager (MDI) onboard the Solar and Heliospheric Observatory (SOHO). Empirical relationships were obtained to estimate the probable importance of upcoming space weather disturbances caused by an eruption, which just occurred, without data on the associated coronal mass ejections. In particular, it was possible to estimate the intensity of a non-recurrent geomagnetic storm (GMS) and Forbush decrease (FD), as well as their onset and peak times. After 2010?–?2011, data on solar eruptions are obtained with the Atmospheric Imaging Assembly (AIA) and the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We use relatively short intervals of overlapping EIT–AIA and MDI–HMI detailed observations, and additionally, a number of large eruptions over the next five years with the 12-hour cadence EIT images to adapt the SOHO diagnostic tool to SDO data. We show that the adopted brightness thresholds select practically the same areas of arcades and dimmings from the EIT 195 Å and AIA 193 Å image, with a cross-calibration factor of 3.6?–?5.8 (5.0?–?8.2) for the AIA exposure time of 2.0 s (2.9 s). We also find that for the same photospheric areas, the MDI line-of-sight magnetic flux systematically exceeds the HMI flux by a factor of 1.4. Based on these results, the empirical diagnostic relationships obtained from SOHO data are adjusted to SDO instruments. Examples of a post-diagnostics based on SDO data are presented. As before, the tool is applicable to non-recurrent GMSs and FDs caused by nearly central eruptions from active regions, provided that the southern component of the interplanetary magnetic field near the Earth is predominantly negative, which is not predicted by this tool.     

Determination of magnetic helicity content of solar active regions from SOHO/MDI magnetograms

Chae (2001) first proposed a method of self-consistently determining the rate of change of magnetic helicity using a time series of longitudinal magnetograms only, such as taken by SOHO/MDI. Assuming that magnetic fields in the photosphere are predominantly vertical, he determined the horizontal component of velocity by tracking the displacements of magnetic flux fragments using the technique of local correlation tracking (LCT). In the present paper, after briefly reviewing the recent advance in helicity rate measurement, we argue that the LCT method can be more generally applied even to regions of inclined magnetic fields. We also present some results obtained by applying the LCT method to the active region NOAA 10365 under emergence during the observable period, which are summarized as follows. (1) Strong shearing flows were found near the polarity inversion line that were very effective in helicity injection. (2) Both the magnetic flux and helicity of the active region steadily increased during the observing period, and reached 1.2 × 10²² Mx and 8 ×10⁴² Mx², respectively, 4.5 days after the birth of the active region. (3) The corresponding ratio of the helicity to the square of the magnetic flux, 0.05, is roughly compatible with the values determined by other studies using linear-force-free modeling. (4) A series of flares took place while the rate of helicity injection was high. (5) The choice of a smaller window size or a shorter time interval in the LCT method resulted in a bigger value of the LCT velocity and a bigger value of the temporal fluctuation of the helicity rate. (6) Nevertheless when averaged over a time period of about one hour or longer, the average rate of helicity became about the same within about 10%, almost irrespective of the chosen window size and time interval, indicating that short-lived, fluctuating flows may be insignificant in transferring magnetic helicity. Our results suggest that the LCT method may be applied to 96-minute cadence full-disk MDI magnetograms or other data of similar kind, to provide a practically useful, if not perfect, way of monitoring the magnetic helicity content of active regions as a function of time.     

Determination of magnetic helicity content of solar active regions from SOHO/MDI magnetograms

Chae (2001) first proposed a method of self-consistently determining the rate of change of magnetic helicity using a time series of longitudinal magnetograms only, such as taken by SOHO/MDI. Assuming that magnetic fields in the photosphere are predominantly vertical, he determined the horizontal component of velocity by tracking the displacements of magnetic flux fragments using the technique of local correlation tracking (LCT). In the present paper, after briefly reviewing the recent advance in helicity rate measurement, we argue that the LCT method can be more generally applied even to regions of inclined magnetic fields. We also present some results obtained by applying the LCT method to the active region NOAA 10365 under emergence during the observable period, which are summarized as follows. (1) Strong shearing flows were found near the polarity inversion line that were very effective in helicity injection. (2) Both the magnetic flux and helicity of the active region steadily increased during the observing period, and reached 1.2 × 10²² Mx and 8 ×10⁴² Mx², respectively, 4.5 days after the birth of the active region. (3) The corresponding ratio of the helicity to the square of the magnetic flux, 0.05, is roughly compatible with the values determined by other studies using linear-force-free modeling. (4) A series of flares took place while the rate of helicity injection was high. (5) The choice of a smaller window size or a shorter time interval in the LCT method resulted in a bigger value of the LCT velocity and a bigger value of the temporal fluctuation of the helicity rate. (6) Nevertheless when averaged over a time period of about one hour or longer, the average rate of helicity became about the same within about 10%, almost irrespective of the chosen window size and time interval, indicating that short-lived, fluctuating flows may be insignificant in transferring magnetic helicity. Our results suggest that the LCT method may be applied to 96-minute cadence full-disk MDI magnetograms or other data of similar kind, to provide a practically useful, if not perfect, way of monitoring the magnetic helicity content of active regions as a function of time.     

Properties of High-Frequency Wave Power Halos Around Active Regions: An Analysis of Multi-height Data from HMI and AIA Onboard SDO

We study properties of waves of frequencies above the photospheric acoustic cut-off of ≈5.3 mHz, around four active regions, through spatial maps of their power estimated using data from the Helioseismic and Magnetic Imager (HMI) and Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). The wavelength channels 1600 Å and 1700 Å from AIA are now known to capture clear oscillation signals due to helioseismic p-modes as well as waves propagating up through to the chromosphere. Here we study in detail, in comparison with HMI Doppler data, properties of the power maps, especially the so-called “acoustic halos” seen around active regions, as a function of wave frequencies, inclination, and strength of magnetic field (derived from the vector-field observations by HMI), and observation height. We infer possible signatures of (magneto)acoustic wave refraction from the observation-height-dependent changes, and hence due to changing magnetic strength and geometry, in the dependences of power maps on the photospheric magnetic quantities. We discuss the implications for theories of p-mode absorption and mode conversions by the magnetic field.     

Automatic Detection and Extraction of Coronal Dimmings from SDO/AIA Data

The volume of data anticipated from the Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA) highlights the necessity for the development of automatic-detection methods for various types of solar activity. Initially recognized in the 1970s, it is now well established that coronal dimmings are closely associated with coronal mass ejections (CMEs), and they are particularly noted as a reliable indicator of front-side (halo) CMEs, which can be difficult to detect in white-light coronagraph data. Existing work clearly demonstrates that several properties derived from the analysis of coronal dimmings can give useful information about the associated CME. The development and implementation of an automated coronal-dimming region detection and extraction algorithm removes visual observer bias, however unintentional, from the determination of physical quantities such as spatial location, area, and volume. This allows for reproducible, quantifiable results to be mined from very large data sets. The information derived may facilitate more reliable early space-weather detection, as well as offering the potential for conducting large-sample studies focused on determining the geo-effectiveness of CMEs, coupled with analysis of their associated coronal dimming signatures. In this paper we present examples of both simple and complex dimming events extracted using our algorithm, which will be run as a module for the SDO/Computer Vision Centre. Contrasting and well-studied events at both the minimum and maximum of solar cycle 23 are identified in Solar and Heliospheric Observatory/Extreme ultra-violet Imaging Telescope (SOHO/EIT) data. A more recent example extracted from Solar and Terrestrial Relations Observatory/Extreme Ultra-Violet Imager (STEREO/EUVI) data is also presented, demonstrating the potential for the anticipated application to SDO/AIA data. The detection part of our algorithm is based largely on the principle of operation of the NEMO software, namely the detection of significant variation in the statistics of the EUV image pixels (Podladchikova and Berghmans in Solar Phys. 228, 265?–?284, 2005). As well as running on historic data sets, the presented algorithm is capable of detecting and extracting coronal dimmings in near real-time.     

On the Nature and Genesis of EUV Waves: A Synthesis of Observations from SOHO, STEREO, SDO, and Hinode (Invited Review)

A major, albeit serendipitous, discovery of the SOlar and Heliospheric Observatory mission was the observation by the Extreme Ultraviolet Telescope (EIT) of large-scale extreme ultraviolet (EUV) intensity fronts propagating over a significant fraction of the Sun??s surface. These so-called EIT or EUV waves are associated with eruptive phenomena and have been studied intensely. However, their wave nature has been challenged by non-wave (or pseudo-wave) interpretations and the subject remains under debate. A string of recent solar missions has provided a wealth of detailed EUV observations of these waves bringing us closer to resolving the question of their nature. With this review, we gather the current state-of-the-art knowledge in the field and synthesize it into a picture of an EUV wave driven by the lateral expansion of the CME. This picture can account for both wave and pseudo-wave interpretations of the observations, thus resolving the controversy over the nature of EUV waves to a large degree but not completely. We close with a discussion on several remaining open questions in the field of EUV waves research.     

Stable Sunspot Area Level of Debrecen Photoheliographic Data and Multivariate Correction Factor of SOON Data

We investigate the spatial and temporal variation of sunspot group areas reported by the Greenwich Photoheliographic Results (GPR), the Solar Optical Observing Network (SOON), the Kislovodsk Mountain Astronomical Station (KMAS), and the Debrecen Photoheliographic Data (DPD) databases. We identify improved correction factors for reconciling these individual records to a common scale. Our results show that the DPD sunspot group areas are stable over the studied interval (1974?–?2014). We find an improved fit between GPR and DPD sunspot group areas when using a correction factor such that \(\mathrm{GPR} = 0.975(\pm 0.006) \times \mathrm{DPD}\), independent of the position of the sunspot group on the solar disk. We also find that the scale of KMAS sunspot group areas fits that of DPD well, but has a small position-dependent trend near the limb. However, in order to set SOON sunspot group area records onto the scale of DPD, we find that there is a need for a multivariate correction factor. This multivariate correction factor has a value ranging between 1.1 and 1.9 and is dependent upon the time of the SOON observation, the distance of the group from disk center, and the observatory within the SOON network. Finally, we provide further context to the systematic bias in SOON sunspot group area observations toward lower values relative to those recorded in the GPR and DPD databases that has previously been reported in the literature. We have identified the two main contributors to the SOON area deficit; some penumbral parts are unobserved, and the spot areas are underestimated. Our analysis is vital for studies that require stable, long-term solar activity records such as solar irradiance models that estimate irradiance reduction from records of sunspot group numbers, areas, and locations.     

Sounding Rocket Observations of Active Region Soft X-Ray Spectra Between 0.5 and 2.5 nm Using a Modified SDO/EVE Instrument

Spectrally resolved measurements of individual solar active regions (ARs) in the soft X-ray (SXR) range are important for studying dynamic processes in the solar corona and their associated effects on the Earth’s upper atmosphere. They are also a means of evaluating atomic data and elemental abundances used in physics-based solar spectral models. However, very few such measurements are available. We present spectral measurements of two individual ARs in the 0.5 to 2.5 nm range obtained on the NASA 36.290 sounding rocket flight of 21 October 2013 (at about 18:30 UT) using the Solar Aspect Monitor (SAM), a channel of the Extreme Ultaviolet Variability Experiment (EVE) payload designed for underflight calibrations of the orbital EVE on the Solar Dynamics Observatory (SDO). The EVE rocket instrument is a duplicate of the EVE on SDO, except the SAM channel on the rocket version was modified in 2012 to include a freestanding transmission grating to provide spectrally resolved images of the solar disk with the best signal to noise ratio for the brightest features, such as ARs. Calibrations of the EVE sounding rocket instrument at the National Institute of Standards and Technology Synchrotron Ultraviolet Radiation Facility (NIST/SURF) have provided a measurement of the SAM absolute spectral response function and a mapping of wavelength separation in the grating diffraction pattern. We discuss techniques (incorporating the NIST/SURF data) for determining SXR spectra from the dispersed AR images as well as the resulting spectra for NOAA ARs 11877 and 11875 observed on the 2013 rocket flight. In comparisons with physics-based spectral models using the CHIANTI v8 atomic database we find that both AR spectra are in good agreement with isothermal spectra (4 MK), as well as spectra based on an AR differential emission measure (DEM) included with the CHIANTI distribution, with the exception of the relative intensities of strong Fe?xvii lines associated with \(2p^{6}\)–\(2p^{5}3{s}\) and \(2p^{6}\)–\(2p^{5}3{d}\) transitions at about 1.7 nm and 1.5 nm, respectively. The ratio of the Fe?xvii lines suggests that the AR 11877 is hotter than the AR 11875. This result is confirmed with analysis of the active regions imaged by X-ray Telescope (XRT) onboard Hinode.