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.     

Variability of Solar Five-Minute Oscillations in the Corona as Observed by the Extreme Ultraviolet Spectrophotometer (ESP) on the Solar Dynamics Observatory/Extreme Ultraviolet Variability Experiment (SDO/EVE)

Solar five-minute oscillations have been detected in the power spectra of two six-day time intervals from soft X-ray measurements of the Sun observed as a star using the Extreme Ultraviolet Spectrophotometer (ESP) onboard the Solar Dynamics Observatory (SDO)/Extreme Ultraviolet Variability Experiment (EVE). The frequencies of the largest amplitude peaks were found to match the known low-degree (?=0?–?3) modes of global acoustic oscillations within 3.7 μHz and can be explained by a leakage of the global modes into the corona. Due to the strong variability of the solar atmosphere between the photosphere and the corona, the frequencies and amplitudes of the coronal oscillations are likely to vary with time. We investigated the variations in the power spectra for individual days and their association with changes of solar activity, e.g. with the mean level of the EUV irradiance, and its short-term variations caused by evolving active regions. Our analysis of samples of one-day oscillation power spectra for a 49-day period of low and intermediate solar activity showed little correlation with the mean EUV irradiance and the short-term variability of the irradiance. We suggest that some other changes in the solar atmosphere, e.g., magnetic fields and/or inter-network configuration may affect the mode leakage to the corona.     

Computer Vision for the Solar Dynamics Observatory (SDO)

In Fall 2008 NASA selected a large international consortium to produce a comprehensive automated feature-recognition system for the Solar Dynamics Observatory (SDO). The SDO data that we consider are all of the Atmospheric Imaging Assembly (AIA) images plus surface magnetic-field images from the Helioseismic and Magnetic Imager (HMI). We produce robust, very efficient, professionally coded software modules that can keep up with the SDO data stream and detect, trace, and analyze numerous phenomena, including flares, sigmoids, filaments, coronal dimmings, polarity inversion lines, sunspots, X-ray bright points, active regions, coronal holes, EIT waves, coronal mass ejections (CMEs), coronal oscillations, and jets. We also track the emergence and evolution of magnetic elements down to the smallest detectable features and will provide at least four full-disk, nonlinear, force-free magnetic field extrapolations per day. The detection of CMEs and filaments is accomplished with Solar and Heliospheric Observatory (SOHO)/Large Angle and Spectrometric Coronagraph (LASCO) and ground-based Hα data, respectively. A?completely new software element is a trainable feature-detection module based on a generalized image-classification algorithm. Such a trainable module can be used to find features that have not yet been discovered (as, for example, sigmoids were in the pre-Yohkoh era). Our codes will produce entries in the Heliophysics Events Knowledgebase (HEK) as well as produce complete catalogs for results that are too numerous for inclusion in the HEK, such as the X-ray bright-point metadata. This will permit users to locate data on individual events as well as carry out statistical studies on large numbers of events, using the interface provided by the Virtual Solar Observatory. The operations concept for our computer vision system is that the data will be analyzed in near real time as soon as they arrive at the SDO Joint Science Operations Center and have undergone basic processing. This will allow the system to produce timely space-weather alerts and to guide the selection and production of quicklook images and movies, in addition to its prime mission of enabling solar science. We briefly describe the complex and unique data-processing pipeline, consisting of the hardware and control software required to handle the SDO data stream and accommodate the computer-vision modules, which has been set up at the Lockheed-Martin Space Astrophysics Laboratory (LMSAL), with an identical copy at the Smithsonian Astrophysical Observatory (SAO).     

Study of Differences Between Sunspot and White Light Facular Area Data Determined from SDO/HMI and SOHO/MDI Observations

