Results of modelled Arctic river discharge trends since 7 ka BP

The Arctic hydrological cycle throughout the Holocene is analyzed based on the results of transient simulations with the coupled atmosphere-ocean circulation model ECHO-G. The results suggest a ~ 2 % increase of mid-Holocene to preindustrial Arctic river discharges for the Eurasian continent. However, rivers of the North America Arctic realm show a moderate runoff decline of approximately 4 to 5 % for the same period. The total river discharge into the Arctic Ocean has remained at an approximately constant preindustrial level since the mid Holocene. The positive discharge trend within Eurasia is caused by a more rapid decrease in local net evaporation compared to a smaller decline in advected moisture and hence precipitation. This effect is neither recognized within the North American Arctic domain nor in the far eastern part of the Eurasian Arctic realm. A detailed comparison of these model findings with a variety of proxy studies is conducted. The collected proxy records show trends of continental surface temperatures and precipitation rates that are consistent with the simulations. A continuation of the transient Holocene runs for the 19th and 20th century with increased greenhouse gases indicates an increase of the total river influx into the Arctic Ocean of up to 7.6 %. The Eurasian river discharges increase by 7.5 %, the North American discharges by up to 8.4 %. The most rapid increases have been detected since the beginning of the 20th century. These results are corroborated by the observed rising of Arctic river discharges during the last century which is attributed to anthropogenic warming. The acceleration of the Arctic hydrological cycle in the 20th century is without precedence in the Holocene.

Sea surface temperatures during Marine Isotope Stage (MIS) 11

The Marine Isotope Stage (MIS) 11 (424-374 ka) was characterized by a protracted deglaciation and an unusually long climatic optimum. It remains unclear to what degree the climate development during this interglacial reflects the unusually weak orbital forcing or greenhouse gas trends. Previously, arguments about the duration and timing of the MIS11 climatic optimum and about the pace of the deglacial warming were based on a small number of key records, which appear to show regional differences. In order to obtain a global signal of climate evolution during MIS11, we compiled a database of 78 sea surface temperature (SST) records from 57 sites spanning MIS11, aligned these individually on the basis of benthic (N = 28) or planktonic (N = 31) stable oxygen isotope curves to a common time frame and subjected 48 of them to an empirical orthogonal function (EOF) analysis. The analysis revealed a high commonality among all records, with the principal SST trend explaining almost 49% of the variability. This trend indicates that on the global scale, the surface ocean underwent rapid deglacial warming during Termination V, in pace with carbon dioxide rise, followed by a broad SST optimum centered at ~410 kyr. The second EOF, which explained ~18% of the variability, revealed the existence of a different SST trend, characterized by a delayed onset of the temperature optimum during MIS11 at ~398 kyr, followed by a prolonged warm period lasting beyond 380 kyr. This trend is most consistently manifested in the mid-latitude North Atlantic and Mediterranean Sea and is here attributed to the strength of the Atlantic meridional overturning circulation. A sensitivity analysis indicates that these results are robust to record selection and to age-model uncertainties of up to 3-6 kyr, but more sensitive to SST seasonal attribution and SST uncertainties >1 °C. In order to validate the CCSM3 (Community Climate System Model, version 3) predictive potential, the annual and seasonal SST anomalies recorded in a total of 74 proxy records were compared with runs for three time slices representing orbital configuration extremes during the peak interglacial of MIS11. The modeled SST anomalies are characterized by a significantly lower variance compared to the reconstructions. Nevertheless, significant correlations between proxy and model data are found in comparisons on the seasonal basis, indicating that the model captures part of the long-term variability induced by astronomical forcing, which appears to have left a detectable signature in SST trends.

Nd bulk isotopes and 40Ar/39Ar ages of individual minerals from Aurora and Wilkes sub-glacial basin sediments

