Tail Reconnection Research Articles

A “Wave‐Like” Evolution of Polar Cap Patches Modulated by Enhanced Magnetopause Reconnection and Extended Magnetotail Reconnection

AbstractA train of large‐scale polar cap patches was observed in the ionosphere during a period of southward interplanetary magnetic field (IMF) during the main phase of a moderate geomagnetic storm. The patches were sequentially detached from a dayside storm enhanced density plume through cusp dynamics due to pulsed enhancement of the dayside magnetic reconnection rate, suggested by the increase of the Kan‐Lee electric field and cross polar cap potential. After formation and IMF By duskward turning, the patches propagated obliquely like “waves” in the polar cap toward the dusk sector, and exited the polar cap over a wide magnetic local time (MLT) region (18–24 MLT). Simulations under a strongly southward IMF Bz and duskward IMF By show that the magnetotail magnetic reconnection occurs over a broad MLT range. Such a tail reconnection modulates the evolution of the patches.

A Multi‐Scale Particle‐In‐Cell Simulation of Plasma Dynamics From Magnetotail Reconnection to the Inner Magnetosphere

AbstractDuring magnetospheric substorms, plasma from magnetic reconnection in the magnetotail is thought to reach the inner magnetosphere and form a partial ring current. We simulate this process using a fully kinetic 3D particle‐in‐cell (PIC) numerical code along with a global magnetohydrodynamics (MHD) model. The PIC simulation extends from the solar wind outside the bow shock to beyond the reconnection region in the tail, while the MHD code extends much further and is run for nominal solar wind parameters and a southward interplanetary magnetic field. By the end of the PIC calculation, ions and electrons from the tail reconnection reach the inner magnetosphere and form a partial ring current and diamagnetic current. The primary source of particles to the inner magnetosphere is bursty bulk flows (BBFs) that originate from a complex pattern of reconnection in the near‐Earth magnetotail at to . Most ion acceleration occurs in this region, gaining from 10 to 50 keV as they traverse the sites of active reconnection. Electrons jet away from the reconnection region much faster than the ions, setting up an ambipolar electric field allowing the ions to catch up after approximately 10 ion inertial lengths. The initial energy flux in the BBFs is mainly kinetic energy flux from the ions, but as they move earthward, the energy flux changes to enthalpy flux at the ring current. The power delivered from the tail reconnection in the simulation to the inner magnetosphere is W, which is consistent with observations.

Global-scale magnetosphere convection driven by dayside magnetic reconnection

Plasma convection on a global scale is a fundamental feature of planetary magnetosphere. The Dungey cycle explains that steady-state convection within the closed part of the magnetosphere relies on magnetic reconnection in the nightside magnetospheric tail. Nevertheless, time-dependent models of the Dungey cycle suggest an alternative scenario where magnetospheric convection can be solely driven by dayside magnetic reconnection. In this study, we provide direct evidence supporting the scenario of dayside-driven magnetosphere convection. The driving process is closely connected to the evolution of Region 1 and Region 2 field-aligned currents. Our global simulations demonstrate that intensified magnetospheric convection and field-aligned currents progress from the dayside to the nightside within 10–20 minutes, following a southward turning of the interplanetary magnetic field. Observational data within this short timescale also reveal enhancements in both magnetosphere convection and the ionosphere’s two-cell convection. These findings provide insights into the mechanisms driving planetary magnetosphere convection, with implications for the upcoming Solar-Wind-Magnetosphere-Ionosphere Link Explorer (SMILE) mission.

Observations of ULF Pulsations During TRINNI Events

AbstractTwo magnetospheric phenomena are commonly observed in both SuperDARN (Dual Auroral Radar Network) and ground‐based magnetometer network data. These are ULF pulsations driven by magnetohydrodynamic waves, and flow bursts in the high latitude ionosphere driven by magnetic reconnection. Possible connections between these two phenomena are investigated. TRINNI (Tail Reconnection during IMF Northward, Non‐substorm Intervals) events, high‐speed ionospheric flows resulting from magnetotail reconnection during IMF Bz‐positive but By‐dominant periods, have been investigated by many researchers. These events produce measurable enhancements of the cross polar cap potential. Four previously reported TRINNI events have been analyzed using SuperDARN data, together with high‐latitude magnetometer data which show simultaneous ULF pulsations. Fourier analysis of the cross polar cap potential measured by SuperDARN and two magnetic field components revealed common spectral peaks in the 1–5 mHz range. The signals were narrow‐band filtered and complex demodulation was performed to construct the analytic signals, allowing determination of amplitude and phase. In all four events the magnetic field and potential signals showed similar wave packet structure with large‐amplitude oscillations at a common frequency for periods of the order of hours. During these periods the phase differences between the field and potential remained constant. It is argued that this phase locking suggests a causal link between the two oscillations, so that the TRINNIs could be the source of the magnetic pulsations, or they at least share a common driver.

