Search

We compile and analyze all available geothermal heat flow measurements collected in and around Greenland into a new database of 419 sites and generate an accompanying spatial map. This database includes 290 sites previously reported by the International Heat Flow Commission (IHFC), for which we now standardize measurement and metadata quality. This database also includes 129 new sites, which have not been previously reported by the IHFC. These new sites consist of 88 offshore measurements and 41 onshore measurements, of which 24 are subglacial. We employ machine learning to synthesize these in situ measurements into a gridded geothermal heat flow model that is consistent across both continental and marine areas in and around Greenland. This model has a native horizontal resolution of 55 km. In comparison to five existing Greenland geothermal heat flow models, our model has the lowest mean geothermal heat flow for Greenland onshore areas. Our modeled heat flow in central North Greenland is highly sensitive to whether the NGRIP (North GReenland Ice core Project) elevated heat flow anomaly is included in the training dataset. Our model's most distinctive spatial feature is pronounced low geothermal heat flow (

It is well known that ice sheet–climate feedbacks are essential for realistically simulating the spatiotemporal evolution of continental ice sheets over glacial–interglacial cycles. However, many of these feedbacks are dependent on the ice sheet thickness, which is poorly constrained by proxy data records. For example, height estimates of the Laurentide Ice Sheet (LIS) topography at the Last Glacial Maximum (LGM; ∼ 21 000 years ago) vary by more than 1 km among different ice sheet reconstructions. In order to better constrain the LIS elevation it is therefore important to understand how the mean climate is influenced by elevation discrepancies of this magnitude. Here we use an atmospheric circulation model coupled to a slab-ocean model to analyze the LGM surface temperature response to a broad range of LIS elevations (from 0 to over 4 km). We find that raising the LIS topography induces a widespread surface warming in the Arctic region, amounting to approximately 1.5 ∘C per km of elevation increase, or about 6.5 ∘C for the highest LIS. The warming is attributed to an increased poleward energy flux by atmospheric stationary waves, amplified by surface albedo and water vapor feedbacks, which account for about two-thirds of the total temperature response. These results suggest a strong feedback between continental-scale ice sheets and the Arctic temperatures that may help constrain LIS elevation estimates for the LGM and explain differences in ice distribution between the LGM and earlier glacial periods.

Coupled climate–ice sheet simulations have been growing in popularity in recent years. Experiments of this type are however challenging as ice sheets evolve over multi-millennial timescales, which is beyond the practical integration limit of most Earth system models. A common method to increase model throughput is to trade resolution for computational efficiency (compromise accuracy for speed). Here we analyze how the resolution of an atmospheric general circulation model (AGCM) influences the simulation quality in a stand-alone ice sheet model. Four identical AGCM simulations of the Last Glacial Maximum (LGM) were run at different horizontal resolutions: T85 (1.4∘), T42 (2.8∘), T31 (3.8∘), and T21 (5.6∘). These simulations were subsequently used as forcing of an ice sheet model. While the T85 climate forcing reproduces the LGM ice sheets to a high accuracy, the intermediate resolution cases (T42 and T31) fail to build the Eurasian ice sheet. The T21 case fails in both Eurasia and North America. Sensitivity experiments using different surface mass balance parameterizations improve the simulations of the Eurasian ice sheet in the T42 case, but the compromise is a substantial ice buildup in Siberia. The T31 and T21 cases do not improve in the same way in Eurasia, though the latter simulates the continent-wide Laurentide ice sheet in North America. The difficulty to reproduce the LGM ice sheets in the T21 case is in broad agreement with previous studies using low-resolution atmospheric models, and is caused by a substantial deterioration of the model climate between the T31 and T21 resolutions. It is speculated that this deficiency may demonstrate a fundamental problem with using low-resolution atmospheric models in these types of experiments.

A recent data campaign in the East Siberian Sea has revealed evidence of grounded and floating ice dynamics in regions of up to 1000 m water depth, and which are attributed to glaciations older than the Last Glacial Maximum (21 kyrs BP). The main hypothesis based on this evidence is that a small ice cap developed over Beringia and expanded over the East Siberian continental margin during some of the Late Pleistocene glaciations. Other similar evidence of ice dynamics that have been previously collected on the shallow continental shelves of the Arctic Ocean have been attributed to the penultimate glaciation, i.e. Marine Isotopes Stage 6 (≈140 kyrs BP). We use an ice sheet model, forced by two previously simulated MIS 6 glacial maximum climates, to carry out a series of sensitivity experiments testing the impact of dynamics and mass-balance related parameters on the geometry of the East Siberian ice cap and ice shelf. Results show that the ice cap developing over Beringia connects to the Eurasian ice sheet in all simulations and that its volume ranges between 6 and 14 m SLE, depending on the climate forcing. This ice cap generates an ice shelf of dimensions comparable with or larger than the present-day Ross ice shelf in West Antarctica. Although the ice shelf extent strongly depends on the ice flux through the grounding line, it is particularly sensitive to the choice of the calving and basal melting parameters. Finally, inhibiting a merging of the Beringia ice cap with the Eurasian ice sheet affects the expansion of the ice shelf only in the simulations where the ice cap fluxes are not large enough to compensate for the fluxes coming from the Eurasian ice sheet.

