Exploring Interpolation Methods (Kriging, IDW, etc) on 30-day NDAWN Temperature Datasets in ArcGIS Pro

Abstract

Temperature data is dynamic and can be published frequently. Tools that can retrieve temperature data can be utilized on the fly to update and add-on to datasets. This type data is also collected at stations like NDAWN stations and due to the nature of that, is missing values where stations are not located. Interpolation can be used to solve this problem by creating and inferring the values in between based on the existing sample values collected. There are different methods of interpolation such as IDW, various Kriging methods, and RBF. There are differences such as exact and non-exact consideration of sample values between these different methods. They are all able to produce an interpolated surface. Overall, there is no best method, but in some studies, various types of Kriging methods have produced the “best” interpolated surfaces with considerations to the scenario variables.


Problem Statement

Temperature data is dynamic and changes over time. NDAWN provides temperature data on the web for different temporal resolutions. To examine data on a monthly frequency, 30 days of temperature data must be extracted on the fly. Temperature data like maximum and minimum daily temperatures can also be extracted on the fly to examine the high and low temperatures. The high temperature is defined in this lab as the largest value of the maximum daily temperatures recorded for the 30-day period. The low temperature is defined as the lowest value of the minimum daily temperature recorded for the 30-day period. Furthermore, NDAWN stations are dispersed unevenly. Temperatures in areas between station collection sites must be inferred through interpolation, which can be done with different methods.


Results

There was a total of 131 NDAWN stations. The daily temperatures collected by the NDAWN stations varied. The minimum daily temperatures for all stations had a range from -40.864 to 34.704 and a maximum temperature range of -21.388 to 62.24. The average temperature for the 30-day period for all stations ranged from -0.37 to 11.72 F (Fig. 2).

The interpolated results for the 30-day period high and low temperatures from Feb 01, 2021 were different for each algorithm used. Each output raster was symbolized using natural break classification with 10 classes. The IDW low interpolated results ranged from -40.86 to -24.41 and the high interpolated results ranged from 38.23 to 62.23. The kriging¬¬¬ low interpolated results ranged from -36.23 to -25.97, and the high interpolated results ranged from 43.43 to 58.91. The EBK low interpolated results ranged from -42.18 to -24.88, and the high interpolated results ranged from 38.46 to 63.53. The RBF low interpolated results ranged from -40.89 to -24.41, and the high interpolated results ranged from 38.06 to 63.27 (Fig 3).

The highest maximum daily temperatures appears to be recorded in the south west portion of North Dakota and the lowest minimum daily temperatures appears to be in areas outside of the red river valley of the north river. In terms of smoothness, EBK produced the most “smooth” interpolated surface and Kriging had the least smooth with edges/straight line sin the interpolated results (Fig.3).


References

  • CRONQVIST, F. (2018). Interpolation of temperature data for improved weather forecasts. DEGREE PROJECT IN MATHEMATICS, SECOND CYCLE, 1–98. https://kth.diva-portal.org/smash/get/diva2:1214240/FULLTEXT01.pdf

  • ESRI. (n.d.-a). How inverse distance weighted interpolation works—ArcGIS Pro | Documentation. ArcGIS Pro. Retrieved April 21, 2021, from https://pro.arcgis.com/en/pro-app/latest/help/analysis/geostatistical-analyst/how-inverse-distance-weighted-interpolation-works.htm

  • ESRI. (n.d.-b). How radial basis functions work—ArcGIS Pro | Documentation. ArcGIS Pro. Retrieved April 21, 2021, from https://pro.arcgis.com/en/pro-app/latest/help/analysis/geostatistical-analyst/how-radial-basis-functions-work.htm#:%7E:text=Radial%20basis%20functions%20(RBFs)%20are,Spline%20with%20tension

  • ESRI. (n.d.-c). Understanding ordinary kriging—ArcMap | Documentation. ArcGIS Desktop. Retrieved April 21, 2021, from https://desktop.arcgis.com/en/arcmap/latest/extensions/geostatistical-analyst/understanding-ordinary-kriging.htm

  • ESRI. (n.d.-d). What is empirical Bayesian kriging?—ArcGIS Pro | Documentation. ArcGIS Pro. Retrieved April 21, 2021, from https://pro.arcgis.com/en/pro-app/latest/help/analysis/geostatistical-analyst/what-is-empirical-bayesian-kriging-.htm

  • Hofstra, N., Haylock, M., New, M., Jones, P., & Frei, C. (2008). Comparison of six methods for the interpolation of daily, European climate data. Journal of Geophysical Research, 113(D21), 1–19. https://doi.org/10.1029/2008jd010100