ibpmodel 1.3.0

Ionospheric Bubble Probability

The ionospheric bubble probability statistical model is a Swarm L2 product, named IBP_CLI. The output of the Ionospheric Bubble Probability (IBP) product is an index, that depends on the day of year or the month of the year, geographic longitude, local time and solar flux index.

The output floating point index ranges 0-1 and characterizes the percentage probability of low latitude bubble occurence at the specified time, location and solar flux.

This empirical IBP model has been derived from magnetic observations obtained by the CHAMP and Swarm missions. The model is representative for the altitude range 350 - 500 km and low geographic latitudes of +/- 45 degree.

Documentation

Detailed documentation can be found at: https://ibp-model.readthedocs.io

Quick Start

Installation

Using pip:

$ pip install ibpmodel

Dependencies:

• numpy

• pandas

• matplotlib

• scipy

• cdflib

Usage

The return value of the function ibpmodel.calculateIBPindex() is of type pandas.DataFrame.

>>> import ibpmodel as ibp
>>> ibp.calculateIBPindex(day_month=15,           # Day of Year or Month
              longitude=0,                        # Longitude in degree
              local_time=20.9,                    # Local time in hours
              f107=150)                           # F10.7 cm Solar Flux index
   Doy  Month  Lon    LT  F10.7     IBP
0   15      1    0  20.9    150  0.3433

 >>> ibp.calculateIBPindex(day_month=['Jan','Feb','Mar'], local_time=22)
      Doy  Month  Lon  LT  F10.7     IBP
 0     15      1 -180  22    150  0.0655
 1     15      1 -175  22    150  0.0667
 2     15      1 -170  22    150  0.0685
 3     15      1 -165  22    150  0.0710
 4     15      1 -160  22    150  0.0754
 ..   ...    ...  ...  ..    ...     ...
 211   74      3  155  22    150  0.2146
 212   74      3  160  22    150  0.2123
 213   74      3  165  22    150  0.2117
 214   74      3  170  22    150  0.2119
 215   74      3  175  22    150  0.2123

[216 rows x 6 columns]

>>> ibp.plotIBPindex(doy=349)
>>>

The IBP model reproduces the high occurrence probability of EPDs ranging between 50-90% over the South American (75-25°W) sector and low occurrence probability over the Pacific sector during the period around December solstice.

>>> ibp.plotButterflyData(f107=150)
>>>

The monthly global occurrence rate of EPDs from the IBP model, is derived for a fixed value of F10.7=150 s.f.u for all integer longitudes at a resolution of 5° at the middle of each month and averaged between 19 and 1 LT. The seasonal and longitudinal variations of the EPD occurrence rates are particularly well-characterized by the IBP model as compared to its climatology with highest rates seen around the equinoxes and winter solstice in the America-Atlantic-Africa region and lowest rates during November-February in the Pacific sector and during May-July in the America-Atlantic and Indian sectors.

References

Stolle, C., Siddiqui, T. A., Schreiter, L., Das, S. K., Rusch, I., Rother, M., & Doornbos, E. (2024). An empirical model of the occurrence rate of low latitude post‐sunset plasma irregularities derived from CHAMP and Swarm magnetic observations. Space Weather, 22, e2023SW003809. https://doi.org/10.1029/2023SW003809

Lucas Schreiter, Anwendungsorientierte Modellierung der Auftretenswahrscheinlichkeit und relativen Häufigkeit von äquatorialen Plasmabubbles, Master’s thesis, Institute of Mathematics, University of Potsdam, 2016. (in German only.)

Information for developers

Setup environment

$ git clone https://igit.iap-kborn.de/ibp/ibp-model.git
$ cd ibp-model
$ pip install -r requirements-dev.txt
$ pip install -e .

Test of package using doctest

$ python src/ibpmodel/ibpcalc.py

No error should occur.

Test run of the documentation

$ cd docs
$ make clean && make html

The docs/build/html/ directory contains the html files. Open index.html in browser. The results of the code examples on the usage page are generated automatically. Therefore the ibpmodel package must be installed (pip install -e .).