Landscape response to tipping points in granite weathering: The case of stepped topography in the Southern Sierra Critical Zone Observatory

The dynamics of granitic landscapes are modulated by bimodal weathering, which produces patchy granular soils and expanses of bare rock ranging from meter-scale boulders to mountain-scale domes. We used terrain analysis and with cosmogenic nuclide measurements of erosion rates to quantitatively explore Wahrhaftig’s decades-old hypothesis for the development of “stepped topography” by differential weathering of bare and soil-mantled granite. According to Wahrhaftig’s hypothesis, bare granite weathers slower than soil-mantled granite; thus random erosional exposure of bare rock leads to an alternating sequence of steep, slowly weathering bedrock “steps” and gently sloped, but rapidly weathering, soil-mantled “treads.” Our investigation focused on the terrain surrounding the Southern Sierra Critical Zone Observatory (CZO), which is underlain by granitic bedrock and lies outside the limits of recent glaciation, in the heart of the stepped topography described by Wahrhaftig. Our digital terrain analysis confirms that steep steps often grade into gentle treads, consistent with Wahrhaftig’s hypothesis. However, we observe a mix-and-match of soil and bare rock on treads and steps, contrary to one of the hypothesis’ major underpinnings – that bare rock should be much more common on steps than on treads. Moreover, the data show that bare rock is not as common as expected at step tops; Wahrhaftig’s hypothesis dictates that step tops should act as slowly eroding base levels for the treads above them. The data indicate that, within each landscape class (i.e., the steps and treads), bare rock erodes more slowly than surrounding soil. This suggests that the coupling between soil production and denudation in granitic landscapes harbors a tipping point wherein erosion rates decrease when soils are stripped to bedrock. Although broadly consistent with the differential weathering invoked by Wahrhaftig, the data also show that steps are eroding faster than treads, undermining Wahrhaftig’s explanation for the origins of the steps. The revised interpretation proposed here is that the landscape evolves by back-wearing of steps in addition to differential erosion due to differences in weathering of bare and soil-mantled granite.     

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.     

Variability of stream extents controlled by flow regime and network hydraulic scaling

Stream networks expand and contract through time, impacting chemical export, aquatic habitat, and water quality. Although recent advances improve prediction of the extent of the wetted channel network (L ) based on discharge at the catchment outlet (Q ), controls on the temporal variability of L remain poorly understood and unquantified. Here we develop a quantitative, conceptual framework to explore how flow regime and stream network hydraulic scaling factors co-determine the relative temporal variability in L (denoted here as the total wetted channel drainage density). Network hydraulic scaling determines how much L changes for a change in Q , while the flow regime describes how Q changes in time. We compiled datasets of co-located dynamic stream extent mapping and discharge to analyze all globally available empirical data using the presented framework. We found that although variability in L is universally damped relative to variability in Q (i.e., streamflow is relatively more variable in time than network extent), the relationship is elastic, meaning that for a given increase in the variability in Q , headwater catchments will experience greater-than-proportional increases in the variability of L . Thus, under anticipated climatic shifts towards more volatile precipitation, relative variability in headwater stream network extents can be expected to increase even more than the relative variability of discharge itself. Comparison between network extents inferred from the L -Q relationship and blue lines on USGS topographic maps shows widespread underestimation of the wetted channel network by the blue line network.     

Picocyanobacterial abundances and diversity in surface water of the northwestern Pacific Ocean

To understand picocyanobacterial distribution patterns in the northwestern Pacific Ocean, their abundances and genetic diversity were studied using flow cytometry and a barcoded amplicon pyrosequencing approach. At open ocean stations affected by the North Equatorial Current, Prochlorococcus was the predominant picocyanobacteria, and a high-light-adapted ecotype (HLII) made up most of the population. In contrast, at stations in shelf areas of the East China Sea (ECS) and South Sea, Synechococcus was the predominant picocyanobacteria and clade II was dominant. At other ECS stations affected by the Kuroshio Current, both Prochlorococcus and Synechococcus made up similar proportions of the picocyanobacterial community. These results indicate that picocyanobacterial diversity differs among oceanic regions, and that physicochemical properties related to dominant water masses, seem to be important in determining picocyanobacterial diversity.     

