Timing of glacier advances and climate in the High Tatra Mountains (Western Carpathians) during the Last Glacial Maximum

During the Last Glacial Maximum (LGM), long valley glaciers developed on the northern and southern sides of the High Tatra Mountains, Poland and Slovakia. Chlorine-36 exposure dating of moraine boulders suggests two major phases of moraine stabilization, at 26–21 ka (LGM I — maximum) and at 18 ka (LGM II). The dates suggest a significantly earlier maximum advance on the southern side of the range. Reconstructing the geometry of four glaciers in the Sucha Woda, Pańszczyca, Mlynicka and Velicka valleys allowed determining their equilibrium-line altitudes (ELAs) at 1460, 1460, 1650 and 1700 m asl, respectively. Based on a positive degree-day model, the mass balance and climatic parameter anomaly (temperature and precipitation) has been constrained for LGM I advance. Modeling results indicate slightly different conditions between northern and southern slopes. The N–S ELA gradient finds confirmation in slightly higher temperature (at least 1 °C) or lower precipitation (15%) on the south-facing glaciers during LGM I. The precipitation distribution over the High Tatra Mountains indicates potentially different LGM atmospheric circulation than at the present day, with reduced northwesterly inflow and increased southerly and westerly inflows of moist air masses.     

Late Pleistocene equilibrium-line reconstructions in the northern Peruvian Andes

Equilibrium-line-altitude (ELA) reconstructions using the toe-to-headwall-altitude ratio method for paleoglaciers in the Cordilleras Blanca and Oriental, northern Peruvian Andes (7–10°S; 77°20'–77°35'W), indicate that ELAs during the last glacial maximum (LGM; marine isotope stage 2) were c . 4300 m in the Cordillera Blanca, c . 3900–3600 m on the west side of the Cordillera Oriental, and c . 3200 m on the east (Amazon Basin) side of the Cordillera Oriental. Comparison with estimated modern ELAs and glaciation thresholds indicate that ELA depression ranged from c . 700 m in the Cordillera Blanca to c . 1200 m on the east side of the Cordillera Oriental. This augments data from many mountain ranges in middle- and low-latitude regions that indicate that ELAs during the LGM were depressed by c . 1000 m. Published palynological evidence for drier conditions during the LGM in the tropical Andes suggests that ELA depression of this amount involved a temperature reduction (> 5–6°C) that greatly exceeded the tropical sea-surface temperature depression estimates of CLIMAP (< 2°C). The west to east increase in ELA depression during the LGM indicates that the steep modern precipitation gradients may have been even steeper during the LGM.     

Late Quaternary glaciation and equilibrium-line altitudes of the Mayan Ice Cap, Guatemala, Central America

The Sierra los Cuchumatanes (3837 m), Guatemala, supported a plateau ice cap and valley glaciers around Montaña San Juan (3784 m) that totaled ∼ 43 km² in area during the last local glacial maximum. Former ice limits are defined by sharp-crested lateral and terminal moraines that extend to elevations of ∼ 3450 m along the ice cap margin, and to ca. 3000-3300 m for the valley glaciers. Equilibrium-line altitudes (ELAs) estimated using the area-altitude balance ratio method for the maximum late Quaternary glaciation reached as low as 3470 m for the valley glaciers and 3670 m for the Mayan Ice Cap. Relative to the modern altitude of the 0°C isotherm of ∼ 4840 m, we determined ELA depressions of 1110-1436 m. If interpreted in terms of a depression of the freezing level during maximal glaciation along the modern lapse rate of − 5.3°C km^− 1, this ΔELA indicates tropical highland cooling of ∼ 5.9 to 7.6 ± 1.2°C. Our data support greater glacial highland cooling than at sea level, implying a high tropical sensitivity to global climate changes. The large magnitude of ELA depression in Guatemala may have been partially forced by enhanced wetness associated with southward excursions of the boreal winter polar air mass.     

A latest Pleistocene and Holocene glacial history and paleoclimate reconstruction at Three Sisters and Broken Top Volcanoes, Oregon, U.S.A.

