Run‐off processes from mountains to foothills: The role of soil stratigraphy and structure in influencing run‐off characteristics across high to low relief landscapes

The critical zone features that control run‐off generation, specifically at the regional watershed scale, are not well understood. Here, we addressed this knowledge gap by quantitatively and conceptually linking regional watershed‐scale run‐off regimes with critical zone structure and climate gradients across two physiographic provinces in the Southeastern United States. We characterized long‐term (~20 years) discharge and precipitation regimes for 73 watersheds with United States Geological Survey in‐stream gaging stations across the Appalachian Mountain and Piedmont physiographic provinces of North Carolina. Watersheds included in this analysis had <10% developed land and ranged in size from 14.1–4,390 km². Thirty‐four watersheds were located in the Piedmont physiographic province, which is typically classified as a low relief landscape with deep, highly weathered soils and regolith. Thirty‐nine watersheds were located in the Appalachian Mountain physiographic province, which is typically classified as a steeper landscape with highly weathered, but shallower soils and regolith. From the United States Geological Survey daily mean run‐off time series, we calculated annual and seasonal baseflow indices (BFI), minimum, mean, and maximum daily run‐off, and Pearson's correlation coefficients between precipitation and baseflow. Our results showed that Appalachian Mountain watersheds systematically had higher minimum daily flows and BFI values. Piedmont watersheds displayed much larger deviations from mean annual BFI in response to year‐to‐year variability in precipitation. A series of linear regression models between 21 landscape metrics and annual BFIs showed non‐linear and complex terrestrial–hydrological relationships across the two provinces. From these results, we discuss how distinct features of critical zone architecture, with specific focus on soil depth and stratigraphy, may be dominating the regulation of hydrological processes and run‐off regimes across these provinces.     

Application of a coil‐type TDR probe for measuring the volumetric water content in weathered granitic bedrock

As a first step toward describing water flow processes in bedrock, a coil‐type time domain reflectometry (TDR) probe capable of measuring volumetric water content, θ, in weathered bedrock at three depths was prepared. Because the coil‐type TDR probe is large in diameter (19 mm), it can be installed even in highly weathered bedrock more easily and appropriately than conventional TDR probes that consists of two or three rods of small diameter (5‐8 mm). The probe calibrations suggest that the values measured by the probe are very sensitive to changes in θ. Using the calibrated probe together with commercially available profile soil moisture sensors, the θ profile was monitored for 1 year. Even rainfall events with relatively small cumulative rainfall of 15 mm increased the bedrock θ, and the increments were comparable to those in the soil. After the end of the rainfall events, the bedrock θ displayed a more rapid drop than the soil, and varied little during the period of no rainfall. The water storage showed similar tendencies. These observations suggest that the bedrock θ is controlled by clearly distinguishable macropores and micropores within the bedrock. It is concluded that the coil‐type TDR probe is very effective in determining θ in weathered bedrock, and that bedrock, conventionally defined by conducting cone penetration tests and treated as impermeable, does conduct and hold substantial amounts of water, and therefore contribute greatly to hydrological processes in headwater catchments. Copyright © 2007 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.     

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.     

Quantifying weathering on variable rocks,an extension of geochemical mass balance: Critical zone and landscape evolution

We present a statistical model of soil and rock weathering in deep profiles to expand the capacity to assess weathering to heterogeneous bedrock types, which are common at the Earth's surface. We developed the Weathering Trends (WT) model by extending the fractional mass change calculation (tau) of the geochemical mass balance model in two important ways. First, WT log transforms the elemental ratio data, to discern the log‐linear patterns that naturally develop from thermodynamic and kinetic laws of chemistry. Second, WT statistically fits log‐transformed element concentration ratio data – log(c_j/c_i), the only depth‐varying term in tau – as a function of depth to determine characteristic depths of transitions in weathering processes, along with confidence intervals. With no prior assumptions, WT estimates average parent material composition, average composition of the upper weathered zone and mean fractional mass change of each element over the entire weathering profile. WT displays the mean shape of weathering profiles of log‐transformed geochemical data bounded by calculated confidence intervals. We share the WT model code as an open‐source R package ( https://github.com/fisherba/WeatheringTrends ). The WT model was designed to interpret two 21 m cores from the Laurels Schist bedrock in the Christina River Basin Critical Zone Observatory in the Pennsylvania Piedmont, where our morphological and elemental data provided inconclusive estimates of bedrock depth. The WT model differentiated between rock variability and weathering to delineate the maximum extent of weathering at 12.3 m (CI 95% [9.2, 21.3]) in Ridge Well 1 and 7.2 m (CI 95% [4.3, 13.0]) in Interfluve Well 2. The water table was 5–8 m below fresh rock at Ridge Well 1, but at the same depth as fresh rock at the lower elevation interfluve. We assess statistical approaches to identify the best immobile element for use in WT and tau calculations. Copyright © 2017 John Wiley & Sons, Ltd.     

