A SUPER-OB FOR THE WSR-88D RADAR RADIAL WINDS FOR USE IN THE NCEP OPERATI ONAL ASSIMILATION SYSTEM. NOAA/NWS/NCEP, Camp Springs, Maryland 20746

Similar documents
13A. 4 Analysis and Impact of Super-obbed Doppler Radial Velocity in the NCEP Grid-point Statistical Interpolation (GSI) Analysis System

Implementation and Evaluation of WSR-88D Radial Velocity Data Assimilation for WRF-NMM via GSI

The use of WSR-88D radar data at NCEP

9.4 INITIAL DEPLOYMENT OF THE TERMINAL DOPPLER WEATHER RADAR SUPPLEMENTAL PRODUCT GENERATOR FOR NWS OPERATIONS

P4.11 SINGLE-DOPPLER RADAR WIND-FIELD RETRIEVAL EXPERIMENT ON A QUALIFIED VELOCITY-AZIMUTH PROCESSING TECHNIQUE

Robert E. Saffle * Mitretek Systems, Inc., Falls Church, VA

Robert E. Saffle * Noblis, Inc., Falls Church, VA. Michael J. Istok National Weather Service, Office of Science and Technology, Silver Spring, MD

Multi-Sensor Precipitation Reanalysis

P r o c e. d i n g s. 1st Workshop. Madrid, Spain September 2002

7 AN ADAPTIVE PEDESTAL CONTROL ALGORITHM FOR THE NATIONAL WEATHER RADAR TESTBED PHASED ARRAY RADAR

results, azimuthal over-sampling and rapid update Y. Umemoto 1, 2, T. Lei 3, T. -Y. Yu 1, 2 and M. Xue 3, 4

P10R.10 NWS Use of FAA Radar Data Status and Operational Considerations

4.3. David E. Rudack*, Meteorological Development Laboratory Office of Science and Technology National Weather Service, NOAA 1.

7 WSR-88D OBSERVATIONS OF AN EXTREME HAIL EVENT IMPACTING ABILENE, TX ON 12 JUNE 2014

THE IMPACT OF SATELLITE-DERIVED WINDS ON GFDL HURRICANE MODEL FORECASTS

The use of WSR-88D radar data at NCEP

IMPACT OF GROUND-BASED GPS PRECIPITABLE WATER VAPOR AND COSMIC GPS REFRACTIVITY PROFILE ON HURRICANE DEAN FORECAST. (a) (b) (c)

Joseph C. Lang * Unisys Weather Information Services, Kennett Square, Pennsylvania

P3.1 Development of MOS Thunderstorm and Severe Thunderstorm Forecast Equations with Multiple Data Sources

Estimating the Impact of a 3-dB Sensitivity Loss on WSR-88D Data

P4.8 PERFORMANCE OF A NEW VELOCITY DEALIASING ALGORITHM FOR THE WSR-88D. Arthur Witt* and Rodger A. Brown

LEVERAGING MULTIPLE FAA RADARS FOR NWS OPERATIONS

5B.2 AN UPDATE ON THE NWS WSR-88D LEVEL II DATA COLLECTION AND DISTRIBUTION NETWORK AND PLANS FOR CHANGES

1. INTRODUCTION 3. VERIFYING ANALYSES

THE DETECTABILITY OF TORNADIC SIGNATURES WITH DOPPLER RADAR: A RADAR EMULATOR STUDY

USE OF SURFACE MESONET DATA IN THE NCEP REGIONAL GSI SYSTEM

Allison Monarski, University of Maryland Masters Scholarly Paper, December 6, Department of Atmospheric and Oceanic Science

Meteorology 311. RADAR Fall 2016

Assimilation of Satellite Infrared Brightness Temperatures and Doppler Radar Observations in a High-Resolution OSSE

USING GRIDDED MOS TECHNIQUES TO DERIVE SNOWFALL CLIMATOLOGIES

P5.4 WSR-88D REFLECTIVITY QUALITY CONTROL USING HORIZONTAL AND VERTICAL REFLECTIVITY STRUCTURE

ASSESMENT OF THE SEVERE WEATHER ENVIROMENT IN NORTH AMERICA SIMULATED BY A GLOBAL CLIMATE MODEL

A Realtime Forecasting System for the Chesapeake Bay Tributaries: Interoperability and Metadata for Automation