At least three sets of moraines mark distinct glacial stands since the last glacial maximum (LGM) in the Three Sisters region of the Oregon Cascade Range. The oldest stand predates 8.1 ka (defined here as post-LGM), followed by a second between ∼ 2 and 8 ka (Neoglacial) and a third from the Little Ice Age (LIA) advance of the last 300 years. The post-LGM equilibrium line altitudes were 260 ± 100 m lower than that of modern glaciers, requiring 23 ± 9% increased winter snowfall and 1.4 ± 0.5°C cooler summer temperatures than at present. The LIA advance had equilibrium line altitudes 110 ± 40 m lower than at present, implying 10 ± 4% greater winter snowfall and 0.6 ± 0.2°C cooler summer temperatures.     

Effect of sea-level lowering on ELA depression during the LGM

Decreases in equilibrium-line altitudes (ELAs) varied geographically during the last glacial maximum (LGM), with a mid-range value of ~ 900 m commonly deduced from altitude ratio and accumulation–area ratio calculations. Sea level, however, was 120 m lower during the LGM, so the ELA lowering relative to sea level would only be 780 m for a 900-m absolute lowering. With a lapse rate of 0.006°C m⁻¹, this implies a 4.7°C lowering of global temperature. It has been argued that this correction for sea-level change is unnecessary, but the logic on which this is based requires adiabatic compression to apply over much longer time scales than is typically invoked. We find that the correction is necessary. In addition, geometric changes in the atmosphere during the LGM, pointed out by Osmaston (2006), could lead to 0.4°C decrease in the average temperature of the troposphere. Additionally, orographic effects could significantly change the snow distribution on mountain masses near sea level.     

Late Pleistocene glaciers in Darhad Basin, northern Mongolia

Late Pleistocene glaciers around Darhad Basin advanced to near their maximum positions at least three times, twice during the Zyrianka glaciation (at ∼ 17-19 ka and ∼ 35-53 ka), and at least once earlier. The Zyrianka glaciers were smaller than their predecessors, but the equilibrium-line altitude (ELA) difference was < 75 m. End moraines of the Zyrianka glaciers were ∼ 1600 m asl; ELAs were 2100-2400 m asl. ¹⁴C and luminescence dating of lake sediments confirm the existence of paleolake highstands in Darhad Basin before ∼ 35 ka. Geologic evidence and ¹⁰Be cosmic-ray exposure dating of drift suggests that at ∼ 17-19 ka the basin was filled at least briefly by a glacier-dammed lake ∼ 140 m deep. However, lake sediments from that time have not yet been recognized in the region. A shallower paleolake briefly occupied the basin at ∼ 11 ka, but between ∼ 11 and 17 ka and after ∼ 10 ka the basin was probably largely dry. The timing of maximum glacier advances in Darhad appears to be approximately synchronous across northern Mongolia, but different from Siberia and western Central Asia, supporting the inference that paleoclimate in Central Asia differed among regions.     

Last Glacial Maximum equilibrium-line altitude trends and precipitation patterns in the Sangre de Cristo Mountains, southern Colorado, USA

Precipitation patterns during the Last Glacial Maximum (LGM) in the Rocky Mountains varied due to the influence of the continental ice sheets and pluvial lakes. However, no constraints have been placed on potential changes of southeasterly Gulf of Mexico-derived moisture that today contributes considerable precipitation to the easternmost ranges of the southern and middle Rocky Mountains. The Sangre de Cristo Mountains of southern Colorado are ideally situated to assess the relative importance of westerly and southeasterly-derived moisture during the LGM. Based on reconstructions of 30 palaeoglaciers in the Sangre de Cristo Mountains, we find that LGM equilibrium-line altitudes (ELAs) on the east side of the range were systematically 100–200 m lower than ELAs on the west side. The observed ELA pattern is strikingly similar to modern precipitation patterns in the study area, suggesting that southeasterly-derived precipitation had a significant influence on the mass balances of LGM glaciers.     

Spatial patterns of Holocene glacier advance and retreat in Central Asia

