Modelling calving front dynamics using a level-set method: application to Jakobshavn Isbræ, West Greenland, links model results in MATLAB format

Calving is a major mechanism of ice discharge of the Antarctic and Greenland ice sheets, and a change in calving front position affects the entire stress regime of marine terminating glaciers. The representation of calving front dynamics in a 2-D or 3-D ice sheet model remains non-trivial. Here, we present the theoretical and technical framework for a level-set method, an implicit boundary tracking scheme, which we implement into the Ice Sheet System Model (ISSM). This scheme allows us to study the dynamic response of a drainage basin to user-defined calving rates. We apply the method to Jakobshavn Isbræ, a major marine terminating outlet glacier of the West Greenland Ice Sheet. The model robustly reproduces the high sensitivity of the glacier to calving, and we find that enhanced calving triggers significant acceleration of the ice stream. Upstream acceleration is sustained through a combination of mechanisms. However, both lateral stress and ice influx stabilize the ice stream. This study provides new insights into the ongoing changes occurring at Jakobshavn Isbræ and emphasizes that the incorporation of moving boundaries and dynamic lateral effects, not captured in flow-line models, is key for realistic model projections of sea level rise on centennial timescales.

Calcareous nannofossils across the Eocene-Oligocene boundary

This work presents the stratigraphic distribution of several species of calcareous nannofossil in the middle Eocene early-Oligocene from four Ocean Drilling Program (ODP) sites located between 60° and 65°S paleolatitude in the Southern Atlantic and Indian Oceans. Useful nannofossil datums that should facilitate construction of age-models and contribute to an integrated chronology for the upper Paleogene Southern Ocean sediments from ~42 to 33 Ma are summarized. The distribution patterns of calcareous nannofossils, studied by means of quantitative and semiquantitative methods, provide an improvement of the classical Southern Ocean biozonations, introducing new biostratigraphically useful biohorizons, and testing their reproducibility within and outside the region.

Larcform 1 - single column model intercomparison of Arctic air formation, link to model results in NetCDF format

Weather and climate models struggle to represent lower tropospheric temperature and moisture profiles and surface fluxes in Arctic winter, not least because they lack or misrepresent physical processes that are specific to high latitudes. The Arctic boundary layer in winter has been observed to be in either a radiatively clear or cloudy state: The radiatively clear state is characterized by strong surface radiative cooling leading to the build-up of surface-based temperature inversions, whereas the cloudy state occurs when cloud liquid water is present in the atmospheric column, allowing little or no surface radiative cooling and leading to weaker and typically elevated temperature inversions. Many large-scale models have been shown to lack the cloudy state, and some do substantially underestimate stability in the clear state. We here present results from the first Lagrangian ARCtic air FORMation experiment (Larcform 1), a GASS (Global atmospheric system studies) single-column model intercomparison which reproduces these biases of large-scale models in an idealised setup.

Simulated Bromoform emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by 6-hourly NCEP data

Bromoform (CHBr3) is one important precursor of atmospheric reactive bromine species that are involved in ozone depletion in the troposphere and stratosphere. In the open ocean bromoform production is linked to phytoplankton that contains the enzyme bromoperoxidase. Coastal sources of bromoform are higher than open ocean sources. However, open ocean emissions are important because the transfer of tracers into higher altitude in the air, i.e. into the ozone layer, strongly depends on the location of emissions. For example, emissions in the tropics are more rapidly transported into the upper atmosphere than emissions from higher latitudes. Global spatio-temporal features of bromoform emissions are poorly constrained. Here, a global three-dimensional ocean biogeochemistry model (MPIOM-HAMOCC) is used to simulate bromoform cycling in the ocean and emissions into the atmosphere using recently published data of global atmospheric concentrations (Ziska et al., 2013) as upper boundary conditions. Our simulated surface concentrations of CHBr3 match the observations well. Simulated global annual emissions based on monthly mean model output are lower than previous estimates, including the estimate by Ziska et al. (2013), because the gas exchange reverses when less bromoform is produced in non-blooming seasons. This is the case for higher latitudes, i.e. the polar regions and northern North Atlantic. Further model experiments show that future model studies may need to distinguish different bromoform-producing phytoplankton species and reveal that the transport of CHBr3 from the coast considerably alters open ocean bromoform concentrations, in particular in the northern sub-polar and polar regions.

Age model, raw data and interpolated raw data from proxy records of sediment cores GeoB6007-1 and GeoB6007-2

Here we present a 1200 yr long benthic foraminiferal Mg/Ca based temperature and oxygen isotope record from a ~900 m deep sediment core off northwest Africa to show that atmosphere-ocean interactions in the eastern subpolar gyre are transferred at central water depth into the eastern boundary of the subtropical gyre. Further we link the variability of the NAO (over the past 165 yrs) and solar irradiance (Late Holocene) and their control on subpolar mode water formation to the multidecadal variability observed at mid-depth in the eastern subtropical gyre. Our results show that eastern North Atlantic central waters cooled by up to ~0.8± 0.7 °C and densities decreased by Sigma theta=0.3±0.2 during positive NAO years and during minima in solar irradiance during the Late Holocene. The presented records demonstrate the sensitivity of central water formation to enhanced atmospheric forcing and ice/freshwater fluxes into the eastern subpolar gyre and the importance of central water circulation for cross-gyre climate signal propagation during the Late Holocene.