WARNING DECISION SUPPORT SYSTEM INTEGRATED INFORMATION (WDSS-II). PART I: MULTIPLE-SENSOR SEVERE WEATHER APPLICATIONS DEVELOPMENT AT NSSL DURING 2002

Prediction of Snow Water Equivalent in the Snake River Basin

Doppler radial wind spatially correlated observation error: operational implementation and initial results

1. INTRODUCTION 2. QPF

*Corresponding author address: Charles Barrere, Weather Decision Technologies, 1818 W Lindsey St, Norman, OK

BUILDING A GRIDDED CLIMATOLOGICAL DATASET FOR USE IN THE STATISTICAL INTERPRETATION OF NUMERICAL WEATHER PREDICTION MODELS

Improved analyses and forecasts with AIRS retrievals using the Local Ensemble Transform Kalman Filter

P6.18 THE IMPACTS OF THUNDERSTORM GEOMETRY AND WSR-88D BEAM CHARACTERISTICS ON DIAGNOSING SUPERCELL TORNADOES

IOP-2 Summary of Operations 03 June UTC 04 June UTC Authors: Market, Kastman

2.7 A PROTOTYPE VERIFICATION SYSTEM FOR EXAMINING NDFD FORECASTS

P 5.16 Documentation of Convective Activity in the North-eastern Italian Region of Veneto

5.2 PRE-PROCESSING OF ATMOSPHERIC FORCING FOR ENSEMBLE STREAMFLOW PREDICTION

Ensemble-variational assimilation with NEMOVAR Part 2: experiments with the ECMWF system

TEMPORAL DISTIRUBTION OF PMP RAINFALL AS A FUNCTION OF AREA SIZE. Introduction

5B.1 OBSERVATIONS OF TORNADOGENESIS USING A MOBILE, PHASED-ARRAY, DOPPLER RADAR

Ionospheric Scintillation Impact Report: South African SKA Site

GIS in Water Resources. Fall Homework #1

Utilizing Ensemble-based Anomalies to Anticipate Significant Events A New Display Tool

P7.1 STUDY ON THE OPTIMAL SCANNING STRATEGIES OF PHASE-ARRAY RADAR THROUGH ENSEMBLE KALMAN FILTER ASSIMILATION OF SIMULATED DATA

Fundamentals of Radar Display. Atmospheric Instrumentation

Enhancing information transfer from observations to unobserved state variables for mesoscale radar data assimilation

Impact of METOP ASCAT Ocean Surface Winds in the NCEP GDAS/GFS and NRL NAVDAS

Scott A. True Project Scientist Geodesy & Geophysics Division Basic and Applied Research Office InnoVision Directorate

J16.1 PRELIMINARY ASSESSMENT OF ASCAT OCEAN SURFACE VECTOR WIND (OSVW) RETRIEVALS AT NOAA OCEAN PREDICTION CENTER

Town of Chino Valley. Survey Control Network Report. mgfneerhg mc N. Willow Creek Road Prescott AZ

Single-Doppler EnKF assimilation of high-resolution data from the 29 May 2004 OKC supercell: Comparisons with dual-doppler analyses.

P1.3 STATISTICAL ANALYSIS OF NUMERICAL MODEL OUTPUT USED FOR FORECASTING DURING JOINT URBAN 2003

RSMC WASHINGTON USER'S INTERPRETATION GUIDELINES ATMOSPHERIC TRANSPORT MODEL OUTPUTS

UNIT 1C. USING TOPOGRAPHIC MAPS WHERE IN THE WORLD... ARE YOU?

RADAR- BASED ANALYSIS OF THE DALLAS COWBOYS PRACTICE FACILITY COLLAPSE OF 2 MAY 2009

Maps. 3. Which altitude of Polaris could be observed in New York State?

Analysis and Forecasting of the Low-Level Wind during the Sydney 2000 Forecast Demonstration Project

Topographic Maps and Profiles

1C.4 TROPICAL CYCLONE TORNADOES: SYNOPTIC SCALE INFLUENCES AND FORECASTING APPLICATIONS

P1.6 Simulation of the impact of new aircraft and satellite-based ocean surface wind measurements on H*Wind analyses

Error Covariance Estimation and Representation for Mesoscale Data Assimilation

CHARACTERIZATION AND MITIGATION OF WIND TURBINE CLUTTER ON THE WSR-88D NETWORK

