Lunar surface potential and electric field

The Moon has no significant atmosphere, thus its surface is exposed to solar ultraviolet radiation and the solar wind. Photoemission and collection of the solar wind electrons and ions may result in lunar surface charging. On the dayside, the surface potential is mainly determined by photoelectrons, modulated by the solar wind;while the nightside surface potential is a function of the plasma distribution in the lunar wake. Taking the plasma observations in the lunar environment as inputs, the global potential distribution is calculated according to the plasma sheath theory, assuming Maxwellian distributions for the surface emitted photoelectrons and the solar wind electrons. Results show that the lunar surface potential and sheath scale length change versus the solar zenith angle, which implies that the electric field has a horizontal component in addition to the vertical one. By differentiating the potential vertically and horizontally, we obtain the global electric field. It is found that the vertical electric field component is strongest at the subsolar point,which has a magnitude of 1 V m-1. The horizontal component is much weaker, and mainly appears near the terminator and on the nightside, with a magnitude of several mV m-1. The horizontal electric field component on the nightside is rotationally symmetric around the wake axis and is strongly determined by the plasma parameters in the lunar wake.     

The three-dimensional shape of the magnetopause

Nearly 1000 magnetopause crossings from HEOS-2, HEOS-1, OGO-5 and 5 IMP space-craft covering most of the northern and part of the southern dayside and near-Earth tail magnetopause (X >?15 R_E) have been used to perform a detailed study of the three-dimensional shape and location of the magnetopause. The long-term influence of the solar wind conditions on the average magnetopause geometry has been reduced by normalising the radial distances of the observed magnetopause crossings to an average dynamical solar wind pressure. Best-fit ellipsoids have been obtained to represent the average magnetopause surface in geocentric solar ecliptic (GSE) and (as a function of tilt angle) in solar magnetic (SM) coordinates. Average geocentric distances to the magnetopause for the 1972–1973 solar wind conditions (density 9.4 cm^?3, velocity 450 km s^?1) are 8.8 R_E in the sunward direction, 14.7 R_E in the dusk direction, 13.4 R_E in the dawn direction and 13.7 R_E in the direction normal to the ecliptic plane. The magnetopause surface is tilted by 6.6° ± 2° in a direction consistent with that expected from the aberration effect of the radial solar wind. Our data suggest that the solar wind plasma density and the interplanetary magnetic field (IMF) orientation affect the distance to the polar magnetopause, larger distances corresponding to higher plasma density and southward fields. Our best-fit magnetopause surface shows larger geocentric distances than predicted by the model of Choe et al. [Planet Space Sci. 21, 485 (1973).] normalised to the same solar wind pressure.     

Signatures of Slow Solar Wind Streams from Active Regions in the Inner Corona

The identification of solar-wind sources is an important question in solar physics. The existing solar-wind models (e.g., the Wang–Sheeley–Arge model) provide the approximate locations of the solar wind sources based on magnetic field extrapolations. It has been suggested recently that plasma outflows observed at the edges of active regions may be a source of the slow solar wind. To explore this we analyze an isolated active region (AR) adjacent to small coronal hole (CH) in July/August 2009. On 1 August, Hinode/EUV Imaging Spectrometer observations showed two compact outflow regions in the corona. Coronal rays were observed above the active-region coronal hole (ARCH) region on the eastern limb on 31 July by STEREO-A/EUVI and at the western limb on 7 August by CORONAS-Photon/TESIS telescopes. In both cases the coronal rays were co-aligned with open magnetic-field lines given by the potential field source surface model, which expanded into the streamer. The solar-wind parameters measured by STEREO-B, ACE, Wind, and STEREO-A confirmed the identification of the ARCH as a source region of the slow solar wind. The results of the study support the suggestion that coronal rays can represent signatures of outflows from ARs propagating in the inner corona along open field lines into the heliosphere.     

Thin layered field aligned current in the solar wind

The distribution of thin layered field aligned currents (current layers) in the leading edge of the fast stream at 1 a.u. is studied by comparing solar wind plasma parameters (bulk velocity, density proton temperature and magnetic field intensity) and occurrence frequency of the current layer. Each leading edge studied is either a stream interface or an interplanetary shock. It is found that the occurrence frequency is related best to the proton temperature variations; the occurrence frequency increases with the rise of proton temperature. Possible mechanisms which cause high occurrence frequency of current layers in the leading edge are discussed.     