This study examines the evolution of a continental-scale ice sheet on a triangular representation of North America, with and without the influence of the Cordilleran region. Simulations are conducted using a comprehensive atmospheric general circulation model asynchronously coupled to a three-dimensional thermomechanical ice-sheet model. The atmospheric state is updated for every 2 × 106 km3 increase in ice volume, and the coupled model is integrated to steady state. In the first experiment a flat continent with no background topography is used. The ice sheet evolves fairly zonally symmetric, and the equilibrium state is continent-wide and has the highest point in the center of the continent. This equilibrium ice sheet forces an anticyclonic circulation that results in relatively warmer (cooler) summer surface temperatures in the northwest (southeast), owing to warm (cold) air advection and radiative heating due to reduced cloudiness. The second experiment includes a simplified representation of the Cordilleran region. The ice sheet’s equilibrium state is here structurally different from the flat continent case; the center of mass is strongly shifted to the east and the interior of the continent remains ice free—an outline broadly resembling the geologically determined ice margin in Marine Isotope Stage 4. The limited glaciation in the continental interior is the result of warm summer surface temperatures primarily due to stationary waves and radiative feedbacks.

Coupled climate–ice-sheet simulations have been growing in popularity in recent years. Experiments of this type are however challenging as ice sheets evolve over multi-millennial time scales, which is beyond the practical integration limit for most Earth-system models. A common method to increase model throughput is to trade resolution for computational efficiency (compromises accuracy for speed). Here, we analyze how the resolution of an atmospheric general circulation model (AGCM) influences the simulation quality of a standalone ice-sheet model. Four identical AGCM simulations of the Last Glacial Maximum (LGM) were run at different horizontal resolutions: T85 (1.4°), T42 (2.8°), T31 (3.8°), and T21 (5.6°). These simulations were subsequently used as forcing of an ice-sheet model. While the T85 climate forcing reproduces the LGM ice sheets to a high accuracy, the intermediate resolution cases (T42 and T31) fail to build the Eurasian Ice Sheet. The T21 case fails in both Eurasia and North America. Sensitivity experiments using different surface mass balance parameterizations improve the simulations of the Eurasian ice-sheet in the T42 case, but the compromise is a substantial ice buildup in Siberia. The T31 and T21 cases are not improving in the same way in Eurasia, though the latter simulates the continent-wide Laurentide Ice Sheet in North America. The difficulty to reproduce the LGM ice sheets in the T21 case is in broad agreement with previous studies using low-resolution atmospheric models, and is caused by a deterioration of the atmospheric climate between the T31 and T21 resolutions. It is speculated that this deficiency may demonstrate a fundamental problem using low-resolution atmospheric models in these types of experiments.

The usefulness of sedimentary charcoal records to document centennial to millennial scale trends in aspects of fire regimes (frequency, severity) is widely acknowledged, yet the long-term variability in these regimes is poorly understood. Here, we use a high-resolution, multi-proxy analysis of a lacustrine sequence located in the lowlands of Transylvania (NW Romania), alongside global climate simulations in order to disentangle the drivers of fire regimes in this dry climatic region of Central-Eastern Europe. Periods of greater fire activity and frequency occurred between 10,700 and 7100 cal yr BP (mean Fire Interval ¼ mFI 112 yr), and between 3300 and 700 cal yr BP (mFI 150 yr), whereas intervals of lower fire activity were recorded between 12,000 and 10,700 cal yr BP (mFI 217 yr), 7100 and 3300 cal yr BP (mFI 317 yr), and over last 700 years (no fire events detected).We found good correlations between simulated early summer (June, July) soil moisture content and near-surface air temperature with fire activity, particularly for the early to mid Holocene. A climateefire relationship is further supported by local hydrological changes, i.e., lake level and runoff fluctuations. Fuel limitation, as a result of arid and strongly seasonal climatic conditions, led to low fire activity before 10,700 cal yr BP. However, fires were most frequent during climatically drier phases for the remaining, fuel-sufficient, part of the Holocene. Our results also suggest that the occurrence of more frequent fires in the early Holocene has kept woodlands open, promoted grassland abundance and sustained a more flammable ecosystem (mFI < 150 years) whereas the decline in fire risk under cooler and wetter climate conditions (mFI ¼ 317 years) favoured woodland development. From 3300 cal yr BP, human impacts clearly were partly responsible for changes in fire activity, first increasing fire frequency and severity in periods with fire-favourable climatic conditions (halving the mFI from 300 years to about 150 years), then effectively suppressing fires over the last several centuries. Given the projected future temperature increase and moisture decline and the biomass accumulation due to the agricultural land abandonment in the region,natural fire frequency would be expected to return to

