Ch02 Radiat

Copyright 2015 by Roland Stull.
Practical Meteorology: An Algebra-based Survey of Atmospheric Science
2 Solar & Infrared Radiation
Contents Solar energy powers the atmosphere. This ener-

gy warms the air and drives the air motion you feel
Orbital Factors 27 as winds. The seasonal distribution of this energy
Planetary Orbits 27 depends on the orbital characteristics of the Earth
Orbit of the Earth 27 around the sun.
Seasonal Effects 30 The Earths rotation about its axis causes a daily
Daily Effects 32
cycle of sunrise, increasing solar radiation until so-
Sunrise, Sunset & Twilight 33
lar noon, then decreasing solar radiation, and finally
Flux 34 sunset. Some of this solar radiation is absorbed at
Radiation principles 36 the Earths surface, and provides the energy for pho-
Propagation 36 tosynthesis and life.
Emission 36 Downward infrared (IR) radiation from the at-
Distribution 39
mosphere to the Earth is usually slightly less than
Average Daily Insolation 40
Absorption, Reflection & Transmission 41
upward IR radiation from the Earth, causing net
Beers Law 43 cooling at the Earths surface both day and night.
The combination of daytime solar heating and con-
Surface Radiation Budget 44
Solar 44
tinuous IR cooling yields a diurnal (daily) cycle of
Longwave (IR) 45 net radiation.
Net Radiation 45
Actinometers 45
Review 47
Homework Exercises 47 Orbital Factors
Broaden Knowledge & Comprehension 47
Apply 48
Evaluate & Analyze 50
Planetary Orbits
Synthesize 51 Johannes Kepler, the 17th century astronomer,
discovered that planets in the solar system have el-
liptical orbits around the sun. For most planets in
the solar system, the eccentricity (deviation from
circular) is relatively small, meaning the orbits are
nearly circular. For circular orbits, he also found
that the time period Y of each orbit is related to the
distance R of the planet from the sun by:
Y = a1 R 3/2 (2.1)
Parameter a1 0.1996 d(Gm)3/2, where d is Earth

days, and Gm is gigameters = 106 km.
Figs. 2.1a & b show the orbital periods vs. dis-
tances for the planets in our solar system. These
figures show the duration of a year for each planet,
which affect the seasons experienced on the planet.
Practical Meteorology: An Algebra-based Survey

Orbit of the Earth
of Atmospheric Science by Roland Stull is licensed The Earth and the moon rotate with a sidereal
under a Creative Commons Attribution-NonCom- (relative to the stars) period of 27.32 days around
mercial-ShareAlike 4.0 International License. View this license at
http://creativecommons.org/licenses/by-nc-sa/4.0/ . This work is their common center of gravity, called the Earth-
available at http://www.eos.ubc.ca/books/Practical_Meteorology/ . moon barycenter. (Relative to the moving Earth,
27
28 chapter 2 Solar & Infrared Radiation

1MVUP &SJT
(a) (b) 1MVUP
.FSDVSZ 7FOVT &BSUI .BST /FQUVOF
6SBOVT
: &BSUI%BZT
/FQUVOF 4BUVSO

+VQJUFS
: &BSUI%BZT

.BST
6SBOVT
&BSUI
+VQJUFS
7FOVT
4BUVSO
.FSDVSZ

