Search

The EU funded the FLARECAST project, that ran from Jan 2015 until Feb 2018. FLARECAST had a R2O focus, and introduced several innovations into the discipline of solar flare forecasting. FLARECAST innovations were: first, the treatment of hundreds of physical properties viewed as promising flare predictors on equal footing, extending multiple previous works; second, the use of fourteen (14) different ML techniques, also on equal footing, to optimize the immense Big Data parameter space created by these many predictors; third, the establishment of a robust, three-pronged communication effort oriented toward policy makers, space-weather stakeholders and the wider public. FLARECAST pledged to make all its data, codes and infrastructure openly available worldwide. The combined use of 170+ properties (a total of 209 predictors are now available) in multiple ML algorithms, some of which were designed exclusively for the project, gave rise to changing sets of best-performing predictors for the forecasting of different flaring levels. At the same time, FLARECAST reaffirmed the importance of rigorous training and testing practices to avoid overly optimistic pre-operational prediction performance. In addition, the project has (a) tested new and revisited physically intuitive flare predictors and (b) provided meaningful clues toward the transition from flares to eruptive flares, namely, events associated with coronal mass ejections (CMEs). These leads, along with the FLARECAST data, algorithms and infrastructure, could help facilitate integrated space-weather forecasting efforts that take steps to avoid effort duplication. In spite of being one of the most intensive and systematic flare forecasting efforts to-date, FLARECAST has not managed to convincingly lift the barrier of stochasticity in solar flare occurrence and forecasting: solar flare prediction thus remains inherently probabilistic., Comment: 67 pages, 14 figures; submitted

This paper shows that compressed sensing realized by means of regularized deconvolution and the Finite Isotropic Wavelet Transform is effective and reliable in hard X-ray solar imaging. The method utilizes the Finite Isotropic Wavelet Transform with Meyer function as the mother wavelet. Further, compressed sensing is realized by optimizing a sparsity-promoting regularized objective function by means of the Fast Iterative Shrinkage-Thresholding Algorithm. Eventually, the regularization parameter is selected by means of the Miller criterion. The method is applied against both synthetic data mimicking the Spectrometer/Telescope Imaging X-rays (STIX) measurements and experimental observations provided by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The performances of the method are compared with the results provided by standard visibility-based reconstruction methods. The results show that the application of the sparsity constraint and the use of a continuous, isotropic framework for the wavelet transform provide a notable spatial accuracy and significantly reduce the ringing effects due to the instrument point spread functions.

In this paper we present a hybrid approach to numerically solve two-dimensional electromagnetic inverse scattering problems, whereby the unknown scatterer is hosted by a possibly inhomogeneous background. The approach is `hybrid' in that it merges a qualitative and a quantitative method to optimize the way of exploiting the a priori information on the background within the inversion procedure, thus improving the quality of the reconstruction and reducing the data amount necessary for a satisfactory result. In the qualitative step, this a priori knowledge is utilized to implement the linear sampling method in its near-field formulation for an inhomogeneous background, in order to identify the region where the scatterer is located. On the other hand, the same a priori information is also encoded in the quantitative step by extending and applying the contrast source inversion method to what we call the `inhomogeneous Lippmann-Schwinger equation': the latter is a generalization of the classical Lippmann-Schwinger equation to the case of an inhomogeneous background, and in our paper is deduced from the differential formulation of the direct scattering problem to provide the reconstruction algorithm with an appropriate theoretical basis. Then, the point values of the refractive index are computed only in the region identified by the linear sampling method at the previous step. The effectiveness of this hybrid approach is supported by numerical simulations presented at the end of the paper., Comment: accepted in SIAM Journal on Applied Mathematics

