The response of AUSM to precession forcing and its relation to EASM and EAWM

The response of AUSM to precession forcing and its relation to EASM and EAWM | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article The response of AUSM to precession forcing and its relation to EASM and EAWM Mi Yan, Zhengyu Liu, Bin Wang, Xingwei Kong, Liang Ning, Qin Wen, and 1 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-5092629/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 20 May, 2025 Read the published version in npj Climate and Atmospheric Science → Version 1 posted 8 You are reading this latest preprint version Abstract The East Asian monsoon (EAM) and the Australian monsoon (AUM), which constitute the two subsystems of the Asian-Australian monsoon system, are dynamically linked through cross-equatorial flows while also responding distinctly to external forcings. Although previous studies have examined their relationships across various timescales, their connection at the orbital timescale remains underexplored. Through a series of simulations, we demonstrate that Australian Summer Monsoon (AUSM) precipitation varies out of phase with East Asian Summer Monsoon (EASM) precipitation at the precession timescale due to local insolation changes. At the same timescale, the relationship between AUSM circulation and East Asian winter monsoon (EAWM) circulation is weak. This weak linkage arises from the precession-induced inter-hemispheric insolation gradient, which diminishes cross-equatorial flow and consequently weakens the dynamical links between AUSM and EAWM circulations. Earth and environmental sciences/Climate sciences/Climate change Earth and environmental sciences/Climate sciences/Palaeoclimate Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Introduction The Asian-Australian monsoon system comprises three sub-monsoon systems: the Australian monsoon (AUM), the South Asian monsoon (SAM), and the East Asian monsoon (EAM). Each sub-monsoon system can show distinct changes in response to global climate forcing and sea surface temperature (SST) influences 1 . 2 . Furthermore, the EAM and AUM can interact with each other through cross-equatorial flow around the Maritime Continent 3 . Proxy reconstructions and climate model simulations have demonstrated that the relationship between East Asian summer monsoon (EASM) precipitation and Australian summer monsoon (AUSM) precipitation is not consistent across timescales, due to their differing responses to external forcings and internal variations 4 – 6 . At the precession timescale, there is an ongoing debate about whether the strengths of AUSM and EASM vary “in-phase” 1 , 7 – 8 or “out-of-phase” 9 – 10 . This debate underscores a long-standing uncertainty regarding the causes of changes in AUSM 7 : Are the changes in AUSM driven by local austral summer insolation over Australia or by remote boreal winter insolation over Asia? Theoretically, the seasonality of the precession-induced insolation changes in opposite directions between the Northern Hemisphere (NH) and Southern Hemisphere (SH). For instance, at a minimum precession, the NH and SH receive more insolation for boreal summer (austral winter) and less insolation in boreal winter (astral summer). The increased insolation in boreal summer will lead to a strengthened EASM 11 (Wen et al. 2016), while the decreased insolation in austral summer tends to weaken the AUSM. Therefore, if the monsoon systems are determined predominantly by local insolation forcing, the opposite change of local summer monsoon over NH and SH should force an out-of-phase variation at the precession scale, i.e., a negative (or “out-of-phase”) relationship between the AUSM and the EASM. This precession-induced interhemispheric out-of-phase monsoon intensity was supported by some proxies from the Asian cave speleothem record, the North African dust record, and the South American cave speleothem record. Records from Australia, however, have not yet provided absolutely dated evidence for the out-of-phasing to the NH monsoons 12 . On the other hand, the AUSM is believed to be connected to the East Asian winter monsoon (EAWM) through cross-equatorial flow 13 . A strong EAWM may enhance the cross-equatorial outflow, the convergence over the Australian monsoon region, and consequently, the AUSM precipitation. The reconstructions by Clemens et al. 14 and Sun et al. 15 indicate that the EASM and EAWM variations were in phase at precession timescales, prior to the significant Northern Hemisphere glacial cycles approximately 3 million years ago. This in-phase relationship between the EAWM and the EASM at precession timescales has been observed in some recent modeling studies, which were solely influenced by orbital forcing 11 , 16 – 17 . If the local insolation-induced AUSM variation is dominated by the remote dynamical forcing from the EAWM, we should anticipate a positive (or “in-phase”) relationship between AUSM precipitation and EASM precipitation at precession timescales. However, the sensitivity of the AUSM to the interhemispheric interaction with the EAWM remains a topic of intense debate due to a limited understanding of past variability and conflicting results from simulations and reconstructions 1 , 8 , 18 – 19 . Recently, some new records have been reconstructed to reflect the Australian climate evolution during the past 300–400 kyrs 20 – 22 , which provide opportunities to better understand the EAM-AUM relationship. However, they have not been used to quantitatively investigate the EASM-AUSM relationship at the precession timescale, which, therefore, remains unsolved. We then raise the following questions: What is the relationship between the AUSM and the EASM at the precession timescale forced by the precession? Is the AUSM predominated by the local insolation or dynamically linked to the EAWM? In this study, we investigate the relationship between the EAM and AUM at the precessional timescale over the past 300,000 years (300 kyrs) in a numerical model simulation. Specifically, we mainly discuss the monsoon precipitation in local summer for the summer monsoons’ relationship and the monsoon circulation in austral summer to find out the dynamical link. It was found that local precession-induced insolation plays a dominant role in modulating the relationships between the EAM and the AUM. Local sea surface temperature (SST) and remote insolation may also have modulation effects. We first compare the relationship between the summer monsoon precipitation of AUSM and EASM, along with the possible mechanisms. Then, we investigate the simultaneous dynamic linkage between the AUSM circulation and the EAWM circulation and its possible mechanisms. Thirdly, we discuss the other processes that may have modulation roles and provide some model-model and model-data comparisons. The last section provides the conclusions. Results The out-of-phase relationship between AUSM and EASM Here, we show the time series of the AUSM precipitation and EASM precipitation in the pure orbital experiment, which is only forced by the transient orbital parameters (ORB), along with those in the realistic experiment, which is forced by orbital, greenhouse gases and ice sheet (OGG) (Fig. 1 a, 1 b). The experiments are conducted using the Community Climate System Model version 3 (CCSM3). The details of the experiments can be found in Lu and Liu 23 and the Method section. The simulated EASM precipitation is defined as the area averaged June-July-August (JJA) mean precipitation over northern China (35°N-50°N, 105°E-120°E). The simulated AUSM precipitation is defined as the area averaged December-January-February (DJF) mean precipitation over northern Australia over land (20°S-10°S, 120°E-150°E). Two proxy-based monsoon indices representing EASM 24 and AUSM 21 (Fig. 1 c) are used to evaluate our simulated results. The composite δ 18 O records from Chinese caves reconstructed by Cheng et al. 24 have been widely used to represent the EASM. Sarim et al. 21 collected high-resolution clay mineral assemblages and siliciclastic grain-size records from the International Ocean Discovery Program (IODP) site U1483, in which the Kaolinite/smectite ratios were adopted as the proxy for the AUSM. A high kaolinite/smectite ratio indicates a strong AUSM. Another record representing the AUSM intensity reconstructed by Pei et al. 20 is also adopted. In Pei et al. 20 , the AUSM intensity was represented by the monsoonal discharge (Log(Terr/Ca)) in a composite record from Core SO157-18548 and IODP site U1482. A large Log(Terr/Ca) value indicated a high runoff input by monsoonal discharge, i.e., a stronger AUSM. It is clear that the simulated EASM and AUSM precipitation are predominated with precession signal and vary oppositely (Fig. 1 ). They show significant anti-phase coherence at the precession band with a ~ 180° phase shift. The correlation coefficient at lag-0 is -0.73 at the precession timescale between the filtered simulated EASM precipitation and the AUSM precipitation in the ORB run. The out-of-phase relationship between the AUSM and the EASM is robust in the pure orbital run (ORB, Fig. 1 a) and in the realistic run (OGG, Fig. 1 b). However, there is uncertainty in the reconstructed AUSM-EASM relationship. The AUSM intensity reconstructed by Sarim et al. 21 exhibits a clear anti-phase variation compared to the reconstructed EASM (blue line vs. orange line in Fig. 1 c), while the AUSM intensity reconstructed by Pei et al. 20 demonstrates a roughly in-phase variation with the reconstructed EASM (green line vs. orange line in Fig. 1 c). We then calculated the correlation coefficient at lag-0 between the AUSM reconstructed by Pei et al. 20 and the EASM reconstructed by Cheng et al. 24 at their precession band, which is positive with a value of 0.55 (not significant). It is important to note that while the AUSM indicated by Sarim et al. 21 shows an anti-phase variation relative to the EASM indicated by Cheng et al. 24 , their correlation coefficient at lag-0 on the precession timescale is -0.06, which is not significant. The inconsistency between the two AUSM proxies reflects a similar proxy inconsistency noted in previous works 8 . This highlights the uncertainty of proxy records in accurately representing the AUSM. Here, we endeavor to determine whether the proxies and the inferred monsoon responses can be understood physically within the context of our model. The composite analysis of maps (Fig. 2 ) and seasonal cycles (Fig. S1 ) shows that at the low precession stages, the AUSM precipitation over the land is decreased (Fig. 2 a, Fig. S1 a), whereas the EASM precipitation is increased (Fig. 2 b, Fig. S1 b). This also illustrates an “out-of-phase” relationship between the AUSM and the EASM. Note that the AUSM precipitation over the adjacent ocean is increased. We then suggest that if the reconstructed AUSM represents the signal over the ocean, it could be expected a positive relationship with the EASM, shown, for example, in the relationship between Pei et al. 20 and Cheng et al. 24 . We noticed that, dominated by the JJA negative anomaly, the annual mean precipitation over the adjacent ocean also decreased (Fig. S2c). The opposite precipitation changes over land and ocean near northwest Australia are also presented in a previous simulation work by Liu et al. 1 using another coupled ocean-atmosphere model (FOAM). However, the relative magnitudes of the opposite changes differ potentially due to model dependence. Previous reconstructions have suggested that the tropical precipitation pattern could be explained by an interhemispheric anti-phase pattern at the precession timescale, which was related to the monsoon intensity modulated by the precession-induced inter-hemispheric contrast insolation between NH and SH 25 – 26 . This agrees with our simulated result that the AUSM and EASM precipitation have out-of-phase variations at the precession timescale. Such an out-of-phase relationship indicates that the AUSM and the EASM might be dominated by the local precession-induced insolation, and there might be a weak dynamical linkage between the EAWM and the AUSM. Local insolation dominates the out-of-phase relationship between AUSM and EASM Reconstruction works have proposed that the monsoon evolution at orbital timescales followed the orbital-induced insolation closely 27 . Our previous simulation works also revealed the important role of the local and remote orbital-induced insolation on the Asian and Australian monsoon changes in the Holocene 28 – 29 . The precession-induced insolation tends to warm/cool the NH and SH in the same direction (Fig. 3 c, d). From the energy point of view, the moist static energy (MSE) at low levels, including near surface (Fig. S3) and 850 hPa (Fig. 4 a, d), reduced/increased over most parts of the globe in DJF/JJA in response to the reduced/increased insolation at low precession stages. The reduced MSE over the Australian continent in DJF (Fig. 4 a) directly leads to the reduced AUSM precipitation (Fig. 2 a). In contrast, the increased MSE over the East Asia in JJA (Fig. S3b, Fig. 4 d) leads to the increased EASM precipitation (Fig. 2 b). Note the peculiar increased surface temperature over the northern Australian continent in DJF (Fig. 4 b), and it doesn’t seem to match the increased SLP and the anomalous anti-cyclone wind at 850 hPa (Fig. 4 c, similar to surface wind), and, in turn, the reduced monsoon precipitation (Fig. 2 a). Instead, it is consistent with the minimum MSE (Fig. 4 a, S3a), implying that a dry condition accompanies the warming. Decreased precipitation may lead to a drier surface and warmer temperature 30 . Thus, this land warming might be the result, rather than the cause, of the reduced AUSM precipitation over land. Therefore, we argue that the precession-induced reduction of local insolation directly dominates the reduced AUSM precipitation and the resulting low MSE. Regarding the strengthened EASM at low precession stages, besides the increased MSE (Fig. 4 d, S3b), the surface temperature increases more over the land than over the ocean, leading to increased land-sea thermal contrast (Fig. 4 e) and thus the increased land-sea SLP contrast and, in turn, the strengthened EASM southerly (Fig. 4 f) and northward moisture transport. Such local processes in the EASM region play a dominant role in enhancing the intensity of EASM with enhanced rainfall, especially in northern China. Meanwhile, the northward cross-equatorial flow is also strengthened. This may imply an additional remote impact that enhances the EASM from the SH. Quantification of the simultaneous dynamical link between AUSM and EAWM The composite results revealed by our transient simulation is consistent with the previous results derived from a series of equilibrium simulations forced by different obliquity-precession configurations 31 . In their experiments, the reduction in local summer insolation associated with a “high precession” (their “high precession” indicates NH precession, i.e., our “low precession”) led to a lowered precipitation rate over most monsoon regions of northern Australia, accompanied by an increased precipitation rate over Maritime Continent, even though the EAWM wind towards the equator is enhanced. These modeling simulations lead us to hypothesize that the simultaneous dynamical link between the EAWM and the AUSM can be weakened, or even overwhelmed by the local response to precessional forcing. To quantitatively assess the simultaneous dynamical links between the EAM and the AUM, we defined three circulation indices: the meridional wind over EAM region (30 \(\:^\circ\:\) -40 \(\:^\circ\:\) N, 105 \(\:^\circ\:\) -120 \(\:^\circ\:\) E), the meridional wind over Maritime Continent (5 \(\:^\circ\:\) S-5 \(\:^\circ\:\) N, 100 \(\:^\circ\:\) -110 \(\:^\circ\:\) E) (green boxes in Fig. 2 ) and the zonal wind over AUM region (15 \(\:^\circ\:\) S-5 \(\:^\circ\:\) S, 100 \(\:^\circ\:\) -130 \(\:^\circ\:\) E) (blue boxes in Fig. 2 ). The meridional wind and zonal wind at 850 hPa geopotential height are used. The JJA mean and DJF mean meridional wind over the EAM region are used to represent the EASM and EAWM circulation intensities, respectively (EASM_v850, EAWM_v850). The meridional wind over the Maritime Continent is used as a proxy of cross-equatorial flow. The DJF mean zonal wind over the AUM region is used to illustrate the AUSM circulation intensity (AUSM_u850). Since the DJF mean zonal wind corresponding to the AUSM circulation is a westerly wind, we refer to it as “AUSM westerly” in the following part. The rainfall indices are the above-mentioned summer monsoon precipitation of EASM (EASM_P) and AUSM (AUSM_P, over land only). In this study, the EAWM and the AUSM are regarded as dynamically linked if they both have a positive (or “in-phase”) simultaneous relationship with the cross-equatorial flow. The correlation coefficient at lag-0 and the phase shift at the precession timescale between each index are listed in Table 1 . Table 1 The correlation coefficient at lag-0 and the phase shift between each index. The results are derived from the ORB run. The correlation coefficient between each index is conducted based on their filtered precession band. Index Lag-0 correlation coefficient Phase shift (°) EAWM_v850 vs. Cross-equatorial flow 0.70* 10# AUSM_u850 vs. -Cross-equatorial flow 0.16 -100 AUSM_u850 vs. -EAWM_v850 -0.13 -122 AUSM_P vs. -EAWM_v850 -0.40 138# An asterisk (*) indicates the correlation coefficient being above the 90% confidence level, as determined by applying the Monte Carlo test. A pound (#) indicates the phase in which the coherence square passed the 90% confidence level in the coherence analysis. In boreal winter (DJF), although the EAWM northerly is significantly positively correlated with the cross-equatorial flow, with a near 0 \(\:^\circ\:\) phase shift, the correlation coefficients are small and insignificant between the AUSM westerly and the cross-equatorial flow and between the AUSM westerly and the EAWM northerly. This confirms our hypothesis of the weak dynamical link between the EAWM and the AUSM. This hypothesis is also consistent with the significant nearly out-of-phase coherence between the AUSM_P and the EAWM_v850, along with the insignificant negative correlation coefficient (-0.40 in Table 1 ). This suggests that the AUSM circulation and precipitation might be influenced predominantly by other mechanisms instead of the dynamical process from the NH. In other words, the remote AUSM-EAWM dynamic linkage might be overwhelmed by the local insolation-induced effect for the precessional forcing. To summarize, our model simulations suggest that the dynamic link between the EAWM and the AUSM is weak at the precession timescale (insignificant − 0.13 correlation coefficient in Table 1 ). The relationships between each index still hold after being imposed with GHGs and ice sheets (figures not shown), suggesting limited roles of GHGs and ice sheets on the precessional time scale in our model. The roles of GHGs and ice sheets in modulating the monsoon dynamics are beyond the scope of this paper and will be discussed elsewhere. The hemispheric insolation gradient regulates the interhemispheric dynamical links Although the precession induced insolation cooling/warming is of the same sign, or in the same direction, in both the NH and SH, the magnitude of the change is different across latitudes. At low precession stages, the insolation in DJF reduces more over the SH, especially the mid-low latitude, than over the NH (Fig. 3 c). The SH-NH insolation gradient over mid-low latitude (from 30S to 30N) is reduced by 5%, from 2.27 \(\:\times\:{10}^{-2}\) W/m 2 /km to 2.16 \(\:\times\:{10}^{-2}\) W/m 2 /km. The reduced SH-NH insolation gradient leads to weakened interhemispheric contrast in temperature and SLP over the low latitudes (shading in Fig. 4 c, and Fig. S4). The weakened north-south SLP gradient (Fig. S4c) then weakens the tropical easterly and, in turn, the cross-equatorial flow from NH to SH. As such, the remote dynamical impact of the EAWM on AUSM is weakened. The weakened AUSM westerly and the reduced AUSM precipitation is caused by the local south-north SLP gradient between the Australian continent and the nearby ocean to the north, which is reversed and induces an anomalous easterly over the Australian monsoon region, exporting moisture away from the Australian continent into the nearby ocean to the northeast, decreasing rainfall over land and increasing rainfall over the ocean (vector in Fig. 4 c). In JJA, on the other hand, although the insolation increases globally at low precession stages, the mid-low latitudes over the NH receives more insolation than other regions (Fig. 3 d). The NH-SH insolation gradient over mid-low latitude (from 30S to 30N) increases by 6% from 3.13 \(\:\times\:{10}^{-2}\) W/m 2 /km to 3.32 \(\:\times\:{10}^{-2}\) W/m 2 /km. The increased NH-SH insolation gradient strengthens the SH-NH SLP gradient (Fig. S4d), the northward cross-equatorial flow, and the remote dynamic link between the SH circulation and the EASM. In short, caused by the precession-induced changes of the insolation gradient between SH and NH over mid-low latitude, the dynamical link between the EAWM and the AUSM is weakened or even “cut off” at the precession timescale. Discussions Wen et al. 28 studied the responses of EAWM and AUSM to local and remote orbital forcing during the Holocene in model simulations that apply anomalous insolation forcing in different latitude bands. They found that the local forcing for AUSM was overwhelmed by remote forcing from both NH low latitude and SH high latitude. Our composite results of the seasonal cycle of local and remote insolation and precipitation in monsoon regions (Fig. 5 ) suggest that the AUSM precipitation is influenced by both local (15 \(\:^\circ\:\) S-5 \(\:^\circ\:\) S) and remote extratropical SH (60 \(\:^\circ\:\) S-30 \(\:^\circ\:\) S) insolation (Fig. 5 c). The EASM precipitation is also affected by local (30 \(\:^\circ\:\) N-40 \(\:^\circ\:\) N) and remote extratropical NH (45 \(\:^\circ\:\) N-60 \(\:^\circ\:\) N) insolation, which is consistent with our previous work 29 . We notice that the EASM precipitation peaks in August, two months after the insolation peak; this might be due to the selected region that defines the monsoon precipitation. We use the precipitation over northern China, where the precipitation peaks in August, to represent the EASM precipitation, and the August peak is related to the seasonal migration of the EASM circulation. Here, we note that the impacts of local and remote forcing are assessed by comparing the seasonal cycle of insolation at different latitude bands to that of the precipitation over the monsoon region. The insolation represents the local insolation for EASM (30 \(\:^\circ\:\) N-40 \(\:^\circ\:\) N), and the local insolation for AUSM is represented by the insolation over (15 \(\:^\circ\:\) S-5 \(\:^\circ\:\) S). The remote insolation for EASM is represented by the insolation over (45 \(\:^\circ\:\) N-60 \(\:^\circ\:\) N), and the remote insolation for AUSM is represented by the insolation over (60 \(\:^\circ\:\) S-30 \(\:^\circ\:\) S). This assessment may not be accurate because, first, the seasonal cycles of insolation at different latitudes are not completely independent, and second, the monsoon responses to localized insolation forcing may not be linear. Therefore, transient sensitivity experiments with insolation forcing in different latitude bands, similar to Wen et al. 28 , are required in the future to quantitatively isolate the roles of the local and remote insolation on the monsoon evolutions more precisely. By comparing the evolution of global monsoons during Holocene in both a coupled climate model and a coupled AGCM-slab ocean model, Liu et al. 1 found that the direct insolation forcing of the EAWM alone was insufficient to enhance the Australian monsoon. Instead, the local oceanic dynamic response also plays an important role, particularly in AUSM. Here, we cannot explicitly isolate the role of oceanic dynamics, as we do not have a parallel transient simulation in an uncoupled model. Alternatively, we compare the changes in monsoon activities during transition seasons when the orbital-induced insolation is minimal. The variations in precession-induced insolation in March and October are relatively small (Fig. 6 a, b). However, the changes in monsoon circulation and precipitation in March (October) (Fig. S5) resemble those in DJF (JJA) as shown in Fig. 2 . The differing responses between March and October may be related to the observation that the SST in March (October) (Fig. 6 c, d) is similar to that in DJF (JJA) (figures not shown). This confirms the role of SST in modulating the responses of EAM and AUM to orbital-induced insolation, even though these SST effects encompass both thermo and dynamic responses. The precise impacts of SST, and separately, the thermo and dynamical effects on modulating the AUM must be evaluated in future studies through sensitive experiments using AGCM with prescribed SST and AGCM coupled with a slab ocean. To assess the robustness of our model response, we further compared our simulation with two equilibrium experiments of CMIP6/PMIP4, the mid-Holocene (6,000 years before present, 6 ka BP, MH) and the pre-industrial (AD 1950, PI) experiments. The main external forcing during the MH that differs from the PI is the latitudinal and seasonal distribution of insolation caused by the orbital parameters, especially the precession 32 . Therefore, the difference between the MH and the PI can be considered to be a good approximation of the response to the precessional forcing. Six models are randomly selected in this study to conduct the inter-model comparison (Table 2 ). Table 2 CMIP6/PMIP4 models used in this study Model Institution Spatial resolution for atmospheric module Lon*Lat (grids) FGOALS-g3 Institute of Atmospheric Physics, Chinese Academy of Sciences (IAP CAS), China 180*80 IPSL-CM6A-LR Institute Pierre-Simon Laplace (IPSL), France 143*144 MIROC-ES2L Japan Agency for Marine-Earth Science and Technology (JAMSTEC), Atmosphere and Ocean Research Institute (AORI), National Institute for Environmental Studies (NIES), and RIKEN Center for Computational Science (R-CCS), Japan 128*64 MRI-ESM2-0 Meteorological Research Institute (MRI), Japan 320*160 MPI-ESM1-2-LR Max Planck Institute for Meteorology (MPI), Germany 192*96 NESM3 Nanjing University of Information Science and Technology (NUIST), China 192*96 The EASM southerly is stronger and the EASM precipitation is increased in the MH than in PI, robust across the 6 models (Fig. 7 lower panels, Fig. 8 a). However, most models simulate a weaker AUSM precipitation over the northern Australia. This suggests an anti-phase response to insolation between the EASM precipitation and the AUSM precipitation, consistent with our simulations. On the other hand, all the 6 models simulate strengthened EAWM northerly (Fig. 7 upper panels, Fig. 8 b). While 4 out of the 6 models simulate weakened AUSM westerly in MH, indicating a weakened dynamical link between the EAWM and the AUSM, consistent with our simulation results. In comparison with the robust EAM changes across different models, the modeled AUSM here shows some discrepancies from some previous models 1 , 31 . This illustrates some uncertainties in simulating the AUM. The uncertainties might come from the model performance in simulating the above-mentioned local oceanic processes, the dynamical vegetation processes, as well as their competition 8 . For example, all six models simulate opposite summer precipitation responses between the Australian continent changes and its adjacent ocean in the northwest (Fig. 7 upper panels), which is consistent with our simulation (Fig. 2 a) and other previous simulations 1 . It is, therefore, conceivable that differences in the relative importance of land and oceanic rainfall processes may affect the rainfall in the coastal region of AUSM. Improved models and further experiments are required to fully understand the responses of AUM to orbital forcing. Our simulated out-of-phase relationship between the EASM_P and the AUSM_P at the precession timescale (Fig. 1 a, 1 b) is consistent with that between the reconstructed EASM by Cheng et al. 24 and the newly reconstructed AUSM by Sarim et al. 21 (Fig. 1 c). The reconstructed AUSM by Pei et al. 20 , however, doesn’t show anti-phase covariances to the reconstructed EASM by Cheng et al. 24 . We have previously noticed that the AUSM precipitation over the main continent varies oppositely to that over the adjacent ocean at the precession timescale (Fig. 2 a, Fig. 7 upper panels). Therefore, we suggest that the reconstructed AUSM by Pei et al. 20 might indicate the signals over the ocean and that by Sarim et al. 21 might indicate the signals over the land. However, linking the record of Pei et al. 20 with oceanic precipitation here is somewhat arbitrary because if it reflects the signal of land river runoff, it will not be influenced by rainfall over the ocean. Why Pei et al. 20 shows in-phase relationship with Cheng et al. 24 , as indicated by the relationship between the modeled oceanic precipitation and EASM precipitation, requires further investigations in the future. Furthermore, as mentioned above, although the proxy from Sarim et al. 21 shows an anti-phase relationship with Cheng et al. 24 , their correlation coefficient in the precession band is very small and insignificant. We must acknowledge that the proxy records are currently insufficient to support our simulated results. More reconstruction and simulation works are strongly needed to further confirm our hypothesis. The local insolation and the south-north insolation gradient are considered the primary mechanisms influencing the evolutions of the EAM and AUM, and thereby, the relationship and dynamical linkage between the EAM and AUM. However, changes in remote insolation and local oceanic dynamics may modulate AUM variations and, consequently, impact the EAM-AUM relationship and their dynamic connection. This study finds that AUSM precipitation is affected by local and remote insolation from high latitudes in the Southern Hemisphere, which differs slightly from our previous work 28 , which indicated significant remote influences from both low latitudes in the Northern Hemisphere and high latitudes in the Southern Hemisphere. Wen et al. 28 isolated the impacts of local and remote orbital-induced insolation changes on the EASM and AUSM through four experiments that examined anomalous insolation across the latitude ranges of 90N-45N, 45N-0, 0-45S, and 45S-90S. The experiments were conducted in equilibrium, potentially overlooking the lead/lag effects that are expected given the gradual changes in orbital forcing over time 33 . Therefore, further transient sensitivity experiments, following the methodology of Wen et al. 28 to isolate local and remote orbital-induced insolation effects but using transient orbital forcing, are necessary in the future to quantify their influences on AUM evolution at orbital timescales. We should also notice that the correlation coefficient or the coherence does not necessarily illustrate the dynamical links between the EAM and the AUM. The relationship between the NH and SH monsoons mainly indicates their responses to local forcings. More importantly, the summer monsoon precipitation over the Australian land changes differently from that over the adjacent ocean. We should be very careful when comparing the EAM and AUM evolutions at orbital timescales in reconstruction works. Our work may have implications to better understand the monsoon evolutions at orbital timescale. In summary, this study discussed the evolutionary relationships between the AUSM precipitation and the EASM precipitation at the precession timescale, along with the dynamical links between the AUSM and EAWM circulations. Due to the seasonal and the hemispheric asymmetries, the relationships and the links show different responses in different seasons. The precession-induced local insolation change plays a dominant role in deciding the “out-of-phase” relationship between the AUSM and the EASM. And the AUSM - EAWM dynamical linkage can be “weakened” or even “cut off” by the south-north insolation gradient change in response to precession. Future studies are still needed to reduce discrepancies on identifying the relationship between AUSM and EASM/EAWM from both model and proxy perspectives. Methods Simulation Data The experiments used in this study are a set of accelerated transient simulations spanning the past 300,000 years (300 kyrs), conducted with the National Center for Atmospheric Research Community Climate System Model version 3 (NCAR-CCSM3) (Lu and Liu, 2018). The results are derived mainly from the pure orbital run (ORB), which is forced only by orbital parameters, with a 100-acceleration factor. The other experiment (OGG) is also used to evaluate the model performance. It is forced by the transient orbital parameters, the greenhouse gases, and the ice sheet. Proxy Data To evaluate the model results, three proxy-based monsoon indices are collected. The composite δ 18 O from Chinese caves over the last 640 kyrs is used to represent the EASM intensity (Cheng et al., 2016). The kaolinite/smectite ratios from the IODP high sedimentation Site U1483 off the Scott Plateau of the northwest Australian margin are adopted as the proxy for the AUSM (Sarim et al., 2023). Another terrigenous sedimentation along the northwestern Australian margin, the monsoonal discharge (Log(Terr/Ca)) in the composite record from core SO157-18548 and IODP site U1482 (Pei et al., 2021), is also collected to represent AUSM intensity. Declarations Data Availability The accelerated ORB and OGG output are available from the corresponding author upon reasonable request. The CMIP6-PMIP4 output dataset can be accessed through the official CMIP6 channels via the Earth System Grid Federation (ESGF; https://earthsystemcog.org/projects/wip/CMIP6DataRequest). The composite δ 18 O reconstructed by Cheng et al. (2016) is available via the National Centers for Environmental Information (NCEI), National Oceanic and Atmospheric Administration (NOAA; https://www.ncei.noaa.gov/access/paleo-search/). The clay minerals data collected by Sarim et al. (2023) is available via their Zenodo database (https://doi.org/10.5281/zenodo.7418088). The monsoonal discharge (Log(Terr/Ca)) data is downloaded from their online version (https://doi.org/10.1016/j.palaeo.2021.110280). Code availability The source codes for the analysis of this study are available from the corresponding author upon reasonable request. Acknowledgements This work is jointly supported by the Science and Technology Innovation Project of Laoshan Laboratory (Grant LSKJ202203303), the National Natural Science Foundation of China (Grants 42130604, 42475051, 42075049, 41971021), and the Priority Academic Program Development of Jiangsu Higher Education Institutions (Grant 164320H116). Author Contributions M. Y. & Z. L.: Conceptualization, Data, Investigation, Analysis, Writing Original draft, Reviewing, Editing. B. W.: Analysis, Reviewing, Editing. X. K.: Methodology, Visualization. N. L. & Q. W.: Reviewing, Editing. J. L.: Reviewing, Editing, Supervision. Competing interests The authors declare that there are no conflicts of interest regarding the publication of this article. References Liu, Z., B. Otto-Bliesner, J. Kutzbach, et al. (2003). "Coupled climate simulation of the evolution of global monsoon in the Holocene." Journal of Climate 16: 2472–2490. Ding, R., K.-J. Ha, J. Li (2010). Interdecadal shift in the relationship between the East Asian summer monsoon and the tropical Indian Ocean. Climate Dynamics, 34: 1059–1071. Matsumoto, J. (1992). "The Seasonal Changes in Asian and Australian Monsoon Regions." Journal of the Meteorological Society of Japan. Ser. II 70(1B): 257–273. Eroglu, D., F. H. McRobie, I. Ozken, et al. (2016). "See-saw relationship of the Holocene East Asian-Australian summer monsoon." Nat Commun 7: 12929. Zhang, J., M.-Q. Liang, T.-Y. Li, et al. (2022). "Asian-Australian monsoon evolution over the last millennium linked to ENSO in composite stalagmite δ18O records." Quaternary Science Reviews 281: 107420. Xu, C., S. Y. S. Wang, K. Borhara, et al. (2023). "Asian-Australian summer monsoons linkage to ENSO strengthened by global warming." npj Climate and Atmospheric Science 6(1):8. Magee, J. W., G. H. Miller, N. A. Spooner and D. Questiaux (2004). "Continuous 150 k.y. monsoon record from Lake Eyre, Australia: Insolation-forcing implications and unexpected Holocene failure." Geology 32(10): 885. Miller, G., J. Mangan, D. Pollard, et al. (2005). "Sensitivity of the Australian Monsoon to insolation and vegetation: Implications for human impact on continental moisture balance." Geology 33(1): 65. Wang, P. X., B. Wang, H. Cheng, et al. (2014). "The global monsoon across timescales: coherent variability of regional monsoons." Climate of the Past 10(6): 2007–2052. Wang, P. X., B. Wang, H. Cheng, et al. (2017). "The global monsoon across time scales: Mechanisms and outstanding issues." Earth-Science Reviews 174: 84–121. Wen, X., Z. Liu, S. Wang, et al. (2016). "Correlation and anti-correlation of the East Asian summer and winter monsoons during the last 21,000 years." Nat Commun 7: 11999. Müller, U. C., W. J. Fletcher, A. M. Milner and S. Scheiter (2013). "Interhemispheric anti-phasing of orbitally driven monsoon intensity: Implications for ice-volume forcing in the high latitudes." Earth and Planetary Science Letters 377–378: 34–42. Webster, P. J., V. O. Magaña, T. N. Palmer, et al. (1998). "Monsoons: Processes, predictability, and the prospects for prediction." Journal of Geophysical Research: Oceans 103(C7): 14451–14510. Clemens, S. C., and Coauthors, (2008) “Southern Hemisphere forcing of Pliocene δ18O and the evolution of Indo-Asian monsoons.” Paleoceanography, 23, PA4210. Sun, Y. B., and Coauthors, (2010). “Seven million years of wind and precipitation variability on the Chinese Loess Plateau.” Earth and Planetary Science Letters, 297: 525–535. Yan, M., Z. Y. Liu, L. Ning and J. Liu (2020). "Holocene EASM-EAWM Relationship Across Different Timescales in CCSM3." Geophysical Research Letters 47(17): e2020GL088451. Yan, M., Z. Liu, J. Han, et al. (2023). "Relationship between the East Asian Summer and Winter Monsoons at Obliquity Time Scales." Journal of Climate 36(12): 3993–4003. Jalihal, C., J. H. C. Bosmans, J. Srinivasan and A. Chakraborty (2019). "The response of tropical precipitation to Earth's precession: the role of energy fluxes and vertical stability." Climate of the Past 15(2): 449–462. Gong, L., A. Holbourn, W. Kuhnt, et al. (2023). "Middle Pleistocene re-organization of Australian Monsoon." Nature Communications 14(1): 2002. Pei, R., W. Kuhnt, A. Holbourn, et al. (2021). "Monitoring Australian Monsoon variability over the past four glacial cycles." Palaeogeography, Palaeoclimatology, Palaeoecology 568: 110280. Sarim, M., J. Xu, P. Zhang, M. U. Rahman, et al. (2023). "Late quaternary clay mineral and grain-size records from northwest Australia and their implications for paleoclimate, ocean currents, and paleodrainage of the Bonaparte basin." Palaeogeography, Palaeoclimatology, Palaeoecology 610: 111353. Weij, R., J. M. K. Sniderman, J. D. Woodhead, et al. (2024). "Elevated Southern Hemisphere moisture availability during glacial periods." Nature 626(7998): 319–326 Lu, Z. and Z. Liu (2018). "Orbital modulation of ENSO seasonal phase locking." Climate Dynamics 52(7–8): 4329–4350. Cheng, H., R. L. Edwards, A. Sinha, C. Spotl, et al. (2016). "The Asian monsoon over the past 640,000 years and ice age terminations." Nature 534(7609): 640–646. Wang, X., A. S. Auler, R. L. Edwards, et al. (2007). "Millennial-scale precipitation changes in southern Brazil over the past 90,000 years." Geophysical Research Letters 34(23): L23701. Tachikawa, K., O. Cartapanis, L. Vidal, et al. (2011). "The precession phase of hydrological variability in the Western Pacific Warm Pool during the past 400 ka." Quaternary Science Reviews 30(25–26): 3716–3727. Wang, Y., H. Cheng, R. L. Edwards, et al. (2008). "Millennial- and orbital-scale changes in the East Asian monsoon over the past 224,000 years." Nature 451(7182): 1090–1093. Wen, Q., M. Yan, Z. Liu and J. Liu (2022a). "Responses of East Asian Winter Monsoon-Australian Summer Monsoon to Local and Remote Orbital Forcing During Holocene." Geophysical Research Letters 49(10): e2022GL098865. Wen, Q., Z. Liu, J. Zhu, et al. (2022b). "Local Insolation Drives Afro-Asian Monsoon at Orbital‐Scale in Holocene." Geophysical Research Letters 49(6): e2021GL097661. Dai, A., Zhao, T., Chen, J (2018). “Climate Change and Drought: a Precipitation and Evaporation Perspective.” Curr Clim Change Rep, 4, 301–312. Wyrwoll, K.-H., Z. Liu, G. Chen, et al. (2007). "Sensitivity of the Australian summer monsoon to tilt and precession forcing." Quaternary Science Reviews 26(25–28): 3043–3057. Otto-Bliesner, B. L., P. Braconnot, S. P. Harrison, et al. (2017). "The PMIP4 contribution to CMIP6 – Part 2: Two interglacials, scientific objective and experimental design for Holocene and Last Interglacial simulations." Geoscientific Model Development 10(11): 3979–4003. Kutzbach, J. E., X. Liu, Z. Liu and G. Chen (2008). "Simulation of the evolutionary response of global summer monsoons to orbital forcing over the past 280,000 years." Climate Dynamics 30(6): 567–579. Additional Declarations No competing interests reported. Supplementary Files EAMAUMrevision1SI20250418.pdf Cite Share Download PDF Status: Published Journal Publication published 20 May, 2025 Read the published version in npj Climate and Atmospheric Science → Version 1 posted Editorial decision: Revision requested 01 May, 2025 Reviews received at journal 01 May, 2025 Reviews received at journal 26 Apr, 2025 Reviewers agreed at journal 21 Apr, 2025 Reviewers agreed at journal 20 Apr, 2025 Reviewers invited by journal 20 Apr, 2025 Submission checks completed at journal 20 Apr, 2025 First submitted to journal 16 Apr, 2025 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-5092629","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":445613920,"identity":"4967484e-fe7f-4d0f-88be-dc809141c599","order_by":0,"name":"Mi Yan","email":"","orcid":"","institution":"Nanjing Normal University","correspondingAuthor":false,"prefix":"","firstName":"Mi","middleName":"","lastName":"Yan","suffix":""},{"id":445613921,"identity":"0b832046-a5cd-4a36-b95b-4813bf2531d5","order_by":1,"name":"Zhengyu Liu","email":"","orcid":"","institution":"The Ohio State University","correspondingAuthor":false,"prefix":"","firstName":"Zhengyu","middleName":"","lastName":"Liu","suffix":""},{"id":445613922,"identity":"04449b07-79bb-410a-b13c-782a3e5cfb0f","order_by":2,"name":"Bin Wang","email":"","orcid":"","institution":"University of Hawaiʻi at Mānoa","correspondingAuthor":false,"prefix":"","firstName":"Bin","middleName":"","lastName":"Wang","suffix":""},{"id":445613923,"identity":"3c909eaa-cad9-49ae-8a24-ccac11103758","order_by":3,"name":"Xingwei Kong","email":"","orcid":"","institution":"Nanjing Normal University","correspondingAuthor":false,"prefix":"","firstName":"Xingwei","middleName":"","lastName":"Kong","suffix":""},{"id":445613927,"identity":"38ada537-edf9-4758-9b78-43fd8984e954","order_by":4,"name":"Liang Ning","email":"","orcid":"","institution":"Nanjing Normal University","correspondingAuthor":false,"prefix":"","firstName":"Liang","middleName":"","lastName":"Ning","suffix":""},{"id":445613929,"identity":"255019e6-d08c-4dbe-bc76-6cc87d98aa0b","order_by":5,"name":"Qin Wen","email":"","orcid":"","institution":"Nanjing Normal University","correspondingAuthor":false,"prefix":"","firstName":"Qin","middleName":"","lastName":"Wen","suffix":""},{"id":445613930,"identity":"3ee38632-1aa4-41af-bc91-bef138e49b3f","order_by":6,"name":"Jian Liu","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAsElEQVRIiWNgGAWjYBACxnbGBoYPDGxgjgRxWpoZGxhnkKSFgRmIeKBs4rQwNzO3Sdvu4MszOMB88DYPg10eMQ5rk849w1ZscIAt2ZqHIbmYSC1tbIkbDvCYSfMwHEhsIEqLJVgL/zcStDBCbGEjWkuzZS9Qy8zDbMaWcwySCWsxbG9/eONn27HEvuPND2+8qbAjQgtExTFwBDEwGBBSDwTyEKqGCKWjYBSMglEwYgEAjbU1n7MihiMAAAAASUVORK5CYII=","orcid":"","institution":"Nanjing Normal University","correspondingAuthor":true,"prefix":"","firstName":"Jian","middleName":"","lastName":"Liu","suffix":""}],"badges":[],"createdAt":"2024-09-15 12:15:53","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-5092629/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-5092629/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41612-025-01084-y","type":"published","date":"2025-05-20T15:58:35+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":81204764,"identity":"bfd26b2e-1b19-485e-9dae-d8759a459ccb","added_by":"auto","created_at":"2025-04-23 12:02:55","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":552976,"visible":true,"origin":"","legend":"\u003cp\u003eTime series of (a, b, d, e) simulated and (c, f) reconstructed EASM (orange) and AUSM (blue) precipitation over the past 300,000 years (300 kyr). The EASM proxy (orange line in c, f) is the composite δ\u003csup\u003e18\u003c/sup\u003eO records from Chinese caves reconstructed by Cheng et al.\u003csup\u003e24\u003c/sup\u003e. The AUSM proxies are the kaolinite/smectite ratio from IODP site U1483 reconstructed by Sarim et al.\u003csup\u003e21\u003c/sup\u003e (blue line in c, f), and the monsoonal discharge (Log(Terr/Ca)) record integrated from Core SO257-18548 and IODP site U1482 by Pei et al.\u003csup\u003e20\u003c/sup\u003e (green line in c, f). Shown are the time series with 1-kyr resolution (left panels) and in their precession band (right panels). The simulated results are derived from pure orbital run (ORB, a, d) and three forcings run (OGG, b, e).\u003c/p\u003e","description":"","filename":"F1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-5092629/v1/38245d8f441c016aa8c33ebd.jpg"},{"id":81206291,"identity":"5ba7c32a-5903-4256-90a5-dc041a32c7f6","added_by":"auto","created_at":"2025-04-23 12:18:55","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":517771,"visible":true,"origin":"","legend":"\u003cp\u003eThe composite maps of total precipitation rate (PRECT) anomalies in DJF (a) and in JJA (b). The differences are obtained between the low and high precession states (P_low - P_high). Only those values above 95 % confidence level are shown. The color boxes illustrate the regions used to define the indices for the following section. In green boxes, the meridional wind is used. In blue boxes, the zonal wind is used. In red boxes, the precipitation is used.\u003c/p\u003e","description":"","filename":"F2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-5092629/v1/72bb0140c2c11f9825d29772.jpg"},{"id":81205742,"identity":"3502db8f-e0aa-489c-9803-7e9efa037eda","added_by":"auto","created_at":"2025-04-23 12:10:55","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":413482,"visible":true,"origin":"","legend":"\u003cp\u003eComposite maps of zonal mean insolation at the top of atmosphere (SOLIN) climatology at high(solid line) and low (dashed line) precession stages (a, b) and its anomalies (c, d) in DJF (a, c) and JJA (b, d). The composite anomalies are the differences between low and high precession states (L-H).\u003c/p\u003e","description":"","filename":"F3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-5092629/v1/4d45990f5527087d151f6dd7.jpg"},{"id":81205746,"identity":"f8f1a450-c127-42b0-84ff-141a26eda3f8","added_by":"auto","created_at":"2025-04-23 12:10:55","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":860381,"visible":true,"origin":"","legend":"\u003cp\u003eComposite maps of (a, d) moist static energy at 850 hPa (MSE850) and (b, e) surface temperature (TS, shading) along with sea level pressure (SLP, contour), and (c, f) SLP (shading) along with wind at 850 hPa (vector) anomalies in DJF (a-c) and JJA (d-f). The composite anomalies are the differences between low and high precession states. Only those values above 95 % confidence level are shown. The results are derived from ORB run.\u003c/p\u003e","description":"","filename":"F4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-5092629/v1/ab744ed36721bb7544905426.jpg"},{"id":81204770,"identity":"dfc13ac5-d222-48a1-8a8b-c4e924c99a3a","added_by":"auto","created_at":"2025-04-23 12:02:55","extension":"jpg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":303965,"visible":true,"origin":"","legend":"\u003cp\u003eComposite seasonal cycles of insolation at top of atmosphere anomalies (a, b) and precipitation anomalies in monsoon regions (c). The composite anomalies are differences between low and high precession stages. (a) shows the related local insolation for EAM (30°N-40°N) and for AUM (15°S-5°S). (b) shows the related remote insolation for EAM (45°N-60°N) and for AUM (60°S-30°S).\u003c/p\u003e","description":"","filename":"F5.jpg","url":"https://assets-eu.researchsquare.com/files/rs-5092629/v1/e64dd9d1d16880342bc44689.jpg"},{"id":81204766,"identity":"55e1c5a5-a2b1-4f7c-90c6-3f0bdf5833da","added_by":"auto","created_at":"2025-04-23 12:02:55","extension":"jpg","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":631729,"visible":true,"origin":"","legend":"\u003cp\u003eComposite maps of (a, b) SOLIN and (c, f) sea surface temperature (SST) anomalies in March (a, c) and October (b, d). The composite anomalies are the differences between low and high precession states. Only those values above 95 % confidence level are shown.\u003c/p\u003e","description":"","filename":"F6.jpg","url":"https://assets-eu.researchsquare.com/files/rs-5092629/v1/e42cc6d46b6a2c0eadfac22c.jpg"},{"id":81204772,"identity":"d82642de-826f-4d3f-bd43-4502d14811fc","added_by":"auto","created_at":"2025-04-23 12:02:55","extension":"jpg","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":1559929,"visible":true,"origin":"","legend":"\u003cp\u003eThe differences of the precipitation (shading) and wind field at 850 hPa (vector) in DJF (a-f) and JJA (g-l) between MH and PI. The results are derived from the 6 models contributed to CMIP6/PMIP4.\u003c/p\u003e","description":"","filename":"F7.jpg","url":"https://assets-eu.researchsquare.com/files/rs-5092629/v1/cf6cdc7576c0af540b0016bd.jpg"},{"id":81204773,"identity":"f41af262-13a3-4378-9a68-a04695a1fa10","added_by":"auto","created_at":"2025-04-23 12:02:55","extension":"jpg","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":251611,"visible":true,"origin":"","legend":"\u003cp\u003eScatter diagrams of the monsoon precipitation changes (a, EASM_P vs. AUSM_P) and the monsoon circulation changes (b, EAWM_v850 vs. AUSM_u850) between MH and PI across the six CMIP6/PMIP4 models. The horizontal black lines are the zero lines for AUSM. The positive values indicate increasing and/or strengthening, and vice versa.\u003c/p\u003e","description":"","filename":"F8.jpg","url":"https://assets-eu.researchsquare.com/files/rs-5092629/v1/729ff657edf7d2aea4bc9194.jpg"},{"id":83460665,"identity":"90b6ac0d-625f-4dd7-9c3f-ba7d672b4fd8","added_by":"auto","created_at":"2025-05-26 16:13:24","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":5758042,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-5092629/v1/e854c3b9-dd58-4541-b9fe-e0f2811f42b3.pdf"},{"id":81205745,"identity":"5bcf6ca0-8d6e-4c51-951d-27be9afa8311","added_by":"auto","created_at":"2025-04-23 12:10:55","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":1376993,"visible":true,"origin":"","legend":"","description":"","filename":"EAMAUMrevision1SI20250418.pdf","url":"https://assets-eu.researchsquare.com/files/rs-5092629/v1/148bf45eea1c177e2ef62e83.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"The response of AUSM to precession forcing and its relation to EASM and EAWM","fulltext":[{"header":"Introduction","content":"\u003cp\u003eThe Asian-Australian monsoon system comprises three sub-monsoon systems: the Australian monsoon (AUM), the South Asian monsoon (SAM), and the East Asian monsoon (EAM). Each sub-monsoon system can show distinct changes in response to global climate forcing and sea surface temperature (SST) influences\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e.\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e. Furthermore, the EAM and AUM can interact with each other through cross-equatorial flow around the Maritime Continent\u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eProxy reconstructions and climate model simulations have demonstrated that the relationship between East Asian summer monsoon (EASM) precipitation and Australian summer monsoon (AUSM) precipitation is not consistent across timescales, due to their differing responses to external forcings and internal variations\u003csup\u003e\u003cspan additionalcitationids=\"CR5\" citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e\u003c/sup\u003e. At the precession timescale, there is an ongoing debate about whether the strengths of AUSM and EASM vary \u0026ldquo;in-phase\u0026rdquo;\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e,\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u003c/sup\u003e or \u0026ldquo;out-of-phase\u0026rdquo;\u003csup\u003e\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e\u003c/sup\u003e. This debate underscores a long-standing uncertainty regarding the causes of changes in AUSM\u003csup\u003e\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e\u003c/sup\u003e: Are the changes in AUSM driven by local austral summer insolation over Australia or by remote boreal winter insolation over Asia?\u003c/p\u003e \u003cp\u003eTheoretically, the seasonality of the precession-induced insolation changes in opposite directions between the Northern Hemisphere (NH) and Southern Hemisphere (SH). For instance, at a minimum precession, the NH and SH receive more insolation for boreal summer (austral winter) and less insolation in boreal winter (astral summer). The increased insolation in boreal summer will lead to a strengthened EASM\u003csup\u003e\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e\u003c/sup\u003e (Wen et al. 2016), while the decreased insolation in austral summer tends to weaken the AUSM. Therefore, if the monsoon systems are determined predominantly by local insolation forcing, the opposite change of local summer monsoon over NH and SH should force an out-of-phase variation at the precession scale, i.e., a negative (or \u0026ldquo;out-of-phase\u0026rdquo;) relationship between the AUSM and the EASM. This precession-induced interhemispheric out-of-phase monsoon intensity was supported by some proxies from the Asian cave speleothem record, the North African dust record, and the South American cave speleothem record. Records from Australia, however, have not yet provided absolutely dated evidence for the out-of-phasing to the NH monsoons\u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eOn the other hand, the AUSM is believed to be connected to the East Asian winter monsoon (EAWM) through cross-equatorial flow\u003csup\u003e\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u003c/sup\u003e. A strong EAWM may enhance the cross-equatorial outflow, the convergence over the Australian monsoon region, and consequently, the AUSM precipitation. The reconstructions by Clemens et al.\u003csup\u003e\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e\u003c/sup\u003e and Sun et al.\u003csup\u003e\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e\u003c/sup\u003e indicate that the EASM and EAWM variations were in phase at precession timescales, prior to the significant Northern Hemisphere glacial cycles approximately 3\u0026nbsp;million years ago. This in-phase relationship between the EAWM and the EASM at precession timescales has been observed in some recent modeling studies, which were solely influenced by orbital forcing\u003csup\u003e\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e,\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e\u003c/sup\u003e. If the local insolation-induced AUSM variation is dominated by the remote dynamical forcing from the EAWM, we should anticipate a positive (or \u0026ldquo;in-phase\u0026rdquo;) relationship between AUSM precipitation and EASM precipitation at precession timescales. However, the sensitivity of the AUSM to the interhemispheric interaction with the EAWM remains a topic of intense debate due to a limited understanding of past variability and conflicting results from simulations and reconstructions\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e,\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e,\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eRecently, some new records have been reconstructed to reflect the Australian climate evolution during the past 300\u0026ndash;400 kyrs\u003csup\u003e\u003cspan additionalcitationids=\"CR21\" citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u003c/sup\u003e, which provide opportunities to better understand the EAM-AUM relationship. However, they have not been used to quantitatively investigate the EASM-AUSM relationship at the precession timescale, which, therefore, remains unsolved. We then raise the following questions: What is the relationship between the AUSM and the EASM at the precession timescale forced by the precession? Is the AUSM predominated by the local insolation or dynamically linked to the EAWM? In this study, we investigate the relationship between the EAM and AUM at the precessional timescale over the past 300,000 years (300 kyrs) in a numerical model simulation. Specifically, we mainly discuss the monsoon precipitation in local summer for the summer monsoons\u0026rsquo; relationship and the monsoon circulation in austral summer to find out the dynamical link. It was found that local precession-induced insolation plays a dominant role in modulating the relationships between the EAM and the AUM. Local sea surface temperature (SST) and remote insolation may also have modulation effects. We first compare the relationship between the summer monsoon precipitation of AUSM and EASM, along with the possible mechanisms. Then, we investigate the simultaneous dynamic linkage between the AUSM circulation and the EAWM circulation and its possible mechanisms. Thirdly, we discuss the other processes that may have modulation roles and provide some model-model and model-data comparisons. The last section provides the conclusions.\u003c/p\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eThe out-of-phase relationship between AUSM and EASM\u003c/h2\u003e \u003cp\u003eHere, we show the time series of the AUSM precipitation and EASM precipitation in the pure orbital experiment, which is only forced by the transient orbital parameters (ORB), along with those in the realistic experiment, which is forced by orbital, greenhouse gases and ice sheet (OGG) (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ea, \u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb). The experiments are conducted using the Community Climate System Model version 3 (CCSM3). The details of the experiments can be found in Lu and Liu\u003csup\u003e\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e\u003c/sup\u003e and the Method section. The simulated EASM precipitation is defined as the area averaged June-July-August (JJA) mean precipitation over northern China (35\u0026deg;N-50\u0026deg;N, 105\u0026deg;E-120\u0026deg;E). The simulated AUSM precipitation is defined as the area averaged December-January-February (DJF) mean precipitation over northern Australia over land (20\u0026deg;S-10\u0026deg;S, 120\u0026deg;E-150\u0026deg;E).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eTwo proxy-based monsoon indices representing EASM\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e and AUSM\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec) are used to evaluate our simulated results. The composite δ\u003csup\u003e\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e\u003c/sup\u003eO records from Chinese caves reconstructed by Cheng et al.\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e have been widely used to represent the EASM. Sarim et al.\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e collected high-resolution clay mineral assemblages and siliciclastic grain-size records from the International Ocean Discovery Program (IODP) site U1483, in which the Kaolinite/smectite ratios were adopted as the proxy for the AUSM. A high kaolinite/smectite ratio indicates a strong AUSM. Another record representing the AUSM intensity reconstructed by Pei et al.\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e is also adopted. In Pei et al.\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e, the AUSM intensity was represented by the monsoonal discharge (Log(Terr/Ca)) in a composite record from Core SO157-18548 and IODP site U1482. A large Log(Terr/Ca) value indicated a high runoff input by monsoonal discharge, i.e., a stronger AUSM.\u003c/p\u003e \u003cp\u003eIt is clear that the simulated EASM and AUSM precipitation are predominated with precession signal and vary oppositely (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). They show significant anti-phase coherence at the precession band with a\u0026thinsp;~\u0026thinsp;180\u0026deg; phase shift. The correlation coefficient at lag-0 is -0.73 at the precession timescale between the filtered simulated EASM precipitation and the AUSM precipitation in the ORB run. The out-of-phase relationship between the AUSM and the EASM is robust in the pure orbital run (ORB, Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ea) and in the realistic run (OGG, Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb).\u003c/p\u003e \u003cp\u003eHowever, there is uncertainty in the reconstructed AUSM-EASM relationship. The AUSM intensity reconstructed by Sarim et al.\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e exhibits a clear anti-phase variation compared to the reconstructed EASM (blue line vs. orange line in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec), while the AUSM intensity reconstructed by Pei et al.\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e demonstrates a roughly in-phase variation with the reconstructed EASM (green line vs. orange line in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec). We then calculated the correlation coefficient at lag-0 between the AUSM reconstructed by Pei et al.\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e and the EASM reconstructed by Cheng et al.\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e at their precession band, which is positive with a value of 0.55 (not significant). It is important to note that while the AUSM indicated by Sarim et al.\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e shows an anti-phase variation relative to the EASM indicated by Cheng et al.\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e, their correlation coefficient at lag-0 on the precession timescale is -0.06, which is not significant. The inconsistency between the two AUSM proxies reflects a similar proxy inconsistency noted in previous works\u003csup\u003e\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u003c/sup\u003e. This highlights the uncertainty of proxy records in accurately representing the AUSM. Here, we endeavor to determine whether the proxies and the inferred monsoon responses can be understood physically within the context of our model.\u003c/p\u003e \u003cp\u003eThe composite analysis of maps (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e) and seasonal cycles (Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e) shows that at the low precession stages, the AUSM precipitation over the land is decreased (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea, Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003ea), whereas the EASM precipitation is increased (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb, Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003eb). This also illustrates an \u0026ldquo;out-of-phase\u0026rdquo; relationship between the AUSM and the EASM. Note that the AUSM precipitation over the adjacent ocean is increased. We then suggest that if the reconstructed AUSM represents the signal over the ocean, it could be expected a positive relationship with the EASM, shown, for example, in the relationship between Pei et al.\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e and Cheng et al.\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e. We noticed that, dominated by the JJA negative anomaly, the annual mean precipitation over the adjacent ocean also decreased (Fig. S2c). The opposite precipitation changes over land and ocean near northwest Australia are also presented in a previous simulation work by Liu et al.\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u003c/sup\u003e using another coupled ocean-atmosphere model (FOAM). However, the relative magnitudes of the opposite changes differ potentially due to model dependence.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003ePrevious reconstructions have suggested that the tropical precipitation pattern could be explained by an interhemispheric anti-phase pattern at the precession timescale, which was related to the monsoon intensity modulated by the precession-induced inter-hemispheric contrast insolation between NH and SH\u003csup\u003e\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e\u003c/sup\u003e. This agrees with our simulated result that the AUSM and EASM precipitation have out-of-phase variations at the precession timescale. Such an out-of-phase relationship indicates that the AUSM and the EASM might be dominated by the local precession-induced insolation, and there might be a weak dynamical linkage between the EAWM and the AUSM.\u003c/p\u003e \u003c/div\u003e\n\u003ch3\u003eLocal insolation dominates the out-of-phase relationship between AUSM and EASM\u003c/h3\u003e\n\u003cp\u003eReconstruction works have proposed that the monsoon evolution at orbital timescales followed the orbital-induced insolation closely\u003csup\u003e\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e\u003c/sup\u003e. Our previous simulation works also revealed the important role of the local and remote orbital-induced insolation on the Asian and Australian monsoon changes in the Holocene\u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eThe precession-induced insolation tends to warm/cool the NH and SH in the same direction (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ec, d). From the energy point of view, the moist static energy (MSE) at low levels, including near surface (Fig. S3) and 850 hPa (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea, d), reduced/increased over most parts of the globe in DJF/JJA in response to the reduced/increased insolation at low precession stages. The reduced MSE over the Australian continent in DJF (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea) directly leads to the reduced AUSM precipitation (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea). In contrast, the increased MSE over the East Asia in JJA (Fig. S3b, Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ed) leads to the increased EASM precipitation (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eNote the peculiar increased surface temperature over the northern Australian continent in DJF (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eb), and it doesn\u0026rsquo;t seem to match the increased SLP and the anomalous anti-cyclone wind at 850 hPa (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ec, similar to surface wind), and, in turn, the reduced monsoon precipitation (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea). Instead, it is consistent with the minimum MSE (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea, S3a), implying that a dry condition accompanies the warming. Decreased precipitation may lead to a drier surface and warmer temperature\u003csup\u003e\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e\u003c/sup\u003e. Thus, this land warming might be the result, rather than the cause, of the reduced AUSM precipitation over land. Therefore, we argue that the precession-induced reduction of local insolation directly dominates the reduced AUSM precipitation and the resulting low MSE.\u003c/p\u003e \u003cp\u003eRegarding the strengthened EASM at low precession stages, besides the increased MSE (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ed, S3b), the surface temperature increases more over the land than over the ocean, leading to increased land-sea thermal contrast (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ee) and thus the increased land-sea SLP contrast and, in turn, the strengthened EASM southerly (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ef) and northward moisture transport. Such local processes in the EASM region play a dominant role in enhancing the intensity of EASM with enhanced rainfall, especially in northern China. Meanwhile, the northward cross-equatorial flow is also strengthened. This may imply an additional remote impact that enhances the EASM from the SH.\u003c/p\u003e\n\u003ch3\u003eQuantification of the simultaneous dynamical link between AUSM and EAWM\u003c/h3\u003e\n\u003cp\u003eThe composite results revealed by our transient simulation is consistent with the previous results derived from a series of equilibrium simulations forced by different obliquity-precession configurations\u003csup\u003e\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e\u003c/sup\u003e. In their experiments, the reduction in local summer insolation associated with a \u0026ldquo;high precession\u0026rdquo; (their \u0026ldquo;high precession\u0026rdquo; indicates NH precession, i.e., our \u0026ldquo;low precession\u0026rdquo;) led to a lowered precipitation rate over most monsoon regions of northern Australia, accompanied by an increased precipitation rate over Maritime Continent, even though the EAWM wind towards the equator is enhanced. These modeling simulations lead us to hypothesize that the simultaneous dynamical link between the EAWM and the AUSM can be weakened, or even overwhelmed by the local response to precessional forcing.\u003c/p\u003e \u003cp\u003eTo quantitatively assess the simultaneous dynamical links between the EAM and the AUM, we defined three circulation indices: the meridional wind over EAM region (30\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003e-40\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eN, 105\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003e-120\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eE), the meridional wind over Maritime Continent (5\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eS-5\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eN, 100\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003e-110\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eE) (green boxes in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e) and the zonal wind over AUM region (15\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eS-5\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eS, 100\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003e-130\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eE) (blue boxes in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). The meridional wind and zonal wind at 850 hPa geopotential height are used. The JJA mean and DJF mean meridional wind over the EAM region are used to represent the EASM and EAWM circulation intensities, respectively (EASM_v850, EAWM_v850). The meridional wind over the Maritime Continent is used as a proxy of cross-equatorial flow. The DJF mean zonal wind over the AUM region is used to illustrate the AUSM circulation intensity (AUSM_u850). Since the DJF mean zonal wind corresponding to the AUSM circulation is a westerly wind, we refer to it as \u0026ldquo;AUSM westerly\u0026rdquo; in the following part. The rainfall indices are the above-mentioned summer monsoon precipitation of EASM (EASM_P) and AUSM (AUSM_P, over land only).\u003c/p\u003e \u003cp\u003eIn this study, the EAWM and the AUSM are regarded as dynamically linked if they both have a positive (or \u0026ldquo;in-phase\u0026rdquo;) simultaneous relationship with the cross-equatorial flow. The correlation coefficient at lag-0 and the phase shift at the precession timescale between each index are listed in Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eThe correlation coefficient at lag-0 and the phase shift between each index. The results are derived from the ORB run. The correlation coefficient between each index is conducted based on their filtered precession band.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"3\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eIndex\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eLag-0 correlation coefficient\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003ePhase shift (\u0026deg;)\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eEAWM_v850 vs. Cross-equatorial flow\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.70*\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e10#\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAUSM_u850 vs. -Cross-equatorial flow\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.16\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e-100\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAUSM_u850 vs. -EAWM_v850\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e-0.13\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e-122\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAUSM_P vs. -EAWM_v850\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e-0.40\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e138#\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"3\"\u003eAn asterisk (*) indicates the correlation coefficient being above the 90% confidence level, as determined by applying the Monte Carlo test.\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"3\"\u003eA pound (#) indicates the phase in which the coherence square passed the 90% confidence level in the coherence analysis.\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eIn boreal winter (DJF), although the EAWM northerly is significantly positively correlated with the cross-equatorial flow, with a near 0\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003e phase shift, the correlation coefficients are small and insignificant between the AUSM westerly and the cross-equatorial flow and between the AUSM westerly and the EAWM northerly. This confirms our hypothesis of the weak dynamical link between the EAWM and the AUSM.\u003c/p\u003e \u003cp\u003eThis hypothesis is also consistent with the significant nearly out-of-phase coherence between the AUSM_P and the EAWM_v850, along with the insignificant negative correlation coefficient (-0.40 in Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). This suggests that the AUSM circulation and precipitation might be influenced predominantly by other mechanisms instead of the dynamical process from the NH. In other words, the remote AUSM-EAWM dynamic linkage might be overwhelmed by the local insolation-induced effect for the precessional forcing.\u003c/p\u003e \u003cp\u003eTo summarize, our model simulations suggest that the dynamic link between the EAWM and the AUSM is weak at the precession timescale (insignificant \u0026minus;\u0026thinsp;0.13 correlation coefficient in Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). The relationships between each index still hold after being imposed with GHGs and ice sheets (figures not shown), suggesting limited roles of GHGs and ice sheets on the precessional time scale in our model. The roles of GHGs and ice sheets in modulating the monsoon dynamics are beyond the scope of this paper and will be discussed elsewhere.\u003c/p\u003e\n\u003ch3\u003eThe hemispheric insolation gradient regulates the interhemispheric dynamical links\u003c/h3\u003e\n\u003cp\u003eAlthough the precession induced insolation cooling/warming is of the same sign, or in the same direction, in both the NH and SH, the magnitude of the change is different across latitudes. At low precession stages, the insolation in DJF reduces more over the SH, especially the mid-low latitude, than over the NH (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ec). The SH-NH insolation gradient over mid-low latitude (from 30S to 30N) is reduced by 5%, from 2.27 \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\times\\:{10}^{-2}\\)\u003c/span\u003e\u003c/span\u003e W/m\u003csup\u003e2\u003c/sup\u003e/km to 2.16 \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\times\\:{10}^{-2}\\)\u003c/span\u003e\u003c/span\u003e W/m\u003csup\u003e2\u003c/sup\u003e/km. The reduced SH-NH insolation gradient leads to weakened interhemispheric contrast in temperature and SLP over the low latitudes (shading in Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ec, and Fig. S4). The weakened north-south SLP gradient (Fig. S4c) then weakens the tropical easterly and, in turn, the cross-equatorial flow from NH to SH. As such, the remote dynamical impact of the EAWM on AUSM is weakened. The weakened AUSM westerly and the reduced AUSM precipitation is caused by the local south-north SLP gradient between the Australian continent and the nearby ocean to the north, which is reversed and induces an anomalous easterly over the Australian monsoon region, exporting moisture away from the Australian continent into the nearby ocean to the northeast, decreasing rainfall over land and increasing rainfall over the ocean (vector in Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ec).\u003c/p\u003e \u003cp\u003eIn JJA, on the other hand, although the insolation increases globally at low precession stages, the mid-low latitudes over the NH receives more insolation than other regions (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ed). The NH-SH insolation gradient over mid-low latitude (from 30S to 30N) increases by 6% from 3.13 \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\times\\:{10}^{-2}\\)\u003c/span\u003e\u003c/span\u003e W/m\u003csup\u003e2\u003c/sup\u003e/km to 3.32 \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\times\\:{10}^{-2}\\)\u003c/span\u003e\u003c/span\u003e W/m\u003csup\u003e2\u003c/sup\u003e/km. The increased NH-SH insolation gradient strengthens the SH-NH SLP gradient (Fig. S4d), the northward cross-equatorial flow, and the remote dynamic link between the SH circulation and the EASM.\u003c/p\u003e \u003cp\u003eIn short, caused by the precession-induced changes of the insolation gradient between SH and NH over mid-low latitude, the dynamical link between the EAWM and the AUSM is weakened or even \u0026ldquo;cut off\u0026rdquo; at the precession timescale.\u003c/p\u003e"},{"header":"Discussions","content":"\u003cp\u003eWen et al.\u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e studied the responses of EAWM and AUSM to local and remote orbital forcing during the Holocene in model simulations that apply anomalous insolation forcing in different latitude bands. They found that the local forcing for AUSM was overwhelmed by remote forcing from both NH low latitude and SH high latitude. Our composite results of the seasonal cycle of local and remote insolation and precipitation in monsoon regions (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e) suggest that the AUSM precipitation is influenced by both local (15\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eS-5\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eS) and remote extratropical SH (60\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eS-30\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eS) insolation (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003ec). The EASM precipitation is also affected by local (30\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eN-40\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eN) and remote extratropical NH (45\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eN-60\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eN) insolation, which is consistent with our previous work\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. We notice that the EASM precipitation peaks in August, two months after the insolation peak; this might be due to the selected region that defines the monsoon precipitation. We use the precipitation over northern China, where the precipitation peaks in August, to represent the EASM precipitation, and the August peak is related to the seasonal migration of the EASM circulation.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eHere, we note that the impacts of local and remote forcing are assessed by comparing the seasonal cycle of insolation at different latitude bands to that of the precipitation over the monsoon region. The insolation represents the local insolation for EASM (30\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eN-40\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eN), and the local insolation for AUSM is represented by the insolation over (15\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eS-5\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eS). The remote insolation for EASM is represented by the insolation over (45\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eN-60\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eN), and the remote insolation for AUSM is represented by the insolation over (60\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eS-30\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:^\\circ\\:\\)\u003c/span\u003e\u003c/span\u003eS). This assessment may not be accurate because, first, the seasonal cycles of insolation at different latitudes are not completely independent, and second, the monsoon responses to localized insolation forcing may not be linear. Therefore, transient sensitivity experiments with insolation forcing in different latitude bands, similar to Wen et al.\u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e, are required in the future to quantitatively isolate the roles of the local and remote insolation on the monsoon evolutions more precisely.\u003c/p\u003e \u003cp\u003eBy comparing the evolution of global monsoons during Holocene in both a coupled climate model and a coupled AGCM-slab ocean model, Liu et al.\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u003c/sup\u003e found that the direct insolation forcing of the EAWM alone was insufficient to enhance the Australian monsoon. Instead, the local oceanic dynamic response also plays an important role, particularly in AUSM.\u003c/p\u003e \u003cp\u003eHere, we cannot explicitly isolate the role of oceanic dynamics, as we do not have a parallel transient simulation in an uncoupled model. Alternatively, we compare the changes in monsoon activities during transition seasons when the orbital-induced insolation is minimal. The variations in precession-induced insolation in March and October are relatively small (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ea, b). However, the changes in monsoon circulation and precipitation in March (October) (Fig. S5) resemble those in DJF (JJA) as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e. The differing responses between March and October may be related to the observation that the SST in March (October) (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ec, d) is similar to that in DJF (JJA) (figures not shown). This confirms the role of SST in modulating the responses of EAM and AUM to orbital-induced insolation, even though these SST effects encompass both thermo and dynamic responses. The precise impacts of SST, and separately, the thermo and dynamical effects on modulating the AUM must be evaluated in future studies through sensitive experiments using AGCM with prescribed SST and AGCM coupled with a slab ocean.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eTo assess the robustness of our model response, we further compared our simulation with two equilibrium experiments of CMIP6/PMIP4, the mid-Holocene (6,000 years before present, 6 ka BP, MH) and the pre-industrial (AD 1950, PI) experiments. The main external forcing during the MH that differs from the PI is the latitudinal and seasonal distribution of insolation caused by the orbital parameters, especially the precession\u003csup\u003e\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e. Therefore, the difference between the MH and the PI can be considered to be a good approximation of the response to the precessional forcing. Six models are randomly selected in this study to conduct the inter-model comparison (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eCMIP6/PMIP4 models used in this study\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"3\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eModel\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eInstitution\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSpatial resolution for atmospheric module\u003c/p\u003e \u003cp\u003eLon*Lat\u003c/p\u003e \u003cp\u003e(grids)\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eFGOALS-g3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eInstitute of Atmospheric Physics, Chinese Academy of Sciences (IAP CAS), China\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e180*80\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eIPSL-CM6A-LR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eInstitute Pierre-Simon Laplace (IPSL), France\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e143*144\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eMIROC-ES2L\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eJapan Agency for Marine-Earth Science and Technology (JAMSTEC), Atmosphere and Ocean Research Institute (AORI), National Institute for Environmental Studies (NIES), and RIKEN Center for Computational Science (R-CCS), Japan\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e128*64\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eMRI-ESM2-0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMeteorological Research Institute (MRI), Japan\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e320*160\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eMPI-ESM1-2-LR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMax Planck Institute for Meteorology (MPI), Germany\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e192*96\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNESM3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eNanjing University of Information Science and Technology (NUIST), China\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e192*96\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eThe EASM southerly is stronger and the EASM precipitation is increased in the MH than in PI, robust across the 6 models (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e lower panels, Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003ea). However, most models simulate a weaker AUSM precipitation over the northern Australia. This suggests an anti-phase response to insolation between the EASM precipitation and the AUSM precipitation, consistent with our simulations. On the other hand, all the 6 models simulate strengthened EAWM northerly (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e upper panels, Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003eb). While 4 out of the 6 models simulate weakened AUSM westerly in MH, indicating a weakened dynamical link between the EAWM and the AUSM, consistent with our simulation results.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eIn comparison with the robust EAM changes across different models, the modeled AUSM here shows some discrepancies from some previous models\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e,\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e\u003c/sup\u003e. This illustrates some uncertainties in simulating the AUM. The uncertainties might come from the model performance in simulating the above-mentioned local oceanic processes, the dynamical vegetation processes, as well as their competition\u003csup\u003e\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u003c/sup\u003e. For example, all six models simulate opposite summer precipitation responses between the Australian continent changes and its adjacent ocean in the northwest (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e upper panels), which is consistent with our simulation (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea) and other previous simulations\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u003c/sup\u003e. It is, therefore, conceivable that differences in the relative importance of land and oceanic rainfall processes may affect the rainfall in the coastal region of AUSM. Improved models and further experiments are required to fully understand the responses of AUM to orbital forcing.\u003c/p\u003e \u003cp\u003eOur simulated out-of-phase relationship between the EASM_P and the AUSM_P at the precession timescale (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ea, \u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb) is consistent with that between the reconstructed EASM by Cheng et al.\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e and the newly reconstructed AUSM by Sarim et al.\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec). The reconstructed AUSM by Pei et al.\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e, however, doesn\u0026rsquo;t show anti-phase covariances to the reconstructed EASM by Cheng et al.\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e. We have previously noticed that the AUSM precipitation over the main continent varies oppositely to that over the adjacent ocean at the precession timescale (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea, Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e upper panels). Therefore, we suggest that the reconstructed AUSM by Pei et al.\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e might indicate the signals over the ocean and that by Sarim et al.\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e might indicate the signals over the land.\u003c/p\u003e \u003cp\u003eHowever, linking the record of Pei et al.\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e with oceanic precipitation here is somewhat arbitrary because if it reflects the signal of land river runoff, it will not be influenced by rainfall over the ocean. Why Pei et al.\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e shows in-phase relationship with Cheng et al.\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e, as indicated by the relationship between the modeled oceanic precipitation and EASM precipitation, requires further investigations in the future. Furthermore, as mentioned above, although the proxy from Sarim et al.\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e shows an anti-phase relationship with Cheng et al.\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e, their correlation coefficient in the precession band is very small and insignificant. We must acknowledge that the proxy records are currently insufficient to support our simulated results. More reconstruction and simulation works are strongly needed to further confirm our hypothesis.\u003c/p\u003e \u003cp\u003eThe local insolation and the south-north insolation gradient are considered the primary mechanisms influencing the evolutions of the EAM and AUM, and thereby, the relationship and dynamical linkage between the EAM and AUM. However, changes in remote insolation and local oceanic dynamics may modulate AUM variations and, consequently, impact the EAM-AUM relationship and their dynamic connection. This study finds that AUSM precipitation is affected by local and remote insolation from high latitudes in the Southern Hemisphere, which differs slightly from our previous work\u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e, which indicated significant remote influences from both low latitudes in the Northern Hemisphere and high latitudes in the Southern Hemisphere. Wen et al.\u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e isolated the impacts of local and remote orbital-induced insolation changes on the EASM and AUSM through four experiments that examined anomalous insolation across the latitude ranges of 90N-45N, 45N-0, 0-45S, and 45S-90S. The experiments were conducted in equilibrium, potentially overlooking the lead/lag effects that are expected given the gradual changes in orbital forcing over time\u003csup\u003e\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e\u003c/sup\u003e. Therefore, further transient sensitivity experiments, following the methodology of Wen et al.\u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e to isolate local and remote orbital-induced insolation effects but using transient orbital forcing, are necessary in the future to quantify their influences on AUM evolution at orbital timescales.\u003c/p\u003e \u003cp\u003eWe should also notice that the correlation coefficient or the coherence does not necessarily illustrate the dynamical links between the EAM and the AUM. The relationship between the NH and SH monsoons mainly indicates their responses to local forcings. More importantly, the summer monsoon precipitation over the Australian land changes differently from that over the adjacent ocean. We should be very careful when comparing the EAM and AUM evolutions at orbital timescales in reconstruction works. Our work may have implications to better understand the monsoon evolutions at orbital timescale.\u003c/p\u003e \u003cp\u003eIn summary, this study discussed the evolutionary relationships between the AUSM precipitation and the EASM precipitation at the precession timescale, along with the dynamical links between the AUSM and EAWM circulations. Due to the seasonal and the hemispheric asymmetries, the relationships and the links show different responses in different seasons. The precession-induced local insolation change plays a dominant role in deciding the \u0026ldquo;out-of-phase\u0026rdquo; relationship between the AUSM and the EASM. And the AUSM - EAWM dynamical linkage can be \u0026ldquo;weakened\u0026rdquo; or even \u0026ldquo;cut off\u0026rdquo; by the south-north insolation gradient change in response to precession. Future studies are still needed to reduce discrepancies on identifying the relationship between AUSM and EASM/EAWM from both model and proxy perspectives.\u003c/p\u003e "},{"header":"Methods","content":"\u003cdiv id=\"Sec8\" class=\"Section2\"\u003e\u003cdiv id=\"Sec9\" class=\"Section3\"\u003e \u003ch2\u003eSimulation Data\u003c/h2\u003e \u003cp\u003eThe experiments used in this study are a set of accelerated transient simulations spanning the past 300,000 years (300 kyrs), conducted with the National Center for Atmospheric Research Community Climate System Model version 3 (NCAR-CCSM3) (Lu and Liu, 2018). The results are derived mainly from the pure orbital run (ORB), which is forced only by orbital parameters, with a 100-acceleration factor. The other experiment (OGG) is also used to evaluate the model performance. It is forced by the transient orbital parameters, the greenhouse gases, and the ice sheet.\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e\n\u003ch3\u003eProxy Data\u003c/h3\u003e\n\u003cp\u003eTo evaluate the model results, three proxy-based monsoon indices are collected. The composite δ\u003csup\u003e\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e\u003c/sup\u003eO from Chinese caves over the last 640 kyrs is used to represent the EASM intensity (Cheng et al., 2016). The kaolinite/smectite ratios from the IODP high sedimentation Site U1483 off the Scott Plateau of the northwest Australian margin are adopted as the proxy for the AUSM (Sarim et al., 2023). Another terrigenous sedimentation along the northwestern Australian margin, the monsoonal discharge (Log(Terr/Ca)) in the composite record from core SO157-18548 and IODP site U1482 (Pei et al., 2021), is also collected to represent AUSM intensity.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eData Availability\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe accelerated ORB and OGG output are available from the corresponding author upon reasonable request.\u003c/p\u003e\n\u003cp\u003eThe CMIP6-PMIP4 output dataset can be accessed through the official CMIP6 channels via the Earth System Grid Federation (ESGF; https://earthsystemcog.org/projects/wip/CMIP6DataRequest).\u003c/p\u003e\n\u003cp\u003eThe composite δ\u003csup\u003e18\u003c/sup\u003eO reconstructed by Cheng et al. (2016) is available via the National Centers for Environmental Information (NCEI), National Oceanic and Atmospheric Administration (NOAA; https://www.ncei.noaa.gov/access/paleo-search/).\u003c/p\u003e\n\u003cp\u003eThe clay minerals data collected by Sarim et al. (2023) is available via their Zenodo database (https://doi.org/10.5281/zenodo.7418088).\u003c/p\u003e\n\u003cp\u003eThe monsoonal discharge (Log(Terr/Ca)) data is downloaded from their online version (https://doi.org/10.1016/j.palaeo.2021.110280).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCode availability\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe source codes for the analysis of this study are available from the corresponding author upon reasonable request.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work is jointly supported by the Science and Technology Innovation Project of Laoshan Laboratory (Grant LSKJ202203303), the National Natural Science Foundation of China (Grants 42130604, 42475051, 42075049, 41971021), and the Priority Academic Program Development of Jiangsu Higher Education Institutions (Grant 164320H116).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor Contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eM. Y. \u0026amp; Z. L.: Conceptualization, Data, Investigation, Analysis, Writing Original draft, Reviewing, Editing. B. W.: Analysis, Reviewing, Editing. X. K.: Methodology, Visualization. N. L. \u0026amp; Q. W.: Reviewing, Editing. J. L.: Reviewing, Editing, Supervision.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that there are no conflicts of interest regarding the publication of this article.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eLiu, Z., B. Otto-Bliesner, J. Kutzbach, et al. (2003). \"Coupled climate simulation of the evolution of global monsoon in the Holocene.\" Journal of Climate 16: 2472\u0026ndash;2490.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDing, R., K.-J. Ha, J. Li (2010). Interdecadal shift in the relationship between the East Asian summer monsoon and the tropical Indian Ocean. Climate Dynamics, 34: 1059\u0026ndash;1071.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMatsumoto, J. (1992). \"The Seasonal Changes in Asian and Australian Monsoon Regions.\" Journal of the Meteorological Society of Japan. Ser. II 70(1B): 257\u0026ndash;273.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEroglu, D., F. H. McRobie, I. Ozken, et al. (2016). \"See-saw relationship of the Holocene East Asian-Australian summer monsoon.\" Nat Commun 7: 12929.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang, J., M.-Q. Liang, T.-Y. Li, et al. (2022). \"Asian-Australian monsoon evolution over the last millennium linked to ENSO in composite stalagmite δ18O records.\" Quaternary Science Reviews 281: 107420.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXu, C., S. Y. S. Wang, K. Borhara, et al. (2023). \"Asian-Australian summer monsoons linkage to ENSO strengthened by global warming.\" npj Climate and Atmospheric Science 6(1):8.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMagee, J. W., G. H. Miller, N. A. Spooner and D. Questiaux (2004). \"Continuous 150 k.y. monsoon record from Lake Eyre, Australia: Insolation-forcing implications and unexpected Holocene failure.\" \u003cem\u003eGeology\u003c/em\u003e 32(10): 885.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMiller, G., J. Mangan, D. Pollard, et al. (2005). \"Sensitivity of the Australian Monsoon to insolation and vegetation: Implications for human impact on continental moisture balance.\" \u003cem\u003eGeology\u003c/em\u003e 33(1): 65.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang, P. X., B. Wang, H. Cheng, et al. (2014). \"The global monsoon across timescales: coherent variability of regional monsoons.\" Climate of the Past 10(6): 2007\u0026ndash;2052.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang, P. X., B. Wang, H. Cheng, et al. (2017). \"The global monsoon across time scales: Mechanisms and outstanding issues.\" Earth-Science Reviews 174: 84\u0026ndash;121.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWen, X., Z. Liu, S. Wang, et al. (2016). \"Correlation and anti-correlation of the East Asian summer and winter monsoons during the last 21,000 years.\" Nat Commun 7: 11999.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eM\u0026uuml;ller, U. C., W. J. Fletcher, A. M. Milner and S. Scheiter (2013). \"Interhemispheric anti-phasing of orbitally driven monsoon intensity: Implications for ice-volume forcing in the high latitudes.\" Earth and Planetary Science Letters 377\u0026ndash;378: 34\u0026ndash;42.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWebster, P. J., V. O. Maga\u0026ntilde;a, T. N. Palmer, et al. (1998). \"Monsoons: Processes, predictability, and the prospects for prediction.\" Journal of Geophysical Research: Oceans 103(C7): 14451\u0026ndash;14510.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eClemens, S. C., and Coauthors, (2008) \u0026ldquo;Southern Hemisphere forcing of Pliocene δ18O and the evolution of Indo-Asian monsoons.\u0026rdquo; Paleoceanography, 23, PA4210.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSun, Y. B., and Coauthors, (2010). \u0026ldquo;Seven million years of wind and precipitation variability on the Chinese Loess Plateau.\u0026rdquo; Earth and Planetary Science Letters, 297: 525\u0026ndash;535.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYan, M., Z. Y. Liu, L. Ning and J. Liu (2020). \"Holocene EASM-EAWM Relationship Across Different Timescales in CCSM3.\" \u003cem\u003eGeophysical Research Letters\u003c/em\u003e 47(17): e2020GL088451.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYan, M., Z. Liu, J. Han, et al. (2023). \"Relationship between the East Asian Summer and Winter Monsoons at Obliquity Time Scales.\" Journal of Climate 36(12): 3993\u0026ndash;4003.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJalihal, C., J. H. C. Bosmans, J. Srinivasan and A. Chakraborty (2019). \"The response of tropical precipitation to Earth's precession: the role of energy fluxes and vertical stability.\" Climate of the Past 15(2): 449\u0026ndash;462.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGong, L., A. Holbourn, W. Kuhnt, et al. (2023). \"Middle Pleistocene re-organization of Australian Monsoon.\" \u003cem\u003eNature Communications\u003c/em\u003e 14(1): 2002.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePei, R., W. Kuhnt, A. Holbourn, et al. (2021). \"Monitoring Australian Monsoon variability over the past four glacial cycles.\" Palaeogeography, Palaeoclimatology, Palaeoecology 568: 110280.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSarim, M., J. Xu, P. Zhang, M. U. Rahman, et al. (2023). \"Late quaternary clay mineral and grain-size records from northwest Australia and their implications for paleoclimate, ocean currents, and paleodrainage of the Bonaparte basin.\" Palaeogeography, Palaeoclimatology, Palaeoecology 610: 111353.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWeij, R., J. M. K. Sniderman, J. D. Woodhead, et al. (2024). \"Elevated Southern Hemisphere moisture availability during glacial periods.\" Nature 626(7998): 319\u0026ndash;326\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLu, Z. and Z. Liu (2018). \"Orbital modulation of ENSO seasonal phase locking.\" Climate Dynamics 52(7\u0026ndash;8): 4329\u0026ndash;4350.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCheng, H., R. L. Edwards, A. Sinha, C. Spotl, et al. (2016). \"The Asian monsoon over the past 640,000 years and ice age terminations.\" Nature 534(7609): 640\u0026ndash;646.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang, X., A. S. Auler, R. L. Edwards, et al. (2007). \"Millennial-scale precipitation changes in southern Brazil over the past 90,000 years.\" Geophysical Research Letters 34(23): L23701.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTachikawa, K., O. Cartapanis, L. Vidal, et al. (2011). \"The precession phase of hydrological variability in the Western Pacific Warm Pool during the past 400 ka.\" Quaternary Science Reviews 30(25\u0026ndash;26): 3716\u0026ndash;3727.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang, Y., H. Cheng, R. L. Edwards, et al. (2008). \"Millennial- and orbital-scale changes in the East Asian monsoon over the past 224,000 years.\" Nature 451(7182): 1090\u0026ndash;1093.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWen, Q., M. Yan, Z. Liu and J. Liu (2022a). \"Responses of East Asian Winter Monsoon-Australian Summer Monsoon to Local and Remote Orbital Forcing During Holocene.\" \u003cem\u003eGeophysical Research Letters\u003c/em\u003e 49(10): e2022GL098865.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWen, Q., Z. Liu, J. Zhu, et al. (2022b). \"Local Insolation Drives Afro-Asian Monsoon at Orbital‐Scale in Holocene.\" \u003cem\u003eGeophysical Research Letters\u003c/em\u003e 49(6): e2021GL097661.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDai, A., Zhao, T., Chen, J (2018). \u0026ldquo;Climate Change and Drought: a Precipitation and Evaporation Perspective.\u0026rdquo; Curr Clim Change Rep, 4, 301\u0026ndash;312.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWyrwoll, K.-H., Z. Liu, G. Chen, et al. (2007). \"Sensitivity of the Australian summer monsoon to tilt and precession forcing.\" Quaternary Science Reviews 26(25\u0026ndash;28): 3043\u0026ndash;3057.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOtto-Bliesner, B. L., P. Braconnot, S. P. Harrison, et al. (2017). \"The PMIP4 contribution to CMIP6 \u0026ndash; Part 2: Two interglacials, scientific objective and experimental design for Holocene and Last Interglacial simulations.\" Geoscientific Model Development 10(11): 3979\u0026ndash;4003.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKutzbach, J. E., X. Liu, Z. Liu and G. Chen (2008). \"Simulation of the evolutionary response of global summer monsoons to orbital forcing over the past 280,000 years.\" Climate Dynamics 30(6): 567\u0026ndash;579.\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"npj-climate-and-atmospheric-science","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"npjclimatsci","sideBox":"Learn more about [npj Climate and Atmospheric Science](http://www.nature.com/npjclimatsci/)","snPcode":"41612","submissionUrl":"https://submission.springernature.com/new-submission/41612/3","title":"npj Climate and Atmospheric Science","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"NPJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"","lastPublishedDoi":"10.21203/rs.3.rs-5092629/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-5092629/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eThe East Asian monsoon (EAM) and the Australian monsoon (AUM), which constitute the two subsystems of the Asian-Australian monsoon system, are dynamically linked through cross-equatorial flows while also responding distinctly to external forcings. Although previous studies have examined their relationships across various timescales, their connection at the orbital timescale remains underexplored. Through a series of simulations, we demonstrate that Australian Summer Monsoon (AUSM) precipitation varies out of phase with East Asian Summer Monsoon (EASM) precipitation at the precession timescale due to local insolation changes. At the same timescale, the relationship between AUSM circulation and East Asian winter monsoon (EAWM) circulation is weak. This weak linkage arises from the precession-induced inter-hemispheric insolation gradient, which diminishes cross-equatorial flow and consequently weakens the dynamical links between AUSM and EAWM circulations.\u003c/p\u003e","manuscriptTitle":"The response of AUSM to precession forcing and its relation to EASM and EAWM","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-04-23 12:02:50","doi":"10.21203/rs.3.rs-5092629/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2025-05-01T18:16:05+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-05-01T16:41:49+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-04-26T06:18:23+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"193381776509941567969722959516788035636","date":"2025-04-21T07:34:47+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"2465760583735297879328519622108829290","date":"2025-04-21T03:37:21+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2025-04-20T23:13:28+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2025-04-20T12:14:07+00:00","index":"","fulltext":""},{"type":"submitted","content":"npj Climate and Atmospheric Science","date":"2025-04-16T07:07:52+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"npj-climate-and-atmospheric-science","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"npjclimatsci","sideBox":"Learn more about [npj Climate and Atmospheric Science](http://www.nature.com/npjclimatsci/)","snPcode":"41612","submissionUrl":"https://submission.springernature.com/new-submission/41612/3","title":"npj Climate and Atmospheric Science","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"NPJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"54d46ddd-e420-4a6d-aa1b-abedfb998922","owner":[],"postedDate":"April 23rd, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":47435171,"name":"Earth and environmental sciences/Climate sciences/Climate change"},{"id":47435172,"name":"Earth and environmental sciences/Climate sciences/Palaeoclimate"}],"tags":[],"updatedAt":"2025-05-26T16:10:59+00:00","versionOfRecord":{"articleIdentity":"rs-5092629","link":"https://doi.org/10.1038/s41612-025-01084-y","journal":{"identity":"npj-climate-and-atmospheric-science","isVorOnly":false,"title":"npj Climate and Atmospheric Science"},"publishedOn":"2025-05-20 15:58:35","publishedOnDateReadable":"May 20th, 2025"},"versionCreatedAt":"2025-04-23 12:02:50","video":"","vorDoi":"10.1038/s41612-025-01084-y","vorDoiUrl":"https://doi.org/10.1038/s41612-025-01084-y","workflowStages":[]},"version":"v1","identity":"rs-5092629","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-5092629","identity":"rs-5092629","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}