3 (N

3 (N
Figure 2.1 a (linear) & b (log-log)

Planetary orbital periods versus distance from sun. Eris (136199) and Pluto (134340) are dwarf planets. Eris, larger than Pluto, has
a very elliptical orbit. (a) Linear graph. (b) Log-log graph. (See Appendix A for comparison of various graph formats.)
the time between new moons is 29.5 days.) Because

Sample Application
the mass of the moon (7.35x1022 kg) is only 1.23% of
Verify that eq. (2.1) gives the correct orbital period
of one Earth year. the mass of the Earth (Earth mass is 5.9726x1024 kg),
the barycenter is much closer to the center of the
Find the Answer: Earth than to the center of the moon. This barycen-
Given: R = 149.6 Gm avg. distance sun to Earth. ter is 4671 km from the center of the Earth, which is
Find: Y = ? days, the orbital period for Earth below the Earths surface (Earth radius is 6371 km).
To a first approximation, the Earth-moon bary-
Use eq. (2.1): center orbits around the sun in an elliptical orbit
Y = (0.1996 d(Gm)3/2) [(149.6 Gm)1.5] (Fig. 2.2, thin-line ellipse) with sidereal period of P
= 365.2 days. = 365.256363 days. Length of the semi-major axis
(half of the longest axis) of the ellipse is a = 149.598
Check: Units OK. Sketch OK. Almost 1 yr. Gm. This is almost equal to an astronomical unit
Exposition: In 365.0 days, the Earth does not quite (au), where 1 au = 149,597,870.691 km.
finish a complete orbit. After four years this shortfall
Semi-minor axis (half the shortest axis) length
accumulates to nearly a day, which we correct using a
is b = 149.090 Gm. The center of the sun is at one of
leap year with an extra day (see Chapter 1).
the foci of the ellipse, and half the distance between
the two foci is c = 2.5 Gm, where a2 = b2 + c2. The or-
bit is close to circular, with an eccentricity of only
about e c/a = 0.0167 (a circle has zero eccentricity).
NPPO The closest distance (perihelion) along the ma-
FBSUI jor axis between the Earth and sun is a c = 146.96
Gm and occurs at about dp 4 January. The farthest
C 3 distance (aphelion) is a + c = 151.96 Gm and oc-
TVO O QFSJIFMJPO curs at about 5 July. The dates for the perihelion
BQIFMJPO D and aphelion jump a day or two from year to year
B because the orbital period is not exactly 365 days.
Figs. 2.3 show these dates at the prime meridian
(Greenwich), but the dates will be slightly different
in your own time zone. Also, the dates of the peri-
helion and aphelion gradually become later by 1 day
every 58 years, due to precession (shifting of the
Figure 2.2
Geometry of the Earths orbit, as viewed from above Earths
location of the major and minor axes) of the Earths
North Pole. (Not to scale.) Dark wavy line traces the Earth- orbit around the sun (see the Climate chapter).
center path, while the thin smooth ellipse traces the barycenter
path.
R. Stull Practical Meteorology 29

%BZJO+BOVBSZ 65$
B
1FSJIFMJPO C
"QIFMJPO

%BZJO+VMZ 65$

:FBS :FBS
Figure 2.3
Dates (UTC) of the (a) perihelion and (b) aphelion, & their
Because the Earth is rotating around the Earth- trends (thick line). From the US Naval Observatory. http://
moon barycenter while this barycenter is revolv- aa.usno.navy.mil/data/docs/EarthSeasons.php
ing around the sun, the location of the center of the
Earth traces a slightly wiggly path as it orbits the
sun. This path is exaggerated in Fig. 2.2 (thick line). Sample Application()
Define a relative Julian Day, d, as the day of the Use a spreadsheet to find the true anomaly and
year. For example, 15 January corresponds to d = sun-Earth distance for several days during the year.
15. For 5 February, d = 36 (which includes 31 days in
January plus 5 days in February). Find the Answer
The angle at the sun between the perihelion and Given: dp = 4 Jan. P = 365.25 days.
the location of the Earth (actually to the Earth-moon Find: = ? and R = ? Gm.
barycenter) is called the true anomaly (see Fig.
Sketch: (same as Fig 2.2)
2.2). This angle increases during the year as the day
For example, for 15 Feb, d = 46
d increases from the perihelion day (about dp = 4; Use eq. (2.2):
namely, about 4 January). According to Keplers M = (2 3.14159)(464)/ 365.256363 = 0.722 radians
second law, the angle increases more slowly when Use eq. (2.3b):
the Earth is further from the sun, such that a line =0.722 +0.0333988sin(0.722) + 0.0003486sin(1.444)
connecting the Earth and the sun will sweep out + 0.000005sin(2.166) = 0.745 radians
equal areas in equal time intervals. Use eq. (2.4):
An angle called the mean anomaly M is a sim- R = (149.457Gm)(1 0.01672)/[1+0.0167cos(0.745)]
ple, but good approximation to . It is defined by: = (149.457Gm) 0.99972 / 1.012527 = 147.60 Gm
d dp Repeating this calculation on a spreadsheet for several

(2.2)
M = C days, and comparing M vs. and R(M) vs. R() , gives:
P
Date d M R(M) R()
where P = 365.256363 days is the (sidereal) orbital (rad) (rad) (Gm) (Gm)
period and C is the angle of a full circle (C = 2 ra- 4 Jan 4 0 0 146.96 146.96
dians = 360. Use whichever is appropriate for your 18 Jan 18 0.241 0.249 147.03 147.04
calculator, spreadsheet, or computer program.) 1 Feb 32 0.482 0.497 147.24 147.25
Because the Earths orbit is nearly circular, M. 15 Feb 46 0.722 0.745 147.57 147.60
A better approximation to the true anomaly for the 1 Mar 60 0.963 0.991 148.00 148.06
elliptical Earth orbit is 15 Mar 74 1.204 1.236 148.53 148.60
29 Mar 88 1.445 1.487 149.10 149.18
M + [2 e (e 3 / 4)]sin( M ) + [(5 / 4)e 2 ]sin(2 M ) 12 Apr 102 1.686 1.719 149.70 149.78
26 Apr 116 1.927 1.958 150.29 150.36
+[(13 / 12)e 3 ]sin(3 M )
(2.3a) 21 Jun 172 2.890 2.898 151.87 151.88
or 23 Sep 266 4.507 4.474 149.93 150.01
22 Dec 356 6.055 6.047 147.02 147.03
M + 0.0333988 sin( M ) + 0.0003486 sin(2 M )
+ 0.0000050 sin(3 M ) (2.3b) Check: Units OK. Physics OK.
Exposition: Because M and are nearly equal, you
for both and M in radians, and e = 0.0167. can use M instead of in eq. (2.4), with good accuracy.
/PSUI1PMF WFSOBMFRVJOPY
TVNNFSTPMTUJDF
.BSDI
+VOF
"SDUJD$
JSDMF
BQIFMJPO FBSUI QFSJIFMJPO
5SPQJD +VMZ
TVO +BOVBSZ
PG$B
ODFS
'S XJOUFSTPMTUJDF
DFOUFSPGTVOUP BVUVNOBMFRVJOPY %FDFNCFS
&RVBUP
S DFOUFSPGFBSUI 4FQUFNCFS
5SPQJD ET
PG$BQ
SJDPSO
"OUBSD Figure 2.5
UJD$JSD Dates (UTC) of northern-hemisphere seasons relative to Earths
MF
orbit.
4PVUI1PMF
The distance R between the sun and Earth (actu-
ally to the Earth-moon barycenter) as a function of
Figure 2.4
Relationship of declination angle s to tilt of the Earths axis, for
time is
a day near northern-hemisphere summer. 1 e2
R = a (2.4)
1 + e cos( )
where e = 0.0167 is eccentricity, and a = 149.457 Gm

is the semi-major axis length. If the simple approxi-
mation of M is used, then angle errors are less
than 2 and distance errors are less than 0.06%.
Seasonal Effects
The tilt of the Earths axis relative to a line per-
Table 2-1 . Dates and times (UTC) for northern hemi-
pendicular to the ecliptic (i.e., the orbital plane of
sphere equinoxes and solstices. Format: dd hhmm
gives day (dd), hour (hh) and minutes (mm). From the the Earth around the sun) is presently r = 23.44 =
US Naval Observatory. http://aa.usno.navy.mil/data/docs/ 0.409 radians. The direction of tilt of the Earths axis
EarthSeasons.php is not fixed relative to the stars, but wobbles or pre-
cesses like a top with a period of about 25,781 years.
Year Spring Summer Fall Winter
Although this is important over millennia (see the
Equinox Solstice Equinox Solstice
Climate chapter), for shorter time intervals (up to a
(March) (June) (Sept.) (Dec.)
century) the precession is negligible, and you can as-
2015 20 2245 21 1638 23 0820 22 0448 sume a fixed tilt.
2016 20 0430 20 2234 22 1421 21 1044 The solar declination angle s is defined as
2017 20 1028 21 0424 22 2002 21 1628 the angle between the ecliptic and the plane of the
2018 20 1615 21 1007 23 0154 21 2222 Earths equator (Fig. 2.4). Assuming a fixed orienta-
tion (tilt) of the Earths axis as the Earth orbits the
2019 20 2158 21 1554 23 0750 22 0419
sun, the solar declination angle varies smoothly as
2020 20 0349 20 2143 22 1330 21 1002 the year progresses. The north pole points partially
2021 20 0937 21 0332 22 1921 21 1559 toward the sun in summer, and gradually changes
2022 20 1533 21 0914 23 0104 21 2148 to point partially away during winter (Fig. 2.5).
2023 20 2124 21 1458 23 0650 22 0327 Although the Earth is slightly closer to the sun
in winter (near the perihelion) and receives slightly
2024 20 0306 20 2051 22 1244 21 0921
more solar radiation then, the Northern Hemisphere
2025 20 0901 21 0242 22 1819 21 1503 receives substantially less sunlight in winter be-
cause of the tilt of the Earths axis. Thus, summers
are warmer than winters, due to Earth-axis tilt.
The solar declination angle is greatest (+23.44) on
about 20 to 21 June (summer solstice in the North-
ern Hemisphere, when daytime is longest) and is
most negative ( 23.44) on about 21 to 22 December
(winter solstice, when nighttime is longest). The

vernal equinox (or spring equinox, near 20 to 21
March) and autumnal equinox (or fall equinox,
near 22 to 23 September) are the dates when daylight
hours equal nighttime hours (eqi nox literally
translates to equal night), and the solar declina-
tion angle is zero.
Astronomers define a tropical year (= 365.242190
days) as the time from vernal equinox to the next
vernal equinox. Because the tropical year is not an
integral number of days, the Gregorian calen-
dar (the calendar adopted by much of the western
world) must be corrected periodically to prevent the
seasons (dates of summer solstice, etc.) from shifting
into different months.
To make this correction, a leap day (29 Feb) is
added every 4th year (i.e., leap years, are years di-
visible by 4), except that years divisible by 100 dont
have a leap day. However, years divisible by 400 do
have a leap day (for example, year 2000). Because of
all these factors, the dates of the solstices, equinoxes,
perihelion, and aphelion jump around a few days on
the Gregorian calendar (Table 2-1 and Figs. 2.3), but Sample Application
remain in their assigned months. Find the solar declination angle on 5 March.
The solar declination angle for any day of the
year is given by Find the Answer
Assume: Not a leap year.
C (d d ) Given: d = 31 Jan + 28 Feb + 5 Mar = 64.
s r cos r (2.5) Find: s = ?
dy Sketch: /PSUI1PMF

where d is Julian day, and dr is the Julian day of the $BODFS

summer solstice, and dy = 365 (or = 366 on a leap ET 'S
year) is the number of days per year. For years when &RVBUPS
FDMJQUJD
the summer solstice is on 21 June, dr = 172. This
$BQSJDPSO
equation is only approximate, because it assumes
the Earths orbit is circular rather than elliptical. As
before, C = 2 radians = 360 (use radians or degrees 4PVUI1PMF
depending on what your calculator, spreadsheet, or
computer program expects). Use eq. (2.5):
By definition, Earth-tilt angle (r = 23.44) equals s = 23.44 cos[360(64172)/365]
the latitude of the Tropic of Cancer in the North-
= 23.44 cos[106.521] = 6.67
ern Hemisphere (Fig. 2.4). Latitudes are defined to
be positive in the Northern Hemisphere. The Trop-
Check: Units OK. Sketch OK. Physics OK.
ic of Capricorn in the Southern Hemisphere is the Exposition: On the vernal equinox (21 March),
same angle, but with a negative sign. The Arctic the angle should be zero. Before that date, it is winter
Circle is at latitude 90 r = 66.56, and the Ant- and the declination angle should be negative. Namely,
arctic Circle is at latitude 66.56 (i.e., 66.56S). the ecliptic is below the plane of the equator. In spring
During the Northern-Hemisphere summer solstice: and summer, the angle is positive. Because 5 March
the solar declination angle equals r ; the sun at is near the end of winter, we expect the answer to be a
noon is directly overhead (90 elevation angle) at the small negative angle.
Tropic of Cancer; and the sun never sets that day at
the Arctic Circle.
Sample Application Daily Effects

Find the local elevation angle of the sun at 3 PM As the Earth rotates about its axis, the local el-
local time on 5 March in Vancouver, Canada. evation angle of the sun above the local horizon
rises and falls. This angle depends on the latitude
Find the Answer and longitude e of the location:
Assume: Not a leap year. Not daylight savings time.
Also, Vancouver is in the Pacific time zone, (2.6)
where tUTC = t + 8 h during standard time. sin( ) = sin() sin( s )
Given: t = 3 PM = 15 h. Thus, tUTC = 23 h. Ct
= 49.25N, e = 123.1W for Vancouver. cos() cos( s ) cos UTC e
s = 6.665 from previous Sample Application. td
Find: = ?
Use eq. (2.6): sin() =

Time of day tUTC is in UTC, C = 2 radians = 360
= sin(49.25)sin(6.665) as before, and the length of the day is td. For tUTC in
cos(49.25)cos(6.665)cos[360(23h/24h) 123.1] hours, then td = 24 h. Latitudes are positive north
= 0.7576 (0.1161) 0.6528 0.9932 cos(221.9) of the equator, and longitudes are positive west of
= 0.08793 + 0.4826 = 0.3946 the prime meridian. The sin() relationship is used
= arcsin(0.3946) = 23.24 later in this chapter to calculate the daily cycle of
solar energy reaching any point on Earth.
Check: Units OK. Physics OK. [CAUTION: Dont forget to convert angles to radians if
Exposition: The sun is above the local horizon, as required by your spreadsheet or programming language.]
expected for mid afternoon.
Beware of other situations such as night that give
negative elevation angle.
Sample Application() Sample Application (continuation)

Use a spreadsheet to plot elevation angle vs. time at
Vancouver, for 21 Dec, 23 Mar, and 22 Jun. Plot these
three curves on the same graph.
Find the Answer

Given: Same as previous Sample Application, except +VO
4PMBS&MFWBUJPO"OHMF
d = 355, 82, 173.

Find: = ?. Assume not a leap year.

A portion of the tabulated results are shown below, as
.BS
well as the full graph.
()
t (h) 21 Dec 23 Mar 22 Jun
3 0.0 0.0 0.0
4 0.0 0.0 0.0 %FD
5 0.0 0.0 6.6

6 0.0 0.0 15.6

7 0.0 7.8 25.1
-PDBM4UBOEBSE5JNF I
8 0.0 17.3 34.9

9 5.7 25.9 44.5
10 11.5 33.2 53.4
Check: Units OK. Physics OK. Graph OK.
11 15.5 38.4 60.5
Exposition: Summers are pleasant with long days.
12 17.2 40.8 64.1
The peak elevation does not happen precisely at local
13 16.5 39.8 62.6 noon, because Vancouver is not centered within its
time zone.
(continues in next column)
The local azimuth angle of the sun relative to

MPDBMTUBOEBSEUJNF
north is

+VO
sin( s ) sin() cos()
cos( ) = (2.7)
&MFWBUJPO"OHMF

cos() sin()

where = C/4 is the zenith angle. After noon, .BS
the azimuth angle might need to be corrected to be
= C , so that the sun sets in the west instead

of the east. Fig. 2.6 shows an example of the eleva-

tion and azimuth angles for Vancouver (latitude = %FD
49.25N, longitude = 123.1W) during the solstices

and equinoxes.

Sunrise, Sunset & Twilight "[JNVUI"OHMF
Geometric sunrise and sunset occur when the

center of the sun has zero elevation angle. Appar- Figure 2.6
ent sunrise/set are defined as when the top of the Position (solid lines) of the sun for Vancouver, Canada for vari-
sun crosses the horizon, as viewed by an observer on ous seasons. September 21 and March 23 nearly coincide. Iso-
the surface. The sun has a finite radius correspond- chrones are dashed. All times are Pacific standard time.
ing to an angle of 0.267 when viewed from Earth.
Also, refraction (bending) of light through the at-
mosphere allows the top of the sun to be seen even Sample Application()
when it is really 0.567 below the horizon. Thus, ap- Use a spreadsheet to find the Pacific standard time
parent sunrise/set occurs when the center of the sun (PST) for all the events of Table 2-2, for Vancouver,
has an elevation angle of 0.833. Canada, during 22 Dec, 23 Mar, and 22 Jun.
When the apparent top of the sun is slightly be-
low the horizon, the surface of the Earth is not re- Find the Answer
ceiving direct sunlight. However, the surface can Given: Julian dates 355, 82, & 173.
still receive indirect light scattered from air mole- Find: t = ? h (local standard time)
cules higher in the atmosphere that are illuminated Assume: Pacific time zone: tUTC = t + 8 h.
by the sun. The interval during which scattered
Use eq. (2.8a) and Table 2-2:
light is present at the surface is called twilight.
Because twilight gradually fades as the sun moves 22Dec 23Mar 22Jun
lower below the horizon, there is no precise defini- Morning: PST (h)
tion of the start of sunrise twilight or the end of sun- geometric sunrise 8.22 6.20 4.19
set twilight. apparent sunrise 8.11 6.11 4.09
Arbitrary definitions have been adopted by dif-
civil twilight starts 7.49 5.58 3.36
ferent organizations to define twilight. Civil twi-
light occurs while the sun center is no lower than military twilight starts 6.80 4.96 2.32
6, and is based on the ability of humans to see astron. twilight starts 6.16 4.31 n/a
objects on the ground. Military twilight occurs Evening:
while the sun is no lower than 12. Astronomi-
geometric sunset 16.19 18.21 20.22
cal twilight ends when the skylight becomes suf-
ficiently dark to view certain stars, at solar elevation apparent sunset 16.30 18.30 20.33
angle 18. civil twilight ends 16.93 18.83 21.05
Table 2-2 summarizes the solar elevation angle military twilight ends 17.61 19.45 22.09
definitions used for sunrise, sunset and twilight. astron. twilight ends 18.26 20.10 n/a
All of these angles are at or below the horizon.
Check: Units OK. Physics OK.
Exposition: During the summer solstice (22 June),
the sun never gets below 18. Hence, it is astronomi-
cal twilight all night in Vancouver in mid summer.
During June, Vancouver is on daylight time, so the
actual local time would be one hour later.
Table 2-2. Elevation angles for diurnal events. Approximate (sundial) time-of-day correspond-
ing to these events can be found by rearranging eq.
Event () (radians) (2.6): (2.8a)
Sunrise & Sunset:
Geometric 0 0 td sin sin s sin
tUTC = e arccos
Apparent 0.833 0.01454 C cos cos s
Twilight: where the appropriate elevation angle is used from
Civil 6 0.10472 Table 2-2. Where the sign appears, use + for sun-
Military 12 0.20944 rise and for sunset. If any of the answers are nega-
tive, add 24 h to the result.
Astronomical 18 0.31416
To correct the time for the tilted, elliptical orbit of
the earth, use the approximate Equation of Time:
INFO Astronomical Values for Time
ta a sin( M ) + bsin(2 M + c) (2.8b)
The constants a = 2e/E , b = [tan2(/2)]/E , and
c = 2 in the Equation of Time are based on the fol-
lowing astronomical values: E = 2/24h = 0.0043633 where a = 7.659 minutes, b = 9.863 minutes, c = 3.588
radians/minute is Earths rotation rate about its axis, radians = 205.58, and where the mean anomaly M
e = 0.01671 is the eccentricity of Earths orbit around from eq. (2.2) varies with day of the year. This time
the sun, = 0.40909 radians = 23.439 is the obliquity correction is plotted in the Sample Application. The
(tilt of Earths axis), = 4.9358 radians = 282.8 is the corrected (mechanical-clock) time teUTC is:
angle (in the direction of Earths orbit around the sun)
between a line from the sun to the vernal (Spring) teUTC = tUTC ta (2.8c)
equinox and a line drawn from the sun to the moving
perihelion (see Fig. 2.5, and Fig. 21.10 in Chapter 21). which corrects sundial time to mechanical-clock
time. Dont forget to convert the answer from UTC
to your local time zone.
Sample Application ()
Plot time correction vs. day of the year.
Find the Answer:

Given: d = 0 to 365 Flux
Find: ta = ? minutes
A flux density, F, called a flux in this book, is
Use a spreadsheet. For example, for d = 45 (which is 14
the transfer of a quantity per unit area per unit time.
Feb), first use eq. (2.2) find the mean anomaly:
M = 2(45-4)/365.25 = 0.705 rad = 40.4
The area is taken perpendicular (normal) to the di-
Then use the Equation of Time (2.8b) in rection of flux movement. Examples with metric (SI)
ta = (7.659 min)sin(0.705 rad) units are mass flux (kg m2s1) and heat flux, ( Jm
2s1). Using the definition of a watt (1 W = 1 Js1),
+ (9.836 min)sin(2 0.705 rad + 3.588 rad)
ta = 4.96 9.46 = 14.4 minutes the heat flux can also be given in units of (Wm2). A
+'.".++"40/% flux is a measure of the amount of inflow or outflow
such as through the side of a fixed volume, and thus
is frequently used in Eulerian frameworks (Fig. 2.7).
/PW
Because flow is associated with a direction, so is

UB
flux associated with a direction. You must account
NJOVUFT
for fluxes Fx, Fy, and Fz in the x, y, and z directions,

+VM
.BZ

respectively. A flux in the positive x-direction (east-
'FC
m ward) is written with a positive value of Fx, while

m
a flux towards the opposite direction (westward) is
mNJOVUFT negative.
m The total amount of heat or mass flowing through
a plane of area A during time interval t is given
EBZPGZFBS
by:
Check: Units OK. Magnitude OK.
Exposition: Because of the negative sign in eq. (2.8c), Amount = F A t (2.9)
you need to ADD 14.4 minutes to the results of eq.
(2.8a) to get the correct sunrise and sunset times. For heat, Amount QH by definition.
Fluxes are sometimes written in kinematic

form, F, by dividing by the air density, air :
'PS'
F
F= (2.10)
air
Kinematic mass flux equals the wind speed, M. Ki-

nematic fluxes can also be in the 3 Cartesian direc-
tions: Fx, Fy, and Fz.
Figure 2.7
[CAUTION: Do not confuse the usage of the symbol
Flux F through an area A into one side of a volume.
M. Here it means wind speed. Earlier it meant mean
anomaly. Throughout this book, many symbols will be
re-used to represent different things, due to the limited Sample Application
number of symbols. Even if a symbol is not defined, you The mass flux of air is 1 kgm2s1 through a door
can usually determine its meaning from context.] opening that is 1 m wide by 2.5 m tall. What amount of
Heat fluxes FH can be put into kinematic form by mass of air passes through the door each minute, and
dividing by both air density air and the specific heat what is the value of kinematic mass flux?
for air Cp, which yields a quantity having the same
units as temperature times wind speed (Kms1). Find the Answer
Given: area A = (1 m) (2.5 m) = 2.5 m2
F = 1 kgm2s1
FH Find: (a) Amount = ? kg, and (b) F = ? ms1
FH = for heat only (2.11)
air Cp Sketch: (see Fig. 2.7)
(a) Use eq. (2.9):

For dry air (subscript d) at sea level: Amount = (1 kgm2s1)(2.5 m2)
(1 min)(60 s min1) = 150 kg.
air Cpd = 1231 (Wm2) / (Kms1)
(b) Assume: = 1.225 kgm3 at sea-level
= 12.31 mbK1 Use eq. (2.10):
F = (1 kgm2s1) / (1.225 kgm3)
= 1.231 kPaK1 . = 0.82 ms 1.
The reason for sometimes expressing fluxes in Check: Units OK. Sketch OK. Physics OK.
kinematic form is that the result is given in terms Exposition: The kinematic flux is equivalent to a very
of easily measured quantities. For example, while light wind speed of less than 1 m s1 blowing through
the doorway, yet it transports quite a large amount of
most people do not have Watt meters to measure
mass each minute.
the normal dynamic heat flux, they do have ther-
mometers and anemometers. The resulting temper-
ature times wind speed has units of a kinematic heat
flux (Kms1). Similarly, for mass flux it is easier to Sample Application
measure wind speed than kilograms of air per area If the heat flux from the sun that reaches the Earths
per time. surface is 600 W m2, find the flux in kinematic units.
Heat fluxes can be caused by a variety of pro-
cesses. Radiative fluxes are radiant energy (elec- Find the Answer:
Given: FH = 600 W m2
tromagnetic waves or photons) per unit area per
Find: FH = ? Km s1
unit time. This flux can travel through a vacuum.
Assume: sea level.
Advective flux is caused by wind blowing through
an area, and carrying with it warmer or colder Use eq. (2.11)
temperatures. For example a warm wind blowing FH = (600 W m2) / [1231 (Wm2) / (Kms1) ]
toward the east causes a positive heat-flux compo- = 0.487 Km s1
nent FHx in the x-direction. A cold wind blowing
toward the west also gives positive FHx. Turbulent Check: Units OK. Physics OK.
fluxes are caused by eddy motions in the air, while Exposition: This amount of radiative heat flux is
conductive fluxes are caused by molecules bounc- equivalent to an advective flux of a 1 m s1 wind blow-
ing into each other. ing air with temperature excess of about 0.5C.
Sample Application
Red light has a wavelength of 0.7 m. Find its Radiation principles
frequency, circular frequency, and wavenumber in a
vacuum.
Propagation
Find the Answer Radiation can be modeled as electromagnetic
Given: co = 299,792,458 m s1, = 0.7 m waves or as photons. Radiation propagates through
Find: = ? Hz, = ? s1 , = ? m1 . a vacuum at a constant speed: co = 299,792,458 ms1.
For practical purposes, you can approximate this
Use eq. (2.12), solving for : speed of light as co 3x108 ms1. Light travels slight-
= co/ = (3x108 m s1) / (0.7x10 6 m cycle 1) ly slower through air, at roughly c = 299,710,000 ms1
= 4.28x1014 Hz. at standard sea-level pressure and temperature, but
the speed varies slightly with thermodynamic state
= 2 = 2(3.14159)(4.28x1014 Hz)
of the air (see the Optics chapter).
= 2.69x1015 s1.
Using the wave model of radiation, the wave-
= 1/ = 1 / (0.7x10 6 m cycle 1) length (mcycle1) is related to the frequency,
= 1.43x10 6 m1. (Hz = cycless1) by:
Check: Units OK. Physics reasonable. = co (2.12)

Exposition: Wavelength, wavenumber, frequency,
and circular frequency are all equivalent ways to ex- Wavelength units are sometimes abbreviated as (m).
press the color of radiation. Because the wavelengths of light are so short, they
are often expressed in units of micrometers (m).
Wavenumber is the number of waves per me-
ter: (cyclesm1) = 1/. Its units are sometimes
abbreviated as (m1). Circular frequency or an-
gular frequency is (radians s1) = 2. Its units
are sometimes abbreviated as (s1).
Emission
Objects warmer than absolute zero can emit ra-
diation. An object that emits the maximum possible
radiation for its temperature is called a blackbody.
Sample Application Plancks law gives the amount of blackbody mono-
Find the blackbody monochromatic radiant exi-
chromatic (single wavelength or color) radiant flux
tance of green light of wavelength 0.53 m from an
leaving an area, called emittance or radiant exi-
object of temperature 3000 K.
tance, E*:
Find the Answer c1
E * = (2.13)
Given: = 0.53 m, T = 3000 K [ exp(c2 / ( T )) 1]
5
Find: E* = ? Wm2 m1 .
Use eq. (2.13): where T is absolute temperature, and the asterisk in-
dicates blackbody. The two constants are:
c1
E * = 5 c1 = 3.74 x 108 Wm2 m4 , and
[ exp(c2 /( T )) 1] c2 = 1.44 x 104 mK.
(3.74x108 Wm -2 m 4 )/(0.53m)5 Eq. (2.13) and constant c1 already include all di-
E * = rections of exiting radiation from the area. E* has
exp (1.44x10 4 Km)/(0.53m 3000K) 1
units of Wm-2 m1 ; namely, flux per unit wave-
= 1.04 x 10 6 Wm2 m1 length. For radiation approaching an area rather
than leaving it, the radiant flux is called irradi-
Check: Units OK. Physics reasonable. ance.
Exposition: Because 3000 K is cooler than the sun, Actual objects can emit less than the theoreti-
about 50 times less green light is emitted. The answer cal blackbody value: E = eE* , where 0 e 1 is
is the watts emitted from each square meter of surface emissivity, a measure of emission efficiency.
per m wavelength increment.

Y 4VO
&BSUI
5,
5,
&M 8pNmpNm
&M 8pNmpNm
Y
Y

Y

8BWFMFOHUI N
8BWFMFOHUI N
Figure 2.8 Figure 2.9

Planck radiant exitance, E*, from a blackbody approximately Planck radiant exitance, E*, from a blackbody approximately
the same temperature as the sun. the same temperature as the Earth.
The Planck curve (eq. 2.13) for emission from a Table 2-3. Ranges of wavelengths of visible colors.
blackbody the same temperature as the sun (T = Approximate center wavelength is in boldface. For
5780 K) is plotted in Fig. 2.8. Peak emissions from more info, see Chapter 22 on Atmospheric Optics.
the sun are in the visible range of wavelengths (0.38
Color (m)
0.74 m, see Table 2-3). Radiation from the sun is
called solar radiation or short-wave radiation.
Red 0.625 - 0.650 - 0.740
The Planck curve for emission from a blackbody Orange 0.590 - 0.600 - 0.625
that is approximately the same temperature as the Yellow 0.565 - 0.577 - 0.590
whole Earth-atmosphere system (T 255 K) is plot- Green 0.500 - 0.510 - 0.565
ted in Fig. 2.9. Peak emissions from this idealized Cyan 0.485 - 0.490 - 0.500
average Earth system are in the infrared range 8 to Blue 0.460 - 0.475 - 0.485
18 m. This radiation is called terrestrial radia- Indigo 0.425 - 0.445 - 0.460
tion, long-wave radiation, or infrared (IR) radia- Violet 0.380 - 0.400 - 0.425
tion.
The wavelength of the peak emission is given by
Wiens law: Sample Application
What is the total radiant exitance (radiative flux)
a emitted from a blackbody Earth at T = 255 K, and what
max = (2.14)
T is the wavelength of peak emission?
where a = 2897 mK.
Find the Answer
The total amount of emission (= area under
Given: T = 255 K &BSUI
Planck curve = total emittance) is given by the Ste- Find: E* = ? Wm2 , max = ? m
fan-Boltzmann law: Sketch:
Use eq. (2.15):
E*= (5.67x10 8 Wm2K4)(255 K)4 = 240 Wm2.
E* = SB T 4 (2.15) Use eq. (2.14):
max = (2897 mK) / (255 K) = 11.36 m
where SB = 5.67x10 8 Wm2K4 is the Stefan-
Boltzmann constant, and E* has units of Wm2. Check: Units OK. agrees with peak in Fig. 2.9.
More details about radiation emission are in the Sat- Exposition: You could create the same flux by plac-
ellites & Radar chapter in the sections on weather ing a perfectly-efficient 240 W light bulb in front of a
satellites. parabolic mirror that reflects the light and IR radiation
into a beam that is 1.13 m in diameter.
For comparison, the surface area of the Earth is
about 5.1x1014 m2, which when multiplied by the flux
gives the total emission of 1.22x1017 W. Thus, the Earth
acts like a large-wattage IR light bulb.
HIGHER MATH Incremental Changes A SCIENTIFIC PERSPECTIVE Scientific

Laws The Myth
What happens to the total radiative exitance given There are no scientific laws. Some theories or
by eq. (2.15) if temperature increases by 1C? Such a models have succeeded for every case tested so far,
question is important for climate change. yet they may fail for other situations. Newtons Laws
of Motion were accepted as laws for centuries, until
Find the Answer using calculus: they were found to fail in quantum mechanical, rela-
Calculus allows a simple, elegant way to find the tivistic, and galaxy-size situations. It is better to use
solution. By taking the derivative of both the left and the word relationship instead of law. In this text-
right sides of eq. (2.15), one finds that: book, the word law is used for sake of tradition.
Einstein said No amount of experimentation can
dE* = 4 SB T 3dT ever prove me right; a single experiment can prove me
wrong.
assuming that SB is constant. It can be written as Because a single experiment can prove a relation-
ship wrong, it behooves us as scientists to test theories
E* = 4 SB T 3 T and equations not only for reasonable values of vari-
ables, but also in the limit of extreme values, such as
for small T. when the variables approach zero or infinity. These
Thus, a fixed increase in temperature of T = 1C are often the most stringent tests of a relationship.
causes a much larger radiative exitance increase at
high temperatures than at cold, because of the T3 fac- Example
tor on the right side of the equation.
Query: The following is an approximation to
Find the Answer using algebra: Plancks law.
This particular problem could also have been
solved using algebra, but with a more tedious and less E* = c1 5 exp[ c2 / ( T) ] (a)
elegant solution. First let
Why is it not a perfect substitute?
E2 = SB T24 and E1 = SB T14
Find the Answer: If you compare the numeri-
Next, take the difference between these two eqs:
cal answers from eqs. (2.13) and (a), you find that they
agree very closely over the range of temperatures of
E = E2 E1 = T24 T14
the sun and the Earth, and over a wide range of wave-
lengths. But...
Recall from algebra that (a2 b2) = (a b)(a + b)
a) What happens as temperature approaches
( )( )

