Postglacial to Holocene landscape evolution and process rates in steep alpine catchments

Climate change and high magnitude mass wasting events pose adverse societal effects and hazards, especially in alpine regions. Quantification of such geomorphic processes and their rates is therefore critical but is often hampered by the lack of appropriate techniques and the various spatiotemporal scales involved in these studies. Here we exploit both in situ cosmogenic beryllium-10 (¹⁰Be) and carbon-14 (¹⁴C) nuclide concentrations for deducing exposure ages and tracing of sediment through small alpine debris flow catchments in central Switzerland. The sediment cascade and modern processes we track from the source areas, through debris flow torrents to their final export out into sink regions with cosmogenic nuclides over an unprecedented five-year time series with seasonal resolution. Data from a seismic survey and a 90 m core revealed a glacially overdeepened basin, filled with glacial and paraglacial sediments. Surface exposure dating of fan boulders and radiocarbon ages constrain the valley fill from the last deglaciation until the Holocene and show that most of the fan existed in early Holocene times already. Current fan processes are controlled by episodic debris flow activity, snow (firn) and rock avalanches. Field investigations, digital elevation models (DEMs) of difference and geomorphic analysis agree with sediment fingerprinting with cosmogenic nuclides, highlighting that the bulk of material exported today at the outlet of the subcatchments derives from the lower fans. Cosmogenic nuclide concentrations steadily decrease from headwater sources to distal fan channels due to the incorporation of material with lower nuclide concentrations. Further downstream the admixture of sediment from catchments with less frequent debris flow activity can dilute the cosmogenic nuclide signals from debris flow dominated catchments but may also reach thresholds where buffering is limited. Consequently, careful assessment of boundary conditions and driving forces is required when apparent denudation rates derived from cosmogenic nuclide analysis are upscaled to larger regions. © 2018 John Wiley & Sons, Ltd.     

Erosion in northwest Tibet from in‐situ‐produced cosmogenic 10Be and 26Al in bedrock

Concentrations of in‐situ‐produced cosmogenic nuclides ¹⁰Be and ²⁶Al in quartz were measured by accelerator mass spectrometry for bedrock basalts and sandstones located in northwest Tibet. The effective exposure ages range between 23 and 134 ka (¹⁰Be) and erosion rates between 4·0 and 24 mm ka^?1. The erosion rates are significantly higher than those in similarly arid Antarctica and Australia, ranging between 0·1 and 1 mm ka^?1, suggesting that precipitation is not the major control of erosion of landforms. Comparison of erosion rates in arid regions with contrasting tectonic activities suggests that tectonic activity plays a more important role in controlling long‐term erosion rates. The obtained erosion rates are, however, significantly lower than the denudation rate of 3000–6000 mm ka^?1 beginning at c. 5‐3 Ma in the nearby Godwin Austen (K2) determined by apatite fission‐track thermochronology. It appears that the difference in erosion rates within different time intervals is indicative of increased tectonic activity at c. 5–3 Ma in northwest Tibet. We explain the low erosion rates determined in this study as reflecting reduced tectonic activity in the last million years. A model of localized thinning of the mantle beneath northwest Tibet may account for the sudden increased tectonic activity at c. 5–3 Ma and the later decrease. Copyright © 2006 John Wiley & Sons, Ltd.     

Estimation of episodic exfoliation rates of rock sheets on a granite dome in Korea from cosmogenic nuclide analysis

We have measured concentrations of cosmogenic ¹⁰Be and ²⁶Al produced in situ at bare bedrock surfaces of successive sheets developing on a granite dome in Korea and calculated the exfoliation rate of sheeting joints. The exfoliation rate was obtained using a simple model in which the sheeting joints experience intermittent denudation, i.e. peeling off along the bedrock face. We find that the average exfoliation (erosion) rate of the episodic peeling‐off process is 5·6 cm/ka^?1. The analysis is useful for understanding the evolution of granite sheeting structures on this dome in Korea. Copyright © 2006 John Wiley & Sons, Ltd.     

