Drivers and mechanisms for enhanced summer monsoon precipitation over East Asia during the mid-Pliocene in the IPSL-CM5A

A comparative analysis of East Asian summer monsoon (EASM) precipitation is performed to reveal the drivers and mechanisms controlling the similarities of the mid-Pliocene EASM precipitation changes compared to the corresponding pre-industrial (PI) experiments derived from atmosphere-only (i.e. AGCM) and fully coupled (i.e. CGCM) simulations, as well as the large simulated differences in the mid-Pliocene EASM precipitation between the two simulations. The area-averaged precipitation over the EASM domain is enhanced in the mid-Pliocene compared to the corresponding PI experiments performed by both the AGCM (LMDZ5A) and the CGCM (IPSL-CM5A). Moisture budget analysis reveals that it is the surface warming over East Asia that drives the area-averaged EASM precipitation increase in the mid-Pliocene in both simulations. The surface warming increases the atmospheric moisture content, as revealed by an increase in the thermodynamic component of vertical moisture advection, resulting in enhanced mid-Pliocene EASM precipitation compared to PI in both simulations. Moist static energy diagnosis identifies the combined effect of enhanced zonal thermal contrast and column-integrated meridional stationary eddy velocity v∗¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{{v^{*} }}$$\end{document} and its convergence ∂v∗¯∂y\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{{\overline{{\partial v^{*} }} }}{\partial y}$$\end{document} as the physical mechanisms that sustain the enhancement of mid-Pliocene EASM precipitation in both simulations compared to the PI experiments. This takes place through a strengthening of the EASM circulation and moisture transport into the EASM domain associated with an increase in local moisture convergence in the mid-Pliocene in both simulations. Moisture budget analysis also reveals that the larger area-averaged mid-Pliocene EASM precipitation increase in the CGCM compared to its AGCM component is mainly caused by the dynamical component contributing more to the vertical moisture advection in the CGCM (i.e. IPSL-CM5A) compared to its AGCM (LMDZ5). The large simulated differences in the spatial pattern of the mid-Pliocene EASM precipitation between the two simulations result from the combined effect of enhanced meridional thermal contrast over the EASM domain and increased v∗¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{{v^{*} }}$$\end{document} convergence over South China in the CGCM simulation compared to the AGCM simulation.


Introduction
The mid-Pliocene warm period (~3.3-3.0 Ma) is widely considered as the most recent warm period in Earth's history when the global average surface air temperature (SAT) Abstract A comparative analysis of East Asian summer monsoon (EASM) precipitation is performed to reveal the drivers and mechanisms controlling the similarities of the mid-Pliocene EASM precipitation changes compared to the corresponding pre-industrial (PI) experiments derived from atmosphere-only (i.e. AGCM) and fully coupled (i.e. CGCM) simulations, as well as the large simulated differences in the mid-Pliocene EASM precipitation between the two simulations. The area-averaged precipitation over the EASM domain is enhanced in the mid-Pliocene compared to the corresponding PI experiments performed by both the AGCM (LMDZ5A) and the CGCM (IPSL-CM5A). Moisture budget analysis reveals that it is the surface warming over East Asia that drives the area-averaged EASM precipitation increase in the mid-Pliocene in both simulations. The surface warming increases the atmospheric moisture content, as revealed by an increase in the thermodynamic was higher than the present day (Haywood et al. 2000;Haywood and Valdes 2004). More importantly, this period was the same as or quite similar in many respects to the modern climate, such as the geographic distribution of the continents and oceans, the atmospheric CO 2 concentration (estimated to have been slightly higher than the present day at 360-440 ppm), and the sea level (estimated to have been 22 ± 10 m higher than the modern level) (Miller et al. 2012). Given the abundance of proxy data and such large similarities with the modern climate, the mid-Pliocene warm period has been recommended as a potential analogue for future climate predictions by the end of the twenty first century (Thompson and Fleming 1996;Jansen et al. 2007;Sun et al. 2013). In particular, Sun et al. (2013) compared the simulated mid-Pliocene climate to that under a Representative Concentration Pathway (RCP) future scenario (RCP4.5) and found that Hadley circulation changes were similar in both simulations, strengthening the idea that the mid-Pliocene could be a useful tool for understanding atmospheric changes in the future.
Given the importance of the mid-Pliocene warm period, data-model comparisons have been performed to better understand the dynamical cause of this warm climate in Earth's history. The PRISM (Pliocene Research, Interpretation and Synoptic Mapping) project was initially devised to reconstruct the sea surface conditions of the mid-Pliocene by utilizing multi-proxy records, the aim being to create the oceanic boundary conditions necessary to drive atmosphere-only models for mid-Pliocene simulation. Recently, the first phase of the Pliocene Modeling Intercomparison Project (PlioMIP1) has been conducted to validate and assess the ability of state-of-the-art climate models to simulate the warm climate of the past . Two types of experiments have been conducted, proposed by the PlioMIP guidelines, to investigate the possible causes of the mid-Pliocene warm climate. The first (Experiment 1) is performed with atmosphere-only climate models (atmospheric general circulation models; AGCMs) to simulate the mid-Pliocene and pre-industrial (PI) climate. The second (Experiment 2) utilizes fully coupled oceanatmosphere climate models [atmosphere-ocean general circulation models (AOGCMs) or coupled general circulation models (CGCMs)] to perform the mid-Pliocene and PI simulations (Haywood et al. 2010(Haywood et al. , 2011. All the PlioMIP1 participating models show good performance in reproducing the global SAT warming in the mid-Pliocene (Chan et al. 2011;Bragg et al. 2012;Contoux et al. 2012;Kamae and Ueda 2012;Koenig et al. 2012;Stepanek and Lohmann 2012;Yan et al. 2012a, b;Chandler et al. 2013;Rosenbloom et al. 2013;Zheng et al. 2013), despite the mechanisms underpinning the trend remaining controversial (Haywood and Valdes 2004;Hill et al. 2014).
East Asian climate during the mid-Pliocene warm period has also been investigated (Yan et al. 2012a, b), revealing a stronger than present East Asian summer monsoon (EASM) using the Community Atmosphere Model version 3.1 (CAM3.1). This result was subsequently confirmed by multi-model comparison analysis, in both atmosphere-only and coupled model experiments from the PlioMIP . However, compared to the EASM, the comparison showed that there is greater uncertainty regarding the East Asian winter monsoon (EAWM) in the mid-Pliocene. Despite the features of the EASM and EAWM during the mid-Pliocene warm period having been revealed well by this multi-model comparison study (Zhang et al. , 2015, a comprehensive comparison of mid-Pliocene EASM precipitation simulated by AGCMs and CGCMs has never been addressed. Changes in EASM precipitation associated with EASM circulation variations have direct and wide societal and economic impacts over East Asia, and thus it is important to investigate the possible factors determining the position of the rain belt and its intensity over East Asia. Sampe and Xie (2010) argued that temperature advection from the southeastern flank of the Tibetan Plateau at 500 hPa is important because of its role in inducing vertical motion and because of its control over the position of the rain belt over East Asia. Moist static energy (MSE) has been proposed to gauge the prominent feature of present day EASM precipitation (Chen and Bordoni 2014a). However, no study has investigated the mechanisms controlling the position of the EASM precipitation during the mid-Pliocene. This is an important topic because it can help provide an insight into the possible status of EASM precipitation in the future, which of course has important social implications.
In this study, an analysis of the similarities and differences of the mid-Pliocene EASM precipitation simulated by an atmosphere-only and a coupled model is carried out. Diagnosis of the moisture budget equation is used to explain the cause of the mid-Pliocene EASM precipitation changes compared to PI in both simulations, as well as the differences between the mid-Pliocene EASM precipitation simulated by the two models.
The paper is organized as follows: The models and experimental design are described in Sect. 2. In Sect. 3, a comparison of the mid-Pliocene EASM precipitation is made using the moisture budget equation. In Sect. 4, we investigate the factors determining the horizontal distribution of the EASM precipitation and then explain the mechanism involved in the simulated mid-Pliocene EASM precipitation changes relative to the PI experiment. The differences between the mid-Pliocene EASM precipitation simulated by the CGCM and AGCM are discussed in Sect. 5. In Sect. 6, the drivers and mechanisms for enhanced mid-Pliocene EASM precipitation, as revealed by IPSL-CM5A-LR and its atmospheric component (LMDZ5), are further confirmed by performing multimodel comparisons. Section 7 provides further discussion and a summary.

Climate model description
The climate model used in the present study is IPSL-CM5A-LR, which is a state-of-the art coupled GCM developed by Institute Pierre Simon Laplace . This model has been widely used for paleoclimate simulations and future projections. IPSL-CM5A-LR is a coupled climate model including four components: the LMDZ5A atmospheric model, which has a horizontal resolution of 1.875° latitude × 3.75° longitude and 39 vertical levels (R96 × 95L39) ); the NEMOv3.2 ocean model, which has a mean grid spacing about 2° (latitudinal resolution of 0.5° near the equator and 1° in the Mediterranean Sea) and 31 vertical levels in the ocean (10 levels in the top 100 m) (R182 × 149L31) (Madec 2008); the LIM2 sea-ice model (Louvain-la-Neuve Sea Ice Model) (Fichefet and Morales-Maqueda 1997); and the ORCHIDEE (Organizing Carbon and Hydrology In Dynamic Ecosystems) land surface model (Krinner et al. 2005). IPSL-CM5A-LR couples the ocean, atmosphere, land and sea ice through the OASIS3 coupler (Valcke 2006). For the atmosphere-only simulations, LMDZ5A is used together with the ORCHIDEE land surface model, with fixed vegetation (see Contoux et al. 2012).

