Dynamical mechanisms linking Indian monsoon precipitation and the circumglobal teleconnection

The circumglobal teleconnection (CGT) is an important mode of circulation variability, with an influence across many parts of the northern hemisphere. Here, we examine the excitation mechanisms of the CGT in the ECMWF seasonal forecast model, and the relationship between the Indian summer monsoon (ISM), the CGT and the extratropical northern hemisphere circulation. Results from relaxation experiments, in which the model is corrected to reanalysis in specific regions, suggest that errors over northwest Europe are more important in inhibiting the model skill at representing the CGT, in addition to northern hemisphere skill more widely, than west-central Asia and the ISM region, although the link between ISM precipitation and the extratropical circulation is weak in all experiments. Thermal forcing experiments in the ECMWF model, in which a heating is applied over India, suggest that the ISM does force an extratropical Rossby wave train, with upper tropospheric anticyclonic anomalies over east Asia, the North Pacific and North America associated with increased ISM heating. However, this eastward-propagating branch of the wave train does not project into Europe, and the response there occurs largely through westward-propagating Rossby waves. Results from barotropic model experiments show a response that is highly consistent with the seasonal forecast model, with similar eastward- and westward-propagating Rossby waves. This westward-propagating response is shown to be important in the downstream reinforcement of the wave train between Asia and North America.


Introduction
The circumglobal teleconnection (CGT) is a prominent wavenumber-5 teleconnection mechanism that manifests in the northern hemisphere summer. First identified by Ding and Wang (2005), it has centres of action over west-central Asia, east Asia, the North Pacific, North America and northwest Europe, where the upper-level pressure variations are nearly all in phase (Fig. 1). It is associated with significant temperature and precipitation anomalies in these regions (see Figures 7 and 9 from Ding and Wang 2005), and so accurate representation of this mechanism in seasonal forecasting systems could be an important source of subseasonal to seasonal forecast skill.
The convective heat source over the Bay of Bengal and the Indian summer monsoon (ISM) region and its associated upper-tropospheric divergence can serve as a Rossby wave source. Wang (2005, 2007) hypothesised that the interaction of ISM precipitation with the extratropical circulation plays an important role in the maintenance of the CGT. A number of studies have found that strong ISM precipitation is associated with anticyclonic anomalies over west-central Asia and east Asia, near Japan (Wang et al. 2001;Wu and Wang 2002). The propagation of Rossby waves from the jet exit region over northwest Europe to Asia results in the development of an anomalous upper-level high over west-central Asia. The easterly anomalies to the south of this high strengthen the monsoon easterly vertical shear over northern India, which enhances the convection in this region. This convection and associated diabatic heating in turn excites a Rossby wave response, similar to those of Gill (1980) and Rodwell and Hoskins (1996), which reinforces the west-central Asian high. This then enhances downstream circulation anomalies through the propagation of these waves along the waveguide (the westerly jet) (Liu and Wang 2013). Linked to this mechanism, Lau and Weng (2002) found that North America summer precipitation variability is connected to a wave train that originates in Asia. Lau et al. (2004) named this the "Tokyo-Chicago Express".
Related to these mechanisms is the Silk Road Pattern (SRP), which is associated with stationary Rossby waves trapped in the Eurasian jet stream (Lu et al. 2002;Enomoto et al. 2003). Stephan et al. (2019) examined the decadal variability of the SRP and identified a positive feedback loop between the SRP and vertical motion over India and the Mediterranean, whereby more coherent monsoon precipitation leads to increased descent over the Mediterranean region, via the monsoon-desert mechanism (Rodwell and Hoskins 1996), which subsequently induces ascent over India following the propagation of Rossby waves along the Asian jet. Ding and Wang (2005) showed that the Tokyo-Chicago Express and the SRP can be viewed as local manifestations of the CGT (Chen and Huang 2012; Kosaka et al. 2012;Hong and Lu 2016).
Despite being a prominent mode of northern hemisphere summertime circulation variability, the ECMWF seasonal forecast model struggles to accurately represent some aspects of the CGT. Beverley et al. (2019) showed that the ECMWF model has a generally weak representation of the CGT, along with poor skill for geopotential height in the location of the centres of action of the CGT, including over west-central Asia and northwest Europe. Several potential causes of this weak representation were proposed-these include propagation errors, such as a northward shifted model jet bias, which will impact on Rossby wave propagation in the model, and forcing errors, such as errors in the Rossby wave source (RWS, Sardeshmukh and Hoskins 1988) over west-central Asia and a  Table 2) weak link between ISM precipitation and the extratropical circulation.
In this paper, the following questions will be addressed: • What impact does heating associated with the ISM have on the extratropical circulation, and what is its role in forcing the CGT? • How does the Indian monsoon influence the circulation over Europe? • What are the main reasons for the lack of skill in the representation of the CGT in the ECMWF model?
To answer these questions, a number of different types of experiment will be presented. Firstly, to further understand the evolution and potential causes of the above errors, results from several relaxation experiments are shown, in which the circulation over west-central Asia, northwest Europe and the ISM region is corrected towards reanalysis. Based on results from these experiments, a further set of experiments in the ECMWF model were carried out in which a thermal forcing was applied in the ISM region to explore the northern hemisphere response to variations in heating associated with the monsoon. In addition to experiments in the coupled seasonal forecast model, we also present results using a linear barotropic model. This type of idealised model can be a useful tool to help understand the propagation of Rossby waves under different conditions. For example, Shaman and Tziperman (2007) showed that El Niño can affect the strength of the Indian monsoon via a westward-propagating Rossby wave response, and O' Reilly et al. (2018) demonstrated that the summertime Euro-Atlantic circulation is largely forced by precipitation anomalies in the eastern tropical Pacific, with contributions from both eastward-and westward-propagating Rossby waves.
The rest of the paper is arranged as follows. In Sect. 2 we outline the data, models and methods used, including an overview of the ECMWF and barotropic models, as well as details of the relaxation and thermal forcing techniques. Results from relaxation experiments are presented in Sect. 3, followed by thermal forcing experiments in Sect. 4 and barotropic model experiments in Sect. 5. A summary and the main conclusions can be found in Sect. 6.

