Main Cause of the Poloidal Plasma Motion Inside a Magnetic Cloud Inferred from Multiple-Spacecraft Observations

Although the dynamical evolution of magnetic clouds (MCs) has been one of the foci of interplanetary physics for decades, only few studies focus on the internal properties of large-scale MCs. Recent work by Wang et al. (J. Geophys. Res.120, 1543, 2015) suggested the existence of the poloidal plasma motion in MCs. However, the main cause of this motion is not clear. In order to find it, we identify and reconstruct the MC observed by the Solar Terrestrial Relations Observatory (STEREO)-A, Wind, and STEREO-B spacecraft during 19 – 20 November 2007 with the aid of the velocity-modified cylindrical force-free flux-rope model. We analyze the plasma velocity in the plane perpendicular to the MC axis. It is found that there was evident poloidal motion at Wind and STEREO-B, but this was not clear at STEREO-A, which suggests a local cause rather than a global cause for the poloidal plasma motion inside the MC. The rotational directions of the solar wind and MC plasma at the two sides of the MC boundary are found to be consistent, and the values of the rotational speeds of the solar wind and MC plasma at the three spacecraft show a rough correlation. All of these results illustrate that the interaction with ambient solar wind through viscosity might be one of the local causes of the poloidal motion. Additionally, we propose another possible local cause: the existence of a pressure gradient in the MC. The significant difference in the total pressure at the three spacecraft suggests that this speculation is perhaps correct.


Introduction
Magnetic clouds (MCs) as a considerable subset of interplanetary manifestations of coronal mass ejections (CMEs) have been studied for many years, since the first identification by Burlaga et al. (1981). The main reasons that we pay so much attention to MCs are that MCs have a relatively well-defined magnetic topology and are the main drivers for many spaceweather events (e.g. Tsurutani et al., 1988;Huttunen, Koskinen, and Schwenn, 2002;Wu and Lepping, 2002;Cane and Richardson, 2003;Zhang et al., 2007). The dynamical evolution of MCs in interplanetary space is one of the important topics that need to be determined to improve the precision of space-weather forecasts.
Abundant observational evidence has shown that MCs may experience expansion (Burlaga and Behannon, 1982;Klein and Burlaga, 1982;Farrugia et al., 1992Farrugia et al., , 1993Lepping et al., 2002;Wang, Du, and Richardson, 2005;Gulisano et al., 2010), latitudinal and longitudinal deflection (Wang et al., 2002(Wang et al., , 2004(Wang et al., , 2006a(Wang et al., , 2014Kilpua et al., 2009;Rodriguez et al., 2011;Isavnin, Vourlidas, and Kilpua, 2013), and/or rotation (Wang et al., 2006b;Yurchyshyn, 2008;Yurchyshyn, Abramenko, and Tripathi, 2009;Vourlidas et al., 2011;Nieves-Chinchilla et al., 2012) when they propagate away from the Sun in the interplanetary medium. In addition to these motions, the "poloidal plasma motion" here refers to the motion of a velocity component along the poloidal direction, i.e., around the MC axis, that has recently been found in the moving MC frame for many MCs by Wang et al. (2015). For this motion, the authors raised three possible explanations: i) it might be locally generated through the interaction with the solar wind; ii) it might be internally generated by the expansion of MCs, during which magnetic energy is converted into kinetic energy, including the rotational component; and iii) it might be initially generated at the eruption of the corresponding CME, and carried all the way to 1 AU.
In order to determine the main cause of the poloidal plasma motion in MCs, we employ multi-spacecraft observations to study the same MC. The Solar Terrestrial Relations Observatory (STEREO) has twin spacecraft, one (STEREO-A) leading and the other (STEREO-B) trailing the Earth, and the separation between them increases by about 45 • per year (Kaiser et al., 2008). The in-situ data from STEREO-A (ST-A), STEREO-B (ST-B), and the Wind spacecraft, which is at the Sun-Earth L 1 point, are often used to perform multi-spacecraft analyses. Nevertheless, the number of MCs that are suitable for a multispacecraft analysis is very limited, because when the distance between two spacecraft is too small (e.g. Wind and the Advanced Composition Explorer (ACE)), the two spacecraft observe almost identical structures; and when the distance between two spacecraft is too great, in addition to the finite longitudinal extent of the MC, which is about 60 • (Wang et al., 2011), the solar wind may also be strongly distorted, which means that it is difficult to recognize whether the structures belong to the same event. The separation angle between two STEREO spacecraft for a multipoint study is suggested to be no larger than 50 -60 • . We examined the MCs that occurred during February 2007 to June 2008, when the STEREO spacecraft were separated by about 1 • to 58 • , and we found that the best case for our study is the MC on 19 November 2007, when the spacecraft were separated by 41 • .
In the next section, we identify and reconstruct the MC at the three spacecraft ST-A, Wind, and ST-B. In Section 3, we show evidence of the plasma poloidal motion in this MC and discuss the possible main cause of this motion. A conclusion is reached in Section 4.

