1 Introduction

In the UK both precipitation and temperature extremes are projected to increase significantly with anthropogenic climate change (Hanlon et al. 2021; Christidis et al. 2020; Cotterill et al. 2021; Kendon et al. 2014). However, less attention has been given to how much atmospheric circulation changes are contributing to these long-term projected changes in addition to the thermodynamic changes, particularly within the transition seasons. Atmospheric circulation impacts the probability of extreme events over Europe significantly, with many attribution studies of extreme events examining changes in probability of the event for given atmospheric conditions (Christidis et al. 2015; Schaller et al. 2016).

Atmospheric circulation and weather types have been examined for a long time for the UK including the seven Lamb weather types (Lamb 1972) which have been developed and used further (Jenkinson 1977; Jones et al. 1993). Since then, there has a been a shift to creating circulation types using k-means clustering. Fereday et al. (2008) created 10 weather patterns for each season using mean sea level pressure (psl), Ferranti et al. (2015) created 4 weather patterns using geopotential height and most recently Neal et al. (2016) created the 30 weather patterns which are used at the Met Office as a forecasting tool. The 30 weather patterns cover the UK and the north-west European domain. The benefit of using the 30 weather patterns is that they don’t only provide information on broadscale circulation, but variations within it, including seasonality (Neal et al. 2016). Studies have started using the 30 weather patterns for research into future and current circulation changes and their relation to climate hazards using the UK Climate Projections (UKCP). Pope et al. (2021) finds that circulation changes will lead to more dry summers and more westerly and cyclonic conditions in the winter and Kendon et al. (2020) examine which patterns are most likely to lead to record breaking warm winter events such as February of 2019 in the UK.

There is still a lot of uncertainty on atmospheric changes over the European region, especially due to their link to the Atlantic Meridional Overturning Circulation (AMOC). There are large differences in the projected weakening of the AMOC in models, the strength thereof would impact the changes in European climate through the uncertainty of heat transport (Haarsma et al. 2015). Rahmstorf et al. (2015) provide evidence of the slowdown of AMOC in the twentieth century leading to a cooling of a region in the North Atlantic, with this subpolar cooling much weaker in the CMIP5 models than seen in the observations.

Some studies are already finding differences in atmospheric circulation over Northern Europe. Schaller et al. (2016) using current and counterfactual climate conditions find that an increase in low-pressure systems northwest of the UK in winter and change in zonal flows is already occurring because of anthropogenic influence. Another study finds that changes in atmospheric conditions using clustering are already being seen in the observations with the summer extending into early autumn and late spring since 1948 (Vrac et al. 2013). Seasonal changes are likely to have a big impact on the Agriculture, Energy and Tourism sectors (Vrac et al. 2013). An analysis using the UKCP projections finds that the summer dry season is longer in the future with less soil moisture in the summer and autumn (Pirret et al. 2020). The UK has already experienced impacts from droughts, a dry period from 2016 to 2019 caused a number of impacts including agriculture, wildfires and freshwater ecosystems, especially in the South-East of England (Turner et al. 2021).

This paper will look at whether circulation changes through weather patterns may be contributing to these changes and by how much. This work will address three main scientific questions. Firstly, using attribution runs, has anthropogenic influence noticeably changed the frequency of weather patterns over the UK already? Secondly, are there future changes in seasonality due to circulation, and if so, how does this change vary between emissions scenarios and models? Finally, what is the impact of these changes in seasonality on mean and extreme rainfall for the UK? This last question will be addressed by isolating the atmospheric circulation signal from the overall signal.

2 Data and methods

2.1 Data

2.1.1 ERA5 weather patterns

To evaluate how well each of the models capture the frequency of daily weather patterns, the ERA5 reanalysis data from the European Centre for Medium Range Weather Forecasting (ECMWF) covering the period 1979–2010 was used. This data with a horizontal resolution of 31 km and 137 vertical levels up to 0.01 hPa uses Observations, Models and a 4D-Var data assimilation scheme to produce high resolution global reanalyses (Hersbach et al. 2020). The weather patterns were assigned to the ERA5 daily mean sea level pressure data using the method described in Neal et al (2016).

