Located in the Sunbelt of the Southwestern United States, Phoenix Arizona finds itself in one of the hottest, driest places in the world. Thankfully, Phoenix has the Salt River, Gila River, Verde River, and a vast aquifer to meet the water demands of the municipal, industrial, and agricultural sectors. However, rampant groundwater pumping and over-allocation of these water supplies based on unprecedented, high flows of the Colorado River have created challenges for water managers to ensure adequate water supply for the future. Combined with the current 17-year drought and the warming and drying projections of climate change, the future of water availability in Phoenix will depend on the strength of water management laws, educating the public, developing a strong sense of community, and using development to manage population and support sustainability. As the prevalence of agriculture declines in and around Phoenix, a substantial amount of water is saved. Instead of storing this saved water, Phoenix is using it to support further development. Despite uncertainty regarding the abundant and continuous availability of Phoenix's water resources, development has hardly slowed and barely shifted directions to support sustainability. Phoenix was made to grow until it legally cannot expand anymore. In order to develop solutions, we must first understand the push for development in water-stressed Phoenix, Arizona.

There is a need to understand spatio-temporal variation of slip in active fault zones, both for the advancement of physics-based earthquake simulation and for improved probabilistic seismic hazard assessments. One challenge in the study of seismic hazards is producing a viable earthquake rupture forecast—a model that specifies the expected frequency and magnitude of events for a fault system. Time-independent earthquake forecasts can produce a mismatch among observed earthquake recurrence intervals, slip-per-event estimates, and implied slip rates. In this thesis, I developed an approach to refine several key geologic inputs to rupture forecasts by focusing on the San Andreas Fault in the Carrizo Plain, California. I use topographic forms, sub-surface excavations, and high-precision geochronology to understand the generation and preservation of slip markers at several spatial and temporal scales—from offset in a single earthquake to offset accumulated over thousands of years. This work results in a comparison of slip rate estimates in the Carrizo Plain for the last ~15 kyr that reduces ambiguity and enriches rupture forecast parameters. I analyzed a catalog of slip measurements and surveyed earth scientists with varying amounts of experience to validate high-resolution topography as a supplement to field-based active fault studies. The investigation revealed that (for both field and remote studies) epistemic uncertainties associated with measuring offset landforms can present greater limitations than the aleatoric limitations of the measurement process itself. I pursued the age and origin of small-scale fault-offset fluvial features at Van Matre Ranch, where topographic depressions were previously interpreted as single-event tectonic offsets. I provide new estimates of slip in the most recent earthquake, refine the centennial-scale fault slip rate, and formulate a new understanding of the formation of small-scale fault-offset fluvial channels from small catchments (<7,000 m2). At Phelan Creeks, I confirm the constancy of strain release for the ~15,000 years in the Carrizo Plain by reconstructing a multistage offset landform evolutionary history. I update and explicate a simplified model to interpret the geomorphic response of stream channels to strike-slip faulting. Lastly, I re-excavate and re-interpret paleoseismic catalogs along an intra-continental strike-slip fault (Altyn Tagh, China) to assess consistency of earthquake recurrence.

For this dissertation, three separate papers explore the study areas of the western Grand Canyon, the Grand Staircase (as related to Grand Canyon) and Desolation Canyon on the Green River in Utah.

In western Grand Canyon, I use comparative geomorphology between the Grand Canyon and the Grand Wash Cliffs (GWC). We propose the onset of erosion of the GWC is caused by slip on the Grand Wash Fault that formed between 18 and 12 million years ago. Hillslope angle and channel steepness are higher in Grand Canyon than along the Grand Wash Cliffs despite similar rock types, climate and base level fall magnitude. These experimental controls allow inference that the Grand Canyon is younger and eroding at a faster rate than the Grand Wash Cliffs.

