Geodetic strain rate and seismicity rates along the Apennines

Gli Appennini sono una regione tettonicamente attiva che è stata luogo di forti terremoti storici in passato, tra i quali il recente L’Aquila (2009) e le sequenze di Amatrice-Visso-Norcia (2016-2017). Una densità di stazioni GPS sempre maggiore permette di stimare il tasso di deformazione orizzontale con più risoluzione spaziale che in passato. In questa tesi studio la correlazione tra tassi di deformazione geodetica orizzontale e tassi di sismicità negli Appennini, esaminando possibili implicazioni per valutazioni del potenziale sismogenico. Uso un nuovo dataset di velocità GPS analizzando lo stato di deformazione attraverso il calcolo dello strain rate geodetico. Per questo calcolo utilizzo tre algoritmi, tra quelli presenti in letteratura scientifica, e discuto punti di forza e limitazioni di ciascuno. Seleziono poi una mappa di strain rate tra le altre e sviluppo un’analisi seguendo l’approccio di uno studio recente per investigare la relazione tra tassi di sismicità e tassi di deformazione lungo la catena Appenninica. Per farlo, verifico se le ipotesi che sostengono l’approccio sono verificate. Calcolo la sismicità di background a partire da un nuovo catalogo sismico per l’Italia e testo se gli eventi indipendenti seguono un processo di Poisson nel tempo. Infine, divido la regione di studio in aree a diversa deformazione sulla base dei tassi di deformazionee analizzado la distribuzione della sismicità di background in questi eventi e analizzando la relazione tra tassi di deformazione e tassi sismicità. Mostro che vi è una relazione è approssimativamente lineare ma con una dipendenza importante dalla risoluzione della mappa di strain rate. Descrivo poi alcune implicazioni che questa relazione avrebbe su stime di due quantità che concorrono a definire il potenziale sismogenico analizzando quali sono i parametri all’interno dell’approccio che possono influenzare maggiormente dette stime.

Present-day deformation of the Pyrenees revealed by GPS surveying and earthquake focal mechanisms until 2011

The Pyrenean mountain range is a slowly deforming belt with continuous and moderate seismic activity. To quantify its deformation field, we present the velocity field estimated from a GPS survey of the Pyrenees spanning 18 yr. The PotSis and ResPyr networks, including a total of 85 GPS sites, were installed and first measured in 1992 and 1995 1997, respectively, and remeasured in 2008 and 2010. We obtain a deformation field with velocities less than 1 mm yr−1 across the range. The estimated velocities for individual stations do not differ significantly from zero with 95 per cent confidence. Even so, we estimate a maximum extensional horizontal strain rate of 2.0 ± 1.7 nanostrain per year in a N S direction in the western part of the range. We do not interpret the vertical displacements due to their large uncertainties. In order to compare the horizontal strain rates with the seismic activity, we analyse a set of 194 focal mechanisms using three methods: (i) the 'r' factor relating their P and T axes, (ii) the stress tensors obtained by fault slip inversion and (iii) the strain-rate tensors. Stress and strain-rate tensors are estimated for: (i) the whole data set, (ii) the eastern and western parts of the range separately, and (iii) eight zones, which are defined based on the seismicity and the tectonic patterns of the Pyrenees. Each of these analyses reveals a lateral variation of the deformation style from compression and extension in the east to extension and strike-slip in the west of the range. Although the horizontal components of the strain-rate tensors estimated from the seismic data are slightly smaller in magnitude than those computed from the GPS velocity field, they are consistent within the 2σ uncertainties. Furthermore, the orientations of their principal axes agree with the mapped active faults.

An evaluation of seismic hazard in La Hispaniola, after the 2010 Haiti earthquake

