This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Intrusion of Magnetic Peninsula toward the Neighboring Opposite-polarity Region That Triggers the Largest Solar Flare in Solar Cycle 24

Yumi Bamba Institute for Advanced Research (IAR)/Nagoya University
Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
Institute for Space-Earth Environmental Research (ISEE)/Nagoya University
Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
Satoshi Inoue Institute for Space-Earth Environmental Research (ISEE)/Nagoya University
Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
Shinsuke Imada Institute for Space-Earth Environmental Research (ISEE)/Nagoya University
Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
(Received MMMM DD, YYYY; Accepted MMMM DD, YYYY)
Abstract

The largest X9.3 solar flare in solar cycle 24 and the preceding X2.2 flare occurred on September 6, 2017, in the solar active region NOAA 12673. This study aims to understand the onset mechanism of these flares via analysis of multiple observational datasets from the Hinode and Solar Dynamics Observatory and results from a non-linear force-free field extrapolation. The most noticeable feature is the intrusion of a major negative-polarity region, appearing similar to a peninsula, oriented northwest into a neighboring opposite-polarity region. We also observe proxies of magnetic reconnection caused by related to the intrusion of the negative peninsula: rapid changes of the magnetic field around the intruding negative peninsula; precursor brightening at the tip of the negative peninsula, including a cusp-shaped brightening that shows a transient but significant downflow (100km/s\sim 100km/s) at a leg of the cusp; a dark tube-like structure that appears to be a magnetic flux rope that erupted with the X9.3 flare; and coronal brightening along the dark tube-like structure that appears to represent the electric current generated under the flux rope. Based on these observational features, we propose that (1) the intrusion of the negative peninsula was critical in promoting the push-mode magnetic reconnection that forms and grows a twisted magnetic flux rope that erupted with the X2.2 flare, (2) the continuing intrusion progressing even beyond the X2.2 flare is further promoted to disrupt the equilibrium that leads the reinforcement of the magnetic flux rope that erupted with the X9.3 flare.

Sun: flares — Sun: magnetic field — Sun: sunspots — Sun: activity
journal: ApJfacilities: Hinode (SOT, EIS, and XRT), SDO (HMI and AIA)software: SSW, SPEDAS, VAPOR

1 Introduction

Solar eruptions, such as solar flares (Carrington, 1859) and coronal mass ejections (CMEs; Hansen et al., 1971; Tousey, 1973; Gopalswamy, 2016), explosively release magnetic energy that is stored in the solar corona. These energetic phenomena are mainly driven by the topological changes of magnetic fields, such as magnetic reconnection (Priest, & Forbes, 2002); however, the detailed mechanisms of these phenomena are not yet fully known. One of the key considerations regarding solar eruptions is the magnetic flux rope, a bundle of helically twisted magnetic field lines and is sometimes represented by filaments on the solar surface (Labrosse et al., 2010; Mackay et al., 2010). Several studies have proposed the onset mechanisms of the magnetic flux rope eruptions i.e., solar flares and CMEs from both observational and numerical standpoints. Ideal magnetohydrodynamics (MHD) instability provides a strong basis for explaining the onset of solar eruptions as flares and CMEs. The torus instability (Bateman, 1978; Kliem & Török, 2006; Démoulin & Aulanier, 2010; Kliem et al., 2014) is proposed as one of the ideal cases of MHD instability that causes solar eruptions. The threshold is determined by the decay index, n=(z/Bex)(Bex/z)n=-(z/B_{ex})(\partial{B_{ex}}/\partial{z}), where BexB_{ex} denotes the horizontal component of the external field (Kliem & Török, 2006). Although the value of n=1.5n=1.5 is well known as the critical value for the onset of the torus instability, other studies have proposed that eruptions can occur even when nn is lesser or greater than 1.5 (Zuccarello et al., 2015; Syntelis et al., 2017). Some observational and computational studies (Savcheva et al., 2015; Vemareddy, & Demóulin, 2018; Roudier et al., 2018) employing nonlinear force-free field (NLFFF) modeling (Wiegelmann & Sakurai, 2012; Inoue, 2016; Cheng et al., 2017; Guo et al., 2017) support the proposal regarding the torus instability. The helical kink instability (Gerrard et al., 2001; Fan & Gibson, 2003; Török et al., 2004; Török, & Kliem, 2005) is also considered to be an ideal case of MHD instability that causes solar eruptions. This type of instability can grow when the magnetic twist becomes strong enough, and it can grow even if a twisted flux rope does not attain the critical height of the torus instability. Some studies show that there exists a consistency between the helical kink instability mechanism and confined flares (Hassanin & Kliem, 2016; Liu et al., 2016). Other non-ideal MHD models propose that certain kinds of magnetic reconnection, such as magnetic break-out (Antiochos et al., 1999; Karpen et al., 2012) and tether-cutting reconnection (Moore et al., 2001; Moore & Sterling, 2006; Ishiguro & Kusano, 2017), causes the onset of solar eruptions (Schmieder et al., 2015). These magnetic reconnection models emphasize that the magnetic flux rope is formed by magnetic reconnection, while the ideal MHD models emphasize the evolution of the flux rope to an eruption by the MHD instabilities. However, in the cases of both the ideal and non-ideal MHD models, the factor responsible for triggering the initiation of the flux rope formation and eruption is not yet fully understood.

Solar active region (AR) NOAA 12673 appeared on the solar disk on August 28, 2017. It was a small, quiet, and single sunspot containing dominant positive polarity and low negative polarity and regarded as a β\beta-type sunspot according to the Mount Wilson (Hale) classification (Hale et al., 1919). The region did not appear to produce flares until September 2; however, sudden flux emergence was observed in the southeast and northern regions of the original bipole from September 3, 2017, and various M- and C-class flares occurred thereafter. By September 5, 2017, this rapid flux emergence complicated the βδγ\beta\delta\gamma-type magnetic field structure that contains a penumbra enclosing the umbrae of both positive and negative polarities. The βδγ\beta\delta\gamma-type spot is reportedly the magnetic configuration that exhibits the highest correlation with large flares, according to Sammis et al. (2000); in fact, the AR 12673 produced four X-class flares on September 6, 7, and 10, 2017. Consecutive X2.2 and X9.3 flares (Figure 1) occurred on September 6, 2017, with rapid growth of the AR (Yan et al., 2018; Hou et al., 2018). The X9.3 flare is the largest produced flare in solar cycle 24 and caused an increase in high-energy protons and disturbances in the Earth’s ionosphere and magnetosphere. According to Redmon et al. (2018), these geo-magnetic and ionospheric disturbances disrupted HF communications in the Caribbean region during the hurricane response. It is important to understand how and why such a large and geo-effective flare occurs, from the perspectives of both solar physics and space weather. Therefore, in this study, we aim to understand the onset process of the geo-effective X9.3 flare and preceding X2.2 flare.

The remainder of this paper is organized as follows. The observational data, analysis method, and NLFFF modeling are reviewed in Section 2. The results of the observational data analysis are presented in Section 3. The results are interpreted and a conceivable flare onset scenario of the X2.2 and X9.3 flares is discussed in Section 4. We finally summarize and conclude the study in Section 5.

2 Data Analysis and Modeling

2.1 Data Description

The X2.2 and X9.3 flares that are the focus of this study commenced at 08:57 and 11:53 UT, and peaked at 09:10 and 12:02 UT, respectively, on September 6, 2017. Hinode (Kosugi et al., 2007) had tracked the AR constantly from approximately 11:00 UT on September 5 until it disappeared from sight on September 14, and perfectly observed the X2.2 and X9.3 flares on September 6 using its three instruments: Solar Optical Telescope (SOT, Tsuneta et al., 2008), Extreme-ultraviolet Imaging Spectrometer (EIS, Culhane et al., 2007), and X-ray Telescope (XRT, Golub et al., 2007).

Hinode/SOT determined the full polarization states (Stokes-I, Q, U, and V) at 6301.5 Å and 6302.5 Å (Fe I line) with a sampling of 21.5 mÅ using a Spectro-Polarimeter (SP, Lites et al., 2013). We used an SP map, which was scanned from 06:14 UT on September 6, 2017, to analyze three components of the photospheric magnetic field prior to the onset of the two X-class flares. The SP map was obtained with the fast-mapping mode observation that covered a field-of-view of 164×′′{}^{\prime\prime}\times 164′′. The sampling was 0.3′′ along both the slit and slit-scan.

Hinode/EIS successfully observed both the X2.2 and X9.3 flares. We investigated the local Doppler velocities of plasma related to the precursor brightening using the data from EIS. The raster scan was repeated from 06:20 to 15:35 UT on September 6, 2017, for 10 different wavelengths. In this study, we mainly used Fe XIV lines corresponding to 264.79 Å and 274.20 Å, Fe XVI/Fe XXIII lines corresponding to 263.3 Å, and He II lines corresponding to 256.3 Å. The Fe and He II lines are sensitive to emissions from log(T)6.36.4log(T)\sim 6.3-6.4 plasmas and log(T)4.9log(T)\sim 4.9 plasmas, respectively (Culhane et al., 2007). The 2′′×120′′2^{\prime\prime}\times 120^{\prime\prime} slit scanned 30 positions with sparse 4′′4^{\prime\prime} steps, where 4′′4^{\prime\prime} is the distance between the center of the slits; hence, the scanning field-of-view covered 120′′×120′′120^{\prime\prime}\times 120^{\prime\prime}. The exposure time at each position was 4s, and the spectral resolution was 22 mÅ.