E = SB T2 2 T12 T2 2 + T12
absolute zero? For eq. (2.13), T is in the denominator
of the argument of an exponential, which itself is in
E = SB (T2 T1 ) (T2 + T1 ) ( T2 + T12
2
) the denominator of eq. (2.13). If T = 0, then 1/T = .
Exp() = , and 1/ = 0. Thus, E* = 0, as it should.
Since (T2 T1) / T1 << 1, then (T2 T1 ) = T, Namely, no radiation is emitted at absolute zero (ac-
but T2 + T1 2T, and T22 + T12 2T2. This gives: cording to classical theory).
For eq. (a), if T = 0, then 1/T = , and exp() = 0.
E = SB T 2T 2T 2 So it also agrees in the limit of absolute zero. Thus,
both equations give the expected answer.
or
b) What happens as temperature approaches
E SB 4T 3 T
infinity? For eq. (2.13), if T = , then 1/T = 0 , and
which is identical to the answer from calculus. exp(0) = 1. Then 1 1 = 0 in the denominator, and 1/0
We were lucky this time, but it is not always pos- = . Thus, E* = , as it should. Namely, infinite
sible to use algebra where calculus is needed. radiation is emitted at infinite temperature.
For eq. (a), if T = , then 1/T = 0 , and exp(0) = 1.
Thus, E* = c1 5 , which is not infinity. Hence, this
approaches the wrong answer in the T limit.
Conclusion: Eq. (a) is not a perfect relationship,

