We carry out a study of the global three-dimensional (3D) structure of the electron density and temperature of the quiescent inner solar corona ($r<1.25 R_\odot$) by means of tomographic reconstructions and magnetohydrodynamic simulations. We use differential emission measure tomography (DEMT) and the Alfv\'en Wave Solar Model (AWSoM), in their latest versions. Two target rotations were selected from the solar minimum between solar cycles (SCs) 23 and 24 and the declining phase of SC 24. We report in quantitative detail on the 3D thermodynamic structure of the core and outer layers of the streamer belt, and of the high latitude coronal holes (CH), as revealed by the DEMT analysis. We report on the presence of two types of structures within the streamer belt, loops with temperature decreasing/increasing with height (dubbed down/up loops), as reported first in previous DEMT studies. We also estimate the heating energy flux required at the coronal base to keep these structures stable, found to be or order $10^5 erg\, cm^{-2} s^{-1}$, consistently with previous DEMT and spectroscopic studies. We discuss how these findings are consistent with coronal dissipation of Alfv\'en waves. We compare the 3D results of DEMT and AWSoM in distinct magnetic structures. We show that the agreement between the products of both techniques is the best so far, with an overall agreement $\lesssim 20\%$, depending on the target rotation and the specific coronal region. In its current implementation the ASWsoM model can not reproduce down loops though. Also, in the source region of the fast and slow components of the solar wind, the electron density of the AWSoM model increases with latitude, opposite to the trend observed in DEMT reconstructions.

Coronal mass ejections (CMEs) on stars other than the Sun have proven very difficult to detect. One promising pathway lies in the detection of type II radio bursts. Their appearance and distinctive properties are associated with the development of an outward propagating CME-driven shock. However, dedicated radio searches have not been able to identify these transient features in other stars. Large Alfv\'en speeds and the magnetic suppression of CMEs in active stars have been proposed to render stellar eruptions "radio-quiet". Employing 3D magnetohydrodynamic simulations, we study here the distribution of the coronal Alfv\'en speed, focusing on two cases representative of a young Sun-like star and a mid-activity M-dwarf (Proxima Centauri). These results are compared with a standard solar simulation and used to characterize the shock-prone regions in the stellar corona and wind. Furthermore, using a flux-rope eruption model, we drive realistic CME events within our M-dwarf simulation. We consider eruptions with different energies to probe the regimes of weak and partial CME magnetic confinement. While these CMEs are able to generate shocks in the corona, those are pushed much farther out compared to their solar counterparts. This drastically reduces the resulting type II radio burst frequencies down to the ionospheric cutoff, which impedes their detection with ground-based instrumentation.

We conduct a post hoc analysis of solar flare predictions made by a Long Short Term Memory (LSTM) model employing data in the form of Space-weather HMI Active Region Patches (SHARP) parameters calculated from data in proximity to the magnetic polarity inversion line where the flares originate. We train the the LSTM model for binary classification to provide a prediction score for the probability of M/X class flares to occur in next hour. We then develop a dimension-reduction technique to reduce the dimensions of SHARP parameter (LSTM inputs) and demonstrate the different patterns of SHARP parameters corresponding to the transition from low to high prediction score. Our work shows that a subset of SHARP parameters contain the key signals that strong solar flare eruptions are imminent. The dynamics of these parameters have a highly uniform trajectory for many events whose LSTM prediction scores for M/X class flares transition from very low to very high. The results demonstrate the existence of a few threshold values of SHARP parameters that when surpassed indicate a high probability of the eruption of a strong flare. Our method has distilled the knowledge of solar flare eruption learnt by deep learning model and provides a more interpretable approximation, which provides physical insight to processes driving solar flares.

We develop a mixed Long Short Term Memory (LSTM) regression model to predict the maximum solar flare intensity within a 24-hour time window 0$\sim$24, 6$\sim$30, 12$\sim$36 and 24$\sim$48 hours ahead of time using 6, 12, 24 and 48 hours of data (predictors) for each Helioseismic and Magnetic Imager (HMI) Active Region Patch (HARP). The model makes use of (1) the Space-weather HMI Active Region Patch (SHARP) parameters as predictors and (2) the exact flare intensities instead of class labels recorded in the Geostationary Operational Environmental Satellites (GOES) data set, which serves as the source of the response variables. Compared to solar flare classification, the model offers us more detailed information about the exact maximum flux level, i.e. intensity, for each occurrence of a flare. We also consider classification models built on top of the regression model and obtain better results in solar flare classifications. Our results suggest that the most efficient time period for predicting the solar activity is within 24 hours before the prediction time using the SHARP parameters and the LSTM model.

