Relativity and the Solar Wind : The Maxwell-Equation Origins of the Solar-Wind Motional Electric Field

The motional electric field of the solar wind as seen by the Earth is examined theoretically and with spacecraft measurements. As it flows outward from the sun, the solar-wind plasma carries a spatially structured magnetic field with it. To calculate the motional electric field of the solar wind the spatially structured magnetic field is Lorentz transformed; for a full physical understanding, it is also necessary to Lorentz transform the current densities and charge densities in the solar wind. Referring to Maxwell’s equations, two related questions are asked: 1) Is the source of the solarwind motional electric field charge density in the solar wind, time derivatives of current densities in the solar wind, or both? 2) Is the solar-wind motional electric field at Earth an electrostatic field, an induction field, or a superposition of the two? A Helmholtz decomposition of the motional electric field of the solar wind is made into a divergence-origin (electrostatic) and a curl-origin (induction) electric field. The global electric field associated with the outward advection of the global Parker-spiral magnetic field is found to be electrostatic with its origin being a distributed charge density in the solar-wind plasma. The electrostatic versus induction nature of the time-varying electric field associated with the advection of mesoscale magnetic structure varies with time as differently shaped magnetic structures in the solar-wind plasma pass the Earth; the mesoscale structure of the solar-wind plasma contains sheets of space charge and sheets wherein the current density has nonzero time derivatives.


Introduction
The solar wind is a magnetized collisionless plasma comprised of ions (mostly protons) and electrons.(Collisionless means that inter-particle collisions that would disrupt the ion or electron orbits in the magnetic and electric fields are rare.)The solar-wind plasma flows at a high velocity radially outward from the Sun, in the range 250 km/s to 1000 km/s.Typical parameters for the solar wind at Earth appear in Table 1.The magnetic field B in the solar-wind plasma originates at the Sun and is carried outward with the supersonic plasma flow.Because of the rotation of the Sun and the radial outflow of the wind, the magnetic field in the solar wind has on average a "Parker-spiral" orientation [1].The magnetic-field vector has strong localized variations about the Parker-spiral direction: these spatial variations in B are advected outward with the solar-wind flow [2].Much of the spatial variation of B is owed to the presence of current sheets embedded in the solar-wind plasma [3] [4], with the value of the magnetic-field vector B jumping across the current sheets.
The motional electric field E of a plasma is a Lorentz transformation of the magnetic field B in the plasma [5].For a velocity vector v between two reference frames, the Lorentz transformation of the components of the electric and magnetic fields parallel to and perpendicular to v is [6] E E′ = (1a) ( ) ( ) and where we will take the primed frame to be at rest in the solar-wind plasma and the unprimed frame to be the rest frame of the Earth or the Sun.In the rest frame of the solar-wind plasma, 0 E ⊥ ′ = .This is because (in the absence of non-electromagnetic forces such as gravity) in a collisionless plasma particle orbit theory says that the ions and the electrons E-cross-B drift with an averaged velocity of ( ) (cf. eq. ( 12.43) of Reference [6]).In the rest frame of the plasma there can be no E-cross-B drift velocity of the ions and electrons, otherwise it would not be the rest frame of the plasma.Therefore, in the rest frame of the solar-wind plasma 0 E ⊥ ′ ≡ .In the reference frame of the Earth or the Sun, the plasma flow is of the order of v/c ~ 10 −3 , so γ ≅ 1 can be taken.For 0 E ⊥ ′ = expression (1d) gives B B ⊥ ⊥ ′ = .For γ ≅ 1, 0 E ⊥ ′ = , and B B ⊥ ⊥ ′ = , expression (1b) yields the motional electric field E for a plasma moving with a velocity v in a magnetic field B to be ( ) which is the low-velocity Lorentz transformation common for a moving conducting fluid [6] [7].The solar-wind motional electric field as seen by the Earth (expression (2)) is often taken to be the driver of geomagnetic activity in the Earth's magnetosphere-ionosphere system [8] (but see Reference [9] for a contrary interpretation).
In this paper, we wish to gain an understanding of the solar-wind motional electric field from first principles.The motional electric field of the solar wind will be examined via the Lorentz transformation of the magnetic-field structure in the solar wind.Because of the spatial variations of the magnetic-field vector B, the motional electric field will have spatial variations.Via Maxwell's equations, the sources of the motional electric field will be sought: charge densities in the solar wind and time derivatives of current densities in the solar wind.The question will be asked whether the motional electric field of the solar wind is an electrostatic electric field, an inductive electric field, or both.We wish to gain an understanding of the solar-wind motional electric field from first principles.Identifying the electrostatic versus induction portions of electric fields within the Earth's magnetosphere has been of interest because the identification yields an understanding of the origin of the electric fields and a capability to add them to numerical models of the magnetosphere [10] [11].
To avoid confusion, the reader should note that there is also a non-motional electric field in the solar wind associated with the "interplanetary electrostatic potential" or "heliospheric electric potential" [12].This is a global parallel-to-B electric field associated with the differences in the collisionless transport of the solar-wind electrons and the solar-wind ions along the magnetic field.The sign of this interplanetary electric field is outward along the magnetic field lines, accelerating ions away from the Sun and retarding the outward electron flux.This interplanetary electrostatic potential is responsible for the observed statistical increase of the bulk speed of the solar wind v sw between 0.3 AU and 1 AU.This interplanetary electrostatic potential is not the subject of the present study.

