Background
I'm posting just to try to start a discussion about uncertainty
fields in the URMA. We're still using the Obs database, but look
forward to switching to URMA since it is much less work for us to
maintain. Many forecasters in the office object to the URMA when they
see the big differences between the Obs database and the URMA. And
then if they find the uncertainty grids in URMA, they balk at it. How
can temperature uncertainty be 2-3 degrees (as it ALWAYS is) when the
observation and grid differ by 10 degrees at times? Of course I agree
completely with that point. But I also point out that we have no idea
what the uncertainty of the Obs database analysis is, since it doesn't
have this field at all.
I think good uncertainty fields would be a big help in selling the
URMA to the rest of the forecast staff, and I think it would be more
useful for doing grid verification. For example, the current defacto
standard for temperature forecasts in BOIVerify is to be within 3
degrees, but in reality, several analysis grid points may have much
larger uncertainty. I would argue that a forecast error near to or
less than the analysis uncertainty is a hit, even if that error is 7
degrees. Along those lines, I started thinking about uncertainty in
our grids, and I came up with a few sources of uncertainty.
- distance from observations that affect that point
- sensor accuracy, roughly +-2 degrees for ASOS temperatures I think
- terrain variation within a grid cell of the analysis
As I've looked in depth at terrain variation in my CWA (MSO) before,
I thought I would take a closer look at that aspect of uncertainty alone.
Methods
First I downloaded the 30m elevation data from the USGS, then
converted it to a text format that I could read in python. I wrote a
procedure to read the elevation for each 30m topo data point and
assign it to one of my GFE grid points. From that I was able to
calculate the maximum, minimum, range, average, and standard deviation
of hi-res topography values that correspond to each GFE grid point in
my CWA.
There was roughly 8,400 30m elevation points per 2.5x2.5 km grid
point in GFE. It wasn't the same for each cell due to differences in
the map projections. It took well over 4-hours to ingest and analyze
all the data points in GFE.
Next, I assumed an average lapse rate (-3.5F / 1000ft), well mixed
lapse rate (-5.5F / 1000ft), and strong inversion (+10F / 1000ft) to
calculate the variation of temperatures across a grid point just due
to the varied elevation across it.
Results
Topograhpy
The first image below is the range (max-min) of elevations for each
grid point in the MSO CWA. The next image (using the same colorscale)
is the standard deviation of elevation values for each grid point.
Above the graphic shows percent of grid points with an elevation
range greater than or equal to the x-axis values. So 79% of grid cells
in the MSO CWA (not the whole image, just the MSO CWA portion) have
cover an area with an elevation range greater than 1,000 feet! Over a
third of the area has grid cells with an elevation range over 2,000 feet.
Average lapse rate
Assuming an average lapse rate of -3.5F per 1,000 feet, the range
(and standard deviation) of temperature values by grid cell is shown
below. The pixels in yellow have a range 5-10 degrees of variation
across a grid cell, but note that the standard deviation is much less.
Gray areas are less than 3 degrees, and green is 3-5 degrees of variation.

Range |

Standard Deviation |
Well mixed, dry adiabatic lapse rate
Assuming a lapse rate of -5.5F per 1,000 feet, the range (and
standard deviation) of temperature values by grid cell is shown below.
The pixels in yellow have a range 5-10 degrees of variation across a
grid cell, and the red pixels have more than 10 degrees of variation
across the grid cell.

Range |

Standard Deviation |
Strong inversion
Assuming a lapse rate of +10F per 1,000 feet, the range (and standard
deviation) of temperature values by grid cell is shown below. Strong
inversions like this frequently occur in mountain cold pool areas of
the west, especially in the winter.

Range |

Standard Deviation |
Discussion
So this only looked at one aspect of uncertainty, adding other
sources as mentioned above would make the uncertainty greater, plus
considering the uncertainty associated with the lapse rates used in
any such analysis. This was a pretty simplistic analysis, I'm hoping
someone with much better statistics can comment as well.
The uncertainty associated with the terrain variations across a grid
point could be considered a source of uncertainty that is intrinsic to
the analysis grid. In this case, I wonder if this could be utilized
during the analysis process as one consideration for how strongly to
weight an observation compared to the background. For instance, most
of our ASOS stations are associated with grid points that have
elevation ranges such that the standard deviation of temperature
across the grid cell for all cases above is 0-2 degrees. I would hope
that this, combined with the quality of the ASOS data, would cause
them to be weighted very heavily. On the other hand, SNOTEL and RAWS
weather stations are often located on slopes or in very hilly areas
and might be weighted less.
Thanks for reading,
-Ryan