The Grand Staircase is the headwaters of some of the streams that flow into Grand Canyon. A space-for-time substitution of erosion rates, supported by landscape simulations, implies that the Grand Canyon is the result of an increase in base level fall rate, with the older, slower base level fall rate preserved in the Grand Staircase. Our data and analyses also support a younger, ~6-million-year estimate of the age of Grand Canyon that is likely related to the integration of the Colorado River from the Colorado Plateau to the Basin and Range. Complicated cliff-band erosion and its effect on cosmogenic erosion rates are also explored, guiding interpretation of isotopic data in landscapes with stratigraphic variation in quartz and rock strength.

Several hypotheses for the erosion of Desolation Canyon are tested and refuted, leaving one plausible conclusion. I infer that the Uinta Basin north of Desolation Canyon is eroding slowly and that its form represents a slow, stable base level fall rate. Downstream of Desolation Canyon, the Colorado River is inferred to have established itself in the exhumed region of Canyonlands and to have incised to near modern depths prior to the integration of the Green River and the production of relief in Desolation Canyon. Analysis of incision and erosion rates in the region suggests integration is relatively recent.

The morphology of mountainous areas is strongly influenced by stream bed incision rates, but most studies of landscape evolution consider erosion at basin scales or larger. The research here attempts to understand the smaller-scale mechanics of erosion on exposed bedrock channels in the conceptual framework of an established saltation-abrasion model by Sklar and Dietrich [2004]. The recirculating flume used in this experiment allows independent control of bed slope, water discharge rate, sediment flux, and sediment grain size – all factors often bundled together in simple models of river incision and typically cross-correlated in natural settings. This study investigates the mechanics of erosion on exposed bedrock channels caused by abrasion of transported particles. Of particular interest are saltating particles, as well as sediment near the threshold between saltation and suspension - sediment vigorously transported but with significant interaction with the bed. The size of these erosive tools are varied over an order of magnitude in mean grain diameter, including a sand of D¬50 = 0.56 mm, and three gravel sizes of 3.39, 4.63, and 5.88 mm. Special consideration was taken to prevent any flow conditions that created a persistent alluvial cover. The erodible concrete substrate is fully exposed at all times during experiments reported here. Rates of erosion into the concrete substrate (a bedrock proxy) were measured by comparing topographic data before and after each experimental run, made possible by a precision laser mounted on a high speed computer-controlled cart. The experimental flume was able to produce flow discharge as high as 75 liters per second, sediment fluxes (of many varieties) up to 215 grams per second, and bed slopes up to 10%. I find a general positive correlation is found between erosion rate and bed slope, shear stress, grain size, and sediment flux.

ABSTRACT

The Spirit landing site in Gusev Crater has been imaged by the Mars Reconnaissance Orbiter High Resolution Imaging Science Experiment (HiRISE) camera more than thirty times since 2006. The breadth of this image set allowed a study of changes to surface features, covering four Mars years.

Small fields of bedforms comprised of dark material, and dark dust devil tracks are among the features revealed in the images. The bedforms are constrained within craters on the plains, and unconstrained in depressions less than 200m wide within the topography of the Columbia Hills, a ~120m-high structure in center of Gusev. Dust devil tracks appear in many images of the bedforms.

Within the Columbia Hills, three bedform fields approximately 180m2 and composed of fine dark basaltic sand were studied, using five HiRISE images taken from 2006 to 2014. Both bedform crests and the dust devil tracks superimposed on them were evaluated for change to azimuth and length, and for correlation between the features. The linear to slightly sinuous transverse crests ranging from less than 1m to 113m in length and two to three meters in wavelength, are primary bedforms. During the study they shifted as much as 33 degrees in azimuth, and individual crests moved on the surface as much as 0.75m. The greatest changes corresponded to a global dust storm in 2007. Average crest movement was documented at the rate of 0.25m per year. Rather than moving progressively, the crests eventually returned to near their original orientation after the storm. The dust devil tracks, reflecting a more complex wind regime, including vortex development during diurnal heating, maintained predominantly NW-SE orientations but also reflected the effects of the storm.

