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.

Emanating from the base of the Sun's corona, the solar wind fills the interplanetary medium with a magnetized stream of charged particles whose interaction with the Earth's magnetosphere has space-weather consequences such as geomagnetic storms. Accurately predicting the solar wind through measurements of the spatio-temporally evolving conditions in the solar atmosphere is important but remains an unsolved problem in heliophysics and space-weather research. In this work, we use deep learning for prediction of solar wind (SW) properties. We use Extreme Ultraviolet images of the solar corona from space based observations to predict the SW speed from the NASA OMNIWEB dataset, measured at Lagragian point 1. We evaluate our model against autoregressive and naive models, and find that our model outperforms the benchmark models, obtaining a best-fit correlation of 0.55 $\pm$ 0.03 with the observed data. Upon visualization and investigation of how the model uses data to make predictions, we find higher activation at the coronal holes for fast wind prediction ($\approx$ 3 to 4 days prior to prediction), and at the active regions for slow wind prediction. These trends bear an uncanny similarity to the influence of regions potentially being the sources of fast and slow wind, as reported in literature. This suggests that our model was able to learn some of the salient associations between coronal and solar wind structure without built-in physics knowledge. Such an approach may help us discover hitherto unknown relationships in heliophysics datasets.

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.

Solar coronal dimmings have been observed extensively in the past two decades and are believed to have close association with coronal mass ejections (CMEs). Recent study found that coronal dimming is the only signature that could differentiate powerful ares that have CMEs from those that do not. Therefore, dimming might be one of the best candidates to observe the stellar CMEs on distant Sun-like stars. In this study, we investigate the possibility of using coronal dimming as a proxy to diagnose stellar CMEs. By simulating a realistic solar CME event and corresponding coronal dimming using a global magnetohydrodynamics model (AWSoM: Alfven-wave Solar Model), we first demonstrate the capability of the model to reproduce solar observations. We then extend the model for simulating stellar CMEs by modifying the input magnetic flux density as well as the initial magnetic energy of the CME flux rope. Our result suggests that with improved instrument sensitivity, it is possible to detect the coronal dimming signals induced by the stellar CMEs.

For a better understanding of magnetic field in the solar corona and dynamic activities such as flares and coronal mass ejections, it is crucial to measure the time-evolving coronal field and accurately estimate the magnetic energy. Recently, a new modeling technique called the data-driven coronal field model, in which the time evolution of magnetic field is driven by a sequence of photospheric magnetic and velocity field maps, has been developed and revealed the dynamics of flare-productive active regions. Here we report on the first qualitative and quantitative assessment of different data-driven models using a magnetic flux emergence simulation as a ground-truth (GT) data set. We compare the GT field with those reconstructed from the GT photospheric field by four data-driven algorithms. It is found that, at least, the flux rope structure is reproduced in all coronal field models. Quantitatively, however, the results show a certain degree of model dependence. In most cases, the magnetic energies and relative magnetic helicity are comparable to or at most twice of the GT values. The reproduced flux ropes have a sigmoidal shape (consistent with GT) of various sizes, a vertically-standing magnetic torus, or a packed structure. The observed discrepancies can be attributed to the highly non-force-free input photospheric field, from which the coronal field is reconstructed, and to the modeling constraints such as the treatment of background atmosphere, the bottom boundary setting, and the spatial resolution.

Understanding and monitoring the complex and dynamic processes of the Sun is important for a number of human activities on Earth and in space. For this reason, NASA's Solar Dynamics Observatory (SDO) has been continuously monitoring the multi-layered Sun's atmosphere in high-resolution since its launch in 2010, generating terabytes of observational data every day. The synergy between machine learning and this enormous amount of data has the potential, still largely unexploited, to advance our understanding of the Sun and extend the capabilities of heliophysics missions. In the present work, we show that deep learning applied to SDO data can be successfully used to create a high-fidelity virtual telescope that generates synthetic observations of the solar corona by image translation. Towards this end we developed a deep neural network, structured as an encoder-decoder with skip connections (U-Net), that reconstructs the Sun's image of one instrument channel given temporally aligned images in three other channels. The approach we present has the potential to reduce the telemetry needs of SDO, enhance the capabilities of missions that have less observing channels, and transform the concept development of future missions.

As a part of NASA's Heliophysics System Observatory (HSO) fleet of satellites,the Solar Dynamics Observatory (SDO) has continuously monitored the Sun since2010. Ultraviolet (UV) and Extreme UV (EUV) instruments in orbit, such asSDO's Atmospheric Imaging Assembly (AIA) instrument, suffer time-dependent degradation which reduces instrument sensitivity. Accurate calibration for (E)UV instruments currently depends on periodic sounding rockets, which are infrequent and not practical for heliophysics missions in deep space. In the present work, we develop a Convolutional Neural Network (CNN) that auto-calibrates SDO/AIA channels and corrects sensitivity degradation by exploiting spatial patterns in multi-wavelength observations to arrive at a self-calibration of (E)UV imaging instruments. Our results remove a major impediment to developing future HSOmissions of the same scientific caliber as SDO but in deep space, able to observe the Sun from more vantage points than just SDO's current geosynchronous orbit.This approach can be adopted to perform autocalibration of other imaging systems exhibiting similar forms of degradation

