OBSERVED

MODEL (M)}

ALTITUDE (KM)

OBSERVED

CURVE YR MO DY TIME SSNO
358 68 111
REFERENCE (R) 354 68 10 25 1200 110

(0)

1200 111

LAT LONG METH MAG REF

74.7 265.1
74.7 265.1

1 127
1 1az

FIG SPC
147
14.0

wi

N

~

7

ma

oe

it))

CURVE YR MO

382

REFERENCE (R) 0

MODEL (M)

ALTITUDE (KM)

69
Q

1
ie}

1

oY
3
0

T
5
ie)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

IME SSNO
10S
ti)

LAT LONG METH MAG REF

$8.8 265.8
0.0 60.0

FIG SPC

2 2138 27.5 8

o 00

0

|

V

2 3 4
LOG ELECTRON CONCENTRATION (N/CC)

SOLAR PROTO!
16

IN EVENT



OBSERVED (0)

CURVE YR
482 72

REFERENCE (R] O a

MODEL (M)

ALTITUDE (KM)

OBSERVED (0)

us

410

s

MO
5
it

D
1

a

Q

T
2
ie)

1

13

‘ME SSNO LAT
oO 0.0

LONG HETH HAG REF
73 ~-69.0 39.0

0.0

FI6 SPC
2 1100 1 it
o oo ()

g

yi

}

CURVE YR
480 72

REFERENCE (R} O 0

MODEL (M)

ALTITUDE (KM)

ite}

MO
S
ie)

oY
1
Q

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

TIME SSNO LAT

0

6 202 73

0 0.0

LONG METH MAG REF
-69.0 39.0

0.0

FIG SPC
2 1100 1 ou
o 00 it)

a 3 4
LOS ELECTRON CONCENTRATION (N/CC)

POLAR SUBSTORM
37



YY TIME SSNO LAT LONG METH MAG REF FIG SPC

OBSERVED (0) 225 66 2 275 09

REFERENCE (R) 0

MODEL (M)

ALTITUDE (KM)

OBSERVED

CURVE YR MO
3
Q

ia)

0

a 0.0

o
10 2326 «34 69.3 16.0
(9)

0.0

o 00

Qo

[|

(0)

CURVE
224

REFERENCE (R} O

MODEL (M)

ALTITUDE (KM)

YR
6
ie}

“Oo
3
Q

oY
9
fe)

T

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

IME SSNO LAT

2320 34 69.3

fa)

0 0.0

LONG METH MAG REF

16.0 2 275

0.0

o 00

FIG SPC

09
ie)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
AURORAL ABSORPTION EVENT

33



21

24
19

124
125

1108

66

129

71

APPENDIX A (Continued)

FIND CORRESPONDING REFERENCE PROFILE
DO 21 I=1,64

LI=I

IF (KURVE.EQ.IK(I)) GO TO 24

CONTINUE

Go TO 19

IKR=IR(II)

IKRL (L) =IKR

CONTINUE

IYL(L)=1Y

60.-FLON

IRFL (L) =IREF

WRITE OUT THE HEADER CARD INFORMATION

WRITE (8,100) KURVE,1Y,1M,1ID,LH,LM,ISSNO,FLAT,FLON,METH,
A MAG, KOL, IREF,IFIG,ISPEC,IKR

PRINT 109, KURVE,1Y,1M,ID,LH,LM,ISSNO, FLAT, FLON,METH,
A MAG, KOL, IREF,IFIG,ISPEC,IKR

GO TO 15

ALL DATA READ IN

REWIND 5

L
KK=9999
WRITE(8,100) KK

WRITE OUT HEADER CARDS FOR REFERENCE PROFILES
READ (5,100) KURVE,IY,IM,ID,LH,LM,ISSNO,FLAT,FLON, METH,
A MAG, KOL, IREF,IFIG,ISPEC

IF (KURVE.EQ.9999) GO TO 125

DO 124 J=1,LL

IF(KURVE.NE.IKRL(J)) GO TO 124

WRITE(11,100) KURVE,IY,IM,ID,LH,LM,ISSNO,FLAT,FLON, METH,
A MAG ,KOL, IREF,IFIG,ISPEC

PRINT 160, KURVE,IY,IM,ID,LH,LM,ISSNO,FLAT,FLON, METH,
A MAG, KOL,IREF,IFIG,ISPEC

GO TO 115

CONTINUE

GO TO 115

KK=9999

WRITE(11,108) KK

REWIND 5

FOR EACH REQUIRED PROFILE, CALCULATE STATISTICAL MODEL AND
FIND REFERENCE PROFILE

PRINT 210, LL

FORMAT (/1X,*THERE ARE *,13,* SUCH PROFILES*)
IF(LL.EQ.8) GO TO 5

READ IN THE DIGITIZED PROFILES ONE BY ONE

READ (6,118) KURVE

FORMAT (14)

IF (KURVE.EQ.9999) GO TO 74

Tl==5

J1=J1+6

J2=J1+5

READ (6,126) (FLN(J) ,H(J) ,J=J1,J2)

FORMAT (12F6, 2)

IF (H(J2).NE.@.) GO TO 69

SEE IF THIS CURVE IS A REQUIRED SPECIAL CASE

DO 71 L=1,LL

IF (KURVE.NE.NUM(L)) GO TO 71

CURVE IS REQUIRED .. WRITE IT OUT

WRITE (9,110) KURVE

WRITE(9,120) (FLN(J) ,H(J) ,J=1,32)

GO TO 40

CONTINUE

42

|
|
|



CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 48. 6112181200 49 45.4 2841 1 192 8s
REFERENCE (R) 49 6112181200 49 45.4 2841 1 192 8 0
MODEL (¥) |

2 i
s LZ |

ALTITUDE (KM)

-+_§41

