Investigation of Field-aligned Currents Onboard of Interkosmos Bulgaria-1300 Satellite

1. Investigation of Field-aligned Currents Onboard of Interkosmos Bulgaria-1300 Satellite Dimitar Danov Solar-Terrestrial Influences Laboratory, Bulgarian academy of sciences

2. Introduction BULGARIA-1300 SATELLITE

3. Introduction BULGARIA-1300 SATELLITE Triaxial Fluxgate Magnetometer

4. Introduction Reason for … The statistical picture of large-scale field-aligned current (FAC) distribution has been obtained by Iijima and Potemra [1976] and by many satellite observations including latest Iridium results (Waters et al. [2001]). At the same time it is well known that large-scale field aligned current sheets are frequently stratified on multiple sheets (see Kamide and Rostoker [1977]). Antonova et al. [1998] discussed the mechanisms of such stratification.

5. Introduction Reason for … Investigations of the stratifications of large-scale field-aligned current bands on multiple field-aligned current sheets were limited by the assumption of the existence of infinite homogeneous current sheets. However, this assumption is valid when the width of the current sheet in longitudinal direction exceeds significantly its size along the meridian. This allowed reducing the three-dimensional problem to one-dimensional one and deriving FAC’s structure from one-satellite measurements. We try to evaluate the corrections considering a finite size of the sheet in the longitudinal direction and show the efficiency of our methodology using Interkosmos–Bulgaria-1300 magnetic field measurements.

6. We identify current sheets in magnetic field data “manually” by the following algorithm: • Process the telemetry and obtain three magnetic components in satellite coordinate system • Subtract the model IGRF geomagnetic field from the measured field • Transform the magnetic vector to three components in spherical Solar Magnetic (SM) coordinate system • Analyse the graphs of these components to identify the sheets • The value (and sign) of the current density is evaluated by a simple program In the past we have done all these steps automatically, but we obtained many false currents.

7. Sources of errors ?? Whenever the current density is estimated from thecurl(B)the following approximation is very popular curl(B) ~ By/x ~ By/Vxt

8. MODELS SELECTED • We use two models of current sheets, to understand possible errors : • Model of a finite plane current sheet • current sheet in X,Y directions with a rectangular form and constant current density inside the rectangular: • j=const if a1 x a2, b1 y b2 and • j=0 outside this interval. • Conical model (used by Tsyganenko 2001) • current sheet is defined inside the spherical cone : • j=const if a1  a2 and • j=0 outside this cone.

9. Model of a finite current sheet The normal (meridional) component of the magnetic field disturbance for such configuration is symmetric with respect to the tangential surface of symmetry and antysymmetric with respect to the normal one. The tangential (zonal) component is antysymmetrical with respect to the tangential surface of symmetry andsymmetrical with respect to the normal one.

10. Model of a finite current sheet One sheet, s/c trajectories are normal to it As seen from graphs on this slide,By/xapproximation and the possible limited size of the current in Y-directioncan be a cause for the appearance of ‘false’ currents. Trace - 3 near the ‘termination’ of current sheet Trace - 1 near the ‘center’ of current sheet Changes in the tangential and normal components of the magnetic field corresponding to the field-aligned current sheet shown at the top, the satellite orbit is perpendicular to the sheet.

11. Model of a finite current sheet One sheet, s/c trajectories are normal to it  the current sheet  False current sheet with By/x >0 Another false current sheet By/x >0 Trace - 3 near the ‘termination’ of current sheet Trace - 1 near the ‘center’ of current sheet

12. Model of a finite current sheet Two sheets, trajectories are normal to it The value of the current density is equal in both sheets, but their directions are opposite. In case of TRACE 1 & 2 one can estimate equivalent (and opposite) current densities. NO false currents Trace - 2 Trace - 1 In case of TRACE 3 one can estimate different (and opposite) current densities. AS WELL false currents Trace - 3

13. Model of a finite current sheet Very complicated case Cross with big angle Cross with fair approach to normal Cross with fair approach to tangent It is good that B-1300 satellite has such Inclination, and cuts across FAC sheets near it’s normal.

14. Obtaining of ‘false’ currents sheets can be demonstrated also with the conikal model.

15. Conical model (used by Tsyganenko in his work from 2001) This model is defined in spherical coordinate system. It is complicated but self consistent. outside the interval a1    a2 and where in the interval a1    a2

16. Conical model (used by Tsiganenko 2001) Formulas for meridional and zonal components of the magnetic field are to long, for that reason we will show only graphs.  the current sheet  zonal False current sheet with By/x <0 meridional False current sheet with By/x <0 Trajectory provisionally parallel to 5 MLT

17. In the model of finite current sheet we evaluate the difference between the value of the modeled tangential component of the magnetic field at both edges of the sheet. The X-axis represents the distance (in percents) from satellite trajectory tosymmetry surface (100%) • The calculated current density is always less than the real one; the By/xapproximation works well only in case when s/c crosses the sheet near its middle; • There are cases when the satellite does not cross the sheet, but the change in the tangential component is significant – one can obtain false current.

18. Magnetic data & Field Aligned Currents The procedure we use to avoid obtaining false currents: We calculate the currents only for those measurements which: 1. The graph of the longitudinal component Bis linear; 2. DB between both edges of the linear interval is greater than 50nT; 3.  = max(B) – min(B) < 30nT (for the meridional component); 4. The correlation coefficient between B and  is greater 0.9 5. j = 0 -1B/  - coefficient of the linear regression B to 

19. Magnetic data & Field Aligned Currents Tsyganenko-2001 model field is plotted with symbols. Tsy-2001 c/ sheet  Tsy-2001 c/ sheet  Seance 217. Midnight satellite trajectory. Series of FACs are present.

20. Magnetic data & Field Aligned Currents Tsyganenko-2001 model field is plotted with symbols. Tsy-2001 c/ sheet  Tsy-2001 c/ sheet  Seance 183. Midnight satellite trajectory. FAC is directed into ionosphere.

21. Statistics of current sheets FAC density thickness distribution measured aboard ICB-1300 satellite It is possible to see that more than 60% of observed current sheets have thickness less than 1.50.

22. Conclusions and discussion • The conducted analysis shows that both components (meridional and zonal) of the magnetic field must be taken into account when estimating the FAC density and position, to eliminate the false groups of sheets. • Estimated from measurements FAC density is always smaller than the real one, up to 25% at the end of the sheet. • The FAC density is greater and thickness is smaller than usually accepted. We compared the results of our analysis with Tsyganenko [2002] and Papitashvili et al. [2001] models and results of Lukianova et al. [2001]. It is possible to mention that Intercosmos-Bulgaria-1300 observations show that FAC density is greater and thickness is usually smaller, than in these works. Unfortunately we have no sufficient data to investigate the dependence of FAC's density or thickness on other parameters like Kp, Dst or Solar Wind parameters

23. THANK YOU