On the convection of ionospheric density features

. We investigate whether the boundaries of an ionospheric region of different density than its surroundings will drift relative to the background E × B drift and, if so, how the drift depends on the degree of density enhancement and the altitude. We ﬁnd analytic solutions for discrete circular features in a 2-D magnetised plasma. The relative drift is proportional to the density difference, which suggests that where density gradients occur they should tend to steepen on one side of a patch while 25 they are weakened on the other. This may have relevance to the morphology of polar ionospheric patches and auroral arcs, since the result is scale-invariant. There is also an altitude dependence which enters through the ion-neutral collision frequency. We discuss how the 2-D analytic result can be applied to the real ionosphere.


Introduction
We are investigating the fundamental transport properties of cold, magnetised plasmas.There appears to be a characteristic property of E × B drift which has not been previously elucidated and which we examine here.The expression for this drift field v a (subscript a for ambipolar) is which does not involve either mass or number density, or indeed any of the intrinsic properties of the plasma.The phenomenon of ambipolar drift can be understood either from the trochoidal trajectory of individual particles, or from the Lorenz transformation, which shows that for a magnetised plasma there is a "preferred" frame of reference (this v a ) in which the perpendicular electric field E ⊥ vanishes.
The E × B drift is independent of the plasma properties -for a given E. However we may ask how the electric field will be arranged in and around a plasma density feature.We show that E can and will likely be structured in such a way that E × B drift does depend inversely on the plasma mass density in most situations.The conditions under which our initial assumptions prevail in the ionosphere and magnetosphere will also be discussed.
The analysis we present entails certain simplifying assumptions, the chief of which is uniformity along magnetic field lines, also called a 2-D plasma.More precisely, it means that the location of plasma along a field line is considered unimportant, so that we need to consider only height-integrated (or field-line averaged) plasma properties.This assumption is clearly a significant one in the ionosphere, especially in the E-region.However the high parallel conductivity at all altitudes causes magnetic field lines to have a consistent electric potential, which forces the E × B drift to be mapped along the entire flux tube (propagated at the Alfvén speed).We can also show that on closed field lines our results are still applicable to an ionospheric situation.We believe that numerical modelling with realistic field-line gradients will also substantiate this analysis.
We discuss electric field rearrangements, denoted δE ⊥ , which require finite time to propagate through the 2-D domain perpendicular to B. However the low-frequency (DC) limit for perpendicular EM propagation in the ionosphere is also the Alfvén speed, v A , as shown e.g. by Baumjohann and Treumann (1996, Eq. 9.142).Typical ambipolar drifts are much slower, so one can assume that δE ⊥ propagates effectively instantaneously.
We call the average electric field far from the density features we consider the "background" electric field, E 0 .As we shall see below, the background field can only be uniform if the background density is as well.
If we were to address collisionless plasma, i.e. with no significant ion-neutral collisions, there is a train of argument we can take to show how the electric field becomes structured merely by propagating through density features.The result one obtains is the same as we obtain below, if we then look at the limit of collision frequency approaching zero; the time scale is finite and dependent on propagation, not collisions.However this case entails as much work as this already lengthy article, and will hopefully be the subject of a follow-on paper.The results below are more than adequate for any ionospheric conditions.
In this treatment, we assume that there are initally no field-aligned currents (FACs) within our domain of interest.Instead the electric field, and associated plasma drift, occur because of driving forces outside our domain.In a region of open field lines, or where the plasma is being forced locally by the neutral dynamo, the drift speed could remain independent of local density, e.g. the polar-cap potential field is "mapped" there more or less directly from the interplanetary electric field.But it still creates an electric field on closed field lines in the polar regions.The scenarios we study are an idealisation of conditions in the polar region, but as long as there is an electric field in the frame of the neutrals our results apply in some measure.
We shall examine a circular plasma density feature, either an enhancement or a depletion.For ease of reference we use the word "patch" in this paper, but with that term we don't mean only features that are two times or greater in density than the background, which is its conventional definition (e.g.Carlson, 2012).
Let the density of the patch be n times the background density.Thus a depletion is a density feature with 0 < n < 1.For ease of later notation, we find it useful to define another dimensionless quantity: Thus η serves as an alternate parameterisation of relative density such that −1 < η < 1, the bounds corresponding to an extreme depletion (η → −1 as n → 0) or an extremely dense patch (η → 1 as n → ∞).