a 3 4
LOG ELECTRON CONCENTRATION (N/CC!
CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC |

OBSERVED (0) 221 661 101200 28 37.9 284.5 2 190 15 I
REFERENCE (R) 220 65 12 15 1200 24 37.9 284.5 2 190 1 0 I
MODEL (M)

105

: 4 |
Ie

ALTITUDE (KM)

9 1 a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

WINTER ANOMALY DAY

{
|
I
|
|

10



UAG-48A

UAG~49

UAG-50

UAG-51

UAG-52

UAG-53

UAG-54

UAG-55

UAG-56

UAG~57

UAG-58

UAG-59

UAG-60

UAG-61

UAG-62

UAG-63

UAG-64

"synoptic Observations of the Solar Corona during Carrington Rotations 1580-1596 (11 October 1971 - 15 January 1973)",
(Reissue with quality images] by R. A. Howard, M. J. Koomen, D. J, Michels, R. Tousey, C. R. Detwiler, D. E.
Roberts, R. T. Seal and J. D. Whitney, E. 0. Hulbert Center for Space Research, NRL, Washington, D. C. 20375
and R. T. and S. F. Hansen, C. J. Garcia and E. Yasukawa, High Altitude Observatory, NCAR, Boulder, Colorado
80303, February 1976, 200 pages, price $4.27.

“Catalog of Standard Geomagnetic Variation Data", prepared by Environmental Data Service, NOAA, Boulder, Colorado,
August 1975, 125 pages, price $1.85.

“High-Latitude Supplement to the URSI Handbook on Ionogram Interpretation and Reduction", by W. R. Piggott, British
Antarctic Survey, c/o SRC, Appleton Laboratory, Ditton Park, Slough, England, October 1975, 292 pages. price $4.00.

“Synoptic Maps of Solar Coronal Hole Boundaries Derived from He II 304A Spectroheliograms from the Manned Skylab
Missions", by J. D. Bohlin and D. M. Rubenstein, E. 0. Hulbert Center for Space Research, Naval Research Lab-
oratory, Washington, D. C. 20375 U.S.A., November 1975, 30 pages, price 54 cents.

“Experimental Comprehensive Solar Flare Indices for Certain Flares, 1970-1974", compiled by Helen W. Dodson and
E. Ruth Hedeman, McMath-Hulbert Observatory, The University of Michigan, 895 Lake Angelus Road North, Pontiac,
Michigan 48055 U.S.A., November 1975, 27 pages, price 60 cents.

“Description and Catalog of Ionospheric F-Region Data, Jicamarca Radar Observatory (November 1966 - April 1969)", by
WoL. Clark and T. E. Van Zandt, Aeronomy Laboratory, NOAA, Boulder, Colorado 80302 and J. P. McClure, University
of Texas at Dallas, Dallas, Texas 75230, April 1976, 10 pages, price 33 cents.

“Catalog of Ionosphere Vertical Soundings Data", prepared by Environmental Data Service, NOAA, Boulder, Colorado
80302, April 1976, 130 pages, price $2.10.

“Equivalent Tonospheric Current Representations by a New Method, Illustrated for 8-9 November 1969 Magnetic Disturb-
ances", by Y. Kamide, Cooperative Institute for Research in Environmental Sciences, University of Colorado,
Boulder, Colorado 80302 and Geophysical Institute, University of Alaska, Fairbanks, Alaska 99701, H. W. Kroehl,
Data Studies Division, NOAA/EDS/NGSDC, Boulder, Colorado 80302, M. Kanamitsu, Advanced Study Program, National
Center for Atmospheric Research, Boulder, Colorado 80303, J. H. Allen, Data Studies Division, NOAA/EDS/NGSDC,
Boulder, Colorado 80302, and S.-1. Akasofu, Geophysical Institute, University of Alaska, Fairbanks, Alaska
99701, April 1976, 91 pages, price $1.60.

“Iso-intensity Contours of Ground Magnetic H Perturbations for the December 16-18, 1971 Geomagnetic Storm", by
X. Kamide, Cooperative Institute for Research in Environmental Sciences, University of Colorado, Boulder,
Colorado 80302 and Geophysical Institute, University of Alaska, Fairbanks, Alaska 99701 (currently Guest worker
at Data Studies Division, NOAA/EDS/NGSDC, Boulder, Colorado 80302), April 1976, 37 pages, price $1.39.

“Wanual on Ionospheric Absorption Measurements", edited by K. Rawer, Institut fiir Physikalische Weltraumforschung,
Freiburg, G.F.R., June 1976, 202 pages, price $4.27.

"ATS6 Radio Beacon Electron Content Measurements at Boulder, July 1974 - May 1975", by R. B. Fritz, Space Environment
Laboratory (currently with Wave Propagation Laboratory), NOAA, Boulder, Colorado 80302 USA, September 1976,
61 pages, price $1.04.

“Auroral Electrojet Magnetic Activity Indices AE(11) for 1974", by Joe Haskell Allen, Carl C. Abston and Leslie D.
Morris, National Geophysical and Solar-Terrestrial Data Center, Environmental Data Service, December 1976,
144 pages, price $2.16.

“Geomagnetic Data for January 1976 (AE(7) Indices and Stacked Magnetograms)" by J. H. Allen, C. A. Abston and
L. R. Morris, NGSDC/EDS/NOAA, July 1977, 57 pages, price $1.07.

“Collected Data Reports for STIP Interval IT 20 March - § May 1976", edited by Helen E. Coffey and John A.

Nckinnon, National Geophysical and Solar-Terrestrial Data Center, Environmental Data Service, August 1977,
313 pages, price $2.95.

“Geomagnetic Data for February 1976 (AE(7) Indices and Stacked Magnetograms)" by J. H. Allen, C. A. Abston and
L. R. Morris, NGSDC/EDS/NOAA, September 1977, 55 pages, price $1.11. "

“Geomagnetic Data for March 1976 (AE(7) Indices and Stacked Magnetograms)" by J. H. Allen, ¢. A. Abston and
L. R. Morris, NGSDC/EDS/NOAA, September 1977, 57 pages, price eth y an

"Geomagnetic Data for April 1976 (AE(8) Indices and Stacked Ma: inetograms)" by
LR. Morris, NGSDC/EDS/NOAA, February 1978, ignetograms)" by J. H. Allen, C. A. Abston and

USCOMM-NOAA-ASHEVILLE ,NC-10/78/550



15

APPENDIX B (Continued)

CONTINUE

RETURN

END

SUBROUTINE PICTUR(I)

DIMENSION X1(5) ,¥1(5) ,X2(9) ,Y2(9)

DATA X1/-.035,-.035,0. 5.035, .035/
DATA Y1/-.05,.85,-.95,.05,-.85/

DATA X2/-.035,-.035, .020,.035, 035,
X.620,-.05,0+,-.035/

DATA Y2/-.05,.05,-05,-04, D406,
XO. Dep 05/

GoTo (1,2), I

CALL PICTRC(X1,Y1,5)
RETURN

CALL PICTRC(X2,¥2,9)

RETURN
END

50



APPENDIX B

The program DREGPLT plots the D-region profiles as indicated in the diagrams. The version listed
below is suitable for a height range of 50 to 105 km. Slight changes are necessary for profiles con-
taining points outside this range.

au

15

Listing of DREGPLT

-THIS PROGRAM PLOTS ELECTRON DENSITY VERSUS HEIGHT
FOR D REGION PROFILES.

MAKCUS ERTLE PROGRAMMER, JUN. 1978.

PROGRAM DREGPLT (INPUT, OUTPUT, TAPE1=INPUT, TAPE2=OUTPUT,
XTAPE3)

COMMON A(2,6) ,X(199) ,¥(100) , INPT(1)

IEND=1

IEND=9

CALL TEKTRN (469)

CALLBGNPL(IEND)

KPT=¢
N

N=O

CALL TITLE(" $",-100,"LOG ELECTRON CONCENTRATION (N/CC) $",
X190,"ALTITUDE (KM)$",190,8.,6.)

CALL YINTAX

CALL XINTAX

CALL YAXANG (@.@)

CALL GRAF(@,1,6,58,5,105)

CALL GRAF (8,1,6,35,5,98)

CALL GRAF(6,1,6,55,5,115)

CALL GRID(1,1)

CALL CURVENO(IDIST, IYR, IMO, IDY,ITIME,ISUN,STLAT,

XSTLONG, IMETH, IMAG , ICOL, [REF ,UFIG,1SC,JREF,
XJDIST,JYR,JMO,JDY ,JTIME,JSUN, SUTLAT,

XSJTLONG ,JMETH, JMAG ,JCOL, JJREF ,UJFIG,JSC,JJJREF ,N,NH,MODE,I6,IC)
IF(N.EQ.2) 46,5

IF(NH.EQ.9) 20,6

CALL MARKER(MODE)

CALL EPOINTS(N)

CALL FILLUP(M,KPT,N)

IF(X(M) .EQ.9.08) 15,4

CALL CURVE(X,Y,KPT,1)

KPT=9

M

Ig=6

IF(IC.NE.3) GOTO 4

CALL MESSAG("CURVE YR MO DY TIME SSNO LAT LONG", 39,8.,7+25)

CALL MESSAG(" METH MAG REF FIG SPC",23,"ABUT", "ABUT")
CALL MESSAG("OBSERVED (0) ",14,-1.5,7.)
CALL INTNO(IDIST,"ABUT" ,"ABUT")

CALL INTNO(IYR,.75,"ABUT")

CALL INTNO(IMO,1.1,"ABUT")

CALL INTNO(IDY,1.45,"ABUT")

CALL INTNO(ITIME,1.8,"ABUT")

CALL INTNO(ISUN,2.6,"ABUT")

CALL REALNO(STLAT,1,3.25,"ABUT")

CALL REALNO(STLONG,1,4.,"ABUT")

CALL INTNO(IMETH,5.2,"ABUT")

CALL INTNO(IMAG,5.65,"ABUT")

CALL INTNO(IREF,5.85,“ABUT")

CALL MESSAG (UFIG,5,6.5,"ABUT")

CALL INTNO(ISC,7.25,"ABUT"

CALL MESSAG ("REFERENCE (R) ",14,-1.5,6.75)
CALL INTNO(JDIST,.2,"ABUT")

CALL INTNO(JYR,.75,"ABUT")

CALL INTNO(JMO,1.1,"ABUT")

CALL INTNO(JDY,1.45,"ABUT")

CALL INTNO(JTIME,1.8,"ABUT")

48



CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 804 69113 100 105 58.8 8 2 2154 3-3 8
REFERENCE (R) 0 O 0 O 0 a 0.0 60.0 o 00 is)
MODEL (M)

ALTITUDE (KM)

s oo

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 805 69113 628 105 «58.8 265.8 2 2154 3-3 8
REFERENCE (R) O oaoo 0 0.0 0.0 Qo oo ie)
MODEL (M) |
105
100
s
0
_
=
Fd
~ 0
w
8
Ss
as
=
5
“7
Ss
: f
. eval
cy

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
SOLAR PROTON EVENT

20



CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 995 699 302241 105 57.2 352.9 2 2151 3
REFERENCE (R) 994 728 312 65 $7.2 352.9 2 1151 30
MODEL (M)

ALTITUDE (KM)

2 3 4
LOG ELECTRON CONCENTRATION (N/CCI
CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC

OBSERVED (0) 996 69101 420 104 57.2 352.9 2 2151 3 ul
REFERENCE (R) 993 728 30 2202 65 57.2 352.9 2 1151 30
MODEL (M4)

105

1 i
« >

‘

ALTITUDE (KM)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
POLAR SUBSTORM

38



SELECTED DISTURBED D-REGION ELECTRON DENSITY PROFILES
Their relation to the undisturbed D region

by

L.F. McNamara
lonospheric Prediction Service
Sydney, Australia*

ABSTRACT. A large collection of experimental D-region electron density profiles has
been searched to find those profiles obtained under disturbed conditions, such as
eclipses and solar flares. Where possible, each disturbed profile has been associ-
ated with a suitable reference profile corresponding to the undisturbed ionosphere,
in order to determine the effect of the particular kind of disturbance. A model
electron density profile corresponding to the undisturbed ionosphere has also been
calculated in each case, and may be used as a reference profile when no suitable
experimental profile is available. All profiles have been plotted in a self-evident
format which makes it easy to see the effect of the disturbance considered.

Method

This report presents a collection of observed and model D-region electron density profiles that is
designed to illustrate the effects on the profiles of various kinds of phenomena. Detailed descriptions
of the phenomena are not given because these are available either in standard texts or in the journal
articles from which the profiles have been extracted. The main emphasis is one of illustrating the ef~
fects on the profile of the different phenomena. References are given in Part III of McNamara [1978].

All observed profiles have been extracted by computer from the D-region electron density data base

established by McNamara [1978], using the "special case parameter" associated with each profile. These
parameters and the associated physical phenomena are listed in Table 1.

Table 1. Phenomena Corresponding to the "Special Case" Parameter.

Parameter Physical Description

1 Solar eclipse

2 Solar flare

3 Solar X-ray event

4 Sudden ionospheric disturbance

5 Winter anomaly day

6 Mid-latitude particle precipitation

7 Solar proton event/solar cosmic ray event
8 Polar cap absorption

9 Auroral absorption

10 Daytime absorption event
11 Polar substorm

12 Polar radio blackout

13 Auroral arc conditions

NOTE: The descriptions are those given by the individual authors
and some are equivalent or have the same causal mechanism.

Each of the observed profiles has associated with it a model profile derived from all undisturbed
profiles in the complete data base [McNamara, 1979]. Many of the profiles also have associated with
then a reference profile corresponding to the undisturbed ionosphere and one usually observed just
before or just after the event in question.

*Guest worker at World Data Center A for Solar-Terrestrial Physics, NOAA, Boulder, Colorado.
1


OBSERVED

REFERENCE (R} O 0

MODEL (M)

ALTITUDE (KM)

OBSERVED

(0)

CURVE YR MO DY TIME SSNO LAT LONG ME:
143 64 10 27 1200 10 $8.8 265.8
g00 0 0.0 60.0

TH MAG REF
2 2140
o 00

FIG SPC
309
ie)

Z

L

a

Z
Va

(Oo)

2 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF

168 653 21 659 13 69.3 16.0

REFERENCE (R) 0 0 0 O O ie) 0.0 60.0

MODEL (M)

ALTITUDE (KM)

2 275
o 00

FIG SPC
09
i)

va

Pa

Lee

