Vertical Fault Research Articles

Application of electrical resistivity tomography for mapping gold mineralization potential in Iperindo, Ilesha schist belt, southwestern Nigeria

This study applied Electrical Resistivity Tomography geophysical investigation technique to evaluate the gold mineralization potential of Iperindo in Ilesha Schist Belt, southwestern Nigeria, where commercial exploitation capable of generating revenue and employment for the inhabitants has been challenged by lack / inadequate subsurface geological/geophysical information. The filtered and inverted electrical resistivity data acquired through five (5) 336 m long E – W trending profiles, established 10 m apart from each other, delineate isolated near surface but thick (> 30 m) low resistivity zones, especially at the eastern and western ends of the study area. Some of the delineated low resistivity zones (3 – 200 𝛀m) present vertical sharp edges, likely created by vertical faults that flank the zones on both sides. The low resistivity of these zones could be attributed to the occurrence of conductive material such as gold and associated base metals which probably exist in pegmatitic veins within the zones.  Â

3D Finite-Element Modeling of Dynamic Rupture and Aseismic Slip over Earthquake Cycles on Geometrically Complex Faults

ABSTRACTWe develop a new dynamic earthquake simulator to numerically simulate both spontaneous rupture and aseismic slip over earthquake cycles on geometrically complex fault systems governed by rate- and state-dependent friction. The method is based on the dynamic finite-element method (FEM) EQdyna, which is directly used in the simulator for modeling 3D spontaneous rupture. We apply an adaptive dynamic relaxation technique and a variable time stepping scheme to EQdyna to model the quasi-static processes of an earthquake cycle, including the postseismic, interseismic, and nucleation processes. Therefore, the dynamic and quasi-static processes of an earthquake cycle are modeled in one FEM framework. Tests on a vertical strike-slip fault verify the correctness of the dynamic simulator. We apply the simulator to thrust faults with various dipping angles, which can be considered as the simplest case of geometrically complex faults by breaking symmetry, compared with vertical faults, to examine effects of dipping fault geometry on earthquake cycle behaviors. We find that shallower dipping thrust faults produce larger seismic slip and longer recurrence time over earthquake cycles with the same rupture area. In addition, we find an empirically linear scaling relation between the recurrence interval (and the seismic moment) and the sinusoidal function of the dip angle. The dip-angle dependence is likely due to the free-surface effect, because of broken symmetry. These results suggest dynamic earthquake simulators that can handle nonvertical dipping fault geometry are needed for subduction-zone earthquake studies.

Main faults of the Betic Cordillera

In the Betic Cordillera the faults correspond to thrusts forming tectonic units and to three important sets of vertical or quasi-vertical faults. In the Internal Zone the thrusts formed mainly during the Oligocene, under an approximately N-S compression. Later, in the Burdigalian, the Internal Zone collided with the External Zone, which underwent major deformations. Then N60ºE-S60ºW to E-W dextral strike-slip faults formed. These faults formed under nearly E-W compression. From the end of the middle Miocene, the near N-S compression was reestablished and then formed two conjugated sets of faults whose directions are: NNE-SSW (with sinistral character) and NNW-SSE to NW-SE (with dextral character). Not only are the aforementioned major faults, but the rocks also have many other small faults and joints that formed under the same stresses.

On the morphology and amplitude of 2D and 3D thermal anomalies induced by buoyancy-driven flow within and around fault zones

