White Dwarf stars are the stellar remnants of low mass stars which have completed their evolution. Nearly all stars will become white dwarfs. The interior of a white dwarf encapsulates its evolution history: unraveling a white dwarf’s internal structure constrains the physical events which occurred to construct its composition. Variable, or pulsating, white dwarfs emit pulsations which are sensitive to their internal stratification. Just as seismology reveals Earth’s interior, asteroseismology can reveal stellar interiors. The standard approach to construe an observed white dwarf’s chemical makeup is to match observed pulsation properties to theoretical stellar models. Observed white dwarf pulsation data has reached 6-7 significant digits of precision. As such, it is important for computational modeling to consider systematic offsets from initial conditions and theoretical uncertainties that are within the detectable threshold. By analyzing the magnitude of pulsation differences among various uncertainties from white dwarf models, one can place constraints on important theoretical uncertainties. In this thesis, I explore impacts on white dwarf pulsations that result from accounting for various uncertainties in computational models. I start by showing the importance of 22Ne, and its impact on the pulsations in Helium atmosphere white dwarfs. Next, I discuss how certain trapped modes of white dwarfs may yield a signal for the 12C(α,γ)16O reaction rate probability distribution function. This reaction occurs during the Helium core burning phase in stellar evolution, and chiefly determines the Carbon and Oxygen abundance of white dwarfs. Following this work, I show how overshooting impacts the pulsation signatures of the 12C(α, γ)16O reaction rate. I then touch on the analytical work I’ve done regarding educational research in the HabWorlds course offered at Arizona State University (ASU). I then summarize my conclusions from these efforts.

Hydrogen is the main constituent of stars, and thus dominates the protoplanetary disc from which planets are born. Many planets may at some point in their growth have a high-pressure interface between refractory planetary materials and ahydrogen-dominated atmosphere. However, little experimental data for these materials at the relevant pressure-temperature conditions exists. I have experimentally explored the interactions between planetary materials and hydrogen at high P-T conditions utilizing the pulsed laser-heated diamond-anvil cell. First, I found that ferric/ferrous iron (as Fe2O3 hematite and (Mg,Fe)O ferropericlase) are reduced to metal by hydrogen: Fe2O3 + 4H2 → 2FeO + H2O + 3H2 → 2FeH + 3H2O and (Mg1−xFex) O + 3/2 xH2 → xFeH + (1 − x) MgO + xH2O respectively. This reduction of iron by hydrogen is important because it produces iron metal and water from iron oxide. This can partition H into the core (as FeH) or mantle (as H2O/OH−) of a growing planet. Next, I expanded my starting materials to silicates. I conducted experiments on San Carlos Olivine at pressures of 5-42 GPa. In the presence hydrogen, I observed the breakdown of molten magnesium silicate and the reduction of both iron and silicon to metal, forming alloys of both Fe-H and Fe-Si: Mg2SiO4 + 2H2 + 3Fe → 2MgO + FeSi + 2FeH + 2H2O. Similar experiments using natural fayalite (Fe2SiO4) as a starting material at pressures of 5-21 GPa yielded similar results. Hydrogen reduced iron to metal as it did in experiments with iron oxides. Unlike with San Carlos olivine, above 10 GPa silicon remained oxidized, implying the following reaction: Fe2SiO4 + 3H2 → 2FeH+2H2O +SiO2. However, below 7 GPa, silicon reduces and alloys with iron. The formation of Fe-Si alloys from silicates facilitated by hydrogen could have important effects for core composition in growing planets. I also observed at low pressures (<10 GPa), quenched iron melt can trap more hydrogen than previously thought (H/Fe nearly 2 instead of 1). This may have important effects for the chemical sequestration of a hydrogen atmosphere at shallow depths in an early magma ocean. All of the experimental work presented herein show that the composition, chemical partitioning, and phase stability of the condensed portion of growing planets can be modified via interaction with overlaying or ingassed volatile species.

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 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.

The pace of exoplanet discoveries has rapidly accelerated in the past few decades and the number of planets with measured mass and radius is expected to pick up in the coming years. Many more planets with a size similar to earth are expected to be found. Currently, software for characterizing rocky planet interiors is lacking. There is no doubt that a planet’s interior plays a key role in determining surface conditions including atmosphere composition and land area. Comparing data with diagrams of mass vs. radius for terrestrial planets provides only a first-order estimate of the internal structure and composition of planets [e.g. Seager et al 2007]. This thesis will present a new Python library, ExoPlex, which has routines to create a forward model of rocky exoplanets between 0.1 and 5 Earth masses. The ExoPlex code offers users the ability to model planets of arbitrary composition of Fe-Si-Mg-Al-Ca-O in addition to a water layer. This is achieved by modeling rocky planets after the earth and other known terrestrial planets. The three distinct layers which make up the Earth's internal structure are: core, mantle, and water. Terrestrial planet cores will be dominated by iron however, like earth, there may be some quantity of light element inclusion which can serve to enhance expected core volumes. In ExoPlex, these light element inclusions are S-Si-O and are included as iron-alloys. Mantles will have a more diverse mineralogy than planet cores. Unlike most other rocky planet models, ExoPlex remains unbiased in its treatment of the mantle in terms of composition. Si-Mg-Al-Ca oxide components are combined by predicting the mantle mineralogy using a Gibbs free energy minimization software package called Perple\_X [Connolly 2009]. By allowing an arbitrary composition, ExoPlex can uniquely model planets using their host star’s composition as an indicator of planet composition. This is a proven technique [Dorn et al 2015] which has not yet been widely utilized, possibly due to the lack of availability of easy to use software. I present a model sensitivity analysis to indicate the most important parameters to constrain in future observing missions. ExoPlex is currently available on PyPI so it may be installed using pip or conda on Mac OS or Linux based operating systems. It requires a specific scripting environment which is explained in the documentation currently stored on the ExoPlex GitHub page.