Image‐Based Classification of Intense Radio Bursts From Spectrograms: An Application to Saturn Kilometric Radiation

AbstractSaturn Kilometric Radiation (SKR) is a non‐thermal auroral emission with peak emission occurring at 100–400 kHz. Its properties have been extensively studied since Cassini's arrival at Saturn until mission end with its Radio and Plasma Wave Science (RPWS) experiment. Low Frequency Extensions (LFEs) of SKR which consist of global intensifications of SKR accompanied by extensions of the main SKR band down to lower frequencies have been studied in particular. Low Frequency Extensions result from internally driven tail reconnection and from solar wind compressions of the magnetosphere, which also trigger tail reconnection. They have been cataloged through visual inspection with two approaches, using an intensity threshold for LFEs in 2006 (Reed et al., 2018, https://doi.org/10.1002/2017ja024499) and more recently O’Dwyer et al. (2023a, https://doi.org/10.25546/103103) produced a sample of LFEs detected by Cassini/RPWS by fitting their exact frequency‐time coordinates with polygons. In this study we use the latter catalog of LFEs as a training set for an image based machine learning algorithm to classify all LFEs detected by Cassini/RPWS. The inputs to the model are multi‐channel images consisting of spectrogram images in flux density and degree of circular polarization. The outputs of the model are binary masks showing the exact location of the LFE in frequency‐time space. The median Intersection Over Union across the testing and training set were calculated to be 0.97 and 0.98, respectively. The output of this study is a list of all 4,874 LFEs detected using this method. The list of LFE frequency‐time coordinates is available for use amongst the scientific community.

Closed field line vortices in planetary magnetospheres

ABSTRACTIn a rotation-dominated magnetosphere, there is a region where closed field lines rotate around the planet, and also a region where the open field lines stretch away from the planet, forming the lobes of the magnetotail. This paper shows that there could be a third, significantly different region, where the closed field lines form twisted vortex structures anchored in the magnetotail. Such patterns form when there are significant plasma sources inside the magnetosphere and the time-scale of the plasmoid formation process is substantially larger than the planetary rotation period. In the presence of vortices, the Dungey and Vasyliunas cycles act differently. The Dungey flow does not penetrate the central region of the polar cap. Tail reconnection events are rare, thus leaving the plasma time enough to participate in the essentially 3D vortex-forming plasma motion. The above conditions are fulfilled for Saturn. We discovered vortex-like patterns in the plasma and magnetic field data measured by the Cassini spacecraft in the nightside magnetosphere of Saturn. The plasma whirling around in these vortices never reaches the dayside, instead, it performs a retrograde motion in the high latitude regions of the magnetotail. Low-energy plasma data suggest that the observed patterns correspond to the closed field line vortices.

Applying Magnetic Curvature to MMS Data to Identify Thin Current Sheets Relative to Tail Reconnection.

Thin current sheets (TCSs) have been postulated to be a necessary precondition for reconnection onset. Magnetic reconnection X-lines in the magnetotail have been observed to be more common duskward of midnight. We take advantage of the MMS tetrahedral formation during the 2017-2020 MMS tail seasons to calculate the thickness of the cross-tail neutral sheet relative to ion gyroradius. While a similar technique was applied to Cluster data, current sheet thickness over a broader range of radial distances has not been robustly explored before this study. We compare our analysis to recent theories regarding mechanisms of tail current sheet thinning and to recent simulations. We find MMS spent more than twice as long in ion-scale TCSs in the pre-midnight sector than post-midnight, despite nearly even plasma sheet dwell time. The dawn-dusk asymmetry in the distribution of Ion Diffusion Regions, as previously reported in relation to regions of TCSs, is also analyzed.

The Role of Magnetic Flux Rope in Ion Acceleration: MHD Simulations and Test-particle Tracing