Abstract. In the first kilometers of the subsurface, temperature anomalies due to heat conduction processes rarely exceed 20–30 ∘C. When fault zones are sufficiently permeable, fluid flow may lead to much larger thermal anomalies, as evidenced by the emergence of thermal springs or by fault-related geothermal reservoirs. Hydrothermal convection triggered by buoyancy effects creates thermal anomalies whose morphology and amplitude are not well known, especially when depth- and time-dependent permeability is considered. Exploitation of shallow thermal anomalies for heat and power production partly depends on the volume and temperature of the hydrothermal reservoir. This study presents a non-exhaustive numerical investigation of fluid flow models within and around simplified fault zones, wherein realistic fluid and rock properties are accounted for, as are appropriate boundary conditions. 2D simplified models point out relevant physical mechanisms for geological problems, such as “thermal inheritance” or pulsating plumes. When permeability is increased, the classic “finger-like” upwellings evolve towards a “bulb-like” geometry, resulting in a large volume of hot fluid at shallow depth. In simplified 3D models wherein the fault zone dip angle and fault zone thickness are varied, the anomalously hot reservoir exhibits a kilometer-sized “hot air balloon” morphology or, when permeability is depth-dependent, a “funnel-shaped” geometry. For thick faults, the number of thermal anomalies increases but not the amplitude. The largest amplitude (up to 80–90 ∘C) is obtained for vertical fault zones. At the top of a vertical, 100 m wide fault zone, temperature anomalies greater than 30 ∘C may extend laterally over more than 1 km from the fault boundary. These preliminary results should motivate further geothermal investigations of more elaborated models wherein topography and fault intersections would be accounted for.

Evolution of Faulting Induced by Deep Fluid Injection, Paradox Valley, Colorado

ABSTRACT High-pressure fluid injection into a subhorizontal confined aquifer at 4.3–4.6 km depth induced >7000 earthquakes between 1991 and 2012 within once seismically quiescent Paradox Valley in Colorado, with magnitudes up to Mw 3.9. Earthquake hypocenters expanded laterally away from the well with time, defining the margins of the aquifer pressurized by injection at the well. Within 5 km of the well, alignment of earthquake hypocenters defines strikes of nine vertical fault zones. Previous studies show that these fault zones predate injection, producing left-stepping offsets in the normal faults of the Wray-Mesa fault system that cradles Paradox Valley. Hypocenters, rakes, and strikes of 2041 well-constrained focal mechanisms show that most injection-related earthquakes occur where these vertical faults intersect the pressurized aquifer. Well-defined focal mechanisms show that this induced seismicity consists of Riedel shear faults at acute angles to the strikes of these fault zones. These small faults develop an anastomosing fault structure of focal planes along each planar fault zone, as fluid injection continues, even as their hypocenters define a single planar fault zone. Failure conditions at each hypocenter are found using a fully coupled poroelastic analysis of stress induced by fluid injection, and this analysis indicates a minimum Coulomb failure condition of 0.1 MPa. This failure condition is primarily a result of aquifer pore-fluid pressurization, as almost all well-located seismicity is within the pressurized aquifer. Reducing the rate of injection and frequent well shutdowns in the second decade nearly eliminated induced seismicity, except very near the well where gradients in pressurization are the largest. Despite these decreases in failure conditions and seismicity, some fault zones continued to produce earthquakes larger than M 3 as injection continued.

Automated Microseismic Processing and Integrated Interpretation of Induced Seismicity during a Multistage Hydraulic-Fracturing Stimulation, Alberta, Canada

ABSTRACT The development of organic-rich, low-permeability formations for hydrocarbon production requires the use of unconventional techniques such as multiwell pad drilling of horizontal wells and massive multistage hydraulic-fracturing stimulations. However, proliferation of these unconventional development methods has been linked to localized cases of fault reactivation during or shortly after hydraulic fracturing. In the Duvernay formation, located in Alberta, Canada, induced seismicity from hydraulic fracturing has occurred on nearly vertical strike-slip faults that are difficult to detect with conventional seismic exploration methods. In such cases, faults may only be discernible from seismic events with precise and accurate locations, which generally requires dense seismic monitoring arrays deployed near the stimulated wells. In this study, we introduce a new, semiautomated workflow for processing passive seismic data from a dense array and then integrate it with a 3D seismic dataset to characterize seismicity clusters related to hydraulic fracturing and pre-existing faults. The reactivated faults inferred from the distribution of the microseismic events directly overlie a system of incised, middle Devonian channels below the Duvernay formation observed in time slices extracted from the 3D seismic data. The channel system exhibits a set of lateral offsets, interpreted as ancient strike-slip fault displacements, the detection of which is further enhanced by use of a similarity attribute calculated from the 3D seismic data. Taken together, integrated interpretation of induced seismicity and 3D seismic data support a model of a regional left-lateral strike-slip fault system that was active during the middle Devonian and reactivated in a reverse sense (right-lateral strike slip) during hydraulic-fracturing operations.