2.1.2 HadGEM3-A

The HadGEM3-A is a global atmospheric climate model created in the UK Met Office Hadley Centre with a horizontal resolution of 60 km at the mid-latitudes. This contains 15 ensemble members for each ALL forcings (historical) ensemble and NAT forcings (historical natural-only) ensemble from 1960 to 2013 and larger sets of ensembles beyond that until present. The historical natural-only runs contain anthropogenic forcings held at 1850 levels with changes in Natural forcings only, and the historical runs represent all forcings as they actually evolved (Vautard et al. 2018). The boundary conditions for both ALL and NAT simulations use observed sea surface temperatures and sea ice, which are adjusted in the NAT ensemble by having an estimate of the changes due to human induced climate change removed. More information can be found in Christidis et al. (2013) and Ciavarella et al. (2018).

2.1.3 UKCP PPE

The UK Climate Projections contains a 28 global climate model ensemble (a 15 member perturbed parameter ensemble (PPE) based on the Met Office Hadley Centre model and 13 models sub-selected from CMIP5) for the period 1899–2099 for two emission scenarios: Representative Concentration Pathways (RCP) RCP 2.6 and RCP 8.5 (Lowe et al. 2019). For this work we use the 15 Hadley Centre PPE members. The 15 PPE members are based on the HadGEM model at GC 3.05 configuration, and the construction of the PPE is discussed in Sexton et al. (2021) and Yamazaki et al. (2021). A summary of the simulations is available in Pope et al. (2021). This paper uses daily precipitation and pressure at sea level (psl) data for both emission scenarios. These RCP scenarios were created to account for a range of outcomes of future greenhouse gas emissions, including assumptions on mitigation, economic factors and population changes (Lowe et al. 2019). The UKCP PPE weather pattern data for RCP 2.6 and RCP 8.5 between 1899 and 2099 was the same data created and used in Pope et al. (2021) from daily psl data. These two RCP scenarios represent the low-end and high-end emissions scenarios respectively of the four outcomes created.

2.1.4 HadUK- grid observations

The HadUK-Grid gridded observations dataset is used to evaluate climate model data for regional precipitation in the UK. This dataset produced by the Met Office Hadley Centre is based on land surface observations and is available both at a number of gridded resolutions and a regional level for 16 administrative regions (Hollis et al. 2019). This paper uses the regional daily precipitation dataset over the time period 1970–2000.

2.1.5 CMIP6 models

The following four CMIP6 models were chosen because they were the only ones that had three or more ensemble members with daily psl data for both a high future emissions scenario and historical period with a horizontal resolution no lower than 100 km.

  1. (i)

    MRI-ESM2-0 (CMIP6): Data from CMIP6 from the Earth System Model MRI-ESM2-0, created by the Meteorological Research Institute in Japan, with a horizontal resolution of 100 km (Yukimoto et al. 2019). The data used for this work came from two of the experiments, historical natural-only (containing only changes in Natural forcings) and a shared socio-economic pathway (SSP) SSP5-8.5. There are five SSP scenarios which represent anthropogenic emissions from five different projected socio-economic pathways, with SSP5-8.5 being the high-end scenario (Meinshausen et al. 2020). The high-end scenario SSP5-8.5 differs slightly from the high-end RCP scenario (RCP 8.5), with higher future levels of Carbon Dioxide, but lower levels of Methane in SSP5-8.5 by the end of the twenty-first century (Meinshausen et al. 2020).

  2. (ii)

    HadGEM3-GC31-MM (CMIP6): Data from CMIP6 from the Global Climate Model HadGEM3-GC31 created by the UK Met Office Hadley Centre with a horizontal resolution of 60 km at the mid-latitudes (Ridley et al. 2019). The data used from this work came from two experiments, historical and a shared socio-economic pathway SSP5-8.5.

  3. (iii)

    NCAR-CESM2-WACCM (CMIP6): Data from CMIP6 from the Global Climate Model CESM2-WACCM created by the National Center for Atmospheric Research in the USA with 100 km horizontal and atmospheric resolution (Danabasoglu 2019). The data used from this work came from two experiments, historical and a shared socio-economic pathway SSP5-8.5.

  4. (iv)

    EC-Earth3 (CMIP6): Data from CMIP6 from the Global Climate Model EC-Earth 3.3 created by a consortium of European meteorological and research institutions, with 100 km resolution for both atmosphere and ocean/land (EC-Earth 2019). The data used from this work came from two experiments, historical and a shared socio-economic pathway SSP5-8.5.

