Southern Tibetan Plateau ice core δ18O reflects abrupt shifts in atmospheric circulation in the late 1970s

Ice cores from the Tibetan Plateau provide high-resolution records of changes in the snow and ice isotopic composition. In the monsoon sector of southern Tibetan Plateau, their climatic interpretation has been controversial. Here, we present a new high-resolution δ18O record obtained from 2206 measurements performed at 2–3 cm depth resolution along a 55.1 m depth ice core retrieved from the Noijinkansang glacier (NK, 5950 m a.s.l.) that spans the period from 1864 to 2006 AD. The data are characterized by high δ18O values in the nineteenth century, 1910s and 1960s, followed by a drop in the late 1970s and a recent increasing trend. The comparison with regional meteorological data and with a simulation performed with the LMDZiso general circulation model leads to the attribution of the abrupt shift in the late 1970s predominantly to changes in regional atmospheric circulation, together with the impact of atmospheric temperature change. Correlation analyses suggest that the large-scale modes of variability (PDO and ENSO, i.e. Pacific Decadal Oscillation and El Nino-Southern Oscillation) play important roles in modulating NK δ18O changes. The NK δ18O minimum at the end of the 1970s coincides with a PDO phase shift, an inflexion point of the zonal index (representing the overall intensity of the surface westerly anomalies over middle latitudes) as well as ENSO, implying interdecadal modulation of the influence of the PDO/ENSO on the Indian monsoon on southern TP precipitation δ18O. While convective activity above North India controls the intra-seasonal variability of precipitation δ18O in southern TP, other processes associated with changes in large-scale atmospheric circulation act at the inter-annual scale.