Direct observation of frost wedging in alpine bedrock

Width and temperature of rock joints were automatically monitored in the Japanese Alps. Three years of monitoring on a sandstone rock face shows two seasonal peaks of joint widening in autumn and spring. The autumn events are associated with short‐term freeze–thaw cycles, and the magnitude of widening reflects the freezing intensity and water availability. The short‐term freezing can produce wedging to a depth of at least 20 cm. The spring events follow a rise in the rock surface temperature to 0 °C beneath the seasonal snowcover, and likely originate from refreezing of meltwater entering the joint. Some of these events contribute to permanent enlargement of the joint. Two other joints on nearby rock faces experience only sporadic widening accompanying freeze–thaw cycles and insignificant permanent enlargement. Observations indicate that no single thermal criterion can explain frost weathering. The temperature range at which wedging occurs varies with the bedrock conditions, water availability and duration of freezing. Copyright © 2001 John Wiley & Sons, Ltd.     

Landscape preservation under Fennoscandian ice sheets determined from in situ produced Be and Al

Some areas within ice sheet boundaries retain pre-existing landforms and thus either remained as ice free islands (nunataks) during glaciation, or were preserved under ice. Differentiating between these alternatives has significant implications for paleoenvironment, ice sheet surface elevation, and ice volume reconstructions. In the northern Swedish mountains, in situ cosmogenic ¹⁰Be and ²⁶Al concentrations from glacial erratics on relict surfaces as well as glacially eroded bedrock adjacent to these surfaces, provide consistent last deglaciation exposure ages (∼8-13 kyr), confirming ice sheet overriding as opposed to ice free conditions. However, these ages contrast with exposure ages of 34-61 kyr on bedrock surfaces in these same relict areas, demonstrating that relict areas were preserved with little erosion through multiple glacial cycles. Based on the difference in radioactive decay between ²⁶Al and ¹⁰Be, the measured nuclide concentration in one of these bedrock surfaces suggests that it remained largely unmodified for a minimum period of 845₋₄₁₈⁺⁴⁶¹ kyr. These results indicate that relict areas need to be accounted for as frozen bed patches in basal boundary conditions for ice sheet models, and in landscape development models. Subglacial preservation also implies that source areas for glacial sediments in ocean cores are considerably smaller than the total area covered by ice sheets. These relict areas also have significance as potential long-term subglacial biologic refugia.     

The critical role of climate and saprolite weathering in landscape evolution

Landscapes evolve in response to external forces, such as tectonics and climate, that influence surface processes of erosion and weathering. Internal feedbacks between erosion and weathering also play an integral role in regulating the landscapes response. Our understanding of these internal and external feedbacks is limited to a handful of field‐based studies, only a few of which have explicitly examined saprolite weathering. Here, we report rates of erosion and weathering in saprolite and soil to quantify how climate influences denudation, by focusing on an elevation transect in the western Sierra Nevada Mountains, California. We use an adapted mass balance approach and couple soil‐production rates from the cosmogenic radionuclide (CRN) ¹⁰Be with zirconium concentrations in rock, saprolite and soil. Our approach includes deep saprolite weathering and suggests that previous studies may have underestimated denudation rates across similar landscapes. Along the studied climate gradient, chemical weathering rates peak at middle elevations (1200–2000 m), averaging 112·3 ± 9·7 t km^–2 y^–1 compared to high and low elevation sites (46·8 ± 5·2 t km^?2 y^?1). Measured weathering rates follow similar patterns with climate as those of predicted silica fluxes, modeled using an Arrhenius temperature relationship and a linear relationship between flux and precipitation. Furthermore, chemical weathering and erosion are tightly correlated across our sites, and physical erosion rates increase with both saprolite weathering rates and intensity. Unexpectedly, saprolite and soil weathering intensities are inversely related, such that more weathered saprolites are overlain by weakly weathered soils. These data quantify exciting links between climate, weathering and erosion, and together suggest that climate controls chemical weathering via temperature and moisture control on chemical reaction rates. Our results also suggest that saprolite weathering reduces bedrock coherence, leading to faster rates of soil transport that, in turn, decrease material residence times in the soil column and limit soil weathering. Copyright © 2009 John Wiley & Sons, Ltd.     