AURORAL ABSORPTION EVENT
31

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)



OBSERVED

REFERENCE (R) O

MODEL (M)

ALTITUDE (KM)

OBSERVED

REFERENCE (R} O
MODEL (M)

ALTITUDE (KM)

(0)

CURVE YR MO OY TIME SSNO

385 69113 752 10S

oa00

LAT

) 0.0

LONG METH MAG REF

58.8 265.8

0.0

2
ie)

2 138
i)

FIG SPC
27.6 8
ie)

|
\
}

WA

£

(0)

oaao0

2 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO DY TIME SSNO
390 69114 1708 105
te} 0.0

LAT

LONG METH MAG REF
58.8 265.8

0.0

2
ie]

2 138
oo

s

FIG SPC
27.8 8
0

Z

s

~

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
SOLAR PROTON EVENT

7

i]
|
{

t
;

f
f


OBSERVED

REFERENCE (R} O

MODEL (M)

ALTITUDE (KM)

OBSERVED

1
REFERENCE (R} 0 O O

MODEL (M)

ALTITUDE (KN)

(0)

CURVE YR MO DY TIME SSNO LAT LONG METH
384 69113 730 105 58.8 265.8 2

oa0a00

tt) 0.0 0.0 Q

MAG REF = =FIG SPC

244 20.3 8

oo

ie)

7

J

L

ZA

|

4

0 1

(0)

105

CURVE YR MO
388 69 11

oY
4
is)

(0)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

TIME SSNO LAT LONG METH
1530 105 S8.8 265.8 2

0 0.0 0.0 ie)

MAG REF FIG SPC

244 20.6 8

00

to)

if
7
yi

Wa

f
}

i

\f

vA

a 3 4
LOG ELECTRON CONCENTRATION (N/CCI

SOLAR PROTON EVENT
23

|
— —|


OBSERVED

(0)

CURVE YR
369 69

REFERENCE (R) 0 O

MODEL (M)

ALTITUDE (KM)

OBSERVED

MO
4
0

OY TIME SSNO
132319 «106
ao ie)

LAT
58.8
0.0

LONG METH MAG REF
265.8 2 294
0.0 Qo 00

FIG SPC
6 8
is)

LZ
V

(0)

CURVE YR
368 69

REFERENCE (R) O 0

MODEL (M)

ALTITUDE (KM)

TMO
4
Q

a 3 a
LOG ELECTRON CONCENTRATION (N/CC)

106

DY TIME SSNO
1
00 tt)

LAT
58.8
0.0

LONG METH MAG REF
265.8 2 294
0.0 o oo

FIG SPC
6 8
0

1

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
SOLAR PROTON EVENT

19

|

}
|



OBSERVED

(0)

CURVE YR
222 66

REFERENCE (R} O 0

MODEL (M)

ALTITUDE (KM)

OBSERVED

“Oo
2
Q

oY
22
0

TIME SSNO LAT LONG METH MAG REF
1929 31 69.3 16.0 2° 275
fo) 0 0.0 0.0 o 00

FIG SPC
a)
ie)

7

(0)

CURVE YR
223 66

REFERENCE (R} O ie)

MODEL (M)

ALTITUDE (KM)

105

“Oo
2
Qa

oY

2
ia)

4

a 3 4
LOG ELECTRON CONCENTRATION (N/CC}

TIME SSNO LAT LONG METH MAG REF
1S 31 69.3 16.0 2 275
ie) Q 0.0 60.0 o 00

FIG SPC
09
0

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
AURORAL ABSORPTION EVENT

32

|



CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 158 652 6 1200 12 74.7 265.1 1 168 27
REFERENCE (R) 160 652 11 1200 12 74.7 265.1 1 168 2.0 7
MODEL (M}

|
t
}
|

ALTITUDE (KM)

=o
mee
f

z

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 159 652 8 1200 12 74.7 265.1 1 168 27 i
REFERENCE (R) 160 652 11 1200 12 74.7 265.1 1 168 20 1
MODEL (M)

ALTITUDE (KM)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC}
SOLAR PROTON EVENT

4



CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 518 7210251208 61 37.1 353.3 2 148 4A 4 |
REFERENCE (R) 517 7210181200 61 37.1 353.3 2 148 4A 0 |
MODEL (M) |

ALTITUDE (KM)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 920 681 161514 103 37.9 284.5 2.12 43
REFERENCE (R) 0 0 0 0 0 ie) 0.0 60.0 0 00 0
MODEL (M)

ALTITUDE (KM)

2 3 4
LOG ELECTRON CONCENTRATION (N/CCI
SOLAR X-RAY EVENT
6



CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 69 621211427 30 69.3 16.0 20277 10 9
REFERENCE (R) 0 0 0 0 0 ie) 0.0 0.0 Oo 00 tt)
MODEL (M)

|
|

ALTITUDE (KM)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC!

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 70 62 12 14 2152 30 69.3 16.0 2 277 10 9
REFERENCE (R} 0 0 000 tt) 0.0 60.0 o 00 it)
MODEL (M)

ALTITUDE (KM)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
AURORAL ABSORPTION EVENT

29



CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 487 622 23140 42 64.9 212.1 S 1104 6 9
REFERENCE (R} 0 0 0 0 0 0 0.0 0.0 o 00 ie)

MODEL (H) '
= |

ALTITUDE (KM)

Z
ez | |

0 lo a 3 4
LOS ELECTRON CONCENTRATION (N/CC)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 282 67 1120137 97 69.3 16.0 2 1112 15 9
REFERENCE (R} O 0 00 0 0.0 0.0 o 00 tt)
MODEL (M)

1
ie}

uo

i. 2
. wa |

&
ph e"|

ALTITUDE (KM)

2 3 4
LOB ELECTRON CONCENTRATION (N/CC)
AURORAL ABSORPTION EVENT
26



OBSERVED

(0)

OBSERVED AND MODEL D-REGION ELECTRON DENSITY PROFILES

CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF
947 637 201503 28 58.6 285.6 2 1125

REFERENCE (R) 948 637 201610 28 58.8 265.8 2 $1125

MODEL (M)

ALTITUDE (KM)

OBSERVED

us

uo

FIG SPC
1301

$s 9

H
=

(0)

2 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF
B50 665 201045 41 36.8 21.9 2 5178

REFERENCE (R) 234 665 201130 41 36.8 21.9 2 178

MODEL (M)

ALTITUDE (KM)

FIG SPC
3 1
300

a

L

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
SOLAR ECLIPSE. #947 AND 850 DURING TOTALITY

3

|



Consideration of each of the following diagrams must include the suitability of the reference pro~
file and the error bars associated with each of the observed, reference and model profiles, as well as
the effect of the phenomena being considered. The error bars have been discussed by McNamara [1978;
1979]. A reasonable and simple estimate is +0.5 in logioN, or #1/2 division along the X axis of each
diagram. This corresponds to a factor of 3 in density. Differences between profiles which are less
than this value cannot be classed as "real". On some occasions, the model profile will agree better
with the disturbed profile than with the undisturbed one. Under these circumstances the model profile
should be ignored. Such failures of the model are only to be expected, especially when the differences
between profiles are of the order of the error bar.

The first figure may be used to illustrate the presentation of the data. In general, there will be
three profiles: an observed (disturbed) profile (0), a reference profile (R), and a model profile (M).
In some cases, no suitable reference profile was available in the data base, so the observed profile can
be compared only with the model profile. The reliability of the model profile is thus an important
consideration.

Each observed or reference profile is identified by a unique three-digit curve identification
number. In the first figure these are 947 and 948, respectively. These numbers are sufficient to
extract the corresponding data from the computer files DB1 and DB3, which comprise part of the data base
[McNamara, 1978], including the parameters listed at the top of the figure. These parameters give the
relevant conditions of the experiment--the year, month, day, local time, monthly average sunspot number,
latitude and longitude of the observing station, method of observation, and a magnetic index. The METH
parameter may assume one of the following five values: (1) partial reflection; (2) rocket; (3) wave
interaction; (4) LFVLF reflection; and (5) other, @.g., incoherent scatter. The magnetic index MAG is
1 for magnetically quiet days (approximately Ap <25) and 2 for disturbed days. The parameter REF is a
literature reference number, which may be used to find the corresponding journal or report reference in
the computer file DB2 of the data base or in Part III of McNamara [1978]. The figure number FIG is the
figure number in that reference.

The last parameter in each line is SPC, the special case parameter described in Table 1. It is
zero for the reference profile. The model profile is calculated using the year, month, day, local time,
sunspot number, latitude and longitude corresponding to the observed profile. The model for magnetical-
ly quiet days differs slightly from that for magnetically disturbed days. It is most likely to fail at
night and at high latitudes, because the data base includes fewer undisturbed nighttime, high-latitude
profiles than daytime mid- and low-latitude profiles.

The computer program DISPPLT used to prepare the data for plotting is listed in Appendix A. The
actual plots were made using the software package of the Integrated Software Systems Corporation called
"DISPLAA" and the program DREGPLT listed in Appendix B.

Some of the phenomena corresponding to different special case parameters have in fact the same
causal mechanism. The associated plots have therefore been grouped together and described by the causal
mechanism. Thus solar flare (2), solar x-ray event (3) and sudden ionospheric disturbance (4) are
grouped under solar x-ray event; solar proton event (7), polar cap absorption (8) and polar radio black-
out (12) are grouped under solar proton event.

Not all available plots are included, since some are very similar to each other and contain no
further useful information. In some cases, consecutive plots show the time development of the event.

Acknowledgment
The plotting program DREGPLT was written and run by Marcus 0. Ertle, NOAA/EDIS.

REFERENCES

MCNAMARA, L.F. 1978 lonospheric D-Region Profile Data Base--A Collection of
Computer-Accessible Experimental Profiles of the D and
Lower E Region. Report UAG-67, World Data Center A for
Solar-Terrestrial Physics, Boulder, Colorado.

MCNAMARA, L.F. 1979 Statistical Model of the D-Region. Submitted to Radio
Science.



CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC |
OBSERVED (0) 938 667 17 1200 50 74.7 265.1 1 168 4 7 |
REFERENCE (R) 250 667 15 1200 50 74.7 265.1 1 168 40 |
MODEL (M)

