Anticipating responses to climate change and planning for resilience in California's freshwater ecosystems.

As human-caused climate changes accelerate, California will experience hydrologic and temperature conditions different than any encountered in recorded history. How will these changes affect the state's freshwater ecosystems? Rivers, lakes, and wetlands are managed as a water resource, but they also support a complex web of life, ranging from bacteria, fungi, and algae to macrophytes, woody plants, invertebrates, fish, amphibians, reptiles, birds, and mammals. In much of the state, native freshwater organisms already struggle to survive massive water diversions and dams, deteriorating water quality, extensive land cover modification for agriculture and urban development, and invasions of exotic species. In the face of climate change, we need to expand efforts to recover degraded ecosystems and to protect the resilience, health, and viability of existing ecosystems. For this, more process-based understanding of river, lake, and wetlands ecosystems is needed to forecast how systems will respond to future climate change and to our interventions. This will require 1) expanding our ability to model mechanistically how freshwater biota and ecosystems respond to environmental change; 2) hypothesis-driven monitoring and field studies; 3) education and training to build research, practitioner, stewardship, and policy capabilities; and 4) developing tools and policies for building resilient ecosystems. A goals-driven, hypothesis-informed collaboration among tribes, state (and federal) agencies, nongovernmental organizations, academicians, and consultants is needed to accomplish these goals and to advance the skills and knowledge of the future workforce of practitioners, regulators, and researchers who must live with the climate changes that are already upon us and will intensify.

Controls on the size distributions of shallow landslides.

Rainfall-triggered shallow landslides are destructive hazards and play an important role in landscape processes. A theory explaining the size distributions of such features remains elusive. Prior work connects size distributions to topography, but field-mapped inventories reveal pronounced similarities in the form, mode, and spread of distributions from diverse landscapes. We analyze nearly identical distributions occurring in the Oregon Coast Range and the English Lake District, two regions of strikingly different topography, lithology, and vegetation. Similarity in minimum sizes at these sites is partly explained by theory that accounts for the interplay of mechanical soil strength controls resisting failure. Maximum sizes, however, are not explained by current theory. We develop a generalized framework to account for the entire size distribution by unifying a mechanistic slope stability model with a flexible spatial-statistical description for the variability of hillslope strength. Using hillslope-scale numerical experiments, we find that landslides can occur not only in individual low strength areas but also across multiple smaller patches that coalesce. We show that reproducing observed size distributions requires spatial strength variations to be strongly localized, of large amplitude, and a consequence of multiple interacting factors. Such constraints can act together with the mechanical determinants of landslide initiation to produce size distributions of broadly similar character in widely different landscapes, as found in our examples. We propose that size distributions reflect the systematic scale dependence of the spatially averaged strength. Our results highlight the critical need to constrain the form, amplitude, and wavelength of spatial variability in material strength properties of hillslopes.

Direct observations of rock moisture, a hidden component of the hydrologic cycle.

Recent theory and field observations suggest that a systematically varying weathering zone, that can be tens of meters thick, commonly develops in the bedrock underlying hillslopes. Weathering turns otherwise poorly conductive bedrock into a dynamic water storage reservoir. Infiltrating precipitation typically will pass through unsaturated weathered bedrock before reaching groundwater and running off to streams. This invisible and difficult to access unsaturated zone is virtually unexplored compared with the surface soil mantle. We have proposed the term "rock moisture" to describe the exchangeable water stored in the unsaturated zone in weathered bedrock, purposely choosing a term parallel to, but distinct from, soil moisture, because weathered bedrock is a distinctly different material that is distributed across landscapes independently of soil thickness. Here, we report a multiyear intensive campaign of quantifying rock moisture across a hillslope underlain by a thick weathered bedrock zone using repeat neutron probe measurements in a suite of boreholes. Rock moisture storage accumulates in the wet season, reaches a characteristic upper value, and rapidly passes any additional rainfall downward to groundwater. Hence, rock moisture storage mediates the initiation and magnitude of recharge and runoff. In the dry season, rock moisture storage is gradually depleted by trees for transpiration, leading to a common lower value at the end of the dry season. Up to 27% of the annual rainfall is seasonally stored as rock moisture. Significant rock moisture storage is likely common, and yet it is missing from hydrologic and land-surface models used to predict regional and global climate.

Inverted channel variations identified on a distal portion of a bajada in the central Atacama Desert, Chile.

