RICO Sounding Composite Highest Resolution Upper-Air Data Set

1.0 General Description

This dataset contains the highest possible resolution (i.e. the 'native' resolution) upper-air sounding data in the National Center for Atmospheric Research (NCAR) Earth Observing Laboratory (EOL) Sounding Composite Format from stations in the following networks:

Data for the Rain in Cumulus Over the Ocean (RICO) 2004-2005 domain (15°N to 21°N latitude and 75°W to 72°W longitude) are contained within this dataset. The dataset covers the time period 1 December 2004 through 31 January 2005. This RICO 2004-2005 Sounding Composite Highest Resolution Upper-Air Dataset contains 554 soundings from 4 sites.

Section 2.0 contains a detailed description of the source network that collected the data. Section 2.1 contains a detailed description of the format of the composite dataset. See Section 3.0 below for the quality control processing performed by NCAR/EOL on this dataset. Section 4.0 contains references.

2.0 Detailed Data Description


This dataset contains 144 high vertical resolution (1 second) soundings from various locations around Barbuda. Soundings were typically released twice per day from 07 December 2004 to 08 January 2005 and four times per day from 09 January to 24 January 2005. Further information can be found in the NCAR GAUS README file. Note that this composite data set includes only the ascending data files.

2.0.2 NCAR/NSF C-130 Dropsonde

This dataset contains 196 high vertical resolution (1/2 second) soundings from the NCAR/NSF C-130 aircraft. During RICO the C-130 conducted eighteen missions and one ferry flight completed near Antigua. The duration of each flight was approximately nine hours long and each started and ended with three circles at 15,000 ft, 1500 ft and 300 ft. The dropsondes were deployed from the highest altitudes. Further information can be found in the NCAR/NSF C-130 dropsonde README file.

More information on the NCAR GPS Dropsonde system can be found at the following site: http://www.atd.ucar.edu/rtf/facilities/dropsonde (NCAR/EOL, 2006).

2.0.3 NOAA National Weather Service (NOAA/NWS)

This dataset contains 126 high vertical resolution (6 second) soundings from the San Juan, Puerto Rico National Weather Service (NWS) rawinsonde station (WMO 78526) site. The soundins were typically release twice a day (0000 and 1200 UTC).

The detailed description of NWS sounding collection and instrumentation is located in NWS, 1991.

2.0.4 R/V Seward Johnson

This dataset contains 88 high vertical resolution (2 second) soundings from the Research Vessel (R/V) Seward Johnson. Soundings were typically released four to six times per day.

2.1 Detailed Format Description

All upper air soundings were converted to National Center for Atmospheric Research/Earth Observing Laboratory (NCAR/EOL) Sounding Composite Format (ESCF). ESCF is a version of the National Center for Atmospheric Research (NCAR) CLASS format and is an ASCII format consisting of 15 header records for each sounding followed by the data records with associated QC information.

Header Records

The header records (15 total records) contain data type, project ID, site ID, site location, release time, sonde type, meteorological and wind data processors, and the operator's name and comments. The first five header lines contain information identifying the sounding, and have a rigidly defined form. The following 7 header lines are used for auxiliary information and comments about the sounding, and may vary from dataset to dataset. The last 3 header records contain header information for the data columns. Line 13 holds the field names, line 14 the field units, and line 15 contains dashes ('-' characters) delineating the extent of the field.

The five standard header lines are as follows:

LineLabel (padded to 35 char)Contents
1Data Type:Description of type and resolution of data.
2Project ID:Id of weather project.
3Release Site Type/Site ID:Description of the release site.
4Release Location (lon,lat,alt):Position of release site in the format described below.
5UTC Release Time (y,m,d,h,m,s):Time of release, in the format: yyyy, mm, dd, hh:mm:ss

The release location is given as: lon (deg min), lat (deg min), lon (dec. deg), lat (dec. deg), alt (m)

Longitude in deg min is in the format: ddd mm.mm'W where ddd is the number of degrees from True North (with leading zeros if necessary), mm.mm is the decimal number of minutes, and W represents W or E for west or east longitude, respectively. Latitude has the same format as longitude, except there are only two digits for degrees and N or S for north/south latitude. The decimal equivalent of longitude and latitude and station elevation follow.

The seven nonstandard header lines may contain any label and contents. The labels are padded to 35 characters to match the standard header lines. Records for this dataset include the following three nonstandard header lines.

LineLabel (padded to 35 char)Contents
12Nominal Release Time (y,m,d,h,m,s):Nominal time of release, in the format: yyyy, mm, dd, hh:mm:ss

