Dynamic characterisation of four nine-story large-panel R.C. buildings in Bishkek (Kyrgyzstan): A comparison between experimental ambient vibration analysis and numerical finite element modeling

With the outlook of improving seismic vulnerability assessment for the city of Bishkek (Kyrgyzstan), the global dynamic behaviour of four nine-storey r.c. large-panel buildings in elastic regime is studied. The four buildings were built during the Soviet era within a serial production system. Since they all belong to the same series, they have very similar geometries both in plan and in height. Firstly, ambient vibration measurements are performed in the four buildings. The data analysis composed of discrete Fourier transform, modal analysis (frequency domain decomposition) and deconvolution interferometry, yields the modal characteristics and an estimate of the linear impulse response function for the structures of the four buildings. Then, finite element models are set up for all four buildings and the results of the numerical modal analysis are compared with the experimental ones. The numerical models are finally calibrated considering the first three global modes and their results match the experimental ones with an error of less then 20%.

From Craton to Rift: Empirically Based Ground-Truth Criteria for Local Events Recorded on Regional Networks

Region-specific empirically based ground-truth (EBGT) criteria used to estimate the epicentral-location accuracy of seismic events have been developed for the Main Ethiopian Rift and the Tibetan plateau. Explosions recorded during the Ethiopia-Afar Geoscientific Lithospheric Experiment (EAGLE), the International Deep Profiling of Tibet, and the Himalaya (INDEPTH III) experiment provided the necessary GT0 reference events. In each case, the local crustal structure is well known and handpicked arrival times were available, facilitating the establishment of the location accuracy criteria through the stochastic forward modeling of arrival times for epicentral locations. In the vicinity of the Main Ethiopian Rift, a seismic event is required to be recorded on at least 8 stations within the local Pg/Pn crossover distance and to yield a network-quality metric of less than 0.43 in order to be classified as EBGT5(95%) (GT5 with 95% confidence). These criteria were subsequently used to identify 10 new GT5 events with magnitudes greater than 2.1 recorded on the Ethiopian Broadband Seismic Experiment (EBSE) network and 24 events with magnitudes greater than 2.4 recorded on the EAGLE broadband network. The criteria for the Tibetan plateau are similar to the Ethiopia criteria, yet slightly less restrictive as the network-quality metric needs to be less than 0.45. Twenty-seven seismic events with magnitudes greater than 2.5 recorded on the INDEPTH III network were identified as GT5 based on the derived criteria. When considered in conjunction with criteria developed previously for the Kaapvaal craton in southern Africa, it is apparent that increasing restrictions on the network-quality metric mirror increases in the complexity of geologic structure from craton to plateau to rift. Accession Number: WOS:000322569200012

Geophysical Analysis of Seasonal Montandon Gravel Ridge Water Table Fluctuation and Moisture Gradient Variation due to Storm Events

The purpose of this research project is to continue exploring the Montandon Long-Term Hydrologic Research Site(LTHR) by using multiple geophysical methods to obtain more accurate and precise information regarding subsurface hydrologic properties of a local gravel ridge,which are important to both the health of surrounding ecosystems and local agriculture. Through using non-invasive geophysical methods such as seismic refraction, Direct Current resistivity and ground penetrating radar (GPR) instead of invasive methods such as boreholedrilling which displace sediment and may alter water flow, data collection is less likely to bias the data itself. In addition to imaging the gravel ridge subsurface, another important researchpurpose is to observe how both water table elevation and the moisture gradient (moisture content of the unsaturated zone) change over a seasonal time period and directly after storm events. The combination of three types of data collection allows the strengths of each method combine together and provide a relatively strongly supported conclusions compared to previous research. Precipitation and geophysical data suggest that an overall increase in precipitation during the summer months causes a sharp decrease in subsurface resistivity within the unsaturated zone. GPR velocity data indicate significant immediate increase in moisture content within the shallow vadose zone (< 1m), suggesting that rain water was infiltrating into the shallow subsurface. Furthermore, the combination of resistivity and GPR results suggest that the decreased resistivity within the shallow layers is due to increased ion content within groundwater. This is unexpected as rainwater is assumed to have a DC resistivity value of 3.33*105 ohm-m. These results may suggest that ions within the sediment must beincorporated into the infiltrating water.

Monitoring Stormwater Redistribution Into Low Resistivity Soils Using Noninvasive Geophysical Techniques