An evaluation of the seismic hazard in La Hispaniola Island has been carried out, as part of the cooperative project SISMO-HAITI, supported by the Technical University of Madrid (UPM) and developed by several Spanish Universities, the National Observatory of Environment and Vulnerability) ONEV of Haiti, and with contributions from the Puerto Rico Seismic Network (PRSN) and University Seismological Institute of Dominican Republic (ISU). The study was aimed at obtaining results suitable for seismic design purposes. It started with the elaboration of a seismic catalogue for the Hispaniola Island, requiring an exhaustive revision of data reported by more than 20 seismic agencies, apart from these from the PRSN and ISU. The final catalogue contains 96 historical earthquakes and 1690 instrumental events, and it was homogenized to moment magnitude, Mw. Seismotectonic models proposed for the region were revised and a new regional zonation was proposed, taking into account geological andtectonic data, seismicity, focal mechanisms, and GPS observations. In parallel, attenuation models for subduction and crustal zones were revised in previous projects and the most suitable for the Caribbean plate were selected. Then, a seismic hazard analysis was developed in terms of peak ground acceleration, PGA, and spectral accelerations, SA (T), for periods of 0.1, 0.2, 0.5, 1 and 2s, using the Probabilistic Seismic Hazard Assessment (PSHA) methodology. As a result, different hazard maps were obtained for the quoted parameters, together with Uniform Hazard Spectra for Port au Prince and the main cities in the country. Hazard deaggregation was also carried out in these towns, for the target motion given by the PGA and SA (1s) obtained for return periods of 475, 975 and 2475 years. Therefore, the controlling earthquakes for short- and long-period target motions were derived. This study was started a few months after the 2010 earthquake, as a response to an aid request from the Haitian government to the UPM, and the results are available for the definition of the first building code in Haiti.

Intraplate earthquake swarm in Belo Jardim, NE Brazil: reactivation of a major Neoproterozoic shear zone (Pernambuco Lineament)

Intraplate earthquakes in stable continental areas have been explained basically by reactivation of pre-existing zones of weakness, stress concentration, or both. Zones of weakness are usually identified as sites of the last major orogeny, provinces of recent alkaline intrusions, or stretched crust in ancient rifts. However, it is difficult to identify specific zones of weakness and intraplate fault zones are not always easily correlated with known geological features. Although Northeastern Brazil is one of the most seismically active areas in the country (magnitudes 5 roughly every 5 yr), with hypocentral depths shallower than similar to 10 km and seismic zones as long as 30-40 km, no clear relationship with the known surface geology can be usually established with confidence, and a clear identification of zones of weakness has not yet been possible. Here we present the first clear case of seismic activity occurring as reactivation of an old structure in Brazil: the Pernambuco Lineament, a major Neoproterozoic shear zone. The 2004 earthquake swarm of Belo Jardim (magnitudes up to 3.1) and the recurrent activities in the nearby towns of Sao Caetano and Caruaru (magnitudes up to 4.0 and 3.8), show that the Pernambuco Lineament is a weak zone. A local seismic network showed that the Belo Jardim swarm of 2004 November occurred by normal faulting on a North dipping, E-W oriented fault plane in close agreement with the E-W trending structures within the Pernambuco Lineament. The Belo Jardim activity was concentrated in a 1.5 km (E-W) by 2 km (downdip) fault area, and average depth of 4.5 km. The nearby Caruaru activity occurs as both strike-slip and normal faulting, also consistent with local structures of the Pernambuco Lineament. The focal mechanisms of Belo Jardim, Caruaru and S. Caetano, indicate E-W compressional and N-S extensional principal stresses. The NS extension of this stress field is larger than that predicted by numerical models such as those of Coblentz & Richardson and we propose that additional factors such as flexural stresses from the nearby Sergipe-Alagoas marginal basin could also affect the current stress field in the Pernambuco Lineament.

The intraplate Porto dos Gauchos seismic zone in the Amazon craton - Brazil

