Suzanne P. Anderson

Glaciers dominate eustatic sea-level rise in the 21st century

Ice loss to the sea currently accounts for virtually all of the sea-level rise that is not attributable to ocean warming, and about 60% of the ice loss is from glaciers and ice caps rather than from the two ice sheets. The contribution of these smaller glaciers has accelerated over the past decade, in part due to marked thinning and retreat of marine-terminating glaciers associated with a dynamic instability that is generally not considered in mass-balance and climate modeling. This acceleration of glacier melt may cause 0.1 to 0.25 meter of additional sea-level rise by 2100.

HYDROLOGIC RESPONSE OF A STEEP, UNCHANNELED VALLEY TO NATURAL AND APPLIED RAINFALL

Observations from natural rain storms and sprinkling experiments at a steep zero-order catchment in the Oregon Coast Range demonstrate the importance offlow through near-surface bedrock on runoff generation and pore pressure development in shallow colluvial soils. Sprinkling experiments, involving irrigation of the entire 860 m 2 catchment at average intensities of 1.5 and 3.0 mm/h, permitted detailed observation of runoff and the development of subsurface saturation under controlled conditions. A weir installed to collectflow through the colluvium at the base of the catchment recovered runoff equal to one third to one half of the precipitation rate during quasi-steady irrigation. Three key observations demonstrate that a significant proportion of storm runoffflows through near-surface bedrock and illustrate the importance of shallow bedrockflow in pore pressure development in the overlying colluvial soil: (1) greater discharge recovery during both the experiments and natural rainfall at a weir installed approximately 15 m downslope of the weir at the base of the catchment, (2) spatially discontinuous patterns of positive pressure head in the colluvium during steady sprinkling, and (3) local development of upward head gradients associated withflow from weathered rock into the overlying colluvium during high-intensity rainfall. Data from natural storms also show that smaller storms produce no significant runoff or piezometric response and point to a critical intensity-duration rainfall to overcome vadose zone storage. Together these observations highlight the role of interaction betweenflow in colluvium and near- surface bedrock in governing patterns of soil saturation, runoff production, and positive pore pressures.

Weathering profiles, mass-balance analysis, and rates of solute loss: Linkages between weathering and erosion in a small, steep catchment

In a headwater catchment in the Oregon Coast Range, we find that solid-phase mass losses due to chemical weathering are equivalent in the bedrock and the soil. However, the long-term rate of mass loss per unit volume of parent rock is greater in the soil than in the rock. We attribute this finding to the effects of biotic processes in the soil and to hydrologic conditions that maximize contact time and water flux through the mineral matrix in the soil. This result stems both from earlier work in which we demonstrated that rock and soil contribute equally to the solute flux and from arguments presented here that the basin is in dynamic equilibrium with respect to erosion and uplift. The silica flux of 10.7 ± 7.1 t·km−2·yr−1 from the basin is several times larger than the flux from older soils elsewhere, but comparable to the flux from sites with similar physical erosion rates. This result argues that physical denudation or uplift rates play an important role in setting the chemical denudation rate. Physical processes appear to influence chemical-weathering rates in several ways. First, they limit chemical evolution by removing material, thus setting the residence time within the weathered rock and the soil. Second, bioturbation mixes rock fragments into the more reactive soil and maintains high soil porosity, allowing free circulation of water. Because the weathering in the soil is more intense than in the rock, we argue that the chemical denudation rate will diminish where uplift rates—and, hence, physical-denudation rates—are great enough to lead to a bedrock-dominated landscape. Chemical denudation rates will increase with physical-denudation rates, but only as long as the landscape remains mantled by soil.

Chemical weathering in glacial environments

Do glaciers enhance or inhibit chemical weathering rates relative to other environments? The importance of glaciers in the global carbon cycle and climate change hinges on the answer. We show that catchments occupied by active alpine glaciers yield cation denudation rates greater than the global mean rate but do not exceed rates in nonglacial catchments with similar water discharge. Silica denudation rates are distinctly lower in glacier-covered catchments than in their nonglacial counterparts. Because sediment yields are high from glaciers, this suggests that water flux, rather than physical erosion, exerts the primary control on chemical erosion by glaciers. Potassium and calcium concentrations are high relative to other cations in glacial water, probably due to dissolution of soluble trace phases, such as carbonates, exposed by comminution, and cation leaching from biotite. Preferential weathering of biotite may result in higher 87 Sr/ 86 Sr in glacial runoff than expected from whole-rock compositions. Thus, although glaciers do not influence total chemical denudation rates at a given runoff, they may yield compositionally distinctive chemical fluxes to the oceans. Disruption of mineral lattices by grinding increases dissolution rates; this and high surface area should make glacial sediments exceptionally weatherable. Weathering of glacial erosion products in environments beyond the glacier margin deserves attention because it may figure prominently in global chemical cycles.

