Download
1 / 33
330 likes | 337 Views
This presentation presents a reliable method for estimating snow water equivalent (SWE) at ad-hoc locations in Alaska using ArcMap's Geostatistical Analyst extension. The method utilizes a SWE estimation model based on the Helsinki University of Technology's snow emission model. The model is applied to 30 years of European Space Agency satellite observations covering Alaska. The SWE estimates are based on data from the Global Snow Monitoring Network and are stored in NETCDF files by the National Snow and Ice Data Center. The estimation technique uses Kriging to predict SWE values at satellite sample points.
E N D
Estimating Alaska Snow Loads By Russ Frith University of Alaska at Anchorage Civil Engineering Department
Problem Statement • Snow loads are an important aspect of structural design in cold climates. In some areas, the weight of accumulated snow on a structure represents the most significant load that the structure will endure in its lifetime. • In Alaska, snow load data are calculated from snow water equivalent (SWE) measurements at weather/climate stations. • There are relatively few climate stations in Alaska that reliably report SWE measurements. As a result, SWE measurements or estimates are unknown for most of the land area in Alaska. • The presentation presents a method for reliably estimating SWE at ad-hoc locations in Alaska. The estimation technique uses ArcMap’s Geostatistical Analyst extension.
Front Matter : Model • A SWE estimation model was selected. Snow depth was estimated using the Helsinki University of Technology snow emission model-based iterative inversion algorithm. • Pullianen, J.T., Grandell, J. & Hallikainen, M. (1999). HUT snow emission model and its applicability to snow water equivalent retrieval. IEEE Transactions on Geoscience and Remote Sensing 37(3), 1378-1390 • Results of this model are readily available in NETCDF files from the National Snow and Ice Data Center. This model was applied to 30 years of European Space Agency satellite observations covering Alaska. The data produced from this model have rigorous quality control checks. • This is one of the few models applied to physical observations from a long period of record and which outputs are publicly available. • Other comparable models exist but many of these models’ outputs are not publicly available or are not processed for such a lengthy period of record.
Front Matter : Data • Global Snow Monitoring Network makes daily satellite observations. The observational data is fed into the SWE model which then estimates SWE for given latitude and longitude on a given grid that covers Alaska’s land area. • SWE estimates from the model are bundled into NETCDF. • The NETCDF files are held by the National Snow and Ice Data Center, www.nsidc.org • The NETCDF locations are converted into ArcInfo shape point shape files. • NETCDF (Network Common Data Form) is a product of Unidata Corporation is a set of software libraries and self-describing, machine-independent data formats that support creation, access, and sharing of array-oriented scientific data.
Front Matter : Data • Each green dot represents a satellite location for a modeled SWE estimate for a 25 km grid cell starting around 1981. • Each blue dot represent a Global Historical Climate Network SWE monitoring site. • At each location, a 50 year mean recurrence interval (MRI) value of SWE was calculated. • The values were stored in a geodatabase.
Use Kriging to Estimate 50 Year-MRI Predicting SWE at Satellite Sample Points 1. Kriging Fits function to • Specified number of points, • All points within specified radius. 2. Kriging is based on the idea that one can make inferences regarding a random function Z(x), given data points Z(x1), Z(x2), …, Z(xn). 3. The basis of this technique is the rate at which the variance between points changes over space. 4. The variance is expressed in a semivariogram which shows how the average difference between values at points changes with distance between points. 5. Z(x) = m(x) + g(h) + e • m(x) is the constant mean • g(h) is a random spatially correlated component • e is the residual error
Variogram The amount and form of spatial autocorrelation can be described by a semivariogram that shows how differences in values increase with geographical separation. Semi-Variogram A semi-variogram is a function that relates semi-variance (or dissimilarity) of data points to the distance that separates them.
Calculating the Semi-Variogram • To compute a variogram, need to determine how variance increases with distance. • Begin be dividing the range of distance into a set of discrete intervals. • For every pair of points, compute distance and the squared difference in z values. • Assign each pair to one of the distance ranges, and accumulate the total variance in each range. • After every pair has been used (or a sample of pairs in a large dataset) compute the average variance in each distance range. • Plot this value at the midpoint distance of each range.
Scatter Plots For each h-scatterplot, a value γ is calculated. These values are used to fit a variogram model.
Kriging Example Estimate a value for point 0 (65E, 137N), based on the seven surrounding sample points. The table indicates the (x,y) coordinates of the seven sample points, their corresponding values of V, the variable of interest and their distances from point 0.
Spatial Continuity Parameters
First, find the distance matrix Covariances are calculated based on the distance between points using the model:
ArcGIS & Kriging Using the Statistical Analyst Extension
Recall that the only places where SWE (i.e. snow load) are measured directly are at GHCN stations. GHCN SWE Monitoring Site The green dots are locations where SWE is estimated by the model. How good are these estimates? Using Kriging, predict the GHCN measurements from the model estimates.
Predicting Measured Values from Modeled Calculations – Results!
