Architecture of the deep critical zone in the Río Icacos watershed (Luquillo Critical Zone Observatory,Puerto Rico) inferred from drilling and ground penetrating radar (GPR)

In the critical zone, surficial bedrock interactions result in the formation of a mantle of chemically‐ and physically‐altered material defined here as regolith. In the watershed of the Río Icacos, an upland river draining the Luquillo Mountains in tropical Puerto Rico, we explored the influence of lithology (quartz diorite versus hornfels‐facies volcaniclastic rock) on weathering. Regolith profiles were studied by drilling boreholes and imaging the subsurface using ground penetrating radar (GPR). Overall, the regolith structure is not laterally continuous but rather is punctuated by zones of deep fractures that host in situ weathering, corestones, and colluvial material. GPR images of these vertical zones show reflectors at 15–20 m depth. Thus, the architecture of the critical zone in the upper Luquillo Mountains is highly dependent on lithology and its influence on fracture development. At the highest elevations where hornfels overlies quartz diorite, positive feedbacks occur when the water table drops so that oxidative weathering of biotite in the more felsic rock creates microfractures and allows deeper infiltration of meteoric waters. Such exposure results in some of the fastest weathering rocks in the world and may contribute to formation of the knickpoint in the Río Icacos watershed. This work represents the first study combining GPR and drilling to look at the structure of the deep critical zone and demonstrates: (1) the importance of combining direct methods (such as drilling) with indirect methods (such as GPR) to understand the architecture of the critical zone in tropical systems; (2) the interplay of the surficial stress regime, lithology and climate in dictating the architecture of weathering. Copyright © 2016 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.     

Mass‐balance modeling of mineral weathering rates and CO2 consumption in the forested,metabasaltic Hauver Branch watershed,Catoctin Mountain,Maryland, USA

Mineral weathering rates and a forest macronutrient uptake stoichiometry were determined for the forested, metabasaltic Hauver Branch watershed in north‐central Maryland, USA. Previous studies of Hauver Branch have had an insufficient number of analytes to permit determination of rates of all the minerals involved in chemical weathering, including biomass. More equations in the mass‐balance matrix were added using existing mineralogic information. The stoichiometry of a deciduous biomass term was determined using multi‐year weekly to biweekly stream‐water chemistry for a nearby watershed, which drains relatively unreactive quartzite bedrock. At Hauver Branch, calcite hosts ~38 mol% of the calcium ion (Ca²⁺) contained in weathering minerals, but its weathering provides ~90% of the stream water Ca²⁺. This occurs in a landscape with a regolith residence time of more than several Ka (kiloannum). Previous studies indicate that such old regolith does not typically contain dissolving calcite that affects stream Ca²⁺/Na⁺ ratios. The relatively high calcite dissolution rate likely reflects dissolution of calcite in fractures of the deep critical zone. Of the carbon dioxide (CO₂) consumed by mineral weathering, calcite is responsible for approximately 27%, with the silicate weathering consumption rate far exceeding that of the global average. The chemical weathering of mafic terrains in decaying orogens thus may be capable of influencing global geochemical cycles, and therefore, climate, on geological timescales. Based on carbon‐balance calculations, atmospheric‐derived sulfuric acid is responsible for approximately 22% of the mineral weathering occurring in the watershed. Our results suggest that rising air temperatures, driven by global warming and resulting in higher precipitation, will cause the rate of chemical weathering in the Hauver Branch watershed to increase until a threshold temperature is reached. Beyond the threshold temperature, increased recharge would produce a shallower groundwater table and reduced chemical weathering rates. Copyright © 2013 John Wiley & Sons, Ltd.     

Controls on deep critical zone architecture: a historical review and four testable hypotheses