Electrostatic Fields versus Induction Fields and Charge Density versus Time-Varying Current Density
According to the Helmholtz decomposition [13] [14], any vector field can be decomposed into the sum of an irrotational (curl-free) field plus a solenoidal (divergence-free) field.Hence, any vector electric field E can be written as the superposition where E es is an electrostatic electric field that is everywhere curl-free and E ind is an induction electric field that is everywhere divergence-free.Further, the field E can be determined everywhere in space from knowledge of the values of the divergence of the field E and the curl of the field E everywhere in space [14].≠ one can obtain E es everywhere in space and from all of the regions where 0 ind E ∇ × ≠ one can obtain E ind everywhere in space.Knowing E es everywhere in space and E ind everywhere in space, one can obtain E everywhere in space.
For this paper the following equivalence of adjectives for each of the two types of electric fields is declared.First type: inductive = divergence-free = solenoidal = non-potential = transverse.Second type: electrostatic = curl-free = irrotational = potential = longitudinal.For an electrostatic field E es , the electric-field lines end; for an induction field E ind , the electric-field lines do not end.
Here we will define a purely electrostatic electric field E es as an irrotational field having its origin in regions where the divergence E ∇ ⋅ is nonzero.Likewise, we will define a purely induction electric field E ind as a solenoidal field having its origin in regions where the curl E ∇ × is nonzero.From this and from expression (3) it follows that it is seen that the electrostatic field originates from regions of nonzero charge density ρ q .From Faraday's law of induction ( ) it is seen that the induction field originates from regions of nonzero ∂B/∂t.It is important to explore expression (5) a little further.Taking the curl of Faraday's law ( ) ( ) with the displacement current neglected (negligible in the v c limit [15]) to replace B ∇ × , and taking 0 (with 0 ind E ∇ ⋅ = ) for the induction field, which indicates that the induction electric field E ind has its origin in regions where the time derivative of the current density ∂j/∂t is nonzero.
The motional electric field ( ) of the solar wind is sometimes called the induced or inductive electric field of the solar wind because it is "induced" by the motion of the plasma.In this paper more care will be taken with the "induction" nomenclature.

Observations of the Solar-Wind Motional Electric Field at Earth
In Figure 1, the solar-wind speed and two components transverse to radial (from the Sun) of solar-wind motional electric field near Earth are plotted as functions of time for 100 days in the year 2004.The GSE (Geocentric Solar Ecliptic) right-handed coordinate system is used where X points from the Earth toward the Sun, Y points transverse to the Sun-Earth line in the ecliptic plane opposite to the direction of the Earth's motion about the Sun, and Z is normal to the ecliptic plane pointing in the northward direction.In the top panel the solar-wind radial velocity v sw is plotted, in the middle panel and bottom panel the north-south component E z and the dawn-dusk component E y of the solar-wind motional electric field ( ) × are plotted.The black points are 64-s time resolution values and the red points are 24-hr running averages.The electric field is calculated from measurements of the bulk-flow vector v sw by the SWEPAM instrument [16] on the ACE spacecraft in the solar wind and from measurements of the magnetic-field vector B from the MAG instrument [17] on ACE.The value of the observed motional electric field depends on the reference frame of the observer, in this case the ACE spacecraft in approximately the rest frame of the Earth.Note the reversals in the sign of the north-south component E z (e.g.Days 111, 122, 139, 149, 165, 177, and 191): these reversals occur at sector reversals of the solar-wind magnetic field wherein the into-the-Sun versus out-of-the-Sun polarity of the solar-wind magnetic field changes.These magnetic sector reversals typically occur prior to intervals of high-speed wind.A higher-time-resolution look at the motional electric field of the solar wind appears in Figure 2. Here, plasma [18] and magnetic-field [19] measurements from the WIND spacecraft with 3-s time resolution are used to calculate ( ) × .In the top panel E y and E z are plotted for two hours of fast solar wind and in the bottom panel E y and E z are plotted for two hours of slow solar wind.Note in Figure 2 that the temporal behavior of the motional electric field is dominated by sudden changes of the values of E y and E z ; these sudden changes occur when a current sheet in the solar wind passes the spacecraft.This is the "mesoscale structure" of the solar wind, dominated by magnetic flux tubes bounded by current sheets.
In the upper two panels of Figure 3 the vector tips of the solar-wind motional electric field ( ) × are plotted in the E y -versus-E z plane.The left panel is for 24-hr averages of the electric field (from the multispacecraft OMNI2 data set [20]) and the right panel is for 64-s values of the electric field (from the ACE spacecraft).The 24-hr averages are much more dominated by a north-south orientation to E, whereas the 64-s values in the right panel are more isotropic.Basically, the deviation of ( ) × from purely north-south is owed to the deviation of the solar-wind magnetic field B from the Parker-spiral direction.In the lower two panels of Figure 3 the vector tips of the solar-wind magnetic field B are plotted in the B y -versus-B z plane, left panel being 24-hr averages of the magnetic field (from OMNI2) and the right panel being 64-s values (from ACE).The 24-hr averages of the electric and magnetic field (left panels) are dominated by a Parker-spiral orientation with B y > 0 and E z > 0 for out-of-the-Sun (away) magnetic fields and B y < 0 and E z < 0 for into-the-Sun (toward) magnetic fields.