2.1.6 The Met Office 30 UK weather patterns

The Met Office 30 Weather Patterns (Neal et al. 2016) are a set of patterns used in operational medium- to long-term weather forecasting aiming to cover the large range of circulation types over the UK and the European Domain (the North Atlantic-European (NAE) Domain) (Fig. 1). These patterns verified by UK Meteorologists were created using the EMULATE (European and North Atlantic daily to multidecadal climate variability- Ansell et al. 2006) sea level pressure reanalysis data from 1850 to 2003 in combination with a k-means clustering algorithm (Philipp et al. 2007), over the NAE domain (30°W–20°E; 35°–70°N). Each of the weather patterns has an associated climatology for meteorological variables (such as surface temperature or wind direction). These weather patterns have been used in understanding potential impacts (such as flooding or volcanic ash – Neal et al. 2018; Harrison et al. 2022) and in understanding future changes in British Isles climate (Pope et al. 2021). The lower numbers have lower psl anomalies, which occur more in the summer and will be referred to as summer-types and the high numbered weather patterns which have very high psl anomalies occurring more in the winter, will be referred to as winter-types (Neal et al. 2016).

Fig. 1
figure 1

Definition maps for the set of 30 predefined weather patterns from Neal et al. (2016) which have been used in this study (and also Pope et al. 2021). Mean sea level pressure (MSLP) anomalies plotted as filled contours (hPa) and MSLP mean values plotted in foreground (2 hPa intervals)

2.2 Methods

2.2.1 Allocating weather patterns to climate model data

To assign weather pattern values to climate model data, we first calculated the daily psl anomalies for each model day and a baseline anomaly period for the corresponding model. Fields of daily anomalous psl are regridded to the 5° × 5° grid of the EMULATE dataset (which is used by the weather patterns) and then each day is independently assigned to its closest match from of the 30 weather patterns. This is done through the area-weighted sum of squares difference (referred to as the distance) calculated at each grid point, with the weather pattern with the smallest distance value being assigned to that day in the model.. The methodology used here is described in more detail in Neal et al. (2016) and is consistent with the methodology applied in Pope et al. (2021).

The climate model data simulated psl anomalies are created by taking a climatological average over a baseline period for each day of the year and comparing that to the respective day of the year for the daily climate model data. The climatologic average is the mean psl value for each grid box on a given day of the year over the baseline period. The climatological average was created for each ensemble member individually, with the choice of baseline period the same for each ensemble member. This was applied to the four CMIP6 models used and the HadGEM3-A model runs. All five models used the historical natural-only/historical runs from 1960/70 to 2000/10 as the baseline period. The details of the choice of baseline period for each model can be found in Supplementary Table 1. The UKCP weather pattern data used was created in the same way but had a longer baseline period of 1900–1999 using historical runs (Pope et al. 2021).

2.2.2 Measuring confidence in the results

To examine the strength of changes between two different time-slices or emission scenarios in the frequency of weather patterns for a given weather pattern, wp we use a bootstrapping method. Weather patterns are compared individually instead of all 30 at once, because consecutive patterns are not directly related to each other and are unique.

