Johan Liakka

Last Millennium hydro-climate variability in Central–Eastern Europe (Northern Carpathians, Romania)

Proxy-based reconstructions of climate variability over the last millennium provide important insights for understanding current climate change within a long-term context. Past hydrological changes are particularly difficult to reconstruct, yet rainfall patterns and variability are among the most critical environmental variables. Ombrotrophic bogs, entirely dependent on water from precipitation and sensitive to changes in the balance between precipitation and evapotranspiration, are highly suitable for such hydro-climate reconstructions. We present a multi-proxy analysis (testate amoebae, plant macrofossils, stable carbon isotopes in Sphagnum, pollen, spores and macroscopic charcoal) from an ombrotrophic peat profile from the Rodna Mountains (northern Romania) to establish a quantitative record of hydro-climatic changes. We identify five main stages: wet surface mire conditions between AD 800 and 1150 and AD 1800 and 1950, and drying of the mire surface between AD 1300 and 1450, AD 1550 and 1750 and AD 1950 and 2012. Our multi-proxy reconstructions suggest that conditions during the Medieval Climate Anomaly (MCA) period (AD 900–1150) were considerably wetter than today, while during most of the ‘Little Ice Age’ (LIA; AD 1500–1850), they were dry. Mire surface conditions in the Rodna Mountains have dried markedly over the last 40 years mainly as a result of anthropogenic climate change approaching the driest conditions seen over the last 1000 years. There is a marked difference between current hydro-climatic conditions (dry mire) and those of the MCA (wet mire). This implies that for the study region, the MCA cannot provide analogous climatic conditions to the contemporary situation. Our reconstructions are in partial agreement with water table estimates elsewhere in central and eastern Europe but generally contrast with those from NW Europe, especially during LIA. We suggest that these distinctive regional differences result from fluctuations in large-scale atmospheric circulation, which determine the relative influences of continental and oceanic air masses.

Interactions between stationary waves and ice sheets: linear versus nonlinear atmospheric response

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 as well as a fully-nonlinear dry atmospheric primitive equation model. The focus is on how the stationary-wave induced ablation feeds back on the ice sheet. Simulations are conducted in which an embryonal ice mass, on an idealised “North American” continent, evolves to an equilibrium ice sheet. Under the coupling to the linear atmospheric model, the equilibrium ice sheet is primarily controlled by the ratio between the wavelength of the stationary waves and the zonal continental extent. When this ratio is near two, the ice sheet has its center of mass shifted far eastward and its shape is broadly reminiscent of the Laurentide ice sheet at LGM. For wavelengths comparable to the continental extent, however, the ice margin extends far equatorward on the central continent but is displaced poleward near the eastern coast. Remarkably, the coupling to the nonlinear atmospheric model yields equilibrium ice sheets that are virtually identical to the ones obtained in uncoupled simulations, i.e. a symmetric ice sheet with a zonal southern margin. Thus, the degree of linearity of the atmospheric response should control to what extent topographically-forced stationary waves can reorganise the structure of ice sheets. If the stationary-wave response is linear, the present results suggest that spatial reconstructions of past ice sheets can provide some information on the zonal-mean atmospheric circulation that prevailed.

The North American Cordillera—An Impediment to Growing the Continent-Wide Laurentide Ice Sheet

AbstractThis 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 t...

The impact of climate‐vegetation interactions on the onset of the Antarctic ice sheet

A global coupled atmosphere/vegetation model and a dynamic ice sheet model were employed to study the impact of climate-vegetation interactions on the onset of the Antarctic ice sheet during the Eocene-Oligocene transition. We found that the CO2 threshold for Antarctic glaciation is highly sensitive to the prevailing vegetation. In our experiments, the CO2 threshold is less than 280 ppm if the Antarctic vegetation is dominated by forests and between 560 and 1120 ppm for tundra and bare ground conditions. The large impact of vegetation on inception is attributed to the ability of canopies to shade the snow-covered ground, which leads to a weaker snow albedo feedback and higher summer temperatures. However, the overall effect of canopy shading on the Antarctic climate also depends on features like local cloudiness and atmospheric meridional heat transport. Our results suggest that vegetation feedbacks on climate are crucial for the timing of the Antarctic glaciation.

