Climate driven coevolution of weathering profiles and hillslope topography generates dramatic differences in critical zone architecture

Considerable debate revolves around the relative importance of rock type, tectonics, and climate in creating the architecture of the critical zone. We demonstrate the importance of climate and in particular the rate of water recharge to the subsurface, using numerical models that incorporate hydrologic flowpaths, chemical weathering, and geomorphic rules for soil production and transport. We track alterations in both solid phase (plagioclase to clay) and water chemistry along hydrologic flowpaths that include lateral flow beneath the water table. To isolate the role of recharge, we simulate dry and wet cases and prescribe identical landscape evolution rules. The weathering patterns that develop differ dramatically beneath the resulting parabolic interfluves. In the dry case, incomplete weathering is shallow and surface parallel, whereas in the wet case, intense weathering occurs to depths approximating the base of the bounding channels, well below the water table. Exploration of intermediate cases reveals that the weathering state of the subsurface is strongly governed by the ratio of the rate of advance of the weathering front itself controlled by the water input rate, and the rate of erosion of the landscape. The system transitions between these end‐member behaviours rather abruptly at a weathering front speed ‐ erosion rate ratio of approximately 1. Although there are undoubtedly direct roles for tectonics and rock type in critical zone architecture, and yet more likely feedbacks between these and climate, we show here that differences in hillslope‐scale weathering patterns can be strongly controlled by climate.     

Economics,Helium, and the U.S. Federal Helium Reserve: Summary and Outlook

In 2017, disruptions in the global supply of helium reminded consumers, distributors, and policy makers that the global helium supply chain lacks flexibility, and that attempts to increase production from the U.S. Federal Helium Reserve (the FHR) may not be able to compensate for the loss of one of the few major producers in the world. Issues with U.S. and global markets for helium include inelastic demand, economic availability of helium only as a byproduct, only 4–5 major producers, helium’s propensity to escape earth’s crust, an ongoing absence of storage facilities comparable to the FHR, and a lack of consequences for the venting of helium. The complex combination of these economic, physical, and regulatory issues is unique to helium, and determining helium’s practical availability goes far beyond estimating the technically accessible volume of underground resources. Although most of these issues have been analyzed since helium was recognized to be a valuable mineral commodity in the early 1900s, very few economic models have been developed that adequately consider the unique characteristics of helium and helium markets. In particular, there is a notable lack of recent empirical work to estimate the responsiveness of helium demand, supply, prices, and trade patterns to the ongoing drawdown and sale of helium reserves stored in the FHR. In general, existing models of helium either do not account for an oligopoly controlling supply, or they do not evaluate potential helium extraction and storage programs based on an intertemporal maximization of the value of the resource. Such models could be of very limited use to decision makers. This review found only one working paper with a helium market model that has incorporated both of these vital considerations. That and other economic studies along similar lines could be very useful in helping inform current helium policy discussions and decisions.     

Conditionally Averaged Large-Scale Motions in the Neutral Atmospheric Boundary Layer: Insights for Aeolian Processes

Aeolian erosion of flat, arid landscapes is induced (and sustained) by the aerodynamic surface stress imposed by flow in the atmospheric surface layer. Conceptual models typically indicate that sediment mass flux, Q (via saltation or drift), scales with imposed aerodynamic stress raised to some exponent, n, where \(n > 1\). This scaling demonstrates the importance of turbulent fluctuations in driving aeolian processes. In order to illustrate the importance of surface-stress intermittency in aeolian processes, and to elucidate the role of turbulence, conditional averaging predicated on aerodynamic surface stress has been used within large-eddy simulation of atmospheric boundary-layer flow over an arid, flat landscape. The conditional-sampling thresholds are defined based on probability distribution functions of surface stress. The simulations have been performed for a computational domain with \(\approx 25 H\) streamwise extent, where H is the prescribed depth of the neutrally-stratified boundary layer. Thus, the full hierarchy of spatial scales are captured, from surface-layer turbulence to large- and very-large-scale outer-layer coherent motions. Spectrograms are used to support this argument, and also to illustrate how turbulent energy is distributed across wavelengths with elevation. Conditional averaging provides an ensemble-mean visualization of flow structures responsible for erosion ‘events’. Results indicate that surface-stress peaks are associated with the passage of inclined, high-momentum regions flanked by adjacent low-momentum regions. Fluid in the interfacial shear layers between these adjacent quasi-uniform momentum regions exhibits high streamwise and vertical vorticity.     

Millisecond pulsar velocities

We present improved timing parameters for 13 millisecond pulsars (MSPs), including nine new proper motion measurements. These new proper motions bring to 23 the number of MSPs with measured transverse velocities. In light of these new results we present and compare the kinematic properties of MSPs with those of ordinary pulsars. The mean transverse velocity of MSPs was found to be 85±13 km s⁻¹, a value consistent with most models for the origin and evolution of MSPs and approximately a factor of 4 lower than that of ordinary pulsars. We also find that, in contrast to young ordinary pulsars, the vast majority of which are moving away from the Galactic plane, almost half of the MSPs are moving towards the plane. This near-isotropy would be expected of a population that has reached dynamic equilibrium. Accurate measurements of MSP velocities have allowed us to correct their measured spin-down rates for Doppler acceleration effects, and thereby derive their intrinsic magnetic field strengths and characteristic ages. We find that close to half of our sample of MSPs have a characteristic age comparable to or greater than the age of the Galactic disc.     

Gonadal development in Japanese medaka (Oryzias latipes) exposed to 17 β-estradiol