Water held in the unsaturated zone is important for agriculture and construction and is replenished by infiltrating rainwater. Monitoring the soil water content of clay soils using ground-penetrating radar (GPR) has not been researched, as clay soils cause attenuation of GPR signal. In this study, GPR common-midpoint soundings (CMPs) are used in the clayey soils of the Miller Run floodplain to monitor changes in the soil water content (SWC) before and after rainfall events. GPR accomplishes this task because increases in water content will increase the dielectric constant of the subsurface material, and decrease the velocity of the GPR wave. Using an empirical relationship between dielectric constant and SWC, the Topp relation, we are able to calculate a SWC from these velocity measurements. Non-invasive electromagnetics, resistivity, and seismic were performed, and from these surveys, the layering at the field site was delineated. EM characterized the horizontal variation of the soil, allowing us to target the most clay rich area. At the CMP location, resistivity indicates the vertical structure of the subsurface consists of a 40 cm thick layer with a resistivity of 100 ohm*m. Between 40 cm and 1.5 m is a layer with a resistivity of 40 ohm*m. The thickness estimates were confirmed with invasive auger and trenching methods away from the CMP location. GPR CMPs were collected relative to a July 2013 and September 2013 storm. The velocity observations from the CMPs had a precision of +/- 0.001 m/ns as assessed by repeat analysis. In the case of both storms, the GPR data showed the expected relationship between the rainstorms and calculated SWC, with the SWC increasing sharply after the rainstorm and decreasing as time passed. We compared these data to auger core samples collected at the same time as the CMPs were taken, and the volumetric analysis of the cores confirmed the trend seen in the GPR, with SWC values between 3 and 5 percent lower than the GPR estimates. Our data shows that we can, with good precision, monitor changes in the SWC of conductive soils in response to rainfall events, despite the attenuation induced by the clay.

Fastener-Based Computational Models of Cold-Formed Steel Shear Walls

Cold-formed steel (CFS) combined with wood sheathing, such as oriented strand board (OSB), forms shear walls that can provide lateral resistance to seismic forces. The ability to accurately predict building deformations in damaged states under seismic excitations is a must for modern performance-based seismic design. However, few static or dynamic tests have been conducted on the non-linear behavior of CFS shear walls. Thus, the purpose of this research work is to provide and demonstrate a fastener-based computational model of CFS wall models that incorporates essential nonlinearities that may eventually lead to improvement of the current seismic design requirements. The approach is based on the understanding that complex interaction of the fasteners with the sheathing is an important factor in the non-linear behavior of the shear wall. The computational model consists of beam-column elements for the CFS framing and a rigid diaphragm for the sheathing. The framing and sheathing are connected with non-linear zero-length fastener elements to capture the OSB sheathing damage surrounding the fastener area. Employing computational programs such as OpenSees and MATLAB, 4 ft. x 9 ft., 8 ft. x 9 ft. and 12 ft. x 9 ft. shear wall models are created, and monotonic lateral forces are applied to the computer models. The output data are then compared and analyzed with the available results of physical testing. The results indicate that the OpenSees model can accurately capture the initial stiffness, strength and non-linear behavior of the shear walls.

Characterization of harmonic tremor at Santiaguito volcano and its implications for eruption mechanisms

We observed Santiaguito volcano in southwestern Guatemala from March 2008 - March 2010. Seismic and infrasound data collected between January and March of 2009 contain records of many diverse processes occurring at the dacitic dome complex, including the recurrence of short lived (30-200 seconds in duration) harmonic tremor concurrent with ash poor gas emissions from the volcano. We employ several different analytical techniques to examine different portions of the tremor and source mechanisms. We use the parameters derived by this analysis to compare the feasibility of several suggested models of eruption mechanisms, and determine that this type of harmonic tremor is most justifiably generated by the flow of gas through crack networks generated by shear fracture along the magma conduit margin.

Dynamics and kinematics of eruptive activity at Fuego volcano, Guatemala 2005-2009