Erosion of a high‐altitude,low‐relief area on the Korean Peninsula: implications for its development processes and evolution

The processes involved in the development of high‐altitude, low‐relief areas (HLAs) are still poorly understood. Although cosmogenic nuclides have provided insights into the evolution of HLAs interpreted as paleo‐surfaces, most studies focus on estimating how slowly they erode and thereby their relative stability. To understand actual development processes of HLAs, we applied several techniques of cosmogenic nuclides in the Daegwanryeong Plateau, a well‐known HLA in the Korean Peninsula. Our denudation data from strath terraces, riverine sediments, soils, and tors provide the following conclusions: (1) bedrock incision rate in the plateau (~127 m Myr^?1) is controlled by the incision rate of the western part of the Korean Peninsula, and is similar to the catchment‐wide denudation rate of the plateau (~93 m Myr^?1); (2) the soil production function we observed shows weak depth dependency that may result from highly weathered bedrock coupled with frequent frost action driven by alpine climate; (3) a discrepancy between the soil production and catchment‐wide denudation rates implies morphological disequilibrium in the plateau; (4) the tors once regarded as fossil landforms of the Tertiary do not reflect Tertiary processes; and (5) when compared with those of global paleo‐surfaces (<20 m Myr^?1), our rapid denudation rates suggest that the plateau cannot have maintained its probable initial paleo landscape, and thus is not a paleo‐surface. Our data contribute to understanding the surface processes of actively eroding upland landscapes as well as call into question conventional interpretations of supposed paleo‐surfaces around the world. Copyright © 2015 John Wiley & Sons, Ltd.     

Climate and soil moisture environment during development of the fifth palaeosol in Guanzhong Plain

Based on weathering characteristics of the fifth palaeosol layer (S5) of four sections in Guanzhong Plain, the thickness of the weathered profile of the paleosol is determined to be greater than the ordi- nary soil, a weathered and leached loess layer thicker than 2 m. The distribution depth of the red argil- lans, the weathered and leached loess layer, Fe2O3, CaCO3 and Sr content under the S5 all indicate that the precipitation in Guanzhong Plain was over 900 mm at that time. The distribution depth of gravity water zone reached 4.2 m at least, and the soil moisture content was generally more than 20% within the range of 4.2 m. At that time there was sufficient soil moisture and no dried earth layer developed in Guanzhong Plain, suitable for the forest to develop. When this soil developed, the mean annual pre- cipitation was more than the annual soil moisture evaporation. The value of soil moisture balance was positive and the atmospheric precipitation could supply the underground water normally. Soil water was weak acidic in the middle and late stages when S5 developed in Guanzhong Plain. It was a kind of subtropical climate and even more humid and warmer than the northern edge of the subtropical climate zone in Guanzhong Plain when the S5 developed. At that time the subtropical climate was prevailing over the northern side and southern side of Qingling Mountains, showing the Mountains no longer to be the boundary between the subtropical zone and the temperate zone in China. The summer monsoon acted intensely and could go over Qingling Mountains frequently bring abundant precipitation.     

Infiltration pattern in a regolith–fractured bedrock profile: field observation of a dye stain pattern

