We introduce a novel technique to bias correct climate model output for impact simulations that preserves its physical consistency and multivariate structure. The methodology considerably improves the representation of extremes in climatic variables relative to conventional bias correction strategies. Illustrative simulations of biosphere–atmosphere carbon and water fluxes with a biosphere model (LPJmL) show that the novel technique can be usefully applied to drive climate impact models.