Modern sediment yield compared to geologic rates of sediment production in a semi‐arid basin,New Mexico: assessing the human impact

An Erratum has been published for this article in Earth Surface Processes and Landforms 29(13) 2004, 1707. In the semi‐arid Arroyo Chavez basin of New Mexico, a 2·28 km² sub‐basin of the Rio Puerco, we contrasted short‐term rates (3 years) of sediment yield measured with sediment traps and dams with long‐term, geologic rates (～10 000 years) of sediment production measured using ¹⁰Be. Examination of erosion rates at different time‐scales provides the opportunity to contrast the human impact on erosion with background or geologic rates of sediment production. Arroyo Chavez is grazed and we were interested in whether differences in erosion rates observed at the two time‐scales are due to grazing. The geologic rate of sediment production, 0·27 kg m^?2 a^?1 is similar to the modern sediment yields measured for geomorphic surfaces including colluvial slopes, gently sloping hillslopes, and the mesa top which ranged from 0·12 to 1·03 kg m^?2 a^?1. The differences between modern sediment yield and geologic rates of sediment production were most noticeable for the alluvial valley ?oor, which had modern sediment yields as high as 3·35 kg m^?2 a^?1. The hydraulic state of the arroyo determines whether the alluvial valley ?oor is aggrading or degrading. Arroyo Chavez is incised and the alluvial valley ?oor is gullied and piped and is a source of sediment. The alluvial valley ?oor is also the portion of the basin most modi?ed by human disturbance including grazing and gas pipeline activity, both of which serve to increase erosion rates. Copyright © 2004 John Wiley & Sons, Ltd.     

Study of the erosion rates in the upper Maracujá Basin (Quadrilátero Ferrífero/MG,Brazil) by the in situ produced cosmogenic 10Be method

The present work quantifies the erosive processes in the two main substrates (schists–phyllites and granites–gneisses) of the upper Maracujá Basin in the Quadrilátero Ferrífero/MG, Brazil, a region of semi‐humid tropical climate. Two measuring methods of concentration were used: (i) in situ produced ¹⁰Be in quartz veins (surface erosion rates) and (ii) ¹⁰Be in fluvial sediments (basin erosion rates). The results confirm that (i) erosion tends to be more aggressive close to the headwaters than in the lower parts of the basin and (ii) the region is now affected by dissection. Copyright © 2006 John Wiley & Sons, Ltd.     

Ice cap erosion patterns from bedrock 10Be and 26Al,southeastern Tibetan Plateau

Quantifying glacial erosion contributes to our understanding of landscape evolution and topographic relief production in high altitude and high latitude areas. Combining in situ ¹⁰Be and ²⁶Al analysis of bedrock, boulder, and river sand samples, geomorphological mapping, and field investigations, we examine glacial erosion patterns of former ice caps in the Shaluli Shan of the southeastern Tibetan Plateau. The general landform pattern shows a zonal pattern of landscape modification produced by ice caps of up to 4000 km² during pre-LGM (Last Glacial Maximum) glaciations, while the dating results and landforms on the plateau surface imply that the LGM ice cap further modified the scoured terrain into different zones. Modeled glacial erosion depth of 0–0.38 m per 100 ka bedrock sample located close to the western margin of the LGM ice cap, indicates limited erosion prior to LGM and Late Glacial moraine deposition. A strong erosion zone exists proximal to the LGM ice cap marginal zone, indicated by modeled glacial erosion depth >2.23 m per 100 ka from bedrock samples. Modeled glacial erosion depths of 0–1.77 m per 100 ka from samples collected along the edge of a central upland, confirm the presence of a zone of intermediate erosion in-between the central upland and the strong erosion zone. Significant nuclide inheritance in river sand samples from basins on the scoured plateau surface also indicate restricted glacial erosion during the last glaciation. Our study, for the first time, shows clear evidence for preservation of glacial landforms formed during previous glaciations under non-erosive ice on the Tibetan Plateau. As patterns of glacial erosion intensity are largely driven by the basal thermal regime, our results confirm earlier inferences from geomorphology for a concentric basal thermal pattern for the Haizishan ice cap during the LGM. © 2018 John Wiley & Sons, Ltd.     