The Wilkes and Aurora basins are large, low-lying sub-glacial basins that may cause areas of weakness in the overlying East Antarctic ice sheet. Previous work based on ice-rafted debris (IRD) provenance analyses found evidence for massive iceberg discharges from these areas during the late Miocene and Pliocene. Here we characterize the sediments shed from the inferred areas of weakness along this margin (94°E to 165°E) by measuring40Ar/39Ar ages of 292 individual detrital hornblende grains from eight marine sediment core locations off East Antarctica and Nd isotopic compositions of the bulk fine fraction from the same sediments. We further expand the toolbox for Antarctic IRD provenance analyses by exploring the application of 40Ar/39Ar ages of detrital biotites; biotite as an IRD tracer eliminates lithological biases imposed by only analyzing hornblendes and allows for characterization of samples with low IRD concentrations. Our data quadruples the number of detrital 40Ar/39Ar ages from this margin of East Antarctica and leads to the following conclusions: (1) Four main sectors between the Ross Sea and Prydz Bay, separated by ice drainage divides, are distinguishable based upon the combination of 40Ar/39Ar ages of detrital hornblende and biotite grains and the e-Nd of the bulk fine fraction; (2) 40Ar/39Ar biotite ages can be used as a robust provenance tracer for this part of East Antarctica; and (3) sediments shed from the coastal areas of the Aurora and Wilkes sub-glacial basins can be clearly distinguished from one another based upon their isotopic fingerprints.

Glacier mass balances of Stubacher Sonnblickkees, Hohe Tauern Range, Eastern Alps, Austria, 1958/1959 to 2012/2013