Experimental designs for the mid-Pliocene and PI simulations
Following PlioMIP guidelines, two types of simulations for mid-Pliocene warm period are conducted using the atmosphere-only climate model (LMDZ5A) and the fully coupled climate model (IPSL-CM5A-LR). The boundary conditions derived from PRISM3D are set for the mid-Pliocene simulation, as proposed by PlioMIP experimental designs (Contoux et al. 2012). The modern coastline is employed in these simulations to avoid changing the land-sea mask in the ocean model. Differences between the mid-Pliocene and modern topographies are added as anomalies to the IPSL-CM5A-LR model topography (Edwards et al. 1992;Sohl et al. 2009), but include few changes apart from the ice sheets. Imposed ice-sheet and vegetation reconstructions are derived from PRISM3D (Hill et al. 2007;Salzmann et al. 2008), and include a Greenland Ice Sheet reduced by 50 % in volume, and an Antarctic Ice Sheet reduced by 33 %. Prescribed vegetation includes a reduction of most present-day desert, and a northward shift of temperate forest biomes. The AGCM simulation for the mid-Pliocene is forced by prescribed monthly SST derived from the PRISM reconstruction. The orbital configuration, solar constant, greenhouse gases and aerosols are specified as similar to in the PI experiment except that the atmospheric CO 2 concentration is prescribed at 405 ppm. More details on the experimental design, the implementation of PlioMIP boundary conditions, as well as the general features of the mid-Pliocene climate simulated by LMDZ5A and IPSL-CM5A can be found in Contoux et al. (2012).
The PI experiments are performed by both the AGCM and the CGCM to enable direct comparisons with the corresponding mid-Pliocene simulations. Greenhouse gas concentrations, the solar constant and orbital parameters are set to the PI values recommended by the Coupled Model Intercomparison Project Phase 3 and Phase 5 (CMIP3 and CMIP5). The solar constant is 1365 W m −2 . The CO 2 and CH 4 concentrations are fixed at 280 ppm and 760 ppb, respectively.
The IPSL-CM5A-LR model and its atmospheric component (LMDZ5A) performance as well as other PlioMIP models in estimating the global annual mean SAT for the mid-Pliocene, i.e. at around 2 °C higher than present (Contoux et al. 2012;Haywood et al. 2013). Here, we focus mainly on the similarities and differences between the EASM precipitation changes in the context of the mid-Pliocene warm period simulated by IPSL-CM5A-LR and LMDZ5A.

Moisture budget analysis
Moisture budget analysis is widely used to understand the global and regional precipitation changes (Yoon and Chen 2005;Chou and Lan 2012;Hsu et al. 2012). Here, to estimate the processes that are important for the EASM precipitation, the moisture budget equation is analyzed (Seo et al. 2013). The vertically integrated moisture budget equation can be expressed as where q is specific humidity, V is three-dimensional wind, and 〈 〉 represents a vertical integration from 1000 to 100 hPa. E is surface evaporation and P is precipitation. Res is the residual term, including transient eddies and nonlinear effects, and is relatively small compared with other terms.
The term on the left-hand side in Eq. (1), moisture flux convergence (〈∇ · q � V〉), is further divided into three terms: vertical moisture advection −�ω∂ P q� (ω is vertical velocity in the pressure coordinate); horizontal moisture advection −〈 � V H · ∇q〉 ( V H is horizontal wind); and moisture divergence −〈q · ∇ · � V H 〉. Thus, Eq. (1) can be expressed as Changes of vertical moisture advection can be further divided into thermodynamic and dynamic components, to assess the relative contribution of specific humidity and atmospheric circulation to regional precipitation variability (Chou and Lan 2012). Thus, where − is climatology and 〈 〉′ is the departure of change from the climatology, and − ω∂ p q ′ and − ω ′ ∂ pq denote the thermodynamic and dynamic contribution to the changes of vertical moisture advection, respectively. The thermodynamic contribution is directly affected by specific humidity changes, which is closely tied to temperature changes roughly following Clausius-Clapeyron (C-C) equation (Chou et al. 2013). The dynamic contribution is directly linked to changes in large-scale vertical motion (Chou et al. 2013).

