A Multi-wavelength Analysis of Active Regions and Sunspots by Comparison of Automatic Detection Algorithms

Since the Solar Dynamics Observatory (SDO) began recording ≈?1 TB of data per day, there has been an increased need to automatically extract features and events for further analysis. Here we compare the overall detection performance, correlations between extracted properties, and usability for feature tracking of four solar feature-detection algorithms: the Solar Monitor Active Region Tracker (SMART) detects active regions in line-of-sight magnetograms; the Automated Solar Activity Prediction code (ASAP) detects sunspots and pores in white-light continuum images; the Sunspot Tracking And Recognition Algorithm (STARA) detects sunspots in white-light continuum images; the Spatial Possibilistic Clustering Algorithm (SPoCA) automatically segments solar EUV images into active regions (AR), coronal holes (CH), and quiet Sun (QS). One month of data from the Solar and Heliospheric Observatory (SOHO)/Michelson Doppler Imager (MDI) and SOHO/Extreme Ultraviolet Imaging Telescope (EIT) instruments during 12 May?–?23 June 2003 is analysed. The overall detection performance of each algorithm is benchmarked against National Oceanic and Atmospheric Administration (NOAA) and Solar Influences Data Analysis Center (SIDC) catalogues using various feature properties such as total sunspot area, which shows good agreement, and the number of features detected, which shows poor agreement. Principal Component Analysis indicates a clear distinction between photospheric properties, which are highly correlated to the first component and account for 52.86% of variability in the data set, and coronal properties, which are moderately correlated to both the first and second principal components. Finally, case studies of NOAA 10377 and 10365 are conducted to determine algorithm stability for tracking the evolution of individual features. We find that magnetic flux and total sunspot area are the best indicators of active-region emergence. Additionally, for NOAA 10365, it is shown that the onset of flaring occurs during both periods of magnetic-flux emergence and complexity development.     

Evaluation of a Revised Interplanetary Shock Prediction Model: 1D CESE-HD-2 Solar-Wind Model

We modified the one-dimensional conservation element and solution element (CESE) hydrodynamic (HD) model into a new version [1D CESE-HD-2], by considering the direction of the shock propagation. The real-time performance of the 1D CESE-HD-2 model during Solar Cycle 23 (February 1997?–?December 2006) is investigated and compared with those of the Shock Time of Arrival Model (STOA), the Interplanetary-Shock-Propagation Model (ISPM), and the Hakamada–Akasofu–Fry version 2 (HAFv.2). Of the total of 584 flare events, 173 occurred during the rising phase, 166 events during the maximum phase, and 245 events during the declining phase. The statistical results show that the success rates of the predictions by the 1D CESE-HD-2 model for the rising, maximum, declining, and composite periods are 64 %, 62 %, 57 %, and 61 %, respectively, with a hit window of ±?24 hours. The results demonstrate that the 1D CESE-HD-2 model shows the highest success rates when the background solar-wind speed is relatively fast. Thus, when the background solar-wind speed at the time of shock initiation is enhanced, the forecasts will provide potential values to the customers. A high value (27.08) of χ ² and low p-value (<?0.0001) for the 1D CESE-HD-2 model give considerable confidence for real-time forecasts by using this new model. Furthermore, the effects of various shock characteristics (initial speed, shock duration, background solar wind, longitude, etc.) and background solar wind on the forecast are also investigated statistically.     

The Helioseismic and Magnetic Imager (HMI) Vector Magnetic Field Pipeline: SHARPs – Space-Weather HMI Active Region Patches

A new data product from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) called Space-weather HMI Active Region Patches (SHARPs) is now available. SDO/HMI is the first space-based instrument to map the full-disk photospheric vector magnetic field with high cadence and continuity. The SHARP data series provide maps in patches that encompass automatically tracked magnetic concentrations for their entire lifetime; map quantities include the photospheric vector magnetic field and its uncertainty, along with Doppler velocity, continuum intensity, and line-of-sight magnetic field. Furthermore, keywords in the SHARP data series provide several parameters that concisely characterize the magnetic-field distribution and its deviation from a potential-field configuration. These indices may be useful for active-region event forecasting and for identifying regions of interest. The indices are calculated per patch and are available on a twelve-minute cadence. Quick-look data are available within approximately three hours of observation; definitive science products are produced approximately five weeks later. SHARP data are available at jsoc.stanford.edu and maps are available in either of two different coordinate systems. This article describes the SHARP data products and presents examples of SHARP data and parameters.     