Glaciers in the southern Himalayas advanced in the early Holocene despite an increase in incoming summer solar insolation at the top of the atmosphere. These glacier advances are in contrast to the smaller alpine glaciers in the western and northern regions of Central Asia. Two different glacier mass-balance models are used to reconcile this Holocene glacier history with climate by quantifying the change in equilibrium-line altitudes (ELA) for simulated changes in Holocene climate. Both ELA models clearly show that the lowering of ELAs in the southern Himalayas is largely due to a decrease in summer temperatures, and that an increase in monsoonal precipitation accounts for less than 30% of the total ELA changes. The decrease in summer temperatures is a dynamic response to the changes in solar insolation, resulting in both a decrease in incoming shortwave radiation at the surface due to an increase in cloudiness and an increase in evaporative cooling. In the western and northern zones of Central Asia, both ELA models show a rise in ELAs in response to a general increase in summer temperatures. This increase in temperatures in the more northern regions is a direct radiative response to the increase in summer solar insolation.     

Late Quaternary glacial chronology on Nevado Illimani, Bolivia, and the implications for paleoclimatic reconstructions across the Andes

¹⁰Be terrestrial cosmogenic nuclide surface exposure ages from moraines on Nevado Illimani, Cordillera Real, Bolivia suggest that glaciers retreated from moraines during the periods 15.5-13.0 ka, 10.0-8.5 ka, and 3.5-2.0 ka. Late glacial moraines at Illimani are associated with an ELA depression of 400-600 m, which is consistent with other local reconstructions of late glacial ELAs in the Eastern Cordillera of the central Andes. A comparison of late glacial ELAs between the Eastern Cordillera and Western Cordillera indicates a marked change toward flattening of the east-to-west regional ELA gradient. This flattening is consistent with increased precipitation from the Pacific during the late glacial period.     

Timing of glaciation and last glacial maximum paleoclimate estimates from the Fish Lake Plateau, Utah

The High Plateaus of Utah include seven separate mountain ranges that supported glaciers during the Pleistocene. The Fish Lake Plateau, located on the eastern edge of the High Plateaus, preserves evidence of at least two glacial advances. Four cosmogenic ³He exposure ages of boulders in an older moraine range from 79 to 159 ka with a mean age of 129 ± 39 ka and oldest ages of 152 ± 3 and 159 ± 5 ka. These ages suggest deposition during the type Bull Lake glaciation and Marine Oxygen Isotope Stage (MIS) 6. Twenty boulder exposure ages from four different younger moraines indicate a local last glacial maximum (LGM) of ~ 21.1 ka, coincident with the type Pinedale glaciation and MIS 2. Reconstructed Pinedale-age glaciers from the Fish Lake Plateau have equilibrium-line altitudes ranging from 2950 to 3190 m. LGM summer temperature depressions for the Fish Lake Plateau range from −10.7 to −8.2°C, assuming no change in precipitation. Comparison of the Fish Lake summer temperature depressions to a regional dataset suggests that the Fish Lake Plateau may have had a slight increase (~ 1.5× modern) in precipitation during the LGM. A series of submerged ridges in Fish Lake were identified during a bathymetric survey and are likely Bull Lake age moraines.     

Sensitivities of the equilibrium line altitude to temperature and precipitation changes along the Andes

Equilibrium line altitudes (ELAs) of alpine glaciers are sensitive indicators of climate change and have been commonly used to reconstruct paleoclimates at different temporal and spatial scales. However, accurate interpretations of ELA fluctuations rely on a quantitative understanding of the sensitivity of ELAs to changes in climate. We applied a full surface energy- and mass-balance model to quantify ELA sensitivity to temperature and precipitation changes across the range of climate conditions found in the Andes. Model results show that ELA response has a strong spatial variability across the glaciated regions of South America. This spatial variability correlates with the distribution of the present-day mean climate conditions observed along the Andes. We find that ELAs respond linearly to changes in temperature, with the magnitude of the response being prescribed by the local lapse rates. ELA sensitivities to precipitation changes are nearly linear and are inversely correlated with the emissivity of the atmosphere. Temperature sensitivities are greatest in the inner tropics; precipitation becomes more important in the subtropics and northernmost mid-latitudes. These results can be considered an important step towards developing a framework for understanding past episodes of glacial fluctuations and ultimately for predicting glacier response to future climate changes.     

Last glacial maximum climate inferences from cosmogenic dating and glacier modeling of the western Uinta ice field, Uinta Mountains, Utah

