National Weather Service United States Department of Commerce
To be presented at International Association for Hydraulic Research (IAHR)
XXVII Congress
San Francisco, California
August 10-15, 1997

Comparison of National Weather Service Operational
Mean Areal Precipitation Estimates Derived from
NEXRAD Radar vs. Rain Gage Networks


Bryce Finnerty
Dennis Johnson
Office of Hydrology
NOAA/National Weather Service
1325 East-West Highway
Silver Spring, Maryland 20910

ABSTRACT

Historically, hydrologic forecasts and calibration of the National Weather Service (NWS) hydrologic models are typically prepared using inputs of 6-hour mean areal precipitation (MAP) estimates derived from rain gage networks. However, operational hydrologic forecasts can now be prepared using 1- or 6-hour MAPs derived from high resolution gridded precipitation estimates from the NWS Next Generation Radar system, NEXRAD. An initial analysis of 7 months (May 93 - December 93) of 1- and 6-hour operational MAP time series for nine basins revealed differences between the gage and radar MAPs in the long term accumulations, individual storm totals, and timing of events. The 7-month radar MAP accumulations were 10-25 percent less than the gage MAP accumulations. The 6-hour radar and gage MAPs had similar estimates of the timing of events, but the 1-hour radar and gage MAPs showed more discrepancies in the timing of events. The radar MAPs captured more of the variability in the precipitation fields than the gage MAPS at the 1-hour time step. However, the variability of the radar and gage MAPs was nearly equal at the 6-hour time step.

INTRODUCTION

New challenges in operational hydrologic forecasting have materialized because of the installation of an advanced national system of radars called Weather Surveillance Radar 1988-Doppler (WSR88-D), or NEXRAD (Hudlow, 1988). NEXRAD produces high resolution gridded precipitation estimates at a 4x4 km2 spatial scale and a 1-hour time scale. These NEXRAD data provide spatial and temporal estimates of highly variable precipitation fields that are not possible from point measurements of precipitation from rain gage networks.

Hydrologic forecasting at the NWS generally uses a lumped modeling approach in which hydrologic model parameters are averaged over a river basin. The Sacramento Soil Moisture Accounting (SAC-SMA) model (Burnash et al.,1973), which is the most commonly used hydrologic model in the NWS, uses a mean areal precipitation (MAP) input and outputs discharge. The SAC-SMA model parameters are calibrated using 6-hour model input MAPs derived from historical point gage precipitation measurements. These parameters are then used operationally with 6-hour gage MAPs to produce river forecasts with a 6-hour time step. River forecasters now have the option of using 1- or 6-hour radar-derived mean areal precipitation (referred to as MAPXs) for operational forecasting. However, the model parameters have not been recalibrated to the MAPX data because only a short 7-month period exists where historical gage and discharge data overlap with the NWS archive of radar data. Therefore, a detailed analysis of model input MAPs is required to understand how these model inputs will impact model performance. The objective of this paper is to analyze the differences in operational methods of estimating 1- and 6-hour MAPs derived from radar and gage networks.

ANALYSIS

The study area is nine basins in the region near the Oklahoma-Arkansas-Missouri state boundaries. This region was analyzed because of its dense gage network, six overlapping radar umbrellas, and one of the longest available periods of archived NEXRAD radar products in the United States to date. MAP and MAPX time series were derived for the period from May 7, 1993 through December 31, 1993. This is the period in which the historical gage and discharge data overlap with archived radar products.

The radar data used for the study are the 4x4 km2 resolution Stage III data. Stage III is a merged radar-gage precipitation field design to provide the spatial resolution of radar data while preserving the precipitation accumulations measured by gages. It assumes the gage data as "ground truth" and scales the radar accumulation estimates to match the gages (Shedd and Smith, 1991). The influence of the gages on the Stage III data diminish with distance from the gage. The gage MAPs were calculated using both the hourly and daily gage accumulations obtained from the National Climatic Data Center's 3200 precipitation data.

Summary statistics were generated for the 1- and 6-hour case, for all nine basins, and for the entire 7-month duration of the MAP and MAPX time series.

Statistics include the mean, standard deviation, and coefficient of variation, all of which are conditioned upon the occurrence of precipitation. The conditional statistics provide information about the ability of radar and gage networks to detect precipitation, as well as their ability to estimate precipitation rates. Cumulative sums for each basin and the percent bias of accumulations from radar vs. gage networks were also tabulated.

RESULTS

Tables 1 and 2 show that for the 7-month period, the Stage III MAPX cumulative totals have a negative bias ranging from -10 to -25 percent as compared to gage MAPs. This under catch of the MAPXs exists for both the 1- and 6-hour cases. The fact that the 6-hour conditional mean of the MAPX is larger than MAP for some basins is attributed to the spatial and temporal averaging of point gage measurements to time increments and areas where it is not raining. This produces rain for time intervals where no rain was occurring, reduces the conditional mean, and hides the actual cumulative bias that exist at the 6-hour case. MAPX is approximately 25 percent greater at detecting the variability in the precipitation at the 1-hour time step than the gage MAPS. However, the normalized mean (coefficient of variation) shows that the variability of the MAPs and MAPXs is nearly equal at the 6-hour time step. This indicates that the radar's ability to capture the variability in precipitation at a 1-hour time step is largely lost when averaging to a 6-hour time step.