The Parker-Spiral Motional Electric Field
Owing to the 27-day rotation of the Sun, the magnetic field in the solar-wind plasma has a global Archimidesspiral pattern.In (r, θ, ϕ) spherical coordinates the Parker-spiral magnetic field B at the radius r from the Sun and colatitude θ in the heliosphere is given by [1] ( ) where B o is a reference magnetic-field strength near the Sun (with B o > 0 for an out-of-the-Sun (away) field and B o < 0 for an into-the-Sun (toward) field); r o is the reference radius where B o is measured; ω ≈ 2.9 × 10 −6 radians/s is the angular rotation rate of the solar surface at the solar equator; and v sw is the solar-wind radial velocity.The angle θ = 0˚ at the north pole of the Sun and θ = 180˚ at the south pole.In general, the Parker-spiral magnetic field of the heliosphere is of one sector (away or toward) in the northern hemisphere of the heliosphere and is of the opposite sector (toward or away) in the southern hemisphere, with a heliospheric current sheet separating the two morphologies [21].In Figure 4 the Parker-spiral motional electric field of the solar wind will be sketched with the heliospheric current sheet lying in the equatorial plane bounding an away sector (B o > 0) in the northern hemisphere and a toward sector (B o < 0) in the southern hemisphere.
For the Parker-spiral magnetic field given by expression (7), the motional electric field of expression (2) for a radial solar-wind velocity v sw = (v sw , 0, 0) is ( ) this electric field is sketched as the red field lines in Figure 4.The motional electric field E for the Parker spiral is purely poloidal: E is northward for an away-from-the-Sun magnetic-field orientation and E is southward for a toward-the-Sun magnetic-field orientation.Note in expression ( 8) that E θ does not depend on v sw .
Examining expression (8) in spherical coordinates one finds 0 E ∇ × = everywhere for the Parker-spiral electric field.The distributed charge density in the solar-wind plasma that supports this global electrostatic electric field is given by Coulomb's law 4π q E ρ ∇ ⋅ = : taking the divergence of expression ( 8) in spherical coordinates yields ( ) a distributed charge density throughout the heliosphere.This charge density is sketched as the blue "+" and "−" symbols in Figure 4.According to expression (9) the charge density is strong near the polar regions of the heliosphere (θ ~ 0˚ and θ ~ 180˚) and weak near the equatorial regions (θ ~ 90˚).For an away-from-the-Sun magnetic-field orientation (B o > 0) the distributed charge density ρ q is negative in the northern hemisphere and positive in the southern hemisphere and for a toward orientation (B o < 0) the distributed charge density ρ q is positive in the northern hemisphere and negative in the southern hemisphere.Note in expression ( 9) that ρ q does not depend on the solar-wind velocity v sw .With 0 E ∇ × = everywhere, the global Parker-spiral electric field is purely electrostatic with its origin in regions where 0 E ∇ ⋅ ≠ .Note in Figure 4 that the electric-field lines of the global Parker-spiral motional electric field end.
The distributed current density j in the solar-wind plasma that supports the global Parker-spiral magnetic-field geometry is calculated by applying Ampere's law ( ) to the Parker-spiral magnetic field given by expression (7): taking the curl of expression (7) in spherical coordinates yields ( ) , 0, 0 This distributed current density is purely radial and it is strong in the polar regions of the heliosphere and weak in the equatorial regions, with j r negative (toward the Sun) when the magnetic field is away from the Sun (B o > 0) and j r is positive (away from the Sun) when the magnetic field is toward the Sun (B o < 0).Taking the heliospheric current sheet to lie in the equatorial plane of the heliosphere there is a radial component of the current in the sheet associated with the jump in the azimuthal component of the Parker-spiral magnetic field across the current sheet and there is an azimuthal current in the sheet associated with the reversal of the radial component of the magnetic field across the current sheet.Taking the thickness of the heliospheric current sheet to be d, Ampere's law gives ( ) ( ) for the magnitudes of the radial and azimuthal current densities in the heliospheric current sheet, note that the thickness d of the heliospheric current sheet could be a function of r, so the current density in the heliospheric current sheet could vary with r differently than the right-hand sides of expressions (11a) and (11b).Note that the poloidal motional electric field of the solar wind reverses across the current sheet owing to the reversal of the direction of B across the sheet as v sw remains radially outward.This electric-field reversal across the current sheet is sketched in Figure 4.There is a nonzero charge density in the equatorial current sheet associated with the reversal of the poloidal electric field: Coulomb's law gives this charge density ρ q as ( ) ( ) where again d is the thickness of the current sheet, and where again d could be a function of r.The charge density of the heliospheric current sheet is positive if the magnetic field is outward (from the Sun) in the northward hemisphere and inward in the southern hemisphere (as drawn in Figure 4): the charge density of the heliospheric current sheet is negative if the magnetic field is inward in the northern hemisphere and outward in the southern hemisphere.The source of this charge density is the relativistic transformation of current density into charge density [22], which will be discussed in Section 5.1.
Owing to the orbital motion of the Earth around the Sun, the radial component B r of the Parker-spiral magnetic field gives rise to an additional north-south poloidal electric field as seen by the Earth.The orbital velocity of the Earth around the Sun is v orbit ≈ 29.8 km/s in the dawnward direction, which is small compared with v sw and so the orbital motional electric field is small compared with the Parker-spiral motional electric field.The motional electric field ( ) in the reference frame of the Earth is southward when the Earth is in an away magnetic sector and northward when the Earth is in a toward magnetic sector, always opposite to the sign of the Parker-spiral motional electric field.
In the real solar wind there are several complications to this ideal-Parker-spiral picture.First, the global solar wind contains solar ejecta, which are impulsively ejected large-scale plasmas with non-Parker-spiral-oriented magnetic fields.Second, the solar-wind velocity is not uniform everywhere: rather there are regions on the rotating Sun that emit fast wind and regions that emit slower wind, disturbing the uniform-Parker-spiral picture with compressions and rarefactions.Third, there is a tilt to the solar magnetic dipole, plus magnetic multipoles that invalidate the uniform-northern-hemisphere uniform-southern-hemisphere picture of Figure 4. Fifth, the heliospheric current sheet was drawn to be planar in Figure 4 but in reality it is tilted and warped.Finally, there are large-amplitude fluctuations in the magnetic-field orientation about the Parker-spiral direction representing a mesoscale spatial structure of the solar-wind magnetic field.This mesoscale structure is considered in the following section.

