We have compiled a catalog of solar flares as observed by the Extreme ultraviolet Imaging Telescope (EIT) aboard the Solar and Heliospheric Observatory (SOHO) spacecraft and the GOES spacecraft over a span from 1997 to 2010. During mid-1998, the cadence of EIT images was revised from two images per day to 12 minutes. However, the low temporal resolution causes significant data gaps in capturing much of the flaring phenomenon. Therefore, we monitor possible errors in flare detection by flare parameters such as temporal overlap, observational wavelength, and considering full field of view (FOV) images. We consider the GOES flare catalog as the primary source. We describe the technique used to enhance the GOES detected flares using the Extreme Ultraviolet (EUV) image captured by the EIT instrument. In order to detect brightenings, we subtract the images with a maximum cadence of 25 minutes. We have downloaded and analyzed the EIT data via the Virtual Solar Observatory (VSO). This flare dataset from the SOHO/EIT period proves indispensable to the process of the solar flare predictions as the instrument has covered most of Solar Cycle 23.

The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.

The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.

It is well known that solar filaments are features in the solar atmosphere which show a hemispheric preference in their chirality. The hemispheric preference is such that the dextral chirality dominates in the northern hemisphere while the sinistral chirality dominates in the southern. Determining the strength and cyclic variation of the degree of this hemispheric preference, however, is challenging and tedious and thus needs to be automated. In this paper, we follow Dr. Pietro Bernasconi's algorithm (Bernasconi et al. 2005) to detect filament chirality with two parallel channels of application. The algorithm is applied to H-alpha images by "Advanced Automated Filament Detection and Characterization Code" (AAFDCC) (Bernasconi et al. 2005) or the algorithm is explained to a human and the human determines the chirality of the solar filament. We have conducted this exercise on data during the month of August from years 2000 to 2016 and we found that 83 % of our visually determined filaments follow the hemispheric chirality preference, while 58 % of automatically determined filaments follow it. Our visually compiled results have over 90 % agreement with those of Pevtsov et al. (2003), yet the visually determined chiralities of filaments disagree with automated determinations significantly. We find that the hemispheric preference remained the same between solar cycles 23 and 24 but the preference is very difficult to determine during the solar minimum of 2008-2010 primarily due to the absence of filaments.

In 1844 Schwabe discovered that the number of sunspots increased and decreased over a period of about 11 years, that variation became known as the sunspot cycle. Almost eighty years later, Hale described the nature of the Sun's magnetic field, identifying that it takes about 22 years for the Sun's magnetic polarity to cycle. It was also identified that the latitudinal distribution of sunspots resembles the wings of a butterfly showing migration of sunspots in each hemisphere that abruptly start at mid-latitudes towards the Sun's equator over the next 11 years. These sunspot patterns were shown to be asymmetric across the equator. In intervening years, it was deduced that the Sun (and sun-like stars) possess magnetic activity cycles that are assumed to be the physical manifestation of a dynamo process that results from complex circulatory transport processes in the star's interior. Understanding the Sun's magnetism, its origin and its variation, has become a fundamental scientific objective \-- the distribution of magnetism, and its interaction with convective processes, drives various plasma processes in the outer atmosphere. In the past few decades, a range of diagnostic techniques have been employed to systematically study finer scale magnetized objects, and associated phenomena. The patterns discerned became known as the ``Extended Solar Cycle'' (ESC). The patterns of the ESC appeared to extend the wings of the activity butterfly back in time, nearly a decade before the formation of the sunspot pattern, and to much higher solar latitudes. In this short review, we describe their observational patterns of the ESC and discuss possible connections to the solar dynamo as we depart on a multi-national collaboration to investigate the origins of solar magnetism through a blend of archived and contemporary data analysis with the goal of improving solar dynamo understanding and modeling.

Solar energetic particles are a result of intense solar events such as solar flares and Coronal Mass Ejections (CMEs). These latter events all together can cause major disruptions to spacecraft that are in Earth's orbit and outside of the magnetosphere. In this work we are interested in establishing the necessary conditions for a major geo-effective solar particle storm immediately after a major flare, namely the existence of a direct magnetic connection. To our knowledge, this is the first work that explores not only the correlations of GOES X-ray and proton channels, but also the correlations that happen across all the proton channels. We found that proton channels auto-correlations and cross-correlations may also be precursors to the occurrence of an SEP event. In this paper, we tackle the problem of predicting >100 MeV SEP events from a multivariate time series perspective using easily interpretable decision tree models.