ALTITUDE (KM)

|

|
|
|
|
|

- Z

2 3 4
LOG ELECTRON CONCENTRATION (N/CC)

MO DY TIME SSNO LAT LONG METH MAG REF FIG SPC
1 291030 91 ~-66.7 140.0 2°10 19.6 7
aoo0o0 ie) 0.0 0.0 9 00 0

CURVE YR
OBSERVED (0) 932 67
REFERENCE (R} OQ 0
MODEL (M)

ALTITUDE (KN)
8
hs.
‘Se

a 3 4 s
LOG ELECTRON CONCENTRATION (N/CC) |
SOLAR PROTON EVENT |
15


CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 409 703 7 1309 106 40.8 2682.1 3 114 21
REFERENCE (R} 417 703 111310 106 40.8 262.1 3 114 20
MODEL (M)

|
|

ALTITUDE (KM)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 410 703 7 1328 106 40.8 282.1 3 1114 21
REFERENCE (R) 417 703 111310 106 40.8 282.1 3 114 2.0
MODEL (M)

ALTITUDE (KM)

8 1 a 3 4 s i
LOG ELECTRON CONCENTRATION (N/CC! |

SOLAR ECLIPSE. MAXIMUM OBSERVATION 90% AT 1338 LT |

5 |



OBSERVED

REFERENCE (R) 925 69112 400 105

MODEL (M)

ALTITUDE (KM)

OBSERVED

REFERENCE (R} 925 69 11 2 400 10S
MODEL (mM)

ALTITUDE (KM)

(0)

CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF
58.8 265.8
58.8 265.8

375 69 11 2 500 105

1
1

217
217

FIG SPC

16.6
16.6

8
0

Le

=a

(0)

376 69112 600 105

1
1

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF
$8.8 265.8
58.8 265.8

217
217

FIG SPC

16.6
16.6

8
0

a

VA

SOLAR PROTON EVENT
21

a 3 4
LOS ELECTRON CONCENTRATION (N/CC)



WORLD DATA CENTER A
or

Solar-Terrestrial Physics

REPORT UAG-69

SELECTED DISTURBED D-REGION
ELECTRON DENSITY PROFILES

Their relation to the undisturbed D-region

Mi

OCTOBER 1978



APPENDIX A (Continued)

c LOOK FOR REFERENCE PROFILE

81 READ (9,119) KKW
IF (EOF (9) .EQ.8.) GO TO 87

c NO REFERENCE PROFILE |
DO 86 I=1,12 t
FLN(I)=0. i

86 H(L

KZ=! |
WRITE(19,1190) KZ |
WRITE(10,120) (FLN(I) ,H(I) ,I=1,6)
REWIND 9
GO TO 88

87 CONTINUE

82

READ (9,126) (FLN(J) ,H(J) ,J=J1,J2)
IF(H(J2).NE.@.) GO TO 82 |
IF (KKW.NE.IKR) GO TO 81 f
WRITE(19,119) KKW }
WRITE(10,120) (FLN(J) ,H(J) ,J=1,32) i
Cc LOOK FOR MODEL PROFILE .... NO NEED TO REWIND |
88 CONTINUE
KKX=1006+KKV
83 READ (9,116) KKK

84

J2=J145

READ (9,120) (FLN(J) ,H(J) ,»J=31,J2)
IF(H(J2).NE.@.) GO TO 84

IF (KKK.NE.KKX) GO TO 83
WRITE(19,118) KKX

WRITE(1@,126) (FLN(J) ,H(J) ,J=1,J2)

REWIND 9
GO TO 78

9

CONTINUE

c
9
Cc
Cc
c

RETURN FOR NEXT SPECIAL CASE PARAMETER
Go TO 5

END

SUBROUTINE PROFILE (IY,IM,ID,LH,LM,ISSNO,FLAT,WLON, DEN)

S/R TO CALCULATE A D-REGION PROFILE USING THE STATISTICAL MODEL
WITH COEFFICIENTS A,B,C,D,E

MAGNETICALLY QUIET DAYS ONLY. DENSITY INCREASES BY A FACTOR
OF ABOUT TWO ON DISTURBED DAYS

DIMENSION A(8) ,B(8) ,C(8) ,D(8) ,E(8) -F(8) ,T(8) ,DEN(8)

A = CONSTANT TERM 90(5)55 KM

DATA A / 3.21,2.32,1.87,1.64,1.49,1.25,-81,.82 /

Qa aaaa

Cc B = COS KI TERM
DATA B / .80,.86,1.93,1.08,.86,.57, 266,477 /
c C = GEOMAGNETIC LATITUDE TERM

DATA C / O.,+23,.16,.69,.17,+58,.78,0. /
c D = SUNSPOT TERM
DATA D / 6440134420, 625,613,226,426,0. /

c E = SEASONAL TERM
DATA E / .15,.16,-17,0.,~.05,04,+12,0. /
c RESIDUAL STANDARD DEVIATION
DATA F / 653, 65404644397 233,032, 440,047 /
c
PI=3,141592654
Cc PREPARE TO CHANGE TO GEOMAGNETIC LATITUDE AT WLON DEGREES (WEST)
DEGRAD=PI/186.
1.5*DEGRAD
OS (THETGM)
SINGM=SIN (THETGM)
PHIGM=291.*DEGRAD
c

44



CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 324 683 251010 105 45.4 284.1 11127 01 2
REFERENCE (R) 323 683 25 945 105 45.4 284.1 11127 01 0
MODEL (M)

}

8

f

ALTITUDE (KM)

zx

8 1

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIL

& SPC
OBSERVED (0) 64 6810211239 110 40.8 282.1 3 1S 1
1

2
REFERENCE (R) 63 68 10 21 1210 110 40.8 262.1 3 1s to)
MODEL (M}

DN

ALTITUDE (KM)
a

Z|

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
SOLAR X-RAY EVENT
8



CURVE YR MO DY TIME SSNO LAT

LONG METH MAG REF

OBSERVED (0) S41 721 191147 71 37.9 284.5 1 143
REFERENCE (R) 540 721 191109 71 37.9 284.5 1 143

MODEL (M)

FIG SPC
S.12 2
$.12 0

ALTITUDE (KM)

CURVE YR MO DY TIME SSNO LAT

a 3 4
LOS ELECTRON CONCENTRATION (N/CC)

LONG METH MAG REF

OBSERVED (0) 922 623 1 1147 40 45.4 284.1 1 06114
REFERENCE (R) 54 623 0 1200 40 45.4 284.1 10115

MODEL (M)

FIG SPC
8 4
3.6 0

ALTITUDE {KN}

SOLAR X-RAY EVENT
7

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)



CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 969 733 24 2250 44 65.3 212.2 1 2153) 4.3 13
REFERENCE (R} 966 733 242100 44 65.3 212.2 1 2153 4.3 0
MODEL (M)

ALTITUDE (KN)

wa

3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 970 733 26 2238 44 65.3 212.2 1 2183) 4.7 13
REFERENCE (R} 0 0 0 0 O 0 0.0 0.0 Oo 00 it)
MODEL (TM}

. =a

ALTITUDE (KM)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

AURORAL ARC CONDITIONS. #966 PRE-ARC; #969 SUBSEQUENT AURORAL ACTIVITY (SEE
PREVIOUS PAGE).
40



CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF FIG SPC

OBSERVED (0) 12 577 4 1216 191 588 265.8 2 214 99 |
REFERENCE (R) 41 616 0 1200 56 45.4 284.1 1 114 90 |
MODEL (1) |
|

105 i

° '

ALTITUDE (KM)
a 8
i—e}

a
pa

: = ?

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 60 628 18809 35 69.3 16.0 2 277 6 9
REFERENCE (R} O oa ie} (0) 0.0 0.0 Oo oOo ie)
MODEL (TM)

Ui
ie)

105

ALTITUDE (KM)

« —

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
AURORAL ABSORPTION EVENT

28



Cc

15
20

c

Coser
Cc

Cases

aanaaa

fo)

APPENDIX A (Continued)

CONVERT TO GEOMAGNETIC LATITUDE
COLAT=(9.-FLAT) *DEGRAD

ALON=WLON*DEGRAD

COSA=COS (COLAT) *COSGM+SIN (COLAT) *SINGM*COS (ALON=PHIGM)
GEOML=99 .~ACOS (COSA) /DEGRAD
TLAT=SIN(2.*GEOML*PI/189.)

1+9.9* (SIN (GEOML*PI/186.)) **4

FMON=(FMON~9.5) /12.*2.*PI

FMON=COS (FMON)

IF (FLAT.LT.@.) FMON=~FMON

SEASON=FMON

IID=ID

IF(ID.EQ.9) IID=15

CALL COSKI(IY,IM,IID,LH,LM,FLAT,WLON, COSX)

DO 20 I=1,8

IHT=99~(I~1) *5

T(1)=A(I)

(2) =COSX*B (I)

XLAT=TLAT

IF(I.GE.5) XLAT=ALAT

T(3)=XLAT*C (I)

T (4) =ISSNO*.91*D (I)

T(5) =SEASON*E (I)

T(6)=F(I)

ED=0.

DO 15 J=1,5

ED=ED+T (J)

DEN(I)=ED

CONTINUE

RETURN

END

SUBROUTINE COSKI(LLY,LMON, LDAY ,LH,LM,DLAT,DLON, COSX)
MODIFIED VERSION OF NEWBERN SMITH'S COS(SOLAR ZENITH ANGLE) S/R
IT IS ASSUMED THAT THE GIVEN LOCAL TIME WAS MEASURED AT THE
NEAREST STANDARD MERIDIAN I.E. THE NEAREST MULTIPLE OF 15 DEGREES
«LY = LOCAL YEAR LAST TWO DIGITS ONLY

LMON = LOCAL MONTH

LDAY = LOCAL DAY

LH = LOCAL HOUR

LM = LOCAL MINUTE

DLAT = NORTHERN LATITUDE IN DEGREES

DLON = WEST LONGITUDE IN DEGREES

DIMENSION NDO(12) ,NDL(12)

DATA DTR / 9.G1745329 /