The Mesoscale Electric Field Structure of the Solar Wind
In this section, the mesoscale spatial structure of the solar-wind motional electric field in the reference frame of the Earth is examined.That electric-field spatial structure is associated with the mesoscale magnetic-field spatial structure of the solar-wind plasma, which is largely flux tubes separated by current sheets [3] [4].The advection of these magnetic spatial structures past the Earth results in temporal variations of E and B at the Earth mani-fested as sudden large changes in the vectors E and B associated with each current sheet passing the Earth (cf. Figure 2).

Idealized Cases
The question of whether the mesoscale solar-wind electric field is an electrostatic field originating from charge density or an induction field originating from time-varying current density is first examined in two idealized cases drawn in the ecliptic plane in Figure 5.The first case (top panel) has the current sheets of the solar wind Figure 5. Two idealized sketches in the ecliptic plane of the solar wind plasma advecting past the Earth with current sheets of infinite extent.In the top panel the current sheets of the solar wind (black) are oriented parallel to the solar wind velocity vector and in the bottom panel the current sheets of the solar wind are oriented with their normal parallel to the solar wind velocity vector.The direction of the solar-wind velocity is indicated in purple, the direction of the magnetic field is indicated in blue, and the direction of the motional electric field (in the Earth's frame) is indicated in red.The motional electric field in the top panel is purely electrostatic (divergence origin, curl-free) and the motional electric field in the bottom panel is purely inductive (curl origin, divergence-free).
oriented along the Sun-Earth line and the second case (bottom panel) has the current sheets of the solar wind oriented normal to the Sun-Earth line.In both cases the plasma in the vicinity of the Earth is divided by two current sheets (black lines) into three plasmas, with the Earth within plasma #2.In both cases the solar-wind plasma is taken to have the uniform velocity v sw (purple vector) everywhere, with v sw parallel to the Sun-Earth line, with the Sun off the bottom of the drawing in the distance.In both cases plasma #2 is taken to have a strong uniform magnetic field B and plasmas #1 and #3 are taken to have weaker uniform magnetic fields B; the B vector (blue) is taken everywhere to be vertical out of the plane of the sketches with the strength of B changing across the current sheets.
In the top panel of Figure 5 the direction of the current density j in each current sheet is indicated by a green arrow, with the jump in the magnetic-field vector ∆B across the current sheet being related to the current density.In the reference frame of the Earth, the motional electric field in the solar wind is given by ( ) × in all three plasmas with the same value of v sw : the motional electric field ( ) ) is uniform and strong within plasma #2 and ( ) × is uniform and weak in plasmas #1 and #3.In the top panel of Figure 5 the normal component of the motional electric field changes across the current sheets, which is of the nature of 0 E ∇ ⋅ ≠ in the sheets; hence according to expression (2) there are nonzero charge densities ρ q in the current sheets of the top panel.These charge densities are indicated by the red "−" and "+" symbols in the sketch.This charge density may seem at first paradoxical to the reader because in the reference frame of the solar-wind plasma (where v = 0 everywhere and E = 0 everywhere) there is 0 E ∇ ⋅ = and no ρ q at the current sheets.The nonzero ρ q seen in the Earth's reference frame has its origin as a Lorentz transformation ( ) of the 3-vector current density j in a current density 4-vector (j, cρ q ) [22] [23].(Of course the total charge is conserved in the Lorentz transformation, but the spatial volumes containing the ions and the electrons moving at different velocities in the nonzero current are Lorentz contracted differently, resulting in a nonzero charge density out of equal amounts of positive and negative charge [14] [23].)As can be seen by examining the geometry of the electric field drawn in the top panel of Figure 5, there is no variation in E in the direction perpendicular to E, so 0 E ∇ × = .Because 0 E ∇ × = , in the top panel of Figure 5 the motional electric field of the mesoscale structure of the solar wind given by expression (2) as seen by the Earth is purely electrostatic E es and has its origin in regions of nonzero charge density where the divergence of the electric field is nonzero.
As an aside, note that if one applies the Lorentz transformation of expression (13) to the radial current density j r that supports the Parker-spiral magnetic field (expression (10)), then one gets the distributed charged density ρ q (expression (9)) of the Parker-spiral electric field.Likewise, a Lorentz transformation of the current density in the heliospheric current sheet that corresponds to the reversal of the Parker-spiral magnetic field (expression (11a)) yields the charge density of the heliospheric current sheet that corresponds to the reversal of the Parker-spiral electric field (expression (12)).Again, those charge densities arise from a difference in the Lorentz contraction of the volumes containing the differently moving ions and electrons in the currents.
The second idealized case (bottom panel of Figure 5) has the current sheets of the solar wind oriented with the current-sheet normal vectors parallel to the solar-wind velocity vector v sw .In the sketch two current sheets (black) divide the region into three plasmas, with the magnetic-field strength of plasma #2 uniform and strong and the field strength in plasmas #1 and #3 uniform and weak.The velocity v sw (purple vector) is everywhere the same and so the motional electric field ( ) vector) as seen by the Earth is uniform and strong in plasma #2 and uniform and weak in each of plasmas #1 and #3.In the bottom panel of Figure 5 the tangential component of the motional electric field E changes across the current sheets, which is of the nature of 0 E ∇ × ≠ , with a nonzero curl.At the locations of the current sheets there is a nonzero ∂j/∂t in the reference frame of the Earth owing to the motion of the current sheet as seen by the Earth.Hence, in the bottom panel of Figure 5 the motional electric field of the solar wind given by expression (2) as seen by the Earth is purely inductive E ind and has its origin at regions of ∂j/∂t (e.g.∂B/∂t) where the curl of the electric field is nonzero.

