The Coronal Global Evolutionary Model (CGEM) provides data-driven simulations of the magnetic field in the solar corona to better understand the build-up of magnetic energy that leads to eruptive events. The CGEM project has developed six capabilities. CGEM modules (1) prepare time series of full-disk vector magnetic field observations to (2) derive the changing electric field in the solar photosphere over active-region scales. This local electric field is (3) incorporated into a surface flux transport model that reconstructs a global electric field that evolves magnetic flux in a consistent way. These electric fields drive a (4) 3D spherical magneto-frictional (SMF) model, either at high-resolution over a restricted range of solid angle or at lower resolution over a global domain, to determine the magnetic field and current density in the low corona. An SMF-generated initial field above an active region and the evolving electric field at the photosphere are used to drive (5) detailed magneto-hydrodynamic (MHD) simulations of active regions in the low corona. SMF or MHD solutions are then used to compute emissivity proxies that can be compared with coronal observations. Finally, a lower-resolution SMF magnetic field is used to initialize (6) a global MHD model that is driven by an SMF electric-field time series to simulate the outer corona and heliosphere, ultimately connecting Sun to Earth. As a demonstration, this report features results of CGEM applied to observations of the evolution of NOAA Active Region 11158 in February 2011.

We describe the PDFI_SS software library, which is designed to find the electric field at the Sun's photosphere from a sequence of vector magnetogram and Doppler velocity measurements, and estimates of horizontal velocities obtained from local correlation tracking using the recently upgraded FLCT code. The library, a collection of Fortran subroutines, uses the "PDFI" technique described by Kazachenko et al. (2014), but modified for use in spherical, Plate-Carr\'ee geometry on a staggered grid. The domain over which solutions are found is a subset of the global spherical surface, defined by user-specified limits of colatitude and longitude. Our staggered-grid approach, based on that of Yee (1966), is more conservative and self-consistent compared to the centered, Cartesian grid used by Kazachenko et al. (2014). The library can be used to compute an end-to-end solution for electric fields from data taken by the HMI instrument aboard NASA's SDO Mission. This capability has been incorporated into the HMI pipeline processing system operating at SDO's JSOC. The library is written in a general and modular way so that the calculations can be customized to modify or delete electric field contributions, or used with other data sets. Other applications include "nudging" numerical models of the solar atmosphere to facilitate assimilative simulations. The library includes an ability to compute "global" (whole-Sun) electric field solutions. The library also includes an ability to compute Potential Magnetic Field solutions in spherical coordinates. This distribution includes a number of test programs which allow the user to test the software.

Observations from the Kepler mission have revealed frequent superflares on young and active solar-like stars. Superflares result from the large-scale restructuring of stellar magnetic fields, and are associated with the eruption of coronal material (a coronal mass ejection, or CME) and energy release that can be orders of magnitude greater than those observed in the largest solar flares. These catastrophic events, if frequent, can significantly impact the potential habitability of terrestrial exoplanets through atmospheric erosion or intense radiation exposure at the surface. We present results from numerical modeling designed to understand how an eruptive superflare from a young solar-type star, $\kappa^{1}Cet$, could occur and would impact its astrospheric environment. Our data-inspired, three-dimensional magnetohydrodynamic modeling shows that global-scale shear concentrated near the radial-field polarity inversion line can energize the closed-field stellar corona sufficiently to power a global, eruptive superflare that releases approximately the same energy as the extreme 1859 Carrington event from the Sun. We examine proxy measures of synthetic emission during the flare and estimate the observational signatures of our CME-driven shock, both of which could have extreme space-weather impacts on the habitability of any Earth-like exoplanets. We also speculate that the observed 1986 Robinson-Bopp superflare from $\kappa^{1}Cet$ was perhaps as extreme for that star as the Carrington flare was for the Sun.