The Multi-slit Solar Explorer (MUSE) is a proposed mission aimed at understanding the physical mechanisms driving the heating of the solar corona and the eruptions that are at the foundation of space weather. MUSE contains two instruments, a multi-slit EUV spectrograph and a context imager. It will simultaneously obtain EUV spectra (along 37 slits) and context images with the highest resolution in space (0.33-0.4 arcsec) and time (1-4 s) ever achieved for the transition region and corona. The MUSE science investigation will exploit major advances in numerical modeling, and observe at the spatial and temporal scales on which competing models make testable and distinguishable predictions, thereby leading to a breakthrough in our understanding of coronal heating and the drivers of space weather. By obtaining spectra in 4 bright EUV lines (Fe IX 171A, Fe XV 284A, Fe XIX-XXI 108A) covering a wide range of transition region and coronal temperatures along 37 slits simultaneously, MUSE will be able to "freeze" the evolution of the dynamic coronal plasma. We describe MUSE's multi-slit approach and show that the optimization of the design minimizes the impact of spectral lines from neighboring slits, generally allowing line parameters to be accurately determined. We also describe a Spectral Disambiguation Code to resolve multi-slit ambiguity in locations where secondary lines are bright. We use simulations of the corona and eruptions to perform validation tests and show that the multi-slit disambiguation approach allows accurate determination of MUSE observables in locations where significant multi-slit contamination occurs.

Because of the complex physics that governs the formation of chromospheric lines, interpretation of solar chromospheric observations is difficult. The origin and characteristics of many chromospheric features are, because of this, unresolved. We focus here on studying two prominent features: long fibrils and flare ribbons. To model them, we use a 3D MHD simulation of an active region which self-consistently reproduces both of them. We model the H$\alpha$, Mg II k, Ca II K, and Ca II 8542 {\AA} lines using the 3D non-LTE radiative transfer code Multi3D. This simulation reproduces long fibrils that span between the opposite-polarity sunspots and go up to 4 Mm in height. They can be traced in all lines due to density corrugation. Opposite to previous studies, H$\alpha$, Mg II h&k, and Ca II H&K, are formed at similar height in this model. Magnetic field lines are aligned with the H$\alpha$ fibrils, but the latter holds to a lesser extent for the Ca II 8542 {\AA} line. The simulation shows structures in the H$\alpha$ line core that look like flare ribbons. The emission in the ribbons is caused by a dense chromosphere and a transition region at high column mass. The ribbons are visible in all chromospheric lines, but least prominent in Ca II 8542 {\AA} line. In some pixels, broad asymmetric profiles with a single emission peak are produced, similar to the profiles observed in flare ribbons. They are caused by a deep onset of the chromospheric temperature rise and large velocity gradients. The simulation produces long fibrils similar to what is seen in observations. It also produces structures similar to flare ribbons despite the lack of non-thermal electrons in the simulation. The latter suggests that thermal conduction might be a significant agent in transporting flare energy to the chromosphere in addition to non-thermal electrons.

Solar flares - bursts of high-energy radiation responsible for severe space-weather effects - are a consequence of the occasional destabilization of magnetic fields rooted in active regions (ARs). The complexity of AR evolution is a barrier to a comprehensive understanding of flaring processes and accurate prediction. Though machine learning (ML) has been used to improve flare predictions, the potential for revealing precursors and associated physics has been underexploited. Here, we train ML algorithms to classify between vector-magnetic-field observations from flaring ARs, producing at least one M-/X-class flare, and non-flaring ARs. Analysis of magnetic-field observations accurately classified by the machine presents statistical evidence for (1) ARs persisting in flare-productive states --- characterized by AR area --- for days, before and after M- and X-class flare events, (2) systematic pre-flare build-up of free energy in the form of electric currents, suggesting that associated subsurface magnetic field is twisted, (3) intensification of Maxwell stresses in the corona above newly emerging ARs, days before first flares. These results provide new insights into flare physics and improving flare forecasting.

In this paper we present a curated dataset from the NASA Solar Dynamics Observatory (SDO) mission in a format suitable for machine learning research. Beginning from level 1 scientific products we have processed various instrumental corrections, downsampled to manageable spatial and temporal resolutions, and synchronized observations spatially and temporally. We illustrate the use of this dataset with two example applications: forecasting future EVE irradiance from present EVE irradiance and translating HMI observations into AIA observations. For each application we provide metrics and baselines for future model comparison. We anticipate this curated dataset will facilitate machine learning research in heliophysics and the physical sciences generally, increasing the scientific return of the SDO mission. This work is a direct result of the 2018 NASA Frontier Development Laboratory Program. Please see the appendix for access to the dataset.