The prediction of solar flares, eruptions, and high energy particle storms is of great societal importance. The data mining approach to forecasting has been shown to be very promising. Benchmark datasets are a key element in the further development of data-driven forecasting. With one or more benchmark data sets established, judicious use of both the data themselves and the selection of prediction algorithms is key to developing a high quality and robust method for the prediction of geo-effective solar activity. We review here briefly the process of generating benchmark datasets and developing prediction algorithms.

The purpose of this paper is to explore a resolution for the Faint Young Sun Paradox that has been mostly rejected by the community, namely the possibility of a somewhat more massive young Sun with a large mass loss rate sustained for two to three billion years. This would make the young Sun bright enough to keep both the terrestrial and Martian oceans from freezing, and thus resolve the paradox. It is found that a large and sustained mass loss is consistent with the well observed spin-down rate of Sun-like stars, and indeed may be required for it. It is concluded that a more massive young Sun must be considered a plausible hypothesis.

In this work we introduce a new way of binning sunspot group data with the purpose of better understanding the impact of the solar cycle on sunspot properties and how this defined the characteristics of the extended minimum of cycle 23. Our approach assumes that the statistical properties of sunspots are completely determined by the strength of the underlying large-scale field and have no additional time dependencies. We use the amplitude of the cycle at any given moment (something we refer to as activity level) as a proxy for the strength of this deep-seated magnetic field. We find that the sunspot size distribution is composed of two populations: one population of groups and active regions and a second population of pores and ephemeral regions. When fits are performed at periods of different activity level, only the statistical properties of the former population, the active regions, is found to vary. Finally, we study the relative contribution of each component (small-scale versus large-scale) to solar magnetism. We find that when hemispheres are treated separately, almost every one of the past 12 solar minima reaches a point where the main contribution to magnetism comes from the small-scale component. However, due to asymmetries in cycle phase, this state is very rarely reached by both hemispheres at the same time. From this we infer that even though each hemisphere did reach the magnetic baseline, from a heliospheric point of view the minimum of cycle 23 was not as deep as it could have been.

In this work we take advantage of eleven different sunspot group, sunspot, and active region databases to characterize the area and flux distributions of photospheric magnetic structures. We find that, when taken separately, different databases are better fitted by different distributions (as has been reported previously in the literature). However, we find that all our databases can be reconciled by the simple application of a proportionality constant, and that, in reality, different databases are sampling different parts of a composite distribution. This composite distribution is made up by linear combination of Weibull and log-normal distributions -- where a pure Weibull (log-normal) characterizes the distribution of structures with fluxes below (above) $10^{21}$Mx ($10^{22}$Mx). We propose that this is evidence of two separate mechanisms giving rise to visible structures on the photosphere: one directly connected to the global component of the dynamo (and the generation of bipolar active regions), and the other with the small-scale component of the dynamo (and the fragmentation of magnetic structures due to their interaction with turbulent convection). Additionally, we demonstrate that the Weibull distribution shows the expected linear behavior of a power-law distribution (when extended into smaller fluxes), making our results compatible with the results of Parnell et al. (2009).

Direct observations over the past four centuries show that the number of sunspots observed on the Sun's surface vary periodically, going through successive maxima and minima. Following sunspot cycle 23, the Sun went into a prolonged minimum characterized by a very weak polar magnetic field and an unusually large number of days without sunspots. Sunspots are strongly magnetized regions and are generated by a dynamo mechanism which recreates the solar polar field mediated via plasma flows. Here we report results from kinematic dynamo simulations which demonstrate that a fast meridional flow in the early half of a cycle, followed by a slower flow in the latter half, reproduces both the characteristics of the minimum of sunspot cycle 23 - a large number of spotless days and a relatively weak polar field. Our model predicts that, in general, very deep minima are associated with weak polar fields. Sunspots govern the solar radiative energy and radio flux, and in conjunction with the polar field, modulate the solar wind, heliospheric open flux and consequently cosmic ray flux at Earth; our results therefore provide a causal link between solar internal dynamics and the atypical values of these heliospheric forcing parameters during the recently concluded solar minimum.

We present a fast method for reconstructing Differential Emission Measures (DEMs) using solar coronal data. On average, the method computes over 1000 DEMs per second for a sample active region observed by the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO), and achieves reduced chi-squared of order unity with no negative emission in all but a few test cases. The high performance of this method is especially relevant in the context of AIA, which images of order one million solar pixels per second. This paper describes the method, analyzes its fidelity, compares its performance and results with other DEM methods, and applies it to an active region and loop observed by AIA and by the Extreme-ultraviolet Imaging Spectrometer (EIS) on Hinode.

