Quantifying the impact of former glaciation on catchment‐wide denudation rates derived from cosmogenic 10Be

Glacial denudation can significantly perturb terrestrial cosmogenic nuclide depth profiles and, if this is not corrected for, derived apparent denudation rates will overestimate the actual denudation rates. Here we determine how much ¹⁰Be‐derived denudation rates – calculated under the assumption of steady state – deviate from actual denudation rates as a function of three parameters: (1) the total amount of glacial denudation, (2) the post‐glacial denudation rate and (3) the time elapsed since deglaciation. We provide correction lines for the full parameter space explored (glacial denudation: 0.01–100 m; post‐glacial denudation rate: 1–1000 mm/ka; deglaciation: 1–100 ka before present), to evaluate and, if necessary, correct denudation rates for the impact of glacial denudation. Applied to ¹⁰Be‐derived catchment‐averaged denudation rates for formerly glaciated catchments in the Black Forest, Germany, we find that uncorrected denudation rates overestimate actual rates by up to a factor of three.     

Erosion rates on different timescales derived from cosmogenic 10Be and river loads: implications for landscape evolution in the Rhenish Massif, Germany

We determined erosion rates on timescales of 10¹–10⁴ years for two catchments in the northeastern Rhenish Massif, in order to unravel the Quaternary landscape evolution in a Variscan mountain range typical of central Europe. Spatially averaged erosion rates derived from in situ produced ¹⁰Be concentrations in stream sediment of the Aabach and M?hne watersheds range from 47 ± 6 to 65 ± 14 mm/ka and integrate over the last 9–13 ka. These erosion rates are similar to local rates of river incision and rock uplift in the Quaternary and to average denudation rates since the Mesozoic derived from fission track data. This suggests that rock uplift is balanced by denudation, i.e., the landscape is in a steady state. Short-term erosion rates were derived from suspended and dissolved river loads subsequent to (1) correcting for atmospheric and anthropogenic inputs, (2) establishing calibration curves that relate the amount of suspended load to discharge, and (3) estimating the amount of bedload. The resulting solid mass fluxes (suspended and bedload) agree with those derived from the sediment volume trapped in three reservoirs. However, resulting geogenic short-term erosion rates range from 9 to 25 mm/ka and are only about one-third of the rates derived from ¹⁰Be. Model simulations in combination with published sediment yield data suggest that this discrepancy is caused by at least three factors: (1) phases with higher precipitation and/or lower evapotranspiration, (2) rare flood events not captured in the short-term records, and (3) prolonged periods of climatic deterioration with increased erosion and sediment transport on hillslopes.     

Short‐time (<10 ka) denudation rates as a marker of active folding in the Zagros Fold Belt (Iran)

Coupling between tectonics and surface processes is usually ill‐quantified as other factors such as climate and lithology affect the later. We provide catchment‐wide ¹⁰Be denudation rates of the Mand catchment in the Zagros Fold Belt (Iran) to infer correlations between these rates and ongoing tectonic shortening in the region. Denudation rates are generally low (~0.05–0.1 mm/a) but increase to ~1 mm/a near the Halikan anticline, where changes in precipitation, lithology or hillslope gradient are insignificant. The denudation rates upstream and downstream of the Halikan anticline are consistent with the GPS convergence rates in these areas. The sharp increase in denudation rates over the Halikan anticline denotes its growth as previously detected from terrace incision. It also reveals small wavelength coupling between crustal deformation and erosion. Denudation rates are therefore a useful and sensitive tool that helps constraining non‐brittle active tectonics such as folding of a sedimentary cover.     

The Lastglacial and Holocene seismostratigraphy and sediment distribution of Lake Bolshoye Shchuchye,Polar Ural Mountains,Arctic Russia