Map Projections. Displaying the earth on 2 dimensional maps

A VAD-Based Dealiasing Method for Radar Velocity Data Quality Control

Current Limited Area Applications

Asymmetry in Wind Field of Typhoon 0115 analyzed by Triple Doppler Radar Observation

4.6 NWS WSR-88D AND TDWR-SPG DATA COLLECTION AND DISTRIBUTION NETWORK STATUS AND PLANS

10.5 PROBABLISTIC LIGHTNING FORECASTS AND FUEL DRYNESS LEVEL FORECASTS IN THE GRAPHICAL FOREAST EDITOR: EXPANDED DOMAIN AND DISTRIBUTION FOR 2009

Hurricane Wave Topography and Directional Wave Spectra in Near Real-Time

Meteorological Studies with the Phased Array Weather Radar and Data Assimilation Using the Ensemble Kalman Filter

SNOWFALL RATE RETRIEVAL USING AMSU/MHS PASSIVE MICROWAVE DATA

P1.12 MESOSCALE VARIATIONAL ASSIMILATION OF PROFILING RADIOMETER DATA. Thomas Nehrkorn and Christopher Grassotti *

Development of Methodology and Tools for Determining the Impact of Cloud-Cover on Satellite Sensors

Ed Tomlinson, PhD Bill Kappel Applied Weather Associates LLC. Tye Parzybok Metstat Inc. Bryan Rappolt Genesis Weather Solutions LLC

SAMPLE. SITE SPECIFIC WEATHER ANALYSIS Wind Report. Robinson, Smith & Walsh. John Smith. July 1, 2017 REFERENCE: 1 Maple Street, Houston, TX 77034

1.2 DEVELOPMENT OF THE NWS PROBABILISTIC EXTRA-TROPICAL STORM SURGE MODEL AND POST PROCESSING METHODOLOGY

Synoptic Meteorology II: Frontogenesis Examples Figure 1

NAM-WRF Verification of Subtropical Jet and Turbulence

CSU RAMS. Standard Input Formats for Pressure Coordinate and Observational Data (Data Prep Stage)

Localized Aviation Model Output Statistics Program (LAMP): Improvements to convective forecasts in response to user feedback

Inter-tropical Convergence Zone (ITCZ) analysis using AIRWAVE retrievals of TCWV from (A)ATSR series and potential extension of AIRWAVE to SLSTR

THE CURRENT STAGE OF DEVELOPMENT OF A METHOD OF PRODUCING MOTION VECTORS AT HIGH LATITUDES FROM NOAA SATELLITES. Leroy D. Herman

A Climatology of supercells in Romania

WindNinja Tutorial 3: Point Initialization

Forecasting of Optical Turbulence in Support of Realtime Optical Imaging and Communication Systems

Recent Data Assimilation Activities at Environment Canada

Long term analysis of convective storm tracks based on C-band radar reflectivity measurements

Assimilation of Doppler radar observations for high-resolution numerical weather prediction

Simon Tschannett*, B. Bica, B. Chimani, M. Dorninger, W. Gepp, S. Schneider, R. Steinacker University of Vienna Vienna, Austria

The Local Ensemble Transform Kalman Filter (LETKF) Eric Kostelich. Main topics

Developmental Testbed Center Visitor Program Final Report: Development and Application of 3-Dimensional Object Algorithms to High Resolution Forecasts

Can the Convective Temperature from the 12UTC Sounding be a Good Predictor for the Maximum Temperature, During the Summer Months?

Transcription:

