Northern Hemisphere Ice Sheets Research Articles

Orbital- and millennial-scale Asian winter monsoon variability across the Pliocene–Pleistocene glacial intensification

Intensification of northern hemisphere glaciation (iNHG), ~2.7 million years ago (Ma), led to establishment of the Pleistocene to present-day bipolar icehouse state. Here we document evolution of orbital- and millennial-scale Asian winter monsoon (AWM) variability across the iNHG using a palaeomagnetically dated centennial-resolution grain size record between 3.6 and 1.9 Ma from a previously undescribed loess-palaeosol/red clay section on the central Chinese Loess Plateau. We find that the late Pliocene–early Pleistocene AWM was characterized by combined 41-kyr and ~100-kyr cycles, in response to ice volume and atmospheric CO2 forcing. Northern hemisphere ice sheet expansion, which was accompanied by an atmospheric CO2 concentration decline, substantially increased glacial AWM intensity and its orbitally oscillating amplitudes across the iNHG. Superposed on orbital variability, we find that millennial AWM intensity fluctuations persisted during both the warmer (higher-CO2) late Pliocene and colder (lower-CO2) early Pleistocene, in response to both external astronomical forcing and internal climate dynamics.

Glacial isostatic adjustment reduces past and future Arctic subsea permafrost

Sea-level rise submerges terrestrial permafrost in the Arctic, turning it into subsea permafrost. Subsea permafrost underlies ~ 1.8 million km2 of Arctic continental shelf, with thicknesses in places exceeding 700 m. Sea-level variations over glacial-interglacial cycles control subsea permafrost distribution and thickness, yet no permafrost model has accounted for glacial isostatic adjustment (GIA), which deviates local sea level from the global mean due to changes in ice and ocean loading. Here we incorporate GIA into a pan-Arctic model of subsea permafrost over the last 400,000 years. Including GIA significantly reduces present-day subsea permafrost thickness, chiefly because of hydro-isostatic effects as well as deformation related to Northern Hemisphere ice sheets. Additionally, we extend the simulation 1000 years into the future for emissions scenarios outlined in the Intergovernmental Panel on Climate Change’s sixth assessment report. We find that subsea permafrost is preserved under a low emissions scenario but mostly disappears under a high emissions scenario.

Glacial inception through rapid ice area increase driven by albedo and vegetation feedbacks

Abstract. We present transient simulations of the last glacial inception using the Earth system model CLIMBER-X with dynamic vegetation, interactive ice sheets, and visco-elastic solid Earth responses. The simulations are initialized at the middle of the Eemian interglacial (125 kiloyears before present, ka) and run until 100 ka, driven by prescribed changes in Earth's orbital parameters and greenhouse gas concentrations from ice core data. CLIMBER-X simulates a rapid increase in Northern Hemisphere ice sheet area through MIS5d, with ice sheets expanding over northern North America and Scandinavia, in broad agreement with proxy reconstructions. While most of the increase in ice sheet area occurs over a relatively short period between 119 and 117 ka, the larger part of the increase in ice volume occurs afterwards with an almost constant ice sheet extent. We show that the vegetation feedback plays a fundamental role in controlling the ice sheet expansion during the last glacial inception. In particular, with prescribed present-day vegetation the model simulates a global sea level drop of only ∼ 20 m, compared with the ∼ 35 m decrease in sea level with dynamic vegetation response. The ice sheet and carbon cycle feedbacks play only a minor role during the ice sheet expansion phase prior to ∼ 115 ka but are important in limiting the deglaciation during the following phase characterized by increasing summer insolation. The model results are sensitive to climate model biases and to the parameterization of snow albedo, while they show only a weak dependence on changes in the ice sheet model resolution and the acceleration factor used to speed up the climate component. Overall, our simulations confirm and refine previous results showing that climate–vegetation–cryosphere feedbacks play a fundamental role in the transition from interglacial to glacial states characterizing Quaternary glacial cycles.

Surface mass balance and climate of the Last Glacial Maximum Northern Hemisphere ice sheets: simulations with CESM2.1