A Parker-Spiral-Aligned Flux Tube
The constructions in the two panels of Figure 5 were idealizations to produce a purely electrostatic motional electric field (top panel) and a purely induction motional electric field (bottom panel).In actuality the current sheets of the solar wind have a wide range of orientations [2] [4] [24], and the current sheets are not infinite in extent [2] [24].Usually, the magnetic field of the solar wind at Earth is approximately in the Parker-spiral orientation with the individual magnetic-flux tubes aligned approximately with the Parker-spiral direction [25].In the top panel of Figure 6 a single Parker-spiral-oriented magnetic flux tube (pink) is sketched in the ecliptic plane at Earth.The magnetic field orientation (blue vectors) is away from the Sun inside the flux tube and the motional electric field in the reference frame of the Earth (red vectors) is northward directed up out of the ecliptic plane.The magnetic-field strength is sketched strong within the tube and weak outside of the tube so the motional electric field is strong within the tube and weak outside of the tube.In the bottom panel of Figure 6 the Parker-spiral-oriented flux tube of the top panel is viewed along the axis of the tube looking back along the spiral toward the Sun with the magnetic field pointed at the viewer.The direction in the sketch of the bottom panel is the Earth's north-south direction and the flux tube is drawn to have a square cross section.The direction of the current density in the current sheets (black) surrounding the pink flux tube is indicated with green vectors.In the reference frame with the Earth at rest, the tube and its surrounding plasma advects toward the right with the solar-wind velocity v sw away from the Sun; the projection of v sw is drawn as the purple vector in the bottom panel.The motional electric field ( ) ).The top and bottom current sheets make a contribution to an electrostatic electric field in the tube and the leading and trailing current sheets make a contribution to an induction electric field in the tube.
In  ( ) The superposition E = E es + E ind of the two electric fields is calculated in the bottom panel: this superposed electric field E is the field that would be measured by a detector.In Figure 7 the tube is taken to have a uniform magnetic field B within the tube and zero magnetic field outside the tube.The current density j in each of the four current sheets bounding the tube is given by Ampere's law to be jw = (c/4π)B, where w is the thickness of a current sheet.As is the case for the flux tube in the bottom panel of Figure 6, the flux tube in Figure 7 is taken to advect in the Earth's frame from left to right with the solar-wind velocity and the magnetic field inside the tube points out of the page of the sketch.The charge density in the top and bottom current sheets is given by the relativistic transformation (taking γ = 1) ( ) 13)), so the charge density is ( ) ( ) in the leading and trailing current sheets of the tube are taken to have , with ∂B/∂t positive (up out of the page of the figure) in the leading edge and negative (into the page) in the trailing edge.For the top and bottom current sheets of the flux tube the electrostatic electric field E es everywhere in space is numerically calculated from expression (14a) by the Coulomb integral over the charge densities of the top and bottom sheets (e.g.eq.(5.15) of [6]) A sampling of electric-field lines E es (x) are numerically traced in the top panel of Figure 7.Note that these electrostatic-electric-field lines end.For the right and left current sheets of the flux tube the induction electric field E ind everywhere in space is numerically calculated from the Biot-Savart integrals applied to Faraday's law of induction (expression (14b)) with the integration over ∂B/∂t in the right and left current sheets (e.g.eq. ( A sampling of the induction electric-field lines are numerically traced and plotted in the second panel of Figure 7.Note that these induction-electric-field lines do not end.The vector sum (superposition) of the electrostatic field E es and the induction field E ind is numerically calculated everywhere in space and a sampling of the field lines of this superposed field are plotted in the third panel of Figure 7.Note that the superposition E es + E ind of the electrostatic field and the induction field is a uniform electric field within the flux tube and zero electric field outside of the flux tube (cf.expression (2)).Note also in the bottom panel of Figure 7 that the electric field lines begin and end on the top boundaries (where there is nonzero charge density) and that the tangential component of the electric field jumps across the left and right moving boundaries (where there is nonzero ∂j/∂t).
In the top panel of Figure 8 the ratio ( ) ( ) − formed from the magnitudes of the electrostatic field E es and the induction field E ind is calculated at every point inside the square flux tube of Figure 7.The ratio R is +1 if the electric field is purely electrostatic and R is −1 of the electric field is purely inductive.In Figure 8 colors approaching R = +1 are red and colors approaching R = −1 are violet.As can be seen in the top panel of Figure 8, the field has an electrostatic dominance near the top and bottom of the tube (where the charge density resides) and the field has an induction dominance near the leading and trailing edges (where the ∂B/∂t or ∂j/∂t reside).The green curve in Figure 9 is a binning of the value of R for all points inside the square flux tube of the top panel of Figure 8; the mean value and standard deviation of the R values in the interior of the square flux tube is <R> = 0.0 ± 0.14.
In the middle and bottom panels of Figure 8 the ratio ( ) ( ) + is plotted in color for a flux tube with a cross section that is tall in the north-south direction (middle) and a flux tube with a cross section that is elongated in the sunward-antisunward direction (bottom).Red colors represent regions where the value of R approaches +1 (electrostatic dominance) and violet colors represent regions where the value of R approaches −1 (induction dominance).The tall flux tube in the middle panel has a motional electric field that is dominantly induction, with the mean (and standard deviation) value <R> = −0.53± 0.19 for the interior of the tube.The flattened flux tube in the bottom panel has a motional electric field that is dominantly electrostatic, with the mean (and standard deviation) value <R> = +0.53± 0.19 for the interior of the tube.In Figure 9 the

