Landscape response to Pleistocene-Holocene precipitation change in the Western Cordillera, Peru: 10Be concentrations in modern sediments and terrace fills

The landscape response to climate change is frequently investigated with models because natural experiments on geologic timescales are rare. In Quebrada Veladera, in the western Andes Mountains, the formation of alluvial terraces during periods of high precipitation presents opportunities for such an experiment. We compare drainage-average erosion rates during Pleistocene terrace deposition with Holocene rates, using cosmogenic 10Be samples for seven pairs of quartz sand taken from the trunk and tributaries of Quebrada Veladera and adjacent terraces. Each pair consists of sediment collected from the modern channel and excavated from an adjacent fill terrace. The terrace fill was deposited at ~16 ka and preserved an isotopic record of paleoerosion rates in the Late Pleistocene. Modern sands yield 10Be concentrations between 1.68 × 105 and 2.28 × 105 atoms/g, corresponding to Holocene erosion rates between 43 ± 3 and 58 ± 4 mm/kyr. The 10Be concentrations in terrace sands range from 9.46 × 104 to 3.73 × 105 atoms/g, corresponding to paleoerosion rates from 27 ± 2 to 103 ± 8 mm/kyr. Smaller, upstream tributaries show a substantial decline in erosion rate following the transition from a wet to dry climate, but larger drainage areas show no change. We interpret this trend to indicate that the wetter climate drove landscape dissection, which ceased with the return to dry conditions. As channel heads propagated upslope, erosion accelerated in low-order drainages before higher-order ones. This contrast disappeared when the drainage network ceased to expand; at that point, erosion rates became spatially uniform, consistent with the uniformity of modern hillslope gradients. Key Points Landscape response to climate change evaluated with 10Be erosion rates Wetter climate associated with rapid erosion in smaller, upstream drainages Drier, Holocene climate associated with spatially uniform erosion rates ©2013. American Geophysical Union. All Rights Reserved

A late Holocene onset of Aboriginal burning in southeastern Australia

The extent to which Aboriginal Australians used fire to modify their environment has been debated for decades and is generally based on charcoal and pollen records rather than landscape responses to land-use change. Here we investigate the sensitivity of in-situ–produced 10Be, an isotope commonly used in geomorphological contexts, to anthropogenic perturbations in the southeastern Australian Tablelands. Comparing 10Be-derived erosion rates from fluvial sediment (8.7 ± 0.9 mm k.y.–1; 1 standard error, SE; n = 11) and rock outcrops (5.3 ± 1.4 mm k.y.–1; 1 SE; n = 6) confirms that landscape lowering rates integrating over 104–105 yr are consistent with rates previously derived from studies integrating over 104 to >107 yr. We then model an expected 10Be inventory in fluvial sediment if background erosion rates were perturbed by a low-intensity, high-frequency Aboriginal burning regime. When we run the model using the average erosion rate derived from 10Be in fluvial sediment (8.7 mm k.y.–1), measured and modeled 10Be concentrations overlap between ca. 3 ka and 1 ka. Our modeling is consistent with intensified Aboriginal use of fire in the late Holocene, a time when Aboriginal population growth is widely recognized

Long-term erosion rates of Panamanian drainage basins determined using in situ 10Be

Erosion rates of tropical landscapes are poorly known. Using measurements of in situ-produced 10Be in quartz extracted from river and landslide sediment samples, we calculate long-term erosion rates for many physiographic regions of Panama. We collected river sediment samples from a wide variety of watersheds (n = 35), and then quantified 24 landscape-scale variables (physiographic, climatic, seismic, geologic, and land-use proxies) for each watershed before determining the relationship between these variables and long-term erosion rates using linear regression, multiple regression, and analysis of variance (ANOVA). We also used grain-size-specific 10Be analysis to infer the effect of landslides on the concentration of 10Be in fluvial sediment and thus on erosion rates. Cosmogenic 10Be-inferred, background erosion rates in Panama range from 26 to 595 m My− 1, with an arithmetic average of 201 m My− 1, and an area-weighted average of 144 m My− 1. The strongest and most significant relationship in the dataset was between erosion rate and silicate weathering rate, the mass of material leaving the basin in solution. None of the topographic variables showed a significant relationship with erosion rate at the 95% significance level; we observed weak but significant correlation between erosion rates and several climatic variables related to precipitation and temperature. On average, erosion rates in Panama are higher than other cosmogenically-derived erosion rates in tropical climates including those from Puerto Rico, Madagascar, Australia and Sri Lanka, likely the result of Panama's active tectonic setting and thus high rates of seismicity and uplift. Contemporary sediment yield and cosmogenically-derived erosion rates for three of the rivers we studied are similar, suggesting that human activities are not increasing sediment yield above long-term erosion rate averages in Panama. 10Be concentration is inversely proportional to grain size in landslide and fluvial samples from Panama; finer grain sizes from landslide material have lower 10Be concentration than fine-grained fluvial sediment. Large grains from both landslide and stream sediments have similarly low 10Be concentrations. These data suggest that fluvial gravel is delivered to the channel by landslides whereas sand is preferentially delivered by soil creep and bank collapse. Furthermore, the difference in 10Be concentration in sand-sized material delivers by soil creep and that delivered by landsliding suggests that the frequency and intensity of landslides influence basin scale erosion rates