Sunspot and white light facular areas are important data for solar activity and are used, for example, in the study of the evolution of sunspots and their effect on solar irradiance. Solar Dynamic Observatory??s Helioseismic and Magnetic Imager (SDO/HMI) solar images have much higher resolution (??0.5????pixel^?1) than Solar and Heliospheric Observatory??s Michelson Doppler Imager (SOHO/MDI) solar images (??2????pixel^?1). This difference in image resolution has a significant impact on the sunspot and white light facular areas measured in the two datasets. We compare the area of sunspots and white light faculae derived from SDO/HMI and SOHO/MDI observations. This comparison helps the calibration of the SOHO sunspot and facular area to those in SDO observations. We also find a 0.22 degree difference between the North direction in SDO/HMI and SOHO/MDI images.     

EUV SpectroPhotometer (ESP) in Extreme Ultraviolet Variability Experiment (EVE): Algorithms and?Calibrations

The Extreme ultraviolet SpectroPhotometer (ESP) is one of five channels of the Extreme ultraviolet Variability Experiment (EVE) onboard the NASA Solar Dynamics Observatory (SDO). The ESP channel design is based on a highly stable diffraction transmission grating and is an advanced version of the Solar Extreme ultraviolet Monitor (SEM), which has been successfully observing solar irradiance onboard the Solar and Heliospheric Observatory (SOHO) since December 1995. ESP is designed to measure solar Extreme UltraViolet (EUV) irradiance in four first-order bands of the diffraction grating centered around 19 nm, 25 nm, 30 nm, and 36 nm, and in a soft X-ray band from 0.1 to 7.0?nm in?the?zeroth-order of the grating. Each band’s detector system converts the photo-current into a count rate (frequency). The count rates are integrated over 0.25-second increments and transmitted to the EVE Science and Operations Center for data processing. An algorithm for converting the measured count rates into solar irradiance and the ESP calibration parameters are described. The ESP pre-flight calibration was performed at the Synchrotron Ultraviolet Radiation Facility of the National Institute of Standards and Technology. Calibration parameters were used to calculate absolute solar irradiance from the sounding-rocket flight measurements on 14 April 2008. These irradiances for the ESP bands closely match the irradiance determined for two other EUV channels flown simultaneously: EVE’s Multiple EUV Grating Spectrograph (MEGS) and SOHO’s Charge, Element and Isotope Analysis System/Solar EUV Monitor (CELIAS/SEM).     

Automated Identification of Coronal Holes from Synoptic EUV Maps

Coronal holes (CHs) are regions of open magnetic field lines in the solar corona and the source of the fast solar wind. Understanding the evolution of coronal holes is critical for solar magnetism as well as for accurate space weather forecasts. We study the extreme ultraviolet (EUV) synoptic maps at three wavelengths (195 Å/193 Å, 171 Å and 304 Å) measured by the Solar and Heliospheric Observatory/Extreme Ultraviolet Imaging Telescope (SOHO/EIT) and the Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA) instruments. The two datasets are first homogenized by scaling the SDO/AIA data to the SOHO/EIT level by means of histogram equalization. We then develop a novel automated method to identify CHs from these homogenized maps by determining the intensity threshold of CH regions separately for each synoptic map. This is done by identifying the best location and size of an image segment, which optimally contains portions of coronal holes and the surrounding quiet Sun allowing us to detect the momentary intensity threshold. Our method is thus able to adjust itself to the changing scale size of coronal holes and to temporally varying intensities. To make full use of the information in the three wavelengths we construct a composite CH distribution, which is more robust than distributions based on one wavelength. Using the composite CH dataset we discuss the temporal evolution of CHs during the Solar Cycles 23 and 24.     

Modeling coronal magnetic field using spherical geometry: cases with several active regions