because it fails for this one case. But otherwise it is
a good approximation over a wide range of normal
temperatures.
Distribution
Radiation emitted from a spherical source de- 4PMBSSBEJBUJPOSFBDIJOH
creases with the square of the distance from the cen- &BSUITPSCJU
&M 8pNmpNm
ter of the sphere:

2
R1 &BSUIT
* * (2.16)
E2 = E1 SBEJBUJPO
R2
where R is the radius from the center of the sphere,

and the subscripts denote two different distances
from the center. This is called the inverse square

law.
The reasoning behind eq. (2.16) is that as ra-
8BWFMFOHUI N
diation from a small sphere spreads out radially, it

passes through ever-larger conceptual spheres. If
no energy is lost during propagation, then the to- Figure 2.10
tal energy passing across the surface of each sphere Blackbody radiance E* reaching top of Earths atmosphere from
must be conserved. Because the surface areas of the the sun and radiance of terrestrial radiation leaving the top of
spheres increase with the square of the radius, this the atmosphere, plotted on a log-log graph.
implies that the energy flux density must decrease
at the same rate; i.e., inversely to the square of the
radius.
From eq. (2.16) we expect that the radiative flux
reaching the Earths orbit is greatly reduced from
that at the surface of the sun. The solar emissions
of Fig. 2.8 must be reduced by a factor of 2.167x10 5, Sample Application
based on the square of the ratio of solar radius to Estimate the value of the solar irradiance reaching
the orbit of the Earth, given a sun surface temperature
Earth-orbital radius from eq. (2.16). This result is
(5780 K), sun radius (6.96x105 km), and orbital radius
compared to the emission from Earth in Fig. 2.10.
(1.495x108 km) of the Earth from the sun.
The area under the solar-radiation curve in Fig.
2.10 is the total (all wavelengths) solar irradiance Find the Answer
(TSI), So, reaching the Earths orbit. We call it an ir- Given: Tsun = 5780 K
radiance here, instead of an emittance, because rela- Rsun = 6.96x105 km = solar radius
tive to the Earth it is an incoming radiant flux. This REarth = 1.495x108 km = Earth orbit radius
quantity was formerly called the solar constant Find: So = ? Wm2
but we now know that it varies slightly. The aver- Sketch:
age value of solar irradiance measured at the Earths &BSUI
orbit by satellites is about
So = 1366 (7) Wm2 . (2.17) 4VO
In kinematic units (based on sea-level density), the

solar irradiance is roughly So = 1.11 Km s1.
About 4 Wm2 of the 7 Wm2 error bars are
due to radiometer calibration errors between the First, use eq. (2.15):
various satellites. Also, during the 11-year sunspot E1* = (5.67x10 8 Wm2K4)(5780 K)4
cycle the solar irradiance normally varies by about = 6.328x107 Wm2.
1.4 Wm2, with peak values during sunspot maxima. Next, use eq. (2.16), with R1 = Rsun & R 2 = REarth
One example of a longer term variation in solar ac- So = E2*=(6.328x107 Wm2)
tivity is the Maunder Minimum in the late 1600s, (6.96x105 km/1.495x108 km)2 = 1372 Wm2.
during which the solar irradiance was perhaps 2.7
to 3.7 Wm2 less than values during modern solar Check: Units OK, Sketch OK. Physics OK.
minima. Such irradiance changes could cause sub- Exposition: Answer is nearly equal to that measured
by satellites, as given in eq. (2.17). The error is due to a
tle changes (0.3 to 0.4C) in global climate. See the
poor estimate of effective sun-surface temperature.
Climate chapter for more info on solar variability.
According to the inverse-square law, variations

Sample Application()
of distance between Earth and sun cause changes of
Using the results from an earlier Sample Applica-
tion that calculated the true anomaly and sun-Earth the solar radiative forcing, S, that reaches the top
distance for several days during the year, find the solar of the atmosphere:
radiative forcing for those days.
2
R
Find the Answer S = So (2.18)
Given: R values from previously Sample Application R
Find: S = ? Wm2
where So = 1366 Wm2 is the average total solar ir-
Sketch: (same as Fig 2.2) radiance measured at an average distance R = 149.6
Use eq. (2.18). Gm between the sun and Earth, and R is the actual
Date d R() S distance between Earth and the sun as given by eq.
(Gm) (W/m2) (2.4). Remember that the solar irradiance and the
4 Jan 4 146.96 1418 solar radiative forcing are the fluxes across a surface
18 Jan 18 147.04 1416 that is perpendicular to the solar beam, measured
1 Feb 32 147.25 1412 above the Earths atmosphere.
15 Feb 46 147.60 1405 Let irradiance E be any radiative flux crossing
1 Mar 60 148.06 1397 a unit area that is perpendicular to the path of the
15 Mar 74 148.60 1386 radiation. If this radiation strikes a surface that is
29 Mar 88 149.18 1376 not perpendicular to the radiation, then the radia-
12 Apr 102 149.78 1365 tion per unit surface area is reduced according to the
26 Apr 116 150.36 1354 sine law. The resulting flux, Frad, at this surface is:
21 Jun 172 151.88 1327
23 Sep 266 150.01 1361 (2.19)
Frad = E sin( )
22 Dec 356 147.03 1416
Check: Units OK. Physics OK. where is the elevation angle (the angle of the
Exposition: During N. Hemisphere winter, solar ra- sun above the surface). In kinematic form, this is
diative forcing is up to 50 Wm2 larger than average.
E
Frad = sin( )
Cp (2.20)
Sample Application where Cp is given under eq. (2.11).

During the equinox at noon at latitude =60, the
solar elevation angle is = 90 60 = 30. If the at-
mosphere is perfectly transparent, then how much ra- Average Daily Insolation
diative flux is absorbed into a perfectly black asphalt The acronym insolation means incoming
parking lot? solar radiation at the top of the atmosphere. The
average daily insolation E takes into account both
Find the Answer the solar elevation angle (which varies with season
Given: = 30 = elevation angle and time of day) and the duration of daylight. For
E = So = 1366 Wm2. solar irradiance example, there is more total insolation at the poles
Find: Frad = ? Wm2 in summer than at the equator, because the low sun
Sketch:
TVOMJHIU angle near the poles is more than compensated by

the long periods of daylight.

