Topographic correlations with soil and regolith thickness from shallow‐seismic refraction constraints across upland hillslopes in the Valles Caldera,New Mexico

How rock is weathered physically and chemically into transportable material is a fundamental question in critical‐zone science. In addition, the distribution of weathered material (soil and intact regolith) across upland landscapes exerts a first‐order control on the hydrology of watersheds. In this paper we present the results of six shallow seismic‐refraction surveys in the Redondo Mountain region of the Valles Caldera, New Mexico. The P‐wave velocities corresponding to soil (≤ 0.6 km s^?1) were inferred from a seventh seismic survey where soil‐thickness data were determined by pit excavation. Using multivariable regression, we quantified the relationships among slope gradient, aspect, and topographic wetness index (TWI) on soil and regolith (soil plus intact regolith) thicknesses. Our results show that both soil and regolith thicknesses vary inversely with TWI in all six survey areas while varying directly with slope aspect (i.e. thicker beneath north‐facing slopes) and inversely with slope gradient (i.e. thinner beneath steep slopes) in the majority of the survey areas. An empirical model based on power‐law relationships between regolith thickness and its correlative variables can fit our inferred thicknesses with R² ‐values up to 0.880 for soil and 0.831 for regolith in areas with significant topographic variations. These results further demonstrate the efficacy of shallow seismic refraction for mapping and determining how soil and regolith variations correlate with topography across upland landscapes. Copyright © 2016 John Wiley & Sons, Ltd.     

Local topography and erosion rate control regolith thickness along a ridgeline in the Sierra Nevada,California

The ridgelines of mountain ranges are a source of geomorphic information unadulterated by the arrival of sediment from upslope. Studies along ridgecrests, therefore, can help identify and isolate the controls on important regolith properties such as thickness and texture. A 1.5 km section of ridgeline in the Sierra Nevada (CA) with a tenfold decrease in erosion rate (inferred from ridgetop convexity) provided an opportunity to conduct a high‐resolution survey of regolith properties and investigate their controls. We found that regolith along the most quickly eroding section of the ridge was the rockiest and had the lowest clay concentrations. Furthermore, a general increase in regolith thickness with a slowing of erosion rate was accompanied by an increase in biomass, changes in vegetation community, broader ridgeline profiles, and an apparent increase in total available moisture. The greatest source of variation in regolith thickness at the 10–100 m scale, however, was the local topography along the ridgeline, with the deepest regolith in the saddles and the thinnest on the knobs. Because regolith in the saddles had higher surface soil moisture than the knobs, we conclude that the hydrological conditions primarily driven by local topography (i.e. rapid vs. slow drainage and water‐storage potential) provide the fundamental controls on regolith thickness through feedbacks incorporating physical weathering by the biota and chemical weathering. Moreover, because the ridgeline saddles are the uppermost extensions of first‐order valleys, we propose that the fluvial network affects regolith properties in the furthest reaches of the watershed. Copyright © 2015 John Wiley & Sons, Ltd.     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 1. Hydraulic theory of reactive transport

This is the first of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory (also known as Dupuit or Boussinesq theory). Here, we examine the effect of lateral flow on the downward fluxes of water and solutes through perched groundwater at steady state. We derive analytical expressions describing the decline in the downward flux rate with depth. Using these, we obtain analytical expressions for water age in a number of cases. The results show that when the permeability field is homogeneous, the spatial structure of water age depends qualitatively on a single dimensionless number, Hi. This number captures the relative contributions to the lateral hydraulic potential gradient of the relief of the lower‐most impermeable boundary (which may be below the weathering front within permeable or incipiently weathered bedrock) and the water table. A “scaled lateral symmetry” exists when Hi is low: age varies primarily in the vertical dimension, and variations in the horizontal dimension x almost disappear when the vertical dimension z is expressed as a fraction z/H(x) of the laterally flowing system thickness H(x). Taking advantage of this symmetry, we show how the lateral dimension of the advection–diffusion‐reaction equation can be collapsed, yielding a 1‐D vertical equation in which the advective flux downward declines with depth. The equation holds even when the permeability field is not homogeneous, as long as the variations in permeability have the same scaled lateral symmetry structure. This new 1‐D approximation is used in the accompanying paper to extend chemical weathering models derived for 1‐D columns to hillslope domains.     