Magnetic flux ropes (MFRs), playing a crucial role in particle energization and energy transport in the solar–terrestrial space, are helical structures produced by magnetic reconnection. It has been both theoretically predicted and observationally confirmed that MFRs and associated processes are inherently three-dimensional in space. Although such structures have been widely suggested as a favorable place for electron acceleration, whether large-scale MFRs can lead to ion acceleration has been rarely investigated. In this study, an MHD model is used to examine the evolution of large-scale MFRs in the magnetotail, and a test-particle simulation is further employed to study the associated ion energization. Results show that magnetic reconnections take place at multiple X-lines in the magnetotail current sheet, generating a twisted MFR with a scale of about 10 R e in azimuth. It propagates earthward following the tail reconnection but its east and west wings are significantly distorted azimuthally. Test-particle tracing reveals that ions (0.1–100 keV) can be trapped within the rope while being effectively accelerated. The rope therefore brings in energetic plasma sources into the inner magnetosphere as it transports earthward. These results demonstrate that the MFR is an important source carrier for the ring-current formation in the inner magnetosphere.

Evolution of IMF By induced asymmetries during substorms: Superposed epoch analysis at geosynchronous orbit

The By component of the magnetic field inside the magnetosphere is positively correlated with the By component of the Interplanetary Magnetic Field (IMF). This leads to asymmetries in aurora, plasma convection and electric currents between the northern and southern hemispheres It has been demonstrated that magnetic conjugate locations in the northern and southern ionosphere become less displaced during magnetospheric substorms, which are associated with enhanced reconnection in the near-Earth tail. Here we directly address how the average By component in the magnetotail evolves relative to substorm onset by performing a superposed epoch analysis of the magnetic field observed at nightside geosynchronous orbit during periods with dominant IMF By. The observations demonstrate that the average |By| in the magnetotail increases during the loading phase prior to onset. |By| maximizes in the expansion phase and is subsequently reduced during the remaining unloading phase. The observed trends become more pronounced using substorm onset lists that on average identify stronger substorms. Since dayside reconnection dominates over tail reconnection during the loading phase, whereas tail reconnection dominates during the unloading phase, the results demonstrate how asymmetries build up during periods with low tail reconnection and are reduced during periods with enhanced tail reconnection in agreement with previous case studies of conjugate auroral substorm features.

Global Magnetohydrodynamic Magnetosphere Simulation With an Adaptively Embedded Particle‐In‐Cell Model

AbstractWe perform a geomagnetic event simulation using a newly developed magnetohydrodynamic with adaptively embedded particle‐in‐cell (MHD‐AEPIC) model. We have developed effective criteria to identify reconnection sites in the magnetotail and cover them with the PIC model. The MHD‐AEPIC simulation results are compared with Hall MHD and ideal MHD simulations to study the impacts of kinetic reconnection at multiple physical scales. At the global scale, the three models produce very similar SYM‐H and SuperMag Electrojet indexes, which indicates that the global magnetic field configurations from the three models are very close to each other. We also compare the ionospheric solver results and all three models generate similar polar cap potentials and field‐aligned currents. At the mesoscale, we compare the simulations with in situ Geotail observations in the tail. All three models produce reasonable agreement with the Geotail observations. At the kinetic scales, the MHD‐AEPIC simulation can produce a crescent shape distribution of the electron velocity space at the electron diffusion region, which agrees very well with MMS observations near a tail reconnection site. These electron scale kinetic features are not available in either the Hall MHD or ideal MHD models. Overall, the MHD‐AEPIC model compares well with observations at all scales, it works robustly, and the computational cost is acceptable due to the adaptive adjustment of the PIC domain. It remains to be determined whether kinetic physics can play a more significant role in other types of events, including but not limited to substorms.

The Response of Saturn's Dawn Field‐Aligned Currents to Magnetospheric and Ring Current Conditions During Cassini's Proximal Orbits: Evidence for a Region 2 Response at Saturn

