We would like to thank you for reviewing our manuscript and providing comments about our work. We are happy to adopt many of the comments and will revise the manuscript as requested. 

—Line 51: A brief explanation is that the highlighted sentences say that simulated mPWP climate change does not include aerosol forcing, while our analysis mainly focuses on aerosol effects in mPWP simulations. The first highlighted sentence is trying to say that traditional mPWP simulations usually set prescribed aerosol concentrations same as the control runs, which means that they do not consider aerosol forcings. Our work uses three simulations: a pre-industrial control simulation (PI), a simulation with mPWP boundary condition plus pre-industrial aerosol (Plio_Pristine) and one with mPWP boundary condition plus industrial pollutants (Plio_Polluted). The simulated mPWP change from the highlighted sentence is equivalent to "Plio_Pristine - PI" in our study (Sect. 3.1). We then compares the simulated mPWP change (Plio_Pristine - PI) with the reconstructions in this section, in purpose of firstly evaluating our simulated mPWP climate and secondly pointing out the existence of data-model mismatch. Plio_Pristine - Plio_Polluted considers effects of aerosol forcing in the mPWP simulations (Sect. 3.2 and 3.3). Idealised aerosol scenarios were chosen, because the realistic mPWP aerosol concentrations are not available. Plio_Polluted uses industrial pollutants (Lamarque et al., 2010), which is a real existed aerosol forcing though not happened during the mPWP. Analysis related to aerosol effects does not touch data-model comparison. We will rewrite line 51 to "The discrepancy between mPWP simulations and reconstructions (Haywood et al., 2013, 2020) implies that models might miss some important mechanisms (Fedorov et al., 2013) or prescribed forcing could be a source of uncertainty (Feng et al., 2019). The mPWP simulations now use modern-day or pre-industrial aerosol concentration same as the control runs that may differ from the conditions during mPWP. Lack of including the effect of aerosol forcings and the usage of unrealistic prescribed aerosol concentration in mPWP simulations implies that aerosol effects may be one of the possible explanation for the mismatch between reconstructions and simulations."

—Aerosol distribution: Plio_Pristine and Plio_Polluted simulations in our study use prescribed aerosols, which give aerosol distribution instantaneously at the beginning and then left the model to respond to it. Aerosols do not distribute globally as the lifetime of aerosols is only a few weeks. The atmosphere is not well-mixed and the North-South difference emissions persists into concentrations. 

—Lines 169-171: The mPWP boundary conditions cause more warming than removing anthropogenic aerosols. That’s what conclude in line 171 as "...This suggests that changes in boundary conditions are relative more important than aerosol effects in mPWP warming." We are trying to say that removing aerosols only causes more over Northeastern Pacific and high-latitude North Atlantic Ocean than the mPWP boundaries. We will rewrite the sentence to "However, the mPWP boundary conditions causes more warming than removing anthropogenic aerosols nearly global except Northeastern Pacific and high-latitude North Atlantic Ocean (Fig. 5a,c,e). More warming lead by the mPWP boundaries than removal of aerosols suggests ..." to emphasise that the mPWP boundaries lead to more warming than removal of aerosols. 

—Line 176: We will add a figure and discussion about changes in the Hadley circulation.

—Lines 180-183: The temperature along the Equator in the Pacific Ocean is greater than surroundings. The warming band is largely overlap with the precipitation increase belt. 

—Lines 193-197: The evidence (the citations in sentences) in the second half says that the cooling effect induced by in- creased aerosol during 1950s-1980s and over the 20th century lead to monsoon reductions. Logically, monsoon would be enhanced in the opposite situation (i.e. reducing aerosols), which matches the first half sentence saying that removing aerosols result in monsoon enhancement. We will add this logic in text. 

—Lines 218-220: Sorry we made a mistake here. We will rewrite the sentence to "Over tropical and subtropical regions, 40 the overall uniform warming induce by mPWP boundaries (Plio_Pristine - PI) result in general decrease in sea level pressure and have relatively small effects on surface wind (Fig. 9a,c,e). Though removal of anthropogenic aerosols also show warming in tropical and subtropical regions, it rises the sea level pressure and have much stronger effects on surface wind that shows seasonal variances (Fig. 9b,d,f)."

—Technical corrections: We will make corrections as required. 

References 

Fedorov, A. V., Brierley, C. M., Lawrence, K. T., Liu, Z., Dekens, P. S., and Ravelo, A. C.: Patterns and mechanisms of early Pliocene warmth, Nature, 496, 43–49, https://doi.org/10.1038/nature12003, 2013. 

Feng, R., Otto-Bliesner, B. L., Xu, Y., Brady, E., Fletcher, T., and Ballantyne, A.: Contributions of aerosol-cloud in- teractions to mid-Piacenzian seasonally sea ice-free Arctic Ocean, Geophysical Research Letters, 46, 9920–9929, https://doi.org/https://doi.org/10.1029/2019GL083960, 2019. 

Haywood, A. M., Hill, D. J., Dolan, A. M., Otto-Bliesner, B. L., Bragg, F., Chan, W.-L., Chandler, M. A., Contoux, C., Dowsett, H. J., Jost, A., Kamae, Y., Lohmann, G., Lunt, D. J., Abe-Ouchi, A., Pickering, S. J., Ramstein, G., Rosenbloom, N. A., Salzmann, U., Sohl, L., Stepanek, C., Ueda, H., Yan, Q., and Zhang, Z.: Large-scale features of Pliocene climate: results from the Pliocene Model Intercomparison Project, Climate of the Past, 9, 191–209, https://doi.org/10.5194/cp-9-191-2013, 2013. 

Haywood, A. M., Tindall, J. C., Dowsett, H. J., Dolan, A. M., Foley, K. M., Hunter, S. J., Hill, D. J., Chan, W.-L., Abe-Ouchi, A., Stepanek, C., Lohmann, G., Chandan, D., Peltier, W. R., Tan, N., Contoux, C., Ramstein, G., Li, X., Zhang, Z., Guo, C., Nisancioglu, K. H., Zhang, Q., Li, Q., Kamae, Y., Chandler, M. A., Sohl, L. E., Otto-Bliesner, B. L., Feng, R., Brady, E. C., von der Heydt, A. S., Baatsen, M. L. J., and Lunt, D. J.: The Pliocene Model Intercomparison Project Phase 2: large-scale climate features and climate sensitivity, Climate of the Past, 16, 2095–2123, https://doi.org/10.5194/cp-16-2095-2020, 2020. 

Lamarque, J. F., Bond, T. C., Eyring, V., Granier, C., Heil, A., Klimont, Z., Lee, D., Liousse, C., Mieville, A., Owen, B., Schultz, M. G., Shindell, D., Smith, S. J., Stehfest, E., Van Aardenne, J., Cooper, O. R., Kainuma, M., Mahowald, N., McConnell, J. R., Naik, V., Riahi, K., and Van Vuuren, D. P.: Historical (1850-2000) gridded anthropogenic and biomass burning emissions of reactive gases and aerosols: Methodology and application, Atmospheric Chemistry and Physics, 10, 7017–7039, https://doi.org/10.5194/acp-10-7017-2010, 2010. 

Please see attached for a response.