In this work we study how the input data cadence affects the photospheric energy and helicity injection estimates in eruptive NOAA active region 11158. We sample the novel 2.25-minute vector magnetogram and Dopplergram data from the \emph{Helioseismic and Magnetic Imager} (HMI) instrument onboard the \emph{Solar Dynamics Observatory} (SDO) spacecraft to create input datasets of variable cadences ranging from 2.25 minutes to 24 hours. We employ state-of-the-art data processing, velocity and electric field inversion methods for deriving estimates of the energy and helicity injections from these datasets. We find that the electric field inversion methods that reproduce the observed magnetic field evolution through the use of Faraday's law are more stable against variable cadence: the PDFI (PTD-Doppler-FLCT-Ideal) electric field inversion method produces consistent injection estimates for cadences from 2.25 minutes up to 2 hours, implying that the photospheric processes acting on time scales below 2 hours contribute little to the injections, or that they are below the sensitivity of the input data and the PDFI method. On other hand, the electric field estimate derived from the output of DAVE4VM (Differential Affine Velocity Estimator for Vector Magnetograms), which does not fulfil Faraday's law exactly, produces significant variations in the energy and helicity injection estimates in the 2.25-minute to 2-hour cadence range. We present also a third, novel DAVE4VM-based electric field estimate, which corrects the poor inductivity of the raw DAVE4VM estimate. This method is less sensitive to the changes of cadence, but still faces significant issues for the lowest of considered cadences ($\geq$2 hours). We find several potential problems in both PDFI- and DAVE4VM-based injection estimates and conclude that the quality of both should be surveyed further in controlled environments.

Active Regions (ARs) that exhibit compact Polarity Inversion Lines (PILs) are known to be very flare-productive. However, the physical mechanisms behind this statistical inference have not been demonstrated conclusively. We show that such PILs can occur due to the collision between two emerging flux tubes nested within the same AR. In such multipolar ARs, the flux tubes may emerge simultaneously or sequentially, each initially producing a bipolar magnetic region (BMR) at the surface. During each flux tube's emergence phase, the magnetic polarities can migrate such that opposite polarities belonging to different BMRs collide, resulting in shearing and cancellation of magnetic flux. We name this process 'collisional shearing' to emphasize that the shearing and flux cancellation develops due to the collision. Collisional shearing is a process different from the known concept of flux cancellation occurring between polarities of a single bipole, a process that has been commonly used in many numerical models. High spatial and temporal resolution observations from the Solar Dynamics Observatory for two emerging ARs, AR11158 and AR12017, show the continuous cancellation of up to 40% of the unsigned magnetic flux of the smallest BMR, which occurs at the collisional PIL for as long as the collision persists. The flux cancellation is accompanied by a succession of solar flares and CMEs, products of magnetic reconnection along the collisional PIL. Our results suggest that the quantification of magnetic cancellation driven by collisional shearing needs to be taken into consideration in order to improve the prediction of solar energetic events and space weather.

We present a database of 3137 solar flare ribbon events corresponding to every flare of GOES class C1.0 and greater within 45 degrees from the central meridian, from April 2010 until April 2016, observed by the \emph{Solar Dynamics Observatory}. For every event in the database, we compare the GOES peak X-ray flux with corresponding active-region and flare-ribbon properties. We find that while the peak X-ray flux is not correlated with the active region unsigned magnetic flux, it is strongly correlated with the flare ribbon reconnection flux, flare ribbon area, and the fraction of active region flux that undergoes reconnection. We find the relationship between the peak X-ray flux and the flare ribbon reconnection flux to be $I_\mathrm{X,peak} \propto \Phi_\mathrm{ribbon}^{1.5}$. This scaling law is consistent with earlier hydrodynamic simulations of impulsively heated flare loops. Using the flare reconnection flux as a proxy for the total released flare energy $E$, we find that the occurrence frequency of flare energies follows a power-law dependence: $dN/dE \propto E^{-1.6}$ for $10^{31}<E<10^{33}$ erg, consistent with earlier studies of solar and stellar flares. The database is available online and can be used for future quantitative studies of flares.