Boundary condition on moving, sharp interface
In studying a scenario with varying density and drift speed, we must maintain conservation of particle number by species.This leads to a constraint on the speed at which the boundary drifts.
Let n be an outward normal, and let ρ s be the density of species s.If we call the drift of the sharp, step boundary v b then conservation of particle number requires that at the interface, for each species, which we can rearrange to get Also, the component of v b tangential to the boundary is arbitrary.So if we obtain a result like Eq. ( 4) where the bracketed expression on the RHS has no angular dependence, then we can drop n• from both sides of the last equation and choose the vector v b to be (nv s,int − v s,ext )/(n − 1).Now, Eq. ( 4) in this model does not include terms for ionisation and recombination, which are significant in the E region.
At the boundary itself, these terms might not be important compared to the flux of particles through the sharp boundary in an arbitrarily short time, but n is different on each side, so the ionisation/recombination balance must be different; and since the production mechanism cannot be assumed to drift with the patch, the results we obtain are valid only on the timescale in which the original density difference remains significant under ambient production.
There are further objections that could be raised to studying sharp boundaries, but we address those challenges in the Discussion.It might also be questioned why particle number but not momentum density appears to be conserved at the boundary, and we explain this in App. C.
We also wish to draw the reader's attention to the fact that we present only a steady-state picture: the generation of smallscale potential features in the ionosphere that are different from the large-scale convection pattern must trigger Alfvén waves that transport small-scale stresses back up from the ionosphere to the magnetosphere.The propagation timescales and the momentum reserves in the magnetosphere are beyond the scope of this study.We only posit that both are finite, so that there must be some conditions under which the present model has at least limited application.This caveat is essentially the same as the earlier warning that our results are only strictly appropriate for a 2-D plasma, and only applicable to the extent that a plasma has a 2-D character.
The first-order effect of currents and charge accumulations in the ionosphere is a δB which launches an Alfvén wave that establishes an FAC.But the net effect is a return to current closure, i.e. ∇ • J = 0, with a modified E. Our analysis strides over those processes and looks at the net result.Vasyliūnas (2005) offers an explanation in his Sec.3.2, especially its last paragraph.Both the large-scale FACs that we assume are present outside our domain, and the small-scale FACs that our model would implicitly trigger but which we ignore, couple the magnetosphere to the ionosphere.The net, time-dependent convection would be quite complicated.Here we are endeavouring to show only what form the ionosphere's side of the forcing terms might take.

