Magnetic flux generated and intensified by the solar dynamo emerges into the solar atmosphere, forming active regions (ARs) including sunspots. Existing theories of flux emergence suggest that the magnetic flux can rise buoyantly through the convection zone but is trapped at the photosphere, while its further rising into the atmosphere resorts to the Parker buoyancy instability. To trigger such an instability, the Lorentz force in the photosphere needs to be as large as the gas pressure gradient to hold up an extra amount of mass against gravity. This naturally results in a strongly non-force-free photosphere, which is indeed shown in typical idealized numerical simulations of flux tube buoyancy from below the photosphere into the corona. Here we conduct a statistical study of the extents of normalized Lorentz forces and torques in the emerging photospheric magnetic field with a substantially large sample of SDO/HMI vector magnetograms. We found that the photospheric field has a rather small Lorentz force and torque on average, and thus is very close to a force-free state, which is not consistent with theories as well as idealized simulations of flux emergence. Furthermore, the small extents of forces and torques seem not to be influenced by the emerging AR's size, the emergence rate, or the non-potentiality of the field. This result puts an important constraint on future development of theories and simulations of flux emergence.

Inspired by recent observations suggesting that photospheric magnetic flux cancellation occurs much more frequently than previously thought, we analytically estimated the energy released from reconnection driven by photospheric flux cancellation, and proposed that it can act as a mechanism for chromospheric and coronal heating (Priest et al., 2018). Using two-dimensional simulations we validated the analytical estimates and studied the resulting atmospheric response (Syntelis et al. 2019). In the present work, we set up three-dimensional resistive MHD simulations of two cancelling polarities in a stratified atmosphere with a horizontal external field to further validate and improve upon the analytical estimates. The computational evaluation of the parameters associated with the energy release are in good qualitative agreement with the analytical estimates. The computational Poynting energy flux into the current sheet is in good qualitative agreement with the analytical estimates, after correcting the analytical expression to better account for the horizontal extent of the current sheet. The atmospheric response to the cancellation is the formation of hot ejections, cool ejections, or a combination of both hot and cool ejections, which can appear with a time difference and/or be spatially offset, depending on the properties of the cancelling region and the resulting height of the reconnection. Therefore, during the cancellation, a wide spectrum of ejections can be formed, which can account for the variety of multi-thermal ejections associated with Ellerman bombs, UV bursts and IRIS bombs, and also other ejections associated with small-scale cancelling regions and spicules.

We report on a three-dimensional MHD numerical experiment of a small scale coronal mass ejection (CME) -like eruption propagating though a non-magnetized solar atmosphere. We find that the Kelvin-Helmholtz instability (KHI) develops at various but specific locations at the boundary layer between the erupting field and the background atmosphere, depending on the relative angle between the velocity and magnetic field. KHI develops at the front and at two of the four sides of the eruption. KHI is suppressed at the other two sides of the eruption. We also find the development of Alfv\'enic vortex shedding flows at the wake of the developing CME due to the 3D geometry of the field. Forward modeling reveals that the observational detectability of the KHI in solar eruptions is confined to a narrow $\approx10^{\circ}$ range when observing off-limb, and therefore its occurrence could be underestimated due to projection effects. The new findings can have significant implications for observations, for heating and for particle acceleration by turbulence from flow-driven instabilities associated with solar eruptions of all scales.