The base of Earth's critical zone (CZ) is commonly shielded from study by many meters of overlying rock and regolith. Though deep CZ processes may seem far removed from the surface, they are vital in shaping it, preparing rock for infusion into the biosphere and breaking Earth materials down for transport across landscapes. This special issue highlights outstanding challenges and recent advances of deep CZ research in a series of articles that we introduce here in the context of relevant literature dating back to the 1500s. Building on several contributions to the special issue, we highlight four exciting new hypotheses about factors that drive deep CZ weathering and thus influence the evolution of life‐sustaining CZ architecture. These hypotheses have emerged from recently developed process‐based models of subsurface phenomena including: fracturing related to subsurface stress fields; weathering related to drainage of bedrock under hydraulic head gradients; rock damage from frost cracking due to subsurface temperature gradients; and mineral reactions with reactive fluids in subsurface chemical potential gradients. The models predict distinct patterns of subsurface weathering and CZ thickness that can be compared with observations from drilling, sampling and geophysical imaging. We synthesize the four hypotheses into an overarching conceptual model of fracturing and weathering that occurs as Earth materials are exhumed to the surface across subsurface gradients in stress, hydraulic head, temperature, and chemical potential. We conclude with a call for a coordinated measurement campaign designed to comprehensively test the four hypotheses across a range of climatic, tectonic and geologic conditions. Copyright © 2016 John Wiley & Sons, Ltd.     

Strong slope-aspect control of regolith thickness by bedrock foliation

The porous near-surface layer of the Earth's crust – the critical zone – constitutes a vital reservoir of water for ecosystems, provides baseflow to streams, guides recharge to deep aquifers, filters contaminants from groundwater, and regulates the long-term evolution of landscapes. Recent work suggests that the controls on regolith thickness include climate, tectonics, lithology, and vegetation. However, the relative paucity of observations of regolith structure and properties at landscape scales means that theoretical models of critical zone structure are incompletely tested. Here we present seismic refraction and electrical resistivity surveys that thoroughly characterize subsurface structure in a small catchment in the Santa Catalina Mountains, Arizona, USA, where slope-aspect effects on regolith structure are expected based on differences in vegetation. Our results show a stark contrast in physical properties and inferred regolith thickness on opposing slopes, but in the opposite sense of that expected from environmental models and observed vegetation patterns. Although vegetation (as expressed by normalized difference vegetation index [NDVI]) is denser on the north-facing slope, regolith on the south-facing slope is four times thicker (as indicated by lower seismic velocities and resistivities). This contrast cannot be explained by variations in topographic stress or conventional hillslope morphology models. Instead, regolith thickness appears to be controlled by metamorphic foliation: regolith is thicker where foliation dips into the topography, and thinner where foliation is nearly parallel to the surface. We hypothesize that, in this catchment, hydraulic conductivity and infiltration capacity control weathering: infiltration is hindered and regolith is thin where foliation is parallel to the surface topography, whereas water infiltrates deeper and regolith is thicker where foliation intersects topography at a substantial angle. These results suggest that bedrock foliation, and perhaps by extension sedimentary layering, can control regolith thickness and must be accounted for in models of critical zone development. © 2020 John Wiley & Sons, Ltd.     

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.     

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.     

In situ stress measurements in NIED boreholes in and around the fault zone near the 1995 Hyogo-ken Nanbu earthquake, Japan

Abstract The 1995 Hyogo-ken Nanbu (Kobe) earthquake, M 7.2, occurred along the north-east–south-west trending Rokko–Awaji Fault system. Three boreholes of 1001 m, 1313 m and 1838 m deep were drilled in the vicinity of the epicenter of the earthquake. Each borehole is located at characteristic sites in relation to active faults and the aftershock distribution. In particular, the Nojima–Hirabayashi borehole [Hirabayashi National Research Institute for Earth Science and Disaster Prevention (NIED) drilling] in Awaji Island was drilled to a depth of 1838 m, approximately 320 m southeast from the surface rupture of the Nojima Fault, and it crosses fracture zones below a depth of 1140 m. In situ stress measurements by the hydraulic fracturing method were conducted in these boreholes within 1.5 years after the earthquake. Measurement results suggest the following: (i) Differential stress values are very small, approximately 10 MPa at a depth of 1000 m at each site; (ii) the orientation of maximum horizontal compression is almost the same in the boreholes, perpendicular to the surface trace of the faults, north-west–south-east; (iii) fault types estimated from the state of stress differ among these sites; and (iv) the differential stress value just beneath the fault fracture zone decreases abruptly to one-half of that above the fault zone in the Hirabayashi NIED drilling. These features support the idea that the shear stress along the Rokko–Awaji Fault system decreased to a low level just after the earthquake.     

Heterogeneity in ground water–surface water interactions in the hyporheic zone of a salmonid spawning stream

