-
Forecasting SEP Events During Solar Cycles 23 and 24 Using Interpretable Machine Learning
Authors:
Spiridon Kasapis,
Irina N. Kitiashvili,
Paul Kosovich,
Alexander G. Kosovichev,
Viacheslav M. Sadykov,
Patrick O'Keefe,
Vincent Wang
Abstract:
Prediction of the Solar Energetic Particle (SEP) events garner increasing interest as space missions extend beyond Earth's protective magnetosphere. These events, which are, in most cases, products of magnetic reconnection-driven processes during solar flares or fast coronal-mass-ejection-driven shock waves, pose significant radiation hazards to aviation, space-based electronics, and particularly,…
▽ More
Prediction of the Solar Energetic Particle (SEP) events garner increasing interest as space missions extend beyond Earth's protective magnetosphere. These events, which are, in most cases, products of magnetic reconnection-driven processes during solar flares or fast coronal-mass-ejection-driven shock waves, pose significant radiation hazards to aviation, space-based electronics, and particularly, space exploration. In this work, we utilize the recently developed dataset that combines the Solar Dynamics Observatory/Helioseismic and Magnetic Imager's (SDO/HMI) Space weather HMI Active Region Patches (SHARP) and the Solar and Heliospheric Observatory/Michelson Doppler Imager's (SoHO/MDI) Space Weather MDI Active Region Patches (SMARP). We employ a suite of machine learning strategies, including Support Vector Machines (SVM) and regression models, to evaluate the predictive potential of this new data product for a forecast of post-solar flare SEP events. Our study indicates that despite the augmented volume of data, the prediction accuracy reaches 0.7 +- 0.1, which aligns with but does not exceed these published benchmarks. A linear SVM model with training and testing configurations that mimic an operational setting (positive-negative imbalance) reveals a slight increase (+ 0.04 +- 0.05) in the accuracy of a 14-hour SEP forecast compared to previous studies. This outcome emphasizes the imperative for more sophisticated, physics-informed models to better understand the underlying processes leading to SEP events.
△ Less
Submitted 4 March, 2024;
originally announced March 2024.
-
Investigating Performance Trends of Simulated Real-time Solar Flare Predictions: The Impacts of Training Windows, Data Volumes, and the Solar Cycle
Authors:
Griffin T. Goodwin,
Viacheslav M. Sadykov,
Petrus C. Martens
Abstract:
This study explores the behavior of machine learning-based flare forecasting models deployed in a simulated operational environment. Using Georgia State University's Space Weather Analytics for Solar Flares benchmark dataset (Angryk et al. 2020a,b), we examine the impacts of training methodology and the solar cycle on decision tree, support vector machine, and multilayer perceptron performance. We…
▽ More
This study explores the behavior of machine learning-based flare forecasting models deployed in a simulated operational environment. Using Georgia State University's Space Weather Analytics for Solar Flares benchmark dataset (Angryk et al. 2020a,b), we examine the impacts of training methodology and the solar cycle on decision tree, support vector machine, and multilayer perceptron performance. We implement our classifiers using three temporal training windows: stationary, rolling, and expanding. The stationary window trains models using a single set of data available before the first forecasting instance, which remains constant throughout the solar cycle. The rolling window trains models using data from a constant time interval before the forecasting instance, which moves with the solar cycle. Finally, the expanding window trains models using all available data before the forecasting instance. For each window, a number of input features (1, 5, 10, 25, 50, 120) and temporal sizes (5, 8, 11, 14, 17, 20 months) were tested. To our surprise, we found that for a 20-month window, skill scores were comparable regardless of the window type, feature count, and classifier selected. Furthermore, reducing the size of this window only marginally decreased stationary and rolling window performance. This implies that, given enough data, a stationary window can be chosen over other window types, eliminating the need for model retraining. Lastly, a moderately strong positive correlation was found to exist between a model's false positive rate and the solar X-ray background flux. This suggests that the solar cycle phase has a considerable influence on forecasting.
△ Less
Submitted 7 February, 2024;
originally announced February 2024.
-
Subsurface Flows Associated with Formation and Flaring Activity of Solar Active Regions
Authors:
Alexander G. Kosovichev,
Viacheslav M. Sadykov
Abstract:
We investigate the evolution of subsurface flows during the emergence and the active phase of sunspot regions using the time-distance helioseismology analysis of the full-disk Dopplergrams from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We present an analysis of emerging active regions of various types, including delta-type active regions and regions w…
▽ More
We investigate the evolution of subsurface flows during the emergence and the active phase of sunspot regions using the time-distance helioseismology analysis of the full-disk Dopplergrams from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We present an analysis of emerging active regions of various types, including delta-type active regions and regions with the reverse polarity order (`anti-Hale active regions'). The results reveal strong vortical and shearing flows during the emergence of magnetic flux, as well as the process of formation of large-scale converging flow patterns around developing active regions, predominantly in the top 6 Mm deep layers of the convection zone. Our analysis revealed a significant correlation between the flow divergence and helicity in the active regions with their flaring activity, indicating that measuring characteristics of subsurface flows can contribute to flare forecasting.
△ Less
Submitted 30 January, 2024;
originally announced January 2024.
-
Time Series of Magnetic Field Parameters of Merged MDI and HMI Space-Weather Active Region Patches as Potential Tool for Solar Flare Forecasting
Authors:
Paul A. Kosovich,
Viacheslav M. Sadykov,
Alexander G. Kosovichev,
Spiridon Kasapis,
Irina N. Kitiashvili,
Patrick M. O'Keefe,
Aatiya Ali,
Vincent Oria,
Samuel Granovsky,
Chun Jie Chong,
Gelu M. Nita
Abstract:
Solar flare prediction studies have been recently conducted with the use of Space-Weather MDI (Michelson Doppler Imager onboard Solar and Heliospheric Observatory) Active Region Patches (SMARP) and Space-Weather HMI (Helioseismic and Magnetic Imager onboard Solar Dynamics Observatory) Active Region Patches (SHARP), which are two currently available data products containing magnetic field character…
▽ More
Solar flare prediction studies have been recently conducted with the use of Space-Weather MDI (Michelson Doppler Imager onboard Solar and Heliospheric Observatory) Active Region Patches (SMARP) and Space-Weather HMI (Helioseismic and Magnetic Imager onboard Solar Dynamics Observatory) Active Region Patches (SHARP), which are two currently available data products containing magnetic field characteristics of solar active regions. The present work is an effort to combine them into one data product, and perform some initial statistical analyses in order to further expand their application in space weather forecasting. The combined data are derived by filtering, rescaling, and merging the SMARP with SHARP parameters, which can then be spatially reduced to create uniform multivariate time series. The resulting combined MDI-HMI dataset currently spans the period between April 4, 1996, and December 13, 2022, and may be extended to a more recent date. This provides an opportunity to correlate and compare it with other space weather time series, such as the daily solar flare index or the statistical properties of the soft X-ray flux measured by the Geostationary Operational Environmental Satellites (GOES). Time-lagged cross-correlation indicates that a relationship may exist, where some magnetic field properties of active regions lead the flare index in time. Applying the rolling window technique makes it possible to see how this leader-follower dynamic varies with time. Preliminary results indicate that areas of high correlation generally correspond to increased flare activity during the peak solar cycle.
△ Less
Submitted 22 February, 2024; v1 submitted 10 January, 2024;
originally announced January 2024.
-
Spectro-Polarimetric Properties of Sunquake Sources in X1.5 Flare and Evidence for Electron and Proton Beam Impacts
Authors:
Alexander G. Kosovichev,
Viacheslav M. Sadykov,
John T. Stefan
Abstract:
The first significant sunquake event of Solar Cycle 25 was observed during the X1.5 flare of May 10, 2022, by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory. We perform a detailed spectro-polarimetric analysis of the sunquake photospheric sources, using the Stokes profiles of the FeI 6173A line, reconstructed from the HMI linear and circular polarized filtergrams…
▽ More
The first significant sunquake event of Solar Cycle 25 was observed during the X1.5 flare of May 10, 2022, by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory. We perform a detailed spectro-polarimetric analysis of the sunquake photospheric sources, using the Stokes profiles of the FeI 6173A line, reconstructed from the HMI linear and circular polarized filtergrams. The results show fast variations of the continuum emission with rapid growth and slower decay lasting 3-4 min, coinciding in time with the hard X-ray impulses observed by the Konus instrument onboard the Wind spacecraft. The variations in the line core appeared slightly ahead of the variations in the line wings, showing that the heating started in the higher atmospheric layers and propagated downward. The most significant feature of the line profile variations is the transient emission in the line core in three of the four sources, indicating intense, impulsive heating in the lower chromosphere and photosphere. In addition, the observed variations of the Stokes profiles reflect transient and permanent changes in the magnetic field strength and geometry in the sunquake sources. Comparison with the radiative hydrodynamics models shows that the physical processes in the impulsive flare phase are substantially more complex than those predicted by proton and electron beam flare models currently presented in the literature.
△ Less
Submitted 13 September, 2023;
originally announced September 2023.
-
Solar Flare Ribbon Fronts I: Constraining flare energy deposition with IRIS spectroscopy
Authors:
Vanessa Polito,
Graham S. Kerr,
Yan Xu,
Viacheslav M. Sadykov,
Juraj Lorincik
Abstract:
Lower atmospheric lines show peculiar profiles at the leading edge of ribbons during solar flares. In particular, increased absorption of the BBSO/GST \hei~10830~Å line \citep[e.g.][]{Xu2016}, as well as broad and centrally reversed profiles in the spectra of the \mgii~and \cii~lines observed by the \iris~satellite \citep[e.g.][]{Panos2018,Panos2021a} have been reported. In this work, we aim to un…
▽ More
Lower atmospheric lines show peculiar profiles at the leading edge of ribbons during solar flares. In particular, increased absorption of the BBSO/GST \hei~10830~Å line \citep[e.g.][]{Xu2016}, as well as broad and centrally reversed profiles in the spectra of the \mgii~and \cii~lines observed by the \iris~satellite \citep[e.g.][]{Panos2018,Panos2021a} have been reported. In this work, we aim to understand the physical origin of the \iris\ ribbon front line profiles, which seem to be common of many, if not all, flares. To achieve this, we quantify the spectral properties of the \iris~\mgii~ribbon front profiles during four large flares and perform a detailed comparison with a grid of radiative hydrodynamic models using the \radynfp~code. We also studied their transition region counterparts, finding that these ribbon front locations are regions where transition region emission and chromospheric evaporation are considerably weaker compared to other parts of the ribbons. Based on our comparison between the \iris~observations and modelling, our interpretation is that there are different heating regimes at play in the leading and trailing regions of the ribbons. More specifically, we suggest that bombardment of the chromosphere by more gradual and modest non-thermal electron energy fluxes can qualitatively explain the \iris~observations at the ribbon front, while stronger and more impulsive energy fluxes are required to drive chromospheric evaporation and more intense TR emission. Our results provide a possible physical origin for the peculiar behaviour of the \iris~chromospheric lines in the ribbon leading edge and new constraints for the flare models.
△ Less
Submitted 9 November, 2022;
originally announced November 2022.
-
Leptocline as a Shallow Substructure of Near-Surface Shear Layer in 3D Radiative Hydrodynamic Simulations
Authors:
Irina N. Kitiashvili,
Alexander G. Kosovichev,
Alan A. Wray,
Viacheslav M. Sadykov,
Gustavo Guerrero
Abstract:
Understanding effects driven by rotation in the solar convection zone is essential for many problems related to solar activity, such as the formation of differential rotation, meridional circulation, and others. We analyze realistic 3D radiative hydrodynamics simulations of solar subsurface dynamics in the presence of rotation in a local domain 80 Mm wide and 25 Mm deep, located at 30 degrees lati…
▽ More
Understanding effects driven by rotation in the solar convection zone is essential for many problems related to solar activity, such as the formation of differential rotation, meridional circulation, and others. We analyze realistic 3D radiative hydrodynamics simulations of solar subsurface dynamics in the presence of rotation in a local domain 80 Mm wide and 25 Mm deep, located at 30 degrees latitude. The simulation results reveal the development of a shallow 10-Mm deep substructure of the Near-Surface Shear Layer (NSSL), characterized by a strong radial rotational gradient and self-organized meridional flows. This shallow layer ("leptocline") is located in the hydrogen ionization zone associated with enhanced anisotropic overshooting-type flows into a less unstable layer between the H and HeII ionization zones. We discuss current observational evidence of the presence of the leptocline and show that the radial variations of the differential rotation and meridional flow profiles obtained from the simulations in this layer qualitatively agree with helioseismic observations.
△ Less
Submitted 27 August, 2022; v1 submitted 2 March, 2022;
originally announced March 2022.
-
Multi-Passband Observations of A Solar Flare over the He I 10830 Å line
Authors:
Yan Xu,
Xu Yang,
Graham S. Kerr,
Vanessa Polito,
Viacheslav M. Sadykov,
Ju Jing,
Wenda Cao,
Haimin Wang
Abstract:
This study presents a C3.0 flare observed by the BBSO/GST and IRIS, on 2018-May-28 around 17:10 UT. The Near Infrared Imaging Spectropolarimeter (NIRIS) of GST was set to spectral imaging mode to scan five spectral positions at $\pm$ 0.8 Å, $\pm$ 0.4 Åand line center of He I 10830. At the flare ribbon's leading edge the line is observed to undergo enhanced absorption, while the rest of the ribbon…
▽ More
This study presents a C3.0 flare observed by the BBSO/GST and IRIS, on 2018-May-28 around 17:10 UT. The Near Infrared Imaging Spectropolarimeter (NIRIS) of GST was set to spectral imaging mode to scan five spectral positions at $\pm$ 0.8 Å, $\pm$ 0.4 Åand line center of He I 10830. At the flare ribbon's leading edge the line is observed to undergo enhanced absorption, while the rest of the ribbon is observed to be in emission. When in emission, the contrast compared to the pre-flare ranges from about $30~\%$ to nearly $100~\%$ at different spectral positions. Two types of spectra, "convex" shape with higher intensity at line core and "concave" shape with higher emission in the line wings, are found at the trailing and peak flaring areas, respectively. On the ribbon front, negative contrasts, or enhanced absorption, of about $\sim 10\% - 20\%$ appear in all five wavelengths. This observation strongly suggests that the negative flares observed in He I 10830 with mono-filtergram previously were not caused by pure Doppler shifts of this spectral line. Instead, the enhanced absorption appears to be a consequence of flare energy injection, namely non-thermal collisional ionization of helium caused by the precipitation of high energy electrons, as found in our recent numerical modeling results. In addition, though not strictly simultaneous, observations of Mg II from the IRIS spacecraft, show an obvious central reversal pattern at the locations where enhanced absorption of He I 10830 is seen, which is in consistent with previous observations.
△ Less
Submitted 18 December, 2021;
originally announced December 2021.
-
Feature Selection on a Flare Forecasting Testbed: A Comparative Study of 24 Methods
Authors:
Atharv Yeoleka,
Sagar Patel,
Shreejaa Talla,
Krishna Rukmini Puthucode,
Azim Ahmadzadeh,
Viacheslav M. Sadykov,
Rafal A. Angryk
Abstract:
The Space-Weather ANalytics for Solar Flares (SWAN-SF) is a multivariate time series benchmark dataset recently created to serve the heliophysics community as a testbed for solar flare forecasting models. SWAN-SF contains 54 unique features, with 24 quantitative features computed from the photospheric magnetic field maps of active regions, describing their precedent flare activity. In this study,…
▽ More
The Space-Weather ANalytics for Solar Flares (SWAN-SF) is a multivariate time series benchmark dataset recently created to serve the heliophysics community as a testbed for solar flare forecasting models. SWAN-SF contains 54 unique features, with 24 quantitative features computed from the photospheric magnetic field maps of active regions, describing their precedent flare activity. In this study, for the first time, we systematically attacked the problem of quantifying the relevance of these features to the ambitious task of flare forecasting. We implemented an end-to-end pipeline for preprocessing, feature selection, and evaluation phases. We incorporated 24 Feature Subset Selection (FSS) algorithms, including multivariate and univariate, supervised and unsupervised, wrappers and filters. We methodologically compared the results of different FSS algorithms, both on the multivariate time series and vectorized formats, and tested their correlation and reliability, to the extent possible, by using the selected features for flare forecasting on unseen data, in univariate and multivariate fashions. We concluded our investigation with a report of the best FSS methods in terms of their top-k features, and the analysis of the findings. We wish the reproducibility of our study and the availability of the data allow the future attempts be comparable with our findings and themselves.
△ Less
Submitted 29 September, 2021;
originally announced September 2021.
-
He I 10830 Å Dimming During Solar Flares, I. The Crucial Role of Non-Thermal Collisional Ionisations
Authors:
Graham S. Kerr,
Yan Xu,
Joel C. Allred,
Vanessa Polito,
Viacheslav M. Sadykov,
Nengyi Huang,
Haimin Wang
Abstract:
While solar flares are predominantly characterised by an intense broadband enhancement to the solar radiative output, certain spectral lines and continua will, in theory, exhibit flare-induced dimmings. Observations of transitions of orthohelium He I $λλ10830$Å and the He I D3 lines have shown evidence of such dimming, usually followed by enhanced emission. It has been suggested that non-thermal c…
▽ More
While solar flares are predominantly characterised by an intense broadband enhancement to the solar radiative output, certain spectral lines and continua will, in theory, exhibit flare-induced dimmings. Observations of transitions of orthohelium He I $λλ10830$Å and the He I D3 lines have shown evidence of such dimming, usually followed by enhanced emission. It has been suggested that non-thermal collisional ionisation of helium by an electron beam, followed by recombinations to orthohelium, is responsible for overpopulating the those levels, leading to stronger absorption. However it has not been possible observationally to preclude the possibility of overpopulating orthohelium via enhanced photoionisation of He I by EUV irradiance from the flaring corona followed by recombinations. Here we present radiation hydrodynamics simulations of non-thermal electron beam-driven flares where (1) both non-thermal collisional ionisation of Helium and coronal irradiance are included, and (2) only coronal irradiance is included. A grid of simulations covering a range of total energies deposited by the electron beam, and a range of non-thermal electron beam low-energy cutoff values, were simulated. In order to obtain flare-induced dimming of the He I 10830Å line it was necessary for non-thermal collisional ionisations to be present. The effect was more prominent in flares with larger low-energy cutoff values and longer lived in weaker flares and flares with a more gradual energy deposition timescale. These results demonstrate the usefulness of orthohelium line emission as a diagnostic of flare energy transport.
△ Less
Submitted 23 May, 2021; v1 submitted 30 March, 2021;
originally announced March 2021.
-
Radiation Data Portal: Integration of Radiation Measurements at the Aviation Altitudes and Solar-Terrestrial Environment Observations
Authors:
Viacheslav M Sadykov,
Irina N Kitiashvili,
W Kent Tobiska,
Madhulika Guhathakurta
Abstract:
The impact of radiation dramatically increases at high altitudes in the Earth's atmosphere and in space. Therefore, monitoring and access to radiation environment measurements are critical for estimating the radiation exposure risks of aircraft and spacecraft crews and the impact of space weather disturbances on electronics. Addressing these needs requires convenient access to multi-source radiati…
▽ More
The impact of radiation dramatically increases at high altitudes in the Earth's atmosphere and in space. Therefore, monitoring and access to radiation environment measurements are critical for estimating the radiation exposure risks of aircraft and spacecraft crews and the impact of space weather disturbances on electronics. Addressing these needs requires convenient access to multi-source radiation environment data and enhancement of visualization and search capabilities. The Radiation Data Portal represents an interactive web-based application for search and visualization of in-flight radiation measurements. The Portal enhances the exploration capabilities of various properties of the radiation environment and provides measurements of dose rates along with information on space weather-related conditions. The Radiation Data Portal back-end is a MySQL relational database that contains the radiation measurements obtained from the Automated Radiation Measurements for Aerospace Safety (ARMAS) device and the soft X-ray and proton flux measurements from the Geostationary Operational Environmental Satellite (GOES). The implemented Application Programming Interface (API) and Python routines allow a user to retrieve the database records without interaction with the web interface. As a use case of the Radiation Data Portal, we examine ARMAS measurements during an enhancement of the Solar Proton (SP) fluxes, known as Solar Proton Events (SPEs), and compare them to measurements during SP-quiet periods.
△ Less
Submitted 12 March, 2021;
originally announced March 2021.
-
Compression of Solar Spectroscopic Observations: a Case Study of Mg II k Spectral Line Profiles Observed by NASA's IRIS Satellite
Authors:
Viacheslav M Sadykov,
Irina N Kitiashvili,
Alberto Sainz Dalda,
Vincent Oria,
Alexander G Kosovichev,
Egor Illarionov
Abstract:
In this study we extract the deep features and investigate the compression of the Mg II k spectral line profiles observed in quiet Sun regions by NASA's IRIS satellite. The data set of line profiles used for the analysis was obtained on April 20th, 2020, at the center of the solar disc, and contains almost 300,000 individual Mg II k line profiles after data cleaning. The data are separated into tr…
▽ More
In this study we extract the deep features and investigate the compression of the Mg II k spectral line profiles observed in quiet Sun regions by NASA's IRIS satellite. The data set of line profiles used for the analysis was obtained on April 20th, 2020, at the center of the solar disc, and contains almost 300,000 individual Mg II k line profiles after data cleaning. The data are separated into train and test subsets. The train subset was used to train the autoencoder of the varying embedding layer size. The early stopping criterion was implemented on the test subset to prevent the model from overfitting. Our results indicate that it is possible to compress the spectral line profiles more than 27 times (which corresponds to the reduction of the data dimensionality from 110 to 4) while having a 4 DN average reconstruction error, which is comparable to the variations in the line continuum. The mean squared error and the reconstruction error of even statistical moments sharply decrease when the dimensionality of the embedding layer increases from 1 to 4 and almost stop decreasing for higher numbers. The observed occasional improvements in training for values higher than 4 indicate that a better compact embedding may potentially be obtained if other training strategies and longer training times are used. The features learned for the critical four-dimensional case can be interpreted. In particular, three of these four features mainly control the line width, line asymmetry, and line dip formation respectively. The presented results are the first attempt to obtain a compact embedding for spectroscopic line profiles and confirm the value of this approach, in particular for feature extraction, data compression, and denoising.
△ Less
Submitted 16 June, 2021; v1 submitted 12 March, 2021;
originally announced March 2021.
-
Connecting Atmospheric Properties and Synthetic Emission of Shock Waves Using 3D RMHD Simulations of Quiet Sun
Authors:
Viacheslav M. Sadykov,
Irina N. Kitiashvili,
Alexander G. Kosovichev,
Alan A. Wray
Abstract:
We analyze the evolution of shock waves in high-resolution 3D radiative MHD simulations of the quiet Sun and their synthetic emission characteristics. The simulations model the dynamics of a 12.8x12.8x15.2 Mm quiet-Sun region (including a 5.2 Mm layer of the upper convection zone and a 10 Mm atmosphere from the photosphere to corona) with an initially uniform vertical magnetic field of 10 G, natur…
▽ More
We analyze the evolution of shock waves in high-resolution 3D radiative MHD simulations of the quiet Sun and their synthetic emission characteristics. The simulations model the dynamics of a 12.8x12.8x15.2 Mm quiet-Sun region (including a 5.2 Mm layer of the upper convection zone and a 10 Mm atmosphere from the photosphere to corona) with an initially uniform vertical magnetic field of 10 G, naturally driven by convective flows. We synthesize the Mg II and C II spectral lines observed by the IRIS satellite and EUV emission observed by the SDO/AIA telescope. Synthetic observations are obtained using the RH1.5D radiative transfer code and temperature response functions at both the numerical and instrumental resolutions. We found that the Doppler velocity jumps of the C II 1334.5 A IRIS line and a relative enhancement of the emission in the 335 A SDO/AIA channel are the best proxies for the enthalpy deposited by shock waves into the corona (with Kendall's $τ$ correlation coefficients of 0.59 and 0.38, respectively). The synthetic emission of the lines and extreme ultraviolet passbands are correlated with each other during the shock wave propagation. All studied shocks are mostly hydrodynamic (i.e., the magnetic energy carried by horizontal fields is $\leq{}$2.6% of the enthalpy for all events) and have Mach numbers > 1.0-1.2 in the low corona. The study reveals the possibility of diagnosing energy transport by shock waves into the solar corona, as well as their other properties, by using IRIS and SDO/AIA sensing observations.
△ Less
Submitted 12 March, 2021; v1 submitted 13 August, 2020;
originally announced August 2020.
-
Comparison of Enhanced Absorption in He I 10830 Å in Observations and Modeling During the Early Phase of a Solar Flare
Authors:
Nengyi Huang,
Viacheslav M. Sadykov,
Yan Xu,
Ju Jing,
Haimin Wang
Abstract:
The He I 10830 Å triplet is a very informative indicator of chromospheric activities as the helium is the second most abundant element in the solar atmosphere. Taking advantage of the high resolution of the 1.6 m Goode Solar Telescope (GST) at Big Bear Solar Observatory (BBSO), previous observations have shown clear evidence of the enhanced absorption, instead of typically-observed emission, for t…
▽ More
The He I 10830 Å triplet is a very informative indicator of chromospheric activities as the helium is the second most abundant element in the solar atmosphere. Taking advantage of the high resolution of the 1.6 m Goode Solar Telescope (GST) at Big Bear Solar Observatory (BBSO), previous observations have shown clear evidence of the enhanced absorption, instead of typically-observed emission, for two M-class flares. In this study, we analyze the evolution of the He I 10830 10830 Å emission in numerical models and compare it with observations. The models represent the RADYN simulation results obtained from the F-CHROMA database. We consider the models with the injected electron spectra parameters close to observational estimates for the 2013-August-17 flare event ($δ=8$, $E_c = \{15,20\}$ keV, $F=\{1\times 10^{11}, 3\times{}10^{11}\}$ erg cm$^{-2}$) in detail, as well as other available models. The modeling results agree well with observations, in the sense of both the maximum intensity decrease (-17.1%, compared to the observed value of -13.7%) and the trend of temporal variation (initial absorption phase followed by the emission). All models demonstrate the increased number densities and decreased ratio of the upper and lower level populations of He I 10830 10830 Å transition in the initial phase, which enhances the opacity and forms an absorption feature. Models suggest that the temperatures and free electron densities at heights of 1.3-1.5 Mm should be larger than $\sim 10^4$ K and $6\times 10^{11}$ cm$^{-3}$ thresholds for the line to start being in emission.
△ Less
Submitted 17 June, 2020;
originally announced June 2020.
-
Response of SDO/HMI observables to heating of the solar atmosphere by precipitating high-energy electrons
Authors:
Viacheslav M Sadykov,
Alexander G Kosovichev,
Irina N Kitiashvili,
Graham S Kerr
Abstract:
We perform an analysis of the line-of-sight (LOS) observables of the SDO/HMI for models of the solar atmosphere heated by precipitating high-energy electrons during solar flares. The radiative hydrodynamic (RADYN) flare models are obtained from the F-CHROMA database. The Stokes profiles for the Fe 6173 A line observed by SDO/HMI are calculated using the radiative transfer code RH1.5D assuming stat…
▽ More
We perform an analysis of the line-of-sight (LOS) observables of the SDO/HMI for models of the solar atmosphere heated by precipitating high-energy electrons during solar flares. The radiative hydrodynamic (RADYN) flare models are obtained from the F-CHROMA database. The Stokes profiles for the Fe 6173 A line observed by SDO/HMI are calculated using the radiative transfer code RH1.5D assuming statistical equilibrium for atomic level populations and imposing uniform background vertical magnetic field of various strength. The SDO/HMI observing sequence and LOS data processing pipeline algorithm are applied to derive the observables (continuum intensity, line depth, Doppler velocity, LOS magnetic field). Our results reveal that the strongest deviations of the observables from the actual spectroscopic line parameters are found for the model with a total energy deposited of $E_{total}=1.0\times{}10^{12}$ erg cm$^{-2}$, injected with power-law spectral index of $δ=3$ above a low-energy cutoff of $E_{c}=25$ keV. The magnitudes of the velocity and magnetic field deviations depend on the imposed magnetic field, and can reach 0.35 km/s for LOS velocities, 90 G for LOS magnetic field, and 3\% for continuum enhancement for the 1000 G imposed LOS magnetic field setup. For $E_{total}\geq{}3.0\times{}10^{11}$ erg cm$^{-2}$ models, the velocity and magnetic field deviations are most strongly correlated with the energy flux carried by $\sim$50 keV electrons, and the continuum enhancement is correlated with the synthesized $\sim$55-60 keV hard X-ray photon flux. The relatively low magnitudes of perturbations of the observables and absence of magnetic field sign reversals suggest that the considered radiative hydrodynamic beam heating models augmented with the uniform vertical magnetic field setups cannot explain strong transient changes found in the SDO/HMI observations.
△ Less
Submitted 1 March, 2020; v1 submitted 25 June, 2019;
originally announced June 2019.
-
Statistical Properties of Soft X-ray emission of solar flares
Authors:
Viacheslav M Sadykov,
Alexander G Kosovichev,
Irina N Kitiashvili,
Alexander Frolov
Abstract:
We present a statistical analysis of properties of Soft X-Ray (SXR) emission, plasma temperature (T), and emission measure (EM), derived from GOES observations of flares in 2002-2017. The temperature and emission measures are obtained using the TEBBS algorithm (Ryan et al. 2012), which delivers reliable results together with uncertainties even for weak B-class flare events. More than 96% of flares…
▽ More
We present a statistical analysis of properties of Soft X-Ray (SXR) emission, plasma temperature (T), and emission measure (EM), derived from GOES observations of flares in 2002-2017. The temperature and emission measures are obtained using the TEBBS algorithm (Ryan et al. 2012), which delivers reliable results together with uncertainties even for weak B-class flare events. More than 96% of flares demonstrate a sequential appearance of T, SXR, and EM maxima, in agreement with the expected behavior of the chromospheric evaporation process. The relative number of such flares increases with increasing the SXR flux maximum. The SXR maximum is closer in time to the T maximum for B-class flares than for >=C-class flares, while it is very close to the EM maximum for M- and X-class flares. We define flares as "T-controlled" if the time interval between the SXR and T maxima is at least two times shorter than the interval between the EM and SXR maxima, and as "EM-controlled" if the time interval between the EM and SXR maxima is at least two times shorter than the interval between the SXR and T maxima. For any considered flare class range, the T-controlled events compared to EM-controlled events have: a) higher EM but lower T; b) longer durations and shorter relative growth times; c) longer FWHM and characteristic decay times. Interpretation of these statistical results based on analysis of a single loop dynamics suggests that for flares of the same class range, the T-controlled events can be developed in longer loops than the EM-controlled events.
△ Less
Submitted 12 February, 2019; v1 submitted 12 October, 2018;
originally announced October 2018.
-
Statistical Study of Chromospheric Evaporation in Impulsive Phase of Solar Flares
Authors:
Viacheslav M Sadykov,
Alexander G Kosovichev,
Ivan N Sharykin,
Graham S Kerr
Abstract:
We present a statistical study of chromospheric evaporation in solar flares using simultaneous observations by the RHESSI X-ray telescope and the IRIS UV spectrograph. The results are compared with radiation hydrodynamic flare models from the F-CHROMA RADYN database. For each event, we study spatially-resolved Doppler shifts of spectral lines formed in the transition region (C\,II\,1334.5\,Å) and…
▽ More
We present a statistical study of chromospheric evaporation in solar flares using simultaneous observations by the RHESSI X-ray telescope and the IRIS UV spectrograph. The results are compared with radiation hydrodynamic flare models from the F-CHROMA RADYN database. For each event, we study spatially-resolved Doppler shifts of spectral lines formed in the transition region (C\,II\,1334.5\,Å) and hot coronal plasma (Fe\,XXI\,1354.1\,Å) to investigate the dynamics of the solar atmosphere during the flare impulsive phase. We estimate the energy fluxes deposited by high-energy electrons using X-ray imaging spectroscopy and assuming the standard thick-target model. Using the RADYN flare models, the RH 1.5D radiative transfer code and the Chianti atomic line database, we calculate C\,II and Fe\,XXI line profiles and compare with the observations. While the RADYN models predict a correlation between the Doppler shifts and deposited energy flux for both lines, this was only observed in the C II data. Several quantitative discrepancies are found between the observations and models: the Fe\,XXI Doppler shifts are substantially stronger in the models than in the data; the C\,II mean blueshifts are absent in the observations but predicted by the models. The transition energy between `gentle' and `explosive' evaporation regimes estimated from the observations ($2-8\times{}10^{9}$~erg~cm$^{-2}$~s$^{-1}$) and derived from the models ($2.2-10.1\times{}10^{9}$~erg~cm$^{-2}$~s$^{-1}$) are comparable with each other. The results illustrate relationships among the processes of chromospheric evaporation, response of the colder layers, and the flare energy flux deposited by high-energy electrons, although demonstrating discrepancy between analyzed observations and RADYN models.
△ Less
Submitted 4 December, 2018; v1 submitted 27 May, 2018;
originally announced May 2018.
-
Flare Energy Release in the Lower Solar Atmosphere Near the Magnetic Field Polarity Inversion Line
Authors:
I. N. Sharykin,
V. M. Sadykov,
A. G. Kosovichev,
S. Vargas-Dominguez,
I. V. Zimovets
Abstract:
We study flare processes in the solar atmosphere using observational data for a M1-class flare of June 12, 2014, obtained by New Solar Telescope (NST/BBSO) and Helioseismic Magnetic Imager (HMI/SDO). The main goal is to understand triggers and manifestations of the flare energy release in the photosphere and chromosphere using high-resolution optical observations and magnetic field measurements. W…
▽ More
We study flare processes in the solar atmosphere using observational data for a M1-class flare of June 12, 2014, obtained by New Solar Telescope (NST/BBSO) and Helioseismic Magnetic Imager (HMI/SDO). The main goal is to understand triggers and manifestations of the flare energy release in the photosphere and chromosphere using high-resolution optical observations and magnetic field measurements. We analyze optical images, HMI Dopplergrams and vector magnetograms, and use Non-Linear Force-Free Field (NLFFF) extrapolations for reconstruction of the magnetic topology and electric currents. The NLFFF modelling reveals interaction of two magnetic flux ropes with oppositely directed magnetic field in the PIL. These flux ropes are observed as a compact sheared arcade along the PIL in the high-resolution broad-band continuum images from NST. In the vicinity of PIL, the NST H alpha observations reveal formation of a thin three-ribbon structure corresponding to a small-scale photospheric magnetic arcade. The observational results evidence in favor of location of the primary energy release site in the chromospheric plasma with strong electric currents concentrated near the polarity inversion line. In this case, magnetic reconnection is triggered by the interacting magnetic flux ropes forming a current sheet elongated along the PIL.
△ Less
Submitted 13 January, 2018;
originally announced January 2018.
-
Relationships Between Characteristics of the Line-of-Sight Magnetic Field and Solar Flare Forecasts
Authors:
Viacheslav M Sadykov,
Alexander G Kosovichev
Abstract:
We analyze the relationship between the flare X-ray peak flux, and characteristics of the Polarity Inversion Line (PIL) and Active Regions (AR), derived from line-of-sight (LOS) magnetograms. The PIL detection algorithm based on a magnetogram segmentation procedure is applied for each AR with 1 hour cadence. The PIL and AR characteristics are associated with the AR flare history and divided into f…
▽ More
We analyze the relationship between the flare X-ray peak flux, and characteristics of the Polarity Inversion Line (PIL) and Active Regions (AR), derived from line-of-sight (LOS) magnetograms. The PIL detection algorithm based on a magnetogram segmentation procedure is applied for each AR with 1 hour cadence. The PIL and AR characteristics are associated with the AR flare history and divided into flaring and non-flaring cases. Effectiveness of the derived characteristics for flare forecasting is determined by the number of non-flaring cases separated from flaring cases by a certain threshold, and by their Fisher ranking score. The Support Vector Machine (SVM) classifier trained only on the PIL characteristics is used for the flare prediction. We have obtained the following results: (1) the PIL characteristics are more effective than global characteristics of ARs, (2) the highest True Skill Statistics (TSS) values of 0.76$\pm$0.03 for $\geq$M1.0 flares and 0.84$\pm$0.07 for $\geq$X1.0 flares are obtained using the "Sigmoid" SVM kernel, (3) the TSS scores obtained using only the LOS magnetograms are slightly lower than the scores obtained using vector magnetograms, but significantly better than current expert-based predictions, (4) for prediction of $\geq$M1.0 class flares 74.4 % of all cases, and 91.2 % for $\geq$X1.0 class, can be pre-classified as negative with no significant effect on the results, (5) the inclusion of global AR characteristics does not improve the forecast. The study confirms the unique role of the PIL region characteristics in the flare initiation process, and demonstrate possibilities of flare forecasting using only the line-of-sight magnetograms.
△ Less
Submitted 27 May, 2018; v1 submitted 11 April, 2017;
originally announced April 2017.
-
Investigation of Relationship Between High-Energy X-ray Sources and Photospheric and Helioseismic Impacts of X1.8 Solar Flare of October 23, 2012
Authors:
Ivan N Sharykin,
Alexander G Kosovichev,
Viacheslav M Sadykov,
Ivan V Zimovets,
Ivan I Myshyakov
Abstract:
The X-class solar flare of October 23, 2012, generated continuum photospheric emission and a strong helioseismic wave ("sunquake") that points to an intensive energy release in the dense part of the solar atmosphere. We study properties of the energy release with high temporal and spatial resolutions, using photospheric data from the Helioseismic Magnetic Imager (HMI) onboard Solar Dynamics Observ…
▽ More
The X-class solar flare of October 23, 2012, generated continuum photospheric emission and a strong helioseismic wave ("sunquake") that points to an intensive energy release in the dense part of the solar atmosphere. We study properties of the energy release with high temporal and spatial resolutions, using photospheric data from the Helioseismic Magnetic Imager (HMI) onboard Solar Dynamics Observatory (SDO), and hard X-ray observations made by the Ramaty High-Energy Solar Spectroscopic Imager (RHESSI). For this analysis we use level-1 HMI data (filtergrams), obtained by scanning the Fe I line (6731~Å) with the time cadence of $\sim 3.6$ s and spatial resolution of $\sim 0.5^{\prime\prime}$ per pixel. It is found that the photospheric disturbances caused by the flare spatially coincide with the region of hard X-ray emission, but are delayed by $\lesssim 4$ seconds. This delay is consistent with predictions of the flare hydrodynamics RADYN models. However, the models fail to explain the magnitude of variations observed by the HMI. The data indicate that the photospheric impact and helioseismic wave might be caused by the electron energy flux substantially higher than that in the current flare radiative hydrodynamic models.
△ Less
Submitted 8 June, 2017; v1 submitted 10 March, 2017;
originally announced March 2017.
-
Interactive Multi-Instrument Database of Solar Flares
Authors:
Viacheslav M Sadykov,
Alexander G Kosovichev,
Vincent Oria,
Gelu M Nita
Abstract:
Solar flares are complicated physical phenomena that are observable in a broad range of the electromagnetic spectrum, from radiowaves to $γ$-rays. For a more comprehensive understanding of flares, it is necessary to perform a combined multi-wavelength analysis using observations from many satellites and ground-based observatories. For efficient data search, integration of different flare lists and…
▽ More
Solar flares are complicated physical phenomena that are observable in a broad range of the electromagnetic spectrum, from radiowaves to $γ$-rays. For a more comprehensive understanding of flares, it is necessary to perform a combined multi-wavelength analysis using observations from many satellites and ground-based observatories. For efficient data search, integration of different flare lists and representation of observational data, we have developed an Interactive Multi-Instrument Database of Solar Flares (https://solarflare.njit.edu/). The web accessible database is fully functional and allows the user to search for uniquely-identified flare events based on their physical descriptors and availability of observations by a particular set of instruments. Currently, the data from three primary flare lists (GOES, RHESSI and HEK) and a variety of other event catalogs (Hinode, Fermi GBM, Konus-Wind, OVSA flare catalogs, CACTus CME catalog, Filament eruption catalog) and observing logs (IRIS and Nobeyama coverage) are integrated, and an additional set of physical descriptors (temperature and emission measure) is provided along with an observing summary, data links, and multi-wavelength light curves for each flare event since January, 2002. We envision that this new tool will allow researchers to significantly speed up the search of events of interest for statistical and case studies.
△ Less
Submitted 13 June, 2017; v1 submitted 9 February, 2017;
originally announced February 2017.
-
Observational Investigation of Energy Release in the Lower Solar Atmosphere of a Solar Flare
Authors:
I. N. Sharykin,
V. M. Sadykov,
A. G. Kosovichev,
S. Vargas-Dominguez,
I. V. Zimovets
Abstract:
We study flare processes in the lower solar atmosphere using observational data for a M1-class flare of June 12, 2014, obtained by New Solar Telescope (NST/BBSO) and Helioseismic Magnetic Imager (HMI/SDO). The main goal is to understand triggers and manifestations of the flare energy release in the lower layers of the solar atmosphere (the photosphere and chromosphere) using high-resolution optica…
▽ More
We study flare processes in the lower solar atmosphere using observational data for a M1-class flare of June 12, 2014, obtained by New Solar Telescope (NST/BBSO) and Helioseismic Magnetic Imager (HMI/SDO). The main goal is to understand triggers and manifestations of the flare energy release in the lower layers of the solar atmosphere (the photosphere and chromosphere) using high-resolution optical observations and magnetic field measurements. We analyze optical images, HMI Dopplergrams and vector magnetograms, and use Non-Linear Force-Free Field (NLFFF) extrapolations for reconstruction of the magnetic topology. The NLFFF modelling reveals interaction of oppositely directed magnetic flux-tubes in the PIL. These two interacting magnetic flux tubes are observed as a compact sheared arcade along the PIL in the high-resolution broad-band continuum images from NST. In the vicinity of the PIL, the NST H alpha observations reveal formation of a thin three-ribbon structure corresponding to the small-scale photospheric magnetic arcade. Presented observational results evidence in favor of location of the primary energy release site in the dense chromosphere where plasma is partially ionized in the region of strong electric currents concentrated near the polarity inversion line. Magnetic reconnection may be triggered by two interacting magnetic flux tubes with forming current sheet elongated along the PIL.
△ Less
Submitted 18 April, 2016;
originally announced April 2016.
-
Relationship between chromospheric evaporation and magnetic field topology in M-class solar flare
Authors:
V. M. Sadykov,
A. G. Kosovichev,
I. N. Sharykin,
I. V. Zimovets,
S. Vargas Dominguez
Abstract:
Chromospheric evaporation is observed as Doppler blueshift during solar flares. It plays one of key roles in dynamics and energetics of solar flares, however, its mechanism is still unknown. In this paper we present a detailed analysis of spatially-resolved multi-wavelength observations of chromospheric evaporation during an M 1.0 class solar flare (SOL2014-06-12T21:12) using data from the NASA's…
▽ More
Chromospheric evaporation is observed as Doppler blueshift during solar flares. It plays one of key roles in dynamics and energetics of solar flares, however, its mechanism is still unknown. In this paper we present a detailed analysis of spatially-resolved multi-wavelength observations of chromospheric evaporation during an M 1.0 class solar flare (SOL2014-06-12T21:12) using data from the NASA's IRIS (Interface Region Imaging Spectrograph) and HMI/SDO (Helioseismic and Magnetic Imager onboard Solar Dynamics Observatory) telescopes, and VIS/NST (Visible Imaging Spectrometer at New Solar Telescope) high-resolution observations, covering the temperature range from 10^4 K to 10^7 K. The results show that the averaged over the region Fe XXI blueshift of the hot evaporating plasma is delayed relative to the C II redshift of the relatively cold chromospheric plasma by about 1 min. The spatial distribution of the delays is not uniform across the region and can be as long as 2 min in several zones. Using vector magnetograms from HMI we reconstruct the magnetic field topology and the quasi-separatrix layer (QSL) and find that the blueshift delay regions as well as the H-alpha flare ribbons are connected to the region of magnetic polarity inversion line (PIL) and an expanding flux rope via a system of low-lying loop arcades with height < ~4.5 Mm. This allows us to propose an interpretation of the chromospheric evaporation based on the geometry of local magnetic fields, and the primary energy source associated with the PIL.
△ Less
Submitted 18 April, 2016;
originally announced April 2016.
-
Using MHD simulations to model H-alpha and UV spectral lines for interpretation of IRIS and NST data
Authors:
Viacheslav M. Sadykov,
Alexander G. Kosovichev
Abstract:
We present results of non-LTE modeling of H-alpha 6563 A and Mg II k&h 2796 A and 2803 A lines. This modeling is important for interpretation of coordinated observations from the recently launched NASA's IRIS mission and from the New Solar Telescope at Big Bear Solar Observatory. Among available codes for the non-LTE modeling, the RH code is chosen as the most appropriate for modeling of the line…
▽ More
We present results of non-LTE modeling of H-alpha 6563 A and Mg II k&h 2796 A and 2803 A lines. This modeling is important for interpretation of coordinated observations from the recently launched NASA's IRIS mission and from the New Solar Telescope at Big Bear Solar Observatory. Among available codes for the non-LTE modeling, the RH code is chosen as the most appropriate for modeling of the line profiles. The most suitable Hydrogen and Magnesium atomic models are selected by performing several tests of the code. The influence of the ionization degree on the line profiles is also studied. Radiative-MHD simulations of the solar atmosphere, obtained with the Bifrost code, are used as input data for calculation of synthetic spectra of the H-alpha and Mg II h&k lines for particular locations evolving with time. The spectral line variations reveal the presence of MHD waves in the simulation results. We construct oscillation power spectra of the line intensity for different wavelength, and compare these with the corresponding height-dependent power spectra of atmospheric parameters from the simulations. We find correlations between the power spectra of intensities of the line profiles at certain wavelengths and the power spectra of the atmospheric parameters at the tau-unity heights for these wavelengths. These results provide a new diagnostic method of chromospheric oscillations; however, larger amounts of data are needed to confirm these correlations.
△ Less
Submitted 30 November, 2014;
originally announced December 2014.
-
Properties of chromospheric evaporation and plasma dynamics of a solar flare from IRIS observations
Authors:
Viacheslav M. Sadykov,
Santiago Vargas Dominguez,
Alexander G. Kosovichev,
Ivan N. Sharykin,
Alexei B. Struminsky,
Ivan V. Zimovets
Abstract:
Dynamics of hot chromospheric plasma of solar flares is a key to understanding of mechanisms of flare energy release and particle acceleration. A moderate M1.0 class flare of 12 June, 2014 (SOL2014-06-12T21:12) was simultaneously observed by NASA's Interface Region Imaging Spectrograph (IRIS), other spacecraft, and also by New Solar Telescope (NST) at the BBSO. This paper presents the first part o…
▽ More
Dynamics of hot chromospheric plasma of solar flares is a key to understanding of mechanisms of flare energy release and particle acceleration. A moderate M1.0 class flare of 12 June, 2014 (SOL2014-06-12T21:12) was simultaneously observed by NASA's Interface Region Imaging Spectrograph (IRIS), other spacecraft, and also by New Solar Telescope (NST) at the BBSO. This paper presents the first part of our investigation focused on analysis of the IRIS data. Our analysis of the IRIS data in different spectral lines reveals strong redshifted jet-like flow with the speed of ~100 km/s of the chromospheric material before the flare. Strong nonthermal emission of the C II k 1334.5 A line, formed in the chromosphere-corona transition region, is observed at the beginning of the impulsive phase in several small (with a size of ~1 arcsec) points. It is also found that the C II k line is redshifted across the flaring region before, during and after the impulsive phase. A peak of integrated emission of the hot (1.1 MK) plasma in the Fe XXI 1354.1 A line is detected approximately 5 minutes after the integrated emission peak of the lower temperature C II k. A strong blueshift of the Fe XXI line across the flaring region corresponds to evaporation flows of the hot chromospheric plasma with a speed of 50 km/s. Additional analysis of the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) data supports the idea that the upper chromospheric dynamics observed by IRIS has features of "gentle" evaporation driven by heating of the solar chromosphere by accelerated electrons and by a heat flux from the flare energy release site.
△ Less
Submitted 9 April, 2015; v1 submitted 29 November, 2014;
originally announced December 2014.