AbstractCassini's 2017 proximal orbits provided the opportunity to examine the auroral field‐aligned currents in the northern hemisphere dawn sector in relation to wider magnetospheric conditions. We combine three recent studies to examine the response of the dawn region auroral field‐aligned currents and the azimuthal ring currents to compressions and expansions of the Saturnian magnetosphere. For compressions of Saturn's magnetosphere resulting in tail reconnection, the currents within the downward current sheet, located equatorward of the main auroral oval, increases in strength with increasing total ring current and location of the peak downwards current moves inwards toward Saturn. While the inverse relation occurs during intervals of quiet or expanded magnetospheric conditions. During compression events there is an increase in the energetic particle intensities, in particular in the protons (35–506 keV), within the downward current region. This current system is akin to an Earth‐like “region 2” field aligned current within Saturn's magnetosphere, with tail reconnection occurring when the magnetosphere is compressed resulting in a partial nightside ring current closed by a downward current near to dawn. Within the upward current sheet, mapping to Saturn's main auroral oval, both non‐rotating subcorotating current and the rotating Planetary Period Oscillations (PPOs) currents flow. The upward current is strongly modulated by the PPOs but also increases in strength, with enhanced high‐energy protons, during intervals of magnetospheric compressions and tail reconnection. We conclude that the enhanced plasma injected into the midnight‐dawn sector during tail reconnection events results in an enhanced subcorotation current system.

Ionospheric plasma flows associated with the formation of the distorted nightside end of a transpolar arc

Abstract. We investigate ionospheric flow patterns occurring on 28 January 2002 associated with the development of the nightside distorted end of a J-shaped transpolar arc (nightside distorted TPA). Based on the nightside ionospheric flows near to the TPA, detected by the SuperDARN (Super Dual Auroral Radar Network) radars, we discuss how the distortion of the nightside end toward the pre-midnight sector is produced. The J-shaped TPA was seen under southward interplanetary magnetic field (IMF) conditions, in the presence of a dominant dawnward IMF-By component. At the onset time of the nightside distorted TPA, particular equatorward plasma flows at the TPA growth point were observed in the post-midnight sector, flowing out of the polar cap and then turning toward the pre-midnight sector of the main auroral oval along the distorted nightside part of the TPA. We suggest that these plasma flows play a key role in causing the nightside distortion of the TPA. SuperDARN also found ionospheric flows typically associated with Tail Reconnection during IMF Northward Non-substorm Intervals (TRINNIs) on the nightside main auroral oval, before and during the TPA interval, indicating that nightside magnetic reconnection is an integral process to the formation of the nightside distorted TPA. During the TPA growth, SuperDARN also detected anti-sunward flows across the open–closed field line boundary on the dayside that indicate the occurrence of low-latitude dayside reconnection and ongoing Dungey cycle driving. This suggests that nightside distorted TPA can grow even in Dungey-cycle-driven plasma flow patterns.

The EDR inflow region of a reconnecting current sheet in the geomagnetic tail

On 6 July 2017, the four Magnetospheric Multiscale spacecrafts were positioned within an electron diffusion region (EDR) just northward of a reconnection X line. The EDR was identified by electron crescent distributions, out-of-plane current, and energy conversion. From this position, the three spacecrafts closest to the X line (within about three electron inertial lengths) were able to accurately measure the reconnection electric field and the electron inflow velocity. The reconnection rates derived from the electric field and inflow velocity measurements agree with theoretical estimates (0.11–0.17) and a previous measurement of EM in a tail reconnection event on 11 July 2017.

Time‐Varying Magnetopause Reconnection During Sudden Commencement: Global MHD Simulations

AbstractIn response to a solar wind dynamic pressure enhancement, the compression of the magnetosphere generates strong ionospheric signatures and a sharp variation in the ground magnetic field, termed sudden commencement (SC). While such compressions have also been associated with a contraction of the ionospheric polar cap due to the triggering of reconnection in the magnetotail, the effect of any changes in dayside reconnection is less clear and is a key component in fully understanding the system response. In this study, we explore the time‐dependent nature of dayside coupling during SC by performing global simulations using the Gorgon magnetohydrodynamic code and impact the magnetosphere with a series of interplanetary shocks with different parameters. We identify the location and evolution of the reconnection region in each case as the shock propagates through the magnetosphere, finding strong enhancement in the dayside reconnection rate and prompt expansion of the dayside polar cap prior to the eventual triggering of tail reconnection. This effect pervades for a variety of interplanetary magnetic field orientations, and the reconnection rate is most enhanced for events with higher dynamic pressure. We explain this by repeating the simulations with a large explicit resistivity, showing that compression of the magnetosheath plasma near the propagating shock front allows for reconnection of much greater intensity and at different locations on the dayside magnetopause than during typical solar wind conditions. The results indicate that the dynamic behavior of dayside coupling may render steady models of reconnection inaccurate during the onset of a severe space weather event.

