The presence of ices (H2O and CO2) and liquid water is key to the evolution ofmartian geology, with implications for the potential for past or extant life, and the future of robotic and human exploration on Mars. In this dissertation, I present the first direct evidence that the smooth deposits covering mid-latitude, martian pole-facing slopes are composed of shallow dusty H2O ice covered by desiccated material. To analyze this H2O ice, I developed the first validated radiative transfer model for dusty martian snow and glacier ice. I found that these ice exposures have < 1% dust in them, and discovered the lowest latitude detection of H2O ice on Mars, at 32.9°S. After observing the ice disappear, and new gully channels form, I proposed a model for gully formation. In this model, dusty ice gets exposed by slumping, leading to melting in the subsurface and channels eroding within the ice and the wall rock beneath. Access to liquid water within this ice could provide potential abodes for any extant life. Next, I developed novel methodology to search for CO2 frosts within the entire Thermal Emission Imaging System (THEMIS) infrared dataset and found that about half of all gullies overlap with CO2 frost detections. I also used the Thermal Emission Spectrometer (TES) water vapor retrievals to assess the formation and distribution of H2O frosts on Mars. Additionally, I used radar data from the Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) instrument to investigate Mars’ ice-rich South Polar Layered Deposits (SPLD). I discovered radar signals similar to those proposed to be caused by a subglacial lake throughout the martian SPLD. Finally, I mapped martian polygonal ridge networks thought to represent fossilized remnants of ancient groundwater near the Perseverance rover landing site with the help of citizen scientists across a fifth of Mars’ total surface area and analyzed their thermophysical properties. All these studies highlight the key role that ices and liquid water have played in shaping Mars’ landscape through time, and provide an intriguing path forward in martian exploration and the search for alien life.

Subsolidus convection in the mantle of Earth is the driving mechanism behind plate tectonics and provides a central framework linking geophysical, geochemical, petrological, hydrological, and biological processes within the system. Seismic observations have revealed mantle heterogeneities in wide-ranging scales from less than tens of to thousands of kilometers. Understanding the origins and dynamics of these anomalies is critical to advance our knowledge on how mantle convection operates and coevolves with the surface system. This dissertation attempts to constrain the past, present and future of mantle dynamics with lines of evidence from seismology, geodynamics, petrology, geochemistry, and astrophysics. Above Earth’s core, two continent-sized large low shear velocity provinces (LLSVPs) beneath Africa and the Pacific Ocean were seismically detected decades ago. Yet their origin, composition, detailed morphology and influence over mantle convection remain elusive. First, I propose the two LLSVPs may represent the mantle remnants of the Moon-forming impactor Theia. I show that the mantle of Theia is intrinsically denser than Earth’s mantle and would have sunk and accumulated into LLSVP-like structures in the deepest mantle after 4.5 billion years. Second, I examined the maximum height of the two LLSVPs and determined that the African LLSVP is ~1,000 km higher than the Pacific counterpart. Using geodynamic simulations, I find the height of a stable LLSVP is mainly controlled by its density and the ambient mantle viscosity. With ~1,000 numerical experiments, I conclude that the origin of the great height difference between the LLSVPs is that the African LLSVP is less dense, and thus less stable than the Pacific LLSVP. Next, I numerically identified another dynamic scenario accounting for the vastly different height of the two LLSVPs, which is caused by catastrophic sinking of accumulated subducted slabs at the 660-km boundary. Last, targeting one ancient carbonatite above the African LLSVP, I show that lithium isotopes in humite measured by nanoscale secondary ion mass spectrometry was able to uncover the signature of a subducted oceanic crust in its magma source, which may return from the interior to the surface by mantle plumes.

