A method, based on climate pattern-scaling, has been developed to expand a small number of projections of fields of a selected climate variable ( X ) into an ensemble that encapsulates a wide range of model structural uncertainties. The method described in this paper is referred to as the Ensemble Projections Incorporating Climate model uncertainty (EPIC) method. Each ensemble member is constructed by adding contributions from (1) a climatology derived from observations that represents the time invariant part of the signal, (2) a contribution from forced changes in X where those changes can be statistically related to changes in global mean surface temperature ( T global ), and (3) a contribution from unforced variability that is generated by a stochastic weather generator. The patterns of unforced variability are also allowed to respond to changes in T global . The statistical relationships between changes in X (and its patterns of variability) with T global are obtained in a "training" phase. Then, in an "implementation" phase, 190 simulations of T global are generated using a simple climate model tuned to emulate 19 different Global Climate Models (GCMs) and 10 different carbon cycle models. Using the generated T global time series and the correlation between the forced changes in X and T global , obtained in the "training" phase, the forced change in the X field can be generated many times using Monte Carlo analysis. A stochastic weather generator model is used to generate realistic representations of weather which include spatial coherence. Because GCMs and Regional Climate Models (RCMs) are less likely to correctly represent unforced variability compared to observations, the stochastic weather generator takes as input measures of variability derived from observations, but also responds to forced changes in climate in a way that is consistent with the RCM projections. This approach to generating a large ensemble of projections is many orders of magnitude more computationally efficient than running multiple GCM or RCM simulations. Such a large ensemble of projections permits a description of a Probability Density Function (PDF) of future climate states rather than a small number of individual story lines within that PDF which may not be representative of the PDF as a whole; the EPIC method corrects for such potential sampling biases. The method is useful for providing projections of changes in climate to users wishing to investigate the impacts and implications of climate change in a probabilistic way. A web-based tool, using the EPIC method to provide probabilistic projections of changes in daily maximum and minimum temperatures for New Zealand, has been developed and is described in this paper.