The Mesoscale Structure of the Solar Wind
In reality the mesoscale magnetic structure of the solar-wind plasma is quite complicated, resembling a "spaghetti" of tangled tubes of various cross-sectional shapes, various cross-sectional sizes, and various orientations [3] [4].Instead of the idealized rectangular cross sections utilized in Section 5.2, the cross-sectional shapes of solar-wind flux tubes are more likely a Voronoi pattern of tubes packed against each other, each with a complicated polygonal shape (cf. Figure 1 of [4]).Portions of current sheets that are obliquely aligned with the solar-wind velocity vector will have both ( ) and 0 E ∇ × ≠ ), giving rise to electrostatic field and to induction field.Examining the panels of Figure 8, it is seen that if the cross-sectional shape of a flux tube is tall in the north-south direction, then the leading-edge and trailing-edge current sheets will dominate and the electric field within the tube will have dominantly an induction ∂j/∂t source; on the contrary if the cross-sectional shape of a flux tube is flat in the ecliptic plane then the top and bottom current sheets will dominate and the electric field within the tube will have dominantly an electrostatic ρ q source.In general, at Earth the electrostatic versus induction dominance of the time variations of the solar-wind motional electric field will vary with time as magnetic structures with varying shapes pass the Earth.However, regions of the solar wind that are compressed or rarefacted show systematic changes to the measured current-sheet orientations indicating flux tubes that are flattened in the radial direction (in compression regions) or stretched in the radial direction (in rarefaction regions) [26].Therefore, in regions where the solar-wind plasma has undergone large-scale radial compressions (e.g.behind interplanetary shocks or in corotating interaction regions) the solar-wind motional electric field of the mesoscale structure should be dominantly induction and in regions where the solar-wind plasma has undergone large-scale radial rarefactions (e.g. in the trailing edges of high-speed streams and behind fast magnetic clouds) the solar-wind motional electric field of the mesoscale structure should be dominantly electrostatic.

