2017-10-16

Channel Energy Ranges in the LANL GPS Charged Particle Dataset

The documentation for the CXD instrument mentions eleven electron channels and five proton channels. Some details are given in a referenced paper. Unfortunately, that paper is reproduced in black-and-white, and the response curves contained therein require colour in order to be interpreted correctly; in addition, a lot of details are still absent from the paper.

So here is the information I have been able to gather so far regarding the channels of the CXD instrument. Thank to John Sullivan of LANL for this information.

I note that there are minor inconsistencies in the publicly available documentation. These inconsistencies, however, would seem likely to be less than the accuracy and reproducibility limits for the instrumentation. What follows is my best interpretation of the information I have been able to gather.

This paper and this paper provide nominal energy ranges for the 11 electron channels. Adding the equivalent values of γ to that information gives us this table:

CXD Electron Channel Energy Ranges (MeV)
Channel Detector Min Energy Min γ Max Energy Max γ
E1 LEP 0.14 1.27 0.23 1.45
E2 LEP 0.23 1.45 0.41 1.80
E3 LEP 0.41 1.80 0.77 2.51
E4 LEP 0.77 2.51 1.25 3.45
E5 LEP 1.26 3.46 68 134
E6 HEP 1.33.54 1.7 4.33
E7 HEP 1.7 4.33 2.2 5.30
E8 HEP 2.2 5.30 3.0 6.87
E9 HEP 3.0 6.87 4.1 9.02
E10 HEP 4.1 9.02 5.8 12.35
E11 HEP 5.8 12.35

There are two distinct detectors on the instrument, the low-energy detector (often denoted "LEP") and the high-energy detector ("HEP"). As shown above, electron channels 1 to 5 are from the LEP, channels 6 to 11 are from the HEP.

The LEP and HEP detectors respond to both electrons and protons (and photons, which we will ignore for now). For protons the LEP detector has two channels, with threshold energies of "about" 6 MeV and 10 MeV, with an upper limit of 70 MeV. This paper gives the complete ranges for the proton channels as:

CXD Proton Channel Energy Ranges (MeV)
Channel Detector Minimum Energy Min γ Maximum Energy Max γ
P1 LEP 6 1.01 10 1.01
P2 LEP 10 1.01 50 1.05
P3 HEP 16 1.02 128 1.14
P4 HEP 57 1.06 75 1.08
P5 HEP 75 1.08

I am informed that the lower limits in these tables are reasonably accurate; however, the upper limits are rather soft, as the channels typically have some response to particles of higher energies.

The detailed transfer function between actual particle energy and the measured flux values is currently available only for the eleven electron channels (i.e., not the proton channels), and only for the satellites carrying the CXD instrument. There are two sets of transfer functions, one for SVN 53 through 61, and one for SVN 62 through 73. The two sets of numerical coefficients that define the transfer functions are available in a spreadsheet file in OpenDocument format here (the  coefficients for SVN 53 to 61 are on the first sheet, the coefficients for the remaining satellites on the second).

Plotting the transfer functions graphically, as below, gives us a better feel for the responses in the various channels.



 







 In the same way, we can plot the response curves for the remaining satellites:












These curves are a far cry from the ideal curves of a perfect instrument: in particular, note that all channels, even the low-energy ones, have a larger response to high energy electrons than to low energy ones: the chief difference between the high-energy channels and the low-energy ones being that the high-energy channels effectively suppress any response to low-energy electrons.

Thus, a notional stream of high-energy particles would be detected by all channels. In a more realistic stream with a low-energy component, the high-energy component might well swamp the contribution from the low-energy particles, even in the (notionally) low-energy channels; but the high-energy components could be determined from the high-energy channels (which are effectively immune to contamination from low-energy electrons), and then removed prior to the analysis of the low-energy channels.

Obtaining actual flux values from the instrument is therefore not a trivial task, and will be examined in more detail in a subsequent post.



2017-10-09

Long-Term Charged Particle Rates from LANL GPS data

Prior posts in this series:
This post contains more substantial plots from the stage 35 dataset of charged particle detections from the CXD instrument aboard the GPS satellites. As in the earlier posts, the event count data are used as-is, with no attempt to impose corrections or calibrations. (Neither do I present herein any interpretations of the several interesting features that are apparent in the figures below.)

A description of the style of each plot follows the first one.


Several items of interest are plotted in a single figure. Each year is divided into twelve segments of equal duration, and data are averaged and plotted for each such bin, with time increasing in the conventional manner from left to right. The date marked on the x axis corresponds to the start of the given year.

  1. The upper panel shows data from the eleven electron channels, E1 being the lowest-energy channel and E11 the highest-energy channel. The key to the colours used is at the extreme right hand of the figure, with the logarithmic scale running from 0.1 event per second to 10,000 events per second.
  2. The lower panel shows data from the five proton channels, P1 being the lowest-energy channel and P5 the highest-energy channel. The key to the colours used is to the left of the rainbow gradient at the right of the figure, with the logarithmic scale running from 0.1 event per second to 1,000 events per second.
  3. The red wiggly line on the upper panel is the mean total event rate for electrons (i.e., over all electron channels, from E1 to E11). The scale, in events per second, is provided by the red numbers in the upper part of the grey margin that surrounds the main plotting region of the figure. The logarithmic scale runs from 10 events per second (at the black line dividing electrons from protons) to 100,000 events per second (at the top of the plotting region).
  4. The red wiggly line on the lower panel is the mean total event rate for protons (i.e., over all proton channels, from P1 to P5). The scale, in events per second, is provided by the red numbers in the lower part of the grey margin that surrounds the main plotting region of the figure. The logarithmic scale runs from 1 event per second (at the bottom of the plotting region) to 1,000 events per second (at the black line dividing electrons from protons).
  5. The purple wiggly lines on the lower and upper panels show the 10.7 cm solar flux as recorded at the Dominion Radio Astrophysical Observatory station at Penticton, British Columbia. The same data are presented in each panel: only the vertical scale differs. For each panel, the scale is given in purple at the left edge of the plotting region, and runs linearly from a flux of 50 units (which is purely notional, as the flux index never reaches such a low value), to a flux of 250 units.
Each of the following figures is constructed similarly, one figure for each satellite carrying a CXD instrument.




















2017-10-02

Improved Plots of Charged Particle Events on GPS Satellites, from LANL Public Data

This is a successor to this earlier post.

The following plots show the same information as the prior post, but with the electron and proton colour gradients scaled independently for each type of particle and each spacecraft, as well as a larger number of plottable colours, so as to provide more visible detail.

As a reminder, the plotted data are raw, taken directly from the stage 35 public dataset, with no attempt to make any corrections or calibrations.






















The electron data appear much more consistent from satellite to satellite on this time scale than the proton data. Also, it is clear that some channels are more or less "sensitive" than their neighbouring channels, insofar as the raw counts reported in the dataset are concerned. Presumably, once proper calibration has been performed, the latter phenomenon will no longer be apparent.