Computer modeling of test particle acceleration at oblique shocks

We review the basic techniques and results of numerical codes used to model the acceleration of charged particles at oblique, fast-mode, collisionless shocks. The emphasis is upon models in which accelerated particles (ions) are treated as test particles, and particle dynamics is calculated by numerically integrating along exact phase-space orbits. We first review the case where ions are sufficiently energetic so that the shock can be approximated by a planar discontinuity, and where the electromagnetic fields on both sides of the shock are defined at the outset of each computer run. When the fields are uniform and static, particles are accelerated by the scatter-free drift acceleration process at a single shock encounter. We review the characteristics of scatter-free drift acceleration by considering how an incident particle distribution is modified by interacting with a shock. Next we discuss drift acceleration when magnetic fluctuations are introduced on both sides of the shock, and compare these results with those obtained under scatter-free conditions. We describe the modeling of multiple shock encounters, discuss specific applications, and compare the model predictions with theory. Finally, we review some recent numerical simulations that illustrate the importance of shock structure to both the ion injection process and to the acceleration of ions to high energies at quasi-perpendicular shocks.     

The Solar System d/h Ratio: Observations and Theories

The measured D/H ratios in interstellar environments and in the solar system are reviewed. The two extreme D/H ratios in solar system water - (720±120)×10⁻⁶ in clay minerals and (88±11)×10⁻⁶ in chondrules, both from LL3 chondritic meteorites - are interpreted as the result of a progressive isotopic exchange in the solar nebula between deuterium-rich interstellar water and protosolar H₂. According to a turbulent model describing the evolution of the nebula (Drouart et al., 1999), water in the solar system cannot be a product of thermal (neutral) reactions occurring in the solar nebula. Taking 720×10⁻⁶ as a face value for the isotopic composition of the interstellar water that predates the formation of the solar nebula, numerical simulations show that the water D/H ratio decreases via an isotopic exchange with H₂. During the course of this process, a D/H gradient was established in the nebula. This gradient was smoothed with time and the isotopic homogenization of the solar nebula was completed in 10⁶ years, reaching a D/H ratio of 88×10⁻⁶. In this model, cometary water should have also suffered a partial isotopic re-equilibration with H₂. The isotopic heterogeneity observed in chondrites result from the turbulent mixing of grains, condensed at different epochs and locations in the solar nebula. Recent isotopic determinations of water ice in cold interstellar clouds are in agreement with these chondritic data and their interpretation (Texeira et al., 1999). This revised version was published online in June 2006 with corrections to the Cover Date.     

Can Surface Flux Transport Account for the Weak Polar Field in Cycle 23?

To reproduce the weak magnetic field on the polar caps of the Sun observed during the declining phase of cycle 23 poses a challenge to surface flux transport models since this cycle has not been particularly weak. We use a well-calibrated model to evaluate the parameter changes required to obtain simulated polar fields and open flux that are consistent with the observations. We find that the low polar field of cycle 23 could be reproduced by an increase of the meridional flow by 55% in the last cycle. Alternatively, a decrease of the mean tilt angle of sunspot groups by 28% would also lead to a similarly low polar field, but cause a delay of the polar field reversals by 1.5 years in comparison to the observations.     

Understanding magnetotail current sheet meso-scale structures using MHD simulations

Recent Cluster observations have strongly supported the existence of meso-scale structure in the magnetotail current sheet. In our study, a magnetohydrodynamic simulation event study exhibited current sheet behavior comparable to that seen in the Cluster observations. Geotail and DoubleStar observations also show that the simulation is providing a realistic representation of the magnetosphere during the period of interest; that is, when the current sheet evidently becomes bifurcated. The magnetohydrodynamic simulation allows us to place the local observations into a global contest. It shows that the observations can be explained in terms of localized reconnection tailward of the Cluster location and the formation of a flux rope nearby. The simulation also features wave-like structure across the current sheet.     

Atmospheric reentry hemisphere prediction for prograde orbits using logical disjunction

A unique logic-based algorithm for atmospheric reentry hemisphere prediction is presented for spacecraft in low-eccentricity, prograde low Earth orbits at altitudes of 300 km and lower. Using two-line element (TLE) data for initial orbit conditions, coupled with coarse estimates for spacecraft aerodynamic characteristics, the algorithm relies on logical disjunction operations based on a dual analysis of histogram and two-weighted Gaussian probability density function (PDF) fits of predicted reentry latitude data. The algorithm requires the execution of a series of parametric simulations to determine the reentry hemisphere for variations in spacecraft aerodynamic coefficients and drag reference area. When implemented, the algorithm yields accurate hemisphere predictions on average 15 days from reentry as demonstrated by historical reentry cases from 1979 to 2018. All reentry cases were selected to demonstrate the algorithm’s ability to deliver accurate reentry hemisphere predictions for spacecraft with varying physical size and mass, and reentering during different periods of solar cycle activity.     

Investigation of GPS-TEC measurements using ANN method indicating seismo-ionospheric anomalies around the time of the Chile (Mw = 8.2) earthquake of 01 April 2014

In 94 km NW of Iquique in Chile (19.610°S, 70.776°W) a powerful earthquake of M_w = 8.2 took place at 23:46:47 UTC (20:46:47 LT) on April 01, 2014. Using GPS-TEC (Total Electron Content) measurements, potential unusual variations around the time and location of the Chile earthquake have been detected based on the median and Artificial Neural Network (ANN) methods. The indices D_st, K_p, A_p and F10.7 were used to distinguish pre-earthquake anomalies from the other anomalies related to the solar-geomagnetic activities. Using the median method, striking anomalies in time series of TEC data are observed 4 days before the earthquake at 14:00 and 16:00 UTC. The ANN method detected a number of anomalies, 4 (02:00 and 16:00 UTC) and 13 (24:00 UTC) days preceding the earthquake. The results indicate that the ANN method due to its capability of non linear learning is quite promising and deserves serious attention as a robust predictor tool for seismo-ionospheric anomalies detection.