The following steps are used to examine whether weather pattern changes between two time-slices named ts1 and ts2 for a specific weather pattern wp in a specific season is significant:

  • Firstly, for each year of data in the time-slice calculate the number of days (day-count) with pattern wp over the season of interest. For a time-slice of length n years and m ensemble members there are n x m season-year day-count values for each weather pattern.

  • For each time-slice these season-year day-count values are bootstrapped with replacement 10,000 times.

  • After every bootstrap, the day-count values are summed across all years in the time-slice giving totals counts for time-slices ts1 and ts2, count_ts1 and count_ts2 respectively.

  • The change for a given pattern is significant if count_ts1 > count_ts2 for 97.5% or count_ts1 < count_ts2 for 2.5% of the bootstrapped values for that pattern. This is a measure of how much the ranges overlap which we will refer to as a significant range difference.

2.2.3 Isolating the circulation signal using weather pattern frequencies

To isolate the climate signal for circulation, we use both the changes in daily frequencies of the 30 weather patterns and climatology profiles of each pattern for the chosen variable. The method follows as such for examining the impact of a circulation signal change in variable X, which for this example is extreme precipitation days (99.5th percentile) for the season of interest.

Climatology Profile: Firstly, a baseline period for the climatology was chosen, which is the UKCP Global data between 1900 and 1950. This is because there are a large number of ensemble members (15) over a long time period and hence a sufficient number of extreme events to produce profiles for each pattern. The time period also occurs far in the past where signals of anthropogenic climate change have not emerged significantly yet. Using daily rainfall data associated with each pattern using the UKCP data, the probability of an extreme rainfall day (p_wp) on any day given its weather pattern (wp), is calculated. The definition of an extreme rainfall day for a particular season, is a day where the daily precip exceeds the 99.5th percentile over the baseline period (1900–1950). This was also carried out for mean rainfall, where in the climatology profile instead of it being a probability of an extreme rainfall day- it was the mean rainfall for a day of a given pattern.

Frequency profile: For each of the models examined, the pattern frequency is calculated for the two time-slices to be compared individually. In each case, a historical/ historical natural-only time-slice (ts_past) is compared to a future time time-slice (ts_future) under a given emissions scenario over the following time periods:

  1. (i)

    Historical Period/historical natural-only period- January 1970 to December 1999.

  2. (ii)

    Future Period (RCP 2.6/8.5 and SSP5-85)- January 2070/71 to November 2099.

Based on the weather pattern frequencies in the two time-slices being compared, the expected number of days for each weather pattern wp (number_dayswp) over a 20-year period for both ts_past and ts_future for the chosen season is calculated. Then using p_wp from the climatology profile and number_dayswp, the projected percentage increase in the number of extreme days between time slices ts_past and ts_future based off the change in frequency of weather patterns alone can be calculated using formulas 1 and 2.

$${\text{Expected extreme rainfall days}}_{ts} = \mathop \sum \limits_{wp = 1}^{30} {\text{number}}\_{\text{days}}_{wp} *{\text{prob extreme day}}_{wp}$$
(1)
$$\Delta \,{\text{change in number of extreme days}}\;( \% ) = \left( \frac{{{\text{Expected extreme rainfall days}}_{{{\text{ts}}}\_{{{\text{future}}}}} }} {{{\text{Expected extreme rainfall days}}_ {{{\text{ts}}}\_{{{\text{past}}}}} }} - 1 \right)*100\% $$
(2)

To produce the same calculation, but for mean precipitation instead of extreme precipitation, formulas 3 and 4 are used where mean_precipitationwp is the mean daily precipitation for weather pattern wp. The uncertainty in the results from formulas (14) for the circulation signal is estimated using the model range. These results give an estimate of the change in the variable of interest due to the isolated Atmospheric Circulation Signal (ACS).

$${\text{Expected 20 year rainfall}}_{ts} = \mathop \sum \limits_{wp = 1}^{30} {\text{number}}\_{\text{days}}_{wp} *{\text{mean}}\_{\text{precipiation}}_{wp}$$
(3)
$$\Delta \user2{ }{\text{change in mean rainfall}}\user2{ }\left( \user2{\%} \right) = \left( {\frac{{{\text{Expected 20 year rainfall}}_{{{\text{ts}}}\_{{{\text{future}}}}} }} {{{\text{Expected 20 year rainfall}}_ {{{\text{ts}}}\_{{{\text{past}}}}} }} - 1} \right)*100\%$$
(4)