The observed modifications were neither progressive, nor strictly seasonal. The apparent stability of the bedform geometry over four seasons supports the predictions of the Mars Regional Atmospheric Modeling System (MRAMS): low speed (1-7.5 ms-1), daily alternating winds of relatively equal force. Crest profiles were found to be nearly symmetrical, without slipfaces to indicate a preferential wind direction; this finding also is supported by the MRAMS model.

On Mars, sedimentary deposits reveal a complex history of water- and wind-related geologic processes. Central mounds – kilometer-scale stacks of sediment located within craters – occur across Mars, but the specific processes responsible for mound formation and subsequent modification are still uncertain. A survey of central mounds within large craters was conducted. Mound locations, mound offsets within their host craters, and relative mound heights were used to address various mound formation hypotheses. The results suggest that mound sediments once filled their host craters and were later eroded into the features observed today. Mounds offsets from the center of their host crater imply that wind caused the erosion of central mounds. An in depth study of a single central mound (Mt. Sharp within Gale crater) was also conducted. Thermal Emission Imaging System Visible Imaging Subsystem (THEMIS-VIS) mosaics in grayscale and false color were used to characterize the morphology and color variations in and around Gale crater. One result of this study is that dunes within Gale crater vary in false color composites from blue to purple, and that these color differences may be due to changes in dust cover, grain size, and/or composition. To further investigate dune fields on Mars, albedo variations at eight dune fields were studied based on the hypothesis that a dune’s ripple migration rate is correlated to its albedo. This study concluded that a dune’s minimum albedo does not have a simple correlation with its ripple migration rate. Instead, dust devils remove dust on slow-moving and immobile dunes, whereas saltating sand caused by strong winds removes dust on faster-moving dunes.

On the Moon, explosive volcanic deposits within Oppenheimer crater that were emplaced ballistically were investigated. Lunar Reconnaissance Orbiter (LRO) Diviner Radiometer mid-infrared data, LRO Camera images, and Chandrayaan-1 orbiter Moon Mineralogy Mapper near-infrared spectra were used to test the hypothesis that the pyroclastic deposits in Oppenheimer crater were emplaced via Vulcanian activity by constraining their composition and mineralogy. The mineralogy and iron-content of the pyroclastic deposits vary significantly (including examples of potentially very high iron compositions), which indicates variability in eruption style. These results suggest that localized lunar pyroclastic deposits may have a more complex origin and mode of emplacement than previously thought.

I investigate the Black Canyon City landslide (BCC landslide), a prominent deep-seated landslide located northeast of Black Canyon City, Arizona. Although the landslide does not appear to pose a significant hazard to structures, its prominent features and high topographic relief make it an excellent site to study the geologic setting under which such features develop. This study has the potential to contribute toward understanding the landscape evolution in similar geologic and topographic settings, and for characterizing the underlying structural processes of this deep-seated feature. We use field and remotely-based surface geology and geomorphological mapping to characterize the landslide geometry and its surface displacement. We use the Structure from Motion (SfM) method to generate a 0.2 m resolution digital elevation model and rectified ortho-photo imagery from unmanned aerial vehicle (UAV) - and balloon-based images and used them as the base map for our mapping. The ~0.6 km2 landslide is easily identified through remotely-sensed imagery and in the field because of the prominent east-west trending fractures defining its upper extensional portion. The landslide displaces a series of Early and Middle Miocene volcanic and sedimentary rocks. The main head scarp is ~600 m long and oriented E-W with some NW-SE oriented minor scarps. Numerous fractures varying from millimeters to meters in opening were identified throughout the landslide body (mostly with longitudinal orientation). The occurrence of a distinctive layer of dark reddish basalt presents a key displaced marker to estimate the long-term deformation of the slide mass. Using this marker, the total vertical displacement is estimated to be ~70 m, with maximum movement of ~95 m to the SE. This study indicates that the landslide motion is translational with a slight rotational character. We estimate the rate of the slide motion by resurvey of monuments on and off the slide, and examination of disturbed vegetation located along the fractures. The analysis indicates a slow integrated average landslide velocity of 10-60 mm/yr. The slide motion is probably driven during annual wet periods when increased saturation of the slide mass weakens the basal slip surface and the overall mass of the slide is increased. Results from our study suggest that the slide is stable and does not pose significant hazard for the surrounding area given no extreme changes in the environmental condition. Although the landslide is categorized as very slow (according to Cruden and Varnes, 1996), monitoring the landslide is still necessary.