The signal measured by an astronomical spectrometer may be due to radiation from a multi-component mixture of plasmas with a range of physical properties (e.g. temperature, Doppler velocity). Confusion between multiple components may be exacerbated if the spectrometer sensor is illuminated by overlapping spectra dispersed from different slits, with each slit being exposed to radiation from a different portion of an extended astrophysical object. We use a compressed sensing method to robustly retrieve the different components. This method can be adopted for a variety of spectrometer configurations, including single-slit, multi-slit (e.g., the proposed MUlti-slit Solar Explorer mission; MUSE) and slot spectrometers (which produce overlappograms).

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.

The term "ultraviolet (UV) burst" is introduced to describe small, intense, transient brightenings in ultraviolet images of solar active regions. We inventorize their properties and provide a definition based on image sequences in transition-region lines. Coronal signatures are rare, and most bursts are associated with small-scale, canceling opposite-polarity fields in the photosphere that occur in emerging flux regions, moving magnetic features in sunspot moats, and sunspot light bridges. We also compare UV bursts with similar transition-region phenomena found previously in solar ultraviolet spectrometry and with similar phenomena at optical wavelengths, in particular Ellerman bombs. Akin to the latter, UV bursts are probably small-scale magnetic reconnection events occurring in the low atmosphere, at photospheric and/or chromospheric heights. Their intense emission in lines with optically thin formation gives unique diagnostic opportunities for studying the physics of magnetic reconnection in the low solar atmosphere. This paper is a review report from an International Space Science Institute team that met in 2016-2017.

We report SDO/AIA observations of an extraordinary global extreme ultraviolet (EUV) wave triggered by the X8.2+ flare-CME eruption on 2017 September 10. This was one of the best EUV waves ever observed with modern instruments, yet it was likely the last one of such magnitudes of Solar Cycle 24 as the Sun heads toward the minimum. Its remarkable characteristics include the following. (1) The wave was observed, for the first time, to traverse the full-Sun corona over the entire visible solar disk and off-limb circumference, manifesting a truly global nature, owing to its exceptionally large amplitude, e.g., with EUV enhancements by up to 300% at 1.1 Rsun from the eruption. (2) This leads to strong transmissions (in addition to commonly observed reflections) in and out of both polar coronal holes, which are usually devoid of EUV waves. It has elevated wave speeds >2000 km/s within them, consistent with the expected higher fast-mode magnetosonic wave speeds. The coronal holes essentially serve as new "radiation centers" for the waves being refracted out of them, which then travel toward the equator and collide head-on, causing additional EUV enhancements. (3) The wave produces significant compressional heating to local plasma upon its impact, indicated by long-lasting EUV intensity changes and differential emission measure increases at higher temperatures (e.g., log T=6.2) accompanied by decreases at lower temperatures (e.g., log T=6.0). These characteristics signify the potential of such EUV waves for novel magnetic and thermal diagnostics of the solar corona {\it on global scales}.

Coronal mass ejections (CMEs) are large-scale eruptions of plasma from the coronae of stars. Understanding the plasma processes involved in CME initiation has applications to space weather forecasting and laboratory plasma experiments. James et al. (Sol. Phys. 292, 71, 2017) used EUV observations to conclude that a magnetic flux rope formed in the solar corona above NOAA Active Region 11504 before it erupted on 14 June 2012 (SOL2012-06-14). In this work, we use data from the Solar Dynamics Observatory to model the coronal magnetic field of the active region one hour prior to eruption using a nonlinear force-free field extrapolation, and find a flux rope reaching a maximum height of 150 Mm above the photosphere. Estimations of the average twist of the strongly asymmetric extrapolated flux rope are between 1.35 and 1.88 turns, depending on the choice of axis, although the erupting structure was not observed to kink. The decay index near the apex of the axis of the extrapolated flux rope is comparable to typical critical values required for the onset of the torus instability, so we suggest that the torus instability drove the eruption.

Self-organization is a property of dissipative nonlinear processes that are governed by an internal driver and a positive feedback mechanism, which creates regular geometric and/or temporal patterns and decreases the entropy, in contrast to random processes. Here we investigate for the first time a comprehensive number of 16 self-organization processes that operate in planetary physics, solar physics, stellar physics, galactic physics, and cosmology. Self-organizing systems create spontaneous {\sl order out of chaos}, during the evolution from an initially disordered system to an ordered stationary system, via quasi-periodic limit-cycle dynamics, harmonic mechanical resonances, or gyromagnetic resonances. The internal driver can be gravity, rotation, thermal pressure, or acceleration of nonthermal particles, while the positive feedback mechanism is often an instability, such as the magneto-rotational instability, the Rayleigh-B\'enard convection instability, turbulence, vortex attraction, magnetic reconnection, plasma condensation, or loss-cone instability. Physical models of astrophysical self-organization processes involve hydrodynamic, MHD, and N-body formulations of Lotka-Volterra equation systems.