Model details
The seasonal forecast model used is version Cy41r1 of the European Centre for Medium-Range Weather Forecasts (ECMWF)'s Integrated Forecasting System model (IFS), coupled to the Nucleus for European Modelling of the Ocean model (NEMO). This version of the atmospheric model is a more recent one than in the former operational model, System 4 (Cy36r4, Molteni et al. 2011), but older than in the current seasonal forecast model, SEAS5 (Cy43r1, Johnson et al. 2019). The horizontal spectral resolution of the atmospheric model (T255) is the same as System 4 and corresponds to a grid length of approximately 80 km . The atmospheric model was run with 60 vertical levels, with the model top at 0.01 hPa , while the ocean model has a resolution of approximately 1 • with 42 vertical levels (Beverley et al. 2019).
The different number of vertical levels used here (60 compared to 91 in Beverley et al. 2019) is a requirement of performing relaxation experiments, as there needs to be the same number of levels in the model run as there are in the reference dataset (ERA-Interim), to enable the relaxation to take place on all levels. This is the only difference between the model used in these experiments and the model analysed in Beverley et al. (2019), and the control integration used here has similar performance to the 91 level model of Beverley et al. (2019).
Experiments were performed using the ECMWF ERA-Interim (Dee et al. 2011) and ORAS4 (Balmaseda et al. 2013) reanalyses for initialisation of the atmosphere and ocean, respectively. They are initialised on 1st May for the period 1981-2014 and run for four months, therefore covering the boreal summer season of June-August (JJA) and much of the ISM season. The relaxation experiments presented here have 25 ensemble members, whereas the thermal forcing experiments have 5 members.

Relaxation technique
In Sect. 3, relaxation experiments in the ECMWF model are performed in order to further understand the CGT mechanism and errors associated with its representation in the model. Relaxation of the atmosphere in a model is a well-established technique for analysing model performance and deficiencies. The technique (also known as nudging) involves relaxing the model temperature, winds and humidity-based variables towards a reference state (usually reanalysis) throughout the length of a forecast in a pre-defined region. It has previously been used to help understand the impact of tropical biases on the extratropical circulation (Ferranti et al. 1990;Jung et al. 2010a), to identify possible sources of predictability (Klinker 1990) as well as to investigate the origin of atmospheric circulation anomalies in the northern hemisphere in specific extreme seasons (Jung et al. 2010b;Douville et al. 2011;Watson et al. 2016).
The relaxation of variable X is carried out through the addition of an extra term, X relax , to the ECMWF model:

3
where S represents the advective tendencies and sources and sinks due to physical parametrisations in the model, and where X relax is of the form: where X old is the model variable being relaxed, X ref is the reference state towards which the model is drawn (in this case ERA-Interim), ( , ) is a function that smooths the transition from relaxed to non-relaxed regions and is the relaxation timescale which has a value of 2.75 h. The relaxation is applied every time step, with X ref updated by linear interpolation between the 6-hourly data from ERA-Interim. Further details about the relaxation technique can be found in the "Appendix" and details of the individual experiments in Sect. 3.

Thermal forcing technique
In Sect. 4, we perform thermal forcing experiments in the ECMWF model in which we impose a heating in a specific region through the addition of an extra temperature tendency, H , at each model time step, so that: where S is the original model forecast temperature tendency. In these experiments, as before, the model is run from the beginning of May for four months, but the heating is only applied from the 1st July onwards. Further details about the definition of the heating term, H , can be found in the "Appendix".
The horizontal structure of the heating which is applied is given by a Gaussian function which is set to zero outside a certain region. The imposed heating has a vertical structure similar to a "typical" tropical convective heating profile, such as those shown in Figure 3 of Schumacher et al. (2004) or Figure 6 of Li et al. (2009). It increases fairly steadily to a peak at around 400 hPa , before dropping away sharply near the tropopause. Figure 2 shows an example heating profile as used in our experiments, for a surface pressure of 1012 hPa . As the heating is applied on model levels, the actual heating applied is concentrated slightly higher in the atmosphere at lower levels due to the adjustment of the model levels to account for the surface orography. However, the impact of this is only seen in the lower part of the atmosphere, and the heating rate at 400 hPa is almost unchanged. (1)