Argument for 2-D assumption
Let us consider a 2-D plasma with collisions with a neutral gas, like the E-region but isolated from any parallel connections with a higher, collisionless region.It is uniform in the parallel direction, but we shall relate it to a realistic scenario in the Discussion.There is a background electric field in the +y direction.The Pedersen and Hall currents will be assumed to be driven by -and to close with -parallel currents far off in the ±y directions.
The FACs have their origin in plasma drifts in the magnetosphere or, in the case of open field lines, in the solar wind.FACs will create or adjust the ionospheric electric field whenever the ionospheric drift is not coherent with the magnetospheric drift: that is, the electric potential is mapped along field lines.
If E-region conductivity features begin to structure E ⊥ , then initially the F-region will create FACs that reduce the trend.
However in App.B we show that the F-region's drift kinetic energy would be absorbed by the E-region in a time on the order of a second or less, leading the F-region to adopt the E-region's E ⊥ structure also.
Then, once the whole height of the ionosphere has developed a structured E ⊥ , the magnetosphere's much greater drift kinetic energy reserve will continue to feed FACs structured in such a way as to restore the ionosphere's drift to the pattern of the magnetosphere's drift.
However, the magnetosphere is far enough from the ionosphere that it cannot provide FACs immediately, once the drift kinetic energy of the F-region has been used up by E-region conductivity.The time for an Alfvén wave to travel to the magnetic equator and back is of order 3 to 5 minutes, and there is only a finite energy available even there.Therefore it is reasonable to consider that most of the auroral ionosphere is typically in a condition where the electric field re-arrangements necessary for the E-region to be in steady state (∇ • J = 0) have been more or less mapped up throughout the F-region as well.(And the E-region would ultimately impose its electric-field structuring upwards onto the remainder of any closed field line, given sufficient time.) Moreover, the E-region's conductivity structure will generate E ⊥ structure within the time scales of the ion-neutral collision frequency ν in and of the perpendicular EM propagation, which is also at the Alfvén speed.
Therefore we feel confident in proceeding with an analysis that excludes FACs within the domain of interest, and in which the electric field in the E-region is determined only by the "background" E 0 and by E-region conductivity.

E-region steady state without local FAC
In a uniform 2-D plasma the current density is where σ P and σ H are the Pedersen and Hall conductivities.(We use the letter σ elsewhere for a surface-charge density, so to avoid confusion we write the conductivity matrix as [σ].)The components are equal to In order for the Pedersen current across the boundaries to be conserved, E1 within the slab will be lower than E0 outside.
where the sum is over all charged species; q s , n s and κ s are the charge, number density and magnetisation ratio of the species s.The latter is defined as κ s = ω s /ν s , where ω s is the cyclotron frequency of the species and ν s is its momentum transfer collision frequency with the neutrals.(We consider both ω e and κ e to be negative.) We denote quantities within the enhancement with a prime.If the patch is a factor n times the number density of the background, assuming similar composition, then 2 A "slab" feature in the E-region Let us first consider a slab geometry.Assume there is a curtain of plasma as in Fig. 1 with different density than its surroundings.
(By slab we mean something like different L shells, not horizontal strata unless we were at the equator.)We assume Cartesian geometry with B = B ẑ.
We say E-region because we are taking account of the Hall current density.The results for both the slab and the circular geometry are equally valid for F-region patches, with the considerations explained in the Discussion.
If E were uniform, the higher conductivity inside the slab would build up charge on one boundary and deplete it on the other, creating a δE oriented towards −ŷ and reducing E within the slab.In steady state, δE will be the value that restores We can see by inspection that E = E 0 + δE must be inversely proportional to the density ratio n, since the Pedersen current density has to be equal on both sides of the boundary.Provided ν in is constant, the Hall current density remains uniform as well, since the product E σ H is concomitantly fixed.The Pedersen component of ion drift is in the +y direction, faster outside the layer than within it.This means that ions are piling up on the incoming side of the layer and being peeled away on the other.The net result is that the boundary between the more and less dense regions does not move, and that the ions within the slab are transient while the electrons remain in it.