Abstract. The Last Glacial Maximum (LGM, from ∼26 to 20 ka BP) was the most recent period with large ice sheets in Eurasia and North America. At that time, global temperatures were 5–7 ∘C lower than today, and sea level ∼125 m lower. LGM simulations are useful to understand earth system dynamics, including climate–ice sheet interactions, and to evaluate and improve the models representing those dynamics. Here, we present two simulations of the Northern Hemisphere ice sheet climate and surface mass balance (SMB) with the Community Earth System Model v2.1 (CESM2.1) using the Community Atmosphere Model v5 (CAM5) with prescribed ice sheets for two time periods that bracket the LGM period: 26 and 21 ka BP. CESM2.1 includes an explicit simulation of snow/firn compaction, albedo, refreezing, and direct coupling of the ice sheet surface energy fluxes with the atmosphere. The simulated mean snow accumulation is lowest for the Greenland and Barents–Kara Sea ice sheets (GrIS, BKIS) and highest for British and Irish (BIIS) and Icelandic (IcIS) ice sheets. Melt rates are negligible for the dry BKIS and GrIS, and relatively large for the BIIS, North American ice sheet complex (NAISC; i.e. Laurentide, Cordilleran, and Innuitian), Scandinavian ice sheet (SIS), and IcIS, and are reduced by almost a third in the colder (lower temperature) 26 ka BP climate compared with 21 ka BP. The SMB is positive for the GrIS, BKIS, SIS, and IcIS during the LGM (26 and 21 ka BP) and negative for the NAISC and BIIS. Relatively wide ablation areas are simulated along the southern (terrestrial), Pacific and Atlantic margins of the NAISC, across the majority of the BIIS, and along the terrestrial southern margin of the SIS. The integrated SMB substantially increases for the NAISC and BIIS in the 26 ka BP climate, but it does not reverse the negative sign. Summer incoming surface solar radiation is largest over the high interior of the NAISC and GrIS, and minimum over the BIIS and southern margin of NAISC. Summer net radiation is maximum over the ablation areas and minimum where the albedo is highest, namely in the interior of the GrIS, northern NAISC, and all of the BKIS. Summer sensible and latent heat fluxes are highest over the ablation areas, positively contributing to melt energy. Refreezing is largest along the equilibrium line altitude for all ice sheets and prevents 40 %–50 % of meltwater entering the ocean. The large simulated melt for the NAISC suggests potential biases in the climate simulation, ice sheet reconstruction, and/or highly non-equilibrated climate and ice sheet at the LGM time.

Equilibrium line altitudes of alpine glaciers in Alaska suggest Last Glacial Maximum summer temperature was 2–5 °C lower than during the pre-industrial

Abstract. The lack of continental ice sheets in Alaska during the Last Glacial Maximum (LGM; 26–19 ka) has long been attributed to extensive aridity in the western Arctic. More recently, climate model outputs, a few isolated paleoclimate studies, and global paleoclimate synthesis products show mild summer temperature depressions in Alaska compared to much of the high northern latitudes. This suggests the importance of limited summer temperature depressions in controlling the relatively limited glacier growth during the LGM. To explore this further, we present a new statewide map of LGM alpine glacier equilibrium line altitudes (ELAs), LGM ΔELAs (LGM ELA anomalies relative to the Little Ice Age, LIA), and ΔELA-based estimates of temperature depressions across Alaska to assess paleoclimate conditions. We reconstructed paleoglacier surfaces in ArcGIS to calculate ELAs using an accumulation area ratio (AAR) of 0.58 and an area–altitude balance ratio (AABR) of 1.56. We calculated LGM ELAs (n= 480) in glaciated massifs in the state, excluding areas in southern Alaska that were covered by the Cordilleran Ice Sheet. The data show a trend of increasing ELAs from the southwest to the northeast during both the LGM and the LIA, indicating a consistent southern Bering Sea and northernmost Pacific Ocean precipitation source. Our LGM–LIA ΔELAs from the Alaska Range, supported with limited LGM–LIA ΔELAs from the Brooks Range and the Kigluaik Mountains, average to −355 ± 176 m. This value is much greater than the global LGM average of ca. −1000 m. Using a range of atmospheric lapse rates, LGM–LIA ΔELAs in Alaska translate to summer cooling of < 2–5 ∘C. Our results are consistent with a growing number of local climate proxy reconstructions and global data assimilation syntheses that indicate mild summer temperature across Beringia during the LGM. Limited LGM summer temperature depressions could be explained by the influence of Northern Hemisphere ice sheets on atmospheric circulation.

Variability of East Asian Summer Monsoon precipitation during the Early Pleistocene based on pollen analysis from the western Loess Plateau, China