Shallow earthquakes in the upper part of the overriding plate of subduction zones can be devastating due to their proximity to population centers despite the smaller rupture extents than commonly occur on subduction megathrusts that produce the largest earthquakes. Damaging effects can be greater in volcanic arcs like Java because ground shaking is amplified by surficial deposits of uncompacted volcaniclastic sediments. Identifying the upper-plate structures and their potential hazards is key for minimizing the dangers they pose. In particular, the knowledge of the regional stress field and deformation pattern in this region will help us to better understand how subduction and collision affects deformation in this part of the overriding plate. The majority of the upper plate deformation studies have been focused on the deformation in the main thrusts of the fore-arc region. Study of deformation within volcanic arc is limited despite the associated earthquake hazards. In this study, I use maps of active upper-plate structures, earthquake moment tensor data and stress orientation deduced from volcano morphology analysis to characterize the strain field of Java arc. In addition, I use sandbox analog modeling to evaluate the mechanical factors that may be important in controlling deformation. My field- and remotely-based mapping of active faults and folds, supplemented by results from my paleoseismic studies and physical models of the system, suggest that Java’s deformation is distributed over broad areas along small-scale structures. Java is segmented into three main zones based on their distinctive structural patterns and stress orientation. East Java is characterized by NW-SE normal and strike-slip faults, Central Java has E-W folds and thrust faults, and NE-SW strike-slip faults dominate West Java. The sandbox analog models indicate that the strain in response to collision is partitioned into thrusting and strike-slip faulting, with the dominance of margin-normal thrust faulting. My models test the effects of convergence obliquity, geometry, preexisting weaknesses, asperities, and lateral strength contrast. The result suggest that slight variations in convergence obliquity do not affect the deformation pattern significantly, while the margin shape, lateral strength contrast, and perturbation of deformation from asperities each have a greater impact on deformation.

The collision between the Indian and Eurasian tectonic plates marked the onset of the rise of the Himalayan-Tibetan orogen, but also brought about profound changes to the Earth's oceans and climate. The exact sequence of events that occurred during this collision is poorly understood, leading to a wide range of estimates of its age. The Indus and Yarlung sutures are generally considered to represent the final collision between India and Eurasia, and together form a mostly continuous belt that can be traced over 2000 km along strike. In the western portions of the orogen the Karakoram Fault introduces a key complexity to the study of timing of collision by offsetting the Indus and Yarlung sutures. Recent work has used the complexities introduced by the Karakoram Fault to suggest that the more northerly Shyok suture, not the Indus suture, represents the India-Eurasia collision zone. Estimates for timing of the India-Eurasia collision fall into one of three groups: 40-34 Ma, 55-50 Ma, and 66-60 Ma. Attempts to reconcile these models have thus far been unsuccessful. In order to provide additional data that might further clarify the timing and location of collision, studies have been performed along the Shyok suture in India and along the Yarlung suture in Tibet at Sangsang. A study along the Shyok suture argues that the suture formed between 92-85 Ma. This timing precludes an interpretation that the Shyok suture marks the location of the India-Eurasia collision. A second study demonstrates the utility of two new geochronometers, (U-Th)/Pb joaquinite and 40Ar/39Ar neptunite, that play an important role in unraveling the tectonic history of the Yarlung suture. A third study is an investigation of the structure and geochronology of the Sangsang ophiolite complex. Here, multiple (U-Th)/Pb and 40Ar/39Ar systems record magmatism and metamorphism spanning ca. 125-52 Ma. By tying these chronometers to tectonic process, a history is reconstructed of the southern margin of Tibet that includes Early Cretaceous to Late Cretaceous forearc rifting associated with mid ocean ridge subduction, Paleocene accretionary wedge uplift and erosion, and finally Eocene metasomatism and collision.