Decaying active region 10942 is investigated from 4:00-16:00 UT on February 24, 2007 using a suite of EUV observing instruments. Results from Hinode/EIS, STEREO and TRACE show that although the active region has decayed and no sunspot is present, the physical mechanisms that produce distinguishable loop structures, spectral line broadening, and plasma flows still occur. A coronal loop that appears as a blue-shifted structure in Doppler maps is apparent in intensity images of log(T) = 6.0-6.3 ions. The loop structure is found to be anti-correlated with spectral line broadening generally attributed to nonthermal velocities. This coronal loop structure is investigated physically (temperature, density, geometry) and temporally. Lightcurves created from imaging instruments show brightening and dimming of the loop structure on two different time scales; short pulses of 10-20 min and long duration dimming of 2-4 hours until its disappearance. The coronal loop structure, formed from relatively blue-shifted material that is anti-correlated with spectral line broadening, shows a density of 10^10 to 10^9.3 cm-3 and is visible for longer than characteristic cooling times. The maximum nonthermal spectral line broadenings are found to be adjacent to the footpoint of the coronal loop structure.

The emergence of tilted bipolar active regions and the dispersal of their flux, mediated via processes such as diffusion, differential rotation and meridional circulation is believed to be responsible for the reversal of the Sun's polar field. This process (commonly known as the Babcock-Leighton mechanism) is usually modeled as a near-surface, spatially distributed $\alpha$-effect in kinematic mean-field dynamo models. However, this formulation leads to a relationship between polar field strength and meridional flow speed which is opposite to that suggested by physical insight and predicted by surface flux-transport simulations. With this in mind, we present an improved double-ring algorithm for modeling the Babcock-Leighton mechanism based on active region eruption, within the framework of an axisymmetric dynamo model. Using surface flux-transport simulations we first show that an axisymmetric formulation -- which is usually invoked in kinematic dynamo models -- can reasonably approximate the surface flux dynamics. Finally, we demonstrate that our treatment of the Babcock-Leighton mechanism through double-ring eruption leads to an inverse relationship between polar field strength and meridional flow speed as expected, reconciling the discrepancy between surface flux-transport simulations and kinematic dynamo models.

The turbulent magnetic diffusivity in the solar convection zone is one of the most poorly constrained ingredients of mean-field dynamo models. This lack of constraint has previously led to controversy regarding the most appropriate set of parameters, as different assumptions on the value of turbulent diffusivity lead to radically different solar cycle predictions. Typically, the dynamo community uses double step diffusivity profiles characterized by low values of diffusivity in the bulk of the convection zone. However, these low diffusivity values are not consistent with theoretical estimates based on mixing-length theory -- which suggest much higher values for turbulent diffusivity. To make matters worse, kinematic dynamo simulations cannot yield sustainable magnetic cycles using these theoretical estimates. In this work we show that magnetic cycles become viable if we combine the theoretically estimated diffusivity profile with magnetic quenching of the diffusivity. Furthermore, we find that the main features of this solution can be reproduced by a dynamo simulation using a prescribed (kinematic) diffusivity profile that is based on the spatiotemporal geometric-average of the dynamically quenched diffusivity. Here, we provide an analytic fit to the dynamically quenched diffusivity profile, which can be used in kinematic dynamo simulations. Having successfully reconciled the mixing-length theory estimated diffusivity profile with kinematic dynamo models, we argue that they remain a viable tool for understanding the solar magnetic cycle.

We study the relationship between full-disk solar radiative flux at different wavelengths and average solar photospheric magnetic-flux density, using daily measurements from the Kitt Peak magnetograph and other instruments extending over one or more solar cycles. We use two different statistical methods to determine the underlying nature of these flux-flux relationships. First, we use statistical correlation and regression analysis and show that the relationships are not monotonic for total solar irradiance and for continuum radiation from the photosphere, but are approximately linear for chromospheric and coronal radiation. Second, we use signal theory to examine the flux-flux relationships for a temporal component. We find that a well-defined temporal component exists and accounts for some of the variance in the data. This temporal component arises because active regions with high magnetic field strength evolve, breaking up into small-scale magnetic elements with low field strength, and radiative and magnetic fluxes are sensitive to different active-region components. We generate empirical models that relate radiative flux to magnetic flux, allowing us to predict spectral-irradiance variations from observations of disk-averaged magnetic-flux density. In most cases, the model reconstructions can account for 85-90% of the variability of the radiative flux from the chromosphere and corona. Our results are important for understanding the relationship between magnetic and radiative measures of solar and stellar variability.

