Wmo-Td 149
Wmo-Td 149
l
Claus Fr6hlich & Julius London
Editors
REFERENCES 45
4. RADIATION MEASUREMENTS 49
4.2.2 Pyranometers 62
4.2.2.1 Most Widely Used Pyranometers 63
4.2.2.2 Installation of Pyranometers 69
4.2.2.3 Maintenance of Pyranometers 69
4.2.2.4 Calibration of Pyranometers 70
4.2.2.5 Measurement of Diffuse Radiation 71
4.2.2.6 Measurement of Solar Radiation
on Inclined Surfaces 73
4.2.2.7 Special Pyranometers 73
REFERENCES 88
TABLE OF CONTENTS
5.1 INTRODUCTION 93
REFERENCES 114
TABLE OF CONTENTS
INDEX 137
FOREWORD
1
1. INTRODUCTION
Before the advent of rocket and satellite observations it was only pos-
sible to infer the mean irradiance at the top of the atmosphere (the solar
constant) by extrapolation of measurements from ground-based, balloonborne or
aircraft mounted instruments. However, it is now feasible to observe directly
the total solar irradiance and its spectral distribution from platforms above
3
CHAPTER 1 INTRODUCTION
the earth's atmosphere. Observations reported during the past five years
indicate an average solar constant of about 1367 Wm- 2 with a standard devia-
tion of the observed values of about ±0.3 per cent (Willson, 1984). Approxi-
mately 95 per cent of this energy is in the visible and near-infrared spec-
tral range (0.325 tlm - 2.8 tlm). About 30 per cent of the solar irradiance
received on a horizontal surface at the top of the atmosphere is reflected
back to space, primarily by clouds, but also by the clear atmosphere and the
earth's surface. Slightly more than 20 per cent is absorbed by atmospheric
water vapour, carbon dioxide, ozone and by clouds. ·The remaining incoming
solar energy, just under 50 per cent, is received at the earth's surface and
utilized to evaporate large quantities of ocean-surface water and to heat the
land and ocean surfaces. A small but essential amount of this energy is
converted for use in sustenance and growth of the world's biological system.
The earth radiates in the infrared (4 tlm - 100 tlm), but most of this
upward radiation is absorbed by clouds, atmospheric water vapour and carbon
dioxide. As a result, less than 10 per cent of the earth's emitted radiation
directly escapes to space. The remaining outgoing radiative energy emanates
from cloud tops and atmospheric· gases. On an average, downward radiation,
chiefly from the base of clouds, and from water vapour and carbon dioxide in
the lower atmosphere almost balances the upward surface radiation. Details of
the radiation characteristics of the earth and its atmosphere are discussed
in Chapters 2 and 3 of this Manual.
4
CHAPTER 1 INTRODUCTION
the net infrared-radiation balance at the ground has led to some effort to-
wards improving the surface network of infrared observations. Such measure-
ments are also required for detection of long-term changes of upward- and
downward-directed radiative components at the earth's surface to provide
useful information of the mechanisms of climate change. Since there are no
methods other than ground-based observations to measure the downward infrared
radiation at the surf ace, it has been proposed by the Radiation Com mission
that the present, very limited network of net-atmospheric-radiation measure-
ments be extended. Detailed descriptions of radiation instruments, both solar
and terrestrial, that are presently operational, including principles of
measurements and techniques for their use, are described in Chapter 4.
5
CHAPTER 1 INTRODUCTION
REFERENCES
Coulson, K.L. (1975) Solar and Terrestrial Radiation, Academic Press, New
York.
CSAGI (1958) Radiation Instruments and Measurements, Part VI, IGY Instruction
Manual, p. 371. Pergamon Press; London, New York, Paris.
Kendall, J .M., Sr. and C.M. Berdahl (1970) Two blackbody radiometers of high
accuracy. App.l. opt. 9: 1082.
Moller, F. (1980) A Brief History of the Radiation Commissions of the IMO and
of the IUGG, Radiation Com mission (lAMAP) p. 43.
Raschke, E. and K. Ya. Kondratyev: Radiation Budget of the Earth and Its
Atmosphere, WCP-40, World Climate Research Programme.
Stephens, G.L., C.G. Campbell and T.H. Vonder Haar (1981) Earth radiation
budgets measurements from satellites and their interpretation for cli-
mate modeling and studies. J. G'eopl1ys. Res. 86: 9739.
Willson, R.C. (1984) Measurements of Solar Total Irradiance and its Variabi-
lity. Space Sel. Rev. 38: 203.
6
2. RADIATION IN THE ATMOSPHERE
2.1 INTRODUCTION
The interaction of radiation with the atmosphere and the underlying sur-
face, such as land, water and ice, involves:
The most important gaseous atmospheric constituents of dry air are lis-
ted· in Table 2.1. Water vapour (H 2 0) has a molecular weight of 18.02 and is
present in variable amounts of 0-0.04 fraction of total molecules.
Table 2.1 Composition of the earth's atmosphere below 100 km (modified after
U.S. Standard Atmosphere, 1976).
Molecular Content
Constituent weight (fraction of total molecules)
These gases are uniformly mixed, except as noted below, up to the meso-
pause (Le. to about 85 km), where photo-dissociation and gravitational dif-
fusive separation becomes increasingly important with increasing height. Car-
bon dioxide has a seasonal variation in the lower troposphere of a few parts
7
CHAPTER 2 RADIATION IN THE ATMOSPHERE
per million, and ozone, abundant mainly in the lower and middle stratosphere
(15 - 35 km) also has large seasonal and geographic variations. Atmospheric
water is found in the atmosphere in all three states: vapour, fluid and
solid.
2.2.2 Aerosols
Soil particles, dust and sea salt form aerosols, which are radiatively
active components of the atmosphere. They are active in a direct way by:
8
10
6 MARITIME
...
'0
-10
AEROSOL MODE L
--"...
Z
10
4
2
CONTINENTAL /
-
'0
10 URIGINN
~
0 ~~, \
0
c 10
CD / / SEA SPRAY \
'0 -2 / PRODUCED \
... 10 I
CD
.0 I
/ \
E 10 -4 I
\
::::J
Z \
-6
,10 \
\
-8
10 -J -2 -1: 0 1
10 10 10 10 10
Radius <lAm)
8
CHAPTER 2 RADIATION IN THR ATMOSPHERE
change in their size distribution and hence their radiative properties. The
largest aerosol concentrations are found in the planetary boundary layer over
industrial areas and over the oceans. After strong volcanic eruptions a
different, much smaller, aerosol maximum may be found in the lower strato-
sphere (Fig. 2.2).
70
60 STANDARD AEROSOL MODEL
~ VERTICAL PROFILE S
50 ~
40 - ' HAZY
---CLEAR
30 "\
CHANGE \ ....... RAYLEIGH
IN \
-E
JiIi:
..,....
25
SCALE \
\
\
\
\ .
'
RAYLE IGH
SCATTERING
CD 20 BACKGROUND \ . . . VOLCANIC
"C STRATOSPHERIC I '.
.....-
-
~
".\I
15
et I
I w
I 0:
UJ
\
10 \
I
Cl..
,,
\ I./)
\ 0
Cl..
'. 0
5 \. 0:>-
~o:
\
'.',
'.
.....
oc{
0
Z
-::>
-4 -3 -2 -1 0 0
cD
10 10 10 10 10
Radius (pm)
Figure 2.2 Vertical profile of the extinction coefficients at 550 nm of
a tmospheric aerosol profiles (Shet tle and Fenn, 1976).
West
150 120 90 60 30 0
o
('I')
-
.c
: :l
o
(f)
9
CHAPTER 2 RADIATION IN THE ATMOSPHERE
2.2.3 Clouds
Clouds are found in all tropospheric layers and, occasionally, in the
lower stratosphere. Clouds contain water particles in one or two states:
liquid and/or solid. Although estimates of the global cloud distribution vary
somewhat, it is likely that the global cloud cover is around 60 per cent
(e.g. Berlyand et al., 1980). Fig. 2.3 (from National Academy of Sciences,
1975) shows the global distribution of fractional cloud cover for the month
of January.
The usual size of cloud droplets and particles ranges from 2 to 20 tim.
Larger raindrops may have diameters of up to 10 mm. The largest "particles"
are found in convective tropical and extra tropical cloud systems as hail-
stones with diameters of up to 10 cm. Typical size distributions are shown in
Fig. 2.4. Fig. 2.5 shows a vertical profile of the size range of droplets in
a low stratocumulus cloud. This figure also shows a profile of the liquid-
water content (LWC in g of liquid water per m- J air); on occasion the LWC may
reach even larger values. The content of solid water in dense cirrus may be
as high as 0.03 gm- J (Griffith and Cox, 1977). .
2
n(r) 10
1
10
0
10
-1
10
-2
10 0 2· 4 6 10 12 14
r(pm)
Figure 2.4 Size distribution of cloud droplets (Deirmendjian, 1969).
8.AUG.1978
Sc -3 8.Aug.197B_ 1000
LWC=0.4gm 3
Sc-1 cloud to~level
Sc-2 mid cloud leve I E800
Sc-3 cloud base level _ -a
Gi 600 ./
:c adiabatic LWC
Clouds primarily reflect and scatter incident solar radiation. They are
almost black for terrestrial thermal radiation and absorb practically all
incident IR radiation. The emission is thus determined mainly by the tempe-
rature at their surface.
10
CHAPTER 2 RADIATION IN THE ATMOSPHERE
~=o
!p=
1 dL
a. = L --
ds
(m- 1 )
'
(2.2)
11
CHAPTER 2 RADIATION IN THE ATMOSPHERE
particle density N(m- 3 ), the single particle cross section (m 2 ) is given by:'
= .GE'A'S.
(2.3)
N
For a single scattering process within the volume element, the single scat-
tering albedo w is defined as the ratio of the part scattered to the total
extinction:
(2.4)
The "optical thickness", OE, between the points Xi and X2 is defined as:
X2
OE ::: JaE' ds (2.5)
Xi
where 12 ::: -1, and the subscripts rand i mean the real and imaginary
parts respectively. The reflection and refraction processes (including scat-
tering) are due to the real part m.. , while the imaginary part m\ stands
for absorption ("damping") of radiation. In general, the index of refraction
depends on the radiation wavelength.
Radiation is scattered and absorbed along its path through the atmo-
sphere - or, in more general terms - through any gaseous, fluid or solid
matter. Also, radiation may be added to a given beam, either due to scat-
tering into the direction of the beam from any other direction or due to
thermal emission within the volume element. Thus the transfer of a beam of
radiation through a volume element can be given' by a differential equation
which describes the energy budget along the path of the beam through this
volume element.
12
CHAPTER 2 RADIATION IN THE ATMOSPHERE
(2.7)
The first part of the right-hand side of Eq. 2.7 is the law of Bouguer -
Lambert as described above. The second part, J(j1,ijI), is called the source
function and describes the direct and scattered solar radiation and thermal
radiation emitted from and scattered into and out of this volume element.
Each contribution of the source function depends on the spectral range consi-
dered and on the properties of the atmosphere.
dt
(2.8)
j1
A similar solution exists for the downward radiance. The first term in
Eq. 2.8 describes the radiance transmitted through the layer. The second
term describes the light scattered and radiation emitted from the layer in
the direction (j1, ijI) and integrated from 00 to 0 where t represents the dummy
integration variable. Many algorithms have been developed to solve this equa-
tion numerically for a great number of practical applications to realistic
atmospheric conditions. A summary of the algorithms and a comparison of
their performance has been published by the Radiation Com mission of IAMAP
(Lenoble, 1977 and Fouquart et aI., 1980). In terrestrial applications radia-
tion·from the sun dominates the source term at wavelengths less than about
4 j1m, whereas thermal radiation is dominant for A > 4 j1m.
The term 'scattering' is. used to explain all processes occurring within
geometrical. dimensions of the order of, or smaller than, the wavelength of
the radiation involved. The spherical waves reflected in all directions to-
gether with the refracted part affect scattered light through subtle inter-
ference.
Cs • Lo
L(8,r) = 4nr 2 • p(8,r) , (2.9)
13
CHAPTER 2 RADIATION IN THR ATMOSPHRRE
with N as the total number of particles per unit volume. Similar formulas
have to be used for the phase function, the single scattering albedo, and
other optical properties.
32 1
aSR =3 N • (m-1)2. (2.11)
where m is the refractive index of medium (e.g. air), and l)p is the so-
called depolarization factor. For air l)p = 0.0279 (Young, 1980). The phase
function p(8 s ) is given by:
p (8 s ) - -4
3
. (1 + cos2 8s ) , (2.12)
where 8s is the scattering angle, Le. the angle between the direction of the
incident beam and the scattered ray (8 s = 0 for forward scattering, 8s = 180
for backward scattering). .
o
It) 0 = 1.50 + O.OOi
o
.
o·
('I)
o
C'I
o
X =604 om
o
o-1-- ....-:~~~;d::::::--_.__--,--~~--,--,..---.,~-__r___'-...__~
0.01 ·0.02 0.05 0.1 0.2 0.5 1.0 2.0 5.0 10.0
Radius (pm)
N = 100cm'3
1.0
.6
...,.._.~~. / I
/
Lt ,\
. C3
.......
E
I
--...
CJ
10
2
\~
....>-
.2 III
c 10"
Gl
.. .1
..
'C
GI
~P3
\') .0
E 10 5
.06 L :::l
Z
.04
\ 10&
.02 ~
.01
H- N p4
_ 10.7
.4 .8 1 2 4 .01 .1 1 10
A (.-n) Radius (IJm)
-.o...
() 2
10
-Sc-1
C ---Sc-3
....::I 550nm
,,
-2
10
-3
1 0 ....T"'T",.,-..,....,..,....",.....'T""I"........TT'"r..-r---y----i
2 6 10 30 50 70 90 110130150170175 180
Scattering angle
Figure 2.9 Phase functions for the drop-size distributions Sc.-1 and Sc-3 in
Fig. 2.5.
15
CHAPTER 2 RADIATION IN THE ATMOSPHERE
There are three principal processes in the natural environment that are
responsible for polarization. First, very strong polarization fields are
16
CHAPTER 2 RADIATION IN THE ATMOSPHERE
,
;'
a: .... \ \ \ \. 0.50
I
w
'~0.75
Cl .
w 1&.1
0
20
~\
A-O ~\
\ \ \. \ ~ 20
1&.1
0
,\,\ \
\
10 0.25~ \ 10 '\
0.50.:=::--", \. '\
-.;:~
0.75 "
0 0
-10 -10
17
CHAPTER 2 RADIATION IN THE ATMOSPHERE
The maxima of P are considerably lower here than for 0.80 Mm, and the
negative portion of the polarization field is much more pronounced. The
neutral points are distinct in both cases, their angular distance from the
sun is greater, and the Brewster point of the previous diagram gives way to
the Arago point here. All of these features are explained by the greater
amount of multiple scattering which occurs at 0.365 Mm.
(2.13)
(2.14)
18
CHAPTER 2 RADIATION IN THE ATMOSPHERE
REFERENCES
Berlyand, T.G., L.A. Stroknia and L.E. Gresnnikova (1980) Zonal Distribution
of Cloud Amount on the Earth. Het.lwd G.idr. 3: 15, Allerton Press.
Coulson, K.L., J. V. Da ve and Z. Sek era (1960) Tables Rela ted to Radia tion
Emerging from a Planetary Atmosphere with Rayleigh Scattering.' Univer-
sity of California Press, Berkeley and Los Angeles
Coulson, K.L. (1974) The polarization of light in the environment, pp. 444-
471, . in Gehrels, T. (Ed.) Planets stars and Nebulae Studied with Photo-
polarimetry, University of Arizona Press, Tucson, Arizona.
Fouquart, Y., W.M. Irvine and J. Lenoble ed. (980) Standard Procedures to
Corn pute Atmospheric Radia ti ve Tra n sfer in a Sca t tering Atmosphere.
Vol.II. Radiation Corn mission, CIAMAP), Boulder, Colorado.
Griffith, K.T. and S.K. Cox (1977) Infrared radiative properties of tropical
cirrus clouds inferred from broadband measurements. Atmos. ScL Paper
No. 269, Colorado State University, Fort Collins, (NTIS PB 26853).
Kneizys, F.X., E.P. Shettle, W.O. Gallery, J.H. Chetwynd, Jr., L.W. Abreu,
J.E.A. Selby, R.W. Fenn and R.A. McClatchey (1980) Atmospheric Transmit-
tance/Radiance: Computer Code LOWTRAN 5, AFGL - TR - 80 - 0067, Environ~
mental Research Papers, No. 697, Optical Physics Division, Air Force
Geophysics Laboratory, Hanscom, Massachusetts.
Mie, G. (1908) Bei tdige' zur Optik tri..iber Medien, speziell kolloidaler
Metalllosungen. Anna.len der Phys.ik 25: 377.
19
CHAPTER 2 RADIATION IN THE ATMOSPHERE
Shettle, E.P. and R.W. Fenn (1976) Models of Atmospheric Aerosols and their
Optical Properties. In: AGARD Conference Proceedings No. 183. Opt.ical
Propagation in the Atmosphere. Presented at the Electromagnetic Wave
Propagation Panel Symposium, Lyngby, Denmark, 27-31 October 1975, AGARD-
CP -18 3 , available from U•S • Na t i on a I T e ch n ic a I Inform a t ion Se r vie,e
(No. AD-A028-615).
3.1 INTRODUCTION
The radiant energy output of the sun is 3.94-10 26 W, of which 1.8 . 1017
W is intercepted by the disk of the earth. R6ughly 30 per cent of the inter-
cepted energy is scattered back to space, but the remainder is absorbed and
tends to heat the earth-atmosphere system. Since in fact the annual mean
temp~rature of the planet is fairly constant, the absorption of solar energy
must be bal~nced by an equal energy output. The planet has adopted a tempe-
rature and temperature distribution such that its average output of thermal-
infrared radiant-energy equals the solar input.
100 6 21 ,4 6 37 26
ATMOSPHERE
Net Emission by ~
Absorbed Water Vapour, CO 2 Emission
18
by Water by Clouds
Vapour, Dust, O.
Net Gain by Absorption
Absorbed 10 by Water Vapour, CO 2
by Clouds
Latent
Heat Flux
Net Surface Sensible'
Emission of H~at Flux
Long-wave Radiation
OCEAN,LAND 48 16 10. 22
Fig. 3.1 gives some detail of what happens to the energy stream on its
way into and out of the earth-atmosphere system. Note that only 48 per cent
of the solar radiation impinging on the top of the atmosphere actually
reaches the bottom, and it is this 48 per cent which is the primary drive for
all processes on or involving the earth's surface. Most of the physical sur-
face processes (evaporation of water, surface heating and so on) involve the
total spectrum of solar radiation which, at the ground, covers the wavelength
band from 0.29 to about 2.5 ,urn. Biological processes, on the other hand,
often depend on the amount of solar radiation in specific wavelength regions
of the total spectrum. For instance, sunburning of skin, DNA breakdown, and
the degradation of paints are all primarily dependent on the shortest-wave-
length (ultra-violet) portion of solar radiation, and photosynthesis involves
mainly the visible part of the spectrum (see Fig. 3.2). The surface values of
solar radiation in different wavelength bands are complicated- by the fact
that absorption and scatter of solar radiation by the atmosphere itself is
highly wavelength-dependent.
21
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
;4
:', n
11
"
Cl) "
I,
0CO ""
C 0. I'
It
0 1\
:
Q. , \\
UiU) : \ 2
Cl) • I \;
~O I
I \
\
./ \
" \
Q)'lit
,
I
> . I I
_.
I
;;0 I I
...
I
(Q.. I I
1
, I
I
G)N \ ; I
\
, f I
o \ .' I
,
.' \ I
\ ' \
" \,
O~ __ ~..,)".. ......._~....-_ ...... ....... ..- ...........--~_
o 300 400 5.00 600 700
Wavelength (nm)
Figure 3.2 Relative "response spectra" to solar radiation for various bio-'
logical processes. Curve 1 = photosynthesis; curve 2 = erythema; curve 3
= product-ion of vitamin D; curve 4 = bactericidal effects.
The infrared radiative output from the surface and from objects on the
surface is also governed by the state of the atmosphere and is also spec-
trally dependent.
The flux density (solar irradiance) impinging on a unit area at the top
of the atmoshere directly facing the sun, ~, depends on the sun - earth dis-
tance R and thus varies significantly throughout the year. The solar con-
stant, S, is then defined as So when the earth is at mean distance from the
sun; that is, when Ra = 149.6 X 108 km. The average value of the solar con-
stant determined from balloon, rocket and satellite measurements over the
period 1976 - 1981 was 1367 Wm- 2 with an uncertainty of about ±0.3 per cent
(Frohlich and Brusa, 1982 and World Climate Programme, 1984). Since the· solar
irradiance ~ is proportional to the inverse square of R and the day-to-day
variation of R is known very accurately, the daily value of So can easily be
computed. The following formula gives an empirical determination of So accu-
rate to 0.01 per cent. (Spencer, 1971),
2:2
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
R~
S = So . R2 = so' (1.000110 + 0.034221·cos 8 0
( 3.1 )
+ 0.001280 • sin 80 + 0.000719 • cos 28 0 + 0.000077 . sin 28 0 )
where Ra is the mean sun - earth distance and 8 0 is defined in terms of the
Julian day of the year - that is, of day d counted from 1 to 365 or 366,
starting on January 1st. Thus,
2nd
= 365
(3.2)
S presently ranges from So +3.3 per cent on about January 3rd (perihelion) to
So -3.3 per cent on about July 5th (aphelion).
For most purposes the value of S can be regarded as constant with time,
although satellite measurements have revealed that it varies at the level of
a few tenths of a per cent on time scales of days to weeks (see Fig. 3.3). It
may be that there are similar variations over the period of a solar cycle,
but these have not yet been verified.
0.15
MEAN
0.05
1367.7
-0.05 %
Wm- 2
0.15
,--"-~--'--.L..-..L-"""""""'---'---'---!.-~-'--L..-...J-....L."...,J,~~~~~--,L-J---'--.L..-..L--'----'~....I.....&-O.
2 5
100 200 300 100 200 300 100
1980 DAY 1981 DAY '1982 DAY
Figure 3.3 Solar "Constant" as measured with ACRIM on the Solar Maxim urn
Sa telli te from 1980 to 1983 (Willson, 1984).
where again 8 0 is in terms of Julian day according to Eq. 3.2. The hour
angle (0 to 2n:) is counted from midday of local apparent or true solar time
23
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
LAT and can be obtained from the clock time of a given place from the rela-
tion:
Here clock time is the time referred to any of the standard meridians
about the world; the longitude correction is 4 minutes for every degree of
the given place east (positive) or west (negative) of the standard meridian;
and Eq is the equation of time for which an empirical relation is as follows
(Spencer, 1971):
This quantity mR is known as the relative air mass, and ranges from 1.0
for an overhead sun, through 2.0 at a zenith angle of 60·, to infinity with
the sun on the horizon. This simple formula is valid for a plane parallel
atmosphere. Taking into account a spherical earth atmosphere with refraction,
mR never really reaches infinity, but is limited to about 36.5 as given by
Bemporad's formula (see, for instance, List, 1968). If more accurate values
are needed, e.g. for sunphotometry at high air masses, the formula of Kasten
(1966) may be used:
1
(3.9)
sin-8-+-0~15(8+3~885)~~-
The solar input at the top of the atmosphere is entirely direct radia-
tion so that the instantaneous input per unit area at the top of the atmo-
s p her e is si mply So' cos 8, wher e 8 is t he I 0 c a Iso I a r zen i t h an g I e a t a
particular time. Calculation of the total daily input involves integration of
So . cos 8 over a 24-hour period. Calculations of this sort for all locations
and all times of year lead to the contour plots of daily input as a function
of latifude and time given in Fig. 3.4. It is notable that the maximum daily
total input occurs at the poles in midsummer. This is because there. is 24
hours of daylight per day at the poles during that time. Note also that, at
present, the solar irradiance during the summer is slightly higher at polar
latitudes in the Southern Hemisphere than in the Northern Hemisphere. This is
because the Southern-Hemisphere summer solstice currently occurs at a time
close to perihelion.
24
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
MJ m- 2
EQ.
J F A M J J A s o
Figure 3.4 Contours of daily total input of solar radiation on a horizontal
surface at the top of the atmosphere as a function of latitude and time
of year.
1I'l...- ...,
-
N
E
C
....... 1I'l
(lE ...:
.......
::
-
25
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
There are three major mechanisms which modify solar radiation on its way
down through the non-cloudy atmosphere.
Ozone exists at all atmospheric levels, but has its maximum concentra-
tion at about 20 - 25 km. It absorbs solar ultraviolet and visible radiation.
The absorption is concentrated in three wavelength bands, the Hartley and
Huggins bands in the ultra-violet (UV) up to about 0.36 tim, and the Chappuis
band from about 0.42 tim to 0.74 Jim in the visible region of the spectrum.
Within each band the absorption is a strong function of wavelength. The UV
absorption by ozone ensures that no radiation of wavelength shorter than
0.29 tim reaches the ground. The absorptance a of each band (the fraction
absorbed of the total spectrum of solar radiation impinging on the top of the
atmosphere) can be expressed to fair accuracy as an empirical function of the
ozone path, x, (in cm STP) traversed by the solar beam (Lacis and Hanson,
1974). Thus
a 0>
YIS (3.10)
and
1082x 0.0658x
a UY
0>
-- ---------------- + -------------
(1 + 1 38.. 6 x ) 0.8Il5
(3.11)
1 + (103.6xP
Since the ultra-violet and visible bands are spectrally separate, the
total absorptance is simply the sum of the absorptance in both regions. Fur-
ther, it is a fairly good approximation to assume that most of the ozone
absorption takes place in the high atmosphere - that is before scattering
processes have modified the direct nature of the solar irradiance.
26
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
S = So (1 - aw v - ao Z - OR) . (3.14)
(3.15)
S
-
So
<3.16)
27
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
Lt')
,
C\I
Extraterrestrial
,..
I
-""'
C\I
,
E
c
N
I
E Lt')
,
;: ~ Elevation
H2 O
(1) 0,
(J ~
CO 2
C
<0
.- - 03
"'0
<0 - O2
-
"-
"-
Lt')
0
,
o,
o 300 500 1000 2500
A (nm)
where A is the wavelength, and ~R,"'6(A) the optical depth for Rayleigh
scat tering (R), aerosol extinction (A) and gaseous absorption (G) respective-
ly. The Rayleigh optical depth can be calculated for the total atmosphere
from Eq. 2.11 and can be expressed by the following approximation (Frohlich
and Shaw, 1980 and Young, 1980).
28
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
ble (for the technique of these measurements see 4.2.4.1). Angstrom assumed
the following behaviour of ~A:
which means that /3, the Angstrom coefficient, is the aerosol optical depth at
a wavelength of 1 /-lm and the extinction follows a power law as shown in
Fig. 2.8 in Chapter 2. SchUepp replaced the exponent 1.3 by parameter a and
changed from natural to decadic logarithmS 5. Furthermore, his coefficient is
referenced to the optical depth at 500 nm:
<3.20)
SchUepp's B is only valid for wavelengths below about 850 nm, whereas
Angstrom's coefficient has no wavelength limitations. From Eq. 3.18 and 3.19
it is easy to calculate /3 from B and a and vice versa:
/3 = B·(ln10-0.51.3) = 0.935-B
B = /3/ (l n 10·2-«) = 0.434 -/3 - 2« (3.21)
3. 5 • 3 Solar Aureole
The downward diffuse flux density for a dust-free atmosphere over a non-
reflecting ground can be written as:
29
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
~~~~o~... ~
~ ~
Figure 3.7 Distribution of the sky radiance at 368 nm over the hemisphere
from measurements in Switzerland at an elevation of 400 m a.s.l. norma-
lized to the zenith radiance (Heimo, 1985).
In the real situation, how'ever, one must also take into account that the
direct-beam radiation, which is initially reflected by the ground, is, in
turn, scattered back to the ground by Rayleigh scatter so as to augment the
diffuse flux. Fig. 3.8 demonstrates that the actual total flux reflected up
from or down from such a multiple reflection system is (in either case) the
sum of an infinite series. If Et and Et are the normalized flux densities
upwelling and downwelling from the layer, extension of Fig. 3.8 shows that
Et (1-p)2 P6 Et (1 - p)
= P + --------- and = --------- (3.26)
E co (1 - PP6 ) E co (1 - PP6 )
p (1-p)2p (1-p)2p 2p
o o
( 1-p)
30
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
These relations yield the flux densities in Wm- 2 • They are derived from
data of a typical mid-latitude station with average turbidity coefficient f3
equal to about 0.04, and could be in error for an instantaneous individual
value by as much as 20 per cent. The spectral distribution for clear sky
conditions are shown in Fig. 3.9.
Cl)
(J
C
.-ca
"......
2.5
ca
--...ca 2.0
...
(J 15
Cl)
Q.
un
en 10
Cl)
...ca
>
-a::
(1)
500 600 700 om
Wavelength A
Figure 3.9 Spectra of sun and clear sky (E g {-) and of northern sky (E d )
(i.e. diffuse) radiation (each normalized to the irradiance at 550 nm).
31
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
features are more pronounced here than in cases of less extreme aerosol con-
centrations, they are indicative of the effects produced by aerosols in gene-
ral and demonstrate the diagnostic potential of polarization measurements.
8o,..-------'------------,
70 !-
...-...•
•
~ :
•• ••
•• o.70~m
~ 60-
z
o
•
•
• .............•
....
....
'
i=
< 50 ~ •
!:::! •
a:
<
..J 40~
•
• .. ..
~
&I-
o 30 .... :
'.'.
~
w
w •
a:
<:)
w
•••
•
o
10 - '.•'.
•, 9
..
,.. ..
••
.'
••
.~ of .'
ol--------'---.,---~~.-~*---t
BA/"" ""'BR
I I I I I I I. I I
-10 L-aJ..0--6.LO--4.LO--2..J.O--I.0---:-L20:-----:J40'::-""'""'=6'::'0--:S"=0.....
ZENITH ANGLE le)
60..--""'T"'--.,.....----,..-----..---r---,.----,..---"'-,---""T"""--t
50
;z40
o
.-
«
~O
0:
«
...I
o
Q.20
I&.
o
~IO
a:
(!)
w
o 0
-10---=8~0----=6'-::-0----J4..."..P-........."2L,,O--0~--:2~O~-4,....,O:----,6:"-:O~-8~0~.....
ZENITH ANGLE (.)
32
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
(3.30)
where a is the Stefan - Boltzmann constant (5.67 . iO- 8 Wm- 2 K-4)· and T is
the .absolute temperature (K). The actual distribution of this radiation as a
function of wavelength is also a function ·0£ temperature (see Fig. 3.12), If
a body absorbs only a fraction f; of thermal radiation (and f; being the same
at all wavelengths), it is said to be a grey body of emittance f;, and for
such cases E = f;aT 4 •
FREQUENCY (Hz)
15 1
,-.. J 10 , 10'\ 4 10 ' • 4
, ,(--;--"rj;'1Irr-I'I-1--.--;"--",,;,.",,-,,'-:;-,--.-,---:r--""T-'" " .-,(--"
E 10·
::L
, '0
,
6000b..L
~OOo':::'J
-----
I 1 ,
I! I
~ II
4000
~
(/) V/ ,/
,/ -~ooo
• 10
N
I
E '0• / I
/
V
/
,
20~'- ~ ~
--..... , ._-_.- i, l!
~
•
1/ /
1
I1
/
11
~~~ K- ~ ~~ -, ...
I
UJ
U '0
JI / / 1/ I
!/
/'
IOO<L::".
800~ ---
,
- ~
~
Z
<C • / / I
1
1// ' 60~
~
~"~
...... " ~
-
Cl
<C
Cl:
I()
10
I , 1/
fll
I
1
I
~--
I / ,
400
27
~ ~~
- ~ ............. + I
2~---: -......::.
I
I /
..J
<C .0
I 1 1/ I 1/
~ ,,~~ '.
"
0:
I- ,I IO- ,
I
I 1/ 1 / 1
/ / ,
"" ' '" , ,
'
U ,1
w-
0.1()
,I 1 I
I
I I 1/ I
/ I
I
"'"
C/)
, I I
I
1/
/
I
I
1
""" "-
<
.6
{
I
I
I
I
/ I
I
I "
...J 16"
01
o. 0.4 0.1 0.' 1.0 6 • 10 '0
WAVELENGTH (lJm)
I" , , I " I I I I I I! I I I
~4
•• ~J
•
WAVENUMBER (cm-')
33
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY.
very little emission of thermal infrared by the gases of the atmosphere and,
by Kirchhoff's law, in which there is little absorption. In the case of clear
skies more of the up welling IR flux from the ground in this spectral region
can be radiated directly out to space.
On the one hand, clouds scatter much solar radiation directly back to
space and thereby reduce the total solar input to the ground. There can be
rare and very specific exceptions to this - as, for instance, when two adja-
cent clouds funnel solar radiation downwards. In terms of reasonable ave-
rages, however, it seems that the 55 per cent cloud cover of the world re-
duces the solar energy absorbed by the planet by about 20 per cent. On the
other hand, the scattering of the direct solar beam by clouds is normally so
great as to ensure that for cloudy skies the total global input of radiation
to the ground (E.. ) is entirely made up of diffuse radiation Ed •
34
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
major cloud types. The mathematics of the calculation of the albedo of "bro-
ken" clouds is enormously difficult and only now being attempted. No simple
and practically useful rules have yet been devised.
. 10000 r----,----y-~--r----,r----o,...,._____,..-T'"""1
Cloud
N
I
E
Ci1
1000
L-
W
13
3
20 40 60 80
Solar Zenith Angle (degrees)
Table 3.1 The average short-wave albedo p and long-wave ~mittance £ as-
sumed to be associated with the major cloud types. These ·values show
relatively good agreement between observations and th.eoretical models.
They are appropriate for average sun height· and cloud thickness.
Cloud Type p £
stratus 60 1.0
Cumulonimbus 70 1.0
Nimbostratus 70 1.0
35
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
Clouds also absorb solar radiation, since liquid water and ice have ab-
sorption bands in the near infrared. However, these bands to a large extent
overlap those of water vapour, so the change in absorption by the whole
atmosphere is not very great. The main effect is to concentrate the absorp-
tion in the upper levels of the cloud (primarily because scatter in the
cloud increases the effective optical paths of water vapour).
This very simple concept can then be expressed in the follow.ing formula
for downcoming IR flux density under cloudy skies:
(3.32)
3.8..1 Albedo
The main difficulty is that, except over the ocean and over certain
highly u'niform land surfaces, such as extensive snow cover on flat country,
the variation of surface character leads to highly variable albedo. Even if
specific spot measurements are a:vailable, it is rare that these measurements
will be generally useful for any purpose other than that for which they were
taken.
To a certain extent this is true even for the ocean. Oceanic surface
albedo depends heavily on wind speed among other things, and is presumably
36
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
80 .-.-.-.-.- .........
--- --
/'-
,,'-.
".
/'
~ --..... '.
"
70
60 "-
p%
50
40
-·-1
--2
---3
l'"
.. -.-
'""".
-"-4
:: ~ ------~-
--5
30 ----6
---7
.-----
;: / ~~r-:" - - - --
··········8 J 'I
---9
20 - - - --10
10
37
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
11.12
PD - 2 J p(8) cos8 sin 8 d8 . (3.33)
o
,,
,,
--- --- --~oULTRAVIOLET
,
,, (FITCH)
,,
........ ',,/N<0.7} CAR ROLL AND FITCH"
......... ,N::::O.7
..... .X
90 ~ ",
\. .......... "
. \. .... ........ '"
\.. ..... ,
o \.. ANNUA .........." ...
o
UI
" " , AVERAGE ,..-~,-=.-:.-:.-:.-=.
ID
...J
et
" ..... ..... .....
80 ............
.......
.......
.......
.....
......
............ KORFF ET AL
10 20 30 40
SUN ELEVATION (0)
Figure 3.16 The albedo of snow as a function of sun elevation. The curve for
Kuhn et a1. (1977) represents' the average of four Antarctic stations,
that of Korff et a1. (1974) is for a snow field in the Rocky Mountains,
and the data of Carroll and Pitch (1981) and of Fitch (1976) are for the
South Pole.
38
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
Although the absolute values of snow albedo show some variability, the
change with sun elevation is quite similar at the different locations. From
an extensive ,series of measurements at the South Pole, Carroll and Fitch
(1981) developed the dependencies of albedo on cloudiness N and sun elevation
r given by:
peN < 0.70) = 0.98 - 0.007~j (3.34)
for cloud cover N < 0.70, and
Contrary to the situation for most other types of surfaces, the albedo
of snow is higher in the ultraviolet than at visible and near-infrared wave-
lengths. The data of Fitch (976) in Fig. 3.16 show a difference of 10 to 11
per cent between the albedo for the region 0.30 - 0.385 tLm and that for total
solar radiation. Spectral measurements of Kuhn and Siogas (978), summarized
in Table 3.2, show that surfaces composed of frozen water in some form have a
general decrease of albedo from visible to the near infrared. One exception
to this trend is a secondary maximum at about 1.4 tLm, which the authors
ascribe to inclusions of bubbles and snow inclusions in sea ice and to den-
dritic forms of crystals in fresh snow.
Albedo (%)
0.4 88. 47 22
0.5 88 47 19
0.6 88 40 17
0.7 87 30 10
0.8 87 20 6
0.9 86 16 5
1,0 77 11 4
1,1 74 11 5
1,2 61 9 4
1,3 55 8 5
1,4 66 15 -
1.5 25 4 -
1, 55 12 4 -
39
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
ultraviolet, blue, and red wavelengths only weakly, making for very lO~o1 re-
flectance in these regions. These features are shown in the corresponding
reflectance spectrum given in Fig. 3.15. A slight increased reflectance in
the green, combined with the spectral distribution of sunlight and the spec-
tral sensitivity of the eye, produce the green colour of most growing plants.
The near-infrared wavelengths are absorbed very little by most plant pig-
ments, a fact which permits this type of radiation to enter the main cellular
structure of plant tissues, where large differences of refractive index bet-
ween cells and cell walls are effective in producing a high reflectance at
infrared wavelengths. High values persist throughout the near infrared ex-
cept in the main absorption bands of liquid water. These bands are' strong
enough to cause a significant reduction in total reflectance of the plants,
as well as within the ba nds them selves. In addition to those shown in
Fig. 3.15, mean values of the reflectance of various agricultural crops at
visible and near-infrared wavelengths, as measured by Kondratyev et al.
(1964), are given in T'able 3.4.
Plant type .55 .:60 .65 .67 .70 .71 .75 .80 .85 .90
, ;
40
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
The albeqo of soils, sands, rocks, and other mineral-type surfaces has a
wide variation from as low as 2 to 3 per cent for wet black soil to as much
as 70 to 80 per cent for gypsum or quartz sand, certain clays, limestone, and
other light-coloured materials. More usual values for bare soils in farming
areas are 15 to 30 per cent when dry and 8 to 15 per cent when wet. Deserts
in general are higly reflective, as shown by the measurements of Otterman and
Fr a se r (1 976) l i s t e din Tab I e 3 . 5.
Table 3.6 Spectral reflectance (per cent) of sands and soils (Data from
. Condit, 1970).
Wavelength (#m)
Type state 0.31 0.40 0.50 0.60 0.70 0.80 0.90 1.00
41
CHAP.TER 3 RADIATION QUANTITIES IN METEOROLOGY
Very high reflectances are exhibited -by the small particles of many
clays and by the naturally light-coloured sands of quartz and gypsum, whereas
the dark-coloured soils prevalent in the midwestern part of the United
States, southern parts of the Soviet Union, and other rich agricultural areas
generally show low reflectances. Typical spectral reflectance values for
some of the different types of soil in the region 0.31 to 1.0 J.lm as measured
by Condi t (1970), are listed in Table 3.6.
42
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
The polarization of light reflected from natural surfaces has been stu-
died much less extensively than has that scattered by the atmosphere, and the
theory of reflection from rough surfaces is not as well developed as is the
theory of scattering. Polarization data are available for a few surfaces,
however, from measurements taken in laboratory conditions with artificial
light. The degree of polarization varies with type of surface, angle of inci-
dence, angle of reflection, azimuth angle from the position of the source,
and wavelength of radiation. In addition, the degree of polarization strong-
ly depends on the particle characteristics and roughness of the surface, and
a given surface polarizes reflected light much more strongly when it is wet
than when dry. A highly reflecting surface generally polarizes light less
strongly than does dark-coloured material.
3.8.4 Emittance
43
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
particle size may decrease spectral contrast and even introduce some absorp-
tion features of its own. The infrared reflectance of some typical minerals,
as determined by Hovis (1966), is shown in Fig. 3.18. The reflectance of
Oregon beach sand is particularly variable, but there is better regularity in
the features for the purer minerals. For instance, the curves for serpentine
and basalt show low· and relatively constant values at A < 8 jlm, whereas
those for other minerals have considerable variations in this region.
The one feature prevalent for all samples is the reststrahlen enhance-
ment in the 9-11-jlm region, and most samples have a similar reststrahlen
feature at about 20 jlm. There is evidence that the reststrahlen peak shifts
toward shorter wavelengths with increasing acidity of the mineral.
r.
/ i
I
,, - - SERPENTINE
I
60 I ,
I - - - BASALT
- - - DUNITE
\,
I _ ••- BEACH SAND
'\ (, ..- ...., GRANITE
I , ,1\ ,
,"\I
. ;. i I \ f: I '
<-!V\ \} \
I \
'\' ,
j : H\ ~
,,
.1\ \ r, •
i
I
I \.
{\I \ IV:1
~i I{\~i(',
\i I i \ , :'
\.; i ,/f'\ I I i i\
• •/ ., I I "
'.~: ~
\::
- "~
I 11
44
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
REFERENCES
Carroll, J.J. and B.W. Fitch (1981) Effects of solar elevation and cloudiness
on snow albedo at the South Pole. J. Geophys. Res. 86: 5271.
Coulson, K.L. and D.W. Reynolds (1971) The spectral reflectance of natural
surfaces. J. Appl. Het. 10: 1285.
Fitch, B.W. (1976) Albedo of the snow-covered Antarctic Plateau, M.S. Thesis,
Dept. of Land, Air, and Water Resources, Univ. of California, Davis, 59
pages.
Frohlich, C. and R.W. Brusa (1982) Solar Radiation and its Variations in
Time. Solar Physics 74: 209.
Gates, D.M. and W. Tantraporn (1952) The reflectivity of deciduous trees and
herbaceous plants in the infrared to 25 microns. Sc.ience 115: 613.
Hovis, W.A. and W.R. Gallahan (1966) Infrared reflectance spectra of igneous
rocks, tuffs, and red sandstone from 0.5 to 22.um. J. opt. Soc. Am. 56:
639.
Kasten, F. (1966) A new table and approximation formula for the relative air
mass. Archiv. Het. Geoph, Biokl. 206: B14.
Kondratyev, K.Ya., Z.F. Mironova and A.N. otto (1964) Spectral albedo of
natural surfaces, PAG£OPH 59: 207.
Korff, H.C., J.J. Gailiun and T.H. Vonder Haar (1974) Radiation Measurements
over a Snow Field at an Elevated Site, Pap.221, Dept. of Atm. ScL Colo.
State Univ., Fort Collins, (NTIS N74-3187813GD.
45
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
Kuhn, M., L.S. Ku ndla and L. A. Stroschei n (1977) The radia tion budget at
Plateau Station, Antarctica, 1966-1967, Antarctic Research Series, Am.
Ceop.!Jys. t/nJon 25: 41.
Lacis, A.A. and J.E. Hanson (1974) A parameterization for the absorption of
solar radiation in the earth's atmosphere. J. Atmos. Sc.l~ 31: 118.
LinkeF. and K. Boda (1922) Vorschlage zur Berechnung des TrUbungsgrades der
Atmosphare aus den Messungen der Intensitat der Sonnenstrahlung. Heteor.
ZeJtschr. 39: 161.
Li st, R.J. (1968) Smi thsonia n Met'eorological Tables, Smi thsonia n Institution,
Washington, D.C., 6th Revised Edition, p. 442.
Otterman, J. and R.S. Fraser (1976) Earth-Atmosphere system and surface re-
flectivities in arid regions from LANDSAT MSS data. Remote Sens. of Env.
5: 247.
Payne, R.E. (1972) Albedo of the sea surface. J. A tm. ScL 29: 959.
Willis, J. (1971) Some high values of the albedo of the sea. J. Appl. Het.
10: 1296.
Willson, R.C. (1984) Measurements of Solar Total Irradiance and its Variabi-
lity. Space Sci. Rev. 38: 203.
46
CHAPTER 3 RADIATION QUANTITIES IN METEOROLOGY
47
4. RADIATION MEASUREMENTS
a = =---- (4•1)
49
CHAPTER 4 RAD I AT 1ON MEASUREMENTS
40
1/ I
35
I
30
-
-
e1< 25
Q)
o
c
ca 20
I1 r r h ./
IV
i
~
o
If\ 1/
-
Q)
Q)
15
~
CC
10
I h .-. 1/
\~
'-' ! \ 11
I
I ,\ 1'-' _J p\\ / ~~
1
5
~/ V -,
11 '1/
/--
L\ff '-.-. - ~
_. -'-'
_2
HV ~
, I' --~
\ t--
"
"-
-....... ~
.b. ::.
/
\1
\
.- ...... r;:'.c':=-.
\
~~
---
I
t._-=- " ::::::.- --~":. -~~:::
::.~....
::'::.::: --- ~
l!
't-, t'- -- ',: ~. f:,-.: ..... .-._. -~
..--'-
o
OA 1.0 3.0 4,0 6.0 8.0 10.0 12.0 14.0 16.0
Wavelength (IJm)
Figure 4.1 Optical properties of black coatings. For identification of cur-
ves see text.
50
CHAPTER 4 RADIATION MEASUREMENTS
51
CHAPTER 4 RADIATION MEASUREMENTS
Before 1956 two "scales"· were used as references for meteorological ra-
diation measurements: the original Angstrom scale and the Smithsonian
scale 1913 revised. The Angstrom scale is based on the calculated perfor-
mance of the Angstrom compensation pyrheliometer A 70. The reference of the
Smithsonian scale is the Waterflow instrument with a self-calibrating cavity-
type. receiver and the Silverdisc instruments as transfer· standards.
The 1956 International Radiation Conference at Davos established the Interna-
tional Pyrheliometric Scale OPS 1956) by correcting the Angstrom and the
Smithsonian scales as follows:
These corrections were based on the assumption that the difference between
the Angstrom scale and the Smithsonian scale was 3.5 per cent as determined
by Angstrom (919), in spite of the fact that Kimball (1910), found this
correction to be 5.4 per cent, 5.3 per cent and 5.1 per cent respectively,
(Morikofer, 1937 and Hinzpeter, 1956).
The IPS 1956 was defined and used during the First International Pyrhe-
liometric Comparisons OPC 0, organized by the World Meteorological Organi-
zation at Davos 1959. The reference instrument was A 158, the Stockholm
standard, (traceable to A 70) corrected to IPS 1956 by changing its original
calibration factor by 1.5 per cent. According to the results of the compari-
sons, the calibration factors of all participating instruments were adjusted
in order to give the same readings as the reference. At IPC II in 1964 the
same procedure was applied. During the Regional Comparison 1969 of the Re-
gional Association VI at Carpentras, a substantial change in the readings of
the Stockholm instrument was noticed. This change was confirmed during
IPC I I I in 1970. A 158 measured 1.2 per cent higher than during IPC I and
IPC 11. Later, some of the discrepancies were explained by an irregularity
in the auxiliary equipment used during IPC I and 11. Therefore, the
"IPS 1956" used after 1959 did not represent the scale as defined by the
Radia tion Com mission in 1956.
Table 4.1 Differences between different scales in use for meteorological ra-
diation measurements
These results demonstrate that the assumed differences for: the defini-'
tion of the IPS were wrong, and that· it is difficult to establish such a
reference with the traditional pyrheliometers. Moreover, it suggests that
52
CHAPTER 4 RADIATION MEASUREMENTS
Figure 4.2 Sum mary of the relations between different scales in use for
meteorological radiation measurement s (Frohlich 1973a).
WR R
= 1. 022. (4.3)
"IPS 1956"
With this value and those given in Table 4.1, results of previous stu-
dies can now easily and clearly be reduced to the WRR.
Mandatory use of the WRR has been put into the regulations of the World
Meteorological Organization in 1981. Its absolute uncertainty is estimated
to be ±0.3 per cent, which means that the true value lies within this range.
53
CHAPTER 4 RADIATION MEASUREMENTS
ACR 310
ACR 311
ACR 701
CROM
ECR 10
ECPR NS21-11A
~~-ljEPAC 11402
EPAC 12843
EPAC 13617
PMO 2
ifii'.imR*t--_..-------. PMO 3
PMO 5
PVS-5
TMI 67502
PACRAD III
4.2.1 Pyrheliometers
54
CHAPTER 4 RADIATION M~ASUREMENTS
Front Aperture
Full viewAngle
Slope Angle
Receiver
S = kP, (4.4)
where S is the solar irradiance, i the measured current and k is the in-
strument constant. In principle k can be calculated according to
k = cr/ab, (4.5)
55
CHAPTER 4· RADIATION MEASUREMENTS
56
CHAPTER 4 RADIATION MEASUREMENTS
ciated with each type of instrument, due to deviations from ideal behaviour,
must be applied in order to convert the observed electrical power to a true
irradiance quantity. The most important items cover the accurate measurement
of the area of the precision aperture, errors in the electrical power mea~
surement, the emission from the cavity through the entrance aperture, heat
losses through the power or sensor leads, and a number of items which are
generally grouped under the term "non-equivalence" and which relate to the
differences in the temperature distribution and, thus, the thermal paths of
the radiatively and electrically produced heat in the cavity. The absorp-
tance of the cavity receiver over the spectral range of the incident solar
radiation must also be accurately assessed.
Identification Sensor and mode Cavi ty type Cav i ty-heat fl Oil Aperature Character i zat ionl )
The drawings ill ustrate the cavi ty-heat f 10" JIleter arrangements:
- cavity. =·heat-flowaeter. '7777777 heat sink
Each individual instrument has its own set of correction factors just as
the older pyrheliometers had their own calibration constants. The absolute
57
CHAPTER 4 RADIATION MEASUREMENTS
accuracy for all described instruments is between ±O.3 and ±O.2 per cent,
depending on the realiabilityof the characterization. Some instruments, such
as the PACRAD and ACR, have been built by different manufacturers. The TMI,
H-F, and PMO radiometers are also in commercial production. The ACR, how-
ever, is no longer commercially available (19·84).
The terms. "active" and "passive" have been employed with regard to the
operation of absolute radiometers. Generally speaking, "passive" devices are
operated in a sequential manner by alternately exposing the receiver to the
sun and then electrically calibrating the heat-flow transducer with no elec-
trical power being applied during t.he sunlit phase. "Active" operation gene-
rally refers to the maintenance of a constant temperature difference along
the heat.-flow resistor during both exposure and reference phases of a se-
quence controlled by a servo or other control system which automaticalJy
adjust.s t.he power to the heater. None of the inst.ruments, except CROM, is
truly bifurcated as was the case with the Angstrcim and the double-tube water--
flow instruments. But Toast are employed with a compensating receiver to mini-
mize the effects of changes in ambient conditions on the instrument perfor--
mance. A number of automatic cont.rol mechanisms and computer interfaces have
been developed for operation of these devices.
Contrary to the i nstrum ents described under 4.2.1.1 and 4.2.1.2, which
are mainly used as reference standards for calibrating the instruments de-'
scribed below, the inst.ruments discussed here are used operationally. The
basic sensing element in these instruments is a thermopile used as a thermal
flux meter mount.ed on a heat sink and painted black on the surface exposed to
the incident radiation. They have to be calibrated by comparison with a
reference standard pyrheliometer.
58
CHAPTER 4 RADIATION MEASUREMENTS
The Kipp & Zonen Pyrheliomet.er corn monly referred t.o as t.he Linke - Feus-
sner Act.inomet.er ("Panzer--Akt.inomet.er"), is shown in Fig. 4.6. It. has a
Moll-t.ype t.hermopile element. wit.h fort.y junct.ions of manganin-const.ant.an ar-
ranged in two groups of twent.y each ina. circle. This arrangement. offers
compensat.ion for inst.rument.-t.emperat.ure fluct.uat.ions. The unit. is fitt.ed with
a Eil t.er wheel for t.he selection of either no iil ter, a quartz glass, or one
of t.he classical Schott filters. There is no t.emperat.ure compensat.ion, but
the calibration constant can be corrected by the factor 1 +a·(t.-20), with t
being t.he instrument temperature and a = O.002K-1. A thermometer is mounted
on the instrument body for establishing the temperatur~at the time of mea-
surement. This pyrheliomet.er is usually supplied on a mount allowing an ad-
just.ment for solar elevation and azimut.h to be made by hand .• The performance
specifications are summarized in Table 4.3.
59
CHAPT'ER 4 R AD I AT ION MEASURE ME'NTS
The pyrheliometer used in the· Soviet Union and E'astern European coun-
tril;!s is the actinometer A80 developed at the Main Geophysical Observatory in
Leningrad CYane.shevky, 1957), The receiver is a thermopile, soldered in the
form of a st.ar, of manganine and constantan bands CFig.4.7a:J. The cold
junctions of the thermopile are fixed on a, copper ring with a thin layer of
electric insulator inbetween. This ensures g.ood thermal contact with the
a·ctinometer body. A thin silver-foil disk of 11 mm in diameter is glued to
the top of the hot junctions, which is painted with a mat black lacquer.
aJ
Figure 4.7 a) Construction of the copper-m.anganin thermopile of the A80
actinometer. 1: hot junctions, 2: cold junctions, 3: copper ring, 4:
silver foil disc. b) A80 actinometer. 1: actinometer, 2: coarse a.nd 3:
fine adjustment of hour angle 4: support, 5: cover, 6: base plate.
The general view of the actinometer is given in Fig. 4.7b and the per-
formance specifications are summarized in Table 4.3. The instrument is nor-
mally fixed on a clock-driven equatorial mount.
Sensivity 8 16 10 /lV/Wm-'
View limiting
geometry:
full view angle 5.7' 9.6 10 . deg
slope- angle 0.8·' 2.4 2;3 deg
60
CHAPTER 4 RADIATION MEASUREMENTS
Since pYl;'heliometers need view only a small field centered on the solar
disk, the mounting site requirements are less restrictive than for most types
of radiation instruments. The tracker carrying the instrument should be
mounted on a stable platform from which a view of the sun is not restricted
by buildings, towers, cables, or other obstruction$ throughout the sunrise -
sunset period. In the choice of a mounting site, the change of - solar azimuth
with season must, of course, be taken into account. But obstructions of the
horizon have no other significance for measurement of the direct solar beam.
Step 1: Set the mount on a stable and approximately level surface with the
hour-angle axis in approximately the north-south plane. Level the mount
by means of a spirit level.
Step 30: . Orient the pyrheliometer to view the zenith by adjust'ing the decli-
nation and the hour angle.
step 4: Precisely at solar noon, rotate the instrument around the decli-
nation axis to approximately view the sun. By small adjustments of the
instrument around the declination axis and the mount in azimuth but not
around the hour axis, orient the instrument precisely to the solar posi-
tion. The mount should now be in the proper orient a tion, but it is well
to check the incl~nation angle of the hour-angle axis again to be sure
it has not been disturbed.
61
GHAPTER 4 RADIATION MEASUREMENTS
4.2.2 Pyranometers
62
CHAPTER 4 RADIATION MEASUREMENTS
over an azimuth range of 360 degrees about the zenith vector. Pyranometers
are usually used for routine measurements and left permanently outside. For
this reason, they must be able to withstand the effects of all types of
meteorological condit.ions. Most. inst.rument.s are of t.he thermoelectric typE'
with the detector protected by glass hemispheres, which also define the spec-
tral range of the radiat.ion being measured. In the case of a black receiver
two domes -are required in order to shield the infrared radiation from the
outer dome, the temperat.ure of which may vary quite rapidly with changing
wind speed and air temperature. Pyranomet.ers have recently been developed
which employ sensors of various t.ypes of silicon diodes (see 4.2.2.8).
A short description of the represent a ti ves of the mo~.t. widely used in-
strument t.ypes: Eppley PSP, Kipp & Zonen eM5 and CM10, the Schenk stern
Pyranometer and the Russian instrument M-BO will be given below; This list is
by no means 'exhaustive, but it is representative. Furthermore, these intru-
ments have been tested unde~ a great variety of conditfons and were compared
many times with each other and with other types during int.ernational trials.
Therefore, a great deal of useful information is available to illustrate the
behaviour of these instrument types. This has helped considerably in their
selection for field use. The general specifications of these instruments are
sum marized in Table 4.4.
Schenk 15 -0.1%IK 6
M - 80M 11 -0.15%IK 8
--------_._-
63
CHAPTER 4 HAD IATI ON MEASUREMENTS
The Eppley PSP (Precision Spectral Pyranometer, Fig. 4.8) has a circular
multijunction t.hermopile of t.he t.ype wound wit.h copper-plated const.ant.an
wire. The receiver is coat.ed with Parson's optical black lacquer. It has
t.wo precision-ground and polished hemispheres of Schott glass WG295 (short-
wave cutoff, 50 percent transmission at 295 ,um),which may be replaced bya
coloured glass dome, hence the name "precision spectral" (see 4.2.4.1).
The Kipp & Zonen CM5 comprises a Moll-type thermopile of very simple and
rugged construction. The domes are made fromSchott glass K5" which has a
short-wave cutoff at about. 310 nm. The thermopile is coated with lamp-black.
The recentlydevelopedCM10 has a thermopile of hundred thermocouples printed
on a ceramic substrate, using Cl thick-film technique. The blacj{ coating is
the 3M velvet black. The material of the domes is the same as that for the
CM5. A CM-iD is shown in Fig. 4.9.
64
CHAPTEf~ 4 RADIATION MEASUREMENTS
The Russian instrument is also of the black-and-white type and was deve-
loped at the Main Geophysical Observatory in Leningrad (Fig. 4.11). The re-
ceiver is a flat square thermopile built up of pieces of manganine and con-
stantan strip soldered together. The black and the white paints have the same
factor 0.£ absorption in the long-wave spectral region, which makes the in-·
strument sensitive .only to short-wave radiation. The thermopile is protected
by a thin hemispherical glass dome.
KIPP &
KIPP &
KIPP &
KJPP&
KIPP &
ZONEN
ZONEN
ZONEN
ZONEN
ZONEN
CM10
CM10
CM10
CM10
CM10
790059
810119
810120
·810121
810122
5.58
4.58
4.54
4.66
4.24
5.65
4.59
4.52
4.62
4.22 ..
'\.
SCHENK STAR 1626 14.26 14.49
SCHENK STAR 2186 14.94 15.15
SCHENK STAR 2209 15.36 15.29
SCHENK STAR 2217 14.16 14.17
SCHENK STAR 2221 15.24 14.97
65
Outdoor and indoor tests have been performed in order to provide infor-
mation of the performance of pyranometers. A very exhaustive investigation
has been completed recently (TEA Task III (984)), during which the perfor-
mance of more than 30 instruments of eight different types was studied. The
outdoor tests covered a wide range of conditions: winter, summer, clear
skies, cloudy skies, and overcast skies. The temperature dependence, lineari-
ty, influence of tilt and directional response (cosine 'law) were meas'ured in
the laboratory for all instruments. The ultimate goal of the trials was to
determine hot\! -a.ccurate or how variablepyranometers ,can be under different
conditions and i f the accuracy can be substantially improved by introducing
corrections based on the laboratory results.
't. " ',JuI2,9 6 .Jul 30 Jon·13 Jan 14 '" Jan ~8; ",Ju129 6JuI 30 + Jan13 x Jan 14 " Jan18
. .:.'::,:'"=~:'I
.,0 + x
-10'"
15.00
6 q-C%"
' - .~.'..
0- . ~
1.2~OO
• ..
V·,
( ",,~
~.
.14.00
'" ".
118.00
-~ ' & ~' 'ii
18;00
']==,:"
-6:00
'£1:"":=4<:'"
10iOO
.~= 1'8.00
+10,'" Of- - - - - - - - - - - - - - - -
c
j .. ":.~. 15 20 25 30 '0 25
6
30
10 15 .20
.2
'1Q
:; Instrument Temperature {O'C) Instrument Temperature (OC)
III
o
"°:l 6
1!l
x~""~l ~--,----'
","10'1.1--1~-~x- _ _- - _ _-"-~_~
o 200 400 600 ::-800 1200
2
:Intensity {Wm c ) --hilensity(Wm '2)
'20 -40 aD 80
Figure 4.13 Comparison of the performance of an Eppley PSP (lef.t) and a Kipp
& Zonen eMS (right) as a function of time of day, temperature, intensity
and' incident angle.
From these results it follows that the Eppley PSP and Kipp & Zonen CMiO
perform quite well and give measurements within 2-3 per cent of each other
over a very large range of conditions. On the other hand, the GMS and the
Schenk stern show systematic differences of lIP to 6 to 9 per cent. This
deviation is partly due to the temperature compensation in the former two in-
struments, hut for the latter instruments the linearity and directional re-
66
CHAPTER 4 RADIATION MEASUREMENTS
----- ----
PSP
~o+---=="'~~:---===-------
w.
a.
o-l
Ul
<l -1
-2
-4
-J +---Ir----+--+---f---!----+--+--~I- -s-J-t-----i--+---f---!---+--r---l-
100 200 JOO ~OO ~OO 600 700 BOO -JO -20 -10 0 " ·-10 20 30 40
67
CHAPTER 4 RADIATION MEASUREMENTS
Clear, summer -0.5 +0.3 -0.2 +0.7 ~2.2 -0.3 +0.4 + 2.2
Cl ear, wint,er -2.4 +0.8 -0:.7 +1.8 -5.7 -3.5 +7.7 +10.4
C.loudy, summer -0.5 +0.1 +0.2 +0.5 -1.7 +0.5 +0.8 + 1. 9
Cloudy, winter -3.0 -1.2 +1.1 +2.0 -1. 8 +0.7 +4.4 + 8.2
Overcast, winter -1.3 +0.6 +1.2 +2 •. 0 +0.6 +3.9 +3.0 + 3.,8
Up to now npthing has been said about long-t.erm stability. and about how
often the different types of instruments need to be calibrated. As shot.ln in
Fig. 4.15, the long-term variability is not very important compared with the
short-term variations due to changing conditions during. the measurements. As
to the calibration stra.tegy it is importantt.hat they are always performed
under similar conditions, otherwise the results will be inconsistent. This
holds true mainly for the less accurate types. Ag.ain, it is ast.onishing how
well the CM-5 performs at the level of daily sums.
::]~._=_~_=_~_~__: : : = = :=-:~:.-q. . .: :.~ CM~6~':-= - -: : :f=>:~" " " ': ~.~ ~'~:·~:=I
: : : .: .. .:.. . ::K:I:P:P:": .•
68
CHAPTER 4 RADIATION MEASUREMENTS
(4.7)
N
AE' = ~ l1E n , (4.8)
n=1
69
CHAPTER 4 RADIATION MEASUREHENTS
Pyranometers should be calibrated at lea.st once per year, and more fre-
quently if there are reasons to sus.pect unusual changes of instrument perfor-
mance. The most precise calibration, and the one g.enerally recommended, is
provided by anyone of several national or regional radiation centres which
have special equipment and expertise in radiometer calibrations. This way of
calibration, of course, requires transport of the instrument to and from the
centre, and its use in the field is interrupted for a period of time. A
second, more appropriate method is to transport a secondary standard pyra-
nometer to the field for a direct comparison with the operational instrument.
In this method the secondary standard is installed alongside the operational
instrument and a statistical sample of their records is obtained under a
variety of radiation conditions, including clear and cloudy skies, over a
period of about two weeks. As for the case of calibration of pyrheliometers"
a relatively large number (at least 100) of simultaneous samples of thesig-
nals Ai and A2 aTe- taken, and the average R of the individual ratios A2 /Ai
computed. If Ci and, C2 are the calibration factors of the operational and
standard pyranometers, respectively, then we determine the calibration factor
for the operational instrument as Ct = C2 • R. This should be checked
against previous calibration factors for the instrument, and any discrepan-
cies not,ed. In case there has been a significant change since the previous
calibration, the process should be repeated within a. feww'eeks to detect any
progressive deterioration of the instrument. If such deterioration conti-
nues, the instrument should be replaced and a new calibration performed.
70
CHAPTER 4 RADIATION MEASUREMENTS
ings of the shaded pyranometer Vs and the pyrheliometer Vp and their calibra-
tion factors Cs and Cp the global radiation incident on a horizontal surface
can be calculated as a function of ro the solar elevation by:
(4.9)
From this value and the reading Vc of the pyranometer under test its calibra-
tion factor can be calculated:
Cc = E';'/V c . (4.10)
With this expression and some assumptions on the sky radiance, a correction
factor f can be determined by:
f = (1 - Dv I ED ~ ) - t (4.12)
71
CHAPTER 4 RADIATION MEASUREMENTS
...0 1.10
u
-IV
-......
0
u
11 1.05
0
()
1. 00 +--r-------r----~---__,...----_.__--+
-23.5°.20°
Declination
72
CHAPTER 4 RADIATION MEASUREMENTS
between disk and sensor will be constant, thereby making the solid angle
subtended by the disk independent of solar declination.
For practical reason, the best approach when choosing a pyranometer for
use in a tilted position would seem to consider the problem and calibrate the
instrument in the tilted position with a shading technique.
calibrated und&r average climatic conditions at th~ station where the instru-
ment will be used. But this is very difficult to quantify and, therefore,
these instruments are not· recom mended for accurate climatological measure-
ments.
When interest is mainly in the net radiation, Le. the difference be-
tween the incoming and the outgoing radiation fluxes, rather than in the two
quantiti&s themselves~ one heat-flow transducer can b& used by exposing both
sides, the "cold" and the "hot", to the radiation. In this case the instru-
ment is called a net pyrradiometer.
The instruments in use are of two different types: the unshielded and
the shielded types. It is obvious that the un shielded radiometer has an
inherently better optical performance, but as the sensors are bas&d on ther-
mal heat-flow meters, varying losses to the surrounding air introduces er-
rors. Courvoisier (1950) and others (e.g. Geir and Dunkle, 1951) tried to
overcome this problem by blowing an "air curtain" of constant speed over the
sensors. Although the influence of wind is drastically reduced, it is not
eliminated. The construction of such devices is rather involved and they are
normally not used for operational measurements, but as a reference for cali-
bration. The Yanishevsky net pyrradiometer (Fig. 4.17a), which is also un-
shielded, but not artificially ventilated, uses individually determined cor-
rection factors for the influence of wind. These typically range from about
1.03 at lms- i to 1.22 at 10ms- i • The sensing element is a thermopile of
silver and constantan which consists of ten sections connected in series. One
section consists of a rectangular copper bar (Fig. 4.17b) around which a thin
ribbon of constantan is wound over a thin insulating layer. One half of each
loop is electroplated with silver. Thus the junctions between the silvered
and the non-silvered parts, the hot and cold junctions, are situated on oppo-
si te sides of the bar.
Ten bars - side by side - form the complete thermopile. A copper foil is
glued onto each of the two surfaces with an insulating spacer and painted mat
black. The sensivity is about 12 J./V/Wm- 2 and the temperature difference
between the hot and cold junctions at 1 kW/m- 2 is only about lK. Under' such
conditions the infrared radiation of both sides of the receiver can be consi-
dered the same. It is obvious that this instrument can only be used occasio-
nally and not continuously. Th& wavelength dependence of the absorptance of
the black coating is one of the sources of errors of the unshielded instru-
ments. In case of most coatings, the absorptance decreases with increasing
wavelength and may eveD have some line features at wavelengths where the
74
CHAPTER 4 RADIATION MEASUREMENTS
a) b)
n
Figure 4.17 a) External view of the balance meter M-l0M. 1: receiving sur-
face, 2: frame, 3: knob, 4: case. b) One section of the thermopile of
t.he balance meter M-l0M. 1: constantan strip, 2: copper bar, 3: electric
insulator, 4: silvered sections
75
CHJ'lPTER 4 RADIATION MEASUREMENTS
organic structure. As in the case· of"" the black coating, this introduces
further uncertainties which are difficult .to quantify. In order to minimize
this effect very thin polye.thylene sheets have to be used ~h:i.ch are no longer
self-supporting. Therefore the domes have to be ~ept in shape by pressurizing
the interior. The sa me technique is used for the net pyrradiometer, developed
by Funk (959) and manufactured by Middleton, Australia, and Swissteco,
Switzerland. 'An exanrPle of such an instrument is shown in Fig. 4.18 and ,the
specifications of the-'two types are summarized in Table 4.5. ,.
Table 4.5 Specification of Funk type pyrradiometers and the Eppley pyrgeo-
'-meter'
Time constant 15 s 6 s 2 S
(1 le)
Temp. coefficient
of sensitivit.y )+0.012%/K +1% for -20 to 40· C
76
CHAPTER 4 RADIATION MEASUREMENTS
to have a nearly grey transmission from 4 to 50 ~m. These domes are used in
the Eppley pyrgeometer shown in Fig. 4.19 The sensor and design of the in-
strument is the same as for the Precision-Spectral Pyranometer and the gene-
ral specifications of the pyrgeometer are listed in Table 4.5. In contrast to
the net pyrradiometers a single instrument needs to take its own infrared
emission into account. The corresponding signal can be calculated from the
instrument's temperature according to the Stefan - Boltzmann law or is, as in
the Eppley pyrgeometer, automatically eliminated from the signal by applying
a battery voltage which is precisely controlled by a thermostat sensing the
temperature bf the instrument.
A final and very si mpIe met hod, s u g ges ted by Id so (1971), utilizes a
large temperature-controlled flat black plate. The sensing element is posi-
tioned directly over or under the plate. If surrounding laboratory condi-
tions are constant, then the net irradiance at the sensor is E =Ee - ~€aT~.
Here Ec is the irradiance (considered constant) from the surroundings, € is
the plate emittance, a is the Stefan-Boltzmann constant, T is plate tempera-
ture, and ~ is the angle factor given by Idso (1971). The output signal V is
V = C~ - K'€a~, where C and K are the calibration factors for the radia-
tion from the room and from the plate, respectively. In performing the cali-
bration, values of V ~re obtained as a function of T. The value of K is
given by the slope of the curve of V plotted versus the quantity ~€aT~. Idso
estimated the accuracy obtainable by the method to be ±5 to ±3 per cent.
-The first quantum detectors in use were the cadmium cells for the mea-
surement of ultraviolet radiation. As the operation was very cumbersome, they
77
CHAPTER 4 RADIATION MEASURRMENTS
were never used operationally. Selenium cells were later introduced for de-
termination of the spectral extinction of the direct solar beam. The modern
sunphotometers were developed from these later instruments.
Table 4.6 Transmittance values at various wavelengths for Schott glas's fil-'
ters (2 mm thickness, +25'C) as given by the Eppley Laboratory, Inc.
(CC: 50 per cent point of cutoff; FF: filter factor; wavelength in #m)
wavelength band covered by the filter in Wm- 2 • If the cut-offs are less
steep and the transmittance is lower than one, or even varies in the
transmitting band, the interpretation of the reading is more difficult. Be-
cause the result will depend on the spectral distribution of the incoming
radiation, the so-called filter factor, FF, has been introduced. This compen-
sates for the attenuation of radiation in the transmission band. The FF of
a specific filter can be calculated by convoluting the spectral
transmission of the filter with the spectral irradiance of the incoming
solar radiation and comparing the result with that found with an ideal fil-
ter. Such a FF, however, will obviously depend on the state of the atmo-
sphere as does the spectral distribution. But if a typical condition is cho-
sen for the calculation, the FF may be representative for some range of
conditions as long as the variation is small enough for the purpose of
the measurements. The first FF was calculated as described above, by using
measured filter transmissions and a typical solar spectrum as mea-
sured, for instance, by the Smithsonian Institution from mountain tops in
their programme to determine the solar constant. This is a very cumbersome
way to calibrate individual filters. The FP's were therefore normally deter-
mined by comparison with standard filters, which had their FF's derived
from calculations. This method was applied at the Davos Observatory for
many filters as a service. One severe drawback to this method is that
the intrinsic transmittance curve of the filters to be calibrated may be
different from those of the standard. This may happen if the filters are from
different batches, e.g. the standard filters, manufactured before
World War Il, have a substantially different transmittance towards the in-
frared cut-off than those delivered later. But they have still been used for
calibrations. Thus most of the more recent FP's should be applied with some
caution and a calculated redetermination should be made from the individua-I
transmittance curve and a representative- terrestrial solar spectrum. Another
problem is the shift of the short-wave cut-off with the temperature
of the filter. The FF changes by 0.02 %/K for OG530, by 0.03 for RG630 and
by 0.04 for .RG695 (e.g. Angstrom and Drummond, 1961). As it is very diffi-
cuI t to measure the temperature of a filter accurately, corrections are
difficult to apply.
4.3.2 Sunphotometers
79
CHAPTER 4' RADIATION MBASUREMEN'TS
Measurements at 610 and 675 nm can be used together with the 500-nm
channel for the in situ determination of ozone from absorption in the Ghapuis
hand, To improve the information in ozone-free spectral regions, channels at
412 nm and 450 nm may be added. These channels may be used as reference for
N02 determinations and together with the 778-nm channel as reference for
ozone determination. Channels at 719, 817 and 946 nm have been used to
determine the columnar content of water vapour by comparing their optical
depth with those at the 675-, 778-, 862- and 1024-nm channels (e.g. Koepke
and Quenzel, 1978).
A fully equiped set-up would therefore have twelve channels, but if only
aerosol optical depth is sought, the three mandatory WMO-channels are
sufficient. Table 4.7 lists all the wavelengths which can be used in sunpho~
tometry available with choices for different applications and possible
combinations for three-, six- and nine-channel instruments. Fig. 4.20 shows
the extraterrestrial spectrum together with the distribution of the channels
and a curve indicating the sig.nal received from a typical silicon cell.
Sll
CHAPTER 4 RADIATION MEASUREMENTS
111
N
r--------------------------,~
rl
E
c
IS)
N
,..,
a:
"<
(\j
E
--
111
:J
L..J
111
"
~
L..J
IS)
+'
C
ID
L
ID L
tJ IS)
:J
C tJ
I'd 0
+'
"U 0
Itl ..!:
L
L
a...
H
ISl
ISl
III -
ISl
ID
U1
J'.
ID - -
lTl
J'.
aI
J'.
J'.
J'.
aI
N
ID
aI
ID
'Of"
lTl
Aerosols
minimum
extended * * *
Ozone * * * * * *
Water vapour
* * * * *
weak * * *
* *
strong
Nitrous Oxide
* * *
* * *
3 Channels
6 Channels * * *
9 Channels
* * * * * *
12 Channels
* * * * * * * * *
* * * * * * * * * * * *
81
CHAPTER 4 RkDIATIOW MEASUREMENTS
/ /29405 /
/ //
~ Electronics
III
,''''-
Shutter Motor Hea~Wln~
(1 of 3) Temperature Sensor
82
CHAPTER 4 RADIATION MEASUREMENTS
The Langley method is based on the fact that the logarithm of the irra-
diance is a linear function of the air mass (see also Eq. 3.15). Thus V. can
be obtained by extrapolating the values at different air masses to air mass
zero. This method, although very simple in principle, relies on the stabili-
ty of the atmosphere during at least half the day, which is very difficult to
check. Thus, this method can only be applied from high-altitude stations and
under conditions with very clear sky. Moreover, a great deal of experience is
needed to judge the quality of the results obtained even under more or less
ideal condition s.
In recent years much progress has been made in sunphotometry and its
calibration from high-altitude balloons and rockets (Brusa et aI., 1983) so
that reliable reference instruments which are directly calibrated with extra-
terrestrial sunlight may be available in the near future. In any case, it is
recom mended that comparisons be made with a stable reference instrument to
standardize field instruments so that the V. can be traced back, whenever a
reference of higher accuracy becomes available.
83
CHAPTER 4 RADIATION HEASUREMENTS
Two different techniques have been developed for measuring the albedo of
natural surfaces. The prevailing method involves the use of an inverted
pyranometer, whereas the use of an integrating sphere has been mainly for
special studies. In the following only the pyranometer method. is described.
85
CHAPTER. 4 RAD IA T I ON MEASUREMENTS
86
CHAPTER 4 RADIATION MEASUREMENTS
Blevin and Geist (974) have developed a very simple reflectometer with
a large-area pyroelectric detector made of Polyvenylfluoride and mounted
inside a cone. For the measurements, a thin cone is placed on the surface,
the reflectance of which has to be determined. Chopped infrared radiation
from a spectrometer illuminates the sample through a small hole on the top of
the cone. The reflected radiation is then collected by the large-area pyro-
electric detector. Thus a direct determination of the reflectance is al-
lowed. As pyroelectric detectors are only sensitive to chopped radiation,
the thermal radiation from the sample will not be measured.
where X is the angle of the plane of polarization of the incident light. The
degree of linear polarization P is given by the relation:
(4.16)
L = La + Lr • (4.17)
87
CHAPTER 4 RADIATION MEASUREMENTS
The use of such a system is outside the scope of the present discussion.
One such system has been described byCoulson (1975), and an extensive set of
measurements using the system was made by Sekera et al. (1955).
Abbot, C.G. (1922) The silver disk pyrheliometer. Sm.ith. H.isc. -Coll. 56: 19.
Aldrich, L.B. and C.G. Abbot (1948) Slllithsonian pyrheliometry and the
standard scale of solar radiation. Smit1J.H.isc. Coll. 110~ 5.
Angstrom, A. and A.J. Drummond (1961) Bas.ic Concepts Concerning Cutoff Glass
Filters used in Radiation Measurements. J. Het. 18: 360.
Beckman, W.A., S.A. Klein and J.A. Duffie (1977) Solar Heating Design by the
f-Chart Method.Wiley and Sons, New York.
Berger, D.S. (1976) The Sunburning Ultraviolet Meter: Design and Performance.
Photochem.istryand Photo.b.iology 24: 587.
Blevin, W.R. and W.!. Brown (1966) Black coatings for absolute radiometers.
Hetrologia 2: 140.
Brusa, R.W., C. Frohlich and Ch. Wehrli (1983) Solar Radiometry from High-
Altitude Balloons. 6th ESA Symposium on European rocket & balloon pro-
88
CHAPTER 4 RADIATION MEASUREMENTS
Coblentz, W.W. (1913) Instruments and methods used in radiometry. If. Bul.
Ha t. Bur. std. 9: 7.
Coulson, K.L. (1975) Solar and Terrestrial Radiation: Methods and Measure-
ments, Academic Press, New York.
Dave, J.V. (1978) Performance of a tilted solar cell under various atmosphe-
ric conditions. Solar £nerqy 21: 263.
Drummond, A.J. (1956) On the Measurement of Sky Radiation. Arch. Het. Ceo-
phys. B.io.kl. Ser.ie B, 7: 413.
Flowers, E.C. (1978) Solar Radiation Facility - Report for 1978, National
Oceanic and Atmospheric Administration, Boulder, Colorado (unpublished).
Flowers, E.C. (1984) Longterm behaviours of· pyranometers. Proc. lEA Meeting,
Norrkopping, Sweden, January 1984.
Frohlich, C.· (1973a) The relation between the IPS now in use and ·Smithsonian
Scale 1913, Angstrom Scale and Absolute Scale. Proc. Sym. Solar Radia-
tion, Smithsonia Institution, Rockville MD, Nov. 13-15.
Geir, J.T. and R.V. Dunkle (1951) Total hemispherical radiometer.s. Trlins.
Am. Inst. £lee. £nq. 70: 339.
Geir, J.T., R.V. Dunkle and J.T. Bevans (1954) Measurement ,of absolute spec-
tral reflectivity from 1.0 to 15 microns. J. apt. Soc. Am. 44: 558.
Goldberg, B. and IV.H. Klein (974) Radiometer to Monitor Low Levels of Ultra-
violet Irradiance. Appl. opt. 13: 493.
89
CHNP'fER 4 RADIAT ION MEASUREMENTS
Idso, S.B. (1971) A simple technique for the calibration of long-wave radia-
tion probes. Aq.r..ic~ . Heteor. 8: 235.
lEA Task HI, (1985) Results of an Outdoor and Indoor Pyranometer Comparison.
In press.
Karoli, A.R., J.R. Hickey and R.. G. Frieden (983). Self calibrating Cavity
Radiometers at. the Eppley Laboratory: Capabilities and Applications.
Froc. SPIE, East Arkington, April 19'83.
KendaU, J.M., Sr.' and C.M. Berdahl (1970)· Two black body radiome.ters of high
accuracy. Appl. opt. 9: 1082.
Kimbal1, H.H. (1910) Solar radiation, atmospheric absorption, and sky polari-
zation at Washington, D.C. B.uIl. /fount W'eatl1er Ohs. 3: 69.
Kmito, A.A., V.A. Parfinsky, M.M. Seredenko and V.A. Klevantsova (1976) Spec-
tral coefficient of blackened surface reflection. "HGO Proceed.inqs" 370:
39.
Koepke, P. and H. Quenzel (1978) Water vapor: . spectral transmission a.t wave-
lengfits between 0.7 J.lm and 1 J.lm. Appl. opt. 17: 2114.
Le Baron, B.A., W.A. Feterson and 1. Dirmhirn (1980) Corrections for diffuse
irradiance measured with shadowbands. Solar Enerqy 25: 1..
Uu, B.Y.H. and R.C. Jordan (1963) Daily insolation on surfaces tilted toward
the equator. ASJlRAE Journal 3: 53.
FeterS'on, J .T., E.C. Flowers and J.H. Rudisill (973) Dew and fros·t deposi-
tion on pyranometers.. J. App.l. Heteor. 12: 1231.
Reid, M.S., C.M. Berdahl and J.M. Kendall (1978) Calibration standards and
field instruments for precision measurements of insolation. Solar li'nerqy'
20: 357.
Sekera, Z., K.L. Coulson, D. Deirm.endjian, R.S. Eraser and C. Seaman (1955)
Investigation of polarization of skylight. Final Rep., Contr. AF19(122)-
239, Univ. o·f Calif., Los Angeles, California.
CHAPTER 4 RADIATION MEASUREMENTS
Temps, R.C. and K.L. Coulson (1977) Solar radiation incident upon slopes of
different orienta tions. ~(.701ar E'nerg-y 19: 179.
Willson, R.C. (1973) Active cavity radiometer. AppJ. opt. 12: 810.
Willson, R.e. (1980) Active cavity radiometer type V. AppJ. opt. 19: 3256.
Zalewski, E.F. and J. Geist (1980) Silicon Photodiode Absolute Spectral Re-
sponse Self-calibrating. Appl. opt. 19: 1214.
91
5. - DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
5.1 INTRODUCTION
93
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
reference to published values for the specific location and time or by refe-
rence to a long series of observed data, in which near-maximum values are
usually connected with clear skies (see for example Fritz. 1949). Kondratyev
(1969) has tabulated representative clear-sky values of global solar radia-
tion, expressing these as a function of latitude. Sample values are included
here as Table 5.1.
LatH ude. Jan. Feb. Marc. Apr. May June Suly Au-g. Sep I. Qcl. Nov. Decc.
de-s,·
; Nor Ih
90 0.0 0.0 0.2 14.0 30.7 36.6 33.3 18. 1 3.3 0.0 0'.0 0.0
85 0.0 0.0 1.0 14'.3 30.6 38·. 1 32.9 18.4 4.3 0.0 0.0 0.0
80 0.0 0.0 2.9 15. 1 30.1 35.4 32.2 18.7 6.0 0.. 6 0.0 O. O·
75 0 .. 0 O. 8 5.6 16.4 29.5 34.4· 31.0 19.4 8.2 1.9 0.0 0.0
70 0.0 2.2 8.5 18.4 28.8 33.0 29.9 20.5 10.6 3.8 0.7 0.0
65 1.0 3.9 11.3 20.4 28.7 32. 1 29.5 28.2 13.3 6.1 1 .9 0.3
80 2.5 6. 1 13.9 22.5 29.2 32.2 30.0 23.5 15.8 8.5 3.8 1.6
55 4.4 8.7 18.4 24.3 30.2 32.8 30.8 25.2 18.1 11.0 5.7 3.0
50 6.8 11 .5 18.7 26.0 31 . 1 33.3 31.7 26.8 20.2 13.8 8.1 5.8
45 9.4 14.5 16.6 27.4 31.9 33.8 32.1 . 28.3 22.2 14.4 10.9 8'.2
4·0 12.4 n:.2 23.0 28.5 32.4 33.7 33.0 29.0 23.9 18.5 13.6 11.1
35 15.0 19.6 24.8 29.4 32.6 33.8 33.1 30.1 25.4 20.6 16.0 13.7
30 17.5 21.7 28.2 30.0 32.8 33.3 32.9 30.8 26.8 22.6 18.4 16. 1
25 19.8 23.8 27.3 30.3 32.2 32.8 32.5 30.7 27.9 24.4 20.6 18.4
20 21.8 25.2 28.3 30.3 31.8 32.0 31.7 30.6 28.7 26.0 22.6 20.7
15 23-.7 26.6 29. 1 30. 1 30.8 30.9 30.8 30.3 29.4 27.2 24.4 22.6
10 2.5.4 27.8 29.7 29.8 29.7 29.5 29.6 29.8 29.8 28.2 26.0 24.6
5 27.1 28.7 30.1 29.4 28.5 28.0 28.3 29.0 29.9 29.1 27. 5 26.4
0 28.4 29.4 30.2 28.7 27.1 26.4 26.8 28.2 29.8 29.7 28.7 28.0
5 29.7 30.1 30.1 27.9 25.7 24.8 25.2 27. 1 29.3 30.1 29.8 29.6
10 30.8 30.5 29.6 26.9 2.4.2 22.8 23.5 25.7 28.6 30.1 30.6 31.0
15 31.8 30.8 29.0 25.7 22.5 21.0 21 .6 24.3 27.7 29.9 31.3 32.2
20 32.5 31.0 28.2 24.2 36.1 18.9 1-9.8 22.7 26.6 29.5 31.9 33.0
25 33.1 31.0 27.3 22.5 18.7 16.9 17.8 21.0 22.2 29.0 32.2 33.9
30 33.6 30.7 26.1 20.9 16.7 14.9 15.6' 19. 1 23.8 28.3 32.5 34.6
35 33.8 29.9 24.7 19. 1 14.6 12.7 13.5 17.0 22.3 27 .4 32.5 35.1
4'0 33.8 29.0 23. 1 17.1 12.3 10.3 11.3 14.9 20.6 26.3 32.1 35.4
45 33.7 27.9 21.3 15. 1 9-.9 7.7 8.8- 12.7 18,.7 24.9 31.5 35.5
50 33.3 26.6 19.4 12.9 7.6 5.3 6.4 10.3 16.8 23.4 30.7 35.2
55 32.5 25.0 17.3 10.5 5.3 3.2 4.2 7.7 14.5 21.6 30.0 34.6
60 3'1.7 23.4 15. 1 7.9 3.4 1 .4 2.2 5.3 12.0 19.8 29.5 34.3
65 31 .4 21.6 '12.7 5.4 1.5 0.1 0.6 3.2 9.4 17.8 29.2 34.5
70 31.7 20.0 10.2 3.2 0.3 0.0 0.0 1.4 7.0 16.0 29.4 35.0
75 32.6 18.9 8.0 1.3 0.0 0.0 0.0 0.2 4.9 14.6 29.9 35.8
80 33.8 17.9 6.0 0.0 0.0 0.0 0.0 0.0 2.9 13.6 30.8 36.6
85 34.6 17.4 4. 1 0.0 0.0 0.0 0.0 0.0 1.3 12.9 31.4 37. 3
"0. 35.0 17.3 2.4 0.0 0.0 0.0 O. o· 0.0 0.1 12.8 31.7 31.8
Sou-Ih
Under conditions with clear sky the ratio (Ed..ji I Eg..ji) will generally be
a smooth function of solar zenith angle allowing unusual measurements to be
readily recognized. However, if direct (S) and global (Eg..ji) data are avail-
able, a preferable approach is to verify the consistency of the data, using
the rela'tionship:
(5.1')
Differences should be less than 5 per cent after the appropriate corrections
have been ma.de to the diffuse reading (see 4.2.2.6). Schmid (1976) demon-
s:trates the effectiveness of this approach (Fig. 5.1).
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
- N
I
E
200
-
3:
UJ
0
lOO
I" 0
UJ
-100
0 5 10 15 20 25 ao as .0 u 10 15 10 15
Figure 5.1 Difference between global and diffuse solar radiation for a
horizontal surface when there is no bright sunshine. Data are plotted
against the solar elevation. (Schmid, 1976).
'" r=-+--!-+-~f-----"'----t---'--'--t-----j!
I I
150 I i "
!
- 125 ~_+__--_-+I----+------+-----+--I
~ tower .ffect
\.A
-"
3=
W
100 1==-=JJrf=====:::t====:t====+::::::=::=fc=t=""'1
~
75l--4--+__---_l_---'----+---____1----t___;
tower .ffect
50 ~-H-+__---_l_---_+---____1----t___;
25 ~..JI-~---_+_---_+---____1~.---t___;
10 11 12 13 HI
Figure 5.2 Temporal variation over one day of the ratio of -global radiation
determined by combining of Sand Ed to the global radiation Eq deter-
mined by the pyra nometer. (Wendler and Ea ton, 1980).
95
CHAPTER 5 ' DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
- cP.'9 0 ! - - - - - - 1 l - - - - - + - - - - - - - - - - - t
_
0, 80 I'-~-~_I+_----+------.--
"C:
CD· 7 0 1 ' - - - - - 1 + - - - - + - - - - - - - - - - ;
.a
_. 80
~. . =='N~====:::::::::::;;:::::::::::::;====i
t=:===j1~fi.=
9h
w: 50' 1'-'- - - - - - - = - - - - - - t , - - - - - - - - -
al); H,
Ffgure 5.3 Two hours of alhedo (upper curve), global (m.iddle curve) and
reflected (lower curve) radIation values. (Wendler, 1981.)...
are less than about 10 per cent (Dav.ies et al., 1970). E*' is the net and Et
the downward solar Eg and terrestrial El irradiance.
and e'nsuring that differences are less than about 10 per cent.
IS less than about 15 per cent (Davies et al.,. 1970). Wendler and Eaton
(1980) also utilized this me·thod.
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
(5.5)
Numerous studies (e.g. Liu and Jordan (1960), Collares-Pereira and Rabl
(1979), Erbs et a1. (1982)) have demonstrated a relationship between
Ed t I Eg {o and Eg {O, which has potential for the verification of diffuse radia-
tion data. However, on an individual hourly and even daily basis the scatter
in the relationship is' large and consequently this represents little more
than a crude verification method.
Tabulated values of albedo and the knowledge of the global solar radia-
tion would normally make it possible to verify the measured reflected solar
radiation. On the other hand, both the considerable variability in the albe-
do for a given cover type and the zenith angle dependence (see 3.8.2) suggest
97
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
Sections 3.6 and 3.7.2 describe equations which may be used to calculate
this irradiance. However, the accuracy of these equations is frequently not
sufficient to estimate Ee t.
(5.6)
have been described in the literature (e.g. Davies (1967), Viswanadham and
Ramanadham (1967)). Potentially these represent a viable verification tool
and they can be used to a·dvantage for time periods of a day or more. At
shorter time intervals the phase shift between the short- and long-wave quan-
ti ties could lead to invalid conclusions about the quality of the observed
data.
If the global and other component solar irradiances are known for a
horizontal surface, it is possible to perform some relatively rigorous veri-
fications of the slope radiation data. The individual components are con-
sidered in the following sections, except for the reflected solar and upward-
terrestrial irradiances, which may be treated in the same way as in the
earlier horizontal-surface cases.
98
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
~
:)
1.4
0
..c:
co
'E
1.2 • -E
..., 1.0
~
z 0.8 ·s
0
i= 0.6
N::~
ci
0
ci
a:
a:
0.4
0.2
t
ci
..J
0 0.0
en 0 30 60 90
SLOPE ANGLE
Figure 5.4 Measured global solar radiation for the hour ending 0800 Local
Apparent Time, June 11, 1980 at Vancouver Canada, Data for given slope
angles and orienta tions (E=East, W=West, S=South, N=North) under partly
cloudy skies. (Hay, pri va te com munication),
where
and s is the slope angle, 8 is the solar zenith angle and a and bare
the solar and slope azimuths, respectively (see Fig 5,5). It is desirable to
use values valid for an hourly interval or less due to non-linearities in the
above expressions. Agreement between calculated and observed values should
be within 2 and 5 per cent, with the higher accuracy being associated with
the use of a pyrheliometer to measure the direct radiation.
NORMAL TO THE
INCLINED PLANE
WEST_.
HORIZONTAL
SURFACE _ _ /
PROJECTION ON A HORIZONTAL
PLANE OF NORMAL TO THE
INCLINED SURFACE
99
CHAPTER 5' DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
5.0;r---..,......---r------.---,----,r---~,
3:0:
• .'
• ..
0.0.
0.0 l'.0 2:0 3,0 4.0 5~0
E~,1~' 91, 0.5 Et~ (i-cos s) + O.S 0-[;) Ett (l+cos s) + 0.5 Et,t(i-cos. s) C5.10)
100
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
The observed net radiation should agree with the sum of the observed or
calculated component fluxes, namely:
(5.11)
However, the accumulated effects of errors in each of these terms limit the
applicability of this relationship.
The most common and effective method for verifying the accuracy of spec-
tr~l measurements is based on the relatively conservative nature of the ratio
of the radiant intensities between two given spectral intervals. Wendler and
Eaton (1980) lean on this fact to assess the quality of data generated by
pyranometers equipped with OG1, RG2 and RG8 filters (see Fig. 5.7). Stewart
(1981) has shown that over short time-intervals numerically modelled values
may be used both for verifying the quality of data and replacing erroneous or
filling in missing data associated with instrument malfunction (Fig. 5.8).
Less success will be achieved with quality-control of the ultraviolet-radia-
tion data by the use of the ratio of UV to global solar radiation. The
wavelength dependence of the at mospheric at tenua tion processes results in
large variations for this ratio, precluding the application of strict quali-
ty-control limits (Wendler and Eaton, 1980).
15.0ct 1978
300 --~---- - ~ - ~ -
250
-C'l
I
E
200
-==
ISO
~ 100
W
50
01 IIH 09 10 11 17 13 11 16
Time Ch)
101
CHAPTER 5 DATA PROCESSING, VALIDATION,. AND QUALITY-CONTROL
£Total 'insolaUon
..0
..
.... -2
.......
·111
...
.
Time 1(ll1'h-12h) .
5.3.1 Interpolation
This procedure involves the use of data observed both before and after
the break in the record to estimate the missing data values. Many methods can
be used, the simplest being to use linearly weighted averages of the two
observations adjacent to the data gap. More sophisticated techniques make
use of advanced curve-fitting procedures (e.g. cubic splines) and a longer
data record on either side of the dis.continuity. All interpolation proce-
10~ ....
:5 06\
:.;i .-\
.~ D.6,.1-;-\-------------'1
·0 '
'u .0" I--'\-,-------------'1
.. '~02; ;.-1\__
(
.'....."
<: Flrst-orderregresslon
..< . . '
;
;
o .~- •
'0·2·~_""""'--2~-~.3--.•~-""'~.....J
;D1s.placement. days
Figure. 5. 9 Averag:e autocorrelation for U1'lited Kingdom global solar-radia-
tiondata. (BriIlkw<?rth, 1977).
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
dures invoke assumptions about the temporal changes in the variable under
consideration (e.g. that it varies linearly with time). However, only under
clear-sky or heavily overcast conditions does radiation have simple temporal
variations. Thus interpolation schemes must be chosen with reference to the
temporal characteristics of the radiation prevailing at the time of the break
in observations. The temporal variability (primarily the auto-correlation)
of the data series and the method used to approximate the time series (e.g.
autoregressive-moving average, maximum entropy, maximum likely hood models)
will also strongly influence the length of the data gap, which may be filled
by interpolation. Fig. 5.9 and the results of Armstrong (1980) and Valko
(1980) indicate that it is of very little value to use interpolation tech-
niques for time periods as long as a day. In an analysis of data from a
single site in continental North America (st. Paul, Minnesota) Skaggs et a1.
(1982) found that the characteristic time between uncorrelated observations
ranged between slightly more than one day to just over 3 days.
960
~)
..6f"1"'bb01.tOfd
~-AirpOfl
o I Km 10
l' , , !: : '
• BMihot
5.3.2 Calculations Using Variables Coincidently Observed vith the Data Break
103
CHAFTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
joint distribution, these being different variables observed at the same site
or identical variables observed at adjacent sites. Fig. 5.10 demonstrates
the high interstation correlations for solar radiation that exist in Vancou-
ver, Canada, and which have made it possible to estimate missing data for a.ny
of the station pairings (Hay, 1984).
Since archive.d data are the basis for all subsequent studies, including
those of both the originator and others who may may have access to the ar-
chive, it is necessary that the archive be accurate, accessible. and unambi-
g.UOliS. The quality-control procedures outlined in 5.2 will do much to ensure
IS I VALUE If.; I
U I I11 LJ
repeated J!' times.
~. VALUE ~.
. I I I I .•
repeated: 12' t; mes
Abbrev i at,i",,'
Fig'ure 5.11 Format of archived hourly, daily and monthly d.ata at the Cana-
dian At mosp:heric' Environm,ent. Service.
104
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
that the first objective is met. The other desirable attributes require that
the archive be organized in a systematic manner. An example of a convenient
format is shown in Fig. 5.11. The flag fields may be used to designate such
information as specific sensor and recording systems, techniques applied for
adjusting or interpolating data or to indicate the quality of the data in
case the quality has deteriorated through such factors, as, for instance,
moisture on the sensor. The archived programme of the World Radiation Data
Centre at the Main Geophysical Observatory, Leningrad, USSR (Fig. 5.12) is
associated with considerable expertise, Additional information on archiving
is contained in World Meteorological Organization <1960, 1981), In 1985 the
archiving at the WRDC will be fully computerized. Thus the data sent to the
WRDC may be delivered on standard magnetic tapes and the data sets can be
requested on tape.
REGION VII PEn10H VI EUROPE! EBPonA
CO\..OOn, iUlD
~D TlJRJ<n
CT!'A.... iIIl&.llt'I T1P'Jitq
nATIC:>HS
CTAHUM!<
U=ne
Kxpyxa
1AI.lal
J1,.190
Stockbola
CToJtronM "no"
.woe
1.000=0
JOIltaIlfO
8aaou.n
CawcfJI
btanbul
~T8»6y.l
4nkarli
AHRap8.
hair
H3l<Zp
Detle,,o,a T '
I ,80 '80 9? ''''2 . 872 920 Wl 868 -
2 "3 Z09 536 953 812 9}3 895
..... 1162 -
3 ,84 131 3"- 849 '90 8}1
,., 850 -
•5 158 16 68 553 244 111 190 -
54 25 86 soe 792 '/13 }64 805 -
6 ,08 '44 252 805 81. 573 ~oo 2.} -
7 50 '7 240 888 2'1/' 75 206 2}3 -
8 112 50 }1} 900 783 72 150 3'5 -
9 90 160 '4C 87' 759 139 1e:6 '64 -
10 ,01 ,6<; 2}<o
.,. 879 77' 21' ,"" 1t}1 -
11 90 '56 810 78' 609 }08 792 -
12
13
"15
94
6'
}2
50
.,
'69
""
&,
263
382
m
2,,",
806
875
263
8.18
630
7&
'86
56'
159
&50
001
858
'·9
507
696
760
525
10!-5
'079
849
-
-
-
-
16 &, '2. 22~ 795 '95 1b; ,80 ;05 -
17 32 79 '32 828 309 32} Ho 09} -
18 50 97 ''/1 185 &39 '08 ,7] 32' -
19 12 83 80Z 61~ ,.} .2} -
""25 Wl
-
20
21 }&
22
58
2"
205 ,.,
105 570
&10
20}
...'"
&<2
&&5
.,8
""" -
2'2
23
25
22
'01
&5
29Z
281
151
63}
584
501 44
-.0
699
158
699
--
2'
25
25 50 '5' - 544 521 113 800 -
-
26
22
22
}6
zz
58
03
259
15'
}58
Z9Z
196
190
8}<o
no .,.
9'lO
-
27 22 90 '22 80 68 151 8}5 651 -
28 22 '01 '48 '71 589 184 1161 590
,., -
29 2.2 25 ,"" 261 6'6 '52 '168 -
JO
31
22 22 2,2 208 621 190 }O} '59 -
IT
-
..
'969 2681 6829 20109 '58" 122&0 15119 '7}52
T 66 228 612 -
..
90 528 '09 50" 518
III
SS
}O
'.0 , .5
15·
2.5
1"-2-·
0.1
,.0
0.1
5'
,.1
>.
5''.8
82
2.1
115
}
Figure 5.12 Sample table for November 1983 from the published USSR radiation
data archive (USSR State Committee). Global radiation is given in
10- 2 MJm- 2 , and LT: 'monthly total, T: monthly meq.n, Lss: monthly total of
sunshine duration in hours, SS:, monthly mean.
105
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
not the case and Fig. 5.13 gives one example which indicates the departure
from normality. Therefore, the distributional characteristics are best de-
scribed using non-parametric statistics including the median, quartiles and
percentiles. Differences between the two measures of central tendency (the
mean and median) indicate the skewness shown in Fig. 5.15 which is a typical,
feature of the radiation data.
b) 05
a)
OA ~---.-----r----y----' OA Vancouver, January. 1977
Vancouver, U. B. C., K160 degrees, /
0.3 0.3
South facing KI
0.1 0.1
0.0 0.0
>-
u
~
::J
ow
"C i f~r;u: ."~ I
0.3
0.2
0.1
~_
. " ]
. "/
: 0.0 _.
0: 0.0
u.
::C=.-I ~ ~~lt:::
::.L.'.· ~~L:"
0.0
0.1. z . - .
O.Ci.obeo
. 0.1
/
,'. .
I
1
0.0 '. . . '.
0.0 1.0 2.0 3.0 , 4.0 0.00.0 1.0 20 3:0 4,0
HOURLY SOLAR RADIATION CMJ m· 1 hour·) HOURLY SOLAR RADIATION (MJ ni ., hOUf 1 1
Figure 5.13 Frequency dist'ributions 'of observed hourly values of global so-
lar radiation at Vancouver, Canada for a) the horizontal and three
south-facing surfaces in January 1977 and b) the 60-degree south-facing
surface in January, April, July and October, 1977. (Hay, 1979a.).
100
~ 80
>-
u
z
w
0 60
w
a:
u.
w
> 40
i=
<{
..J
:;)
~ 20
:;)
u
2 3 4 6
106
CHAPTER 5 DATA PROCESSING, VALIDATION; AND QUALITY-CONTROL
VANCOUVER January
Figure 5.15 Percentiles and quartiles of measured global solar radiation for
a horizontal surface at Vancouver in January for the period 1959-1975
and for averaging times of one to ten days. (Hay, 1979a).
Before some uses of the solar energy (e.g. solar heating) are possible,
the radiant intensity must be above a minimum level. In many cases the mini-
mum threshold is also a function of temperature. For this reason tabulations
which consider radiative-energy inputs above specified thresholds and for
specified temperatures can be of immense use. Table 5.2 is one such example
and provides the mean monthly value of global solar radiation, which is
greater than the energy-threshold level given at the top of the column and
which was received while the ambient temperature was within the temperature
range noted at the left of the table. To illustrate this, the average values
in January were 27.37 MJ m- 2 of global solar radiation above a threshold
intensity of 300 W m- 2 that occurred at the ambient temperature between -15
and -10 degrees Celsius.
107
CHAPTER 5 DATA PROCESSING~ VALIDATION, AND QUALITY-tONTR6L
-50 V-'-"T
. 1.
10
'/'
1"'\1=l~UIY"
!50'/' 1,10'10 . '
0,0 I, -11-~-4----1
-, I ,
X / Stintis I ]
5,0 J, ,~L:'/_+
C-~ I
/ +
I
-+-;I
--1
','I'
I,
o
.. oJ?I.
U
10
"
~
~ 0
.o
CL
I
r:
. -+---=-/-+--.I-
E
t-
15,0
3.9,0";-,..L-l-l-"'-!-'---J'--"Y'If---.'L.'L... 1-'
o 8 12 16
•• +' • ~
Table 5.2 An example of the analysis of radiation data for selected thresh-
olds of global solar radiation intensity and air-temperature classes.
(McLaren et aL. 1979).
SUDBURY (OMT> (60HUO)
Month: .Jan-
-3S
30.-: 3S
2S-: 30
20-: 2S
IS-: 2D-
to-: IS
S.: 10
0.02
0.-: ~ 4.43 4.43 1.24 0.36 0.13 0.05
-5+:
-10+: -S
0 n."
38.19
22.6'-
38.19
10.82
25.08
5.94
17,28
3.47
11.93
2.03
7.94
1.13
4.94
0.54
2.65
0.21
1.14
0.07
0.25
0.01
51. 57 38.13 27.37 \9.03 12.48 7.15 3.37 1.08 0.14 0.01
-15+ : -10 68.89 68.89
82.51 64.13 49.49 36.82 26.35 17.59 10.11 4.75 1.37 0.15 0.00
-20+ : -IS 82.51
57.82 47.24 37.69 28.75 20.93 14.35 8.74 4.45 1.39 0.07
-25+ : -20 57.82
37.28 37.28 30.79 24.81 \9;08 13.90 9.81 6.16 3.40 1.28 0.11
-25-
228.40 228.40 167 • .,. 1\7. SI 74.93 47.88 25.10 9.64 I. 85 0.02
TolaJ 90 PCTL \,16 0.01
Hean 192.93 191.93 137.62 94.93 61.42 36.87 18.3£ 6.88
Direct
10 PCTL 135.57 135.51 91.61 61.47 39.09 22.81 11. 09 4.37 0.71
Rad
(HJ/Sq "etre)
293.88 221. 67 17'.24 153." 120.40 '1.91 50. " 10.n 1.34
Tot_J 90 pcn IOU.'" 377.82 2S•• " , .80 0.2'
336.86 248.23 200.68 170.17 134 .04 114.13 '2.32 £2.30 36.84
De9ree Hean 964.20 2.53
889.74 306.83 204.23 152.24 127.93 103.29 82.3S 60.79 39." 23.50
Day. 10 PCTL
108
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
N = e
(5.12)
Figure 5.17 Map of the southern United states showing the number of years of
observation required for an estimate within 5 per cent of the true mean
95 per cent of the time. (Based on data for July). (Black, 1961).
Table 5.3 Sample size needed to estimate the true monthly mean solar radia-
tion within 1 MJm- 2 day-t at a confidence level of 95 per cent. "Random"
means the persistence has been removed from the data, "Real" means the
real data has been used.(after Skaggset a1.1982)
109
CHAPTER 5 DATA PROCESSING; ·VALIDATION, AND QlIALITY-CON'fROL
In many cases the radiation data are needed for a site other than the
measurement location. This requires a knowledge of the spatial coherency of
the radiation field. Variations in this coherency will have an influence on
the distance for which data may be extrapolated or interpolated, and hence on
the density requirements for an observational network. The spatial coherency
is a function of the integration period of observation (e.g. hour, day,
month) as well as of locational and seasonal factors. These effects have
been demonstrated in a number of papers including those of the World Meteoro-
logical Organization, 1969 and 1981; Gandin, 1970; Alaka, 1970; Pivovarova,
1978; Wilson and Petzold, 1976; Hay and Suckling, 1979; Latimer and Won,
1981; Pivovarova and Stadnik, 1969 and Hay, 1984. The results presented in
Figs. 5.18a, band 5.19 are from meso- and macroscale networks and show that
an increase in observational integration period increases the spatial cohe-
rency and hence decreases interpolation and extrapolation errors. Fig. 5.20
illustrates the change in extrapolation error, as distance increases, while
Figs. 5.21 and 5.22 provide an indication of the considerable impact that
topography has on the coherency of the solar radiation field. Distance-
dependent correlation variations such as those discussed above, are an essen
I- I-
zQ Zc::J
W •
W •
&.;.-...t ..... a) .-. .... I f , __ b)
U
.-. . I f 111
U
.-. .,lJlf"II . . **
IfIlfs
IJ...
IJ...
.~.
•
5
III • .. IJ...
IJ... -+ t-
5 , "If
11 i'J'-ltw
W
0
..1II;ifl lf ,
11
W
0 -+
t-
T • If
't'+
III It
-+ 11 If
~
t.-~ u~ T+
u~ + t- Will 111 11
z::Q III~ •
zc::J t-
o
.-. TT t- t- + + -+
0
.-.
T
I- + I-
a: a:
.-J .-J
WlQ WlQ
0:: • 0::- .•
0::0
o 0.0 20.0 40,0 60,0 80.0 ~';Jo.o 20.0 40.0 60.0 80.0
u STRT ION SEPRRRT10N LK/1) Ll STRTlON SEP8RRTION (KM)
Figure 5.18 a) Correlation 'of hourly global solar-radiation data for indi-
vidual station pairs versus the distance between that station pair. Data
involving Grouse Mountain (see Fig. 5.8) shown by (+). All other station
pairs shown by (*). Data from June 1979 to May 1980. b) Same as a) but
for daily data.
,,
,
- ' .... ....
..........
........
-"---3
0.6
0.4
0.2
Figure 5.19 Distance (::orrelation function for (1)" monthly and for (2) daily
global solar radiation for Jllne-Julyand for (3) samples for the entire
year. From World Meteorological Orga niza tion, 1981.
110
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
-.'" 8
...
•
'\lI
~
.,E 5
/'
-
::E
Z
0
4
.
,,
I
I
I "
j: I
~
3 I
t
> /
w I
C
2
c ,
a:
« 1
,
I
C t
Z
«
I- 0
f/l
0 200 400 800 800 1000 1200
DISTANCE BETWEEN STATIONS (km)
Figure 5.20 The standard deviation of the daily differences in' global solar
radiation forgiven station pairs in British Columbia and Alberta, Cana-
da, as a function of the distance between the station pair. Broken line
for data from Southern Canada by Wilson and Petzold, 1972. (Suckling and
Hay, 1976).
Strait
of
Georgia
•: ' , •
, ';'
It"' 10
Pit"
111
CH~PTER 5 DATA PROCESSING, VALIDATION, AND QUAL.ITY-CONTROL
; r
ij
U 0.7
CQHIS.WGadt
'_n 1
Wo,,~""
1 ---==1: - - IL-
• Central Swit.",land
--II_. -J
o 50 100 ISO DJ 2SO 300
Stollon •• porollan di.lonce 0 (km!
Figure 5.22 Correlation coefficients for monthly sunshine totals for given
station pairs in Switzerland. For some' categories only the regression
curves are drawn. Also shown (dashed' .line) is the regression curve for
Southern -Canadian global solar -radiation "data (see Fig. 5.21). (Valko,
1980). .
Table 5.4 Errors (E) associated with the spatial interpolation of hourly and
daily values of global solar" radiation at Vancouver, Canada. Based on
the technique of Gandin (1970) and Alaka (1970).. Interpolation error (E)
in per cent of mean, kJ m- 2 h- 1 (for hourly data) and MJ m- 2 day-1 (for
daily da tal, .
e--
". -----
a'
Figure 5.23 The relationship between the optimum interpolation error in the
. centre of a square and of a triangle and distance for different values
of the observation error. (Gandin, 1970).
112
CHAPTER 5 DATA PROCESSfNG, ",VALIDATION, AND QUALITY-C'ONTROL
Table 5.5 RMS errors of linear interpolation between anomalies of daily ra-
diation totals for given mean distances between stations (Pivovarova,
1978).
Component Distance in km
January July
Units: Kilometres
113
,DA'FA P'ROCE&SDfG, llALIDA'tlON,' ANI) OUALITY-~OlfTROt. '
5.8 COIfCLUSIOltS
, REFERENCES
Arnfield, A.J. (1975) A note on the diurnal, latitudinal and se<jlsonal varia-
tion of the surface reflection coefficient. .1;, AjJp.l. He{. 14: 1603.
Atwater.' M.A. and J.T. Ball (1978) A nUinerical solar ra:diat'ioni model based
on standard meteorological observations. Sohfr-Bnerqy 21: ,163.
Bahm, R.J. (1978:) Reg,iqnaL 'differences in' solar radiation availability. So-
lar Age, June 1978, p. 34.
Berland, T.G. (1962) Geographical principles of the: solar radiation regime.
Proceedings All-Union Scientific Meteorology Conference, 1961: Vol. IV.
114
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
Black, J.N. (1961) Some aspects of the climatology of solar radiation. Paper
E/CONF. 35/S13, U.N. Conference on New Sources of Energy, Rome, p. 311.
Budyko, M.I. and O.A. Drozdov (1966) On the use of means in climatological
studies. Het. and G.idr. 10: 3.
Da vies, J .A. (1981) Models for estim a ting incoming solar irradiance. Unpub-
lished Manuscript. Canadian Climate Centre Report No. 81-2, Atmospheric
Environment Service, Downsview, Ontario, p. 120.
Davies, J.A. and J.E. Hay (1980) Calculation of solar radiation incident on
a horizontal surface. In Proceedings First Canadian Solar Radiation
Data Workshop, J.E. Hay and T.K. Won eds., Canadian Atmospheric Environ-
ment Service, Downsview, Ontario, p. 32.
Davies, J.A., P.J. Robinson and M. Nunez (1970) Radiation Measurement Over
Lake Ontario and the Determination of Emissivity. First Report, March
1970, Department of Geography, McMaster University, Hamilton, Ontario,
Ca nada, p. SS.
Da vies, J .A., W. Schertzer and M. Nunez (1975) . Estimating global solar ra-
diation. Boundary layer Heteoro.loqy 9: 33.
Erbs, D.G., S.A. Klein and J.A. Duffie (1982) Estimation of the diffuse
radiation fraction for hourly, daily and monthly average global ra-
diation. So.lar £nerqy 28: 293.
Foulkes, C.W. (1980) Controlling the accuracy of solar radiation data from a
-low cost, 30 station network. Proceedings of the Sixth National Passive
Solar Conference, American Section, International Solar Energy Society,
Newark, Delaware, p. 377.
Hay, J.E. (1976) A rev'ised method for determining the direct and diffuse
components of the total short wave radiation. Atmosphere 14: 278.
Hay, J.E. (1979b) An analysis of solar radiation data for selected loca-
tions in Canada. Climatological Studies No. 32, Atmospheric Environment
Service, Downview, Ontario, p. 158.
115
~
/ -
CHAPI:ER 5 DATA PROCESSING, VALIDATION, AND QUALIT'l-CONTRftL
Hay, J.B. and J.A. Davies(1980) Calculation of solar radiation incident on-
an inclined surface. In Proceedings of the First Canadian Solar Radia-
tion Data Workshop, J.E. Hay and T.K. Won, eds., Canadian Atmospheric
Environment Service, Downsview, Ontario, p. 59.
Hay, J.E. and P.W. Suckling (1979) An assessment of the networks for
measuring and modelling solar radiation in British _Columbia and adjacent
areas of Western Canada. Canad.ian G'eoqrapker 23: 222.
Hoyt, D.V. (1978) A model for the calculation of solar global insolation.
So.lar £nerqy 21: 27.
Kimball, H.H. (1928) Amount of solar radiation that reaches the surface of
the earth on the land and on the sea, and the methods by which it is
measured. Honth.ly fiea ther Rev. 56: 393.
Liu, B.Y.H. and R.C. Jordan (1960) The interrelations ship and characteristic
distribution of direct, diffuse and total solar radiation. So.lar £nerqy
4: 1. -
McLaren, J.F., J.E. Hay and J.A. Davies, (1979) Define, Develope and Estab-
lisha Merged Solar and Meteorological Computer Data Base. Report to the
Atmospheric Environment Service, Downsview, Ontario, p. 250. '~
Machta, L. (1980) -An ,approximate method for quality control-of solar radia-
tion in'struments. In an Introduction to Meteorological Measurements and
Data Handling for Solar Energy Applications; International Energy Agen-
cy, Task IV, United States Department of Energy, DOE/ER-0084, Washing-
ton, D.C., A-II-1 to A-II-28. -
Pal tridge, G. W. and C.M .R. PIa t t (1976) Radia ti ve Processes in Meteorology
and Climatology, Elsevier, Amsterdam, p. 318.
Pivovarov..a, Z.I. (1978) Accuracy of solar radiation data and density of net-
work of solar radiation stations. Report to the World Meteorological
Organization, Solar Energy Meeting, Geneva, October 1978, Appendix F
(b), WMO, Geneva, p. 7.
116
CHAPTER 5 DATA PROCESSING, VALIDATION, AND QUALITY-CONTROL
Skaggs, R.H., D.G. Baker and J.E. Ljungkull (198~) The influence of per si-
stance and variability on the required solar radiation record length.
Solar £nerqy 28: 281.
Suckling, P.W. and J.E. Hay (1976) The spatial variability of daily values of
solar radiation for British Columbia and Alberta, Canada. Climatological
Bulletin, No. 20, McGill University, Montreal, Quebec, p. 1.
USSR State Committee for Hydrometeorology and Control of the Natural En-
vironment: Solar Radiation and Radiation Balance Data (The World Net-
work), Leningrad.
Wilson, R.G. and G.E. Petzold (1976) A solar radiation network evaluation for
the Churchill Rive'r Basin. Papers in Climatology - The Cam AlIen Memo-
rial Volume, Department of Geography, McMaster University, Hamilton,
Ontario, p. 84.
Wolfe, W.J., G. Clark and J.E. Rudzki (1980) On-site quality control of
insolation measurements. Proceedings of the American Society Mechanical
Engineers, 1980 Winter Annual Meeting, Chicago, Illinois, p. 1.
117
A.1. EXTRATERRESTRIAL SOLAR SPECTRUM
200 .. 310 nm: G.Brasseur and P.C.Simon (1981) L~ Ceophys. Res. 86: 7343
310 •. 330 nm: J.C.Arvesen, R.N.Griffin and B.D.Pearson (1969) Appl. opt.
8: 2215
330 .. 869 nm: H. Neckel an D.Labs (1984) Solar PhysJcs 90: 205
870 .. IR nm: E.V.P.Smith and D.M.Gottlieb (1974) Space Sc.i.Rev. 16: 771
At the wavelength limits between the different spectra the values are forced
to be equal. After these adjustments the total spectrum is scaled in order to
yield a solar constant of 1367 Wm- 2 •
A(nm) A(nm)
119
APPENDIX 1 EXTRATERRESTRIAL SOLAR SPECTRUM
A(nm) SA (Wm- 2 nm- t ) So -' A(Wm- 2 ) A (nm) . SA (Wm- 2 nm- t:) So _A(Wm- 2 )
120
APPENDIX 1 EXTRATERRESTRIAL SOLAR SPECTRUM
121
APPENDIX 1 EXTRATERRESTRIAL SOLAR SPECTRUM
12,2
APPENDIX r EXTRATERRESTRIAL SOLAR SPECTRUM
123
APPENDIX 1 EXTRATERRESTRIAL SOLAR SPECTRUM
124
APPENDIX 1 EXTRATERRESTRIAL SOLAR SPECTRUM
125
APPENPIX 1 EXTRATERRESTRIAL SOLAR SPECTRUM
126
A.2. NOMENCLATURE OF RADIATION QUANTITIES
cr
8
• J
411'-0
Y (fl ,Il ) dO
8 i 8 •
[1I- i ) P(ll i ) - f Yr(lli'rlr) cos"r d!l r T(ll ) -
i f Yt (Ili ,ll t )C08 "tdll t
2TT ·Il 211'rl
0 0 0
J (H )
s s
• ..!..
0 f Ys(lli,lls)dEr'lli) J (Il )
r- r
. f y/rli,rlr)cos"idF.i(rl i ) Jt(rl ) •
t f y (ll i ,llt)cos"idEi (P. )
t i
s 4T1'rl 21T·1l 2T1 ·Il
0 0 0
[\01 -2
m sr-i) [\01 -2 -h
m sr J [\01 m-2 sr -i)
directional reflectance factor directional transmittance factor
for hemispherical incide~ for hemispherical incidtmc-e- -
P(ll ) •
r
J (ll )'TI'rl
r r
Et
o T(rl )
t
. Jt~llt)'1T'rlo
E
i
where Ei - f
21T'1l
cos "idEi (rl i ) [\01 m-2)
127
APPENDIX 2 NOMENCLATURE OF RADIATION QUANTITIES
radiant
energy Q. (W) J . W 8 ,45-05-130
radiant
flux tIl'. (1') w tI>'. =-
d
dt
Cl
pow-er 45-05-135
-1
2 The radiance is a
-2 d t1>
radiance L Wm, sr L • ""'dO-=-'',---=d;"'A--=-,c-o-s~" . conservative quantity fn
an optical system
Me:
'emittance £
£ - M £ - 1 for a black body
£... 1
~.r and
. ~i are the
If!
r
reflectance P • f"""" reflected: and' incident
1 radiant flux,
respectively
129
A.3. CRITERIA FOR THE CHOICE OF INSTRUMENTS
A.J.t Accuracy
131
APPENDIX 3 CRITERIA FOR THE CHOICE OF INSTRUMENTS
A.3.7Cost
The cost factor tends to be greatly overemphasized in the choice of
instruments. The purchase cost is usually so s,ma.ll relative to the. costoE
long term operations that it pales to insignificance in the overall program.
Only for networks of instruments does the first cost assume a major signifi-
cance, and,even then, economies in the cost of basic sensors are likely to
be mistaken economies. It is normally good judgment to g.et the best quality
instrument available for a particular task, regardless of the first cost of
the device. .
CM ,Cs ,Ca single particle cross section for extinction, scattering and
absorbtion respectively
E· irradiance
E* net radiation
J source function
K degrees Kelvin
L radiance
P degree of polarization
d Julian day
a absorption coefficient
133
APPENDIX 4 GLOSSARY OF SYMBOLS
A electromagnetic wavelength
p reflectance
P6 ground albedo
a standard deviation
transmittance
latitude
134
A.5. CONSTANTS AND CONVERSION FACTORS
135
INDEX
A
Angstrom
pyrheliometer, 54f
scale, 52
A
Absolute radiometers, 53
active, 58
passive, 58
Absorption, 11, 18
cloud, 36
coefficient, 11, 18
Actinometer, 59
Aerosols
absorption, 8
extinction, 9, 27, 28
scattering, 8
Air mass, 23f
Albedo, 17
cloud, 34f
measurements, 85
Rayleigh scattering, 29
sea ice, 39
sea surface, 42
snow, 38
soils, 41
solar elevation, 37
vegetated surfaces, 39
Atmospheric constituents, 7
B
Black body radiance, 33
Bouguer - Lambert law, 11, 13
C
Carbon dioxide, 7, 33
Cavity receiver, 56
Characterization of a radiometer, 57
Clouds, 10
solar radiation, 34
albedo, 34
thermal radiation, 36
emittance, 35
Criteria for instrument selection, 131
D
Data gaps, 102f
Data statistics, 105f
frequency distributions, 107
threshold analyses, 107
Declination, 23
Detectors
quantum, 51
thermal, 49
Diffuse radiation, 29
measurements, 71
137
INDEX
E
Emittance, 11, 33, 43
measurements, 86
Equation of time, 24
Equatorial mount. 61
Extinction coefficient, 11, 15
F
Filter factor, 79
G
Gaseous absorption, 28
Global radiation, 29
measurements, 62
I
Index of refraction, 12
International Pyrheliometric Scale, 52
H
Methane-, 8
N
Net radiation, 5
measurements, 74
o
Optical thickness, 12, 28
Oxygen, 8
Ozone, 8, 26
p
Phase function, 13, 15
Planck function, 18
Polarization, 16
aerosols, 18, 32
degree of, 16f, 31
measurements, 87
surface reflection, 43
Pyranometer, 62
calibration, 70
cut-off glass filters, 79
installation, 69
maintenance, 69
performance, 66
self-calibrating, 73
shading disk, 72
shadow ring, 71
silicon detectors, 73
Pyrgeometer, 74, 76
Pyrheliometer, 54
Angstrom, 54f
maintenance, 62
operational, 58
self-calibrating, 56
Silver~Disk, 55f
spectral measurements, 78
Pyrradiometer, 74
calibration, 77
maintenance, 77
138
INDEX
Q
Quality-control. 93ff
broad-band spectral data, 101
consistency checks based on calculations, 97
consistency checks based on observations, 93
consistency checks for inclined-surface irradiances, 98
R
Radiation, 74
budget, 3, 21
diffuse, 29
global, 29
solar, 8
thermal, 8
scales, 51
transfer, 12
Rayleigh-Scattering, 14, 17, 26, 28
Reflectance, 11, 13
mineral surfaces, 43
natural, 37
soils, 41
surfaces, 37, 44
Representativness of data
spatial, 110f
temporal, 109f
S
Scattering, 13
coefficients, 11
efficiency, 14
indicatrix, 13
Schott glass filters, 78
Silver-Disk pyrheliometer, 55£
Single scattering albedo, 12
Size distribution, 9, 14
cloud droplets, 10
aerosols, 8
Smithsonian scale, 52
Solar constant, 4, 23
Solar radiation, 7
clear skies, 26
clouds, 34
daily total input, 25
inclin&d surfaces, 73
measurement, 54, 73, 77
spectral distribution; 24, 28, 31, 119ff
Source function, 13
stratosphere, 8
Sun - earth distance, 23
Sunphotometer, 77
calibration, 83
recommended wavelengths, 79
T
Terminology,S, 127ff
radiative quantities, 11
radiative transfer, 11
Thermal radiation, 4, 7
clear skies, 33
clouds, 36
measurement, 74
Tracker, 61
Transmittance, 11
Troposphere, 7
139
INDEX
Turbidity, 27
coefficients, 28
If
Ultraviolet radiation, 26
measurem'ent, 8'4
V
View-limiting geometry, 55
W
Water vapour, 7, 33
absorptance, 2,6
World Radiometric Reference, 53f
140