We examined the infiltration pattern of water in a regolith–bedrock profile consisting of two overburdens (OB1 and OB2), a buried rice paddy soil (PS), two texturally distinctive weathered materials (WM1 and WM2) and a fractured sedimentary rock (BR), using a Brilliant Blue FCF dye tracer. A black‐coloured coating in conducting fractures in WM1, WM2 and BR was analysed by X‐ray diffraction and scanning electron microscopy. The dye tracer penetrated to greater than 2 m depth in the profile. The macropore flow and saturated interflow were the major infiltration patterns in the profile. Macropore flow and saturated interflow were observed along fractures in WM1, WM2 and BR and at the dipping interfaces of PS–WM1, PS–WM2 and PS–BR respectively. Heterogeneous matrix flow occurred in upper overburden (OB1) and PS. Compared with OB1, the coarser textured OB2 acted as a physical barrier for vertical flow of water. The PS with low bulk density and many fine roots was another major conducting route of water in the profile. Manganese oxide and iron oxide were positively identified in the black coating material and had low crystallinity and high surface area, indicating their high reactivity with conducting contaminants. Copyright © 2006 John Wiley & Sons, Ltd.     

Modelling investigation of water partitioning at a semiarid ponderosa pine hillslope

The effects of vegetation root distribution on near‐surface water partitioning can be two‐fold. On the one hand, the roots facilitate deep percolation by root‐induced macropore flow; on the other hand, they reduce the potential for deep percolation by root‐water‐uptake processes. Whether the roots impede or facilitate deep percolation depends on various conditions, including climate, soil, and vegetation characteristics. This paper examines the effects of root distribution on deep percolation into the underlying permeable bedrock for a given soil profile and climate condition using HYDRUS modelling. The simulations were based on previously field experiments on a semiarid ponderosa pine (Pinus ponderosa) hillslope. An equivalent single continuum model for simulating root macropore flow on hillslopes is presented, with root macropore hydraulic parameterization estimated based on observed root distribution. The sensitivity analysis results indicate that the root macropore effect dominates saturated soil water flow in low conductivity soils (K_matrix below 10^?7 m/s), while it is insignificant in soils with a K_matrix larger than 10^?5 m/s, consistent with observations in this and other studies. At the ponderosa pine site, the model with simple root‐macropore parameterization reasonably well reproduces soil moisture distribution and some major runoff events. The results indicate that the clay‐rich soil layer without root‐induced macropores acts as an impeding layer for potential groundwater recharge. This impeding layer results in a bedrock percolation of less than 1% of the annual precipitation. Without this impeding layer, percolation into the underlying permeable bedrock could be as much as 20% of the annual precipitation. This suggests that at a surface with low‐permeability soil overlying permeable bedrock, the root penetration depth in the soil is critical condition for whether or not significant percolation occurs. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

An Excel®‐Based Visualization Tool of Two‐Dimensional Soil Gas Concentration Profiles in Petroleum Vapor Intrusion

In this study, we present a petroleum vapor intrusion (PVI) tool implemented in Microsoft^® Excel^® using Visual Basic for Applications and integrated within a graphical interface. The latter helps users easily visualize two‐dimensional soil gas concentration profiles and indoor concentrations as a function of site‐specific conditions such as source strength and depth, biodegradation reaction rate constant, soil characteristics and building features. This tool is based on a two‐dimensional explicit analytical model that combines steady‐state diffusion‐dominated vapor transport in a homogeneous soil with a piecewise first‐order aerobic biodegradation model, in which rate is limited by oxygen availability. As recommended in the recently released United States Environmental Protection Agency's final PVI guidance, a sensitivity analysis and a simplified Monte Carlo uncertainty analysis are also included in the spreadsheet.     

Spatial Variability of the Depth of Weathered and Engineering Bedrock using Multichannel Analysis of Surface Wave Method

In this paper an attempt has been made to evaluate the spatial variability of the depth of weathered and engineering bedrock in Bangalore, south India using Multichannel Analysis of Surface Wave (MASW) survey. One-dimensional MASW survey has been carried out at 58 locations and shear-wave velocities are measured. Using velocity profiles, the depth of weathered rock and engineering rock surface levels has been determined. Based on the literature, shear-wave velocity of 330 ± 30 m/s for weathered rock or soft rock and 760 ± 60 m/s for engineering rock or hard rock has been considered. Depths corresponding to these velocity ranges are evaluated with respect to ground contour levels and top surface levels have been mapped with an interpolation technique using natural neighborhood. The depth of weathered rock varies from 1 m to about 21 m. In 58 testing locations, only 42 locations reached the depths which have a shear-wave velocity of more than 760 ± 60 m/s. The depth of engineering rock is evaluated from these data and it varies from 1 m to about 50 m. Further, these rock depths have been compared with a subsurface profile obtained from a two-dimensional (2-D) MASW survey at 20 locations and a few selected available bore logs from the deep geotechnical boreholes.     