The redox conditions of Earth have been changing since proto-Earth’s accretion from the solar nebula. These changes have influenced the distribution and partitioning of volatile elements between the atmosphere and the mantle (Righter et al., 2020; Stagno and Fei, 2020. Though oxygen fugacity fO2 is arguably not the main factor for phase stability at certain pressure-temperature conditions (McCammon, 2005), it can influence which phases are stable, especially within a closed system such as the ones presented in this study. Despite the importance of controlling fO2 for interpreting the history of planetary bodies, there have been no methods to control the redox conditions in the laser-heated diamond anvil cell (LHDAC). This thesis has examined the feasibility for controlling redox conditions in the LHDAC using a mixture of Ar and H2 for insulation media. The experiments of this study were carried out at the GSECARS sector of the Advanced Photon Source at Argonne National Laboratory. In this study, α-Fe2O3 (hematite), ε-FeOOH (CaCl2-type), and Fe3O4 (magnetite) starting materials were used for probing changes of redox conditions. Experiments were also conducted with a pure Ar-medium for ε-FeOOH at the same pressure-temperature conditions of the hydrogen-bearing medium in order to provide a reference point for data which has uncontrolled redox conditions for an initially Fe(2+)-free material. The results for the ε-FeOOH starting material in Ar show transformation to ι-Fe2O3 (Rh2O3(II)-type) at 30.0 GPa and 1900 K, while in Ar + H2 it transformed to Fe5O7 with minor FeH (dhcp) at 30.0 GPa and 1850 K. For α-Fe2O3 in Ar + H2, it was found to convert to ε-FeOOH, Fe5O7, Fe5O6, and FeH (dhcp) at 36.5 GPa and 1800 K. For Fe3O4 in Ar + H2, it was found to convert to Fe4O5 (CaFe3O5-type), Fe5O6, and minor FeH (fcc) at 26.0 GPa and 1800 K. These results demonstrate that H in an Ar medium can promote the conversion of some Fe(3+) to Fe(2+) and Fe(0). However, the formation of ε-FeOOH in the α-Fe2O3 starting material suggests that H may participate in the chemical reaction of iron oxides.

Astrobiology is premised on the idea that life beyond Earth can exist. Yet, everything known about life is derivative from life on Earth. To understand life beyond Earth, then, requires a definition of life that is abstracted beyond a particular geophysical context. To do this requires a formal understanding of the physical mechanisms by which matter is animated into life. At current, such descriptions are completely lacking for the emergence of life, but do exist for the emergence of consciousness. Namely, contemporary neuroscience offers definitions for universal physical processes that are in one-to-one correspondence with conscious experience. Since consciousness is a sufficient condition for life, these universal definitions of consciousness offer an interesting way forward in terms of the search for life in the cosmos. In this work, I systematically examine Integrated Information Theory (IIT), a well-established theory of consciousness, with the aim of applying it in both biological and astrobiological settings. Surprisingly, I discover major problems with Integrated Information Theory on two fronts: mathematical and epistemological. On the mathematical side, I show how degeneracies buried deep within the theory render it mathematically ill-defined, while on the epistemological side, I prove that the postulates of IIT are scientifically unfalsifiable and inherently metaphysical. Given that IIT is the preeminent theory of consciousness in modern neuroscience, these results have far-reaching implications in this field. In addition, I show that the epistemic issues of falsifiability that hamstring IIT apply quite generally to all contemporary theories of consciousness, which suggests a major reframing of the problem is necessary. The problems that I reveal in regard to defining consciousness offer an important parallel in regard to defining life, as both fields seek to define their topic of study in absence of an existing theoretical framework. To avoid metaphysical problems related to falsifiability, universal theories of both life and consciousness must be framed with respect to independent empirical observations that can be used to benchmark predictions from the theory. In this regard, I argue that the epistemic debate over scientific theories of consciousness should be used to inform the discussion regarding theoretical definitions of life.

Volcanic eruptions can be serious geologic hazards, and have the potential to effect human life, infrastructure, and climate. Therefore, an understanding of the evolution and conditions of the magmas stored beneath volcanoes prior to their eruption is crucial for the ability to monitor such systems and develop effective hazard mitigation plans. This dissertation combines classic petrologic tools such as mineral chemistry and thermometry with novel techniques such as diffusion chronometry and statistical modeling in order to better understand the processes and timing associated with volcanic eruptions. By examining zoned crystals from the fallout ash of Yellowstone’s most recent supereruption, my work shows that the rejuvenation of magma has the ability to trigger a catastrophic supereruption at Yellowstone caldera in the years (decades at most) prior to eruption. This provides one of the first studies to thoroughly identify a specific eruption trigger of a past eruption using the crystal record. Additionally, through experimental investigation, I created a novel diffusion chronometer with application to determine magmatic timescales in silicic volcanic systems (i.e., rhyolite/dacite). My results show that Mg-in-sanidine diffusion operates simultaneously by both a fast and slow diffusion path suggesting that experimentally-derived diffusion chronometers may be more complex than previously thought. When applying Mg-in-sanidine chronometry to zoned sanidine from the same supereruption at Yellowstone, the timing between rejuvenation and eruption is further resolved to as short as five months, providing a greater understanding of the timing of supereruption triggers. Additionally, I developed a new statistical model to examine the controls on a single volcano’s distribution of eruptions through time, therefore the controls on the timing between successive eruptions, or repose time. When examining six Cascade volcanoes with variable distribution patterns through time, my model shows these distributions are not result of sampling bias, rather may represent geologic processes. There is a robust negative correlation between average repose time and average magma composition (i.e., SiO2), suggesting this may be a controlling factor of long-term repose time at Cascade volcanoes. Together, my work provides a better vision for forecasting models to mitigate potential destruction.

The mineralogy of the deep mantle is one of the key factors for the chemical evolution of the Earth. The constituent minerals of the mantle rock control the physical properties of the mantle, which have significant impacts on the large-scale processes occurring in the Earth's interior. In my PhD research, I adopted experimental approaches to investigate the mineralogy and the physical properties of the Earth's materials in the deep mantle by using the diamond anvil cells (DACs) combined with in-situ X-ray diffraction techniques.

First, I found that Ca-bearing bridgmanite can be stable in the deeper part of the Earth's lower mantle where temperature is sufficiently high. The dissolution of calcium into bridgmanite can change the physical properties of the mantle, such as compressibility and viscosity. This study suggests a new mineralogical model for the lower mantle, which is composed of the two layers depending on whether calcium dissolves in bridgmanite or forms CaSiO3 perovskite as a separate phase.

Second, I investigated the mineralogy and density of the subducting materials in the Archean at the P-T conditions near 670 km-depth. The experiments suggest that the major phases of Archean volcanic crust is majoritic garnet and ringwoodite in the P-T conditions of the deep transition zone, which become bridgmanite with increasing pressure. The density model showed that Archean volcanic crust would have been denser than the surrounding mantle, promoting sinking into the lower mantle regardless of the style of the transportation in the Archean.

Lastly, I further investigated the mineralogies and densities of the ancient volcanic crusts for the Archean and Proterozoic at the P-T conditions of the lower mantle. The experiments suggest that the mineralogy of the ancient volcanic crusts is composed mostly of bridgmanite, which is systemically denser than the surrounding lower mantle. This implies that the ancient volcanic crusts would have accumulated at the base of the mantle because of their large density and thickness. Therefore, the distinctive chemistry of the ancient volcanic crusts from the surrounding mantle would have given a rise to the chemical heterogeneities in the region for billions of years.

The interior of Earth is stratified due to gravity. Therefore, the lateral heterogeneities observed as seismic anomalies by seismologists are extremely interesting: they hold the key to understand the composition, thermal status and evolution of the Earth. This work investigates seismic anomalies inside Earth’s lowermost mantle and focuses on patch-like ultra-low velocity zones (ULVZs) found on Earth’s core-mantle boundary (CMB). Firstly, all previous ULVZ studies are compiled and ULVZ locations on the CMB are digitized. The result is a database, which is publicly available online. A key finding is that there is not a simple mapping between the locations of the observed ULVZs and the large low velocities provinces (LLVPs). Instead, ULVZs are more likely to occur near LLVP boundaries. This spatial correlation study supports a compositionally distinct origin for at least some ULVZs. Next, the seismic structure of the basal mantle beneath the Central America is investigated. This region hosts present and past subducted slabs, which could have brought compositionally distinct oceanic basalt all the way down to the CMB. The waveform distortions of a core-reflected seismic phase and a forward modeling method are used to constrain the causes of the CMB structures. In addition to ULVZ structures, isolated patches of thin zones with shear velocity increased by over 10% relative to background mantle are found for the first time. Ultra-high velocity zones (UHVZs) are interspersed with ULVZs and could be caused by subducted mid-ocean ridge basalt (MORB) that undergoes partial melting and melt segregation. Fe-rich partial melt of MORB can form ULVZs, and silica polymorphs (SiO2) and calcium-perovskite (CaPv) rich solid residue can explain the UHVZs. Finally, large-scale heterogeneities in the lowermost mantle are investigated using S waveform broadening observations. Several basal layer models are case-studied via synthetic calculations. S wave arrivals received at a distance larger than 80˚ in a global dataset from large earthquakes between the years 1994 and 2017 are examined and S waveform broadenings are documented. This approach exploits large distance data for the first time, and therefore is complementary to previous studies in terms of sampling locations. One possible explanation of S waveform broadening is velocity discontinuity inside the D″ layer due to the temperature controlled Bm-pPv phase transition.

The transport of hydrogen to the Earth’s deep interior remains uncertain. The upper mantle minerals have very low hydrogen solubilities (hundreds of ppm). The hydrogen storage capability in the transition zone minerals (2 wt%) is high compared to those of the upper mantle. The hydrogen storage in the lower mantle is not well known. The main minerals in the lower mantle bridgmanite and ferropericlase have very low hydrogen storage capacities (less than 20 ppm). In order to further understand the hydrogen storage in the lower mantle, a series of experiments had been conducted to simulate the environment similar to the Earth’s mantle. The experiments with hydrous Mg2SiO4 ringwoodite (Rw) show that it converts to crystalline dense hydrous silica, stishovite (Stv) or CaCl2-type SiO2(mStv), containing ∼1 wt% H2O together with bridgmanite (Brd) and MgO at the pressure-temperature conditions expected for lower mantle depths between approximately 660 to 1600 km. Brd would break down partially to dense hydrous silica (6–25 mol%) and(Mg,Fe)O in mid-mantle regions with 0.05–0.27 wt% H2O. The hydrous stishovite has a CaCl2 structure, which is common among hydrous minerals in the lower mantle. Based on this observation, I hypothesize the existence of hydrous phases in the lower mantle. The experiments found a new hexagonal iron hydroxide (η-Fe12O18+x/2Hx) between the stability fields of the epsilon and pyrite-type FeOOH at 60–80 GPa and high temperature. The new phase contains less H2O, limiting the H2O transport from the shallow to the deep mantle in the Fe–O–H system. Possible hydrogen storage in Ca-perovskite was studied. CaPv could contain 0.5–1 wt% water and the water in CaPv could distort the crystal structure of CaPv from cubic to tetragonal structure. In conclusion, hydrogen can be stored in hydrous stishovite in the shallower depth of the lower mantle. At greater depth, the new η phase and pyrite-type phase would take over the hydrogen storage. The role of CaPv in deep water storage needs to be considered in future studies.

This research investigates the fine scale structure in Earth's mantle, especially for the lowermost mantle, where strong heterogeneity exists. Recent seismic tomography models have resolved large-scale features in the lower mantle, such as the large low shear velocity provinces (LLSVPs). However, differences are present between different models, especially at shorter length scales. Fine scale structures both within and outside LLSVPs are still poorly constrained. The drastic growth of global seismic networks presents densely sampled seismic data in unprecedented quality and quantity. In this work, the Empirical Wavelet construction method has been developed to document seismic travel time and waveform information for a global shear wave seismic dataset. A dataset of 250K high-quality seismic records with comprehensive measurements is documented and made publicly available. To more accurately classify high quality seismic signal from the noise, 1.4 million manually labeled seismic records have been used to train a supervised classification model. The constructed model performed better than the empirical model deployed in the Empirical Wavelet method, with 87% in precision and 83% in recall. To utilize lower amplitude phases such as higher multiples of S and ScS waves, we have developed a geographic bin stacking method to improve signal-to-noise ratio. It is then applied to Sn waves up to n=6 and ScSn wave up to n=5 for both minor and major arc phases. The virtual stations constructed provide unique path sampling and coverage, vastly improving sampling in the Southern Hemisphere. With the high-quality dataset we have gathered, ray-based layer stripping iterative forward tomography is implemented to update a starting tomography model by mapping the travel time residuals along the ray from the surface down to the core mantle boundary. Final updated models with different starting tomography models show consistent updates, suggesting a convergent solution. The final updated models show higher resolution results than the starting tomography models, especially on intermediate-scale structures. The combined analyses and results in this work provide new tools and new datasets to image the fine-scale heterogeneous structures in the lower mantle, which advances our understanding of the dynamics and evolution of the Earth's mantle.

A mineral’s helium content reflects a balance between two competing processes: accumulation by radioactive decay and temperature-dependent diffusive loss. (U-Th)/He dating of zircon and other uranium and thorium-bearing minerals provides insight into the temperature histories of rocks at or near Earth’s surface that informs geoscientists’ understanding of tectonic and climate-driven exhumation, magmatic activity, and other thermal events. The crystal structure and chemistry of minerals affect helium diffusion kinetics, recorded closure temperatures, and interpretations of (U-Th)/He datasets. I used empirical and experimental methods to investigate helium systematics in two minerals chronometers: zircon and xenotime.

The same radioactivity that makes zircon a valuable chronometer damages its crystal structure over time and changes zircon helium kinetics. I used a zircon, titanite, and apatite (U-Th)/He dataset combined with previously published data and a new thermal model to place empirical constraints on the closure temperature for helium in a suite of variably damaged zircon crystals from the McClure Mountain syenite of Colorado. Results of this study suggest that the widely-used zircon damage accumulation and annealing model (ZRDAAM) does not accurately predict helium closure temperatures for a majority of the dated zircons. Detailed Raman maps of Proterozoic zircon crystals from the Lyon Mountain Granite of New York document complex radiation damage zoning. Models based on these results suggest that most ancient zircons are likely to exhibit intracrystalline variations in helium diffusivity due to radiation damage zoning, which may, in part, explain discrepancies between my empirical findings and ZRDAAM.

Zircon crystallography suggests that helium diffusion should be fastest along the crystallographic c-axis. I used laser depth profiling to show that diffusion is more strongly anisotropic than previously recognized. These findings imply that crystal morphology affects the closure temperature for helium in crystalline zircon. Diffusivity and the magnitude of diffusive anisotropy decrease with low doses of radiation damage.

Xenotime would make a promising (U-Th)/He thermochronometer if its helium kinetics were better known. I performed classic step-wise degassing experiments to characterize helium diffusion in xenotime FPX-1. Results suggest that this xenotime sample is sensitive to exceptionally low temperatures (∼50 °C) and produces consistent (U-Th)/He dates.