A hybrid GCM paleo ice-sheet model, ANICE2.1 – [email protected]: set up and benchmark experiments

Fully coupled ice-sheet-climate modelling over 10,000 – 100,000-year time scales on high spatial and temporal resolution remains beyond the capability of current computational systems. Hybrid GCM-ice-sheet modelling offers a middle ground, balancing the need to accurately capture both long-term processes, in particular circulation driven changes in precipitation, and processes requiring a high spatial resolution like ablation. Here, we present and evaluate a model set-up that forces the ANICE 3D thermodynamic ice-sheet-shelf model calculating the four large continental ice-sheets (Antarctica, Greenland, North America and Eurasia), with pre-calculated output from two steady-state simulations with the HadCM3 general circulation model (GCM), using a so-called matrix method of coupling both components, where simulations with various levels of pCO2 and ice-sheet configuration are combined to form a time-continuous transient climate forcing consistent with the modelled ice-sheets. We address the difficulties in downscaling low-resolution GCM output to the higher-resolution grid of an ice-sheet model, and account for differences between GCM and ice-sheet model surface topography ranging from interglacial to glacial conditions. Although the approach presented here can be applied to a matrix with any number of GCM snapshots, we limited our experiments to a matrix of only two snapshots. As a benchmark experiment to assess the validity of this model set-up, we perform a simulation of the entire last glacial cycle, from 120 kyr ago to present-day. The simulated eustatic sea-level drop at the Last Glacial maximum (LGM) for the combined Antarctic, Greenland, Eurasian and North-American ice-sheets amounts to 100 m, in line with many other studies. The simulated ice-sheets at LGM agree well with the ICE-5G reconstruction and the more recent DATED-1 reconstruction in terms of total volume and geographical location of the ice sheets. Moreover, modelled benthic oxygen isotope abundance and the relative contributions from global ice volume and deep-water temperature agree well with available data, as do surface temperature histories for the Greenland and Antarctic ice-sheets. This model strategy can be used to create time-continuous ice-sheet distribution and sea-level reconstructions for geological periods up to several millions of years in duration, capturing climate model driven variations in the mass balance of the ice sheet

Antarctic ice sheet and oceanographic response to eccentricity forcing during the early Miocene

Stable isotope records of benthic foraminifera from ODP Site 1264 in the southeastern Atlantic Ocean are presented which resolve the latest Oligocene to early Miocene (~24–19 Ma) climate changes at high temporal resolution (<3 kyr). Using an inverse modelling technique, we decomposed the oxygen isotope record into temperature and ice volume and found that the Antarctic ice sheet expanded episodically during the declining phase of the long-term (~400 kyr) eccentricity cycle and subsequent low short-term (~100 kyr) eccentricity cycle. The largest glaciations are separated by multiple long-term eccentricity cycles, indicating the involvement of a non-linear response mechanism. Our modelling results suggest that during the largest (Mi-1) event, Antarctic ice sheet volume expanded up to its present-day configuration. In addition, we found that distinct ~100 kyr variability occurs during the termination phases of the major Antarctic glaciations, suggesting that climate and ice-sheet response was more susceptible to short-term eccentricity forcing at these times. During two of these termination-phases, ?18O bottom water gradients in the Atlantic ceased to exist, indicating a direct link between global climate, enhanced ice-sheet instability and major oceanographic reorganisations

The mid-Cretaceous North Atlantic nutrient trap: Black shales and OAEs

Organic-rich sediments are the salient marine sedimentation product in the mid-Cretaceous of the ocean basins formed in the Mesozoic. Oceanic anoxic events (OAEs) are discrete and particularly organic-rich intervals within these mid-Cretaceous organic-rich sequences and are defined by pronounced carbon isotope excursions. Marine productivity during OAEs appears to have been enhanced by the increased availability of biolimiting nutrients in seawater due to hydrothermal alteration of submarine basalts in the Pacific and proto-Indian oceans. The exact mechanisms behind the deposition of organic-rich sediments in the mid-Cretaceous are still a matter of discussion, but a hypothesis which is often put forward is that their deposition was a consequence of the coupling of a particular paleogeography with changes in ocean circulation and nutrient supply. In this study, we used a global coupled climate model to investigate oceanic processes that affect the interbasinal exchange of nutrients as well as their spatial distribution and bioavailability. We conclude that the mid-Cretaceous North Atlantic was a nutrient trap as a consequence of an estuarine circulation with respect to the Pacific. Organic-rich sediments in the North Atlantic were deposited below regions of intense upwelling. We suggest that enhanced productivity during OAEs was a consequence of upwelling of Pacific-derived nutrient-rich seawater associated with submarine igneous events

