Weather Station Siting: Effects On Phenological Models
Weather Station Siting: Effects On Phenological Models
9
WEATHER STATION SITING
Effects on Phenological Models
1. INTRODUCTION
100 degrees as the ice point. However, shortly after Celsius proposed his
temperature scale, others began to use degrees centigrade for the scale
defined with 100 gradations between the ice point and boiling point (i.e., the
reverse of the original Celsius scale). Because there were several
temperature scales in usage during the late 1800’s, Callendar recommended
a single practical temperature using the ice and boiling points of water with
standard platinum resistance thermometers for interpolation.
During the 20th century, several conferences on the international
temperature scale were held in 1927, 1948, 1968, and 1990 (i.e., ITS-27,
ITS-48, ITS-68, and ITS-90). The ITS-27 adopted the concepts of Callender,
with the addition of thermocouples being the standard instrument. The ITS-
48 changed the name of the scale from centigrade to Celsius to credit Celsius
for developing the scale and because the word centigrade represents a unit of
angle measure in French. The ITS-48 also defined the absolute temperature
scale as having exactly 100 gradations between ice point and the boiling
point of water. The ITS-68 fixed the Kelvin temperature scale to absolute
zero by setting 1.0 K = 1/273.16 of the triple point (i.e., the equilibrium
between the phases of ice, liquid water, and water vapor). The triple point of
water was set at exactly 273.16 K and the ice point was given the value of
273.15 K = 0°C. The boiling point was still set at 373.15 K = 100°C. They
also eliminated the use of the term degrees for the absolute temperature scale
(i.e., the unit is Kelvin rather than degrees Kelvin). Later, the ITS-90
reported the boiling point was really 99.975°C, so the definition of the
boiling point was eliminated from the absolute temperature scale. However,
at standard temperature and pressure, the Celsius scale still is defined as
going from 0°C at the freezing point to 100°C at the boiling point, and the
Fahrenheit scale goes from 32°F at the freezing point to 212°F at the boiling
point. Today, absolute temperature (K) and temperature in Celsius, where
(0°C = 273.15 K), are commonly used in science.
3. MEASUREMENT THEORY
3.1 Sensors
distance from the fluid reservoir. The liquid expands past the constriction,
and expands to a maximum temperature. On cooling, the liquid above the
constriction contracts leaving a vacuum barrier, while the liquid below the
constriction contracts into the reservoir. This leaves the upper end of the
liquid at the maximum temperature. The minimum thermometer has a small
rider, which is pushed by the surface forces of the meniscus of the liquid
within the capillary bore. The upper end of the rider is left at a minimum
temperature position when the liquid rises again.
Thermistors are made of sintered semiconductor materials (e.g.,
manganese, nickel, copper, iron, cobalt, and uranium oxides) that are pressed
into the thermistor form and aged to promote stability. The resistance to
electrical current decreases as the temperature rises and simple electronics
are used to monitor the resistance. The physics involved in resistance sensor
operation is described in Quinn (1983). With a thermistor, a curvilinear
relationship between resistance and temperature is needed to calibrate the
sensors. Older thermistors were somewhat unstable and, during their
lifetime, the resistance would drift requiring frequent calibration. However,
now thermistors are aged before calibration to minimize drift. Because the
bead thermistors are manufactured with their leads inside the powdered
material before they are pressed, whereas disk thermistors are manufactured
prior to the leads being sprayed or printed on the material, thermistors made
in the shape of beads are more stable than thermistors pressed into disk
shapes. A major advantage of thermistors is that they are small sensors.
The wire/metal resistors may be made very thin, but they still must have an
appreciable length (several millimeters at the minimum) to have a practical
level of resistance (tens of ohms) for temperature measurement. Thermistors
can be made as small as 0.03 mm.
3.2 Shielding
and painted white to reflect away most solar radiation. For electronic
sensors, white plastic shields like the Gill shield (Gill 1983) are often used
(Figure 2). Long wave radiation is unavoidable, because it is emitted by all
surfaces at terrestrial temperatures.
Figure 4.9-2. A standard Stevenson Screen, image by Jack Kelly Clark, © Regents of the
University of California, used with permission (left), and a Gill radiation shield (right).
rH Ts rRT (1)
Tt =
rR + rH
by assuming the net radiation and convective heat transfer are equal ((Rn = C)
C
as illustrated in Figure 3. When estimating air temperature, the goal is to
have Tt approach T by making rH much smaller than rR and/or by making Ts
very close to T T. To make rH much smaller than rR, one can maximize
convective heat transfer by making the sensor very small or by ventilating
the sensor. Generally, ventilation of a few meters per second is adequate.
When ventilating, the air should be pulled over the sensor to avoid heating
from the fans. In naturally ventilated Stevenson Screens or Gill shields,
Chapter 4.9: Weather Station Siting 351
errors can result when wind speed is low. Using double walls and roof and
painting the shelter white to reflect solar radiation makes Ts closer to T.
T
Figure 4.9-3. Drawing of the energy balance between a thermometer and its surrounding
radiation shield where Rn is the net radiation, C is convective heat transfer, σ is the Stefan-
Boltzmann constant, Ts, Tt, and T are the shield, thermometer, and air temperatures (K), ρ is
the air density, Cp is the specific heat of air at constant pressure, and rR and rH are resistances
to radiative and sensible heat transfer.
3.3 Height
develop in fewer days because they are exposed to the greater heat
accumulation than organisms grown under cooler conditions. The
accumulation of heat is called “physiological time” and “heat units” are a
measure of physiological time, which typically are better than “calendar
Figure 4.9-4. Temperature profile data from a walnut orchard at 0.5, 1.0, 2.0, 6.0, and 10.0 m
height prior to and during a radiation freeze night in a northern California valley.
Figure 4.9-5. Days and cumulative degree-days to mature cotton bollworm from neonate
larvae to adult stages for different temperature ranges (after Wilson and Barnett 1983).
Chapter 4.9: Weather Station Siting 353
time” for predicting days between phenological stages. For example, the
cotton bollworm development from neonate larvae to adult (Figure 5) is
closely related to the number of accumulated heat units regardless of the
large differences in number of days between the stages due to temperature
differences during growth. The cumulative number of degree-days (oD) to
develop is about the same, whereas the number of days to develop is greatly
different for the different growth temperatures.
Most organisms show little growth or development below a lower
threshold temperature (T TL) and often there is little increase in development
rate above an upper threshold (T TU), but the number of hours between TL and
TU provide a measure of heat units for estimating the days between
phenological stages. Each hour of time that the organism is exposed to
temperature between the two thresholds is called a “degree-hour”. Dividing
the accumulated degree-hours by 24 converts the heat units to degree-days,
which are commonly used in phenological models.
There are numerous papers on methods to estimate degree-days from
daily maximum and minimum temperature data and on how to determine
threshold temperatures and the cumulative degree-day requirement for
various crops and pests (Baskerville and Emin 1969; Allen 1976; Johnson
and Fitzpatrick 1977; Parton and Logan 1981; Kline et al. 1982; Snyder
1985; Wann et al. 1985; Reicosky 1989; de Gaetano and Knapp 1993;
Kramer 1994; Yin et al. 1995; Pellizzaro et al. 1996; Roltsch et al. 1999;
Snyder et al. 1999, 2001; Cesaraccio et al. 2001). In general, researchers
have attempted to model diurnal temperature changes using mathematical
models with a range of complexity (Zalom et al. 1983; Snyder et al. 1999).
Calculating the number of degree-days between the two phenological
stages of interest using selected lower and upper threshold temperatures and
all appropriate data sets best identifies threshold temperatures. Then the
mean cumulative degree-days (oDm) over all of the data sets are computed.
The same thresholds are used to calculate degree-days for each data set until
the cumulative degree-days adds up to oDm. Then the mean absolute
difference or root mean square error between the predicted and observed
number of days over all of the data sets is calculated and recorded. This
process is repeated for a range of threshold temperatures until the smallest
mean absolute difference or root mean square error is identified. The
threshold temperatures with the smallest difference between predicted and
observed days between stages and the corresponding oDm provide the best
prediction model.
354 Phenology: An Integrative Environmental Science
es − e r γ *( n ) (2)
To = T + − a
Δ + γ * ρC p Δ + γ *
where T is the air temperature (°C), es is the saturation vapor pressure (kPa)
at T,
T e is the vapor pressure (kPa), ra is the aerodynamic resistance (s/m), Rn
is the net radiation (W/m2), G is the soil heat flux density (W/m2), ρ is the air
density (g/m3), Cp is the specific heat of air at constant pressure (J/g K), Δ is
the slope of the saturation vapor pressure at air temperature (kPa/°C), γ is the
psychrometric constant (kPa/°C), and γ* = γ (rs/ra), where rs is the surface
resistance (s/m). The adiabatic part (i.e., the middle term on the right) is
directly proportional to the vapor pressure deficit (es – e), which is
approximately an exponential function of the air temperature. The diabatic
term (i.e., the right-hand term) depends on the external energy supply
(mainly net radiation). Therefore, a wet surface temperature will be warmer
than air temperature when the vapor pressure deficit is high enough that the
adiabatic term is bigger than the diabatic term.
6. SITING EFFECTS
The main indirect source of energy for heating the air is the sum of direct
and diffuse short-wave (solar) radiation minus that which is reflected away
Chapter 4.9: Weather Station Siting 355
Rn = G + λE + H + M (3)
balance over bare soil that is similar to the grass during the rainy season but
can be quite different when the soil dries and the grass continues to transpire.
Also, it is likely that differences in vegetation around the weather stations
contributed to differences in the degree-day accumulation.
The type of underlying surface around a weather station should be
selected depending on the purpose of the data. If the purpose is to
characterize phenological development of natural vegetation, the stations
should be sited in a natural setting that represents the conditions of the
vegetation being studied. However, if the data are used to model phenology
of irrigated crops, it is best to site the stations over an irrigated grass surface
that removes variability due to intermittent rainfall.
6.2 Fetch
Even when temperatures are recorded over the same standard surface,
sometimes big differences are recorded due to the surrounding environment.
For example, Figure 7 shows the corresponding temperatures recorded at
Torrey Pines and Miramar Naval Air Station near San Diego, California.
Torrey Pines is on the coast and Miramar is located about 30 km inland.
Clearly, although these sites are not far apart, there were big differences in
the temperature data with the coastal site having warmer temperatures when
values are low and the inland site having warmer temperatures when values
are high. If phenological models were to be applied in this region, many
stations would be needed to account the temperature effects of the ocean and
how they change with distance from the coast.
358 Phenology: An Integrative Environmental Science
Figure 4.9-6. Corresponding hourly temperature data (May-June 2000) measured on the west
edge of a grass field and 177 m east of the west edge of a large grass field surrounded by
desert near Indio, California. The prevailing wind was from the west.
Figure 4.9-7. Corresponding hourly temperature data (June-November 2000) from Torrey
Pines on the coast and Miramar Naval Air Station, which is 30 km inland.
In areas that are far from the ocean, there are often big differences in
temperature over similar underlying surfaces. For example, Figure 8 shows
the corresponding temperatures measured at the Indio CIMIS station and the
Vintage Country Club, which are only about 30 km apart in the below sea
Chapter 4.9: Weather Station Siting 359
Figure 4.9-8. Corresponding hourly temperature data (June - November 2000) measured at
the Vintage Country Club near Indio, California and the Indio CIMIS station. The wind speed
at the Vintage Country Club averaged about 8% of the Indio CIMIS site.
7. CONCLUSIONS
REFERENCES CITED
Allen, J. C., A modified sine wave method for calculating degree days, Environ. Ent., 5, 388-
396, 1976
Baskerville, G. L., and P. Emin, Rapid estimation of heat accumulation from maximum and
minimum temperatures, Ecology, 50, 514-517, 1969.
Cesaraccio, C., D. Spano, P. Duce, and R. L. Snyder, An improved model for degree-days
from temperature data, Int. J. Biometeorol., 45, 161-169, 2001.
de Gaetano, A., and W. W. Knapp, Standardization of weekly growing degree day
accumulations based on differences in temperature observation and method, Agric. For.
Meteorol., 66, 1-19, 1993.
Gill, G. C., Comparison testing of selected naturally ventilated solar radiation shields, Final
Report Contract # NA-82-OA-A-266, NOAA, St. Louis, 15 pp., 1983.
Horstmeyer, S., Building blocks – What goes on in a cubic meter of air?, Weatherwise, 54(5),
20-27, 2001.
Idso, S. B., R. D. Jackson, P. J. Pinter, R. J. Reginato, and J. L. Hatfield, Normalizing the
stress-degree-day parameter for environmental variability, Agric. Meteorol., 24, 45-55,
1981.
Jackson, R. D., S. B. Idso, R. J. Reginato, and P. J. Pinter, Canopy temperature as a crop
water stress indicator, Water Resour. Res., 13, 651-656, 1981.
Johnson, M. E., and E. A. Fitzpatrick, A comparison of methods of estimating a mean diurnal
temperature curve during the daylight hours, Arch. Meteorol. Geophys. Bioklimatol.,
Series B, 25, 251-263, 1977.
Chapter 4.9: Weather Station Siting 361