Moist static energy diagnosis
Following Chen and Bordoni (2014a), MSE is selected to quantify the relative roles of temperature, moisture, and radiative processes in sustaining the horizontal distribution of the mid-Pliocene EASM precipitation. The highly simplified MSE balance in the atmospheric column is expressed as where MSE is the atmospheric moist static energy, E is the atmospheric moist enthalpy, and F net is the net energy flux into the atmosphere. The atmospheric moist static energy and moist enthalpy are expressed as MSE = C p T + gz + L v q and E = C p T + L v q, where C p is specific heat at constant pressure, T is temperature, z is geopotential height, g is gravitational acceleration, L v is latent heat of vaporization, and q is specific humidity. The MSE stratification in the atmosphere is negative, hence regions of positive (negative) vertical MSE advection correspond to ascending (descending) vertical motion. The time-mean advection of the atmospheric energy 〈V · ∇E〉 can be written as where (·)′ and (·)* denote deviations from the time mean (·) and zonal mean [·], respectively. The first term on the right-hand side of Eq. (5) is the zonal-mean energy advection by the zonal mean flow; the second term is the advection of the stationary eddy energy by the zonal-mean flow; the third term is the advection of the zonal-mean energy by the stationary eddy velocity; the fourth term is the advection of the stationary eddy energy by the stationary eddy velocity; the fifth term is the advection of the transient eddy energy by the transient eddies. The first term is very small compared with the other terms and can be neglected. A more detailed description of the MSE used to analyze the EASM can be found in Chen and Bordoni (2014a).

Simulated spatial distribution of the EASM precipitation changes in the mid-Pliocene
The EASM precipitation in the mid-Pliocene simulated by the coupled model and the atmosphere-only model share similarities and differences ( Fig. 1). Figure 1a, b show the maps of the simulated differences in precipitation between the mid-Pliocene and PI experiments. The red boxes in Fig. 1a, b denote the selected EASM domain for the present study. The EASM precipitation changes in the mid-Pliocene derived from the AGCM simulation are nonuniform, with a decrease in precipitation along the southeastern coast of China and increased precipitation on the northwestern flank ( Fig. 1a) compared to PI. In contrast, a general increase in summer monsoon precipitation over East Asia is seen in the coupled simulation for the mid-Pliocene warm period compared to PI, associated with a maximal increase of precipitation confined to southern China ( Fig. 1b). Despite the area-averaged mean summer monsoon precipitation over East Asia increasing in both simulations when compared with the corresponding PI experiments (Fig. 3b), important differences exist between the AGCM and CGCM simulations for the mid-Pliocene. Compared with the AGCM simulation, the increased mid-Pliocene EASM precipitation in the CGCM simulation is confined to the south of 25°N, while the decreased precipitation anomalies are northward of 25°N (Fig. 1c).
What caused the similarities and large differences in the horizontal distributions of EASM precipitation in the mid-Pliocene between the two simulations? There are differences in SAT and SST between the mid-Pliocene and PI (Fig. 2), and between the CGCM and AGCM for the mid-Pliocene (Fig. 2c, d). What is the impact of these differences in SST and SAT on the precipitation patterns of the EASM? What are the dominant mechanisms through which the SST and SAT differences impact upon the EASM precipitation? To answer these questions, we analyze the moisture budget and diagnose the MSE in the following section.

Moisture budget analysis of the area-averaged EASM precipitation changes in the mid-Pliocene
Moisture budget analysis is carried out to examine the simulated area-averaged EASM precipitation in the mid-Pliocene warm period. Climatologically, local surface evaporation and vertical moisture advection rather than horizontal moisture advection are the main contributors to the precipitation over East Asia (20°-45°N, 105°-135°E) (Fig. 3a). Furthermore, moisture divergence contributes a smaller amount to sustain the mean state of the EASM precipitation when compared with local surface evaporation and vertical moisture advection, while the contribution of the residual term is negative (Fig. 3a). Moisture budget analysis additionally reveals the dominant factor determining the simulated similarities of the EASM precipitation changes in the mid-Pliocene between the two simulations. Figure 3b shows the AGCM-and CGCM-simulated moisture budget differences between the mid-Pliocene and PI. The simulated enhancement of mid-Pliocene precipitation over the EASM domain in both simulations is mainly attributed to increased vertical moisture advection (−〈ω∂ p q〉′). The vertical moisture advection (−〈ω∂ p q〉′) is further divided into thermodynamic and dynamic components to assess the relative contribution of specific humidity (q) and atmospheric circulation (ω) to the enhanced mid-Pliocene EASM precipitation. The results show that the simulated increase of vertical moisture advection (−〈ω∂ p q〉′) over East Asia results from an increase in the thermodynamic component (−〈ω∂ p q ′ 〉), thereby leading to enhanced EASM precipitation in the mid-Pliocene AGCM and CGCM simulations (Fig. 3b). The simulated increase in the thermodynamic component's contribution to the enhanced EASM precipitation in the mid-Pliocene essentially reflects the increasing moisture in the atmospheric column in response to surface warming. For instance, a 1.2 K (1.8 K) increase in summer (June-July-August; JJA) mean SAT over East Asia (Fig. 2) derived from the AGCM (CGCM) simulation for the mid-Pliocene produces about a 6.3 % (7.8 %) increase in the atmospheric moisture content.
Moisture budget analysis is also used to reveal what causes the large aforementioned differences of the mid-Pliocene EASM precipitation between the AGCM and CGCM simulations. The mean JJA EASM precipitation in the mid-Pliocene is larger in the coupled simulation than in the atmosphere-only simulation (Fig. 3a,c). There is also a larger increase in the mid-Pliocene EASM precipitation relative to the PI experiment in the coupled model (IPSL-CM5A-LR) compared to that in its atmospheric component alone (LMDZ5A) (Fig. 3b). The simulated larger increase in the mid-Pliocene EASM precipitation in the coupled model results mainly from the larger vertical moisture advection (−〈ω∂ p q〉′), which is due to the larger increase in the dynamic component (−〈ω ′ ∂ pq 〉), itself closely related to the enhanced vertical motion in the coupled simulation compared to the atmosphere-only simulation (Fig. 3c).
Moisture budget analysis provides a useful way to identify the driver of the enhanced summer monsoon precipitation over East Asia during the mid-Pliocene in the atmosphere-only and coupled simulations with IPSL-CM5A. In essence, we are able to explain that it is the increase in SAT over East Asia that drives the increased EASM precipitation in the mid-Pliocene simulations compared to PI: The increase in SAT produces an increase in atmospheric a b c Fig. 1 The horizontal distributions of JJA precipitation (units: mm/ day) differences between the mid-Pliocene (MP) and PI experiments in the a AGCM and b CGCM. c The simulated JJA MP precipitation difference between the CGCM and AGCM. The red box denotes the EASM domain (20°-45°N, 105°-135°E) moisture content, revealed by the increase of the thermodynamic component of the vertical moisture advection. Additionally, the dynamic component contributes more to the vertical moisture advection in the coupled simulation than in the atmosphere-only simulation (Fig. 3c). That is why there is a larger increase in the mid-Pliocene precipitation in the coupled model simulation than in the atmosphereonly simulation. This essentially reveals an enhanced role of vertical motion in the mid-Pliocene EASM precipitation increase in the coupled model.

