Solar radiation drives many environmental processes; however, needs to be estimated indirectly from more commonly measured meteorological variables since these data are not readily available from most climate stations. A geo-referenced dataset from 159 climate stations of 17 variables (maximum possible sunshine duration, mean, minimum and maximum air temperature, soil temperature, mean and maximum relative humidity, precipitation, cloudiness, evapotranspiration, extraterrestrial radiation, day length, declination angle, day of the year, latitude, longitude, and altitude) was used to model spatio-temporal dynamics of solar radiation over Turkey. A total of 78 empirical models of different mathematical functions with a different combination of 17 explanatory variables were compared based on the error statistics of the Jackknifing validation. The empirical models had adjusted coefficient of determination (R adj 2 ) values of 22.7–96.5% based on the parameterization dataset (P 0.05), while the performance of the hybrid models was worse than that of the linear, quadratic and cubic models in terms of maximum relative percentage error (e) (P < 0.01). In comparing the interpolation methods of inverse distance weighting and universal co-kriging, anisotropic spherical semi-variogram model of universal co-kriging was found to provide the best description of spatial autocorrelation and variability latent in these data based on the spatial leave-one-out cross-validation.