The solar active region photospheric magnetic field evolves rapidly during major eruptive events, suggesting appreciable feedback from the corona. Previous studies of these "magnetic imprints" are mostly based on line-of-sight only or lower-cadence vector observations; a temporally resolved depiction of the vector field evolution is hitherto lacking. Here, we introduce the high-cadence (90~s or 135~s) vector magnetogram dataset from the Helioseismic and Magnetic Imager (HMI) that is well suited for investigating the phenomenon. These observations allow quantitative characterization of the permanent, step-like changes that are most pronounced in the horizontal field component ($B_h$). A highly structured pattern emerges from analysis of an archetypical event, \texttt{SOL2011-02-15T01:56}, where $B_h$ near the main polarity inversion line increases significantly during the earlier phase of the associated flare with a time scale of several minutes, while $B_h$ in the periphery decreases at later times with smaller magnitudes and a slightly longer time scale. The dataset also allows effective identification of the "magnetic transient" artifact, where enhanced flare emission alters the Stokes profiles and the inferred magnetic field becomes unreliable. Our results provide insights on the momentum processes in solar eruptions. The dataset may also be useful to the study of sunquakes and data-driven modeling of the corona.

We present a detailed analysis of the properties of magnetic reconnection at large-scale current sheets in a high cadence version of the Lynch & Edmondson (2013) 2.5D MHD simulation of sympathetic magnetic breakout eruptions from a pseudostreamer source region. We examine the resistive tearing and breakup of the three main current sheets into chains of X- and O-type null points and follow the dynamics of magnetic island growth, their merging, transit, and ejection with the reconnection exhaust. For each current sheet, we quantify the evolution of the length-to-width aspect ratio (up to $\sim$100:1), Lundquist number ($\sim$10$^3$), and reconnection rate (inflow-to-outflow ratios reaching $\sim$0.40). We examine the statistical and spectral properties of the fluctuations in the current sheets resulting from the plasmoid instability, including the distribution of magnetic island area, mass, and flux content. We show that the temporal evolution of the spectral index of the reconnection-generated magnetic energy density fluctuations appear to reflect global properties of the current sheet evolution. Our results are in excellent agreement with recent, high resolution reconnection-in-a-box simulations even though our current sheets' formation, growth, and dynamics are intrinsically coupled to the global evolution of sequential sympathetic CME eruptions.

How much electromagnetic energy crosses the photosphere in evolving solar active regions? With the advent of high-cadence vector magnetic field observations, addressing this fundamental question has become tractable. In this paper, we apply the "PTD-Doppler-FLCT-Ideal" (PDFI) electric field inversion technique of Kazachenko et al. (2014) to a 6-day HMI/SDO vector magnetogram and Doppler velocity sequence, to find the electric field and Poynting flux evolution in active region NOAA 11158, which produced an X2.2 flare early on 2011 February 15. We find photospheric electric fields ranging up to $2$ V/cm. The Poynting fluxes range from $[-0.6$ to $2.3]\times10^{10}$ ergs$\cdot$cm$^{-2}$s$^{-1}$, mostly positive, with the largest contribution to the energy budget in the range of $[10^9$-$10^{10}]$ ergs$\cdot$cm$^{-2}$s$^{-1}$. Integrating the instantaneous energy flux over space and time, we find that the total magnetic energy accumulated above the photosphere from the initial emergence to the moment before the X2.2 flare to be $E=10.6\times10^{32}$ ergs, which is partitioned as $2.0$ and $8.6\times10^{32}$ ergs, respectively, between free and potential energies. Those estimates are consistent with estimates from preflare non-linear force-free field (NLFFF) extrapolations and the Minimum Current Corona estimates (MCC), in spite of our very different approach. This study of photospheric electric fields demonstrates the potential of the PDFI approach for estimating Poynting fluxes and opens the door to more quantitative studies of the solar photosphere and more realistic data-driven simulations of coronal magnetic field evolution.

