Space Physics Seminar spring-2015

Accurate evaluation of the physical processes during the substorm growth phase, including formation of field-aligned currents (FACs), isotropization by current sheet scattering, instabilities, and ionosphere-magnetosphere connection relies on knowing the realistic 3 dimensional (3D) magnetic field configuration, which cannot be reliably provided by current available empirical models. We have first developed a 2D empirical plasma sheet pressure model using the Support Vector Regression Machine (SVRM) with observational data from THEMIS and Geotail. The model predicts the plasma sheet pressure accurately with median errors of 5%, and predicted pressure gradients agree reasonably well with observed gradients obtained from two-probe measurements. Then we established a 3D substorm growth phase magnetic field model, which is uniquely constructed from empirical plasma sheet pressures under the constraint of force balance. We investigated the evolution of model pressure and magnetic field responding to increasing energy loading, and their configurations under different solar wind dynamic pressure (PSW) and sunspot number. Our model reproduces the typical growth phase evolution signatures: plasma pressure increases, magnetic field lines become more stretched, westward perpendicular current is intensified and moves earthward, and the Region-2 FACs are enhanced. The model magnetic fields agree quantitatively well with observed fields. The magnetic field is substantially more stretched under higher PSW while the dependence on sunspot number is non-linear and less substantial. By applying our modeling to a substorm event, we found that (1) the equatorward movement of proton aurora during the growth phase is mainly due to continuous stretching of magnetic field lines, (2) the ballooning instability is more favorable during late growth phase around midnight tail where there is a localized plasma beta peak, and (3) the equatorial mapping of the breakup auroral arc is at X ~ –14 RE near midnight, coinciding with the location of the maximum growth rate for the ballooning instability.