Observation of bottom water renewal and export production in the japan basin, east sea using tritium and helium isotopes

Tritium (₃H or T) has been produced mostly by atmospheric nuclear weapon tests, and entered the ocean in the form of water (HTO). As tritium exists as water itself, it has been regarded as an ideal tool to study the transport of water masses. In April 2001 we collected water samples in the western Japan Basin (WJB) for tritium and helium measurement. The timely sampling provided direct evidence of the bottom water formation, resulting in the drastic increase in tritium concentration from 0.3 TU in 2000 to 0.67 TU in 2001. Considering that the new bottom waters were found mostly in the WJB, it implies that maximum 1% of the whole bottom layer below 2600 m should be replaced with the surface water during the severely cold winter 2000—2001.₃H-₃He age, showing the elapsed time since the water left from the surface, can be used to calculate oxygen utilization rate by dividing AOU by the age. Under the condition of 90% oxygen saturation in the surface water, the integration of OUR in the water column below 200 m yields net oxygen consumption of 12 mol (O₂) m^-2 yr^-1, which corresponds to the export production of 99 g C m^-2 yr^-1 . This estimate is comparable to a previous estimate based on satellite data and implies that the ratio of export to primary production (f -ratio) is as high as 0.5 in the WJB.     

Quantification of the seasonal hillslope water storage that does not drive streamflow

The relationship between seasonal catchment water storage and discharge is typically nonunique due to water storage that is not directly hydraulically connected to streams. Hydraulically disconnected water volumes are often ecologically and hydrologically important but cannot be explicitly estimated using current storage–discharge techniques. Here, we propose that discharge is explicitly sensitive to changes in only some fraction of seasonally dynamic storage that we call “direct storage,” whereas the remaining storage (“indirect storage”) varies without directly influencing discharge. We use a coupled mass balance and storage–discharge function approach to partition seasonally dynamic storage between these 2 pools in the Northern California Coast Ranges. We find that indirect storage constitutes the vast majority of dynamic catchment storage, even at the wettest times of the year. Indirect storage exhibits lower variability over the course of the wet season (and in successive winter periods) than does direct storage. Predicted indirect storage volumes and dynamics match field observations. Comparison of 2 neighbouring field sites reveals that indirect storage volumes can occur as unsaturated storage held under tension in soils and weathered bedrock and as near‐surface saturated storage that remains on hillslopes (and is eventually evapotranspired). Indirect storage volumes (including moisture in the weathered bedrock) may support plant transpiration, and our method indicates that this important water source could be quantified from precipitation and stream discharge records.     

Determining hypocentral parameters for local earthquakes under ill conditions using genetic algorithm

We demonstrate that GA-MHYPO determines accurate hypocentral parameters for local earthquakes under ill conditions, such as limited number of stations (phase data), large azimuthal gap, and noisy data. The genetic algorithm (GA) in GA-MHYPO searches for the optimal 1-D velocity structure which provides the minimum traveltime differences between observed (true) and calculated P and S arrivals within prescribed ranges. GA-MHYPO is able to determine hypocentral parameters more accurately in many circumstances than conventional methods which rely on an a priori (and possibly incorrect) 1-D velocity model. In our synthetic tests, the accuracy of hypocentral parameters obtained by GA-MHYPO given ill conditions is improved by more than a factor of 20 for error-free data, and by a factor of five for data with errors, compared to that obtained by conventional methods such as HYPOINVERSE. In the case of error-free data, GA-MHYPO yields less than 0.1 km errors in focal depths and hypocenters without strong dependence on azimuthal coverage up to 45°. Errors are less than 1 km for data with errors of a 0.1-s standard deviation. To test the performance using real data, a well-recorded earthquake in the New Madrid seismic zone and earthquakes recorded under ill conditions in the High Himalaya are relocated by GA-MHYPO. The hypocentral parameters determined by GA-MHYPO under both good and ill conditions show similar computational results, which suggest that GA-MHYPO is robust and yields more reliable hypocentral parameters than standard methods under ill conditions for natural earthquakes.