Solar eruptions are the most powerful drivers of space weather. To understand their cause and nature, it is crucial to know how the coronal magnetic field evolves before eruption. Here we study the formation process of a relatively large-scale magnetic flux rope (MFR) in active region NOAA~12371 that erupts with a major flare and coronal mass ejection on 2015 June 21. A data-driven numerical magnetohydrodynamic model is employed to simulate three-dimensional coronal magnetic field evolution of one-day duration before the eruption. Comparison between the observed features and our modeled magnetic field discloses how the pre-eruption MFR forms. Initially, the magnetic field lines were weakly twisted as being simple sheared arcades. Then a long MFR was formed along the polarity inversion line due to the complex photospheric motion, which is mainly shearing rather than twisting. The presence of the MFR is evidenced by a coherent set of magnetic field lines with twist number above unity. Below the MFR a current sheet is shown in the model, suggesting that tether-cutting reconnection plays a key role in the MFR formation. The MFR's flux grows as more and more field lines are twisted due to continuous injection of magnetic helicity by the photospheric motions. Meanwhile, the height of the MFR's axis increases monotonely from its formation. By an analysis of the decay index of its overlying field, we suggest that it is because the MFR runs into the torus instability regime and becomes unstable that finally triggers the eruption.

Two X-class solar flares occurred on 2017 September 6 from active region NOAA 12673: the first one is a confined X2.2 flare, and it is followed only $\sim 3$ hours later by the second one, which is the strongest flare in solar cycle 24, reaching X9.3 class and accompanied with a coronal mass ejection. Why these two X-class flares occurred in the same position with similar magnetic configurations, but one is eruptive while the other is not? Here we track the coronal magnetic field evolution via nonlinear force-free field extrapolations from a time sequence of vector magnetograms with high cadence. A detailed analysis of the magnetic field shows that a magnetic flux rope (MFR) forms and grows gradually before the first flare, and shortly afterwards, the MFR's growth is significantly enhanced with a much faster rise in height, from far below the threshold of torus instability to above it, while the magnetic twist only increases mildly. Combining EUV observations and the magnetic field extrapolation, we found that overlying the MFR is a null-point magnetic topology, where recurrent brightening is seen after the first flare. We thus suggest a scenario to interpret the occurrence of the two flares. The first flare occurred since the MFR reached a high enough height to activate the null point, and its continuous expansion forces the null-point reconnection recurrently. Such reconnection weakens the overlying field, allowing the MFR to rise faster, which eventually crosses the threshold of torus instability and triggers the second, eruptive flare.

Soar filament eruptions are often associated with solar flares and coronal mass ejections (CMEs), which are the major impacts on space weather. However, the fine structures and the trigger mechanisms of solar filaments are still unclear. To address these issues, we studied a failed solar active-region filament eruption associated with a C-class flare by using high-resolution H$\alpha$ images from the New Vacuum Solar Telescope (NVST), supplemented by EUV observations of the Solar Dynamical Observatory (SDO). Before the filament eruption, a small bi-pole magnetic field emerged below the filament. And then magnetic reconnection between the filament and the emerging bi-pole magnetic field triggered the filament eruption. During the filament eruption, the untwisting motion of the filament can be clearly traced by the eruptive threads. Moreover, the foot-points of the eruptive threads are determined by tracing the descending filament materials. Note that the filament twisted structure and the right part of the eruptive filament threads cannot be seen before the filament eruption. These eruptive threads at the right part of the filament are found to be rooting in the weak negative polarities near the main negative sunspot. Moreover, a new filament formed in the filament channel due to material injection from the eruptive filament. The above observations and the potential field extrapolations are inclined to support that the filament materials were transferred into the overlying magnetic loops and the nearby filament channel by magnetic reconnection. These observations shed light on better understanding on the complexity of filament eruptions.

The Kelvin-Helmholtz (KH) instability is commonly found in many astrophysical, laboratory, and space plasmas. It could mix plasma components of different properties and convert dynamic fluid energy from large scale structure to smaller ones. In this study, we combined the ground-based New Vacuum Solar Telescope (NVST) and the Solar Dynamic Observatories (SDO) / Atmospheric Imaging Assembly (AIA) to observe the plasma dynamics associated with active region 12673 on 09 September 2017. In this multi-temperature view, we identified three adjacent layers of plasma flowing at different speeds, and detected KH instabilities at their interfaces. We could unambiguously track a typical KH vortex and measure its motion. We found that the speed of this vortex suddenly tripled at a certain stage. This acceleration was synchronized with the enhancements in emission measure and average intensity of the 193 \AA{} data. We interpret this as evidence that KH instability triggers plasma heating. The intriguing feature in this event is that the KH instability observed in the NVST channel was nearly complementary to that in the AIA 193 \AA{}. Such a multi-thermal energy exchange process is easily overlooked in previous studies, as the cold plasma component is usually not visible in the extreme ultraviolet channels that are only sensitive to high temperature plasma emissions. Our finding indicates that embedded cold layers could interact with hot plasma as invisible matters. We speculate that this process could occur at a variety of length scales and could contribute to plasma heating.

Magnetic flux ropes (MFRs) are thought to be the central structure of solar eruptions, and their ideal MHD instabilities can trigger the eruption. Here we performed a study of all the MFR configurations that lead to major solar flares, either eruptive or confined, from 2011 to 2017 near the solar disk center. The coronal magnetic field is reconstructed from observed magnetograms, and based on magnetic twist distribution, we identified the MFR, which is defined as a coherent group of magnetic field lines winding an axis with more than one turn. It is found that 90% of the events possess pre-flare MFRs, and their three-dimensional structures are much more complex in details than theoretical MFR models. We further constructed a diagram based on two parameters, the magnetic twist number which controls the kink instability (KI), and the decay index which controls the torus instability (TI). It clearly shows lower limits for TI and KI thresholds, which are $n_{\rm crit} = 1.3$ and $|T_w|_{\rm crit} = 2$, respectively, as all the events above $n_{\rm crit}$ and nearly 90% of the events above $|T_w|_{\rm crit}$ erupted. Furthermore, by such criterion, over 70% of the events can be discriminated between eruptive and confined flares, and KI seems to play a nearly equally important role as TI in discriminating between the two types of flare. There are more than half of events with both parameters below the lower limits, and 29% are eruptive. These events might be triggered by magnetic reconnection rather than MHD instabilities.