Volcanoes are the surficial expressions of complex pathways that vent magma and gasses generated deep in the Earth. Geophysical data record at least the partial history of magma and gas movement in the conduit and venting to the atmosphere. This work focuses on developing a more comprehensive understanding of explosive degassing at Fuego volcano, Guatemala through observations and analysis of geophysical data collected in 2005 – 2009. A pattern of eruptive activity was observed during 2005 – 2007 and quantified with seismic and infrasound, satellite thermal and gas measurements, and lava flow lengths. Eruptive styles are related to variable magma flux and accumulation of gas. Explosive degassing was recorded on broadband seismic and infrasound sensors in 2008 and 2009. Explosion energy partitioning between the ground and the atmosphere shows an increase in acoustic energy from 2008 to 2009, indicating a shift toward increased gas pressure in the conduit. Very-long-period (VLP) seismic signals are associated with the strongest explosions recorded in 2009 and waveform modeling in the 10 – 30 s band produces a best-fit source location 300 m west and 300 m below the summit crater. The calculated moment tensor indicates a volumetric source, which is modeled as a dike feeding a SW-dipping (35°) sill. The sill is the dominant component and its projection to the surface nearly intersects the summit crater. The deformation history of the sill is interpreted as: 1) an initial inflation due to pressurization, followed by 2) a rapid deflation as overpressure is explosively release, and finally 3) a reinflation as fresh magma flows into the sill and degasses. Tilt signals are derived from the horizontal components of the seismometer and show repetitive inflation deflation cycles with a 20 minute period coincident with strong explosions. These cycles represent the pressurization of the shallow conduit and explosive venting of overpressure that develops beneath a partially crystallized plug of magma. The energy released during the strong explosions has allowed for imaging of Fuego’s shallow conduit, which appears to have migrated west of the summit crater. In summary, Fuego is becoming more gas charged and its summit centered vent is shifting to the west - serious hazard consequences are likely.

Risk analysis of light-frame wood construction due to multiple hazards

Light-frame wood buildings are widely built in the United States (U.S.). Natural hazards cause huge losses to light-frame wood construction. This study proposes methodologies and a framework to evaluate the performance and risk of light-frame wood construction. Performance-based engineering (PBE) aims to ensure that a building achieves the desired performance objectives when subjected to hazard loads. In this study, the collapse risk of a typical one-story light-frame wood building is determined using the Incremental Dynamic Analysis method. The collapse risks of buildings at four sites in the Eastern, Western, and Central regions of U.S. are evaluated. Various sources of uncertainties are considered in the collapse risk assessment so that the influence of uncertainties on the collapse risk of lightframe wood construction is evaluated. The collapse risks of the same building subjected to maximum considered earthquakes at different seismic zones are found to be non-uniform. In certain areas in the U.S., the snow accumulation is significant and causes huge economic losses and threatens life safety. Limited study has been performed to investigate the snow hazard when combined with a seismic hazard. A Filtered Poisson Process (FPP) model is developed in this study, overcoming the shortcomings of the typically used Bernoulli model. The FPP model is validated by comparing the simulation results to weather records obtained from the National Climatic Data Center. The FPP model is applied in the proposed framework to assess the risk of a light-frame wood building subjected to combined snow and earthquake loads. The snow accumulation has a significant influence on the seismic losses of the building. The Bernoulli snow model underestimates the seismic loss of buildings in areas with snow accumulation. An object-oriented framework is proposed in this study to performrisk assessment for lightframe wood construction. For home owners and stake holders, risks in terms of economic losses is much easier to understand than engineering parameters (e.g., inter story drift). The proposed framework is used in two applications. One is to assess the loss of the building subjected to mainshock-aftershock sequences. Aftershock and downtime costs are found to be important factors in the assessment of seismic losses. The framework is also applied to a wood building in the state of Washington to assess the loss of the building subjected to combined earthquake and snow loads. The proposed framework is proven to be an appropriate tool for risk assessment of buildings subjected to multiple hazards. Limitations and future works are also identified.

Triggering of volcanic activity by large earthquakes