We report the intriguing large-scale dynamic phenomena associated with the M6.5 flare~(SOL2015-06-22T18:23) in NOAA active region 12371, observed by RHESSI, Fermi, and the Atmospheric Image Assembly (AIA) and Magnetic Imager (HMI) on the Solar Dynamic Observatory (SDO). The most interesting feature of this event is a third ribbon (R3) arising in the decay phase, propagating along a dimming channel (seen in EUV passbands) towards a neighboring sunspot. The propagation of R3 occurs in the presence of hard X-ray footpoint emission, and is broadly visible at temperatures from 0.6~MK to over 10~MK through the Differential Emission Measure (DEM) analysis. The coronal loops then undergo an apparent slipping motion following the same path of R3, after a $\sim$ 80~min delay. To understand the underlying physics, we investigate the magnetic configuration and the thermal structure of the flaring region. Our results are in favor of a slipping-type reconnection followed by the thermodynamic evolution of coronal loops. In comparison with those previously reported slipping reconnection events, this one proceeds across a particularly long distance ($\sim$60 Mm) over a long period of time ($\sim$50 min), and shows two clearly distinguished phases: the propagation of the footpoint brightening driven by nonthermal particle injection and the apparent slippage of loops governed by plasma heating and subsequent cooling.

Emerging flux regions (EFRs) are known to exhibit various sporadic local heating events in the lower atmosphere. To investigate the characteristics of these events, especially to link the photospheric magnetic fields and atmospheric dynamics, we analyze Hinode, IRIS, and SDO data of a new EFR in NOAA AR 12401. Out of 151 bright points (BPs) identified in Hinode/SOT Ca images, 29 are overlapped by an SOT/SP scan. Seven BPs in the EFR center possess mixed-polarity magnetic backgrounds in the photosphere. Their IRIS UV spectra (e.g., Si IV 1402.8 A) are strongly enhanced and red- or blue-shifted with tails reaching +/- 150 km/s, which is highly suggestive of bi-directional jets, and each brightening lasts for 10 - 15 minutes leaving flare-like light curves. Most of this group show bald patches, the U-shaped photospheric magnetic loops. Another 10 BPs are found in unipolar regions at the EFR edges. They are generally weaker in UV intensities and exhibit systematic redshifts with Doppler speeds up to 40 km/s, which could exceed the local sound speed in the transition region. Both types of BPs show signs of strong temperature increase in the low chromosphere. These observational results support the physical picture that heating events in the EFR center are due to magnetic reconnection within cancelling undular fields like Ellerman bombs, while the peripheral heating events are due to shocks or strong compressions caused by fast downflows along the overlying arch filament system.

With the global view and high-cadence observations from SDO/AIA and STEREO, many spatially separated solar eruptive events appear to be coupled. However, the mechanisms for "sympathetic" events are still largely unknown. In this study, we investigate the impact of an erupting flux rope on surrounding solar structures through large-scale magnetic coupling. We build a realistic environment of the solar corona on 2011 February 15 using a global magnetohydrodynamics (MHD) model and initiate coronal mass ejections (CMEs) in active region (AR) 11158 by inserting Gibson-Low analytical flux ropes. We show that a CME's impact on the surrounding structures depends not only on the magnetic strength of these structures and their distance to the source region, but also on the interaction between the CME with the large-scale magnetic field. Within the CME expansion domain where the flux rope field directly interacts with the solar structures, expansion-induced reconnection often modifies the overlying field, thereby increasing the decay index. This effect may provide a primary coupling mechanism underlying the sympathetic eruptions. The magnitude of the impact is found to depend on the orientation of the erupting flux rope, with the largest impacts occurring when the flux rope is favorably oriented for reconnecting with the surrounding regions. Outside the CME expansion domain, the influence of the CME is mainly through field line compression or post-eruption relaxation. Based on our numerical experiments, we discuss a way to quantify the eruption impact, which could be useful for forecasting purposes.

We carried out a magnetohydrodynamics simulation where a subsurface twisted kink-unstable flux tube emerges from the solar interior to the corona. Unlike the previous expectations based on the bodily emergence of a knotted tube, we found that the kinked tube can spontaneously form a complex quadrupole structure at the photosphere. Due to the development of the kink instability before the emergence, the magnetic twist at the kinked apex of the tube is greatly reduced, although the other parts of the tube is still strongly twisted. This leads to the formation of a complex quadrupole structure: a pair of the coherent, strongly twisted spots and a narrow complex bipolar pair between it. The quadrupole is formed by the submergence of a portion of emerged magnetic fields. This result is relevant for understanding of the origin of the complex multipolar $\delta$-spot regions that have a strong magnetic shear and emerge with polarity orientations not following Hale-Nicholson and Joy Laws.

Light bridges, the bright structure dividing umbrae in sunspot regions, show various activity events. In Paper I, we reported on analysis of multi-wavelength observations of a light bridge in a developing active region (AR) and concluded that the activity events are caused by magnetic reconnection driven by magnetconvective evolution. The aim of this second paper is to investigate the detailed magnetic and velocity structures and the formation mechanism of light bridges. For this purpose, we analyze numerical simulation data from a radiative magnetohydrodynamics model of an emerging AR. We find that a weakly-magnetized plasma upflow in the near-surface layers of the convection zone is entrained between the emerging magnetic bundles that appear as pores at the solar surface. This convective upflow continuously transports horizontal fields to the surface layer and creates a light bridge structure. Due to the magnetic shear between the horizontal fields of the bridge and the vertical fields of the ambient pores, an elongated cusp-shaped current layer is formed above the bridge, which may be favorable for magnetic reconnection. The striking correspondence between the observational results of Paper I and the numerical results of this paper provides a consistent physical picture of light bridges. The dynamic activity phenomena occur as a natural result of the bridge formation and its convective nature, which has much in common with those of umbral dots and penumbral filaments.

