Geology Reports⌕ Search

Geology topics

J. Ole Kaven

Publications and source records attributed to J. Ole Kaven.

16 recordsLinked to original sources

Long-term trends in microseismicity during operational shut-ins at the Coso Geothermal Field, California

Pausing injection and production can lead to induced seismicity in a variety of settings, with some of the largest events occurring during these so-called shut-ins. In geothermal fields, shut-ins are periodically conducted for maintenance on wells and surface infrastructure, thereby offering recurring means of estimating stress changes in the subsurface that lead to increased seismicity rates. Here, we examine daily production and microseismicity data from the Coso Geothermal Field (CGF) in California between 1996 and 2010 to investigate the repetitive trends of operational shut-in microseismicity. Using 24 local seismic stations, we first analyze spatial and temporal trends of over 60,000 earthquakes with magnitudes between -0.4 to 3.8. We find that the northern region exhibits no significant seismicity changes during shut-ins, whereas the rest of the field experiences induced seismicity during almost every shut-in with an increasing intensity towards the southern and eastern portions of the field, highlighting local differences in stress within the CGF. Additionally, we cluster the seismicity using waveform cross-correlation, revealing several earthquake clusters primarily occurring during shut-in periods. These observations suggest that certain fracture and fault sections respond quicker to changes in pore pressure and poroelastic stresses within the geothermal system, possibly highlighting main fluid pathways.

California↗

Induced seismicity strategic vision

Executive Summary The U.S. Geological Survey has a long history of contributions to the understanding and resolution of various scientific questions related to earthquakes associated with human activities, referred to as induced seismicity. Work started with the Rocky Mountain Arsenal studies in the 1960’s (Healy and others, 1968) when it was first discovered that fluid waste-disposal operations can cause earthquakes. U.S. Geological Survey work on induced seismicity continued in the intervening years but expanded dramatically in the early 2010s. Disposal of large volumes of wastewater, a byproduct of oil and gas production, into the subsurface caused an exponential increase in earthquakes in the central United States, including earthquakes that caused damage to buildings and infrastructure in nearby communities. Established in 2012 within the Earthquake Hazards Program (EHP), the Induced Seismicity Project (ISP) examines earthquakes that are induced by human activities to assess and mitigate the hazards associated with induced earthquakes as well as understand the conditions and processes controlling induced and natural earthquake generation and recurrence. The ISP examines instances of suspected induced earthquakes, in real-time and retrospectively, to assess the probabilistic hazard and investigate possibilities for reducing that hazard, something currently not possible with natural earthquake activity. The ISP has many synergies with work across multiple areas of EHP, including earthquake monitoring, hazard mitigation, and fundamental research into earthquake processes. The overarching questions that guided the development of our research strategies and approaches are as follows: How can seismicity that is induced by industrial activities be better distinguished from natural processes, particularly in regions with higher rates of natural seismicity? Can faults that may be at increased hazard for induced seismicity be identified based on existing fault structures, past earthquake activity, local stress conditions, or other information? What geologic and operational conditions control the occurrence, spatial and temporal evolution, maximum magnitude, ground shaking, and other characteristics of induced seismic sequences? What roles do surface and subsurface deformation, aseismic slip, and other processes play in the occurrence of injection-induced earthquakes? How can induced seismicity be better forecasted to support the Nation’s energy production and climate mitigation goals while limiting related earthquake hazards? In this report, background on induced seismicity and an overview of the state of knowledge is provided in three topical task areas: (1) oil and gas, (2) geothermal, and (3) geologic CO2 sequestration. This report then presents the EHP’s goals, strategies, and approaches—primarily focused on work in the ISP—for understanding the conditions that lead to induced seismicity, mitigating the associated hazards, and assisting the Nation in meeting its future needs in energy production and mitigation of climate change. Finally, the key limitations and opportunities for innovation are outlined at the end of the report.

Circular↗

Intelligent monitoring system for real-time geologic storage, optimization, and reservoir management