Figure 1 illustrates how the cumulative sum of both the MAP and MAPX time series behaves over the 7-month duration for the TIFM7 basin. Only this basin is shown because all nine basins exhibited similar behavior. The figure shows how the gages and radar have different estimates of the volume and timing of individual precipitation events. Figure 1 also shows a plot of the cumulative differences in the MAP and MAPX time series, which highlights the timing and storm-by-storm catch of the two rainfall measurements. Both methods appear to be similar in their estimation of event timing, but a spike appears in the cumulative difference plot when their timing is off. MAPXs have a lower accumulation even in the summer, which may be due to averaging point gage measurements over areas where it is not raining, which overestimates precipitation. Or, the radar may simply be underestimating precipitation rates. Radar generally performs better than gage networks during high intensity convective storms because of the radar's high resolution spatial coverage (Seo and Smith, 1996; Smith et al., 1996).

Table 1. Conditional Summary Statistics for 9 Basins Using 1-hour MAP and MAPX Time Series
Basin ID SUM
(mm)
HOURLY MEAN
(mm)
HOURLY STD. DV.
(mm)
COEF. of
VARIATION
% BIAS
of
SUMS
MAP3 MAPX MAP3 MAPX MAP3 MAPX MAP3 MAPX
JOPM7 1207 1040 2.07 1.51 3.26 3.26 1.575 2.160 -13.0
TIFM7 1104 996 1.48 1.32 2.69 2.66 1.820 2.020 -9.8
WTTO2 977 798 1.27 1.04 2.68 2.20 2.107 2.111 -18.3
KNSO2 948 837 1.50 1.41 2.67 2.83 1.782 2.000 -11.7
ELDO2 1056 930 1.51 1.34 2.55 2.90 1.680 2.176 -11.9
TALO2 933 839 1.47 1.33 2.77 3.28 1.880 2.460 -10.1
TENO2 979 784 1.62 1.18 2.85 2.90 1.755 2.460 -19.9
VLBA4 1055 798 1.73 1.07 2.99 2.29 1.730 2.128 -24.4
MLBA4 864 726 1.523 1.11 2.37 2.48 1.547 2.240 -15.9
% BIAS = [MAPX - MAP) / MAP] X 100

Table 2. Conditional Summary Statistics for 9 Basins Using 6-hour MAP and MAPX Time Series
Basin ID SUM
(mm)
6-HOUR MEAN
(mm)
6-HOUR STD. DV.
(mm)
COEF. of
VARIATION
% BIAS
of
SUMS
MAP3 MAPX MAP3 MAPX MAP3 MAPX MAP3 MAPX
JOPM7 1207 1040 5.51 5.10 9.29 9.52 1.69 1.87 -13.9
TIFM7 1104 996 4.30 4.70 7.90 8.55 1.84 1.82 -9.8
WTTO2 977 798 3.60 3.82 6.94 6.89 1.93 1.80 -18.3
KNSO2 948 837 4.35 4.65 7.75 8.32 1.78 1.79 -11.7
ELDO2 1056 930 4.19 4.72 7.68 8.79 1.83 1.86 -11.9
TALO2 933 839 4.15 4.49 8.34 9.34 2.01 2.08 -10.1
TENO2 979 784 4.45 4.15 8.64 9.02 1.94 2.17 -19.9
VLBA4 1055 798 4.67 3.80 7.69 6.20 1.65 1.63 -24.4
MLBA4 864 726 4.23 3.65 6.58 6.77 1.56 1.85 -15.9


Figure 1. 7-month cumulative sum of 1-hour MAP and MAPX time series for the TIFM7 basin. Time=0 is May 7, 1993, and time=5736 is December 31, 1993. The radar MAPX has a lower cumulative sum of precipitation for the period analyzed. The bottom of the figure shows a plot of the difference between gage and radar mean areal precipitation estimates (MAP - MAPX).

CONCLUSIONS

The most significant result of this preliminary study is that MAPX has a negative bias in the 7-month accumulations for all nine basins analyzed. Any bias in the hydrologic model inputs may impact operational flood forecasts, soil moisture accounting, and long-term water supply forecasts. The 1-hour MAPX detected more precipitation variability than the gage MAPs, but this effect was minimized when averaging to a 6-hour time step. There were clear differences in each method's ability to estimate the timing and volume of the same particular event; however, no clear trend existed for all events. This indicated that hydrologic forecasting can benefit from having both measurement systems available to the forecasters, from which they can choose. Further analysis is required to evaluate what impact these differences between MAP and MAPX model inputs may have on operational hydrologic forecasting.

REFERENCES

Burnash, R.J.C., Ferral, R.L., and McGuire, R.A., 1973. "A generalizedstreamflow simulation system - conceptual modeling for digital computers," U.S. Department of Commerce, National Weather Service and State of California, Department of Water Resources.

Seo, D.J., and Smith, J.A., 1996. "Characterization of climatological variability of mean areal rainfall through fractional coverage," Water Resources Research, 33(7): 2087-2095.

Shedd, R.C., and Smith, J.A., 1991. "Interactive precipitation processing for the modernized National Weather Service," Preprints, Seventh International Conference on Interactive Information and Processing Systems for Meteorology, Oceanography, and Hydrology, New Orleans, Louisiana, American Meteorological Society, pp. 320-323.

Smith, J.A., Seo, D.J., Baeck, M.L., and Hudlow, M.D., 1996. "An intercomparison study of NEXRAD precipitation estimates," Water Resources Research, 32(7): 2035-2045.



Main Link Categories:
Home | OHD | NWS