A deep learning network, Long-Short Term Memory (LSTM) network, is used in this work to predict whether the maximum flare class an active region (AR) will produce in the next 24 hours is class $\Gamma$. We considered $\Gamma$ are $\ge M$, $\ge C$ and any flare class. The essence of using LSTM, which is a recurrent neural network, is its capability to capture temporal information of the data samples. The input features are time sequences of 20 magnetic parameters from SHARPs - Space-weather HMI Active Region Patches. We analyzed active regions from June 2010 to Dec 2018, using the Geostationary Operational Environmental Satellite (GOES) X-ray flare catalogs and label the data samples with identified ARs in the GOES X-ray flare catalogs. Our results (i) shows consistent skill scores with recently published results using LSTMs and better than the previous work using single time input (eg. DeFN) (ii) The skill scores from the model show essential differences when different years of data was chosen for training and testing.

We perform a validation study of the latest version of the Alfv\'{e}n Wave Solar atmosphere Model (AWSoM) within the Space Weather Modeling Framework (SWMF). To do so, we compare the simulation results of the model with a comprehensive suite of observations for Carrington rotations representative of the solar minimum conditions extending from the solar corona to the heliosphere up to the Earth. In the low corona ($r < 1.25$ \Rs), we compare with EUV images from both STEREO-A/EUVI and SDO/AIA and to three-dimensional (3-D) tomographic reconstructions of the electron temperature and density based on these same data. We also compare the model to tomographic reconstructions of the electron density from SOHO/LASCO observations ($2.55 < r < 6.0$\Rs). In the heliosphere, we compare model predictions of solar wind speed with velocity reconstructions from InterPlanetary Scintillation (IPS) observations. For comparison with observations near the Earth, we use OMNI data. Our results show that the improved AWSoM model performs well in quantitative agreement with the observations between the inner corona and 1 AU. The model now reproduces the fast solar wind speed in the polar regions. Near the Earth, our model shows good agreement with observations of solar wind velocity, proton temperature and density. AWSoM offers an extensive application to study the solar corona and larger heliosphere in concert with current and future solar missions as well as being well suited for space weather predictions.

We present several methods towards construction of precursors, which show great promise towards early predictions, of solar flare events in this paper. A data pre-processing pipeline is built to extract useful data from multiple sources, Geostationary Operational Environmental Satellites (GOES) and Solar Dynamics Observatory (SDO)/Helioseismic and Magnetic Imager (HMI), to prepare inputs for machine learning algorithms. Two classification models are presented: classification of flares from quiet times for active regions and classification of strong versus weak flare events. We adopt deep learning algorithms to capture both the spatial and temporal information from HMI magnetogram data. Effective feature extraction and feature selection with raw magnetogram data using deep learning and statistical algorithms enable us to train classification models to achieve almost as good performance as using active region parameters provided in HMI/Space-Weather HMI-Active Region Patch (SHARP) data files. Case studies show a significant increase in the prediction score around 20 hours before strong solar flare events.

Ground-based solar observations provide key contextual data (i.e., the 'big picture') to produce a complete description of the only astrosphere we can study in situ: our Sun's heliosphere. The next decade will see the beginning of operations of the Daniel K. Inouye Solar Telescope (DKIST). DKIST will join NASA's Parker Solar Probe and the NASA/ESA Solar Orbital mission, which together will study our Sun's atmosphere with unprecedented detail. This white paper outlines the current paradigm for ground-based solar synoptic observations, and indicates those areas that will benefit from focused attention.

We examine Alfv\'en Wave Solar atmosphere Model (AWSoM) predictions of the first Parker Solar Probe (PSP) encounter. We focus on the 12-day closest approach centered on the 1st perihelion. AWSoM (van der Holst et al., 2014) allows us to interpret the PSP data in the context of coronal heating via Alfv\'en wave turbulence. The coronal heating and acceleration is addressed via outward-propagating low-frequency Alfv\'en waves that are partially reflected by Alfv\'en speed gradients. The nonlinear interaction of these counter-propagating waves results in a turbulent energy cascade. To apportion the wave dissipation to the electron and anisotropic proton temperatures, we employ the results of the theories of linear wave damping and nonlinear stochastic heating as described by Chandran et al. (2011). We find that during the first encounter, PSP was in close proximity to the heliospheric current sheet (HCS) and in the slow wind. PSP crossed the HCS two times, namely at 2018/11/03 UT 01:02 and 2018/11/08 UT 19:09 with perihelion occuring on the south of side of the HCS. We predict the plasma state along the PSP trajectory, which shows a dominant proton parallel temperature causing the plasma to be firehose unstable.