In deserts, the interplay between occasional fluvial events and persistent aeolian erosion can form composite modern and relict surfaces, especially on the distal portion of alluvial fans. There, relief inversion of alluvial deposits by differential erosion can form longitudinal ridges. We identified two distinct ridge types formed by relief inversion on converging alluvial fans in the hyperarid Chilean Atacama Desert. Although they are co-located and similar in scale, the ridge types have different ages and formation histories that apparently correspond to minor paleoclimate variations. Gravel-armored ridges are remnants of deflated alluvial deposits with a bimodal sediment distribution (gravel and sand) dated to a minor pluvial phase at the end of the Late Pleistocene (~12 kyr). In contrast, younger (~9 kyr) sulfate-capped ridges formed during a minor arid phase with evaporite deposition in a pre-existing channel that armored the underlying deposits. Collectively, inverted channels at Salar de Llamara resulted from multiple episodes of surface overland flow and standing water spanning several thousand years. Based on ridge relief and age, the minimum long-term deflation rate is 0.1-0.2 m/kyr, driven primarily by wind erosion. This case study is an example of the equifinality concept whereby different processes lead to similar landforms. The complex history of the two ridge types can only be generally constrained in remotely sensed data. In situ observations are required to discern the specifics of the aqueous history, including the flow type, magnitude, sequence, and paleoenvironment. These findings have relevance for interpreting similar landforms on Mars.

A bottom-up control on fresh-bedrock topography under landscapes.

The depth to unweathered bedrock beneath landscapes influences subsurface runoff paths, erosional processes, moisture availability to biota, and water flux to the atmosphere. Here we propose a quantitative model to predict the vertical extent of weathered rock underlying soil-mantled hillslopes. We hypothesize that once fresh bedrock, saturated with nearly stagnant fluid, is advected into the near surface through uplift and erosion, channel incision produces a lateral head gradient within the fresh bedrock inducing drainage toward the channel. Drainage of the fresh bedrock causes weathering through drying and permits the introduction of atmospheric and biotically controlled acids and oxidants such that the boundary between weathered and unweathered bedrock is set by the uppermost elevation of undrained fresh bedrock, Zb. The slow drainage of fresh bedrock exerts a "bottom up" control on the advance of the weathering front. The thickness of the weathered zone is calculated as the difference between the predicted topographic surface profile (driven by erosion) and the predicted groundwater profile (driven by drainage of fresh bedrock). For the steady-state, soil-mantled case, a coupled analytical solution arises in which both profiles are driven by channel incision. The model predicts a thickening of the weathered zone upslope and, consequently, a progressive upslope increase in the residence time of bedrock in the weathered zone. Two nondimensional numbers corresponding to the mean hillslope gradient and mean groundwater-table gradient emerge and their ratio defines the proportion of the hillslope relief that is unweathered. Field data from three field sites are consistent with model predictions.

The search for a topographic signature of life.

Landscapes are shaped by the uplift, deformation and breakdown of bedrock and the erosion, transport and deposition of sediment. Life is important in all of these processes. Over short timescales, the impact of life is quite apparent: rock weathering, soil formation and erosion, slope stability and river dynamics are directly influenced by biotic processes that mediate chemical reactions, dilate soil, disrupt the ground surface and add strength with a weave of roots. Over geologic time, biotic effects are less obvious but equally important: biota affect climate, and climatic conditions dictate the mechanisms and rates of erosion that control topographic evolution. Apart from the obvious influence of humans, does the resulting landscape bear an unmistakable stamp of life? The influence of life on topography is a topic that has remained largely unexplored. Erosion laws that explicitly include biotic effects are needed to explore how intrinsically small-scale biotic processes can influence the form of entire landscapes, and to determine whether these processes create a distinctive topography.

Experimental evidence for the conditions necessary to sustain meandering in coarse-bedded rivers.