The objective of the subtask was to develop a near-real-time monitoring system for seismic data at the Decatur, IL, geologic carbon sequestration (GCS) site and specifically include fiber-optic cable derived distributed acoustic signal (DAS) data in the process. Owing to the large volumes of data, we opted to utilize existing deep borehole conventional seismic sensors for detection and pull DAS and shallow borehole seismic data once a detection has been made. Unfortunately, the horizontal fiber-optic cables did not yield microseismic signals for use in locating events near the GCS site. Various stacking and filtering approaches were tested without any coherent detection becoming apparent. We attribute the insensitivity to local microseismic events to a lack of coupling in downhole, vertical cable and the non-ideal alignment of the horizontal fiber-optic cable to the vertically polarized seismic energy. Despite the inability to detect and utilize the DAS data from the fiber-optic cables, we developed a general processing framework that enables easy adaptation for future deployment of fiber-optic cables with better suited alignment at Decatur or elsewhere.

Technical Report↗

Seismic hazard assessment at the Fallon, Nevada, Frontier Observatory for Research in Geothermal Energy site

Enhanced geothermal systems (EGS) technology aims to engineer a productive geothermal reservoir in regions of hot, but low permeability, rocks. In any EGS operation, the rock mass requires stimulation by high pressure injection of fluids, which has the potential to induce seismicity. To address the seismic hazard specifically, a probabilistic seismic hazard assessment (PSHA) is often required and is generally part of an induced seismicity mitigation plan (ISMP). A specific PSHA for the proposed Fallon, Nev., FORGE site is outlined below that relies solely on hypothetical stimulation scenarios and analog sites to assess the hazard of induced seismicity in the absence of local microseismicity. Partially due to the lack of existing seismicity at the site and partially to arrive at conservative estimates of the hazard, the PSHA is calculated for a range of b-values. Results indicate that the conservative estimates of seismic hazard at locations having significant, sensitive infrastructure near the proposed site are very low.

Nevada↗

Microseismic events associated with the Oroville Dam spillway

On 14 February 2017, two small (equivalent M D 0.8 and 1.0) seismic events occurred in proximity to the Oroville Dam in the Sierra Nevada foothills, California. To examine possible causal relationships between these events and reservoir operations, including the spillway failure starting prior to these events, we applied a new optimized template matching approach to seismic data between May 1993 - April 2018. We identified more than 19,000 smaller-magnitude events that were similar in character to the February 14 events. These events are located in proximity to the Oroville spillway and occurred in tight temporal clusters that strongly correlate with periods of spillway discharge. Seismic source inversion is inconclusive, but we suggest that these events might be induced by rapid changes in pore pressure along a fracture (or fractures) near the spillway. Cavitation cannot be ruled out, but it is unlikely to be the primary cause of the signals observed because these events are intermittent, impulsive and of short duration. The inferred repetitive opening and closing of the fracture(s) occurred long before any damage to the spillway and is thus probably not directly associated with spillway failure in February 2017. These events were not related to the 1975 M L 5.7 earthquake sequence that may have been induced by the filling of the Oroville reservoir.

California↗

Physical properties of sidewall cores from Decatur, Illinois

To better assess the reservoir conditions influencing the induced seismicity hazard near a carbon dioxide sequestration demonstration site in Decatur, Ill., core samples from three deep drill holes were tested to determine a suite of physical properties including bulk density, porosity, permeability, Young’s modulus, Poisson’s ratio, and failure strength. Representative samples of the shale cap rock, the sandstone reservoir, and the Precambrian basement were selected for comparison. Physical properties were strongly dependent on lithology. Bulk density was inversely related to porosity, with the cap rock and basement samples being both least porous (<3 percent) and densest (~2.6 grams per cubic centimeter [g/cc]). Permeability was highest in the reservoir sandstones (10-15 to 10-18 meters squared [m2]) relative to the cap rock and basement rocks (<10-21 m2). Young’s modulus was distinctly higher in the basement rocks (45 to 80 gigapascal [GPa]) compared to the cap rock and sandstones (19 to 57 GPa). Poisson’s ratio for the sandstones varied widely (0.14 to 0.27), but the highest values were similar to the cap rock and basement rocks (0.24 to 0.28). These physical properties reflect the layered structure of the reservoir and adjacent rocks at the Decatur site. However, within the sandstone there is a great deal of lithologic variety, accounting for the large range in physical parameters for this geologic unit.