Dependence of DOLP on Coronal Electron Temperature,Speed, and Structure

We study the predictive capabilities of magnetic-feature properties (MF) generated by the Solar Monitor Active Region Tracker (SMART: Higgins et al. in Adv. Space Res. 47, 2105, 2011) for solar-flare forecasting from two datasets: the full dataset of SMART detections from 1996 to 2010 which has been previously studied by Ahmed et al. (Solar Phys. 283, 157, 2013) and a subset of that dataset that only includes detections that are NOAA active regions (ARs). The main contributions of this work are: we use marginal relevance as a filter feature selection method to identify the most useful SMART MF properties for separating flaring from non-flaring detections and logistic regression to derive classification rules to predict future observations. For comparison, we employ a Random Forest, Support Vector Machine, and a set of Deep Neural Network models, as well as lasso for feature selection. Using the linear model with three features we obtain significantly better results (True Skill Score: TSS = 0.84) than those reported by Ahmed et al. (Solar Phys. 283, 157, 2013) for the full dataset of SMART detections. The same model produced competitive results (TSS = 0.67) for the dataset of SMART detections that are NOAA ARs, which can be compared to a broader section of flare-forecasting literature. We show that more complex models are not required for this data.

     

Transient Artifacts in a Flare Observed by the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory

The Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) provides a new tool for the systematic observation of white-light flares, including Doppler and magnetic information as well as continuum. In our initial analysis of the highly impulsive $\mathrm{\gamma}$ -ray flare SOL2010-06-12T00:57 (Martínez Oliveros et al., Solar Phys. 269, 269, 2011), we reported the signature of a strong blueshift in the two footpoint sources. Concerned that this might be an artifact due to aliasing peculiar to the HMI instrument, we undertook a comparative analysis of Global Oscillation Network Group (GONG++) observations of the same flare, using the PArametric Smearing Correction ALgorithm (PASCAL) algorithm to correct for artifacts caused by variations in atmospheric smearing. This analysis confirms the artifactual nature of the apparent blueshift in the HMI observations, finding weak redshifts at the footpoints instead. We describe the use of PASCAL with GONG++ observations as a complement to the SDO observations and discuss constraints imposed by the use of HMI far from its design conditions. With proper precautions, these data provide rich information on flares and transients.     

EIDOSCOPE: particle acceleration at plasma boundaries

We describe the mission concept of how ESA can make a major contribution to the Japanese Canadian multi-spacecraft mission SCOPE by adding one cost-effective spacecraft EIDO (Electron and Ion Dynamics Observatory), which has a comprehensive and optimized plasma payload to address the physics of particle acceleration. The combined mission EIDOSCOPE will distinguish amongst and quantify the governing processes of particle acceleration at several important plasma boundaries and their associated boundary layers: collisionless shocks, plasma jet fronts, thin current sheets and turbulent boundary layers. Particle acceleration and associated cross-scale coupling is one of the key outstanding topics to be addressed in the Plasma Universe. The very important science questions that only the combined EIDOSCOPE mission will be able to tackle are: 1) Quantitatively, what are the processes and efficiencies with which both electrons and ions are selectively injected and subsequently accelerated by collisionless shocks? 2) How does small-scale electron and ion acceleration at jet fronts due to kinetic processes couple simultaneously to large scale acceleration due to fluid (MHD) mechanisms? 3) How does multi-scale coupling govern acceleration mechanisms at electron, ion and fluid scales in thin current sheets? 4) How do particle acceleration processes inside turbulent boundary layers depend on turbulence properties at ion/electron scales? EIDO particle instruments are capable of resolving full 3D particle distribution functions in both thermal and suprathermal regimes and at high enough temporal resolution to resolve the relevant scales even in very dynamic plasma processes. The EIDO spin axis is designed to be sun-pointing, allowing EIDO to carry out the most sensitive electric field measurements ever accomplished in the outer magnetosphere. Combined with a nearby SCOPE Far Daughter satellite, EIDO will form a second pair (in addition to SCOPE Mother-Near Daughter) of closely separated satellites that provides the unique capability to measure the 3D electric field with high accuracy and sensitivity. All EIDO instrumentation are state-of-the-art technology with heritage from many recent missions. The EIDOSCOPE orbit will be close to equatorial with apogee 25-30 RE and perigee 8-10 RE. In the course of one year the orbit will cross all the major plasma boundaries in the outer magnetosphere; bow shock, magnetopause and magnetotail current sheets, jet fronts and turbulent boundary layers. EIDO offers excellent cost/benefits for ESA, as for only a fraction of an M-class mission cost ESA can become an integral part of a major multi-agency L-class level mission that addresses outstanding science questions for the benefit of the European science community.     