Recent detections of high-energy $\gamma$-rays from behind-the-limb (BTL) solar flares by the \emph{Fermi $\gamma$-ray Space Telescope} pose a puzzle and challenge on the particle acceleration and transport mechanisms. In such events, the $\gamma$-ray emission region is located away from the BTL flare site by up to tens of degrees in heliogrpahic longitude. It is thus hypothesized that particles are accelerated at the shock driven by the coronal mass ejection (CME) and then travel from the shock downstream back to the front side of the Sun to produce the observed $\gamma$-rays. To test this scenario, we performed data-driven, global magnetohydrodynamics simulations of the CME associated with a well-observed BTL flare on 2014 September 1. We found that part of the CME-driven shock develops magnetic connectivity with the $\gamma$-ray emission region, facilitating transport of particles back to the Sun. Moreover, the observed increase in $\gamma$-ray flux is temporally correlated with (1) the increase of the shock compression ratio and (2) the presence of a quasi-perpendicular shock over the area that is magnetically connected to the $\gamma$-ray emitting region, both conditions favoring the diffusive shock acceleration (DSA) of particles. These results support the above hypothesis and can help resolve another puzzle, i.e., long-duration (up to 20 hours) $\gamma$-rays flares. We suggest that, in addition to DSA, stochastic acceleration by plasma turbulence may also play a role, especially in the shock downstream region and during the early stage when the shock Alfv\'{e}n Mach number is small.

The history and present state of large-scale magnetohydrodynamic (MHD) modeling of the solar corona and the solar wind with steady or quasi-steady coronal physics is reviewed. We put the evolution of ideas leading to the recognition of the existence of an expanding solar atmosphere into historical context. The development and main features of the first generation of global corona and solar wind models are described in detail. This historical perspective is also applied to the present suite of global corona and solar wind models. We discuss the evolution of new ideas and their implementation into numerical simulation codes. We point out the scientific and computational challenges facing these models and discuss the ways various groups tried to overcome these challenges. Next, we discuss the latest, state-of-the art models and point to the expected next steps in modeling the corona and the interplanetary medium.

Coronal mass ejections (CMEs) were discovered in the early 1970s when space-borne coronagraphs revealed that eruptions of plasma are ejected from the Sun. Today, it is known that the Sun produces eruptive flares, filament eruptions, coronal mass ejections and failed eruptions; all thought to be due to a release of energy stored in the coronal magnetic field during its drastic reconfiguration. This review discusses the observations and physical mechanisms behind this eruptive activity, with a view to making an assessment of the current capability of forecasting these events for space weather risk and impact mitigation. Whilst a wealth of observations exist, and detailed models have been developed, there still exists a need to draw these approaches together. In particular more realistic models are encouraged in order to asses the full range of complexity of the solar atmosphere and the criteria for which an eruption is formed. From the observational side, a more detailed understanding of the role of photospheric flows and reconnection is needed in order to identify the evolutionary path that ultimately means a magnetic structure will erupt.

We present an updated global model of the solar corona, including the transition region. We simulate the realistic tree-dimensional (3D) magnetic field using the data from the photospheric magnetic field measurements and assume the magnetohydrodynamic (MHD) Alfv\'en wave turbulence and its non-linear dissipation to be the only source for heating the coronal plasma and driving the solar wind. In closed field regions the dissipation efficiency in a balanced turbulence is enhanced. In the coronal holes we account for a reflection of the outward propagating waves, which is accompanied by generation of weaker counter-propagating waves. The non-linear cascade rate degrades in strongly imbalanced turbulence, thus resulting in colder coronal holes. The distinctive feature of the presented model is the description of the low corona as almost-steady-state low-beta plasma motion and heat flux transfer along the magnetic field lines. We trace the magnetic field lines through each grid point of the lower boundary of the global corona model, chosen at some heliocentric distance, $R=R_{b}\sim1.1\ R_\odot$ well above the transition region. One can readily solve the plasma parameters along the magnetic field line from 1D equations for the plasma motion and heat transport together with the Alfv\'en wave propagation, which adequately describe physics within the heliocentric distances range, $R_{\odot}<R<R_{b}$, in the low solar corona. By interfacing this threaded-field-lines model with the full MHD global corona model at $r=R_{b}$, we find the global solution and achieve a faster-than-real-time performance of the model on $\sim200$ cores.