A circular, E-region patch
Suppose next that there is a circular patch of higher density, as suggested in Fig. 2. The electric field strength E must be lower inside the patch, or else charge would continually build up at the boundary.In fact in steady state there must be a net charge on the interface in order for E to be lower inside.By inspection or simple argument we can see that it will have a cylindrical dipole arrangement.
This problem has been solved already by Hysell and Drexler (2006) using complex analysis, and they have even obtained the solution for a more general elliptical problem using a conformal mapping.Up to Eq. ( 27) we provide an alternate derivation of the same result.Our method for the circular patch is somewhat less abstract, and perhaps simpler, because it does not use the double shell required for their elliptical result.
The cylindrical (or 2-D) dipole has a distinct character from the spherical dipole that is more familiar in space physics contexts.Using cylindrical polar coordinates ρ and θ, the components of a 2-D dipole aligned with the x axis are (Mallinson, 1981, Eq. 8) where µ is a constant and a = (a ρ , a θ ) represents the dipole field.Both the field lines and equipotential surfaces of a 2-D dipole are circular.So the disturbed velocity field is also a dipole rotated by 90 • as seen in Fig. 3.In App.A we provide a Cartesian expression for the dipole and details of the algebraic steps.
In Figs. 2 and 3 the dipole is aligned with the y axis, and the net charge along the circular boundary has the form σ = σ 0 sin θ; the sense of θ is shown in Fig. 2. Let E dip be a vector oriented parallel to the charge dipole and representing its maximum strength; then σ 0 = 2ε 0 E dip .The electric field in and around the patch has the form where R is the radius of the patch and D is a matrix defined in App. A.
When we take account of the Hall conductivity, we find that the dipole is no longer oriented exactly in the +y direction.
Nor can E and J inside the patch be simply parallel, but scaled-down versions of their values outside the patch as they are for the slab, because the Hall current would then accumulate on the interface as a dipole oriented towards −x.The angle which J makes with E means that the positive pole of the dipole will be shifted from 6 o'clock in that figure toward 7 o'clock.So the field inside the patch, E int , will be oriented somewhat towards 5 o'clock.See Fig. 4. If we assume a dipole orientation angle as a variable we can solve for both it and the dipole strength by requiring the current density expressions on the inner and outer sides of the circular boundary to be identical.(It is something like getting a doughnut and its hole to travel at the same speed.) Using the two mentioned parameters, this condition can be satisfied.
With E 0 = E 0 ŷ we assume the perimeter of the patch develops a surface charge density in the form of a rotated dipole: We solve for σ x and σ y under the condition that the current J into the boundary from the inside must balance the current on the outside.The outward normal is n = [cos θ, sin θ] T .

Steady state currents around a circular, E-region patch
Using the result of App.A and the definition in Eq. ( 9), δE int = −E dip .For sake of brevity we let k = E dip for the coming passage up to Eq. (20).
The current into the inner side of the boundary is Using another result of App.A, at ρ = R, which after some straightforward steps yields a current out of the boundary of Setting the cosine terms in J int and J ext equal, and setting the sine terms equal, These two linear equations in k x and k y are The LHS represents the rate of loss of net charge from the dipole σ = 2ε 0 E dip .On the RHS, the difference in [σ] between the patch and its surroundings is forcing its polarisation.Hence we can write where the matrix on the LHS is which is nearly scalar for n ≈ 1, but let us define H by writing Solving for E dip , a polarisation with yields a divergence-free current field.
The steady-state electric field in and around the patch, using Eqs.(A7) and (A8), is From this expression for E int , one can verify that it is at an angle relative to E 0 whose tangent is ησ H /σ P .This agrees with Hysell and Drexler's Eq. ( 9), which gives us confidence in our results, although we focus below on the boundary's drift, rather than that of the ions inside.

Simplifying assumptions
We shall deal here with a single ion species, and assume that electrons are fully magnetised.These assumptions are not necessary for a unique solution, but will greatly simplify the algebra.Under these assumptions, and using the ion magnetisation parameter κ i = ω i /ν in , we have Thus σ P = κ i σ H , and |κ e | is much larger than both κ i and unity.
We also assume no neutral drift, or equivalently that the electric field and all of the species' drifts are expressed in the neutrals' frame of reference.Finally, we assume a vertical magnetic field, but this is addressed in the Discussion.