A Statistical Study on Wide‐Amplitude Kinetic Alfvénic Pulse at 8–12RE in the Magnetotail by THEMIS Spacecraft From 2008 to 2010

AbstractUsing THEMIS spacecraft observations from 2008 to 2010, we perform statistical study on the wide‐amplitude kinetic Hall‐like Alfvénic wave (KAW) pulse in the near‐Earth tail region of 8–12 RE. The KAW pulse is identified by Ez‐By pulse with parallel‐predominantly Poynting flux (|P///PT| > 0.7) and large Ez/By ratio (|∆Ez/∆By | > VA, VA is Alfvénic velocity). Three typical events are presented. Therein, the KAW pulse on 24 March 2009 at the magnetic dipolarization is accompanied by the ion/electron dispersed structure (IDS/EDS; time of flight effect from reconnection source). Statistical study shows that the near‐Earth KAW pulse has a typical temporal scale of ∼15–25 s. The KAW pulses are mainly distributed in the plasma sheet boundary layer (PSBL) and tail lobe regions of β < 0.1 (β is the ratio of thermal pressure to magnetic pressure). The KAW pule at the PSBL/lobe has a wide distribution in the dawn‐dusk direction (−10 RE < |Y| < 10 RE), with the maximum occurrence at the midnight (|Y| < 4 RE). The KAW pulses at the PSBL/lobe are predominantly earthward‐propagating. Only a small number of the KAW pulses are tailward‐propagating. The earthward‐propagating pulse has a higher Poynting flux (PT) than tailward‐propagating pulse. The KAW pulse occurs during quiet time (AE < 200 nT) as well as substorm interval (AE > 200 nT). The substorm‐interval pulse has a higher Poynting flux than the quiet‐time pulse. Finally, tail reconnection is likely the significant source of the near‐Earth KAW pulses at 8–12 RE.

Global MHD Simulation of the Weak Southward IMF Condition for Different Time Resolutions

We performed high-resolution three-dimensional global MHD simulations to determine the impact of weak southward interplanetary magnetic field (IMF) (Bz = −2 nT) and slow solar wind to the Earth’s magnetosphere and ionosphere. We considered two cases of differing, uniform time resolution with the same grid spacing simulation to find any possible differences in the simulation results. The simulation results show that dayside magnetic reconnection and tail reconnection continuously occur even during the weak and steady southward IMF conditions. A plasmoid is generated on closed plasma sheet field lines. Vortices are formed in the inner side of the magnetopause due to the viscous-like interaction, which is strengthened by dayside magnetic reconnection. We estimated the dayside magnetic reconnection which occurred in relation to the electric field at the magnetopause and confirmed that the enhanced electric field is caused by the reconnection and the twisted structure of the electric field is due to the vortex. The simulation results of the magnetic field and the plasma properties show quasi-periodic variations with a period of 9–11 min between the appearances of vortices. Also the peak values of the cross-polar cap potential are both approximately 50 kV, the occurrence time of dayside reconnections are the same, and the polar cap potential patterns are the same in both cases. Thus, there are no significant differences in outcome between the two cases.

Simultaneous UV Images and High‐Latitude Particle and Field Measurements During an Auroral Dawn Storm at Jupiter

AbstractWe present multi‐instrument Juno observations on day‐of‐year 86, 2017 that link particles and fields in Jupiter's polar magnetosphere to transient UV emissions in Jupiter's northern auroral region known as dawn storms. Juno ranged from 42°N to 51°N in magnetic latitude and 5.8–7.8 Jovian radii (1 RJ = 71,492 km) during this period. These dawn storm emissions consisted of two separate, elongated structures which extended into the nightside, rotated with the planet, had enhanced brightness (up to at least 1.4 megaRayleigh) and high color ratios. The color ratio is a proxy for the atmospheric penetration depth and therefore the energy of the electrons that produce the UV emissions. Juno observed electrons and ions on magnetic field lines mapping to these emissions. The electrons were primarily field‐aligned, bidirectional, and, at times, exhibited sudden intensity decreases below ∼10 keV coincident with intensity enhancements up to energies of ∼1,000 keV, consistent with the high color ratio observations. The more energetic electron distributions had characteristic energies of ∼160–280 keV and downward energy fluxes (∼70–135 mW m−2) that were a significant fraction needed to produce the UV emissions for this event. Magnetic field perturbations up to ∼0.7% of the local magnetic field showing evidence of upward and downward field‐aligned currents, whistler mode waves, and broadband kilometric radio emissions were also observed along Juno's trajectory during this time frame. These high‐latitude observations show similarities to those in the equatorial magnetosphere associated with dynamics processes such as interchange events, plasma injections, and/or tail reconnection.

Evolution of IMF By Induced Asymmetries: The Role of Tail Reconnection

AbstractNorth‐south asymmetries arise in the magnetosphere‐ionosphere system when a significant east‐west (By) component is present in the interplanetary magnetic field (IMF). During such conditions, a By component with the same sign as the IMF By component is induced in the magnetosphere, and the locations of conjugate magnetic footpoints are displaced between the two hemispheres. It has been suggested that these asymmetries are introduced into the closed magnetosphere by tail reconnection. However, recent studies instead suggest that asymmetric lobe pressure induces the asymmetries, which are then reduced during periods of enhanced tail reconnection. To address this, we use the Lyon‐Fedder‐Mobarry (LFM) model and initiate a loading‐unloading cycle in multiple runs by changing the IMF. Asymmetries are induced during the loading phase and reduced during the unloading phase. The model results thus suggest that asymmetries arise during periods with low tail reconnection and are reduced during periods with enhanced tail reconnection.

Changes in the Magnetic Field Topology and the Dayside/Nightside Reconnection Rates in Response to a Solar Wind Dynamic Pressure Front: A Case Study

AbstractOne of the most significant observations associated with a sharp enhancement in solar wind dynamic pressure, , is the poleward expansion of the auroral oval and the closing of the polar cap. The polar cap shrinking over a wide range of magnetic local times (MLTs), in connection with an observed increase in ionospheric convection and the transpolar potential, led to the conclusion that the nightside reconnection rate is significantly enhanced after a pressure front impact. However, this enhanced tail reconnection has never been directly measured. We demonstrate the effect of a solar wind dynamic pressure front on the polar cap closure, and for the first time, measure the enhanced reconnection rate in the magnetotail, for a case occurring during southward background Interplanetary Magnetic Field (IMF) conditions. We use Polar Ultra‐Violet Imager (UVI) measurements to detect the location of the open‐closed field line boundary, and combine them with Assimilative Mapping of Ionospheric Electrodynamics (AMIE) potentials to calculate the ionospheric electric field along the polar cap boundary, and thus evaluate the variation of the dayside/nightside reconnection rates. We find a strong response of the polar cap boundary at all available MLTs, exhibiting a significant reduction of the open flux content. We also observe an immediate response of the dayside reconnection rate, plus a phased response, delayed by ∼15–20 min, of the nightside reconnection rate. Finally, we provide comparison of the observations with the results of the Open Geospace General Circulation Model (OpenGGCM), elucidating significant agreements and disagreements.

Acceleration of Oxygen Ions in Dipolarization Events: 1. CPS Distributions

AbstractUsing an MHD simulation of near tail reconnection associated with a flow burst and the collapse (dipolarization) of the inner tail in combination with test particle tracing we study the acceleration and flux increases of energetic oxygen ions (O+). The characteristic orbits, distributions, and acceleration mechanisms are governed by the dimensionless parameter σ = ωcitn, where ωci is the ion gyro frequency and tn a characteristic Alfvén time of the MHD simulation. For σ < 1, central plasma sheet (CPS) populations after the passage of the dipolarization front are found to resemble half‐shells in velocity space oriented toward dusk. They originate from within the CPS and are energized typically by a single encounter of the region of enhanced cross‐tail electric field associated with the flow burst. For larger σ values (σ > 1) the O+ distributions resemble more closely those of protons, consisting of two counter‐streaming field‐aligned beams and an, albeit more tenuous and irregular, ring population perpendicular to the magnetic field. The existence of the beams, however, depends on suitable earthward moving source populations in the plasma sheet boundary layer or the adjacent lobes. The acceleration to higher energies is found to indicate a charge dependence, consistent with a dominance of more highly charged ions at energies of a few hundred keV. As in earlier simulations, the simulated fluxes show large anisotropies and nongyrotropic effects, phase bunching, and spatially and temporally localized beams.