Context. Some of the most dynamic active regions are associated with complex photospheric magnetic configurations such as quadrupolar regions, and especially ones with a $\delta$-spot configuration and a strong Polarity Inversion Line (PIL). Aims. We study the formation and eruption of magnetic flux ropes in quadrupolar regions. Methods. 3D MHD simulations of the partial emergence of a highly twisted flux tube from the solar interior into a non-magnetized stratified atmosphere. Results. We emerge two $\Omega$-shaped loops forming a quadrupolar region. The emerging flux forms two initially separated bipoles, that later come in contact creating a $\delta$-spot central region. Above the two bipoles, two magnetic lobes expand and interact through a series of current sheets at the interface between them. Two recurrent confined eruptions are produced. In both cases, reconnection between sheared low-lying field lines forms a flux rope. Reconnection between the two lobes high in the atmosphere forms retracting field lines that push against the flux rope, creating a current sheet between them. It also forms a third magnetic lobe between the two emerged ones, that later acts as a strapping field. The flux rope eruptions are triggered when the reconnection between the flux ropes and the field above them becomes efficient enough to remove the tension of the overlying field. These reconnection events occur internally in the system. The first erupting flux rope almost fully reconnects with the overlying field. The second eruption is confined by the overlying field. During the confined eruption, the flux rope is enhanced in size, flux and twist, similar to confined-flare-to-flux-rope observations. Proxies of the emission reveal the two erupting filaments channels. A flare arcade is formed only in the second eruption due to the longer-lasting and more efficient reconnection below the flux rope.

A plethora of solar dynamic events, such as the formation of active regions, the emission of jets and the occurrence of eruptions is often associated to the emergence of magnetic flux from the interior of the Sun to the surface and above. Here, we present a short review on the onset, driving and/or triggering of such events by magnetic flux emergence. We briefly describe some key observational examples, theoretical aspects and numerical simulations, towards revealing the mechanisms that govern solar dynamics and activity related to flux emergence. We show that the combination of important physical processes like shearing and reconnection of magnetic fieldlines in emerging flux regions or at their vicinity, can power some of the most dynamic phenomena in the Sun on various temporal and spatial scales. Based on previous and recent observational and numerical studies, we highlight that, in most cases, none of these processes alone can drive and also trigger explosive phenomena releasing considerable amount of energy towards the outer solar atmosphere and space, such as flares, jets and large-scale eruptions (e.g. CMEs). In addition, one has to take into account the physical properties of the emerging field (e.g. strength, amount of flux, relative orientation to neighbouring and pre-existing magnetic fields, etc.) in order to better understand the exact role of magnetic flux emergence on the onset of solar dynamic events.

We present results of three-dimensional MHD simulations of recurrent eruptions in emerging flux regions. The initial numerical setup is the same with the work by Syntelis et al. 2017 (hereafter, Paper I). Here, we perform a parametric study on the magnetic field strength ($B_0$) of the emerging field. The kinetic energy of the produced ejective eruptions in the emerging flux region ranges from $10^{26}-10^{28}$ erg, reaching up to the energies of small Coronal Mass Ejections (CMEs). The kinetic and magnetic energies of the eruptions scale linearly in a logarithmic plot. We find that the eruptions are triggered earlier for higher $B_0$ and that $B_0$ is not directly correlated to the frequency of occurrence of the eruptions. Using large numerical domains, we show the initial stage of the partial merging of two colliding erupting fields. The partial merging occurs partly by the reconnection between the field lines of the following and the leading eruption at the interface between them. We also find that tether-cutting reconnection of the field lines of the leading eruption underneath the following eruption magnetically links the two eruptions. Shocks develop inside the leading eruption during the collision.

We report on our three-dimensional (3D) magnetohydrodynamic (MHD) simulations of cylindrical weakly twisted flux tubes emerging from 18 Mm below the photosphere. We perform a parametric study, by varying the initial magnetic field strength ($B_0$), radius ($R$), twist ($\alpha)$ and length of the emerging part of the flux tube ($\lambda$) to investigate how these parameters affect the transfer of the magnetic field from the convection zone to the photosphere. We show that the efficiency of emergence at the photosphere (i.e. how strong the photospheric field will be in comparison to $B_0$) depends not only on the $B_0$ but also the morphology of the emerging field and the twist. We show that parameters such as $B_0$ and magnetic flux cannot alone determine whether a flux tube will emerge to the solar surface. For instance, high-$B_0$ (weak-$B_0$) fields may fail (succeed) to emerge at the photosphere, depending on their geometrical properties. We also show that the photospheric magnetic field strength can vary greatly for flux tubes with the same $B_0$ but different geometric properties. Moreover, in some cases we have found scaling laws, whereby the magnetic field strength scales with the local density as $B\propto \rho^\kappa$, where $\kappa \approx 1$ deeper in the convection zone and $\kappa <1$, close to the photosphere. The transition between the two values occurs approximately when the local pressure scale ($H_p$) becomes comparable to the diameter of the flux tube ($H_p\approx2R$). We derive forms to explain how and when these scaling laws appear and compare them with the numerical simulations.