Light bridges, the bright structures that divide the umbra of sunspots and pores into smaller pieces, are known to produce wide variety of activity events in solar active regions (ARs). It is also known that the light bridges appear in the assembling process of nascent sunspots. The ultimate goal of this series of papers is to reveal the nature of light bridges in developing ARs and the occurrence of activity events associated with the light bridge structures from both observational and numerical approaches. In this first paper, exploiting the observational data obtained by Hinode, IRIS, and Solar Dynamics Observatory (SDO), we investigate the detailed structure of the light bridge in NOAA AR 11974 and its dynamic activity phenomena. As a result, we find that the light bridge has a weak, horizontal magnetic field, which is transported from the interior by large-scale convective upflow and is surrounded by strong, vertical fields of adjacent pores. In the chromosphere above the bridge, a transient brightening occurs repeatedly and intermittently, followed by a recurrent dark surge ejection into higher altitudes. Our analysis indicates that the brightening is the plasma heating due to magnetic reconnection at lower altitudes, while the dark surge is the cool, dense plasma ejected from the reconnection region. From the observational results, we conclude that the dynamic activity observed in a light bridge structure such as chromospheric brightenings and dark surge ejections are driven by magnetoconvective evolution within the light bridge and its interaction with surrounding magnetic fields.

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.

We study the buoyant rise of magnetic flux tubes embedded in an adiabatic stratification using two-and three-dimensional, MHD simulations. We analyze the dependence of the tube evolution on the field line twist and on the curvature of the tube axis in different diffusion regimes. To be able to achieve a comparatively high spatial resolution we use the FLASH code, which has a built-in Adaptive Mesh Refinement (AMR) capability. Our 3D experiments reach Reynolds numbers that permit a reasonable comparison of the results with those of previous 2D simulations. When the experiments are run without AMR, hence with a comparatively large diffusivity, the amount of longitudinal magnetic flux retained inside the tube increases with the curvature of the tube axis. However, when a low-diffusion regime is reached by using the AMR algorithms, the magnetic twist is able to prevent the splitting of the magnetic loop into vortex tubes and the loop curvature does not play any significant role. We detect the generation of vorticity in the main body of the tube of opposite sign on the opposite sides of the apex. This is a consequence of the inhomogeneity of the azimuthal component of the field on the flux surfaces. The lift force associated with this global vorticity makes the flanks of the tube move away from their initial vertical plane in an antisymmetric fashion. The trajectories have an oscillatory motion superimposed, due to the shedding of vortex rolls to the wake, which creates a Von Karman street.

The outer solar atmosphere, the corona, contains plasma at temperatures of more than a million K, more than 100 times hotter that solar surface. How this gas is heated is a fundamental question tightly interwoven with the structure of the magnetic field in the upper atmosphere. Conducting numerical experiments based on magnetohydrodynamics we account for both the evolving three-dimensional structure of the atmosphere and the complex interaction of magnetic field and plasma. Together this defines the formation and evolution of coronal loops, the basic building block prominently seen in X-rays and extreme ultraviolet (EUV) images. The structures seen as coronal loops in the EUV can evolve quite differently from the magnetic field. While the magnetic field continuously expands as new magnetic flux emerges through the solar surface, the plasma gets heated on successively emerging fieldlines creating an EUV loop that remains roughly at the same place. For each snapshot the EUV images outline the magnetic field, but in contrast to the traditional view, the temporal evolution of the magnetic field and the EUV loops can be different. Through this we show that the thermal and the magnetic evolution in the outer atmosphere of a cool star has to be treated together, and cannot be simply separated as done mostly so far.

We present a new method for performing differential emission measure (DEM) inversions on narrow-band EUV images from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). The method yields positive definite DEM solutions by solving a linear program. This method has been validated against a diverse set of thermal models of varying complexity and realism. These include (1) idealized gaussian DEM distributions, (2) 3D models of NOAA Active Region 11158 comprising quasi-steady loop atmospheres in a non-linear force-free field, and (3) thermodynamic models from a fully-compressible, 3D MHD simulation of AR corona formation following magnetic flux emergence. We then present results from the application of the method to AIA observations of Active Region 11158, comparing the region's thermal structure on two successive solar rotations. Additionally, we show how the DEM inversion method can be adapted to simultaneously invert AIA and XRT data, and how supplementing AIA data with the latter improves the inversion result. The speed of the method allows for routine production of DEM maps, thus facilitating science studies that require tracking of the thermal structure of the solar corona in time and space.

