Impact of Realistic Detector Response and Azimuthal Drift on the Quantification of Energetic Electron Precipitation using POES/MetOp

Abstract

Drift-Diffusion model outputs and its derivative products used in this study are included here as a part of data avaliability for the publication to appear in Geophysical Research Letters. For this study, analaysis was performed in Python, and the data output are in formats of .txt and .pickle, which can be easily parsed into Python. The ".pickle" extension is a Python specific format, and can be loaded using the Pickle Python library. Please refer to the newest version of this repository. Contained in this repository is a compressed folder, Fig1&2.zip, and 5 seperate data files. Descriptions for each of these files can their usage is described below: Fig1&2.zip: This compressed folder contains the necessary data to reproduce Figure 1b1, Figure 1b2, Figure 2b1, and Figure 2b2. Inside, you will find "a0bavg.txt", "energy_keV.txt", and "phi.txt" - these are the equatorial pitch angle, energy (in keV), and longtitudinal grids for the model distribution function. The equatorial distribution function are contained in 4 .pickle files, corresponding to: f_Dw1e-3_withDrift.pickle: Distribution function simulated with Dw = 1e-3 per second with azimuthal drift f_Dw1e-3_noDrift.pickle: Distribution function simulated with Dw = 1e-3 per second with no azimuthal drift f_Dw1e-5_withDrift.pickle: Distribution function simulated with Dw = 1e-5 per second with azimuthal drift f_Dw1e-7_withDrift.pickle: Distribution function simulated with Dw = 1e-7 per second with azimuthal drift f_Dw1e-7_noDrift.pickle: Distribution function simulated with Dw = 1e-7 per second with no azimuthal drift The indexing order for each distribution function is [pitch angle, phi ,energy]. The distribution at each energy is normalized, and requires energy scaling parameters (e.g., a known energy spectra). response_contribution_example.pickle: Specifically contains the data shown in Figure 1c1 - Figure 1c4. Distribution functions at the equator are mapped to POES/MetOp's location, and is stored inside structure as "f_highdiff" and "f_lowdiff" for Figure 1c1 and Figure 1c2 respectively. The computed E3 detector response at various incident angle are stored inside the structure as "E3_90_lowdiff" for 90 degree telescope and under low diffusion levels. The high diffusion and 0 degrees telescope response cases are named similarly. Lastly, the local pitch angle grid and energy grids are stored as "aL" and "Ek" respectively. ratio_Dw_longitude_noDrift_fullresponse.pickle & ratio_Dw_longitude_noDrift_inFOV.pickle & ratio_Dw_longitude_withDrift_fullresponse.pickle: Stored within these structures are the the simulated count rate ratio parameterized by Dw and longitude for the North hemisphere, "CR_Ratio_North", and the South hemisphere, "CR_Ratio_South". The longitude grids are stored under "phi_grid", and Dw grids as "Dw_grid". These count rate ratio products are directly used to produce Figure 3b, Figure 4a1, and Figure 4a2. ionS_lon126_Dw1p00e-5.txt & ionS_lon126_Dw5p62e-5.txt: These text files contain the data required to produce Figure 3c - computed with our model distribution function using the G4EPP method (see paper). There are two roles in these files: row 1 is altitude in kilometers, and row 2 is the ionization rate. Fig3d.txt: Text file containing the computed data used to produce Figure 3d. First colomn is E3 CR90/CR0 ratio, and second column is the ratio of the different ionization peaks computed between more and less accurate methods.