DAY NUMBERS FOR ORDINARY AND LEAP YEARS

DATA NDO,NDL / 1,32,60,91,121,152,182,213,244,274,305,335,
A 1,32,61,92,122,153,183,214,245,275,366,336/

AMPLITUDES OF FOURIER EXPANSION COEFFICIENTS 1955 EPOCH
DATA P1,P2,P3,P4,P6 /
A 6.917203534,0.634407068 ,0.051610602,0.068814136,8.193221204 /
IF(MOD(LY,4) .EQ.6) LD=NDL(LMON) +LDAY~1

IF (MOD(LY,4) .NE.@) LD=NDO(LMON)+LDAY~1

LY=LLY-68

MER=INT ( (DLON+7.5) /15) *15.

FIND TIME OF THE EQUINOX FOR THAT YEAR

IF(LY) 7,7,9

IL=INT((3.1-LY) /4.)

Go TO 11
INT( (LY+.1)/4.)
78194+0.24225*LY+IL

D+ (LH+LM/6@.+MER/15.) /24.
TE=TD-DT
DECLINATION

DCL=23.256*SIN(P1* (TE-~82.242) )+.381*SIN(P2* (TE~44.855) )+.167*SIN
1(P3* (TE~23.355) )~.613*SIN(P4* (TE+11.97) )+.O11*SIN (P6* (TE~10.41) )
2+ .339137

45

t
,
|

7
f



OBSERVED

(0)

CURVE YR
248 66

REFERENCE (R} O 0

MODEL (4)

ALTITUDE (KM)

OBSERVED

MO
6
ie)

1200 45 69.3 16.0
ie) 0 0.0 60.0

Oy TIME SSNO LAT LONG METH MAG REF
26
ie)

2 1112
a 00

i

FIG SPC
89
tt)

L

{

Lv

[>

4

(0)

CURVE YR
272 67

REFERENCE (R} O 0

MODEL (mM)

ALTITUDE (KN)

“oO
3
Q

oY
1
fe)

2 3 4
LOG ELECTRON CONCENTRATION (N/CCI

TIME SSNO LAT LONG METH HAG REF

46 82 69.3 16.0
ie) 0 0.0 60.0

2 5175
o oOo

FIG SPC

09
Q

Z

if

eal

Zo

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

AURORAL ABSORPTION EVENT
27



CURVE YR MQ OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 809 66 11121040 70 -32.0 308.0 1 1155 7 1
REFERENCE (R) 808 66 110 1045 70 ~-32.0 308.0 1 1155 170
MODEL (M)

s ea |

ALTITUDE (KM)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC}

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 944 66 11121040 70 -32.0 308.0 2 189 41
REFERENCE (R) 943 66 11121230 70 ~-32.0 308.0 2 189 4 0
MODEL (mM)

ALTITUDE (KN)

: 7
—S

a a 4
LOG ELECTRON CONCENTRATION (N/CC)
SOLAR ECLIPSE. #809 AND 944 DURING TOTALITY
4



CURVE

YR MQ OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 366 692 25 1 1100 50.0 15.0 2 278 1 8
REFERENCE (R} 0 O08 0 0 0 0 0.0 0.0 Oo oo ie)
MODEL (M)
105

ALTITUDE (KM)
a

Ss =|

a 3 4
LOG ELECTRON CONCENTRATION (N/CC}

OY TIME SSNO LAT LONG METH MAG REF FIG SPC
25 2237 «110 =50.0 15.0 2 278 1 8
00 Q 0.0 0.0 6 oO ie)

Z |
: |
t

OBSERVED (0) 365 69
REFERENCE (R} O 0
MODEL (M)

CURVE YR MO
2
0

ALTITUDE (KN)
a

cal

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
SOLAR PROTON EVENT
18



CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 353 68 10 242150 110 67.8 20.4 2 247 19
REFERENCE (R) O O 0 0 0 0 0.0 60.0 o 00 ie)
MODEL (HM)

ALTITUDE (KM)

|
|
|

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 362 691 182024 110 69.3 16.0 2 247 1 9
REFERENCE (R} 0 O 0 0 0 0 0.0 60.0 Oo 00 oO
MODEL (M) |

105

s +7
0 |

N

ALTITUDE (KM)
a

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
AURORAL ABSORPTION EVENT

35

|


OBSERVED (0)
REFERENCE (R) O

MODEL (13 5

ALTITUDE (KM)

OBSERVED

110

105

(0)

CURVE YR MO DY TIME SSNO

483.728 7 445 66
oo0600 ie)

LAT
-69.0
0.0

LONG METH
33.0 2
0.0 ie}

MAG REF
1 100
oo

FIG SPC
tou
ia)

REFERENCE (RI O
MODEL (M1 ys

ALTITUDE (KM)

U9

105

ooo

1

23354 36
oo i)

LAT
69.0
o.0

LONG METH
39.0 2
0.0 oO

2 3 4
LOG ELECTRON CONCENTRATION (N/CC}

CURVE YR MO DY TIME SSNO
488 738

MAG REF
1 100
oo

FIG SPC
7
a

q

2

x

man

a 3 4
LOG ELECTRON CONCENTRATION (N/CC}

POLAR SUBSTORM

36

|
|
}
|
|
|
|
|

1
i}


UAG Series of Reports

Prepared by World Data Center A for Solar-Terrestrial Physics, NOAA, Boulder, Colorado, U.S.A.

These reports are for sale through the National Climatic Center, Federal Building, Asheville, NC 28801, Attn: Publications.
Subscription price: $25.20 a year; $12.00 additional for foreign mailing; single copy price varies. These reports are issued
on an irregular basis with 6 to 12 reports being issued each year. Therefore, in some years the single copy rate will be less
than the subscription price, and in some years the single copy rate will be more than the subscription price. Make check or
money order payable to: Department of Commerce, NOAA.

Some issues are now out of print and are available only on microfiche as indicated. Requests for microfiche should be sent
to World Data Center A for Solar-Terrestrial Physics, NOAA, Boulder, CO 80302, with check or money order made payable to
Department of Commerce, NOAA.

UAG-1 “IQSY Night Airglow Data", price $1.75.

UAG-2 "A Reevaluation of Solar Flares, 1964-1966", price 30 cents.

UAG-3 “Observations of Jupiter's Sporadic Radio Emission in the Range 7.6-41 MHz, 6 July 1966 through 8 September 1968",
microfiche only, price 45 cents.

UAG-4 "Abbreviated Calendar Record 1966-1967", price $1.25.

UAG-5 “Data on Solar Event of May 23, 1967 and its Geophysical Effects", price 65 cents.

UAG-6 "International Geophysical Calendars 1957-1969", price 30 cents.

UAG-7 “Observations of the Solar Electron Corona: February 1964-January 1968", price 15 cents.

UAG-8 "Data on Solar-Geophysical Activity October 24-November 6, 1968", price (includes Parts 1 & 2) $1.75.
UAG-9 "Data on Cosmic Ray Event of November 18, 1968 and Associated Phenomena", price 55 cents.

UAG-10 “Atlas of Ionograms", price $1.50.

UAG-11 "Catalogue of Data on Solar-Terrestrial Physics" (now obsolete).
UAG-12 “solan-Geophysical Activity Associated with the Major Geomagnetic Storm of March 8, 1970", price (includes Parts
1-3) $3.00.

UAG-13 "Data on the Solar Proton Event of November 2, 1969 through the Geomagnetic Storm of November 8-10, 1969, price 50 cents.
UAG-14 “An Experimental, Comprehensive Flare Index and Its Derivation for ‘Major’ Flares, 1955-1969", price 30 cents.
UAG-15 “Catalogue of Data on Solar-Terrestrial Physics" (now obsolete).

UAG-16 "Temporal Development of the Geographical Distribution of Auroral Absorption for 30 Substorm Events in each of IQSY
(1964-65) and IASY (1969)", price 70 cents.

UAG-17 "Ionospheric Drift Velocity Measurements at Jicamarca, Peru (July 1967-Narch 1970)", microfiche only, price 45 cents.

UAG-18 "A Study of Polar Cap and Auroral Zone Magnetic Variations", price 20 cents.

UAG-19 "Reevaluation of Solar Flares 1967", price 15 cents.

UAG-20 "Catalogue of Data on Solar-Terrestrial. Physics" (now obsolete).

UAG-21 “Preliminary Compilation of Data for Retrospective World Interval duly 26 - August 14, 1972", price 70 cents.

UAG-22 “Auroral Electrojet Magnetic Activity Indices (AE) for 1970", price 75 cents.

UAG-23 "U.R.S.I, Handbook of Ionogram Interpretation and Reduction", price $1.75.

UAG-24 “Data on Solar-Geophysical Activity Associated with the Major Ground Level Cosmic Ray Events of 24 January and
1 September 1971", price (includes Parts 1 and 2) $2.00.

UAG-25 “Observations of Jupiter's Sporadic Radio Emission in the Range 7.6-41 MHz, 9 September 1968 through 9 December 1971",
price 35 cents.

UAG-26 “Data Compilation for the Magnetospherically Quiet Periods February 19-23 and November 29 - December 3, 1970",
price 70 cents.

UAG-27 “High Speed Streams in the Solar Wind", price 15 cents.

UAG-28 “Collected Data Reports on August 1972 Solar-Terrestrial Events", price (includes Parts 1-3) $4.50.

UAG-29 “Auroral Electrojet Magnetic Activity Indices AE (11) for 1968", price 75 cents.

UAG-30 "Catalogue of Data on Solar-Terrestrial Physics", price $1.75.

UAG-31 "Aurgral Electrojet Magnetic Activity Indices AE (11) for 1969", price 75 cents.

UAG-32 "Synoptic Radio Naps of the Sun at 3.3 mm for the Years 1967-1969", price 35 cents.

UAG-33 “Auroral Electrojet Magnetic Activity Indices AE (10) for 1967", price 75 cents.

UAG-34 "Absorption Data for the IGY/IGC and IQSY", price $2.00.

UAG-35 "Catalogue of Digital Geomagnetic Variation Data at World Data Center A for Solar-Terrestrial Physics", price 20 cents.