Denudation rates on cratonic landscapes: comparison between suspended and dissolved fluxes,and 10Be analysis in the Nyong and Sanaga River basins,south Cameroon

South Cameroon is located in a tropical and tectonically quiescent region, with landscapes characterized by thick highly weathered regolith, indicative of the long‐term predominance of chemical weathering over erosion. Currently this region undergoes huge changes due to accelerated mutations related to a growing population and economical developments with associated needs and increasing pressures on land and natural resources. We analysed two of the main south Cameroon rivers: the Nyong River and Sanaga River. The Sanaga catchment undergoes a contrasted tropical climate from sub‐humid mountainous and humid climate and is impacted by deforestation, agriculture, damming, mining and urbanization, especially in the Mbam sub‐basin, draining the highly populated volcanic highlands. By contrast, the Nyong catchment, only under humid tropical climate, is preserved from anthropogenic disturbance with low population except in the region of Yaoundé (Méfou sub‐basin). Moreover the Nyong basin is dam‐free and less impacted by agriculture and logging. We explore both denudation temporal variability and the ratio between chemical and physical denudation through two catchment‐averaged erosion and denudation datasets. The first one consists of an 11‐year long gauging dataset, while the second one comes from cosmogenic radionuclides [CRNs, here beryllium‐10 (¹⁰Be)] from sand sampled in the river mainstreams (timescale of tens to hundreds of thousands of years). Modern fluxes estimated from gauging data range from 5 to 100 m/Ma (10 to 200 t/km²/yr); our calculations indicate that the usual relative contribution of chemical versus physical denudation is 60% and 40%, respectively, of the total denudation. Beryllium‐10 denudation rates and sediment fluxes range from 4.8 to 40.3 m/Ma or 13 to 109 t/km²/yr, respectively, after correction for quartz enrichment. These fluxes are slightly less than the modern fluxes observed in Cameroon and other stable tropical areas. The highest ¹⁰Be‐derived fluxes and the highest physical versus chemical denudation ratios are attributed to anthropogenic impact. Copyright © 2016 John Wiley & Sons, Ltd.     

Erosion processes and rates for gullies in granitic soils bathurst,new south wales,australia

Soil movement rates from fluted gully sides in New South Wales were measured over three years by erosion pins and survey. Erosion rates varied considerably over short distances from 0 to 53 mm yr^?1 with an average removal rate of 19-8 mm yr^?1.     

26Al and 10Be dating of late pleistocene and holocene fill terraces: a record of fluvial deposition and incision,Colorado front range

Cosmogenic ²⁶Al, ¹⁰Be, and ¹⁴C dating of fluvial fill terraces in steep canyons of the Colorado Front Range provides a temporal framework for analysing episodic aggradation and incision. Results from Boulder Canyon show that terrace heights above the modern channel (grade) can be divided into: (1) Bull Lake (≳100 ka; 20–15 m above grade); (2) Pinedale (32–10 ka; 15–4 m above grade); and (3) Holocene age (<4 m above grade). No pre‐Bull Lake deposits are preserved along Boulder Canyon, and only three small remnants >15 m above grade record Bull Lake deposition. Well‐preserved terraces of Pinedale age suggest that the range of terrace height above grade reflects short‐term fluctuations in the river profile during periods of rapidly changing stream load and power. Net river incision apparently occurred during transitions to interglacial periods. Soil development and stratigraphic position, along with limited cosmogenic and ¹⁴C dating, suggest that ∼130 ka terraces in Boulder Canyon correlate with the Louviers Alluvium, and that 32 to 10 ka fills in the canyon correlate with the Broadway Alluvium on the adjacent High Plains. Late Pleistocene incision rates (∼0·15 m ka⁻¹) along Boulder Canyon exceed pre‐late Pleistocene incision rates, and are higher than middle to late Pleistocene incision rates (∼0·04 m ka⁻¹) on the High Plains. This study provides an example of how modern geochronologic techniques allow us to understand better rivers that drain glaciated catchments. Copyright © 2002 John Wiley & Sons, Ltd.     