S a 2 (2.21)
&BSUI E = o ho sin()sin( s ) +
R
Use eq. (2.19): cos()) cos( s )sin( ho )
Frad = (1366 Wm2)sin(30) = 683 Wm2.
Check: Units OK. Sketch OK. Physics OK. where So =1366 W m2 is the solar irradiance, a =
Exposition: Because the solar radiation is striking 149.457 Gm is Earths semi-major axis length, R is
the parking lot at an angle, the radiative flux into the the actual distance for any day of the year, from eq.
parking lot is half of the solar irradiance.
(2.4). In eq. (2.21), ho is the sunset and sunrise hour +BO 'FC .BS "QS .BZ +VO +VM "VH 4FQ 0DU /PW %FD
angle in radians.
The hour angle ho at sunrise and sunset can be

found using the following steps:

= tan()tan( s ) (2.22a)

= min[1, (max( 1, )] (2.22b)
-BUJUVEF
8N
ho = arccos() (2.22c)

Eq. (2.22b) truncates the argument of the arccos to
be between 1 and 1, in order to account for high m

latitudes where there are certain days when the sun

never sets, and other days when it never rises.
m
[CAUTION. When finding the arccos, your answer
might be in degrees or radians, depending on your calcu-
lator, spreadsheet, or computer program. Determine the m
units by experimenting first with the arccos(0), which will
either give 90 or /2 radians. If necessary, convert the 3FMBUJWF+VMJBO%BZ
hour angle to units of radians, the result of which is ho.]
Fig. 2.11 shows the average incoming solar radia- Figure 2.11
tion vs. latitude and day of the year, found using Average daily insolation E (W m2) over the globe.
eq. (2.21). For any one hemisphere, E has greater
difference between equator and pole during winter
than during summer. This causes stronger winds
and more active extratropical cyclones in the winter Sample Application
hemisphere than in the summer hemisphere. Find the average daily insolation over Vancouver
during the summer solstice.
Absorption, Reflection & Transmission Find the Answer:

The emissivity, e, is the fraction of blackbody Given: d = dr = 173 at the solstice,
radiation that is actually emitted (see Table 2-4) at = 49.25N, e = 123.1W for Vancouver.
any wavelength . The absorptivity, a, is the frac- Find: E = ? W m2
tion of radiation striking a surface of an object that
is absorbed (i.e., stays in the object as a different Use eq. (2.5): s = r = 23.45
form of energy). Kirchhoffs law states that the Use eq. (2.2): M = 167.55, and assume M.
absorptivity equals the emissivity of a substance Use eq. (2.4): R = 151.892 Gm.
Use eq. (2.22):
at each wavelength, . Thus,
ho = arccos[tan(49.25)tan(23.45)] = 120.23
ho = ho2/360 = 2.098 radians
a = e (2.23) Use eq. (2.21):
Some substances such as dark glass are semi- (1368 Wm -2 ) 149Gm 2

E=
transparent (i.e., some radiation passes through). 151.892 Gm
A fraction of the incoming (incident) radiation might [2.098sin( 49.25)sin(23.45) +
also be reflected (bounced back), and another por- cos( 49.25)cos(23.45)sin(120.23)]
tion might be absorbed. Thus, you can define the
efficiencies of reflection, absorption, and transmis- E = (1327 W m2)[2.098(0.3016)+0.5174]
sion as:
= 486 W m2
E reflected
r = = reflectivity (2.24)
E incident Check: Units OK. Physics OK. Agrees with Fig. 2.11.
Exposition: At the equator on this same day, the av-
erage daily insolation is less than 400 W m2 .
E
a = absorbed = absorptivity (2.25)
E incident
Table 2-4. Typical infrared emissivities. E transmitted

t = = transmissivity (2.26)
Surface e Surface e E incident
alfalfa 0.95 iron, galvan. 0.13-0.28
aluminum 0.01-0.05 leaf 0.8 m 0.05-0.53 Values of e, a, r, and t are between 0 and 1.
asphalt 0.95 leaf 1 m 0.05-0.6 The sum of the last three fractions must total 1,
bricks, red 0.92 leaf 2.4 m 0.7-0.97 as 100% of the radiation at any wavelength must be
cloud, cirrus 0.3 leaf 10 m 0.97-0.98 accounted for:
cloud, alto 0.9 lumber, oak 0.9
cloud, low 1.0 paper 0.89-0.95 1 = a + r + t (2.27)
concrete 0.71-0.9 plaster, white 0.91 or
desert 0.84-0.91 sand, wet 0.98 (2.28)
forest, conif. 0.97 sandstone 0.98 E incoming = E absorbed +E reflected +E transmitted
forest, decid. 0.95 shrubs 0.9
glass 0.87-0.94 silver 0.02 For opaque (t = 0) substances such as the Earths
grass 0.9-0.95 snow, fresh 0.99 surface, you find: a = 1 r.
grass lawn 0.97 snow, old 0.82 The reflectivity, absorptivity, and transmissivity
gravel 0.92 soils 0.9-0.98 usually vary with wavelength. For example, clean
human skin 0.95 soil, peat 0.97-0.98 snow reflects about 90% of incoming solar radiation,
ice 0.96 urban 0.85-0.95 but reflects almost 0% of IR radiation. Thus, snow
is white in visible light, and black in IR. Such
behavior is crucial to surface temperature forecasts.
Instead of considering a single wavelength, it is
Sample Application
also possible to examine the net effect over a range
If 500 W m2 of visible light strikes a translucent
object that allows 100 W m2 to shine through and 150 of wavelengths. The ratio of total reflected to total
W m2 to bounce off, find the transmissivity, reflectiv- incoming solar radiation (i.e., averaged over all solar
ity, absorptivity, and emissivity. wavelengths) is called the albedo, A :
Ereflected (2.29)
Find the Answer A=
Given: E incoming = 500 W m2 , Eincoming
E transmitted = 100 W m2 , E reflected = 150 W m2
Find: a = ? , e = ? , r = ? , and t = ? The average global albedo for solar radiation reflect-
ed from Earth is A = 30% (see the Climate chapter).
Use eq. (2.26): t = (100 W m2) / (500 W m2) = 0.2 The actual global albedo at any instant varies with
Use eq. (2.24): r = (150 W m2) / (500 W m2) = 0.3
ice cover, snow cover, cloud cover, soil moisture, to-
Use eq. (2.27): a = 1 0.2 0.3 = 0.5
pography, and vegetation (Table 2-5). The Moons
Use eq. (2.23): e = a = 0.5
albedo is only 7%.
Check: Units dimensionless. Physics reasonable. The surface of the Earth (land and sea) is a very
Exposition: By definition, translucent means part- strong absorber and emitter of radiation.
ly transparent, and partly absorbing.
Table 2-5. Typical albedos (%) for sunlight. Table 2-5 (continuation). Typical albedos (%).
Surface A (%) Surface A (%) Surface A (%) Surface A (%)
alfalfa 23-32 forest, decid. 10-25 rice paddy 12 soil, red 17
buildings 9 granite 12-18 road, asphalt 5-15 soil, sandy 20-25
clay, wet 16 grass, green 26 road, dirt 18-35 sorghum 20
clay, dry 23 gypsum 55 rye winter 18-23 steppe 20
cloud, thick 70-95 ice, gray 60 sand dune 20-45 stones 20-30
cloud, thin 20-65 lava 10 savanna 15 sugar cane 15
concrete 15-37 lime 45 snow, fresh 75-95 tobacco 19
corn 18 loam, wet 16 snow, old 35-70 tundra 15-20
cotton 20-22 loam, dry 23 soil, dark wet 6-8 urban, mean 15
field, fallow 5-12 meadow, green 10-20 soil, light dry 16-18 water, deep 5-20
forest, conif. 5-15 potatoes 19 soil, peat 5-15 wheat 10-23
Within the air, however, the process is a bit more

complicated. One approach is to treat the whole at- &JODJEFOU
mospheric thickness as a single object. Namely, you
can compare the radiation at the top versus bottom of
the atmosphere to examine the total emissivity, ab-
sorptivity, and reflectivity of the whole atmosphere. TQBUIMFOHUI
Over some wavelengths called windows there is
little absorption, allowing the radiation to shine
through. In other wavelength ranges there is partial
&USBOTNJUUFE QBSUJDMFPS
or total absorption. Thus, the atmosphere acts as a NPMFDVMF
filter. Atmospheric windows and transmissivity are
discussed in detail in the Satellites & Radar and Cli-
mate chapters.
Figure 2.12
Reduction of radiation across an air path due to absorption and
Beers Law scattering by particles such as air-pollution aerosols or cloud
Sometimes you must examine radiative extinc- droplets, illustrating Beers law.
tion (reduction of radiative flux) across a short path
length s within the atmosphere (Fig. 2.12). Let n be
the number density of radiatively important par- Sample Application
ticles in the air (particles m3), and b be the extinc- Suppose that soot from a burning automobile tire
tion cross section of each particle (m2 particle1), has number density n = 107 m3, and an extinction
where this latter quantity gives the area of the shad- cross section of b = 10 9 m2. (a) Find the radiative flux
that was not attenuated (i.e., reduction from the solar
ow cast by each particle.
constant) across the 20 m diameter smoke plume. As-
Extinction can be caused by absorption and scat- sume the incident flux equals the solar constant. (b)
tering of radiation. If the change in radiation is due Find the optical depth. (c) If this soot fills the air, what
only to absorption, then the absorptivity across this is the visual range?
layer is
E Etransmitted Find the Answer
a = incident
E (2.30) Given: n = 107 m3 , b = 10 9 m2 , s = 20 m
incident
Eincident = S = 1366 Wm2 solar constant
Beers law gives the relationship between inci- Find: Etran =? Wm2, =? (dimensionless), V =? km
dent radiative flux, Eincident, and transmitted radia- Sketch:
TVOMJHIU
tive flux, Etransmitted, as
Etransmitted = Eincident e n b s (2.31a)
Beers law can be written using an extinction coef-

ficient, k:
(2.31b)
Etransmitted = Eincident e k s
where is the density of air, and k has units of Use eq. (2.32): = (107 m3)(10 9 m2) = 0.01 m1
m2 gair1 . The total extinction across the whole path
can be quantified by a dimensionless optical thick- (a) Use eq. (2.31a)
ness (or optical depth in the vertical), , allowing Etran = 1366 Wm2 exp[ (0.01 m1)(20 m)]
Beers law to be rewritten as: = 1366 Wm2 exp[0.2] = 1118 Wm2
E =E e (2.31c)
transmitted incident (b) Rearrange eq. (2.31c):
= ln(Ein/Etransmitted) = ln(1366/1118) = 0.20
To simplify these equations, sometimes a volume
extinction coefficient is defined by (c) Rearrange eq. (2.31a): s = [ln(Eincident/Etrans)]/ .
V = s = [ln(1/0.02)]/0.01 m1 = 0.391 km
= nb = k (2.32)
Check: Units OK. Sketch OK. Physics OK.
Visual range (V, one definition of visibility) is Exposition: Not much attenuation through this small
the distance where the intensity of transmitted light smoke plume; namely, 1366 1118 = 248 Wm2 was ab-
has decreased to 2% of the incident light. It estimates sorbed by the smoke. However, if this smoke fills the
the max distance s (km) you can see through air. air, then visibility is very poor.

*m Surface Radiation Budget
'MVY 8N
' ,m Define F* as the net radiative flux (positive

upward) perpendicular to the Earths surface. This
net flux has contributions (Fig. 2.13) from down-
welling solar radiation K , reflected upwelling
m *o solar K , downwelling longwave (IR) radiation
emitted from the atmosphere I , and upwelling
longwave emitted from the Earth I:
m ,o
F* = K + K + I + I (2.33)
where K and I are negative because they are