Solar flares are often associated with coronal eruptions, but there are confined ones without eruption, even for some X-class flares. How such large flares occurred and why they are confined are still not well understood. Here we studied a confined X2.2 flare in NOAA 12673 on 2017 September 6. It exhibits two episodes of flare brightening with rather complex, atypical ribbons. Based on topology analysis of extrapolated coronal magnetic field, we revealed that there is a two-step magnetic reconnection process during the flare. Prior to the flare, there is a magnetic flux rope (MFR) with one leg rooted in a rotating sunspot. Neighboring to the leg is a magnetic null-point structure. The sunspot drives the MFR to expand, pushing magnetic flux to the null point, and reconnection is first triggered there. The disturbance from the null-point reconnection triggers the second reconnection, i.e., a tether-cutting reconnection below the rope. However, these two reconnections failed to produce an eruption, because the rope is firmly held by its strapping flux. Furthermore, we compared this flare with an eruptive X9.3 flare in the same region with 2 hours later, which has a similar MFR configuration. The key difference between them is that, for the confined flare, the MFR is fully below the threshold of torus instability, while for the eruptive one, the MFR reaches entirely above the threshold. This study provides a good evidence supporting that reconnection alone may not be able to trigger eruption, rather, MHD instability plays a more important role.

We present the Data Release 4&5 quasar catalog from the quasar survey by Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), which includes quasars observed between September 2015 and June 2017. There are a total of 19,253 quasars identified by visual inspections of the spectra. Among them, 11,458 are independently discovered by LAMOST, in which 3296 were reported by SDSS DR12 and DR14 quasar catalog after our survey began, while the rest 8162 are new discoveries of LAMOST. We provide the emission line measurements for the Halpha, Hbeta, MgII and/or CIV for 18100 quasars. Since LAMOST does not have absolute flux calibration information, we obtain the monochromatic continuum luminosities by fitting the SDSS photometric data using the quasar spectra, and then estimate the black hole masses. The catalog and spectra for these quasars are available online. This is the third installment in the series of LAMOST quasar survey which has released spectra for totally ~43,000 quasars hitherto. There are 24,772 independently discovered quasars, 17,128 of which are newly discovered. In addition to the great supplement to the new quasar discoveries, LAMOST has also provided a large database (overlapped with SDSS) for investigating the quasar spectral variability and discovering unusual quasars, including changing-look quasars, with ongoing and upcoming large surveys.

Three-dimensional magnetic topology is crucial to understanding the explosive release of magnetic energy in the corona during solar flares. Much attention has been given to the pre-flare magnetic topology to identify candidate sites of magnetic reconnection, yet it is unclear how the magnetic reconnection and its attendant topological changes shape the eruptive structure and how the topology evolves during the eruption. Here we employed a realistic, data-constrained magnetohydrodynamic simulation to study the evolving magnetic topology for an X9.3 eruptive flare that occurred on 2017 September 6. The simulation successfully reproduces the eruptive features and processes in unprecedented detail. The numerical results reveal that the pre-flare corona contains multiple twisted flux systems with different connections, and during the eruption, these twisted fluxes form a coherent flux rope through tether-cutting-like magnetic reconnection below the rope. Topological analysis shows that the rising flux rope is wrapped by a quasi-separatrix layer, which intersects itself below the rope, forming a topological structure known as hyperbolic flux tube, where a current sheet develops, triggering the reconnection. By mapping footpoints of the newly-reconnected field lines, we are able to reproduce both the spatial location and, for the first time, the temporal separation of the observed flare ribbons, as well as the dynamic boundary of the flux rope's feet. Futhermore, the temporal profile of the total reconnection flux is comparable to the soft X-ray light curve. Such a sophisticated characterization of the evolving magnetic topology provides important insight into the eventual understanding and forecast of solar eruptions.

Seven different models are applied to the same problem of simulating the Sun's coronal magnetic field during the solar eclipse on 2015 March 20. All of the models are non-potential, allowing for free magnetic energy, but the associated electric currents are developed in significantly different ways. This is not a direct comparison of the coronal modelling techniques, in that the different models also use different photospheric boundary conditions, reflecting the range of approaches currently used in the community. Despite the significant differences, the results show broad agreement in the overall magnetic topology. Among those models with significant volume currents in much of the corona, there is general agreement that the ratio of total to potential magnetic energy should be approximately 1.4. However, there are significant differences in the electric current distributions; while static extrapolations are best able to reproduce active regions, they are unable to recover sheared magnetic fields in filament channels using currently available vector magnetogram data. By contrast, time-evolving simulations can recover the filament channel fields at the expense of not matching the observed vector magnetic fields within active regions. We suggest that, at present, the best approach may be a hybrid model using static extrapolations but with additional energization informed by simplified evolution models. This is demonstrated by one of the models.

The rare case of changing-look (CL) AGNs, with the appearance or disappearance of broad Balmer emission lines within a few years, challenges our understanding of the AGN unified model. We present a sample of 21 new CL AGNs at $0.08<z<0.58$, which doubles the number of such objects known to date. These new CL AGNs were discovered by several ways, from (1) repeat spectra in the SDSS, (2) repeat spectra in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and SDSS, and (3) photometric variability and new spectroscopic observations. We use the photometric data from surveys, including the SDSS imaging survey, the Pan-STARRS1, the DESI Legacy imaging survey, the Wide-field Infrared Survey Explorer (WISE), the Catalina Real-time Transient Survey, and the Palomar Transient Factory. The estimated upper limits of transition timescale of the CL AGNs in this sample spans from 0.9 to 13 years in the rest frame. The continuum flux in the optical and mid-infrared becomes brighter when the CL AGNs turn on, or vice versa. Variations of more than 0.2 mag in $W1$ band were detected in 15 CL AGNs during the transition. The optical and mid-infrared variability is not consistent with the scenario of variable obscuration in 10 CL AGNs at more than $3\sigma$ confidence level. We confirm a bluer-when-brighter trend in the optical. However, the mid-infrared WISE colors $W1-W2$ become redder when the objects become brighter in the $W1$ band, possibly due to a stronger hot dust contribution in the $W2$ band when the AGN activity becomes stronger. The physical mechanism of type transition is important for understanding the evolution of AGNs.