Meandering rivers are common on Earth and other planetary surfaces, yet the conditions necessary to maintain meandering channels are unclear. As a consequence, self-maintaining meandering channels with cutoffs have not been reproduced in the laboratory. Such experimental channels are needed to explore mechanisms controlling migration rate, sinuosity, floodplain formation, and planform morphodynamics and to test theories for wavelength and bend propagation. Here we report an experiment in which meandering with near-constant width was maintained during repeated cutoff and regeneration of meander bends. We found that elevated bank strength (provided by alfalfa sprouts) relative to the cohesionless bed material and the blocking of troughs (chutes) in the lee of point bars via suspended sediment deposition were the necessary ingredients to successful meandering. Varying flood discharge was not necessary. Scaling analysis shows that the experimental meander migration was fast compared to most natural channels. This high migration rate caused nearly all of the bedload sediment to exchange laterally, such that bar growth was primarily dependent on bank sediment supplied from upstream lateral migration. The high migration rate may have contributed to the relatively low sinuosity of 1.19, and this suggests that to obtain much higher sinuosity experiments at this scale may have to be conducted for several years. Although patience is required to evolve them, these experimental channels offer the opportunity to explore several fundamental issues about river morphodynamics. Our results also suggest that sand supply may be an essential control in restoring self-maintaining, actively shifting gravel-bedded meanders.

Orbital and In-Situ Investigation of Periodic Bedrock Ridges in Glen Torridon, Gale Crater, Mars.

Gale crater, the field site for NASA's Mars Science Laboratory Curiosity rover, contains a diverse and extensive record of aeolian deposition and erosion. This study focuses on a series of regularly spaced, curvilinear, and sometimes branching bedrock ridges that occur within the Glen Torridon region on the lower northwest flank of Aeolis Mons, the central mound within Gale crater. During Curiosity's exploration of Glen Torridon between sols ∼2300-3080, the rover drove through this field of ridges, providing the opportunity for in situ observation of these features. This study uses orbiter and rover data to characterize ridge morphology, spatial distribution, compositional and material properties, and association with other aeolian features in the area. Based on these observations, we find that the Glen Torridon ridges are consistent with an origin as wind-eroded bedrock ridges, carved during the exhumation of Mount Sharp. Erosional features like the Glen Torridon ridges observed elsewhere on Mars, termed periodic bedrock ridges (PBRs), have been interpreted to form transverse to the dominant wind direction. The size and morphology of the Glen Torridon PBRs are consistent with transverse formative winds, but the orientation of nearby aeolian bedforms and bedrock erosional features raise the possibility of PBR formation by a net northeasterly wind regime. Although several formation models for the Glen Torridon PBRs are still under consideration, and questions persist about the nature of PBR-forming paleowinds, the presence of PBRs at this site provides important constraints on the depositional and erosional history of Gale crater.

Origin and composition of three heterolithic boulder- and cobble-bearing deposits overlying the Murray and Stimson formations, Gale Crater, Mars.

Heterolithic, boulder-containing, pebble-strewn surfaces occur along the lower slopes of Aeolis Mons ("Mt. Sharp") in Gale crater, Mars. They were observed in HiRISE images acquired from orbit prior to the landing of the Curiosity rover. The rover was used to investigate three of these units named Blackfoot, Brandberg, and Bimbe between sols 1099 and 1410. These unconsolidated units overlie the lower Murray formation that forms the base of Mt. Sharp, and consist of pebbles, cobbles and boulders. Blackfoot also overlies portions of the Stimson formation, which consists of eolian sandstone that is understood to significantly postdate the dominantly lacustrine deposition of the Murray formation. Blackfoot is elliptical in shape (62 × 26 m), while Brandberg is nearly circular (50 × 55 m), and Bimbe is irregular in shape, covering about ten times the area of the other two. The largest boulders are 1.5-2.5 m in size and are interpreted to be sandstones. As seen from orbit, some boulders are light-toned and others are dark-toned. Rover-based observations show that both have the same gray appearance from the ground and their apparently different albedos in orbital observations result from relatively flat sky-facing surfaces. Chemical observations show that two clasts of fine sandstone at Bimbe have similar compositions and morphologies to nine ChemCam targets observed early in the mission, near Yellowknife Bay, including the Bathurst Inlet outcrop, and to at least one target (Pyramid Hills, Sol 692) and possibly a cap rock unit just north of Hidden Valley, locations that are several kilometers apart in distance and tens of meters in elevation. These findings may suggest the earlier existence of draping strata, like the Stimson formation, that would have overlain the current surface from Bimbe to Yellowknife Bay. Compositionally these extinct strata could be related to the Siccar Point group to which the Stimson formation belongs. Dark, massive sandstone blocks at Bimbe are chemically distinct from blocks of similar morphology at Bradbury Rise, except for a single float block, Oscar (Sol 516). Conglomerates observed along a low, sinuous ridge at Bimbe consist of matrix and clasts with compositions similar to the Stimson formation, suggesting that stream beds likely existed nearly contemporaneously with the dunes that eventually formed the Stimson formation, or that they had the same source material. In either case, they represent a later pulse of fluvial activity relative to the lakes associated with the Murray formation. These three units may be local remnants of infilled impact craters (especially circular-shaped Brandberg), decayed buttes, patches of unconsolidated fluvial deposits, or residual mass-movement debris. Their incorporation of Stimson and Murray rocks, the lack of lithification, and appearance of being erosional remnants suggest that they record erosion and deposition events that post-date the exposure of the Stimson formation.