Seismostratigraphical studies of the 11.8‐km²‐large and ~140‐m‐deep Lake Bolshoye Shchuchye, Polar Ural Mountains, reveal up to 160‐m‐thick acoustically laminated sediments in the lake basin. Using a dense grid of seismic lines, the spatial and temporal distributions of the sedimentary history have been reconstructed. Three regional seismic horizons have been identified and correlated with the well‐dated 24‐m‐long sediment core retrieved from the lake. Isopach maps constructed from the seismic data show four phases of sedimentation. A contour map of the deepest regional seismic reflector represents the earliest hemipelagic sedimentation in the lake. Three contour maps represent time intervals covering the last 23 cal. ka based on the well‐dated core stratigraphy from the lake. The detailed time constraints on the upper stratigraphical units in the lake allow calculation of the lake's development in terms of sediment fluxes and the denudation rates from the Last Glacial Maximum (LGM) to the present. The sedimentation in Lake Bolshoye Shchuchye has been dominated by hemipelagic processes during at least the last 24 cal. ka BP only locally interrupted by delta progradation and slope processes. A major shift in the sediment accumulation at c. 18.7 cal. ka BP is interpreted to mark the end of the local glacial maximum, greatly reduced denudation and the onset of the deglaciation period; this also demonstrates how fast the glaciers melted and possibly disappeared at the end of the LGM. The denudation rate during the Holocene is only a fifth of the LGM rate. The age of the oldest stratified sediments in Lake Bolshoye Shchuchye is not well constrained, but estimated as c. 50–60 ka.     

Late Quaternary denudation of the Alps,valley and lake fillings and modern river loads

Erosional denudation of the Alps and their role as sediment source underwent major changes throughout the Quaternary, by repeated glaciation and deglaciation. The sediment fluxes of 16 major Alpine drainage basins were quantified by determining the sediment volumes which have been trapped in valleys and lake basins. These became sedimentologically closed after the last glacier retreat around 17 000 cal. BP. The sediment volumes distributed over their provenance areas yield mean mechanical denudation rates between 250 to 1060 mm ka^–1. In contrast, modern denudation rates, derived from river loads and delta surveys, range from 30 to 360 mm ka^–1. Relief, such as mean elevation and slope, turned out to be the primary control of both modern and Late Glacial mechanical denudation. Rock types seem to be responsible for some scatter of the data, but their role is masked by other factors. Modern denudation rates increase with higher proportions of bare rocks and glaciated area, but decrease with forest cover. An area-weighted extrapolation of the studied drainage basins to the entire Alps on the basis of major morphotectonic zones yields a mean denudation rate of 620 mm ka^–1 over the last 17 000 years. This rate clearly exceeds the modern rate of 125 mm ka^–1. Lake sediments and palaeoclimatic reconstructions confirm that the sediment yield of the Alps reached a maximum during deglaciation when large masses of unconsolidated materials were available, vegetation was scarse, and transport capacities were high. During the early Holocene sediment yield declined to a minimum before some climate deterioration and human activities again accelerated erosional processes. Assuming a denudation rate in the early Holocene half of the modern one, the Late Glacial denudation rates must have been in the order of 1100 to 2900 mm ka^–1. Consequently, denudation rates during a glacial/interglacial cycle probably varied by a factor of 14, which lies well within the range of other studies in central Europe, Scandinavia and North America. From large scale sediment budgets of perialpine sedimentary basins the overall denudation rate of the Alps during the Quaternary has been c. 400 mm ka^–1, i.e. about one third lower than the estimate for the last 17000 years. This can be well explained by the outstanding role which deglaciation played in the time span studied here.     

Sediment production and transport from in situ-produced cosmogenic 10Be and river loads in the Napo River basin,an upper Amazon tributary of Ecuador and Peru