We present a clear case study on the occurrence of two successive X-class flares including a decade-class flare (X9.3) and two coronal mass ejections (CMEs) triggered by shearing motion and sunspot rotation in active region NOAA 12673 on 2017 September 6. A shearing motion between the main sunspots with opposite polarities started on September 5 and even lasted after the second X-class flare on September 6. Moreover, the main sunspot with negative polarity rotated around its umbral center and another main sunspot with positive polarity also exhibited a slow rotation. The sunspot with negative polarity at the northwest of active region also began to rotate counter-clockwise before the onset of the first X-class flare. The successive formation and eruption of two S-shaped structures were closely related to the counter-clockwise rotation of three sunspots. It is also found that the rotation of sunspots is faster during four hours prior to the onset of the flares than the period before. The existence of a flux rope is found prior to the onset of two flares by using non-linear force free field extrapolation based on the vector magnetograms observed by SDO/HMI. These results suggest that shearing motion and sunspot rotation play an important role in the buildup of the free energy and the formation of flux ropes in the corona which produces solar flares and CMEs.

Three-dimensional magnetic topology of solar flare plays a crucial role in understanding its explosive release of magnetic energy in the corona. However, such three-dimensional coronal magnetic field is still elusive in direct observation. Here we realistically simulate the magnetic evolution during the eruptive process of a great flare, using a numerical magnetohydrodynamic model constrained by observed solar vector magnetogram. The numerical results reveal that the pre-flare corona contains multi-set twisted magnetic flux, which forms a coherent rope during the eruption. The rising flux rope is wrapped by a quasi-separatrix layer, which intersects itself below the rope, forming a hyperbolic flux tube and magnetic reconnection is triggered there. By tracing the footprint of the newly-reconnected field lines, we reproduce both the spatial location and its temporal evolution of flare ribbons with an expected accuracy in comparison of observed images. This scenario strongly confirms the three-dimensional version of standard flare model.

Further studies on the corner singularity of GS reconstruction are are compiled in this paper. It's focused on solution of the Data Completion (DC) problem with the Extended Hilbert Transform (EHT) over plane rectangular region. Optimal selections of the regular parameters in{\it Tikhonov} solution of corresponding DC problem are developed in this study. The 6-parameter regular solutions and the {\it Jacobian} matrix and one {\it Hessian} tensor to the regular parameters are derived in this work. A concise formula for EHT in the near field of corners, which shows property for EHT near the corner, are also provided. It serves as the additional constraints for our parameter optimization problem (OP). Third, a nonlinear convex function defined by the regular solution and the corner constraints is introduced and is taken as the object function for the OP of the 6 regular parameters on half-space $(\mathbf{p}>0)$. Given an initial guess of $\mathbf{p}_0$, the optimal parameters are solved from the OP through a well known constrained nonlinear optimization method. Last, the benchmark tests to the proposed solution approach are carried out, and detailed results from totally 9 different bench-cases are tabulated. In contrast to solutions with given regular parameters, our bench results demonstrate that an objective way for selection of the optimal $\mathbf{p}$ is successfully laid out here. Robustness and efficiency of the suggested new approach are also highlighted in this study.

As a fundamental magnetic structure in the solar corona, electric current sheets (CSs) can form either prior to or during solar flare, and they are essential for magnetic energy dissipation in the solar corona by enabling magnetic reconnection. However static reconstruction of CS is rare, possibly due to limitation inherent in available coronal field extrapolation codes. Here we present the reconstruction of a large-scale pre-flare CS in solar active region 11967 using an MHD-relaxation model constrained by SDO/HMI vector magnetogram. The CS is found to be associated with a set of peculiar homologous flares that exhibit unique X-shaped ribbons and loops occurring in a quadrupolar magnetic configuration. This is evidenced by that the field lines traced from the CS to the photosphere form an X shape which nearly precisely reproduces the shape of the observed flare ribbons, suggesting that the flare is a product of the dissipation of the CS through reconnection. The CS forms in a hyperbolic flux tube, which is an intersection of two quasi-separatrix layers. The recurrence of the X-shaped flares might be attributed to the repetitive formation and dissipation of the CS, as driven by the photospheric footpoint motions. These results demonstrate the power of data-constrained MHD model in reproducing CS in the corona as well as providing insight into the magnetic mechanism of solar flares.

With SDO observations and a data-constrained MHD model, we identify a confined multi-ribbon flare occurred on 2010 October 25 in solar active region 11117 as a magnetic bald patch (BP) flare with strong evidences. From the photospheric magnetic field observed by SDO/HMI, we find there is indeed magnetic BPs on the PILs which match parts of the flare ribbons. From the 3D coronal magnetic field derived from a MHD relaxation model constrained by the vector magnetograms, we find strikingly good agreement of the BP separatrix surface (BPSS) footpoints with the flare ribbons, and the BPSS itself with the hot flaring loop system. Moreover, the triggering of the BP flare can be attributed to a small flux emergence under the lobe of the BPSS, and the relevant change of the coronal magnetic field through the flare is well reproduced by the pre-flare and post-flare MHD solutions, which match the corresponding pre and post-flare AIA observations, respectively. Our work contributes to the study of non-typical flares that constitute the majority of solar flares but cannot be explained by the standard flare model.