Spatial and temporal variability in ground water–surface water interactions in the hyporheic zone of a salmonid spawning stream was investigated. Four locations in a 150‐m reach of the stream were studied using hydrometric and hydrochemical tracing techniques. A high degree of hydrological connectivity between the riparian hillslope and the stream channel was indicated at two locations, where hydrochemical changes and hydraulic gradients indicated that the hyporheic zone was dominated by upwelling ground water. The chemistry of ground water reflected relatively long residence times and reducing conditions with high levels of alkalinity and conductivity, low dissolved oxygen (DO) and nitrate. At the other locations, connectivity was less evident and, at most times, the hyporheic zone was dominated by downwelling stream water characterized by high DO, low alkalinity and conductivity. Substantial variability in hyporheic chemistry was evident at fine (<10 m) spatial scales and changed rapidly over the course of hydrological events. The nature of the hydrochemical response varied among locations depending on the strength of local ground water influence. It is suggested that greater emphasis on spatial and temporal heterogeneity in ground water–surface water interactions in the hyporheic zone is necessary for a consideration of hydrochemical effects on many aspects of stream ecology. For example, the survival of salmonid eggs in hyporheic gravels varied considerably among the locations studied and was shown to be associated with variation in interstitial chemistry. River restoration schemes and watershed management strategies based only on the surface expression of catchment characteristics risk excluding consideration of potentially critical subsurface processes. Copyright © 2002 John Wiley & Sons, Ltd.     

Quantifying aggregation and change in runoff source in accordance with catchment area increase in a forested headwater catchment

There has been a great deal of research interest regarding changes in flow path/runoff source with increases in catchment area. However, there have been very few quantitative studies taking subscale variability and convergence of flow path/runoff source into account, especially in relation to headwater catchments. This study was performed to elucidate how the contributions and discharge rates of subsurface water (water in the soil layer) and groundwater (water in fractured bedrock) aggregate and change with catchment area increase, and to elucidate whether the spatial variability of the discharge rate of groundwater determines the spatial variability of stream discharge or groundwater contribution. The study area was a 5‐km² forested headwater catchment in Japan. We measured stream discharge at 113 points and water chemistry at 159 points under base flow conditions. End‐member mixing analysis was used to separate stream water into subsurface water and groundwater. The contributions of both subsurface water and groundwater had large variability below 1 km². The contribution of subsurface water decreased markedly, while that of groundwater increased markedly, with increases in catchment area. The specific discharge of subsurface water showed a large degree of variability and decreased with catchment area below 0.1 km², becoming almost constant above 0.1 km². The specific discharge of groundwater showed large variability below 1 km² and increased with catchment area. These results indicated that the variabilities of stream discharge and groundwater contribution corresponded well with the variability of the discharge rate of groundwater. However, below 0.1 km², it was necessary to consider variations in the discharge rates of both subsurface water and groundwater. Copyright © 2016 John Wiley & Sons, Ltd.     

Lagged rejuvenation of groundwater indicates internal flow structures and hydrological connectivity

Large proportions of rainwater and snowmelt infiltrate into the subsurface before contributing to stream flow and stream water quality. Subsurface flow dynamics steer the transport and transformation of contaminants, carbon, weathering products and other biogeochemistry. The distribution of groundwater ages with depth is a key feature of these flow dynamics. Predicting these ages are a strong test of hypotheses about subsurface structures and time-varying processes. Chlorofluorocarbon (CFC)-based groundwater ages revealed an unexpected groundwater age stratification in a 0.47 km² forested catchment called Svartberget in northern Sweden. An overall groundwater age stratification, representative for the Svartberget site, was derived by measuring CFCs from nine different wells with depths of 2–18 m close to the stream network. Immediately below the water table, CFC-based groundwater ages of already 30 years that increased with depth were found. Using complementary groundwater flow models, we could reproduce the observed groundwater age stratification and show that the 30 year lag in rejuvenation comes from return flow of groundwater at a subsurface discharge zone that evolves along the interface between two soil types. By comparing the observed groundwater age stratification with a simple analytical approximation, we show that the observed lag in rejuvenation can be a powerful indicator of the extent and structure of the subsurface discharge zone, while the vertical gradient of the age-depth-relationship can still be used as a proxy of the overall aquifer recharge even when sampled in the discharge zone. The single age stratification profile measured in the discharge zone, close to the aquifer outlet, can reveal the main structure of the groundwater flow pattern from recharge to discharge. This groundwater flow pattern provides information on the participation of groundwater in the hydrological cycle and indicates the lower boundary of hydrological connectivity.     

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.     