Cosmogenic nuclide-based denudation rates and published erosion rates from recent river gauging in the Napo River basin (Peruvian Amazonia) are used to decipher erosion and sedimentation processes along a 600 km long transect from the headwaters to the lowlands. The sediment-producing headwaters to the Napo floodplain are the volcanically active Ecuadorian Andes, which discharge sediment at a cosmogenic nuclide-based denudation rate of 0.49 ± 0.12 mm/yr. This denudation rate was calculated from an average ¹⁰Be nuclide concentration of 2.2 ± 0.5 × 10⁴ at/g_(Qz) that was measured in bedload-derived quartz. Within the Napo lowlands, a significant drop in trunk stream ¹⁰Be nuclide concentrations relative to the Andean hinterland is recorded, with an average concentration of 1.2 ± 0.5 × 10⁴ at/g_(Qz). This nuclide concentration represents a mixture between the ¹⁰Be nuclide concentration of eroded floodplain deposits, and that of sediment eroded from the Andean hinterland that is now carried in the trunk stream. Evidence for addition of sediment from the floodplain to the trunk stream is provided by published decadal-scale sediment flux measurements from gauging stations operated in the Napo basin, from which an increase from 12 × 10⁶ t/yr at the outflow of the Andes to ～47 × 10⁶ t/yr at the confluence with the Solimões (upper Amazon River) is recorded. Therefore, approximately 35 × 10⁶ t of floodplain sediment are added annually to the active Napo trunk stream. Combined with our nuclide concentration measurements, we can estimate that the eroded floodplain deposits yield a nuclide concentration of ～0.9 × 10⁴ at/g_(Qz) only. Under steady state surface erosion conditions, this nuclide concentration would translate to a denudation rate of the floodplain of ～0.47 mm/yr. However, we have no geomorphologic explanation for this high denudation rate within the low relief floodplain and thus suggest that this low-nuclide concentrated sediment is Andean-derived and would have been deposited in the floodplain at a time when erosion rates of the Andes were elevated. Therefore, the recently eroded floodplain sediment provides an Andean “paleo denudation rate” of 1.2 mm/yr that was calculated for high Andean production rates. A likely period for elevated erosion rates is the LGM, where climate and vegetation cover of the Andes differed from that of the Holocene. A possible cause for the erosion of the floodplain is the tectonic uplift of the Eastern Andes, which progressively shifts the Napo River northwards. Hence, the river cuts into ancient lowland sediment, which is admixed to the Andean sediment carried in the main Napo River.     

Carbonate drifts as marine archives of aeolian dust (Santaren Channel,Bahamas)

Sediment data from the Bahamian Santaren carbonate drift reveal the variability of trans‐Atlantic Saharan dust transport back to about 100 ka bp (Marine Isotope Stage 5·3) and demonstrate that carbonate drifts are a valuable pelagic archive of aeolian dust flux. Carbonate drift bodies are common around tropical carbonate platforms; they represent large‐scale accumulations of ocean‐current transported material, which originates from the adjacent shallow‐water carbonate factory as well as from pelagic production, i.e. periplatform ooze. Subordinately, there is a clay‐size to silt‐size non‐carbonate fraction, which typically amounts to less than 10% of the sediment volume and originates from aeolian and fluvial input. Sedimentation rates in the 5·42 m long core GeoHH‐M95‐524 recovered 25 km west of Great Bahama Bank in the Santaren Channel ranges from 1·5 to 24·5 cm ka^?1 with lowest values during the last glacial lowstand and highest values following platform re‐flooding around 8 ka bp . These sedimentation rates imply that carbonate drifts have not only the potential to resolve long‐term environmental changes on orbital timescales, but also millennial to centennial fluctuations during interglacials. The sediment core has been investigated with the aim of characterizing the lithogenic dust fraction. Laboratory analyses included X‐ray fluorescence core scanning, determination of carbonate content and grain‐size analyses (of bulk and terrigenous fraction), as well as visual inspections of the lithogenic residue; the age model is based on oxygen isotopes and radiocarbon ages. Data show that the input of aeolian dust in the periplatform ooze as indicated by Ti/Al and Fe/Al element ratios abruptly increases at 57 ka bp , stays elevated during glacial times, and reaches a Holocene minimum around 6·5 ka bp , contemporary to the African Humid Period. Subsequently, there is a gradual increase in dust flux which almost reaches glacial levels during the last centuries. Grain‐size data show that the majority of dust particles fall into the fine silt range (below 10 μm); however, there is a pronounced coarse dust fraction in the size range up to 63 μm and individual ‘giant’ dust particles are up to 515 μm in size. Total dust flux and the relative amounts of fine and coarse dust are decoupled. The time‐variable composition of the grain‐size spectrum is interpreted to reflect different dust transport mechanisms: fine dust particles are delivered by the trade winds and the geostrophic winds of the Saharan Air Layer, whereas coarse dust particles travel with convective storm systems. This mode of transport ensures continuous re‐suspension of large particles and results in a prolonged transport. In this context, grain‐size data from the terrigenous fraction of carbonate drifts provide a measure for past coarse dust transport, and consequently for the frequency of convective storm systems over the dust source areas and the tropical Atlantic.     