Role of “Core” and “Halo” solar electrons in ionization of the interstellar medium

The probability of the interstellar wind atoms (H and He) to survive ionization by solar wind electrons is presented. For the first time a dual temperature electron distribution is used to model the effects of “core” (10 eV) and “halo” (60 eV) solar electrons on the probabilities. Survival probability distributions as a function of helicocentric distance were calculated for variations in the electron temperature, solar radiation force, and the interstellar wind flow velocity. These probabilities are important in determining the radial density distributions of the interstellar atoms. It has been found that the interstellar wind has a distinctively higher probability of surviving “halo” rather than “core” electron ionization only at heliocentric distances, ρ, smaller than about 0.5 a.u. For distances larger than 0.5 a.u., the probabilities of surviving “halo” electrons are close to the probabilities of surviving “core” electrons. Also, the probabilities for both “core” and “halo” electrons are relatively insensitive to changes in u_∞ (interstellar wind velocity at infinity), μ (the solar ratio of radiation to gravitational force) and α (a model parameter for solar electron temperature) for ρ > 0.5. For distances smaller than that, the sensitivity increases significantly.     

GEOS-2 measurements of cold ions in the magnetosheath

The Suprathermal Plasma Analysers on GEOS-2 are able to make differential energy measurements of plasma particles down to sub-eV energies because the entire sensor package can be biased relative to the spacecraft. When the package is biased negatively with respect to space potential, low energy positive ions are sucked in and are more easily detected against the background. Large fluxes of ions with temperatures of the order of 1 eV or less were consistently detected at space potential when the spacecraft was in the magnetosheath though not when it was in the nearby magnetosphere. This apparent geophysical correlation, suggesting that the ions were part of the magnetosheath ion population, was contradicted by the fact that the ions showed no signs of the large drift velocity associated with the electric field in the magnetosheath. We conclude, after further investigation, that the observed ions were probably sputtered as neutrals from the spacecraft surface by the impact of solar wind ions and subsequently ionized by sunlight or electron impact. The effect of sputtering by solar wind ions has not been previously observed, although it could have consequences for the long-term stability of spacecraft surfaces.     

Sources of sodium in the lunar exosphere: Modeling using ground-based observations of sodium emission and spacecraft data of the plasma

Observations of the equatorial lunar sodium emission are examined to quantify the effect of precipitating ions on source rates for the Moon’s exospheric volatile species. Using a model of exospheric sodium transport under lunar gravity forces, the measured emission intensity is normalized to a constant lunar phase angle to minimize the effect of different viewing geometries. Daily averages of the solar Lyman α flux and ion flux are used as the input variables for photon-stimulated desorption (PSD) and ion sputtering, respectively, while impact vaporization due to the micrometeoritic influx is assumed constant. Additionally, a proxy term proportional to both the Lyman α and to the ion flux is introduced to assess the importance of ion-enhanced diffusion and/or chemical sputtering. The combination of particle transport and constrained regression models demonstrates that, assuming sputtering yields that are typical of protons incident on lunar soils, the primary effect of ion impact on the surface of the Moon is not direct sputtering but rather an enhancement of the PSD efficiency. It is inferred that the ion-induced effects must double the PSD efficiency for flux typical of the solar wind at 1 AU. The enhancement in relative efficiency of PSD due to the bombardment of the lunar surface by the plasma sheet ions during passages through the Earth’s magnetotail is shown to be approximately two times higher than when it is due to solar wind ions. This leads to the conclusion that the priming of the surface is more efficiently carried out by the energetic plasma sheet ions.     

Tracking Back the Solar Wind to Its Photospheric Footpoints from Wind Observations – A Statistical Study

It is of great importance to track the solar wind back to its photospheric source region and identify the related current sheets; this will provide key information for investigating the origin and predictions of the solar wind. We report a statistical study relating the photospheric footpoint motion and in-situ observation of current sheets in the solar wind. We used the potential force-free source–surface (PFSS) model and the daily synoptic charts to trace the solar wind back from 1 AU, as observed by the Wind spacecraft, to the solar surface. As the footpoints move along the solar surface we obtain a time series of the jump times between different points. These jumps can be within a cell and between adjacent cells. We obtained the distribution of the jump times and the distribution for a subset of the jump times in which only jumps between adjacent cells were counted. For both cases, the distributions clearly show two populations. These distributions are compared with the distribution of in-situ current sheets reported in an earlier work of Miao, Peng, and Li (Ann. Geophys. 29, 237, 2011). Its implications on the origin of the current sheets are discussed.     