The spatial variation of weathering and soil depth across a Triassic sandstone outcrop

In this study, we used an archive of borehole logs from the British Geological Survey to collect information on the spatial structure of weathering that extends from the surface to competent bedrock across the Triassic Sherwood Sandstone Group outcrop (750 km²), in the East Midlands, UK. The borehole logs were used to estimate the thickness of the soil (n = 280) and soil and saprolite (S&S) to competent rock (n = 500). The weathering profile of the sandstone consisted of soil (median thickness ～ 1·5 m) overlying a transition zone of compacted and weakly cemented weathered sandstone saprolite over bedrock. Topographic analysis using a NEXTMAP 5 m × 5 m digital elevation model (DEM) revealed no significant relationships between slope properties (relief, flow length, flow accumulation or slope angle) and soil or S&S thickness. A weak, but statistically significant correlation was found between the thickness of the soil and S&S (r_s = 0·25, p < 0·001, n = 192). The variation in soil thickness may be related to changes in current and historic and land‐use, variation in sandstone properties and the influence of glacial/peri‐glacial processes. The thickness of the saprolite was more variable towards the southern part of the study area, where it increased to a maximum 40 m. We hypothesize and provide evidence that the greater weathering thickness is related to the occurrence of increased faulting in this part of the study region, allowing increased access to meteoric waters. A possible source of increased water supply is meltwater from Quaternary ice sheets; the overburden of ice may have increased sub‐glacial pore water pressure, with the fractures and faults acting as a drainage system for the removal of dissolved weathering products. British Geological Survey © NERC 2010     

Enhanced application of root‐reinforcement algorithms for bank‐stability modeling

Riparian vegetation is known to exert a number of mechanical and hydrologic controls on bank stability. In particular, plant roots provide mechanical reinforcement to a soil matrix due to the different responses of soils and roots to stress. Root reinforcement is largely a function of the strength of the roots crossing potential shear planes, and the number and diameter of such roots. However, previous bank stability models have been constrained by limited field data pertaining to the spatial and temporal variability of root networks within stream banks. In this paper, a method is developed to use root‐architecture data to derive parameters required for modeling temporal and spatial changes in root reinforcement. Changes in root numbers over time were assumed to follow a sigmoidal curve, which commonly represents the growth rates of organisms. Regressions for numbers of roots crossing potential shear planes over time showed small variations between species during the juvenile growth phase, but extrapolation led to large variations in root numbers by the time the senescent phase of the sigmoidal growth curve had been reached. In light of potential variability in the field data, the mean number of roots crossing a potential shear plane at each year of tree growth was also calculated using data from all species and an additional sigmoidal regression was run. After 30 years the mean number of roots predicted to cross a 1 m shear plane was 484, compared with species‐specific curves whose values ranged from 240 roots for black willow trees to 890 roots for western cottonwood trees. In addition, the effect of spatial variations in rooting density with depth on stream‐bank stability was modeled using the bank stability and toe erosion model (BSTEM). Three root distributions, all approximating the same average root reinforcement (5 kPa) over the top 1 m of the bank profile, were modeled, but with differing vertical distributions (concentrated near surface, non‐linear decline with depth, uniform over top meter). It was found that stream‐bank F_S varied the most when the proportion of the failure plane length to the depth of the rooting zone was greatest. Copyright © 2008 John Wiley & Sons, Ltd.     

Snowmelt‐driven macropore flow and soil saturation in a semiarid forest

Lateral subsurface flow is generally assumed to occur as a result of the development of a saturated zone above a low‐permeability interface such as at the soil–bedrock contact, and it is often augmented by macropore flow. Our objective was to evaluate the development of lateral subsurface flow and soil saturation at a semiarid ponderosa pine forest in New Mexico with respect to the conceptual model of saturation building above the soil–bedrock contact. At this site, we have long‐term observations of the water budget components, including lateral flow. A 1·5 m deep by 7 m long trench was constructed to observe lateral subsurface flow and development of saturation directly. Our observations are based on flow resulting from a melting snowdrift. The edge of the drift was about 7 m upslope from the trench. Lateral subsurface flow only occurred from root macropores in the Bt soil horizon. Saturation developed and grew outward from flowing root macropores, rather than growing upward from the soil–bedrock interface. This macropore‐centred saturation resulted in a highly heterogeneous distribution of water content until enough macropores began flowing and individual macropore saturated zones grew large enough to coalesce and saturate large volumes of the soil. Our observations are based on one snowmelt event and a relatively short hillslope flow path, and thus do not represent a full range of hydrologic conditions. Nevertheless, the observed behaviour did not conform to the traditional model of soil–bedrock control of saturation and lateral flow. Copyright © 2004 John Wiley & Sons, Ltd.     