Hinode/XRT continuously observed the AR using a standard AR observation program; this program captures coronal images in multiple wavelengths with high time cadence. We investigated the high-temperature coronal loops that appeared to relate to the magnetic reconnection among sheared magnetic fields using the XRT images that were captures with a thin-Be/Al-poly filter covering a 384′′×384′′384^{\prime\prime}\times 384^{\prime\prime} field-of-view. The spatial resolution was 2′′2^{\prime\prime}, and the data cadence of each thin-Be/Al-poly images was approximately 90s.

The Solar Dynamics Observatory (SDO, Pesnell et al., 2012) observed the X2.2 and X9.3 flares using the Helioseismic and Magnetic Imager (HMI, Schou et al., 2012) and the Atmospheric Imaging Assembly (AIA, Lemen et al., 2012). We used the SDO data to investigate the temporal evolution of the photospheric magnetic fields and the precursor activities in the solar atmosphere. Moreover, we primarily used vector magnetic field data referred to as the Space-weather HMI AR Patch (SHARP111hmi.sharp_720s” series downloaded from Joint Science Operations Center (JSOC) at Stanford University (http://jsoc.stanford.edu/), Bobra et al., 2014) data series. The SHARP data series has been calibrated assuming the Milne-Eddington atmosphere; its 180180^{\circ} ambiguity in the horizontal component was resolved via the minimum energy method (Metcalf, 1994; Leka et al., 2009b, 2012). The SHARP data was continually obtained every 12 min, and we considered the data obtained from 00:00 UT on September 6, to 00:00 UT on September 7, 2017, for use in this study. We also used AIA level 1 data222aia.lev1_euv_12s” and “aia.lev1_uv_24s” series that was obtained in the time window of 00:00 UT on September 6, to 00:00 UT on September 7, 2017. AIA observed the AR with 10 wavelengths, but we mainly analyzed the images captured corresponding to 1600Å and 171Å. We used AIA 211Å and 335 Å data to spatially co-align the HMI/SHARP magnetic field data and Hinode/EIS or XRT data (see details in Section 2.2). The data cadence was 24s for 1600Å and 12s for 171Å, 211Å, and 335Å. In addition, we used HMI level 1.5 line-of-sight magnetic field data333hmi.M_45s” series that was observed via HMI filtergram observation for the Fe I line (6173Å) with a 45s cadence in the same time window as that of AIA data to superpose SHARP vector magnetograms and AIA or Hinode/EIS or XRT images. The spatial resolution was 1.0′′1.0^{\prime\prime} for SHARP and HMI/line-of-sight magnetograms, and 1.5′′1.5^{\prime\prime} for AIA images. Note that we avoided using the SDO data obtained from 06:22:34 UT to 07:53:11 UT on September 6, 2017, because these data were not of good quality owing to the lunar transit.

2.2 Analysis Methods for Hinode Data

The Hinode/SOT-SP data was calibrated via the sp_prep procedure (Lites, & Ichimoto, 2013) using the Solar-SoftWare (SSW) package (Freeland, & Handy, 1998). Subsequently, vector magnetic fields were derived from the calibrated Stokes profiles based on the assumption of the Milne-Eddington atmosphere using the inversion code MERLIN at HAO/CSAC (DOI: 10.5065/D6JH3J8D). The minimum energy code ME0 (Leka et al., 2009a, b, 2012), which is based on the minimum energy method of Metcalf (1994), was adopted to resolve the 180 ambiguity in the tangential component of the magnetic fields. As highlighted by Leka et al. (2017); Bamba, & Kusano (2018), the three components of the photospheric magnetic field are influenced by a projection-effect; that is, the line-of-sight component of the magnetic field can be treated as the radial component only when the observing angle θ=0\theta=0. The AR 12673 was located on September 6, 2017, at S09W42, which is close to the south-east limb. Therefore, to investigate locations and shapes of magnetic polarity inversion lines, we transformed the vector magnetic field data obtained via SOT-SP, which is in image plane on the Charge Coupled Device (CCD), to its corresponding heliocentric coordinates. To transform the magnetic field vectors, we used the transformation matrix presented by equation (1) of Gary, & Hagyard (1990). In this manner, we obtained the radial component (BzB_{z}) and the horizontal components (Bx,ByB_{x},B_{y}) of the photospheric magnetic field vectors.

Hinode/EIS raster scan data was calibrated via the eis_prep procedure using the SSW package. For reasons concerning thermal factors, there exists an orbital variation of the line position causing an artificial Doppler shift of ±20km/s\pm 20km/s, which follows a sinusoidal behavior. This orbital variation and wavelength calibration of the line position were corrected using the house keeping data (Kamio et al., 2010). Subsequently, we first performed wavelength calibration using the Fe XIV (264.787Å) line, following which we finalized the central value λobs\lambda^{\prime}_{obs} for the observed Fe XIV line profiles via the single Gaussian fitting. Herein, we integrated the observed line profile over the region that was relatively quiet as per the EIS scanning field-of-view. We then applied the single Gaussian fitting to the integrated line profile and repeated the process for the data that were scanned from 06:21:03444 UT to 06:58:03444It means scan start time. UT on September 6, 2017. We averaged the value of fitted line center over the abovementioned time period and obtained the observed line center λobs_FeXIV\lambda_{obs\_FeXIV}. Subsequently, we determined dw by subtracting λobs_FeXIV\lambda_{obs\_FeXIV} from λlit_FeXIV\lambda_{lit\_FeXIV}, which is the literature wavelength of the Fe XIV line as per literature (Brown et al., 2008): dw=λlit_FeXIVλobsdw=\lambda_{lit\_FeXIV}-\lambda_{obs}. Using dw, we calibrated the other spectral lines as per λ0=λlit𝑑𝑤\lambda_{0}=\lambda_{lit}-{\it dw}, where λlit\lambda_{lit} is the wavelength of the other lines based on literature (Brown et al., 2008), and λ0\lambda_{0} is the reference wavelength of those lines as per the EIS raster scan observation. Following this, we identified the center of observed spectral lines λobs\lambda_{obs} via the eis_auto_fit procedure using the SSW package, which applies single Gaussian fitting to the line profiles. We then generated intensity, Doppler velocity, and line width maps for each of the observed lines. The Doppler velocity was calculated as per v=Δλc/λ0v={\Delta\lambda}c/{\lambda_{0}}, where Δλ=λ0λobs\Delta\lambda=\lambda_{0}-\lambda_{obs}, and cc is the speed of light. Using the intensity, Doppler velocity, and line width maps, we checked the locations and timings of precursor brightening occurrences relative to the two X-class flares. We further applied double Gaussian fitting to the averaged line profile of the Fe XIV lines (264.7Å and 274.2Å) at the single-slit position of the precursor brightening and measured the Doppler velocity of a local and transient flow, which is described in Section 3.2. Note that the blend from the other lines is negligibly small for Fe XIV (264.7Å) and Fe XVI (262.9Å) (Young et al., 2007).

The Hinode/XRT data was calibrated via the xrt_prep procedure (Kobelski et al., 2014) using the SSW package. Subsequently, we first spatially co-aligned the XRT Al-poly images and SDO/AIA 335Å, which show structures similar to each other (Yoshimura & McKenzie, 2015). The AIA 335Å images have been spatially well aligned with the HMI/line-of-sight images via the aia_prep procedure, as we described in Section 2.3. Thus, we aligned the radial magnetic field (BzB_{z}) images of Hinode/SOT-SP or SDO/HMI with the line-of-sight magnetic field images. In this manner, we finally overplotted the XRT intensities onto the BzB_{z} field images or vice versa (Figure 6).

2.3 Analysis Methods for SDO Data

HMI/line-of-sight magnetograms and AIA data were calibrated via the aia_prep procedure using SSW, whereby can be rotated such that the solar EW and NS axes are aligned with the horizontal and vertical axes of the image. The pixel scales, which are different between HMI and AIA, were resampled to the same size. Thus, the positions of the line-of-sight magnetograms and AIA images positions were aligned. Next, we transformed SHARP vector magnetograms to the heliocentric coordinates in the same manner for Hinode/SOT-SP. Therefore, we obtained the radial component (BzB_{z}) and horizontal components (Bx,ByB_{x},B_{y}) of the photospheric magnetic field vectors. We then extracted an line-of-sight magnetogram, SHARP vector magnetogram, and AIA image, which were the closest in terms of time, and superposed the SHARP vector magnetic field and AIA image via the line-of-sight magnetogram. Moreover, we superposed the magnetic polarity inversion lines of the radial component (BzB_{z}) of the SHARP data onto the Hinode/EIS maps corresponding to the intensity, Doppler velocity, and line width maps. For this superposition, the Hinode/EIS maps corresponding to the Fe XIV (264.787Å) line and AIA 211Å images were superposed. The AIA 211Å images and the SHARP magnetic field data had already been co-aligned via the line-of-sight magnetic field data; thus the magnetic polarity inversion lines of BzB_{z} from the SHARP data were overplotted onto the Hinode/EIS maps. We showed the magnetic polarity inversion lines on the EIS maps, which were observed during the precursor brightening (Figure 4). Then magnetic polarity inversion lines of BzB_{z} fields from the SHARP were also overplotted on the Hinode/XRT images (Figure 6) in the same manner, although herein we employed AIA 335Å instead of 211Å, as we mentioned in Section 2.2.

2.4 Nonlinear Force-free Field Modeling

We extrapolated three-dimensional (3D) coronal magnetic fields based on the photospheric magnetic fields before the two X-class flares. We used the modeling of the 3D coronal magnetic field using an NLFFF approximation, as described in Wiegelmann & Sakurai (2012), and employed the MHD relaxation method developed in Inoue et al. (2014); Inoue (2016). The SHARP vector magnetic fields (Bx,By,BzB_{x},B_{y},B_{z}) observed at 08:36 and 11:00 UT on September 6, 2017, were set as the bottom boundary condition. We solved the MHD equations in 3D space starting with a potential field, which was extrapolated from the radial component BzB_{z} of the photospheric magnetic field using the Green function method (Sakurai, 1982). Note that pressure and gravity were neglected here because the low β\beta approximation was well satisfied in the solar corona (Gary, 2001). Therefore, we solved the MHD equations simplified in the zero β\beta approximation. The basic equations are as follows:

ρ=|𝐁|,\rho=|{\bf B}|, (1)
𝐯t=(𝐯)𝐯+1ρ𝐉×𝐁+ν2𝐯,\frac{\partial{\bf v}}{\partial t}=-({\bf v}\cdot{\bf\nabla}){\bf v}+\frac{1}{\rho}{\bf J}\times{\bf B}+\nu{\bf\nabla}^{2}{\bf v}, (2)
𝐁t=×(𝐯×𝐁η𝐉)ϕ,\frac{\partial{\bf B}}{\partial t}={\bf\nabla}\times({\bf v}\times{\bf B}-\eta{\bf J})-{\bf\nabla}\phi, (3)
𝐉=×𝐁,{\bf J}={\bf\nabla}\times{\bf B}, (4)
ϕt+ch2𝐁=ch2cp2ϕ,\frac{\partial\phi}{\partial t}+c^{2}_{h}{\bf\nabla}\cdot{\bf B}=-\frac{c^{2}_{h}}{c^{2}_{p}}\phi, (5)

where ρ\rho, 𝐁{\bf B}, 𝐯{\bf v}, 𝐉{\bf J}, and ϕ\phi are plasma density, magnetic flux density, velocity, electric current density, and a convenient potential to remove errors derived from 𝐁{\bf\nabla}\cdot{\bf B} (Dedner et al., 2002), respectively. In these equations, the length, magnetic field, density, velocity, time, and electric current density are normalized by L=244.8MmL^{*}=244.8Mm, B=2,500GB^{*}=2,500G, ρ\rho^{*} = |B||B^{*}|, VAB/(μ0ρ)1/2V_{\rm A}^{*}\equiv B^{*}/(\mu_{0}\rho^{*})^{1/2}, where μ0\mu_{0} is the magnetic permeability, τAL/VA\tau_{\rm A}^{*}\equiv L^{*}/V_{\rm A}^{*}, and J=B/μ0LJ^{*}=B^{*}/\mu_{0}L^{*}, respectively.

The pseudo density is assumed to be proportional to |𝐁||{\bf B}| to ease the relaxation by equalizing the Alfve´\acute{e}n speed in space. ν\nu is the viscosity fixed at 1.0×1031.0\times 10^{-3}, and the coefficients ch2c_{h}^{2}, cp2c_{p}^{2} in Equation (5) are also fixed with the constant values, 0.04 and 0.1, respectively. The initial density condition is given by ρ=|𝐁|\rho=|{\bf B}|, and the velocity is set to zero in the entire space. The resistivity is given by η=η0+η1|𝐉×𝐁||𝐯|/|𝐁|\eta=\eta_{0}+\eta_{1}|{\bf J}\times{\bf B}||{\bf v}|/|{\bf B}| where η0=5.0×105\eta_{0}=5.0\times 10^{-5} and η1=1.0×103\eta_{1}=1.0\times 10^{-3} in non-dimensional units. The second term is introduced to accelerate the relaxation to the force free field, particularly in weak field regions. At the boundaries, all of velocities are fixed to zero, and the magnetic fields are fixed to the potential fields except at the bottom boundary. More details are presented in Inoue et al. (2018).

3 Results

3.1 Evolution Steps of the X2.2 and X9.3 Flares

The X2.2 and X9.3 flares occurred in an interval of approximately 3h as shown in GOES555Geostationary Operational Environmental Satellite: We employed GOES-13 because GOES-15, which was the latest GOES series as of September 2017, had brief outage of soft X-ray due to lunar transit near the onset time of the X2.2 flare. soft X-ray light curves of Figure 1. The onset time was 08:57 UT for the X2.2 flare and 11:53 UT for the X9.3 flare, on September 6, 2017. Another C2.7 flare occurred before the two X-class flares, and its onset time was 07:29 UT, on September 6, 2017. Figure 2(a-c) shows the temporal evolution of the radial magnetic field BzB_{z} from SDO/HMI SHARP. Panels (a), (b), and (c) are snapshots at 04:58:42 UT (approximately 4h before the X2.2 flare), 08:46:42 UT (approximately 10min before the X2.2 flare) and 11:22:42 UT (approximately 30min before the X9.3 flare). We labeled the major polarities as N1 and N2 for the negative region, and P1 and P2 for the positive region, as shown in panel (a). The most remarkable magnetic field change was the intrusion of a part of the N1, which is indicated by the red circle, toward the northwest direction as indicated by the red arrow in panel (a). The region was shaped like a peninsula and intruded into the neighboring positive polarity region P1. Henceforth, this region is called ”intruding negative peninsula.” It was originally shaped as shown in panel (a), and it seemed to start the intrusion between 06:50 and 07:50 UT, when the HMI observation was not available, or the data were not of good quality because of the lunar transit. Then, the intrusion became rapid from around 08:00 UT, i.e., from approximately 1h before the X2.2 flare onset. The intrusion continued throughout the two X-flares (panel (b)). Then, the intruding negative peninsula finally merged with the N2, dividing the P1, as shown in panel (c).

Precursor brightening666We defined precursor brightening as a small brightening intermittently observed at the same site on a magnetic polarity inversion lines where initial flare ribbons appeared on both sides. A physical interpretation of precursor brightening in terms of flare trigger process is discussed in Bamba, & Kusano (2018). and initial flare ribbons of the two X-flares are observed near the intruding negative peninsula. Panels (d-i) are snapshots of the last precursor brightening (panels (d) and (g)), initial flare ribbons (panels (e) and (h)), and enhanced flare ribbons (panels (f) and (i)) for the X2.2 and X9.3 flares, observed by SDO/AIA 1600Å. Several noticeable features were obtained by comparing the X2.2 and X9.3 flares. First, the last precursor brightening was observed at the tip of the intruding negative peninsula before the X2.2 flare, while the location of the last precursor for the X9.3 flare was slightly south from the tip of the intruding negative peninsula, as indicated by the red arrows in panels (d) and (g). In addition, the shape of the last precursor brightening was dot-like for the X2.2 flare, while it was elongated along the magnetic polarity inversion line between P1 and N1 for the X9.3 flare. Second, the initial flare ribbons of the X2.2 flare appeared near the site where the last precursor brightening was observed. The ribbons over the P1 were labeled as R1 and R2, and those over the N1 and N2 were labeled as R3, R4, and R5, as indicated by the white arrows in panel (e). The ribbons R1-R3 were locating around the the tip of the intruding negative peninsula, although the R4 and R5 were at a slight distance from it. The initial flare ribbons of the X9.3 flare also appeared on both sides of the location of the its precursor brightening, although the ribbons were saturated immediately, as seen in panel (h). Third, the enhanced flare ribbons were extended further south for the X9.3 flare, compared with those of the X2.2 flare. Both the X-flares showed north-south (along the magnetic polarity inversion line between P1 and N1) and east-west (along the magnetic polarity inversion line between P1-N2) ribbons, as shown in panels (f) and (i). In the X2.2 flare, the east-west ribbons R4 and R5 appeared almost simultaneously with the north-south ribbons R1-R3; R5 was fainter than the others. Even after the ribbons were enhanced, their shape and location were not very different in panels (e) and (f). In contrast, the X9.3 flare ribbons were extended more toward west and south directions in panel (i) than those in panel (h). These temporal variations in precursor brightening and flare ribbons can be clearly seen in the supplemental animation of the SDO/AIA 1600 Å.

3.2 Local Structures in the X2.2 and X9.3 Flaring Site

The intruding negative peninsula is likely to related to the onset of the two X-class flares, as described in Section 3.1. Therefore, we further investigated the magnetic field structure and corresponding plasma dynamics near the intruding negative peninsula. Figure 3 (a) shows the three components of magnetic fields in the AR observed by Hinode/SOT-SP, approximately 3h before the X2.2 flare. Hinode/SOT-SP successfully scanned the central part of the AR including the intruding negative peninsula, although it missed the southern half of P2. The magnetic polarity inversion line between P1 and N1, including the intruding negative peninsula, had a strong magnetic gradient in the radial component Bz{B_{z}} and strong magnetic shear in the horizontal components (Bx,By)(B_{x},B_{y}), as shown in panel (a). The magnetic shear distribution is more clearly shown in panel (b), where a strong magnetic shear (<70<-70^{\circ} of shear angles) is distributed along the magnetic polarity inversion line between P1 and N1. In particular, the tip of the intruding negative peninsula was close to 90-90^{\circ} magnetic shear as the vectors of the horizontal fields were almost parallel to the magnetic polarity inversion line in panel (a).

Interestingly, a cusp-shaped brightening was observed at the tip of the intruding negative peninsula approximately 1.5h before the X2.2 flare onset. It was observed as a transient enhancement of soft X-ray as indicated by the red arrow in Figure 1. Figure 4 shows quicklook maps of intensity, Doppler velocity, and line width from Hinode/EIS spectroscopic observations of Fe XIV (264.7Å) line, which is sensitive to emissions from plasma at a temperature of log(T)6.3log(T)\sim 6.3 with overlying magnetic polarity inversion lines of the radial magnetic field BzB_{z} from SDO/HMI SHARP. Hinode/EIS field-of-view successfully covered both the north-south and east-west magnetic polarity inversion lines, and the cusp-shaped brightening can be seen in the intensity map of panel (b). The cusp-shaped brightening bridged over the north-south magnetic polarity inversion lines between P1 and N1 at the tip of the intruding negative peninsula. Hinode/EIS scanned the AR from west to east (from right to left in the maps) in the field-of-view for 3min, and the time in the Figure 4 indicates the scan start time. Therefore, the cusp-shaped brightening was observed approximately at 07:19 UT. The brightening was simultaneously observed by SDO/AIA in 171Å, which is sensitive to emissions from plasma at a temperature of log(T)5.8log(T)\sim 5.8 (Lemen et al., 2012), with considerably higher 12min temporal resolution, as shown in panels (e-i)777Note that data of the AIA 171 Å on panels (e-i) are not good quality because of lunar transit. We could still clearly identify the cusp-shaped brightening and its evolution as shown in panels (e-i).. The brightening started in a small region indicated by the white arrow in panel (e), and propagated to both sides; then, the cusp-shaped brightening appeared (see panels (e-i)). This suggests that there was an energy input at the top of the cusp, and plasma fell down along the magnetic field lines.

Moreover, transient but significant downflow was observed to the west of the intruding negative peninsula, as shown in panel (c). Because Hinode/EIS performed observations with a 2′′2^{\prime\prime} slit, the red-shifted Doppler signal was only seen in one slit position. It can be identified that the red-shifted signal corresponds to the west (right in the image, roots in P1) foot point of the cusp-shaped brightening. The Doppler velocity was more than 80km/s80km/s in the Fe XIV line. Significant line broadening was also observed for the cusp-shape in panel (d). One of the reasons of line broadening is the presence of bi-directional plasma jets (e.g., Innes et al., 1997) or turbulent motion (e.g., Imada et al., 2008) caused by magnetic reconnection. There has been considerable discussion on line broadening associated with magnetic reconnection. The shape of the observed spectral line in our study is similar to the spectral line associated with magnetic reconnection discussed in a previous study (see Figure1 in Innes et al. (1997)). Furthermore, coronal dimming and strong blue-shifted signals were observed in the northeast of the field-of-view in Figure 4(b-d). At 06:21 UT, this coronal dimming had already started and was further enhanced after the X2.2 and X9.3 flares. It is inferred that the coronal dimming related to the expanding coronal loops and corresponding upflows were observed a few days or a day before the flare onset, as suggested by Imada et al. (2014).

We analyzed the spectral line profiles at the site of the cusp-shaped brightening via double Gaussian fitting. Figure 5 shows the spectral line profiles and spectral images of Fe XIV (264.7 Å and 274.2 Å), Fe XVI/Fe XXIII (263.3 Å), and He II (256.3 Å) lines at the single site of the red-shifted signal, shown in Figure 4(c). The vertical axis of the spectral images (a-2, b-2, c-2, and d-2) corresponding to the Y-axis of Figure 4(c), and the white horizontal line-cut denotes the Y-coordinate of the strongest red-shifted signal. The intensity profiles along the white line-cut are plotted in the left for each spectral image (a-1, b-1, c-1, and d-1). The horizontal axis for both intensity profiles and spectral images are converted to denote velocity unit by single (in panels (c) and (d)) or double (in panels (a) and (b)) Gaussian fitting. Evidently, both Fe XIV lines have significant red-shifted components, as shown in panels (a-1) and (b-1). The Doppler velocity of the second components (red-colored profiles) relative to the first components (blue-colored profiles) was approximately 103±8km/s103\pm 8km/s for 264.7 Å and 115±16km/s115\pm 16km/s for 274.2 Å, respectively. Note that the errors are defined as the summation of errors due to photon noise for each component of the double Double Gaussian fitting. Fe XVI/Fe XXIII (263.3 Å) and He II (256.3 Å) lines in panels (c-1) and (d-1) also show red-shifted profiles although it was not well fitted via double Gaussian fitting. These profiles Fe XVI/Fe XXIII and He II lines were also seen to have red-shifted components with a velocity of approximately 100km/s100km/s.

3.3 Global Structures in AR 12673

The global coronal structure of the AR 12673 was observed using Hinode/XRT and SDO/AIA. Figure 6 shows the temporal variation of high-temperature coronal loops observed using Hinode/XRT, in which the magnetic polarity inversion lines of the radial magnetic fields BzB_{z} from SDO/HMI SHARP or Hinode/SOT-SP are overplotted. The red global coronal loops connect the major bipoles P1-N1, P1-N2, and P2-N1 in panel (a). In particular, P1-N1 loops were highly sheared and almost parallel to the magnetic polarity inversion lines. The small loop that is indicated by the white arrow in panel (a) existed in the north of the intruding negative peninsula at least from 00:00 UT on September 6, 2017. A cusp-shaped brightening then appeared just before the C2.7 flare as seen in panel (b). This cusp-shaped brightening (c.f. Tsuneta et al., 1992) corresponded to the one that was observed by Hinode/EIS and SDO/AIA 171 Å in Figures 4 and 5, and the transient enhancement of soft X-ray as denoted by the red arrow in Figure 1. The cusp-shaped brightening remained until the X2.2 flare onset, following which, a dark tube-like feature was observed just before the X9.3 flare onset, as denoted by the sky-ble arrows in panels (c) and (d). This tube-like dark feature was elongated north-south and was along the magnetic polarity inversion lines between P1 and N1.

Figure 7 shows the coronal loops of lower temperature than those shown in Figure 6, observed by SDO/AIA 171 Å. Panels (a) and (b) are snapshots of HMI radial magnetic field and AIA 171 Å image, respectively, taken approximately 45min before the X2.2 flare onset. There was a small loop that is bridging over the magnetic polarity inversion line between P1 and N1 at the tip of the intruding negative peninsula, as indicated by the red arrow in panel (b). This small loop corresponded to the cusp-shaped brightening seen in Figures 4(e-i) and 6(b). Moreover, another bright loop L1 is denoted by the white arrow in panel (b). The loop L1 remained until the X2.2 flare onset, as seen in panel (c). The initial flare ribbons, shown in Figure 2(e), are also seen in panel (c), although the loop L1 only can be seen in the SDO/AIA 171 Å image. The loop L1 seemed to correspond the coronal loop that was rooted in P1 and N1. The loop L2, which has the same connectivity as L1, appeared after the X2.2 flare, as denoted by the white arrow in panel (d). The intensity of L2 gradually decayed, meanwhile a long north-south loop that is connecting P2-N2, denoted by a white arrow in panel (e), gradually appeared. The long north-south loop expanded as is outlined by white dashed line in panel (f) and erupted simultaneously with intense emission of the X9.3 flare onset.

Interestingly, the dark tube-like structure that is shown in Figure 6 (c) and (d) is also observed in the AIA 171 Å image, as denoted by the sky-ble arrow in Figure 7 (e). This tube-like structure is along the magnetic polarity inversion lines between P1 and N1 and is also observed under the bright loop L2 in panel (d). It was not observed in AIA 171 Å image before the X2.2 flare; however, it started to be seen as gradually rise from around 10:45 UT, when the brightening of the X2.2 flare ribbons was decaying as shown in the supplemental animation. The dark tube-like structure then became invisible because of intense emission from the X9.3 flare ribbons in panel (f). It seemed to be erupted by the X9.3 flare because there was no clear structure such as dark tube after the X9.3 flare emission started decaying.

4 Discussion

4.1 Conceivable Scenario of sequential X2.2 and X9.3 Flares

The X-class flares occurred while the intruding negative peninsula was moving toward northwest as we described in Section 3.1. Our results suggest that the intrusion possibly causes variations in magnetic topologies owing to magnetic reconnection along the magnetic polarity inversion line between P1 and N1. Figure 8(a) shows the magnetic field lines that produce the X2.2 flare as blue and magenta tubes on the radial magnetic field image. The positive/negative foot point of magenta (blue) tubes is labeled as F1/F2 (F3/F4), as shown in panel (b), in which the back ground is the image of the initial flare ribbons in AIA 1600 Å. The negative foot points of magenta (F2) and blue (F4) tubes are rooted in intruding the negative peninsula (i.e., N1) and N2 respectively, while the positive foot points of both tubes (F1 and F3) are rooted in P1. Note that a part of the positive foot point of the blue tubes is rooted slightly closer to the intruding negative peninsula although still in P1; thus, we labeled the foot point as F’3.

The negative foot point F2 of the magenta tubes is pressed against the blue tubes as the intruding negative peninsula moves toward northwest. This motion possibly causes “push-mode” magnetic reconnection (Yamada, 1999) that forms a magnetic flux rope erupting with the X2.2 flare. Figures 9(a) and (b) illustrate the before and after versions of push-mode magnetic reconnection that varies magnetic topologies from F1-F2 /F3-F4 (the magenta/blue tubes) to F2-F3/F1-F4 (yellow tubes). It is inferred that F1-F4 develops as a flux rope and the bright loop L1 in Figures 7 (b) and (c) corresponds to overlying magnetic fields on the magnetic flux rope. The overlying loop L1 becomes bright as the flux rope grows and rises toward its eruption, i.e., the X2.2 flare, owing to the intrusion of the intruding negative peninsula. The initial flare ribbons then appear at each foot point as shown in Figure 2 (e). The foot points F1 and F3 in Figure 8 (b) correspond to the initial flare ribbons R1 and R2 on the positive polarity P1. The other ribbons R3, R4, and R5 on the negative polarities N1 and N2 appear at the foot points of F2 and F4. Here, the foot point F4 is rooted by both the flare ribbons R4 and R5.

Moreover, F’3-F4, which is a part of the pre-existing blue tubes and is denoted by the blue dashed line in Figure 9(a), is also reconnected with the magenta tubes F1-F2. F’3-F4 reconnects with F1-F2 earlier than F3-F4 because it is originally located closer to F1-F2 than F3-F4. This magnetic reconnection with F’3-F4 and F1-F2 forms a small loop F2-F’3, which is denoted by the red dashed line in Figure 9(b), and corresponds to the small cusp-shaped brightening that is observed in Hinode/EIS, XRT, and SDO/AIA 171 Å (cf. Figures 4, 6, and 7) before the X2.2 flare.

Figure 8(c) displays the magnetic field lines before the X9.3 flare. The intrusion of the negative peninsula continues even after the X2.2 flare, and it further causes magnetic reconnection under the yellow tubes. A part of yellow tubes then becomes unstable and is lifted. Thus, the initial flare ribbons of the X9.3 flare appear under the yellow tubes as shown in Figure 8 (d). The yellow tubes further reconnect with the green tubes that are rooted in P2. This magnetic reconnection causes variations in the magnetic topologies from F1-F4/F5-F6 (the yellow/green loops) to F4-F6/F1-F5 (thick and thin red loops), as illustrated in Figures 9(c) and (d). Furthermore, the initial flare ribbons propagate to the P2 region, where F6 is rooted in, and are elongated toward the south more than that of the X2.2 flare, as shown in Figure 2(i).

The magnetic flux rope F1-F4, which is indicated by the thick yellow line in Figure 9(c), might correspond to the dark tube-like structure that was observed in the Hinode/XRT and SDO/AIA 171 Å images between the X2.2 and X9.3 flares. The tube-like structure started to observed in SDO/AIA 171 Å from 10:45 UT (approximately 1.5h before the X9.3 flare onset), while it was not seen in SDO/AIA 171 Å and Hinode/XRT before and just after the X2.2 flare. Furthermore, we observed precursor brightening that is elongated along the magnetic polarity inversion line between P1 and N1 before the X9.3 flare onset, as denoted by the red arrow in Figures 2(g) and 7(e). Those observational results suggest that successive magnetic reconnection under the flux rope F1-F4 occur and reinforce the flux rope (cf. Fan, 2010).

Therefore, according to the results obtained from observational data analysis, the intrusion of the negative peninsula that presses the magenta magnetic field lines to the blue magnetic field lines promote push-mode magnetic reconnection that causes the X2.2 flare. Moreover, the successive intrusion of the negative peninsula further plays an important role to cause successive magnetic reconnection and to hasten the eruption of the yellow magnetic flux rope by an MHD instability. Zou et al. (2020) also stated the importance of the first X2.2 flare reconnection to facilitate the rise of the flux rope and its eruption, i.e., the X9.3 flare, by torus instability. The magnetic polarity inversion line between P1 and N1 has strong magnetic gradient, resulting in strong electric current. Therefore, there is a possibility that magnetic reconnection between sheared magnetic fields over the P1-N1 magnetic polarity inversion lines occurs and forms a magnetic flux rope by consuming a considerable amount of time, even without the intrusion of the negative peninsula. Nevertheless, our results suggest that the intrusion of the negative peninsula encourages the formation and destabilization of magnetic flux ropes both in the X2.2 and X9.3 flares.

4.2 Interpretation of the Local and Transient Downflow

In this study, we found local and transient downflow related to the cusp-shaped brightening as we described in Section 3.2. Generally, Lorentz force driven flows such as reconnection out flow do not exhibit clear temperature dependence while pressure gradient driven flows such as chromospheric evaporation exhibit the strong temperature dependence (e.g., Fisher et al., 1985; Watanabe et al., 2010; Imada et al., 2015). We observed the downflow of 100km/s\sim 100km/s at both high temperature Fe lines (log(T)6.36.4log(T)\sim 6.3-6.4) and low temperature He II line (log(T)4.9log(T)\sim 4.9). The temperature is believed to rises proportionately to height in the solar atmosphere in the region where is higher than the chromosphere (e.g., Vernazza et al., 1981). Therefore, each ion’s sensitive temperature is suggested to correspond to the formation height of each emission line. Further, the downflow illustrate in Figure 4(c) is driven by Lorentz force and that there is an energy input, such as magnetic reconnection, that occurred in a higher atmosphere than formation height of the Hinode/EIS Fe lines shown in Figure 5. The velocity of the downflow was approximately 100km/s100km/s, which is significantly lower than the Alfve´\acute{e}n velocity (order of 1,000km/s1,000km/s) in the corona, in all the Fe and He lines shown in Figure 5. The small cusp-shaped high-temperature loop that was observed using Hinode/XRT in Figure 6(b) and magnetic topology in Figure 8 also supports this interpretation.

As we discussed in Section 4.1, the push-mode magnetic reconnection between the magenta and blue sheared magnetic fields due to the intrusion of the negative peninsula can form a small magnetic loop connecting P1 and N1 at the tip of the intruding negative peninsula, as well as the cusp-shaped brightening. Moreover, downflow was observed in the west leg of the cusp rather than the top of the cusp, as seen in Figure 4(b, c). Therefore, it is suggested that the downflow indicates a secondary flow in which coronal plasma decelerate and fall down along the coronal loop (e.g., Imada et al., 2013; Polito et al., 2018), rather than the reconnection out-flow itself (e.g., Yokoyama, & Shibata, 1998).

Bamba et al. (2017) found similar local and transient upflow using lower temperature lines (log(T)4.04.8log(T)\sim 4.0-4.8), which were observed by the IRIS888Interface Region Imaging Spectrograph satellite (De Pontieu et al., 2014), than the Hinode/EIS observed lines in this study. In their case, temperature independent strong upflow with a velocity of up to 100km/s100km/s was observed, corresponding to a brightening in the chromosphere before an X-class flare. They interpreted that the local and transient upflow suggest the presence of an energy input from magnetic reconnection in a lower temperature region lower than the middle chromosphere. Unfortunately, data from IRIS, which can observe emission lines from a lower temperature atmosphere such as that used in Bamba et al. (2017), was not available for this AR 12673. Even so, it is suggested that our downflow is driven by Lorentz force because both the high-temperature Fe lines and low-temperature He II line showed red-shifted signals with the velocities of approximately 100km/s100km/s.

Moreover, in Bamba et al. (2017), the local and transient upflow, which was observed approximately 30min before an X-class flare onset by IRIS, is reported as a proxy of precursory phenomenon that directly represents the trigger of an X-class flare. Their upflow was very transient but the corresponding brightening in the chromosphere continued until 10min before an X-class flare onset. In our case, it seems that the location of the cusp-shaped brightening and downflow corresponds to the location of the push-mode magnetic reconnection, which plays a crucial role in triggering the X2.2 flare. However, it is unreasonable to assume that those features represent the direct trigger of the X2.2 flare because those were observed more than 1.5h before the flare onset. Therefore, it is suggested that the cusp-shaped brightening and downflow are proxies of the push-mode magnetic reconnection that gradually occurs at the tip of the intruding negative peninsula due to its intrusion from before the X2.2 flare. As we described in Section 4.1, the magnetic connectivity of the cusp-shaped brightening seems to be formed by push-mode magnetic reconnection between the magenta and blue magnetic arcades in Figure 8 (a) and (b). This topological change possibly forms a magnetic flux rope corresponding to the F1-F4 loop and is erupted with the X2.2 flare.

5 Summary

In the present study, we examined the onset process of the consecutive X2.2 and X9.3 flares because the X9.3 flare is the largest flare in the solar cycle 24 and impacted the space weather of the Earth. We analyzed observational data from the Hinode and SDO satellites and interpreted the results using supplemental NLFFF extrapolation. As a result, we found that the key feature to trigger the consecutive X-class flares is the intrusion of the negative peninsula towards the neighboring opposite polarity region in the northwest.

The intrusion is crucial to promote the “push-mode” magnetic reconnection (Yamada, 1999) between the sheared magnetic field along the magnetic polarity inversion line between the major bipoles P1 and N1 in the AR, and to form a twisted magnetic flux rope. The “push-mode” magnetic reconnection is introduced based on reconnection experiments performed in laboratories, as well as another “pull-mode” magnetic reconnection. Both the push-mode and pull-mode reconnection can occur on a vertical electric current sheet formed between large-scale sheared magnetic fields by converging or diverging flows, respectively. Flux cancellation (van Ballegooijen & Martens, 1989; Green et al., 2011) on the photosphere corresponds to the push-mode reconnection.

Therefore, we propose an onset scenario of the consecutive X-class flares as follows: the flux rope that is formed and evolved by the intrusion of the negative peninsula erupts and causes the X2.2 flare; the successive intrusion of the negative peninsula promotes further breaking of the magnetic equilibrium, leading to reinforcement of the magnetic flux rope; the flux rope then erupts due to destabilization of the magnetic field and causes the X9.3 flare.

We presented observational evidences that support the above scenario. First, magnetic field data showed the intrusion of a part of a major negative sunspot (i.e., negative peninsula) towards the opposite polarity region in the northwest, and precursor brightening was correspondingly observed at the tip of the intruding negative peninsula. Second, the cusp-shaped brightening and corresponding downflow that can be proxies of push-mode magnetic reconnection in the presented magnetic topology during the intrusion of the negative peninsula were detected by both imaging and spectroscopic observations. Third, the magnetic flux rope was observed in extreme-ultraviolet and X-rays between the X2.2 and X9.3 flares, and it erupted as the X9.3 flare. However, it was difficult to clearly identify whether the flux rope formed by X2.2 flare was necessary for the X9.3 flare occurrence. We require further studies to examine reveal this problem.

We further demonstrated the usefulness of complementary observational data analysis by combining observations, such as Hinode and SDO data. Those satellites are unique owing to their high spatial resolution, high temporal cadence, and multiple observational methods such as imaging, spectroscopy, and spectropolarimetry. It would be useful to investigate the detailed physical mechanisms of flare triggering through multiple observations and MHD simulations. Therefore, further understanding of the physical mechanism of flare triggers with a combination of observational and numerical studies is expected.

The authors appreciate all the people working on Hinode for the perfect operation to obtain such high-quality data of AR 12673. Moreover, the authors thank to the proposers of the observations that were postponed owing to the AR 12673 flare observation for their understanding. The HMI and AIA data have been used courtesy of NASA/SDO and the AIA and HMI science teams. The authors are grateful for the science teams. This work was partly carried out at the Hinode Science Center at NAOJ and Nagoya University, Japan, and the Solar Data Analysis System (SDAS) operated by the Astronomy Data Center in cooperation with the Hinode Science Center of the NAOJ. A part of this study was carried by using the computational resource of the Center for Integrated Data Science, Institute for Space-Earth Environmental Research, Nagoya University through the joint research program. The authors and part of this work were supported by MEXT/JSPS KAKENHI Grant Numbers JP15H05812, JP15H05814, JP15H05816, and JP15K21709. YB deeply appreciate the many instructive comments from Dr. D. H. Brooks and Dr. K. D. Leka, discussion with Dr. K. Kusano and Dr. J. Muhamad, and the kind technical support from Dr. T. Hara. The authors thank the anonymous referee for their valuable comments that improved the clarity of the manuscript. We would like to thank Editage (www.editage.com) for English language editing. We evaluate the NLFFF model extrapolated from the photospheric magnetic field, focusing on data taken at 02:30 UT on September 6, 2017. The temporal evolution of residual Lorentz-force R=|J×B|2𝑑VR=\int|\vec{J}\times\vec{B}|^{2}dV and divergence of magnetic field D=|B|2𝑑VD=\int|\vec{\nabla}\cdot\vec{B}|^{2}dV during the NLFFF calculation are plotted in Figures 10(a) and (b), where dVdV is a volume element. Both profiles decrease as the evolution progresses. We change the boundary magnetic field during the calculation; in particular, the horizontal component is according to a linear combination, Bt=γBt(obs)+(1γ)Bt(pot)\vec{B}_{t}=\gamma\vec{B}_{t(obs)}+(1-\gamma)\vec{B}_{t(pot)}, where Bt(obs)\vec{B}_{t(obs)} and Bt(pot)\vec{B}_{t(pot)} correspond to the horizontal component of the observed and potential magnetic fields measured at the bottom surface, respectively. γ\gamma increases as per the following formulation: γ=γ+δγ\gamma=\gamma+\delta\gamma falls below RminR_{min}, where RminR_{min} is set at 1.0×1021.0\times 10^{-2} as a given parameter. If γ\gamma satisfies 1.01.0, Bt\vec{B}_{t} is consistent with observed magnetic field Bt(obs)\vec{B_{t(obs)}}. As shown in Figure 10(a), the value of RR saturates around 1.0×1021.0\times 10^{-2} at a very early stage of the calculation because the horizontal magnetic field at the bottom boundary is changing from a potential field to the observed magnetic field. When γ\gamma is 11, the magnetic field is relaxed towards to the NLFFF. Since the value of RR slightly increases after 8,0008,000 iterations, we took the magnetic field at 12,00012,000 iterations, where the value of RR falls below 1.0×1061.0\times 10^{-6}. Figures 10(c) and (d) show the extreme ultra violet (EUV) image taken by AIA at 131Å and the NLFFFF superimposed on the image, respectively. The colors of the field lines correspond to the value of |J|/|B||\vec{J}|/|\vec{B}|, and the strong brightening region shown in the EUV image is well captured by the NLFFF.

References

  • Antiochos et al. (1999) Antiochos, S. K., DeVore, C. R., & Klimchuk, J. A. 1999, ApJ, 510, 485
  • Bamba et al. (2017) Bamba, Y., Inoue, S., Kusano, K., & Shiota, D. 2017, ApJ, 838, 134
  • Bamba et al. (2017) Bamba, Y., Lee, K.-S., Imada, S., & Kusano, K. 2017, ApJ, 840, 116
  • Bamba, & Kusano (2018) Bamba, Y., & Kusano, K. 2018, ApJ, 856, 43.
  • Bateman (1978) Bateman, G. 1978, Cambridge, Mass., MIT Press, 1978. 270 p.,
  • Bobra et al. (2014) Bobra, M. G., Sun, X., Hoeksema, J. T., et al. 2014, Sol. Phys., 289, 3549.
  • Brown et al. (2008) Brown, C. M., Feldman, U., Seely, J. F., et al. 2008, The Astrophysical Journal Supplement Series, 176, 511.
  • Culhane et al. (2007) Culhane, J. L., Harra, L. K., James, A. M., et al. 2007, Sol. Phys., 243, 19
  • Carrington (1859) Carrington, R. C. 1859, MNRAS, 20, 13
  • Culhane et al. (2007) Culhane, J. L., Harra, L. K., James, A. M., et al. 2007, Sol. Phys., 243, 19
  • Cheng et al. (2017) Cheng, X., Guo, Y., & Ding, M., 2017, ScChE, 60, 1383
  • Dedner et al. (2002) Dedner, A., Kemm, F., Kröner, D., et al. 2002, Journal of Computational Physics, 175, 645
  • Démoulin & Aulanier (2010) Démoulin, P., & Aulanier, G. 2010, ApJ, 718, 1388
  • De Pontieu et al. (2014) De Pontieu, B., Title, A. M., Lemen, J. R., et al. 2014, Solar Phys., 289, 2733
  • Fan & Gibson (2003) Fan, Y., & Gibson, S. E. 2003, ApJL, 589, L10
  • Fan (2010) Fan, Y. 2010, ApJ, 719, 728
  • Fisher et al. (1985) Fisher, G. H., Canfield, R. C., & McClymont, A. N. 1985, ApJ, 289, 425
  • Freeland, & Handy (1998) Freeland, S. L., & Handy, B. N. 1998, Sol. Phys., 182, 497
  • Gerrard et al. (2001) Gerrard, C. L., Arber, T. D., Hood, A. W., & Van der Linden, R. A. M. 2001, A&A, 373, 1089
  • Gary, & Hagyard (1990) Gary, G. A., & Hagyard, M. J. 1990, Sol. Phys., 126, 21.
  • Gary (2001) gary, G. A. 2001, Sol. Phys., 203, 71
  • Golub et al. (2007) Golub, L., Deluca, E., Austin, G., et al. 2007, Sol. Phys., 243, 63
  • Gopalswamy (2016) Gopalswamy, N. 2016, Geoscience Letters, 3, 8
  • Green et al. (2011) Green, L. M., Kliem, B., & Wallace, A. J. 2011, A&A, 526, A2
  • Guo et al. (2017) Guo, Y., Cheng, X., & Ding, M. 2017, Science in China Earth Sciences, 60, 1408
  • Hale et al. (1919) Hale, G. E., Ellerman, F., Nicholson, S. B., & Joy, A. H. 1919, ApJ, 49, 153
  • Hansen et al. (1971) Hansen, R. T., Garcia, C. J., Grognard, R. J.-M., et al. 1971, Proceedings of the Astronomical Society of Australia, 2, 57
  • Hassanin & Kliem (2016) Hassanin, A., & Kliem, B. 2016, ApJ, 832, 106
  • Hou et al. (2018) Hou, Y. J., Zhang, J., Li, T., et al. 2018, A&A, 619, A100
  • Imada et al. (2008) Imada, S., Hara, H., Watanabe, T., et al. 2008, ApJ, 679, L155
  • Imada et al. (2013) Imada, S., Aoki, K., Hara, H., et al. 2013, ApJ, 776, L11
  • Imada et al. (2014) Imada, S., Bamba, Y., & Kusano, K. 2014, Publications of the Astronomical Society of Japan, 66, S17
  • Imada et al. (2015) Imada, S., Murakami, I.,& Watanabe, T. 2015, Phys. Plasmas, 22, 101206
  • Innes et al. (1997) Innes, D. E., Inhester, B., Axford, W. I., et al. 1997, Nature, 386, 811
  • Inoue et al. (2014) Inoue, S., Magara, T., Pandey, V. S., et al. 2014, ApJ, 780, 101
  • Inoue (2016) Inoue, S. 2016, Progress in Earth and Planetary Science, 3, 19
  • Inoue et al. (2018) Inoue, S., Shiota, D., Bamba, Y., et al. 2018, ApJ, 867, 83.
  • Ishiguro & Kusano (2017) Ishiguro, N., & Kusano, K. 2017, ApJ, 843, 101
  • Kamio et al. (2010) Kamio, S., Hara, H., Watanabe, T., Fredvik, T., Hansteen, V. H., Solar Phys., 266, 209 (2010)
  • Karpen et al. (2012) Karpen, J. T., Antiochos, S. K., & DeVore, C. R. 2012, ApJ, 760, 81
  • Kliem & Török (2006) Kliem, B., Török, T. 2006, Physical Review Letters, 96, 255002
  • Kliem et al. (2014) Kliem, B., Lin, J., Forbes, T. G., Priest, E. R., & Török, T. 2014, ApJ, 789, 46
  • Kobelski et al. (2014) Kobelski, A. R., Saar, S. H., Weber, M. A., et al. 2014, Sol. Phys., 289, 2781
  • Kosugi et al. (2007) Kosugi, T., Matsuzaki, K., Sakao, T., et al. 2007, Solar Phys., 243, 3
  • Leka et al. (2009a) Leka, K. D., Barnes, G., & Crouch, A. 2009, The Second Hinode Science Meeting: Beyond Discovery-Toward Understanding, 415, 365
  • Labrosse et al. (2010) Labrosse, N., Heinzel, P., Vial, J.-C., et al. 2010, Space Sci. Rev., 151, 243
  • Leka et al. (2009b) Leka, K. D., Barnes, G., Crouch, A. D., et al. 2009, Sol. Phys., 260, 83
  • Leka et al. (2012) Leka, K. D., Barnes, G., Gary, G. A., Crouch, A. D., & Liu, Y. 2012, Sol. Phys., 276, 441
  • Leka et al. (2017) Leka, K. D., Barnes, G., & Wagner, E. L. 2017, Sol. Phys., 292, 36.
  • Lemen et al. (2012) Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Solar Phys., 275, 17
  • Lites et al. (2013) Lites, B. W., Akin, D. L., Card, G., et al. 2013, Sol. Phys., 283, 579.
  • Lites, & Ichimoto (2013) Lites, B. W., & Ichimoto, K. 2013, Sol. Phys., 283, 601.
  • Liu et al. (2016) Liu, R., Kliem, B., Titov, V. S., et al. 2016, ApJ, 818, 148
  • Mackay et al. (2010) Mackay, D. H., Karpen, J. T., Ballester, J. L., et al. 2010, Space Sci. Rev., 151, 333
  • Metcalf (1994) Metcalf, T. R. 1994, Sol. Phys., 155, 235
  • Moore et al. (2001) Moore, R. L., Sterling, A. C., Hudson, H. S., & Lemen, J. R. 2001, ApJ, 552, 833
  • Moore & Sterling (2006) Moore, R. L., & Sterling, A. C. 2006, Washington DC American Geophysical Union Geophysical Monograph Series, 165, 43
  • Pesnell et al. (2012) Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, Solar Phys., 275, 3
  • Polito et al. (2018) Polito, V., Testa, P., Allred, J., et al. 2018, ApJ, 856, 178
  • Priest, & Forbes (2002) Priest, E. R., & Forbes, T. G. 2002, A&A Rev., 10, 313
  • Redmon et al. (2018) Redmon, R. J., Seaton, D. B., Steenburgh, R., et al. 2018, Space Weather, 16, 1190
  • Redmon et al. (2018) Redmon, R. J., Seaton, D. B., Steenburgh, R., et al. 2018, Space Weather, 16, 1190
  • Roudier et al. (2018) Roudier, T., Schmieder, B., Filippov, B., et al. 2018, A&A, 618, A43
  • Tousey (1973) Tousey, R. 1973, Space Research Conference, 713
  • Young et al. (2007) Young, P. R., Del Zanna, G., Mason, H. E., et al. 2007, PASJ, 59, S857
  • Sakurai (1982) Sakurai, T. 1982, Sol. Phys., 76, 301.
  • Sammis et al. (2000) Sammis, I., Tang, F., & Zirin, H. 2000, ApJ, 540, 583.
  • Savcheva et al. (2015) Savcheva, A., Pariat, E., McKillop, S., et al. 2015, ApJ, 810, 96
  • Schmieder et al. (2015) Schmieder, B., Aulanier, G., & Vršnak, B. 2015, Sol. Phys., 290, 3457
  • Schou et al. (2012) Schou, J., Scherrer, P. H., Bush, R. I., et al. 2012, Solar Phys., 275, 229
  • Syntelis et al. (2017) Syntelis, P., Archontis, V., & Tsinganos, K. 2017, ApJ, 850, 95
  • Török et al. (2004) Török, T., Kliem, B., & Titov, V. S. 2004, A&A, 413, L27
  • Török, & Kliem (2005) Török, T., & Kliem, B. 2005, ApJ, 630, L97
  • Tsuneta et al. (1992) Tsuneta, S., Hara, H., Shimizu, T., et al. 1992, PASJ, 44, L63
  • Tsuneta et al. (2008) Tsuneta, S., Ichimoto, K., Katsukawa, Y., et al. 2008, Sol. Phys., 249, 167
  • van Ballegooijen & Martens (1989) van Ballegooijen, A. A., & Martens, P. C. H. 1989, ApJ, 343, 971
  • Vemareddy, & Demóulin (2018) Vemareddy, P., & Demóulin, P. 2018, ApJ, 857, 90
  • Vernazza et al. (1981) Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635
  • Watanabe et al. (2010) Watanabe, T., Hara, H., Sterling, A. C., & Harra, L. K. 2010, ApJ, 719, 213
  • Wiegelmann & Sakurai (2012) Wiegelmann, T., & Sakurai, T. 2012, Living Reviews in Solar Physics, 9, 5
  • Yamada (1999) Yamada, M. 1999, J. Geophys. Res., 104, 14529
  • Yan et al. (2018) Yan, X. L., Wang, J. C., Pan, G. M., et al. 2018, ApJ, 856, 79
  • Yokoyama, & Shibata (1998) Yokoyama, T., & Shibata, K. 1998, ApJ, 494, L113
  • Yoshimura & McKenzie (2015) Yoshimura, K., & McKenzie, D. E. 2015, Sol. Phys., 290, 2355
  • Young et al. (2007) Young, P. R., Del Zanna, G., Mason, H. E., et al. 2007, PASJ, 59, S857
  • Zou et al. (2020) Zou, P., Jiang, C., Wei, F., et al. 2020, ApJ, 890, 10
  • Zuccarello et al. (2015) Zuccarello, F. P., Aulanier, G., & Gilchrist, S. A. 2015, ApJ, 814, 126
Refer to caption
Figure 1: GOES soft X-ray light curve from 00:00 UT on September 6 to 00:00 UT on September 7, 2017. The solid/broken lines represent 1-8 Å/0.5-4 Å light curves. The onset times of the X2.2 and X9.3 flares, which are indicated by the black arrows, are 08:57 and 11:53 UT on September 6, 2017, respectively. A small C2.7 flare occurs at 07:29 UT before the X-class flares, as indicated by the third black arrow. The red arrow indicates soft X-ray intensity enhancement corresponding to a transient brightening and downflow observed by Hinode/EIS. The brief data gap near 05:00 UT is due to the lunar transit.
Refer to caption
Figure 2: (a-c) Temporal evolution of the radial magnetic field BzB_{z} observed by SDO/HMI SHARP. The background white/black areas indicate the positive/negative polarities of BzB_{z} in the range of ±1,000G{\pm 1,000}G, and the magnetic polarity inversion lines are indicated by the green lines. The major polarities are labeled as N1 and N2 for the negative regions and P1 and P2 for the positive regions. The intruding negative peninsula is denoted by the red circle in panel (a). (d-i) Temporal variation of brightening in SDO/AIA 1600Å in the range of 0-1,000 DN. The left, middle, and right columns depict snapshots of the precursor brightening, initial flare ribbons, and enhanced flare ribbons, respectively. Panels (d-f) and (g-i) correspond to the X2.2 and X9.3 flares, respectively. The magnetic polarity inversion lines of BzB_{z} are overplotted with the green lines. The red arrows in panels (d, g) indicate the location of the last precursor brightening that is observed just before the X2.2 and X9.3 flare onsets, and the white arrows in panel (e) indicate the X2.2 initial flare ribbons. The field-of-views of panels (d-i) are slightly larger than those of panels (a-c). The supplemental animation of the AIA 1600 Å (without the magnetic polarity inversion lines) is available.
Refer to caption
Figure 3: (a) Distribution of three components of the magnetic field (BxB_{x}, ByB_{y}, BzB_{z}) in AR 12673, which was observed by Hinode/SOT-SP fast-map scanning from 06:14:05 UT on September 6, 2017. The background white/black and green lines are similar to those in Figure 2(a-c), except that the background BzB_{z} level is in the range of ±1,500G{\pm 1,500}G. The small arrows indicate the orientation and strength of the horizontal components (BxB_{x}, ByB_{y}) of the photospheric magnetic fields, whose magnitudes are greater than 100G. (b) Distribution of the magnetic shear angles relative to the potential field. The red/blue regions indicates counter-clockwise/clockwise shear in the range of ±180{\pm 180^{\circ}}. The black lines represent the magnetic polarity inversion lines, similar to the green lines in panel (a). The red/black circles in panel (a)/(b) denote the intruding negative peninsula. Additionally, the southern half of the P2 is located outside of the SP scanning field-of-view.
Refer to caption
Figure 4: Intensity (upper-right, panel (b)), Doppler velocity (lower-left, panel (c)), and line width (lower-right, panel (d)) in Fe XIV (264.7Å) line observed by Hinode/EIS. The start time of the scanning is 07:17:59 UT on September 6, 2017, and it takes approximately 3min to scan from the west to east (right to left) in the field-of-view. The color scales for intensity, Doppler velocity, and line width are depicted on the right. Panel (a) is a reference Bz{B_{z}} magnetogram observed by SDO/HMI SHARP at 06:00 UT on September 6, 2017, which is closest in time to the Hinode/EIS observation of panels (b-d). The intensity, Doppler velocity, and line width maps and Bz{B_{z}} magnetogram are co-aligned with each other, and the green colored magnetic polarity inversion lines in panel (a) are overplotted on panels (b-d) denoted in red, black, and white lines, respectively. Temporal evolution of the small cusp-shaped brightening in panel (b) is depicted in panels (e-i), which are coronal images captured by SDO/AIA 171Å in the field-of-view of the red rectangle in panel (a). The cusp-shaped brightening in panel (b) starts as a tiny brightening that is indicated by the white arrow in panel (e).
Refer to caption
Figure 5: Spectral line profiles and spectral images at the single site of the red-shifted signal in Figure 4(c). Each panel is corresponds to the Fe XIV (264.7 Å), Fe XIV (274.2 Å), Fe XVI/Fe XXIII (263.3 Å), and He II (256.3 Å) lines, respectively. The intensity profiles along the horizontal white line-cut in the spectral images on the right are plotted in the images on the left. The vertical axes of the profiles represent intensity, and the vertical axes of the spectral images are represented in spatial arcseconds along the slit. The horizontal axes represent the velocities resulting from single (in panels (c) and (d)) and double (in panels (a) and (b)) Gaussian fitting. The intensity range of the spectral images is depicted in the color-scale on the right of each panel. The histogram with the black solid line depicts the intensity of each spectral line. The black broken lines are the result of Gaussian fitting, and the blue/red curves denote the first/second components obtained from double Gaussian fitting for the intensity profiles in panels (a-1) and (b-1). The Doppler velocities of the second components relative to the first components are approximately 103±8km/s103\pm 8km/s for 264.7 Å and 115±16km/s115\pm 16km/s for 274.2 Å, respectively.
Refer to caption
Figure 6: Temporal variation of high-temperature coronal structures observed by Hinode/XRT. The red-colored X-ray images are displayed over the radial magnetic field (BzB_{z}) images that are closest in time to the X-ray images in panels (a) and (b). The radial magnetic field image is obtained from the SDO/HMI SHARP series in panel (a), whereas the radial field image in panels (b) is obtained by Hinode/SOT-SP, because of the data gap in the SDO observation. The white arrow in panel (a) indicates a small bright loop in the north of the intruding negative peninsula. The background Hinode/XRT images in panels (c) and (d) are the same, whereas the magnetic polarity inversion lines of the radial magnetic field (BzB_{z} obtained from SDO/HMI SHARP) are overlaid only in panel (c), depicted as green lines. The sky-ble arrows in panels (c) and (d) indicate a dark tube-like structure along the magnetic polarity inversion line between P1 and N1.
Refer to caption
Figure 7: Temporal variation of the SDO/AIA 171 Å images, whose characteristic emission temperature is log(T)5.8log(T)\sim 5.8. Panel (a) depicts a reference magnetic field image at 08:11 UT on September 6, 2017, the same time to panel (b). The magnetic polarity inversion lines of the radial magnetic field (BzB_{z}) obtained from SDO/HMI SHARP that are closest in time to the background SDO/AIA 171 Å images are overlaid in panels (b-d). The intruding negative peninsula is indicated by the red circle in panel (a). The red arrow in panel (b) indicates the cusp-shaped bright loop that are seen in Figures 4(b) and 6(b). The white arrows in panels (b) and (c) are indicate the same bright loop L1 at different times before the X2.2 flare onset. Another white arrows in panel (c) indicates the initial flare ribbons of the X2.2 flare that is seen in Figure 2(e). The white arrow in panel (d) points to the bright loop L2 that appears between the X2.2 and X9.3 flares. In panel (e), a dark tube-like structure along the magnetic polarity inversion line between P1 and N1, which is seen in Figure 6(c, d), is indicated by the sky-blue arrow. The white arrows in panels (e) and (f) indicate the same long north-south loop before and after the X9.3 flare onset, and the white dashed line in panel (f) outlines the expanded north-south loop. The supplemental animation of the SDO/AIA 171 Å (without the magnetic polarity inversion lines) is available.
Refer to caption
Figure 8: Three-dimensional coronal magnetic field structures before the onset of the X2.2 and X 9.3 flares, which are reconstructed by the NLFFF extrapolation. All the images depict AR as seen from the southeast (i.e., from the bottom-left in the images of Figure 2). The upper two panels (a) and (b) are correspond to the magnetic field structure before onset of the X2.2 flare (calculated from SDO/HMI SHARP on 08:36 UT on September 6, 2017), whereas the bottom two panels (c) and (d) are correspond to that before the onset of the X9.3 flare (calculated from SDO/HMI SHARP on 11:00 UT on September 6, 2017). The different colored tubes indicate magnetic field lines with different magnetic connectivities. Panels (a) and (c) depict the side view of the magnetic fields with a bottom image of the radial magnetic field (BzB_{z} obtained from SDO/HMI SHARP). Panels (b) and (d) depict the same view with the background of AIA 1600 Å images at the time of the initial flare ribbons at 08:59:27 and 11:55:27 UT on September 6, 2017, which are the same images as Figures 2 (e) and (h), respectively.
Refer to caption
Figure 9: Schematics of magnetic field lines before and after the onset of the X2.2 and X9.3 flares. All the images depict AR as seen from the southeast (i.e., from the bottom-left in the images of Figure 2), similar to that in Figure 8. The colored lines in panels (a) and (c) represent the same colored tubes in Figure 8(a) and (c), and are correspond to the snapshots before the onset of the X2.2 and X9.3 flares, respectively. The red and yellow solid lines in panels (b) and (d) indicate the magnetic field lines after the onset of the X2.2 and X9.3 flares. The small red dashed line loop corresponds to the cusp-shaped brightening depicted in Figures 4 and 6.
Refer to caption
Figure 10: (a)-(b) Temporal evolution of residual Lorentz-force and divergence of magnetic field, which are defined as R=|J×B|2dVR=\int|\vec{J}\times\vec{B|}^{2}dV and D=|B|2D=\int|\vec{\nabla}\cdot\vec{B}|^{2}dV, respectively, where dVdV is a volume element. We considered the magnetic field as the NLFFF at 12,00012,000 iterations denoted by the red circle. (c) EUV image obtained from AIA 131 Å observed at 08:36 UT on September 6, 2017. (d) The NLFFF extrapolated from the photospheric magnetic field observed at 08:30 UT on September 6, 2017, is superimposed on the image depicted in (c).