Illinois↗

Micro-seismicity within the Coso Geothermal field, California, from 1996-2012

We extend our previous catalog of seismicity within the Coso Geothermal field by adding over two and a half years of additional data to prior results. In total, we locate over 16 years of seismicity spanning from April 1996 to May of 2012 using a refined velocity model, apply it to all events and utilize differential travel times in relocations to improve the accuracy of event locations. The improved locations elucidate major structural features within the reservoir that we interpret to be faults that contribute to heat and fluid flow within the reservoir. Much of the relocated seismicity remains diffuse between these major structural features, suggesting that a large volume of accessible and distributed fracture porosity is maintained within the geothermal reservoir through ongoing brittle failure. We further track changes in b value and seismic moment release within the reservoir as a whole through time. We find that b values decrease significantly during 2009 and 2010, coincident with the occurrence of a greater number of moderate magnitude earthquakes (3.0 ≤ ML < 4.5). Analysis of spatial variations in seismic moment release between years reveals that localized seismicity tends to spread from regions of high moment release into regions with previously low moment release, akin to aftershock sequences. These results indicate that the Coso reservoir is comprised of a network of fractures at a variety of spatial scales that evolves dynamically over time, with progressive changes in characteristics of microseismicity and inferred fractures and faults that are only evident from a long period of seismic monitoring analyzed using self-consistent methods.

California↗

Geologic setting of the proposed West Flank Forge Site, California: Suitability for EGS research and development

The proposed West Flank FORGE site is within the China Lake Naval Air Weapons Station (NAWS), China Lake, CA. The West Flank is west of the Coso geothermal field, an area of China Lake NAWS dominated by the Quaternary Coso volcanic field largely comprised of rhyolite domes and their volcaniclastic and epiclastic horizons. The largest dome flow complex, Sugarloaf Mountain, marks the northwestern margin of the geothermal field. The West Flank is situated due west of Sugarloaf. The geologic setting of the West Flank was determined from one deep well (83-11) drilled as a potential production hole in 2009. The bottom-hole temperature (BHT) of well 83-11 approaches 600 oF (315˚C), but flow tests demonstrate very low, non-commercial permeabilities. With the exception of the upper 600 feet of volcaniclastic alluvium, well 83-11 is completed in granitic basement. The West Flank possesses the primary attributes of a FORGE site: non-commercial permeability ( <10 -16 m 2 ), a 175˚ to 225˚C temperature range in crystalline rocks, and a location outside an existing geothermal field The Coso Mountains host the Coso volcanic field and are within a right-releasing stepover between the dextral Airport Lake (ALF) and Little Lake fault zones (LLFZ) and the Wild Horse Mesa and Owens Valley faults. Two distinct fault populations have been identified at Coso: WNW-trending and antithetical, NE-trending strike-slip faults and N- to NNE-trending normal faults. These faults are both high permeability drilling targets at depth within the main (productive) geothermal field and they locally segment the field into distinct hydrothermal regimes. The West Flank may be segmented from the rest of the field by one such northerly trending fault. The overall minimum principal stress orientation in the main geothermal field varies from 103˚ to 108˚; however, the minimum horizontal principal stress in 83-11 is rotated to 081˚.

California↗

Compartmentalization of the Coso East Flank geothermal field imaged by 3-D full-tensor MT inversion