The dominant factors determining the EASM precipitation changes in the mid-Pliocene
The cause of the enhanced area-averaged EASM precipitation in both the mid-Pliocene simulations compared to the corresponding PI simulations, as well as the simulated differences of the mid-Pliocene EASM precipitation between the two simulations, has been revealed in the previous section by performing a moisture budget analysis in the EASM domain. But what is the mechanism that sustains the spatial distribution of this enhanced mid-Pliocene EASM precipitation in both simulations, and the large differences in the mid-Pliocene EASM precipitation between the two simulations? MSE has recently been proposed as an approach to reveal the physical mechanism underlying the observed maintenance of EASM precipitation (Chen and Bordoni 2014a), and has also been used to assess model performance in simulating EASM precipitation (Chen and Bordoni 2014b).

Explaining the EASM precipitation changes in the mid-Pliocene using MSE budget analysis
Following Chen and Bordoni (2014a), MSE budget analysis is used to reveal the possible factor determining the EASM precipitation changes in the mid-Pliocene derived from the AGCM and CGCM simulations. We first examine the JJA-averaged MSE budget integrated from the surface to 100 hPa derived from the AGCM simulation for the PI and mid-Pliocene climate. The spatial distributions of vertical MSE advection (〈ω∂ p MSE〉) derived from the AGCM simulation for the PI and mid-Pliocene experiments resemble the corresponding large-scale summertime precipitation pattern over East Asia (Fig. 4a, d, contours). Given that vertical MSE advection (〈ω∂ p MSE〉) is computed as the vertically integrated product of vertical velocity in the p-coordinates ω and ∂ p MSE, associated with the column-integrated negative ∂ p MSE in the troposphere (i.e. 〈∂ p MSE〉 < 0), the regions where the vertical MSE advection is positive (〈ω∂ p MSE〉 > 0) thus correspond to the ascending motion over East Asia, and vice versa (Chen and Bordoni 2014a, b). Therefore, the vertical MSE advection (〈ω∂ p MSE〉) can be considered as a good gauge to identify the vertical motion over East Asia. These results also provide an alternative way to understand the mechanism sustaining the rain belt position over East Asia by revealing the dominant factor determining the horizontal distribution of column-integrated vertical MSE advection (〈ω∂ p MSE〉). As shown in Fig. 4, the vertically integrated JJA-averaged vertical MSE advection (〈ω∂ p MSE〉) derived from the PI and mid-Pliocene experiments simulated by the AGCM are further partitioned into two terms: net energy flux into the atmospheric column F net . and vertically integrated horizontal moist enthalpy advection (−〈V · ∇E〉). The spatial distribution of the horizontal moist enthalpy advection (−〈V · ∇E〉) over East Asia derived from the PI and mid-Pliocene experiments is similar to the spatial distribution of the vertical MSE advection (〈ω∂ p MSE〉). (Figure 4c, f), while the net energy flux into the atmospheric column (F net ) over East Asia exhibits a horizontal distribution that is distinctly different to that of the vertical MSE advection (〈ω∂ p MSE〉), with positive values over the China mainland and minimal net energy flux into the atmospheric column (F net ) extending from the Bohai Sea to Japan sea (Fig. 4b, e). Thus, the spatial patterns of JJA-averaged vertical MSE advection (〈ω∂ p MSE〉) over East Asia derived from the PI and mid-Pliocene experiments simulated by the AGCM are largely determined by moist enthalpy advection (−〈V · ∇E〉) rather than the energy flux in the atmospheric column (F net ).
We also examine the AGCM-simulated differences of the following three terms between the mid-Pliocene and PI experiments: 〈ω∂ p MSE〉, F net , −〈V · ∇E〉. There are large similarities in the AGCM-simulated spatial patterns between the vertical MSE advection (〈ω∂ p MSE〉) and EASM precipitation changes in the mid-Pliocene (Fig. 4g), indicating that vertical MSE advection (〈ω∂ p MSE〉) is a useful gauge to measure the horizontal distributions of the AGCM-simulated EASM precipitation changes in the mid-Pliocene (Figs. 1a, 4g). The simulated vertical MSE advection (〈ω∂ p MSE〉) decreases along the southeastern coast of China, and increased vertical MSE advection (〈ω∂ p MSE〉) on the northwestern flank largely results from the horizontal moist enthalpy advection (−〈V · ∇E〉) changes in the mid-Pliocene compared to PI (Fig. 4i). In contrast, changes in the AGCM-simulated net energy flux into the atmospheric column (F net ) in the mid-Pliocene are much smaller compared to those of horizontal moist enthalpy advection (−〈V · ∇E〉), and thus its contribution to the changes in vertical MSE advection (〈ω∂ p MSE〉) can be neglected.
Following the same analysis method, we also examine the MSE budget-related variables of vertical MSE advection (〈ω∂ p MSE〉), net energy flux into the atmospheric column (F net ), and horizontal moist enthalpy advection (−〈V · ∇E〉), as derived from the CGCM simulations. The spatial patterns of JJA-averaged vertical MSE advection (〈ω∂ p MSE〉) derived from the CGCM simulations for the PI and mid-Pliocene experiments resemble the horizontal distributions of the corresponding JJA-averaged EASM precipitation derived from the PI and mid-Pliocene experiments, respectively (Fig. 5a, d, contours). This result suggests that vertical MSE advection (〈ω∂ p MSE〉) is also a good indicator of the spatial pattern of the CGCM-simulated summer monsoon precipitation over East Asia. The net energy flux into the atmospheric column (F net ) and horizontal moist enthalpy advection (−〈V · ∇E〉) derived from the CGCM simulation are also shown in order to assess the relative importance of these terms in sustaining the spatial distributions of JJA-averaged vertical MSE advection (〈ω∂ p MSE〉) and its changes in the mid-Pliocene (Fig. 5). The CGCM-simulated spatial patterns of vertical MSE advection (〈ω∂ p MSE〉) are also largely determined by the horizontal distributions of horizontal moist enthalpy advection (−〈V · ∇E〉) with respect to the mean state of vertical MSE advection (〈ω∂ p MSE〉) and its changes in the mid-Pliocene (Fig. 5).
In brief, the spatial patterns of the summer monsoon precipitation over East Asia with respect to the mean state and its changes in the mid-Pliocene, as revealed by the coupled model and its atmospheric component alone, are explained well by the spatial patterns of vertical MSE advection (〈ω∂ p MSE〉), which are largely determined by the horizontal distributions of moist enthalpy advection (−〈V · ∇E〉) in both simulations. However, the mechanisms controlling the spatial pattern of moist enthalpy advection (−〈V · ∇E〉) remain to be further analyzed, so as to provide a reasonable explanation for the enhanced EASM precipitation in the mid-Pliocene. This is the focus of the following section.