Understanding Early Pleistocene variability in precipitation associated with the East Asian Summer Monsoon (EASM) may help improve understanding of the effects of future global warming on this system. In this study, a quantitative EASM precipitation record with ∼2.7 kyr resolution was reconstructed for the Early Pleistocene (∼2.2–1.7 Ma) based on a pollen record from a loess-paleosol sequence in Lanzhou Basin, on the western Chinese Loess Plateau (CLP). The results show that the average annual precipitation (Pann) value was ∼505 mm, almost ∼1.6 times as much as today (∼310 mm), implying much stronger EASM intensity from ∼2.2–1.7 Ma. The EASM experienced a period of stronger intensity from ∼2.2 to 2.07 Ma, with Pann up to ∼525 mm and vegetation dominated by forest-steppe, as well as a period of weaker intensity between ∼2.07 and 1.99 Ma with lower average Pann (∼481 mm) and steppe dominated vegetation. A gradual intensification followed from ∼1.99 to 1.7 Ma, with increasing Pann (∼503 mm) accompanied by a significant increase in arboreal abundance. Our analysis suggests that the west to east temperature gradients across the Equatorial Pacific region influenced the Early Pleistocene EASM evolution over geological timescales. This finding differs from the traditional view that high-latitude ice sheet fluctuations were the primary driving force of Quaternary EASM. On orbital timescales, the Pann exhibits a dominant ∼41 kyr obliquity periodicity between ∼2.2–1.99 Ma and a ∼18 kyr precession periodicity between ∼1.99–1.7 Ma, indicating that Northern Hemisphere ice sheet and low-latitude solar insolation were the dominant drivers of EASM precipitation cycles in these two periods.

Ice core evidence for an orbital-scale climate transition on the Northwest Tibetan Plateau

The influences on climate in the Northwest Tibetan Plateau (NwTP) have changed on millennial to precessional timescales and have been dependent on the size of the Northern Hemisphere ice sheets, Northern Hemisphere summer insolation, and the migration of the Intertropical Convergence Zone (ITCZ). All these influences control the position and intensity of the westerlies over Central Asia and the Tibetan Plateau and of the Asian Monsoon over South Asia. The top 187.4 m of a 309.7-m ice core (2015GP) drilled on the plateau of the Guliya ice cap contains a 41-kyr climate history of the NwTP, a region where information on past climate and environment is limited. The oxygen isotope ratio (δ18O), and ammonium (NH4+) and dust concentration records from 2015GP show that during the glacial (41–17.5 ka BP) temperature and precipitation in the NwTP were primarily influenced on the precessional timescale by summer insolation, while at millennial resolution the climate was linked to the North Atlantic temperature via the westerlies. During the deglaciation (17.5–12 ka BP) summer insolation remained an important temperature forcing, but the influence of the North Atlantic climate on the NwTP climate weakened as the westerlies shifted northward. The Guliya Holocene δ18O record shows that NwTP climate was no longer in phase with decreasing summer insolation or with North Atlantic climate, perhaps as the moisture source and pathways were more determinative factors in isotopic fractionation than temperature, and/or incoming solar radiation (insolation) forcing was replaced by rising greenhouse gas concentrations as the primary driver of warming in NwTP.

The discovery of the southernmost ultra-high-resolution Holocene paleoclimate sedimentary record in Antarctica

The response of the Antarctic ice sheet to climate warming is the main source of uncertainty regarding future global sea level rise, since little is known about its present and past dynamics. The last deglaciation is the most recent interval of large-scale climate warming, during which the Northern and Southern Hemisphere ice sheets retreated, and sea level rose globally, although at a non-uniform rate. Geologic records from the polar regions are fundamental in determining the factors that caused the major changes in ice sheets during the last deglacial under different boundary conditions. Here, we combine morpho-bathymetric and seismic data with sediment cores and oceanographic measurements to reconstruct the processes that influenced the deposition of the southernmost, most extensive, ultrahigh-resolution record of the Holocene in Edisto Inlet fjord (Ross Sea, Antarctica). We find that post-glacial sedimentation resulted in a layered diatom mud up to 110 m thick that was locally redistributed by bottom currents over confined drifts-moats in the central part of the fjord. After the Holocene climatic optimum, the fjord was not carved by ground ice, and there continued to be internal fjord water circulation associated with Ross Sea circulation. These results support a retreat of coastal glaciers by about 11 kiloyears ago (ka) from the continental shelf of North Victoria Land.