The Mars Science Laboratory (MSL) Mast cameras and Descent imager: Investigation and instrument descriptions.

The Mars Science Laboratory Mast camera and Descent Imager investigations were designed, built, and operated by Malin Space Science Systems of San Diego, CA. They share common electronics and focal plane designs but have different optics. There are two Mastcams of dissimilar focal length. The Mastcam-34 has an f/8, 34 mm focal length lens, and the M-100 an f/10, 100 mm focal length lens. The M-34 field of view is about 20° × 15° with an instantaneous field of view (IFOV) of 218 µrad; the M-100 field of view (FOV) is 6.8° × 5.1° with an IFOV of 74 µrad. The M-34 can focus from 0.5 m to infinity, and the M-100 from ~1.6 m to infinity. All three cameras can acquire color images through a Bayer color filter array, and the Mastcams can also acquire images through seven science filters. Images are ≤1600 pixels wide by 1200 pixels tall. The Mastcams, mounted on the ~2 m tall Remote Sensing Mast, have a 360° azimuth and ~180° elevation field of regard. Mars Descent Imager is fixed-mounted to the bottom left front side of the rover at ~66 cm above the surface. Its fixed focus lens is in focus from ~2 m to infinity, but out of focus at 66 cm. The f/3 lens has a FOV of ~70° by 52° across and along the direction of motion, with an IFOV of 0.76 mrad. All cameras can acquire video at 4 frames/second for full frames or 720p HD at 6 fps. Images can be processed using lossy Joint Photographic Experts Group and predictive lossless compression.

A multidimensional stability model for predicting shallow landslide size and shape across landscapes.

The size of a shallow landslide is a fundamental control on both its hazard and geomorphic importance. Existing models are either unable to predict landslide size or are computationally intensive such that they cannot practically be applied across landscapes. We derive a model appropriate for natural slopes that is capable of predicting shallow landslide size but simple enough to be applied over entire watersheds. It accounts for lateral resistance by representing the forces acting on each margin of potential landslides using earth pressure theory and by representing root reinforcement as an exponential function of soil depth. We test our model's ability to predict failure of an observed landslide where the relevant parameters are well constrained by field data. The model predicts failure for the observed scar geometry and finds that larger or smaller conformal shapes are more stable. Numerical experiments demonstrate that friction on the boundaries of a potential landslide increases considerably the magnitude of lateral reinforcement, relative to that due to root cohesion alone. We find that there is a critical depth in both cohesive and cohesionless soils, resulting in a minimum size for failure, which is consistent with observed size-frequency distributions. Furthermore, the differential resistance on the boundaries of a potential landslide is responsible for a critical landslide shape which is longer than it is wide, consistent with observed aspect ratios. Finally, our results show that minimum size increases as approximately the square of failure surface depth, consistent with observed landslide depth-area data.

Formation of Box Canyon, Idaho, by megaflood: implications for seepage erosion on Earth and Mars.

Amphitheater-headed canyons have been used as diagnostic indicators of erosion by groundwater seepage, which has important implications for landscape evolution on Earth and astrobiology on Mars. Of perhaps any canyon studied, Box Canyon, Idaho, most strongly meets the proposed morphologic criteria for groundwater sapping because it is incised into a basaltic plain with no drainage network upstream, and approximately 10 cubic meters per second of seepage emanates from its vertical headwall. However, sediment transport constraints, 4He and 14C dates, plunge pools, and scoured rock indicate that a megaflood (greater than 220 cubic meters per second) carved the canyon about 45,000 years ago. These results add to a growing recognition of Quaternary catastrophic flooding in the American northwest, and may imply that similar features on Mars also formed by floods rather than seepage erosion.