We report on observations of recurrent jets by instruments onboard the Interface Region Imaging Spectrograph (IRIS), Solar Dynamics Observatory (SDO) and Hinode spacecrafts. Over a 4-hour period on July 21st 2013, recurrent coronal jets were observed to emanate from NOAA Active Region 11793. FUV spectra probing plasma at transition region temperatures show evidence of oppositely directed flows with components reaching Doppler velocities of +/- 100 km/s. Raster Doppler maps using a Si IV transition region line show all four jets to have helical motion of the same sense. Simultaneous observations of the region by SDO and Hinode show that the jets emanate from a source region comprising a pore embedded in the interior of a supergranule. The parasitic pore has opposite polarity flux compared to the surrounding network field. This leads to a spine-fan magnetic topology in the coronal field that is amenable to jet formation. Time-dependent data-driven simulations are used to investigate the underlying drivers for the jets. These numerical experiments show that the emergence of current-carrying magnetic field in the vicinity of the pore supplies the magnetic twist needed for recurrent helical jet formation.

We present the first results of sunspot oscillations from observations by the Interface Region Imaging Spectrograph. The strongly nonlinear oscillation is identified in both the slit-jaw images and the spectra of several emission lines formed in the transition region and chromosphere. We first apply a single Gaussian fit to the profiles of the Mgii 2796.35 {\AA}, Cii 1335.71 {\AA}, and Si iv 1393.76 {\AA} lines in the sunspot. The intensity change is about 30%. The Doppler shift oscillation reveals a sawtooth pattern with an amplitude of about 10 km/s in Si iv. In the umbra the Si iv oscillation lags those of Cii and Mgii by about 3 and 12 s, respectively. The line width suddenly increases as the Doppler shift changes from redshift to blueshift. However, we demonstrate that this increase is caused by the superposition of two emission components. We then perform detailed analysis of the line profiles at a few selected locations on the slit. The temporal evolution of the line core is dominated by the following behavior: a rapid excursion to the blue side, accompanied by an intensity increase, followed by a linear decrease of the velocity to the red side. The maximum intensity slightly lags the maximum blueshift in Si iv, whereas the intensity enhancement slightly precedes the maximum blueshift in Mgii. We find a positive correlation between the maximum velocity and deceleration, a result that is consistent with numerical simulations of upward propagating magnetoacoustic shock waves.

We present the first model that couples the formation of the corona of a solar active region to a model of the emergence of a sunspot pair. This allows us to study when, where, and why active region loops form, and how they evolve. We use a 3D radiation MHD simulation of the emergence of an active region through the upper convection zone and the photosphere as a lower boundary for a 3D MHD coronal model. The latter accounts for the braiding of the magnetic fieldlines, which induces currents in the corona heating up the plasma. We synthesize the coronal emission for a direct comparison to observations. Starting with a basically field-free atmosphere we follow the filling of the corona with magnetic field and plasma. Numerous individually identifiable hot coronal loops form, and reach temperatures well above 1 MK with densities comparable to observations. The footpoints of these loops are found where small patches of magnetic flux concentrations move into the sunspots. The loop formation is triggered by an increase of upwards-directed Poynting flux at their footpoints in the photosphere. In the synthesized EUV emission these loops develop within a few minutes. The first EUV loop appears as a thin tube, then rises and expands significantly in the horizontal direction. Later, the spatially inhomogeneous heat input leads to a fragmented system of multiple loops or strands in a growing envelope.

We present numerical simulations of active region scale flux emergence covering a time span of up to 6 days. Flux emergence is driven by a bottom boundary condition that advects a semi-torus of magnetic field with $1.7 \times 10^{22}$ Mx flux into the computational domain. The simulations show that, even in the absense of twist, the magnetic flux is able the rise through the upper $15.5$ Mm of the convection zone and emerge into the photosphere to form spots. We find that spot formation is sensitive to the persistence of upflows at the bottom boundary footpoints, i.e. a continuing upflow would prevent spot formation. In addition, the presence of a torus-aligned flow (such flow into the retrograde direction is expected from angular momentum conservation during the rise of flux ropes through the convection zone) leads to a significant asymmetry between the pair of spots, with the spot corresponding to the leading spot on the Sun being more axisymmetric and coherent, but also forming with a delay relative to the following spot. The spot formation phase transitions directly into a decay phase. Subsurface flows fragment the magnetic field and lead to intrusions of almost field free plasma underneath the photosphere. When such intrusions reach photospheric layers, the spot fragments. The time scale for spot decay is comparable to the longest convective time scales present in the simulation domain. We find that the dispersal of flux from a simulated spot in the first two days of the decay phase is consistent with self-similar decay by turbulent diffusion.