Nature and timing of large landslides in the Himalaya and Transhimalaya of northern India

Four large landslides, each with a debris volume >10⁶ m³, in the Himalaya and Transhimalaya of northern India were examined, mapped, and dated using ¹⁰Be terrestrial cosmogenic radionuclide surface exposure dating. The landslides date to 7.7±1.0 ka (Darcha), 7.9±0.8 ka (Patseo), 6.6±0.4 ka (Kelang Serai), and 8.5±0.5 ka (Chilam). Comparison of slip surface dips and physically reasonable angles of internal friction suggests that the landslides may have been triggered by increased pore water pressure, seismic shaking, or a combination of these two processes. However, the steepness of discontinuities in the Darcha rock-slope, suggests that it was more likely to have started as a consequence of gravitationally-induced buckling of planar slabs. Deglaciation of the region occurred more than 2000 years before the Darcha, Patseo, and Kelang Serai landslides; it is unlikely that glacial debuttressing was responsible for triggering the landslides. The four landslides, their causes, potential triggers and mechanisms, and their ages are compared to 12 previously dated large landslides in the region. Fourteen of the 16 dated landslides occurred during periods of intensified monsoons. Seismic shaking, however, cannot be ruled out as a mechanism for landslide initiation, because the Himalaya has experienced great earthquakes on centennial to millennial timescales. The average Holocene landscape lowering due to large landslides for the Lahul region, which contains the Darcha, Patseo, and Kelang Serai landslides, is ～0.12 mm/yr. Previously published large-landslide landscape-lowering rates for the Himalaya differ significantly. Furthermore, regional glacial and fluvial denudation rates for the Himalaya are more than an order of magnitude greater. This difference highlights the lack of large-landslide data, lack of chronology, problems associated with single catchment/large landslide-based calculations, and the need for regional landscape-lowering determinations over a standardized time period.     

Glacial retreat history of Nanhuta Shan (north-east Taiwan) from preserved glacial features: the cosmic ray exposure perspective.

Taiwan, located at the junction of the Pacific Ocean, the Eurasian continent, and the South China marginal Sea, is of particular interest for reconstructing paleoclimate periods in Eastern Asia. This study reports the first cosmic ray exposure dating (CRE) of glacial features in Taiwan. Among the areas where glacial relicts have been described in Taiwan, the Nanhuta Shan range is probably the place where glacial landforms are best preserved. We consequently focused on this area combining glacial geomorphology observations together with CRE dating using in situ produced ¹⁰Be of erratic boulders and ice-sculpted surfaces. When combined with the geomorphic characteristics of the sampled areas, the obtained minimum CRE ages suggest that the glacial retreat in the Nanhuta Shan commenced about 10±3 ka ago and retreat was complete by 7±1 ka ago. This is consistent with the Holocene warming trend deduced from other biological and physico-chemical paleoclimatic records for the region. Estimates of local bedrock surface denudation rates either directly from in situ produced ¹⁰Be measurements or from geomorphic considerations are employed to determine the preservation of such glacial features within the highly dynamic setting of Taiwan.     

Environmental controls on 10Be‐based catchment‐averaged denudation rates along the western margin of the Peruvian Andes

We present ¹⁰Be‐based basin‐averaged denudation rates for the entire western margin of the Peruvian Andes. Denudation rates range from c. 9 mm ka^?1 to 190 mm ka^?1 and are related neither to the subduction of the Nazca plate nor to the current seismicity along the Pacific coast and the occurrence of raised Quaternary marine terraces. Therefore, we exclude a tectonic control on denudation on a millennial time‐scale. Instead, we explain >60% of the observed denudation rates with a model where erosion rates increase either with mean basin slope angles or with mean annual water discharge. These relationships suggest a strong environmental control on denudation.     

