In order to perform continuous simulation in GSSHA hydrometeorological data are required for the entire period of the continuous simulation. The required data are hourly values of:
- barometric pressure,
- relative humidity,
- total sky cover,
- wind speed,
- dry bulb temperature,
- direct radiation, and
- global radiation.
These data are available from a wide variety of sources and three different input file formats can be used. The GSSHA model also contains provisions to synthetically produce solar radiation estimates and fill in short term data gaps.
9.3.1 Hydrometeorological Input File Formats
The three file formats that GSSHA supports are discussed below in order of preference. Simple file formats are preferred because the likelihood of a data format error is reduced. The WES format is the simplest because it uses no alphabetical characters. Users are encouraged to produce and use hydrometeorological data files in the WES format for continuous simulations. The other formats contain data quality flags, extraneous data, and occasionally errors. Such errors can be very difficult to locate. Reducing data down to the WES format through independent software (such as a spreadsheet) ensures that the data are properly formatted and contains no extraneous values.
126.96.36.199 HMET_WES (Recommended format)
This file contains only hourly values of the required hydrometeorological parameters: barometric pressure, relative humidity, total sky cover, wind speed, temperature, direct radiation, and global radiation. The HMET_WES format contains the following numbers in columns 1-11 (Table 15):
|Col.||Variable||Units||No Data Flags||Type|
|1||Year (4 digit)||integer|
|5||Barometric Pressure||in Hg||No Data (ND) = 99.999||real|
|7||Total Sky Cover||%||ND=999||integer|
|9||Dry Bulb Temperature||°F||ND=999||integer|
|10||Direct Radiation||W h m-2||ND=9999.99||real|
|11||Global Radiation||W h m-2||ND=9999.99||real|
Table 15 – HMET_WES file format
Data format is important; the year is stored as a 4 digit number; the hour varies between 0 and 23. The data in columns 1-4, 6-9 are stored as integers, while the data in columns 5, 10, and 11 are stored as real numbers. The data in this file should be space delimited. The number of spaces between columns does not matter, provided that the total width of the line is less than 256 characters. There must be 1 hourly record for each hour of the intended duration of the continuous simulation, even if there are no data available for that hour. Gaps in the hourly record will cause the GSSHA simulation to terminate prematurely.
The following are three example lines from a HMET_WES formatted hydrometeorological data file (note the distinct presence of both real and integer numbers):
1995 4 14 23 29.625 86 100 000 53 9999.99 9999.99
1995 4 15 0 29.625 83 100 003 53 9999.99 9999.99
1995 4 15 1 29.615 90 80 003 52 9999.99 9999.99
This data input format is strongly recommended because it is compact, easy to read and write and contains no extraneous information. The no data flags are important. GSSHA has a set procedure for dealing with missing data.
188.8.131.52 HMET_SAMSON (Recommended over surface airways format)
In 1993 the National Climatic Data Center (NCDC) published a CD-ROM containing surface meteorological observations from 1961-1990. This CD-ROM is the result of collaboration between the National Oceanic and Atmospheric Administration (NOAA) and the U.S. Department of Energy, National Renewable Energy Laboratory. The SAMSON CD-ROM represents an excellent nationwide dataset for the 1961-1990 period. Unfortunately, this CD-ROM is not being updated as new data become available.
The following is an example of data from the SAMSON CD-ROM. A reader program provided with the SAMSON data produces this data format. The ASCII Samson data format includes a one line station header (in this case the data are from the Memphis, Tennessee station), and a one line descriptive header with column ID’s. These two header lines are followed by a user selectable number of hourly data points (in this example two hourly lines are shown).
~13893 MEMPHIS TN -6 N35 03 W089 59 87
~YR MO DA HR I 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
82 4 1 1 0 0 0 0 ?0 0 ?0 0 ?0 0 0 12.8 6.1 64 1011 100 2.1 24.1 77777 999999999 16 99999. 0 26
82 4 1 2 0 0 0 0 ?0 0 ?0 0 ?0 0 0 12.2 6.1 67 1011 80 2.1 24.1 77777 999999999 16 99999. 0 26
The first header line consist of the Weather Bureau Army Navy (WBAN) number for the station (13893), followed by 26 character descriptive string (MEMPHIS), the state name where the station is located (TN), the time difference between local standard and UTC in hours (-6), degrees and minutes of latitude, N|S, (N35 03), degrees and minutes of longitude ,E|W, (W089 59) and the station elevation in meters above sea level (87). The first header record containing the station information was written by NCDC with the following FORTRAN format statement:
The second line is for the identifier record, which shows the numerical identifier for the data elements in the National Solar Radiation Data Base (NSRDB) Synoptic Format (see: http://rredc.nrel.gov/solar/pubs/NSRDB/).
The observed data occupy the third record and all ensuing records (lines) in the file until a new year of data occurs. When a new year of data is encountered the header record and identifier record are repeated, followed by more data records. The data records include year, month, day, hour (1-24) in local standard time. The number below the “I” column is an “observation indicator” for which 0 means the weather observation was made, and a 9 indicates that the observation is missing. The following table describes the contents of each data record (Table 16):
|No.||Name||Value Range||Description and Notes:|
|YR||Year||(61-90)||Year of observation|
|MO||Month||(1-12)||Month of observation|
|DA||Day||(1-31)||Day of month|
|HR||Hour||(1-24)||Hour of the day, local standard time.|
|I||Observation Indicator|| 0
| 0= weather observation made|
9= no observation made
NOTE: if this field=9 OR if field 13 (wind speed) = missing (9999. or 99.0) then fields 6,7,8,10,11,17, and 18 were all modeled and not actually observed.
|1||Extraterrestrial Horizontal Radiation||0-1415||Amount of solar radiation in W-h m-2 on a horizontal surface at the top of the atmosphere during the 60 minutes preceding the hour indicated.|
|2||Extraterrestrial Direct Normal Radiation||0-1415||Amount of solar radiation in W-h m-2 received on a surface normal to the sun at the top of the atmosphere during the 60 minutes preceding the hour indicated.|
|3||Global Horizontal Radiation|| 0-1415
flag for data 0-9
|Total amount of direct and diffuse solar radiation in W-h m-2 on a horizontal surface during the 60 minutes preceding the hour indicated.|
|4||Direct Normal Radiation|| 0-1415
flag for data 0-9
|Amount of solar radiation in W-h m-2 received within 5.7 degrees field of view centered on the sun during the 60 minutes preceding the hour indicated.|
|5||Diffuse Horizontal Radiation|| 0-1415
flag for data 0-9
|Amount of solar radiation in W-h m-2 received from the sky (excluding the solar disk) on a horizontal surface, during the 60 minutes preceding the hour indicated.|
|6||Total Sky Cover|| 1-10
|Amount of sky dome (in tenths) covered by clouds|
|7||Opaque Sky Cover|| 1-10
|Amount of sky dome (in tenths) covered by clouds that prevent observing the sky or higher cloud layers.|
|8||Dry Bulb Temperature|| -70.0 to +60.0
|Dry bulb air temperature in degrees C.|
|9||Dew point|| -70.0 to +60.0
|Dew point temperature in degrees C.|
|10||Relative Humidity|| 0-100
|Relative humidity in percent.|
|11||Station pressure|| 700-1100
|Station barometric pressure in mb.|
|12||Wind Direction|| 0-360
|Wind direction in degrees. N=0 or 360, E=90, S=180, W=270|
|13||Wind Speed|| 0.0-99.0
9999. or 99.0= missing data
|Wind speed in m s-1|
|Horizontal visibility in kilometers.|
|15||Ceiling height|| 0-30450
|Ceiling height in meters|
|16||Present weather||Table denoted by 9 indicators.||Present weather conditions. (See NOAA document for codes).|
|17||Precipitable water|| 0-100
|Precipitable water in millimeters|
|18||Broadband aerosol optical depth|| 0.0-0.900
|Broadband aerosol optical depth (broadband turbidity) on the day indicated.|
|19||Snow depth|| 1-100
|Snow depth in centimeters on the day indicated.|
|20||Days since last snowfall|| 0-88
88=88 or > days
|Number of days since last snowfall|
|21||Hourly precipitation||000000-099999||In inches and hundredths.|
Table 16 – HMET_SAMPSON file format
- The SAMSON CD-ROM is available by contacting the NCDC at:
- Climatic Services Branch
- National Climatic Data Center
- Federal Building
- Asheville, NC 28801-2696
- Phone - (704) 271-4800
- Fax - (704) 271-4876
184.108.40.206 HMET_SUFAWAYS - NOAA/NCDC Surface Airways Format (Not recommended)
The Surface Airways data format is included in the list of supported formats because surface meteorological data ordered from NCDC are delivered in this format. In the Surface Airways data format each data element (wind, temperature, radiation, etc.) is written on a separate line in the file. Each line contains 24 hourly values. Missing data are not flagged.
The details of this data format are discussed in “TD-3280, Hourly Surface Airways Observations”, published by the National Oceanic and Atmospheric Administration, National Climatic Data Center. GSSHA users are strongly encouraged to reformat Surface Airways observations into the WES data format. There is a strong likelihood that data in the Surface Airways format will produce problems during the GSSHA simulation. This occurs not as consequence of faulty data, but because of the large number of data flags in the Surface Airways output format make the Surface Airways data format difficult to read in a robust fashion.
220.127.116.11 HMET_ASCII Gridded Data (Recommended format)
(Version 7 and higher) This format is for data that comes from other models or is otherwise pre-processed into gridded data. The grid shape and size of the hmet data does not need to match the GSSHA grid (and usually does not), but it does need to match itself (be consistent across data types and over time.)
The data are in individual grid files for each data type and each time step. They are usually placed in a folder by themselves as they will be quite numerous. They have a specific naming convention that GSSHA relies on: YYYYMMDDHH_TYPE.asc
The data for the ASCII Gridded format follows the HMET_WES list of data.
The type flags and units for each class of data are as follows:
|Data||File name code||Example file name||Units|
|Cloud Cover||Clod||2023010512_Clod.asc||0-100, Fraction of the sky obscured by clouds|
|Atmospheric Pressure||Pres||2023010512_Pres.asc||in Hg, Atmospheric pressure near the ground surface|
|Relative Humidity||RlHm||2023010512_RlHm.asc||0-100, Relative humidity|
|Temperature||Temp||2023010512_Temp.asc||deg F, air temperature near the ground surface|
|Wind Speed||WndS||2023010512_WndS.asc||kts, wind speed near the ground surface|
|Direct Radiation||Drad||2023010512_Drad.asc||W h m-2, Direct Radiation|
|Global Radiation||Grad||2023010512_Grad.asc||W h m-2, Global Radiation|
The hmet file format is straightforward: it is simply the list of paths and date/time stamps for each hour. GSSHA will take the path and date/time info and append the underscore and type name to open the individual data files.
In this example, the files in the HMET_DATA folder would look like this:
To use gridded ascii data in your model, in the project file use the card HMET_ASCII "filename.hmet".
9.2.2 Missing Hydrometeorological Data
Frequently, hydrometeorological records have periods of missing data because observational instruments are inoperative, the data are not reliable, or certain observations (e.g. direct radiation, global radiation) are not available for a given station. Because hourly hydrometeorological inputs are required for the entire simulation period, missing data are synthetically produced within the GSSHA model.
18.104.22.168 Radiation Data
Many hydrometeorological stations do not record solar radiation measurements. Missing solar radiation data are simulated in GSSHA based on the location of the watershed, day of the year, and time of day, as described in Section 9.1. When only a portion of the solar radiation data is missing, it is replaced as described in the next section.
22.214.171.124 Short Duration Gaps in Hydrometeorological Measurements
When required data are missing, GSSHA creates the needed data by either persistence or persistence adjusted for the time of day. Missing variables without strong diurnal variability are filled in by “persistence” estimates. These variables include:
- barometric pressure,
- wind speed,
- dew point temperature and,
- percent cloud cover.
Persistence means that the value is held constant until a new observation becomes available. Missing variables with strong diurnal components:
- air temperature,
- global radiation,
- and direct radiation.
are replaced with the last good reading from the same time of day. It is the responsibility of the user to check for missing data. Extended periods of missing record (many days) will cause GSSHA to simulate conditions for the last day of good record over and over again. This scheme requires that the hydrometeorological data file begin with at least one day of good observations of all required variables.
- 9 Continuous
- 9.1 Computation of Evaporation and Evapo-transpiration
- 9.2 Computation of Soil Moisture
- 9.3 Hydrometeorological Data
- 9.4 Snowfall Accumulation and Melting
- 9.5 Sequence of Events During Long-Term Simulations