Observations of solar wind turbulence indicate the existence of multi-scale pressure-balanced structures (PBSs) in the solar wind. In this work, we conduct a numerical simulation to investigate multi-scale PBSs and in particular their formation in compressive MHD turbulence. By the use of a higher order Godunov code Athena,a driven compressible turbulence with an imposed uniform guide field is simulated. The simulation results show that both the magnetic pressure and the thermal pressure exhibit a turbulent spectrum with a Kolmogorov-like power law, and that in many regions of the simulation domain they are anti-correlated. The computed wavelet cross-coherence spectrum of the magnetic pressure and the thermal pressure, as well as their space series, indicate the existence of multi-scale PBSs, with the small PBSs being embedded in the large ones. These multi-scale PBSs are likely to be related with the highly oblique-propagating slow-mode waves, as the traced multi-scale PBS is found to be traveling in a certain direction at a speed consistent with that predicted theoretically for a slow-mode wave propagating in the same direction.

In the solar atmosphere, the jets are ubiquitous and found to be at various spatia-temporal scales. They are significant to understand energy and mass transport in the solar atmosphere. Recently, the high-speed transition region jets are reported from the observation. Here we conduct a numerical simulation to investigate the mechanism in their formation. Driven by the supergranular convection motion, the magnetic reconnection between the magnetic loop and the background open flux occurring in the transition region is simulated with a two-dimensional magnetohydrodynamics model. The simulation results show that not only a fast hot jet, much resemble the found transition region jets, but also a adjacent slow cool jet, mostly like classical spicules, is launched. The force analysis shows that the fast hot jet is continually driven by the Lorentz force around the reconnection region, while the slow cool jet is induced by an initial kick through the Lorentz force associated with the emerging magnetic flux. Also, the features of the driven jets change with the amount of the emerging magnetic flux, giving the varieties of both jets. These results will inspire our understanding of the formation of the prevalence of both the fast hot jet and slow cool jet from the solar transition region and chromosphere.

In the solar atmosphere, jets are prevalent and they are significant for the mass and energy transport. Here we conduct numerical simulations to investigate the mass and energy contributions of the recently observed high-speed jets to the solar wind. With a one-dimensional hydrodynamic solar wind model, the time-dependent pulses are imposed at the bottom to simulate the jets. The simulation results show that without other energy source, the injected plasmas are accelerated effectively to be a transonic wind with a substantial mass flux. The rapid acceleration occurs close to the Sun, and the resulting asymptotic speed, number density at 0.3 AU, as well as mass flux normalized to 1 AU are compatible with in situ observations. As a result of the high speed, the imposed pulses generate a train of shocks traveling upward. By tracing the motions of the injected plasma, it is found that these shocks heat and accelerate the injected plasmas successively step by step to push them upward and eventually allow them to escape. The parametric studies show that increasing the speed of the imposed pulses or their temperature gives a considerably faster, and hotter solar wind, while increasing their number density or decreasing their recurring period only bring a denser solar wind. These studies provide a possibility that the ubiquitous high-speed jets are a substantial mass and energy contributions to the solar wind.

The prediction of the arrival time for fast coronal mass ejections (CMEs) and their associated shocks is highly desirable in space weather studies. In this paper, we use two shock propagation models, i.e. Data Guided Shock Time Of Arrival (DGSTOA) and Data Guided Shock Propagation Model (DGSPM), to predict the kinematical evolution of interplanetary shocks associated with fast CMEs. DGSTOA is based on the similarity theory of shock waves in the solar wind reference frame, and DGSPM on the non-similarity theory in the stationary reference frame. The inputs are the kinematics of the CME front at the maximum speed moment obtained from the geometric triangulation method applied to STEREO imaging observations together with the Harmonic Mean approximation. The outputs provide the subsequent propagation of the associated shock. We apply these models to the CMEs on 2012 January 19, January 23, and March 7. We find that the shock models predict reasonably well the shock's propagation after the impulsive acceleration. The shock's arrival time and local propagation speed at Earth predicted by these models are consistent with in situ measurements of WIND. We also employ the Drag-Based Model (DBM) as a comparison, and find that it predicts a steeper deceleration than the shock models after the rapid deceleration phase. The predictions of DBM at 1 AU agree with the following ICME or sheath structure, not the preceding shock. These results demonstrate the applicability of the shock models used here for future arrival time prediction of interplanetary shocks associated with fast CMEs.

We study the physical mechanism of a major X-class solar flare that occurred in the super NOAA active region (AR) 12192 using a data-driven numerical magnetohydrodynamic (MHD) modeling complemented with observations. With the evolving magnetic fields observed at the solar surface as bottom boundary input, we drive an MHD system to evolve self-consistently in correspondence with the realistic coronal evolution. During a two-day time interval, the modeled coronal field has been slowly stressed by the photospheric field evolution,which gradually created a large-scale coronal current sheet, i.e., a narrow layer with intense current, in the core of the AR. The current layer was successively enhanced until it became so thin that a tether-cutting reconnection between the sheared magnetic arcades was set in, which led to a flare. The modeled reconnecting field lines and their footpoints match well the observed hot flaring loops and the flare ribbons, respectively, suggesting that the model has successfully "reproduced" the macroscopic magnetic process of the flare. In particular, with simulation, we explained why this event is a confined eruption-the consequent of the reconnection is the shared arcade instead of a newly formed flux rope. We also found much weaker magnetic implosion effect comparing to many other X-class flares

Solar post-flare loops (PFLs) are arcade-like loop systems that appear during the gradual phases of eruptive flares. The extreme ultraviolet (EUV) observations from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) allow us to investigate the fine structures in PFLs. In this work, we focus on studying the dark post-flare loops (DPFLs) during X-class flares, which are more evident in SDO/AIA data than in previous EUV data. We identify and analyze the DPFLs observed by SDO and find that: (1) the DPFLs of an X5.4 flare have an average lifetime of 10.0 $\pm$ 5.5 minutes, an average width of 1022 $\pm$ 339 km, and an average maximum length of 33 $\pm$ 10 Mm, (2) blob-like falling features with a size close to the resolution of SDO/AIA are identified in the DPFLs and they have an average velocity of 76 $\pm$ 19 km s$^{-1}$, and (3) the average widths of the DPFLs slightly increase with the characteristic temperatures in AIA 304, 171, 193, and 211 {\AA} channels. Our investigation shows that DPFLs are found in all of the 20 cases within this study, which suggests that they are a common phenomenon in X-class flares and are probably produced by the same mechanism that creates coronal rain.