We present a first-principles-based coronal mass ejection (CME) model suitable for both scientific and operational purposes by combining a global magnetohydrodynamics (MHD) solar wind model with a flux rope-driven CME model. Realistic CME events are simulated self-consistently with high fidelity and forecasting capability by constraining initial flux rope parameters with observational data from GONG, SOHO/LASCO, and STEREO/COR. We automate this process so that minimum manual intervention is required in specifying the CME initial state. With the newly developed data-driven Eruptive Event Generator Gibson-Low (EEGGL), we present a method to derive Gibson-Low (GL) flux rope parameters through a handful of observational quantities so that the modeled CMEs can propagate with the desired CME speeds near the Sun. A test result with CMEs launched with different Carrington rotation magnetograms are shown. Our study shows a promising result for using the first-principles-based MHD global model as a forecasting tool, which is capable of predicting the CME direction of propagation, arrival time, and ICME magnetic field at 1 AU (see companion paper by Jin et al. 2016b).

We perform and analyze results of a global magnetohydrodyanmic (MHD) simulation of the fast coronal mass ejection (CME) that occurred on 2011 March 7. The simulation is made using the newly developed Alfv\'en Wave Solar Model (AWSoM), which describes the background solar wind starting from the upper chromosphere and extends to 24 R$_{\odot}$. Coupling AWSoM to an inner heliosphere (IH) model with the Space Weather Modeling Framework (SWMF) extends the total domain beyond the orbit of Earth. Physical processes included in the model are multi-species thermodynamics, electron heat conduction (both collisional and collisionless formulations), optically thin radiative cooling, and Alfv\'en-wave turbulence that accelerates and heats the solar wind. The Alfv\'en-wave description is physically self-consistent, including non-Wentzel-Kramers-Brillouin (WKB) reflection and physics-based apportioning of turbulent dissipative heating to both electrons and protons. Within this model, we initiate the CME by using the Gibson-Low (GL) analytical flux rope model and follow its evolution for days, in which time it propagates beyond STEREO A. A detailed comparison study is performed using remote as well as \textit{in situ} observations. Although the flux rope structure is not compared directly due to lack of relevant ejecta observation at 1 AU in this event, our results show that the new model can reproduce many of the observed features near the Sun (e.g., CME-driven extreme ultraviolet (EUV) waves, deflection of the flux rope from the coronal hole, "double-front" in the white light images) and in the heliosphere (e.g., shock propagation direction, shock properties at STEREO A).