Rock damage and regolith transport by frost: an example of climate modulation of the geomorphology of the critical zone

In this article we craft process‐specific algorithms that capture climate control of hillslope evolution in order to elucidate the legacy of past climate on present critical zone architecture and topography. Models of hillslope evolution traditionally comprise rock detachment into the mobile layer, mobile regolith transport, and a channel incision or aggradation boundary condition. We extend this system into the deep critical zone by considering a weathering damage zone below the mobile regolith in which rock strength is diminished; the degree of damage conditions the rate of mobile regolith production. We first discuss generic damage profiles in which appropriate length and damage scales govern profile shapes, and examine their dependence upon exhumation rate. We then introduce climate control through the example of rock damage by frost‐generated crack growth. We augment existing frost cracking models by incorporating damage rate limitations for long transport distances for water to the freezing front. Finally we link the frost cracking damage model, a mobile regolith production rule in which rock entrainment is conditioned by the damage state of the rock, and a frost creep transport model, to examine the evolution of an interfluve under oscillating climate. Aspect‐related differences in mean annual surface temperatures result in differences in bedrock damage rate and mobile regolith transport efficiency, which in turn lead to asymmetries in critical zone architecture and hillslope form (divide migration). In a quasi‐steady state hillslope, the lowering rate is uniform, and the damage profile is better developed on north‐facing slopes where the frost damage process is most intense. Because the residence times of mobile regolith and weathered bedrock in such landscapes are on the order of 10 to 100 ka, climate cycles over similar timescales result in modulation of transport and damage efficiencies. These lead to temporal variation in mobile regolith thickness, and to corresponding changes in sediment delivery to bounding streams. Copyright © 2012 John Wiley & Sons, Ltd.     

On the relationship between bedrock lowering and regolith thickness

A general equation is proposed for the relationship between bedrock lowering and regolith thickness which includes, as special cases, hypotheses by Ahnert, Armstrong, Culling, Kirkby and Young, at least approximately.     

Strong variation in weathering of layered rock maintains hillslope‐scale strength under high precipitation

The evolution of volcanic landscapes and their landslide potential are both dependent upon the weathering of layered volcanic rock sequences. We characterize critical zone structure using shallow seismic V_p and V_s profiles and vertical exposures of rock across a basaltic climosequence on Kohala peninsula, Hawai’i, and exploit the dramatic gradient in mean annual precipitation (MAP) across the peninsula as a proxy for weathering intensity. Seismic velocity increases rapidly with depth and the velocity–depth gradient is uniform across three sites with 500–600 mm/yr MAP, where the transition to unaltered bedrock occurs at a depth of 4 to 10 m. In contrast, velocity increases with depth less rapidly at wetter sites, but this gradient remains constant across increasing MAP from 1000 to 3000 mm/yr and the transition to unaltered bedrock is near the maximum depth of investigation (15–25 m). In detail, the profiles of seismic velocity and of weathering at wet sites are nowhere monotonic functions of depth. The uniform average velocity gradient and the greater depths of low velocities may be explained by the averaging of velocities over intercalated highly weathered sites with less weathered layers at sites where MAP > 1000 mm/yr. Hence, the main effect of climate is not the progressive deepening of a near‐surface altered layer, but rather the rapid weathering of high permeability zones within rock subjected to precipitation greater than ~1000 mm/yr. Although weathering suggests mechanical weakening, the nearly horizontal orientation of alternating weathered and unweathered horizons with respect to topography also plays a role in the slope stability of these heterogeneous rock masses. We speculate that where steep, rapidly evolving hillslopes exist, the sub‐horizontal orientation of weak/strong horizons allows such sites to remain nearly as strong as their less weathered counterparts at drier sites, as is exemplified by the 50°–60° slopes maintained in the amphitheater canyons on the northwest flank of the island. Copyright © 2017 John Wiley & Sons, Ltd.     