Turbulence is a chaotic flow regime filled by irregular flows. The dissipation of turbulence is a fundamental problem in the realm of physics. Theoretically, dissipation cannot be ultimately achieved without collisions, and so how turbulent kinetic energy is dissipated in the nearly collisionless solar wind is a challenging problem. Wave particle interactions and magnetic reconnection are two possible dissipation mechanisms, but which mechanism dominates is still a controversial topic. Here we analyze the dissipation region scaling around a solar wind magnetic reconnection region. We find that the magnetic reconnection region shows a unique multifractal scaling in the dissipation range, while the ambient solar wind turbulence reveals a monofractal dissipation process for most of the time. These results provide the first observational evidences for the intermittent multifractal dissipation region scaling around a magnetic reconnection site, and they also have significant implications for the fundamental energy dissipation process.

Solar active region (AR) 11283 is a very magnetically complex region and it has produced many eruptions. However, there exists a non-eruptive filament in the plage region just next to an eruptive one in the AR, which gives us an opportunity to perform a comparison analysis of these two filaments. The coronal magnetic field extrapolated using a CESE-MHD-NLFFF code (Jiang & Feng 2013) reveals that two magnetic flux ropes (MFRs) exist in the same extrapolation box supporting these two filaments, respectively. Analysis of the magnetic field shows that the eruptive MFR contains a bald-patch separatrix surface (BPSS) co-spatial very well with a pre-eruptive EUV sigmoid, which is consistent with the BPSS model for coronal sigmoids. The magnetic dips of the non-eruptive MFRs match H{\alpha} observation of the non-eruptive filament strikingly well, which strongly supports the MFR-dip model for filaments. Compared with the non-eruptive MFR/filament (with a length of about 200 Mm), the eruptive MFR/filament is much smaller (with a length of about 20 Mm), but it contains most of the magnetic free energy in the extrapolation box and holds a much higher free energy density than the non-eruptive one. Both the MFRs are weakly twisted and cannot trigger kink instability. The AR eruptive MFR is unstable because its axis reaches above a critical height for torus instability, at which the overlying closed arcades can no longer confine the MFR stably. On the contrary, the quiescent MFR is very firmly held by its overlying field, as its axis apex is far below the torus-instability threshold height. Overall, this comparison investigation supports that MFR can exist prior to eruption and the ideal MHD instability can trigger MFR eruption.

In the solar corona, magnetic flux rope is believed to be a fundamental structure accounts for magnetic free energy storage and solar eruptions. Up to the present, the extrapolation of magnetic field from boundary data is the primary way to obtain fully three-dimensional magnetic information of the corona. As a result, the ability of reliable recovering coronal magnetic flux rope is important for coronal field extrapolation. In this paper, our coronal field extrapolation code (CESE-MHD-NLFFF, Jiang & Feng 2012) is examined with an analytical magnetic flux rope model proposed by Titov & Demoulin (1999), which consists of a bipolar magnetic configuration holding an semi-circular line-tied flux rope in force-free equilibrium. By using only the vector field in the bottom boundary as input, we test our code with the model in a representative range of parameter space and find that the model field is reconstructed with high accuracy. Especially, the magnetic topological interfaces formed between the flux rope and the surrounding arcade, i.e., the "hyperbolic flux tube" and "bald patch separatrix surface", are also reliably reproduced. By this test, we demonstrate that our CESE-MHD-NLFFF code can be applied to recovering magnetic flux rope in the solar corona as long as the vector magnetogram satisfies the force-free constraints.

Measurements of solar wind turbulence reveal the ubiquity of discontinuities. In this study, we investigate how the discontinuities, especially rotational discontinuities (RDs), are formed in magnetohydrodynamic (MHD) turbulence. In a simulation of the decaying compressive three-dimensional (3-D) MHD turbulence with an imposed uniform background magnetic field, we detect RDs with sharp field rotations and little variations of magnetic field intensity as well as mass density. At the same time, in the de Hoffman-Teller (HT) frame, the plasma velocity is nearly in agreement with the Alfv\'{e}n speed, and is field-aligned on both sides of the discontinuity. We take one of the identified RDs to analyze in details its 3-D structure and temporal evolution. By checking the magnetic field and plasma parameters, we find that the identified RD evolves from the steepening of the Alfv\'{e}n wave with moderate amplitude, and that steepening is caused by the nonuniformity of the Alfv\'{e}n speed in the ambient turbulence.

The dynamic process of coronal mass ejections (CMEs) in the heliosphere provides us the key information for evaluating CMEs' geo-effectiveness and improving the accurate prediction of CME induced Shock Arrival Time (SAT) at the Earth. We present a data constrained three dimensional (3D) magnetohydrodynamic (MHD) simulation of the evolution of the CME in a realistic ambient solar wind for the July 12-16, 2012 event by using the 3D COIN-TVD MHD code. A detailed comparison of the kinematic evolution of the CME between the observations and the simulation is carried out, including the usage of the time-elongation maps from the perspectives of both Stereo A and Stereo B. In this case study, we find that our 3D COIN-TVD MHD model, with the magnetized plasma blob as the driver, is able to re-produce relatively well the real 3D nature of the CME in morphology and their evolution from the Sun to Earth. The simulation also provides a relatively satisfactory comparison with the in-situ plasma data from the Wind spacecraft.

Though coronal mass ejections (CMEs) are magnetized fully-ionized gases, a recent observational study of a CME collision event in 2008 November has suggested that their behavior in the heliosphere is like elastic balls, and their collision is probably super-elastic \citep{Shen_etal_2012}. If this is true, this finding has an obvious impact on the space weather forecasting because the direction and veliocity of CMEs may change. To verify it, we numerically study the event through three-dimensional MHD simulations. The nature of CMEs' collision is examined by comparing two cases. In one case the two CMEs collide as observed, but in the other, they do not. Results show that the collision leads to extra kinetic energy gain by 3%--4% of the initial kinetic energy of the two CMEs. It firmly proves that the collision of CMEs could be super-elastic.