The higher charge states found in slow ($<$400km s$^{-1}$) solar wind streams compared to fast streams have supported the hypothesis that the slow wind originates in closed coronal loops, and released intermittently through reconnection. Here we examine whether a highly ionized slow wind can also form along steady and open magnetic field lines. We model the steady-state solar atmosphere using AWSoM, a global magnetohydrodynamic model driven by Alfv{\'e}n waves, and apply an ionization code to calculate the charge state evolution along modeled open field lines. This constitutes the first charge states calculation covering all latitudes in a realistic magnetic field. The ratios $O^{+7}/O^{+6}$ and $C^{+6}/C^{+5}$ are compared to in-situ Ulysses observations, and are found to be higher in the slow wind, as observed; however, they are under-predicted in both wind types. The modeled ion fractions of S, Si, and Fe are used to calculate line-of-sight intensities, which are compared to EIS observations above a coronal hole. The agreement is partial, and suggests that all ionization rates are under-predicted. Assuming the presence of suprathermal electrons improved the agreement with both EIS and Ulysses observations; importantly, the trend of higher ionization in the slow wind was maintained. The results suggest there can be a sub-class of slow wind that is steady and highly ionized. Further analysis shows it originates from coronal hole boundaries (CHB), where the modeled electron density and temperature are higher than inside the hole, leading to faster ionization. This property of CHBs is global, and observationally supported by EUV tomography.

We present a new version of the Alfven Wave Solar Model (AWSoM), a global model from the upper chromosphere to the corona and the heliosphere. The coronal heating and solar wind acceleration are addressed with low-frequency Alfven wave turbulence. The injection of Alfven wave energy at the inner boundary is such that the Poynting flux is proportional to the magnetic field strength. The three-dimensional magnetic field topology is simulated using data from photospheric magnetic field measurements. This model does not impose open-closed magnetic field boundaries; those develop self-consistently. The physics includes: (1) The model employs three different temperatures, namely the isotropic electron temperature and the parallel and perpendicular ion temperatures. The firehose, mirror, and ion-cyclotron instabilities due to the developing ion temperature anisotropy are accounted for. (2) The Alfven waves are partially reflected by the Alfven speed gradient and the vorticity along the field lines. The resulting counter-propagating waves are responsible for the nonlinear turbulent cascade. The balanced turbulence due to uncorrelated waves near the apex of the closed field lines and the resulting elevated temperatures are addressed. (3) To apportion the wave dissipation to the three temperatures, we employ the results of the theories of linear wave damping and nonlinear stochastic heating. (4) We have incorporated the collisional and collisionless electron heat conduction. We compare the simulated multi-wavelength EUV images of CR2107 with the observations from STEREO/EUVI and SDO/AIA instruments. We demonstrate that the reflection due to strong magnetic fields in proximity of active regions intensifies the dissipation and observable emission sufficiently.

We report on a numerical investigation of two coronal mass ejections (CMEs) which interact as they propagate in the inner heliosphere. We focus on the effect of the orientation of the CMEs relative to each other by performing four different simulations with the axis of the second CME rotated by 90 degrees from one simulation to the next. Each magneto-hydrodynamic (MHD) simulation is performed in three dimensions (3-D) with the Space Weather Modeling Framework (SWMF) in an idealized setting reminiscent of solar minimum conditions. We extract synthetic satellite measurements during and after the interaction and compare the different cases. We also analyze the kinematics of the two CMEs, including the evolution of their widths and aspect ratios. We find that the first CME contracts radially as a result of the interaction in all cases, but the amount of subsequent radial expansion depends on the relative orientation of the two CMEs. Reconnection between the two ejecta and between the ejecta and the interplanetary magnetic field (IMF) determines the type of structure resulting from the interaction. When a CME with a high inclination with respect to the ecliptic overtakes one with a low inclination, it is possible to create a compound event with a smooth rotation in the magnetic field vector over more than 180 degrees. Due to reconnection, the second CME only appears as an extended "tail", and the event may be mistaken for a glancing encounter with an isolated CME. This configuration differs significantly from the one usually studied of a multiple-magnetic cloud event, which we found to be associated with the interaction of two CMEs with the same orientation.

We present a new global model of the solar corona, including the low corona, the transition region and the top of chromosphere. The realistic 3D magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating the MHD Alfven wave turbulence. We assume this turbulence and its non-linear dissipation to be the only momentum and energy source for heating the coronal plasma and driving the solar wind. The difference between the turbulence dissipation efficiency in coronal holes and that in closed field regions is because the non-linear cascade rate degrades in strongly anisotropic turbulence in coronal holes (no inward propagating wave), thus resulting in colder coronal holes with the bi-modal solar wind originating from them. The detailed presentation of the theoretical model is illustrated with the synthetic images for multi-wavelength EUV emission compared with the observations from SDO AIA and Stereo EUVI instruments for the Carrington rotation 2107.

Magnetic loops are building blocks of the closed-field corona. While active region loops are readily seen in images taken at EUV and X-ray wavelengths, quiet Sun loops are seldom identifiable and therefore difficult to study on an individual basis. The first analysis of solar minimum (Carrington Rotation 2077) quiet Sun (QS) coronal loops utilizing a novel technique called the Michigan Loop Diagnostic Technique (MLDT) is presented. This technique combines Differential Emission Measure Tomography (DEMT) and a potential field source surface (PFSS) model, and consists of tracing PFSS field lines through the tomographic grid on which the Local Differential Emission Measure (LDEM) is determined. As a result, the electron temperature Te and density Ne at each point along each individual field line can be obtained. Using data from STEREO/EUVI and SOHO/MDI, the MLDT identifies two types of QS loops in the corona: so-called "up" loops in which the temperature increases with height, and so-called "down" loops in which the temperature decreases with height. Up loops are expected, however, down loops are a surprise, and furthermore, they are ubiquitous in the low-latitude corona. Up loops dominate the QS at higher latitudes. The MLDT allows independent determination of the empirical pressure and density scale heights, and the differences between the two remain to be explained. The down loops appear to be a newly discovered property of the solar minimum corona that may shed light on the physics of coronal heating. The results are shown to be robust to the calibration uncertainties of the EUVI instrument.

We examine a simulation of flux emergence and cancellation, which shows a complex sequence of processes that accumulate free magnetic energy in the solar corona essential for the eruptive events such as coronal mass ejections (CMEs), filament eruptions and flares. The flow velocity at the surface and in the corona shows a consistent shearing pattern along the polarity inversion line (PIL), which together with the rotation of the magnetic polarities, builds up the magnetic shear. Tether-cutting reconnection above the PIL then produces longer sheared magnetic field lines that extend higher into the corona, where a sigmoidal structure forms. Most significantly, reconnection and upward energy-flux transfer are found to occur even as magnetic flux is submerging and appears to cancel at the photosphere. A comparison of the simulated coronal field with the corresponding coronal potential field graphically shows the development of nonpotential fields during the emergence of the magnetic flux and formation of sunspots.

We simulate the buoyant rise of a magnetic flux rope from the solar convection zone into the corona to better understand the energetic coupling of the solar interior to the corona. The magnetohydrodynamic model addresses the physics of radiative cooling, coronal heating and ionization, which allow us to produce a more realistic model of the solar atmosphere. The simulation illustrates the process by which magnetic flux emerges at the photosphere and coalesces to form two large concentrations of opposite polarities. We find that the large-scale convective motion in the convection zone is critical to form and maintain sunspots, while the horizontal converging flows in the near surface layer prevent the concentrated polarities from separating. The foot points of the sunspots in the convection zone exhibit a coherent rotation motion, resulting in the increasing helicity of the coronal field. Here, the local configuration of the convection causes the convergence of opposite polarities of magnetic flux with a shearing flow along the polarity inversion line. During the rising of the flux rope, the magnetic energy is first injected through the photosphere by the emergence, followed by energy transport by horizontal flows, after which the energy is subducted back to the convection zone by the submerging flows.

A three dimensional (3D) tomographic reconstruction of the local differential emission measure (LDEM) of the global solar corona during the whole heliosphere interval (WHI, Carrington rotation CR-2068) is presented, based on STEREO/EUVI images. We determine the 3D distribution of the electron density, mean temperature, and temperature spread, in the range of heliocentric heights 1.03 to 1.23 Rsun. The reconstruction is complemented with a potential field source surface (PFSS) magnetic-field model. The streamer core, streamer legs, and subpolar regions are analyzed and compared to a similar analysis previously performed for CR-2077, very near the absolute minimum of the Solar Cycle 23. In each region, the typical values of density and temperature are similar in both periods. The WHI corona exhibits a streamer structure of relatively smaller volume and latitudinal extension than during CR-2077, with a global closed-to-open density contrast about 6% lower, and a somewhat more complex morphology. The average basal electron density is found to be about 2.23 and 1.08 x 10^8 cm^-3, in the streamer core and subpolar regions, respectively. The electron temperature is quite uniform over the analyzed height range, with average values of about 1.13 and 0.93 MK, in the streamer core and subpolar regions, respectively. Within the streamer closed region, both periods show higher temperatures at mid-latitudes and lower temperatures near the equator. Both periods show beta>1 in the streamer core and beta<1 in the surrounding open regions, with CR-2077 exhibiting a stronger contrast. Hydrostatic fits to the electron density are performed, and the scale height is compared to the LDEM mean electron temperature. Within the streamer core, the results are consistent with an isothermal hydrostatic plasma regime, with the temperatures of ions and electrons differing by up to about 10% .. (continues)..

On 13 February 2009, a coronal wave -- CME -- dimming event was observed in quadrature by the STEREO spacecraft. We analyze this event using a three-dimensional, global magnetohydrodynamic (MHD) model for the solar corona. The numerical simulation is driven and constrained by the observations, and indicates where magnetic reconnection occurs between the expanding CME core and surrounding environment. We focus primarily on the lower corona, extending out to $3R_{\odot}$; this range allows simultaneous comparison with both EUVI and COR1 data. Our simulation produces a diffuse coronal bright front remarkably similar to that observed by STEREO/EUVI at 195 \AA. It is made up of \emph{two} components, and is the result of a combination of both wave and non-wave mechanisms. The CME becomes large-scale quite low ($<$ 200 Mm) in the corona. It is not, however, an inherently large-scale event; rather, the expansion is facilitated by magnetic reconnection between the expanding CME core and the surrounding magnetic environment. In support of this, we also find numerous secondary dimmings, many far from the initial CME source region. Relating such dimmings to reconnecting field lines within the simulation provides further evidence that CME expansion leads to the "opening" of coronal field lines on a global scale. Throughout the CME expansion, the coronal wave maps directly to the CME footprint. Our results suggest that the ongoing debate over the "true" nature of diffuse coronal waves may be mischaracterized. It appears that \emph{both} wave and non-wave models are required to explain the observations and understand the complex nature of these events.

Since the first discovery of an extrasolar planetary system more than a decade ago, hundreds more have been discovered. Surprisingly, many of these systems harbor Jupiter-class gas giants located close to the central star, at distances of 0.1 AU or less. Observations of chromospheric 'hot spots' that rotate in phase with the planetary orbit, and elevated stellar X-ray luminosities,suggest that these close-in planets significantly affect the structure of the outer atmosphere of the star through interactions between the stellar magnetic field and the planetary magnetosphere. Here we carry out the first detailed three-dimensional MagnetoHydroHynamics (MHD) simulation containing the two magnetic bodies and explore the consequences of such interactions on the steady-state coronal structure. The simulations reproduce the observable features of 1) increase in the total X-ray luminosity, 2) appearance of coronal hot spots, and 3) phase shift of these spots with respect to the direction of the planet. The proximate cause of these is an increase in the density of coronal plasma in the direction of the planet, which prevents the corona from expanding and leaking away this plasma via a stellar wind. The simulations produce significant low temperature heating. By including dynamical effects, such as the planetary orbital motion, the simulation should better reproduce the observed coronal heating.