Three-dimensional computer modeling of dynamic reconnection in the magnetotail

Magnetospheric substorm in the magnetotail region is studied numerically by means of a three-dimensional MHD code. The analytic solution for the quiet magnetotail is emloyed as an initial configuration. The localized solar wind is modeled to enter the simulation domain through the boundaries located in the magnetotail lobe region. As a result of the interaction between the solar wind and the magnetosphere, the magnetic field lines are stretched, and the plasma sheet becomes thinner and thinner. When the current-driven resistivity is generated, magnetic reconnection is triggered by this resistivity. The resulting plasma jetting is found to be super-magnetosonic. Although the plasmoid formation and its tailward motion is not quite clear as in the two-dimensional simulation, which is mainly because of the numerical model chosen for the present simulation, the rarification of the plasmas near thex-point is observed. Field-aligned currents are observed in the late expansive stage of the magnetospheric substorm. These field-aligned currents flow from the tail toward the ionosphere on the dawn side and from the ionosphere toward the tail on the dusk side, namely in the same sence of the region 1 current. As the field-aligned currents develop, it is found that the cross-tail current in the Earthside midnight section of the magneticx-point is reduced.     

Lofted charged dust distribution above the Moon surface

We developed kinetic theory for the charging processes of small dust grains near the lunar surface due to interaction with the anisotropic solar wind plasma. Once charged, these dust grains, which are exposed to the electric field in the sheath region near the lunar surface, could loft and distribute around such heights off the surface where they reach equilibrium with the local gravitational force. Analytical solutions were derived for the charging time, grain floating potential, and grain charge, characterizing the charging processes of small dust grains in a two-component and in a multi-component solar wind plasma, and further highlighting the unique features presented by the high streaming plasma velocity. We have also formulated a novel kinetic theory of sheath formation around an absorbing planar surface immersed in the anisotropic solar wind plasma in the case of a negligible photoelectric effect and presented solutions for the sheath structure. In this study we combined the results from these analyses and provided estimates for the size distribution function of dust that is expected to be lofted in regions dominated by the solar wind plasma, such as near the terminator and in nearby shadowed craters. Corresponding to the two dominant streaming velocity peaks of 300 and 800 km/s, mean dust diameters of 500 and 350 nm, respectively, are expected to be found at equilibrium at heights of relevance to exploration operations, e.g., around 1.5 m height off the lunar surface. In shadowed craters near the terminator region, where isotropic plasma should be dominating, we estimate mean lofted dust diameter of 800 nm around the same 1.5 m height off the lunar surface. The generally applicable solutions could be used to readily calculate the expected lofted size distribution near the lunar surface as a function of plasma parameters, dust grain composition, and other parameters of interest.     

The solar wind interaction with Venus through the eyes of the Pioneer Venus Orbiter

Venus has no internal magnetic dynamo and thus its ionosphere and hot oxygen exosphere dominate the interaction with the solar wind. The solar wind at 0.72 AU has a dynamic pressure that ranges from 4.5 nPa (at solar max) to 6.6 nPa (at solar min), and its flow past the planet produces a shock of typical magnetosonic Mach number 5 at the subsolar point. At solar maximum the pressure in the ionospheric plasma is sufficient to hold off the solar wind at an altitude of 400 km above the surface at the subsolar point, and 1000 km above the terminators. The deflection of the solar wind occurs through the formation of a magnetic barrier on the inner edge of the magnetosheath, or shocked solar wind. Under typical solar wind conditions the time scale for diffusion of the magnetic field into the ionosphere is so long that the ionosphere remains field free and the barrier deflects almost all the incoming solar wind. Any neutral atoms of the hot oxygen exosphere that reach the altitude of the magnetosheath are accelerated by the electric field of the flowing magnetized plasma and swept along cycloidal paths in the antisolar direction. This pickup process, while important for the loss of the Venus atmosphere, plays a minor role in the deceleration and deflection of the solar wind. Like at magnetized planets, the Venus shock and magnetosheath generate hot electrons and ions that flow back along magnetic field lines into the solar wind to form a foreshock. A magnetic tail is created by the magnetic flux that is slowed in the interaction and becomes mass-loaded with thermal ions.The structure of the ionosphere is very much dependent on solar activity and the dynamic pressure of the solar wind. At solar maximum under typical solar wind conditions, the ionosphere is unmagnetized except for the presence of thin magnetic flux ropes. The ionospheric plasma flows freely to the nightside forming a well-developed night ionosphere. When the solar wind pressure dominates over the ionospheric pressure the ionosphere becomes completely magnetized, the flow to the nightside diminishes, and the night ionosphere weakens. Even at solar maximum the night ionosphere has a very irregular density structure. The electromagnetic environment of Venus has not been well surveyed. At ELF and VLF frequencies there is noise generated in the foreshock and shock. At low altitude in the night ionosphere noise, presumably generated by lightning, can be detected. This paper reviews the plasma environment at Venus and the physics of the solar wind interaction on the threshold of a new series of Venus exploration missions.     