Stubacher Sonnblickkees (SSK) is located in the Hohe Tauern Range (Eastern Alps) in the south of Salzburg Province (Austria) in the region of Oberpinzgau in the upper Stubach Valley. The glacier is situated at the main Alpine crest and faces east, starting at elevations close to 3050 m and in the 1980s terminated at 2500 m a.s.l. It had an area of 1.7 km² at that time, compared with 1 km² in 2013. The glacier type can be classified as a slope glacier, i.e. the relief is covered by a relatively thin ice sheet and there is no regular glacier tongue. The rough subglacial topography makes for a complex shape in the surface topography, with various concave and convex patterns. The main reason for selecting this glacier for mass balance observations (as early as 1963) was to verify on a complex glacier how the mass balance methods and the conclusions - derived during the more or less pioneer phase of glaciological investigations in the 1950s and 1960s - could be applied to the SSK glacier. The decision was influenced by the fact that close to the SSK there was the Rudolfshütte, a hostel of the Austrian Alpine Club (OeAV), newly constructed in the 1950s to replace the old hut dating from 1874. The new Alpenhotel Rudolfshütte, which was run by the Slupetzky family from 1958 to 1970, was the base station for the long-term observation; the cable car to Rudolfshütte, operated by the Austrian Federal Railways (ÖBB), was a logistic advantage. Another factor for choosing SSK as a glaciological research site was the availability of discharge records of the catchment area from the Austrian Federal Railways who had turned the nearby lake Weißsee ('White Lake') - a former natural lake - into a reservoir for their hydroelectric power plants. In terms of regional climatic differences between the Central Alps in Tyrol and those of the Hohe Tauern, the latter experienced significantly higher precipitation , so one could expect new insights in the different response of the two glaciers SSK and Hintereisferner (Ötztal Alps) - where a mass balance series went back to 1952. In 1966 another mass balance series with an additional focus on runoff recordings was initiated at Vernagtfener, near Hintereisferner, by the Commission of the Bavarian Academy of Sciences in Munich. The usual and necessary link to climate and climate change was given by a newly founded weather station (by Heinz and Werner Slupetzky) at the Rudolfshütte in 1961, which ran until 1967. Along with an extension and enlargement to the so-called Alpine Center Rudolfshütte of the OeAV, a climate observatory (suggested by Heinz Slupetzky) has been operating without interruption since 1980 under the responsibility of ZAMG and the Hydrological Service of Salzburg, providing long-term met observations. The weather station is supported by the Berghotel Rudolfshütte (in 2004 the OeAV sold the hotel to a private owner) with accommodation and facilities. Direct yearly mass balance measurements were started in 1963, first for 3 years as part of a thesis project. In 1965 the project was incorporated into the Austrian glacier measurement sites within the International Hydrological Decade (IHD) 1965 - 1974 and was afterwards extended via the International Hydrological Program (IHP) 1975 - 1981. During both periods the main financial support came from the Hydrological Survey of Austria. After 1981 funds were provided by the Hydrological Service of the Federal Government of Salzburg. The research was conducted from 1965 onwards by Heinz Slupetzky from the (former) Department of Geography of the University of Salzburg. These activities received better recognition when the High Alpine Research Station of the University of Salzburg was founded in 1982 and brought in additional funding from the University. With recent changes concerning Rudolfshütte, however, it became unfeasible to keep the research station going. Fortunately, at least the weather station at Rudolfshütte is still operating. In the pioneer years of the mass balance recordings at SSK, the main goal was to understand the influence of the complicated topography on the ablation and accumulation processes. With frequent strong southerly winds (foehn) on the one hand, and precipitation coming in with storms from the north to northwest, the snow drift is an important factor on the undulating glacier surface. This results in less snow cover in convex zones and in more or a maximum accumulation in concave or flat areas. As a consequence of the accentuated topography, certain characteristic ablation and accumulation patterns can be observed during the summer season every year, which have been regularly observed for many decades . The process of snow depletion (Ausaperung) runs through a series of stages (described by the AAR) every year. The sequence of stages until the end of the ablation season depends on the weather conditions in a balance year. One needs a strong negative mass balance year at the beginning of glacier measurements to find out the regularities; 1965, the second year of observation resulted in a very positive mass balance with very little ablation but heavy accumulation. To date it is the year with the absolute maximum positive balance in the entire mass balance series since 1959, probably since 1950. The highly complex ablation patterns required a high number of ablation stakes at the beginning of the research and it took several years to develop a clearer idea of the necessary density of measurement points to ensure high accuracy. A great number of snow pits and probing profiles (and additional measurements at crevasses) were necessary to map the accumulation area/patterns. Mapping the snow depletion, especially at the end of the ablation season, which coincides with the equilibrium line, is one of the main basic data for drawing contour lines of mass balance and to calculate the total mass balance (on a regular-shaped valley glacier there might be an equilibrium line following a contour line of elevation separating the accumulation area and the ablation area, but not at SSK). - An example: in 1969/70, 54 ablation stakes and 22 snow pits were used on the 1.77 km² glacier surface. In the course of the study the consistency of the accumulation and ablation patterns could be used to reduce the number of measurement points. - At the SSK the stratigraphic system, i.e. the natural balance year, is used instead the usual hydrological year. From 1964 to 1981, the yearly mass balance was calculated by direct measurements. Based on these records of 17 years, a regression analysis between the specific net mass balance and the ratio of ablation area to total area (AAR) has been used since then. The basic requirement was mapping the maximum snow depletion at the end of each balance year. There was the advantage of Heinz Slupetzky's detailed local and long-term experience, which ensured homogeneity of the series on individual influences of the mass balance calculations. Verifications took place as often as possible by means of independent geodetic methods, i.e. monoplotting , aerial and terrestrial photogrammetry, more recently also the application of PHOTOMODELLER and laser scans. The semi-direct mass balance determinations used at SSK were tentatively compared with data from periods of mass/volume change, resulting in promising first results on the reliability of the method. In recent years re-analyses of the mass balance series have been conducted by the World Glacier Monitoring Service and will be done at SSK too. - The methods developed at SSK also add to another objective, much discussed in the 1960s within the community, namely to achieve time- and labour-saving methods to ensure continuation of long-term mass balance series. The regression relations were used to extrapolate the mass balance series back to 1959, the maximum depletion could be reconstructed by means of photographs for those years. R. Günther (1982) calculated the mass balance series of SSK back to 1950 by analysing the correlation between meteorological data and the mass balance; he found a high statistical relation between measured and determined mass balance figures for SSK. In spite of the complex glacier topography, interesting empirical experiences were gained from the mass balance data sets, giving a better understanding of the characteristics of the glacier type, mass balance and mass exchange. It turned out that there are distinct relations between the specific net balance, net accumulation (defined as Bc/S) and net ablation (Ba/S) to the AAR, resulting in characteristic so-called 'turnover curves'. The diagram of SSK represents the type of a glacier without a glacier tongue. Between 1964 and 1966, a basic method was developed, starting from the idea that instead of measuring years to cover the range between extreme positive and extreme negative yearly balances one could record the AAR/snow depletion/Ausaperung during one or two summers. The new method was applied on Cathedral Massif Glacier, a cirque glacier with the same area as the Stubacher Sonnblickkees, in British Columbia, Canada. during the summers of 1977 and 1978. It returned exactly the expected relations, e.g. mass turnover curves, as found on SSK. The SSK was mapped several times on a scale of 1:5000 to 1:10000. Length variations have been measured since 1960 within the OeAV glacier length measurement programme. Between 1965 and 1981, there was a mass gain of 10 million cubic metres. With a time lag of 10 years, this resulted in an advance until the mid-1980s. Since 1982 there has been a distinct mass loss of 35 million cubic metres by 2013. In recent years, the glacier has disintegrated faster, forced by the formation of a periglacial lake at the glacier terminus and also by the outcrops of rocks (typical for the slope glacier type), which have accelerated the meltdown. The formation of this lake is well documented. The glacier has retreated by some 600 m since 1981. - Since August 2002, a runoff gauge installed by the Hydrographical Service of Salzburg has recorded the discharge of the main part of SSK at the outlet of the new Unterer Eisboden See. The annual reports - submitted from 1982 on as a contractual obligation to the Hydrological Service of Salzburg - document the ongoing processes on the one hand, and emphasize the mass balance of SSK and outline the climatological reasons, mainly based on the met-data of the observatory Rudolfshütte, on the other. There is an additional focus on estimating the annual water balance in the catchment area of the lake. There are certain preconditions for the water balance equation in the area. Runoff is recorded by the ÖBB power stations, the mass balance of the now approx. 20% glaciated area (mainly the Sonnblickkees) is measured andthe change of the snow and firn patches/the water content is estimated as well as possible. (Nowadays laserscanning and ground radar are available to measure the snow pack). There is a net of three precipitation gauges plus the recordings at Rudolfshütte. The evaporation is of minor importance. The long-term annual mean runoff depth in the catchment area is around 3.000 mm/year. The precipitation gauges have measured deficits between 10% and 35%, on average probably 25% to 30%. That means that the real precipitation in the catchment area Weißsee (at elevations between 2,250 and 3,000 m) is in an order of 3,200 to 3,400 mm a year. The mass balance record of SSK was the first one established in the Hohe Tauern region (and now since the Hohe Tauern National Park was founded in 1983 in Salzburg) and is one of the longest measurement series worldwide. Great efforts are under way to continue the series, to safeguard against interruption and to guarantee a long-term monitoring of the mass balance and volume change of SSK (until the glacier is completely gone, which seems to be realistic in the near future as a result of the ongoing global warming). Heinz Slupetzky, March 2014