The largest earthquake observed in the stable continental interior of the South American plate occurred in Serra do Tombador, Mato Grosso state - Brazil, on January 31,1955 with a magnitude of 6.2 m(b). Since then no other earthquake has been located near the 1955 epicentre. However, in Porto dos Gauchos, 100 km northeast of Serra do Tombador, a recurrent seismicity has been observed since 1959. Both Serra do Tombador and Porto dos Gauchos are located in the Phanerozoic Parecis basin. Two magnitude 5 earthquakes occurred in Porto dos Gauchos, in 1998 and 2005, with intensities up to VI and V, respectively. These two main shocks were followed by aftershock sequences lasting more than three years each. Local seismic stations have been deployed by the Seismological Observatory of the University of Brasilia since 1998 to study the ""Porto dos Gauchos"" seismic zone (PGSZ). A local seismic refraction survey was carried out with two explosions to help define the seismic velocity model. Both the 1998 and 2005 earthquake sequences occurred in the same WSW-ENE oriented fault zone with right-lateral strike-slip mechanisms. The epicentral zone is in the Parecis basin, near its northern border where there are buried grabens, generally trending WNW-ESE, such as the deep Mesoproterozoic Caiabis graben which lies partly beneath the Parecis basin. However, the epicentral distribution indicates that the 1998 and 2005 sequences are related to a N60 degrees E fault which probably crosses the entire Caiabis graben. The 1955 earthquake, despite the uncertainty in its epicentre, does not seem to be directly related to any buried graben either. The seismicity in the Porto dos Gauchos seismic zone, therefore, is not directly related to rifted crust. The probable direction of the maximum horizontal stress near Porto dos Gauchos is roughly E-W, consistent with other focal mechanisms further south in the Pantanal basin and Paraguay. but seems to be different from the NW-SE direction observed further north in the Amazon basin. The recurrent seismicity observed in Porto dos Gauchos, and the large 1955 earthquake nearby, make this area of the Parecis basin one of the most important seismic zones of Brazil. (C) 2009 Elsevier B.V. All rights reserved.

Sismicidade e correlação com feições geológicas: o caso do lineamento Pernambuco e seu entorno

This thesis presents and discusses the results of the various seismic areas in the State of Pernambuco, with the aim of having a vision of regional seismicity and its causes. To the papers published in journals were added two new original works submitted to international journals, dealing with seismic areas located in the counties of São Caetano, Cupira, and Agrestina. All seismic areas mentioned in this thesis are located on the Pernambuco Lineament and its surroundings (both in branches or single faults within 40 km of it). The Pernambuco Lineament is a Neoproterozoic shear zone of continental-scale that deformed the Borborema Province, and presents as branches, shear zones with NE-SW direction. The new submitted papers are from the analysis of data collected by three local networks of stations that operated in the following areas: network SO07 (seismicity in the district of Santa Luzia - São Caetano, 2007), network BM10 (data from seismic areas of Serra Verde ( Cupira) and Barra do Chata (Agrestina), in 2010), network SO10 (seismicity near the urban center of São Caetano in 2010). These data were used for determining the hypocenters and focal mechanisms in order to discuss the relationship between the seismicity and geological features of the area. The new mechanisms obtained, as well as the previously published allowed the determination of the direction of the average stress in the region. The direction of stress in the region involving the various seismic areas, now or previously studied, is quite stable and approximate EW direction (SHmax). The correlation between seismicity and geological features is observed on the lineament and north of it. In the south (Cupira and Agrestina), in seismic areas nearby shear zones NE-SW, there is no correlation and seismogenic EW normal faults are active and its motion is compatible with regional stresses. It is probable that these active faults are more recent than the Neoproterozoic, probably of the Cretaceous period, when the last great movement of the Pernambuco Lineament occurred

Estudo sísmico no açude Tucunduba, Senador Sá, CE