Uncertainty assessment for 3D geologic modeling of fault zones based on geologic inputs and prior knowledge

Abstract. Characterizing the zone of damaged and altered rock surrounding a fault surface is highly relevant to geotechnical and geo-environmental engineering works in the subsurface. Evaluating the uncertainty associated with 3D geologic modeling of these fault zones is made possible using the popular and flexible input-based uncertainty propagation approach to geologic model uncertainty assessment – termed probabilistic geomodeling. To satisfy the automation requirements of probabilistic geomodeling while still preserving the key geometry of fault zones in the subsurface, a clear and straightforward modeling approach is developed based on four geologic inputs used in implicit geologic modeling algorithms (surface trace, structural orientation, vertical termination depth and fault zone thickness). The rationale applied to identifying and characterizing the various sources of uncertainty affecting each input are explored and provided using open-source codes. In considering these sources of uncertainty, a novel model formulation is implemented using prior geologic knowledge (i.e., empirical and theoretical relationships) to parameterize modeling inputs which are typically subjectively interpreted by the modeler (e.g., vertical termination depth of fault zones). Additionally, the application of anisotropic spherical distributions to modeling disparate levels of information available regarding a fault zone's dip azimuth and dip angle is demonstrated, providing improved control over the structural orientation uncertainty envelope. The probabilistic geomodeling approach developed is applied to a simple fault zone geologic model built from historically available geologic mapping data, allowing for a visual comparison of the independent contributions of each modeling input on the combined model uncertainty, revealing that vertical termination depth and structural orientation uncertainty dominate model uncertainty at depth, while surface trace uncertainty dominates model uncertainty near the ground surface. The method is also successfully applied to a more complex fault network model containing intersecting major and minor fault zones. The impacts of the model parameterization choices, the fault zone modeling approach and the effects of fault zone interactions on the final geologic model uncertainty assessment are discussed.

Minerals and ore deposits exploration using meta-heuristic based optimization on magnetic data

The difficulties in unravelling the tectonic structures, in some cases, prevent the understanding of the ore bodies' geometry, leading to mistakes in mineral exploration, mine planning, evaluation of ore deposits, and even mineral exploitation. For that reason, many geophysical techniques are introduced to reveal the type, dimension, and geometry of these structures. Among them, electric methods, self-potential, electromagnetic, magnetic and gravity methods. Global meta-heuristic technique using Whale Optimization Algorithm (WOA) has been utilized for assessing model parameters from magnetic anomalies due to a thin dike, a dipping dike, and a vertical fault like/shear zone geological structure. These structures are commonly associated with mineralization. This modern algorithm was firstly applied on a free-noise synthetic data and to a noisy data with three different levels of random noise to simulate natural and artificial anomaly disturbances. Good results obtained through the inversion of such synthetic examples prove the validity and applicability of our algorithm. Thereafter, the method is applied to real case studies taken from different ore mineralization resembling different geologic conditions. Data are taken from Canada, United States, Sweden, Peru, India, and Australia. The obtained results revealed good correlation with previous interpretations of these real field examples.

Tectonic Inheritance With Dipping Faults and Deformation Fabric in the Brittle and Ductile Southern California Crust

