P-spline smoothing for spatial data collected worldwide