Linear barotropic model
In Sect. 5, an idealised model is employed to further explore the relationship between ISM heating and the extratropical circulation. The barotropic model used integrates the barotropic vorticity equation on a sphere, following Hoskins and Ambrizzi (1993): where u is the rotational velocity field, is the absolute vorticity, is the relative vorticity, F = F + F � is a constant forcing, is a linear damping with a timescale of 10 days and = 2.4 × 10 16 m 4 s −1 is a diffusion coefficient. Here, F is chosen to keep the model stationary in the absence of any additional forcing so that it exactly maintains the basic state: The equation is solved using spectral harmonics with triangular truncation at wavenumber 42 (T42). The model is initialised from the basic state and integrated forward for 50 days. The anomalous forcing, F ′ , is scaled to be small so that the response can be interpreted as linear.

Verifying data
Observational atmospheric monthly mean fields come from ERA-Interim reanalysis (Dee et al. 2011) at 80 km resolution for the period 1981-2014. Observations of precipitation come from the Global Precipitation Climatology Project (GPCP) dataset (Adler et al. 2003). As in Beverley et al. (2019), in order to focus on year-to-year variations, we only use the interannual component of the ERA-Interim and GPCP data. We performed Fourier harmonic analysis on the monthly mean data to remove the long term trend and decadal variations with a period of longer than 8.5 years. For both datasets, the period used is 1981-2014, to align with the model hindcast period. A Student's t test was used to assess the statistical significance of calculated correlation values (Wilks 2011).

Ensemble sampling
For the relaxation experiments in Sect. 3, the ECMWF forecast system provides us with multiple realisations in the form of 25 ensemble members. When analysing model output, a common technique is to use the ensemble mean. However, using the ensemble mean reduces the noise and the contribution of the forced variability is increased relative to the total variability. Therefore, to ensure that we do not mistake noise in the observations for forced variability, we need to compare individual ensemble members from the relaxation experiments to the observations. To do this, as in Beverley et al. (2019), we follow the method of Johnson et al. (2016) and construct many time series of monthly averaged variables by randomly selecting an ensemble member from each year and repeating this until we have 2000 sets of time series, of 34 years each. We are then able to compare the single realisation of the observed system to multiple realisations of the simulated system.

