Spherical Triangulation

After 5 days on Heron Island and a narrow escape from Cyclone Debbie, I am now back to spherical gridding!

The most elegant method for spatial integration of irregular temperature data must surely be spherical triangulation over  the earth’s surface. This is because it treats each measurement equally by covering the earth’s surface with a triangular mesh of station & SST nodes. Unlike linear triangulation (described previously), spherical triangulation also spans polar regions. There is no need for any ‘kriging’ or linear interpolation into sparse polar regions, since they are naturally included. I finally deciphered IDL’s spherical triangulation output,  thanks to Nick Stokes.  Here then is the result in 3D for temperatures in January 2016. The shading for each triangle is the average of each node’s temperature anomaly (-5C – +5C).

Figure 1: Triangular mesh with temperature shading showing US stations in fine mesh. Transparency of 50% also shows ghosted the grid on the hidden hemisphere. Image currently takes 72 hours to render, so it is a little difficult to optimise!

I now need to find a better visualisation method as this one takes way too long, however at least it shows how triangulation now covers both poles.

Each triangular area shown in Figure 1 is calculated in 3D cartesian coordinates to derive the area weighting used for averaging. Figure 2 shows  how the final spherical results compare to the 2D (lat,lon) triangulation results.

Figure 2: Comparison of spherical triangulation temperature results with 2-D triangulation. Main difference is the handling of polar regions. In blue is shown  Cowtan & Way who use kriging to extrapolate data into polar regions

There is really very little difference in the annual temperature anomalies between the spherical results and the 2-D triangulation results. Based on these results it would seem that Cowtan & Way have exaggerated polar warming effects between 2005 and 2013. Figure 3 shows the monthly comparison and just how remarkably similar the 2-D and 3-D results are despite completely independent methods of integration.

Figure 3: Comparison of global monthly tempearture anomalies calculated using spherical triangulation and 2-d (lat,lon) triangulation. The difference proves that there are only small differences in final result.

Aesthetically the spherical triangulation grid is my favourite. Unfortunately, the extra effort makes only a tiny difference compared to the more easy 2D triangulation solution based on lat,lon coordinates. Despite this, both methods, in my opinion, are better than simple rectangular gridding as used by Hadcrut4 and (partially)  by GISS and NOAA. Furthermore they avoid interpolation.

I will post the code and the data soon.

This entry was posted in climate science, Science and tagged . Bookmark the permalink.

2 Responses to Spherical Triangulation

  1. Olof R says:

    Nice work Clive,

    I guess that the latest version is virtually the same as Nick Stokes TempLSmesh, but with adjusted GHCNv3 and HADSST3 instead.

    I think that the main difference with C&W is the use of GHCN v3 adjusted, which has a known Arctic cooling bias unlike CRUTEM4 land data. (the PHA distrust the rapidly warming arctic stations and adjust them down)

    To test this, and see if the fit with C&W improves you could use:
    – unadjusted GHCNv3
    – adjusted GHCNv4 (which hasn’t the cooling bias, because more arctic stations convince the PHA that the warming is real)
    – HadCRUT4 gridded data (probably the ultimate apples-to-apples test of spherical triangulation vs C&W kriging)

    • Clive Best says:


      Yes I think this version in effect uses the same algorithm as Nick Stokes’s TempLSmesh, but I didn’t know that beforehand, and in addition we are using completely different software as well.

      My plan is to eventually repeat this using the latest CRUTEM4 station data and HadSST3. That would then be an exact comparison with Cowtan & Way as you say.

Leave a Reply