The Tucunduba Dam, is west of Fortaleza, Ceará State. The seismic monitoring of the area, with an analogical station and seven digital stations, had beginning on June 11, 1997. The digital stations, operated from June to November 1997. The data collected in the period of digital monitoring was analyzed for determination of hypocenters, focal mechanisms, and shear-wave anisotropy analysis. For determination of hypocenters, it was possible to find an active zone of nearly 1 km in length, with depth between 4.5 and 5.2 km. A 60AZ/88SE fault plane was determined using the least-squares method and hypocenters of a selected set of 16 earthquakes recorded. Focal mechanisms were determined, in the composite fault plane solution, a strike-slip fault, trending nearly E-W, was found. Single fault plane solutions were obteined to some earthquakes presented mean values of 65 (azimuth), and 80 (dip). Shear-wave anisotropy was found in the data. Polarization directions and travel time delays, between S spliting waves, were determined. It was not possible to obtain any conclusion on the cause of the observed anisotropy. It is not clear if there is correlation between seismicity and mapped faults in the area, although the directions obtained starting from the hipocentros and focal mechanism are they are consistent with directions, observed in the area, photo, topographic and fractures directions observed in the area

Estudoda sismicidade na região de Sobral - CE em 2008

The town of Sobral is located at the northwest part of the Ceará State, 250 km away from its capital, Fortaleza. In January 2008, an intense seismic activity began near Sobral with one event with magnitude 4,2mb on May 21. Since the start of its seismic activity, all events were recorded by the SBBR station (located on EMBRAPA Caprinos Farm), which operates in the region since August 2007. After this event, monitoring the seismic activity was carried out with the deployment of a local three component digital seismographic network, from June 5, 2008 until September 24. Initially, this network was composed of six seismographic stations. Later additional five stations were deployed until August 26 2008. This local network detected approximately 2,800 earthquakes. In this study we analyzed 581 earthquakes recorded by at least three stations for hypocentral and focal mechanism determination, and to contribute to a better explanation of the seismicity which in this region. To determine the hypocenters, we used a half-space model, with vP = 6,00 km/s and vP/vS = 1,71. From the hypocentral determination, it was revealed an active seismic zone with depth ranging between 1 and 8 km, 6 km long in E - W direction. The determination of fault planes and focal mechanism was obtained using the programs FPFIT and PLAN, which allowed comparison between their respective results in order to obtain more accurate results. A set of 24 earthquakes were selected to determine fault using PLAN planes and focal mechanisms using FPFIT. With the aid of detailed map of hypocenters this set, it was possible to identify three structures. Therefore, the set of 24 earthquakes were divided into three subsets. The type of mechanism was predominantly strike-slip with a dextral direction. Although the region has two tectonic structures near the site of the study area: the Café- Ipueiras Fault (normal fault) and the Sobral-Pedro II Lineament (dextral strike-slip fault) it was not possible to correlate the seismicity founded with those structures

Orientação e magnitude de tensões na bacia potiguar: implicações para evolução de bacias em margens passivas

This study presents new stress orientations and magnitudes from the Potiguar basin in the continental margin of Brazil. We analyzed breakout and drilled induced fractures derived from resistivity image logs run in ten oil wells. We also used direct Shmin measurements determined from hydraulic fractures and rock strength laboratory analysis. In addition, we compared these results with 19 earthquake focal mechanisms located in the crystalline basement. We observed that stress directions and magnitudes change across the basin and its basement. In the basin, the SHmax gradient of 20.0 MPa/km and the SHmax/Shmin ratio of 1.154 indicate a normal stress regime from 0.5 to 2.0 km, whereas the SHmax gradient of 24.5MPa/km and the SHmax/Shmin ratio of 1.396 indicate a strike slip stress regime from 2.5 to 4.0 km. The deeper strike-slip stress regime in the basin is similar to the regime in the basement at 1-12 km deep. This stress regime transition is consistent with an incipient tectonic inversion process in the basin. We also noted that the SHmax direction rotates from NW SE in the western part of the Potiguar basin to E W in its central and eastern part, following roughly the shoreline geometry. It indicates that local factors, as density contrast between continental and oceanic crust and sediment loading at the continental shelf influence the stress field. The concentration of fluid pressure in faults of the lowpermeability crystalline basement and its implications to establish a critically stressed fault regime in the basement is also discussed

