Matching Items (9)
Filtering by
- All Subjects: Geophysics
- Creators: Shirzaei, Manoochehr
- Creators: Desch, Steven
Description
Aquifers host the largest accessible freshwater resource in the world. However, groundwater reserves are declining in many places. Often coincident with drought, high extraction rates and inadequate replenishment result in groundwater overdraft and permanent land subsidence. Land subsidence is the cause of aquifer storage capacity reduction, altered topographic gradients which can exacerbate floods, and differential displacement that can lead to earth fissures and infrastructure damage. Improving understanding of the sources and mechanisms driving aquifer deformation is important for resource management planning and hazard mitigation.
Poroelastic theory describes the coupling of differential stress, strain, and pore pressure, which are modulated by material properties. To model these relationships, displacement time series are estimated via satellite interferometry and hydraulic head levels from observation wells provide an in-situ dataset. In combination, the deconstruction and isolation of selected time-frequency components allow for estimating aquifer parameters, including the elastic and inelastic storage coefficients, compaction time constants, and vertical hydraulic conductivity. Together these parameters describe the storage response of an aquifer system to changes in hydraulic head and surface elevation. Understanding aquifer parameters is useful for the ongoing management of groundwater resources.
Case studies in Phoenix and Tucson, Arizona, focus on land subsidence from groundwater withdrawal as well as distinct responses to artificial recharge efforts. In Christchurch, New Zealand, possible changes to aquifer properties due to earthquakes are investigated. In Houston, Texas, flood severity during Hurricane Harvey is linked to subsidence, which modifies base flood elevations and topographic gradients.
Poroelastic theory describes the coupling of differential stress, strain, and pore pressure, which are modulated by material properties. To model these relationships, displacement time series are estimated via satellite interferometry and hydraulic head levels from observation wells provide an in-situ dataset. In combination, the deconstruction and isolation of selected time-frequency components allow for estimating aquifer parameters, including the elastic and inelastic storage coefficients, compaction time constants, and vertical hydraulic conductivity. Together these parameters describe the storage response of an aquifer system to changes in hydraulic head and surface elevation. Understanding aquifer parameters is useful for the ongoing management of groundwater resources.
Case studies in Phoenix and Tucson, Arizona, focus on land subsidence from groundwater withdrawal as well as distinct responses to artificial recharge efforts. In Christchurch, New Zealand, possible changes to aquifer properties due to earthquakes are investigated. In Houston, Texas, flood severity during Hurricane Harvey is linked to subsidence, which modifies base flood elevations and topographic gradients.
ContributorsMiller, Megan Marie (Author) / Shirzaei, Manoochehr (Thesis advisor) / Reynolds, Stephen (Committee member) / Tyburczy, James (Committee member) / Semken, Steven (Committee member) / Werth, Susanna (Committee member) / Arizona State University (Publisher)
Created2018
Description
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.
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.
ContributorsZhai, Guang (Author) / Shirzaei, Manoochehr (Thesis advisor) / Garnero, Edward (Committee member) / Clarke, Amanda (Committee member) / Tyburczy, James (Committee member) / Li, Mingming (Committee member) / Arizona State University (Publisher)
Created2018
Description
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.
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.
ContributorsKhoshmanesh, Mostafa (Author) / Shirzaei, Manoochehr (Thesis advisor) / Arrowsmith, Ramon (Committee member) / Garnero, Edward (Committee member) / Tyburczy, James (Committee member) / Whipple, Kelin (Committee member) / Arizona State University (Publisher)
Created2018
Description
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.
ContributorsLorenzo, Alejandro M., Jr (Author) / Desch, Steven (Thesis advisor) / Shim, Dan S.-H. (Committee member) / Line, Michael (Committee member) / Li, Mingming (Committee member) / Arizona State University (Publisher)
Created2018
Description
The Kuiper Belt Object Haumea is one of the most fascinating objects in the solar system. Spectral reflectance observations reveal a surface of almost pure water ice, yet it has a mass of 4.006 × 1021 kg, measured from orbits of its moons, along with an inferred mean radius of 715 km, and these imply a mean density of around 2600 kg m−3. Thus the surface ice must be a veneer over a rocky core. This model is supported by observations of Haumea's light curve, which shows large photometric variations over an anomalously rapid 3.9154-hour rotational period. Haumea's surface composition is uniform, therefore the light curve must be due to a varying area presented to the observer, implying that Haumea has an oblong, ellipsoidal shape. If Haumea's rotation axis is normal to our line of sight, and Haumea reflects with a lunar-like scattering function, then its axis ratios are p = b/a = 0.80 (in the equatorial cross section) and q = c/a = 0.52 (in the polar cross section). In this work, I assume that Haumea is in hydrostatic equilibrium, and I model it as a two-phase ellipsoid with an ice mantle and a rocky core. I model the core assuming it has a given density in the range between 2700–3300 kg m−3 with axis ratios that are free to vary. The metric which my code uses calculates the angle between the gravity vector and the surface normal, then averages this over both the outer surface and the core-mantle boundary. When this fit angle is minimized, it allows an interpretation of the size and shape of the core, as well as the thickness of the ice mantle. Results of my calculations show that Haumea's most likely core density is 2700–2800 kg m−3, with ice thicknesses anywhere from 12–32 km over the poles and as thin as 4–18 km over the equator.
ContributorsProbst, Luke (Author) / Desch, Steven (Thesis advisor) / Asphaug, Erik (Committee member) / Bell, James (Committee member) / Arizona State University (Publisher)
Created2015
Description
Silicic volcanoes produce many styles of activity over a range of timescales. Eruptions vary from slow effusion of viscous lava over many years to violent explosions lasting several hours. Hazards from these eruptions can be far-reaching and persistent, and are compounded by the dense populations often surrounding active volcanoes. I apply and develop satellite and ground-based remote sensing techniques to document eruptions at Merapi and Sinabung Volcanoes in Indonesia. I use numerical models of volcanic activity in combination with my observational data to describe the processes driving different eruption styles, including lava dome growth and collapse, lava flow emplacement, and transitions between effusive and explosive activity.
Both effusive and explosive eruptions have occurred recently at Merapi volcano. I use satellite thermal images to identify variations during the 2006 effusive eruption and a numerical model of magma ascent to explain the mechanisms that controlled those variations. I show that a nearby tectonic earthquake may have triggered the peak phase of the eruption by increasing the overpressure and bubble content of the magma and that the frequency of pyroclastic flows is correlated with eruption rate. In 2010, Merapi erupted explosively but also shifted between rapid dome-building and explosive phases. I explain these variations by the heterogeneous addition of CO2 to the melt from bedrock under conditions favorable to transitions between effusive and explosive styles.
At Sinabung, I use photogrammetry and satellite images to describe the emplacement of a viscous lava flow. I calculate the flow volume (0.1 km3) and average effusion rate (4.4 m3 s-1) and identify active regions of collapse and advance. Advance rate was controlled by the effusion rate and the flow’s yield strength. Pyroclastic flow activity was initially correlated to the decreasing flow advance rate, but was later affected by the underlying topography as the flow inflated and collapsed near the vent, leading to renewed pyroclastic flow activity.
This work describes previously poorly understood mechanisms of silicic lava emplacement, including multiple causes of pyroclastic flows, and improves the understanding, monitoring capability, and hazard assessment of silicic volcanic eruptions.
Both effusive and explosive eruptions have occurred recently at Merapi volcano. I use satellite thermal images to identify variations during the 2006 effusive eruption and a numerical model of magma ascent to explain the mechanisms that controlled those variations. I show that a nearby tectonic earthquake may have triggered the peak phase of the eruption by increasing the overpressure and bubble content of the magma and that the frequency of pyroclastic flows is correlated with eruption rate. In 2010, Merapi erupted explosively but also shifted between rapid dome-building and explosive phases. I explain these variations by the heterogeneous addition of CO2 to the melt from bedrock under conditions favorable to transitions between effusive and explosive styles.
At Sinabung, I use photogrammetry and satellite images to describe the emplacement of a viscous lava flow. I calculate the flow volume (0.1 km3) and average effusion rate (4.4 m3 s-1) and identify active regions of collapse and advance. Advance rate was controlled by the effusion rate and the flow’s yield strength. Pyroclastic flow activity was initially correlated to the decreasing flow advance rate, but was later affected by the underlying topography as the flow inflated and collapsed near the vent, leading to renewed pyroclastic flow activity.
This work describes previously poorly understood mechanisms of silicic lava emplacement, including multiple causes of pyroclastic flows, and improves the understanding, monitoring capability, and hazard assessment of silicic volcanic eruptions.
ContributorsCarr, Brett B (Author) / Clarke, Amanda B (Thesis advisor) / Arrowsmith, Ramon (Committee member) / Mcnamara, Allen (Committee member) / Shirzaei, Manoochehr (Committee member) / Williams, Stanley (Committee member) / Arizona State University (Publisher)
Created2016
Description
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.
ContributorsAllu Peddinti, Divya (Author) / McNamara, Allen Keith (Thesis advisor) / Garnero, Edward (Committee member) / Desch, Steven (Committee member) / Zolotov, Mikhail (Committee member) / Clarke, Amanda (Committee member) / Arizona State University (Publisher)
Created2017
Description
Hydraulic fracturing, or fracking, has become a common practice in United States oil fields for enhancing their productivity. Among the concerns regarding fracking, however, is the possibility that it could trigger shallow earthquakes. The brine that results from fracking is injected into the subsurface for disposal. This brine causes a pore pressure gradient that is commonly believed to trigger failure along critically stressed subsurface faults. In Timpson, a small city in eastern Texas, earthquakes have become much more common since two injection wells were installed in 2007. 16 events of M_W > 2 have been detected since 2008 and are believed to be associated with failure along a subsurface fault. Applying interferometric synthetic aperture radar, we analyzed 3 sets of SAR images from the Advanced Land Observing Satellite (ALOS) from May 2007 to December 2010. Using these data sets, XX interferograms were generated. From these interferograms, it was possible to determine the spatial and temporal evolution of the crustal deformation in the line-of-sight of the satellite. The results show strong evidence of uplift in the region adjacent to the injection wells. While previous studies have established a strong connection between fluid injection and increased seismicity, this is to our knowledge the first observed case of crustal deformation that has been observed as a result of hydraulic fracturing fluid disposal.
ContributorsSedlak, Alexander Paul Woodward (Author) / Shirzaei, Manoochehr (Thesis director) / Vivoni, Enrique (Committee member) / Barrett, The Honors College (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Department of Chemistry and Biochemistry (Contributor) / School of Geographical Sciences and Urban Planning (Contributor) / School of Earth and Space Exploration (Contributor)
Created2014-05
Description
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.
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.
ContributorsKo, Byeongkwan (Author) / Shim, Sang-Heon (Thesis advisor) / Garnero, Edward (Committee member) / Leinenweber, Kurt (Committee member) / Li, Mingming (Committee member) / Desch, Steven (Committee member) / Arizona State University (Publisher)
Created2020