Modeling the Effect of Vapor Wall Deposition on the Formation of Secondary Organic Aerosol in Chamber Studies

Laboratory chamber experiments are used to investigate formation of secondary organic aerosol (SOA) from biogenic and anthropogenic precursors under a variety of environmental conditions. Simulations of these experiments test our understanding of the prevailing chemistry of SOA formation as well as the dynamic processes occurring in the chamber itself. One dynamic process occurring in the chamber that was only recently recognized is the deposition of vapor species to the Teflon walls of the chamber. Low-volatility products formed from the oxidation of volatile organic compounds (VOCs) deposit on the walls rather than forming SOA, decreasing the amount of SOA formed (quantified as the SOA yield: mass of SOA formed per mass of VOC reacted). In this work, several modeling studies are presented that address the effect of vapor wall deposition on SOA formation in chambers.

A coupled vapor-particle dynamics model is used to examine the competition among the rates of gas-phase oxidation to low volatility products, wall deposition of these products, and mass transfer to the particle phase. The relative time scales of these rates control the amount of SOA formed by affecting the influence of vapor wall deposition. Simulations show that an effect on SOA yield of changing the vapor-particle mass transfer rate is only observed when SOA formation is kinetically limited. For systems with kinetically limited SOA formation, increasing the rate of vapor-particle mass transfer by increasing the concentration of seed particles is an effective way to minimize the effect of vapor wall deposition.

This coupled vapor-particle dynamics model is then applied to α-pinene ozonolysis SOA experiments. Experiments show that the SOA yield is affected when changing the oxidation rate but not when changing the rate of gas-particle mass transfer by changing the concentration of seed particles. Model simulations show that the absence of an effect of changing the seed particle concentration is consistent with SOA formation being governed by quasi-equilibrium growth, in which gas-particle equilibrium is established much faster than the rate of change of the gas-phase concentration. The observed effect of oxidation rate on SOA yield arises due to the presence of vapor wall deposition: gas-phase oxidation products are produced more quickly and condense preferentially onto seed particles before being lost to the walls. Therefore, for α-pinene ozonolysis, increasing the oxidation rate is the most effective way to mitigate the influence of vapor wall deposition.

Finally, the detailed model GECKO-A (Generator for Explicit Chemistry and Kinetics of Organics in the Atmosphere) is used to simulate α-pinene photooxidation SOA experiments. Unexpectedly, α-pinene OH oxidation experiments show no effect when changing either the oxidation rate or the vapor-particle mass transfer rate, whereas GECKO-A predicts that changing the oxidation rate should drastically affect the SOA yield. Sensitivity studies show that the assumed magnitude of the vapor wall deposition rate can greatly affect conclusions drawn from comparisons between simulations and experiments. If vapor wall loss in the Caltech chamber is of order 10^-5 s^-1, GECKO-A greatly overpredicts SOA during high UV experiments, likely due to an overprediction of second-generation products. However, if instead vapor wall loss in the Caltech chamber is of order 10^-3 s^-1, GECKO-A greatly underpredicts SOA during low UV experiments, possibly due to missing autoxidation pathways in the α-pinene mechanism.

Probabilistic Imaging and Dynamic Modeling of Earthquake Source Processes

Investigation of large, destructive earthquakes is challenged by their infrequent occurrence and the remote nature of geophysical observations. This thesis sheds light on the source processes of large earthquakes from two perspectives: robust and quantitative observational constraints through Bayesian inference for earthquake source models, and physical insights on the interconnections of seismic and aseismic fault behavior from elastodynamic modeling of earthquake ruptures and aseismic processes.

To constrain the shallow deformation during megathrust events, we develop semi-analytical and numerical Bayesian approaches to explore the maximum resolution of the tsunami data, with a focus on incorporating the uncertainty in the forward modeling. These methodologies are then applied to invert for the coseismic seafloor displacement field in the 2011 Mw 9.0 Tohoku-Oki earthquake using near-field tsunami waveforms and for the coseismic fault slip models in the 2010 Mw 8.8 Maule earthquake with complementary tsunami and geodetic observations. From posterior estimates of model parameters and their uncertainties, we are able to quantitatively constrain the near-trench profiles of seafloor displacement and fault slip. Similar characteristic patterns emerge during both events, featuring the peak of uplift near the edge of the accretionary wedge with a decay toward the trench axis, with implications for fault failure and tsunamigenic mechanisms of megathrust earthquakes.

To understand the behavior of earthquakes at the base of the seismogenic zone on continental strike-slip faults, we simulate the interactions of dynamic earthquake rupture, aseismic slip, and heterogeneity in rate-and-state fault models coupled with shear heating. Our study explains the long-standing enigma of seismic quiescence on major fault segments known to have hosted large earthquakes by deeper penetration of large earthquakes below the seismogenic zone, where mature faults have well-localized creeping extensions. This conclusion is supported by the simulated relationship between seismicity and large earthquakes as well as by observations from recent large events. We also use the modeling to connect the geodetic observables of fault locking with the behavior of seismicity in numerical models, investigating how a combination of interseismic geodetic and seismological estimates could constrain the locked-creeping transition of faults and potentially their co- and post-seismic behavior.