AbstractLarge uncertainty exists in the sign of long‐term changes in regional scale mean precipitation across the current generation of global climate models. To explore the physical drivers of this uncertainty for New Zealand, here we adopt a storyline approach applying cluster analysis to spatial patterns of future projected seasonal mean precipitation change across CMIP6 models (n = 43). For the winter precipitation change signal, the models split roughly into two main groups: both groups have a very robust wet signal across the west coast of the South Island but differ notably in terms of the sign of precipitation change across the north of the North Island. These far north winter precipitation differences appear related to how far the Hadley cell edge and regional eddy‐driven jet shift across the models relative to their historical positions. In contrast, for summer, most models have a markedly weaker and spatially non‐uniform response, where internal variability often plays a large role. However, a small group of models predict a robust wet signal across most of the country in summer. This “wet model” group is characterized by a regional La Niña‐like increase in high pressure shifted further to the south‐east of New Zealand, associated with more frequent north‐easterly flow over the country and accompanied by significant warming of local sea surface temperatures. This regional circulation response appears related to changes in stationary Rossby wave paths as opposed to changes in La Niña occurrence frequency itself.