Investigating the CGT mechanism using relaxation experiments
In Beverley et al. (2019) it was shown that regions of reduced geopotential height skill tend to be located in the same regions as the centres of action of the CGT, and that the ECMWF model CGT is generally weaker than in observations. In addition, it was hypothesised by Wang (2005, 2007) Ding and Wang (2005) [hereafter the D&W region]) is an important region in maintaining the CGT. It was also shown in Beverley et al. (2019) (and also in Fig. 3a-c) that errors in geopotential height are present in this region in July and August. Therefore, to investigate whether model errors in the CGT pattern arise as a result of the development of these errors, the first relaxation experiment ("DW_RELAX") relaxed the circulation in a region centred approximately over this area (Table 1). The region chosen is centred slightly to the north of the D&W region, in order to minimise the impact of the relaxation on the ISM circulation. Analysis of this experiment will determine whether correcting the circulation in this area results in an improved model representation of the CGT.
The weak link between ISM precipitation and west-central Asia in the ECMWF model was also identified by Beverley et al. (2019) as a potential source of the weak model representation of the CGT, due to a weak reinforcement of the wave train by monsoon heating Wang 2005, 2007). The second relaxation experiment ("ISM_RELAX") was therefore designed to investigate the impact of errors in the representation of the ISM circulation on the simulation of the CGT and on extratropical skill.
The third relaxation experiment aimed to investigate the impact of the observed atmospheric variability over northwest Europe on the atmospheric evolution elsewhere, as well as the hemispheric impact of any model errors in this region. In the control experiment, an area of reduced model skill first appears over northwest Europe in June, before further areas of reduced skill develop elsewhere in the northern hemisphere. Therefore, to determine whether the errors seen elsewhere (including in the D&W region) occur as a result of errors propagating from northwest Europe, a third experiment (hereafter "NWEUR_RELAX") was carried out in which a region over northwest Europe was relaxed. . This is defined such that a positive value indicates that the skill in the relaxation experiment is higher than in the control, and vice versa, and the boxes represent the relaxation regions used. As expected, in DW_RELAX there are large improvements in skill in and around the relaxation region. Away from this region, improvements are largely confined to central and east Asia, downstream of the relaxation region. There is very little change in skill over the North Pacific, North America or Europe. Indeed, in June and July in particular, there is actually a reduction in skill in some parts of Europe. This implies that the errors that appear over west-central Asia in July and August (Fig. 3b, c) are not the dominant cause of reductions in skill elsewhere in the northern hemisphere in these months.
The change in skill in ISM_RELAX is shown in Fig. 3g-i. The changes in the relaxation region in this experiment are smaller compared to the other relaxation experiments. This is because the overall skill for 200 hPa geopotential height in the tropics tends to be much better than for the extratropics. The overall pattern of skill change in the extratropics is similar to DW_RELAX. There are some large increases in skill to the north of the relaxation region, including in the D&W region where there was poor skill in the control experiment. Away from the relaxation region, the patterns of skill change vary from month-to-month. The skill over Europe in June is largely unchanged and is slightly reduced in July. In August,  Table 1 however, the changes in skill over Europe are largely positive, which may suggest that the ISM is a potential source of skill in this month. There are also large increases in skill over North Africa in June, associated with the monsoondesert mechanism (Rodwell and Hoskins 1996).
When compared to the other two relaxation experiments, it can be seen that relaxing over northwest Europe (Fig. 3j-l) results in a more widespread and larger hemispheric improvement in skill. In NWEUR_RELAX, improvements in skill are made across much of Eurasia. In particular, the skill in west-central Asia (in the D&W region) is improved, suggesting that errors from northwest Europe propagate to this region. However, similar to DW_RELAX, this relaxation has caused a reduction in skill upstream of the relaxation region. In general, the relaxation in NWEUR_RELAX appears to have resulted in a larger improvement in skill across the northern hemisphere than in DW_RELAX and ISM_RELAX, implying that the observed variability over northwest Europe has a larger influence on other parts of the northern hemisphere. This may be associated with the propagation of quasi-stationary Rossby waves from the North Atlantic region along the Asian jet stream (Bollasina and Messori 2018).
A measure of the change in the model representation of the CGT in each of the relaxation experiments compared to the control is shown in the histograms in Fig. 4. For these, the CGT pattern (the correlation between the D&W Index and 200 hPa geopotential height elsewhere) was calculated for each of the 2000 artificial time series for each experiment, and the pattern correlations for these maps between 30 • -70 • N compared to the equivalent from ERA-Interim ( Fig. 1) were calculated. These are shown as green bars for the control and blue bars for the three relaxation experiments: DW_RELAX (Fig. 4a), ISM_RELAX (Fig. 4b) and NWEUR_RELAX (Fig. 4c).
It can be seen that relaxing in west-central Asia has actually resulted in a worsening of the representation of the CGT in the model. The range of pattern correlations in DW_ RELAX is around −0.1 to 0.7, with a much lower median value of 0.35 compared to 0.59 for the control. Given that in this experiment the D&W Index has been corrected to ERA-Interim, this may suggest that this region is not actually forcing the CGT pattern. This is consistent with the findings of Ding and Wang (2007), who hypothesised that variations in Indian summer monsoon (ISM) precipitation reinforce the west-central Asian high and re-energise the further propagation of the wave train but do not force it directly.
From Fig. 4b it can be seen that the representation of the CGT in ISM_RELAX is very similar to the control. The two distributions are very alike, with similar median values (0.55 for ISM_RELAX and 0.59 for the control), suggesting that the simulation of the CGT has not been improved by correcting the monsoon circulation. This suggests that either the ISM is not a driver of the CGT in the ECMWF model, or that the pathway that connects the ISM to the CGT has not been improved by the relaxation.
Contrary to DW_RELAX, relaxing over northwest Europe results in an improvement in the representation of the CGT in the model compared to the control, with a median pattern correlation of 0.67. The distribution is also narrower Fig. 4 Histograms of the pattern correlations between ERA-Interim and the ECMWF model CGT correlations for August for the control (green bars) and for a DW_RELAX, b ISM_RELAX and c NWEUR_RELAX (blue bars on each panel), calculated using 2000 time series created from the 25 ensemble members. The pattern correlations were calculated for all longitudes between 30 • -70 • N. The median value for the control is 0.59, for DW_RELAX is 0.35, for ISM_RELAX is 0.55 and for NWEUR_RELAX is 0.67. In all three cases, the medians of the relaxation experiment distributions are significantly different from the control distribution median value at the 5% level than in the control and DW_RELAX, indicating that more members have a good pattern correlation, and the range of pattern correlations is around 0.35-0.85. The improvements seen in this experiment are largely due to an improvement of the representation of the wave train between northwest Europe and east Asia. Histograms of the pattern correlation for both the eastern and western hemispheres separately have also been produced (not shown), and the median value for the eastern hemisphere (0.77) is much greater than that for the western hemisphere (0.46). This suggests that northwest Europe is more likely to be forcing west-central Asia, rather than the other way round. This wave train was also identified in the composite analysis in Ding and Wang (2007), where a Rossby wave is seen propagating from northwest Europe from eight days before a maximum in the west-central Asian high. This is also in agreement with the results of Di Capua et al. (2020), who used causal effect network analysis to show that a robust link exists on intraseasonal timescales between the North Atlantic Oscillation and their CGT Index (weekly mean 200 hPa geopotential height in the D&W region). We now compare the ensemble correlations of the D&W Index against the other centres of action of the CGT, as well as the relationship between the D&W Index, ISM precipitation and east Asia (EASIA), for each relaxation experiment and for the control for August (Fig. 5). These correlations were calculated using the 2000 time series created from the 25 ensemble members. The CGT regions used are the same as those defined in Beverley et al. (2019), based on the August observed CGT pattern and are shown as boxes on Fig. 1 and are listed in Table 2. A common theme between the relaxation experiments for the correlations between the D&W Index and all other centres of action is that DW_RELAX is poorer than NWEUR_ RELAX and the control. This is perhaps unexpected, given that in DW_RELAX the D&W Index is the same (or nearly the same) as in ERA-Interim. What this may imply is that while the D&W region is an important centre of action, the CGT is not forced directly from this region in the ECMWF model, and the improvements to the representation of the CGT seen in NWEUR_RELAX suggest that errors over northwest Europe are much more important. This is further emphasised when looking at the correlation between the D&W Index and NWEUR. In NWEUR_RELAX, the relationship between these regions is improved relative to the control, which provides further evidence that the correlation between D&W and NWEUR in observations is more likely to be acting from Europe to Asia, rather than the other way round. This result was one motivation for the heating experiments presented in Sect. 4. Given the difficulty in understanding how the ISM, EASIA and D&W region are related to each other, we are not going to focus on the interpretation of these parts of the figure, but we note that relaxing in the D&W region significantly reduces the correlation between the D&W region and NWEUR. This might be because when we force the D&W region the forcing region is relatively close to NWEUR and we may be interfering with the propagation of waves into the relaxation region, which generates errors in the circulation.
It is perhaps surprising that the correlation between the D&W Index and EASIA region in DW_RELAX is further from ERA-Interim than in the control and NWEUR_ RELAX, given that the skill of 200 hPa geopotential height in EASIA in DW_RELAX is increased. A possible explanation for this lies in the relationship between these two regions and ISM precipitation. The correlations between an Indian summer monsoon heating index (ISMH), defined as the precipitation averaged between 20 • -27.5 • N, 70 • -85 • E, and the D&W Index in DW_RELAX are weaker than in both the control and NWEUR_RELAX, suggesting that the link between these two regions has been partially broken by relaxing in the D&W region. However, the link between ISMH and the EASIA Index, also shown on Fig. 5, is improved in DW_RELAX compared to the control and NWEUR_RELAX. This means that ISMH is still driving the variations in the EASIA Index in DW_RELAX, which results in an improvement in geopotential height skill in EASIA in this experiment. However, as the D&W Index explains less of the variance in EASIA, the correlation between the two regions is reduced in DW_RELAX.
It can also be seen that the control and the relaxation experiments have a similar weak representation of the link between the D&W Index and ISMH, which may partially explain the lack of improvement in the representation of the CGT in these experiments. This provided a motivation for more model experiments, in which the relationship between ISM heating and the extratropical circulation is further explored, results from which can be found in the next section.

