Ionospheric Data
The ionosphere is that part of the upper atmosphere where free electrons occur in sufficient density to have an appreciable influence on the propagation of radio frequency electromagnetic waves. This ionization depends primarily on the Sun and its activity. Ionospheric structures and peak densities in the ionosphere vary greatly with time (sunspot cycle, seasonally, and diurnally), with geographical location (polar, auroral zones, mid-latitudes, and equatorial regions), and with certain solar-related ionospheric disturbances.
The major part of the ionization is produced by solar X-ray andultraviolet radiation and by corpuscular radiation from the Sun. The most noticeable effect is seen as the Earth rotates with respect to the Sun; ionization increases in the sunlit atmosphere and decreases on the shadowed side. Although the Sun is the largest contributor toward the ionization, cosmic rays make a small contribution. Any atmospheric disturbance effects the distribution of the ionization.
The ionosphere is a dynamic system controlled by many parameters including acoustic motions of the atmosphere, electromagnetic emissions, and variations in the geomagnetic field. Because of its extreme sensitivity to atmospheric changes, the ionosphere is a very sensitive monitor of atmospheric events.
In some circles it is thought that there is persuasive evidence of an ionospheric precursor to large earthquakes that can be used a predictor. Besides the obvious acoustic waves generated before and during an earthquake, a part of the preparation process of large earthquakes is the generation of electromagnetic emissions (EMEs). These EMEs have been detected in the ionosphere up to six days prior to a large earthquake, such as with the May 1960, Chilean 8.3 earthquake.
The most accurate way of measuring the ionosphere is with a ground-based ionosonde, which records data as ionograms
Ionograms are recorded tracings of reflected high frequency radio pulses generated by an ionosonde. Unique relationships exist between the sounding frequency and the ionization densities which can reflect it. As the sounder sweeps from lower to higher frequencies, the signal rises above the noise of commercial radio sources and records the return signal reflected from the different layers of the ionosphere. These echoes form characteristic patterns of "traces" that comprise the ionogram. Radio pulses travel more slowly within the ionosphere than in free space, therefore, the apparent or "virtual" height is recorded instead of a true height. For frequencies approaching the level of maximum plasma frequency in a layer, the virtual height tends to infinity, because the pulse must travel a finite distance at effectively zero speed. The frequencies at which this occurs are called the critical frequencies. Characteristic values of virtual heights (designated as h'E, h'F, and h'F2, etc.) and critical frequencies (designated as foE, foF1, and foF2, etc.) of each layer are scaled, manually or by computer, from these ionograms. Typically, an ionosonde station obtains one ionogram recording every 15 minutes. When the scaling is done manually only the hourly recordings are routinely reduced to numerical data. Modern ionosondes with computer-driven automatic scaling procedures routinely scale all the ionograms recorded. The resulting numerical values, along with the original ionograms and station reports, are archived at five World Data Centers (WDCs) for Ionosphere.
DEFINITIONS OF THE IONOSPHERIC REGIONS
The ionosphere is divided into four broad regions called D,E, F, and topside. These regions may be further divided into several regularly occurring layers, such as F1 or F2.D Region: The region between about 75 and 95km above the Earth in which the relatively weak) ionization is mainly responsible for absorption of high-frequency radio waves. E Region: The region between about 95 and 150km above the Earth that marks the height of the regular daytime E layer. Other subdivisions isolating separate layers of irregular occurrence within this region are also labeled with an E prefix, such as the thick layer, E2, and a highly variable thin layer, Sporadic E. Ions in this region are mainly O2+. F Region: The region above about 150km in which the important reflec-ting layer, F2, is found. Other layers within this region are also described using the prefix F, such as a temperate-latitude regular stratification, F1, and a low-latitude, semi-regular stratification, F1.5. Ions in the lower part of the F layer are mainly NO+ and are predominantly O+ in the upper part. The F layer is the region of primary interest for radio communications.
Topside: This part of the ionosphere starts at the height of the maximum density of the F2 layer of the ionosphere and extends upward with decreasing density to a transition height where O+ ions become less numerous than H+ and He+ ions. The transition height varies but seldom drops below 500km at night or 800km in the daytime, although it may lie above 1000km. Above the transition height, the weak ionization has little influence on transionospheric radio signals.
IIWG Data Exchange Format
The International Union of Radio Science (URSI) Ionospheric Informatic Working Group (IIWG) has produced a data exchange format to address problems such as scaled time resolutions greater than the traditional hour-based usage found in previously accepted formats. For a full description of the data format and the philosophy behind its creation see "Ionogram Characteristics at Uneven Data Rates," by Gamache and Reinisch [1989].
The IDD software can access data only in the IIWG format. This format allows data from all sampling rates to be included in the database. Even though most of the data on these CDs contain only hourly values, data taken on autoscaling equipment produce fully scaled data for every ionogram. Most future data will come from equipment which automatically scales the ionograms. The IIWG format allows the archiving of all data regardless of the resolution.
IIWG Format Description
|
Data Group 1: Station headers (informative and encoding data)
|
|
FORTRAN
RECORD # |
FORMAT |
DESCRIPTION |
|
1 |
A30 |
Station name |
|
1 |
A5 |
Station code |
|
1 |
I4 |
Meridian time used by station (in degrees) |
|
1 |
F5.1 |
Latitude +North, -South |
|
1 |
F5.1 |
Longitude, East |
|
1 |
A10 |
Scaling type: Manual/Automatic |
|
1 |
A10 |
Data editing: Edited/Non-edited/Mixed |
|
1 |
A30 |
Ionosonde system name
|
|
Data Group 2: Data headers (measurement times for the month)
|
|
FORTRAN
RECORD # |
FORMAT |
DESCRIPTION |
|
2,3 |
30I4* |
Year |
|
2,3 |
30I4* |
Month |
|
2,3 |
30I4* |
Number of days in the month, M |
|
2,3 |
30I4* |
Number of Parameters |
|
2,3 |
30I4* |
Numbers of measurements total |
|
2,3 |
30I4* |
Numbers of measurements for each of the M days, NM |
|
4,i |
12A10* |
List of parameters |
|
i+1,j |
12A10* |
Dimensions |
|
j+1,k |
60A2* |
List of corresponding URSI codes |
|
Data Group 3: The data (actual values of the parameters and the corresponding hourly medians and statistics)
|
|
FORTRAN
RECORD # |
FORMAT |
DESCRIPTION |
|
k+1,l |
20(3I2)* |
The NM sample times Hh:Mm:Ss (repeated for each of the M days) |
|
l+1,m |
24(I3,2A1)* |
The N1 values of characteristic 1 for day 1 (repeated for each of the M days) |
|
m+1 |
24(I3,2A1) |
Hourly Medians for characteristic 1 |
|
m+2 |
24(I2,2A1) |
The counts for the hourly medians, Range |
|
m+3 |
24(I3,2A1) |
Upper quartile |
|
m+4 |
24(I3,2A1) |
Lower quartile |
|
m+5 |
24(I3,2A1) |
Upper decile |
|
m+6 |
24(I3,2A1) |
Lower decile |
|
m+7,n |
24(I3,2A1)* |
The N2 values of characteristic 2 for day 1 (and so on, repeated for each characteristic) |
|
* Format is repeated as many times as necessary to read/write the data.
|
URSI Qualifying and Descriptive Letters
Parameter Coding by International Union of Radio Science (URSI)
Five positions are reserved for scaled values: 3 numeric for the value, 1 alpha for qualifier, and 1 alpha for descriptor. Valid combinations are: a full weight value (123__), an unqualified, but described value (123_A), a qualified and described value (123UA), or a replaced value(____A).
QUALIFICATION Letters
A less than (used only fbEs in case of total blanketing).
D greater than.
E less than.
I interpolated.
J deduced from x component.
M mode uncertain.
T smoothed from sequence.
U uncertain.
Z deduced from z component.
DESCRIPTIVE Letters
A blanketing.
B absorption.
C non-ionospheric (equipment).
D above upper frequency range of equipment.
E below lower frequency range of equipment.
F frequency spread.
G ionization density too small.
H stratification.
K particle E layer.
L no sufficiently defined cusp between F layers.
M mode uncertain.
N superimposed layers.
O measurement refers to o component.
P man-made perturbation of parameters.
Q range spread.
R attenuation near critical frequency.
S interference.
T interpolated.
V forked trace.
W above height range.
X measurement refers to x component.
Y lacuna (tilt).
Z measurement refers to z component.
When qualifier and descriptive letters are not used with the data, it must be clear that their absence does not indicate full weight values. The IIWG has developed a code to designate this condition properly, as indicated in the table below:
|
Scaling Method
|
Editing |
QD entry |
|
Manual |
Edited
(QD entered)
---------------------------
Edited
(no QD entered) |
?QD?
---------------------------
?/ ?
|
|
Automatic |
Edited
(no QD entered)
---------------------------
Non-edited
(no QD entered) |
?/ ?
---------------------------
?// ? |
Note QD = Qualification and Descriptive Letters
For more details on qualification and description letters, see W.R. Piggott and K. Rawer, editors, URSI Handbook of Ionogram Interpretation and Reduction, 2nd edition (Report UAG-23, World Data Center A for Solar-Terrestrial Physics, Boulder, CO, 1972).
F Layer Parameters
|
Parameter |
URSI code # |
Dimension |
Description
|
|
foF2 |
00 |
.1 MHz |
F2 layer o-mode (ordinary) critical frequency |
|
fxF2 |
01 |
.1 MHz |
F2 layer x-mode (extraordinary) critical frequency |
|
fzF2 |
02 |
.1 MHz |
F2 layer z-mode critical frequency |
|
M3000F2 |
03 |
.01 MHz |
F2 layer M factor (the ratio of the maximum usable frequency divided by the critical frequency) |
|
h'F2 |
04 |
km |
F2 layer o-mode minimum virtual height |
|
hpF2 |
05 |
km |
An estimate of the true height of the F2 layer (measurement of the ordinary mode virtual height at a frequency of 83.4% of foF2) |
|
h'Ox |
06 |
km |
F layer minimum virtual height of the x-mode trace at a frequency equal to the foF2 |
|
MUF3000F2 |
07 |
.1 MHz |
F2 layer maximum usable frequency for a 3000km path |
|
hc |
08 |
km |
The height of the maximum obtained by fitting a theoretical h'F curve for the parabola of best fit to the observed ordinary mode trace near foF2 and correcting for under-lying ionization |
|
qc |
09 |
km |
F layer scale height |
|
foF1 |
10 |
.01 MHz |
F1 layer o-mode critical frequency |
|
fxF1 |
11 |
.01 MHz |
F1 layer x-mode critical frequency |
|
M3000F1 |
13 |
.01 MHz |
F1 layer M factor (see URSI code 03) |
|
h'F1 |
14 |
km |
F1 layer o-mode minimum virtual height |
|
h'F |
16 |
km |
F layer o-mode minimum virtual height |
|
MUF3000F1 |
17 |
.1 MHz |
F1 layer maximum usable frequency (see URSI code 07) |
E Layer Parameters
|
Parameter |
URSI code # |
Dimension |
Description |
|
foE |
20 |
.01 MHz |
E layer o-mode critical frequency |
|
foE2 |
22 |
.01 MHz |
E2 layer o-mode critical frequency (when it occurs it is between the normal E and F1 layers) |
|
h'E |
24 |
km |
E layer o-mode minimum virtual height |
|
h'E2 |
26 |
km |
E2 layer o-mode minimum virtual height |
Es Layer Parameters
|
Parameter |
URSI code # |
Dimension |
Description |
|
foEs |
30 |
.1 MHz |
Es layer highest o-mode frequency at which a mainly continuous Es trace is observed |
|
fxE |
31 |
.1 MHz |
Es layer highest x-mode frequency at which a mainly continuous Es trace is observed |
|
fbEs |
32 |
.1 MHz |
The blanketing frequency of layer used to derive foEs |
|
ftEs |
33 |
.1 MHz |
Top frequency of the Es trace (any mode) |
|
h'Es |
34 |
km |
The minimum virtual height of the layer used to derive foEs |
|
Type Es |
36 |
|
A characterization of the shape of the Es trace |
Other Parameters
|
Parameter |
URSI code # |
Dimension |
Description |
|
foF1.5 |
40 |
.01 MHz |
The o-mode critical frequency of the F1.5 intermediate stratification (between F1 and F2) |
|
fmin |
42 |
.1 MHz |
The lowest frequency at which an o-mode echo is observed on the ionogram |
|
M3000 F1.5 |
43 |
.01 MHz |
F1.5 layer M factor (see URSI code 03) |
|
h'F1.5 |
44 |
km |
F1.5 layer o-mode minimum virtual height |
|
fm2 |
47 |
.1 MHz |
The fmin for the second order o-mode trace |
|
hm |
48 |
km |
The height of the maximum electron density of the F2 layer calculated by the Titheridge method |
|
fm3 |
49 |
.1 MHz |
The fmin for the third order o-mode trace |
Spread E/Oblique Parameters
|
Parameter |
URSI code # |
Dimension |
Description |
|
foI |
50 |
.1 MHz |
The highest o-mode frequency of spread F |
|
fxI |
51 |
.1 MHz |
The highest frequency of spread F traces (any mode) |
|
fmI |
52 |
.1 MHz |
The lowest o-mode frequency at which spread traces are observed for the F layer |
|
M3000I |
53 |
.01 MHz |
M Factor deduced from upper frequency edge of spread traces and fxI (see URSI code 03) |
|
h'I |
54 |
km |
Minimum slant range of the spread F trace |
|
dfs |
57 |
.1 MHz |
Frequency range of the spread |
N(h) Parameters
|
Parameter |
URSI code # |
Dimension |
Description |
|
fh'F2 |
60 |
.1 MHz |
The frequency at which h'F2 is measured |
|
fh'F |
61 |
.1 MHz |
The frequency at which h'F is measured |
|
h'mF1 |
63 |
km |
The maximum virtual height in the o-mode F1 cusp |
|
h1 |
64 |
km |
True height at f1 Titheridge method |
|
h2 |
65 |
km |
True height at f2 Titheridge method |
|
h3 |
66 |
km |
True height at f3 Titheridge method |
|
h4 |
67 |
km |
True height at f4 Titheridge method |
|
h5 |
68 |
km |
True height at f5 Titheridge method |
|
H |
69 |
km |
Effective scale height at hmF2 Titheridge method |
REFERENCE: Definition of Characteristic extracted from UAG23 (URSI Handbook of Ionogram Interpretation and Reduction, November, 1972)
|