Superhalo electrons appear to be continuously present in the interplanetary medium, even at very quiet times, with a power-law spectrum at energies above $\sim$2 keV. Here we numerically investigate the generation of superhalo electrons by magnetic reconnection in the solar wind source region, using the MHD and test particle simulations for both single X-line reconnection and multiple X-line reconnection. We find that the direct current electric field, produced in the magnetic reconnection region, can accelerate electrons from an initial thermal energy of T $\sim10^5$ K up to hundreds of keV. After acceleration, some of the accelerated electrons, together with the nascent solar wind flow driven by the reconnection, propagate upwards along the newly-opened magnetic field lines into the interplanetary space, while the rest move downwards into the lower atmosphere. Similar to the observed superhalo electrons at 1 AU, the flux of the upward-traveling accelerated electrons versus energy displays a power-law distribution at $\sim$ 2 $-$ 100 keV, $f(E) \sim E^{-\delta}$, with a $\delta$ of $\sim$ 1.5 $-$ 2.4. For single (multiple) X-line reconnection, the spectrum becomes harder (softer) as the anomalous resistivity parameter $\alpha$ (uniform resistivity $\eta$) increases. These modeling results suggest that the acceleration in the solar wind source region may contribute to superhalo electrons.

Solar filament are commonly thought to be supported in magnetic dips, in particular, of magnetic flux ropes (FRs). In this Letter, from the observed photospheric vector magnetogram, we implement a nonlinear force-free field (NLFFF) extrapolation of a coronal magnetic FR that supports a large-scale intermediate filament between an active region and a weak polarity region. This result is the first in that current NLFFF extrapolations with presence of FRs are limited to relatively small-scale filaments that are close to sunspots and along main polarity inversion line (PIL) with strong transverse field and magnetic shear, and the existence of a FR is usually predictable. In contrast, the present filament lies along the weak-field region (photospheric field strength $\lesssim 100$ G), where the PIL is very fragmented due to small parasitic polarities on both side of the PIL and the transverse field has a low value of signal-to-noise ratio. Thus it represents a far more difficult challenge to extrapolate a large-scale FR in such case. We demonstrate that our CESE--MHD--NLFFF code is competent for the challenge. The numerically reproduced magnetic dips of the extrapolated FR match observations of the filament and its barbs very well, which supports strongly the FR-dip model for filaments. The filament is stably sustained because the FR is weakly twisted and strongly confined by the overlying closed arcades.

We present a magnetic analysis of the formation and eruption of an active region (AR) sigmoid in AR 11283 from 2011 September 4 to 6. To follow the quasi-static evolution of the coronal magnetic field, we reconstruct a time sequence of static fields using a recently developed nonlinear force-free field model constrained by the SDO/HMI vector magnetograms. A detailed analysis of the fields compared with the SDO/AIA observations suggests the following scenario for the evolution of the region. Initially, a new bipole emerges into the negative polarity of a pre-existing bipolar AR, forming a null point topology between the two flux systems. A weakly twisted flux rope (FR) is then built up slowly in the embedded core region, largely through flux-cancellation photospheric reconnections, forming a bald patch separatrix surface (BPSS) separating the FR from its ambient field. The FR grows gradually until its axis runs into a torus instability (TI) domain near the end of the third day, and the BPSS also develops a fully S-shape. Unlike in the case of standard TI, the FR does not erupt instantly since it is still attached at the photosphere along the bald patch (BP) portion of the polarity inversion line. The combined effects of the TI-driven expansion of the FR and the line-tying at the BP tear the FR into two parts with the upper portion freely expelled and the lower portion remaining behind the post-flare arcades. This process dynamically perturbs the BPSS and results in the transient enhanced brightening of the sigmoid. The accelerated expansion of the upper portion of the FR strongly pushes its envelope flux near the null point and triggers breakout reconnection at the null, as evidenced by a remarkable circular flare ribbon, which further facilitates the eruption. We discuss the important implications of these results for the formation and disruption of sigmoid region with FR.

Coronal mass ejections (CMEs) can be continuously tracked through a large portion of the inner heliosphere by direct imaging in visible and radio wavebands. White-light (WL) signatures of solar wind transients, such as CMEs, result from Thomson scattering of sunlight by free electrons, and therefore depend on both the viewing geometry and the electron density. The Faraday rotation (FR) of radio waves from extragalactic pulsars and quasars, which arises due to the presence of such solar wind features, depends on the line-of-sight magnetic field component $B_\parallel$, and the electron density. To understand coordinated WL and FR observations of CMEs, we perform forward magnetohydrodynamic modelling of an Earth-directed shock and synthesise the signatures that would be remotely sensed at a number of widely distributed vantage points in the inner heliosphere. Removal of the background solar wind contribution reveals the shock-associated enhancements in WL and FR. While the efficiency of Thomson scattering depends on scattering angle, WL radiance $I$ decreases with heliocentric distance $r$ roughly according to the expression $I \propto r^{-3}$. The sheath region downstream of the Earth-directed shock is well viewed from the L4 and L5 Lagrangian points, demonstrating the benefits of these points in terms of space weather forecasting. The spatial position of the main scattering site $\mathbf{r}_{\rm sheath}$ and the mass of plasma at that position $M_{\rm sheath}$ can be inferred from the polarisation of the shock-associated enhancement in WL radiance. From the FR measurements, the local $B_{\parallel {\rm sheath}}$ at $\mathbf{r}_{\rm sheath}$ can then be estimated. Simultaneous observations in polarised WL and FR can not only be used to detect CMEs, but also to diagnose their plasma and magnetic field properties.

Numerical reconstruction/extrapolation of coronal nonlinear force-free magnetic field (NLFFF) usually takes the photospheric vector magnetogram as input at the bottom boundary. Magnetic field observed at the photosphere, however, contains force which is in conflict with the fundamental assumption of the force-free model and measurement noise which is unfavorable for practical computation. Preprocessing of the raw magnetogram has been proposed by Wiegelmann, Inhester, and Sakurai (2006) to remove the force and noise for providing better input for NLFFF modeling. In this paper we develop a new code of magnetogram preprocessing which is consistent with our extrapolation method CESE-MHD-NLFFF (Jiang, Feng, and Xiang, 2012; Jiang and Feng, 2012). Basing on a magnetic-splitting rule that a magnetic field can be split into a potential field part and a non-potential part, we split the magnetogram and deal with the two parts separately. Preprocessing of the magnetogram's potential part is based on a numerical potential field model, and the non-potential part is preprocessed using the similar optimization method of Wiegelmann et al (2006). The code is applied to the SDO/HMI data and results show that the method can remove efficiently the force and noise and improve the quality of extrapolation.