UAG-36 “An Atlas of Extreme Ultraviolet Flashes of Solar Flares Observed Via Sudden Frequency Deviations During the ATM-SKYLAB
Missions", price 55 cents.

UAG-37 “Auroral Electrojet Magnetic Activity Indices AE (10) for 1966", price 75 cents.

UAG-38 "Master Station List for Solar-Terrestrial Physics Data at WOC-A for Solar-Terrestrial Physics", price $1.60.

UAG-39 "Auroral Electrojet Magnetic Activity Indices AE (11) for 1971", by Joe Haskell Allen, Carl C. Abston and Leslie D.
Morris, National Geophysical nd Solar-Terrestrial Data Center, Environmental Data Service, February 1975,
144 pages, price $2.05.

UAG-40 "H-Alpha Synoptic Charts of Solar Activity For the Period of Skylab Observations, May, 1973-March, 1974", by Patrick
S. McIntosh, NOAA Environmental Research Laboratories, February 1975, 32 pages, price 56 cents.

UAG-41 "H-Alpha Synoptic Charts of Solar Activity During the First Year of Solar Cycle 20, October, 1964 - August, 1965",
by Patrick S. McIntosh, NOAA Environmental Research Laboratories, and Jerome T. Nolte, American Science and
Engineering, Cambridge, Massachusetts, March 1975, 25 pages, price 48 cents.

UAG-42 "Observations of Jupiter's Sporadic Radio Emission in the Range 7.6-80 “Hz 10 December 1971 through 21 March 1975",
by James W. Warwick, George A. Dulk, and Anthony C. Riddle, Department of Astro-Geophysics, University of
Colorado, Boulder, Colorado 80302, April 1975, 49 pages, price $1.15.

UAG-43 “Catalog of Observation Times of Ground-Based Skylab-Coordinated Solar Observing Programs", compiled by Helen E.
Coffey, World Data Center A for Solar-Terrestrial Physics, May 1975, 159 pages, price $3.00.

UAG-44 “Synoptic Maps of Solar 9.1 cm Microwave Emission from June 1962 to.August 1973", by Werner Graf and Ronald N,
Bracewell, Radio Astronomy Institute, Stanford University, Stanford, California 94305, May 1975, 183 pages,
price $2.55.

UAG-45 "Auroral Electrojet Magnetic Activity Indices AE (11) for 1972", by Joe Haskell Allen, Carl C. Abston and Leslie
D. Morris, National Geophysical and Solar-Terrestrial Data Center, Environmental Data Service, May 1975,
144 pages, price $2.10.

UAG-46 "Interplanetary Magnetic Field Data 1963-1974", by Joseph H, King, National Space Science Data Center, NASA Goddard
Space Flight Center, Greenbelt, Maryland 20771, June 1975, 382 pages, price $2,95.

UAG-47 “Auroral Electrojet Magnetic Activity Indices AE (11) for 1973", by Joe Haskell Allen, Carl C, Abston and Leslie
D. Morris, National Geophysical and Solar-Terrestrial Data Center, Environmental Data Service, June 1975,
144 pages, price $2.10.



CURVE YR MG DY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 16 5811 262100 181 69.7 19.0 3.4172 6 12
REFERENCE (R} 0 0 0 0 O ie) 0.0 0.0 0 00 ie)

MODEL (M) E

|
FE
|

ALTITUDE (KM)

a. |

s wo

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 20 8 26 2120 151 69.7 19.0 3.4172 ib 12
REFERENCE (R} 0 O 0 0 O 0 0.0 60.0 o 00 ie) q
MODEL (M) |
105
}
100 |
6
20
5s

ALTITUDE (KM)
a 8

}
A
so a |

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

SOLAR PROTON EVENT
24

|
i
;
}
|
1
|
1
t
|
|
|
|


TABLE OF CONTENTS

Page

Abstract... . 2... 2.200. Fe ee ee eee ee 1
Method... ee ee ee wee Be eee sed
Acknowledgment... . ee ee 4
OBSERVED AND MODEL D-REGION ELECTRON DENSITY PROFILES ........ 3
Plots During Solar Eclipse ...........00- eee ee 3B
Plots During Solar X-ray Event ....... Se ee ~. 6
Plots During Winter Anomaly Day. .........2.000 » + 10
Plots During Solar Proton Event. ...... se ee wees 1
Plots During Auroral Absorption Event. ......... wee 25
Plots During Polar Substorm. ...... eee see 36
Plots During Auroral Arc... 1... Se ee eee » 39
APPENDIXA............ eee bee eee wee ane 41
Program DISPPLT. 2... ee ee ee ee ee ee ee Al
File PRODAT. . 2... 0... Bee ee ee 46

Subroutine PROFILE for magnetically disturbed days ....... 47
APPENDIXB. 2... . ee. ee ae Se ee ee eee eee eee 48



APPENDIX A (Continued)

c SEE IF CURVE IS A REQUIRED REFERENCE CURVE
DO 73 L=1,LL
IF (KURVE.NE.IKRL(L)) GO TO 73
c CURVE IS REQUIRED
WRITE(9,110) KURVE
WRITE(9,126) (FLN(J) ,H(J) ,J=1,32)
GO TO 49
73 CONTINUE
GO TO 48
c CALCULATE STATISTICAL MODEL PROFILE
74 CONTINUE
DO 75 L=1,LL
ABSLAT=ABS (FFL (L) )
FLON=36.~FWL (L)
CALL PROFILE(IYL(L) ,IML(L) ,IDL(L) ,LHL(L) ,LML(L) ,ISL(L) ,
A_ABSLAT, FWL(L) ,DEN)
: ARRANGE THE POINTS IN ASCENDING HEIGHT ORDER :
DO 77 I=1,8
J=8*I+1
ADEN (I) =DEN (J)
77 AHT(I)=HT (J)
c WRITE OUT THE STATISTICAL PROFILE
NUML=1999+NUM (L)
WRITE(9,116) NUML
WRITE(9,126) (ADEN(I) ,AHT(I) ,I=1,8)

c PRINT 300,IYL(L) ,IML(L) ,IDL(L) ,LHL(L) ,LML(L) ,ISL(L) ,
c A PFL(L) ,FLON

75 CONTINUE

300 FORMAT (//* 19*,12,313,12,15,2F6.1)

c

c PRINT 105,NUM(L) ,HT

c PRINT 156,DEN

c PRINT 151,0BS

165 FORMAT (1X,*PROFILE NUMBER = *,14,/1X,* HEIGHTS *,8F6.1)
158 FORMAT (1X,*CALCULATED* , 8F6.1)
isl FORMAT(1X,* OBSERVED *,8F6.1)

c
REWIND 8
REWIND 9
REWIND 11
c WRITE OUT THE DATA IN A MORE SUITABLE FORM
78 READ (8,198) KURVE,IY,IM,ID,LH,LM,ISSNO,FLAT,FLON, METH,

A _MAG,KOL,IREF,IFIG,ISPEC,IKR
IF (KURVE.EQ.9999) GO TO 99
79 READ(9,118) KKV
5

80 J1=I1+6
J2=J1+5
READ (9,128) (FLN(J) ,H(J) ,J=J1,32)
IF(H(J2) .NE.B.) GO TO 80
c SEE IF THIS IS THE REQUIRED DISTURBED PROFILE
IF (KKV.NE.KURVE) GO TO 79
WRITE(10,100) KURVE,IY,IM,ID,LH,LM,ISSNO, FLAT, FLON, METH,
A MAG ,KOL,IREF,IFIG, ISPEC,IKR
c FIND HEADER CARD FOR REFERENCE PROFILE AND WRITE IT OUT
93 READ (11,100) KURVE,IY,IM,ID,LH,LM,ISSNO, FLAT, FLON, METH,
A MAG ,KOL, IREF,IFIG, ISPEC
IF (KURVE.EQ.9999) GO TO 94
IF (KURVE.NE.IKR) GO TO 93
WRITE (1,160) KURVE,IY,IM,ID,LH,LM,ISSNO,FLAT,FLON, METH,
A MAG, KOL, IREF,IFIG, ISPEC
REWIND 11
GO TO 95
94 REWIND 11
WRITE(19,432) (BLANKS(II) ,II=1,20)
432. FORMAT (2@A4)
95 CONTINUE
WRITE(19,110) KKV
WRITE(1@,128) (FLN(J) ,H(J) ,J=1,J2)
REWIND 9

43



APPENDIX A

The computer program DISPPLT selects all D-region profiles obtained under specified disturbed
conditions, and prepares the data for the plotting program DREGPLT. The input files TAPES and TAPE6 are
the “header card file," DB1, and the “digitized data file," DB3, described by McNamara [1978]. The
input file TAPE7 is the file PRODAT which contains the corresponding values of curve identification
number for the disturbed and reference profiles. The program expects the special case parameter to be
entered on the keyboard in FORMAT (12). The output file, TAPE10, should be saved for plotting.

The program listing of DISPPLT below must be changed slightly for magnetically disturbed days to be
considered. The statement IF (MAG.NE.1) GO TO 15 must be changed to IF (MAG.NE.2) GO TO 15. The sub-
routine PROFILE must also be replaced by the version appropriate to magnetically disturbed days, and the
first six lines of TAPE7 deleted.

PROGRAM DISPPLT (INPUT, OUTPUT, TAPES , TAPE6 , TAPE7 , TAPES,
A TAPES ,TAPE1@,TAPE11)
PROGRAM TO SELECT "UNUSUAL CONDITION" PROFILES AND THE CORRESPONDING
REFERENCE PROFILE (IF ANY) AND CALCULATE STATISTICAL MODEL PROFILE
THIS VERSION MAG QUIET DAYS ONLY
TAPES = DBL
TAPE6 = DB3
TAPE7 = TABLE OF CORRESPONDING REFERENCE PROFILES
TAPES = OUTPUT FILE .. HEADER CARDS
TAPE11 = OUTPUT FILE . HEADER CARDS FOR REFERENCE PROFILES
TAPES = OUTPUT FILE .. DIGITIZED DATA
TAPE1@ = OUTPUT ARRANGED IN MORE SUITABLE ORDER
DIMENSION NUM(10@) ,IYL(18G) ,IML(10@) ,IDL(10) ,LHL (190) ,LML(190),
A ISL (10) ,FFL(190) ,FWL(10@) ,DEN(8) ,HT(8) ,H(205),
A_FLN(200) ,OBS(8) ,IK(106) ,IR(199) , IRFL (100) , IKRL(180)
DIMENSION ADEN (8) ,AHT(8)
DIMENSION BLANKS (20)
DATA BLANKS / 26*4H

c

c READ IN THE REFERENCE PROFILES CORRESPONDING TO EACH "UNUSUAL PROFILE"
READ (7,408) (IK(L),IR(L) ,L=1,64)

400 FORMAT (2613)

Cc
Cc HEIGHTS AT WHICH DENSITIES ARE TO BE FOUND
DO 3 I=l1,
3 HT(I)=96.-5.* (I~1)
c
c*** WHAT SPECIAL CASE PARAMETER IS TO BE CONSIDERED??????????
5 READ 16, ISPECL

IF(ISPECL.EQ.%) STOP

PRINT 260,ISPECL
208 FORMAT (///1X,*SPECIAL CASE PARAMETER =*,13)
16 FORMAT (12)

REWIND 5

REWIND 6

REWIND 7

REWIND 8

REWIND 9

REWIND 11

c
c READ IN THE HEADER CARD FILE TO FIND ALL REQUIRED PROFILES

L=@
15 READ (5,108) KURVE,IY,IM,ID,LH,LM,ISSNO, FLAT, FLON, METH,

A _MAG,KOL, IREF ,IFIG,ISPEC
108 FORMAT (15,413,12,14,2F6.1,212,13,14,A5,213)

IF (KURVE.EQ.9999) GO TO 25
C=~~ MAGNETICALLY QUIET DAYS ONLY

IF(MAG.NE.1) GO TO 15

IF(ISPEC.EQ.ISPECL) GO TO 26

GO TO 15
20 L=L+1

IF(L.LE.109) GO TO 22

PRINT 23
23 FORMAT(* TOO MANY CASES*)

STOP
22 NUM(L) =KURVE

41



OBSERVED

MODEL (M)

ALTITUDE (KM)

OBSERVED

(0)

CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF

268 67 1 31 1200 75
REFERENCE (R) 220 65 12 15 1200 24

37.9 284.5
37.9 284.5

2
2

1119
190

FIG SPC
3.5
1 0

eA
eo

(0)

288 681 8

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
CURVE YR MOQ OY TIME SSNO LAT LONG METH MAG REF

1500 102 40.8 282.1

REFERENCE (R) 287 67 12 14 900 102 40.8 282.1

MODEL (M)

ALTITUDE (KM)

3
3

1116
1116

FIG SPC
3.05
300

v4
A

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

WINTER ANOMALY DAY
1

i
f

i
i



OBSERVED (0)

CURVE YR MO DY TIME SSNO
992 71123 1258 69

REPERENCE (R) 991 71 121 1334 69
MODEL (M4)

ALTITUDE (KM)

OBSERVED

ie)

LAT LONG METH MAG REF

57.2 352.9 2
57.2 352.9 2

115)
1151