Data Source
The positions of the spacecraft at the occurrence of this event are shown in Figure 1. For ST-A and ST-B, the in-situ magnetic-field and solar-wind plasma data are obtained from the In situ Measurements of Particles And CME Transients (IMPACT: Luhmann et al., 2008) and Plasma and Suprathermal Ion Composition (PLASTIC: Galvin et al., 2008) instruments, and these are shown in Figures 2 and 4, respectively. Meanwhile, for Wind, the in-situ magnetic field and solar-wind plasma data from the Magnetic Field Investigation (MFI: Lepping et al., 1995) and Solar Wind Experiment (SWE: Ogilvie et al., 1995) instruments are shown in Figure 3. All magnetic-field strengths and plasma properties are plotted in one-minute temporal cadences, and the charge states of iron ions are averaged for two hours. The data of the average iron-ion charge state from the Solar Wind Ion Composition Spectrometer (SWICS: Gloeckler et al., 1998) onboard ACE are used as a complement to the Wind data, because there is no instrument to measure the charge state of iron ions onboard Wind.

Identification of the MC at Three Spacecraft
This event was studied by Farrugia et al. (2011), Kilpua et al. (2011), andRuffenach et al. (2012). The MC has different features than the ambient solar-wind properties when a spacecraft passes through the MC, including i) an enhanced magnetic field, ii) a smooth rotation over a large angle in the direction of the magnetic field, and iii) low proton β and temperature (Burlaga et al., 1981;Klein and Burlaga, 1982). According to these three criteria, the MC was identified to last from 19 November, 22:30 UT to 20 November, 20:00 UT at ST-A, 20 November, 00:20 -09:00 UT at Wind, and 19 November, 23:00 UT to 20 November, 07:00 of the magnetic-field orientation, three components of the bulk velocity in RTN coordinate system, proton density [N p ], proton temperature [T p ], proton β p , total pressure [P t ], and the average charge state of iron ions [F e ]. The ratio of the proton temperature to the expected proton temperature [T p /T exp ] that is calculated based on the empirical formula from Lopez and Freeman (1986) is shown in the eighth panel by the red curve.
UT at ST-B in Farrugia et al. (2011). By comparing with the data at the three spacecraft, the authors found that the data at ST-B and Wind show a gradual transition from the slow to the fast solar wind, and ST-A may be located at the trailing edge of the previous high-speed    Ruffenach et al., 2012).
By combining the above three classic signatures of an MC and the features in the density and iron charge state, we make our own identification of the boundaries of the MC at the three spacecraft. Our identified boundaries of the MC at ST-A and Wind are quite similar to those by Farrugia et al. (2011). At ST-A, the MC in our view started at 19 November, 22:25 UT and ended at 19:00 UT on the next day (see Figure 2), which is five minutes and one hour earlier than the leading and trailing edge identified by Farrugia et al. (2011). At Wind, the MC we identified started on 20 November, 00:30 UT and ended at 09:20 UT (see Figure 3), which is 10 minutes and 20 minutes later than the leading and trailing edge identified by Farrugia et al. (2011). At ST-B, the MC boundaries suggested by Farrugia et al. (2011) are marked by the first solid line and the dotted line (see Figure 4). There is a highdensity region in the middle of the MC, which was regarded as filament material contained within the MC (Farrugia et al., 2011). However, we note that a similar high-density region appeared after the trailing edge of the MC in both ST-A and Wind observations. We therefore propose here a different scenario: the MC at ST-B was located between the two solid lines, as shown in Figure 4, in which the high-density region followed the MC, just like those observed at ST-A and Wind. Moreover, the profile of the average charge state of iron ions at ST-B shows that the charge state first reached a peak inside the MC and then decreased to a low outside the MC in our identification. This is consistent with the measurements at both ST-A and Wind. Although we proposed a rear boundary different from that by Farrugia et al. (2011) for the MC at ST-B, it is difficult to judge which identification is more correct. We therefore consider both possibilities in the following analysis. Hereafter, we use ST-B for our identification and ST-B for the identification of Farrugia et al. (2011).
According to the velocity (about 450 km s −1 ) of the in-situ MC at Wind, the source CME can be roughly estimated to have launched on 15 November 2007. After searching in the online LASCO CME catalog (Yashiro et al., 2004), we found that there is no suitable CME from 14 to 16 November except for the faint partial halo CME that erupted around 18:50 UT on 15 November, as shown in Figure 5a. This CME can also be seen from ST-B, as shown in Figure 5b. We conclude that it is the source of the MC. However, the exact source region cannot be identified confidently because of the weak activity on the solar surface.