X-radiation from energetic electrons is the prime diagnostic of flare-accelerated electrons. The observed X-ray flux (and polarization state) is fundamentally a convolution of the cross-section for the hard X-ray emission process(es) in question with the electron distribution function, which is in turn a function of energy, direction, spatial location and time. To address the problems of particle propagation and acceleration one needs to infer as much information as possible on this electron distribution function, through a deconvolution of this fundamental relationship. This review presents recent progress toward this goal using spectroscopic, imaging and polarization measurements, primarily from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). Previous conclusions regarding the energy, angular (pitch angle) and spatial distributions of energetic electrons in solar flares are critically reviewed. We discuss the role and the observational evidence of several radiation processes: free-free electron-ion, free-free electron-electron, free-bound electron-ion bremsstrahlung, photoelectric absorption and Compton back-scatter (albedo), using both spectroscopic and imaging techniques. This unprecedented quality of data allows for the first time inference of the angular distributions of the X-ray-emitting electrons using albedo, improved model-independent inference of electron energy spectra and emission measures of thermal plasma. Moreover, imaging spectroscopy has revealed hitherto unknown details of solar flare morphology and detailed spectroscopy of coronal, footpoint and extended sources in flaring regions. Additional attempts to measure hard X-ray polarization were not sufficient to put constraints on the degree of anisotropy of electrons, but point to the importance of obtaining good quality polarization data., Comment: This is an article for a monograph on the physics of solar flares, inspired by RHESSI observations. The individual articles are to appear in Space Science Reviews (2011); 10.1007/s11214-011-9804-x

Hard X-Ray bremsstrahlung continuum spectra, such as from solar flares, are commonly described in terms of power-law fits, either to the photon spectra themselves or to the electron spectra responsible for them. In applications various approximate relations between electron and photon spectral indices are often used for energies both above and below electron low-energy cutoffs. We examine the form of the exact relationships in various situations, and for various cross-sections, showing that empirical relations sometimes used can be highly misleading and consider how to improve fitting procedures. We obtain expressions for photon spectra from single, double and truncated power-law electron spectra for a variety of cross-sections and for the thin and thick target models and simple analytic expressions for the Bethe-Heitler cases. We show that above a low-energy cutoff the Kramers and Bethe-Heitler results match reasonably well with results for exact cross-sections up to energies around 100 keV; that below the low-energy cutoff, Kramers and other constant spectral index forms commonly used are very poor approximations to accurate results; but that our analytical forms are a very good match. Analytical forms of the Bethe-Heitler photon spectra from general power-law electron spectra are an excellent match to exact results for both thin and thick targets and they enable much faster spectral fitting than evaluation of the full spectral integrations., Comment: 7 pages, 4 figures, submitted to Astronomy and Astrophysics

Hard X-ray spectra in solar flares provide knowledge of the electron spectrum that results from acceleration and propagation in the solar atmosphere. However, the inference of the electron spectra from solar X-ray spectra is an ill-posed inverse problem. Here we develop and apply an enhanced regularization algorithm for this process making use of physical constraints on the form of the electron spectrum. The algorithm incorporates various features not heretofore employed in the solar flare context: Generalized Singular Value Decomposition (GSVD) to deal with different orders of constraints; rectangular form of the cross-section matrix to extend the solution energy range; regularization with various forms of the smoothing operator; and "preconditioning" of the problem. We show by simulations that this technique yields electron spectra with considerably more information and higher quality than previous algorithms., Comment: 21 pages, 8 fugures, accepted to Solar Physics

Kontar et al. (2004) have shown how to recover mean source electron spectra in solar flares through a physical constraint regularization analysis of the bremsstrahlung photon spectra that they produce. They emphasize the use of non-square inversion techniques, and preconditioning combined with physical properties of the spectra to achieve the most meaningful solution to the problem. Higher-order regularization techniques may be used to generate ${overline F}(E)$ forms with certain desirable properties (e.g., higher order derivatives). They further note that such analyses may be used to infer properties of the electron energy spectra at energies well above the maximum photon energy observed. In this paper we apply these techniques to data from a solar flare observed by RHESSI on 26 February, 2002. Results using different orders of regularization are presented and compared for various time intervals. Clear evidence is presented for a change in the value of the high-energy cutoff in the mean source electron spectrum with time. We also show how the construction of the injected (accelerated) electron spectrum $F_0(E_0)$ (assuming that Coulomb collisions in a cold target dominate the electron energetics) is facilitated by the use of higher-order regularization methods., Comment: 10 pages, 7 figures, accepted to Solar Physics