The magnetic fields in the solar atmosphere structure the plasma, store free magnetic energy and produce a wide variety of active solar phenomena, like flare and coronal mass ejections (CMEs). The distribution and strength of magnetic fields are routinely measured in the solar surface (photosphere). Therefore, there is considerable interest in accurately modeling the 3D structure of the coronal magnetic field using photospheric vector magnetograms. Knowledge of the 3D structure of magnetic field lines also help us to interpret other coronal observations, e.g., EUV images of the radiating coronal plasma. Nonlinear force-free field (NLFFF) models are thought to be viable tools for those task. Usually those models use Cartesian geometry. However, the spherical nature of the solar surface cannot be neglected when the field of view is large. In this work, we model the coronal magnetic field above multiple active regions using NLFFF extrapolation code using vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS)/Vector Spectromagnetograph (VSM) as a boundary conditions. We compare projections of the resulting magnetic field lines solutions with their respective coronal EUV-images from the Atmospheric Imaging Assembly (SDO/AIA) observed on October 15, 2011 and November 13, 2012. This study has found that the NLFFF model in spherical geometry reconstructs the magnetic configurations for several active regions which agrees to some extent with observations. During October 15, 2011 observation, there are substantial number of trans-equatorial loops carrying electric current.     

Coronal Temperature Maps from Solar EUV Images: A Blind Source Separation Approach

Multi-wavelength solar images in the extreme ultraviolet (EUV) are routinely used for analysing solar features such as coronal holes, filaments, and flares. However, images taken in different bands often look remarkably similar, as each band receives contributions coming from regions with a range of different temperatures. This has motivated the search for empirical techniques that may unmix these contributions and concentrate salient morphological features of the corona in a smaller set of less redundant source images. Blind Source Separation (BSS) does precisely this. Here we show how this novel concept also provides new insight into the physics of the solar corona, using observations made by SDO/AIA. The source images are extracted using a Bayesian positive source-separation technique. We show how observations made in six spectral bands, corresponding to optically thin emissions, can be reconstructed by a linear combination of three sources. These sources have a narrower temperature response and allow for considerable data reduction, since the pertinent information from all six bands can be condensed into a single composite picture. In addition, they give access to empirical temperature maps of the corona. The limitations of the BSS technique and some applications are briefly discussed.     

Coronal Holes and Open Magnetic Flux over Cycles 23 and 24

As the observational signature of the footprints of solar magnetic field lines open into the heliosphere, coronal holes provide a critical measure of the structure and evolution of these lines. Using a combination of Solar and Heliospheric Observatory/Extreme ultraviolet Imaging Telescope (SOHO/EIT), Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA), and Solar Terrestrial Relations Observatory/Extreme Ultraviolet Imager (STEREO/EUVI A/B) extreme ultraviolet (EUV) observations spanning 1996?–?2015 (nearly two solar cycles), coronal holes are automatically detected and characterized. Coronal hole area distributions show distinct behavior in latitude, defining the domain of polar and low-latitude coronal holes. The northern and southern polar regions show a clear asymmetry, with a lag between hemispheres in the appearance and disappearance of polar coronal holes.     

Automated Coronal Hole Detection Using Local Intensity Thresholding Techniques

We identify coronal holes using a histogram-based intensity thresholding technique and compare their properties to fast solar wind streams at three different points in the heliosphere. The thresholding technique was tested on EUV and X-ray images obtained using instruments onboard STEREO, SOHO and Hinode. The full-disk images were transformed into Lambert equal-area projection maps and partitioned into a series of overlapping sub-images from which local histograms were extracted. The histograms were used to determine the threshold for the low intensity regions, which were then classified as coronal holes or filaments using magnetograms from the SOHO/MDI. For all three instruments, the local thresholding algorithm was found to successfully determine coronal hole boundaries in a consistent manner. Coronal hole properties extracted using the segmentation algorithm were then compared with in situ measurements of the solar wind at ～?1 AU from ACE and STEREO. Our results indicate that flux tubes rooted in coronal holes expand super-radially within 1 AU and that larger (smaller) coronal holes result in longer (shorter) duration high-speed solar wind streams.     

Using Polar Coronal Hole Area Measurements to Determine the Solar Polar Magnetic Field Reversal in Solar Cycle 24