Recent observations at high spatial resolution have shown that magnetic flux cancellation occurs on the solar surface much more frequently than previously thought, and so this led Priest et al 2018 to propose magnetic reconnection driven by photospheric flux cancellation as a mechanism for chromospheric and coronal heating. In particular, they estimated analytically the amount of energy released as heat and the height of the energy release during flux cancellation. In the present work, we take the next step in the theory by setting up a two-dimensional resistive MHD simulation of two canceling polarities in the presence of a horizontal external field and a stratified atmosphere in order to check and improve upon the analytical estimates. Computational evaluation of the energy release during reconnection is found to be in good qualitative agreement with the analytical estimates. In addition, we go further and undertake an initial study of the atmospheric response to reconnection. We find that, during the cancellation, either hot ejections or cool ones or a combination of both hot and cool ejections can be formed, depending on the height of the reconnection location. The hot structures can have the density and temperature of coronal loops, while the cooler structures are suggestive of surges and large spicules.

Nanoflare models for heating the solar corona usually assume magnetic braiding and reconnection as the source of the energy. However, recent observations at record spatial resolution from the {\sc Sunrise} balloon mission suggest that photospheric magnetic flux cancellation is much more common than previously realised. We therefore examine the possibility of three-dimensional reconnection driven by flux cancellation as a cause of chromospheric and coronal heating. In particular, we estimate how the heights and amount of energy release produced by flux cancellation depend on flux size, flux cancellation speed and overlying field strength.

We report on three-dimensional (3D) Magnetohydrodynamic (MHD) simulations of recurrent eruptions in emerging flux regions. We find that reconnection of sheared fieldlines, along the polarity inversion line of an emerging bipolar region, leads to the formation of a new magnetic structure, which adopts the shape of a magnetic flux rope during its rising motion. Initially, the flux rope undergoes a slow-rise phase and, eventually, it experiences a fast-rise phase and ejective eruption towards the outer solar atmosphere. In total, four eruptions occur during the evolution of the system. For the first eruption, our analysis indicates that the torus instability initiates the eruption and that tether-cutting reconnection of the fieldlines, which envelope the flux rope, triggers the rapid acceleration of the eruptive field. For the following eruptions, we conjecture that it is the interplay between tether-cutting reconnection and torus instability, which causes the onset of the various phases. We show the 3D shape of the erupting fields, focusing more on how magnetic fieldlines reconnect during the eruptions. We find that when the envelope fieldlines reconnect mainly with themselves, hot and dense plasma is transferred closer to the core of the erupting flux rope. The same area appears to be cooler and less dense when the envelope fieldlines reconnect with neighboring sheared fieldlines. The plasma density and temperature distribution, together with the rising speeds, the energies and the size of the erupting fields indicate that they may account for small-scale (mini) Coronal Mass Ejections (CMEs).

We present a spectroscopic analysis of the pre-eruptive configuration of active region NOAA 11429, prior to two very fast coronal mass ejections (CMEs) on March 7, 2012 that are associated with this active region. We study the thermal components and the dynamics associated with the ejected flux ropes. Using differential emission measure (DEM) analysis of Hinode/EIS and SDO/AIA observations, we identify the emission components of both the flux rope and the host active region. We then follow the time evolution of the flux rope emission components by using AIA observations. The plasma density and the Doppler and non-thermal velocities associated with the flux ropes are also calculated from the EIS data. The eastern and western parts of the active region, in which the two different fast CMEs originated during two X-class flares, were studied separately. In both regions we identified an emission component in the temperature range of $\log T=6.8 - 7.1$ associated with the presence of flux ropes. The time evolution of the eastern region showed an increase in the mean DEM in this temperature range by an order of magnitude, 5 hours prior to the first CME. This was associated with a gradual rise and heating of the flux rope as manifested by blue-shifts and increased non-thermal velocities in \ion{Ca}{XV}~200.97\AA, respectively. An overall upward motion of the flux ropes was measured (relative blue-shifts of $\sim12$~\kms). The measured electron density was found to be $4\times 10^9 - 2 \times10^{10}$ cm$^{-3}$ (using the ratio of \ion{Ca}{XV}~181.90\AA\ over \ion{Ca}{XV}~200.97\AA). We compare our findings with other works on the same AR to provide a unified picture of its evolution.

