Elevation Models of Antarctica and Greenland derived from CryoSat-2 in the period 2011 to 2013

This study focuses on the present-day surface elevation of the Greenland and Antarctic ice sheets. Based on 3 years of CryoSat-2 data acquisition we derived new elevation models (DEMs) as well as elevation change maps and volume change estimates for both ice sheets. Here we present the new DEMs and their corresponding error maps. The accuracy of the derived DEMs for Greenland and Antarctica is similar to those of previous DEMs obtained by satellite-based laser and radar altimeters. Comparisons with ICESat data show that 80% of the CryoSat-2 DEMs have an uncertainty of less than 3 m ± 15 m. The surface elevation change rates between January 2011 and January 2014 are presented for both ice sheets. We compared our results to elevation change rates obtained from ICESat data covering the time period from 2003 to 2009. The comparison reveals that in West Antarctica the volume loss has increased by a factor of 3. It also shows an anomalous thickening in Dronning Maud Land, East Antarctica which represents a known large-scale accumulation event. This anomaly partly compensates for the observed increased volume loss of the Antarctic Peninsula and West Antarctica. For Greenland we find a volume loss increased by a factor of 2.5 compared to the ICESat period with large negative elevation changes concentrated at the west and southeast coasts. The combined volume change of Greenland and Antarctica for the observation period is estimated to be -503 ± 107 km**3/yr. Greenland contributes nearly 75% to the total volume change with -375 ± 24 km**3/yr.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (models indicating extremely rapid release of isotopically light carbon, possibly from seafloor methane hydrate, as the proximal cause of the event. However, the 3He[ET] technique indicates a previously unrecognized and extreme increase in sedimentation rate coincident with the return of climate proxies to pre-event values. The 3He[ET]-based age model thus suggests a far more rapid recovery from the climatic perturbation than previously proposed or predicted on the basis of the modern carbon cycle, and so may indicate additional or accelerated mechanisms of carbon removal from the ocean-atmosphere system during this period. 3He[ET] was also measured at ODP Site 1051 to test the validity of the Site 690 chronology. Comparison of these data sets seems to require removal of several tens of kyr of sediment within the climatic excursion at Site 1051, an observation consistent with sediment structures and previous age modeling efforts. The Site 1051 age model shows a ~30 kyr period in which climate proxies return toward pre-event values, after which they remain invariant for ~80 kyr. If this rise represents the recovery interval identified at Site 690, then the 3HeET-based age models of the two sites are in good agreement. However, alternative interpretations are possible, and work on less disrupted sites is required to evaluate the reliability of the proposed new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.

Three species distribution models for the marine diatom Fragilariopsis kerguelensis of the Southern Ocean, links to NetCDF files

Fragilariopsis kerguelensis, a dominant diatom species throughout the Antarctic Circumpolar Current, is coined to be one of the main drivers of the biological silicate pump. Here, we study the distribution of this important species and expected consequences of climate change upon it, using correlative species distribution modeling and publicly available presence-only data. As experience with SDM is scarce for marine phytoplankton, this also serves as a pilot study for this organism group. We used the maximum entropy method to calculate distribution models for the diatom F. kerguelensis based on yearly and monthly environmental data (sea surface temperature, salinity, nitrate and silicate concentrations). Observation data were harvested from GBIF and the Global Diatom Database, and for further analyses also from the Hustedt Diatom Collection (BRM). The models were projected on current yearly and seasonal environmental data to study current distribution and its seasonality. Furthermore, we projected the seasonal model on future environmental data obtained from climate models for the year 2100. Projected on current yearly averaged environmental data, all models showed similar distribution patterns for F. kerguelensis. The monthly model showed seasonality, for example, a shift of the southern distribution boundary toward the north in the winter. Projections on future scenarios resulted in a moderately to negligibly shrinking distribution area and a change in seasonality. We found a substantial bias in the publicly available observation datasets, which could be reduced by additional observation records we obtained from the Hustedt Diatom Collection. Present-day distribution patterns inferred from the models coincided well with background knowledge and previous reports about F. kerguelensis distribution, showing that maximum entropy-based distribution models are suitable to map distribution patterns for oceanic planktonic organisms. Our scenario projections indicate moderate effects of climate change upon the biogeography of F. kerguelensis.

Temperature observation of CORK installation during ODP Leg 168

Four boreholes, drilled a few tens of meters into igneous basement on the eastern flank of the Juan de Fuca Ridge during ODP Leg 168, were sealed and instrumented for long-term monitoring to observe the hydrologic state of young sediment-sealed oceanic crust. The thermal regime is dominated by the effects of rapid fluid circulation in uppermost igneous basement driven by very small non-hydrostatic pressure gradients. Upper basement temperatures are uniform laterally between pairs of holes over distances of hundreds of meters to kilometers. In the case of two holes drilled into a sediment-buried basement ridge and adjacent valley, basement temperatures differ by less than 2 K despite the 2.2 km lateral separation of the sites and the 2.5:1 contrast in sediment cover thickness. Under conductive conditions, upper basement temperatures would differ by roughly 50 K. By comparison with modeling results, the observed degree of isothermality suggests a fluid flux of at least 10**-6 m/s (30 m/yr), and an effective permeability in the range of 10**-10 to 10**-9 m**2 in the uppermost igneous crust. The pressure difference available to drive this rapid flux between the ridge and valley, estimated by comparing the observed pressures via the isothermal upper basement hydrostat that is inferred to connect the two sites, is small (~2 kPa) and also suggests high permeability. Relative to the hydrostats defined by the local conductive sediment geotherms, substantial super-hydrostatic pressure (+18 kPa) is present within the buried basement ridge, and sub-hydrostatic pressure is present in the adjacent valley (-26 kPa). Such pressure differentials are the direct consequence of the advection-dominated thermal regime and small pressure losses in high-permeability basement, and are available to drive fluid seepage through sediment sections vertically up above and horizontally away from buried ridges, and down above valleys. No constraints are provided by any of the observations on the depth in the crust to which thermally or chemically significant flow might extend, although just as in the overlying sediments, the pattern of deep flow may be affected by the near-isothermal and near-hydrostatic conditions present in the permeable uppermost crustal section.

Age models for sediment cores off Peru

Planktonic foraminiferal faunas of the southeast Pacific indicate that sea surface temperatures (SST) have varied by as much as 8-10°C in the Peru Current, and by ?5-7°C along the equator, over the past 150,000 years. Changes in SST at times such as the Last Glacial Maximum reflect incursion of high-latitude species Globorotalia inflata and Neogloboquadrina pachyderma into the eastern boundary current and as far north as the equator. A simple heat budget model of the equatorial Pacific shows that observed changes in Peru Current advection can account for about half of the total variability in equatorial SSTs. The remaining changes in equatorial SST, which are likely related to local changes in upwelling or pycnocline depth, precede changes in polar climates as recorded by d18O. This partitioning of processes in eastern equatorial Pacific SST reveals that net ice-age cooling here reflects first a rapid response of equatorial upwelling to insolation, followed by a later response to changes in the eastern boundary current associated with high-latitude climate (which closely resembles variations in atmospheric CO2 as recorded in the Vostok ice core). Although precise mechanisms responsible for the equatorial upwelling component of climate change remain uncertain, one likely candidate that may operate independently of the ice sheets is insolation-driven changes in El Niño/Southern Oscillation (ENSO) frequency. Early responses of equatorial SST detected both here and elsewhere highlight the sensitivity of tropical systems to small changes in seasonal insolation. The scale of tropical changes we have observed are substantially greater than model predictions, suggesting a need for further quantitative assessment of processes associated with long-term climate change.

Age models and stable isotope record of sediment cores from the Atlantic Ocean

The high-resolution delta18O and delta13C records of benthic foraminifera from a 150,000-year long core from the Caribbean Sea indicate that there was generally high delta13C during glaciations and low delta13C during interglaciations. Due to its 1800-m sill depth, the properties of deep water in the Caribbean Sea are similar to those of middepth tropical Atlantic water. During interglaciations, the water filling the deep Caribbean Sea is an admixture of low delta13C Upper Circumpolar Water (UCPW) and high delta13C Upper North Atlantic Deep Water (UNADW). By contrast, only high delta13C UNADW enters during glaciations. Deep ocean circulation changes can influence atmospheric CO2 levels (Broecker and Takahashi, 1985; Boyle, 1988 doi:10.1029/JC093iC12p15701; Keir, 1988 doi:10.1029/PA003i004p00413; Broecker and Peng, 1989 doi:10.1029/GB003i003p00215). By comparing delta13C records of benthic foraminifera from cores lying in Southern Ocean Water, the Caribbean Sea, and at several other Atlantic Ocean sites, the thermohaline state of the Atlantic Ocean (how close it was to a full glacial or full interglacial configuration) is characterized. A continuum of circulation patterns between the glacial and interglacial extremes appears to have existed in the past. Subtracting the deep Pacific (~mean ocean water) delta13C record from the Caribbean delta13C record yields a record which describes large changes in the Atlantic Ocean thermohaline circulation. The delta13C difference varies as the vertical nutrient distribution changes. This new proxy record bears a striking resemblance to the 150,000-year-long atmospheric CO2 record (Barnola et al., 1987 doi:10.1038/329408a0). This favorable comparison between the new proxy record and the atmospheric CO2 record is consistent with Boyle's (1988a) model that vertical nutrient redistribution has driven large atmospheric CO2 changes in the past. Changes in the relative contribution of NADW and Pacific outflow water to the Southern Ocean are also consistent with Broecker and Peng's (1989) recent model for atmospheric CO2 changes.

Age models, accumulation rates, and grain size distributions of sediment cores from the northern subtropical Atlantic

The terrigenous sediment proportion of the deep sea sediments from off Northwest Africa has been studied in order to distinguish between the aeolian and the fluvial sediment supply. The present and fossil Saharan dust trajectories were recognized from the distribution patterns of the aeolian sediment. The following timeslices have been investigated: Present, 6,000, 12,000 and 18,000 y. B. P. Furthermore, the quantity of dust deposited off the Saharan coast has been estimated. For this purpose, 80 surface sediment samples and 34 sediment cores have been analysed. The stratigraphy of the cores has been achieved from oxygen isotopic curves, 14C-dating, foraminiferal transfer temperatures, and carbonate contents. Silt sized biogenic opal generally accounts for less than 2 % of the total insoluble sediment proportion. Only under productive upwelling waters and off river mouths, the opal proportion exceeds 2 % significantly. The modern terrigenous sediment from off the Saharan coast is generally characterized by intensely stained quartz grains. They indicate an origin from southern Saharan and Sahelian laterites, and a zonal aeolian transport in midtropospheric levels, between 1.5 an 5.5 km, by 'Harmattan' Winds. The dust particles follow large outbreaks of Saharan air across the African coast between 15° and 21° N. Their trajectories are centered at about 18° N and continue further into a clockwise gyre situated south of the Canary Islands. This course is indicated by a sickle-shaped tongue of coarser grain sizes in the deep-sea sediment. Such loess-sized terrigenous particles only settle within a zone extending to 700 km offshore. Fine silt and clay sized particles, with grain sizes smaller than 10- 15 µm, drift still further west and can be traced up to more than 4,000 km distance from their source areas. Additional terrigenous silt which is poor in stained quartz occurs within a narrow zone off the western Sahara between 20° and 27° N only. It depicts the present dust supply by the trade winds close to the surface. The dust load originates from the northwestern Sahara, the Atlas Mountains and coastal areas, which contain a particularly low amount of stained quartz. The distribution pattern of these pale quartz sediments reveals a SSW-dispersal of dust being consistent with the present trade wind direction from the NNE. In comparison to the sediments from off the Sahara and the deeper subtropical Atlantic, the sediments off river mouths, in particular off the Senegal river, are characterized by an additional input of fine grained terrigenous particles (< 6 µm). This is due to fluvial suspension load. The fluvial discharge leads to a relative excess of fine grained particles and is observed in a correlation diagram of the modal grain sizes of terrigenous silt with the proportion of fine fraction (< 6 µm). The aeolian sediment contribution by the Harmattan Winds strongly decreased during the Climatic Optimum at 6,000 y. B. P. The dust discharge of the trade winds is hardly detectable in the deep-sea sediments. This probably indicates a weakened atmospheric circulation. In contrast, the fluvial sediment supply reached a maximum, and can be traced to beyond Cape Blanc. Thus, the Saharan climate was more humid at 6,000 y B. P. A latitudinal shift of the Harmattan driven dust outbreaks cannot be observed. Also during the Glacial, 18,000 y. B. P., Harmattan dust transport crossed the African coast at latitudes of 15°-20° N. Its sediment load increased intensively, and markedly coarser grains spread further into the Atlantic Ocean. An expanded zone of pale-quart sediments indicates an enhanced dust supply by the trade winds blowing from the NE. No synglacial fluvial sediment contribution can be recognized between 12° and 30° N. This indicates a dry glacial climate and a strengthened stmospheric circulation over the Sahelian and Saharan region. The climatic transition pahes, at 12, 000 y. B. P., between the last Glacial and the Intergalcial, which is compareable to the Alerod in Europe, is characterized by an intermediate supply of terrigenous particles. The Harmattan dust transport wa weaker than during the Glacial. The northeasterly trade winds were still intensive. River supply reached a first postglacial maximum seaward of the Senegal river mouth. This indicates increasing humidity over the southern Sahara and a weaker atmospheric circulation as compared to the glacial. The accumulation rates of the terrigenous silt proportion (> 6 µm) decrcase exponentially with increasing distance from the Saharan coast. Those of the terrigenous fine fraction (< 6 µm) follow the same trend and show almost similar gradients. Accordingly, also the terrigenous fine fraction is believed to result predominantly from aeolian transport. In the Atlantic deep-sea sediments, the annual terrigenous sediment accumulation has fluctuated, from about 60 million tons p. a. during the Late Glacial (13,500-18,000 y. B. P, aeolian supply only) to about 33 million tons p. a. during the Holocene Climatic Optimum (6,000-9,000 y. B. P, mainly fluvial supply), when the river supply has reached a maximum, and to about 45 million tons p. a. during the last 4,000 years B. P. (fluvial supply only south of 18° N).

Surface velocity fields, digital elevation models, ice front positions and grounding line derived from remote sensing data at Dinsmoor-Bombardier-Edgeworth glacier system, Antarctic Peninsula (1992-2014)

The northern Antarctic Peninsula is one of the fastest changing regions on Earth. The disintegration of the Larsen-A Ice Shelf in 1995 caused tributary glaciers to adjust by speeding up, surface lowering, and overall increased ice-mass discharge. In this study, we investigate the temporal variation of these changes at the Dinsmoor-Bombardier-Edgeworth glacier system by analyzing dense time series from various spaceborne and airborne Earth observation missions. Precollapse ice shelf conditions and subsequent adjustments through 2014 were covered. Our results show a response of the glacier system some months after the breakup, reaching maximum surface velocities at the glacier front of up to 8.8 m/d in 1999 and a subsequent decrease to ~1.5 m/d in 2014. Using a dense time series of interferometrically derived TanDEM-X digital elevation models and photogrammetric data, an exponential function was fitted for the decrease in surface elevation. Elevation changes in areas below 1000 m a.s.l. amounted to at least 130±15 m130±15 m between 1995 and 2014, with change rates of ~3.15 m/a between 2003 and 2008. Current change rates (2010-2014) are in the range of 1.7 m/a. Mass imbalances were computed with different scenarios of boundary conditions. The most plausible results amount to -40.7±3.9 Gt-40.7±3.9 Gt. The contribution to sea level rise was estimated to be 18.8±1.8 Gt18.8±1.8 Gt, corresponding to a 0.052±0.005 mm0.052±0.005 mm sea level equivalent, for the period 1995-2014. Our analysis and scenario considerations revealed that major uncertainties still exist due to insufficiently accurate ice-thickness information. The second largest uncertainty in the computations was the glacier surface mass balance, which is still poorly known. Our time series analysis facilitates an improved comparison with GRACE data and as input to modeling of glacio-isostatic uplift in this region. The study contributed to a better understanding of how glacier systems adjust to ice shelf disintegration.