Nutrients and oxygen concentrations measured in water samples from the North Atlantic during the James Cook cruise JC087 in March 2013

Analysis for micro-molar concentrations of nitrate and nitrite, nitrite, phosphate, silicate and ammonia was undertaken on a SEAL Analytical UK Ltd, AA3 segmented flow autoanalyser following methods described by Kirkwood (1996). Samples were drawn from Niskin bottles on the CTD into 15ml polycarbonate centrifuge tubes and kept refrigerated at approximately 4oC until analysis, which generally commenced within 30 minutes. Overall 23 runs with 597 samples were analysed. This is a total of 502 CTD samples, 69 underway samples and 26 from other sources. An artificial seawater matrix (ASW) of 40g/litre sodium chloride was used as the inter-sample wash and standard matrix. The nutrient free status of this solution was checked by running Ocean Scientific International (OSI) low nutrient seawater (LNS) on every run. A single set of mixed standards were made up by diluting 5mM solutions made from weighed dried salts in 1litre of ASW into plastic 250ml volumetric flasks that had been cleaned by washing in MilliQ water (MQ). Data processing was undertaken using SEAL Analytical UK Ltd proprietary software (AACE 6.07) and was performed within a few hours of the run being finished. The sample time was 60 seconds and the wash time was 30 seconds. The lines were washed daily with wash solutions specific for each chemistry, but comprised of MQ, MQ and SDS, MQ and Triton-X, or MQ and Brij-35. Three times during the cruise the phosphate and silicate channels were washed with a weak sodium hypochlorite solution.