Previous magnetotelluric (MT) studies of the high-temperature Coso geothermal system in California identified a subvertical feature of low resistivity (2–5 Ohm m) and appreciable lateral extent (>1 km) in the producing zone of the East Flank field. However, these models could not reproduce gross 3-D effects in the recorded data. We perform 3-D full-tensor inversion and retrieve a resistivity model that out-performs previous 2-D and 3-D off-diagonal models in terms of its fit to the complete 3-D MT data set as well as the degree of modelling bias. Inclusion of secondary Z xx and Z yy data components leads to a robust east-dip (60†) to the previously identified conductive East Flank reservoir feature, which correlates strongly with recently mapped surface faults, downhole well temperatures, 3-D seismic reflection data, and local microseismicity. We perform synthetic forward modelling to test the best-fit dip of this conductor using the response at a nearby MT station. We interpret the dipping conductor as a fractured and fluidized compartment, which is structurally controlled by an unmapped blind East Flank fault zone.

Coso Range↗

Geologic setting of the West Flank, a FORGE site adjacent to the Coso geothermal field

The West Flank FORGE (Frontier Observatory for Research in Geothermal Energy) site is located immediately west and outside of the Coso geothermal field, eastern California. Coso is a fluid-dominated, high temperature geothermal system that has been producing power continuously since 1987. The reservoir is composed of highly faulted, fractured and hydrothermally altered Cretaceous and Jurassic plutonic basement rocks. The heat source is a shallow, silicic magma chamber associated with overlying, Quaternary rhyolites and basalts of the Coso volcanic field. Over 30 years of development drilling and associated investigations demonstrates that several well-defined boundaries exist at Coso beyond which there is no commercial, hydrothermal geothermal potential. The West Flank is just west of one of these boundaries and it meets all required FORGE temperature (175-225˚C), host rock (crystalline rock), and depth (1.5-4.0 km) criteria. It has been hypothesized that the Coso volcanic field exists within a right-releasing step-over between two NW-striking, dextral fault systems. Two distinct fault populations yield high permeability drilling targets in the geothermal field. The first population contains WNW-trending with antithetical, NE-trending strike-slip faults and the second includes N- to NNE-trending normal faults. The West Flank appears to be separated from the geothermal field by one such northerly trending fault bounded by a felsic dike swarm. The West Flank’s 83-11 well suggests that the West Flank is comprised of weakly altered, leucogranite and diorites typical of the Jurassic to Cretaceous basement. The stress field in the West Flank has been rotated to 81˚ in contrast to the minimum principal stress orientations within the geothermal field which range from 103˚ to 108˚. Limited drilling, geophysics and other data sets are being assessed and synthesized to develop a working, 3d conceptual geologic model of the West Flank for FORGE.

California↗

Surface monitoring of microseismicity at the Decatur, Illinois, CO2 sequestration demonstration site