To show these changes from the ACS in context, the overall signal change is also calculated. The overall signal change in these variables between a past and future time-slice for a given emissions scenario is calculated separately using data from the UKCP Global in each time-slice. The uncertainty on the overall signal is estimated using bootstrapping with replacement for the daily rainfall data over the region. After bootstrapping the daily precipitation data 10 000 times, the 99.5th percentile of daily rainfall and mean yearly rainfall for the chosen season are calculated 10 000 times. The 95% confidence intervals calculated for each are then compared between time-slices to produce upper and lower estimates of the overall change in each variable between the time-slices.

3 Model evaluation

3.1 Attribution runs

The model used for the attribution of weather patterns, HadGEM3-A is validated for weather pattern frequencies over the time period 1979–2010, using the ERA5 reanalysis data. The 15 ensemble members from the ALL forcings runs of the HadGEM3-A are shown as boxplots in Fig. 2, where strong agreement can be seen for all seasons between the HadGEM3-A model and ERA5. The ERA5 reanalyses fall within the model range for 29/26/26/28 out of 30 patterns for Winter/Spring/Summer/Autumn (djf/mam/jja/son) respectively. For the few patterns where the ERA5 value lies outside the model range, the ERA5 value is no more than 23% outside of the ensemble range for the weather pattern-season combination.

Fig. 2
figure 2

Seasonal daily frequencies of the 30 weather patterns over the UK in both the ERA5 reanalyses and the HadGEM3-A model ALL forcings runs between 1979 and 2010. The points making up the boxplot for each of the patterns represent the 15 individual ensemble members, where the shaded box covers the interquartile range, the whiskers are the min and max values, and the dots represent outliers

3.2 Model validation of UKCP and CMIP6 model data for weather patterns

Pope et al. (2021) evaluates both the annual frequency of the 30 weather patterns and persistence within the UKCP Global climate projections, finding that there is good agreement between the UKCP Global members and ERA5. A further evaluation of seasonal frequencies for UKCP Global 15 PPE members used in this paper, shows the strong agreement between the model and ERA5, with the ERA5 data falling within the UCKP 15 PPE models range for 21/29/25/26 out of the 30 patterns for djf/mam/jja/son respectively. The four CMIP6 models chosen HadGEM3-GC31-MM (4 ensemble members), MRI-ESM2-0 (5 ensemble members), NCAR-CESM2-WACCM (3 ensemble members) and EC-Earth3 (7 ensemble members) when combined and compared against the ERA5 data also show relatively strong agreement. The ERA5 data falls within the CMIP6 models range for 24/25/29/29 out of the 30 patterns for djf/mam/jja/son respectively.

3.3 Model validation of UKCP Global for rainfall

The UKCP Global 15 PPE between 1900 and 1950 is used as the baseline for the rainfall climatology for each of the 30 weather patterns. The model validation is chosen over a later period for the model, where the quality of observations is improved (1970–2000).

The UKCP Global projections show regional variation in its ability to capture UK rainfall in Autumn (Fig. 3). The model agrees more with the observations for UK regions without mountains for mean and extreme rainfall (99th percentile). Therefore, for the climatological examination of the changes in UK rainfall in Autumn, this work focusses on three distinctly located regions where the observations for both mean yearly rainfall and the 99th percentile of daily rainfall in Autumn fall within the model ensemble range. These are chosen as South East England, the West Midlands and Yorkshire and Humber with the model spatial coverage of that region shown in Fig. 4.

Fig. 3
figure 3

Validation of the UKCP Global (60 km) 15 PPE for mean yearly and extreme daily rainfall (99th percentile) over the administrative regions of the UK when compared to the HadUK-gridded observations (red dots) between 1970 and 2000 in Autumn

4 Results

4.1 Attribution of weather patterns

