Tmy 2
Tmy 2
•
I
This report was prepared as an account of work sponsored by an agency of the United States government.
Neither the United States government nor any agency thereof, nor any of their employees, makes any
warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness,
or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would
not infringe privately owned rights. Reference herein to any specific commercial product, process, or
service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its
endorsement, recommendation, or favoring by the United States government or any agency thereof. The
views and opinions of authors expressed herein do not necessarily state or reflect those of the United
States government or any agency thereof.
Available from:
National Technical Information Service
U.S. Department of Commerce
5285 Port Royal Road
Springfield, VA 22161
(703) 487-4650
Information pertaining to the pricing codes can be found in the current issue of the following publications
which are generally available in most libraries: Government Reports Announcements and Index ( GRA
and~; Scientific and Technical Abstract Reports (STAR); and publication NTIS-PR-360 available from NTIS
at the above address .
.,..
f.J Printed with a renewable source ink on paper containing at least 50% wastepaper, including 20% postconsumer waste
DISCLAIMER
To distinguish between the old and new TMY data sets, the new TMY data sets
are referred to as TMY2s. TMY and TMY2 data sets cannot be used
interchangeably because of differences in time (solar versus local), formats,
elements, and units. Unless they are revised, computer programs designed for
TMY data will not work with TMY2 data.
The TMY2s are data sets of hourly values of solar radiation and meteorological
elements for a 1-year period. Their intended use is for computer simulations of
solar energy conversion systems and building systems to facilitate performance
comparisons of different system types, configurations, and locations in the United
States and its territories. Because they represent typical rather than extreme
conditions, they are not suited for designing systems to meet the worst-case
conditions occurring at a location.
The TMY2 data sets and this manual were produced by the National Renewable
Energy Laboratory's (NREL's) Analytic Studies Division under the Resource
Assessment Program, which is funded and monitored by the U.S. Department of
Energy's Office of Solar Energy Conversion.
iii
Acknowledgments
The authors would like to acknowledge NREL staff members Steve Wilcox,
Martin Rymes, Mary Anne Dunlap, Lisa Shertz, Craig Robben, Carol Riordan,
Dave Renne, Gene Maxwell, Tom Stoffel, Chris Cornwall, Eric Hammond,
Roland Hulstrom, Daryl Myers, and Chet Wells and National Climatic Data
Center staff member Marc Plantico for their contributions to the development of
the TMY2 data sets and user's manual.
Other individuals also reviewed NREL's plans to generate the TMY2 data sets
and provided valuable recommendations. This feedback early in the project
permitted efforts to be focused on maximizing the benefits of the TMY2s for
users. We are thankful for the efforts of these individuals, whose names and
affiliations are: Raymond Bahm (Raymond J. Bahm and Associates), William
Beckman (University of Wisconsin), Larry Degelman (Texas A&M University),
Nolan Doesken (Colorado State University), Randy Gee (Industrial Solar
Technology Corporation), Chris Gueymard (Florida Solar Energy Center), Doug
Hittle (Colorado State University), Michael Holtz (Architectural Energy
Corporation), Michael Kennedy (Ecotope), Ed Kern (Ascension Technology,
Inc.), Sandy Klein (University of Wisconsin), Jan Kreider (University of
Colorado), Hans Lund (Technical University of Denmark), Ken May (Industrial
Solar Technology Corporation), Dave Menicucci (Sandia National Laboratories),
John Schaefer (Consultant), Arvid Skartveit (Geophysical Institute, Norway),
Veronica Soebarto (Texas A&M University), Didier Thevenard (Watsun
Simulation Laboratory), Mike Thomas (Sandia National Laboratories), and Frank
Vignola (University of Oregon).
IV
Contents
Preface iii
Acknowledgments iv
Section 1 Overview 1
Typical Meteorological Year-A Description 1
NSRDB-Source of Data for the TMY2s 2
Methodology 4
TMY2 Station Classification 4
Data Elements 5
Where to Order 6
References 7
Section 2 Stations 9
Locations 9
NSRDB Classification 9
TMY2 Classification 9
V
Appendix A (Continued)
Quality Control 41
Calculation of Illuminance Data 42
Assignment of Source and Uncertainty Flags 42
References 44
vi
SECTION 1
Overview
This user's manual describes typical meteorological year (TMY) data sets derived
from the 1961-1990 National Solar Radiation Data Base (NSRDB). Based on
more recent and accurate data, these data sets are recommended for use in place of
earlier TMY data sets (NCDC 1981) that were derived from the 1952-1975
SOLMET/ERSATZ data base (SOLMET-Vol. 1 1978 and SOLMET-Vol. 2
1979). To distinguish between the two TMY data sets, the new TMY data sets are
referred to as TMY2s.
TMY and TMY2 data sets cannot be used interchangeably because of differences
in time (solar versus local), formats, elements, and units. Unless they are revised,
programs designed for TMY data will not work with TMY2 data.
Section 1 of the manual provides general information about the TMY2s and how
they were developed; Section 2 lists the stations and provides station identifying
information and classification; Section 3 details the contents of the TMY2 files
and provides the location in the hourly records of data values and their source and
uncertainty flags; Section 4 compares the TMY2s with 30-year data sets;
Appendix A provides a description of the procedures used to develop the TMY2s;
Appendix B provides a key for present weather elements; and Appendix C
contains a table of unit conversion factors for converting SI data to other units.
A TMY provides a standard for hourly data for solar radiation and other
meteorological elements that permit performance comparisons of system types
and configurations for one or more locations. A TMY is not necessarily a good
indicator of conditions over the next year, or even the next 5 years. Rather, it
represents conditions judged to be typical over a long period of time, such as 30
years. Because they represent typical rather than extreme conditions, they are not
suited for designing systems and their components to meet the worst-case
conditions occurring at a location.
Overview - 1
NSRDB-Source of Data for the TMY2s
The TMY2s were derived from the NSRDB, Version 1.1, which was completed in
March 1994 by the National Renewable Energy Laboratory (NREL). The NSRDB
contains hourly values of measured or modeled solar radiation and meteorological
data for 239 stations for the 30-year period from 1961-1990. A complete
description of the NSRDB and how it was produced is presented in its user's
manual (NSRDB-Vol. 1 1992) and the final technical report (NSRDB-Vol. 2
1995). The original version of the NSRDB, Version 1.0, was completed in August
1992. Version 1.1 corrects two types of minor errors in Version 1. 0 that affected
about 10% of the stations (Rymes 1994).
There are two types of stations in the NSRDB: primary (denoted by asterisks in
the station map in Figure 1-1) and secondary (denoted by dots in the station map
in Figure 1-1). The 56 primary stations measured solar radiation for a part (from 1
to 27 years) of the 30-year period. The remaining 183 stations, designated as
secondary stations, made no solar radiation measurements and therefore use
modeled solar radiation data that are derived from meteorological data, such as
cloud cover. Both primary and secondary stations are National Weather Service
stations that collected meteorological data for the period 1961-1990.
• Better model for estimating values (More than 90% of the solar radiation
data in both data bases are modeled.)
• More measured data, some of which is direct normal radiation
• lmproved instrument calibration methods
• Rigorous procedures for assessing quality of data.
A comparison of the old and new data bases provided an incentive for developing
the TMY2s. On an annual basis, 40% of the NSRDB and SOLMET/ERSATZ
stations are in disagreement for global horizontal radiation by more than 5%, with
some stations showing disagreement of up to 18% (Marion and Myers 1992). For
direct normal radiation, 60% of the NSRDB and SOLMET/ERSATZ stations are
in disagreement by more than 5%, with some stations showing disagreement of
up to 33%. Disagreement between the two data bases is even greater when
compared on a monthly basis.
An analysis of cloud cover data indicated little or no change for the two periods;
consequently, most of the disagreement for NSRDB and SOLMET/ERSATZ data
is attributed to differences in reconstructing the instrument calibrations and
differences in the solar radiation models (NSRDB-Vol. 2 1995).
Overview - 2
•• tf'.J>··
Figure 1-1. Map showing the 239 stations in the National Solar Radiation Data Base, whose data were used to derive the TMY2s
Because of differences in the data bases from which they were derived, the old
TMYs and the new TMY2s will differ. For some stations, the differences may be
minor, but other stations will have large differences.
Methodology
Except for a few changes to the weighting criteria, which accounts for the relative
importance of the solar radiation and meteorological elements, the TMY2s were
created using similar procedures that were developed by Sandia National
Laboratories (Hall et al. 1978) to create the original TMYs from the 1952-1975
SOLMET/ERSATZ data. Studies by Freeman (1979), Siurna, D' Andrea, and
Hollands (1984), and Menicucci and Fernandez (1988) have shown that this
procedure gives reasonable results. Sandia's procedure has also been adopted by
Siurna, D' Andrea, and Hollands (1984) for developing TMYs for Canada.
The Sandia method is an empirical approach that selects individual months from
different years from the period of record. For example, in the case of the NSRDB
that contains 30 years of data, all 30 Januarys are examined, and the one judged
most typical is selected to be included in the TMY. The other months of the year
are treated in a like manner, and then the 12 selected typical months are
concatenated to form a complete year.
The 12 selected typical months for each station were chosen from statistics
determined by using five elements: global horizontal radiation, direct normal
radiation, dry bulb temperature, dew point temperature, and wind speed. These
elements are considered the most important for simulation of solar energy
conversion systems and building systems.
For other elements in the TMY2s, the selected months may or may not be typical.
Cloud cover, which correlates well with solar radiation, is probably reasonably
typical. Other elements, such as snow depth, are not related to the elements used
for selection; consequently, their values may not be typical. Even though wind
speed was used in the selection of the typical months, its relatively low weighting
with respect to the other weighted elements prevents it from being sufficiently
typical for simulation of wind energy conversion systems.
Overview -4
achieve a minimum of 15 candidate months, data filling for periods of up to 47
hours were required. For some elements not required for the selection of the
typical meteorological months, the data are unfilled in the TMY2 data files. The
elements horizontal visibility, ceiling height, and present weather may be missing
for up to 2 consecutive hours for Class A stations and for up to 47 hours for Class
B stations. No data are missing for more than 47 hours, except for snow depth and
days since last snowfall for Colorado Springs, Colorado.
Data Elements
Table 1-1 shows the data elements in the TMY2 data files. These are the same
elements as for the 30-year NSRDB, except that illuminance and luminance
elements were added to support building energy analysis. The table includes
information by element and station classification to alert the user to the possibility
of missing data. Definitions of the elements and their units are provided in Table
3-2 of Section 3.
Data Completeness
Element Class A Class B
Extraterrestrial Horizontal Radiation 1 1
Extraterrestrial Direct Normal Radiation 1 1
Global Horizontal Radiation 1 1
Direct Normal Radiation 1 1
Diffuse Horizontal Radiation 1 1
Global Horizontal Illuminance 1 1
Direct Normal Illuminance 1 1
Diffuse Horizontal Illuminance 1 1
Zenith Luminance 1 1
Total Sky Cover 1 1
Opaque Sky Cover 1 1
Dry Bulb Temperature 1 1
Dew Point Temperature 1 1
Relative Humidity 1 1
Atmospheric Pressure 1 1
Wind Direction 1 1
Wind Speed 1 1
Horizontal Visibility 2 2, 3,4
Ceiling Height 2 2,3,4
Present Weather 2 2,3,4
Precipitable Water 1 1
Broadband Aerosol Optical Depth 1 1
Snow Depth 1 5
Days Since Last Snowfall 1 5
Notes:
1. Serially complete, no missing data.
2. Data may be present only every third hour.
3. Nighttime data may be missing.
4. Data may be missing for up to 47 hours.
5. Serially complete, except for Colorado Springs, CO.
Overview - 5
Where to Order
TMY2 data sets are available over Internet from NREL' s Renewable Resource
Data Center (RReDC). The Universal Resource Locator (URL) address of the
RReDC is "http://rredc.nrel.gov." Users should have World Wide Web (WWW)
browsing software, such as Mosaic or Netscape, to access the RReDC.
TMY2 data sets for all 239 stations may also be obtained on a CD-ROM. A
"Readme" file, which describes the contents, is included on the CD-ROM. The
CD-ROM may be ordered from:
Overview - 6
References
Freeman, T. L. (1979). Evaluation of the "Typical Meteorological Years" for
Solar Heating and Cooling System Studies. SERI/TR-8150-1. Golden, CO: Solar
Energy Research Institute.
Hall, I.; Prairie, R.; Anderson, H.; Boes, E. (1978). Generation of Typical
Meteorological Years for 26 SOLMET Stations. SAND78-1601. Albuquerque,
NM: Sandia National Laboratories.
Rymes, M. (1994). "Beyond the NSRDB Version 1.0." Proceedings of the 1994
Annual American Solar Energy Society Conference; June 25-30, 1994, San Jose,
CA; pp. 445-450.
Overview - 7
SECTION 2
Stations
There are 239 TMY2 stations for the United States and its territories. These are
the same stations as for the NSRDB, from which the TMY2 data sets were
derived. The stations are National Weather Service stations that collected
meteorological data for the period of 1961-1990. Table 2-1 lists the stations by
state or territory and provides information describing the station location and the
NSRDB and TMY2 classifications.
Locations
The station locations are described in Table 2-1 by the city and state name, the
station Weather Bureau Army Navy (WBAN) identification number, the latitude
and longitude in degrees and minutes, and the elevation in meters.
NSRDB Classification
Stations are classified with respect to being NSRDB primary (P) or secondary (S)
stations. The 56 primary stations measured solar radiation for a part (from 1 to 27
years) of the 30-year period of 1961-1990. The remaining 183 secondary stations
made no solar radiation measurements and therefore use modeled solar radiation
data that are derived from meteorological data such as cloud cover.
TMY2 Classification
This classification pertains to the amount of measured meteorological data
available for a station to select typical months to form the typical meteorological
year. Class A stations, of which there are 216, had a minimum of 15 candidate
months without more than 2 consecutive hours of missing data. For the 23 Class B
stations to achieve a minimum of 15 candidate months, data filling for periods of
up to 4 7 hours were required. For some elements not required for the selection of
the typical meteorological months, the data are unfilled in the TMY2 data files.
The elements horizontal visibility, ceiling height, and present weather may be
missing for up to 2 consecutive hours for Class A stations and for up to 47 hours
for Class B stations. No data are missing for more than 47 hours, except for snow
depth and days since last snowfall for Colorado Springs, Colorado.
Stations - 9
Table 2-1. Station Locations and Classifications
Stations - 10
Table 2-1. Station Locations and Classifications (Continued)
Stations - 11
Table 2-1. Station Locations and Classifications (Continued)
Stations - 12
Table 2-1. Station Locations and Classifications (Continued)
Stations - 13
Table 2-1. Station Locations and Classifications (Continued)
Stations - 14
Table 2-1. Station Locations and Classifications (Continued)
Stations - 15
SECTION 3
Each hourly record in the file contains values for solar radiation, illuminance, and
meteorological elements. A two-character source and uncertainty flag is attached
to each data value to indicate whether the data value was measured, modeled, or
missing, and to provide an estimate of the uncertainty of the data value.
Users should be aware that the format of the TMY2 data files is different from the
format used for the NSRDB and the original TMY data files.
File Convention
File naming convention uses the WBAN number as the file prefix, with the
characters TM2 as the file extension. For example, 13876.TM2 is the TMY2 file
name for Birmingham, Alabama. The TMY2 files contain computer readable
ASCII characters and have a file size of 1.26 MB.
File Header
The first record of each file is the file header that describes the station. The file
header contains the WBAN number, city, state, time zone, latitude, longitude, and
elevation. The field positions and definitions of these header elements are given in
Table 3-1, along with sample FORTRAN and C formats for reading the header. A
sample of a file header and data for January 1 is shown in Figure 3-1.
Hourly Records
Following the file header, 8760 hourly data records provide 1 year of solar
radiation, illuminance, and meteorological data, along with their source and
uncertainty flags. Table 3-2 provides field positions, element definitions, and
sample FORTRAN and C formats for reading the hourly records.
Each hourly record begins with the year (field positions 2-3) from which the
typical month was chosen, followed by the month, day, and hour information in
field positions 4-9. The times are in local standard time (previous TMYs based on
SOLMETIERSATZ data are in solar time).
Fi,gure 3-1. Sample file header and data in the TMY2 format for January 1
I I II ili'·
l,·, ,·i',·
'
1 ·1, 1-:
11 ,1 I
Table 3-1. Header Elements in the TMY2 Format
(For First Record of Each File)
Field
Position Element Definition
002- 006 WBANNumber Station's Weather Bureau Army Navy number (see Table 2-1)
008 - 029 City City where the station is located (maximum of22 characters)
031- 032 State State where the station is located (abbreviated to two letters)
034- 036 Time Zone Time zone is the number of hours by which the local standard
time is ahead of or behind Universal Time. For example,
Mountain Standard Time is designated -7 because it is 7 hours
behind Universal Time.
038- 044 Latitude Latitude of the station
038 N == North of equator
040- 041 Degrees
043 - 044 Minutes
046 - 053 Longitude Longitude of the station
046 W == West, E == East
048 - 050 Degrees
052 - 053 Minutes
056- 059 Elevation Elevation of station in meters above sea level
FORTRAN Sample Format:
( lX,A5,lX,A22,lX,A2,lX,I3,lX,Al,lX,I2,lX,I2,lX,Al,lX,I3,lX,I2,2X,I4 )
C Sample Format:
( %s %s %s %d %s %d %d %s %d %d %d)
Field
Position Element Values Definition
002- 009 Local Standard Time
002- 003 Year 61- 90 Year, 1961-1990
004 - 005 Month 1 - 12 Month
006 - 007 Day 1 - 31 Day of month
008 - 009 Hour 1 - 24 Hour of day in local standard time
010- 013 Extraterrestrial Horizontal 0 - 1415Amount of solar radiation in Wh/m2
Radiation received on a horizontal surface at the
top of the atmosphere during the 60
minutes preceding the hour indicated
014 - 017 Extraterrestrial Direct 0 - 1415 Amount of solar radiation in Wh/m2
Normal Radiation received on a surface normal to the
sun at the top of the atmosphere
during the 60 minutes preceding the
hour indicated
018 - 023 Global Horizontal Radiation Total amount of direct and diffuse
018 - 021 Data Value 0 - 1200 solar radiation in Wh/m 2 received on
022 Flag for Data Source A-H,? a horizontal surface during the 60
023 Flag for Data Uncertainty 0-9 minutes preceding the hour indicated
024 - 029 Direct Normal Radiation Amount of solar radiation in Wh/m2
024 - 027 Data Value 0- 1100 received within a 5.7° field of view
028 Flag for Data Source A-H,? centered on the sun during the 60
029 Flag for Data Uncertainty 0-9 minutes preceding the hour indicated
Field
Position Element Values Definition
030 - 035 Diffuse Horizontal Radiation Amount of solar radiation in Wh/m2
030- 033 Data Value 0- 700 received from the sky (excluding the
034 Flag for Data Source A-H,? solar disk) on a horizontal surface
035 Flag for Data Uncertainty 0-9 during the 60 minutes preceding the
hour indicated
036 - 041 Global Horiz. Illuminance Average total amount of direct and
036- 039 Data Value 0- 1300 diffuse illuminance in hundreds of lux
040 Flag for Data Source I, ? received on a horizontal surface
041 Flag for Data Uncertainty 0-9 during the 60 minutes preceding the
hour indicated.
0 to 1300 = 0 to 130,000 lux
042-047 Direct Normal Illurninance Average amount of direct normal
042 - 045 Data Value 0 - 1100 illuminance in hundreds of lux
046 Flag for Data Source I, ? received within a 5.7° field of view
047 Flag for Data Uncertainty 0-9 centered on the sun during the 60
minutes preceding the hour indicated.
0 to 1100 = 0 to 110,000 lux
048 - 053 Diffuse Horiz. Illuminance Average amount of illuminance in
048 - 051 Data Value 0- 800 hundreds of lux received from the sky
052 Flag for Data Source I,? (excluding the solar disk) on a
053 Flag for Data Uncertainty 0-9 horizontal surface during the 60
minutes preceding the hour indicated.
0 to 800 = 0 to 80,000 lux
054 - 059 Zenith Luminance Average amount of luminance at the
054 - 057 Data Value 0- 7000 sky's zenith in tens of Cd/m2 during
058 Flag for Data Source I, ? the 60 minutes preceding the hour
059 Flag for Data Uncertainty 0-9 indicated.
0 to 7000 =0 to 70,000 Cd/m2
060 - 063 Total Sky Cover Amount of sky dome in tenths
060 - 061 Data Value 0- 10 covered by clouds or obscuring
062 Flag for Data Source A-F phenomena at the hour indicated
063 Flag for Data Uncertainty 0-9
064 - 067 Opaque Sky Cover Amount of sky dome in tenths
064 - 065 Data Value 0-10 covered by clouds or obscuring
066 Flag for Data Source A-F phenomena that prevent observing the
067 Flag for Data Uncertainty 0-9 sky or higher cloud layers at the hour
indicated
068 - 073 Dry Bulb Temperature Dry bulb temperature in tenths of °C
068 - 071 Data Value -500 to 500 at the hour indicated.
072 Flag for Data Source A-F -500 to 500 =-50.0 to 50.0°C
073 Flag for Data Uncertainty 0-9
074 - 079 Dew Point Temperature Dew point temperature in tenths of
074 - 077 Data Value -600 to 300 °C at the hour indicated.
078 Flag for Data Source A-F -600 to 300 =-60.0 to 30.0°C
079 Flag for Data Uncertainty 0-9
080 - 084 Relative Humidity Relative humidity in percent at the
080- 082 Data Value 0- 100 hour indicated
083 Flag for Data Source A-F
084 Flag for Data Uncertainty 0-9
Field
Position Element Values Definition
085 - 090 Atmospheric Pressure Atmospheric pressure at station in
085 - 088 Data Value 700- 1100 millibars at the hour indicated
089 Flag for Data Source A-F
090 Flag for Data Uncertainty 0-9
091- 095 Wind Direction Wind direction in degrees at the hour
091 - 093 Data Value 0- 360 indicated. ( N = 0 or 360, E = 90,
094 Flag for Data Source A-F S = 180,W = 270 ). For calm winds,
095 Flag for Data Uncertainty 0-9 wind direction equals zero.
096 -100 Wind Speed Wind speed in tenths of meters per
096- 98 Data Value 0-400 second at the hour indicated.
99 Flag for Data Source A-F 0 to 400 = 0 to 40.0 mis
100 Flag for Data Uncertainty 0-9
101-106 Visibility Horizontal visibility in tenths of
101 - 104 Data Value 0 - 1609 kilometers at the hour indicated.
105 Flag for Data Source A-F,? 7777 = unlimited visibility
106 Flag for Data Uncertainty 0-9 0 to 1609 = 0.0 to 160.9 km
9999 = missing data
107 -113 Ceiling Height Ceiling height in meters at the hour
107 - 111 Data Value 0- 30450 indicated.
112 Flag for Data Source A-F,? 77777 = unlimited ceiling height
113 Flag for Data Uncertainty 0-9 88888 = cirroform
99999 = missing data
114 - 123 Present Weather See Present weather conditions denoted by
Appendix B a 10-digit number. See Appendix B
for key to present weather elements.
124 -128 Precipitable Water Precipitable water in millimeters at
124 - 126 Data Value 0-100 the hour indicated
127 Flag for Data Source A-F
128 Flag for Data Uncertainty 0-9
129 -133 Aerosol Optical Depth Broadband aerosol optical depth
129 - 131 Data Value 0- 240 (broad-band turbidity) in thousandths
132 Flag for Data Source A-F on the day indicated.
133 Flag for Data Uncertainty 0-9 0 to 240 = 0.0 to 0.240
134 -138 Snow Depth Snow depth in centimeters on the day
134 - 136 Data Value 0- 150 indicated.
137 Flag for Data Source A-F,? 999 = missing data
138 Flag for Data Uncertainty 0-9
139 -142 Days Since Last Snowfall Number of days since last snowfall.
139 - 140 Data Value 0- 88 88 = 88 or greater days
141 Flag for Data Source A-F,? 99 = missing data
142 Flag for Data Uncertainty 0-9
FORTRAN Sample Format:
(1X,4I2,2I4,7(I4,A1,I1l ,2(I2,Al,Ill,2(I4,Al,Ill,l(I3,Al,Ill,
1 (I4,Al, Il), 2 (I3 ,Al, Il), 1 (I4,Al, Ill, l (IS,Al, Ill, lOil, 3 (I3 ,Al, Ill,
1 (I2 ,Al, Ill l
C Sample Format:
(%2d%2d%2d%2d%4d%4d%4d%ls%ld%4d%ls%ld%4d%1s%ld%4d%1s%1d%4d%1s%1d%4d%1s
%1d%4d%ls%1d%2d%1s%1d%2d%1s%1d%4d%1s%1d%4d%1s%1d%3d%1s%1d%4d%1s%1d%3d
%1s%1d%3d%1s%ld%4d%1s%ld%51d%1s%1d%1d%1d%1d%1d%1d%1d%1d%1d%1d%1d%3d%1s
%1d%3d%ls%1d%3d%1s%1d%2d%1s%1dl
Note: For ceiling height data, integer variable should accept data values as large as 99999.
Missing Data
Data for some stations, times, and elements are missing. The causes for missing
data include such things as equipment problems, some stations not operating at
night, and a NOAA cost-saving effort from 1965 to 1981 that digitized data for
only every third hour.
Although both the NSRDB and the TMY2 data sets used methods to fill data
where possible, some elements, because of their discontinuous nature, did not
lend themselves to interpolation or other data-filling methods. Consequently, data
in the TMY2 data files may be missing for horizontal visibility, ceiling height, and
present weather for up to 2 consecutive hours for Class A stations and for up to 47
hours for Class B stations. For Colorado Springs, Colorado, snow depth and days
since last snowfall may also be missing. No data are missing for more than 47
hours, except for snow depth and days since last snowfall for Colorado Springs,
Colorado. As indicated in Table 3-2, missing data values are represented by 9's
and the appropriate source and uncertainty flags.
For the most part, the source and uncertainty flags in the TMY2 data files are the
same as the ones in NSRDB, from which the TMY2 files were derived. However,
differences do exist for data that were missing in the NSRDB, but then filled
while developing the TMY2 data sets. Uncertainty values apply to the data with
respect to when the data were measured, and not as to how "typical" a particular
hour is for a future month and day. More information on data filling and the
assignment of source and uncertainty flags is found in Appendix A.
Tables 3-3 through 3-6 define the source and uncertainty flags for the solar
radiation, illuminance, and meteorological elements.
Fla_g Definition
A Post-1976 measured solar radiation data as received from NCDC or
other sources
B Same as "A" except the global horizontal data underwent a
calibration correction
C Pre-1976 measured global horizontal data (direct and diffuse were
not measured before 1976), adjusted from solar to local time, usually
with a calibration correction
D Data derived from the other two elements of solar radiation using the
relationship, _global= diffuse + direct x cosine(zenith)
E Modeled solar radiation data using inputs of observed sky cover
(cloud amount) and aerosol optical depths derived from direct
normal data collected at the same location
F Modeled solar radiation data using interpolated sky cover and
aerosol optical depths derived from direct normal data collected at
the same location
G Modeled solar radiation data using observed sky cover and aerosol
optical depths estimated from _geo_graphical relationships
H Modeled solar radiation data using interpolated sky cover and
estimated aerosol optical depths
I Modeled illuminance or luminance data derived from measured or
modeled solar radiation data
? Source does not fit any of the above categories. Used for nighttime
values and missin_g data
Flag Definition
A Data as received from NCDC, converted to SI units
B Linearly interpolated
C Non-linearly interpolated to fill data gaps from 6 to 47 hours in
length
D Not used
E Modeled or estimated, except: precipitable water, calculated from
radiosonde data; dew point temperature calculated from dry bulb
temperature and relative humidity; and relative humidity calculated
from dry bulb temperature and dew point temperature
F Precipitable water, calculated from surface vapor pressure; aerosol
optical depth, estimated from geographic correlation
? Source does not fit any of the above. Used mostly for missing data
Flag Definition
1-6 Not used
7 Uncertainty consistent with NWS practices and the instrument or
observation used to obtain the data
8 Greater uncertainty than 7 because values were interpolated or
estimated
9 Greater uncertainty than 8 or unknown
0 Not definable
Results of these comparisons are shown in Figures 4-1 through 4-6. TMY2 values
for all stations are plotted against their respective 30-year average from the 1961-
1990 NSRDB. As indicated by the scatter of the data and the statistical
information at the top of the figures, agreement is better on an annual basis than
on a monthly basis. This is a consequence of cancellation of some of the monthly
differences when the monthly values are summed for the annual value. The
statistical information presented is the mean difference between the TMY2 value
and the 1961-1990 average and the standard deviation of the differences.
Table 4-1 provides 95% confidence intervals, determined as twice the standard
deviation of the differences between TMY2 and NSRDB values, for TMY2
monthly and annual solar radiation. The confidence intervals are given in units of
kWh/m2/day. Differences between TMY2 and NSRDB 30-year values should be
within the confidence interval 95% of the time.
Table 4-1. 95% Confidence Intervals for Monthly and Annual Solar Radiation
Comparison - 25
DIFFERENCE= TMY2 minus 1961-1990 Average
2
Mean = 0.026 kWh/m /day
2
Std Dev= 0.101 kWh/m /day
8
ca
+-'
C
2
·;::
6
0
I
<ti
~ 4
CJ
>,
:c
+-'
§ 2
~
C\J
>-
~
I- 0
0 2 4 6 8 10
10
DIFFERENCE= TMY2 minus 1961-1990 Average
2
Mean = 0.026 kWh/m /day
2
Std Dev= 0.030 kWh/m /day
8
"ctj.
+-'
C
0 6
N
·;::
0
I
ca
.0 4
0
a
ca:l
C
C 2
<(
~
~
I- 0
0 2 4 6 8 10
Comparison - 26
10
DIFFERENCE= TMY2 minus 1961-1990 Average
2
Mean = 0.057 kWh/m /day
2
Std Dev= 0.250 kWh/m /day
8
0 2 4 6 8 10
10
DIFFERENCE= TMY2 minus 1961-1990 Average
Mean = 0.057 kWh/m 2/day
Std Dev= 0.082 kWh/m 2/day
8
cij
E
,._ 6
0
z
t5
~ 4
0
ca::::,
C
C
<C 2
C\I
>-
~
I-
0
0 2 4 6 8 10
Comparison - 27
10
DIFFERENCE= TMY2 minus 1961-1990 Average
Mean = 0.040 kWh/m 2/day
Std Dev= 0.145 kWh/m 2/day
Q)
-
"O
::J
:g
_J 4
1!-
-
.c
C
0
~ 2
~
~
l-
o
0 2 4 6 8 10
10
DIFFERENCE= TMY2 minus 1961-1990 Average
Mean = 0.040 kWh/m 2/day
Std Dev= 0.046 kWh/m 2/day
8
0 2 4 6 8 10
Comparison - 28
Heating and Cooling Degree Day Comparisons
Degree days are the difference between the average temperature for the day and a
base temperature. If the average for the day (calculated by averaging the
maximum and minimum temperature for the day) is less than the base value, then
the difference is designated as heating degree days. If the average for the day is
greater than the base value, the difference is designated as cooling degree days.
Monthly and annual heating and cooling degree days (base 18.3°C) calculated
from the TMY2 data sets were compared with those for the same stations from
NCDC's data tape, "1961-1990 Monthly Station Normals All Elements." This
data tape includes temperature and degree day normals for about 477 5 stations in
the United States and its territories. The normals are averages computed by NCDC
for the period 1961-1990.
Results of these comparisons are shown in Figures 4-7 through 4-10. TMY2
values for all stations are plotted against their respective 30-year average from
NCDC's data tape. As seen for solar radiation, agreement is better on an annual
basis than on a monthly basis.
Table 4-2 provides 95% confidence intervals, determined as twice the standard
deviation of the differences between TMY2 and NCDC values, for TMY2
monthly and annual heating and cooling degree days. The confidence intervals are
given in units of degree days. Differences between TMY2 and NCDC 30-year
values should be within the confidence interval 95% of the time.
Table 4-2. 95% Confidence Intervals for Monthly and Annual Degree Days
Comparison - 29
1500
DIFFERENCE= TMY2 minus 1961-1990 Average
Cl) Mean = 1 .6 degree days
~ Std Dev= 21.2 degree days
0 1200
Q)
~
0)
Q)
0
0) 900
C:
~
Q)
::r:
>. 600
-
.c:
~
C:
0
~ 300
~
I-
Figure 4-7. Comparison of monthly heating degree days for NCDC and TMY2 data
12500
DIFFERENCE= TMY2 minus 1961-1990 Average
Mean = 19.3 degree days
Std Dev= 90.7 degree days
10000
7500
5000
2500
Figure 4-8. Comparison of annual heating degree days for NCDC and TMY2 data
Comparison - 30
DIFFERENCE= TMY2 minus 1961-1990 Average
CJ) Mean = -2.1 degree days
~ Std Dev = 11.0 degree days
Cl 400
Q)
Q)
~
C)
Q)
Cl
C) 300
C
0
0
()
~ 200
-
.c
~
C
0
~ 100
~
I-
Figure 4-9. Comparison of monthly cooling degree days for NCDC and TMY2 data
5000
DIFFERENCE= TMY2 minus 1961-1990 Average
Mean = -25.3 degree days
~ Std Dev = 48.8 degree days
~ 4000
Q)
~
C)
Q)
Cl
C)
3000
.!::
0
0
()
cii 2000
::i
C
C
<(
~ 1000
~
I-
Figure 4-10. Comparison of annual cooling degree days for NCDC and TMY2 data
Comparison - 31
References
Marion, W.; Wilcox, S. (1994). Solar Radiation Data Manual for Flat-Plate and
Concentrating Collectors. NRELffP-463-5607. Golden, CO: National Renewable
Energy Laboratory.
Comparison - 32
APPENDIX A
Sandia Method
The Sandia method is an empirical approach that selects individual months from
different years of the period of record. For example, in the case of the NSRDB
that contains 30 years of data, all 30 Januarys are examined and the one judged
most typical is selected to be included in the TMY. The other months of the year
are treated in a like manner, and then the 12 selected typical months are
concatenated to form a complete year. Because adjacent months in the TMY may
be selected from different years, discontinuities at the month interfaces are
smoothed for 6 hours on each side.
The Sandia method selects a typical month based on nine daily indices consisting
of the maximum, minimum, and mean dry bulb and dew point temperatures; the
maximum and mean wind velocity; and the total global horizontal solar radiation.
Final selection of a month includes consideration of the monthly mean and
median and the persistence of weather patterns. The process may be considered a
series of steps.
Step 1-For each month of the calendar year, five candidate months with
cumulative distribution functions (CDFs) for the daily indices that are closest to
the long-term (30 years for the NSRDB) CDFs are selected. The CDF gives the
proportion of values that are less than or equal to a specified value of an index.
Candidate monthly CDFs are compared to the long-term CDFs by using the
following Finkelstein-Schafer (FS) statistics (Finkelstein and Schafer 1971) for
each index.
Appendix A Procedures - 33
n
FS = (1/n) Loi
i =1
where
oi = absolute difference between the long-term CDF and the candidate
month CDF at xi
n = the number of daily readings in a month.
Four CDFs for global horizontal solar radiation for the month of June are shown
in Figure A-1. Compared to the long-term CDF by using FS statistics, the CDF
for June of 1981 compared the best and the CDF for June of 1989 compared the
worst. Even though it was not the best month with respect to the long-term CDF,
June of 1962 was selected for the TMY2. This was a consequence of additional
selection steps described in the following paragraphs.
Because some of the indices are judged more important than others, a weighted
sum (WS) of the FS statistics is used to select the 5 candidate months that have
the lowest weighted sums.
1.0
Long-Term (1961-1990)
C
0
:;:;
(.)
0.8 -e-- Best (1981)
C
::J
LL
---8- TMY2 Selected (1962)
-
.0
·;::
en
ci 0.4
Q)
>
~
::J
E
::J 0.2
0
0 2 4 6 8 10
Appendix A Procedures - 34
WS=lw-FS·
1 1
where
wi = weighting for index
Fsi = FS statistic for index.
Step 2-The 5 candidate months are ranked with respect to closeness of the
month to the long-term mean and median.
Step 3-The persistence of mean dry bulb temperature and daily global horizontal
radiation are evaluated by determining the frequency and run length above and
below fixed long-term percentiles. For mean daily dry bulb temperature, the
frequency and run length above the 67th percentile (consecutive warm days) and
below the 33rd percentile (consecutive cool days) were determined. For global
horizontal radiation, the frequency and run length below the 33rd percentile
(consecutive low radiation days) were determined.
The persistence data are used to select from the five candidate months the month
to be used in the TMY. The highest ranked candidate month from step 2 that
meets the persistence criteria is used in the TMY. The persistence criteria
excludes the month with the longest run, the month with the most runs, and the
month with zero runs.
Step 4--The 12 selected months were concatenated to make a complete year and
smooth discontinuities at the month interfaces for 6 hours each side using curve-
fitting techniques.
For the TMY2s, an index for direct normal radiation was added. This improves
the comparison between annual direct normal radiation for the TMY2s and the 30-
year annual average by about a factor of 2 (based on 20 geographically
representative NSRDB stations). When only global horizontal radiation is used for
the solar index, the TMY annual direct radiation values for the 20 stations were
within 4% (95% confidence level) of the 30-year annual average. Using both
global horizontal and direct radiation indices reduced the differences to 2%, with
no adverse effect on global horizontal radiation comparisons.
Appendix A Procedures - 35
Table A-1. Weightings for FS Statistics
Sandia NSRDB
Index Method TMY2s
Max Dry Bulb Temp 1/24 1/20
Min Dry Bulb Temp 1/24 1/20
Mean Dry Bulb Temp 2/24 2/20
Max Dew Point Temp 1/24 1/20
Min Dew Point Temp 1/24 1/20
Mean Dew Point Temp 2/24 2/20
Max Wind Velocity 2/24 1/20
Mean Wind Velocity 2/24 1/20
Global Radiation 12/24 5/20
Direct Radiation Not Used 5/20
Weightings for dry bulb and dew point temperature were changed slightly to give
more emphasis to dry bulb and dew point temperatures and less to wind velocity,
which is of less importance for solar energy conversion systems and buildings.
Neither of the TMY weightings is appropriate for wind energy conversion
systems.
The relative weights between solar and the other elements were not found to be
particularly sensitive. As an indicator, annual heating and cooling degree days
(base 18.3°C) were compared for the TMY2s and the 30-year period for the 20
stations. With the selected solar weighting of 50% (global and direct), annual
heating degree days for the TMY2s were within 5% (95% confidence level) of the
30-year annual average. As an extreme, reducing the solar weighting to zero only
reduced the differences to within 2Y2%. Differences between the TMY2 annual
averages and the 30-year averages for cooling degree days were within 9%, for
both 0% and 50% solar weightings.
As a consequence of adding the index for direct normal radiation, the persistence
check in Step 3 was modified to determine the frequency and run length below the
33rd percentile (consecutive low radiation days) for daily values of direct normal
radiation. This information, along with that for the other persistence indices, was
then used to select the month satisfying the persistence criteria.
El Chichon Years
The volcanic eruption of El Chichon in Mexico in March 1982 spewed large
amounts of aerosols into the stratosphere. The aerosols spread northward and
circulated around the earth. This phenomenon noticeably decreased the amount of
solar radiation reaching the United States during May 1982 until December 1984,
when the effects of the aerosols had diminished. Consequently, these months were
not used in any of the TMY2 procedures because they were considered not typical.
Appendix A Procedures - 36
Leap Years
TMY2 files do not include data for February 29. Consequently, data for February
29 were not used in leap year Februarys to determine their candidate month CDFs.
However, to maximize the use of available data, data for February 29 were
included for determining the long-term CDFs.
Class A stations are those stations whose 30-year meteorological data records
were the most complete and that had an adequate number ( 15) of candidate
months after eliminating any months with data missing for more than 2
consecutive hours. The minimum of 15 candidate months permitted completion
of 90% of the stations without extensive data filling. As indicated in Figure A-2,
as few as 15 candidate months yielded typical months that were within the range
of differences established by 25 or more candidate months when comparing
monthly values of direct normal for TMY2 months with monthly averages of
direct normal for the 1961-1990 period. This relationship was also found to be
true for global horizontal radiation and heating and cooling degree days.
Appendix A Procedures - 37
..-..
>,
C1l
1.0 - u
--E
"'O
C\I
0.8 -
()
8
fl
0
§
0
0
-- 0 0 0 8
0 § §
~
0
.c 0 0
0.6 - 8 8
3: 0
--
€l 0
i
0
.:.t:. 0 0
0
0 8
I
0
0 8 0 @ 0
u 0.4 -
~
0
() 8
~
Q)
,._ 0
€l fl fl
0 fl 8
i:5 0.2 -
I- I I
0
,._ 0
C1l - 8 0 0
Q) 0 R 8 "'
>- 0.0 fl
9
II I
(J
0
I
- 0 !l
8
0
§
0 0
(")
CJ)
-0.2 - 0
8
0
0 8 0
::I - 0 0 8 0
8 ~ §
€l
0
l,l
C 0 8 €1
-0.4 - €1 0 0
0
8 0
~ ~
-u
,._
Q) -0.6 -
- 0
0
8
0
0
0
0
0
0
0
0
€1
0
0
€1 g ~
0
8
Cl
- 0
0
0
0
C\J -0.8 - 0
0
8 0
>-
~
-
I- -1.0 I I
I I
I I I
I I
I I I I
10 12 14 16 18 20 22 24 26 28 30
Number of Candidate Months
Class B stations had more missing data than Class A stations, and the data were
filled for the index elements used to select the TMY2s. Other elements in Class B
TMY2s were not filled and may be missing. Table 1-1 on page 5 shows elements
that may have missing data values in TMY2 files for Class A and B stations.
Class A Stations. There are 216 Class A stations. Missing data for these
stations were accounted for in the following fashion:
Class B Stations. The NSRDB data from which the 23 Class B stations were
derived have substantially more missing data than the NSRDB data from which
the Class A stations were derived. This situation required filling missing data to
have sufficient candidate months from which to select typical months. The
additional missing data for the Class B stations resulted from such things as
equipment problems and the fact that some stations did not operate at night for
Appendix A Procedures - 38
some or all of the 30-year period. Criteria were relaxed for Class B stations to
permit filled data for periods of up to 47 hours to be used in determining the long-
term CDFs, and months were eligible to be candidate months if they had no
missing or filled data for periods greater than 47 hours. For Colorado Springs,
Colorado, the criteria were further relaxed to permit missing data for snow depth
and days since last snowfall.
Data-Filling Methods
The TMY2 data sets required filling some missing data that were not filled during
the development of the NSRDB. The NSRDB was made complete with respect to
solar radiation elements (NSRDB-Vol. 1 1992). This required NSRDB filling of
missing data, at least for daylight hours, for elements used to model solar
radiation, such as total and opaque sky cover, dry bulb temperature, relative
humidity, and atmospheric pressure.
For other meteorological elements, data were not filled in the NSRDB.
Consequently, to develop the TMY2s, missing data for dry bulb temperature
(nighttime), dew point temperature, and wind speed required data filling to
complete the selection of typical months. These elements, along with global
horizontal and direct normal radiation, were used to generate statistics to
determine the appropriate selection of typical months.
Data filling for TMY2 Class B stations was more extensive than for the Class A
stations. TMY2s for Class A stations were restricted to the selection of typical
months that had no more than 2 consecutive hours of data missing, whereas
Class B stations could have up to 47 consecutive hours of data missing.
Two-hour gaps in data records for Class A and Class B stations were filled by
linear interpolation, except for relative humidity, which was calculated based on
psychometric relationships (ASHRAE 1993) using measured or filled dry bulb
temperature and dew point temperature. For Class B stations, longer gaps from 3
to 47 hours were filled using filled data from the NSRDB if available; otherwise
TMY2 data filling-methods were used.
The NSRDB contains filled data for total and opaque sky cover, dry bulb
temperature, relative humidity, and atmospheric pressure. NSRDB data gaps up to
5 hours were filled by linear interpolation. Gaps from 6 to 47 hours were filled for
the above elements by using data from adjacent days for identical hours and then
by adjusting the data so that there were no abrupt changes in data values between
the filled and measured data. Many Class B stations did not operate for parts of
Appendix A Procedures - 39
the night and/or early morning and late afternoon. For these stations, NSRDB data
were filled from sunrise to sunset to allow model estimates of solar radiation.
However, nighttime data were not necessarily filled.
The TMY2 data sets used procedures to fill nighttime data and other data not
filled in the NSRDB. These procedures were used for total and opaque sky cover,
atmospheric pressure, dry bulb temperatures, dew point temperatures, relative
humidity, wind speed, precipitable water, broadband aerosol optical depth, snow
depth, and days since last snowfall. Data elements not filled are horizontal
visibility, ceiling height, and present weather.
Total and opaque sky cover, and atmospheric pressure were linearly interpolated
over any missing nighttime periods.
Nighttime dry bulb temperatures were linearly interpolated, and then the filled
values were adjusted to preserve nonlinearities, such as more rapid changes in
temperature near sunrise and sunset. These adjustments were based on average
diurnal profiles determined for each calendar month and appropriately scaled to
match the endpoints of the interpolation interval.
Missing wind speed data, for up to 47 hour gaps, were filled by the procedure used
to fill nighttime missing dry bulb temperatures-linear interpolation and then
adjustment of filled values based on average diurnal profiles determined for each
calendar month.
Missing wind direction and precipitable water, for up to 47 hour gaps, were
linearly interpolated. For calm winds, wind direction was set to zero (north).
Appendix A Procedures - 40
Broadband aerosol optical depth values in the TMY2s are daily values provided
by seasonal functions derived during the development of the NSRDB. The
seasonal functions are sinusoidal with respect to the day of the year and have peak
values occurring in the summer.
Snow depth and days since last snowfall data were available from the NSRDB for
all but Colorado Springs and a few stations at southern latitudes, such as Guam
and Puerto Rico. So much data were missing for Colorado Springs that no attempt
was made to fill the data, and missing data for the elements snow depth and days
since last snowfall were flagged as missing. For the southern latitude sites that do
not receive snow, snow depth was set to zero and days since last snowfall was set
to 88, meaning 88 or more days.
Quality Control
Data were checked before and after processing to ensure that data were
reasonable. NCDC provided information identifying some erroneous dew point
temperature data in Version 1.1 of the NSRDB, where dew point temperatures
exceeded dry bulb temperatures. During processing of the NSRDB data to
generate the TMY2s, dew point temperatures were checked to make sure they did
not exceed dry bulb temperatures. If they did, the dew point temperature was
calculated using relative humidity and dry bulb temperature, if available;
otherwise, the data were considered missing.
Post-processing checks revealed that some of the selected TMY2 months had
solar radiation values with obvious errors (diffuse radiation values were zero even
though global horizontal and direct normal radiation were a few hundred watt
hours). Consequently, these stations were reprocessed with the affected data
being excluded. The stations with months excluded during the reprocessing
because of erroneous solar data are: Boulder, Colorado (2/88, 3/85, 5/85, and
10/85);· Lake Charles, Louisiana (2/80); Caribou, Maine (4/78, 7/85, and 7/72);
Great Falls, Montana (10/89); Omaha, Nebraska (5/85, 5/89, and 11/81); Ely,
Nevada (6/89 and 9/88); Guam, Pacific Islands (1/88, 9/79, and 9/88); El Paso,
Texas (12/88); Midland, Texas (5/80 and 12/79); Salt Lake City, Utah (5/88, 8/80,
and 10/89); Lander, Wyoming (3/88 and 8/80).
Appendix A Procedures - 41
Calculation of llluminance Data
To facilitate lighting and energy analysis of buildings, hourly values for global
horizontal illuminance, direct normal illuminance, diffuse horizontal illuminance,
and zenith luminance were added to the TMY2 data sets. These elements were
calculated using luminous efficacy models developed by Perez et al. (1990).
Inputs to the models are global horizontal radiation, direct normal radiation,
diffuse horizontal radiation, and dew point temperature. The luminous efficacy in
terms of lumens per watt is determined as a function of sky clearness, sky
brightness, and zenith angle.
Usually, the source and uncertainty flags in the TMY2 data files are the same as
the ones in the NSRDB, from which the TMY2 files were derived. However,
differences do exist for data that were flagged missing in the NSRDB, but then
filled while developing the TMY2 data sets. Differences are also present for
illuminance and luminance data values that were not included in the NSRDB.
Uncertainty values apply to the data with respect to the time stamp of the data, and
not as to how "typical" a particular hour is for a future month and day. The
uncertainty values represent the plus or minus interval about the data value that
contains the true value 95% of the time.
The uncertainty assigned to modeled solar radiation data includes only the bias
error in the model and not the random error component, which could be several
times larger for partly cloudy skies. For partly cloudy skies, an hour can be
composed of large or small amounts of sunshine, depending on whether the sun is
mostly free of the clouds or occluded by the clouds. Consequently, modeled
hourly values may depart significantly from true values for partly cloudy skies.
The uncertainty assigned to modeled solar radiation data represents the average
uncertainty for a large number of model estimates (such as for a month). When
averaging large data sets, random errors tend to cancel, leaving only the bias error.
Appendix A Procedures - 42
The uncertainty of the model estimates are based on the evaluation presented by
Perez et al. (1990) for six test stations. To be conservative, the following model
mean bias errors for the stations with the largest errors were used:
The uncertainty of the illuminance data value was then determined as the root-
sum-square of the model uncertainty and solar radiation element uncertainty.
The use of the bias error, instead of bias and random error, is consistent with the
approach in the above paragraph concerning the assignment of uncertainty values
to modeled solar radiation elements. Consequently, it also has the same
implications. The assigned uncertainty is representative of the average uncertainty
for a large number of model estimates (such as for a month), but the actual
uncertainty of the individual modeled illuminance and luminance values is greater
than indicated.
For meteorological elements, relative uncertainties from the NSRDB were used.
These uncertainties do not portray a quantitative evaluation of the uncertainty of
the meteorological elements, but rather give relative uncertainties based on the
data and the manner in which they were derived (NSRDB-Vol. 1 1992).
The source and uncertainty flags for the solar radiation, illuminance, and
meteorological elements are presented in Tables 3-3 through 3-6 on pages 21 and
22.
Appendix A Procedures - 43
References
ASHRAE (1993). 1993 ASHRAE Handbook: Fundamentals. Atlanta, GA:
American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.
Hall, I.; Prairie, R.; Anderson, H.; Boes, E. (1978). Generation of Typical
Meteorological Years for 26 SOLMET Stations. SAND78-1601. Albuquerque,
NM: Sandia National Laboratories.
Perez, R.; Ineichen, P.; Seals, R.; Michalsky, J.; Stewart, R. (1990). "Modeling
Daylight Availability and Irradiance Components from Direct and Global
Irradiance." Solar Energy, 44(5), pp. 271-289.
Appendix A Procedures - 44
APPENDIX B
Field
Position Element Values Definition
114 Observation Indicator Oor9 0 = Weather observation made
9 = Weather observation not made, or missing
115 Occurrence of 0- 2, 4, 0 = Thunderstorm-lightning and thunder.
Thunderstorm, 6-9 Wind gusts less than 25.7 mis, and hail, if
Tornado, or Squall any, less than 1.9 cm diameter
1 = Heavy or severe thunderstorm-frequent
intense lightning and thunder. Wind gusts
greater than 25.7 mis and hail, if any, 1.9
cm or greater diameter
2 = Report of tornado or waterspout
4 = Moderate squall-sudden increase of wind
speed by at least 8.2 mis, reaching 11.3 mis
or more and lasting for at least 1 minute
6 =Waterspout ( beginning January 1984)
7 = Funnel cloud ( beginning January 1984 )
8 = Tornado ( beginning January 1984)
9 = None if Observation Indicator element
equals 0, or else unknown or missing if
Observation Indicator element equals 9
116 Occurrence of Rain, 0-9 0 = Light rain
Rain Showers, or 1 = Moderate rain
Freezing Rain 2 = Heavy rain
3 = Light rain showers
4 = Moderate rain showers
5 = Heavy rain showers
6 = Light freezing rain
7 = Moderate freezing rain
8 = Heavy freezing rain
9 = None if Observation Indicator element
equals 0, or else unknown or missing if
Observation Indicator element equals 9
Notes:
Light = up to 0.25 cm per hour
Moderate= 0.28 to 0.76 cm per hour
Heavy= greater than 0.76 cm per hour
117 Occurrence of Rain 0, 1, 0 = Light rain squalls
Squalls, Drizzle, or 3-9 1 = Moderate rain squalls
Freezing Drizzle 3 = Light drizzle
4 = Moderate drizzle
5 = Heavy drizzle
6 = Light freezing drizzle
7 = Moderate freezing drizzle
8 = Heavy freezing drizzle
9 = None if Observation Indicator element
equals 0, or else unknown or missing if
Observation Indicator element equals 9
Field
Position Element Values Definition
Occurrence of Rain Notes:
Squalls, Drizzle, or When drizzle or freezing drizzle occurs with
Freezing Drizzle other weather phenomena:
( continued ) Light = up to 0.025 cm per hour
Moderate= 0.025 to 0.051 cm per hour
Heavy= greater than 0.051 cm per hour
When drizzle or freezing drizzle occurs alone:
Light = visibility 1 km or greater
Moderate = visibility between 0.5 and 1 km
Heavy = visibility 0.5 km or less
118 Occurrence of Snow, 0-9 0 = Light snow
Snow Pellets, or Ice 1 = Moderate snow
Crystals 2 = Heavy snow
3 = Light snow pellets
4 = Moderate snow pellets
5 = Heavy snow pellets
6 = Light ice crystals
7 = Moderate ice crystals
8 = Heavy ice crystals
9 = None if Observation Indicator element
equals 0, or else unknown or missing if
Observation Indicator element equals 9
Notes:
Beginning in April 1963, any occurrence of ice
crystals is recorded as a 7.
119 Occurrence of Snow 0- 7, 9 0 = Light snow
Showers, Snow 1 = Moderate snow showers
Squalls, or Snow 2 = Heavy snow showers
Grains 3 = Light snow squall
4 = Moderate snow squall
5 = Heavy snow squall
6 = Light snow grains
7 = Moderate snow grains
9 = None if Observation Indicator element
equals 0, or else unknown or missing if
Observation Indicator element equals 9
120 Occurrence of Sleet, 0- 2, 0 = Light ice pellet showers
Sleet Showers, or Hail 4,9 1 = Moderate ice pellet showers
2 = Heavy ice pellet showers
4 = Hail
9 = None if Observation Indicator element
equals 0, or else unknown or missing if
Observation Indicator element equals 9
Notes:
Prior to April 1970, ice pellets were coded as
sleet. Beginning in April 1970, sleet and small
hail were redefined as ice pellets and are coded
as 0, 1, or 2.
Field
Position Element Values Definition
121 Occurrence of Fog, 0-9 O=Fog
Blowing Dust, or 1 = Ice fog
Blowing Sand 2 = Ground fog
3 = Blowing dust
4 = Blowing sand
5 = Heavy fog
6 = Glaze (beginning 1984)
7 = Heavy ice fog (beginning 1984)
8 = Heavy ground fog (beginning 1984)
9 = None if Observation Indicator element
equals 0, or else unknown or missing if
Observation Indicator element equals 9
Notes:
These values recorded only when visibility is
less than 11 km.
122 Occurrence of Smoke, 0- 7, 9 O=Smoke
Haze, Smoke and 1 = Haze
Haze, Blowing Snow, 2 = Smoke and haze
Blowing Spray, or 3 = Dust
Dust 4 = Blowing snow
5 = Blowing spray
6 = Dust storm (beginning 1984)
7 = Volcanic ash
9 = None if Observation Indicator element
equals 0, or else unknown or missing if
Observation Indicator element equals 9
Notes:
These values recorded only when visibility is
less than 11 km.
123 Occurrence of Ice 0- 2, 9 0 = Light ice pellets
Pellets 1 = Moderate ice pellets
2 = Heavy ice pellets
9 = None if Observation Indicator element
equals 0, or else unknown or missing if
Observation Indicator element equals 9