The impact of diamond extraction on natural denudation rates in the Diamantina Plateau (Minas Gerais,Brazil)

The ¹⁰Be method was used to investigate the effect of mining activities on the natural denudation rates in alluvial sediments from catchments of the Southern Espinhaço Range (SER) in Minas Gerais State (Brazil). In this region, which is predominantly composed of quartzites, the ¹⁰Be concentrations were measured in alluvial sediments from catchments in a preserved natural area of the Serra do Cipó National Park and on the Diamantina Plateau, which was subjected to diamond extraction from beginnings of XVIII century until the end of the XX. Two types of drainage were identified in the Diamantina Plateau area: (i) reworked drainage (alluvial sediments reworked by panning) and (ii) overloaded drainage (alluvial sediments originating from panning processes on saprolites located upstream). The mean denudation value for the natural drainages (∼4.4 m.My-1) is similar to that of the reworked drainages (∼4.3 m.My-1) However, the denudation rates obtained for eleven samples from three sites in overloaded basins range from ∼6.4 m My-1 to ∼22.8 m My-1 and are thus higher than those determined for the reworked and natural basins. These results show that despite the alluvium deposits have been intensely reworked by panning, the values of denudation rates were not changed, they are similar to denudation rates from the natural drainages. However, the natural rates are lower than those affected by panning processes on saprolites.     

Solar influence and hydrological variability during the Holocene from a speleothem annual record (Molinos Cave,NE Spain)

We present a multi‐proxy approach to reconstructing Holocene climate conditions in northeastern Spain based on an excellent correlation among the lamina thickness, colour parameters and isotope (δ¹⁸O and δ¹³C) variations recorded in a speleothem. An age model constructed from five U/Th dates and annual lamina counting suggests that the uppermost 14.7 cm of the MO‐7 stalagmite grew between 7.2 and 2.5 ka before present but experienced a growth hiatus from 4.9 to 4.3 ka. Three spectral analysis methods were applied to 11 time series. The results reveal common solar periodicities on decennial (Gleissberg cycle) and centennial (De Vries‐Suess cycle) scales. The onset of Holocene carbonate precipitation in the MO‐7 stalagmite appears to be associated with a cold, wet period, whereas the hiatus and the end of growth are related to warm, dry periods. This environmental trend fits well within the regional Holocene climate.     

Climatic and eustatic signals in a global compilation of shallow marine carbonate accumulation rates

Two of the most important factors that control the accumulation rate of material in carbonate platform environments on geological time scales are climate and eustasy. Accurately assessing the importance of these inter‐related factors through the study of both modern and ancient carbonate facies, however, is problematic. These difficulties arise from both the complexities inherent in carbonate depositional systems and the demonstrable incompleteness of the stratigraphic record. Here, a new compilation of more than 19 000 global Phanerozoic shallow marine carbonate accumulation rates derived from nearly 300 individual stratigraphic sections is presented. These data provide the first global holistic view of changes in shallow marine carbonate production in response to climate and eustasy on geological time scales. Notably, a clear latitudinal dependence on carbonate accumulation rates is recognized in the data. Moreover, it can also be demonstrated that rates calculated across the last glacial maximum and Holocene track changes in sea‐level. In detail, the data show that globally averaged changes in carbonate accumulation rates lagged changes in sea‐level by ca 3 kyr, reflecting the commonly observed delay in the response of individual carbonate successions to sea‐level rise. Differences between the rates of carbonate accumulation and sea‐level change over the past 25 kyr ostensibly reflect changing accumulation mode, with platform drowning (give‐up mode) pervasive during peak Early Holocene sea‐level rise, followed by a switch to catch‐up mode accumulation from ca 9 ka to the present. Carbonate accumulation rates older than the Quaternary are typically calculated over time spans much greater than 100 kyr, and at these time spans, rates primarily reflect long‐term tectonically mediated accommodation space changes rather than shorter term changes in climate/eustasy. This finding, coupled with issues of stratigraphic incompleteness and data abundance, tempers the utility of this and other compilations for assessing accurately the role of climate and eustasy in mediating carbonate accumulation rates through geological time.     

