Filtering by
- Creators: Shirzaei, Manoochehr
- Creators: Barrett, The Honors College
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.
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.
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.
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.
Although nitrogen is the dominant element in Earth’s atmosphere, it is depleted in the bulk silicate Earth (relative to expected volatile abundances established by carbonaceous chondrites). To resolve this inconsistency, it has been hypothesized that this “missing nitrogen” may actually be stored within the Earth’s deep interior. In this work, we use multi-anvil press experiments to synthesize solid solution mixtures of the mantle transition zone mineral wadsleyite (Mg2SiO4) and silicon nitride (Si3N4). Successful synthesis of a 90% Si3N4, 10% Mg2SiO4 solid solution implies that nitrogen may not be sequestered within the most abundant mineral phases in the Earth’s mantle. Instead, nitrogen-rich accessory phases may hold the key to studying nitrogen storage within the deep interior. Ultimately, quantifying the amount of nitrogen within the mantle will further our understanding of the N cycle, which is vital to maintaining planetary habitability. Similar N cycling processes may be occurring on other rocky bodies; therefore, studying nitrogen storage may be an important part of determining habitability conditions on other worlds, both within in our solar system and beyond.