ages840424 yield models: an application to north dakota · 2017-02-24 · numerals with numeric...
TRANSCRIPT
United StatesDepartment ofAgriculture
StatisticalReportingService
StatisticalResearchDivision
Apri I 1984AGES840424
Methods for theEvaluation ofReal-Time WeatherData for use in CropYield Models:An Application toNorth DakotaIN COOPERATION WITH THE STATISTICSDEPARTMENT OF THE UNIVERSITY OF MISSOURI
Mark KaiserandJeanne L. Sebaugh
METHODS FOR THE EVALUATION OF REAL-TIME WEATHER DATA FOR USE IN CROP YIELDMODELS: AN APPLICATION TO NORTH DAKOTA. By Mark Kaiser and Jeanne L. Sebaugh,Ph.D., Statistical Research Division, Statistical Reporting Service, U. S.Department of Agriculture, Columbia, Missouri 65201; April 1984
ABSTRACT
A problem with the operational use of many crop yield models is a time delay inobtaining current values for the required weather data. Real-time estimates ofsubsequently published weather values (considered the standard) are often used.A procedure for evaluating these real-time estimates is developed through twoapproaches. One approach directly compares real-time data with the standard setusing statistical analysis. The other approach uses each set of data in a cropyield model and compares the results. Strengths and weaknesses of theseapproaches are illustrated using weather data and barley yield models for NorthDakota.
Key Words: Weather data, real-time data, data evaluation, crop yield models.
******************************************************************* This paper was prepared for limited distribution to the ** research community outside the U. S. Department of Agricul- ** ture. The views expressed herein are not necessarily those ** of SRS or USDA. *******************************************************************
ACKNOWLEDGMENTS
The authors wish to thank Asit Basu, Vikki French, Sharon LeDuc, Charles Perry,Clarence Sakamoto, Gary Sullivan and Wendell Wilson for their helpful commentsand criticism, Jean Sparks for typing numerous drafts of this report andJerry Wright for preparation of the cover and Figure 3.
i
METHODS FOR THE EVALUATION OF REAL-TIME WEATHER DATA FOR USEIN CROP YIELD MODELS: AN APPLICATION TO NORTH DAKOTA
Mark S. KaiserGraduate Research Assistant
andJeanne L. Sebaugh, Ph.D.
Mathematical Statistician
This research was conducted as a part of the AgRISTARS*Program and is part of the task identified as SRS-YES-IXin the 1984 AgRISTARS Program Plan.
*AgRISTARS is an acronym for Agriculture and ResourcesInventory Surveys Through Aerospace Remote Sensing. Itis a multi-agency research program to meet some currentand new information needs of the U. S. Department ofAgriculture.
AgRISTARS Staff ReportSRS-YES-IX (A) (84-04.1)
ii
TABLE OF CONTENTS
INTRODUCTION •••....•••.••EVALUATION OF A REAL-TIME DATA SOURCE
Tests of Central Tendency ••Measures of Reliability ••.••• 0 0
Investigation of DifferencesAdjusting CAB Estimates 0 • 0 0 ••• 0
EFFECT OF ESTIMATED WEATHER VALUES ON MODEL PREDICTIONCONCLUSIONS OF CAB EVALUATIONDISCUSSION OF EVALUATION PROCEDURE 0
LITERATURE CITED • 0 • 0 •• 0 •• 0
APPENDIX A. REAL-TIME WEATHER DATA SOURCES IN NORTH DAKOTA
1
237
........ 12•••••• 0 19
• 19
• • 20
• • 24
o • 26
28
APPENDIX BoPart 10
Part II 0
Part III.
DATA FOR CRDS 10, 20, 40, 50, 60,Corresponds to Figure 1 in Text 0
Corresponds to Figure 2 in Text 0
Corresponds to Table 4 in Text
80 AND 90 ...... 3535
• • 42• • 49
APPENDIX C. OBSERVED SIGNIFICANCE LEVELS FOR STATISTICALTESTS REFERRED TO IN TEXT . 0 0 •••••••••• 0 •• 0 •• 0 • 56
APPENDIX D. FORMULAE FOR MEASURES OF RELIABILITY FOR REAL-TIMEWEATHER VALUE ESTIMATES • 0 0 • 0 •• 0 0 ••••• 0 •• o .• 0 •• 57
APPENDIX E. ESTIMATES OF EQUATION PARAMETERS FOR REGRESSIONS OF CABPRECIPITATION ESTIMATES ON NCDC PRECIPITATION VALUES FOR 47MONTHS IN NORTH DAKOTA . 0 ••• 0 • 0 • 0 0 0 • 0 0 • 0 0 •• 0 • 0 •• 58
iii
1.
2.
3.
4.5.
6.
7.
8.
9.
LIST OF TABLES
North Dakota CRD values for average percent of statewideproduction (barley and spring wheat), 1970-1979, and averageJanuary temperature and annual precipitation, 1951-1980 ••••
Observed significance levels for parametric and non-parametrictests of data source effects on weather values •.•.•.••
Mean differences in temperature and precipitation databetween CAB and NCDC sources, 1979-1982
Measures of reliability for CAB weather value estimates .
Average monthly value differences (CAB values-NCDC values) foreach CRD, in descending order (n=44) •.......••
Measures of reliability for transformed CAB temperature estimates
Measures of reliability for transformed CAB precipitation estimates
Description of variables included in state level CEAS and YESbarley models for North Dakota .•.•..•••.•.
Comparison of CEAS and YES estimated barley yields with actualyields, 1979-1982 (bushels/acre) •.•....•.•..•.
Page
4
8
9
13
15
17
18
21
22
10. Measures of yield reliability for CEAS and YES barley modelsusing NCDC and CAB weather values over 4 years
A1. Potential real-time weather data sources for North Dakota.
23
29A2.A3.
A4.
North Dakota stations reported on DMONTH program
Potential methods for conversion of real-time station datato CD level values •••••..•.••.•••..
Dates CAB and WAOB weather value estimates were receivedat YES, 1981-1982 .••••.•...•••.•.••.•
iv
31
33
34
1.
2.
3.
LIST OF FIGURES
Plots of reported temperature and precipitation valuesover time by CRD •••••••••••••••Plots of differences between CAB and NCDC temperature andprecipitation values ••••••••••.
Flow chart for suggested paradigm in evaluation ofreal-time data sources .•••••.•••••••
v
Page
5
10
25
METHODS FOR THE EVALUATION OF REAL-TIME WEATHER DATA FOR USE INCROP YIELD MODELS: AN APPLICATION TO NORTH DAKOTA
by
Mark Kaiserand
Jeanne L. Sebaugh, Ph.D.
INTRODUCTION
A major objective of the AgRISTARS Program has been to develop and evaluatemodels which predict seasonal yields for a variety of crops and geographicalareas. These are often regression models which consist of one or more lineartrend terms (used as surrogate variables for technological advances), and anumber of weather related variables. For simple models the weather variablesare typically derived from monthly average temperature and total precipitationvalues. As many as twenty-five years of historical values may be required todevelop models and estimate regression coefficients for use in the current year.Often these historical weather data have been climate division (CD) monthlyvalues available from the National Climatic Data Center (NCDC) in Asheville,North Carolina. NCDC values are averages over many reporting weather stationswithin each CD. A problem with real-time use of these CD values has been adelay in obtaining published data from NCDC. However, there are sources ofreal-time weather data for some individual stations, and a number of proceduresfor converting these station data to climatic division values could potentiallybe used to provide input variables for operational crop-yield models. Given asource of real-time station data and a process by which those data are convertedto CD level values, a fundamental concern is the compatibility of weather valuesproduced by that process with the historical data base used to develop themodel. More specifically, the question becomes 'will the use of a real-timeprocess to estimate weather values in the current year affect the performance ofa crop-yield prediction modelr' Several approaches to addressing this questionhave been taken by different authors. Sakamoto et al. (1977) used plots andhistograms to compare monthly temperature and precipitation data from twosources, and computed correlations between crop-yield predictions produced byusing those data sources. Perry, Rogers and Ritchie (1982) used components ofvariance and other techniques to compare satellite-derived estimates of dailyminimum and maximum temperatures, precipitation, and solar radiation values toground measured values. Willmott and Wicks (1980) partitioned the root meansquare error obtained from regressing estimated monthly precipitation values onobserved values for a set of weather stations in California.
The objective of this study was to develop methods and a generalized procedurefor evaluating real-time weather data. Two approaches were used. One approachinvolved the direct comparison of real-time data with a standard set of values.The other approach involved an empirical investigation of the effect of esti-mated weather values on the predictive ability of two crop yield models.
1
Conducting an analysis of real-time data sources for one state was a convenientway to limit the project to a manageable geographic region. North Dakota wasselected because (1) Climatic Division (CD) and Crop Reporting District (CRD)boundaries coincide in North Dakota, (2) barley and spring wheat are major cropsin North Dakota, and models for these crops have been developed and evaluated aspart of the AgRISTARS Program (Motha, 1980; Barnett, 1981; LeDuc, 1981; Sebaugh,1981), and (3) data from both a real-time source and the historical data baseused to develop models are readily available for North Dakota. A description ofreal-time point-source data available in North Dakota and methods of convertingthose data to CD level values may be found in Appendix A. The criteria used toevaluate real-time data were:
1. Timeliness. Are the data available on a time scale that permitsweekly or monthly use of prediction models?
2. Availability and cost. Are data available in a usable form thatpermits economic derivation of model weather variables?
3. Accuracy and precision. How well do data estimate historical valuesused in model development?
The first two criteria are discussed in Appendix A. The remainder of thisreport focuses on criterion 3.
Historical weather values calculated by NCOC and used to develop crop yieldmodels are considered to be "true" values in this report. The use of estimatesin the current year which may deviate from these true values has a potentiallygreat, but unknown, impact on the crop yield forecast of a model. In the firstportion of this report we adopt the approach that any deviation of an estimatedweather value from the "true" values subsequently reported by NCDC is undesir-able. In the second portion of the report, we address the question of howdeviations from NCDC values affect the predictive quality of barley models forNorth Dakota.
EVALUATION OF A REAL-TIME DATA SOURCE
A detailed evaluation was conducted on North Dakota temperature and precip-itation values derived by NOAA, Climatic Impact Assessment Division, ClimateAssessment Division, Climate Assessment Branch (CAB) in Washington, D.C. Valuesderived by CAB are subjective estimates based on weather maps containing iso-lines (see Appendix A). This particular data source was selected for furtherstudy because the subjective estimates were similar to NCOC values in form(Le., one a verage temperature and total precipitation val ue per month per CRD).Thus CAB values could be used directly in crop yield models if they proved suf-ficiently accurate and precise. CAB estimates were available for the periodFebruary 1979 to December 1982. Differences between CRDs in the number ofvalues reported in subsequent tables are due to missing values. To facilitatepresentation of charts and figures, only CRDs 30 and 70 are included in the bodyof this report. Charts and figures for the other CRDs are included in AppendixB. In recent years CRD 30 reported the highest average area harvested for both
2
barley and spring wheat, while CRD 70 provided a geographic and climatic con-trast to CRD 30 (Table 1).
Tests of Central TendencyTemperature values followed a well defined annual cycle in all North DakotaCRDs, but precipitation values did not appear to follow a consistent pattern(Fig. 1). Also, the range in temperature values reported since 1979 by both CABand NCDC was consistent across CRDs, while the range in precipitation valueschanged between CRDs. Because CAB and NCDC values are reported for the sametime periods, they are correlated and should be compared in a pairwise manner.For this purpose a variable called YR_MONTH was constructed by juxtaposing yearnumerals with numeric month values. Thus 7901 was the value assigned to Janu-ary, 1979, 7902 was assigned to February, 1979, etc. This variable matches theCAB and NCDC values according to time of observation and was used to account fordifferences between years as well as differences between the months within ayear. If only the month were taken into consideration, values for the differentyears would be averaged for each month and these averages compared instead ofthe individual monthly observations. Because of the unstable annual cycle ofprecipitation, the use of the variable YR_MONTH as a blocking factor was prob-ably more important for precipitation values than for temperature values.
Because it was desired to compare NCDC and CAB values in a pairwise fashion, thedistribution of difference values (CAB-NCDC) was important. Assumptions ofnormality necessary for most parametric tests of location with pairwise datacould be met in only 2 of 9 CRDs, CRD 10 and CRD 40, for temperature and nonefor precipitation differences when tested with a Shapiro-Wilk statistic (Appen-dix C). Although the parametric procedure used is sometimes considered robustenough to ignore such conditions (Box 1953), both parametric and nonparametricanalysis of variance (ANOVA) procedures were used to test for differencesbetween data sources to protect against possible errors in the resulting infer-ence (Geary 1947; Tukey 1948). The variable YR MONTH was used as a blockingfactor, resulting in a randomized complete block design with two treatments perblock. The model tested for each CRD was:
Y·· 11 + 0.1, + B, + E: ••1J J 1J
where grand mean over both data sources
o.i = main fixed effects for data source
Bj = main effects for year and month (YR_MONTH), and
E:•. = independent error effects (N(O, 02).1J
The two-way ANOVA by ranks, or Friedman Ranks test, is the direct nonparam~tricanalog of the parametric randomized complete block design and was used to testthe null hypothesis that the data source effects are identical.
For only two data sources, these parametric and nonparametric procedures areequivalent to paired "t" and Wilcoxon matched-pairs signed-rank tests, respec-tively. However, they were employed because of their applicability when morethan two data sources are to be compared.
3
Table 1. North Dakota CRD values for average percent of statewideproduction (barley and spring wheat), 1970-1979, and average January
temperature and annual precipitation, 1951-1980
Average Percent of N.D. Average AverageProduction January Annual
CRD Bar1e S rin Wheat Tem . (C) Preci itation (mm)
10: Northwest 5.7 16.5 -14.9 400
20: North central 10.3 10.6 -16.3 427
30: Northeast 34.3 21.7 -16.4 458
40: West central 2.5 7.6 -10.9 418
50: Central 7.2 10.3 -14.8 444
60: East central 24.6 12.6 -15.2 490
70: Southwest 3.4 7.2 -11. 6 406
80: South central 2.2 4.8 -13.0 412
90: Southeast 9.9 8.7 -11. 2 494
4
Figure 1. Plots of reported temperature and precipitation values over time byCRn. Values from NCDC (Asheville) are represented by N; values from CAB by C.
CIltD=30II+III+III+III+III+ eII,+I---+++++++++++++++++++t+++++++++++++++++++++++++++---
1117111111188888888~~888888888888888888888888889999999999900000000UOOOll1ll111l1112222222222220000000011100000000wIII0000000001110000000001112345618901212345618~0121234S6789012123456189n12YEAR AND MONTH CYR_MONTH)
20
T 10EMPE 0RATu -10RE
-20
-30
CN
c
CN
cc ecCN
NN c
cc
c
cc CN c
CN
c
cCCN
c
cCN
C
c c
e
c
CN
CCNCC CN N
c
CCN
CIltD=30
IIISO +III+III90 +III60 +I NI N eI N30 + eN cI Neece N N N N NI N I\JNNCI C N NN ~ eN C eo + CC Ncecc ee ec c---+++++++++++++++++++t+++++++++++++++++++++++++++---~~~~~~~~~~~g8gggg88:gggYr?rrr~?r??r~~~~~~~~~~~~
0000000011100000000WIII000000000III0000000001112345618901212345618~0121234561890121234S6189012
YEAR AND MONTH CYR_MONTH)
PRE 120CIpITATIoN
NNN NN
N
cN
CN
eCNC
C CN CNC NC C C
N
ec
N
CCC
r.N
CNC
5
Figure 1. Continued
Ct.lO=70
+II eI+III+I---+++++++++++++++++++++++++++++++++++++++++++++++---
77777777777A8888888~~8888888888888888888888BAB89999999999900000000W0001111111111112222222222220000000011100000000W111000000000111000000000111234567890121234S678~0121234S67890121234S67890 2
YEAR AND ~ONTH CYR_MONTH)
II C?O + NC, C CII
10 + C CI NII C
0 +I CN NNII
eN ee NN l;;e e
'" cN
~MPE~ATU -10RE
-20
-30
eN
eN
CN
cce N
e
ee eN eee e
e
CNCe e
NeeN
NeNe e
e
e
ee
C1o(0=70
III
150 +I CI NI
120 + NIII C
90 + N NI N, N C, N C
pRECIpITATI 60 + N Co I NC N NNN I e
I C N N eee N e30 + C N C N CI N NC ce N N NC e ee
I e ~,(,. CNN C NI C NN NNNN C ee C C N eco + C Nceee ceee e---++~++++++++++++++++++++++++++++++++++++++++++++---~~~~~~~~~~~6gg3S8gg:3Sgr~~~~~~~~?~?~~~~~~~~~~~~
0000000011100000000WIII000000000111000000000111234S67890121234S678~0121234S6789012123456789012
YEAR AND ~ONTH (YR_MONTH)
6
Both the parametric and nonparametric procedures were run for each CRD, anddifferences between data sources were tested at the 0.05 and 0.10 significancelevels, respectively. Significant data source effects were found in 8 out of 9CRDs for temperature and in all CRDs for precipitation with both the parametricand nonparametric procedures (Table 2).
Difference values were calculated for each month (by CRD) by subtracting theNCDC value from the CAB value. The mean difference in temperature values rangedover CRDs from 0.12 Co to 0.98 Co, while mean precipitation differences rangedfrom -19.9 rom to -5.3 rom(Table 3). Values reported by CAB appeared to over-estimate true temperature values (mean differences positive for all CRDs), andunderestimate true precipitation values (all mean differences negative). Plotsof temperature and precipitation differences (Figure 2) also showed a fairlyconsistent overestimation of temperature values by CAB. Precipitation valuesgenerally were slightly underestimated by CAB early and late in the year andgrossly underestimated during the summer months. In addition, difference valuesfarthest from 0 (thus representing the worst predictions by CAB) were consis-tently the same year and month across CRDs for temperature but varied forprecipitation; the least accurate CAB temperature estimate in the 47 months ofdata was the June 1981 value, and this was true for all CRDs (see Figure 2 andAppendix B). The least accurate CAB precipitation estimates varied from CRD toCRD but generally occurred between May and August.
Measures of Reliability
The two ANOVA procedures (parametric and nonparametric) used to test for differ-ences between CAB and NCDC values were both tests of central tendency, that is,they determined if the CAB and NCDC values were centered around the same pointafter taking into account differences between years and montp' problem withthese procedures is that overestimates and underestimates may counteract eachother such that the mean difference between CAB and NCDC values is close to zerobut the CAB value is not a good estimate of the NCDC value in any given year andmonth. Thus, a number of poor estimates may "cancel" each other's effect on theoutcome of tests of location, yet alter the performance of a crop-yield model ifthey are used as real-time weather values. It is desirable that CAB weathervalues be good estimates of the true value not only on the average but also inany given year and month. The closer a set of real-time weather values is tothe set of true values when compared retrospectively in a pairwise fashion, themore reliable that real-time data source is considered to be. A set of descrip-tive statistics, originally developed for the evaluation of crop yield models,is also appropriate for the comparison of a real-time data source to a standard.These statistics are called indicators of reliability, are described in Wilsonand Sebaugh (1981), and formulas are given in Appendix D. In this report, thefollowing indicators of reliability are used:
BiasRelative BiasRoot Mean Square ErrorRelative Root Mean Square ErrorStandard DeviationRelative Standard DeviationPercent of Absolute Difference Values Greater than 1°C or 10 romLargest Positive Difference
7
Table 2. Observed significance levels for parametric and non-parametrictests of data source effects on weather values
Tern PreclCRD Parametric Parametric
10 .0001 .0001 .0407 .0065
20 .0001 .0001 .0001 .0001
30 .0001 .0001 .0001 .0001
40 .0001 .0023 .0001 .0001
50 .0008 .0007 .0001 .0001
60 .0172 .0007 .0001 .0001
70 .0001 .0001 .0001 .0001
80 .5445 .5430 .0001 .0002
90 .0001 .0001 .0011 .0166
8
Table 3. Mean differences in temperature and precipitationdata between CAB and NCDC sources, 1979-1982
Mean Difference (CAB-NCDC) Standard DeviationCRD n Temperature n Precipitation Temperature Precipitation
(Co) (mm) Co mm10 47 0.98 47 -5.3 .787 17.07
20 45 0.93 45 -19.9 .688 25.11
30 46 0.82 47 -18.7 .803 22.56
40 46 0.56 47 -10.9 1.114 18.71
50 47 0.43 47 -17.3 .791 22.07
60 47 0.36 47 -10.4 .968 17.08
70 47 0.70 47 -10.4 1.091 15.44
80 47 0.12 47 -6.9 1.15 10.37
90 44 0.67 45 -8.2 1.125 16.93
9
Figure 2. Plots of differences between CAB and NCDC temperature and precipitationvalues (CAB-NCDC) over time by CRD.
CII/()=30
TE"'1p
•DIF'F'EQF.~CES
5.0
2.5
0.0
-2.5
-5.0
II+II AII+I A A A AI A A AI AA A A A AA A A A AA AAI A A A A A A A A A A A A+------A-A-------A--------------A---A-------A---A----I A A AIII+IIII+I---+++++++++++++++++++++++++++++++++++++++++++++++---
777777777778888888~H888888888888888888888888888999999999990000000W000011111l1l111122222222?22200000000111000000010111000000000111000000000111234567890121234567a~012123456789012123456789012yEAR AND ~ONTH (YR~~ONTH)
CII/O=30
I AI A A A-75 + AI
I-100 +
PRECIp•DIF'F'ERENCES
25 +II AI A A A
o +----------------A---~---A------A---AA-A--A-----A----I A A A AA AI A A A AA AAAAI A A A A-25 + AI A AI AI A A-50 + A A
I---+++++++++++++++++++++++++++++++++++++++++++++++---
777777777778888888~~88888888888888888888888q888999999999990000000VOOOOlll1111111l1222222222222000000001110000000WOll1000000000111000000000111234567890121234567~~012123456789012123456789012
yEAR AND ~ONTH (YR_~ONTH)10
Figure 2. Continued
ChtD=70A
TEMp•ofF"ERENCE5
5.0
2.5
0.0
-2.5
-5.0
+IIII A+II A A AA A AI AA A A AA A AA A A A A AAA A AAAAA.--__-AA-----_-~ --~---A---A---A----AA---A---------I A A A AI A AII+IIII+I---+++++++++++++++++++++++++++++++++++++++++++++++---
777777777778888A88~H888888888A8888888888A88A888999999999990000000wOOOOlllllllllll1222222222222OOOOOOOOlllOOOOOOOWOIIlOOOOOOOOOlllOOOOOOOOOlll234567890121234567~~012123456789012123456789012yEAR AND MONTH (YR_MONTH)
ChtD=70I•+II A AI A A+---AA---------------~-AA-AAA---------A-------A-AA---I A A AAAA AAA A AAA AI AAI+III+III+III+I---+++++++++++++++++++++++++++++++++++++++++++++++---
777777177778888888~H88888888888888888888888A888999999999990000000wUOOOllllllllll11222222222222OOOOOOOOll10000000WOIIIOOOOOOOOOl1100000000011123456789012123456709012123456789012123456789012yEAR AND MONTH (YR~MONTH)
p 25RECI 0p•D -25IF"FE -50RENC -75ES
-100
AA AAA
A
A
A
A A
A
A
11
Largest Negative DifferenceDirection of Change from Previous ValuePearson Correlation Coefficient
Measures of reliability are presented in Table 4. Overestimation of temperaturevalues and underestimation of precipitation values by CAB are again seen in theconsistently positive and negative bias terms, respectively. It should be notedthat many of the quantities presented in Table 4 were calculated on the differ-ence values (CAB-NCDC), not the original CAB estimates. Thus while the varianceof CAB precipitation estimates should be large to reflect the true variance inprecipitation, the variance of the set of precipitation differences should besmall for a high level of reliability. The calculated measures indicate thatCAB temperature estimates were more reliable predictors of the true values thanwere CAB precipitation estimates. This difference in reliability is seen in thelower relative bias, relative RMSE, and relative standard deviation of tempera-ture differences, as well as the percent of YR MONTHS for which the direction ofchange was correct (Table 4). The percent of YR MONTHS for which absolutedifference values were greater than the specified limits of 1° C or 10 mm aver-aged 43.4% and 35.4% across the 9 CRDs respectively. Although the Pearson cor-relation coefficients between CAB and NCDC values also indicated that CAB tem-perature estimates were better predictors of NCDC values than were CAB precipi-tation estimates, the correlation coefficient is not a sensitive measure ofreliability; the correlation coefficient measures the tendency of two sets ofvalues to vary in the same fashion, but not necessarily on the same scale. Thusthe correlation coefficient for CAB prcipitation estimates is fairly high eventhough these estimates have extremely large relative standard deviations and areoff by more than 10 mm in about 45% of the observations (Table 4, Appendix B).
Investigation of Differences
An important aspect of studying temperature and precipitation difference values(CAB value-NCDC value) was to determine if the observed statistical differencebetween CAB and NCDC estimates could be explained. If so, these explanationsmight indicate consistent flaws in the CAB process. Consistent flaws may, ofcourse, be more easily corrected than inconsistent flaws. The search, then, wasfor systematic patterns in the dev iation of CAB val ues from the "true" NCDCvalues.
In order to investigate whether the degree of difference in the weather valueswas related to geographic location, a randomized complete block ANOVA was usedto test for significant CRD effects. The model tested was:
y .. == 11 + C\ + S. + E: •. ,1J J 1Jwhere 11 == grand mean difference over all CRDs,
ai == main fixed effects for CRD,
S' == main effects for year and month (YR_MONTH) , andJ
E:' • == the independent error effects N(O, rl-).1J
12
Table 4. Measures of reliability for CAB weather values estimates
Measure
Relative Bias = RB (%)
Relative Standard Deviation = RSD (%)
itation
0.82 -18.69 0.70 -10.36
19.2 -44.8 11.2 -30.51.14 29.12 1.29 18.46
26.8 69.8 20.6 54.3
0.79 22.33 LOa 15.28
15.7 96.9 15.6 64.6
36.9 63.0 29.8 36.2
-0.7 -73.3 -1.1 -59.4
3.9 12.3 6.1 11.1
*or (mm)
Bias = B (Co) or (mm)
oRoot Mean Square Error = RMSE (C ) or (mm)
Relative RMSE = RRMSE (%)Standard Deviation = SD (Co)
•......w Percent of YR MONTHS Absolute Differenceo -
> (1 C or 10 mm)
Largest Negative Difference (Co or mm)
Largest Positive Difference (Co or mm)
Percent of YR~ONTHS direction of changefrom the previous YR MONTH in the CABvalue agrees with the actual (NCDC)values (%) 98 63 96 80
Pearson correlation coefficient betweenCAB and NCDC values 1.00 0.76 1.00 0.88
*Differences in standard deviation between Tables 3 and 4 are due to the use of an unbiased estimator ofsample variance in Table 3 and the Maximum Likelihood estimator of sample variance in Table 4.
There was a significant difference among CRDs (P < 0.0001) for both temperatureand precipitation values in the amount of difference between CAB and NCDCvalues. Table 5 shows the average differences for the forty-four year/monthswith complete data for all CRDs. The values are presented in descending orderand one can see that the ordering of CRDs was not similar for temperature andprecipitation differences nor were consistent patterns of N-S, E-W gradientsseen, nor was there division by biotic region (Stewart 1975).Next, an investigation was made of whether the differences might be related to aconsistent bias. Plots of temperature differences indicated that the techniqueemployed by CAB generally overestimated the NCDC value, and that this overesti-mation was fairly consistent throughout the year (Fig. 2). A simple transforma-tion was made to remove the bias in CAB estimates by subtracting the meandifference (from Table 3) from each CAB temperature value for a particular CRD.This transformed data set was then tested against the NCDC values with thenonparametric ANOVA described above. Significant differences between CAB andNCDC values were found in only 2 CRDs (CRD 20 and CRD 90, see Appendix C). Thenonparametric ANOVA alone was run because the transformation (subtraction ofmean differences) rendered any parametric test for differences in means useless.Thus, it seems that the CAB process for estimating temperature values is biasedin a consistent and simple manner, varying primarily with geographic location(CRD).Plots of precipitation differences indicated that the CAB process generallyunderestimated the NCDC values, but that this underestimation was not of a con-sistent magnitude throughout the year (Fig. 2). The magnitude with which CABprecipitation values underestimated the true value appears to be the greatestduring the warm season months, May through September. The influence of thesemonths on CAB estimate bias may be due to greater spatial variability andamounts of rainfall during the summer. Plots of raw precipitation values (Fig.1) show that high amounts of precipitation occurred during the summer months.Precipitation patterns during this period tend to be discontinuous and patchy(Huff and Shipp 1969). Since the average distance between stations used by CABis much greater than the average distance between stations in the denser networkused by NCDC, the probability of recording patchy precipitation may be lowerwith the stations used by CAB. It is also possible that the particular (non-randomly chosen) stations used by CAB generally have less rainfall than the sur-rounding areas. However, since storm system patterns tend to vary from year toyear, it is doubtful whether this underrepresentation by itself could explainthe underestimation of rainfall over the four summers. The bias in CAB precip-itation estimates appears to vary not only geographically, but temporally and inproportion to the true amount of precipitation also. This conclusion is sup-ported by the effect of different data transformations on the testing procedureoutlined earlier. A simple additive adjustment to remove the bias was made toCAB precipitation estimates by subtracting the mean differences by CRD fromTable 3) from each value, just as was done for temperature estimates. Thetransformed CAB data were then tested against the NCDC data using the nonpara-metric statistic. The simple, additive, transformation produced precipitationvalues significantly different than the NCDC values in all but two CRDs(Appendix C).
14
Table 5. Average monthly value differences (CAB va1ues-NCDC values) for eachCRD, in descending order (n=44)
A. Temperature differences (MSE for 344 d.f. = 0.5779)
CRD Mean
10 1.02520 0.94130 0.79670 0.67390 0.67140 0.58650 0.43460 0.40280 0.121
B. Precipitation differences (MSE for 344 d.£. = 232.6970)CRD Mean
10 -5.44580 -7.19890 -8.31860 -10.32070 -10.84340 -11. 36850 -17.41430 -19.17720 -20.120
15
A second transformation of CAB precipitation estimates was then made as follows:
(i) A simple linear regression of CAB values on NCDC values was conductedfor each CRD. The NCDC values were used as the independent variablein this regression because they were assumed to be the '~rue" precip-itation values. Regression coefficients are given in Appendix E.
(ii) The equations obtained in (i) were then solved to give predictionequations for NCDC values of the form
NCDC =
where bo and b1 are the estimated regression coefficients of interceptand slope, respectively, from the regressions conducted in (i).
(iii) CAB precipitation estimates were then plugged into these predictionequations and a set of transformed precipitation values obtained.Negative estimates were set equal to zero. This process is commonlyreferred to as calibration or inverse regression (Krutchhoff 1967,Halperin 1970, Martinelle 1970). The transformation made in thismanner was basically a proportional transformation, having a greatereffect on CAB values when actual precipitation was high. This secondset of precipitation estimates was also tested against NCDC valueswith the nonparametric ANOVA. The proportional transformation pro-duced precipitation values that were not significantly different fromNCDC values in any CRD (Appendix C).
The preceding analysis of differences in CAB and NCDC weather values is subjectto the same criticism as the original testing procedure; that is, it givesinformation about the central tendency of the distribution of differences butprovides little insight to the reliability of CAB in producing values that areclose to NCDC values. For example, it seems that the error in CAB temperatureestimates does not depend on the magnitude of the temperature value. As aresult, the average deviation of CAB from NCDC values in any given CRD may beexplained by a simple additive effect, as shown by testing the transformed CABdata against the NCDC values. To determine if the bias of CAB temperatureestimates was consistent across all years and months, reliability measures werecalculated from the set of differences between transformed CAB and NCDC data(Table 6). Because the transformation removed the bias in CAB estimates, theroot mean square error and standard deviation for these difference values areequal. Comparing the reliability measures in Table 6 to those calculated foruntransformed CAB temperature data (Table 4), it can be seen that the reliabil-ity of transformed CAB temperature estimates was somewhat better than untrans-formed CAB estimates. The RMSE was smaller for transformed values, and thepercent of observations which were greater than 10 C decreased by about 15% whenaveraged over CRDs.
Measures of reliability for proportionally transformed CAB precipitation esti-mates are presented in Table 7. The bias of these transformed values wasshifted from large negative values (see Table 4) to small positive values orzero. The bias associated with transformed values was not taken to zero in all
16
Table 6. Measures of reliability for transformed CAB temperature estimates
CRDReliabilit Measure 10 20 30 40 50 60 70 80 90
Bias (Co) 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
Root MSE (Co) 0.78 0.68 0.79 1.10 0.78 0.96 1.08 1.14 1.11
Relative RMSE (%) 17.2 15.6 18.7 19.7 16.1 20.0 17.3 19.3 19.2
Percent of YR MONTHS AbsoluteDifference > 1°C 19.1 15.6 15.2 34.8 19.1 17.0 17.0 25.5 22.7
Largest Positive Difference (Co) 1.6 1.5 3.1 3.0 3.3 4.2 5.4 4.3 4.1
Largest Negative Difference (CO) -2.2 -1.3 -1.5 -2.8 -1.2 -2.0 -1.8 -1.6 -1.8I-'....,
Pearson correlationcoefficient betweentransformed CAB and NCDCvalues 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00
Table 7. Measures of reliability for transformed CAB precipitation estimates
CRDReliabilit Measure 10 20 30 40 50 60 70 80 90
Bias (rom) 0.45 0.23 0.00 1.40 1.07 0.00 0.00 0.00 0.63
Root MSE (mm) 20.98 35.71 29.2!f 20.84 26.29 19.21 17.33 9.75 19.02
Relative RHSE (%) 67.0 87.8 70.1 57.2 65.8 47.0 51.0 27.3 48.3
Percent of YR MONTHS AbsoluteDifference > 10 rom 34.0 64.4 65.2 45.6 46.8 55.3 34.0 21. 3 52.3
Largest Positive Difference (mm) 62.0 98.3 84.2 48.8 74.7 56.4 56.5 26.2 44.0
f-" Largest Negative Difference (rom) -91. 7 -82.0 -64.2 -59.0 -74.4 -48.4 -46.3 -23.8 -44.3CP
Pearson ccrrelationcoefficient betweentransformed CAB andNCDC values 0.75 0.69 0.76 0.83 0.79 0.88 0.88 0.96 0.86
CRDs because when the transformation yielded a negative estimate of precipita-tion, that estimate was set equal to zero for the sake of logic. Thus, althoughnot strictly equal for CRDs with positive bias, the RMSE and RRMSE were extreme-ly close in value to the standard deviation and relative standard deviation(these later quantities not reported). In contrast to the situation for theadditive temperature transformation, the RMSE of differences between NCDC andproportionally transformed CAB precipitation values (Table 7) showed an increaseover the RMSE of differences calculated from untransformed CAB values (Table 4).Increases were also seen in relative RMSE, and the percent of observations forwhich the absolute difference was greater than 10 IDID. These increases indicatethat although the proportional transformation decreased the average differencebetween CAB and NCDC precipitation estimates, it provided little informationabout consistent patterns in the CAB process for estimating precipitation; thetransformed CAB values are closer to NCDC values "on the average," but thereliability of such values is even worse than that of untransformed values.
Adjusting CAB Estimates
Although it might be possible to use the transformations made to CAB weatherestimates as adjustments to improve the original values, such adjustment is notwholly appropriate. The transformations were made to investigate the way inwhich CAB values differed from NCDC values and thus utilized the entire data setavailable. Adjustments made in this manner cannot be independently tested andmay not be applicable in the future. Also, the proportional transformation madeto CAB precipitation estimates is clearly undesirable as an adjustment, since itfailed to improve the reliability as compared to untransformed CAB values; suchan adjustment might even yield worse precipitation estimates than the originalvalues. Any adjustments made to improve CAB weather values should be adjust-ments to the CAB process itself, not to the set of data that results from thatprocess. Finally, it should be noted that the ability to make an appropriateadjustment depends in part on a well-defined and consistent method for convert-ing point source data to CD level values. The process used by CAB to derive CDweather values is largely subjective, thus adjustments made to a set of CABvalues may be inappropriate for further use if that subjective process changesin any way.
EFFECT OF ESTIMATED WEATHER VALUES ON CROP YIELD MODEL PREDICTIONS
The focus of this study was the direct comparison of a real-time weather valuesource to the standard NCDC source commonly used in the development of crop-yield prediction models. In the case investigated, statistically significantdifferences were found to exist between the real-time and NCDC data. An obviousquestion is whether differences between data sources have a substantial effecton the performance of predictive models. If not, evaluation of real-time datamay become a moot point.
To determine if prediction models are sensitive to statistical differences inweather input data, barley yield prediction models for North Dakota were runwith both NCDC and CAB weather values as input data. Two models were run, onedeveloped by NOAA, Center for Environmental Assessment Services (CEAS), now AISC(Motha 1980), which has been slightly revised by USDA, Yield Evaluation Section
19
(YES) and a different, more recent model developed by YES. The CEAS model con-tains 2 trend terms and 4 weather-derived variables, while the YES model con-tains the same two trend terms and four different weather-derived variables(Table 8). The CEAS model has been evaluated by several authors (Barnett 1981,LeDuc 1982) and has been found to perform better in years with high yields thanin years with poor yields (LeDuc 1982, p. 42). The YES model has been recentlydeveloped and not been fully evaluated. Using a bootstrap procedure, bothmodels were run at the state level for the years 1979-1982, and the predictedyields compared to actual yields for each year (Table 9). Both models provedsensitive to differences in input weather data and performed better when NCDCvalues were used. In only one instance, the CEAS model in 1980, was a betterestimate of actual yield obtained when CAB data were used. This is, perhaps,not surprising when one considers the circumstances. Barley yield in NorthDakota was extremely low in 1980, and the CEAS model is insensitive to badyears. Thus, provided accurate weather input, one would expect the CEAS modelto overestimate yield in this bad year, and such is the result when NCDC dataare used. But CAB values generally are gross underestimates of precipitationfor the months included in the CEAS model (Table 8), and this tends to pullyield estimates down; CEAS yield estimates using CAB values were much lowerthan the actual yield in 1979, 1981 and 1982. Thus, in 1980, the insensitivityof the CEAS model to poor years was countered by the underestimation of precipi-tation by CAB. Consequently, the two inaccuracies produced a more nearly cor-rect estimate. The YES model appears to be sensitive to poor years as it pro-duced an extremely low estimate of barley yield in 1980 when CAB data were used.In addition, the YES model provided a closer estimate to the observed yield thanthe CEAS model in all years when NCDC values were used (Table 9). Thus, it isimportant to consider the characteristics of a model used in a study of real-time data sources. Despite the interesting results for 1980 and the higheraccuracy of the YES model, statistical differences between CAB and NCDC weathervalues were reflected in the yield predictions of both models, and both per-formed better when given NCDC values as input. This conclusion is furtherstrengthened by measures of yield reliability (Table 10). Bias and Root MeanSquare Error were both substantially greater for models using CAB data than formodels using NCDC data. Although calculated using only four years of data,these measures indicate that the reliability of yield estimates was greater whenNCDC values were used in the models.
CONCLUSIONS OF CAB EVALUATIONThe CAB procedure of deriving CD level temperature and precipitation estimatesuses point-source data that are available on a real-time basis. The time framein which the estimates from CAB are sent to AISC is currently poorly defined anderratic but could be improved to provide CAB estimates in a time frame compat-ible with operational use of crop-yield prediction models. CAB weather esti-mates are not equivalent to NCDC values. The differences between CAB and NCDCvalues can be explained, in part, by systematic bias in the CAB procedure, par-ticularly for temperature estimates. Thus modification of the CAB process maybe possible for temperature estimation, but it would be more difficult forprecipitation estimation. Because CAB weather values, as they are now derived,are not equivalent to NCDC values and because differences between CAB and NCDC
20
Table 8. Description of variables included in state levelCEAS and YES barley models for North Dakota
Model
CEAS
YES
Variable
Trend 1Trend 2
Cumulative Precip.
Precip. SDFN
Precip./PET
Temp. DFN
Trend 1Trend 2
PPS~DEV
AWT DEV
TAHD DEV
MSTR DEV
Description
Trend terms are linear or quadratic functionsused as surrogate variables for technologicaladvances.
Total precipitation from September throughMay.
Departure, squared, of June precipitation fromlong-term normal.
June precipitation divided by potential evapo-transpiration for June.
Departure of July temperature from long-termnormal.
See above.
Preseason precipitation from September toplanting, deviation from average, squared.
Average winter temperature over December,January, and February, deviation from average.
Temperature around heading (1 week before and2 after), deviation from average
Moisture stress, deviation from average.
21
Table 9. Comparison of CEAS and YES estimated barley yieldswith actual yieldst 1979-1982 (bushels/acre)
CEAS YESYear Actual NCDC CAB NCDC CAB
1979 46.0 40.0 32.9 45.0 42.3
1980 32.0 39.2 26.2 32.5 9.9
1981 48.0 45.4 35.2 47.2 28.1
1982 53.0 46.9 43.3 51.8 46.4
22
Table 10. Measures of yield reliability for CEAS and YES barleymodels using NCDC and CAB weather values over 4 years
CEAS YESMeasure of Reliability NCDC CAB NCDC CAB
Bias (Bu/acre) -1.88 -10.35 -0.63 -13.082 32.95 115.78 0.83 235.32MSE (Bu/acre)
RMSE (Bu/acre) 5.74 10.76 0.91 15.34
Largest IRDI (%) 23 -28 ±2 -69
Pearson r 0.82 0.94 1.00 0.90
23
precipitation values are complex and inconsistent, it would be desirable to finda more objective method to obtain CD level weather values from real-time point-source data.
DISCUSSION OF EVALUATION PROCEDUREA number of complicating factors influence the objective evaluation of real-timeweather data. The characteristics of the model or models in which the data willbe used may mask effects of different data sources. The spatial scale on whichdata are evaluated may dictate the types of models for which the evaluation ispertinent (Strand 1981). Because of these factors, it is desirable to evaluatereal-time data sources directly against a standard set of values before anattempt is made to use them in crop-yield models. An evaluation of predictedyields obtained from using real-time data as input variables is primarily anevaluation of the prediction model, not the data source. More work is needed onthe effects of errors in input weather variables on the reliability and accuracyof yield predictions. Understanding the way in which different sources of errorinfluence the final product of a prediction equation would help eliminate theneed for re-evaluation of data sources as new models are developed and broughtinto operational use.Direct comparison of real-time data with standard values has been incorporatedinto a generalized procedure for evaluation of real-time data sources (Fig. 3).This procedure considers the three evaluation criteria listed in this documentand is designed to provide a concrete conclusion without unnecessarily con-straining statistical methodology. The authors hope it will prove useful toother investigators evaluating real-time weather data for use in crop-yieldprediction models.
In the direct comparison of real-time data to a set of standard values, statis-tical procedures which test for differences in average values or central ten-dency are not sufficient to infer equivalence. Real-time weather values shouldbe good estimates of true values on the average and distributed in the same wayas true values, but they must also be consistent estimates at all points intime. Indicators of reliability can be used to measure this consistency. A setof measures of reliability is necessary because no one measure alone adequatelydescribes the behavior of real-time estimates. The correlation coefficient, inparticular, is not a sensitive indicator of reliability. A set of real-timedata may show a strong linear relationship with true values, yet have a highpercentage of values for which the relative difference is greater than acritical limit (see Table 7) or have a high bias (see Table 4). These types ofdifferences can affect the prediction of crop yield models, particularly modelswhich are sensitive to changes in weather values.
24
Figure 3. Flow chart for suggested paradigm in evaluation of real-time data sources.
®
Compare predictedyields from oper-at ional IIIOdels.
lmpleoent use ofreal-time datasource.
Search for bettermethod of obtainingCD level values.
25
Modify methodto cor:rect
inaccuracies
Reject DataSource
Barnett, T. L., 1981.Dakota and Minnesota.YMD-1-4-1(81-12.1).
LITERATURE CITED
Evaluation of the CEAS model ,-barley yie.Lds in NorthAgRISTARS Yield Model Develol-'u.~ntProject, Document
Box, G. E. P., 1953. Non-normality and tests on variances. Biometrika 40:318-335.
Cochrane, M. A., Jr., 1981. Soi 1 moisture and Agromet mode Is. USAF Env iron-mental Technical Applications Center, Technical Note, Report. No. USDAETAC/TN-81/001.
Geary, R. C., 1947. Testing for normality. Biometrika 34:209-242.
Halperin, M., 1970. On inverse estimation in linear regression. Technometrics12:727-736.
Huff, A. A., and W. L. Shipp, 1969. Spatial correlations of storm, monthly andseasonal precipitation. ~ Applied Meteorology 8:542-550.
Krutchkoff, R. G., 1967. Classical and inverse regression methods of calibra-tion. Technometrics 9:425-439.
LeDuc, S. K., 1981. Spring wheat models for North Dakota and Minnesota.NOAA/CEAS/Models Branch, Columbia, Missouri.
LeDuc, S. K., 1982. Comparison of the CEAS and Williams-type barley yieldmodels for North Dakota and Minnesota. AgRISTARS Yield Model DevelopmentProject, Document YMD-1-4-2(82-1.2).
Martinelle, S., 1970. On the choice of regression in linear calibration.Comments on a paper by R. G. Krutchkoff. Technometrics 12:157-161.
Motha, R. P. 1980. Barley models for North Dakota and Minnesota. NOAA/CEAS/Models Branch, Columbia, Missouri.
Perry, C. R., J. L. Rogers, and Joe T. Ritchie, 1982. A comparison of measuredand estimated meteorological data for use in crop growth modeling. Proceedingsof Business and Economics Statistics Section of the American Statistical Associ-ation, pp. 77-82.
Sakamoto, C. M., N. D. Strommen, and S. K. LeDuc, 1977. Synoptic vs.cooperative climate data as inputs for wheat yield model estimates. In 13thAmer. Meteor. Soc. Conference on Agriculture and Forest Meteorology. PurdueUniversity, April 4-6, 1977.
Sebaugh, J. L., 1981. Evaluation of the CEAS trend and monthly weather datamodels for spring wheat yields in North Dakota and Minnesota. USDA/SRS/SRDStaff Report No. AGES8112l4.
Stewart, R. E., 1975. Breeding birds of North Dakota. Tri-Coll ege Center forEnvironmental Studies. Fargo, North Dakota. 295 pp.
26
Strand, B. W., 1981. Spatial scale of crop-yield models: A review of the rela-tionship between scale of models and accuracy. USDA/ESS/SRS Staff Report No.AGES810320.Tukey, J. W., 1948. Some elementary problems of importance to small samplepractice. Human Biology 20:205-214.
Willmott, C. J. and D. E. Wicks, 1980. An empirical method for the spatialinterpolation of monthly precipitation within California. Physical Geography1:59-73.
Wilson, W. W. and J. L. Sebaugh, 1981. Established criteria and selectedmethods for evaluating crop yield models in the AgRISTARS program. 1981 Proc.American Statistical Association: 24-31.
27
APPENDIX A. REAL-TIME WEATHER DATA SOURCES IN NORTH DAKOTA
Potential sources of weather values for North Dakota were identified throughcontact with personnel from the U. S. Department of Agriculture (USDA) YieldEvaluation Section (YES), the U. S. National Oceanic and Atmospheric Adminis-tration (NOAA) Assessment and Information Services Center (AISC) and NationalWeather Service (NWS), and the North Dakota Crop and Livestock ReportingService.
Daily weather data for the U. S. (which may be processed to provide monthlyvalues) are available from several types of systems, including (1) cooperativesurface networks, (2) synoptic surface networks, and (3) satellite observations.The distinction between cooperative and synoptic networks may appear vaguebecause some weather stations provide information to both systems. However, thedivision indicates different approaches to the collection and compilation ofweather station data. Cooperative networks are characterized by the collectionof basic weather data from as many weather stations as possible in a given areaand make use of many varieties of weather stations. Synoptic networks attemptto provide an efficient summary of large scale weather patterns from stationsdeemed to prov ide re liab Ie and time ly data. Synoptic networks are general Iysparse, making use of principal weather stations, and rely on automated report-ing techniques.
Types of Collection Systems
Four sources of real-time weather data were identified from the three types ofcollection systems listed above (Table AI):
1. Cooperative Networks
Cooperative weather observers in North Dakota record temperature and precipita-tion values daily on forms provided by the NWS office in Bismarck, N.D.Although the number of stations reporting varies, records received by YES indi-cate that about 100 cooperative stations report in any given month. The NWSforms are designed to contain one month of data and are generally received inBismarck between the 4th and 10th days of the following month (Car" e personalcommunication). Data are collated in Bismarck and sent to NCDC in Asheville,North Carolina. In addition, seven "major" stations send monthl' summariesdirectly to Asheville. These data are the basis for the monthly averages pub-lished by NCDC. The strong point of this system is that data are available froma large number of stations throughout the state. However, because all data arenot available until the 10th day of the following month, current procedures arenot sufficient for real-time models.
Another cooperative program in North Dakota provides weather values on a weeklybasis. From April to October a subset of the NWS cooperators report weeklytemperature and precipitation values to the NWS office in Bismarck by telephoneeach Monday morning. These data are then transmitted to the offices of theNorth Dakota Crop and Livestock Reporting Service in Fargo, N. D., also bytelephone, where they are compiled and published in the weekly Crop-WeatherReport. This program involves about 30 reporting stations, and the Crop-WeatherReport is released on the first work day of each week.
28
Table A1. Potential real-time weather data sourcesfor North Dakota
Type of ReportingSystem
1. Cooperative Network
2. Synoptic Network
3. Satellite Data
Source
NWS monthly forms
N.D. Crop-WeatherWeekly Report
WMO DMONTH
USAF Agromet
Type of Data
Point Source,100 stations
Point Source,30 stations
Point Source,13 stations
Estimated PointSource on grid
NWS: National Weather ServiceWMO: World Meteorological OrganizationUSAF: United States Air Force
29
2. Synoptic NetworkThe Global Telecommunications System (GTS), administered by the WMO, is a largescale network that reports values from first order and automated weather sta-tions throughout the world. These values are processed weekly by the ClimateAnalysis Center (CAC) of the NWS, and a summary produced weekly in the form of acomputer printout known as DMONDAY. A monthly summary known as DMONTH is alsoproduced by CAC and contains the following weather data for each reportingstation: number of days precipitation readings taken; percent of observationsfor which precipitation was recorded; total precipitation; normal precipitation;departure of recorded precipitation from normal; percent of normal precipitationrecorded; number of years used to calculate normal precipitation; number of daystemperature readings taken; mean temperature; normal mean temperature; departureof recorded mean temperature from normal; number of years used to calculatenormal mean temperature; highest daily maximum temperature for month; date whenhighest daily maximum occurred; number of days maximum temperature recorded;lowest daily minimum temperature for month; date when lowest daily minimumoccurred; and number of days minimum temperature recorded. Precipitation valuesare reported in millimeters and temperature values in degrees Celsius. TheDMONTH summary currently reports values for thirteen stations in North Dakota; 3National Weather Service (NWS) stations, 4 Federal Aviation Administration (FAA)stations, 2 stations from the Continental U. S. Meteorological Data System(COMEDS), U. S. Air Force, 3 Automatic or Remote Automatic MeteorologicalObserving Stations (AMOS/RAMOS), and 1 municipal airport station (Table A2).The DMONTH summary is available from CAC within several days of the end of amonth of recorded data, making it a timely source of data for real-time models.Because the GTS is a large scale network, however, the number of stationsreported for North Dakota is small, and some type of data manipulation is neces-sary to derive values suitable for us with models requiring weather values on aCD or CRD scale.
3. Satellite DataThe U. S. Air Force has developed a procedure known as Agromet to estimatemeteorological values from data provided by polar orbiting satellites and theWMO network. The variables used by Agromet are reflectance, radiance, andtemperature from satellites, and temperature and precipitation from surfacereports. Agromet produces point estimates of maximum and minimum temperatures,precipitation, solar radiation, and evaporation with an Air Force 25-NM gridpoint spacing at 60° N. Agromet is run every three hours at Global WeatherCentral (GWC), Offutt Air Force Base in Nebraska. Daily summaries are preparedand compiled as weekly data sets each Monday. Weekly data sets are currentlyflown to the USDA, Foreign Agricultural Service (FAS), Foreign Crop ConditionAssessment Division (FCAD) in Houston, Texas where they are p.ocessed and madeavailable for use within 60 hours (Perry and Rogers, 1982). A complete descrip-tion of Agromet is provided in USAF ETAC/TN-81-001 (Cochrane, Jr., 1981). IfAgromet data prove appropriate for use with real-time models, data collectioncosts could be reduced (Perry and Rogers, 1982). Unfortunately, Agromet datafor the U. S. have been available only since June of 1981.
30
Table A2. North Dakota Stations Reported on DMONTH Program
Station CallNo. No. Latitude Longitude CD/CRD Location Type
72764 BIS 46.8 100 .8 80 Bismarck NWS72753 FAR 46.9 96.8 60 Fargo NWS72768 ISN 48.2 103.6 10 Williston NWS
MOT 48.3 101.3 20 Minot FAAGFK 47.9 97.2 30 Grand Forks FAADIK 46.8 102.8 70 Dickinson FAAJMS 46.9 98.7 50 James town FAAMIB 48.4 101.3 20 Mino t COMEDSRDR 47.9 97.4 30 Grand Forks COMEDS
72758 DVL 48.1 98.9 30 Devil's Lake MUNI ARPT72757 Pll 48.1 98.9 30 Devi1's Lake AIDS/RAMOS
P24 47.8 101.8 40 Roseg1en AMOS/RAMOSP67 46.1 97.1 90 Lidgerwood AMOS/RAMOS
31
Conversion of Point Source Data !£CD Level Values
Since crop-yield models use weather variables calculated at the CD level, pointsource data must be manipulated in some manner. Three methods currently showpotential for converting raw data to CD level estimates (Table A3):
1. Simp le Averages.
Monthly weather values are simply averaged across all reporting stations for aCD by NCDC. Thus taking simple averages of real-time station data is an obviouscandidate when searching for operational methods of obtaining timely CD levelvalues. The potential of this method lies in the possibility that a sample ofcooperative network stations could (1) provide sufficient information for accu-rate estimates and (2) be obtained on a real-time basis.
2. Subjective Estimation.
Two sources of real-time weather data make use of subjective estimation proce-dures to derive temperature and precipitation values for CDs. These sources are(1) the USDA World Agricultural Outlook Board (WAOB) and (2) NOAA, ClimaticImpact Assessment Division, Climate Assessment Branch (CAB), both located inWashington, D.C. Both receive station data from the WMO GTS (DMONTH program runin Suitland, Maryland), plot temperature and precipitation values on weathermaps, draw isotherms and isohyets (lines connecting points of equal temperatureand precipitation, respectively), and make subjective estimates of averagetemperature and local precipitation for each CD. Because these are the sameweather variables reported by NCDC, and data are available on a CD level, littlemanipulation of the data would be needed to use these estimates in formulae forderiving model weather variables. In addition, the cost of obtaining weatherdata from these sources is low since the estimates are derived for purposesother than USDA crop yield modeling. Weather estimates from CAB have beenobtained by YES since January 1979 through NOAA, Assessment and InformationServices Center (AISC) personnel. Estimates from WAOB have been received sinceFebruary 1982. CAB and WAOB data were received by AISC an average of 26 and 7days, respectively, past the end of the month for which data were pertinent (seeTable A4). The procedure for transfer of values from CAB and WAOB to AISC hasbeen poorly defined on both ends. The most rapid turn around time (5 days)indicates that figures could be reported in a time frame consistent with monthlyuse of crop yield models, but the reporting method is currently too erratic forsuch use.
3. Objective Analysis or Spatial Interpolation
A number of techniques have been developed in the past 15 years to estimatemeteorological values for an area, given a minimum of input data. Often asurface is fitted to a small number of points, and estimates of weather vari-ables produced for many locations over the same geographic area. An appealingquality of these techniques is that they remove the subjective bias inherent inWAOB and CAB procedures, and may be useful with data provided by the DMONTH orAgromet programs. Whether objective methods can provide more accurate informa-tion than subjective procedures or data obtained from a sample of cooperativestations, however, remains open to question.
32
Table A3. Potential methods for conversion of real-timestation data to CD level values
Type of Method Source
l. Simple average Sample of coop.network stations
2. Subjective estimation USDA WAOB
NOAA CAB
3. Spatial interpolation Any of a variety oftechniques
Data Form
CD level avg.
CD level avg.
CD level avg.
Yields estimatedpoint source datawhich may beaveraged, somemethods give areaavg. directly
33
Table A4. Dates CAB and WAOB weather valueestimates were received at YES, 1981-1982
Source No. days betweenof Year and month last day of month
Estimate of data Date received at YES and date receivedCAB June 1981 6 August 1981 37CAB July 1981 14 August 1981 14CAB August 1981 26 September 1981 26CAB September 1981 missing recordCAB October 1981 missing recordCAB November 1981 missing recordCAB December 1981 5 March 1982 63CAB January 1982 missing recordCAB February 1982 24 March 1982 24CAB March 1982 28 April 1982 28CAB April 1982 6 May 1982 6CAB May 1982 16 June 1982 16CAB June 1982 27 July 1982 27CAB July 1982 25 August 1982 25CAB August 1982 21 September 1982 21CAB September 1982 3 November 1982 34CAB October 1982 22 November 1982 22CAB November 1982 3 January 1983 34CAB December 1982 26 January 1983 26
WAOB February 1982 5 March 1982 5
WAOB March 1982 6 April 1982 6WAOB April 1982 6 May 1982 6WAOB May 1982 9 June 1982 9WAOB June 1982 6 July 1982 6WAOB July 1982 5 August 1982 5WAOB August 1982 9 September 1982 9WAOB September 1982 6 October 1982
34
APPENDIX B, Part 1: Corresponds to Figure 1 in text.
Cllto=10
III+III C+III+,---+++++++++++++++++++++++++++++++++++++++++++++++---
7777777117788888888~~88888888888888888888888A8A9999999999900000000WOOOl1ll11111111222222222?22OOOOOOOOll100000000WII10000000001110000000001112345618901212345678~0121234567890121234S6789012
YEAR AND MONTH (YR_MONTH)
-30
10 +
-20
CCN
CN
C
C
C
CCNC
C
CN
C
CC Ce CC NN NC
e C e N C C NC N C eN C N C
N N C NC e C NN N C N C C
C '" N C NC C NCC
C
o +
II20 +III
TEMPERATU -10RE
CIttI)=lOII
150 I+III120 +I NII
PRECrpT 90 + C CT I ,\jA ITIC CI 60 + C C NN Co I NN N NC '"N I C C '" CI C N t) N N30 + CC C N N CC e CI NN N C NC NN C CI N NC~ C CC C CI CCN N NC C C CCN C C Co + CN ec N C N
---+++++++++++++++++++++++++++++++++++++++++++++++---771777171778888888d~H8888888888888888888888AR889999999999900000000WOOOll1111111111222222222?220000000011100000000wIIIOOOOOOOOOIIIOOOOOOOOOlll234567890121234S678~012123456189012123456789012
YEAR AND MONTH (YR_MONTH)
35
Appendix B, Part 1. Continued
CI<I)=20
o +
II!
If) +
II20 +
TEMPERATlJ -10 +R 'E
C
C
C~J
CNC C
C
CC NN
I\J
N eeN
C
eC C~
l.
cC
eCCN
C
C
COJe
N CC N
CN
C~)
CN e"I
ccNC Ce "I
rC NN
c e~I C
-20 + CCN
-10 •.I---+++++++++++++++++++++++++++++++++++++++++++++++---
7777777777788888888~H8888A8888888888888888e-q88999999Q999900000000wOOOlllllllllll1222222222?22OOOOOOOOlllOOOOOOOO~lllOOOOOOOOOlllOOOOOOOOOlll214S678901212345678~o121234567890121234S6789012YEAR AND MONTH (YR_MONTH)
ehiD=20N
f\JN ~IN N NeeCe e C NC
N NNN N C C C N
C r "J N "Ie N CN Nee NCC C NNN NNN C eCNNC e e e ceo + eee eee ee e ec c---+++++++++++++++++++++++++++++++++++++++++++++++---
7777777777788888888dH888888888888888888888888889999999999900000000~(lOOllllll]11111222222222?22OOOOOOOOlllOOOOOOOOijlll00000000011loooooonOOlll234S678901212345678~()121234567890121234S6789n12
YFAR AND ~ONTH (YR_MONTH)
ISo ..DqE 120 ..C IT Ip II 90 ..T I
" IT II ~o +0~J
10 ..
C
N
C
NN ~J
Ncet-'
~'J
36
Appendix B, Part 1. Continued
CII/0=40
I) .•.
-20 +
20 +
If) +II,
CNC
C
c
CNCC
C
cCN
CN
C CN
cC
ecN
CeN CN
CN
c
eNCC N~.J b
N
eeCI\jN
cC
c
',.J CN N
CN
C
C
NCNC CC
C
c
C
TfMPEQATU -10 +R IF: I
I
-10 +I---+++++++++++++++++++++++++++++++++++++++++++++++---
77777777777g8888888~H888888888888888888888888889999999999900000000WOOOlllllllllll1222222222?22OOOOOOOOlllOOOOOOOOWIIIOOOOOOOOOIIIOOOOOOOOOlll2145678901212345678~012123456789012123456789n12YEAR AND MONTrl (YR_MONTH)
CII/0=40
C30 •. N N C N C N CC C CN N e C C CN C NN C N C NN N NC NCCCCN ec~c C C e C N0+ C N CC C
---+++++++++++++++++++++++++++++++++++++++++++++++---7777777777788888888~~888888888888888888888888889999999999900000000wOOOlllllllllll1222222222?220000000011100000000wIIIOOOOOOOOOIIIOOOOOOOOOlll234567~90121234567~~1I12123456789012123456789012
YEAR AND ~nNTH (YR_MONTH)
ISO +
ppE 120 .•.C IT Ip II go +T IA IT II ~O ..ntJ
t\1
NN
NN CN ~.J C
I\J C NCNNC C NN C
C C N e CCN N e e
37
Appendix B, Part 1. Continued
CIl(I)=50
II C N C C2() + NC C CNC NC
C N C C I\JC I\J (,; C C
N C C CT 10 + C CE I C C C C"'1 I C CP I C C C '"E 0 + C I\J r-tR I NA I N N C C CCT I C C N NU -10 + C CCR I N C CE I r
I-20 + C CI I\JiI-30 +I---+++++++++++++++++++++++++++++++++++++++++++++++---
77777777777A8888888~M8888A8888A8888888888888A889999999999900000000wOOOllllllll1111222222222?22OOOOOOOOl1100000000WIII0000000001110000000001112345678901212345618'li012123456789012123456789n12
YFAR AND ~ONTH (YR_MONTHl
C••••')=50
N "CJ CC C N C C
I\J N (,;N NN C C N C CNNN NN N~CCC C C C No + cc ceee eee c---+++++++++++++++++++++++++++++++++++++++++++++++---
7777777777788888888oM8888888888888888888888BA88999999999g900000000woOOlll111111111222222222?22OOOOOOOOlllOOOOOOOO~111000000000lI10000000001112345678901212345678~o12123456789012123456789n12
YFAR AND ~ONTH (YR_MONTH)
N150 +
Ip IR I~ 120 +C IJ Ip IJ gO +TA.TI 60 +0 II\J I
I10 +
N
NI\Je C
"I CC
I\JC c '"N
'"N
I\J NC
I\JC '"C
C C CN C C CI\J
I\J C
38
Appendix B, Part 1. Continued
CIl(II=60
II20 +
T In +
fMPf (')+Q
ATI) -10 +o IE I
I-~o + c
c
Ct\1
C
CCC t.J
C
Po.J
N CC
NN CN
CC
NCCC
rJ LN
C
c c
CCc
C
rCNC'" CC N
C
c
C C
C
C
C
('
~iC
CC CI\J ~I
('
CCt.J
-10 +I---+++++++++++++++++++++++++++++++++++++++++++++++---
77777777777q8888888~~8888A88AAq8888888889A88A889999999999900000000WIlOOll1111111111222222222?2?0000000011100000000w1110000000001110000000001112145678901212345~78~11121234567R9012123456789n12YEAR AND MONTH (YR_MONTH)
('
N
C Po. 1 ••••JI\j NN·C N CC
I\J N e CCI\jCo + ce ceee C---+++++++++++++++++++++++++++++++++++++++++++++++---
77777777777q88R8R88~H8888P88AA8888888q88888qA889999999999900000000w000111111111111222222?22?220000000011100000000Wl11000000000111000000000111234567890121234~678~o12123456789012123456789n12
YEAR AND MONTH (YR_MONTH)
ISO +
pqE 120 +C IT Ip II 90 + 1\1TA ".JT CT 60 + C NC0 NPo.J CC NNN N10 + C I\j
NNC
"J CC NNC I\JI\JC CC
N C "ICbC C CN
CC C C
'"
C CC CN I\J
C
39
Appendix B, Part 1. Continued
elo(l)=80
eccN
(lO +
T If) +J='Mp
F n +p0.TlJ -If) +
R~
-20 +
-30 +
e
r
NC
c
cce '\JC
c
N COJ
er
c
NNCeC N1\1 C
rc r
CC
CC CN e
c
e
C
c
c
cr
ce
I---+++++++++++++++++++++++++++++++++++++++++++++++---
7777777777788888888bH88888888888888888888888A889999999999900000000WOOOlllllllllll1222222222222OOOOOOOOlllOOOnOOOOWIIIOOOOOOOOOIIIOOOOOOOOOlll234567A90121234S678~1)121234567B9012123456789n12
YEAR AND MnNT~ (YR_MONTH)
C••.0=80
C NCCC "J "J
C C C N b NCN cec e c c
f) + OIN C C C---+++++++++++++++++++++++++++++++++++++++++++++++---
77777777777R8888R88b~A8888888888888888888888A889999999999900000000wOOOlllllllllll1222222222222OOOOOOOOlllOOOOnOOOUIIIOOOOOOnoOlllOOOOOOOOOlll2345678901212345678~012123456789012123456789n12
YEAR AND ~ONTH (YR_MONTH)
IS!) +
ppF 120 +CrDr 91) +TATI 60 +0 I
'"10 +
c
c
N
e,..v
N
e "" NNCC
CNC
c eCN Ne e
CNN
eN CCC
NC cC CN
40
Appendix B, Part 1. Continued
CIo/O=90
Ne
20 +
T 10 +E IM IP If 0 +R IA IT Iu -10 +Rf
-20 +e
e1\1
C
CeC eeN
"J N Ce
C
c
eN eC
11.1 L
C
C
Cec
I\l
e
CC
eee"J
N e
C
C
eC CN N
e
C
cC
ec e
e
ee
·30 +I---+++++++++++++++++++++++++++++++++++++++++++++++---
77777777777g8888888dH8888A8A8AA888888888888g~889999999999900000000wOOOlll1l111l11122222222222200000000111nOOooOOO~111000000000111000000000111234567890121234,678~0121234567g9012123456789n12
YEAR AND MONTH (YR_MONTH)
ClttO=90
PRFeIp CIT N NA eNTr 60 + N ef) I CCN INN NI C N N C NC C N C
30 + N C NC eI C N L eN C Nec N CI '" N NN N C CC N CINN cce CNe c eo + cc e eNC---+++++++++++++++++++++++++++++++++++++++++++++++---
777777777778888888~OH88888888888888888888888R889999999999900000000~oOOI11111111111222222222?22OOOOOOOOll100000000~111000000000l11000000000l112345678901212345678~012123456789012123456789n12
YEAR AND MONTH (yq_MONTH)
III150 .•.III120 .•.III90 .•.
NN
N C
CN
C Ce N eC"JNe N
41
APPENDIX B, Part 2: Corresponds to Figure 2 in text.
C~U=lO
T~,~D
•ntFFF~FI\j
CEc;
;::>Q
FCTp
•nrF.-FPf'.1CFc;
IIS.O +IIII2.5 + AI It A AA A AI A A A AA A AA AAI A A A ~AA A AA A A AAI A A AA A A0.0 +-------A--------A-AA~--------------------------A----I .A AI 4II-?5 +IIII
-"1.0 +I---+++++++++++++++++++++++++++++++++++++++++++++++---
777777777778d8RAd8~H88b8A888A8888888888B8HAA8R8q9999~9~99qOOOOOOO~OOOOlllllllllll122222222?22~nOOOOOOOlllOOOOOOO~OlllOOOOOOOOOlllOOOOOOOonlll~34S67A9n)?1234567~~01212345A7890121234561Bqn12vEAR AND ~ONTH (YP_~ONTH)
C~iJ=lO
II A
25 + AII A A r,I t.. A A
a +------A------A-------------AA--A-A--AAAAAA-----AA---I A A A A A A AA AAI A A A A AI A A A A
-?5 + A AIII
-50 +III-75 +III r:.-In!) +I---+++++++++++++++++++++++++++++++++++++++++++++++---
77777777777888ARdd~HA888R88888888888R88d888P888q99~999999qOOOOOOO~uOOOlllllllllll12222222~?222nOOOOOOOlllOOO()OOO~ulllO()OOOOOOOlllOOoooooonlll?34~67g9nl?1234S67"~01212345A789012123456789012yEAR A~O ~ONTH (YR~~ONTH)
42
Appendix B, Part 2. Continued
p ;;>5qfCI ap
r) -25IFFF -soRF~I
C -7SFc;
-100
•nTFFFRF~JCFC;
IS.O +IIII
2.5 + AI A A AAAAI A A AAA A AI A A A A A AA A AAA AAI A A AA A A A A A A
0.0 +--4-A---A-----------~A-A----A-----------------------I AIII-2.5 +IIII-S.O +I---+++++++++++++++++++++++++++++++++++++++++++++++---
77777777777888AA8d~H8888A888A88888888888888A888999999999990000000ij0000111111111111222222222222OOOOOOOOlllOOOOOOOijOll1000000000111000000000111?3456789nl?1234S67~~0121?34SA789012123456789nI2yEAR ANO ""'0NTH (YR_~lONTH)
CI(L)=20II+ AI AI aI+--A----------A--------------------AAA-AA-------A----I A AAAA A AAI 1\ A A A AI A A A A+ AI A A A AI A A fJ.I A A+I A AII AA+I ~II A+I---+++++++++++++++++++++++++++++++++++++++++++++++---
77777777777888A8dddH8888888888888888A888A88A8889999999999900000001ll000011111111111122222222?222OOOOOOOOl110000000wol11000000000111000000000111?3456789()121234S67d~OI21c34S67890121234567A9012
YEAR AND ~ONTH (YR_~ONTH)
43
Appendix B, Part 2. Continued
CIt(()=40II+tIII+ AI A A A A AI A A AA A A A AAA AI A A AA AA A A+---------A----------~--------------A----A-A--AA-----I A A A A A
A A A AII+II
I
TE"'1p
•otFFEqF.:~CES
c:;.o
2.5
0.0
-2.5
A
4A
-5.0 +I---+++++++++++++++++++++++++++++++++++++++++++++++---
777777777778888888bH88888888888888888888888A8A8999999999990000000w0000111111111111222222222222000000001110000000WOII1000000000111000000000111234567890121234567b~01212345678901212345678q012
YEAR AND ~ONTH (YR_~ONTH)
C""D=40
+II AI A A+---------------A----~--AAA-A-----A-AAA-A-----A------I A AA AA AA A A A A A AI A A A AAI A AA+III+I
I+II
p 25RECI 0p
•0 -25IFFE -50RF:t-,jC -75ES
AA
A
AA
AA
A
A
-100 +I---+++++++++++++++++++++++++++++++++++++++++++++++---
777777777778888888~M88888888888888888888888A8A8999999999990000000w0000111111111111222222222222n00000001110000000WOII10000000001110000000001112345678Y0121234567~~0121234S6789012123456789012
yEAR AND ~ONTH (YR~~ONTH)
44
Appendix B. Part 2. Continued
C"'I)= I;) 0IIS.O +III AI
2.5 +•[)
rFFF.qE'JCEc:;
0.0
-2.5
-s.o
p 25pFCr 0p
•n -25TFFF -SOREfl.)C -7SFc:;
-100
I A A AI AA A A A A A A A A A AI A A AA A A A A A+--AA---------A----A-~-A---------A-AA------A---AA----11\ A A A A AI A AII+I1I1+I---+++++++++++++++++++++++++++++++++++++++++++++++---777777777778888888bM88~88888~888888888888888888
999999999990000000WOOOOll1l1111l1l1222222222222OOOOOOOOlllOOOOOOOW0111000000000II1000000000111?34567890121234567~~0121?345~7890121?345~78Q012vEAP AND ~ONTH (YR_MONTH)
CIolI)=50II+I1I A A A A+----------------A-------A----------AA-AAA----A-A----I A AA AAA AA AI A AA A AI A A A A+ A AI A A AI A A A AI+II AI A+ A1I AI+I---+++++++++++++++++++++++++++++++++++++++++++++++---
777777777778BA8R8BbH8888A8888A88A8888888888A8889999999999qOOOOOOO~oOOOlllllllllll1222222222222000000001110000000wOI11000000000111000000000111~34567890121234S67~~0121234S~78901212345~78QOI2
yEAR AND ~ONTH (YR~MONTH)
45
Appendix B, Part 2. Continued
C~I)=60
TE'-1P
•r)TFFEp
Ft\JCEc:;
5.0
2.5
0.0
+I AIII A+II A AI A A A AAI A A A A A A AA AA AAAA A AA A A+--_A----A-----------A-A----A--------A-A---A---------I A A AA AI AI A AI+IIII-5.0 +I---++++++++++++++++++.++++++++++++++++++++++++++++---
7111711111188RB888oH8888A888RR8888888888888P888999999999990000000wUOOOlllllllllll1222222222222OOOOOOOOlllOOOOOOOWOIIIOOOOOOOOOIIIOOOOOOOOOlll?34567B90121234567d~0121234S678901212345618Q012vEAD AND ~ONTH (YR_MONTH)
(1011)=60
D ?5RE(I 0D
•[) -25IFFE -50RFNC -15ES
-100
II+II f1..A A A A t\A+-- __ -- ---A----AAAAAAA-A __ ---AA-AA----AA--A---I A A AAI A AA A fJ.I+ A A A AI A AI A A AI+ A AI AII+III+I---+++++++++++++++++++++++++++++++++++++++++++++++---
777711171718888888dH88888888888888888888888A888999999999990000000WOOOOlllllllllll1222222222222OOOOOOOOlllOOOOOOOWOIIIOOOOOOOOOIIIOOOOOOOOOlll?34567B901?1234567~~0121234S6789012123456789012vEAR AND ~ONTH (YR_MONTH)
46
Appendix B, Part 2. Continued
C~[)=80
+IIII+II A AA A AI A A AAA A A A A+-------------------------A-AA--AAA---AAA-------A----I AAA A A A A A AAAAAI AA A AI+IIII+I---+++++++++++++++++++++++++++++++++++++++++++++++---
777777777778888888MH8888A88888888888888888AAAA89999999999QOOOOOOO~oOOOlllllllllll122222222?222OOOOOOOOlllOOOOOOOWDlllOOOOOOOOOlllOOOOOOOOOlll234567890121234567u~012123456189n12123456789012yEAR AND ~ONTH (YR_~ONTH)
TEMP•DIFFEREt'IJCES
5.0
0.0
-2.5
-5.0
A
A
p 25RECT 0p•D -25TFFE -50RENC -75tS
-100
CIt([)=80
II+III A. A A+---A-A-AA-A--A-A---A---AAA-A-----AAAA-A-----A--AA---I A A A A A A A AI A A A AI A A A+ hI A A A A AII+tIt+III+I---+++++++++++++++++++++++++++++++++++++++++++++++---
717777777178888A88MH8888A888888888888888888A8889999999999QOOOOOOOWOOOOlllllllllll122222222?222OOOOOOOOlllOOOOOOO~OlllOOOOOOOOOlllOOOOOOOonlll234567890121234567o~01212345618901212345678Q012vEAR A.ND ~ON1H (YR_~ONTH)
47
Appendix B, Part 2. Continued
ChiD=9U
TF.:~p
I)IFFF'-?F1\1
CFS
5.0
0.0
-2.5
-5.0
II+ AI AIII+I A AI A AAI A AA A AAAA AA AI A A A A AAA A A AAA+--A---A-A----------A~--------A---------AA-----------I A AAft.I A A AII+IIII+
---+++++++++++++++++++++++++++++++++++++++++++++++---777777777778888R88~M8888A8888A88888~888888AA888999999999990000000WUOOOlllllllllll122222222?222OOOOOOOOlllOOOOOOO~Ul1100000000011100000000nlll?345678901?1234567d~012123456789012123456789012
yEAR AND MONTH (YR_~O~TH)
-50 + AI AII-75 +III-100 +
poF
CTp
•nTFFF:RF~.JCF
II?5 + Af t,I AI A A A A AAA A
o +----------------A------AAAA-----A---A---A----A-A----I A A AA A AI AA AA AA AI A A
-25 + AI A A A AIt A
I---+++++++++++++++++++++++++++++++++++++++++++++++---
7777777777788AAA8ddM88888888A88888888888888A888999999999990000000~oOOOlllllllllll122222222?222OOOOOOOOlllOOOOOOOijOlllOOOOOOOOOlllOOOOOOOOOlll~345678901?1234567~~0121234567~901212)45678q012yEAQ AND ~ONTH (YP_~O~TH)
48
APPENDIX B, Part 3: Corresponds to Table 4 in text.
CRD 10
Bias = B (Co) or (rom)
Relative Bias = RB (%)oRoot Mean Square Error = RMSE (C ) or (rom)
Relative Root Mean Square Error = RRMSE (%)
Standard Deviation = SD (Co) or (rom)
Relative Standard Deviation = RSD (%)
Percent of YR MONTHS Absolute Differenceo -> (1 C or 10 rom)
Largest Negative Difference (Co) or (rom)
Largest Positive Difference (Co) or (rom)
Temperature Precipitation
0.99 -5.27
21.8 -16.8
1.26 17.69
27.7 56.5
14.1 64.9
49 81
59.6 29.8
-1. 2 -94.8
2.6 22.7
Percent of YR_MONTHS direction of changefrom the previous YR_MONTH in the CABvalues agrees with the actual (NCDC)values (%)
Pearson correlation coefficient betweenCAB and NCDC values
49
98
1.00
76
0.76
Appendix B, Part 3. Continued
CRD 20
Bias = B (Co) or (rom)
Relative Bias = RB (%)
Root Mean Square Error = RMSE (Co) or (mm)
Relative Root Mean Square Error = RRMSE (%)
Standard Deviation = SD (Co) or (rom)
Relative Standard Deviation = RSD (%)
Percent of YR MONTHS Absolute Difference> (1°C or 10 rom)
Largest Negative Difference (Co) or (rom)
Largest Positive Difference (Co) or (rom)
Temperature Precipitation
0.93 -19.89
21.3 -48.9
1.15 31. 82
26.4 78.2
0.68 24.84
12.9 119.5
51.1 57.8
-0.4 -91. 9
2.4 23.6
Percent of YR_MONTHS direction of changefrom the previous YR MONTH in the CABvalues agrees with the actual (NCDC)values (%)
Pearson correlation coefficient betweenCAB and NCDC values
50
100
1.00
75
0.69
Appendix B, Part 3. Continued
CRD 40
Bias = B (Co) or (mm)
Relative Bias = RB (%)
°Root Mean Square Error = fu~SE (C ) or (rom)
Relative Root Mean Square Error = RRMSE (%)
Standard Deviation = SD (Co) or (rom)
Relative Standard Deviation = RSD (%)
Percent of YR MONTHS Absolute Difference° -> (1 C or 10 rom)
Largest Negative Difference (Co) or (rom)
Largest Positive Difference (Co) or (rom)
Temperature Precipitation
0.56 -10.90
10.0 -29.9
1.24 21. 48
22.1 59.0
1.10 18.51
17.9 72.5
47.8 39.1
-2.3 -89.7
3.6 20.8
Percent of YR_MONTHS direction of changefrom the previous YR_MONTH in the CABvalues agrees with the actual (NCDC)values (%)
Pearson correlation coefficient betweenCAB and NCDC values
51
98
1.00
83
0.83
Appendix B, Part 3. Continued
CRD 50
Bias = B (Co) or (mm)
Relative Bias = RB (%)oRoot Mean Square Error = RMSE (C ) or (mm)
Relative Root Mean Square Error = RRMSE (%)
Standard Deviation = SD (Co) or (mm)
Relative Standard Deviation = RSD (%)
Percent of YR MONTHS Absolute Differenceo -> (1 C or 10 mm)
Largest Negative Difference (Co) or (mm)
Largest Positive Difference (Co) or (mm)
Temperature Precipitation
0.43 -17. 32
8.9 -43.3
0.89 27.87
18.4 69.7
0.78 21. 84
14.8 96.4
17.0 46.8
-0.8 -89.0
3.7 8.3
Percent of YRMONTHS direction of changefrom the previous YR_MONTH in the CABvalues agrees with the actual (NCDC)values (%)
Pearson correlation coefficient betweenCAB and NCDC values
52
98
1.00
89
0.79
Appendix B, Part 3. Continued
CRD 60
Bias = B (Co) or (mm)
Relative Bias = RB (%)oRoot Mean Square Error = RMSE (C ) or (mm)
Relative Root Mean Square Error = RRMSE (%)oStandard Deviation = SD (C ) or (mm)
Relative Standard Deviation = RSD (%)
Percent of YR MONTHS Absolute Differenceo -> (1 C or 10 mm)
Largest Negative Difference (Co) or (mm)
Largest Positive Difference (Co) or (mm)
Temperature Precipitation
0.36 -10.44
7.6 -25.6
1.02 19.87
21.4 48.6
0.96 16.90
18.6 55.5
17.0 46.8
-1. 6 -54.3
4.6 11. 9
Percent of YR_MONTHS direction of changefrom the previous YR_MONTH in the CABvalues agrees with the actual (NCDC)values (%)
Pearson correlation coefficient betweenCAB and NCDC values
53
96
1.00
65
0.88
Appendix B, Part 3. Continued
CRD 80
Bias = B (Co) or (mm)
Relative Bias = RB (%)oRoot Mean Square Error = RMSE (C ) or (mm)
Relative Root Mean Square Error = RRMSE (%)oStandard Deviation = SD (C ) or (mm)
Relative Standard Deviation = RSD (%)
Percent of YR MONTHS Absolute Differenceo -
> (1 C or 10 mm)
Largest Negative Difference (Co) or (mm)
Largest Positive Difference (Co) or (mm)
Temperature Precipitation
0.12 -6.93
2.0 -19.4
1.14 12.38
19.4 34.7
1.14 10.26
18.9 35.7
25.5 27.7
-1.5 -31. 3
4.4 3.7
Percent of YR_MONTHS direction of changefrom the previous YR_MONTH in the CABvalues agrees with the actual (NCDC)values (%)
Pearson correlation coefficient betweenCAB and NCDC values
54
96
1.00
89
0.96
Appendix B, Part 3. ContinuedCRD 90
Bias = B (Co) or (rom)
Relative Bias = RB (%)
°Root Mean Square Error = RMSE (C ) or (mm)
Relative Root Mean Square Error = RRMSE (%)
Standard Deviation = SD (Co) or (mm)
Relative Standard Deviation = RSD (%)
Percent of YR MONTHS Absolute Difference° -> (1 C or 10 rom)
Largest Negative Difference (Co) or (mm)
Largest Positive Difference (Co) or (mm)
Temperature Precipitation
0.67 -8.2611.5 -21.a
1.30 18.6722.4 47.4
1.11 16.74
17.2 53.8
34.1 43.2
-1.1 -58.4
4.8 22.6Percent of YR_MONTHS direction of change
from the previous YR_MONTH in the CABvalues agrees with the actual (NCDC)values (%)
Pearson correlation coefficient betweenCAB and NCDC values
55
98
1.00
77
0.86
Appendix C. Observed Significance Levels for Statistical Tests Referred to in Text
CRD Referredto on
10 20 30 40 50 60 70 80 90 page no.lo Shapiro-Wi1k tests of
normality on tempera-ture differences 0.849 <0.01 < 0.01 0.209 <0.01 <0.01 <0.01 <0.01 <0.01 3
2. Shapiro-Wi1k tests ofnormality on precipi-tation differences <0.01 <0.01 <0.01 <0.01 <0.01 <0.01 <0.01 <0.01 0.015 3
3. Nonparametric ANOVAon temperature valuesafter additive trans-
\.J1 formation 0.3822 0.2424 0.0204 0.5611 0.5611 0.2424 0.5611 0.1421 0.0236 140"-
4. Nonparametric ANOVA onprecipitation valuesafter additive trans-formation 0.0166 0.0765 0.0375 0.0375 0.1421 0.0375 0.0166 0.0375 0.1421 14
5. Nonparametric ANOVA onprecipitation valuesafter proportionaltransformation 0.3822 0.1421 0.7717 0.3822 0.2424 0.2424 0.7717 0.5611 0.2424 16
APPENDIX D
Formulae for measures of reliability for real-timeweather value estimates
Yi = Weather value as reported by NCDC for YR_MONTH i ("true" value)~Yi Weather value as estimated by CAB for YR MONTH i
d.1.
difference between CAB and NCDC values for YR-MONTH i
= 100 d·/Y· = relative difference for YR MONTH i1. 1.
i = 1, ... ,n
1Y = n I y.
1.
Bias = B
nnumber of YR MONTHs and I = i~l = summation over the YR MONTHs
l/n I di = d
Relative Bias = RB = 100 B/Y
Mean Square Error = MSE = l/n I d.21.
Root Mean Square Error = RMSE 1(MSE)2
Relative RMSE = RRMSE = 100 RMSE/Y-2Variance = Var = l/n I (di - d)
1Standard Deviation = SD = (Var)2
Relative SD = RSD = 100 SD/(Y + d)
(LY. ) (IY .) J+ B2 [ IT.y. 1. 1.MSE = Var 1. 1. n
~Pearson r between y. and Yi: r~1.
O:y .)2(L y2 J[d/ - ] ~y/1.
n
57
APPENDIX E
Estimates of equation parameters for regressions of CAB precipitationestimates on NCDC precipitation values for 47 months in North Dakota
CRD bO b1
10 2.8863 0.7395
20 0.6486 0.4950
30 -1. 0837 0.5782
40 4.6596 0.5728
50 2.7011 0.4992
60 -1. 6531 0.7849
70 -2.0513 0.7557
80 0.2572 0.7985
90 4.0424 0.6876
58
*u.s. GOVERIiIIElIT PRII'ITIliG OFPICE : 1984 D-42Q-931/SRS-727