Aims. We study the emergence of a non-twisted flux tube from the solar interior into the solar atmosphere. We investigate whether the length of the buoyant part of the flux tube (i.e. {\lambda}) affects the emergence of the field and the dynamics of the evolving magnetic flux system. Methods. We perform three-dimensional (3D), time-dependent, resistive, compressible MHD simulations using the Lare3D code. Results. We find that there are considerable differences in the dynamics of the emergence of a magnetic flux tube when {\lambda} is varied. In the solar interior, for larger values of {\lambda}, the rising magnetic field emerges faster and expands more due to its lower magnetic tension. As a result, its field strength decreases and its emergence above the photosphere occurs later than in the smaller {\lambda} case. However, in both cases, the emerging field at the photosphere becomes unstable in two places, forming two magnetic bipoles that interact dynamically during the evolution of the system. Most of the dynamic phenomena occur at the current layer, which is formed at the interface between the interacting bipoles. We find the formation and ejection of plasmoids, the onset of successive jets from the interface, and the impulsive heating of the plasma in the solar atmosphere. We discuss the triggering mechanism of the jets and the atmospheric response to the emergence of magnetic flux in the two cases.

We study plasma flows along selected coronal loops in NOAA Active Region 10926, observed on 3 December 2006 with Hinode's EUV Imaging Spectrograph (EIS). From the shape of the loops traced on intensity images and the Doppler shifts measured along their length we compute their three-dimensional (3D) shape and plasma flow velocity using a simple geometrical model. This calculation was performed for loops visible in the Fe VIII 185 Ang., Fe X 184 Ang., Fe XII 195 Ang., Fe XIII 202 Ang., and Fe XV 284 Ang. spectral lines. In most cases the flow is unidirectional from one footpoint to the other but there are also cases of draining motions from the top of the loops to their footpoints. Our results indicate that the same loop may show different flow patterns when observed in different spectral lines, suggesting a dynamically complex rather than a monolithic structure. We have also carried out magnetic extrapolations in the linear force-free field approximation using SOHO/MDI magnetograms, aiming toward a first-order identification of extrapolated magnetic field lines corresponding to the reconstructed loops. In all cases, the best-fit extrapolated lines exhibit left-handed twist (alpha < 0), in agreement with the dominant twist of the region.

We investigate the initiation and formation of Coronal Mass Ejections (CMEs) via detailed two-viewpoint analysis of low corona observations of a relatively fast CME acquired by the SECCHI instruments aboard the STEREO mission. The event which occurred on January 2, 2008, was chosen because of several unique characteristics. It shows upward motions for at least four hours before the flare peak. Its speed and acceleration profiles exhibit a number of inflections which seem to have a direct counterpart in the GOES light curves. We detect and measure, in 3D, loops that collapse toward the erupting channel while the CME is increasing in size and accelerates. We suggest that these collapsing loops are our first evidence of magnetic evacuation behind the forming CME flux rope. We report the detection of a hot structure which becomes the core of the white light CME. We observe and measure unidirectional flows along the erupting filament channel which may be associated with the eruption process. Finally, we compare these observations to the predictions from the standard flare-CME model and find a very satisfactory agreement. We conclude that the standard flare-CME concept is a reliable representation of the initial stages of CMEs and that multi-viewpoint, high cadence EUV observations can be extremely useful in understanding the formation of CMEs.