Age model, lipids, and iron/calcium-ratios of sediment core GeoB9307-3

Hide Intense debate persists about the climatic mechanisms governing hydrologic changes in tropical and subtropical southeast Africa since the Last Glacial Maximum, about 20,000 years ago. In particular, the relative importance of atmospheric and oceanic processes is not firmly established. Southward shifts of the intertropical convergence zone (ITCZ) driven by high-latitude climate changes have been suggested as a primary forcing, whereas other studies infer a predominant influence of Indian Ocean sea surface temperatures on regional rainfall changes. To address this question, a continuous record representing an integrated signal of regional climate variability is required, but has until now been missing. Here we show that remote atmospheric forcing by cold events in the northern high latitudes appears to have been the main driver of hydro-climatology in southeast Africa during rapid climate changes over the past 17,000 years. Our results are based on a reconstruction of precipitation and river discharge changes, as recorded in a marine sediment core off the mouth of the Zambezi River, near the southern boundary of the modern seasonal ITCZ migration. Indian Ocean sea surface temperatures did not exert a primary control over southeast African hydrologic variability. Instead, phases of high precipitation and terrestrial discharge occurred when the ITCZ was forced southwards during Northern Hemisphere cold events, such as Heinrich stadial 1 (around 16,000 years ago) and the Younger Dryas (around 12,000 years ago), or when local summer insolation was high in the late Holocene, i.e., during the last 4,000 years.

Bulk stable carbon and oxygen isotopes from ODP Site 1258, and Fe intensities from ODP Site 1262, and minimal tuning age model options for ODP Site 1258 and 1262

Timing is crucial to understanding the causes and consequences of events in Earth history. The calibration of geological time relies heavily on the accuracy of radioisotopic and astronomical dating. Uncertainties in the computations of Earth's orbital parameters and in radioisotopic dating have hampered the construction of a reliable astronomically calibrated time scale beyond 40 Ma. Attempts to construct a robust astronomically tuned time scale for the early Paleogene by integrating radioisotopic and astronomical dating are only partially consistent. Here, using the new La2010 and La2011 orbital solutions, we present the first accurate astronomically calibrated time scale for the early Paleogene (47-65 Ma) uniquely based on astronomical tuning and thus independent of the radioisotopic determination of the Fish Canyon standard. Comparison with geological data confirms the stability of the new La2011 solution back to ~54 Ma. Subsequent anchoring of floating chronologies to the La2011 solution using the very long eccentricity nodes provides an absolute age of 55.530 {plus minus} 0.05 Ma for the onset of the Paleocene/Eocene Thermal Maximum (PETM), 54.850 {plus minus} 0.05 Ma for the early Eocene ash -17, and 65.250 {plus minus} 0.06 Ma for the K/Pg boundary. The new astrochronology presented here indicates that the intercalibration and synchronization of U/Pb and 40Ar/39Ar radiometric geochronology is much more challenging than previously thought.

Clay mineralogical assemblages, Sm-Nd concentrations, Sm and Nd isotopic ratios and oxygen isotope ratios, OIS number and the age model of ODP Site 105-646

We present 40 Sm-Nd isotope measurements of the clay-size (<2 µm) fractions of sediments from the Southern Greenland rise (ODP-646) that span the last 365 kyr. These data track changes in the relative supply of fine particles carried into the deep Labrador Sea by the Western Boundary Under Current (WBUC) back to the fourth glacial-interglacial cycles. Earlier studies revealed three general sources of particles to the core site: (i) Precambrian crustal material from Canada, Greenland, and/or Scandinavia (North American Shield - NAS), (ii) Palaeozoic or younger crustal material from East Greenland, NW Europe, and/or western Scandinavia (Young Crust - YC) and (iii) volcanic material from Iceland and the Mid-Atlantic Ridge (MAR). Clay-size fractions from glacial sediments have the lowest Nd isotopic ratios. Supplies of young crustal particles were similar during glacial oxygen isotope stages (OIS) 2, 6, and 10. In contrast the mean volcanic contributions decreased relative to old craton material from OIS 10 to OIS 6 and then from OIS 6 to OIS 2. The glacial OIS 8 interval displays a mean Sm/Nd ratio similar to those of interglacials OIS 1, 5, and 9. Compared with other interglacials, OIS 7 was marked by a higher YC contribution but a similar ~30% MAR supply. The overall NAS contribution dropped by a factor of 2 during each glacial/interglacial transition, with the MAR contribution broadly replacing it during interglacials. To decipher between higher supplies and/or dilution, particle fluxes from each end member were estimated. Glacial NAS fluxes were systematically higher than interglacial fluxes. During the time interval examined, fine particle supplies to the Labrador Sea were strongly controlled by proximal ice-margin erosion and thus echoed the glacial stage intensity. In contrast, the WBUC-carried MAR supplies from the eastern basins did not change significantly throughout the last 365 kyr, except for a marked increase in surface-sediments that suggests unique modern conditions. Distal WBUC-controlled inputs from the Northern and NE North Atlantic seem to have been less variable than proximal supplies linked with glacial erosion rate.