Unsaturated zone processes and the hydrologic response of a steep, unchanneled catchment

As part of a larger, collaborative study, we conducted field experiments to investigate how rainfall signals propagate through an unsaturated soil profile, leading to a rapid pore pressure response and slope instability. We sprinkler-irrigated an entire, unchanneled headwater basin in the steep, humid Oregon Coast Range, and we drove the system to quasi steady state as indicated by tensiometers, piezometers, and discharge. During initial wetting some of the deeper tensiometers responded before the arrival of an advancing head gradient front. With continued irrigation most tensiometers attained near- zero pressure heads before most piezometers responded fully, and a stable unsaturated flow field preceded the development of a stable saturated flow field. Steady discharge occurred after the last piezometer reached steady state. With the onset of steady discharge the unsaturated zone, saturated zone, and discharge became delicately linked, and a spike increase in rain intensity led to a response in the saturated zone and discharge much faster than could have happened through advection alone. We propose that the rain spike produced a slight pressure wave that traveled relatively rapidly through the unsaturated zone, where it caused a large change in hydraulic conductivity and the rapid effusion of stored soil water. An important control on the hydrologic response of this catchment lies with the soil-water retention curve. In general, below pressure heads of about 20.05 m, soil-water contents change slightly with changes in pressure head, but above 20.05 m the soil-water content is highly variable. Minor rainstorms upon a wet soil can produce slight changes in pressure head and corresponding large changes in soil-water content, giving rise to the passage of pressure waves in response to increased rain intensity and a relatively rapid response in the unsaturated zone. This rapid unsaturated zone response led to a rapid rise in the saturated zone, and it may be the underlying mechanism enabling short bursts of rain to cause slope instability.

Subsurface flow paths in a steep, unchanneled catchment

Tracer studies during catchment-scale sprinkler experiments illuminate the pathways of subsurface flow in a small, steep catchment in the Oregon Coast Range. Bromide point injections into saturated materials showed rapid flow in bedrock to the catchment outlet. Bedrock flow returned to the colluvium, sustaining shallow subsurface flow there. The bromide peak velocity of ;10 23 ms 21 exceeded the saturated hydraulic conductivity of intact bedrock. This, and the peak shapes, verify that fractures provide important avenues for saturated flow in the catchment. Deuterium added to the sprinkler water moved through the vadose zone as plug flow controlled by rainfall rate and water content. Ninety-two percent of the labeled water remained in the vadose zone after 3 days (;140 mm) of sprinkling. Preferential flow of new water was not observed during either low-intensity irrigation or natural storms; however, labeled preevent water was mobile in shallow colluvium during a storm following our spiking experiment. In response to rainfall, waters from the deeper bedrock pathway, which have traveled through the catchment, exfiltrate into the colluvium mantle and mix with relatively young vadose zone water, derived locally, creating an area of subsurface saturation near the channel head. This effectively becomes a subsurface variable source area, which, depending on its size and the delivery of water from the vadose zone, dictates the apportioning of old and new water in the runoff and, correspondingly, the runoff chemistry. The slow movement of water through the vadose zone allows for chemical modification and limits the amount of new water in the runoff. Moreover, it suggests that travel time of new rain water does not control the timing of runoff generation.

Microbial Community Succession in an Unvegetated, Recently Deglaciated Soil

Primary succession is a fundamental process in macroecosystems; however, if and how soil development influences microbial community structure is poorly understood. Thus, we investigated changes in the bacterial community along a chronosequence of three unvegetated, early successional soils (∼20-year age gradient) from a receding glacier in southeastern Peru using molecular phylogenetic techniques. We found that evenness, phylogenetic diversity, and the number of phylotypes were lowest in the youngest soils, increased in the intermediate aged soils, and plateaued in the oldest soils. This increase in diversity was commensurate with an increase in the number of sequences related to common soil bacteria in the older soils, including members of the divisions Acidobacteria, Bacteroidetes, and Verrucomicrobia. Sequences related to the Comamonadaceae clade of the Betaproteobacteria were dominant in the youngest soil, decreased in abundance in the intermediate age soil, and were not detected in the oldest soil. These sequences are closely related to culturable heterotrophs from rock and ice environments, suggesting that they originated from organisms living within or below the glacier. Sequences related to a variety of nitrogen (N)-fixing clades within the Cyanobacteria were abundant along the chronosequence, comprising 6–40% of phylotypes along the age gradient. Although there was no obvious change in the overall abundance of cyanobacterial sequences along the chronosequence, there was a dramatic shift in the abundance of specific cyanobacterial phylotypes, with the intermediate aged soils containing the greatest diversity of these sequences. Most soil biogeochemical characteristics showed little change along this ∼20-year soil age gradient; however, soil N pools significantly increased with soil age, perhaps as a result of the activity of the N-fixing Cyanobacteria. Our results suggest that, like macrobial communities, soil microbial communities are structured by substrate age, and that they, too, undergo predictable changes through time.