Data Records

The data records each contain time from release, pressure, temperature, dew point, relative humidity, U and V wind components, wind speed and direction, ascent rate, balloon position data, altitude, and quality control flags (see the QC code description). Each data line contains 21 fields, separated by spaces, with a total width of 130 characters. The data are right-justified within the fields. All fields have one decimal place of precision, with the exception of latitude and longitude, which have three decimal places of precision. The contents and sizes of the 21 fields that appear in each data record are as follows:

Field No.WidthFormatParameterUnitsMissing Value
35F5.1Dry-bulb TemperatureDegrees C999.0
45F5.1Dew Point TemperatureDegrees C999.0
55F5.1Relative HumidityPercent999.0
66F6.1U Wind ComponentMeters / Second9999.0
76F6.1V Wind ComponentMeters / Second9999.0
85F5.1Wind SpeedMeters / Second999.0
95F5.1Wind DirectionDegrees999.0
105F5.1Ascent RateMeters / Second999.0
135F5.1Elevation AngleDegrees999.0
145F5.1Azimuth AngleDegrees999.0
164F4.1QC for PressureCode (see below)99.0
174F4.1QC for TemperatureCode (see below)99.0
184F4.1QC for HumidityCode (see below)99.0
194F4.1QC for U Wind ComponentCode (see below)99.0
204F4.1QC for V Wind ComponentCode (see below)99.0
214F4.1QC for Ascension RateCode (see below)99.0

Fields 16 through 21 contain the Quality Control information derived at the NCAR Earth Observing Laboratory (NCAR/EOL). Any QC information from the original sounding is replaced by the following EOL Sounding Composite Format codes:

99.0Unchecked (QC information is "missing.") ("UNCHECKED")
1.0Checked, datum seems physically reasonable. ("GOOD")
2.0Checked, datum seems questionable on physical basis. ("MAYBE")
3.0Checked, datum seems to be in error. ("BAD")
4.0Checked, datum is interpolated. ("ESTIMATED")
9.0Checked, datum was missing in original file. ("MISSING")

Sample Data

The following is a sample record of a RICO high resolution rawinsonde data in EOL Sounding Composite format. The data portion is much longer than 80 characters and, therefore, wraps around to a second line. See section 2.1 for an exact format specification

Data Type:                         High Resolution Sounding
Project ID:                        RICO
Release Site Type/Site ID:         R/V Seward Johnson SWD
Release Location (lon,lat,alt):     74 21.00'W, 21 34.20'N,  -74.35,  21.57,    10.0
UTC Release Time (y,m,d,h,m,s):    2004, 12, 31, 19:34:00
Radiosonde Serial Number:          4117302
Nominal Release Time (y,m,d,h,m,s): 2004, 12, 31, 21:00:00
 Time  Press  Temp  Dewpt  RH    Ucmp   Vcmp   spd   dir   Wcmp     Lon     Lat   Ele   Azim   Alt    Qp   Qt   Qrh  Qu   Qv   QdZ
  sec    mb     C     C     %     m/s    m/s   m/s   deg   m/s      deg     deg   deg   deg     m    code code code code code code
------ ------ ----- ----- ----- ------ ------ ----- ----- ----- -------- ------- ----- ----- ------- ---- ---- ---- ---- ---- ----
   0.0 1019.0  24.6  17.9  66.0   -9.7   -2.2  10.0  77.0 999.0  -74.350  21.570 999.0 999.0    10.0  1.0  1.0  1.0  1.0  1.0  9.0
   2.0 1018.2  24.7  16.7  61.0   -9.3   -6.3  11.2  56.0   3.0 9999.000 999.000 999.0 999.0    16.0  1.0  1.0  1.0  1.0  1.0 99.0
   4.0 1017.4  24.6  16.8  62.0  -10.1   -6.3  11.9  58.0   3.5 9999.000 999.000 999.0 999.0    23.0  1.0  1.0  1.0  1.0  1.0 99.0
   6.0 1016.2  24.5  16.8  62.0  -10.6   -6.4  12.4  59.0   5.0 9999.000 999.000 999.0 999.0    33.0  1.0  1.0  1.0  1.0  1.0 99.0
   8.0 1015.2  24.4  16.9  63.0  -11.2   -6.5  12.9  60.0   4.5 9999.000 999.000 999.0 999.0    42.0  1.0  1.0  1.0  1.0  1.0 99.0
  10.0 1014.2  24.3  17.1  64.0  -11.5   -6.4  13.2  61.0   4.0 9999.000 999.000 999.0 999.0    50.0  1.0  1.0  1.0  1.0  1.0 99.0