Identifying the dominant factors determining the EASM precipitation changes in the mid-Pliocene
As mentioned above, column-integrated horizontal moist enthalpy advection (−〈V · ∇E〉) is an effective metric to describe the spatial distributions of AGCM-and CGCMsimulated summer monsoon precipitation over East Asia (Figs. 4,5). According to Eq. (5), this term (−〈V · ∇E〉) can be decomposed into five terms. Here, we have only focused on the second term (−〈[V ] · ∇E * 〉), which is the advection of the stationary eddy energy by the zonalmean flow, and the third term (−〈V * · [∇E]〉), which is the advection of the zonal-mean energy by the stationary eddy velocity, while the other terms on the right-hand side in Eq. (5) are neglected, as suggested by Chen and Bodorni (2014a, b). Given that atmospheric moist enthalpy (E) can be expressed as a sum of dry enthalpy C p T and latent heating (L v q), each of the aforementioned E-relevant advection terms (−〈V · ∇E〉, −〈[V ] · ∇E * 〉, −〈V * · [∇E]〉) can be further represented by the sum of the corresponding temperature-relevant advection (Figs. 6, 7, middle panels) and specific humidity-relevant advection (Figs. 6, 7, right-hand panels). As shown in Fig. 6a, d and g, the JJA-averaged variables of −〈V · ∇E〉, −〈[V ] · ∇E * 〉 and −〈V * · [∇E]〉 derived from the AGCM simulation for the PI (contours) and mid-Pliocene (shading) experiments, associated with separate contributions from dry enthalpy (Fig. 6b, e, h) and latent heating (Fig. 6c, f, (Fig. 6d) and the advection of the zonalmean energy by the stationary eddy velocity (−〈V * · [∇E]〉) (Fig. 6g). The moist enthalpy advections (−〈V · ∇E〉) (Fig. 6a) derived from the PI and mid-Pliocene experiments are both largely determined by the total dry enthalpy advection (−〈−CpV · ∇T 〉, Fig. 6b), while the AGCMsimulated total latent heating advection (−〈−LvV · ∇q〉, Fig. 6c) derived from the PI and mid-Pliocene experiments both contribute less to the moist enthalpy advections (−〈V · ∇E〉) and are thus neglected in the current study. The total dry enthalpy advections derived from the AGCM simulation for the PI and mid-Pliocene experiments shown in Fig. 6b are largely attributed to the advection of the stationary eddy dry enthalpy by zonal-mean flow (−〈[CpV ] · ∇T * 〉, Fig. 6e) and the advection of zonal-mean dry enthalpy by the stationary eddy velocity (−〈V * · [Cp∇T ]〉, Fig. 6h).
Following the same analysis method, we also examine the moist enthalpy advection (−〈V · ∇E〉) and its associated components, the total dry enthalpy advection (−〈CpV · ∇T 〉), and the total latent heating advection (−〈−LvV · ∇q〉) derived from the coupled simulations for the PI and mid-Pliocene experiments. The moist enthalpy advections (−〈V · ∇E〉, Fig. 7a) in the coupled simulations for the PI and mid-Pliocene experiments are also both largely determined by the total dry enthalpy advection (−〈−CpV · ∇T 〉, Fig. 7b). The advection of the stationary eddy dry enthalpy by zonal-mean flow (−〈[CpV ] · ∇T * 〉, Fig. 7e) and the advection of the zonal-mean dry enthalpy by the stationary eddy velocity (−〈V * · [Cp∇T ]〉, Fig. 7h) are both dominant factors determining the spatial patterns of the total dry enthalpy advection (−〈−CpV · ∇T 〉, Fig. 7b) derived from the coupled simulations for the PI and mid-Pliocene experiments.
As discussed above, the total dry enthalpy advection (−〈−CpV · ∇T 〉) associated with its components (−〈[CpV ] · ∇T * 〉 and −〈V * · [Cp∇T ]〉) in both simulations are the dominant factors determining the spatial patterns of the JJA-averaged moist enthalpy advection (−〈V · ∇E〉) over East Asia. Whether these aforementioned dominant factors used to identify the JJA-averaged moist enthalpy advection (−〈V · ∇E〉) enable us to explain the moist enthalpy advection (−〈V · ∇E〉) changes in the mid-Pliocene in the EASM domain derived from the AGCM and CGCM simulations is further discussed in the following paragraph. As shown in Fig. 8a, the AGCM-simulated moist enthalpy advection (−〈V · ∇E〉) changes in the mid-Pliocene largely result from the total dry enthalpy advection (−〈−CpV · ∇T 〉, Fig. 8b) changes in the mid-Pliocene simulation, while the AGCM-simulated latent heat advection (−〈−LvV · ∇q〉, Fig. 8c) changes in the mid-Pliocene contribute less to the moist enthalpy advection (−〈V · ∇E〉) changes in mid-Pliocene. The two components of the total dry enthalpy advection (−〈−CpV · ∇T 〉), the advection of the stationary eddy dry enthalpy by zonal-mean flow (−〈[CpV ] · ∇T * 〉), and the advection of zonal-mean dry enthalpy by the stationary eddy velocity (−〈V * · [Cp∇T ]〉) both contribute to the AGCM-simulated non-uniform distribution of the total dry enthalpy advection (−〈−CpV · ∇T 〉) changes in the mid-Pliocene over East Asia (Fig. 8e-h). This means that the negative anomalies of the total dry enthalpy advection (−〈−CpV · ∇T 〉) in the mid-Pliocene along the southeastern coast of China can be largely attributed to the decreased advection of the stationary eddy dry enthalpy by zonal mean flow (−〈[CpV ] · ∇T * 〉) in that region (Fig. 8e). Meanwhile, the simulated total dry enthalpy advection (−〈−CpV · ∇T 〉) increases on the northwestern flank of the aforementioned negative anomalies as a result of the combined contribution of the increased advection of the stationary eddy dry enthalpy by zonal-mean flow (−〈[CpV ] · ∇T * 〉) extending from Northeast China to the Korean peninsula (Fig. 8e), and the advection of the zonal-mean dry enthalpy by the stationary eddy velocity (−〈V * · [Cp∇T ]〉) from Southwest China to the Shandong peninsula (Fig. 8h).
The CGCM-simulated moist enthalpy advection (−〈V · ∇E〉) changes associated with the two components −〈[CpV ] · ∇T * 〉 and −〈V * · [Cp∇T ]〉 the mid-Pliocene are also examined to identify the dominant factors determining the spatial pattern of moist enthalpy advection (−〈V · ∇E〉) changes in the mid-Pliocene over East Asia (Fig. 9a, d,  g). The enhanced moist enthalpy advection (−〈V · ∇E〉, Fig. 9a) in the CGCM mid-Pliocene simulation is also largely dominated by the general increase in the total dry enthalpy advection (−〈−CpV · ∇T 〉, Fig. 9b), which is also attributed to the combined increases of advection of the stationary eddy dry enthalpy by zonal-mean flow As discussed above, the moist enthalpy advection largely attributed to the dry enthalpy advection is a good gauge to describe the EASM precipitation with respect to its mean state and changes in the mid-Pliocene. As such, the JJA-averaged dry enthalpy advections in both simulations have been examined first to confirm the two dominant components of the dry enthalpy advection (advection of the stationary eddy dry enthalpy by zonal-mean flow −〈[CpV ] · ∇T * 〉 and advection of the zonal-mean dry enthalpy by the stationary eddy velocity −〈V * · [Cp∇T ]〉) sustaining the spatial pattern of JJA-averaged dry moisture advection. These two components are also considered to be the two dominant factors determining the horizontal distribution of dry enthalpy advection changes in the mid-Pliocene. This means that the two factors of (1) advection of the stationary eddy dry enthalpy by zonal-mean flow (−〈[CpV ] · ∇T * 〉) and (2)  A recent modeling study confirmed that the spatial pattern of JJA-averaged advection of the stationary eddy dry enthalpy by zonal-mean zonal flow (−〈[CpU] · ∂T * x 〉) over East Asia indeed reflects the thermal effect of the Tibetan Plateau on the EASM circulation through enhancing the land-sea thermal contrast and transporting warm and moist air advection to downstream of the Tibetan Plateau into East Asia (Chen and Bodorni 2014a). In the current study, we use this term's (−〈[CpU] · ∂T * x 〉) changes in the mid-Pliocene compared to the PI experiment to explain the enhanced mid-Pliocene EASM circulation and associated enhanced moisture transport from the tropical ocean to the EASM region (Fig. 10e, f). As shown in Fig. 2, the simulated zonal thermal contrasts are enhanced over East Asia in the mid-Pliocene derived from the AGCM and CGCM compared to the corresponding PI experiments, thereby strengthening the EASM circulation and moisture transport into the EASM region in both simulations (Fig. 10e,  f), resulting in the enhanced mid-Pliocene EASM precipitation in both simulations (Fig. 3b).
Column-integrated meridional stationary eddy velocity (〈v * 〉) as a component of the advection of the zonal-mean dry enthalpy by the meridional stationary eddy velocity ((−〈−CpV * · [∂T y ]〉) plays the fundamental role in determining the spatial pattern of the advection of the zonalmean dry enthalpy by the meridional stationary eddy velocity (−〈−CpV * · [∂T y ]〉) with respect to its mean state and changes in the mid-Pliocene (not shown). Nevertheless, the way that the column-integrated meridional stationary eddy velocity (〈v * 〉) effects the summer monsoon precipitation over East Asia remains to be further clarified. Zonally averaged meridional stationary eddy velocity (v * ) over East Asia can capture well the characteristics of the JJA-averaged EASM circulation derived from the AGCM and CGCM simulations, as evidenced by the prevailing southerlies (northerlies) in the lower (upper) troposphere associated with the lower-level (upper-level) convergence (divergence) induced by the non-uniform distribution of meridional stationary eddy velocity (v * ) over the EASM domain ( Fig. 11a, b, d, e). The simulated increase of meridional stationary eddy velocity (v * ) and its convergence ∂v * ∂y over the EASM domain are evident in the mid-Pliocene compared to the PI experiments in both simulations. These simulated increases are beneficial to the enhanced moisture transport into the EASM domain and increased local moisture convergence in the mid-Pliocene in both simulations, respectively (Fig. 11c, f). Therefore, these simulated increases of meridional stationary eddy velocity (v * ) and its convergence ∂v * ∂y over the EASM domain during the mid-Pliocene in both simulations contribute to the enhanced EASM precipitation in the mid-Pliocene simulated by both the AGCM and the CGCM experiments (Fig. 3a, b).

Explaining the differences in the mid-Pliocene EASM precipitation between CGCM and AGCM simulations
Having identified the possible factors-advection of the stationary eddy dry enthalpy by zonal-mean zonal flow (−〈[CpU] · ∂T * x 〉) and column-integrated meridional stationary eddy velocity (〈v * 〉)-determining the spatial patterns of JJA-averaged EASM precipitation and its changes in the mid-Pliocene relative to the corresponding PI experiments derived from the AGCM and CGCM simulations (Fig. 1a, b), we now turn our attention to the key issue of whether or not these factors can be used to explain the differences in EASM precipitation between the CGCM and AGCM simulations (Fig. 1c). The differences in the mid-Pliocene moist enthalpy advection between the CGCM and AGCM simulations in terms of spatial distribution resemble the non-uniform distribution of the EASM precipitation differences between the CGCM and AGCM simulations (Figs. 1c, 12a). Therefore, moist enthalpy advection can also be considered as the dominant factor determining the simulated mid-Pliocene EASM precipitation differences between the CGCM and AGCM. Identifying the simulated differences of the dominant factors determining the moist enthalpy advection thus provides an alternative way to understand the mechanism controlling the simulated EASM precipitation difference between the two simulations.
The enhanced meridional thermal contrast derived from CGCM compared with the AGCM (Fig. 2c) leads to a weakened EASM circulation in the CGCM simulation compared to the AGCM simulation, associated with a b c d e f g h i Fig. 9 The CGCM-simulated differences of the corresponding variables shown in Fig. 7 between the mid-Pliocene and PI experiments (units: decreased moisture transport into East Asia (Fig. 12b), which results in decreased EASM precipitation over the northwest of South China (Fig. 1c). These results provide a reasonable explanation as to why the thermodynamic component of vertical moisture advection contributes negatively to the enhanced mean EASM precipitation in the mid-Pliocene CGCM simulation compared to the AGCM simulation (Fig. 3c). The weakened EASM in the mid-Pliocene CGCM simulation compared to the AGCM simulation is closely associated with the northeasterly anomalies on the northwestern flank of the anomalous cyclonic circulation dominant over South China and the South China Sea (Fig. 12b). This anomalous cyclone is closely linked to the warmer SST anomalies in the CGCM simulation than the AGCM simulation (Fig. 2d). Note that the weakened mid-Pliocene EASM circulation derived from CGCM simulation compared with the AGCM simulation features non-uniform spatial distributions, as evidenced by the simulated differences of meridional stationary eddy velocity convergence between the two simulations (Fig. 12c). This non-uniform weakening of the mid-Pliocene EASM circulation in the CGCM simulation compared to the AGCM simulation is highly consistent with the non-uniform distribution of the mid-Pliocene EASM precipitation differences between the CGCM and AGCM simulations (Fig. 1c).
Relative to the AGCM, the CGCM-simulated increase of stationary eddy convergence over South China associated with increased vertical motion (Fig. 12f) thereby leads to the enhanced mid-Pliocene EASM precipitation over South China (Fig. 1c).
In brief, the enhanced meridional thermal contrast over the EASM domain and the increased stationary eddy convergence over South China both contribute to sustain the non-uniform distribution of the mid-Pliocene EASM precipitation difference between the CGCM and AGCM simulations. The former makes a positive contribution to the weakened precipitation in the CGCM simulation on the northwestern flank of South China, while the latter can explain well the enhanced precipitation over South China in the CGCM simulation relative to the AGCM simulation. Thus, the enhancement of area-averaged EASM precipitation in the CGCM relative to the AGCM is attributed to the increased dynamic component of vertical advection rather than the thermodynamic component (Fig. 3c).  (Fig. 13a). These models used to validate the robustness of enhanced mid-Pliocene EASM precipitation based on IPSL-CM5A have been well documented in previous study . The enhanced zonal thermal contrast has been identified as a key mechanism controlling the enhancement of mid-Pliocene EASM precipitation by performing a comparative analysis of mid-Pliocene warm climate between IPSL-CM5A-LR and its stand-alone AGCM (LMDZ5) simulations. This enhanced zonal thermal contrast associated with column-integrated v * increase is also evident in  Fig. 13 The PlioMIP models simulated a mid-Pliocene EASM precipitation changes relative to PI experiment (by percentage)versus surface warming (units: °C) over EASM domain; b zonal thermal contrast changes (units: °C) in mid-Pliocene; c simulated columnintegrated v * (units: m s −1 ) changes in mid-Pliocene by PlioMIP models and the corresponding stand-alone AGCMs; d same as Fig. 3; except for the simulated differences of mid-Pliocene EASM precipi-tation between CGCM and AGCM as represented by other PlioMIP models; e simulated differences of mid-Pliocene ∂v * ∂y at 1000 hPa (units: 10 −6 s −1 ) averaged over (105°-135°E) between CGCM and AGCM derived from other PlioMIP models; f the SST bias (units: °C) between multi-CGCMs ensemble mean and PRISM3 SST imposed in AGCMs (Fig. 13b, c), which are highly consistent with enhanced mid-Pliocene EASM precipitation revealed by these models. Note that a decrease of mid-Pliocene EASM precipitation revealed by other six models shown in Fig. 13a is also closely related to a weaker zonal thermal contrast in these models (Fig. 13b, thick purple: ensemble mean).
Multi-model comparisons are also carried out to investigate the simulated differences of mid-Pliocene EASM precipitation between CGCMs and corresponding atmospheric components. The simulated EASM precipitation is larger in these coupled models than their atmospheric components (Fig. 13d), which is consistent with the results from IPSL-CM5A model. The multi-model simulated JJA-averaged convergence of meridional stationary eddy velocity ( ∂v * ∂y ) at lower tropopause can well capture the lower-level convergence over the EASM domain (not shown), thus ∂v * ∂y has been considered a good metrics to gauge the mid-Pliocene EASM precipitation differences between the two simulations (Fig. 13e). The larger mid-Pliocene EASM precipitation in these CGCM simulations than their AGCM simulations is closely associated with CGCM-simulated increase in the ∂v * ∂y over EASM domain compared to that of AGCM, despite that there are large inter-model spread of this convergence center over EASM domain (Fig. 13e). This increase of convergence of meridional stationary eddy velocity ( ∂v * ∂y ) over EASM domain in IPSL-CM5A simulation is driven by warmer SST anomalies over South China (Fig. 2d). There are also robust SST warmer biases over south China and tropical ocean between CGCM and PRISM3 SST imposed in stand-alone AGCM (Fig. 13f), which is consistent with the large overestimation of SST at low latitude in coupled PlioMIP1 models than PRISM3 SST in previous studies Salzmann et al. 2013). Thus the larger simulated mid-Pliocene EASM precipitation in coupled models than their atmospheric components is essentially driven by this warmer SST biases between CGCM and AGCM over South China through inducing larger increase in ∂v * ∂y over EASM domain in coupled models than corresponding AGCMs. A data/model comparison on this area should be carried out in near future to investigate whether the SST in couple models or PRISM3 SST are realistic.
In brief, multi-model comparisons demonstrate that the relationships to explain the differences between CGCM and AGCM in the case of IPSL-CM5A have been confirmed when using the model database of PlioMIP. This result pinpoints the essential impact to have reliable SST reconstruction to infer the EASM evolution during mid-Pliocene.

Summary and conclusions
The drivers and mechanisms of enhanced summer monsoon precipitation over East Asia during the mid-Pliocene in atmosphere-only and coupled simulations with IPSL-CM5A, as well as the simulated differences between the two simulations, are revealed in the present paper by performing a comparative analysis. This analysis demonstrates that the increase in area-averaged EASM summer precipitation in the mid-Pliocene simulations compared to the PI simulations is mostly driven by the increase in SAT over East Asia in both the mid-Pliocene simulations. The surface warming increases atmospheric moisture, thereby leading to an increased thermodynamic component of the vertical moisture advection. MSE diagnosis provides an effective way to identify the combined effect of enhanced thermal contrast and increased column-integrated meridional stationary eddy velocity (〈v * 〉) as the principal mechanism controlling the large similarities and differences of the spatial patterns of EASM precipitation changes in the mid-Pliocene between the two simulations. A simple schematic illustration is provided in Fig. 14 to summarize these drivers and mechanisms using moisture and MSE budget analysis, respectively. The study and its main conclusions can be summarized as follows: By performing a moisture budget analysis in the EASM domain, we first examine the causes of the similarities and differences in area-averaged mid-Pliocene EASM precipitation changes compared to the PI experiments in both the simulations, as well as the differences in the mid-Pliocene EASM precipitation between the two simulations. The simulated mid-Pliocene precipitation over the EASM domain is enhanced in both simulations due to the increased vertical moisture advection in both simulations. The increased contribution of the thermodynamic component to the enhanced vertical moisture advection essentially reflects the response of atmospheric moisture to the mid-Pliocene global surface warming in both simulations. This result indicates that the enhanced mid-Pliocene EASM could be used as a test bed to further validate the projection of an intensified EASM in the future (Bao 2012). The dynamic component rather than thermodynamic component contributes to the increased vertical moisture advection in the CGCM simulation than the AGCM simulation, thereby leading to larger mid-Pliocene EASM precipitation in the CGCM simulation compared to the AGCM simulation.
The MSE budget analysis provides an effective way to identify the dominant factors determining the spatial distributions of JJA-averaged EASM precipitation and its changes in the mid-Pliocene derived from both simulations. The spatial patterns of the EASM precipitation with respect to the mean state and its changes in the mid-Pliocene derived from both simulations are explained well by the dry enthalpy advection. The advection of the stationary eddy dry enthalpy by zonal-mean zonal flow (−〈[CpU] · ∂T * x 〉) and the advection of the zonal-mean dry enthalpy by the meridional stationary eddy velocity (−〈−CpV * · [∂T y ]〉) are identified as the two dominant components of the dry enthalpy advection, and thus they are selected to explain the EASM precipitation changes in the mid-Pliocene derived from both simulations. The simulated zonal thermal contrasts are enhanced over East Asia in the mid-Pliocene derived from both simulations, as evidenced by the advection of the stationary eddy dry enthalpy by zonal-mean zonal flow (−〈[CpU] · ∂T * x 〉), thereby strengthening the EASM circulation and moisture transport into the EASM domain, and thus contributing to the enhancement of the mid-Pliocene EASM precipitation derived from the AGCM and CGCM simulations compared with the corresponding PI experiments. Column-integrated meridional stationary eddy velocity (〈v * 〉) as a dominant component of the advection of the zonal-mean dry enthalpy by the meridional stationary eddy velocity (−〈−CpV * · [∂T y ]〉) can reasonably capture the characteristics of the JJA-averaged EASM circulation derived from the AGCM and CGCM simulations, and thus the enhanced column-integrated meridional stationary eddy velocity (〈v * 〉) associated with the increase of lower-level convergence and upper-level divergence over the EASM domain also contribute to the enhancement of EASM precipitation in the mid-Pliocene derived from both simulations.
Simulated differences of the mid-Pliocene EASM precipitation between the CGCM and AGCM are non-uniform, despite the simulated EASM precipitation in the mid-Pliocene compared with the PI experiment being enhanced in both the Fig. 14 Schematic diagram summarizing the drivers and mechanisms of the enhanced mid-Pliocene EASM precipitation compared to the PI experiment in both simulations, as well as the differences in the mid-Pliocene EASM precipitation between the two simulations, using moisture budget analysis and MSE budget analysis, respectively AGCM and CGCM simulations. The enhanced meridional thermal contrast over the EASM domain and increased stationary eddy convergence over South China both contribute to sustain the non-uniform distribution of the mid-Pliocene EASM precipitation difference between the CGCM and AGCM simulations. The enhanced meridional thermal contrast contributes to the decreased EASM precipitation through weakening the EASM circulation and associated moisture transport. The weakening of EASM circulation is linked to the anomalous cyclone over South China induced by the warmer SST bias between the CGCM and the PRISM3 SST imposed in the AGCM. The increased stationary eddy convergence over South China provides a reasonable explanation for the enhanced precipitation associated with increased vertical motion over South China in the CGCM compared to that in the AGCM simulation. Therefore, the larger mean mid-Pliocene precipitation over the EASM domain in the CGCM simulation compared to the AGCM simulation is attributed to the enhanced vertical moisture advection; and it is the dynamic component rather than the thermodynamic component of the vertical moisture advection that contributes to the enhancement of vertical moisture advection.
Multi-model comparisons are carried out to validate the drivers and mechanisms of enhanced mid-Pliocene EASM precipitation, as revealed by IPSL-CM5A and its stand-alone LMZD5 simulations. The enhancement of mid-Pliocene EASM precipitation relative to the corresponding PI experiments in two simulations driven by surface warming over EASM domain is also reproduced by another eight climate models (FGOALS-g2, GISS-E2-R, HadCM3, HadAM3, MRI-CGCM2.3-CGCM, MRI-CGCM2.3-AGCM, CCSM4, COSMOS). The enhanced zonal thermal contrast simulated by these models has also considered to be a key mechanism controlling these enhanced mid-Pliocene EASM precipitation in both simulations through strengthening the moisture transport into EASM domain and inducing an increase in the local convergence of meridional stationary eddy velocity.
The mid-Pliocene EASM precipitation produced by other coupled models (HadCM3, MIROC4 m-CGCM, MRI-CGCM2.3-CGCM, NorESM-L and COSMOS) is also larger than those in stand-alone AGCM (HadAM3, MIROC4 m-AGCM, MRI-CGCM2.3-AGCM, CAM4 and ECHAM5) simulations. Here we explain the lager simulated EASM precipitation in coupled models than their atmospheric components by examining the dominator factors in determining these differences and their connection to warmer SST anomalies. The ∂v * ∂y has been considered a good metric to gauge the mid-Pliocene EASM precipitation differences between the two simulations. The larger simulated EASM precipitation is associated with an increase in the ∂v * ∂y to south of 20°N, despite there are large inter-model spread of the convergence center over the EASM domain. The simulated ∂v * ∂y increase matches well with the warmer SST anomalies between multi-CGCM ensemble mean and PRISM3 SST imposed in stand-alone AGCM. In present study, we emphasize the primary role of the local warmer SST anomalies over Southern China in the simulated ∂v * ∂y increase; despite that larger warm SST bias also exist in the tropical ocean. In addition, we should note the warmer SST in other oceans such as the Indian Ocean. In present climate, a warmer tropical Indian Ocean at both interannual and inter-decadal time scales would affect the EASM by forcing an Indian Ocean-Western Pacific anticyclone tele-connection pattern (Zhou et al. 2008;Wu et al. 2010;Song and Zhou 2014a, b), whether a similar mechanism works in mid-Pliocene warrants further study.