During the last glacial maximum (LGM), the western Uinta Mountains of northeastern Utah were occupied by the Western Uinta Ice Field. Cosmogenic ¹⁰Be surface-exposure ages from the terminal moraine in the North Fork Provo Valley and paired ²⁶Al and ¹⁰Be ages from striated bedrock at Bald Mountain Pass set limits on the timing of the local LGM. Moraine boulder ages suggest that ice reached its maximum extent by 17.4 ± 0.5 ka (± 2σ). ¹⁰Be and ²⁶Al measurements on striated bedrock from Bald Mountain Pass, situated near the former center of the ice field, yield a mean ²⁶Al/¹⁰Be ratio of 5.7 ± 0.8 and a mean exposure age of 14.0 ± 0.5 ka, which places a minimum-limiting age on when the ice field melted completely. We also applied a mass/energy-balance and ice-flow model to investigate the LGM climate of the western Uinta Mountains. Results suggest that temperatures were likely 5 to 7°C cooler than present and precipitation was 2 to 3.5 times greater than modern, and the western-most glaciers in the range generally received more precipitation when expanding to their maximum extent than glaciers farther east. This scenario is consistent with the hypothesis that precipitation in the western Uintas was enhanced by pluvial Lake Bonneville during the last glaciation.     

Should Quaternary sea-level changes be used to correct glacier ELAs, vegetation belt altitudes and sea level temperatures for inferring climate changes?

Changes in the altitudes of glacier snowlines (ELAs) and the altitudes of montane vegetation belts (VBAs) measure Quaternary climatic change. An accepted ‘correction’ to such changes by deducting the amount of contemporary sea level fall is wrong, since the air displaced by the ice sheets approximately fills the space left by the falling sea level and so there is no overall downward movement of the troposphere. This also causes a reduced cooling at the lowered sea level relative to that at the former inter-glacial sea level, about 1°C at the Las Glacial Maximum, which reduces the discrepancies previously noted by others between terrestrial and marine estimates of sea-level cooling. The change in temperature is indicated by the product of the ELA or VBA lowering and the environmental lapse rate (ELR). Prior estimates of ΔELA (−900 ± 135 m) and ELR (−6° ± 0.1°C km⁻¹) would indicate a cooling of −5.4°C at interglacial sea level and −4.4°C at glacial sea level, although glacial-period ELRs are not known reliably. Established ELA corrections for local epeirogenic uplift or subsidence are appropriate.     

长白山现代理论雪线和古雪线高度

根据冰川地貌和地形特征、岩性、冰川沉积物的风化程度以及OSL测年结果,认为长白山地区发育两期冰川作用,即末次冰盛期和晚冰期,测年结果分别为20.0±2.1ka和11.3±1.2ka。根据平衡线（ELA）处6～8月多年平均气温（T）和年降水量（P）的关系,计算长白山现代理论雪线高度为3380±100m。通过积累区面积比率AAR（accumulation-arearatio）、冰川末端到山顶高度TSAM（the terminal to summit altitudinal）,冰川末端至分水岭平均高度Hofer（the terminal to average elevation of the catchment area）、末端至冰斗后壁比率THAR（toe-to headwall altitude ratios）、冰斗底部高程CF（cirque-floor altitudes method）、侧碛堤最大高度法MELM（maximum elevation of lateral moraines）等方法计算该区末次冰盛期雪线高度为2250～2383m,平均值2320±20m。考虑到末次冰盛期后地壳上升20m,当时雪线的实际高度为2300±20m,冰盛期的雪线降低值为10⁸0±100m。晚冰期北坡和西坡的雪线高度分别为2490m和2440m,平均值2465m,考虑新构造运动后的雪线实际高度2454m,降低值926±100m。长白山新构造运动（LGM上升约20m,晚冰期上升约11m）在末次冰盛期以来对冰川发育的影响不明显。     

中国现代冰川平衡线分布特征与末次冰期平衡线下降值研究