Understanding fracture distribution and its relation to knickpoint evolution in the Rio Icacos watershed (Luquillo Critical Zone Observatory,Puerto Rico) using landscape-scale hydrogeophysics

The Rio Icacos watershed in the Luquillo Mountains (Puerto Rico) is unique due to its extremely rapid weathering rates. The watershed is incised into a quartz diorite that has developed a large knickzone defining the river profile. Regolith thickness within the watershed generally decreases from 20 to 30 m at the ridges to several meters in the quartz diorite-dominated valley to tens of centimeters near the major river knickpoint, as determined from previous studies. Above the knickzone, we observe spheroidal corestones, but below this weathering is much less apparent. Measured erosion rates from previous studies are also high in the knickzone compared with upper elevations within the river profile. A suite of near-surface geophysical methods (i.e. ground penetrating radar and terrain conductivity) capable of fast data acquisition in rugged landscapes, was deployed at kilometer scales to characterize critical zone structure. Concentrations of chaotic ground penetrating radar (GPR) reflections and diffraction hyperbolas with low electrical conductivity were observed in vertical zones that outcrop at the land surface as areas of intense fracturing and spheroidally weathered corestones. The width of these fractured and weathered zones showed an increase with proximity to the knickpoint, and was attributed to dilation of these sub-vertical fractures near the knickpoint, as postulated theoretically by a stress model calculated for the topographic variability across the knickzone in the Rio Icacos, and that shows a release of compressive stress near the knickpoint. We hypothesize that erosion rates increase in the knickzone because of this inferred dilation of fractures. Specifically, opened fractures could enhance access of water and in turn promote spalling, erosion, and spheroidal weathering. This study shows that ground-based hydrogeophysical methods used at the landscape-scale (traditionally applied at smaller scales) can be used to explore critical zone architecture at the scales needed to explain the extreme variability in erosion rates across river profiles. © 2018 John Wiley & Sons, Ltd.     

Probing the deep critical zone beneath the Luquillo Experimental Forest,Puerto Rico

Recent work has suggested that weathering processes occurring in the subsurface produce the majority of silicate weathering products discharged to the world's oceans, thereby exerting a primary control on global temperature via the well‐known positive feedback between silicate weathering and CO₂. In addition, chemical and physical weathering processes deep within the critical zone create aquifers and control groundwater chemistry, watershed geometry and regolith formation rates. Despite this, most weathering studies are restricted to the shallow critical zone (e.g. soils, outcrops). Here we investigate the chemical weathering, fracturing and geomorphology of the deep critical zone in the Bisley watershed in the Luquillo Critical Zone Observatory, Puerto Rico, from two boreholes drilled to 37.2 and 27.0 m depth, from which continuous core samples were taken. Corestones exposed aboveground were also sampled. Weathered rinds developed on exposed corestones and along fracture surfaces on subsurface rocks slough off of exposed corestones once rinds attain a thickness up to ~1 cm, preventing the corestones from rounding due to diffusion limitation. Such corestones at the land surface are assumed to be what remains after exhumation of similar, fractured bedrock pieces that were observed in the drilled cores between thick layers of regolith. Some of these subsurface corestones are massive and others are highly fractured, whereas aboveground corestones are generally massive with little to no apparent fracturing. Subsurface corestones are larger and less fractured in the borehole drilled on a road where it crosses a ridge compared with the borehole drilled where the road crosses the stream channel. Both borehole profiles indicate that the weathering zone extends to well below the stream channel in this upland catchment; hence weathering depth is not controlled by the stream level within the catchment and not all of the water in the watershed is discharged to the stream. Copyright © 2013 John Wiley & Sons, Ltd.     

