Chemical and physical interactions of flowing ice and rock have inexorably shaped planetary surfaces. Weathering in glacial environments is a significant link in biogeochemical cycles – carbon and strontium – on Earth, and may have once played an important role in altering Mars’ surface. Despite growing recognition of the importance of low-temperature chemical weathering, these processes are still not well understood. Debris-coated glaciers are also present on Mars, emphasizing the need to study ice-related processes in the evolution of planetary surfaces. During Earth’s history, subglacial environments are thought to have sheltered communities of microorganisms from extreme climate variations. On Amazonian Mars, glaciers such as lobate debris aprons (LDA) could have hosted chemolithotrophic communities, making Mars’ present glaciers candidates for life preservation. This study characterizes glacial processes on both Earth and Mars.

Chemical weathering at Robertson Glacier, a small alpine glacier in the Canadian Rocky Mountains, is examined with a multidisciplinary approach. The relative proportions of differing dissolution reactions at various stages in the glacial system are empirically determined using aqueous geochemistry. Synthesis of laboratory and orbital thermal infrared spectroscopy allows identification of dissolution rinds on hand samples and characterization of carbonate dissolution signals at orbital scales, while chemical and morphological evidence for thin, discontinuous weathering rinds at microscales are evident from electron microscopy. Subglacial dissolution rates are found to outpace those of the proglacial till plain; biologically-mediated pyrite oxidation drives the bulk of this acidic weathering.

Second, the area-elevation relationship, or hypsometry, of LDA in the midlatitudes of Mars is characterized. These glaciers are believed to have formed ~500 Ma during a climate excursion. Hypsometric measurements of these debris-covered glaciers enable insight into past flow regimes and drive predictions about past climate scenarios. The LDA in this study fall into three major groups, strongly dependent on basal elevation, implying regional and climatic controls on ice formation and flow.

I show that biologically-mediated mineral reactions drive high subglacial dissolution rates, such that variations within the valley can be detected with remote sensing techniques. In future work, these insights can be applied to examining Mars’ glacial regions for signs of chemical alteration and biosignatures.

Sedimentary basins are defined by extensional tectonics. Rugged mountain ranges stand in stark relief adjacent to muted structural basins filled with sediment. In simplest terms, this topography is the result of ranges uplifted along normal faults, and this uplift drives erosion within upland drainages, shedding sediment into subsiding basins. In southeastern Arizona's Basin and Range province extensional tectonics waned at approximately 3-5 Myr, and the region's structural basins began transitioning from internal to external drainage, forming the modern Gila River fluvial network. In the Atacama Desert of northern Chile, some basins of the Central Depression remain internally drained while others have integrated to the Pacific Ocean. In northern Chile, rates of landscape evolution are some of the slowest on Earth due to the region's hyperarid climate. While the magnitude of upland erosion driven by extensional tectonics is largely recorded in the stratigraphy of the structural basins, the landscape's response to post-tectonic forcings is unknown.

I employ the full suite of modern geomorphic tools provided by terrestrial cosmogenic nuclides - surface exposure dating, conventional burial dating, isochron burial dating, quantifying millennial-scale upland erosion rates using detrital TCN, quantifying paleo-erosion rates using multiple TCN such as Ne-21/Be-10 and Al-26l/Be-10, and assessing sediment recycling and complex exposure using multiple TCN - to quantify the rates of landscape evolution in southeastern Arizona and northern Chile during the Late Cenozoic. In Arizona, I also use modern remnants of the pre-incision landscape and digital terrain analyses to reconstruct the landscape, allowing the quantification of incision and erosion rates that supplement detrital TCN-derived erosion rates. A new chronology for key basin high stand remnants (Frye Mesa) and a flight of Gila River terraces in Safford basin provides a record of incision rates from the Pliocene through the Quaternary, and I assess how significantly regional incision is driving erosion rates. Paired nuclide analyses in the Atacama Desert of northern Chile reveal complex exposure histories resulting from several rounds of transport and burial by fluvial systems. These results support a growing understanding that geomorphic processes in the Atacama Desert are more active than previously thought despite the region's hyperarid climate.