青藏高原及周边山地拥有地球上最高大且最广阔的高山高原,是除两极外最大的现代冰川作用中心,这也使得中国成为中低纬度地区现代冰川最发育的国家之一. 现代冰川平衡线分布具有纬度地带性特征,在青藏高原上还呈不对称的环状. 根据相关研究资料估算,中国末次冰期最盛期时的冰川面积约为50×10⁴ km²,是现代的8.4倍. 基于平衡线处年降水量和夏季平均气温（6-8月）之间的相关关系重建的中国西部（105° E以西）末次冰期最盛期时的平衡线分布图与现代的相似. 在青藏高原内部与西北部,平衡线下降值在500 m以内,小的仅为200~300 m;在青藏高原东南边缘下降值约800 m,最大可达1 000~1 200 m. 天山与阿尔泰山平衡线下降值均在500 m左右. 中国东部（105° E以东）没有发育现代冰川,仅有数处中高山地,如贺兰山、太白山、长白山与台湾山地保存有确切的末次冰期冰川地形,末次冰期最盛期时的平衡线下降800~900 m,大于青藏高原、天山与阿尔泰山地区的下降值. 根据中国东部末次冰期的平衡线分布图以及相关的古气候与古环境研究资料,海拔2 000 m以下的中低山地在第四纪期间任何一次冰川作用中都不具备冰川发育所需的地势条件.     

Progressive glacial retreat in the Southern Altiplano (Uturuncu volcano, 22°S) between 65 and 14 ka constrained by cosmogenic He dating

This work presents the first reconstruction of late Pleistocene glacier fluctuations on Uturuncu volcano, in the Southern Tropical Andes. Cosmogenic ³He dating of glacial landforms provides constraints on ancient glacier position between 65 and 14 ka. Despite important scatter in the exposure ages on the oldest moraines, probably resulting from pre-exposure, these ³He data constrain the timing of the moraine deposits and subsequent glacier recessions: the Uturuncu glacier may have reached its maximum extent much before the global LGM, maybe as early as 65 ka, with an equilibrium line altitude (ELA) at 5280 m. Then, the glacier remained close to its maximum position, with a main stillstand identified around 40 ka, and another one between 35 and 17 ka, followed by a limited recession at 17 ka. Then, another glacial stillstand is identified upstream during the late glacial period, probably between 16 and 14 ka, with an ELA standing at 5350 m. This stillstand is synchronous with the paleolake Tauca highstand. This result indicates that this regionally wet and cold episode, during the Heinrich 1 event, also impacted the Southern Altiplano. The ELA rose above 5450 m after 14 ka, synchronously with the Bolling–Allerod.     

Late pleistocene glacial equilibrium-line altitudes in the Colorado Front Range: A comparison of methods

Six methods for approximating late Pleistocene (Pinedale) equilibrium-line altitudes (ELAs) are compared for rapidity of data collection and error (RMSE) from first-order trend surfaces, using the Colorado Front Range. Trend surfaces computed from rapidly applied techniques, such as glaciation threshold, median altitude of small reconstructed glaciers, and altitude of lowest cirque floors have relatively high RMSEs (97–186 m) because they are subjectively derived and are based on small glaciers sensitive to microclimatic variability. Surfaces computed for accumulation-area ratios (AARs) and toe-to-headwall altitude ratios (THARs) of large reconstructed glaciers show that an AAR of 0.65 and a THAR of 0.40 have the lowest RMSEs (about 80 m) and provide the same mean ELA estimate (about 3160 m) as that of the more subjectively derived maximum altitudes of Pinedale lateral moraines (RMSE = 149 m). Second-order trend surfaces demonstrate low ELAs in the latitudinal center of the Front Range, perhaps due to higher winter accumulation there. The mountains do not presently reach the ELA for large glaciers, and small Front Range cirque glaciers are not comparable to small glaciers existing during Pinedale time. Therefore, Pleistocene ELA depression and consequent temperature depression cannot reliably be ascertained from the calculated ELA surfaces.     

A 40,000-yr record of environmental change from Burial Lake in Northwest Alaska