19th Conference on IIPS P151 A SUPER-OB FOR THE WSR-88D RADAR RADIAL WINDS FOR USE IN THE NCEP OPERATI ONAL ASSIMILATION SYSTEM Jordan C Alpert 1*, Peter Pickard 2, Yukuan Song 2, Melanie Taylor 2, William B Facey 1, Michael Istok 2 and David F Parrish 1 1 NOAA/NWS/NCEP, Camp Springs, Maryland 20746 2 NWS/OST/SEC, Silver Spring, Maryland 20910 1 INTRODUCTION The spatial and temporal densities of WSR-88D raw radar radial wind represent a rich source of high resolution observations A characteristic of these observations is the presence of a significant degree of redundant information present in the radar returns Redundant observations impose a burden on an operational assimilation system since each datum is processed with repetitive interpolations from the analysis grid to its location and back again This effort is carried out for each datum regardless of the information that can be attributed to it in the overall assimilation The time and storage expended on mutually redundant data could be better spent on improving other aspects of an assimilation (Purser, et al, (2000) Therefore, it is desirable to effect whatever data compression the ensemble of fresh observations allow while minimizing any degradation of the information content The term for a surrogate datum which replaces several partially redundant actual data is a "super observation" or "super-ob" Super-obs have been applied to operational analysis at the National Center for Environmental Prediction (NCEP) Operational assimilation system for subsets of the WSR-88D radar radial wind observations from the NEXRAD Information Dissemination Service (NIDS) However, the precision and information content of the NIDS radial wind can be improved if data at each radar site is directly utilized at the resolution and precision of the WSR-88D radar to construct the super-ob and then delivered to a central site We report a new product for the Open RPG software system called the Radar Superob which takes the result of radar scans and averages data points within a prescribed time and spatial 3- dimensional volume before transmitting reports to a central collection of base data The new product will be super-obs of radar radial winds * Corresponding author address: Jordan C Alpert, DOC/NOAA/NWS/NCEP/Environmental Modeling Center, 5200 Auth Road, Rm 207, Camp Springs, MD 20746; e-mail: jordanalpert@noaagov 2 DESCRIPTION The full-resolution WSR-88D base radial wind data provides sufficient data am ounts for statistically significant averaging The super-ob product is programmed using the Open Systems RPG (Radar Product Generator) to control all aspects of the calculation at WSR-88D radar sites The Open RPG is the system which operates between the Open Radar Data Acquisition (RDA) system and more sophisticated display devices, such as the National Weather Service's Advanced Weather Interactive Processing System (AWIPS) The Open RDA collects data from the WSR-88D radar and forwards base data products to the Open RPG These base data products consist of reflectivity, radial wind, and spectrum width The Open RPG creates the special purpose products from the base data and forwards them to other systems for display or for further processing Super-ob is one of the newly enabled, enhanced products under the Open RPG Adaptable parameters for super-ob are the Time window, Cell Range Size, Cell Azimuth Size, Maximum Range and Minimum Number of Points The values of these parameters, which define the super-ob averaging, are shown in Table 1 The default settings indicate that at each elevation angle, a wedge shape volume of 6 Table 1 Adaptable parameters for the Super-ob product Parameter Default Range Time Window 60 minutes [5-90 min] Cell Range Size 5 km [1-10 km ] Cell Azimuth Size 6 degrees [2-12 deg] Maximum Range 100 km [60-230km] Minimum Number of Points required 50 unitless [20-200] azimuth degrees by 5 km along a radius, averaged over a time of 60 minutes define each super-ob cell Each of the super-obs contain no less than 50 points and no cell would extend past 100 km, as the radar beam width becomes too wide and returns become less certain at

larger distances from the radar These adaptable parameter values are programmed through the Open System RPG (ORPG) The range of possible values will allow for the super-ob product to adjust to different analysis resolution requirements All the WSR-88D installations will contribute super-ob data to initialize cycling analysis systems, with identical settings through the ORPG The data transmission precision of the super-ob (mean) radial wind is ±001 m/s or 32767 possible levels This is an improvement over the current transmission of 15 levels of data transmission precision for winds which can span ±100 m/s The standard deviation of the mean radial wind super-ob is transmitted at 15 levels of precision 3 RESULTS A number of tests have been made of the radar system within the Systems Engineering Center at the NOAA/NWS Office of Science and Technology An example of the Radial wind scalar is shown on a vector plot in Fig 1 where the arrow direction is toward or away from the radar and colors, as well as arrow length, indicate the speed of the return The Data for this example is from KATX, the Seattle, WA radar from May 28, 2002 at 1800Z and represents a typical return The characteristic pattern of radar radial returns is apparent and there are large areas that do not have returns At larger distances from the radar the radial wind is seen to increase southwest and northeast of the radar Fig 2 Super-ob height [m] for each observation beginning at 0 az imuth (north ) and e xten ding th rough larger radii and 360 degrees (clockwise) for the same time as in Fig 1 Fig 3 Super-ob radial wind Standard deviation [m/s] for each obser vation b eginning at 0 azimuth (north) and extending through larger radii and 360 degrees (clockwise) for the same time as in Fig 1 Fig 1 Radar radial (scalar) wind super-obs speed and location from station KATX in Seattle, WA on M ay 28, 2002 Arrow s are toward or away form the radar with the color and arrow length indicating the spe ed of the r eturn Fig 4 Super-ob radar radial (mean) wind [m/s] for each observation beginning at 0 azimuth (north) and extending through larger radii and 360 degrees (clockwise) for the same time a s in Fig 1

due to synoptic conditions However, a number of superob wind values in Fig 4 extend beyond the normal range as defined by its neighbors One can scan this data set and find values outside the range of the super-ob wind and compare with the standard deviation, height and the location of a return in question Fig 5 Super-ob radar radial wind location [longitude & latitude] for each observation for the same as in Fig 1 In a clear air test at station KCRI, the Norman, OK test radar, (not shown) the coverage was more complete The number of super-ob returns near the radar has not been reduced at this time although there is code to accomplish this and may be done in the future As the radar beam extends to larger radial distance from the radar, the height of the super-ob increases In the case of the radial wind shown in Fig 1, at the lowest elevation angle, the height ranges from about 250 m to 2 km Higher elevation angles have heights that extend from 1 km to throughout the entire troposphere One way to show the height variation is to plot the super-ob data set according to the first observation which is along 0 degrees north closest to the radar extending out to larger distance from the radar, and then clockwise around in azimuth for the remaining returns in the report, in this case over 700 super-obs The height in meters is displayed this way in Fig 2 Each point is the result of the super-ob at the next 5 km radius increasing its height as the radial distance increases and then clockwise from north around with increasing azimuth Fig 4 shows each super-ob radial wind plotted in this fashion where a sinusoidal pattern of the radial wind results due to the radial wind definition Each successive super-ob moves out alone the radii with increasing beam height causing a range of super-ob wind magnitudes along each azimuth The super-ob returns labeled mean wind in Fig 4 range from a few m/s near point 200 to 10 m/s near point 650 This range of wind values is from varying height and location of the return and the result of natural variance in the winds in the radar vicinity Radial wind speeds at the top of the height range, for each radar beam elevation angle, in Fig 2 (elevation angle = 04), sometimes appear less correlated with their neighbors (Fig 4) For example, point 108 at 4885 at (lat, lon), (4885, -121632) have a super-ob wind of 121 m/s at a height of 1549 m Points 458 and 474 at (47263, -123238) and (47579, -123110) respectively, located on Fig 1, show radial wind magnitudes of -1412 m/s and -1258m/s on Fig 4 which is outside the range of common values The returns for these super-obs are expected to have larger wind magnitudes with increased elevation The variance for these points is low, 2 m/s, for each point An example of a high standard deviation return is shown for point 589 at (48039, -123789) with a standard deviation of 12 m/s and wind magnitude of -1033 m/s The height for point 589 is 1572 m which makes it among the larger distances from the radar and higher height for this radar elevation angle Fig 3 shows a moderate spike at point 589 in line with the increased height in Fig 2 and the increased negative super-ob wind In this case, the value is suspect because of the large standard deviation in the radar radial winds making up the super-ob One may expect a larger standard deviation of radial wind values because of the widening of the radar beam width with increasing distance from the radar as well as with increased wind speeds at higher altitudes When the radar beam grows wider with increasing distance from the radar, the vertical position of the return can not be determined with precision The standard deviation, although transmitted as a 4 bit marker (15 levels), can be useful in determining the quality of the returns The quality control of the super-ob data in the analysis program may also be able to utilize the reflectivity to determine locations of severe convection (and vertical motion) which could cause errors in the radial wind returns The standard deviation can also be large at points that have less than 1 m/s winds such as points 630, 638, and 639 and 634 In this case these were among the nearly calm values located north and west of the radar (Fig 1) Fig 5 illustrates the location of each super-ob The centroid of WSR-88D returns making up each superob is calculated and made part of the transmitted header of each super-ob value If actual returns are weighted toward the longitude and latitude position then the superob position will reflect this This is the reason the superob locations, shown by a + in Fig 5 or the tails of the vectors in Fig 1, are not orientated in exact concentric circles This is also done in time as several scans are averaged to construct the super-ob (Table 1) The time

deviation, accounting for the number of raw radar returns differing in each scan are transmit ted in the super-ob header Close to the radar, the raw return numbers are sufficient to create super-obs under the criterion of Table 1 even though the super-ob area grows smaller At higher resolution (smaller super-ob area) the near radar super-obs will not have enough returns to satisfy the minimum number of points required to form the super-ob This criterion may be adjusted so the number of super-obs will be thinned out or the large number of super-obs in close to the radar shown in Fig 5 can be decreased by directly removing super-obs as a function of radial distance from the radar 4 ANALYSIS AND VERIFICATION The super-ob examples shown herein are typical of a number of examined preliminary data sets Fig 6 shows analyzed regular winds at a height of 04 km from the NCEP Global Forecast System (GFS) Spectral Statistical Interpolation (SSI) Analysis of May 5, 2002 at 18Z The wind barbs are displayed over the same domain and time as the super-ob radial wind observations in Fig 1 The operational SSI uses all available remote and conventional observations to create the analysis at a resolution (spectral triangular truncation T170) of less than 1 degree of latitude and longitude It is plotted at 1 degree in the figure We can use the SSI analysis as a proxy for observed winds in the radar area as the number of in conventional wind observations in the area depicted around the radar are few The black cross hair in Fig 6 is located at the Seattle, WA radar location The regular winds shown in Fig 6 are at a height of 04 km which was interpolated vertically by the SSI analysis The super-ob radial winds located close to the radar are at or near a height of 04 km as shown in Fig 2 In viewing the analyzed regular winds, composed of wind components (u,v), in Fig 6, one must keep in mind the definition of radial wind: u radial = u cos + v sin is a scaler, where is the angle between a line from the radar directed at 0 degrees azimuth (north) and the radial wind azimuth There is one direction where the regular wind will be identical to the radial wind and regular winds perpendicular to this will have zero radial wind The direction of the regular wind north and west of the radar in Fig 6 is southwest at 15 knots The radial wind along a line northwest of the radar line up in the same direction The radial wind speeds in the northeast Fig 6 Analyzed regular Wind barbs (blue) at 04 km height from the NCEP Global Forecast System Spectral Statistical Interpolation Analysis of May 5, 2002 at 18Z The black cross hair is at the Seattle, WA radar location quadrant vary from 5-15 knots close in as shown by the blue and green wind arrows in Fig 1 To the east of the radar, the wind magnitudes are far less, 0-5 knots This is because of the radial wind definition and the winds in this region being more southerly and of lower magnitude (Fig 6) so their component in the radial direction will be small The analyzed winds in the southwest quadrant (Fig 6) tend toward larger magnitudes as do the radial winds shown on Fig 1 Plots of 1 km and 15 km regular analyzed winds (not shown) show that winds at higher elevations in this area have increased magnitudes even though the angle of the wind (southerly) and the radial wind are at large angles The regular winds in the southeast quadrant are SSE with similar magnitudes to the radial wind That we show consistent magnitudes of the radial wind along the SE - NW direction as the observed winds is a necessary condition for acceptance of the data by the analysis system The analyzed regular winds in the southeast and northwest quadrants show changes that are reflected in the radial wind pattern as a turning of the regular wind direction from south to southwesterly This is an example of the high resolution potential of the radar radial wind data set

4 SUMMARY AND FUTURE WORK Implementation of the new super-ob product is scheduled for the fall of 2002 Once the radial wind super-obs are received at NCEP they will be used in the analysis cycling system to become a component, with the other remotely sensed and conventional observations, providing a high resolution wind component In addition to tests run with and without the presence of the radial winds from the full complement of WSR-88D radars, Purser, et al, (2000), have shown an objective means of quantifying the information lost by the conversion of many raw data into few super-obs Refinements about the decisions made on the size of clusters of raw data in space and time can be rationalized using methods of information theory as described in the Purser, et al (2000) approach A comparison with the NIDS form of the data will be also done Work continues to extend the super-ob to distances beyond 100 km from the radar Data and plots of other elevations and radar stations are available at the web site: http://sgi62wwb noaagov:8080/superob The raw data for this report is located at: http://sgi62wwbnoaagov:8080/superob/read_superob The pre-print figures shown in this report are located at: http://sgi62wwbnoaagov:8080/superob/amsfigures 4 REFERENCES Purser, R J,, D Parrish and M Masutani, 2000: Meteorological observational data compression; an alternative to conventional super-obbing, Office Note 430, DOC/NOOA/NWS/NCEP, 5200 Auth Road, Camp Springs, MD 20746-4304,12pp