The Coupled Ice Sheet–Earth System Model Bern3D v3.0

Abstract Understanding climate variability from millennial to glacial–interglacial time scales remains challenging due to the complex and nonlinear feedbacks between ice, ocean, sediments, biosphere, and atmosphere. Complex climate models generally struggle to dynamically and comprehensively simulate such long time periods as a result of the large computational costs. Here, we therefore coupled a dynamical ice sheet model to the Bern3D Earth system model of intermediate complexity, which allows for simulating multiple glacial–interglacial cycles. The performance of the model is first validated against modern observations and its response to abrupt perturbations, such as atmospheric CO2 changes and North Atlantic freshwater hosing, is investigated. To further test the fully coupled model, the climate evolution over the entire last glacial cycle is explored in a transient simulation forced by variations in the orbital configuration and greenhouse gases and aerosols. The model simulates global mean surface temperature in fair agreement with reconstructions, exhibiting a gradual cooling trend since the last interglacial that is interrupted by two more rapid cooling events during the early Marine Isotope Stage (MIS) 4 and Last Glacial Maximum (LGM). Simulated Northern Hemispheric ice sheets show pronounced variability on orbital time scales, and ice volume more than doubles from MIS3 to the LGM in good agreement with recent sea level reconstructions. At the LGM, the Atlantic overturning has a strength of about 14 Sv (1 Sv ≡ 106 m3 s−1), which is a reduction by about one-quarter compared to the preindustrial. We thus demonstrate that the new coupled model is able to simulate large-scale aspects of glacial–interglacial cycles.

Estimating Plio‐Pleistocene North African Monsoon Runoff Into the Mediterranean Sea and Temperature Impacts

AbstractSapropels are dark, organic‐rich layers found in Mediterranean sediments that formed during periods of bottom water anoxia. While various mechanisms have been proposed to have caused anoxic conditions, a primary factor is considered to be water column stratification induced by freshwater runoff related to intensified North African monsoon precipitation during precession minima. Monsoon intensification also induced Green Sahara Periods that may have impacted North African hominin dispersal. In this study, we present a novel regression‐based deconvolution of a high‐resolution planktonic foraminiferal oxygen isotope record to estimate the combination of freshwater runoff reaching the eastern Mediterranean and associated surface warming of the water column over the past 5 million years. Sapropels are known to occur in clusters associated with periods of high orbital eccentricity. Our analysis reveals a consistent influence of orbital eccentricity in modulating the North African monsoon, and a possible shift in runoff source area induced by the initiation of Northern Hemisphere ice sheets. Our findings provide important insights into the role of the North African monsoon in shaping Mediterranean environmental changes over the past 5 million years.

Holocene summer temperature record based on branched tetraethers in Northeast China

The temperature evolution during the Holocene shows conflicting patterns and yields the controversy regarding whether the Holocene followed a warming or cooling trend. Resolving this controversy is critical for understanding the underlying mechanisms of climate change and evaluating global warming on a longer timescale. Here we present a quantitative summer temperature reconstruction based on branched glycerol dialkyl glycerol tetraethers analyzed from a sedimentary sequence retrieved from the volcanic lake Kielguo crater lake in Northeast (NE) China. Our record revealed that the summer temperature in NE China followed a rough cooling trend during ~10–2.2 cal ka BP and increased after ~2.2 cal ka BP. The warmest period appeared at ~10–8.8 cal ka BP with ~15.8°C, and the coldest time occurred at ~2.3 cal ka BP with ~13.4°C. Compared with other temperature records from NE China, we conclude that the summer and winter temperature change trends in NE China were probably controlled by the summer and winter insolation respectively, and the ice volume. The mean annual temperature changes in NE China resulted from the superposition of summer cooling and winter warming during the Holocene, and the decrease in ice volume forcing the increase in mean annual temperatures before ~7–6 cal ka BP. The switch in mean annual temperature changes was probably linked with the collapse of Northern Hemisphere ice sheets.

Quantitative reconstruction of precipitation changes in the Mongolian Altai Mountains since 13.7 ka