Reconstruction of the MCs
We used the velocity-modified cylindrical force-free flux-rope model that was developed by Wang et al. (2015) to fit the MC. For ST-A and ST-B, we show the measurements in the RTN coordinate system, in which R points from the Sun to the spacecraft, T is perpendicular to R and points in the direction of planetary motion, and N completes the right-hand coordinate system. For Wind, the GSE coordinate system was used. Meanwhile, a cylindrical coordinate system (r, ϕ, z) and a Cartesian coordinate system (x , y , z ) were used to present the data in the MC frame, which were illustrated in detail in Wang et al. (2015). The axis orientation of the flux-rope is indicated by the elevation and azimuthal angles θ and φ. When θ = 0, the axis is in the ecliptic plane, and θ = 0 and φ = 0 means that the axis points to the Sun in GSE coordinates or to the spacecraft in RTN coordinates. By considering the linear propagating motion [v c ], the expanding motion [v e ], and the poloidal motion [v p ], the model not only uses three components of the measured magnetic field, but also three components of the measured velocity to constrain the fitting parameters.
In the light of the boundary of the MC at three spacecraft that we identified, we use our model to fit the five-minute-averaged data of the MC. All of the fitting parameters obtained by the model are shown in Table 1, and the fitting results are indicated by the blue lines in Figures 2 -4. We also fit to the MC at ST-B according to the boundary of Farrugia et al. (2011) and show this as the red line in Figure 4.
Our fitting results of the MC at the three spacecraft are roughly consistent with those given by Farrugia et al. (2011). The largest difference occurs in the orientation of the MC axis. At Wind, the orientation difference is about 40 -50 • , while at ST-A and ST-B the difference is smaller than 20 • . The orientation of the MC at the three spacecraft is plotted in Figure 1. From this figure, it is clear that the axes of the MC at the three spacecraft roughly align along a loop, in agreement with the common picture of an MC that is a globally looplike structure with two ends rooted on the Sun.
From Table 1, the parameter d, i.e. the closest approach, at ST-B and Wind is about 0.75 -0.85, which means that the spacecraft crossed the edge of the MC, which perhaps affects the reliability of our model results. Riley et al. (2004) performed "blind tests" for a magnetohydrodynamically (MHD) simulated MC using five different fitting techniques, and they found that a large deviation for the orientation of the MC axis existed, especially when the closest approach of the spacecraft to the MC axis is large. However, in the three types of linear force-free models, the orientation of the MC axis from Lepping's model (Lepping, Burlaga, and Jones, 1990) is closest to the value that the MHD model gave, in which the difference is less than 60 • . Our model has been demonstrated statistically to yield similar fitting results as Lepping's linear force-free flux-rope model (Wang et al., 2015). In this study, we assume that the modeled orientations are close to the real orientations, and we investigate the poloidal motion of the MC plasma in the following sections.