The structure of electric current and magnetic helicity in the solar corona is closely linked to solar activity over the 11-year cycle, yet is poorly understood. As an alternative to traditional current-free "potential field" extrapolations, we investigate a model for the global coronal magnetic field which is non-potential and time-dependent, following the build-up and transport of magnetic helicity due to flux emergence and large-scale photospheric motions. This helicity concentrates into twisted magnetic flux ropes, which may lose equilibrium and be ejected. Here, we consider how the magnetic structure predicted by this model-in particular the flux ropes-varies over the solar activity cycle, based on photospheric input data from six periods of cycle 23. The number of flux ropes doubles from minimum to maximum, following the total length of photospheric polarity inversion lines. However, the number of flux rope ejections increases by a factor of eight, following the emergence rate of active regions. This is broadly consistent with the observed cycle modulation of coronal mass ejections, although the actual rate of ejections in the simulation is about a fifth of the rate of observed events. The model predicts that, even at minimum, differential rotation will produce sheared, non-potential, magnetic structure at all latitudes.

The relative importance of different initiation mechanisms for coronal mass ejections (CMEs) on the Sun is uncertain. One possible mechanism is the loss of equilibrium of coronal magnetic flux ropes formed gradually by large-scale surface motions. In this paper, the locations of flux rope ejections in a recently-developed quasi-static global evolution model are compared with observed CME source locations over a 4.5-month period in 1999. Using EUV data, the low-coronal source locations are determined unambiguously for 98 out of 330 CMEs. Despite the incomplete observations, positive correlation (with coefficient up to 0.49) is found between the distributions of observed and simulated ejections, but only when binned into periods of one month or longer. This binning timescale corresponds to the time interval at which magnetogram data are assimilated into the coronal simulations, and the correlation arises primarily from the large-scale surface magnetic field distribution; only a weak dependence is found on the magnetic helicity imparted to the emerging active regions. The simulations are limited in two main ways: they produce fewer ejections, and they do not reproduce the strong clustering of observed CME sources into active regions. Due to this clustering, the horizontal gradient of radial photospheric magnetic field is better correlated with the observed CME source distribution (coefficient 0.67). Our results suggest that, while the gradual formation of magnetic flux ropes over weeks can account for many observed CMEs, especially at higher latitudes, there exists a second class of CMEs (at least half) for which dynamic active region flux emergence on shorter timescales must be the dominant factor.

An essential ingredient in kinematic dynamo models is the velocity field within the solar convection zone. In particular, the differential rotation is now well constrained by helioseismic observations. Helioseismology also gives us information about the depth-dependence of the meridional circulation in the near-surface layers. The typical velocity inputs used in solar dynamo models, however, continue to be an analytic fit to the observed differential rotation and a theoretically constructed meridional flow profile that matches only the peak flow speed at the surface. Here we take the first steps towards realistic helioseismic data assimilation, by presenting methodologies for constructing differential rotation and meridional circulation profiles that more closely conform to the observational constraints currently available. We also present simulations driven by the assimilated rotation and four plausible profiles for the internal meridional circulation -- all of which match the helioseismically inferred near-surface depth-dependence, but whose magnitudes are made to vary. We discuss how the results compare with those that are driven by purely analytic fits. Our results indicate that the latitudinal shear of the rotation in the bulk of the solar convection zone plays a more important role, than either the tachocline or surface radial shear, in the induction of toroidal field. We also find that it is the speed of the equatorward counter-flow in the meridional flow at the base of the convection zone, and not how far into the radiative interior it penetrates, that primarily determines the dynamo cycle period. Given that improved helioseismic constraints are expected to be available in the future, our analysis lays the basis for assimilating these data within dynamo models.

