Numerical Weather Prediction (NWP) models over limited areas enable the simulation of local atmospheric process in more detail and with a higher degree of accuracy when compared to global models. Limited-area NWP models can out-perform their global counterparts due to higher resolution (ability to explicitly simulate processes) and tailored physics (global models, unless run as a physics ensemble, have one set of parameterisation schemes for the whole globe). However, increased accuracy from an NWP model is not guaranteed and can vary based on location and variable of interest. In this paper we present a method for combining the output of a limited-area NWP model, the Weather Research and Forecasting model (WRF) and its global model-the European Centre for Medium Range Weather Forecasting (ECMWF) deterministic model. We simulate day-ahead Global Horizontal Irradiance (GHI) for a location in Qinghai, China. WRF model configurations optimized by the type of day (cloud amount) are then implemented based on the ECMWF model forecast of cloud amount. A regression model to combined ECMWF and WRF model forecasts is then trained. The optimised coefficients (weights) of ECMWF and WRF show increasing WRF importance with higher cloud amounts and the combination out-performs the ECMWF input by 5.2% and the best WRF configuration by 7.2% on a 2.5-month testing set. The performance of the combined model increased with observed cloud amount where the combined model out-performed the ECMWF model by 12.6% for cloudy days indicating the relative importance of physical downscaling for the simulation of clouds.