Holocene climate changes in the Altai Mountains and the surrounding areas are well documented, but the timing of the Holocene precipitation maximum remains controversial. We obtained 17 accelerator mass spectrometry (AMS) 14C dates and 310 fossil pollen assemblages from sediments in Tolbo Lake, which we used to reconstruct the history of vegetation and climate change in the Mongolian Altai Mountains over the past 13,750 years. The results suggest that mean annual precipitation (Pann) was the most significant control on the fossil pollen record at the study site. Based on the fossil pollen record from Tolbo Lake and 469 surface pollen samples within the radius of 1300 km, the reconstructed Pann was lowest during 13.7–11.3 ka (1 ka = 1,000 cal yr BP); relatively low but increasing Pann occurred in the early Holocene (11.3–9.4 ka); and maximum Pann occurred during the middle to late Holocene (after 6.8 ka), possibly because of the high temperatures. A significant climatic and ecological transition occurred in the Altai Mountains at ∼7 ka: steppe and desert steppe vegetation with generally low Pann was present before ∼7 ka, while subsequently there was the expansion of forest vegetation and a marked increase in Pann, possibly in response to the minimum in Northern Hemisphere ice sheets and the increased radiative forcing caused by rising atmospheric greenhouse gases.

Climate Evolution Through the Onset and Intensification of Northern Hemisphere Glaciation

AbstractThe Pliocene Epoch (∼5.3–2.6 million years ago, Ma) was characterized by a warmer than present climate with smaller Northern Hemisphere ice sheets, and offers an example of a climate system in long‐term equilibrium with current or predicted near‐future atmospheric CO2 concentrations (pCO2). A long‐term trend of ice‐sheet expansion led to more pronounced glacial (cold) stages by the end of the Pliocene (∼2.6 Ma), known as the “intensification of Northern Hemisphere Glaciation” (iNHG). We assessed the spatial and temporal variability of ocean temperatures and ice‐volume indicators through the late Pliocene and early Pleistocene (from 3.3 to 2.4 Ma) to determine the character of this climate transition. We identified asynchronous shifts in long‐term means and the pacing and amplitude of shorter‐term climate variability, between regions and between climate proxies. Early changes in Antarctic glaciation and Southern Hemisphere ocean properties occurred even during the mid‐Piacenzian warm period (∼3.264–3.025 Ma) which has been used as an analog for future warming. Increased climate variability subsequently developed alongside signatures of larger Northern Hemisphere ice sheets (iNHG). Yet, some regions of the ocean felt no impact of iNHG, particularly in lower latitudes. Our analysis has demonstrated the complex, non‐uniform and globally asynchronous nature of climate changes associated with the iNHG. Shifting ocean gateways and ocean circulation changes may have pre‐conditioned the later evolution of ice sheets with falling atmospheric pCO2. Further development of high‐resolution, multi‐proxy reconstructions of climate is required so that the full potential of the rich and detailed geological records can be realized.

Moist and warm conditions in Eurasia during the last glacial of the Middle Pleistocene Transition

The end of the Middle Pleistocene Transition (MPT, ~ 800-670 thousand years before present, ka) was characterised by the emergence of large glacial ice-sheets associated with anomalously warm North Atlantic sea surface temperatures enhancing moisture production. Still, the direction and intensity of moisture transport across Eurasia towards potential ice-sheets is poorly constrained. To reconstruct late MPT moisture production and dispersal, we combine records of upper ocean temperature and pollen-based Mediterranean forest cover, a tracer of westerlies and precipitation, from a subtropical drill-core collected off South-West Iberia, with records of East Asia summer monsoon (EASM) strength and West Pacific surface temperatures, and model simulations. Here we show that south-western European winter precipitation and EASM strength reached high levels during the Marine Isotope Stage 18 glacial. This anomalous situation was caused by nearly-continuous moisture supply from both oceans and its transport to higher latitudes through the westerlies, likely fuelling the accelerated expansion of northern hemisphere ice-sheets during the late MPT.

Timing and Consequences of Bering Strait Opening: New Insights From40Ar/39Ar Dating of the Barmur Group (Tjörnes Beds), Northern Iceland