This study examines the tail disconnection event on April 20, 2007 on comet 2P/Encke, caused by a coronal mass ejection (CME) at a heliocentric distance of 0.34 AU. During their interaction, both the CME and the comet are visible with high temporal and spatial resolution by the STEREO-A spacecraft. Previously, only current sheets or shocks have been accepted as possible reasons for comet tail disconnections, so it is puzzling that the CME caused this event. The MHD simulation presented in this work reproduces the interaction process and demonstrates how the CME triggered a tail disconnection in the April 20 event. It is found that the CME disturbs the comet with a combination of a $180^\circ$ sudden rotation of the interplanetary magnetic field (IMF), followed by a $90^\circ$ gradual rotation. Such an interpretation applies our understanding of solar wind-comet interactions to determine the \textit{in situ} IMF orientation of the CME encountering Encke.

Shear flows have been prescribed in numerical models of coronal mass ejections and flares for decades as a way of energizing magnetic fields to erupt. While such shear flows have long been observed in the solar atmosphere, until recently, there was no compelling physical explanation for them. This paper will discuss the discovery that such shear flows are readily explained as a response to the Lorentz force that naturally occurs as bipolar magnetic fields emerge and expand in a gravitationally stratified atmosphere. It will be shown that shearing motions transport axial flux, and magnetic energy from the submerged portion of the field to the expanding portion, strongly coupling the solar interior to the corona. This physical process explains active region shear flows and why the magnetic field is found to be nearly parallel to photospheric polarity inversion lines where prominences form. Finally, shear flows driven by the Lorentz force are shown to produce a loss of equilibrium and eruption in magnetic arcades and flux ropes offering a convincing explanation for CMEs and flares.