Deformation and failure of the ice bridge on the Wilkins Ice Shelf, Antarctica

A narrow bridge of floating ice that connected the Wilkins Ice Shelf, Antarctica, to two confining islands eventually collapsed in early April 2009. In the month preceding the collapse, we observed deformation of the ice bridge by means of satellite imagery and from an in situ GPS station. TerraSAR-X images (acquired in stripmap mode) were used to compile a time series. The ice bridge bent most strongly in its narrowest part (westerly), while the northern end (near Charcot Island) shifted in a northeasterly direction. In the south, the ice bridge experienced compressive strain parallel to its long axis. GPS position data were acquired a little south of the narrowest part of the ice bridge from 19 January 2009. Analysis of these data showed both cyclic and monotonic components of motion. Meteorological data and re-analysis of the output of weather-prediction models indicated that easterly winds were responsible for the cyclic motion component. In particular, wind stress on the rough ice melange that occupied the area to the east exerted significant pressure on the ice bridge. The collapse of the ice bridge began with crack formation in the southern section parallel to the long axis of the ice bridge and led to shattering of the southern part. Ultimately, the narrowest part, only 900 m wide, ruptured. The formation of many small icebergs released energy of >125 × 106 J

Constraints on the delta H-2 diffusion rate in firn from field measurements at Summit, Greenland

We performed detailed 2H isotope diffusion measurements in the upper 3 m of firn at Summit, Greenland. Using a small snow gun, a thin snow layer was formed from 2H-enriched water over a 6 × 6 m2 area. We followed the diffusion process, quantified as the increase of the δ2H diffusion length, over a 4-year period, by retrieving the layer once per year by drilling a firn core and slicing it into 1 cm layers and measuring the δ2H signal of these layers. We compared our experimental findings to calculations based on the model by Johnsen et al. (2000) and found substantial differences. The diffusion length in our experiments increased much less over the years than in the model. We discuss the possible causes for this discrepancy and conclude that several aspects of the diffusion process in firn are still poorly constrained, in particular the tortuosity

Будівельна лихоманка на Київських схилах

Forty marine-terminating glaciers have been surveyed daily since 2000 using cloud-free MODIS visible imagery (Box and Decker 2011; http://bprc. osu.edu/MODIS/). The net area change of the 40 glaciers during the period of observation has been -1775 km2, with the 18 northernmost (>72°N) glaciers alone contributing to half of the net area change. In 2012, the northernmost glaciers lost a collective area of 255 km2, or 86% of the total net area change of the 40 glaciers surveyed. The six glaciers with the largest net area loss in 2012 were Petermann (-141 km2), 79 glacier (-27 km2), Zachariae (-26 km2), Steenstrup (-19 km2), Steensby (-16 km2, the greatest retreat since observations began), and Jakobshavn (-13 km2). While the total area change was negative in 2012, the area of four of the forty glaciers did increase relative to the end of the 2011 melt season. The anomalous advance of these four glaciers is not easily explained, as the mechanisms controlling the behavior of individual glaciers are uncertain due to their often unique geographic!settings

A state-dependent quantification of climate sensitivity based in paleo data of the last 2.1 million years

The evidence from both data and models indicates that specific equilibrium climate sensitivity S[X] — the global annual mean surface temperature change (DTg) as a response to a change in radiative forcing X (DR[X]) — is state-dependent. Such a state dependency implies that the best fit in the scatter plot of (DTg versus DR[X] is not a linear regression, but can be some non-linear or even non-smooth function. While for the conventional linear case the slope (gradient) of the regression is correctly interpreted as the specific equilibrium climate sensitivity S[X], the interpretation is not straightforward in the non-linear case. We here explain how such a state-dependent scatter plot needs to be interpreted, and provide a theoretical understanding — or generalization — how to quantify S[X] in the non-linear case. Finally, from data covering the last 2.1 Myr we show that — due to state dependency — the specific equilibrium climate sensitivity which considers radiative forcing of CO2 and land ice sheet (LI) albedo, S[CO2;LI], is larger during interglacial states than during glacial conditions by more than a factor two