Magnetic states of the ionosphere of Venus observed by Venus Express

Strong ultraviolet radiation from the Sun ionizes the upper atmosphere of Venus, creating a dense ionosphere on the dayside of the planet. In contrast to Earth, the ionosphere of Venus is not protected against the solar wind by a magnetic field. However, the interaction between charged ionospheric particles and the solar wind dynamic and magnetic pressure creates a pseudo-magnetosphere which deflects the solar wind flow around the planet (Schunk and Nagy, 1980). The combination of changing solar radiation and solar wind intensities leads to a highly variable structure and plasma composition of the ionosphere. The instrumentation of the Venus Express spacecraft allows to measure the magnetic field (MAG experiment) as well as the electron energy spectrum and the ion composition (ASPERA-4 experiment) of the upper ionosphere and ionopause. In contrast to the earlier Pioneer Venus Orbiter (PVO) measurements which were conducted during solar maximum, the solar activity was very low in the period 2006-2009. A comparison with PVO allows for an investigation of ionospheric properties under different solar wind and EUV radiation conditions. Observations of MAG and ASPERA have been analyzed to determine the positions of the photoelectron boundary (PEB) and the “magnetopause” and their dependence on the solar zenith angle (SZA). The PEB was determined using the ELS observations of ionospheric photoelectrons, which can be identified by their specific energy range. It is of particular interest to explore the different magnetic states of the ionosphere, since these influence the local plasma conductivity, currents and probably the escape of electrons and ions. The penetration of magnetic fields into the ionosphere depends on the external conditions as well as on the ionospheric properties. By analyzing a large number of orbits, using a combination of two different methods, we define criteria to distinguish between the so-called magnetized and unmagnetized ionospheric states. Furthermore, we confirm that the average magnetic field inside the ionosphere shows a linear dependence on the magnetic field in the region directly above the PEB.     

Solar wind energy transfer regions inside the dayside magnetopause—II. Evidence for an MHD generator process

In this paper a quantitative analysis of magnetosheath injection regions observed by PROGNOZ-7 in the dayside high latitude boundary layer is performed. Particular emphasis is laid on describing the consequences of the observed excess transverse momentum of solar wind ions (H⁺ and He²⁺) as compared to the magnetospheric ions (e.g. He⁺ and O⁺) in the magnetosheath injection regions, hereafter referred to as energy transfer regions.An important result of this study is that the observed excess drift velocity of the solar wind ions as compared to the magnetospheric ions can be interpreted as a negative inertia current being present in the boundary layer. This means that the inertia current goes against the local electric field and that particle kinetic energy is converted into electric energy there. The dayside high-latitude boundary layer therefore constitutes a voltage generator (at least with respect to the injected magnetosheath plasma).The MHD-theory predicts a strong coupling of the energy transfer process in the boundary layer and the ionosphere, both regions being connected by field aligned currents. The rate of decay of the inertia current in the injected plasma element is in the range of a few minutes, a value which is directly proportional to the ionospheric resistance. By taking into account both the Hall and the Pedersen conductivities in the ionosphere, the theory also predicts a strong coupling between ionospheric East/West and North/South currents. A considerable part of the inertia current may actually flow in the tangential (East/West) direction due to this coupling. Thus, a consequence of the boundary layer energy transfer process is that it may generate currents, powering other magnetospheric plasma processes, down to ionospheric heights.     

Alfven Wave Reflection model of field-aligned currents at Mercury