fe

FIG SPC

25
20

(0)

CURVE YR MO OY TIME SSNO
545 71 12 12 1200 69

2 3 4
LOG ELECTRON CONCENTRATION (N/CCI

LAT LONG METH MAG REF

37.9 284.5 1

REFERENCE (R) 544 71129 1200 69 40.1 271.8 1

MODEL (M)

ALTITUDE (KM)

10

143
143

FIs
S.7
5.7

sec
5
0

va

WINTER ANOMALY DAY

13

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)



WORLD DATA CENTER A
for
Solar-Terrestrial Physics

REPORT UAG-69

SELECTED DISTURBED D-REGION
ELECTRON DENSITY PROFILES

Their relation to the undisturbed D-region

by

LF. McNamara
lonospheric Prediction Service
Sydney, Australia

OCTOBER 1978

Published by World Data Center A for Si,
BZ a a

Solar-Terrestrial Physics, NOAA, Boulder, Colorado
and printed by

U.S. DEPARTMENT OF COMMERCE WS
NATIONAL OCEANIC AND ATMOSPHERIC ADMINISTRATION \\) .
ENVIRONMENTAL DATA AND INFORMATION SERVICE
Asheville, North Carolina, USA 28801

SUBSCRIPTION PRICE: $25.20 a year; $17.30 additional for foreign mailing; single copy price varies."
Checks and money orders should be made payable to the Department of Commerce, NOAA/NGSDC.

Remittance and correspondence regarding subscriptions should be sent to the National Geophysical and
Solar—Terrestrial Data Center, NOAA, Boulder, CO 80303.

PRICE THIS ISSUE $1.29



CURVE YR MO OY
OBSERVED (0) 265 66
REFERENCE (R) O tt)
MODEL (M)

TIME SSNO LAT LONG METH MAG REF FIG SPC
120138 70 69.3 16.0 2 275 29
oo

1
(0) ie) 0.0 0.0 o oOo ie) t

Z |

ALTITUDE (KM)

a

ste |

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC |
OBSERVED (0) 226 663 224 34 69.3 16.0 2 275 03 |
REFERENCE (R} 0 0 0 0 0 0 0.0 0.0 0 00 0 |
MODEL (HM) |

ALTITUDE (KM)
a

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)
AURORAL ABSORPTION EVENT

34


OBSERVED

(0)

CURVE YR MO DY T.
12 700
oa

377

REFERENCE (R) O

MODEL (M)

ALTITUDE (KM)

OBSERVED

69
QO

1
Q

IME SSNO
105
io)

LAT LONG METH MAG REF

$8.8 265.8 1
0.0 60.0

217

Oo oo

FIG SPC

16.6

8
ie)

Va

ra

(0)

CURVE YR MO

381

69

REFERENCE (R) 923 69

MODEL (M)

ALTITUDE (KM)

o
ug
11.

ry

T
it)
ie)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

IME SSNO
105
105

LAT LONG METH MAG REF

$8.8 265.8 1
$8.8 265.8 1

217
117

FIG SPC

16.6
16.7

8
is)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

SOLAR PROTO!
22

IN EVENT



WORLD DATA CENTER A
National Academy of Sciences
2101 Constitution Avenue, N.W.
Washington, D.C., U.S.A., 20418

World Data Center A consists of the Coordination Office
and seven Subcenters:

World Data Center A
Coordination Office
National Academy of Sciences
2101 Constitution Avenue, N.W.
Washington, D.C., U.S.A., 20418
[Telephone: (202) 389-6478]

Glaciology: Rotation of the Earth:
World Data Center A: Glaciology World Data Center A: Rotation
[Snow and Ice] of the Earth
Inst. of Arctic & Alpine Research U.S. Naval Observatory
University of Colorado Washington, D.C., U.S.A. 20390
Boulder, Colorado, U.S.A. 80309 [Telephone: (202) 254-4023]

[Telephone: (303) 492-5171]
Solar-Terrestrial Physics (Solar and

Meteorology (and Nuclear Radiation) : Interplanetary Phenomena, Ionospheric
Phenomena, Flare-Associated Events,
World Data Center A: Meteorology Geomagnetic Variations, Magnetospheric
National Climatic Center and Interplanetary Magnetic Phenomena,
Federal Building Aurora, Cosmic Rays, Airglow):
Asheville, North Carolina, U.S.A. 28801
[Telephone: (704) 258-2850] World Data Center A
for Solar-Terrestrial Physics
Oceanography: Environmental Data and Information
Service, NOAA
World Data Center A: Oceanography Boulder, Colorado, U.S.A. 80303
National Oceanic and Atmospheric {[Telephone: (303) 499-1000, Ext. 6467]
Administration
Washington, D.C., U.S.A. 20235 Solid-Earth Geophysics (Seismology,
[Telephone: (202) 634-7249] Tsunamis, Gravimetry, Earth Tides,
Recent Movements of the Earth's
Rockets and Satellites: Crust, Magnetic Measurements,
Paleomagnetism and Archeomagnetism,
World Data Center A: Rockets and Volcanology, Geothermics):
Satellites
Goddard Space Flight Center World Data Center A
Code 601 for Solid-Earth Geophysics
Greenbelt, Maryland, U.S.A. 20771 Environmental Data and Information
[Telephone: (301) 982-6695] Service, NOAA

Boulder, Colorado, U.S.A. 80303
[Telephone: (303) 499-1000, Ext. 6521]

NOTES:

1. World Data Centers conduct international exchange of geophysical observations
in accordance with the principles set forth by the International Council of Scientific
Unions. WDC-A is established in the United States under the auspices of the National
Academy of Sciences.

2. Communications regarding data interchange matters in general and World Data
Center A as a whole should be addressed to: World Data Center A, Coordination Office
(see address above).

3. Inquiries and communications concerning data in specific disciplines should
be addressed to the appropriate subcenter listed above.



APPENDIX A (Continued)

c THE EQUATION OF TIME

EQT=~7,38*SIN(P1* (TF~4.))~9.87*SIN(P2* (TF+9.))+.27*SIN (P3* (TF-53.
2*COS (P4* (TF-17.))
QT*DTR/4.

FC=MER*DTR

PHO=FC~FB+ET

T=LH+LM/60.
PHI=,26179939* (T-12.) +PHO
A=SIN(FA) *SIN (DC)

B=COS (PA) *COS (DC)
COSK=A+B*COS (PHI)

RETURN

END

Listing of File PRODAT

246090850234409417419417414417944943809808971972947948541540542549543549 64 63
3243239 20000945946496495407405518517922 54548583579583580583581583221220450451
455451288287268229545544 48 49992991527 G158169159169357354358354938250933
932 6807 42 41282 6248 269 6272 9484 6485 G486 6487 6559 O560 @
488900480000481000482090483000562000 16006 17000 18090 19909 29000 21900
DELETE FIRST 6 LINES IF DOING MAG DISTURBED DAYS

523524275900356354 368000369000 3600003819233879233759 2537692537 7080378000379009
38 309038000038 4900 386090388000 389000 365000 36690038 2000 390000385000804000805000
89600935300036290026500015015114300014400H 86000 99094100000 60000 69090 76000
222600223090224900225000226000567 41168608 12 4125389625490099599499699397B009
967966968966969966

46



Oy TIME SSNO LAT LONG METH MAG REF FIG SPC
127200 27 69.3 16.0 2 274 369
ao Q 0.0 0.0 o 00 ie)

OBSERVED (0) 86 63
REFERENCE (R) 0 0
MODEL (M)

|
|
|

CURVE YR MO
9
ie)

ALTITUDE (KM)

« ||

a 3 4
LOS ELECTRON CONCENTRATION (N/CC)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC

OBSERVED (0) 99 643 1 15 69.3 16.0 2 274 3.9
REFERENCE (R} 0 O08 0 0 0 0 0.0 60.0 o oo ie}
MODEL (mM)