AbstractPlate motions in Southern California have undergone a transition from compressional and extensional regimes to a dominantly strike‐slip regime in the Miocene. Strike‐slip motion is most easily accommodated on vertical faults, and major transform fault strands in the region are typically mapped as near vertical on the surface. However, some previous work suggests that these faults have a dipping geometry at depth. We analyze receiver function arrivals that vary harmonically with back azimuth at all available broadband stations in the region. The results show a dominant signal from contrasts in dipping foliation as well as dipping isotropic velocity contrasts from all crustal depths, including from the ductile middle to lower crust. We interpret these receiver function observations as a dipping fault‐parallel structural fabric that is pervasive throughout the region. The strike of these structures and fabrics is parallel to that of nearby fault surface traces. We also plot microseismicity on depth profiles perpendicular to major strike‐slip faults and find consistently NE dipping features in seismicity changing from near vertical (80–85°) on the Elsinore Fault in the Peninsular Ranges to 60–65° slightly further inland on the San Jacinto Fault to 50–55° on the San Andreas Fault. Taken together, the dipping features in seismicity and in rock fabric suggest that preexisting fabrics and faults may have acted as strain guides in the modern slip regime, with reactivation and growth of strike‐slip faults along northeast dipping fabrics both above and below the brittle‐ductile transition.

Seismicity in Strike‐Slip Foreland Faults (Central Betic Cordillera Front): Evidence of Indentation Tectonics

AbstractUnexpected seismicity has been detected in the past few years along the northern Guadiana Menor River (Guadalquivir foreland basin) at the northern boundary of the Betic Cordillera. Earthquake focal mechanisms evidence the activity of N‐S to NNE‐SSW sinistral faults in the basement. Yet continuous GPS (CGPS) data show a westward movement in both the Prebetic Arc and the eastern Guadalquivir basin infill, which disagrees with the strike‐slip faults. To more precisely describe the structure and evolution of the area, new Bouguer anomaly data from the southeastern basin combined with seismicity data, electric resistivity tomography profiles (ERT), and surface studies are provided. Given the soft consistency of the sedimentary infill in the area, surface evidence of faulting is scarce, limited to elongated channel and minor vertical faults affecting Quaternary sediments, where ERT profiles suggested the presence of faults. These results suggest that the Guadalquivir basin infill, the Prebetic Arc and the central Betic Cordillera move towards the west, independent from the Iberian crust, but indicating to some extent an upward propagation of the basement deformation. However, the basement is coupled with the Alboran Domain and undergoes NNW‐SSE Eurasia‐Nubia convergence and orthogonal extension. This scenario is suggestive of the initial stages of indentation tectonics, better developed in the eastern Betic Cordillera and the central Alboran Sea.

On the impact of alternative seismic time imaging methods on subsurface fault mapping in the northernmost Molasse Basin, Switzerland

The identification and characterization of tectonic faults in the subsurface represent key aspects of geologic exploration activities across the world. We have evaluated the impact of alternative seismic time imaging methods on initial subsurface fault mapping in three dimensions in the form of a case study situated in the most external foreland of the European Central Alps (the northernmost Molasse Basin). Four different seismic amplitude volumes of one and the same 3D seismic data set, differing in imaging technologies and parameterizations applied, were considered for the interpretation of a fault zone dissecting a Mesozoic sedimentary sequence that is characterized by a pronounced mechanical stratigraphy and has witnessed a multiphase tectonic evolution. For this purpose, we interpreted each seismic amplitude volume separately. In addition, we computed a series of seismic attributes individually for each volume. Comparison of the different data interpretations revealed consistent results concerning the mapping of the seismic marker horizons and main fault segments. Deviations concern the apparent degree of vertical and lateral fault zone segmentation and the occurrence of small-scale fault strands that may be regarded as important fault kinematic indicators. The compilation of all fault interpretations in map form allows the critical assessment of the robustness of the initial seismic fault mapping, highlighting well-constrained from poorly defined fault zone elements. We conclude that the consideration of multiple seismic processing products for subsurface fault mapping is advisable to evaluate general imaging uncertainties and potentially guide the development of fault zone model variants to tackle previously discussed aspects of conceptual interpretation uncertainties.

FD3D_TSN: A Fast and Simple Code for Dynamic Rupture Simulations with GPU Acceleration