The dynamic Earth involves feedbacks between the solid crust and both natural and anthropogenic fluid flows. Fluid-rock interactions drive many Earth phenomena, including volcanic unrest, seismic activities, and hydrological responses. Mitigating the hazards associated with these activities requires fundamental understanding of the underlying physical processes. Therefore, geophysical monitoring in combination with modeling provides valuable tools, suitable for hazard mitigation and risk management efforts. Magmatic activities and induced seismicity linked to fluid injection are two natural and anthropogenic processes discussed in this dissertation.

Successful forecasting of the timing, style, and intensity of a volcanic eruption is made possible by improved understanding of the volcano life cycle as well as building quantitative models incorporating the processes that govern rock melting, melt ascending, magma storage, eruption initiation, and interaction between magma and surrounding host rocks at different spatial extent and time scale. One key part of such models is the shallow magma chamber, which is generally directly linked to volcano’s eruptive behaviors. However, its actual shape, size, and temporal evolution are often not entirely known. To address this issue, I use space-based geodetic data with high spatiotemporal resolution to measure surface deformation at Kilauea volcano. The obtained maps of InSAR (Interferometric Synthetic Aperture Radar) deformation time series are exploited with two novel modeling schemes to investigate Kilauea’s shallow magmatic system. Both models can explain the same observation, leading to a new compartment model of magma chamber. Such models significantly advance the understanding of the physical processes associated with Kilauea’s summit plumbing system with potential applications for volcanoes around the world.

The unprecedented increase in the number of earthquakes in the Central and Eastern United States since 2008 is attributed to massive deep subsurface injection of saltwater. The elevated chance of moderate-large damaging earthquakes stemming from increased seismicity rate causes broad societal concerns among industry, regulators, and the public. Thus, quantifying the time-dependent seismic hazard associated with the fluid injection is of great importance. To this end, I investigate the large-scale seismic, hydrogeologic, and injection data in northern Texas for period of 2007-2015 and in northern-central Oklahoma for period of 1995-2017. An effective induced earthquake forecasting model is developed, considering a complex relationship between injection operations and consequent seismicity. I find that the timing and magnitude of regional induced earthquakes are fully controlled by the process of fluid diffusion in a poroelastic medium and thus can be successfully forecasted. The obtained time-dependent seismic hazard model is spatiotemporally heterogeneous and decreasing injection rates does not immediately reduce the probability of an earthquake. The presented framework can be used for operational induced earthquake forecasting. Information about the associated fundamental processes, inducing conditions, and probabilistic seismic hazards has broad benefits to the society.

Energy harvesting from ambient is important to configuring Wireless Sensor Networks (WSN) for environmental data collecting. In this work, highly flexible thermoelectric generators (TEGs) have been studied and fabricated to supply power to the wireless sensor notes used for data collecting in hot spring environment. The fabricated flexible TEGs can be easily deployed on the uneven surface of heated rocks at the rim of hot springs. By employing the temperature gradient between the hot rock surface and the air, these TEGs can generate power to extend the battery lifetime of the sensor notes and therefore reduce multiple batteries changes where the environment is usually harsh in hot springs. Also, they show great promise for self-powered wireless sensor notes. Traditional thermoelectric material bismuth telluride (Bi2Te3) and advanced MEMS (Microelectromechanical systems) thin film techniques were used for the fabrication. Test results show that when a flexible TEG array with an area of 3.4cm2 was placed on the hot plate surface of 80°C in the air under room temperature, it had an open circuit voltage output of 17.6mV and a short circuit current output of 0.53mA. The generated power was approximately 7mW/m2.

On the other hand, high pressure, temperatures that can reach boiling, and the pH of different hot springs ranging from <2 to >9 make hot spring ecosystem a unique environment that is difficult to study. WSN allows many scientific studies in harsh environments that are not feasible with traditional instrumentation. However, wireless pH sensing for long time in situ data collection is still challenging for two reasons. First, the existing commercial-off-the-shelf pH meters are frequent calibration dependent; second, biofouling causes significant measurement error and drift. In this work, 2-dimentional graphene pH sensors were studied and calibration free graphene pH sensor prototypes were fabricated. Test result shows the resistance of the fabricated device changes linearly with the pH values (in the range of 3-11) in the surrounding liquid environment. Field tests show graphene layer greatly prevented the microbial fouling. Therefore, graphene pH sensors are promising candidates that can be effectively used for wireless pH sensing in exploration of hot spring ecosystems.