A key-predictand and key-station approach was employed in downscaling general circulation model outputs to monthly evaporation, minimum temperature (Tmin) and maximum temperature (Tmax) at five observation stations concurrently. Tmax was highly correlated (magnitudes above 0.80 at p ≤ 0.05) with evaporation and Tmin at each individual station, hence Tmax was identified as the key predictand. One station was selected as the key station, as Tmax at that station showed high correlations with evaporation, Tmin and Tmax at all stations. Linear regression relationships were developed between the key predictand at the key station and evaporation, Tmin and Tmax at all stations using observations. A downscaling model was developed at the key station for Tmax. Then, outputs of this downscaling model at the key station were introduced to the linear regression relationships to produce projections of monthly evaporation, Tmin and Tmax at all stations. This key-predictand and key-station approach was proved to be effective as the statistics of the predictands simulated by this approach were in close agreement with those of observations. This simple multi-station multivariate downscaling approach enabled the preservation of the cross-correlation structures of each individual predictand among the stations and also the cross-correlation structures between different predictands at individual stations.