This study examines mutualinteractions between stationary waves and ice sheets using a dry atmospheric primitive-equation model coupled to a three-dimensional thermomechanical ice-sheet model. The emphasis is on how non-linear interactions between thermal and topographical forcing of the stationary waves influence the ice-sheet evolution by changing the ablation. Simulations are conducted in which a small ice cap, on an idealised Northern Hemisphere continent, evolves to an equilibrium continental-scale ice sheet. In the absence of stationary waves, the equilibrium ice sheet arrives at symmetric shape with a zonal equatorward margin. In isolation, the topographically induced stationary waves have essentially no impact on the equilibrium features of the ice sheet. The reason is that the temperature anomalies are located far from the equatorward ice margin. When forcing due to thermal cooling is added to the topographical forcing, thermally induced perturbation winds amplify the topographically induced stationary-wave response, which that serves to increase both the equatorward extent and the volume of the ice sheet. Roughly, a 10% increase in the ice volume is reported here. Hence, the present study suggests that the topographically induced stationary-wave response can be substantially enhanced by the high albedo of ice sheets. Keywords: stationary waves; ice sheets; topographical forcing; thermal forcing (Published: 2 January 2012) Citation: Tellus A 2012, 64 , 11088, DOI: 10.3402/tellusa.v64i0.11088

This study examines the mutual interaction between topographically-forced atmospheric stationary waves and continental-scale ice sheets using a thermomechanical ice-sheet model coupled to a linear ...

A linear two-level atmospheric model is employed to study the influence of ice-sheet topography on atmospheric stationary waves. In particular, the stationary-wave-induced temperature anomaly is considered locally over a single ice-sheet topography, which is computed using the plastic approximation. It is found that stationary waves induce a local cooling which increases linearly with the ice volume for ice sheets of horizontal extents smaller than ∼1400 km. Beyond this horizontal scale, the dependence of stationary-wave-induced cooling on the ice volume becomes gradually weaker. For a certain ice-sheet size, and for small changes of the surface zonal wind, it is further shown that the strength of the local stationary-wave-induced cooling is proportional to the basic state meridional temperature gradient multiplied by the vertical stratification in the atmosphere. These results are of importance for the nature of the feedback between ice sheets and stationary waves, and may also serve as a basis for parameterizing this feedback in ice-sheet model simulations (e.g. through the Pleistocene glacial/interglacial cycles).

Modeling studies show that the massive ice sheet expanding over the North American and Eurasian continents in the last glacial cycle has a large impact on the atmospheric stationary waves and thus yielded a glacial climate distinctly different from the present. However, to what extent the two ice sheets influenced each others growth trajectories remains largely unexplored. In this study we investigate how ice sheets in North America influence the downstream evolution of the Eurasian ice sheet, using a thermomechanical ice-sheet model forced by climate data from snapshot simulations of three distinctly different phases of the last glacial cycle: the Marine Isotope Stages 5b, 4 and 2 (LGM). Our results suggest that changes in the North American paleo-topography may have had a large influence on evolution of the Eurasian ice sheet. In the MIS4 and LGM experiments, the Eurasian ice sheet migrates westward towards the Atlantic sector – largely consistent with geological data and contemporary ice-sheet reconstructions – due to a low wavenumber stationary wave response, which yields a cooling in Europe and a warming in northeastern Siberia. The expansion of the North American ice sheet between MIS4 and LGM amplifies the Siberian warm anomaly, which limits the glaciation there and may therefore help to explain the progressive westward migration of the Eurasian ice sheet over this time period. While the Eurasian ice sheet in the MIS4 and LGM experiments appears to be in equilibrium with the simulated climate conditions, the MIS5b climate forcing is too warm to grow an ice sheet. First-order sensitivity experiments suggest that most of the MIS5b ice sheet was established during preceding colder stages.

International audience; This work focuses on the Late Saalian (140 ka) Eurasian ice sheets' surface mass balance (SMB) sensitivity to changes in sea surface temperatures (SST). An Atmospheric General Circulation Model (AGCM), forced with two preexisting Last Glacial Maximum (LGM, 21 ka) SST reconstructions, is used to compute climate at 140 and 21 ka (reference glaciation). Contrary to the LGM, the ablation almost stopped at 140 ka due to the climatic cooling effect from the large ice sheet topography. Late Saalian SST are simulated using an AGCM coupled with a mixed layer ocean. Compared to the LGM, these 140 ka SST show an inter-hemispheric asymmetry caused by the larger ice-albedo feedback, cooling climate. The resulting Late Saalian ice sheet SMB is smaller due to the extensive simulated sea ice reducing the precipitation. In conclusion, SST are important for the stability and growth of the Late Saalian Eurasian ice sheet.