Engineering topics
Song, P.
Publications and source records attributed to Song, P..
Surface Penetrating Radar Simulations for Europa
The space environment above the icy surface of Europa is a source of radio noise in this frequency range from natural sources in the Jovian magnetosphere. The ionospheric and magnetospheric plasma environment of Europa affects propagation of transmitted and return signals between the spacecraft and the solid surface in a frequency-dependent manner. The ultimate resolution of the subsurface sounding measurements will be determined, in part, by a capability to mitigate these effects. We discuss an integrated multi-frequency approach to active radio sounding of the Europa ionospheric and local magnetospheric environments, based on operational experience from the Radio Plasma Imaging @PI) experiment on the IMAGE spacecraft in Earth orbit, in support of the subsurface measurement objectives.
Radio Sounding Science at High Powers
Future space missions like the Jupiter Icy Moons Orbiter (JIMO) planned to orbit Callisto, Ganymede, and Europa can fully utilize a variable power radio sounder instrument. Radio sounding at 1 kHz to 10 MHz at medium power levels (10 W to kW) will provide long-range magnetospheric sounding (several Jovian radii) like those first pioneered by the radio plasma imager instrument on IMAGE at low power (less than l0 W) and much shorter distances (less than 5 R(sub E)). A radio sounder orbiting a Jovian icy moon would be able to globally measure time-variable electron densities in the moon ionosphere and the local magnetospheric environment. Near-spacecraft resonance and guided echoes respectively allow measurements of local field magnitude and local field line geometry, perturbed both by direct magnetospheric interactions and by induced components from subsurface oceans. JIMO would allow radio sounding transmissions at much higher powers (approx. 10 kW) making subsurface sounding of the Jovian icy moons possible at frequencies above the ionosphere peak plasma frequency. Subsurface variations in dielectric properties, can be probed for detection of dense and solid-liquid phase boundaries associated with oceans and related structures in overlying ice crusts.
Empirical Models of the Plasma Density in the Inner Magnetosphere
The radio plasma imager (RPI) on the IMAGE satellite performs radio sounding in the magnetosphere, transmitting coded signals stepping through the frequency range of interest and receiving the returned echoes. The measurements provide the echo amplitude as a function of frequency and echo delay time on a so-called plasmagram. A newly developed algorithm inverts THE echo traces on a plasmagram to electron density spatial distributions. Rased on these observed density distributions, an empirical model is constructed to describe the two-dimensional density distribution in the inner magnetosphere.
Remote Radio Sounding Science for JIMO
Radio sounding of the Earth's top side ionosphere and magnetosphere is a proven technique from geospace missions such as the International Satellites for Ionospheric Studies (ISIS) and the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE). Application of this technique to the Jupiter Icy Moons Orbiter (JIMO) mission will provide unique remote sensing observations of the plasma and magnetic field environments, and the subsurface conductivities, of Europa, Ganymede, and Callisto. Spatial structures of ionospheric plasma above the moon surfaces vary in response to magnetic field perturbations from (1) magnetospheric plasma flows, (2) ionospheric currents from ionization of sputtered surface material, and (3) induced electric currents in salty subsurface oceans. Radio sounding at 3 kHz to 10 MHz can provide globally-determined electron densities necessary for the extraction of the oceanic current signals and supplements in-situ plasma and magnetic field measurements. Subsurface variations in conductivity, can be investigated by radio sounding from 10 MHz to 40 MHz allowing the determination of the presence of dense and solid-liquid phase boundaries associated with oceans and related structures in overlying ice crusts.
Plasmaspheric Mass Loss and Refilling as a Result of a Magnetic Storm
Using the sounding measurements from the radio plasma imager on IMAGE, and a plasma density inversion algorithm, we derive the plasma density profiles along the magnetic field in a few L-shells every 14 hours at magnetic local noon before, during, and after the March 31,2001 magnetic storm. An empirical model of the plasmaspheric plasma density distribution is derived as a reference using the measurements before the storm. During the storm, the equatorial plasma was substantially depleted in a range of L-shells. The flux tubes were refilled after the storm. The filling ratio, the equatorial plasma density normalized by its quiet-time value before the storm, is introduced to assess the time evolution of the depletion and refilling processes. The depletion, more than two thirds of the quiet time content, appeared to occur rather quickly after the storm onset, as determined by the limited temporal resolution of the measurements. The refilling proceeded, although more slowly than the depletion process, significantly faster than the theoretical prediction of a 3-day time scale. Dynamic structures are observed in situ and confirmed by the extreme ultraviolet imager (EUV) measurements.
Surface Penetrating Radar Simulations for Jupiter's Icy Moons
The icy moons of Jupiter (Europa, Callisto, and Ganymede) are of similar overall composition but show different surface features as a result of different sub-surface processes. Furthermore, each of these moons could have a liquid ocean of water buried underneath the icy crust, but their depth can only be speculated. For Europa, estimates put the thickness of the ice shell anywhere between 2-30 km, with'a few models predicting up to 100 km. Much of the uncertainties are due to the largely unknown temperature gradients and levels of water impurities across different surface layers. One of the most important geological processes is the possible transportation of heat by ice convection. If the ice is convecting, then an upper limit of about 20 km is set for the depth of the ocean underneath. Convection leads to a sharp increase in temperature followed by a thick region of nearly constant temperature. If ice is not convecting, then an exponentially increasing temperature profile is expected. The crust is thought to be a mixture of ice and rock, and although the exact percentage of rock is not known, it is expected to be low. Additionally, the ice crust could contain salt, similar to sea ice on Earth. The exact amount of salt and how that amount changes with depth is also unknown. In preparation for the Jupiter Icy Moons Orbiter (JIMO) mission, we performed simulations for a surface-penetrating radar investigating signatures for different possible surface and sub-surface structures of these moons in order to estimate the applicability of using radar with a frequency range between 1 and 50 MHz. This includes simulations of power requirements, attenuation losses, layer resolutions for scenarios with and without the presence of a liquid ocean underneath the ice, cases of convecting and non-convecting ice, different impurities within the ice, and different surface roughnesses.
POLAR Magnetosheath Observations on May 4,1998
The unusually high solar wind pressure and strongly southward IMF on May 4, 1998, pushed the magnetopause well into the geosynchronous orbit which exposed the POLAR satellite to the magnetosheath and solar wind. We use a gasdynamic convected field model to predict the magnetosheath quantities and then compare them with the in situ observations. The model prediction helps to reduce the uncertainty in the timing of the solar wind arrival time and provides a reference value for each physical parameter. It also helps to resolve the location of the satellite during strong magnetic fluctuations near the magnetopause. The plasma measurements from the TIDE instrument, in conjunction with the magnetometer measurements, indicate that there is a magnetospheric boundary layer during the event. There are also transient signatures near the magnetopause which may be caused by magnetospheric flux transfer events.
Anisotropic shock jump conditions: Theory and observations
The MHD Rankine-Hugoniot (RH) relations for shock waves in a collisionless plasma with bi-Maxwellian distribution functions are considered. While by introducing the pressure anisotropy parameter xi in the RH relations, the number of unknowns -- B, V, n, p and xi (a total of 9) -- becomes one more than the total number of the conservation equations, it is possible to use the observed quantities on both sides of the shock to study the anisotropy changes across the shock. A simple relation for the anisotropy change across the shock is derived as a function of the ratio of magnetic fields m(= B'/B), the shock normal angle theta(sub Bn) and the plasma beta and beta' (primes are downstream values). Since m and theta(sub Bn) can be determined accurately in observation, the reliability of the anisotropy change deduced is mostly dependent on the accuracy of the measurements beta and beta'. We have applied the results to six low-beta quasi-perpendicular (Q perpendicular) laminar bow shock crossings with temperature anisotropy measured in the magnetosheath. In the six test cases, it is found that the predicted pressure anisotropies agree well with those observed in the magnetosheath.
Identification of low-frequency fluctuations in the terrestrial magnetosheath
On the basis of magnetohydrodynamic (MHD) theory we develop a scheme for distinguishing among the four low-frequency modes which may propagate in a high-beta anisotropic plasma such as the magnetosheath: the fast and slow magnetosonic, the Alfven, and mirror modes. We use four parameters: the ratio of transverse to compressional powers in the magnetic field, the ratio of the wave powers in the thermal pressure and in the magnetic field, the ratio of the perturbations in the thermal and magnetic pressures, and the ratio of the wave powers in the velocity and in the magnetic field. In the test case of an Active Magnetospheric Particle Tracer Explorers/Ion Release Module (AMPTE/IRM) magnetosheath pass near the Sun-Earth line downstream of a quasi-perpendicular shock, the four modes can be clearly distinguished both spatially and spectrally. Near the bow shock, the waves are Alfvenic in a large frequency range, 1 to 100 mHz. In the middle and inner magnetosheath, the waves below 10 mHz are Alfvenic. The fast mode waves occur in the higher-frequency end of the enhanced spectrum, 80 mHz for the middle magnetosheath and 55 mHz for the inner sheath. The wave enhancement in the intermediate frequencies is slow modes in the inner sheath and mirror modes in the middle sheath. This confirms the earlier report of the existence of the slow mode waves near the magnetopause. These slow waves provide evidence that the magnetopause is an active source of the waves in the sheath. We also show that the measured frequency of a wave is close to an invariant if the magnetosheath flow is in a steady state. Therefore changes in the frequencies of enhanced waves indicate emergence, or damping, or mode conversion of the waves.
Modelling the low-latitude boundary layer with reconnection entry
We develop a one-dimensional Low Latitude Boundary Layer (LLBL) model for northward interplanetary magnetic field (IMF). The boundary layer in this model is uniform in the direction normal to the magnetopause, a 'plateau-type' boundary layer. The boundary layer motion is decoupled from the magnetosheath motion and driven by the plasma pressure associated with the incoming solar wind plasma near local noon, which has become entrained on closed field lines as a result of reconnection in the cusp region. Dissipation in the ionosphere at the feet of the boundary layer field lines opposes this motion. There are two physical solutions for the model. In one, the boundary layer reaches a terminal velocity in the tail as the boundary layer plasma effectively joins the solar wind flow. In the other solution, the flow is nearly stopped in the far tail. In combination with other mechanisms, this latter solution may correspond to the case in which the boundary layer plasma participates in magnetospheric convection and returns sunward. The density, velocity, and thickness as functions of distance from local noon are studied, assuming that the magnetopause hasa elliptical shape and the magnetospheric field is dipolar.
Structure and properties of the subsolar magnetopause for northward interplanetary magnetic field - Multiple-instrument particle observations
The paper examines the structure and properties of the subsolar magnetopause for northward IMF on the basis of measurements from 10 different instrument for three ISEE crossings. It is shown that the overall structure and properties are similar for the three crossings, indicating that the magnetopause is relatively well determined in the subsolar region for strongly northward IMF. The combined data set suggests that the magnetopause region is best organized by defining a sheath transition layer and steplike boundary layers. The electron flux enhancements in the lowest energies in the boundary layers and magnetosphere are found to be ionospheric electrons and not photoelectrons from the spacecraft. For northward IMF, they are photoelectrons, but for southward IMF they may be secondary electrons. The density measurements from differential and integral techniques are similar, leaving no room for a significant 'invisible' population.
Wave properties near the subsolar magnetopause - Pc 1 waves in the sheath transition layer
We study the waves in the frequency range of Pc 1 in the sheath transition layer of the magnetopause from the ISEE 1 and 2 observations. The waves are enhanced in the sheath transition layer, although they are scattered into the magnetosheath when the outer edge of the sheath transition layer is not sharp. The wave frequency is proportional to, and equal to, about 44 percent of the ion gyrofrequency. The waves are left-handed polarized for southward IMF, but linearly polarized for northward IMF. The direction of maximum variation is perpendicular to both the background field and the gradients of the field and density for northward IMF; for southward IMF, the waves are more turbulent. Wave generation mechanisms may depend on IMF orientations rather than the shock geometry. To investigate the free energy generating the waves for northward IMF, a method is developed combining the measurements from the fast plasma experiment and Lepedea to obtain a high time resolution estimate of the temperature anisotropy for strongly northward IMF. The estimated ion temperature anisotropy is enhanced, up to a factor of 2, within the sheath transition layer for northward IMF.
Wave properties near the subsolar magnetopause - Pc 3-4 energy coupling for northward interplanetary magnetic field
Strong slow mode waves in the Pc 3-4 frequency range are found in the magnetosheath close to the magnetopause. We have studied these waves at one of the ISEE subsolar magnetopause crossings using the magnetic field, electric field, and plasma measurements. We use the pressure balance at the magnetopause to calibrate the Fast Plasma Experiment data versus the magnetometer data. When we perform such a calibration and renormalization, we find that the slow mode structures are not in pressure balance and small scale fluctuations in the total pressure still remain in the Pc 3-4 range. Energy in the total pressure fluctuations can be transmitted through the magnetopause by boundary motions. The Poynting flux calculated from the electric and magnetic field measurements suggests that a net Poynting flux is transmitted into the magnetopause. The two independent measurements show a similar energy transmission coefficient. The transmitted energy flux is about 18 percent of the magnetic energy flux of the waves in the magnetosheath. Part of this transmitted energy is lost in the sheath transition layer before it enters the closed field line region. The waves reaching the boundary layer decay rapidly. Little wave power is transmitted into the magnetosphere.
Waves in the inner magnetosheath - A case study
We study the waves within, and upstream of, the density enhancements in front of the magnetopause using magnetic field and plasma measurements from ISEE 1 and 2. The waves upstream of the density enhancements are most likely to be mirror modes and are convected with the magnetosheath flow. Within the density enhancements there is additional wave power at lower frequencies. These low frequency waves appear to be slow modes propagating sunward and quasi-standing in the flow. These observations lend support to the contention that the outer edge of the density enhancement is a slow mode wave front which marks the greatest distance that slow modes can propagate upstream into the solar wind.
Slow mode transition in the frontside magnetosheath
Three magnetosheath passes with density enhancements in front of the magnetopause are studied with data from ISEE 1, 2, and 3. The density structure appears to be locally generated and slow mode in nature. In one pass when ISEE 1 and 2 were well separated, the motion of the density structure can be determined. The density structure appears to stand in the magnetosheath flow. Thus, it propagates upstream in the rest frame of the flow. The flow in and near the density structure appears to be closer to isothermal than adiabatic. The flow velocity decreases from super-slow to being close to the intermediate and slow mode velocities at the outer edge of the density structure. This study provides additional evidence that the density structure in front of the magnetopause is a slow mode transition in which the flow velocity decreases to the MHD slow mode velocity. The slow mode transition may consist of two wave fronts and a region with strong slow mode waves. This slow mode transition may play an important role in establishing the flow and field pattern near the magnetopause.
Model of the formation of the low-latitude boundary layer for strongly northward interplanetary magnetic field
A model for the formation of the low-boundary layer in which a magnetosheath flux tube reconnects in the north and south beyond the cusp when the IMF is strongly northward is presented. For northward IMF the geomagnetic field captures solar wind flux tubes through intermittent reconnection at the cusp region. These newly captured flux tubes shorten and sink into the magnetosphere while the flux tubes reorient themselves as they become assimilated into the magnetosphere. There is no significant acceleration for the particles within the flux tubes. In the magnetosphere the interchange instability disperses the flux tube azimuthally along the magnetopause to form a boundary layer. Subsequent reconnection forms sublayers of the boundary layer and different sublayers represent different ages after reconnection. The interchange instability is stable radially, which keeps sharp boundaries between sublayers and between the boundary layer and magnetosphere.
Solar cycle variations in the size and shape of the magnetopause
The 10 years of the ISEE 1 and 2 mission covering much of solar cycle 21 and the beginning of solar cycle 22 make it possible to study the position, shape, and motion of the magnetopause throughout the course of changing solar activity. The size and shape of the magnetopause were determined for each observing season using the ISEE 1 and 2 magnetometer data IMP 8 data were used to monitor the solar wind changes with the solar cycle. During the 1979-1980 season, at solar maximum, the solar wind dynamic pressure was at its lowest values, and, at solar minimum, the solar wind pressure was at its largest values, more than double the value in the 1979-1980 season. During this solar cycle, the magnetopause was about 0.5 R(E) farther when the interplanetary magnetic field (IMF) was strongly northward, than when strongly southward. Both standoff distance values are fround to be smaller than the value found by Fairfield (1971). The standoff distance of the magnetopause for northward IMF is anticorrelated with the solar wind pressure. However, the standoff distance for southward IMF seems relatively insensitive to solar wind dynamic pressure.