Species drift
The drift of a species s (ion or electron) in the plane perpendicular to B can be described by where and are the ambipolar and Pedersen mobilities, respectively.If we let φ s be an angle defined by κ s = cot(φ s ) then we can write the mobility matrix as where R is a rotation matrix, anticlockwise if looking from positive z.We use R s below as shorthand for R( where v a is the ambipolar drift.We see that φ s is the angle which the species' drift makes w.r.t. the ambipolar drift, however we shall use the previous expression in order to get v s as a function of E.

Patch drift speed
Applying Eq. ( 4) to our dipole field, We now make use of Eq. ( A11), the fact that nT R s D = nT R T s .It can also be shown that Let us call that last matrix M s .Simplifying, The factor in square brackets on the RHS of the last equation has no θ dependence, and the tangential component of v b is arbitrary.So we may drop the n from both sides and choose where the first term on the RHS is the background drift of the species, and we have used J = [σ]E.We now use J = q i n i (v i − v e ); σ P = q i n i κ i /B(1 + κ 2 i ); and κ i = cot φ i to get (dropping the subscript 0 now since all quantities except η are background values) This equation ought to yield the same answer for either ions or electrons.We first demonstrate this for |η| 1.
For ions (s = i) we obtain The matrix is approximately a rotation, anticlockwise if we are looking parallel to B, by a small angle where in the last step we use n + 1 ≈ 2. See Fig. 5 for a sketch of this construction of v b .
For electrons (s = e) we obtain, using κ e −1  The magnitude of this drift, relative to the electron (≈ ambipolar) drift is and the orientation of the patch's drift is at right angles to J 0 .It is slower than v a for the case of density enhancements and faster for density depletions.Curiously, the drift of an enhancement's boundary has a component against the electric field.

General solution for arbitrary relative density
Still restricting ourselves to the three simplifying assumptions (one ion species, magnetised electrons, no neutral drift) we can determine the drift velocity of a patch or depletion with n much smaller or larger than unity.
Using Eq. ( 38), with either species, we can show that in the limit of large n, v b → 0. In the limit of n → 0 (a deep hole), we get v b → 2v i .And of course for n = 1 we have v b = v a .Looking at the trend for n ∼ 1 suggested by Fig. 5, and using the intuition obtained with a bit of numerical experimentation, one can appreciate that the range of values which v b can take on as a function of n, for a fixed value of κ i , describes a circular arc passing through the origin, v a and 2v i .
Since the origin and v a both lie on the circular arc, the arc's centre must lie on the line v x = 1 2 v a .Furthermore, the trend around n ∼ 1 being perpendicular to J shows that the centre lies along the line of J extended from v a .Thus the centre is at So we can separate the centre from the rest of Eq. ( 38) and, using s = e and repeating the same O(κ −1 e ) approximations used to get Eq.( 41), write Let β be an angle defined by tan β = η/κ i (this is the angle between E int and E 0 ).Then Fig. 6 illustrates the portion of the arc along which the boundary drift v b can occur for various density ratios n, for one particular value of κ i .
If we look at the electron drift velocity inside the patch using Eqs.( 26) and (32), we find that the electrons drift along with the boundary.So the patch (or depletion) keeps its original electrons (so to speak) whereas the ions are transiently within it, as in the slab geometry.
In Fig. 7 we show a series of possibilities for different ion magnetisation ratios from 0.1 to 10.Each coloured arc shows, for a fixed value of κ i , how patches or depletions (varying n, or η) should drift.An enhancement always drifts slower than ambipolar, with a component against E, whereas a depletion always has a component parallel to E. Usually a depletion drifts faster than ambipolar, but with highly demagnetised ions an enhancement's drift can be slower.
In the frame of reference of the background drift, the current vector J separates the two cases, with enhancements drifting towards its right side (boreal pole), and depletions towards its left.

Discussion
We have derived the drift behaviour of a circular density enhancement (or depletion) in a 2-D magnetised plasma.Our study grew from a 2-D auroral modeling effort in the meridional plane (de Boer et al., 2010) which led to thinking about appropriate upper boundary conditions (∼ 1000 km) for the electric potential, and the time scales over which various effects should assert themselves.The eventual conclusion was that J = 0 was appropriate over the time scales of interest, and this inquiry led to the arguments presented above in justifying the 2-D analysis in a surface perpendicular to B.
The idea that charge accumulation from a non-uniform conductivity [σ] would set up an electric field disturbance in a way that would force the system towards steady-state (∇ • J = 0) is simple, although the algebra turned out to be more complicated than expected.Perhaps a shorter derivation is possible.But we found a unique result that we are confident in.
Still, one might challenge the relevance of our result in the limit as κ i → ∞, since the time required to approach an equilibrium condition also tends to infinity.However by following other lines of argument not included in this paper, and considering the effective, perpendicular permittivity of a magnetised plasma with κ i → ∞, it can be shown that the steady-state towards which the E-region is driving the convection pattern is the same as the E ⊥ structure obtained from considering a circular patch (or depletion) with E initially zero, and increasing E 0 up to some steady value.
One should also reasonably question how these 2-D results -specifically Eq. ( 44) -can be extended to the real ionosphere, where κ i and φ i vary continuously with height.Each altitude layer will be forcing a different patch drift speed, while parallel conductivity is trying to enforce a coherent flux-tube drift.A conjecture (which awaits analytic justification) is that the effective ion magnetisation of a flux tube can be obtained from the ratio of height-integrated conductivities: and that values of φ i and β determined from this κ eff.can be used in Eq. ( 44) to find the theoretical drift direction and speed of a patch (or a depletion).Of course such a 3-D drift structure would entail dipolar structures of FACs closed within the ionosphere.There would be a special height at which κ eff.= κ i , and the dipolar FAC would be oppositely oriented above and below this height.

Addressing some idealisations
If we are considering F-region structure, our results are still applicable.We presented the analysis in the context of the E-region because we are treating the Hall current, whereas the F-region alone, with negligible Hall current, would constitute a narrower problem.Also any E-region structure is weighted much more highly than the F-region's in determining electric field structure, due to its stronger contribution to Σ P .
We assumed a vertical B field.In the polar region this is not a large approximation, but the E × B drift acts perpendicular to B so it can have a vertical component.In the polar cap, between the dayside open-close boundary (OCB) and the line across 06-18 MLT, the convection adds an upward component to the plasma's vertical momentum balance, while on the midnight side of the cap it is driving the plasma downwards.This adds a layer of complexity, but also ensures that even initially purely F-region patches should create some conductivity structure by the time the plasma reaches the nightside OCB.
It may be possible to extend our analysis to elliptical patches.However it has some complexities that might only yield to complex analysis, as Hysell and Drexler (2006) have accomplished.At least this study reinforces their result for circular patches.ionosphere developing very short gradient scale lengths and generating FAC over equally small scales, smaller than any initial density structure, without the requirement for precipitation to initiate the fine-scale structure.
This effect could explain why in auroral phenomena we see structures cascading to smaller and smaller gradient scale lengths.The initial density gradients generated in the auroral oval by precipitation are stronger to begin with than at other latitudes.Also the strong electric fields found there reduce the time scale for features to cascade to smaller spatial scales.This combination of conditions yields an increased chance that this cascade time scale might prevail over the erasure of structure on the scale of the ion chemical lifetime.

Search for observational support
The predicted structuring of the electric field around density features, and the relative drift of those features, is independent of scale.Hysell and Drexler were motivated by Farley-Buneman waves.We began with the goal of understanding small-scale auroral structure, however an observational test of our results for individual features would appear to be very difficult.The most promising avenue for testing our hypotheses will be in the observation of large-scale patches in the polar cap and auroral oval.Both Carlson (2012) and Zhang et al. (2013) have provided analyses of events using, respectively, EISCAT for Ne data and GPS receiver arrays for TEC.The latter dataset is overlaid with SuperDARN convection patterns.These papers show incontrovertibly that patches of ionisation do cross the polar cap, and can return in the sunward return flow.However within the spatial and temporal resolution of their available data, it is challenging to see either confirmation or refutation of the electric-field structuring that we posit.
There are however two ways that the effect we argue could bias patch motion statistically.Fig. 7 shows that a patch in the boreal polar cap might be expected to drift somewhat to the left (towards post-midnight) relative to the average convection, and that it should travel more slowly than the mean drift.But if it is a predominantly F-region feature, then the leftward deviation may be negligible.It is curious that in their Fig. 6 the highest speed shown is inferred from the SuperDARN convection map, and is clearly higher than the patch's speed measured by radar backscatter.Also Gillies et al. (2009) have examined the factor of about 0.75 by which SuperDARN Doppler velocities are lower than those obtained for convection from the DMSP satellites.Their work shows that about a third of the discrepancy is accounted for by taking the index of refraction into account.But it is intriguing to speculate that the remaining factor may arise because the DMSP data yield mean convection, while the Doppler speeds may be biased towards plasma with higher density and therefore stronger radar return.Perhaps the remaining systematic difference is due to an inverse correlation between plasma density and local electric field strength.
In some convection maps, it seems like it might even be possible to see the effect we posit.For example in Fig. 2

Modeling
We should address how our work is relevant to ionospheric modeling.Our result shows that the electric field cannot be simply prescribed for some region, but that it will have structure implicitly determined by the plasma density structure.Most models, including our own cited work, assume an E or potential field that is prescribed in some way.For example Schunk and Sojka (1987) used an E field that remained fixed despite the introduction of very strong density features.
A numerical model intended to address E field structuring might begin with an initial electric potential map, but the actual Pedersen and Hall currents will generate FACs wherever they converge or diverge in the ionosphere.These FACs cannot be driven immediately or indefinitely by magnetospheric processes.Charge accumulations will then force a structuring of the ionospheric and magnetospheric potential towards a situation where FACs are no longer required to maintain current closure, i.e. exactly the sort of structure we have identified.Where the plasma density is higher, the field will be lower, and vice versa.
In the limit of closed field lines, this will amount to solving the Laplace equation in 2-D.

Conclusions
The following characteristics of E × B drift in 2-D, magnetised plasma, shown by Hysell and Drexler (2006), have been confirmed through an alternate analysis: 1.While plasma on an open flux tube may have a uniform electric field more-or-less "imposed" on it regardless of density structure, plasma on closed flux tubes will experience a structuring of the steady-state electric field that depends on density features -weighted towards dependence on E-region density.
2. For a circular density feature, the assumption of a dipolar net charge with appropriate magnitude and orientation can yield a divergence-free current field.
3. A density feature does not "own" a particular parcel of ions -the ions both inside and out can convect through the boundary -nevertheless the boundary of a circular density feature retains a circular shape, and the electrons convect with the density feature.
We have also shown that: 4. The boundary of a circular feature should convect with a velocity given by Eq. ( 44) and shown in Fig. 7 -always slower than ambipolar for an enhancement and usually faster for a depletion; and with a component against or with, respectively, the background electric field.
5. An obvious extension of the result for a sharp, circular feature is that features with density gradients will see gradients on one side steepened and gradients on the other side weakened.
6.The E-region can therefore generate smaller-scale structure than its initial structure, at any length scale and even without instability present, and the time scale for this to occur is inversely related to the electric field strength.
As well, we have provided some arguments as to why and how these 2-D results are still applicable to the real ionosphere with its altitude-dependence of plasma properties.Moreover, we wish to show in a future paper that the structuring described at point 1., in a non-conducting plasma, will depend on plasma mass density features.

Figure 1 .
Figure1.A slab of denser plasma extending perpendicular to the background electric field E0 and extending parallel to B (e.g. between two L-shells).In order for the Pedersen current across the boundaries to be conserved, E1 within the slab will be lower than E0 outside.

Figure 2 .
Figure 2. A circular patch of denser plasma will acquire a polarisation that reduces the electric field strength inside compared to the background field strength.This sets up a cylindrical (2-D) dipole.The + and − show net surface charge.The sense of the angle θ used in the analysis is shown (the conventional sense, but appearing clockwise because ẑ is into the page).This figure ignores Hall current, but it is included in our analysis, and it is shown in Fig. 4. The sense of plasma drift around such a feature is sketched in Fig. 3.

Figure 3 .
Figure 3.The polarisation of the patch in Fig. 2 establishes a dipole disturbance in the E × B drift.Note that plasma flows through the patch: while the density enhancement should retain its circular shape, the constituent ions are transient.

Figure 4 .
Figure 4.The Hall current, not yet shown in Fig. 2, adds a rotation to the polarisation of the patch, however the exact solution for a sharp circular boundary is still a cylindrical dipole.

Fig. 5
Fig. 5 also sketches this construction of v b , and one can see it is coïncident with the value of v b obtained using the ions.

Figure 5 .
Figure 5.The dashed semicircle shows the locus of vi for various values of κi.The blue arrows show the construction of vb, the drift of the density enhancement, relative to either the ion or electron background drift, for a modest enhancement (n 1).

Figure 6 .
Figure 6.The construction of the circular arc (blue dotted curve) along which the drift of the boundary of a circular patch (vb) can occur for a given ion magnetisation κi.The blue text shows the variation of the drift along the arc for various values of the relative density n, which determines the angle β.In this example κi ∼ 1.5.

Figure 7 .
Figure 7.The variation w.r.t.ion magnetisation κi of the arc along which the drift of a circular patch can lie.All of the arcs coïncide at the origin (the limit of high n) and at the ambipolar drift (for the case of n = 1).An enhancement always drifts slower than ambipolar, with a component against E, whereas a depletion always has a component parallel to E.

Figure 8 .
Figure 8. Equipotential contours that satisfy ∇ • J = 0 for two concentric patches of density 2× and 5× the background density, and a uniform κ = 1.Although the E field and the boundary shapes will become quite complicated as the inner patch approaches the boundary of the outer one, we can appreciate that the step boundaries will approach each other around their 10 o'clock position, while growing farther apart around 4 o'clock.
Moen et al. (2007) did a statistical study of patch exit times and found a high degree of symmetry around midnight.There was a small bias toward pre-midnight exits, the distribution being centred on 23:25 MLT, which might only reflect the small bias of IMF By in their sample.This offset is the opposite of what we are positing, albeit in a region where convection is expected to be the least structured by conductivity.Moen et al. (2015) examined the statistical variation of patch exit times with IMF By and Bz more closely, and found that the pre-midnight shift for By positive is mirrored by a post-midnight shift for By negative.So this test is inconclusive.Next we look at drift speed.Oksavik et al. (2010) studied two particular polar cap patch events in 2001, of which we examine no. 2. Their Fig.3.e) is interesting because it shows that this patch's velocity was distinctly lower than the plasma ahead and behind it.The SuperDARN flow plotted in their Fig.2.d) shows a speed of approximately 470 m/s on the patch's right (east) side and 620 m/s on its left.(The patch undergoes a clockwise rotation of 90 • .)The SuperDARN maps for that day show electric field strengths of about 40 and 30 mV/m at 08 and 09 UT, respectively, or roughly 700 m/s at 0830.Yet the patch progressed through the east and west beams of the radar at 226 and 566 m/s, respectively.So there is some evidence that patches convect more slowly than the average convective flow around them.Hosokawa et al. (2010) looked at an event where two patches were pulled away from each other by a shear in the convection.
in Zhang et al., looking around the terminator in the polar cap in panels C, D & E, and in the return flow in panel G, it appears that the contours of potential are slightly spread out (weaker E) around the stronger TEC structures.