This paper presents the development, testing and validation of a novel generic type universal model consisting of a set of sine and cosine harmonics in the temporal and spatial domain suitably parameterized for the prediction of the mean expected global solar radiation H(n,phi) on the horizontal for a day, n, at any latitude phi. Its prediction power is further enhanced with the introduction of a correction term for the site altitude taking into account the phi dependent atmospheric height. Solar radiation data from 53 stations around the earth were obtained from GEBA database to train the model. H(n,phi) is expressed by a Fourier series of compact form with the zero frequency component dependent on phi providing the main spatial dependence and two n dependent harmonics in the form of cosine functions giving the time dependence. The phi dependent model parameters follow symmetry rules and are expressed by Fourier series up to the 3rd order harmonic. The 3D spatiotemporal profile of the model is in agreement to the extraterrestrial one. The model was validated using GEBA data from additional 28 sites and compared with NASA, PVGIS and SoDa data, showing the robustness, reliability and prediction accuracy of the proposed model. (C) 2018 Elsevier Ltd. All rights reserved.