Plateau reduction by drainage divide migration in the Eastern Cordillera of Colombia defined by morphometry and 10Be terrestrial cosmogenic nuclides

Catchment‐wide erosion rates were defined using ¹⁰Be terrestrial cosmogenic nuclides for the Eastern Cordillera of the Colombian Andes to help determine the nature of drainage development and landscape evolution. The Eastern Cordillera, characterized by a smooth axial plateau bordered by steep flanks, has a mean erosion rate of 11 ± 1 mm/ka across the plateau and 70 ± 10 mm/ka on its flanks, with local high rates >400 mm/ka. The erosional contrast between the plateau and its flanks was produced by the increase in the orogen regional slope, derived from the progressive shortening and thickening of the Eastern Cordillera. The erosion rates together with digital topographic analysis show that the drainage network is dynamic and confirms the view that drainage divides in the Eastern Cordillera are migrating towards the interior of the mountain belt resulting in progressive drainage reorganization from longitudinal to transverse‐dominated rivers and areal reduction of the Sabana de Bogotá plateau. Copyright © 2016 John Wiley & Sons, Ltd.     

26Al/10Be ratios reveal the source of river sediments in the Kimberley,NW Australia

We use cosmogenic ¹⁰Be and ²⁶Al in both bedrock and fluvial sediments to investigate controls on erosion rates and sediment supply to river basins at the regional scale in the Kimberley, NW Australia. The area is characterised by lithologically controlled morphologies such as cuestas, isolated mesas and extensive plateaus made of slightly dipping, extensively jointed sandstones. All sampled bedrock surfaces at plateau tops, ridgelines, and in the broader floodplain of major rivers over the region show similar slow lowering rates between 0.17 and 4.88 m.Myr^-1, with a mean value of 1.0 ± 0.6 m.Myr^-1 (n=15), whilst two bedrock samples collected directly within river-beds record rates that are one to two orders of magnitude higher (14.4 ± 1.5 and 20.9 ± 2.5 m.Myr^-1, respectively). Bedrock ²⁶Al/¹⁰Be ratios are all compatible with simple, continuous sub-aerial exposure histories. Modern river sediment yield lower ¹⁰Be and ²⁶Al concentrations, apparent ¹⁰Be basin-wide denudation rates ranging between 1.8 and 7.7 m.Myr^-1, with a median value of 2.6 m.Myr^-1, more than double the magnitude of bedrock erosion rates. ²⁶Al/¹⁰Be ratios of the sediment samples are lower than those obtained for bedrock samples. We propose that these depleted ²⁶Al/¹⁰Be ratios can largely be explained by the supply of sediment to river basins from the slab fragmentation and chemical weathering of channel gorge walls and plateau escarpments that result in diluting the cosmogenic nuclide concentration in river sediments measured at the basin outlets. The results of a mass-balance model suggest that ~60–90% of river sediment in the Kimberley results from the breakdown and chemical weathering of retreating vertical sandstone rock-walls in contrast to sediment generated by bedrock weathering and erosion on the plateau tops. This study emphasises the value of analysing two or more isotopes in basin-scale studies using cosmogenic nuclides, especially in slowly eroding post-orogenic settings. © 2019 John Wiley & Sons, Ltd.     

Improved discrimination of subglacial and periglacial erosion using 10Be concentration measurements in subglacial and supraglacial sediment load of the Bossons glacier (Mont Blanc massif,France)