Summary, Refinements, and Conclusion
A Helmholtz decomposition of the motional electric field in the solar wind was made and the electrostatic versus induction nature of the field was analyzed as was the charge-density origin versus the time-varying-currentdensity origin of the field.To calculate the motional electric fields, Lorentz transformations were used; when Lorentz transforming the electric and magnetic fields, the currents and charges also had to be Lorentz transformed.Just as the transformation of the motional electric field E from B is important at non-relativistic so-lar-wind velocities, the transformation of the charge density ρ q from the current density j is important at non-relativistic velocities.
The large-scale time-averaged electric field of the solar wind associated with the outward advection of the Parker-spiral magnetic field in the solar-wind plasma is an electrostatic electric field having its source in the distributed charge density in the solar-wind plasma and the more-concentrated charge density of the heliospheric current sheet.This large-scale electric field is in the poloidal direction, northward at Earth in away magnetic sectors and southward at Earth in toward magnetic sectors.The origin of the large-scale electric field is a relativistic transformation of the large-scale magnetic field in the solar-wind plasma and the origin of the distributed charge density is the same relativistic transformation of the distributed current density in the solar-wind plasma that supports the Parker-spiral magnetic field.
The mesoscale structure of the solar-wind electric field spatially varies from being dominantly electrostatic to being dominantly inductive.At the Earth the time-varying electric field will be dominantly inductive near current sheets that have their normal vectors nearly aligned with the Sun-Earth line and the time-varying electric field will be dominantly electrostatic near current sheets that have their normal vectors nearly perpendicular to the Sun-Earth line.In general the mesoscale electric field will be a combination of electrostatic and induction electric fields, varying with time as differently shaped magnetic structures pass over the Earth.In compression regions of the solar wind time-varying induction fields will tend to dominate, and in rarefaction regions of the solar wind time-varying electrostatic fields will tend to dominate.
Note that in the actual solar wind the solar-wind velocity vector v sw also has spatial structure, on both large scales and mesoscales.On the mesoscale, there is usually a jump in the plasma velocity vector v as well as the magnetic-field vector B across the current sheets in the solar wind [27].The jump in v across the current sheet amounts to a vorticity v ω = ∇ × in the sheets; in a plasma there is a charge density ρ q associated with vorticity [28] with the charge density in a vorticity layer given by eq. ( 11) of [9] ( ) (Note that for charge density associated with vorticity there is no reference frame in which the charge density is zero.)In general the change in velocity ∆v across a solar-wind current sheet is much smaller in magnitude than the bulk flow velocity v sw of the solar-wind plasma, so the inhomogeneity of the flow-velocity field is typically a less-than-10% effect on magnitude of the motional electric fields of the solar wind.
To conclude, a first-principles understanding of the motional electric field of the solar wind associated with the advection of magnetic structure in the solar-wind plasma has been developed.This understanding is based on the Lorentz transformations of electric and magnetic fields, the Lorentz transformations of current densities and charge densities, and the application of Maxwell's equations.

Discussion: Charge Density in the Solar-Wind Plasma
As seen in Sections 4 and 5, there is a nonzero charge density ρ q in the solar-wind plasma associated with the motional electric field of the solar-wind: a distributed charge density associated with the advection of the global Parker-spiral magnetic field and more-concentrated charge densities in the current sheets of the solar wind.The physical origin of that charge density is the difference in the Lorentz contractions of volumes containing the positive ions of the solar wind and the negative electrons of the solar wind when the ions and electrons move at (slightly) different velocities owing to the presence of an electrical current.The differences in ion versus electron velocity are on the order of 1 km/s in the solar-wind current sheets.
The charge density in the solar-wind plasma associated with the motional electric field is frame dependent, being zero in the reference frame of the plasma and nonzero in other reference frames such as the frame of the Earth or the Sun.
In the directions perpendicular to B, the solar-wind plasma behaves on large scales as an ideal-MHD fluid: in the parallel-to-B direction the MHD description of the collisionless solar wind can fail (cf.Sect. 1 of [29] for MHD tests that are failed by the solar wind).Sometimes in space physics textbooks it is declared that MHD has 0 E ∇ ⋅ = (e.g.eq.(5.22) of Parks [30], eq.(4.82) of Gombosi [31], eq.(6.1.31) of Gurnett and Bhattacharjee [32], or Table 6-6 of Knipp [33]) or equivalently textbooks declare that ρ q = 0 (e.g.eq.(5.23) of Parks [30] or Sect.6.1 of Gurnett and Bhattacharjee [32]), which cannot be true relativistically (cf.expressions ( 13) and ( 17)) and which is not true in MHD computer simulations of the solar-wind flow (cf.Sect.2.3 of reference [9]).
In deriving ideal MHD, a quasineutrality condition n p ≈ n e (ion number density approximately equals the electron number density) is often assumed.The assumption of quasineutrality does not mean that Coulomb's law 4π q E ρ ∇ ⋅ = is irrelevant or inconsistent with the application of ideal MHD.For the charge density in the solar wind, one can assess just how quasineutral the solar-wind plasma is by calculating the fractional charge imbalance n q /n o , where n q = ρ q /e is the number density of excess charges and where n o ≈ n p ≈ n e is the number density of the plasma.A current sheet of the mesoscale structure of the solar wind plasma is considered.If the magnetic-field vector jumps by ∆B across the current sheet of thickness w, then the current density in the sheet is 4π j c B w = ∆ .The charge density ρ q in the current sheet from the transformation ( ) .If the magnetic-field vector across the current sheet jumps by ∆B = 5 × 10 −5 Gauss, if the width of the current sheet is w = 2500 km, and if the speed of the solar wind is v sw = 438 km/s (cf.Table 1), then the estimate for the number density of charges is n q = 4.8 × 10 −8 cm −3 .A typical number density of the solar-wind plasma is n o = 6.8 cm −3 (cf.Table 1).Thus n q /n o ~ 10 −8 for the charge density of the solar wind.
This amount of charge imbalance is small, but it is not zero.Even with n q /n o ~ 10 −8 in the solar wind, the limit ρ q = 0 cannot be taken and 0 E ∇ ⋅ ≠ : this is demonstrated with the following example.One can compare the solar-wind degree of quasineutrality to the degree of quasineutrality in a lightning bolt.The linear charge density Q/L in a conducting lightning channel is about Q/L = 10 −4 C/m [34].Dividing by e, this is a number of excess charges per unit length N/L = 6.2 × 10 6 cm −1 .A conducting lightning channel has a radius r of about r = 1 cm [35], so the net number density of charges n q in the channel is n q = 2 × 10 6 cm −3  .The number density n o of free electrons in the highly ionized lightning channel is about n o = 10 18 cm −3 [36].Hence, the fractional charge density in the channel is n q /n o ~ 2 × 10 −12 .(In actuality, the charge density n q inside the channel is lower than 2 × 10 6 cm −3 since the strong radial electrostatic electric field of the charged lightning channel pushes excess charge out of the channel where it migrates radially outward into the unionized air, into a "corona sheath" of radius of a meter or so [37]).Note that the lightning channel, with an electrostatic potential of ~10 8 V [34], is several orders of magnitude more charge neutral than the n q /n o ~ 10 −8 of the solar-wind plasma.Lightning is an electrodynamic problem in which Coulomb's law (Gauss' law, Poisson's equation, capacitance calculations) cannot be ignored (e.g.[38] [39]).
Further, the condition of quasineutrality does not mean that the electric field is negligibly small.Certainly in calculating test-particle orbits in a plasma (or MHD fluid) both E and B are needed [40].In Gaussian-cgs units the electric field E and the magnetic field B have the same units (1 gauss = 1 Stat Volt/cm = 1 gm 1/2 cm −1/2 s −1 ) so their magnitudes are directly comparable: even though the magnitude of the motional electric field is v/c times smaller than the magnitude of the magnetic field (cf.expression (2)), the electric force on a particle qE is c/v times larger than the magnetic force q(1/c)v × B on a particle.