-PDBM5JNF I
downward.
Figure 2.13
Typical diurnal variation of radiative fluxes at the surface. Solar
Fluxes are positive upward. Recall that the solar irradiance (i.e., the solar con-
stant) is So 1366 7 Wm2 (equivalent to 1.11 Km
s1 in kinematic form after dividing by Cp ) at the
top of the atmosphere. Some of this radiation is at-
DMPVET
tenuated between the top of the atmosphere and the
4 ,o surface (Fig. 2.14). Also, the sine law (eq. 2.19) must
QBUI be used to find the component of downwelling solar
TVOMJHIU : flux K that is perpendicular to the surface. The re-
sult for daytime is
IFSF
K = So Tr sin( ) (2.34)
&BSUI
BUNPTQ
where Tr is a net sky transmissivity. A negative

sign is incorporated into eq. (2.34) because K is a
Figure 2.14 downward flux. Eq. (2.6) can be used to find sin().
Fate of sunlight en route to the Earths surface. At night, the downwelling solar flux is zero.
Net transmissivity depends on path length
through the atmosphere, atmospheric absorption
Sample Application characteristics, and cloudiness. One approximation
Downwelling sunlight shines through an atmo-
for the net transmissivity of solar radiation is
sphere with 0.8 net sky transmissivity, and hits a
ground surface of albedo 0.5, at a time when sin() =
0.3. The surface emits 400 W m2 IR upward into the Tr = (0.6 + 0.2 sin )(1 0.4 H )(1 0.7 M )(1 0.4 L )
atmosphere, and absorbs 350 W m2 IR coming down
from the atmosphere. Find the net radiative flux. (2.35)
Find the Answer where cloud-cover fractions for high, middle, and
Given: Tr = 0.8, So = 1366 W m2 , A = 0.5. low clouds are H, M, and L, respectively. These
I = 400 W m2 , I = 350 W m2 cloud fractions vary between 0 and 1, and the
Find: F* = ? W m2 transmissivity also varies between 0 and 1.
Of the sunlight reaching the surface, a portion
Use eq. (2.34): might be reflected:
K =(1366 W m2)(0.8) (0.3) = 381 W m2
Use eq. (2.36): K= (0.5)(381 W m2)= 164 W m2
K = A K (2.36)
Use eq. (2.33): F* =(381)+(164)+(350)+(400) W m2
F* = 167 W m2
where the surface albedo is A.
Check: Units OK. Magnitude and sign OK.
Exposition: Negative sign means net inflow to the
surface, such as would cause daytime warming.
Longwave (IR) Sample Application

Upward emission of IR radiation from the Earths Find the net radiation at the surface in Vancouver,
surface can be found from the Stefan-Boltzmann re- Canada, at noon (standard time) on 22 Jun. Low clouds
lationship: are present with 30% coverage.
I = eIR SB T 4 (2.37) Find the Answer

Assume: Grass lawns with albedo A = 0.2 .
where eIR is the surface emissivity in the IR por- No other clouds.
tion of the spectrum (eIR = 0.9 to 0.99 for most sur- Given: L = 0.3
faces), and SB is the Stefan-Boltzmann constant (= Find: F* = ? Wm2
5.67x10 8 Wm2K4).
However, downward IR radiation from the at- Use = 64.1 from an earlier Sample Application.
Use eq. (2.35) to find the transmissivity:
mosphere is much more difficult to calculate. As an
Tr = [0.6+0.2sin](10.4L)
alternative, sometimes a net longwave flux is de- = [0.6 + 0.2sin(64.1)][1(0.40.3)]
fined by = [0.80](0.88) = 0.686
I* = I + I
(2.38) Use eq. (2.39) to find net IR contribution:
I* = b(1 0.6L)
One approximation for this flux is = (98.5 Wm2 )[1(0.60.3)] = 80.77 Wm2
Use eq. (2.40a):

I * = b (1 0.1 H 0.3 M 0.6 L ) (2.39)
F* = (1A) S Tr sin() + I*
= (10.2) (1366 Wm2) 0.686 sin(64.1)
where parameter b = 98.5 Wm2 , or b = 0.08 Kms1 + 80.77 Wm2
=( 674.36 + 80.77) Wm2 = 593.59 Wm2
in kinematic units.
Check: Units OK. Physics OK.
Net Radiation Exposition: The surface flux is only about 43% of that
Combining eqs. (2.33), (2.34), (2.35), (2.36) and at the top of the atmosphere, for this case. The negative
(2.39) gives the net radiation (F*, defined positive sign indicates a net inflow of radiation to the surface,
upward): such as can cause warming during daytime.
F * = (1 A) S Tr sin( ) + I * daytime (2.40a)
= I* nighttime (2.40b)
4PMBSSBEJBUJPOSFBDIJOH
&BSUITPSCJU

Actinometers
&M 8pNmpNm

Sensors designed to measure electromagnetic &BSUIT
SBEJBUJPO
radiative flux are generically called actinometers
or radiometers. In meteorology, actinometers are
usually oriented to measure downwelling or up-
welling radiation. Sensors that measure the differ-
ence between down- and up-welling radiation are

called net actinometers.
Special categories of actinometers are designed M N
to measure different wavelength bands: TIPSUXBWF MPOHXBWF

N N
USBOTQBSFODZ USBOTQBSFODZ
pyranometer broadband solar (short-wave) ir- PGHMBTT PGTJMJDPO
radiance, viewing a hemisphere of solid angle, XJOEPX XJOEPX
with the radiation striking a flat, horizontal plate
(Fig. 4.15). Figure 4.15
net pyranometer difference between top and The wavelength bands observed by pyranometers (shortwave
bottom hemispheres for short-wave radiation. radiation) and pyrgeometers (longwave radiation) depends on
the transparency of the windows used in those instruments.
pyrheliometer solar (short wave) direct-beam

A SCI. PERSPECTIVE Seek Solutions
radiation normal to a flat surface (and shielded
Most differential equations describing meteo- from diffuse radiation).
rological phenomena cannot be solved analytically. diffusometer a pyranometer that measures
They cannot be integrated; they do not appear in a only diffuse solar radiation scattered from air,
table of integrals; and they are not covered by the particles, and clouds in the sky, by using a device
handful of mathematical tricks that you learned in
that shades the sensor from direct sunlight.
math class.
But there is nothing magical about an analytical
pyrgeometer infrared (long-wave) radiation
solution. Any reasonable solution is better than
from a hemisphere that strikes a flat, horizontal
no solution. Be creative. surface (Fig. 4.15).
While thinking of creative solutions, also think
net pyrgeometer difference between top and
of ways to check your answer. Is it the right order bottom hemispheres for infrared (long-wave) ra-
of magnitude, right sign, right units, does it approach diation.
a known answer in some limit, must it satisfy some radiometer measure all wavelengths of radia-
other physical constraint or law or budget? tion (short, long, and other bands).
net radiometer difference between top and bot-
Example tom hemispheres of radiation at all wavelengths.
Find the irradiance that can pass through an atmo- spectrometers measures radiation as a func-
spheric window between wavelengths 1 and 2. tion of wavelength, to determine the spectrum of
radiation.
Find the Answer:
Approach: Integrate Plancks law between the Inside many radiation sensors is a bolometer,
specified wavelengths. This is the area under a por- which works as follows. Radiation strikes an object
tion of the Planck curve.
such as a metal plate, the surface of which has a coat-
Check: The area under the whole spectral curve
ing that absorbs radiation mostly in the wavelength
should yield the Stefan-Boltzmann (SB) law. Namely, band to be measured. By measuring the tempera-
the answer should be smaller than the SB answer, but ture of the radiatively heated plate relative to a non-
should increase and converge to the SB answer as the irradiated reference, the radiation intensity can be
lower and upper limits approach 0 and , respec- inferred for that wavelength band. The metal plate
tively. is usually enclosed in a glass or plastic hemispheric
chamber to reduce error caused by heat conduction
Methods: with the surrounding air.
Pay someone else to get the answer (Dont do this in Inside other radiation sensors are photometers.
school!), but be sure to check it yourself. Some photometers use the photoelectric effect,
Look up the answer in a Table of Integrals. where certain materials release electrons when
Integrate it using the tricks you learned in math
struck by electromagnetic radiation. One type of
class.
Integrate it using a symbolic equation solver on a
photometer uses photovoltaic cells (also called
computer, such as Mathematica or Maple. solar cells), where the amount of electrical energy
Find an approximate solution to the full equation. generated can be related to the incident radiation.
For example, integrate it numerically on a computer. Another photometric method uses photoresistor,
(Trapezoid method, Gaussian integration, finite dif- which is a high-resistance semiconductor that be-
ference iteration, etc.) comes more conductive when irradiated by light.
Find an exact solution for an approximation to the Other photometers use charge-coupled devic-
eq., such as a model or idealization of the physics. es (CCDs) similar to the image sensors in digital
Most eqs. in this textbook have used this approach. cameras. These are semiconductor integrated cir-
Draw the Planck curve on graph paper. Count the cuits with an array of tiny capacitors that can gain
squares under the curve between the wavelength
their initial charge by the photoelectric effect, and
bands, and compare to the value of each square, or to
the area under the whole curve. (We will use this ap-
can then transfer their charge to neighboring capaci-
proach extensively in the Thunderstorm chapter.) tors to eventually be read by the surrounding cir-
Draw the curve, and measure area with a planim- cuits.
eter. Simple spectrometers use different filters in front
Draw the Planck curve on cardboard or thick pa- of bolometers or photometers to measure narrow
per. Cut out the whole area under the curve. Weigh it. wavelength bands. Higher spectral-resolution spec-
Then cut the portion between wavelengths, & weigh trometers use interferometry (similar to the Mi-
again. chelson inteferometer described in physics books),
...and there are probably many more methods. where the fringes of an interference pattern can be
measured and related to the spectral intensities. B2. Access web cam camera images from a city,
These are also sometimes called Fourier-trans- town, ski area, mountain pass, or highway near you.
form spectrometers, because of the mathematics Use visible shadows on sunny days, along with your
used to extract the spectral information from the knowledge of solar azimuth angles, to determine the
spacing of the fringes. direction that the camera is looking.
You can learn more about radiation, including
the radiative transfer equation, in the weather-satel- B3. Access from the web the exact time from mili-
lite section of the Satellites & Radar chapter. Satel- tary (US Navy) or civilian (National Institute of Stan-
lites use radiometers and spectrometers to remotely dards and Technology) atomic clocks. Synchronize
observe the Earth-atmosphere system. Other satel- your clocks at home or school, utilizing the proper
lite-borne radiometers are used to measure the glob- time zone for your location. What is the time dif-
al radiation budget (see the Climate chapter). ference between local solar noon (the time when the
sun is directly overhead) and the official noon ac-
cording to your time zone. Use this time difference
to determine the number of degrees of longitude that
you are away from the center of your time zone.
Review
B4. Access orbital information about one planet
The variations of temperature and humidity that
(other than Earth) that most interests you (or a plan-
you feel near the ground are driven by the diurnal
et assigned by the instructor). How elliptical is the
cycle of solar heating during the day and infrared
orbit of the planet? Also, enjoy imagery of the planet
cooling at night. Both diurnal and seasonal heating
if available.
cycles can be determined from the geometry of the
Earths rotation and orbit around the sun. The same
B5. Access runway visual range reports from sur-
orbital mechanics describes weather-satellite orbits,
face weather observations (METARs) from the web.
as is discussed in the Satellites & Radar chapter.
Compare two different locations (or times) having
Short-wave radiation is emitted from the sun and
different visibilities, and calculate the appropriate
propagates through space. It illuminates a hemi-
volume extinction coefficients and optical thickness.
sphere of Earth. The portion of this radiation that is
Also search the web to learn how runway visual
absorbed is the heat input to the Earth-atmosphere
range (RVR) is measured.
system that drives Earths weather.
IR radiation from the atmosphere is absorbed at
B6. Access both visible and infrared satellite pho-
the ground, and IR radiation is also emitted from the
tos from the web, and discuss why they look differ-
ground. The IR and short-wave radiative fluxes do
ent. If you can access water-vapor satellite photos,
not balance, leaving a net radiation term acting on
include them in your comparison.
the surface at any one location. But when averaged
over the whole globe, the earth-atmosphere system
B7. Search the web for information about the sun.
is approximately in radiative equilibrium.
Examine satellite-based observations of the sun
Instruments to measure radiation are called ac-
made at different wavelengths. Discuss the struc-
tinometers or radiometers. Radiometers and spec-
ture of the sun. Do any of the web pages give the
trometers can be used in remote sensors such as
current value of the solar irradiance (i.e., the solar
weather satellites.
constant)? If so, how has it varied recently?
B8. Access from the web daytime visible photos of