Earliest Holocene deglaciation of the central Uummannaq Fjord system,West Greenland

Uummannaq Fjord, West Greenland, held the Uummannaq Ice Stream system that drained an estimated ~6% of the Greenland Ice Sheet (GrIS) during the Last Glacial Maximum. Published ages for the final deglaciation in Uummannaq Fjord vary from as early as c. 9.8 ka to as late as c. 5.3 ka. Assessing this variability requires additional chronological controls to improve the deglaciation history of central West Greenland. Here, we combine ¹⁴C dating of lake sediment cores with cosmogenic ¹⁰Be exposure dating at sites adjacent to the present GrIS margin in the central‐inland sector of the Uummannaq Fjord system. We find that ice retreated to or within the present GrIS margin at 10.8±0.2 ka (n = 6). Although this ‘final deglaciation’ to or within the present GrIS margin across the Uummannaq Fjord system varies from c. 10.8 to 5.3 ka, all chronologies indicate collapse from the continental shelf to the inner fjords at c. 11.0 ka, which occurred at a net retreat rate of 300–1100 m a⁻¹. The Uummannaq Fjord system deglaciated c. 1000 years earlier than the major fjord system to the south, Disko Bugt. However, similarly rapid retreat rates of the two palaeo‐ice stream systems suggest that their collapse may have been aided by high calving rates. The asynchronous deglaciation of the GrIS throughout the Uummannaq Fjord system probably relates to the influence of varying fjord geometry on marine glacier behaviour.     

Insights into Wasatch fault vertical slip rates using the age of sediments in Timpanogos Cave, Utah

Timpanogos Cave, located near the Wasatch fault, is about 357 m above the American Fork River. Fluvial cave sediments and an interbedded carbonate flowstone yield a paleomagnetic and U–Th depositional age of 350 to 780 ka. Fault vertical slip rates, inferred from calculated river downcutting rates, range between 1.02 and 0.46 mm yr^− 1. These slip rates are in the range of the 0–12 Ma Wasatch Range exhumation rate ( 0.5–0.7 mm yr^− 1), suggesting that the long-term vertical slip rate remained stable through mid-Pleistocene time. However, the late Pleistocene (0–250 ka) decelerated slip rate ( 0.2–0.3 mm yr^− 1) and the accelerated Holocene slip rate ( 1.2 mm yr^− 1) are consistent with episodic fault activity. Assuming that the late Pleistocene vertical slip rate represents an episodic slowing of fault movement and the long-term (0–12 Ma) average vertical slip rate, including the late Pleistocene and Holocene, should be  0.6 mm yr^− 1, there is a net late Pleistocene vertical slip deficit of  50–75 m. The Holocene and late Pleistocene slip rates may be typical for episodes of accelerated and slowed fault movement, respectively. The calculated late Pleistocene slip deficit may mean that the current accelerated Wasatch fault slip rate will extend well into the future.     

U‐Th Zircon Dating by Laser Ablation Single Collector Inductively Coupled Plasma‐Mass Spectrometry (LA‐ICP‐MS)

