The role of the Tibetan Plateau (TP) in maintaining the large-scale overturning circulation in the Atlantic and Pacific is investigated using a coupled atmosphere–ocean model. For the present day with a realistic topography, model simulation shows a strong Atlantic meridional overturning circulation (AMOC) but a near absence of the Pacific meridional overturning circulation (PMOC), which are in good agreement with the present observations. In contrast, the simulation without the TP depicts a collapsed AMOC and a strong PMOC that dominates deep-water formation. The switch in deep-water formation between the two basins results from changes in the large-scale atmospheric circulation and atmosphere–ocean feedback over the Atlantic and Pacific. The intensified westerly winds and increased freshwater flux over the North Atlantic cause an initial slowdown of the AMOC, while the weakened East Asian monsoon circulation and associated decreased freshwater flux over the North Pacific give rise to the initial intensification of the PMOC. The further decreased heat flux and the associated increase in sea-ice fraction promote the final AMOC collapse over the Atlantic, while the further increased heat flux leads to the final PMOC establishment over the Pacific. Although the simulations were performed in a cold world, it still importantly implicates that the uplift of the TP alone could have been a potential driver for the reorganization of PMOC–AMOC between the late Eocene and early Oligocene.
The uplift of the Tibetan Plateau (TP) was a major tectonic event that had occurred throughout the Cenozoic, and its gradual growth had exerted a strong influence on the atmospheric circulation and climate (Molnar et al., 2010). Since the pioneering work of Bolin (1950), the impacts of mountain uplift on regional and global climate have been extensively investigated. Nevertheless, most studies have emphasized the role of mountain ranges on atmosphere dynamics, while quantifications of the associated impact on ocean dynamics have been rare. For example, most previous works have used atmospheric general circulation models to address regional climate effects, notably the Asian monsoon and arid environment evolutions (e.g., Ruddiman and Kutzbach, 1989; Ramstein et al., 1997; An et al., 2001; Liu and Yin, 2002; Jiang et al., 2008; Zhang et al., 2015). However, simulations have recently been applied to investigate the effect of mountain uplift in the context of the atmosphere–ocean system, and a few studies have proposed that the uplift of the Andes (Sepulchre et al., 2009) and Rocky Mountains (Seager et al., 2002) is closely linked to the evolution of oceanic circulations, including the Gulf Stream and Humboldt Current, and the El Niño–Southern Oscillation system (Feng and Poulsen, 2014). Although it has been indicated that the TP uplift affects sea surface temperatures, sea surface salinity, precipitation, and trade winds for both the Pacific and equatorial Indian Ocean (Abe et al., 2003; Kitoh, 2004; Okajima and Xie, 2007), the influence of the TP uplift on the high-latitude oceanic circulations, particularly in the North Atlantic, has rarely been explored.
Two topographic height configurations used in experiments:
The potential importance of mountain uplift in modifying the oceanic thermohaline circulation has previously been investigated. Ruddiman and Kutzbach (1989) indicated that mountain-uplift-induced changes in the North Atlantic surface circulation are expected to increase the North Atlantic Deep Water formation. In addition, Rind et al. (1997) performed the coupled model simulations with and without the TP and proposed that the TP may have a considerable impact on the large-scale meridional overturning circulation (MOC). However, the integration time used in these pioneering simulations was too short to fully evaluate the deep oceanic circulation response, and thus more studies are still needed to evaluate the possible role of the TP in modulating the MOC.
On a geological timescale, remarkable reorganization and evolution of the large-scale oceanic overturning circulation, from the Southern Ocean deep-water dominating mode to the modern-like North Atlantic deep-water mode, have been evidenced through the late Eocene to the early Oligocene (Wright and Miller, 1993; Davies et al., 2001; Via and Thomas, 2006). This dramatic shift is possibly associated with major rearrangements in the ocean seaways and other tectonic changes, although the ultimate trigger is still being debated (Zhang et al., 2011). In addition, it is suggested that the regional surface of the TP had reached a high elevation of more than 4000 m around 40 Myr ago (Dupont-Nivet et al., 2008; Wang et al., 2008), although debates regarding paleoaltitude reconstructions remain (Botsyun et al., 2016). Given this timing of TP uplift, it is important to quantify the contribution of TP uplift on the meridional oceanic circulation of the Northern Hemisphere. In this study, therefore, two coupled atmosphere–ocean numerical integrations, with and without the TP, are designed to investigate the role of the TP on the Atlantic MOC (AMOC) and Pacific MOC (PMOC).
The Community Earth System Model (CESM) version 1.0.5 of the National Center
for Atmospheric Research is a widely used, well-validated coupled model with
dynamic atmosphere, land, ocean, and sea-ice components (Gent et al., 2011).
It is applied to this study at a low-resolution configuration that is
computationally efficient and well described (Shields et al., 2012) and
employs an atmospheric horizontal grid of roughly
3.75
Two experiments are conducted; firstly, a control run with the modern
topography (MTP, Fig. 1a), and secondly a sensitivity run in which topography
within the region of 20–60
Because the major aim of this paper is to analyze changes in the
meridional oceanic circulation, we decided to focus on the density flux
parameter, which is appropriate to diagnose these oceanic circulation
changes. The dense deep-water masses are formed in the area with relatively
high surface density achieved by cooling or increasing salinity. To better
understand which processes dominated the MOC changes in simulations, it is
instructive to further analyze the time evolution of density flux budget.
Therefore, a density flux analysis method, in which the total density flux
decomposes into the haline contribution due to freshwater flux and the
thermal contribution due to heat flux (Schmitt et al., 1989), is adopted in
our study. The total density flux is calculated from a linearized state
equation of seawater as
There are evident changes in the AMOC and PMOC indices in response to the TP
uplift (Fig. 1d). With MTP, the AMOC stabilizes at around 17 Sv
(Sv
The modified AMOC and PMOC are linked to the large-scale atmospheric circulation changes. In terms of model results of the NTP run relative to the MTP run, the surface air temperature over and around the TP and in the North Pacific increased but decreased over the North Atlantic (Fig. 2a), which agrees with the previous simulations (Broccoli and Manabe, 1992; Kutzbach et al., 1993). In addition, there are intensified westerlies over the North Atlantic and weakened subtropical anticyclones and trade winds over the North Pacific (Fig. 2c); the former results from a significant increase in the meridional pressure gradient driven by a large-scale equatorward shift of air mass occupying the current position of the TP (Fig. 2c) and from a reduced drag of the orographically induced gravity waves associated with the absence of the TP (Palmer et al., 1986; Sinha et al., 2012), while the latter is derived from the weakening of zonal Eurasia–Pacific thermal contrast in the middle troposphere (not shown), especially in boreal summertime, in relation to removal of the mountains (Ruddiman and Kutzbach, 1989; Rodwell and Hoskins, 2001; Kitoh, 2004).
Regional annual mean across the North Atlantic basin at
40–70
When the TP is removed, the atmospheric moisture transport between the Pacific and Atlantic oceans undergoes a basin–basin asymmetric redistribution as a response to the large-scale circulation anomalies (Fig. 2c). In comparison with MTP results, the NTP simulation shows large amounts of anomalous westerly moisture flux transported through the lowlands of Central America to the North Atlantic, causing weak moisture convergence therein (Fig. 2d). In addition, removal of the TP leads to a significant divergence of moisture over East Asia and the western North Pacific marginal seas (Fig. 2d), which is linked to a weakened monsoon circulation and is consistent with the previous simulations using both atmospheric and coupled ocean–atmosphere general circulation models (Liu and Yin, 2002; Kitoh, 2004; Molnar et al., 2010). The anomalies of atmospheric circulation shown above are derived partially from the positive feedback caused by the changed sea surface temperature due to the removal of the TP. In particular, the weakening of both the Asian monsoon and North Pacific subtropical anticyclone in association with the TP removal are shown to be greater in the context of a coupled model as compared to that in an atmosphere-only model due to the additional ocean–atmosphere feedback (Kitoh, 2004). Thus, the atmosphere–ocean feedbacks also play an important role in maintaining the inter-basin atmospheric moisture asymmetric redistribution.
The North Atlantic and Pacific region features of annual mean
sea-ice formation rate (shading; positive stands for formation, and negative
stands for melting) and sea-ice velocity (vectors, cm s
The above changes in the large-scale atmospheric circulation markedly
decrease the total ocean density flux in the North
Atlantic (Fig. 3a, brown), supporting the trend of the AMOC (Fig. 1d). The increases in both net freshwater and wind-driven sea-ice expansion are
responsible for the initial reduction of total ocean density flux and further
induce a gradual weakening of the AMOC. In more detail, the
anomalous atmospheric circulation associated with the removal of the TP
transports more water vapor northward (Fig. 2d) over the North Atlantic
Ocean, causing more precipitation at the beginning of NTP simulation
(Fig. 3b, ca. 1101–1200 years, red line). Correspondingly, the net
freshwater flux (precipitation plus runoff minus evaporation) convergence
into the North Atlantic basin at 40–70
There is also a significant increase in the area-averaged sea-ice coverage over the North Atlantic through wind-driven processes (Fig. 3c, green). With the TP, the annual mean sea ice forms mainly in the northern and western regions of the subpolar North Atlantic, and it shifts southward and eastward when driven by cyclonic wind stress associated with the Icelandic Low, and melts in the Labrador Sea (subpolar gyre) due to a warm condition (Fig. 4a). By comparison, after removal of the TP, anomalously intensified cyclonic winds induce an anomalous eastward sea-ice velocity (Fig. 4c) and also a rapid eastward shift of the sea-ice margin (Fig. 4c). Meanwhile, the locally melted sea ice due to thermodynamics processes reduces in the southeast of Greenland (red shading, Fig. 4c), but increases in the south of Greenland (blue shading, Fig. 4c). It suggests that there is more sea ice transporting from the high latitudes into the subpolar gyre region, and the anomalous expansion of the sea-ice margin in this region primarily originates from the wind-driven eastward transportation (dynamic processes), but not the local formation (thermodynamic processes). Because of this increased sea ice through thermodynamically insulating the sea water from the freezing air, the release of sensible and latent heat into the atmosphere decreases and the density of sea water finally reduces; these processes have also been previously elucidated by Zhu et al. (2014).
The North Atlantic annual mean
Moreover, the total ocean density flux increases in the North Pacific in
response to the removal of the TP. Due to the weakened Asian monsoon
circulation and associated decrease in rainfall and runoff after lowering the
topography, the net freshwater flux received by the North Pacific decreases
by 0.08 (
The aforementioned weakening of the AMOC due to the atmospheric processes further triggers a positive atmosphere–ocean feedback loop through a reduction of northward heat transport, and a subsequent decrease in sea surface temperatures, which then allow sea ice to expand, suppressing the release of evaporating latent and sensible heat, reducing the sea water density, and further weakening the AMOC, as previously shown in Jayne and Marotzke (1999) and Zhu et al. (2014). Note that the negative effect of net freshwater becomes increasingly unimportant in comparison to the heat flux feedback associated with the latent–sensible heat changes (Fig. 3a). Finally, the thermal density flux decreases by 49 % relative to the MTP run, which substantially dominates the total density flux changes (Fig. 3a). To be specific, the annual mean total density flux and mixed-layer depth over the North Atlantic, especially around Iceland where the collapse of deep-water formation occurs, is dramatically decreased in NTP (Fig. 5d, the maximum mixed-layer depth is approximately 100 m) in comparison to that in MTP (Fig. 5a, the maximum mixed-layer depth is approximately 900 m). Moreover, this reduced total density flux over the North Atlantic is attributed more to the decreased thermal density flux associated with less latent and sensible heat release (Fig. 5e) than the changed haline density flux (Fig. 5f).
Atmosphere–ocean feedbacks also strengthen the PMOC. Due to the initial development of the PMOC mentioned in Sect. 3.1, a positive feedback (as pointed out in Warren, 1983) is initiated by the intensifying meridional oceanic circulation, which transports warmer subtropical water northward and leads to the buoyancy loss and evaporation increase (Fig. 6b). This feedback is also able to re-trigger PMOC enhancement. By comparison to the changes in the North Atlantic, both the regionally averaged sea-ice coverage (Fig. 6c) and February sea-ice margin (Fig. 4f) over the North Pacific experience a slightly northward retreat and have a relatively smaller effect on the simulated strengthening of the PMOC. Over a longer time, the thermal density flux, which is due to the loss of total heat, contributes more to the total density flux than the haline flux in relation to reduction in the net freshwater discharge (Fig. 6b). Spatially, both increased total density flux and mixed-layer depth in the North Pacific Ocean show opposite change characteristics from the North Atlantic (Fig. 7d). Correspondingly, in comparison to the MTP run, there is a widespread increase in the thermally induced density flux in the subpolar North Pacific in NTP (Fig. 7f), but with little spatial change in the haline density flux (Fig. 7f). Thus, in contrast to the results shown in the North Atlantic, the increased total heat exchange between the atmosphere and ocean due to the processes of sensible and latent heat releases (Fig. 6b) ultimately becomes a dominant factor in maintaining a vigorous PMOC by controlling the increased total density flux (Fig. 6a).
As in Fig. 3, but for the North Pacific basin at 30–70
This study investigates the effect of TP uplift on the large-scale oceanic circulation using a low-resolution version of CESM. Results show that the removal of the TP initially changes the wind-driven atmospheric moisture transport process and the wind-driven sea-ice coverage expansion process, which are responsible for the initial weakening of the AMOC. Meanwhile, the suppressed monsoonal circulation in East Asia and the western Pacific marginal seas induces the decrease in rainfall and runoff and further causes the initially increased PMOC. Moreover, the positive feedback further changes the AMOC and PMOC. In particular, the AMOC weakening can further decrease the North Atlantic sea surface temperatures, ocean–atmosphere temperature contrast, evaporation, and precipitation, and subsequently increase sea-ice coverage. These processes together cause the final changes of the AMOC and PMOC (Fig. 8).
The same as Fig. 5, but for the North Pacific (30–80
Schematic diagram about the influence of the removal of TP on the AMOC and PMOC. Vectors in gray denote the climate responses in relation to the increased wind-induced and decreased monsoonal-driven net precipitation–evaporation and wind-driven sea-ice processes. The black vectors denote the feedback processes. The bold characters A and P stand for the physical processes occurring over the North Atlantic and Pacific, respectively.
A previous study demonstrated the role of Rocky Mountain uplift on heat transport and Gulf Stream patterns in the North Atlantic (Seager et al., 2002). In this study, we focus on the most prominent long-term orogenesis occurring since the Eocene: the TP and Himalayan uplift and associated impacts on the MOC. Our results can be compared with those derived from the earlier simulations, although experimental configurations differ somewhat. It has been indicated that the removal of global mountains triggers the collapse of deep water in the North Atlantic but enables formation in the North Pacific in three different coupled models (Schmittner et al., 2011; Sinha et al., 2012; Maffre et al., 2018). The simulated weakening of the AMOC is also qualitatively consistent with recent experiments using a decreased elevation of the TP and central Asia (Fallah et al., 2016). However, only TP topography is reduced in our study, but our results are comparable with those of the earlier studies, therefore highlighting the key role that TP has played in forming the current large-scale deep oceanic circulation pattern. Nevertheless, given that all existing simulations (including ours) have used a rather coarse resolution of the coupled model configuration, it is considered that a finer-resolution model may provide a better representation of the western boundary currents and allow for a more accurate and realistic resolving of the ocean eddies, which are believed to be critically important oceanic processes that should be taken into account in the realistic simulations of the AMOC (Spence et al., 2008). In addition, the low-resolution CESM is also found to generally have a cold bias with the underestimated ocean heat transport and excessive Arctic sea ice (Shields et al., 2012), which could potentially exert modulations on the AMOC weakening. It is thus considered that investigating the response of the PMOC and AMOC to the TP uplift using an atmosphere–ocean general circulation model with a higher spatial resolution would be useful. In addition, the robust changes in the AMOC and PMOC and the associated mechanisms due to the TP uplift can be evaluated through multi-model comparison.
Based on a comprehensive analysis of modern climatological data, Warren (1983) and Emile-Geay et al. (2003) hypothesized that the present MOC (mainly occurring in the Atlantic but not in the Pacific) is determined by the large mountains, namely the Himalayas and Rockies, which induce an asymmetric distribution of wind stress and moisture transport features between the Atlantic and Pacific basins. However, previous studies have also demonstrated that the asymmetric continental extents and basin widths (basin geometries) between the two basins (Weaver et al., 1999; Nilsson et al., 2013) also play a possible key role in maintaining the present-day AMOC. Our simulations support the hypothesis proposed by Warren (1983) and highlight the significant role of the TP alone in maintaining the modern AMOC. Moreover, the similar PMOC–AMOC seesaw dynamics have also been seen in simulations (Saenko et al., 2004; Chikamoto et al., 2012; Hu et al., 2012) as well as in reconstructions (Okazaki et al., 2010; Menviel et al., 2014; Freeman et al., 2015) for the last deglaciation. Such studies have also suggested that large PMOC–AMOC seesaw modulations can be triggered by slight changes in the freshwater–salinity redistributions between the Pacific and Atlantic. Furthermore, we provide an insight that the maintenance mechanism of the PMOC in the simulation without the TP, to some extent, is the same as the AMOC in the present day. In particular, for the current North Atlantic, there is a persistently northward movement of warm and salty water mass from the tropical–subtropical Gulf Stream region into the North Atlantic and farther poleward into the Norwegian and Greenland seas, where it is exposed to very cold atmospheric temperatures and followed by a gradual cooling and in turn a higher density due to the release of substantial sensible and latent heat into the overlying cold atmosphere, which is the same as the PMOC, before eventually forming the North Atlantic Deep Water.
Our simulations have potential implications for understanding
paleotemperature reconstructions and paleoceanographic circulation
reorganization. The Earth has experienced a long-term cooling trend
throughout the Cenozoic as testified by many proxies and stacked records
(Zachos et al., 2001, 2008), in association with an increased Equator-to-pole
thermal gradient. A very important contribution to understanding the large
cooling during the Cenozoic has been determined as the drastic decrease in
atmospheric CO
Finally, this simulation is performed with the constant atmospheric CO
The model output data used in the figures are provided in the Supplement.
The supplement related to this article is available online at:
The authors declare that they have no conflict of interest.
We sincerely thank the two anonymous reviewers for their insightful comments and suggestions to improve this paper. We also thank Jiang Zhu for discussions and technical support during the writing of the paper. This work was supported by the National Natural Science Foundation of China (41421004, 41572159, and 41625018). Edited by: Arne Winguth Reviewed by: two anonymous referees