Sequestration of CO2 into subsurface reservoirs can play an important role in limiting future emission of CO2 into the atmosphere (e.g., Benson and Cole, 2008). For geologic sequestration to become a viable option to reduce greenhouse gas emissions, large-volume injection of supercritical CO2 into deep sedimentary formations is required. These formations offer large pore volumes and good pore connectivity and are abundant (Bachu, 2003; U.S. Geological Survey Geologic Carbon Dioxide Storage Resources Assessment Team, 2013). However, hazards associated with injection of CO2 into deep formations require evaluation before widespread sequestration can be adopted safely (Zoback and Gorelick, 2012). One of these hazards is the potential to induce seismicity on pre-existing faults or fractures. If these faults or fractures are large and critically stressed, seismic events can occur with magnitudes large enough to pose a hazard to surface installations and, possibly more critical, the seal integrity of the cap rock. The Decatur, Illinois, carbon capture and storage (CCS) demonstration site is the first, and to date, only CCS project in the United States that injects a large volume of supercritical CO2 into a regionally extensive, undisturbed saline formation. The first phase of the Decatur CCS project was completed in November 2014 after injecting a million metric tons of supercritical CO2 over three years. This phase was led by the Illinois State Geological Survey (ISGS) and included seismic monitoring using deep borehole sensors, with a few sensors installed within the injection horizon. Although the deep borehole network provides a more comprehensive seismic catalog than is presented in this paper, these deep data are not publicly available. We contend that for monitoring induced microseismicity as a possible seismic hazard and to elucidate the general patterns of microseismicity, the U.S. Geological Survey (USGS) surface and shallow borehole network described below provides an adequate event detection threshold. The formation targeted for injection is the Mount Simon Sandstone, which is laterally extensive, has high porosity and permeability and has the potential to host future CCS projects due to its favorable hydrologic characteristics and proximity to industrial sources of CO2 (Birkholzer and Zhou, 2009). At Decatur, CO2, a byproduct of ethanol production at the Archer Daniels Midland (ADM) facility, is compressed to supercritical state and injected at 2.1 km depth into the 460 m thick Mount Simon Sandstone. This sandstone has varying properties, ranging from the lower, fine- to coarse-grained sandstone with high permeability and porosity, to the middle and upper Mount Simon, which consist of planar, cross-bedded layers of varied permeability and porosity (Leetaru and Freiburg, 2014). The changes in permeability and porosity within the Mount Simon Sandstone, due to depositional and diagenetic differences, create horizontal baffles, which inhibit vertical flow and restrict the injected CO2 to remain near the injection horizon (Bowen et al., 2011). The lowest portion of the Mount Simon Sandstone overlying the Precambrian rhyolite basement is the Pre-Mount Simon interval, generally < 15 m in thickness and composed of fine- to medium-grain size sandstone that is highly deformed (Leetaru and Freiburg, 2014). The basement rhyolite has a clayrich matrix and is fractured, with significant alterations within the fractures. The primary sealing cap rock is the Eau Claire Formation, a 100–150 m thick unit at a depth of roughly 1.69 km (Leetaru and Freiburg, 2014). The Maquoketa Shale Group and the New Albany Shale serve as secondary and tertiary seals at shallower depths of ∼820 and ∼650 m, respectively. The ISGS managed the Illinois Basin–Decatur Project (IBDP), a three-year project beginning in November 2011, during which carbon dioxide was injected at a rate of ∼1000 metric tons per day until November 2014 (Finley et al., 2011, 2013). ADM manages the Illinois Industrial CCS (ICCS) project, which will inject ∼3000 metric tons/day into a second injection well starting in the summer of 2015. The USGS began monitoring microseismicity with a 13- station seismic network at Decatur in July 2013 (Fig. 1). This network provides good detection capabilities and azimuthal (focal sphere) coverage for microseismicity with moment magnitudes (Mw) above about −0:5. Here, we report on 19 months of microseismicity monitoring at the Decatur CO2 sequestration site, which permits a detailed look at the evolution and character of injection-induced seismicity.

Illinois↗

Natural or Induced: Identifying Natural and Induced Swarms from Pre-production and Co-production Microseismic Catalogs at the Coso Geothermal Field

Increased levels of seismicity coinciding with injection of reservoir fluids have prompted interest in methods to distinguish induced from natural seismicity. Discrimination between induced and natural seismicity is especially difficult in areas that have high levels of natural seismicity, such as the geothermal fields at the Salton Sea and Coso, both in California. Both areas show swarm-like sequences that could be related to natural, deep fluid migration as part of the natural hydrothermal system. Therefore, swarms often have spatio-temporal patterns that resemble fluid-induced seismicity, and might possibly share other characteristics. The Coso Geothermal Field and its surroundings is one of the most seismically active areas in California with a large proportion of its activity occurring as seismic swarms. Here we analyze clustered seismicity in and surrounding the currently produced reservoir comparatively for pre-production and co-production periods. We perform a cluster analysis, based on the inter-event distance in a space-time-energy domain to identify notable earthquake sequences. For each event j, the closest previous event i is identified and their relationship categorized. If this nearest neighbor&rsquo;s distance is below a threshold based on the local minimum of the bimodal distribution of nearest neighbor distances, then the event j is included in the cluster as a child to this parent event i. If it is above the threshold, event j begins a new cluster. This process identifies subsets of events whose nearest neighbor distances and relative timing qualify as a cluster as well as a characterizing the parent-child relationships among events in the cluster. We apply this method to three different catalogs: (1) a two-year microseismic survey of the Coso geothermal area that was acquired before exploration drilling in the area began; (2) the HYS_catalog_2013 that contains 52,000 double-difference relocated events and covers the years 1981 to 2013; and (3) a catalog of 57,000 events with absolute locations from the local network recorded between 2002 and 2007. Using this method we identify 10 clusters of more than 20 events each in the pre-production survey and more than 200 distinct seismicity clusters that each contain at least 20 and up to more than 1000 earthquakes in the more extensive catalogs. The cluster identification method used yields a hierarchy of links between multiple generations of parent and offspring events. We analyze different topological parameters of this hierarchy to better characterize and thus differentiate natural swarms from induced clustered seismicity and also to identify aftershock sequences of notable mainshocks. We find that the branching characteristic given by the average number of child events per parent event is significantly different for clusters below than for clusters around the produced field.