Earth's topographic surface forms an interface across which the geodynamic and geomorphic engines interact. This interaction is best observed along crustal margins where topography is created by active faulting and sculpted by geomorphic processes. Crustal deformation manifests as earthquakes at centennial to millennial timescales. Given that nearly half of Earth's human population lives along active fault zones, a quantitative understanding of the mechanics of earthquakes and faulting is necessary to build accurate earthquake forecasts. My research relies on the quantitative documentation of the geomorphic expression of large earthquakes and the physical processes that control their spatiotemporal distributions. The first part of my research uses high-resolution topographic lidar data to quantitatively document the geomorphic expression of historic and prehistoric large earthquakes. Lidar data allow for enhanced visualization and reconstruction of structures and stratigraphy exposed by paleoseismic trenches. Lidar surveys of fault scarps formed by the 1992 Landers earthquake document the centimeter-scale erosional landforms developed by repeated winter storm-driven erosion. The second part of my research employs a quasi-static numerical earthquake simulator to explore the effects of fault roughness, friction, and structural complexities on earthquake-generated deformation. My experiments show that fault roughness plays a critical role in determining fault-to-fault rupture jumping probabilities. These results corroborate the accepted 3-5 km rupture jumping distance for smooth faults. However, my simulations show that the rupture jumping threshold distance is highly variable for rough faults due to heterogeneous elastic strain energies. Furthermore, fault roughness controls spatiotemporal variations in slip rates such that rough faults exhibit lower slip rates relative to their smooth counterparts. The central implication of these results lies in guiding the interpretation of paleoseismically derived slip rates that are used to form earthquake forecasts. The final part of my research evaluates a set of Earth science-themed lesson plans that I designed for elementary-level learning-disabled students. My findings show that a combination of concept delivery techniques is most effective for learning-disabled students and should incorporate interactive slide presentations, tactile manipulatives, teacher-assisted concept sketches, and student-led teaching to help learning-disabled students grasp Earth science concepts.

The southern Tibetan Plateau margin between ~ 83E and 86.5E is defined by an abrupt change from the low-relief Tibetan Plateau to the rugged topography and deep gorges of the Himalaya. This physiographic transition lies well to the north of active thrusting, and thus, the mechanism responsible for the distinct topographic break remains the focus of much debate. While numerous studies have utilized thermochronology to examine the exhumation history of the Himalaya, few have done so with respect to variations across the Himalaya-Tibetan Plateau transition. In this work, we examine the nature of the transition where it is accessible and well-defined in the Nyalam valley of south-central Tibet. We employ several new and previously published thermochronologic datasets (with a closure temperature range of ~ 70C–300C) in conjunction with river incision patterns inferred by the longitudinal profile of the Bhote Kosi River.

The results reveal a sharp change in cooling rate at ~ 3.5 Ma at a location corresponding to a pronounced river knickpoint representing a sharp increase in river gradient and presumably incision rate (a proxy for rock uplift). Margin retreat models for the physiographic transition are inconsistent with the cooling pattern revealed by low-temperature thermochronologic data. Models invoking passive uplift of the upper crust over a midcrustal ramp and associated duplex to account for the physiographic transition do not explain the observed break in cooling rate there, although they may explain a suggesting in the thermochronologic data of progressively increasing exhumation rates south of the transition. The simplest model consistent with all observations is that passive uplift is augmented by contemporaneous differential uplift across a young (Pliocene-Quaternary) normal fault at the physiographic transition. Drawing on observations elsewhere, we hypothesize that similar structural relationships may be characteristic of the Tibetan Plateau-Himalaya transition from ~83E – 86.5E.

Many shallow craters near the Spirit Mars Exploration Rover landing site contain asymmetric deposits of windblown sediments which could indicate the predominant local wind direction at the time of deposition or redistribution. Wind tunnel simulations and field studies of terrestrial craters were used to determine trends in deposition as a function of crater morphometry and wind direction. Terrestrial analog field work at the Amboy lava field, Mojave Desert, California, included real-time wind measurements and assessments of active sediment deposition in four small (<100 m) craters. Preliminary results indicate that reverse flow or stagnant wind and deposition on the upwind side of the crater floor occurs in craters with depth-to-diameter (d/D) ratios ≥0.05. Measurements taken within a crater of d/D of ~0.02 do not indicate reverse flow. Therefore, reverse flow is expected to cease within a d/D range of 0.02 to 0.05, resulting in wind movement directly over the crater floor in the downwind direction with no asymmetric sediment deposition. Wind tunnel simulations using six crater models, including a scaled model of a crater from the Amboy lava field, were completed to assess the wind flow in and around craters as a function of crater morphometry (depth, diameter). Reverse flow occurred in craters with d/D ratios ≥0.033, resulting in sediment deposition in the upwind portion of the crater floor. Visual observations of a crater with a d/D of ~0.020 did not indicate reverse flow, similar to the results of field studies; therefore, reverse flow appears to cease within a d/D range of 0.020 to 0.033. Craters with asymmetric aeolian deposits near the Mars Spirit landing site have d/D ratios of 0.034 to 0.076, suggesting that reverse flow occurs in these craters. Thus, the position of windblown sediments in the northwest parts of the crater floors would indicate prevailing winds from the northwest to the southeast, consistent with late afternoon winds as predicted by the Mars Regional Atmospheric Modeling System (MRAMS) circulation model.