105

100 Fi

s

0

s

ALTITUDE (KM)
a 8

a 4
LOG ELECTRON CONCENTRATION (N/CC)
AURORAL ABSORPTION EVENT
30



CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 967 733 242115 44 65.3 212.2 1 21530 4.3 13
REFERENCE (R] 966 733 24 2100 44 65.3 212.2 1 2153 4.3 0
MODEL (TM) i

: a,

ALTITUDE (KM)

3 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC }
OBSERVED (0) 968 733 242135 44 65.3 212.2 1 2153 4.3 13 |
REFERENCE (R) 966 733 24 2100 44 65.3 212.2 1 2153 4.3 0 i]
MODEL (M)

ALTITUDE (KM)

° a

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

AURORAL ARC CONDITIONS. #966 PRE-ARC; #967 ARC FORMATION; #968 ARC BREAK UP
(SEE NEXT PAGE).

39



CURVE YR MQ OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 42 617 141100 53 45.4 284.1 106114 99
REFERENCE (R) 41 616 QO 1200 S6 45.4 284.1 10114 9 0
MODEL (M)

|
|

° x
. yA

7 |
, |
as a

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 484 622 18915 42 64.9 212.1 S 1104 6 9
REFERENCE (R} 0 O 0 00 tt) 0.0 «60.0 o 00 oO
MODEL (M)

ALTITUDE (KM)

a

\
Mm

ALTITUDE (KM)
a

0 ele

9 1 a 3 4
LOG ELECTRON CONCENTRATION (N/CC}
AURORAL ABSORPTION EVENT

25


OBSERVED (0)

MODEL (M)

ALTITUDE (KM)

OBSERVED

(0)

CURVE YR MO DY TIME SSNO LAT LONG METH MAG REF

945 69 11 18 1250
REFERENCE (RI 946 69 11 18 1430

105 45.4 284.1 1 198
105 45.4 284.1 1 198

FIG SPC
2 3
20

ZI

VA

406 703

REFERENCE (R) 40S 70 3

MODEL (M)

ALTITUDE (KM}

1
1

Te

2 3 4
LOG ELECTRON CONCENTRATION (N/CC)
CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF

10091
1003 1

06 45.4 284.1 1 198
06 45.4 284.1 1 198

FIG SPC
3.3
3.0

fal

4

3 4
LOG ELECTRON CONCENTRATION (N/CC)

SOLAI

R X-RAY EVENT
9



aaa

aa

aaa

aa

15

26

APPENDIX A (Continued)
Listing of Subroutine PROFILE for Magnetically Disturbed Days

SUBROUTINE PROFILE (IY,IM,ID,LH,LM,ISSNO,FLAT,WLON,DEN)

S/R TO CALCULATE A D~REGION PROFILE USING THE STATISTICAL

WITH COEFFICIENTS A,B,C,D,E,F
MAGNETICALLY DISTURBED DAYS ONLY

DIMENSION A(8) ,B(8) ,C(8) ,D(8) ,E(8) -F(8) »T(8) ,DEN(8) ,G(8)
A = CONSTANT TERM 99(5)55 KM

DATA A / 3.12,2.36,1.88,1.66,1.49,1.25,0.81,0.82 /

B = COS KI TERM

DATA B / .82,.86,.98,.95,+84,.58,.68,.77 /

C = GEOMAGNETIC LATITUDE TERM

DATA C / 05,024,215, 209,018,257, +78,0+ /

D = SUNSPOT TERM

DATA D / .15,.15,.22,.26,+13,.26,.25,0. /

E = SEASONAL TERM

DATA E / .17,.18,.18,0.,~0.06,0+,412,0. /

MAGNETIC TERM

DATA F / 017,634,631, +32, +32,-21,+45,0./

RESIDUAL STANDARD DEVIATION

DATA G / 153,251,247, 240,233, +32, 240,047 /

PI=3.141592654

PREPARE TO CHANGE TO GEOMAGNETIC LATITUDE AT WLON DEGREES (WEST)

1.5*DEGRAD
COSGM=COS (THETGM)
SINGM=SIN (THETGM)
PHIGM=291.*DEGRAD

CONVERT TO GEOMAGNETIC LATITUDE

COLAT=(99.~FLAT) *DEGRAD

LON* DEGRAD

‘OS (COLAT) *COSGM+SIN (COLAT) *SINGM*COS (ALON-PHIGM)
GEOML=99 .~ACOS (COSA) /DEGRAD

TLAT=SIN (2.*GEOML*PI/180.)

214.9% (SIN (GEOML*PI/186.)) **4

M

FMON~6.5) /12.*2.*PI

0S (FMON)

IF(FLAT.LT.@.) PMON=~FMON

SEASON=FMON

IID=ID

IF(ID.EQ.8) TID=15

CALL COSKI(IY,IM,IID,LH,LM,FLAT,WLON, COSX)
DO 29 I=1,8

IHT=99~(I~1) *5

XLAT=TLAT

IF(I.GE.5) XLAT=ALAT
LAT*C (I)
SSNO*.G1*D(I)
EASON*E (I)

ED=ED+T (J)
DEN (I)=ED
CONTINUE
RETURN

END

47

MODEL,



APPENDIX B (Continued)

CALL INTNO(JSUN,2.6,"ABUT")

CALL REALNO(SJTLAT,1,3.25,"ABUT")

CALL REALNO(SJTLONG,1,4.,"ABUT")

CALL INTNO(JMETH,5.2,"ABUT") |
CALL INTNO(JMAG,5.65,"ABUT") I
CALL INTNO(JJREF,5.85,"ABUT")
CALL MESSAG (UJFIG,5,6.5,"ABUT")
CALL INTNO(JSC,7.25,"ABUT") i

CALL MESSAG ("MODEL (M) ",14,-1.5,6.5)
20 CALL ENDPL(IEND)
c IEND=IEND+1

GOTO 2

40 CALL ENDPL (IEND)
CALL DONEPL
END
SUBROUTINE CURVENO(IDIST,IYR,IMO,IDY,ITIME,ISUN,STLAT,
XSTLONG , IMETH, IMAG, ICOL, [REF ,UFIG,1ISC,JREF,
XJDIST,JYR,JMO,JDY,JTIME,JSUN,SJTLAT,
XSJTLONG , JMETH , JMAG ,JCOL,JJREF ,UJFIG , JSC, JJJREF ,N,NH, MODE ,10,1C)
COMMON A(2,6) ,X(100) ,¥ (100) ,INPT(1)
c PRINT 103,NH
193 FORMAT(I3)
IF(NH.GE.1) GOTO 3
1 READ(3,199) IDIST,1YR,IMO,IDY,ITIME,ISUN, STLAT, STLONG,
XIMETH, IMAG, COL, LREF ,UFIG, ISC ,JREF
IF(EOF(3)) 2,308 |
300 READ(3,190) JDIST,JYR,JMO,JDY,JTIME ,JSUN, SJTLAT,SJTLONG,
XIMETH, JMAG ,JCOL, JOREF ,UJFIG ,JSC ,JJJREF
IF(EOF(3)) 2,3
192 FORMAT(I3,13)
2 N=2
RETURN
3 NH=NH+1
IF(IG.EQ.1) 17,4 |
4 READ (3,191) IMODEL,ICURV i
c PRINT 191,IMODEL, ICURV
IF(EOP(3)) 2,5 F
5 10=1 f
IC=IC+1
IF(ICURV.EQ.IDIST) 6,7

6 IF(IMODEL.EQ.1) 8,9 |

7 IF(ICURV.EQ.JREF) 12,17

8 MODE=-1 |
GOTO 17 |

9 MODE=1 |
GoTo 17

12. MODE==2

17 RETURN

108 FORMAT(2X,13,313,15,14,2F6.1,212,13,14,A5,13,13)
191  FORMAT(I1,13)
END
SUBROUTINE EPOINTS(N)
COMMON A(2,6) ,X(10B) ,¥ (190) ,INPT(1)
READ (3,160) ((A(I,J) ,1=1,2) ,J=1,6)
IF(EOF(3)) 5,10
N=2
vi) RETURN
90 FORMAT(2X,6(P4.2,F6.2,2X) )
END
SUBROUTINE FILLUP(M,KPT,N)
COMMON A(2,6) ,X(16) ,¥ (190) ,INPT(1)
18 DO 15 K=1,6
M=M+1
X(M)=A(1,K)
¥(M) =A(2,K)
IF(X(M) .EQ.0.89) 12,14
12 RETURN
14 KPT=KPT+1

BREW

49

|


OBSERVED (01 548 68125 1200 110 57.2 352.7 2 132 105
REFERENCE (R) 583 71 10 14 1255 66 57.4 352.6 2 °=133 4 0
MODEL (M)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC |
|

° BP |
a) La |
7 |

ALTITUDE (KM)

a 3 4
LOG ELECTRON CONCENTRATION (N/CC)

CURVE YR MO OY TIME SSNO LAT LONG METH MAG REF FIG SPC
OBSERVED (0) 450 71 1 221200 80 32.3 253.5 2 196 7 5
REFERENCE (R) 451 712 1 1215 78 32.3 253.5 2 196 7 0
MODEL (M4)

: ‘

4
LOG ELECTRON CONCENTRATION (N/CC)
WINTER ANOMALY DAY

12

ALTITUDE (KM)

t
}
i
|
1
|
|
{