The Interface Region Imaging Spectrograph (IRIS) small explorer spacecraft provides simultaneous spectra and images of the photosphere, chromosphere, transition region, and corona with 0.33-0.4 arcsec spatial resolution, 2 s temporal resolution and 1 km/s velocity resolution over a field-of-view of up to 175 arcsec x 175 arcsec. IRIS was launched into a Sun-synchronous orbit on 27 June 2013 using a Pegasus-XL rocket and consists of a 19-cm UV telescope that feeds a slit-based dual-bandpass imaging spectrograph. IRIS obtains spectra in passbands from 1332-1358, 1389-1407 and 2783-2834 Angstrom including bright spectral lines formed in the chromosphere (Mg II h 2803 Angstrom and Mg II k 2796 Angstrom) and transition region (C II 1334/1335 Angstrom and Si IV 1394/1403 Angstrom). Slit-jaw images in four different passbands (C II 1330, Si IV 1400, Mg II k 2796 and Mg II wing 2830 Angstrom) can be taken simultaneously with spectral rasters that sample regions up to 130 arcsec x 175 arcsec at a variety of spatial samplings (from 0.33 arcsec and up). IRIS is sensitive to emission from plasma at temperatures between 5000 K and 10 MK and will advance our understanding of the flow of mass and energy through an interface region, formed by the chromosphere and transition region, between the photosphere and corona. This highly structured and dynamic region not only acts as the conduit of all mass and energy feeding into the corona and solar wind, it also requires an order of magnitude more energy to heat than the corona and solar wind combined. The IRIS investigation includes a strong numerical modeling component based on advanced radiative-MHD codes to facilitate interpretation of observations of this complex region. Approximately eight Gbytes of data (after compression) are acquired by IRIS each day and made available for unrestricted use within a few days of the observation.

We present a method for performing data-driven simulations of solar active region formation and evolution. The approach is based on magnetofriction, which evolves the induction equation assuming the plasma velocity is proportional to the Lorentz force. The simulations of active region coronal field are driven by temporal sequences of photospheric magnetograms from the Helioseismic Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO). Under certain conditions, the data-driven simulations produce flux ropes that are ejected from the modeled active region due to loss of equilibrium. Following the ejection of flux ropes, we find an enhancement of the photospheric horizontal field near the polarity inversion line. We also present a method for the synthesis of mock coronal images based on a proxy emissivity calculated from the current density distribution in the model. This method yields mock coronal images that are somewhat reminiscent of images of active regions taken by instruments such as SDO's Atmospheric Imaging Assembly (AIA) at extreme ultraviolet wavelengths.

We investigate the importance of ambipolar diffusion and Hall currents for high-resolution comprehensive ('realistic') photospheric simulations. To do so we extended the radiative magnetohydrodynamics code \emph{MURaM} to use the generalized Ohm's law under the assumption of local thermodynamic equilibrium. We present test cases comparing analytical solutions with numerical simulations for validation of the code. Furthermore, we carried out a number of numerical experiments to investigate the impact of these neutral-ion effects in the photosphere. We find that, at the spatial resolutions currently used (5-20 km per grid point), the Hall currents and ambipolar diffusion begin to become significant -- with flows of 100 m/s in sunspot light bridges, and changes of a few percent in the thermodynamic structure of quiet-Sun magnetic features. The magnitude of the effects is expected to increase rapidly as smaller-scale variations are resolved by the simulations.

The immense volume of data generated by the suite of instruments on SDO requires new tools for efficient identifying and accessing data that is most relevant to research investigations. We have developed the Heliophysics Events Knowledgebase (HEK) to fill this need. The HEK system combines automated data mining using feature-detection methods and high-performance visualization systems for data markup. In addition, web services and clients are provided for searching the resulting metadata, reviewing results, and efficiently accessing the data. We review these components and present examples of their use with SDO data.

We present a radiative magnetohydrodynamics simulation of the formation of an Active Region on the solar surface. The simulation models the rise of a buoyant magnetic flux bundle from a depth of 7.5 Mm in the convection zone up into the solar photosphere. The rise of the magnetic plasma in the convection zone is accompanied by predominantly horizontal expansion. Such an expansion leads to a scaling relation between the plasma density and the magnetic field strength such that $B\propto\varrho^{1/2}$. The emergence of magnetic flux into the photosphere appears as a complex magnetic pattern, which results from the interaction of the rising magnetic field with the turbulent convective flows. Small-scale magnetic elements at the surface first appear, followed by their gradual coalescence into larger magnetic concentrations, which eventually results in the formation of a pair of opposite polarity spots. Although the mean flow pattern in the vicinity of the developing spots is directed radially outward, correlations between the magnetic field and velocity field fluctuations allow the spots to accumulate flux. Such correlations result from the Lorentz-force driven, counter-streaming motion of opposite-polarity fragments. The formation of the simulated Active Region is accompanied by transient light bridges between umbrae and umbral dots. Together with recent sunspot modeling, this work highlights the common magnetoconvective origin of umbral dots, light bridges and penumbral filaments.

We study the observational signature of flux emergence in the photosphere using synthetic data from a 3D MHD simulation of the emergence of a twisted flux tube. Several stages in the emergence process are considered. At every stage we compute synthetic Stokes spectra of the two iron lines Fe I 6301.5 {\AA} and Fe I 6302.5 {\AA} and degrade the data to the spatial and spectral resolution of Hinode's SOT/SP. Then, following observational practice, we apply Milne-Eddington-type inversions to the synthetic spectra in order to retrieve various atmospheric parameters and compare the results with recent Hinode observations. During the emergence sequence, the spectral lines sample different parts of the rising flux tube, revealing its twisted structure. The horizontal component of the magnetic field retrieved from the simulations is close to the observed values. The flattening of the flux tube in the photosphere is caused by radiative cooling, which slows down the ascent of the tube to the upper solar atmosphere. Consistent with the observations, the rising magnetized plasma produces a blue shift of the spectral lines during a large part of the emergence sequence.