Thermal forcing experiments
As seen in the previous section, the link between ISM heating and the extratropical circulation (the D&W region) in the ECMWF model is too weak compared to observations. Therefore, in order to test the response of the model to heating in the monsoon region, a set of thermal forcing experiments was carried out. We may not expect the model to respond to the applied heating in the same way as it would in the real world, due to the variety of real-world external factors and the idealised nature of the additional forcing. Nevertheless, these experiments can still be used to provide insight into the relationship between monsoon heating and the extratropical circulation in the model.
In these experiments, a heating in the form of an extra temperature tendency is applied at each model time step (further details in Sect. 2.1.3). As with the relaxation experiments presented in Sect. 3, the model is run from the beginning of May, but the heating is only applied from the 1st July onwards and is then constant until the end of the model run at the end of August. Each experiment has five ensemble members and the model is run for start dates between 1981-2014.
In order to find the most suitable region to apply the heating to, the regression of August global precipitation against the D&W Index was calculated to determine which regions of tropical precipitation (particularly in the vicinity of the ISM) are most closely associated with variations of the D&W Index. This is shown in the upper panel of Fig. 6, which shows that the strongest regression is located over northwest India, indicating that precipitation in this region has a close relationship with variations of geopotential height in the D&W region. This is in agreement with other studies which have shown that ISM precipitation is associated with an anticyclonic anomaly over west-central Asia, which is generated through a Rossby wave response to the diabatic heating (e.g. Wang et al. 2001;Wu and Wang 2002).
We then calculated the correlation between precipitation in a box centred over the area of strong regression seen in Fig. 6a ( 20 • -27.5 • N, 70 • -85 • E, the same region (ISMH) as used in Fig. 5) and 200 hPa geopotential height (Fig. 6b). From this it can be seen that precipitation over northern India is not only associated with variations of the D&W Index, but also the CGT, with a similar correlation pattern to the observed August CGT pattern (Fig. 1). This suggests that heating in this region is strongly related to the strength of the CGT. Therefore, this region was chosen as the location in which to apply the thermal forcing in the ECMWF model.
In these experiments we look at the impact of applying both a positive and a negative heating over the ISM region. The amount of additional heating was specified to have a peak value which is equivalent to approximately 4 mm day −1 of extra precipitation in the positive heating experiment, and vice versa in the negative heating experiment. This corresponds approximately to a two standard deviation change in the D&W Index (Fig. 6a), and equates to a peak thermal forcing of ±2.02 K day −1 .
When calculating the total heating difference in these experiments compared to the control, we have to consider both the heating that is applied directly, and also the diabatic heating associated with precipitation changes. This was calculated using Eqs. 13-15 given in the "Appendix". Using these equations, the heating applied in the positive thermal forcing experiment was calculated to have a peak of around 120 W m −2 , and the combined effect of this imposed thermal forcing and the resulting precipitation changes on the overall differences in heating in the thermal forcing experiments are shown in Fig. 7. These are with respect to the control experiment.
From Figs. 7a and 7b it can be seen that the maximum increase in heating in the positive heating experiment is centred slightly to the west of the region in which the heating was imposed. This is in the region of the largest increases in precipitation (not shown). This precipitation increase is in a similar location to the region of observed positive regression of precipitation against the D&W Index in Fig. 6a.
The increase is also accompanied by a slight decrease in heating further south, over central parts of India and the Western Ghats, associated with a reduction in precipitation in these areas. From the 850 hPa winds it can be seen that the increase in heating in the positive heating experiment is associated with an extension of the monsoon circulation further north, with stronger winds in this region bringing increased moisture.
The heating differences seen in the negative heating experiment (Fig. 7c, d) are slightly stronger in magnitude than for the positive heating experiment. They are also positioned more towards the centre of the forcing region, associated with a reduction in precipitation over central and northern parts of India. There is also a slight increase in heating over the Arabian Sea, to the west of India.
We now look at the monthly mean response of the model across the northern hemisphere. Figure 8a, b show the difference in ensemble mean 200 hPa geopotential height for the positive minus negative heating experiments in July and August. It can be seen that the geopotential height in several of the centres of action of the CGT has been affected by the Fig. 6  heating. The geopotential height in east Asia and the North Pacific, and also to a lesser extent North America, is up to 50 m higher in the positive than in the negative heating experiment by August, and the location of the geopotential height change maxima in east Asia and the North Pacific are close to the centres of action of the CGT (white crosses on Fig 8). The differences over east Asia and the North Pacific are approximately equivalent to those associated with a two standard deviation change in the D&W Index, and the North America anomalies are equivalent to a one standard deviation change (not shown). This suggests that heating associated with the ISM is a driver of the variability in these parts of the northern hemisphere through driving a wave-like response.
However, the response over the North Atlantic and the NWEUR region is very different to the CGT. Over central and northern Europe and the far North Atlantic there are broadly cyclonic anomalies, which are the opposite of what would be expected in the NWEUR region due to the CGT. There is a strong response visible over North Africa and the Mediterranean, as a result of the Rodwell and Hoskins (1996) monsoon-desert mechanism whereby heating associated with the monsoon triggers westward-propagating Rossby waves. The response over North Africa and the Mediterranean is very different to that seen in Fig. 6b. This may suggest that the precipitation in the heating region is a response to forcing from the CGT, rather than a driver of it. This is further supported by the different response in the The difference in 200 hPa geopotential height relative to the standard deviation of ERA-Interim is shown in Fig. 8c, d. The differences seen associated with the monsoon-desert mechanism are equivalent to up to 3 times the observed standard deviation, indicating that these differences represent a large change relative to usual variability. The difference over east Asia is between 1 and 1.5 observed standard deviations, although the changes over the North Pacific, North America and Europe are equivalent to less than 1 standard deviation.
To examine how the response to the thermal forcing develops over time, Fig. 9 shows the difference in 200 hPa geopotential height between the two experiments averaged over five day periods after the heating is applied. In this figure, Day 0 is the first day after the heating is turned on (1st July). In the first 5 days, the only anomalies that are visible occur as part of the monsoon-desert mechanism, with anticyclonic anomalies appearing in and to the west of the heating region, and a weak cyclonic anomaly to the north of this, which grows in magnitude by Days 6-10. By this point, anticyclonic anomalies have also appeared over east Asia and have begun to grow over the North Pacific as an eastward-propagating Rossby wave response develops. The westward-propagating anomalies have also extended across North Africa to the North Atlantic, and by Days  16-20 these anomalies have reached their peak magnitude, from which point they remain almost constant. The east Asia anomalies also reach their peak magnitude by Days 16-20, while the North Pacific and North America nodes continue to grow until Days 21-25. The North America anomalies seem to arise in part from the interference between the westward-and eastward-propagating waves. From Days 21-25 onwards, the response is largely stationary over the Mediterranean, east Asia and the North Pacific, but the node over North America is a little more transient, shifting in location between west and east North America.
The European response is set up within 10-15 days of the heating being applied, and it is clear that the main differences seen over Europe occur not as a result of the propagation of Rossby waves eastwards from the heating region but as a consequence of the westward-propagating Rossby waves. Throughout the period shown in Fig. 9 the differences in the NWEUR CGT centre of action are cyclonic, unlike the anticyclonic differences seen elsewhere in the hemisphere, and they appear to be the opposite of what would be expected for a CGT-like response. To further understand the development of the response to the heating, we now perform experiments in a linear barotropic model, which are presented in the next section.