An Alfven Wave Reflection (AWR) model is proposed that provides closure for strong field-aligned currents (FACs) driven by the magnetopause reconnection in the magnetospheres of planets having no significant ionospheric and surface electrical conductance. The model is based on properties of the Alfven waves, generated at high altitudes and reflected from the low-conductivity surface of the planet. When magnetospheric convection is very slow, the incident and reflected Alfven waves propagate along approximately the same path. In this case, the net field-aligned currents will be small. However, as the convection speed increases, the reflected wave is displaced relatively to the incident wave so that the incident and reflected waves no longer compensate each other. In this case, the net field-aligned current may be large despite the lack of significant ionospheric and surface conductivity. Our estimate shows that for typical solar wind conditions at Mercury, the magnitude of Region 1-type FACs in Mercury’s magnetosphere may reach hundreds of kilo-Amperes. This AWR model of field-aligned currents may provide a solution to the long-standing problem of the closure of FACs in the Mercury’s magnetosphere.     

Coronal sources of high-speed plasma streams in the solar wind during the declining phase of solar cycle 20

A correlative study is made between inferred solar sources of high-speed solar wind streams and extended white-light coronal features. The solar wind data used in the study consists of 110 co-rotating high-speed plasma streams observed from spacecraft at 1 AU in the period February 1971-December 1974; the coronal data consists of 144 equatorward extensions of polar coronal holes and 15 equatorial coronal holes, derived fromK-coronometer maps of the white-light corona during the same period. Of 110 observed solar wind streams 88 could directly be associated with an equatorward extension of a polar-cap coronal hole and 14 could be associated with a low-latitude equatorial coronal hole. In 8 cases no visible coronal feature was identified. Of 144 identified polar-cap extensions 102 were associated with a high-speed stream observed at 1 AU; 19 coronal features were related in time to data gaps in the solar wind measurements, while 38 features did not give rise to solar wind streams observed at Earth orbit. The probability of an association depended on the heliographic co-latitude of a polar hole extension, being 50% for a polar lobe extending down to 45° co-latitude and 100% for a polar coronal hole extending to 80° co-latitude or more.Paper presented at the 11th European Regional Astronomical Meeting of the IAU on New Windows to the Univese, held 3–8 July, 1989, Tenerife, Canary Islands, Spain.     

Geoactive zones in the solar wind

Three parameters of the solar wind, proton number density n, Z-component of frozen-in magnetic field, in solar ecliptic coordinates and magnetic field variability ΔB, may be called geoactive parameters since each of them is responsible for a certain phase or stage of a geomagnetic storm.An undisturbed solar corpuscular stream differs from the quiet solar wind mainly in higher bulk velocity v; other parameters, in particular, n, Z and ΔB, are not enhanced in the stream. However, the examination of a number of geomagnetic storms shows that v is not a geoactive parameter. Hence the corpuscular stream itself is not more geoactive than the quiet solar wind.The retarding of corpuscular stream by the quiet solar wind results in various plasma deformations (compression, torsion, shear). This, in turn, leads to the creation, in the stream and ambient quiet solar wind, of geoactive zones. Each zone is characterized by the enhancement of some geoactive parameter. The entry of the Earth into a geoactive zone causes a corresponding phase or stage of a geomagnetic storm.The concept of geoactive zones is applied to the analysis of the geomagnetic storm of 8–10 July 1966.     

Distributions of tangential discontinuity in the corotating solar wind structure

Distributions of the tangential discontinuity (TD) in the solar wind sector structure are investigated on the basis of the magnetic field data and the ion plasma parameters from the Explorer 33 satellite from 23 January to 23 March 1968. The TD is separated from the observed field fluctuations by calculating the direction of the plasma flow and also the direction of the minimum field fluctuation with respect to the ambient magnetic field direction.It is found that the TD is formed by the thin layered field-aligned currents (the current sheets), and that the TD is predominantly built up in the leading edge of the solar wind where the compression of the plasma and the magnetic field takes place.It is suggested that the current sheets might be locally generated in the leading edge in the turbulent conditions arising from collisions between the fast- and the slow-stream of the solar sector structure.     

Modulation of Jovian middle magnetosphere currents and auroral precipitation by solar wind-induced compressions and expansions of the magnetosphere: initial response and steady state