Nonlinear force-free field (NLFFF) models are thought to be viable tools for investigating the structure, dynamics and evolution of the coronae of solar active regions. In a series of NLFFF modeling studies, we have found that NLFFF models are successful in application to analytic test cases, and relatively successful when applied to numerically constructed Sun-like test cases, but they are less successful in application to real solar data. Different NLFFF models have been found to have markedly different field line configurations and to provide widely varying estimates of the magnetic free energy in the coronal volume, when applied to solar data. NLFFF models require consistent, force-free vector magnetic boundary data. However, vector magnetogram observations sampling the photosphere, which is dynamic and contains significant Lorentz and buoyancy forces, do not satisfy this requirement, thus creating several major problems for force-free coronal modeling efforts. In this article, we discuss NLFFF modeling of NOAA Active Region 10953 using Hinode/SOT-SP, Hinode/XRT, STEREO/SECCHI-EUVI, and SOHO/MDI observations, and in the process illustrate the three such issues we judge to be critical to the success of NLFFF modeling: (1) vector magnetic field data covering larger areas are needed so that more electric currents associated with the full active regions of interest are measured, (2) the modeling algorithms need a way to accommodate the various uncertainties in the boundary data, and (3) a more realistic physical model is needed to approximate the photosphere-to-corona interface in order to better transform the forced photospheric magnetograms into adequate approximations of nearly force-free fields at the base of the corona. We make recommendations for future modeling efforts to overcome these as yet unsolved problems.

We present results from numerical modeling of emerging flux regions on the solar surface. The modeling was carried out by means of 3D radiative MHD simulations of the rise of buoyant magnetic flux tubes through the convection zone and into the photosphere. Due to the strong stratification of the convection zone, the rise results in a lateral expansion of the tube into a magnetic sheet, which acts as a reservoir for small-scale flux emergence events at the scale of granulation. The interaction of the convective downflows and the rising magnetic flux undulates it to form serpentine field lines emerging into the photosphere. Observational characteristics including the pattern of emerging flux regions, the cancellation of surface flux and associated high speed downflows, the convective collapse of photospheric flux tubes, the appearance of anomalous darkenings, the formation of bright points and the possible existence of transient kilogauss horizontal fields are discussed in the context of new observations from the Hinode Solar Optical Telescope. Implications for the local helioseismology of emerging flux regions are also discussed.

We study the emergence of magnetic flux from the near-surface layers of the solar convection zone into the photosphere. To model magnetic flux emergence, we carried out a set of numerical radiative magnetohydrodynamics simulations. Our simulations take into account the effects of compressibility, energy exchange via radiative transfer, and partial ionization in the equation of state. All these physical ingredients are essential for a proper treatment of the problem. Furthermore, the inclusion of radiative transfer allows us to directly compare the simulation results with actual observations of emerging flux. We find that the interaction between the magnetic flux tube and the external flow field has an important influence on the emergent morphology of the magnetic field. Depending on the initial properties of the flux tube (e.g. field strength, twist, entropy etc.), the emergence process can also modify the local granulation pattern. The emergence of magnetic flux tubes with a flux of $10^{19}$ Mx disturbs the granulation and leads to the transient appearance of a dark lane, which is coincident with upflowing material. These results are consistent with observed properties of emerging magnetic flux.

We study the structure and reveal the physical nature of the reversed granulation pattern in the solar photosphere by means of 3-dimensional radiative hydrodynamics simulations. We used the MURaM code to obtain a realistic model of the near-surface layers of the convection zone and the photosphere. The pattern of horizontal temperature fluctuations at the base of the photosphere consists of relatively hot granular cells bounded by the cooler intergranular downflow network. With increasing height in the photosphere, the amplitude of the temperature fluctuations diminishes. At a height of z=130-140 km in the photosphere, the pattern of horizontal temperature fluctuations reverses so that granular regions become relatively cool compared to the intergranular network. Detailed analysis of the trajectories of fluid elements through the photosphere reveal that the motion of the fluid is non-adiabatic, owing to strong radiative cooling when approaching the surface of optical depth unity followed by reheating by the radiation field from below. The temperature structure of the photosphere results from the competition between expansion of rising fluid elements and radiative heating. The former acts to lower the temperature of the fluid whereas the latter acts to increase it towards the radiative equilibrium temperature with a net entropy gain. After the fluid overturns and descends towards the convection zone, radiative energy loss again decreases the entropy of the fluid. Radiative heating and cooling of fluid elements that penetrate into the photosphere and overturn do not occur in equal amounts. The imbalance in the cumulative heating and cooling of these fluid elements is responsible for the reversal of temperature fluctuations with respect to height in the photosphere.