To assess whether there has been an anthropogenic influence on the patterns seen over the UK, seasonal frequencies of weather patterns are compared for HadGEM3-A, which contains both 15 members ensembles with Natural Forcings only and with All Forcings (Fig. 5). The figure shows that for all seasons between 1980 and 2010 there are very few noticeable differences in pattern frequencies between ALL and NAT with there being no significant range differences for 28/26/23/25 of the 30 weather patterns for djf/mam/jja/son respectively. The same analysis was carried out for MRI-ESM2-0 for the 5 NAT and 5 ALL ensemble members between 1980–2010. This showed a few extra noticeable differences with no significant range differences for 22/27/28/19 out of the 30 patterns for djf/mam/jja/son respectively. The patterns out of the 30 that showed significant range differences were compared for both models and there was very little agreement between the two models on these weather pattern changes with no overarching trends (Supplementary Table 2). The only weather patterns that showed the same change in both models were pattern 20 (increasing in frequency in SON), pattern 5 (decreasing in frequency in SON), and pattern 29 (increasing in frequency in MAM). Overall, based on these models, no clear signal has appeared suggesting that anthropogenic climate change has not influenced the frequency of weather patterns significantly before 2010.

Fig. 4
figure 4

The grid boxes used for the 3 Regions of the UK analysed for rainfall changes using the N216 grid for UKCP Global

4.2 Future climate signal in weather patterns

To assess the influence of anthropogenic climate change on weather patterns in the future, the UKCP Global is examined for a past period (1960–2010), as well as a future period (2071–2100) for both a low (RCP 2.6) and high (RCP 8.5) emissions scenario. Given very few significant range differences were found between ALL and NAT forcings from 1980 to 2010, the historical period is used as a baseline for weather pattern frequencies when looking at the influence of anthropogenic climate change on the ACS.

Fig. 5
figure 5

Seasonal frequencies of the 30 weather patterns over the UK of 15 ensemble members of HadGEM3-A ALL forcings runs (red) and NAT forcings runs (blue) from 1980 to 2010. The points making up the boxplots for each of the 30 patterns represent the frequency values for the 15 individual ensemble members, where the shaded box covers the interquartile range and the whiskers are the min and max values

The results show that there are significant range differences for the majority of the 30 weather patterns within all seasons between UKCP Global historical and future RCP 8.5. Figure 6 shows the differences for Autumn with the boxplots for future RCP 2.6 falling between the past and future RCP 8.5 time slices, indicating that anthropogenic influence is behind these changes in weather pattern frequencies. The main pattern seen in this figure for Autumn, is that there are significant increases in the lowest numbered patterns (summer-types) and a decrease in the highest numbered patterns (winter-types). This is also seen for the season of Summer, but not Winter or Spring in the UKCP projections (Supp. Figures 4, 5, 6). These differences seen between the lowest and highest emission scenarios in the future are much greater than the differences seen between the past and the future low emissions scenario.

Fig. 6
figure 6

Plot showing Autumn seasonal frequencies of daily pattern numbers for the 30 patterns for different time-slices and emissions scenarios using UKCP Global. The points making up the boxplots for each of the 30 patterns represent the frequency values for the 15 individual ensemble members, where the shaded box covers the interquartile range and the whiskers are the min and max values

The CMIP6 models also show the same trends in Autumn but to a lesser extent, with a large increase in the lower numbered summer-type weather patterns and a decrease in the higher-numbered winter-type ones (Fig. 7). This change can also be seen in the Summer in the CMIP6 models (Supplementary Table 3), however such a trend is not seen in Spring or Winter. This extension of the summer seen in the models could have big impacts on a number of weather variables including rainfall- which is examined in detail in the next section. It is also important to note that these changes have not appeared to occur in weather patterns seen when comparing ALL and NAT forcings before 2010 in the attribution section.

Fig. 7
figure 7