Poloidal Motion Inside the MC at Different Spacecraft
After we obtained the orientation of the MC axis, we converted the magnetic field and velocity into the MC frame, (x , y , z ), in which z is along the axis of the MC that is perpendicular to the x -y plane. The thick black arrows in Figures 6d -8d show the paths along which the spacecraft pass through the MC projected on the x -y plane. Along the path from the beginning to the end of the MC, the observed x -and y -components of the magnetic field and plasma velocity are shown by color-coded dots in Figures 6a, b -8a, b. A clear and regular Table 1 The fitting parameters of the MC obtained by the velocity-modified cylindrical force-free flux-rope model.  arc-structure is found in the magnetic-field data (see Figure 6a -8a), confirming the helical magnetic structure in the MC. In Figures 7b -8b, we can also find an arc-like distribution of the velocity in the x -y plane, which indicates the poloidal motion of the plasma flow inside the MC. For the MC observed by Wind, this motion rotates anticlockwise, with the velocity changing from +v y to +v x . The signature of the poloidal motion for the MC observed by ST-B is not as clear as that observed by the Wind spacecraft, but it can still be recognized as an anticlockwise rotation in the x -y plane. For the MC observed by ST-A, there is no clear plasma rotation. We also analyzed the poloidal motion inside the MC at ST-B with the boundary identified by Farrugia et al. (2011) (Figure 9), and the poloidal motion was also very weak. The modeled overall poloidal velocities are listed in Table 1.
To further check the reliability of our model results, we used the equation where c p and c e are two free parameters, B t is the measured magnetic field in the x -y plane to fit the measured velocity [v t ] in the x -y plane. Here, we consider that the magnetic field is a well-organized helical structure and the plasma is frozen on the magnetic-field lines, therefore v t is decomposed into two components, i.e. poloidal motion and expansion, which are assumed to be parallel and perpendicular to the local magnetic-field line, respectively. Figures 6c -9c show the results. The horizontal axis indicates the fitted velocity, which is calculated from the above formula based on the measured B after c p and c e are obtained from the fitting to all the data points, and the vertical axis represents the measured velocity. The blue dots represent the v x -component of the velocity, and the red dots represent the v y -component. Significant correlations are found between the fitting and measured velocities at all the three spacecraft. The correlation coefficients are 0.6189, 0.6446, 0.6359, and 0.6959 at ST-A, Wind, ST-B, and ST-B , respectively. This suggests that the velocity model incorporated in our flux-rope forward-modeling is reasonable and acceptable, and the significant deviation might result from fluctuations in the measured magnetic field and velocity, and it reveals why this poloidal motion of the plasma inside MCs was not noted before.
We have given three hypotheses for the causes of the poloidal motion (Wang et al. 2015; see also the Introduction). Two of them are global reasons, which should cause a consistent poloidal motion at the three spacecraft. The above analysis has suggested anticlockwise rotation at Wind and ST-B, but no clear rotation at ST-A. This implies that the main cause of the poloidal motion cannot be global, but must be local. We note that even if we had replaced ST-B with ST-B , the conclusion would not be changed as the poloidal motion was quite different at the three spacecraft.
What type of local process may cause this poloidal motion inside the MC? One possible process is the interaction with the ambient solar wind. Viscosity is an intrinsic property of the solar wind (e.g. Whang, Liu, and Chang, 1966;Hundhausen, 1968;Hollweg, 1978;Bruno and Carbone, 2005). As long as there is a velocity difference between the MC plasma and the ambient solar wind at the boundary of the MC, the poloidal motion probably has a natural cause. To check this speculation, we first calculated the average solar-wind velocity every 15 minutes before the leading and after the trailing edge in the x -y plane. The values of the velocity component tangential to the modeled MC boundary (labeled v pswli and v pswti , i = 1, 2 and 3, with the smaller number indicating a place closer to the MC boundary) are listed in Table 2, in which the positive number means the anticlockwise rotation and the negative number the clockwise rotation. The uncertainty is the standard deviation of the speed in the 15-minute interval. Second, we calculated the average tangential plasma velocities within the MC every 15 minutes after the leading and before the trailing edge [v pl1 , v pl2 , v pl3 , v pt3 , v pt2 , and v pt1 ] which are also listed in Table 2.
The tangential velocities of the solar wind closest to the MC boundary are marked by the red arrows in Figures 6d -9d. From Figures 7d and 8d, we find that the rotational directions of the tangential velocities [v pswl1 and v pswt1 ] are consistent at Wind and ST-B, where the plasma poloidal motions inside the MC are significant. However, at ST-A the rotations of the ambient solar wind near the front and rear boundaries were opposite, as illustrated in Figure 6d. The rotations of the MC plasmas near the front and rear boundaries were also opposite (see Table 2). This is probably the reason why the modeled overall poloidal speed is not evident. We note, however, that the rotational direction inside the MC is consistent with the rotational direction of the ambient solar wind at both the front and rear boundaries. The situation at ST-B is similar. This apparent relationship between the rotational directions of the solar-wind plasma and the MC plasma qualitatively supports the idea that the interaction between the solar wind and the MC might be a cause of the poloidal plasma motion inside the MC.
To further quasi-quantitatively investigate the interaction between the solar wind and MC, we show the relation between the (v pswli + v pswti )/2 and (v pli + v pti )/2, where i is equal to 1, 2, and 3 in Figure 10. We use the sum of the tangential speeds at the front and real boundaries to reduce possible local bias and uncertainties and show an overall relation of the plasma velocities at the two sides of the MC boundary. The error bars in Figure 10 are the mean value of the errors of the velocities at the front and rear boundaries. From the left panel to the right panel, the data points used in the plot move farther from the boundary of the MC. These panels show that the tangential velocities at the two sides of the MC boundary are roughly correlated. Considering the possible uncertainty in identifying the boundary of the MC, the correlations presented in the plots are expected to become more and more reliable In particular, the middle and right panels show that the rotational speeds of the ambient solar wind are higher overall than those of the MC plasma at ST-A and Wind, suggesting that the viscosity may affect the poloidal motion. However, the data point at ST-B does not follow this pattern. Furthermore, the tangential velocity does not show a clear pattern that the value decreases from the edge of the MC toward the axis of the MC. These inconsistencies imply that the interaction with the solar wind through viscosity might be only one of the local causes of the poloidal motion. Here we propose another possible cause: a pressure gradient in the MC. An MC is a large-scale structure in interplanetary space. Its interaction with ambient solar wind cannot be uniform everywhere. In other words, the different parts of the MC may be affected by the solar wind to different degrees, and therefore the local pressures inside the MC may be non-uniform, which might also lead to local poloidal motion following the helical magnetic-field lines. However, this scenario is difficult to examine. One preliminary preparation to aid this verification is to check the total pressure at the three spacecraft, which is plotted in the second panels from the bottom of Figures 2 -4. They show that the total pressures in the MC at the three spacecraft are different, perhaps confirming this speculation.