Statistical analyses of temporal relationships between large earthquakes and volcanic eruptions suggest seismic waves may trigger eruptions even over great (>1000 km) distances, although the causative mechanism is not well constrained. In this study the relationship between large earthquakes and subtle changes in volcanic activity was investigated in order to gain greater insight into the relationship between dynamic stresses propagated by surface waves and volcanic response. Daily measurements from the Ozone Monitoring Instrument (OMI), onboard the Aura satellite, provide constraints on volcanic sulfur-dioxide (SO2) emission rates as a measure of subtle changes in activity. Time series of SO2 emission rates were produced from OMI data for thirteen persistently active volcanoes from 1 October 2004 to 30 September 2010. In order to quantify the affect of earthquakes at teleseismic distances, we modeled surface-wave amplitudes from the source mechanisms of moment magnitude (Mw) ≥7 earthquakes, and calculated the Peak Dynamic Stress (PDS). We assessed the influence of earthquakes on volcanic activity in two ways: 1) by identifying increases in the SO2 time series data and looking for causative earthquakes and 2) by examining the average emission rate before and after each earthquake. In the first, the SO2 time series for each volcano was used to calculate a baseline threshold for comparison with post-earthquake emission. Next, we generated a catalog of responses based on sustained SO2 emission increases above this baseline. Delay times between each SO2 response and each prior earthquake were analyzed using both the actual earthquake catalog, and a randomly generated catalog of earthquakes. This process was repeated for each volcano. Despite varying multiple parameters, this analysis did not demonstrate a clear relationship between earthquake-generated PDS and SO2 emission. However, the second analysis, which was based on the occurrence of large earthquakes indicated a response at most volcanoes. Using the PDS calculations as a filtering criterion for the earthquake catalog, the SO2 mass for each volcano was analyzed in 28-day windows centered on the earthquake origin time. If the average SO2 mass after the earthquake was greater than an arbitrary percentage of pre-earthquake mass, we identified the volcano as having a response to the event. This window analysis provided insight on what type of volcanic activity is more susceptible to triggering by dynamic stress. The volcanoes with very open systems included in this study, Ambrym, Gaua, Villarrica, Erta Ale and, Turrialba, showed a clear response to dynamic stress while the volcanoes with more closed systems, Merapi, Semeru, Fuego, Pacaya, and Bagana, showed no response.

Ultraviolet digital imaging of volcanic plumes : implementation and application to magmatic processes at basaltic volcanoes

Magmatic volatiles play a crucial role in volcanism, from magma production at depth to generation of seismic phenomena to control of eruption style. Accordingly, many models of volcano dynamics rely heavily on behavior of such volatiles. Yet measurements of emission rates of volcanic gases have historically been limited, which has restricted model verification to processes on the order of days or longer. UV cameras are a recent advancement in the field of remote sensing of volcanic SO2 emissions. They offer enhanced temporal and spatial resolution over previous measurement techniques, but need development before they can be widely adopted and achieve the promise of integration with other geophysical datasets. Large datasets require a means by which to quickly and efficiently use imagery to calculate emission rates. We present a suite of programs designed to semi-automatically determine emission rates of SO2 from series of UV images. Extraction of high temporal resolution SO2 emission rates via this software facilitates comparison of gas data to geophysical data for the purposes of evaluating models of volcanic activity and has already proven useful at several volcanoes. Integrated UV camera and seismic measurements recorded in January 2009 at Fuego volcano, Guatemala, provide new insight into the system’s shallow conduit processes. High temporal resolution SO2 data reveal patterns of SO2 emission rate relative to explosions and seismic tremor that indicate tremor and degassing share a common source process. Progressive decreases in emission rate appear to represent inhibition of gas loss from magma as a result of rheological stiffening in the upper conduit. Measurements of emission rate from two closely-spaced vents, made possible by the high spatial resolution of the camera, help constrain this model. UV camera measurements at Kilauea volcano, Hawaii, in May of 2010 captured two occurrences of lava filling and draining within the summit vent. Accompanying high lava stands were diminished SO2 emission rates, decreased seismic and infrasonic tremor, minor deflation, and slowed lava lake surface velocity. Incorporation of UV camera data into the multi-parameter dataset gives credence to the likelihood of shallow gas accumulation as the cause of such events.

Understanding volcanic processes using UV camera measurements of sulfur dioxide and coincident infrasound and seismicity