The Astro-WISE approach to quality control for astronomical data

We present a novel approach to quality control during the processing of astronomical data. Quality control in the Astro-WISE Information System is integral to all aspects of data handing and provides transparent access to quality estimators for all stages of data reduction from the raw image to the final catalog. The implementation of quality control mechanisms relies on the core features in this Astro-WISE Environment (AWE): an object-oriented framework, full data lineage, and both forward and backward chaining. Quality control information can be accessed via the command-line awe-prompt and the web-based Quality-WISE service. The quality control system is described and qualified using archive data from the 8-CCD Wide Field Imager (WFI) instrument (http://www.eso.org/lasilla/instruments/wfi/) on the 2.2-m MPG/ESO telescope at La Silla and (pre-)survey data from the 32-CCD OmegaCAM instrument (http://www.astro-wise.org/~omegacam/) on the VST telescope at Paranal.     

Observation of the impulsive phase of a simple flare

We present a broad range of complementary observations of the onset and impulsive phase of a fairly large (1B, M1.2) but simple two-ribbon flare. The observations consist of hard X-ray flux measured by the SMM HXRBS, high-sensitivity measurements of microwave flux at 22 GHz from Itapetinga Radio Observatory, sequences of spectroheliograms in UV emission lines from Ov (T ≈ 2 × 10⁵ K) and Fexxi (T ≈ 1 × 10⁷ K) from the SMM UVSP, Hα and Hei D₃ cine-filtergrams from Big Bear Solar Observatory, and a magnetogram of the flare region from the MSFC Solar Observatory. From these data we conclude:

1. The overall magnetic field configuration in which the flare occurred was a fairly simple, closed arch containing nonpotential substructure.
2. The flare occurred spontaneously within the arch; it was not triggered by emerging magnetic flux.
3. The impulsive energy release occurred in two major spikes. The second spike took place within the flare arch heated in the first spike, but was concentrated on a different subset of field lines. The ratio of Ov emission to hard X-ray emission decreased by at least a factor of 2 from the first spike to the second, probably because the plasma density in the flare arch had increased by chromospheric evaporation.
4. The impulsive energy release most likely occurred in the upper part of the arch; it had three immediate products:

1. An increase in the plasma pressure throughout the flare arch of at least a factor of 10. This is required because the Fexxi emission was confined to the feet of the flare arch for at least the first minute of the impulsive phase.
2. Nonthermal energetic (～ 25 keV) electrons which impacted the feet of the arch to produce the hard X-ray burst and impulsive brightening in Ov and D₃. The evidence for this is the simultaneity, within ± 2 s, of the peak Ov and hard X-ray emissions.
3. Another population of high-energy (～100keV) electrons (decoupled from the population that produced the hard X-rays) that produced the impulsive microwave emission at 22 GHz. This conclusion is drawn because the microwave peak was 6 ± 3 s later than the hard X-ray peak.

     

CorPITA: An Automated Algorithm for the Identification and Analysis of Coronal “EIT Waves”

The continuous stream of data available from the Atmospheric Imaging Assembly (AIA) telescopes onboard the Solar Dynamics Observatory (SDO) spacecraft has allowed a deeper understanding of the Sun. However, the sheer volume of data has necessitated the development of automated techniques to identify and analyse various phenomena. In this article, we describe the Coronal Pulse Identification and Tracking Algorithm (CorPITA) for the identification and analysis of coronal “EIT waves”. CorPITA uses an intensity-profile technique to identify the propagating pulse, tracking it throughout its evolution before returning estimates of its kinematics. The algorithm is applied here to a data set from February 2011, allowing its capabilities to be examined and critiqued. This algorithm forms part of the SDO Feature Finding Team initiative and will be implemented as part of the Heliophysics Event Knowledgebase (HEK). This is the first fully automated algorithm to identify and track the propagating “EIT wave” rather than any associated phenomenon and will allow a deeper understanding of this controversial phenomenon.     

Interpreting the Helioseismic and Magnetic Imager (HMI) Multi-Height Velocity Measurements

The Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) filtergrams, taken at six wavelengths around the Fe i 6173.3 Å line, contain information about the line-of-sight velocity over a range of heights in the solar atmosphere. Multi-height velocity inferences from these observations can be exploited to study wave motions and energy transport in the atmosphere. Using realistic convection-simulation datasets provided by the STAGGER and MURaM codes, we generate synthetic filtergrams and explore several methods for estimating Dopplergrams. We investigate at which height each synthetic Dopplergram correlates most strongly with the vertical velocity in the model atmospheres. On the basis of the investigation, we propose two Dopplergrams other than the standard HMI-algorithm Dopplergram produced from HMI filtergrams: a line-center Dopplergram and an average-wing Dopplergram. These two Dopplergrams correlate most strongly with vertical velocities at the heights of 30?–?40 km above (line center) and 30?–?40 km below (average wing) the effective height of the HMI-algorithm Dopplergram. Therefore, we can obtain velocity information from two layers separated by about a half of a scale height in the atmosphere, at best. The phase shifts between these multi-height Dopplergrams from observational data as well as those from the simulated data are also consistent with the height-difference estimates in the frequency range above the photospheric acoustic-cutoff frequency.     