AbstractWe introduce FD3D_TSN—an open-source Fortran code for 3D dynamic earthquake rupture modeling based on the staggered grid fourth-order finite-difference method employing a regular cubical spatial discretization. Slip-weakening and fast-velocity-weakening rate-and-state fault friction laws are combined with vertical planar fault geometry, orthogonal to a planar free surface. FD3D_TSN demonstrates good agreement with other methods in a range of benchmark exercises of the Southern California Earthquake Center and U.S. Geological Survey dynamic rupture code verification project. Efficient graphic processing units (GPU) acceleration using the OpenACC framework yields a factor of 10 speed-up in terms of time to solution compared to a single-core solution for current hardware (Intel i9-9900K and Nvidia RTX 2070). The software is fast and easy-to-use and suitable explicitly for data-driven applications requiring a large number of forward simulations such as dynamic source inversion or probabilistic ground-motion modeling. The code is freely available for the scientific community and may be incorporated in physics-based earthquake source imaging and seismic hazard assessment, or for teaching purposes.

A comparison study using particle swarm optimization inversion algorithm for gravity anomaly interpretation due to a 2D vertical fault structure

A new approach to the inversion of gravity data utilizing the Particle Swarm Optimization (PSO) algorithm is used to model 2D vertical faults. The PSO algorithm is stochastic in nature; its development was motivated by the communal in-flight performance of birds looking for food. The birds are represented by particles (or models). Individual particles have a location and a velocity vector. The location vectors represent the parameter value. PSO is adjusted with random particles (models) and searches for targets by updating generations.Herein, the PSO algorithm is applied to three synthetic data sets (residual only with and without noise, residual plus regional, residual plus anomaly generated by a buried cylinder structure) and two field gravity data sets acquired across known faults in Egypt. Assessment of the synthetic data demonstrates that the PSO algorithm generates superior results if a first horizontal gradient (FHG) filter is applied first. The robustness of the PSO inversion algorithm is demonstrated for both synthetic and field gravity data.

Potential of Surface-to-Tunnel Seismic Tomography to Detect Vertical Faults: Application to the Tournemire Underground Research Laboratory, France

The Tournemire Underground Research Laboratory (Aveyron, France), developed by the French Institute for Radiological Protection and Nuclear Safety, is composed of a tunnel excavated in an argillaceous geological layer belonging to a limestone–argillite–limestone sub-horizontal sedimentary sequence. Sub-vertical fault zones are intercepted by drifts and boreholes excavated from the tunnel in the argillaceous layer at a depth of about 250 m. These fault zones are characterized by velocities lower than those of the embedding rocks. In this study, we assess the capacity of first arrival traveltime tomography to detect and characterize such fault zones in a surface-to-depth transmission configuration with a limited angle aperture. The source at the surface is a sledgehammer and turns out to be very repetitive. The source trigger is transmitted to the acquisition bay via a VHF transmitter. Receivers are located at depth in underground works. First arrivals are clear and traveltime picking is straightforward. Before inverting real data, we conduct synthetic tests with the real acquisition geometry to assess limits and validity of first arrival traveltime tomography to detect a lower velocity fault zone using different initial configurations and inversion constraints during the inversion. We show that introducing a priori information on velocities in the shallow subsurface, limestones and argillaceous layers and constraints on maximum velocities greatly helps improving results and subsequent interpretation. Velocities obtained from the real data display a main lower velocity zone that extends from the galleries in the argillaceous layer right up to the surface in limestones. The location and extension of this lower velocity zone are consistent with fault and fracture zones mapped from geological observations.

Use of prismatic waves in full-waveform inversion with the exact Hessian