The most violent space weather events (eruptive solar flares and coronal mass ejections) are driven by the release of free magnetic energy stored in the solar corona. Energy can build up on timescales of hours to days, and then may be suddenly released in the form of a magnetic eruption, which then propagates through interplanetary space, possibly impacting the Earth's space environment. Can we use the observed evolution of the magnetic and velocity fields in the solar photosphere to model the evolution of the overlying solar coronal field, including the storage and release of magnetic energy in such eruptions? The objective of CGEM, the Coronal Global Evolutionary Model, funded by the NASA/NSF Space Weather Modeling program, is to develop and evaluate such a model for the evolution of the coronal magnetic field. The evolving coronal magnetic field can then be used as a starting point for magnetohydrodynamic (MHD) models of the corona, which can then be used to drive models of heliospheric evolution and predictions of magnetic field and plasma density conditions at 1AU.

Solar Active Region NOAA 11158 has hosted a number of strong flares, including one X2.2 event. The complexity of current density and current helicity are studied through cancellation analysis of their sign-singular measure, which features power-law scaling. Spectral analysis is also performed, revealing the presence of two separate scaling ranges with different spectral index. The time evolution of parameters is discussed. Sudden changes of the cancellation exponents at the time of large flares, and the presence of correlation with EUV and X-ray flux, suggest that eruption of large flares can be linked to the small scale properties of the current structures.

Photospheric electric fields, estimated from sequences of vector magnetic field and Doppler measurements, can be used to estimate the flux of magnetic energy (the Poynting flux) into the corona and as time-dependent boundary conditions for dynamic models of the coronal magnetic field. We have modified and extended an existing method to estimate photospheric electric fields that combines a poloidal-toroidal (PTD) decomposition of the evolving magnetic field vector with Doppler and horizontal plasma velocities. Our current, more comprehensive method, which we dub the "{\bf P}TD-{\bf D}oppler-{\bf F}LCT {\bf I}deal" (PDFI) technique, can now incorporate Doppler velocities from non-normal viewing angles. It uses the \texttt{FISHPACK} software package to solve several two-dimensional Poisson equations, a faster and more robust approach than our previous implementations. Here, we describe systematic, quantitative tests of the accuracy and robustness of the PDFI technique using synthetic data from anelastic MHD (\texttt{ANMHD}) simulations, which have been used in similar tests in the past. We find that the PDFI method has less than $1%$ error in the total Poynting flux and a $10%$ error in the helicity flux rate at a normal viewing angle $(\theta=0$) and less than $25%$ and $10%$ errors respectively at large viewing angles ($\theta<60^\circ$). We compare our results with other inversion methods at zero viewing angle, and find that our method's estimates of the fluxes of magnetic energy and helicity are comparable to or more accurate than other methods. We also discuss the limitations of the PDFI method and its uncertainties.

n order to better understand the solar genesis of interplanetary magnetic clouds (MCs) we model the magnetic and topological properties of four large eruptive solar flares and relate them to observations. We use the three-dimensional Minimum Current Corona model \cite{Longcope1996d} and observations of pre-flare photospheric magnetic field and flare ribbons to derive values of reconnected magnetic flux, flare energy, flux rope helicity and orientation of the flux rope poloidal field. We compare model predictions of those quantities to flare and MC observations and within the estimated uncertainties of the methods used find the following. The predicted model reconnection fluxes are equal to or lower than the reconnection fluxes inferred from the observed ribbon motions. Both observed and model reconnection fluxes match the MC poloidal fluxes. The predicted flux rope helicities match the MC helicities. The predicted free energies lie between the observed energies and the estimated total flare luminosities. The direction of the leading edge of the MC's poloidal field is aligned with the poloidal field of the flux rope in the AR rather than the global dipole field. These findings compel us to believe that magnetic clouds associated with these four solar flares are formed by low-corona magnetic reconnection during the eruption, rather than eruption of pre-existing structures in the corona or formation in the upper corona with participation of the global magnetic field. We also note that since all four flares occurred in active regions without significant pre-flare flux emergence and cancellation, the energy and helicity we find are stored by shearing and rotating motions, which are sufficient to account for the observed radiative flare energy and MC helicity.