Deciphering the complex interplays between climate, uplift and erosion is not straightforward and estimating present‐day erosion rates can provide useful insights. Glaciers are thought to be powerful erosional agents, but most published ‘glacial’ erosion rates combine periglacial, subglacial and proglacial erosion processes. Within a glaciated catchment, sediments found in subglacial streams originate either from glacial erosion of substratum or from the rock walls above the glacier that contribute to the supraglacial load. Terrestrial cosmogenic nuclides (TCN) are produced by interactions between cosmic ray particles and element targets at the surface of the Earth, but their concentration becomes negligible under 15 m of ice. Measuring TCN concentrations in quartz sand sampled in subglacial streams and in supraglacial channels is statistically compliant with stochastic processes (e.g. rockfalls) and may be used to discriminate subglacial and periglacial erosion. Results for two subglacial streams of the Bossons glacier (Mont Blanc massif, France) show that the proportion of sediments originating from glacially eroded bedrock is not constant: it varies from 50% to 90% (n = 6). The difference between the two streams is probably linked to the presence or absence of supraglacial channels and sinkholes, which are common features of alpine glaciers. Therefore, most of the published mean catchment glacial erosion rates should not be directly interpreted as subglacial erosion rates. In the case of catchments with efficient periglacial erosion and particularly rockfalls, the proportion of sediments in the subglacial stream originating from the supraglacial load could be considerable and the subglacial erosion rate overestimated. Here, we estimate warm‐based subglacial and periglacial erosion rates to be of the same order of magnitude: 0.39 ± 0.33 and 0.29 ± 0.17 mm a^?1, respectively. Copyright © 2015 John Wiley & Sons, Ltd.     

Spatial and temporal replicability of meteoric and in situ 10Be concentrations in fluvial sediment

Cosmogenic isotopes, short‐lived radionuclides, elemental concentrations and thermochronometric indicators are measured in river sand to quantify erosion rates and trace sediment sources, and/or infer erosional processes. Interpretations of detrital sediment analyses are often based on the rarely tested assumption of time‐invariant tracer concentration. A better understanding of when and where this assumption breaks down and what sampling strategies minimize temporal and small‐scale spatial variance will improve science done using detrital river sediment. Here, we present new and previously published spatial and temporal replicates measured for in situ and meteoric ¹⁰Be (¹⁰Be_i and ¹⁰Be_m, respectively). Our new data include 113 replicate pairs, taken from agricultural and/or tectonically active watersheds in China months to millennia apart and spatial replicates taken up to 2 km apart on the same day. The mean percentage difference is 10% (?122% to 150%) for both systems considered together; the mode is close to 0% for both systems; and 36% of pairs of samples replicate within our analytical accuracy at 2σ . We find that ¹⁰Be_i replicates better than ¹⁰Be_m (p < 0.01). ¹⁰Be_i replicability is worse in steeper basins, suggesting that stochastic processes (i.e. landslides) affect reproducibility. ¹⁰Be_m replicability is worse in larger basins, suggesting non‐conservative behavior of ¹⁰Be_m as sediment moves downstream. Our results are consistent with the few previously published replicate studies. Considering all replicate data in a wide range of landscapes, in areas with deep erosional processes, replicability is poor; in other areas, replicability is good. This suggests that, in steep, tectonically active, and/or agricultural landscapes, individual detrital sediment measurements do not represent upstream rates as well as larger populations of samples. To ensure that measurements are representative of the upstream watershed, our data suggest that samples be amalgamated either over time or from several places close by in the same channel. Copyright © 2017 John Wiley & Sons, Ltd.     

Cosmogenic nuclide‐derived rates of diffusive and episodic erosion in the glacially sculpted upper Rhone Valley,Swiss Alps