We numerically model the coronal mass ejection (CME) event of October 28, 2003 that erupted from active region 10486 and propagated to Earth in less than 20 hours causing severe geomagnetic storms. The magnetohydrodynamic (MHD) model is formulated by first arriving at a steady state corona and solar wind employing synoptic magnetograms. We initiate two CMEs from the same active region, one approximately a day earlier that preconditions the solar wind for the much faster CME on the 28th. This second CME travels through the corona at a rate of over 2500 km s$^{-1}$ driving a strong forward shock. We clearly identify this shock in an image produced by the Large Angle Spectrometric Coronagraph (LASCO) C3, and reproduce the shock and its appearance in synthetic white light images from the simulation. We find excellent agreement with both the general morphology and the quantitative brightness of the model CME with LASCO observations. These results demonstrate that the CME shape is largely determined by its interaction with the ambient solar wind and may not be sensitive to the initiation process. We then show how the CME would appear as observed by wide-angle coronagraphs onboard the Solar Terrestrial Relations Observatory (STEREO) spacecraft. We find complex time evolution of the white-light images as a result of the way in which the density structures pass through the Thomson sphere. The simulation is performed with the Space Weather Modeling Framework (SWMF).

We describe a method to measure the magnetic field orientation of coronal mass ejections (CMEs) using Faraday rotation (FR). Two basic FR profiles, Gaussian-shaped with a single polarity or "N"-like with polarity reversals, are produced by a radio source occulted by a moving flux rope depending on its orientation. These curves are consistent with the Helios observations, providing evidence for the flux-rope geometry of CMEs. Many background radio sources can map CMEs in FR onto the sky. We demonstrate with a simple flux rope that the magnetic field orientation and helicity of the flux rope can be determined 2-3 days before it reaches Earth, which is of crucial importance for space weather forecasting. An FR calculation based on global magnetohydrodynamic (MHD) simulations of CMEs in a background heliosphere shows that FR mapping can also resolve a CME geometry curved back to the Sun. We discuss implementation of the method using data from the Mileura Widefield Array (MWA).