Delineating depth to bedrock beneath shallow unconfined aquifers: a gravity transect across the Palmer River Basin

In view of the increasing demand on ground water supplies in the northeastern United States, it is imperative to develop appropriate methods to geophysically characterize the most widely used sources of ground water in the region: shallow unconfined aquifers consisting of well-sorted, stratified glacial deposits laid down in bedrock valleys and channels. The gravity method, despite its proven value in delineating buried bedrock valleys elsewhere, is seldom used by geophysical contractors in this region. To demonstrate the method's effectiveness for evaluating such aquifers, a pilot study was undertaken in the Palmer River Basin in southeastern Massachusetts. Because bedrock is so shallow beneath this aquifer (maximum depth is 30 m), the depth-integrated mass deficiency of the overlying unconsolidated material was small, so that the observed gravity anomaly was on the order of 1 milligal (mGal) or less. Thus data uncertainties were significant. Moreover, unlike previous gravity studies elsewhere, we had no a priori information on the density of the sediment. Under such circumstances, it is essential to include model constraints and weighted least-squares in the inversion procedure. Among the model constraints were water table configuration, bedrock outcrops, and depth to bedrock from five water wells. Our procedure allowed us to delineate depth to bedrock along a 3.5 km profile with a confidence interval of 1.8 m at a nominal depth of 17 m. Moreover, we obtained a porosity estimate in the range of 39% to 44%. Thus the gravity method, with appropriate refinements, is an effective tool for the reconnaissance of shallow unconfined aquifers.     

Role of bedrock groundwater in the rainfall–runoff process in a small headwater catchment underlain by volcanic rock

The role of bedrock groundwater in rainfall–runoff processes is poorly understood. Hydrometric, tracer and subsurface water potential observations were conducted to study the role of bedrock groundwater and subsurface flow in the rainfall–runoff process in a small headwater catchment in Shiranui, Kumamoto prefecture, south‐west Japan. The catchment bedrock consists of a strongly weathered, fractured andesite layer and a relatively fresh continuous layer. Major chemical constituents and stable isotopic ratios of δ¹⁸O and δD were analysed for spring water, rainwater, soil water and bedrock groundwater. Temporal and spatial variation in SiO₂ showed that stream flow under the base flow condition was maintained by bedrock groundwater. Time series of three components of the rainstorm hydrograph (rainwater, soil water and bedrock groundwater) separated by end member mixing analysis showed that each component fluctuated during rainstorm, and their patterns and magnitudes differed between events. During a typical mid‐magnitude storm event, a delayed secondary runoff peak with 1·0 l s⁻¹ was caused by increase in the bedrock groundwater component, whereas during a large rainstorm event the bedrock groundwater component increased to ≈ 2·5 l s⁻¹. This research shows that the contribution of bedrock groundwater and soil water depends strongly on the location of the groundwater table, i.e. whether or not it rises above the soil–bedrock interface. Copyright © 2010 John Wiley & Sons, Ltd.     

Spatial and temporal variability of 0‐ to 5‐m soil–water storage at the watershed scale

Dynamic relationships among rainfall patterns, soil water distribution, and plant growth are crucial for sustainable conservation of soil and water resources in water‐limited ecosystems. Spatial and temporal variation in deep soil water content at a watershed scale have not yet been characterized adequately due to the lack of deep soil water data. Deep soil–water storage (SWS) up to a depth of 5 m (n = 73) was measured at 19 sampling occasions at the LaoYeManQu watershed on the Chinese Loess Plateau (CLP). At a depth of 0–1.5 m, the annual mean SWS was highly correlated with rain intensity, and the correlation decreased with depth, but within the layers at 1.5–5.0 m, the changes in SWS indicated a lag between precipitation and the replenishment of soil water. Geostatistical parameters of SWS were also highly dependent on depth, and the mean SWS presented similar spatial structures in two adjacent layers. Temporal stability of SWS as indicated by mean relative difference, standard deviation of the relative difference (SDRD), and mean absolute bias error (MABE) was significantly weaker at the shallow than at deeper layers. Soil separates and organic carbon content controlled the spatial pattern of SWS at the watershed scale. One representative location (Site 57) was identified to estimate the mean SWS in the 1‐ to 5‐m layer of the watershed. Semivariograms of the SDRD and MABE were best fitted by an isotropic spherical model, and their spatial distributions were depth‐dependent. Both temporal stability and spatial variability of SWS increased over depth. This study is helpful for deep SWS estimation and sustainable management of soil and water on the CLP, and for other similar regions around the world.     