Burial Lake in northwest Alaska records changes in water level and regional vegetation since ∼ 39,000 cal yr BP based on terrestrial macrofossil AMS radiocarbon dates. A sedimentary unconformity is dated between 34,800 and 23,200 cal yr BP. During all or some of this period there was a hiatus in deposition indicating a major drop in lake level and deflation of lacustrine sediments. MIS 3 vegetation was herb-shrub tundra; more xeric graminoid-herb tundra developed after 23,200 cal yr BP. The tundra gradually became more mesic after 17,000 cal yr BP. Expansions of Salix then Betula, at 15,000 and 14,000 cal yr BP, respectively, are coincident with a major rise in lake level marked by increasing fine-grained sediment and higher organic matter content. Several sites in the region display disrupted sedimentation and probable hiatuses during the last glacial maximum (LGM); together regional data indicate an arid interval prior to and during the LGM and continued low moisture levels until ∼ 15,000 cal yr BP. AMS ¹⁴C dates from Burial Lake are approximately synchronous with AMS ¹⁴C dates reported for the Betula expansion at nearby sites and sites across northern Alaska, but 1000-2000 yr younger than bulk-sediment dates.     

Chemical weathering, river geochemistry and atmospheric carbon fluxes from volcanic and ultramafic regions on Luzon Island, the Philippines

We investigated rates of chemical weathering of volcanic and ophiolitic rocks on Luzon Island, the Philippines. Luzon has a tropical climate and is volcanically and tectonically very active, all factors that should enhance chemical weathering. Seventy-five rivers and streams (10 draining ophiolites, 65 draining volcanic bedrock) and two volcanic hot springs were sampled and analyzed for major elements, alkalinity and ⁸⁷Sr/⁸⁶Sr. Cationic fluxes from the volcanic basins are dominated by Ca²⁺ and Mg²⁺ and dissolved silica concentrations are high (500-1900 μM). Silica concentrations in streams draining ophiolites are lower (400-900 μM), and the cationic charge is mostly Mg²⁺. The areally weighted average CO₂ export flux from our study area is 3.89 ± 0.21 × 10⁶ mol/km²/yr, or 5.99 ± 0.64 × 10⁶ mol/km²/yr from ophiolites and 3.58 ± 0.23 × 10⁶ mol/km²/yr from volcanic areas (uncertainty given as ±1 standard error, s.e.). This is ∼6-10 times higher than the current best estimate of areally averaged global CO₂ export by basalt chemical weathering and ∼2-3 times higher than the current best estimate of CO₂ export by basalt chemical weathering in the tropics. Extrapolating our findings to all tropical arcs, we estimate that around one tenth of all atmospheric carbon exported via silicate weathering to the oceans annually is processed in these environments, which amount to ∼1% of the global exorheic drainage area. Chemical weathering of volcanic terranes in the tropics appears to make a disproportionately large impact on the long-term carbon cycle.     

O/O and D/H in goethite from a North American Oxisol of the Early Eocene climatic optimum

An Early Eocene Oxisol in the Ione Formation of California formed in a coastal continental weathering environment at a paleolatitude of ∼38°N. The dominant minerals in the Oxisol are goethite, quartz, and kaolinite. Material balance calculations were applied to new measurements of chemical composition, D/H, and ¹⁸O/¹⁶O ratios of Oxisol samples to determine the δD (−150 ± 3‰) and δ¹⁸O (−2.4 ± 0.3‰) values of the goethite (α-FeOOH). These data, in combination with the global meteoric water line (MWL), yielded an isotopic temperature of 21(±4) °C. The nominal value of 21 °C contrasts with the modern mean annual temperature (MAT) of 16 °C in that area. The warmer temperature is consistent with formation of the goethite during the Early Eocene climatic optimum. The isotopic composition of the goethite and a temperature of 21 °C imply ancient water with a δD value of −61(±4)‰ and a δ¹⁸O value of −8.9(±0.5)‰. This Early Eocene δ¹⁸O (or δD) value is more negative than values in the range of isotopic scatter observed for modern global precipitation at sites with a MAT of 21 °C.At times of warm global climates, the location of a near-surface atmospheric isotherm would generally shift relative to its location under modern climatic conditions. A simple Rayleigh-type condensation model indicates that, if one “follows the isotherm”, the associated scatter in δD and δ¹⁸O of precipitation in very warm global climates should shift (for a given isotherm) to more negative values that may be detectable in proxy records. The isotopic results from the goethite of the Early Eocene Oxisol appear to add to evidence in support of this idea.