Crack-filling clays and weathered cracks in the DPRI 1800 m core near the Nojima Fault, Japan: Evidence for deep surface-water circulation near an active fault

Abstract Crack-filling clays and weathered cracks were observed in the Disaster Prevention Research Institute, Kyoto University (DPRI) 1800 m cores drilled from the Nojima Fault Zone, which was activated during the 1995 Hyogo-ken Nanbu earthquake (Kobe earthquake). The crack-filling clays consist mainly of unconsolidated fine-grained materials that fill opening cracks with no shear textures. Most of the cracks observed in the DPRI 1800 m cores are yellow-brown to brown in color due to weathering. Powder X-ray diffraction analyses show that the crack-filling clays are composed mainly of clay minerals and carbonates such as siderite and calcite. Given that the top of the borehole is approximately 45 m above sea level, most of the core is far below the stable groundwater table. Hence, it is suggested that the crack-filling clays and weathered cracks in the cores taken at depths of 1800 m were formed by the flow of surface water down to the deep fractured zone of the Nojima Fault Zone during seismic faulting.     

The role of wetland coverage within the near‐stream zone in predicting of seasonal stream export chemistry from forested headwater catchments

Stream chemistry is often used to infer catchment‐scale biogeochemical processes. However, biogeochemical cycling in the near‐stream zone or hydrologically connected areas may exert a stronger influence on stream chemistry compared with cycling processes occurring in more distal parts of the catchment, particularly in dry seasons and in dry years. In this study, we tested the hypotheses that near‐stream wetland proportion is a better predictor of seasonal (winter, spring, summer, and fall) stream chemistry compared with whole‐catchment averages and that these relationships are stronger in dryer periods with lower hydrologic connectivity. We evaluated relationships between catchment wetland proportion and 16‐year average seasonal flow‐weighted concentrations of both biogeochemically active nutrients, dissolved organic carbon (DOC), nitrate (NO₃‐N), total phosphorus (TP), as well as weathering products, calcium (Ca), magnesium (Mg), at ten headwater (<200 ha) forested catchments in south‐central Ontario, Canada. Wetland proportion across the entire catchment was the best predictor of DOC and TP in all seasons and years, whereas predictions of NO₃‐N concentrations improved when only the proportion of wetland within the near‐stream zone was considered. This was particularly the case during dry years and dry seasons such as summer. In contrast, Ca and Mg showed no relationship with catchment wetland proportion at any scale or in any season. In forested headwater catchments, variable hydrologic connectivity of source areas to streams alters the role of the near‐stream zone environment, particularly during dry periods. The results also suggest that extent of riparian zone control may vary under changing patterns of hydrological connectivity. Predictions of biogeochemically active nutrients, particularly NO₃‐N, can be improved by including near‐stream zone catchment morphology in landscape models.     

Bulk density optimization to determine subsurface hydraulic properties in Rocky Mountain catchments using the GEOtop model

Integrated watershed models can be used to calculate streamflow generation in snow‐dominated mountainous catchments. Parameterization of water flow is often complicated by the lack of information on subsurface hydraulic properties. In this study, bulk density optimization was used to determine hydraulic parameters for the upper and lower regolith in the GEOtop model. The methodology was tested in two small catchments in the Dry Creek Watershed in Idaho and the Libby Creek Watershed in Wyoming. Modelling efficiencies for profile‐average soil–water content for the two catchments were between 0.52 and 0.64. Modelling efficiencies for stream discharge (cumulative stream discharge) were 0.45 (0.91) and 0.54 (0.94) for the Idaho and Wyoming catchments, respectively. The calculated hydraulic properties suggest that lateral flow across the upper–lower regolith interface is an important driver of streamflow in both the Idaho and Wyoming watersheds. The overall calibration procedure is computationally efficient because only two bulk density values are optimized. The two‐parameter calibration procedure was complicated by uncertainty in hydraulic conductivity anisotropy. Different upper regolith hydraulic conductivity anisotropy factors had to be tested in order to describe streamflow in both catchments.     

Weathering and erosion of fractured bedrock systems