Denudation rates of small tributary valleys in the upper Rhone valley of the Swiss Central Alps vary by more than an order of magnitude within a very small distance (tens of kilometers). Morphometric data indicate two distinct erosion processes operate in these steep mountain valleys. We determined the rates of these processes using cosmogenic beryllium‐10 (¹⁰Be) in pooled soil and stream sediment samples. Denudation in deep, glacially scoured valleys is characterized by rapid, non‐uniform processes, such as debris flows and rock falls. In these steep valleys denudation rates are 760–2100 mm kyr^?1. In those basins which show minimal previous glacial modification denudation rates are low with 60–560 mm kyr^?1. The denudation rate in each basin represents a binary mixture between the rapid, non‐uniform processes, and soil creep. The soil production rate measured with cosmogenic ¹⁰Be in soil samples averages at 60 mm kyr^?1. Mixing calculations suggest that the debris flows and rock falls are occurring at rates up to 3000–7000 mm kyr^?1. These very high rates occur in the absence of baselevel lowering, since the tributaries drain into the Rhone trunk stream up‐stream of a knickzone. The flux‐weighted spatial average of denudation rates for the upper Rhone valley is 1400 mm kyr^?1, which is similar to rock uplift rates determined in this area from leveling. The pace and location of erosion processes are determined by the oscillation between a glacial and a non‐glacial state, preventing the landscape from reaching equilibrium. Copyright © 2010 John Wiley & Sons, Ltd.     

Determination of predevelopment denudation rates of an agricultural watershed (Cayaguás River, Puerto Rico) using in-situ-produced Be in river-borne quartz

Accurate estimates of watershed denudation absent anthropogenic effects are required to develop strategies for mitigating accelerated physical erosion resulting from human activities, to model global geochemical cycles, and to examine interactions among climate, weathering, and uplift. We present a simple approach to estimate predevelopment denudation rates using in-situ-produced cosmogenic ¹⁰Be in fluvial sediments. Denudation processes in an agricultural watershed (Cayaguás River Basin, Puerto Rico) and a matched undisturbed watershed (Icacos River Basin) were compared using ¹⁰Be concentrations in quartz for various size fractions of bed material. The coarse fractions in both watersheds bear the imprint of long subsurface residence times. Fine material from old shallow soils contributes little, however, to the present-day sediment output of the Cayaguás. This confirms the recent and presumably anthropogenic origin of the modern high denudation rate in the Cayaguás Basin and suggests that pre-agricultural erosional conditions were comparable to those of the present-day Icacos.     

Temporal and spatial variability of cation and silica export in an alpine watershed,Emerald Lake,California

A reaction set of possible mineral weathering reactions is proposed to explain observed cation and silica export for the Emerald Lake watershed, a small Sierra Nevada, California catchment. The reaction set was calculated through a stoichiometric mole‐balance method, using a multiyear record of stream flow and snowpack chemical analyses and site‐specific mineral compositions. Reaction‐set calculations were intended to explore how the processes controlling stream cation and silica export depend on differing bedrock mineralogy across the catchment as snowmelt and runoff patterns change over the year. Different regions within the watershed can be differentiated by lake inflow subdrainages, each exhibiting different stream‐flow chemistry and calculated weathering stoichiometry, indicating that different silica and cation generation processes are dominant in wet steep portions of the catchment. Short‐term differences in stream concentrations were assumed to reflect ion exchange equilibria and rapid biological processes, whereas long‐term persistent stream concentration differences in different areas of the catchment were assumed to reflect spatial variability in mineral weathering stoichiometry. Mineralogical analyses of rock samples from the watershed provided site‐specific chemical compositions of major mineral species for reaction calculations. Reaction sets were evaluated by linear regression of calculated versus observed differences between snowmelt and stream‐flow chemistry and by a combined measure. Initially, single weathering reactions were balanced and evaluated to determine the reactions that best explained observed stream chemical export. Next, reactions were combined, using mineral compositions from different rock types to estimate the dependence of ion fluxes on lithology. The seasonal variability of major solute calculated fluxes is low, approximately one order of magnitude, relative to the observed three orders of magnitude variability in basin discharge. Reaction sets using basin‐averaged lithology and Aplite lithologies gave superior explanations of stream chemical composition. Copyright © 2004 John Wiley & Sons, Ltd.