Experimental petrology of basement basaltic rocks from ODP Leg 127/128 samples

The relatively fresh basement basaltic rocks cored at Sites 794 and 797 during ODP Legs 127 and 128 show compositional variations suggesting the following: (1) the aphyric rocks might be differentiated from compositional equivalents of the aphyric sample with the lowest FeO*/MgO (Sample 127-797C-12R-4, 35-37 cm); and (2) the plagioclase-phyric rocks (i.e., another constituent of the basement basaltic rocks from the sites) may be derivatives from the same parents; in this case, however, crystallized plagioclase was not effectively removed. Melting experiments were conducted for Sample 127-797C-12R-4, 35-37 cm, and the differentiation processes for the basement basaltic rocks were assessed. The high-pressure melting-phase relation can not account for the compositional variation of the aphyric rocks, suggesting that the variation was developed at relatively low pressure where olivine and plagioclase fractionation was followed by Ca-rich clinopyroxene fractionation. The density of Sample 127-797C-12R-4,35-37 cm, is comparable to that of plagioclase at some depth, but at still relatively low pressure, making it possible that the liquidus plagioclase was retained in the successive liquids to produce the plagioclase-phyric rocks. According to backtrack calculation assuming the olivine maximum fractionation, Sample 127-797C-12R-4, 35-37 cm, was differentiated from primary picritic high-Al basalt magma. The estimated primary magma composition was experimentally proved to coexist with harzburgite mantle at about 14 kbar, suggesting relatively shallow production (approximately 40-50 km below surface) of the rifting-related primary magma.

Geologic-hydrologic setting and surface sedimentology in the Persian Gulf