Aspect control of water movement on hillslopes near the rain–snow transition of the Colorado Front Range

In the Colorado Front Range, forested catchments near the rain–snow transition are likely to experience changes in snowmelt delivery and subsurface water transport with climate warming and associated shifts in precipitation patterns. Snowpack dynamics are strongly affected by aspect: Lodgepole pine forested north‐facing slopes develop a seasonal snowpack, whereas Ponderosa pine‐dotted south‐facing slopes experience intermittent snow accumulation throughout winter and spring. We tested the degree to which these contrasting water input patterns cause different near‐surface hydrologic response on north‐facing and south‐facing hillslopes during the snowmelt period. During spring snowmelt, we applied lithium bromide (LiBr) tracer to instrumented plots along a north–south catchment transect. Bromide broke through immediately at 10‐ and 30‐cm depths on the north‐facing slope and was transported out of soil waters within 40 days. On the south‐facing slope, Br^? was transported to significant depths only during spring storms and remained above the detection limit throughout the study. Modelling of unsaturated zone hydrologic response using Hydrus‐1D corroborated these aspect‐driven differences in subsurface transport. Our multiple lines of evidence suggest that north‐facing slopes are dominated by connected flow through the soil matrix, whereas south‐facing slope soils experience brief periods of rapid vertical transport following snowmelt events and are drier overall than north‐facing slopes. These differences in hydrologic response were largely a function of energy‐driven differences in water supply, emphasizing the importance of aspect and climate forcing when considering contributions of water and solutes to streamflow in catchments near the snow line. Copyright © 2012 John Wiley & Sons, Ltd.     

The role of peri‐glacial active layer development in determining soil‐regolith thickness across a Triassic sandstone outcrop in the UK

This paper examines the weathering processes that have combined to produce the distribution of soil‐regolith (SR) thickness across the Triassic Sherwood Sandstone Group outcrop (750 km²) in Nottinghamshire, UK. Archive borehole logs (n = 282) taken across the outcrop showed that SR thickness had mean and median depths of ~1·8 and 1·5 m, respectively. Cores were taken from a forested site to depths ~3 m for geochemical analysis. At this site the SR thickness was ~1·7 m. Analysis of the loss of elements, compared to bedrock using mass balance calculations (τ) showed that all the calcite and gypsum cement had been removed to depths of >3 m. Thus the major difference between the SR and the underlying saprolite was that the former exists as loose sand as opposed to a semi‐durable rock. Scanning electron microscopy (SEM) analysis of core samples suggested that the non‐durable rock or saprolite had greater cementation of clay particles. We propose that the mechanism through which the clay cement (and other interlocking grain bonds) was eased apart was through freeze–thaw processes associated with the summer ‘active layer development (ALD)’ during the last glacial activity in the UK. We tested this theory by developing a Monte Carlo simulation based on a simplified version of the Stefan equation. Current Arctic datasets of air and ground temperatures were obtained to provide reasonable starting conditions for input variables. These were combined with known data for thermal conductivity, bulk density and moisture content of the Sherwood Sandstone regolith. Model predictions (n = 1000) of the distribution of SR thickness accurately reflect the observed distribution thickness from the borehole logs. This is strong evidence that freeze–thaw and ‘ALD’ processes are major factors in determining the thickness of SR across this outcrop. British Geological Survey © NERC 2012     

Detachment and infiltration variations as consequence of regolith development in a Pyrenean badland system