Barotropic model experiments
We now use the barotropic model described in Sect. 2.2. The anomalous forcing, F ′ , is taken to be the Rossby wave source (RWS) calculated using the difference in divergence between the positive thermal forcing experiment and the control (a negative RWS, "BT_ISMN"). The vorticity is taken from the ECMWF model control experiment, and so the RWS equation becomes: where is the absolute vorticity, D is the divergence and v is the divergent part of the wind field. The terms with subscript "control" are taken from the ECMWF control experiment and with "exp" are calculated from the difference in divergence between the positive thermal forcing experiment and the control. The calculated RWS forcing used in this experiment is shown in Fig. 10. The basic state used is the July-August average 200 hPa relative vorticity field from the ECMWF model control experiment. Figure 11 shows the evolution of the streamfunction anomalies in this experiment every other day up to Day 12. There is clear eastward propagation of Rossby waves from the forcing region across the North Pacific to North America. Associated with these are centres of action (positive, anticyclonic streamfunction anomalies) in east Asia and the North Pacific, although these are shifted westwards relative to the location of the centres of action of the CGT (the white crosses on Fig. 11). There is also a westward-propagating signal which influences the circulation over the Mediterranean and North Africa, which is consistent with that seen in the thermal forcing experiments. The wave signal becomes less clear over North America as it interacts with the westward-propagating waves, and this interaction means that there is no clear centre of action over North America. The propagation of the waves in this experiment is highly consistent with the results from the thermal forcing experiments shown earlier, albeit with the slight westward shift of the anomalies. The evolution of the signal over Europe in these experiments is also similar to the thermal forcing experiments, whereby the circulation anomalies in this region appear to be influenced by a westward-or northwestwardpropagating signal from the forcing region, rather than from the west via a CGT-like response.
To explore the development of the European response in greater detail, we now follow the method of Shaman andTziperman (2007) andO'Reilly et al. (2018) and add sponge layers at specific longitudes which damp the nearby vorticity anomalies and prevent the zonal propagation of Rossby waves. This allows us to determine the contribution of the westward-and eastward-propagating waves to the development of vorticity anomalies in the northern hemisphere.
The sponge layers, which have a Gaussian longitudinal profile, were implemented through the addition of a term on the right-hand side of the vorticity equation (Eq. 4) of the form: where Here, is the longitude and DAMP is the longitude at which the damping is centred. The first of these sponge layers was added at 60 • W, from pole-to-pole, and the evolution of the vorticity anomalies in this experiment ("BT_NAM_DAMP") and the equivalent experiment with no damping (BT_ISMN) are shown in Fig. 12 (left and centre columns). It can be seen that the midlatitude response in these two experiments develops in a very similar manner, with an obvious wave train propagating eastwards from the forcing region. The westward propagating signals in these two experiments are also extremely similar, such that in the day 40-50 average the anomalies over Europe are almost identical in the two experiments, providing further evidence that the response in this region to forcing from the ISM occurs via westwardpropagating Rossby waves. These results are similar to those of Shaman (2014), who found that the response over Europe to ENSO-related forcing in late summer (JAS) is dominated by a westward-propagating signal, and also the results from the tropical Pacific forcings of O'Reilly et al. (2018), in which a similar signal was observed. Stephan et al. (2019) proposed a positive feedback loop between ascent over India and descent over the Mediterranean, whereby the ISM-induced Mediterranean descent (which occurs through the monsoon-desert mechanism) can subsequently induce rising motion over South Asia following the propagation of anomalies along the Asian jet, associated with the Silk Road Pattern (Lu et al. 2002;Enomoto et al. 2003). To investigate this mechanism, and its relationship with the development of downstream anomalies over North America, a further experiment was carried out with the sponge layer centred at 30 • E ("BT_EUR_DAMP"). The vorticity anomalies in this experiment are shown in the righthand column of Fig. 12. Comparing the damped experiment to the undamped experiment, the effect of preventing the westward-propagating waves can be clearly seen. Early in the simulation, the eastward-propagating response is similar in both experiments. However, by Day 40-50, the anomalies in BT_EUR_DAMP over the North Pacific and North America are weaker. This suggests that the propagation of waves between Asia and North America, associated with the "Tokyo-Chicago Express" wave train (Lau and Weng 2002;Lau et al. 2004), does occur when the westward-propagating waves are damped, but they are much weaker. This agrees with the hypothesis of Stephan et al. (2019), and implies that the westward-moving Rossby wave response is crucial in amplifying the subsequent eastward propagation of the wave train between Asia and North America. This is also similar to the results of Di Capua et al. (2020), who showed that rainfall in the monsoon trough region does exert a forcing on the CGT (the D&W region), although this link is weaker than the reverse. The interference between the eastwardand westward-propagating signals over North America may also be important in the development of the North America anomalies, and this agrees with results from the ECMWF thermal forcing experiments.
To further understand this mechanism, we now examine the response near the forcing region in BT_EUR_DAMP in greater detail. Figure 13 shows the vorticity and wind anomalies in BT_ISMN and BT_EUR_DAMP, focussing this time near the forcing region in the first few days of the simulation. Fig. 12, the vorticity and wind anomalies early in the simulation are very similar in the two experiments, with little difference between them at hour 42. However, from this point on, the anomalies in the two experiments begin to differ. In BT_ISMN (left-hand column of Fig. 13) a strong westward-propagating negative vorticity anomaly is able to develop, and this extends across the northern Arabian Sea and the Middle East by hour 138. In BT_EUR_DAMP, these anomalies are prevented from extending much beyond the eastern parts of the Middle East. As a consequence of this, there is a much greater southerly component to the anomalous wind over the Middle East in the damped experiment, and these southerly wind anomalies will advect low vorticity air to higher latitudes. This prevents the development of the area of positive vorticity anomalies near the Caspian Sea that can be seen in the experiment with no damping, which are hypothesised to be important in the subsequent downstream reinforcement of the wave train through the propagation of these anomalies along the Asian jet stream. When no damping is applied, the negative vorticity anomalies associated with the westward-moving response effectively cut off the anomalous southerlies, allowing the positive vorticity anomalies to develop to the north.