Conference Paper↗

Seismic monitoring at the Decatur, Ill., CO2 sequestration demonstration site

The viability of carbon capture and storage (CCS) to reduce emissions of greenhouse gases depends on the ability to safely sequester large quantities of CO2 over geologic time scales. One concern with CCS is the potential of induced seismicity. We report on ongoing seismic monitoring by the U.S. Geological Survey (USGS) at a CCS demonstration site in Decatur, IL, in an effort to understand the potential hazards posed by injection-induced seismicity associated with geologic CO2 sequestration. At Decatur, super-critical CO2 is injected at 2.1 km depth into the 550-m-thick Mt. Simon Sandstone, which directly overlies granitic basement. The primary sealing cap rock is the Eau Claire Shale, a 100- to 150-m-thick unit at a depth of roughly 1.5 km. The USGS seismic network consists of 12 stations, three of which have surface accelerometers and three-component borehole geophones. We derived a one-dimensional velocity models from a vertical seismic profile acquired by Archer-Daniels-Midland (ADM) and the Illinois State Geological Survey (ISGS) to a depth of 2.2 km, tied into shallow acoustic logs from our borehole stations and assuming a 6 km/sec P-wave velocity for granite below 2.2 km. We further assume a constant ratio of P- to S-wave velocities of 1.83, as derived from velocity model inversions. We use this velocity model to locate seismic events, all of which are within the footprint of our network. So far magnitudes of locatable events range from Mw = -1.52 to 1.07. We further improved the hypocentral precision of microseismic events when travel times and waveforms are sufficiently similar by employing double-difference relocation techniques, with relative location errors less than 80 m horizontally and 100 m vertically. We observe tend to group in three distinct clusters: ∼0.4 to 1.0 km NE, 1.6 to 2.4 km N, and ∼1.8 to 2.6 km WNW from the injection well. The first cluster of microseismicity forms a roughly linear trend, which may represent a pre-existing geologic structure. Most of these microearthquakes occur in the granitic basement at depths greater than 2.2 km, well below the caprock, and likely do not compromise the integrity of the seal. We conclude that because the observed microseismicity is occurring in the granitic basement, the integrity of the caprock seal has not been compromised by CCS activities.

Illinois↗

Micro-seismicity and seismic moment release within the Coso Geothermal Field, California