The exsolution of volatiles from magma maintains an important control on volcanic eruption styles. The nucleation, growth, and connectivity of bubbles during magma ascent provide the driving force behind eruptions, and the rate, volume, and ease of gas exsolution can affect eruptive activity. Volcanic plumes are the observable consequence of this magmatic degassing, and remote sensing techniques allow us to quantify changes in gas exsolution. However, until recently the methods used to measure volcanic plumes did not have the capability of detecting rapid changes in degassing on the scale of standard geophysical observations. The advent of the UV camera now makes high sample rate gas measurements possible. This type of dataset can then be compared to other volcanic observations to provide an in depth picture of degassing mechanisms in the shallow conduit. The goals of this research are to develop a robust methodology for UV camera field measurements of volcanic plumes, and utilize this data in conjunction with seismoacoustic records to illuminate degassing processes. Field and laboratory experiments were conducted to determine the effects of imaging conditions, vignetting, exposure time, calibration technique, and filter usage on the UV camera sulfur dioxide measurements. Using the best practices determined from these studies, a field campaign was undertaken at Volcán de Pacaya, Guatemala. Coincident plume sulfur dioxide measurements, acoustic recordings, and seismic observations were collected and analyzed jointly. The results provide insight into the small explosive features, variations in degassing rate, and plumbing system of this complex volcanic system. This research provides useful information for determining volcanic hazard at Pacaya, and demonstrates the potential of the UV camera in multiparameter studies.

Evaluating stress patterns for Plio-Quaternary brittle deformation along the Calama-Olacapato-El Toro fault zone, Central Andes

Understanding the geometry and kinematics of the major structures of an orogen is important to elucidate its style of deformation, as well as its tectonic evolution. We describe the temporal and spatial changes in the state of stress of the trans-orogen area of the Calama-Olacapato-El Toro (COT) Fault Zone in the Central Andes, at about 24°S within the northern portion of the Puna Plateau between the Argentina-Chile border. The importance of the COT derives principally from the Quaternary-Holocene activity recognized on some segments, which may shed new light on its possible control on Quaternary volcanism and on the seismic hazard evaluation of the area. Field geological surveys along with kinematic analysis and numerical inversion of ∼140 new fault-slip measurements have revealed that this portion of the COT zone, previously considered a continuous, long-lived lineament, in reality has been subjected to three different kinematic regimes: 1) a Miocene transpressional phase with the maximum principal stress (σ1) chiefly trending NNE-SSW; 2) an extensional phase that started by 9 Ma, with a horizontal NW-SE-striking minimum principal stress (σ3) – permutations between σ2 and σ3 axes have been recognized at two sites – and 3) a left-lateral strike-slip phase with a horizontal ∼E-W &sigma1 and ∼N-S σ3 dating to the Late Pliocene-Quaternary. Spatially, in the Quaternary, the left-lateral component decreases toward the westernmost tip of the COT, where it transitions to extension; this produced to a N-S horst and graben structure. Hence, even if transcurrence is still active in the eastern portion of the COT, as focal mechanisms of crustal earthquakes indicate, our study demonstrates that extension is becoming the predominant structural style of deformation, at least in the western region. These major temporal and spatial changes in the tectonic regimes are attributed in part to changes in the magnitude of the boundary forces due to subduction processes. The overall orogen-perpendicular extension might be the result of vertical stress larger than both the horizontal stresses induced by gravitational effect of a thickened crust.

TECHNOLOGICAL ADVANCEMENT IN REMOTE SENSING AND INNOVATIVE APPLICATION TO VOLCANOLOGY

The integration of remote monitoring techniques at different scales is of crucial importance for monitoring of volcanoes and assessment of the associated hazard. In this optic, technological advancement and collaboration between research groups also play a key role. Vhub is a community cyberinfrastructure platform designed for collaboration in volcanology research. Within the Vhub framework, this dissertation focuses on two research themes, both representing novel applications of remotely sensed data in volcanology: advancement in the acquisition of topographic data via active techniques and application of passive multi-spectral satellite data to monitoring of vegetated volcanoes. Measuring surface deformation is a critical issue in analogue modelling of Earth science phenomena. I present a novel application of the Microsoft Kinect sensor to measurement of vertical and horizontal displacements in analogue models. Specifically, I quantified vertical displacement in a scaled analogue model of Nisyros volcano, Greece, simulating magmatic deflation and inflation and related surface deformation, and included the horizontal component to reconstruct 3D models of pit crater formation. The detection of active faults around volcanoes is of importance for seismic and volcanic hazard assessment, but not a simple task to be achieved using analogue models. I present new evidence of neotectonic deformation along a north-south trending fault from the Mt Shasta debris avalanche deposit (DAD), northern California. The fault was identified on an airborne LiDAR campaign of part of the region interested by the DAD and then confirmed in the field. High resolution LiDAR can be utilized also for geomorphological assessment of DADs, and I describe a size-distance analysis to document geomorphological aspects of hummock in the Shasta DAD. Relating the remote observations of volcanic passive degassing to conditions and impacts on the ground provides an increased understanding of volcanic degassing and how satellite-based monitoring can be used to inform hazard management strategies in nearreal time. Combining a variety of satellite-based spectral time series I aim to perform the first space-based assessment of the impacts of sulfur dioxide emissions from Turrialba volcano, Costa Rica, on vegetation in the surrounding environment, and establish whether vegetation indices could be used more broadly to detect volcanic unrest.