We present results from a theoretical model which has been used to investigate the modulation of the magnetosphere-ionosphere coupling currents in the Jovian middle magnetosphere by solar wind-induced compressions and expansions of the magnetosphere. We consider an initial system in which the current sheet field lines extend to 50R_J in the equatorial plane, and where the iogenic plasma in the current sheet undergoes steady outward radial diffusion under the influence of the ionospheric torque which tends to maintain corotation with the planet. We show using typical Jovian parameters that the upward-directed field-aligned currents flowing throughout the middle magnetosphere region in this system peak at values requiring the existence of significant field-aligned voltages to drive them, resulting in large precipitating energy fluxes of accelerated electrons and bright ‘main oval’ UV auroras. We then consider the changes in these parameters which take place due to sudden expansions or compressions of the magnetosphere, resulting from changes in the solar wind dynamic pressure. Two cases are considered and compared, these being first the initial response of the system to the change, determined approximately from conservation of angular momentum of the radially displaced plasma and frozen-in field lines, and second the subsequent steady state of steady outward radial diffusion applied to the compressed or expanded system. We show that moderate inward compressions of the outer boundary of the current sheet field lines, e.g. from 50 to 40R_J, are effective in significantly reducing the coupling currents and precipitation in the initial state, the latter then recovering, but only partly so, during the evolution to the steady state. Strong inward compressions, e.g. to 30R_J cause significant super-corotation of the plasma and a reversal in sense of the current system in the initial state, such that bright auroras may then be formed poleward of the usual ‘main auroral oval’ due to the ‘return’ currents. The sense of the currents subsequently reverts back to the usual direction as steady-state conditions are restored, but they are weak, and so is the consequent electron precipitation. For outward expansions of the current sheet, however, the field-aligned currents and electron precipitation are strongly enhanced, particularly at the poleward border mapping to the outer weak field region of the current sheet. In this case there is little evolution of the parameters between the initial expansion and the subsequent steady state. Overall, the results suggest that the Jovian middle magnetosphere coupling currents and resulting ‘main oval’ auroral acceleration and precipitation will be strongly modulated by the solar wind dynamic pressure in the sense of anti-correlation, through the resulting compressions and expansions in the size of the magnetosphere.     

Evolution of the Martian atmosphere and hydrosphere: Solar wind erosion studied by ASPERA-3 on Mars Express

The evolution of the Martian atmosphere and the potential existence of a past hydrosphere is a scientific issue of great interest in planetary research. Although the first missions to Mars had a focus on surface features and atmospheric properties, some of the missions (e.g., The Soviet Mars 2, 3 and 5) also carried instruments addressing the solar wind interaction with the Martian atmosphere and ionosphere and the potential existence of an intrinsic magnetic field on Mars. However, it took until 1989 before a spacecraft, Phobos-2, was able to carry out a more detailed investigation of the solar wind interaction with Mars. Phobos-2 gave valuable data on the Solar wind interaction with Mars during about 2 months of operations, leading to a better understanding of the solar wind impact on a weakly magnetized planet. However, Phobos-2 also raised a number of critical issues that has left science without adequate data since 1989.Investigations planned for Mars Express will cast new light on important aspects of the solar wind interaction with Mars. ASPERA-3 (Analyzer of Space Plasma and Energetic Atoms) on Mars Express will focus on the overall plasma outflow and monitor remotely the outflow and inflow of energetic neutral atoms produced by charge exchange processes. This report will discuss some of the unsolved issues about the solar wind interaction with Mars and how we plan to address these issues with Mars Express.     

Mapping Solar Wind Streams from the Sun to 1 AU: A Comparison of Techniques

A variety of techniques exist for mapping solar wind plasma and magnetic field measurements from one location to another in the heliosphere. Such methods are either applied to extrapolate solar data or coronal model results from near the Sun to 1 AU (or elsewhere), or to map in-situ observations back to the Sun. In this study, we estimate the sensitivity of four models for evolving solar wind streams from the Sun to 1 AU. In order of increasing complexity, these are: i) ballistic extrapolation; ii) ad hoc kinematic mapping; iii) 1-D upwinding propagation; and iv) global heliospheric MHD modeling. We also consider the effects of the interplanetary magnetic field on the evolution of the stream structure. The upwinding technique is a new, simplified method that bridges the extremes of ballistic extrapolation and global heliospheric MHD modeling. It can match the dynamical evolution captured by global models, but is almost as simple to implement and as fast to run as the ballistic approximation.