2.2 Station List

EOL_GAUS EOL_GAUS     9999.0 999.0 99999.0
N130AR Lockheed C-130     9999.0 999.0 99999.0
SJU San Juan   PR -66.0 18.0 3.0
SWD R/V Seward Johnson SWD   US 9999.0 999.0 99999.0

3.0 Quality Control Processing

This dataset underwent an automated QC process. The dataset underwent internal consistency checks which included two types of checks, gross limit checks on all parameters and rate-of-change checks on temperature, pressure and ascension rate. Some further information on the QC processing conducted by EOL can be found in Loehrer et al. (1996) and Loehrer et al. (1998).

3.1 Gross Limit Checks

These checks were conducted on each sounding and data were automatically flagged as appropriate. Only the data point under examination was flagged. EOL conducted the following gross limit checks on the RICO sounding datasets. In the table P = pressure, T = temperature, RH = relative humidity, U = U wind component, V = V wind component, B = bad, and Q = questionable.

ParameterGross Limit CheckParameter(s) FlaggedFlag Applied
Pressure < 0 mb or > 1050 Mb PB
Altitude < 0 m or > 40000 m P, T, RHQ
Temperature < -90 °C or > 45 °C TT
Dew Point < -99.9 °C or > 33 °C
> Temperature
Relative Humidity < 0% or > 100% RHB
Wind Speed < 0 m/s or > 100 m/s
> 150 m/s
U, V
U, V
U Wind Component < 0 m/s or > 100 m/s
> 150 m/s
V Wind Component < 0 m/s or > 100 m/s
> 150 m/s
Wind Direction < 0° or > 360° U, VB
Ascent Rate < -10 m/s or > 10 m/s P, T, RHQ

3.2 Vertical Consistency Checks