Spatial and Temporal Evolution of the Volcanics and Sediments of the Kenya Rift

New volumetric and mass flux estimates have been calculated for the Kenya Rift. Spatial and temporal histories for volcanic eruptions, lacustrine deposition, and hominin fossil sites are presented, aided by the compilation of a new digital geologic map. Distribution of volcanism over time indicates several periods of southward expansion followed by relative positional stasis. Volcanism occurs throughout the activated rift length, with no obvious abandonment as the rift system migrated. The main exception is a period of volcanic concentration around 10 Ma, when activity was constrained within 2° of the equator. Volumes derived from seismic data indicate a total volume of c. 310,000 km3 (2.47 x 1010 kg/yr ), which is significantly more than the map-derived volumes found here or published previously. Map-based estimates are likely affected by a bias against recognizing small volume events in the older record. Such events are, however, the main driver of erupted volume over the last 5 Ma. A technique developed here to counter this bias results in convergence of the two volume estimation techniques. Relative erupted composition over time is variable. Overall, the erupted material has a mafic to silicic ratio of 0.9:1. Basalts are distinctly more common in the Turkana region, which previously experienced Mesozoic rifting. Despite the near equal ratio of mafic to silicic products, the Kenya Rift otherwise fits the definition of a SLIP. It is proposed that the compositions would better fit the published definition if the Turkana region was not twice-rifted. Lacustrine sedimentation post-dates initial volcanism by about 5 million years, and follows the same volcanic trends, showing south and eastward migration over time. This sedimentation delay is likely related to timing of fault displacements. Evidence of hominin habitation is distinctly abundant in the northern and southern sections of the Kenya Rift, but there is an observed gap in the equatorial rift between 4 and 0.5 million years ago. After 0.5 Ma, sites appear to progress towards the equator. The pattern and timing of hominid site distributions suggests that the equatorial gap in habitation may be the result of active volcanic avoidance.

EVALUATION OF THE EVOLVING STRESS FIELD OF THE YELLOWSTONE VOLCANIC PLATEAU FROM 1988 TO 2010 FROM EARTHQUAKE FIRST MOTIONS INVERSION

Within the Yellowstone National Park, Wyoming, the silicic Yellowstone volcanic field is one of the most active volcanic systems all over the world. Although the last rhyolite eruption occurred around 70,000 years ago, Yellowstone is still believed to be volcanically active, due to high hydrothermal and seismic activity. The earthquake data used in this study cover the period of time between 1988 and 2010. Earthquake relocations and a set of 369 well-constrained, double-couple, focal mechanism solutions were computed. Events were grouped according to location and time to investigate trends in faulting. The majority of the events has oblique, normal-faulting solutions. The overall direction of extension throughout the 0.64 Ma Yellowstone caldera looks nearly ENE, consistently with the direction of alignments of volcanic vents within the caldera, but detailed study revealed spatial and temporal variations. Stress-field solutions for different areas and time periods were calculated from earthquake focal mechanism inversion. A well-resolved rotation of σ3 was found, from NNE-SSW near the Hebgen Lake fault zone, to ENE-WSW near Norris Junction. In particular, the σ3 direction changed throughout the years in the Norris Junction area, from being ENE-WSW, as calculated in the study by Waite and Smith (2004), to NNE-SSW, while the other σ3 directions are mostly unchanged over time. The Yellowstone caldera was subject to periods of net uplift and subsidence over the past century, explained in previous studies as caused by expanding or contracting sills, at different depths. Based on the models used to explain these deformation periods, we investigated the relationship between variability in aseismic deformation and seismic activity and faulting styles. Focal mechanisms and P and T axes were divided into temporal and depth intervals, in order to identify spatial or temporal trends in deformation. The presence of “chocolate tablet” structures, with composite dilational faults, was identified in many stages of the deformation history both in the Norris Geyser Basin area and inside the caldera. Strike-slip component movement was found in a depth interval below a contracting sill, indicating the movement of magma towards the caldera.