An analysis of solar polar coronal hole (PCH) areas since the launch of the Solar Dynamics Observatory (SDO) shows how the polar regions have evolved during Solar Cycle 24. We present PCH areas from mid-2010 through 2013 using data from the Atmospheric Imager Assembly (AIA) and Helioseismic and Magnetic Imager (HMI) instruments onboard SDO. Our analysis shows that both the northern and southern PCH areas have decreased significantly in size since 2010. Linear fits to the areas derived from the magnetic-field properties indicate that, although the northern hemisphere went through polar-field reversal and reached solar-maximum conditions in mid-2012, the southern hemisphere had not reached solar-maximum conditions in the polar regions by the end of 2013. Our results show that solar-maximum conditions in each hemisphere, as measured by the area of the polar coronal holes and polar magnetic field, will be offset in time.     

Rotation of solar structures in the upper chromosphere. II. Time variations in the latitudinal distribution of the rotation of active regions and coronal holes

The time variations in the latitudinal distribution of the rotation of active regions and coronal holes are investigated. The synoptic maps obtained from observations in the He I 1083 nm line at Kitt Peak Observatory over almost three solar cycles are used as observational data. A Fourier analysis of the time series constructed from synoptic maps has yielded the following results. The rotation of active regions differs significantly from the rotation of coronal holes in all parameters: the set of the most significant rotation periods, their latitudinal distribution, and time variations. The rotation of active regions and coronal holes is characterized by variations from cycle to cycle, a time-varying north-south asymmetry. The power spectra for consecutive cycles of solar activity differ significantly for both epochs of high activity and minima. Analysis of the total power of the spectra within four selected intervals of periods from 21 to 33 days has shown that the total power is highest in the intervals of periods 24–27 and 27–30 days. This is valid for both active regions and coronal holes. The correlation between the total powers in the above intervals of periods changes noticeably with time. Long-lived or successively appearing active regions with rotation periods in the range 24–30 days are typical of the time of a sharp decrease in the total equivalent width of active regions. This includes not only the decline time of the 11-year cycles, but also the minima between recurrent activity maxima during one cycle. A predominance of long-lived coronal holes as their total equivalent width decreases is noticeable for coronal holes with rotation periods in the interval 30–33 days. All of the above results suggest that the rotation of solar structures is determined mainly by the subphotospheric sources of specific structures, not by the rotation of the main volumes of solar plasma of the quiet Sun.     

Solar Si XI Line Ratios Observed by the Normal Incidence Spectrometer on SOHO CDS

New measurements of line intensity ratios in the Be-like ion Si xi are presented for observations of the quiet Sun, active regions, coronal holes and above-limb regions obtained using the Coronal Diagnostic Spectrometer on the Solar and Heliospheric Observatory. A model ion, constructed using the best available atomic data, is used to predict the line intensity ratios for a wide range of electron temperatures and densities. Comparisons of the theoretical ratios with the new intensity ratios as well as with those from previous solar observations and laboratory measurements are given. The usefulness of the ratios for electron temperature and density diagnostics, as well as for spectrometer calibration, is discussed.     

Synoptic charts of solar 9.1 cm and coronal hole data

Synoptic charts for Carrington rotations 1601–1605 (May–August, 1973) were prepared using the central meridian column of the daily 9.1 cm Stanford solar radio maps. These charts were especially contoured to emphasize temperatures near the quiet solar disk level. Synoptic charts of coronal holes from the ATM-Skylab were superimposed on the radio data to investigate the ability of the radio charts to show coronal holes. This brief period is unfortunately the only interval for which both sets of data are available. The conclusion reached is that in spite of certain problems due to active regions, side-lobe effects and a rather large beamwidth, the 9.1 cm synoptic charts can be of substantial value in identifying large coronal holes, especially during periods of low solar activity. Such synoptic charts, therefore, for the years 1962–1973 that Stanford data are available, could enhance significantly the meagre data pool for coronal holes prior to the Skylab mission.     

Structure and Dynamics of Interconnecting Loops and Coronal Holes in Active Longitudes

Using SOHO/MDI and SOHO/EIT data we study properties and dynamics of interconnected active regions, and the relations between the photospheric magnetic fields and coronal structures in active longitudes during the beginning of solar cycle 23. The emergence of new magnetic flux results in appearance of new interconnecting loops. The existence of stable coronal structures strongly depends on the photospheric magnetic fluxes and their variations. We present some initial results for a complex of solar activity observed in April 1997, and discuss the role of reconnection in the formation of the interconnected loops and coronal holes.     

