HR: 0800h
AN: SM41A-1102    [Abstracts]
TI: Preliminary Simulation Results for Stormtime Ring Current in a Self-Consistent Magnetic Field Model
AU: * Liu, S
EM: hanzo@atmos.ucla.edu
AF: Dept. of Atmospheric and Oceanic Sciences, UCLA, 405 Hilgard Ave., Los Angeles, CA 90095 United States
AU: Chen, M W
EM: margaret.w.chen@aero.org
AF: Space Science Applications Laboratory, Aerospace Corporation, 2350 E. El Segundo Blvd., El Segundo, CA 90245 United States
AU: Schulz, M
EM: mike.schulz@lmco.com
AF: Space Sciences Department, Lockheed Palo Alto Research Laboratory, 3251 Hanover St., Palo Alto, CA 94304 United States
AU: Lyons, L R
EM: larry@atmos.ucla.edu
AF: Dept. of Atmospheric and Oceanic Sciences, UCLA, 405 Hilgard Ave., Los Angeles, CA 90095 United States
AB: The stormtime ring current generates a strong and time-dependent perturbation of the magnetospheric $\vec{B}$ field, and this magnetic-field perturbation can have important feedback on the dynamics of ring current particles themselves. In particular, the modification of $\vec{B}$ can significantly alter the gradient-curvature drifts of ring current particles, and the induced electric field associated with $\partial\vec{B}/\partial t$ can inhibit ring current particle injection and energization. Thus, in order to accurately simulate the stormtime ring current, we need a self-consistent magnetic field model that takes into account effects of the ring current on the particles that produce it. This study is our first attempt to address this issue. We assume for simplicity a model for $\vec{B}$ ($= \nabla\alpha\times\nabla\beta$) such that magnetic field lines lie in meridional planes and satisfy the generic equation $r = La(1+0.5r^3/b^3)\sin^2\theta$, where $r$ is the radial distance from the point dipole, $\theta$ is the magnetic colatitude, $a$ is the radius of the Earth, $L$ is a dimensionless field-line label inversely proportional to the Euler potential $\alpha = -\mu_E/La$, $\beta$ is the magnetic local time, and the parameter $b$ (a function of $L$ and $\beta$) controls the amount by which a field line is stretched. The special case of constant $b$ yields Dungey's model magnetosphere (dipole field plus uniform southward $\Delta B$), in which the limit $b\rightarrow\infty$ corresponds to a purely dipolar $B$ field ($\Delta B = 0$). More generally, we now let the value of $b$ varies from field line to field line so as to account also for the ring current's contribution to equatorial $\Delta B$ as a function of $r$ and $\beta$. The self-consistent magnetic field should satisfy the force balance, as specified by the equation $\mu_0^{-1}\vec{B}\times(\nabla\times\vec{B}) = -\nabla\cdot\vec{\vec{P}}$, where $\mu_0$ is the permeability of free space and $\vec{\vec{P}}$ is the pressure tensor. Under these assumptions, the radial component of the left-hand side of the force balance equation in the equatorial plane can be expressed via an ordinary differential equation in $L$ or in $r_0$, where $r_0$ is the geocentric radial distance in the equatorial plane (related to $L$ by solving the equation of a field line for $r$ at $\theta = \pi/2$). Given the plasma pressure in the equatorial plane as a function of $r_0$ and $\beta$ as a result of our bounce-averaged guiding center simulations of representative ring-current particles, we can thus solve the pressure-balance equation to obtain the self-consistent magnetic field. We will show some preliminary results found by applying this method to actual ring-current plasma simulations.
DE: 2730 Magnetosphere--inner
DE: 2740 Magnetospheric configuration and dynamics
DE: 2753 Numerical modeling
DE: 2778 Ring current
DE: 2788 Storms and substorms
SC: SPA-Magnetospheric Physics [SM]
MN: 2004 AGU Fall Meeting