We explore the contribution of fractures (joints) in controlling the rate of weathering advance for a low‐porosity rock by using methods of homogenization to create averaged weathering equations. The rate of advance of the weathering front can be expressed as the same rate observed in non‐fractured media (or in an individual block) divided by the volume fraction of non‐fractured blocks in the fractured parent material. In the model, the parent has fractures that are filled with a more porous material that contains only inert or completely weathered material. The low‐porosity rock weathers by reaction‐transport processes. As observed in field systems, the model shows that the weathering advance rate is greater for the fractured as compared to the analogous non‐fractured system because the volume fraction of blocks is < 1. The increase in advance rate is attributed both to the increase in weathered material that accompanies higher fracture density, and to the increase in exposure of surface of low‐porosity rock to reaction‐transport. For constant fracture aperture, the weathering advance rate increases when the fracture spacing decreases. Equations describing weathering advance rate are summarized in the ‘List of selected equations’. If erosion is imposed at a constant rate, the weathering systems with fracture‐bounded bedrock blocks attain a steady state. In the erosional transport‐limited regime, bedrock blocks no longer emerge at the air‐regolith boundary because they weather away. In the weathering‐limited (or kinetic) regime, blocks of various size become exhumed at the surface and the average size of these exposed blocks increases with the erosion rate. For convex hillslopes, the block size exposed at the surface increases downslope. This model can explain observations of exhumed rocks weathering in the Luquillo mountains of Puerto Rico. Published 2017. This article is a U.S. Government work and is in the public domain in the USA     

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.     

Characterization of seasonal local recharge using electrical resistivity tomography and magnetic resonance sounding

A groundwater recharge process of heterogeneous hard rock aquifer in the Moole Hole experimental watershed, south India, is being studied to understand the groundwater flow behaviour. Significant seasonal variations in groundwater level are observed in boreholes located at the outlet area indicating that the recharge process is probably taking place below intermittent streams. In order to localize groundwater recharge zones and to optimize implementation of boreholes, a geophysical survey was carried out during and after the 2004 monsoon across the outlet zone. Magnetic resonance soundings (MRS) have been performed to characterize the aquifer and measure groundwater level depletion. The results of MRS are consistent with the observation in boreholes, but it suffers from degraded lateral resolution. A better resolution of the regolith/bedrock interface is achieved using electrical resistivity tomography (ERT). ERT results are confirmed by resistivity logging in the boreholes. ERT surveys have been carried out twice—before and during the monsoon—across the stream area. The major feature of recharge is revealed below the stream with a decrease by 80% of the calculated resistivity. The time‐lapse ERT also shows unexpected variations at a depth of 20 m below the slopes that could have been interpreted as a consequence of a deep seasonal water flow. However, in this area time‐lapse ERT does not match with borehole data. Numerical modelling shows that in the presence of a shallow water infiltration, an inversion artefact may take place thus limiting the reliability of time‐lapse ERT. A combination of ERT with MRS provides valuable information on structure and aquifer properties respectively, giving a clue for a conceptual model of the recharge process: infiltration takes place in the conductive fractured‐fissured part of the bedrock underlying the stream and clayey material present on both sides slows down its lateral dissipation. Copyright © 2007 John Wiley & Sons, Ltd.     

Using hysteretic behaviour and hydrograph classification to identify hydrological function across the “hillslope–depression–stream” continuum in a karst catchment

In cockpit karst landscapes, fluxes from upland areas contribute large volumes of water to low-lying depressions and stream flow. Hydrograph hysteresis and similarity between monitoring sites is important for understanding the space–time variability of hydrologic responses across the “hillslope–depression–stream” continuum. In this study, the hysteretic feature of hydrographs was assessed by characterizing the loop-like relationships between responses at upstream sites relative to subsurface discharge at the outlet of a small karst catchment. A classification of hydrograph responses based on the multi-scale smoothing Kernel -derived distance classifies the hydrograph responses on the basis of similarities between hillslope and depression sites, and those at the catchment outlet. Results demonstrate that the temporal and spatial variability of hydrograph hysteresis and similarity between hillslope flow and outlet stream flow can be explained by the local heterogeneity of depression aquifer. Large depression storage deficits emerging in the highly heterogeneous aquifer produce strong hysteresis and multiple relationships of upstream hydrographs relative to the outlet subsurface discharge. In contrast, when depression storage deficits are filled during consecutive rainfall events, depression hydrographs at the high permeability sites are almost synchronous or exhibit a monotonous function with the hydrographs at the outlet. This reduced hydrograph hysteresis enhances preferential flow paths in fractured rocks and conduits that can accelerate the hillslope flow to the outlet. Therefore, classification of hydrograph similarities between any upstream sites and the catchment outlet can help to identify the dominant hydrological functions in the heterogeneous karst catchment.