Figure 1 .
Figure 1. 100 days of the motional electric field E from the ACE spacecraft upstream of the Earth are shown.In the top panel the solar-wind speed v sw is plotted with 64-s time resolution.In the middle panel the north-south component of E is plotted: black points are 64-s time resolution and the red points are 24-hr running averages of the black points.In the bottom panel the dawn-dusk component of E is plotted: black points are 64-s time resolution and the red points are 24-hr running averages of the black points.The GSE coordinate system is used.

Figure 2 .
Figure 2. The motional electric field of the solar wind is plotted for 2 hours of fast solar wind (top panel) and for 2 hours of slow solar wind (bottom panel).The 3-s time-resolution values are from the WIND spacecraft upstream of the Earth.E y is in red and E z is in blue.The GSE coordinate system is used.

Figure 3 .
Figure 3.In the top two panels the vector tips of solar-wind motional electric field at Earth are plotted in the (GSE) E z -versus-E z plane.The plane is viewed from the Sun looking at the Earth.In the left panel 24-hr averages of E from OMNI2 are used and in the right panel 64-s resolution values of E from ACE are used.In the bottom two panels the vector tips of the solar-wind magnetic field are plotted, 24-hr averages from OMNI2 on the left and 64-s values from ACE on the right.

Figure 4 .
Figure 4.As viewed from an observer in the ecliptic plane far from the Sun (yellow circle), Parker-spiral electric field lines are sketched in red and the charge density supporting that electric field is sketched in blue.North is toward the top.In this sketch the Parker-spiral magnetic field is away from the Sun in the northern hemisphere and toward the Sun in the southern hemisphere.The positively charged heliospheric current sheet is seen horizontally across the center of the sketch.

Figure 6 .
Figure 6.In the top panel a Parker-spiral-oriented flux tube passing the Earth is sketched in the ecliptic plane; the magnetic field in the tube has an "away from the Sun" orientation.In the bottom panel the same flux tube is viewed looking down the axis of the tube back toward the Sun.In both panels the electric fields are sketched in the reference frame of the Earth.

Figure 7 ,
the Earth's-reference-frame charge density the current sheets of a square flux tube are numerically integrated to calculate the divergence-origin electric field E es (top panel) and the curl-origin electric field E ind (middle panel) in and around the flux tube:

Figure 7 .
Figure 7.For the flux tube drawn in the bottom panel of Figure 6 (with the magnetic field pointing up out of the plane of the figure) the electric fields in the Earth's reference frame are numerically calculated.In the top panel the 4π q E ρ ∇ ⋅ = electrostatic electric field is calculated everywhere from the Coulomb's integral over the charge density in the top and bottom current sheets, in the middle panel the ( ) 1 E c B t ∇ × = − ∂ ∂ induction electric field is calculated everywhere from the Biot-Savart integral over the ∂B/∂t in the left and right current sheets, and in the bottom panel the superposition of the two fields is shown.

Figure 8 .
Figure 8.For flux tubes similar to the one in Figure 7, the local value of the ratio

Figure 9 .
Figure 9.For the three flux tubes of Figure 8, the local value of ( ) ( ) es ind es ind E E E E − +

Table 1 .
Typical values (mean values) of solar-wind parameters at Earth.