Earth system models (ESMs) produce the highest-quality future climate data available, but they are costly to run, so only a few runs from each model are publicly available. What is needed are emulators that tell us what would have happened, if we had been able to perform as many ESM runs as we might have liked. Much of the existing work on emulators has focused on deterministic projections of average values. Here we present a way to imbue emulators with the variability seen in ESM runs.