The truncated Newton method uses information contained in the exact Hessian in full-waveform inversion (FWI). The exact Hessian physically contains information regarding doubly scattered waves, especially prismatic events. These waves are mainly caused by the scattering at steeply dipping structures, such as salt flanks and vertical or nearly vertical faults. We have systematically investigated the properties and applications of the exact Hessian. We begin by giving the formulas for computing each term in the exact Hessian and numerically analyzing their characteristics. We show that the second term in the exact Hessian may be comparable in magnitude to the first term. In particular, when there are apparent doubly scattered waves in the observed data, the influence of the second term may be dominant in the exact Hessian and the second term cannot be neglected. Next, we adopt a migration/demigration approach to compute the Gauss-Newton-descent direction and the Newton-descent direction using the approximate Hessian and the exact Hessian, respectively. In addition, we determine from the forward and the inverse perspectives that the second term in the exact Hessian not only contributes to the use of doubly scattered waves, but it also compensates for the use of single-scattering waves in FWI. Finally, we use three numerical examples to prove that by considering the second term in the exact Hessian, the role of prismatic waves in the observed data can be effectively revealed and steeply dipping structures can be reconstructed with higher accuracy.

Normal fault networks and their spatial relationships in Plio-Quaternary sedimentary series: A case study in the Zanjan Depression, NW Iran

Well-exposed normal fault networks in the Plio-Quaternary sedimentary series of the Zanjan Depression, Iran, were studied in order to analyse the spatial relationships and evolution of fault and damage zone structures in low-burial (~200 m) conditions. Despite the poorly-consolidated, geologically young nature of sediment, fault networks show several characteristics more typical of faulting in lithified strata. Key controls on the evolution of fault networks are found to be fault interactions, segment linkage, and the mechanical stratigraphy resulting in dip refraction, fault segmentation, the development of fault bends and steps, vertical fault overlap zone breaching, and bed-parallel slip. The data show a decrease of fracture frequency away from the fault cores. The obtained data show that damage zones width increases with fault displacement in a non-linear manner. Furthermore, damage zones appear more deformed and wider in the hanging-wall than in the footwall. A plot of damage zone width versus total fault throw reveals a positive but non-linear correlation, and the best-fit function corresponds better to a power-law relationship. The relationship between displacement and length shows a sub-vertical growth path with large displacement gradients typical of faulting in partially-lithified strata, due to shale ductility accommodating fault tip-related strains. A model for the evolution of the normal fault system is proposed based on spatial distribution of deformation, particularly fault segmentation, segment linkage, and damage zone development, as influenced by initial fault spacing, mechanical stratigraphy and the low burial conditions.

Structural characterization and hydrocarbon prediction for the SB5M strike-slip fault zone in the Shuntuo Low Uplift, Tarim Basin

Newly discovered superdeep marine carbonate reservoirs with burial depths of more than 7000 m are found in the Shuntuo Low Uplift, Tarim Basin, China. Previous exploration experience indicates that strike-slip faults in these deposits played a key role in reservoir formation and hydrocarbon accumulation. This paper describes the structural and hydrocarbon characteristics of the SB5M strike-slip fault zone in the Shuntuo Low Uplift based on the study of high-quality 3-D seismic volumes and drilling data. The SB5M fault is a repeatedly reactivated strike-slip vertical fault zone. In the Paleozoic strata, the SB5M strike-slip fault zone is composed of sub-vertical strike-slip faults in the deep layer, a graben structure in the middle layer, and en-echelon normal faults in the shallow layer. The deep, middle and shallow faults developed sequentially at different times and formed as a result of three-stage evolution: the Middle Caledonian, Late Caledonian, and Middle-Late Hercynian. Fracture networks formed by repeated fault activity provided pathways for karst flow and hydrocarbons, which are conducive to the formation of fracture–cavity reservoirs. Due to the fault damage zone-fault core transition and the complex structural features, the reservoirs controlled by the SB5M strike-slip fault zone feature strong heterogeneity. The data sets presented and evolution models established in this paper provide important data for the exploration and development of the SB5M strike-slip fault zone and represent references for structural analysis and hydrocarbon exploration of strike-slip fault zones in the Tarim Basin and elsewhere.

Development of railway track condition monitoring from multi-train in-service vehicles