Plot showing Autumn seasonal frequencies of daily pattern numbers for the 30 weather patterns as boxplots for a past period (blue) and a future period with SSP5-8.5 (red) for the CMIP6 models used. The 19 points using CMIP6 models come from MRI-ESM2-0 (5 ensemble members), HadGEM3-GC31 MM (4 ensemble members), NCAR-CESM2-WACCM (3 ensemble members) and EC-Earth3 (7 ensemble members). The points making up the boxplots for each of the 30 patterns represent the frequency values for the 19 individual ensemble members, where the shaded box covers the interquartile range and the whiskers are the min and max values

To examine when this signal emerges, the combined frequency of the predominantly summer type patterns (1–6) is plotted against the combined frequency of the predominantly winter type patterns (25–30) for the lowest and highest emissions scenarios for UKCP Global RCP 2.6 and RCP 8.5 respectively (Fig. 8). This shows that the signal could emerge as early as 2025 under these scenarios with the two emissions scenarios deviating to outside the 1900–2000 range from around 2035 onwards. By Autumn 2095 the frequency of winter-type patterns in SON could decrease by over a third under the highest emissions scenarios and summer-types increasing by a quarter. Under low emissions this change in seasonality in Autumn is much reduced as can be seen in Fig. 8.

Fig. 8
figure 8

The transient evolution of daily weather pattern frequencies in SON for summer-types (Patterns 1–6) and winter-types (Patterns 25–30) between 1900 and 2100 in 10-year time-slices for the 15 Global UKCP PPE under emissions scenarios RCP 2.6 (dashed lines) and RCP 8.5. The shaded regions responding to summer and winter types represent the 10-year time-slice frequency ranges between 1900–2000

4.3 Impacts of weather pattern changes on UK rainfall

To investigate whether the increase in summer-type patterns impacts rainfall, the ACS is isolated and compared to the overall change. The main focus of this analysis is for Autumn, given the strong signal seen in Sect. 4.2. To assess the relative impact of the ACS against the overall signal, the method used in Sect. 2.3 is applied for a past time period (1971–2000) and the future time period (2071–2100) under different emissions scenarios for Autumn, using the three regions chosen (Fig. 4). The climatology profile for each weather pattern for both the probability of an extreme rainfall day in Autumn (SON) (a day exceeding 99.5th percentile) and mean SON daily rainfall is calculated from the UKCP Global dataset over a baseline period running from 1900 to 1950.

Future circulation changes in the models when isolated show a 9–12% decrease in mean rainfall in SON from 1985 to 2085 over the three English regions (Table 1). The signal of this change is true for all five models examined under high emissions scenarios, with the strongest signals seen in the two Hadley Centre models used; UCKP Global RCP 8.5 and HadGEM3-GC31. The overall signal for Autumn mean rainfall, a decrease of 9–11% using the UKCP Global RCP 8.5 simulations is almost identical to the ACS, a decrease of 9–12%. This suggests that the overall mean rainfall signal leading to drier Autumns in the future is primarily driven by atmospheric circulation changes. Future ACS changes in the models when isolated also shows a 20–21% decrease in extreme daily rainfall events in Autumn (Table 1). The overall signal over the same time interval, however, shows a 71–103% increase in extreme rainfall days in SON using the UKCP Global RCP 8.5 simulations. This would likely be even higher if it wasn't moderated by the 20–21% decrease due to the ACS.

Table 1 The change in rainfall indices between 1985 and 2085 for both the ACS and overall signal in Autumn

There are also big differences for rainfall indices between emission scenarios RCP 2.6 and RCP 8.5 for UKCP Global, both for the ACS and overall signal (Table 2). Overall, under a low emissions scenario the increase in Autumn extreme rainfall days is around 21–47% depending on the region, compared to 71–103% under the highest emissions scenario at the end of this century. The difference is even more pronounced for mean rainfall where on average there is only a 5% decrease in mean rainfall compared to a 10% decrease, when comparing RCP 2.6 and RCP 8.5 respectively. The ACS is also significantly reduced for both indices under a lower emissions scenario, with the signal between 2–3 times greater for the higher emissions scenario.