Introduction
The Tibetan Plateau (TP) is a unique area where the topography plays a key role in the regional atmospheric circulation and the hydrological cycle, which encompasses the largest number of glaciers outside of polar regions ). Short instrumental records from inhabited areas depict a recent large warming trend (Liu and Chen 2000;Wang et al. 2008;Qin et al. 2009). Natural archives allow us to place recent changes in the broader context of natural climate variability. Among them, ice cores provide high-resolution proxy records of climate variability, archived through e.g. changes in water stable isotopes (usually expressed in a delta notation, such as δ 18 O used hereafter). However, their climatic interpretation requires us to decipher the climatic controls on high-altitude precipitation δ 18 O. While northern TP δ 18 O is mostly controlled by temperature-driven distillation, as in mid and high latitudes, the situation is more complex in the monsoon sector of the southern TP .
A number of studies have attributed recent increasing trends in south TP ice core δ 18 O to regional warming (Thompson et al. 2000;Yao et al. 2006a;Zhao et al. 2012), while process-based studies have stressed the importance of changes in precipitation amount and Indian monsoon intensity (Vuille et al. 2005;Kaspari et al. 2007;Joswiak et al. 2010). Some authors have also stressed that spatial patterns are related to changes in air mass trajectories and moisture sources (Tian et al. 2003;Johnson and Ingram 2004;Breitenbach et al. 2010).
The southern TP climate is characterized by strong seasonal shifts in atmospheric circulation, as winter air masses are advected by the westerlies while summer atmospheric circulation is dominated by the Indian monsoon (Fig. 1a). Large-scale circulation indeed plays a key role in the spatial and seasonal distribution of southern TP precipitation δ 18 O at intra-annual scales (Tian et al. 2001;Gao et al. 2011;Yao et al. 2013). Several studies have therefore suggested that TP ice cores record past variations of the regional atmospheric circulation (Qin et al. 2000;Kang et al. 2002Kang et al. , 2007b which transports moisture from different sources (Tian et al. 2006), as also proposed for other locations (Meeker and Mayewski 2002;Kreutz et al. 2000;Hoffmann et al. 2003b;Xiao et al. 2004). So far, the exact processes relating atmospheric circulation and southern TP δ 18 O have not been precisely identified and quantified.
Here, we present a new high-resolution isotopic record from an ice core recently drilled on the Noijinkansang (NK) glacier, southwest of Lhasa, which covers 2.7 km 2 in the southern TP (Fig. 1a) and extends from 5540 to 6674 m a.s.l.. During May 2007, our ice core was drilled to bedrock (55.1 m depth) at 5950 m a.s.l. in a flat basin of the accumulation zone of the NK glacier (29°2′8.68″N, 90°11′54.10″E, Fig. 1a). The nearest meteorological stations, Nagarze (29.96°N, 90.33°E, 4430 m a.s.l.) and Lhasa (29.70°N, 91.13°E, 3650 m a.s.l.), show that more than 80 % of annual precipitation occurs during the summer Indian monsoon season (June-September, hereafter JJAS, Fig. 1b), with a large inter-annual variability (Gao et al. 2009;Tian et al. 2001) (Fig. 1b). Monthly precipitation δ 18 O from these stations displays a clear seasonal cycle, typical of southern TP, with maximum values in spring and minimum values in summer (Fig. 1b, Gao et al. 2009Yao et al. 2013). This seasonal pattern is identified in the NK ice core and provides the basis for its age scale (Sect. 2).
This paper focuses on the climatic interpretation of the NK ice core δ 18 O record (hereafter NK δ 18 O). In Sect. 3, we explore the statistical relationships between interannual δ 18 O variations and local or regional temperature and precipitation data. Large-scale atmospheric teleconnections with NK δ 18 O are also investigated. Atmospheric reanalyses, climate indices and an atmospheric general circulation model (LMDZiso) are used to identify and quantify the role of changes in large-scale atmospheric circulation on NK δ 18 O. Section 4 reports our conclusions and perspectives.    of surface air temperature (T, °C, circles), precipitation amount (P, mm per month, vertical bars) at Lhasa (LS) and Nagarze (LKZ) meteorological stations, and monthly weighted precipitation δ 18 O (‰, 1994-2006, triangles) at Lhasa. Error bars show the interannual standard deviation of air temperature and precipitation amount

Data
The NK ice core was cut into 2206 samples with intervals of 2 or 3-cm for analyses of stable isotopes, major soluble ions and carbonaceous aerosols, and 55 samples (~1 m resolution) for total beta radioactivity (hereafter β-activity) measurements. All the samples were transported to laboratory in frozen condition and melted at room temperature prior to analyses. Oxygen isotopes were measured using a MAT-253 mass spectrometer (with an analytical precision of 0.05 ‰) at the Key Laboratory of CAS for Tibetan Environment and Land Surface Processes, China. The β-activity was measured at the Analytical Laboratory of the Beijing Research Institute of Uranium Geology, China. The ice core is drilled on a vast and flat platform in the accumulation zone of the glacier, close to the ice divide. Field observations report no sign of surface wind mixing. Zhao et al. (2012) describes in details the ice core drilling and analytical procedures.
Below 16 m depth, the ice core density gradually increases to 0.9 g/cm 3 , with no indication of melt or stratigraphic perturbation. The chronology is based on the identification of seasonal cycles in δ 18 O and major soluble ions and verified using β-activity. The seasonal cycle of precipitation δ 18 O is used to identify annual layers from one maximum of ice core δ 18 O to the next one, therefore from spring (April-May) of 1 year to spring of the next year. The mean accumulation rate of the NK ice core of 470 mm w.e. is larger than the climatological precipitation amount at Lhasa and Nagarze (441.3 and 359.3 mm, respectively). Peaks in β-activity are detected at ~13.73 m (1987 Chernobyl event) and ~23.59 m (1963 nuclear bomb testing). The 55.1 m ice core is estimated to span 142 years (Fig. 2c). From repeated  1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 1900 1910 1920 1930 1940 1950 1960 1970 1980 1990   . From low-resolution measurements of δ 18 O and δD (220 merged samples), Zhao et al. (2012) observed an increasing trend of deuterium excess. They attributed this signal to a decrease of Indian monsoon precipitation on the southern TP through the reduced fraction of annual mean accumulation provided by summer precipitation (associated with a low deuterium excess, Zhao et al. 2012). Here, we expand the investigation of the NK water stable isotope record through our new high-resolution δ 18 O data.
This new record will be compared in the next section with local and regional meteorological data. For this purpose, surface air temperature data are available from Lhasa and Nagarze meteorological stations (National Meteorological Information Center of the China Meteorological Administration). We also use 2.5° × 2.5° regional monthly precipitation amount data spanning 1979-2006 from the GPCP (Global Precipitation Climatology Project) (http:// www.esrl.noaa.gov/psd/data/gridded/data.gpcp.html). Finally, we also use monthly regional surface air temperature and vertical winds data for the period of interest from NCEP reanalysis data (http://www.esrl.noaa.gov/psd/data/ reanalysis/reanalysis.shtml).
To explore the impacts of El Nino/Southern Oscillation (ENSO) on NK δ 18 O, we use the Nino 3.4 sea surface temperature (SST) index derived from the HadISST1 data set (Rayner et al. 2003) for the period 1870-2006. The relationships between our NK δ 18 O and Pacific Decadal Oscillation for the period of 1900-2006 (PDO, http://jisao.washington.edu/pdo/PDO.latest, Zhang et al. 1997;Mantua et al. 1997) as well as Zonal index (defined by the normalized difference in zonal-averaged sea level pressure anomalies between 35 N and 65 N, representing the overall intensity of the surface westerly anomalies over middle latitudes) for the period 1864-2006 (http://ljp.lasg.ac.cn/dct/page/65592, Li and Wang 2003) are also analyzed.
LMDZiso is an atmospheric general circulation model which explicitly includes water stable isotopes (Risi et al. 2010). It captures the seasonal and intra-seasonal variability of TP precipitation δ 18 O (Gao et al. 2011;Yao et al. 2012;Gao et al. 2013), but with a cold and dry bias in the southern TP (Gao et al. 2011). Here, we use the interannual precipitation and vapor δ 18 O outputs from a simulation performed at 2.5° × 3.75° resolution, run from 1958 to 2006 in a "nudged" mode. The Atmospheric Model Intercomparison Project (AMIP) protocol (Gates 1992) is used to prescribe monthly sea surface boundary conditions. The European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis results (Uppala et al. 2005) are used to nudge the horizontal wind fields with a relaxation time of 1 h. This warrants that the large-scale LMDZiso simulated wind fields are consistent with the observations, and allows to use of the LMDZiso framework to understand and quantify the role of changes in atmospheric circulation in the observed δ 18 O variability. A zoomed simulation using the zoom functionality of the LMDZ model, stretching the grid to the horizontal resolution of 50-60 km between 0° and 55°N and 60° and 130°E, are also used here to detect the inter-annual variability.

Ice core δ 18 O variation
The NK δ 18 O record shows large seasonal variations, from −7.0 ‰ at 24.19 m depth to −24.7 ‰ at 16.27 m depth (not shown), with an annual mean value of −16.9 ± 1.7 ‰ from 1864 to 2006 (Fig. 2c). This record contains remarkable multi-decadal variations. Utilizing the regime shift detection algorithm of Rodionov (2004), six different regimes are identified with the regime shift years of 1911, 1921, 1957, 1973 and 1980 and in 1976-1977 in the Guliya, Geladaidong and Tanggula ice cores (not shown). Because the dating uncertainties of these ice cores are usually around ±1 year (Hou et al. 2003;Yao et al. 2006a, b;Joswiak et al. 2010;Kang et al. 2007a;Yao et al. 1997;, we conclude that there is a widespread negative δ 18 O anomaly in late 1970s in the Tibetan ice cores, most strongly expressed in NK and Tanggula ice cores.

Relationships with regional meteorological data
These results motivate further exploration of the climatic controls on NK δ 18 O, with a focus on the period from 1960 to 2006, where the age scale is well established. We first investigate statistical relationship of NK δ 18 O with regional meteorological data. Since more than 80 % of precipitation occurs in JJAS, we focus on the relationships between JJAS temperature and precipitation amount and NK δ 18 O. Since the 1960s, NK δ 18 O shows decadal variations which parallel those recorded in Nagarze and Lhasa JJAS temperature.

Comparison with LMDZiso simulations and large scale climate fields
We now compare these data with LMDZiso outputs. The simulated temperature is in better agreement with the observed temperature at Nagarze station than that at Lhasa (Fig. 2e, f), especially for the cool periods (1968,(1977)(1978)(1979)(1992)(1993)(1994)(1995)(1996)(1997). However, the model overestimates the variability of precipitation amount, compared with observations at Lhasa and Nagarze stations (Fig. 2g, h), especially before 1980s. We observe a better model-data agreement between NK δ 18 O and the simulated water vapor δ 18 O than with the simulated precipitation δ 18 O for the period of 1958-2006 (Fig. 2a-c); this could be related to the intermittency of precipitation, or to post-deposition processes (Steen-Larsen et al. 2014). Some inter-annual changes in the simulated annual temperature, precipitation amount and water vapor δ 18 O coincide with changes recorded in NK δ 18 O (Fig. 2a-c, f, h). Within LMDZiso, the percentage of the simulated precipitation δ 18 O which can be explained by a linear correlation with local simulated temperature and precipitation is weak and can only explain less than 10 % of precipitation δ 18 O variance (Fig. 2b, f, h), even when considering a multiple regression that includes the temperature and precipitation amount together. The lowest NK δ 18 O is observed and simulated in 1978, and it is associated with the most negative anomalies of simulated annual temperature (Fig. 2e), precipitation amount (Fig. 2h) and vapor δ 18 O (Fig. 2a).
NK δ 18 O is weakly positively correlated with JJAS simulated precipitation amount in the tropical Indian Ocean (R < 0.6, p < 0.1) and, as in GPCP, anti-correlated with LMDZiso precipitation in the southeast TP (R > −0.3, p < 0.1) (Fig. 3c). NK δ 18 O is most strongly correlated with vapor δ 18 O at 425 hPa for the period of 1979-2006 (pressure level closest to NK ice core drilling site elevation) in the tropical Indian Ocean (10°S-10°N, 85°E-130°E) and the west Pacific Ocean (Fig. 3d). In addition, NK δ 18 O is significantly positive correlated with JJAS zonal wind (U wind) at 500 hPa in 32°N-40°N and 50°E-80°E from 1979 to 2006 (Fig. 3e). Meanwhile, NK δ 18 O is significantly negatively correlated with JJAS U wind at 500 hPa over the tropical Indo-Pacific Ocean during the same period (Fig. 3e). After 1979, the significant positive correlation between NK δ 18 O and meridional wind (V wind) exists in northwest TP (Fig. 3f). NK δ 18 O is also significantly positively correlated with the observed SST since 1982 in the tropical Indo-Pacific Ocean (between 5°N and 20°S, Fig. 3i). In winter, the significantly positive correlation only exists in 0°-10°N and 95°E-120°E with U wind, and the Persian Gulf as well as a small part of Arabian Sea with V wind (Fig. 3f, h). NK δ 18 O show significantly positive correlation with winter SST in the western and northern Arabian Sea and the Philippines (Fig. 3j).
These statistical analyses suggest large-scale atmospheric circulation drivers of NK δ 18 O inter-annual variability. The mechanisms relating large-scale atmospheric circulation and tropical precipitation isotopic composition are complex, and have been explored in atmospheric models. Following earlier works of Hoffmann (2003) and Vuille et al. (2003) for the Andes, Brown et al. (2006) identified strong imprints of monsoons, meridional migration of the ITCZ and ENSO in tropical precipitation isotopic composition, through processes affecting rainout intensity and convection height. Other ice cores in the central TP, Arctic, Antarctic, and subpolar regions were also interpreted to depict rapid changes in atmospheric circulation after the mid-1970s (Kang et al. 2002;Kaspari et al. 2007;Schneider et al. 2012;Kelsey et al. 2012).
The statistical analyses show that local temperature or precipitation account for less than 50 % of NK δ 18 O variability in both observations and LMDZiso, and point to processes associated with moisture transport and largescale atmospheric circulation, at the inter-annual scale. At the intra-summer scale, different controls were identified for Lhasa precipitation δ 18 O variability . Combining daily precipitation isotopic composition measurements and satellite information on convection and water vapor isotopic composition, He et al. (2015) showed a dominant control through convective activity above North India: enhanced convection (and increased precipitation above N. India) depletes surface water δ 18 O, and this depletion is further amplified along the orographic uplift of moister air, and finally by enhanced subsidence of depleted vapor towards Lhasa. Our study suggests that controls at the intra-seasonal time scale cannot be extrapolated at the inter-annual scale. Indeed, the statistical analyses show that, at the inter-annual scale, enhanced JJAS precipitation above the southern Bay of Bengal and enriched vapor δ 18 O Fig. 3 Spatial distribution of linear correlation coefficients between NK ice core δ 18 O and regional air temperature and precipitation from observations and LMDZiso simulations, as well as winds and SST. a NCEP gridded JJAS surface air temperature data , b GPCP gridded JJAS precipitation data , c JJAS precipitation amount simulated from LMDZiso . d Correlation coefficients between NK ice core δ 18 O and LMDZiso JJAS vapor δ 18 O at 425 hPa hight . e Correlation coefficients between NK ice core δ 18 O and NCEP gridded JJAS zonal wind (U wind) at 500 hPa height . f Correlation coefficients between NK ice core δ 18 O and NCEP gridded JJAS meridional wind (V wind) at 500 hPa height . g Correlation coefficients between NK ice core δ 18 O and NCEP gridded DJF zonal wind at 500 hPa height . h Correlation coefficients between NK ice core δ 18 O and NCEP gridded DJF meridional wind at 500 hPa height . i Correlation coefficients between NK ice core δ 18 O and NCEP gridded JJAS SST (1982SST ( -2006. j Correlation coefficients between NK ice core δ 18 O and NCEP gridded DJF SST (1982SST ( -2006. On all maps, only correlations significant at the 95 % confidence level are shown in the tropical ocean and west Pacific Ocean at 425 hPa level are consistent with an δ 18 O enrichment at NK high elevation site, responding to a common large-scale mode of variability. We now discuss the mechanisms which may affect the recent NK δ 18 O trend, since 1980. While warming is expected to induce enrichment in the NK δ 18 O, one process is expected to compensate for temperature effects and weaken the NK δ 18 O increase after 1980. Enhanced westerlies (Yuan et al. 2014), accounting for a larger proportion of annual precipitation, also provide depleted moisture though long distance transport.

Causes of the late 1970s NK δ 18 O shift: statistical analysis
In order to further explore the processes controlling interannual variations of NK δ 18 O, we focus on the sharp decrease observed in the late 1970s (Fig. 2c). Meteorological data show cold conditions in 1978 at both stations ( Fig. 2e) as well as in the central TP (Kang et al. 2007b), with the coldest summer encountered in 1976 at Nagarze and Lhasa (Fig. 2d).
During the mid-1970s, shifts in North Pacific climate regime were associated with a deepening of the Aleutian Low and a strengthening of the mid-latitude westerlies (Miller et al. 1994;Nakamura et al. 1997;Mantua et al. 1997;Zhang et al. 1997;Wu et al. 2005), which affected Indian monsoon strength (Trenberth and Hurrell 1994;Wang and Ding 2006). This regime shift also involved changes in the PDO and ENSO (Guilderson and Schrag 1998;Mantua and Hare 2002;Alexander et al. 2002). The ENSO influence on TP moisture and ice core δ 18 O are also found at Dasuopu (Brown et al. 2006) and other tropical ice cores Hoffmann 2003a;Bradley et al. 2003). This motivates an investigation of the relationships between NK δ 18 O and Pacific modes of variability.
The spectral analysis of NK δ 18 O, performed with a wavelet and Multi-Taper method (not shown) reveals a significant periodicity at 2.7 years, also suggesting an ENSO influence on NK δ 18 O. During several intervals, NK δ 18 O is significantly correlated with Nino 3.4 SST index for different seasons. During the cold PDO phase of 1943-1976 (Fig. 4e), NK δ 18 O is positively correlated with December-February (DJF) Nino 3.4 SST anomalies (R = 0.32, p = 0.06 and up to 0.44, p = 0.01 when excluding 1973, Fig. 4a, b). Since 1980, during the warm phase of PDO, NK δ 18 O is anti-correlated with May-September (MJJAS) Nino 3.4 SST anomalies (R = −0.32, p = 0.11 and up to R = −0.45, p = 0.02 excluding 1984) (Fig. 4a, c, e). The NK δ 18 O minimum at the end of the 1970s coincides with a PDO phase shift (from negative phase to positive phase, Grigholm et al. 2009), an inflexion point of the zonal index (onset of a fast increase, Fig. 4a, d) as well as ENSO year (increase in frequency and intensity). This again indicates that the most remarkable NK δ 18 O negative anomaly is associated with changes of the large-scale atmospheric circulations, although a partial role of temperature cannot be excluded. The mechanism by which the PDO/ENSO could affect the southern TP ice core δ 18 O was hypothesized through the Indian monsoon and westerlies changes. An "atmospheric bridge" connects the ENSO/PDO-driven large-scale atmospheric circulation changes with the SST, wind, and humidity from the North Pacific to Indian Oceans (Alexander et al. 2002), which then alter the strength and intensity of the Indian monsoon. In the decade prior to the climate shift, SST changed in the western subtropical Pacific, leading to a shift of atmospheric circulation and changes of SST in other parts of Pacific, then driving ENSO variability by modulating the trade winds (Barnett et al. 1999;Vimont et al. 2001Vimont et al. , 2003Alexander et al. 2002). After 1976, the cooling of SST of central Pacific and subtropics associated with central Pacific El Niño forced changes in the extra-tropical atmospheric circulation (Lorenzo et al. 2010), and surface wind anomalies switched from easterly to westerly. In return, tropical SST changes associated with the recurrence of El Niño from the mid 1970s to the late 1980s reinforced north Pacific changes and westerly wind anomalies (Wu et al. 2005). The warm phase of PDO during the boreal winter season affects subtropical and tropical SST persisting into the next monsoon season (Krishnamurthy and Krishnamurthy 2013). This pattern explains the anti-correlation between NK δ 18 O and May-September (MJJAS) Nino 3.4 SST anomalies since 1980. They drive changes in equatorial trade winds that reinforce the equatorial Walker circulation and lead to enhanced ascending motion in the central Pacific and the equatorial Indian Ocean, weakening the Indian monsoon (Krishnamurthy and Krishnamurthy 2013).
With the increase of SST and zonal wind in the tropical Indo-Pacific Ocean, the evaporation of surface ocean water strengthens, providing local vapor which is not depleted in δ 18 O. The increased local evaporation replenishes the monsoon flow and may explain the reduced depletion of NK δ 18 O. In the southern TP, these processes can be compensated by the fact that a weaker Indian monsoon also leads to an increase in moisture transported from western remote sources (e.g. Persian Gulf), which is characterized by low δ 18 O. El Niño events further reduce the strength of Indian monsoon through the descending Hadley branch, while the situation is more ambiguous during La Niña events (Mantua and Hare 2002;Krishnamurthy and Krishnamurthy 2013;Sano et al. 2013;Miyasaka et al. 2014). The transport of vapor over the tropical Indo-Pacific Ocean changes with winds corresponding to the PDO/ENSO phase. In El Nino winter, the westerlies decreases due to the warmer SST in tropical Indian Ocean, resulting in more water vapor transport to north; after 1980, during the warm phase of PDO, cyclones prevail in the west of Philippines in La Nina winter, helping stronger moisture transport from Indian Ocean to continent; before 1980, the condition is opposite (Yuan et al. 2014). The influences of El Nino/La Nina on anticyclone/cyclone close to the Philippines however are inconsistent and asymmetric . Meanwhile, the deepening of the Aleutian Low and strengthening of the winter mid-latitude westerlies are associated with colder conditions (Miller et al. 1994;Alexander et al. 2002;Mantua and Hare 2002), further depleting precipitation δ 18 O.

Causes of the late 1970s NK δ 18 O shift: analysis using LMDZiso
The role of large-scale circulation is now explored using LMDZiso JJAS precipitation δ 18 O. The model simulates too enriched precipitation δ 18 O, a feature common to other atmospheric models ), but also underestimates the magnitude of inter-annual δ 18 O variations, and does not produce decadal variability. We note that the model caveats could be linked to the atmospheric model resolution, as it cannot resolve small-scale topographic features which may play a key role on the inter-annual variability of high elevation precipitation δ 18 O (Gao et al. 2011). The zoomed simulation obviously improves the intraannual precipitation δ 18 O, but not its inter-annual variability. We also stress the fact that the ice core record can be associated with a noticeable noise level, due to post deposition processes such as wind scouring, and does include a signal from the non-monsoon season. However, LMDZiso does produce a sharp decrease in NK precipitation δ 18 O in the late 1970s, and minimum values in 1978. We therefore use the framework provided by the LMDZiso simulation to understand the processes controlling the JJAS δ 18 O depletion in 1978-1979, compared to 1982-1983). The associated precipitation δ 18 O anomaly has a magnitude of ~1.5 ‰ (Fig. 5a), extends over all the TP and is also depicted in the vapor δ 18 O (~1.7 ‰, Fig. 5c).
We use a decomposition of the associated δ 18 O anomaly into several terms, which account for the effects of rainvapor disequilibrium, fractionation coefficients, local temperature, temperature at the source, local relative humidity, and residual effects including convection and large-scale transport (See Appendix in Supplementary Material). This approach shows a key role of changes in the intensity of Rayleigh distillation due to changes in temperature. In 1978-1979, LMDZ simulates TP regional cooling, consistent with meteorological data, which explains a depletion of about 0.7 ‰. Residual effects including large-scale transport however account for an additional depletion of about 1 ‰.
The pattern of changes in atmospheric circulation (Fig. 5c) depicts a reduced northward monsoon flow, and a large North-Easterly anomaly of −0.5 m/s (Fig. 5b, c). Combined with the average JJAS vapor δ 18 O South West-North East gradient of about −2.4 × 10 −3 ‰ per km (Fig. 5b) and a residence time of water vapor of about 10 days, this anomalous circulation pattern can account for the residual anomaly of low-level water vapor δ 18 O anomaly of ~−1 ‰. The precipitation δ 18 O composition then follows tightly that of the low-level water vapor, variations in precipitation-vapor differences accounting for less than 0.2 ‰ of precipitation variations.
The large-scale circulation effects are clearly visible in the latitudinal and vertical structure of JJAS zonal and meridional winds from NCEP reanalyses along a transect at 90°E for 1978 and 1983 (Fig. 5d, e). Regarding westerly wind (U, Fig. 5d), the data depict a decrease of low level winds (<500 hPa) at 0°-20°N, and an increase of high level winds (500 hPa) at 20°N-30°N in 1983 compared to 1978. The meridional profile displays a decrease of the low atmosphere southerly winds in 1983 compared to 1978 ( Fig. 5e). This illustrates the decrease of Indian monsoon and increase of westerlies in 1983 compared to 1978, consistent with the statistical analysis and the tele-connections evidenced in Sect. 3.4.
In summary, the abrupt depletion of NK δ 18 O in 1978 is captured by the LMDZiso simulation. It is mostly attributed to a depletion of water vapor δ 18 O due to the change in temperature and in large-scale atmospheric circulation. The effect of large-scale atmospheric circulation on vapor δ 18 O is found to be larger than that of temperature.

Conclusions and perspectives
In this manuscript, we have reported the δ 18 O signal recorded by a new Southern Tibetan ice core. In contrast with earlier studies, our record does not show unprecedented warming during the last decades. It is dominated by decadal variations, including a remarkable isotopic depletion during the late 1970s. Statistical analyses stress a significant anti-correlation with precipitation amount  1978-1979 to 1982-1983 using the LMDZiso model. a Decomposition of the precipitation δ 18 O anomaly into a sum of five contributions due to rainvapor disequilibrium (deseq), fractionation coefficients (α), temperature (T), local relative humidity (h), and residual effects (residual) including convection and large-scale transport (see Appendix in Supplementary Material). b Map showing the 1982-1983 pattern of JJAS vapor δ 18 O at the lowest atmospheric model level (colors) and pattern of JJAS winds at the lowest atmospheric level (vectors). c Same as (b) but for the difference between 1982and 1983and 1978-1979. Vertical profiles of difference between 1983and 1978 and meridional (e) winds along 90°E longitude from 30°S to 60°N, from NCEP reanalyses, downloaded from http://www. esrl.noaa.gov/psd/data/gridded/ data.ncep.reanalysis.derived. pressure.html. The grey shading depicts the TP topography. The blue square shows the location of NK ice core (p < 0.05) in the surrounding areas, as well as zonal wind at 500 hPa in the tropical Indo-Pacific oceans, and a positive correlation with precipitation amount in the Bay of Bengal, as well as with surface temperature in TP, India and summer tropical SSTs. Significant relationships are detected with large-scale modes of variability (ENSO and PDO). The importance of large-scale atmospheric circulation is further supported by the ability of the atmospheric general circulation model LMDZiso to resolve the 1978 isotopic anomaly. The model shows that this anomaly extends over all of TP, and is positively related to an anomaly in water vapor δ 18 O in the tropical Indo-Pacific region. Within the model framework, we isolated the driver of this anomaly, which is mostly due to an abnormal atmospheric circulation, marked by enhanced westerlies, and a reduced monsoon flow, while air temperature also affects this anomaly through changes in Rayleigh distillation, but with a smaller effect than that of circulation alone.
This work demonstrates the ability of TP annual ice core δ 18 O to record abrupt shifts in atmospheric circulation and provides a methodology to evaluate the magnitude of local versus large-scale controls on ice core δ 18 O. This will motivate further integration of high-resolution precipitation δ 18 O records from multiple archives (ice cores, but also tree ring cellulose, speleothems and lake sediment records) in order to map the spatial and temporal variability of δ 18 O, and to characterize past changes in regional atmospheric circulation.