Open magnetic fields and the solar cycle

Models of open magnetic structures on the Sun are presented for periods near solar minimum (CR 1626–1634) and near solar maximum (CR 1668–1678). Together with previous models of open magnetic structures during the declining phase (CR 1601–1611) these calculations provide clues to the relations between open structures, coronal holes, and active regions at different times of the solar cycle. Near solar minimum the close relation between active regions and open structures does not exist. It is suggested that near solar minimum the systematic emergence of new flux with the proper polarity imbalance to maintain open magnetic structures may occur primarily at very small spatial scales. Near solar maximum the role of active regions in maintaining open structures and coronal holes is strong, with large active regions emerging in the proper location and orientation to maintain open structures longer than typical active region lifetimes. Although the use of He I 10830 Å spectroheliograms as a coronal hole indicator is shown to be subject to significant ambiguity, the agreement between calculated open structures and coronal holes determined from He I 10830 Å spectroheliograms is very good. The rotation properties of calculated open structures near solar maximum strongly suggest two classes of features: one that rotates differentially similar to sunspots and active regions and a separate class that rotates more rigidly, as was the case for single large coronal holes during Skylab.     

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.     

Fast Segmentation of Solar Extreme Ultraviolet Images

A segmentation scheme for identifying large-scale structures (coronal holes, active regions, etc.) in solar extreme ultraviolet images, is presented. Unlike standard approaches, both the image intensity and the relative contribution of different wavelengths are used. Spectral information is important for compensating luminosity changes. The approach is illustrated with images taken in the extreme ultraviolet by the EIT telescope onboard SOHO. This supervised segmentation scheme, which incorporates a Bayesian classifier, is computationally simple, and can easily be used to track in near-real time structures, such as coronal holes.     

The relationship between solar activity and coronal hole evolution

We examine the relationship between coronal hole evolution and solar active regions during the Skylab period. We find a tendency for holes to grow or remain stable when the activity nearby, seen as calcium plages and bright regions in X-rays, is predominantly large, long-lived regions. This is consistent with results of previous studies, using somewhat different methods. We also find that there is a significantly higher number of small, short-lived active regions, as indicated by X-ray bright points, in the vicinity of decaying holes than there is near other holes. We interpret this to mean that holes disappear at least in part because they become filled with many small scale, magnetically closed, X-ray emitting features. This interpretation, together with the previously reported observation that the number of X-ray bright points was much larger near solar minimum than it was during the Skylab period, provides a possible explanation for the disappearance of the large, near-equatorial coronal holes at the time of solar minimum.     

How Should One Optimize Nonlinear Force-Free Coronal Magnetic Field Extrapolations from SDO/HMI Vector Magnetograms?

The Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) provides photospheric vector magnetograms with a high spatial and temporal resolution. Our intention is to model the coronal magnetic field above active regions with the help of a nonlinear force-free extrapolation code. Our code is based on an optimization principle and has been tested extensively with semianalytic and numeric equilibria and applied to vector magnetograms from Hinode and ground-based observations. Recently we implemented a new version which takes into account measurement errors in photospheric vector magnetograms. Photospheric field measurements are often affected by measurement errors and finite nonmagnetic forces inconsistent for use as a boundary for a force-free field in the corona. To deal with these uncertainties, we developed two improvements: i) preprocessing of the surface measurements to make them compatible with a force-free field, and ii) new code which keeps a balance between the force-free constraint and deviation from the photospheric field measurements. Both methods contain free parameters, which must be optimized for use with data from SDO/HMI. In this work we describe the corresponding analysis method and evaluate the force-free equilibria by how well force-freeness and solenoidal conditions are fulfilled, by the angle between magnetic field and electric current, and by comparing projections of magnetic field lines with coronal images from the Atmospheric Imaging Assembly (SDO/AIA). We also compute the available free magnetic energy and discuss the potential influence of control parameters.