Automated Solar Flare Detection and Feature Extraction in High-Resolution and Full-Disk H$\upalpha$ Images

In this article, an automated solar flare detection method applied to both full-disk and local high-resolution H\(\upalpha\) images is proposed. An adaptive gray threshold and an area threshold are used to segment the flare region. Features of each detected flare event are extracted, e.g. the start, peak, and end time, the importance class, and the brightness class. Experimental results have verified that the proposed method can obtain more stable and accurate segmentation results than previous works on full-disk images from Big Bear Solar Observatory (BBSO) and Kanzelhöhe Observatory for Solar and Environmental Research (KSO), and satisfying segmentation results on high-resolution images from the Goode Solar Telescope (GST). Moreover, the extracted flare features correlate well with the data given by KSO. The method may be able to implement a more complicated statistical analysis of H\(\upalpha\) solar flares.     

The Astro-WISE optical image pipeline

We have designed and implemented a novel way to process wide-field astronomical data within a distributed environment of hardware resources and humanpower. The system is characterized by integration of archiving, calibration, and post-calibration analysis of data from raw, through intermediate, to final data products. It is a true integration thanks to complete linking of data lineage from the final catalogs back to the raw data. This paper describes the pipeline processing of optical wide-field astronomical data from the WFI (http://www.eso.org/lasilla/instruments/wfi/) and OmegaCAM (http://www.astro-wise.org/~omegacam/) instruments using the Astro-WISE information system (the Astro-WISE Environment or simply AWE). This information system is an environment of hardware resources and humanpower distributed over Europe. AWE is characterized by integration of archiving, data calibration, post-calibration analysis, and archiving of raw, intermediate, and final data products. The true integration enables a complete data processing cycle from the raw data up to the publication of science-ready catalogs. The advantages of this system for very large datasets are in the areas of: survey operations management, quality control, calibration analyses, and massive processing.     

On the Formation Height of the SDO/HMI Fe 6173 ? Doppler Signal

The Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) is designed to study oscillations and the magnetic field in the solar photosphere. It observes the full solar disk in the Fe?i absorption line at 6173 Å. We use the output of a high-resolution, 3D, time-dependent, radiation-hydrodynamic simulation based on the CO ⁵ BOLD code to calculate profiles F(??,x,y,t) for the Fe?i 6173 Å line. The emerging profiles F(??,x,y,t) are multiplied by a representative set of HMI filter-transmission profiles R _i (??, 1??i??6) and filtergrams I _i (x,y,t; 1??i??6) are constructed for six wavelengths. Doppler velocities V _HMI(x,y,t) are determined from these filtergrams using a simplified version of the HMI pipeline. The Doppler velocities are correlated with the original velocities in the simulated atmosphere. The cross-correlation peaks near 100 km, suggesting that the HMI Doppler velocity signal is formed rather low in the solar atmosphere. The same analysis is performed for the SOHO/MDI Ni?i line at 6768 Å. The MDI Doppler signal is formed slightly higher at around 125 km. Taking into account the limited spatial resolution of the instruments, the apparent formation height of both the HMI and MDI Doppler signal increases by 40 to 50 km. We also study how uncertainties in the HMI filter-transmission profiles affect the calculated velocities.     