the whole disk of the Earth, taken from geostation-
Homework Exercises ary weather satellites. Discuss how variations in the
apparent brightness at different locations (different
latitudes; land vs. ocean, etc.) might be related to re-
Broaden Knowledge & Comprehension flectivity and other factors.
(Dont forget to cite each web address you use.)
B9. Some weather stations and research stations re-
B1. Access a full-disk visible satellite photo image of
port hourly observations on the web. Some of these
Earth from the web. What visible clues can you use
stations include radiative fluxes near the surface.
to determine the current solar declination angle?
Use this information to create surface net radiation
How does your answer compare with that expected
graphs.
for your latitude and time of year.
B10. Access information from the web about how

color relates to wavelength. Also, how does the A4. For the date assigned from exercise A3, find:
range of colors that can be perceived by eye compare (i) mean anomaly
to the range of colors that can be created on a com- (ii) true anomaly
puter screen? (iii) distance between the sun and the Earth
(iv) solar declination angle
B11. Search the web for information about albedos (v) average daily insolation
and IR emissivities for substances or surfaces that
are not already listed in the tables in this chapter. A5(). Plot the local solar elevation angle vs. local
time for 22 December, 23 March, and 22 June for the
B12. Find on the web satellite images of either forest- following city:
fire smoke plumes or volcanic ash plumes. Compare a. Seattle, WA, USA
the intensity of reflected radiation from the Earths b. Corvallis, OR, USA
surface as it shines through these plumes with ear- c. Boulder, CO, USA
lier satellite photos when the plumes were not there. d. Norman, OK, USA
Use these data to estimate the extinction coefficient. e. Madison, WI, USA
f. Toronto, Canada
B13. Access photos and diagrams from the web that g. Montreal, Canada
describe how different actinometers are constructed h. Boston, MA, USA
and how they work. Also, list any limitations of i. New York City, NY, USA
these instruments that are described in the web. j. University Park, PA, USA
k. Princeton, NJ, USA
l. Washington, DC, USA
Apply m. Raleigh, NC, USA
(Students, dont forget to put a box around each an-
n. Tallahassee, FL, USA
swer.)
o. Reading, England
p. Toulouse, France
A1. Given distances R between the sun and planets
q. Mnchen, Germany
compute the orbital periods (Y) of:
r. Bergen, Norway
a. Mercury (R = 58 Gm)
s. Uppsala, Sweden
b. Venus (R = 108 Gm)
t. DeBilt, The Netherlands
c. Mars (R = 228 Gm)
u. Paris, France
d. Jupiter (R = 778 Gm)
v. Tokyo, Japan
e. Saturn (R = 1,427 Gm)
w. Beijing, China
f. Uranus (R = 2,869 Gm)
x. Warsaw, Poland
g. Neptune (R = 4,498 Gm)
y. Madrid, Spain
h. Pluto (R = 5,900 Gm)
z. Melbourne, Australia
i. Eris: given Y = 557 Earth years, estimate the
aa. Your location today.
distance R from the sun assuming a circular orbit.
bb. A location assigned by your instructor
(Note: Eris orbit is highly eccentric and steeply tilt-
ed at 44 relative to the plane of the rest of the solar
A6(). Plot the local solar azimuth angle vs. local
system, so our assumption of a circular orbit was
time for 22 December, 23 March, and 22 June, for the
made here only to simplify the exercise.)
location from exercise A5.
A2. This year, what is the date and time of the:
A7(). Plot the local solar elevation angle vs. azimuth
a. perihelion b. vernal equinox
angle (similar to Fig. 2.6) for the location from exer-
c. summer solstice d. aphelion
cise A5. Be sure to add tic marks along the resulting
e. autumnal equinox f. winter solstice
curve and label them with the local standard times.
A3. What is the relative Julian day for:
A8(). Plot local solar elevation angle vs. azimuth
a. 10 Jan b. 25 Jan c. 10 Feb d. 25 Feb
angle (such as in Fig. 2.6) for the following location:
e. 10 Mar f. 25 Mar g. 10 Apr h. 25 Apr
a. Arctic Circle b. 75N c. 85N d. North Pole
i. 10 May j. 25 May k. 10 Jun l. 25 Jun
e. Antarctic Circle f. 70S g. 80S h. South Pole
m. 10 Jul n. 25 Jul o. 10 Aug p. 25 Aug
for each of the following dates:
q. 10 Sep r. 25 Sep s. 10 Oct t. 25 Oct
(i) 22 Dec (ii) 23 Mar (iii) 22 Jun
u. 10 Nov v. 25 Nov w. 10 Dec x. 25 Dec
y. todays date z. date assigned by instructor
A9(). Plot the duration of evening civil twilight (dif- b. Using the average solar irradiance, plot the ra-
ference between end of twilight and sunset times) diative flux (reaching the Earths surface through a
vs. latitude between the south and north poles, for perfectly clear atmosphere) vs. latitude. Assume lo-
the following date: cal noon.
a. 22 Dec b. 5 Feb c. 21 Mar d. 5 May
e. 21 Jun f. 5 Aug g. 23 Sep h. 5 Nov A19(). For the city of exercise A1, plot the average
daily insolation vs. Julian day.
A10. On 15 March for the city listed from exercise
A5, at what local standard time is: A20. What is the value of IR absorptivity of:
a. geometric sunrise a. aluminum b. asphalt c. cirrus cloud
b. apparent sunrise d. conifer forest e. grass lawn f. ice
c. start of civil twilight g. oak h. silver i. old snow
d. start of military twilight j. urban average k. concrete average
e. start of astronomical twilight l. desert average m. shrubs n. soils average
f. geometric sunset
g. apparent sunset A21. Suppose polluted air reflects 30% of the incom-
h. end of civil twilight ing solar radiation. How much (W m2) is absorbed,
i. end of military twilight emitted, reflected, and transmitted? Assume an in-
j. end of astronomical twilight cident radiative flux equal to the solar irradiance,
given a transmissivity of:
A11. Calculate the Eq. of Time correction for: a. 0 b. 0.05 c. 0.1 d. 0.15 e. 0.2 f. 0.25
a. 1 Jan b. 15 Jan c. 1 Feb d. 15 Feb e. 1 Mar g. 0.3 h. 0.35 i. 0.4 j. 0.45 k. 0.5 l. 0.55
f. 15 Mar g. 1 Apr h. 15 Apr i. 1 May j. 15 May m. 0.6 n. 0.65 o. 0.7
A12. Find the mass flux (kgm2s1) at sea-level, giv- A22. What is the value of albedo for the following
en a kinematic mass flux (m s1) of: land use?
a. 2 b. 5 c. 7 d. 10 e. 14 f. 18 g. 21 a. buildings b. dry clay c. corn d. green grass
h. 25 i. 30 j. 33 k. 47 l. 59 m. 62 n. 75 e. ice f. potatoes g. rice paddy h. savanna
i. red soil j. sorghum k. sugar cane
A13. Find the kinematic heat fluxes at sea level, giv- l. tobacco
en these regular fluxes (Wm2):
a. 1000 b. 900 c. 800 d. 700 e. 600 A23. What product of number density times absorp-
f. 500 g. 400 h. 300 i. 200 j. 100 tion cross section is needed in order for 50% of the
k. 43 l. 50 m. 250 n. 325 o. 533 incident radiation to be absorbed by airborne volca-
nic ash over the following path length (km)?
A14. Find the frequency, circular frequency, and a. 0.2 b. 0.4 c. 0.6 d. 0.8 e. 1.0 f. 1.5 g. 2
wavenumber for light of color: h. 2.5 i. 3 j. 3.5 k. 4 l. 4.5 m. 5 n. 7
a. red b. orange c. yellow d. green
e. cyan f. blue f. indigo g. violet A24. What fraction of incident radiation is transmit-
ted through a volcanic ash cloud of optical depth:
A15(). Plot Planck curves for the following black- a. 0.2 b. 0.5 c. 0.7 d. 1.0 e. 1.5 f. 2 g. 3
body temperatures (K): h. 4 j. 5 k. 6 l. 7 m. 10 n. 15 o. 20
a. 6000 b. 5000 c. 4000 d. 3000 e. 2500
f. 2000 g. 1500 h. 1000 i. 750 j. 500 A25. What is the visual range (km) for polluted air
k. 300 l. 273 m. 260 n. 250 h. 240 that has volume extinction coefficient (m-1) of:
a. 0.00001 b. 0.00002 c. 0.00005 d. 0.0001
A16. For the temperature of exercise A15, find: e. 0.0002 f. 0.0005 g. 0.001 h. 0.002
(i) wavelength of peak emissions i. 0.005 j. 0.01 k. 0.02 l. 0.05
(ii) total emittance (i.e., total amount of emis-
sions) A26. (i) What is the value of solar downward direct
radiative flux reaching the surface at the city from
A17. Estimate the value of solar irradiance reaching exercise A5 at noon on 4 July, given 20% coverage of
the orbit of the planet from exercise A1. cumulus (low) clouds.
(ii) If the albedo is 0.5 in your town, what is the
A18(). a. Plot the value of solar irradiance reaching reflected solar flux at that same time?
Earths orbit as a function of relative Julian day.
(iii) What is the approximate value of net E10. a. On a clear day at your location, observe and
longwave radiation at that time? record actual sunrise and sunset times, and the du-
(iv) What is the net radiation at that time, given ration of twilight.
all the info from parts (i) - (iii)? b. Use that information to determine the day of
the year.
A27. For a surface temperature of 20C, find the c. Based on your personal determination of the
emitted upwelling IR radiation (W m2) over the length of twilight, and based on your latitude and
surface-type from exercise A20. season, is your personal twilight most like civil, mil-
itary, or astronomical twilight?
Evaluate & Analyze E11. Given a flux of the following units, convert it to
E1. At what time of year does the true anomaly
a kinematic flux, and discuss the meaning and/or
equal:
advantages of this form of flux.
a. 45 b. 90 c. 135 d. 180 e. 225
a. Moisture flux: g waterm2s1
f. 270 g. 315 h. 360
b. Momentum flux: (kgair ms1)m2s1
c. Pollutant flux: gpollutantm2s1
E2() a. Calculate and plot the position (true anom-
aly and distance) of the Earth around the sun for the
E12. a. What solar temperature is needed for the
first day of each month.
peak intensity of radiation to occur at 0.2 microm-
b. Verify Keplers second law.
eters?
c. Compare the elliptical orbit to a circular orbit.
b. Remembering that humans can see light only
between 0.38 and 0.74 microns, would the sun look
E3. What is the optimum angle for solar collectors
brighter or dimmer at this new temperature?
at your town?
E13. A perfectly black asphalt road absorbs 100% of
E4. Design a device to measure the angular diam-
the incident solar radiation. Suppose that its result-
eter of the sun when viewed from Earth. (Hint, one
ing temperature is 50C. How much visible light
approach is to allow the sun to shine through a pin
does it emit?
hole on to a flat surface. Then measure the width
of the projected image of the sun on this surface di-
E14. If the Earth were to cool 5C from its present
vided by the distance between the surface and the
radiative equilibrium temperature, by what percent-
pin hole. What could cause errors in this device?)
age would the total emitted IR change?
E5. For your city, plot the azimuth angle for appar-
E15(). Evaluate the quality of the approximation to
ent sunrise vs. relative Julian day. This is the direc-
Plancks Law [see eq. (a) in the A Scientific Perspec-
tion you need to point your camera if your want to
tive Scientific Laws the Myth box] against the
photograph the sunrise.
exact Planck equation (2.13) by plotting both curves
for a variety of typical sun and Earth temperatures.
E6. a. Compare the length of daylight in Fairbanks,
AK, vs Miami, FL, USA.
E16. Find the solar irradiance that can pass through
b. Why do vegetables grow so large in Alaska?
an atmospheric window between 1 = 0.3 m and
c. Why are few fruits grown in Alaska?
2 = 0.8 m. (See the A Scientific Perspective
Seek Solutions box in this chapter for ways to do
E7. How would Fig. 2.6 be different if daylight (sum-
this without using calculus.)
mer) time were used in place of standard time dur-
ing the appropriate months?
E17. How much variation in Earth orbital distance
from the sun is needed to alter the solar irradiance
E8(). Plot a diagram of geometric sunrise times
by 10%?
and of sunset times vs. day of the year, for your loca-
tion.
E18. Solar radiation is a diffuse source of energy,
meaning that it is spread over the whole Earth rath-
E9(). Using apparent sunrise and sunset, calculate
er than being concentrated in a small region. It has
and plot the hours of daylight vs. Julian day for your
been proposed to get around the problem of the
city.
inverse square law of radiation by deploying very
large mirrors closer to the sun to focus the light as
collimated rays toward the Earth. Assuming that
all the structural and space-launch issues could be Synthesize