Using the 2004 Venus transit of the Sun to constrain a semi-empirical point-spread function for the TRACE EUV solar telescope, we have measured the effect of stray light in that telescope. We find that 43% of 171A EUV light that enters TRACE is scattered, either through diffraction off the entrance filter grid or through other nonspecular effects. We carry this result forward, via known-PSF deconvolution of TRACE images, to identify its effect on analysis of TRACE data. Known-PSF deconvolution by this derived PSF greatly reduces the effect of visible haze in the TRACE 171A images, enhances bright features, and reveals that the smooth background component of the corona is considerably less bright (and hence much more rarefied) than commonly supposed. Deconvolution reveals that some prior conlclusions about the Sun appear to have been based on stray light in the images. In particular, the diffuse background "quiet corona" becomes consistent with hydrostatic support of the coronal plasma; feature contrast is greatly increased, possibly affecting derived parameters such as the form of the coronal heating function; and essentially all existing differential emission measure studies of small features appear to be affected by contamination from nearby features. We speculate on further implications of stray light for interpretation of EUV images from TRACE and similar instruments, and advocate deconvolution as a standard tool for image analysis with future instruments such as SDO/AIA.

The bulk of solar coronal radiative loss consists of soft X-ray emission from quasi-static loops at the cores of Active Regions. In order to develop diagnostics for determining the heating mechanism of these loops from observations by coronal imaging instruments, I have developed analytical solutions for the temperature structure and scaling laws of loop strands for a wide range of heating functions, including footpoint heating, uniform heating, and heating concentrated at the loop apex. Key results are that the temperature profile depends only weakly on the heating distribution -- not sufficiently to be of significant diagnostic value -- and that the scaling laws survive for this wide range of heating distributions, but with the constant of proportionality in the RTV scaling law ($P_{0}L \thicksim T_{max}^3$) depending on the specific heating function. Furthermore, quasi-static analytical solutions do not exist for an excessive concentration of heating near the loop footpoints, a result in agreement with recent numerical simulations. It is demonstrated that a generalization of the solutions to the case of a strand with a variable diameter leads to only relatively small correction factors in the scaling laws and temperature profiles for constant diameter loop strands. A quintet of leading theoretical coronal heating mechanisms is shown to be captured by the formalism of this paper, and the differences in thermal structure between them may be verified through observations. Preliminary results from full numerical simulations demonstrate that, despite the simplifying assumptions, the analytical solutions from this paper are stable and accurate.

The observations from 1998 April 20 taken with the Coronal Diagnostics Spectrometer CDS on SOHO of a coronal loop on the limb have shown that the plasma was multi-thermal along each line of sight investigated, both before and after background subtraction. The latter result relied on Emission Measure Loci plots, but in this Letter, we used a forward folding technique to produce Differential Emission Measure curves. We also calculate DEM-weighted temperatures for the chosen pixels and find a gradient in temperature along the loop as a function of height that is not compatible with the flat profiles reported by numerous authors for loops observed with EIT on SOHO and TRACE. We also find discrepancies in excess of the mathematical expectation between some of the observed and predicted CDS line intensities. We demonstrate that these differences result from well-known limitations in our knowledge of the atomic data and are to be expected. We further show that the precision of the DEM is limited by the intrinsic width of the ion emissivity functions that are used to calculate the DEM. Hence we conclude that peaks and valleys in the DEM, while in principle not impossible, cannot be confirmed from the data.

Coronal Mass Ejections continuously drag closed magnetic field lines away from the Sun, adding new flux to the interplanetary magnetic field (IMF). We propose that the outward-moving blobs that have been observed in helmet streamers are evidence of ongoing, small-scale reconnection in streamer current sheets, which may play an important role in the prevention of an indefinite buildup of the IMF. Reconnection between two open field lines from both sides of a streamer current sheet creates a new closed field line, which becomes part of the helmet, and a disconnected field line, which moves outward. The blobs are formed by plasma from the streamer that is swept up in the trough of the outward moving field line. We show that this mechanism is supported by observations from SOHO/LASCO. Additionally, we propose a thorough statistical study to quantify the contribution of blob formation to the reduction of the IMF, and indicate how this mechanism may be verified by observations with SOHO/UVCS and the proposed NASA STEREO and ESA Polar Orbiter missions.

A self-consistent mathematical description of the magnetic field of fluted sunspot penumbrae is presented. This description is based on an expansion of the nonlinear force-free magnetohydrostatic equations written in cylindrical coordinates. The lowest order solutions are mathematically equivalent to laminated force-free equilibria in Cartesian geometry. The lowest order solutions have no toroidal component of the magnetic field and the magnetic pressure does not vary with azimuth but the solutions allow arbitrary variations of the magnetic field components with azimuth. Explicit solutions are presented which have a realistic radial profile of the magnetic field strength and reproduce the basic features of the observations.