Current magnetohydrodynamic (MHD) simulations of the initiation of solar eruptions are still commonly carried out with idealized magnetic field models, whereas the realistic coronal field prior to eruptions can possibly be reconstructed from the observable photospheric field. Using a nonlinear force-free field extrapolation prior to a sigmoid eruption in AR 11283 as the initial condition in a MHD model, we successfully simulate the realistic initiation process of the eruption event, as is confirmed by a remarkable resemblance to the SDO/AIA observations. Analysis of the pre-eruption field reveals that the envelope flux of the sigmoidal core contains a coronal null and furthermore the flux rope is prone to a torus instability. Observations suggest that reconnection at the null cuts overlying tethers and likely triggers the torus instability of the flux rope, which results in the eruption. This kind of simulation demonstrates the capability of modeling the realistic solar eruptions to provide the initiation process.

Due to the absence of direct measurement, the magnetic field in the solar corona is usually extrapolated from the photosphere in numerical way. At the moment, the nonlinear force-free field (NLFFF) model dominates the physical models for field extrapolation in the low corona. Recently we have developed a new NLFFF model with MHD relaxation to reconstruct the coronal magnetic field. This method is based on CESE--MHD model with the conservation-element/solution-element (CESE) spacetime scheme. In this paper, we report the application of the CESE--MHD--NLFFF code to \SDO/HMI data with magnetograms sampled for two active regions (ARs), NOAA AR 11158 and 11283, both of which were very non-potential, producing X-class flares and eruptions. The raw magnetograms are preprocessed to remove the force and then inputted into the extrapolation code. Qualitative comparison of the results with the \SDO/AIA images shows that our code can reconstruct magnetic field lines resembling the EUV-observed coronal loops. Most important structures of the active regions are reproduced excellently, like the highly-sheared field lines that suspend filaments in AR 11158 and twisted flux rope which corresponds to a sigmoid in AR 11283. Quantitative assess of the results shows that the force-free constraint is fulfilled very well in the strong-field regions but apparently not that well in the weak-field regions because of data noise and numerical errors in the small currents.

In recent studies (Feng et al., 2012; Kong et al., 2012), we proposed that source properties of type II radio bursts can be inferred through a causal relationship between the special shape of the type II dynamic spectrum (e.g., bump or break) and simultaneous extreme ultraviolet (EUV)/white light imaging observations (e.g., CME-shock crossing streamer structures). As a further extension of these studies, in this paper we examine the CME event dated on December 31 2007 associated with a multiple type II radio burst. We identify the presence of two spectral bump features on the observed dynamic spectrum. By combining observational analyses of the radio spectral observations and the EUV-white light imaging data, we conclude that the two spectral bumps are resulted from CME-shock propagating across dense streamers on the southern and northern sides of the CME, respectively. It is inferred that the corresponding two type II emissions originate separately from the two CME-shock flanks where the shock geometries are likely quasi-perpendicular or oblique. Since the emission lanes are bumped as a whole within a relatively short time, it suggests that the type IIs with bumps of the study are emitted from spatially confined sources (with a projected lateral dimension smaller than 0.05-0.1 R$_\odot$ at a fundamental frequency level of 20-30 MHz).

We apply a data-driven MHD model to investigate the three-dimensional (3D) magnetic field of NOAA active region (AR) 11117 around the time of a C-class confined flare occurred on 2010 October 25. The MHD model, based on the spacetime conservation-element and solution-element (CESE) scheme, is designed to focus on the magnetic-field evolution and to consider a simplified solar atomsphere with finite plasma $\beta$. Magnetic vector-field data derived from the observations at the photoshpere is inputted directly to constrain the model. Assuming that the dynamic evolution of the coronal magnetic field can be approximated by successive equilibria, we solve a time sequence of MHD equilibria basing on a set of vector magnetograms for AR 11117 taken by the Helioseismic and Magnetic Imager (HMI) on board the {\it Solar Dynamic Observatory (SDO)} around the time of flare. The model qualitatively reproduces the basic structures of the 3D magnetic field, as supported by the visual similarity between the field lines and the coronal loops observed by the Atmospheric Imaging Assembly (AIA), which shows that the coronal field can indeed be well characterized by the MHD equilibrium in most time. The magnetic configuration changes very limited during the studied time interval of two hours. A topological analysis reveals that the small flare is correlated with a bald patch (BP, where the magnetic field is tangent to the photoshpere), suggesting that the energy release of the flare can be understood by magnetic reconnection associated with the BP separatrices. The total magnetic flux and energy keep increasing slightly in spite of the flare, while the computed magnetic free energy drops during the flare with an amount of $\sim 10^{30}$ erg, which seems to be adequate to provide the energy budget of the minor C-class confined flare.

Here is reported in situ observation of energetic electrons (~100-500 keV) associated with magnetic reconnection in the solar wind by the ACE and Wind spacecraft. The properties of this magnetic cloud driving reconnection and the associated energetic electron acceleration problem are discussed. Further analyses indicate that the electric field acceleration and Fermi type mechanism are two fundamental elements in the electron acceleration processes and the trapping effect of the specific magnetic field configuration maintains the acceleration status that increases the totally gained energy.