solved, would this be a viable method of increasing (Dont forget to state and justify all assumptions.)
energy on Earth?
S1. What if the eccentricity of the Earths orbit
E19. The sine law for radiation striking a surface around the sun changed to 0.2 ? How would the
at an angle is sometimes written as a cosine law, seasons and climate be different than now?
but using the zenith angle instead of the elevation
angle. Use trig to show that the two equations are S2. What if the tilt of the Earths axis relative to the
physically identical. ecliptic changed to 45 ? How would the seasons
and climate be different than now?
E20. Explain the meaning of each term in eq. (2.21).
S3. What if the tilt of the Earths axis relative to the
E21. a. Examine the figure showing average daily ecliptic changed to 90 ? How would the seasons
insolation. In the summer hemisphere during the and climate be different than now?
few months nearest the summer solstice, explain
why the incoming solar radiation over the pole is S4. What if the tilt of the Earths axis relative to the
nearly equal to that over the equator. ecliptic changed to 0 ? How would the seasons and
b. Why are not the surface temperatures near the climate be different than now?
pole nearly equal to the temperatures near the equa-
tor during the same months? S5. What if the rotation of the Earth about its axis
matched its orbital period around the sun, so that
E22. Using Table 2-5 for the typical albedos, specu- one side of the Earth always faced the sun and the
late on the following: other side was always away. How would weather
a. How would the average albedo will change if and climate be different, if at all?
a pasture is developed into a residential neighbor-
hood. S6. What if the Earth diameter decreased to half of
b. How would the changes in affect the net radia- its present value? How would sunrise and sunset
tion budget? time, and solar elevation angles change?
E23. Use Beers law to determine the relationship S7. Derive eq. (2.6) from basic principles of geometry
between visual range (km) and volume extinction and trigonometry. This is quite complicated. It can
coefficient (m1). (Note that extinction coefficient be done using plane geometry, but is easier if you
can be related to concentration of pollutants and use spherical geometry. Show your work.
relative humidity.)
S8. What if the perihelion of Earths orbit happened
E24(). For your city, calculate and plot the noontime at the summer solstice, rather than near the winter
downwelling solar radiation every day of the year, solstice. How would noontime, clear-sky values of
assuming no clouds, and considering the change in insolation change at the solstices compared to now?
solar irradiance due to changing distance between
the Earth and sun. S9. What if radiative heating was caused by the
magnitude of the radiative flux, rather than by the
E25. Consider cloud-free skies at your town. If 50% radiative flux divergence. How would the weather
coverage of low clouds moves over your town, how or the atmospheric state be different, if at all?
does net radiation change at noon? How does it
change at midnight? S10. Linearize Plancks Law in the vicinity of one
temperature. Namely, derive an equation that gives
E26. To determine the values of terms in the surface a straight line that is tangent to any point on the
net-radiation budget, what actinometers would you Planck curve. (Hint: If you have calculus skills, try
use, and how would you deploy them (i.e., which di- using a Taylors series expansion.) Determine over
rections does each one need to look to get the data what range of temperatures your equation gives rea-
you need)? sonable answers. Such linearization is sometimes
used retrieving temperature soundings from satel-
lite observations.
S11. Suppose that Kirchhoffs law were to change

such that a = 1 e. What are the implications?
S21. Consider Beers law. If there are n particles per

S12. What if Wiens law were to be repealed, be- cubic meter of air, and if a vertical path length in
cause it was found instead that the wavelength of air is s, then multiplying the two gives the number
peak emissions increases as temperature increases. of particles over each square meter of ground. If
a. Write an equation that would describe this. You the absorption cross section b is the area of shadow
may name this equation after yourself. cast by each particle, then multiplying this times the
b. What types of radiation from what sources would previous product would give the shadowed area di-
affect the radiation budget of Earth? vided by the total area of ground. This ratio is just
the absorptivity a. Namely, by this reasoning, one
S13. What if the solar constant were even less con- would expect that a = nbs.
stant than it is now. Suppose the solar constant ran- However, Beers law is an exponential function.
domly varies within a range of 50% of its present Why? What was wrong with the reasoning in the
value, with each new value lasting for a few years previous paragraph?
before changing again. How would weather and cli-
mate be different, if at all? S22. What if the atmosphere were completely trans-
parent to IR radiation. How would the surface net
S14. What if the distance between sun and Earth radiation budget be different?
was half what it is now. How would weather and
climate be different, if at all? S23. Existing radiometers are based on a bolometer,
photovoltaic cell, or charge-coupled device. Design
S15. What if the distance between sun and Earth a new type of actinometer that is based on a differ-
was double what it is now. How would weather and ent principle. Hint, think about what is affected in
climate be different, if at all? any way by radiation or sunlight, and then use that
effect to measure the radiation.
S16. Suppose the Earth was shaped like a cube, with
the axis of rotation perpendicular to the ecliptic, and
with the axis passing through the middle of the top
and bottom faces of the cube. How would weather
and climate be different, if at all?
S17. Suppose the Earth was shaped like a narrow

cylinder, with the axis of rotation perpendicular to
the ecliptic, and with the axis passing through the
middle of the top and bottom faces of the cylinder.
How would weather & climate be different, if at all?
S18. Derive eq. (2.21) from the other equations in

this chapter. Show your work. Discuss the physical
interpretation of the hour angle, and what the effect
of truncating it is.
S19. Suppose that the Earths surface was perfectly

reflective everywhere to short-wave radiation, but
that the atmosphere absorbed 50% of the sunlight
passing through it without reflecting any. What
percentage of the insolation would be reflected to
space? Also, how would the weather and climate be
different, if at all?
S20. Suppose that the atmosphere totally absorbed

all short wave radiation that was incident on it, but
also emitted an exactly equal amount of short wave
radiation as it absorbed. How would weather and
climate be different, if at all?

Ch02 Radiat

Uploaded by

Copyright:

Available Formats

Ch02 Radiat

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Ch02 Radiat

Uploaded by

Copyright:

Available Formats

Copyright 2015 by Roland Stull.

Practical Meteorology: An Algebra-based Survey of Atmospheric Science

2 Solar & Infrared Radiation

Contents Solar energy powers the atmosphere. This ener-

Parameter a1 0.1996 d(Gm)3/2, where d is Earth

Practical Meteorology: An Algebra-based Survey

Figure 2.1 a (linear) & b (log-log)

the time between new moons is 29.5 days.) Because

d dp Repeating this calculation on a spreadsheet for several

where e = 0.0167 is eccentricity, and a = 149.457 Gm

(winter solstice, when nighttime is longest). The

where d is Julian day, and dr is the Julian day of the $BODFS

Sample Application Daily Effects

Use eq. (2.6): sin() =

Sample Application() Sample Application (continuation)

Find the Answer  

d = 355, 82, 173.

8 0.0 17.3 34.9

The local azimuth angle of the sun relative to

Geometric sunrise and sunset occur when the

Find the Answer:

 Because flow is associated with a direction, so is

for fluxes Fx, Fy, and Fz in the x, y, and z directions,

m ward) is written with a positive value of Fx, while

Fluxes are sometimes written in kinematic

Kinematic mass flux equals the wind speed, M. Ki-

(a) Use eq. (2.9):

Check: Units OK. Physics reasonable. = co (2.12)

Figure 2.8 Figure 2.9

HIGHER MATH Incremental Changes A SCIENTIFIC PERSPECTIVE Scientific

Conclusion: Eq. (a) is not a perfect relationship,

ter of the sphere:

where R is the radius from the center of the sphere, 

diation from a small sphere spreads out radially, it

So = 1366 (7) Wm2 . (2.17) 4VO

In kinematic units (based on sea-level density), the

According to the inverse-square law, variations

Sample Application where Cp is given under eq. (2.11).

Absorption, Reflection & Transmission Find the Answer:

Some substances such as dark glass are semi- (1368 Wm -2 ) 149Gm 2

Table 2-4. Typical infrared emissivities. E transmitted

Within the air, however, the process is a bit more

Etransmitted = Eincident e n b s (2.31a)

Beers law can be written using an extinction coef-

' ,m Define F* as the net radiative flux (positive

     where K and I are negative because they are

where Tr is a net sky transmissivity. A negative

Longwave (IR) Sample Application

I = eIR SB T 4 (2.37) Find the Answer

Use eq. (2.40a):

F * = (1 A) S Tr sin( ) + I * daytime (2.40a)

to measure different wavelength bands: TIPSUXBWF MPOHXBWF

pyrheliometer solar (short wave) direct-beam

B8. Access from the web daytime visible photos of

B10. Access information from the web about how

all the structural and space-launch issues could be Synthesize

S11. Suppose that Kirchhoffs law were to change

S21. Consider Beers law. If there are n particles per

S17. Suppose the Earth was shaped like a narrow

Find the Answer

Because flow is associated with a direction, so is

m ward) is written with a positive value of Fx, while

where R is the radius from the center of the sphere,

where K and I are negative because they are