Conclusions
The circumglobal teleconnection (CGT) is a wavenumber-5 wave train with centres of action over west-central Asia, east Asia, the North Pacific, North America and northwest Europe (Fig. 1). Here, we further examined the representation of the CGT in the ECMWF seasonal forecast model, and also the mechanisms linking the Indian summer monsoon (ISM) and the extratropical circulation. This work was performed with the aim of answering the following questions, as posed in Sect. 1: • What impact does heating associated with the ISM have on the extratropical circulation, and what is its role in forcing the CGT?
To explore the impact of ISM heating on the northern hemisphere circulation, two thermal forcing experiments in the ECMWF model were carried out. In these, forcing over the ISM region leads to a strong signal over North Africa, associated with the monsoon-desert mechanism (Rodwell and Hoskins 1996), and also a clear eastwardpropagating wave train response, with areas of increased geopotential height over east Asia, the North Pacific and North America associated with increased monsoon heating, with maxima in similar locations to the centres of action of the CGT in these regions. However, the wave train does not project eastwards into the North Atlantic/ Europe, suggesting that the CGT centre of action in this region is not forced directly by the ISM. The anomalies over northern Europe are cyclonic, and appear to result from a westward-and northwestward-propagating signal from the ISM region. These anomalies develop before the eastward-propagating Rossby waves have travelled around the hemisphere.
Barotropic model experiments, forced using a Rossby wave source applied over India, show a response which is very similar to that from the full ECMWF model, although with a slight westward shift of the circulation anomalies, which supports the idea that the mechanism is through Rossby wave propagation from the heating region. These experiments have similar eastward-and westward-propagating signals, and the development of the response over North America through the interference of these is also very similar in the thermal forcing and barotropic model experiments. • How does the Indian monsoon influence the circulation over Europe? Following Shaman and Tziperman (2007) and O'Reilly et al. (2018), further barotropic model experiments were carried out in which sponge layers were added at certain longitudes to prevent the zonal propagation of Rossby waves through the layer. When a sponge layer was added at 60 • W, the 40-50 day average vorticity anomalies throughout the midlatitudes were very similar to the equivalent experiment with no damping. In particular, the response over Europe was very similar, which further indicates that the influence of the ISM on Europe occurs largely as a result of westward-propagating Rossby waves from the monsoon region, rather than an eastward-moving wave train. Stephan et al. (2019) hypothesised that the Silk Road Pattern (SRP) wave train relies on a positive feedback loop between ascent over India and descent over the Mediterranean. In order to further investigate this mechanism, and its importance in the propagation of the Rossby wave train across the Pacific to North America, another experiment was carried out in which the damping was applied at 30 • E. In this experiment, the vorticity anomalies over the eastern North Pacific and North America have a weaker magnitude by Day 40-50 when compared to the experiment without damping. When examining the response near the forcing region in closer detail, it was found that in the damped experiment there is a much greater southerly component to the anomalous wind over the Middle East. This acts to advect low vorticity air to higher latitudes and suppresses the positive vorticity anomaly that develops near the Caspian Sea in the undamped experiment, which is important in the subsequent downstream reinforcement of the wave train. This reinforcement occurs through the mechanism proposed by Stephan et al. (2019) as follows. Ascent over India leads to westward-propagating Rossby waves, associated with the Rodwell and Hoskins (1996) monsoondesert mechanism. This in turn results in descent over the eastern Mediterranean, which subsequently induces downstream anomalies which propagate along the Asian jet stream and reinforce ISM precipitation in a positive feedback loop. By adding the damping at 30 • E, we have cut off part of this feedback loop and the downstream wave train is weaker as a result. This suggests that the westward-propagating Rossby wave response to forcing in the ISM region is crucial in the downstream propagation of the SRP/Tokyo-Chicago Express wave train. • What are the main reasons for lack of skill in the representation of the CGT in the ECMWF model? Beverley et al. (2019) showed that ECMWF model errors in geopotential height develop in several of the centres of action of the CGT. Relaxation experiments were carried out to explore different aspects of the CGT mechanism and the relationship between these model errors and the weak representation of the CGT in the ECMWF model. Three experiments were performed, in which the circulation was corrected towards ERA-Interim in each of west-central Asia, the ISM region and northwest Europe in turn.
Of these experiments, the largest improvements in skill in the northern hemisphere were seen in the northwest Europe relaxation, suggesting that this region is a more important source of poor skill in the northern hemisphere than the ISM or west-central Asia. The representation of the CGT was also improved in the northwest Europe relaxation, but was unchanged or slightly worsened in the ISM and west-central Asia relaxations. However, the link between monsoon precipitation and the extratropical circulation was found to be weaker than observed in the control and all relaxation experiments. This may mean that the CGT wave train is not being correctly reinforced by monsoon precipitation, which may impact on its subsequent downstream propagation and limit any potential improvements in its representation in the ECMWF model. Alternatively, Ding et al. (2011) suggested that one explanation of the structure of the CGT may lie in the barotropic instability of the midlatitude basic state, and therefore tropical sources may not be responsible for the excitation of the CGT. They found that this barotropic instability is significant over the North Atlantic region, in approximately the same location as the relaxation region in NWEUR_RELAX, which may explain the improvements seen in this experiment.
Through the utilisation of a range of different types of model experiment, including relaxation and thermal forcing experiments in the ECMWF model and simple damped barotropic model experiments, we have shown that ISM heating drives both an eastward-and westward-propagating response, with the westward response important in reinforcing the wave train between Asia and North America. A potential future avenue of research may lie in the nature of thermal forcing that is applied-in our experiments, the forcing is constant from the start of July onwards, but it would be interesting to repeat the experiments with a shorter forcing (e.g. heating for 1 week) or with pulses of heating, as in the experiments of Branstator (2014). Although not part of a CGT-like wave train, the northern European circulation is still influenced by variations in ISM heating and further investigation is required to determine if this mechanism can be used as a potential source of skill for summer weather in this region.