We relocate 16 years of seismicity in the Coso Geothermal Field (CGF) using differential travel times and simultaneously invert for seismic velocities to improve our knowledge of the subsurface geologic and hydrologic structure. We expand on our previous results by doubling the number of relocated events from April 1996 through May 2012 using a new field-wide 3-D velocity model. Relocated micro-seismicity sharpens in many portions of the active geothermal reservoir, likely defining large-scale fault zones and fluid pressure compartment boundaries. However, a significant fraction of seismicity remains diffuse and does not cluster into sharply defined structures, suggesting that permeability is maintained within the reservoir through distributed brittle failure. The seismic velocity structure reveals heterogeneous distributions of compressional (Vp) and shear (Vs) wave speed, with Vs generally higher in the Main Field and East Flank and Vp remaining relatively uniform across the CGF, but with significant local variations. The Vp/Vs ratio appears to outline the two main producing compartments of the reservoir at depths below mean ground level of approximately 1 to 2.5 km, with a ridge of relatively high Vp/Vs separating the Main Field from the East Flank. Detailed analyses of spatial and temporal variations in earthquake relocations and cumulative seismic moment release in the East Flank reveal three regions with persistently high rates of seismic activity. Two of these regions exhibit sharp, stationary boundaries at the margins of the East Flank that likely represent barriers to fluid flow and advective heat transport. However, seismicity and moment release in a third region at the northern end of the East Flank spread over time to form an elongated NE to SW structure, roughly parallel both to an elongated cluster of seismicity at the southern end of the East Flank and to regional fault traces mapped at the surface. Our results indicate that high-precision relocations of micro-seismicity and simultaneous velocity inversions in conjunction with mapping of seismic moment release can provide useful insights into subsurface structural features and hydrologic compartmentalization within the Coso Geothermal Field.

California↗

Linear complementarity formulation for 3D frictional sliding problems

Frictional sliding on quasi-statically deforming faults and fractures can be modeled efficiently using a linear complementarity formulation. We review the formulation in two dimensions and expand the formulation to three-dimensional problems including problems of orthotropic friction. This formulation accurately reproduces analytical solutions to static Coulomb friction sliding problems. The formulation accounts for opening displacements that can occur near regions of non-planarity even under large confining pressures. Such problems are difficult to solve owing to the coupling of relative displacements and tractions; thus, many geomechanical problems tend to neglect these effects. Simple test cases highlight the importance of including friction and allowing for opening when solving quasi-static fault mechanics models. These results also underscore the importance of considering the effects of non-planarity in modeling processes associated with crustal faulting.

Computational Geosciences↗

Mechanical analysis of fault slip data: Implications for paleostress analysis

Stress inversions are a useful and popular tool for structural geologist and seismologist alike. These methods were first introduced by Wallace (1951) and Bott (1959) and subsequent studies continue to be based on their assumptions: the remote stress tensor is spatially uniform for the rock mass containing the faults and temporally constant over the history of faulting in that region, and the slip on each fault surface has the same direction and sense as the maximum shear stress resolved on that surface from the remote stress tensor. Furthermore, successful implementation requires that slip accumulates on faults of diverse orientation. Many studies employ these methods on isolated faults or on fault systems with limited ranges of orientations, which can lead to erroneous results. We propose a new method that incorporates the effects of mechanical interaction of the entire fault or fault system, and solves the complete mechanical problem rather than employing empirical relationships between slip and stress or strain (or strain rate). The method requires knowledge of the fault geometry and information on at least one slip vector component along portions of the known fault geometry. For example, if throw is known, the strike-slip component can be solved for. We test the method using a single synthetic fault with anisotropic roughness similar to that measured at fault outcrops. While the orientation of remote stress may be determined precisely, the lack of diverse fault orientations introduces a systematic error in the remote stress ratio. We further test the effect of diversity of fault orientations and find that Wallace–Bott type inversions do not perform as well for limited ranges of orientations when compared to the proposed method. Finally, we use published data from the 1999 Chi-Chi, Taiwan, earthquake, and find that the method using surface data only, and surface data with subsurface focal mechanisms, produce similar results. The resulting stress orientations are in good agreement with results from Wallace–Bott inversions. Furthermore, the slip distribution is in general agreement with kinematic slip inversions using coseismic surface deformation. Stress inversion methods using fault slip data can thus be improved upon, significantly in some cases, by solving a mechanical boundary value problem that takes into account the geometry of faults or fault systems. As a bonus, the solution provides the stress, strain, and displacement fields throughout the region and the slip distributions on the faults.

Journal of Structural Geology↗