Summary and Conclusion
In this article, multi-spacecraft observations were employed to study the possible main reason that causes the poloidal plasma motion inside the MC that was observed by ST We then checked the possible local cause, which we hypothesize is the interaction between the solar-wind plasma. By calculating the measured tangential velocities of the ambient solar wind and the MC plasma at the two sides of the MC boundary in the MC frame, we find that the rotational direction of the solar wind is consistent with the modeled overall poloidal motion. Furthermore, we compared the values of the tangential velocities, Table 2 The poloidal component of the averaged velocity in the x -y plane of the MC frame.
Spacecraft  (v pswli + v pswti )/2 and (v pli + v pti )/2 (i = 1, 2 and 3), at the three spacecraft and derived a rough correlation between them (see Figure 10). These results illustrate that the solar wind might be able to cause the poloidal motion in the MC through viscosity. The strength of the viscosity mainly depends on the gradient of the velocity and may be modified by the magnetic field in magnetized plasmas (e.g. Hollweg, 1985). The importance of viscosity in the solar wind has been demonstrated in the studies of CME's dynamical evolution, e.g. the aerodynamic drag force is a manifestation of the viscosity, which may dominant the dynamic evolution of ICMEs in interplanetary space (e.g. Subramanian, Lara, and Borgazzi, 2012;Vršnak et al., 2013;Hess and Zhang, 2015). Here our study suggests, from another point of view, that the viscosity plays a significant role in some circumstances. The question now is how the viscosity takes effect in a collisionless plasma. More work is necessary to address this issue.
However, in the above picture, the tangential speed of the solar wind should be higher than that inside the MC, which is roughly true at ST-A and Wind, but not at ST-B. Moreover, the viscosity should cause the tangential velocity to be higher at the edge than in the inner part of the MC. However, no such pattern is evident in our study. These inconsistencies imply that there are probably other local reasons causing the poloidal motion. By briefly checking the total pressure inside the MC at the three spacecraft, we propose that the nonuniform pressure in the large-scale MC structure might be another local cause.
We note that the conclusions obtained in this article are based on only one case, which needs to be verified by more events in future work.