Temporal solar variability significantly affects the integration of solar power systems into the grid. It is thus essential to predict temporal solar variability, particularly given the increasing popularity of solar power generation globally. In this paper, the daily variability of solar irradiance at four sites across Australia is quantified using observed time series of global horizontal irradiance for 2003-2012. It is shown that the daily variability strongly depends on sky clearness with generally low values under a clear or overcast condition and high values under an intermittent cloudiness condition. Various statistical techniques are adopted to model the daily variability using meteorological variables selected from the ERA-Interim reanalysis as predictors. The nonlinear regression technique (i.e. random forest) is demonstrated to perform the best while the performance of the simple analog method is only slightly worse. Among the four sites, Alice Springs has the lowest daily variability index on average and Rockhampton has the highest daily variability index on average. The modelling results of the four sites produced by random forest have a correlation coefficient of above 0.7 and a median relative error around 40%. While the approach of statistical downscaling from a large spatial domain has been applied for other problems, it is shown in this study that it generally suffices to use only the predictors at a single near point for the problem of solar variability. The relative importance of the involved meteorological variables and the effects of clearness on the modelling of the daily variability are also explored. Crown Copyright (C) 2014 Published by Elsevier Ltd. All rights reserved.