AbstractThe Barmur Group (informally Tjörnes beds) sedimentary succession of northern Iceland is key to reconstructing the opening of the Bering Strait oceanic gateway because these rocks record migration of bivalve molluscs from the Pacific to the Atlantic via the Arctic. However, the timing of the migration event is poorly constrained owing to a lack of reliable absolute ages. To address this problem, we present the first Ar‐Ar radiometric dates from four basaltic lavas that underlie, are intercalated with, and overlie the Barmur Group, and integrate them with existing paleomagnetic records. We show that the Barmur Group has a latest Miocene to early Pliocene age range (c. 6.0–4.4 Ma; C3r–C3n.2n), older than all previous age models. Thus, the Barmur Group does not record the mid‐Piacenzian Warm Period, contra some previous suggestions. Abundant Pacific bivalve molluscs appeared in the Barmur Group during subchrons C3n.4n–C3n.3r at 5.235–4.896 Ma, over 1.3 million years earlier than previously suggested. Appearance of Pacific bivalves in the northern Atlantic occurred shortly after the 5.6–5.4 Ma age previously inferred for first appearance of Arctic bivalves in the Pacific. Thus, our data suggest that first opening of the Bering Strait gateway by the latest Miocene (c. 5.5 Ma) was soon followed by bidirectional trans‐Arctic faunal exchange, and argue against a hypothesized two‐stage faunal exchange process spanning c. 2 million years. Our results also confirm that first opening of the Bering Strait gateway was not directly associated with the growth of large northern hemisphere icesheets, which occurred several million years later.

Optically stimulated luminescence dating and paleoclimatic implications of the Holocene dune sands in the Hunshandake Sandy Land, Northeast China

Aeolian sediments are important geological archives of paleoclimate changes. Due to the fragile ecosystem of desert, sand dune mobility is particularly sensitive to climate change, and is thus important for revealing the stabilization and mobilization processes of dunes during the geological past. In northeastern China, several sandy lands are located in the East Asian monsoonal region, and desert evolution has been used to infer the waxing and waning of monsoonal climate. However, the relationship between dynamic aeolian landscapes and forcing mechanisms is still under debate due to the complex nature of the aeolian dune systems. In this study, we analyse two aeolian sand dune sections in the Hunshandake Sandy Land, northeast China, located at the northern edge of the East Asian summer monsoon, explore regional paleoclimatic changes during the Holocene. We carried out high-resolution optically stimulated luminescence (OSL) dating to develop a well-constrained chronology for the dune sands, as well as analyses of multiple climatic proxies including magnetic susceptibility, particle size, and color parameters, to obtain knowledge about dune mobility and climatic changes. Our results indicate that the early Holocene was characterized by the accumulation of aeolian dune sands, implying an arid climate, driven by an enhanced East Asian winter monsoon. Paleosols were generally developed in the mid-Holocene (∼8 and 4 ka) when the sand dunes were stabilized by vegetation under a relatively warm and humid climate driven by an enhanced East Asian summer monsoon. In the late Holocene, from ∼4 ka to the present, climate returned to an arid and semi-arid phase, marked by aeolian dune sands with some weakly-developed paleosols. Through comparison of the aeolian record with the other East Asian monsoonal records and the sea-surface temperature of the North Pacific, we suggest that the evolution of the East Asian monsoon, which was paced by both the melting and expansion from the Northern Hemisphere ice sheet as well as the changes of solar radiation, played an important role in driving climate change in the Hunshandake Sandy Land during the Holocene.

Northern hemisphere ice sheet expansion intensified Asian aridification and the winter monsoon across the mid-Pleistocene transition

The mid-Pleistocene transition 1.25 to 0.6 million years ago marked a major shift in global climate periodicity from 41,000 to around 100,000 years without a concomitant orbital forcing shift. Here, we investigate Asian climate dynamics associated with two extreme glacial loess coarsening events at the onset and middle of the mid-Pleistocene transition by combining new and existing grain size and magnetic susceptibility records from the Chinese Loess Plateau spanning the last 1.6 million years with general circulation model simulations. We find that the two extreme glacial events reflect exceptionally enhanced Asian aridification and winter monsoon activity. They coincided with notable Northern Hemisphere glacial ice sheet expansion at 1.25 and 0.9 million years ago when the 100,000-year periodicity initiated and intensified, respectively. Our results indicate that these anomalously dry and windy Asian glacials were probably driven by an amplified terrestrial climate response to the coincident Northern Hemisphere ice sheet expansion.

Modelling feedbacks between the Northern Hemisphere ice sheets and climate during the last glacial cycle