The Araguás experimental catchment has been monitored to study badland dynamics in the Central Pyrenees. Previous studies of weathering processes within the catchment reported strong regolith dynamics associated with seasonal variations in the temperature and moisture regimes. A preliminary analysis of hydrological response and suspended sediment transport data recorded at a gauging station also demonstrated seasonal trends. The main objective of the present study is to understand the effect of regolith dynamics on sediment detachment and infiltration processes, based on field studies using simulated rainfall. The experiment design was based on seasonal differences in the physical conditions of surface regolith and the general trends of hydro‐sedimentological responses. Rainfall simulations were conducted on small plots using a pressure nozzle. Similar experimental rainfall conditions were set for all plots (rainfall intensity around 45 mm h^–1). The results showed strong variations in the infiltration and detachment responses closely associated with the regolith conditions and crusting development. Infiltration showed seasonal differences in time lag and intensity: average infiltration rates ranged from very low (2·05 mm h^–1) to moderated high values (44·04 mm h^–1) associated to regolith development conditions. Maximum sediment concentration, as an indicator of particles produced by detachment, also ranged from moderate (3 g l^–1) to extreme values (145 g l^–1). Mean and minimum infiltration rates showed negative correlations with initial moisture content. Sediment concentration showed a positive correlation with time lag, ponding, and sealing time, and a negative correlation with initial moisture. In terms of seasonal trends, infiltration and erosion responses were relatively stable during spring and autumn, whereas wide variations were recorded in infiltration rates and sediment detachment during summer and winter. As a general conclusion, the obtained results indicate that seasonal differences in detachment and infiltration depend on the nature of regolith development. Copyright © 2009 John Wiley & Sons, Ltd.     

Nano-texture of penetrative foliation in metamorphic rocks

By observing four samples obtained from Jiangxi Province, China, under the scanning electron mi-croscope (SEM), we discovered that nano-particle layers were commonly formed on sliding planes of the penetrative foliation in metamorphic rocks. We also successfully reproduced this phenomenon with a tri-axial pressure experiment. Having gone through the granulitization-alienation-partition in the shear sliding process, the nano-particles (40-95 nm in diameter) display different individual shapes and dis-tinct layered textures. This nano-confinement layer is essentially a frictional-viscous stripe with vis-cous-elastic deformation. In the micro-domain stripe, activities in structural stress field-rheological physical field-geochemical field are very dynamic, corresponding to the three stages (i.e., shear sliding strengthening-weakening-exfoliating) of the foliation development in metamorphism rocks. As such, the viscous-elastic deformation behavior helps shed light on the understanding of the micro-dynamic mechanism of the structural shearing.     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 2. Controls on weathering and permeability profiles,drainage hydraulics,and solute export pathways

The advance of a chemical weathering front into the bedrock of a hillslope is often limited by the rate weathering products that can be carried away, maintaining chemical disequilibrium. If the weathering front is within the saturated zone, groundwater flow downslope may affect the rate of transport and weathering—however, weathering also modifies the rock permeability and the subsurface potential gradient that drives lateral groundwater flow. This feedback may help explain why there tends to be neither “runaway weathering” to great depth nor exposed bedrock covering much of the earth and may provide a mechanism for weathering front advance to keep pace with incision of adjacent streams into bedrock. This is the second of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory. Here, we show how a simplified kinetic model of 1‐D rock weathering can be extended to consider lateral flow in a 2‐D hillslope. Exact and approximate analytical solutions for the location and thickness of weathering within the hillslope are obtained for a number of cases. A location for the weathering front can be found such that lateral flow is able to export weathering products at the rate required to keep pace with stream incision at steady state. Three pathways of solute export are identified: “diffusing up,” where solutes diffuse up and away from the weathering front into the laterally flowing aquifer; “draining down,” where solutes are advected primarily downward into the unweathered bedrock; and “draining along,” where solutes travel laterally within the weathering zone. For each pathway, a different subsurface topography and overall relief of unweathered bedrock within the hillslope is needed to remove solutes at steady state. The relief each pathway requires depends on the rate of stream incision raised to a different power, such that at a given incision rate, one pathway requires minimal relief and, therefore, likely determines the steady‐state hillslope profile.     

Quartz luminescence sensitivity from sediment versus bedrock in highly weathered soils of the Piedmont of North Carolina,south-eastern USA