Atmospherically produced beryllium-10 in annually laminated late-glacial sediments of the North American Varve Chronology

We attempt to synchronize the North American Varve Chronology (NAVC) with ice core and calendar year timescales by comparing records of atmospherically produced 10Be fallout in the NAVC and in ice cores. The North American Varve Chronology (NAVC) is a sequence of 5659 varves deposited in a series of proglacial lakes adjacent to the southeast margin of the retreating Laurentide Ice Sheet between approximately 18 200 and 12 500 years before present. Because properties of NAVC varves are related to climate, the NAVC is also a climate proxy record with annual resolution, and our overall goal is to place the NAVC and ice core records on the same timescale to facilitate high-resolution correlation of climate proxy variations in both. Total 10Be concentrations in NAVC sediments are within the range of those observed in other lacustrine records of 10Be fallout, but 9Be and 10Be concentrations considered together show that the majority of 10Be is present in glacial sediment when it enters the lake, and only a minority of total 10Be derives from atmospheric fallout at the time of sediment deposition. Because of this, an initial experiment to determine whether or not 10Be fallout variations were recorded in NAVC sediments by attempting to observe the characteristic 11-year solar cycle in short varve sections sampled at high resolution was inconclusive: short-period variations at the expected magnitude of this cycle were not distinguishable from measurement scatter. On the other hand, longer varve sequences sampled at decadal resolution display centennial-period variations in reconstructed 10Be fallout that have similar properties as coeval 10Be fallout variations recorded in ice core records. These are most prominent in glacial sections of the NAVC that were deposited in proglacial lakes and are suppressed in paraglacial sections of the NAVC that were deposited in lakes lacking direct glacial sediment input. We attribute this difference to the fact that buffering of 10Be fallout by soil adsorption can filter out short-period variations in an entirely deglaciated watershed, but such buffering cannot occur in the ablation zone of an ice sheet. This implies that proglacial lakes whose watershed is mostly glacial may effectively record 10Be fallout variations. We attempted to match centennial-period variations in reconstructed 10Be fallout flux from two segments of the NAVC with ice core fallout records. For both records, it is possible to obtain matches that result in acceptable correlation between NAVC and ice core 10Be fallout records, but the best-fitting matches for the two segments disagree, and only one of them is consistent with independent calendar year calibrations of the NAVC and therefore potentially valid. This leaves several remaining ambiguities in whether or not 10Be fallout variations can, in fact, be used for synchronizing NAVC and ice core timescales, but these could most likely be resolved by higher-resolution and replicate 10Be measurements on targeted sections of the NAVC

Cosmogenic 26Al/10Be surface production ratio in Greenland

The assumed value for the cosmogenic 26Al/10Be surface production rate ratio in quartz is an important parameter for studies investigating the burial or subaerial erosion of long-lived surfaces and sediments. Recent models and data suggest that the production ratio is spatially variable and may be greater than originally thought. Here we present measured 26Al/10Be ratios for 24 continuously exposed bedrock and boulder surfaces spanning ~61–77°N in Greenland. Empirical measurements, such as ours, include nuclides produced predominately by neutron-induced spallation with percent-level contributions by muon interactions. The slope of a York regression line fit to our data is 7.3 ± 0.3 (1σ), suggesting that the 26Al/10Be surface production ratio exceeds the commonly used value of 6.75, at least in the Arctic. A higher 26Al/10Be production ratio has implications for multinuclide cosmogenic isotope studies because it results in greater modeled burial durations and erosion rates

Evidence for the stability of the West Antarctic Ice Sheet divide for 1.4 million years

Past fluctuations of the West Antarctic Ice Sheet (WAIS) are of fundamental interest because of the possibility of WAIS collapse in the future and a consequent rise in global sea level. However, the configuration and stability of the ice sheet during past interglacial periods remains uncertain. Here we present geomorphological evidence and multiple cosmogenic nuclide data from the southern Ellsworth Mountains to suggest that the divide of the WAIS has fluctuated only modestly in location and thickness for at least the last 1.4 million years. Fluctuations during glacial–interglacial cycles appear superimposed on a long-term trajectory of ice-surface lowering relative to the mountains. This implies that as a minimum, a regional ice sheet centred on the Ellsworth-Whitmore uplands may have survived Pleistocene warm periods. If so, it constrains the WAIS contribution to global sea level rise during interglacials to about 3.3 m above present