The purpose of this research was to evaluate the Japanese medaka (Oryzias latipes) as an ecotoxicological model for the rapid evaluation of environmental estrogens. A novel short-term (48-h) exposure to 17 β-estradiol is proposed in development of a positive control for disruption of gonadal development. Recently hatched medaka fry (30 fry per dose group) with undifferentiated gonads were exposed to 4.0, 29.4, and 115.6 μg/liter of 17 β-estradiol (acetone carrier) for 48 h in a water bath at 25 °C. The fry were then grown-out in spring water for 2 weeks, killed, and processed for histological evaluation. High lethality was encountered during the grow-out period in the 115.6 μg/liter dose group. Fry in the spring water and acetone (carrier) control groups developed into females or males. Fry exposed to 17 β-estradiol developed primarily into females or had testis-ova.     

Influence of m-cresol purple indicator additions on the pH of seawater samples: correction factors evaluated from a chemical speciation model

The perturbation of the indicator m-cresol purple on the pH in seawater is illustrated in diagrams, representing measurements in 1-cm and 5-cm cells. The diagrams apply to a measured pH interval of 7.4–8.4 using a 2-mM stock solution of m-cresol purple sodium salt dissolved in seawater. The magnitude of the perturbation is described as correction values, i.e., the change in seawater pH caused by the indicator. The diagrams are based on calculations made by using the equilibrium speciation programme, MARINHALT. From these calculations, and least squares fitting methods, pH correction values are described in terms of the pH difference between each seawater sample and the pH of an indicator stock solution. Calculations are performed for a typical high latitude water and a north Pacific deep water. Diagrams are presented for a salinity of 35 and a temperature of 15°C. Responses to salinities between 32 and 36 and temperatures 15–25°C are illustrated as well. A ±0.05 pH difference between a seawater sample and an indicator stock solution gives a correction of less than 0.001 pH unit for a 1-cm cell. For a 5-cm cell, pH differences between the indicator stock solution and a seawater sample as large as ±0.3 cause corrections smaller than ±0.001 pH unit. Calculations demonstrate that the five-fold lower indicator concentration used with 5-cm cells decreases the perturbation effect by approximately a factor of five relative to 1-cm cells.     

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.     

Long-Term Evaluation of Mulch Biowall Performance to Treat Chlorinated Solvents

Permeable reactive barriers (PRBs), such as mulch biowalls, have been installed at numerous groundwater cleanup sites, and laboratory and field studies have demonstrated biotic and abiotic processes that degrade chlorinated volatile organic compounds (CVOCs) in groundwater passing through these engineered remedies. However, the longevity of mulch biowalls remains a fundamental research question. Soil and groundwater sampling at seven mulch biowalls at Altus Air Force Base (AFB) approximately 10 years after installation demonstrated the ongoing degradation of CVOCs. Trichloroethene was not detected in five of seven groundwater samples collected from the biowall despite upgradient detections above federal drinking water standards. Microbial sampling established the presence of key dechlorinating bacteria and the abundance of genes encoding specific enzymes for degradation, high methane concentrations, low sulfate concentrations, and negative oxidation-reduction potential, all indicative of highly reducing conditions within the biowalls and favorable conditions for CVOC destruction via microbial reductive dechlorination. High cellulose content (>79%) of the mulch, elevated total organic carbon (TOC) content in groundwater, and elevated potentially bioavailable organic carbon (PBOC) measurements in soil samples further supports an ongoing, long-lived source of carbon. These results demonstrate the ongoing and long-term efficacy of the mulch biowalls at Altus AFB. In addition, concentrations of bacteria, TOC, PBOC, and other geochemical parameters suggest a modest impact of the biowalls downgradient. The continued presence of CVOCs downgradient may be attributable to back diffusion from low-permeability shale. However, the biowalls continue to provide benefits by removing CVOCs in groundwater, thus reducing further CVOC loading to the downgradient, low-permeability strata.     

Field Demonstration of the Horizontal Treatment Well (HRX Well®) for Passive In Situ Remediation

The horizontal reactive media treatment well (HRX Well®) uses directionally drilled horizontal wells filled with a treatment media to induce flow-focusing behavior created by the well-to-aquifer permeability contrast to passively capture proportionally large volumes of groundwater. Groundwater is treated in situ as it flows through the HRX Well and downgradient portions of the aquifer are cleaned via elution as these zones are flushed with clean water discharging from the HRX Well. The HRX Well concept is particularly well suited for sites where long-term mass discharge control is a primary performance objective. This concept is appropriate for recalcitrant and difficult-to-treat constituents, including chlorinated solvents, per- and polyfluoroalkyl substances (PFAS), 1,4-dioxane, and metals. A full-scale HRX Well was installed and operated to treat trichloroethene (TCE) with zero valent iron (ZVI). The model-predicted enhanced flow through the HRX Well (compared to the flow in and equivalent cross-sectional area orthogonal to flow in the natural formation before HRX Well installation) and treatment zone width was consistent with flows and widths estimated independently by point velocity probe (PVP) testing, HRX Well tracer testing, and observed treatment in downgradient monitoring wells. The actual average capture zone width was estimated to be between 45 and 69 feet. Total TCE mass discharge reduction was maintained through the duration of the performance monitoring period and exceeded 99.99% (%). Decreases in TCE concentrations were observed at all four downgradient monitoring wells within the treatment zone (ranging from 50 to 74% at day 436), and the first arrival of treated water was consistent with model predictions. The field demonstration confirmed the HRX Well technology is best suited for long-term mass discharge control, can be installed under active infrastructure, requires limited ongoing operation and maintenance, and has low life cycle energy and water requirements.