1. Morphology and sedimentation The deepest parts of the Persian Gulf lie off the Iranian coast. Several swells separate the Persian Gulf into the Western Basin, the Central Basin and the Strait of Hormuz, which leads without noticeable morphological interruption onto the Biaban Shelf; the latter gradually drops off towards the continental slope, which itself has a strongly subdivided morphology. The sediment distribution in the Western Basin runs parallel to the basin's axis to a depth of 50 -60 m. This is caused by the shallow and uniform slope of the Iranian coast into the Western Basin, by clear exposure of the area to the Shamal-Winds and by tidal currents parallel to the basin's axis. Most other parameters also show isolines parallel to the coast line. Data from the sediment analyses show a net transport which extends out along the Central Swell: coarse fraction > 63 µ, total carbonate content, carbonate in fine fractions < 2 µ, 2-6 µ and 20-63 µ, calcite-aragonite ratios in the fine fractions 2-6 µ and 20-63 µ and quartz-dolomite ratios in fine fraction 2-6 µ. At least the uppermost 10-40 m of this sediment is late Holocene. This implies sedimentation rates of several meters per 1000 years. The slope from the Iranian coast into the Central Basin (max. depth 100 m) is generally steeper, with interspersed islands and flats. Both facts tend to disturb a sediment dustribition parallel to the basin's axis over extensive areas and may preclude any such trend from being detected by the methods and sample net used. The spatial distribution of the coarse fraction, however, seems to indicate sediment transport at greater water depths perpendicular to the basin's long axis and along the steepest gradients well into the Central Basin. The flats of the Central Basin have a sediment cover distinctly different from those of the deeper basin areas. Characteristic parameters are the extremely high percentages of coarse grained sediments, total content of carbonate CO2 over 40, low total organic carbon content, (however values are high if calculated on the basis of the < 63 µ fraction), low total N-content, and low C/N ratios. These characteristics probably result from the absence of any terrigenous material being brought in as well as from exposure to wave action. Finest terrigenous material is deposited in the innermost protected part of the Hormuz Bay. In the deep channel cut into the Biaban Shelf which carries the Persian Gulf out-flow water to the Indian Ocean, no fine grained sediment is deposited as shown by grain size data. 2. Geographic settings and sedimentation Flat lands border the Arabian coast of the Persian Gulf except for the Oman region. The high and steep Zagros Mountains form the Iranian coastline. Flat topography in combination with generally low precipitation precludes fluviatile sediment being added from the South. Inorganic and biogenic carbonates accumulating under low sedimentation rates are dominant on the shallow Arabic Shelf and the slopes into the Western and Central Basins. The fluviatile sediment brought in from the Iranian side, however decisively determine the composition of the Holocene sediment cover in the Persian Gulf and on the Biaban Shelf. Holocene sediments extend 20-30 km seaward into the Western Basin and about 25 km on to the Biaban Shelf. As mentioned before, sedimentation rates are of several meters/1000 years. The rocks exposed in the hinterland influence the sediments. According to our data the Redbeds of the Zagros Mountains determine the colour of the very fine grained sediments near the Iranian Coast of the Persian Gulf. To the West of Hormuz, addition of carbonate minerals is particularly high. Dolomite and protodolomite, deposited only in this area, as well as palygorskite, have proven to be excellent trace minerals. To the East of Hormuz, the supply of terrigenous carbonates is considerably lower. Clay minerals appear to bring in inorganically bound nitrogen thus lowering the C/N ratio in these sediments especially off river mouths. 3. Climate and sedimentation The Persian Gulf is located in a climatically arid region. This directly affects sedimentation through increased wind action and the infrequent but heavy rainfalls which cause flash floods. Such flash floods could be responsible for transporting sedheats into the Central Basin in a direction perpendicular to the Gulf's axis. Eolian influx is difficult to asses from our data; however, it probably is of minor importance from the Iranian side and may add, at the most, a few centimeters of fine sediment per 1000 years. 4. Hydrology and sedimentation High water temperatures favor inorganic carbonate precipitation in southern margin of the Gulf, and probably on the flats, as well as biogenic carbonate production in general. High evaporation plus low water inflow through rivers and precipitation cause a circulation pattern that is typical for epicontinental seas within the arid climate region. Surface water flows in from the adjoining ocean, in this case the Indian Ocean and sinks to the bottom of the Persian Gulf mainly in the northern part of the Western Basin, on the "Mesopotamischer Flachschelf" ard probably in the area of the "Arabischer Flachschelf". This sinking water continually rejuvenates the bottom out-flow water. The inflowing surface water from the Indian Ocean brings organic matter into the Persian Gulf, additional nutrients are added by the "fresh" upwelling waters of the Gulf of Oman. Both nutrients and organic matter diminish very rapidly as the water moves into the Persian Gulf. This depletion of nutrients and organic matter is the reasonfor generally low organic carbon contents of the Persian Gulf sediments. The Central Swell represents a distinct boundary, to the west of which the organic carbon content are lower than to the east when sediment samples of similar grain size distribution are compared. The outflow carries well oxygenated water over the bottom of the Persian Gulf and the resulting oxidation further decreases the content of organic matter. In the Masandam-Channel and in the Biaban-Shelf channel, the outflowing water prevents deposition of fine material and transports sediment particles well beyond the shelf margin. The outflowing water remains at a depth of 200-300 m depending on its density and releases ist suspending sediment load to the ocean floor, irrespectative of the bottom morphology. This is reflected in several parameters in which the sediments from beneath the outflow differ from nearby sediments not affected by the outflowing water. High carbonate content of total samples and of the individual size fraction as well as high aragonite and dolomite contents of individual size fractions characterize the sediment beneath the outflowing water. The tidal currents, which avt more or less parallel to the Gulf's axis, favor mixing of the water masses, they rework sediments at velocities reported here. This fact enlarges to a certain degree the extent of our interfaces which are based on only a few sample points (Persian Gulf and Biaban Shelf one sample per 620 km**2, continental slope one sample per 1000 km**2). The water on the continental slope shows and oxygen minimum at 200-1200 m which favors preservation of organically-bound carbon in the sediment. The low pH-values may even permit dissolution of carbonate minerals.