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

The movement between tectonic plates is accommodated through brittle (elastic) displacement on the plate boundary faults and ductile permanent deformation on the fault borderland. The elastic displacement along the fault can occur in the form of either large seismic events or aseismic slip, known as fault creep. Fault creep mainly occurs at the deep ductile portion of the crust, where the temperature is high. Nonetheless, aseismic creep can also occur on the shallow brittle portion of the fault segments that are characterized by frictionally weak material, elevated pore fluid pressure, or geometrical complexity. Creeping segments are assumed to safely release the accumulated strain(Kodaira et al., 2004; Rice, 1992)(Kodaira et al., 2004; Rice, 1992)(Kodaira et al., 2004; Rice, 1992)(Kodaira et al., 2004; Rice, 1992)(Kodaira et al., 2004; Rice, 1992) on the fault and also impede propagation of the seismic rupture. The rate of aseismic slip on creeping faults, however, might not be steady in time and instead consist of successive periods of acceleration and deceleration, known as slow slip events (SSEs). SSEs, which aseismically release the strain energy over a period of days to months, rather than the seconds to minutes characteristic of a typical earthquake, have been interpreted as earthquake precursors and as possible triggering factor for major earthquakes. Therefore, understanding the partitioning of seismic and aseismic fault slip and evolution of creep is fundamental to constraining the fault earthquake potential and improving operational seismic hazard models. Thanks to advances in tectonic geodesy, it is now possible to detect the fault movement in high spatiotemporal resolution and develop kinematic models of the creep evolution on the fault to determine the budget of seismic and aseismic slip.

In this dissertation, I measure the decades-long time evolution of fault-related crustal deformation along the San Andrea Fault in California and the northeast Japan subduction zone using space-borne geodetic techniques, such as Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar (InSAR). The surface observation of deformation combined with seismic data set allow constraining the time series of creep distribution on the fault surface at seismogenic depth. The obtained time-dependent kinematic models reveal that creep in both study areas evolves through a series of SSEs, each lasting for several months. Using physics-based models informed by laboratory experiments, I show that the transient elevation of pore fluid pressure is the driving mechanism of SSEs. I further investigate the link between SSEs and evolution of seismicity on neighboring locked segments, which has implications for seismic hazard models and also provides insights into the pattern of microstructure on the fault surface. I conclude that while creeping segments act as seismic rupture barriers, SSEs on these zones might promote seismicity on adjacent seismogenic segments, thus change the short-term earthquake forecast.

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.

The Jovian moon Europa's putative subsurface ocean offers one of the closest astrobiological targets for future exploration. It’s geologically young surface with a wide array of surface features aligned with distinct surface composition suggests past/present geophysical activity with implications for habitability. In this body of work, I propose a hypothesis for material transport from the ocean towards the surface via a convecting ice-shell. Geodynamical modeling is used to perform numerical experiments on a two-phase water-ice system to test the hypotheses. From these models, I conclude that it is possible for trace oceanic chemistry, entrapped into the newly forming ice at the ice-ocean phase interface, to reach near-surface. This new ice is advected across the ice-shell and towards the surface affirming a dynamical possibility for material transport across the ice-ocean system, of significance to astrobiological prospecting. Next, I use these self-consistent ice-ocean models to study the thickening of ice-shell over time. Europa is subject to the immense gravity field of Jupiter that generates tidal heating within the moon. Analysis of cases with uniform and localized internal tidal heating reveal that as the ice-shell grows from a warm initial ocean, there is an increase in the size of convection cells which causes a dramatic increase in the growth rate of the ice-shell. Addition of sufficient amount of heat also results in an ice-shell at an equilibrium thickness. Localization of tidal heating as a function of viscosity controls the equilibrium thickness. These models are then used to understand how compositional heterogeneity can be created in a growing ice-shell. Impurities (e.g. salts on the surface) that enter the ice-shell get trapped in the thickening ice-shell by freezing. I show the distribution pattern of heterogeneities that can form within the ice-shell at different times. This may be of potential application in identifying the longevity and mobility of brine pockets in Europa's ice-shell which are thought to be potential habitable niches.

For the geoscience community to continue to grow, students need to be attracted to the field. Here we examine the Incorporated Research Institutions for Seismology (IRIS) Research Experience for Undergraduates (REU) program to understand how the participants' experiences' affects their interest in geoscience and educational and career goals. Eleven interns over two years (2013-2014) were interviewed prior to the start of their internship, after their internship, and after presenting their research at the American Geophysical Union annual meeting. This internship program is of particular interest because many of the interns come into the REU with non-geoscience or geophysics backgrounds (e.g., physics, mathematics, chemistry, engineering). Both a priori and emergent codes are used to convert interview transcripts into quantitative data, which is analyzed alongside demographic information to understand how the REU influences their decisions. Increases in self-efficacy and exposure to multiple facets of geoscience research are expressed as primary factors that help shape their future educational and career goals. Other factors such as networking opportunities and connections during the REU also can play a role in their decision. Overall, REU participants who identified as geosciences majors solidified their decisions to pursue a career in geosciences, while participants who identified as non-geosciences majors were inclined to change majors, pursue geosciences in graduate school, or explore other job opportunities in the geosciences.

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.