Deeply weathered soils cover most of the Piedmont physiographic province of the south-eastern United States of America (USA). These soils have traditionally been inferred to derive from weathered bedrock, but recent work (e.g. Ferguson et al., 2019) suggests that deposited sediments are more prevalent than recognized. Distinguishing sediment from weathered bedrock is integral to understanding critical-zone processes and overall Quaternary landscape evolution, yet the well-developed, red, clay-dominated Ultisols of this temperate and humid region mask differences between transported from non-transported material. Our goal is to determine if optically stimulated luminescence (OSL) methods can distinguish quartz sand from allochthonous (e.g. transported sediment) versus autochthonous (e.g. in situ weathered bedrock) material in soil-profile and core samples from the Redlair Observatory in southwestern North Carolina, USA. Here, we turn to OSL sensitivity and linear-modulated OSL (LM-OSL) to observe the intensity or lack thereof of the fast-decay luminescence component (most light-sensitive signal) in quartz grains from soil horizons and crystalline bedrock-derived saprolite. We find that quartz grains sampled from in situ weathered bedrock as well as from saprolotized clasts of rock have weak luminescence properties and are not dominated by a fast-decay luminescence component. In contrast, quartz grains from transported sediment (e.g. mobile regolith; colluvium; alluvium) contain sensitive grains with more dominant fast components. These results suggest that quartz luminescence sensitivity can be a tool to differentiate between in situ weathered bedrock and similar looking mobile regolith and colluvium over-printed by soil development.     

Identification of Strong Earthquake Sites in Hebei Region by Making Use of Cluster Method

The Cluster and Hamming methods are used in this paper for a comprehensive study on geology,geomorphology,geophysical field,crustal deformation,active faults,regional stress axes and their relation in Hebei region.Fourteen potential seismic zones in which shocks with M≥6 may happen have been identified.Shocks with M≥6 have occurred in seven of them,and the others have been considered as a future strong earthquake areas.Both the K value and testing of deleting nodes show the stability of results obtained in this paper.The potential seismic zones identified in the paper fall into the areas of marked risk areas within 10 years in North China,but the scale of the identified zones is smaller.The Datong-Yanggao earthquake with M-6.1 occurred in October 1989 precisely in the 14th potential seismic zone mentioned above.     

Analytical solution for enhanced recharge around a bedrock exposure caused by deep-aquifer dewatering through a variable thickness aquitard

In this study an analytical solution was developed to predict steady radially-symmetric percolation rates from an aquifer underlain by a variable thickness aquitard. The solutions consider an aquitard with constant thickness and with radial-symmetrically increasing thickness outward from the center. The solution was used to predict the percolation rate from a peat layer around a bedrock outcrop in the James Bay Lowland near the De Beers Victor diamond mine. In this case the marine sediment layer limited the direct connection between the peat layer and the bedrock as an aquitard. Our zero order solution with constant marine sediment thickness showed the best fit to the steady state water level data of June 2012. It was found that the enhanced recharge around bioherms (i.e., at rates greater than the regional average of 0.7 mm/day) will only occur in marine sediments less than 4.3 m thick, for extreme depressurization of 30 m.     

Asymmetry of weathering‐limited hillslopes: the importance of diurnal covariation in solar insolation and temperature

