Plasmapause location

The lpp provides functionalities for calculating the plasmapause location.

Models:

  • Carpenter and Anderson (1992)

  • O’Brien and Moldwin (2003)

  • Moldwin et al., (2002)

CA1992(time, kp)

Calculates plasmapause location Lpp following Carpenter and Anderson [1992], eq (7).

M2002(time, kp)

Calculates plasmapause location (lpp) following Moldwin et al. [2002] eq.

OBM2003(time, index, index_type)

Calculates plasmapause location (lpp) following O'Brien and Moldwin [2003] eq.

Functions

rbamlib.models.lpp.CA1992(time, kp)

Calculates plasmapause location Lpp following Carpenter and Anderson [1992], eq (7).

Parameters:

  • time (ndarray) – Time in days as a numpy array.

  • kp (ndarray) – Kp-index, vector of geomagnetic activity indices corresponding to the times.

Returns:

lpp – Plasmapause location Lpp as a vector.

Return type:

ndarray

Notes

The plasmapause location Lpp is calculated using the formula:

\[L_{pp} = 5.6 - 0.46 \cdot Kp_{24}\]

where \(Kp_{24}\) is the maximum Kp-index value in the last 24 hours.

rbamlib.models.lpp.OBM2003(time, index, index_type)

Calculates plasmapause location (lpp) following O’Brien and Moldwin [2003] eq. (1).

Parameters:

  • time (ndarray) – Time in days as a numpy array.

  • index (ndarray) – Vector of the index corresponding to the type, same size as time.

  • index_type (str) – A string representing the type of index (‘Kp’, ‘Ae’, ‘Dst’).

Returns:

lpp – Vector plasmapause location, same size as time.

Return type:

ndarray

Notes

The plasmapause location Lpp is calculated based on the empirical model by O’Brien and Moldwin (2003).

\[L_{pp} = a \cdot Q + b\]

where ( Q ) is a representation of one of the indicators:

\[\begin{split}Q_{Kp} = \max_{-36,-2} Kp, \\ Q_{AE} = \log_{10} \left( \max_{-36,0} AE \right), \\ Q_{Dst} = \log_{10} \left| \min_{-24,0} Dst \right|\end{split}\]

When using Kp index, the first value will always be np.nan because maximum Kp is calculated between -36 and -2 hours prior current time

rbamlib.models.lpp.M2002(time, kp)

Calculates plasmapause location (lpp) following Moldwin et al. [2002] eq. (2).

Parameters:

  • time (ndarray) – Vector of time in days as a numpy array.

  • kp (ndarray) – Vector of Kp index, same size as time.

Returns:

lpp – Vector of plasmapause location in Re, same size as time.

Return type:

ndarray

Notes

The plasmapause location Lpp is calculated using the formula:

\[L_{pp} = 5.39 - 0.382 \cdot Kp_{12}\]

where \(Kp_{24}\) is the maximum Kp-index value in the last 12 hours.