Table 2 The change in rainfall indices between 1985 and 2085 for both the ACS and the overall signal for different UKCP Global emissions scenarios

5 Discussion and conclusions

In this paper we have found a climate circulation signal showing a predicted increase in the frequency of predominantly summer-type and a decrease in predominantly winter type weather patterns in both Autumn and Summer over the UK. This climate signal could emerge as early as the 2020s based on the UKCP projections. The estimated impact of what this signal means for Autumn mean and extreme rainfall is examined for multiple English regions. We find that these changes in the frequency of weather patterns could result in up to a 4–14% decrease in mean Autumn rainfall for English regions by 2085, dependent on the emissions scenario (RCP 2.6 – RCP 8.5). The overall mean reduction of rainfall in Autumn seen in the UKCP Global is of almost the same magnitude as the reduction in rainfall due to changes in the frequency of weather patterns, suggesting that circulation is playing a dominant role in these changes. This combined with drier summers predicted by the UKCP projections (Lowe et al. 2019) has big implications for future risks of drought over the UK, especially given parts of England have only recently experienced some of the impacts from drought.

There are however both strong differences between models and emission scenarios for the magnitude of this change. The decrease in Autumn rainfall due to circulation in the UKCP Global is three times smaller for the lowest emissions scenario RCP 2.6 than the highest emissions scenario RCP 8.5. Furthermore, this circulation climate signal seen in Autumn producing drier weather is around 1.5–2 times smaller in MRI-ESM2-0 and EC-Earth3 from CMIP6 compared to HadGEM3-GC31 from CMIP6 under the same future emissions scenario, with the difference between the models being the main uncertainty in these results.

Despite drier Autumns, there is a large increase in the number of extreme rainfall days in SON of up to 47–103% (RCP 2.6- RCP 8.5) by the end of this century. We do find that changes in circulation for Autumn moderate this increase and the circulation signal when isolated from the overall signal acts to reduce the number of extreme rainfall days over parts of England by 20–21%. This could explain why we see a larger increase in winter daily rainfall extremes than in Autumn by the end of this century in the UKCP projections (Cotterill et al. 2021).

Despite these future changes, no obvious detectable changes have been found when comparing the NAT forcings runs to the ALL forcings runs for any of the seasons before 2010. Even for those of the 30 patterns that show small differences between ALL and NAT for models, there is very little agreement between the two models examined HadGEM3-A and MRI-ESM2-0. One could argue that a greater number of simulations than the 20 each from ALL/NAT are required to examine changes for weak signals, but these results suggest that this Extension of Summer signal is not visible before 2010. The UKCP Global suggests that under both the low and high emissions scenarios, this Autumn circulation signal with more summer-type and less winter-type patterns will emerge in the next decade or so (Fig. 8). The magnitude of this change by the end of this century, however, is heavily dependent on the emissions scenario.

Despite the detail of the 30 weather patterns the climate models used had strong agreement with the Observations for all seasons, especially the HadGEM3-A. This may be expected as the model uses SST and Arctic Sea Ice as boundary conditions, therefore constraining it more to the global patterns over the time period examined. However, there are limitations to this method at only looking at the daily frequencies of weather patterns. Firstly, the pattern assigned from the daily psl data for a given day could be a combination of two weather patterns if there is a transition between patterns that day. Secondly, the weather patterns based on observations from 1850 to 2010, may evolve and be less suited to future patterns and associated impacts. Multidecadal variability may have also impacted the results when comparing the two time slices. But given the atmospheric climate signal found in this work in Autumn is seen across all models, there can be confidence in the results.

This paper mainly focusses on the impacts of the circulation signal on rainfall for Autumn, but there are a number of other possible impacts as a result of this climate signal including heat stress and impacts for agriculture. This could be a good area for future research with the number of sectors and industries that these long-term changes could impact, with large differences seen between the lowest and highest emissions scenarios. The impacts of these could also be examined using the UKCP18 2.2 km convective permitting model, to look at the regional impacts of these long-term changes.