Hillslope asymmetry, i.e. variation in hillslope form as a function of slope aspect and/or mean solar insolation, has been documented in many climates and geologic contexts. Such patterns have the potential to help us better understand the hydrologic, ecologic, and geomorphologic processes and feedbacks operating on hillslopes. Here we document asymmetry in the fraction of hillslope relief accommodated by cliffs in weathering‐limited hillslopes of drainage basins incised into the East Kaibab Monocline (northern Arizona) and Raplee Ridge Monocline (southern Utah) of the southern Colorado Plateau. We document that south‐ and west‐facing hillslopes have a larger proportion of hillslope relief accommodated by cliffs compared with north‐ and east‐facing hillslopes. Cliff abundance correlates positively with mean solar insolation and, by inference, negatively with soil/rock moisture. Solar insolation control of hillslope asymmetry is an incomplete explanation, however, because it cannot account for the fact that the greatest asymmetry occurs between southwest‐ and northeast‐facing hillslopes rather than between south‐ and north‐facing hillslopes in the study sites. Modeling results suggest that southwest‐facing hillslopes are more cliff‐dominated than southeast‐facing hillslopes of the same mean solar insolation in part because potential evapotranspiration rates, which control the soil/rock moisture that drives weathering, are controlled by the product of solar insolation and a nonlinear function of surface temperature, together with the fact that southwest‐facing hillslopes receive peak solar insolation during warmer times of day compared with southeast‐facing hillslopes. The dependence of water availability on both solar insolation and surface temperature highlights the importance of the diurnal cycle in controlling water availability, and it provides a general explanation for the fact that vegetation cover tends to exhibit the greatest difference between northeast‐ and southwest‐facing hillslopes in the Northern Hemisphere and between southeast‐ and northwest‐facing hillslopes in the Southern Hemisphere. Copyright © 2017 John Wiley & Sons, Ltd.     

Subsurface architecture of the Boulder Creek Critical Zone Observatory from electrical resistivity tomography

The architecture of the critical zone includes the distribution, thickness, and contacts of various types of slope deposits and weathering products such as saprolite and weathered bedrock resting on solid bedrock. A quantitative analysis of architecture is necessary for many model‐driven approaches used by pedologic, geomorphic, hydrologic or biologic studies. We have used electrical resistivity tomography, a well‐established geophysical technique causing minimum surficial disturbance, to portray the subsurface electrical resistivity differences at three study sites (Green Lakes Valley; Gordon Gulch; Betasso) at the Boulder Creek Critical Zone Observatory (BcCZO). Possible limitations of the technique are discussed. Interpretation of the specific resistivity values using natural outcrops, pits, roadcuts and drilling data as ground truth information allows us to image the critical zone architecture of each site. Green Lakes Valley (3700 MASL), a glacially eroded alpine basin, shows a rather simple, split configuration with coarse blockfields and sediments, partly containing permafrost above bedrock. The critical zone in Gordon Gulch (2650 MASL), a montane basin with rolling hills, and Betasso (1925 MASL), a lower montane basin with v‐shaped valleys, is more variable due to a complex Quaternary geomorphic history. Boundaries between overlying stratified slope deposits and saprolite were identified at mean depths of 3.0 ± 2.2 m and 4.1 ± 3.6 m in the respective sites. The boundary between saprolite and weathered bedrock is deeper in Betasso at 5.8 ± 3.7 m, compared with 4.3 ± 3.0 m in Gordon Gulch. In general, the data are consistent with results from seismic studies, but electrical resistivity tomography documents a 0.5–1.5 m shallower critical zone above the weathered bedrock on average. Additionally, we document high lateral variability, which results from the weathering and sedimentation history and seems to be a consistent aspect of critical zone architecture within the BcCZO. Copyright © 2013 John Wiley & Sons, Ltd.     

兰桥排土场边坡失稳模式及其稳定性数值分析

以兰桥矿区排土场边坡工程为例,在掌握排土场边坡工程地质条件的基础上,详细分析该多台阶排土场边坡失稳的破坏模式和主要影响因素,采用极限平衡法和有限元强度折减法对其稳定性进行分析。计算结果表明,该边坡可能沿着排土场下卧腐殖土层滑动,塑性区首先在腐殖土层贯通,接着主要向上产生2条塑性区变化带,其中沿着中部台阶边坡发展的塑性区较早贯通,而塑性区朝坡顶的贯通相对滞后;随着地下水位的升高,竖向拉应力区有扩展趋势,其产生的附加位移越大,对边坡的安全越不利;以竖向增量位移的变化趋势作为监测和分析边坡稳定性的依据时,应避开位移为零或变化很小的坡段;水平增量位移集中于坡脚,并且在边坡中部偏下的位置开始产生剧变,这种不连续的过渡,更加剧了边坡失稳的可能性。