The Magnetic cloud boundary layer (BL) is a dynamic region formed by the interaction of the magnetic cloud (MC) and the ambient solar wind. In the present study, we comparatively investigate the proton and electron mean flux variations in the BL, in the interplanetary reconnection exhaust (RE) and across the MC-driven shock by using the Wind 3DP and MFI data from 1995 to 2006. In general, the proton flux has higher increments at lower energy bands compared with the ambient solar wind. Inside the BL, the core electron flux increases quasi-isotropically and the increments decrease monotonously with energy from ~30% (at 18 eV) to ~10% (at 70 eV); the suprathermal electron flux usually increases in either parallel or antiparallel direction; the correlation coefficient of electron flux variations in parallel and antiparallel directions changes sharply from ~0.8 below 70 eV to ~0 above 70 eV. Similar results are also found for RE. However, different phenomena are found across the shock where the electron flux variations first increase and then decrease with a peak increment (>200%) near 100 eV. The correlation coefficient of electron flux variations in parallel and antiparallel directions is always around 0.8. The similar behavior of flux variations in BL and RE suggests that reconnection may commonly occur in BL. Our work also implies that the strong energy dependence and direction selectivity of electron flux variations, which are previously thought to have no enough relevance to magnetic reconnection, could be considered as an important signature of solar wind reconnection in the statistical point of view.

Reliable measurements of the solar magnetic field are still restricted to the photosphere, and our present knowledge of the three-dimensional coronal magnetic field is largely based on extrapolation from photospheric magnetogram using physical models, e.g., the nonlinear force-free field (NLFFF) model as usually adopted. Most of the currently available NLFFF codes have been developed with computational volume like Cartesian box or spherical wedge while a global full-sphere extrapolation is still under developing. A high-performance global extrapolation code is in particular urgently needed considering that Solar Dynamics Observatory (SDO) can provide full-disk magnetogram with resolution up to $4096\times 4096$. In this work, we present a new parallelized code for global NLFFF extrapolation with the photosphere magnetogram as input. The method is based on magnetohydrodynamics relaxation approach, the CESE-MHD numerical scheme and a Yin-Yang spherical grid that is used to overcome the polar problems of the standard spherical grid. The code is validated by two full-sphere force-free solutions from Low & Lou's semi-analytic force-free field model. The code shows high accuracy and fast convergence, and can be ready for future practical application if combined with an adaptive mesh refinement technique.

It has been suggested that type II radio bursts are due to energetic electrons accelerated at coronal shocks. Radio observations, however, have poor or no spatial resolutions to pinpoint the exact acceleration locations of these electrons. In this paper, we discuss a promising approach to infer the electron acceleration location by combining radio and white light observations. The key assumption is to relate specific morphological features (e.g. spectral bumps) of the dynamic spectra of type II radio bursts, to imaging features (e.g. CME going into a streamer) along the CME (and its driven shock) propagation. In this study, we examine the CME-streamer interaction for the solar eruption dated on 2003 November 1. The presence of spectral bump in the relevant type II radio burst is identified, which is interpreted as a natural result of the shock-radio emitting region entering the dense streamer structure. The study is useful for further determinations of the location of type II radio burst and the associated electron acceleration by CME-driven shock.

Magnetic field in the solar corona is usually extrapolated from photospheric vector magnetogram using a nonlinear force-free field (NLFFF) model. NLFFF extrapolation needs a considerable effort to be devoted for its numerical realization. In this paper we present a new implementation of the magnetohydrodynamics (MHD)-relaxation method for NLFFF extrapolation. The magneto-frictional approach which is introduced for speeding the relaxation of the MHD system is novelly realized by the spacetime conservation-element and solution-element (CESE) scheme. A magnetic field splitting method is used to further improve the computational accuracy. The bottom boundary condition is prescribed by changing the transverse field incrementally to match the magnetogram, and all other artificial boundaries of the computational box are simply fixed. We examine the code by two types of NLFFF benchmark tests, the Low & Lou (1990) semi-analytic force-free solutions and a more realistic solar-like case constructed by van Ballegooijen et al. (2007). The results show that our implementation are successful and versatile for extrapolations of either the relatively simple cases or the rather complex cases which need significant rebuilding of the magnetic topology, e.g., a flux rope. We also compute a suite of metrics to quantitatively analyze the results and demonstrate that the performance of our code in extrapolation accuracy basically reaches the same level of the present best-performing code, e.g., that developed by Wiegelmann (2004).

In this paper we conduct a data survey searching for well-defined streamer wave events observed by the Large Angle and Spectrometric Coronagraph (LASCO) on-board the Solar and Heliospheric Observatory (SOHO) throughout Solar Cycle 23. As a result, 8 candidate events are found and presented here. We compare different events and find that in most of them the driving CMEs ejecta are characterized by a high speed and a wide angular span, and the CME-streamer interactions occur generally along the flank of the streamer structure at an altitude no higher than the bottom of the field of view of LASCO C2. In addition, all front-side CMEs have accompanying flares. These common observational features shed light on the excitation conditions of streamer wave events. We also conduct a further analysis on one specific streamer wave event on 5 June 2003. The heliocentric distances of 4 wave troughs/crests at various exposure times are determined; they are then used to deduce the wave properties like period, wavelength, and phase speeds. It is found that both the period and wavelength increase gradually with the wave propagation along the streamer plasma sheet, and the phase speed of the preceding wave is generally faster than that of the trailing ones. The associated coronal seismological study yields the radial profiles of the Alfv\'en speed and magnetic field strength in the region surrounding the streamer plasma sheet. Both quantities show a general declining trend with time. This is interpreted as an observational manifestation of the recovering process of the CME-disturbed corona. It is also found that the Alfv\'enic critical point is at about 10 R$_\odot$ where the flow speed, which equals the Alfv\'en speed, is $\sim$ 200 km s$^{-1}$.

Solar magnetic synoptic charts obtained by NSO/Kitt Peak and SOHO/MDI are analyzed for studying the appearance of bipolar magnetic regions (BMRs) during solar minima. As a result, we find the emergence of long-lived BMRs has three typical features. (1) BMRs' emerging rates of the new cycles increase about 3 times faster than those of the old cycles decrease. (2) Two consecutive solar cycles have an overlapping period of near 10 Carrington rotations. During this very short overlapping time interval, BMRs of two cycles tend to concentrate in the same longitudes. (3) About 53% BMRs distribute with a longitudinal distance of 1/8 solar rotation. Such phenomenon suggests a longitudinal mode of m=8 existing during solar minima.