New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology
Disappearance of the last tropical glaciers in the Western Pacific Warm Pool (Papua, Indonesia) appears imminent
Edited by Michael L. Bender, Princeton University, Princeton, NJ, and approved November 5, 2019 (received for review December 27, 2018)

Significance
The glaciers near Puncak Jaya, Papua, Indonesia, the last tropical glaciers in the Western Pacific Warm Pool, have recently undergone a rapid pace of loss of ice cover and a 5.4-fold increase in the rate of thinning, augmented by the strong 2015–2016 El Niño. Ice cores recovered in 2010 record approximately the past half-century of tropical Pacific climate variability and reveal the effects of El Niño–Southern Oscillation (ENSO). It appears that the regional warming has passed a threshold such that the next very strong ENSO event, which typically exacerbates the rising freezing levels and associated feedbacks such as reduced snow cover, could lead to the demise of the only remaining tropical glaciers between the Himalayas and the Andes.
Abstract
The glaciers near Puncak Jaya in Papua, Indonesia, the highest peak between the Himalayas and the Andes, are the last remaining tropical glaciers in the West Pacific Warm Pool (WPWP). Here, we report the recent, rapid retreat of the glaciers near Puncak Jaya by quantifying the loss of ice coverage and reduction of ice thickness over the last 8 y. Photographs and measurements of a 30-m accumulation stake anchored to bedrock on the summit of one of these glaciers document a rapid pace in the loss of ice cover and a ∼5.4-fold increase in the thinning rate, which was augmented by the strong 2015–2016 El Niño. At the current rate of ice loss, these glaciers will likely disappear within the next decade. To further understand the mechanisms driving the observed retreat of these glaciers, 2 ∼32-m-long ice cores to bedrock recovered in mid-2010 are used to reconstruct the tropical Pacific climate variability over approximately the past half-century on a quasi-interannual timescale. The ice core oxygen isotopic ratios show a significant positive linear trend since 1964 CE (0.018 ± 0.008‰ per year; P < 0.03) and also suggest that the glaciers’ retreat is augmented by El Niño–Southern Oscillation processes, such as convection and warming of the atmosphere and sea surface. These Papua glaciers provide the only tropical records of ice core-derived climate variability for the WPWP.
Tropical glaciers are highly sensitive indicators and recorders of climate changes (1). Unfortunately, most tropical glaciers are currently in retreat (2⇓⇓⇓⇓–7) due primarily to recent anthropogenic atmospheric warming, although strong El Niños play an intermittent role in many regions by increasing air temperature and decreasing precipitation (8⇓–10). The warmest global sea surface and upper level atmospheric temperatures occur in the Western Pacific Warm Pool (WPWP) where the associated energy drives intense and deep convective precipitation (11). Here, we present the recent, rapid retreat of the last tropical glaciers in the WPWP by quantifying the loss of ice coverage and reduction of ice thickness over the last 8 y. Total ice loss was assessed by implementing a glacier mass balance model and future regional climate model projections for the region under different climate scenarios. We also present an ice core-derived tropical Pacific climate reconstruction from these glaciers on an interannual [El Niño–Southern Oscillation (ENSO)] timescale from 1964 to ∼2010 CE in order to help understand the mechanisms driving their observed retreat.
Study Site
The glaciers near Puncak Jaya in Papua, Indonesia (4.083°S; 137.167°E; 4,884 m above sea level), which prior to 2002 was also known as Irian Jaya, are located within an east–west-oriented mountain range on the western half of the island of New Guinea (Fig. 1A). These glaciers are the last remaining tropical glaciers in the WPWP, the heat engine for Earth’s climate system, and sit near the highest peak (Carstensz Pyramid) between the Himalayas and the Andes. At lower elevations, the island is covered mainly by tropical lowland and montane rainforest. As one of Earth’s wettest regions, many sites in Papua receive 2,500 to 4,500 mm of precipitation annually, with up to ∼12,500 mm falling on the southern slope of the central Maoke Mountains (12, 13). In general, Papua experiences wet conditions during austral summer (December to March) and relatively drier conditions, particularly on the southern part of the island, during austral winter (May to October). Interannual precipitation variability is strongly linked to ENSO, with El Niño events generally associated with dry conditions, enhanced dust, and increased biomass burning, with opposite conditions dominating during La Niña events (14). The stable oxygen (δ18O) and hydrogen (δD) isotopic ratios in precipitation on Papua are controlled by regional convective activity at daily to interannual timescales with higher δ18O values (e.g., enriched in 18O) during El Niño years (13).
Geographical and physical setting of the Papua ice fields. (A) Topography of New Guinea Island, showing the location of Puncak Jaya in Papua. (B) Papua glaciers and locations of stations with available instrumental records (DISP, GRS, and ALP) (background image: https://eoimages.gsfc.nasa.gov/images/imagerecords/91000/91716/puncak_oli_2017339_lrg.jpg). (C) Oblique aerial photo taken in June 2010 of the ENF, showing the locations of the drill sites on the west and east domes. The distance between the summits of the domes is about 600 m.
Most mountain glaciers around the world are retreating (15), and for many sites atmospheric warming has been identified as the main factor affecting ice loss (16). Likewise, the glaciers near Puncak Jaya have been retreating since the end of the most recent neoglacial period ∼1850 CE (17). Total glacier area has decreased from ∼19 km2 in ∼1850 (18, 19) to ∼1.8 km2 in 2005 (20). An ice core drilling campaign was conducted in Papua in May to June 2010 during which 3 cores were recovered from the East Northwall Firn (ENF) (Fig. 1 B and C). Two ice cores, measuring 32.13 m (D1) and 31.25 m (D1B) in depth, were drilled to bedrock and extracted from the west dome (Sumantri Peak) and a third core (D2), measuring 26.19 m, was extracted from the east dome (Soekarno Peak/Ngga Pulu) (1, 17). At the same time, a flexible accumulation “stake” composed of 15, 2-m-long polyvinyl chloride (PVC) pipes linearly connected by rope was placed in the D1 borehole (4°03′ S; 137°11′ E) (Fig. 2A and SI Appendix, Fig. S1A).
Ice thinning and retreat of Papua glaciers. (A) Schematic of connected PVC pipes illustrating ice thickness changes on the ENF in June 2010, November 2015 (gray), May 2016 (red), and November 2016 (blue). Oblique aerial photographs of the ENF taken in (B) June 2010, (C) November 2015, (D) November 2016, and (E) March 2018 show the bifurcation of the ice mass.
Results and Discussion
Ice Thickness Reduction.
Changes in ice thickness since June 2010 were determined by measuring the accumulation stake during subsequent visits to the drill site, which revealed progressive exposure of PVC sections due to surface ice loss (SI Appendix, Fig. S1 B–D). During the first site visit on November 3, 2015, 2 pipe segments and 1.49 m of rope were observed lying on the surface, indicating that the ice had thinned a total of ∼5 m at a rate of ∼1.05 m/y since 2010 (Fig. 2A and SI Appendix, Fig. S1B). During a monitoring visit on May 31, 2016, 4 pipe segments were exposed at the surface (Fig. 2A and SI Appendix, Fig. S1C), indicating an additional ice thinning of ∼4.26 m within just over 6 mo. This was likely due to the effects of the very strong 2015–2016 El Niño. The most recent measurement on November 23, 2016, revealed an additional 1.43 m of ice thinning since May 2016, confirming that, although the annual rate of thinning had slowed, it remained 2.7 times greater than the average rate from 2010 to 2015 (Fig. 2A and SI Appendix, Fig. S1D). This thinning accompanied the progressive contraction of the glacier surface area resulting in the separation of the ENF into 2 ice masses (Fig. 2 B–E). These measurements illustrate the overall rapid rate of ice thinning as well as the particularly strong impact of the temperature increase and precipitation decrease associated with the 2015–2016 El Niño.
Ice Coverage Loss.
The 2 remaining large glaciers near Puncak Jaya (Fig. 3A) experienced a dramatic loss of ice coverage from 2002 to 2018 that intensified during the 2015–2016 El Niño (Fig. 3 B and C). The new data from 2015, 2016, and 2018 are integrated with those from previous observations (17) to show that the total glacier area near Puncak Jaya has decreased linearly from ∼1850 to 2018 (Fig. 3D). From 2002 to 2015, the ice cover diminished by 1.45 km2, equivalent to a loss rate of 0.11 km2 per y. During the 2015–2016 El Niño, an additional 0.11 km2 of ice cover was lost from 2015 to 2016. In March 2018, the ice cover was only 0.458 ± 0.036 km2 or a loss of 0.09 km2 since 2016. During the period from 2002 to 2018, the glacier retreated at a rapid but steady pace (Fig. 3D, Inset).
Changes of total ice area and impact of the 2015–2016 El Niño on Papua glaciers. (A) Landsat-5 satellite image of the glaciated area on October 9, 2009. The ENF and the Carstensz Glacier are outlined in black boxes. (B) Planet’s satellite image of ENF taken on May 16, 2016, showing outlines of ice extent on June 11, 2002 (as in ref. 6); March 16, 2015; May 16, 2016; and March 28, 2018. (C) Planet’s satellite image of Carstensz Glacier taken on May 16, 2016, showing outlines of ice extent on same dates as in B. (D) The changes of total glacier area near Puncak Jaya from ∼1850 to 2018. The Inset highlights the ice area changes from 2002 to 2018. Color-coded triangles depict the coverage outlined in B and C; the gray circles and uncertainty bars are derived from other studies (17). Additional information on the satellite images in B and C are provided in SI Appendix, Table S1.
The closest available instrumental records of temperature and precipitation are located 50 to 450 m below the minimum glacierized elevation (Grasberg [GRS], Dispatch Tower [DISP], and Alpine [ALP] stations, Fig. 1B) and cover the period 1997–2016 (SI Appendix, Table S2). Although these records have a number of data gaps, they indicate a near-zero probability that average daily temperatures have fallen below the freezing point since 1997 (SI Appendix, Fig. S2). This assumption is supported by calculating freezing levels from the average temperature at the nearest station (ALP) using 2 realistic lapse rates, i.e., the mean lapse rate between GRS-DISP and GRS-ALP (4.6 °C/km) (SI Appendix, Fig. S3), and the environmental lapse rate (6.5 °C/km), to identify the remaining ice surface that lies above the freezing line. Coupling these freezing levels with 3 different hypsometric curves of the glacierized area indicate that between 0 and 10% of the remaining ice surface lies above the freezing line (SI Appendix, Fig. S4). Thus, even with steep lapse rates, this hypsometric distribution implies that prevailing temperatures are generally above freezing over the glaciers, which promotes glacier melting. Furthermore, accumulation during these warm conditions is restricted to a few hours during the night and only over the highest elevations where frozen precipitation covers a small fraction of the glacierized area (SI Appendix, Fig. S5).
Glacier Mass Balance Model and Shrinkage Projections.
The rapid rate of area loss and dramatic ∼5.4-fold increase in the rate of thinning in 2015–2016 present a scenario of unstoppable glacier demise. Depending on the interannual variability of temperature and precipitation over the next few years, the remaining ice fields of Papua, Indonesia, are very likely to disappear during the next decade. To provide a quantitative basis for this scenario, a glacier mass balance model was implemented by combining the available instrumental climatic records and dynamical/statistical downscaling of future climate model projections for the region under different climate scenarios (SI Appendix, Figs. S6 and S7 and Tables S3 and S4). We performed simulations of glacier volumetric change for the period 2017–2030 using a spatially distributed mass balance model. The model computes the difference between ablation and accumulation at daily time steps and transforms that change into volume relative to initial volume. Surface area reported in the Randolph Glacier Inventory (RGI) (21) is utilized to compute an initial volume of 0.13 km3 according to an area–volume scaling (22, 23). Ablation is calculated from daily positive temperatures while accumulation utilizes temperature thresholds to determine the amount of solid precipitation falling over the glacierized area. The model was evaluated and calibrated to reproduce the surface thinning based on the available stake measurements by using optimal combinations of lapse rates derived from the closest instrumental records of temperature and precipitation. Simulation of volumetric changes was driven by climate model projections over the Southeast Asia region (Southeast Asia Regional Climate Downscaling [SEACLID]/Coordinated Regional Climate Downscaling Experiment–Southeast Asia [CORDEX-SEA] project; http://www.ukm.edu.my/seaclid-cordex). The model predicts that under practically all scenarios glacier shrinkage will continue unabated, leading to total ice loss no later than 2026 (SI Appendix, Fig. S8).
Ice Core Records.
The ice cores recovered in 2010 from the glaciers near Puncak Jaya, the first to bedrock, provide both an absolute constraint on ice thickness and proxy evidence to examine the role of atmospheric and sea surface warming and other potential mechanisms driving glacier loss in the region. To further understand the mechanisms driving the recently observed retreat of these glaciers, analyses were conducted on the Papua ice cores to reconstruct the past tropical Pacific climate variability. The 2 longest cores drilled on the ENF (D1, 32.13 m, and D1B, 31.25 m long; located ∼5 m apart) were analyzed for δ18O, δD, and concentrations of insoluble dust and major ions (17). The D2 (26.19 m) core was not analyzed as the drilling was terminated when vertical ice layers, suggesting contorted ice flow, were encountered. Basal melting was confirmed by the 0 °C temperature at the bottom of the D1 and D1B boreholes and by meltwater encountered 80 cm above the bottom of the D2 borehole. Reproducibility of the stable isotopes in both cores is illustrated (Fig. 4 A and B and SI Appendix, Table S5). The presence of the 1964 tritium (3H) peak is detectable in a record of Papua precipitation at the Global Network of Isotopes in Precipitation [GNIP (24)] Jayapura station from 1957 to 1991 (SI Appendix, Fig. S9 A and B). This tritium peak is also recorded at 23.4-m depth in the D1 core, although its value (2.98 ± 0.42 TU) (Fig. 5A) is reduced due to radioactive decay over the 47 intervening years. This serves as the sole chronological marker for the timescale reconstruction. Modest precipitation seasonality in the Papua highlands (13) and post-depositional processes minimize seasonal aerosol variations in the cores (Fig. 4 C–E). Dating of the ice core records was challenging due to 1) this lack of seasonal variability, 2) the presence of only one chronological marker, and 3) melting of the glacier at both the surface and the base. However, construction of post-1964 chronologies for both cores was attempted by δ18O reference matching with the NINO3 extended reconstructed sea surface temperatures [ERSST (25)]. The rationale for the matching is based on a positive, albeit weak, correlation between rainfall δ18O at GNIP stations in Papua and NINO3 ERSST (SI Appendix, Figs. S10 and S11 and Table S6). The resulting timescales demonstrate that the ENF cores encompass at most 47 y (1964–2010 CE) in the top 23.4 m (Fig. 5B). The correlations between the 2 ice cores’ δ18O values and the NINO3 sea surface temperatures (SSTs) are statistically robust with r = 0.58 (P < 0.001) for the D1 core and r = 0.52 (P < 0.001) for the D1B core. The correlations are just slightly stronger for the 1964–1999 period with r = 0.61 (P < 0.001) for the D1 core and r = 0.53 (P < 0.001) for the D1B core. However, ablation at the surface of the glacier, which was observed during the 2010 drilling operation and subsequently documented by stake measurements (Fig. 2), may have truncated the top of the time series. The very high concentrations of soluble and insoluble aerosols in the upper ∼4 m of the ice cores (Fig. 4 C–E) illustrate the extent of this ablation, which may have accelerated in the early 21st century (∼2005) (Fig. 5 C–E). Therefore, this climate record since 1964 should not be regarded as precise at an annual resolution. Since melting is also occurring at the base of the glacier, the bottom ages of the D1 and D1B ice cores cannot be fixed, and thus timescale reconstruction for the cores was not attempted below the depth of the 3H signal.
Ice core records. The (A) δ18O and (B) δD in the D1 (black) and D1B (red) cores show high reproducibility and positive trends toward the surface. Comparison of (C) mineral dust, (D) chloride, and (E) potassium concentrations among the ENF cores shows reproducibility and increased concentrations toward the top as the stable isotope variability decreases and the trend of isotopic enrichment increases (gray shading).
Annual ice core records since 1964. (A) The δ18O and 3H data from 68 samples in D1 core (red closed circles with error bars). The 1964 horizon is marked by the 3H peak (2.98 ± 0.42 TU) at 23.4-m depth. (B) The D1 and D1B δ18O time series are compared, and 12 points in D1 (green closed circles) are paired with corresponding points in 13-mo running means of NINO3 ERSST (24) after assigning the 1964 3H peak at 23.4 m and assuming the top layer is May 2010. The correlations between the 2 ice cores’ δ18O values and the NINO3 SSTs are shown. (C) The logged values of annual averages of dust concentration (size of >0.63 μm) in the D1 and D1B cores post-1964. (D) As in C, but for chloride concentrations. (E) As in C, but for potassium concentrations.
An independent check on the timescale is possible by comparing the core D1 tritium data (Fig. 6A) with the record of annual tritium in precipitation from the GNIP stations in Jayapura, Indonesia, and Madang, Papua New Guinea (SI Appendix, Fig. S9A), from 1957 to 1991 (Fig. 6B). Unfortunately, tritium data from these stations do not exist after 1991. However, the comparison of these records over 35 y is comparable with the ice core δ18O matching with the NINO3 SST time series.
Comparison of core D1 and GNIP station 3H records. (A) The 3H record from core D1 is shown in depth next to (B) the annual 3H in precipitation collected from GNIP stations in Jayapura, Indonesia (blue curve), and Madang, Papua New Guinea (red curve), from 1957 to 1991 (Bottom). The green curve depicts the average 3H values from the 2 stations where the records overlap (1974–1982). The thin lines connect similar features in the ice core and precipitation curves. Note that the 3H concentration in precipitation collected at the stations has been calculated using the decay half-life of 12.3 y in order to show that the value of the peak in the sample collected in 1964 is compatible with that in the ice core record.
In the Tropics, precipitation δ18O is negatively correlated with precipitation amount at lower altitudes and in coastal regions through the amount effect (26, 27). However, this relationship is much more complex, since stable isotopes in tropical precipitation are affected by source conditions such as SSTs, along with large-scale atmospheric circulation, transport pathways, cloud-top temperature and cloud-top pressure in convective clouds, and prevailing low-latitude upper tropospheric temperatures (10, 13, 28⇓⇓⇓⇓–33). The precipitation type also appears to influence δ18O, i.e., a higher proportion of convective versus stratiform precipitation results in isotopic enrichment (34). Moreover, convective precipitation is more responsive to increasing temperature than stratiform precipitation (35). The influences of coupled ocean–atmosphere processes on stable isotopes are particularly pronounced during ENSO events, which is evident in the correlation between SSTs in the equatorial Pacific and δ18O in rainfall at Jayapura (SI Appendix, Fig. S10), which facilitated the development of the ENF ice core timescales.
Since the middle of the 20th century, the increasing tropical Pacific sea surface and upper level temperatures have been augmented by the effects of El Niño (36). The ENF ice core δ18O records show increasing trends (SI Appendix, Fig. S12A), with an average slope of 0.018 ± 0.008‰ per year (P < 0.03), which reflects the trends in global and regional surface and upper level (550 mb) temperature anomalies (SI Appendix, Fig. S12 B and C). However, from the mid-20th century to 2010, seasonal trends in precipitation anomalies were constant (SI Appendix, Fig. S12D). This also implies that the relationship between δ18O and the amount effect, which is obvious on a seasonal basis in the Tropics, is indirect on intraannual timescales. Precipitation amount and temperature are components of the complex atmospheric and oceanic processes that control the stable isotopes in the high-altitude precipitation and hence in the glacier ice.
Although the ENF ice core records contain less than a century of climate history, they serve as the only ice core-derived, high-altitude climate reconstruction from the WPWP. These ice core records not only demonstrate the deteriorating condition of these few surviving glaciers, but also provide information on the changing conditions of the atmospheric and oceanic conditions that will drive the loss of the only tropical ice fields between the Himalayas and the Andes over the next few years, especially if the loss is augmented by strong El Niños.
The glaciers near Puncak Jaya are remnants of glaciers that have existed for ∼5,000 y (37, 38), and at the current rate of ice loss these glaciers will disappear in a very short time. The precise ice coverage over several years provides a reasonable determination of the variations in the rate of ice loss between El Niño and non-El Niño years.
The climate records from the ENF ice cores establish that the current glaciers have undergone significant mass loss, are melting from the surface and at the base, and have been doing so for many years, and that the rate of melting and hence the destruction of information regarding interannual climate variability is increasing. There is a strong probability that the ENF now contains fewer years than it did in 2010 as time (ice) is being removed by ablation from the bottom and the top. Along with the Furtwängler Glacier on the top of Kilimanjaro (5), these cores illustrate what happens to a climate record contained in a glacier as it melts away.
Materials and Methods
Ice Surface Area and Thickness Changes.
The changes in ice surface area were determined using 5 satellite images to outline and quantify the surface area of the glaciers near Puncak Jaya (Fig. 3 B and C and SI Appendix, Table S1). The 2002 flat image was extracted from Google Maps (stereoscopic 3D mode turned off), and all images were obtained as orthorectified products. Changes in ice thickness since June 2010 were determined by measuring the accumulation stake during subsequent visits (Fig. 2A and SI Appendix, Fig. S1 B–D).
Mass Balance Model and Future Projections.
Volumetric changes from 2017 to 2030 were simulated using a spatially distributed model that calculates mass balance at daily time steps on each grid-cell flagged as ice or snow, according to the RGI (21). Ablation was computed from daily positive temperatures, while accumulation was calculated using temperature thresholds to determine the amount of solid precipitation falling over the glacierized area. The model was calibrated by simulating periods with available stake measurements, and identifying optimal combinations of lapse rates and degree day factors (DDFs) that reproduced the surface thinning. One thousand combinations were tested using the lapse rates shown in SI Appendix, Fig. S3 and published DDFs (39). Precipitation and temperature records around the study area (SI Appendix, Table S2) were combined to construct a relatively continuous time series for the highest elevation weather station (ALP). Simulation of future volume change was performed for the period 2017–2030 using input from 5 regional climate models over Southeast Asia [CORDEX-SEA (40)] and for 3 Coupled Model Intercomparison Project (Phase 5) experiments (historical, Representative Concentration Pathway 4.5 [RCP4.5], and RCP8.5) (SI Appendix, Table S4). Time series for the Puncak Jaya ice fields were obtained from bilinear interpolation of the gridded model output, and then downscaled using the cumulative distribution function-transform method (41).
Ice Core Analysis and Timescale Reconstruction.
The D1 and D1B ice cores were cut into discrete samples (each ∼2 to 3 cm long), melted, and analyzed for oxygen and hydrogen isotope ratios (δ18O, δD, respectively), and insoluble dust and major anion and cation concentrations. The annual means of δ18O of precipitation from the GNIP station at Jayapura, Papua (1961–1991), were compared with NINO3 ERSST. The timescale reconstruction was based on identification of the 1964 peak in tritium and matching of the core D1 tritium record with tritium in precipitation from GNIP stations at Jayapura and Madang, and on δ18O reference matching with NINO3 ERSST.
Detailed explanations of the methods, analysis, model simulations, and associated references are in SI Appendix.
Acknowledgments
D.S.P., L.G.T., E.M.-T., M.E.D., P.-N.L., J.P.N., B.W.B., V.N.M., P.G., V.Z., K.R.M., U.S., and B.G.M. were supported by National Science Foundation (NSF) Grant ATM-0823586. R.D.S. was supported by NSF Grant ATM-0823617. A.F. was supported by National Fund for Scientific and Technological Development (Fondo Nacional de Desarrollo Científico y Tecnológico) Grant 11160454. D.S.P., W.H., M.N.H., and D.G. were supported by Agency for Meteorology Climatology and Geophysics (Badan Meteorologi, Klimatologi, dan Geofisika) 2009–2010 funds. Y.K. and G.S. were supported by PT Freeport Indonesia (PTFI) funds. J.F.B., V.Z., and K.R.M. are retired. We gratefully acknowledge the considerable staff and logistical support of PTFI in Indonesia. This manuscript is Byrd Polar and Climate Research Center contribution number C-1558.
Footnotes
- ↵1To whom correspondence may be addressed. Email: donaldi.permana{at}bmkg.go.id or thompson.3{at}osu.edu.
Author contributions: L.G.T., E.M.-T., and R.D.S. designed research; D.S.P., L.G.T., E.M.-T., M.E.D., P.-N.L., B.W.B., V.N.M., P.G., V.Z., K.R.M., W.H., M.N.H., Y.K., D.G., G.S., and R.D.S. performed research; D.S.P., L.G.T., M.E.D., P.-N.L., J.P.N., J.F.B., U.S., A.F., and B.G.M. analyzed data; D.S.P., L.G.T., E.M.-T., and M.E.D. wrote the paper; D.S.P., L.G.T., V.N.M., P.G., V.Z., K.R.M., W.H., Y.K., D.G., G.S., and R.D.S. supported the ice core drilling project and collected ice core samples; D.S.P., L.G.T., K.R.M., W.H., M.N.H., Y.K., and G.S. measured the stake accumulation; M.E.D., and P.-N.L. conducted the ice core stable isotope, dust, and chemical analyses; and U.S. conducted the ice core tritium analysis.
The authors declare no competing interest.
This article is a PNAS Direct Submission.
Data deposition: The data reported in this paper have been archived at the National Center for Environmental Information (NCEI) National Oceanic and Atmospheric Administration (NOAA) World Data Center for Paleoclimatology: https://www.ncdc.noaa.gov/data-access/paleoclimatology-data/datasets/ice-core; https://www.ncdc.noaa.gov/paleo/study/24351.
This article contains supporting information online at https://www.pnas.org/lookup/suppl/doi:10.1073/pnas.1822037116/-/DCSupplemental.
Published under the PNAS license.
References
- ↵
- ↵
- L. G. Thompson et al
- ↵
- ↵
- ↵
- L. G. Thompson et al
- ↵
- ↵
- S. Haberle,
- J. Stevenson,
- M. Prebble
- M. L. Prentice,
- S. Glidden
- ↵
- ↵
- B. Francou,
- M. Vuille,
- V. Favier,
- B. Cáceres
- ↵
- L. G. Thompson et al
- ↵
- G. Chen,
- H. Qin
- ↵
- A. J. Marshall,
- B. M. Beehler
- M. L. Prentice,
- G. S. Hope
- ↵
- D. S. Permana,
- L. G. Thompson,
- G. Setyadi
- ↵
- ↵
- ↵
- ↵
- D. S. Permana
- ↵
- J. A. Peterson,
- L. F. Peterson
- ↵
- Q. van Ufford,
- P. Sedgwick
- ↵
- J. L. Kincaid
- ↵
- ↵
- D. B. Bahr,
- M. F. Meier,
- S. D. Peckham
- ↵
- I. Gärtner-Roer et al
- ↵
- IAEA/WMO
- ↵
- B. Huang et al
- ↵
- ↵
- P. K. Swart,
- K. C. Lohmann,
- J. McKenzie,
- S. Savin
- K. Rozanski,
- L. Araguás-Araguás,
- R. Gonfiantini
- ↵
- ↵
- Z. Cai,
- L. Tian
- ↵
- K. E. Samuels-Crow et al
- ↵
- M. A. Scholl,
- J. B. Shanley,
- J. P. Zegarra,
- T. B. Coplen
- ↵
- M. Vuille et al
- ↵
- J. V. Hurley,
- M. Vuille,
- D. R. Hardy,
- S. J. Burns,
- L. G. Thompson
- ↵
- ↵
- P. Berg,
- C. Moseley,
- J. O. Haerter
- ↵
- R. Seager et al
- ↵
- G. S. Hope,
- J. A. Peterson
- ↵
- E. Loffler
- ↵
- R. Hock
- ↵
- J. X. Chung,
- J. Liew,
- T. Fredolin,
- A. F. Jamaluddin
- ↵
Citation Manager Formats
Sign up for Article Alerts
Article Classifications
- Physical Sciences
- Earth, Atmospheric, and Planetary Sciences