Zircon crystals in the age range of ca. 10–300 ka can be dated by ²³⁰Th/²³⁸U (U‐Th) disequilibrium methods because of the strong fractionation between Th and U during crystallisation of zircon from melts. Laser ablation inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) analysis of nine commonly used reference zircons (at secular equilibrium) and a synthetic zircon indicates that corrections for abundance sensitivity and dizirconium trioxide molecular ions (Zr₂O₃⁺) are critical for reliable determination of ²³⁰Th abundances in zircon. When corrected for abundance sensitivity and interferences, mean activity ratios of (²³⁰Th)/(²³⁸U) for nine reference zircons analysed on five different days averaged 0.995 ± 0.023 (95% confidence weighted by data‐point uncertainty only, MSWD = 1.6; n = 9), consistent with their U‐Pb ages > 4 Ma that imply equilibrium for all intermediate daughter isotopes (including ²³⁰Th) within the ²³⁸U decay chain. U‐Th zircon ages generated by LA‐ICP‐MS without mitigating (e.g., by high mass resolution) or correcting for abundance sensitivity and molecular interferences on ²³⁰Th are potentially unreliable. To validate the applicability of LA‐ICP‐MS to this dating method, we acquired data from three late Quaternary volcanic units: the 41 ka Campanian Ignimbrite (plutonic clasts), the 161 ka Kos Plateau Tuff (juvenile clasts) and the 12 ka Puy de Dôme trachyte lava (all eruption ages by Ar/Ar, with zircon U‐Th ages being of equal or slightly older). A comparison of the corrected LA‐ICP‐MS results with previously published secondary ion mass spectrometry (SIMS) data for these rocks shows comparable ages with equivalent precision for LA‐ICP‐MS and SIMS, but much shorter analysis durations (~ 2 min vs. ~ 15 min) per spot with LA‐ICP‐MS and much simpler sample preparation. Previously undated zircons from the Yali eruption (Kos‐Nisyros volcanic centre, Greece) were analysed using this method. This yielded a large age spread (~ 45 to > 300 ka), suggesting significant antecryst recycling. The youngest zircon age (~ 45 ± 10 ka) provides a reasonable maximum estimate for the eruption age, in agreement with the previously published age using oxygen isotope stratigraphy (~ 31 ka).     

Gypsum karst intensification as a consequence of sulphur mining activity (Jaziv field,Western Ukraine)

The scope of this study is to evaluate the parameters of the gypsum–anhydrite strata karstification under the influence of the quarry exploitation of Jaziv sulphur field (West Ukraine) accompanying by drainage. The quarry drainage provoked the enormous depression cone forming with 100 km² in area and the corresponding intensification of chemical denudation of the gypsum–anhydrite strata. The gypsum karstification rate for the 29-year period of the quarry drainage was 0.06% with the dissolved rock volume of 1,328,507 m³ that is about 80 times higher than in natural conditions. For the drainage period (29 years) the karst denudation rate was 17,952 m³/km² against the natural background of 231.3 m³/km². The absolute value of artificial denudation for the studied massif area is 1.79 cm for 29 years or 0.062 mm/year while the denudation rate under natural conditions would be 0.231 cm for 29 years or 0.0008 mm/year. The forecasted volumes of the surface-evident collapses are evaluated using the relationship between the calculated amount of dissolved sulphate rocks and volumes of the current surface-evident collapses. According to this evaluation, the current collapses correspond to about 34% only from the dissolved rocks volume calculated for the drainage period.     

Retreat of the Smith Sound Ice Stream in the Early Holocene

Nares Strait, a major connection between the Arctic Ocean and Baffin Bay, was blocked by coalescent Innuitian and Greenland ice sheets during the last glaciation. This paper focuses on the events and processes leading to the opening of the strait and the environmental response to establishment of the Arctic‐Atlantic throughflow. The study is based on sedimentological, mineralogical and foraminiferal analyses of radiocarbon‐dated cores 2001LSSL‐0014PC and TC from northern Baffin Bay. Radiocarbon dates on benthic foraminifera were calibrated with ΔR = 220±20 years. Basal compact pebbly mud is interpreted as a subglacial deposit formed by glacial overriding of unconsolidated marine sediments. It is overlain by ice‐proximal (red/grey laminated, ice‐proximal glaciomarine unit barren of foraminifera and containing >2 mm clasts interpreted as ice‐rafted debris) to ice‐distal (calcareous, grey pebbly mud with foraminifera indicative of a stratified water column with chilled Atlantic Water fauna and species associated with perennial and then seasonal sea ice cover) glacial marine sediment units. The age model indicates ice retreat into Smith Sound as early as c. 11.7 and as late as c. 11.2 cal. ka BP followed by progressively more distal glaciomarine conditions as the ice margin retreated toward the Kennedy Channel. We hypothesize that a distinct IRD layer deposited between 9.3 and 9 (9.4–8.9 1σ) cal. ka BP marks the break‐up of ice in Kennedy Channel resulting in the opening of Nares Strait as an Arctic‐Atlantic throughflow. Overlying foraminiferal assemblages indicate enhanced marine productivity consistent with entry of nutrient‐rich Arctic Surface Water. A pronounced rise in agglutinated foraminifers and sand‐sized diatoms, and loss of detrital calcite characterize the uppermost bioturbated mud, which was deposited after 4.8 (3.67–5.55 1σ) cal. ka BP. The timing of the transition is poorly resolved as it coincides with the slow sedimentation rates that ensued after the ice margins retreated onto land.     