Abstract. During the last glacial cycle (LGC), ice sheets covered large parts of Eurasia and North America, which resulted in ∼120 m of sea level change. Ice sheet–climate interactions have considerable influence on temperature and precipitation patterns and therefore need to be included when simulating this time period. Ideally, ice sheet–climate interactions are simulated by a high-resolution Earth system model. While these models are capable of simulating climates at a certain point in time, such as the pre-industrial (PI) or the Last Glacial Maximum (LGM; 21 000 years ago), a full transient glacial cycle is currently computationally unfeasible as it requires a too-large amount of computation time. Nevertheless, ice sheet models require forcing that captures the gradual change in climate over time to calculate the accumulation and melt of ice and its effect on ice sheet extent and volume changes. Here we simulate the LGC using an ice sheet model forced by LGM and PI climates. The gradual change in climate is modelled by transiently interpolating between pre-calculated results from a climate model for the LGM and the PI. To assess the influence of ice sheet–climate interactions, we use two different interpolation methods: the climate matrix method, which includes a temperature–albedo and precipitation–topography feedback, and the glacial index method, which does not. To investigate the sensitivity of the results to the prescribed climate forcing, we use the output of several models that are part of the Paleoclimate Modelling Intercomparison Project Phase III (PMIP3). In these simulations, ice volume is prescribed, and the climate is reconstructed with a general circulation model (GCM). Here we test those models by using their climate to drive an ice sheet model over the LGC. We find that the ice volume differences caused by the climate forcing exceed the differences caused by the interpolation method. Some GCMs produced unrealistic LGM volumes, and only four resulted in reasonable ice sheets, with LGM Northern Hemisphere sea level contribution ranging between 74–113 m with respect to the present day. The glacial index and climate matrix methods result in similar ice volumes at the LGM but yield a different ice evolution with different ice domes during the inception phase of the glacial cycle and different sea level rates during the deglaciation phase. The temperature–albedo feedback is the main cause of differences between the glacial index and climate matrix methods.

Ice Volume and Insolation Forcing of Abrupt Strengthening of East Asian Winter Monsoon During Glacial Inceptions

AbstractIt is generally accepted that the glacial‐interglacial variations of the East Asian winter monsoon (EAWM) are controlled by the volume of the Northern Hemisphere ice sheets (NHIS), but the fact that they exhibit different evolution patterns during glacial inceptions is often overlooked. By generating an independent chronology framework and integrating multi‐proxy records from the loess sections on the central Chinese Loess Plateau, here we show that the rapid intensifying of the EAWM during glacial inceptions reflects millennial perturbations in the large‐scale atmospheric circulation in East Asia in response to insolation‐triggered abrupt North Atlantic cooling. This climate teleconnection between North Atlantic and East Asia is found to be particularly effective only when the NHIS reaches a critical large size. Our integrated multi‐proxy records thus highlight the key role of ice volume in modulating the response of the EAWM to insolation‐triggered North Atlantic cooling during the interglacial‐glacial transitions.

The Mid-Pleistocene Climate Transition

The timing of ice ages over the past ∼2,600 thousand years (kyr) follows pacing by cyclical changes in three aspects of Earth's orbit that influence the solar energy received as a function of latitude and season. Explaining the large magnitude of the climate changes is challenging, particularly so across the period of time from ∼1,250 to 750 ka—the Mid-Pleistocene Transition or MPT. The average repeat time of ice age cycles changed from an earlier 41-kyr rhythm to longer and more intense glaciations at a spacing of about 100 kyr. Explaining this change is very difficult because there was no corresponding change in the orbital pacing that would trigger a change in timing. While the first generation of hypotheses looked largely to changes in the behavior of Northern Hemisphere ice sheets, more recent work integrates ice behavior with new data capturing the evolution of other important aspects of past climate. A full explanation is still lacking, but attention increasingly focuses on the ocean carbon cycle and atmospheric CO2levels as the crucial agents involved in the MPT. ▪ The pattern of climate changes connected to the ice ages of the past few million years changed radically between about 1,250 and 750 thousand years ago, a time known as the Mid-Pleistocene Transition or MPT. ▪ While the glacial cycles were ultimately triggered by cyclical changes in Earth's orbit, the changes across the MPT came from changes in the Earth system itself, most likely in the form of a change in the carbon cycle. ▪ The dramatic change in many essential aspects of climate—ice volume, temperature, rainfall on land, and many others—in the absence of an external change suggests how important feedbacks are to the climate system.