These checks were conducted on each sounding and data were automatically flagged as appropriate. These checks were started at the lowest level of the sounding and compared neighboring 2-sec data points (except at pressures less than 100 Mb where 30-sec average values were used. In the case of checks ensuring that the values increased/decreased as expected, only the data point under examination was flagged. However, for the other checks, all of the data points used in the examination were flagged. All items within the table are as previously defined.

Parameter Vertical Consistency Check Parameter(s) FlaggedFlag Applied
Time decreasing / equal NoneNone
Altitude decreasing / equal P, T, RHQ
Pressure increasing / equal
> 1 Mb/s or < -1 Mb/s
> 2 Mb/s or < -2 Mb/s
P, T, RH
P, T, RH
P, T, RH
Temperature < -15 °C/km
< -30 °C/km
> 50 °C/km (not applied at p < 250 Mb)
>100 °C/km (not applied at p < 250 Mb)
P, T, RH
P, T, RH
P, T, RH
P, T, RH
Ascent Rate change of > 3 m/s or & lt; -3 m/s
change of > 5 m/s or < -5 m/s

3.3 Data Quality Issues


The surface elevations utilized in these files were derived from the Shuttle Radar Topography Mission (SRTM) 90-meter (3 arc-second) resolution elevation data set. The altitudes above the surface were then derived utilizing the hydrostatic equation.

3.3.2 NCAR/NSF C-130 Dropsonde

Data File Specifics

The files contain data calculated at half-second intervals. The variables pressure, temperature, and relative humidity are calibrated values from measurements made by the dropsonde. The dew point is calculated from the relative humidity. The altitude value is calculated from the hydrostatic equation using first available pressure, temperature, and dew point. For the dropsondes specifically, if the sonde is launched over water and transmits data to the surface, the height is calculated by integrating from the surface (sea level) upward. However, if the sonde failed to transmit data to the surface or if the dropsonde is launched over land, because of unknown surface elevations, we integrate from the flight level down. The descent rate of the dropsonde is computed using the time-differentiated hydrostatic equation. The position (lat, lon) are integrated from the initial launch position using the horizontal winds. If a time gap of greater than 60 seconds occurs in the wind measurements, the integration terminates.

All wind data are computed from GPS navigation signals received from the sonde. The raw wind values are calculated at a one half second data rate by a commercial processing card. The resolution of the data is half second.

3.3.3 NOAA National Weather Service (NOAA/NWS)

NWS soundings during RICO utilized either the VIZ type radiosonde produced by Sippican Inc. (http://www.sippican.com/meteorological.html) or the Vaisala RS-80 radiosonde (http://www.vaisala.com, under "Product and Systems/Radiosonde,dropsondes,rockectsondes").

The use of the raw 6-sec resolution elevation and azimuth angle data to derive the winds sometimes led to large oscillations in wind speed, due to the presence of oscillations in the elevation angle data, particularly at low elevation angles. The general approach to correct this problem was to remove the out-lier radiosonde position data before computing the wind components (Williams et al. 1993). For both the azimuth and elevation angles from 360 sec to the end of the sounding, a ninth order polynomial was fit to the curve. The residuals were calculated and compared to the observed values. The outliers of the residuals were then removed.

Then to help correct the more extensive problems at low elevation angles within 10 degrees of the limiting angles (LA) some additional smoothing was applied. If the elevation angle was between (LA + 7.5) and (LA + 10), the new elevation angle was computed with a 2 min linear fit. If the elevation angle was between (LA + 5) and (LA + 7.5), the new elevation angle was computed with a 3 min linear fit. If the elevation angle was less than (LA + 5), the new elevation angle was computed with a 4 min linear fit. If the number of observations with low elevation angles was greater than 20% of the total number of observations for the sounding no frequency smoothing occurred.

Then, for the elevation angle only, a finite Fourier analysis was performed on the residuals. Periods from 90-190 sec were removed and those below 30 sec were flattened.

Finally, a 2 min second order polynomial was then fit to the position to derive the u and v wind components, except for the beginning and end minute (or 1.5 minutes if over 50 Mb) which used a 3 min fit. If there were less than 15% of the total number of points, not counting the beginning or end of the flight, on one side of the point for which the wind value was being computed, a linear fit was used.

For further information on this methodology and its changes since Williams et al. (1993) please see Williams, et al. (1998).

Near Surface Winds

A common problem in near surface wind speed values calculated from the 6-second position data is that the first radiosonde wind speed is much higher than the independently measured surface value. The calculated radiosonde winds then decrease rapidly so that within about 60 s (20-30 Mb) after release the wind speeds are more realistic. The cause of this appears to be the acceptance of radiosonde position data prior to a "good lock" being achieved on the radiosonde by the tracking system. Thus there appear to be rapid positional shifts of the radiosonde while the tracking system "searches" for the radiosonde.

Wind Oscillations

Despite the extensive efforts to remove oscillations in wind speeds caused by oscillations in elevation angles (see Section 2.2) there are occasional cases with remaining oscillations. Most of the remaining oscillations have periods just slightly longer than the 190 s maximum point of our notch filter.

3.3.4 R/V Seward Johnson

Visual Checks

All soundings were examined via skewt plots and where needed changes were made to the flags applied by the automated checks listed above.

4.0 References

Loehrer, S. M., T. A. Edmands, and J. A. Moore, 1996: TOGA COARE upper-air sounding data archive: development and quality control procedures. Bull. Amer. Meteor. Soc., 77, 2651-2671.

Loehrer, S. M., S. F. Williams, and J. A. Moore, 1998: Results from UCAR/JOSS quality control of atmospheric soundings from field projects. Preprints, Tenth Symposium on Meteorological Observations and Instrumentation, Phoenix, AZ, Amer. Meteor. Soc., 1-6.

NCAR/EOL, cited 2006: GPS Dropsonde [Available online from http://www.atd.ucar.edu/rtf/facilities/dropsonde].

NWS, 1991: Micro-ART Observation and Rework Programs Technical Document, National Weather Service, National Oceanic and Atmospheric Administration, Washington, D.C., March 1991.

Wade, C. G., 1995: Calibration and data reduction problems affecting National Weather Service radiosonde humidity measurements. Preprints, Ninth Symposium on Meteorological Observations and Instrumentation, Charlotte, NC, Amer. Meteor. Soc., 37-42.

Williams, S. F., C. G. Wade, and C. Morel, 1993: A comparison of high resolution radiosonde winds: 6-second Micro-ART winds versus 10-second CLASS LORAN winds. Preprints, Eighth Symposium on Meteorological Observations and Instrumentation, Anaheim, California, Amer. Meteor. Soc., 60-65.

Williams, S. F., S. M. Loehrer, and D. R. Gallant, 1998: Computation of high-resolution National Weather Service rawinsonde winds. Preprints, Tenth Symposium on Meteorological Observations and Instrumentation, Phoenix, AZ, Amer. Meteor. Soc., 387-391.