Deglaciation and landscape history around Annapurna,Nepal, based on 10Be surface exposure dating

The High Himalaya is a key area for tectonic, geomorphological and climate studies, because of its extreme relief and location at the transition zone between areas with abundant monsoonal precipitation and the arid/semiarid Tibetan Plateau. We present ¹⁰Be surface exposure ages on 22 boulders from the Annapurna area in Nepal. The ages improve understanding of the Late Quaternary landscape history and the geomorphological processes operating in this part of the Himalaya.Although our study is reconnaissance in nature, it highlights the importance of catastrophic events, such as landslides and debris flows, for denudation of high mountains. Holocene exposure ages for the Dhampu–Chooya landslide (～4.1 ± 0.6 ka) and for 600 m of alluviation in Kali Gandaki Valley (～2.1 ± 0.6 ka), for example, indicate the frequent occurrence and extent of catastrophic events and their implications for natural hazards. We also offer an explanation for the differences in Late Quaternary glacial chronologies at closely spaced study sites in the Nepal Himalaya. Topographically controlled and spatially variable precipitation in the Himalaya determines the sensitivity of glaciers to changes in temperature and precipitation. Accordingly, some glaciers advanced in-phase with Northern Hemisphere ice sheets, whereas others reached their maximum extent at times of increased monsoonal precipitation during Marine Isotope Stage 3 and the early Holocene.     

The last Eurasian ice sheets – a chronological database and time‐slice reconstruction,DATED‐1

We present a new time‐slice reconstruction of the Eurasian ice sheets (British–Irish, Svalbard–Barents–Kara Seas and Scandinavian) documenting the spatial evolution of these interconnected ice sheets every 1000 years from 25 to 10 ka, and at four selected time periods back to 40 ka. The time‐slice maps of ice‐sheet extent are based on a new Geographical Information System (GIS) database, where we have collected published numerical dates constraining the timing of ice‐sheet advance and retreat, and additionally geomorphological and geological evidence contained within the existing literature. We integrate all uncertainty estimates into three ice‐margin lines for each time‐slice; a most‐credible line, derived from our assessment of all available evidence, with bounding maximum and minimum limits allowed by existing data. This approach was motivated by the demands of glaciological, isostatic and climate modelling and to clearly display limitations in knowledge. The timing of advance and retreat were both remarkably spatially variable across the ice‐sheet area. According to our compilation the westernmost limit along the British–Irish and Norwegian continental shelf was reached up to 7000 years earlier (at c. 27–26 ka) than the eastern limit on the Russian Plain (at c. 20–19 ka). The Eurasian ice sheet complex as a whole attained its maximum extent (5.5 Mkm²) and volume (~24 m Sea Level Equivalent) at c. 21 ka. Our continental‐scale approach highlights instances of conflicting evidence and gaps in the ice‐sheet chronology where uncertainties remain large and should be a focus for future research. Largest uncertainties coincide with locations presently below sea level and where contradicting evidence exists. This first version of the database and time‐slices (DATED‐1) has a census date of 1 January 2013 and both are available to download via the Bjerknes Climate Data Centre and PANGAEA ( www.bcdc.no ; http://doi.pangaea.de/10.1594/PANGAEA.848117 ).