Chemical weathering in the foreland of a retreating glacier

Chemical denudation rates and strontium isotope ratios in streams vary substantially and systematically in the foreland of the retreating Bench Glacier in south-central Alaska. To study weathering of young glacier sediments, we sampled 12 streams draining a chronosequence of till and moraine soils derived from Cretaceous metagraywacke–metapelite bedrock. Both sediment age and vegetation cover increase with distance from the glacier. Cation denudation rates decline with increasing distance from the glacier, whereas silica denudation rates increase. Carbonate dissolution and sulfide oxidation account for roughly 90% of the solute flux from the youngest sediments. Biotite alteration accounts for 5–11% of the solute flux; its peak contribution is found in the glacier outlet stream. Silicate weathering is the dominant reaction only in the oldest sediments. In a laboratory dissolution experiment using fresh glacial sediment, carbonate dissolution dominated the solute flux during the first 700 hours, paralleling the behavior of young sediments in the field. In contrast to trends in the field, the silica flux did not increase after the carbonate was exhausted from the reactor. A possible reason for this difference is that establishment of vegetation causes an increase in silicate weathering. The 87Sr/86Sr ratio in the glacier outlet stream is greater than that in proglacial streams and in bulk rock, due to a greater contribution of biotite weathering in the outlet than in proglacial streams. Strontium isotope ratios decline with sediment exposure age in the proglacial streams, and are consistent with a carbonate source. Because the dominant weathering reactions in the young sediments are of carbonate and sulfide rather than silicate minerals, weathering at the glacier margin is not an important long-term sink for atmospheric CO2.

Comparison of Microbial Community Compositions of Two Subglacial Environments Reveals a Possible Role for Microbes in Chemical Weathering Processes

ABSTRACT Viable microbes have been detected beneath several geographically distant glaciers underlain by different lithologies, but comparisons of their microbial communities have not previously been made. This study compared the microbial community compositions of samples from two glaciers overlying differing bedrock. Bulk meltwater chemistry indicates that sulfide oxidation and carbonate dissolution account for 90% of the solute flux from Bench Glacier, Alaska, whereas gypsum/anhydrite and carbonate dissolution accounts for the majority of the flux from John Evans Glacier, Ellesmere Island, Nunavut, Canada. The microbial communities were examined using two techniques: clone libraries and dot blot hybridization of 16S rRNA genes. Two hundred twenty-seven clones containing amplified 16S rRNA genes were prepared from subglacial samples, and the gene sequences were analyzed phylogenetically. Although some phylogenetic groups, including the Betaproteobacteria, were abundant in clone libraries from both glaciers, other well-represented groups were found at only one glacier. Group-specific oligonucleotide probes were developed for two phylogenetic clusters that were of particular interest because of their abundance or inferred biochemical capabilities. These probes were used in quantitative dot blot hybridization assays with a range of samples from the two glaciers. In addition to shared phyla at both glaciers, each glacier also harbored a subglacial microbial population that correlated with the observed aqueous geochemistry. These results are consistent with the hypothesis that microbial activity is an important contributor to the solute flux from glaciers.

Numerical simulations of glacial-valley longitudinal profile evolution

Glaciers shape alpine landscapes. They broaden valley bottoms, enhance local valley relief, generate multiple steps, overdeepen valley floors, and cause tributary valleys to hang. These distinctive glacial signatures result from 10 4 –10 5 yr of erosion, during which swings in climate drive advances and retreats of alpine glaciers. We use a numerical model of glacial erosion to explore the development of the longitudinal profiles of glaciated valleys. The model is driven by the past 400 k.y. of variable climate. Because both sliding speed, which dictates abrasion rate, and water-pressure fluctuations, which strongly modulate quarrying rate, should peak at the equilibrium-line altitude (ELA), we expect the locus of most rapid erosion to follow the transient ELA. Simulations of a single glacial valley show rapid flattening of the longitudinal profile. Inclusion of a tributary glacier creates a step immediately downvalley of the tributary junction that persists over multiple glaciations and commonly leaves the tributary valley hanging. Steps and overdeepenings result from an increase in ice discharge immediately below the tributary junction, which is accommodated primarily by increased ice thickness and hence sliding rate. The size of the step increases with the ratio of tributary to trunk ice discharge, while the height of a hanging valley reflects the difference in the time-integrated ice discharge in tributary and trunk valleys and therefore increases as the discharge ratio decreases.