Geophysical constraints on deep weathering and water storage potential in the Southern Sierra Critical Zone Observatory

The conversion of bedrock to regolith marks the inception of critical zone processes, but the factors that regulate it remain poorly understood. Although the thickness and degree of weathering of regolith are widely thought to be important regulators of the development of regolith and its water‐storage potential, the functional relationships between regolith properties and the processes that generate it remain poorly documented. This is due in part to the fact that regolith is difficult to characterize by direct observations over the broad scales needed for process‐based understanding of the critical zone. Here we use seismic refraction and resistivity imaging techniques to estimate variations in regolith thickness and porosity across a forested slope and swampy meadow in the Southern Sierra Critical Zone Observatory (SSCZO). Inferred seismic velocities and electrical resistivities image a weathering zone ranging in thickness from 10 to 35 m (average = 23 m) along one intensively studied transect. The inferred weathering zone consists of roughly equal thicknesses of saprolite (P‐velocity < 2 km s^?1) and moderately weathered bedrock (P‐velocity = 2–4 km s^?1). A minimum‐porosity model assuming dry pore space shows porosities as high as 50% near the surface, decreasing to near zero at the base of weathered rock. Physical properties of saprolite samples from hand augering and push cores are consistent with our rock physics model when variations in pore saturation are taken into account. Our results indicate that saprolite is a crucial reservoir of water, potentially storing an average of 3 m³ m^?2 of water along a forested slope in the headwaters of the SSCZO. When coupled with published erosion rates from cosmogenic nuclides, our geophysical estimates of weathering zone thickness imply regolith residence times on the order of 10⁵ years. Thus, soils at the surface today may integrate weathering over glacial–interglacial fluctuations in climate. Copyright © 2013 John Wiley & Sons, Ltd.     

Joint Audio-Magnetotelluric and Passive Seismic Imaging of the Cerdanya Basin

The structure of Cerdanya Basin (north-east of Iberian Peninsula) is partly known from geological cross sections, geological maps and vintage geophysical data. However, these data do not have the necessary resolution to characterize some parts of Cerdanya Basin such as the thickness of soft soil, geometry of bedrock or geometry of geological units and associated faults. For all these reasons, the main objective of this work is to improve this deficiency carrying out a detailed study in this Neogene basin applying jointly the combination of passive seismic methods (H/V spectral ratio and seismic array) and electromagnetic methods (audio-magnetotelluric and magnetotelluric method). The passive seismic techniques provide valuable information of geometry of basement along the profile. The maximum depth is located near Alp village with a bedrock depth of 500 m. The bedrock is located in surface at both sites of profile. The Neogene sediments present a shear-wave velocity between 400 and 1000 m/s, and the bedrock basement presents a shear-wave velocity values between 1700 and 2200 m/s. These results are used as a priori information to create a 2D resistivity initial model which constraints the inversion process of electromagnetic data. We have obtained a 2D resistivity model which is characterized by (1) a heterogeneous conductivity zone (<40 Ohm m) that corresponds to shallow part of the model up to 500 m depth in the centre of the profile. These values have been associated with Quaternary and Neogene sediments formed by silts, clays, conglomerates, sandstones and gravels, and (2) a deeper resistive zone (1000–3000 Ohm m) interpreted as Palaeozoic basement (sandstones, limestones and slates at NW and conglomerates and microconglomerates at SE). The resistive zone is truncated by a discontinuity at the south-east of the profile which is interpreted as the Alp-La Tet Fault. This discontinuity is represented by a more conductive zone (600 Ohm m approx.) and is explained as a combination of fractured rock and a fluid network. The result highlights that the support between different geophysical methods is essential in producing geophysical meaningful models.