TomograPy: A Fast, Instrument-Independent, Solar Tomography Software

Solar tomography has progressed rapidly in recent years thanks to the development of robust algorithms and the availability of more powerful computers. It can today provide crucial insights in solving issues related to the line-of-sight integration present in the data of solar imagers and coronagraphs. However, there remain challenges such as the increase of the available volume of data, the handling of the temporal evolution of the observed structures, and the heterogeneity of the data in multi-spacecraft studies. We present a generic software package that can perform fast tomographic inversions that scales linearly with the number of measurements, linearly with the length of the reconstruction cube (and not the number of voxels), and linearly with the number of cores and can use data from different sources and with a variety of physical models: TomograPy ( http://nbarbey.github.com/TomograPy/ ), an open-source software freely available on the Python Package Index. For performance, TomograPy uses a parallelized-projection algorithm. It relies on the World Coordinate System standard to manage various data sources. A variety of inversion algorithms are provided to perform the tomographic-map estimation. A test suite is provided along with the code to ensure software quality. Since it makes use of the Siddon algorithm it is restricted to rectangular parallelepiped voxels but the spherical geometry of the corona can be handled through proper use of priors. We describe the main features of the code and show three practical examples of multi-spacecraft tomographic inversions using STEREO/EUVI and STEREO/COR1 data. Static and smoothly varying temporal evolution models are presented.     

Lifetimes of some excited states in the rare earths: La ii,Ce ii,Pr ii,Nd ii,Sm ii,Yb i,Yb ii,and Lu ii

Radiative lifetimes for excited states in La ii, Ce ii, Pr ii, Nd ii, Sm ii, Yb i, Yb ii, and Lu ii have been determined by means of the beam-foil technique or the zero-field level-crossing method. The lifetimes for La, Ce, Pr, Nd, Sm, and Yb are shorter than those computed by summing the transition probabilities of Corliss and Bozman by a factor of up to ~5. The large discrepancies between the abundance of La, Ce, Pr, Nd, and Sm in the solar photosphere and in meteorites are eliminated or greatly reduced if the abundance determinations of the solar photosphere are based upon the gf values of Corliss and Bozman, corrected for by the present lifetimes.     

An evaluation of the possibility of studying flare plasma turbulence using the satellites of Hei line forbidden components

Using the Baranger-Mozer method, we explore the possibility of diagnosing the flare plasma of forbidden Hei lines, that permits the determination of the plasma oscillation frequency and noise level. Examination of the Hei lines observed in solar flare has led us to conclude that:

1. the appearance of satellites of forbidden components in the flares spectrum, due to turbulent electric fields, is the most probable for Hei 3819.606 Å lines;
2. the Baranger-Mozer method is more sensitive to the high-frequency component of turbulent fields than to the low-frequency ones;
3. the upper limit of the turbulent oscillation level in flares is evaluated.

In the spectrum of the solar flare of 26 September, 1963 we detected satellites of the forbidden component of the 3820 Å line and used its relative intensity to derive the level of low-frequency oscillations (～1.5 kVcm^-1).     

Variability of Fe ii in Two NLS1s, PG 1700+815 and NGC 4051

We analyze the spectral variability for two narrow line Seyfert 1 galaxies, PG 1700+518 and NGC 4051 using the spectral decomposition method. We focus on their optical Fe ii variability to investigate the origin of Fe ii in AGNs. For PG 1700+518, we find that the Fe ii size is about 200 light-days, which is consistent with the Hβ size derived from the empirical R–L relation. For NGC 4051, the [O iii] 5007 Å flux is strongly correlated with continuum flux, suggesting that we should recalibrate the spectral flux on a scale defined by [O iii] flux. The corrected light curves of Fe ii, Hβ, He ii, f _λ (5100 Å) are given here. A detailed analysis will be given in the near future.     

A Statistical Study on Photospheric Magnetic Nonpotentiality of Active Regions and Its Relationship with Flares During Solar Cycles 22?–?23

A statistical study is carried out on the photospheric magnetic nonpotentiality in solar active regions and its relationship with associated flares. We select 2173 photospheric vector magnetograms from 1106 active regions observed by the Solar Magnetic Field Telescope at Huairou Solar Observing Station, National Astronomical Observatories of China, in the period of 1988??C?2008, which covers most of the 22nd and 23rd solar cycles. We have computed the mean planar magnetic shear angle ( $\overline{\Delta\phi}$ ), mean shear angle of the vector magnetic field ( $\overline{\Delta\psi}$ ), mean absolute vertical current density ( $\overline{|J_{z}|}$ ), mean absolute current helicity density ( $\overline{|h_{\mathrm{c}}|}$ ), absolute twist parameter (|?? _av|), mean free magnetic energy density ( $\overline{\rho_{\mathrm{free}}}$ ), effective distance of the longitudinal magnetic field (d _E), and modified effective distance (d _Em) of each photospheric vector magnetogram. Parameters $\overline{|h_{\mathrm{c}}|}$ , $\overline{\rho_{\mathrm{free}}}$ , and d _Em show higher correlations with the evolution of the solar cycle. The Pearson linear correlation coefficients between these three parameters and the yearly mean sunspot number are all larger than 0.59. Parameters $\overline {\Delta\phi}$ , $\overline{\Delta\psi}$ , $\overline{|J_{z}|}$ , |?? _av|, and d _E show only weak correlations with the solar cycle, though the nonpotentiality and the complexity of active regions are greater in the activity maximum periods than in the minimum periods. All of the eight parameters show positive correlations with the flare productivity of active regions, and the combination of different nonpotentiality parameters may be effective in predicting the flaring probability of active regions.     

Optical spectral decomposition of NGC 4051

Considering the host galaxy contribution, a spectral decomposition method is used to reanalyzed the archive data of optical spectra for a narrow line Seyfert 1 galaxy, NGC 4051. The light curves of the continuum f _λ (5100 Å), and Hβ, He ii, Fe ii emission lines are given. We find strong flux correlations between line emissions of Hβ, He ii, Fe ii and the continuum f _λ (5100 Å). These low-ionization lines (Hβ, Fe ii, He ii) have “inverse” intrinsic Baldwin effects. Using the methods of the cross-correlation function and the Monte Carlo simulation, we find the time delays, with respect to the continuum, are $3.45^{+12.0}_{-0.5}~\mbox{days}$ with the probability of 34 % for the intermediate component of Hβ, $6.45^{+13.0}_{-1.0}~\mbox{days}$ with the probability of 65 % for the intermediate component of He ii. From these intermediate components of Hβ and He ii, the calculated central black hole masses are $0.86^{+4.35}_{-0.33}\times 10^{6}$ and $0.82^{+3.12}_{-0.45}\times 10^{6}~M_{\odot }$ . We also find that the time delays for Fe ii are $9.7^{+3.0}_{-5.0}~\mbox{days}$ with the probability of 36 %, $8.45^{+1.0}_{-2.0}~\mbox{days}$ with the probability of 18 % for the total epochs and “subset 1” data, respectively. It seems that the Fe ii emission region is outside of the Hβ emission region.