Caracterização tectonoestratigráfica da porção adjacente ao epicentro do sismo São Vicente na bacia de Santos

Pós-graduação em Geologia Regional - IGCE

Tectonic geomorphology and active strain of the Northern Apennines mountain front

The Northern Apennines (NA) chain is the expression of the active plate margin between Europe and Adria. Given the low convergence rates and the moderate seismic activity, ambiguities still occur in defining a seismotectonic framework and many different scenarios have been proposed for the mountain front evolution. Differently from older models that indicate the mountain front as an active thrust at the surface, a recently proposed scenario describes the latter as the frontal limb of a long-wavelength fold (> 150 km) formed by a thrust fault tipped around 17 km at depth, and considered as the active subduction boundary. East of Bologna, this frontal limb is remarkably very straight and its surface is riddled with small, but pervasive high- angle normal faults. However, west of Bologna, some recesses are visible along strike of the mountain front: these perturbations seem due to the presence of shorter wavelength (15 to 25 km along strike) structures showing both NE and NW-vergence. The Pleistocene activity of these structures was already suggested, but not quantitative reconstructions are available in literature. This research investigates the tectonic geomorphology of the NA mountain front with the specific aim to quantify active deformations and infer possible deep causes of both short- and long-wavelength structures. This study documents the presence of a network of active extensional faults, in the foothills south and east of Bologna. For these structures, the strain rate has been measured to find a constant throw-to-length relationship and the slip rates have been compared with measured rates of erosion. Fluvial geomorphology and quantitative analysis of the topography document in detail the active tectonics of two growing domal structures (Castelvetro - Vignola foothills and the Ghiardo plateau) embedded in the mountain front west of Bologna. Here, tilting and river incision rates (interpreted as that long-term uplift rates) have been measured respectively at the mountain front and in the Enza and Panaro valleys, using a well defined stratigraphy of Pleistocene to Holocene river terraces and alluvial fan deposits as growth strata, and seismic reflection profiles relationships. The geometry and uplift rates of the anticlines constrain a simple trishear fault propagation folding model that inverts for blind thrust ramp depth, dip, and slip. Topographic swath profiles and the steepness index of river longitudinal profiles that traverse the anti- clines are consistent with stratigraphy, structures, aquifer geometry, and seismic reflection profiles. Available focal mechanisms of earthquakes with magnitude between Mw 4.1 to 5.4, obtained from a dataset of the instrumental seismicity for the last 30 years, evidence a clear vertical separation at around 15 km between shallow extensional and deeper compressional hypocenters along the mountain front and adjacent foothills. In summary, the studied anticlines appear to grow at rates slower than the growing rate of the longer- wavelength structure that defines the mountain front of the NA. The domal structures show evidences of NW-verging deformation and reactivations of older (late Neogene) thrusts. The reconstructed river incision rates together with rates coming from several other rivers along a 250 km wide stretch of the NA mountain front and recently available in the literature, all indicate a general increase from Middle to Late Pleistocene. This suggests focusing of deformation along a deep structure, as confirmed by the deep compressional seismicity. The maximum rate is however not constant along the mountain front, but varies from 0.2 mm/yr in the west to more than 2.2 mm/yr in the eastern sector, suggesting a similar (eastward-increasing) trend of the apenninic subduction.

Detailed study of the seismotectonic setting in the Lucanian Apennines and surrounding areas (Southern Italy)