The impact of topographically forced stationary waves on local ice-sheet climate

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 co ...

Effect of changing vegetation on denudation (part 1): Predicted vegetation composition and cover over the last 21 thousand years along the Coastal Cordillera of Chile

Vegetation is crucial for modulating rates of denudation and landscape evolution, as it stabilizes and protects hillslopes and intercepts rainfall. Climate conditions and the atmospheric CO2 concentration, hereafter [CO2], influence the establishment and performance of plants; thus, these factors have a direct influence on vegetation cover. In addition, vegetation dynamics (competition for space, light, nutrients, and water) and stochastic events (mortality and fires) determine the state of vegetation, response times to environmental perturbations and successional development. In spite of this, state-of-the-art reconstructions of past transient vegetation changes have not been accounted for in landscape evolution models. Here, a widely used dynamic vegetation model (LPJ-GUESS) was used to simulate vegetation composition/cover and surface runoff in Chile for the Last Glacial Maximum (LGM), the mid-Holocene (MH) and the present day (PD). In addition, transient vegetation simulations were carried out from the LGM to PD for four sites in the Coastal Cordillera of Chile at a spatial and temporal resolution adequate for coupling with landscape evolution models. A new landform mode was introduced to LPJ-GUESS to enable a better simulation of vegetation dynamics and state at a sub-pixel resolution and to allow for future coupling with landscape evolution models operating at different spatial scales. Using a regionally adapted parameterization, LPJ-GUESS was capable of reproducing PD potential natural vegetation along the strong climatic gradients of Chile, and simulated vegetation cover was also in line with satellite-based observations. Simulated vegetation during the LGM differed markedly from PD conditions. Coastal cold temperate rainforests were displaced northward by about 5 and the tree line and vegetation zones were at lower elevations than PD. Transient vegetation simulations indicate a marked shift in vegetation composition starting with the past glacial warming that coincides with a rise in [CO2]. Vegetation cover between the sites ranged from 13 % (LGM: 8 %) to 81 % (LGM: 73 %) for the northern Pan de Azúcar and southern Nahuelbuta sites, respectively, but did not vary by more than 10 % over the 21 000 year simulation. A sensitivity study suggests that [CO2] is an important driver of vegetation changes and, thereby, potentially Published by Copernicus Publications on behalf of the European Geosciences Union. 830 C. Werner et al.: Effect of changing vegetation and precipitation on denudation landscape evolution. Comparisons with other paleoclimate model drivers highlight the importance of model input on simulated vegetation. In the near future, we will directly couple LPJ-GUESS to a landscape evolution model (see companion paper) to build a fully coupled dynamic-vegetation/landscape evolution model that is forced with paleoclimate data from atmospheric general circulation models.

Arctic warming induced by the Laurentide Ice Sheet topography

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 twothirds 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.

On the limited ice intrusion in Alaska at the LGM: ON THE ICE-FREE LGM ALASKA

The Last Glacial Maximum (LGM) Laurentide Ice Sheet covered most of the North American continent poleward of 40∘N, with the exception of Alaska that remained relatively warm, dry, and largely ice free. Experiments with a global atmospheric circulation model are in broad agreement with proxies: the Alaskan summer temperatures are comparable to the preindustrial, and the annual precipitation is reduced by 30–50%. The warm conditions are attributed to a lowering of the local planetary albedo—due to a decreased cloudiness in response to the cold LGM sea surface temperatures (SSTs) and a stationary anticyclone forced by the ice sheet—that allows more shortwave radiation to reach the surface. Stationary waves are shown to counteract the shortwave cloud feedback by converging less heat over the target region. The LGM SST field also yields an equatorward shifted Pacific stormtrack, which results in drier conditions in Alaska and abundant precipitation at the southern margin of the Laurentide Ice Sheet.