PASJ: Publ. Astron. Soc. Japan 51, 553-563 (1999) A Quadruple Magnetic Source Model for Arcade Flares and X-Ray Arcade Formations outside Active Regions I. Dark Filament Suspension and the Magnetic Structure in the Pre-Event Regions Yutaka UCHIDA, 1 Shigenobu Masaya T O R I I , 1 Shuhei 1 H I R O S E , 1 Samuel CABLE, 2 Satoshi M O R I T A , 1 U E M U R A , 1 and Tomotaka YAMAGUCHI 1 Department of Physics, Science University of Tokyo, Shinjuku-ku, Tokyo 162-8601 2 Physics Department, Auburn University, Alabama, USA E-mail (YU): [email protected] (Received 1999 April 9; accepted 1999 June 3) Abstract The high-sensitivity, wide dynamic-range observations by the Soft X-ray Telescope (SXT) aboard Yohkoh has enabled us to look into the faint pre-event structures of arcade flares, and of even fainter X-ray arcade formation events outside active regions. What we have found in the pre-event structure of the latter, however, was not a sagged simple (bipolar) arcade, as expected in the classical model, but a "dual-arcades" type structure in which the inside legs of each "arcade" cross with the other's, landing at the closer part of the domain of the other. Similar features, together with some other features inexplicable in the classical arcade flare model, were also found in strong arcade flares in active regions seen axis-on at the limb. These features raised a severe problem with the classical "reclosing of the once opened simple arcade" model(s). In the present paper, we propose interpretations of what we discovered by Yohkoh-SXT by reviving a quadruple source model proposed by one of the authors (YU) years ago, pointing out a serious difficulty in the classical model(s). This model, based on the quadruple magnetic sources in the photosphere, has a "neutral sheet" already in the pre-event phase in the corona above the field polarity-reversal line in the photosphere, and turns out to explain quite nicely the structures of both faint pre-event corona before arcade formation events, and that of arcade flares discovered by Yohkoh-SXT. A dynamic model of arcade flares and arcade formation events based on this dark filament model will be discussed in the forthcoming Paper II of this series. K e y words: Sun: flares — Sun: magnetic field — Sun: prominences — Sun: X-rays 1. Introduction The magnetic structure surrounding a dark filament is of intrinsic relevance to the understanding of the physics of arcade flarings (energetic arcade flares in active regions, and large-scale weaker X-ray arcade formation events outside active regions). A model of mass suspension in a dipped magneticfield configuration has been discussed by Kippenhahn and Schluter (1957), and many authors have assumed that such a mass suspension could occur due to sagging of the field lines, even at the top of a convex magnetic field arching over the field-polarity reversal lines in the photosphere (see Tandberg-Hanssen 1974; Pickelner 1971; Priest 1989). A model for arcade-type flares (Carmichael 1964; Sturrock 1966; Hirayama 1974; Kopp, Pneuman 1976) was developed based on this type dark-filament model: When a dark filament is somehow destabilized and flies away, it drags the arcade field with it into interplanetary space. Magnetic reconnection occurs between the antiparallel fields at the leg parts of this stretched fieid. A reclosed arcade in the low corona and a detached c i o s e d loop field surrounding a flying blob will be formed as the result of magnetic reconnection. The energy released in the flare is interpreted as being the difference 'm the energies in the stretched and reclosed magnetic n e lds, a substantial part of which is liberated in the rec i o s e d arcade during the process. The magnetic energy [s converted into kinetic and thermal energies of the gas ejected downward from the reconnection site and the energy of high-energy particles accelerated in the reconnect ion process, and transported down to the chromosphere a i o n g the magnetic fields. The soft X -ray emitting mass in the reclosed arcade is considered to be chromospheric © Astronomical Society of Japan • Provided by the NASA Astrophysics Data System 554 Y. Uchida et al. mass "evaporated" by being given excess energy via heat conduction and/or electron bombardment from above. Some researchers (e.g., Shibata et al. 1995) refer to this model as the CSHKP model, taking the initials of the above authors. A difficulty with this model of dark filament came to light in the observations by Leroy et al. (1983). Using the Hanle effect to observe the magnetic field in dark filaments, they showed that the direction of the field perpendicular to the dark filament is opposite to the direction expected from the simple connection of the bipolar field below. A solution of this dilemma has been proposed, for example, by introducing the Kuperus and Raadu (1974) type model with a field configuration having an inverted circular sub-connection at the top of the convex arcade (cf., Priest 1989). However, no sufficiently persuasive physical explanation has been given for the creation and stability of such an intricate configuration. A more basic question about this classical model of arcade flares was raised (Uchida 1980) concerning the kinetic energy of the rising dark filament. In order to stretch the strong field (whose magnetic energy ultimately should explain the large energy of flares in the reconnection process) and eventually cut it open, the energy in the rise of the dark filament should be larger than the energy of the flare itself. [This point, which was the motivation for the work of Uchida and Jockers (1979), and Uchida (1980), was later formulated in a clearer form by Aly (1991).] It should be pointed out here that there would be no magnetic neutral sheet in the classical model if the process of cutting the field open could not actually take place due to the energy problem. This is a more severe dilemma. If the classical picture is correct, what requires a more serious investigation is the dark filament rise, not the flare itself, because the flare is then merely a repairing process of a more energetic break-up caused by the rise of the dark filament. Some researchers, accordingly, have tried to explain the dark filament as a high-stress structure having a highly sheared field (cf., Moore, Roumeliotis 1992). The dark filament in this type model is an extremely energetic entity, because if, for example, an initially current-free arcade with a 10 G field in a high latitude dark filament case, is sheared to produce a structure (a dark filament) with an aspect ratio of 10-100, the field strength in the long dark filament should become 10 2 - 3 G in the corona by a simple argument, and this is difficult to hold down. (This is a simplified argument. In actuality, the coronal part of the field may expand, but then it does not explain the thin partition-like shape of the dark filaments). In this connection, it should also be noted that there is no evidence for an actual very long and systematic antiparallel excursion of the photospheric footpoints along the field-polarity reversal line to consistently shear a simple arcade structure to produce an elongated dark filament [Vol. 51, with an aspect ratio of 10-100. Also, there is no reason for the photospheric motion to follow the polarity reversal line, if the field is thought to be a passive entity down in the photosphere. Comparing this type model with the observations, the greatest difficulty is that the observed dark filaments are an intrinsically passive entity easily shaken by the influence of some other disturbance, and, on the other hand, are quite stable as a whole until they are disturbed by some newly emerged magnetic patches appearing nearby. The highly sheared appearance is better interpreted as in our model discussed later in sections 2 and 4. In section 2, some explanation is given about a quadruple photospheric magnetic source model that one of us proposed years ago after being motivated by the difficulty of the classical models as mentioned above. 2. Quadruple Magnetic Source Model to Solve the Difficulty One suggestion that explains the elongated magnetic. structure along the "polarity reversal line" in a more natural way was proposed by Uchida (1980), based on Uchida and Jockers (1979), to avoid the energy paradox mentioned above, with a related proposal for a mechanism to produce arcade-type flares from such a configuration. Their proposal was that the magnetic structure involving a dark filament may be due to a quadruple array distribution of the magnetic field (elongated regions of +,—,+, — polarities side by side, referred to as A, B, C, and D) in the photosphere below. In that case, there appears a "neutral line" (the locus of the "neutral points", which will become a vertical "neutral sheet" if the photospheric footpoints are squeezed towards the central line) in the planar component of the magnetic field, B± (the field components in the xz-plane perpendicular to the central polarity-reversal line in the photosphere which is taken to be the y-axis), in the corona above the central polarity reversal line. We refer to this symbolically as the "neutral sheet" in the following, although it is actually a current sheet with the y-component of the field lying in it. The longitudinal component along the y-axis, Byi is dominant in the "neutral sheet", because B± is close to zero there. In our model, By, which may be relatively week initially, but will become of the same order of magnitude as B± at certain distance on both sides when squeezed, is what is suspending the dark filament gas lying in the thin, high and long vertical "neutral sheet", and is supported, in turn, by the gradient in B± a^ciiscussed later in section 4. j This model proposed by Uchida (1980) is a natural one for the suspension of dark-filament material, compared with the Kippenhahn-Schluter model. When the Kippenhahn-Schluter model is applied to the top of a convex arcade, the material initially loaded at the top of © Astronomical Society of Japan • Provided by the NASA Astrophysics Data System No. 4] Quadruple Source Model for Arcade Flarings a convex loop will readily slip down, and it would be hard to produce a dip at the top of the arcade. If one looks into Kippenhahn-Schluter's treatment, it would be readily seen that their treatment is for a concave field with current flowing on the side boundaries in an appropriate direction to support it, and does not apply to the convex loop top situation. The quadruple-array source model, in contrast, can sustain the dark filament mass in the thin neutral sheet, thus explaining the observed very thin vertical partitiontype structure of the dark filament. The quadruple source model, however, was not accepted too well when proposed by Uchida (1980) when there were no sufficiently precise X-ray observations to compare with it. Many years after our proposal (Uchida 1980), the quadruple source characteristics have attracted increasingly greater attention, and have been discussed by Dahlburg and Antiochos (1991), Demoulin, Priest (1993), Hundhausen, Low (1994), Antiochos et al. (1994), and Sciffer (1997). In the following, we first summarize what we obtained from the observation by Yohkoh in section 3, and describe our theoretical model in section 4. Comparisons between the observation and the model will be given in section 5. Finally, some discussion will be given in section 6. 3. Coronal Structure Surrounding a Dark Filament, or Features in the Pre-Event Region of Arcade Formation and Arcade Flares, Observed by Yohkoh It is not necessary to say that the information about the structure and its change of the magnetic field (represented by coronal X-ray loops) in the still faint preflare and initial stages is vital for clarifying the flare mechanism from the point of view of causality. But this information could not be obtained before Yohkoh. We first concentrate on the pre-event structure of X-ray arcade formation outside of the active regions in subsection 3.1, since the details are more clearly observable in those large-scale, slowly developing fainter version of arcade flares outside of the active regions, and then come to the case of arcade flares in active regions in subsection 3.2. 3.1. Large Scale Arcade Formation outside Active Regions McAllister et al. (1992) reported on a detailed observation of this type event obtained on 1991 September 28 by Yohkoh. The observed behavior, however, was quite different from what was expected from the classical model, and people thought that it might have been an exceptional event. It turned out, however, from analyses of several more events (Uchida 1996a, 1996b, 1996c; Uchida et al. 1999a; Fujisaki et al. 1999) that what was ob- 555 served in this event was common for arcade formation events. Namely, a bundle of bright threads (we referred to this as a "spine") rises together with the arcade of loops, sometime after the dark filament has disappeared. The "spine" is definitely not the locus of the reconnecting points in the reclosing process of the once-opened arcade, as suggested by the classical model, because the "spine" consists of multiple longitudinal threads coming up from below! The locus of the reconnecting points can not form multiple threads. The "spine" eventually balloons up on one end, and forms a large cusp-like shape. Considering that these X-ray arcade formations are larger scale fainter versions of arcade flares occurring outside active regions with a weaker magnetic field, we examined their pre-event structure in detail. An important result of this investigation reported separetely by Uchida et al. (1999b) is that the faint pre-event structure of these X-ray arcade formations were found to be "dual-arcades", one on respective sides of the "field polarity reversal line" in the photosphere (figure la). The inside legs of each arcade land in the domain of the other arcade in an intermingled way, and the dark-filament lies in the region in which the legs are crossing, quite unlike what was expected from the classical model. Upon examining the Kitt Peak magnetogram, rotated precisely to the time of the X-ray observation, we found that the "field polarity reversal line" is not at all a clearly definable line as assumed by previous reporters, but is a belt of mixed polarity regions extending to both sides of such an imaginary line (figure lb. Jack Harvey immediately endorsed what we found as true, and said he always noted it. Also see Uchida et al. 1999b). The assumption in the classical model that the parts closer to the polarity reversal lines are preferentially connected to the corresponding closer parts of the opposite polarity side, while the distant parts are preferentially connected to the distant parts of the opposite polarity side, and form a familiar simple arcade, does not actually hold. Distant parts are connected to the close parts of opposite polarity region, and vice versa, quite unlike the classical model. The observed coronal structure, "dual-arcades", whose inside legs cross, just looks like the flux separatrix surfaces seen in the quadruple source model proposed by Uchida and Jockers (1979) and Uchida (1980) (see figures 3 and 4 later). The situation with a considerable number of flux patches transported into the opposite polarity region can be approximated by four belt of regions with different field intensities; especially when the density of the transported opposite polarity flux exceeds that of the normal flux in that locality, the situation may be represented by belts of sources +,—,+,—, quadruple array sources, if approximated in 2.5 D. © Astronomical Society of Japan • Provided by the NASA Astrophysics Data System 556 Y. Fig. 1. Coronal structure surrounding a dark filament: (a) Faint structure-enhanced SXT image around a high-latitude dark filament on southern hemisphere (16:44 UT, 1992 January 18). The white curve is the locus of the top of the corresponding dark filament (16:40 UT, 1992 January 16) taken at BBSO, rotated to the exact time of the SXT observation by using a technique described in Uchida et al. (1999b). (b) Kitt Peak magnetogram of 16:16 UT, 1992 Jan 18 (not rotated because the time is close enough to the time of the SXT observation). The white curve is the same as in (a), but down-projected to the photosphere to compare with the photospheric magnetic field, (c) The relevant Ha dark filament observed at BBSO, 16:40 UT, 1992 January 16. The height of the top of the dark filament is determined as 1.2 x 10 5 km in the method of matching the down-projected white locus to the border of the opposite polarity fields in the photosphere which is not a well-defined line as presumed before (Uchida et al. 1999b). 3.2. Arcade Flares in Active Regions The first example of the case of an arcade flare in active regions suggesting quadrupolarity, was, to our surprise, the famous flare of 1992 February 21 (GOES class M2). This flare, seen at the east limb, showed us for the first time the clear structure of an arcade flare seen axis-on (Tsuneta et al. 1992). The first report said that Yohkoh confirmed what the magnetic reconnection model predicted. What the young Yohkoh colleagues had in mind in writing their reports was the CSHKP model mentioned in section 1 (the "reclosing of the once opened simple bipolar arcade" model), without knowing about the debates on the energy difficulty of that model. Their con- a et al. [Vol. 51, elusion was derived a bit too hastily, because the quadruple source model, an opponent of the CSHKP model, can also explain the observed behavior of the cusp, for example, that their foot point distance increases with time, etc., equally well. The difference can be found (1) in the structures above the top of the cusp, which are faint, and could not be seen easily before, but possibly can be seen by Yohkoh, and (2) in the energy aspect that we claimed to be fatal, but did not attract much general attention when first claimed. Those could not be made clear due to the absence of precise enough observations, and we noted that Yohkoh would be able to contribute essentially to this point. We therefore pursued observational evidence that could discriminate those models by using Yohkoh-SXT. We tried to look into the faint pre-event phase of the 1992 February 21 event, and actually found that (a) there exist overlying structures connecting the top of the "pre-flare core" (which developed into the flare cusp later) back to the photosphere on both sides at several times of the width of the core structure (figure 2a). We also found (b) a bright structure existed near the axis of the dark tunnel, together with a vertical partition-like structure in the middle of the dark tunnel, below the pre-flare cusp (figure 2a). There was in this event a faint rising blob ejected from the region of the pre-flare core before the flare started. It is worth noting that this did not pull up and stretch the arcade, but rather, seemed to be pulled out by the rising "dual-loops" (back connections on both sides mentioned above) from the valley right above the pre-flare core (S. Morita et al. 1999, in preparation), and continued to be pulled up rather than making the loop into a sharply tipped shape by pulling it. All of these are not in agreement with what the classical model predicts, but in good agreement with what a quadruple source model suggests. Exactly similar features were found in another event, the 1991 December 2 flare of GOES class Ml.5 (Tsuneta 1993). Namely, there appeared clear quadruple connections in the flare site with the central X-shape (figure 2b). There was also a brightening right at the axis of the dark tunnel below the "candle flame" shaped cusp of the flare in this case. There was a rising blob from the valley part of the "dual arcades" type structure above the top of the cusp, leaving the X-shape behind, without pulling the loop into a sharp-tipped shape. These examples clearly showed similar structures connecting the top of the pre-flare core back to the photosphere on both sides, together with the low-lying activities near the axis of the dark tunnel, and a rising blob from the top part of the cusp ejected without pulling the loop into a sharply tipped shape. The classical model(s) have difficulties in providing the model counterparts for these newly discovered features, in addition to the energy paradox. We are convinced © Astronomical Society of Japan • Provided by the NASA Astrophysics Data System Quadruple Source Model for Arcade Flarings No. 4] 557 Fig. 2. Pre-flare structure above the core of arcade flares seen at the limb: (a) Time development of the 1992 February 21 flare. There are clear high loops connecting the top of the pre-flare core back to the photosphere on both sides, as well as a bright structure at the axis of the dark tunnel below the "candle flame" type flare cusp, that are not explicable in the classical models. A blob rises prior to the start of the flare, but it is being pulled up by the expanding "dual-loops", rather than pulling the arcade into a tipped shape as expected from the classical model, (b) The same of the 1991 December 2 flare. Corresponding features that can not be explained by the classical model, and supporting the quadruple source model, are seen clearly. that the situations for arcade flares and arcade formations are better represented by our quadruple magnetic source model to be shown in section 4. 4. Magnet ohydrost at ic Configuration due to Quadruple Array Sources in the Photosphere (2.5 D Model) We next discuss a 2.5 D solution of magnetohydrostatics in order to look into some basic properties of a model having a neutral point in B±, in order to better understand what occur in flares. We deal with the equation of magnetohydrostatics, 1 . j x B -gmdp + pg = 0, (1) where B is the magnetic field, j = (c/47r)rot B is the current flowing in the plasma, p is the pressure and p is the density of the plasma, and g is the gravity directed downward to the — z direction. B is assumed to have the form B = dz'"y' dx)' (2) where ip is the magnetic stream function in the xz-plane, and this, as well as Byj the field component in the ydirection (taken along the polarity-reversal line in the photospheric plane), is assumed to be independent of y (so-called 2.5 D approximation). Equation (1) is then written as 'dtp dx 2 ±<P + J_ /dipdBy An \ dx dz l^\ 2 dx I + dP dx 0, d(p dBy dz dx dp 1- ^vi^ + i ^ + ~d~z 4TT I dx ^ 2 dz (3) (4) -pa, (5) where V^_ = d2/dx2 + d2/dz2. Equation (4) can be written as B± - \7±By — 0, where B± — (BX,BZ), and indicates that By is constant along the field line JB_L, or By is constant along the curve ip =constant in the xzplane. By assuming that the temperature T is constant and that p and p depend on (p and z in a separable way like © Astronomical Society of Japan • Provided by the NASA Astrophysics Data System H 2a(z + l) [x2 + (z + l)2Y i p = f(<p)F(z), then the equations (3) and (5) reduce respec5 tively to S [Vol. 51, Y. Uchida et al. 5 5 8 / I 2 ^ 1 dB2y 2 dip Equation (9) is accordingly transformed into ox d<p ox (6) and V# w%% + «rfi; + ™- (7) If we further assume that the dependence of p on z fulfills a hydrostatic distribution with a scaleheight z s = ?RT/g where 5ft is the gas constant, or dlnF dz (14) 1 zs' 2a 2 d V| c ^ + [xi2 + (C - a) 2 ] 2 d<p [£>)] 4a(C - a) + l)/zt +8-rrp0((p)exp 2 £ + (C - a)2 0, (15) where V | c = d2/d£2 + d2/dC,2 . We solved this transformed equation by extending the Jockers'(1977) numerical code, which is a three-level sc heme of relaxation by Marder and Weitzner (1970). As f o r t h e boundary conditions, Jockers (1976, 1977) chose the following for a simple bipolar array situation: a 9 U=o= x 2 + l , ( a > 0 ) , (16) we have from equations (6) and (7) 2 2 V ^ + \^[B V{<P) + 87rp 0 (^)e- z/z '] = 0, (9) which is a two-dimensional Poisson's equation with a source distribution that depends on ip itself. We solve this equation by giving the distribution of ip on the boundary, together with the dependence of By and po on (p. Jockers (1977) had developed a numerical scheme to solve the equation for the force-free field case for which equation (1) simplifies itself into ( l / c ) j X B — 0 with a simple arcade-like geometry. Uchida and Jockers (1979) extended the scheme and took into account the pressure gradient force and gravity terms as in equations (1) through (9), and further, dealt with a different type boundary distribution of the field (quadruple array) for which a neutral sheet structure in B± appears in the corona above. By a conformal transformation from the half plane xz to the inside of a circle on the ££-plane by the use of a complex mapping function, W = ai^A, Z + i' (10) where Z x + iz, (11) and W = £ + iC. (12) The photospheric boundary in the xz-plane is mapped to a circle of radius a with the origin mapped to (£ = 0, ( — —a) and the infinity to (£ = 0, ( = a) by means of equation (10) through the relations C _ 2ax " [X2 + (Z + l) 2 ] ' (13) 3v(w) 'dip 8TT = -P<pn (/?,n>0), (17) in his calculation of a force-free field satisfying ( l / c ) j X B = 0, in place of our (1). Equation (16) has a positivesloped region in x < 0 and a negative-sloped region in x > 0, implying that he considered extended negative and positive polarity regions lying side by side divided by x — 0. He treated a sheared simple arcade. In the present paper in which we deal with the nonforce-free field having a neutral sheet structure in JB_L, we use a more general boundary distribution for ip (x, z — 0), which results in a neutral sheet structure in JB_L, and takes a dependence of j y on ip more suitable to our situation corresponding to the physical circumstance of the dark filament to be considered. We consider a magnetoplasma system in equilibrium which fulfills equation (1). We must note that the time history of the system before attaining the equilibrium under consideration is relevant in specifying the boundary condition on ip and jy(ip). In our case, we consider two systems of closed magnetic flux initially independent from each other like two pairs of extended magnetic regions (for the case of a quiescent dark filament), or two pairs of sunspots or two sets of spot-plage combinations (for the case of an active region dark filament), and assume that these are brought into contact by a systematic motion of the photospheric footpoints of lines of force approaching each other (Tang 1987). Corresponding to this process, we assume that initially *.<*-<»~S z=0 =£*o, exp (x - x 0 J 2 (18) and thus <P £^Bb,{Erf (x - x0i) + 1}, © Astronomical Society of Japan • Provided by the NASA Astrophysics Data System (19) No. 4] Quadruple Source Model for Arcade Flarings where B^ is the local peak of the field in each of the two-dimensional magnetic belts at x — xoi with a Gaussian widths c^, and Erf(x) is the Gauss' error function, Erf(x)= (2/V5F) /0* exp(-t2)dt. We assume that the pairs of magnetic fluxes were initially self-closed within them and that the net total flux vanishes at ±00. The condition is ip (x — +00, z = 0) = 0, or ] P ViBoi (20) 0 with (p(x = - 0 0 , z = 0) = 0 automatically fulfilled. A simplest example fulfilling equation (20) is a pair of extended field regions antisymmetric with respect to the central field polarity-reversal line x = 0 in the photosphere. By giving the boundary values for ip(x) on the lower boundary, we can solve equation (15) numerically by a relaxation method, as mentioned above, if we specify po and By as functions of <p, as in the next section. 5. A Solution and Its Characteristic Features to Compare with Observations If we consider a case with p0 and By independent of ip, the source term in equations (9) and (15) is identically zero and the equations reduce to the two-dimensional Laplace's equation. The solution is then a potential field for the given distribution of the sources on the boundary, equation (18) for our quadruple sources. We, however, consider here the source term for equations (9) and (15), for example, with P o M = Po c exp (<P ~ Vc) 21 + P00 (21) distributed around the critical lines of force, <p — tpc, passing through the neutral point in B±. A constant background pressure does not affect the solution. As for the dependence of By on y?, we also assume a function of the form By(ip) =BVcexp (<P - Vcf (22) with a2B — 2a2. This may well be the case if the weak background field with plasma which pre-existed between the flux regions is squeezed together with the plasma frozen to it by the approaching motion of the two systems towards each other, and some part of it remains confined in the intervening region within the sheet-like structure. In order to reproduce the process of the squeezing of the two field systems, we here move the distribution of the footpoints of the potential field lines towards x — 0, for example, by dtanh f — J , (23) 559 with appropriate values of d and h. The magnetic lines of force in initially different regions can not penetrate into each other's region if there are plasmas frozen-in to them. Therefore, a sheet current is formed on the interface of these regions by this squeezing in order to maintain the identity of the field lines in each region. An example of the solution for the case in which the photospheric footpoints are moved from x to x' is given. Figure 3a shows the field with the thus-formed current sheet, and figure 3b is the same one, but seen from a different angle in order to show the 3 D structures of the field. In this squeezed state, energy is stored in the distorted magnetic field, and we propose this to be the source of energy to be released in the process of flaring when the relaxation to a lower energy configuration is allowed via magnetic reconnection at the "neutral points". The first thing to note from figure 3b is a pattern similar to the characteristic pattern of the observed dark filaments as seen from above (Martin, Levi 1992; Martin et al. 1994). This pattern was commonly attributed to the sheared arcade (the By field produced from a simple arcade by a shear due to a footpoint motion to a large aspect ratio, however, should be unreasonably strong, and also no observational evidence for such a shearing motion exists as noted in section 1). It is now quite naturally explained by a reasonable By imbedded along the neutral sheet of B± in our model. In the neutral sheet in JB_L, even a By of a modest strength can dominate, and the field pattern is elongated into the y-direction without actually shearing the structure. By in the thin neutral sheet may be considered to be the original weak field which pre-existed in the intervening region, strengthened by being squeezed between the two approaching field regions, and has a considerable value to render magnetohydrostatic equilibrium together with B± and the pressure of the dark-filament gas. The plasma may initially be the plasma frozen-in to the squeezed field, but it will be supplemented by a plasma flowing into the region by a mass-supply mechanism, which will be mentioned later. Secondly, the thin partition-like structure of the Ha dark filament, which is difficult to explain in the previous models, because of rather slight sagging of strong fields due to the relatively small gravitational force on the dark filament material (Low 1975) for the classical model, is now naturally explicable because the thin partition-like structure is basically the vertical shape of the squeezed neutral sheet in B± in our model, as can be seen in the lower-middle of figure 3a. If we look at the field structure from the sides, it is seen to be nested in this thin sheet, because the field at the center of the sheet is horizontal, while the field in the region on both sides are gradualy tilted toward the vertical with distance according to the nature of the solution. The line-of-sight field component in the dark filament in this case is nicely opposite from what the polarities of the "bipolar regions" predict. Also, © Astronomical Society of Japan • Provided by the NASA Astrophysics Data System 560 Y. Uchida et al. [Vol. 51, (a) Fig. 3. A solution of magnetohydrostatic equation (15): (a) A mildly squeezed solution described in the text, seen almost axis-on. An X-point-like structure is seen near the center, and the magnetic fluxes from each one of the sources are (i) fluxes connecting the source A to B and D, (ii) fluxes connecting the source B to C and A, (iii) fluxes connecting the source C to B and D, and (iv) fluxes connecting the source D to C and A. The part in the low central part will show up as an X-ray cusp when the reconnect ion takes place at the "neutral point" right above it, and various characteristic behavior, like the widening of the footpoints of the cusped arcade etc., take place just as observed. The difference from the classical model is that the strong part of the field is already "opened" by the effect of the second pair of magnetic fields, (b) A bird's-eye view of the same solution. The B y component dominates in the neutral region in -Bj_, and the structure appear as if highly sheared in the middle into which the lines of force coming from A and D flow into, and then leave from it after running some distance along the valley where the neutral line in B± is located. (In 3 D, B y is likely due to the presence of a flux pair at skewed positions outside of the finite region in which the elongated inner pair of sources exist.) there can be vertically flattened helices in the neutral sheet region corresponding to the "island" between the oppositely directed fields on both sides, explaining the observed helical field in the eruptive prominences. Thirdly, our model has simply-connecting high loops over this neutral sheet region of B±. The presence of such overlying loops are noted in soft X-ray observations (Vaiana et al. 1973), and thought of as evidence for the simply-connected configuration in favor of the classical picture. It is interesting to note that in some cases the overlying loops are not affected by the disappearance of the dark filament (Serio et al. CfA Preprint No. 935, 1978). In our present picture, the mass in the dark filament has been pumped in along the field lines which connect the dark filament to the hot coronal regions on both sides in skew (see figure 3b) by the "syphon" mechanism (Meyer, Schmidt 1969) due to the decrease of the pressure caused by the radiative cooling in the high density dark filament. The direction of the flow can be inverted without interfering the overlying loops if heating, rather than cooling, takes place in the dark filament. If, however, the heating is too violent, or some mechanism of eruption operates, for example, due to some magnetic instability, as in large flares exerting a force on the dark filament in the perpendicular direction to By, the dark filament as well as the overlying loops may also be disturbed, and pushed upward. We will present in a following paper a magnetohydrodynamic treatment of the 2.5 D dark filament model by dealing with the magnetodynamical simulation of the processes following a trigger destabilizing the system (Hirose et al. 1999 in preparation). The expansion of such large scale loops is actually observed by Yohkoh, but a noticeable point is that they are not made to be sharp-tipped as they would be if pulled by the rising dark filament as in the classical model, but rise without changing their initial rounded shape of the overlying loops. A simulation of the quadruple source model shows the loop behavior just as observed (Hirose et al. 1999 in preparation), and provides an inclusive interpretation of the phenomena associated with flarings much better. Fourth, we have a longitudinal field lying along the neutral sheet (mainly By) supporting the radiativelycooled dense plasma held down by gravity; this massloaded By) in turn, is supported in the gradient in the structure of JB_L- Our configuration is equipped with a structure supporting these mass-loaded field lines in the sheet from below. This is an array of small-scale loop structures lying below, directly bridging both sides of the "central polarity reversal line" in the photosphere, as can be seen in figure 3a (strong field connecting the region B to C). Such a structure is actually observed in the low corona, and was often erroneously thought to © Astronomical Society of Japan • Provided by the NASA Astrophysics Data System Quadruple Source Model for Arcade Flarings evidence against the existence of a physical neutral structure in association with the "polarity reversal . The neutral sheet in B± (with By lying in it) exght above this in the corona in our model. The lowlying structure plays an important role in supporting the mass-loaded field lines lying inside the neutral sheet by providing the gradient in B±, providing the model counterpart of "barbs" (Martin et al. 1994), whereas none of the previous models supporting the gas above the convex arcade could explain the "barbs". The lower part of the dark filament is considered to show up as a heated S-shaped feature pressed down on the structure below, while the upper part of it is squeezed out and flies away as an erupting dark filament before arcade flarings. The major antiparallel vertical fields can come into contact and are allowed to reconnect after the dark filament with By is expelled out. It will be shown in Paper II that the presence of the dark filament containing By prevents the antiparallel field from contacting, and stabilizes the system before the dark filament is squeezed out. Anomalous resistivity can only set in when the current density that drives plasma turbulence can get high after the antiparallel fields can be squeezed together after the dark filament is gone. Finally, it should be mentioned that the structure in the central part of figure 3a below the "neutral sheet" shows up as an X-ray cusp when the reconnection takes place in the "neutral point" right above it, and various characteristic behavior like the widening of the footpoints of the cusped arcade etc., take place just as observed. The essential difference from the classical model is (a) that the strong part of the field is already "opened" by the effect of the second pair of magnetic field, removing the difficulty in energy in our model, while there does not exist any neutral sheet in the classical model unless the dark filament opens up the strong part of the closed magnetic arcade. This causes the energy difficulty. And (b) the actual structure above the "neutral sheet" revealed by the observations by Yohkoh, favoring our model, as mentioned in section 3. All the features from the quadruple source model summarized above seem to explain the characteristics of the observed dark filaments, as well as the characteristic preevent X-ray coronal structures very well. 6. Conclusion and Discussion We have revisited the quadruple magnetic source model for the magnetic field surrounding dark filaments (Uchida, Jockers 1979). The model was not well accepted at the time of its proposal, and was not published. Some figures and explanations, however, are found in the Proceedings of the Skylab Workshop as a part of the Reports of the Working Groups, reported by Uchida (1980). It is now revived because the results of the new observations 561 of the pre-event coronal structure around the dark filaments from Yohkoh seem to suggest that quadruple photospheric magnetic source model explains the observation quite well (Uchida et al. 1994; Uchida 1996; Uchida et al. 1999a). The long-conceived bipolar arcade model of the dark filament suspension, which provided the basis of the " reclosing of the once opened simple bipolar arcade" model (so-called CSHKP model) for the arcade type flarings, not only have a difficulty with respect to energy, but can not provide model counterparts for various new features discovered by Yohkoh. 6.1. Separatric Surfaces for a More Realistic Case with Magnetic Patches Exchanged to the Opposite Polarity Regions We claim that what are observed in X-rays in the preevent state (figure la) are the flux separatrix surfaces in the quadruple source model. This is plausible because the contact line of the two branches of the separatrix surfaces is the magnetic neutral line in that case, and it is quite possible that the neutral line provides the separatrix surfaces with some mass and heat already in the pre-event phase of the flaring event that starts there one or two days later. We first show a set of two separatrix surfaces for the quadruple array magnetic sources in the photosphere, for the cases of a potential field in 2.5 D, as a reference, in figure 4a. If we name the four sources A, B, C, and D as before, it is seen that the near-critical lines of force coming from source A, for example, go either to source B or source D, but together with the part of other critical lines connecting D to C, they seem to land at source C of the same polarity as A. This explains the seemingly paradoxical fact that the observed pre-event loops in the "dual-arcades" seem to connect a region to the patches of the same polarity on the other side (precise comparison of figures la and b suggest this), though not very clear, due to the lowering of the temperature of the loop near the footpoints. Although the elongated field in the ydirection will give an impression of a highly sheared field, we now know that this is due to the dominance of By in the "neutral sheet" of B±, and no actual mechanical shearing is necessary. Next, in a more realistic situation, in which the exchanged patches of one polarity are transported into the opposite side, and vice versa, are taken into account, the magnetic configuration is no longer 2.5 D but 3 D. The critical lines of force in such a case in the potential field state are shown in figure 4b. Here, we can see in figure 4b that the configuration explains even better the observed features of "dual-arcades with crossed legs", seemingly connecting to the same polarity parts, as seen in figure la. © Astronomical Society of Japan • Provided by the NASA Astrophysics Data System 562 [Vol. 51, Y. Uchida et al. (b) ^^^^^^^fe^^j ifflm / /1' "^flM ^^^WjF^ & ^''S^TBH Sni^l Fig. 4. Flux separatrix surfaces for 2.5 D and 3 D cases: (a) The separatrix surfaces are in the form of "dual-arcades" shape as seen from figures 3a and 3b with critical lines of force spanning the surfaces. The contact line of the two surfaces is the neutral line (figures 4a,b are the potential field cases, for simplicity), (b) Critical lines of force for a case with exchanged patches of oppsite polarity sources. These critical lines of force, corresponding to the critical lines in figure 3b, seemingly connecting source A to C and B to D spanning on the separatrix surfaces, represent the observed "dual-arcades" quite well (cf., Uchida et al. 1999a). 6.2. Relation of Our Global Model to the Local X-Point Reconnection Models Finally, we would like to point out that the local treatment of the reconnection problem, started by Sweet (1958) and treated by many authors, has not really been solved yet. The huge discrepancy in the time scales between the rise time of flares, on the order of 102 s, and the magnetic diffusion time scale, on the order of 10 12 s, has annoyed researchers, and various kinds of trials have been made to resolve the discrepancy without any real success. For example, an introduction of anomalous resistivity, or an enhancement of the resistivity by an order of 10 6 , still leaves a 4-5 orders of magnitude discrepancy. In the model proposed by Uchida (1980), the rapid flare rise is attributed to a rapid dynamical collapse to the intervening lower energy interchanged state through an interchange instability after the longitudinal field (By with dark filament) which stabilized the configuration against this instability in the neutral sheet with some asymmetry in general, is lost with the dark filament rises. The long-endured energy release in the later phase of arcade flares for hours to one day has been identified with dissipation enhanced through an increase in the contact surface of the opposite-polarity fields by a large factor as the result of this interchange collapse, thus producing a highly interleaved state of opposite polarity fields. One of the most important findings from the dynamical simulations is that the energy release comes from a relax- ation of the magnetic stress all over the squeezed magnetic field below the separatrix surfaces. This relaxation is allowed by the occurrence of magnetic reconnection at the "neutral points". The release of the magnetic stress and mass stored within this structure contained below the separatrix surfaces can exceed by considerable factors the energy and mass directly released in the locality of the critical point. This may solve the riddle that the total release of mass and energy in CME's (corresponding to the mass and energy released dynamically from the stressed region under the separatrix surfaces in our model) exceeds those liberated in the flare itself (corresponding to the direct release of the mass and energy in the locality of reconnection site), and they seem to provide a consistent explanations for this. Actual results of our MHD simulations of the dynamical energy and mass releases will be discussed in Paper II of the series (Hirose et al. 1999 in preparation). We acknowledge that this work was made possible through the results of observations by Yohkoh, and thank all the team members who made the successful satellite possible through hardware and software construction, continued successful operation of the satellite, and the collaborative data analysis and scientific discussions. References Aly J.J. 1991, ApJ 375, L61 © Astronomical Society of Japan • Provided by the NASA Astrophysics Data System No. 4] Quadruple Source Model for Arcade Flarings Antiochos S.K Dahlburg R.B., Klimchuk J.A. 1994, ApJ 420, L41 Carmichael H. 1964, in Proc. NASA Symp. on the Physics of Solar Flares, ed W.N. Hess (NASA SP-50) p451 Dahlburg R.B., Antiochos S. 1991, ApJ 383, 420 Demoulin P., Priest E.R. 1993, Sol. Phys. 144, 283 Pujisaki K., Uchida Y., Morita S., Tsuneta S., Hirose S., Cable S. 1999, PASJ submitted Hirayama T. 1974, Sol. Phys. 34, 323 Hundhausen A., Low B.C. 1994, ApJ 429, 876 Jockers K. 1976, Sol. Phys. 50, 405 Jockers K. 1977, Sol. Phys. 56, 35 Kippenhahn R., Schluter A. 1957, Z. Astrophys 43, 36 Kopp R.A., Pneuman G.W. 1976, Sol. Phys. 50 85 Kuperus M., Raadu M. 1974, A&A 31, 189 Leroy J.-L., Bommier V., Sahal-Brechot S. 1983, Sol. Phys 83, 135 Low B.C. 1975, ApJ 198, 211 Marder M.B., Weizner H . 1970, Plasma Phys. 12, 435 Martin S.F., Bilimoria R , Tracadas P.W. 1994, in Solar Surface Magnetism, ed R.J. Rut ten, C.J Schrijver (Springer, New York) p303 Martin S., Levi S. 1992, in Eruptive Solar Flares, ed Z. Svestka, B. Jackson, M. Machado (Springer, Berlin) p33 McAllister A.H., Uchida Y., Tsuneta S., Strong K., Acton L., Hiei E., Bruner E., Watanabe Ta. et al. 1992, PASJ 44, L205 Meyer F., Schmidt H.U. 1968, Z. Angew. Math. Mech. 48, 218 Moore R., Roumeliotis G. 1992, in Eruptive Solar Flares, ed Z. Svestka, B. Jackson, M. Machado (Springer, Berlin) p69 Pickelner S.B. 1971, Sol. Phys. 17, 44 Priest E.R. 1989, Dynamics and Structure of Quiescent Solar Prominences (Kluwer, Dordrecht) Sciffer M. 1997, Thesis, University of New Castle, Australia Shibata K., Masuda S., Shimojo M., Hara H., Yokoyama H., Tsuneta S., Kosugi T., Ogawara Y. 1995, ApJ 451, L83 563 Sturrock P A . 1966, Nature 221, 695 Sweet P A . 1958, in Electromagnetic Phenomena in Cosmic physics, ed B. Lehnert (Cambridge University Press) pi23 Tandberg-Hanssen E. 1974, Solar Prominences (D. Reidel, Dordrecht) Tang F. 1987, Sol. Phys. 107, 203 Tsuneta S. 1993, in The Magnetic and Velocity Fields of Solar Active regions, ed H. Zirin, G.-X. Ai, H. Wang, ASP Conf. Ser. 46, p239 Tsuneta S., Hara H., Shimizu T., Acton L.W., Strong K.T., Hudson H.S., Ogawara Y. 1992, PASJ 44, L63 Uchida Y. 1980, contributions by Uchida, printed in the Reports of Working Groups, in Sky lab Workshop, Solar Flares, ed P.A. Sturrock (University of Colorado Press) pp67, 110 Uchida Y. 1996a, Adv. Space Res. 17, 19 Uchida Y. 1996b, in The Structure of the Sun, ed T.R. Cortes, F. Sanchez (Cambridge University Press) p355 Uchida Y. 1996c, in Magnetodynamic Phenomena in the Solar Atmosphere-Prototypes of Stellar Magnetic Activity, ed Y. Uchida, T. Kosugi, H.S. Hudson (Kluwer, Dordrecht) p295 Uchida Y., Fujisaki K., Morita S., Torii M., Hirose S., Cable S. 1999b, PASJ 51, 53 Uchida Y., Hirose S., Morita S., Fujisaki K., Torii M., Tanaka T., Yabiku T., Miyagoshi T., Uemura S., Yamaguchi T. 1999a, Ap&SS in press Uchida Y., Jockers K. 1979, presented at the Skylab Workshop in 1980 (Uchida 1980) Uchida Y., McAllister A., Khan J., Sakurai T., Jockers K. 1994, in X-ray Solar Physics from Yohkoh, ed Y. Uchida, T. Watanabe, K. Shibata, H. Hudson (Universal Academy Press, Tokyo) pl61 Vaiana G.S., Krieger A.S., Timothy A.F. 1973, Sol. Phys. 32, 81 Webb D.F., Krieger A.S., Rust D.M. 1976, Sol. Phys. 48, 159 © Astronomical Society of Japan • Provided by the NASA Astrophysics Data System
© Copyright 2026 Paperzz