In this research work I analyzed the instrumental seismicity of Southern Italy in the area including the Lucanian Apennines and Bradano foredeep, making use of the most recent seismological database available so far. I examined the seismicity occurred during the period between 2001 and 2006, considering 514 events with magnitudes M ≥ 2.0. In the first part of the work, P- and S-wave arrival times, recorded by the Italian National Seismic Network (RSNC) operated by the Istituto Nazionale di Geofisica e Vulcanologia (INGV), were re-picked along with those of the SAPTEX temporary array (2001–2004). For some events located in the Upper Val d'Agri, I also used data from the Eni-Agip oil company seismic network. I computed the VP/VS ratio obtaining a value of 1.83 and I carried out an analysis for the one-dimensional (1D) velocity model that approximates the seismic structure of the study area. After this preliminary analysis, making use of the records obtained in the SeSCAL experiment, I incremented the database by handpicking new arrival times. My final dataset consists of 15,666 P- and 9228 S-arrival times associated to 1047 earthquakes with magnitude ML ≥ 1.5. I computed 162 fault-plane solutions and composite focal mechanisms for closely located events. I investigated stress field orientation inverting focal mechanism belonging to the Lucanian Apennine and the Pollino Range, both areas characterized by more concentrated background seismicity. Moreover, I applied the double difference technique (DD) to improve the earthquake locations. Considering these results and different datasets available in the literature, I carried out a detailed analysis of single sub-areas and of a swarm (November 2008) recorded by SeSCAL array. The relocated seismicity appears more concentrated within the upper crust and it is mostly clustered along the Lucanian Apennine chain. In particular, two well-defined clusters were located in the Potentino and in the Abriola-Pietrapertosa sector (central Lucanian region). Their hypocentral depths are slightly deeper than those observed beneath the chain. I suggest that these two seismic features are representative of the transition from the inner portion of the chain with NE-SW extension to the external margin characterized by dextral strike-slip kinematics. In the easternmost part of the study area, below the Bradano foredeep and the Apulia foreland, the seismicity is generally deeper and more scattered and is associated to the Murge uplift and to the small structures present in the area. I also observed a small structure NE-SW oriented in the Abriola-Pietrapertosa area (activated with a swarm in November 2008) that could be considered to act as a barrier to the propagation of a potential rupture of an active NW-SE striking faults system. Focal mechanisms computed in this study are in large part normal and strike-slip solutions and their tensional axes (T-axes) have a generalized NE-SW orientation. Thanks to denser coverage of seismic stations and the detailed analysis, this study is a further contribution to the comprehension of the seismogenesis and state of stress of the Southern Apennines region, giving important contributions to seismotectonic zoning and seismic hazard assessment.

Fault delineation and stress orientations from the analysis of background, low magnitude seismicity in Southern Apennines (Italy)

The aim of this work was to show that refined analyses of background, low magnitude seismicity allow to delineate the main active faults and to accurately estimate the directions of the regional tectonic stress that characterize the Southern Apennines (Italy), a structurally complex area with high seismic potential. Thanks the presence in the area of an integrated dense and wide dynamic network, was possible to analyzed an high quality microearthquake data-set consisting of 1312 events that occurred from August 2005 to April 2011 by integrating the data recorded at 42 seismic stations of various networks. The refined seismicity location and focal mechanisms well delineate a system of NW-SE striking normal faults along the Apenninic chain and an approximately E-W oriented, strike-slip fault, transversely cutting the belt. The seismicity along the chain does not occur on a single fault but in a volume, delimited by the faults activated during the 1980 Irpinia M 6.9 earthquake, on sub-parallel predominant normal faults. Results show that the recent low magnitude earthquakes belongs to the background seismicity and they are likely generated along the major fault segments activated during the most recent earthquakes, suggesting that they are still active today thirty years after the mainshock occurrences. In this sense, this study gives a new perspective to the application of the high quality records of low magnitude background seismicity for the identification and characterization of active fault systems. The analysis of the stress tensor inversion provides two equivalent models to explain the microearthquake generation along both the NW-SE striking normal faults and the E- W oriented fault with a dominant dextral strike-slip motion, but having different geological interpretations. We suggest that the NW-SE-striking Africa-Eurasia convergence acts in the background of all these structures, playing a primary and unifying role in the seismotectonics of the whole region.

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.

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.