To model the interaction between the solar wind and the interstellar wind, magnetic fields must be included. Recently Opher et al. 2003 found that, by including the solar magnetic field in a 3D high resolution simulation using the University of Michigan BATS-R-US code, a jet-sheet structure forms beyond the solar wind Termination Shock. Here we present an even higher resolution three-dimensional case where the jet extends for $150AU$ beyond the Termination Shock. We discuss the formation of the jet due to a de Laval nozzle effect and it's su bsequent large period oscillation due to magnetohydrodynamic instabilities. To verify the source of the instability, we also perform a simplified two dimensional-geometry magnetohydrodynamic calculation of a plane fluid jet embedded in a neutral sheet with the profiles taken from our 3D simulation. We find remarkable agreement with the full three-dimensional evolution. We compare both simulations and the temporal evolution of the jet showing that the sinuous mode is the dominant mode that develops into a velocity-shear-instability with a growth rate of $5 \times 10^{-9} sec^{-1}=0.027 years^{-1}$. As a result, the outer edge of the heliosphere presents remarkable dynamics, such as turbulent flows caused by the motion of the jet. Further study, e.g., including neutrals and the tilt of the solar rotation from the magnetic axis, is required before we can definitively address how this outer boundary behaves. Already, however, we can say that the magnetic field effects are a major player in this region changing our previous notion of how the solar system ends.

There is currently a controversy as to whether Voyager 1 has already crossed the Termination Shock, the first boundary of the Heliosphere. The region between the Termination Shock and the Heliopause, the Helisheath, is one of the most unknown regions theoretically. In the Heliosheath magnetic effects are crucial, as the solar magnetic field is compressed at the Termination Shock by the slowing flow. Recently, our simulations showed that the Heliosheath presents remarkable dynamics, with turbulent flows and the presence of a jet flow at the current sheet that is unstable due to magnetohydrodynamic instabilities \cite{opher,opher1}. In this paper we review these recent results, and present an additional simulation with constant neutral atom background. In this case the jet is still present but with reduced intensity. Further study, e.g., including neutrals and the tilt of the solar rotation from the magnetic axis, is required before we can definitively address how the Heliosheath behaves. Already we can say that this region presents remarkable dynamics, with turbulent flows, indicating that the Heliosheath might be very different from what we previously thought.

The Voyager spacecraft is now approaching the edge of the solar system. Near the boundary between the solar system and the interstellar medium we find that an unstable ``jet-sheet'' forms. The jet-sheet oscillates up and down due to a velocity shear instability. This result is due to a novel application of a state-of-art 3D Magnetohydrodynamic (MHD) code with a highly refined grid. We assume as a first approximation that the solar magnetic and rotation axes are aligned. The effect of a tilt of the magnetic axis with respect to the rotation axis remains to be seen. We include in the model self-consistently magnetic field effects in the interaction between the solar and interstellar winds. Previous studies of this interaction had poorer spatial resolution and did not include the solar magnetic field. This instability can affect the entry of energetic particles into the solar system and the intermixing of solar and interstellar material. The same effect found here is predicted for the interaction of rotating magnetized stars possessing supersonic winds and moving with respect to the interstellar medium, such as O stars.