A cab-based track monitoring system has been developed which makes use of the existing on-board GSM-R cab radio present in the majority of trains operating in the UK. With the addition of a low-cost sensor, type, location and severity of the track defects are reported using the system. The system improves safety and network performance by efficiently directing maintenance crews to the location of defects, minimising time spent on maintenance and inspection. Initially, vehicle dynamic simulation was used to test the feasibility of the system for defect monitoring and to develop compensation factors for vehicle type and operating speed. Novel on-board signal processing techniques are also presented through comparison of vibration response from sites with known defects and outputs from Network Rail’s (NR) New Measurement Train (NMT). Good agreement was reported for track faults in relation to vertical and lateral alignment and dip faults. Statistically, good agreement has been demonstrated, suggesting that the data acquired could be used to provide an indication of track quality thereby improving network performance, reducing rough ride and leading to improved passenger comfort. Improvements in the measured and statistical correlation are anticipated through the use, of multi-train / multi-journey and machine learning methods.

Conjugate ruptures and seismotectonic implications of the 2019 Mindanao earthquake sequence inferred from Sentinel-1 InSAR data

In 2019, four strong earthquakes of Mw>6.4 occurred successively in Mindanao, Philippines. Based on the reports from the USGS and PHIVOLCS, these earthquakes were dominated by strike-slip ruptures. Whether these earthquakes are temporally and spatially related remained unknown. We characterized the coseismic displacement fields during the earthquake sequence using an InSAR technique with Sentinel-1 SAR data. The InSAR deformation measurements convincingly reveal that the four earthquakes produced distinct coseismic displacement patterns. We estimated the source parameters of the earthquakes with a two-step inversion strategy. The optimal model suggests that the earthquake sequence resulted from the reactivation of a conjugate fault structure that involves two nearly vertical left-lateral strike-slip faults and two high-angle right-lateral strike-slip faults. We calculated Coulomb stress changes from the earthquake sequence, suggesting that the previous strong earthquakes had significant stress-encouraging effects on the following events. The regional velocities based on the GPS analysis suggest that the formation of this conjugate structure is mainly due to the westward movement of the subducting Philippine Sea Plate. This earthquake sequence provides a seismotectonic background for subsequent strong earthquakes and helps to better understand the formation mechanisms and seismotectonic implications of conjugate structure rupturing.

Interpretation of resistivity data using 3D Euler deconvolution and Radially Averaged Power Spectrum

Different electrode configuration in resistivity measurements over the same geologic structures generally produce different anomaly patterns. This is related to the position of the structure concerning the electrodes. Given the above, a unified approach has been proposed using the concept of analytical signal to interpret the resistivity data as electric potential follows Laplace’s equation. We interpret the data sets using Euler deconvolution and Radially Averaged Power Spectrum (RAPS) to determine both lithological boundaries and units and compare the results with inverted resistivity section. We analysed the resistivity data using electrical resistivity tomography (ERT) technique over conducting dyke and vertical fault for Wenner and dipole–dipole arrays, respectively. The obtained structural indices for dyke and fault are 1.0 and 0.6, respectively. The results from Euler depth solutions match well with the inverted resistivity section. Subsequently, two field examples one each over ground water and mineral exploration were analysed. The delineated lineaments over ground water exploration matches fairly well with the available results. In addition, some additional lineaments are also mapped. These new features could be a zone of interest for a detailed survey. The depths from Euler depth solutions and RAPS are in agreement. However, no depth information was available earlier. The data for mineral exploration have been acquired by ERT technique with a profile length of ~500 m over Dhanjori Basin, Jharkhand, India. The Euler depth solution and RAPS indicated the presence of two interfaces at an average depth of 8 and 20 m and 7 and 21 m, respectively. The first interface is present all along the profile, whereas the second interface lies between a part of the profile. The location and depth of these are in broad agreement with the resistivity sections obtained from ERT/AMT and borehole data. A borehole in the vicinity of the survey area indicated the presence of two interfaces. The first interface coincides with the soil depth and the second interface coincides with the sulphide mineralization.